ATF3 as a potential diagnostic marker of early-stage osteoarthritis and its correlation with immune infiltration through bioinformatics analysis

A comparative metabolomic analysis reveals the metabolic variations among cartilage of Kashin-Beck disease and osteoarthritis

Aims

The metabolic variations between the cartilage of osteoarthritis (OA) and Kashin-Beck disease (KBD) remain largely unknown. Our study aimed to address this by conducting a comparative analysis of the metabolic profiles present in the cartilage of KBD and OA.

Methods

Cartilage samples from patients with KBD (n = 10) and patients with OA (n = 10) were collected during total knee arthroplasty surgery. An untargeted metabolomics approach using liquid chromatography coupled with mass spectrometry (LC-MS) was conducted to investigate the metabolomics profiles of KBD and OA. LC-MS raw data files were converted into mzXML format and then processed by the XCMS, CAMERA, and metaX toolbox implemented with R software. The online Kyoto Encyclopedia of Genes and Genomes (KEGG) database was used to annotate the metabolites by matching the exact molecular mass data of samples with those from the database.

Results

A total of 807 ion features were identified for KBD and OA, including 577 positive (240 for upregulated and 337 for downregulated) and 230 negative (107 for upregulated and 123 for downregulated) ions. After annotation, LC-MS identified significant expressions of ten upregulated and eight downregulated second-level metabolites, and 183 upregulated and 162 downregulated first-level metabolites between KBD and OA. We identified differentially expressed second-level metabolites that are highly associated with cartilage damage, including dimethyl sulfoxide, uric acid, and betaine. These metabolites exist in sulphur metabolism, purine metabolism, and glycine, serine, and threonine metabolism.

Conclusion

This comprehensive comparative analysis of metabolism in OA and KBD cartilage provides new evidence of differences in the pathogenetic mechanisms underlying cartilage damage in these two conditions.

ATF3 affects osteogenic differentiation in inflammatory hPDLSCs by mediating ferroptosis via regulating the Nrf2/HO-1 signaling pathway

Activating transcription factor 3 (ATF3) has been identified as a regulator associated with osteoblast differentiation. However, the effects of ATF3 on the osteogenic differentiation and proliferation of human periodontal stem cells (hPDLSCs) in periodontitis have not been reported. With the purpose of establishing an in vitro model of periodontitis, hPDLSCs were challenged with lipopolysaccharide (LPS). The Cell Counting Kit-8 assay was applied to assess cell viability, while reverse transcription-quantitative PCR and western blotting were employed to detect ATF3 expression. Inflammatory release was assessed using ELISA, together with western blotting. Lipid peroxidation was explored using the C11 BODIPY 581/591 probe, biochemical kits, thiobarbituric acid reactive substances (TBARS) assay and DCFH-DA staining. Iron and Fe2+ levels, and the expression levels of ferroptosis-related proteins were measured using corresponding kits and western blotting. Osteogenic differentiative capability was evaluated using alkaline phosphatase staining, Alizarin red staining and western blotting. The expression levels of proteins associated with Nrf2/HO-1 signaling were identified using western blotting. The results indicated that ATF3 expression was upregulated in LPS-induced hPDLSCs. The knockdown of ATF3 alleviated the LPS-induced inflammatory response in hPDLSCs, together with increased levels of TNF-α, IL-6, IL-1β, Cox-2 and iNOS, and decreased levels of IL-10. ATF3 silencing also led to lower TBARS production rate, and reduced levels of reactive oxygen species, iron, Fe2+, ACSL4 and TFR1, whereas it elevated the levels of SLC7A11 and GPX4. In addition, ATF3 silencing promoted hPDLSC mineralization and cell differentiation, and elevated the levels of OCN2, RUNX2 and BMP2. Additionally, ATF3 depletion upregulated the expression levels of proteins related with Nrf2/HO-1 signaling. The Nrf2 inhibitor ML385 partially counteracted the effects of ATF3 interference on the LPS-challenged inflammatory response, lipid peroxidation, ferroptosis as well as osteogenic differentiative capability in hPDLSCs. In summary, the results revealed that ATF3 silencing suppressed inflammation and ferroptosis, while it improved osteogenic differentiation in LPS-induced hPDLSCs by regulating Nrf2/HO-1 signaling, which may provide promising therapeutic targets for the treatment of periodontitis.

Cartilage oligomeric matrix protein as a potential biomarker for knee osteoarthritis

Aims

This study aimed to determine the expression and clinical significance of a cartilage protein, cartilage oligomeric matrix protein (COMP), in knee osteoarthritis (OA) patients.

Methods

A total of 270 knee OA patients and 93 healthy controls were recruited. COMP messenger RNA (mRNA) and protein levels in serum, synovial fluid, synovial tissue, and fibroblast-like synoviocytes (FLSs) of knee OA patients were determined using enzyme-linked immunosorbent assay, real-time polymerase chain reaction, and immunohistochemistry.

Results

COMP protein levels were significantly elevated in serum and synovial fluid of knee OA patients, especially those in the advanced stages of the disease. Serum COMP was significantly correlated with radiological severity as well as measures of body composition, physical performance, knee pain, and disability. Receiver operating characteristic curve analysis unveiled a diagnostic value of serum COMP as a biomarker of knee OA (41.64 ng/ml, area under the curve (AUC) = 1.00), with a sensitivity of 99.6% and a specificity of 100.0%. Further analysis uncovered that COMP mRNA expression was markedly upregulated in the inflamed synovium of knee OA, consistent with immunohistochemical staining revealing localization of COMP protein in the lining and sub-lining layers of knee OA inflamed synovium. Most notably, relative COMP mRNA expression in knee OA synovium was positively associated with its protein levels in serum and synovial fluid of knee OA patients. In human knee OA FLSs activated with tumour necrosis factor-alpha, COMP mRNA expression was considerably up-regulated in a time-dependent manner.

Conclusion

All results indicate that COMP might serve as a supportive diagnostic marker for knee OA in conjunction with the standard diagnostic methods.

Identification and Verification of Endoplasmic Reticulum Stress-Related Genes as Novel Signatures for Osteoarthritis Diagnosis and Therapy: A Bioinformatics Analysis-Oriented Pilot Study

BACKGROUND AND PURPOSE

Endoplasmic reticulum stress (ERS) has been reported to be closely associated with the development of osteoarthritis (OA), but the underlying mechanisms are not fully delineated. The present study was designed to investigate the involvement of ERS-related genes in regulating OA progression.

METHODS

The expression profiles of OA patients and normal people were downloaded from the gene expression omnibus (GEO) database. The differentially expressed genes (DEGs) in datasets GSE55457 and GSE55235 were screened and identified by R software with the construction of the protein-protein interaction (PPI) networks. Through the STRING and Venn diagram analysis, hub ERS-related genes were obtained. Gene ontology (GO) and kyoto encyclopedia of genes and genomes (KEGG) enrichment analyses were performed. Biomarkers with high diagnostic values of osteoarthritis (OA) were studied. The hematoxylin and eosin (H&E) staining and micro-CT were applied to evaluate the establishment of the OA model. The expression levels of biomarkers were validated with the use of reverse transcription‑quantitative polymerase chain reaction (RT-qPCR) and western blot. Finally, we evaluated the correlations of hub ERS-related genes with the immune infiltration cells via the CIBERSORT algorithm.

RESULTS

A total of 60 downregulated and 52 upregulated DEGs were identified, and the following GO and KEGG pathway analyses verified that those DEGs were mainly enriched in biological process (BP), cellular component (CC), molecular function (MF), and inflammation-associated signal pathways. Interestingly, among all the DEGs, six ER stress-associated genes, including activating transcription factor 3 (ATF3), DEAD-Box Helicase 3 X-Linked (DDX3X), AP-1 transcription factor subunit (JUN), eukaryotic initiation factor 4 (EIF4A1), KDEL endoplasmic reticulum protein retention receptor 3 (KDELR3), and vascular endothelial growth factor A (VEGFA), were found to be closely associated with OA progression, and the following RT-qPCR and Western Blot analysis confirmed that DDX3X, JUN, and VEGFA were upregulated, whereas KDELR3, EIF4A1, and ATF3 were downregulated in OA rats tissues compared to the normal tissues, which were in accordance with our bioinformatics findings. Furthermore, our receiver operating characteristic (ROC) curve analysis verified that the above six ER stress-associated genes could be used as ideal biomarkers for OA diagnosis and those genes also potentially regulated immune responses by influencing the biological functions of mast cells and macrophages.

CONCLUSION

Collectively, the present study firstly identified six ER stress-associated genes (ATF3, DDX3X, JUN, EIF4A1, KDELR3, and VEGFA) that may play critical role in regulating the progression of OA.

Lipolysis products from triglyceride-rich lipoproteins induce stress protein ATF3 in osteoblasts

High fat diets overwhelm the physiological mechanisms for absorption, storage, and utilization of triglycerides (TG); consequently TG, TG-rich lipoproteins (TGRL), and TGRL remnants accumulate, circulate systemically, producing dyslipidemia. This associates with, or is causative for increased atherosclerotic cardiovascular risk, ischemic stroke, fatty liver disease, and pancreatitis. TGRL hydrolysis by endothelial surface-bound lipoprotein lipase (LPL) generates metabolites like free fatty acids which have proinflammatory properties. While osteoblasts utilize fatty acids as an energy source, dyslipidemia is associated with negative effects on the skeleton. In this study we investigated the effects of TGRL lipolysis products (TGRL-LP) on expression of a stress responsive transcription factor, termed activating transcription factor 3 (ATF3), reactive oxygen species (ROS), ATF3 target genes, and angiopoietin-like 4 (Angptl4) in osteoblasts. As ATF3 negatively associates with osteoblast differentiation, we also investigated the skeletal effects of global ATF3 deletion in mice. TGRL-LP increased expression of Atf3, proinflammatory proteins Ptgs2 and IL-6, and induced ROS in MC3T3-E1 osteoblastic cells. Angptl4 is an endogenous inhibitor of LPL which was transcriptionally induced by TGRL-LP, while recombinant Angptl4 prevented TG-driven Atf3 induction. Atf3 global knockout male mice demonstrated increased trabecular and cortical microarchitectural parameters. In summary, we find that TGRL-LP induce osteoblastic cell stress as evidenced by expression of ATF3, which may contribute to the negative impact of dyslipidemia in the skeleton. Further, concomitant induction of Angptl4 in osteoblasts might play a protective role by reducing local lipolysis.

Insights into the Roles of Natural Killer Cells in Osteoarthritis

Osteoarthritis (OA) is now widely acknowledged as a low-grade inflammatory condition, in which the intrinsic immune system plays a significant role in its pathogenesis. While the involvement of macrophages and T cells in the development of OA has been extensively reviewed, recent research has provided mounting evidence supporting the crucial contribution of NK cells in both the initiation and advancement of OA. Accumulated evidence has emerged in recent years indicating that NK cells play a critical role in OA development and progression. This review will outline the ongoing understanding of the utility of NK cells in the etiology of OA, focusing on how NK cells interact with chondrocytes, synoviocytes, osteoclasts, and other immune cells to influence the course of OA disease.

Overexpression of RAD54L attenuates osteoarthritis by suppressing the HIF-1α/VEGF signaling pathway: Bioinformatics analysis and experimental validation

Osteoarthritis (OA) is a widespread chronic, progressive, degenerative joint disease that causes pain and disability. Current treatments for OA have limited effectiveness and new biomarkers need to be identified. Bioinformatics analysis was conducted to explore differentially expressed genes and DNA repair/recombination protein 54 L (RAD54L) was selected. We firstly overexpressed RAD54L in interleukin-1β (IL-1β)-induced human articular chondrocytes or in OA rats to investigate its effect on OA. Chondrocyte viability and apoptotic rate were measured by Cell Counting Kit-8 and flow cytometry, respectively. Then we evaluated OA severity in vivo by Hematoxylin-eosin staining and Osteoarthritis Research Society International standards. The expression of inflammatory mediators was tested by enzyme-linked immunosorbent assay. Finally, western blot was performed to determine the relative expression level of hypoxia-inducible factors 1α (HIF-1α) and vascular endothelial growth factor (VEGF). Overexpression of RAD54L promoted cell viability and attenuated apoptosis in IL-1β-induced human chondrocytes. A lower Osteoarthritis Research Society International score and a remarkable alleviation of chondrocyte disordering and infiltration of inflammatory cells were found in cartilage tissues of OA rats after overexpressing RAD54L. The inflammatory response induced by OA was decreased by RAD54L overexpression in vitro and in vivo. In addition, RAD54L overexpression decreased the relative expression level of HIF-1α and VEGF. Overexpression of RAD54L could attenuate OA by suppressing the HIF-1α/VEGF signaling pathway, indicating that RAD54L may be a potential treatment target for OA.

Transcriptomic analyses and machine-learning methods reveal dysregulated key genes and potential pathogenesis in human osteoarthritic cartilage

Aims

This study aimed to explore the biological and clinical importance of dysregulated key genes in osteoarthritis (OA) patients at the cartilage level to find potential biomarkers and targets for diagnosing and treating OA.

Methods

Six sets of gene expression profiles were obtained from the Gene Expression Omnibus database. Differential expression analysis, weighted gene coexpression network analysis (WGCNA), and multiple machine-learning algorithms were used to screen crucial genes in osteoarthritic cartilage, and genome enrichment and functional annotation analyses were used to decipher the related categories of gene function. Single-sample gene set enrichment analysis was performed to analyze immune cell infiltration. Correlation analysis was used to explore the relationship among the hub genes and immune cells, as well as markers related to articular cartilage degradation and bone mineralization.

Results

A total of 46 genes were obtained from the intersection of significantly upregulated genes in osteoarthritic cartilage and the key module genes screened by WGCNA. Functional annotation analysis revealed that these genes were closely related to pathological responses associated with OA, such as inflammation and immunity. Four key dysregulated genes (cartilage acidic protein 1 (CRTAC1), iodothyronine deiodinase 2 (DIO2), angiopoietin-related protein 2 (ANGPTL2), and MAGE family member D1 (MAGED1)) were identified after using machine-learning algorithms. These genes had high diagnostic value in both the training cohort and external validation cohort (receiver operating characteristic > 0.8). The upregulated expression of these hub genes in osteoarthritic cartilage signified higher levels of immune infiltration as well as the expression of metalloproteinases and mineralization markers, suggesting harmful biological alterations and indicating that these hub genes play an important role in the pathogenesis of OA. A competing endogenous RNA network was constructed to reveal the underlying post-transcriptional regulatory mechanisms.

Conclusion

The current study explores and validates a dysregulated key gene set in osteoarthritic cartilage that is capable of accurately diagnosing OA and characterizing the biological alterations in osteoarthritic cartilage; this may become a promising indicator in clinical decision-making. This study indicates that dysregulated key genes play an important role in the development and progression of OA, and may be potential therapeutic targets.

Bioinformatics analysis of hub genes as osteoarthritis prognostic biomarkers

Osteoarthritis (OA) is a progressive cartilage degradation disease, concomitant with synovitis, osteophyte formation, and subchondral bone sclerosis. Over 37% of the elderly population is affected by OA, and the number of cases is increasing as the global population ages. Therefore, the objective of this study was to identify and analyze the hub genes of OA combining with comprehensive bioinformatics analysis tools to provide theoretical basis in further OA effective therapies. Two sample sets of GSE46750 contained 12 pairs OA synovial membrane and normal samples harvested from patients as well as GSE98918 including 12 OA and non-OA patients were downloaded from the Gene Expression Omnibus database (GEO) database. Differentially expressed genes (DEGs) were identified using Gene Expression Omnibus 2R (GEO2R), followed by functional enrichment analysis, protein-protein interaction networks construction. The hub genes were identified and evaluated. An OA rat model was constructed, hematoxylin and eosin staining, safranin O/fast green staining, cytokines concentrations of serum were used to verify the model. The hub genes expression level in the knee OA samples were verified using RT-qPCR. The top 20 significantly up-regulated and down-regulated DEGs were screened out from the two datasets, respectively. The top 18 GO terms and 10 KEGG pathways were enriched. Eight hub genes were identified, namely MS4A6A, C1QB, C1QC, CD74, CSF1R, HLA-DPA1, HLA-DRA and ITGB2. Among them, the hub genes were all up-regulated in in vivo OA rat model, compared with healthy controls. The eight hub genes identified (MS4A6A, C1QB, C1QC, CD74, CSF1R, HLA-DPA1, HLA-DRA and ITGB2) were shown to be associated with OA. These genes can serve as disease markers to discriminate OA patients from healthy controls.

Differences between soluble and insoluble undenatured type II collagen in improving osteoarthritis in rats and their potential mechanisms

Our previous research showed that soluble (SC II) and insoluble (IC II) undenatured type II collagen had significant differences during gastrointestinal digestion in vitro, and SC II exposed more type II collagen with triple helix structure. However, the differences in their in vivo digestive characteristics, improvement on osteoarthritis (OA), and possible mechanisms have not been elucidated. The aim of this study was to explore these issues. After oral administration of SC II and IC II, the joint swelling of OA rats significantly reduced, and the weight bearing ratio of right hind limb significantly increased, especially in SC II group (raised to 48%). The Mankin and OARSI scores decreased by 35% and 48% in SC II group, respectively. SC II and IC II increased the mRNA expression of anti-inflammatory factors and the proportion of regulatory T cells (Treg). Importantly, type II collagen released by IC II during in vivo gastrointestinal digestion was far less than SC II, which explained the higher ability of SC II to induce immune tolerance in small intestine than IC II. Bioinformatics analysis showed that the differential genes between model and control were significantly enriched in PI3K/AKT, PPAR and AMPK signalling pathways, and 24 hub genes were analyzed. SC II significantly down-regulated the mRNA expression of Il6, Ccl7, NF-κB, AKT and up-regulated the mRNA expression of Scd1. These results showed that SC II was superior to IC II in improving OA by inducing immune tolerance and could regulate key biomarkers and signalling pathways in OA rats.

Undenatured type II collagen and its role in improving osteoarthritis

Osteoarthritis (OA) is a degenerative joint disease, affecting 32.5 million US adults or 242 million people worldwide. There is no cure for OA. Many animal and clinical trials showed that oral administration of undenatured type II collagen could significantly reduce the incidence of OA or alleviate the symptoms of articular cartilage. Type II collagen is an important component of cartilage matrix. This article reviewed research progress of undenatured type II collagen including its methods of extraction and preparation, structure and characterization, solubility, thermal stability, gastrointestinal digestive stability, its role in improving OA, and the mechanism of its action in improving OA. Type II collagen has been extensively explored for its potential in improving arthritis. Methods of extraction of type II collagen are inefficient and tedious. The method of limited enzymatic hydrolysis is mainly used to prepare soluble undenatured type II collagen (SC II). The solubility, thermal and gastrointestinal digestive stability of SC II are affected by the sources of raw material, pH, salt ions, and temperature. Oral administration of undenatured type II collagen improves OA, whereas its activity is affected by the sources, degree of denaturalization, intervention methods and doses. However, the influence of the structure of undenatured type II collagen on its activity and the mechanism are unclear. The findings in this review support that undenatured type II collagen can be used in the intervention or auxiliary intervention of patients with OA.

Identification of Ferroptosis-related potential biomarkers and immunocyte characteristics in Chronic Thromboembolic Pulmonary Hypertension via bioinformatics analysis

BACKGROUND

Chronic Thromboembolic Pulmonary Hypertension (CTEPH) is a form of pulmonary hypertension with a high mortality rate. A new type of iron-mediated cell death is Ferroptosis, which is characterized by the accumulation of lethal iron ions and lipid peroxidation leading to mitochondrial atrophy and increased mitochondrial membrane density. Now, there is a lack of Ferroptosis-related biomarkers (FRBs) associated with pathogenic process of CTEPH.

METHODS

The differentially expressed genes (DEGs) of CTEPH were obtained by GEO2R. Genes related to Ferroptosis were obtained from FerrDb database. The intersection of Ferroptosis and DEGs results in FRBs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were performed in Database for Annotation, Visualization and Integrated Discovery (DAVID) database. The optimal potential biomarkers for CTEPH were analyzed by least absolute shrinkage and selection operator (LASSO) and support vector machine-recursive feature elimination (SVM-RFE) machine learning. The four hub genes were verified from the Gene Expression Omnibus (GEO) dataset GSE188938. Immune infiltration was analyzed by CIBERSORT. SPSS software was used to analyze the Spearman rank correlation between FRBs identified and infiltration-related immune cells, and p < 0.05 was considered as statistically significant.

RESULTS

In this study, potential genetic biomarkers associated with Ferroptosis in CTEPH were investigated and explored their role in immune infiltration. In total, we identified 17 differentially expressed Ferroptosis-associated genes by GEOquery package. The key FRBs including ARRDC3, HMOX1, BRD4, and YWHAE were screened using Lasso and SVM-RFE machine learning methods.Through gene set GSE188938 verification, only upregulation of gene ARRDC3 showed statistical difference. In addition, immune infiltration analysis using the CIBERSORT algorithm revealed the infiltration of Eosinophils and Neutrophils in CTEPH samples was less than that in the control group. And correlation analysis revealed that ARRDC3 was positively correlated with T cells follicular helper (r = 0.554, p = 0.017) and negatively correlated with Neutrophils (r = -0.47, p = 0.049).

CONCLUSIONS

In conclusion, ARRDC3 upregulation with different immune cell infiltration were involved in the development of CTEPH. ARRDC3 might a potential Ferroptosis-related biomarker for CTEPH treatment. This study provided a new insight into pathogenesis CTEPH.

Appraising causal risk and protective factors for rheumatoid arthritis

Aims

Mendelian randomization (MR) is considered to overcome the bias of observational studies, but there is no current meta-analysis of MR studies on rheumatoid arthritis (RA). The purpose of this study was to summarize the relationship between potential pathogenic factors and RA risk based on existing MR studies.

Methods

PubMed, Web of Science, and Embase were searched for MR studies on influencing factors in relation to RA up to October 2022. Meta-analyses of MR studies assessing correlations between various potential pathogenic factors and RA were conducted. Random-effect and fixed-effect models were used to synthesize the odds ratios of various pathogenic factors and RA. The quality of the study was assessed using the Strengthening the Reporting of Observational Studies in Epidemiology using Mendelian Randomization (STROBE-MR) guidelines.

Results

A total of 517 potentially relevant articles were screened, 35 studies were included in the systematic review, and 19 studies were eligible to be included in the meta-analysis. Pooled estimates of 19 included studies (causality between 15 different risk factors and RA) revealed that obesity, smoking, coffee intake, lower education attainment, and Graves' disease (GD) were related to the increased risk of RA. In contrast, the causality contribution from serum mineral levels (calcium, iron, copper, zinc, magnesium, selenium), alcohol intake, and chronic periodontitis to RA is not significant.

Conclusion

Obesity, smoking, education attainment, and GD have real causal effects on the occurrence and development of RA. These results may provide insights into the genetic susceptibility and potential biological pathways of RA.

Dihydrocaffeic acid improves IL-1β-induced inflammation and cartilage degradation via inhibiting NF-κB and MAPK signalling pathways

Aims

Osteoarthritis (OA) is a prevalent joint disorder with inflammatory response and cartilage deterioration as its main features. Dihydrocaffeic acid (DHCA), a bioactive component extracted from natural plant (gynura bicolor), has demonstrated anti-inflammatory properties in various diseases. We aimed to explore the chondroprotective effect of DHCA on OA and its potential mechanism.

Methods

In vitro, interleukin-1 beta (IL-1β) was used to establish the mice OA chondrocytes. Cell counting kit-8 evaluated chondrocyte viability. Western blotting analyzed the expression levels of collagen II, aggrecan, SOX9, inducible nitric oxide synthase (iNOS), IL-6, matrix metalloproteinases (MMPs: MMP1, MMP3, and MMP13), and signalling molecules associated with nuclear factor-kappa B (NF-κB) and mitogen-activated protein kinase (MAPK) pathways. Immunofluorescence analysis assessed the expression of aggrecan, collagen II, MMP13, and p-P65. In vivo, a destabilized medial meniscus (DMM) surgery was used to induce mice OA knee joints. After injection of DHCA or a vehicle into the injured joints, histological staining gauged the severity of cartilage damage.

Results

DHCA prevented iNOS and IL-6 from being upregulated by IL-1β. Moreover, the IL-1β-induced upregulation of MMPs could be inhibited by DHCA. Additionally, the administration of DHCA counteracted IL-1β-induced downregulation of aggrecan, collagen II, and SOX9. DHCA protected articular cartilage by blocking the NF-κB and MAPK pathways. Furthermore, DHCA mitigated the destruction of articular cartilage in vivo.

Conclusion

We present evidence that DHCA alleviates inflammation and cartilage degradation in OA chondrocytes via suppressing the NF-κB and MAPK pathways, indicating that DHCA may be a potential agent for OA treatment.

Cite this article: Bone Joint Res 2023;12(4):259–273.

A two-gene random forest model to diagnose osteoarthritis based on RNA-binding protein-related genes in knee cartilage tissue

Osteoarthritis (OA) is one of the most common diseases in the orthopedic clinic, characterized by progressive cartilage degradation. RNA-binding proteins (RBPs) are capable of binding to RNAs at transcription and translation levels, playing an important role in the pathogenesis of OA. This study aims to investigate the diagnosis values of RBP-related genes in OA. The RBPs were collected from previous studies, and the GSE114007 dataset (control = 18, OA = 20) was downloaded from the Gene Expression Omnibus (GEO) as the training cohort. Through various bioinformatical and machine learning methods, including genomic difference detection, protein-protein interaction network analyses, Lasso regression, univariate logistic regression, Boruta algorithm, and SVM-RFE, RNMT and RBM24 were identified and then included into the random forest (RF) diagnosis model. GSE117999 dataset (control = 10, OA = 10) and clinical samples collected from local hospital (control = 10, OA = 11) were used for external validation. The RF model was a promising tool to diagnose OA in the training dataset (area under curve [AUC] = 1.000, 95% confidence interval [CI] = 1.000-1.000), the GSE117999 cohort (AUC = 0.900, 95% CI = 0.769-1.000), and local samples (AUC = 0.759, 95% CI = 0.568-0.951). Besides, qPCR and Western Blotting experiments showed that RNMT (P < 0.05) and RBM24 (P < 0.01) were both down-regulated in CHON-001 cells with IL-1β treatment. In all, an RF model to diagnose OA based on RNMT and RBM24 in cartilage tissue was constructed, providing a promising clinical tool and possible cut-in points in molecular mechanism clarification.

Characterization and role exploration of ferroptosis-related genes in osteoarthritis

Osteoarthritis (OA), viewing as a degenerative aseptic inflammatory disease, is characterized by joint pain and inflammation that significantly affects the quality of patients' life, especially for the elder. Although rapid progress has been achieved in elucidating the underlying mechanisms of OA occurrence and progression, there is still a lack of effective clinical therapeutics for OA patients. Currently the most common treatments including drug therapy and surgical operations are not very satisfactory in majority of cases, so it is worthy to explore new remedies. During the past few decades, a number of novel forms of regulated cell death have been reported widely, typified by ferroptosis, with its prominent features including reactive oxygen species (ROS) elevation, lipid peroxidation, iron accumulation and glutathione deprivation. Our study was designed to identify the functional roles of differentially expressed ferroptosis-related genes in OA, which were screened out by referring to GEO database via bioinformatics analyses. Human chondrocytes were applied to validate the above findings in the scenario of ferroptosis inhibitors administration. Results partially proved the consistency with bioinformatics analyses that ATF3 and TFRC were highly expressed in interleukin-1β (IL-1β)-stimulated chondrocytes whereas CXCL2 and JUN were downregulated. Besides, TFRC was firstly validated to be upregulated in IL-1β-stimulated chondrocytes, which could be reversed by ferroptosis inhibitors. In conclusion, our study reported two prominent ferroptosis-related genes, ATF3 and TFRC are upregulated in IL-1β-stimulated chondrocytes while CXCL2 and JUN are downregulated. And preliminary results demonstrated that TFRC might serve as an accomplice of ferroptosis process in IL-1β-stimulated chondrocytes and ferroptosis inhibitors have the potential to inhibit ROS in IL-1β-stimulated chondrocytes.