Phenotypic dissection of bone mineral density reveals skeletal site specificity and facilitates the identification of novel loci in the genetic regulation of bone mass attainment

Genetic Evaluation for Monogenic Disorders of Low Bone Mass and Increased Bone Fragility: What Clinicians Need to Know

PURPOSE OF REVIEW

The purpose of this review is to outline the principles of clinical genetic testing and to provide practical guidance to clinicians in navigating genetic testing for patients with suspected monogenic forms of osteoporosis.

RECENT FINDINGS

Heritability assessments and genome-wide association studies have clearly shown the significant contributions of genetic variations to the pathogenesis of osteoporosis. Currently, over 50 monogenic disorders that present primarily with low bone mass and increased risk of fractures have been described. The widespread availability of clinical genetic testing offers a valuable opportunity to correctly diagnose individuals with monogenic forms of osteoporosis, thus instituting appropriate surveillance and treatment. Clinical genetic testing may identify the appropriate diagnosis in a subset of patients with low bone mass, multiple or unusual fractures, and severe or early-onset osteoporosis, and thus clinicians should be aware of how to incorporate such testing into their clinical practices.

Bone health index in the assessment of bone health: The Generation R Study

Bone Health Index (BHI) has been proposed as a useful instrument for assessing bone health in children. However, its relationship with fracture risk remains unknown. We aimed to investigate whether BHI is associated with bone mineral density (BMD) and prevalent fracture odds in children from the Generation R Study. We also implemented genome-wide association study (GWAS) and polygenic score (PGS) approaches to improve our understanding of BHI and its potential. In total, 4150 children (49.4 % boys; aged 9.8 years) with genotyped data and bone assessments were included in this study. BMD was measured across the total body (less head following ISCD guidelines) using a GE-Lunar iDXA densitometer; and BHI was determined from the hand DXA scans using BoneXpert®. Fractures were self-reported collected with home questionnaires. The association of BHI with BMD and fractures was evaluated using linear models corrected for age, sex, ethnicity, height, and weight. We observed a positive correlation between BHI and BMD (ρ = 0.32, p-value<0.0001). Further, every SD decrease in BHI was associated with an 11 % increased risk of prevalent fractures (OR:1.11, 95 % CI 1.00-1.24, p-value = 0.05). Our BHI GWAS identified variants (lead SNP rs1404264-A, p-value = 2.61 × 10-14) mapping to the ING3/CPED1/WNT16 locus. Children in the extreme tails of the BMD PGS presented a difference in BHI values of -0.10 standard deviations (95% CI -0.14 to -0.07; p-value<0.0001). On top of the demonstrated epidemiological association of BHI with both BMD and fracture risk, our results reveal a partially shared biological background between BHI and BMD. These findings highlight the potential value of using BHI to screen children at risk of fracture.

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.

Revealing chronic disease progression patterns using Gaussian process for stage inference

OBJECTIVE

The early stages of chronic disease typically progress slowly, so symptoms are usually only noticed until the disease is advanced. Slow progression and heterogeneous manifestations make it challenging to model the transition from normal to disease status. As patient conditions are only observed at discrete timestamps with varying intervals, an incomplete understanding of disease progression and heterogeneity affects clinical practice and drug development.

MATERIALS AND METHODS

We developed the Gaussian Process for Stage Inference (GPSI) approach to uncover chronic disease progression patterns and assess the dynamic contribution of clinical features. We tested the ability of the GPSI to reliably stratify synthetic and real-world data for osteoarthritis (OA) in the Osteoarthritis Initiative (OAI), bipolar disorder (BP) in the Adolescent Brain Cognitive Development Study (ABCD), and hepatocellular carcinoma (HCC) in the UTHealth and The Cancer Genome Atlas (TCGA).

RESULTS

First, GPSI identified two subgroups of OA based on image features, where these subgroups corresponded to different genotypes, indicating the bone-remodeling and overweight-related pathways. Second, GPSI differentiated BP into two distinct developmental patterns and defined the contribution of specific brain region atrophy from early to advanced disease stages, demonstrating the ability of the GPSI to identify diagnostic subgroups. Third, HCC progression patterns were well reproduced in the two independent UTHealth and TCGA datasets.

CONCLUSION

Our study demonstrated that an unsupervised approach can disentangle temporal and phenotypic heterogeneity and identify population subgroups with common patterns of disease progression. Based on the differences in these features across stages, physicians can better tailor treatment plans and medications to individual patients.

Zebrafish as a Model for Osteoporosis: Functional Validations of Genome-Wide Association Studies

PURPOSE OF REVIEW

GWAS, as a largely correlational analysis, requires in vitro or in vivo validation. Zebrafish (Danio rerio) have many advantages for studying the genetics of human diseases. Since gene editing in zebrafish has been highly valuable for studying embryonic skeletal developmental processes that are prenatally or perinatally lethal in mammalian models, we are reviewing pros and cons of this model.

RECENT FINDINGS

The true power for the use of zebrafish is the ease by which the genome can be edited, especially using the CRISPR/Cas9 system. Gene editing, followed by phenotyping, for complex traits such as BMD, is beneficial, but the major physiological differences between the fish and mammals must be considered. Like mammals, zebrafish do have main bone cells; thus, both in vivo stem cell analyses and in vivo imaging are doable. Yet, the "long" bones of fish are peculiar, and their bone cavities do not contain bone marrow. Partial duplication of the zebrafish genome should be taken into account. Overall, small fish toolkit can provide unmatched opportunities for genetic modifications and morphological investigation as a follow-up to human-first discovery.

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.

Lowering of Circulating Sclerostin May Increase Risk of Atherosclerosis and Its Risk Factors: Evidence From a Genome-Wide Association Meta-Analysis Followed by Mendelian Randomization

OBJECTIVE

In this study, we aimed to establish the causal effects of lowering sclerostin, target of the antiosteoporosis drug romosozumab, on atherosclerosis and its risk factors.

METHODS

A genome-wide association study meta-analysis was performed of circulating sclerostin levels in 33,961 European individuals. Mendelian randomization (MR) was used to predict the causal effects of sclerostin lowering on 15 atherosclerosis-related diseases and risk factors.

RESULTS

We found that 18 conditionally independent variants were associated with circulating sclerostin. Of these, 1 cis signal in SOST and 3 trans signals in B4GALNT3, RIN3, and SERPINA1 regions showed directionally opposite signals for sclerostin levels and estimated bone mineral density. Variants with these 4 regions were selected as genetic instruments. MR using 5 correlated cis-SNPs suggested that lower sclerostin increased the risk of type 2 diabetes mellitus (DM) (odds ratio [OR] 1.32 [95% confidence interval (95% CI) 1.03-1.69]) and myocardial infarction (MI) (OR 1.35 [95% CI 1.01-1.79]); sclerostin lowering was also suggested to increase the extent of coronary artery calcification (CAC) (β = 0.24 [95% CI 0.02-0.45]). MR using both cis and trans instruments suggested that lower sclerostin increased hypertension risk (OR 1.09 [95% CI 1.04-1.15]), but otherwise had attenuated effects.

CONCLUSION

This study provides genetic evidence to suggest that lower levels of sclerostin may increase the risk of hypertension, type 2 DM, MI, and the extent of CAC. Taken together, these findings underscore the requirement for strategies to mitigate potential adverse effects of romosozumab treatment on atherosclerosis and its related risk factors.

Identification of an intronic enhancer regulating RANKL expression in osteocytic cells

The bony skeleton is continuously renewed throughout adult life by the bone remodeling process, in which old or damaged bone is removed by osteoclasts via largely unknown mechanisms. Osteocytes regulate bone remodeling by producing the osteoclast differentiation factor RANKL (encoded by the TNFSF11 gene). However, the precise mechanisms underlying RANKL expression in osteocytes are still elusive. Here, we explored the epigenomic landscape of osteocytic cells and identified a hitherto-undescribed osteocytic cell-specific intronic enhancer in the TNFSF11 gene locus. Bioinformatics analyses showed that transcription factors involved in cell death and senescence act on this intronic enhancer region. Single-cell transcriptomic data analysis demonstrated that cell death signaling increased RANKL expression in osteocytic cells. Genetic deletion of the intronic enhancer led to a high-bone-mass phenotype with decreased levels of RANKL in osteocytic cells and osteoclastogenesis in the adult stage, while RANKL expression was not affected in osteoblasts or lymphocytes. These data suggest that osteocytes may utilize a specialized regulatory element to facilitate osteoclast formation at the bone surface to be resorbed by linking signals from cellular senescence/death and RANKL expression.

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.

SEESAW: detecting isoform-level allelic imbalance accounting for inferential uncertainty

Detecting allelic imbalance at the isoform level requires accounting for inferential uncertainty, caused by multi-mapping of RNA-seq reads. Our proposed method, SEESAW, uses Salmon and Swish to offer analysis at various levels of resolution, including gene, isoform, and aggregating isoforms to groups by transcription start site. The aggregation strategies strengthen the signal for transcripts with high uncertainty. The SEESAW suite of methods is shown to have higher power than other allelic imbalance methods when there is isoform-level allelic imbalance. We also introduce a new test for detecting imbalance that varies across a covariate, such as time.

Bone mineral density loci specific to the skull portray potential pleiotropic effects on craniosynostosis

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

Computed tomography Osteoabsorptiometry: Review of bone density, mechanical strength of material and clinical application

Computed Tomography (CT) imaging is an effective non-invasive examination. It is widely used in the diagnosis of fractures, arthritis, tumor, and some anatomical characteristics of patients. The density value (Hounsfield unit, HU) of a material in computed tomography can be the same for materials with varying elemental compositions. This value depends on the mass density of the material and the degree of X-ray attenuation. Computed Tomography Osteoabsorptiometry (CTOAM) imaging technology is developed on the basis of CT imaging technology. By applying pseudo-color image processing to the articular surface, it is used to analyze the distribution of bone mineralization under the articular cartilage, evaluate the position of prosthesis implantation, track the progression of osteoarthritis, and determine the joint injury prognosis. Furthermore, this technique was combined with indentation testing to discuss the relationship between the high bone density area of the articular surface, the mechanical strength of the bone, and the anchorage stability of the implant, in addition to the study of the relationship between mechanical strength and bone density. This narrative study discusses the pre- and postoperative evaluation of medical device implantation position, orthopedic surgery, and the clinical treatment of bone injury and degeneration. It also discusses the research status of CTOAM technology in image post-processing engineering and the relationship between bone material and mechanical strength.

Genetic variation in WNT16 and its association with bone mineral density, fractures and osteoporosis in children with bone fragility

Several genome-wide association studies (GWAS), GWAS meta-analyses, and mouse studies have demonstrated that wingless-related integration site 16 (WNT16) gene is associated with bone mineral density (BMD), cortical bone thickness, bone strength and fracture risk. Practically no data exist regarding the significance of WNT16 in childhood-onset osteoporosis and related fractures. We hypothesized that pathogenic variants and genetic variations in WNT16 could explain skeletal fragility in affected children. We screened the WNT16 gene by Sanger sequencing in three pediatric cohorts: 35 with primary osteoporosis, 59 with multiple fractures, and in 95 healthy controls. Altogether, we identified 12 variants in WNT16. Of them one was a rare 5′UTR variant rs1386898215 in genome aggregate and medical trans-omic databases (GnomAD, TOPMED; minor allele frequency (MAF) 0.00 and 0.000008, respectively). One variant rs1554366753, overrepresented in children with osteoporosis (MAF = 0.06 vs healthy controls MAF = 0.01), was significantly associated with lower BMD. This variant was found associated with increased WNT16 gene expression at mRNA level in fibroblast cultures. None of the other identified variants were rare (MAF < 0.001) or deemed pathogenic by predictor programs. WNT16 may play a role in childhood osteoporosis but genetic WNT16 variation is not a common cause of skeletal fragility in childhood.

•

No pathogenic WNT16 variants were found associated with pediatric osteoporosis or fracture-prone patients

•

Altogether, twelve WNT16 variants were found in pediatric osteoporosis or fracture-prone patients

•

The genetic variation rs1554366753 in the WNT16 gene is associated with bone mineral density and primary osteoporosis

Ras and Rab Interactor 3: From Cellular Mechanisms to Human Diseases

Ras and Rab interactor 3 (RIN3) functions as a Guanine nucleotide Exchange Factor (GEF) for some members of the Rab family of small GTPase. By promoting the activation of Rab5, RIN3 plays an important role in regulating endocytosis and endocytic trafficking. In addition, RIN3 activates Ras, another small GTPase, that controls multiple signaling pathways to regulate cellular function. Increasing evidence suggests that dysregulation of RIN3 activity may contribute to the pathogenesis of several disease conditions ranging from Paget’s Disease of the Bone (PDB), Alzheimer’s Disease (AD), Chronic Obstructive Pulmonary Disease (COPD) and to obesity. Recent genome-wide association studies (GWAS) identified variants in the RIN3 gene to be linked with these disease conditions. Interestingly, some variants appear to be missense mutations in the functional domains of the RIN3 protein while most variants are located in the noncoding regions of the RIN3 gene, potentially altering its gene expression. However, neither the protein structure of RIN3 nor its exact function(s) (except for its GEF activity) has been fully defined. Furthermore, how the polymorphisms/variants contribute to disease pathogenesis remain to be understood. Herein, we examine, and review published studies in an attempt to provide a better understanding of the physiological function of RIN3; More importantly, we construct a framework linking the polymorphisms/variants of RIN3 to altered cell signaling and endocytic traffic, and to potential disease mechanism(s).

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.

Estrogen receptor alpha and NFATc1 bind to a bone mineral density-associated SNP to repress WNT5B in osteoblasts

Genetic factors and estrogen deficiency contribute to the development of osteoporosis. The single-nucleotide polymorphism (SNP) rs2887571 is predicted from genome-wide association studies (GWASs) to associate with osteoporosis but has had an unknown mechanism. Analysis of osteoblasts from 110 different individuals who underwent joint replacement revealed that the genotype of rs2887571 correlates with WNT5B expression. Analysis of our ChIP-sequencing data revealed that SNP rs2887571 overlaps with an estrogen receptor alpha (ERα) binding site. Here we show that 17β-estradiol (E2) suppresses WNT5B expression and further demonstrate the mechanism of ERα binding at the enhancer containing rs2887571 to suppress WNT5B expression differentially in each genotype. ERα interacts with NFATc1, which is predicted to bind directly at rs2887571. CRISPR-Cas9 and ChIP-qPCR experiments confirm differential regulation of WNT5B between each allele. Homozygous GG has a higher binding affinity for ERα than homozygous AA and results in greater suppression of WNT5B expression. Functionally, WNT5B represses alkaline phosphatase expression and activity, decreasing osteoblast differentiation and mineralization. Furthermore, WNT5B increases interleukin-6 expression and suppresses E2-induced expression of alkaline phosphatase during osteoblast differentiation. We show that WNT5B suppresses the differentiation of osteoblasts via receptor tyrosine kinase-like orphan receptor 1/2 (ROR1/2), which activates DVL2/3/RAC1/CDC42/JNK/SIN3A signaling and inhibits β-catenin activity. Together, our data provide mechanistic insight into how ERα and NFATc1 regulate the non-coding SNP rs2887571, as well as the function of WNT5B on osteoblasts, which could provide alternative therapeutic targets for 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.

Identification of PDXDC1 as a novel pleiotropic susceptibility locus shared between lumbar spine bone mineral density and birth weight

An increasing number of epidemiological studies have suggested that birth weight (BW) may be a determinant of bone health later in life, although the underlying genetic mechanism remains unclear. Here, we applied a pleiotropic conditional false discovery rate (cFDR) approach to the genome-wide association study (GWAS) summary statistics for lumbar spine bone mineral density (LS BMD) and BW, aiming to identify novel susceptibility variants shared between these two traits. We detected 5 novel potential pleiotropic loci which are located at or near 7 different genes (NTAN1, PDXDC1, CACNA1G, JAG1, FAT1P1, CCDC170, ESR1), among which PDXDC1 and FAT1P1 have not previously been linked to these phenotypes. To partially validate the findings, we demonstrated that the expression of PDXDC1 was dramatically reduced in ovariectomized (OVX) mice in comparison with sham-operated (SHAM) mice in both the growth plate and trabecula bone. Furthermore, immunohistochemistry assay with serial sections showed that both osteoclasts and osteoblasts express PDXDC1, supporting its potential role in bone metabolism. In conclusion, our study provides insights into some shared genetic mechanisms for BMD and BW as well as a novel potential therapeutic target for the prevention of OP in the early stages of the disease development. KEY MESSAGES : We investigated pleiotropy-informed enrichment between LS BMD and BW. We identified genetic variants related to both LS BMD and BW by utilizing a cFDR approach. PDXDC1 is a novel pleiotropic gene which may be related to both LS BMD and BW. Elevated expression of PDXDC1 is related to higher BMD and lower ratio n-6/n-3 PUFA indicating a bone protective effect of PDXDC1.

Contrasting effects of <i>Ksr2</i>, an obesity gene, on trabecular bone volume and bone marrow adiposity

Pathological obesity and its complications are associated with an increased propensity for bone fractures. Humans with certain genetic polymorphisms at the kinase suppressor of ras2 (KSR2) locus develop severe early-onset obesity and type 2 diabetes. Both conditions are phenocopied in mice with <i>Ksr2</i> deleted, but whether this affects bone health remains unknown. Here we studied the bones of global <i>Ksr2</i> null mice and found that <i>Ksr2</i> negatively regulates femoral, but not vertebral, bone mass in two genetic backgrounds, while the paralogous gene, <i>Ksr1</i>, was dispensable for bone homeostasis. Mechanistically, KSR2 regulates bone formation by influencing adipocyte differentiation at the expense of osteoblasts in the bone marrow. Compared with <i>Ksr2</i>'s known role as a regulator of feeding by its function in the hypothalamus, pair-feeding and osteoblast-specific conditional deletion of <i>Ksr2</i> reveals that <i>Ksr2</i> can regulate bone formation autonomously. Despite the gains in appendicular bone mass observed in the absence of <i>Ksr2</i>, bone strength, as well as fracture healing response, remains compromised in these mice. This study highlights the interrelationship between adiposity and bone health and provides mechanistic insights into how <i>Ksr2</i>, an adiposity and diabetic gene, regulates bone metabolism.

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.

A transcriptome-wide association study to detect novel genes for volumetric bone mineral density

Transcriptome-wide association studies (TWAS) systematically investigate the association of genetically predicted gene expression with disease risk, providing an effective approach to identify novel susceptibility genes. Osteoporosis is the most common metabolic bone disease, associated with reduced bone mineral density (BMD) and increased risk of osteoporotic fractures, whereas genetic factors explain approximately 70% of the variance in phenotypes associated with bone. BMD is commonly assessed using dual-energy X-ray absorptiometry (DXA) to obtain measurements (g/cm2) of areal BMD. However, quantitative computed tomography (QCT) measured 3D volumetric BMD (vBMD) (g/cm3) has important advantages compared with DXA since it can evaluate cortical and trabecular microstructural features of bone quality, which can be used to directly predict fracture risk. Here, we performed the first TWAS for volumetric BMD (vBMD) by integrating genome-wide association studies (GWAS) data from two independent cohorts, namely the Framingham Heart Study (FHS, n = 3298) and the Osteoporotic Fractures in Men (MrOS, n = 4641), with tissue-specific gene expression data from the Genotype-Tissue Expression (GTEx) project. We first used stratified linkage disequilibrium (LD) score regression approach to identify 12 vBMD-relevant tissues, for which vBMD heritability is enriched in tissue-specific genes of the given tissue. Focusing on these tissues, we subsequently leveraged GTEx expression reference panels to predict tissue-specific gene expression levels based on the genotype data from FHS and MrOS. The associations between predicted gene expression levels and vBMD variation were then tested by MultiXcan, an innovative TWAS method that integrates information available across multiple tissues. We identified 70 significant genes associated with vBMD, including some previously identified osteoporosis-related genes such as LYRM2 and NME8, as well as some novel loci such as DNAAF2 and SPAG16. Our findings provide novel insights into the pathophysiological mechanisms of osteoporosis and highlight several novel vBMD-associated genes that warrant further investigation.

A multiethnic whole genome sequencing study to identify novel loci for bone mineral density

At present, there have only been a few DNA sequencing-based studies to explore the genetic determinants of bone mineral density (BMD). We carried out the largest whole genome sequencing analysis to date for femoral neck and spine BMD (n = 4981), with one of the highest average sequencing depths implemented thus far at 22×, in a multiethnic sample (58% Caucasian and 42% African American) from the Louisiana Osteoporosis Study (LOS). The LOS samples were combined with summary statistics from the GEFOS consortium and several independent samples of various ethnicities to perform GWAS meta-analysis (n = 44 506). We identified 31 and 30 genomic risk loci for femoral neck and spine BMD, respectively. The findings substantiate many previously reported susceptibility loci (e.g. WNT16 and ESR1) and reveal several others that are either novel or have not been widely replicated in GWAS for BMD, including two for femoral neck (IGF2 and ZNF423) and one for spine (SIPA1). Although we were not able to uncover ethnicity specific differences in the genetic determinants of BMD, we did identify several loci which demonstrated sex-specific associations, including two for women (PDE4D and PIGN) and three for men (TRAF3IP2, NFIB and LYSMD4). Gene-based rare variant association testing detected MAML2, a regulator of the Notch signaling pathway, which has not previously been suggested, for association with spine BMD. The findings provide novel insights into the pathophysiological mechanisms of osteoporosis.

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.

Genome-wide functional screen of 3'UTR variants uncovers causal variants for human disease and evolution

3' untranslated region (3'UTR) variants are strongly associated with human traits and diseases, yet few have been causally identified. We developed the massively parallel reporter assay for 3'UTRs (MPRAu) to sensitively assay 12,173 3'UTR variants. We applied MPRAu to six human cell lines, focusing on genetic variants associated with genome-wide association studies (GWAS) and human evolutionary adaptation. MPRAu expands our understanding of 3'UTR function, suggesting that simple sequences predominately explain 3'UTR regulatory activity. We adapt MPRAu to uncover diverse molecular mechanisms at base pair resolution, including an adenylate-uridylate (AU)-rich element of LEPR linked to potential metabolic evolutionary adaptations in East Asians. We nominate hundreds of 3'UTR causal variants with genetically fine-mapped phenotype associations. Using endogenous allelic replacements, we characterize one variant that disrupts a miRNA site regulating the viral defense gene TRIM14 and one that alters PILRB abundance, nominating a causal variant underlying transcriptional changes in age-related macular degeneration.

Targeted Inactivation of Rin3 Increases Trabecular Bone Mass by Reducing Bone Resorption and Favouring Bone Formation

Common genetic variants at the RIN3 locus on chromosome 14q32 predispose to Paget's disease of bone (PDB) but the mechanisms by which they do so are unknown. Here, we analysed the skeletal phenotype of female mice with targeted inactivation of the mouse Rin3 gene (Rin3^-/-) as compared with wild-type littermates. The Rin3^-/- mice had higher trabecular bone volume (BV/TV%) compared with wild type. Mean ± standard deviation values at the distal femur at 8 weeks were 9.0 ± 2.5 vs. 7.0 ± 1.5 (p = 0.002) and at 52 weeks were 15.8 ± 9.5 vs. 8.5 ± 4.2 (p = 0.002). No differences were observed in femoral cortical bone parameters with the exception of marrow diameter which was significantly smaller in 52-week-old Rin3^-/- mice compared to wild type: (0.43 mm ± 0.1 vs. 0.57 mm ± 0.2 (p = 0.001). Bone histomorphometry showed a lower osteoclast surface / bone surface (Oc.S/BS%) at 8 weeks in Rin3^-/- mice compared to wild type (24.1 ± 4.7 vs. 29.7 ± 6.6; p = 0.025) but there were no significant differences in markers of bone formation at this time. At 52 weeks, Oc.S/BS did not differ between genotypes but single labelled perimeter (SL.Pm/B.Pm (%)) was significantly higher in Rin3^-/- mice (24.4 ± 6.4 vs. 16.5 ± 3.8, p = 0.003). We conclude that Rin3 negatively regulates trabecular bone mass in mice by inhibiting osteoclastic bone resorption and favouring bone formation. Our observations also suggest that the variants that predispose to PDB in humans probably do so by causing a gain-in-function of RIN3.

WNT5B in Physiology and Disease

WNT5B, a member of the WNT family of proteins that is closely related to WNT5A, is required for cell migration, cell proliferation, or cell differentiation in many cell types. WNT5B signals through the non-canonical β-catenin-independent signaling pathway and often functions as an antagonist of canonical WNT signaling. Although WNT5B has a high amino acid identity with WNT5A and is often assumed to have similar activities, WNT5B often exhibits unique expression patterns and functions. Here, we describe the distinct effects and mechanisms of WNT5B on development, bone, adipose tissue, cardiac tissue, the nervous system, the mammary gland, the lung and hematopoietic cells, compared to WNT5A. We also highlight aberrances in non-canonical WNT5B signaling contributing to diseases such as osteoarthritis, osteoporosis, obesity, type 2 diabetes mellitus, neuropathology, and chronic diseases associated with aging, as well as various cancers.

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 and physical activity interact to affect bone density in Hispanic children

Background

Our aim was to investigate if moderate to vigorous physical activity (MVPA), calcium intake interacts with bone mineral density (BMD)-related single nucleotide polymorphisms (SNPs) to influence BMD in 750 Hispanic children (4-19y) of the cross-sectional Viva La Familia Study.

Methods

Physical activity and dietary intake were measured by accelerometers and multiple-pass 24 h dietary recalls, respectively. Total body and lumbar spine BMD were measured by dual energy X-ray absorptiometry. A polygenic risk score (PRS) was computed based on SNPs identified in published literature. Regression analysis was conducted with PRSs, MVPA and calcium intake with total body and lumbar spine BMD.

Results

We found evidence of statistically significant interaction effects between the PRS and MVPA on total body BMD and lumbar spine BMD (p < 0.05). Higher PRS was associated with a lower total body BMD (β = − 0.040 ± 0.009, p = 1.1 × 10^− 5) and lumbar spine BMD (β = − 0.042 ± 0.013, p = 0.0016) in low MVPA group, as compared to high MVPA group (β = − 0.015 ± 0.006, p = 0.02; β = 0.008 ± 0.01, p = 0.4, respectively).

Discussion

The study indicated that calcium intake does not modify the relationship between genetic variants and BMD, while it implied physical activity interacts with genetic variants to affect BMD in Hispanic children. Due to limited sample size of our study, future research on gene by environment interaction on bone health and functional studies to provide biological insights are needed.

Conclusions

Bone health in Hispanic children with high genetic risk for low BMD is benefitted more by MVPA than children with low genetic risk. Our results may be useful to predict disease risk and tailor dietary and physical activity advice delivery to people, especially children.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12887-021-02537-y.

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.

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 study implicates novel loci and reveals candidate effector genes for longitudinal pediatric bone accrual

BACKGROUND

Bone accrual impacts lifelong skeletal health, but genetic discovery has been primarily limited to cross-sectional study designs and hampered by uncertainty about target effector genes. Here, we capture this dynamic phenotype by modeling longitudinal bone accrual across 11,000 bone scans in a cohort of healthy children and adolescents, followed by genome-wide association studies (GWAS) and variant-to-gene mapping with functional follow-up.

RESULTS

We identify 40 loci, 35 not previously reported, with various degrees of supportive evidence, half residing in topological associated domains harboring known bone genes. Of several loci potentially associated with later-life fracture risk, a candidate SNP lookup provides the most compelling evidence for rs11195210 (SMC3). Variant-to-gene mapping combining ATAC-seq to assay open chromatin with high-resolution promoter-focused Capture C identifies contacts between GWAS loci and nearby gene promoters. siRNA knockdown of gene expression supports the putative effector gene at three specific loci in two osteoblast cell models. Finally, using CRISPR-Cas9 genome editing, we confirm that the immediate genomic region harboring the putative causal SNP influences PRPF38A expression, a location which is predicted to coincide with a set of binding sites for relevant transcription factors.

CONCLUSIONS

Using a new longitudinal approach, we expand the number of genetic loci putatively associated with pediatric bone gain. Functional follow-up in appropriate cell models finds novel candidate genes impacting bone accrual. Our data also raise the possibility that the cell fate decision between osteogenic and adipogenic lineages is important in normal bone accrual.

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.

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

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

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

Gene Expression and RNA Splicing Imputation Identifies Novel Candidate Genes Associated with Osteoporosis

CONTEXT

Though genome-wide association studies (GWASs) have identified hundreds of genetic variants associated with osteoporosis related traits, such as bone mineral density (BMD) and fracture, it remains a challenge to interpret their biological functions and underlying biological mechanisms.

OBJECTIVE

Integrate diverse expression quantitative trait loci and splicing quantitative trait loci data with several powerful GWAS datasets to identify novel candidate genes associated with osteoporosis.

DESIGN, SETTING, AND PARTICIPANTS

Here, we conducted a transcriptome-wide association study (TWAS) for total body BMD (TB-BMD) (n = 66 628 for discovery and 7697 for validation) and fracture (53 184 fracture cases and 373 611 controls for discovery and 37 857 cases and 227 116 controls for validation), respectively. We also conducted multi-SNP-based summarized mendelian randomization analysis to further validate our findings.

RESULTS

In total, we detected 88 genes significantly associated with TB-BMD or fracture through expression or ribonucleic acid splicing. Summarized mendelian randomization analysis revealed that 78 of the significant genes may have potential causal effects on TB-BMD or fracture in at least 1 specific tissue. Among them, 64 genes have been reported in previous GWASs or TWASs for osteoporosis, such as ING3, CPED1, and WNT16, as well as 14 novel genes, such as DBF4B, GRN, TMUB2, and UNC93B1.

CONCLUSIONS

Overall, our findings provide novel insights into the pathogenesis mechanisms of osteoporosis and highlight the power of a TWAS to identify and prioritize potential causal genes.

Gene regulatory landscape in osteoblast differentiation

The development of osteoblasts, a bone-forming cell population, occurs in conjunction with development of the skeleton, which creates our physical framework and shapes the body. In the past two decades, genetic studies have uncovered the molecular framework of this process-namely, transcriptional regulators and signaling pathways coordinate the cell fate determination and differentiation of osteoblasts in a spatial and temporal manner. Recently emerging genome-wide studies provide additional layers of understanding of the gene regulatory landscape during osteoblast differentiation, allowing us to gain novel insight into the modes of action of the key regulators, functional interaction among the regulator-bound enhancers, epigenetic regulations, and the complex nature of regulatory inputs. In this review, we summarize current understanding of the transcriptional regulation in osteoblasts, in terms of the gene regulatory landscape.

Pregnancy and lactation, a challenge for the skeleton

In this review we discuss skeletal adaptations to the demanding situation of pregnancy and lactation. Calcium demands are increased during pregnancy and lactation, and this is effectuated by a complex series of hormonal changes. The changes in bone structure at the tissue and whole bone level observed during pregnancy and lactation appear to largely recover over time. The magnitude of the changes observed during lactation may relate to the volume and duration of breastfeeding and return to regular menses. Studies examining long-term consequences of pregnancy and lactation suggest that there are small, site-specific benefits to bone density and that bone geometry may also be affected. Pregnancy- and lactation-induced osteoporosis (PLO) is a rare disease for which the pathophysiological mechanism is as yet incompletely known; here, we discuss and speculate on the possible roles of genetics, oxytocin, sympathetic tone and bone marrow fat. Finally, we discuss fracture healing during pregnancy and lactation and the effects of estrogen on this process.

LOX gene polymorphisms are associated with osteoporotic vertebral compression fracture in postmenopausal Chinese women

INTRODUCTION

Collagen cross-linking, which is regulated by lysyl oxidase (LOX), plays critical roles in bone mechanical strength. LOX can influence bone remodeling by modulating osteoblast and osteoclast activity. This study aimed to explore the effect of LOX gene polymorphisms on osteoporotic fractures susceptibility in postmenopausal Chinese women.

METHODS

This was a prospective study of postmenopausal women who visited the outpatient and community clinics of the local Hospital. Five tagging single nucleotide polymorphisms (SNPs) in the LOX gene were determined. Bone mineral density (BMD) was measured at the lumbar spine, femoral neck, and hip using dual-energy X-ray absorptiometry. Fractures were confirmed by X-ray and divided into: vertebral compression fracture (OVCF) and non-OVCF (all other fractures).

RESULTS

This study included 602 patients with non-traumatic fractures and 1343 healthy volunteers. The rs1800449 was significantly associated with vertebral compression fracture (OVCF) after adjusting for age and BMI (P = 0.012). Compared with subjects with the GG genotype, the risk of having OVCF was 1.28 and 1.74, respectively for subjects with the GA and AA genotypes (P = 0.043 and P = 0.018). A recessive genetic model showed that carriers of the AA genotype had higher fracture risk compared to G carriers (GA and GG genotypes) (P = 0.015). The rs2288393 SNP exhibited marginally significant association with OVCF (P = 0.051). Haplotype analyses corroborated our single SNP results: both haplotype CGA and CCG contained rs10519694, rs2288393, and rs1800449, and were significant associated with OVCF (P = 0.048 and P = 0.032, respectively). On the other hand, we found no evidence of an association of LOX gene allelic variants with either BMD or non-OVCF (all P > 0.05).

CONCLUSION

The results suggest that genetic polymorphisms in LOX may contribute to susceptibility to OVCF in Chinese postmenopausal women.

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.

Influence of Genetic Variation in PDE3A on Endothelial Function and Stroke

We aimed to characterize the genetics of endothelial function and how this influences risk for cardiovascular diseases such as ischemic stroke. We integrated genetic data from a study of ultrasound flow-mediated dilatation of brachial artery in adolescents from ALSPAC (Avon Longitudinal Study of Parents and Children; n=5214) with a study of ischemic stroke (MEGASTROKE: n=60 341 cases and 452 969 controls) to identify variants that confer risk of ischemic stroke through altered endothelial function. We identified a variant in PDE3A (Phosphodiesterase 3A), encoding phosphodiesterase 3A, which was associated with flow-mediated dilatation in adolescents (9-12 years of age; β[SE], 0.38 [0.070]; P=3.8×10-8) and confers risk of ischemic stroke (odds ratio, 1.04 [95% CI, 1.02-1.06]; P=5.2×10-6). Bayesian colocalization analyses showed the same underlying variation is likely to lead to both associations (posterior probability, 97%). The same variant was associated with flow-mediated dilatation in a second study in young adults (age, 24-27 years; β[SE], 0.47 [0.23]; P=0.047) but not in older adults (β[SE], -0.012 [0.13]; P=0.89). We conclude that a genetic variant in PDE3A influences endothelial function in early life and leads to increased risk of ischemic stroke. Subtle, measurable changes to the vasculature that are influenced by genetics also influence risk of ischemic stroke.

Influence of Estrogen Receptor Alpha Polymorphism on Bone Mineral Density in Iranian Children

BACKGROUND

Bone mass acquisition in childhood is directly linked to adult bone mineral density (BMD) and fracture risk. BMD is a heritable trait, more than 70% of its variability among a population is affected by genetic factors.

OBJECTIVES

In the present study, we wanted to investigate the association between estrogen receptor alpha (ESR1) polymorphisms, PvuII (rs2234693) and XbaI (rs9340799), and bone area, bone mineral content (BMC), and BMD of the lumbar spine, femoral neck, and also of the total body less the head in Iranian children.

METHODS

The ESR1 gene PvuII and XbaI genotypes were determined by polymerase chain reaction-restriction fragment length polymorphism. Bone area, BMC, BMD, and bone mineral apparent density (BMAD) were assessed by dual-energy X-ray absorptiometry (DEXA). Linear regression was carried out to examine the effects of the ESR1 (PvuII and XbaI) polymorphisms on DEXA outputs when adjusted for confounding factors (i.e., age, sex, BMI, and pubertal stage) in 3 models.

RESULTS

ESR1 (PvuII) gene polymorphisms (CT vs. CC) showed significant effects on the BMC of the total body less the head in all 3 models. For ESR1 (XbaI), individuals with the AG genotype had higher lumbar spine BMD and lumbar spine BMAD compared to other genotypes.

CONCLUSIONS

It seems that the PvuII and XbaI polymorphisms of ESR1 could be associated with BMC and BMD variation in Iranian children and adolescents.

Pilot GWAS of caries in African-Americans shows genetic heterogeneity

BACKGROUND

Dental caries is the most common chronic disease in the US and disproportionately affects racial/ethnic minorities. Caries is heritable, and though genetic heterogeneity exists between ancestries for a substantial portion of loci associated with complex disease, a genome-wide association study (GWAS) of caries specifically in African Americans has not been performed previously.

METHODS

We performed exploratory GWAS of dental caries in 109 African American adults (age > 18) and 96 children (age 3-12) from the Center for Oral Health Research in Appalachia (COHRA1 cohort). Caries phenotypes (DMFS, DMFT, dft, and dfs indices) assessed by dental exams were tested for association with 5 million genotyped or imputed single nucleotide polymorphisms (SNPs), separately in the two age groups. The GWAS was performed using linear regression with adjustment for age, sex, and two principal components of ancestry. A maximum of 1 million adaptive permutations were run to determine empirical significance.

RESULTS

No loci met the threshold for genome-wide significance, though some of the strongest signals were near genes previously implicated in caries such as antimicrobial peptide DEFB1 (rs2515501; p = 4.54 × 10- 6) and TUFT1 (rs11805632; p = 5.15 × 10- 6). Effect estimates of lead SNPs at suggestive loci were compared between African Americans and Caucasians (adults N = 918; children N = 983). Significant (p < 5 × 10- 8) genetic heterogeneity for caries risk was found between racial groups for 50% of the suggestive loci in children, and 12-18% of the suggestive loci in adults.

CONCLUSIONS

The genetic heterogeneity results suggest that there may be differences in the contributions of genetic variants to caries across racial groups, and highlight the critical need for the inclusion of minorities in subsequent and larger genetic studies of caries in order to meet the goals of precision medicine and to reduce oral health disparities.

Meta-Analysis of Genome-Wide Association Studies Identifies Three Loci Associated With Stiffness Index of the Calcaneus

The stiffness index (SI) from quantitative ultrasound measurements is a good indicator of BMD and may be used to predict the risk of osteoporotic fracture. We conducted a genomewide association study (GWAS) for SI using 7742 individuals from the Taiwan Biobank, followed by a replication study in a Korean population (n = 2955). Approximately 6.1 million SNPs were subjected to association analysis, and SI-associated variants were identified. We further conducted a meta-analysis of Taiwan Biobank significant SNPs with a Korean population-based cohort. Candidate genes were prioritized according to epigenetic annotations, gene ontology, protein-protein interaction, GWAS catalog, and expression quantitative trait loci analyses. Our results revealed seven significant single-nucleotide polymorphisms (SNPs) within three loci: 7q31.31, 17p13.3, and 11q14.2. Conditional analysis showed that three SNPs, rs2536195 (CPED1/WNT16), rs1231207 (SMG6), and rs4944661 (LOC10050636/TMEM135), were the most important signals within these regions. The associations for the three SNPs were confirmed in a UK Biobank estimated BMD GWAS; these three cytobands were replicated successfully after a meta-analysis with a Korean population cohort as well. However, two SNPs were not replicated. After prioritization, we identified two novel genes, RAB15 and FNTB, as strong candidates for association with SI. Our study identified three SI-associated SNPs and two novel SI-related genes. Overall, these results provide further insight into the genetic architecture of osteoporosis. Further studies in larger East Asian populations are needed. © 2019 American Society for Bone and Mineral Research.

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.

Genetic Variation in RIN3 in the Belgian Population Supports Its Involvement in the Pathogenesis of Paget's Disease of Bone and Modifies the Age of Onset

Paget's disease of bone (PDB) is a common, late-onset bone disorder characterized by focal increase of bone turnover. Mutations in the SQSTM1 gene are found in up to 40% of patients and recent GWAS have led to novel associations with several loci. RIN3, the candidate gene located at the associated 14q32 locus, has recently been studied in a British cohort to elucidate its contribution to the pathogenesis. In this study, we performed a genetic screening of RIN3 in an unrelated cohort to validate these findings and to further explore genetic variation in this gene in the context of PDB. In our screening, we examined the 5' untranslated region (UTR), the exonic regions and the intron-exon boundaries of the gene in a control cohort and a patient cohort. Our findings show clustering of variation similar to the British cohort and support a protective role for common genetic variation (rs117068593, p.R279C) in the proline-rich region and a functionally relevant role for rare genetic variation in the domains that mediate binding and activation of its interaction partner, Rab5. Additive regression models, fitted for the common variants, validated the association of the rs117068593 variant with the disease (OR^+/+ 0.315; OR^+/- 0.562). In addition, our analyses revealed a potentially modifying effect of this variant on the age of onset of the disease. In conclusion, our findings support the involvement of genetic variation in RIN3 in PDB and suggest a role for RIN3 as a potential modifier of the age of onset of the disease.

High-density genotyping reveals signatures of selection related to acclimation and economically important traits in 15 local sheep breeds from Russia

Background

Domestication and centuries of selective breeding have changed genomes of sheep breeds to respond to environmental challenges and human needs. The genomes of local breeds, therefore, are valuable sources of genomic variants to be used to understand mechanisms of response to adaptation and artificial selection. As a step toward this we performed a high-density genotyping and comprehensive scans for signatures of selection in the genomes from 15 local sheep breeds reared across Russia.

Results

Results demonstrated that the genomes of Russian sheep breeds contain multiple regions under putative selection. More than 50% of these regions matched with intervals identified in previous scans for selective sweeps in sheep genomes. These regions contain well-known candidate genes related to morphology, adaptation, and domestication (e.g., KITLG, KIT, MITF, and MC1R), wool quality and quantity (e.g., DSG@, DSC@, and KRT@), growth and feed intake (e.g., HOXA@, HOXC@, LCORL, NCAPG, LAP3, and CCSER1), reproduction (e.g., CMTM6, HTRA1, GNAQ, UBQLN1, and IFT88), and milk-related traits (e.g., ABCG2, SPP1, ACSS1, and ACSS2). In addition, multiple genes that are putatively related to environmental adaptations were top-ranked in selected intervals (e.g., EGFR, HSPH1, NMUR1, EDNRB, PRL, TSHR, and ADAMTS5). Moreover, we observed that multiple key genes involved in human hereditary sensory and autonomic neuropathies, and genetic disorders accompanied with an inability to feel pain and environmental temperatures, were top-ranked in multiple or individual sheep breeds from Russia pointing to a possible mechanism of adaptation to harsh climatic conditions.

Conclusions

Our work represents the first comprehensive scan for signatures of selection in genomes of local sheep breeds from the Russian Federation of both European and Asian origins. We confirmed that the genomes of Russian sheep contain previously identified signatures of selection, demonstrating the robustness of our integrative approach. Multiple novel signatures of selection were found near genes which could be related to adaptation to the harsh environments of Russia. Our study forms a basis for future work on using Russian sheep genomes to spot specific genetic variants or haplotypes to be used in efforts on developing next-generation highly productive breeds, better suited to diverse Eurasian environments.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5537-0) contains supplementary material, which is available to authorized users.

Genome-scale Capture C promoter interactions implicate effector genes at GWAS loci for bone mineral density

Osteoporosis is a devastating disease with an essential genetic component. GWAS have discovered genetic signals robustly associated with bone mineral density (BMD), but not the precise localization of effector genes. Here, we carry out physical and direct variant to gene mapping in human mesenchymal progenitor cell-derived osteoblasts employing a massively parallel, high resolution Capture C based method in order to simultaneously characterize the genome-wide interactions of all human promoters. By intersecting our Capture C and ATAC-seq data, we observe consistent contacts between candidate causal variants and putative target gene promoters in open chromatin for ~ 17% of the 273 BMD loci investigated. Knockdown of two novel implicated genes, ING3 at ‘CPED1-WNT16’ and EPDR1 at ‘STARD3NL’, inhibits osteoblastogenesis, while promoting adipogenesis. This approach therefore aids target discovery in osteoporosis, here on the example of two relevant genes involved in the fate determination of mesenchymal progenitors, and can be applied to other common genetic diseases.

GWAS have identified numerous genetic loci for bone mineral density (BMD) and fracture risk. Here, the authors map these variants to putative target genes using ATAC-seq and Capture C of human osteoblasts and confirm ING3 and EPDR1 as BMD genes in in vitro osteoblast differentiation experiments.

LncRNA ZBTB40-IT1 modulated by osteoporosis GWAS risk SNPs suppresses osteogenesis

Previous genome-wide linkage and association studies have identified an osteoporosis-associated locus at 1p36 that harbors SNPs rs34920465 and rs6426749. The 1p36 locus also comprises the WNT4 gene with known role in bone metabolism and functionally unknown ZBTB40/lncRNA ZBTB40-IT1 genes. How these might interact to contribute to osteoporosis susceptibility is not known. In this study, we show that lncRNA ZBTB40-IT1 is able to suppress osteogenesis and promote osteoclastogenesis by regulating the expression of WNT4, RUNX2, OSX, ALP, COL1A1, OPG and RANKL in U-2OS and hFOB1.19 cell lines, whereas ZBTB40 plays an opposite role in bone metabolism. Treatment with parathyroid hormone significantly downregulates the expression of ZBTB40-IT1 in U-2OS cell lines. ZBTB40 can suppress ZBTB40-IT1 expression but has no response to parathyroid hormone treatment. Dual-luciferase assay and biotin pull-down assay demonstrate that osteoporosis GWAS lead SNPs rs34920465-G and rs6426749-C alleles can respectively bind transcription factors JUN::FOS and CREB1, and upregulate ZBTB40 and ZBTB40-IT1 expression. Our study discovers the critical role of ZBTB40 and lncRNA ZBTB40-IT1 in bone metabolism, and provides a mechanistic basis for osteoporosis GWAS lead SNPs rs34920465 and rs6426749.