Recent evidence has indicated that long non-coding RNAs (lncRNAs) can function as competing endogenous RNAs (ceRNAs) to modulate mRNAs expression by sponging microRNAs (miRNAs). However, the specific mechanism and function of lncRNA-miRNA-mRNA regulatory network in non-small cell lung cancer (NSCLC) remains unclear.
We constructed a lung cancer related lncRNA-mRNA network (LCLMN) by integrating differentially expressed genes (DEGs) with miRNA-target interactions. We further performed topological feature analysis and random walk with restart (RWR) analysis of LCLMN. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed to investigate the target DEGs in LCLMN. The expression levels of significant lncRNAs in NSCLC were validated by quantitative real-time PCR (RT-qPCR). The prognostic value of the potential lncRNA was evaluated by Kaplan-Meier analysis.
A total of 33 lncRNA nodes, 580 mRNA nodes and 2105 edges were identified from LCLMN. Based on functional enrichment analysis and co-expression analysis, lncRNA EPB41L4A-AS1 was demonstrated to be correlated with the tumorigenesis of NSCLC. RT-qPCR results confirmed that the expression levels of lncRNA EPB41L4A-AS1 in NSCLC tissues were downregulated compared with adjacent non-cancerous tissues. Kaplan-Meier analysis showed that high expression of lncRNA EPB41L4A-AS1 was associated with better overall survival (OS) in NSCLC patients. Further investigation identified that high expression levels of COL4A3BP, CDS2, PURA, PDCD6IP, and TMEM245 were also correlated with better OS in NSCLC patients.
In this study, we constructed a lncRNA-miRNA-mRNA ceRNA network to investigate potential prognostic biomarkers for NSCLC. We found that lncRNA EPB41L4A-AS1 could function as a regulator in the pathogenesis of NSCLC.