Journal of Jilin University(Medicine Edition) ›› 2022, Vol. 48 ›› Issue (1): 154-162.doi: 10.13481/j.1671-587X.20220119

• Research in clinical medicine • Previous Articles     Next Articles

Bioinformatics analysis on miRNA-mRNA regulatory networks based on fusion genes acting in rhabdomyosarcoma

Zhijuan ZHAO1,Lian MENG1,Chunxia LIU1,2()   

  1. 1.Department of Pathology,First Affiliated Hospital,College of Medical Sciences,Shihezi University,Shihezi 832002,China
    2.Department of Pathology,Second Affiliated Hospital,Guangzhou Medical University,Guangzhou 510260,China
  • Received:2021-05-21 Online:2022-01-28 Published:2022-01-17
  • Contact: Chunxia LIU E-mail:liuliu2239@sina.com

Abstract: Objective

To construct the potential miRNA-mRNA regulatory network of fusion genes in the pathogenesis of rhabdomyosarcoma (RMS)with bioinformatics analysis, and to provide a new direction for the study of RMS.

Methods

GEO2R analysis tool was used to screen the differentially expressed genes (DEGs) and differentially expressed miRNA in fusion-gene-positive and -negative RMS tissues; the target genes of differentially expressed miRNA were predicted; and the target genes on the basis of overlapping DEGs and target genes were screened out. DAVID database was implemented for the GO and KEGG enrichment analysis of target genes. STRING and Cytoscape software were utilized to construct the protein protein interation (PPI) network, the Top10 hub genes were screened out and the molecular regulation networks of hub genes and miRNA were constructed, and the Kaplan-Meier survival curves of the sarcoma patients were drawn.

Results

A total of 891 DEGs (P<0.01,|logFC|≥1) and 14 differentially expressed miRNAs (P<0.05, |logFC|≥3) were screened out. Moreover, 1 654 target genes differentially expressing miRNAs were predicted, and 115 target genes were identified by overlapping the target genes and DEGs. The GO enrichment analysis revealed that the target genes were mainly enriched in the positive regulation of cell proliferation, cell surface, protein binding, and other biological processes. The KEGG signaling pathway analysis revealed that the target genes were mainly enriched in the extracellular matrix-receptor interaction, proteoglycan in cancer, and other pathways. The Top10 hub genes screened out from the PPI network were ERBB4, PTPRD, IRS1, ADAM10, YAP1, TFAP2A, CADM1, ELAVL2, SNAI1, and ERRFI1. The survival analysis showed that the increased expression levels of PTPRD (HR=1.79, P=0.005),ADAM10 (HR =1.94,P=0.010), ELAVL2 (HR = 1.56,P=0.031),and ERRFI1 (HR=2.05,P=0.005) were related to the bad prognosis of the sarcoma patients.

Conclusion

The miRNA-mRNA network is constructed by selecting the hub genes and corresponding miRNA involed in the pathogenesis of RMS, and the study provides a theoretical basis for further study of the relationship between fusion genes and RMS.

Key words: Rhabdomyosarcoma, Fusion gene, Bioinformatics, MiRNA-mRNA regulatory network, Differentially expressed genes

CLC Number: 

  • R738.6