Pierobon, Mariaelena

Of the thousands of genes and substrates identified in KRAS-mutant signaling networks in pancreatic ductal adenocarcinoma (PDAC), more than 200 are transcription factors, implying extensive and complex transcriptional regulation. However, we observed that genetic suppression of the transcription factor MYC alone was sufficient to phenocopy the effect of KRAS suppression in signaling, growth, and metabolic processes in PDAC cells. We determined the gene transcription changes caused by acute suppression of MYC function in KRAS -mutant PDAC cell lines and performed dependency map and pathway analyses on the affected gene sets. The expression of 1685 genes was increased upon suppression of MYC, and this gene set may comprise the bulk of the MYC-regulated genes essential for PDAC growth. In contrast, the 1325 genes whose expression was inhibited may comprise a compensatory response to oncogenic stress, mediated in part by the GTPase RHO. MYC-dependent transcriptional activity was largely ERK dependent, and almost one-third of ERK-regulated genes were also regulated by MYC in PDAC cells. Furthermore, chemical proteomic profiling revealed MYC-regulated protein kinases that can be targeted therapeutically. Together, these data provide a molecular portrait of MYC-dependent signaling that encompasses potentially exploitable mechanisms for treating PDAC.

This prospective phase II clinical trial (Side Out 2) explored the clinical benefits of treatment selection informed by multi‐omic molecular profiling (MoMP) in refractory metastatic breast cancers (MBCs). Core needle biopsies were collected from 32 patients with MBC at trial enrollment. Patients had received an average of 3.94 previous lines of treatment in the metastatic setting before enrollment in this study. Samples underwent MoMP, including exome sequencing, RNA sequencing (RNA‐Seq), immunohistochemistry, and quantitative protein pathway activation mapping by Reverse Phase Protein Microarray (RPPA). Clinical benefit was assessed using the previously published growth modulation index (GMI) under the hypothesis that MoMP‐selected therapy would warrant further investigation for GMI ≥ 1.3 in ≥ 35% of the patients. Of the 32 patients enrolled, 29 received treatment based on their MoMP and 25 met the follow‐up criteria established by the trial protocol. Molecular information was delivered to the tumor board in a median time frame of 14 days (11–22 days), and targetable alterations for commercially available agents were found in 23/25 patients (92%). Of the 25 patients, 14 (56%) reached GMI ≥ 1.3. A high level of DNA topoisomerase I (TOPO1) led to the selection of irinotecan‐based treatments in 48% (12/25) of the patients. A pooled analysis suggested clinical benefit in patients with high TOPO1 expression receiving irinotecan‐based regimens (GMI ≥ 1.3 in 66.7% of cases). These results confirmed previous observations that MoMP increases the frequency of identifiable actionable alterations (92% of patients). The MoMP proposed allows the identification of biomarkers that are frequently expressed in MBCs and the evaluation of their role as predictors of response to commercially available agents. Lastly, this study confirmed the role of MoMP for informing treatment selection in refractory MBC patients: more than half of the enrolled patients reached a GMI ≥ 1.3 even after multiple lines of previous therapies for metastatic disease.

Abstract Pancreatic ductal adenocarcinoma (PDAC) is characterized by KRAS- and autophagy-dependent tumorigenic growth. We determined the role of mutationally activated KRAS, found in ~95% of PDAC, in supporting autophagy. Surprisingly, acute KRAS suppression, which blocks proliferation, was associated with increased rather than decreased autophagic flux. Pharmacologic inhibition of ERK MAPK phenocopied the genetic silencing of KRAS and also increased autophagic flux. We speculated that the loss of ERK-driven metabolic processes may induce compensatory mechanisms to increase autophagy. Addressing a mechanism for ERK suppression-increased autophagy, we describe three mechanisms: 1) ERK inhibition-mediated increased transcription of autophagy and lysosomal genes, 2) ERK inhibitor-induced AMPK activation and suppression of mTOR signaling, and 3) ERK inhibition-facilitated decreased glycolytic flux. We then addressed whether ERK inhibition increased PDAC dependence on autophagy. Supporting this possibility, we found that cotreatment with the autophagy inhibitor chloroquine (CQ) synergistically enhanced ERK inhibitor-mediated antiproliferative activity. Similarly, genetic or pharmacologic inhibition of specific regulators of autophagy also enhanced ERK inhibitor activity. Encouraged by the synergistic relationship between ERK and autophagy inhibition, we performed a CRISPR/Cas-9 mediated genetic loss of function screen in the presence of CQ to determine additional sensitizers as well as mediators of resistance to autophagy inhibition. Top sensitizers included multiple mediators of the DNA damage response. One such sensitizing gene was CHEK1, which encodes the CHK1 serine/threonine kinase, required for checkpoint-mediated cell cycle arrest and activation of DNA repair in response to DNA damage. Accordingly, we showed that treatment with prexasertib, a clinical candidate CHK1 inhibitor, increased autophagic flux in PDAC cells and synergized with CQ to decrease PDAC cell proliferation and increase apoptosis. We conclude that concurrent suppression of multiple metabolic processes, to block compensatory rebound activities, will be needed for effective PDAC treatment. This abstract is also being presented as Poster B07. Citation Format: Kirsten L. Bryant, Jennifer E. Klomp, Ye S. Lee, Clint A. Stalnecker, Kajal R. Grover, A. Cole Edwards, Sen Peng, Mariaelena Pierobon, Emanuel F. Petricoin III, Nhan Tran, Alec C. Kimmelman, Adrienne D. Cox, Channing J. Der. Enhancing the effect of autophagy inhibition for pancreatic cancer treatment [abstract]. In: Proceedings of the AACR Special Conference on Pancreatic Cancer: Advances in Science and Clinical Care; 2019 Sept 6-9; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2019;79(24 Suppl):Abstract nr PR10.

Abstract Background The dearth of relevant tumor models reflecting the heterogeneity of human central nervous system metastasis (CM) has hindered development of novel therapies. Methods We established 39 CM patient-derived xenograft (PDX) models representing the histological spectrum, and performed phenotypic and multi-omic characterization of PDXs and their original patient tumors. PDX clonal evolution was also reconstructed using allele-specific copy number and somatic variants. Results PDXs retained their metastatic potential, with flank-implanted PDXs forming spontaneous metastases in multiple organs, including brain, and CM subsequent to intracardiac injection. PDXs also retained the histological and molecular profiles of the original patient tumors, including retention of genomic aberrations and signaling pathways. Novel modes of clonal evolution involving rapid expansion by a minor clone were identified in 2 PDXs, including CM13, which was highly aggressive in vivo forming multiple spontaneous metastases, including to brain. These PDXs had little molecular resemblance to the patient donor tumor, including reversion to a copy number neutral genome, no shared nonsynonymous mutations, and no correlation by gene expression. Conclusions We generated a diverse and novel repertoire of PDXs that provides a new set of tools to enhance our knowledge of CM biology and improve preclinical testing. Furthermore, our study suggests that minor clone succession may confer tumor aggressiveness and potentiate brain metastasis.

Abstract Despite recent advances in cancer therapy, options for central nervous system metastasis (CM) remain limited and patient outcomes are dismal. The lack of tumor models that reflect the heterogeneity and complexity of human CM tumors is a barrier to understanding the underlying biology and to developing novel therapies. Here, we report the development and characterization of 39 patient derived xenograft (PDX) models of CM tumors with histological subtypes representing the disease spectrum. The PDX models successfully grew both in the flank and brains of mice, and retained their histology regardless of the site of implantation. Brain metastasis also formed in 5/12 models after intra-cardiac injection. Multi-omic analysis of the PDX tumor and the patient-matched CM tumor showed that PDXs largely retained the molecular profiles of the tumor from which they were derived, including the retention of aberrations in key driver genes and signaling pathways. PDXs differed primarily from the ungrafted patient tumors by a downregulation in immunity-related pathways, which was expected due to the use of an immuno-compromised host, arguing for the validity of using the flank as a site of implantation. Integrated analysis of mutations and gene expression changes in CM revealed an upregulation of neuronal signaling and DNA damage response, and downregulation of calcium signaling and inflammation. We also assessed the clonal evolution of the PDX from the original patient CM by tracing variant allele frequencies in copy number neutral regions. We found that the PDX tumors were heterogeneous and comprised of multiple tumor clones that were either unique or in common with the original tumor, and displayed evidence of clonal expansion. Two PDXs that underwent phylogeny assessment showed evidence of minor clone expansion. One such PDX was a tumor from a Her2+ breast cancer (CM13), which largely did not molecularly resemble its matched patient tumor but was confirmed to have some shared expressed variants. The PDX was nearly copy number neutral, which included loss of the Her2 amplicon, and had few mutations compared with its parent tumor. Most noteworthy, is that CM13 displayed unique molecular patterns, aggressive growth and was highly metastatic in vivo. From the flank, it led to multiple metastases, including micrometastasis in the brain. This suggests that minor clone takeover may confer tumor aggressiveness and potentiate brain metastasis. The large and diverse repertoire of PDXs developed and characterized in this study provides a new set of tools that will help deepen our understanding of CM and improve preclinical testing for CM therapies. Our data also show that expansion of minor clones could confer aggressive tumor behavior and may be important in tumor evolution in PDX models and human cancers. Citation Format: Ben Yi Tew, Christophe Legendre, Tim Triche, Gerald C. Gooden, Kyle N. Johnson, Rae Anne Martinez, Emanuel F. Petricoin, Mariaelena Pierobon, Joyce O'Shaughnessy, Cindy Osborne, Mark A. Schroeder, Daniel J. Ma, Mark Bernstein, Jann N. Sarkaria, Steven A. Toms, Bodour Salhia. Development and characterization of patient derived xenografts from central nervous system metastasis reveal minor clone expansion linked with aggressive tumor behavior [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 991.

Abstract Pancreatic ductal adenocarcinoma (PDAC) cell growth is dependent on altered biochemical and cellular activities that facilitate increased metabolic dependencies. We determined the role of mutationally activated KRAS, found in ~95% of PDAC, in supporting autophagy. Surprisingly, acute KRAS suppression, which blocks proliferation, was associated with increased rather than decreased autophagic flux. Similarly, KRAS transformation of immortalized human pancreatic epithelial cells showed decreased autophagic flux. Pharmacologic inhibition of ERK MAPK phenocopied the genetic silencing of KRAS and also increased autophagic flux. Addressing a mechanism for ERK suppression increased autophagy, we found increased AMPK activity, decreased mTOR activity, and increased expression of autophagic proteins. We speculated that increased autophagy may be a compensatory mechanism caused by the loss of KRAS- and ERK-dependent metabolic processes. To address this possibility, we performed RNA-seq analyses to monitor gene expression changes in response to ERK inhibition. We observed increased transcription of genes responsible for the metabolic processes of autophagy and β-oxidation, with concurrent decreased transcription of glycolysis and mitochondrial-associated genes. These changes were validated in metabolomic analyses following either genetic silencing of KRAS or ERK inhibition. Thus, the loss of ERK-driven metabolic processes may induce compensatory mechanisms to increase autophagy. We then addressed whether ERK inhibition increased PDAC dependence on autophagy. Supporting this possibility, we found that cotreatment with the autophagy inhibitor chloroquine synergistically enhanced ERK inhibitor mediated antiproliferative activity. Similarly, genetic or pharmacologic inhibition of specific regulators of autophagy also enhanced ERK inhibitor activity. We conclude that concurrent suppression of multiple metabolic processes, to block compensatory rebound activities, will be needed for effective PDAC treatment. Citation Format: Kirsten Bryant, Sen Peng, Andrey Tikunov, Mariaelena Pierobon, Venugopal Gunda, Garima Tomar, Pankaj Singh, Emanuel Petricoin, Jeffrey Macdonald, Nhan Tran, Alec Kimmelman, Channing Der. Inhibition of ERK MAPK signaling increases pancreatic cancer dependency on autophagy [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 5492.

Abstract Background: Although the efficacy of different immunotherapeutic agents varies greatly across tumor types, little is known about the role of different metastatic tumor microenvironments in modulating the expression of immune-checkpoints like PD-L1. Using metastatic breast cancer as a model, this work explored whether metastatic lesions derived from different patients, but colonizing the same host target organ, showed organ-specific immune signatures. Methods: Snap frozen metastatic lesions collected from 30 breast cancer patients enrolled in a prospective phase II multi-omic precision medicine clinical trial (“Side-Out II”) were used for this analysis. Sites of metastasis were: liver (LI) (n=10), skin/chest wall (SC) (n=10), and other lesions (OL) (n=10) [including lymph nodes (LN) (n=4), lung (LU) (n=3), and intra-abdominal lesions (IA) (n=3)]. Pure tumor epithelia were isolated from each specimen via Laser Capture Microdissection (LCM). PD-L1 expression (Antibody Clone E1L3N) was quantitatively measured using the Reverse Phase Protein Microarray (RPPA), a high-throughput, multiplex immunoassay that provides operator-independent quantitative data starting from a small amount of biological material. Proportion of patients with PD-L1 expression above the median and the 75th percentile of the population was compared across metastatic sites. Results: PD-L1 expression had a >5 fold dynamic range across breast cancer-derived metastatic lesions with different distribution across metastatic sites. Specifically, PD-L1 expression of 6/10 LIs, 5/10 OLs (including 2 LU and 3 LN positive lesions), and 4/10 SC metastases was greater than the population median of all metastases combined. LI metastases were the most represented group above the 75th percentile of the population (4/10), followed by OS (2/10) (including 1 LU and 1 LN positive lesions), and SC (1/10). Of the 15 lesions above the population median, 6 were LI (40%), 5 were OS (33%), and 4 were SC (27%); of the 7 metastases above the 75th percentile, 4 were LI (57%), 2 were OS (29%), and 1 was SC (14%). Conclusions: The LCM-RPPA workflow can capture PD-L1 protein expression on a continuous quantitative scale and provide a broader understanding of the dynamic range of expression. Our data suggest that organ specific microenvironments in which metastatic lesions develop may strongly influence overall PD-L1 tumor cell expression. If different metastatic host organs can modulate PD-L1 expression, location of the metastatic lesion should be kept in consideration when selecting patients that are candidates for immunotherapy. Such hypotheses should be further validated in prospective clinical trial where quantitative PD-L1 expression of different metastatic lesions is evaluated along with response to targeted immuno-oncology compounds. Citation Format: Mariaelena Pierobon, K Alex Hodge, Elisa Baldelli, Steven Anthony, Nicholas J. Robert, Donald W. Northfelt, Mohammad J. Jahanzeb, Linda Vocila, Lance A. Liotta, Bryant Dunetz, Emanuel F. Petricoin. PD-L1 expression differs across cancer metastatic sites from breast tumors [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 3129.

Abstract Background: The efficacy of immunotherapy, including therapeutic strategies capable of modulating innate/adaptive immune resistance, varies greatly across tumor types. As the number of available immunotherapies accelerates, the study of predictive markers by IHC (e.g. PDL1 expression) is under intense investigation. However, staining protocol inconsistency, variation across scoring systems, and subjective interpretation of the immunostaining have produced conflicting results thus far. This work explored the role of Laser Capture Microdissection (LCM) coupled with Reverse Phase Protein Microarray (RPPA) as an alternative high throughput, quantitative, operator independent platform for measuring PDL1 expression across tumor types. Material and Methods: Pure epithelial cells were isolated via LCM from 178 samples including: 72 ovarian cancers (OC), 57 lung adenocarcinomas (LC), 30 metastatic breast cancers (MBC), and 19 pancreatic cancers (PC). PDL1 expression was measured on a continuous scale using quantitative RPPA based analysis. Each tumor type was processed and arrayed independently. Experimental samples and reference standards used for inter-assay normalization were printed in triplicates. Results: PDL1 expression varied greatly across tumor types. LC were characterized by the greatest intra-tumor fold dynamic range (> 35-fold), followed by OC (< 13-fold), MBC (< 4-fold), and PC (< 2-fold). PDL1 expression of 46/57 (80.7%) LC, 17/30 (56.7%) MBC, 6/19 (31.5%) PC, and 20/72 (27.8%) OC was greater than the population median of all tumors combined. Within the LC samples with PDL1 expression equal to the top quartile of the population, 10 (71.4%) were KRAS mutant lesions and 4 (28.6%) were WT tumors. Finally, amongst LC and PC harboring a KRAS mutation, PC showed an overall lower expression of PDL1 with only 2/19 (10.5%) cases been above the population median and none within the top quartile of the population. Conclusions: The LCM-RPPA workflow has the unique ability to capture immune checkpoint expression on a continuous quantitative scale as well as capture its broad dynamic range. Because RPPA is unconstrained by antigen retrieval issues as well as subjectivity of IHC interpretation, this approach may generate more accurate cut-point of therapeutic response prediction. Overall the dynamic range of PDL1 was broader in LC compared to other solid tumors, and LC had a much higher proportion of patients with tumors expressing high levels of PDL1. These quantitative differences may explain therapeutic efficacy of PDL1 inhibition across tumor types. Such speculative hypothesis should be further validated in prospective clinical trials. Finally, these preliminary data suggest that organ specific microenvironments more than specific driving mutations (e.g. KRAS) may strongly influence PDL1 expression in malignant lesions. Citation Format: Elisa Baldelli, Valerie Calvert, K. Alex Hodge, Maria Isabella Sereni, Guido Gambara, Eric B. Haura, Lucio Crino', Bryant Dunetz, Sergio Pecorelli, David J. Perry, Stephen P. Anthony, Nicholas Robert, Donald W. Northfelt, Mohammad Jahanzeb, Emanuel F. Petricoin, Mariaelena Pierobon. Quantitative measurement of PDL1 expression across tumor types using laser capture microdissection and reverse phase protein microarray [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 5656. doi:10.1158/1538-7445.AM2017-5656