Purpose EMT has been associated with metastatic spread and EGFR inhibitor resistance. We developed and validated a robust 76-gene EMT signature using gene expression profiles from four platforms using NSCLC cell lines and patients treated in the BATTLE study. Methods We conducted an integrated gene expression, proteomic, and drug response analysis using cell lines and tumors from NSCLC patients. A 76-gene EMT signature was developed and validated using gene expression profiles from four microarray platforms of NSCLC cell lines and patients treated in the BATTLE (Biomarker-integrated Approaches of Targeted Therapy for Lung Cancer Elimination) study, and potential therapeutic targets associated with EMT were identified. Results Compared with epithelial cells, mesenchymal cells demonstrated significantly greater resistance to EGFR and PI3K/Akt pathway inhibitors, independent of EGFR mutation status, but more sensitivity to certain chemotherapies. Mesenchymal cells also expressed increased levels of the receptor tyrosine kinase Axl and showed a trend towards greater sensitivity to the Axl inhibitor SGI-7079, while the combination of SGI-7079 with erlotinib reversed erlotinib resistance in mesenchymal lines expressing Axl and in a xenograft model of mesenchymal NSCLC. In NSCLC patients, the EMT signature predicted 8-week disease control in patients receiving erlotinib, but not other therapies. Conclusion We have developed a robust EMT signature that predicts resistance to EGFR and PI3K/Akt inhibitors, highlights different patterns of drug responsiveness for epithelial and mesenchymal cells, and identifies Axl as a potential therapeutic target for overcoming EGFR inhibitor resistance associated with the mesenchymal phenotype

The molecular underpinnings that drive the heterogeneity of KRAS-mutant lung adenocarcinoma (LUAC) are poorly characterized. We performed an integrative analysis of genomic, transcriptomic and proteomic data from early-stage and chemo-refractory LUAC and identified three robust subsets of KRAS-mutant LUAC dominated, respectively, by co-occurring genetic events in STK11/LKB1 (the KL subgroup), TP53 (KP) and CDKN2A/B inactivation coupled with low expression of the NKX2-1 (TTF1) transcription factor (KC). We further reveal biologically and therapeutically relevant differences between the subgroups. KC tumors frequently exhibited mucinous histology and suppressed mTORC1 signaling. KL tumors had high rates of KEAP1 mutational inactivation and expressed lower levels of immune markers, including PD-L1. KP tumors demonstrated higher levels of somatic mutations, inflammatory markers, immune checkpoint effector molecules and improved relapse-free survival. Differences in drug sensitivity patterns were also observed; notably, KL cells showed increased vulnerability to HSP90-inhibitor therapy. This work provides evidence that co-occurring genomic alterations identify subgroups of KRAS-mutant LUAC with distinct biology and therapeutic vulnerabilities.

Systematic efforts are underway to decipher the genetic changes associated with tumor initiation and progression. However, widespread clinical application of this information is hampered by an inability to identify critical genetic events across the spectrum of human tumors with adequate sensitivity and scalability. Here, we have adapted high-throughput genotyping to query 238 known oncogene mutations across 1,000 human tumor samples. This approach established robust mutation distributions spanning 17 cancer types. Of 17 oncogenes analyzed, we found 14 to be mutated at least once, and 298 (30%) samples carried at least one mutation. Moreover, we identified previously unrecognized oncogene mutations in several tumor types and observed an unexpectedly high number of co-occurring mutations. These results offer a new dimension in tumor genetics, where mutations involving multiple cancer genes may be interrogated simultaneously and in 'real time' to guide cancer classification and rational therapeutic intervention.