Immune and malignant cell phenotypes of ovarian cancer are determined by distinct mutational processes

High-grade serous ovarian cancer (HGSOC) is an archetypal cancer of genomic instability patterned by distinct mutational processes, intratumoral heterogeneity and intraperitoneal spread. We investigated determinants of immune recognition and evasion in HGSOC to elucidate co- evolutionary processes underlying malignant progression and tumor immunity. Mutational processes and anatomic sites of tumor foci were key determinants of tumor microenvironment cellular phenotypes, inferred from whole genome sequencing, single-cell RNA sequencing, digital histopathology and multiplexed immunofluorescence of 160 tumor sites from 42 treatment-naive HGSOC patients. Homologous recombination-deficient (HRD)-Dup (BRCA1 mutant-like) and HRD- Del (BRCA2 mutant-like) tumors harbored increased neoantigen burden, inflammatory signaling and ongoing immunoediting, reflected in loss of HLA diversity and tumor infiltration with highly- differentiated dysfunctional CD8+ T cells. Foldback inversion (FBI, non-HRD) tumors exhibited elevated TGFβ signaling and immune exclusion, with predominantly naive/stem-like and memory T cells. Our findings implicate distinct immune resistance mechanisms across HGSOC subtypes which can inform future immunotherapeutic strategies. HIGHLIGHTS Multi-region, multi-modal profiling of malignant and immune cell phenotypes in ovarian cancer Anatomic site specificity is a determinant of cancer cell and intratumoral immune phenotypes Tumor mutational processes impact mechanisms of immune control and immune evasion Spatial topology of HR-deficient tumors is defined by immune interactions absent from immune inert HR-proficient subtypes

[1]  Marc J. Williams,et al.  Evolutionary tracking of cancer haplotypes at single-cell resolution , 2021, bioRxiv.

[2]  Marc J. Williams,et al.  The impact of mutational processes on structural genomic plasticity in cancer cells , 2021, bioRxiv.

[3]  K. Robison,et al.  Atezolizumab, Bevacizumab, and Chemotherapy for Newly Diagnosed Stage III or IV Ovarian Cancer: Placebo-Controlled Randomized Phase III Trial (IMagyn050/GOG 3015/ENGOT-OV39) , 2021, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[4]  Chandler D. Gatenbee,et al.  IgA transcytosis and antigen recognition govern ovarian cancer immunity , 2021, Nature.

[5]  Benjamin E. Gross,et al.  OncoTree: A Cancer Classification System for Precision Oncology. , 2021, JCO clinical cancer informatics.

[6]  Yuanhua Huang,et al.  Cellsnp-lite: an efficient tool for genotyping single cells , 2021, bioRxiv.

[7]  J. Reis-Filho,et al.  Mutations in BRCA1 and BRCA2 differentially affect the tumor microenvironment and response to checkpoint blockade immunotherapy , 2020, Nature Cancer.

[8]  R. Weinberg,et al.  Genetically defined syngeneic mouse models of ovarian cancer as tools for the discovery of combination immunotherapy. , 2020, Cancer discovery.

[9]  J. Reis-Filho,et al.  The genetic landscape of metaplastic breast cancers and uterine carcinosarcomas , 2020, Molecular oncology.

[10]  P. Gimotty,et al.  Rapid tumor vaccine using Toll-like receptor-activated ovarian cancer ascites monocytes , 2020, Journal for ImmunoTherapy of Cancer.

[11]  R. Weinberg,et al.  Genetically Defined, Syngeneic Organoid Platform for Developing Combination Therapies for Ovarian Cancer , 2020, bioRxiv.

[12]  B. Monk,et al.  Avelumab in combination with and/or following chemotherapy vs chemotherapy alone in patients with previously untreated epithelial ovarian cancer: Results from the phase 3 javelin ovarian 100 trial , 2020 .

[13]  T. Schumacher,et al.  CD8+ T cell states in human cancer: insights from single-cell analysis , 2020, Nature Reviews Cancer.

[14]  Malachi Griffith,et al.  pVACtools: A Computational Toolkit to Identify and Visualize Cancer Neoantigens , 2020, Cancer Immunology Research.

[15]  M. Baum,et al.  A Hybrid Approach To Hierarchical Density-based Cluster Selection , 2019, 2020 IEEE International Conference on Multisensor Fusion and Integration for Intelligent Systems (MFI).

[16]  Kieran R. Campbell,et al.  Probabilistic cell-type assignment of single-cell RNA-seq for tumor microenvironment profiling , 2019, Nature Methods.

[17]  L. Saal,et al.  Whole-genome-sequencing of triple negative breast cancers in a population-based clinical study , 2019, Nature Medicine.

[18]  E. Sala,et al.  Computed Tomography-Derived Radiomic Metrics Can Identify Responders to Immunotherapy in Ovarian Cancer. , 2019, JCO precision oncology.

[19]  Paul J. Hoffman,et al.  Comprehensive Integration of Single-Cell Data , 2018, Cell.

[20]  Ali Bashashati,et al.  Integrated structural variation and point mutation signatures in cancer genomes using correlated topic models , 2019, PLoS Comput. Biol..

[21]  Kieran R. Campbell,et al.  Probabilistic cell type assignment of single-cell transcriptomic data reveals spatiotemporal microenvironment dynamics in human cancers , 2019, bioRxiv.

[22]  G. Freeman,et al.  PARP Inhibition Elicits STING-Dependent Antitumor Immunity in Brca1-Deficient Ovarian Cancer , 2018, Cell reports.

[23]  Fan Zhang,et al.  Fast, sensitive, and accurate integration of single cell data with Harmony , 2018, bioRxiv.

[24]  L. Cantley,et al.  The Multifaceted Role of Chromosomal Instability in Cancer and Its Microenvironment , 2018, Cell.

[25]  Gregory M. Chen,et al.  Consensus on Molecular Subtypes of High-Grade Serous Ovarian Carcinoma , 2018, Clinical Cancer Research.

[26]  Eugene W. Myers,et al.  Cell Detection with Star-convex Polygons , 2018, MICCAI.

[27]  Wyeth W. Wasserman,et al.  Interfaces of Malignant and Immunologic Clonal Dynamics in Ovarian Cancer , 2018, Cell.

[28]  Ville Mustonen,et al.  The repertoire of mutational signatures in human cancer , 2018, Nature.

[29]  L. Vermeulen,et al.  Turning Cold Tumors Hot by Blocking TGF-β. , 2018, Trends in cancer.

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

[31]  Camille Stephan-Otto Attolini,et al.  TGFβ drives immune evasion in genetically reconstituted colon cancer metastasis , 2018, Nature.

[32]  R. Bourgon,et al.  TGF-β attenuates tumour response to PD-L1 blockade by contributing to exclusion of T cells , 2018, Nature.

[33]  Angela E. Leek,et al.  Allele-Specific HLA Loss and Immune Escape in Lung Cancer Evolution , 2017, Cell.

[34]  R. Scharpf,et al.  High grade serous ovarian carcinomas originate in the fallopian tube , 2017, Nature Communications.

[35]  L. Cope,et al.  Molecular analysis of high-grade serous ovarian carcinoma with and without associated serous tubal intra-epithelial carcinoma , 2017, Nature Communications.

[36]  Oliver Hofmann,et al.  Copy-number signatures and mutational processes in ovarian carcinoma , 2017, Nature Genetics.

[37]  Benjamin J. Raphael,et al.  The evolutionary history of 2,658 cancers , 2017, Nature.

[38]  D. Bajorin,et al.  Global Cancer Transcriptome Quantifies Repeat Element Polarization between Immunotherapy Responsive and T Cell Suppressive Classes , 2017, bioRxiv.

[39]  Ali Bashashati,et al.  Genomic consequences of aberrant DNA repair mechanisms stratify ovarian cancer histotypes , 2017, Nature Genetics.

[40]  Donavan T. Cheng,et al.  Mutational Landscape of Metastatic Cancer Revealed from Prospective Clinical Sequencing of 10,000 Patients , 2017, Nature Medicine.

[41]  Sohrab P. Shah,et al.  deStruct: Accurate Rearrangement Detection using Breakpoint Specific Realignment , 2017, bioRxiv.

[42]  S. Elledge,et al.  Tumor aneuploidy correlates with markers of immune evasion and with reduced response to immunotherapy , 2017, Science.

[43]  Peter Bankhead,et al.  QuPath: Open source software for digital pathology image analysis , 2017, Scientific Reports.

[44]  V. Seshan,et al.  FACETS: allele-specific copy number and clonal heterogeneity analysis tool for high-throughput DNA sequencing , 2016, Nucleic acids research.

[45]  David C. Jones,et al.  Landscape of somatic mutations in 560 breast cancer whole genome sequences , 2016, Nature.

[46]  Ali Bashashati,et al.  Divergent modes of clonal spread and intraperitoneal mixing in high-grade serous ovarian cancer , 2016, Nature Genetics.

[47]  Charles H. Yoon,et al.  Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq , 2016, Science.

[48]  D. Levine,et al.  Invasion Patterns of Metastatic Extrauterine High-grade Serous Carcinoma With BRCA Germline Mutation and Correlation With Clinical Outcomes , 2016, The American journal of surgical pathology.

[49]  G. Shapiro,et al.  Homologous Recombination Deficiency: Exploiting the Fundamental Vulnerability of Ovarian Cancer. , 2015, Cancer discovery.

[50]  K. Cibulskis,et al.  Comprehensive analysis of cancer-associated somatic mutations in class I HLA genes , 2015, Nature Biotechnology.

[51]  Donavan T. Cheng,et al.  Memorial Sloan Kettering-Integrated Mutation Profiling of Actionable Cancer Targets (MSK-IMPACT): A Hybridization Capture-Based Next-Generation Sequencing Clinical Assay for Solid Tumor Molecular Oncology. , 2015, The Journal of molecular diagnostics : JMD.

[52]  Ryan M. Layer,et al.  SpeedSeq: Ultra-fast personal genome analysis and interpretation , 2014, Nature Methods.

[53]  S. Bohlson,et al.  Complement, C1q, and C1q-Related Molecules Regulate Macrophage Polarization , 2014, Front. Immunol..

[54]  Lovelace J. Luquette,et al.  Diverse Mechanisms of Somatic Structural Variations in Human Cancer Genomes , 2013, Cell.

[55]  R. Schreiber,et al.  Persistent LCMV Infection Is Controlled by Blockade of Type I Interferon Signaling , 2013, Science.

[56]  K. Cibulskis,et al.  Prognostically relevant gene signatures of high-grade serous ovarian carcinoma. , 2012, The Journal of clinical investigation.

[57]  Ryan M. Layer,et al.  LUMPY: a probabilistic framework for structural variant discovery , 2012, Genome Biology.

[58]  Wendy S. W. Wong,et al.  Strelka: accurate somatic small-variant calling from sequenced tumor-normal sample pairs , 2012, Bioinform..

[59]  Narciso Olvera,et al.  Morphologic patterns associated with BRCA1 and BRCA2 genotype in ovarian carcinoma , 2012, Modern Pathology.

[60]  Benjamin J. Raphael,et al.  Integrated Genomic Analyses of Ovarian Carcinoma , 2011, Nature.

[61]  Süleyman Cenk Sahinalp,et al.  deFuse: An Algorithm for Gene Fusion Discovery in Tumor RNA-Seq Data , 2011, PLoS Comput. Biol..

[62]  T. Okazaki,et al.  IFN-α Directly Promotes Programmed Cell Death-1 Transcription and Limits the Duration of T Cell-Mediated Immunity , 2011, The Journal of Immunology.

[63]  Carlos Caldas,et al.  Driver mutations in TP53 are ubiquitous in high grade serous carcinoma of the ovary , 2010, The Journal of pathology.

[64]  R. Tothill,et al.  Novel Molecular Subtypes of Serous and Endometrioid Ovarian Cancer Linked to Clinical Outcome , 2008, Clinical Cancer Research.

[65]  P. Angel,et al.  S100A8 and S100A9 in inflammation and cancer. , 2006, Biochemical pharmacology.

[66]  D. Zamarin,et al.  Preparation of single cells from tumors for single-cell RNA sequencing. , 2020, Methods in enzymology.