Identification of ceRNA network based on a RNA-seq shows prognostic lncRNA biomarkers in human lung adenocarcinoma

Previous studies have emphasized the significant functions of long non-coding RNAs (lncRNAs) as competing endogenous RNAs (ceRNAs) in tumor biology. However, the functions of certain cancer lncRNAs in the lncRNA-related ceRNA network in lung adenocarcinoma (LUAD) are unknown. A systematic and integrative survey of RNA-seq data from The Cancer Genome Atlas (TCGA) was performed to identify candidate lncRNAs for the prognosis of LUAD. In total, 20,502 genes that contain 181 lncRNAs were evaluated in a cohort of 570 LUAD cases. Initially, 6,280 differentially expressed genes (fold-change >2, P<0.05) were obtained using R package, which includes 75 lncRNAs. Next, by univariate regression and multivariate Cox proportional hazards analysis, 32 genes were associated with survival in LUAD. Using these 29 mRNAs and 3 lncRNAs, a prognosis index (PI) was calculated to accurately estimate the survival in LUAD: PI=∑exp_{risk gene} × HRrisk gene. Furthermore, the 32-gene signature was an independent prognostic indicator for LUAD (HR >1; P<0.05, by multivariate analysis). Weighted gene co-expression network analysis (WGCNA) of three risk lncRNAs-FAM138B, NHEG1 and TLX1NB-was performed, based on the P-values of the associated genes, and the top 27 miRNAs that bound to these lncRNAs were predicted by Miranda as target miRNAs. Next, these target miRNAs were transferred to the TarBase, miRTarBase, miRecards and starBase v2.0 databases to obtain their target genes. According to the previous miRNA-mRNA and miRNA-lncRNA data, three lncRNA-miRNA-mRNA ceRNA networks were established, based on the 29 prognostic mRNAs, forming a regulatory network in LUAD. The present study provided insight into the lncRNA-related ceRNA network in LUAD and has identified potential diagnostic and prognostic biomarkers.