Co-Expression Network Analysis Unveiled lncRNA-mRNA Links Correlated to Epidermal Growth Factor Receptor-Tyrosine Kinase Inhibitor Resistance and/or Intermediate Epithelial-to-Mesenchymal Transition Phenotypes in a Human Non-Small Cell Lung Cancer Cellular Model System

We investigated mRNA-lncRNA co-expression patterns in a cellular model system of non-small cell lung cancer (NSCLC) sensitive and resistant to the epithelial growth factor receptor (EGFR) tyrosine kinase inhibitors (TKIs) erlotinib/gefitinib. The aim of this study was to unveil insights into the complex mechanisms of NSCLC targeted therapy resistance and epithelial-to-mesenchymal transition (EMT). Genome-wide RNA expression was quantified for weighted gene co-expression network analysis (WGCNA) to correlate the expression levels of mRNAs and lncRNAs. Functional enrichment analysis and identification of lncRNAs were conducted on modules associated with the EGFR-TKI response and/or intermediate EMT phenotypes. We constructed lncRNA-mRNA co-expression networks and identified key modules and their enriched biological functions. Processes enriched in the selected modules included RHO (A, B, C) GTPase and regulatory signaling pathways, apoptosis, inflammatory and interleukin signaling pathways, cell adhesion, cell migration, cell and extracellular matrix organization, metabolism, and lipid metabolism. Interestingly, several lncRNAs, already shown to be dysregulated in cancer, are connected to a small number of mRNAs, and several lncRNAs are interlinked with each other in the co-expression network.