Identification of EGFR-Related LINC00460/mir-338-3p/MCM4 Regulatory Axis as Diagnostic and Prognostic Biomarker of Lung Adenocarcinoma Based on Comprehensive Bioinformatics Analysis and Experimental Validation

SIMPLE SUMMARY: While the epidermal growth factor receptor (EGFR) is an important target for lung adenocarcinoma (LUAD) therapy, acquired resistance is still inevitable. A comprehensive bioinformatics analysis strongly suggested that the closely related LINC00460-mir-338-3p-MCM4 ceRNA network of EGFR plays an important role in the diagnosis and prognosis of LUAD. High expression of LINC00460 and MCM4 predicts shorter patient survival. Univariate and multivariate regression analyses demonstrated that higher expression of LINC00460 and MCM4 was significantly associated with tumor size, lymph node metastasis, distant metastasis and TNM stage. A multi-gene regulation model analysis revealed that the LINC00460 (downregulation)—mir-338-3p (upregulation))—MCM4 (downregulation) pattern significantly improved the overall survival of LUAD patients. We also verified the expression of these genes in LUAD cell lines and tumor tissues by RT-PCR and immunohistochemistry. Finally, the possible targeted drugs of MCM4 were queried through the drug database platform, hoping to solve its drug resistance problem by targeting EGFR-related genes. ABSTRACT: Background: Lung adenocarcinoma (LUAD) is one of the most aggressive and lethal tumor types and requires effective diagnostic and therapeutic targets. Though the epidermal growth factor receptor (EGFR) is an important target for LUAD therapy, acquired resistance is still inevitable. In recent years, the regulation of the EGFR by competing endogenous RNAs (ceRNAs) has been extensively studied and significant progress has been made. Therefore, we aim to find new targets for the diagnosis and treatment of LUAD by analyzing the EGFR-related ceRNA network in LUAD and expect to address the problem of EGFR resistance. Methods: We identified differentially expressed lncRNAs, miRNAs and mRNAs closely associated with the EGFR by analyzing transcriptome data from LUAD samples. Comprehensive bioinformatics analysis strongly suggests that the LINC00460—mir-338-3p—MCM4 ceRNA network plays an important role in the diagnosis and prognosis of LUAD. The effects of different patterns of the LINC00460/MCM4 axis on the overall survival of patients with LUAD were analyzed by a polygene regulation model. We also verified the expression of these genes in LUAD cell lines and tumor tissues by RT-PCR and immunohistochemistry. The functional enrichment analysis and targeted drug prediction of the MCM4 gene were performed. Results: Survival analysis indicated that high expressions of LINC00460 and MCM4 predict a shorter survival period for patients. Univariate and multivariate regression analyses demonstrated that higher expressions of LINC00460 and MCM4 were significantly associated with tumor size, lymph node metastasis, distant metastasis and TNM stage. A multi-gene regulation model analysis revealed that the LINC00460 (downregulation)—mir-338-3p (upregulation)—MCM4 (downregulation) pattern significantly improved the overall survival of LUAD patients (p = 0.0093). RT-PCR and immunohistochemical experiments confirmed our analytical results. In addition, the functional enrichment analysis indicated that MCM4-related genes were mainly enriched in the cell cycle and cell division. A functional association network analysis showed that MCM4 was closely related to the EGFR. Finally, the possible targeted drugs of MCM4 were queried through the drug database platform, hoping to solve its drug resistance problem by targeting EGFR-related genes. Conclusions: In summary, the LINC00460/MCM4 axis can be used as a potential new perspective for targeting EGFR genes in precision medicine and is expected to serve as a diagnostic, prognostic and drug target for LUAD.