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

The GWAS candidate far upstream element binding protein 3 (FUBP3) is required for normal skeletal growth, and adult bone mass and strength in mice

Bone mineral density (BMD) and height are highly heritable traits for which hundreds of genetic loci have been linked through genome wide association studies (GWAS). FUBP3 is a DNA and RNA binding protein best characterised as a transcriptional regulator of c-Myc, but little is known about its role in vivo. Single nucleotide polymorphisms in FUBP3 at the 9q34.11 locus have been associated with BMD, fracture and height in multiple GWAS, but FUBP3 has no previously established role in the skeleton. We analysed Fubp3-deficient mice to determine the consequence of FUBP3 deficiency in vivo. Mice lacking Fubp3 had reduced survival to adulthood and impaired skeletal growth. Bone mass was decreased, most strikingly in the vertebrae, with altered trabecular micro-architecture. Fubp3 deficient bones were also weak. These data provide the first functional demonstration that Fubp3 is required for normal skeletal growth and development and maintenance of adult bone structure and strength, indicating that FUBP3 contributes to the GWAS association of 9q34.11 with variation in height, BMD and fracture.

Causal association between female reproductive factors and periodontitis: evidence from Mendelian randomization study

OBJECTIVES

To evaluate the causal relationship between female reproductive factors and periodontitis, and identify potential mediators underlying these associations.

MATERIALS AND METHODS

The genome-wide association study (GWAS) summary statistics from different consortiums was used to search for single nucleotide polymorphisms (SNPs) strongly associated with exposures as instrumental variables (IVs). The causality between female reproductive factors and periodontitis was measured by univariable mendelian randomization (UVMR). The mediation analysis was conducted by two-step mendelian randomization (Two-Step MR). The reverse MR was utilized to examine potential reverse causality.

RESULTS

The age of menarche (β 0.1328; 95%CI 0.0454-0.2203; p 0.0029) and the age of first birth (β 0.1222; 95%CI 0.171-0.2247; p 0.0227) were positively correlated with chronic periodontitis. And the mediation analysis suggested that the heel bone mineral density (BMD) mediated the correlation between the age of menarche and chronic periodontitis.

CONCLUSIONS

The age of menarche and the age of first birth have causal relationships with the infection risk of chronic periodontitis. Further clinical and laboratory researches are needed to clarify the mechanism behind the newly found causality.

CLINICAL RELEVANCE

This study suggested that the later the age of menarche and first labor, the higher the susceptibility to chronic periodontitis, which may help the clinicians strengthen the prevention and treatment of periodontitis in women.

Causal Associations of Inflammatory Cytokines With Osteosarcopenia: Insights From Mendelian Randomization and Single Cell Analysis

Background: Osteosarcopenia, the coexistence of osteoporosis and sarcopenia, poses significant challenges in aging populations due to its dual impact on bone and muscle health. Inflammation, mediated by specific cytokines, is thought to play a crucial role in the development of osteosarcopenia, though the underlying mechanisms are not fully understood. Objective: This study aimed to clarify the causal role of circulating cytokines in the pathogenesis of osteosarcopenia by employing mendelian randomization (MR) and single-cell RNA sequencing (scRNA-seq) to identify cell-specific cytokine expression patterns. The ultimate objective was to uncover potential pathological mechanisms and therapeutic targets for treating osteosarcopenia. Methods: A two-sample MR approach was employed, leveraging publicly available genome-wide association study (GWAS) data from multiple cohorts. A total of 91 circulating cytokines were examined using genetic instruments, and their causal effects on traits related to osteoporosis and sarcopenia were evaluated. Various complementary and sensitivity analyses were performed to ensure robust findings. Additionally, scRNA-seq datasets from human muscle and bone marrow were analyzed to validate the single-cell expression profiles of candidate cytokines. Results: MR analysis identified several cytokines with causal effects on osteosarcopenia traits, including LTA, CD40, CXCL6, CXCL10, DNER (delta and notch-like epidermal growth factor-related receptor), and VEGFA (vascular endothelial growth factor A). LTA and CD40 were protective for both bone and muscle, while VEGFA posed a risk. Other cytokines demonstrated opposite effects on bone and muscle. Single cell analysis revealed distinct expression patterns, with LTA highly expressed in lymphocytes, CD40 in immune cells, and VEGFA in various musculoskeletal cell types. Age-related differences in cytokine expression were also noted, with LTA more highly expressed in younger individuals, and VEGFA in older individuals. Conclusion: This study offers preliminary insights into the inflammatory mechanisms potentially driving osteosarcopenia, identifying key cytokines that may be involved in its pathogenesis. By integrating MR and scRNA-seq data, we highlight potential therapeutic targets, though further research is needed to confirm these findings and their implications for musculoskeletal health.

Identifying rare variants in genes related to bone phenotypes in a cohort of postmenopausal women

Rare genetic variants in genes previously described to be involved in bone monogenic disorders were identified in postmenopausal women split into two groups according to extreme bone mineral density (BMD) values and lumbar spine Z-scores. A pathogenic variant in COL1A2 gene found in a woman with low BMD highlights the overlap between osteogenesis imperfecta and osteoporosis, which may share their genetic etiology. Other variants were not clearly associated with the extreme BMD, suggesting that there is little contribution of rare variants to postmenopausal osteoporosis.

PURPOSE

We aimed to evaluate whether extreme values of bone mineral density (BMD) in a population-based cohort of postmenopausal women (BARCOS) could be determined by rare genetic variants in genes related to monogenic bone disorders.

METHODS

A panel of 127 genes related to different skeletal phenotypes was designed. Massive sequencing by targeted capture of these genes was performed in 104 DNA samples from those women of the BARCOS cohort that exhibited the highest (HZ group) and lowest (LZ group) LS Z-scores, ranging from + 0.70 to + 3.80 and from - 2.35 to - 4.26, respectively. 5'UTR, 3'UTR, splice region, missense, nonsense, and short indel variants with MAF < 0.01 were annotated with CADD version 1.6 and considered in the analysis.

RESULTS

After filtering those variants with CADD > 25 and present only in one of the groups (either LZ or HZ), six variants were detected, most of which (5/6) were in the LZ group in TCIRG1, COL1A2, SEC24D, LRP4, and ANO5 genes, while only one, in the LMNA gene, was in the HZ group. According to the ClinVar database, the COL1A2 variant, causative of a recessive form of osteogenesis imperfecta, is described as pathogenic, while the other variants are considered of uncertain significance (VUS).

CONCLUSION

The variant identified in COL1A2 in a woman from the LZ group highlights the genetic overlap between monogenic diseases such as osteogenesis imperfecta and complex diseases like osteoporosis. However, the other variants were not clearly associated with the extreme BMD, suggesting that there is little contribution of rare variants to postmenopausal osteoporosis.

Regulation of WNT16 in bone may involve upstream enhancers within CPED1

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

Development and validation of a novel prediction model for osteoporosis : from serotonin to fat-soluble vitamins

Aims

We aimed to develop and validate a novel prediction model for osteoporosis based on serotonin, fat-soluble vitamins, and bone turnover markers to improve prediction accuracy of osteoporosis.

Methods

Postmenopausal women aged 55 to 65 years were recruited and divided into three groups based on DXA (normal, osteopenia, and osteoporosis). A total of 109 participants were included in this study and split into healthy (39/109, 35.8%), osteopenia (35/109, 32.1%), and osteoporosis groups (35/109, 32.1%). Serum concentrations of serotonin, fat-soluble vitamins, and bone turnover markers of participants were measured. Stepwise discriminant analysis was performed to identify efficient predictors for osteoporosis. The prediction model was developed based on Bayes and Fisher's discriminant functions, and validated via leave-one-out cross-validation. Normal and empirical volume under the receiver operating characteristic (ROC) surface (VUS) tests were used to evaluate predictive effects of variables in the prediction model.

Results

Significant variables including oestrogen (E2), total procollagen type 1 amino-terminal propeptide (TP1NP), parathyroid hormone (PTH), BMI, vitamin K, serotonin, osteocalcin (OSTEOC), vitamin A, and vitamin D3 were used for the development of the prediction model. The training accuracy for normal, osteopenia, and osteoporosis is 74.4% (29/39), 80.0% (28/35), and 85.7% (30/35), respectively, while the total training accuracy is 79.8% (87/109). The internal validation showed excellent performance with 72.5% testing accuracy (72/109). Among these variables, serotonin and vitamin K exert important roles in the prediction of osteoporosis.

Conclusion

We successfully developed and validated a novel prediction model for osteoporosis based on serum concentrations of serotonin, fat-soluble vitamins, and bone turnover markers. In addition, interactive communication between serotonin and fat-soluble vitamins was observed to be critical for bone health in this study.

Identification of a global gene expression signature associated with the genetic risk of catastrophic fracture in iPSC-derived osteoblasts from Thoroughbred horses

Bone fractures are a significant problem in Thoroughbred racehorses. The risk of fracture is influenced by both genetic and environmental factors. To determine the biological processes that are affected in genetically susceptible horses, we utilised polygenic risk scoring to establish induced pluripotent stem cells (iPSCs) from horses at high and low genetic risk. RNA-sequencing on iPSC-derived osteoblasts revealed 112 genes that were significantly differentially expressed. Forty-three of these genes have known roles in bone, 27 are not yet annotated in the equine genome and 42 currently have no described role in bone. However, many of the proteins encoded by the known and unknown genes have reported interactions. Functional enrichment analyses revealed that the differentially expressed genes were overrepresented in processes regulating the extracellular matrix and pathways known to be involved in bone remodelling and bone diseases. Gene set enrichment analysis also detected numerous biological processes and pathways involved in glycolysis with the associated genes having a higher expression in the iPSC-osteoblasts from horses with low polygenic risk scores for fracture. Therefore, the differentially expressed genes may be relevant for maintaining bone homeostasis and contribute to fracture risk. A deeper understanding of the consequences of mis-regulation of these genes and the identification of the DNA variants which underpin their differential expression may reveal more about the molecular mechanisms which are involved in equine bone health and fracture risk.

Family-based whole-exome sequencing implicates a variant in lysyl oxidase like 4 in atypical femur fractures

Atypical femur fractures (AFFs) are rare adverse events associated with bisphosphonate use, having unclear pathophysiology. AFFs also cluster in families and have occurred in patients with monogenetic bone diseases sometimes without bisphosphonate use, suggesting an underlying genetic susceptibility. Our aim was to identify a genetic cause for AFF in a Caucasian family with 7 members affected by osteoporosis, including 3 siblings with bisphosphonate-associated AFFs. Using whole-exome sequencing, we identified a rare pathogenic variant c.G1063A (p.Gly355Ser) in lysyl oxidase like 4 (LOXL4) among 64 heterozygous rare, protein-altering variants shared by the 3 siblings with AFFs. The same variant was also found in a fourth sibling with a low-trauma femur fracture above the knee, not fulfilling all the ASBMR criteria of AFF and in 1 of 73 unrelated European AFF patients. LOXL4 is involved in collagen cross-linking and may be relevant for microcrack formation and bone repair mechanisms. Preliminary functional analysis showed that skin fibroblast-derived osteoblasts from the unrelated patient with the LOXL4 variant expressed less collagen type I and elastin, while osteogenic differentiation and mineralization were enhanced compared with 2 controls. In conclusion, this LOXL4 variant may underlie AFF susceptibility possibly due to abnormal collagen metabolism, leading to increased formation of microdamage or compromised healing of microcracks in the femur.

Hypoglycemic agents and bone health; an umbrella systematic review of the clinical trials' meta-analysis studies

BACKGROUND

No clear consensus exists regarding the safest anti-diabetic drugs with the least adverse events on bone health. This umbrella systematic review therefore aims to assess the published meta-analysis studies of randomized controlled trials (RCTs) conducted in this field.

METHODS

All relevant meta-analysis studies of RCTs assessing the effects of anti-diabetic agents on bone health in patients with diabetes mellitus (DM) were collected in line with the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA). English articles published until 15 March 2023 were collected through the search of Cochrane Library, Scopus, ISI Web of Sciences, PubMed, and Embase using the terms "Diabetes mellitus", "anti-diabetic drugs", "Bone biomarker", "Bone fracture, "Bone mineral density" and their equivalents. The methodological and evidence quality assessments were performed for all included studies.

RESULTS

From among 2220 potentially eligible studies, 71 meta-analyses on diabetic patients were included. Sodium-glucose cotransporter-2 inhibitors (SGLT-is) showed no or equivalent effect on the risk of fracture. Dipeptidyl peptidase-4 inhibitors (DPP-4is) and Glucagon-like peptide-1 receptor agonists (GLP-1Ras) were reported to have controversial effects on bone fracture, with some RCTs pointing out the bone protective effects of certain members of these two medication classes. Thiazolidinediones (TZDs) were linked with increased fracture risk as well as higher concentrations of C-terminal telopeptide of type I collagen (CTx), a bone resorption marker.

CONCLUSION

The present systematic umbrella review observed varied results on the association between the use of anti-diabetic drugs and DM-related fracture risk. The clinical efficacy of various anti-diabetic drugs, therefore, should be weighed against their risks and benefits in each patient.

Genetic variants associated with longevity in long-living Indians

Genetic factors play a significant role in determining an individual's longevity. The present study was aimed at identifying genetic variants associated with longevity in Indian population. Long living individuals (LLIs), aged 85+, were compared with younger controls, aged 18-49 years, using data from GenomegaDB, a genetic database of Indians living in India. An in-house developed custom chip, having variants associated with various cancers, cardiovascular, neurological, gastro-intestinal, metabolic and auto-immune disorders, was used to generate genotype data. Logistic regression analysis with sex and top three genetic principal components as covariates resulted in 9 variants to be significantly associated with longevity at a p-value threshold of 5 × 10-4. Alleles associated with slower heart rate (rs365990, MYH6), decreased risk of osteoporosis and short body height (rs2982570, ESR1), decreased risk of schizophrenia (rs1339227, RIMS1-KCNQ5) and decreased risk of anxiety and neuroticism (rs391957, HSPA5) were found to have higher frequency in LLIs. Alleles associated with increased risk of atrial fibrillation (rs3903239, GORAB-PRRX1) and biliary disorders (rs2002042, ABCC2) were found to have lower frequency. The G allele of rs2802292 from FOXO3A gene, associated with longevity in Japanese, German and French centenarians, was also found to be significant in this population (P = 0.032). Pathway enrichment analysis revealed that the genes involved in oxidative stress, apoptosis, DNA damage repair, glucose metabolism and energy metabolism were significantly involved in affecting the longevity. Results of our study demonstrate the genetic basis of healthy aging and longevity in the population.

Valid inference for machine learning-assisted genome-wide association studies

Machine learning (ML) has become increasingly popular in almost all scientific disciplines, including human genetics. Owing to challenges related to sample collection and precise phenotyping, ML-assisted genome-wide association study (GWAS), which uses sophisticated ML techniques to impute phenotypes and then performs GWAS on the imputed outcomes, have become increasingly common in complex trait genetics research. However, the validity of ML-assisted GWAS associations has not been carefully evaluated. Here, we report pervasive risks for false-positive associations in ML-assisted GWAS and introduce Post-Prediction GWAS (POP-GWAS), a statistical framework that redesigns GWAS on ML-imputed outcomes. POP-GWAS ensures valid and powerful statistical inference irrespective of imputation quality and choice of algorithm, requiring only GWAS summary statistics as input. We employed POP-GWAS to perform a GWAS of bone mineral density derived from dual-energy X-ray absorptiometry imaging at 14 skeletal sites, identifying 89 new loci and revealing skeletal site-specific genetic architecture. Our framework offers a robust analytic solution for future ML-assisted GWAS.

Causal inference in health and disease: a review of the principles and applications of Mendelian randomization

Mendelian randomization (MR) is a genetic epidemiological technique that uses genetic variation to infer causal relationships between modifiable exposures and outcome variables. Conventional observational epidemiological studies are subject to bias from a range of sources; MR analyses can offer an advantage in that they are less prone to bias as they use genetic variants inherited at conception as "instrumental variables", which are proxies of an exposure. However, as with all research tools, MR studies must be carefully designed to yield valuable insights into causal relationships between exposures and outcomes, and to avoid biased or misleading results that undermine the validity of the causal inferences drawn from the study. In this review, we outline Mendel's laws of inheritance, the assumptions and principles that underlie MR, MR study designs and methods, and how MR analyses can be applied and reported. Using the example of serum phosphate concentrations on liability to kidney stone disease we illustrate how MR estimates may be visualized and, finally, we contextualize MR in bone and mineral research including exemplifying how this technique could be employed to inform clinical studies and future guidelines concerning BMD and fracture risk. This review provides a framework to enhance understanding of how MR may be used to triangulate evidence and progress research in bone and mineral metabolism as we strive to infer causal effects in health and disease.

Cystatin C is associated with osteoporosis and fractures: An observational study based on Mendelian randomization analysis

OBJECTIVES

Osteoporosis is characterized by decreased bone mass and damaged bone microstructure, often leading to fragility fractures. Low bone mineral density is a key risk factor for fractures. Serum cystatin C (CysC), an endogenous marker of glomerular filtration rate, is negatively correlated with bone mineral density and may be a potential risk factor for osteoporosis. This study aims to investigate the association and potential pathogenic mechanisms between CysC and osteoporosis and fractures in the general population by combining cohort analysis and Mendelian randomization (MR) analysis.

METHODS

Large-scale prospective cohort data from the UK Biobank and summary statistics from genome-wide association study (GWAS) in European populations were utilized, with strict exclusion criteria applied (excluding non-white individuals, those with thyroid diseases, gastrointestinal dysfunction, kidney diseases, rheumatoid diseases, malignant tumors, chronic infections or inflammatory diseases, diabetes, hypertension, and individuals taking medications that affect bone metabolism). Multivariable linear regression, logistic regression, and Cox proportional hazards models were used to analyze the relationship between CysC and bone mineral density, osteoporosis, and fracture risk. All analyses were performed using three sequential models to adjust for confounding factors: Model 1 adjusted for demographic characteristics and lifestyle factors; Model 2 further adjusted for renal function based on Model 1; and Model 3 further adjusted for physical activity based on Model 2. Restricted cubic spline models were used to explore non-linear relationships, and MR analysis was conducted to assess the causal associations between CysC and osteoporosis and fractures.

RESULTS

Multivariate analysis showed that after adjusting for basic variables (Model 1), there was no correlation between CysC and estimated bone mineral density (eBMD) in the overall study population; however, when stratified by gender, both males and females exhibited a significant negative correlation (P<0.001). After further adjustment for renal function (Model 2) and physical activity level (Model 3), CysC became negatively correlated with eBMD in the overall population (P<0.001). Moreover, multivariable logistic regression consistently demonstrated that CysC concentration was significantly positively associated with osteoporosis risk (P<0.01), and this association remained stable across all models. In all populations and models, multivariate Cox regression analysis indicated that subjects in the highest quartile (Q4) of CysC had a significantly increased risk of developing osteoporosis (P<0.001). In the overall population, the positive association between Q4 CysC levels and fractures was observed only in Models 2 and 3, with a hazard ratio of 1.118 (both P<0.001); however, after gender stratification, this association disappeared in males (P>0.05). Additionally, restricted cubic spline regression analyses revealed a significant non-linear relationship between CysC and the incidence of osteoporosis and fractures (P<0.05). MR analysis, using 167 single nucleotide polymorphisms (SNPs) as instrumental variables, showed no direct causal relationship between CysC and osteoporosis or fractures (P≥0.05), a finding that differs from previous studies in special populations.

CONCLUSIONS

Elevated levels of CysC are significantly associated with an increased risk of osteoporosis and fractures, and this association is more pronounced in females. Renal function and physical activity levels may be important factors influencing this relationship. The link between CysC and osteoporosis and fractures may be mediated by several mechanisms: Eelevated CysC may lead to abnormalities in vitamin D and mineral metabolism, thereby inhibiting bone formation; renal dysfunction may exacerbate inflammation, affecting bone resorption; or in the osteoporosis state, increased osteoclast differentiation may result in elevated CysC levels. These findings support the potential use of CysC as a biomarker for predicting the risk of osteoporosis.

A GWAS for grip strength in cohorts of children-Advantages of analysing young participants for this trait

Grip strength (GS) is a proxy measure for muscular strength and a predictor for bone fracture risk among other diseases. Previous genome-wide association studies (GWASs) have been conducted in large cohorts of adults focusing on scores collected for the dominant hand, therefore increasing the likelihood of confounding effects by environmental factors. Here, we perform the first GWAS meta-analyses on maximal GS with the dominant (GSD) and non-dominant (GSND) hand in two cohorts of children (ALSPAC, N = 5450; age range = 10.65-13.61; Raine Study, N = 1162, age range: 9.42-12.38 years). We identified a novel significant association for GSND (rs9546244, LINC02465, p = 3.43e-08) and replicated associations previously reported in adults including with a HOXB3 gene marker that shows an expression quantitative trait locus (eQTL) effect. Despite a much smaller sample (~3%) compared with the UK Biobank we replicated correlation analyses previously reported in this much larger adult cohort, such as a negative correlation with coronary artery disease. Although the results from the polygenic risk score (PRS) analyses did not survive multiple testing correction, we observed nominally significant associations between GS and risk of overall fracture, as previously reported, as well ADHD which will require further investigations. Finally, we observed a higher SNP-heritability (24%-41%) compared with previous studies (4%-24%) in adults. Overall, our results suggest that cohorts of children might be better suited for genetic studies of grip strength, possibly due to the shorter exposure to confounding environmental factors compared with adults.

The potential protective effect of 3-Hydroxybutyrate against aortic dissection: a mendelian randomization analysis

BACKGROUND

3-Hydroxybutyrate, also called β-hydroxybutyrate, is a significant constituent of ketone bodies. Previous observational and experimental studies have suggested that ketogenic diet, especially 3-hydroxybutyrate, may have a protective effect against cardiovascular disease. However, the relationship between ketone bodies, especially 3-hydroxybutyrate, and aortic dissection remains uncertain.

MATERIALS AND METHODS

Publicly accessible data from genome-wide association study (GWAS) was utilized to obtain information on ketone bodies, including 3-hydroxybutyrate, acetoacetate and acetone as exposure respectively, while GWAS data on aortic dissection was used as outcome. Subsequently, two-sample Mendelian randomization (MR) analysis was conducted to examine the potential relationship between ketone bodies and aortic dissection. Then, reverse and multivariate Mendelian randomization analyses were performed. Additionally, sensitivity tests were conducted to assess the robustness of MR study.

RESULTS

The inverse-variance weighted (IVW) method of Mendelian randomization analysis of gene prediction observed a negative correlation between 3-hydroxybutyrate and risk of aortic dissection (OR 0.147, 95% CI 0.053-0.410). Furthermore, consistent findings were obtained through the implementation of the weighted median, simple mode, Mendelian randomization-Egger (MR-Egger), and weighted mode methods. After adjusting acetoacetate (OR 0.143, 95% CI 0.023-0.900) or acetone (OR 0.100, 95% CI 0.025-0.398), MR analysis of gene prediction still observed a negative correlation between 3-hydroxybutyrate and risk of aortic dissection. No indications of heterogeneity or pleiotropy among the SNPs were detected.

CONCLUSION

The findings from the MR analysis demonstrated that genetically predicted 3-hydroxybutyrate exhibits a protective effect against aortic dissection.

Mining Candidate Genes and Identifying Risk Factors for Leg Disease in Broilers: A Mendelian Randomization Study

Clinical investigations have highlighted disruptions in bone metabolic processes and abnormal fluctuations in serum indicator levels during the onset of leg disease (LD) in broilers. However, the presence of a genetic causal relationship for this association remains undetermined. Therefore, the aim of this study is to discern the risk factors underlying LD development using 1235 sequenced white-feathered broilers. We employed Mendelian randomization (MR) analysis to assess the associations of bone strength (BS), bone mineral density (BMD), tibial bone weight (TBW), tibial bone length (TBL), tibial bone diameter (TBD), bone ash (BA), ash calcium (Ash Ca), ash phosphorus (Ash P), serum calcium (Ca), serum phosphorus (P), serum alkaline phosphatase (ALP), and serum osteoprotegerin (OPG) with the incidence of LD. Compelling evidence underscores a causal link between the risk of developing LD and decreased BMD (odds ratio (OR) = 0.998; 95% CI: 0.983, 0.993; P < 0.001) and narrower TBD (OR = 0.985, 95% CI: 0.975, 0.994, P = 0.002). Additionally, serum OPG concentrations (OR: 0.995, 95% CI: 0.992, 0.999, P = 0.008) were associated with BMD (OR = 0.0078, 95% CI = 0.0043 to 0.0140, P < 0.001), indicating a robust genetic relationship between ALP concentrations (OR: 0.988, 95% CI: 0.984, 0.993, P < 0.001) and TBD (OR = 0.0046, 95% CI = 0.0026, 0.0083, P < 0.001). Moreover, elevated serum Ca (OR: 0.564, 95% CI: 0.487, 0.655, P < 0.001) and P (OR: 0.614, 95% CI: 0.539, 0.699, P < 0.001) levels were associated with a narrower TBD. Elevated serum levels of Ca, P, ALP, and OPG contribute to disturbances in bone metabolism, while decreased BMD and narrower TBD are associated with a greater risk of developing LD in broilers. This discovery elucidates the metabolic risk factors for LD in broilers and could provide information on LDs, such as osteoporosis, in humans.

Associations between gut microbiota and incident fractures in the FINRISK cohort

The gut microbiota (GM) can regulate bone mass, but its association with incident fractures is unknown. We used Cox regression models to determine whether the GM composition is associated with incident fractures in the large FINRISK 2002 cohort (n = 7043, 1092 incident fracture cases, median follow-up time 18 years) with information on GM composition and functionality from shotgun metagenome sequencing. Higher alpha diversity was associated with decreased fracture risk (hazard ratio [HR] 0.92 per standard deviation increase in Shannon index, 95% confidence interval 0.87-0.96). For beta diversity, the first principal component was associated with fracture risk (Aitchison distance, HR 0.90, 0.85-0.96). In predefined phyla analyses, we observed that the relative abundance of Proteobacteria was associated with increased fracture risk (HR 1.14, 1.07-1.20), while the relative abundance of Tenericutes was associated with decreased fracture risk (HR 0.90, 0.85-0.96). Explorative sub-analyses within the Proteobacteria phylum showed that higher relative abundance of Gammaproteobacteria was associated with increased fracture risk. Functionality analyses showed that pathways related to amino acid metabolism and lipopolysaccharide biosynthesis associated with fracture risk. The relative abundance of Proteobacteria correlated with pathways for amino acid metabolism, while the relative abundance of Tenericutes correlated with pathways for butyrate synthesis. In conclusion, the overall GM composition was associated with incident fractures. The relative abundance of Proteobacteria, especially Gammaproteobacteria, was associated with increased fracture risk, while the relative abundance of Tenericutes was associated with decreased fracture risk. Functionality analyses demonstrated that pathways known to regulate bone health may underlie these associations.

Genetic and Modifiable Risk Factors for Postoperative Complications of Total Joint Arthroplasty: A Genome-Wide Association and Mendelian Randomization Study

Background: Total joint arthroplasty (TJA) is an orthopedic procedure commonly used to treat damaged joints. Despite the efficacy of TJA, postoperative complications, including aseptic prosthesis loosening and infections, are common. Moreover, the effects of individual genetic susceptibility and modifiable risk factors on these complications are unclear. This study analyzed these effects to enhance patient prognosis and postoperative management. Methods: We conducted an extensive genome-wide association study (GWAS) and Mendelian randomization (MR) study using UK Biobank data. The cohort included 2964 patients with mechanical complications post-TJA, 957 with periprosthetic joint infection (PJI), and a control group of 398,708 individuals. Genetic loci associated with postoperative complications were identified by a GWAS analysis, and the causal relationships of 11 modifiable risk factors with complications were assessed using MR. Results: The GWAS analysis identified nine loci associated with post-TJA complications. Two loci near the PPP1R3B and RBM26 genes were significantly linked to mechanical complications and PJI, respectively. The MR analysis demonstrated that body mass index was positively associated with the risk of mechanical complications (odds ratio [OR]: 1.42; p < 0.001). Higher educational attainment was associated with a decreased risk of mechanical complications (OR: 0.55; p < 0.001) and PJI (OR: 0.43; p = 0.001). Type 2 diabetes was suggestively associated with mechanical complications (OR, 1.18, p = 0.02), and hypertension was suggestively associated with PJI (OR, 1.41, p = 0.008). Other lifestyle factors, including smoking and alcohol consumption, were not causally related to postoperative complications. Conclusions: The genetic loci near PPP1R3B and RBM26 influenced the risk of post-TJA mechanical complications and infections, respectively. The effects of genetic and modifiable risk factors, including body mass index and educational attainment, underscore the need to perform personalized preoperative assessments and the postoperative management of surgical patients. These results indicate that integrating genetic screening and lifestyle interventions into patient care can improve the outcomes of TJA and patient quality of life.

A plasma protein-based risk score to predict hip fractures

As there are effective treatments to reduce hip fractures, identification of patients at high risk of hip fracture is important to inform efficient intervention strategies. To obtain a new tool for hip fracture prediction, we developed a protein-based risk score in the Cardiovascular Health Study using an aptamer-based proteomic platform. The proteomic risk score predicted incident hip fractures and improved hip fracture discrimination in two Trøndelag Health Study validation cohorts using the same aptamer-based platform. When transferred to an antibody-based proteomic platform in a UK Biobank validation cohort, the proteomic risk score was strongly associated with hip fractures (hazard ratio per s.d. increase, 1.64; 95% confidence interval 1.53-1.77). The proteomic risk score, but not available polygenic risk scores for fractures or bone mineral density, improved the C-index beyond the fracture risk assessment tool (FRAX), which integrates information from clinical risk factors (C-index, FRAX 0.735 versus FRAX + proteomic risk score 0.776). The developed proteomic risk score constitutes a new tool for stratifying patients according to hip fracture risk; however, its improvement in hip fracture discrimination is modest and its clinical utility beyond FRAX with information on femoral neck bone mineral density remains to be determined.

Genome-wide association study reveals the candidate genes of humerus quality in laying duck

Bone plays a crucial role in poultry's health and production. However, during the selection and cage farming, there has been a decline in bone quality. As the development of breeding theory, researchers find that it's possible to enhance bone quality through selective breeding.This study measure 8 humerus quality in 260 samples of the 350-day-old female duck. By descripting the basic characteristic traits, mechanical property traits we found that all the bone quality traits had a large variable coefficient, especially mechanical properties trait (20-70%), indicating that there was a large difference in bone health status among laying ducks. The phenotypic correlations showed a high correlation between weight and density, diameter and perimeter, breaking and toughness (r = 0.52-0.68). And then, we performed the Genome-wide association study (GWAS) to reveal the candidate genes of humerus quality in ducks. Seven candidate protein-coding genes were identified with perimeter trait, and 52 protein-coding genes were associated with toughness trait. We also analysed the candidate region and performed KEGG and GO analyse for 75 candidate genes. Furthermore, the expression analyse of the above candidate genes in different stage of humerus and different tissues were performed. Finally, AP2A2, SMAD3, SMNDC1, NFIA, EPHB2, PMEPA1, UNC5C, ESR1, VAV3, NFATC2 deserve further focus. The obtained results can contribute to new insight into bone quality and provide new genetic biomarkers for application in duck breeding programs.

Vitamin D and vitamin K status in postmenopausal women with normal and low bone mineral density

OBJECTIVES

Vitamin D and K are believed to promote bone health, but existing evidence is controversial. This study aimed to measure several metabolites of both vitamins by liquid chromatography tandem mass spectrometry (LC-MS/MS) in a cohort of postmenopausal women with low and normal bone mineral density (BMD).

METHODS

Vitamin metabolites (25-hydroxyvitamin D (25[OH]D), 24,25-dihydroxyvitamin D (24,25(OH)2D), phylloquinone (K1), menaquinone-4 (MK-4) and MK-7) were measured in 131 serum samples by LC-MS/MS. The vitamin D metabolite ratio (VMR) was calculated. Parathyroid hormone (PTH), type I procollagen-N-terminal-peptide (PINP) and C-terminal telopeptides of type I collagen (CTX-I) were measured by immunoassay. Dual X-ray absorptiometry was performed to identify participants with normal (T-score>-1) and low (T-score<-1) BMD.

RESULTS

Mean age was 58.2±8.5 years. BMD was normal in 68 and low in 63 women. Median (interquartile range) for 25(OH)D and total vitamin K concentrations were 53.5 (39.6-65.9) nmol/L and 1.33 (0.99-2.39) nmol/L. All vitamin metabolites were comparable in individuals with normal and low BMD. Furthermore, BMD and trabecular bone score were comparable in participants with adequate and inadequate vitamin status (at least one criterion was met: 25(OH)D <50 nmol/L, 24,25(OH)2D <3 nmol/L, VMR <4 %, total vitamin K <0.91 nmol/L). PTH, but not PINP or CTX-I, was inversely correlated with 25(OH)D, 24,25(OH)2D and VMR. Synergistic effects between vitamin D and K were not observed.

CONCLUSIONS

Vitamin D and K status is not related to BMD and trabecular bone quality in postmenopausal women. Inverse associations were only seen between vitamin D metabolites and PTH.

GREM2 inactivation increases trabecular bone mass in mice

Osteoporosis is a common skeletal disease affecting millions of individuals world-wide, with an increased risk of fracture, and a decreased quality of life. Despite its well-known consequences, the etiology of osteoporosis and optimal treatment methods are not fully understood. Human genetic studies have identified genetic variants within the FMN2/GREM2 locus to be associated with trabecular volumetric bone mineral density (vBMD) and vertebral and forearm fractures, but not with cortical bone parameters. GREM2 is a bone morphogenetic protein (BMP) antagonist. In this study, we employed Grem2-deficient mice to investigate whether GREM2 serves as the plausible causal gene for the fracture signal at the FMN2/GREM2 locus. We observed that Grem2 is moderately expressed in bone tissue and particularly in osteoblasts. Complete Grem2 gene deletion impacted mouse survival and body growth. Partial Grem2 inactivation in Grem2+/- female mice led to increased trabecular BMD of femur and increased trabecular bone mass in tibia due to increased trabecular thickness, with an unchanged cortical thickness, as compared with wildtype littermates. Furthermore, Grem2 inactivation stimulated osteoblast differentiation, as evidenced by higher alkaline phosphatase (Alp), osteocalcin (Bglap), and osterix (Sp7) mRNA expression after BMP-2 stimulation in calvarial osteoblasts and osteoblasts from the long bones of Grem2-/- mice compared to wildtype littermates. These findings suggest that GREM2 is a possible target for novel osteoporotic treatments, to increase trabecular bone mass and prevent osteoporotic fractures.

Association of serum trimethylamine N-oxide levels and bone mineral density in type 2 diabetes mellitus

PURPOSE

The relationship between trimethylamine N-oxide (TMAO) and bone mineral density (BMD) in type 2 diabetes mellitus (T2DM) is unclear. We explore the relationship between TMAO levels and BMD in T2DM.

METHODS

This is a cross-sectional study. 254 T2DM patients were enrolled and divided into three groups by TMAO tertiles, and the clinical data were collected. BMD was determined by dual-energy X-ray absorptiometry (DXA) and serum TMAO levels was determined by stable isotope dilution high-performance liquid chromatography tandem mass spectrometry (HPLC-MS/MS).

RESULTS

Patients in the highest tertile of TMAO levels (TMAO > 6.72 μmol/L) showed relatively low BMD and a higher number of fracture history, osteoporosis (OP) than those in the lower tertiles. Spearman correlation analysis showed that serum TMAO was negatively correlated with BMD of whole body (WB), lumbar spine (LS) and femoral neck (FN), while TMAO was positive correlated with osteoporotic fracture (p < 0.05). Logistic regression models showed that TMAO was an independent influencing factor of fracture history after adjusting for confounders in TMAO > 6.72 μmol/L group.

CONCLUSIONS

There is a significant linear correlation between TMAO levels and BMD in T2DM patients. Especially in TMAO > 6.72 μmol/L group, TMAO was negatively correlated with WB, LS, and FN BMD, and was positive correlated with osteoporotic fracture in T2DM patients. The findings suggest that elevated TMAO levels are associated with OP and osteoporotic fracture in T2DM patients.

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.

Multimodal cell atlas of the ageing human skeletal muscle

Muscle atrophy and functional decline (sarcopenia) are common manifestations of frailty and are critical contributors to morbidity and mortality in older people1. Deciphering the molecular mechanisms underlying sarcopenia has major implications for understanding human ageing2. Yet, progress has been slow, partly due to the difficulties of characterizing skeletal muscle niche heterogeneity (whereby myofibres are the most abundant) and obtaining well-characterized human samples3,4. Here we generate a single-cell/single-nucleus transcriptomic and chromatin accessibility map of human limb skeletal muscles encompassing over 387,000 cells/nuclei from individuals aged 15 to 99 years with distinct fitness and frailty levels. We describe how cell populations change during ageing, including the emergence of new populations in older people, and the cell-specific and multicellular network features (at the transcriptomic and epigenetic levels) associated with these changes. On the basis of cross-comparison with genetic data, we also identify key elements of chromatin architecture that mark susceptibility to sarcopenia. Our study provides a basis for identifying targets in the skeletal muscle that are amenable to medical, pharmacological and lifestyle interventions in late life.

Vitamin D and human health: evidence from Mendelian randomization studies

We summarized the current evidence on vitamin D and major health outcomes from Mendelian randomization (MR) studies. PubMed and Embase were searched for original MR studies on vitamin D in relation to any health outcome from inception to September 1, 2022. Nonlinear MR findings were excluded due to concerns about the validity of the statistical methods used. A meta-analysis was preformed to synthesize study-specific estimates after excluding overlapping samples, where applicable. The methodological quality of the included studies was evaluated according to the STROBE-MR checklist. A total of 133 MR publications were eligible for inclusion in the analyses. The causal association between vitamin D status and 275 individual outcomes was examined. Linear MR analyses showed genetically high 25-hydroxyvitamin D (25(OH)D) concentrations were associated with reduced risk of multiple sclerosis incidence and relapse, non-infectious uveitis and scleritis, psoriasis, femur fracture, leg fracture, amyotrophic lateral sclerosis, anorexia nervosa, delirium, heart failure, ovarian cancer, non-alcoholic fatty liver disease, dyslipidemia, and bacterial pneumonia, but increased risk of Behçet's disease, Graves' disease, kidney stone disease, fracture of radium/ulna, basal cell carcinoma, and overall cataracts. Stratified analyses showed that the inverse association between genetically predisposed 25(OH)D concentrations and multiple sclerosis risk was significant and consistent regardless of the genetic instruments GIs selected. However, the associations with most of the other outcomes were only pronounced when using genetic variants not limited to those in the vitamin D pathway as GIs. The methodological quality of the included MR studies was substantially heterogeneous. Current evidence from linear MR studies strongly supports a causal role of vitamin D in the development of multiple sclerosis. Suggestive support for a number of other health conditions could help prioritize conditions where vitamin D may be beneficial or harmful.

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.

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

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

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.

Causal effect of education on bone mineral density: A Mendelian randomization study

Education level may have some association with the incidence of osteoporosis, but it is elusive if this association is causal. This two-sample Mendelian randomization analysis focused on the causal effect of education level on femoral neck bone mineral density (FN-BMD), forearm BMD, lumbar spine BMD, and heel BMD. Twelve single nucleotide polymorphisms were used as instrumental variables. The results suggested that high education level was associated with improved FN-BMD (beta-estimate: 0.406, 95% confidence interval: 0.061 to 0.751, standard error: 0.176, P-value = .021). There were null association between education and other sites of bone mineral density. Our results found the causal effect of high education level on improved FN-BMD, and improved educational attainment may be beneficial to prevent osteoporosis.

Integrative high-throughput enhancer surveying and functional verification divulges a YY2-condensed regulatory axis conferring risk for osteoporosis

The precise roles of chromatin organization at osteoporosis risk loci remain largely elusive. Here, we combined chromatin interaction conformation (Hi-C) profiling and self-transcribing active regulatory region sequencing (STARR-seq) to qualify enhancer activities of prioritized osteoporosis-associated single-nucleotide polymorphisms (SNPs). We identified 319 SNPs with biased allelic enhancer activity effect (baaSNPs) that linked to hundreds of candidate target genes through chromatin interactions across 146 loci. Functional characterizations revealed active epigenetic enrichment for baaSNPs and prevailing osteoporosis-relevant regulatory roles for their chromatin interaction genes. Further motif enrichment and network mapping prioritized several putative, key transcription factors (TFs) controlling osteoporosis binding to baaSNPs. Specifically, we selected one top-ranked TF and deciphered that an intronic baaSNP (rs11202530) could allele-preferentially bind to YY2 to augment PAPSS2 expression through chromatin interactions and promote osteoblast differentiation. Our results underline the roles of TF-mediated enhancer-promoter contacts for osteoporosis, which may help to better understand the intricate molecular regulatory mechanisms underlying osteoporosis risk loci.

Multi-modal molecular determinants of clinically relevant osteoporosis subtypes

Due to a rapidly aging global population, osteoporosis and the associated risk of bone fractures have become a wide-spread public health problem. However, osteoporosis is very heterogeneous, and the existing standard diagnostic measure is not sufficient to accurately identify all patients at risk of osteoporotic fractures and to guide therapy. Here, we constructed the first prospective multi-omics atlas of the largest osteoporosis cohort to date (longitudinal data from 366 participants at three time points), and also implemented an explainable data-intensive analysis framework (DLSF: Deep Latent Space Fusion) for an omnigenic model based on a multi-modal approach that can capture the multi-modal molecular signatures (M3S) as explicit functional representations of hidden genotypes. Accordingly, through DLSF, we identified two subtypes of the osteoporosis population in Chinese individuals with corresponding molecular phenotypes, i.e., clinical intervention relevant subtypes (CISs), in which bone mineral density benefits response to calcium supplements in 2-year follow-up samples. Many snpGenes associated with these molecular phenotypes reveal diverse candidate biological mechanisms underlying osteoporosis, with xQTL preferences of osteoporosis and its subtypes indicating an omnigenic effect on different biological domains. Finally, these two subtypes were found to have different relevance to prior fracture and different fracture risk according to 4-year follow-up data. Thus, in clinical application, M3S could help us further develop improved diagnostic and treatment strategies for osteoporosis and identify a new composite index for fracture prediction, which were remarkably validated in an independent cohort (166 participants).

Bone fragility and osteoporosis in children and young adults

Osteoporosis is a metabolic bone disorder which increases fragility fracture risk. Elderly individuals, especially postmenopausal women, are particularly susceptible to osteoporosis. Although rare, osteoporosis in children and young adults is becoming increasingly evident, highlighting the need for timely diagnosis, management and follow-up. Early-onset osteoporosis is defined as the presence of a low BMD (Z-score of ≤ -2.0 in individuals aged < 20 years; T-score of ≤ -2.5 in those aged between 20 to 50 years) accompanied by a clinically significant fracture history, or the presence of low-energy vertebral compression fractures even in the absence of osteoporosis. Affected children and young adults should undergo a thorough diagnostic workup, including collection of clinical history, radiography, biochemical investigation and possibly bone biopsy. Once secondary factors and comorbidities are excluded, genetic testing should be considered to determine the possibility of an underlying monogenic cause. Defects in genes related to type I collagen biosynthesis are the commonest contributors of primary osteoporosis, followed by loss-of-function variants in genes encoding key regulatory proteins of canonical WNT signalling (specifically LRP5 and WNT1), the actin-binding plastin-3 protein (encoded by PLS3) resulting in X-linked osteoporosis, and the more recent sphingomyelin synthase 2 (encoded by SGMS2) which is critical for signal transduction affecting sphingomyelin metabolism. Despite these discoveries, genetic causes and underlying mechanisms in early-onset osteoporosis remain largely unknown, and if no causal gene is identified, early-onset osteoporosis is deemed idiopathic. This calls for further research to unravel the molecular mechanisms driving early-onset osteoporosis that consequently will aid in patient management and individualised targeted therapy.

Causal relationship between glycemic traits and bone mineral density in different age groups and skeletal sites: a Mendelian randomization analysis

INTRODUCTION

Previous research has confirmed that patients with type 2 diabetes mellitus tend to have higher bone mineral density (BMD), but it is unknown whether this pattern holds true for individuals without diabetes. This Mendelian randomization (MR) study aims to investigate the potential causal relationship between various glycemic trait (including fasting glucose, fasting insulin, 2-h postprandial glucose, and glycated hemoglobin) and BMD in non-diabetic individuals. The investigation focuses on different age groups (15-30, 30-45, 45-60, and 60 + years) and various skeletal sites (forearm, lumbar spine, and hip).

MATERIALS AND METHODS

We utilized genome-wide association study data from large population-based cohorts to identify robust instrumental variables for each glycemic traits parameter. Our primary analysis employed the inverse-variance weighted method, with sensitivity analyses conducted using MR-Egger, weighted median, MR-PRESSO, and multivariable MR methods to assess the robustness and potential horizontal pleiotropy of the study results.

RESULTS

Fasting insulin showed a negative modulating relationship on both lumbar spine and forearm. However, these associations were only nominally significant. No significant causal association was observed between blood glucose traits and BMD across the different age groups. The direction of fasting insulin's causal effects on BMD showed inconsistency between genders, with potentially decreased BMD in women with high fasting insulin levels and an increasing trend in BMD in men.

CONCLUSIONS

In the non-diabetic population, currently available evidence does not support a causal relationship between glycemic traits and BMD. However, further investigation is warranted considering the observed gender differences.

Fracture risk prediction in old Chinese people-a narrative review

With aging, the burden of osteoporotic fracture (OF) increases substantially, while China is expected to carry the greatest part in the future. The risk of fracture varies greatly across racial groups and geographic regions, and systematically organized evidence on the potential predictors for fracture risk is needed for Chinese. This review briefly introduces the epidemiology of OF and expands on the predictors and predictive tools for the risk of OF, as well as the challenges for their potential translation in the old Chinese population. There are regional differences of fracture incidence among China. The fracture incidences in Hong Kong and Taiwan have decreased in recent years, while it is still increasing in mainland China. Although the application of dual-energy X-ray absorptiometry (DXA) is limited among old Chinese in the mainland, bone mineral density (BMD) by DXA has a predictive value similar to that worldwide. Other non-DXA modalities, especially heel QUS, are helpful in assessing bone health. The fracture risk assessment tool (FRAX) has a good discrimination ability for OFs, especially the FRAX with BMD. And some clinical factors have added value to FRAX, which has been verified in old Chinese. In addition, although the application of the osteoporosis self-assessment tool for Asians (OSTA) in Chinese needs further validation, it may help identify high-risk populations in areas with limited resources. Moreover, the translation use of the muscle quality and genetic or serum biomarkers in fracture prediction needs further works. More applicable and targeted fracture risk predictors and tools are still needed for the old Chinese population.

Targeting Longevity Gene SLC13A5: A Novel Approach to Prevent Age-Related Bone Fragility and Osteoporosis

Reduced expression of the plasma membrane citrate transporter SLC13A5, also known as INDY, has been linked to increased longevity and mitigated age-related cardiovascular and metabolic diseases. Citrate, a vital component of the tricarboxylic acid cycle, constitutes 1-5% of bone weight, binding to mineral apatite surfaces. Our previous research highlighted osteoblasts' specialized metabolic pathway facilitated by SLC13A5 regulating citrate uptake, production, and deposition within bones. Disrupting this pathway impairs bone mineralization in young mice. New Mendelian randomization analysis using UK Biobank data indicated that SNPs linked to reduced SLC13A5 function lowered osteoporosis risk. Comparative studies of young (10 weeks) and middle-aged (52 weeks) osteocalcin-cre-driven osteoblast-specific Slc13a5 knockout mice (Slc13a5cKO) showed a sexual dimorphism: while middle-aged females exhibited improved elasticity, middle-aged males demonstrated enhanced bone strength due to reduced SLC13A5 function. These findings suggest reduced SLC13A5 function could attenuate age-related bone fragility, advocating for SLC13A5 inhibition as a potential osteoporosis treatment.

High-throughput screening of glucocorticoid-induced enhancer activity reveals mechanisms of stress-related psychiatric disorders

Exposure to stressful life events increases the risk for psychiatric disorders. Mechanistic insight into the genetic factors moderating the impact of stress can increase our understanding of disease processes. Here, we test 3,662 single nucleotide polymorphisms (SNPs) from preselected expression quantitative trait loci in massively parallel reporter assays to identify genetic variants that modulate the activity of regulatory elements sensitive to glucocorticoids, important mediators of the stress response. Of the tested SNP sequences, 547 were located in glucocorticoid-responsive regulatory elements of which 233 showed allele-dependent activity. Transcripts regulated by these functional variants were enriched for those differentially expressed in psychiatric disorders in the postmortem brain. Phenome-wide Mendelian randomization analysis in 4,439 phenotypes revealed potentially causal associations specifically in neurobehavioral traits, including major depression and other psychiatric disorders. Finally, a functional gene score derived from these variants was significantly associated with differences in the physiological stress response, suggesting that these variants may alter disease risk by moderating the individual set point of the stress response.

Assessing the Casual Association between Sex Hormone Levels and Fracture Risk: A Two-Sample Mendelian Randomization Study

OBJECTIVE

Prior observational studies have reported that levels of sex hormones constitute a risk factor for the fracture. The aim of this study was to ascertain whether there is a causal relationship between the levels of sex hormones and the risk of fracture through Mendelian randomization (MR).

METHODS

Single-nucleotide polymorphisms (SNPs) associated with two indicators of sex hormone levels, circulating sex hormone-binding globulin (SHBG) and bioavailable testosterone levels, as exposures were selected from a large genome-wide association study (GWAS) from UK Biobank. The summary statistics for 11 different types of fracture as outcomes from the FinnGen consortium. This study employed the two-sample MR approach. For the main analysis, the inverse-variance-weighted (IVW) method was utilized. To assess the heterogeneity of MR results, the IVW method and MR-Egger method were utilized. To evaluate potential pleiotropy, MR-Egger regression was conducted. Additionally, a leave-one-SNP-out test was performed to assess the robustness of MR results to the exclusion of any individual SNP.

RESULTS

The MR analyses demonstrated a conspicuous impact of SHBG on the risk of pathological fracture with osteoporosis (OP). We found that an increase of one standard deviation (SD) in SHBG correspondingly increased the risk of pathological fracture with OP [odds ratio (OR) 2.42, 95% confidence interval (CI), 1.52-3.85; p = 1.93 × 10-4 ]. The bioavailable testosterone showed the negative casual genetic associations with fractures of foot and forearm. An increase of one SD in the genetically predetermined bioavailable testosterone was associated with a reduction of 37% in the risk of fracture of foot (OR 0.63, 95% Cl 0.49 to 0.81; p = 3.37 × 10-4 ), as well as a 39% decrease in the risk of fracture of forearm (OR 0.61, 95% Cl 0.50 to 0.76; p = 5.40 × 10-6 ).

CONCLUSIONS

Our study confirms that individuals experiencing elevated SHBG concentrations showed a major causal effect on pathological fracture with OP. High bioavailable testosterone levels play an important role in preventing the fractures of foot and forearm. Although increasing bioavailable testosterone and decreasing SHBG levels had no casual effect on most fractures in the general population, they are likely to have the most clinically relevant effect on certain fracture risk reduction.

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.

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.

An atlas of genetic determinants of forearm fracture

Osteoporotic fracture is among the most common and costly of diseases. While reasonably heritable, its genetic determinants have remained elusive. Forearm fractures are the most common clinically recognized osteoporotic fractures with a relatively high heritability. To establish an atlas of the genetic determinants of forearm fractures, we performed genome-wide association analyses including 100,026 forearm fracture cases. We identified 43 loci, including 26 new fracture loci. Although most fracture loci associated with bone mineral density, we also identified loci that primarily regulate bone quality parameters. Functional studies of one such locus, at TAC4, revealed that Tac4-/- mice have reduced mechanical bone strength. The strongest forearm fracture signal, at WNT16, displayed remarkable bone-site-specificity with no association with hip fractures. Tall stature and low body mass index were identified as new causal risk factors for fractures. The insights from this atlas may improve fracture prediction and enable therapeutic development to prevent fractures.

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.

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.

Prevention of Hip Fractures: Trade-off between Minor Benefits to Individuals and Large Benefits to the Community

Goeffrey Rose postulated that a population-based measure bringing a small benefit to each individual can yield large benefits to the community. We aimed to test this axiom by quantifying the relationship between change in bone mineral density (BMD) and hip fracture incidence between two prospective cohorts separated by ~10 years. In this prospective population-based Dubbo Osteoporosis Epidemiology Study (DOES), the participants aged 60+ were recruited in two waves: the initial cohort (1311 women, 842 men) in 1989 to 1992 and the second cohort (974 women, 544 men) in 1999 to 2001. The incident hip fracture was radiologically ascertained. Femoral neck BMD was measured biannually. Multivariable-adjusted Cox's proportional hazards models were adjusted for the predefined covariates such as age, BMI, lifestyle factors, falls, and prior fracture. Compared with the initial cohort, the second cohort had a higher femoral neck BMD by ~0.04 g/cm2 in women and 0.03 g/cm2 in men. However, the prevalence of osteoporosis in the second cohort was halved (prevalence ratio 0.51, 95% CI 0.36 to 0.73 in women; 0.45, 0.24 to 0.84 in men), and its hip fracture incidence was significantly reduced (hazard ratio 0.54, 95% CI, 0.38 to 0.78 in women; 0.39, 0.19 to 0.80 in men). Sensitivity analyses indicated that the "effect" was unlikely due to unmeasured confounders. These findings suggest that a population-wide strategy aimed at enhancing BMD across the entire population could lead to a substantial decrease in the incidence of hip fractures. © 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).

Native lamin A/C proteomes and novel partners from heart and skeletal muscle in a mouse chronic inflammation model of human frailty

Clinical frailty affects ∼10% of people over age 65 and is studied in a chronically inflamed (Interleukin-10 knockout; "IL10-KO") mouse model. Frailty phenotypes overlap the spectrum of diseases ("laminopathies") caused by mutations in LMNA. LMNA encodes nuclear intermediate filament proteins lamin A and lamin C ("lamin A/C"), important for tissue-specific signaling, metabolism and chromatin regulation. We hypothesized that wildtype lamin A/C associations with tissue-specific partners are perturbed by chronic inflammation, potentially contributing to dysfunction in frailty. To test this idea we immunoprecipitated native lamin A/C and associated proteins from skeletal muscle, hearts and brains of old (21-22 months) IL10-KO versus control C57Bl/6 female mice, and labeled with Tandem Mass Tags for identification and quantitation by mass spectrometry. We identified 502 candidate lamin-binding proteins from skeletal muscle, and 340 from heart, including 62 proteins identified in both tissues. Candidates included frailty phenotype-relevant proteins Perm1 and Fam210a, and nuclear membrane protein Tmem38a, required for muscle-specific genome organization. These and most other candidates were unaffected by IL10-KO, but still important as potential lamin A/C-binding proteins in native heart or muscle. A subset of candidates (21 in skeletal muscle, 30 in heart) showed significantly different lamin A/C-association in an IL10-KO tissue (p < 0.05), including AldoA and Gins3 affected in heart, and Lmcd1 and Fabp4 affected in skeletal muscle. To screen for binding, eleven candidates plus prelamin A and emerin controls were arrayed as synthetic 20-mer peptides (7-residue stagger) and incubated with recombinant purified lamin A "tail" residues 385-646 under relatively stringent conditions. We detected strong lamin A binding to peptides solvent exposed in Lmcd1, AldoA, Perm1, and Tmem38a, and plausible binding to Csrp3 (muscle LIM protein). These results validated both proteomes as sources for native lamin A/C-binding proteins in heart and muscle, identified four candidate genes for Emery-Dreifuss muscular dystrophy (CSRP3, LMCD1, ALDOA, and PERM1), support a lamin A-interactive molecular role for Tmem38A, and supported the hypothesis that lamin A/C interactions with at least two partners (AldoA in heart, transcription factor Lmcd1 in muscle) are altered in the IL10-KO model of frailty.

Review of mendelian randomization studies on age at natural menopause

Menopause marks the end of the reproductive phase of life. Based on epidemiological studies, abnormal age at natural menopause (ANM) is thought to contribute to a number of adverse outcomes, such as osteoporosis, cardiovascular disease, and cancer. However, the causality of these associations remains unclear. A powerful epidemiological method known as Mendelian randomization (MR) can be used to clarify the causality between ANM and other diseases or traits. The present review describes MR studies that included ANM as an exposure, outcome and mediator. The findings of MR analyses on ANM have revealed that higher body mass index, poor educational level, early age at menarche, early age at first live birth, early age at first sexual intercourse, and autoimmune thyroid disease appear to be involved in early ANM etiology. The etiology of late ANM appears to be influenced by higher free thyroxine 4 and methylene tetrahydrofolate reductase gene mutations. Furthermore, early ANM has been found to be causally associated with an increased risk of osteoporosis, fracture, type 2 diabetes mellitus, glycosylated hemoglobin, and the homeostasis model of insulin resistance level. In addition, late ANM has been found to be causally associated with an increased systolic blood pressure, higher risk of breast cancer, endometrial cancer, endometrioid ovarian carcinoma, lung cancer, longevity, airflow obstruction, and lower risk of Parkinson's disease. ANM is also a mediator for breast cancer caused by birth weight and childhood body size. However, due to the different instrumental variables used, some results of studies are inconsistent. Future studies with more valid genetic variants are needed for traits with discrepancies between MRs or between MR and other types of epidemiological studies.

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.

Long-term and sequential treatment for osteoporosis

Osteoporosis is a skeletal disorder that causes impairment of bone structure and strength, leading to a progressively increased risk of fragility fractures. The global prevalence of osteoporosis is increasing in the ageing population. Owing to the chronic character of osteoporosis, years or even decades of preventive measures or therapy are required. The long-term use of bone-specific pharmacological treatment options, including antiresorptive and/or osteoanabolic approaches, has raised concerns around adverse effects or potential rebound phenomena after treatment discontinuation. Imaging options, risk scores and the assessment of bone turnover during initiation and monitoring of such therapies could help to inform individualized treatment strategies. Combination therapies are currently used less often than 'sequential' treatments. However, all patients with osteoporosis, including those with secondary and rare causes of osteoporosis, as well as specific patient populations (for example, young adults, men and pregnant women) require new approaches for long-term therapy and disease monitoring. New pathophysiological aspects of bone metabolism might therefore help to inform and revolutionize the diagnosis and treatment of osteoporosis.

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.