Integrated analysis of lncRNA-miRNA-mRNA ceRNA network and the potential prognosis indicators in sarcomas
Abstract The ceRNA network has been demonstrated to play crucial roles in multiple biological processes and the development of neoplasms, which have the potential to become diagnostic and prognosis markers and therapeutic targets. In this work, we comparing the expression profiles between sarcoma identified differentially expressed genes (DEGs), lncRNAs (DELs) and miRNAs (DEMs) in sarcomas and normal tissue samples in GEO datasets. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were applied to investigate the major functions of the overlapping DEGs. Then, lncRNA-miRNA interactions and miRNA-mRNA interactions were predicted, and a ceRNA regulatory network was constructed. In addition, the mRNAs included in ceRNA network were used to construct the protein-protein interactions network, and the survival analysis of sarcomas was performed according to the biomarkers included in the ceRNA network. According to the RNA sequence data from GEO dataset, 1296 DEGs were identified in sarcoma samples by combining the GO and Pathway enrichment analysis, 338 DELs were discovered after re-annotating the probes, and 36 DEGs were ascertained through intersecting two different expression miRNAs sets. Further, 448 miRNA-mRNA interactions and 454 miRNA-lncRNA interactions were obtained through target gene prediction, and then, we constructed a lncRNA-miRNA-mRNA ceRNA network containing 9 miRNAs, 69 lncRNAs and 113 mRNAs. PPI network showed that the hub up-regulated nodes include IGF1, PRKCB and GNAI3, and the hub down-regulated nodes include AR, CYCS and PPP1CB. Survival analysis revealed that the expression levels of 12 RNAs involved in the ceRNA network were associated with overall survival of sarcoma patients. Our study showed that the ceRNA network in sarcomas based on that lncRNA could serve as ceRNA and discovered the potential indicators for prognosis of sarcoma patients.