New roles for DNA cytosine modification, eRNA, anchors, and superanchors in developing B cell progenitors

The Function of E2A in B-Cell Development

Helix-loop-helix (HLH) transcription factors (TFs) play a key role in various cellular differentiation and function through the regulation of enhancer activity. E2A, a member of the mammalian E-protein family (class I HLH protein), is well known to play an important role in hematopoiesis, especially in adaptive lymphocyte development. E2A instructs B- and T-cell lineage development through the regulation of enhancer activity for B- or T-cell signature gene expression, including Rag1 and Rag2 (Rag1/2) genes. In this chapter, we mainly focus on the function of E2A in B-cell development and on the roles of E2A in establishing the enhancer landscape through the recruitment of EP300/KAT3B, chromatin remodeling complex, mediator, cohesion, and TET proteins. Finally, we demonstrate how E2A orchestrates the assembly of the Rag1/2 gene super-enhancer (SE) formation by changing the chromatin conformation across the Rag gene locus.

Intra- and interchromosomal contact mapping reveals the Igh locus has extensive conformational heterogeneity and interacts with B-lineage genes

To produce a diverse antibody repertoire, immunoglobulin heavy-chain (Igh) loci undergo large-scale alterations in structure to facilitate juxtaposition and recombination of spatially separated variable (VH), diversity (DH), and joining (JH) genes. These chromosomal alterations are poorly understood. Uncovering their patterns shows how chromosome dynamics underpins antibody diversity. Using tiled Capture Hi-C, we produce a comprehensive map of chromatin interactions throughout the 2.8-Mb Igh locus in progenitor B cells. We find that the Igh locus folds into semi-rigid subdomains and undergoes flexible looping of the VH genes to its 3' end, reconciling two views of locus organization. Deconvolution of single Igh locus conformations using polymer simulations identifies thousands of different structures. This heterogeneity may underpin the diversity of V(D)J recombination events. All three immunoglobulin loci also participate in a highly specific, developmentally regulated network of interchromosomal interactions with genes encoding B cell-lineage factors. This suggests a model of interchromosomal coordination of B cell development.

Insight into the regulatory mechanism of dynamic chromatin 3D interactions during cardiomyocyte differentiation in human

Cardiogenesis is an extremely complicated process involved with DNA regulatory elements, and trans factors regulate gene expression pattern spatiotemporally. Enhancers, as the well-known DNA elements, activate target gene expression by transcription factors (TFs) occupied to organize dynamic three-dimensional (3D) interactions, which when affected or interrupted might cause heart defects or diseases. In this study, we integrated transcriptome, 3D genome, and regulatome to reorganize the global 3D genome in cardiomyogenesis, showing a gradually decreased trend of both chromatin interactions and topological associating domains (TADs) during cardiomyocyte differentiation. And almost all of the chromatin interactions occurred within the same or between adjacent TADs involved with enhancers, indicating that dynamical rewiring of enhancer-related chromatin interactions in the continuous expansive TADs is closely correlated to cardiogenesis. Moreover, we found stage-specific interactions activate stage-specific expression to be involved within corresponding biological functions, and the stage-specific combined regulations of enhancers and binding TFs form connected networks to control stage-specific expression and biological processes, which promote cardiomyocyte differentiation. Finally, we identified markers based on regulatory networks, which might drive cardiac development. This study demonstrates the power of enhancer interactome combined with active TFs to reveal insights into transcriptional regulatory networks during cardiomyogenesis.

Contribution of the IGCR1 regulatory element and the 3'Igh CTCF-binding elements to regulation of Igh V(D)J recombination

Immunoglobulin heavy chain variable region exons are assembled in progenitor-B cells, from VH, D, and JH gene segments located in separate clusters across the Igh locus. RAG endonuclease initiates V(D)J recombination from a JH-based recombination center (RC). Cohesin-mediated extrusion of upstream chromatin past RC-bound RAG presents Ds for joining to JHs to form a DJH-RC. Igh has a provocative number and organization of CTCF-binding elements (CBEs) that can impede loop extrusion. Thus, Igh has two divergently oriented CBEs (CBE1 and CBE2) in the IGCR1 element between the VH and D/JH domains, over 100 CBEs across the VH domain convergent to CBE1, and 10 clustered 3'Igh-CBEs convergent to CBE2 and VH CBEs. IGCR1 CBEs segregate D/JH and VH domains by impeding loop extrusion-mediated RAG-scanning. Downregulation of WAPL, a cohesin unloader, in progenitor-B cells neutralizes CBEs, allowing DJH-RC-bound RAG to scan the VH domain and perform VH-to-DJH rearrangements. To elucidate potential roles of IGCR1-based CBEs and 3'Igh-CBEs in regulating RAG-scanning and elucidate the mechanism of the ordered transition from D-to-JH to VH-to-DJH recombination, we tested effects of inverting and/or deleting IGCR1 or 3'Igh-CBEs in mice and/or progenitor-B cell lines. These studies revealed that normal IGCR1 CBE orientation augments RAG-scanning impediment activity and suggest that 3'Igh-CBEs reinforce ability of the RC to function as a dynamic loop extrusion impediment to promote optimal RAG scanning activity. Finally, our findings indicate that ordered V(D)J recombination can be explained by a gradual WAPL downregulation mechanism in progenitor-B cells as opposed to a strict developmental switch.

Contribution of the IGCR1 regulatory element and the 3 'Igh CBEs to Regulation of Igh V(D)J Recombination

Immunoglobulin heavy chain variable region exons are assembled in progenitor-B cells, from V _H , D, and J _H gene segments located in separate clusters across the Igh locus. RAG endonuclease initiates V(D)J recombination from a J _H -based recombination center (RC). Cohesin-mediated extrusion of upstream chromatin past RC-bound RAG presents Ds for joining to J _H s to form a DJ _H -RC. Igh has a provocative number and organization of CTCF-binding-elements (CBEs) that can impede loop extrusion. Thus, Igh has two divergently oriented CBEs (CBE1 and CBE2) in the IGCR1 element between the V _H and D/J _H domains, over 100 CBEs across the V _H domain convergent to CBE1, and 10 clustered 3' Igh -CBEs convergent to CBE2 and V _H CBEs. IGCR1 CBEs segregate D/J _H and V _H domains by impeding loop extrusion-mediated RAG-scanning. Down-regulation of WAPL, a cohesin unloader, in progenitor-B cells neutralizes CBEs, allowing DJ _H -RC-bound RAG to scan the VH domain and perform VH-to-DJH rearrangements. To elucidate potential roles of IGCR1-based CBEs and 3' Igh -CBEs in regulating RAG-scanning and elucidate the mechanism of the "ordered" transition from D-to-J _H to V _H -to-DJ _H recombination, we tested effects of deleting or inverting IGCR1 or 3' Igh -CBEs in mice and/or progenitor-B cell lines. These studies revealed that normal IGCR1 CBE orientation augments RAG-scanning impediment activity and suggest that 3' Igh -CBEs reinforce ability of the RC to function as a dynamic loop extrusion impediment to promote optimal RAG scanning activity. Finally, our findings indicate that ordered V(D)J recombination can be explained by a gradual WAPL down-regulation mechanism in progenitor B cells as opposed to a strict developmental switch.

SIGNIFICANCE STATEMENT

To counteract diverse pathogens, vertebrates evolved adaptive immunity to generate diverse antibody repertoires through a B lymphocyte-specific somatic gene rearrangement process termed V(D)J recombination. Tight regulation of the V(D)J recombination process is vital to generating antibody diversity and preventing off-target activities that can predispose the oncogenic translocations. Recent studies have demonstrated V(D)J rearrangement is driven by cohesin-mediated chromatin loop extrusion, a process that establishes genomic loop domains by extruding chromatin, predominantly, between convergently-oriented CTCF looping factor-binding elements (CBEs). By deleting and inverting CBEs within a critical antibody heavy chain gene locus developmental control region and a loop extrusion chromatin-anchor at the downstream end of this locus, we reveal how these elements developmentally contribute to generation of diverse antibody repertoires.

Profiling the epigenetic landscape of the antigen receptor repertoire: the missing epi-immunogenomics data

High-resolution sequencing methods that capture the epigenetic landscape within the T cell receptor (TCR) gene loci are pivotal for a fundamental understanding of the epigenetic regulatory mechanisms of the TCR repertoire. In our opinion, filling the gaps in our understanding of the epigenetic mechanisms regulating the TCR repertoire will benefit the development of strategies that can modulate the TCR repertoire composition by leveraging the dynamic nature of epigenetic modifications.

An Igh distal enhancer modulates antigen receptor diversity by determining locus conformation

The mouse Igh locus is organized into a developmentally regulated topologically associated domain (TAD) that is divided into subTADs. Here we identify a series of distal V_H enhancers (E_VHs) that collaborate to configure the locus. E_VHs engage in a network of long-range interactions that interconnect the subTADs and the recombination center at the D_HJ_H gene cluster. Deletion of E_VH1 reduces V gene rearrangement in its vicinity and alters discrete chromatin loops and higher order locus conformation. Reduction in the rearrangement of the V_H11 gene used in anti-PtC responses is a likely cause of the observed reduced splenic B1 B cell compartment. E_VH1 appears to block long-range loop extrusion that in turn contributes to locus contraction and determines the proximity of distant V_H genes to the recombination center. E_VH1 is a critical architectural and regulatory element that coordinates chromatin conformational states that favor V(D)J rearrangement.

Roles of super enhancers and enhancer RNAs in skeletal muscle development and disease

Skeletal muscle development is a multistep biological process regulated by a variety of myogenic regulatory factors, including MyoG, MyoD, Myf5, and Myf6 (also known as MRF4), as well as members of the FoxO subfamily. Differentiation and regeneration during skeletal muscle myogenesis contribute to the physiological function of muscles. Super enhancers (SEs) and enhancer RNAs (eRNAs) are involved in the regulation of development and diseases. Few studies have identified the roles of SEs and eRNAs in muscle development and pathophysiology. To develop approaches to enhance skeletal muscle mass and function, a more comprehensive understanding of the key processes underlying muscular diseases is needed. In this review, we summarize the roles of SEs and eRNAs in muscle development and disease through affecting of DNA methylation, FoxO subfamily, RAS-MEK signaling, chromatin modifications and accessibility, MyoD and cis regulating target genes. The summary could inform strategies to increase muscle mass and treat muscle-related diseases.

Locus architecture and RAG scanning determine antibody diversity

Diverse mammalian antibody repertoires are produced via distant genomic contacts involving immunoglobulin Igh variable (V), diversity (D), and joining (J) gene segments and result in V(D)J recombination. How such interactions determine V gene usage remains unclear. The recombination-activating gene (RAG) chromatin scanning model posits that RAG recombinase bound to the recombination center (RC) linearly tracks along chromatin by means of cohesin-mediated loop extrusion; a proposition supported by cohesin depletion studies. A mechanistic role for chromatin loop extrusion has also been implicated for Igh locus contraction. In this opinion, we provide perspective on how loop extrusion interfaces with the 3D conformation of the Igh locus and newly identified enhancers that regionally regulate VH gene usage during V(D)J recombination, shaping the preselected repertoire.

RANKL-responsive epigenetic mechanism reprograms macrophages into bone-resorbing osteoclasts

Monocyte/macrophage lineage cells are highly plastic and can differentiate into various cells under different environmental stimuli. Bone-resorbing osteoclasts are derived from the monocyte/macrophage lineage in response to receptor activator of NF-κB ligand (RANKL). However, the epigenetic signature contributing to the fate commitment of monocyte/macrophage lineage differentiation into human osteoclasts is largely unknown. In this study, we identified RANKL-responsive human osteoclast-specific superenhancers (SEs) and SE-associated enhancer RNAs (SE-eRNAs) by integrating data obtained from ChIP-seq, ATAC-seq, nuclear RNA-seq and PRO-seq analyses. RANKL induced the formation of 200 SEs, which are large clusters of enhancers, while suppressing 148 SEs in macrophages. RANKL-responsive SEs were strongly correlated with genes in the osteoclastogenic program and were selectively increased in human osteoclasts but marginally presented in osteoblasts, CD4+ T cells, and CD34+ cells. In addition to the major transcription factors identified in osteoclasts, we found that BATF binding motifs were highly enriched in RANKL-responsive SEs. The depletion of BATF1/3 inhibited RANKL-induced osteoclast differentiation. Furthermore, we found increased chromatin accessibility in SE regions, where RNA polymerase II was significantly recruited to induce the extragenic transcription of SE-eRNAs, in human osteoclasts. Knocking down SE-eRNAs in the vicinity of the NFATc1 gene diminished the expression of NFATc1, a major regulator of osteoclasts, and osteoclast differentiation. Inhibiting BET proteins suppressed the formation of some RANKL-responsive SEs and NFATc1-associated SEs, and the expression of SE-eRNA:NFATc1. Moreover, SE-eRNA:NFATc1 was highly expressed in the synovial macrophages of rheumatoid arthritis patients exhibiting high-osteoclastogenic potential. Our genome-wide analysis revealed RANKL-inducible SEs and SE-eRNAs as osteoclast-specific signatures, which may contribute to the development of osteoclast-specific therapeutic interventions.

The role of chromatin loop extrusion in antibody diversification

Cohesin mediates chromatin loop formation across the genome by extruding chromatin between convergently oriented CTCF-binding elements. Recent studies indicate that cohesin-mediated loop extrusion in developing B cells presents immunoglobulin heavy chain (Igh) variable (V), diversity (D) and joining (J) gene segments to RAG endonuclease through a process referred to as RAG chromatin scanning. RAG initiates V(D)J recombinational joining of these gene segments to generate the large number of different Igh variable region exons that are required for immune responses to diverse pathogens. Antigen-activated mature B cells also use chromatin loop extrusion to mediate the synapsis, breakage and end joining of switch regions flanking Igh constant region exons during class-switch recombination, which allows for the expression of different antibody constant region isotypes that optimize the functions of antigen-specific antibodies to eliminate pathogens. Here, we review recent advances in our understanding of chromatin loop extrusion during V(D)J recombination and class-switch recombination at the Igh locus.

The HDAC7-TET2 epigenetic axis is essential during early B lymphocyte development

Correct B cell identity at each stage of cellular differentiation during B lymphocyte development is critically dependent on a tightly controlled epigenomic landscape. We previously identified HDAC7 as an essential regulator of early B cell development and its absence leads to a drastic block at the pro-B to pre-B cell transition. More recently, we demonstrated that HDAC7 loss in pro-B-ALL in infants associates with a worse prognosis. Here we delineate the molecular mechanisms by which HDAC7 modulates early B cell development. We find that HDAC7 deficiency drives global chromatin de-condensation, histone marks deposition and deregulates other epigenetic regulators and mobile elements. Specifically, the absence of HDAC7 induces TET2 expression, which promotes DNA 5-hydroxymethylation and chromatin de-condensation. HDAC7 deficiency also results in the aberrant expression of microRNAs and LINE-1 transposable elements. These findings shed light on the mechanisms by which HDAC7 loss or misregulation may lead to B cell–based hematological malignancies.

Super-Enhancers, Phase-Separated Condensates, and 3D Genome Organization in Cancer

3D chromatin organization plays an important role in transcription regulation and gene expression. The 3D genome is highly maintained by several architectural proteins, such as CTCF, Yin Yang 1, and cohesin complex. This structural organization brings regulatory DNA elements in close proximity to their target promoters. In this review, we discuss the 3D chromatin organization of super-enhancers and their relationship to phase-separated condensates. Super-enhancers are large clusters of DNA elements. They can physically contact with their target promoters by chromatin looping during transcription. Multiple transcription factors can bind to enhancer and promoter sequences and recruit a complex array of transcriptional co-activators and RNA polymerase II to effect transcriptional activation. Phase-separated condensates of transcription factors and transcriptional co-activators have been implicated in assembling the transcription machinery at particular enhancers. Cancer cells can hijack super-enhancers to drive oncogenic transcription to promote cell survival and proliferation. These dysregulated transcriptional programs can cause cancer cells to become highly dependent on transcriptional regulators, such as Mediator and BRD4. Moreover, the expression of oncogenes that are driven by super-enhancers is sensitive to transcriptional perturbation and often occurs in phase-separated condensates, supporting therapeutic rationales of targeting SE components, 3D genome organization, or dysregulated condensates in cancer.

Enhancer methylation dynamics drive core transcriptional regulatory circuitry in pan-cancer

Accumulating evidence has demonstrated that enhancer methylation has strong and dynamic regulatory effects on gene expression. Some transcription factors (TFs) can auto- and cross-regulate in a feed-forward manner, and cooperate with their enhancers to form core transcriptional regulatory circuitries (CRCs). However, the elaborated regulatory mechanism between enhancer methylation and CRC remains the tip of the iceberg. Here, we revealed that DNA methylation could drive the tissue-specific enhancer basal transcription and target gene expression in human cancers. By integrating methylome, transcriptome, and 3D genomic data, we identified enhancer methylation triplets (enhancer methylation-enhancer transcription-target gene expression) and dissected potential regulatory patterns within them. Moreover, we observed that cancer-specific core TFs regulated by enhancers were able to shape their enhancer methylation forming the enhancer methylation-driven CRCs (emCRCs). Further parsing of clinical implications showed rewired emCRCs could serve as druggable targets and prognostic risk markers. In summary, the integrative analysis of enhancer methylation regulome would facilitate portraying the cancer epigenomics landscape and developing the epigenetic anti-cancer approaches.

DNMT3A and TET2; Potential Estimates of Generic DNA Methylation in Children and Adolescents with Obesity; Relation to Metabolic Dysregulation

INTRODUCTION

The role of DNA methylation in metabolic dysregulation is emerging. However, the functional role of methylation in obesity and metabolic dysregulation is poorly understood.

AIM

The aim of this study was to compare DNA methyltransferase-3A (DNMT3A) and ten-eleven translocase-2 (TET2) levels in children and adolescents with obesity to normal-weighed children and adolescents and to correlate them to various metabolic parameters.

METHODS

Fifty children and adolescents with obesity were compared to 50 matched normal-weighed children and adolescents. Participants underwent assessment for anthropometric measurements, Tanner staging, acanthosis nigricans, and mean blood pressure percentile on three different occasions. TET2, DNMT3A, fasting lipids, and insulin were measured with calculation of the homeostatic model assessment insulin resistance (HOMA-IR).

RESULTS

The median BMI SDS of the studied children and adolescents with obesity was 3.40, their mean TET2 was 178.40 ng/mL, and their mean DNMT3A was 2.18 ng/mL. TET2 is significantly lower (p = 0.009), while DNMT3A is significantly higher (p < 0.001) in children and adolescents with obesity than controls. Children and adolescents with obesity and insulin resistance have significantly lower TET2 (p = 0.012) and significantly higher DNMT3A (p = 0.013) than those without insulin resistance. Diastolic blood pressure percentile and HOMA-IR are positively correlated to DNMT3A (p < 0.001) and negatively correlated to TET-2 (p < 0.001). Multivariate logistic regression analysis revealed that TET2 and DNMT3A are independently associated with diastolic blood pressure percentile (p = 0.03 and p = 0.014, respectively) and HOMA-IR (p = 0.003 and p = 0.001, respectively).

CONCLUSIONS

Children and adolescents with obesity have significantly higher DNMT3A and significantly lower TET2 than controls. This is more evident in those having insulin resistance than those without. DNMT3A and TET2 are independently associated with systemic hypertension and insulin resistance in children with obesity.

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.

Exploring the role of eRNA in regulating gene expression

eRNAs as the products of enhancers can regulate gene expression via various possible ways, but which regulation way is more reasonable is debatable in biology, and in particular, how eRNAs impact gene expression remains unclear. Here we introduce a mechanistic model of gene expression to address these issues. This model considers three possible regulation ways of eRNA: Type-I by which eRNA regulates transcriptional activity by facilitating the formation of enhancer-promoter (E-P) loop, Type-II by which eRNA directly promotes the mRNA production rate, and mixed regulation (i.e., the combination of Type-I and Type-II). We show that with the increase of the E-P loop length, mRNA distribution can transition from unimodality to bimodality or vice versa in all the three regulation cases. However, in contrast to the other two regulations, Type-II regulation can lead to the highest mean mRNA level and the lowest mRNA noise, independent of the E-P loop length. These results would not only reveal the essential mechanism of how eRNA regulates gene expression, but also imply a new mechanism for phenotypic switching, namely the E-P loop can induce phenotypic switching.

Loop extrusion promotes an alternate pathway for isotype switching

Class-switch recombination (CSR) involves replacement of the Cμ constant region with another downstream C_H region. CSR is initiated by activation-induced cytidine deaminase (AID)-mediated DNA breaks that are targeted to transcriptionally active switch (S) regions. S region promoters (Prs) direct synapsis by associating with the Eμ and 3′Eα enhancers that jointly anchor a chromatin loop. We report that asymmetric loop extrusion allows 3′Eα to track along the locus and form Pr-Pr-E interactions that mediate CSR between downstream S regions, followed by switching to donor Sμ. This alternative pathway bypasses sequential switching and creates immunoglobulin (Ig)E⁺ B cells in the absence of IgG1 expression. Based on the analysis of diagnostic CSR products in B cell subsets, we identify a BCR-negative cell intermediate that is pivotal to efficient CSR.

Shen et al. report that 3′Eα tracks along the Igh locus via unidirectional loop extrusion to form germline transcript promoter (Pr)-Pr-E interactions that mediate an alternative CSR pathway. B cell intermediates of CSR are identified, which are AID-dependent, surface BCR-negative, and in the G₁ phase of the cell cycle.

DNA methylation signatures reveal that distinct combinations of transcription factors specify human immune cell epigenetic identity

Epigenetic reprogramming underlies specification of immune cell lineages, but patterns that uniquely define immune cell types and the mechanisms by which they are established remain unclear. Here, we identified lineage-specific DNA methylation signatures of six immune cell types from human peripheral blood and determined their relationship to other epigenetic and transcriptomic patterns. Sites of lineage-specific hypomethylation were associated with distinct combinations of transcription factors in each cell type. By contrast, sites of lineage-specific hypermethylation were restricted mostly to adaptive immune cells. PU.1 binding sites were associated with lineage-specific hypo- and hypermethylation in different cell types, suggesting that it regulates DNA methylation in a context-dependent manner. These observations indicate that innate and adaptive immune lineages are specified by distinct epigenetic mechanisms via combinatorial and context-dependent use of key transcription factors. The cell-specific epigenomics and transcriptional patterns identified serve as a foundation for future studies on immune dysregulation in diseases and aging.

Long-Range Control of Class Switch Recombination by Transcriptional Regulatory Elements

Immunoglobulin class switch recombination (CSR) plays a crucial role in adaptive immune responses through a change of the effector functions of antibodies and is triggered by T-cell-dependent as well as T-cell-independent antigens. Signals generated following encounter with each type of antigen direct CSR to different isotypes. At the genomic level, CSR occurs between highly repetitive switch sequences located upstream of the constant gene exons of the immunoglobulin heavy chain locus. Transcription of switch sequences is mandatory for CSR and is induced in a stimulation-dependent manner. Switch transcription takes place within dynamic chromatin domains and is regulated by long-range regulatory elements which promote alignment of partner switch regions in CSR centers. Here, we review recent work and models that account for the function of long-range transcriptional regulatory elements and the chromatin-based mechanisms involved in the control of CSR.

B1a and B2 cells are characterized by distinct CpG modification states at DNMT3A-maintained enhancers

The B1 and B2 lineages of B cells contribute to protection from pathogens in distinct ways. The role of the DNA CpG methylome in specifying these two B-cell fates is still unclear. Here we profile the CpG modifications and transcriptomes of peritoneal B1a and follicular B2 cells, as well as their respective proB cell precursors in the fetal liver and adult bone marrow from wild-type and CD19-Cre Dnmt3a floxed mice lacking DNMT3A in the B lineage. We show that an underlying foundational CpG methylome is stably established during B lineage commitment and is overlaid with a DNMT3A-maintained dynamic methylome that is sculpted in distinct ways in B1a and B2 cells. This dynamic DNMT3A-maintained methylome is composed of novel enhancers that are closely linked to lineage-specific genes. While DNMT3A maintains the methylation state of these enhancers in both B1a and B2 cells, the dynamic methylome undergoes a prominent programmed demethylation event during B1a but not B2 cell development. We propose that the methylation pattern of DNMT3A-maintained enhancers is determined by the coincident recruitment of DNMT3A and TET enzymes, which regulate the developmental expression of B1a and B2 lineage-specific genes.

The Interplay Between Chromatin Architecture and Lineage-Specific Transcription Factors and the Regulation of Rag Gene Expression

Cell type-specific gene expression is driven through the interplay between lineage-specific transcription factors (TFs) and the chromatin architecture, such as topologically associating domains (TADs), and enhancer-promoter interactions. To elucidate the molecular mechanisms of the cell fate decisions and cell type-specific functions, it is important to understand the interplay between chromatin architectures and TFs. Among enhancers, super-enhancers (SEs) play key roles in establishing cell identity. Adaptive immunity depends on the RAG-mediated assembly of antigen recognition receptors. Hence, regulation of the Rag1 and Rag2 (Rag1/2) genes is a hallmark of adaptive lymphoid lineage commitment. Here, we review the current knowledge of 3D genome organization, SE formation, and Rag1/2 gene regulation during B cell and T cell differentiation.

Dynamic 3D Locus Organization and Its Drivers Underpin Immunoglobulin Recombination

A functional adaptive immune system must generate enormously diverse antigen receptor (AgR) repertoires from a limited number of AgR genes, using a common mechanism, V(D)J recombination. The AgR loci are among the largest in the genome, and individual genes must overcome huge spatial and temporal challenges to co-localize with optimum variability. Our understanding of the complex mechanisms involved has increased enormously, due in part to new technologies for high resolution mapping of AgR structure and dynamic movement, underpinning mechanisms, and resulting repertoires. This review will examine these advances using the paradigm of the mouse immunoglobulin heavy chain (Igh) locus. We will discuss the key regulatory elements implicated in Igh locus structure. Recent next generation repertoire sequencing methods have shown that local chromatin state at V genes contribute to recombination efficiency. Next on the multidimensional scale, we will describe imaging studies that provided the first picture of the large-scale dynamic looping and contraction the Igh locus undergoes during recombination. We will discuss chromosome conformation capture (3C)-based technologies that have provided higher resolution pictures of Igh locus structure, including the different models that have evolved. We will consider the key transcription factors (PAX5, YY1, E2A, Ikaros), and architectural factors, CTCF and cohesin, that regulate these processes. Lastly, we will discuss a plethora of recent exciting mechanistic findings. These include Rag recombinase scanning for convergent RSS sequences within DNA loops; identification of Igh loop extrusion, and its putative role in Rag scanning; the roles of CTCF, cohesin and cohesin loading factor, WAPL therein; a new phase separation model for Igh locus compartmentalization. We will draw these together and conclude with some horizon-scanning and unresolved questions.

The transcription factor E2A activates multiple enhancers that drive Rag expression in developing T and B cells

Cell type-specific gene expression is driven by the interplay between lineage-specific transcription factors and cis-regulatory elements to which they bind. Adaptive immunity relies on RAG-mediated assembly of T cell receptor (TCR) and immunoglobulin (Ig) genes. Although Rag1 and Rag2 expression is largely restricted to adaptive lymphoid lineage cells, it remains unclear how Rag gene expression is regulated in a cell lineage-specific manner. Here, we identified three distinct cis-regulatory elements, a T cell lineage-specific enhancer (R-TEn) and the two B cell-specific elements, R1B and R2B By generating mice lacking either R-TEn or R1B and R2B, we demonstrate that these distinct sets of regulatory elements drive the expression of Rag genes in developing T and B cells. What these elements have in common is their ability to bind the transcription factor E2A. By generating a mouse strain that carries a mutation within the E2A binding site of R-TEn, we demonstrate that recruitment of E2A to this site is essential for orchestrating changes in chromatin conformation that drive expression of Rag genes in T cells. By mapping cis-regulatory elements and generating multiple mouse strains lacking distinct enhancer elements, we demonstrate expression of Rag genes in developing T and B cells to be driven by distinct sets of E2A-dependent cis-regulatory modules.

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.

Enhancers as regulators of antigen receptor loci three-dimensional chromatin structure

Enhancers are defined as regulatory elements that control transcription in a cell-type and developmental stage-specific manner. They achieve this by physically interacting with their cognate gene promoters. Significantly, these interactions can occur through long genomic distances since enhancers may not be near their cognate promoters. The optimal coordination of enhancer-regulated transcription is essential for the function and identity of the cell. Although great efforts to fully understand the principles of this type of regulation are ongoing, other potential functions of the long-range chromatin interactions (LRCIs) involving enhancers are largely unexplored. We recently uncovered a new role for enhancer elements in determining the three-dimensional (3D) structure of the immunoglobulin kappa (Igκ) light chain receptor locus suggesting a structural function for these DNA elements. This enhancer-mediated locus configuration shapes the resulting Igκ repertoire. We also propose a role for enhancers as critical components of sub-topologically associating domain (subTAD) formation and nuclear spatial localization.

The fundamental role of chromatin loop extrusion in physiological V(D)J recombination

RAG endonuclease initiates IgH locus (Igh) V(D)J assembly in progenitor (pro)-B cells by joining Ds to J_Hs, before joining upstream V_Hs to DJ_H intermediates¹. In mouse pro-B cells, the CTCF-binding element (CBE)-anchored chromatin loop domain² at the 3’end of Igh contains an internal sub-domain spanning the 5’CBE anchor (IGCR1)³, the D_Hs, and a RAG-bound recombination center (RC)⁴. The RC comprises J_H-proximal D (DQ52), 4 J_Hs, and the intronic enhancer (“iEμ”)⁵. Robust RAG cleavage is restricted to paired V(D)J segments flanked by complementary recombination signal sequences (12RSSs and 23RSSs)⁶. Ds are flanked downstream and upstream by 12RSSs that, respectively, mediate deletional joining with convergently-oriented J_H-23RSSs and V_H-23RSSs⁶. Despite 12/23 compatibility, inversional D to J_H joining via upstream D-12RSSs is rare^7,8. Plasmid-based assays attributed lack of inversional D to J_H joining to sequence-based preference for downstream D-12RSSs⁹, as opposed to putative linear scanning mechanisms^10,11. Given recent findings that RAG linearly scans convergent CBE-anchored chromatin loops^4,12-14, potentially formed by cohesin-mediated loop extrusion^15-18, we revisited a scanning role. Here, we report that J_H-23RSS chromosomal orientation programs RC-bound RAG to linearly scan upstream chromatin in the 3’Igh sub-domain for convergently-oriented D-12RSSs and, thereby, to mediate deletional joining of all Ds, except RC-based DQ52 that joins by a diffusion-related mechanism. In a DQ52-based RC, formed in the absence of J_Hs, RAG bound by the downstream DQ52-RSS scans the downstream constant region exon-containing 3’Igh sub-domain in which scanning can be impeded by targeted nuclease-dead Cas9 (dCas9) binding, by transcription through repetitive Igh switch sequences, and by the 3’Igh CBE-based loop anchor. Notably, each scanning impediment focally increases RAG activity on potential substrate sequences within the impeded region. High resolution mapping of RC chromatin interactions reveals that such focal RAG targeting is associated with corresponding impediments to the loop extrusion process that drives chromatin past RC-bound RAG.

Integrative characterization of G-Quadruplexes in the three-dimensional chromatin structure

DNA molecules are highly compacted in the eukaryotic nucleus where distal regulatory elements reach their targets through three-dimensional chromosomal interactions. G-quadruplexes, stable four-stranded non-canonical DNA structures, can change local chromatin organization through the exclusion of nucleosomes. However, the relationship between G-quadruplexes and higher-order genome organization remains unknown. Here, we found that G-quadruplexes are significantly enriched at boundaries of topological associated domains (TADs). Architectural protein occupancy, which plays critical roles in the formation of TADs, was highly correlated with the content of G-quadruplexes at TAD boundaries. Moreover, adjacent boundaries containing G-quadruplexes frequently interacted with each other because of the high enrichment of architectural protein binding sites. Similar to CCCTC-binding factor (CTCF) binding sites, G-quadruplexes also showed strong insulation ability in the separation of adjacent regions. Additionally, the insulation ability of CTCF binding sites and TAD boundaries was significantly reinforced by G-quadruplexes. Furthermore, G-quadruplex motifs on different strands were associated with the orientation of CTCF binding sites. These findings suggest a potential role for G-quadruplexes in loop extrusion. The enrichment of transcription factor binding sites (TFBSs) around regulatory elements containing G-quadruplexes led to frequent interactions between regulatory elements containing G-quadruplexes. Intriguingly, more than 99% of G-quadruplexes overlapped with TFBSs. The binding sites of CTCF and cohesin proteins were preferentially located surrounding G-quadruplexes. Accordingly, we proposed a new mechanism of long-distance gene regulation in which G-quadruplexes are involved in distal interactions between enhancers and promoters.

Noncoding RNA transcription at enhancers and genome folding in cancer

Changes of nuclear localization of lineage-specific genes from a transcriptionally inert to permissive environment are a crucial step in establishing the identity of a cell. Noncoding RNA transcription-mediated genome folding and activation of target gene expression have been found in a variety of cell types. Noncoding RNA ThymoD (thymocyte differentiation factor) transcription at superenhancers is essential for mouse T-cell lineage commitment. The cessation of ThymoD transcription abolishes transcription-mediated demethylation, recruiting looping factors such as the cohesin complex, CCCTC-binding factor (CTCF), ultimately leading to the phenotype of severe combined immunodeficiency and T-cell leukemia/lymphoma. In this review, we describe the functional role of RNA polymerase II-mediated transcription at enhancers and in genome folding. We also highlight the involvement of faulty activation or suppression of enhancer transcription and enhancer-promoter interaction in cancer development.

The murine IgH locus contains a distinct DNA sequence motif for the chromatin regulatory factor CTCF

Antigen receptor assembly in lymphocytes involves stringently-regulated coordination of specific DNA rearrangement events across several large chromosomal domains. Previous studies indicate that transcription factors such as paired box 5 (PAX5), Yin Yang 1 (YY1), and CCCTC-binding factor (CTCF) play a role in regulating the accessibility of the antigen receptor loci to the V(D)J recombinase, which is required for these rearrangements. To gain clues about the role of CTCF binding at the murine immunoglobulin heavy chain (IgH) locus, we utilized a computational approach that identified 144 putative CTCF-binding sites within this locus. We found that these CTCF sites share a consensus motif distinct from other CTCF sites in the mouse genome. Additionally, we could divide these CTCF sites into three categories: intergenic sites remote from any coding element, upstream sites present within 8 kb of the V_H-leader exon, and recombination signal sequence (RSS)-associated sites characteristically located at a fixed distance (∼18 bp) downstream of the RSS. We noted that the intergenic and upstream sites are located in the distal portion of the V_H locus, whereas the RSS-associated sites are located in the D_H-proximal region. Computational analysis indicated that the prevalence of CTCF-binding sites at the IgH locus is evolutionarily conserved. In all species analyzed, these sites exhibit a striking strand-orientation bias, with >98% of the murine sites being present in one orientation with respect to V_H gene transcription. Electrophoretic mobility shift and enhancer-blocking assays and ChIP–chip analysis confirmed CTCF binding to these sites both in vitro and in vivo.

Chromosome dynamics near the sol-gel phase transition dictate the timing of remote genomic interactions

Diverse antibody repertoires are generated through remote genomic interactions involving immunoglobulin variable (VH), diversity (DH) and joining (JH) gene segments. How such interactions are orchestrated remains unknown. Here we develop a strategy to track VH-DHJH motion in B-lymphocytes. We find that VH and DHJH segments are trapped in configurations that allow only local motion, such that spatially proximal segments remain in proximity, while spatially remote segments remain remote. Within a subset of cells, however, abrupt changes in VH-DHJH motion are observed, plausibly caused by temporal alterations in chromatin configurations. Comparison of experimental and simulated data suggests that constrained motion is imposed by a network of cross-linked chromatin chains characteristic of a gel phase, yet poised near the sol phase, a solution of independent chromatin chains. These results suggest that chromosome organization near the sol-gel phase transition dictates the timing of genomic interactions to orchestrate gene expression and somatic recombination.

Functional Implications of DNA Methylation in Adipose Biology

The twin epidemics of obesity and type 2 diabetes (T2D) are a serious health, social, and economic issue. The dysregulation of adipose tissue biology is central to the development of these two metabolic disorders, as adipose tissue plays a pivotal role in regulating whole-body metabolism and energy homeostasis (1). Accumulating evidence indicates that multiple aspects of adipose biology are regulated, in part, by epigenetic mechanisms. The precise and comprehensive understanding of the epigenetic control of adipose tissue biology is crucial to identifying novel therapeutic interventions that target epigenetic issues. Here, we review the recent findings on DNA methylation events and machinery in regulating the developmental processes and metabolic function of adipocytes. We highlight the following points: 1) DNA methylation is a key epigenetic regulator of adipose development and gene regulation, 2) emerging evidence suggests that DNA methylation is involved in the transgenerational passage of obesity and other metabolic disorders, 3) DNA methylation is involved in regulating the altered transcriptional landscape of dysfunctional adipose tissue, 4) genome-wide studies reveal specific DNA methylation events that associate with obesity and T2D, and 5) the enzymatic effectors of DNA methylation have physiological functions in adipose development and metabolic function.

Developmental regulation of DNA cytosine methylation at the immunoglobulin heavy chain constant locus

DNA cytosine methylation is involved in the regulation of gene expression during development and its deregulation is often associated with disease. Mammalian genomes are predominantly methylated at CpG dinucleotides. Unmethylated CpGs are often associated with active regulatory sequences while methylated CpGs are often linked to transcriptional silencing. Previous studies on CpG methylation led to the notion that transcription initiation is more sensitive to CpG methylation than transcriptional elongation. The immunoglobulin heavy chain (IgH) constant locus comprises multiple inducible constant genes and is expressed exclusively in B lymphocytes. The developmental B cell stage at which methylation patterns of the IgH constant genes are established, and the role of CpG methylation in their expression, are unknown. Here, we find that methylation patterns at most cis-acting elements of the IgH constant genes are established and maintained independently of B cell activation or promoter activity. Moreover, one of the promoters, but not the enhancers, is hypomethylated in sperm and early embryonic cells, and is targeted by different demethylation pathways, including AID, UNG, and ATM pathways. Combined, the data suggest that, rather than being prominently involved in the regulation of the IgH constant locus expression, DNA methylation may primarily contribute to its epigenetic pre-marking.

Helix-loop-helix proteins and the advent of cellular diversity: 30 years of discovery

In this review from Murre, the evolution of HLH genes, the structures of HLH domains, and the elaborate activities of HLH proteins in multicellular life are discussed.

Helix–loop–helix (HLH) proteins are dimeric transcription factors that control lineage- and developmental-specific gene programs. Genes encoding for HLH proteins arose in unicellular organisms >600 million years ago and then duplicated and diversified from ancestral genes across the metazoan and plant kingdoms to establish multicellularity. Hundreds of HLH proteins have been identified with diverse functions in a wide variety of cell types. HLH proteins orchestrate lineage specification, commitment, self-renewal, proliferation, differentiation, and homing. HLH proteins also regulate circadian clocks, protect against hypoxic stress, promote antigen receptor locus assembly, and program transdifferentiation. HLH proteins deposit or erase epigenetic marks, activate noncoding transcription, and sequester chromatin remodelers across the chromatin landscape to dictate enhancer–promoter communication and somatic recombination. Here the evolution of HLH genes, the structures of HLH domains, and the elaborate activities of HLH proteins in multicellular life are discussed.

TET2 facilitates PPARγ agonist-mediated gene regulation and insulin sensitization in adipocytes

Emerging evidence indicates that epigenetic mechanisms like DNA methylation directly contribute to metabolic regulation. For example, we previously demonstrated that de novo DNA methyltransferase Dnmt3a plays a causal role in the development of adipocyte insulin resistance. Recent studies suggest that DNA demethylation plays an important role in the developmental process of adipocytes. However, little is known about whether DNA demethylase ten-eleven translocation (TET) proteins regulate the metabolic functions of adipocytes.

METHODS

The expression of Tet genes was assessed in the fractionated adipocytes of chow- and high fat diet-fed C57/Bl6 mice using qPCR and western blotting. The effect of Tet2 gain- or loss-of-function in fully mature 3T3-L1 adipocytes in the presence/absence of Rosiglitazone (Rosi) and TNF-α on insulin sensitivity was using the insulin-stimulated glucose uptake and insulin signaling assays. Gene expression and DNA methylation analyses of PPARγ target genes was performed in the same setting. In addition, PPARγ reporter assays, co-immunoprecipitation assays, PPARγ ChIP-PCR analyses were performed.

RESULTS

We found that adipose expression of TET2, alone among its family members, was significantly reduced in diet-induced insulin resistance. TET2 gain-of-function was sufficient to promote insulin sensitivity while loss-of-function was necessary to facilitate insulin sensitization in response to the PPARγ agonist Rosiglitazone (Rosi) in cultured adipocytes. Consistent with this, TET2 was required for Rosi-dependent gene activation of certain PPARγ targets accompanied by changes in DNA demethylation at the promoter regions. Furthermore, TET2 was necessary to sustain PPARγ binding to target loci upon activation with Rosi via physical interaction with PPARγ.

CONCLUSIONS

Our data demonstrate that TET2 works as an epigenetic regulator of Rosi-mediated insulin sensitization and transcriptional regulation in adipocytes.

RAG Chromatin Scanning During V(D)J Recombination and Chromatin Loop Extrusion are Related Processes

An effective adaptive immune system depends on the ability of developing B and T cells to generate diverse immunoglobulin (Ig) and T cell receptor repertoires, respectively. Such diversity is achieved through a programmed somatic recombination process whereby germline V, D, and J segments of antigen receptor loci are assembled to form the variable region V(D)J exons of Ig and TCRs. Studies of this process, termed V(D)J recombination, have provided key insights into our understanding of a variety of general gene regulatory and DNA repair processes over the last several decades. V(D)J recombination is initiated by the RAG endonuclease which generates DNA double-stranded breaks at the borders of V, D, and J segments. In this review, we cover recent work that has elucidated RAG structure and work that revealed that RAG has a novel chromatin scanning activity, likely mediated by chromatin loop extrusion, that contributes to its ability to locate V, D, J gene segment substrates within large chromosomal loop domains bounded by CTCF-binding elements (CBEs). This latter function, coupled with the role CBE-based chromatin loop domains and subdomains within them play in focusing V(D)J recombination activity within antigen receptor loci, provide mechanistic explanations for long-standing questions regarding V(D)J segment usage diversification and in limiting potentially deleterious off-target RAG-initiated recombination events genome-wide. This review will focus mainly on studies of the mouse Ig heavy chain locus, but the principles described also apply to other Ig loci and to TCR loci in mice and humans.

CTCF-Binding Elements Mediate Accessibility of RAG Substrates During Chromatin Scanning

RAG endonuclease initiates antibody heavy chain variable region exon assembly from V, D, and J segments within a chromosomal V(D)J recombination center (RC) by cleaving between paired gene segments and flanking recombination signal sequences (RSSs). The IGCR1 control region promotes DJH intermediate formation by isolating Ds, JHs, and RCs from upstream VHs in a chromatin loop anchored by CTCF-binding elements (CBEs). How VHs access the DJHRC for VH to DJH rearrangement was unknown. We report that CBEs immediately downstream of frequently rearranged VH-RSSs increase recombination potential of their associated VH far beyond that provided by RSSs alone. This CBE activity becomes particularly striking upon IGCR1 inactivation, which allows RAG, likely via loop extrusion, to linearly scan chromatin far upstream. VH-associated CBEs stabilize interactions of D-proximal VHs first encountered by the DJHRC during linear RAG scanning and thereby promote dominant rearrangement of these VHs by an unanticipated chromatin accessibility-enhancing CBE function.

The Energetics and Physiological Impact of Cohesin Extrusion

Cohesin extrusion is thought to play a central role in establishing the architecture of mammalian genomes. However, extrusion has not been visualized in vivo, and thus, its functional impact and energetics are unknown. Using ultra-deep Hi-C, we show that loop domains form by a process that requires cohesin ATPases. Once formed, however, loops and compartments are maintained for hours without energy input. Strikingly, without ATP, we observe the emergence of hundreds of CTCF-independent loops that link regulatory DNA. We also identify architectural "stripes," where a loop anchor interacts with entire domains at high frequency. Stripes often tether super-enhancers to cognate promoters, and in B cells, they facilitate Igh transcription and recombination. Stripe anchors represent major hotspots for topoisomerase-mediated lesions, which promote chromosomal translocations and cancer. In plasmacytomas, stripes can deregulate Igh-translocated oncogenes. We propose that higher organisms have coopted cohesin extrusion to enhance transcription and recombination, with implications for tumor development.

Metabolites, genome organization, and cellular differentiation gene programs

The nutrient environment and metabolism play a dynamic role in cellular differentiation and research is elucidating the mechanisms that contribute to this process. Metabolites serve as an effective bridge that helps to translate information about nutrient states into specific interpretations of the genome. Part of this activity relates to the role for metabolites in regulating epigenetic processes as well as a newly appreciated role for metabolites in the regulation of genome organization. In this review, we will highlight recent research that has defined roles for metabolism in the organization and interpretation of the genome and how this influences cellular differentiation decisions. We will integrate information about how nutrients, such as glutamine, regulate metabolites, such as alpha-ketoglutarate, and highlight how these pathways influence epigenetic states as well as CTCF association and genome organization. We will also discuss mechanistic similarities and differences between normal differentiation states associated with embryonic stem (ES) cells and T cells and how this might relate to dysregulated states such as those associated with tumor infiltrating lymphocytes.

Variable Extent of Lineage-Specificity and Developmental Stage-Specificity of Cohesin and CCCTC-Binding Factor Binding Within the Immunoglobulin and T Cell Receptor Loci

CCCTC-binding factor (CTCF) is largely responsible for the 3D architecture of the genome, in concert with the action of cohesin, through the creation of long-range chromatin loops. Cohesin is hypothesized to be the main driver of these long-range chromatin interactions by the process of loop extrusion. Here, we performed ChIP-seq for CTCF and cohesin in two stages each of T and B cell differentiation and examined the binding pattern in all six antigen receptor (AgR) loci in these lymphocyte progenitors and in mature T and B cells, ES cells, and fibroblasts. The four large AgR loci have many bound CTCF sites, most of which are only occupied in lymphocytes, while only the CTCF sites at the end of each locus near the enhancers or J genes tend to be bound in non-lymphoid cells also. However, despite the generalized lymphocyte restriction of CTCF binding in AgR loci, the Igκ locus is the only locus that also shows significant lineage-specificity (T vs. B cells) and developmental stage-specificity (pre-B vs. pro-B) in CTCF binding. We show that cohesin binding shows greater lineage- and stage-specificity than CTCF at most AgR loci, providing more specificity to the loops. We also show that the culture of pro-B cells in IL7, a common practice to expand the number of cells before ChIP-seq, results in a CTCF-binding pattern resembling pre-B cells, as well as other epigenetic and transcriptional characteristics of pre-B cells. Analysis of the orientation of the CTCF sites show that all sites within the large V portions of the Igh and TCRβ loci have the same orientation. This suggests either a lack of requirement for convergent CTCF sites creating loops, or indicates an absence of any loops between CTCF sites within the V region portion of those loci but only loops to the convergent sites at the D-J-enhancer end of each locus. The V region portions of the Igκ and TCRα/δ loci, by contrast, have CTCF sites in both orientations, providing many options for creating CTCF-mediated convergent loops throughout the loci. CTCF/cohesin loops, along with transcription factors, drives contraction of AgR loci to facilitate the creation of a diverse repertoire of antibodies and T cell receptors.

Dynamic EBF1 occupancy directs sequential epigenetic and transcriptional events in B-cell programming

B-cell fate determination requires the action of transcription factors that operate in a regulatory network to activate B-lineage genes and repress lineage-inappropriate genes. However, the dynamics and hierarchy of events in B-cell programming remain obscure. To uncouple the dynamics of transcription factor expression from functional consequences, we generated induction systems in developmentally arrested Ebf1-/- pre-pro-B cells to allow precise experimental control of EBF1 expression in the genomic context of progenitor cells. Consistent with the described role of EBF1 as a pioneer transcription factor, we show in a time-resolved analysis that EBF1 occupancy coincides with EBF1 expression and precedes the formation of chromatin accessibility. We observed dynamic patterns of EBF1 target gene expression and sequential up-regulation of transcription factors that expand the regulatory network at the pro-B-cell stage. A continuous EBF1 function was found to be required for Cd79a promoter activity and for the maintenance of an accessible chromatin domain that is permissive for binding of other transcription factors. Notably, transient EBF1 occupancy was detected at lineage-inappropriate genes prior to their silencing in pro-B cells. Thus, persistent and transient functions of EBF1 allow for an ordered sequence of epigenetic and transcriptional events in B-cell programming.

Sexual-lineage-specific DNA methylation regulates meiosis in Arabidopsis

DNA methylation regulates eukaryotic gene expression and is extensively reprogrammed during animal development. However, whether developmental methylation reprogramming during the sporophytic life cycle of flowering plants regulates genes is presently unknown. Here we report a distinctive gene-targeted RNA-directed DNA methylation (RdDM) activity in the Arabidopsis thaliana male sexual lineage that regulates gene expression in meiocytes. Loss of sexual-lineage-specific RdDM causes mis-splicing of the MPS1 gene (also known as PRD2), thereby disrupting meiosis. Our results establish a regulatory paradigm in which de novo methylation creates a cell-lineage-specific epigenetic signature that controls gene expression and contributes to cellular function in flowering plants.

Dnmt3a is an epigenetic mediator of adipose insulin resistance

Insulin resistance results from an intricate interaction between genetic make-up and environment, and thus may be orchestrated by epigenetic mechanisms like DNA methylation. Here, we demonstrate that DNA methyltransferase 3a (Dnmt3a) is both necessary and sufficient to mediate insulin resistance in cultured mouse and human adipocytes. Furthermore, adipose-specific Dnmt3a knock-out mice are protected from diet-induced insulin resistance and glucose intolerance without accompanying changes in adiposity. Unbiased gene profiling studies revealed Fgf21 as a key negatively regulated Dnmt3a target gene in adipocytes with concordant changes in DNA methylation at the Fgf21 promoter region. Consistent with this, Fgf21 can rescue Dnmt3a-mediated insulin resistance, and DNA methylation at the FGF21 locus was elevated in human subjects with diabetes and correlated negatively with expression of FGF21 in human adipose tissue. Taken together, our data demonstrate that adipose Dnmt3a is a novel epigenetic mediator of insulin resistance in vitro and in vivo.

Non-coding Transcription Instructs Chromatin Folding and Compartmentalization to Dictate Enhancer-Promoter Communication and T Cell Fate

It is now established that Bcl11b specifies T cell fate. Here, we show that in developing T cells, the Bcl11b enhancer repositioned from the lamina to the nuclear interior. Our search for factors that relocalized the Bcl11b enhancer identified a non-coding RNA named ThymoD (thymocyte differentiation factor). ThymoD-deficient mice displayed a block at the onset of T cell development and developed lymphoid malignancies. We found that ThymoD transcription promoted demethylation at CTCF bound sites and activated cohesin-dependent looping to reposition the Bcl11b enhancer from the lamina to the nuclear interior and to juxtapose the Bcl11b enhancer and promoter into a single-loop domain. These large-scale changes in nuclear architecture were associated with the deposition of activating epigenetic marks across the loop domain, plausibly facilitating phase separation. These data indicate how, during developmental progression and tumor suppression, non-coding transcription orchestrates chromatin folding and compartmentalization to direct with high precision enhancer-promoter communication.

CCCTC-Binding Factor Translates Interleukin 2- and α-Ketoglutarate-Sensitive Metabolic Changes in T Cells into Context-Dependent Gene Programs

Despite considerable research connecting cellular metabolism with differentiation decisions, the underlying mechanisms that translate metabolite-sensitive activities into unique gene programs are still unclear. We found that aspects of the interleukin-2 (IL-2)-sensitive effector gene program in CD4⁺ and CD8⁺ T cells in type 1 conditions (Th1) were regulated by glutamine and alpha-ketoglutarate (αKG)-induced events, in part through changes in DNA and histone methylation states. We further identified a mechanism by which IL-2- and αKG-sensitive metabolic changes regulated the association of CCCTC-binding factor (CTCF) with select genomic sites. αKG-sensitive CTCF sites were often associated with loci containing IL-2- and αKG-sensitive genome organization patterns and gene expression in T cells. IL-2- and αKG-sensitive CTCF sites in T cells were also associated with genes from developmental pathways that had αKG-sensitive expression in embryonic stem cells. The data collectively support a mechanism wherein CTCF serves to translate αKG-sensitive metabolic changes into context-dependent differentiation gene programs.

Inducible CTCF insulator delays the IgH 3' regulatory region-mediated activation of germline promoters and alters class switching

Class switch recombination (CSR) plays an important role in adaptive immune response by enabling mature B cells to switch from IgM expression to the expression of downstream isotypes. CSR is preceded by inducible germline (GL) transcription of the constant genes and is controlled by the 3' regulatory region (3'RR) in a stimulus-dependent manner. Why the 3'RR-mediated up-regulation of GL transcription is delayed to the mature B-cell stage is presently unknown. Here we show that mice devoid of an inducible CTCF binding element, located in the α constant gene, display a marked isotype-specific increase of GL transcription in developing and resting splenic B cells and altered CSR in activated B cells. Moreover, insertion of a GL promoter downstream of the CTCF insulator led to premature activation of the ectopic promoter. This study provides functional evidence that the 3'RR has a developmentally controlled potential to constitutively activate GL promoters but that this activity is delayed, at least in part, by the CTCF insulator, which borders a transcriptionally active domain established by the 3'RR in developing B cells.

Biogenesis and Transcriptional Regulation of Long Noncoding RNAs in the Human Immune System

The central dogma of molecular biology states that DNA makes RNA makes protein. Discoveries over the last quarter of a century found that the process of DNA transcription into RNA gives rise to a diverse array of functional RNA species, including genes that code for protein and noncoding RNAs. For decades, the focus has been on understanding how protein-coding genes are regulated to influence protein expression. However, with the completion of the Human Genome Project and follow-up ENCODE data, it is now appreciated that only 2-3% of the genome codes for protein-coding gene exons and that the bulk of the transcribed genome, apart from ribosomal RNAs, is at the level of noncoding RNA genes. In this article, we focus on the biogenesis and regulation of a distinct class of noncoding RNA molecules termed long, noncoding RNAs in the context of the immune system.

Tet2 and Tet3 cooperate with B-lineage transcription factors to regulate DNA modification and chromatin accessibility

Ten-eleven translocation (TET) enzymes oxidize 5-methylcytosine, facilitating DNA demethylation and generating new epigenetic marks. Here we show that concomitant loss of Tet2 and Tet3 in mice at early B cell stage blocked the pro- to pre-B cell transition in the bone marrow, decreased Irf4 expression and impaired the germline transcription and rearrangement of the Igκ locus. Tet2/3-deficient pro-B cells showed increased CpG methylation at the Igκ 3' and distal enhancers that was mimicked by depletion of E2A or PU.1, as well as a global decrease in chromatin accessibility at enhancers. Importantly, re-expression of the Tet2 catalytic domain in Tet2/3-deficient B cells resulted in demethylation of the Igκ enhancers and restored their chromatin accessibility. Our data suggest that TET proteins and lineage-specific transcription factors cooperate to influence chromatin accessibility and Igκ enhancer function by modulating the modification status of DNA.

DNA Breaks and End Resection Measured Genome-wide by End Sequencing

DNA double-strand breaks (DSBs) arise during physiological transcription, DNA replication, and antigen receptor diversification. Mistargeting or misprocessing of DSBs can result in pathological structural variation and mutation. Here we describe a sensitive method (END-seq) to monitor DNA end resection and DSBs genome-wide at base-pair resolution in vivo. We utilized END-seq to determine the frequency and spectrum of restriction-enzyme-, zinc-finger-nuclease-, and RAG-induced DSBs. Beyond sequence preference, chromatin features dictate the repertoire of these genome-modifying enzymes. END-seq can detect at least one DSB per cell among 10,000 cells not harboring DSBs, and we estimate that up to one out of 60 cells contains off-target RAG cleavage. In addition to site-specific cleavage, we detect DSBs distributed over extended regions during immunoglobulin class-switch recombination. Thus, END-seq provides a snapshot of DNA ends genome-wide, which can be utilized for understanding genome-editing specificities and the influence of chromatin on DSB pathway choice.