Romi Biswas

Identifying novel targets of oncogenic EGF receptor signaling in lung cancer through global phosphoproteomics

Mutations in the epidermal growth factor receptor (EGFR) kinase domain occur in 10–30% of lung adenocarcinoma and are associated with tyrosine kinase inhibitor (TKI) sensitivity. We sought to identify the immediate direct and indirect phosphorylation targets of mutant EGFRs in lung adenocarcinoma. We undertook SILAC strategy, phosphopeptide enrichment, and quantitative MS to identify dynamic changes of phosphorylation downstream of mutant EGFRs in lung adenocarcinoma cells harboring EGFRL858R and EGFRL858R/T790M, the TKI‐sensitive, and TKI‐resistant mutations, respectively. Top canonical pathways that were inhibited upon erlotinib treatment in sensitive cells, but not in the resistant cells include EGFR, insulin receptor, hepatocyte growth factor, mitogen‐activated protein kinase, mechanistic target of rapamycin, ribosomal protein S6 kinase beta 1, and Janus kinase/signal transducer and activator of transcription signaling. We identified phosphosites in proteins of the autophagy network, such as ULK1 (S623) that is constitutively phosphorylated in these lung adenocarcinoma cells; phosphorylation is inhibited upon erlotinib treatment in sensitive cells, but not in resistant cells. Finally, kinase–substrate prediction analysis from our data indicated that substrates of basophilic kinases from, AGC and Calcium and calmodulin‐dependent kinase groups, as well as STE group kinases were significantly enriched and those of proline‐directed kinases from, CMGC and Casein kinase groups were significantly depleted among substrates that exhibited increased phosphorylation upon EGF stimulation and reduced phosphorylation upon TKI inhibition. This is the first study to date to examine global phosphorylation changes upon erlotinib treatment of lung adenocarcinoma cells and results from this study provide new insights into signaling downstream of mutant EGFRs in lung adenocarcinoma. All MS data have been deposited in the ProteomeXchange with identifier PXD001101 (http://proteomecentral.proteomexchange.org/dataset/PXD001101).

Quantitative tyrosine phosphoproteomics of Epidermal growth factor receptor (EGFR) tyrosine kinase inhibitor-treated lung adenocarcinoma cells reveals potential novel biomarkers of therapeutic response

Mutations in the Epidermal growth factor receptor (EGFR) kinase domain, such as the L858R missense mutation and deletions spanning the conserved sequence 747LREA750, are sensitive to tyrosine kinase inhibitors (TKIs). The gatekeeper site residue mutation, T790M accounts for around 60% of acquired resistance to EGFR TKIs. The first generation EGFR TKIs, erlotinib and gefitinib, and the second generation inhibitor, afatinib are FDA approved for initial treatment of EGFR mutated lung adenocarcinoma. The predominant biomarker of EGFR TKI responsiveness is the presence of EGFR TKI-sensitizing mutations. However, 30–40% of patients with EGFR mutations exhibit primary resistance to these TKIs, underscoring the unmet need of identifying additional biomarkers of treatment response. Here, we sought to characterize the dynamics of tyrosine phosphorylation upon EGFR TKI treatment of mutant EGFR-driven human lung adenocarcinoma cell lines with varying sensitivity to EGFR TKIs, erlotinib and afatinib. We employed stable isotope labeling with amino acids in cell culture (SILAC)-based quantitative mass spectrometry to identify and quantify tyrosine phosphorylated peptides. The proportion of tyrosine phosphorylated sites that had reduced phosphorylation upon erlotinib or afatinib treatment correlated with the degree of TKI-sensitivity. Afatinib, an irreversible EGFR TKI, more effectively inhibited tyrosine phosphorylation of a majority of the substrates. The phosphosites with phosphorylation SILAC ratios that correlated with the TKI-sensitivity of the cell lines include sites on kinases, such as EGFR-Y1197 and MAPK7-Y221, and adaptor proteins, such as SHC1-Y349/350, ERRFI1-Y394, GAB1-Y689, STAT5A-Y694, DLG3-Y705, and DAPP1-Y139, suggesting these are potential biomarkers of TKI sensitivity. DAPP1, is a novel target of mutant EGFR signaling and Y-139 is the major site of DAPP1 tyrosine phosphorylation. We also uncovered several off-target effects of these TKIs, such as MST1R-Y1238/Y1239 and MET-Y1252/1253. This study provides unique insight into the TKI-mediated modulation of mutant EGFR signaling, which can be applied to the development of biomarkers of EGFR TKI response.

Genomic profiling of multiple sequentially acquired tumor metastatic sites from an “exceptional responder” lung adenocarcinoma patient reveals extensive genomic heterogeneity and novel somatic variants driving treatment response

We used next-generation sequencing to identify somatic alterations in multiple metastatic sites from an “exceptional responder” lung adenocarcinoma patient during his 7-yr course of ERBB2-directed therapies. The degree of heterogeneity was unprecedented, with ∼1% similarity between somatic alterations of the lung and lymph nodes. One novel translocation, PLAG1-ACTA2, present in both sites, up-regulated ACTA2 expression. ERBB2, the predominant driver oncogene, was amplified in both sites, more pronounced in the lung, and harbored an L869R mutation in the lymph node. Functional studies showed increased proliferation, migration, metastasis, and resistance to ERBB2-directed therapy because of L869R mutation and increased migration because of ACTA2 overexpression. Within the lung, a nonfunctional CDK12, due to a novel G879V mutation, correlated with down-regulation of DNA damage response genes, causing genomic instability, and sensitivity to chemotherapy. We propose a model whereby a subclone metastasized early from the primary site and evolved independently in lymph nodes.

Quantitative mass spectrometry to interrogate proteomic heterogeneity in metastatic lung adenocarcinoma and validate a novel somatic mutation CDK12-G879V

Lung cancer is the leading cause of cancer death both in men and women. Tumor heterogeneity is an impediment to targeted treatment of all cancers, including lung cancer. Here, we sought to characterize changes in tumor proteome and phosphoproteome by longitudinal, prospective collection of tumor tissue of an exceptional responder lung adenocarcinoma patient who survived with metastatic lung adenocarcinoma for more than seven years with HER2-directed therapy in combination with chemotherapy. We employed “Super-SILAC” and TMT labeling strategies to quantify the proteome and phosphoproteome of a lung metastatic site and ten different metastatic progressive lymph nodes collected across a span of seven years, including five lymph nodes procured at autopsy. We identified specific signaling networks enriched in lung compared to the lymph node metastatic sites. We correlated the changes in protein abundance with changes in copy number alteration (CNA) and transcript expression. To further interrogate the mass spectrometry data, patient-specific database was built incorporating all the somatic variants identified by whole genome sequencing (WGS) of genomic DNA from the lung, one lymph node metastatic site and blood. An extensive validation pipeline was built for confirmation of variant peptides. We validated 360 spectra corresponding to 55 germline and 6 somatic variant peptides. Targeted MRM assays demonstrated expression of two novel variant somatic peptides, CDK12-G879V and FASN-R1439Q, with expression in lung and lymph node metastatic sites, respectively. CDK12 G879V mutation likely results in a nonfunctional CDK12 kinase and chemotherapy susceptibility in lung metastatic sites. Knockdown of CDK12 in lung adenocarcinoma cells results in increased chemotherapy sensitivity, explaining the complete resolution of the lung metastatic sites in this patient.

Integrated proteogenomic analysis of metastatic thoracic tumors identifies APOBEC mutagenesis and copy number alterations as drivers of proteogenomic tumor evolution and heterogeneity

Elucidation of the proteogenomic evolution of metastatic tumors may offer insight into the poor prognosis of patients harboring metastatic disease. We performed whole-exome and transcriptome sequencing, copy number alterations (CNA) and mass spectrometry-based quantitative proteomics of 37 lung adenocarcinoma (LUAD) and thymic carcinoma (TC) metastases obtained by rapid autopsy and found evidence of patient-specific, multi-dimensional heterogeneity. Extreme mutational heterogeneity was evident in a subset of patients whose tumors showed increased APOBEC-signature mutations and expression of APOBEC3 region transcripts compared to patients with lesser mutational heterogeneity. TP53 mutation status was associated with APOBEC hypermutators in our cohort and in three independent LUAD datasets. In a thymic carcinoma patient, extreme heterogeneity and increased APOBEC3AB expression was associated with a high-risk germline APOBEC3AB variant allele. Patients with CNA occurring late in tumor evolution had corresponding changes in gene expression and protein abundance indicating genomic instability as a mechanism of downstream transcriptomic and proteomic heterogeneity between metastases. Across all tumors, proteomic heterogeneity was greater than copy number and transcriptomic heterogeneity. Enrichment of interferon pathways was evident both in the transcriptome and proteome of the tumors enriched for APOBEC mutagenesis despite a heterogeneous immune microenvironment across metastases suggesting a role for the immune microenvironment in the expression of APOBEC transcripts and generation of mutational heterogeneity. The evolving, heterogeneous nature of LUAD and TC, through APOBEC-mutagenesis and CNA illustrate the challenges facing treatment outcomes.

Abstract 4759: Whole genome sequencing of sequentially acquired lung and lymph node metastatic sites from a never smoker lung adenocarcinoma patient revealed extensive genomic heterogeneity

Proceedings: AACR 106th Annual Meeting 2015; April 18-22, 2015; Philadelphia, PA Recent large-scale cancer genome sequencing studies have uncovered extensive diverse mutational landscapes in lung cancer patients. Furthermore, tumor heterogeneity has been widely recognized and it has significant clinical implications in selection of targeted treatment strategies as well as treatment response. Using whole genome sequencing, we demonstrate an unprecedented genomic heterogeneity between sequentially acquired lung and lymph node metastatic sites from an African American never-smoker lung adenocarcinoma patient who has survived with metastatic disease for over seven years while being treated with single or combination HER2-directed therapies. We determined that less than 1% of somatic variants were common between the two tumor sites. Copy number variations were more intense in the lung tumor than in the metastatic lymph node. We identified several novel somatic mutations in key cancer genes in both sites. Interestingly, one novel translocation, PLAG1-ACTA2 was highly expressed in both the lung and lymph node metastases resulting in overexpression of ACTA2, which has been suggested to increase the metastatic potential in lung adenocarcinoma. Using ultra deep targeted re-sequencing, we validated all non-synonymous variants, and approximately 80% of those identified in the metastatic lymph node were also present in a second lymph node biopsied two years after the first one. Although this degree of tumor heterogeneity was surprising, somatic variants affected key hallmarks of tumorigenesis in both sites. These findings suggest a model of early metastatic spread and parallel clonal evolution in disparate metastatic sites. Citation Format: Shaojian Gao, Constance Cultraro, Romi BIswas, Corey A. Carter, Tapan K. Maity, Anish Thomas, Arun Rajan, Paul Meltzer, David Schrump, Giuseppe Giaccone, Javed Khan, Udayan Guha. Whole genome sequencing of sequentially acquired lung and lymph node metastatic sites from a never smoker lung adenocarcinoma patient revealed extensive genomic heterogeneity. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 4759. doi:10.1158/1538-7445.AM2015-4759

Abstract 5260: Identification of substrates of lung cancer-specific mutant EGFR kinases.

Proceedings: AACR 104th Annual Meeting 2013; Apr 6-10, 2013; Washington, DC Among various cancers, lung cancer is the leading cause of cancer-related mortality in both men and women in the US. The major driver mutations in lung cancer are in KRAS and EGFR. Mutations in EGFR account for 10-30% of lung adenocarcinoma. Somatic activating mutations in the kinase domain of EGFR occur in a subset of patients with lung adenocarcinoma and are associated with sensitivity to EGFR-directed tyrosine kinase inhibitors (TKIs). Unfortunately, all patients who respond dramatically to TKIs (gefitinib and erlotinib) eventually develop secondary resistance. A mutation in the gatekeeper residue in the ATP binding pocket of EGFR (T790M) occurs in around 50% of such patients. We hypothesized that mutant EGFR kinases have altered substrate specificities in vivo; these substrates play intermediary roles in transducing mutant EGFR signaling to downstream components to promote tumorigenesis. We have employed a screening method using an in vitro kinase assay on protein arrays containing around 17,000 proteins to compare the substrate specificity of the kinase domains of wild type EGFR, L858R EGFR (the most common activating mutant), and the EGFR harboring both the L858R and T790M mutations (TKI resistant mutant). We have identified several potential substrates such as Crk, Crk-like protein, Grb2, Stat3, Nck1, PIK3R3, etc., which have altered phosphorylation status in response to mutant EGFR kinase activity. We also performed a SILAC (Stable isotope labeling by amino acids in cell culture)-based global phosphoproteomic analysis in human lung adenocarcinoma cells using mass spectrometry and observed that phosphorylation of some of the potential substrates is modulated by EGF stimulation and TKI inhibition. We are currently validating promising substrates from the above assay by performing in vitro kinase assays in presence or absence of TKIs. Furthermore, to identify specific phosphorylation sites of the validated targets in vivo, we have adapted a SILAC-based mass spectrometry approach to immunoprecipitate specific substrates from isogenic human bronchial epithelial cells expressing various EGFR mutants followed by identification and quantification of phosphorylation. We have used siRNA-mediated gene silencing to interrogate the role of several of these substrates in survival of human lung adenocarcinoma cells. Finally, we are examining the phosphorylation status of identified substrates in tumor samples from patients of lung adenocarcinoma to correlate phosphorylation with EGFR mutation status. Citation Format: Tapan K. Maity, Hee-Sool Rho, Natalya Belinka, Zhi Xie, Xu Zhang, Abhilash Venugopalan, Romi Biswas, Constance Cultraro, Heng Zhu, Udayan Guha. Identification of substrates of lung cancer-specific mutant EGFR kinases. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 5260. doi:10.1158/1538-7445.AM2013-5260