Dataset Information

Bioinformatics analysis of miRNA and mRNA expression profiles to reveal the key miRNAs and genes in osteoarthritis.

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.

SUBMITTER: Huang PY

PROVIDER: S-EPMC7814623 | biostudies-literature | 2021 Jan

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Bioinformatics analysis of miRNA and mRNA expression profiles to reveal the key miRNAs and genes in osteoarthritis.

Huang Pei-Yan PY Wu Jun-Guo JG Gu Jun J Zhang Tie-Qi TQ Li Ling-Feng LF Wang Si-Qun SQ Wang Minghai M

Journal of orthopaedic surgery and research 20210119 1

<h4>Background</h4>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.<h4>Methods</h4>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 ge ...[more]

PMID: 33468167

Similar Datasets

Project description:Background: Osteoarthritis (OA) is one of the main causes of disability in the elderly population, accompanied by a series of underlying pathologic changes, such as cartilage degradation, synovitis, subchondral bone sclerosis, and meniscus injury. The present study aimed to identify key genes, signaling pathways, and miRNAs in knee OA associated with the entire joint components, and to explain the potential mechanisms using computational analysis. Methods: The differentially expressed genes (DEGs) in cartilage, synovium, subchondral bone, and meniscus were identified using the Gene Expression Omnibus 2R (GEO2R) analysis based on dataset from GSE43923, GSE12021, GSE98918, and GSE51588, respectively and visualized in Volcano Plot. Venn diagram analyses were performed to identify the overlapping DEGs (overlapping DEGs) that expressed in at least two types of tissues mentioned above. Gene Ontology (GO) enrichment analysis, Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, protein-protein interaction (PPI) analysis, and module analysis were conducted. Furthermore, qRT-PCR was performed to validate above results using our clinical specimens. Results: As a result, a total of 236 overlapping DEGs were identified, of which 160 were upregulated and 76 were downregulated. Through enrichment analysis and constructing the PPI network and miRNA-mRNA network, knee OA-related key genes, such as HEY1, AHR, VEGFA, MYC, and CXCL12 were identified. Clinical validation by qRT-PCR experiments further supported above computational results. In addition, knee OA-related key miRNAs such as miR-101, miR-181a, miR-29, miR-9, and miR-221, and pathways such as Wnt signaling, HIF-1 signaling, PI3K-Akt signaling, and axon guidance pathways were also identified. Among above identified knee OA-related key genes, pathways and miRNAs, genes such as AHR, HEY1, MYC, GAP43, and PTN, pathways like axon guidance, and miRNAs such as miR-17, miR-21, miR-155, miR-185, and miR-1 are lack of research and worthy for future investigation. Conclusion: The present informatic study for the first time provides insight to the potential therapeutic targets of knee OA by comprehensively analyzing the overlapping genes differentially expressed in multiple joint components and their relevant signaling pathways and interactive miRNAs.

Project description:Gastric cancer (GC) is one of the most common malignancies in the world, with morbidity and mortality ranking second among all cancers. Accumulating evidences indicate that circular RNAs (circRNAs) are closely correlated with tumorigenesis. However, the mechanisms of circRNAs still remain unclear. This study is aimed at determining hub genes and circRNAs and analyzing their potential biological functions in GC. Expression profiles of mRNAs and circRNAs were downloaded from the Gene Expression Omnibus (GEO) data sets of GC and paracancer tissues. Differentially expressed genes (DEGs) and differentially expressed circRNAs (DE-circRNAs) were identified. The target miRNAs of DE-circRNAs and the bidirectional interaction between target miRNAs and DEGs were predicted. Functional analysis was performed, and the protein-protein interaction (PPI) network and the circRNA-miRNA-mRNA network were established. A total of 456 DEGs and 2 DE-circRNAs were identified with 3 mRNA expression profiles and 2 circRNA expression profiles. GO analysis indicated that DEGs were mainly enriched in extracellular matrix and cell adhesion, and KEGG confirmed that DEGs were mainly associated with focal adhesion, the PI3K-Akt signaling pathway, extracellular matrix- (ECM)- receptor interaction, and gastric acid secretion. 15 hub DEGs (BGN, COL1A1, COL1A2, FBN1, FN1, SPARC, SPP1, TIMP1, UBE2C, CCNB1, CD44, CXCL8, COL3A1, COL5A2, and THBS1) were identified from the PPI network. Furthermore, the survival analysis indicate that GC patients with a high expression of the following 9 hub DEGs, namely, BGN, COL1A1, COL1A2, FBN1, FN1, SPARC, SPP1, TIMP1, and UBE2C, had significantly worse overall survival. The circRNA-miRNA-mRNA network was constructed based on 1 circRNA, 15 miRNAs, and 45 DEGs. In addition, the 45 DEGs included 5 hub DEGs. These results suggested that hub DEGs and circRNAs could be implicated in the pathogenesis and development of GC. Our findings provide novel evidence on the circRNA-miRNA-mRNA network and lay the foundation for future research of circRNAs in GC.

Project description:IntroductionCircRNA is closely correlated with a wide variety of processes of acute myeloid leukemia (AML), whereas the novel circRNAs, their molecular mechanism and the specific function they played in AML should be explored in depth.MethodsThe microarray chip data of AML patients and normal samples in the Gene Expression Omnibus (GEO) database were selected to differentially expressed (DE) circRNA, miRNA, and mRNA genes. The miRNA gene was the intersection of the circRNA target gene predicted using CSCD and the miRNA gene screened from AML patients, while the mRNA gene was the intersection of the target gene mRNA of miRNA predicted using miRanda and miRTarBase software and the mRNA gene screened from AML patients. The hub mRNAs related to survival were further screened through Cox proportional hazard regression. CircRNA/miRNA/mRNA interaction network was constructed by using Cytoscape software.10 circRNAs and 6 miRNAs in bone marrow mononuclear cells (BMMNCs) of AML patients (n=43) and healthy controls (n=35) were determined by RT-qPCR. Correlations between them were analyzed by Pearson correlation coefficient.Results10 circRNAs, 6 miRNAs, and 33 mRNAs were identified. Subsequently, the network of circRNAs, miRNAs, and hub genes was built using Cystoscope. Four key circRNAs, seven hub genes and their regulatory pathways were identified. The result of RT-qPCRs showed that hsa_circ_0009581 and hsa_circ_0005273 were significantly upregulated in AML patients while hsa_circ_0000497 and hsa_circ_0001947 were significantly downregulated. Hsa-miR-150-5p was significantly downregulated; hsa-miR-454-3p was upregulated in AML patients. Hsa_circ_0009581 and hsa-miR-150-5p; hsa_ circ_0001947 and hsa-miR-454-3p were inversely correlated using Pearson's correlation coefficient.ConclusionThis study suggests that differentially expressed circRNAs take on a critical significance to AML development and may be the effective therapeutic targets. We suppose that hsa_circ_0009581 promotes leukemia development through hsa-miR-150-5p and hsa_circ_0001947 through hsa-miR-454-3p. hsa_circ_0001947 and hsa_circ_0009581 may provide new directions in the pathogenesis of AML.

Dataset Information

Bioinformatics analysis of miRNA and mRNA expression profiles to reveal the key miRNAs and genes in osteoarthritis.

Background

Methods

Results

Conclusion

Publications

Bioinformatics analysis of miRNA and mRNA expression profiles to reveal the key miRNAs and genes in osteoarthritis.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets