Dynamic recruitment of Ets1 to both nucleosome-occupied and -depleted enhancer regions mediates a transcriptional program switch during early T-cell differentiation

Regulation of T helper cell differentiation by the interplay between histone modification and chromatin interaction

The development and function of the immune system are controlled by temporospatial gene expression programs, which are regulated by cis-regulatory elements, chromatin structure, and trans-acting factors. In this study, we cataloged the dynamic histone modifications and chromatin interactions at regulatory regions during T helper (Th) cell differentiation. Our data revealed that the H3K4me1 landscape established by MLL4 in naive CD4+ T cells is critical for restructuring the regulatory interaction network and orchestrating gene expression during the early phase of Th differentiation. GATA3 plays a crucial role in further configuring H3K4me1 modification and the chromatin interaction network during Th2 differentiation. Furthermore, we demonstrated that HSS3-anchored chromatin loops function to restrict the activity of the Th2 locus control region (LCR), thus coordinating the expression of Th2 cytokines. Our results provide insights into the mechanisms of how the interplay between histone modifications, chromatin looping, and trans-acting factors contributes to the differentiation of Th cells.

Multimodal stimulation screens reveal unique and shared genes limiting T cell fitness

Genes limiting T cell antitumor activity may serve as therapeutic targets. It has not been systematically studied whether there are regulators that uniquely or broadly contribute to T cell fitness. We perform genome-scale CRISPR-Cas9 knockout screens in primary CD8 T cells to uncover genes negatively impacting fitness upon three modes of stimulation: (1) intense, triggering activation-induced cell death (AICD); (2) acute, triggering expansion; (3) chronic, causing dysfunction. Besides established regulators, we uncover genes controlling T cell fitness either specifically or commonly upon differential stimulation. Dap5 ablation, ranking highly in all three screens, increases translation while enhancing tumor killing. Loss of Icam1-mediated homotypic T cell clustering amplifies cell expansion and effector functions after both acute and intense stimulation. Lastly, Ctbp1 inactivation induces functional T cell persistence exclusively upon chronic stimulation. Our results functionally annotate fitness regulators based on their unique or shared contribution to traits limiting T cell antitumor activity.

Investigating open access new approach methods (NAM) to assess biological points of departure: A case study with 4 neurotoxic pesticides

Open access new approach methods (NAM) in the US EPA ToxCast program and NTP Integrated Chemical Environment (ICE) were used to investigate activities of four neurotoxic pesticides: endosulfan, fipronil, propyzamide and carbaryl. Concordance of in vivo regulatory points of departure (POD) adjusted for interspecies extrapolation (AdjPOD) to modelled human Administered Equivalent Dose (AEDHuman) was assessed using 3-compartment or Adult/Fetal PBTK in vitro to in vivo extrapolation. Model inputs were from Tier 1 (High throughput transcriptomics: HTTr, high throughput phenotypic profiling: HTPP) and Tier 2 (single target: ToxCast) assays. HTTr identified gene expression signatures associated with potential neurotoxicity for endosulfan, propyzamide and carbaryl in non-neuronal MCF-7 and HepaRG cells. The HTPP assay in U-2 OS cells detected potent effects on DNA endpoints for endosulfan and carbaryl, and mitochondria with fipronil (propyzamide was inactive). The most potent ToxCast assays were concordant with specific components of each chemical mode of action (MOA). Predictive adult IVIVE models produced fold differences (FD) < 10 between the AEDHuman and the measured in vivo AdjPOD. The 3-compartment model was concordant (i.e., smallest FD) for endosulfan, fipronil and carbaryl, and PBTK was concordant for propyzamide. The most potent AEDHuman predictions for each chemical showed HTTr, HTPP and ToxCast were mainly concordant with in vivo AdjPODs but assays were less concordant with MOAs. This was likely due to the cell types used for testing and/or lack of metabolic capabilities and pathways available in vivo. The Fetal PBTK model had larger FDs than adult models and was less predictive overall.

Transcriptional reprogramming by mutated IRF4 in lymphoma

Disease-causing mutations in genes encoding transcription factors (TFs) can affect TF interactions with their cognate DNA-binding motifs. Whether and how TF mutations impact upon the binding to TF composite elements (CE) and the interaction with other TFs is unclear. Here, we report a distinct mechanism of TF alteration in human lymphomas with perturbed B cell identity, in particular classic Hodgkin lymphoma. It is caused by a recurrent somatic missense mutation c.295 T > C (p.Cys99Arg; p.C99R) targeting the center of the DNA-binding domain of Interferon Regulatory Factor 4 (IRF4), a key TF in immune cells. IRF4-C99R fundamentally alters IRF4 DNA-binding, with loss-of-binding to canonical IRF motifs and neomorphic gain-of-binding to canonical and non-canonical IRF CEs. IRF4-C99R thoroughly modifies IRF4 function by blocking IRF4-dependent plasma cell induction, and up-regulates disease-specific genes in a non-canonical Activator Protein-1 (AP-1)-IRF-CE (AICE)-dependent manner. Our data explain how a single mutation causes a complex switch of TF specificity and gene regulation and open the perspective to specifically block the neomorphic DNA-binding activities of a mutant TF.

Recruitment of transcription factor ETS1 to activated accessible regions promotes the transcriptional program of cilia genes

Defects in cilia genes, which are critical for cilia formation and function, can cause complicated ciliopathy syndromes involving multiple organs and tissues; however, the underlying regulatory mechanisms of the networks of cilia genes in ciliopathies remain enigmatic. Herein, we have uncovered the genome-wide redistribution of accessible chromatin regions and extensive alterations of expression of cilia genes during Ellis-van Creveld syndrome (EVC) ciliopathy pathogenesis. Mechanistically, the distinct EVC ciliopathy-activated accessible regions (CAAs) are shown to positively regulate robust changes in flanking cilia genes, which are a key requirement for cilia transcription in response to developmental signals. Moreover, a single transcription factor, ETS1, can be recruited to CAAs, leading to prominent chromatin accessibility reconstruction in EVC ciliopathy patients. In zebrafish, the collapse of CAAs driven by ets1 suppression subsequently causes defective cilia proteins, resulting in body curvature and pericardial oedema. Our results depict a dynamic landscape of chromatin accessibility in EVC ciliopathy patients, and uncover an insightful role for ETS1 in controlling the global transcriptional program of cilia genes by reprogramming the widespread chromatin state.

Batf stabilizes Th17 cell development via impaired Stat5 recruitment of Ets1-Runx1 complexes

Although the activator protein-1 (AP-1) factor Batf is required for Th17 cell development, its mechanisms of action to underpin the Th17 program are incompletely understood. Here, we find that Batf ensures Th17 cell identity in part by restricting alternative gene programs through its actions to restrain IL-2 expression and IL-2-induced Stat5 activation. This, in turn, limits Stat5-dependent recruitment of Ets1-Runx1 factors to Th1- and Treg-cell-specific gene loci. Thus, in addition to pioneering regulatory elements in Th17-specific loci, Batf acts indirectly to inhibit the assembly of a Stat5-Ets1-Runx1 complex that enhances the transcription of Th1- and Treg-cell-specific genes. These findings unveil an important role for Stat5-Ets1-Runx1 interactions in transcriptional networks that define alternate T cell fates and indicate that Batf plays an indispensable role in both inducing and maintaining the Th17 program through its actions to regulate the competing actions of Stat5-assembled enhanceosomes that promote Th1- and Treg-cell developmental programs.

Heterogeneity of T cells and macrophages in chlorine-induced acute lung injury in mice using single-cell RNA sequencing

OBJECTIVE

Chlorine (Cl₂), as an asphyxiant toxicant, induced poisoning incidents and acute lung injury (ALI) occur frequently. The specific pathogenesis of Cl₂-induced ALI remains unclear. Immune cells play an important role in the process of lung damage. We used single-cell RNA sequencing (scRNA-seq) technology to explore T cells and macrophages molecular mechanism.

METHODS

Female BALB/c mice were exposed to 400 ppm Cl₂ for 15 min. scRNA-seq technology was used to observe the heterogeneity of T cells and macrophages. Hematoxylin-eosin (H&E) staining was used to evaluate the degree of lung injury. Immunofluorescence was used to verify the highly expressed genes of our interest.

RESULTS

A total of 5316 to 7742 cells were classified into eight different cell types. Several new highly expressed anti-inflammatory and pro-inflammatory genes were found in T cells and macrophages, which were further verified in vitro. Through the pseudotime analysis of macrophages, it was found that the expression of pro-inflammatory and anti-inflammatory genes showed opposite trends in the development of Cl₂-induced ALI. This study also mapped T cells-macrophage communication and identified the development of several important receptor-ligand complexes in Cl₂-induced ALI.

CONCLUSIONS

These findings are worthy of further exploration and provide new resources and directions for the study of Cl₂-induced ALI in mice, especially in immune and inflammation mechanisms.

Transcriptional Regulation of Early T-Lymphocyte Development in Thymus

T-lymphocytes play crucial roles for maintaining immune homeostasis by fighting against various pathogenic microorganisms and establishing self-antigen tolerance. They will go through several stages and checkpoints in the thymus from progenitors to mature T cells, from CD4^-CD8^- double negative (DN) cells to CD4⁺CD8⁺ double positive (DP) cells, finally become CD4⁺ or CD8⁺ single positive (SP) cells. The mature SP cells then emigrate out of the thymus and further differentiate into distinct subsets under different environment signals to perform specific functions. Each step is regulated by various transcriptional regulators downstream of T cell receptors (TCRs) that have been extensively studied both in vivo and vitro via multiple mouse models and advanced techniques, such as single cell RNA sequencing (scRNA-seq) and Chromatin Immunoprecipitation sequencing (ChIP-seq). This review will summarize the transcriptional regulators participating in the early stage of T cell development reported in the past decade, trying to figure out cascade networks in each process and provide possible research directions in the future.

Inflammatory Immune-Associated eRNA: Mechanisms, Functions and Therapeutic Prospects

The rapid development of multiple high-throughput sequencing technologies has made it possible to explore the critical roles and mechanisms of functional enhancers and enhancer RNAs (eRNAs). The inflammatory immune response, as a fundamental pathological process in infectious diseases, cancers and immune disorders, coordinates the balance between the internal and external environment of the organism. It has been shown that both active enhancers and intranuclear eRNAs are preferentially expressed over inflammation-related genes in response to inflammatory stimuli, suggesting that enhancer transcription events and their products influence the expression and function of inflammatory genes. Therefore, in this review, we summarize and discuss the relevant inflammatory roles and regulatory mechanisms of eRNAs in inflammatory immune cells, non-inflammatory immune cells, inflammatory immune diseases and tumors, and explore the potential therapeutic effects of enhancer inhibitors affecting eRNA production for diseases with inflammatory immune responses.

NKG7 Is a T-cell-Intrinsic Therapeutic Target for Improving Antitumor Cytotoxicity and Cancer Immunotherapy

Cytotoxic CD8⁺ T cells (CTL) are a crucial component of the immune system notable for their ability to eliminate rapidly proliferating malignant cells. However, the T-cell intrinsic factors required for human CTLs to accomplish highly efficient antitumor cytotoxicity are not well defined. By evaluating human CD8⁺ T cells from responders versus nonresponders to treatment with immune checkpoint inhibitors, we sought to identify key factors associated with effective CTL function. Single-cell RNA-sequencing analysis of peripheral CD8⁺ T cells from patients treated with anti-PD-1 therapy showed that cells from nonresponders exhibited decreased expression of the cytolytic granule-associated molecule natural killer cell granule protein-7 (NKG7). Functional assays revealed that reduced NKG7 expression altered cytolytic granule number, trafficking, and calcium release, resulting in decreased CD8⁺ T-cell-mediated killing of tumor cells. Transfection of T cells with NKG7 mRNA was sufficient to improve the tumor-cell killing ability of human T cells isolated from nonresponders and increase their response to anti-PD-1 or anti-PD-L1 therapy in vitro. NKG7 mRNA therapy also improved the antitumor activity of murine tumor antigen-specific CD8⁺ T cells in an in vivo model of adoptive cell therapy. Finally, we showed that the transcription factor ETS1 played a role in regulating NKG7 expression. Together, our results identify NKG7 as a necessary component for the cytotoxic function of CD8⁺ T cells and establish NKG7 as a T-cell-intrinsic therapeutic target for enhancing cancer immunotherapy.See related article by Li et al., p. 154.

Multilayer omics analysis reveals a non-classical retinoic acid signaling axis that regulates hematopoietic stem cell identity

Hematopoietic stem cells (HSCs) rely on complex regulatory networks to preserve stemness. Due to the scarcity of HSCs, technical challenges have limited our insights into the interplay between metabolites, transcription, and the epigenome. In this study, we generated low-input metabolomics, transcriptomics, chromatin accessibility, and chromatin immunoprecipitation data, revealing distinct metabolic hubs that are enriched in HSCs and their downstream multipotent progenitors. Mechanistically, we uncover a non-classical retinoic acid (RA) signaling axis that regulates HSC function. We show that HSCs rely on Cyp26b1, an enzyme conventionally considered to limit RA effects in the cell. In contrast to the traditional view, we demonstrate that Cyp26b1 is indispensable for production of the active metabolite 4-oxo-RA. Further, RA receptor beta (Rarb) is required for complete transmission of 4-oxo-RA-mediated signaling to maintain stem cells. Our findings emphasize that a single metabolite controls stem cell fate by instructing epigenetic and transcriptional attributes.

Sequences to Differences in Gene Expression: Analysis of RNA-Seq Data

RNA-Seq is now a routinely employed assay to measure gene expression. As the technique matured over the last decade, so have dedicated analytic tools. In this chapter, we first describe the mainstream as well as the most up-to-date protocols and their implications on downstream analysis. We then detail the steps entailing RNA-Seq analysis in three main stages: (i) preprocessing and data preparation, (ii) upstream processing, and (iii) high-level analyses. We review the most recent and relevant tools as one workflow following a stepwise order. The chapter further encompasses in-depth features of these tools. Details of the required code are made available throughout the chapter, as well as of the underlying statistics. We illustrate these steps with analysis of publicly available RNA-Seq data.

Arid1a promotes thymocyte development through β-selection-dependent and β-selection-independent mechanisms

Early T-cell development from CD4^-  CD8^- double-negative (DN) stage to CD4⁺  CD8⁺ double-positive (DP) stage in the thymus is regulated through multiple steps involving a batch of sequentially expressed factors. Our preliminary data and a recent report showed that AT-rich interaction domain 1A (Arid1a) is required for the transition from DN to DP stages, but the mechanism is not fully understood. In this study, we consolidated that conditional deletion of Arid1a in T-cell lineage intrinsically caused developmental blocks from DN3 to DN4 stages, as well as from DN4 to DP stages using both in vivo adoptive T-cell transfer model and in vitro culture system. The expression of intracellular TCRβ is significantly decreased in Arid1a-deficient DN4 cells compared with WT cells. OT1 transgenic TCR can rescue the defect in the transition from DN3 to DN4 stages, but not from DN to DP stages. Furthermore, we observed a comparable or stronger proliferation capacity accompanied by a significant increase in cell death in Arid1a^-/- DP cells compared with that in WT controls. RNA-Seq analysis shows a significant enrichment of apoptotic pathway within differentially expressed genes between Arid1a^-/- and WT DP cells, including the upregulation of Bim, Casp3 and Trp53 and the downregulation of Rorc, Bcl-XL and Mcl1. Therefore, our study reveals a novel mechanism that Arid1a controls early T-cell development by maintaining intracellular TCRβ expression-mediated β-selection and activating parallel cell survival pathways.

Global patterns of enhancer activity during sea urchin embryogenesis assessed by eRNA profiling

We used capped analysis of gene expression with sequencing (CAGE-seq) to profile eRNA expression and enhancer activity during embryogenesis of a model echinoderm: the sea urchin, Strongylocentrotus purpuratus. We identified more than 18,000 enhancers that were active in mature oocytes and developing embryos and documented a burst of enhancer activation during cleavage and early blastula stages. We found that a large fraction (73.8%) of all enhancers active during the first 48 h of embryogenesis were hyperaccessible no later than the 128-cell stage and possibly even earlier. Most enhancers were located near gene bodies, and temporal patterns of eRNA expression tended to parallel those of nearby genes. Furthermore, enhancers near lineage-specific genes contained signatures of inputs from developmental gene regulatory networks deployed in those lineages. A large fraction (60%) of sea urchin enhancers previously shown to be active in transgenic reporter assays was associated with eRNA expression. Moreover, a large fraction (50%) of a representative subset of enhancers identified by eRNA profiling drove tissue-specific gene expression in isolation when tested by reporter assays. Our findings provide an atlas of developmental enhancers in a model sea urchin and support the utility of eRNA profiling as a tool for enhancer discovery and regulatory biology. The data generated in this study are available at Echinobase, the public database of information related to echinoderm genomics.

ZFP451-mediated SUMOylation of SATB2 drives embryonic stem cell differentiation

Here, Urrutia et al. set out to study the mechanism that regulates the choice between pluripotency and differentiation in embryonic stem cells (ESCs). Using biochemical and genomic analyses, the authors identify SUMO2 modification of Satb2 by the E3 ligase Zfp451 as a driver of ESC differentiation.

The establishment of cell fates involves alterations of transcription factor repertoires and repurposing of transcription factors by post-translational modifications. In embryonic stem cells (ESCs), the chromatin organizers SATB2 and SATB1 balance pluripotency and differentiation by activating and repressing pluripotency genes, respectively. Here, we show that conditional Satb2 gene inactivation weakens ESC pluripotency, and we identify SUMO2 modification of SATB2 by the E3 ligase ZFP451 as a potential driver of ESC differentiation. Mutations of two SUMO-acceptor lysines of Satb2 (Satb2^{K →R}) or knockout of Zfp451 impair the ability of ESCs to silence pluripotency genes and activate differentiation-associated genes in response to retinoic acid (RA) treatment. Notably, the forced expression of a SUMO2-SATB2 fusion protein in either Satb2^{K →R} or Zfp451^−/− ESCs rescues, in part, their impaired differentiation potential and enhances the down-regulation of Nanog. The differentiation defect of Satb2^{K →R} ESCs correlates with altered higher-order chromatin interactions relative to Satb2^wt ESCs. Upon RA treatment of Satb2^wt ESCs, SATB2 interacts with ZFP451 and the LSD1/CoREST complex and gains binding at differentiation genes, which is not observed in RA-treated Satb2^{K →R} cells. Thus, SATB2 SUMOylation may contribute to the rewiring of transcriptional networks and the chromatin interactome of ESCs in the transition of pluripotency to differentiation.

Blueprint of human thymopoiesis reveals molecular mechanisms of stage-specific TCR enhancer activation

Cell differentiation is accompanied by epigenetic changes leading to precise lineage definition and cell identity. Here we present a comprehensive resource of epigenomic data of human T cell precursors along with an integrative analysis of other hematopoietic populations. Although T cell commitment is accompanied by large scale epigenetic changes, we observed that the majority of distal regulatory elements are constitutively unmethylated throughout T cell differentiation, irrespective of their activation status. Among these, the TCRA gene enhancer (Eα) is in an open and unmethylated chromatin structure well before activation. Integrative analyses revealed that the HOXA5-9 transcription factors repress the Eα enhancer at early stages of T cell differentiation, while their decommission is required for TCRA locus activation and enforced αβ T lineage differentiation. Remarkably, the HOXA-mediated repression of Eα is paralleled by the ectopic expression of homeodomain-related oncogenes in T cell acute lymphoblastic leukemia. These results highlight an analogous enhancer repression mechanism at play in normal and cancer conditions, but imposing distinct developmental constraints.

Regulation of T-cell Receptor Gene Expression by Three-Dimensional Locus Conformation and Enhancer Function

The adaptive immune response in vertebrates depends on the expression of antigen-specific receptors in lymphocytes. T-cell receptor (TCR) gene expression is exquisitely regulated during thymocyte development to drive the generation of αβ and γδ T lymphocytes. The TCRα, TCRβ, TCRγ, and TCRδ genes exist in two different configurations, unrearranged and rearranged. A correctly rearranged configuration is required for expression of a functional TCR chain. TCRs can take the form of one of three possible heterodimers, pre-TCR, TCRαβ, or TCRγδ which drive thymocyte maturation into αβ or γδ T lymphocytes. To pass from an unrearranged to a rearranged configuration, global and local three dimensional (3D) chromatin changes must occur during thymocyte development to regulate gene segment accessibility for V(D)J recombination. During this process, enhancers play a critical role by modifying the chromatin conformation and triggering noncoding germline transcription that promotes the recruitment of the recombination machinery. The different signaling that thymocytes receive during their development controls enhancer activity. Here, we summarize the dynamics of long-distance interactions established through chromatin regulatory elements that drive transcription and V(D)J recombination and how different signaling pathways are orchestrated to regulate the activity of enhancers to precisely control TCR gene expression during T-cell maturation.

Microanatomy of the Human Atherosclerotic Plaque by Single-Cell Transcriptomics

RATIONALE

Atherosclerotic lesions are known for their cellular heterogeneity, yet the molecular complexity within the cells of human plaques has not been fully assessed.

OBJECTIVE

Using single-cell transcriptomics and chromatin accessibility, we gained a better understanding of the pathophysiology underlying human atherosclerosis.

METHODS AND RESULTS

We performed single-cell RNA and single-cell ATAC sequencing on human carotid atherosclerotic plaques to define the cells at play and determine their transcriptomic and epigenomic characteristics. We identified 14 distinct cell populations including endothelial cells, smooth muscle cells, mast cells, B cells, myeloid cells, and T cells and identified multiple cellular activation states and suggested cellular interconversions. Within the endothelial cell population, we defined subsets with angiogenic capacity plus clear signs of endothelial to mesenchymal transition. CD4+ and CD8+ T cells showed activation-based subclasses, each with a gradual decline from a cytotoxic to a more quiescent phenotype. Myeloid cells included 2 populations of proinflammatory macrophages showing IL (interleukin) 1B or TNF (tumor necrosis factor) expression as well as a foam cell-like population expressing TREM2 (triggering receptor expressed on myeloid cells 2) and displaying a fibrosis-promoting phenotype. ATACseq data identified specific transcription factors associated with the myeloid subpopulation and T cell cytokine profiles underlying mutual activation between both cell types. Finally, cardiovascular disease susceptibility genes identified using public genome-wide association studies data were particularly enriched in lesional macrophages, endothelial, and smooth muscle cells.

CONCLUSIONS

This study provides a transcriptome-based cellular landscape of human atherosclerotic plaques and highlights cellular plasticity and intercellular communication at the site of disease. This detailed definition of cell communities at play in atherosclerosis will facilitate cell-based mapping of novel interventional targets with direct functional relevance for the treatment of human disease.

Implication of a new function of human tDNAs in chromatin organization

Transfer RNA genes (tDNAs) are essential genes that encode tRNAs in all species. To understand new functions of tDNAs, other than that of encoding tRNAs, we used ENCODE data to examine binding characteristics of transcription factors (TFs) for all tDNA regions (489 loci) in the human genome. We divided the tDNAs into three groups based on the number of TFs that bound to them. At the two extremes were tDNAs to which many TFs bound (Group 1) and those to which no TFs bound (Group 3). Several TFs involved in chromatin remodeling such as ATF3, EP300 and TBL1XR1 bound to almost all Group 1 tDNAs. Furthermore, almost all Group 1 tDNAs included DNase I hypersensitivity sites and may thus interact with other chromatin regions through their bound TFs, and they showed highly conserved synteny across tetrapods. In contrast, Group 3 tDNAs did not possess these characteristics. These data suggest the presence of a previously uncharacterized function of these tDNAs. We also examined binding of CTCF to tDNAs and their involvement in topologically associating domains (TADs) and lamina-associated domains (LADs), which suggest a new perspective on the evolution and function of tDNAs.

Batf Pioneers the Reorganization of Chromatin in Developing Effector T Cells via Ets1-Dependent Recruitment of Ctcf

The basic leucine zipper transcription factor activating transcription factor-like (Batf) contributes to transcriptional programming of multiple effector T cells and is required for T helper 17 (Th17) and T follicular helper (Tfh) cell development. Here, we examine mechanisms by which Batf initiates gene transcription in developing effector CD4 T cells. We find that, in addition to its pioneering function, Batf controls developmentally regulated recruitment of the architectural factor Ctcf to promote chromatin looping that is associated with lineage-specific gene transcription. The chromatin-organizing actions of Batf are largely dependent on Ets1, which appears to be indispensable for the Batf-dependent recruitment of Ctcf. Moreover, most of the Batf-dependent sites to which Ctcf is recruited lie outside of activating protein-1-interferon regulatory factor (Ap-1-Irf) composite elements (AICEs), indicating that direct involvement of Batf-Irf complexes is not required. These results identify a cooperative role for Batf, Ets1, and Ctcf in chromatin reorganization that underpins the transcriptional programming of effector T cells.

Alternative Enhancer Usage and Targeted Polycomb Marking Hallmark Promoter Choice during T Cell Differentiation

During thymic development and upon peripheral activation, T cells undergo extensive phenotypic and functional changes coordinated by lineage-specific developmental programs. To characterize the regulatory landscape controlling T cell identity, we perform a wide epigenomic and transcriptional analysis of mouse thymocytes and naive CD4 differentiated T helper cells. Our investigations reveal a dynamic putative enhancer landscape, and we could validate many of the enhancers using the high-throughput CapStarr sequencing (CapStarr-seq) approach. We find that genes using multiple promoters display increased enhancer usage, suggesting that apparent "enhancer redundancy" might relate to isoform selection. Furthermore, we can show that two Runx3 promoters display long-range interactions with specific enhancers. Finally, our analyses suggest a novel function for the PRC2 complex in the control of alternative promoter usage. Altogether, our study has allowed for the mapping of an exhaustive set of active enhancers and provides new insights into their function and that of PRC2 in controlling promoter choice during T cell differentiation.

Diversity and Emerging Roles of Enhancer RNA in Regulation of Gene Expression and Cell Fate

Enhancers are cis-regulatory elements in the genome that cooperate with promoters to control target gene transcription. Unlike promoters, enhancers are not necessarily adjacent to target genes and can exert their functions regardless of enhancer orientations, positions and spatial segregations from target genes. Thus, for a long time, the question as to how enhancers act in a temporal and spatial manner attracted considerable attention. The recent discovery that enhancers are also abundantly transcribed raises interesting questions about the exact roles of enhancer RNA (eRNA) in gene regulation. In this review, we highlight the process of enhancer transcription and the diverse features of eRNA. We review eRNA functions, which include enhancer-promoter looping, chromatin modifying, and transcription regulating. As eRNA are transcribed from active enhancers, they exhibit tissue and lineage specificity, and serve as markers of cell state and function. Finally, we discuss the unique relationship between eRNA and super enhancers in phase separation wherein eRNA may contribute significantly to cell fate decisions.

Programming for T-lymphocyte fates: modularity and mechanisms

In this review, Rothenburg discusses the gene regulatory network and chromatin-based kinetic constraints that determine activities of transcription factors in the primary establishment of T-cell identity.

T-cell development in mammals is a model for lineage choice and differentiation from multipotent stem cells. Although T-cell fate choice is promoted by signaling in the thymus through one dominant pathway, the Notch pathway, it entails a complex set of gene regulatory network and chromatin state changes even before the cells begin to express their signature feature, the clonal-specific T-cell receptors (TCRs) for antigen. This review distinguishes three developmental modules for T-cell development, which correspond to cell type specification, TCR expression and selection, and the assignment of cells to different effector types. The first is based on transcriptional regulatory network events, the second is dominated by somatic gene rearrangement and mutation and cell selection, and the third corresponds to establishing a poised state of latent regulator priming through an unknown mechanism. Interestingly, in different lineages, the third module can be deployed at variable times relative to the completion of the first two modules. This review focuses on the gene regulatory network and chromatin-based kinetic constraints that determine activities of transcription factors TCF1, GATA3, PU.1, Bcl11b, Runx1, and E proteins in the primary establishment of T-cell identity.

An integrated transcriptional switch at the β-selection checkpoint determines T cell survival, development and leukaemogenesis

In T cell development, a pivotal decision-making stage, termed β-selection, integrates a TCRβ checkpoint to coordinate survival, proliferation and differentiation to an αβ T cell. Here, we review how transcriptional regulation coordinates fate determination in early T cell development to enable β-selection. Errors in this transcription control can trigger T cell acute lymphoblastic leukaemia. We describe how the β-selection checkpoint goes awry in leukaemic transformation.

Transcriptional regulation of Hepatic Stellate Cell activation in NASH

Non-alcoholic steatohepatitis (NASH) signified by hepatic steatosis, inflammation, hepatocellular injury, and fibrosis is a growing cause of chronic liver disease, cirrhosis, and hepatocellular carcinoma. Hepatic fibrosis resulting from accumulation of extracellular matrix proteins secreted by hepatic myofibroblasts plays an important role in disease progression. Activated hepatic stellate cells (HSCs) have been identified as the primary source of myofibroblasts in animal models of hepatotoxic liver injury; however, so far HSC activation and plasticity have not been thoroughly investigated in the context of NASH-related fibrogenesis. Here we have determined the time-resolved changes in the HSC transcriptome during development of Western diet- and fructose-induced NASH in mice, a NASH model recapitulating human disease. Intriguingly, HSC transcriptional dynamics are highly similar across disease models pointing to HSC activation as a point of convergence in the development of fibrotic liver disease. Bioinformatic interrogation of the promoter sequences of activated genes combined with loss-of-function experiments indicates that the transcriptional regulators ETS1 and RUNX1 act as drivers of NASH-associated HSC plasticity. Taken together, our results implicate HSC activation and transcriptional plasticity as key aspects of NASH pathophysiology.

Genetic determinants of co-accessible chromatin regions in activated T cells across humans

Over 90% of genetic variants associated with complex human traits map to non-coding regions, but little is understood about how they modulate gene regulation in health and disease. One possible mechanism is that genetic variants affect the activity of one or more cis-regulatory elements leading to gene expression variation in specific cell types. To identify such cases, we analyzed ATAC-seq and RNA-seq profiles from stimulated primary CD4⁺ T cells in up to 105 healthy donors. We found that regions of accessible chromatin (ATAC-peaks) are co-accessible at kilobase and megabase resolution, consistent with the three-dimensional chromatin organization measured by in situ Hi-C in T cells. Fifteen percent of genetic variants located within ATAC-peaks affected the accessibility of the corresponding peak (local-ATAC-QTLs). Local-ATAC-QTLs have the largest effects on co-accessible peaks, are associated with gene expression and are enriched for autoimmune disease variants. Our results provide insights into how natural genetic variants modulate cis-regulatory elements, in isolation or in concert, to influence gene expression.

Cytokines, Transcription Factors, and the Initiation of T-Cell Development

Multipotent blood progenitor cells migrate into the thymus and initiate the T-cell differentiation program. T-cell progenitor cells gradually acquire T-cell characteristics while shedding their multipotentiality for alternative fates. This process is supported by extracellular signaling molecules, including Notch ligands and cytokines, provided by the thymic microenvironment. T-cell development is associated with dynamic change of gene regulatory networks of transcription factors, which interact with these environmental signals. Together with Notch or pre-T-cell-receptor (TCR) signaling, cytokines always control proliferation, survival, and differentiation of early T cells, but little is known regarding their cross talk with transcription factors. However, recent results suggest ways that cytokines expressed in distinct intrathymic niches can specifically modulate key transcription factors. This review discusses how stage-specific roles of cytokines and transcription factors can jointly guide development of early T cells.

Stochastics of Cellular Differentiation Explained by Epigenetics: The Case of T-Cell Differentiation and Functional Plasticity

Epigenetic marks including histone modifications and DNA methylation are associated with the regulation of gene expression and activity. In addition, an increasing number of non-coding RNAs with regulatory activity on gene expression have been identified. Alongside, technological advancements allow for the analysis of these mechanisms with high resolution up to the single-cell level. For instance, the assay for transposase-accessible chromatin using sequencing (ATAC-seq) simultaneously probes for chromatin accessibility and nucleosome positioning. Thus, it provides information on two levels of epigenetic regulation. Development and differentiation of T cells into functional subset cells including memory T cells are dynamic processes driven by environmental signals. Here, we briefly review the current knowledge of how epigenetic regulation contributes to subset specification, differentiation and memory development in T cells. Specifically, we focus on epigenetic mechanisms differentially active in the two distinct T cell populations expressing αβ or γδ T cell receptors. We also discuss examples of epigenetic alterations of T cells in autoimmune diseases. DNA methylation and histone acetylation are subject to modification by several classes of 'epigenetic modifiers', some of which are in clinical use or in preclinical development. Therefore, we address the impact of some epigenetic modifiers on T-cell activation and differentiation, and discuss possible synergies with T cell-based immunotherapeutic strategies.

Activation of Mouse Tcrb: Uncoupling RUNX1 Function from Its Cooperative Binding with ETS1

T lineage commitment requires the coordination of key transcription factors (TFs) in multipotent progenitors that transition them away from other lineages and cement T cell identity. Two important TFs for the multipotent progenitors to T lineage transition are RUNX1 and ETS1, which bind cooperatively to composite sites throughout the genome, especially in regulatory elements for genes involved in T lymphopoiesis. Activation of the TCR β (Tcrb) locus in committed thymocytes is a critical process for continued development of these cells, and is mediated by an enhancer, Eβ, which harbors two RUNX-ETS composite sites. An outstanding issue in understanding T cell gene expression programs is whether RUNX1 and ETS1 have independent functions in enhancer activation that can be dissected from cooperative binding. We now show that RUNX1 is sufficient to activate the endogenous mouse Eβ element and its neighboring 25 kb region by independently tethering this TF without coincidental ETS1 binding. Moreover, RUNX1 is sufficient for long-range promoter-Eβ looping, nucleosome clearance, and robust transcription throughout the Tcrb recombination center, spanning both DβJβ clusters. We also find that a RUNX1 domain, termed the negative regulatory domain for DNA binding, can compensate for the loss of ETS1 binding at adjacent sites. Thus, we have defined independent roles for RUNX1 in the activation of a T cell developmental enhancer, as well as its ability to mediate specific changes in chromatin landscapes that accompany long-range induction of recombination center promoters.

Chromatin Immunoprecipitation (ChIP) in Mouse T-cell Lines

Signaling pathways regulate gene expression programs via the modulation of the chromatin structure at different levels, such as by post-translational modifications (PTMs) of histone tails, the exchange of canonical histones with histone variants, and nucleosome eviction. Such regulation requires the binding of signal-sensitive transcription factors (TFs) that recruit chromatin-modifying enzymes at regulatory elements defined as enhancers. Understanding how signaling cascades regulate enhancer activity requires a comprehensive analysis of the binding of TFs, chromatin modifying enzymes, and the occupancy of specific histone marks and histone variants. Chromatin immunoprecipitation (ChIP) assays utilize highly specific antibodies to immunoprecipitate specific protein/DNA complexes. The subsequent analysis of the purified DNA allows for the identification the region occupied by the protein recognized by the antibody. This work describes a protocol to efficiently perform ChIP of histone proteins in a mature mouse T-cell line. The presented protocol allows for the performance of ChIP assays in a reasonable timeframe and with high reproducibility.

Signatures of DNA target selectivity by ETS transcription factors

The ETS family of transcription factors is a functionally heterogeneous group of gene regulators that share a structurally conserved, eponymous DNA-binding domain. DNA target specificity derives from combinatorial interactions with other proteins as well as intrinsic heterogeneity among ETS domains. Emerging evidence suggests molecular hydration as a fundamental feature that defines the intrinsic heterogeneity in DNA target selection and susceptibility to epigenetic DNA modification. This perspective invokes novel hypotheses in the regulation of ETS proteins in physiologic osmotic stress, their pioneering potential in heterochromatin, and the effects of passive and pharmacologic DNA demethylation on ETS regulation.

Two possible modes of pioneering associated with combinations of H2A.Z and p300/CBP at nucleosome-occupied enhancers

Pioneer transcription factors are defined by their ability to bind nucleosome-occupied regions. Here, we discuss the properties of nucleosomes bound by pioneers at enhancer regions. We describe how select pioneers bind nucleosome-occupied or -depleted enhancer sites. Importantly, by revisiting and expanding existing data sets, we show differential H2A.Z and p300/CBP association at bound enhancers, highlighting two possible pioneering modes.

Hematopoiesis and T-cell specification as a model developmental system

The pathway to generate T cells from hematopoietic stem cells guides progenitors through a succession of fate choices while balancing differentiation progression against proliferation, stage to stage. Many elements of the regulatory system that controls this process are known, but the requirement for multiple, functionally distinct transcription factors needs clarification in terms of gene network architecture. Here, we compare the features of the T-cell specification system with the rule sets underlying two other influential types of gene network models: first, the combinatorial, hierarchical regulatory systems that generate the orderly, synchronized increases in complexity in most invertebrate embryos; second, the dueling 'master regulator' systems that are commonly used to explain bistability in microbial systems and in many fate choices in terminal differentiation. The T-cell specification process shares certain features with each of these prevalent models but differs from both of them in central respects. The T-cell system is highly combinatorial but also highly dose-sensitive in its use of crucial regulatory factors. The roles of these factors are not always T-lineage-specific, but they balance and modulate each other's activities long before any mutually exclusive silencing occurs. T-cell specification may provide a new hybrid model for gene networks in vertebrate developmental systems.

Impact of the gut microbiota on enhancer accessibility in gut intraepithelial lymphocytes

The gut microbiota impacts many aspects of host biology including immune function. One hypothesis is that microbial communities induce epigenetic changes with accompanying alterations in chromatin accessibility, providing a mechanism that allows a community to have sustained host effects even in the face of its structural or functional variation. We used Assay for Transposase-Accessible Chromatin with high-throughput sequencing (ATAC-seq) to define chromatin accessibility in predicted enhancer regions of intestinal αβ⁺ and γδ⁺ intraepithelial lymphocytes purified from germ-free mice, their conventionally raised (CONV-R) counterparts, and mice reared germ free and then colonized with CONV-R gut microbiota at the end of the suckling-weaning transition. Characterizing genes adjacent to traditional enhancers and super-enhancers revealed signaling networks, metabolic pathways, and enhancer-associated transcription factors affected by the microbiota. Our results support the notion that epigenetic modifications help define microbial community-affiliated functional features of host immune cell lineages.