Abstract:Background Lung adenocarcinoma (LUAD) is the most common histological subtype of lung cancer worldwide. Until now, the molecular mechanisms underlying LUAD progression have not been fully explained. This study aimed to identify a competing endogenous RNA (ceRNA) network in LUAD. Methods Differentially expressed lncRNAs (DELs), miRNAs (DEMs), and mRNAs (DEGs) were identified from The Cancer Genome Atlas (TCGA) database with a |log2FC| > 1.0 and a false discovery rate (FDR) < 0.05. Then, these DELs, DEMs, and DEGs were used to construct the initial ceRNA network. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), protein-protein interaction (PPI) network, and survival analyses were performed to analyse these DEGs involved in the ceRNA network. Subsequently, the drug-gene interaction database (DGIdb) was utilized to select candidate LUAD drugs interacting with significant DEGs. Then, lasso-penalized Cox regression and multivariate Cox regression models were used to construct the risk score system. Kaplan-Meier (K-M) survival curves and receiver operating characteristic (ROC) curves were utilized to validate the reliability of the risk score system. Finally, based on the correlations between DELs and DEGs involved in the risk score system, the final ceRNA network was identified. Results A total of 340 DELs, 29 DEMs, and 218 DEGs were selected to construct the initial ceRNA network. Functional enrichment analyses indicated that 218 DEGs were significantly enriched in the GO terms “nucleoplasm”, “transcription factor complex”, “protein binding”, and “metal ion binding”, whereas these DEGs were associated with the KEGG pathway terms “microRNAs in cancer”, “pathways in cancer”, “cell cycle”, “HTLV-1 infection”, and the “PI3K-Akt signalling pathway”. K-M survival analysis of all differentially expressed genes involved in the ceRNA network identified 24 DELs, 4 DEMs, and 29 DEGs, all of which were significantly correlated with LUAD progression (P < 0.05). Furthermore, 15 LUAD drugs interacting with 29 DEGs were selected. After lasso-penalized Cox regression and multivariate Cox regression modelling, 4 DEGs, PRKCE, DLC1, LATS2, and DPY19L1, were incorporated into the risk score system. The area under the curve (AUC) values of the time-dependent ROC curves at 3 years and 5 years were both higher than 0.5. Finally, the correlation coefficients between these 4 DEGs and their corresponding DELs involved in the ceRNA network suggested that there were 2 DEL-DEG pairs, NAV2-AS2 – PRKCE (r = 0.430, P < 0.001) and NAV2-AS2 – LATS2 (r = 0.338, P < 0.001). Considering the previously constructed ceRNA network, NAV2-AS2 – mir-31 – PRKCE and NAV2-SA2 – mir-31 – LATS2 were identified. Conclusions The lncRNA-miRNA-mRNA ceRNA network plays an essential role in LUAD. These results may improve our understanding and provide novel mechanistic insights to explore diagnostics, tumourigenesis, prognosis, and therapeutic drugs for LUAD patients.

Identification of Prognostic Factors Related to Super Enhancer-Regulated ceRNA Network in Metastatic Lung Adenocarcinoma

Identification of Key Genes and Evaluation of Clinical Outcomes in Lung Squamous Cell Carcinoma Using Integrated Bioinformatics Analysis

An integrated analysis of the competing endogenous RNA network associated of prognosis of stage I lung adenocarcinoma

Predictions of the dysregulated competing endogenous RNA signature involved in the progression of human lung adenocarcinoma.

Functional characterization and clinical significance of super-enhancers in lung adenocarcinoma

Comprehensive analysis of the putative competing endogenous RNA network in human lung adenocarcinoma

Identification of Cerna Network Based on a RNA-seq Shows Prognostic Lncrna Biomarkers in Human Lung Adenocarcinoma

Comprehensive analysis to identify a novel PTEN-associated ceRNA regulatory network as a prognostic biomarker for lung adenocarcinoma

Integrated Analysis of Long Non-Coding RNA-associated Cerna Network Reveals Potential Lncrna Biomarkers in Human Lung Adenocarcinoma

Identification of an Exosome-relevant SNHG6-hsa-miR-429- CHRDL1/CCNA2 Axis for Lung Adenocarcinoma Prognosis Evaluation

Identification of Three Lncrnas As Potential Predictive Biomarkers of Lung Adenocarcinoma

Identification of LINC02310 As an Enhancer in Lung Adenocarcinoma and Investigation of Its Regulatory Network Via Comprehensive Analyses

Identification of a ZC3H12D-regulated competing endogenous RNA network for prognosis of lung adenocarcinoma at single-cell level

Identification of Hub Driving Genes and Regulators of Lung Adenocarcinoma Based on the Gene Co-expression Network.

Identification of a Putative Competitive Endogenous RNA Network for Lung Adenocarcinoma Using TCGA Datasets

Identification of a Novel Cancer Stemness-Associated Cerna Axis in Lung Adenocarcinoma Via Stemness Indices Analysis.

Comprehensive Analysis of Aberrantly Expressed Profiles of Lncrnas and Mirnas with Associated Cerna Network in Lung Adenocarcinoma and Lung Squamous Cell Carcinoma

The Regulatory Network and Potential Role of LINC00973-miRNA-mRNA ceRNA in the Progression of Non-Small-Cell Lung Cancer

Identification key genes, key miRNAs and key transcription factors of lung adenocarcinoma

The super-enhancer-driven lncRNA LINC00880 acts as a scaffold between CDK1 and PRDX1 to sustain the malignance of lung adenocarcinoma

Prognostic Value of Lncrnas Related to Fatty Acid Metabolism in Lung Adenocarcinoma and Their Correlation with Tumor Microenvironment Based on Bioinformatics Analysis.