Meta-analysis of genome-wide scans for total body BMD in children and adults reveals allelic heterogeneity and age-specific effects at the WNT16 locus

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.

Amoli M, Fana SE, Dimai HP, Obermayer-Pietsch B, Luegger B, Rivadeneira F, Nabipour I, Larijani B, Khashayar P. Protocol for preliminary, multicenteric validation of "PoCOsteo device": A point of care tool for proteomic and genomic study of osteoporosis

One of the goals of the HORIZON 2020 project PoCOsteo was to develop a medical device, which would measure and/or quantify proteomic as well as genomic factors as present in whole blood samples collected through finger prick. After validating the tool in the clinical setting, the next step would be its clinical validation based on the existing guidelines. This article presents the protocol of a validation study to be carried out independently at two different centers (Division of Endocrinology and Diabetology at the Medical University of Graz as a clinic-based cohort, and the Endocrinology and Metabolism Research Institute at the Tehran University of Medical Sciences as a population-based cohort). It aims to assess the tool according to the Clinical & Laboratory Standards Institute guidelines, confirming if the proteomics and genomics measurements provided by the tool are accurate and reproducible compared with the existing state-of-the-art tests. This is the first time that such a detailed protocol for lab validation of a medical tool for proteomics and genomic measurement is designed based on the existing guidelines and thus could be used as a template for clinical validation of future point-of-care tools. Moreover, the multicentric cohort design will allow the study of a large number of diverse individuals, which will improve the validity and generalizability of the results for different settings.

A longitudinal genome-wide association study of bone mineral density mean and variability in the UK Biobank

Bone mineral density (BMD) is an essential predictor of osteoporosis and fracture. We conducted a genome-wide trajectory analysis of BMD and analyzed the BMD change.

PURPOSE

This study aimed to identify the genetic architecture and potential biomarkers of BMD.

METHODS

Our analysis included 141,261 white participants from the UK Biobank with heel BMD phenotype data. We used a genome-wide trajectory analysis tool, TrajGWAS, to conduct a genome-wide association study (GWAS) of BMD. Then, we validated our findings in previously reported BMD genetic associations and performed replication analysis in the Asian participants. Finally, gene-set enrichment analysis (GSEA) of the identified candidate genes was conducted using the FUMA platform.

RESULTS

A total of 52 genes associated with BMD trajectory mean were identified, of which the top three significant genes were WNT16 (P = 1.31 × 10-126), FAM3C (P = 4.18 × 10-108), and CPED1 (P = 8.48 × 10-106). In addition, 114 genes associated with BMD within-subject variability were also identified, such as AC092079.1 (P = 2.72 × 10-13) and RGS7 (P = 4.72 × 10-10). The associations for these candidate genes were confirmed in the previous GWASs and replicated successfully in the Asian participants. GSEA results of BMD change identified multiple GO terms related to skeletal development, such as SKELETAL SYSTEM DEVELOPMENT (Padjusted = 2.45 × 10-3) and REGULATION OF OSSIFICATION (Padjusted = 2.45 × 10-3). KEGG enrichment analysis showed that these genes were mainly enriched in WNT SIGNALING PATHWAY.

CONCLUSIONS

Our findings indicated that the CPED1-WNT16-FAM3C locus plays a significant role in BMD mean trajectories and identified several novel candidate genes contributing to BMD within-subject variability, facilitating the understanding of the genetic architecture of BMD.

A study of the influence of genetic variance and sex on the density and thickness of the calvarial bone in collaborative cross mice

BACKGROUND

Bone microarchitecture is affected by multiple genes, each having a small effect on the external appearance. It is thus challenging to characterize the genes and their specific effect on bone thickness and porosity. The purpose of this study was to assess the heritability and the genetic variation effect, as well as the sex effect on the calvarial bone thickness (Ca.Th) and calvarial porosity (%PoV) using the Collaborative Cross (CC) mouse population.

METHODS

In the study we examined the parietal bones of 56 mice from 9 lines of CC mice. Morphometric parameters were evaluated using microcomputed tomography (μCT) and included Ca.Th and %PoV. We then evaluated heritability, genetic versus environmental variance and the sex effect for these parameters.

RESULTS

Our morphometric analysis showed that Ca.Th and %PoV are both significantly different among the CC lines with a broad sense heritability of 0.78 and 0.90, respectively. The sex effect within the lines was significant in line IL111 and showed higher values of Ca.Th and %PoV in females compared to males. In line IL19 there was a borderline sex effect in Ca.Th in which males showed higher values than females.

CONCLUSIONS

These results stress the complexity of sex and genotype interactions controlling Ca.Th and %PoV, as the skeletal sexual dimorphism was dependent on the genetic background. This study also shows that the CC population is a powerful tool for establishing the genetic effect on these traits.

Examining craniofacial variation among crispant and mutant zebrafish models of human skeletal diseases

Genetic diseases affecting the skeletal system present with a wide range of symptoms that make diagnosis and treatment difficult. Genome-wide association and sequencing studies have identified genes linked to human skeletal diseases. Gene editing of zebrafish models allows researchers to further examine the link between genotype and phenotype, with the long-term goal of improving diagnosis and treatment. While current automated tools enable rapid and in-depth phenotyping of the axial skeleton, characterizing the effects of mutations on the craniofacial skeleton has been more challenging. The objective of this study was to evaluate a semi-automated screening tool can be used to quantify craniofacial variations in zebrafish models using four genes that have been associated with human skeletal diseases (meox1, plod2, sost, and wnt16) as test cases. We used traditional landmarks to ground truth our dataset and pseudolandmarks to quantify variation across the 3D cranial skeleton between the groups (somatic crispant, germline mutant, and control fish). The proposed pipeline identified variation between the crispant or mutant fish and control fish for four genes. Variation in phenotypes parallel human craniofacial symptoms for two of the four genes tested. This study demonstrates the potential as well as the limitations of our pipeline as a screening tool to examine multi-dimensional phenotypes associated with the zebrafish craniofacial skeleton.

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.

Associations Between Body Mass Index, WNT16 rs2908004 and Osteoporosis: Findings from Taiwan Biobank

Purpose

Osteoporosis is a degenerative disease that affects women and men of all races. We studied the association between body mass index (BMI), rs2908004 polymorphism of the WNT16 gene, and osteoporosis using data from Taiwan Biobank (TWB).

Patients and Methods

We analyzed data from 10,942 subjects aged 30 to 70. We defined osteoporosis based on a mean T-score of -2.5 and below in the hip. Body mass index was classified following the guidelines of the Health Promotion Administration. Imputation was carried out using the IMPUTE2 (v2.3.1) program. Multiple logistic regression was used for analysis. The odds ratios (ORs) and 95% confidence interval (CI) for osteoporosis were determined.

Results

In the multivariate regression model, variant rs2908004 had a significant association with osteoporosis. That is, the rs2908004-GA+AA genotype was associated with lower osteoporosis risk than the GG genotype (OR, 0.651; 95% CI = 0.544 to 0.780). Compared to normal-weight, underweight was significantly associated with a higher risk of osteoporosis (OR, 6.517; 95% CI = 4.624 to 9.186) while overweight and obesity were protective (OR, 0.176; 95% CI = 0.140 to 0.221 and 0.057; 95% CI = 0.039 to 0.083, respectively). There was an interaction between rs2908004 and BMI (p = 0.0148). Subgroup analyses (using rs2908004-GG/normal-weight as the reference group) indicated ORs of 7.66 (95% CI = 5.153 to 11.394) in the rs2908004-GG/underweight group and 3.002 (95% CI = 1.509 to 5.974) in the rs2908004-GA+AA/underweight group (95% CI = 1.509 to 5.974). Odds ratios were substantially lower in rs2908004-GG/obese, rs2908004-GG/overweight, GA+AA/normal-weight, rs2908004-GA+AA/overweight, and rs2908004-GA+AA/obese groups, respectively.

Conclusion

According to our study, underweight was associated with an increased risk of osteoporosis irrespective of WNT16 rs2908004 genotypes, while overweight and obesity were associated with a lower risk.

wnt16 regulates spine and muscle morphogenesis through parallel signals from notochord and dermomyotome

Bone and muscle are coupled through developmental, mechanical, paracrine, and autocrine signals. Genetic variants at the CPED1-WNT16 locus are dually associated with bone- and muscle-related traits. While Wnt16 is necessary for bone mass and strength, this fails to explain pleiotropy at this locus. Here, we show wnt16 is required for spine and muscle morphogenesis in zebrafish. In embryos, wnt16 is expressed in dermomyotome and developing notochord, and contributes to larval myotome morphology and notochord elongation. Later, wnt16 is expressed at the ventral midline of the notochord sheath, and contributes to spine mineralization and osteoblast recruitment. Morphological changes in wnt16 mutant larvae are mirrored in adults, indicating that wnt16 impacts bone and muscle morphology throughout the lifespan. Finally, we show that wnt16 is a gene of major effect on lean mass at the CPED1-WNT16 locus. Our findings indicate that Wnt16 is secreted in structures adjacent to developing bone (notochord) and muscle (dermomyotome) where it affects the morphogenesis of each tissue, thereby rendering wnt16 expression into dual effects on bone and muscle morphology. This work expands our understanding of wnt16 in musculoskeletal development and supports the potential for variants to act through WNT16 to influence bone and muscle via parallel morphogenetic processes.

In humans, genetic variants (DNA sequences that vary amongst individuals) have been identified that appear to influence two tissues, bone and skeletal muscle. However, how single genes and genetic variants exert dual influence on both tissues is not well understood. In this study, we found that the wnt16 gene is necessary for specifying the size and shape of both muscle and bone during development in zebrafish. We also disentangled how wnt16 affects both tissues: distinct cellular populations adjacent to muscle and bone secrete Wnt16, where it acts as a signal guiding the size and shape of each tissue. This is important because in humans, genetic variants near the WNT16 gene have effects on both bone- and muscle-related traits. This study expands our understanding of the role of WNT16 in bone and muscle development, and helps to explain how genetic variants near WNT16 affect traits for both tissues. Moreover, WNT16 is actively being explored as a target for osteoporosis therapies, thus our study could have implications with regard to the potential of targeting WNT16 to treat bone and muscle simultaneously.

Where is bone science taking us?

Bone science has over the last decades unraveled many important pathways in bone and mineral metabolism and the interplay between genetic factors and the environment. Some of these discoveries have led to the development of pharmacological treatments of osteoporosis and rare bone diseases. Other scientific avenues have uncovered a role for the gut microbiome in regulating bone mass, which have led to investigations on the possible therapeutic role of probiotics in the prevention of osteoporosis. Huge advances have been made in identifying the genes that cause rare bone diseases, which in some cases have led to therapeutic interventions. Advances have also been made in understanding the genetic basis of the more common polygenic bone diseases, including osteoporosis and Paget's disease of bone (PDB). Polygenic profiles are used for establishing genetic risk scores aiming at early diagnosis and intervention, but also in Mendelian randomization (MR) studies to investigate both desired and undesired effects of targets for drug design.

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.

Early-Onset Osteoporosis

Osteoporosis is a skeletal disorder with enhanced bone fragility, usually affecting the elderly. It is very rare in children and young adults and the definition is not only based on a low BMD (a Z-score < - 2.0 in growing children and a Z-score ≤ - 2.0 or a T-score ≤ - 2.5 in young adults) but also on the occurrence of fragility fractures and/or the existence of underlying chronic diseases or secondary factors such as use of glucocorticoids. In the absence of a known chronic disease, fragility fractures and low BMD should prompt extensive screening for secondary causes, which can be found in up to 90% of cases. When fragility fractures occur in childhood or young adulthood without an evident secondary cause, investigations should explore the possibility of an underlying monogenetic bone disease, where bone fragility is caused by a single variant in a gene that has a major role in the skeleton. Several monogenic forms relate to type I collagen, but other forms also exist. Loss-of-function variants in LRP5 and WNT1 may lead to early-onset osteoporosis. The X-chromosomal osteoporosis caused by PLS3 gene mutations affects especially males. Another recently discovered form relates to disturbed sphingolipid metabolism due to SGMS2 mutations, underscoring the complexity of molecular pathology in monogenic early-onset osteoporosis. Management of young patients consists of treatment of secondary factors, optimizing lifestyle factors including calcium and vitamin D and physical exercise. Treatment with bone-active medication should be discussed on a personalized basis, considering the severity of osteoporosis and underlying disease versus the absence of evidence on anti-fracture efficacy and potential harmful effects in pregnancy.

Competitive blocking of LRP4-sclerostin binding interface strongly promotes bone anabolic functions

Induction of bone formation by Wnt ligands is inhibited when sclerostin (Scl), an osteocyte-produced antagonist, binds to its receptors, the low-density lipoprotein receptor-related proteins 5 or 6 (LRP5/6). Recently, it was shown that enhanced inhibition is achieved by Scl binding to the co-receptor LRP4. However, it is not clear if the binding of Scl to LRP4 facilitates Scl binding to LRP5/6 or inhibits the Wnt pathway in an LRP5/6-independent manner. Here, using the yeast display system, we demonstrate that Scl exhibits a stronger binding affinity for LRP4 than for LRP6. Moreover, we found stronger Scl binding to LRP6 in the presence of LRP4. We further show that a Scl mutant (SclN93A), which tightly binds LRP4 but not LRP6, does not inhibit the Wnt pathway on its own. We demonstrate that SclN93A competes with Scl for a common binding site on LRP4 and antagonizes Scl inhibition of the Wnt signaling pathway in osteoblasts in vitro. Finally, we demonstrate that 2 weeks of bi-weekly subcutaneous injections of SclN93A fused to the fragment crystallizable (Fc) domain of immunoglobulin (SclN93AFc), which retains the antagonistic activity of the mutant, significantly increases bone formation rate and enhances trabecular volumetric bone fraction, trabecular number, and bone length in developing mice. Our data show that LRP4 serves as an anchor that facilitates Scl-LRP6 binding and that inhibition of the Wnt pathway by Scl depends on its prior binding to LRP4. We further provide evidence that compounds that inhibit Scl-LRP4 interactions offer a potential strategy to promote anabolic bone functions.

Functional Validation of Osteoporosis Genetic Findings Using Small Fish Models

The advancement of human genomics has revolutionized our understanding of the genetic architecture of many skeletal diseases, including osteoporosis. However, interpreting results from human association studies remains a challenge, since index variants often reside in non-coding regions of the genome and do not possess an obvious regulatory function. To bridge the gap between genetic association and causality, a systematic functional investigation is necessary, such as the one offered by animal models. These models enable us to identify causal mechanisms, clarify the underlying biology, and apply interventions. Over the past several decades, small teleost fishes, mostly zebrafish and medaka, have emerged as powerful systems for modeling the genetics of human diseases. Due to their amenability to genetic intervention and the highly conserved genetic and physiological features, fish have become indispensable for skeletal genomic studies. The goal of this review is to summarize the evidence supporting the utility of Zebrafish (Danio rerio) for accelerating our understanding of human skeletal genomics and outlining the remaining gaps in knowledge. We provide an overview of zebrafish skeletal morphophysiology and gene homology, shedding light on the advantages of human skeletal genomic exploration and validation. Knowledge of the biology underlying osteoporosis through animal models will lead to the translation into new, better and more effective therapeutic approaches.

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.

The genetic overlap between osteoporosis and craniosynostosis

Osteoporosis is the most prevalent bone condition in the ageing population. This systemic disease is characterized by microarchitectural deterioration of bone, leading to increased fracture risk. In the past 15 years, genome-wide association studies (GWAS), have pinpointed hundreds of loci associated with bone mineral density (BMD), helping elucidate the underlying molecular mechanisms and genetic architecture of fracture risk. However, the challenge remains in pinpointing causative genes driving GWAS signals as a pivotal step to drawing the translational therapeutic roadmap. Recently, a skull BMD-GWAS uncovered an intriguing intersection with craniosynostosis, a congenital anomaly due to premature suture fusion in the skull. Here, we recapitulate the genetic contribution to both osteoporosis and craniosynostosis, describing the biological underpinnings of this overlap and using zebrafish models to leverage the functional investigation of genes associated with skull development and systemic skeletal homeostasis.

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.

R-spondin 3 deletion induces Erk phosphorylation to enhance Wnt signaling and promote bone formation in the appendicular skeleton

Activation of Wnt signaling leads to high bone density. The R-spondin family of four secreted glycoproteins (Rspo1-4) amplifies Wnt signaling. In humans, RSPO3 variants are strongly associated with bone density. Here, we investigated the role of Rspo3 in skeletal homeostasis in mice. Using a comprehensive set of mouse genetic and mechanistic studies, we show that in the appendicular skeleton, Rspo3 haplo-insufficiency and Rspo3 targeted deletion in Runx2⁺ osteoprogenitors lead to an increase in trabecular bone mass, with increased number of osteoblasts and bone formation. In contrast and highlighting the complexity of Wnt signaling in the regulation of skeletal homeostasis, we show that Rspo3 deletion in osteoprogenitors results in the opposite phenotype in the axial skeleton, i.e., low vertebral trabecular bone mass. Mechanistically, Rspo3 deficiency impairs the inhibitory effect of Dkk1 on Wnt signaling activation and bone mass. We demonstrate that Rspo3 deficiency leads to activation of Erk signaling which in turn, stabilizes β-catenin and Wnt signaling activation. Our data demonstrate that Rspo3 haplo-insufficiency/deficiency boosts canonical Wnt signaling by activating Erk signaling, to favor osteoblastogenesis, bone formation, and bone mass.

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.

The Polygenic and Monogenic Basis of Paediatric Fractures

PURPOSE OF REVIEW

Fractures are frequently encountered in paediatric practice. Although recurrent fractures in children usually unveil a monogenic syndrome, paediatric fracture risk could be shaped by the individual genetic background influencing the acquisition of bone mineral density, and therefore, the skeletal fragility as shown in adults. Here, we examine paediatric fractures from the perspective of monogenic and complex trait genetics.

RECENT FINDINGS

Large-scale genome-wide studies in children have identified ~44 genetic loci associated with fracture or bone traits whereas ~35 monogenic diseases characterized by paediatric fractures have been described. Genetic variation can predispose to paediatric fractures through monogenic risk variants with a large effect and polygenic risk involving many variants of small effects. Studying genetic factors influencing peak bone attainment might help in identifying individuals at higher risk of developing early-onset osteoporosis and discovering drug targets to be used as bone restorative pharmacotherapies to prevent, or even reverse, bone loss later in life.

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.

Loss of Wnt16 Leads to Skeletal Deformities and Downregulation of Bone Developmental Pathway in Zebrafish

Wingless-type MMTV integration site family, member 16 (wnt16), is a wnt ligand that participates in the regulation of vertebrate skeletal development. Studies have shown that wnt16 can regulate bone metabolism, but its molecular mechanism remains largely undefined. We obtained the wnt16^−/− zebrafish model using the CRISPR-Cas9-mediated gene knockout screen with 11 bp deletion in wnt16, which led to the premature termination of amino acid translation and significantly reduced wnt16 expression, thus obtaining the wnt16^−/− zebrafish model. The expression of wnt16 in bone-related parts was detected via in situ hybridization. The head, spine, and tail exhibited significant deformities, and the bone mineral density and trabecular bone decreased in wnt16^−/− using light microscopy and micro-CT analysis. RNA sequencing was performed to explore the differentially expressed genes (DEGs). Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis found that the down-regulated DEGs are mainly concentrated in mTOR, FoxO, and VEGF pathways. Protein–protein interaction (PPI) network analysis was performed with the detected DEGs. Eight down-regulated DEGs including akt1, bnip4, ptena, vegfaa, twsg1b, prkab1a, prkab1b, and pla2g4f.2 were validated by qRT-PCR and the results were consistent with the RNA-seq data. Overall, our work provides key insights into the influence of wnt16 gene on skeletal development.

Approaches to osteoporosis in paleopathology: How did methodology shape bone loss research?

OBJECTIVE

This paper will review how different methods employed to study bone loss in the past were used to explore different questions and aspects of bone loss, how methodology has changed over time, and how these different approaches have informed our understanding of bone loss in the past.

MATERIALS AND METHODS

A review and discussion is conducted on research protocols and results of 84 paleopathology publications on bone loss in archaeological skeletal collections published between 1969 and 2021.

CONCLUSIONS

The variety in research protocols confounds accurate meta-analysis of previously published research; however, more recent publications incorporate a combination of bone mass and bone quality based methods. Biased sample selection has resulted in a predominance of European and Medieval publications, limiting more general observations on bone loss in the past. Collection of dietary or paleopathological covariables is underemployed in the effort to interpret bone loss patterns.

SIGNIFICANCE

Paleopathology publications have demonstrated differences in bone loss between distinct archaeological populations, between sex and age groups, and have suggested factors underlying observed differences. However, a lack of a gold standard has encouraged the use of a wide range of methods. Understanding how this array of methods effects results is crucial in contextualizing our knowledge of bone loss in the past.

LIMITATIONS

The development of a research protocol is also influenced by available expertise, available equipment, restrictions imposed by the curator, and site-specific taphonomic aspects. These factors will likely continue to cause (minor) biases even if a best practice can be established.

SUGGESTIONS FOR FUTURE RESEARCH

Greater effort to develop uniform terminology and operational definitions of osteoporosis in skeletal remains, as well as the expansion of time scale and geographical areas studied. The Next-Generation Sequencing revolution has also opened up the possibility of ancient DNA analyses to study genetic predisposition to bone loss in the past.

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.

Genetic variants in the FAM3C gene are associated with lipid traits in Chinese children

BACKGROUND

Previous studies have related FAM3C gene with childhood bone health, and the regulation of lipid metabolism in hepatocytes. The present case-control study aimed to analyze the association of FAM3C genetic variants with overweight/obesity and lipid traits among Chinese children.

METHODS

Two genetic variants (rs7776725 and rs7793554) within the FAM3C gene were genotyped in 3305 Chinese children aged 6-18 years.

RESULTS

In the whole study population, the T-allele of rs7776725 and A-allele of rs7793554 within the FAM3C gene were associated with 40.2% (95% CI: 11.6-76.1%; P = 0.004) and 29.1% (6.9-56.0%; P = 0.008) increased risk of dyslipidemia, higher triglyceride (P = 0.014 and P = 0.001) and lower HDL-C (P = 0.015 and P = 0.003). In addition, we found that rs7776725 interacted with sex on dyslipidemia (P_{for interaction} = 0.004), and sex-stratified analyses showed that it was significantly associated with dyslipidemia only in girls (P = 8.78 × 10^-5). The variant also showed nominally significant interactions with sex on total cholesterol and LDL-C (P_{for interaction} = 0.012 and 0.008).

CONCLUSION

We found that FAM3C genetic variants were associated with dyslipidemia and lipid traits among Chinese children. In addition, we found significant gene-by-sex interactions. Our findings provided evidence supporting the role of FAM3C gene in regulating lipid metabolism in humans.

IMPACT

FAM3C genetic variants were associated with dyslipidemia and lipid traits among Chinese children. In addition, we found significant gene-by-sex interactions. FAM3C/rs7776725 was associated with dyslipidemia and lipid traits only in girls. Our findings provided evidence supporting the role of FAM3C gene in regulating lipid metabolism in humans.

Genome-Wide Association Analysis of Longitudinal Bone Mineral Content Data From the Iowa Bone Development Study

The foundation for osteoporosis risk is, in part, established during childhood, adolescence, and young adulthood, all periods of development when bone mass is acquired rapidly. The relative quantity of bone mass accrued is influenced by both lifestyle and genetic factors, although the genetic component is not yet well understood. The purpose of this study was to use a genome-wide association (GWA) analysis to discover single nucleotide polymorphisms (SNPs) associated with: (1) the sex-specific hip bone mineral content at approximately the age of 19 when the amount of bone accrued is near its peak; and (2) the sex-specific rate of hip bone mineral content accrual during the adolescent growth spurt. The Iowa Bone Development Study, a longitudinal cohort study exploring bone health in children, adolescents, and young adults was the source of data. From this cohort, n = 364 (190 females, 174 males) participants were included in GWA analyses to address (1) and n = 258 participants (125 females and 133 males) were included in GWA analyses to address (2). Twenty SNPS were detected having p < 1.0 × 10-5. Of most biologic relevance were 2 suggestive SNPs (rs2051756 and rs2866908) detected in an intron of the DKK2 gene through the GWA analysis that explored peak bone mass in females.

USF3 modulates osteoporosis risk by targeting WNT16, RANKL, RUNX2, and two GWAS lead SNPs rs2908007 and rs4531631

Osteoporotic fractures cause major morbidity and mortality in the aging population. Genome-wide association studies (GWAS) have identified USF3 as the novel susceptibility gene of osteoporosis. However, the functional role in bone metabolism and the target gene of the basic helix-loop-helix transcription factor USF3 are unclear. Here, we show that USF3 enhances osteoblast differentiation and suppresses osteoclastogenesis in cultured human osteoblast-like U-2OS cells. Mechanistic studies revealed that transcription factor USF3 antagonistically interacts with anti-osteogenic TWIST1/TCF12 heterodimer in the WNT16 and RUNX2 promoter, and counteracts CREB1 and JUN/FOS in the RANKL promoter. Importantly, the osteoporosis GWAS variant g.1744A>G (rs2908007A>G) located in the WNT16 promoter confers G-allele-specific transcriptional modulation by USF3, TWIST1/TCF12 and TBX5/TBX15, and USF3 transactivates the osteoclastogenesis suppressor WNT16 promoter activity and antagonizes the repression of WNT16 by TWIST1 and TCF12. The risk G allele of osteoporosis GWAS variant g.3260A>G (rs4531631A>G) in the RANKL promoter facilitates the binding of CREB1 and JUN/FOS and enhances transactivation of the osteoclastogenesis contributor RANKL that is inhibited by USF3. Our findings uncovered the functional mechanisms of osteoporosis novel GWAS-associated gene USF3 and lead single nucleotide polymorphisms rs2908007 and rs4531631 in the regulation of bone formation and resorption.

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.

Whole exome sequencing of a Saudi family and systems biology analysis identifies CPED1 as a putative causative gene to Celiac Disease

Celiac disease (CD) is a gastrointestinal disorder whose genetic basis is not fully understood. Therefore, we studied a Saudi family with two CD affected siblings to discover the causal genetic defect. Through whole exome sequencing (WES), we identified that both siblings have inherited an extremely rare and deleterious CPED1 genetic variant (c.241 A > G; p.Thr81Ala) segregating as autosomal recessive mutation, suggesting its putative causal role in the CD. Saudi population specific minor allele frequency (MAF) analysis has confirmed its extremely rare prevalence in homozygous condition (MAF is 0.0004). The Sanger sequencing analysis confirmed the absence of this homozygous variant in 100 sporadic Saudi CD cases. Genotype-Tissue Expression (GTEx) data has revealed that CPED1 is abundantly expressed in gastrointestinal mucosa. By using a combination of systems biology approaches like protein 3D modeling, stability analysis and nucleotide sequence conservation analysis, we have further established that this variant is deleterious to the structural and functional aspects of CPED1 protein. To the best of our knowledge, this variant has not been previously reported in CD or any other gastrointestinal disease. The cell culture and animal model studies could provide further insight into the exact role of CPED1 p.Thr81Ala variant in the pathophysiology of CD. In conclusion, by using WES and systems biology analysis, present study for the first-time reports CPED1 as a potential causative gene for CD in a Saudi family with potential implications to both disease diagnosis and genetic counseling.

Genetic variants affecting bone mineral density and bone mineral content at multiple skeletal sites in Hispanic children

CONTEXT

Osteoporosis is a major public health burden with significant economic costs. However, the correlates of bone health in Hispanic children are understudied.

OBJECTIVE

We aimed to identify genetic variants associated with bone mineral density (BMD) and bone mineral content (BMC) at multiple skeletal sites in Hispanic children.

METHODS

We conducted a cross-sectional genome-wide linkage analysis, genome-wide and exome-wide association analysis of BMD and BMC. The Viva La Familia Study is a family-based cohort with a total of 1030 Hispanic children (4-19 years old at baseline) conducted in Houston, TX. BMD and BMC were measured by Dual-energy X-ray absorptiometry.

RESULTS

Significant heritability were observed for BMC and BMD at multiple skeletal sites ranging between 44 and 68% (P < 2.8 × 10-9). Significant evidence for linkage was found for BMD of pelvis and left leg on chromosome 7p14, lumbar spine on 20q13 and left rib on 6p21, and BMC of pelvis on chromosome 20q12 and total body on 14q22-23 (logarithm of odds score > 3). We found genome-wide significant association between BMC of right arm and rs762920 at PVALB (P = 4.6 × 10-8), and between pelvis BMD and rs7000615 at PTK2B (P = 7.4 × 10-8). Exome-wide association analysis revealed novel association of variants at MEGF10 and ABRAXAS2 with left arm and lumber spine BMC, respectively (P < 9 × 10-7).

CONCLUSIONS

We identified novel loci associated with BMC and BMD in Hispanic children, with strongest evidence for PTK2B. These findings provide better understanding of bone genetics and shed light on biological mechanisms underlying BMD and BMC variation.

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.

Inducible expression of Wnt7b promotes bone formation in aged mice and enhances fracture healing

There remain unmet clinical needs for safe and effective bone anabolic therapies to treat aging-related osteoporosis and to improve fracture healing in cases of nonunion or delayed union. Wnt signaling has emerged as a promising target pathway for developing novel bone anabolic drugs. Although neutralizing antibodies against the Wnt antagonist sclerostin have been tested, Wnt ligands themselves have not been fully explored as a potential therapy. Previous work has demonstrated Wnt7b as an endogenous ligand upregulated during osteoblast differentiation, and that Wnt7b overexpression potently stimulates bone accrual in the mouse. The earlier studies however did not address whether Wnt7b could promote bone formation when specifically applied to aged or fractured bones. Here we have developed a doxycycline-inducible strategy where Wnt7b is temporally induced in the bones of aged mice or during fracture healing. We report that forced expression of Wnt7b for 1 month starting at 15 months of age greatly stimulated trabecular and endosteal bone formation, resulting in a marked increase in bone mass. We further tested the effect of Wnt7b on bone healing in a murine closed femur fracture model. Induced expression of Wnt7b at the onset of fracture did not affect the initial cartilage formation but promoted mineralization of the subsequent bone callus. Thus, targeted delivery of Wnt7b to aged bones or fracture sites may be explored as a potential therapy.

Coordination of WNT signaling and ciliogenesis during odontogenesis by piezo type mechanosensitive ion channel component 1

Signal transmission from the mechanical forces to the various intracellular activities is a fundamental process during tissue development. Despite their critical role, the mechanism of mechanical forces in the biological process is poorly understood. In this study, we demonstrated that in the response to hydrostatic pressure (HP), the piezo type mechanosensitive ion channel component 1 (PIEZO1) is a primary mechanosensing receptor for odontoblast differentiation through coordination of the WNT expression and ciliogenesis. In stem cells from human exfoliated deciduous teeth (SHED), HP significantly promoted calcium deposition as well as the expression of odontogenic marker genes, PANX3 and DSPP, and WNT related-genes including WNT5b and WNT16, whereas HP inhibited cell proliferation and enhanced primary cilia expression. WNT signaling inhibitor XAV939 and primary cilia inhibitor chloral hydrate blocked the HP-induced calcium deposition. The PIEZO1 activator Yoda1 inhibited cell proliferation but induced ciliogenesis and WNT16 expression. Interestingly, HP and Yoda1 promoted nuclear translocation of RUNX2, whereas siRNA-mediated silencing of PIEZO1 decreased HP-induced nuclear translocation of RUNX2. Taken together, these results suggest that PIEZO1 functions as a mechanotransducer that connects HP signal to the intracellular signalings during odontoblast differentiation.

Using zebrafish to study skeletal genomics

While genome-wide association studies (GWAS) have revolutionized our understanding of the genetic architecture of skeletal diseases, animal models are required to identify causal mechanisms and to translate underlying biology into new therapies. Despite large-scale knockout mouse phenotyping efforts, the skeletal functions of most genes residing at GWAS-identified loci remain unknown, highlighting a need for complementary model systems to accelerate gene discovery. Over the past several decades, zebrafish (Danio rerio) has emerged as a powerful system for modeling the genetics of human diseases. In this review, our goal is to outline evidence supporting the utility of zebrafish for accelerating our understanding of human skeletal genomics, as well as gaps in knowledge that need to be filled for this purpose. We do this by providing a basic foundation of the zebrafish skeletal morphophysiology and phenotypes, and surveying evidence of skeletal gene homology and the use of zebrafish for post-GWAS analysis in other tissues and organs. We also outline challenges in translating zebrafish mutant phenotypes. Finally, we conclude with recommendations of future directions and how to leverage the large body of tools and knowledge of skeletal genetics in zebrafish for the needs of human skeletal genomic exploration. Due to their amenability to rapid genetic approaches, as well as the large number of conserved genetic and phenotypic features, there is a strong rationale supporting the use of zebrafish for human skeletal genomic studies.

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.

Genotype-by-environment interactions inferred from genetic effects on phenotypic variability in the UK Biobank

Genotype-by-environment interaction (GEI) is a fundamental component in understanding complex trait variation. However, it remains challenging to identify genetic variants with GEI effects in humans largely because of the small effect sizes and the difficulty of monitoring environmental fluctuations. Here, we demonstrate that GEI can be inferred from genetic variants associated with phenotypic variability in a large sample without the need of measuring environmental factors. We performed a genome-wide variance quantitative trait locus (vQTL) analysis of ~5.6 million variants on 348,501 unrelated individuals of European ancestry for 13 quantitative traits in the UK Biobank and identified 75 significant vQTLs with P < 2.0 × 10-9 for 9 traits, especially for those related to obesity. Direct GEI analysis with five environmental factors showed that the vQTLs were strongly enriched with GEI effects. Our results indicate pervasive GEI effects for obesity-related traits and demonstrate the detection of GEI without environmental data.

Predicting human disease mutations and identifying drug targets from mouse gene knockout phenotyping campaigns

Two large-scale mouse gene knockout phenotyping campaigns have provided extensive data on the functions of thousands of mammalian genes. The ongoing International Mouse Phenotyping Consortium (IMPC), with the goal of examining all ∼20,000 mouse genes, has examined 5115 genes since 2011, and phenotypic data from several analyses are available on the IMPC website (www.mousephenotype.org). Mutant mice having at least one human genetic disease-associated phenotype are available for 185 IMPC genes. Lexicon Pharmaceuticals' Genome5000™ campaign performed similar analyses between 2000 and the end of 2008 focusing on the druggable genome, including enzymes, receptors, transporters, channels and secreted proteins. Mutants (4654 genes, with 3762 viable adult homozygous lines) with therapeutically interesting phenotypes were studied extensively. Importantly, phenotypes for 29 Lexicon mouse gene knockouts were published prior to observations of similar phenotypes resulting from homologous mutations in human genetic disorders. Knockout mouse phenotypes for an additional 30 genes mimicked previously published human genetic disorders. Several of these models have helped develop effective treatments for human diseases. For example, studying Tph1 knockout mice (lacking peripheral serotonin) aided the development of telotristat ethyl, an approved treatment for carcinoid syndrome. Sglt1 (also known as Slc5a1) and Sglt2 (also known as Slc5a2) knockout mice were employed to develop sotagliflozin, a dual SGLT1/SGLT2 inhibitor having success in clinical trials for diabetes. Clinical trials evaluating inhibitors of AAK1 (neuropathic pain) and SGLT1 (diabetes) are underway. The research community can take advantage of these unbiased analyses of gene function in mice, including the minimally studied 'ignorome' genes.

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.

Increased glycolysis mediates Wnt7b-induced bone formation

Wingless/integrated (Wnt) signaling has emerged as a major mechanism for promoting bone formation and a target pathway for developing bone anabolic agents against osteoporosis. However, the downstream events mediating the potential therapeutic effect of Wnt proteins are not fully understood. Previous studies have indicated that increased glycolysis is associated with osteoblast differentiation in response to Wnt signaling, but direct genetic evidence for the importance of glucose metabolism in Wnt-induced bone formation is lacking. Here, we have generated compound transgenic mice to overexpress Wnt family member 7B (Wnt7b) transiently in the osteoblast lineage of postnatal mice, with or without concurrent deletion of the glucose transporter 1 (Glut1), also known as solute carrier family 2, facilitated glucose transporter member 1. Overexpression of Wnt7b in 1-mo-old mice for 1 wk markedly stimulated bone formation, but the effect was essentially abolished without Glut1, even though transient deletion of Glut1 itself did not affect normal bone accrual. Consistent with the in vivo results, Wnt7b increased Glut1 expression and glucose consumption in the primary culture of osteoblast lineage cells, and deletion of Glut1 diminished osteoblast differentiation in vitro. Thus, Wnt7b promotes bone formation in part through stimulating glucose metabolism in osteoblast lineage cells.-Chen, H., Ji, X., Lee, W.-C., Shi, Y., Li, B., Abel, E. D., Jiang, D., Huang, W., Long, F. Increased glycolysis mediates Wnt7b-induced bone formation.

The interactions of physical activity, exercise and genetics and their associations with bone mineral density: implications for injury risk in elite athletes

Low bone mineral density (BMD) is established as a primary predictor of osteoporotic risk and can also have substantial implications for athlete health and injury risk in the elite sporting environment. BMD is a highly multi-factorial phenotype influenced by diet, hormonal characteristics and physical activity. The interrelationships between such factors, and a strong genetic component, suggested to be around 50-85% at various anatomical sites, determine skeletal health throughout life. Genome-wide association studies and case-control designs have revealed many loci associated with variation in BMD. However, a number of the candidate genes identified at these loci have no known associated biological function or have yet to be replicated in subsequent investigations. Furthermore, few investigations have considered gene-environment interactions-in particular, whether specific genes may be sensitive to mechanical loading from physical activity and the outcome of such an interaction for BMD and potential injury risk. Therefore, this review considers the importance of physical activity on BMD, genetic associations with BMD and how subsequent investigation requires consideration of the interaction between these determinants. Future research using well-defined independent cohorts such as elite athletes, who experience much greater mechanical stress than most, to study such phenotypes, can provide a greater understanding of these factors as well as the biological underpinnings of such a physiologically "extreme" population. Subsequently, modification of training, exercise or rehabilitation programmes based on genetic characteristics could have substantial implications in both the sporting and public health domains once the fundamental research has been conducted successfully.

Mechanoactivation of Wnt/β-catenin pathways in health and disease

Mechanical forces play an important role in regulating tissue development and homeostasis in multiple cell types including bone, joint, epithelial and vascular cells, and are also implicated in the development of diseases, e.g. osteoporosis, cardiovascular disease and osteoarthritis. Defining the mechanisms by which cells sense and respond to mechanical forces therefore has important implications for our understanding of tissue function in health and disease and may lead to the identification of targets for therapeutic intervention. Mechanoactivation of the Wnt signalling pathway was first identified in osteoblasts with a key role for β-catenin demonstrated in loading-induced osteogenesis. Since then, mechanoregulation of the Wnt pathway has also been observed in stem cells, epithelium, chondrocytes and vascular and lymphatic endothelium. Wnt can signal through both canonical and non-canonical pathways, and evidence suggests that both can mediate responses to mechanical strain, stretch and shear stress. This review will discuss our current understanding of the activation of the Wnt pathway in response to mechanical forces.

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.

Genome-wide association studies and CRISPR/Cas9-mediated gene editing identify regulatory variants influencing eyebrow thickness in humans

Hair plays an important role in primates and is clearly subject to adaptive selection. While humans have lost most facial hair, eyebrows are a notable exception. Eyebrow thickness is heritable and widely believed to be subject to sexual selection. Nevertheless, few genomic studies have explored its genetic basis. Here, we performed a genome-wide scan for eyebrow thickness in 2961 Han Chinese. We identified two new loci of genome-wide significance, at 3q26.33 near SOX2 (rs1345417: P = 6.51×10(-10)) and at 5q13.2 near FOXD1 (rs12651896: P = 1.73×10(-8)). We further replicated our findings in the Uyghurs, a population from China characterized by East Asian-European admixture (N = 721), the CANDELA cohort from five Latin American countries (N = 2301), and the Rotterdam Study cohort of Dutch Europeans (N = 4411). A meta-analysis combining the full GWAS results from the three cohorts of full or partial Asian descent (Han Chinese, Uyghur and Latin Americans, N = 5983) highlighted a third signal of genome-wide significance at 2q12.3 (rs1866188: P = 5.81×10(-11)) near EDAR. We performed fine-mapping and prioritized four variants for further experimental verification. CRISPR/Cas9-mediated gene editing provided evidence that rs1345417 and rs12651896 affect the transcriptional activity of the nearby SOX2 and FOXD1 genes, which are both involved in hair development. Finally, suitable statistical analyses revealed that none of the associated variants showed clear signals of selection in any of the populations tested. Contrary to popular speculation, we found no evidence that eyebrow thickness is subject to strong selective pressure.

Assessment of the genetic and clinical determinants of fracture risk: genome wide association and mendelian randomisation study

OBJECTIVES

To identify the genetic determinants of fracture risk and assess the role of 15 clinical risk factors on osteoporotic fracture risk.

DESIGN

Meta-analysis of genome wide association studies (GWAS) and a two-sample mendelian randomisation approach.

SETTING

25 cohorts from Europe, United States, east Asia, and Australia with genome wide genotyping and fracture data.

PARTICIPANTS

A discovery set of 37 857 fracture cases and 227 116 controls; with replication in up to 147 200 fracture cases and 150 085 controls. Fracture cases were defined as individuals (>18 years old) who had fractures at any skeletal site confirmed by medical, radiological, or questionnaire reports. Instrumental variable analyses were performed to estimate effects of 15 selected clinical risk factors for fracture in a two-sample mendelian randomisation framework, using the largest previously published GWAS meta-analysis of each risk factor.

RESULTS

Of 15 fracture associated loci identified, all were also associated with bone mineral density and mapped to genes clustering in pathways known to be critical to bone biology (eg, SOST, WNT16, and ESR1) or novel pathways (FAM210A, GRB10, and ETS2). Mendelian randomisation analyses showed a clear effect of bone mineral density on fracture risk. One standard deviation decrease in genetically determined bone mineral density of the femoral neck was associated with a 55% increase in fracture risk (odds ratio 1.55 (95% confidence interval 1.48 to 1.63; P=1.5×10-68). Hand grip strength was inversely associated with fracture risk, but this result was not significant after multiple testing correction. The remaining clinical risk factors (including vitamin D levels) showed no evidence for an effect on fracture.

CONCLUSIONS

This large scale GWAS meta-analysis for fracture identified 15 genetic determinants of fracture, all of which also influenced bone mineral density. Among the clinical risk factors for fracture assessed, only bone mineral density showed a major causal effect on fracture. Genetic predisposition to lower levels of vitamin D and estimated calcium intake from dairy sources were not associated with fracture risk.

Common and rare variants of WNT16, DKK1 and SOST and their relationship with bone mineral density

Numerous GWAS and candidate gene studies have highlighted the role of the Wnt pathway in bone biology. Our objective has been to study in detail the allelic architecture of three Wnt pathway genes: WNT16, DKK1 and SOST, in the context of osteoporosis. We have resequenced the coding and some regulatory regions of these three genes in two groups with extreme bone mineral density (BMD) (n = ∼50, each) from the BARCOS cohort. No interesting novel variants were identified. Thirteen predicted functional variants have been genotyped in the full cohort (n = 1490), and for ten of them (with MAF > 0.01), the association with BMD has been studied. We have found six variants nominally associated with BMD, of which 2 WNT16 variants predicted to be eQTLs for FAM3C (rs55710688, in the Kozak sequence and rs142005327, within a putative enhancer) withstood multiple-testing correction. In addition, two rare variants in functional regions (rs190011371 in WNT16b 3′UTR and rs570754792 in the SOST TATA box) were found only present in three women each, all with BMD below the mean of the cohort. Our results reinforce the higher importance of regulatory versus coding variants in these Wnt pathway genes and open new ways for functional studies of the relevant variants.

Global gene expression analysis identifies Mef2c as a potential player in Wnt16-mediated transcriptional regulation

Wnt16 is a major Wnt ligand involved in the regulation of postnatal bone homeostasis. Previous studies have shown that Wnt16 promotes bone formation and inhibits bone resorption, suggesting that this molecule could be targeted for therapeutic interventions to treat bone thinning disorders such as osteoporosis. However, the molecular mechanisms by which Wnt16 regulates bone metabolism is not yet fully understood. To better understand the molecular mechanisms by which Wnt16 promotes bone formation and to identify the target genes regulated by Wnt16 in osteoblasts, we treated calvarial osteoblasts purified from C57Bl/6 mice with recombinant Wnt16 and profiled the gene expression changes by RNA-seq at 24 h post-treatment. We also compared gene expression profiles of Wnt16-treated osteoblasts to canonical Wnt3a- and non-canonical Wnt5a-treated osteoblasts. This study identified 576 genes differentially expressed in Wnt16-treated osteoblasts compared to sham-treated controls; these included several members of Wnt pathway (Wnt2b, Wnt7b, Wnt11, Axin2, Sfrp2, Sfrp4, Fzd5 etc.) and TGF-β/BMP signaling pathway (Bmp7, Inhba, Inhbb, Tgfb2 etc.). Wnt16 also regulated a large number of genes with known bone phenotypes. We also found that about 37% (215/576) of the Wnt16 targets overlapped with Wnt3a targets and ~15% (86/576) overlapped with Wnt5a targets, suggesting that Wnt16 activates both canonical and non-canonical Wnt signaling targets in osteoblasts. Transcription factor binding motif enrichment analysis in the promoter regions of Wnt16 targets identified noncanonical Wnt/JNK pathway activated transcription factors Fosl2 and Fosl1 as two of the most significantly enriched transcription factors associated with genes activated by Wnt16 while Mef2c was the most significantly enriched transcription factor associated with genes repressed by Wnt16. We also found that a large number of Mef2c targets overlapped with genes down-regulated by Wnt16 and Mef2c itself was transcriptionally repressed by Wnt16 suggesting that Mef2c plays a role in Wnt16-mediated transcriptional regulation.

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.