|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

* The Jackson Laboratory, 600 Main St., Bar Harbor, ME 04609
Wake Forest University, Winston-Salem, NC
Published, JLR Papers in Press, March 1, 2004. DOI 10.1194/jlr.M300460-JLR200
1 To whom correspondence should be addressed. e-mail: rkorstan{at}jax.org
| ABSTRACT |
|---|
|
|
|---|
The dependence of QTL effects on sex suggests an important role for the sex hormones in HDL regulation. We recommend that sex should be explicitly accounted for in future studies in the genetics of HDL regulation in both mice and humans.
Abbreviations: apoB, apolipoprotein B; ER, estrogen receptors; QTL, quantitative trait locus; LDLR, LDL receptor; LOD, logarithm of the odds; PLTP, phospholipid transfer protein; SR-BI, scavenger receptor class B type I; VLDLR, VLDL receptor
Supplementary key words diet high density lipoprotein mice sex
| INTRODUCTION |
|---|
|
|
|---|
One approach to identify genes involved in HDL regulation is using quantitative trait locus (QTL) analysis in animal models. By crossing inbred strains that significantly differ in HDL levels, and subsequently testing for association between HDL levels and genetic markers in the progeny, many QTL involved in HDL levels have been identified in human, mouse, rat, and rabbit (513). A recent review summarized the QTL for HDL that have been found in mouse and human (5). At least 12 murine QTL are located in regions homologous to human QTL, suggesting that the underlying candidate genes may be the same. Several QTL found in mice are found repeatedly in different crosses, suggesting that a limited set of important genes control HDL levels in the mouse.
Despite the many studies involving HDL QTL in mice, not much attention has been paid to the differences that we observe between males and females and between different diets. Mouse inbred strains NZB/BlNJ (NZB) and SM/J (SM) differ in their HDL levels [3.63 and 1.42 mmol/l in females and 4.33 and 1.42 mmol/l in males, respectively, according to the Mouse Phenome Database (14)]. Previously, an intercross between these two strains identified a QTL on Chr 1 (92.3 cM) and two QTL, one each on Chr 5 at 62 cM and 72 cM (11). Another study using the same two strains, but in a backcross to NZB, confirmed the loci on Chr 5 (15). The current study, which also uses the NZB and SM inbred strains, differs from previous QTL studies in several ways. First, the cross is considerably larger than in previous studies. Second, HDL was measured prior to and at different time points following the introduction of an atherogenic diet. Lastly, we included both sexes in this study and explicitly tested for dependence of QTL effects on sex. The study confirms the QTL on Chr 1 and Chr 5 found previously between the two inbred strains and reports new QTL on chromosomes 3, 6, 8, 10, 15, 16, 17, 18, and 19. More importantly, this study explicitly characterizes the interaction of QTL with sex.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Lipid measurements
At 0, 6, 12, and 16 weeks of diet consumption, mice were fasted for 4 h, blood was collected by retroorbital bleeding into EDTA-coated tubes, and plasma was separated by centrifugation at 1,500 rpm for 5 min at 4°C. Plasma total cholesterol concentrations were measured by commercial colorimetric enzymatic assay as described previously (18). HDL was measured after selective precipitation of apolipoprotein B (apoB)-containing lipoproteins with polyethylene glycol (19). The results are expressed as mean ± SE in mmol/l.
DNA isolation and genotyping
DNA was isolated from a centimeter of tail digested overnight in 500 µl of 1x digestion buffer (50 mM Tris-Cl pH 8.0, 100 mM EDTA pH 8.0, 100 mM NaCl, 1% SDS, 1 mg/ml proteinase K) in a 55°C waterbath. Digests were mixed with 1 vol of 25:24:1 phenol-chloroform-isoamyl alcohol and centrifuged for 5 min at 14,000 rpm at room temperature. DNA was precipitated by adding 2 vol of 100% ethanol to the isolated aqueous phase. The dried DNA pellets were resuspended in 1 ml TE (10 mM Tris-HCl, 1 mM EDTA, pH 7.58.0). Genotyping by PCR using mouse MIT MapPairs primers (Research Genetics, Huntsville, AL) was carried out under standard conditions at an annealing temperature of 55°C. A list of markers used for this cross is available as an online supplement. Polymorphisms were detected either by capillary electrophoresis or by standard electrophoresis using 4% agarose (Nusieve) gels.
Statistical analyses
Genome scans
We performed genomewide scans and multiple QTL modeling using the method of Sen and Churchill (20). For simple genome scans, this approach is similar to the interval mapping procedure of Lander and Botstein (21) but is based on an imputation algorithm. An advantage of the imputation approach is that it can be generalized to incorporate covariates, such as sex, as well as include multiple QTL in genome scan analyses. Logarithm of the odds (LOD) ratio scores for genome scans were computed at 2 cM intervals across the genome, and significance was determined by permutation testing (22). Significant and suggestive QTL meet or exceed the 95% and 90% genomewide thresholds, respectively. We ran multiple scans on each phenotype, and the usual suggestive threshold [37%, as reviewed in ref. (23)] was deemed too liberal. The software package used in this study was PSEUDOMARKER release version 1.0 (24).
We carried out single-locus genome scans to detect QTL with main effects. These scans included sex as an additive covariate to account for overall differences in HDL levels between the sexes. We then carried out a second set of single-locus scans that included an interaction between sex and the putative QTL at each locus. These scans allow the QTL effect to differ between the two sexes and can be used to identify sex-specific QTL. The difference in LOD scores between these two scans constitutes a test for sex by QTL 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. Lastly, we subdivided the population by sex and scanned male and female data separately for the main effect QTL.
We carried out simultaneous genome scans of all pairs of markers in the search for epistatic interactions. The search strategy is described by Sen and Churchill (20) and Sugiyama et al. (25). We failed to detect any significant QTL by QTL interactions using stringent genomewide adjusted criteria.
For some of the QTL detected in the single-locus genome scans, the shape of the LOD curve suggested that there might be two or more QTL present. In order to resolve multiple-linked QTL, we carried out multi-locus scans on a per chromosome basis. The change in the LOD score between one QTL and two QTL models was used as a test statistic. In one case, we also tested two QTL versus three QTL models. Thresholds for these chromosome-wide searches were established by running simulations.
Multiple regression We selected all significant and suggestive loci and interactions with sex for incorporation into multiple regression analysis and followed a backward elimination strategy to retain only individually significant terms (P < 0.01). The percent of variance explained for each QTL is based on the adjusted (type III) sums of squares.
| RESULTS |
|---|
|
|
|---|
|
|
|
|
|
|
LOD = 2.54) for two QTL on the chow diet (Table 3), with one QTL interacting with sex (Table 4). Evidence was also found for two QTL (both interacting with sex) on Chr 17 on the atherogenic diet. Furthermore, the tests suggest two QTL each for Chr 16 and Chr 19 on chow and for Chr 18 on the atherogenic diet.
|
|
| DISCUSSION |
|---|
|
|
|---|
Second, the results from this study imply that it is important to differentiate between the two sexes in QTL analysis because four QTL were found to have effects only in female mice, and one has effects only in male mice. The phenomenon of differences between sexes is not restricted to HDL (e.g., atherosclerosis, gallstones, albuminuria), and it is not restricted to phenotypes that show an overall difference between the two sexes. To identify sex-specific QTL, it is not sufficient to carry out separate analyses of males and females. This approach can be misleading, as the subdivided populations will be much smaller than the combined population, and this increases the likelihood that real QTL effects would fail to be detected in one or both subsets of the data. Analysis of the entire population with and without a sex-by-QTL interaction provides a more appropriate basis for interpretation. Failure to detect a significant interaction suggests a common QTL, the null hypothesis, and a sex-specific QTL is declared only when there is substantial evidence to support the claim. The shortcoming of the subdivided analysis approach is highlighted by our analysis of Hdlq26 on chromosome 10, where the individual scans do not detect a QTL, but the analysis with sex as an interactive covariate does. This strategy also applies for humans, where QTL analysis is commonly carried out on mixed-sex populations.
It would be interesting to examine the QTL in female mice that do not have estradiol. Will the QTL that are present only in females disappear? Will QTL that are normally present only in males now appear in females as well? A recently developed method for chemically-induced menopause in mice by follicle depletion using vinylcyclohexene diepoxide (28) might be a good way to address these questions.
A strong gender effect can be seen not only in the F2 population, but also in the analysis of the reciprocal F1 populations. In females, the SMxNZB F1 animals resemble the SM parents, while the reciprocal NZBxSM F1 animals are intermediate between SM and NZB. In males, the SMxNZB F1 animals also resemble the SM parents (although HDL levels are significantly higher), and the reciprocal NZBxSM F1 animals resemble the NZB parents. These observations do seem to indicate that an interaction between the X-chromosome (animals with the X-chromosome from SM have lower HDL levels) and the mother (animals with an SM mother have lower HDL levels) may be acting in crosses between the strains SM and NZB. This nurture effect has been observed in NZO mice, which are closely related to NZB (29).
Third, the most important changes happen in the first weeks on the high fat diet. We observed different QTL between chow and 6 weeks of high-fat diet, which suggests that different genes are regulating HDL, but after 6 weeks, the HDL levels remain relatively stable.
Fourth, male mice have significantly higher HDL levels than female mice. In contrast, premenopausal human females have higher HDL levels and reduced cardiovascular disease compared with men (30). Several studies have shown that the relatively higher estrogen levels in women cause this difference. Estrogen may protect women against atherosclerosis by increasing HDL levels and stimulating reverse cholesterol transport. One way is through the binding of the estrogen receptors (ER), as homo- or heterodimers, to an estrogen response element found in the promoter of estrogen sensitive genes. Another way is through the binding of the ER to other transcription factors, resulting in enhancement or repression of transcription (31). In the human male population, higher levels of testosterone are associated with higher HDL levels (32), and it has been shown that testosterone treatment causes female mice to have HDL levels that resemble that of the male (27).
Most QTL that we describe have been reported previously; only five out of the 15 QTL found in this study are novel, while two of the novel QTL (Hdlq26 and Hdlq28) show overlap with the QTL recently found for phospholipid transfer protein (PLTP) activity in the same cross (33). PLTP is known to be involved in HDL metabolism, and it is possible that genes involved in PLTP activity are indirectly involved in HDL levels. These findings confirm our earlier statement that a limited set of genes control HDL levels in the mouse (5). Of special interest are the QTL found on Chr 1, 5, 15, and 16, which are located in regions homologous to human QTL for HDL on chromosomes 1q21-25, 13q12, 8q24, and 3q26, respectively (5). The identification of the underlying genes may provide therapeutic targets to increase HDL levels in humans. The gene underlying the Chr 1 QTL (Hdlq20) is most likely Apoa2 [Wang and Paigen (5), submitted] and is caused by an amino acid difference at position 61 in the protein (valine in NZB, alanine in SM). The QTL on Chr 3 (Hdlq21) overlaps with QTL found in previous crosses (10, 13) and a QTL for PLTP activity in this cross (33). We have tested the gene encoding the apoA-I binding protein (Apoa1bp) as a candidate for the PLTP-activity QTL by sequence comparison and expression analysis in the liver but found no evidence to support the candidacy of the gene (33). The gene encoding the LDL receptor (LDLR)-related protein associated protein 1 (Lrpap1) is a good candidate for Hdlq22, the most proximal QTL (at 18 cM) on Chr 5. An overexpression of this gene in the liver is associated with decreased apoA-I and HDL levels (34). The gene encoding the HDL receptor (HDL-R) scavenger receptor class B type I (SR-BI) has been suggested as a candidate gene for Hdlq1, the distal QTL (at 70 cM) on Chr 5 (10, 15). However, both studies found no significant differences in Scarb1 (formerly Srb1) expression, in sequence variation between CAST and B6, or between SM and NZB mice. Tcf1, the gene encoding the transcription factor HNF1A, also maps within the confidence interval of the QTL, and it has been demonstrated that HNF1A is an essential transcriptional regulator of HDL metabolism (35). Hdlq27 on Chr 15 overlaps with Pltpq4, a QTL for PLTP activity. Ppara expression in the liver was shown to be significantly higher (
4-fold) in male SM mice compared with male NZB mice fed a chow diet (33). The apoD gene (Apod) is located close to the peak of Hdlq28 on Chr 16. Polymorphism in this gene has been associated with HDL levels in the African population (36). A good candidate gene for Hdlq31 is the gene encoding endothelial lipase (Lipg). Overexpression of endothelial lipase in mice reduces the plasma concentration of HDL (37). Of more interest in this study are the sex-specific QTL. In considering candidate genes for these QTL, we have to keep in mind that these are probably genes that are either regulated at the hormone level or act in combination with sex-specific factors. The gene for apoF (Apof) is located near the peak of Hdlq26 on Chr 10. This gene has also been suggested as a candidate for Pltpq2, a QTL for PLTP activity. Interestingly, Hdlq26 has only been found in females, while Pltpq2 has been observed in males. Sequence comparison of Apof between SM and NZB shows a basepair change leading to an amino acid difference in the protein (33). However, it is difficult to find a sex-specific link with this amino acid difference. Hdlq23, which was specific for females on the chow diet, was mapped to Chr 6. Searching the MGI and Ensembl databases did not lead to any obvious candidate genes.
The QTL on Chr 17 (Hdlq29) overlaps with an unnamed QTL found between CAST/Ei and C57BL/6J (10). Abcg5/Abcg8 have been suggested as candidate genes for this QTL (10).
The locus on Chr 19 (Hdlq32) has been reported previously in a cross between MRL/lpr and BALB/cJ (8). Also, a QTL for cholesterol absorption (Chab5) has been mapped to this region (38). The first report (8) suggests that Fas underlies this QTL. MRL/lpr has a defective Fas gene, which leads to a breakdown of the central and/or peripheral tolerance, resulting in the failure to properly clear CD4/CD8 negative cells. The Fas gene is identified as the major defect underlying the autoimmune manifestation, especially autoantibody levels, which strongly correlated with HDL cholesterol levels (8). Toda, Wickham, and Sullivan (37) have shown gender differences in Fas expression and the alteration of Fas expression by testosterone treatment in MRL/lpr mice. The gene encoding the VLDL receptor (Vldlr) is also located in this region.
In a human population (Alberta Hutteries), genetic variation of VLDLR is associated with variation in HDL cholesterol levels (40). However, Vldlr/ mice have normal plasma lipoprotein levels (41).
In summary, we have identified 15 significant QTL involved in HDL cholesterol levels. Most of these loci have been found in crosses using other mouse inbred strains and in homologous human regions. There are currently 64 HDL QTL mapped in the mouse using 16 different crosses. Most QTL have been found multiple times, resulting in 23 confirmed locations and four locations found only once. Recurrence of the same QTL in different crosses suggests a limited number of key genes involved in the phenotype. Furthermore, concordance between species shows that 22 QTL found in the mouse are also found at the homologous regions in human, and one QTL is concordant between mouse, rat, and rabbit (5). The combining of all this data provides an excellent basis for gene identification.
After applying statistical analyses that investigate sex dependence of these QTL, we found five QTL to be sex-specific (four specific for female and one for male). A likely explanation is the influence of sex hormones (estrogen and testosterone) on the regulation of the genes that underlie these QTL. Our findings highlight the importance of taking gender into account in the analysis of HDL QTL in both animal models and human studies.
| ACKNOWLEDGMENTS |
|---|
Manuscript received November 3, 2003 and in revised form December 23, 2003.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
C. Doorenbos, S.-W. Tsaih, S. Sheehan, N. Ishimori, G. Navis, G. Churchill, K. DiPetrillo, and R. Korstanje Quantitative Trait Loci for Urinary Albumin in Crosses Between C57BL/6J and A/J Inbred Mice in the Presence and Absence of Apoe Genetics, May 1, 2008; 179(1): 693 - 699. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. M. Stylianou, J. P. Affourtit, K. R. Shockley, R. Y. Wilpan, F. A. Abdi, S. Bhardwaj, J. Rollins, G. A Churchill, and B. Paigen Applying Gene Expression, Proteomics and Single-Nucleotide Polymorphism Analysis for Complex Trait Gene Identification Genetics, March 1, 2008; 178(3): 1795 - 1805. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Llamas, S. Belanger, S. Picard, and C. F. Deschepper Cardiac mass and cardiomyocyte size are governed by different genetic loci on either autosomes or chromosome Y in recombinant inbred mice Physiol Genomics, October 19, 2007; 31(2): 176 - 182. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. S. Wang, E. E. Schadt, H. Wang, X. Wang, L. Ingram-Drake, W. Shi, T. A. Drake, and A. J. Lusis Identification of Pathways for Atherosclerosis in Mice: Integration of Quantitative Trait Locus Analysis and Global Gene Expression Data Circ. Res., August 3, 2007; 101(3): e11 - e30. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. E. Wergedal, C. L. Ackert-Bicknell, W. G. Beamer, S. Mohan, D. J. Baylink, and A. K. Srivastava 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 J. Lipid Res., August 1, 2007; 48(8): 1724 - 1734. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. L. Svenson, R. Von Smith, P. A. Magnani, H. R. Suetin, B. Paigen, J. K. Naggert, R. Li, G. A. Churchill, and L. L. Peters Multiple trait measurements in 43 inbred mouse strains capture the phenotypic diversity characteristic of human populations J Appl Physiol, June 1, 2007; 102(6): 2369 - 2378. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. A. Raby, M. E. Soto-Quiros, L. Avila, S. L. Lake, A. Murphy, C. Liang, E. Fournier, M. Spesny, J. S. Sylvia, A. Verner, et al. Sex-specific linkage to total serum immunoglobulin E in families of children with asthma in Costa Rica Hum. Mol. Genet., February 1, 2007; 16(3): 243 - 253. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Wittenburg, M. A. Lyons, R. Li, U. Kurtz, X. Wang, J. Mossner, G. A. Churchill, M. C. Carey, and B. Paigen QTL mapping for genetic determinants of lipoprotein cholesterol levels in combined crosses of inbred mouse strains, J. Lipid Res., August 1, 2006; 47(8): 1780 - 1790. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. K. Srivastava, S. Mohan, G. L. Masinde, H. Yu, and D. J. Baylink Identification of quantitative trait loci that regulate obesity and serum lipid levels in MRL/MpJ x SJL/J inbred mice J. Lipid Res., January 1, 2006; 47(1): 123 - 133. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Wang and B. Paigen Genetics of Variation in HDL Cholesterol in Humans and Mice Circ. Res., January 7, 2005; 96(1): 27 - 42. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Wang, R. Korstanje, D. Higgins, and B. Paigen Haplotype Analysis in Multiple Crosses to Identify a QTL Gene Genome Res., September 1, 2004; 14(9): 1767 - 1772. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||