Whole-transcriptome sequencing of knee joint cartilage from osteoarthritis patients

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.

Whole-Transcriptome Sequencing of Knee Joint Cartilage from Kashin-Beck Disease and Osteoarthritis Patients

The aim of this study was to provide a comprehensive understanding of similarities and differences in mRNAs, lncRNAs, and circRNAs within cartilage for Kashin-Beck disease (KBD) compared to osteoarthritis (OA). We conducted a comparison of the expression profiles of mRNAs, lncRNAs, and circRNAs via whole-transcriptome sequencing in eight KBD and ten OA individuals. To facilitate functional annotation-enriched analysis for differentially expressed (DE) genes, DE lncRNAs, and DE circRNAs, we employed bioinformatic analysis utilizing Gene Ontology (GO) and KEGG. Additionally, using quantitative reverse transcriptase polymerase chain reaction (qRT-PCR), we validated the expression levels of four cartilage-related genes in chondrocytes. We identified a total of 43 DE mRNAs, 1451 DE lncRNAs, and 305 DE circRNAs in KBD cartilage tissue compared to OA (q value < 0.05; |log2FC| > 1). We also performed competing endogenous RNA network analysis, which identified a total of 65 lncRNA-mRNA interactions and 4714 miRNA-circRNA interactions. In particular, we observed that circRNA12218 had binding sites for three miRNAs targeting ACAN, while circRNA12487 had binding sites for seven miRNAs targeting COL2A1. Our results add a novel set of genes and non-coding RNAs that could potentially serve as candidate diagnostic biomarkers or therapeutic targets for KBD patients.

Single-cell RNA sequencing in osteoarthritis

Osteoarthritis is a progressive and heterogeneous joint disease with complex pathogenesis. The various phenotypes associated with each patient suggest that better subgrouping of tissues associated with genotypes in different phases of osteoarthritis may provide new insights into the onset and progression of the disease. Recently, single-cell RNA sequencing was used to describe osteoarthritis pathogenesis on a high-resolution view surpassing traditional technologies. Herein, this review summarizes the microstructural changes in articular cartilage, meniscus, synovium and subchondral bone that are mainly due to crosstalk amongst chondrocytes, osteoblasts, fibroblasts and endothelial cells during osteoarthritis progression. Next, we focus on the promising targets discovered by single-cell RNA sequencing and its potential applications in target drugs and tissue engineering. Additionally, the limited amount of research on the evaluation of bone-related biomaterials is reviewed. Based on the pre-clinical findings, we elaborate on the potential clinical values of single-cell RNA sequencing for the therapeutic strategies of osteoarthritis. Finally, a perspective on the future development of patient-centred medicine for osteoarthritis therapy combining other single-cell multi-omics technologies is discussed. This review will provide new insights into osteoarthritis pathogenesis on a cellular level and the field of applications of single-cell RNA sequencing in personalized therapeutics for osteoarthritis in the future.

MicroRNA expression in osteoarthritis: a meta-analysis

Osteoarthritis (OA) is one of the most prevalent musculoskeletal diseases globally, leading to chronic disability and poor prognosis. One of the approaches for optimizing OA treatment is to find early effective diagnostic biomarkers. The contribution of microRNAs (miRNAs) in OA progression is now being increasingly recognized. This review provides a comprehensive summary on studies reporting the expression profiling of miRNAs in OA and associated signaling pathways. We performed a systematic search of the Embase, Web of Science, PubMed, and Cochrane library databases. This systematic review is reported according to the PRISMA checklist. Studies which identified miRNAs with aberrant expression compared to controls during OA progression were included, and a meta-analysis was performed. Results from the random effects model were provided as log10 odds ratios (logORs) and 95% confidence intervals. Sensitivity analysis was conducted to confirm the accuracy of the results. Subgroup analysis was conducted based on tissue source. The target genes of miRNAs identified in this study were extracted from the MiRWalk database, and these target genes were enriched in Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathways. A total of 191 studies reporting 162 miRNAs were included in our meta-analysis. Among them, 36 miRNAs distributed across 96 studies were expressed in the same direction in at least two studies (13 up-regulated and 23 down-regulated). Subgroup analysis of tissue source revealed that the highest number of studies was conducted using articular cartilage, where the most up-regulated miRNAs were miR-146a-5p (logOR 7.355; P < 0.001) and miR-34a-5p (logOR 6.955; P < 0.001), and the most down-regulated miRNAs were miR-127-5p (logOR 6.586; P < 0.001) and miR-140-5p (logOR 6.373; P < 0.001). Enrichment analysis of 752 downstream target genes of all identified miRNAs was performed, and the regulatory relationships among them were displayed. Mesenchymal stem cells and transforming growth factor-β were found to be the most important downstream effectors regulated by miRNA in OA. This study highlighted the importance of miRNA signaling in OA progression and identified a number of prominent miRNAs including miR-146a-5p, miR-34a-5p, miR-127-5p, and miR-140-5p which might be considered as potential biomarkers for OA.

Paeoniflorin shows chondroprotective effects under IL-1β stress by regulating circ-PREX1/miR-140-3p/WNT5B axis

BACKGROUND

Osteoarthritis (OA) is a chronic and degenerative bone and joint disease, and paeoniflorin shows anti-arthritis role in OA. This study planned to investigate the mechanism related to chondroprotective role of paeoniflorin in OA.

METHODS

Real-time quantitative PCR and western blotting were performed to measure expression levels of circ-PREX1, microRNA (miR)-140-3p, Wingless-type MMTV integration site family, member 5B (WNT5B), B cell lymphoma (Bcl)-2, and Bcl-2 Associated X Protein (Bax). MTT assay, EdU assay, flow cytometry and enzyme-linked immunosorbent assay evaluated cell viability, proliferation, apoptosis and inflammatory response, respectively. Dual-luciferase reporter assay and RNA immunoprecipitation assay identified the relationship among circ-PREX1, miR-140-3p, and WNT5B.

RESULTS

IL-1β highly induced apoptosis rate, Bax expression and TNF-α product, accompanied with decreased cell viability, cell proliferation and IL-10 secretion, whereas these effects were partially reversed after paeoniflorin pretreatment. Expression of circ-PREX1 was upregulated and miR-140-3p was downregulated in cartilage tissues of patients with knee OA (KOA) and IL-1β-induced human chondrocytes (C28/I2). Circ-PREX1 overexpression and miR-140-3p silencing attenuated the suppressive effect of paeoniflorin in IL-1β-induced C28/I2 cells. Furthermore, miR-140-3p was negatively regulated by circ-PREX1. WNT5B was a downstream target of miR-140-3p and could be modulated by the circ-PREX1/miR-140-3p pathway in IL-1β-induced C28/I2 cells.

CONCLUSION

Paeoniflorin might protect human chondrocytes from IL-1β-induced inflammatory injury via circ-PREX1-miR-140-3p-WNT5B pathway, suggesting a potential preventative agent and a novel target for the treatment of KOA.

CircSFMBT2-OA alleviates chondrocyte apoptosis and extracellular matrix degradation through repressing NF-κB/NLRP3 inflammasome activation

Background

Intra-articular inflammation and cartilage degradation are the major pathological characteristics of osteoarthritis (OA). Mounting studies have revealed that circular RNAs (circRNAs) act as an important regulatory role in inflammatory diseases and are frequently dys-expressed in OA cartilage tissues.

Objective

Here, a dys-regulated cicrRNA (has_circ_0017636, termed circSFMBT2-OA) was identified, and its role in regulating lipopolysaccharide (LPS)-induced chondrocyte injury was next investigated.

Methods

CHON-001 chondrocytes were treated with LPS, and then the levels of circSFMBT2-OA, cartilage-related genes, and pro-inflammatory cytokines were measured using quantitative real-time PCR (qRT-PCR) and Western blot analysis. CHON-001 cell viability, proliferation, and apoptosis were assayed using Cell Counting Kit-8 (CCK-8), 5-Ethynyl-2'-deoxyuridine (EDU), and terminal deoxynucleotidyl transferase (TdT) dUTP nick-end labeling (TUNEL) assay, respectively.

Results

CircSFMBT2-OA level was significantly down-regulated in OA cartilage tissues and LPS-treated CHON-001 cells. Functionally, circSFMBT2-OA overexpression accelerated cell proliferation, and suppressed cell apoptosis, pro-inflammatory cytokines production, matrix-degrading enzymes expression, and ECM degradation in CHON-001 cells. Inversely, circSFMBT2-OA depletion decreased cell viability and increased matrix-degrading enzymes expression and ECM degradation. Mechanistically, circSFMBT2-OA inhibited LPS-induced NF-κB/NOD-like receptor family pyrin domain containing protein 3 (NLRP3) inflammasome activation in CHON-001 cells. Consequently, NLRP3 activator reversed the effect of circSFMBT2-OA on repressing LPS-induced CHON-001 cell injury.

Conclusion

These data reveal a vital effect of a novel circSFMBT2-OA on repressing OA progression and provide a promising target to treat OA.

CircRNA-mediated ceRNA mechanism in Osteoarthritis: special emphasis on circRNAs in exosomes and the crosstalk of circRNAs and RNA methylation

Osteoarthritis (OA) is an age-related joint disease with chronic inflammation, progressive articular cartilage destruction and subchondral bone sclerosis. CircRNAs (circRNAs) are a class of non-coding RNA with a circular structure that participate in a series of important pathophysiological processes of OA, especially its ceRNA mechanisms, and play an important role in OA. CircRNAs may be potential biomarkers for the diagnosis and prognosis of OA. Additionally, differentially expressed circRNAs were found in patients with OA, indicating that circRNAs are involved in the pathogenesis of OA. Experiments have shown that the intra-articular injection of modified circRNAs can effectively relieve OA. Exosomal circRNAs and methylated circRNAs also provide new ideas for the treatment of OA. Clarifying the important roles of circRNAs in OA will deepen people's understanding of the pathogenesis of OA. CircRNAs may be developed as new biomarkers or drug targets for the diagnosis of OA and provide new methods for the treatment of OA.

Implications and theragnostic potentials of circular RNAs in rheumatic diseases

Circular RNAs (circRNAs), a class of non-coding RNAs (ncRNAs), are highly stable and ubiquitous molecules that exhibit tissue-specific expression. Accumulating evidence has shown that aberrant expression of circRNAs can play a role in the pathogenesis of several diseases. Rheumatic diseases are a varied group of autoimmune and inflammatory disorders affecting mainly the musculoskeletal system. Notably, circRNAs, which are essential immune system gene modulators, are strongly linked to the occurrence and progression of autoimmune disorders. Here, we present and discuss the current findings concerning the roles, implications and theragnostic potentials of circRNAs in common rheumatic diseases, including ankylosing spondylitis (AS), osteoarthritis (OA), osteoporosis (OP), rheumatoid arthritis (RA), systemic lupus erythematosus (SLE), Crohn's disease (CD), and gout. This review aims to provide new insights to support the development of novel diagnostic and therapeutic strategies for these disabling diseases.

Whole-transcriptome sequencing and ceRNA interaction network of temporomandibular joint osteoarthritis

Purpose: The aim of this study was to conduct a comprehensive transcriptomic analysis to explore the potential biological functions of noncoding RNA (ncRNAs) in temporomandibular joint osteoarthritis (TMJOA). Methods: Whole transcriptome sequencing was performed to identify differentially expressed genes (DEGs) profiles between the TMJOA and normal groups. The functions and pathways of the DEGs were analyzed using Metascape, and a competitive endogenous RNA (ceRNA) network was constructed using Cytoscape software. Results: A total of 137 DEmRNAs, 65 DEmiRNAs, 132 DElncRNAs, and 29 DEcircRNAs were identified between the TMJOA and normal groups. Functional annotation of the DEmRNAs revealed that immune response and apoptosis are closely related to TMJOA and also suggested key signaling pathways related to TMJOA, including chronic depression and PPAR signaling pathways. We identified vital mRNAs, including Klrk1, Adipoq, Cryab, and Hspa1b. Notably, Adipoq expression in cartilage was significantly upregulated in TMJOA compared with normal groups (10-fold, p < 0.001). According to the functional analysis of DEmRNAs regulated by the ceRNA network, we found that ncRNAs are involved in the regulation of autophagy and apoptosis. In addition, significantly DEncRNAs (lncRNA-COX7A1, lncRNA-CHTOP, lncRNA-UFM1, ciRNA166 and circRNA1531) were verified, and among these, circRNA1531 (14.5-fold, p < 0.001) and lncRNA-CHTOP (14.8-fold, p < 0.001) were the most significantly downregulated ncRNAs. Conclusion: This study showed the potential of lncRNAs, circRNAs, miRNAs, and mRNAs may as clinical biomarkers and provides transcriptomic insights into their functional roles in TMJOA. This study identified the transcriptomic signatures of mRNAs associated with immunity and apoptosis and the signatures of ncRNAs associated with autophagy and apoptosis and provides insight into ncRNAs in TMJOA.

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

AIMS

This study aimed, through bioinformatics analysis, to identify the potential diagnostic markers of osteoarthritis, and analyze the role of immune infiltration in synovial tissue.

METHODS

The gene expression profiles were downloaded from the Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) were identified by R software. Functional enrichment analyses were performed and protein-protein interaction networks (PPI) were constructed. Then the hub genes were screened. Biomarkers with high value for the diagnosis of early osteoarthritis (OA) were validated by GEO datasets. Finally, the CIBERSORT algorithm was used to evaluate the immune infiltration between early-stage OA and end-stage OA, and the correlation between the diagnostic marker and infiltrating immune cells was analyzed.

RESULTS

A total of 88 DEGs were identified. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses indicated that DEGs were significantly enriched in leucocyte migration and interleukin (IL)-17 signalling pathways. Disease ontology (DO) indicated that DEGs were mostly enriched in rheumatoid arthritis. Six hub genes including FosB proto-oncogene, AP-1 transcription factor subunit (FOSB); C-X-C motif chemokine ligand 2 (CXCL2); CXCL8; IL-6; Jun proto-oncogene, AP-1 transcription factor subunit (JUN); and Activating transcription factor 3 (ATF3) were identified and verified by GEO datasets. ATF3 (area under the curve = 0.975) turned out to be a potential biomarker for the diagnosis of early OA. Several infiltrating immune cells varied significantly between early-stage OA and end-stage OA, such as resting NK cells (p = 0.016), resting dendritic cells (p = 0.043), and plasma cells (p = 0.043). Additionally, ATF3 was significantly correlated with resting NK cells (p = 0.034), resting dendritic cells (p = 0.026), and regulatory T cells (Tregs, p = 0.018).

CONCLUSION

ATF3 may be a potential diagnostic marker for early diagnosis and treatment of OA, and immune cell infiltration provides new perspectives for understanding the mechanism during OA progression.Cite this article: Bone Joint Res 2022;11(9):679-689.

ASPN Is a Potential Biomarker and Associated with Immune Infiltration in Endometriosis

Objective: Endometriosis is a benign gynecological disease characterized by distant metastasis. Previous studies have discovered abnormal numbers and function of immune cells in endometriotic lesions. We aimed to find potential biomarkers of endometriosis and to explore the relationship between ASPN and the immune microenvironment of endometriosis. Methods: We obtained the GSE141549 and GSE7305 datasets containing endometriosis and normal endometrial samples from the Gene Expression Omnibus database (GEO). In the GSE141549 dataset, differentially expressed genes (DEGs) were found. The Least Absolute Shrinkage and Selection Operator (Lasso) regression and generalized linear models (GLMs) were used to screen new biomarkers. The expression levels and diagnostic utility of biomarkers were assessed in GSE7305, and biomarker expression levels were further validated using qRT-PCR and western blot. We identified DEGs between high and low expression groups of key biomarkers. Enrichment analysis was carried out to discover the target gene’s biological function. We analyzed the relationship between key biomarker expression and patient clinical features. Finally, the immune cells that infiltrate endometriosis were assessed using the Microenvironment Cell Population-Counter (MCP-counter), and the correlation of biomarker expression with immune cell infiltration and immune checkpoints genes was studied. Results: There were a total of 38 DEGs discovered. Two machine learning techniques were used to identify 10 genes. Six biomarkers (SCG2, ASPN, SLIT2, GEM, EGR1, and FOS) had good diagnostic efficiency (AUC > 0.7) by internal and external validation. We excluded previously reported related genes (SLIT2, EGR1, and FOS). ASPN was the most significantly differentially expressed biomarker between normal and ectopic endometrial tissues, as verified by qPCR. The western blot assay revealed a significant upregulation of ASPN expression in endometriotic tissues. The investigation for DEGs in the ASPN high- and low-expression groups revealed that the DEGs were particularly enriched in extracellular matrix tissue, vascular smooth muscle contraction, cytokine interactions, the calcium signaling pathway, and the chemokine signaling pathway. High ASPN expression was related to r-AFS stage (p = 0.006), age (p = 0.03), and lesion location (p < 0.001). Univariate and multivariate logistic regression analysis showed that ASPN expression was an independent influencing factor in patients with endometriosis. Immune cell infiltration analysis revealed a significant increase in T-cell, B-cell, and fibroblast infiltration in endometriosis lesions; cytotoxic lymphocyte, NK-cell, and endothelial cell infiltration were reduced. Additionally, the percentage of T cells, B cells, fibroblasts, and endothelial cells was favorably connected with ASPN expression, while the percentage of cytotoxic lymphocytes and NK cells was negatively correlated. Immune checkpoint gene (CTLA4, LAG3, CD27, CD40, and ICOS) expression and ASPN expression were positively associated. Conclusions: Increased expression of ASPN is associated with immune infiltration in endometriosis, and ASPN can be used as a diagnostic biomarker as well as a potential immunotherapeutic target in endometriosis.

The Molecular Mechanism of Long Non-Coding RNA MALAT1-Mediated Regulation of Chondrocyte Pyroptosis in Ankylosing Spondylitis

Long non-coding RNAs (lncRNAs) may be important regulators in the progression of ankylosing spondylitis (AS). The competing endogenous RNA (ceRNA) activity of lncRNAs plays crucial roles in osteogenesis. We identified the mechanism of the differentially expressed lncRNA MALAT1 in AS using bioinformatic analysis and its ceRNA mechanism. The interaction of MALAT1, microRNA-558, and GSDMD was identified using integrated bioinformatics analysis and validated. Loss- and gain-of-function assays evaluated their effects on the viability, apoptosis, pyroptosis and inflammation of chondrocytes in AS. We found elevated MALAT1 and GSDMD but reduced miR-558 in AS cartilage tissues and chondrocytes. MALAT1 contributed to the suppression of cell viability and facilitated apoptosis and pyroptosis in AS chondrocytes. GSDMD was a potential target gene of miR-558. Depletion of MALAT1 expression elevated miR-558 by inhibiting GSDMD to enhance cell viability and inhibit inflammation, apoptosis and pyroptosis of chondrocytes in AS. In summary, our key findings demonstrated that knockdown of MALAT1 served as a potential suppressor of AS by upregulating miR-558 via the downregulation of GSDMD expression.

Cartilage Homeostasis and Osteoarthritis

Healthy limb joints are important for maintaining health and attaining longevity. Endochondral ossification (the replacement of cartilage with bone, occurring during skeletal development) is essential for bone formation, especially in long-axis bones. In contrast to endochondral ossification, chondrocyte populations in articular cartilage persist and maintain joint tissue into adulthood. Articular cartilage, a connective tissue consisting of chondrocytes and their surrounding extracellular matrices, plays an essential role in the mechanical cushioning of joints in postnatal locomotion. Osteoarthritis (OA) pathology relates to disruptions in the balance between anabolic and catabolic signals, that is, the loss of chondrocyte homeostasis due to aging or overuse of cartilages. The onset of OA increases with age, shortening a person’s healthy life expectancy. Although many people with OA experience pain, the mainstay of treatment is symptomatic therapy, and no fundamental treatment has yet been established. To establish regenerative or preventative therapies for cartilage diseases, further understanding of the mechanisms of cartilage development, morphosis, and homeostasis is required. In this review, we describe the general development of cartilage and OA pathology, followed by a discussion on anabolic and catabolic signals in cartilage homeostasis, mainly microRNAs.

RNA-seq and Network Analysis Reveal Unique Chemokine Activity Signatures in the Synovial Tissue of Patients With Rheumatoid Arthritis

Purpose

This study aimed to provide a comprehensive understanding of the genome-wide expression patterns in the synovial tissue samples of patients with rheumatoid arthritis (RA) to investigate the potential mechanisms regulating RA occurrence and development.

Methods

Transcription profiles of the synovial tissue samples from nine patients with RA and 15 patients with osteoarthritis (OA) (control) from the East Asian population were generated using RNA sequencing (RNA-seq). Gene set enrichment analysis (GSEA) was used to analyze all the detected genes and the differentially expressed genes (DEGs) were identified using DESeq. To further analyze the DEGs, the Gene Ontology (GO) functional enrichment and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed. The protein-protein interaction (PPI) network of the DEGs was constructed using the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) and the hub genes were identified by topology clustering with the Molecular Complex Detection (MCODE)-Cytoscape. The most important hub genes were validated using quantitative real-time PCR (qRT-PCR).

Results

Of the 17,736 genes detected, 851 genes were identified as the DEGs (474 upregulated and 377 downregulated genes) using the false discovery rate (FDR) approach. GSEA revealed that the significantly enriched gene sets that positively correlated with RA were CD40 signaling overactivation, Th1 cytotoxic module, overactivation of the immune response, adaptive immune response, effective vs. memory CD8+ T cells (upregulated), and naïve vs. effective CD8+ T cells (downregulated). Biological process enrichment analysis showed that the DEGs were significantly enriched for signal transduction (P = 3.01 × 10-6), immune response (P = 1.65 × 10-24), and inflammatory response (P = 5.76 × 10-10). Molecule function enrichment analysis revealed that the DEGs were enriched in calcium ion binding (P = 1.26 × 10-5), receptor binding (P = 1.26 × 10-5), and cytokine activity (P = 2.01 × 10-3). Cellular component enrichment analysis revealed that the DEGs were significantly enriched in the plasma membrane (P = 1.91 × 10-31), an integral component of the membrane (P = 7.39 × 10-13), and extracellular region (P = 7.63 × 10-11). The KEGG pathway analysis showed that the DEGs were enriched in the cytokine-cytokine receptor interaction (P = 3.05 × 10-17), chemokine signaling (P = 3.50 × 10-7), T-cell receptor signaling (P = 5.17 × 10-4), and RA (P = 5.17 × 10-4) pathways. We confirmed that RA was correlated with the upregulation of the PPI network hub genes, such as CXCL13, CXCL6, CCR5, CXCR5, CCR2, CXCL3, and CXCL10, and the downregulation of the PPI network hub gene such as SSTR1.

Conclusion

This study identified and validated the DEGs in the synovial tissue samples of patients with RA, which highlighted the activity of a subset of chemokine genes, thereby providing novel insights into the molecular mechanisms of RA pathogenesis and identifying potential diagnostic and therapeutic targets for RA.

Silencing of circ_0000205 mitigates interleukin-1β-induced apoptosis and extracellular matrix degradation in chondrocytes via targeting miR-766-3p/ADAMTS5 axis

The aim of this study was to explore the role of hsa_circRNA_0000205 (circ_0000205) in chondrocyte injury in osteoarthritis (OA) and the underlying mechanism. Expression of circ_0000205, microRNA (miR)-766-3p and a disintegrin and metalloproteinase with thrombospondin motif (ADAMTS)-5 was detected by quantitative real time (qRT)-polymerase chain reaction (PCR) and Western blot assays. Cell proliferation, apoptosis, and extracellular matrix (ECM) synthesis were examined by 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide and 5-ethynyl-2-deoxyuridine assays, flow cytometry, and qRT-PCR and Western blot assays. The target relationship between miR-766-3p and circ_0000205 or ADAMTS5 was confirmed by luciferase reporter assay and RNA immunoprecipitation. IL-1β treatment could attenuate cell viability of primary chondrocytes and proliferating cell nuclear antigen (PCNA) and collagen II type alpha-1 (COL2A1) levels, and elevate apoptosis rate and cleaved caspase-3, ADAMTS5 and matrix metalloproteinase-13 (MMP13) levels, suggesting that IL-1β induced chondrocyte apoptosis and ECM degradation. Expression of circ_0000205 was up-regulated in OA tissues and IL-1β-induced primary chondrocytes, accompanied with miR-766-3p down-regulation and ADAMTS5 up-regulation. Knockdown of circ_0000205 could mitigate IL-1β-induced above effects and improve cell proliferation. Moreover, both depleting miR-766-3p and promoting ADAMTS5 could partially counteract circ_0000205 knockdown roles in IL-1β-cultured primary chondrocytes. Notably, circ_0000205 was verified as a sponge for miR-766-3p via targeting, and ADAMTS5 was a direct target for miR-766-3p. Silencing circ_0000205 could protect chondrocytes from IL-1β-induced proliferation reduction, apoptosis, and ECM degradation by targeting miR-766-3p/ADAMTS5 axis.

CircSCAPER contributes to IL-1β-induced osteoarthritis in vitro via miR-140-3p/EZH2 axis

Aims

Circular RNA (circRNA) S-phase cyclin A-associated protein in the endoplasmic reticulum (ER) (circSCAPER, ID: hsa_circ_0104595) has been found to be highly expressed in osteoarthritis (OA) patients and has been associated with the severity of OA. Hence, the role and mechanisms underlying circSCAPER in OA were investigated in this study.

Methods

In vitro cultured human normal chondrocyte C28/I2 was exposed to interleukin (IL)-1β to mimic the microenvironment of OA. The expression of circSCAPER, microRNA (miR)-140-3p, and enhancer of zeste homolog 2 (EZH2) was detected using quantitative real-time polymerase chain reaction and Western blot assays. The extracellular matrix (ECM) degradation, proliferation, and apoptosis of chondrocytes were determined using Western blot, cell counting kit-8, and flow cytometry assays. Targeted relationships were predicted by bioinformatic analysis and verified using dual-luciferase reporter and RNA immunoprecipitation (RIP) assays. The levels of phosphoinositide 3-kinase (PI3K)/protein kinase B (AKT) pathway-related protein were detected using Western blot assays.

Results

CircSCAPER was highly expressed in OA cartilage tissues and IL-1β-induced chondrocytes. Knockdown of circSCAPER reduced IL-1β-evoked ECM degradation, proliferation arrest, and apoptosis enhancement in chondrocytes. Mechanistically, circSCAPER directly bound to miR-140-3p, and miR-140-3p inhibition reversed the effects of circSCAPER knockdown on IL-1β-induced chondrocytes. miR-140-3p was verified to target EZH2, and overexpression of miR-140-3p protected chondrocytes against IL-1β-induced dysfunction via targeting EZH2. Additionally, we confirmed that circSCAPER could regulate EZH2 through sponging miR-140-3p, and the circSCAPER/miR-140-3p/EZH2 axis could activate the PI3K/AKT pathway.

Conclusion

CircSCAPER promoted IL-1β-evoked ECM degradation, proliferation arrest, and apoptosis enhancement in chondrocytes via regulating miR-140-3p/EZH2 axis, which gained a new insight into the pathogenesis of OA.

Cite this article: Bone Joint Res 2022;11(2):61–72.

Circ_0020093 Overexpression Alleviates Interleukin-1 Beta-induced Inflammation, Apoptosis and Extracellular Matrix Degradation in Human Chondrocytes by Targeting the miR-181a-5p/ERG Pathway

Osteoarthritis (OA) is a well-known chronic degenerative joint disease, with multiple changes in the phenotype of chondrocytes. Circular RNAs (circRNAs) have been shown to be involved in various human diseases, including OA. The purpose of this study was to determine the role of circ_0020093 in OA pathological changes in vitro. C28/I2 cells were treated with interleukin-1 beta (IL-1β) to mimic OA pathological conditions. The expression levels of circ_0020093, miR-181a-5p and ETS-related gene (ERG) mRNA were measured by quantitative real-time PCR (qRT-PCR). For functional analyses, cell proliferative capacity was detected using EdU assay and CCK-8 assay. Inflammatory response was assessed by determining the release of pro-inflammatory factors using ELISA kits. Cell apoptosis was examined by flow cytometry assay. The levels of apoptosis-related proteins and extracellular matrix (ECM)-associated proteins were assessed by Western blot. The binding relationship between miR-181a-5p and circ_0020093 or ERG was confirmed by RNA pull-down assay, dual-luciferase reporter assay or RIP assay. The expression level of circ_0020093 was decreased in IL-1β-treated C28/I2 cells. Circ_0020093 overexpression relieved inflammatory responses, cell apoptosis and ECM degradation in IL-1β-induced C28/I2 cells. Circ_0020093 directly targeted miR-181a-5p, and miR-181a-5p bound to the 3' -untranslated region (3'UTR) of ERG to regulate ERG expression. Circ_0020093 overexpression promoted the expression of ERG by sponging miR-181a-5p. Rescue experiments showed that miR-181a-5p overexpression or ERG knockdown could reverse the inhibitory effects of circ_0020093 overexpression on the pathological changes in IL-1β-induced C28/I2 cells. Circ_0020093 overexpression alleviated IL-1β-induced human chondrocyte inflammatory injury, apoptosis and ECM degradation by targeting miR-181a-5p/ERG pathway.

Compatibility of Achyranthes bidentata components in reducing inflammatory response through Arachidonic acid pathway for treatment of Osteoarthritis

Achyranthes bidentate is a common traditional Chinese medicine (TCM) used in treating osteoarthritis (OA). The compatibility between effective components has now become a breakthrough in understanding the mechanism of TCM. This study aimed at determining the optimal compatibility and possible mechanism of Achyranthes bidentate for OA treatment. Results showed that the adhesion score of the OA group is higher than NC group, and showed a trend of down-regulation in the intervention group. The CHI3L1 and IL-1β in joint fluid of the OA group was significantly increased compared to the sham operation group (NC group). Group G, I, and L exhibited significantly down-regulated CHI3L1, while groups C, F, I, K, and L exhibited reduced IL-1β. Joint adhesion, damage in cartilage, and synovial tissue was found in the OA model, cartilage tissue was found recovered in groups I, J, and L, and synovial tissue was recovered in group G, I, and L. Thus, group I and L were chosen for metabolite analysis, and indole-3-propionic acid was slightly up-regulated, while koeiginequinone A, prostaglandin H2, and 1-hydroxy-3-methoxy-10-methylacridonew were down-regulated in group I and L. According to functional analysis, the arachidonic acid (AA) metabolic pathway is enriched. Down-regulated expression of vital proteins in the AA metabolism pathway, such as PGE2 and COX2 in group I and L were verified. In conclusion, Hydroxyecdysone, Oleanolic acid, Achyranthes bidentata polysaccharide at a compatibility of 0.03-μg/mg, 2.0-μg/mg, 20.0-μg/mg or 0.03-μg/mg, 2.0-μg/mg, 10.0-μg/mg, respectively, may be the optimal compatibility of Achyranthes bidentate.

Crosstalk Among circRNA/lncRNA, miRNA, and mRNA in Osteoarthritis

Osteoarthritis (OA) is a joint disease that is pervasive in life, and the incidence and mortality of OA are increasing, causing many adverse effects on people's life. Therefore, it is very vital to identify new biomarkers and therapeutic targets in the clinical diagnosis and treatment of OA. ncRNA is a nonprotein-coding RNA that does not translate into proteins but participates in protein translation. At the RNA level, it can perform biological functions. Many studies have found that miRNA, lncRNA, and circRNA are closely related to the course of OA and play important regulatory roles in transcription, post-transcription, and post-translation, which can be used as biological targets for the prevention, diagnosis, and treatment of OA. In this review, we summarized and described the various roles of different types of miRNA, lncRNA, and circRNA in OA, the roles of different lncRNA/circRNA-miRNA-mRNA axis in OA, and the possible prospects of these ncRNAs in clinical application.

Shared genetic liability between major depressive disorder and osteoarthritis

Aims

Deciphering the genetic relationships between major depressive disorder (MDD) and osteoarthritis (OA) may facilitate an understanding of their biological mechanisms, as well as inform more effective treatment regimens. We aim to investigate the mechanisms underlying relationships between MDD and OA in the context of common genetic variations.

Methods

Linkage disequilibrium score regression was used to test the genetic correlation between MDD and OA. Polygenic analysis was performed to estimate shared genetic variations between the two diseases. Two-sample bidirectional Mendelian randomization analysis was used to investigate causal relationships between MDD and OA. Genomic loci shared between MDD and OA were identified using cross-trait meta-analysis. Fine-mapping of transcriptome-wide associations was used to prioritize putatively causal genes for the two diseases.

Results

MDD has a significant genetic correlation with OA (r_g = 0.29) and the two diseases share a considerable proportion of causal variants. Mendelian randomization analysis indicates that genetic liability to MDD has a causal effect on OA (b_xy = 0.24) and genetic liability to OA conferred a causal effect on MDD (b_xy = 0.20). Cross-trait meta-analyses identified 29 shared genomic loci between MDD and OA. Together with fine-mapping of transcriptome-wide association signals, our results suggest that Estrogen Receptor 1 (ESR1), SRY-Box Transcription Factor 5 (SOX5), and Glutathione Peroxidase 1 (GPX1) may have therapeutic implications for both MDD and OA.

Conclusion

The study reveals substantial shared genetic liability between MDD and OA, which may confer risk for one another. Our findings provide a novel insight into phenotypic relationships between MDD and OA.

Cite this article: Bone Joint Res 2022;11(1):12–22.

Recent Advances in Pharmacological Intervention of Osteoarthritis: A Biological Aspect

Osteoarthritis (OA) is a degenerative joint disease in the musculoskeletal system with a relatively high incidence and disability rate in the elderly. It is characterized by the degradation of articular cartilage, inflammation of the synovial membrane, and abnormal structure in the periarticular and subchondral bones. Although progress has been made in uncovering the molecular mechanism, the etiology of OA is still complicated and unclear. Nevertheless, there is no treatment method that can effectively prevent or reverse the deterioration of cartilage and bone structure. In recent years, in the field of pharmacology, research focus has shifted to disease prevention and early treatment rather than disease modification in OA. Biologic agents become more and more attractive as their direct or indirect intervention effects on the initiation or development of OA. In this review, we will discuss a wide spectrum of biologic agents ranging from DNA, noncoding RNA, exosome, platelet-rich plasma (PRP), to protein. We searched for key words such as OA, DNA, gene, RNA, exosome, PRP, protein, and so on. From the pharmacological aspect, stem cell therapy is a very special technique, which is not included in this review. The literatures ranging from January 2016 to August 2021 were included and summarized. In this review, we aim to help readers have a complete and precise understanding of the current pharmacological research progress in the intervention of OA from the biological aspect and provide an indication for the future translational studies.

Whole-Transcriptome Sequence of Degenerative Meniscus Cells Unveiling Diagnostic Markers and Therapeutic Targets for Osteoarthritis

Meniscus plays an important role in joint homeostasis. Tear or degeneration of meniscus might facilitate the process of knee osteoarthritis (OA). Hence, to investigate the transcriptome change during meniscus degeneration, we reveal the alterations of messenger RNA (mRNA), microRNA (miRNA), long noncoding RNA (lncRNA), and circular RNA (circRNA) in meniscus during OA by whole-transcriptome sequence. A total of 375 mRNAs, 15 miRNAs, 56 lncRNAs, and 90 circRNAs were significantly altered in the degenerative meniscus treated with interleukin-1β (IL-1β). More importantly, highly specific co-expression RNA (ceRNA) networks regulated by lncRNA LOC107986251-miR-212-5p-SESN3 and hsa_circ_0018069-miR-147b-3p-TJP2 were screened out during IL-induced meniscus degeneration, unveiling potential therapeutic targets for meniscus degeneration during the OA process. Furthermore, lipocalin-2 (LCN2) and RAB27B were identified as potential biomarkers in meniscus degeneration by overlapping three previously constructed databases of OA menisci. LCN2 and RAB27B were both upregulated in osteoarthritic menisci and IL-1β-treated menisci and were highly associated with the severity of OA. This could introduce potential novel molecules into the database of clinical diagnostic biomarkers and possible therapeutic targets for early-stage OA treatment.

Transcriptomic and epigenomic analyses uncovered Lrrc15 as a contributing factor to cartilage damage in osteoarthritis

In osteoarthritis (OA), articular chondrocytes display phenotypic and functional changes associated with epigenomic alterations. These changes contribute to the disease progression, which is characterized by dysregulated reparative processes and abnormal extracellular matrix remodeling leading to cartilage degradation. Recent studies using a murine model of posttraumatic OA highlighted the contribution of changes in DNA hydroxymethylation (5hmC) to OA progression. Here, we integrated transcriptomic and epigenomic analyses in cartilage after induction of OA to show that the structural progression of OA is accompanied by early transcriptomic and pronounced DNA methylation (5mC) changes in chondrocytes. These changes accumulate over time and are associated with recapitulation of developmental processes, including cartilage development, chondrocyte hypertrophy, and ossification. Our integrative analyses also uncovered that Lrrc15 is differentially methylated and expressed in OA cartilage, and that it may contribute to the functional and phenotypic alterations of chondrocytes, likely coordinating stress responses and dysregulated extracellular matrix remodeling.

MiR-4303 relieves chondrocyte inflammation by targeting ASPN in osteoarthritis

Background

Osteoarthritis (OA) is a severe articular cartilage disease whose pathogenesis involves the inflammation of chondrocytes. MicroRNAs (miRNAs) are considered to be effective inflammation regulators. However, the regulatory mechanism of miRNAs in osteoarthritis needs to be further elucidated. In this paper, we aim to investigate the underlying mechanisms by which miR-4303 regulates osteoarthritis.

Methods

RT-qPCR is performed to detect the mRNA expression levels of miR-4303, ASPN, PDIA3, PIK3CA, and TRAF3. CCK-8 assay and EdU assay are carried to assess chondrocyte viability. The protein expression levels of ASPN, PCNA, Ki-67, CyclinA1, CyclinB1, CyclinD2, p27, Bax, Bcl-2, cleaved caspase-3, and Cleaved caspase-9 were measured by western blot. FACs is performed to detect the cell cycle and apoptosis of chondrocyte. ELISA is conducted to assess the levels of TNF-β, IL-1β and IL-6 in the supernatant of chondrocytes. The potential binding sites of miR-4303 and ASPN are predicted by the miRDB database and confirmed by the dual-luciferase reporter gene assay.

Results

Our findings illustrated that miR-4303 was down-regulated in arthritic tissues and LPS-induced chondrocytes; miR-4303 overexpression rescued the decrease in cell viability, cell cycle arrest and apoptosis induced by LPS. Furthermore, miR-4303 overexpression inhibited the release of inflammatory factors in LPS-induced chondrocytes, miR-4303 relieved chondrocyte inflammation via targeting ASPN.

Conclusion

MiR-4303 serves as a prognostic biomarker and relieves chondrocyte inflammation via targeting ASPN. Our findings provide novel prognostic biomarkers in predicting the progression and prognosis of osteoarthritis.

Role of circular RNAs in osteoarthritis (Review)

Osteoarthritis (OA) is a chronic bone and joint disease characterized by articular cartilage degeneration and joint inflammation. OA is the most common form of arthritis, and the major clinical manifestations of OA are chronic pain and joint activity disorder, which severely affect patients' quality of life. Circular RNA (circRNA) is a type of non-coding RNA that is ubiquitous in eukaryotic cells. Unlike standard linear RNAs, they form a covalently closed continuous loop without 5' or 3' polarity. They are usually considered as byproducts of mis-splicing or mRNA processing. CircRNAs have been detected and identified in numerous species. Various studies have confirmed that certain circRNAs are differentially expressed in OA cartilage and are closely associated with a variety of pathological processes of OA, including extracellular matrix degradation, inflammation and apoptosis. The present study reviewed the latest research on circRNAs in the pathogenesis of OA, providing a novel direction for the prevention, diagnosis and treatment of OA.

Single-cell RNA-sequencing atlas of bovine caudal intervertebral discs: Discovery of heterogeneous cell populations with distinct roles in homeostasis

Back and neck pain are significant healthcare burdens that are commonly associated with pathologies of the intervertebral disc (IVD). The poor understanding of the cellular heterogeneity within the IVD makes it difficult to develop regenerative IVD therapies. To address this gap, we developed an atlas of bovine (Bos taurus) caudal IVDs using single-cell RNA-sequencing (scRNA-seq). Unsupervised clustering resolved 15 unique clusters, which we grouped into the following annotated partitions: nucleus pulposus (NP), outer annulus fibrosus (oAF), inner AF (iAF), notochord, muscle, endothelial, and immune cells. Analyzing the pooled gene expression profiles of the NP, oAF, and iAF partitions allowed us to identify novel markers for NP (CP, S100B, H2AC18, SNORC, CRELD2, PDIA4, DNAJC3, CHCHD7, and RCN2), oAF (IGFBP6, CTSK, LGALS1, and CCN3), and iAF (MGP, COMP, SPP1, GSN, SOD2, DCN, FN1, TIMP3, WDR73, and GAL) cells. Network analysis on subpopulations of NP and oAF cells determined that clusters NP1, NP2, NP4, and oAF1 displayed gene expression profiles consistent with cell survival, suggesting these clusters may uniquely support viability under the physiological stresses of the IVD. Clusters NP3, NP5, oAF2, and oAF3 expressed various extracellular matrix (ECM)-associated genes, suggesting their role in maintaining IVD structure. Lastly, transcriptional entropy and pseudotime analyses found that clusters NP3 and NP1 had the most stem-like gene expression signatures of the NP partition, implying these clusters may contain IVD progenitor cells. Overall, results highlight cell type diversity within the IVD, and these novel cell phenotypes may enhance our understanding of IVD development, homeostasis, degeneration, and regeneration.

Intra-articular low-dose parathyroid hormone (1-34) improves mobility and articular cartilage quality in a preclinical age-related knee osteoarthritis model

Aims

Osteoarthritis (OA) is prevalent among the elderly and incurable. Intra-articular parathyroid hormone (PTH) ameliorated OA in papain-induced and anterior cruciate ligament transection-induced OA models; therefore, we hypothesized that PTH improved OA in a preclinical age-related OA model.

Methods

Guinea pigs aged between six and seven months of age were randomized into control or treatment groups. Three- or four-month-old guinea pigs served as the young control group. The knees were administered 40 μl intra-articular injections of 10 nM PTH or vehicle once a week for three months. Their endurance as determined from time on the treadmill was evaluated before kill. Their tibial plateaus were analyzed using microcalculated tomography (μCT) and histological studies.

Results

PTH increased the endurance on the treadmill test, preserved glycosaminoglycans, and reduced Osteoarthritis Research Society International score and chondrocyte apoptosis rate. No difference was observed in the subchondral plate bone density or metaphyseal trabecular bone volume and bone morphogenetic 2 protein staining.

Conclusion

Subchondral bone is crucial in the initiation and progression of OA. Although previous studies have shown that subcutaneous PTH alleviates knee OA by improving subchondral and metaphyseal bone mass, we demonstrated that intra-articular PTH injections improved spontaneous OA by directly affecting the cartilage rather than the subchondral or metaphyseal bone in a preclinical age-related OA model.

Cite this article: Bone Joint Res 2021;10(8):514–525.

PR11-364P22.2/ATF3 protein interaction mediates IL-1β-induced catabolic effects in cartilage tissue and chondrocytes

Osteoarthritis (OA) is a degenerative joint disease which lacks effective medical treatment due to ill-defined molecular mechanisms underlying the pathology. Inflammation is a key factor that induces and aggravates OA. Therefore, the current study aims to explore roles of the dysregulated long non-coding RNAs in the pro-inflammatory cytokine IL-1β-mediated catabolic effects in cartilage tissue and chondrocytes. We identified RP11-364P22.2 as dysregulated in OA patient-derived cartilage tissues and highly responsive to IL-1β stimulus. RNA pull-down coupled with mass spectrometry demonstrated that RP11-364P22.2 physically binds to activating transcription factor 3 (ATF3) and thus increases the protein stability and facilitates its nuclear translocation. Loss- and gain-of-function assays indicated that the interaction between RP11-364P22.2 and ATF3 is indispensable for the detrimental effects of IL-1β including growth inhibition, apoptosis induction as well as degradation of the key chondrocyte structural proteins of type II collage and Aggrecan and synthesis of the extracellular matrix-degrading enzyme MMP13 in chondrocytes. In vivo, depletion of the RP11-364P22.2 effector ATF3 drastically prevented OA development in the rats with surgical destabilization of the medial meniscus (DMM). These results highlight the important roles of lncRNAs in the pathogenesis of OA and indicate the RP11-364P22.2/ATF3 regulatory axis as a potential therapeutic target of inflammation-induced OA.

Clinical evaluation after matrix-associated autologous chondrocyte transplantation : a comparison of four different graft types

Aims

The aim of this retrospective study was to determine if there are differences in short-term clinical outcomes among four different types of matrix-associated autologous chondrocyte transplantation (MACT).

Methods

A total of 88 patients (mean age 34 years (SD 10.03), mean BMI 25 kg/m² (SD 3.51)) with full-thickness chondral lesions of the tibiofemoral joint who underwent MACT were included in this study. Clinical examinations were performed preoperatively and 24 months after transplantation. Clinical outcomes were evaluated using the International Knee Documentation Committee (IKDC) Subjective Knee Form, the Brittberg score, the Tegner Activity Scale, and the visual analogue scale (VAS) for pain. The Kruskal-Wallis test by ranks was used to compare the clinical scores of the different transplant types.

Results

The mean defect size of the tibiofemoral joint compartment was 4.28 cm² (SD 1.70). In total, 11 patients (12.6%) underwent transplantation with Chondro-Gide (matrix-associated autologous chondrocyte implantation (MACI)), 40 patients (46.0%) with Hyalograft C (HYAFF), 21 patients (24.1%) with Cartilage Regeneration System (CaReS), and 15 patients (17.2%) with NOVOCART 3D. The mean IKDC Subjective Knee Form score improved from 35.71 (SD 6.44) preoperatively to 75.26 (SD 18.36) after 24 months postoperatively in the Hyalograft group, from 35.94 (SD 10.29) to 71.57 (SD 16.31) in the Chondro-Gide (MACI) group, from 37.06 (SD 5.42) to 71.49 (SD 6.76) in the NOVOCART 3D group, and from 45.05 (SD 15.83) to 70.33 (SD 19.65) in the CaReS group. Similar improvements were observed in the VAS and Brittberg scores.

Conclusion

Two years postoperatively, there were no significant differences in terms of outcomes. Our data demonstrated that MACT, regardless of the implants used, resulted in good clinical improvement two years after transplantation for localized tibiofemoral defects.

Cite this article: Bone Joint Res 2021;10(7):370–379.

MicroRNA-486 promotes a more catabolic phenotype in chondrocyte-like cells by targeting SIRT6 : possible involvement in cartilage degradation in osteoarthritis

AIMS

Osteoarthritis (OA) is characterized by persistent destruction of articular cartilage. It has been found that microRNAs (miRNAs) are closely related to the occurrence and development of OA. The purpose of the present study was to investigate the mechanism of miR-486 in the development and progression of OA.

METHODS

The expression levels of miR-486 in cartilage were determined by quantitative real-time polymerase chain reaction (qRT-PCR). The expression of collagen, type II, alpha 1 (COL2A1), aggrecan (ACAN), matrix metalloproteinase (MMP)-13, and a disintegrin and metalloproteinase with thrombospondin motifs-4 (ADAMTS4) in SW1353 cells at both messenger RNA (mRNA) and protein levels was determined by qRT-PCR, western blot, and enzyme-linked immunosorbent assay (ELISA). Double luciferase reporter gene assay, qRT-PCR, and western blot assay were used to determine whether silencing information regulator 6 (SIRT6) was involved in miR-486 induction of chondrocyte-like cells to a more catabolic phenotype.

RESULTS

Compared with osteonecrosis, the expression of miR-486 was significantly upregulated in cartilage from subjects with severe OA. In addition, overexpressed miR-486 promoted a catabolic phenotype in SW1353 cells by upregulating the expressions of ADAMTS4 and MMP-13 and down-regulating the expressions of COL2A1 and ACAN. Conversely, inhibition of miR-486 had the opposite effect. Furthermore, overexpression of miR-486 significantly inhibited the expression of SIRT6, confirming that SIRT6 is a direct target of miR-486. Moreover, SW1353 cells were transfected with small interfering RNA (si)-SIRT6 and it was found that SIRT6 was involved in and inhibited miR-486-induced changes to SW1353 gene expression.

CONCLUSION

Our results indicate that miR-486 promotes a catabolic phenotype in SW1353 cells in OA by targeting SIRT6. Our findings might provide a potential therapeutic target and theoretical basis for OA. Cite this article: Bone Joint Res 2021;10(7):459-466.

PARP-1 inhibition attenuates the inflammatory response in the cartilage of a rat model of osteoarthritis

Aims

Poly (ADP-ribose) polymerase (PARP) inhibitor has been reported to attenuate inflammatory response in rat models of inflammation. This study was designed to investigate the effect of PARP signalling in osteoarthritis (OA) cartilage inflammatory response in an OA rat model.

Methods

The OA model was established by anterior cruciate ligament transection with medial meniscectomy in Wistar rats. The poly (ADP-ribose) polymerase 1 (PARP-1) shRNA (short hairpin (sh)-PARP-1) and negative control shRNA (sh-NC) were delivered using a lentiviral vector and were intra-articularly injected into rats after surgery. The weight-bearing distribution of the hind limbs and the knee joint width were measured every two weeks. The expression levels of PARP-1, inducible nitric oxide synthase (iNOS), and cyclooxygenase-2 (COX-2) in cartilage were determined using real-time quantitative reverse transcription polymerase chain reaction (RT-qPCR) and Western blot. The serum concentrations of inflammatory cytokines were detected using enzyme-linked immunosorbent assay (ELISA).

Results

PARP-1 expression level significantly increased in the cartilage of the established OA rat model. sh-PARP-1 treatment suppressed PARP-1 levels, decreased the Δ Force (the difference between the weight on ipsilateral limb and contralateral limb) and the knee joint width, inhibited cartilage matrix catabolic enzymes, and ameliorated OA cartilage degradation and attenuated inflammatory response.

Conclusion

PARP-1 inhibition attenuates OA cartilage inflammatory response in the OA rat model.

Cite this article: Bone Joint Res 2021;10(7):401–410.

Exhausting circ_0136474 and Restoring miR-766-3p Attenuate Chondrocyte Oxidative Injury in IL-1β-Induced Osteoarthritis Progression Through Regulating DNMT3A

Circular RNA circ_0136474 is a new contributor of human osteoarthritis (OA) by suppressing chondrocyte proliferation. However, its role and mechanism in OA chondrocyte injury remain ill defined. Herein, we performed real-time quantitative PCR to detect RNA expression of circ_0136474, microRNA (miR)-766-3p, and DNA methyltransferase 3A (DNMT3A) and utilized Western blotting to measure protein expression of DNMT3A, matrix metalloproteinase-1 (MMP1), MMP13, collagen II, proliferating cell nuclear antigen (PCNA) and B cell lymphoma (Bcl)-2, and Bcl-2-associated X protein (Bax). Direct interaction between miR-766-3p and circ_0136474 or DNMT3A was confirmed by bioinformatics algorithms, dual-luciferase reporter assay, and RNA immunoprecipitation. Functional experiments including cell counting kit-8 assay, flow cytometry, and special assay kits were employed to measure oxidative injury in interleukin (IL)-1β-induced OA-like chondrocytes. First, IL-1β administration induced cell viability inhibition, collagen II suppression, and promotion of MMP1 and MMP13 in human chondrocyte CHON-001 cells. Expression of circ_0136474 and DNMT3A was upregulated, and miR-766-3p was downregulated in human OA cartilages and IL-1β-induced CHON-001 cells. Functionally, both blocking circ_0136474 and upregulating miR-766-3p could rescue cell viability and levels of PCNA, Bcl-2, reduced glutathione (GSH), and total superoxide dismutase (SOD), and attenuate apoptosis rate and levels of Bax, reactive oxygen species (ROS), and lipid peroxidation malondialdehyde (MDA). Mechanically, circ_0136474 served as miR-766-3p sponge to govern miR-766-3p-targeted DNMT3A expression. Accidently, restoring DNMT3A counteracted the miR-766-3p upregulation role, and silencing miR-766-3p weakened circ_0136474 knockdown effect in IL-1β-induced CHON-001 cells. In conclusion, exhausting circ_0136474 could mitigate OA chondrocyte oxidative injury through regulating miR-766-3p/DNMT3A axis.

Identification of key regulators responsible for dysregulated networks in osteoarthritis by large-scale expression analysis

BACKGROUND

Osteoarthritis (OA) is a worldwide musculoskeletal disorder. However, disease-modifying therapies for OA are not available. Here, we aimed to characterize the molecular signatures of OA and to identify novel therapeutic targets and strategies to improve the treatment of OA.

METHODS

We collected genome-wide transcriptome data performed on 132 OA and 74 normal human cartilage or synovium tissues from 7 independent datasets. Differential gene expression analysis and functional enrichment were performed to identify genes and pathways that were dysregulated in OA. The computational drug repurposing method was used to uncover drugs that could be repurposed to treat OA.

RESULTS

We identified several pathways associated with the development of OA, such as extracellular matrix organization, inflammation, bone development, and ossification. By protein-protein interaction (PPI) network analysis, we prioritized several hub genes, such as JUN, CDKN1A, VEGFA, and FOXO3. Moreover, we repurposed several FDA-approved drugs, such as cardiac glycosides, that could be used in the treatment of OA.

CONCLUSIONS

We proposed that the hub genes we identified would play a role in cartilage homeostasis and could be important diagnostic and therapeutic targets. Drugs such as cardiac glycosides provided new possibilities for the treatment of OA.

Combining bulk and single-cell RNA-sequencing data to reveal gene expression pattern of chondrocytes in the osteoarthritic knee

Osteoarthritis (OA) occurs mostly in the knees, hips, finger interphalangeal joints, and spinal facet joints, and is characterized by cartilage degeneration. The existing bulk RNA sequencing (bulk RNA-seq) and single-cell sequencing (scRNA-seq) data for chondrocytes in the osteoarthritic knee joint provide the expression profiles of entire cell populations and individual cells, respectively. Here, we aimed to analyze these two types of sequencing data in order to obtain a more comprehensive understanding of OA. We compared the analysis results of bulk RNA-seq and scRNA-seq from the dataset GSE114007 and the dataset GSE104782, respectively, and identified the differentially expressed genes (DEGs). Then, we tried to find the key The transcription factor is a more fomal term (TFs) and long non-coding RNA (lncRNA) regulation. We highlighted 271 genes that were simultaneously suggested by these two types of data and provided their possible expression pattern in OA. Among the 271 genes, we identified 14 TFs, and TWIST2, MYBL2, RELA, JUN, KLF4, and PTTG1 could be the key TFs for the 271 genes. We also found that 8 lncRNAs among the 271 genes and the lncRNA regulation between CYTOR and NRP1 could contribute to the pain and vascularization of cartilage in the osteoarthritic knee. In short, our research combined the analysis results of bulk RNA-seq and scRNA-seq data for OA chondrocytes, which will contribute to further elucidation of the molecular mechanisms of OA pathogenesis.

The function of lncRNAs in the pathogenesis of osteoarthritis

Osteoarthritis (OA), one of the most common motor system disorders, is a degenerative disease involving progressive joint destruction caused by a variety of factors. At present, OA has become the fourth most common cause of disability in the world. However, the pathogenesis of OA is complex and has not yet been clarified. Long non-coding RNA (lncRNA) refers to a group of RNAs more than 200 nucleotides in length with limited protein-coding potential, which have a wide range of biological functions including regulating transcriptional patterns and protein activity, as well as binding to form endogenous small interference RNAs (siRNAs) and natural microRNA (miRNA) molecular sponges. In recent years, a large number of lncRNAs have been found to be differentially expressed in a variety of pathological processes of OA, including extracellular matrix (ECM) degradation, synovial inflammation, chondrocyte apoptosis, and angiogenesis. Obviously, lncRNAs play important roles in regulating gene expression, maintaining the phenotype of cartilage and synovial cells, and the stability of the intra-articular environment. This article reviews the results of the latest research into the role of lncRNAs in a variety of pathological processes of OA, in order to provide a new direction for the study of OA pathogenesis and a new target for prevention and treatment. Cite this article: Bone Joint Res 2021;10(2):122-133.

Current Models for Development of Disease-Modifying Osteoarthritis Drugs

Osteoarthritis (OA) is a painful and disabling disease that affects millions of people worldwide. Symptom-alleviating treatments exist, although none with long-term efficacy. Furthermore, there are currently no disease-modifying OA drugs (DMOADs) with demonstrated efficacy in OA patients, which is, in part, attributed to a lack of full understanding of the pathogenesis of OA. The inability to translate findings from basic research to clinical applications also highlights the deficiencies in the available OA models at simulating the clinically relevant pathologies and responses to treatments in humans. In this review, the current status in the development of DMOADs will be first presented, with special attention to those in Phase II-IV clinical trials. Next, current in vitro, ex vivo, and in vivo OA models are summarized and the respective advantages and disadvantages of each are highlighted. Of note, the development and application of microphysiological or tissue-on-a-chip systems for modeling OA in humans are presented and the issues that need to be addressed in the future are discussed. Microphysiological systems should be given serious consideration for their inclusion in the DMOAD development pipeline, both for their ability to predict drug safety and efficacy in human clinical trials at present, as well as for their potential to serve as a test platform for personalized medicine. Impact statement At present, no disease-modifying osteoarthritis (OA) drugs (DMOADs) have been approved for widespread clinical use by regulatory bodies. The failure of developing effective DMOADs is likely owing to multiple factors, not the least of which are the intrinsic differences between the intact human knee joint and the preclinical models. This work summarizes the current OA models for the development of DMOADs, discusses the advantages/disadvantages of each, and then proposes future model development to aid in the discovery of effective and personalized DMOADs. The review also highlights the microphysiological systems, which are emerging as a new platform for drug development.

Understanding osteoarthritis pathogenesis: a multiomics system-based approach

PURPOSE OF REVIEW

Osteoarthritis is a heterogeneous, multifactorial condition regulated by complex biological interactions at multiple levels. Comprehensive understanding of these regulatory interactions is required to develop feasible advances to improve patient outcomes. Improvements in technology have made extensive genomic, transcriptomic, epigenomic, proteomic, and metabolomic profiling possible. This review summarizes findings over the past 20 months related to omics technologies in osteoarthritis and examines how using a multiomics approach is necessary for advancing our understanding of osteoarthritis as a disease to improve precision osteoarthritis treatments.

RECENT FINDINGS

Using the search terms 'genomics' or 'transcriptomics' or 'epigenomics' or 'proteomics' or 'metabolomics' and 'osteoarthritis' from January 1, 2018 to August 31, 2019, we identified advances in omics approaches applied to osteoarthritis. Trends include untargeted whole genome, transcriptome, proteome, and metabolome analyses leading to identification of novel molecular signatures, cell subpopulations and multiomics validation approaches.

SUMMARY

To address the complexity of osteoarthritis, integration of multitissue analyses by multiomics approaches with the inclusion of longitudinal clinical data is necessary for a comprehensive understanding of the disease process, and for appropriate development of efficacious diagnostics, prognostics, and biotherapeutics.

Circ0083429 Regulates Osteoarthritis Progression via the Mir-346/SMAD3 Axis

Osteoarthritis (OA) is a degenerative joint disease. Currently, apart from symptomatic treatment or joint replacement, no other effective treatments for OA exist. The mechanisms underlying OA remain elusive and require further research. Circular RNAs (circRNAs) are known to be involved in many diseases; however, their function in OA is not yet fully understood. Here, we identified a novel circRNA, Circ0083429. The role of Circ0083429 in OA was confirmed via western blot (WB), quantitative real-time PCR (qRT-PCR), and immunofluorescence (IF) through knockdown and overexpression experiments. The binding of Circ0083429 to downstream miR-346 and its target gene SMAD3 was predicted via bioinformatics analysis and verified using a luciferase reporter assay and RNA pulldown experiments. Finally, the function of Circ0083429 was evaluated in mouse OA models. In our study, we found that Circ0083429 regulates the homeostasis of the extracellular matrix (ECM) in human chondrocytes. Mechanistically, Circ0083429 affects OA by regulating the mRNA level of SMAD3 through the sponging of microRNA (miRNA)-346. Injecting adeno-associated virus Circ0083429 into the intra-junction of the mouse knee alleviated OA. In conclusion, Circ0083429 regulates the ECM via the regulation of the downstream miRNA-346/SMAD3 in human chondrocytes, which provides a new therapeutic strategy for OA.

Circular RNAs in osteoarthritis: indispensable regulators and novel strategies in clinical implications

Over the past decades, circular RNAs (circRNAs) have emerged as a hot spot and sparked intensive interest. Initially considered as the transcriptional noises, further studies have indicated that circRNAs are crucial regulators in multiple cellular biological processes, and thus engage in the development and progression of many diseases including osteoarthritis (OA). OA is a prevalent disease that mainly affects those aging, obese and post-traumatic population, posing as a major source of socioeconomic burden. Recently, numerous circRNAs have been found aberrantly expressed in OA tissues compared with counterparts. More importantly, circRNAs have been demonstrated to interplay with components in OA microenvironments, such as chondrocytes, synoviocytes and macrophages, by regulation of their proliferation, apoptosis, autophagy, inflammation, or extracellular matrix reorganization. Herein, in this review, we extensively summarize the roles of circRNAs in OA microenvironment, progression, and putative treatment, as well as envision the future directions for circRNAs research in OA, with the aim to provide a novel insight into this field.

Dysregulation of USP18/FTO/PYCR1 signaling network promotes bladder cancer development and progression

N6-methyladenosine refers to a methylation of adenosine base at the 6th nitrogen position, which is the dominant methylation modification in both message and non-coding RNAs. Dysregulation of RNA m6A methylation causes tumorigenesis in humans. The key N6-methyladenosine demethylase fat-mass and obesity-associated protein (FTO) is negatively correlated with the overall survival of bladder cancer patients, but the underlying mechanism remains poorly understood. In this study, we demonstrated that the post-translational deubiquitination by USP18 up-regulates the protein but not mRNA of FTO in bladder cancer tissues and cells. As a result, FTO decreased N6-methyladenosine methylation level in PYCR1 through its demethylase enzymatic activity and stabilized PYCR1 transcript to promote bladder cancer initiation and progression. Our work shows the importance of N6-methyladenosine RNA modification in bladder cancer development, and highlights UPS18/FTO/PYCR1 signaling network as potential therapeutic targets of bladder cancer.

Role of adipose tissues in osteoarthritis

PURPOSE OF REVIEW

Epidemiologic studies reveal that the link between obesity and osteoarthritis cannot be uniquely explained by overweight-associated mechanical overload. For this reason, much attention focuses on the endocrine activity of adipose tissues. In addition to the systemic role of visceral and subcutaneous adipose tissues, many arguments highlight the involvement of local adipose tissues in osteoarthritis.

RECENT FINDINGS

Alteration in MRI signal intensity of the infrapatellar fat pad may predict both accelerated knee osteoarthritis and joint replacement. In this context, recent studies show that mesenchymal stromal cells could play a pivotal role in the pathological remodelling of intra-articular adipose tissues (IAATs) in osteoarthritis. In parallel, recent findings underline bone marrow adipose tissue as a major player in the control of the bone microenvironment, suggesting its possible role in osteoarthritis.

SUMMARY

The recent description of adipose tissues of various phenotypes within an osteoarthritic joint allows us to evoke their direct involvement in the initiation and progression of the osteoarthritic process. We can expect in the near future the discovery of novel molecules targeting these tissues.

Circular RNA CSNK1G1 promotes the progression of osteoarthritis by targeting the miR‑4428/FUT2 axis

Osteoarthritis (OA) is a chronic disease that results in chronic arthralgia and functional disability of the affected joint. To date, there is no effective treatment available for this disease. Circular RNAs (circRNAs) are a type of intracellular stable RNA that can regulate the development and progression of OA. However, the function of circCSNK1G1 in OA has not yet been investigated. In the present study, it was found that circCSNK1G1 was upregulated in OA cartilage tissues. The upregulation of circCSNK1G1 was associated with extracellular matrix (ECM) degradation and chondrocyte apoptosis. Moreover, the expression of miR‑4428 was downregulated and that of fucosyltransferase 2 (FUT2) was upregulated in OA‑affected cartilage tissues. Dual‑luciferase reporter assay and RNA immunoprecipitation confirmed that miR‑4428 targeted FUT2 mRNA to inhibit FUT2 expression. circCSNK1G1 and FUT2 induced ECM degradation and chondrocyte apoptosis. The negative effects of circCSNK1G1 and FUT2 were reversed by miR‑4428. On the whole, the present study demonstrates that circCSNK1G1 promotes the development of OA by targeting the miR‑4428/FUT2 axis. Thus, the circCSNK1G1/miR‑4428/FUT2 axis may present a novel target for the treatment of OA in the clinical setting.

Intra-Articular Injection of (-)-Epigallocatechin 3-Gallate to Attenuate Articular Cartilage Degeneration by Enhancing Autophagy in a Post-Traumatic Osteoarthritis Rat Model

(-)-Epigallocatechin 3-gallate (EGCG) is the main active green tea catechin and has a wide variety of benefits for health. Post-traumatic osteoarthritis (PTOA) occurs as a consequence of joint injuries that commonly happen in the young population. In this study, we investigated the effects of EGCG on PTOA prevention by using the anterior cruciate ligament transection (ACLT)–OA model and further investigated the roles of autophagy in OA treatment. Our results showed that intra-articular injection of EGCG significantly improved the functional performances and decreased cartilage degradation. EGCG treatment attenuated the inflammation on synovial tissue and cartilage through less immunostained cyclooxygenase-2 and matrix metalloproteinase-13. We further noted EGCG may modulate the chondrocyte apoptosis by activation of the cytoprotective autophagy through reducing the expression of the mTOR and enhancing the expression of microtubule-associated protein light chain 3, beclin-1, and p62. In conclusion, intra-articular injection of EGCG after ACL injury inhibited the joint inflammation and cartilage degradation, thereby increasing joint function. EGCG treatment also reduced the chondrocyte apoptosis, possibly by activating autophagy. These findings suggested that EGCG may be a potential disease-modifying drug for preventing OA progression.

Synovial fluid peptidase activity as a biomarker for knee osteoarthritis clinical progression

AIMS

To analyze the potential role of synovial fluid peptidase activity as a measure of disease burden and predictive biomarker of progression in knee osteoarthritis (KOA).

METHODS

A cross-sectional study of 39 patients (women 71.8%, men 28.2%; mean age of 72.03 years (SD 1.15) with advanced KOA (Ahlbäck grade ≥ 3 and clinical indications for arthrocentesis) recruited through the (Orthopaedic Department at the Complejo Asistencial Universitario de León, Spain (CAULE)), measuring synovial fluid levels of puromycin-sensitive aminopeptidase (PSA), neutral aminopeptidase (NAP), aminopeptidase B (APB), prolyl endopeptidase (PEP), aspartate aminopeptidase (ASP), glutamyl aminopeptidase (GLU) and pyroglutamyl aminopeptidase (PGAP).

RESULTS

Synovial fluid peptidase activity varied significantly as a function of clinical signs, with differences in levels of PEP (p = 0.020), ASP (p < 0.001), and PGAP (p = 0. 003) associated with knee locking, PEP (p = 0.006), ASP (p = 0.001), GLU (p = 0.037), and PGAP (p = 0.000) with knee failure, and PEP (p = 0.006), ASP (p = 0.001), GLU (p = 0.037), and PGAP (p < 0.001) with knee effusion. Further, patients with the greatest functional impairment had significantly higher levels of APB (p = 0.005), PEP (p = 0.005), ASP (p = 0.006), GLU (p = 0.020), and PGAP (p < 0.001) activity, though not of NAP or PSA, indicating local alterations in the renin-angiotensin system. A binary logistic regression model showed that PSA was protective (p = 0.005; Exp (B) 0.949), whereas PEP (p = 0.005) and GLU were risk factors (p = 0.012).

CONCLUSION

These results suggest synovial fluid peptidase activity could play a role as a measure of disease burden and predictive biomarker of progression in KOA. Cite this article: Bone Joint Res 2020;9(11):789-797.

Knockout of miR-21-5p alleviates cartilage matrix degradation by targeting Gdf5 in temporomandibular joint osteoarthritis

AIMS

The study aimed to determine whether the microRNA miR21-5p (MiR21) mediates temporomandibular joint osteoarthritis (TMJ-OA) by targeting growth differentiation factor 5 (Gdf5).

METHODS

TMJ-OA was induced in MiR21 knockout (KO) mice and wild-type (WT) mice by a unilateral anterior crossbite (UAC) procedure. Mouse tissues exhibited histopathological changes, as assessed by: Safranin O, toluidine blue, and immunohistochemistry staining; western blotting (WB); and quantitative real-time polymerase chain reaction (RT-qPCR). Mouse condylar chondrocytes were transfected with a series of MiR21 mimic, MiR21 inhibitor, Gdf5 siRNA (si-GDF5), and flag-GDF5 constructs. The effects of MiR-21 and Gdf5 on the expression of OA related molecules were evaluated by immunofluorescence, alcian blue staining, WB, and RT-qPCR.

RESULTS

UAC altered the histological structure and extracellular matrix content of cartilage in the temporomandibular joint (TMJ), and KO of MiR21 alleviated this effect (p < 0.05). Upregulation of MiR21 influenced the expression of TMJ-OA related molecules in mandibular condylar chondrocytes via targeting Gdf5 (p < 0.05). Gdf5 overexpression significantly decreased matrix metalloproteinase 13 (MMP13) expression (p < 0.05) and reversed the effects of MiR21 (p < 0.05).

CONCLUSION

MiR21, which acts as a critical regulator of Gdf5 in chondrocytes, regulates TMJ-OA related molecules and is involved in cartilage matrix degradation, contributing to the progression of TMJ-OA. Cite this article: Bone Joint Res 2020;9(10):689-700.

Promising targets and drugs in rheumatoid arthritis: a module-based and cumulatively scoring approach

AIMS

Rheumatoid arthritis (RA) is a systematic autoimmune disorder, characterized by synovial inflammation, bone and cartilage destruction, and disease involvement in multiple organs. Although numerous drugs are employed in RA treatment, some respond little and suffer from severe side effects. This study aimed to screen the candidate therapeutic targets and promising drugs in a novel method.

METHODS

We developed a module-based and cumulatively scoring approach that is a deeper-layer application of weighted gene co-expression network (WGCNA) and connectivity map (CMap) based on the high-throughput datasets.

RESULTS

Four noteworthy RA-related modules were identified, revealing the immune- and infection-related biological processes and pathways involved in RA. HLA-DMA, HLA-DMB, HLA-DPA1, HLA-DPB1, HLA-DQB1, HLA-DRA, HLA-DRB1, BLNK, BTK, CD3D, CD4, IL2RG, INPP5D, LCK, PTPRC, RAC2, SYK, and VAV1 were recognized as the key hub genes with high connectivity in gene regulation networks and gene pathway networks. Moreover, the long noncoding RNAs (lncRNAs) in the RA-related modules, such as FAM30A and NEAT1, were identified as the indispensable interactors with the hub genes. Finally, candidate drugs were screened by developing a cumulatively scoring approach based on the selected modules. Niclosamide and the other compounds of T-type calcium channel blocker, IKK inhibitor, and PKC activator, HIF activator, and proteasome inhibitor, which harbour the similar gene signature with niclosamide, were promising drugs with high specificity and broad coverage for the RA-related modules.

CONCLUSION

This study provides not only the promising targets and drugs for RA but also a novel methodological insight into the target and drug screening.Cite this article: Bone Joint Res 2020;9(8):501-514.

Cell cycle-related lncRNAs and mRNAs in osteoarthritis chondrocytes in a Northwest Chinese Han Population

BACKGROUND

A group of differentially expressed long non-coding RNAs (lncRNAs) have been shown to play key roles in osteoarthritis (OA), although they represented only a small proportion of lncRNAs that may be biologically and physiologically relevant. Since our knowledge of regulatory functions of non-coding RNAs is still limited, it is important to gain better understanding of their relation to the pathogenesis of OA.

METHODS

We performed mRNA and lncRNA microarray analysis to detect differentially expressed RNAs in chondrocytes from three OA patients compared with four healthy controls. Then, enrichment analysis of the differentially expressed mRNAs was carried out to define disease molecular networks, pathways and gene ontology (GO) function. Furthermore, target gene prediction based on the co-expression network was performed to reveal the potential relationships between lncRNAs and mRNAs, contributing an exploration of a role of lncRNAs in OA mechanism. Quantitative RT-PCR analyses were used to demonstrate the reliability of the experimental results.

FINDINGS

Altogether 990 lncRNAs (666 up-regulated and 324 down-regulated) and 546 mRNAs (419 up-regulated and 127 down-regulated) were differentially expressed in OA samples compared with the normal ones. The enrichment analysis revealed a set of genes involved in cell cycle. In total, 854 pairs of mRNA and lncRNA were highly linked, and further target prediction appointed 12 genes specifically for their corresponding lncRNAs. The lncRNAs lncRNA-CTD-2184D3.4, ENST00000564198.1, and ENST00000520562.1 were predicted to regulate SPC24, GALM, and ZNF345 mRNA expressions in OA.

INTERPRETATION

This study uncovered several novel genes potentially important in pathogenesis of OA, and forecast the potential function of lnc-CTD-2184D3.4, especially for the cell cycle in the chondrocytes. These findings may promote additional aspects in studies of OA.

An integrated analysis of the competing endogenous RNA network and co-expression network revealed seven hub long non-coding RNAs in osteoarthritis

Aims

This study aimed to uncover the hub long non-coding RNAs (lncRNAs) differentially expressed in osteoarthritis (OA) cartilage using an integrated analysis of the competing endogenous RNA (ceRNA) network and co-expression network.

Methods

Expression profiles data of ten OA and ten normal tissues of human knee cartilage were obtained from the Gene Expression Omnibus (GEO) database (GSE114007). The differentially expressed messenger RNAs (DEmRNAs) and lncRNAs (DElncRNAs) were identified using the edgeR package. We integrated human microRNA (miRNA)-lncRNA/mRNA interactions with DElncRNA/DEmRNA expression profiles to construct a ceRNA network. Likewise, lncRNA and mRNA expression profiles were used to build a co-expression network with the WGCNA package. Potential hub lncRNAs were identified based on an integrated analysis of the ceRNA network and co-expression network. StarBase and Multi Experiment Matrix databases were used to verify the lncRNAs.

Results

We detected 1,212 DEmRNAs and 49 DElncRNAs in OA and normal knee cartilage. A total of 75 dysregulated lncRNA-miRNA interactions and 711 dysregulated miRNA-mRNA interactions were obtained in the ceRNA network, including ten DElncRNAs, 69 miRNAs, and 72 DEmRNAs. Similarly, 1,330 dysregulated lncRNA-mRNA interactions were used to construct the co-expression network, which included ten lncRNAs and 407 mRNAs. We finally identified seven hub lncRNAs, named MIR210HG, HCP5, LINC00313, LINC00654, LINC00839, TBC1D3P1-DHX40P1, and ISM1-AS1. Subsequent enrichment analysis elucidated that these lncRNAs regulated extracellular matrix organization and enriched in osteoclast differentiation, the FoxO signalling pathway, and the tumour necrosis factor (TNF) signalling pathway in the development of OA.

Conclusion

The integrated analysis of the ceRNA network and co-expression network identified seven hub lncRNAs associated with OA. These lncRNAs may regulate extracellular matrix changes and chondrocyte homeostasis in OA progress.

Cite this article:Bone Joint Res. 2020;9(3):90–98.