|
|
||||||||
1 Departments of Pathology and Laboratory Medicine
2 Orthopedic Surgery
3 Radiology
4 Medicine
5 Microbiology, Immunology and Molecular Genetics
6 Biomathematics, University of California, Los Angeles, California 90095
7 Department of Bioinformatics, Rosetta Inpharmatics, Kirkland, Washington 98034
| ABSTRACT |
|---|
|
|
|---|
adipose tissue; plasma lipids; quantitative trait locus analysis; leptin; biomechanics; genetic pleiotropy
| INTRODUCTION |
|---|
|
|
|---|
The mouse has proved to be a useful model for the genetic dissection of complex traits. Quantitative trait locus (QTL) analysis in intercrosses or backcrosses of selected inbred strains allows survey of the entire genome for loci containing genes that influence a particular trait (26). We have previously used this technique to identify genetic loci affecting myocardial injury and calcification, lipid metabolism, and obesity, among other complex traits in mice fed an atherogenic diet (17, 27, 32). Inbred mouse strains have been shown to have significant heritable differences in a variety of bone phenotypes, including bone mass and susceptibility to osteoporosis (2, 6, 12, 22, 35, 55). Application of the technique of QTL analysis in mice for peak bone mass has begun to identify genetic loci associated with this trait (5, 7, 23, 47, 51).
In the above instances, which are typical of most QTL analyses, single or several closely related traits were examined, with studies designed to minimize potentially complicating variables. In the current study we explore the use of QTL and other statistical analyses applied to an F2 intercross population for examining the phenotypic and potential genetic relationships among multiple complex traits. Of particular interest is the potential for assessing whether the presence of pleiotropic genetic loci may account for phenotypic correlations among traits, by analyzing for colocalization of QTLs for different traits, and applying additional statistical tests. Although these analyses cannot prove the presence of pleiotropy, they can generate important testable hypotheses as to whether pleiotropy may be present. We have applied these concepts to the setting described above of determination of bone density in older animals with hyperlipidemia and atherosclerosis, using an intercross between strains C57BL/6J and DBA/2J in which an atherogenic diet was administered. We report here the phenotypic and genetic relationships observed among the above metabolic and morphometric parameters, atherosclerotic lesion formation, and various measures of bone mass and density, including biomechanical properties, in mature (16 mo old) F2 female mice with diet-induced atherosclerosis. Advanced statistical analyses are employed to examine these multiple traits and assess potential pleiotropic effects of genetic loci identified.
| MATERIALS AND METHODS |
|---|
|
|
|---|
21 days, female mice were selected and placed on a diet of rodent chow (Purina 5001) containing 17.5% protein, 11% fat, 0.8% calcium, and 0.5% phosphorus. At 12 mo of age, this was switched to an atherogenic high-fat, high-cholesterol diet for another 16 wk. This consisted of 75% chow supplemented with 7.5% cocoa butter, 2.5% dextrose, 1.625% each of sucrose and dextrin, 1.25% cholesterol, and 0.5% sodium cholate, with 0.73% calcium and 0.68% phosphorus (diet 90221; Harlan Teklad, Madison, WI). The mice were given free access to food and had a 12:12-h light-dark cycle maintained throughout. Body weight and length (nose to anus) were determined prior to death. At death, blood was collected by retro-orbital sampling into EDTA anticoagulant, spleens were removed for DNA isolation, and fat pads (omental, retroperitoneal, parametrial, and subcutaneous) were removed and weighed. Hearts and ascending aortas were removed and embedded in OCT compound and frozen at -70°C. Carcasses were placed in airtight individual bags and frozen at -20°C. Subsequently, the left femurs were dissected and stripped of soft tissue and placed directly in 10% buffered formalin. The right femurs were similarly obtained and were wrapped in gauze soaked in normal saline, placed in airtight plastic bags, frozen at -70°C, and used for biomechanical analyses.
Femoral bone analyses.
Bone mineral mass and indices for bone mineral density (BMD) were measured radiographically and by computerized tomography. For the former, all samples were collectively exposed on a single X-ray film using a Faxitron and subsequently developed under similar conditions. Separate X-rays were obtained for the anterior-posterior (AP) and lateral orientations. The X-rays were digitized and comparatively analyzed. Image 1/AT image processing software on a 486 IBM compatible computer and a CCD 72 video camera (Dage-MTI, Michigan City, IN) with Fujinon Zoom TV lens were used to digitize the X-rays. Image 1/AT was used to analyze the brightness and contrast level as well as the cortical diameter at a fixed axial location of the mouse femur diaphysis. Indices calculated from the digitized measurements for the AP as well as the lateral views were the average gray scale (integral GS) for the entire femur and for cross sections at fixed locations of the diaphysis and intertrochanteric regions [radiographic bone mineral content (BMC-R)]. Each of these was averaged by the total area or respective bone width (D). This method of normalizing for variation in bone size yielded the best correlation with computed tomography (CT) determined mineral content (r = 0.77) and density (r = 0.71) and is analogous to validated measures used for metacarpal bone density determination in human studies (1, 14, 16, 49). Values given are relative densitometry units. Variation in film and radiation exposure at various points on the film was determined to be less than 2%. The coefficients of variation for the primary measurements (integral GS, D) were less than 2%, based on repeated measures by a single observer on two separately prepared radiographs. Only results for the indices derived from the AP views is presented subsequently, as results for lateral views were strongly correlated with these and in almost all instances yielded results for phenotypic and genotypic analyses that were comparable.
CT determination of bone density was performed using electron beam computed tomography (EBCT). This technology has been used clinically for the past decade to detect and quantitate calcium deposits in coronary arteries and was adapted for determination of BMC in mice femurs. Femurs obtained as described above were embedded in a semisolid reusable gel phantom and scanned along with calibration bar phantoms containing known quantities of calcium hydroxyapatite in an Evolution EBCT scanner (Siemens Medical Systems, Iselin, NJ). Bones were oriented perpendicular to the scanning plane, and a total of 40 axial slices were obtained to image the entire phantom. Imaging employed a 1.5-mm collimation continuous volume scan mode, and images were reconstructed using vendor supplied sharp reconstruction kernel, a 12.7-cm field of view, and a 512 x 512 matrix. The calcium score was quantitated using a modification of the vendor supplied scoring software and a calcium threshold of 130 Hounsfield units (HU). Bone volume was determined by summation of the cross-sectional areas of each axial image multiplied by the 1.5-mm collimation thickness. BMD was assessed as the average HU value for each bone, and a measure of BMC was obtained by multiplying the average HU value by the bone volume. This method was demonstrated to have a highly significant correlation of the BMC determination with ash weight (r2 = 0.909), and a between-run coefficient of variation of 1.19%.
Torsion testing.
The proximal and distal mouse femora were potted into square molds using polymethylmethacrylate cement. Torsional testing was subsequently performed using a Burstein-Frankel torsion tester (8, 33). A deadweight pendulum was used to reach torsional failure in less than 0.1 s. The femora were tested in a random sequence and always in external rotation. A torque/rotation curve was recorded using a storage oscilloscope, and the results were photographed. The curve was analyzed using a computer program written in BASIC to measure the failure torque (N-m), angular deformation at failure (rad), stiffness (N-m/rad), and toughness (N-m·rad). These values were calculated using a linear regression model over the linear aspect of the curve. Reproducibility of the technique determined by testing of left and right femurs from 10 mice was 90%. Bones were excluded from the analysis if they did not demonstrate a clean spiral fracture indicating pure torque failure, or the angular displacement was excessively large with low torque indicating slippage. Satisfactory response curves were obtained on 117 of 140 attempted. Unsatisfactory test results were most commonly due to slippage of the ends in the potting material.
Plasma biochemical measurements and aortic fatty lesion quantitation.
Plasma total cholesterol, high-density and low-density lipoprotein (HDL and LDL) cholesterol, triglycerides, and free fatty acids were measured as previously described (32). Plasma leptin and insulin were measured by immunoassay according to manufacturers instructions, using kits specific for mouse proteins (R&D Systems, Minneapolis, MN; and Crystal Chem, Chicago, IL; respectively). Quantitation of aortic fatty lesion formation was performed as previously described. This entails histomorphometric analyses of Oil red-O stained serial frozen sections of the ascending aorta from the aortic valves through the aortic arch (42).
Linkage and data analyses.
A complete linkage map for all chromosomes except the Y chromosome was constructed at an average density of 13 cM using microsatellite markers as previously described (17). The constructed linkage map with all markers used is available on request. PCR primers for microsatellite markers were purchased from Research Genetics (Huntsville, AL); methods for use were as described. ANOVA, correlation, regression analyses, and hierarchical clustering were performed using StatView version 5.0 (SAS Institute, Berkeley, CA) and SPLUS 2000 (MathSoft). Linkage maps were constructed and QTL analysis was performed using MapMaker QTL (29) and QTL Cartographer (4). "Log of the odds ratio" (LOD) scores were calculated at 1-cM intervals throughout the genome for each of the traits under consideration. Before analysis, the normal distribution assumption underlying parametric tests for linkage was assessed for each trait by examining the histogram, as well as the skewness and kurtosis statistics for the trait values. In addition, despite its low power, the Kolmogorov-Smirnov test was applied to assist in determining whether a given trait was sampled from a normal distribution. Traits demonstrating obvious departures from normality were log-transformed or square root-transformed and then reassessed for normality. Linkage results for traits in which normality was questionable were verified using nonparametric association tests, such as the Wilcoxon rank sum test, in determining marker/trait associations. Trait values falling outside the mean ± 3 SD were considered outliers and excluded from analysis, accounting for the variation in numbers given in Table 1. A LOD score of 2.8 was the threshold for suggestive evidence of linkage, and a LOD score greater than 4.3 was interpreted as significant evidence for linkage, as proposed by Lander and Kruglyak for an F2 intercross with 2 degrees of freedom (25).
|
| RESULTS |
|---|
|
|
|---|
|
|
Analysis of plasma lipid correlations with bone measures showed that a prooxidant lipid profile was associated with reduced bone mass and density. HDL levels were positively correlated with bone measurements, whereas LDL levels were inversely correlated. Expressing these as a ratio of HDL to LDL increased the correlation coefficients, which were significant for almost all the bone measurements (r values for bone density measurements ranged from 0.2020.315). Aortic fatty lesions showed significant inverse correlations with bone mass and density. As plasma lipid levels are themselves correlated with body adiposity and morphometry, stepwise regression analyses were performed to determine whether these were independent factors influencing bone or were secondary. Results showed that among the body morphometric and lipid- and lesion-related parameters discussed, adiposity was the major determinant for all, with body length being the only other contributing independent variable for the bone mass parameters. These measured variables together accounted for
25% of the observed variability in the various bone mass and density parameters.
Identification of genetic loci controlling bone-related traits.
A complete linkage map for all the chromosomes except the Y chromosome was constructed using microsatellite markers, with an average interval of 13 cM. CIM was performed as implemented by QTL Cartographer software, using four background markers identified by forward/backward regression analysis, in addition to two flanking markers for each trait and locus. This analysis identified six loci showing linkage with multiple bone-related traits (Table 3). These six loci were further analyzed using the multi-trait program (JZMAPQTL) of QTL Cartographer, incorporating up to four bone-related traits showing linkage to the respective loci. At each locus, the LOD score obtained by the multi-trait analysis significantly exceeded that obtained by CIM for any one trait and were highly significant (LOD >5) for all six loci (Table 3; Fig. 2).
|
|
In addition to considering multiple traits simultaneously and detecting epistatic effects, it is of interest to determine the additive and dominance effects for a given QTL. The CIM model and joint CIM models provided in the QTL Cartographer software automatically generate test statistics for assessing additive and dominance effects. Table 3 summarizes the additive and dominance effects for each of the six loci.
Colocalization of QTL for bone and nonbone traits: potential pleiotropy vs. close linkage.
QTL analyses were also performed for the morphometric and metabolic parameters that showed phenotypic correlations with the bone-related traits. One or more QTL each were identified for adipose tissue traits, leptin, body weight and length, plasma lipids, and aortic lesions. Three of the six QTL for bone-related traits were found to be adjacent to or overlapping with QTL for nonbone traits (chromosome 2, adipose tissue and plasma HDL; chromosome 3, plasma LDL and body length; chromosome 6, adipose tissue and plasma HDL) (Figs. 35). Among bone and nonbone traits with adjacent or colocalizing QTL, allele effects were similar in each case.
|
|
|
![]() |
![]() |
2 distribution with 1 degree of freedom. The C++ software for testing these particular hypotheses is available upon request. Multi-trait analysis was performed on bone and nonbone traits for QTL reported on chromosomes 2, 3, and 6. Incorporation of the non-bone-related trait(s) led to significant increases in the LOD scores for the loci on chromosomes 2, 3, and 6 (Figs. 3 5). The pleiotropy vs. close linkage test was then applied to each of the specific genome regions defined by the QTL resulting from the multi-trait analyses. This test yielded the following results. For chromosome 2, BMD-R and BMD-Rd were tested against adipose mass and HDL cholesterol to determine whether the QTL on chromosome 2 potentially represented the same locus (pleiotropy) or distinct, closely linked loci. The pleiotropy test could not be rejected at the 0.05 significance level for any of the comparisons. For chromosome 3, the chromosome 3 QTL for BMC-CT was tested against the adjacent chromosome 3 QTL for body length and HDL cholesterol, to determine whether they potentially represented the same QTL or closely linked loci. The test was rejected at the P < 0.05 level, for each comparison, indicating that these loci were distinct. Finally, on chromosome 6, total radiographic intertrochanteric density (BMC-Rit) (at position 68 cM) was compared with plasma HDL (at position 73 cM) and subcutaneous fat pad mass (at position 65 cM), and the pleiotropy test could not be rejected at the 0.05 significance level. It is of note that the pleiotropy test could be rejected at the 0.01 significance level when adipose tissue mass (at position 50 cM) was compared against BMC-Rit (at position 68 cM), indicating there may be two QTL controlling for adipose tissue mass and subcutaneous fat pad mass on chromosome 6.
The specificity and sensitivity of the pleiotropy/close linkage test have not been thoroughly explored. However, the position of the QTLs tested affects the power of this test, since the test is dependent on observing recombination in the testing regions defined by the positions under consideration. Therefore, although insufficient support to reject this test is consistent with pleiotropy, it does not prove conclusively that pleiotropic effects of a single gene are controlling for the multiple traits, since the power to detect multiple QTLs at close distances is greatly reduced and since these small distances may contain a large number of genes.
| DISCUSSION |
|---|
|
|
|---|
Adipose tissue and plasma HDL traits were correlated with bone mass and density as well as biomechanical properties, and loci on distal chromosome 2 and distal chromosome 6 showed colocalization of QTLs for bone traits with both of these. For the chromosome 2 and 6 QTLs, incorporating adipose tissue and plasma HDL traits into the multi-trait QTL analyses significantly raised the strength of linkage, which is consistent with pleiotropy. Additionally, the test for pleiotropy vs. close linkage supported potential pleiotropy for the chromosome 2 and the distal chromosome 6 loci. The chromosome 2 locus falls within a region identified in other studies as showing linkage with obesity (32).
Clinical studies have suggested that obesity confers a protective effect against osteoporosis and fractures that is independent of its contribution to body weight (43, 44, 46). Traditional explanations for this include the increased peripheral conversion of estrogen and increased insulin levels associated with adipose tissue, each having positive effects on bone mass and density (11, 45). However, recent studies of the role of leptin have shown that the situation is likely much more complicated. Clinical studies have demonstrated that plasma leptin levels are correlated with bone density (31), as we observed in this study, and in vitro leptin stimulates osteoblastic differentiation in progenitor stromal cells (52). It is intriguing that adipocytes and osteoblasts can be derived from common marrow stromal precursor cells and thus might share common cell biologic features (41). On the other hand, leptin is now well documented to be a central neuroendocrine mediator causing inhibition of bone formation, among other effects (13). Also, administration of a high-fat diet to mice can impair both central and peripheral sensitivity to leptin (28) and cause reduced bone density (40).
Similarly, there is observational and experimental evidence linking plasma lipids with bone density. Atherosclerotic lesions and the concomitant lower plasma HDL and higher plasma LDL cholesterol concentrations were associated with reduced bone density, a relationship observed in clinical studies in women (3, 19, 54). Although the relationship could simply be parallel changes occurring with age, recent in vitro and in vivo studies showing inhibitory effects of oxidized lipids on osteoblast differentiation and on bone density have provided evidence for a causative relationship between these (3840). Our findings of two loci showing linkage with bone, adipose tissue, and plasma lipid traits are intriguing in light of the above and allow us to hypothesize that genetic pleiotropy may be responsible. Further studies will be needed to definitively establish whether this is the case or simply close linkage of distinct genes is responsible for the observation.
Several recent studies have reported results of QTL analyses for peak bone density in mice, including one by Klein and colleagues (23) that utilized recombinant inbred (RI) strains derived from the same parental strains studied here. Before discussing our findings in relation to these, it should first be noted that the bone density traits examined here are not directly comparable. Rather than reflecting peak bone density, they purposely reflect bone density in the more complicated setting of older age and the concurrent administration of a highly atherogenic diet. Both of these factors could be responsible for some or all of the QTL identified. However, since the traits reflect the cumulative genetic effects over the life of the mice, the QTL may also reflect those affecting peak bone density. In light of the latter, it is of interest that the distal chromosome 2 locus was linked with bone density in both our study and that of Klein et al. (23), suggesting that it is indeed an important locus. Although the same genetic variations should be present in the RI strains and in the F2 intercross mice (since the parental strains are similar), other loci were not found in common, which may be explained by differing study conditions as discussed above and, potentially, the different techniques for measuring bone density. Also, in general, studies with RI strains also have less statistical power than those using intercrosses and are more susceptible to false-positive results. Loci for peak bone density in an intercross between strains C3H/HeJ and C57BL/6J, C57BL/6J and CAST/Ei, and in the osteoporotic SAM-P6 mouse have also been reported recently and do not correspond to loci observed in our study (5, 7, 51). In these cases, strain differences would be expected to yield different loci, in addition to the age and diet effects.
In humans, variations in the genotype for the vitamin D receptor, estrogen receptor, collagen type 1a, and interleukin 6 have been associated with differences in bone density, as has a QTL on chromosome 11, but these are not mapped to any of the loci identified in this study (15, 20, 24, 34, 53). It is important to note that, in studies such as this, lack of linkage is simply noninformative, as there may be no biologically relevant genetic variation between the two strains used. Also, as noted above, the QTL we have detected may be related to age or diet effects rather than peak bone mass. The next steps in the identification of the genes responsible at each of the loci identified here will require finer mapping through the use of strains congenic for each of the loci. Continued identification and mapping of new genes through the various genome projects will also expand the pool of candidate genes to consider.
This study demonstrates the value of combining phenotypic and genetic analysis for elucidating the pathogenic interactions among multiple related traits. Over the past 5 years, a variety of statistical techniques have been developed to resolve coincident QTL, to consider multiple traits simultaneously, and to test whether coincident QTL are closely linked or represent pleiotropic effects from a single QTL. As the number of traits being simultaneously analyzed increases, the power to detect marker/trait associations can potentially dramatically increase if appropriate statistical models are used. The results obtained in this study demonstrate this dramatic increase in power, and the pleiotropy test developed by Jiang and Zeng (18) further demonstrates that an application of more sophisticated statistical tests can go far in elucidating the underlying genetic factors contributing to multiple quantitative traits. We were able to increase the evidence for linkage, given by the LOD score, for six QTL by more than one order of magnitude, on average, and by more than three orders of magnitude in two cases. In addition, we were able to reject the pleiotropy hypothesis for the chromosome 3 loci, but were unable to reject this hypothesis for the chromosome 2 and the distal chromosome 6 QTL, by implementing and then applying the pleiotropy vs. close linkage statistical test. As more phenotype data obtain from a variety of sources (including microarray-measured transcript levels in addition to the more classic traits measured here) there will be increased demand for more sophisticated statistical methods capable of deriving more information from the data, such as those we have employed here. In our own study, we have tried to maximize the available information in our data by making use of the most advanced statistical techniques available to date.
| ACKNOWLEDGMENTS |
|---|
This work was supported in part by National Heart, Lung, and Blood Institute Grants HL-30568 (to A. J. Lusis and T. A. Drake) and HL-28481 (to A. J. Lusis) and by the Laubisch Fund, UCLA. K. Krass was supported by National Research Service Award GM-07014.
| FOOTNOTES |
|---|
Address for reprint requests and other correspondence: T. A. Drake, Dept. of Pathology and Laboratory Medicine, UCLA, Los Angeles, CA 90095-1732 (E-mail: tdrake{at}mednet.ucla.edu).
| REFERENCES |
|---|
|
|
|---|
1 gene to bone density and the risk of osteoporotic fractures in postmenopausal women. N Engl J Med 338: 10161021, 1998.This article has been cited by other articles:
![]() |
B. Lecka-Czernik, C. Ackert-Bicknell, M. L. Adamo, V. Marmolejos, G. A. Churchill, K. R. Shockley, I. R. Reid, A. Grey, and C. J. Rosen Activation of Peroxisome Proliferator-Activated Receptor {gamma} (PPAR{gamma}) by Rosiglitazone Suppresses Components of the Insulin-Like Growth Factor Regulatory System in Vitro and in Vivo Endocrinology, February 1, 2007; 148(2): 903 - 911. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. M. Delahunty, K. L. Shultz, G. A. Gronowicz, B. Koczon-Jaremko, M. L. Adamo, L. G. Horton, J. Lorenzo, L. R. Donahue, C. Ackert-Bicknell, B. E. Kream, et al. Congenic Mice Provide in Vivo Evidence for a Genetic Locus that Modulates Serum Insulin-Like Growth Factor-I and Bone Acquisition Endocrinology, August 1, 2006; 147(8): 3915 - 3923. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Doss, E. E. Schadt, T. A. Drake, and A. J. Lusis Cis-acting expression quantitative trait loci in mice Genome Res., May 1, 2005; 15(5): 681 - 691. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. Nishimura, T. A. Drake, A. J. Lusis, K. M. Lyons, J. H. Nadeau, and J. Zernik ENU LARGE-SCALE MUTAGENESIS AND QUANTITATIVE TRAIT LINKAGE (QTL) ANALYSIS IN MICE: NOVEL TECHNOLOGIES FOR SEARCHING POLYGENETIC DETERMINANTS OF CRANIOFACIAL ABNORMALITIES Crit. Rev. Oral. Biol. Med., September 1, 2003; 14(5): 320 - 330. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Allayee, A. Ghazalpour, and A. J. Lusis Using Mice to Dissect Genetic Factors in Atherosclerosis Arterioscler. Thromb. Vasc. Biol., September 1, 2003; 23(9): 1501 - 1509. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. A. Lyons, H. Wittenburg, R. Li, K. A. Walsh, G. A. Churchill, M. C. Carey, and B. Paigen Quantitative trait loci that determine lipoprotein cholesterol levels in DBA/2J and CAST/Ei inbred mice, J. Lipid Res., May 1, 2003; 44(5): 953 - 967. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. L. Svenson, M. A. Bogue, and L. L. Peters Genetic Models in Applied Physiology: Invited Review: Identifying new mouse models of cardiovascular disease: a review of high-throughput screens of mutagenized and inbred strains J Appl Physiol, April 1, 2003; 94(4): 1650 - 1659. [Abstract] [Full Text] [PDF] |
||||
![]() |
H.-W. Deng, F.-H. Xu, Q.-Y. Huang, H. Shen, H. Deng, T. Conway, Y.-J. Liu, Y.-Z. Liu, J.-L. Li, H.-T. Zhang, et al. A Whole-Genome Linkage Scan Suggests Several Genomic Regions Potentially Containing Quantitative Trait Loci for Osteoporosis J. Clin. Endocrinol. Metab., November 1, 2002; 87(11): 5151 - 5159. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||