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

Pancreatic ductal adenocarcinoma (PDAC) is a devastating malignancy driven by a heterogeneous tumor microenvironment enriched with cancer associated fibroblasts (CAFs) that influence its overall immunosuppressive composition. We examined inflammation in PDAC by leveraging our existing patient-derived organoid (PDO) model and a novel PDO-CAF co-culture. We first identified induction of major histocompatibility complex class II expression following treatment with interferon gamma. In parallel, we collated an atlas of 6 published single-cell RNA-sequencing datasets (174,394 cells) combining 61 PDAC (142,807 cells) and 16 non-malignant samples. By combining in silico modeling and in vitro PDO co-culture, we define a gene expression pattern of inflammatory processes in epithelial tumor cells. Following computational inferences, we examined interactions between epithelial tumor cells and CAFs, focusing on VEGF-A and ITGB1 pathways. This work, integrating computational and biological approaches, highlights the value of convergence to accelerate our understanding of key drivers of PDAC. Statement of Significance We established PDO-CAF co-cultures to model tumor cell interactions and validated discoveries using transfer learning into a single-cell RNA-seq atlas of PDAC tumors. This bidirectional approach from human to experimental systems facilitates interrogation of PDAC biology, including the role of inflammation associated with VEGF-A crosstalk between CAFs and tumor cells.

[1]  E. Jaffee,et al.  Precision medicine in pancreatic cancer: Patient derived organoid pharmacotyping is a predictive biomarker of clinical treatment response. , 2022, Clinical cancer research : an official journal of the American Association for Cancer Research.

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

[3]  Kevin Wei,et al.  Fibroblast pathology in inflammatory diseases. , 2021, The Journal of clinical investigation.

[4]  P. Cahan,et al.  Computational reconstruction of the signalling networks surrounding implanted biomaterials from single-cell transcriptomics , 2021, Nature Biomedical Engineering.

[5]  A. Rao,et al.  Stabilized epithelial phenotype of cancer cells in primary tumors leads to increased colonization of liver metastasis in pancreatic cancer , 2021, Cell reports.

[6]  B. Nicolet,et al.  CD29 Enriches for Cytotoxic Human CD4+ T Cells , 2021, The Journal of Immunology.

[7]  Yinmo Yang,et al.  Single-cell RNA-seq reveals dynamic change in tumor microenvironment during pancreatic ductal adenocarcinoma malignant progression , 2020, EBioMedicine.

[8]  R. Burkhart,et al.  Modeling human pancreatic ductal adenocarcinoma for translational research: current options, challenges, and prospective directions , 2020, Annals of pancreatic cancer.

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

[10]  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.

[11]  C. Conrad,et al.  Single-cell analysis of patient-derived PDAC organoids reveals cell state heterogeneity and a conserved developmental hierarchy , 2020, Nature Communications.

[12]  J. Cameron,et al.  Patient-derived Organoid Pharmacotyping is a Clinically Tractable Strategy for Precision Medicine in Pancreatic Cancer , 2020, Annals of surgery.

[13]  E. Jaffee,et al.  Transfer learning between preclinical models and human tumors identifies a conserved NK cell activation signature in anti-CTLA-4 responsive tumors , 2020, Genome Medicine.

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

[15]  O. Chinot,et al.  Bevacizumab (Avastin®) in cancer treatment: A review of 15 years of clinical experience and future outlook. , 2020, Cancer treatment reviews.

[16]  Erin Helms,et al.  Fibroblast Heterogeneity in the Pancreatic Tumor Microenvironment. , 2020, Cancer discovery.

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

[18]  Maxim N. Artyomov,et al.  MHC-II neoantigens shape tumor immunity and response to immunotherapy , 2019, Nature.

[19]  S. Yamada,et al.  Meflin-positive cancer-associated fibroblasts inhibit pancreatic carcinogenesis. , 2019, Cancer research.

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

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

[22]  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.

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

[24]  Jennifer L Hu,et al.  MULTI-seq: sample multiplexing for single-cell RNA sequencing using lipid-tagged indices , 2019, Nature Methods.

[25]  Martin J. Aryee,et al.  Stromal Microenvironment Shapes the Intratumoral Architecture of Pancreatic Cancer , 2019, Cell.

[26]  A. Maitra,et al.  Cellular heterogeneity during mouse pancreatic ductal adenocarcinoma progression at single-cell resolution , 2019, bioRxiv.

[27]  Hans Clevers,et al.  Generation of Tumor-Reactive T Cells by Co-culture of Peripheral Blood Lymphocytes and Tumor Organoids , 2018, Cell.

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

[29]  S. Byers,et al.  Curation of the Pancreatic Ductal Adenocarcinoma Subset of the Cancer Genome Atlas Is Essential for Accurate Conclusions about Survival-Related Molecular Mechanisms , 2018, Clinical Cancer Research.

[30]  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.

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

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

[33]  J. Aerts,et al.  SCENIC: Single-cell regulatory network inference and clustering , 2017, Nature Methods.

[34]  Hans Clevers,et al.  Distinct populations of inflammatory fibroblasts and myofibroblasts in pancreatic cancer , 2017, The Journal of experimental medicine.

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

[36]  Li Gao,et al.  Proinflammatory Cytokines Increase Vascular Endothelial Growth Factor Expression in Alveolar Epithelial Cells , 2015, Mediators of inflammation.

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

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

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

[40]  P. Muranski,et al.  Naive tumor-specific CD4+ T cells differentiated in vivo eradicate established melanoma , 2010, The Journal of experimental medicine.

[41]  Peter Olson,et al.  Cancer-Associated Fibroblasts Are Activated in Incipient Neoplasia to Orchestrate Tumor-Promoting Inflammation in an NF-kappaB-Dependent Manner. , 2010, Cancer cell.

[42]  E. Birney,et al.  Mapping identifiers for the integration of genomic datasets with the R/Bioconductor package biomaRt , 2009, Nature Protocols.

[43]  H. Pandha,et al.  Loss of expression of antigen‐presenting molecules in human pancreatic cancer and pancreatic cancer cell lines , 2007, Clinical and experimental immunology.

[44]  J. Weitz,et al.  Specific immune recognition of pancreatic carcinoma by patient-derived CD4 and CD8 T cells and its improvement by interferon-gamma. , 2006, International journal of oncology.

[45]  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.

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

[47]  Thomas D. Schmittgen,et al.  Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. , 2001, Methods.

[48]  A. Scarpa,et al.  Expression of MHC class I and class II antigens in pancreatic adenocarcinomas. , 1996, Tissue antigens.