PanIN and CAF transitions in pancreatic carcinogenesis revealed with spatial data integration.

Spatial transcriptomics (ST) is a powerful new approach to characterize the cellular and molecular architecture of the tumor microenvironment. Previous single-cell RNA-sequencing (scRNA-seq) studies of pancreatic ductal adenocarcinoma (PDAC) have revealed a complex immunosuppressive environment characterized by numerous cancer associated fibroblasts (CAFs) subtypes that contributes to poor outcomes. Nonetheless, the evolutionary processes yielding that microenvironment remain unknown. is a premalignant lesion with potential to develop into PDAC, but the formalin-fixed and paraffin-embedded (FFPE) specimens required for PanIN diagnosis preclude scRNA-seq profiling. We developed a new experimental pipeline for FFPE ST analysis of PanINs that preserves clinical specimens for diagnosis. We further developed novel multi-omics analysis methods for three-fold integration of imaging, ST, and scRNA-seq data to analyze the premalignant microenvironment. The integration of ST and imaging enables automated cell type annotation of ST spots at a single-cell resolution, enabling spot selection and deconvolution for unique cellular components of the tumor microenvironment (TME). Overall, this approach demonstrates that PanINs are surrounded by the same subtypes of CAFs present in invasive PDACs, and that the PanIN lesions are predominantly of the classical PDAC subtype. Moreover, this new experimental and computational protocol for ST analysis suggests a biological model in which CAF-PanIN interactions promote inflammatory signaling in neoplastic cells which transitions to proliferative signaling as PanINs progress to PDAC.

[1]  L. Wood,et al.  Inflammatory Signaling and Fibroblast-Cancer Cell Interactions Transfer from a Harmonized Human Single-cell RNA Sequencing Atlas of Pancreatic Ductal Adenocarcinoma to Organoid Co-Culture , 2022, bioRxiv.

[2]  James M. McFarland,et al.  Microenvironment drives cell state, plasticity, and drug response in pancreatic cancer , 2021, Cell.

[3]  J. Lundeberg,et al.  Genome-wide spatial expression profiling in formalin-fixed tissues , 2021, Cell Genomics.

[4]  H. V. van Laarhoven,et al.  The Cellular Origins of Cancer-Associated Fibroblasts and Their Opposing Contributions to Pancreatic Cancer Growth , 2021, Frontiers in Cell and Developmental Biology.

[5]  S. Lowe,et al.  MYC Levels Regulate Metastatic Heterogeneity in Pancreatic Adenocarcinoma , 2021, Cancer discovery.

[6]  Gustavo S. França,et al.  Exploring tissue architecture using spatial transcriptomics , 2021, Nature.

[7]  E. Jaffee,et al.  From bench to bedside: single-cell analysis for cancer immunotherapy. , 2021, Cancer cell.

[8]  S. Hingorani,et al.  Mesenchymal Lineage Heterogeneity Underlies Nonredundant Functions of Pancreatic Cancer–Associated Fibroblasts , 2021, bioRxiv.

[9]  D. Saur,et al.  Visualization of stem cell activity in pancreatic cancer expansion by direct lineage tracing with live imaging , 2021, eLife.

[10]  M. Guttman,et al.  Integrated spatial genomics reveals global architecture of single nuclei , 2020, Nature.

[11]  Raphael Gottardo,et al.  Integrated analysis of multimodal single-cell data , 2020, Cell.

[12]  C. Iacobuzio-Donahue,et al.  Pancreatic cancer stem cells may define tumor stroma characteristics and recurrence patterns in pancreatic ductal adenocarcinoma , 2020, BMC cancer.

[13]  Toby C. Cornish,et al.  In situ characterization of the 3D microanatomy of the pancreas and pancreatic cancer at single cell resolution , 2020, bioRxiv.

[14]  Evan Z. Macosko,et al.  Highly sensitive spatial transcriptomics at near-cellular resolution with Slide-seqV2 , 2020, Nature Biotechnology.

[15]  Mingyao Li,et al.  Integrating gene expression, spatial location and histology to identify spatial domains and spatially variable genes by graph convolutional network , 2020, bioRxiv.

[16]  Jessica S. Yu,et al.  Multimodal mapping of the tumor and peripheral blood immune landscape in human pancreatic cancer , 2020, Nature Cancer.

[17]  Joon-Oh Park,et al.  Single-cell transcriptome analysis of tumor and stromal compartments of pancreatic ductal adenocarcinoma primary tumors and metastatic lesions , 2020, Genome medicine.

[18]  Q. Nguyen,et al.  stLearn: integrating spatial location, tissue morphology and gene expression to find cell types, cell-cell interactions and spatial trajectories within undissociated tissues , 2020, bioRxiv.

[19]  E. Jaffee,et al.  The tumour microenvironment in pancreatic cancer — clinical challenges and opportunities , 2020, Nature Reviews Clinical Oncology.

[20]  Bora Lim,et al.  Advancing Cancer Research and Medicine with Single-Cell Genomics. , 2020, Cancer cell.

[21]  Steven J. M. Jones,et al.  Transcription phenotypes of pancreatic cancer are driven by genomic events during tumor evolution , 2020, Nature Genetics.

[22]  Itai Yanai,et al.  Integrating microarray-based spatial transcriptomics and single-cell RNA-seq reveals tissue architecture in pancreatic ductal adenocarcinomas , 2020, Nature Biotechnology.

[23]  Gennady Korotkevich,et al.  Fast gene set enrichment analysis , 2019, bioRxiv.

[24]  Gaurav Sharma,et al.  projectR: An R/Bioconductor package for transfer learning via PCA, NMF, correlation, and clustering , 2019, bioRxiv.

[25]  P. Robson,et al.  Cross-species single-cell analysis of pancreatic ductal adenocarcinoma reveals antigen-presenting cancer-associated fibroblasts. , 2019, Cancer discovery.

[26]  Thomas Sherman,et al.  CoGAPS 3: Bayesian non-negative matrix factorization for single-cell analysis with asynchronous updates and sparse data structures , 2019, BMC Bioinformatics.

[27]  Yun-Gui Yang,et al.  Single-cell RNA-seq highlights intra-tumoral heterogeneity and malignant progression in pancreatic ductal adenocarcinoma , 2019, Cell Research.

[28]  Evan Z. Macosko,et al.  Slide-seq: A scalable technology for measuring genome-wide expression at high spatial resolution , 2019, Science.

[29]  G. Evan,et al.  Myc instructs and maintains pancreatic adenocarcinoma phenotype , 2019, bioRxiv.

[30]  Andrew J. Hill,et al.  The single cell transcriptional landscape of mammalian organogenesis , 2019, Nature.

[31]  Vincent A. Traag,et al.  From Louvain to Leiden: guaranteeing well-connected communities , 2018, Scientific Reports.

[32]  Christoph Hafemeister,et al.  Comprehensive integration of single cell data , 2018, bioRxiv.

[33]  Carlo Colantuoni,et al.  Decomposing cell identity for transfer learning across cellular measurements, platforms, tissues, and species , 2018, bioRxiv.

[34]  P. Hermann,et al.  Pancreatic cancer stem cells: A state or an entity? , 2018, Seminars in cancer biology.

[35]  T. Ishiwata,et al.  Pancreatic cancer stem cells: features and detection methods , 2018, Pathology & Oncology Research.

[36]  P. Scheet,et al.  Single-Cell Transcriptomics of Pancreatic Cancer Precursors Demonstrates Epithelial and Microenvironmental Heterogeneity as an Early Event in Neoplastic Progression , 2018, Clinical Cancer Research.

[37]  Paul Hoffman,et al.  Integrating single-cell transcriptomic data across different conditions, technologies, and species , 2018, Nature Biotechnology.

[38]  Laleh Haghverdi,et al.  Batch effects in single-cell RNA-sequencing data are corrected by matching mutual nearest neighbors , 2018, Nature Biotechnology.

[39]  Leland McInnes,et al.  UMAP: Uniform Manifold Approximation and Projection for Dimension Reduction , 2018, ArXiv.

[40]  Geoff S. Higgins,et al.  Oxidative Phosphorylation as an Emerging Target in Cancer Therapy , 2018, Clinical Cancer Research.

[41]  B. Sáinz,et al.  The Ever-Evolving Concept of the Cancer Stem Cell in Pancreatic Cancer , 2018, Cancers.

[42]  M. Nagino,et al.  Trefoil factor family 1 expression in the invasion front is a poor prognostic factor associated with lymph node metastasis in pancreatic cancer. , 2017, Pancreatology : official journal of the International Association of Pancreatology (IAP) ... [et al.].

[43]  Hannah A. Pliner,et al.  Reversed graph embedding resolves complex single-cell trajectories , 2017, Nature Methods.

[44]  R. Hruban,et al.  IL2RG, identified as overexpressed by RNA-seq profiling of pancreatic intraepithelial neoplasia, mediates pancreatic cancer growth , 2017, Oncotarget.

[45]  Michael F. Ochs,et al.  PatternMarkers & GWCoGAPS for novel data-driven biomarkers via whole transcriptome NMF , 2016, bioRxiv.

[46]  Long Cai,et al.  seqFISH Accurately Detects Transcripts in Single Cells and Reveals Robust Spatial Organization in the Hippocampus , 2017, Neuron.

[47]  Roland Eils,et al.  Complex heatmaps reveal patterns and correlations in multidimensional genomic data , 2016, Bioinform..

[48]  Jen Jen Yeh,et al.  Virtual microdissection identifies distinct tumor- and stroma-specific subtypes of pancreatic ductal adenocarcinoma , 2015, Nature Genetics.

[49]  Sean C. Bendall,et al.  Data-Driven Phenotypic Dissection of AML Reveals Progenitor-like Cells that Correlate with Prognosis , 2015, Cell.

[50]  P. Linsley,et al.  MAST: a flexible statistical framework for assessing transcriptional changes and characterizing heterogeneity in single-cell RNA sequencing data , 2015, Genome Biology.

[51]  A. Regev,et al.  Spatial reconstruction of single-cell gene expression data , 2015 .

[52]  J. Mesirov,et al.  The Molecular Signatures Database (MSigDB) hallmark gene set collection. , 2015, Cell systems.

[53]  Cole Trapnell,et al.  The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells , 2014, Nature Biotechnology.

[54]  R. Gibbs,et al.  Genetic alterations associated with progression from pancreatic intraepithelial neoplasia to invasive pancreatic tumor. , 2013, Gastroenterology.

[55]  Z. Naito,et al.  Expression of cancer stem cell markers in pancreatic intraepithelial neoplasias and pancreatic ductal adenocarcinomas. , 2012, International journal of oncology.

[56]  M. Kanda,et al.  Presence of somatic mutations in most early-stage pancreatic intraepithelial neoplasia. , 2012, Gastroenterology.

[57]  R. Hwang,et al.  Trefoil Factor 1 Stimulates Both Pancreatic Cancer and Stellate Cells and Increases Metastasis , 2011, Pancreas.

[58]  Jie Ding,et al.  CoGAPS: an R/C++ package to identify patterns and biological process activity in transcriptomic data , 2010, Bioinform..

[59]  F. Schmidt.  Meta-Analysis , 2008 .

[60]  Pablo Tamayo,et al.  Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[61]  Bart De Moor,et al.  BioMart and Bioconductor: a powerful link between biological databases and microarray data analysis , 2005, Bioinform..

[62]  M. Löhr,et al.  Frequency of K-ras mutations in pancreatic intraductal neoplasias associated with pancreatic ductal adenocarcinoma and chronic pancreatitis: a meta-analysis. , 2005, Neoplasia.

[63]  M. Daly,et al.  PGC-1α-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes , 2003, Nature Genetics.

[64]  R H Hruban,et al.  Progression model for pancreatic cancer. , 2000, Clinical cancer research : an official journal of the American Association for Cancer Research.