Meta-analysis of genome-wide studies identifies WNT16 and ESR1 SNPs associated with bone mineral density in premenopausal women

Muscle metabolic resilience and enhanced exercise adaptation by Esr1-induced remodeling of mitochondrial cristae-nucleoid architecture in males

Reduced estrogen action is associated with obesity and insulin resistance. However, the cell and tissue-specific actions of estradiol in maintaining metabolic health remain inadequately understood, especially in men. We observed that skeletal muscle ESR1/Esr1 (encodes estrogen receptor α [ERα]) is positively correlated with insulin sensitivity and metabolic health in humans and mice. Because skeletal muscle is a primary tissue involved in oxidative metabolism and insulin sensitivity, we generated muscle-selective Esr1 loss- and gain-of-expression mouse models. We determined that Esr1 links mitochondrial DNA replication and cristae-nucleoid architecture with metabolic function and insulin action in the skeletal muscle of male mice. Overexpression of human ERα in muscle protected male mice from diet-induced disruption of metabolic health and enhanced mitochondrial adaptation to exercise training intervention. Our findings indicate that muscle expression of Esr1 is critical for the maintenance of mitochondrial function and metabolic health in males and that tissue-selective activation of ERα can be leveraged to combat metabolic-related diseases in both sexes.

Regulation of Skeletogenic Pathways by m6A RNA Modification: A Comprehensive Review

In the complex process of skeletal development, the significance of m6A RNA methylation-a predominant form of RNA modification-has not been fully explored. This review discuss how m6A RNA methylation plays an important, though not yet fully understood, role in regulating skeletal formation. It examines how m6A influences key signaling pathways essential for skeletal development and homeostasis, suggesting various possible interactions between m6A methylation and these critical pathways. While the exact mechanisms for many of these interactions remain to be elucidated, m6A RNA methylation is anticipated to be a key emerging regulator in skeletal structure development across vertebrates. Highlighting the need for further research, this overview provides an in-depth look at the potential regulatory interactions of m6A RNA methylation within skeletal system. Uniquely, this review is the most comprehensive compilation of evidence linking components of m6A RNA methylation to signaling pathways involved in skeletogenesis.

Regulation of WNT16 in bone may involve upstream enhancers within CPED1

WNT16 stands up as an essential gene for bone homeostasis. Here, we present new evidence of the functional role of a particular region within WNT16. Performing 4 C chromatin conformation analysis in three osteoblast-related cells (the human fetal osteoblast hFOB 1.19 cell line, Saos 2 osteosarcoma cell line and mesenchymal Stem Cells -MSC-), we identify physical interactions between the proximal part of WNT16 intron 2, shown here to be an active promoter in Saos 2 osteosarcoma cells, and several putative regulatory regions within CPED1. Analysis of previously published RNA-seq data from hFOB cells disclosed low expression of a region located downstream of this promoter. Our results suggest a novel regulatory mechanism of WNT16 in bone, mediated by physical interaction with various enhancer regions within CPED1.

Loss of cped1 does not affect bone and lean tissue in zebrafish

Human genetic studies have nominated cadherin-like and PC-esterase domain-containing 1 (CPED1) as a candidate target gene mediating bone mineral density (BMD) and fracture risk heritability. Recent efforts to define the role of CPED1 in bone in mouse and human models have revealed complex alternative splicing and inconsistent results arising from gene targeting, making its function in bone difficult to interpret. To better understand the role of CPED1 in adult bone mass and morphology, we conducted a comprehensive genetic and phenotypic analysis of cped1 in zebrafish, an emerging model for bone and mineral research. We analyzed two different cped1 mutant lines and performed deep phenotyping to characterize more than 200 measures of adult vertebral, craniofacial, and lean tissue morphology. We also examined alternative splicing of zebrafish cped1 and gene expression in various cell/tissue types. Our studies fail to support an essential role of cped1 in adult zebrafish bone. Specifically, homozygous mutants for both cped1 mutant alleles, which are expected to result in loss-of-function and impact all cped1 isoforms, exhibited no significant differences in the measures examined when compared to their respective wildtype controls, suggesting that cped1 does not significantly contribute to these traits. We identified sequence differences in critical residues of the catalytic triad between the zebrafish and mouse orthologs of CPED1, suggesting that differences in key residues, as well as distinct alternative splicing, could underlie different functions of CPED1 orthologs in the two species. Our studies fail to support a requirement of cped1 in zebrafish bone and lean tissue, adding to evidence that variants at 7q31.31 can act independently of CPED1 to influence BMD and fracture risk.

The effects of endurance trainability phenotype, sex, and interval running training on bone collagen synthesis in adult rats

Bone is influenced by many factors such as genetics and mechanical loading, but the short-term physiological effects of these factors on bone (re)modelling are not well characterised. This study investigated the effects of endurance trainability phenotype, sex, and interval running training (7-week intervention) on bone collagen formation in rats using a deuterium oxide stable isotope tracer method. Bone samples of the femur diaphysis, proximal tibia, mid-shaft tibia, and distal tibia were collected after necropsy from forty-six 9 ± 3-month male and female rats selectively bred for yielding low (LRT) or high (HRT) responses to endurance training. Bone collagen proteins were isolated and hydrolysed, and fractional synthetic rates (FSRs) were determined by the incorporation of deuterium into protein-bound alanine via GC-pyrolysis-IRMS. There was a significant large main effect of phenotype at the femur site (p < 0.001; η2g = 0.473) with HRT rats showing greater bone collagen FSRs than LRT rats. There was a significant large main effect of phenotype (p = 0.008; η2g = 0.178) and a significant large main effect of sex (p = 0.005; η2g = 0.196) at the proximal site of the tibia with HRT rats showing greater bone collagen FSRs than LRT rats, and male rats showing greater bone collagen FSRs compared to female rats. There was a significant large main effect of training at the mid-shaft site of the tibia (p = 0.012; η2g = 0.159), with rats that underwent interval running training having greater bone collagen FSRs than control rats. Similarly, there was a significant large main effect of training at the distal site of the tibia (p = 0.050; η2g = 0.156), with rats in the interval running training group having greater bone collagen FSRs compared to rats in the control group. Collectively, this evidence highlights that bone responses to physiological effects are site-specific, indicating that interval running training has positive effects on bone collagen synthesis at the tibial mid-shaft and distal sites, whilst genetic factors affect bone collagen synthesis at the femur diaphysis (phenotype) and proximal tibia (phenotype and sex) in rats.

Causal relationship between reproductive factors and female bone density: a univariate and multivariate mendelian randomization study

Objective

Observational studies have found associations between reproductive factors and bone density in women. However, the causal relationships are not well understood. This study aims to investigate whether various reproductive factors are causally related to bone density at different skeletal sites using both univariable and multivariable Mendelian randomization (MR) methods.

Methods

The study incorporated four reproductive factors, namely, age at menarche (AAM), age at first live birth (AFB), age at menopause (ANM), and age at last live birth (ALB), as well as five distinct skeletal sites, including bone mineral density (BMD), heel calcaneus BMD, ultradistal forearm bone mineral density (FA-BMD), lumbar spine bone mineral density (LS-BMD), and femoral neck bone mineral density (FN-BMD). Univariable two-sample MR and multivariable MR analyses were conducted using data from published genome-wide association studies (GWASs). A total of 150 single nucleotide polymorphisms (SNPs) associated with the four reproductive factors were extracted from GWAS databases. The primary statistical analysis method utilized in this study was the inverse variance weighted (IVW) method.

Results

In the univariate MR analysis, we observed causal connections between four reproductive factors and bone density. Specifically, AAM had a significant impact on BMD and heel calcaneus BMD. Age at first live birth was negatively associated with FA-BMD. Age at last live birth showed a negative correlation with BMD and heel calcaneus BMD. ANM exhibited positive associations with BMD, heel calcaneus BMD, FA-BMD, and LS-BMD. Subsequently, we performed a multivariable MR analysis to examine the combined effects of multiple variables, which confirmed the persistence of associations between age at menopause and bone density at various sites. Additionally, we found a negative correlation between age at last live birth and heel calcaneus BMD.

Conclusion

This study offers a fresh perspective on the prevention of osteoporosis in women, explicitly stating that reproductive factors such as early menopause and late childbirth play a significant predictive role in individual bone density decline. Therefore, when developing osteoporosis screening and management protocols, reproductive factors should be included for a more comprehensive guidance of clinical practice.

Loss of cped1 does not affect bone and lean tissue in zebrafish

Human genetic studies have nominated Cadherin-like and PC-esterase Domain-containing 1 (CPED1) as a candidate target gene mediating bone mineral density (BMD) and fracture risk heritability. Recent efforts to define the role of CPED1 in bone in mouse and human models have revealed complex alternative splicing and inconsistent results arising from gene targeting, making its function in bone difficult to interpret. To better understand the role of CPED1 in adult bone mass and morphology, we conducted a comprehensive genetic and phenotypic analysis of cped1 in zebrafish, an emerging model for bone and mineral research. We analyzed two different cped1 mutant lines and performed deep phenotyping to characterize more than 200 measures of adult vertebral, craniofacial, and lean tissue morphology. We also examined alternative splicing of zebrafish cped1 and gene expression in various cell/tissue types. Our studies fail to support an essential role of cped1 in adult zebrafish bone. Specifically, homozygous mutants for both cped1 mutant alleles, which are expected to result in loss-of-function and impact all cped1 isoforms, exhibited no significant differences in the measures examined when compared to their respective wildtype controls, suggesting that cped1 does not significantly contribute to these traits. We identified sequence differences in critical residues of the catalytic triad between the zebrafish and mouse orthologs of CPED1, suggesting that differences in key residues, as well as distinct alternative splicing, could underlie different functions of CPED1 orthologs in the two species. Our studies fail to support a requirement of cped1 in zebrafish bone and lean tissue, adding to evidence that variants at 7q31.31 can act independently of CPED1 to influence BMD and fracture risk.

Large-scale circulating proteome association study (CPAS) meta-analysis identifies circulating proteins and pathways predicting incident hip fractures

Hip fractures are associated with significant disability, high cost, and mortality. However, the exact biological mechanisms underlying susceptibility to hip fractures remain incompletely understood. In an exploratory search of the underlying biology as reflected through the circulating proteome, we performed a comprehensive Circulating Proteome Association Study (CPAS) meta-analysis for incident hip fractures. Analyses included 6430 subjects from two prospective cohort studies (Cardiovascular Health Study and Trøndelag Health Study) with circulating proteomics data (aptamer-based 5 K SomaScan version 4.0 assay; 4979 aptamers). Associations between circulating protein levels and incident hip fractures were estimated for each cohort using age and sex-adjusted Cox regression models. Participants experienced 643 incident hip fractures. Compared with the individual studies, inverse-variance weighted meta-analyses yielded more statistically significant associations, identifying 23 aptamers associated with incident hip fractures (conservative Bonferroni correction 0.05/4979, P < 1.0 × 10-5). The aptamers most strongly associated with hip fracture risk corresponded to two proteins of the growth hormone/insulin growth factor system (GHR and IGFBP2), as well as GDF15 and EGFR. High levels of several inflammation-related proteins (CD14, CXCL12, MMP12, ITIH3) were also associated with increased hip fracture risk. Ingenuity pathway analysis identified reduced LXR/RXR activation and increased acute phase response signaling to be overrepresented among those proteins associated with increased hip fracture risk. These analyses identified several circulating proteins and pathways consistently associated with incident hip fractures. These findings underscore the usefulness of the meta-analytic approach for comprehensive CPAS in a similar manner as has previously been observed for large-scale human genetic studies. Future studies should investigate the underlying biology of these potential novel drug targets.

Potential Interaction between WNT16 and Vitamin D on Bone Qualities in Adolescent Idiopathic Scoliosis Patients and Healthy Controls

Adolescent idiopathic scoliosis (AIS) is a three-dimensional spinal deformity that is associated with low bone mineral density (BMD). Vitamin D (Vit-D) supplementation has been suggested to improve BMD in AIS, and its outcomes may be related to genetic factors. The present study aimed to (a) investigate the synergistic effect between a low BMD-related gene (wingless-related integration site 16, WNT16) and two important Vit-D pathway genes (Vit-D receptor, VDR, and Vit-D binding protein, VDBP) on serum Vit-D and bone qualities in Chinese AIS patients and healthy adolescents, and (b) to further investigate the effect of ablating Wnt16 on the cortical bone quality and whether diets with different dosages of Vit-D would further influence bone quality during the rapid growth phase in mice in the absence of Wnt16. A total of 519 girls (318 AIS vs. 201 controls) were recruited, and three selected single-nucleotide polymorphisms (SNPs) (WNT16 rs3801387, VDBP rs2282679, and VDR rs2228570) were genotyped. The serum 25(OH)Vit-D level was significantly associated with VDBP rs2282679 alleles (OR = -4.844; 95% CI, -7.521 to -2.167, p < 0.001). Significant multi-locus models were identified by generalized multifactor dimensionality reduction (GMDR) analyses on the serum 25(OH)Vit-D level (p = 0.006) and trabecular area (p = 0.044). In the gene-edited animal study, Wnt16 global knockout (KO) and wildtype (WT) male mice were provided with different Vit-D diets (control chow (1000 IU/Kg) vs. Vit-D-deficient chow (Nil in Vit-D) vs. high-dose Vit-D chow (20,000 IU/Kg)) from 4 weeks to 10 weeks old. Wnt16 global KO mice had significantly lower serum 25(OH)Vit-D levels and higher liver Vdbp mRNA expression levels than WT mice. In addition, Wnt16 global KO mice showed a decrease in bone density, cortical thickness and cortical area compared with WT mice. Interestingly, high-dose Vit-D chow led to lower bone density, cortical thickness, and cortical area in WT mice, which were less obvious in Wnt16 global KO mice. In conclusion, WNT16 may regulate the serum 25(OH)Vit-D level and bone qualities, which might be associated with VDBP expression. Further investigations with a larger sample size and wider spectrum of scoliosis severity are required to validate our findings regarding the interaction between WNT16 and Vit-D status in patients with AIS.

Gut microbiota impacts bone via Bacteroides vulgatus-valeric acid-related pathways

Although the gut microbiota has been reported to influence osteoporosis risk, the individual species involved, and underlying mechanisms, remain largely unknown. We performed integrative analyses in a Chinese cohort of peri-/post-menopausal women with metagenomics/targeted metabolomics/whole-genome sequencing to identify novel microbiome-related biomarkers for bone health. Bacteroides vulgatus was found to be negatively associated with bone mineral density (BMD), which was validated in US white people. Serum valeric acid (VA), a microbiota derived metabolite, was positively associated with BMD and causally downregulated by B. vulgatus. Ovariectomized mice fed B. vulgatus demonstrated increased bone resorption and poorer bone micro-structure, while those fed VA demonstrated reduced bone resorption and better bone micro-structure. VA suppressed RELA protein production (pro-inflammatory), and enhanced IL10 mRNA expression (anti-inflammatory), leading to suppressed maturation of osteoclast-like cells and enhanced maturation of osteoblasts in vitro. The findings suggest that B. vulgatus and VA may represent promising targets for osteoporosis prevention/treatment.

Multiple Mechanisms Explain Genetic Effects at the CPED1-WNT16 Bone Mineral Density Locus

PURPOSE OF REVIEW

Chromosome region 7q31.31, also known as the CPED1-WNT16 locus, is robustly associated with BMD and fracture risk. The aim of the review is to highlight experimental studies examining the function of genes at the CPED1-WNT16 locus.

RECENT FINDINGS

Genes that reside at the CPED1-WNT16 locus include WNT16, FAM3C, ING3, CPED1, and TSPAN12. Experimental studies in mice strongly support the notion that Wnt16 is necessary for bone mass and strength. In addition, roles for Fam3c and Ing3 in regulating bone morphology in vivo and/or osteoblast differentiation in vitro have been identified. Finally, a role for wnt16 in dually influencing bone and muscle morphogenesis in zebrafish has recently been discovered, which has brought forth new questions related to whether the influence of WNT16 in muscle may conspire with its influence in bone to alter BMD and fracture risk.

Association of the Estrogen Receptor 1 Polymorphisms rs2046210 and rs9383590 with the Risk, Age at Onset and Prognosis of Breast Cancer

Estrogen receptor α (ERα), encoded by the ESR1 gene, is a key prognostic and predictive biomarker firmly established in routine diagnostics and as a therapeutic target of breast cancer, and it has a central function in breast cancer biology. Genetic variants at 6q25.1, containing the ESR1 gene, were found to be associated with breast cancer susceptibility. The rs2046210 and rs9383590 single nucleotide variants (SNVs) are located in the same putative enhancer region upstream of ESR1 and were separately identified as candidate causal variants responsible for these associations. Here, both SNVs were genotyped in a hospital-based case-control study of 409 female breast cancer patients and 422 female controls of a Central European (Austrian) study population. We analyzed the association of both SNVs with the risk, age at onset, clinically and molecularly relevant characteristics and prognosis of breast cancer. We also assessed the concordances between both SNVs and the associations of each SNV conditional on the other SNV. The minor alleles of both SNVs were found to be non-significantly associated with an increased breast cancer risk. Significant associations were found in specific subpopulations, particularly in patients with an age younger than 55 years. The minor homozygotes of rs2046210 and the minor homozygotes plus heterozygotes of rs9383590 exhibited a several-years-younger age at onset than the common homozygotes, which was more pronounced in ER-positive and luminal patients. Importantly, the observed associations of each SNV were not consistently nullified upon correction for the other SNV nor upon analyses in common homozygotes for the other SNV. We conclude that both SNVs remain independent candidate causal variants.

Transcriptome-wide association study identifies novel genes associated with bone mineral density and lean body mass in children

OBJECTIVE

To identify novel candidate genes whose expression is associated with bone mineral density (BMD) and body lean mass (LM) in children.

METHODS

A tissue-specific transcriptome-wide association study (TWAS) was conducted utilizing a large-scale genome-wide association study (GWAS) dataset associated with BMD and LM and involving 10,414 participants. The measurement of BMD and LM phenotypes was made based on total-body dual-energy X-ray absorptiometry (TB-DXA) scans. TWAS was conducted by using FUSION software. Reference panels for muscle skeleton (MS), peripheral blood (NBL) and whole blood (YBL) were used for TWAS analysis. Functional enrichment and protein-protein interaction (PPI) analyses of the genes identified by TWAS were performed by using the online tool Metascape ( http://metascape.org ).

RESULTS

For BMD, we identified 174 genes with P < 0.05, such as IKZF1 (P = 1.46 × 10^-9) and CHKB (P = 8.31 × 10^-7). For LM, we identified 208 genes with P < 0.05, such as COPS5 (P = 3.03 × 10^-12) and MRPS33 (P = 5.45 × 10^-10). Gene ontology (GO) enrichment analysis of the BMD-associated genes revealed 200 GO terms, such as protein catabolic process (Log P = -5.09) and steroid hormone-mediated signaling pathway (Log P = -3.13). GO enrichment analysis of the LM-associated genes detected 287 GO terms, such as the apoptotic signaling pathway (Log P = -8.08) and lipid storage (Log P = -3.55).

CONCLUSION

This study identified several candidate genes for BMD and LM in children, providing novel clues to the genetic mechanisms underlying the development of childhood BMD and LM.

Identifying Causes of Fracture Beyond Bone Mineral Density: Evidence From Human Genetics

New therapies may help to prevent osteoporotic fractures other than through increasing bone mineral density (BMD). Because fracture risk has an important genetic component, we aim to identify loci increasing fracture risk that do not decrease BMD, using a recently-proposed structural equation model adapted to remove genetic influences of BMD on fracture risk. We used summary statistics of the largest genome-wide association studies (GWASs) for BMD and for fracture in these analyses. We next estimated the genetic correlation between the non-BMD or BMD-related genetic effects and other clinical risk factors for fracture. Last, based on white British participants in the UK Biobank, we conducted genetic risk score analyses to assess whether the aggregated genetic effects conferred increased major osteoporotic fracture risk. We found that only three loci affecting fracture risk exhibited genetic effects not mediated by BMD: SOST, CPED1-WNT16, and RSPO3, while these three loci simultaneously conferred BMD-related effects. No strong genetic associations between non-BMD or BMD-related effects and 16 clinical risk factors were observed. However, non-BMD effects might be genetic correlated with hip bone size. In the UK Biobank, a 1 standard deviation (1-SD) increase in the non-BMD genetic risk score conferred an odds ratio of 1.17 for incident major osteoporotic fracture, compared to 1.29 by a BMD-related genetic risk score. Our study suggests that the majority of common genetic predisposition toward fracture risk acts upon BMD. Although non-BMD genetic effects may exist, they are not strongly correlated with most traditional clinical risk factors. Risk loci harboring non-BMD genetic effects may influence other perspectives of bone quality, or confer effects that existing GWASs fail to capture, but they demonstrate weaker impact on fracture risk than BMD-related genetic effects. These findings suggest that most successful drug development programs for osteoporosis should focus on pathways identified through BMD-associated loci. © 2022 American Society for Bone and Mineral Research (ASBMR).

Associations between WNT signaling pathway-related gene polymorphisms and risks of osteoporosis development in Chinese postmenopausal women: a case-control study

BACKGROUND

The WNT signaling pathway is involved in the regulation of bone homeostasis, and the effect of WNT signaling pathway-related gene (WNT16 and LRP5) polymorphisms on osteoporosis risk among Chinese postmenopausal women is still unknown. Hence, we performed a case-control study to assess the association of WNT signaling pathway-related gene polymorphisms and osteoporosis risk.

METHODS

A total of 1026 women (515 osteoporosis patients and 511 controls) of postmenopausal age who were randomly sampled from Xi'an 630 Hospital (Shaanxi Province, China) were involved in this study. Seven genetic polymorphisms in WNT16 (rs3779381, rs3801387, rs917727 and rs7776725) and LRP5 (rs2291467, rs11228240 and rs12272917) were selected and genotyped using the Agena MassARRAY iPLEX system. The association of the genetic polymorphisms and osteoporosis risk was assessed by odds ratios and 95% confidence intervals. The multifactor dimensionality reduction (MDR) method was conducted to analyze single nucleotide polymorphism (SNP)-SNP interaction.

RESULTS

We found that LRP5 polymorphisms (rs2291467, rs11228240 and rs12272917) were significantly associated with a decreased risk of osteoporosis in homozygote, recessive and additive models (p < 0.05). Stratification analysis showed that LRP5 polymorphisms (rs2291467, rs11228240 and rs12272917) significantly decreased the osteoporosis risk in the subgroup of body mass index (BMI) ≤ 24 (p < 0.05) and that individuals carrying a heterozygote genotype of WNT16 polymorphisms (rs3779381, rs3801387, rs917727 and rs7776725) had a higher osteoporosis risk in the subgroup of BMI > 24 (p < 0.05). Two haplotypes (haplotype 1: rs3779381, rs3801387, rs917727 and rs7776725; haplotype 2: rs2291467 and rs11228240) were observed, yet only Trs2291467Trs11228240 and Crs2291467Crs11228240 had a strong association with a decreased risk of osteoporosis (p < 0.05). Additionally, MDR analysis revealed that LRP5 rs2291467 was the best model in single-locus MDR analysis. A seven-locus model including rs3779381-AG, rs7776725-TC, rs3801387-GA and rs917727-TC in WNT16 and rs11228240-CC, rs12272917-TC and rs2291467-CC in LRP5 was the best model in multiple-loci MDR analysis (p < 0.001). These two best models were the most significantly associated with osteoporosis risk.

CONCLUSIONS

Our findings suggested that WNT16 and LRP5 genetic polymorphisms are associated with osteoporosis risk among Chinese postmenopausal women.

Wnt Pathway Extracellular Components and Their Essential Roles in Bone Homeostasis

The Wnt pathway is involved in several processes essential for bone development and homeostasis. For proper functioning, the Wnt pathway is tightly regulated by numerous extracellular elements that act by both activating and inhibiting the pathway at different moments. This review aims to describe, summarize and update the findings regarding the extracellular modulators of the Wnt pathway, including co-receptors, ligands and inhibitors, in relation to bone homeostasis, with an emphasis on the animal models generated, the diseases associated with each gene and the bone processes in which each member is involved. The precise knowledge of all these elements will help us to identify possible targets that can be used as a therapeutic target for the treatment of bone diseases such as osteoporosis.

Wnt16 signaling in bone homeostasis and osteoarthristis

Wnts are secreted cysteine-rich glycoproteins involved in joint development and skeletal homeostasis and have been implicated in the occurrence of osteoarthritis. Over the past decade, Wnt16, a member of the Wnt family, has received widespread attention for its strong association with bone mineral density, cortical bone thickness, bone strength, and osteoporotic fracture risk. In recent years, further studies have shed light on the role of Wnt16 a positive regulator of bone mass and protective regulator of osteoarthritis progression. Transduction mechanisms and crosstalk involving Wnt16 signaling have also been illustrated. More importantly, local Wnt16 treatment has been shown to ease osteoarthritis, inhibit bone resorption, and promote new bone formation in bone defect models. Thus, Wnt16 is now a potential therapeutic target for skeletal diseases and osteoarthritis. This paper reviews our current understanding of the mechanisms by which Wnt16 signaling regulates bone homeostasis and osteoarthritis.

WNT16 elevation induced cell senescence of osteoblasts in ankylosing spondylitis

BACKGROUND

WNT16 is critical for bone homeostasis, but the effect of WNT16 in ankylosing spondylitis (AS) is still unknown. Here, we investigated whether WNT16 influences bone formation and pathophysiological changes of AS in an in vitro model.

METHODS

The bone tissue from the facet joints was obtained from seven disease control and seven AS patients. Primary osteoprogenitor cells of the facet joints were isolated using an outgrowth method. Isolated osteoprogenitor cells from both control and AS tissues were analyzed by microarray, RT-qPCR, immunoblotting, and immunohistochemistry. The bone-forming activity of osteoprogenitor cells was assessed by various in vitro assays. β-galactosidase staining and senescence-associated secretory phenotype (SASP) using RT-qPCR were used to assess cell senescence.

RESULTS

In microarray analysis, WNT16 expression was significantly elevated in AS osteoprogenitor cells compared to the control. We also validated that WNT16 expression was elevated in AS-osteoprogenitor cells and human AS-bone tissues. WNT16 treatment inhibited bone formation in AS-osteoprogenitor cells but not in the control. Intriguingly, AS-osteoprogenitor cells were stained markedly with β-galactosidase for cell senescence in WNT16 treatment. Furthermore, in an H2O2 stress-induced premature senescence condition, WNT16 treatment increased cell senescence in AS-osteoprogenitor cells and WNT16 treatment under the H2O2 stress condition showed an increase in p21 protein and SASP mRNA expression. The WNT16-induced SASP expression in AS-osteoprogenitor cells was reduced in WNT16 knockdown cultures.

CONCLUSION

WNT16 is highly expressed in AS and WNT16 treatment facilitated cell senescence in AS-osteoprogenitor cells during osteoblast differentiation accompanied by suppression of bone formation. The identified role of WNT16 in AS could influence bone loss in AS patients.

Using "-omics" Data to Inform Genome-wide Association Studies (GWASs) in the Osteoporosis Field

PURPOSE OF REVIEW

Osteoporosis constitutes a major societal health problem. Genome-wide association studies (GWASs) have identified over 1100 loci influencing bone mineral density (BMD); however, few of the causal genes have been identified. Here, we review approaches that use "-omics" data and genetic- and systems genetics-based analytical strategies to facilitate causal gene discovery.

RECENT FINDINGS

The bone field is beginning to adopt approaches that are commonplace in other disease disciplines. The slower progress has been due in part to the lack of large-scale "omics" data on bone and bone cells. This is however changing, and approaches such as eQTL colocalization, transcriptome-wide association studies (TWASs), network, and integrative approaches are beginning to provide significant insight into the genes responsible for BMD GWAS associations. The use of "-omics" data to inform BMD GWASs has increased in recent years, leading to the identification of novel regulators of BMD in humans. The ultimate goal will be to use this information to develop more effective therapies to treat and ultimately prevent osteoporosis.

Level and change in bone microarchitectural parameters and their relationship with previous fracture and established bone mineral density loci

BACKGROUND

Osteoporosis is characterised by a reduction of bone mineral density (BMD) and predisposition to fracture. Bone microarchitecture, measured by high resolution peripheral quantitative computed tomography (HR-pQCT), has been related to fragility fractures and BMD and has been the subject of large-scale genome-wide analysis. We investigated whether fracture was related to baseline values and longitudinal changes in bone microarchitecture and whether bone microarchitecture was associated with established BMD loci.

METHODS

115 males and 99 females (aged 72-81 at baseline) from the Hertfordshire Cohort Study (HCS) were analysed. Fracture history was determined in 2011-2012 by self-report and vertebral fracture assessment. Participants underwent HR-pQCT scans of the distal radius and tibia in 2011-2012 and 2017. Previous fracture in relation to baseline values and changes in tibial HR-pQCT parameters was examined using sex-adjusted logistic regression with and without adjustment for age, sociodemographic, lifestyle and clinical characteristics; baseline values and changes in parameters associated with previous fracture were then examined in relation to four established BMD loci after adjustment for sex and age.

RESULTS

Previous fracture was related to: higher trabecular area (fully-adjusted odds ratio [95% CI] per SD greater baseline value: 2.18 [1.27,3.73], p = 0.005); lower total volumetric BMD (0.53 [0.34,0.84], p = 0.007), cortical area (0.53 [0.30,0.95], p = 0.032), cortical BMD (0.56 [0.36,0.88], p = 0.011) and cortical thickness (0.45 [0.27,0.77], p = 0.004); and greater declines in trabecular BMD (p = 0.001). Associations were robust in sex- and fully-adjusted analysis. Relationships between BMD loci and these HR-pQCT parameters were weak: rs3801387 (WNT16) was related to decline in trabecular BMD (p = 0.011) but no other associations were significant (p > 0.05).

CONCLUSION

Baseline values of HR-pQCT parameters and greater decline in trabecular BMD were associated with fracture. Change in trabecular BMD was associated with WNT16 which has been demonstrated to influence bone health in murine models and human genome-wide association studies (GWAS).

Twelve years of GWAS discoveries for osteoporosis and related traits: advances, challenges and applications

Osteoporosis is a common skeletal disease, affecting ~200 million people around the world. As a complex disease, osteoporosis is influenced by many factors, including diet (e.g. calcium and protein intake), physical activity, endocrine status, coexisting diseases and genetic factors. In this review, we first summarize the discovery from genome-wide association studies (GWASs) in the bone field in the last 12 years. To date, GWASs and meta-analyses have discovered hundreds of loci that are associated with bone mineral density (BMD), osteoporosis, and osteoporotic fractures. However, the GWAS approach has sometimes been criticized because of the small effect size of the discovered variants and the mystery of missing heritability, these two questions could be partially explained by the newly raised conceptual models, such as omnigenic model and natural selection. Finally, we introduce the clinical use of GWAS findings in the bone field, such as the identification of causal clinical risk factors, the development of drug targets and disease prediction. Despite the fruitful GWAS discoveries in the bone field, most of these GWAS participants were of European descent, and more genetic studies should be carried out in other ethnic populations to benefit disease prediction in the corresponding population.

Wnt16 Elicits a Protective Effect Against Fractures and Supports Bone Repair in Zebrafish

Bone homeostasis is a dynamic, multicellular process that is required throughout life to maintain bone integrity, prevent fracture, and respond to skeletal damage. WNT16 has been linked to bone fragility and osteoporosis in human genome wide‐association studies, as well as the functional hematopoiesis of leukocytes in vivo. However, the mechanisms by which WNT16 promotes bone health and repair are not fully understood. In this study, CRISPR‐Cas9 was used to generate mutant zebrafish lacking Wnt16 (wnt16^−/−) to study its effect on bone dynamically. The wnt16 mutants displayed variable tissue mineral density (TMD) and were susceptible to spontaneous fractures and the accumulation of bone calluses at an early age. Fractures were induced in the lepidotrichia of the caudal fins of wnt16^−/− and WT zebrafish; this model was used to probe the mechanisms by which Wnt16 regulates skeletal and immune cell dynamics in vivo. In WT fins, wnt16 expression increased significantly during the early stages for bone repair. Mineralization of bone during fracture repair was significantly delayed in wnt16 mutants compared with WT zebrafish. Surprisingly, there was no evidence that the recruitment of innate immune cells to fractures or soft callus formation was altered in wnt16 mutants. However, osteoblast recruitment was significantly delayed in wnt16 mutants postfracture, coinciding with precocious activation of the canonical Wnt signaling pathway. In situ hybridization suggests that canonical Wnt‐responsive cells within fractures are osteoblast progenitors, and that osteoblast differentiation during bone repair is coordinated by the dynamic expression of runx2a and wnt16. This study highlights zebrafish as an emerging model for functionally validating osteoporosis–associated genes and investigating fracture repair dynamically in vivo. Using this model, it was found that Wnt16 protects against fracture and supports bone repair, likely by modulating canonical Wnt activity via runx2a to facilitate osteoblast differentiation and bone matrix deposition. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC. on behalf of American Society for Bone and Mineral Research.

Factors influencing peak bone mass gain

Bone mass is a key determinant of osteoporosis and fragility fractures. Epidemiologic studies have shown that a 10% increase in peak bone mass (PBM) at the population level reduces the risk of fracture later in life by 50%. Low PBM is possibly due to the bone loss caused by various conditions or processes that occur during adolescence and young adulthood. Race, gender, and family history (genetics) are responsible for the majority of PBM, but other factors, such as physical activity, calcium and vitamin D intake, weight, smoking and alcohol consumption, socioeconomic status, age at menarche, and other secondary causes (diseases and medications), play important roles in PBM gain during childhood and adolescence. Hence, the optimization of lifestyle factors that affect PBM and bone strength is an important strategy to maximize PBM among adolescents and young people, and thus to reduce the low bone mass or osteoporosis risk in later life. This review aims to summarize the available evidence for the common but important factors that influence bone mass gain during growth and development and discuss the advances of developing high PBM.

Osteoblast-Osteoclast Communication and Bone Homeostasis

Bone remodeling is tightly regulated by a cross-talk between bone-forming osteoblasts and bone-resorbing osteoclasts. Osteoblasts and osteoclasts communicate with each other to regulate cellular behavior, survival and differentiation through direct cell-to-cell contact or through secretory proteins. A direct interaction between osteoblasts and osteoclasts allows bidirectional transduction of activation signals through EFNB2-EPHB4, FASL-FAS or SEMA3A-NRP1, regulating differentiation and survival of osteoblasts or osteoclasts. Alternatively, osteoblasts produce a range of different secretory molecules, including M-CSF, RANKL/OPG, WNT5A, and WNT16, that promote or suppress osteoclast differentiation and development. Osteoclasts also influence osteoblast formation and differentiation through secretion of soluble factors, including S1P, SEMA4D, CTHRC1 and C3. Here we review the current knowledge regarding membrane bound- and soluble factors governing cross-talk between osteoblasts and osteoclasts.

WNT16 Requires Gα Subunits as Intracellular Partners for Both Its Canonical and Non-Canonical WNT Signalling Activity in Osteoblasts

In the past years, WNT16 became an interesting target in the field of skeletal research, as it was identified as an essential regulator of the cortical bone compartment, with the ability to increase both cortical and trabecular bone mass and strength in vivo. Even though there are indications that these advantageous effects are coming from canonical and non-canonical WNT-signalling activity, a clear model of WNT signalling by WNT16 is not yet depicted. We, therefore, investigated the modulation of canonical (WNT/β-catenin) and non-canonical [WNT/calcium, WNT/planar cell polarity (PCP)] signalling in human embryonic kidney (HEK) 293 T and SaOS2 cells. Here, we demonstrated that WNT16 activates all WNT-signalling pathways in osteoblasts, whereas only WNT/calcium signalling was activated in HEK293T cells. In osteoblasts, we therefore, additionally investigated the role of Gα subunits as intracellular partners in WNT16's mechanism of action by performing knockdown of Gα12, Gα13 and Gαq. These studies point out that the above-mentioned Gα subunits might be involved in the WNT/β-catenin and WNT/calcium-signalling activity by WNT16 in osteoblasts, and for Gα12 in its WNT/PCP-signalling activity, illustrating a novel possible mechanism of interplay between the different WNT-signalling pathways in osteoblasts. Additional studies are needed to demonstrate whether this mechanism is specific for WNT16 signalling or relevant for all other WNT ligands as well. Altogether, we further defined WNT16's mechanism of action in osteoblasts that might underlie the well-known beneficial effects of WNT16 on skeletal homeostasis. These findings on WNT16 and the activity of specific Gα subunits in osteoblasts could definitely contribute to the development of novel therapeutic approaches for fragility fractures in the future.

A road map for understanding molecular and genetic determinants of osteoporosis

Osteoporosis is a highly prevalent disorder characterized by low bone mineral density and an increased risk of fracture, termed osteoporotic fracture. Notably, bone mineral density, osteoporosis and osteoporotic fracture are highly heritable; however, determining the genetic architecture, and especially the underlying genomic and molecular mechanisms, of osteoporosis in vivo in humans is still challenging. In addition to susceptibility loci identified in genome-wide association studies, advances in various omics technologies, including genomics, transcriptomics, epigenomics, proteomics and metabolomics, have all been applied to dissect the pathogenesis of osteoporosis. However, each technology individually cannot capture the entire view of the disease pathology and thus fails to comprehensively identify the underlying pathological molecular mechanisms, especially the regulatory and signalling mechanisms. A change to the status quo calls for integrative multi-omics and inter-omics analyses with approaches in 'systems genetics and genomics'. In this Review, we highlight findings from genome-wide association studies and studies using various omics technologies individually to identify mechanisms of osteoporosis. Furthermore, we summarize current studies of data integration to understand, diagnose and inform the treatment of osteoporosis. The integration of multiple technologies will provide a road map to illuminate the complex pathogenesis of osteoporosis, especially from molecular functional aspects, in vivo in humans.

The genetic architecture of osteoporosis and fracture risk

Osteoporosis and fracture risk are common complex diseases, caused by an interaction of numerous disease susceptibility genes and environmental factors. With the advances in genomic technologies, large-scale genome-wide association studies (GWAS) have been performed which have broadened our understanding of the genetic architecture and biological mechanisms of complex disease. Currently, more than ~90 loci have been found associated with DXA derived bone mineral density (BMD), over ~500 loci with heel estimated BMD and several others with other less widely available bone parameters such as bone geometry, shape, and microarchitecture. Notably, several of the pathways identified by the GWAS efforts correspond to pathways that are currently targeted for the treatment of osteoporosis. Overall, tremendous progress in the field of the genetics of osteoporosis has been achieved with the discovery of WNT16, EN1, DAAM2, and GPC6 among others. Assessment of the function and biological mechanisms of the remaining genes may further untangle the complex genetic landscape of osteoporosis and fracture risk. With this review we aimed to provide a general overview of the existing GWAS studies on osteoporosis traits and fracture risk.

Meta-Analysis of Genomewide Association Studies Reveals Genetic Variants for Hip Bone Geometry

Hip geometry is an important predictor of fracture. We performed a meta-analysis of GWAS studies in adults to identify genetic variants that are associated with proximal femur geometry phenotypes. We analyzed four phenotypes: (i) femoral neck length; (ii) neck-shaft angle; (iii) femoral neck width, and (iv) femoral neck section modulus, estimated from DXA scans using algorithms of hip structure analysis. In the Discovery stage, 10 cohort studies were included in the fixed-effect meta-analysis, with up to 18,719 men and women ages 16 to 93 years. Association analyses were performed with ∼2.5 million polymorphisms under an additive model adjusted for age, body mass index, and height. Replication analyses of meta-GWAS significant loci (at adjusted genomewide significance [GWS], threshold p ≤ 2.6 × 10-8 ) were performed in seven additional cohorts in silico. We looked up SNPs associated in our analysis, for association with height, bone mineral density (BMD), and fracture. In meta-analysis (combined Discovery and Replication stages), GWS associations were found at 5p15 (IRX1 and ADAMTS16); 5q35 near FGFR4; at 12p11 (in CCDC91); 11q13 (near LRP5 and PPP6R3 (rs7102273)). Several hip geometry signals overlapped with BMD, including LRP5 (chr. 11). Chr. 11 SNP rs7102273 was associated with any-type fracture (p = 7.5 × 10-5 ). We used bone transcriptome data and discovered several significant eQTLs, including rs7102273 and PPP6R3 expression (p = 0.0007), and rs6556301 (intergenic, chr.5 near FGFR4) and PDLIM7 expression (p = 0.005). In conclusion, we found associations between several genes and hip geometry measures that explained 12% to 22% of heritability at different sites. The results provide a defined set of genes related to biological pathways relevant to BMD and etiology of bone fragility. © 2019 American Society for Bone and Mineral Research.

Advanced Genetic Approaches in Discovery and Characterization of Genes Involved With Osteoporosis in Mouse and Human

Osteoporosis is a complex condition with contributions from, and interactions between, multiple genetic loci and environmental factors. This review summarizes key advances in the application of genetic approaches for the identification of osteoporosis susceptibility genes. Genome-wide linkage analysis (GWLA) is the classical approach for identification of genes that cause monogenic diseases; however, it has shown limited success for complex diseases like osteoporosis. In contrast, genome-wide association studies (GWAS) have successfully identified over 200 osteoporosis susceptibility loci with genome-wide significance, and have provided most of the candidate genes identified to date. Phenome-wide association studies (PheWAS) apply a phenotype-to-genotype approach which can be used to complement GWAS. PheWAS is capable of characterizing the association between osteoporosis and uncommon and rare genetic variants. Another alternative approach, whole genome sequencing (WGS), will enable the discovery of uncommon and rare genetic variants in osteoporosis. Meta-analysis with increasing statistical power can offer greater confidence in gene searching through the analysis of combined results across genetic studies. Recently, new approaches to gene discovery include animal phenotype based models such as the Collaborative Cross and ENU mutagenesis. Site-directed mutagenesis and genome editing tools such as CRISPR/Cas9, TALENs and ZNFs have been used in functional analysis of candidate genes in vitro and in vivo. These resources are revolutionizing the identification of osteoporosis susceptibility genes through the use of genetically defined inbred mouse libraries, which are screened for bone phenotypes that are then correlated with known genetic variation. Identification of osteoporosis-related susceptibility genes by genetic approaches enables further characterization of gene function in animal models, with the ultimate aim being the identification of novel therapeutic targets for osteoporosis.

Molecular signaling in bone cells: Regulation of cell differentiation and survival

The achievement of proper bone mass and architecture, and their maintenance throughout life requires the concerted actions of osteoblasts, the bone forming cells, and osteoclasts, the bone resorbing cells. The differentiation and activity of osteoblasts and osteoclasts are regulated by molecules produced by matrix-embedded osteocytes, as well as by cross talk between osteoblasts and osteoclasts through secreted factors. In addition, it is likely that direct contact between osteoblast and osteoclast precursors, and the contact of these cells with osteocytes and cells in the bone marrow, also modulates bone cell differentiation and function. With the advancement of molecular and genetic tools, our comprehension of the intracellular signals activated in bone cells has evolved significantly, from early suggestions that osteoblasts and osteoclasts have common precursors and that osteocytes are inert cells in the bone matrix, to the very sophisticated understanding of a network of receptors, ligands, intracellular kinases/phosphatases, transcription factors, and cell-specific genes that are known today. These advances have allowed the design and FDA-approval of new therapies to preserve and increase bone mass and strength in a wide variety of pathological conditions, improving bone health from early childhood to the elderly. We have summarized here the current knowledge on selected intracellular signal pathways activated in osteoblasts, osteocytes, and osteoclasts.

Integrative genomic analysis predicts novel functional enhancer-SNPs for bone mineral density

Osteoporosis is a skeletal disorder characterized by low bone mineral density (BMD) and deterioration of bone microarchitecture. To identify novel genetic loci underlying osteoporosis, an effective strategy is to focus on scanning of variants with high potential functional impacts. Enhancers play a crucial role in regulating cell-type-specific transcription. Therefore, single-nucleotide polymorphisms (SNPs) located in enhancers (enhancer-SNPs) may represent strong candidate functional variants. Here, we performed a targeted analysis for potential functional enhancer-SNPs that may affect gene expression and biological processes in bone-related cells, specifically, osteoblasts, and peripheral blood monocytes (PBMs), using five independent cohorts (n = 5905) and the genetics factors for osteoporosis summary statistics, followed by comprehensive integrative genomic analyses of chromatin states, transcription, and metabolites. We identified 15 novel enhancer-SNPs associated with femoral neck and lumbar spine BMD, including 5 SNPs mapped to novel genes (e.g., rs10840343 and rs10770081 in IGF2 gene) and 10 novel SNPs mapped to known BMD-associated genes (e.g., rs2941742 in ESR1 gene, and rs10249092 and rs4342522 in SHFM1 gene). Interestingly, enhancer-SNPs rs10249092 and rs4342522 in SHFM1 were tightly linked, but annotated to different enhancers in PBMs and osteoblasts, respectively, suggesting that even tightly linked SNPs may regulate the same target gene and contribute to the phenotype variation in cell-type-specific manners. Importantly, ten enhancer-SNPs may also regulate BMD variation by affecting the serum metabolite levels. Our findings revealed novel susceptibility loci that may regulate BMD variation and provided intriguing insights into the genetic mechanisms of osteoporosis.

The Aging Skeleton

Skeletal aging begins after peak bone mass is reached; progressive bone loss then occurs. Peak bone mass may occur at different ages in different skeletal sites and varies between sexes. Accelerated loss of bone occurs in the perimenopausal period in women, whereas more gradual but progressive loss of bone occurs in aging men. Changes in bone quality as well as bone quantity occur during growth and subsequent aging. These include changes in bone microarchitecture which may differ between cortical and trabecular compartments and in different sites, and may impact on bone size and geometry. Changes in material properties of bone matrix may also occur with aging. Loss of bone quantity and altered bone quality with aging may weaken bones and culminate in osteoporosis with an increased risk of fractures. Both genetic and epigenetic mechanisms may predispose to osteoporosis. Cellular and molecular events underlie the alterations in bone quantity and quality. Osteoclastic bone resorption and osteoblastic bone formation, tightly regulated by hormones, growth factors, and cytokines, are organized in coordinated activities resulting in remodeling and modeling. Malignancies, and anti-neoplastic therapies, may impact on the cellular and molecular events in the aging skeleton and produce focal or diffuse skeletal lesions and fractures.

Genetic Predisposition for Osteoporosis and Fractures in Postmenopausal Women

Osteoporosis is a disease with complex etiology where the genetic factors may account for as much as 50-85% of the risk of its development in postmenopausal women. The polymorphism of estrogen receptor genes (ESR1, ESR2) seems essential among the genetic factors. The goal of this study was to analyze polymorphisms of selected genes in a population of postmenopausal women treated for osteoporosis and to evaluate the influence of genetic and nongenetic factors on the estimated 10-year risk of fracture. The study group consisted of 214 women hospitalized for treatment of postmenopausal osteoporosis. We investigated the presence of ESR1, ESR2, LRP5, and WNT16 genetic polymorphisms and the risk of fracture in each woman. The main finding was that of significant differences in the polymorphisms of the WNT16 rs2908004 genetic variant, notably, the less frequent presence of TC allele in women with a greater risk of osteoporotic fractures. We conclude that the polymorphism of the WNT16 gene seems highly relevant in the pathogenesis of osteoporosis, which makes it a promising object for further research on the genetic background of fracture risk.

Overexpression of WNT16 Does Not Prevent Cortical Bone Loss Due to Glucocorticoid Treatment in Mice

Glucocorticoids (GC) are commonly used for the treatment of a wide variety of autoimmune, pulmonary, gastrointestinal, and malignancy conditions. One of the devastating side effects of GC use is osteoporotic fractures, particularly in the spine and hip. Bisphosphonates (BP) are the most commonly prescribed pharmacological agents for the prevention and treatment of GC-induced osteoporosis (GIO). However, GIO is marked by reduced bone formation and BP serves mainly to decrease bone resorption. The WNT signaling pathway plays a major role in bone and mineral homeostasis. Previously, we demonstrated that overexpression of WNT16 in mice led to higher bone mineral density and improved bone microarchitecture and strength. We hypothesized that WNT16 overexpression would prevent bone loss due to glucocorticoid treatment in mice. To test our hypothesis, we treated adult wild-type and WNT16-transgenic mice with vehicle and GC (prednisolone; 2.1 mg/kg body weight) via slow-release pellets for 28 days. We measured bone mass and microarchitecture by dual-energy X-ray absorptiometry (DXA) and micro-CT, and performed gene expression and serum biochemical analysis. We found that GC treatment compared with the vehicle significantly decreased femoral areal bone mineral density (aBMD), bone mineral content (BMC), and cortical bone area and thickness in both wild-type and transgenic female mice. In contrast, the trabecular bone parameters at distal femur were not significantly changed by GC treatment in male and female mice for both genotypes. Further, we observed significantly lower level of serum P1NP and a tendency of higher level of serum TRAP in wild-type and transgenic mice due to GC treatment in both sexes. Gene expression analysis showed lower mRNA levels of Wnt16, Opg, and Opg/Rankl ratio in GC-treated female mice for both genotypes compared with the sex-matched vehicle-treated mice. These data suggest that although WNT16 overexpression resulted in higher baseline bone mineral density and bone volume per trabecular volume (BV/TV) in the transgenic mice, this was insufficient to prevent bone loss in mice due to glucocorticoid treatment.

Characterization of expression and alternative splicing of the gene cadherin-like and PC esterase domain containing 1 (Cped1)

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.

Meta-Analysis of Genome-Wide Association Studies Identifies Novel Functional CpG-SNPs Associated with Bone Mineral Density at Lumbar Spine

Osteoporosis is a serious public health issue, which is mostly characterized by low bone mineral density (BMD). To search for additional genetic susceptibility loci underlying BMD variation, an effective strategy is to focus on testing of specific variants with high potential of functional effects. Single nucleotide polymorphisms (SNPs) that introduce or disrupt CpG dinucleotides (CpG-SNPs) may alter DNA methylation levels and thus represent strong candidate functional variants. Here, we performed a targeted GWAS for 63,627 potential functional CpG-SNPs that may affect DNA methylation in bone-related cells, in five independent cohorts (n = 5905). By meta-analysis, 9 CpG-SNPs achieved a genome-wide significance level (p < 7.86 × 10-7) for association with lumbar spine BMD and additional 15 CpG-SNPs showed suggestive significant (p < 5.00 × 10-5) association, of which 2 novel SNPs rs7231498 (NFATC1) and rs7455028 (ESR1) also reached a genome-wide significance level in the joint analysis. Several identified CpG-SNPs were mapped to genes that have not been reported for association with BMD in previous GWAS, such as NEK3 and NFATC1 genes, highlighting the enhanced power of targeted association analysis for identification of novel associations that were missed by traditional GWAS. Interestingly, several genomic regions, such as NEK3 and LRP5 regions, contained multiple significant/suggestive CpG-SNPs for lumbar spine BMD, suggesting that multiple neighboring CpG-SNPs may synergistically mediate the DNA methylation level and gene expression pattern of target genes. Furthermore, functional annotation analyses suggested a strong regulatory potential of the identified BMD-associated CpG-SNPs and a significant enrichment in biological processes associated with protein localization and protein signal transduction. Our results provided novel insights into the genetic basis of BMD variation and highlighted the close connections between genetic and epigenetic mechanisms of complex disease.

Glucocorticoids suppress Wnt16 expression in osteoblasts in vitro and in vivo

Glucocorticoid-induced osteoporosis is a frequent complication of systemic glucocorticoid (GC) therapy and mainly characterized by suppressed osteoblast activity. Wnt16 derived from osteogenic cells is a key determinant of bone mass. Here, we assessed whether GC suppress bone formation via inhibiting Wnt16 expression. GC treatment with dexamethasone (DEX) decreased Wnt16 mRNA levels in murine bone marrow stromal cells (mBMSCs) time- and dose-dependently. Similarly, Wnt16 expression was also suppressed after DEX treatment in calvarial organ cultures. Consistently, mice receiving GC-containing slow-release prednisolone pellets showed lower skeletal Wnt16 mRNA levels and bone mineral density than placebo-treated mice. The suppression of Wnt16 by GCs was GC-receptor-dependent as co-treatment of mBMSCs with DEX and the GR antagonist RU-486 abrogated the GC-mediated suppression of Wnt16. Likewise, DEX failed to suppress Wnt16 expression in GR knockout-mBMSCs. In addition, Wnt16 mRNA levels were unaltered in bone tissue of GC-treated GR dimerization-defective GR^dim mice, suggesting that GCs suppress Wnt16 via direct DNA-binding mechanisms. Consistently, DEX treatment reduced Wnt16 promoter activity in MC3T3-E1 cells. Finally, recombinant Wnt16 restored DEX-induced suppression of bone formation in mouse calvaria. Thus, this study identifies Wnt16 as a novel target of GC action in GC-induced suppression of bone formation.

Inducible Wnt16 inactivation: WNT16 regulates cortical bone thickness in adult mice

Substantial progress has been made in the therapeutic reduction of vertebral fracture risk in patients with osteoporosis, but non-vertebral fracture risk has been improved only marginally. Human genetic studies demonstrate that the WNT16 locus is a major determinant of cortical bone thickness and non-vertebral fracture risk and mouse models with life-long Wnt16 inactivation revealed that WNT16 is a key regulator of cortical thickness. These studies, however, could not exclude that the effect of Wnt16 inactivation on cortical thickness might be caused by early developmental and/or growth effects. To determine the effect of WNT16 specifically on adult cortical bone homeostasis, Wnt16 was conditionally ablated in young adult and old mice through tamoxifen-inducible Cre-mediated recombination using CAG-Cre-ER; Wnt16^flox/flox (Cre-Wnt16^flox/flox) mice. First, 10-week-old Cre-Wnt16^flox/flox and Wnt16^flox/flox littermate control mice were treated with tamoxifen. Four weeks later, Wnt16 mRNA levels in cortical bone were reduced and cortical thickness in femur was decreased in Cre-Wnt16^flox/flox mice compared to Wnt16^flox/flox mice. Then, inactivation of Wnt16 in 47-week-old mice (evaluated four weeks later) resulted in a reduction of Wnt16 mRNA levels, cortical thickness and cortical bone strength with no effect on trabecular bone volume fraction. Mechanistic studies demonstrated that the reduced cortical bone thickness was caused by a combination of increased bone resorption and reduced periosteal bone formation. In conclusion, WNT16 is a crucial regulator of cortical bone thickness in young adult and old mice. We propose that new treatment strategies targeting the adult regulation of WNT16 might be useful to reduce fracture risk at cortical bone sites.

Interaction between LRP5 and periostin gene polymorphisms on serum periostin levels and cortical bone microstructure

We investigated the interaction between periostin SNPs and the SNPs of the genes assumed to modulate serum periostin levels and bone microstructure in a cohort of postmenopausal women. We identified an interaction between LRP5 SNP rs648438 and periostin SNP rs9547970 on serum periostin levels and on radial cortical porosity.

PURPOSE

The purpose of this study is to investigate the interaction between periostin gene polymorphisms (SNPs) and other genes potentially responsible for modulating serum periostin levels and bone microstructure in a cohort of postmenopausal women.

METHODS

In 648 postmenopausal women from the Geneva Retirees Cohort, we analyzed 6 periostin SNPs and another 149 SNPs in 14 genes, namely BMP2, CTNNB1, ESR1, ESR2, LRP5, LRP6, PTH, SPTBN1, SOST, TGFb1, TNFRSF11A, TNFSF11, TNFRSF11B and WNT16. Volumetric BMD and bone microstructure were measured by high-resolution peripheral quantitative computed tomography at the distal radius and tibia.

RESULTS

Serum periostin levels were associated with radial cortical porosity, including after adjustment for age, BMI, and years since menopause (p = 0.036). Sixteen SNPs in the ESR1, LRP5, TNFRSF11A, SOST, SPTBN1, TNFRSF11B and TNFSF11 genes were associated with serum periostin levels (p range 0.03-0.001) whereas 26 SNPs in 9 genes were associated with cortical porosity at the radius and/or at the tibia. WNT 16 was the gene with the highest number of SNPs associated with both trabecular and cortical microstructure. The periostin SNP rs9547970 was also associated with cortical porosity (p = 0.04). In particular, SNPs in LRP5, ESR1 and near the TNFRSF11A gene were associated with both cortical porosity and serum periostin levels. Eventually, we identified an interaction between LRP5 SNP rs648438 and periostin SNP rs9547970 on serum periostin levels (interaction p = 0.01) and on radial cortical porosity (interaction p = 0.005).

CONCLUSION

These results suggest that periostin expression is genetically modulated, particularly by polymorphisms in the Wnt pathway, and is thereby implicated in the genetic variation of bone microstructure.

Life-Course Genome-wide Association Study Meta-analysis of Total Body BMD and Assessment of Age-Specific Effects

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.

The RSPO3 gene as genetic markers for bone mass assessed by quantitative ultrasound in a population of young adults

Ultrasound bone mass measurement has been postulated as a valuable bone-health assessment tool for primary care. The aim of this study was to analyse the possible relationship between the SPTBN1, RSPO3, CCDC170, DKK1, GPATCH1, and TMEM135 genes, with calcaneal quantitative ultrasound (QUS) in a population of young adults. These genes were first associated with broadband ultrasound attenuation (BUA) in the GEFOS/GENOMOS study. A cross-sectional study was conducted on 575 individuals (mean age 20.41 ± 2.69). Bone mass at the right calcaneus was estimated by QUS. Six single-nucleotide polymorphisms (SNPs) in SPTBN1 (rs11898505), RSPO3 (rs7741021), CCDC170 (rs4869739), DKK1 (rs7902708), TMEM135 (rs597319), and GPATCH1 (rs10416265) were selected as genetic markers based on their previous association with calcaneal QUS. After adjusting for multiple confounding factors, the only significant association with QUS in our population was found for the rs7741021 SNP in the RSPO3 gene (P = 0.006) using the dominant model of inheritance. This suggests the possible implication of the RSPO3 gene in bone mass acquisition during early adulthood.

Genome-wide association study of coronary artery calcified atherosclerotic plaque in African Americans with type 2 diabetes

Background

Coronary artery calcified atherosclerotic plaque (CAC) predicts cardiovascular disease (CVD). Despite exposure to more severe conventional CVD risk factors, African Americans (AAs) are less likely to develop CAC, and when they do, have markedly lower levels than European Americans. Genetic factors likely contribute to the observed ethnic differences. To identify genes associated with CAC in AAs with type 2 diabetes (T2D), a genome-wide association study (GWAS) was performed using the Illumina 5 M chip in 691 African American-Diabetes Heart Study participants (AA-DHS), with replication in 205 Jackson Heart Study (JHS) participants with T2D. Genetic association tests were performed on the genotyped and 1000 Genomes-imputed markers separately for each study, and combined in a meta-analysis.

Results

Single nucleotide polymorphisms (SNPs), rs11353135 (2q22.1), rs16879003 (6p22.3), rs5014012, rs58071836 and rs10244825 (all on chromosome 7), rs10918777 (9q31.2), rs13331874 (16p13.3) and rs4459623 (18q12.1) were associated with presence and/or quantity of CAC in the AA-DHS and JHS, with meta-analysis p-values ≤8.0 × 10⁻⁷. The strongest result in AA-DHS alone was rs6491315 in the 13q32.1 region (parameter estimate (SE) = −1.14 (0.20); p-value = 9.1 × 10⁻⁹). This GWAS peak replicated a previously reported AA-DHS CAC admixture signal (rs7492028, LOD score 2.8).

Conclusions

Genetic association between SNPs on chromosomes 2, 6, 7, 9, 16 and 18 and CAC were detected in AAs with T2D from AA-DHS and replicated in the JHS. These data support a role for genetic variation on these chromosomes as contributors to CAC in AAs with T2D, as well as to variation in CAC between populations of African and European ancestry.

Electronic supplementary material

The online version of this article (10.1186/s12863-017-0572-9) contains supplementary material, which is available to authorized users.

Analysis of SNP-SNP interactions and bone quantitative ultrasound parameter in early adulthood

Background

Osteoporosis individual susceptibility is determined by the interaction of multiple genetic variants and environmental factors. The aim of this study was to conduct SNP-SNP interaction analyses in candidate genes influencing heel quantitative ultrasound (QUS) parameter in early adulthood to identify novel insights into the mechanism of disease.

Methods

The study population included 575 healthy subjects (mean age 20.41; SD 2.36). To assess bone mass QUS was performed to determine Broadband ultrasound attenuation (BUA, dB/MHz). A total of 32 SNPs mapping to loci that have been characterized as genetic markers for QUS and/or BMD parameters were selected as genetic markers in this study. The association of all possible SNP pairs with QUS was assessed by linear regression and a SNP-SNP interaction was defined as a significant departure from additive effects.

Results

The pairwise SNP-SNP analysis showed multiple interactions. The interaction comprising SNPs rs9340799 and rs3736228 that map in the ESR1 and LRP5 genes respectively, revealed the lowest p value after adjusting for confounding factors (p-value = 0.001, β (95% CI) = 14.289 (5.548, 23.029). In addition, our model reported others such as TMEM135-WNT16 (p = 0.007, β(95%CI) = 9.101 (2.498, 15.704), ESR1-DKK1 (p = 0.012, β(95%CI) = 13.641 (2.959, 24.322) or OPG-LRP5 (p = 0.012, β(95%CI) = 8.724 (1.936, 15.512). However, none of the detected interactions remain significant considering the Bonferroni significance threshold for multiple testing (p<0.0001).

Conclusion

Our analysis of SNP-SNP interaction in candidate genes of QUS in Caucasian young adults reveal several interactions, especially between ESR1 and LRP5 genes, that did not reach statistical significance. Although our results do not support a relevant genetic contribution of SNP-SNP epistatic interactions to QUS in young adults, further studies in larger independent populations would be necessary to support these preliminary findings.

Non-synonymous WNT16 polymorphisms alleles are associated with different osteoarthritis phenotypes

Hereditary factors have a strong influence on osteoarthritis (OA). The Wnt pathway is involved in bone and cartilage homeostasis. Hence, we hypothesized that allelic variations of WNT16 could influence the OA phenotype. We studied 509 Caucasian patients undergoing joint replacement due to severe primary OA. Radiographs were used to classify the OA as atrophic or hypertrophic. Two nonsynonymous polymorphisms of WNT16 (rs2707466 and rs2908004) were analyzed. The association between the genotypes and the OA phenotype was analyzed by logistic regression and adjusted for age and body mass index. A genotype-phenotype association was found in the sex-stratified analysis. Thus, there was a significant difference in the genotypic frequencies of rs2707466 between hypertrophic and atrophic hip OA in males (p = 0.003), with overrepresentation of G alleles in the hypertrophic phenotype (OR 2.08; CI 1.28-3.38). An association in the same direction was observed between these alleles and the type of knee OA, with G alleles being more common in the hypertrophic than in atrophic knee phenotypes (p = 0.008; OR 1.956, CI 1.19-3.19). Similar associations were found for the rs2908004 SNP, but it only reached statistical significance for knee OA (p = 0.017; OR 0.92, CI 0.86-0.989). This is the first study attempting to explore the association of genetic variants with the OA phenotype. These data suggest the need to consider the OA phenotype in future genetic association studies of OA.

The Protein Encoded by the CCDC170 Breast Cancer Gene Functions to Organize the Golgi-Microtubule Network

Genome-Wide Association Studies (GWAS) and subsequent fine-mapping studies (>50) have implicated single nucleotide polymorphisms (SNPs) located at the CCDC170/C6ORF97-ESR1 locus (6q25.1) as being associated with the risk of breast cancer. Surprisingly, our analysis using genome-wide differential allele-specific expression (DASE), an indicator for breast cancer susceptibility, suggested that the genetic alterations of CCDC170, but not ESR1, account for GWAS-associated breast cancer risk at this locus. Breast cancer-associated CCDC170 nonsense mutations and rearrangements have also been detected, with the latter being specifically implicated in driving breast cancer. Here we report that the wild type CCDC170 protein localizes to the region of the Golgi apparatus and binds Golgi-associated microtubules (MTs), and that breast cancer-linked truncations of CCDC170 result in loss of Golgi localization. Overexpression of wild type CCDC170 triggers Golgi reorganization, and enhances Golgi-associated MT stabilization and acetyltransferase ATAT1-dependent α-tubulin acetylation. Golgi-derived MTs regulate cellular polarity and motility, and we provide evidence that dysregulation of CCDC170 affects polarized cell migration. Taken together, our findings demonstrate that CCDC170 plays an essential role in Golgi-associated MT organization and stabilization, and implicate a mechanism for how perturbations in the CCDC170 gene may contribute to the hallmark changes in cell polarity and motility seen in breast cancer.

Bivariate genome-wide association meta-analysis of pediatric musculoskeletal traits reveals pleiotropic effects at the SREBF1/TOM1L2 locus

Bone mineral density is known to be a heritable, polygenic trait whereas genetic variants contributing to lean mass variation remain largely unknown. We estimated the shared SNP heritability and performed a bivariate GWAS meta-analysis of total-body lean mass (TB-LM) and total-body less head bone mineral density (TBLH-BMD) regions in 10,414 children. The estimated SNP heritability is 43% (95% CI: 34–52%) for TBLH-BMD, and 39% (95% CI: 30–48%) for TB-LM, with a shared genetic component of 43% (95% CI: 29–56%). We identify variants with pleiotropic effects in eight loci, including seven established bone mineral density loci: WNT4, GALNT3, MEPE, CPED1/WNT16, TNFSF11, RIN3, and PPP6R3/LRP5. Variants in the TOM1L2/SREBF1 locus exert opposing effects TB-LM and TBLH-BMD, and have a stronger association with the former trait. We show that SREBF1 is expressed in murine and human osteoblasts, as well as in human muscle tissue. This is the first bivariate GWAS meta-analysis to demonstrate genetic factors with pleiotropic effects on bone mineral density and lean mass.

Bone mineral density and lean skeletal mass are heritable traits. Here, Medina-Gomez and colleagues perform bivariate GWAS analyses of total body lean mass and bone mass density in children, and show genetic loci with pleiotropic effects on both traits.

Bone Mass and Strength are Significantly Improved in Mice Overexpressing Human WNT16 in Osteocytes

Recently, we demonstrated that osteoblast-specific overexpression of human WNT16 increased both cortical and trabecular bone mass and structure in mice. To further identify the cell-specific role of Wnt16 in bone homeostasis, we created transgenic (TG) mice overexpressing human WNT16 in osteocytes using Dmp1 promoter (Dmp1-hWNT16 TG) on C57BL/6 (B6) background. We analyzed bone phenotypes and serum bone biomarkers, performed gene expression analysis and measured dynamic bone histomorphometry in Dmp1-hWNT16 TG and wild-type (WT) mice. Compared to WT mice, Dmp1-hWNT16 TG mice exhibited significantly higher whole-body, spine and femoral aBMD, BMC and trabecular (BV/TV, Tb.N, and Tb.Th) and cortical (bone area and thickness) parameters in both male and female at 12 weeks of age. Femur stiffness and ultimate force were also significantly improved in the Dmp1-hWNT16 TG female mice, compared to sex-matched WT littermates. In addition, female Dmp1-hWNT16 TG mice displayed significantly higher MS/BS, MAR and BFR/BS compared to the WT mice. Gene expression analysis demonstrated significantly higher mRNA level of Alp in both male and female Dmp1-hWNT16 TG mice and significantly higher levels of Osteocalcin, Opg and Rankl in the male Dmp1-hWNT16 TG mice in bone tissue compared to sex-matched WT mice. These results indicate that WNT16 plays a critical role for acquisition of both cortical and trabecular bone mass and strength. Strategies designed to use WNT16 as a target for therapeutic interventions will be valuable to treat osteoporosis and other low bone mass conditions.

Association Between Single Gene Polymorphisms and Bone Biomarkers and Response to Calcium and Vitamin D Supplementation in Young Adults Undergoing Military Training

Initial military training (IMT) is associated with increased stress fracture risk. In prior studies, supplemental calcium (Ca) and vitamin D provided daily throughout IMT reduced stress fracture incidence, suppressed parathyroid hormone (PTH), and improved measures of bone health compared with placebo. Data were analyzed from a randomized, double-blind, placebo-controlled trial to determine whether single-nucleotide polymorphisms (SNPs) in Ca and vitamin D-related genes were associated with circulating biomarkers of bone metabolism in young adults entering IMT, and whether responses to Ca and vitamin D supplementation were modulated by genotype. Associations between SNPs, including vitamin D receptor (VDR), vitamin D binding protein (DBP), and 1-alpha-hydroxylase (CYP27B1), and circulating biomarkers were measured in fasting blood samples from volunteers (n = 748) starting IMT. Volunteers were block randomized by race and sex to receive Ca (2000 mg) and vitamin D (1000 IU) or placebo daily throughout Army or Air Force IMT (7 to 9 weeks). Total Ca and vitamin D intakes were calculated as the sum of supplemental intake based on intervention compliance and dietary intake. Relationships between SNPs, Ca, and vitamin D intake tertile and change in biomarkers were evaluated in trial completers (n = 391). At baseline, the minor allele of a DBP SNP (rs7041) was positively associated with both 25OHD (B = 4.46, p = 1.97E-10) and 1,25(OH)₂ D₃ (B = 9.63, p < 0.001). Combined genetic risk score (GRS) for this SNP and a second SNP in the VDR gene (rs1544410) was inversely associated with baseline 25OHD (r = -0.28, p < 0.001) and response to Ca and vitamin D intake differed by GRS (p < 0.05). In addition, presence of the minor allele of a second VDR SNP (rs2228570) was associated with lower P1NP (B = -4.83, p = 0.04) and osteocalcin (B = -0.59, p = 0.03). These data suggest that VDR and DBP SNPs are associated with 25OHD status and bone turnover and those with the highest GRS require the greatest vitamin D intake to improve 25OHD during IMT. © 2016 American Society for Bone and Mineral Research.