Early transcriptional and epigenetic regulation of CD8+ T cell differentiation revealed by single-cell RNA-seq

During microbial infection, responding CD8+ T lymphocytes differentiate into heterogeneous subsets that together provide immediate and durable protection. To elucidate the dynamic transcriptional changes that underlie this process, we applied a single-cell RNA-sequencing approach and analyzed individual CD8+ T lymphocytes sequentially throughout the course of a viral infection in vivo. Our analyses revealed a striking transcriptional divergence among cells that had undergone their first division and identified previously unknown molecular determinants that controlled the fate specification of CD8+ T lymphocytes. Our findings suggest a model for the differentiation of terminal effector cells initiated by an early burst of transcriptional activity and subsequently refined by epigenetic silencing of transcripts associated with memory lymphocytes, which highlights the power and necessity of single-cell approaches.

[1]  C. Glass,et al.  Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. , 2010, Molecular cell.

[2]  R. Wersto,et al.  Histone Acetylation Facilitates Rapid and Robust Memory CD8 T Cell Response through Differential Expression of Effector Molecules (Eomesodermin and Its Targets: Perforin and Granzyme B)1 , 2008, The Journal of Immunology.

[3]  Andrew J. McMichael,et al.  Molecular Signatures Distinguish Human Central Memory from Effector Memory CD8 T Cell Subsets1 , 2005, The Journal of Immunology.

[4]  D. Pellicci,et al.  T-box Transcription Factors Combine with the Cytokines TGF-β and IL-15 to Control Tissue-Resident Memory T Cell Fate. , 2015, Immunity.

[5]  John T. Chang,et al.  Molecular regulation of effector and memory T cell differentiation , 2014, Nature Immunology.

[6]  Rona S. Gertner,et al.  Single-Cell Genomics Unveils Critical Regulators of Th17 Cell Pathogenicity , 2015, Cell.

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

[8]  Sean C. Bendall,et al.  viSNE enables visualization of high dimensional single-cell data and reveals phenotypic heterogeneity of leukemia , 2013, Nature Biotechnology.

[9]  J. Bluestone,et al.  The chromatin-modifying enzyme Ezh2 is critical for the maintenance of regulatory T cell identity after activation. , 2015, Immunity.

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

[11]  M. Ffrench,et al.  BTG1, a member of a new family of antiproliferative genes. , 1992, The EMBO journal.

[12]  E. Yang,et al.  Transcriptional insights into the CD8+ T cell response to infection and memory T cell formation , 2013, Nature Immunology.

[13]  R. Emerson,et al.  Common clonal origin of central and resident memory T cells following skin immunization , 2015, Nature Medicine.

[14]  J. Powell,et al.  Cytosolic Branched Chain Aminotransferase (BCATc) Regulates mTORC1 Signaling and Glycolytic Metabolism in CD4+ T Cells* , 2014, The Journal of Biological Chemistry.

[15]  Sean C. Bendall,et al.  Wishbone identifies bifurcating developmental trajectories from single-cell data , 2016, Nature Biotechnology.

[16]  Gaël Varoquaux,et al.  Scikit-learn: Machine Learning in Python , 2011, J. Mach. Learn. Res..

[17]  E. Zúñiga,et al.  Cell-intrinsic transforming growth factor-beta signaling mediates virus-specific CD8+ T cell deletion and viral persistence in vivo. , 2009, Immunity.

[18]  Sean C. Bendall,et al.  Single-Cell Mass Cytometry of Differential Immune and Drug Responses Across a Human Hematopoietic Continuum , 2011, Science.

[19]  Rustom Antia,et al.  Lineage relationship and protective immunity of memory CD8 T cell subsets , 2003, Nature Immunology.

[20]  S. Jameson,et al.  Transcriptional downregulation of S1pr1 is required for establishment of resident memory CD8+ T cells , 2013, Nature Immunology.

[21]  N. Neff,et al.  Reconstructing lineage hierarchies of the distal lung epithelium using single cell RNA-seq , 2014, Nature.

[22]  T. Speed,et al.  Distinct epigenetic signatures delineate transcriptional programs during virus-specific CD8(+) T cell differentiation. , 2014, Immunity.

[23]  Yutaka Suzuki,et al.  The polycomb protein Ezh2 regulates differentiation and plasticity of CD4(+) T helper type 1 and type 2 cells. , 2013, Immunity.

[24]  W. Huber,et al.  Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 , 2014, Genome Biology.

[25]  Lior Pachter,et al.  Near-optimal probabilistic RNA-seq quantification , 2016, Nature Biotechnology.

[26]  Greg M. Delgoffe,et al.  Asymmetric inheritance of mTORC1 kinase activity during division dictates CD8 T cell differentiation , 2016, Nature Immunology.

[27]  Michael Poidinger,et al.  High-dimensional analysis of the murine myeloid cell system , 2014, Nature Immunology.

[28]  Lior Pachter,et al.  Fast and accurate single-cell RNA-seq analysis by clustering of transcript-compatibility counts , 2016, Genome Biology.

[29]  O. Kretz,et al.  Mitochondrial Dynamics Controls T Cell Fate through Metabolic Programming , 2016, Cell.

[30]  Nu Zhang,et al.  Transforming growth factor-β signaling is constantly shaping memory T-cell population , 2015, Proceedings of the National Academy of Sciences.

[31]  Susan M. Kaech,et al.  Molecular and Functional Profiling of Memory CD8 T Cell Differentiation , 2002, Cell.

[32]  R. Klose,et al.  Targeting Polycomb systems to regulate gene expression: modifications to a complex story , 2015, Nature Reviews Molecular Cell Biology.

[33]  Mouse Genome Sequencing Consortium.  Initial sequencing and comparative analysis of the mouse genome , 2002, Nature.

[34]  Yuka Kanno,et al.  BACH2 regulates CD8+ T cell differentiation by controlling access of AP-1 factors to enhancers , 2016, Nature Immunology.

[35]  S. Varambally,et al.  Cancer mediates effector T cell dysfunction by targeting microRNAs and EZH2 via glycolysis restriction , 2015, Nature Immunology.

[36]  Hendrik Blockeel,et al.  Introduction to the special issue on multi-relational data mining and statistical relational learning , 2006, Machine Learning.

[37]  D. Green,et al.  Metabolic Maintenance of Cell Asymmetry following Division in Activated T Lymphocytes , 2016, Nature.

[38]  P. Kourilsky,et al.  Lineage relationships, homeostasis, and recall capacities of central– and effector–memory CD8 T cells in vivo , 2005, The Journal of experimental medicine.

[39]  John T. Chang,et al.  Asymmetric proteasome segregation as a mechanism for unequal partitioning of the transcription factor T-bet during T lymphocyte division. , 2011, Immunity.

[40]  Scott N. Mueller,et al.  Tissue-resident memory T cells: local specialists in immune defence , 2015, Nature Reviews Immunology.

[41]  John T. Chang,et al.  Asymmetric T Lymphocyte Division in the Initiation of Adaptive Immune Responses , 2007, Science.

[42]  U. Klein,et al.  Asymmetric PI3K Signaling Driving Developmental and Regenerative Cell Fate Bifurcation. , 2015, Cell reports.

[43]  A. Quinlan.  BEDTools: The Swiss‐Army Tool for Genome Feature Analysis , 2014, Current protocols in bioinformatics.

[44]  J. Lieberman,et al.  Effector differentiation is not prerequisite for generation of memory cytotoxic T lymphocytes. , 2001, The Journal of clinical investigation.

[45]  Mark D. Robinson,et al.  edgeR: a Bioconductor package for differential expression analysis of digital gene expression data , 2009, Bioinform..

[46]  M. Stephens.  EDF Statistics for Goodness of Fit and Some Comparisons , 1974 .

[47]  Christoph Wülfing,et al.  Polycomb Group Protein Ezh2 Controls Actin Polymerization and Cell Signaling , 2005, Cell.

[48]  G. V. D. van der Windt,et al.  Mitochondrial respiratory capacity is a critical regulator of CD8+ T cell memory development. , 2012, Immunity.

[49]  S. Jameson,et al.  Expression of the transcription factor lung Krüppel-like factor is regulated by cytokines and correlates with survival of memory T cells in vitro and in vivo. , 1999, Journal of immunology.

[50]  Yang Cheng,et al.  Categorical Analysis of Human T Cell Heterogeneity with One-Dimensional Soli-Expression by Nonlinear Stochastic Embedding , 2016, The Journal of Immunology.

[51]  Takeshi Yamada,et al.  The transcription factor ELF4 controls proliferation and homing of CD8+ T cells via the Krüppel-like factors KLF4 and KLF2 , 2009, Nature Immunology.

[52]  Nikhil S. Joshi,et al.  Inflammation directs memory precursor and short-lived effector CD8(+) T cell fates via the graded expression of T-bet transcription factor. , 2007, Immunity.

[53]  Cathy H. Wu,et al.  Software for pre-processing Illumina next-generation sequencing short read sequences , 2014, Source Code for Biology and Medicine.

[54]  Rona S. Gertner,et al.  Single-cell transcriptomics reveals bimodality in expression and splicing in immune cells , 2013, Nature.

[55]  Pierre Geurts,et al.  Extremely randomized trees , 2006, Machine Learning.

[56]  T. Speed,et al.  Distinct Epigenetic Signatures Delineate Transcriptional Programs during Virus-Specific CD8 T Cell Differentiation , 2014 .

[57]  Robert A. Edwards,et al.  Quality control and preprocessing of metagenomic datasets , 2011, Bioinform..

[58]  R. Sandberg,et al.  Full-Length mRNA-Seq from single cell levels of RNA and individual circulating tumor cells , 2012, Nature Biotechnology.

[59]  John T. Chang,et al.  Early specification of CD8+ T lymphocyte fates during adaptive immunity revealed by single-cell gene expression analyses , 2014, Nature Immunology.

[60]  J. Harty,et al.  Initial T cell receptor transgenic cell precursor frequency dictates critical aspects of the CD8(+) T cell response to infection. , 2007, Immunity.

[61]  Geoffrey E. Hinton,et al.  Visualizing Data using t-SNE , 2008 .

[62]  Diego G. Silva,et al.  Identification of T Cell-Restricted Genes, and Signatures for Different T Cell Responses, Using a Comprehensive Collection of Microarray Datasets1 , 2005, The Journal of Immunology.

[63]  R. Jaenisch,et al.  DNA Methylation by DNA Methyltransferase 1 Is Critical for Effector CD8 T Cell Expansion1 , 2006, The Journal of Immunology.

[64]  John T. Chang,et al.  Regulation of Asymmetric Division and CD8+ T Lymphocyte Fate Specification by Protein Kinase Cζ and Protein Kinase Cλ/ι , 2015, The Journal of Immunology.

[65]  R. Ahmed,et al.  Chronic virus infection enforces demethylation of the locus that encodes PD-1 in antigen-specific CD8(+) T cells. , 2011, Immunity.

[66]  John T. Chang,et al.  Epigenetic landscapes reveal transcription factors regulating CD8+ T cell differentiation , 2017, Nature Immunology.

[67]  Dustin E. Schones,et al.  Genome-wide analysis of histone methylation reveals chromatin state-based regulation of gene transcription and function of memory CD8+ T cells. , 2009, Immunity.

[68]  Robert Tibshirani,et al.  Finding consistent patterns: A nonparametric approach for identifying differential expression in RNA-Seq data , 2013, Statistical methods in medical research.

[69]  F. Sallusto,et al.  Two subsets of memory T lymphocytes with distinct homing potentials and effector functions , 1999, Nature.