Nuclear positioning rather than contraction controls ordered rearrangements of immunoglobulin loci

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.

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.

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.

Contribution of Immunoglobulin Enhancers to B Cell Nuclear Organization

B cells undergo genetic rearrangements at immunoglobulin gene (Ig) loci during B cell maturation. First V(D)J recombination occurs during early B cell stages followed by class switch recombination (CSR) and somatic hypermutation (SHM) which occur during mature B cell stages. Given that RAG1/2 induces DNA double strand breaks (DSBs) during V(D)J recombination and AID (Activation-Induced Deaminase) leads to DNA modifications (mutations during SHM or DNA DSBs during CSR), it is mandatory that IgH rearrangements be tightly regulated to avoid any mutations or translocations within oncogenes. Ig loci contain various cis-regulatory elements that are involved in germline transcription, chromatin modifications or RAG/AID recruitment. Ig cis-regulatory elements are increasingly recognized as being involved in nuclear positioning, heterochromatin addressing and chromosome loop regulation. In this review, we examined multiple data showing the critical interest of studying Ig gene regulation at the whole nucleus scale. In this context, we highlighted the essential function of Ig gene regulatory elements that now have to be considered as nuclear organizers in B lymphocytes.

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.

Genomic Repeats Categorize Genes with Distinct Functions for Orchestrated Regulation

Repetitive elements are abundantly distributed in mammalian genomes. Here, we reveal a striking association between repeat subtypes and gene function. SINE, L1, and low-complexity repeats demarcate distinct functional categories of genes and may dictate the time and level of gene expression by providing binding sites for different regulatory proteins. Importantly, imaging and sequencing analysis show that L1 repeats sequester a large set of genes with specialized functions in nucleolus- and lamina-associated inactive domains that are depleted of SINE repeats. In addition, L1 transcripts bind extensively to its DNA in embryonic stem cells (ESCs). Depletion of L1 RNA in ESCs leads to relocation of L1-enriched chromosomal segments from inactive domains to the nuclear interior and de-repression of L1-associated genes. These results demonstrate a role of L1 DNA and RNA in gene silencing and suggest a general theme of genomic repeats in orchestrating the function, regulation, and expression of their host genes.

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.

A B-Cell-Specific Enhancer Orchestrates Nuclear Architecture to Generate a Diverse Antigen Receptor Repertoire

The genome is organized into topologically associated domains (TADs) that enclose smaller subTADs. Here, we identify and characterize an enhancer that is located in the middle of the V gene region of the immunoglobulin kappa light chain (Igκ) locus that becomes active preceding the stage at which this locus undergoes V(D)J recombination. This enhancer is a hub of long-range chromatin interactions connecting subTADs in the V gene region with the recombination center at the J genes. Deletion of this element results in a highly altered long-range chromatin interaction pattern across the locus and, importantly, affects individual V gene utilization locus-wide. These results indicate the existence of an enhancer-dependent framework in the Igκ locus and further suggest that the composition of the diverse antibody repertoire is regulated in a subTAD-specific manner. This enhancer thus plays a structural role in orchestrating the proper folding of the Igκ locus in preparation for V(D)J recombination.

Epigenetic Enhancer Marks and Transcription Factor Binding Influence Vκ Gene Rearrangement in Pre-B Cells and Pro-B Cells

To date there has not been a study directly comparing relative Igκ rearrangement frequencies obtained from genomic DNA (gDNA) and cDNA and since each approach has potential biases, this is an important issue to clarify. Here we used deep sequencing to compare the unbiased gDNA and RNA Igκ repertoire from the same pre-B cell pool. We find that ~20% of Vκ genes have rearrangement frequencies ≥2-fold up or down in RNA vs. DNA libraries, including many members of the Vκ3, Vκ4, and Vκ6 families. Regression analysis indicates Ikaros and E2A binding are associated with strong promoters. Within the pre-B cell repertoire, we observed that individual Vκ genes rearranged at very different frequencies, and also displayed very different Jκ usage. Regression analysis revealed that the greatly unequal Vκ gene rearrangement frequencies are best predicted by epigenetic marks of enhancers. In particular, the levels of newly arising H3K4me1 peaks associated with many Vκ genes in pre-B cells are most predictive of rearrangement levels. Since H3K4me1 is associated with long range chromatin interactions which are created during locus contraction, our data provides mechanistic insight into unequal rearrangement levels. Comparison of Igκ rearrangements occurring in pro-B cells and pre-B cells from the same mice reveal a pro-B cell bias toward usage of Jκ-distal Vκ genes, particularly Vκ10-96 and Vκ1-135. Regression analysis indicates that PU.1 binding is the highest predictor of Vκ gene rearrangement frequency in pro-B cells. Lastly, the repertoires of iEκ^−/− pre-B cells reveal that iEκ actively influences Vκ gene usage, particularly Vκ3 family genes, overlapping with a zone of iEκ-regulated germline transcription. These represent new roles for iEκ in addition to its critical function in promoting overall Igκ rearrangement. Together, this study provides insight into many aspects of Igκ repertoire formation.

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.