Michael J. T. Stubbington

T cell fate and clonality inference from single-cell transcriptomes

We developed TraCeR, a computational method to reconstruct full-length, paired T cell receptor (TCR) sequences from T lymphocyte single-cell RNA sequence data. TraCeR links T cell specificity with functional response by revealing clonal relationships between cells alongside their transcriptional profiles. We found that T cell clonotypes in a mouse Salmonella infection model span early activated CD4+ T cells as well as mature effector and memory cells.

Science Forum: The Human Cell Atlas

The recent advent of methods for high-throughput single-cell molecular profiling has catalyzed a growing sense in the scientific community that the time is ripe to complete the 150-year-old effort to identify all cell types in the human body. The Human Cell Atlas Project is an international collaborative effort that aims to define all human cell types in terms of distinctive molecular profiles (such as gene expression profiles) and to connect this information with classical cellular descriptions (such as location and morphology). An open comprehensive reference map of the molecular state of cells in healthy human tissues would propel the systematic study of physiological states, developmental trajectories, regulatory circuitry and interactions of cells, and also provide a framework for understanding cellular dysregulation in human disease. Here we describe the idea, its potential utility, early proofs-of-concept, and some design considerations for the Human Cell Atlas, including a commitment to open data, code, and community.

Aging increases cell-to-cell transcriptional variability upon immune stimulation

Single-cell sequencing of mouse immune cells reveals how aging destabilizes a conserved transcriptional activation program. Aging and variability among immune cells How and why the immune system becomes less effective with age are not well understood. Martinez-Jimenez et al. performed single-cell sequencing of CD4+ T cells in old and young mice of two species. In young mice, the gene expression program of early immune activation was tightly regulated and conserved between species. However, as mice aged, the expression of genes involved in pathways responding to immune cell stimulation was not as robust and exhibited increased cell-to-cell variability. Science, this issue p. 1433 Aging is characterized by progressive loss of physiological and cellular functions, but the molecular basis of this decline remains unclear. We explored how aging affects transcriptional dynamics using single-cell RNA sequencing of unstimulated and stimulated naïve and effector memory CD4+ T cells from young and old mice from two divergent species. In young animals, immunological activation drives a conserved transcriptomic switch, resulting in tightly controlled gene expression characterized by a strong up-regulation of a core activation program, coupled with a decrease in cell-to-cell variability. Aging perturbed the activation of this core program and increased expression heterogeneity across populations of cells in both species. These discoveries suggest that increased cell-to-cell transcriptional variability will be a hallmark feature of aging across most, if not all, mammalian tissues.

Single-cell transcriptomics to explore the immune system in health and disease

The immune system varies in cell types, states, and locations. The complex networks, interactions, and responses of immune cells produce diverse cellular ecosystems composed of multiple cell types, accompanied by genetic diversity in antigen receptors. Within this ecosystem, innate and adaptive immune cells maintain and protect tissue function, integrity, and homeostasis upon changes in functional demands and diverse insults. Characterizing this inherent complexity requires studies at single-cell resolution. Recent advances such as massively parallel single-cell RNA sequencing and sophisticated computational methods are catalyzing a revolution in our understanding of immunology. Here we provide an overview of the state of single-cell genomics methods and an outlook on the use of single-cell techniques to decipher the adaptive and innate components of immunity.

An atlas of mouse CD4 + T cell transcriptomes

BackgroundCD4+ T cells are key regulators of the adaptive immune system and can be divided into T helper (Th) cells and regulatory T (Treg) cells. During an immune response Th cells mature from a naive state into one of several effector subtypes that exhibit distinct functions. The transcriptional mechanisms that underlie the specific functional identity of CD4+ T cells are not fully understood.ResultsTo assist investigations into the transcriptional identity and regulatory processes of these cells we performed mRNA-sequencing on three murine T helper subtypes (Th1, Th2 and Th17) as well as on splenic Treg cells and induced Treg (iTreg) cells. Our integrated analysis of this dataset revealed the gene expression changes associated with these related but distinct cellular identities. Each cell subtype differentially expresses a wealth of ‘subtype upregulated’ genes, some of which are well known whilst others promise new insights into signalling processes and transcriptional regulation. We show that hundreds of genes are regulated purely by alternative splicing to extend our knowledge of the role of post-transcriptional regulation in cell differentiation.ConclusionsThis CD4+ transcriptome atlas provides a valuable resource for the study of CD4+ T cell populations. To facilitate its use by others, we have made the data available in an easily accessible online resource at www.th-express.org.ReviewersThis article was reviewed by Wayne Hancock, Christine Wells and Erik van Nimwegen.

Single-cell transcriptome analysis of fish immune cells provides insight into the evolution of vertebrate immune cell types

The immune system of vertebrate species consists of many different cell types that have distinct functional roles and are subject to different evolutionary pressures. Here, we first analyzed conservation of genes specific for all major immune cell types in human and mouse. Our results revealed higher gene turnover and faster evolution of trans-membrane proteins in NK cells compared with other immune cell types, and especially T cells, but similar conservation of nuclear and cytoplasmic protein coding genes. To validate these findings in a distant vertebrate species, we used single-cell RNA sequencing of lck:GFP cells in zebrafish and obtained the first transcriptome of specific immune cell types in a nonmammalian species. Unsupervised clustering and single-cell TCR locus reconstruction identified three cell populations, T cells, a novel type of NK-like cells, and a smaller population of myeloid-like cells. Differential expression analysis uncovered new immune-cell-specific genes, including novel immunoglobulin-like receptors, and neofunctionalization of recently duplicated paralogs. Evolutionary analyses confirmed the higher gene turnover of trans-membrane proteins in NK cells compared with T cells in fish species, suggesting that this is a general property of immune cell types across all vertebrates.

Single-cell analysis at the threshold

1111 Xi Chen, Michael J.T. Stubbington and Sarah A. Teichmann are at the Wellcome Trust Sanger Institute, Cambridge, UK; J. Christopher Love is at the Koch Institute for Integrative Cancer Research, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA; Nicholas E. Navin is in the Department of Genetics and Department of Bioinformatics and Computational Biology, The University of Texas, MD Anderson Cancer Center, Houston, Texas, USA; Lior Pachter is in the Department of Mathematics, University of California, Berkeley, California, USA; Valentine Svensson is at the European Bioinformatics Institute, European Molecular Biology Laboratory, Cambridge, UK; and Jonathan V. Sweedler is in the Department of Chemistry, University of Illinois, Urbana, Illinois, USA. e-mail: lpachter@math.berkeley.edu, sweedler@scs.illinois.edu, clove@mit.edu, nnavin@mdanderson.org or sarah.teichmann@sanger.ac.uk Single-cell analysis at the threshold

The Regulatory T Cell Lineage Factor Foxp3 Regulates Gene Expression through Several Distinct Mechanisms Mostly Independent of Direct DNA Binding.

The lineage factor Foxp3 is essential for the development and maintenance of regulatory T cells, but little is known about the mechanisms involved. Here, we demonstrate that an N-terminal proline-rich interaction region is crucial for Foxp3’s function. Subdomains within this key region link Foxp3 to several independent mechanisms of transcriptional regulation. Our study suggests that Foxp3, even in the absence of its DNA-binding forkhead domain, acts as a bridge between DNA-binding interaction partners and proteins with effector function permitting it to regulate a large number of genes. We show that, in one such mechanism, Foxp3 recruits class I histone deacetylases to the promoters of target genes, counteracting activation-induced histone acetylation and thereby suppressing their expression.

Two Mutually Exclusive Local Chromatin States Drive Efficient V(D)J Recombination.

Summary Variable (V), diversity (D), and joining (J) (V(D)J) recombination is the first determinant of antigen receptor diversity. Understanding how recombination is regulated requires a comprehensive, unbiased readout of V gene usage. We have developed VDJ sequencing (VDJ-seq), a DNA-based next-generation-sequencing technique that quantitatively profiles recombination products. We reveal a 200-fold range of recombination efficiency among recombining V genes in the primary mouse Igh repertoire. We used machine learning to integrate these data with local chromatin profiles to identify combinatorial patterns of epigenetic features that associate with active VH gene recombination. These features localize downstream of VH genes and are excised by recombination, revealing a class of cis-regulatory element that governs recombination, distinct from expression. We detect two mutually exclusive chromatin signatures at these elements, characterized by CTCF/RAD21 and PAX5/IRF4, which segregate with the evolutionary history of associated VH genes. Thus, local chromatin signatures downstream of VH genes provide an essential layer of regulation that determines recombination efficiency.

Circulating and Tissue-Resident CD4+ T Cells With Reactivity to Intestinal Microbiota Are Abundant in Healthy Individuals and Function Is Altered During Inflammation

Background & Aims Interactions between commensal microbes and the immune system are tightly regulated and maintain intestinal homeostasis, but little is known about these interactions in humans. We investigated responses of human CD4+ T cells to the intestinal microbiota. We measured the abundance of T cells in circulation and intestinal tissues that respond to intestinal microbes and determined their clonal diversity. We also assessed their functional phenotypes and effects on intestinal resident cell populations, and studied alterations in microbe-reactive T cells in patients with chronic intestinal inflammation. Methods We collected samples of peripheral blood mononuclear cells and intestinal tissues from healthy individuals (controls, n = 13−30) and patients with inflammatory bowel diseases (n = 119; 59 with ulcerative colitis and 60 with Crohn’s disease). We used 2 independent assays (CD154 detection and carboxy-fluorescein succinimidyl ester dilution assays) and 9 intestinal bacterial species (Escherichia coli, Lactobacillus acidophilus, Bifidobacterium animalis subsp lactis, Faecalibacterium prausnitzii, Bacteroides vulgatus, Roseburia intestinalis, Ruminococcus obeum, Salmonella typhimurium, and Clostridium difficile) to quantify, expand, and characterize microbe-reactive CD4+ T cells. We sequenced T-cell receptor Vβ genes in expanded microbe-reactive T-cell lines to determine their clonal diversity. We examined the effects of microbe-reactive CD4+ T cells on intestinal stromal and epithelial cell lines. Cytokines, chemokines, and gene expression patterns were measured by flow cytometry and quantitative polymerase chain reaction. Results Circulating and gut-resident CD4+ T cells from controls responded to bacteria at frequencies of 40−4000 per million for each bacterial species tested. Microbiota-reactive CD4+ T cells were mainly of a memory phenotype, present in peripheral blood mononuclear cells and intestinal tissue, and had a diverse T-cell receptor Vβ repertoire. These cells were functionally heterogeneous, produced barrier-protective cytokines, and stimulated intestinal stromal and epithelial cells via interleukin 17A, interferon gamma, and tumor necrosis factor. In patients with inflammatory bowel diseases, microbiota-reactive CD4+ T cells were reduced in the blood compared with intestine; T-cell responses that we detected had an increased frequency of interleukin 17A production compared with responses of T cells from blood or intestinal tissues of controls. Conclusions In an analysis of peripheral blood mononuclear cells and intestinal tissues from patients with inflammatory bowel diseases vs controls, we found that reactivity to intestinal bacteria is a normal property of the human CD4+ T-cell repertoire, and does not necessarily indicate disrupted interactions between immune cells and the commensal microbiota. T-cell responses to commensals might support intestinal homeostasis, by producing barrier-protective cytokines and providing a large pool of T cells that react to pathogens.