Interleukin-23 (IL-23) is a pro-inflammatory cytokine required for the pathogenicity of T helper 17 (Th17) cells but the molecular mechanisms governing this process remain unclear. We identified the transcription factor Blimp-1 (Prdm1) as a key IL-23-induced factor that drove the inflammatory function of Th17 cells. In contrast to thymic deletion of Blimp-1, which causes T cell development defects and spontaneous autoimmunity, peripheral deletion of this transcription factor resulted in reduced Th17 activation and reduced severity of autoimmune encephalomyelitis. Furthermore, genome-wide occupancy and overexpression studies in Th17 cells revealed that Blimp-1 co-localized with transcription factors RORγt, STAT-3, and p300 at the Il23r, Il17a/f, and Csf2 cytokine loci to enhance their expression. Blimp-1 also directly bound to and repressed cytokine loci Il2 and Bcl6. Taken together, our results demonstrate that Blimp-1 is an essential transcription factor downstream of IL-23 that acts in concert with RORγt to activate the Th17 inflammatory program.

Interleukin-27 (IL-27) is a key immunosuppressive cytokine that counters T helper 17 (Th17) cell-mediated pathology. To identify mechanisms by which IL-27 might exert its immunosuppressive effect, we analyzed genes in T cells rapidly induced by IL-27. We found that IL-27 priming of naive T cells upregulated expression of programmed death ligand 1 (PD-L1) in a signal transducer and activator of transcription 1 (STAT1)-dependent manner. When cocultured with naive CD4(+) T cells, IL-27-primed T cells inhibited the differentiation of Th17 cells in trans through a PD-1-PD-L1 interaction. In vivo, coadministration of naive TCR transgenic T cells (2D2 T cells) with IL-27-primed T cells expressing PD-L1 inhibited the development of Th17 cells and protected from severe autoimmune encephalomyelitis. Thus, these data identify a suppressive activity of IL-27, by which CD4(+) T cells can restrict differentiation of Th17 cells in trans.

T cell development is orchestrated by transcription factors that regulate the expression of genes initially buried within inaccessible chromatin, but the transcription factors that establish the regulatory landscape of the T cell lineage remain unknown. Profiling chromatin accessibility at eight stages of T cell development revealed the selective enrichment of TCF-1 at genomic regions that became accessible at the earliest stages of development. TCF-1 was further required for the accessibility of these regulatory elements and at the single-cell level, it dictated a coordinate opening of chromatin in T cells. TCF-1 expression in fibroblasts generated de novo chromatin accessibility even at chromatin regions with repressive marks, inducing the expression of T cell-restricted genes. These results indicate that a mechanism by which TCF-1 controls T cell fate is through its widespread ability to target silent chromatin and establish the epigenetic identity of T cells.

Genetics is a major determinant of susceptibility to autoimmune disorders. Here, we examined whether genome organization provides resilience or susceptibility to sequence variations, and how this would contribute to the molecular etiology of an autoimmune disease. We generated high-resolution maps of linear and 3D genome organization in thymocytes of NOD mice, a model of type 1 diabetes (T1D), and the diabetes-resistant C57BL/6 mice. Multi-enhancer interactions formed at genomic regions harboring genes with prominent roles in T cell development in both strains. However, diabetes risk-conferring loci coalesced enhancers and promoters in NOD, but not C57BL/6 thymocytes. 3D genome mapping of NODxC57BL/6 F1 thymocytes revealed that genomic misfolding in NOD mice is mediated in cis. Moreover, immune cells infiltrating the pancreas of humans with T1D exhibited increased expression of genes located on misfolded loci in mice. Thus, genetic variation leads to altered 3D chromatin architecture and associated changes in gene expression that may underlie autoimmune pathology.

Exhausted CD8 T cells (TEX) are a distinct state of T cell differentiation associated with failure to clear chronic viruses and cancer. Immunotherapies such as PD-1 blockade can reinvigorate TEX cells, but reinvigoration is not durable. A major unanswered question is whether TEX cells differentiate into functional durable memory T cells (TMEM) upon antigen clearance. Here, using a mouse model, we found that upon eliminating chronic antigenic stimulation, TEX cells partially (re)acquire phenotypic and transcriptional features of TMEM cells. These 'recovering' TEX cells originated from the T cell factor (TCF-1+) TEX progenitor subset. Nevertheless, the recall capacity of these recovering TEX cells remained compromised as compared to TMEM cells. Chromatin-accessibility profiling revealed a failure to recover core memory epigenetic circuits and maintenance of a largely exhausted open chromatin landscape. Thus, despite some phenotypic and transcriptional recovery upon antigen clearance, exhaustion leaves durable epigenetic scars constraining future immune responses. These results support epigenetic remodeling interventions for TEX cell-targeted immunotherapies.

Chimeric antigen receptor (CAR)-T immunotherapy has yielded impressive results in several B cell malignancies, establishing itself as a powerful means to redirect the natural properties of T lymphocytes. In this strategy, the T cell genome is modified by the integration of lentiviral vectors encoding CAR that direct tumor cell killing. However, this therapeutic approach is often limited by the extent of CAR-T cell expansion in vivo. A major outstanding question is whether or not CAR-T integration itself enhances the proliferative competence of individual T cells by rewiring their regulatory landscape. To address this question, it is critical to define the identity of an individual CAR-T cell and simultaneously chart where the CAR-T vector integrates into the genome. Here, we report the development of a method called EpiVIA (https://github.com/VahediLab/epiVIA) for the joint profiling of the chromatin accessibility and lentiviral integration site analysis at the population and single-cell levels. We validate our technique in clonal cells with previously defined integration sites and further demonstrate the ability to measure lentiviral integration sites and chromatin accessibility of host and viral genomes at the single-cell resolution in CAR-T cells. We anticipate that EpiVIA will enable the single-cell deconstruction of gene regulation during CAR-T therapy, leading to the discovery of cellular factors associated with durable treatment.

Lymphocyte migration is essential for adaptive immune surveillance. However, our current understanding of this process is rudimentary, because most human studies have been restricted to immunological analyses of blood and various tissues. To address this knowledge gap, we used an integrated approach to characterize tissue-emigrant lineages in thoracic duct lymph (TDL). The most prevalent immune cells in human and non-human primate efferent lymph were T cells. Cytolytic CD8+ T cell subsets with effector-like epigenetic and transcriptional signatures were clonotypically skewed and selectively confined to the intravascular circulation, whereas non-cytolytic CD8+ T cell subsets with stem-like epigenetic and transcriptional signatures predominated in tissues and TDL. Moreover, these anatomically distinct gene expression profiles were recapitulated within individual clonotypes, suggesting parallel differentiation programs independent of the expressed antigen receptor. Our collective dataset provides an atlas of the migratory immune system and defines the nature of tissue-emigrant CD8+ T cells that recirculate via TDL.

The generation of high-affinity neutralizing antibodies, the objective of most vaccine strategies, occurs in B cells within germinal centers (GCs) and requires rate-limiting "help" from follicular helper CD4+ T (Tfh) cells. Although Tfh differentiation is an attribute of MHC II-restricted CD4+ T cells, the transcription factors driving Tfh differentiation, notably Bcl6, are not restricted to CD4+ T cells. Here, we identified a requirement for the CD4+-specific transcription factor Thpok during Tfh cell differentiation, GC formation, and antibody maturation. Thpok promoted Bcl6 expression and bound to a Thpok-responsive region in the first intron of Bcl6. Thpok also promoted the expression of Bcl6-independent genes, including the transcription factor Maf, which cooperated with Bcl6 to mediate the effect of Thpok on Tfh cell differentiation. Our findings identify a transcriptional program that links the CD4+ lineage with Tfh differentiation, a limiting factor for efficient B cell responses, and suggest avenues to optimize vaccine generation.

Through their functional diversification, distinct lineages of CD4(+) T cells can act to either drive or constrain immune-mediated pathology. Transcription factors are critical in the generation of cellular diversity, and negative regulators antagonistic to alternate fates often act in conjunction with positive regulators to stabilize lineage commitment. Genetic polymorphisms within a single locus encoding the transcription factor BACH2 are associated with numerous autoimmune and allergic diseases including asthma, Crohn's disease, coeliac disease, vitiligo, multiple sclerosis and type 1 diabetes. Although these associations point to a shared mechanism underlying susceptibility to diverse immune-mediated diseases, a function for BACH2 in the maintenance of immune homeostasis has not been established. Here, by studying mice in which the Bach2 gene is disrupted, we define BACH2 as a broad regulator of immune activation that stabilizes immunoregulatory capacity while repressing the differentiation programs of multiple effector lineages in CD4(+) T cells. BACH2 was required for efficient formation of regulatory (Treg) cells and consequently for suppression of lethal inflammation in a manner that was Treg-cell-dependent. Assessment of the genome-wide function of BACH2, however, revealed that it represses genes associated with effector cell differentiation. Consequently, its absence during Treg polarization resulted in inappropriate diversion to effector lineages. In addition, BACH2 constrained full effector differentiation within TH1, TH2 and TH17 cell lineages. These findings identify BACH2 as a key regulator of CD4(+) T-cell differentiation that prevents inflammatory disease by controlling the balance between tolerance and immunity.

Enhancers regulate spatiotemporal gene expression and impart cell-specific transcriptional outputs that drive cell identity. Super-enhancers (SEs), also known as stretch-enhancers, are a subset of enhancers especially important for genes associated with cell identity and genetic risk of disease. CD4(+) T cells are critical for host defence and autoimmunity. Here we analysed maps of mouse T-cell SEs as a non-biased means of identifying key regulatory nodes involved in cell specification. We found that cytokines and cytokine receptors were the dominant class of genes exhibiting SE architecture in T cells. Nonetheless, the locus encoding Bach2, a key negative regulator of effector differentiation, emerged as the most prominent T-cell SE, revealing a network in which SE-associated genes critical for T-cell biology are repressed by BACH2. Disease-associated single-nucleotide polymorphisms for immune-mediated disorders, including rheumatoid arthritis, were highly enriched for T-cell SEs versus typical enhancers or SEs in other cell lineages. Intriguingly, treatment of T cells with the Janus kinase (JAK) inhibitor tofacitinib disproportionately altered the expression of rheumatoid arthritis risk genes with SE structures. Together, these results indicate that genes with SE architecture in T cells encompass a variety of cytokines and cytokine receptors but are controlled by a 'guardian' transcription factor, itself endowed with an SE. Thus, enumeration of SEs allows the unbiased determination of key regulatory nodes in T cells, which are preferentially modulated by pharmacological intervention.

Mutations in MECP2 cause Rett syndrome (RTT), an X-linked neurological disorder characterized by regressive loss of neurodevelopmental milestones and acquired psychomotor deficits. However, the cellular heterogeneity of the brain impedes an understanding of how MECP2 mutations contribute to RTT. Here we developed a Cre-inducible method for cell-type-specific biotin tagging of MeCP2 in mice. Combining this approach with an allelic series of knock-in mice carrying frequent RTT-associated mutations (encoding T158M and R106W) enabled the selective profiling of RTT-associated nuclear transcriptomes in excitatory and inhibitory cortical neurons. We found that most gene-expression changes were largely specific to each RTT-associated mutation and cell type. Lowly expressed cell-type-enriched genes were preferentially disrupted by MeCP2 mutations, with upregulated and downregulated genes reflecting distinct functional categories. Subcellular RNA analysis in MeCP2-mutant neurons further revealed reductions in the nascent transcription of long genes and uncovered widespread post-transcriptional compensation at the cellular level. Finally, we overcame X-linked cellular mosaicism in female RTT models and identified distinct gene-expression changes between neighboring wild-type and mutant neurons, providing contextual insights into RTT etiology that support personalized therapeutic interventions.

Signaling pathways are intimately involved in cellular differentiation, allowing cells to respond to their environment by regulating gene expression. Although enhancers are recognized as key elements that regulate selective gene expression, the interplay between signaling pathways and actively used enhancer elements is not clear. Here, we use CD4(+) T cells as a model of differentiation, mapping the activity of cell-type-specific enhancer elements in T helper 1 (Th1) and Th2 cells. Our data establish that STAT proteins have a major impact on the activation of lineage-specific enhancers and the suppression of enhancers associated with alternative cell fates. Transcriptome analysis further supports a functional role for enhancers regulated by STATs. Importantly, expression of lineage-defining master regulators in STAT-deficient cells fails to fully recover the chromatin signature of STAT-dependent enhancers. Thus, these findings point to a critical role of STATs as environmental sensors in dynamically molding the specialized enhancer architecture of differentiating cells.

Interleukin-6 (IL-6) and IL-27 signal through a shared receptor subunit and employ the same downstream STAT transcription proteins, but yet are ascribed unique and overlapping functions. To evaluate the specificity and redundancy for these cytokines, we quantified their global transcriptomic changes and determined the relative contributions of STAT1 and STAT3 using genetic models and chromatin immunoprecipitation-sequencing (ChIP-seq) approaches. We found an extensive overlap of the transcriptomes induced by IL-6 and IL-27 and few examples in which the cytokines acted in opposition. Using STAT-deficient cells and T cells from patients with gain-of-function STAT1 mutations, we demonstrated that STAT3 is responsible for the overall transcriptional output driven by both cytokines, whereas STAT1 is the principal driver of specificity. STAT1 cannot compensate in the absence of STAT3 and, in fact, much of STAT1 binding to chromatin is STAT3 dependent. Thus, STAT1 shapes the specific cytokine signature superimposed upon STAT3's action.

Tumour necrosis factor (TNF) superfamily cytokines and their receptors regulate diverse immune system functions through a common set of signalling pathways. Genetic variants in and expression of individual TNF superfamily cytokines, receptors and signalling proteins have been associated with autoimmune and inflammatory diseases, but their interconnected biology has been largely unexplored.

Cellular composition and anatomical organization influence normal and aberrant organ functions. Emerging spatial single-cell proteomic assays such as Image Mass Cytometry (IMC) and Co-Detection by Indexing (CODEX) have facilitated the study of cellular composition and organization by enabling high-throughput measurement of cells and their localization directly in intact tissues. However, annotation of cell types and quantification of their relative localization in tissues remain challenging. To address these unmet needs, we developed AnnoSpat (Annotator and Spatial Pattern Finder) that uses neural network and point process algorithms to automatically identify cell types and quantify cell-cell proximity relationships. Our study of data from IMC and CODEX show the superior performance of AnnoSpat in rapid and accurate annotation of cell types compared to alternative approaches. Moreover, the application of AnnoSpat to type 1 diabetic, non-diabetic autoantibody-positive, and non-diabetic organ donor cohorts recapitulated known islet pathobiology and showed differential dynamics of pancreatic polypeptide (PP) cell abundance and CD8 + T cells infiltration in islets during type 1 diabetes progression.

The roles of EZH2 in various subsets of CD4(+) T cells are controversial and its mechanisms of action are incompletely understood. FOXP3-positive Treg cells are a critical helper T cell subset, and dysregulation of Treg generation or function results in systemic autoimmunity. FOXP3 associates with EZH2 to mediate gene repression and suppressive function. Herein, we demonstrate that deletion of Ezh2 in CD4 T cells resulted in reduced numbers of Treg cells in vivo and differentiation in vitro and an increased proportion of memory CD4 T cells in part due to exaggerated production of effector cytokines. Furthermore, we found that both Ezh2-deficient Treg cells and T effector cells were functionally impaired in vivo: Tregs failed to constrain autoimmune colitis and T effector cells neither provided a protective response to T. gondii infection nor mediated autoimmune colitis. The dichotomous function of EZH2 in regulating differentiation and senescence in effector and regulatory T cells helps to explain the apparent existing contradictions in literature.

Interleukin (IL)-22-producing innate lymphoid cells (ILCs; ILC22) comprise a heterogeneous population of cells that are dependent on the transcription factor retinoid-related orphan γt (RORγt) and are critical for barrier function of the intestinal mucosa. A distinct ILC22 subset expresses the natural cytotoxicity receptor NKp46 (NKp46+ ILC22); however, the factors that contribute to the generation of this population versus other subsets are largely unknown. Herein, we show that T-bet (encoded by Tbx21) was highly expressed in NKp46+ ILC22, a feature shared by all NKp46+ cells present in the intestine but not by other IL-22-producing populations. Accordingly, the absence of T-bet resulted in loss of NKp46+ ILC22 in the intestinal lamina propria. The residual NKp46+ ILC22 present in Tbx21(-/-) mice showed a marked reduction of Rorγt expression and impairment in IL-22 production. Generation and functions of gut NK1.1+ cells were also altered. Bone marrow chimera experiments revealed a cell-intrinsic requirement for T-bet in these subsets and competitive reconstitution experiments revealed roles for T-bet in multiple ILC subsets. Thus, T-bet has a general importance for ILC in the gut and plays a selective and critical role in the generation of NKp46+ ILC22.

Histone post-translational modifications (PTMs) contribute to chromatin accessibility due to their chemical properties and their ability to recruit enzymes responsible for DNA readout and chromatin remodeling. To date, more than 400 different histone PTMs and thousands of combinations of PTMs have been identified, the vast majority with still unknown biological function. Identification and quantification of histone PTMs has become routine in mass spectrometry (MS) but, since raising antibodies for each PTM in a study can be prohibitive, lots of potential is lost from MS datasets when uncharacterized PTMs are found to be significantly regulated. We developed an assay that uses metabolic labeling and MS to associate chromatin accessibility with histone PTMs and their combinations. The labeling is achieved by spiking in the cell media a 5x concentration of stable isotope labeled arginine and allow cells to grow for at least one cell cycle. We quantified the labeling incorporation of about 200 histone peptides with a proteomics workflow, and we confirmed that peptides carrying PTMs with extensively characterized roles in active transcription or gene silencing were in highly or poorly labeled forms, respectively. Data were further validated using next-generation sequencing to assess the transcription rate of chromatin regions modified with five selected PTMs. Furthermore, we quantified the labeling rate of peptides carrying co-existing PTMs, proving that this method is suitable for combinatorial PTMs. We focus on the abundant bivalent mark H3K27me3K36me2, showing that H3K27me3 dominantly represses histone swapping rate even in the presence of the more permissive PTM H3K36me2. Together, we envision this method will help to generate hypotheses regarding histone PTM functions and, potentially, elucidate the role of combinatorial histone codes.

TCF-1 is a key transcription factor in progenitor exhausted CD8 T cells (Tex). Moreover, this Tex cell subset mediates responses to PD-1 checkpoint pathway blockade. However, the role of the transcription factor TCF-1 in early fate decisions and initial generation of Tex cells is unclear. Single-cell RNA sequencing (scRNA-seq) and lineage tracing identified a TCF-1+Ly108+PD-1+ CD8 T cell population that seeds development of mature Tex cells early during chronic infection. TCF-1 mediated the bifurcation between divergent fates, repressing development of terminal KLRG1Hi effectors while fostering KLRG1Lo Tex precursor cells, and PD-1 stabilized this TCF-1+ Tex precursor cell pool. TCF-1 mediated a T-bet-to-Eomes transcription factor transition in Tex precursors by promoting Eomes expression and drove c-Myb expression that controlled Bcl-2 and survival. These data define a role for TCF-1 in early-fate-bifurcation-driving Tex precursor cells and also identify PD-1 as a protector of this early TCF-1 subset.

Type 1 diabetes (T1D) is an autoimmune disease in which immune cells destroy insulin-producing beta cells. The aetiology of this complex disease is dependent on the interplay of multiple heterogeneous cell types in the pancreatic environment. Here, we provide a single-cell atlas of pancreatic islets of 24 T1D, autoantibody-positive and nondiabetic organ donors across multiple quantitative modalities including ~80,000 cells using single-cell transcriptomics, ~7,000,000 cells using cytometry by time of flight and ~1,000,000 cells using in situ imaging mass cytometry. We develop an advanced integrative analytical strategy to assess pancreatic islets and identify canonical cell types. We show that a subset of exocrine ductal cells acquires a signature of tolerogenic dendritic cells in an apparent attempt at immune suppression in T1D donors. Our multimodal analyses delineate cell types and processes that may contribute to T1D immunopathogenesis and provide an integrative procedure for exploration and discovery of human pancreatic function.