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

Associations of sleep behaviors and genetic risk with risk of incident osteoporosis: A prospective cohort study of 293,164 participants

BACKGROUND

Unhealthy sleep behaviors are associated with higher risks of osteoporosis (OP), while prospective evidence is limited. This study aimed to prospectively investigate this association, quantify the attributable burden of OP incidence reduction due to unhealthy sleep behaviors, and explore potential modifications by genetic risk factors.

METHODS

This longitudinal cohort study was conducted utilizing data from the UK Biobank, comprising 293,164 participants initially free of OP and with requisite sleep behaviors data at baseline. We followed the participants after recruitment until November 30, 2022, to ascertain incident OP. We assessed the associations of five sleep behaviors including sleep duration, chronotype, insomnia, daytime napping, and morning wake-up difficulties, as well as sleep behavior patterns identified based on the above sleep behaviors, with the risk of OP, using Cox models adjusted for multiple confounders. The analyses were then performed separately among individuals with different OP susceptibility, indexed by standard polygenetic risk scores(PRS) for OP. Our secondary outcome was OP with pathologic fracture. Subgroup and sensitivity analyses were performed. Additionally, attributable risk percent in the exposed population (AR%) and population attributable fraction (PAF) of sleep behaviors were calculated.

RESULTS

Over a median follow-up of 13.7 years, 8253 new-onset OP cases were documented. Unhealthy sleep behaviors, such as long or short sleep duration, insomnia, daytime napping, morning wake-up difficulties, and unhealthy sleep patterns, were associated with elevated risks of OP (HRs ranging from 1.14 to 1.46, all P-value <0.001) compared to healthy sleep behaviors. Similar associations were observed for OP with pathologic fractures. Insomnia exhibited the largest AR% of 39.98 % (95%CI: 36.46, 43.31) and PAF of 33.25 % (95%CI: 30.00, 36.34) among healthy sleep patterns and components. A statistically significant multiplicative interaction was noted between sleep behaviors and OP PRS on OP risk (all P-interaction <0.001).

CONCLUSIONS

Four unhealthy sleep behaviors and sleep behavior patterns were associated to increased OP risk, with insomnia contributing the most to OP incidence, while genetic risk for OP modified this association. These findings underscore the crucial role of adhering to healthy sleep behaviors for effective OP prevention.

Genetic and Epigenetic Landscape for Drug Development in Polycystic Ovary Syndrome

The treatment of polycystic ovary syndrome (PCOS) faces challenges as all known treatments are merely symptomatic. The US Food and Drug Administration has not approved any drug specifically for treating PCOS. As the significance of genetics and epigenetics rises in drug development, their pivotal insights have greatly enhanced the efficacy and success of drug target discovery and validation, offering promise for guiding the advancement of PCOS treatments. In this context, we outline the genetic and epigenetic advancement in PCOS, which provide novel insights into the pathogenesis of this complex disease. We also delve into the prospective method for harnessing genetic and epigenetic strategies to identify potential drug targets and ensure target safety. Additionally, we shed light on the preliminary evidence and distinctive challenges associated with gene and epigenetic therapies in the context of PCOS.

Unraveling the molecular and immunological landscape: Exploring signaling pathways in osteoporosis

Osteoporosis, characterized by compromised bone density and microarchitecture, represents a significant global health challenge, particularly in aging populations. This comprehensive review delves into the intricate signaling pathways implicated in the pathogenesis of osteoporosis, providing valuable insights into the pivotal role of signal transduction in maintaining bone homeostasis. The exploration encompasses cellular signaling pathways such as Wnt, Notch, JAK/STAT, NF-κB, and TGF-β, all of which play crucial roles in bone remodeling. The dysregulation of these pathways is a contributing factor to osteoporosis, necessitating a profound understanding of their complexities to unveil the molecular mechanisms underlying bone loss. The review highlights the pathological significance of disrupted signaling in osteoporosis, emphasizing how these deviations impact the functionality of osteoblasts and osteoclasts, ultimately resulting in heightened bone resorption and compromised bone formation. A nuanced analysis of the intricate crosstalk between these pathways is provided to underscore their relevance in the pathophysiology of osteoporosis. Furthermore, the study addresses some of the most crucial long non-coding RNAs (lncRNAs) associated with osteoporosis, adding an additional layer of academic depth to the exploration of immune system involvement in various types of osteoporosis. Finally, we propose that SKP1 can serve as a potential biomarker in osteoporosis.

Pre-hypertrophic chondrogenic enhancer landscape of limb and axial skeleton development

Chondrocyte differentiation controls skeleton development and stature. Here we provide a comprehensive map of chondrocyte-specific enhancers and show that they provide a mechanistic framework through which non-coding genetic variants can influence skeletal development and human stature. Working with fetal chondrocytes isolated from mice bearing a Col2a1 fluorescent regulatory sensor, we identify 780 genes and 2'704 putative enhancers specifically active in chondrocytes using a combination of RNA-seq, ATAC-seq and H3K27ac ChIP-seq. Most of these enhancers (74%) show pan-chondrogenic activity, with smaller populations being restricted to limb (18%) or trunk (8%) chondrocytes only. Notably, genetic variations overlapping these enhancers better explain height differences than those overlapping non-chondrogenic enhancers. Finally, targeted deletions of identified enhancers at the Fgfr3, Col2a1, Hhip and, Nkx3-2 loci confirm their role in regulating cognate genes. This enhancer map provides a framework for understanding how genes and non-coding variations influence bone development and diseases.

Combined application of mesenchymal stem cells and different glucocorticoid dosing alleviates osteoporosis in SLE murine models

OBJECTIVE

Bone mesenchymal stem cells (BMSCs) have been tentatively applied in the treatment of glucocorticoid-induced osteoporosis (GIOP) and systemic lupus erythematosus (SLE). However, the effects of BMSCs on osteoporosis within the context of glucocorticoid (GC) application in SLE remain unclear. Our aim was to explore the roles of BMSCs and different doses of GC interventions on osteoporosis in SLE murine models.

METHODS

MRL/MpJ-Faslpr mice were divided into eight groups with BMSC treatment and different dose of GC intervention. Three-dimensional imaging analysis and hematoxylin and eosin (H&E) staining were performed to observe morphological changes. The concentrations of osteoprotegerin (OPG) and receptor activator of nuclear factor κB ligand (RANKL) in serum were measured by enzyme-linked immunosorbent assay (ELISA). The subpopulation of B cells and T cells in bone marrows and spleens were analyzed by flow cytometry. Serum cytokines and chemokines were assessed using Luminex magnetic bead technology.

RESULTS

BMSCs ameliorated osteoporosis in murine SLE models by enhancing bone mass, improving bone structure, and promoting bone formation through increased bone mineral content and optimization of trabecular morphology. BMSC and GC treatments reduced the number of B cells in bone marrows, but the effect was not significant in spleens. BMSCs significantly promoted the expression of IL-10 while reducing IL-18. Moreover, BMSCs exert immunomodulatory effects by reducing Th17 expression and rectifying the Th17/Treg imbalance.

CONCLUSION

BMSCs effectively alleviate osteoporosis induced by SLE itself, as well as osteoporosis resulting from SLE combined with various doses of GC therapy. The therapeutic effects of BMSCs appear to be mediated by their influence on bone marrow B cells, T cell subsets, and associated cytokines. High-dose GC treatment exerts a potent anti-inflammatory effect but may hinder the immunotherapeutic potential of BMSCs. Our research may offer valuable guidance to clinicians regarding the use of BMSC treatment in SLE and provide insights into the judicious use of GCs in clinical practice.

Understanding Osteoporosis: Human Bone Density, Genetic Mechanisms, Gut Microbiota, and Future Prospects

Osteoporosis is a systemic condition of the skeleton that leads to diminished bone mass, a breakdown in the bone tissue's microscopic architecture, and an elevated risk of breaking a bone. The elderly and women particularly after menopause are disproportionately affected, and the condition generally stays undiagnosed until a broken bone causes severe pain and immobility. Causes of osteoporosis include low bone mass, more than normal bone loss, changes in hormone levels (decreased estrogen or testosterone), certain diseases and therapies, and lifestyle factors like smoking and inactivity. The spine, hip, and forearm are particularly vulnerable to osteoporosis-related fractures. The purpose of this article is to present a thorough understanding of osteoporosis, including the disease's connection to bone density in humans, and the major part played by genetic pathways and gut flora. The causes of osteoporosis, the effects of aging on bone density, and why some groups experience a higher incidence of the disease than others are investigated. The paper also includes animal and human experiments investigating the link between gut flora and osteoporosis. Finally, it looks to the future and speculates on possible developments in osteoporosis prevention and therapy.

Common and novel haplotype structures between different types of cancer

BACKGROUND

Background: Genome-wide association studies (GWAS) have identified hundreds of genetic variants associated with cancer risk. GWAS data are important for cancer prevention and understanding the underlying mechanisms of cancer.

AIMS

This study aimed to investigate the genetic association between different types of cancer using GWAS data and a bioinformatics approach.

METHODS AND RESULTS

The significant GWAS variants associated with more than one cancer type were identified. Common linkage disequilibrium (LD) variants between different types of cancer were identified by 1000 genomes phase 3 LD data. Haplotype blocks were identified by analyzing 1000 Genomes phase 3 genotyping data in the GWAS populations. Subsequent analyses included functional SNP analyses and TCGA gene expression. The results associated with significant GWAS variants (P<5E-8) showed the following haplotype associations in European population: GT rs4808075-rs8170 haplotype on BABAM1 with breast and ovarian cancers, GC rs16857609-rs11693806 haplotype on DIRC3 with breast and thyroid cancers, GCG rs380286-rs401681-rs31487 haplotype on CLPTM1L with skin and lung cancers, GGG rs4430796-rs11651052-rs11263763 haplotype on HNF1B with prostate and endometrial cancers, and GT rs10505477-rs6983267 haplotype on CASC8 associated with colorectal and prostate cancers. All these genes had significantly different expressions in tumor tissues (P<1E-3). In addition, the rs11693806 variant is located in the hsa-miR-873-5p binding site and has an enhancing effect on the hsa-miR-873-5p:DIRC3 interaction.

CONCLUSION

These novel haplotype structures and miRNA:lncRNA interactions are important for understanding the common genetic link between cancers. These results can potentially be used in genetic panels.

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.

A comprehensive review and advanced biomolecule-based therapies for osteoporosis

BACKGROUND

The prevalence of osteoporosis (OP) on a global scale is significantly elevated that causes life threatening issues. The potential of groundbreaking biomolecular therapeutics in the field of OP is highly encouraging. The administration of biomolecular agents has the potential to mitigate the process of bone demineralization while concurrently augmenting the regenerative capacity of bone tissue, thereby facilitating a personalized therapeutic approach. Biomolecules-based therapies showed promising results in term of bone mass protection and restoration in OP.

AIM OF REVIEW

We summarized the recent biomolecular therapies with notable progress in clinical, demonstrating the potential to transform illness management. These treatments frequently utilize different biomolecule based strategies. Biomolecular therapeutics has a targeted character, which results in heightened specificity and less off-target effects, ultimately leading to increased patient outcomes. These aspects have the capacity to greatly enhance the management of OP, thus resulting in a major enhancement in the quality of life encountered by individuals affected by this condition.

Enhanced fracture risk prediction: a novel multi-trait genetic approach integrating polygenic scores of fracture-related traits

The novel metaPGS, integrating multiple fracture-related genetic traits, surpasses traditional polygenic scores in predicting fracture risk. Demonstrating a robust association with incident fractures, this metaPGS offers significant potential for enhancing clinical fracture risk assessment and tailoring prevention strategies.

INTRODUCTION

Current polygenic scores (PGS) have limited predictive power for fracture risk. To improve genetic prediction, we developed and evaluated a novel metaPGS combining genetic information from multiple fracture-related traits.

METHODS

We derived individual PGS from genome-wide association studies of 16 fracture-related traits and employed an elastic-net logistic regression model to examine the association between the 16 PGSs and fractures. An optimal metaPGS was constructed by combining 11 significant individual PGSs selected by the elastic regularized regression model. We evaluated the predictive power of the metaPGS alone and in combination with clinical risk factors recommended by guidelines. The discrimination ability of metaPGS was assessed using the concordance index. Reclassification was assessed using net reclassification improvement (NRI) and integrated discrimination improvement (IDI).

RESULTS

The metaPGS had a significant association with incident fractures (HR 1.21, 95% CI 1.18-1.25 per standard deviation of metaPGS), which was stronger than previously developed bone mineral density (BMD)-related individual PGSs. Models with PGS_FNBMD, PGS_TBBMD, and metaPGS had slightly higher but statistically non-significant c-index than the base model (0.640, 0.644, 0.644 vs. 0.638). However, the reclassification analysis showed that compared to the base model, the model with metaPGS improves the reclassification of fracture.

CONCLUSIONS

The metaPGS is a promising approach for stratifying fracture risk in the European population, improving fracture risk prediction by combining genetic information from multiple fracture-related traits.

Enhanced osteoporotic fracture prediction in postmenopausal women using Bayesian optimization of machine learning models with genetic risk score

This study aimed to enhance the fracture risk prediction accuracy in major osteoporotic fractures (MOFs) and hip fractures (HFs) by integrating genetic profiles, machine learning (ML) techniques, and Bayesian optimization. The genetic risk score (GRS), derived from 1,103 risk single nucleotide polymorphisms (SNPs) from genome-wide association studies (GWAS), was formulated for 25,772 postmenopausal women from the Women's Health Initiative dataset. We developed four ML models: Support Vector Machine (SVM), Random Forest, XGBoost, and Artificial Neural Network (ANN) for binary fracture outcome and 10-year fracture risk prediction. GRS and FRAX clinical risk factors (CRFs) were used as predictors. Death as a competing risk was accounted for in ML models for time-to-fracture data. ML models were subsequently fine-tuned through Bayesian optimization, which displayed marked superiority over traditional grid search. Evaluation of the models' performance considered an array of metrics such as accuracy, weighted F1 Score, the area under the precision-recall curve (PRAUC), and the area under the receiver operating characteristic curve (AUC) for binary fracture predictions, and the C-index, Brier score, and dynamic mean AUC over a 10-year follow-up period for fracture risk predictions. We found that GRS-integrated XGBoost with Bayesian optimization is the most effective model, with an accuracy of 91.2% (95% CI: 90.4-92.0%) and an AUC of 0.739 (95% CI: 0.731-0.746) in MOF binary predictions. For 10-year fracture risk modeling, the XGBoost model attained a C-index of 0.795 (95% CI: 0.783-0.806) and a mean dynamic AUC of 0.799 (95% CI: 0.788-0.809). Compared to FRAX, the XGBoost model exhibited a categorical net reclassification improvement (NRI) of 22.6% (P = .004). A sensitivity analysis, which included BMD but lacked GRS, reaffirmed these findings. Furthermore, portability tests in diverse non-European groups, including Asians and African Americans, underscored the model's robustness and adaptability. This study accentuates the potential of combining genetic insights and optimized ML in strengthening fracture predictions, heralding new preventive strategies for postmenopausal women.

Osteoporosis GWAS-implicated DNM3 locus contextually regulates osteoblastic and chondrogenic fate of mesenchymal stem/progenitor cells through oscillating miR-199a-5p levels

Genome wide association study (GWAS)-implicated bone mineral density (BMD) signals have been shown to localize in cis-regulatory regions of distal effector genes using 3D genomic methods. Detailed characterization of such genes can reveal novel causal genes for BMD determination. Here, we elected to characterize the "DNM3" locus on chr1q24, where the long non-coding RNA DNM3OS and the embedded microRNA MIR199A2 (miR-199a-5p) are implicated as effector genes contacted by the region harboring variation in linkage disequilibrium with BMD-associated sentinel single nucleotide polymorphism, rs12041600. During osteoblast differentiation of human mesenchymal stem/progenitor cells (hMSC), miR-199a-5p expression was temporally decreased and correlated with the induction of osteoblastic transcription factors RUNX2 and Osterix. Functional relevance of miR-199a-5p downregulation in osteoblastogenesis was investigated by introducing miR-199a-5p mimic into hMSC. Cells overexpressing miR-199a-5p depicted a cobblestone-like morphological change and failed to produce BMP2-dependent extracellular matrix mineralization. Mechanistically, a miR-199a-5p mimic modified hMSC propagated normal SMAD1/5/9 signaling and expressed osteoblastic transcription factors RUNX2 and Osterix but depicted pronounced upregulation of SOX9 and enhanced expression of essential chondrogenic genes ACAN, COMP, and COL10A1. Mineralization defects, morphological changes, and enhanced chondrogenic gene expression associated with miR-199a-5p mimic over-expression were restored with miR-199a-5p inhibitor suggesting specificity of miR-199a-5p in chondrogenic fate specification. The expression of both the DNM3OS and miR-199a-5p temporally increased and correlated with hMSC chondrogenic differentiation. Although miR-199a-5p overexpression failed to further enhance chondrogenesis, blocking miR-199a-5p activity significantly reduced chondrogenic pellet size, extracellular matrix deposition, and chondrogenic gene expression. Taken together, our results indicate that oscillating miR-199a-5p levels dictate hMSC osteoblast or chondrocyte terminal fate. Our study highlights a functional role of miR-199a-5p as a BMD effector gene at the DNM3 BMD GWAS locus, where patients with cis-regulatory genetic variation which increases miR-199a-5p expression could lead to reduced osteoblast activity.

Identification of an novel genetic variant associated with osteoporosis: insights from the Taiwan Biobank Study

Purpose

The purpose of this study was to identify new independent significant SNPs associated with osteoporosis using data from the Taiwan Biobank (TWBB).

Material and Methods

The dataset was divided into discovery (60%) and replication (40%) subsets. Following data quality control, genome-wide association study (GWAS) analysis was performed, adjusting for sex, age, and the top 5 principal components, employing the Scalable and Accurate Implementation of the Generalized mixed model approach. This was followed by a meta-analysis of TWBB1 and TWBB2. The Functional Mapping and Annotation (FUMA) platform was used to identify osteoporosis-associated loci. Manhattan and quantile-quantile plots were generated using the FUMA platform to visualize the results. Independent significant SNPs were selected based on genome-wide significance (P < 5 × 10-8) and independence from each other (r2 < 0.6) within a 1 Mb window. Positional, eQTL(expression quantitative trait locus), and Chromatin interaction mapping were used to map SNPs to genes.

Results

A total of 29 084 individuals (3154 osteoporosis cases and 25 930 controls) were used for GWAS analysis (TWBB1 data), and 18 918 individuals (1917 cases and 17 001 controls) were utilized for replication studies (TWBB2 data). We identified a new independent significant SNP for osteoporosis in TWBB1, with the lead SNP rs76140829 (minor allele frequency = 0.055, P-value = 1.15 × 10-08). Replication of the association was performed in TWBB2, yielding a P-value of 6.56 × 10-3. The meta-analysis of TWBB1 and TWBB2 data demonstrated a highly significant association for SNP rs76140829 (P-value = 7.52 × 10-10). In the positional mapping of rs76140829, 6 genes (HABP2, RP11-481H12.1, RNU7-165P, RP11-139 K1.2, RP11-57H14.3, and RP11-214 N15.5) were identified through chromatin interaction mapping in mesenchymal stem cells.

Conclusions

Our GWAS analysis using the Taiwan Biobank dataset unveils rs76140829 in the VTI1A gene as a key risk variant associated with osteoporosis. This finding expands our understanding of the genetic basis of osteoporosis and highlights the potential regulatory role of this SNP in mesenchymal stem cells.

Deciphering the complex relationship between type 2 diabetes and fracture risk with both genetic and observational evidence

The 'diabetic bone paradox' suggested that type 2 diabetes (T2D) patients would have higher areal bone mineral density (BMD) but higher fracture risk than individuals without T2D. In this study, we found that the genetically predicted T2D was associated with higher BMD and lower risk of fracture in both weighted genetic risk score (wGRS) and two-sample Mendelian randomization (MR) analyses. We also identified ten genomic loci shared between T2D and fracture, with the top signal at SNP rs4580892 in the intron of gene RSPO3. And the higher expression in adipose subcutaneous and higher protein level in plasma of RSPO3 were associated with increased risk of T2D, but decreased risk of fracture. In the prospective study, T2D was observed to be associated with higher risk of fracture, but BMI mediated 30.2% of the protective effect. However, when stratified by the T2D-related risk factors for fracture, we observed that the effect of T2D on the risk of fracture decreased when the number of T2D-related risk factors decreased, and the association became non-significant if the T2D patients carried none of the risk factors. In conclusion, the genetically determined T2D might not be associated with higher risk of fracture. And the shared genetic architecture between T2D and fracture suggested a top signal around RSPO3 gene. The observed effect size of T2D on fracture risk decreased if the T2D-related risk factors could be eliminated. Therefore, it is important to manage the complications of T2D to prevent the risk of fracture.

The Relationship Between Systemic Lupus Erythematosus and Osteoporosis Based on Different Ethnic Groups: a Two-Sample Mendelian Randomization Analysis

The previous observational studies could not overcome the effects of confounding variables and reverse causality. We aimed to determine whether there is a causal relationship between systemic lupus erythematosus and osteoporosis in East Asian and European populations, respectively, by two-sample Mendelian Randomization analysis. We obtained and downloaded data from publicly available genome-wide association study databases and analyses for East Asian and European populations, including systemic lupus erythematosus (SLE), osteoporosis (OP), multisite bone mineral density (BMD), and OP with fracture. After screening for instrumental single-nucleotide polymorphisms (SNPs) significantly correlated to SLE, the inverse-variance weighted (IVW) method was used for calculating the ratio and 95% confidence interval, besides utilizing MR-Egger, weighted median, and weighted mode to assess the robustness of the primary outcome. Moreover, multiple analyses, including MR-PRESSO, MR-Egger intercept, Cochran's Q test, as well as "leave-one-out" sensitivity, were used for evaluating horizontal pleiotropy, heterogeneity, and stability. Finally, we exchanged exposure and outcome and performed a reverse MR analysis. IVW (OR = 1.05, 95% CI = 1.01-1.09, P = 0.009) indicated a significant positive correlation between genetically predicted SLE and OP in East Asians. Furthermore, neither heterogeneity nor horizontal pleiotropy was observed. In Europe, there was no significant genetically predicted causal relation between SLE and OP. Bi-directional MR analysis showed no reverse causality between SLE and OP. In the East Asian population, genetically predicted SLE may have had a positive causal relationship with OP. In Europe, there is insufficient evidence for a potential causal relation between SLE and OP or BMD and fracture, and the correlations currently observed may be attributed to a variety of confounder variables.

Garcinone C attenuates RANKL-induced osteoclast differentiation and oxidative stress by activating Nrf2/HO-1 and inhibiting the NF-kB signaling pathway

Osteoporosis is the result of osteoclast formation exceeding osteoblast production, and current osteoporosis treatments targeting excessive osteoclast bone resorption have serious adverse effects. There is a need to fully understand the mechanisms of osteoclast-mediated bone resorption, identify new drug targets, and find better drugs to treat osteoporosis. Gar C (Gar C) is a major naturally occurring phytochemical isolated from mangosteen, and is a derivative of the naturally occurring phenolic antioxidant lutein. We used an OP mouse model established by ovariectomy (OVX). We found that treatment with Gar C significantly increased bone mineral density and significantly decreased the expression of TRAP, NFATC1 and CTSK relative to untreated OP mice. We found that Garcinone C could disrupt osteoclast activation and resorption functions by inhibiting RANKL-induced osteoclast differentiation as well as inhibiting the formation of multinucleated osteoclasts. Immunoblotting showed that Gar C downregulated the expression of osteoclast-related proteins. In addition, Gar C significantly inhibited RANKL-induced ROS production and affected NF-κB activity by inhibiting phosphorylation Formylation of P65 and phosphorylation and degradation of ikba. These data suggest that Gar C significantly reduced OVX-induced osteoporosis by inhibiting osteoclastogenesis and oxidative stress in bone tissue. Mechanistically, this effect was associated with inhibition of the ROS-mediated NF-κB pathway.

Association between Fat Mass and Obesity-Related Transcript Polymorphisms and Osteoporosis Phenotypes

BACKGROUND

Common variants in the fat mass and obesity-related transcript (FTO) gene are related to body mass index and obesity, suggesting its potential association with bone mineral density (BMD) and fracture risk. This study sought to define the association between FTO gene variants and the following phenotypes: (1) BMD; (2) bone loss; and (3) fracture risk.

METHODS

This analysis was based on the Dubbo Osteoporosis Epidemiology Study that included 1,277 postmenopausal women aged ≥60 years living in Dubbo, Australia. BMD at the femoral neck and lumbar spine was measured biennially by dual energy X-ray absorptiometry (GE Lunar). Fractures were radiologically ascertained. Six single nucleotide polymorphisms (SNPs; rs1421085, rs1558902, rs1121980, rs17817449, rs9939609, and rs9930506) of the FTO gene were genotyped using TaqMan assay.

RESULTS

Women homozygous for the minor allele (GG) of rs9930506 had a significantly higher risk of hip fracture (adjusted hazard ratio, 1.93; 95% confidence interval, 1.15-3.23) than those homozygous for the major allele (AA) after adjusting for potential confounding effects. Similar associations were also observed for the minor allele of rs1121980. However, there was no significant association between the FTO SNPs and BMD or the rate of bone loss.

CONCLUSIONS

Common variations in the FTO gene are associated with a hip fracture risk in women, and the association is not mediated through BMD or bone loss.

A Missense Variant in TP53 Could Be a Genetic Biomarker Associated with Bone Tissue Alterations

Metabolic bone diseases cover a broad spectrum of disorders that share alterations in bone metabolism that lead to a defective skeleton, which is associated with increasing morbidity, disability, and mortality. There is a close connection between the etiology of metabolic bone diseases and genetic factors, with TP53 being one of the genes associated therewith. The single nucleotide polymorphism (SNP) Arg72Pro of TP53 is a genetic factor associated with several pathologies, including cancer, stroke, and osteoporosis. Here, we aim to analyze the influence of the TP53 Arg72Pro SNP on bone mass in humanized Tp53 Arg72Pro knock-in mice. This work reports on the influence of the TP53 Arg72Pro polymorphism in bone microarchitecture, OPG expression, and apoptosis bone status. The results show that the proline variant of the TP53 Arg72Pro polymorphism (Pro72-p53) is associated with deteriorated bone tissue, lower OPG/RANK ratio, and lower apoptosis in bone tissue. In conclusion, the TP53 Arg72Pro polymorphism modulates bone microarchitecture and may be a genetic biomarker that can be used to identify individuals with an increased risk of suffering metabolic bone alterations.

Systems genetics and bioinformatics analyses using ESR1-correlated genes identify potential candidates underlying female bone development

Estrogen receptor α (ESR1) is involved in E2 signaling and plays a major role in postmenopausal bone loss. However, the molecular network underlying ESR1 has not been explored. We used systems genetics and bioinformatics to identify important genes associated with Esr1 in postmenopausal bone loss. We identified ~2300 Esr1-coexpressed genes in female BXD bone femur, functional analysis of which revealed 'osteoblast signaling' as the most enriched pathway. PPI network led to the identification of 25 'female bone candidates'. The gene-regulatory analysis revealed RUNX2 as a key TF. ANKRD1 and RUNX2 were significantly different between osteoporosis patients and healthy controls. Sp7, Col1a1 and Pth1r correlated with multiple femur bone phenotypes in BXD mice. miR-3121-3p targeted Csf1, Ankrd1, Sp7 and Runx2. β-estradiol treatment markedly increased the expression of these candidates in mouse osteoblast. Our study revealed that Esr1-correlated genes Ankrd1, Runx2, Csf1 and Sp7 may play important roles in female bone development.

Genetic and Gene Expression Resources for Osteoporosis and Bone Biology Research

PURPOSE OF REVIEW

The integration of data from multiple genomic assays from humans and non-human model organisms is an effective approach to identify genes involved in skeletal fragility and fracture risk due to osteoporosis and other conditions. This review summarizes genome-wide genetic variation and gene expression data resources relevant to the discovery of genes contributing to skeletal fragility and fracture risk.

RECENT FINDINGS

Genome-wide association studies (GWAS) of osteoporosis-related traits are summarized, in addition to gene expression in bone tissues in humans and non-human organisms, with a focus on rodent models related to skeletal fragility and fracture risk. Gene discovery approaches using these genomic data resources are described. We also describe the Musculoskeletal Knowledge Portal (MSKKP) that integrates much of the available genomic data relevant to fracture risk. The available genomic resources provide a wealth of knowledge and can be analyzed to identify genes related to fracture risk. Genomic resources that would fill particular scientific gaps are discussed.

Exosome-targeted delivery of METTL14 regulates NFATc1 m6A methylation levels to correct osteoclast-induced bone resorption

Osteoporosis has a profound influence on public health. First-line bisphosphonates often cause osteonecrosis of the jaw meanwhile inhibiting osteoclasts. Therefore, it is important to develop effective treatments. The results of this study showed that the increased level of NFATc1 m6A methylation caused by zoledronic acid (ZOL), with 4249A as the functional site, is highly correlated with the decreased bone resorption of osteoclasts. Upstream, METTL14 regulates osteoclast bone absorption through the methylation functional site of NFATc1. Downstream, YTHDF1 and YTHDF2 show antagonistic effects on the post-transcriptional regulation of NFATc1 after the m6A methylation level is elevated by METTL14. In this study, meRIP-Seq, luciferase reporter assays, meRIP and other methods were used to elucidate the NFATc1 regulatory mechanism of osteoclasts from the perspective of RNA methylation. In addition, EphA2 overexpression on exosomes is an effective biological method for targeted delivery of METTL14 into osteoclasts. Importantly, this study shows that METTL14 released by exosomes can increase the m6A methylation level of NFATc1 to inhibit osteoclasts, help postmenopausal osteoporosis patients preserve bone mass, and avoid triggering osteonecrosis of the jaw, thus becoming a new bioactive molecule for the treatment of osteoporosis.

The role of genetic and epigenetic factors in determining the risk of spinal fragility fractures: new insights in the management of spinal osteoporosis

Osteoporosis predisposes patients to spinal fragility fractures. Imaging plays a key role in the diagnosis and prognostication of these osteoporotic vertebral fractures (OVF). However, the current imaging knowledge base for OVF is lacking sufficient standardisation to enable effective risk prognostication. OVF have been shown to be more prevalent in Caucasian patient cohorts in comparison to the Eastern Asian population. These population-based differences in risk for developing OVF suggest that there could be genetic and epigenetic factors that drive the pathogenesis of osteoporosis, low bone mineral density (BMD) and OVF. Several genetic loci have been associated with a higher vertebral fracture risk, although at varying degrees of significance. The present challenge is clarifying whether these associations are specific to vertebral fractures or osteoporosis more generally. Furthermore, these factors could be exploited for diagnostic interpretation as biomarkers [including novel long non-coding (lnc)RNAs, micro (mi)RNAs and circular (circ)RNAs]. The extent of methylation of genes, alongside post-translational histone modifications, have shown to affect several interlinked pathways that converge on the regulation of bone deposition and resorption, partially through their influence on osteoblast and osteoclast differentiation. Lastly, in addition to biomarkers, several exciting new imaging modalities could add to the established dual-energy X-ray absorptiometry (DXA) method used for BMD assessment. New technologies, and novel sequences within existing imaging modalities, may be able to quantify the quality of bone in addition to the BMD and bone structure; these are making progress through various stages of development from the pre-clinical sphere through to deployment in the clinical setting. In this mini review, we explore the literature to clarify the genetic and epigenetic factors associated with spinal fragility fractures and delineate the causal genes, pathways and interactions which could drive different risk profiles. We also outline the cutting-edge imaging modalities which could transform diagnostic protocols for OVF.

Inhibition of Caspase-1-mediated pyroptosis promotes osteogenic differentiation, offering a therapeutic target for osteoporosis

BACKGROUND

Pyroptosis, an emerging inflammatory form of cell death, has been previously demonstrated to stimulate a massive inflammatory response, thus hindering the osteogenic differentiation of bone marrow mesenchymal stem cells (BMSCs). Nevertheless, the impact of pyroptosis in thwarting osteogenic differentiation and exacerbating the advancement of osteoporosis (OP) remains enigmatic.

METHODS

We evaluated the expression levels of pyroptosis-associated indicators, including NOD-like receptor family pyrin domain-containing protein 3 (NLRP3), CASPASE-1, IL-1β, and IL-18, in specimens obtained from femoral heads of OP patients, as well as in an ovariectomy-induced mouse model of OP. Subsequently, the precise roles of pyroptosis in osteogenic differentiation were investigated using bioinformatics analysis, alongside morphological and biochemical assessments.

RESULTS

The pivotal pyroptotic proteins, including NLRP3, Caspase-1, IL-1β, and IL-18, exhibited significant upregulation within the bone tissue samples of clinical OP cases, as well as in the femoral tissues of ovariectomy (OVX)-induced mouse OP model, displaying a negatively associated with compromised osteogenic capacity, as represented by lessened bone mass, suppressed expression of osteogenic proteins such as Runt-related transcription factor 2 (RUNX2), Alkaline phosphatase (ALP), Osterix (OSX), and Osteopontin (OPN), and increased lipid droplets. Moreover, bioinformatics analysis substantiated shared gene expression patterns between pyroptosis and OP pathology, encompassing NLRP3, Caspase-1, IL-1β, IL-18, etc. Furthermore, our in vitro investigation using ST2 cells revealed that dexamethasone treatment prominently induced pyroptosis while impeding osteogenic differentiation. Notably, gene silencing of Caspase-1 effectively counteracted the inhibitory effects of dexamethasone on osteogenic differentiation, as manifested by increased ALP activity and enhanced expression of RUNX2, ALP, OSX, and OPN.

CONCLUSION

Our findings unequivocally underscore that inhibition of Caspase-1-mediated pyroptosis promotes osteogenic differentiation, providing a promising therapeutic target for managing OP.

Genetic Variants in WNT16 and PKD2L1 Locus Affect Heel Ultrasound Bone Stiffness: Analyses from the General Population and Patients Evaluated for Osteoporosis

Osteoporosis, a complex chronic disease with increasing prevalence, is characterised by reduced bone mineral density (BMD) and increased fracture risk. The high heritability of BMD suggests substantial impact of the individual genetic disposition on bone phenotypes and the development of osteoporosis. In the past years, genome-wide association studies (GWAS) identified hundreds of genetic variants associated with BMD or osteoporosis. Here, we analysed 1103 single nucleotide polymorphisms (SNPs), previously identified as associated with estimated BMD (eBMD) in the UK Biobank. We assessed whether these SNPs are related to heel stiffness index obtained by quantitative ultrasound in 5665 adult participants of the Study of Health in Pomerania (SHIP). We confirmed 45 significant associations after correction for multiple testing. Next, we analysed six selected SNPs in 631 patients evaluated for osteoporosis [rs2707518 (CPED1/WNT16), rs3779381 (WNT16), rs115242848 (LOC101927709/EN1), rs10239787 (JAZF1), rs603424 (PKD2L1) and rs6968704 (JAZF1)]. Differences in minor allele frequencies (MAF) of rs2707518 and rs3779381 between SHIP participants (higher MAF) and patients evaluated for osteoporosis (lower MAF) indicated a protective effect of the minor allele on bone integrity. In contrast, differences in MAF of rs603424 indicated a harmful effect. Co-localisation analyses indicated that the rs603424 effect may be mediated via stearoyl-CoA desaturase (SCD) expression, an enzyme highly expressed in adipose tissue with a crucial role in lipogenesis. Taken together, our results support the role of the WNT16 pathway in the regulation of bone properties and indicate a novel causal role of SCD expression in adipose tissue on bone integrity.

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.

Study of genes polymorphisms in RANK/RANKL/OPG and WNT signaling pathways and their associations with bone parameters in broiler chicken

Limb problems are one of the most common problems with fast-growing meat-type chickens. Various bone abnormalities, which can lead to limping, bone weakness, or even fractures, bring overall discomfort to birds and a loss of production. Genetic aspects are often associated with these side effects on bone stability and are also cited as the dominant cause. These points to a close negative relationship of genetic selection for rapid growth with traits involved in bone integrity. Due to the assumption of an additive genetic background, improvements through genetic tools can be used. Our study is focused on selected genes of important signaling pathways for bone metabolism. We tried to detect polymorphisms that would show associations with selected bone parameters in a total of 48 broilers. Those were fast-growing Ross 308 hybrids and slow-growing Hubbard M22BxJA87A hybrids. The TNFRSF11A and WISP1 genes were tested. A total of fourteen polymorphisms were found, three of them were synonymous and five in the intron. In the case of four polymorphisms found in exons of the TNFRSF11A gene (c.11G > T, c.31G > A, c.37C > G, c.514G > A), associations with the observed bone parameters (bone strength, bone dimensions and bone mass) were demonstrated. The genetic architecture of bone traits is not fully understood, therefore the present study and the knowledge gained can help to increase the potential in poultry breeding processes and thus reduce the death of individuals.

Molecular and Cellular Mechanisms of Osteoporosis

Osteoporosis is a widespread systemic disease characterized by a decrease in bone mass and an imbalance of the microarchitecture of bone tissue. Experimental and clinical studies devoted to investigating the main pathogenetic mechanisms of osteoporosis revealed the important role of estrogen deficiency, inflammation, oxidative stress, cellular senescence, and epigenetic factors in the development of bone resorption due to osteoclastogenesis, and decreased mineralization of bone tissue and bone formation due to reduced function of osteoblasts caused by apoptosis and age-depended differentiation of osteoblast precursors into adipocytes. The current review was conducted to describe the basic mechanisms of the development of osteoporosis at molecular and cellular levels and to elucidate the most promising therapeutic strategies of pathogenetic therapy of osteoporosis based on articles cited in PubMed up to September 2023.

Runt-related Transcription Factors and Gene Regulatory Mechanisms in Skeletal Development and Diseases

PURPOSE OF REVIEW

Runt-related transcription factors (RUNX) play critical roles in skeletal development, metabolism, and diseases. In mammals, three RUNX members, namely RUNX1, RUNX2, and RUNX3, play distinct and redundant roles, although RUNX2 is a dominant factor in skeletal development and several skeletal diseases. This review is to provide an overview of the current understanding of RUNX-mediated transcriptional regulation in different skeletal cell types.

RECENT FINDINGS

Advances in chromatin immunoprecipitation and next-generation sequencing (ChIP-seq) have revealed genome-wide RUNX-mediated gene regulatory mechanisms, including their association with cis-regulatory elements and putative target genes. Further studies with genome-wide analysis and biochemical assays have shed light on RUNX-mediated pioneering action and involvements of RUNX2 in lipid-lipid phase separation. Emerging multi-layered mechanisms of RUNX-mediated gene regulations help us better understanding of skeletal development and diseases, which also provides clues to think how genome-wide studies can help develop therapeutic strategies for skeletal diseases.

Proceedings of the 2023 Santa Fe Bone Symposium: Progress and Controversies in the Management of Patients with Skeletal Diseases

The Santa Fe Bone Symposium (SFBS) held its 23rd annual event on August 5-6, 2023, in Santa Fe, New Mexico, USA. Attendees participated in-person and remotely, representing many states and countries. The program included plenary presentations, panel discussions, satellite symposia, a Project ECHO workshop, and a session on healthcare policy and reimbursement for fracture liaison programs. A broad range of topics were addressed, including transitions of osteoporosis treatments over a lifetime; controversies in vitamin D; update on Official Positions of the International Society for Clinical Densitometry; spine surgery and bone health; clinical applications of bone turnover markers; basic bone biology for clinicians; premenopausal-, pregnancy-, and lactation-associated osteoporosis; cancer treatment induced bone loss in patients with breast cancer and prostate cancer; genetic testing for skeletal diseases; and an update on nutrition and bone health. There were also sessions on rare bone diseases, including managing patients with hypophosphatasia; treatment of X-linked hypophosphatemia; and assessment and treatment of patients with hypoparathyroidism. There were oral presentations of abstracts by endocrinology fellows selected from those who participated in the Santa Fe Fellows Workshop on Metabolic Bone Diseases, held the 2 days prior to the SFBS. These proceedings of the 2023 SFBS present the clinical highlights and insights generated from many formal and informal discussions in Santa Fe.

Ambient PM2.5 Exposure and Bone Homeostasis: Analysis of UK Biobank Data and Experimental Studies in Mice and in Vitro

BACKGROUND

Previous evidence has identified exposure to fine ambient particulate matter (PM 2.5 ) as a leading risk factor for adverse health outcomes. However, to date, only a few studies have examined the potential association between long-term exposure to PM 2.5 and bone homeostasis.

OBJECTIVE

We sought to examine the relationship between long-term PM 2.5 exposure and bone health and explore its potential mechanism.

METHODS

This research included both observational and experimental studies. First, based on human data from UK Biobank, linear regression was used to explore the associations between long-term exposure to PM 2.5 (i.e., annual average PM 2.5 concentration for 2010) and bone mineral density [BMD; i.e., heel BMD (n = 37,440 ) and femur neck and lumbar spine BMD (n = 29,766 )], which were measured during 2014-2020. For the experimental animal study, C57BL/6 male mice were assigned to ambient PM 2.5 or filtered air for 6 months via a whole-body exposure system. Micro-computed tomography analyses were applied to measure BMD and bone microstructures. Biomarkers for bone turnover and inflammation were examined with histological staining, immunohistochemistry staining, and enzyme-linked immunosorbent assay. We also performed tartrate-resistant acid phosphatase (TRAP) staining and bone resorption assay to determine the effect of PM 2.5 exposure on osteoclast activity in vitro. In addition, the potential downstream regulators were assessed by real-time polymerase chain reaction and western blot.

RESULTS

We observed that long-term exposure to PM 2.5 was significantly associated with lower BMD at different anatomical sites, according to the analysis of UK Biobank data. In experimental study, mice exposed long-term to PM 2.5 exhibited excessive osteoclastogenesis, dysregulated osteogenesis, higher tumor necrosis factor-alpha (TNF- α ) expression, and shorter femur length than control mice, but they demonstrated no significant differences in femur structure or BMD. In vitro, cells stimulated with conditional medium of PM 2.5 -stimulated macrophages had aberrant osteoclastogenesis and differences in the protein/mRNA expression of members of the TNF- α / Traf 6 / c -Fos pathway, which could be partially rescued by TNF- α inhibition.

DISCUSSION

Our prospective observational evidence suggested that long-term exposure to PM 2.5 is associated with lower BMD and further experimental results demonstrated exposure to PM 2.5 could disrupt bone homeostasis, which may be mediated by inflammation-induced osteoclastogenesis. https://doi.org/10.1289/EHP11646.

Single-Cell Integration of BMD GWAS Results Prioritize Candidate Genes Influencing Age-Related Bone Loss

The regulation of bone mineral density (BMD) is highly influenced by genetics and age. Although genome-wide association studies (GWAS) for BMD have uncovered many genes through their proximity to associated variants (variant nearest-neighbor [VNN] genes), the cell-specific mechanisms of each VNN gene remain unclear. This is primarily due to the inability to prioritize these genes by cell type and age-related expression. Using age-related transcriptomics, we found that the expression of many VNN genes was upregulated in the bone and marrow from aged mice. Candidate genes from GWAS were investigated using single-cell RNA-sequencing (scRNA-seq) datasets to enrich for cell-specific expression signatures. VNN candidate genes are highly enriched in osteo-lineage cells, osteocytes, hypertrophic chondrocytes, and Lepr+ mesenchymal stem cells. These data were used to generate a "blueprint" for Cre-loxp mouse line selection for functional validation of candidate genes and further investigation of their role in BMD maintenance throughout aging. In VNN-gene-enriched cells, Sparc, encoding the extracellular matrix (ECM) protein osteonectin, was robustly expressed. This, along with expression of numerous other ECM genes, indicates that many VNN genes likely have roles in ECM deposition by osteoblasts. Overall, we provide data supporting streamlined translation of GWAS candidate genes to potential novel therapeutic targets for the treatment of osteoporosis. © 2023 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Identification of ROCK1 as a novel biomarker for postmenopausal osteoporosis and pan-cancer analysis

BACKGROUND

Postmenopausal osteoporosis (PMOP) is a prevalent bone disorder with significant global impact. The elevated risk of osteoporotic fracture in elderly women poses a substantial burden on individuals and society. Unfortunately, the current lack of dependable diagnostic markers and precise therapeutic targets for PMOP remains a major challenge.

METHODS

PMOP-related datasets GSE7429, GSE56814, GSE56815, and GSE147287, were downloaded from the GEO database. The DEGs were identified by "limma" packages. WGCNA and Machine Learning were used to choose key module genes highly related to PMOP. GSEA, DO, GO, and KEGG enrichment analysis was performed on all DEGs and the selected key hub genes. The PPI network was constructed through the GeneMANIA database. ROC curves and AUC values validated the diagnostic values of the hub genes in both training and validation datasets. xCell immune infiltration and single-cell analysis identified the hub genes' function on immune reaction in PMOP. Pan-cancer analysis revealed the role of the hub genes in cancers.

RESULTS

A total of 1278 DEGs were identified between PMOP patients and the healthy controls. The purple module and cyan module were selected as the key modules and 112 common genes were selected after combining the DEGs and module genes. Five Machine Learning algorithms screened three hub genes (KCNJ2, HIPK1, and ROCK1), and a PPI network was constructed for the hub genes. ROC curves validate the diagnostic values of ROCK1 in both the training (AUC = 0.73) and validation datasets of PMOP (AUC = 0.81). GSEA was performed for the low-ROCK1 patients, and the top enriched field included protein binding and immune reaction. DCs and NKT cells were highly expressed in PMOP. Pan-cancer analysis showed a correlation between low ROCK1 expression and SKCM as well as renal tumors (KIRP, KICH, and KIRC).

CONCLUSIONS

ROCK1 was significantly associated with the pathogenesis and immune infiltration of PMOP, and influenced cancer development, progression, and prognosis, which provided a potential therapy target for PMOP and tumors. However, further laboratory and clinical evidence is required before the clinical application of ROCK1 as a therapeutic target.

Single-Cell Transcriptomics of Bone Marrow Stromal Cells in Diversity Outbred Mice: A Model for Population-Level scRNA-Seq Studies

Genome-wide association studies (GWASs) have advanced our understanding of the genetics of osteoporosis; however, the challenge has been converting associations to causal genes. Studies have utilized transcriptomics data to link disease-associated variants to genes, but few population transcriptomics data sets have been generated on bone at the single-cell level. To address this challenge, we profiled the transcriptomes of bone marrow-derived stromal cells (BMSCs) cultured under osteogenic conditions from five diversity outbred (DO) mice using single-cell RNA-seq (scRNA-seq). The goal of the study was to determine if BMSCs could serve as a model to generate cell type-specific transcriptomic profiles of mesenchymal lineage cells from large populations of mice to inform genetic studies. By enriching for mesenchymal lineage cells in vitro, coupled with pooling of multiple samples and downstream genotype deconvolution, we demonstrate the scalability of this model for population-level studies. We demonstrate that dissociation of BMSCs from a heavily mineralized matrix had little effect on viability or their transcriptomic signatures. Furthermore, we show that BMSCs cultured under osteogenic conditions are diverse and consist of cells with characteristics of mesenchymal progenitors, marrow adipogenic lineage precursors (MALPs), osteoblasts, osteocyte-like cells, and immune cells. Importantly, all cells were similar from a transcriptomic perspective to cells isolated in vivo. We employed scRNA-seq analytical tools to confirm the biological identity of profiled cell types. SCENIC was used to reconstruct gene regulatory networks (GRNs), and we observed that cell types show GRNs expected of osteogenic and pre-adipogenic lineage cells. Further, CELLECT analysis showed that osteoblasts, osteocyte-like cells, and MALPs captured a significant component of bone mineral density (BMD) heritability. Together, these data suggest that BMSCs cultured under osteogenic conditions coupled with scRNA-seq can be used as a scalable and biologically informative model to generate cell type-specific transcriptomic profiles of mesenchymal lineage cells in large populations. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Generic Platform for the Multiplexed Targeted Electrochemical Detection of Osteoporosis-Associated Single Nucleotide Polymorphisms Using Recombinase Polymerase Solid-Phase Primer Elongation and Ferrocene-Modified Nucleoside Triphosphates

Osteoporosis is a multifactorial disease influenced by genetic and environmental factors, which contributes to an increased risk of bone fracture, but early diagnosis of this disease cannot be achieved using current techniques. We describe a generic platform for the targeted electrochemical genotyping of SNPs identified by genome-wide association studies to be associated with a genetic predisposition to osteoporosis. The platform exploits isothermal solid-phase primer elongation with ferrocene-labeled nucleoside triphosphates. Thiolated reverse primers designed for each SNP were immobilized on individual gold electrodes of an array. These primers are designed to hybridize to the SNP site at their 3'OH terminal, and primer elongation occurs only where there is 100% complementarity, facilitating the identification and heterozygosity of each SNP under interrogation. The platform was applied to real blood samples, which were thermally lysed and directly used without the need for DNA extraction or purification. The results were validated using Taqman SNP genotyping assays and Sanger sequencing. The assay is complete in just 15 min with a total cost of 0.3€ per electrode. The platform is completely generic and has immense potential for deployment at the point of need in an automated device for targeted SNP genotyping with the only required end-user intervention being sample addition.

Anaphase-Promoting Complex Subunit 1 Associates with Bone Mineral Density in Human Osteoporotic Bone

Genome-wide association studies (GWAS) are one of the most common approaches to identify genetic loci that are associated with bone mineral density (BMD). Such novel genetic loci represent new potential targets for the prevention and treatment of fragility fractures. GWAS have identified hundreds of associations with BMD; however, only a few have been functionally evaluated. A locus significantly associated with femoral neck BMD at the genome-wide level is intronic SNP rs17040773 located in the intronic region of the anaphase-promoting complex subunit 1 (ANAPC1) gene (p = 1.5 × 10-9). Here, we functionally evaluate the role of ANAPC1 in bone remodelling by examining the expression of ANAPC1 in human bone and muscle tissues and during the osteogenic differentiation of human primary mesenchymal stem cells (MSCs). The expression of ANAPC1 was significantly decreased 2.3-fold in bone tissues and 6.2-fold in muscle tissue from osteoporotic patients as compared to the osteoarthritic and control tissues. Next, we show that the expression of ANAPC1 changes during the osteogenic differentiation process of human MSCs. Moreover, the silencing of ANAPC1 in human osteosarcoma (HOS) cells reduced RUNX2 expression, suggesting that ANAPC1 affects osteogenic differentiation through RUNX2. Altogether, our results indicate that ANAPC1 plays a role in bone physiology and in the development of osteoporosis.

Examining Mechanisms for Voltage-Sensitive Calcium Channel-Mediated Secretion Events in Bone Cells

In addition to their well-described functions in cell excitability, voltage-sensitive calcium channels (VSCCs) serve a critical role in calcium (Ca2+)-mediated secretion of pleiotropic paracrine and endocrine factors, including those produced in bone. Influx of Ca2+ through VSCCs activates intracellular signaling pathways to modulate a variety of cellular processes that include cell proliferation, differentiation, and bone adaptation in response to mechanical stimuli. Less well understood is the role of VSCCs in the control of bone and calcium homeostasis mediated through secreted factors. In this review, we discuss the various functions of VSCCs in skeletal cells as regulators of Ca2+ dynamics and detail how these channels might control the release of bioactive factors from bone cells. Because VSCCs are druggable, a better understanding of the multiple functions of these channels in the skeleton offers the opportunity for developing new therapies to enhance and maintain bone and to improve systemic health.

COVID-19 severity: does the genetic landscape of rare variants matter?

Rare variants affecting host defense against pathogens may be involved in COVID-19 severity, but most rare variants are not expected to have a major impact on the course of COVID-19. We hypothesized that the accumulation of weak effects of many rare functional variants throughout the exome may contribute to the overall risk in patients with severe disease. This assumption is consistent with the omnigenic model of the relationship between genetic and phenotypic variation in complex traits, according to which association signals tend to spread across most of the genome through gene regulatory networks from genes outside the major pathways to disease-related genes. We performed whole-exome sequencing and compared the burden of rare variants in 57 patients with severe and 29 patients with mild/moderate COVID-19. At the whole-exome level, we observed an excess of rare, predominantly high-impact (HI) variants in the group with severe COVID-19. Restriction to genes intolerant to HI or damaging missense variants increased enrichment for these classes of variants. Among various sets of genes, an increased signal of rare HI variants was demonstrated predominantly for primary immunodeficiency genes and the entire set of genes associated with immune diseases, as well as for genes associated with respiratory diseases. We advocate taking the ideas of the omnigenic model into account in COVID-19 studies.

Association between AXIN1 gene polymorphism (rs9921222) of WNT signaling pathway and susceptibility to osteoporosis in Egyptian patients: a case-control study

BACKGROUND

Osteoporosis (OP) is the most prevalent metabolic bone disease. Numerous genetic loci are strongly related to OP. AXIN1 is a significant gene that serves an important role in the WNT signaling pathway. The aim of this study was to explore the association between the AXIN1 genetic polymorphism (rs9921222) and OP susceptibility.

METHODS

A total of 101 subjects were enrolled in the study (50 patients with OP and 51 healthy individuals). Genomic DNA was extracted from whole blood using the QIAamp DNA Blood Mini Kit, and the AXIN1 gene polymorphism (rs9921222) was genotyped by TaqMan allelic discrimination assays. A logistic regression analysis was used to assess the association between genotypes and OP risk.

RESULTS

We found that AXIN1 rs9921222 had a significant association with the susceptibility of OP under the homozygote model (TT vs. CC: OR = 16.6, CI = 2.03-136.4, p = 0.009), (CT vs. CC: OR = 6.3, CI = 1.23-31.8, p = 0.027), recessive genetic model (TT vs.TC-CC: OR = 13.6, CI = 1.7-110.4, p = 0.015), and the dominant model (TT-TC vs. CC: OR = 9.7, CI = 2.6-36.3, p < 0.001). Allele T was significantly associated with OP risk (T vs. C: OR = 10.5, CI = 3.5-31.15, p = 0.001). There was a statistically significant difference between genotypes in mean platelet volume (p = 0.004), and platelet distribution width (p = 0.025). In addition, lumbar spine bone density, and femur neck bone density were significantly different between genotypes (p < 0.001).

CONCLUSION

AXIN1 rs9921222 was associated with OP susceptibility in the Egyptian population and should be considered a potential determinant risk for OP.

A Prediction Model for Osteoporosis Risk Using a Machine-Learning Approach and Its Validation in a Large Cohort

BACKGROUND

Osteoporosis develops in the elderly due to decreased bone mineral density (BMD), potentially increasing bone fracture risk. However, the BMD is not regularly measured in a clinical setting. This study aimed to develop a good prediction model for the osteoporosis risk using a machine learning (ML) approach in adults over 40 years in the Ansan/Anseong cohort and the association of predicted osteoporosis risk with a fracture in the Health Examinees (HEXA) cohort.

METHODS

The 109 demographic, anthropometric, biochemical, genetic, nutrient, and lifestyle variables of 8,842 participants were manually selected in an Ansan/Anseong cohort and included in the ML algorithm. The polygenic risk score (PRS) of osteoporosis was generated with a genome-wide association study and added for the genetic impact of osteoporosis. Osteoporosis was defined with < -2.5 T scores of the tibia or radius compared to people in their 20s-30s. They were divided randomly into the training (n = 7,074) and test (n = 1,768) sets-Pearson's correlation between the predicted osteoporosis risk and fracture in the HEXA cohort.

RESULTS

XGBoost, deep neural network, and random forest generated the prediction model with a high area under the curve (AUC, 0.86) of the receiver operating characteristic (ROC) with 10, 15, and 20 features; the prediction model by XGBoost had the highest AUC of ROC, high accuracy and k-fold values (> 0.85) in 15 features among seven ML approaches. The model included the genetic factor, genders, number of children and breastfed children, age, residence area, education, seasons to measure, height, smoking status, hormone replacement therapy, serum albumin, hip circumferences, vitamin B6 intake, and body weight. The prediction models for women alone were similar to those for both genders, with lower accuracy. When the prediction model was applied to the HEXA study, the correlation between the fracture incidence and predicted osteoporosis risk was significant but weak (r = 0.173, P < 0.001).

CONCLUSION

The prediction model for osteoporosis risk generated by XGBoost can be applied to estimate osteoporosis risk. The biomarkers can be considered for enhancing the prevention, detection, and early therapy of osteoporosis risk in Asians.

Discovery of potential biomarkers for osteoporosis diagnosis by individual omics and multi-omics technologies

INTRODUCTION

Global aging has made osteoporosis an increasingly serious public health problem. Osteoporotic fractures seriously affect the quality of life of patients and increase disability and mortality rates. Early diagnosis is important for timely intervention. The continuous development of individual- and multi-omics methods is helpful for the exploration and discovery of biomarkers for the diagnosis of osteoporosis.

AREAS COVERED

In this review, we first introduce the epidemiological status of osteoporosis and then describe the pathogenesis of osteoporosis. Furthermore, the latest progress in individual- and multi-omics technologies for exploring biomarkers for osteoporosis diagnosis is summarized. Moreover, we clarify the advantages and disadvantages of the application of osteoporosis biomarkers obtained using the omics method. Finally, we put forward valuable views on the future research direction of diagnostic biomarkers of osteoporosis.

EXPERT OPINION

Omics methods undoubtedly provide greatly contribute to the exploration of diagnostic biomarkers of osteoporosis; however, in the future, the clinical validity and clinical utility of the obtained potential biomarkers should be thoroughly examined. In addition, the improvement and optimization of the detection methods for different types of biomarkers and standardization of the detection process guarantee the reliability and accuracy of the detection results.

Ethnic disparities in fracture risk assessment using polygenic scores

Whether the PGS developed using data from European ancestry is predictive of fracture risk for minorities remains unclear. This study demonstrated that PGSs based on common BMD-related genetic variants discovered in the European ancestry cohort are predictive of fracture risk in people of Asian but not African ancestry.

PURPOSE

Large-scale genome-wide association studies (GWAS) on bone mineral density (BMD) have been conducted predominantly in European cohorts. Genetic models based on common variants associated with BMD have been evaluated using almost exclusively European data, which could potentially exacerbate health disparities due to different linkage disequilibrium among different ethnic groups.

METHODS

UK Biobank (UKB) is a large-scale population-based observational study starting in 2006 that recruited 502,617 individuals aged between 40 and 69 years with genotypic and phenotypic data available. Based on the summary statistics of two GWAS studies of femoral neck BMD and total body BMD, we derived four PGSs and assessed the association between each PGS and prevalent/incident fractures within each ethnic group separately using Multivariate logistic regressions and Cox proportional hazard models. All models were adjusted for age, sex, and the first four principal components.

RESULTS

We assessed four PGSs derived from European cohorts. Significant associations were observed between PGSs and fracture in European and Asian cohorts but not in the African cohort. Of all four PGSs, [Formula: see text] performed the best. A standard deviation decreases in [Formula: see text] were associated with an increased hazard ratio (HR) of 1.24 (1.22-1.27), 1.28 (0.83-1.99), and 1.34 (1.10-1.64) in European, African, and Asian ancestry, respectively. A low BMD-related PGS is associated with up to 2.35- and 4.31-fold increased fracture risk in European and Asian populations.

CONCLUSIONS

These results showed that PGSs based on common BMD-related genetic variants discovered in the European ancestry cohort are predictive of fracture risk in people of Asian but not African ancestry.

SP7: from Bone Development to Skeletal Disease

PURPOSE OF REVIEW

The purpose of this review is to summarize the different roles of the transcription factor SP7 in regulating bone formation and remodeling, discuss current studies in investigating the causal relationship between SP7 mutations and human skeletal disease, and highlight potential therapeutic treatments that targeting SP7 and the gene networks that it controls.

RECENT FINDINGS

Cell-type and stage-specific functions of SP7 have been identified during bone formation and remodeling. Normal bone development regulated by SP7 is strongly associated with human bone health. Dysfunction of SP7 results in common or rare skeletal diseases, including osteoporosis and osteogenesis imperfecta with different inheritance patterns. SP7-associated signaling pathways, SP7-dependent target genes, and epigenetic regulations of SP7 serve as new therapeutic targets in the treatment of skeletal disorders. This review addresses the importance of SP7-regulated bone development in studying bone health and skeletal disease. Recent advances in whole genome and exome sequencing, GWAS, multi-omics, and CRISPR-mediated activation and inhibition have provided the approaches to investigate the gene-regulatory networks controlled by SP7 in bone and the therapeutic targets to treat skeletal disease.

Copy Number Variation and Osteoporosis

PURPOSE OF REVIEW

The purpose of this review is to summarize recent findings on copy number variations and susceptibility to osteoporosis.

RECENT FINDINGS

Osteoporosis is highly influenced by genetic factors, including copy number variations (CNVs). The development and accessibility of whole genome sequencing methods has accelerated the study of CNVs and osteoporosis. Recent findings include mutations in novel genes and validation of previously known pathogenic CNVs in monogenic skeletal diseases. Identification of CNVs in genes previously associated with osteoporosis (e.g. RUNX2, COL1A2, and PLS3) has confirmed their importance in bone remodelling. This process has been associated also with the ETV1-DGKB, AGBL2, ATM, and GPR68 genes, identified by comparative genomic hybridisation microarray studies. Importantly, studies in patients with bone pathologies have associated bone disease with the long non-coding RNA LINC01260 and enhancer sequences residing in the HDAC9 gene. Further functional investigation of genetic loci harbouring CNVs associated with skeletal phenotypes will reveal their role as molecular drivers of osteoporosis.

Emerging role of mesenchymal stem/stromal cells (MSCs) and MSCs-derived exosomes in bone- and joint-associated musculoskeletal disorders: a new frontier

Exosomes are membranous vesicles with a 30 to 150 nm diameter secreted by mesenchymal stem/stromal cells (MSCs) and other cells, such as immune cells and cancer cells. Exosomes convey proteins, bioactive lipids, and genetic components to recipient cells, such as microRNAs (miRNAs). Consequently, they have been implicated in regulating intercellular communication mediators under physiological and pathological circumstances. Exosomes therapy as a cell-free approach bypasses many concerns regarding the therapeutic application of stem/stromal cells, including undesirable proliferation, heterogeneity, and immunogenic effects. Indeed, exosomes have become a promising strategy to treat human diseases, particularly bone- and joint-associated musculoskeletal disorders, because of their characteristics, such as potentiated stability in circulation, biocompatibility, low immunogenicity, and toxicity. In this light, a diversity of studies have indicated that inhibiting inflammation, inducing angiogenesis, provoking osteoblast and chondrocyte proliferation and migration, and negative regulation of matrix-degrading enzymes result in bone and cartilage recovery upon administration of MSCs-derived exosomes. Notwithstanding, insufficient quantity of isolated exosomes, lack of reliable potency test, and exosomes heterogeneity hurdle their application in clinics. Herein, we will deliver an outline respecting the advantages of MSCs-derived exosomes-based therapy in common bone- and joint-associated musculoskeletal disorders. Moreover, we will have a glimpse the underlying mechanism behind the MSCs-elicited therapeutic merits in these conditions.

Emerging RUNX2-Mediated Gene Regulatory Mechanisms Consisting of Multi-Layered Regulatory Networks in Skeletal Development

Skeletal development is tightly coordinated by chondrocytes and osteoblasts, which are derived from skeletal progenitors, and distinct cell-type gene regulatory programs underlie the specification and differentiation of cells. Runt-related transcription factor 2 (Runx2) is essential to chondrocyte hypertrophy and osteoblast differentiation. Genetic studies have revealed the biological functions of Runx2 and its involvement in skeletal genetic diseases. Meanwhile, molecular biology has provided a framework for our understanding of RUNX2-mediated transactivation at a limited number of cis-regulatory elements. Furthermore, studies using next-generation sequencing (NGS) have provided information on RUNX2-mediated gene regulation at the genome level and novel insights into the multiple layers of gene regulatory mechanisms, including the modes of action of RUNX2, chromatin accessibility, the concept of pioneer factors and phase separation, and three-dimensional chromatin organization. In this review, I summarize the emerging RUNX2-mediated regulatory mechanism from a multi-layer perspective and discuss future perspectives for applications in the treatment of skeletal diseases.

High Bone Mass Disorders: New Insights From Connecting the Clinic and the Bench

Monogenic high bone mass (HBM) disorders are characterized by an increased amount of bone in general, or at specific sites in the skeleton. Here, we describe 59 HBM disorders with 50 known disease-causing genes from the literature, and we provide an overview of the signaling pathways and mechanisms involved in the pathogenesis of these disorders. Based on this, we classify the known HBM genes into HBM (sub)groups according to uniform Gene Ontology (GO) terminology. This classification system may aid in hypothesis generation, for both wet lab experimental design and clinical genetic screening strategies. We discuss how functional genomics can shape discovery of novel HBM genes and/or mechanisms in the future, through implementation of omics assessments in existing and future model systems. Finally, we address strategies to improve gene identification in unsolved HBM cases and highlight the importance for cross-laboratory collaborations encompassing multidisciplinary efforts to transfer knowledge generated at the bench to the clinic. © 2022 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Nutritional and lifestyle management of the aging journey: A narrative review

With age, the physiological responses to occasional or regular stressors from a broad range of functions tend to change and adjust at a different pace and restoring these functions in the normal healthy range becomes increasingly challenging. Even if this natural decline is somehow unavoidable, opportunities exist to slow down and attenuate the impact of advancing age on major physiological processes which, when weakened, constitute the hallmarks of aging. This narrative review revisits the current knowledge related to the aging process and its impact on key metabolic functions including immune, digestive, nervous, musculoskeletal, and cardiovascular functions; and revisits insights into the important biological targets that could inspire effective strategies to promote healthy aging.

Twelve Loci Associated With Bone Density in Middle-aged and Elderly Chinese: The Shanghai Changfeng Study

CONTEXT

Previous genome-wide association studies (GWASs) of bone mineral density (BMD) were mainly conducted in Europeans.

OBJECTIVE

To explore genetic variants that affect BMD and sex differences in a Chinese population.

METHODS

A total of 5428 middle-aged and elderly Chinese were included. Dual-energy X-ray absorptiometry was used to measure BMD at the lumbar spine, and total and specific sites of the hip. A mixed linear model was used to analyze the associations between BMD and autosomal genetic variants, adjusting for age, age squared, sex, and menopausal women (model 1); model 2 was further adjusted for height and weight. A GWAS of osteoporosis in the Biobank Japan (BBJ) project was used for replication. GWAMA software was used to detect the statistical significance of sex differences of estimated effects. Gene annotation and pathway enrichment analysis were performed.

RESULTS

Women lost BMD at earlier ages and faster than men. The 2 models identified a total of 12 loci that were associated with BMD at any site. Single nucleotide polymorphisms rs72354346, rs2024219, rs1463093, rs10037512, and rs5880932 were successfully replicated in the BBJ. Variations of rs79262027 G>A (VKORC1L1) and rs4795209 A>G (DDX52) were associated with BMD only in men, and rs1239055408 G>GA (KCNJ2) was associated with BMD only in women. Gene enrichment analysis showed that BMD in a Chinese elderly population was mainly related to ossification, bone resorption, sex hormones, and kidney physiology.

CONCLUSION

The present GWAS identified 12 loci that were significantly associated with BMD at any site in a Chinese population, and 3 of them showed sex differences in their effects.

Factors associated with bone thickness: Comparison of the cranium and humerus

Cortical bone thickness is important for the mechanical function of bone. Ontogeny, aging, sex, body size, hormone levels, diet, behavior, and genetics potentially cause variations in postcranial cortical robusticity. However, the factors associated with cranial cortical robusticity remain poorly understood. Few studies have examined cortical robusticity in both cranial and postcranial bones jointly. In the present study, we used computed tomography (CT) images to measure cortical bone thicknesses in the cranial vault and humeral diaphysis. This study clearly showed that females have a greater cranial vault thickness and greater age-related increase in cranial vault thickness than males. We found an age-related increase in the full thickness of the temporal cranial vault and the width of the humeral diaphysis, as well as an age-related decrease in the cortical thickness of the frontal cranial vault and the cortical thickness of the humeral diaphysis, suggesting that the mechanisms of bone modeling in cranial and long bones are similar. A positive correlation between cortical indices in the cranial vault and humeral diaphysis also suggested that common factors affect cortical robusticity. We also examined the association of polymorphisms in the WNT16 and TNFSF11 genes with bone thickness. However, no significant associations were observed. The present study provides fundamental knowledge about similarities and differences in the mechanisms of bone modeling between cranial and postcranial bones.