If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
In the past 15 years, the quantitative trait locus (QTL) mapping approach has been applied to crosses between different inbred mouse strains to identify genetic loci associated with plasma HDL cholesterol levels. Although successful, a disadvantage of this method is low mapping resolution, as often several hundred candidate genes fall within the confidence interval for each locus. Methods have been developed to narrow these loci by combining the data from the different crosses, but they rely on the accurate mapping of the QTL and the treatment of the data in a consistent manner. We collected 23 raw datasets used for the mapping of previously published HDL QTL and reanalyzed the data from each cross using a consistent method and the latest mouse genetic map. By utilizing this approach, we identified novel QTL and QTL that were mapped to the wrong part of chromosomes. Our new HDL QTL map allows for reliable combining of QTL data and candidate gene analysis, which we demonstrate by identifying Grin3a and Etv6, as candidate genes for QTL on chromosomes 4 and 6, respectively. In addition, we were able to narrow a QTL on Chr 19 to five candidates.
In the past decade, high-density lipoprotein (HDL) cholesterol has emerged as a new potential therapeutic target for the treatment of cardiovascular disease (CVD). The key role of HDL as a carrier of excess cellular cholesterol in the reverse cholesterol transport pathway is believed to provide protection against atherosclerosis (
High-density lipoprotein cholesterol as a predictor of coronary heart disease risk. The PROCAM experience and pathophysiological implications for reverse cholesterol transport.
Coronary heart disease prediction from lipoprotein cholesterol levels, triglycerides, lipoprotein(a), apolipoproteins A-I and B, and HDL density subfractions: The Atherosclerosis Risk in Communities (ARIC) Study.
). Although much is known about HDL metabolism and reverse cholesterol transport, we still do not have a complete understanding of all the genetic factors involved in the regulation of HDL levels. Also, HDL is not a single entity but can be subdivided into different subpopulations. The relationships among these different subpopulations are not understood, and the genes responsible for the subpopulation differences among individuals are unknown.
Several approaches, including genetic mapping studies, have been used to identify novel genes involved in the regulation of HDL levels. In mice, the technique of quantitative trait locus (QTL) mapping, which uses crosses between different inbred strains of mice, has been heavily employed over the past 15 years. These results have been summarized in several reviews (
). One limitation of traditional QTL analysis is low mapping resolution, which is a result of the limited genetic recombination possible in one-generation backcrosses [i.e., A×(A×B)] and two-generation intercrosses [i.e., (A×B)×(A×B)]. The 95% confidence interval (CI), the interval in which the causative gene is most likely to reside, is usually very broad. For example, one large study of bone mineral density QTL found that the average CI width for traditionally mapped QTL was 32 cM (
Mouse BMD quantitative trait loci show improved concordance with human genome-wide association loci when recalculated on a new, common mouse genetic map.
), the number of candidate genes per QTL can be very large, making the identification of the causative gene very difficult. In the past few years, several methods have been developed that combine accumulated data from the different crosses, allowing for narrowing of the CI for QTL and reducing candidate gene lists (
). However, the success of these methods heavily depends on the accuracy of the QTL mapping.
The current standard genetic map for the mouse is curated and maintained by the Mouse Genome Informatics (MGI) Group at The Jackson Laboratory (www.informatics.jax.org) (
). Cox and colleagues integrated a total of 7,080 standard, simple-sequence length polymorphism (SSLP) markers to this single-nucleotide polymorphism (SNP)-based map, generating a corrected mouse genetic map (
). This new map resolved inconsistencies between the physical and genetic maps and is now the standard MGI genetic map, providing highly accurate genetic distances. A recent mapping study, in which the new and traditional genetic maps were compared, suggests that up to 20% of published QTL may have been mislocalized due to marker order and positioning errors in the old genetic map (
In addition, many advances have been made over the years in the statistics and insights underlying QTL mapping. For example, in the past, males and females were analyzed as one population, disregarding possible sex-associated differences in the phenotype, or as two separate populations, resulting in a loss of power. We now analyze males and females as one population, using sex as a covariate. This method results in power compared with analyzing them separately, and it takes into account the obvious differences between males and females. Another issue is the direction of the cross. While, for most crosses, all F1 animals were produced by crossing a female from strain A with a male from strain B, in some crosses, the F1 animals originated from reciprocal crosses (A×B and B×A) and the F2 progeny were analyzed regardless of the parental grandmother. Now that we are aware of the epigenetic differences that can occur depending on the direction of the cross, we use the paternal grandmother as a covariate when analyzing the data (
The different HDL QTL datasets have been analyzed in various ways using different genetic maps and often without the above-mentioned covariates. This has resulted in misplaced QTL, missed QTL, or QTL that were only suggestive because of the loss of power. In this study, we made a concerted effort to identify and collect as many published, historical HDL QTL datasets as possible and reanalyze them. While most of the data used in this study were generated by our research group, several datasets were kindly provided by other investigators. All data underwent rigorous quality control and were then analyzed using current statistical methods. In all cases, the new corrected genetic mouse map was used. By analyzing the different datasets in the same way and by using the same genetic map, we created a more consistent and complete mouse HDL QTL map. As we demonstrate in our study, this improved map will aid in the search for novel HDL genes.
METHODS
Identification of mapping-cross datasets
First, a literature search to identify published reports of mouse HDL QTL was done using the following key words: “HDL,”“QTL,” and “mouse.” Second, the Mouse Genome Informatics Database (www.informatics.jax.org/) was searched for HDL QTL using the Genes and Markers Query form. Specifically, this database was searched using the keyword “HDL” in the Gene/Marker Symbol/Name field, “QTL” in the Type field, and “Any” in the 'Chromosome field; the No Limit box was checked under the heading of Maximum Returned. The cross in which each HDL QTL was mapped was identified and compared with the list of crosses identified by conventional literature search. In short, the result was a list of mouse mapping crosses in which HDL QTL were mapped regardless of the method used to measure HDL.
Main-effects QTL mapping
Map positions for the markers for all datasets were updated to the new mouse genetic map using a mouse map converter tool (http://cgd.jax.org). All QTL analyses were done using the R/qtl software package (
) (R Version 2.6.2, qtl Library Version 1.09-43, www.rqtl.org/). Data were examined for phenotypic outliers and for genotyping errors, as previously described (
Mouse BMD quantitative trait loci show improved concordance with human genome-wide association loci when recalculated on a new, common mouse genetic map.
). A single locus main-scan for QTL was performed, and LOD scores were calculated at 2 cM intervals across the genome using the EM (or expectation-maximization) method in R/qtl for all datasets. The LOD thresholds for significant and suggestive QTL were determined in a cross-specific manner based on 1,000 permutations of the data. A QTL was considered to be suggestive if the LOD score exceeded the P < 0.63 threshold and significant if it exceeded the P < 0.05 threshold. These thresholds were chosen because they are the widely accepted cutoffs for suggestive and significant QTL.
For 5 of the 23 crosses, data were available for both male and female mice (Table 1). To account for the average differences between the males and females, we carried out scans using sex as an additive covariate. To identify sex-dependent QTL effects, we carried out additional scans using sex as an interactive covariate and computed the differences in LOD scores between these two scans (i.e., the ΔLOD) (
). The interactive scan model identified the most likely position of the sex-specific QTL. Calculating the ΔLOD score at the peak position is the secondary test for the QTL-by-sex interaction. This secondary test is carried out with no adjustment for multiple testing, and the threshold, based on the usual chi-square distribution of the likelihood ratio, is 2.0 on the LOD scale. However, sex specificity of QTL was not confirmed by further analysis, and thus all “sex-specificity” of QTL should be considered putative. For 6 of the 23 crosses, data were available for two different cross directions (i.e., B×A and A×B, Table 1). Putative cross-direction-specific QTL were identified using the strain of paternal grandmother (pgm) for each mouse as both an additive and interactive covariate, as was done for sex-specific QTL. In only one cross, both males and females were examined in crosses generated in a reciprocal fashion. In this last dataset, the interactive term of sex:pgm was examined to identify QTL that interacted with sex and were cross-direction specific; no such QTL, however, were identified in this study.
Uni = all F1 animals were produced using the direction as indicated in the first column, with the first strain being the female; bi = F1 animals were produced using both directions.
Quantitative trait loci analysis for plasma HDL-cholesterol concentrations and atherosclerosis susceptibility between inbred mouse strains C57BL/6J and 129S1/SvImJ.
Mapping genetic loci that regulate lipid levels in a NZB/B1NJxRF/J intercross and a combined intercross involving NZB/B1NJ, RF/J, MRL/MpJ, and SJL/J mouse strains.
a Uni = all F1 animals were produced using the direction as indicated in the first column, with the first strain being the female; bi = F1 animals were produced using both directions.
High-quality genomic DNA from the different inbred mouse strains was obtained from the Jackson Laboratory's DNA Resource (www.jax.org/dnares/). Primers were designed spanning the exon that contained the SNP and ordered from Integrated DNA Technologies. PCR and Sanger sequencing were performed using standard protocols. Sequence data were analyzed using Sequencher 4.9.
RESULTS
Datasets used in this study
From our literature search, we determined that 28 datasets would be of interest to us for this project. However, 5 of these datasets were no longer available due to loss of records or the inability to contact the lead authors. In sum, we obtained a total of 23 independent datasets (summarized in Table 1). Of these datasets, 19 were from our own laboratory. Of the remaining 4 datasets, 2 were kindly provided by Drs. Breslow and Lusis, respectively (
), and the other two were obtained from the QTL archive (www.qtlarchive.org), a free public access repository for QTL mapping datasets.
A large number of main-effects QTL were mapped
The analyses for single main effect QTL resulted in 143 QTL (Table 2, Fig. 1). In many instances, our analyses confirmed QTL that had been previously reported in the literature, but in other instances, our analyses yielded disparate results when compared with the published findings (Figs. 2, 3). For example, we observed that some QTL moved significantly on the same chromosome, the confidence interval width changed for some QTL, and peak LOD scores were different for many QTL. These dissimilarities were caused by the discrepancies between the genetic map we used and those used in the original studies (
Mouse BMD quantitative trait loci show improved concordance with human genome-wide association loci when recalculated on a new, common mouse genetic map.
Fig. 1Mouse chromosome map with the recalculated HDL QTL. Each bar represents the 95% confidence interval, with the black dot indicating the peak position of the QTL.
Fig. 2Chromosome 1 with the old map positions on the right and the new QTL map positions on the left. Each bar represents the confidence interval of a QTL within a cross, with the black dot indicating the peak. Bars with the same number on the left and right are the QTL in the same cross. Careful examination of the QTL mapped to distal Chr 1 in this study shows three discrete clusters of comapping QTL: one cluster of five QTL centered at 70 cM, one cluster of nine QTL centered at 80 cM, and one cluster of four QTL centered at 85 cM.
Fig. 3Chromosomes 4 (A), 6 (B), and 19 (C) with the old map positions on the right and the new QTL map positions on the left. Each bar represents the confidence interval of a QTL within a cross, with the black dot indicating the peak. Bars with the same color on the left and right are QTL in the same cross.
). The most striking example of HDL QTL clustering is observed on the distal end of Chromosome (Chr) 1 (Fig. 2). By using bioinformatics techniques, such as haplotyping, which leverages these comapping QTL clusters, we previously showed that a single nucleotide polymorphism (SNP) in Apoa2 is the most likely causative SNP underlying the QTL mapped to distal Chr 1 in all crosses where the two strains were polymorphic for this SNP (
). Careful examination of the QTL mapped to distal Chr 1 in this study showed three clusters of comapping QTL: one cluster of five QTL centered at 70 cM, one cluster of nine QTL centered at 80 cM, and one cluster of four QTL centered at 85 cM. The Apoa2 gene (80 cM) is the most likely candidate for this middle cluster, whereas Soat1 (67.7 cM) is the most likely candidate for the upper cluster (
). The gene underlying the most distal cluster remains unknown.
Grin3a is a candidate gene for an HDL QTL on mouse Chr 4
We mapped nine QTL to mouse Chr 4, including a cluster of four QTL with peak locations between 26 and 29 cM (Table 2, Fig. 3A). This QTL cluster consisted of QTL mapped in the following crosses (underlined strains indicate the high allele for each QTL): PERA×I (peak at 26 cM), B6×A (peak at 27 cM), D2×129 (peak at 27 cM), and 129×CAST (peak at 29 cM). The minimum shared interval for these QTL extended from 20 to 31 cM (38.6–56.2 Mb), a region containing 343 genes with assigned MGI IDs (http://www.informatics.jax.org), as determined using the Biomart tool available from the Ensembl genome database [http://www.ensembl.org/index.html, Ensembl Genes 69 dataset:Mus musculus genes (GRCm38)]. The raw centimorgan positions for the confidence interval were converted to genome coordinates (NCBI, Bld37) using the Mouse Map Converter tool (http://cgd.jax.org). We then examined this interval to identify SNPs that met the following allele conditions:
Eq. 1
The SNP dataset used for this analysis was the CGD-MDA1 dataset, which is publically available (http://phenome.jax.org); all SNP in the above QTL interval in this database were examined. While this SNP dataset does not contain all known mouse SNP, it is the most complete, publically available SNP dataset containing data for all of the strains listed in Equation 1. Another limitation is that we excluded all SNPs that did not have complete genotype information for our strains of interest. Therefore, until complete genome sequences for all involved strains is available, we cannot exclude the possibility that candidate genes are missed. Of the 3,483 SNP for which data was available for the minimal QTL interval, only 8 SNPs met the allele distribution pattern defined by Equation 1. These SNP were all located in a narrow region between 49.86 and 50.58 Mb. This narrow region does not contain the coding region for any gene. However, the most proximal of the SNPs is located 937 bp distal to the start codon for the Grin3a gene. The Grin3a gene runs on the negative strand, and our SNPs overlap with its predicted regulatory/promoter region. There are no polymorphisms in the coding region of Grin3a that match the strain pattern shown in Equation 1.
Identification of Etv6 as a candidate gene for an HDL QTL on distal mouse Chr 6
On the distal end of Chr 6, we identified a cluster of three comapping QTL with peak locations between 70 and 74 cM (Table 2, Fig. 3B). These QTL were mapped in the following crosses (underlined strains have the high allele): SM×NZB (peak at 70 cM), D2×CAST (peak at 72 cM), and B6×D2 (peak at 74 cM). The B6×D2 QTL was only a suggestive QTL, and the allele effects for this QTL were not conclusive. Very few mice were phenotyped in the B6×D2 cross, severely underpowering the mapping in this study. However, the D2 strain was used in another cross with a QTL mapping to this location, suggesting that this B6×D2 QTL is valid. Using methods similar to that which was used for the Chr 4 locus, we then attempted to narrow this QTL to a single candidate gene. The minimum shared interval for this cluster on distal Chr 6 extended from 64 to 76 cM (131.92–143.65 Mb). This relatively gene-dense region contains 124 known or predicted genes. Using SNP comparison such that the conditions in Equation 2 were met, we were able to narrow this QTL to a single region spanning from 134.00 to 134.09 Mb.
Eq. 2
This region contains a single gene, Etv6. No polymorphisms are found in the coding region of Etv6 that match the strain pattern shown in Equation 2. This gene appears to be moderately expressed in the liver (
) using the reported Q values (http://cgd.jax.org/datasets/expression/10strain.shtml), but it did not show differential expression between strains in the three crosses defining this QTL. Thus, it appears that this gene either does not function directly in the liver to modulate HDL or that the causative mutation(s) in this gene do not affect transcription levels but may still influence protein levels (
Cluster of QTL on distal Chr 19 can be narrowed to a small number of candidate genes
The new QTL map for Chr 19 differs significantly from the previous map (Fig. 3C). In addition to shifts in the peaks of QTL identified in several crosses, we identified several new peaks in other crosses but could not find some of the peaks reported in the literature. We did identify an interesting cluster of three QTL at the distal end of Chr 19 in the following crosses (underlined strains have the high allele): NZB×RF (peak at 49 cM), NZW×B6 (peak at 54 cM), and NZB×SM (peak at 56 cM). The last marker genotyped in the NZB×RF cross is located at 50 cM. The LOD score at this marker is nearly identical to the LOD score at the location of peak association. Thus, while the proximal boundary of the confidence interval for this QTL can be determined, the distal boundary cannot. By using the distal limit of the confidence interval for the B6×NZW cross, we determined that the narrowest common region for this QTL extends from 42 to 56 cM (48.5–59.5 Mb). As for the other QTL clusters described above, we used SNP comparison to narrow the list of candidate genes. The conditions for the comparison were:
Eq. 3
Unlike the two previous QTL clusters, we were unable to narrow this QTL to a single candidate gene. Fifteen known genes (Sorcs3, Sorcs1, Add3, Mxi1, 5830416P10Rik, Rbm20, Pdcd4, Shoc2, Acsl5, Vti1a, Tcf7l2, B230217O12Rik, Fam160b1, Trub1, and Atrnl1) and two predicted genes (1700001K23Rik and Gm16745) were identified as containing at least one SNP for which the alleles met the conditions defined in Equation 3.
To further narrow the list of genes for this QTL, we examined liver expression for strains included in our previously published gene expression strain survey (
). While liver expression data are not available for two of the strains used to map QTL identified in this cluster (NZW and RF), data are available for NZB, B6, and SM. Statistically significant differential expression was noted for Add3, Shoc2, and Fam160b1, and nearly statistically significant differential expression was noted for Rbm20. The differential expression profile for Fam160b1, Rbm20, and Shoc2 showed the expected profile such that the two low allele strains (B6 and SM) were expressed at levels different than that for the high allele NZB strain, but expression was not different when comparing B6 with SM. The expression profile for the Add3 gene did not follow this pattern, as expression was not different when comparing the high allele NZB strain to the low allele B6 strain.
As coding region polymorphisms can affect gene function in a way that is not necessarily reflected in expression profiles, we looked for coding region polymorphisms (nonsynonymous, nonsense) and splice site SNPs in which the allele pattern met the requirements defined in Equation 3; none were found, however, in the SNP databases available for these five strains. We then examined the Sanger resequencing databases to find any known nonsynonymous (Cn) SNPs in our 15 candidate genes in the classical inbred strains (i.e., not wild derived). Of the five strains of interest for this QTL cluster, only B6 has been resequenced to date. Nevertheless, given the number of strains that are resequenced, examination of this database should give a more complete listing of Cn SNP than other available databases. Nonsynonymous SNPs are listed in the Sanger database in Add3, Rbm20, Sorcs3, Sorcs1, Fam160b1, Atrnl1, and Trub1. The only Cn SNP listed in Add3 appeared to be private to the NON/ShiLtJ strain, and the only two in Rbm20 appeared to be private to the C57BL substrains (C57BL/6J and C57BL/6NJ). We then used the SIFT tool (http://sift.jcvi.org/) to determine whether the resulting amino acid changes for the remaining Cn SNPs were predicted to have a functional consequence(s). We found one such Cn SNP (rs37375751) in the Sorcs3 gene, but this SNP was not predicted to affect protein function. Four Cn SNPs are listed in the Sorcs1 gene in the classical inbred strains. One SNP (at 50255976 Mb, no rs number available, NCBI, Bld37) appeared to be a private SNP in the AKR/J strain. Two others, rs37453589 and a SNP at 50752804 Mb (no rs number available), were predicted to be tolerated. A fourth SNP in this gene, at 50752876 Mb (no rs number available), was predicted to result in an A-to-V substitution. This substitution was predicted not to be tolerated and thus had the potential to impact protein function. While this SNP appears quite polymorphic among inbred strains, the alanine amino acid at this location is quite conserved among other mammals. We resequenced this SNP in the five strains of interest for this QTL cluster and determined that the genotypes for the five strains fit Equation 3. Five Cn SNPs are listed in the Fam160b1 gene, for the classical inbred strains (rs36742147, rs31128381, rs38343270, rs37035338, and rs36675440). However, none of these SNPs is predicted to affect protein function. Likewise, a single Cn SNP is listed in the Atrnl1 gene (rs36943069), but again this SNP is expected to be tolerated with regard to protein function. A single Cn SNP is listed in the Trub1 gene. This SNP at 57532595 Mb (no rs number available), results in a change from S to F at amino acid position 40, which is predicted not to be tolerated with regard to protein function. While serine at this position is well conserved among mammals, a phenylalanine is found at this position in opossums. A number of strains are polymorphic for this SNP, including the NZO/HILtJ strain, a close relative of the NZB/B1NJ strain. We then resequenced this SNP in the five strains of interest for this QTL cluster and determined that the genotypes for the five strains fit Equation 3. In sum, the expression data provide support for Fam160b1, Shoc2, and Rbm20, and the presence of Cn SNPs provide support for Sorcs1 and Trub1 as candidates for this Chr 19 QTL. Further work will be required to further narrow this QTL.
DISCUSSION
QTL analysis is a reliable method to identify loci involved in the variation of a phenotype of interest, and since the key publication by Lander and Botstein (
), thousands of QTL studies have been performed for many different phenotypes and utilizing different species. A disadvantage of the method is the low resolution of mapping, as it is dependent on the number of recombinations in the population, which is inherently low in a two-generation cross between inbred strains. Because of this low resolution, a typical QTL spans 10–20 cM and contains hundreds of genes. We and others have developed statistical approaches that combine the data from multiple crosses to increase mapping resolution yielding narrower QTL intervals and bioinformatics methods, such as haplotyping, that allow to decrease the number of possible candidate genes within an interval (
). Using these methods, either in isolation or in combination, has allowed for the identification of causative genes for QTLs. However, these approaches strongly depend on accurate QTL mapping, and inclusion of an inaccurately mapped QTL or exclusion of an accurately mapped QTL can lead to false conclusions.
The literature reports a number of clusters of putatively comapping HDL QTL for which no candidate gene has been identified. For many of these clusters, our analysis suggests the existence of additional comapping QTL (i.e., an addition of QTL to a locus) or inappropriate conclusions regarding comapping (i.e., the eliminated QTL from a locus). In the first situation, the addition of a QTL can provide additional haplotypes, which might allow for further narrowing of the interval. In the second situation, candidate genes may have been inappropriately eliminated by using the “false” strains that were included in the haplotype comparison.
The advantages of our reanalysis are demonstrated in the three examples presented above. In the first example highlighting a QTL on Chr 4, we were able to narrow the interval to one single gene, Glutamate receptor ionotropic NMDA3A (Grin3a). Grin3a, which is also known as NR3A, codes for an inhibitory subunit of NMDA-type glutamate receptors. Traditional NMDA receptors are key players in glutamate-facilitated neurotransmission (
The NMDA receptor may participate in widespread suppression of circuit level neural activity, in addition to a similarly prominent role in circuit level activation.
). Studies have suggested that expression of Grin3a is highest in selected regions of the brain in the postnatal period and in the adult retina, where the GRIN3A/NR3A subunit is thought to inhibit NMDA receptor Ca2+ permeability and has a neuroprotective function (
). Similarly, expression of Grin3a is high in the neonatal kidney, but expression does persist into adulthood. As in the brain, Grin3a appears to serve a protective role for the cells of the collecting duct of the kidney (
). While mechanism by which Grin3a influences HDL levels remains unknown, this gene was recently identified to be associated with HDL in a human genome-wide association study (
In the second example, for a Chr 6 QTL, we were also able to narrow the interval to a single gene (Etv6) using the QTL data. Etv6 is an ubiquitously expressed gene encoding a transcription repressor, and ETV6 is a member of the ETS family of transcription factors (
). Most studies on ETV6 have focused on the fusion proteins that can arise from chromosomal rearrangements between Etv6 and other genes encoding transcription factors, such as Runx1, and which have been implicated as being causative in some forms of leukemia (
). While the mechanism by which ETV6 regulates HDL levels is unclear, studies of gene networks associated with ETV6 action have suggested that this transcription repressor regulates genes associated with cholesterol biogenesis (
Finally, in our third example, Chr 19, combining the newly calculated QTL narrowed the region to a cluster of 15 genes, for which additional evidence exists supporting 5 as candidates for this QTL. This QTL is different from the other two examples, as it seems to be exclusive to NZ strains. Both the NZB/B1NJ and the NZW/LacJ strains can be traced back to a common set of mice brought to Australia in the 1930s, and recent analysis of the relatedness of inbred strains shows that the NZW/LacJ and the NZB/B1NJ strains are indeed genetically very closely related (
). In all three crosses contributing to the QTL cluster on distal Chr 19, one of the NZ strains was the high allele, suggesting that the alleles underlying this QTL are private to the NZ strains. As the modern NZO/H1LtJ is related to both the NZW/LacJ and the NZB/B1NJ strains, it is possible that this strain also carries the private “NZ” allele at this locus. The addition of the NZO/H1LtJ strain to the haplotyping for this locus does eliminate the Add3, 170001K23Rik, Mxi1I, and Tcf7l2 genes. Future work must be conducted, however, to determine whether this QTL is indeed caused by private NZ alleles. The 5 candidate genes for this QTL are Fam160b1, Shoc2, Rbm20, Sorcs1, and Trub1. Of these 5, Sorcs1 is the most likely candidate gene based on the existing evidence. This gene codes for a Vsp10p-D receptor family member transmembrane protein, which functions to sort and traffic proteins (
). These two diseases are correlated with serum lipid levels and genetic variation in another Vsp10p-D receptor gene, Sort1, which has been linked to serum LDL levels (
Evaluation of the gene-age interactions in HDL cholesterol, LDL cholesterol, and triglyceride levels: the impact of the SORT1 polymorphism on LDL cholesterol levels is age dependent.
). The remaining 4 genes have been poorly studied, and thus it is unclear how or whether these genes could function to modulate serum HDL levels. Mutations in SHOC2 can cause Noonan syndrome, an autosomal dominant disorder characterized by short stature, down-slanted eyes, cardiovascular defects, and increased tumor risks (
). The Rbm20 gene codes for RNA binding protein/RNA splicing protein, which shows particularly high expression in muscle, heart, and brain, and recent studies have shown that null mutations in this gene are associated with hereditary cardiomyopathies (
). Little is known about Trub1; it has been identified as a tRNA pseudouridinilation gene.
In conclusion, our reanalysis of the QTL data for 23 crosses results in an improved map for HDL QTL and facilitates the identification of HDL genes in the mouse. We show that using this approach, we can identify previously missed QTL and locate QTL on different regions of a chromosome. The new HDL QTL map allows for reliable combining of QTL data and candidate gene analysis, which we demonstrate by identifying Grin3a, Etv6, and a cluster of five genes as likely candidate genes for QTL on Chr 4, 6, and 19, respectively.
High-density lipoprotein cholesterol as a predictor of coronary heart disease risk. The PROCAM experience and pathophysiological implications for reverse cholesterol transport.
Coronary heart disease prediction from lipoprotein cholesterol levels, triglycerides, lipoprotein(a), apolipoproteins A-I and B, and HDL density subfractions: The Atherosclerosis Risk in Communities (ARIC) Study.
Mouse BMD quantitative trait loci show improved concordance with human genome-wide association loci when recalculated on a new, common mouse genetic map.
The NMDA receptor may participate in widespread suppression of circuit level neural activity, in addition to a similarly prominent role in circuit level activation.
Evaluation of the gene-age interactions in HDL cholesterol, LDL cholesterol, and triglyceride levels: the impact of the SORT1 polymorphism on LDL cholesterol levels is age dependent.
Quantitative trait loci analysis for plasma HDL-cholesterol concentrations and atherosclerosis susceptibility between inbred mouse strains C57BL/6J and 129S1/SvImJ.
Mapping genetic loci that regulate lipid levels in a NZB/B1NJxRF/J intercross and a combined intercross involving NZB/B1NJ, RF/J, MRL/MpJ, and SJL/J mouse strains.