Integrative Analysis of Three Novel Competing Endogenous RNA Biomarkers with a Prognostic Value in Lung Adenocarcinoma

Abstract

Increasing evidence has shown competitive endogenous RNAs (ceRNAs) play key roles in numerous cancers. Nevertheless, the ceRNA network that can predict the prognosis of lung adenocarcinoma (LUAD) is still lacking. The aim of the present study was to identify the prognostic value of key ceRNAs in lung tumorigenesis. Differentially expressed (DE) RNAs were identified between LUAD and adjacent normal samples by limma package in R using The Cancer Genome Atlas database (TCGA). Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway function enrichment analysis was performed using the clusterProfiler package in R. Subsequently, the LUAD ceRNA network was established in three steps based on ceRNA hypothesis. Hub RNAs were identified using degree analysis methods based on Cytoscape plugin cytoHubba. Multivariate Cox regression analysis was implemented to calculate the risk score using the candidate ceRNAs and overall survival information. The survival differences between the high-risk and low-risk ceRNA groups were determined by the Kaplan-Meier and log-rank test using survival and survminer package in R. A total of 2,989 mRNAs, 185 lncRNAs, and 153 miRNAs were identified. GO and KEGG pathway function enrichment analysis showed that DE mRNAs were mainly associated with “sister chromatid segregation,” “regulation of angiogenesis,” “cell adhesion molecules (CAMs),” “cell cycle,” and “ECM-receptor interaction.” LUAD-related ceRNA network was constructed, which comprised of 54 nodes and 78 edges. Top ten hub RNAs (hsa-miR-374a-5p, hsa-miR-374b-5p, hsa-miR-340-5p, hsa-miR-377-3p, hsa-miR-21-5p, hsa-miR-326, SNHG1, RALGPS2, and PITX2) were identified according to their degree. Kaplan-Meier survival analyses demonstrated that hsa-miR-21-5p and RALGPS2 had a significant prognostic value. Finally, we found that a high risk of three novel ceRNA interactions (SNHG1-hsa-miR-21-5p-RALGPS2, SNHG1-hsa-miR-326-RALGPS2, and SNHG1-hsa-miR-377-3p-RALGPS2) was positively associated with worse prognosis. Three novel ceRNAs (SNHG1-hsa-miR-21-5p-RALGPS2, SNHG1-hsa-miR-326-RALGPS2, and SNHG1-hsa-miR-377-3p-RALGPS2) might be potential biomarkers for the prognosis and treatment of LUAD.