Search

Skull bone mineral density (SK-BMD) provides a suitable trait for the discovery of key genes in bone biology, particularly to intramembranous ossification, not captured at other skeletal sites. We perform a genome-wide association meta-analysis (n ~ 43,800) of SK-BMD, identifying 59 loci, collectively explaining 12.5% of the trait variance. Association signals cluster within gene-sets involved in skeletal development and osteoporosis. Among the four novel loci (ZIC1, PRKAR1A, AZIN1/ATP6V1C1, GLRX3), there are factors implicated in intramembranous ossification and as we show, inherent to craniosynostosis processes. Functional follow-up in zebrafish confirms the importance of ZIC1 on cranial suture patterning. Likewise, we observe abnormal cranial bone initiation that culminates in ectopic sutures and reduced BMD in mosaic atp6v1c1 knockouts. Mosaic prkar1a knockouts present asymmetric bone growth and, conversely, elevated BMD. In light of this evidence linking SK-BMD loci to craniofacial abnormalities, our study provides new insight into the pathophysiology, diagnosis and treatment of skeletal diseases., (© 2023. The Author(s).)

The "omnigenic" model of the genetic architecture of complex traits proposed two categories of causal genes: core and peripheral. Core genes are hypothesized to directly regulate disease and may serve as therapeutic targets. Using a cell-type- and time-point-specific gene co-expression network for mineralizing osteoblasts, we identify a co-expression module enriched for genes implicated by bone mineral density (BMD) genome-wide association studies (GWASs), correlated with in vitro osteoblast mineralization and associated with skeletal phenotypes in human monogenic disease and mouse knockouts. Four genes from this module (B4GALNT3, CADM1, DOCK9, and GPR133) are located within the BMD GWAS loci with colocalizing expression quantitative trait loci (eQTL) and exhibit altered BMD in mouse knockouts, suggesting that they are causal genetic drivers of BMD in humans. Our network-based approach identifies a "core" module for BMD and provides a resource for expanding our understanding of the genetics of bone mass., Competing Interests: Declaration of Interests The authors declare no competing interests., (Copyright © 2020 The Authors. Published by Elsevier Inc. All rights reserved.)

Osteoporosis is a genetic disease characterized by progressive reductions in bone mineral density (BMD) leading to an increased risk of fracture. Over the last decade, genome-wide association studies (GWASs) have identified over 1000 associations for BMD. However, as a phenotype BMD is challenging as bone is a multicellular tissue affected by both local and systemic physiology. Here, we focused on a single component of BMD, osteoblast-mediated bone formation in mice, and identified associations influencing osteoblast activity on mouse Chromosomes (Chrs) 1, 4, and 17. The locus on Chr. 4 was in an intergenic region between Wnt4 and Zbtb40, homologous to a locus for BMD in humans. We tested both Wnt4 and Zbtb40 for a role in osteoblast activity and BMD. Knockdown of Zbtb40, but not Wnt4, in osteoblasts drastically reduced mineralization. Additionally, loss-of-function mouse models for both genes exhibited reduced BMD. Our results highlight that investigating the genetic basis of in vitro osteoblast mineralization can be used to identify genes impacting bone formation and BMD., Competing Interests: The authors have declared that no competing interests exist.

Osteoporosis is a common aging-related disease diagnosed primarily using bone mineral density (BMD). We assessed genetic determinants of BMD as estimated by heel quantitative ultrasound in 426,824 individuals, identifying 518 genome-wide significant loci (301 novel), explaining 20% of its variance. We identified 13 bone fracture loci, all associated with estimated BMD (eBMD), in ~1.2 million individuals. We then identified target genes enriched for genes known to influence bone density and strength (maximum odds ratio (OR) = 58, P = 1 × 10 -75 ) from cell-specific features, including chromatin conformation and accessible chromatin sites. We next performed rapid-throughput skeletal phenotyping of 126 knockout mice with disruptions in predicted target genes and found an increased abnormal skeletal phenotype frequency compared to 526 unselected lines (P < 0.0001). In-depth analysis of one gene, DAAM2, showed a disproportionate decrease in bone strength relative to mineralization. This genetic atlas provides evidence linking associated SNPs to causal genes, offers new insight into osteoporosis pathophysiology, and highlights opportunities for drug development.

Bone mineral density (BMD) assessed by DXA is used to evaluate bone health. In children, total body (TB) measurements are commonly used; in older individuals, BMD at the lumbar spine (LS) and femoral neck (FN) is used to diagnose osteoporosis. To date, genetic variants in more than 60 loci have been identified as associated with BMD. To investigate the genetic determinants of TB-BMD variation along the life course and test for age-specific effects, we performed a meta-analysis of 30 genome-wide association studies (GWASs) of TB-BMD including 66,628 individuals overall and divided across five age strata, each spanning 15 years. We identified variants associated with TB-BMD at 80 loci, of which 36 have not been previously identified; overall, they explain approximately 10% of the TB-BMD variance when combining all age groups and influence the risk of fracture. Pathway and enrichment analysis of the association signals showed clustering within gene sets implicated in the regulation of cell growth and SMAD proteins, overexpressed in the musculoskeletal system, and enriched in enhancer and promoter regions. These findings reveal TB-BMD as a relevant trait for genetic studies of osteoporosis, enabling the identification of variants and pathways influencing different bone compartments. Only variants in ESR1 and close proximity to RANKL showed a clear effect dependency on age. This most likely indicates that the majority of genetic variants identified influence BMD early in life and that their effect can be captured throughout the life course., (Copyright © 2017 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Osteoporosis is a common disease diagnosed primarily by measurement of bone mineral density (BMD). We undertook a genome-wide association study (GWAS) in 142,487 individuals from the UK Biobank to identify loci associated with BMD as estimated by quantitative ultrasound of the heel. We identified 307 conditionally independent single-nucleotide polymorphisms (SNPs) that attained genome-wide significance at 203 loci, explaining approximately 12% of the phenotypic variance. These included 153 previously unreported loci, and several rare variants with large effect sizes. To investigate the underlying mechanisms, we undertook (1) bioinformatic, functional genomic annotation and human osteoblast expression studies; (2) gene-function prediction; (3) skeletal phenotyping of 120 knockout mice with deletions of genes adjacent to lead independent SNPs; and (4) analysis of gene expression in mouse osteoblasts, osteocytes and osteoclasts. The results implicate GPC6 as a novel determinant of BMD, and also identify abnormal skeletal phenotypes in knockout mice associated with a further 100 prioritized genes.

Purpose: Familial exudative vitreoretinopathy (FEVR) is caused by mutations in the genes encoding low-density lipoprotein receptor-related protein (LRP5) or its interacting partners, namely frizzled class receptor 4 (FZD4) and norrin cystine knot growth factor (NDP). Mouse models for Lrp5 , Fzd4 , and Ndp have proven to be important for understanding the retinal pathophysiology underlying FEVR and systemic abnormalities related to defective Wnt signaling. Here, we report a new mouse mutant, tvrm111B , which was identified by electroretinogram (ERG) screening of mice generated in the Jackson Laboratory Translational Vision Research Models (TVRM) mutagenesis program., Methods: ERGs were used to examine outer retinal physiology. The retinal vasculature was examined by in vivo retinal imaging, as well as by histology and immunohistochemistry. The tvrm111B locus was identified by genetic mapping of mice generated in a cross to DBA/2J, and subsequent sequencing analysis. Gene expression was examined by real-time PCR of retinal RNA. Bone mineral density (BMD) was examined by peripheral dual-energy X-ray absorptiometry., Results: The tvrm111B allele is inherited as an autosomal recessive trait. Genetic mapping of the decreased ERG b-wave phenotype of tvrm111B mice localized the mutation to a region on chromosome 19 that included Lrp5 . Sequencing of Lrp5 identified the insertion of a cytosine (c.4724&#95;4725insC), which is predicted to cause a frameshift that disrupts the last three of five conserved PPPSPxS motifs in the cytoplasmic domain of LRP5, culminating in a premature termination. In addition to a reduced ERG b-wave, Lrp5 tvrm111B homozygotes have low BMD and abnormal features of the retinal vasculature that have been reported previously in Lrp5 mutant mice, including persistent hyaloid vessels, leakage on fluorescein angiography, and an absence of the deep retinal capillary bed., Conclusions: The phenotype of the Lrp5 tvrm111B mutant includes abnormalities of the retinal vasculature and of BMD. This model may be a useful resource to further our understanding of the biological role of LRP5 and to evaluate experimental therapies for FEVR or other conditions associated with LRP5 dysfunction.

Genome-wide association studies (GWASs) have revealed numerous loci for areal bone mineral density (aBMD). We completed the first GWAS meta-analysis (n = 15,275) of lumbar spine volumetric BMD (vBMD) measured by quantitative computed tomography (QCT), allowing for examination of the trabecular bone compartment. SNPs that were significantly associated with vBMD were also examined in two GWAS meta-analyses to determine associations with morphometric vertebral fracture (n = 21,701) and clinical vertebral fracture (n = 5893). Expression quantitative trait locus (eQTL) analyses of iliac crest biopsies were performed in 84 postmenopausal women, and murine osteoblast expression of genes implicated by eQTL or by proximity to vBMD-associated SNPs was examined. We identified significant vBMD associations with five loci, including: 1p36.12, containing WNT4 and ZBTB40; 8q24, containing TNFRSF11B; and 13q14, containing AKAP11 and TNFSF11. Two loci (5p13 and 1p36.12) also contained associations with radiographic and clinical vertebral fracture, respectively. In 5p13, rs2468531 (minor allele frequency [MAF] = 3%) was associated with higher vBMD (β = 0.22, p = 1.9 × 10 -8 ) and decreased risk of radiographic vertebral fracture (odds ratio [OR] = 0.75; false discovery rate [FDR] p = 0.01). In 1p36.12, rs12742784 (MAF = 21%) was associated with higher vBMD (β = 0.09, p = 1.2 × 10 -10 ) and decreased risk of clinical vertebral fracture (OR = 0.82; FDR p = 7.4 × 10 -4 ). Both SNPs are noncoding and were associated with increased mRNA expression levels in human bone biopsies: rs2468531 with SLC1A3 (β = 0.28, FDR p = 0.01, involved in glutamate signaling and osteogenic response to mechanical loading) and rs12742784 with EPHB2 (β = 0.12, FDR p = 1.7 × 10 -3 , functions in bone-related ephrin signaling). Both genes are expressed in murine osteoblasts. This is the first study to link SLC1A3 and EPHB2 to clinically relevant vertebral osteoporosis phenotypes. These results may help elucidate vertebral bone biology and novel approaches to reducing vertebral fracture incidence. © 2016 American Society for Bone and Mineral Research., (© 2016 American Society for Bone and Mineral Research.)

Bisphosphonates (BPs) have been in use for many years for the treatment of osteoporosis, multiple myeloma, Paget's disease, as well as a variety of other diseases in which there is reduced bone mineral density. Given that bisphosphonates inhibit bone resorption, an important stage of fracture healing; this class of compounds has been widely studied in preclinical models regarding their influence on fracture healing. In animal models, bisphosphonate treatment is associated with a larger fracture callus, coincident with a delay in remodeling from primary woven bone to lamellar bone, but there is no delay in formation of the fracture callus. In humans, de novo use of bisphosphonate therapy after fracture does not appear to have a significant effect on fracture healing. Rarely, patients with long term use of Bisphosphonates may develop an atypical fracture and delay in fracture healing has been observed. In summary, bisphosphonates appear safe for use in the setting of acute fracture management in the upper and lower extremity in humans. While much remains unknown about the effects on healing of long-term bisphosphonates, use prior to "typical" fracture, in the special case of atypical fracture, evidence suggests that bisphosphonates negatively influence healing., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

The extent to which low-frequency (minor allele frequency (MAF) between 1-5%) and rare (MAF ≤ 1%) variants contribute to complex traits and disease in the general population is mainly unknown. Bone mineral density (BMD) is highly heritable, a major predictor of osteoporotic fractures, and has been previously associated with common genetic variants, as well as rare, population-specific, coding variants. Here we identify novel non-coding genetic variants with large effects on BMD (ntotal = 53,236) and fracture (ntotal = 508,253) in individuals of European ancestry from the general population. Associations for BMD were derived from whole-genome sequencing (n = 2,882 from UK10K (ref. 10); a population-based genome sequencing consortium), whole-exome sequencing (n = 3,549), deep imputation of genotyped samples using a combined UK10K/1000 Genomes reference panel (n = 26,534), and de novo replication genotyping (n = 20,271). We identified a low-frequency non-coding variant near a novel locus, EN1, with an effect size fourfold larger than the mean of previously reported common variants for lumbar spine BMD (rs11692564(T), MAF = 1.6%, replication effect size = +0.20 s.d., Pmeta = 2 × 10(-14)), which was also associated with a decreased risk of fracture (odds ratio = 0.85; P = 2 × 10(-11); ncases = 98,742 and ncontrols = 409,511). Using an En1(cre/flox) mouse model, we observed that conditional loss of En1 results in low bone mass, probably as a consequence of high bone turnover. We also identified a novel low-frequency non-coding variant with large effects on BMD near WNT16 (rs148771817(T), MAF = 1.2%, replication effect size = +0.41 s.d., Pmeta = 1 × 10(-11)). In general, there was an excess of association signals arising from deleterious coding and conserved non-coding variants. These findings provide evidence that low-frequency non-coding variants have large effects on BMD and fracture, thereby providing rationale for whole-genome sequencing and improved imputation reference panels to study the genetic architecture of complex traits and disease in the general population.

Heritability of bone mineral density (BMD) varies across skeletal sites, reflecting different relative contributions of genetic and environmental influences. To quantify the degree to which common genetic variants tag and environmental factors influence BMD, at different sites, we estimated the genetic (rg) and residual (re) correlations between BMD measured at the upper limbs (UL-BMD), lower limbs (LL-BMD) and skull (SK-BMD), using total-body DXA scans of ∼ 4,890 participants recruited by the Avon Longitudinal Study of Parents and their Children (ALSPAC). Point estimates of rg indicated that appendicular sites have a greater proportion of shared genetic architecture (LL-/UL-BMD rg = 0.78) between them, than with the skull (UL-/SK-BMD rg = 0.58 and LL-/SK-BMD rg = 0.43). Likewise, the residual correlation between BMD at appendicular sites (r(e) = 0.55) was higher than the residual correlation between SK-BMD and BMD at appendicular sites (r(e) = 0.20-0.24). To explore the basis for the observed differences in rg and re, genome-wide association meta-analyses were performed (n ∼ 9,395), combining data from ALSPAC and the Generation R Study identifying 15 independent signals from 13 loci associated at genome-wide significant level across different skeletal regions. Results suggested that previously identified BMD-associated variants may exert site-specific effects (i.e. differ in the strength of their association and magnitude of effect across different skeletal sites). In particular, variants at CPED1 exerted a larger influence on SK-BMD and UL-BMD when compared to LL-BMD (P = 2.01 × 10(-37)), whilst variants at WNT16 influenced UL-BMD to a greater degree when compared to SK- and LL-BMD (P = 2.31 × 10(-14)). In addition, we report a novel association between RIN3 (previously associated with Paget's disease) and LL-BMD (rs754388: β = 0.13, SE = 0.02, P = 1.4 × 10(-10)). Our results suggest that BMD at different skeletal sites is under a mixture of shared and specific genetic and environmental influences. Allowing for these differences by performing genome-wide association at different skeletal sites may help uncover new genetic influences on BMD.

Trps1 has been proposed as a candidate gene for a mouse bone mineral density (BMD) QTL on Chromosome (Chr) 15, but it remained unclear if this gene was associated with BMD in humans. We used newly available data and advanced bioinformatics techniques to confirm that Trps1 is the most likely candidate gene for the mouse QTL. In short, by combining the raw genetic mapping data from two F2 generation crosses of inbred strains of mice, we narrowed the 95% confidence interval of this QTL down to the Chr 15 region spanning from 6 to 24cM. This region contains 131 annotated genes. Using block haplotyping, all other genes except Trps1 were eliminated as candidates for this QTL. We then examined associations of 208 SNPs within 10kb of TRPS1 with BMD and hip geometry, using human genome-wide association study (GWAS) data from the GEFOS consortium. After correction for multiple testing, six TRPS1 SNPs were significantly associated with femoral neck BMD (P=0.0015-0.0019; adjusted P=0.038-0.048). We also found that three SNPs were highly associated with femoral neck width in women (rs10505257, P=8.6×10(-5), adjusted P=2.15×10(-3); rs7002384, P=5.5×10(-4), adjusted P=01.38×10(-2)). In conclusion, we demonstrated that combining association studies in humans with murine models provides an efficient strategy to identify new candidate genes for bone phenotypes., (Copyright © 2012 Elsevier Inc. All rights reserved.)

Overwhelming evidence has linked cardiovascular disease and osteoporosis, but the shared root cause of these two diseases of the elderly remains unknown. Low levels of high density lipoprotein cholesterol (HDL) and bone mineral density (BMD) are risk factors for cardiovascular disease and osteoporosis respectively. A number of correlation studies have attempted to determine if there is a relationship between serum HDL and BMD but these studies are confounded by a number of variables including age, diet, genetic background, gender and hormonal status. Collectively, these data suggest that there is a relationship between these two phenotypes, but that the nature of this relationship is context specific. Studies in mice plainly demonstrate that genetic loci for BMD and HDL co-map and transgenic mouse models have been used to show that a single gene can affect both serum HDL and BMD. Work completed to date has demonstrated that HDL can interact directly with both osteoblasts and osteoclasts, but no direct evidence links bone back to the regulation of HDL levels. Understanding the genetic relationship between BMD and HDL has huge implications for understanding the clinical relationship between CVD and osteoporosis and for the development of safe treatment options for both diseases., (Copyright © 2011 Elsevier Inc. All rights reserved.)

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., (© 2011 American Society for Bone and Mineral Research.)

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.

Unlabelled: Genetic analyses with mouse congenic strains for distal Chr1 have identified three closely linked QTLs regulating femoral vBMD, mid-diaphyseal cortical thickness, and trabecular microstructure in a sex-dependent fashion. The homologous relationship between distal mouse Chr 1 and human 1q21-24 offers the possibility of finding common regulatory genes for cortical and trabecular bone., Introduction: The distal third of mouse chromosome 1 (Chr 1) has been shown to carry a major quantitative trait locus (QTL) for BMD from several inbred mouse strain crosses. Genetic and functional analyses are essential to identify genes and cellular mechanisms for acquisition of peak bone mass., Materials and Methods: Nested congenic sublines of mice were developed with a C57BL/6J (B6) background carrying <1- to 9-Mbp-sized segments donated from C3H/HeJ (C3H). Isolated femurs from 16-wk-old female and male mice were measured by pQCT and microCT40 for volumetric (v)BMD, mid-diaphyseal cortical thickness, and distal trabecular phenotypes. Static and dynamic histomorphologic data were obtained on selected females and males at 16 wk., Results and Conclusions: We found that the original BMD QTL, Bmd5, mapped to distal Chr 1 consists of three QTLs with different effects on vBMD and trabecular bone in both sexes. Compared with B6 controls, femoral vBMD, BMD, and cortical thickness (p < 0.0001) were significantly increased in congenic subline females, but not in males, carrying C3H alleles at QTL-1. Both females and males carrying C3H alleles at QTL-1 showed marked increases in BV/TV by microCT compared with B6 mice (p < 0.0001). Females increased BV/TV by increasing trabecular thickness, whereas males increased trabecular number. In addition, the microCT40 data showed two unique QTLs for male trabecular bone, QTL-2 and QTL-3, which may interact to regulate trabecular thickness and number. These QTLs are closely linked with and proximal to QTL-1. The histomorphometric data revealed sex-specific differences in cellular and bone formation parameters. Mice and humans share genetic homology between distal mouse Chr 1 and human Chr 1q20-24 that is associated with adult human skeletal regulation. Sex- and compartment-specific regulatory QTLs in the mouse suggest the need to partition human data by sex to improve accuracy of mapping and genetic loci identification.

Unlabelled: Genetic analysis of an NZB/B1NJ x RF/J cross has identified QTLs for femur mechanical, geometric, and densitometric phenotypes. Most mechanical QTLs were associated with geometric QTLs, strongly suggesting common genetic regulation., Introduction: Previous studies have shown that bone architecture and BMD are important factors affecting bone strength, and both are genetically regulated. We conducted genetic analyses for loci regulating femur mechanical properties, geometric properties, and BMD in a cohort of F2 mice derived from intercross matings of (NZB/B1NJ x RF/J)F1 parents., Materials and Methods: Femurs were isolated from 662 10-week-old females. Mechanical properties were determined for a femur from each animal by three-point bending. Geometric properties and volumetric BMD (vBMD) were determined by pQCT. Genotype data were obtained by PCR assays for polymorphic markers carried in the genomic DNA of each mouse. Genome-wide scans were carried out for co-segregation of genetic marker data with values from 23 different phenotypes. Quantitative trait loci (QTLs) were identified for mechanical, geometric, and mineral density phenotypes., Results: QTLs for many phenotypes were significantly refined by covariate analyses using body weight and femur length. Major QTLs for mechanical and geometric phenotypes were found on chromosomes 5, 7, 9, 11, and 12. Nine chromosomal locations were identified with mechanical QTLs and 17 locations with one or more geometric QTLs. The significance of five mechanical and nine geometric QTLs was affected by the inclusion of covariates. These changes included both decreases and increases in significance. The QTLs on chromosomes 5 and 12 were decreased by inclusion of the covariates in the analysis, but QTLs on 7 and 11 were unaffected. Mechanical QTLs were almost always associated with geometric QTLs and less commonly (two of six) with vBMD QTLs., Conclusions: Genetic regulation of mechanical properties in the F(2) mice of this NZB/B1NJ x RF/J cross seems to be caused by genes regulating femur geometry.

The femurs from groups of mice from 29 different inbred strains were characterized to study the genetic variations in bone parameters. For these analyses, we used peripheral quantitative computed tomography to assess bone size and density in addition to three-point bend testing to assess bone mechanical properties. Highly significant differences between inbred strains were found for all size, density, and mechanical parameters measured (P < 0.0001). Correcting femoral cross-sectional geometry values or bone mechanical properties values for body weight or femur length reduced but did not eliminate the variations in bone geometry or bone mechanical properties. Mice of similar body size had as much as a 40% difference in the midshaft total area of the femur. Regression analysis suggested that 50.9% of the variation in maximum load among strains was related to variations in section modulus, i.e., cross-sectional geometry, 21.5% was related to variations in material bone density, and 27.7% to variations in quality. These components were further analyzed to show that 3.9-27.8% of the variation in maximum load was related to adaptation to mechanical stress. These findings indicate that there is a significant genetic variation in the femur cross-sectional area, density, and mechanical properties between inbred mouse strains. These studies identify inbred mouse strains suitable for future studies identifying genes regulating bone geometry and mechanical properties.

Targeted gene studies have demonstrated the importance of insulin-like growth factor-I (IGF-I) for osteoblast (OB) differentiation and the acquisition of peak bone mineral density (BMD). The skeletal response to allelic differences in IGF-I expression can also be measured in vivo, using congenic mice. We created a congenic strain with reduced (approximately 20%) circulating IGF-I (C3H.B6-6T [6T]) by backcrossing a small genomic region (30 cM) of Chromosome 6 (Chr6) from C3H/HeJ (C3H) onto a C57Bl/6J (B6) background. 6T female mice have lower serum IGF-I (P<0.001 vs. B6) but similar growth hormone (GH) and serum IGF binding protein (IGFBP) concentrations as B6. At 16 weeks of age, congenics have greater body fat (P<0.02 vs. B6) despite less total body weight, and exhibit smaller femoral cross-sectional size (P=0.001), reduced cortical thickness (P<0.001) and lower trabecular BV/TV (P<0.05) than B6. 6T mice also have suppressed serum leptin (P<0.01), but compared to B6 have similar markers of bone resorption (i.e., urine CTx and serum TRAP 5B). At 8 weeks of age, skeletal IGF-I mRNA from long bones was reduced by 40% (P<0.05) as were liver mRNA transcripts (i.e., 50%, P<0.01). Osteoblast progenitors from the bone marrow of 6T mice formed less colony forming unit fibroblasts by crystal violet staining than B6 (P<0.007) and had significantly reduced alkaline phosphatase-positive colonies than B6(P<0.0001). In addition, staining of bone marrow with oil red O revealed greater numbers of adipocytes in 6T than B6. Several candidate genes in the Chr6 QTL were excluded by lack of strain-related expression differences in bone, but genes positively regulating adipocyte differentiation including Alox 5 and PPAR-gamma require further study as either "pathway" or candidate genes. In summary, allelic differences in a QTL on Chr6 result in altered IGF-I gene expression, changes in OB lineage allocation, and reduced peak bone mass. Congenic mice are useful models not only for mapping genes related to bone mass but also for elucidating the biology underlying various skeletal phenotypes associated with more subtle manipulation of the mouse genome.

IGF-1 is a growth-promoting polypeptide that is essential for normal growth and development. In serum, the majority of the IGFs exist in a 150-kDa complex including the IGF molecule, IGF binding protein 3 (IGFBP-3), and the acid labile subunit (ALS). This complex prolongs the half-life of serum IGFs and facilitates their endocrine actions. Liver IGF-1-deficient (LID) mice and ALS knockout (ALSKO) mice exhibited relatively normal growth and development, despite having 75% and 65% reductions in serum IGF-1 levels, respectively. Double gene disrupted mice were generated by crossing LID+ALSKO mice. These mice exhibited further reductions in serum IGF-1 levels and a significant reduction in linear growth. The proximal growth plates of the tibiae of LID+ALSKO mice were smaller in total height as well as in the height of the proliferative and hypertrophic zones of chondrocytes. There was also a 10% decrease in bone mineral density and a greater than 35% decrease in periosteal circumference and cortical thickness in these mice. IGF-1 treatment for 4 weeks restored the total height of the proximal growth plate of the tibia. Thus, the double gene disruption LID+ALSKO mouse model demonstrates that a threshold concentration of circulating IGF-1 is necessary for normal bone growth and suggests that IGF-1, IGFBP-3, and ALS play a prominent role in the pathophysiology of osteoporosis.

Bone mineral density (BMD) is a strong predictor of osteoporotic fracture. It is also one of the most heritable disease-associated quantitative traits. As a result, there has been considerable effort focused on dissecting its genetic basis. Here, we performed a genome-wide association study (GWAS) in a panel of inbred strains to identify associations influencing BMD. This analysis identified a significant (P = 3.1 x 10−12) BMD locus on Chromosome 3@52.5 Mbp that replicated in two separate inbred strain panels and overlapped a BMD quantitative trait locus (QTL) previously identified in a F2 intercross. The association mapped to a 300 Kbp region containing four genes; Gm2447, Gm20750, Cog6, and Lhfp. Further analysis found that Lipoma HMGIC Fusion Partner (Lhfp) was highly expressed in bone and osteoblasts. Furthermore, its expression was regulated by a local expression QTL (eQTL), which overlapped the BMD association. A co-expression network analysis revealed that Lhfp was strongly connected to genes involved in osteoblast differentiation. To directly evaluate its role in bone, Lhfp deficient mice (Lhfp-/-) were created using CRISPR/Cas9. Consistent with genetic and network predictions, bone marrow stromal cells (BMSCs) from Lhfp-/- mice displayed increased osteogenic differentiation. Lhfp-/- mice also had elevated BMD due to increased cortical bone mass. Lastly, we identified SNPs in human LHFP that were associated (P = 1.2 x 10−5) with heel BMD. In conclusion, we used GWAS and systems genetics to identify Lhfp as a regulator of osteoblast activity and bone mass. [ABSTRACT FROM AUTHOR]

Cadherin-like and PC-esterase domain containing 1 ( CPED1 ) is an uncharacterized gene with no known function. Human genome wide association studies (GWAS) for bone mineral density (BMD) have repeatedly identified a significant locus on Chromosome 7 that contains the gene CPED1 , but it remains unclear if this gene could be causative. While an open reading frame for this gene has been predicted, there has been no systematic exploration of expression or alternate splicing for CPED1 in humans or mice.Using mouse models, we demonstrate that Cped1 is alternately spliced whereby transcripts are generated with exon 3 or exons 16 and 17 removed. In calvarial-derived pre-osteoblasts, Cped1 utilizes the predicted promoter upstream of exon 1 as well as alternate promoters upstream of exon 3 and exon 12.Lastly, we have determined that some transcripts terminate at the end of exon 10 and therefore do not contain the cadherin like and the PC esterase domains.Together, these data suggest that multiple protein products may be produced by this gene, with some products either lacking or containing both the predicted functional domains. Our data provide a framework upon which future functional studies will be built to understand the role of this gene in bone biology. [ABSTRACT FROM AUTHOR]

Osteoporosis is a complex metabolic bone disorder. Recently it has been appreciated that the "obesity in bone" phenomenon occurs at the expense of bone formation, and that is a key component of the pathology of this disease. Mouse models with altered bone expression levels of peroxisome proliferator-activated receptor gamma (PPARG) impact bone formation, but genetic studies connecting PPARG polymorphisms to skeletal phenotypes in humans have proven to be less than satisfactory. One missense polymorphism in exon one has been linked to low bone mineral density (BMD), but the most studied polymorphism, Pro12Ala, has not yet been examined in the context of skeletal phenotype. The studies to date are a promising start in leading to our understanding of the genetic contribution of PPARG to the phenotypes of BMD and fracture risk. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Genome-wide association studies (GWASs) for bone mineral density (BMD) in humans have identified over 1100 associations to date. However, identifying causal genes implicated by such studies has been challenging. Recent advances in the development of transcriptome reference datasets and computational approaches such as transcriptome-wide association studies (TWASs) and expression quantitative trait loci (eQTL) colocalization have proven to be informative in identifying putatively causal genes underlying GWAS associations. Here, we used TWAS/eQTL colocalization in conjunction with transcriptomic data from the Genotype-Tissue Expression (GTEx) project to identify potentially causal genes for the largest BMD GWAS performed to date. Using this approach, we identified 512 genes as significant using both TWAS and eQTL colocalization. This set of genes was enriched for regulators of BMD and members of bone relevant biological processes. To investigate the significance of our findings, we selected PPP6R3, the gene with the strongest support from our analysis which was not previously implicated in the regulation of BMD, for further investigation. We observed that Ppp6r3 deletion in mice decreased BMD. In this work, we provide an updated resource of putatively causal BMD genes and demonstrate that PPP6R3 is a putatively causal BMD GWAS gene. These data increase our understanding of the genetics of BMD and provide further evidence for the utility of combined TWAS/colocalization approaches in untangling the genetics of complex traits. [ABSTRACT FROM AUTHOR]

Abstract: Hepatic osteodystrophy (HOD) denotes the alterations in bone morphology and metabolism frequently observed in patients with chronic liver diseases, in particular in case of cholestatic conditions. The molecular mechanisms underlying HOD are only partially understood. In the present study, we characterized the bone phenotypes of the ATP-binding cassette transporter B4 knockout mouse (Abcb4 −/− ), a well-established mouse model of chronic cholestatic liver disease, with the aim of identifying and characterizing a mouse model for HOD. Furthermore, we investigated the influence of vitamin D on bone quality in this model. The bone morphology analyses revealed reduced bone mineral contents as well as changes in trabecular bone architecture and decreased cortical bone densities in Abcb4 −/− mice with severe liver fibrosis. We observed dysregulation of genes involved in bone remodeling (osteoprotegerin, osteocalcin, osteopontin) and vitamin D metabolism (7-dehydrocholesterol reductase, Gc-globulin, Cyp2r1, Cyp27a1) as well as alterations in calcium and vitamin D homeostasis. In addition, serum RANKL and TGF-β levels were increased in Abcb4 −/− mice. Vitamin D dietary intervention did not restore the bone phenotypes of Abcb4 −/− animals. We conclude that the Abcb4 −/− mouse provides an experimental framework and a preclinical model to gain further insights into the molecular pathobiology of HOD and to study the systemic effects of therapeutic interventions. [Copyright &y& Elsevier]