Recalculation of 23 mouse HDL QTL datasets improves accuracy and allows for better candidate gene analysis.

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.

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 ( 6,7 ). One limitation of traditional QTL analysis is low mapping resolution, which is a result of the limited genetic recombination possible in onegeneration backcrosses [i.e., A×(A×B)] and two-generation intercrosses [i.e., (A×B)×(A×B)]. The 95% confi dence 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 ( 8 ). As it can be assumed that there are, on average, 20 genes per cM ( 9 ), the number of candidate genes per QTL can be very large, making the identifi cation of the causative gene very diffi cult. 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 ( 10 ). 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) ( 11 ). Mapping QTL requires accurate genetic map information for both the relative order of markers and the distances between them ( 12 ). Recently, Shifman and colleagues published a new genetic map based on a large population of a heterogeneous stock ( 13 ). Cox and colleagues integrated a total of 7,080 standard, simplesequence length polymorphism (SSLP) markers to this single-nucleotide polymorphism (SNP)-based map, generating a corrected mouse genetic map ( 14 ). This new map Abstract 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 confi dence 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 identifi ed 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 fi ve candidates. -Ackert-Bicknell, C., B. Paigen, and R. Korstanje 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 ( 1 ). Several epidemiologic studies have demonstrated that HDL is a strong inverse predictor of CVD risk (2)(3)(4)(5). 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

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 ( 12 ) (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 ( 8 ). 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 signifi cant and suggestive QTL were determined in a cross-specifi c 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 signifi cant if it exceeded the P < 0.05 threshold. These thresholds were chosen because they are the widely accepted cutoffs for suggestive and signifi cant 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) ( 12 ). The interactive scan model identifi ed the most likely position of the sex-specifi c 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 specifi city of QTL was not confi rmed by further analysis, and thus all "sex-specifi city" 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-specifi c QTL were identifi ed using the strain of paternal grandmother (pgm) for each mouse as both an additive and interactive covariate, as was done for sex-specifi c 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 specifi c; no such QTL, however, were identifi ed in this study.

Resequencing of SNPs
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.

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 ( 15,16 ), and the other two were obtained from the QTL archive 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 ( 14 ).
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 higher 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 ( 12 ).
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.

Identifi cation 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. Specifi cally, this database was searched using the keyword "HDL" in the Gene/Marker Symbol/Name fi eld, "QTL" in the Type fi eld, and "Any" in the 'Chromosome fi eld; the No Limit box was checked under the heading of Maximum Returned. The cross in which each HDL QTL was mapped was identifi ed and compared with the list of crosses identifi ed 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. 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 ( 17 ). However, Apoa2 is not the only distal Chr 1 gene to infl uence HDL ( 18 ). Careful examination of the QTL mapped to distal Chr 1 in this study showed three clusters of comapping QTL: one cluster of fi ve 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 ( 19 ). 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), (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 confi rmed QTL that had been previously reported in the literature, but in other instances, our analyses yielded disparate results when compared with the published fi ndings ( Figs. 2 , 3 ). For example, we observed that some QTL moved signifi cantly 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 ( 14 ) and because of methodological differences in the way we conducted our analyses compared with the original ( 8 ).

New HDL QTL map is required to determine the number of independent HDL QTL
As previously reported, many of the QTL comap to the same loci, suggesting a common underlying causative gene ( 7 ). 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, a Uni = all F1 animals were produced using the direction as indicated in the fi rst column, with the fi rst strain being the female; bi = F1 animals were produced using both directions.  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 defi ned 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 .

Identifi cation of Etv6 as a candidate gene for an HDL QTL on distal mouse Chr 6
On the distal end of Chr 6, we identifi ed 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 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 confi dence interval for this QTL can be determined, the distal boundary cannot. By using the distal limit of the confi dence interval for the B6×NZW cross, we determined that the narrowest common region for this QTL extends from 42 to 56 cM . 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: 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 identifi ed as containing at least one SNP for which the alleles met the conditions defi ned 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 ( 21 ). While liver expression data are not available for two of the strains used to map QTL identifi ed 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 signifi cant differential expression was noted for Rbm20 . The differential expression profi le for Fam160b1 , Rbm20 , and Shoc2 showed the expected profi le 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 profi le 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 refl ected in expression profi les, we looked for coding region polymorphisms (nonsynonymous, nonsense) and splice site SNPs in which the allele pattern met the requirements defi ned in Equation 3 ; none were found, however, in the SNP databases available for these fi ve strains. We then examined the Sanger resequencing databases to fi nd any known nonsynonymous (Cn) SNPs in our 15 candidate genes in the classical inbred strains (i.e., not wild derived). Of the fi ve 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, 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.

SM D
NZB CAST B 2 6 ( 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 ( 20 ). Therefore, we examined the expression pattern for this gene in our previously published liver expression strain survey ( 21 ) 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 defi ning 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 infl uence protein levels ( 22 ).

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 signifi cantly from the previous map ( Fig. 3C ). In addition to shifts in the peaks of QTL identifi ed in several crosses, we identifi ed several new peaks in other crosses but could not fi nd 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 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 fi ve strains of interest for this QTL cluster and determined that the genotypes for the fi ve strains fi t 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. 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 fi ve strains of interest for this QTL cluster and determined that the genotypes for the fi ve strains fi t Equation 3 . Five Cn SNPs are listed

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 ( 23 ), 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 ( 10 ). Using these methods, either in isolation or in combination, has allowed for the identifi cation 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 identifi ed. 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 fi rst 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 fi rst 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 ( 24 ). 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 Ca 2+ permeability and has a neuroprotective function ( 25 ). 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 ( 26 ). While mechanism by which Grin3a infl uences HDL levels remains unknown, this gene was recently identifi ed to be associated with HDL in a human genomewide association study ( 27 ).
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 ( 28 ). 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 ( 29 ). 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 ( 30 ).
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 ( 31 ). 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 traffi c proteins ( 32 ). Genetic variation in this gene has been associated with both the development of type II diabetes ( 33,34 ) and risk for development of Alzheimer disease ( 35 ). 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 ( 36,37 ). 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 ( 38 ). 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 ( 39 ). Little is known about Trub1 ; it has been identifi ed 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 identifi cation 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 fi ve genes as likely candidate genes for QTL on Chr 4, 6, and 19, respectively.