Back to Journals » Journal of Inflammation Research » Volume 17
Network Analysis of Osteoarthritis Progression Using a Steiner Minimal Tree Algorithm
Authors Xie Y, Shao F, Ji Y, Feng D, Wang L, Huang Z, Wu S, Sun F, Jiang H, Miyamoto A, Wang H, Zhang C
Received 3 September 2023
Accepted for publication 9 December 2023
Published 18 May 2024 Volume 2024:17 Pages 3201—3209
DOI https://doi.org/10.2147/JIR.S438407
Checked for plagiarism Yes
Review by Single anonymous peer review
Peer reviewer comments 3
Editor who approved publication: Professor Ning Quan
Yujie Xie,1 Fanglin Shao,1 Yuxiu Ji,1 Dechao Feng,1,2 Ling Wang,3 Zonghai Huang,3 Shengjian Wu,1 Fuhua Sun,1 Hong Jiang,4 Akira Miyamoto,5 Haiming Wang,6 Chi Zhang1
1Rehabilitation Medicine Department, The Affiliated Hospital of Southwest Medical University, Luzhou, Sichuan, People’s Republic of China; 2Department of Urology, Institute of Urology, West China Hospital, Sichuan University, Chengdu, Sichuan, People’s Republic of China; 3School of Information and Communication Engineering, University of Electronic Science and Technology, Chengdu, Sichuan, People’s Republic of China; 4Rehabilitation Medicine Department, Xichong County People’s hospital, Nanchong, Sichuan, People’s Republic of China; 5Faculty of Rehabilitation, Nishikyushu University, Kansaitama, Japan; 6Department of Rehabilitation Medicine, the First Affiliated Hospital of Zhengzhou University, Zhengzhou, Henan, People’s Republic of China
Correspondence: Haiming Wang, Department of Rehabilitation Medicine, the First Affiliated Hospital of Zhengzhou University, Zhengzhou, Henan, People’s Republic of China, Tel +868950771, Fax +86 830 8950771, Email [email protected] Chi Zhang, Rehabilitation Medicine Department, The Affiliated Hospital of Southwest Medical University, Luzhou, Sichuan, People’s Republic of China, Email [email protected]
Purpose: To provide a comprehensive analysis of associated genes with osteoarthritis (OA). Here, we reported a network analysis of OA progression by using a Steiner minimal tree algorithm.
Methods: We collected the OA-related genes through screening the publications in MEDLINE. We performed functional analysis to analyze the associated biochemical pathways of the OA-related genes. Pathway crosstalk analysis was constructed to explore interactions of the enriched pathways. Steiner minimal tree algorithm was used to analyze molecular pathway networks. The average clustering coefficient was compared with the corresponding values of the Osteoarthritis-specific network. The new finding RNA was compared with former single-cell RNA-seq analysis results.
Results: A gene set with 177 members reported to be significantly associated with Osteoarthritis was collected from 187 studies. Functional enrichment analysis revealed a specific related-OA gene including skeletal system development, cytokine-mediated signaling pathway, inflammatory response, cartilage development, and extracellular matrix organization. We performed a pathway crosstalk analysis among the 72 significantly enriched pathways. A total of 151 of the 177 genes in the Osteoarthritis gene set were included in the human interactome network. There were 31 genes in the former single-cell RNA-seq analysis results. The CLU, ENO1, SRRM1, UBC, HMGB1, NR3C1, NOTCH2NL, and CBX5 have significantly increased expression in seven molecularly defined populations of OA cartilage.
Conclusion: The Steiner tree-based approach finds new biological molecules associated with OA genes.
Keywords: network analysis, Steiner minimal tree, osteoarthritis, bioinformatics, osteoarthritis genes
Introduction
Osteoarthritis (OA) is the leading cause of disability and source of societal cost in older adults.1 Risk factors associated with OA include joint injuries, aging, and obesity.2 It is projected that by 2040, OA will affect 78 million people.3 The impact of OA on individuals includes pain, loss of mobility, and loss of independence, with 25% of patients unable to perform normal daily activities. OA is characterized by various factors including subchondral bone remodeling,4 degeneration of the meniscus,5,6 inflammation, and fibrosis of the infrapatellar fat pad and synovial membrane.7
OA will remain a large and growing global problem for clinical and public health systems. Although OA is a significant medical condition, there are currently no registered therapies available that can stop the structural damage caused by the disease. Symptom-modifying interventions only offer moderate long-term effects, at best.8 To advance in the treatment of OA, it is crucial to have a deep understanding of the cellular and molecular pathophysiology.
Several genes and pathways associated with mechanical, inflammatory, and synovial factors have been identified in OA.9 The identification of the possible essential genes such as MMP-3, MMP-13, VEGF, TIMP-1, TGFB1, FGF18 has important influence on molecular causes.10–13
Although studies have identified important pathogenic genes in the development of OA, the underlying molecular mechanism remains unclear.14 Lots of genes has effects on OA progress rather than one or two major genes with large effects.15 Network and pathway analysis reveals crosstalk between genes. So, a comprehensive analysis of potentially causal genes can provide many important insights finding.16 We used a comprehensive analysis of genes associated with OA. Then, functional enrichment analyzes explore the significant biological themes within these genetic factors. We intended to analyze the topological characteristics of related genes in the context of a human protein-protein interaction network. To accomplish this, we first inferred the specific molecular network using the Steiner minimal tree algorithm. By isolating the network specific to Osteoarthritis from the larger human interaction group network, we gained insights into the potential pathological molecular network associated with the Osteoarthritis gene set. Additionally, we compared the significant RNA findings with previous results obtained from single-cell RNA-seq analysis.17 This comprehensive analysis provides systematic insights into the mechanisms underlying OA.
Materials and Methods
Identification of OA-Related Genes
Candidate genes associated with OA were curated by retrieving the human genetic association studies deposited in MEDLINE (http://www.ncbi.nlm.nih.gov/pubmed/). We searched OA with the term (osteoarthritis [MeSH]) and (polymorphism [MeSH] or genotype [MeSH] or alleles [MeSH]) not (neoplasms [MeSH]). Non-English articles studies were excluded. Two authors screened the titles and abstracts of the articles, and the full text of all eligible studies was reviewed for inclusion and exclusion criteria. Any discrepancies in literature screening were resolved by the third author. By 2020–12-01, we found 1203 publications for OA. Then, we collected the association studies after reviewing the abstracts, and narrowed the selection publication by focusing on OA. We only included publication that results support the conclusions. The negative or insignificant associations reports were excluded. The genes significantly associated with OA in these studies were selected for this study.
Functional Enrichment Analysis of OA-Related Genes
The functional features of the OA-related genes were examined by ClusterProfiler based on Gene Ontology (GO), BioCarta and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. GO terms of biological processes with a p-adjust value smaller than 0.01 were kept as the significantly enriched ones. All the pathways with one or more genes overlapping the candidate genes were extracted. A p-value was assigned to each pathway by Fisher’s test to indicate the significant overlap between the pathway and the input gene. The pathways were significantly enriched with a p-adjust < 0.01.
Pathway Crosstalk Analysis
We further performed pathway crosstalk following the former way.18 To describe the overlap between any given pair of pathways, two measurements were computed. The Jaccard Coefficient (JC) \ (JC=|\frac {A\cap B {A\cup B} |\) and the Overlap Coefficient (OC) \ (OC=\frac {|A\cap B|} {min (|A|, |B|)} \). Where A and B are the lists of genes included in the two tested pathways. To construct the pathway crosstalk, we implemented the following procedure:
Select a set of pathways for crosstalk analysis. Only the pathways with a p-adjust < 0.01 were used. Meanwhile, the pathways containing less than six candidate genes were removed.
Count the number of shared candidate genes between any pair of pathways. Pathway pairs with less than seven overlapped genes were removed.
Calculate the overlap of all pathway pairs. All the pathway pairs were ranked according to their JC and OC values. The overlapping level between the two pathways was measured by the JC and OC average scores. The pathways could be grouped into one major module by their crosstalk analysis results.
Visualize the selected pathway crosstalk with the packages cytoscape.18
Construction of the Human Interactome
We constructed a comprehensive and reliable human interactome to investigate the correlation between the genes. First, we downloaded the human protein-protein interaction (PPI) data from the Protein Interaction Network Analysis (PINA) platform (release version: 21 May 2014). The database pooled and curated non-redundant physical interaction data from six databases: IntAct, BioGRID, DIP, HPRD, MINT, and MIPS/MPact. After excluding the redundant and self-interacting pairs and using org.Hs.ens.db R packages to map these interactome data onto gene SYMBOL, we constructed a comprehensive human physical interactome by merging the two data sets, which contained 17,457 nodes and 233,333 edges.
Construction of OA-Specific Network
We intended to explore the disease subnetwork extraction. We used the Steiner Net to extract specific potential pathological networks and take the collected OA genes as input seeds. Then, we generated 1000 random networks with the same number of interactions as the OA-specific network using the Erdos-Renyi model in the R igraph package to assess the non-randomness. Next, we calculated the average values of the shortest-path distance and clustering coefficient. The significance level of non-randomness needs two requirements; 1) the number of random networks with average shortest-path distance (ND) was smaller than the OA-specific networks, 2) the number of random networks with an average clustering coefficient (NC) was higher than the observed clustering coefficient. Finally, the empirical p-value was calculated by using ND/1000 and NC/1000.
Single-Cell RNA-Seq Analysis Results Identify
The important RNA was compared with former single-cell RNA-seq analysis results.17 Quanbo Ji has identified seven molecularly defined populations of chondrocytes in the human OA cartilage, including ECs, effector chondrocytes (RegCs), regulatory chondrocytes (ProCs), proliferative chondrocytes (preHTCs), pre-hypertrophic chondrocytes (FCs), fibrocartilage chondrocytes (HTCs), hypertrophic chondrocytes (HomCs) and homeostatic chondrocytes. They also presented gene expression profiles at different OA stages at single-cell resolution. We identified our crosstalk results from these expression profiles.
Results
Identification of Genes Associated with OA
(Supplementary Table 1) provided all genes related to OA. Altogether, a gene set with 177 members significantly associated with OA was collected from 187 studies.
Significant Enrichment Pathways for OA
Biological Functions and Pathway Enrichment Analysis Enriched in Gene Set Some GO terms significantly enriched in the candidate genes (Supplementary Table 2), including those associated with skeletal system development, cytokine-mediated signaling pathway, inflammatory response, cartilage development, extracellular matrix organization, negative regulation of cell population proliferation, ossification, positive regulation of tyrosine phosphorylation of STAT protein and positive regulation of epithelial to mesenchymal transition. We found 72 significant enrichment pathways for OA (Supplementary Table 3) through.
Crosstalk Among Significantly Enriched Pathways
To understand how these 72 pathways interact with each other, we performed a pathway crosstalk analysis. A total of 66 pathways were containing six or more members in the OA gene set. A number of 44 pathways met the criterion for crosstalk analysis. Each pathway shared at least seven genes with one or more other pathways (Figure 1).
Gene Set Network Topological Characteristics
PPI analyzed the topological properties of nodes and interactions. The degree of the constructed network was analyzed in the context of the constructed human interactome. We intended to analyze the two other gene sets network features including the osteoporosis list with 323 members and the cancer genes with 506 genes (Supplementary Table 4). The nodes with specific degrees were scattered in a range from 1 to 836. The nodes were measured by the number of genes connecting with a given gene. The 151 of 177 genes mapped onto the human interactome network with the 52.51 mean degree, which each gene connected with 52.51 other genes on average. The 300 of 323 genes could be mapped onto the human interactome in osteoporosis, with 54.05 average degree. The 485 out of 506 had the corresponding nodes in the cancer genes with the 74.97 average degree. Further, 52.98% (80/151) and 51% (153/300) genes fell in the degree interval of 1–20, respectively for OA and Osteoporosis, while the cancer genes only 30.31%.
OA-Specific Molecular Network Inference
Steiner minimal tree algorithm linked the maximal and minimal members of genes nodes. The OA network contained 182 nodes and 181 edges (Figure 2). For these random subnetworks, the mean shortest-path distance was 6.51, significantly larger than that of the OA-related network (shortest-path distance, 4.06; empirical p=0). The average clustering coefficient of the random networks was 0.01, statistically significantly less than the OA distinctive network (clustering coefficient, 0; empirical p=0.72). As specified, 151 of the 177 genes in the OA gene set were included in the extracted OA-specific network. A total of 85.31% of the genes in the OA gene set and 82.96% of 182 genes in the OA-specific network, indicating a high coverage of the OA gene set in the subnetwork. The remaining 31 genes in the OA-specific network were not excluded, as they were non-input genes screened from the Steiner minimum network. These 31 genes are not part of the Osteoarthritis gene concentration, but they are involved in the biological process and have close interactions with each other. These genes provide potential candidate genes for further exploration in the OA biological processes. In addition, the expression of CLU, ENO1, SRRM1, UBC, HMGB1, NR3C1, NOTCH2NL, and CBX5 were significantly increased in seven molecularly defined populations of OA cartilage, which is consistent with the results of former single-cell RNA-seq analysis, making the consistency more convincing. Figure 3 presents violin plots depicting the expression levels of specific marker genes across seven different populations in human OA cartilage. These populations include effector chondrocytes, regulatory chondrocytes, proliferative chondrocytes, prehypertrophic chondrocytes, fibrocartilage chondrocytes, hypertrophic chondrocytes, and homeostatic chondrocytes.
Discussion
OA is characterized by cartilage degeneration, synovial and infrapatellar fat pad inflammation, and subchondral bone remodeling.19,20 In this study, we aimed to conduct a comprehensive collection of pathogenic associated with OA by including three databases. We found the intersection of all the signaling pathways through the GO and the KEGG functional enrichment analyses. Then, we further analyzed these novel genes in the context of the human OA protein-protein interaction network. Next, the Steiner minimal tree algorithm was used in OA-specific molecular network analysis.
Former research17 has identified seven molecularly defined populations of chondrocytes in the human OA cartilage, including three novel phenotypes with distinct functions. It presented gene expression profiles at different OA stages at single-cell resolution. These genes had close interaction with genes known to be related to the biological processes involved in OA.21 Thus, these genes provided a list of potential candidates for further exploration. After reviewing all these genes, we finally identified 31 important genes. Notably, CLU, ENO1, SRRM1, UBC, HMGB1, NR3C1, NOTCH2NL, and CBX5 showed significantly increased expression in OA chondrocytes.
Clusterin (CLU) has been found as a ubiquitous glycoprotein and its upregulation has been found in OA cartilage.22 Also known as apolipoprotein J, CLU acts as an ATP-independent holdase chaperone, preventing protein aggregation and proteotoxicity.23 It is a disulfide-linked heterodimeric protein of approximately 60 kDa, involved in the clearance of cellular debris and the regulation of apoptosis. Recent studies suggest that the intracellular form of CLU may suppress stress-induced apoptosis and the secreted form of CLU functions as an extracellular chaperone that prevents protein aggregation.24 Increased levels of CLU mRNA and protein in the joint environments of knee OA may reflect the severity of the condition. Csaba Matta’s findings demonstrate that the secretion of clusterin can be regulated by interleukin-1β and tumor necrosis factor-α, influencing cartilage degradation.25 Furthermore, three neopeptides associated with Clusterin show potential as markers for OA disease.26
Enolase 1(ENO1) has been implicated in autoimmune and inflammatory diseases.27 ENO1 has been identified as a novel stimulatory receptor on monocytes.27 ENO1 is a multifunctional glycolytic enzyme expressed in the cell-surface and has been found to be increased rapidly in response to inflammatory loops in rheumatoid arthritis.28 Weijuan Ma investigated the protein changes in Kashin-Beck disease cartilage and founded the candidate proteins α-enolase.29 Additionally, proteomic analysis of human osteoarthritic chondrocytes reveals protein changes, and Western blotting and immunohistochemistry in tissue cartilage were used to identify a decrease in ENO1.30
Evidence supports increased expression of high mobility group box 1 (HMGB1) in OA cartilage.31 HMGB1 protein is associated with chemoattractant properties.32 Chondrocytes release HMGB1, which has migratory effects on chondrogenic progenitor cells.33 Additionally, increased expression of HMGB1 has been observed in the synovium of rat models with anterior cruciate ligament transection-induced KOA.34 HMGB1 is a downstream product of pyroptosis and may have potential inflammatory regulating effects. Chromobox 4 (CBX4), a component of polycomb repressive complex 1, is involved in maintaining cell identity and organ development through gene silencing. The expression of CBX4 has been shown to attenuate the development of OA in mice.35
Ubiquitin C(UBC) is a gene encoding ubiquitin precursor in mammals.32 Conjugation of ubiquitin monomers or polymers can lead to diverse effects within a cell, depending on the residues to which ubiquitin is conjugated.36 A study suggests that UBC could be a promising therapeutic target for ovarian cancer patients with recurrent ubiquitin B(UBB) silencing.37 Hypoxia-inducible factor 1α (HIF-1α) plays a crucial role in preventing apoptosis in chondrocytes. The protein UCHL1 can inhibit apoptosis in chondrocytes by increasing the levels of HIF-1α, which in turn promotes mitophagy and preserves mitochondrial function.38 Bobin Mi reported one of the top 10 central genes of female patients with OA was UBC, and the upregulation of UBC is a therapeutic strategy for maintaining the normal function of cells in the synovium.39
Notch is a single-pass transmembrane cell surface receptor, it regulates differentiation and apoptosis during embryogenesis, as well as in various developmental systems such as neurogenesis and hematopoiesis.40 The proliferation and differentiation of disc cells are dependent on the Notch signaling pathway.41 In a study, miR-485-3p was found to promote proliferation and prevent apoptosis, ECM degradation, inflammation, and oxidative stress in OA chondrocytes by inhibiting Notch2 and the NF-jB pathway.42 The inhibition of Notch2 expression resulted in increased proliferation and reduced apoptosis rate.42 It has been reported that the activation of Notch2 in specific chondrocytes can stimulate the progression of OA.43 Notch2 serves as the receptor for the canonical NF-jB pathway, which is essential for chondrocytes to express MMPs, release ECM components, trigger inflammatory responses, and further contribute to cartilage destruction.44
Serine and arginine repetitive matrix 1(SRRM1), also known as Srm160,45 is an SR-related protein.46 SRRM1 was initially identified as a nuclear matrix antigen highly concentrated in interphase nuclear “speckle” domains, which are enriched in splicing components.47 It functions as a coactivator of both constitutive and exon enhancer-dependent splicing by forming cross-intron interactions with multiple splicing factors bound directly to pre-mRNA.46 SRRM1 has been found to be correlated with tumor malignancy in patients with poor clinical prognosis.48
We first made these systems biology-based analyses to explore complicated OA phenotypes. Steiner minimal tree algorithm results might provide m important insights beyond the conventional single-gene analyses and evaluated. This analysis provided important findings to understand OA progress in a systems biological view.
Although no OA genetic study has yet led to a diagnostic tool or treatment, it has been observed that these genes show a significant increase in expression in seven distinct molecularly defined populations of OA cartilage. Therefore, the current gap between discovery and utility is expected to narrow.
Conclusion
The Steiner tree-based approach identified new biological molecules that are associated with OA genes. Specifically, the CLU, ENO1, SRRM1, UBC, HMGB1, NR3C1, NOTCH2NL, and CBX5 genes were found to be associated with the progression of OA cartilage. In this study, it was observed that these genes exhibited significantly increased expression in seven distinct molecularly defined populations of OA cartilage.
Abbreviations
OA, osteoarthritis; GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; PPI, protein-protein interaction; PINA, Protein Interaction Network Analysis; EC, effector chondrocyte; RegC, regulatory chondrocyte; ProC, proliferative chondrocyte; preHTC, prehypertrophic chondrocyte; FC, fibrocartilage chondrocyte; HTC, hypertrophic chondrocyte; HomC, homeostatic chondrocyte. CLU, clusterin; ENO1, enolase 1; HMGB1, high mobility group box 1; CBX4, chromobox 4; UBC, ubiquitin C; UBB, ubiquitin B; SRRM1, serine and arginine repetitive matrix 1.
Acknowledgments
Thanks to Yuanzhi Cheng from Chengdu Basebiotech Co., Ltd for providing assistance on bioinformatic analysis.
Disclosure
The authors report no conflicts of interest in this work.
References
1. Vos T, Allen C, Arora M, et al. Global regional, and national incidence, prevalence, and years lived with disability for 310 diseases and injuries, 1990–2015: a systematic analysis for the Global Burden of Disease Study 2015. Lancet. 2016;388(10053):1545–1602. doi:10.1016/s0140-6736(16)31678-6
2. Johnson VL, Hunter DJ. The epidemiology of osteoarthritis. Best Pract Res Clin Rheumatol. 2014;28(1):5–15. doi:10.1016/j.berh.2014.01.004
3. Hootman JM, Helmick CG, Barbour KE, et al. Updated projected prevalence of self-reported doctor-diagnosed arthritis and arthritis-attributable activity limitation among US adults, 2015–2040. Arthritis Rheumatol. 2016;68(7):1582–1587. doi:10.1002/art.39692
4. Zhu X, Chan YT, Yung PSH, et al. Subchondral bone remodeling: a therapeutic target for osteoarthritis. Front Cell Dev Biol. 2020;8:607764. doi:10.3389/fcell.2020.607764
5. Ozeki N, Koga H, Sekiya I. Degenerative meniscus in knee osteoarthritis: from pathology to treatment. Life. 2022;12(4). doi:10.3390/life12040603
6. Battistelli M, Favero M, Burini D, et al. Morphological and ultrastructural analysis of normal, injured and osteoarthritic human knee menisci. Eur J Histochem. 2019;63(1). doi:10.4081/ejh.2019.2998
7. Favero M, El-Hadi H, Belluzzi E, et al. Infrapatellar fat pad features in osteoarthritis: a histopathological and molecular study. Rheumatology. 2017;56(10):1784–1793. doi:10.1093/rheumatology/kex287
8. Zhou D, Zhou F, Sheng S, et al. Intra-articular nanodrug delivery strategies for treating osteoarthritis. Drug Discov Today. 2023;28(3):103482. doi:10.1016/j.drudis.2022.103482
9. Felson DT, Lawrence RC, Hochberg MC, et al. Osteoarthritis: new insights. Part 2: treatment approaches. Ann Internal Med. 2000;133(9):726–737. doi:10.7326/0003-4819-133-9-200011070-00015
10. Mazur CM, Woo JJ, Yee CS, et al. Osteocyte dysfunction promotes osteoarthritis through MMP13-dependent suppression of subchondral bone homeostasis. Bone Res. 2019;7:34. doi:10.1038/s41413-019-0070-y
11. Haraden CA, Huebner JL, Hsueh MF, et al. Synovial fluid biomarkers associated with osteoarthritis severity reflect macrophage and neutrophil related inflammation. Arthritis Res Ther. 2019;21(1):146. doi:10.1186/s13075-019-1923-x
12. Tachmazidou I, Hatzikotoulas K, Southam L, et al. Identification of new therapeutic targets for osteoarthritis through genome-wide analyses of UK Biobank data. Nat Genet. 2019;51(2):230–236. doi:10.1038/s41588-018-0327-1
13. Son YO, Kim HE, Choi WS, et al. RNA-binding protein ZFP36L1 regulates osteoarthritis by modulating members of the heat shock protein 70 family. Nat Commun. 2019;10(1):77. doi:10.1038/s41467-018-08035-7
14. Shen P, Serve S, Wu P, et al. NOS inhibition reverses TLR2-induced chondrocyte dysfunction and attenuates age-related osteoarthritis. Proc Natl Acad Sci U S A. 2023;120(29):e2207993120. doi:10.1073/pnas.2207993120
15. Wang T, Liang Y, Li H, et al. Single nucleotide polymorphisms and osteoarthritis: an overview and a meta-analysis. Medicine. 2016;95(7):e2811. doi:10.1097/md.0000000000002811
16. Wang J, Fan Q, Yu T, et al. Identifying the hub genes and immune cell infiltration in synovial tissue between osteoarthritic and rheumatoid arthritic patients by bioinformatic approach. Curr Pharm Des. 2022;28(6):497–509. doi:10.2174/1381612827666211104154459
17. Ji Q, Zheng Y, Zhang G, et al. Single-cell RNA-seq analysis reveals the progression of human osteoarthritis. Ann Rheumatic Dis. 2019;78(1):100–110. doi:10.1136/annrheumdis-2017-212863
18. Shannon P, Markiel A, Ozier O, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–2504. doi:10.1101/gr.1239303
19. Wang Y, Chen Y, Wei Y. Osteoarthritis animal models for biomaterial-assisted osteochondral regeneration. Biomater Transl. 2022;3(4):264–279. doi:10.12336/biomatertransl.2022.04.006
20. Tamaddon M, Gilja H, Wang L, et al. Osteochondral scaffolds for early treatment of cartilage defects in osteoarthritic joints: from bench to clinic. Biomater Transl. 2020;1(1):3–17. doi:10.3877/cma.j.issn.2096-112X.2020.01.002
21. Zhang H, Wang L, Cui J, et al. Maintaining hypoxia environment of subchondral bone alleviates osteoarthritis progression. Sci Adv. 2023;9(14):eabo7868. doi:10.1126/sciadv.abo7868
22. Tarquini C, Pucci S, Scioli MG, et al. Clusterin exerts a cytoprotective and antioxidant effect in human osteoarthritic cartilage. Aging. 2020;12(11):10129–10146. doi:10.18632/aging.103310
23. Poon S, Easterbrook-Smith SB, Rybchyn MS, et al. Clusterin is an ATP-independent chaperone with very broad substrate specificity that stabilizes stressed proteins in a folding-competent state. Biochemistry. 2000;39(51):15953–15960. doi:10.1021/bi002189x
24. Kalvaityte U, Matta C, Bernotiene E, et al. Exploring the translational potential of clusterin as a biomarker of early osteoarthritis. J Orthopa Transl. 2022;32:77–84. doi:10.1016/j.jot.2021.10.001
25. Matta C, Fellows CR, Quasnichka H, et al. Clusterin secretion is attenuated by the proinflammatory cytokines interleukin-1β and tumor necrosis factor-α in models of cartilage degradation. J Orthopaedic Res. 2021;39(5):1017–1029. doi:10.1002/jor.24814
26. Fernández-Puente P, González-Rodríguez L, Calamia V, et al. Analysis of endogenous peptides released from osteoarthritic cartilage unravels novel pathogenic markers. Mol Cell Proteomics. 2019;18(10):2018–2028. doi:10.1074/mcp.RA119.001554
27. Bae S, Kim H, Lee N, et al. α-Enolase expressed on the surfaces of monocytes and macrophages induces robust synovial inflammation in rheumatoid arthritis. J Immunol. 2012;189(1):365–372. doi:10.4049/jimmunol.1102073
28. Lee JY, Kang MJ, Choi JY, et al. Apolipoprotein B binds to enolase-1 and aggravates inflammation in rheumatoid arthritis. Ann Rheumatic Dis. 2018;77(10):1480–1489. doi:10.1136/annrheumdis-2018-213444
29. Ma WJ, Guo X, Liu JT, et al. Proteomic changes in articular cartilage of human endemic osteoarthritis in China. Proteomics. 2011;11(14):2881–2890. doi:10.1002/pmic.201000636
30. Ruiz-Romero C, Carreira V, Rego I, et al. Proteomic analysis of human osteoarthritic chondrocytes reveals protein changes in stress and glycolysis. Proteomics. 2008;8(3):495–507. doi:10.1002/pmic.200700249
31. Wang Y, Shen S, Li Z, et al. MIR-140-5p affects chondrocyte proliferation, apoptosis, and inflammation by targeting HMGB1 in osteoarthritis. Inflammat Res. 2020;69(1):63–73. doi:10.1007/s00011-019-01294-0
32. Lund PK, Moats-Staats BM, Simmons JG, et al. Nucleotide sequence analysis of a cDNA encoding human ubiquitin reveals that ubiquitin is synthesized as a precursor. J Biol Chem. 1985;260(12):7609–7613.
33. Wagner G, Lehmann C, Bode C, et al. High mobility group box 1 protein in osteoarthritic knee tissue and chondrogenic progenitor cells: an ex vivo and in vitro study. Cartilage. 2021;12(4):484–495. doi:10.1177/1947603519835897
34. Xiao Y, Ding L, Yin S, et al. Relationship between the pyroptosis of fibroblast‑like synoviocytes and HMGB1 secretion in knee osteoarthritis. Mole Med Rep. 2021;23(2). doi:10.3892/mmr.2020.11736
35. Ren X, Hu B, Song M, et al. Maintenance of nucleolar homeostasis by CBX4 alleviates senescence and osteoarthritis. Cell Rep. 2019;26(13):3643–56.e7. doi:10.1016/j.celrep.2019.02.088
36. Hershko A, Ciechanover A. The ubiquitin system. Annu. Rev. Biochem. 1998;67:425–479. doi:10.1146/annurev.biochem.67.1.425
37. Kedves AT, Gleim S, Liang X, et al. Recurrent ubiquitin B silencing in gynecological cancers establishes dependence on ubiquitin C. J Clin Invest. 2017;127(12):4554–4568. doi:10.1172/jci92914
38. Yan Q, Shi S, Ge Y, et al. UCHL1 alleviates apoptosis in chondrocytes via upregulation of HIF‑1α‑mediated mitophagy. Int J Mol Med. 2023;52(4). doi:10.3892/ijmm.2023.5302
39. Mi B, Liu G, Zhou W, et al. Identification of genes and pathways in the synovia of women with osteoarthritis by bioinformatics analysis. Mole Med Rep. 2018;17(3):4467–4473. doi:10.3892/mmr.2018.8429
40. Pillai S, Cariappa A. The follicular versus marginal zone B lymphocyte cell fate decision. Nat Rev Immunol. 2009;9(11):767–777. doi:10.1038/nri2656
41. Dong Y, Jesse AM, Kohn A, et al. RBPjkappa-dependent Notch signaling regulates mesenchymal progenitor cell proliferation and differentiation during skeletal development. Development. 2010;137(9):1461–1471. doi:10.1242/dev.042911
42. Zhou Y, Zhao Z, Yan L, et al. MiR-485-3p promotes proliferation of osteoarthritis chondrocytes and inhibits apoptosis via Notch2 and the NF-kappaB pathway. Immunopharmacol Immunotoxicol. 2021;43(3):370–379. doi:10.1080/08923973.2021.1918150
43. Zanotti S, Yu J, Bridgewater D, et al. Mice harboring a Hajdu Cheney Syndrome mutation are sensitized to osteoarthritis. Bone. 2018;114:198–205. doi:10.1016/j.bone.2018.06.020
44. Goldring MB, Otero M. Inflammation in osteoarthritis. Curr Opin Rheumatol. 2011;23(5):471–478. doi:10.1097/BOR.0b013e328349c2b1
45. Chwalenia K, Qin F, Singh S, et al. A cell-based splicing reporter system to identify regulators of cis-splicing between adjacent genes. Nucleic Acids Res. 2019;47(4):e24. doi:10.1093/nar/gky1288
46. McCracken S, Lambermon M, Blencowe BJ. SRm160 splicing coactivator promotes transcript 3’-end cleavage. Mol Cell Biol. 2002;22(1):148–160. doi:10.1128/MCB.22.1.148-160.2002
47. McCracken S, Longman D, Marcon E, et al. Proteomic analysis of SRm160-containing complexes reveals a conserved association with cohesin. J Biol Chem. 2005;280(51):42227–42236. doi:10.1074/jbc.M507410200
48. Cheng C, Sharp PA. Regulation of CD44 alternative splicing by SRm160 and its potential role in tumor cell invasion. Mol Cell Biol. 2006;26(1):362–370. doi:10.1128/MCB.26.1.362-370.2006
© 2024 The Author(s). This work is published and licensed by Dove Medical Press Limited. The
full terms of this license are available at https://www.dovepress.com/terms.php
and incorporate the Creative Commons Attribution
- Non Commercial (unported, v3.0) License.
By accessing the work you hereby accept the Terms. Non-commercial uses of the work are permitted
without any further permission from Dove Medical Press Limited, provided the work is properly
attributed. For permission for commercial use of this work, please see paragraphs 4.2 and 5 of our Terms.
Recommended articles
Bioinformatics and Integrative Experimental Method to Identifying and Validating Co-Expressed Ferroptosis-Related Genes in OA Articular Cartilage and Synovium
Ma J, Yu P, Ma S, Li J, Wang Z, Hu K, Su X, Zhang B, Cheng S, Wang S
Journal of Inflammation Research 2024, 17:957-980
Published Date: 12 February 2024