Lineage tracking reveals dynamic relationships of T cells in colorectal cancer

T cells are key elements of cancer immunotherapy1 but certain fundamental properties, such as the development and migration of T cells within tumours, remain unknown. The enormous T cell receptor (TCR) repertoire, which is required for the recognition of foreign and self-antigens2, could serve as lineage tags to track these T cells in tumours3. Here we obtained transcriptomes of 11,138 single T cells from 12 patients with colorectal cancer, and developed single T cell analysis by RNA sequencing and TCR tracking (STARTRAC) indices to quantitatively analyse the dynamic relationships among 20 identified T cell subsets with distinct functions and clonalities. Although both CD8+ effector and ‘exhausted’ T cells exhibited high clonal expansion, they were independently connected with tumour-resident CD8+ effector memory cells, implicating a TCR-based fate decision. Of the CD4+ T cells, most tumour-infiltrating T regulatory (Treg) cells showed clonal exclusivity, whereas certain Treg cell clones were developmentally linked to several T helper (TH) cell clones. Notably, we identified two IFNG+ TH1-like cell clusters in tumours that were associated with distinct IFNγ-regulating transcription factors —the GZMK+ effector memory T cells, which were associated with EOMES and RUNX3, and CXCL13+BHLHE40+ TH1-like cell clusters, which were associated with BHLHE40. Only CXCL13+BHLHE40+ TH1-like cells were preferentially enriched in patients with microsatellite-instable tumours, and this might explain their favourable responses to immune-checkpoint blockade. Furthermore, IGFLR1 was highly expressed in both CXCL13+BHLHE40+ TH1-like cells and CD8+ exhausted T cells and possessed co-stimulatory functions. Our integrated STARTRAC analyses provide a powerful approach to dissect the T cell properties in colorectal cancer comprehensively, and could provide insights into the dynamic relationships of T cells in other cancers.An integrated RNA-sequencing approach demonstrates that CXCL13+ TH1-like cells are preferentially enriched in microsatellite-instable tumours from patients with colorectal cancer, and IGFLR1 is identified as a co-stimulatory molecule.

[1]  Aviv Regev,et al.  Induction and transcriptional regulation of the co-inhibitory gene module in T cells , 2018, Nature.

[2]  Chih-Chung Lin,et al.  Bhlhe40 is an essential repressor of IL-10 during Mycobacterium tuberculosis infection , 2018, The Journal of experimental medicine.

[3]  Liza Konnikova,et al.  Individual intestinal symbionts induce a distinct population of RORγ+ regulatory T cells , 2015, Science.

[4]  Bert Vogelstein,et al.  PD-1 Blockade in Tumors with Mismatch-Repair Deficiency. , 2015, The New England journal of medicine.

[5]  J. Mesirov,et al.  The Molecular Signatures Database Hallmark Gene Set Collection , 2015 .

[6]  Anjana Rao,et al.  Transcription factors T-bet and Runx3 cooperate to activate Ifng and silence Il4 in T helper type 1 cells , 2007, Nature Immunology.

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

[8]  Alessandro Laio,et al.  Clustering by fast search and find of density peaks , 2014, Science.

[9]  Keji Zhao,et al.  The transcription factor Bhlhe40 is a switch of inflammatory versus antiinflammatory Th1 cell fate determination , 2017, The Journal of experimental medicine.

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

[11]  Jeffrey J Delrow,et al.  Tumor-Specific T Cell Dysfunction Is a Dynamic Antigen-Driven Differentiation Program Initiated Early during Tumorigenesis. , 2016, Immunity.

[12]  Alessandro Sette,et al.  Identifying specificity groups in the T cell receptor repertoire , 2017, Nature.

[13]  Mark M Davis,et al.  Linking T-cell receptor sequence to functional phenotype at the single-cell level , 2014, Nature Biotechnology.

[14]  Hiroyoshi Nishikawa,et al.  Natural and Induced T Regulatory Cells in Cancer , 2013, Front. Immunol..

[15]  Burton E. Barnett,et al.  Progenitor and Terminal Subsets of CD8+ T Cells Cooperate to Contain Chronic Viral Infection , 2012, Science.

[16]  H. Hakonarson,et al.  ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data , 2010, Nucleic acids research.

[17]  Aviv Regev,et al.  A Distinct Gene Module for Dysfunction Uncoupled from Activation in Tumor-Infiltrating T Cells , 2016, Cell.

[18]  Hilde Cheroutre,et al.  The light and dark sides of intestinal intraepithelial lymphocytes , 2011, Nature Reviews Immunology.

[19]  J. Marioni,et al.  Pooling across cells to normalize single-cell RNA sequencing data with many zero counts , 2016, Genome Biology.

[20]  S. Halgamuge,et al.  Inferring copy number and genotype in tumour exome data , 2014, BMC Genomics.

[21]  K. Kinzler,et al.  Cancer Genome Landscapes , 2013, Science.

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

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

[24]  Nargis Khan,et al.  Triggering through Toll-like receptor 2 limits chronically stimulated T-helper type 1 cells from undergoing exhaustion. , 2015, The Journal of infectious diseases.

[25]  Jovana Maksimovic,et al.  Genome-wide DNA methylation analysis identifies hypomethylated genes regulated by FOXP3 in human regulatory T cells. , 2013, Blood.

[26]  Liming Lu,et al.  Metabolic control of regulatory T cell stability and function by TRAF3IP3 at the lysosome , 2018, The Journal of experimental medicine.

[27]  J. Schenkel,et al.  Tissue-resident memory T cells. , 2014, Immunity.

[28]  R. Balderas,et al.  Repression of the genome organizer SATB1 in regulatory T cells is required for suppressive function and inhibition of effector differentiation , 2011, Nature Immunology.

[29]  A. Luster,et al.  Chemokines and chemokine receptors: positioning cells for host defense and immunity. , 2014, Annual review of immunology.

[30]  D. Godfrey,et al.  Raising the NKT cell family , 2010, Nature Immunology.

[31]  Richard Durbin,et al.  Sequence analysis Fast and accurate short read alignment with Burrows – Wheeler transform , 2009 .

[32]  Joydeep Ghosh,et al.  Cluster Ensembles --- A Knowledge Reuse Framework for Combining Multiple Partitions , 2002, J. Mach. Learn. Res..

[33]  Edward C Hutchinson,et al.  MAIT cells are activated during human viral infections , 2016, Nature Communications.

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

[35]  I. Mellman,et al.  Elements of cancer immunity and the cancer–immune set point , 2017, Nature.

[36]  Serban Nacu,et al.  Fast and SNP-tolerant detection of complex variants and splicing in short reads , 2010, Bioinform..

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

[38]  Wenjun Ouyang,et al.  Murine Insulin Growth Factor-like (IGFL) and Human IGFL1 Proteins Are Induced in Inflammatory Skin Conditions and Bind to a Novel Tumor Necrosis Factor Receptor Family Member, IGFLR1 , 2011, The Journal of Biological Chemistry.

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

[40]  M. Vignali,et al.  T‐cell receptor profiling in cancer , 2015, Molecular oncology.

[41]  M. Schaub,et al.  SC3 - consensus clustering of single-cell RNA-Seq data , 2016, Nature Methods.

[42]  Richard P. Junghans,et al.  Second-Generation Anti–Carcinoembryonic Antigen Designer T Cells Resist Activation-Induced Cell Death, Proliferate on Tumor Contact, Secrete Cytokines, and Exhibit Superior Antitumor Activity In vivo: A Preclinical Evaluation , 2008, Clinical Cancer Research.

[43]  C. Klein,et al.  A transcriptionally and functionally distinct PD-1+ CD8+ T cell pool with predictive potential in non-small cell lung cancer treated with PD-1 blockade , 2018, Nature Medicine.

[44]  Steven J. M. Jones,et al.  Comprehensive molecular characterization of human colon and rectal cancer , 2012, Nature.

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

[46]  Ludmila V. Danilova,et al.  Mismatch repair deficiency predicts response of solid tumors to PD-1 blockade , 2017, Science.

[47]  Boxi Kang,et al.  Global characterization of T cells in non-small-cell lung cancer by single-cell sequencing , 2018, Nature Medicine.

[48]  Gavin D. Grant,et al.  Common markers of proliferation , 2006, Nature Reviews Cancer.

[49]  Andrew J. Hill,et al.  Single-cell mRNA quantification and differential analysis with Census , 2017, Nature Methods.

[50]  Dana Pe’er,et al.  Distinct Cellular Mechanisms Underlie Anti-CTLA-4 and Anti-PD-1 Checkpoint Blockade , 2017, Cell.

[51]  Z. Trajanoski,et al.  Integrative Analyses of Colorectal Cancer Show Immunoscore Is a Stronger Predictor of Patient Survival Than Microsatellite Instability. , 2016, Immunity.