Clinical, pathological, genetics and intratumoural immune milieu of micropapillary carcinoma of the colon

AIM

Micropapillary carcinoma (MPC) is a recognised WHO variant of colonic carcinoma (CC), although little is known about its prognosis, immune microenvironment and molecular alterations. We investigated its clinical, pathological and immunological characteristics.

METHODS

We assessed 903 consecutive CCs and used the WHO definition to identify MPC. We recorded serrated and mucinous differentiation and mismatch repair (MMR) status. We performed immunohistochemistry and quantification on tissue microarrays for HLA class I/II proteins, beta-2-microglobulin (B2MG), CD8, CD163, LAG3, PD-L1, FoxP3, PD-L1and BRAF V600E.

RESULTS

We classified 8.6% (N=78) of CC as MPC. Relative to non-MPC, MPC was more often high grade (p=0.03) and showed serrated morphology (p<0.01); however, we found no association with extramural venous invasion (p=0.41) and American Joint Committee on Cancer stage (p=0.95). MPCs showed lower numbers of CD8 positive lymphocytes (p<0.01), lower tumour cell B2MG expression (p=0.04) and lower tumour cell PD-L1 expression (p<0.01). There was no difference in HLA class I/II, LAG3, FOXP3, CD163 and PD-L1 positive histiocytes. There was no association with MMR status or BRAF V600E relative to non-MPC. MPC was not associated with decreased disease-specific survival (p=0.36).

CONCLUSION

MPCs are associated with high-grade differentiation and a less active immune microenvironment than non-MPC. MPC is not associated with inferior disease-specific survival.

Transfer Learning Reveals Cancer-Associated Fibroblasts Are Associated with Epithelial-Mesenchymal Transition and Inflammation in Cancer Cells in Pancreatic Ductal Adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is an aggressive malignancy characterized by an immunosuppressive tumor microenvironment enriched with cancer-associated fibroblasts (CAF). This study used a convergence approach to identify tumor cell and CAF interactions through the integration of single-cell data from human tumors with human organoid coculture experiments. Analysis of a comprehensive atlas of PDAC single-cell RNA sequencing data indicated that CAF density is associated with increased inflammation and epithelial-mesenchymal transition (EMT) in epithelial cells. Transfer learning using transcriptional data from patient-derived organoid and CAF cocultures provided in silico validation of CAF induction of inflammatory and EMT epithelial cell states. Further experimental validation in cocultures demonstrated integrin beta 1 (ITGB1) and vascular endothelial factor A (VEGFA) interactions with neuropilin-1 mediating CAF-epithelial cell cross-talk. Together, this study introduces transfer learning from human single-cell data to organoid coculture analyses for experimental validation of discoveries of cell-cell cross-talk and identifies fibroblast-mediated regulation of EMT and inflammation.

SIGNIFICANCE

Adaptation of transfer learning to relate human single-cell RNA sequencing data to organoid-CAF cocultures facilitates discovery of human pancreatic cancer intercellular interactions and uncovers cross-talk between CAFs and tumor cells through VEGFA and ITGB1.

Mesothelin CAR T Cells Secreting Anti-FAP/Anti-CD3 Molecules Efficiently Target Pancreatic Adenocarcinoma and its Stroma

PURPOSE

Targeting solid tumors with chimeric antigen receptor (CAR) T cells remains challenging due to heterogenous target antigen expression, antigen escape, and the immunosuppressive tumor microenvironment (TME). Pancreatic cancer is characterized by a thick stroma generated by cancer-associated fibroblasts (CAF), which may contribute to the limited efficacy of mesothelin-directed CAR T cells in early-phase clinical trials. To provide a more favorable TME for CAR T cells to target pancreatic ductal adenocarcinoma (PDAC), we generated T cells with an antimesothelin CAR and a secreted T-cell-engaging molecule (TEAM) that targets CAF through fibroblast activation protein (FAP) and engages T cells through CD3 (termed mesoFAP CAR-TEAM cells).

EXPERIMENTAL DESIGN

Using a suite of in vitro, in vivo, and ex vivo patient-derived models containing cancer cells and CAF, we examined the ability of mesoFAP CAR-TEAM cells to target PDAC cells and CAF within the TME. We developed and used patient-derived ex vivo models, including patient-derived organoids with patient-matched CAF and patient-derived organotypic tumor spheroids.

RESULTS

We demonstrated specific and significant binding of the TEAM to its respective antigens (CD3 and FAP) when released from mesothelin-targeting CAR T cells, leading to T-cell activation and cytotoxicity of the target cell. MesoFAP CAR-TEAM cells were superior in eliminating PDAC and CAF compared with T cells engineered to target either antigen alone in our ex vivo patient-derived models and in mouse models of PDAC with primary or metastatic liver tumors.

CONCLUSIONS

CAR-TEAM cells enable modification of tumor stroma, leading to increased elimination of PDAC tumors. This approach represents a promising treatment option for pancreatic cancer.

Interplay between B7-H3 and HLA class I in the clinical course of pancreatic ductal adenocarcinoma

Human leukocyte antigen (HLA) class I defects are associated with cancer progression. However, their prognostic significance is controversial and may be modulated by immune checkpoints. Here, we investigated whether the checkpoint B7-H3 modulates the relationship between HLA class I and pancreatic ductal adenocarcinoma (PDAC) prognosis. PDAC tumors were analyzed for the expression of B7-H3, HLA class I, HLA class II molecules, and for the presence of tumor-infiltrating immune cells. We observed defective HLA class I and HLA class II expressions in 75% and 59% of PDAC samples, respectively. HLA class I and B7-H3 expression were positively related at mRNA and protein level, potentially because of shared regulation by RELA, a sub-unit of NF-kB. High B7-H3 expression and low CD8+ T cell density were indicators of poor survival, while HLA class I was not. Defective HLA class I expression was associated with unfavorable survival only in patients with low B7-H3 expression. Favorable survival was observed only when HLA class I expression was high and B7-H3 expression low. Our results provide the rationale for targeting B7-H3 in patients with PDAC tumors displaying high HLA class I levels.

Tumor cell-based liquid biopsy using high-throughput microfluidic enrichment of entire leukapheresis product

Circulating Tumor Cells (CTCs), interrogated by sampling blood from patients with cancer, contain multiple analytes, including intact RNA, high molecular weight DNA, proteins, and metabolic markers. However, the clinical utility of tumor cell-based liquid biopsy has been limited since CTCs are very rare, and current technologies cannot process the blood volumes required to isolate a sufficient number of tumor cells for in-depth assays. We previously described a high-throughput microfluidic prototype utilizing high-flow channels and amplification of cell sorting forces through magnetic lenses. Here, we apply this technology to analyze patient-derived leukapheresis products, interrogating a mean blood volume of 5.83 liters from patients with metastatic cancer, with a median of 2,799 CTCs purified per patient. Isolation of many CTCs from individual patients enables characterization of their morphological and molecular heterogeneity, including cell and nuclear size and RNA expression. It also allows robust detection of gene copy number variation, a definitive cancer marker with potential diagnostic applications. High-volume microfluidic enrichment of CTCs constitutes a new dimension in liquid biopsies.

Ultrasensitive Detection of Circulating LINE-1 ORF1p as a Specific Multicancer Biomarker

Improved biomarkers are needed for early cancer detection, risk stratification, treatment selection, and monitoring treatment response. Although proteins can be useful blood-based biomarkers, many have limited sensitivity or specificity for these applications. Long INterspersed Element-1 (LINE-1) open reading frame 1 protein (ORF1p) is a transposable element protein overexpressed in carcinomas and high-risk precursors during carcinogenesis with negligible expression in normal tissues, suggesting ORF1p could be a highly specific cancer biomarker. To explore ORF1p as a blood-based biomarker, we engineered ultrasensitive digital immunoassays that detect mid-attomolar (10-17 mol/L) ORF1p concentrations in plasma across multiple cancers with high specificity. Plasma ORF1p shows promise for early detection of ovarian cancer, improves diagnostic performance in a multianalyte panel, provides early therapeutic response monitoring in gastroesophageal cancers, and is prognostic for overall survival in gastroesophageal and colorectal cancers. Together, these observations nominate ORF1p as a multicancer biomarker with potential utility for disease detection and monitoring.

SIGNIFICANCE

The LINE-1 ORF1p transposon protein is pervasively expressed in many cancers and is a highly specific biomarker of multiple common, lethal carcinomas and their high-risk precursors in tissue and blood. Ultrasensitive ORF1p assays from as little as 25 μL plasma are novel, rapid, cost-effective tools in cancer detection and monitoring. See related commentary by Doucet and Cristofari, p. 2502. This article is featured in Selected Articles from This Issue, p. 2489.

Dysregulated Repeat Element Viral-like Immune Response in Hepatocellular Carcinoma

Purpose

Dysregulation of viral-like repeat RNAs are a common feature across many malignancies that are linked with immunological response, but the characterization of these in hepatocellular carcinoma (HCC) is understudied. In this study, we performed RNA in situ hybridization (RNA-ISH) of different repeat RNAs, immunohistochemistry (IHC) for immune cell subpopulations, and spatial transcriptomics to understand the relationship of HCC repeat expression, immune response, and clinical outcomes.

Experimental Design

RNA-ISH for LINE1, HERV-K, HERV-H, and HSATII repeats and IHC for T-cell, Treg, B-cell, macrophage, and immune checkpoint markers were performed on 43 resected HCC specimens. Spatial transcriptomics on tumor and vessel regions of interest was performed on 28 specimens from the same cohort.

Results

High HERV-K and high LINE1 expression were both associated with worse overall survival. There was a positive correlation between LINE1 expression and FOXP3 T-regulatory cells (r = 0.51 p < 0.001) as well as expression of the TIM3 immune checkpoint (r = 0.34, p = 0.03). Spatial transcriptomic profiling of HERV-K high and LINE-1 high tumors identified elevated expression of multiple genes previously associated with epithelial mesenchymal transition, cellular proliferation, and worse overall prognosis in HCC including SSX1, MAGEC2, and SPINK1.

Conclusion

Repeat RNAs may serve as useful prognostic biomarkers in HCC and may also serve as novel therapeutic targets. Additional study is needed to understand the mechanisms by which repeat RNAs impact HCC tumorigenesis.

The landscape of human SVA retrotransposons

SINE-VNTR-Alu (SVA) retrotransposons are evolutionarily young and still-active transposable elements (TEs) in the human genome. Several pathogenic SVA insertions have been identified that directly mutate host genes to cause neurodegenerative and other types of diseases. However, due to their sequence heterogeneity and complex structures as well as limitations in sequencing techniques and analysis, SVA insertions have been less well studied compared to other mobile element insertions. Here, we identified polymorphic SVA insertions from 3646 whole-genome sequencing (WGS) samples of >150 diverse populations and constructed a polymorphic SVA insertion reference catalog. Using 20 long-read samples, we also assembled reference and polymorphic SVA sequences and characterized the internal hexamer/variable-number-tandem-repeat (VNTR) expansions as well as differing SVA activity for SVA subfamilies and human populations. In addition, we developed a module to annotate both reference and polymorphic SVA copies. By characterizing the landscape of both reference and polymorphic SVA retrotransposons, our study enables more accurate genotyping of these elements and facilitate the discovery of pathogenic SVA insertions.

Spatial transcriptomics reveals distinct tissue niches linked with steroid responsiveness in acute gastrointestinal GVHD

Severe acute graft-versus-host disease (aGVHD) is associated with significant mortality and morbidity, especially in steroid-resistant (SR) cases. Spatial transcriptomic technology can elucidate tissue-based interactions in vivo and possibly identify predictors of treatment response. Tissue sections from 32 treatment-naïve patients with biopsy-confirmed lower gastrointestinal (GI) aGVHD were obtained. The GeoMx digital spatial profiler was used to capture transcriptome profiles of >18 000 genes from different foci of immune infiltrates, colonic epithelium, and vascular endothelium. Each tissue compartment sampled showed 2 distinct clusters that were analyzed for differential expression and spatially resolved correlation of gene signatures. Classic cell-mediated immunity signatures, normal differentiated epithelial cells, and inflamed vasculature dominated foci sampled from steroid-sensitive cases. In contrast, a neutrophil predominant noncanonical inflammation with regenerative epithelial cells and some indication of angiogenic endothelial response was overrepresented in areas from SR cases. Evaluation of potential prognostic biomarkers identified ubiquitin specific peptidase 17-like (USP17L) family of genes as being differentially expressed in immune cells from patients with worsened survival. In summary, we demonstrate distinct tissue niches with unique gene expression signatures within lower GI tissue from patients with aGVHD and provide evidence of a potential prognostic biomarker.

Proceedings of the inaugural Dark Genome Symposium: November 2022

In November 2022 the first Dark Genome Symposium was held in Boston, USA. The meeting was hosted by Rome Therapeutics and Enara Bio, two biotechnology companies working on translating our growing understanding of this vast genetic landscape into therapies for human disease. The spirit and ambition of the meeting was one of shared knowledge, looking to strengthen the network of researchers engaged in the field. The meeting opened with a welcome from Rosana Kapeller and Kevin Pojasek followed by a first session of field defining talks from key academics in the space. A series of panels, bringing together academia and industry views, were then convened covering a wide range of pertinent topics. Finally, Richard Young and David Ting gave their views on the future direction and promise for patient impact inherent in the growing understanding of the Dark Genome.

Quantitative p53 immunostaining aids in the detection of prevalent dysplasia

AIMS

The lack of accepted scoring criteria has precluded the use of p53 in routine practice. We evaluate the utility of automated quantitative p53 analysis in risk stratifying Barrett's oesophagus (BE) patients using non-dysplastic BE (NDBE) biopsies in a multicentric cohort of BE progressor (P) and non-progressor (NP) patients.

METHODS

NDBE biopsies prior to the diagnosis of advanced neoplasia from 75 BE-P, and index and last surveillance biopsies from 148 BE-NP were stained for p53, and scored digitally as 1+, 2+ and 3+. A secondary cohort of 30 BE-P was evaluated.

RESULTS

Compared with BE-NP, BE-P was predominantly men (p=0.001), ≥55 years of age (p=0.008), with longer BE segments (71% vs 33%; p<0.001). The mean number of 3+p53 positive cells and 3+ positive glands were significantly more in BE-P versus BE-NP NDBE biopsies (175 vs 9.7, p<0.001; 9.8 vs 0.1; p<0.001, respectively). At a cut-off of ≥10 p53 (3+) positive cells, the sensitivity and specificity of the assay to identify BE-P were 39% and 93%. On multivariate analysis, scoring p53 in NDBE biopsies, age, gender and length of BE were significantly associated with neoplastic progression. 54% of patients classified as prevalent dysplasia showed an abnormal p53 immunohistochemical stain. These findings were validated in the secondary cohort.

CONCLUSIONS

Automated p53 analysis in NDBE biopsies serves as a promising tool for assessing BE neoplastic progression and risk stratification. Our study highlights the practical applicability of p53 assay to routine surveillance practice and its ability to detect prevalent dysplasia.

Therapy-associated remodeling of pancreatic cancer revealed by single-cell spatial transcriptomics and optimal transport analysis

In combination with cell intrinsic properties, interactions in the tumor microenvironment modulate therapeutic response. We leveraged high-plex single-cell spatial transcriptomics to dissect the remodeling of multicellular neighborhoods and cell-cell interactions in human pancreatic cancer associated with specific malignant subtypes and neoadjuvant chemotherapy/radiotherapy. We developed Spatially Constrained Optimal Transport Interaction Analysis (SCOTIA), an optimal transport model with a cost function that includes both spatial distance and ligand-receptor gene expression. Our results uncovered a marked change in ligand-receptor interactions between cancer-associated fibroblasts and malignant cells in response to treatment, which was supported by orthogonal datasets, including an ex vivo tumoroid co-culture system. Overall, this study demonstrates that characterization of the tumor microenvironment using high-plex single-cell spatial transcriptomics allows for identification of molecular interactions that may play a role in the emergence of chemoresistance and establishes a translational spatial biology paradigm that can be broadly applied to other malignancies, diseases, and treatments.

Cancer cells co-evolve with retrotransposons to mitigate viral mimicry

Overexpression of repetitive elements is an emerging hallmark of human cancers ¹ . Diverse repeats can mimic viruses by replicating within the cancer genome through retrotransposition, or presenting pathogen-associated molecular patterns (PAMPs) to the pattern recognition receptors (PRRs) of the innate immune system ^2-5 . Yet, how specific repeats affect tumor evolution and shape the tumor immune microenvironment (TME) in a pro- or anti-tumorigenic manner remains poorly defined. Here, we integrate whole genome and total transcriptome data from a unique autopsy cohort of multiregional samples collected in pancreatic ductal adenocarcinoma (PDAC) patients, into a comprehensive evolutionary analysis. We find that more recently evolved S hort I nterspersed N uclear E lements (SINE), a family of retrotransposable repeats, are more likely to form immunostimulatory double-strand RNAs (dsRNAs). Consequently, younger SINEs are strongly co-regulated with RIG-I like receptor associated type-I interferon genes but anti-correlated with pro-tumorigenic macrophage infiltration. We discover that immunostimulatory SINE expression in tumors is regulated by either L ong I nterspersed N uclear E lements 1 (LINE1/L1) mobility or ADAR1 activity in a TP53 mutation dependent manner. Moreover, L1 retrotransposition activity tracks with tumor evolution and is associated with TP53 mutation status. Altogether, our results suggest pancreatic tumors actively evolve to modulate immunogenic SINE stress and induce pro-tumorigenic inflammation. Our integrative, evolutionary analysis therefore illustrates, for the first time, how dark matter genomic repeats enable tumors to co-evolve with the TME by actively regulating viral mimicry to their selective advantage.

Abstract 2513: Distinct cancer-intrinsic mechanisms mediate nerve recruitment/outgrowth versus perineural invasion

Pancreatic ductal adenocarcinoma (PDAC) features abundant perineural invasion (PNI). Intra-tumoral nerves play critical roles in cancer initiation, progression, recurrence, treatment-resistance, metastasis, and mortality for many malignancies but the diverse molecular mechanisms underlying tumor-nerve crosstalk remain largely unknown—hindering the development of therapies targeting this key pathological process. To address this gap, we performed whole transcriptome digital spatial profiling on twelve custom tissue microarrays (n=288 cores) derived from intratumorally-matched regions with and without PNI in primary PDAC specimens (n=31 patients) and independently measured gene expression from cancer cells, fibroblasts, and nerves. We undertook a differential gene expression analysis comparing malignant cells in PNI-present and -absent regions. We quantified the growth kinetics of dorsal root ganglia (DRG) sensory neurons cultured with exogenous candidate proteins, which validated that some candidates augment (e.g., Lgals1) and others inhibit (e.g., Sema3b) neurite outgrowth. Next, we mapped our previously discovered malignant cell programs onto the epithelial segments and observed a significant enrichment of the mesenchymal, basal-like, and neural-like progenitor (NRP) programs versus depletion of the classical program in PNI-present regions. To determine the effects of malignant subtype on nerve outgrowth, we engineered isogenic KrasG12D/+;Trp53FL/FL;Rosa26-dCas9-VPR (KP;dCas9-VPR) organoids to overexpress the master transcription factors (TFs) for each malignant subtype (e.g., Gata6 for classical, Glis3 for NRP). We then performed the DRG neuronal outgrowth assay using conditioned media from each subtype-specific organoid line versus an off-target control and observed that the classical line suppressed neurite outgrowth, the mesenchymal and basal-like lines were neutral, and the NRP line enhanced neurite outgrowth dynamics comparable to the Ngf positive control. Taken together, our findings suggest that the mechanisms underlying nerve recruitment/outgrowth and perineural invasion may be partly decoupled. To further test this hypothesis, we are performing transwell invasion assays comparing KP;dCas9-VPR cancer cell lines that overexpress each of the candidate PNI-associated genes and malignant subtype TFs. We anticipate that this study will transform our understanding of how cancer cells and the peripheral nervous system collaborate and guide prioritization for therapeutic intervention in the burgeoning cancer neuroscience field.

Citation Format: William L. Hwang, Jennifer Su, Carina Shiau, Peter L. Wang, Jimmy A. Guo, Nicole A. Lester, Jaimie L. Barth, Hannah I. Hoffman, Prajan Divakar, Jason W. Reeves, Eric Miller, Joseph M. Beechem, Andrew J. Aguirre, Lei Zheng, David T. Ting, Mari Mino-Kenudson, Tyler Jacks. Distinct cancer-intrinsic mechanisms mediate nerve recruitment/outgrowth versus perineural invasion [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 2513.

Ultrasensitive detection of circulating LINE-1 ORF1p as a specific multi-cancer biomarker

Improved biomarkers are needed for early cancer detection, risk stratification, treatment selection, and monitoring treatment response. While proteins can be useful blood-based biomarkers, many have limited sensitivity or specificity for these applications. Long INterspersed Element-1 (LINE-1, L1) open reading frame 1 protein (ORF1p) is a transposable element protein overexpressed in carcinomas and high-risk precursors during carcinogenesis with negligible detectable expression in corresponding normal tissues, suggesting ORF1p could be a highly specific cancer biomarker. To explore the potential of ORF1p as a blood-based biomarker, we engineered ultrasensitive digital immunoassays that detect mid-attomolar (10-17 M) ORF1p concentrations in patient plasma samples across multiple cancers with high specificity. Plasma ORF1p shows promise for early detection of ovarian cancer, improves diagnostic performance in a multi-analyte panel, and provides early therapeutic response monitoring in gastric and esophageal cancers. Together, these observations nominate ORF1p as a multi-cancer biomarker with potential utility for disease detection and monitoring.

SULT1A1-dependent sulfonation of alkylators is a lineage-dependent vulnerability of liver cancers

Adult liver malignancies, including intrahepatic cholangiocarcinoma and hepatocellular carcinoma, are the second leading cause of cancer-related deaths worldwide. Most individuals are treated with either combination chemotherapy or immunotherapy, respectively, without specific biomarkers for selection. Here using high-throughput screens, proteomics and in vitro resistance models, we identify the small molecule YC-1 as selectively active against a defined subset of cell lines derived from both liver cancer types. We demonstrate that selectivity is determined by expression of the liver-resident cytosolic sulfotransferase enzyme SULT1A1, which sulfonates YC-1. Sulfonation stimulates covalent binding of YC-1 to lysine residues in protein targets, enriching for RNA-binding factors. Computational analysis defined a wider group of structurally related SULT1A1-activated small molecules with distinct target profiles, which together constitute an untapped small-molecule class. These studies provide a foundation for preclinical development of these agents and point to the broader potential of exploiting SULT1A1 activity for selective targeting strategies.

TGF-β in the microenvironment induces a physiologically occurring immune-suppressive senescent state

TGF-β induces senescence in embryonic tissues. Whether TGF-β in the hypoxic tumor microenvironment (TME) induces senescence in cancer and how the ensuing senescence-associated secretory phenotype (SASP) remodels the cellular TME to influence immune checkpoint inhibitor (ICI) responses are unknown. We show that TGF-β induces a deeper senescent state under hypoxia than under normoxia; deep senescence correlates with the degree of E2F suppression and is marked by multinucleation, reduced reentry into proliferation, and a distinct 14-gene SASP. Suppressing TGF-β signaling in tumors in an immunocompetent mouse lung cancer model abrogates endogenous senescent cells and suppresses the 14-gene SASP and immune infiltration. Untreated human lung cancers with a high 14-gene SASP display immunosuppressive immune infiltration. In a lung cancer clinical trial of ICIs, elevated 14-gene SASP is associated with increased senescence, TGF-β and hypoxia signaling, and poor progression-free survival. Thus, TME-induced senescence may represent a naturally occurring state in cancer, contributing to an immune-suppressive phenotype associated with immune therapy resistance.

FGFR mRNA Expression in Cholangiocarcinoma and Its Correlation with FGFR2 Fusion Status and Immune Signatures

PURPOSE

Selective FGFR inhibitors are effective against cholangiocarcinomas that harbor gene alterations in FGFR2. Clinical trials suggest that expression of wild-type FGFR mRNA can predict sensitivity to FGFR inhibitors, but this biomarker has not been well characterized in cholangiocarcinoma. This study explores the prevalence of FGFR mRNA overexpression in cholangiocarcinoma, its role in predicting sensitivity to FGFR inhibitors, and its association with immune markers.

EXPERIMENTAL DESIGN

Tissue microarrays of intrahepatic (ICC) and extrahepatic cholangiocarcinomas (ECC) resected between 2004 and 2015 were used to evaluate FGFR1-4 mRNA expression levels by RNA in situ hybridization (ISH). Expression levels of FGFR2 mRNA were correlated with FGFR2 fusion status and with patient outcomes. Immune markers expression was assessed by IHC and CSF1 and CSF1 receptor expression were examined by RNA ISH.

RESULTS

Among 94 patients with resected cholangiocarcinoma, the majority had ICC (77%). FGFR2 fusions were identified in 23% of ICCs and 5% of ECCs. High levels of FGFR mRNA in FGFR2 fusion-negative ICC/ECC were seen for: FGFR1 (ICC/ECC: 15%/0%), FGFR2 (ICC/ECC: 57%/0%), FGFR3 (ICC/ECC: 53%/18%), and FGFR4 (ICC/ECC: 32%/0%). Overall, 62% of fusion-negative cholangiocarcinomas showed high levels of FGFR mRNA. In patients with advanced FGFR2 fusion-positive ICC, high levels of FGFR2 mRNA did not correlate with clinical benefit. FGFR2 fusion-positive tumors showed a paucity of PD-L1 on tumor cells.

CONCLUSIONS

FGFR mRNA overexpression occurs frequently in cholangiocarcinoma in the absence of genetic alterations in FGFR. This study identifies a molecular subpopulation in cholangiocarcinoma for which further investigation of FGFR inhibitors is merited outside currently approved indications.

Abstract B016: Quantifying and dissecting pancreatic cancer cell phenotypic plasticity using lineage tracing, single-cell multiomics and CRISPR perturbations reveals novel regulators of plastic states

Pancreatic cancer is a lethal disease in part because tumor cells exist in distinct transcriptional phenotypes (e.g. basal and classical states), each with a selective ability to evade current chemotherapy regimens. Two major mechanisms have been suggested for treatment evasion: 1) intrinsic resistance of certain phenotypes to particular chemotherapy regimens and 2) plasticity of treatment sensitive phenotypes to adopt more resistant phenotypes. However, the relative contribution of these mechanisms to treatment resistance is still poorly understood. Whereas previous work has described the redistribution of tumor cell states under selective treatment pressure, there is no direct evidence that tumor cells exhibit phenotypic plasticity at steady state or with treatment. By leveraging technological advancements in single-cell methods, lineage tracing and functional genomics, we have now shown direct evidence of phenotypic state switching in human pancreatic cancer cell lines. By performing single-cell RNA-seq on 5 barcoded PDAC cell lines over a steady state timecourse and under chemotherapy selective pressure (&gt;600k cells total), we identify unique plasticity phenotypes within these cell lines and infer regulators of these plastic states. We validate the role of several of these regulators using bulk phenotypic CRISPRi screens in these cell lines. We next perform CRISPRi perturbations along with lineage tracing and single-cell multiomics (&gt;300k cells) to dissect the regulatory relationships that underlie these cell states. We identify several novel epithelial and mesenchymal biasing factors, including those with unique roles in the most plastic clones. Collectively, we nominate several regulators that bias PDAC cell states thus posing a paradigm whereby perturbations may be used to homogenize tumor populations towards treatment-sensitive phenotypes. We believe this approach combined with current chemotherapy regimens could benefit pancreatic cancer patients by targeting residual, resistant tumor cells in the localized and metastatic disease settings to improve patient survival.

Citation Format: Arnav Mehta, Lynn Bi, Aziz Al'Khafaji, Martin Jankowiak, Milan Parikh, Mehrtash Babadi, Alex Bloemendal, Marc Schwartz, Glen Munson, Joeseph Chan, Cassandra Burdziak, Elisa Donnard, Ryan Park, Chen Lu, Philippe Rigollet, Andrew Aguirre, Vidya Subramanian, Ray Jones, Eric S. Lander, David T. Ting, Dana Pe'er, Nir Hacohen. Quantifying and dissecting pancreatic cancer cell phenotypic plasticity using lineage tracing, single-cell multiomics and CRISPR perturbations reveals novel regulators of plastic states [abstract]. In: Proceedings of the AACR Special Conference on Pancreatic Cancer; 2022 Sep 13-16; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2022;82(22 Suppl):Abstract nr B016.

Abstract PR004: Characterizing the effects of neoadjuvant therapy in PDAC

Pancreatic ductal adenocarcinoma lethality can be attributed to a combination of rapid metastatic dissemination and intrinsic resistance to conventional therapies. Our prior studies using single cell RNA-seq in pancreatic circulating tumor cells (CTCs) revealed that these “seeds of metastasis” had a biphenotypic state of both epithelial and mesenchymal features. This suggested that CTCs are highly plastic cells that exist in an intermediate state along epithelial (E) to mesenchymal (M) transition (EMT). Analysis of CTC markers in primary PDAC tumors revealed that this subpopulation of cancer cells was concentrated at the interface of tumor cells and stromal cancer associated fibroblasts (CAFs). To evaluate the contribution of CAFs in PDAC EMT plasticity, we utilized patient derived PDAC-CAF co-culture preclinical models demonstrating that PDAC EMT heterogeneity was modulated by the density of CAFs and partially driven by TGF-b secretion. Moreover, analysis of patient derived PDAC tumor spheres treated with FOLFIRINOX (FFX) demonstrated selection or induction of EMT changes that were also found in patient tumors that were resected after neoadjuvant FFX. Altogether, our collective work along with others demonstrates the importance of EMT plasticity in PDAC cell metastatic propensity and the ability to resist cytotoxic chemotherapy. A multi-institutional randomized Phase II trial supported by SU2C-Lustgarten Foundation evaluating the role of the TGF-b modulating activities of losartan on locally advanced PDAC response to chemotherapy and suppression of metastatic dissemination has neared completion. This 3 arm (n=40 per arm) trial of neoadjuvant FFX, FFX + losartan, or FFX+ losartan + nivolumab (anti-PD1). Here, we have performed EMT RNA in situ hybridization in all post-treatment resection specimens to determine if there is a difference in EMT proportions between the different arms of the study. As an orthogonal unbiased approach, we have utilized the NanoString GeoMx Digital Spatial Profiler whole transcriptome assay (18,000+ protein coding genes) on each of the resection specimens with analysis of the tumor cells, CAFs, and immune cell separately for multiple regions of interest in these specimens. This provides deeper molecular insight of the PDAC tumor cells that resisted neoadjuvant therapy, the changes in the surrounding CAFs, and the modulation of immune infiltrates that might differ between each arm of the study. We anticipate the results of the trial will be completed before the conference and initial correlative analysis of spatial transcriptomics and RNA-ISH will be presented.

Citation Format: Amaya Pankaj, Michael J. Raabe, Bidish Patel, Evan R. Lang, Joshua Kocher, Katherine H. Xu, Linda T. Nieman, Alec C. Kimmelman, David P. Ryan, Theodore S. Hong, William L. Hwang, Martin Aryee, David T. Ting. Characterizing the effects of neoadjuvant therapy in PDAC [abstract]. In: Proceedings of the AACR Special Conference on Pancreatic Cancer; 2022 Sep 13-16; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2022;82(22 Suppl):Abstract nr PR004.

Abstract C038: Repeat RNA dysregulation of cellular states in the pancreatic cancer microenvironment

Aberrant transcription of the repeat RNAs is a common feature in epithelial cancers including PDAC, but the function of these non-coding RNAs in cancer development is relatively unexplored. We have found that these repeat RNAs are sensed and replicate like retroviruses, and now have identified the ability of these viral-like elements to be transmitted from cancer cells through extracellular vesicles (EVs). PDAC-derived EVs applied to cancer-associated fibroblasts (CAFs) activates interferon-stimulated genes (ISGs) and is able to drive CAFs towards an inflammatory CAF (iCAF) phenotype with concordant loss of myofibroblast CAF (myCAF) marker genes. Using in-vitro transcription, we demonstrate that individual repeat RNAs (HSATII, HERVK (env), LINE-1 5’UTR and LINE-1 3’UTR) are sufficient to induce ISG response in CAFs with HSATII and HERVK (env) having the most potent ISG response. In contrast, PDAC cells were found to induce epithelial-mesenchymal transition (EMT) with loss of epithelial gene expression. To determine the potential sensor of HSATII repeat RNAs, we utilized CRISPR/Cas9 knockout of the viral RNA sensors RIG-I, MDA5, and MAVS in PDAC and CAF cells. Notably, these sensors were important for PDAC repeat RNA sensing and response, but not in CAF cells. Evaluating the innate immune pathway further downstream, we used genetic knockout of IRF3 with CRISPR/Cas9 knockout and find significant downregulation of key EMT genes that are shared with myCAF markers (ACTA2, FN1, SERPINE1). Interestingly, HSATII RNA activated IRF3 dependent EMT genes in PDAC cells, but induced IRF3 degradation in CAF cells that results in loss of myCAF gene expression. Furthermore, we found that conditioned media from HSATII transfected CAF activates EMT-related gene expression (ACTA2, FN1, SERPINE1) in PDAC cell lines, which indicates an cell extrinsic mechanism to augment EMT induction in PDAC cells. We utilized next generation spatial transcriptomic platforms NanoString GeoMx and CosMx to understand the spatial distribution of repeat RNAs in human PDAC tumors. We find that repeat RNAs can be found as a gradient from PDAC cells to the surround tumor microenvironment consistent with delivery of these RNA species. Analysis of over 300,000 individual cells in 3 PDAC tumor specimens, we find that high repeat PDAC cells have lost epithelial gene expression and high repeat CAFs have lost myCAF gene expression. Altogether, these findings support the “infection” of repeat RNAs disrupts cellular identity in both tumor cells and the CAF microenvironment as a mechanism for tumor progression.

Citation Format: Eunae You, Luli Zou, Patrick Danaher, Ildiko E. Phillips, Michael J. Raabe, Bidish Patel, Amaya Pankaj, Khrystyna North, Sean Kim, Youngmi Kim, Martin Aryee, David T. Ting. Repeat RNA dysregulation of cellular states in the pancreatic cancer microenvironment [abstract]. In: Proceedings of the AACR Special Conference on Pancreatic Cancer; 2022 Sep 13-16; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2022;82(22 Suppl):Abstract nr C038.

Abstract A066: 5z-7-oxozeaenol as an epithelial mesenchymal transition plasticity inhibitor in PDAC

Transcriptional profiling has defined PDAC into distinct classical epithelial (E) or quasi-mesenchymal (QM) subtypes, and these subtypes exist on a continuum of interconverting cell states. This epithelial to mesenchymal transition (EMT) plasticity is thought to be important for the metastatic dissemination and intrinsic resistance to chemotherapy. Our prior work identified TAK1 as an important driver of EMT in mouse circulating tumor cells, and we had identified a compound (5z)-7-oxozeaenol (Oxo) as a small molecule inhibitor of TAK1. Here, we have expanded upon this prior work by evaluating Oxo in 3 human patient derived PDAC cell lines as an EMT plasticity inhibitor. Treatment with Oxo induced a shift of all PDAC cell lines towards an E phenotype in both E and QM PDAC cell lines using RNA-seq and quantitative RNA in situ hybridization. This was associated with increased FFX sensitivity, loss of migratory abilities of all cell lines by Boyden chamber, and decreased invasion in matrigel. Interestingly, genetic TAK1 knockout via CRISPR/Cas9 in these PDAC cell lines showed different effects from Oxo. We found that TAK1 knockout decreased E PDAC cell line migratory behavior and increased FFX sensitivity, but QM PDAC cell lines were not affected by TAK1 knockout. This suggests that QM PDAC cell lines have additional signaling cascades that compensate for TAK1 loss and that Oxo effects across PDAC cell lines are not solely TAK1 dependent. Using phospho-proteomics, in FFX treated PDAC cell lines, we have identified potential candidate pathways in RNA splicing and viral processes that may be additional processes affected by Oxo that merits further investigation. Altogether, we have identified Oxo as a potential broad EMT plasticity inhibitor across PDAC subtypes, and there is a potential role of selective TAK1 inhibitors for E PDAC tumors.

Citation Format: Eunae You, Ildiko Phillips E. Phillips, Richard Y. Ebright, Niyati Desai, Michael J. Raabe, Evan R. Lang, Linda T. Nieman, Soroush Hajizadeh, Johannes Kreuzer, Wilhelm Haas, David T. Ting. 5z-7-oxozeaenol as an epithelial mesenchymal transition plasticity inhibitor in PDAC [abstract]. In: Proceedings of the AACR Special Conference on Pancreatic Cancer; 2022 Sep 13-16; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2022;82(22 Suppl):Abstract nr A066.

Abstract C052: Identifying mediators of perineural invasion in pancreatic cancer using spatial transcriptomics

Intratumoral nerves play important and versatile roles in cancer initiation, progression, recurrence, treatment-resistance, metastasis, morbidity, and mortality for many malignancies but the diverse molecular mechanisms underlying tumor-nerve crosstalk remain largely unknown. One of the differentiating hallmarks of pancreatic ductal adenocarcinoma (PDAC) is an exceptionally high frequency of perineural invasion (PNI), a histopathologic manifestation of tumor-nerve crosstalk whereby cancer cells recruit, migrate towards, and envelop or invade peripheral nerves. Evidence for some neurochemicals/neurotrophins involved in PNI have been uncovered, but most of the underlying work was limited by a lack of cell-type specificity, spatial context, and fragmented focus on individual pathways. To address these shortcomings, we set out to comprehensively identify cell-type specific genes spatially linked to PNI in patient tumors and then dissect the functional roles of these genes through live imaging of dorsal root ganglia (DRG) sensory neurons incubated in conditioned media from cancer cell organoids overexpressing candidate genes via CRISPR activation (CRISPRa). First, we performed whole transcriptome digital spatial profiling (NanoString GeoMx) on twelve custom tissue microarrays (n=288 cores) derived from intratumorally-matched malignant regions with and without PNI in primary resected PDAC specimens (n=31 patients). Differential gene expression (DE) analysis (FDR &lt; 0.001) for PNI demonstrated that for malignant cells there were 271 enriched and 65 depleted genes, and for fibroblasts there were 16 enriched and 27 depleted genes. We further evaluated associations between PNI and expression of malignant subtypes previously identified from single-nucleus RNA-seq applied to 43 primary resected PDAC specimens. We found that malignant cells engaged in PNI were enriched in the mesenchymal, basaloid and neural-like progenitor (NRP) subtypes and depleted in the classical subtype. To test these associations functionally, we generated isogenic murine organoid lines (KrasG12D/+;Trp53FL/FL;R26-dCas9-VPR) overexpressing subtype-driving transcription factors and collected conditioned media. DRG sensory neurons demonstrate enhanced and suppressed growth kinetics when grown in NRP and classical conditioned media, respectively; mesenchymal and basal-like conditioned media do not appear to influence growth kinetics. These results suggest that while mesenchymal, basaloid, and NRP cells likely all play a role in cancer cell invasion of nerves, NRP cells may have an additional role in tumor-nerve tropism. Additional experiments exploring the functional effects of the top enriched and depleted genes from the DE analysis are ongoing. We anticipate that this study will provide a high-resolution understanding of critical intercellular interactions in the PDAC tumor microenvironment that facilitate PNI and tumor-nerve crosstalk more broadly to guide novel therapeutic strategies.

Citation Format: William L. Hwang, Jennifer Su, Jimmy A. Guo, Carina Shiau, Jaimie L. Barth, Hannah I. Hoffman, Prajan Divakar, Jason W. Reeves, Eric Miller, Grissel Cervantes-Jaramillo, William Freed-Pastor, Vanessa Funes, Jennifer Y. Wo, Theodore S. Hong, Carlos Fernandez-del Castillo, Lei Zheng, Andrew J. Aguirre, David T. Ting, Mari Mino-Kenudson, Tyler Jacks. Identifying mediators of perineural invasion in pancreatic cancer using spatial transcriptomics [abstract]. In: Proceedings of the AACR Special Conference on Pancreatic Cancer; 2022 Sep 13-16; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2022;82(22 Suppl):Abstract nr C052.

Abstract C012: Dissecting the reorganization of pancreatic tumor microenvironments after radiation and immunotherapy reveals insights into immunotherapy resistance

Immune checkpoint blockade (ICB) has revolutionized the treatment of many cancers but has been ineffective for the treatment of microsatellite stable (MSS) PDAC. The lack of efficacy of immunotherapies in PDAC is due to: 1) a desmoplastic tumor microenvironment (TME); 2) the presence of suppressive cells, including myeloid derived suppressor cells and regulatory T cells; and 3) the lack of antigen-presenting dendritic cells (DCs) that are important in priming an effective immune response to generate functionally effective tumor antigen-specific T cells. We recently completed a pilot study of dual ICB (Ipilumamab and Nivolumab) with radiation therapy (SBRT 8Gy for 3 fractions) in a cohort of 25 metastatic PDAC patients that had progressed on conventional chemotherapy; this combination conferred an impressive 18% ORR and 29% disease control rate measured on non-irradiated lesions (historical 0% ORR with ICB in PDAC). This led to a phase 2 study in 30 metastatic PDAC patients using this dual modality treatment paradigm. To understand the role of radiation and ICB in altering the PDAC tumor microenvironment we performed single-cell RNA-sequencing and TCR-sequencing (&gt;180k cells), and single-nucleus RNA-sequencing (&gt;300k cells) on 36 tumor biopsies (23 pre-treatment, 13 paired on-treatment between day 10 and 21) from patients undergoing treatment in our phase 2 study. Tumor tissue was taken from distinct tissue sites, including primary tumors in the pancreas, and liver and abdominal wall metastases. We identified distinct tumor cell state distributions within different tissues, and a redistribution of cells from basal/mesenchymal states to classical states after radiation. We identified several state-specific interferon stimulated gene programs thus cataloging distinct responses of epithelial cells with different transcriptional states. Importantly, we found a redistribution of T cells states towards proliferating and exhausted T cells with unique clonality after radiation. Additionally, the myeloid compartment after radiation was enriched for C1QC+ and MHCII+ macrophage subsets, as well as infiltrating CD16/CD16 monocytes and CD14 monocytes, each showing induction of unique sets of interferon stimulated genes (ISGs). We next sought to better understand immunotherapy resistance mechanisms within these PDAC patients despite finding strong ISG induction in several subsets. We analyzed covarying gene programs and identified multicellular communities of cells before and after radiation that underlie interaction networks associated with radiation. Together our data provides the most comprehensive single-cell atlas of paired biopsies to study tumor and immune cell states in the context of radiation and ICB response.

Citation Format: Arnav Mehta, Aparna Parikh, Milan Parikh, Ryan Park, Moshe Sade-Feldman, Lynn Bi, Nicole Carzo, Tarin M. Grillo, Islam Baiev, Olanike Asupoto, Irena Gushterova, Tom LaSalle, Anna Gonye, Emily Blaum, Sebastien Vigneau, Ronan Chaligne, Ana Lako, Thomas Lila, David Nelson, Caroline Porter, Orr Ashenberg, Karthik Jagadesh, William L. Hwang, Christopher Smillie, David P. Ryan, David T. Ting, Theodore Hong, Dana Pe'er, Nir Hacohen. Dissecting the reorganization of pancreatic tumor microenvironments after radiation and immunotherapy reveals insights into immunotherapy resistance [abstract]. In: Proceedings of the AACR Special Conference on Pancreatic Cancer; 2022 Sep 13-16; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2022;82(22 Suppl):Abstract nr C012.

Prospective Phase II Trials Validate the Effect of Neoadjuvant Chemotherapy on Pattern of Recurrence in Pancreatic Adenocarcinoma

OBJECTIVE

The objective of this study was to characterize the patterns of first recurrence after curative-intent resection for pancreatic adenocarcinoma (PDAC).

SUMMARY OF BACKGROUND DATA

We evaluated the first site of recurrence after neoadjuvant treatment as locoregional (LR) or distant metastasis (DM). To validate our findings, we evaluated the pattern from 2 phase II clinical trials evaluating neoadjuvant chemotherapy (NAC) in PDAC.

METHODS

We identified site of first recurrence from a retrospective cohort of patients from 2011 to 2017 treated with NAC followed by chemoradiation and then an operation or an operation first followed by adjuvant therapy, and 2 separate prospective cohorts of patients derived from 2 phase II clinical trials evaluating patients treated with NAC in borderline-resectable and locally advanced PDAC.

RESULTS

In the retrospective cohorts, 160 out of 285 patients (56.1%) recurred after a median disease-free survival (mDFS) of 17.2 months. The pattern of recurrence was DM in 81.9% of patients, versus LR in 11.1%. This pattern was consistent in patients treated with upfront resection and adjuvant chemotherapy (DM 83.0%, LR 16.9%) regardless of margin-involvement (DM 80.1%, LR 19.4%). The use of NAC did not alter pattern of recurrence; 81.7% had DM and 18.3% had LR. This pattern also remained consistent regardless of margin-involvement (DM 94.1%, LR 5.9%). In the Phase II borderline-resectable trial (NCI# 01591733) cohort of 32 patients, the mDFS was 34.2 months. Pattern of recurrence remained predominantly DM (88.9%) versus LR (11.1%). In the Phase II locally-advanced trial (NCI# 01821729) cohort of 34 patients, the mDFS was 30.7 months. Although there was a higher rate of local recurrence in this cohort, pattern of first recurrence remained predominantly DM (66.6%) versus LR (33.3%) and remained consistent independent of margin-status.

CONCLUSIONS

The pattern of recurrence in PDAC is predominantly DM rather than LR, and is consistent regardless of the use of NAC and margin involvement.

Abstract PO016: Spatial transcriptomic profiling to characterize the tumor-vascular interactome of hepatocellular carcinoma

BACKGROUND: The hepatocellular carcinoma (HCC) tumor microenvironment (TME) is composed of a complex ecosystem dominated by cancer cells and the endothelial cells that line tumor blood vessels. Although many genomic drivers have been identified and at least three transcriptional subsets have been proposed, these efforts have not yet led to novel therapies or otherwise significantly impacted management options. Single cell transcriptional profiling has generated deep insights into the multiple heterogeneous cell types within tissues, but the spatial context of these data is lost during single cell processing. Spatial transcriptomic approaches aim to bridge the gap between dissociative single cell technologies and in situ histopathological characterization.METHODS: To gain insight into potential in situ cancer-endothelial crosstalk interactions, we utilized the Nanostring GeoMx spatial transcriptomics platform with the Cancer Transcriptome Atlas ~1800 gene oligonucleotide probe panel to generate tumor (Arginase+) and blood vessel (CD31+) areas of interest (AOI) gene expression profiles from formalin-fixed, paraffin-embedded archival tissue specimens obtained from HCC resection specimens. Oligonucleotides released from each microscopic AOI were then captured, processed by DNA sequencing, and analyzed using custom computational pipelines.RESULTS: Using the 119 ROI containing data from both tumor and vessels that passed quality control filters, we performed unbiased hierarchical clustering of both the tumor and vessel areas of interest (AOI) within each ROI using the most highly variable genes for each AOI set and identified at least 3 clusters within each AOI type (tumor and vessel). Based on gene ontology analysis of the tumor AOIs, the two subsets were distinguished by unique immune and inflammatory-related genes. Analogous ontology-based characterization of the vessel AOIs demonstrated two groups: 1) an interferon-activated, inflamed progenitor, and immune checkpoint-associated cluster; and 2) a TGF-beta and oxidative stress-associated cluster. Notably, both vessel clusters also contained significant numbers of leukocyte genes, concordant with the intimate relationship of the vasculature and immune system. Canonical correlation analysis (CCA) utilizing both the most variable genes within each AOI set showed significant correlated gene sets within tumor AOIs and vessel AOIs, implying biologically significant interactions in multiple signaling pathways.CONCLUSIONS: Spatial transcriptomic profiling enables an understanding of cell-cell interactions in situ that can uncover biologically distinct tumor and blood vessel niches within the HCC microenvironment. Subsequent efforts will be focused on functionally assessing the spatially linked cancer and endothelial cell phenotypes with the goals of developing improved prognostic and predictive biomarkers and generating novel drug targets.

Citation Format: Joseph W Franses, Michael J Raabe, Amaya Pankaj, Bidish Patel, Avril Coley, Irun Bhan, Martin Aryee, David T Ting. Spatial transcriptomic profiling to characterize the tumor-vascular interactome of hepatocellular carcinoma [abstract]. In: Proceedings of the AACR Special Conference: Advances in the Pathogenesis and Molecular Therapies of Liver Cancer; 2022 May 5-8; Boston, MA. Philadelphia (PA): AACR; Clin Cancer Res 2022;28(17_Suppl):Abstract nr PO016.

Satellite repeat RNA expression in epithelial ovarian cancer associates with a tumor-immunosuppressive phenotype

Aberrant expression of viral-like repeat elements is a common feature of epithelial cancers, and the substantial diversity of repeat species provides a distinct view of the cancer transcriptome. Repeatome profiling across ovarian, pancreatic, and colorectal cell lines identifies distinct clustering independent of tissue origin that is seen with coding gene analysis. Deeper analysis of ovarian cancer cell lines demonstrated that human satellite II (HSATII) satellite repeat expression was highly associated with epithelial-mesenchymal transition (EMT) and anticorrelated with IFN-response genes indicative of a more aggressive phenotype. SATII expression - and its correlation with EMT and anticorrelation with IFN-response genes - was also found in ovarian cancer RNA-Seq data and was associated with significantly shorter survival in a second independent cohort of patients with ovarian cancer. Repeat RNAs were enriched in tumor-derived extracellular vesicles capable of stimulating monocyte-derived macrophages, demonstrating a mechanism that alters the tumor microenvironment with these viral-like sequences. Targeting of HSATII with antisense locked nucleic acids stimulated IFN response and induced MHC I expression in ovarian cancer cell lines, highlighting a potential strategy of modulating the repeatome to reestablish antitumor cell immune surveillance.

Precision Medicine in Pancreatic Cancer: Patient-Derived Organoid Pharmacotyping Is a Predictive Biomarker of Clinical Treatment Response

PURPOSE

Patient-derived organoids (PDO) are a promising technology to support precision medicine initiatives for patients with pancreatic ductal adenocarcinoma (PDAC). PDOs may improve clinical next-generation sequencing (NGS) and enable rapid ex vivo chemotherapeutic screening (pharmacotyping).

EXPERIMENTAL DESIGN

PDOs were derived from tissues obtained during surgical resection and endoscopic biopsies and studied with NGS and pharmacotyping. PDO-specific pharmacotype is assessed prospectively as a predictive biomarker of clinical therapeutic response by leveraging data from a randomized controlled clinical trial.

RESULTS

Clinical sequencing pipelines often fail to detect PDAC-associated somatic mutations in surgical specimens that demonstrate a good pathologic response to previously administered chemotherapy. Sequencing the PDOs derived from these surgical specimens, after biomass expansion, improves the detection of somatic mutations and enables quantification of copy number variants. The detection of clinically relevant mutations and structural variants is improved following PDO biomass expansion. On clinical trial, PDOs were derived from biopsies of treatment-naïve patients prior to treatment with FOLFIRINOX (FFX). Ex vivo PDO pharmacotyping with FFX components predicted clinical therapeutic response in these patients with borderline resectable or locally advanced PDAC treated in a neoadjuvant or induction paradigm. PDO pharmacotypes suggesting sensitivity to FFX components were associated with longitudinal declines of tumor marker, carbohydrate-antigen 19-9 (CA-19-9), and favorable RECIST imaging response.

CONCLUSIONS

PDOs established from tissues obtained from patients previously receiving cytotoxic chemotherapies can be accomplished in a clinically certified laboratory. Sequencing PDOs following biomass expansion improves clinical sequencing quality. High in vitro sensitivity to standard-of-care chemotherapeutics predicts good clinical response to systemic chemotherapy in PDAC. See related commentary by Zhang et al., p. 3176.

Molecular Basis of Extramural Vascular Invasion (EMVI) in Colorectal Carcinoma

BACKGROUND

Extramural vascular invasion (EMVI) is a known poor prognostic factor in colorectal carcinoma; however, its molecular basis has not been defined. This study aimed to assess the expression of molecular markers in EMVI positive colorectal carcinoma to understand their tumor microenvironment.

METHODS

Immunohistochemistry was performed on tissue microarrays of surgically resected colorectal cancer specimens for immunological markers, and BRAFV600E mutation (and on the tissue blocks for mismatch repair proteins). Automated quantification was used for CD8, LAG3, FOXP3, PU1, and CD163, and manual quantification was used for PDL1, HLA I markers (beta-2 microglobulin, HC10), and HLA II. The Wilcoxon rank-sum test was used to compare EMVI positive and negative tumors. A logistic regression model was fitted to assess the predictive effect of biomarkers on EMVI.

RESULTS

There were 340 EMVI positive and 678 EMVI negative chemo naïve tumors. PDL1 was barely expressed on tumor cells (median 0) in the entire cohort. We found a significantly lower expression of CD8, LAG3, FOXP3, PU1 cells, PDL1 positive macrophages, and beta-2 microglobulin on tumor cells in the EMVI positive subset (p ≤ 0.001). There was no association of BRAFV600E or deficient mismatch repair proteins (dMMR) with EMVI. PU1 (OR 0.8, 0.7-0.9) and low PDL1 (OR 1.6, 1.1-2.3) independently predicted EMVI on multivariate logistic regression among all biomarkers examined.

CONCLUSION

There is a generalized blunting of immune response in EMVI positive colorectal carcinoma, which may contribute to a worse prognosis. Tumor-associated macrophages seem to play the most significant role in determining EMVI.

Single-nucleus and spatial transcriptome profiling of pancreatic cancer identifies multicellular dynamics associated with neoadjuvant treatment

Pancreatic ductal adenocarcinoma (PDAC) is a highly lethal and treatment-refractory cancer. Molecular stratification in pancreatic cancer remains rudimentary and does not yet inform clinical management or therapeutic development. Here, we construct a high-resolution molecular landscape of the cellular subtypes and spatial communities that compose PDAC using single-nucleus RNA sequencing and whole-transcriptome digital spatial profiling (DSP) of 43 primary PDAC tumor specimens that either received neoadjuvant therapy or were treatment naive. We uncovered recurrent expression programs across malignant cells and fibroblasts, including a newly identified neural-like progenitor malignant cell program that was enriched after chemotherapy and radiotherapy and associated with poor prognosis in independent cohorts. Integrating spatial and cellular profiles revealed three multicellular communities with distinct contributions from malignant, fibroblast and immune subtypes: classical, squamoid-basaloid and treatment enriched. Our refined molecular and cellular taxonomy can provide a framework for stratification in clinical trials and serve as a roadmap for therapeutic targeting of specific cellular phenotypes and multicellular interactions.

Abstract SY12-04: Multicellular spatial community featuring a novel neuronal-like malignant phenotype is enriched in pancreatic cancer after neoadjuvant chemotherapy and radiotherapy

Pancreatic ductal adenocarcinoma (PDAC) is projected to be the second leading cause of cancer mortality in the United States by 2030. Given that resistance to cytotoxic therapy is pervasive, there is a critical need to elucidate salient gene expression programs and spatial relationships among malignant and stromal cells in the tumor microenvironment (TME), particularly in residual disease. We developed and applied a single-nucleus RNA-seq (snRNA-seq) technique to 43 banked frozen primary PDAC specimens that either received neoadjuvant therapy (n=25) or were treatment-naïve (n=18). We discovered expression programs across malignant cell and fibroblast profiles that formed the basis for a refined molecular taxonomy, including a novel neural-like progenitor (NRP) malignant program enriched with neoadjuvant treatment in tumors and organoids, and associated with the worst prognosis in bulk profiles from independent cohorts.

To elucidate how neoadjuvant treatment and cancer cell- and fibroblast-intrinsic programs modulate the composition of multicellular neighborhoods, we performed spatial profiling with the GeoMx[1] platform (NanoString) on 21 formalin-fixed paraffin-embedded sections using the human whole transcriptome atlas (WTA). Each tumor showed intra-tumoral heterogeneity in tissue architecture and regions of interest (ROIs) with diverse patterns of neoplastic cells, cancer-associated fibroblasts (CAFs), and immune cells were selected for profiling. We deconvolved the WTA data with our snRNA-seq cell type signatures and mapped expression programs onto the tumor architecture to reveal three distinct multicellular neighborhoods, which we annotated as classical, squamoid-basaloid, and treatment-enriched. The observed enrichment in post-treatment residual disease of multiple spatially-defined receptor-ligand interactions and a neighborhood featuring the NRP program, neurotropic CAF program, and CD8+ T cells may open new therapeutic opportunities.

Next, we mapped malignant/CAF programs and immune cell subsets at single-cell spatial resolution by performing spatial molecular imaging (SMI[2]; NanoString CosMx) using a panel of 960 RNA targets on a subset of seven tumors (2 untreated, 5 treated) and captured over 200,000 cells with an average of more than 450 transcripts detected per cell. Correlating ROIs from whole-transcriptome DSP to matched fields of view in kiloplex SMI enabled further dissection of PDAC architecture and treatment-associated remodeling of cell type distributions and receptor-ligand interactions.

Ongoing functional studies have begun to elucidate the key regulatory elements underlying the distinct treatment-associated NRP malignant program and its interactions with the TME. Overall, the complementary combination of snRNA-seq, whole-transcriptome DSP, and kiloplex SMI provides a high-resolution molecular framework that can be harnessed to augment precision oncology efforts in pancreatic cancer.

[1] GeoMx DSP is for Research Use Only and not for use in diagnostic procedures. [2] CosMx SMI is for Research Use Only and not for use in diagnostic procedures.

Citation Format: William L. Hwang, Karthik A. Jagadeesh, Jimmy A. Guo, Hannah I. Hoffman, Carina Shiau, Jennifer Su, Payman Yadollahpour, Jason W. Reeves, Youngmi Kim, Sean Kim, Mark Gregory, Prajan Divakar, Eric Miller, Michael Rhodes, Sarah Warren, Erroll Rueckert, Kit Fuhrman, Daniel R. Zollinger, Robin Fropf, Joseph M. Beechem, Arnav Mehta, Toni Delorey, Cristin McCabe, Jaimie L. Barth, Piotr Zelga, Cristina R. Ferrone, Motaz Qadan, Keith D. Lillemoe, Rakesh K. Jain, Jennifer Y. Wo, Theodore S. Hong, Ramnik Xavier, Orit Rozenblatt-Rosen, Andrew J. Aguirre, Carlos Fernandez-Del Castillo, Andrew S. Liss, Mari Mino-Kenudson, David T. Ting, Tyler Jacks, Aviv Regev. Multicellular spatial community featuring a novel neuronal-like malignant phenotype is enriched in pancreatic cancer after neoadjuvant chemotherapy and radiotherapy [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr SY12-04.

Reverse Transcriptase Inhibition Disrupts Repeat Element Life Cycle in Colorectal Cancer

Altered RNA expression of repetitive sequences and retrotransposition are frequently seen in colorectal cancer, implicating a functional importance of repeat activity in cancer progression. We show the nucleoside reverse transcriptase inhibitor 3TC targets activities of these repeat elements in colorectal cancer preclinical models with a preferential effect in p53-mutant cell lines linked with direct binding of p53 to repeat elements. We translate these findings to a human phase II trial of single-agent 3TC treatment in metastatic colorectal cancer with demonstration of clinical benefit in 9 of 32 patients. Analysis of 3TC effects on colorectal cancer tumorspheres demonstrates accumulation of immunogenic RNA:DNA hybrids linked with induction of interferon response genes and DNA damage response. Epigenetic and DNA-damaging agents induce repeat RNAs and have enhanced cytotoxicity with 3TC. These findings identify a vulnerability in colorectal cancer by targeting the viral mimicry of repeat elements.

SIGNIFICANCE

Colorectal cancers express abundant repeat elements that have a viral-like life cycle that can be therapeutically targeted with nucleoside reverse transcriptase inhibitors (NRTI) commonly used for viral diseases. NRTIs induce DNA damage and interferon response that provide a new anticancer therapeutic strategy. This article is highlighted in the In This Issue feature, p. 1397.

ISL2 is a putative tumor suppressor whose epigenetic silencing reprograms the metabolism of pancreatic cancer

Pancreatic ductal adenocarcinoma (PDA) cells reprogram their transcriptional and metabolic programs to survive the nutrient-poor tumor microenvironment. Through in vivo CRISPR screening, we discovered islet-2 (ISL2) as a candidate tumor suppressor that modulates aggressive PDA growth. Notably, ISL2, a nuclear and chromatin-associated transcription factor, is epigenetically silenced in PDA tumors and high promoter DNA methylation or its reduced expression correlates with poor patient survival. The exogenous ISL2 expression or CRISPR-mediated upregulation of the endogenous loci reduces cell proliferation. Mechanistically, ISL2 regulates the expression of metabolic genes, and its depletion increases oxidative phosphorylation (OXPHOS). As such, ISL2-depleted human PDA cells are sensitive to the inhibitors of mitochondrial complex I in vitro and in vivo. Spatial transcriptomic analysis shows heterogeneous intratumoral ISL2 expression, which correlates with the expression of critical metabolic genes. These findings nominate ISL2 as a putative tumor suppressor whose inactivation leads to increased mitochondrial metabolism that may be exploitable therapeutically.

A standardized definition of placental infection by SARS-CoV-2, a consensus statement from the National Institutes of Health/Eunice Kennedy Shriver National Institute of Child Health and Human Development SARS-CoV-2 Placental Infection Workshop

Pregnant individuals infected with SARS-CoV-2 have higher rates of intensive care unit admission, oxygen requirement, need for mechanical ventilation, and death than nonpregnant individuals. Increased COVID-19 disease severity may be associated with an increased risk of viremia and placental infection. Maternal SARS-CoV-2 infection is also associated with pregnancy complications such as preeclampsia and preterm birth, which can be either placentally mediated or reflected in the placenta. Maternal viremia followed by placental infection may lead to maternal-fetal transmission (vertical), which affects 1% to 3% of exposed newborns. However, there is no agreed-upon or standard definition of placental infection. The National Institutes of Health/Eunice Kennedy Shriver National Institute of Child Health and Human Development convened a group of experts to propose a working definition of placental infection to inform ongoing studies of SARS-CoV-2 during pregnancy. Experts recommended that placental infection be defined using techniques that allow virus detection and localization in placental tissue by one or more of the following methods: in situ hybridization with antisense probe (detects replication) or a sense probe (detects viral messenger RNA) or immunohistochemistry to detect viral nucleocapsid or spike proteins. If the abovementioned methods are not possible, reverse transcription polymerase chain reaction detection or quantification of viral RNA in placental homogenates, or electron microscopy are alternative approaches. A graded classification for the likelihood of placental infection as definitive, probable, possible, and unlikely was proposed. Manuscripts reporting placental infection should describe the sampling method (location and number of samples collected), method of preservation of tissue, and detection technique. Recommendations were made for the handling of the placenta, examination, and sampling and the use of validated reagents and sample protocols (included as appendices).

Abstract PR-002: A phase II pilot trial of nivolumab (N) + albumin bound paclitaxel (AP) + paricalcitol (P) + cisplatin (C) + gemcitabine (G) (NAPPCG) in patients with previously untreated metastatic pancreatic ductal adenocarcinoma (PDAC)

Background: Effective therapy for the treatment of PDAC remains one of the greatest unmet oncology clinical needs. The addition of C to G and AP has shown 71% ORR in a previously reported study [JAMA Oncol. 2019 Oct 3;6(1):125-32]. In preclinical work, vitamin D (Vit D) analog therapy decreases myeloid derived suppressor cells and regulatory T cells, turning PDAC into a more immune favorable microenvironment. This trial combines AP/C/G with Vit D analog P and the anti-PD-1 antibody N as a combination therapy for patients with previously untreated metastatic PDAC. This trial evaluates the efficacy and safety of NAPPCG in that patient population (NCT02754726). Methods: Eligibility criteria include Stage IV PDAC, no prior chemotherapy for systemic disease, KPS ≥ 70, and RECIST 1.1 measurable disease. Doses are AP 125 mg/m2, G 1000 mg/m2, each infused over 30 minutes with C 25 mg/m2 infused over 60 minutes on days 1, 8, 22, and 29 of a 42-day cycle. N is given at a fixed dose of 240 mg as a 60 minute infusion on days 1, 15, and 29. P is given at a fixed dose of 25 µg IV twice weekly. Primary objective was to determine the efficacy of the combination for patients with previously untreated metastatic PDAC through determining CR, ORR, PFS, and OS. The secondary objective was to evaluate safety in patients with previously untreated metastatic PDAC. Exploratory endpoints include evaluating tissue molecular profile as it relates to treatment outcomes. Results: Trial was conducted May 2016 with enrollment completed August 2020. 35 patients have been enrolled in the study and 32 are evaluable (baseline and ≥1 follow up CT scan). Most common drug-related grade (Gr) 3-4 adverse events (AE’s), are thrombocytopenia 76% (gr 3 = 34%, gr 4 = 28%) with no serious bleeding events, anemia 37% (gr 3 = 37%, gr 4 = 0%), and CIPN 11% (gr 3 = 11%, gr 4 = 0%). Immune Related Adverse Events &gt;5% were colitis (gr 3=8.6%, gr 4= 0%) and dermatitis (gr 3=8.6%, gr 4= 0%). By RECIST 1.1 criteria, the best response is 1 CR, 26 PR, 4 SD, 1 PD, yielding an 84% ORR (95% CI = (67%, 95%). Median PFS is 6 months (95% CI = (5, 8)). Median OS is 18 months (95% CI = (13, 22)). Conclusions: Although a small study, the high response rate is encouraging. Evaluation of exploratory endpoints is ongoing. Pursuing this regimen in localized PDAC is warranted due to its high ORR. Supported by grants from the Seena Magowitz Foundation, Mattress Firm, Bristol Myers Squibb, and SU2C.

Citation Format: Erkut Borazanci, Gayle S. Jameson, Sunil Sharma, Frank Tsai, Ronald L. Korn, Lana Caldwell, Karen Ansaldo, David T. Ting, Denise Roe, Anna Bermudez, Daniel D. Von Hoff. A phase II pilot trial of nivolumab (N) + albumin bound paclitaxel (AP) + paricalcitol (P) + cisplatin (C) + gemcitabine (G) (NAPPCG) in patients with previously untreated metastatic pancreatic ductal adenocarcinoma (PDAC) [abstract]. In: Proceedings of the AACR Virtual Special Conference on Pancreatic Cancer; 2021 Sep 29-30. Philadelphia (PA): AACR; Cancer Res 2021;81(22 Suppl):Abstract nr PR-002.

LIN28B alters ribosomal dynamics to promote metastasis in MYCN-driven malignancy

High expression of LIN28B is associated with aggressive malignancy and poor survival. Here, probing MYCN-amplified neuroblastoma as a model system, we showed that LIN28B expression was associated with enhanced cell migration in vitro and invasive and metastatic behavior in murine xenografts. Sequence analysis of the polyribosome fraction of LIN28B-expressing neuroblastoma cells revealed let-7-independent enrichment of transcripts encoding components of the translational and ribosomal apparatus and depletion of transcripts of neuronal developmental programs. We further observed that LIN28B utilizes both its cold shock and zinc finger RNA binding domains to preferentially interact with MYCN-induced transcripts of the ribosomal complex, enhancing their translation. These data demonstrated that LIN28B couples the MYCN-driven transcriptional program to enhanced ribosomal translation, thereby implicating LIN28B as a posttranscriptional driver of the metastatic phenotype.

Radiation therapy enhances immunotherapy response in microsatellite stable colorectal and pancreatic adenocarcinoma in a phase II trial

Overcoming intrinsic resistance to immune checkpoint blockade for microsatellite stable (MSS) colorectal cancer (CRC) and pancreatic ductal adenocarcinoma (PDAC) remains challenging. We conducted a single-arm, non-randomized, phase II trial (NCT03104439) combining radiation, ipilimumab and nivolumab to treat patients with metastatic MSS CRC (n = 40) and PDAC (n = 25) with an Eastern Cooperative Oncology Group (ECOG) performance status of 0 or 1. The primary endpoint was disease control rate (DCR) by intention to treat. DCRs were 25% for CRC (ten of 40; 95% confidence interval (CI), 13-41%) and 20% for PDAC (five of 25; 95% CI, 7-41%). In the per-protocol analysis, defined as receipt of radiation, DCR was 37% (ten of 27; 95% CI, 19-58%) in CRC and 29% (five of 17; 95% CI, 10-56%) in PDAC. Pretreatment biopsies revealed low tumor mutational burden for all samples but higher numbers of natural killer (NK) cells and expression of the HERVK repeat RNA in patients with disease control. This study provides proof of concept of combining radiation with immune checkpoint blockade in immunotherapy-resistant cancers.

Selective targeting of MYC mRNA by stabilized antisense oligonucleotides

MYC is a prolific proto-oncogene driving the malignant behaviors of numerous common cancers, yet potent and selective cell-permeable inhibitors of MYC remain elusive. In order to ultimately realize the goal of therapeutic MYC inhibition in cancer, we have initiated discovery chemistry efforts aimed at inhibiting MYC translation. Here we describe a series of conformationally stabilized synthetic antisense oligonucleotides designed to target MYC mRNA (MYCASOs). To support bioactivity, we designed and synthesized this focused library of MYCASOs incorporating locked nucleic acid (LNA) bases at the 5'- and 3'-ends, a phosphorothioate backbone, and internal DNA bases. Treatment of MYC-expressing cancer cells with MYCASOs leads to a potent decrease in MYC mRNA and protein levels. Cleaved MYC mRNA in MYCASO-treated cells is detected with a sensitive 5' Rapid Amplification of cDNA Ends (RACE) assay. MYCASO treatment of cancer cell lines leads to significant inhibition of cellular proliferation while specifically perturbing MYC-driven gene expression signatures. In a MYC-induced model of hepatocellular carcinoma, MYCASO treatment decreases MYC protein levels within tumors, decreases tumor burden, and improves overall survival. MYCASOs represent a new chemical tool for in vitro and in vivo modulation of MYC activity, and promising therapeutic agents for MYC-addicted tumors.

Conditional Survival in Resected Pancreatic Ductal Adenocarcinoma Patients Treated with Total Neoadjuvant Therapy

BACKGROUND

Dynamic survival data based on time already survived are lacking for resected borderline resectable/locally advanced (BR/LA) pancreatic ductal adenocarcinoma (PDAC) patients who received total neoadjuvant therapy (TNT) with FOLFIRINOX followed by chemoradiation. Conditional survival, i.e., the probability of surviving an additional length of time after having already survived an amount of time, offers such information. We aimed to determine actuarial and conditional overall (OS, COS) and disease-free survival (DFS, CDFS) among this cohort.

METHODS

Clinicopathologic data were retrospectively collected for resected BR/LA PDAC patients who received TNT (2011-2019). COS and CDFS rates were calculated for patients being event (death/recurrence)-free at multiple intervals and by recurrence status.

RESULTS

After a median follow-up of 32.1 months, the 183 patients had a median OS and DFS of 39.1 months and 16.8 months, respectively. COS and CDFS increased as a function of time already survived. The probability of surviving an additional 24 months if a patient survived 2 years post-operatively was 70%, whereas the 4-year actuarial OS was 47%. Similarly, the probability of surviving disease-free an additional 24 months after 2 years was 66%, while actuarial 48-month DFS was 27%. COS for disease-free patients increased further over time. For patients remaining disease-free 12 months post-operatively, BR vs. LA status at diagnosis, tumor ≤ 4 cm at diagnosis, and R0 resection were independent predictors of favorable additional OS and DFS.

CONCLUSIONS

For resected TNT-treated BR/LA PDAC patients, the probability of surviving an additional length of time increases as a function of survival already accrued. Dynamic survival estimates may allow personalized follow-up and counseling.

Refining the Molecular Framework for Pancreatic Cancer with Single-cell and Spatial Technologies

Pancreatic ductal adenocarcinoma (PDAC) is a treatment-refractory malignancy in urgent need of a molecular framework for guiding therapeutic strategies. Bulk transcriptomic efforts over the past decade have yielded two broad consensus subtypes: classical pancreatic/epithelial versus basal-like/squamous/quasi-mesenchymal. Although this binary classification enables prognostic stratification, it does not currently inform the administration of treatments uniquely sensitive to either subtype. Furthermore, bulk mRNA studies are challenged by distinguishing contributions from the neoplastic compartment versus other cell types in the microenvironment, which is accentuated in PDAC given that neoplastic cellularity can be low. The application of single-cell transcriptomics to pancreatic tumors has generally lagged behind other cancer types due in part to the difficulty of extracting high-quality RNA from enzymatically degradative tissue, but emerging studies have and will continue to shed light on intratumoral heterogeneity, malignant-stromal interactions, and subtle transcriptional programs previously obscured at the bulk level. In conjunction with insights provided by single-cell/nucleus dissociative techniques, spatially resolved technologies should also facilitate the contextualization of gene programs and inferred cell-cell interactions within the tumor architecture. Finally, given that patients often receive neoadjuvant chemotherapy and/or chemoradiotherapy even in resectable disease, deciphering the gene programs enriched in or induced by cytotoxic therapy will be crucial for developing insights into complementary treatments aimed at eradicating residual cancer cells. Taken together, single-cell and spatial technologies provide an unprecedented opportunity to refine the foundations laid by prior bulk molecular studies and significantly augment precision oncology efforts in pancreatic cancer.

Abstract 94: Multi-compartment reprogramming and spatially-resolved interactions in frozen pancreatic cancer with and without neoadjuvant chemotherapy and radiotherapy at single-cell resolution

A molecular classification of pancreatic ductal adenocarcinoma (PDAC) that informs clinical management remains elusive. Previously identified bulk expression subtypes in the untreated setting were influenced by contaminating stroma whereas single cell RNA-seq (scRNA-seq) of fresh tumors under-represented key cell types. Two consensus subtypes have arisen from these prior efforts: (1) classical-like, and (2) basal-like. Basal-like tumors were associated with worse survival in the metastatic setting but attempts to refine this binary classification have failed to further stratify patient survival. Here, we developed a robust single-nucleus RNA-seq (snRNA-seq) technique for banked frozen PDAC specimens and studied a cohort of untreated resected primary tumors (n ~ 20). Gene expression programs learned across malignant cell and cancer-associated fibroblast (CAF) profiles uncovered a clinically-relevant molecular taxonomy with improved prognostic stratification compared to prior classifications. Digital spatial profiling revealed an association between malignant cells expressing basal-like programs and greater immune infiltration with relatively fewer macrophages, whereas those exhibiting classical-like programs were linked to inflammatory CAFs and macrophage-predominant microniches. Recent clinical trials have supported the increasing adoption of neoadjuvant therapy to aggressively address the risk of micro-metastatic spread and to circumvent concerns of treatment tolerance in the postoperative setting. There is an urgent need to understand how preoperative treatment impacts residual tumor cells and their interactions with other cell types in the tumor microenvironment to identify additional therapeutic vulnerabilities that can be exploited. Towards this end, we performed snRNA-seq on an unmatched cohort of neoadjuvant-treated resected primary tumors (n ~ 25) with most cases involving FOLFIRINOX chemotherapy followed by chemoradiation. Remarkably, the quality of single-nucleus mRNA profiles was comparable between heavily pre-treated and untreated specimens. We identified differentially expressed genes between treated and untreated samples to infer cell-type specific reprogramming in the residual tumor. This analysis revealed that in the neoadjuvant treatment context, there was lower expression of classical-like phenotypes in malignant cells in favor of basal-like phenotypes associated with TNF-NFkB and interferon signaling as well as the presence of novel acinar and neuroendocrine classical-like states. Our refined molecular taxonomy and spatial resolution may help advance precision oncology in PDAC through informative stratification in clinical trials and insights into compartment-specific therapies.

Citation Format: William L. Hwang, Karthik A. Jagadeesh, Jimmy A. Guo, Hannah I. Hoffman, Payman Yadollahpour, Jason Reeves, Eugene Drokhlyansky, Nicholas Van Wittenberghe, Samouil Farhi, Denis Schapiro, George Eng, Jason M. Schenkel, William A. Freed-Pastor, Orr Ashenberg, Clifton Rodrigues, Domenic Abbondanza, Toni Delorey, Devan Phillips, Jorge Roldan, Debora Ciprani, Marina Kern, Jaimie L. Barth, Daniel R. Zollinger, Kit Fuhrman, Robin Fropf, Joseph Beechem, Colin Weekes, Cristina R. Ferrone, Jennifer Y. Wo, Theodore S. Hong, Orit Rozenblatt-Rosen, Andrew J. Aguirre, Mari Mino-Kenudson, Carlos Fernandez-del- Castillo, Andrew S. Liss, David T. Ting, Tyler Jacks, Aviv Regev. Multi-compartment reprogramming and spatially-resolved interactions in frozen pancreatic cancer with and without neoadjuvant chemotherapy and radiotherapy at single-cell resolution [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 94.

Abstract 381: Elevated PARP7 expression in select cancers identifies a target population for RBN-2397 therapy

PARP7 (TIPARP) is a monoPARP which catalyzes the transfer of single units of ADP-ribose onto substrates to change their function. Its expression is upregulated during cellular stress, including viral infection or through the activation of the aryl hydrocarbon receptor after exposure to cigarette smoke. We and others have shown that PARP7 activity suppresses the Type I interferon (IFN) response following activation by cytosolic nucleic acid sensing pathways. RBN-2397 is a first-in class PARP7 inhibitor, which induces cancer cell autonomous and immune stimulatory effects in preclinical models through enhanced Type I IFN signaling in cancer cells. Here we describe the presence of PARP7 genomic amplification with corresponding increased mRNA expression across select cancers. Elevated PARP7 expression or amplification may identify cancer patients that could derive benefit from treatment with RBN-2397.

In characterizing PARP7 copy number and mRNA expression from The Cancer Genome Atlas (TCGA) database, we found the presence of PARP7 copy number amplification in a subset of tumor types, particularly those of squamous histology, as well as ovarian cancer that corresponded to higher mRNA expression levels. High PARP7 expression correlated with poor survival in squamous cancers, while it had no effect on survival in ovarian cancer. Interestingly, tumor types with high PARP7 expression also expressed higher levels of baseline interferon stimulated genes (ISGs). This parallels our findings that cancer cell lines with higher ISGs at baseline are more responsive to PARP7 inhibition.

To explore PARP7 copy number variations (CNVs) in advanced cancer patients, we queried the FoundationCore® (Foundation Medicine, Inc) database. Similar to TCGA, squamous cancers as well as ovarian, breast, and pancreatic ductal adenocarcinoma (PDAC) were among the tumor types with PARP7 amplifications. Moreover, PARP7 was found to be amplified both on the background of chromosome 3q arm-level amplifications as well as focally.

Congruent to our analysis of PARP7 amplifications, we evaluated PARP7 mRNA expression in both squamous and non-squamous non-small cell lung carcinoma (NSCLC), breast cancer, and PDAC primary tumor samples. Using a validated RNAscope ISH probe set, we analyzed over 1,000 patient samples and found that PARP7 was more highly expressed in tumor cells relative to corresponding normal tissues. Overall, there were varying levels of PARP7 expression across samples with higher expression levels of PARP7 in tumor cells, compared to stroma, across all cancers examined.

Our analyses describing the presence of PARP7 amplifications as well as high expression levels in several tumor types including NSCLC, breast, and PDAC provides evidence for the therapeutic relevance of PARP7 inhibition and highlights potential patient selection strategies to identify those patients more likely to benefit from RBN-2397 treatment.

Citation Format: Jodie Wong, Kristy Kuplast-Barr, Ryan P. Abo, Anupriya S. Kulkarni, Linda T. Nieman, Azfar Neyaz, Katherine H. Xu, Radwa Sharaf, Lee A. Albacker, David T. Ting, Heike Keilhack, Kristen A. McEachern. Elevated PARP7 expression in select cancers identifies a target population for RBN-2397 therapy [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 381.

Tumor Microenvironment Immune Response in Pancreatic Ductal Adenocarcinoma Patients Treated With Neoadjuvant Therapy

BACKGROUND

Neoadjuvant folinic acid, fluorouracil, irinotecan, and oxaliplatin (FOLFIRINOX) and chemoradiation have been used to downstage borderline and locally advanced pancreatic ductal adenocarcinoma (PDAC). Whether neoadjuvant therapy-induced tumor immune response contributes to the improved survival is unknown. Therefore, we evaluated whether neoadjuvant therapy induces an immune response towards PDAC.

METHODS

Clinicopathological variables were collected for surgically resected PDACs at the Massachusetts General Hospital (1998-2016). Neoadjuvant regimens included FOLFIRINOX with or without chemoradiation, proton chemoradiation (25 Gy), photon chemoradiation (50.4 Gy), or no neoadjuvant therapy. Human leukocyte antigen (HLA) class I and II expression and immune cell infiltration (CD4+, FoxP3+, CD8+, granzyme B+ cells, and M2 macrophages) were analyzed immunohistochemically and correlated with clinicopathologic variables. The antitumor immune response was compared among neoadjuvant therapy regimens. All statistical tests were 2-sided.

RESULTS

Two hundred forty-eight PDAC patients were included. The median age was 64 years and 50.0% were female. HLA-A defects were less frequent in the FOLFIRINOX cohort (P = .006). HLA class II expression was lowest in photon and highest in proton patients (P = .02). The FOLFIRINOX cohort exhibited the densest CD8+ cell infiltration (P < .001). FOLFIRINOX and proton patients had the highest CD4+ and lowest T regulatory (FoxP3+) cell density, respectively. M2 macrophage density was statistically significantly higher in the treatment-naïve group (P < .001) in which dense M2 macrophage infiltration was an independent predictor of poor overall survival.

CONCLUSIONS

Neoadjuvant FOLFIRINOX with or without chemoradiation may induce immunologically relevant changes in the tumor microenvironment. It may reduce HLA-A defects, increase CD8+ cell density, and decrease T regulatory cell and M2 macrophage density. Therefore, neoadjuvant FOLFIRINOX therapy may benefit from combinations with checkpoint inhibitors, which can enhance patients' antitumor immune response.

Glioma-Derived miRNA-Containing Extracellular Vesicles Induce Angiogenesis by Reprogramming Brain Endothelial Cells

Glioblastoma (GBM) is characterized by aberrant vascularization and a complex tumor microenvironment. The failure of anti-angiogenic therapies suggests pathways of GBM neovascularization, possibly attributable to glioblastoma stem cells (GSCs) and their interplay with the tumor microenvironment. It has been established that GSC-derived extracellular vesicles (GSC-EVs) and their cargoes are proangiogenic in vitro. To further elucidate EV-mediated mechanisms of neovascularization in vitro, we perform RNA-seq and DNA methylation profiling of human brain endothelial cells exposed to GSC-EVs. To correlate these results to tumors in vivo, we perform histoepigenetic analysis of GBM molecular profiles in the TCGA collection. Remarkably, GSC-EVs and normal vascular growth factors stimulate highly distinct gene regulatory responses that converge on angiogenesis. The response to GSC-EVs shows a footprint of post-transcriptional gene silencing by EV-derived miRNAs. Our results provide insights into targetable angiogenesis pathways in GBM and miRNA candidates for liquid biopsy biomarkers.

Transcriptomic Analysis of Laser Capture Microdissected Tumors Reveals Cancer- and Stromal-Specific Molecular Subtypes of Pancreatic Ductal Adenocarcinoma

PURPOSE

Pancreatic ductal adenocarcinoma (PDAC) lethality is multifactorial; although studies have identified transcriptional and genetic subsets of tumors with different prognostic significance, there is limited understanding of features associated with the minority of patients who have durable remission after surgical resection. In this study, we performed laser capture microdissection (LCM) of PDAC samples to define their cancer- and stroma-specific molecular subtypes and identify a prognostic gene expression signature for short-term and long-term survival.

EXPERIMENTAL DESIGN

LCM and RNA sequencing (RNA-seq) analysis of cancer and adjacent stroma of 19 treatment-naïve PDAC tumors was performed. Gene expression signatures were tested for their robustness in a large independent validation set. An RNA-ISH assay with pooled probes for genes associated with disease-free survival (DFS) was developed to probe 111 PDAC tumor samples.

RESULTS

Gene expression profiling identified four subtypes of cancer cells (C1-C4) and three subtypes of cancer-adjacent stroma (S1-S3). These stroma-specific subtypes were associated with DFS (P = 5.55E-07), with S1 associated with better prognoses when paired with C1 and C2. Thirteen genes were found to be predominantly expressed in cancer cells and corresponded with DFS in a validation using existing RNA-seq datasets. A second validation on an independent cohort of patients using RNA-ISH probes to six of these prognostic genes demonstrated significant association with overall survival (median 17 vs. 25 months; P < 0.02).

CONCLUSIONS

Our results identified specific signatures from the epithelial and the stroma components of PDAC, which add clarity to the nature of PDAC molecular subtypes and may help predict survival.

HIF1A signaling selectively supports proliferation of breast cancer in the brain

Blood-borne metastasis to the brain is a major complication of breast cancer, but cellular pathways that enable cancer cells to selectively grow in the brain microenvironment are poorly understood. We find that cultured circulating tumor cells (CTCs), derived from blood samples of women with advanced breast cancer and directly inoculated into the mouse frontal lobe, exhibit striking differences in proliferative potential in the brain. Derivative cell lines generated by serial intracranial injections acquire selectively increased proliferative competency in the brain, with reduced orthotopic tumor growth. Increased Hypoxia Inducible Factor 1A (HIF1A)-associated signaling correlates with enhanced proliferation in the brain, and shRNA-mediated suppression of HIF1A or drug inhibition of HIF-associated glycolytic pathways selectively impairs brain tumor growth while minimally impacting mammary tumor growth. In clinical specimens, brain metastases have elevated HIF1A protein expression, compared with matched primary breast tumors, and in patients with brain metastases, hypoxic signaling within CTCs predicts decreased overall survival. The selective activation of hypoxic signaling by metastatic breast cancer in the brain may have therapeutic implications.