Global transcriptome analysis to identify critical genes involved in the pathology of osteoarthritis

Fosl2 Deficiency Predisposes Mice to Osteopetrosis, Leading to Bone Marrow Failure

Arthritis causes Fos-like 2 (Fosl2) inactivation, and various immune cells contribute to its pathogenesis. However, little is known about the role of Fosl2 in hematopoiesis and the possible pathological role of Fosl2 inactivation in the hematopoietic system in arthritis. In this study, we show that Fosl2 maintains hematopoietic stem cell (HSC) quiescence and differentiation while controlling the inflammatory response via macrophages. Fosl2-specific deletion in the hematopoietic system caused the expansion of HSCs and myeloid cell growth while affecting erythroid and B cell differentiation. Fosl2 inactivation enhanced macrophage M1 polarization and stimulated proinflammatory cytokines and myeloid growth factors, skewing HSCs toward myeloid cell differentiation, similar to hematopoietic alterations in arthritic mice. Loss of Fosl2 mediated by Vav-iCre also displays an unexpected deletion in embryonic erythro-myeloid progenitor-derived osteoclasts, leading to osteopetrosis and anemia. The reduced bone marrow cellularity in Vav-iCreFosl2f/f mice is a consequence of the reduced bone marrow space in osteopetrotic mice rather than a direct role of Fosl2 in hematopoiesis. Thus, Fosl2 is indispensable for erythro-myeloid progenitor-derived osteoclasts to maintain the medullary cavity to ensure normal hematopoiesis. These findings improve our understanding of the pathogenesis of bone-destructive diseases and provide important implications for developing therapeutic approaches for these diseases.

Identification of the shared genes in type 2 diabetes mellitus and osteoarthritis and the role of quercetin

This study investigated the underlying comorbidity mechanism between type 2 diabetes mellitus (T2DM) and osteoarthritis (OA), while also assessing the therapeutic potential of quercetin for early intervention and treatment of these two diseases. The shared genes were obtained through GEO2R, limma and weighted gene co-expression network analysis (WGCNA), and validated using clinical databases and the area under the curves (ROC). Functional enrichment analysis was conducted to elucidate the underlying mechanisms of comorbidity between T2DM and OA. The infiltration of immune cells was analysed using the CIBERSORT algorithm in conjunction with ESTIMATE algorithm. Subsequently, transcriptional regulation analysis, potential chemical prediction, gene-disease association, relationships between the shared genes and ferroptosis as well as immunity-related genes were investigated along with molecular docking. We identified the 12 shared genes (EPHA3, RASIP1, PENK, LRRC17, CEBPB, EFEMP2, UBAP1, PPP1R15A, SPEN, MAFF, GADD45B and KLF4) across the four datasets. Our predictions suggested that targeting these shared genes could potentially serve as therapeutic interventions for both T2DM and OA. Specifically, they are involved in key signalling pathways such as p53, IL-17, NF-kB and MAPK signalling pathways. Furthermore, the regulation of ferroptosis and immunity appears to be interconnected in both diseases. Notably, in this context quercetin emerges as a promising drug candidate for treating T2DM and OA by specifically targeting the shared genes. We conducted a bioinformatics analysis to identify potential therapeutic targets, mechanisms and drugs for T2DM and OA, thereby offering novel insights into molecular therapy for these two diseases.

RNA-binding proteins that are highly expressed and enriched in healthy cartilage but suppressed in osteoarthritis

Objectives: RNA-binding proteins (RBPs) have diverse and essential biological functions, but their role in cartilage health and disease is largely unknown. The objectives of this study were (i) map the global landscape of RBPs expressed and enriched in healthy cartilage and dysregulated in osteoarthritis (OA); (ii) prioritize RBPs for their potential role in cartilage and in OA pathogenesis and as therapeutic targets. Methods: Our published bulk RNA-sequencing (RNA-seq) data of healthy and OA human cartilage, and a census of 1,542 RBPs were utilized to identify RBPs that are expressed in healthy cartilage and differentially expressed (DE) in OA. Next, our comparison of healthy cartilage RNA-seq data to 37 transcriptomes in the Genotype-Tissue Expression (GTEx) database was used to determine RBPs that are enriched in cartilage. Finally, expression of RBPs was analyzed in our single cell RNA-sequencing (scRNA-seq) data from healthy and OA human cartilage. Results: Expression of RBPs was higher than nonRBPs in healthy cartilage. In OA cartilage, 188 RBPs were differentially expressed, with a greater proportion downregulated. Ribosome biogenesis was enriched in the upregulated RBPs, while splicing and transport were enriched in the downregulated. To further prioritize RBPs, we selected the top 10% expressed RBPs in healthy cartilage and those that were cartilage-enriched according to GTEx. Intersecting these criteria, we identified Tetrachlorodibenzodioxin (TCDD) Inducible Poly (ADP-Ribose) Polymerase (TIPARP) as a candidate RBP. TIPARP was downregulated in OA. scRNA-seq data revealed TIPARP was most significantly downregulated in the "pathogenic cluster". Conclusion: Our global analyses reveal expression patterns of RBPs in healthy and OA cartilage. We also identified TIPARP and other RBPs as novel mediators in OA pathogenesis and as potential therapeutic targets.

PGC-1α in osteoarthritic chondrocytes: From mechanism to target of action

Osteoarthritis (OA) is one of the most common degenerative joint diseases, often involving the entire joint. The degeneration of articular cartilage is an important feature of OA, and there is growing evidence that the mitochondrial biogenesis master regulator peroxisome proliferator-activated receptor γ coactivator 1α (PGC-1α) exert a chondroprotective effect. PGC-1α delays the development and progression of OA by affecting mitochondrial biogenesis, oxidative stress, mitophagy and mitochondrial DNA (mtDNA) replication in chondrocytes. In addition, PGC-1α can regulate the metabolic abnormalities of OA chondrocytes and inhibit chondrocyte apoptosis. In this paper, we review the regulatory mechanisms of PGC-1α and its effects on OA chondrocytes, and introduce potential drugs and novel nanohybrid for the treatment of OA which act by affecting the activity of PGC-1α. This information will help to further elucidate the pathogenesis of OA and provide new ideas for the development of therapeutic strategies for OA.

Screening crucial lncRNAs and genes in osteoarthritis by integrated analysis

BACKGROUND

Osteoarthritis (OA) is one of the most frequent chronic diseases with high morbidity worldwide, marked by degradation of the cartilage and bone, joint instability, stiffness, joint space stenosis and subchondral sclerosis. Due to the elusive mechanism of osteoarthritis (OA), we aimed to identify potential markers for OA and explore the molecular mechanisms underlying OA.

METHODS

Expression profiles data of OA were collected from the Gene Expression Omnibus database to identify differentially expressed mRNAs (DEmRNAs) and differentially expressed lncRNAs (DElncRNAs) in OA. Functional annotation and protein-protein interaction (PPI) networks were performed. Then, nearby DEmRNAs of DElncRNAs was obtained. Moreover, GO and KEGG pathway enrichment analysis of nearby DEmRNAs of DElncRNAs was performed. Finally, expression validation of selected mRNAs and lncRNAs was performed by quantitative reverse transcriptase-polymerase chain reaction.

RESULTS

In total, 2080 DEmRNAs and 664 DElncRNAs were determined in OA. PI3K-Akt signaling pathway, Endocytosis and Rap1 signaling pathway were significantly enriched KEGG pathways in OA. YWHAB, HSPA8, NEDD4L and SH3KBP1 were four hub proteins in PPI network. The AC093484.4/TRPV2 interact pair may be involved in the occurrence and development of OA.

CONCLUSION

Our study identified several DEmRNAs and DElncRNAs associated with OA. The molecular characters could provide more information for further study on OA.

RNA-Seq Reveals the mRNAs, miRNAs, and lncRNAs Expression Profile of Knee Joint Synovial Tissue in Osteoarthritis Patients

Osteoarthritis (OA) is a chronic disease common in the elderly population and imposes significant health and economic burden. Total joint replacement is the only currently available treatment but does not prevent cartilage degeneration. The molecular mechanism of OA, especially the role of inflammation in disease progression, is incompletely understood. We collected knee joint synovial tissue samples of eight OA patients and two patients with popliteal cysts (controls), measured the expression levels of lncRNAs, miRNAs, and mRNAs in these tissues by RNA-seq, and identified differentially expressed genes (DEGs) and key pathways. In the OA group, 343 mRNAs, 270 lncRNAs, and 247 miRNAs were significantly upregulated, and 232 mRNAs, 109 lncRNAs, and 157 miRNAs were significantly downregulated. mRNAs potentially targeted by lncRNAs were predicted. Nineteen overlapped miRNAs were screened based on our sample data and GSE 143514 data. Pathway enrichment and functional annotation analyses showed that the inflammation-related transcripts CHST11, ALDH1A2, TREM1, IL-1β, IL-8, CCL5, LIF, miR-146a-5p, miR-335-5p, lncRNA GAS5, LINC02288, and LOC101928134 were differentially expressed. In this study, inflammation-related DEGs and non-coding RNAs were identified in synovial samples, suggesting that competing endogenous RNAs have a role in OA. TREM1, LIF, miR146-5a, and GAS5 were identified to be OA-related genes and potential regulatory pathways. This research helps elucidate the pathogenesis of OA and identify novel therapeutic targets for this disorder.

Multiple Genetic Loci Associated with Pug Dog Thoracolumbar Myelopathy

Pug dogs with thoracolumbar myelopathy (PDM) present with a specific clinical phenotype that includes progressive pelvic limb ataxia and paresis, commonly accompanied by incontinence. Vertebral column malformations and lesions, excessive scar tissue of the meninges, and central nervous system inflammation have been described. PDM has a late onset and affects more male than female dogs. The breed-specific presentation of the disorder suggests that genetic risk factors are involved in the disease development. To perform a genome-wide search for PDM-associated loci, we applied a Bayesian model adapted for mapping complex traits (BayesR) and a cross-population extended haplotype homozygosity test (XP-EHH) in 51 affected and 38 control pugs. Nineteen associated loci (harboring 67 genes in total, including 34 potential candidate genes) and three candidate regions under selection (with four genes within or next to the signal) were identified. The multiple candidate genes identified have implicated functions in bone homeostasis, fibrotic scar tissue, inflammatory responses, or the formation, regulation, and differentiation of cartilage, suggesting the potential relevance of these processes to the pathogenesis of PDM.

Integrated Serum Metabolomics and Network Pharmacology to Reveal the Interventional Effects of Quzhi Decoction against Osteoarthritis Pain

Objectives

Chronic pain, the main symptom of knee osteoarthritis (OA), remains the primary reason for decreased functional capacity. Quzhi decoction, a TCM prescription, is effective in treating chronic pain in OA, but the potential mechanisms require further exploration.

Methods

An anterior cruciate ligament transection (ACLT) rat model was established, and pain-like behavior was evaluated. Metabolomics analysis of serum samples was performed to identify differential metabolites, and network pharmacology was used to identify potential targets of Quzhi decoction for the treatment of OA. Finally, we constructed a comprehensive network of serum metabolomics and network pharmacology. At the same time, the obtained key targets were verified by molecular docking.

Results

Quzhi decoction was shown to attenuate pain-like behavior and joint inflammation in OA rats. Through serum metabolomics, thirty potentially significant metabolites were found to be involved in the therapeutic effects of Quzhi decoction against OA pain. According to network pharmacology, 107 active drug components were matched with 115 disease targets, which was partly consistent with the metabolomics findings. Further analysis focused on 6 key targets, including CYP3A4, PLA2G4A, PTGS1, PTGS2, TYR, and ALOX5, and their associated core metabolites and pathways. Molecular docking results showed that the related targets had high affinity with the active pharmaceutical ingredients in Quzhi decoction.

Conclusion

The effect of Quzhi decoction on OA pain may be related to the inhibition of joint inflammation, mainly through disturbing arachidonic acid metabolism, tyrosine metabolism, and leukotriene metabolism. Further systematic molecular biology experiments are needed to verify the accurate mechanism.

Cartilage debris and osteoarthritis risk factors influence gene expression in the synovium in end stage osteoarthritis

BACKGROUND

Gene expression in healthy synovium remains poorly characterised. Thus, synovial functional activity changes associated with osteoarthritis (OA) are difficult to define. This study sought to identify differentially expressed genes (DEG) of end-stage OA and assess the influence of OA risk factors on these DEG.

METHODS

Anonymised patient clinical data and x-ray images were analysed. Osteoarthritic and non-osteoarthritic patients with soft tissue or traumatic knee injuries were matched for body mass index (BMI) and sex. Tissue samples were partitioned for immunocytochemistry (IHC) and microarray analysis. Multiple bioinformatics applications were utilised to determine changes in functional and canonical pathway activation.

RESULTS

Age, disease-modifying injections and hypertension were confounding factors between patient groups. Inflammation was present in all tissues. Cartilage debris and inflammatory aggregates were noted in many osteoarthritic patient tissues. IHC and expression analyses revealed upregulation of synoviolin 1 (SYVN1) in osteoarthritic synovium. Significant differential expression was noted in 2084 genes. Osteoarthritic synovium displayed a significant upregulation of 95% of DEG coding for proteins, relative to non-osteoarthritic synovium tissues. Unfolded protein response (UPR)-related genes were upregulated in osteoarthritic synovium; gene expression of molecules within many canonical pathways including protein ubiquitination and UPR pathways was modified by BMI and sex.

CONCLUSIONS

The synovium of all three pathologies exhibited elements of an inflammatory response. Cartilage debris, age, BMI and sex influence DEG of osteoarthritic synovium. UPR pathway is the top deregulated canonical pathway identified in osteoarthritic synovium regardless of BMI and sex, while typical OA-associated inflammatory and matrix gene responses were minimal.

Exploring the Mystery of Osteoarthritis using Bioinformatics Analysis of Cartilage Tissue

BACKGROUND

Osteoarthritis (OA) is a kind of chronic disease relating to joints, which seriously affectsthe daily life activities of the elderly and can also lead to disability. However, the pathogenesis of OA is still unclear, which leads to limited treatment and the therapeutic effect far from people's expectations. This study aims to filter out key genes in the pathogenesis of OA and explore their potential role in the occurrence and development of OA.

METHODS

The dataset of GSE117999 was obtained and analyzed in order to identify the differentially expressed genes (DEGs), hub genes and key genes. We also identified potential miRNAs which may play a major role in the pathogenesis of OA, and verified their difference in OA by real-time quantitative PCR (RT-qPCR). DGldb was found to serve as an indicator to identify drugs with potential therapeutic effects on key genes and Receiver Operating Characteristic (ROC) analysis was used for identifying underlying biomarkers of OA.

RESULTS

We identified ten key genes, including MDM2, RB1, EGFR, ESR1, UBE2E3, WWP1, BCL2, OAS2, TYMS and MSH2. Then, we identified hsa-mir-3613-3p, hsa-mir-548e-5p and hsamir- 5692a to be potentially related to key genes. In addition, RT-qPCR confirmed the differential expression of identified genes in mouse cartilage with or without OA. We then identified Etoposide and Everolimus, which were potentially specific to the most key genes. Finally, we speculated that ESR1 might be a potential biomarker of OA.

CONCLUSION

In this study, potential key genes related to OA and their biological functions were identified, and their potential application value in the diagnosis and treatment of OA has been demonstrated, which will help us to improve the therapeutic effect of OA.

软骨细胞线粒体生物发生在骨关节炎发病机制中的作用

OBJECTIVE

To summarize the role of chondrocytes mitochondrial biogenesis in the pathogenesis of osteoarthritis (OA), and analyze the applications in the treatment of OA.

METHODS

A review of recent literature was conducted to summarize the changes in mitochondrial biogenesis in the course of OA, the role of major signaling molecules in OA chondrocytes, and the prospects for OA therapeutic applications.

RESULTS

Recent studies reveales that mitochondria are significant energy metabolic centers in chondrocytes and its dysfunction has been considered as an essential mechanism in the pathogenesis of OA. Mitochondrial biogenesis is one of the key processes maintaining the normal quantity and function of mitochondria, and peroxisome proliferator-activated receptor-gamma coactivator 1 alpha (PGC-1α) is the central regulator of this process. A regulatory network of mitochondrial biogenesis with PGC-1α as the center, adenosine monophosphate-activated protein kinase, sirtuin1/3, and cyclic adenosine monophosphate response element-binding protein as the main upstream regulatory molecules, and nuclear respiratory factor 1, estrogen-related receptor α, and nuclear respiratory factor 2 as the main downstream regulatory molecules has been reported. However, the role of mitochondrial biogenesis in OA chondrocytes still needs further validation and in-depth exploration. It has been demonstrated that substances such as puerarin and omentin-1 can retard the development of OA by activating the damaged mitochondrial biogenesis in OA chondrocytes, which proves the potential to be used in the treatment OA.

CONCLUSION

Mitochondrial biogenesis in chondrocytes plays an important role in the pathogenesis of OA, and further exploring the related mechanisms is of great clinical significance.

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.

Comprehensive analysis of lncRNA and mRNA based on expression microarray profiling reveals different characteristics of osteoarthritis between Tibetan and Han patients

Background

Osteoarthritis (OA) is thought to be the most prevalent chronic joint disease, especially in Tibet of China. Here, we aimed to explore the integrative lncRNA and mRNA landscape between the OA patients of Tibet and Han.

Methods

The lncRNA and mRNA expression microarray profiling was performed by SurePrint G3 Human Gene Expression 8x60K v2 Microarray in articular cartilage samples from OA patients of Han nationality and Tibetans, followed by GO, KEGG, and trans-regulation and cis-regulation analysis of lncRNA and mRNA.

Results

We found a total of 117 lncRNAs and 297 mRNAs differently expressed in the cartilage tissues of Tibetans (n = 5) comparing with those of Chinese Han (n = 3), in which 49 lncRNAs and 158 mRNAs were upregulated, and 68 lncRNAs and 139 mRNAs were downregulated. GO and KEGG analysis showed that several unreported biological processes and signaling pathways were particularly identified. LncRNA-mRNA co-expression analysis revealed a remarkable lncRNA-mRNA relationship, in which OTOA may play a critical role in the different mechanisms of the OA progression between Tibetans and Chinese Han.

Conclusion

This study identified different lncRNA/mRNA expression profiling between OA patients of Tibetans and Han, which were involved in many characteristic biological processes and signaling pathways.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13018-021-02213-y.

Pro-osteogenic Effects of WNT in a Mouse Model of Bone Formation Around Femoral Implants

Wnt signaling maintains homeostasis in the bone marrow cavity: if Wnt signaling is inhibited then bone volume and density would decline. In this study, we identified a population of Wnt-responsive cells as osteoprogenitor in the intact trabecular bone region, which were responsible for bone development and turnover. If an implant was placed into the long bone, this Wnt-responsive population and their progeny contributed to osseointegration. We employed Axin2Cre^CreERT2/+;R26^mTmG/+ transgenic mouse strain in which Axin2-positive, Wnt-responsive cells, and their progeny are permanently labeled by GFP upon exposure to tamoxifen. Each mouse received femoral implants placed into a site prepared solely by drilling, and a single-dose liposomal WNT3A protein was used in the treatment group. A lineage tracing strategy design allowed us to identify cells actively expressing Axin2 in response to Wnt signaling pathway. These tools demonstrated that Wnt-responsive cells and their progeny comprise a quiescent population residing in the trabecular region. In response to an implant placed, this population becomes mitotically active: cells migrated into the peri-implant region, up-regulated the expression of osteogenic proteins. Ultimately, those cells gave rise to osteoblasts that produced significantly more new bone in the peri-implant region. Wnt-responsive cells directly contributed to implant osseointegration. Using a liposomal WNT3A protein therapeutic, we showed that a single application at the time of implant placed was sufficient to accelerate osseointegration. The Wnt-responsive cell population in trabecular bone, activated by injury, ultimately contributes to implant osseointegration. Liposomal WNT3A protein therapeutic accelerates implant osseointegration in the long bone.

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.

Ubiquitin-specific protease 49 attenuates IL-1β-induced rat primary chondrocyte apoptosis by facilitating Axin deubiquitination and subsequent Wnt/β-catenin signaling cascade inhibition

Osteoarthritis (OA) is an age-related chronic joint degenerative disease. Interleukin 1 beta (IL-1β) is considered a marker for the progression of OA. In this study, we found that Ubiquitin-Specific Peptidase 49 (USP49) was significantly less expressed in OA patients compared with healthy individuals. Treating primary rat chondrocytes with different concentrations of IL-1β resulted in decreased Usp49 expression, while Usp49 overexpression could attenuate IL-1β-induced chondrocyte apoptosis by promoting Axin deubiquitination. The deubiquitination of Axin led to the accumulation of the protein, which in turn resulted in β-catenin degradation and Wnt/β-catenin signaling cascade inhibition. Interestingly, we also found that [6]-gingerol, an anti-OA drug, could upregulate the protein level of Usp49 and suppress the Wnt/β-catenin signaling cascade in primary rat chondrocytes. Taken together, our study not only demonstrates that Usp49 can negatively regulate the progression of OA by inhibiting the Wnt/β-catenin signaling cascade, but also elucidates the underlying molecular mechanisms.

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.

Integration of transcriptome-wide association study and messenger RNA expression profile to identify genes associated with osteoarthritis

Aims

Osteoarthritis (OA) is the most prevalent joint disease. However, the specific and definitive genetic mechanisms of OA are still unclear.

Methods

Tissue-related transcriptome-wide association studies (TWAS) of hip OA and knee OA were performed utilizing the genome-wide association study (GWAS) data of hip OA and knee OA (including 2,396 hospital-diagnosed hip OA patients versus 9,593 controls, and 4,462 hospital-diagnosed knee OA patients versus 17,885 controls) and gene expression reference to skeletal muscle and blood. The OA-associated genes identified by TWAS were further compared with the differentially expressed genes detected by the messenger RNA (mRNA) expression profiles of hip OA and knee OA. Functional enrichment and annotation analysis of identified genes was performed by the DAVID and FUMAGWAS tools.

Results

We detected 33 common genes, eight common gene ontology (GO) terms, and one common pathway for hip OA, such as calcium and integrin-binding protein 1 (CIB1) (P_TWAS = 0.025, FC_mRNA = -1.575 for skeletal muscle), adrenomedullin (ADM) (P_TWAS = 0.022, FC_mRNA = -4.644 for blood), Golgi apparatus (P_TWAS <0.001, P_mRNA = 0.012 for blood), and phosphatidylinositol 3' -kinase-protein kinase B (PI3K-Akt) signalling pathway (P_TWAS = 0.033, P_mRNA = 0.005 for blood). For knee OA, we detected 24 common genes, eight common GO terms, and two common pathways, such as histocompatibility complex, class II, DR beta 1 (HLA-DRB1) (P_TWAS = 0.040, FC_mRNA = 4.062 for skeletal muscle), Follistatin-like 1 (FSTL1) (P_TWAS = 0.048, FC_mRNA = 3.000 for blood), cytoplasm (P_TWAS < 0.001, P_mRNA = 0.005 for blood), and complement and coagulation cascades (P_TWAS = 0.017, P_mRNA = 0.001 for skeletal muscle).

Conclusion

We identified a group of OA-associated genes and pathways, providing novel clues for understanding the genetic mechanism of OA.

Cite this article:Bone Joint Res. 2020;9(3):130–138.

Profile of Matrix-Remodeling Proteinases in Osteoarthritis: Impact of Fibronectin

The extracellular matrix (ECM) is a complex and specialized three-dimensional macromolecular network, present in nearly all tissues, that also interacts with cell surface receptors on joint resident cells. Changes in the composition and physical properties of the ECM lead to the development of many diseases, including osteoarthritis (OA). OA is a chronic degenerative rheumatic disease characterized by a progressive loss of synovial joint function as a consequence of the degradation of articular cartilage, also associated with alterations in the synovial membrane and subchondral bone. During OA, ECM-degrading enzymes, including urokinase-type plasminogen activator (uPA), matrix metalloproteinases (MMPs), and a disintegrin and metalloproteinase with thrombospondin motifs (ADAMTSs), cleave ECM components, such as fibronectin (Fn), generating fibronectin fragments (Fn-fs) with catabolic properties. In turn, Fn-fs promote activation of these proteinases, establishing a degradative and inflammatory feedback loop. Thus, the aim of this review is to update the contribution of ECM-degrading proteinases to the physiopathology of OA as well as their modulation by Fn-fs.

Identification of Key Genes and Pathways Associated with Sex Differences in Osteoarthritis Based on Bioinformatics Analysis

Sex differences have been suggested to play critical roles in the pathophysiology of osteoarthritis (OA), resulting in sex-specific prevalence and incidence. However, their roles in the development of OA remain largely unknown. The aim of this study was to screen out key genes and pathways mediating biological differences between OA females after menopause and OA males. First, the gene expression data of GSE36700 and GSE55457 were downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) between sexes were identified using R software, respectively. The overlapping DEGs were obtained. Then, protein-protein interactive (PPI) network was constructed to further analyze interactions between the overlapping DEGs. Finally, enrichment analyses were separately performed using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes tools. In our results, a total of 278 overlapping DEGs were identified between OA females after menopause and OA males, including 219 upregulated and 59 downregulated genes. In the PPI network, seven hub genes were identified, including EGF, ERBB2, CDC42, PIK3R2, LCK, CBL, and STAT1. Functional enrichment analysis revealed that these genes were mainly enriched in PI3K-Akt signaling pathway, osteoclast differentiation, and focal adhesion. In conclusion, the results in the current study suggest that pathways of PI3K-Akt, osteoclast differentiation, and focal adhesion may play important roles in the development of OA females after menopause. EGFR, ERBB2, CDC42, and STAT1 may be key genes related to OA progression in postmenopausal women and may be promising therapeutic targets for OA.

Analysis of gene expression and methylation datasets identified ADAMTS9, FKBP5, and PFKBF3 as biomarkers for osteoarthritis

BACKGROUND

Osteoarthritis (OA) is a kind of chronic osteoarthropathy and degenerative joint disease. Epigenetic regulation in the gene expression dynamics has become increasingly important in OA. We performed a combined analysis of two types of microarray datasets (gene expression and DNA methylation) to identify methylation-based key biomarkers to provide a better understanding of molecular biological mechanisms of OA.

METHODS

We obtained two expression profiling datasets (GSE55235, GSE55457) and one DNA methylation profiling data set (GSE63695) from the Gene Expression Omnibus. First, differentially expressed genes (DEGs) between patients with OA and controls were identified using the Limma package in R(v3.4.4). Then, function enrichment analysis of DEGs was performed using a DAVID database. For DNA methylation datasets, ChAMP methylation analysis package was used to identify differential methylation genes (DMGs). Finally, a comprehensive analysis of DEGs and DMGs was conducted to identify genes that exhibited differential expression and methylation simultaneously.

RESULTS

We identified 112 DEGs and 2,896 DMGs in patients with OA compared with controls. Functional analysis of DEGs obtained that inflammatory responses, immune responses, and positive regulation of apoptosis, tumor necrosis factor (TNF) signaling pathway, and osteoclast differentiation may be involved in the pathogenesis of OA. Cross-analysis revealed 26 genes that exhibited differential expression and methylation in OA. Among them, ADAMTS9, FKBP5, and PFKBF3 were identified as valuable methylation-based biomarkers for OA.

CONCLUSION

In summary, our study identified different molecular features between patients with OA and controls. This may provide new clues for clarifying the pathogenetic mechanisms of OA.