Confirmation and fine mapping of chromosomal regions influencing peak bone mass in mice

Genetic effects for femoral biomechanics, structure, and density in C57BL/6J and C3H/HeJ inbred mouse strains

Genome-wide QTL analysis for bone density, structure, and biomechanical phenotypes was performed in 999 (B6xC3H)F2 mice. Multivariate phenotypes were also derived to test for pleiotropic QTL effects. Highly significant QTLs were detected with pleiotropic effects on many of these phenotypes, and QTLs with unique effects on specific phenotypes were found as well.

INTRODUCTION

The inbred C57BL/6J (B6) and C3H/HeJ (C3H) mouse strains were previously shown to segregate quantitative trait loci (QTLs) for femoral bone density.

MATERIALS AND METHODS

The 999 s filial (F2) mouse progeny were further phenotyped for measures of femoral biomechanics (load to failure, Fu; work to failure, U; stiffness, S), structure (polar moment of inertia, Ip; moment of inertia ratio, Ir), and more specific femoral midshaft bone density measures (cortical and total vBMD). Two novel multivariate phenotypes were computed using principal component analysis, thus aiding in the exploration of pleiotropic effects of the QTLs detected.

RESULTS AND CONCLUSIONS

Results of a genome-wide analysis provided strong evidence of pleiotropic QTL effects on chromosome 4, with six of the seven primary phenotypic measures, representing femoral biomechanics, density, and structure, producing LOD scores greater than 8. Chromosomes 1, 8, 13, and 14 were also identified as harboring QTLs that affect phenotypes in two of the three aspects of bone properties. QTLs uniquely contributing to variability in biomechanical measures were identified on chromosomes 10 and 12, whereas a QTL solely affecting structure was found on chromosome 17. Analysis of the evidence for pleiotropic effects using principal component analysis revealed pleiotropic QTLs on chromosomes 4 and 14, influencing nearly all the bone phenotypes measured and revealed QTLs on chromosomes 1, 8, 13, and 17 with pleiotropic effects restricted to either density or the structure and stiffness phenotypes. The use of multivariate phenotypes has allowed us to identify pleiotropic effects of several QTLs previously linked in studies of other mouse strains and in human studies of bone mineral density and femoral structure, which will provide important insight regarding the importance of allelic variation on the entire skeleton.

Quantitative trait loci for femoral size and shape in a genetically heterogeneous mouse population

The aim of this study was to examine the genetic effects on cortical bone geometry. Genotypes from 487 mice were compared with geometric traits obtained from microCT. We found 14 genetic markers that associate with geometric traits, showing the complexity of genetic control over bone geometry.

INTRODUCTION

Previous studies have shown that genetic background affects bone characteristics, particularly bone mineral density, in both mouse and human populations. Much less is known, however, about the effects of polymorphic genes on bone size, shape, and mechanical integrity. In this study, we investigated the genetic determinants of geometric properties of cortical bone in mice.

MATERIALS AND METHODS

This study used a genetically heterogeneous mouse population, which is denoted UM-HET3 stock and is derived as the progeny of (BALB/cJ X C57BL/6J) F1 females and (C3H/HeJ X DBA/2J) F1 males. The experimental group consisted of 487 female UM-HET3 mice. Genotypic data from 99 polymorphic genetic loci was obtained from the mice at 4 weeks of age. At 18 months of age, the mice were humanely killed, and the right femurs were scanned with microcomputed tomography to assess geometric properties of cortical bone. A permutation-based test was used to detect significant associations between genetic markers and geometric traits. This test generates experiment-wise p values, which account for the effect of testing multiple hypotheses. An experiment-wise p < or = 0.05 was considered statistically significant.

RESULTS

Fourteen genetic markers were found to significantly associate with one or more geometric traits. Two markers (D3Mit62 and D4Mit155) were associated with traits describing bone size; 2 (D12Mit167 and D14Mit170) were linked with traits describing bone shape; and 10 (D1Nds2, D5Mit95, D6Mit216, D7Mit91, D8Mit51, D9Mit110, D11Mit83, D15Mit100, D15Mit171, and D17Mit46) were associated with both size and shape.

CONCLUSIONS

Our results indicate that the genetic control of cortical bone geometry is complex and that femoral size and shape may be influenced by different, although overlapping, groups of polymorphic loci.

Precision and accuracy of peripheral quantitative computed tomography (pQCT) in the mouse skeleton compared with histology and microcomputed tomography (microCT)

pQCT was evaluated for accuracy of phenotypic characterization of mouse bone in vivo. Bones (tibia, femur, spine) of 27 animals were measured ex vivo with pQCT, microCT, and histomorphometry and of 23 mice in vivo (pQCT). pQCT yielded satisfactory in vivo precision and accuracy in skeletal characterization.

INTRODUCTION

Important aspects of modern skeletal research depend on the phenotypic characterization of genetically manipulated mice, with some approaches requiring in vivo measurement. Peripheral quantitative computed tomography (pQCT) is applicable in vivo and provides opportunities to determine a large variety of bone parameters. Here we test the ex vivo and in vivo reproducibility of pQCT, and its accuracy in comparison with histomorphometry and microcomputed tomography (microCT).

MATERIALS AND METHODS

We examined the tibia, femur, and lumbar spine of 27 mice ex vivo with high-resolution pQCT, using two mouse models (wild-type and ob/ob) with known differences in bone density. Measurements were repeated three times at different days in nine animals. In a second experiment, 23 animals (10 wild-type and 13 bGH transgenic mice) were repeatedly measured in vivo at 12 and 13 weeks of age, respectively.

RESULTS

Among metaphyseal sites, the ex vivo precision was highest at the distal femur (RMS CV < 1% for density and < 2% for area). The correlation between density (pQCT) and bone volume fraction (histomorphometry) was r2 = 0.79 (tibia, femur, and spine), and that with microCT was r2 = 0.94 (femur). At the diaphysis, the precision was highest at the femur (< 2% for total and cortical area), and the correlation with microCT was r2 > 0.77. The in vivo precision for bone density (distal femur) was 2.3-5.1%, and that for absolute and relative cortical area (tibia) was 3.1% and 2.2%.

CONCLUSIONS

The results show that pQCT can yield satisfactory precision and accuracy in skeletal characterization of mouse bones, if properly applied. The potential advantage of pQCT is that it provides a large set of parameters on bone properties and that it can be used in vivo, extending the available methodological repertoire for genetic studies.

High bone mass in mice expressing a mutant LRP5 gene

A unique mutation in LRP5 is associated with high bone mass in man. Transgenic mice expressing this LRP5 mutation have a similar phenotype with high bone mass and enhanced strength. These results underscore the importance of LRP5 in skeletal regulation and suggest targets for therapies for bone disease. A mutation (G171V) in the low-density lipoprotein receptor related protein 5 (LRP5) has been associated with high bone mass (HBM) in two independent human kindreds. To validate the role of the mutation, several lines of transgenic mice were created expressing either the human LRP5 G171V substitution or the wildtype LRP5 gene in bone. Volumetric bone mineral density (vBMD) analysis by pQCT showed dramatic increases in both total vBMD (30-55%) and trabecular vBMD (103-250%) of the distal femoral metaphysis and increased cortical size of the femoral diaphysis in mutant G171V transgenics at 5, 9, 17, 26, and 52 weeks of age (p < 0.01 for all). In addition, high-resolution microcomputed tomography (microCT) analysis of the distal femorae and lumbar vertebrae revealed an increase (110-232%) in trabecular bone volume fraction caused by both increased trabecular number (41-74%) and increased trabecular thickness (34-46%; p < 0.01 for all) in the mutant G171V mice. The increased bone mass was associated with significant increases in vertebral compressive strength (80-140%) and the increased cortical size with significant increases in femoral bending strength (50-130%). There were no differences in osteoclast number at 17 weeks of age. However, compared with littermate controls, the mutant G171V transgenic mice showed an increase in actively mineralizing bone surface, enhanced alkaline phosphatase staining in osteoblasts, and a significant reduction in the number of TUNEL-positive osteoblasts and osteocytes. These results suggest that the increased bone mineral density in mutant G171V mice was caused by increased numbers of active osteoblasts, which could in part be because of their increased functional lifespan. While slight bone anabolic activity was observed from overexpression of the wildtype LRP5 gene, it is clear that the G171V mutation, rather than overexpression of the receptor itself, is primarily responsible for the dramatic HBM bone effects. Together, these findings establish the importance of this novel and unexpected role of a lipoprotein receptor in regulating bone mass and afford a new model to explore LRP5 and its recent association with Wnt signaling in bone biology.

Congenic strains of mice for verification and genetic decomposition of quantitative trait loci for femoral bone mineral density

Peak femoral volumetric bone mineral density (femoral bone mineral density) in C57BL/6J (B6) 4-month-old female mice is 50% lower than in C3H/HeJ (C3H) and 34% lower than in CAST/EiJ (CAST) females. Genome-wide analyses of (B6 x C3H)F2 and (B6 x CAST)F2 4-month-old female progeny demonstrated that peak femoral bone mineral density is a complex quantitative trait associated with genetic loci (QTL) on numerous chromosomes (Chrs) and with trait heritabilities of 83% (C3H) and 57% (CAST). To test the effect of each QTL on femoral bone mineral density, two sets of loci (six each from C3H and CAST) were selected to make congenic strains by repeated backcrossing of donor mice carrying a given QTL-containing chromosomal region to recipient mice of the B6 progenitor strain. At the N6F1 generation, each B6.C3H and B6.CAST congenic strain (statistically 98% B6-like in genomic composition) was intercrossed to obtain N6F2 progeny for testing the effect of each QTL on femoral bone mineral density. In addition, the femoral bone mineral density QTL region on Chr 1 of C3H was selected for congenic subline development to facilitate fine mapping of this strong femoral bone mineral density locus. In 11 of 12 congenic strains, 6 B6.C3H and 5 B6.CAST, femoral bone mineral density in mice carrying c3h or cast alleles in the QTL regions was significantly different from that of littermates carrying b6 alleles. Differences also were observed in body weight, femoral length, and mid-diaphyseal periosteal circumference among these 11 congenic strains when compared with control littermates; however, these latter three phenotypes were not consistently correlated with femoral bone mineral density. Analyses of eight sublines derived from the B6.C3H-1T congenic region revealed two QTLs: one located between 36.9 and 49.7 centiMorgans (cM) and the other located between 73.2 and 100.0 cM distal to the centromere. In conclusion, these congenic strains provide proof of principle that many QTLs identified in the F2 analyses for femoral bone mineral density exert independent effects when transferred and expressed in a common genetic background. Furthermore, significant differences in femoral bone mineral density among the congenic strains were not consistently accompanied by changes in body weight, femur length, or periosteal circumference. Finally, decomposition of QTL regions by congenic sublines can reveal additional loci for phenotypes assigned to a QTL region and can markedly refine genomic locations of quantitative trait loci, providing the opportunity for candidate gene testing.

Evidence for a major gene underlying bone size variation in the Chinese

Osteoporosis is a major public health problem defined as a loss of bone strength, of which bone size is an important determinant. In the present study, familial correlation and segregation analyses for the spine and hip bone sizes were performed for the first time in a Chinese sample composed of 393 nuclear families with a total of 1,193 individuals. The results indicate a major gene of codominant inheritance for spine bone size; however, there is no evidence of a major gene influencing hip bone size. Significant familial residual effects are found for both traits, suggesting their polygenic inheritance. Heritability estimates (+/-SE) for spine and hip bone size were 0.62 (0.13) and 0.59 (0.12), respectively. Sex and age differences in genotype-specific average bone size were observed. Compared with our previous study on bone mineral density (BMD) in the same population, this study suggests that genetic determination of bone size may be different from that of BMD, and thus studying bone size as one surrogate phenotype for osteoporotic fractures may be necessary.

Mapping quantitative trait loci that influence femoral cross-sectional area in mice

Size and shape are critical determinants of the mechanical properties of skeletal elements and can be anticipated to be highly heritable. Moreover, the genes responsible may be independent of those that regulate bone mineral density (BMD). To begin to identify the heritable determinants of skeletal geometry, we have examined femoral cross-sectional area (FCSA) in male and female mice from two inbred strains of mice with divergent FCSA (C57BL/6 [B6] and DBA/2 [D2]), a large genetically heterogeneous population (n = 964) of B6D2F2 mice and 18 BXD recombinant inbred (RI) strains derived from their F2 cross. Femora were harvested from 16-week-old mice and FCSA (bone and marrow space enclosed within the periosteum) was measured at the midshaft by digital image analysis. In all mouse populations examined, FCSA was positively correlated with body weight and weight-corrected FCSA (WC-FCSA) values were normally distributed in the BXD-RI and F2 populations, suggesting polygenic control of this trait. Genome-wide quantitative trait locus (QTL) analysis of the B6D2F2 population revealed regions on four different chromosomes that were very strongly linked to WC-FCSA (chromosomes 6, 8, 10, and X) in both genders. Evidence of gender-specific genetic influences on femoral geometry was also identified at three other chromosomal sites (chromosomes 2, 7, and 12). Supporting evidence for the WC-FCSA QTLs on chromosomes 2, 7, 8, 10, and 12 also was present in the RI strains. Interestingly, none of these WC-FCSA QTLs were identified in our previous QTL analysis of whole body BMD in the same B6D2F2 population. Thus, the genetic determinants of bone size appear to be largely, if not entirely, distinct from those that regulate BMD attainment. The identification of the genes responsible for geometric differences in bone development should reveal fundamentally important processes in the control of skeletal integrity.

Gender specificity in the genetic determinants of peak bone mass

Peak bone mass is a major determinant of osteoporotic fracture risk. Gender differences in peak bone mass acquisition are well recognized in humans and may account for a substantial share of the increased prevalence of fragility fractures in women compared with men. Skeletal development is regulated by both heritable and environmental factors. Experimental animal models provide a means to circumvent complicating environmental factors. In this study we examined the heritability of peak bone mineral density (BMD) in genetically distinct laboratory mouse strains raised under strict environmental control and sought to identify genetic loci that may contribute to gender differences in this skeletal phenotype. Peak whole body BMD of male and female mice from a panel of 18 recombinant inbred (RI) strains derived from a cross between C57BL/6 and DBA/2 progenitors (BXD) was measured by dual-energy X-ray absorptiometry (DXA). A highly significant relationship existed between body weight and BMD in the BXD RI mice (r2 = 0.25; p = 1 x 10(-43)). To allow for comparison between male and female RI strains, whole body BMD values were corrected for the influence of body weight. The distribution of weight-corrected BMD (WC-BMD) values among the strains indicated the presence of strong genetic influences in both genders, with an estimated narrow sense heritability of 45% and 22% in male and female mice, respectively. Comparison of RI strain results by two-way analysis of variance (ANOVA) revealed a significant strain-by-gender interaction (F1,17,479 = 6.13; p < 0.0001). Quantitative trait locus (QTL) analysis of the BXD RI strain series provisionally identified nine chromosomal sites linked to peak bone mass development in males and seven regions in females. In two cases, the provisional chromosomal loci were shared between genders, but in most cases they were distinct (five female-specific QTLs and six male-specific QTLs). QTL analysis of a genetically heterogeneous F2 population derived from the B6 and D2 progenitor strains provided additional support for the gender specificity of two loci. A significant phenotype-genotype correlation was only observed in male F2 mice at microsatellite marker D7Mit114 on chromosome 7, and a correlation at D2Mit94 on chromosome 2 was only observed in female F2 mice. The present data highlight the important role of gender in the genetic basis of peak bone mass in laboratory mice. Because the male phenotype is associated with considerable fracture risk reduction, an elucidation of the nature of that effect could provide the basis for novel diagnostic, preventative, or therapeutic approaches.