Bioinformatics Analysis of MiRNA and MRNA Expression Profiles to Reveal the Key MiRNAs and Genes in Osteoarthritis

Overview

Journal J Orthop Surg Res

Publisher Biomed Central

Specialty Orthopedics

Date 2021 Jan 20

PMID 33468167

Citations 13

Authors

Pei-Yan Huang

Jun-Guo Wu

Jun Gu

Tie-Qi Zhang

Ling-Feng Li

Si-Qun Wang

Minghai Wang

Affiliations

Soon will be listed here.

Abstract

Background: Osteoarthritis (OA) is a chronic degenerative joint disease and the most frequent type of arthritis. This study aimed to identify the key miRNAs and genes associated with OA progression.

Methods: The GSE105027 (microRNA [miRNA/miR] expression profile; 12 OA samples and 12 normal samples) and GSE48556 (messenger RNA [mRNA] expression profile; 106 OA samples and 33 normal samples) datasets were selected from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) and miRNAs (DEMs) were analyzed using the limma and ROCR packages in R, respectively. The target genes that negatively correlated with the DEMs were predicted, followed by functional enrichment analysis and construction of the miRNA-gene and miRNA-transcription factor (TF)-gene regulatory networks. Additionally, key miRNAs and genes were screened, and their expression levels were verified by real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR).

Results: A total of 1696 DEGs (739 upregulated and 957 downregulated) and 108 DEMs (56 upregulated and 52 downregulated) were identified in the OA samples. Furthermore, 56 target genes that negatively correlated with the DEMs were predicted and found to be enriched in three functional terms (e.g., positive regulation of intracellular protein transport) and three pathways (e.g., human cytomegalovirus infection). In addition, three key miRNAs (miR-98-5p, miR-7-5p, and miR-182-5p) and six key genes (murine double minute 2, MDM2; glycogen synthase kinase 3-beta, GSK3B; transmembrane P24-trafficking protein 10, TMED10; DDB1 and CUL4-associated factor 12, DCAF12; caspase 3, CASP3; and ring finger protein 44, RNF44) were screened, among which the miR-7-5p → TMED10/DCAF12, miR-98-5p → CASP3/RNF44, and miR-182-5p → GSK3B pairs were observed in the regulatory network. Moreover, the expression levels of TMED10, miR-7-5p, CASP3, miR-98-5p, GSK3B, and miR-182-5p showed a negative correlation with qRT-PCR verification.

Conclusion: MiR-98-5p, miR-7-5p, miR-182-5p, MDM2, GSK3B, TMED10, DCAF12, CASP3, and RNF44 may play critical roles in OA progression.

Citing Articles

Exploration of the mechanism of tetramethoxyflavone in treating osteoarthritis based on network pharmacology and molecular docking.

Chen P, Ye B, Lin C, Zhang C, Chen J, Li L Tzu Chi Med J. 2025; 37(1):99-108.

PMID: 39850388 PMC: 11753527. DOI: 10.4103/tcmj.tcmj_77_24.

STUDY ON THE ROLE AND MECHANISM OF MICRORNA-650/WNT1 IN THE REPAIR OF ARTICULAR CARTILAGE INJURY.

Liu H, Wang Y, Wang S, Yang B, Sun D, Han S Acta Ortop Bras. 2024; 32(4):e278218.

PMID: 39386291 PMC: 11460656. DOI: 10.1590/1413-785220243204e278218.

SP1‑mediated ADAMTS5 transcription promotes IL‑1β‑induced chondrocyte injury via Wnt/β‑catenin pathway in osteoarthritis.

Hao X, Wu X Mol Med Rep. 2024; 30(2).

PMID: 38940327 PMC: 11228694. DOI: 10.3892/mmr.2024.13273.

MiR-98-5p plays suppressive effects on IL-1β-induced chondrocyte injury associated with osteoarthritis by targeting CASP3.

Lv H, Liu P, Hu H, Li X, Li P J Orthop Surg Res. 2024; 19(1):239.

PMID: 38615043 PMC: 11015643. DOI: 10.1186/s13018-024-04628-9.

Hsa_circular RNA_0045474 Facilitates Osteoarthritis Via Modulating microRNA-485-3p and Augmenting Transcription Factor 4.

Zhang Z, Yu P, Bai L Mol Biotechnol. 2024; 66(5):1174-1187.

PMID: 38206529 DOI: 10.1007/s12033-023-01019-z.

References

Nakano N, Tsuchiya Y, Kako K, Umezaki K, Sano K, Ikeno S . TMED10 Protein Interferes with Transforming Growth Factor (TGF)-β Signaling by Disrupting TGF-β Receptor Complex Formation. J Biol Chem. 2017; 292(10):4099-4112. PMC: 5354491. DOI: 10.1074/jbc.M116.769109. View

Wang J, Chen L, Jin S, Lin J, Zheng H, Zhang H . Altered expression of microRNA-98 in IL-1β-induced cartilage degradation and its role in chondrocyte apoptosis. Mol Med Rep. 2017; 16(3):3208-3216. PMC: 5547958. DOI: 10.3892/mmr.2017.7028. View

Soyocak A, Kurt H, Ozgen M, Cosan D, Colak E, Gunes H . miRNA-146a, miRNA-155 and JNK expression levels in peripheral blood mononuclear cells according to grade of knee osteoarthritis. Gene. 2017; 627:207-211. DOI: 10.1016/j.gene.2017.06.027. View

Chou C, Shrestha S, Yang C, Chang N, Lin Y, Liao K . miRTarBase update 2018: a resource for experimentally validated microRNA-target interactions. Nucleic Acids Res. 2017; 46(D1):D296-D302. PMC: 5753222. DOI: 10.1093/nar/gkx1067. View

Li G, Yin J, Gao J, Cheng T, Pavlos N, Zhang C . Subchondral bone in osteoarthritis: insight into risk factors and microstructural changes. Arthritis Res Ther. 2013; 15(6):223. PMC: 4061721. DOI: 10.1186/ar4405. View

Trachana V, Ntoumou E, Anastasopoulou L, Tsezou A . Studying microRNAs in osteoarthritis: Critical overview of different analytical approaches. Mech Ageing Dev. 2018; 171:15-23. DOI: 10.1016/j.mad.2018.02.006. View

Livak K, Schmittgen T . Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2002; 25(4):402-8. DOI: 10.1006/meth.2001.1262. View

Oka S, Furukawa H, Shimada K, Hashimoto A, Komiya A, Fukui N . Plasma miRNA expression profiles in rheumatoid arthritis associated interstitial lung disease. BMC Musculoskelet Disord. 2017; 18(1):21. PMC: 5244611. DOI: 10.1186/s12891-017-1389-4. View

Quinlan A, Hall I . BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010; 26(6):841-2. PMC: 2832824. DOI: 10.1093/bioinformatics/btq033. View

10.

McAlindon T, Bannuru R, Sullivan M, Arden N, Berenbaum F, Bierma-Zeinstra S . OARSI guidelines for the non-surgical management of knee osteoarthritis. Osteoarthritis Cartilage. 2014; 22(3):363-88. DOI: 10.1016/j.joca.2014.01.003. View

11.

Si H, Zeng Y, Zhou Z, Pei F, Lu Y, Cheng J . Expression of miRNA-140 in Chondrocytes and Synovial Fluid of Knee Joints in Patients with Osteoarthritis. Chin Med Sci J. 2017; 31(4):207-212. DOI: 10.1016/s1001-9294(17)30002-0. View

12.

Shen J, Chen D . Recent progress in osteoarthritis research. J Am Acad Orthop Surg. 2014; 22(7):467-8. PMC: 4124725. DOI: 10.5435/JAAOS-22-07-467. View

13.

March L, R Smith E, Hoy D, Cross M, Sanchez-Riera L, Blyth F . Burden of disability due to musculoskeletal (MSK) disorders. Best Pract Res Clin Rheumatol. 2014; 28(3):353-66. DOI: 10.1016/j.berh.2014.08.002. View

14.

Yu G, Wang L, Han Y, He Q . clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012; 16(5):284-7. PMC: 3339379. DOI: 10.1089/omi.2011.0118. View

15.

Fu J, Tang W, Du P, Wang G, Chen W, Li J . Identifying microRNA-mRNA regulatory network in colorectal cancer by a combination of expression profile and bioinformatics analysis. BMC Syst Biol. 2012; 6:68. PMC: 3418553. DOI: 10.1186/1752-0509-6-68. View

16.

Taranto E, Xue J, Lacey D, Hutchinson P, Smith M, Morand E . Detection of the p53 regulator murine double-minute protein 2 in rheumatoid arthritis. J Rheumatol. 2005; 32(3):424-9. View

17.

Giordano L, Della Porta G, Peretti G, Maffulli N . Therapeutic potential of microRNA in tendon injuries. Br Med Bull. 2020; 133(1):79-94. DOI: 10.1093/bmb/ldaa002. View

18.

. Gene Ontology Consortium: going forward. Nucleic Acids Res. 2014; 43(Database issue):D1049-56. PMC: 4383973. DOI: 10.1093/nar/gku1179. View

19.

Sing T, Sander O, Beerenwinkel N, Lengauer T . ROCR: visualizing classifier performance in R. Bioinformatics. 2005; 21(20):3940-1. DOI: 10.1093/bioinformatics/bti623. View

20.

Li J, Callegari S, Masucci M . The Epstein-Barr virus miR-BHRF1-1 targets RNF4 during productive infection to promote the accumulation of SUMO conjugates and the release of infectious virus. PLoS Pathog. 2017; 13(4):e1006338. PMC: 5413087. DOI: 10.1371/journal.ppat.1006338. View