Opposing immune and genetic forces shape oncogenic programs in synovial sarcoma

Synovial sarcoma is an aggressive mesenchymal neoplasm, driven by the SS18-SSX fusion, and characterized by immunogenic antigens expression and exceptionally low T cell infiltration levels. To study the cancer-immune interplay in this disease, we profiled 16,872 cells from 12 human synovial sarcoma tumors using single-cell RNA-sequencing (scRNA-Seq). Synovial sarcoma manifests antitumor immunity, high cellular plasticity and a core oncogenic program, which is predictive of low immune levels and poor clinical outcomes. Using genetic and pharmacological perturbations, we demonstrate that the program is controlled by the SS18-SSX driver and repressed by cytokines secreted by macrophages and T cells in the tumor microenvironment. Network modeling predicted that SS18-SSX promotes the program through HDAC1 and CDK6. Indeed, the combination of HDAC and CDK4/6 inhibitors represses the program, induces immunogenic cell states, and selectively targets synovial sarcoma cells. Our study demonstrates that immune evasion, cellular plasticity, and cell cycle are co-regulated and can be co-targeted in synovial sarcoma and potentially in other malignancies.

[1]  Gordon B. Mills,et al.  High multiplex, digital spatial profiling of proteins and RNA in fixed tissue using genomic detection methods , 2019, bioRxiv.

[2]  Evan Z. Macosko,et al.  Integrative inference of brain cell similarities and differences from single-cell genomics , 2018, bioRxiv.

[3]  Monika S. Kowalczyk,et al.  A Cancer Cell Program Promotes T Cell Exclusion and Resistance to Checkpoint Blockade , 2018, Cell.

[4]  Paul J. Hoover,et al.  Defining T Cell States Associated with Response to Checkpoint Immunotherapy in Melanoma , 2018, Cell.

[5]  Karoly Szuhai,et al.  PRAME and HLA Class I expression patterns make synovial sarcoma a suitable target for PRAME specific T-cell receptor gene therapy , 2018, Oncoimmunology.

[6]  J. Luke,et al.  T Cell–Inflamed versus Non-T Cell–Inflamed Tumors: A Conceptual Framework for Cancer Immunotherapy Drug Development and Combination Therapy Selection , 2018, Cancer Immunology Research.

[7]  Erik Sundström,et al.  RNA velocity of single cells , 2018, Nature.

[8]  P. A. Futreal,et al.  Mutations in the SWI/SNF complex induce a targetable dependence on oxidative phosphorylation in lung cancer , 2018, Nature Medicine.

[9]  D. Templeton,et al.  Macrophage Migration Inhibitory Factor protects cancer cells from immunogenic cell death and impairs anti-tumor immune responses , 2018, PloS one.

[10]  Hannah C. Beird,et al.  The SS18-SSX Fusion Oncoprotein Hijacks BAF Complex Targeting and Function to Drive Synovial Sarcoma. , 2018, Cancer cell.

[11]  Sandro Santagata,et al.  Highly multiplexed immunofluorescence imaging of human tissues and tumors using t-CyCIF and conventional optical microscopes , 2018, eLife.

[12]  Tracy T Batchelor,et al.  Developmental and oncogenic programs in H3K27M gliomas dissected by single-cell RNA-seq , 2018, Science.

[13]  C. Garbutt,et al.  Inhibition of cyclin-dependent kinase 4 as a potential therapeutic strategy for treatment of synovial sarcoma , 2018, Cell Death & Disease.

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

[15]  Darjus F. Tschaharganeh,et al.  The SS18-SSX Oncoprotein Hijacks KDM2B-PRC1.1 to Drive Synovial Sarcoma. , 2018, Cancer cell.

[16]  Mark W. Ball,et al.  Genomic correlates of response to immune checkpoint therapies in clear cell renal cell carcinoma , 2018, Science.

[17]  Henry W. Long,et al.  A major chromatin regulator determines resistance of tumor cells to T cell–mediated killing , 2018, Science.

[18]  S. Pollack.  The potential of the CMB305 vaccine regimen to target NY-ESO-1 and improve outcomes for synovial sarcoma and myxoid/round cell liposarcoma patients , 2017, Expert review of vaccines.

[19]  David J. Arenillas,et al.  JASPAR 2018: update of the open-access database of transcription factor binding profiles and its web framework , 2017, Nucleic acids research.

[20]  Steven J. M. Jones,et al.  Comprehensive and Integrated Genomic Characterization of Adult Soft Tissue Sarcomas , 2017, Cell.

[21]  Shawn M. Gillespie,et al.  Single-Cell Transcriptomic Analysis of Primary and Metastatic Tumor Ecosystems in Head and Neck Cancer , 2017, Cell.

[22]  Boxi Kang,et al.  Landscape of Infiltrating T Cells in Liver Cancer Revealed by Single-Cell Sequencing , 2017, Cell.

[23]  P. Hwu,et al.  Overexpressed PRAME is a potential immunotherapy target in sarcoma subtypes , 2017, Clinical Sarcoma Research.

[24]  Shugang Li,et al.  Evaluation of expression of cancer stem cell markers and fusion gene in synovial sarcoma: Insights into histogenesis and pathogenesis. , 2017, Oncology reports.

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

[26]  Robin L. Jones,et al.  T‐cell infiltration and clonality correlate with programmed cell death protein 1 and programmed death‐ligand 1 expression in patients with soft tissue sarcomas , 2017, Cancer.

[27]  Timothy L. Tickle,et al.  STAR-Fusion: Fast and Accurate Fusion Transcript Detection from RNA-Seq , 2017, bioRxiv.

[28]  K. Harimaya,et al.  Cancer-testis antigen expression in synovial sarcoma: NY-ESO-1, PRAME, MAGEA4, and MAGEA1. , 2017, Human pathology.

[29]  P. Tamayo,et al.  Abstract S1-01: Whole exome and transcriptome sequencing of resistant ER+ metastatic breast cancer , 2017 .

[30]  Mariella G. Filbin,et al.  Decoupling genetics, lineages, and microenvironment in IDH-mutant gliomas by single-cell RNA-seq , 2017, Science.

[31]  Kevin B. Jones,et al.  HDAC and Proteasome Inhibitors Synergize to Activate Pro-Apoptotic Factors in Synovial Sarcoma , 2017, PloS one.

[32]  Ashley V. Menk,et al.  Efficacy of PD-1 Blockade Is Potentiated by Metformin-Induced Reduction of Tumor Hypoxia , 2016, Cancer Immunology Research.

[33]  Mariella G. Filbin,et al.  Single-cell RNA-seq supports a developmental hierarchy in human oligodendroglioma , 2016, Nature.

[34]  Grace X. Y. Zheng,et al.  Massively parallel digital transcriptional profiling of single cells , 2016, Nature Communications.

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

[36]  Anneliese O. Speak,et al.  T cell fate and clonality inference from single cell transcriptomes , 2016, Nature Methods.

[37]  Joseph L. Herman,et al.  Characterizing transcriptional heterogeneity through pathway and gene set overdispersion analysis , 2015, Nature Methods.

[38]  George Michailidis,et al.  A non-negative matrix factorization method for detecting modules in heterogeneous omics multi-modal data , 2015, Bioinform..

[39]  M. Mann,et al.  Functional classification of memory CD8+ T cells by CX3CR1 expression , 2015, Nature Communications.

[40]  Fabian J. Theis,et al.  Diffusion maps for high-dimensional single-cell analysis of differentiation data , 2015, Bioinform..

[41]  Piero Carninci,et al.  A draft network of ligand–receptor-mediated multicellular signalling in human , 2015, Nature Communications.

[42]  Steven J. M. Jones,et al.  Genomic Classification of Cutaneous Melanoma , 2015, Cell.

[43]  Jung Eun Shim,et al.  TRRUST: a reference database of human transcriptional regulatory interactions , 2015, Scientific Reports.

[44]  T. Gajewski,et al.  Melanoma-intrinsic β-catenin signalling prevents anti-tumour immunity , 2015, Nature.

[45]  Trevor J Pugh,et al.  Oncotator: Cancer Variant Annotation Tool , 2015, Human mutation.

[46]  Marc Ladanyi,et al.  Synovial sarcoma: recent discoveries as a roadmap to new avenues for therapy. , 2015, Cancer discovery.

[47]  Stein Aerts,et al.  iRegulon: From a Gene List to a Gene Regulatory Network Using Large Motif and Track Collections , 2014, PLoS Comput. Biol..

[48]  Shawn M. Gillespie,et al.  Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma , 2014, Science.

[49]  Adam Kiezun,et al.  Whole-exome sequencing and clinical interpretation of FFPE tumor samples to guide precision cancer medicine , 2013, Nature Medicine.

[50]  Åsa K. Björklund,et al.  Full-length RNA-seq from single cells using Smart-seq2 , 2014, Nature Protocols.

[51]  Ludo Waltman,et al.  A smart local moving algorithm for large-scale modularity-based community detection , 2013, The European Physical Journal B.

[52]  G. Crabtree,et al.  Reversible Disruption of mSWI/SNF (BAF) Complexes by the SS18-SSX Oncogenic Fusion in Synovial Sarcoma , 2013, Cell.

[53]  Olivier Delattre,et al.  Chromosome instability accounts for reverse metastatic outcomes of pediatric and adult synovial sarcomas. , 2013, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[54]  A. Sivachenko,et al.  Sensitive detection of somatic point mutations in impure and heterogeneous cancer samples , 2013, Nature Biotechnology.

[55]  Trevor J Pugh,et al.  Discovery and characterization of artifactual mutations in deep coverage targeted capture sequencing data due to oxidative DNA damage during sample preparation , 2013, Nucleic acids research.

[56]  Andrew McDavid,et al.  Data exploration, quality control and testing in single-cell qPCR-based gene expression experiments , 2012, Bioinform..

[57]  Thomas R. Gingeras,et al.  STAR: ultrafast universal RNA-seq aligner , 2013, Bioinform..

[58]  Klemens Vierlinger,et al.  Meta-Analysis of Gene Expression Signatures Defining the Epithelial to Mesenchymal Transition during Cancer Progression , 2012, PloS one.

[59]  S. Aerts,et al.  i-cisTarget: an integrative genomics method for the prediction of regulatory features and cis-regulatory modules , 2012, Nucleic acids research.

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

[61]  C. Sautès-Fridman,et al.  The immune contexture in human tumours: impact on clinical outcome , 2012, Nature Reviews Cancer.

[62]  Kevin B. Jones,et al.  Deconstruction of the SS18-SSX fusion oncoprotein complex: insights into disease etiology and therapeutics. , 2012, Cancer cell.

[63]  M. Acencio,et al.  HTRIdb: an open-access database for experimentally verified human transcriptional regulation interactions , 2012, BMC Genomics.

[64]  D. McNeel,et al.  Vaccines Targeting the Cancer-testis Antigen SSX-2 Elicit HLA-A2 Epitope-specific Cytolytic T Cells , 2011, Journal of immunotherapy.

[65]  Kristian Cibulskis,et al.  ContEst: estimating cross-contamination of human samples in next-generation sequencing data , 2011, Bioinform..

[66]  Colin N. Dewey,et al.  RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome , 2011, BMC Bioinformatics.

[67]  Helga Thorvaldsdóttir,et al.  Molecular signatures database (MSigDB) 3.0 , 2011, Bioinform..

[68]  M. DePristo,et al.  A framework for variation discovery and genotyping using next-generation DNA sequencing data , 2011, Nature Genetics.

[69]  Dennis C. Friedrich,et al.  A scalable, fully automated process for construction of sequence-ready human exome targeted capture libraries , 2011, Genome Biology.

[70]  Allan Kuchinsky,et al.  GLay: community structure analysis of biological networks , 2010, Bioinform..

[71]  Teruhiko Yoshida,et al.  Gene Expression Profiling of Synovial Sarcoma: Distinct Signature of Poorly Differentiated Type , 2010, The American journal of surgical pathology.

[72]  D. McNeel,et al.  The SSX Family of Cancer-Testis Antigens as Target Proteins for Tumor Therapy , 2010, Clinical & developmental immunology.

[73]  Kakajan Komurov,et al.  Core epithelial-to-mesenchymal transition interactome gene-expression signature is associated with claudin-low and metaplastic breast cancer subtypes , 2010, Proceedings of the National Academy of Sciences.

[74]  H. Yoshikawa,et al.  Synovial Sarcoma Is a Stem Cell Malignancy , 2010, Stem cells.

[75]  Nicolò Riggi,et al.  EWS-FLI-1 modulates miRNA145 and SOX2 expression to initiate mesenchymal stem cell reprogramming toward Ewing sarcoma cancer stem cells. , 2010, Genes & development.

[76]  Cole Trapnell,et al.  Ultrafast and memory-efficient alignment of short DNA sequences to the human genome , 2009, Genome Biology.

[77]  Alberto Mantovani,et al.  Transcriptional Profiling of the Human Monocyte-to-Macrophage Differentiation and Polarization: New Molecules and Patterns of Gene Expression1 , 2006, The Journal of Immunology.

[78]  H. Inoue,et al.  Significant growth suppression of synovial sarcomas by the histone deacetylase inhibitor FK228 in vitro and in vivo. , 2005, Cancer letters.

[79]  S. Stevanović,et al.  CD4+ T Cell Responses to SSX-4 in Melanoma Patients1 , 2005, The Journal of Immunology.

[80]  K. Lindblad-Toh,et al.  Systematic discovery of regulatory motifs in human promoters and 3′ UTRs by comparison of several mammals , 2005, Nature.

[81]  M. Wigler,et al.  Circular binary segmentation for the analysis of array-based DNA copy number data. , 2004, Biostatistics.

[82]  M. Newman,et al.  Finding and evaluating community structure in networks. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[83]  D. Speiser,et al.  Tumor-reactive, SSX-2-specific CD8+ T cells are selectively expanded during immune responses to antigen-expressing tumors in melanoma patients. , 2003, Cancer research.

[84]  Sergey Brin,et al.  The Anatomy of a Large-Scale Hypertextual Web Search Engine , 1998, Comput. Networks.

[85]  Y. Benjamini,et al.  Controlling the false discovery rate: a practical and powerful approach to multiple testing , 1995 .

[86]  M. Carlson,et al.  Genes affecting the regulation of SUC2 gene expression by glucose repression in Saccharomyces cerevisiae. , 1984, Genetics.