Endochondral ossification pathway genes and postmenopausal osteoporosis: Association and specific allele related serum bone sialoprotein levels in Han Chinese

IL-16 Mediates the Effect of Circulating Metabolites on Postmenopausal Osteoporosis: A Two-Step, Multivariable Mendelian Randomization Study

Objectives: This study aimed to explore the relationship between circulating metabolites and postmenopausal osteoporosis (PMOP) and to assess the mediating role of inflammatory factors. Methods: Utilizing summary-level data from genome-wide association studies (GWAS) and employing a Mendelian Randomization approach, a two-sample MR analysis was conducted to assess the relationship between circulating metabolites and PMOP. Additionally, a two-step MR was used to quantify the mediating impact of inflammatory factors on the effect of circulating metabolites on PMOP. Results: The results revealed a significant association between certain metabolites and the risk of PMOP, notably the ratio of free cholesterol to total lipids in very large VLDL particles (OR: 1.399, 95% CI: 1.002-1.954, p = 0.048) and IL-16 (OR: 0.773, 95% CI: 0.608-0.983, p = 0.036). IL-16 was found to partially mediate the impact of circulating metabolites on PMOP, with a mediation effect of 10.4%. Conclusion: This study underscores the crucial role of circulating metabolites and inflammatory factors in PMOP pathogenesis. A causal relationship between circulating metabolites and PMOP was established, with IL-16 mediating some effects. These findings hold promise for clinical applications in early detection, personalized medicine, and the identification of therapeutic targets for PMOP.

The role of bone sialoprotein in bone healing

Bone sialoprotein (BSP) is a multi-functional extracellular matrix (ECM) protein associated with mineralized tissues, particularly bone and cementum. The amino acid sequence of BSP includes three evolutionarily conserved sequences which contribute to functions of the protein: an N-terminal collagen-binding domain, polyglutamic acid (polyE) sequences involved in hydroxyapatite nucleation and crystal growth, and a C-terminal arginine-glycine-aspartic acid (RGD) integrin-binding domain. BSP promotes attachment and differentiation of osteogenic and osteoclastic cells. Genetic ablation of BSP in mice results in skeletal and dental developmental defects and impaired bone healing in both appendicular bone and alveolar bone of the jaw. Several studies demonstrated positive effects of BSP on bone healing in rodent models, though other experiments show negligible results. Native (harvested from rat bones) BSP cross-linked to collagen induced slight improvements in calvarial bone healing in rats. Recombinant BSP and collagen delivered in a polylactide (PLA) cylinder improved bone defect healing in rat femurs. Both native and recombinant BSP delivered in a collagen gel improved alveolar bone healing in wild-type and BSP-deficient mice. These advances suggest BSP is a new player in bone healing that has potential to be an alternative or complimentary to other bioactive factors. Future studies are necessary to understand mechanisms of how BSP influences bone healing and optimize delivery and dose in different types of bone defects and injuries.

The bone-protective benefits of amino-conjugated calcium in an ovariectomized (OVX) rat model

Low bone density, fragility, and microarchitectural disintegration are the symptoms of osteoporosis. An imbalance between bone growth and resorption can lead to osteoporosis. This study evaluated the effects of amino-calcium (AC) on bone protection in ovariectomized control group (NC) rats. Amino-calcium (AC) was characterized using Fourier-transform infrared spectroscopy (FT-IR), energy-dispersive X-ray spectroscopy (EDS), and nuclear magnetic resonance spectroscopy analyses (NMR). After determining the biocompatibility of amino-calcium (AC) with MC3T3-E1 cells, alkaline phosphatase staining revealed significant changes on day 7. Three of the four groups underwent ovariectomy, whereas one group received a placebo. On micro-computed tomography, in vivo, data showed increased bone volume fraction in the femoral head and shaft areas in the amino-calcium (AC) group. Hematoxylin and eosin staining showed a bone mass and architectural protection in the amino-calcium (AC) group compared with the calcium carbonate and OVX control group. RNA sequencing analysis revealed high expression of osteogenesis-related genes in MC3T3-E1 cells. RNA sequencing revealed a significant fold change in the expression of integrin-binding sialoprotein (IBSP), bone gamma-carboxyglutamate proteins 1 and 2(BGLAP1 and BGLAP2), and periostin (POSTN). The study concluded that supplementing the OVX rats with calcium enhanced bone protection.

Mechanism of Erzhiwan in treating osteoporosis based on molecular docking technology and molecular dynamics simulation

This experiment was a network pharmacology research based on the theoretical system of traditional Chinese medicine. TCMSP database, PubChem database, RCSB database, and SwissTargetPrediction database were used to study the effective chemical constituents of Ligustri lucidi Fructus and Ecliptae Herba in Erzhiwan, a traditional prescription for nourishing the liver and kidney. Then Genecards database, OMIM database, OMIM Gene Map, and Metascape database were used to study the therapeutic targets of osteoporosis. At last, Cytoscape 3.6.0 software, its built-in Bisogenet and CytoNCA, AutoDockTools-1.5.6 software, PYMOL-2.2.0 software, and Gromacs software, by drawing the relationship diagram between chemical components and disease targets, PPI network of disease, semi-flexible molecular docking technology, evaluation and analysis of enrichment pathway, and molecular dynamics simulation, were used to study the therapeutic mechanism of Erzhiwan on osteoporosis. It is found that the intervention and regulation of Erzhiwan on osteoporosis were mainly realized through multiple targets of active ingredients and multiple pathways, which provided support for the continued development of Erzhiwan.

Development and validation of a gene signature predicting the risk of postmenopausal osteoporosis

AIMS

We aimed to develop a gene signature that predicts the occurrence of postmenopausal osteoporosis (PMOP) by studying its genetic mechanism.

METHODS

Five datasets were obtained from the Gene Expression Omnibus database. Unsupervised consensus cluster analysis was used to determine new PMOP subtypes. To determine the central genes and the core modules related to PMOP, the weighted gene co-expression network analysis (WCGNA) was applied. Gene Ontology enrichment analysis was used to explore the biological processes underlying key genes. Logistic regression univariate analysis was used to screen for statistically significant variables. Two algorithms were used to select important PMOP-related genes. A logistic regression model was used to construct the PMOP-related gene profile. The receiver operating characteristic area under the curve, Harrell's concordance index, a calibration chart, and decision curve analysis were used to characterize PMOP-related genes. Then, quantitative real-time polymerase chain reaction (qRT-PCR) was used to verify the expression of the PMOP-related genes in the gene signature.

RESULTS

We identified three PMOP-related subtypes and four core modules. The muscle system process, muscle contraction, and actin filament-based movement were more active in the hub genes. We obtained five feature genes related to PMOP. Our analysis verified that the gene signature had good predictive power and applicability. The outcomes of the GSE56815 cohort were found to be consistent with the results of the earlier studies. qRT-PCR results showed that RAB2A and FYCO1 were amplified in clinical samples.

CONCLUSION

The PMOP-related gene signature we developed and verified can accurately predict the risk of PMOP in patients. These results can elucidate the molecular mechanism of RAB2A and FYCO1 underlying PMOP, and yield new and improved treatment strategies, ultimately helping PMOP monitoring.Cite this article: Bone Joint Res 2022;11(8):548-560.

NFKB2 polymorphisms associate with the risk of developing rheumatoid arthritis and response to TNF inhibitors: Results from the REPAIR consortium

This study sought to evaluate the association of 28 single nucleotide polymorphisms (SNPs) within NFKB and inflammasome pathway genes with the risk of rheumatoid arthritis (RA) and response to TNF inhibitors (TNFi). We conducted a case-control study in a European population of 1194 RA patients and 1328 healthy controls. The association of potentially interesting markers was validated with data from the DANBIO (695 RA patients and 978 healthy controls) and DREAM (882 RA patients) registries. The meta-analysis of our data with those from the DANBIO registry confirmed that anti-citrullinated protein antibodies (ACPA)-positive subjects carrying the NFKB2_rs11574851T allele had a significantly increased risk of developing RA (PMeta_ACPA + = 0.0006) whereas no significant effect was found in ACPA-negative individuals (PMeta_ACPA- = 0.35). An ACPA-stratified haplotype analysis including both cohorts (n = 4210) confirmed that ACPA-positive subjects carrying the NFKB2_TT haplotype had an increased risk of RA (OR = 1.39, P = 0.0042) whereas no effect was found in ACPA-negative subjects (OR = 1.04, P = 0.82). The meta-analysis of our data with those from the DANBIO and DREAM registries also revealed a suggestive association of the NFKB2_rs1056890 SNP with larger changes in DAS28 (OR = 1.18, P = 0.007). Functional experiments showed that peripheral blood mononuclear cells from carriers of the NFKB2_rs1005044C allele (in LD with the rs1056890, r2 = 1.00) showed increased production of IL10 after stimulation with LPS (P = 0.0026). These results provide first evidence of a role of the NFKB2 locus in modulating the risk of RA in an ACPA-dependent manner and suggest its implication in determining the response to TNFi. Additional studies are now warranted to further validate these findings.

Combining information from multiple bone turnover markers as diagnostic indices for osteoporosis using support vector machines

CONTEXT

Osteoporosis (OP) is a progressive systemic bone disease. Dual-energy X-ray absorptiometry (DXA) is routinely employed and is considered the gold standard method for the diagnosis of OP.

OBJECTIVE

We aimed to investigate the potential use of combined information from multiple bone turnover markers (BTMs) as a clinical diagnostic tool for OP.

MATERIALS AND METHODS

A total of 9053 Chinese postmenopausal women (2464 primary OP patients and 6589 healthy controls) were recruited. Serum levels of six common BTMs, including BAP, BSP, CTX, OPG, OST and sRANKL were assayed. Models based on support vector machine (SVM) were constructed to explore the efficiency of different combinations of multiple BTMs for OP diagnosis.

RESULTS

Increasing the number of BTMs used in generating the models increased the predictive power of the SVM models for determining the disease status of study subjects. The highest kappa coefficient for the model with one BTM (BAP) compared to DXA was 0.7783. The full model incorporating all six BTMs resulted in a high kappa coefficient of 0.9786.

CONCLUSION

Our findings showed that although single BTMs were not sufficient for OP diagnosis, appropriate combinations of multiple BTMs incorporated into the SVM models showed almost perfect agreement with the DXA.

Identification of aberrantly expressed long non-coding RNAs in postmenopausal osteoporosis

Postmenopausal osteoporosis (PMOP) is a common skeletal disorder in postmenopausal women. The present study aimed to identify the key long non‑coding RNAs (lncRNAs) in PMOP through RNA sequencing. RNA sequencing was performed to obtain the expression profile of lncRNAs and mRNAs in blood samples of patients with PMOP and normal controls (NCs). Following the identification of differentially expressed mRNAs (DEmRNAs) and differentially expressed lncRNAs (DElncRNAs), the DElncRNA-DEmRNA co‑expression network was constructed. A search was performed for the DEGs transcribed within a 100‑kb window upstream or downstream of DElncRNAs, which served as nearby DEmRNAs of DElncRNAs. Functional annotation of the DEmRNAs co‑expressed with DElncRNAs was performed. The GSE56815 dataset was used to verify the expression of selected DEmRNAs and DElncRNAs. Three blood samples from patients with PMOP and two blood samples from NCs were used for RNA sequencing. Compared with the NC group, a total of 185 DEmRNAs and 51 DElncRNAs were obtained in PMOP. A total of 3,057 co‑expression DElncRNA‑DEmRNA pairs and 97 DElncRNA‑nearby DEmRNA pairs were obtained. Six DEmRNAs [diacylglycerol O‑acyltransferase 2, potassium voltage‑gated channel subfamily S member 1, peptidase inhibitor 3, secretory leukocyte peptidase inhibitor, galectin‑related protein and alkaline phosphatase, liver/bone/kidney (ALPL)] were nearby co‑expressed genes of four DElncRNAs, including LOC105376834, LOC101929866, LOC105374771 and LOC100506113. Three PMOP-associated DEmRNAs, including ALPL, suppressor of cytokine signaling 3 and adrenomedullin, were co‑expressed with the hub DElncRNAs (LINC00963, LOC105378415, LOC105377067, HCG27, LOC101928143 and LINC01094) of the positively and negatively co‑expressed DElncRNA‑DEmRNA interaction network. The expression of selected DEmRNAs and DElncRNAs was consistent with the RNA‑sequencing results. In conclusion, the present study identified the key DEmRNAs and DElncRNAs in PMOP, which may provide clues for understanding the mechanism and developing novel biomarkers for PMOP.