Genetic dissection of mouse distal chromosome 1 reveals three linked BMD QTLs with sex-dependent regulation of bone phenotypes

The role of Zfp467 in mediating the pro-osteogenic and anti-adipogenic effects on bone and bone marrow niche

Conditional deletion of the PTH receptor (Pth1r) in mesenchymal progenitors reduces osteoblast differentiation and bone mass while enhancing adipogenesis and bone marrow adipose tissue. Mechanistically, PTH suppresses the expression of Zfp467, a pro-adipogenic zinc finger transcription factor. Consequently, Pth1r deficiency in mesenchymal progenitors leads to increased Zfp467 expression. Based on these observations, we hypothesized that genetic loss of Zfp467 would lead to a shift in marrow progenitor cell fate towards osteogenesis and increased bone mass. To test this hypothesis, we generated Zfp467-/- mice. Zfp467-/- mice (-/-) were significantly smaller than Zfp467+/+ mice (+/+). μCT showed significantly higher trabecular bone and cortical bone area in -/- vs. +/+, and histomorphometry showed higher structural and dynamic formation parameters in -/- mice vs. +/+. Femoral gene expression including Alpl, Sp7, and Acp5 were increased in -/-mice, whereas Adiponectin, Cebpa, Lepr, and Ppraγ mRNA were lower in -/- mice. Similarly, Fabp4 and Lep in the inguinal depot were also decreased in -/- mice. Moreover, marrow adipocyte numbers were reduced in -/- vs +/+ mice (p<0.007). In vitro, COBs and BMSCs-/- showed more positive ALP and Alizarin Red staining and a decrease in ORO droplets. Pth1r mRNA and protein levels were increased in COBs and BMSCs from -/- mice vs +/+ (p<0.02 for each parameter, -/- vs. +/+). -/- cells also exhibited enhanced endogenous levels of cAMP vs. control cells. Moreover, in an ovariectomy (OVX) mouse model, Zfp467-/- mice had significantly lower fat mass but similar bone mass compared to OVX +/+ mice. In contrast, in a high fat diet (HFD) mouse model, in addition to reduced adipocyte volume and adipogenesis related gene expression in both peripheral and bone marrow fat tissue, greater osteoblast number and higher osteogenesis related gene expression were also observed in -/- HFD mice vs. +/+ HFD mice. Taken together, these results demonstrate that ZFP467 negatively influences skeletal homeostasis and favors adipogenesis. Global deletion of Zfp467 increases PTHR1, cAMP and bone turnover, hence its repression is a component of PTH signaling and its regulation. These data support a critical role for Zfp467 in early lineage allocation and provide a novel potential mechanism by which PTH acts in an anabolic manner on the bone remodeling unit.

Mouse mutant phenotyping at scale reveals novel genes controlling bone mineral density

The genetic landscape of diseases associated with changes in bone mineral density (BMD), such as osteoporosis, is only partially understood. Here, we explored data from 3,823 mutant mouse strains for BMD, a measure that is frequently altered in a range of bone pathologies, including osteoporosis. A total of 200 genes were found to significantly affect BMD. This pool of BMD genes comprised 141 genes with previously unknown functions in bone biology and was complementary to pools derived from recent human studies. Nineteen of the 141 genes also caused skeletal abnormalities. Examination of the BMD genes in osteoclasts and osteoblasts underscored BMD pathways, including vesicle transport, in these cells and together with in silico bone turnover studies resulted in the prioritization of candidate genes for further investigation. Overall, the results add novel pathophysiological and molecular insight into bone health and disease.

Patients affected by osteoporosis frequently present with decreased BMD and increased fracture risk. Genes are known to control the onset and progression of bone diseases such as osteoporosis. Therefore, we aimed to identify osteoporosis-related genes using BMD measures obtained from a large pool of mutant mice genetically modified for deletion of individual genes (knockout mice). In a collaborative endeavor involving several research sites world-wide, we generated and phenotyped 3,823 knockout mice and identified 200 genes which regulated BMD. Of the 200 BMD genes, 141 genes were previously not known to affect BMD. The discovery and study of novel BMD genes will help to better understand the causes and therapeutic options for patients with low BMD. In the long run, this will improve the clinical management of osteoporosis.

Interleukin-6 Knockout Inhibits Senescence of Bone Mesenchymal Stem Cells in High-Fat Diet-Induced Bone Loss

Obesity, a chronic low-grade inflammatory state, not only promotes bone loss, but also accelerates cell senescence. However, little is known about the mechanisms that link obesity, bone loss, and cell senescence. Interleukin-6 (IL-6), a pivotal inflammatory mediator increased during obesity, is a candidate for promoting cell senescence and an important part of senescence-associated secretory phenotype (SASP). Here, wild type (WT) and (IL-6 KO) mice were fed with high-fat diet (HFD) for 12 weeks. The results showed IL-6 KO mice gain less weight on HFD than WT mice. HFD induced trabecular bone loss, enhanced expansion of bone marrow adipose tissue (BMAT), increased adipogenesis in bone marrow (BM), and reduced the bone formation in WT mice, but it failed to do so in IL-6 KO mice. Furthermore, IL-6 KO inhibited HFD-induced clone formation of bone marrow cells (BMCs), and expression of senescence markers (p53 and p21). IL-6 antibody inhibited the activation of STAT3 and the senescence of bone mesenchymal stem cells (BMSCs) from WT mice in vitro, while rescued IL-6 induced senescence of BMSCs from IL-6 KO mice through the STAT3/p53/p21 pathway. In summary, our data demonstrated that IL-6 KO may maintain the balance between osteogenesis and adipogenesis in BM, and restrain senescence of BMSCs in HFD-induced bone loss.

Genetic Dissection of Femoral and Tibial Microarchitecture

Our understanding of the genetic control of bone strength has relied mainly on estimates of bone mineral density. Here we have mapped genetic factors that influence femoral and tibial microarchitecture using high‐resolution x‐ray computed tomography (8‐μm isotropic voxels) across a family of 61 BXD strains of mice, roughly 10 isogenic cases per strain and balanced by sex. We computed heritabilities for 25 cortical and trabecular traits. Males and females have well‐matched heritabilities, ranging from 0.25 to 0.75. We mapped 16 genetic loci most of which were detected only in females. There is also a bias in favor of loci that control cortical rather than trabecular bone. To evaluate candidate genes, we combined well‐established gene ontologies with bone transcriptome data to compute bone‐enrichment scores for all protein‐coding genes. We aligned candidates with those of human genome‐wide association studies. A subset of 50 strong candidates fell into three categories: (1) experimentally validated genes already known to modulate bone function (Adamts4, Ddr2, Darc, Adam12, Fkbp10, E2f6, Adam17, Grem2, Ifi204); (2) candidates without any experimentally validated function in bone (eg, Greb1, Ifi202b), but linked to skeletal phenotypes in human cohorts; and (3) candidates that have high bone‐enrichment scores, but for which there is not yet any functional link to bone biology or skeletal system disease (including Ifi202b, Ly9, Ifi205, Mgmt, F2rl1, Iqgap2). Our results highlight contrasting genetic architecture between sexes and among major bone compartments. The alignment of murine and human data facilitates function analysis and should prove of value for preclinical testing of molecular control of bone structure. © 2019 The Authors. JBMR Plus published by Wiley Periodicals, Inc. on behalf of American Society for Bone and Mineral Research.

Intestinal clock system regulates skeletal homeostasis

The circadian clock network is an evolutionarily conserved system involved in the regulation of metabolic homeostasis; however, its impacts on skeletal metabolism remain largely unknown. We herein demonstrated that the circadian clock network in the intestines plays pivotal roles in skeletal metabolism such that the lack of the Bmal1 gene in the intestines (Bmal1Int-/- mice) caused bone loss, with bone resorption being activated and bone formation suppressed. Mechanistically, Clock protein interaction with the vitamin D receptor (VDR) accelerated its binding to the VDR response element by enhancing histone acetylation in a circadian-dependent manner, and this was lost in Bmal1Int-/- mice because nuclear translocation of Clock required the presence of Bmal1. Accordingly, the rhythmic expression of VDR target genes involved in transcellular calcium (Ca) absorption was created, and this was not observed in Bmal1Int-/- mice. As a result, transcellular Ca absorption was impaired and bone resorption was activated in Bmal1Int-/- mice. Additionally, sympathetic tone, the activation of which suppresses bone formation, was elevated through afferent vagal nerves in Bmal1Int-/- mice, the blockade of which partially recovered bone loss by increasing bone formation and suppressing bone resorption in Bmal1Int-/- mice. These results demonstrate that the intestinal circadian system regulates skeletal bone homeostasis.

Weak Epistasis Generally Stabilizes Phenotypes in a Mouse Intercross

The extent and strength of epistasis is commonly unresolved in genetic studies, and observed epistasis is often difficult to interpret in terms of biological consequences or overall genetic architecture. We investigated the prevalence and consequences of epistasis by analyzing four body composition phenotypes—body weight, body fat percentage, femoral density, and femoral circumference—in a large F2 intercross of B6-lit/lit and C3.B6-lit/lit mice. We used Combined Analysis of Pleiotropy and Epistasis (CAPE) to examine interactions for the four phenotypes simultaneously, which revealed an extensive directed network of genetic loci interacting with each other, circulating IGF1, and sex to influence these phenotypes. The majority of epistatic interactions had small effects relative to additive effects of individual loci, and tended to stabilize phenotypes towards the mean of the population rather than extremes. Interactive effects of two alleles inherited from one parental strain commonly resulted in phenotypes closer to the population mean than the additive effects from the two loci, and often much closer to the mean than either single-locus model. Alternatively, combinations of alleles inherited from different parent strains contribute to more extreme phenotypes not observed in either parental strain. This class of phenotype-stabilizing interactions has effects that are close to additive and are thus difficult to detect except in very large intercrosses. Nevertheless, we found these interactions to be useful in generating hypotheses for functional relationships between genetic loci. Our findings suggest that while epistasis is often weak and unlikely to account for a large proportion of heritable variance, even small-effect genetic interactions can facilitate hypotheses of underlying biology in well-powered studies.

The role of statistical epistasis in the genetic architecture of complex traits has been of great interest to the genetics community since Fisher introduced the concept in 1918. However, assessing epistasis in human and model organism populations has been impeded by limited statistical power. To mitigate this limitation, we analyzed bone and body composition traits in an unusually large mouse intercross population of over 2000 mice, paired with a recently-developed computational approach that leverages information to detect interactions across multiple phenotypes. We discovered a large network of highly significant genetic interactions between variants that influence complex body composition traits. Although epistasis was abundant, the interaction network was dominated by epistasis that stabilizes phenotypes by reducing phenotypic deviation from the parent strains. Nevertheless, the observed network provides an overview of genetic architecture and specific hypotheses of how QTL combine to affect phenotypes. These findings suggest that epistatic effects are generally of lesser magnitude than main QTL effects, and therefore are unlikely to account for major components of variance, but also reinforce genetic interaction analysis as a potent tool for dissecting the biology of complex traits.

Studies of OC-STAMP in Osteoclast Fusion: A New Knockout Mouse Model, Rescue of Cell Fusion, and Transmembrane Topology

The fusion of monocyte/macrophage lineage cells into fully active, multinucleated, bone resorbing osteoclasts is a complex cell biological phenomenon that utilizes specialized proteins. OC-STAMP, a multi-pass transmembrane protein, has been shown to be required for pre-osteoclast fusion and for optimal bone resorption activity. A previously reported knockout mouse model had only mononuclear osteoclasts with markedly reduced resorption activity in vitro, but with paradoxically normal skeletal micro-CT parameters. To further explore this and related questions, we used mouse ES cells carrying a gene trap allele to generate a second OC-STAMP null mouse strain. Bone histology showed overall normal bone form with large numbers of TRAP-positive, mononuclear osteoclasts. Micro-CT parameters were not significantly different between knockout and wild type mice at 2 or 6 weeks old. At 6 weeks, metaphyseal TRAP-positive areas were lower and mean size of the areas were smaller in knockout femora, but bone turnover markers in serum were normal. Bone marrow mononuclear cells became TRAP-positive when cultured with CSF-1 and RANKL, but they did not fuse. Expression levels of other osteoclast markers, such as cathepsin K, carbonic anhydrase II, and NFATc1, were not significantly different compared to wild type. Actin rings were present, but small, and pit assays showed a 3.5-fold decrease in area resorbed. Restoring OC-STAMP in knockout cells by lentiviral transduction rescued fusion and resorption. N- and C-termini of OC-STAMP were intracellular, and a predicted glycosylation site was shown to be utilized and to lie on an extracellular loop. The site is conserved in all terrestrial vertebrates and appears to be required for protein stability, but not for fusion. Based on this and other results, we present a topological model of OC-STAMP as a 6-transmembrane domain protein. We also contrast the osteoclast-specific roles of OC- and DC-STAMP with more generalized cell fusion mechanisms.

Sex and genetic factors determine osteoblastic differentiation potential of murine bone marrow stromal cells

Sex and genetic factors determine skeletal mass, and we tested whether bone histomorphometric parameters were sexually dimorphic in femurs from 1 to 6 month old C57BL/6 mice. Trabecular bone volume declined more rapidly in female mice than in male littermates because of enhanced bone resorption. Although bone formation was not different between sexes, female mice exhibited a higher number of osteoblasts than male littermates, suggesting that osteoblasts from female mice may have a reduced ability to form bone. To determine the impact of sex on osteoblastogenesis, we investigated the potential for osteoblastic differentiation of bone marrow stromal cells from C57BL/6, Friend leukemia virus-B (FVB), C3H/HeJ and BALB/c mice of both sexes. Bone marrow stromal cells from female FVB, C57BL/6 and C3H/HeJ mice exhibited lower Alpl and Osteocalcin expression and alkaline phosphatase activity, and formed fewer mineralized nodules than cells from male littermates. Proliferative capacity was greater in cells from male than female C57BL/6, but not FVB, mice. Sorting of bone marrow stromal cells from mice expressing an α-Smooth muscle actin-green fluorescent protein transgene, revealed a higher yield of mesenchymal stem cells in cultures from male mice than in those from female littermates. Sex had a modest impact on osteoblastic differentiation of mesenchymal stem cells. To determine the influence of sex and genetic factors on osteoblast function, calvarial osteoblasts were harvested from C57BL/6, FVB, C3H/HeJ and BALB/c mice. Alpl expression and activity were lower in osteoblasts from C57BL/6 and C3H/HeJ, but not FVB or BALB/c, female mice than in cells from littermates. Sex had no effect on osteoclastogenesis of bone marrow cultures of C57BL/6 mice, but osteoblasts from female mice exhibited higher Rankl and lower Opg expression than cells from male littermates. In conclusion, osteoblastogenesis is sexually dimorphic and influenced by genetic factors.

Quantitative trait loci for bone mineral density and femoral morphology in an advanced intercross population of mice

Osteoporosis, characterized by low levels of bone mineral density (BMD), is a prevalent medical condition in humans. We investigated its genetic and environmental basis by searching for quantitative trait loci (QTLs) affecting six skeletal (including three BMD) traits in a G10 advanced intercross population produced from crosses of mice from the inbred strain C57BL/6J with mice from a strain selected for high voluntary wheel running. The mice in this population were fed either a high-fat or a matched control diet throughout the study, allowing us to test for QTL by diet interactions for the skeletal traits. Our genome scan uncovered a number of QTLs, the great majority of which were different from QTLs previously found for these same traits in an earlier (G4) generation of the same intercross. Further, the confidence intervals for the skeletal trait QTLs were reduced from an average of 18.5 Mb in the G4 population to an equivalent of about 9 Mb in the G10 population. We uncovered a total of 50 QTLs representing 32 separate genomic sites affecting these traits, with a distal region on chromosome 1 harboring several QTLs with large effects on the BMD traits. One QTL was located on chromosome 5 at 4.0 Mb with a confidence interval spanning from 4.0 to 4.6 Mb. Only three protein coding genes reside in this interval, and one of these, Cyp51, is an attractive candidate as others have shown that developing Cyp51 knockout embryos exhibit shortened and bowed limbs and synotosis of the femur and tibia. Several QTLs showed significant interactions with sex, although only two QTLs interacted with diet, both affecting only mice fed the high-fat diet.

Identification of gender-specific candidate genes that influence bone microarchitecture in chromosome 1

Studies on the identification of the genetic basis for sexual dimorphism in peak bone mass are obviously important for providing novel therapeutic approaches to prevent or treat metabolic bone diseases. Our goal in this study was to identify the bone microstructure that could lead to differences in volumetric bone mineral density (vBMD) and new candidate genes that regulate the gender effect on bone. We used a congenic line of mice that carry the BMD1-4 locus from CAST/EiJ (CAST) mice in a C57BL/6J (B6) background and show greater vBMD in female, but not male, congenics compared to age- and gender-matched B6 mice. To assess the vBMD variations between the two lines of mice, we performed μCT measurements and found no difference in cortical bone volume by tissue volume (BV/TV) between congenics and B6 mice. However, trabecular BV/TV was significantly greater in female, but not male, congenics compared to corresponding B6 mice, which was due to increased trabecular thickness but not reduced trabecular separation, suggesting that bone formation, but not bone resorption, is responsible for the trabecular bone phenotype observed in the female, but not male, congenics. To identify the gender candidate genes, we determined the polymorphisms between B6 and CAST within the BMD1-4 locus and performed gene expression profiling. We identified EF-hand calcium binding domain (Efcab2), consortin, connexin sorting protein (Cnst), and presenilin 2 (Psen2) as potential candidate genes that regulate bone mass by influencing trabecular thickness in a gender-specific manner.

Genetic determinants of trabecular and cortical volumetric bone mineral densities and bone microstructure

Most previous genetic epidemiology studies within the field of osteoporosis have focused on the genetics of the complex trait areal bone mineral density (aBMD), not being able to differentiate genetic determinants of cortical volumetric BMD (vBMD), trabecular vBMD, and bone microstructural traits. The objective of this study was to separately identify genetic determinants of these bone traits as analysed by peripheral quantitative computed tomography (pQCT). Separate GWA meta-analyses for cortical and trabecular vBMDs were performed. The cortical vBMD GWA meta-analysis (n = 5,878) followed by replication (n = 1,052) identified genetic variants in four separate loci reaching genome-wide significance (RANKL, rs1021188, p = 3.6×10⁻¹⁴; LOC285735, rs271170, p = 2.7×10⁻¹²; OPG, rs7839059, p = 1.2×10⁻¹⁰; and ESR1/C6orf97, rs6909279, p = 1.1×10⁻⁹). The trabecular vBMD GWA meta-analysis (n = 2,500) followed by replication (n = 1,022) identified one locus reaching genome-wide significance (FMN2/GREM2, rs9287237, p = 1.9×10⁻⁹). High-resolution pQCT analyses, giving information about bone microstructure, were available in a subset of the GOOD cohort (n = 729). rs1021188 was significantly associated with cortical porosity while rs9287237 was significantly associated with trabecular bone fraction. The genetic variant in the FMN2/GREM2 locus was associated with fracture risk in the MrOS Sweden cohort (HR per extra T allele 0.75, 95% confidence interval 0.60–0.93) and GREM2 expression in human osteoblasts. In conclusion, five genetic loci associated with trabecular or cortical vBMD were identified. Two of these (FMN2/GREM2 and LOC285735) are novel bone-related loci, while the other three have previously been reported to be associated with aBMD. The genetic variants associated with cortical and trabecular bone parameters differed, underscoring the complexity of the genetics of bone parameters. We propose that a genetic variant in the RANKL locus influences cortical vBMD, at least partly, via effects on cortical porosity, and that a genetic variant in the FMN2/GREM2 locus influences GREM2 expression in osteoblasts and thereby trabecular number and thickness as well as fracture risk.

Osteoporosis is a common highly heritable skeletal disease characterized by reduced bone mineral density (BMD) and deteriorated bone microstructure, resulting in an increased risk of fracture. Most previous genetic epidemiology studies have focused on the genetics of the complex trait BMD, not being able to separate genetic determinants of the trabecular and cortical bone compartments and bone microstructure. The trabecular and cortical BMDs can be analysed separately by computed tomography. Therefore, we performed separate genome-wide association studies for trabecular and cortical BMDs, demonstrating that the genetic determinants of cortical and trabecular BMDs differ. Genetic variants in the RANKL, LOC285735, OPG, and ESR1 loci were associated with cortical BMD, while a genetic variant in the FMN2/GREM2 locus was associated with trabecular BMD. Two of these are novel bone-related loci. Follow-up analyses of bone microstructure demonstrated that a genetic variant in the RANKL locus is associated with cortical porosity and that the FMN2/GREM2 locus is associated with trabecular number and thickness. We propose that a genetic variant in the RANKL locus influences cortical BMD via effects on cortical porosity, and that a genetic variant in the FMN2/GREM2 locus influences trabecular BMD and fracture risk via effects on both trabecular number and thickness.

Age and sex affect quantitative genetic parameters for dominance rank and aggression in free-living greylag geese

Knowledge of the genetic and environmental influences on a character is pivotal for understanding evolutionary changes in quantitative traits in natural populations. Dominance and aggression are ubiquitous traits that are selectively advantageous in many animal societies and have the potential to impact the evolutionary trajectory of animal populations. Here we provide age- and sex-specific estimates of additive genetic and environmental components of variance for dominance rank and aggression rate in a free-living, human-habituated bird population subject to natural selection. We use a long-term data set on individually marked greylag geese (Anser anser) and show that phenotypic variation in dominance-related behaviours contains significant additive genetic variance, parental effects and permanent environment effects. The relative importance of these variance components varied between age and sex classes, whereby the most pronounced differences concerned nongenetic components. In particular, parental effects were larger in juveniles of both sexes than in adults. In paired adults, the partner's identity had a larger influence on male dominance rank and aggression rate than in females. In sex- and age-specific estimates, heritabilities did not differ significantly between age and sex classes. Adult dominance rank was only weakly genetically correlated between the sexes, leading to considerably higher heritabilities in sex-specific estimates than across sexes. We discuss these patterns in relation to selection acting on dominance rank and aggression in different life history stages and sexes and suggest that different adaptive optima could be a mechanism for maintaining genetic variation in dominance-related traits in free-living animal populations.

Assessment of gene-by-sex interaction effect on bone mineral density

Sexual dimorphism in various bone phenotypes, including bone mineral density (BMD), is widely observed; however, the extent to which genes explain these sex differences is unclear. To identify variants with different effects by sex, we examined gene-by-sex autosomal interactions genome-wide, and performed expression quantitative trait loci (eQTL) analysis and bioinformatics network analysis. We conducted an autosomal genome-wide meta-analysis of gene-by-sex interaction on lumbar spine (LS) and femoral neck (FN) BMD in 25,353 individuals from 8 cohorts. In a second stage, we followed up the 12 top single-nucleotide polymorphisms (SNPs; p < 1 × 10(-5) ) in an additional set of 24,763 individuals. Gene-by-sex interaction and sex-specific effects were examined in these 12 SNPs. We detected one novel genome-wide significant interaction associated with LS-BMD at the Chr3p26.1-p25.1 locus, near the GRM7 gene (male effect = 0.02 and p = 3.0 × 10(-5) ; female effect = -0.007 and p = 3.3 × 10(-2) ), and 11 suggestive loci associated with either FN- or LS-BMD in discovery cohorts. However, there was no evidence for genome-wide significant (p < 5 × 10(-8) ) gene-by-sex interaction in the joint analysis of discovery and replication cohorts. Despite the large collaborative effort, no genome-wide significant evidence for gene-by-sex interaction was found to influence BMD variation in this screen of autosomal markers. If they exist, gene-by-sex interactions for BMD probably have weak effects, accounting for less than 0.08% of the variation in these traits per implicated SNP. © 2012 American Society for Bone and Mineral Research.

Bone morphogenetic protein receptor type Ia localization causes increased BMP2 signaling in mice exhibiting increased peak bone mass phenotype

Bone morphogenetic protein 2 (BMP2) is a growth factor that initiates osteoblast differentiation. Recent studies show that BMP2 signaling regulates bone mineral density (BMD). BMP2 interacts with BMP receptor type Ia (BMPRIa) and type II receptor leading to the activation of the Smad signaling pathway. BMPRIa must shuttle between distinct plasma membrane domains, enriched of Caveolin-1 alpha and Caveolin-1 beta isoforms, and receptor activation occurs in these domains. Yet it remains unknown whether the molecular mechanism that regulates BMP2 signaling is driving mineralization and BMD. Therefore, the B6.C3H-1-12 congenic mouse model with increased BMD and osteoblast mineralization was utilized in this study. Using the family image correlation spectroscopy, we determined if BMP2 led to a significant re-localization of BMPRIa to caveolae of the alpha/beta isoforms in bone marrow stromal cells (BMSCs) isolated from B6.C3H-1-12 mice compared to the C57BL/6J mice, which served as controls. The control, C57BL/6J mice, was selected due to only 4 Mb of chromosome 1 from the C3H/HeJ mouse was backcrossed to a C57BL/6J background. Using reporter gene assays, the B6.C3H-1-12 BMSCs responded to BMP2 with increased Smad activation. Furthermore, disrupting caveolae reduced the BMP2-induced Smad signaling in BMSCs isolated from B6.C3H-1-12 and C57BL/6J. This study suggests for the first time a regulatory mechanism of BMPRIa signaling at the plasma membrane of BMSCs that (i) associated with genetic differences in the distal Chromosome 1 segment carried by the B6.C3H-1-12 congenic and (ii) contributes to increase BMD of the B6.C3H-1-12 compared to the C57BL/6J control mice.

Weak genetic relationship between trabecular bone morphology and obesity in mice

Obesity, in addition to being associated with metabolic diseases, such as diabetes, has also been found to lower the risk of osteoporotic fractures. The relationship between obesity and bone trabecular structure is complex, involving responses to mechanical loading and the effects of adipocyte-derived hormones, both directly interacting with bone tissue and indirectly through central nervous system signaling. Here we examine the effects of sex, a high fat diet, and genetics on the trabecular density and structure of the lumbar and caudal vertebra and the proximal tibia along with body weight, fat pad weight, and serum leptin levels in a murine obesity model, the LGXSM recombinant inbred (RI) mouse strains. The sample included 481 mice from 16 RI strains. We found that vertebral trabecular density was higher in males while the females had higher tibial trabecular density. The high fat diet led to only slightly higher trabecular density in both sexes despite its extreme effects on obesity and serum leptin levels. Trait heritabilities are moderate to strong and genetic correlations among trabecular features are high. Most genetic variation contrasts strains with large numbers of thick, closely-spaced, highly interconnected, plate-like trabeculae with a high bone volume to total volume ratio against strains displaying small numbers of thin, widely-spaced, sparsely connected, rod-like trabeculae with a low bone volume to total volume ratio. Genetic correlations between trabecular and obesity-related traits were low and not statistically significant. We mapped trabecular properties to 20 genomic locations. Only one-quarter of these locations also had effects on obesity. In this population obesity has a relatively minor effect on trabecular bone morphology.

Contribution of myostatin gene polymorphisms to normal variation in lean mass, fat mass and peak BMD in Chinese male offspring

AIM

Myostatin gene is a member of the transforming growth factor-β (TGF-β) family that negatively regulates skeletal muscle growth. Genetic polymorphisms in Myostatin were found to be associated with the peak bone mineral density (BMD) in Chinese women. The purpose of this study was to investigate whether myostatin played a role in the normal variation in peak BMD, lean mass (LM), and fat mass (FM) of Chinese men.

METHODS

Four hundred male-offspring nuclear families of Chinese Han ethnic group were recruited. Anthropometric measurements, including the peak BMD, body LM and FM were measured using dual-energy X-ray absorptiometry (DXA). The single nucleotide polymorphisms (SNPs) studied were tag-SNPs selected by sequencing. Both rs2293284 and +2278GA were genotyped using TaqMan assay, and rs3791783 was genotyped with PCR-restriction fragment length polymorphism (RFLP) analysis. The associations of the SNPs with anthropometric variations were analyzed using the quantitative transmission disequilibrium test (QTDT).

RESULTS

Using QTDT to detect within-family associations, neither single SNP nor haplotype was found to be associated with peak BMD at any bone site. However, rs3791783 was found to be significantly associated with fat mass of the trunk (P<0.001). Moreover, for within-family associations, haplotypes AGG, AAA, and TGG were found to be significantly associated with the trunk fat mass (all P<0.001).

CONCLUSION

Our results suggest that genetic variation within myostatin may play a role in regulating the variation in fat mass in Chinese males. Additionally, the myostatin gene may be a candidate that determines body fat mass in Chinese men.

A high-fat diet induces bone loss in mice lacking the Alox5 gene

5-Lipoxygenase catalyzes leukotriene generation from arachidonic acid. The gene that encodes 5-lipoxygenase, Alox5, has been identified in genome-wide association and mouse Quantitative Trait Locus studies as a candidate gene for obesity and low bone mass. Thus, we tested the hypothesis that Alox5(-/-) mice would exhibit metabolic and skeletal changes when challenged by a high-fat diet (HFD). On a regular diet, Alox5(-/-) mice did not differ in total body weight, percent fat mass, or bone mineral density compared with wild-type (WT) controls (P < 0.05). However, when placed on a HFD, Alox5(-/-) gained more fat mass and lost greater areal bone mass vs. WT (P < 0.05). Microarchitectural analyses revealed that on a HFD, WT showed increases in cortical area (P < 0.01) and trabecular thickness (P < 0.01), whereas Alox5(-/-) showed no change in cortical parameters but a decrease in trabecular number (P < 0.05) and bone volume fraction compared with WT controls (P < 0.05). By histomorphometry, a HFD did not change bone formation rates of either strain but produced an increase in osteoclast number per bone perimeter in Alox5(-/-) mice (P < 0.03). In vitro, osteoclastogenesis of marrow stromal cells was enhanced in mutant but not WT mice fed a HFD. Gene expression for Rankl, Pparg, and Cox-2 was greater in the femur of Alox5(-/-) than WT mice on a HFD (P < 0.01), but these increases were suppressed in the Alox5(-/-) mice after 8 wk of treatment with celecoxib, a cyclooxygenase-2 inhibitor. In sum, there is a strong gene by environmental interaction for bone mass when mice lacking the Alox5 gene are fed a HFD.

Multiple quantitative trait loci for cortical and trabecular bone regulation map to mid-distal mouse chromosome 4 that shares linkage homology to human chromosome 1p36

The mid-distal region of mouse chromosome 4 (Chr 4) is homologous with human Chr 1p36. Previously, we reported that mouse Chr 4 carries a quantitative trait locus (QTL) with strong regulatory effect on volumetric bone mineral density (vBMD). The intent of this study is to utilize nested congenic strains to decompose the genetic complexity of this gene-rich region. Adult females and males from 18 nested congenic strains carrying discrete C3H sequences were phenotyped for femoral mineral and volume by pQCT and for trabecular bone volume (BV), tissue volume (TV), trabecular number (Trab.no), and trabecular thickness (Trab.thk) by MicroCT 40. Our data show that the mouse Chr 4 region consists of at least 10 regulatory QTL regions that affected either or both pQCT and MicroCT 40 phenotypes. The pQCT phenotypes were typically similar between sexes, whereas the MicroCT 40 phenotypes were divergent. Individual congenic strains contained one to seven QTL regions. These regions conferred large positive or negative effects in some congenic strains, depending on the particular bone phenotype. The QTL regions II to X are syntenic with human 1p36, containing from 1 to 102 known genes. We identified 13 candidate genes that can be linked to bone within these regions. Six of these genes were linked to osteoblasts, three linked to osteoclasts, and two linked to skeletal development. Three of these genes have been identified in Genome Wide Association Studies (GWAS) linked to 1p36. In region III, there is only one gene, Lck, which conferred negative pQCT and MicroCT 40 phenotypes in both sexes. This gene is important to development and functioning of T cells, has been associated with osteoclast activity, and represents a novel bone regulatory gene that merits further experimental evaluation. In summary, congenic strains are powerful tools for identifying regulatory regions that influence bone biology and offer models for testing hypotheses about gene-gene and gene-environment interactions that are not available to experimental work in humans.

Mitogen-activated protein kinase phosphatase 1 regulates bone mass, osteoblast gene expression, and responsiveness to parathyroid hormone

Parathyroid hormone (PTH) signaling via PTH 1 receptor (PTH1R) involves mitogen-activated protein kinase (MAPK) pathways. MAPK phosphatase 1 (MKP1) dephosphorylates and inactivates MAPKs in osteoblasts, the bone-forming cells. We previously showed that PTH1R activation in differentiated osteoblasts upregulates MKP1 and downregulates pERK1/2-MAPK and cyclin D1. In this study, we evaluated the skeletal phenotype of Mkp1 knockout (KO) mice and the effects of PTH in vivo and in vitro. Microcomputed tomography analysis of proximal tibiae and distal femora from 12-week-old Mkp1 KO female mice revealed osteopenic phenotype with significant reduction (8-46%) in bone parameters compared with wild-type (WT) controls. Histomorphometric analysis showed decreased trabecular bone area in KO females. Levels of serum osteocalcin (OCN) were lower and serum tartrate-resistant acid phosphatase 5b (TRAP5b) was higher in KO animals. Treatment of neonatal mice with hPTH (1-34) for 3 weeks showed attenuated anabolic responses in the distal femora of KO mice compared with WT mice. Primary osteoblasts derived from KO mice displayed delayed differentiation determined by alkaline phosphatase activity, and reduced expressions of Ocn and Runx2 genes associated with osteoblast maturation and function. Cells from KO females exhibited attenuated PTH response in mineralized nodule formation in vitro. Remarkably, this observation was correlated with decreased PTH response of matrix Gla protein expression. Expressions of pERK1/2 and cyclin D1 were inhibited dramatically by PTH in differentiated osteoblasts from WT mice but much less in osteoblasts from Mkp1 KO mice. In conclusion, MKP1 is important for bone homeostasis, osteoblast differentiation and skeletal responsiveness to PTH.

Contribution of the sclerostin domain-containing protein 1 (SOSTDC1) gene to normal variation of peak bone mineral density in Chinese women and men

A genome-wide linkage analysis in Chinese families revealed a significant quantitative trait loci on chromosome 7p21.1 for femoral neck bone mineral density (BMD) (LOD = 3.68), and a potential candidate gene, sclerostin domain-containing protein 1 (SOSTDC1), is located in this region. SOSTDC1 belongs to a class of bone morphogenetic protein (BMP) antagonists that bind BMPs and regulate their signaling. We therefore genotyped 6 tag single nucleotide polymorphisms (tag-SNPs) in SOSTDC1 gene using allele-specific PCR method and investigated the association between SOSTDC1 gene polymorphisms and peak BMD variation in 401 Chinese female-offspring nuclear families (including 1260 subjects) and 400 Chinese male-offspring nuclear families (including 1215 subjects), respectively. Using both family-based (quantitative transmission disequilibrium test) and population-based (ANOVA) methods of analyses, BMD values were adjusted for age, height and weight. In female-offspring nuclear families, we found a significant within family association between rs16878759 and the lumbar spine peak BMD (P = 0.003) and rs16878759 accounted for 1.4% of the lumbar spine peak BMD variation. Moreover, haplotype CCC (containing rs12699800, rs16878759, and rs17619769) had a significant within family association with the lumbar spine peak BMD (P = 0.001) and accounted for 1.9% of the peak BMD variation at this bone site. However, in the male-offspring nuclear families, we failed to detect any significant association between any SNP or haplotype and peak BMD at any bone site. In conclusion, our results indicate for the first time that the genetic polymorphisms in SOSTDC1 have an effect on attainment and maintenance of peak bone mass in Chinese women.

An essential role for the association of CD47 to SHPS-1 in skeletal remodeling

Integrin-associated protein (IAP/CD47) has been implicated in macrophage-macrophage fusion. To understand the actions of CD47 on skeletal remodeling, we compared Cd47(-/-) mice with Cd47(+/+) controls. Cd47(-/-) mice weighed less and had decreased areal bone mineral density compared with controls. Cd47(-/-) femurs were shorter in length with thinner cortices and exhibited lower trabecular bone volume owing to decreased trabecular number and thickness. Histomorphometry revealed reduced bone-formation and mineral apposition rates, accompanied by decreased osteoblast numbers. No differences in osteoclast number were observed despite a nonsignificant but 40% decrease in eroded surface/bone surface in Cd47(-/-) mice. In vitro, the number of functional osteoclasts formed by differentiating Cd47(-/-) bone marrow cells was significantly decreased compared with wild-type cultures and was associated with a decrease in bone-resorption capacity. Furthermore, by disrupting the CD47-SHPS-1 association, we found that osteoclastogenesis was markedly impaired. Assays for markers of osteoclast maturation suggested that the defect was at the point of fusion and not differentiation and was associated with a lack of SHPS-1 phosphorylation, SHP-1 phosphatase recruitment, and subsequent dephosphorylation of non-muscle cell myosin IIA. We also demonstrated a significant decrease in osteoblastogenesis in bone marrow stromal cells derived from Cd47(-/-) mice. Our finding of cell-autonomous defects in Cd47(-/-) osteoblast and osteoclast differentiation coupled with the pronounced skeletal phenotype of Cd47(-/-) mice support the conclusion that CD47 plays an important role in regulating skeletal acquisition and maintenance through its actions on both bone formation and bone resorption.

Identification of candidate gene regions in the rat by co-localization of QTLs for bone density, size, structure and strength

Susceptibility to osteoporotic fracture is influenced by genetic factors that can be dissected by whole-genome linkage analysis in experimental animal crosses. The aim of this study was to characterize quantitative trait loci (QTLs) for biomechanical and two-dimensional dual-energy X-ray absorptiometry (DXA) phenotypes in reciprocal F2 crosses between diabetic GK and normo-glycemic F344 rat strains and to identify possible co-localization with previously reported QTLs for bone size and structure. The biomechanical measurements of rat tibia included ultimate force, stiffness and work to failure while DXA was used to characterize tibial area, bone mineral content (BMC) and areal bone mineral density (aBMD). F2 progeny (108 males, 98 females) were genotyped with 192 genome-wide markers followed by sex- and reciprocal cross-separated whole-genome QTL analyses. Significant QTLs were identified on chromosome 8 (tibial area; logarithm of odds (LOD) = 4.7 and BMC; LOD = 4.1) in males and on chromosome 1 (stiffness; LOD = 5.5) in females. No QTLs showed significant sex-specific interactions. In contrast, significant cross-specific interactions were identified on chromosome 2 (aBMD; LOD = 4.7) and chromosome 6 (BMC; LOD = 4.8) for males carrying F344mtDNA, and on chromosome 15 (ultimate force; LOD = 3.9) for males carrying GKmtDNA, confirming the effect of reciprocal cross on osteoporosis-related phenotypes. By combining identified QTLs for biomechanical-, size- and qualitative phenotypes (pQCT and 3D CT) from the same population, overlapping regions were detected on chromosomes 1, 3, 4, 6, 8 and 10. These are strong candidate regions in the search for genetic risk factors for osteoporosis.

Suggestive linkage to chromosome 1q for bone mineral apparent density in Brazilian sister adolescents

OBJECTIVE

To investigate linkage to chromosome 1q and 11q region for lumbar spine, femoral neck and total body BMD and volumetric BMD in Brazilian sister adolescents aged 10-20-year-old and 57 mothers.

METHODS

We evaluated 161 sister pairs (n=329) aged 10-20 years old and 57 of their mothers in this study. Physical traits and lifestyle factors were collected as covariates for lumbar spine (LS), femoral neck (FN) and total body (TB) BMD and bone mineral apparent density (BMAD). We selected nine microsatellite markers in chromosome 1q region (spanning nearly 33cM) and eight in chromosome 11q region (spanning nearly 34cM) to perform linkage analysis.

RESULTS

The highest LOD score values obtained from our data were in sister pairs LS BMAD analysis. Their values were: 1.32 (P<0.006), 2.61 (P<0.0002) and 2.44 (P<0.0004) in D1S218, D1S2640 and D1S2623 markers, respectively. No significant LOD score was found with LS and FN BMD/BMAD in chromosome 11q region. Only TB BMD showed significant linkage higher than 1.0 for chromosome 11q region in the markers D11S4191 and D11S937.

DISCUSSION/CONCLUSIONS

Our results provided suggestive linkage for LS BMAD at D1S2640 marker in adolescent sister pairs and suggest a possible candidate gene (LHX4) related to adolescent LS BMAD in this region. These results reinforce chromosome 1q21-23 as a candidate region to harbor one or more bone formation/maintenance gene. In the other hand, it did not repeat for chromosome 11q12-13 in our population.

BMD regulation on mouse distal chromosome 1, candidate genes, and response to ovariectomy or dietary fat

The distal end of mouse chromosome 1 (Chr 1) harbors quantitative trait loci (QTLs) that regulate bone mineral density (BMD) and share conserved synteny with human chromosome 1q. The objective of this article was to map this mouse distal Chr 1 region and identify gene(s) responsible for BMD regulation in females. We used X-ray densitometry [ie, dual-energy X-ray Absorptiometry (DXA), micro-computed tomography (µCT), and peripheral quantitative computed tomography (pQCT)] to phenotype a set of nested congenic strains constructed from C57BL/6BmJ (B6/Bm) and C3H/HeJ (C3H) mice to map the region associated with the BMD QTL. The critical region has been reduced to an interval of 0.152 Mb that contributes to increased BMD when C3H alleles are present. Histomorphometry and osteoblast cultures indicated that increased osteoblast activity was associated with increased BMD in mouse strains with C3H alleles in this critical region. This region contains two genes, Aim2, which binds with cytoplasmic dsDNA and results in apoptosis, and AC084073.22, a predicted gene of unknown function. Ovariectomy induced bone loss in the B6/Bm progenitor and the three congenic strains regardless of the alleles present in the critical BMD region. High dietary fat treatment (thought to suppress distal Chr 1 QTL for BMD in mice) did not induce bone loss in the congenics carrying C3H alleles in the critical 0.152 Mb carrying the AIM2 and AC084073.22 genes. Gene expression studies in whole bone of key congenics showed differential expression of AC084073.22 for strains carrying B6/Bm versus C3H alleles but not for Aim2. In conclusion, our data suggest that osteoblasts are the cellular target of gene action and that AC084073.22 is the best candidate for female BMD regulation in the distal region of mouse Chr 1.

Genetic loci for bone architecture determined by three-dimensional CT in crosses with the diabetic GK rat

The F344 rat carries alleles contributing to bone fragility while the GK rat spontaneously develops type-2 diabetes. These characteristics make F344×GK crosses well suited for the identification of genes related to bone size and allow for future investigation on the association with type-2 diabetes. The aim of this study was to identify quantitative trait loci (QTLs) for bone size phenotypes measured by a new application of three-dimensional computed tomography (3DCT) and to investigate the effects of sex- and reciprocal cross. Tibia from male and female GK and F344 rats, representing the parental, F1 and F2 generations, were examined with 3DCT and analyzed for: total and cortical volumetric BMD, straight and curved length, peri- and endosteal area at mid-shaft. F2 progeny (108 male and 98 female) were genotyped with 192 genome-wide microsatellite markers (average distance 10 cM). Sex- and reciprocal cross-separated QTL analyses were performed for the identification of QTLs linked to 3DCT phenotypes and true interactions were confirmed by likelihood ratio analysis in all F2 animals. Several genome-wide significant QTLs were found in the sex- and reciprocal cross-separated progeny on chromosomes (chr) 1, 3, 4, 9, 10, 14, and 17. Overlapping QTLs for both males and females in the (GK×F344)F2 progeny were located on chr 1 (39-67 cM). This region confirms previously reported pQCT QTLs and overlaps loci for fasting glucose. Sex separated linkage analysis confirmed a male specific QTL on chr 9 (67-82 cM) for endosteal area at the fibula site. Analyses separating the F2 population both by sex and reciprocal cross identified cross specific QTLs on chr 14 (males) and chr 3 and 4 (females). Two loci, chr 4 and 6, are unique to 3DCT and separate from pQCT generated loci. The 3DCT method was highly reproducible and provided high precision measurements of bone size in the rat enabling identification of new sex- and cross-specific loci. The QTLs on chr 1 indicate potential genetic association between bone-related phenotypes and traits affecting type-2 diabetes. The results illustrate the complexity of the genetic architecture of bone size phenotypes and demonstrate the importance of complementary methods for bone analysis.

Identification of genetic determinants of IGF-1 levels and longevity among mouse inbred strains

The IGF-1 signaling pathway plays an important role in regulating longevity. To identify the genetic loci and genes that regulate plasma IGF-1 levels, we intercrossed MRL/MpJ and SM/J, inbred mouse strains that differ in IGF-1 levels. Quantitative trait loci (QTL) analysis of IGF-1 levels of these F2 mice detected four QTL on chromosomes (Chrs) 9 (48 Mb), 10 (86 Mb), 15 (18 Mb), and 17 (85 Mb). Haplotype association mapping of IGF-1 levels in 28 domesticated inbred strains identified three suggestive loci in females on Chrs 2 (13 Mb), 10 (88 Mb), and 17 (28 Mb) and in four males on Chrs 1 (159 Mb), 3 (52 and 58 Mb), and 16 (74 Mb). Except for the QTL on Chr 9 and 16, all loci co-localized with IGF-1 QTL previously identified in other mouse crosses. The most significant locus was the QTL on Chr 10, which contains the Igf1 gene and which had a LOD score of 31.8. Haplotype analysis among 28 domesticated inbred strains revealed a major QTL on Chr 10 overlapping with the QTL identified in the F2 mice. This locus showed three major haplotypes; strains with haplotype 1 had significantly lower plasma IGF-1 and extended longevity (P < 0.05) than strains with haplotype 2 or 3. Bioinformatic analysis, combined with sequencing and expression studies, showed that Igf1 is the most likely QTL gene, but that other genes may also play a role in this strong QTL.

Mouse BMD quantitative trait loci show improved concordance with human genome-wide association loci when recalculated on a new, common mouse genetic map

Bone mineral density (BMD) is a heritable trait, and in mice, over 100 quantitative trait loci (QTLs) have been reported, but candidate genes have been identified for only a small percentage. Persistent errors in the mouse genetic map have negatively affected QTL localization, spurring the development of a new, corrected map. In this study, QTLs for BMD were remapped in 11 archival mouse data sets using this new genetic map. Since these QTLs all were mapped in a comparable way, direct comparisons of QTLs for concordance would be valid. We then compared human genome-wide association study (GWAS) BMD loci with the mouse QTLs. We found that 26 of the 28 human GWAS loci examined were located within the confidence interval of a mouse QTL. Furthermore, 14 of the GWAS loci mapped to within 3 cM of a mouse QTL peak. Lastly, we demonstrated that these newly remapped mouse QTLs can substantiate a candidate gene for a human GWAS locus, for which the peak single-nucleotide polymorphism (SNP) fell in an intergenic region. Specifically, we suggest that MEF2C (human chromosome 5, mouse chromosome 13) should be considered a candidate gene for the genetic regulation of BMD. In conclusion, use of the new mouse genetic map has improved the localization of mouse BMD QTLs, and these remapped QTLs show high concordance with human GWAS loci. We believe that this is an opportune time for a renewed effort by the genetics community to identify the causal variants regulating BMD using a synergistic mouse-human approach.

A novel spontaneous mutation of Irs1 in mice results in hyperinsulinemia, reduced growth, low bone mass and impaired adipogenesis

A spontaneous mouse mutant, designated 'small' (sml), was recognized by reduced body size suggesting a defect in the IGF1/GH axis. The mutation was mapped to the chromosome 1 region containing Irs1, a viable candidate gene whose sequence revealed a single nucleotide deletion resulting in a premature stop codon. Despite normal mRNA levels in mutant and control littermate livers, western blot analysis revealed no detectable protein in mutant liver lysates. When compared with the control littermates, Irs1(sml)/Irs1(sml) (Irs1(sml/sml)) mice were small, lean, hearing impaired; had 20% less serum IGF1; were hyperinsulinemic; and were mildly insulin resistant. Irs1(sml/sml) mice had low bone mineral density, reduced trabecular and cortical thicknesses, and low bone formation rates, while osteoblast and osteoclast numbers were increased in the females but not different in the males compared with the Irs1(+/+) controls. In vitro, Irs1(sml/sml) bone marrow stromal cell cultures showed decreased alkaline phosphatase-positive colony forming units (pre-osteoblasts; CFU-AP+) and normal numbers of tartrate-resistant acid phosphatase-positive osteoclasts. Irs1(sml/sml) stromal cells treated with IGF1 exhibited a 50% decrease in AKT phosphorylation, indicative of defective downstream signaling. Similarities between engineered knockouts and the spontaneous mutation of Irs1(sml) were identified as well as significant differences with respect to heterozygosity and gender. In sum, we have identified a spontaneous mutation in the Irs1 gene associated with a major skeletal phenotype. Changes in the heterozygous Irs1(+)(/sml) mice raise the possibility that similar mutations in humans are associated with short stature or osteoporosis.

Variation in within-bone stiffness measured by nanoindentation in mice bred for high levels of voluntary wheel running

The hierarchical structure of bone, involving micro-scale organization and interaction of material components, is a critical determinant of macro-scale mechanics. Changes in whole-bone morphology in response to the actions of individual genes, physiological loading during life, or evolutionary processes, may be accompanied by alterations in underlying mineralization or architecture. Here, we used nanoindentation to precisely measure compressive stiffness in the femoral mid-diaphysis of mice that had experienced 37 generations of selective breeding for high levels of voluntary wheel running (HR). Mice (n = 48 total), half from HR lines and half from non-selected control (C) lines, were divided into two experimental groups, one with 13-14 weeks of access to a running wheel and one housed without wheels (n = 12 in each group). At the end of the experiment, gross and micro-computed tomography (microCT)-based morphometric traits were measured, and reduced elastic modulus (E(r)) was estimated separately for four anatomical quadrants of the femoral cortex: anterior, posterior, lateral, and medial. Two-way, mixed-model analysis of covariance (ancova) showed that body mass was a highly significant predictor of all morphometric traits and that structural change is more apparent at the microCT level than in conventional morphometrics of whole bones. Both line type (HR vs. C) and presence of the mini-muscle phenotype (caused by a Mendelian recessive allele and characterized by a approximately 50% reduction in mass of the gastrocnemius muscle complex) were significant predictors of femoral cortical cross-sectional anatomy. Measurement of reduced modulus obtained by nanoindentation was repeatable within a single quadrant and sensitive enough to detect inter-individual differences. Although we found no significant effects of line type (HR vs. C) or physical activity (wheel vs. no wheel) on mean stiffness, anterior and posterior quadrants were significantly stiffer (P < 0.0001) than medial and lateral quadrants (32.67 and 33.09 GPa vs. 29.78 and 30.46 GPa, respectively). Our findings of no significant difference in compressive stiffness in the anterior and posterior quadrants agree with previous results for mice, but differ from those for large mammals. Integrating these results with others from ongoing research on these mice, we hypothesize that the skeletons of female HR mice may be less sensitive to the effects of chronic exercise, due to decreased circulating leptin levels and potentially altered endocannabinoid signaling.

Nephroblastoma overexpressed (Nov) inactivation sensitizes osteoblasts to bone morphogenetic protein-2, but nov is dispensable for skeletal homeostasis

Overexpression of nephroblastoma overexpressed (Nov), a member of the Cyr 61, connective tissue growth factor, Nov family of proteins, inhibits osteoblastogenesis and causes osteopenia. The consequences of Nov inactivation on osteoblastogenesis and the postnatal skeleton are not known. To study the function of Nov, we inactivated Nov by homologous recombination. Nov null mice were maintained in a C57BL/6 genetic background after the removal of the neomycin selection cassette and compared with wild-type controls of identical genetic composition. Nov null mice were identified by genotyping and absent Nov mRNA in calvarial extracts and osteoblast cultures. Nov null mice did not exhibit developmental skeletal abnormalities or postnatal changes in weight, femoral length, body fat, or bone mineral density and appeared normal. Bone volume and trabecular number were decreased only in 1-month-old female mice. In older mice, after 7 months of age, osteoblast surface and bone formation were increased in females, and osteoclast and eroded surfaces were increased in male Nov null mice. Calvarial osteoblasts from Nov null mice displayed enhanced alkaline phosphatase activity, alkaline phosphatase mRNA, and transactivation of a bone morphogenetic protein (BMP)/phosphorylated mothers against decapentaplegic reporter construct in response to BMP-2. Similar results were obtained after the down-regulation of Nov by RNA interference in ST-2 stromal and MC3T3 cells. Osteoclast number was increased in marrow stromal cell cultures from Nov null mice. Surface plasmon resonance demonstrated direct interactions between Nov and BMP-2. In conclusion, Nov sensitizes osteoblasts to BMP-2, but Nov is dispensable for the maintenance of bone mass.

Systems analysis of bone

The genetic variants contributing to variability in skeletal traits has been well studied, and several hundred QTLs have been mapped and several genes contributing to trait variation have been identified. However, many questions remain unanswered. In particular, it is unclear whether variation in a single gene leads to alterations in function. Bone is a highly adaptive system and genetic variants affecting one trait are often accompanied by compensatory changes in other traits. The functional interactions among traits, which is known as phenotypic integration, has been observed in many biological systems, including bone. Phenotypic integration is a property of bone that is critically important for establishing a mechanically functional structure that is capable of supporting the forces imparted during daily activities. In this paper, bone is reviewed as a system and primarily in the context of functionality. A better understanding of the system properties of bone will lead to novel targets for future genetic analyses and the identification of genes that are directly responsible for regulating bone strength. This systems analysis has the added benefit of leaving a trail of valuable information about how the skeletal system works. This information will provide novel approaches to assessing skeletal health during growth and aging and for developing novel treatment strategies to reduce the morbidity and mortality associated with fragility fractures.

Quantitative trait loci for biomechanical performance and femoral geometry in an intercross of recombinant congenic mice: restriction of the Bmd7 candidate interval

Despite steady progress in identifying quantitative trait loci (QTLs) for bone phenotypes, relatively little progress has been made in moving from QTLs to identifying the relevant gene. We exploited the genetic structure of recombinant congenic mouse strains by performing a reciprocal intercross of the strains HcB-8 and HcB-23, phenotyped for body size, femoral biomechanical performance, and femoral diaphyseal geometry and mapped with R/qtl and QTL Cartographer. Significant QTLs are present on chromosomes 1, 2, 3, 4, 6, and 10. We found significant sex x QTL and cross-direction x QTL interactions. The chromosome 4 QTL affects multiple femoral anatomic features and biomechanical properties. The known segregating segment of chromosome 4 contains only 18 genes, among which Ece1, encoding endothelin-converting enzyme 1, stands out as a candidate. Endothelin signaling has been shown to promote the growth of osteoblastic metastases and to potentiate signaling via the Wnt pathway. The colocalizing chromosome 4 QTL Bmd7 (for bone mineral density 7) increases responsiveness to mechanical loading. By exploiting the short informative segment of chromosome 4 and the known biology, we propose that Ece1 is the gene responsible for Bmd7 and that it acts by increasing responsiveness to mechanical loading through modulation of Wnt signaling.

Quantitative trait loci, genes, and polymorphisms that regulate bone mineral density in mouse

This is an in silico analysis of data available from genome-wide scans. Through analysis of QTL, genes and polymorphisms that regulate BMD, we identified 82 BMD QTL, 191 BMD-associated (BMDA) genes, and 83 genes containing known BMD-associated polymorphisms (BMDAP). The catalogue of all BMDA/BMDAP genes and relevant literatures are provided. In total, there are substantially more BMDA/BMDAP genes in regions of the genome where QTL have been identified than in non-QTL regions. Among 191 BMDA genes and 83 BMDAP genes, 133 and 58 are localized in QTL regions, respectively. The difference was still noticeable for the chromosome distribution of these genes between QTL and non-QTL regions. These results have allowed us to generate an integrative profile of QTL, genes, polymorphisms that determine BMD. These data could facilitate more rapid and comprehensive identification of causal genes underlying the determination of BMD in mouse and provide new insights into how BMD is regulated in humans.

Gender- and compartment-specific bone loss in C57BL/6J mice: correlation to season?

Seasonal variation in bone mineral density (BMD) has been documented in humans, and has been attributed to changes in 25-hydroxyvitamin D [25(OH)D] synthesis. To test the hypothesis that seasonal changes in bone mass occur in laboratory mice, we measured body composition, femoral bone phenotypes, and serum bone markers in 16-wk-old male and female C57BL/6 (B6) mice during the summer (June-August) and winter (December-February) months at The Jackson Laboratory in Bar Harbor, Maine. Both male and female B6 mice had higher volumetric BMD in the summer than winter. Females showed reduced trabecular bone, whereas males showed changes in bone volume. Males, but not females, had higher insulin-like growth factor 1 in summer than in winter, and only males showed an increase in body weight during the winter. No seasonal differences in serum TRAP5b, osteocalcin, or 25(OH)D were noted for either sex. We conclude that seasonal variation in skeletal and body composition parameters in B6 mice is significant and must be considered when performing longitudinal phenotyping of the skeleton. Further studies are needed to determine the environmental factors that cue seasonal changes in body composition and the mechanisms that produce these changes.

Quantitative trait loci for BMD in an SM/J by NZB/BlNJ intercross population and identification of Trps1 as a probable candidate gene

Identification of genes that regulate BMD will enhance our understanding of osteoporosis and could provide novel molecular targets for treatment or prevention. We generated a mouse intercross population and carried out a quantitative trait locus (QTL) analysis of 143 female and 124 male F(2) progeny from progenitor strains SM/J and NZB/BlNJ using whole body and vertebral areal BMD (aBMD) as measured by DXA. We found that both whole body and vertebral aBMD was affected by two loci on chromosome 9: one with a significant epistatic interaction on distal chromosome 8 and the other with a sex-specific effect. Two additional significant QTLs were identified on chromosome 12, and several suggestive ones were identified on chromosomes 5, 8, 15, and 19. The chromosome 9, 12, and 15 loci have been previously identified in other crosses. SNP-based haplotype analysis of the progenitor strains identified blocks within the QTL region that distinguish the low allele strains from the high allele strains, significantly narrowing the QTL region and reducing the possible candidate genes to 98 for chromosome 9, 31 for chromosome 12, and only 2 for chromosome 15. Trps1 is the most probable candidate gene for the chromosome 15 QTL. The sex-specific effects may help to elucidate the BMD differences between males and females. This study shows the power of statistical modeling to resolve linked QTLs and the use of haplotype analysis in narrowing the list of candidates.

Genomic expression analysis of rat chromosome 4 for skeletal traits at femoral neck

Hip fracture is the most devastating osteoporotic fracture type with significant morbidity and mortality. Several studies in humans and animal models identified chromosomal regions linked to hip size and bone mass. Previously, we identified that the region of 4q21-q41 on rat chromosome (Chr) 4 harbors multiple femoral neck quantitative trait loci (QTLs) in inbred Fischer 344 (F344) and Lewis (LEW) rats. The purpose of this study is to identify the candidate genes for femoral neck structure and density by correlating gene expression in the proximal femur with the femoral neck phenotypes linked to the QTLs on Chr 4. RNA was extracted from proximal femora of 4-wk-old rats from F344 and LEW strains, and two other strains, Copenhagen 2331 and Dark Agouti, were used as a negative control. Microarray analysis was performed using Affymetrix Rat Genome 230 2.0 arrays. A total of 99 genes in the 4q21-q41 region were differentially expressed (P < 0.05) among all strains of rats with a false discovery rate <10%. These 99 genes were then ranked based on the strength of correlation between femoral neck phenotypes measured in F2 animals, homozygous for a particular strain's allele at the Chr 4 QTL and the expression level of the gene in that strain. A total of 18 candidate genes were strongly correlated (r(2) > 0.50) with femoral neck width and prioritized for further analysis. Quantitative PCR analysis confirmed 14 of 18 of the candidate genes. Ingenuity pathway analysis revealed several direct or indirect relationships among the candidate genes related to angiogenesis (VEGF), bone growth (FGF2), bone formation (IGF2 and IGF2BP3), and resorption (TNF). This study provides a shortened list of genetic determinants of skeletal traits at the hip and may lead to novel approaches for prevention and treatment of hip fracture.

Complex genetic regulation of bone mineral density and insulin-like growth factor-I in C57BL/6J-Chr #A/J/NaJ chromosome substitution strains

Low bone mineral density (BMD) is a phenotype associated with osteoporosis and increased risk of fracture. Since 60-80% of variation in BMD is associated with genetic factors, we used the novel approach of chromosome substitution strains (CSS) to identify chromosomes that harbor genes that regulate BMD. Specifically, we evaluated 24 wk old C57BL/6J-Chr #(A/J)/NaJ CSS (n = 7 to 18) in which each chromosome in the host C57BL/6J strain is replaced by a corresponding chromosome from the donor A/J strain (19 autosomes, X, Y). We determined several A/J chromosomes contribute to body weight (BW), percent body fat (BF), whole body areal BMD (aBMD), and serum insulin-like growth factor (IGF)-I in both a positive and negative manner when compared with C57BL/6J. Specifically, C57BL/6J-Chr 5(A/J)/NaJ (B.A-5) (males) and B.A-13 (females) contributed to increased BW, whereas B.A-3, 4, 8, 9, 12, and 18 (males) and B.A-3, 4, and 11 (females) contributed to reduced BW. B.A-5 (males) and B.A-13 (females) contributed to increased BF, whereas B.A-12 (males) and B.A-3, 4, 10, and 11 (females) contributed to reduced BF. B.A-14 (females) contributed to increased aBMD and B.A-1, 2, 6, 9, 10, and 18 (males) and B.A-8, 9, and 10 (females) contributed to reduced aBMD. To determine if similar chromosomes regulate aBMD and IGF-I, we determined serum concentrations of IGF-I. B.A-14 and Y (males) and B.A-6 (females) contributed to increased IGF-I and male B.A-3 and female B.A-8 contributed to reduced IGF-I. Overall, we identified several sex-dependent and -independent chromosomes that regulate aBMD and IGF-I in adult mice.

Genetic variation in the patterns of skeletal progenitor cell differentiation and progression during endochondral bone formation affects the rate of fracture healing

These studies examined how genetic differences that regulate architectural and bone material properties would be expressed during fracture healing and determine whether any of these features would affect rates of healing as defined by regain of strength. Controlled fractures were generated in three inbred strains of mice: A/J, C57Bl/6J (B6), and C3H/HeJ (C3H). Both the A/J and B6 strains showed faster healing than the C3H strain based on regains in strength and stiffness. Strain-specific architectural features such as moment of inertia, cross-sectional area, and cortical thickness were all recapitulated during the development of the callus tissues. None of these traits were directly relatable to rates of fracture healing. However, rates of healing were related to variations in the temporal patterns of chondrogenic and osteogenic lineage development. The B6 strain expressed the highest percentage of cartilage gene products and had the longest period of chondrocyte maturation and hypertrophy. The slowest healing strain (C3H) had the shortest period of chondrogenic development and earliest initiation of osteogenic development. Although the A/J strain showed an almost identical pattern of chondrogenic development as the C3H strain, A/J initiated osteogenic development several days later than C3H during fracture healing. Long bone growth plates at 28 days after birth showed similar strain-specific variation in cartilage tissue development as seen in fracture healing. Thus, the B6 strain had the largest growth plate heights, cell numbers per column, and the largest cell size, whereas the C3H columns were the shortest, had the smallest number of cells per column, and showed the smallest cell sizes. These results show that (1) different strains of mice express variations of skeletal stem cell lineage differentiation and (2) that these variations affect the rate of fracture healing.

Contribution of gender-specific genetic factors to osteoporosis risk

Common diseases result from the complex relationship between genetic and environmental factors. The aim of this review is to provide perspective for a conceptual framework aimed at studying the interplay of gender-specific genetic and environmental factors in the etiology of complex disease, using osteoporosis as an example. In recent years, gender differences in the heritability of the osteoporosis-related phenotypes have been reported and sex-specific quantitative-trait loci were discovered by linkage studies in humans and mice. Results of numerous allelic association studies also differed by gender. In most cases, it was not clear whether or not this phenomenon should be attributed to the effect of sex-chromosomes, sex hormones, or other intrinsic or extrinsic differences between the genders, such as the level of bioavailable estrogen and of physical activity. We conclude that there is need to consider gender-specific genetic and environmental factors in the planning of future association studies on the etiology of osteoporosis and other complex diseases prevalent in the general population.

Gender-specific changes in bone turnover and skeletal architecture in igfbp-2-null mice

IGF-binding protein-2 (IGFBP-2) is a 36-kDa protein that binds to the IGFs with high affinity. To determine its role in bone turnover, we compared Igfbp2(-/-) mice with Igfbp2(+/+) colony controls. Igfbp2(-/-) males had shorter femurs and were heavier than controls but were not insulin resistant. Serum IGF-I levels in Igfbp2(-/-) mice were 10% higher than Igfbp2(+/+) controls at 8 wk of age; in males, this was accompanied by a 3-fold increase in hepatic Igfbp3 and Igfbp5 mRNA transcripts compared with Igfbp2(+/+) controls. The skeletal phenotype of the Igfbp2(-/-) mice was gender and compartment specific; Igfbp2(-/-) females had increased cortical thickness with a greater periosteal circumference compared with controls, whereas male Igfbp2(-/-) males had reduced cortical bone area and a 20% reduction in the trabecular bone volume fraction due to thinner trabeculae than Igfbp2(+/+) controls. Serum osteocalcin levels were reduced by nearly 40% in Igfbp2(-/-) males, and in vitro, both CFU-ALP(+) preosteoblasts, and tartrate-resistant acid phosphatase-positive osteoclasts were significantly less abundant than in Igfbp2(+/+) male mice. Histomorphometry confirmed fewer osteoblasts and osteoclasts per bone perimeter and reduced bone formation in the Igfbp2(-/-) males. Lysates from both osteoblasts and osteoclasts in the Igfbp2(-/-) males had phosphatase and tensin homolog (PTEN) levels that were significantly higher than Igfbp2(+/+) controls and were suppressed by addition of exogenous IGFBP-2. In summary, there are gender- and compartment-specific changes in Igfbp2(-/-) mice. IGFBP-2 may regulate bone turnover in both an IGF-I-dependent and -independent manner.