Novel Approach for Pancreas Transcriptomics Reveals the Cellular Landscape in Homeostasis and Acute Pancreatitis

BACKGROUND & AIMS

Acinar cells produce digestive enzymes that impede transcriptomic characterization of the exocrine pancreas. Thus, single-cell RNA-sequencing studies of the pancreas underrepresent acinar cells relative to histological expectations, and a robust approach to capture pancreatic cell responses in disease states is needed. We sought to innovate a method that overcomes these challenges to accelerate study of the pancreas in health and disease.

METHODS

We leverage FixNCut, a single-cell RNA-sequencing approach in which tissue is reversibly fixed with dithiobis(succinimidyl propionate) before dissociation and single-cell preparation. We apply FixNCut to an established mouse model of acute pancreatitis, validate findings using GeoMx whole transcriptome atlas profiling, and integrate our data with prior studies to compare our method in both mouse and human pancreas datasets.

RESULTS

FixNCut achieves unprecedented definition of challenging pancreatic cells, including acinar and immune populations in homeostasis and acute pancreatitis, and identifies changes in all major cell types during injury and recovery. We define the acinar transcriptome during homeostasis and acinar-to-ductal metaplasia and establish a unique gene set to measure deviation from normal acinar identity. We characterize pancreatic immune cells, and analysis of T-cell subsets reveals a polarization of the homeostatic pancreas toward type-2 immunity. We report immune responses during acute pancreatitis and recovery, including early neutrophil infiltration, expansion of dendritic cell subsets, and a substantial shift in the transcriptome of macrophages due to both resident macrophage activation and monocyte infiltration.

CONCLUSIONS

FixNCut preserves pancreatic transcriptomes to uncover novel cell states during homeostasis and following pancreatitis, establishing a broadly applicable approach and reference atlas for study of pancreas biology and disease.

scKINETICS: inference of regulatory velocity with single-cell transcriptomics data

MOTIVATION

Transcriptional dynamics are governed by the action of regulatory proteins and are fundamental to systems ranging from normal development to disease. RNA velocity methods for tracking phenotypic dynamics ignore information on the regulatory drivers of gene expression variability through time.

RESULTS

We introduce scKINETICS (Key regulatory Interaction NETwork for Inferring Cell Speed), a dynamical model of gene expression change which is fit with the simultaneous learning of per-cell transcriptional velocities and a governing gene regulatory network. Fitting is accomplished through an expectation-maximization approach designed to learn the impact of each regulator on its target genes, leveraging biologically motivated priors from epigenetic data, gene-gene coexpression, and constraints on cells' future states imposed by the phenotypic manifold. Applying this approach to an acute pancreatitis dataset recapitulates a well-studied axis of acinar-to-ductal transdifferentiation whilst proposing novel regulators of this process, including factors with previously appreciated roles in driving pancreatic tumorigenesis. In benchmarking experiments, we show that scKINETICS successfully extends and improves existing velocity approaches to generate interpretable, mechanistic models of gene regulatory dynamics.

AVAILABILITY AND IMPLEMENTATION

All python code and an accompanying Jupyter notebook with demonstrations are available at http://github.com/dpeerlab/scKINETICS.

Epigenetic plasticity cooperates with cell-cell interactions to direct pancreatic tumorigenesis

The response to tumor-initiating inflammatory and genetic insults can vary among morphologically indistinguishable cells, suggesting as yet uncharacterized roles for epigenetic plasticity during early neoplasia. To investigate the origins and impact of such plasticity, we performed single-cell analyses on normal, inflamed, premalignant, and malignant tissues in autochthonous models of pancreatic cancer. We reproducibly identified heterogeneous cell states that are primed for diverse, late-emerging neoplastic fates and linked these to chromatin remodeling at cell-cell communication loci. Using an inference approach, we revealed signaling gene modules and tissue-level cross-talk, including a neoplasia-driving feedback loop between discrete epithelial and immune cell populations that was functionally validated in mice. Our results uncover a neoplasia-specific tissue-remodeling program that may be exploited for pancreatic cancer interception.

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 (>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 (>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.

Mapping the evolution of T cell states during response and resistance to adoptive cellular therapy

To elucidate mechanisms by which T cells eliminate leukemia, we study donor lymphocyte infusion (DLI), an established immunotherapy for relapsed leukemia. We model T cell dynamics by integrating longitudinal, multimodal data from 94,517 bone marrow-derived single T cell transcriptomes in addition to chromatin accessibility and single T cell receptor sequencing from patients undergoing DLI. We find that responsive tumors are defined by enrichment of late-differentiated T cells before DLI and rapid, durable expansion of early differentiated T cells after treatment, highly similar to "terminal" and "precursor" exhausted subsets, respectively. Resistance, in contrast, is defined by heterogeneous T cell dysfunction. Surprisingly, early differentiated T cells in responders mainly originate from pre-existing and novel clonotypes recruited to the leukemic microenvironment, rather than the infusion. Our work provides a paradigm for analyzing longitudinal single-cell profiling of scenarios beyond adoptive cell therapy and introduces Symphony, a Bayesian approach to infer regulatory circuitry underlying T cell subsets, with broad relevance to exhaustion antagonists across cancers.

Abstract LT008: Mapping the evolution of T cell states during response and resistance to adoptive cellular therapy

Immune therapies have transformed the cancer therapeutic landscape but fail to benefit most patients. To elucidate the underlying mechanisms by which T cells mediate elimination of leukemia, we generated a high-resolution map of longitudinal T cell dynamics within the same tumor microenvironment (TME; bone marrow) during response or resistance to donor lymphocyte infusion (DLI), a widely used immunotherapy for relapsed leukemia. We analyzed 87,939 bone marrow-derived single T cell transcriptomes, along with chromatin accessibility and single T cell receptor clonality profiles, by developing novel machine learning tools for integrating longitudinal and multimodal data. We found that pre-treatment enrichment and post-treatment rapid, durable expansion of ‘terminal’ (TEX) and ‘precursor’ (TPEX) exhausted subsets, respectively, defined DLI response. In contrast to the common, shared pathways marking DLI response, a heterogeneous pattern of T cell dysfunction marked DLI resistance. Unexpectedly, TPEX cells that expanded in responders did not arise from the infusion product but instead from both pre-existing and novel clonotypes recruited to the TME. Further, we introduce a Bayesian method, Symphony, to define the T cell regulatory circuitry and master regulators underlying TEX and TPEX subsets that may be broadly relevant to other exhaustion antagonists across cancers. In conclusion, our data implicate the hierarchy of both TEX and TPEX subsets for immunotherapeutic responses in leukemia, extending the scope of their relevance beyond checkpoint blockade to adoptive cellular therapy. Moreover, our results provocatively suggest that immunologic ‘help’ from DLI, rather than direct transfer of anti-leukemic T cells, drove leukemic remission. Finally, we provide a general analysis paradigm for exploiting temporal single-cell genomic profiling for deep understanding of how immune therapies differentially shape the evolutionary trajectories of the TME in accordance with clinical outcome.

Citation Format: Pavan Bachireddy, Elham Azizi, Cassandra Burdziak, Vinhkhang Nguyen, Christina Ennis, Zi- Ning Choo, Shuqiang Li, Kenneth Livak, Donna Neuberg, Robert Soiffer, Jerome Ritz, Edwin Alyea, Dana Pe'er, Catherine Wu. Mapping the evolution of T cell states during response and resistance to adoptive cellular therapy [abstract]. In: Proceedings of the AACR Virtual Special Conference on the Evolving Tumor Microenvironment in Cancer Progression: Mechanisms and Emerging Therapeutic Opportunities; in association with the Tumor Microenvironment (TME) Working Group; 2021 Jan 11-12. Philadelphia (PA): AACR; Cancer Res 2021;81(5 Suppl):Abstract nr LT008.

A gene-environment-induced epigenetic program initiates tumorigenesis

Tissue damage increases the risk of cancer through poorly understood mechanisms1. In mouse models of pancreatic cancer, pancreatitis associated with tissue injury collaborates with activating mutations in the Kras oncogene to markedly accelerate the formation of early neoplastic lesions and, ultimately, adenocarcinoma2,3. Here, by integrating genomics, single-cell chromatin assays and spatiotemporally controlled functional perturbations in autochthonous mouse models, we show that the combination of Kras mutation and tissue damage promotes a unique chromatin state in the pancreatic epithelium that distinguishes neoplastic transformation from normal regeneration and is selected for throughout malignant evolution. This cancer-associated epigenetic state emerges within 48 hours of pancreatic injury, and involves an 'acinar-to-neoplasia' chromatin switch that contributes to the early dysregulation of genes that define human pancreatic cancer. Among the factors that are most rapidly activated after tissue damage in the pre-malignant pancreatic epithelium is the alarmin cytokine interleukin 33, which recapitulates the effects of injury in cooperating with mutant Kras to unleash the epigenetic remodelling program of early neoplasia and neoplastic transformation. Collectively, our study demonstrates how gene-environment interactions can rapidly produce gene-regulatory programs that dictate early neoplastic commitment, and provides a molecular framework for understanding the interplay between genetic and environmental cues in the initiation of cancer.

Recovering Gene Interactions from Single-Cell Data Using Data Diffusion

Single-cell RNA sequencing technologies suffer from many sources of technical noise, including under-sampling of mRNA molecules, often termed "dropout," which can severely obscure important gene-gene relationships. To address this, we developed MAGIC (Markov affinity-based graph imputation of cells), a method that shares information across similar cells, via data diffusion, to denoise the cell count matrix and fill in missing transcripts. We validate MAGIC on several biological systems and find it effective at recovering gene-gene relationships and additional structures. Applied to the epithilial to mesenchymal transition, MAGIC reveals a phenotypic continuum, with the majority of cells residing in intermediate states that display stem-like signatures, and infers known and previously uncharacterized regulatory interactions, demonstrating that our approach can successfully uncover regulatory relations without perturbations.