Bone mineral density, osteoporosis, and osteoporotic fractures: a genome-wide association study

Genome-wide study identifies novel genes associated with bone toxicities in children with acute lymphoblastic leukaemia

Bone toxicities are common among paediatric patients treated for acute lymphoblastic leukaemia (ALL) with potentially major negative impact on patients' quality of life. To identify the underlying genetic contributors, we conducted a genome-wide association study (GWAS) and a transcriptome-wide association study (TWAS) in 260 patients of European-descent from the DFCI 05-001 ALL trial, with validation in 101 patients of European-descent from the DFCI 11-001 ALL trial. We identified a significant association between rs844882 on chromosome 20 and bone toxicities in the DFCI 05-001 trial (p = 1.7 × 10-8). In DFCI 11-001 trial, we observed a consistent trend of this variant with fracture. The variant was an eQTL for two nearby genes, CD93 and THBD. In TWAS, genetically predicted ACAD9 expression was associated with an increased risk of bone toxicities, which was confirmed by meta-analysis of the two cohorts (meta-p = 2.4 × 10-6). In addition, a polygenic risk score of heel quantitative ultrasound speed of sound was associated with fracture risk in both cohorts (meta-p = 2.3 × 10-3). Our findings highlight the genetic influence on treatment-related bone toxicities in this patient population. The genes we identified in our study provide new biological insights into the development of bone adverse events related to ALL treatment.

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.

Osteoporotic osseointegration: therapeutic hallmarks and engineering strategies

Osteoporosis is a systemic skeletal disease caused by an imbalance between bone resorption and formation. Current treatments primarily involve systemic medication and hormone therapy. However, these systemic treatments lack directionality and are often ineffective for locally severe osteoporosis, with the potential for complex adverse reactions. Consequently, treatment strategies using bioactive materials or external interventions have emerged as the most promising approaches. This review proposes twelve microenvironmental treatment targets for osteoporosis-related pathological changes, including local accumulation of inflammatory factors and reactive oxygen species (ROS), imbalance of mitochondrial dynamics, insulin resistance, disruption of bone cell autophagy, imbalance of bone cell apoptosis, changes in neural secretions, aging of bone cells, increased local bone tissue vascular destruction, and decreased regeneration. Additionally, this review examines the current research status of effective or potential biophysical and biochemical stimuli based on these microenvironmental treatment targets and summarizes the advantages and optimal parameters of different bioengineering stimuli to support preclinical and clinical research on osteoporosis treatment and bone regeneration. Finally, the review addresses ongoing challenges and future research prospects.

Signalling interaction between β-catenin and other signalling molecules during osteoarthritis development

Osteoarthritis (OA) is the most prevalent disorder of synovial joint affecting multiple joints. In the past decade, we have witnessed conceptual switch of OA pathogenesis from a 'wear and tear' disease to a disease affecting entire joint. Extensive studies have been conducted to understand the underlying mechanisms of OA using genetic mouse models and ex vivo joint tissues derived from individuals with OA. These studies revealed that multiple signalling pathways are involved in OA development, including the canonical Wnt/β-catenin signalling and its interaction with other signalling pathways, such as transforming growth factor β (TGF-β), bone morphogenic protein (BMP), Indian Hedgehog (Ihh), nuclear factor κB (NF-κB), fibroblast growth factor (FGF), and Notch. The identification of signalling interaction and underlying mechanisms are currently underway and the specific molecule(s) and key signalling pathway(s) playing a decisive role in OA development need to be evaluated. This review will focus on recent progresses in understanding of the critical role of Wnt/β-catenin signalling in OA pathogenesis and interaction of β-catenin with other pathways, such as TGF-β, BMP, Notch, Ihh, NF-κB, and FGF. Understanding of these novel insights into the interaction of β-catenin with other pathways and its integration into a complex gene regulatory network during OA development will help us identify the key signalling pathway of OA pathogenesis leading to the discovery of novel therapeutic strategies for OA intervention.

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.

The respective and dependent effects of scattering and bone matrix absorption on ultrasound attenuation in cortical bone

Cortical bone is characterized by a dense solid matrix permeated by fluid-filled pores. Ultrasound scattering has potential for the non-invasive evaluation of changes in bone porosity. However, there is an incomplete understanding of the impact of ultrasonic absorption in the solid matrix on ultrasound scattering. In this study, maps were derived from scanning acoustic microscopy images of human femur cross-sections. Finite-difference time domain ultrasound scatter simulations were conducted on these maps. Pore density, diameter distribution of the pores, and nominal absorption values in the solid and fluid matrices were controlled. Ultrasound pulses with a central frequency of 8.2 MHz were propagated, both in through-transmission and backscattering configurations. From these data, the scattering, bone matrix absorption, and attenuation extinction lengths were calculated. The results demonstrated that as absorption in the solid matrix was varied, the scattering, absorption, and attenuation extinction lengths were significantly impacted. It was shown that for lower values of absorption in the solid matrix (less than 2 dB mm-1), attenuation due to scattering dominates, whereas at higher values of absorption (more than 2 dB mm-1), attenuation due to absorption dominates. This will impact how ultrasound attenuation and scattering parameters can be used to extract quantitative information on bone microstructure.

Associations of residential greenness with bone mineral density and osteoporosis: the modifying effect of genetic susceptibility

OBJECTIVES

To investigate the associations of residential greenness with bone mineral density and incident osteoporosis, and further evaluate the potential modifying effect of genetic susceptibility.

METHODS

We used the Normalised Difference Vegetation Index (NDVI) at various buffer distances, including 300 m (NDVI300m), 500 m (NDVI500m), 1000 m (NDVI1000m) and 1500 m (NDVI1500m), to serve as indicators of greenness. We fitted linear regression, logistic regression and Cox proportional hazard models to assess the associations of residential greenness with estimated bone mineral density (eBMD), prevalent osteoporosis and incident osteoporosis, respectively. With the Polygenic Risk Score (PRS) for osteoporosis, we further assessed the joint effects of genetic risk and greenness on the risk of osteoporosis. We conducted causal mediation analyses to explore potential mediators.

RESULTS

Each IQR increase in NDVI300m was associated with 0.0007 (95% CI 0.0002 to 0.0013) increase in eBMD, 6% lower risk of prevalent osteoporosis (OR 0.94; 95% CI 0.92 to 0.97) and 5% lower risk of incident osteoporosis (HR 0.95; 95% CI 0.93 to 0.98). The joint effects of greenness and PRS on the risk of osteoporosis displayed a clear dose-response pattern. Compared with individuals exposed to low NDVI levels and high genetic risk, those exposed to high NDVI levels and low genetic risk had a 56% (95% CI 51% to 61%) lower risk of osteoporosis. The primary mediators in the association between greenness and incident osteoporosis were identified as PM2.5 and NO2.

CONCLUSIONS

Residential greenness was associated with higher bone mineral density and decreased risk of incident osteoporosis.

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.

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.

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.

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.

Associations between epigenome-wide DNA methylation and height-related traits among Sub-Saharan Africans: the RODAM study

Human height and related traits are highly complex, and extensively research has shown that these traits are determined by both genetic and environmental factors. Such factors may partially affect these traits through epigenetic programing. Epigenetic programing is dynamic and plays an important role in controlling gene expression and cell differentiation during (early) development. DNA methylation (DNAm) is the most commonly studied epigenetic feature. In this study we conducted an epigenome-wide DNAm association analysis on height-related traits in a Sub-Saharan African population, in order to detect DNAm biomarkers across four height-related traits. DNAm profiles were acquired in whole blood samples of 704 Ghanaians, sourced from the Research on Obesity and Diabetes among African Migrants study, using the Illumina Infinium HumanMethylation450 BeadChip. Linear models were fitted to detect differentially methylated positions (DMPs) and regions (DMRs) associated with height, leg-to-height ratio (LHR), leg length, and sitting height. No epigenome-wide significant DMPs were recorded. However we did observe among our top DMPs five informative probes associated with the height-related traits: cg26905768 (leg length), cg13268132 (leg length), cg19776793 (height), cg23072383 (LHR), and cg24625894 (sitting height). All five DMPs are annotated to genes whose functions were linked to bone cell regulation and development. DMR analysis identified overlapping DMRs within the gene body of HLA-DPB1 gene, and the HOXA gene cluster. In this first epigenome-wide association studies of these traits, our findings suggest DNAm associations with height-related heights, and might influence development and maintenance of these traits. Further studies are needed to replicate our findings, and to elucidate the molecular mechanism underlying human height-related traits.

Bone Trans-omics: Integrating Omics to Unveil Mechanistic Molecular Networks Regulating Bone Biology and Disease

PURPOSE OF REVIEW

Recent advancements in "omics" technologies and bioinformatics have afforded researchers new tools to study bone biology in an unbiased and holistic way. The purpose of this review is to highlight recent studies integrating multi-omics data gathered from multiple molecular layers (i.e.; trans-omics) to reveal new molecular mechanisms that regulate bone biology and underpin skeletal diseases.

RECENT FINDINGS

Bone biologists have traditionally relied on single-omics technologies (genomics, transcriptomics, proteomics, and metabolomics) to profile measureable differences (both qualitative and quantitative) of individual molecular layers for biological discovery and to investigate mechanisms of disease. Recently, literature has grown on the implementation of integrative multi-omics to study bone biology, which combines computational and informatics support to connect multiple layers of data derived from individual "omic" platforms. This emerging discipline termed "trans-omics" has enabled bone biologists to identify and construct detailed molecular networks, unveiling new pathways and unexpected interactions that have advanced our mechanistic understanding of bone biology and disease. While the era of trans-omics is poised to revolutionize our capacity to answer more complex and diverse questions pertinent to bone pathobiology, it also brings new challenges that are inherent when trying to connect "Big Data" sets. A concerted effort between bone biologists and interdisciplinary scientists will undoubtedly be needed to extract physiologically and clinically meaningful data from bone trans-omics in order to advance its implementation in the field.

Osteoprotegerin genetic polymorphisms and their influence on therapeutic response to ibandronate in postmenopausal osteoporotic females

OBJECTIVES

The present study investigated osteoprotegerin (OPG) genetic polymorphisms and their influence on the therapeutic response to ibandronate in postmenopausal osteoporotic females.

METHODS

This case-control study included 135 postmenopausal females (89 osteoporotic females and 46 non-osteoporotic females). Each osteoporotic patient received a monthly 150 mg ibandronate tablet for six months, and blood samples were taken before and after treatment. Bone mineral density (BMD) was measured using DEXA Scan. Three SNPs (A163G, T245G, and G1181C) of the OPG gene were selected for analysis.

RESULTS

Serum OPG levels were significantly lower in osteoporotic subjects than in the control group. The percentage changes in OPG levels in the osteoporotic group before and after treatment with ibandronate were significant (p < .001). After six months of therapy with ibandronate, the percentage changes in OPG levels with AA, TT, TC, GC, and GG genotypes were significant. Following six months of ibandronate treatment, the AA genotype of rs3134069, TT, TC genotypes of rs3102735, GG, and GC genotypes of rs2073618 SNP showed a significant increase in OPG levels. Age, BMI, and GC polymorphism (rs2073618 (G/C) G1181C) were inversely associated with low BMD. Adjusted odds ratios (OR) showed that BMI, GC, GG polymorphism (rs2073618 (G/C) G1181C) and TC polymorphism (rs3102735 (T/C) A163G) were inversely associated with low BMD.

CONCLUSION

The inverse association of rs2073618 and rs3102735 with low BMD indicates the protective role of these SNPs in our population. More research is needed to replicate these results in another cohort and to determine the molecular processes by which such SNPs may influence BMD.

Hedgehog signaling is a potential therapeutic target for vascular calcification

BACKGROUND

Patients with chronic kidney disease (CKD) suffered from vascular calcification (VC), one major contributor for their increased mortality rate. Hedgehog (Hh) signaling plays a crucial role in physiological bone mineralization and is associated with several cardiovascular diseases. However, the molecular changes underlying VC is ill defined and it remains unclear whether Hh signaling intervention affects VC.

METHODS

We constructed human primary vascular smooth muscle cell (VSMC) calcification model and performed RNA sequencing. Alizarin red staining and calcium content assay were conducted to identify the occurrence of VC. Three different R packages were applied to determine differentially expressed genes (DEGs). Enrichment analysis and protein-protein interaction (PPI) network analysis were carried out to explore the biological roles of DEGs. qRT-PCR assay was then applied to validate the expression of key genes. By using Connectivity Map (CMAP) analysis, several small molecular drugs targeting these key genes were obtained, including SAG (Hedgehog signaling activator) and cyclopamine (CPN) (Hedgehog signaling inhibitor), which were subsequently used to treat VSMC.

RESULTS

Obvious Alizarin red staining and increased calcium content identified the occurrence of VC. By integrating results from three R packages, we totally obtained 166 DEGs (86 up-regulated and 80 down-regulated), which were significantly enriched in ossification, osteoblast differentiation, and Hh signaling. PPI network analysis identified 10 key genes and CMAP analysis predicted several small molecular drugs targeting these key genes including chlorphenamine, isoeugenol, CPN and phenazopyridine. Notably, our in vitro experiment showed that SAG markedly alleviated VSMC calcification, whereas CPN significantly exacerbated VC.

CONCLUSIONS

Our research provided deeper insight to the pathogenesis of VC and indicated that targeting Hh signaling pathway may represent a potential and effective therapy for VC.

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.

Biomineralized Piezoelectrically Active Scaffolds for Inducing Osteogenic Differentiation

There is an endogenous electric field in living organisms, which plays a vital role in the development and regeneration of bone tissue. Therefore, self-powered piezoelectric material for bone repair has become hot research in recent years. However, the current piezoelectric materials for tissue regeneration still have the shortcomings of lack of biological activity and three-dimensional structure. Here, we proposed a three-dimensional polyurethane foam (PUF) scaffold coated with piezoelectric poly (vinylidene fluoride-co-hexafluoropropylene) (PVDF-HFP) and modified by a calcium phosphate (CaP) mineralized coating. The preferred scaffold has an open circuit voltage and short circuit current output of 5 V and 200 nA. Combining the physical and chemical properties of the CaP coating, the piezoelectric signal of PVDF-HFP and the three-dimensional structure of PUF, the scaffold exhibits superior promotion of cell osteogenic differentiation and ectopic bone formation in vivo. The mechanism is attributed to an increase in intracellular Ca²⁺ levels in response to chemical and piezoelectric stimulation with the material. This research not only paves the way for the application of piezoelectric scaffolds to stimulate osteoblasts differentiation in situ, but also lays the foundation for the clinical treatment of long-term osteoporosis.

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Frailty and genetic risk predict fracture after lung transplantation

Fractures negatively impact quality of life and survival. We hypothesized that recipient frailty score and genetic profile measured before transplant would predict risk of fracture after lung transplant. We conducted a retrospective cohort study of bone mineral density (BMD) and fracture among lung transplant recipients at a single center. The association between predictors and outcomes were assessed by multivariable time-dependent Cox models or regression analysis. Among the 284 participants, osteoporosis and fracture were highly prevalent. Approximately 59% of participants had posttransplant osteopenia, and 35% of participants developed at least 1 fracture. Low BMD was associated with a polygenic osteoporosis risk score, and the interaction between genetic score and BMD predicted fracture. Pretransplant frailty was associated with risk for spine and hip fracture, which were not associated with chronic lung allograft dysfunction or death. Chest fractures were the most frequent type of fracture and conferred a 2.2-fold increased risk of chronic lung allograft dysfunction or death (time-dependent P < .001). Pneumonia, pleural effusions, and acute rejection frequently occurred surrounding chest fracture. Pretransplant frailty and recipient genotype may aid clinical risk stratification for fracture after transplant. Fracture carries significant morbidity, underscoring the importance of surveillance and osteoporosis prevention.

Aging and Bone Metabolism

Changes in bone architecture and metabolism with aging increase the likelihood of osteoporosis and fracture. Age-onset osteoporosis is multifactorial, with contributory extrinsic and intrinsic factors including certain medical problems, specific prescription drugs, estrogen loss, secondary hyperparathyroidism, microenvironmental and cellular alterations in bone tissue, and mechanical unloading or immobilization. At the histological level, there are changes in trabecular and cortical bone as well as marrow cellularity, lineage switching of mesenchymal stem cells to an adipogenic fate, inadequate transduction of signals during skeletal loading, and predisposition toward senescent cell accumulation with production of a senescence-associated secretory phenotype. Cumulatively, these changes result in bone remodeling abnormalities that over time cause net bone loss typically seen in older adults. Age-related osteoporosis is a geriatric syndrome due to the multiple etiologies that converge upon the skeleton to produce the ultimate phenotypic changes that manifest as bone fragility. Bone tissue is dynamic but with tendencies toward poor osteoblastic bone formation and relative osteoclastic bone resorption with aging. Interactions with other aging physiologic systems, such as muscle, may also confer detrimental effects on the aging skeleton. Conversely, individuals who maintain their BMD experience a lower risk of fractures, disability, and mortality, suggesting that this phenotype may be a marker of successful aging. © 2023 American Physiological Society. Compr Physiol 13:4355-4386, 2023.

Integrating Common Risk Factors with Polygenic Scores Improves the Prediction of Type 2 Diabetes

We tested associations between 13 established genetic variants and type 2 diabetes (T2D) in 1371 study participants from the Volga-Ural region of the Eurasian continent, and evaluated the predictive ability of the model containing polygenic scores for the variants associated with T2D in our dataset, alone and in combination with other risk factors such as age and sex. Using logistic regression analysis, we found associations with T2D for the CCL20 rs6749704 (OR = 1.68, P_FDR = 3.40 × 10^-5), CCR5 rs333 (OR = 1.99, P_FDR = 0.033), ADIPOQ rs17366743 (OR = 3.17, P_FDR = 2.64 × 10^-4), TCF7L2 rs114758349 (OR = 1.77, P_FDR = 9.37 × 10^-5), and CCL2 rs1024611 (OR = 1.38, P_FDR = 0.033) polymorphisms. We showed that the most informative prognostic model included weighted polygenic scores for these five loci, and non-genetic factors such as age and sex (AUC 85.8%, 95%CI 83.7-87.8%). Compared to the model containing only non-genetic parameters, adding the polygenic score for the five T2D-associated loci showed improved net reclassification (NRI = 37.62%, 1.39 × 10^-6). Inclusion of all 13 tested SNPs to the model with age and sex did not improve the predictive ability compared to the model containing five T2D-associated variants (NRI = -17.86, p = 0.093). The five variants associated with T2D in people from the Volga-Ural region are linked to inflammation (CCR5, CCL2, CCL20) and glucose metabolism regulation (TCF7L, ADIPOQ2). Further studies in independent groups of T2D patients should validate the prognostic value of the model and elucidate the molecular mechanisms of the disease development.

CircZNF367 suppresses osteogenic differentiation of human bone marrow mesenchymal stromal/stem cells via reducing HuR-mediated mRNA stability of LRP5

Osteoporosis is a highly prevalent disease characterized by bone mass loss and structural deterioration. There are evidences that altered differentiation of human bone marrow mesenchymal stromal/stem cells (hBMSCs) is a major cause for osteoporosis. Recent studies suggest that circular RNAs (circRNAs) are dysregulated in osteoporosis patients and involved in the pathogenesis of osteoporosis. In the present study, we are aimed to analyze the circRNA expression profiles in osteoporosis patients and identify potential circRNAs that involved in the differentiation of hBMSCs during osteoporosis. Transcriptome RNA-sequencing was conducted to search for differentially expressed circRNAs. Transwell assay, ARS and ALP staining, and ectopic bone formation model were performed to evaluate osteogenic differentiation of hBMSCs. RNA pull-down assay, RNA immunoprecipitation, western blot, and in vitro binding assay were conducted to evaluate the interaction of circRNAs and RNA-binding protein HuR. We found that hsa_circ_0008842 (designated as circZNF367) was upregulated in osteoporosis patients and decreased in hBMSCs during osteogenic differentiation. CircZNF367 overexpression suppressed migration, invasion and osteogenic differentiation of hBMSCs in vitro and in vivo. In comparison, knockdown of circZNF367 promoted migration, invasion and osteogenic differentiation of hBMSCs. CircZNF367 could interact with the RNA-binding protein HuR, thus reduced the mRNA stability of LRP5. Furthermore, HuR overexpression or LRP5 restoration abrogated the effects of circZNF367 overexpression on osteogenic differentiation of hBMSCs. Our results indicated that circZNF367 played a role in osteogenic differentiation of hBMSCs via reducing HuR-mediated mRNA stability of LRP5.

The association between depression and bone metabolism: a US nationally representative cross-sectional study

This population-based study investigated the association between depression and bone mineral density (BMD), fractures, and osteoporosis in the US population. We found that participants with depression had lower BMD and were more likely to have fractures and osteoporosis.

BACKGROUND

Depression, fractures, and osteoporosis are common in middle-aged and elderly, but their associations remained unclear.

OBJECTIVE

To investigate the association between depression and bone mineral density (BMD), osteoporosis, and fracture in a middle-aged and elderly US population.

METHODS

A nationally representative cross-sectional study used the National Health and Nutrition Examination Survey (NHANES) datasets. Depression was assessed and stratified using the Patient Health Questionnaire (PHQ-9). The multiple logistic regression models and the logistic binary regression models were used to analyze the association between depression and BMD, fractures, and osteoporosis. Gender, age, race, educational level, poverty ratio, body mass index (BMI), smoke, alcohol use, physical activity, and diabetes were included as covariates. Subgroup analysis was also conducted on gender, age, race, and education level.

RESULTS

In total, 9766 participants were included after a series of exclusions, and 4179 (42.79%) had at least mild depressive symptoms. Compared to the participants without depression, those with depression had a lower total femur, femoral neck, and total spine BMD after adjusting multiple covariates. The multivariable-adjusted logistic binary regression models demonstrated that participants with depression more likely have hip fractures (OR = 1.518, 95% CI: 1.377-2.703, P = 0.000), spine fractures (OR = 1.311, 95% CI: 1.022-1.678, P = 0.030), and osteoporosis (OR = 1.621, 95% CI: 1.388-1.890, P = 0.000). Subgroup analysis revealed that depressed participants who were males, non-Hispanic White, ≤ 70 years, and not highly educated had a lower BMD and easily had osteoporosis.

CONCLUSION

Depression was associated with lower BMD, particularly in the spine, males, Hispanic-White, and not highly educated populations. Moreover, people with depression were more likely to have fractures and osteoporosis.

Genetic Diagnostics in Routine Osteological Assessment of Adult Low Bone Mass Disorders

CONTEXT

Many different inherited and acquired conditions can result in premature bone fragility/low bone mass disorders (LBMDs).

OBJECTIVE

We aimed to elucidate the impact of genetic testing on differential diagnosis of adult LBMDs and at defining clinical criteria for predicting monogenic forms.

METHODS

Four clinical centers broadly recruited a cohort of 394 unrelated adult women before menopause and men younger than 55 years with a bone mineral density (BMD) Z-score < -2.0 and/or pathological fractures. After exclusion of secondary causes or unequivocal clinical/biochemical hallmarks of monogenic LBMDs, all participants were genotyped by targeted next-generation sequencing.

RESULTS

In total, 20.8% of the participants carried rare disease-causing variants (DCVs) in genes known to cause osteogenesis imperfecta (COL1A1, COL1A2), hypophosphatasia (ALPL), and early-onset osteoporosis (LRP5, PLS3, and WNT1). In addition, we identified rare DCVs in ENPP1, LMNA, NOTCH2, and ZNF469. Three individuals had autosomal recessive, 75 autosomal dominant, and 4 X-linked disorders. A total of 9.7% of the participants harbored variants of unknown significance. A regression analysis revealed that the likelihood of detecting a DCV correlated with a positive family history of osteoporosis, peripheral fractures (> 2), and a high normal body mass index (BMI). In contrast, mutation frequencies did not correlate with age, prevalent vertebral fractures, BMD, or biochemical parameters. In individuals without monogenic disease-causing rare variants, common variants predisposing for low BMD (eg, in LRP5) were overrepresented.

CONCLUSION

The overlapping spectra of monogenic adult LBMD can be easily disentangled by genetic testing and the proposed clinical criteria can help to maximize the diagnostic yield.

Photoacoustic characterization of bone physico-chemical information

Osteoporosis usually alters the chemical composition and physical microstructure of bone. Currently, most clinical techniques for bone assessment are focused on the either bone microstructure or bone mineral density (BMD). In this study, a novel multi-wavelength photoacoustic time-frequency spectral analysis (MWPA-TFSA) method was introduced based on the optical absorption spectra and photoacoustic effects of biological macromolecules, which evaluates changes in bone chemical composition and microstructure. The results demonstrated that osteoporotic bones had decreased BMD, more lipids, and wider trabecular separation filled with larger marrow clusters, which were consistent with multiple gold-standard results, suggesting that the MWPA-TFSA method has the potential to provide a thorough bone physico-chemical information evaluation noninvasively and nonradiatively.

The surgical menopause

Surgical menopause (iatrogenic menopause) happens when both ovaries are removed before the natural "switching off" of ovarian function; it can cause premature ovarian insufficiency where the menopause occurs in women before the age of 40. Surgical menopause is associated with a sudden reduction of ovarian sex steroid production rather than a gradual one as is the case in natural menopause. In women who have undergone bilateral salpingo-oophorectomy (BSO) before the natural age of menopause, strong consideration should be given to giving hormone replacement therapy (HRT) till the natural age of menopause at least. Sexual function and sexual desire are altered post-BSO, especially in younger women hence part of HRT prescription must include consideration of androgen too.

Prioritization of Genes Relevant to Bone Fragility Through the Unbiased Integration of Aging Mouse Bone Transcriptomics and Human GWAS Analyses

Identifying new genetic determinants of bone mineral density (BMD) and fracture promises to yield improved diagnostics and therapies for bone fragility. However, prioritizing candidate genes from genome-wide screens can be challenging. To overcome this challenge, we prioritized mouse genes that are differentially expressed in aging mouse bone based on whether their human homolog is associated with human BMD and/or fracture. Unbiased RNA-seq analysis of young and old male C57BL/6 mouse cortical bone identified 1499, 1685, and 5525 differentially expressed genes (DEGs) in 1, 2, and 2.5-year-old bone, relative to 2-month-old bone, respectively. Gene-based scores for heel ultrasound bone mineral density (eBMD) and fracture were estimated using published genome-wide association studies (GWAS) results of these traits in the UK Biobank. Enrichment analysis showed that mouse bone DEG sets for all three age groups, relative to young bone, are significantly enriched for eBMD, but only the oldest two DEG sets are enriched for fracture. Using gene-based scores, this approach prioritizes among thousands of DEGs by a factor of 5- to 100-fold, yielding 10 and 21 genes significantly associated with fracture in the two oldest groups of mouse DEGs. Though these genes were not the most differentially expressed, they included Sost, Lrp5, and others with well-established functions in bone. Several others have, as yet, unknown roles in the skeleton. Therefore, this study accelerates identification of new genetic determinants of bone fragility by prioritizing a clinically relevant and experimentally tractable number of candidate genes for functional analysis. Finally, we provide a website (www.mouse2human.org) to enable other researchers to easily apply our strategy. © 2022 American Society for Bone and Mineral Research (ASBMR).

Advances in pathogenesis and therapeutic strategies for osteoporosis

Osteoporosis, is the most common bone disorder worldwide characterized by low bone mineral density, leaving affected bones vulnerable to fracture. Bone homeostasis depends on the precise balance between bone resorption by osteoclasts and bone matrix formation by mesenchymal lineage osteoblasts, and involves a series of complex and highly regulated steps. Bone homeostasis will be disrupted when the speed of bone resorption is faster than bone formation. Based on various regulatory mechanisms of bone homeostasis, a series of drugs targeting osteoporosis have emerged in clinical practice, including bisphosphonates, selective estrogen receptor modulators, calcitonin, molecular-targeted drugs and so on. However, many drugs have major adverse effects or are unsuitable for long-term use. Therefore, it is very urgent to find more effective therapeutic drugs based on the new pathogenesis of osteoporosis. In this review, we summarize novel mechanisms involved in the pathological process of osteoporosis, including the roles of gut microbiome, autophagy, iron balance and cellular senescence. Based on the above pathological mechanism, we found promising drugs for osteoporosis treatment, such as: probiotics, alpha-ketoglutarate, senolytics and hydrogen sulfide. This new finding may provide an important basis for elucidating the complex pathological mechanisms of osteoporosis and provide promising drugs for clinical osteoporosis treatment.

Shared Genetic Architecture Between Rheumatoid Arthritis and Varying Osteoporotic Phenotypes

Rheumatoid arthritis (RA) and low bone mineral density (BMD), an indicator of osteoporosis (OP), appear epidemiologically associated. Shared genetic factors may explain this association. This study aimed to investigate the presence of pleiotropy to clarify the potential genetic association between RA and OP. We examined BMDs at varying skeletal sites reported in UK Biobank as well as OP fracture acquired from the Genetic Factors for Osteoporosis (GEFOS) Consortium and the TwinsUK study. PRSice-2 was used to assess the potential shared genetic overlap between RA and OP. The presence of pleiotropy was examined using colocalization analysis. PRSice-2 revealed that RA was significantly associated with OP fracture (β = 351.6 ± 83.9, p value = 2.76E-05), total BMD (β = -1763.5 ± 612.8, p = 4.00E-03), spine BMD (β = -919.8 ± 264.6, p value = 5.09E-04), and forearm BMD (β = -66.09 ± 31.40, p value = 3.53E-02). Through colocalization analysis, the same causal genetic variants, associated with both RA and OP, were apparent in 12 genes: PLCL1, BOLL, AC011997.1, TNFAIP3, RP11-158I9.1, CDK6, CHCHD4P2, RP11-505C13.1, PHF19, TRAF1, C5, and C11orf49 with moderate posterior probabilities (>50%). Pleiotropy is involved in the association between RA and OP phenotypes. These findings contribute to the understanding of disease mechanisms and provide insight into possible therapeutic advancements and enhanced screening measures. © 2021 American Society for Bone and Mineral Research (ASBMR).

Associations of osteoprotegerin (OPG) TNFRSF11B gene polymorphisms with risk of fractures in older adult populations: meta-analysis of genetic and genome-wide association studies

The meta-analysis of osteoprotegerin (OPG) (TNFRSF11B) polymorphisms from genetic association studies and genome-wide association studies was performed in order to test the hypothesis of association between OPG polymorphisms and fracture. The findings showed a significant 13% to 37% protective effect of OPG on fractures in postmenopausal women (PSM) (rs2073618), overall, ≥ 60y and Western subjects (rs3134069 and rs3134070).

PURPOSE

Fractures in older people usually result from compromised bone integrity. The multifactorial aetiology of fractures includes both genetic and environmental factors. Inconsistency of reported associations of osteoprotegerin (OPG) (TNFRSF11B) polymorphisms with fracture in the older adult population warranted a meta-analysis to determine more precise estimates.

METHODS

We searched for all available literature on OPG (TNFRSF11B) and fracture. Four polymorphisms were examined, one exonic (rs2073618) and three intronic (rs3134069, rs3134070 and rs3102735). The first two intron polymorphisms were combined (OPGI: osteoprotegerin intron) on account of complete linkage disequilibrium. Risks were estimated with odds ratios (ORs) and 95% confidence intervals (CIs) using the allele-genotype model that included variant (var), wild-type (wt) and heterozygote (het). Multiple comparisons were Bonferroni-corrected. We used meta-regression to examine sources of heterogeneity. Zero heterogeneity (homogeneity: I² = 0%) and high significance (P^a < 0.00001) were the criteria for strength of evidence. Significant outcomes were subjected to sensitivity analysis and publication bias assessment.

RESULTS

From 13 articles (11 genetic association and two genome-wide), this meta-analysis generated five significant pooled ORs, all indicating reduced risks (ORs 0.44-0.87). Of the five, four highly significant comparisons (P^a ≤ 0.00001-0.002) survived the Bonferroni correction, one in rs2073618 het model of the postmenopausal women (OR 0.87, 95% CI 0.81-0.92, I² = 0%) and the other three in OPGI wt model of the overall analysis, ≥ 60 y and Western subjects (ORs 0.63-0.71, 95% CI 0.47-0.86, I² = 97-99%). These findings were consistent, had high significance and high statistical power and were robust and without evidence of publication bias. Four covariates (year of publication, study quality, fracture type/site and sample size) were the sources of heterogeneity in the OPGI overall outcomes (P^a = 0.0001-0.03).

CONCLUSION

Evidence showed that the OPG (TNFRSF11B) polymorphisms reduced the risk for fracture in older adults, particularly protective among postmenopausal women, ≥ 60 y and Western subjects.

The genetic polymorphisms of key genes in WNT pathway (LRP5 and AXIN1) was associated with osteoporosis susceptibility in Chinese Han population

BACKGROUND

Genetic factors play a critical role in the pathogenesis of osteoporosis. The imbalance of WNT/β-catenin will cause the occurrence of osteoporosis. LRP5 and AXIN1 play an important role in the classical Wnt/β-catenin signaling pathway. Our study was aimed to determine the association between five candidate single nucleotide polymorphisms (SNPs) of LRP5 or AXIN1 and osteoporosis susceptibility in Chinese Han population.

METHODS

A total of 599 osteoporosis patients and 599 healthy individuals were recruited for this case-control study. Agena MassARRAY was used to genotype SNPs. The association between SNPs and osteoporosis susceptibility in different genetic models was analyzed by PLINK software. We used false-positive report probability (FPRP) analysis to detect whether the positive results were just chance or noteworthy observations. Multifactor dimension reduction (MDR) was used to analyze the interaction of SNP-SNP in the osteoporosis risk. Finally, haplotype analysis was performed by plink1.07 and Haploview software.

RESULTS

We found that LRP5 rs11228240, AXIN1 rs2301522, and rs9921222 were significantly associated with the osteoporosis susceptibility. The results of subgroup analysis showed that LRP5 rs11228240 (protective factor) and AXIN1 rs2301522 (risk factor) were associated with the susceptibility of osteoporosis among participants who were age >60 years, female or BMI ≤ 24; AXIN1 rs9921222 significantly increased the risk of osteoporosis among participants with BMI ≤ 24. The genotype A_rs2301522C_rs9921222 could increase the susceptibility of osteoporosis (p = 0.026). The rs11228219_LPR5, rs11228240 _LPR5, rs2301522_AXIN1, and rs9921222_AXIN1 four-site model was the best model for predicting the osteoporosis risk (test accuracy = 0.541; CVC = 10/10).

CONCLUSIONS

The LRP5-rs11228240, AXIN1-rs2301522, and AXIN1- rs9921222 were associated with osteoporosis susceptibility in Chinese Han population.

Diagnostic yield of bone fragility gene panel sequencing in children and young adults referred for idiopathic primary osteoporosis at a single regional reference centre

Aim

To describe the presenting features, bone characteristics and molecular genetics in a large monocentric cohort of children and young adults with idiopathic primary osteoporosis.

Methods

Sixty-six patients (19 children, 47 adults; 28 males, 38 females; age at referral: 3.8 to 65 years) diagnosed with primary osteoporosis were included in this study; patients with features of osteogenesis imperfecta or other known syndromes associated with osteoporosis were excluded. For each patient, the following data were collected by retrospective chart review: family and personal history of fracture and osteoporosis, mineral homeostasis parameters and markers of bone formation and resorption, bone mineral density (BMD) of the lumbar spine (LS-BMD), the total body less head (TB-BMD), and total hip levels (TH-BMD) measured by DXA. As part of the initial assessment process, a bone fragility gene panel sequencing was performed in all of these patients.

Results

There was a higher predominance of males in the children (63%) and of females in the adults (66%) (p = 0.030). Compared to the adults, the children had a significantly lower frequency of vertebral fractures (26 vs 57%, p = 0.022) and a higher frequency of peripheral fractures (84 vs 53%; p = 0.019). Bone fragility gene panel sequencing allowed the identification of the heterozygous pathogenic variant in 27% of patients (most frequently in LRP5, WNT1 and COL1A1 or 2 genes) and the heterozygous p.(Val667Met) LRP5 variant in 11% of them. The frequency of pathogenic variants tended to be higher in the children compared to the adults without reaching statistical significance (42 vs 19%; p = 0.053). The frequency of the p.(Val667Met) LRP5 variant was similar in children and adults. No significant differences were found regarding the various clinical, biological and radiological characteristics of the patients according to genotype.

Conclusion

In this study, we reported the presenting features and bone characteristics in a large cohort of children and young adults with idiopathic primary osteoporosis. Bone fragility gene panel sequencing allowed the identification of genetic variants in a significant proportion of these patients. Molecular diagnosis in these patients is important in order to be able to offer genetic counselling and organise patient management.

Predictors of decreased bone mineral density in childhood systemic lupus erythematosus: possible role of osteoprotegerin gene polymorphisms

OBJECTIVES

This study aims to explore effects of osteoprotegerin (OPG) gene polymorphisms and other possible factors on bone mineral density (BMD) in children with systemic lupus erythematosus (SLE).

METHODS

Osteoprotegerin gene rs2073617 and rs3134069 were evaluated in 74 SLE patients and 100 controls then genotypes, alleles and haplotypes' frequencies were compared between cases and controls and between patients with BMD z-scores above and below -2 evaluated by dual energy X-ray absorptiometry (DEXA). Disease activity was evaluated by SLE disease activity index (SLEDAI).

RESULTS

The patients aged 14.01 ± 2.6 years and included 57 (77%) females and 27 (36%) patients with BMD z-score below -2. Genotypes, alleles, and haplotypes frequencies did not differ between patients and controls (p>0.05 for all). Rs3134069 GG genotype and G allele (p=0.001, 0.002) and rs2073617 TT genotype and T allele (p=0.01, 0.006) were significantly higher in patients with BMD below -2. Cumulative glucocorticoids dose, disease duration, and SLEDAI scores were higher in patients with BMD below -2 (p=0.01, 0.01, <0.001, respectively). Regression analysis showed T allele of rs2073617, duration of illness (above 36 months), and cumulative SLEDAI (above 10) as independent predictors of decreased BMD (p 0.02, 0.003, and 0.002, respectively).

CONCLUSIONS

This is the first study to demonstrate OPG gene influence on BMD in children with SLE. The studied SNPs are not risk for developing SLE but, rs2073617 T allele is a possible predictor for reduced BMD in SLE. Other predictors include long disease duration and high activity supporting that osteoporosis in SLE is multifactorial.

Wnt Pathway Extracellular Components and Their Essential Roles in Bone Homeostasis

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

The genetic overlap between osteoporosis and craniosynostosis

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

No genetic causal association between iron status and osteoporosis: A two-sample Mendelian randomization

OBJECTIVE

To explore the genetic causal association between osteoporosis (OP) and iron status through Mendelian randomization (MR).

METHODS

Publicly available genome-wide association study (GWAS) summary data were used for MR analysis with four iron status-related indicators (ferritin, iron, total iron binding capacity, and transferrin saturation) as exposures and three different types of OP (OP, OP with pathological fracture, and postmenopausal OP with pathological fracture) as outcomes. The inverse-variance weighted (IVW) method was used to analyze the genetic causal association between the four indicators of iron status and OP. The heterogeneity of MR results was determined using IVW and MR-Egger methods. The pleiotropy of MR results was determined using MR-Egger regression. A leave-one-SNP-out test was performed to determine whether the MR results were affected by a single nucleotide polymorphism (SNP). The weighted median method was conducted to further validate our results.

RESULTS

Based on IVW, MR-Egger and weighted median models, we found no causal association between iron status (ferritin, iron, total iron binding capacity, or transferrin saturation) and OP (P_beta > 0.05 in all models). IVW and MR-Egger analysis of OP with pathological fracture and iron status indicators showed no potential genetic causal association (P_beta> 0.05 in the two analyses). The results of the weighted median were consistent with those of IVW (P_beta> 0.05 in all analyses). There was no potential genetic causal association between iron status and postmenopausal OP with pathological fracture based on serum iron (P_beta>0.05 in all models). No heterogeneity or horizontal pleiotropy was found in any of the analyses. None of the leave-one-out tests in the analyses found any SNP that could affect the results of MR.

CONCLUSION

Our results demonstrate that there is no genetic causal association between OP and iron status, but the effects of other factors were not excluded.

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

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

Wnt7b Inhibits Osteoclastogenesis via AKT Activation and Glucose Metabolic Rewiring

The imbalance between bone formation and bone resorption causes osteoporosis, which leads to severe bone fractures. It is known that increases in osteoclast numbers and activities are the main reasons for increasing bone resorption. Although extensive studies have investigated the regulation of osteoclastogenesis of bone marrow macrophages (BMMs), new pharmacological avenues still need to be unveiled for clinical purpose. Wnt ligands have been widely demonstrated as stimulators of bone formation; however, the inhibitory effect of the Wnt pathway in osteoclastogenesis is largely unknown. Here, we demonstrate that Wnt7b, a potent Wnt ligand that enhances bone formation and increases bone mass, also abolishes osteoclastogenesis in vitro. Importantly, enforced expression of Wnt in bone marrow macrophage lineage cells significantly disrupts osteoclast formation and activity, which leads to a dramatic increase in bone mass. Mechanistically, Wnt7b impacts the glucose metabolic process and AKT activation during osteoclastogenesis. Thus, we demonstrate that Wnt7b diminishes osteoclast formation, which will be beneficial for osteoporosis therapy in the future.

Regulatory SNP of RREB1 is Associated With Bone Mineral Density in Chinese Postmenopausal Osteoporosis Patients

Postmenopausal osteoporosis (PMO) is the most common bone disorder in elderly Chinese women. Although genetic factors have been shown to have a pivotal role in PMO, studies on genetic loci associated with PMO in Chinese individuals are still lacking. We aimed to identify SNPs that contribute to PMO in Chinese individuals by conducting a genome-wide association study (GWAS). Bone mineral density (BMD) of postmenopausal Chinese women was assessed. Participants with T-score < -2.5 standard deviations (n = 341) were recruited and divided into a discovery group (n = 150) and a replication group (n = 191). GWAS was performed, with T-score as the quantitative trait, using linear regression. Our results revealed that an SNP cluster upstream of RREB1 showed a trend of association with BMD in Chinese PMO patients. The leading SNP of the cluster was rs475011 (p _combined = 1.15 × 10^-6, beta = 0.51), which is a splicing quantitative trait locus (sQTL) of RREB1. This association was further supported by data from the UK Biobank (UKBB; p = 9.56 × 10^-12). The high BMD-associated allele G of rs475011 is related to a high intron excision ratio. This SNP may increase BMD by upregulating mature RREB1 mRNA, based on data from the Genotype-Tissue Expression (GTEx) database. We identified BMD-associated SNPs that regulate RREB1 in Chinese PMO patients. Future functional experiments are needed to further link rs475011, RREB1, and PMO in Chinese individuals.

Associations of LRP5 and MTHFR Gene Variants with Osteoarthritis Prevalence in Elderly Women: A Japanese Cohort Survey Randomly Sampled from a Basic Resident Registry

Objective

Osteoarthritis (OA) is a common and degenerative joint disorder in the elderly. A greater importance of understanding the relationship between genetic factors and OA prevalence has emerged with population aging. We therefore investigated the associations of several bone disease-related genetic variants with the prevalence of OA and osteoporosis in Japanese elderly women from the Obuse study cohort, which was randomly sampled from a basic town resident registry.

Methods and Results

In total, 206 female participants (mean ± standard deviation age: 69.7 ± 11.0 years) who completed OA, bone mineral density, and genotype assessments were included. The number of patients diagnosed as having knee/hip OA and osteoporosis was 59 (28.6%) and 30 (14.6%), respectively. Fisher’s exact testing revealed significant relationships between the minor T allele of LDL receptor related protein 5 (LRP5) rs3736228 and the prevalence of knee/hip OA and osteoporosis. The respective odds ratios (ORs) of the TT genotype for knee/hip OA and osteoporosis were 7.28 (95% confidence interval [CI] 2.22–28.08) and 5.24 (95% CI 0.95–26.98). An additional subgroup analysis for knee OA revealed that the frequency of the common C allele of methylenetetrahydrofolate reductase (MTHFR) rs1801133 had a statistically significant protective association with the prevalence of knee OA (OR 0.58, 95% CI 0.35–0.97).

Conclusion

In sum, the present study demonstrated significant associations of LRP5 rs3736228 and MTHFR rs1801133 with knee/hip OA and osteoporosis prevalences and knee OA prevalence, respectively, in Japanese elderly women. These results will help further the understanding of OA pathogenesis and related genetic risk factors.

Bone mineral density in high-level endurance runners: Part B-genotype-dependent characteristics

Purpose

Inter-individual variability in bone mineral density (BMD) exists within and between endurance runners and non-athletes, probably in part due to differing genetic profiles. Certainty is lacking, however, regarding which genetic variants may contribute to BMD in endurance runners and if specific genotypes are sensitive to environmental factors, such as mechanical loading via training.

Method

Ten single-nucleotide polymorphisms (SNPs) were identified from previous genome-wide and/or candidate gene association studies that have a functional effect on bone physiology. The aims of this study were to investigate (1) associations between genotype at those 10 SNPs and bone phenotypes in high-level endurance runners, and (2) interactions between genotype and athlete status on bone phenotypes.

Results

Female runners with P2RX7 rs3751143 AA genotype had 4% higher total-body BMD and 5% higher leg BMD than AC + CC genotypes. Male runners with WNT16 rs3801387 AA genotype had 14% lower lumbar spine BMD than AA genotype non-athletes, whilst AG + GG genotype runners also had 5% higher leg BMD than AG + GG genotype non-athletes.

Conclusion

We report novel associations between P2RX7 rs3751143 genotype and BMD in female runners, whilst differences in BMD between male runners and non-athletes with the same WNT16 rs3801387 genotype existed, highlighting a potential genetic interaction with factors common in endurance runners, such as high levels of mechanical loading. These findings contribute to our knowledge of the genetic associations with BMD and improve our understanding of why some runners have lower BMD than others.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00421-021-04789-z.

Common Variants of the OPG gene Are Associated with Osteoporosis Risk: A Meta-Analysis

Background: The RANKL/RANK/OPG signaling pathway plays a critical role in osteoclastogenesis and bone remodeling. The associations between sequence variants of the osteoprotegerin (OPG) gene and osteoporosis risk have been widely investigated but remain inconclusive. Objective: We performed a meta-analysis to evaluate the associations between OPG gene variants and osteoporosis risk. Methods: We searched electronic databases and included studies meeting inclusion criteria. The genetic associations of four common OPG variants, A163G, T245G, T950C, and G1181C, with osteoporosis risk were explored. Pooled odds ratios (OR) and 95% confidence intervals (95% CI) were calculated for multiple genetic models. Subgroup analyses including ethnicity, gender, menopausal status, sample size, and control source were also performed. Results: Twenty-six studies comprising 4879 osteoporosis cases and 5708 controls were included. The A163G variant was found to be significantly associated with an increased risk of osteoporosis under both the allelic (G vs. A: OR = 1.45, 95% CI 1.29-1.64, p < 0.001) and dominant models (GG+GA vs. AA: OR = 1.48, 95% CI 1.29-1.70, p < 0.001). Significant associations were also found between the T245G variant and osteoporosis risk. In addition, we observed a reduced risk of osteoporosis in women with the CC genotype at T950C (OR = 0.76, 95% CI 0.64-0.89, p = 0.001) and among Caucasians with the GG or CG genotypes at the G1181C locus (OR = 0.78, 95% CI 0.64-0.94, p = 0.010). In postmenopausal women, only the GG/GA genotypes at the A163G variant were more predisposed to osteoporosis (OR = 1.31, 95% CI 1.00-1.71), whereas CC/CG carriers of G1181C locus may have reduced risk (OR = 0.83, 95% CI 0.66-1.03). Conclusions: Common variants of the OPG gene are associated with osteoporosis risk, especially in the Caucasian population and in the female subgroup. These genetic markers could potentially be used as predictive markers for osteoporosis.

Lrp5 Mutant and Crispant Zebrafish Faithfully Model Human Osteoporosis, Establishing the Zebrafish as a Platform for CRISPR-Based Functional Screening of Osteoporosis Candidate Genes

Genomewide association studies (GWAS) have improved our understanding of the genetic architecture of common complex diseases such as osteoporosis. Nevertheless, to attribute functional skeletal contributions of candidate genes to osteoporosis-related traits, there is a need for efficient and cost-effective in vivo functional testing. This can be achieved through CRISPR-based reverse genetic screens, where phenotyping is traditionally performed in stable germline knockout (KO) mutants. Recently it was shown that first-generation (F0) mosaic mutant zebrafish (so-called crispants) recapitulate the phenotype of germline KOs. To demonstrate feasibility of functional validation of osteoporosis candidate genes through crispant screening, we compared a crispant to a stable KO zebrafish model for the lrp5 gene. In humans, recessive loss-of-function mutations in LRP5, a co-receptor in the Wnt signaling pathway, cause osteoporosis-pseudoglioma syndrome. In addition, several GWAS studies identified LRP5 as a major risk locus for osteoporosis-related phenotypes. In this study, we showed that early stage lrp5 KO larvae display decreased notochord mineralization and malformations of the head cartilage. Quantitative micro-computed tomography (micro-CT) scanning and mass-spectrometry element analysis of the adult skeleton revealed decreased vertebral bone volume and bone mineralization, hallmark features of osteoporosis. Furthermore, regenerating fin tissue displayed reduced Wnt signaling activity in lrp5 KO adults. We next compared lrp5 mutants with crispants. Next-generation sequencing analysis of adult crispant tissue revealed a mean out-of-frame mutation rate of 76%, resulting in strongly reduced levels of Lrp5 protein. These crispants generally showed a milder but nonetheless highly comparable skeletal phenotype and a similarly reduced Wnt pathway response compared with lrp5 KO mutants. In conclusion, we show through faithful modeling of LRP5-related primary osteoporosis that crispant screening in zebrafish is a promising approach for rapid functional screening of osteoporosis candidate genes. © 2021 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

The Effect of Individual Musculoskeletal Conditions on Depression: Updated Insights From an Irish Longitudinal Study on Aging

Few longitudinal studies have systematically investigated whether or how individual musculoskeletal conditions (IMCs) convey risks for negative psychological health outcomes, and approaches to assess such risk in the older population are lacking. In this Irish nationally representative longitudinal prospective study of 6,715 individuals aged 50 and above, machine learning algorithms and various models, including mediation models, were employed to elaborate the underlying mechanisms of IMCs leading to depression and to develop an IMC-induced negative psychological risk (IMCPR) classification approach. Resultantly, arthritis [odds ratio (95% confidence interval): 2.233 (1.700-2.927)], osteoporosis [1.681 (1.133-2.421)], and musculoskeletal chronic pain [MCP, 2.404 (1.838-3.151)] were found to increase the risk of depression after 2 years, while fracture and joint replacement did not. Interestingly, mediation models further demonstrated that arthritis per se did not increase the risk of depression; such risk was augmented only when arthritis-induced restrictions of activities (ARA) existed [proportion of mediation: 316.3% (ARA of usual), 213.3% (ARA of social and leisure), and 251.3% (ARA of sleep)]. The random forest algorithm attested that osteoarthritis, not rheumatoid arthritis, contributed the most to depressive symptoms. Moreover, bone mineral density was negatively associated with depressive symptoms. Systemic pain contributed the most to the increased risk of depression, followed by back, knee, hip, and foot pain (mean Gini-Index: 3.778, 2.442, 1.980, 1.438, and 0.879, respectively). Based on the aforementioned findings, the IMCPR classification approach was developed using an interpretable machine learning model, which stratifies participants into three grades. Among the IMCPR grades, patients with a grade of "severe" had higher odds of depression than those with a "mild" [odds ratio (95% confidence interval): 4.055 (2.907-5.498)] or "moderate" [3.584 (2.101-5.883)] grade. Females with a "severe" grade had higher odds of depression by 334.0% relative to those with a "mild" grade, while males had a relative risk of 258.4%. In conclusion, the present data provide systematic insights into the IMC-induced depression risk and updated the related clinical knowledge. Furthermore, the IMCPR classification approach could be used as an effective tool to evaluate this risk.

The association of OPG polymorphisms with risk of osteoporotic fractures: A systematic review and meta-analysis

BACKGROUND

Subjects with low bone mineral density and osteoporosis are more likely to suffer osteoporotic fractures during their lifetime. Polymorphisms in osteoprotegerin (OPG) gene are found to be associated with low bone mineral density and osteoporosis risk but their association with fracture risk is inconclusive. Here, we performed a meta-analysis to investigate the relationship between OPG polymorphisms with susceptibility to osteoporotic fractures.

METHODS

Eligible studies investigating the association between common OPG polymorphisms (A164G, T245G, T950C, and G1181C) and risk of osteoporotic fracture were retrieved from PubMed, EMBASE, Web of Science, and the Cochrane Library. Odds ratio (OR) and the 95% confidence interval (CI) were calculated in the allelic, dominant, recessive, and homozygous model. Subgroup analyses of vertebral fractures, Caucasians, and postmenopausal women were also performed.

RESULTS

A total of 14 studies comprising 5459 fracture cases and 9860 non-fracture controls were included. A163G was associated with fracture risk in dominant (OR = 1.29, 95%CI 1.11-1.50), recessive (OR = 1.64, 95%CI 1.10-2.44), and homozygous model (OR = 1.73, 95%CI 1.16-2.59). T245G was significantly correlated with susceptibility to fractures in all genetic models. Subjects with CC genotype of T950C had a reduced risk of fracture compared to those with CT or TT genotypes (OR = 0.81, 95%CI 0.70-0.94, P = .004). Subgroup analysis showed that A163G and T245G but not T950C and G1181C were associated with vertebral fracture risk.

CONCLUSION

OPG A163G and T245G polymorphisms were risk factors of osteoporotic fractures while T950C had a protective role. These polymorphisms can be used as predictive markers of fractures.

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

OBJECTIVE

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

MATERIALS AND METHODS

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

CONCLUSIONS

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

SIGNIFICANCE

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

LIMITATIONS

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

SUGGESTIONS FOR FUTURE RESEARCH

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

Maternal High-Fat Diet Induces Long-Lasting Defects in Bone Structure in Rat Offspring Through Enhanced Osteoclastogenesis

Maternal stressors during the prenatal and perinatal periods are associated with increased susceptibility for and severity of chronic disease phenotypes in adult offspring. In this study, we used a rat model of maternal high-fat diet (HFD) exposure during pregnancy and lactation to investigate the impact on skeletal homeostasis in offspring. In the distal femur, young male and female offspring (up to 3 weeks of age) from dams fed a HFD exhibited marked increases in trabecular bone volume relative to offspring from dams fed a chow diet, but this was followed by sustained bone loss. By 15 weeks of age, male offspring of HFD fed dams exhibited a 33% reduction in trabecular bone volume fraction that histomorphometric analyses revealed was due to a nearly threefold increase in the abundance of bone-resorbing osteoclasts, while there were no differences between female control and HFD offspring by 15 weeks of age. The osteoblastic differentiation of male offspring-derived bone marrow stromal cells was not affected by maternal diet. However, osteoclastic precursors isolated from the male offspring of HFD fed dams exhibited enhanced differentiation in vitro, forming larger osteoclasts with higher expression of the fusion marker DC-STAMP. This effect appears to be mediated by a cell autonomous increase in the sensitivity of precursors to RANKL. Taken together, these results suggest that maternal stressors like HFD exposure have persistent consequences for the skeletal health of offspring that may ultimately lead to a predisposition for osteopenia/osteoporosis.

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.

Swimming as Treatment for Osteoporosis: A Systematic Review and Meta-analysis

Osteoporosis is a chronic disease that seriously affects human health and quality of life. This study is aimed at determining whether swimming had an effect on the bone mineral density (BMD) of the spine and femoral neck in postmenopausal and premenopausal osteoporosis patients. We retrieved relevant literature and analyzed data from randomized controlled trials to assess the effect of swimming on BMD in postmenopausal and premenopausal women. Relevant studies, with no language restrictions, from inception to September 2019, were retrieved from the PubMed, Cochrane, EMBASE, and EBSCO databases independently by two investigators. The keywords used for the literature search were “osteoporosis” and “swimming.” The main results included BMD and T-score. We searched 256 relevant articles and finally screened five articles, including 263 participants. Lumbar spine density was mentioned in three articles. Although the heterogeneity of lumbar vertebral density is moderate, the analysis of swimmers to nonswimmers shows that the lumbar vertebral density in swimmers is improved [heterogeneity: chi² = 5.16, df = 2 (P = 0.08); I² = 61%]. We analyzed the following heterogeneous subgroups: subgroup 1 (3–6 hours) and subgroup 2 (<3 hours). The BMD in subgroup 1 was significantly higher than that in the placebo, while no effect on BMD was found in subgroup 2 [heterogeneity: chi² = 0.15, df = 3 (P = 0.70); I² = 0%]. According to the current evidence, swimming may improve the BMD of postmenopausal women participants, if the swimming time is between 3 and 6 hours, especially in long-term swimmers. However, the effectiveness of swimming does require further investigation.