Transcriptional stalling in B-lymphocytes: a mechanism for antibody diversification and maintenance of genomic integrity

Somatic hypermutation mechanisms during lymphomagenesis and transformation

B cells undergoing physiologically programmed or aberrant genomic alterations provide an opportune system to study the causes and consequences of genome mutagenesis. Activated B cells in germinal centers express activation-induced cytidine deaminase (AID) to accomplish physiological somatic hypermutation (SHM) of their antibody-encoding genes. In attempting to diversify their immunoglobulin (Ig) heavy- and light-chain genes, several B-cell clones successfully optimize their antigen-binding affinities. However, SHM can sometimes occur at non-Ig loci, causing genetic alternations that lay the foundation for lymphomagenesis, particularly diffuse large B-cell lymphoma. Thus, SHM acts as a double-edged sword, bestowing superb humoral immunity at the potential risk of initiating disease. We refer to off-target, non-Ig AID mutations - that are often but not always associated with disease - as aberrant SHM (aSHM). A key challenge in understanding SHM and aSHM is determining how AID targets and mutates specific DNA sequences in the Ig loci to generate antibody diversity and non-Ig genes to initiate lymphomagenesis. Herein, we discuss some current advances regarding the regulation of AID's DNA mutagenesis activity in B cells.

Noncoding mutations cause super-enhancer retargeting resulting in protein synthesis dysregulation during B cell lymphoma progression

Whole-genome sequencing of longitudinal tumor pairs representing transformation of follicular lymphoma to high-grade B cell lymphoma with MYC and BCL2 rearrangements (double-hit lymphoma) identified coding and noncoding genomic alterations acquired during lymphoma progression. Many of these transformation-associated alterations recurrently and focally occur at topologically associating domain resident regulatory DNA elements, including H3K4me3 promoter marks located within H3K27ac super-enhancer clusters in B cell non-Hodgkin lymphoma. One region found to undergo recurrent alteration upon transformation overlaps a super-enhancer affecting the expression of the PAX5/ZCCHC7 gene pair. ZCCHC7 encodes a subunit of the Trf4/5-Air1/2-Mtr4 polyadenylation-like complex and demonstrated copy number gain, chromosomal translocation and enhancer retargeting-mediated transcriptional upregulation upon lymphoma transformation. Consequently, lymphoma cells demonstrate nucleolar dysregulation via altered noncoding 5.8S ribosomal RNA processing. We find that a noncoding mutation acquired during lymphoma progression affects noncoding rRNA processing, thereby rewiring protein synthesis leading to oncogenic changes in the lymphoma proteome.

How Estrogen, Testosterone, and Sex Differences Influence Serum Immunoglobulin Isotype Patterns in Mice and Humans

Females often exhibit superior immune responses compared to males toward vaccines and pathogens such as influenza viruses and SARS-CoV-2. To help explain these differences, we first studied serum immunoglobulin isotype patterns in C57BL/6 male and female mice. We focused on IgG2b, an isotype that lends to virus control and that has been previously shown to be elevated in murine females compared to males. Improvements in IgG2b serum levels, and/or IgG2b ratios with other non-IgM isotypes, were observed when: (i) wildtype (WT) female mice were compared to estrogen receptor knockout mice (IgG2b, IgG2b/IgG3, IgG2b/IgG1, and IgG2b/IgA were all higher in WT mice), (ii) unmanipulated female mice were compared to ovariectomized mice (IgG2b/IgA was higher in unmanipulated animals), (iii) female mice were supplemented with estrogen in the context of an inflammatory insult (IgG2b and IgG2b/IgG3 were improved by estrogen supplementation), and (iv) male mice were supplemented with testosterone, a hormone that can convert to estrogen in vivo (IgG2b, IgG2b/IgG3, IgG2b/IgG1, and IgG2b/IgA were all improved by supplementation). We next examined data from three sets of previously described male and female human blood samples. In each case, there were higher IgG2 levels, and/or ratios of IgG2 with non-IgM isotypes, in human females compared to males. The effects of sex and sex hormones in the mouse and human studies were subtle, but frequent, suggesting that sex hormones represent only a fraction of the factors that influence isotype patterns. Examination of the gene loci suggested that upregulation of murine IgG2b or human IgG2 could be mediated by estrogen receptor binding to estrogen response elements and cytosine-adenine (CA) repeats upstream of respective Cγ genes. Given that murine IgG2b and human IgG2 lend to virus control, the isotype biases in females may be sufficient to improve outcomes following vaccination or infection. Future attention to sex hormone levels, and consequent immunoglobulin isotype patterns, in clinical trials are encouraged to support the optimization of vaccine and drug products for male and female hosts.

HMCES protects immunoglobulin genes specifically from deletions during somatic hypermutation

Somatic hypermutation (SHM) produces point mutations in immunoglobulin (Ig) genes in B cells when uracils created by the activation-induced deaminase are processed in a mutagenic manner by enzymes of the base excision repair (BER) and mismatch repair (MMR) pathways. Such uracil processing creates DNA strand breaks and is susceptible to the generation of deleterious deletions. Here, we demonstrate that the DNA repair factor HMCES strongly suppresses deletions without significantly affecting other parameters of SHM in mouse and human B cells, thereby facilitating the production of antigen-specific antibodies. The deletion-prone repair pathway suppressed by HMCES operates downstream from the uracil glycosylase UNG and is mediated by the combined action of BER factor APE2 and MMR factors MSH2, MSH6, and EXO1. HMCES's ability to shield against deletions during SHM requires its capacity to form covalent cross-links with abasic sites, in sharp contrast to its DNA end-joining role in class switch recombination but analogous to its genome-stabilizing role during DNA replication. Our findings lead to a novel model for the protection of Ig gene integrity during SHM in which abasic site cross-linking by HMCES intercedes at a critical juncture during processing of vulnerable gapped DNA intermediates by BER and MMR enzymes.

The role of HIRA-dependent H3.3 deposition and its modifications in the somatic hypermutation of immunoglobulin variable regions

The H3.3 histone variant and its chaperone HIRA are involved in active transcription, but their detailed roles in regulating somatic hypermutation (SHM) of immunoglobulin variable regions in human B cells are not yet fully understood. In this study, we show that the knockout (KO) of HIRA significantly decreased SHM and changed the mutation pattern of the variable region of the immunoglobulin heavy chain (IgH) in the human Ramos B cell line without changing the levels of activation-induced deaminase and other major proteins known to be involved in SHM. Except for H3K79me2/3 and Spt5, many factors related to active transcription, including H3.3, were substantively decreased in HIRA KO cells, and this was accompanied by decreased nascent transcription in the IgH locus. The abundance of ZMYND11 that specifically binds to H3.3K36me3 on the IgH locus was also reduced in the HIRA KO. Somewhat surprisingly, HIRA loss increased the chromatin accessibility of the IgH V region locus. Furthermore, stable expression of ectopic H3.3G34V and H3.3G34R mutants that inhibit both the trimethylation of H3.3K36 and the recruitment of ZMYND11 significantly reduced SHM in Ramos cells, while the H3.3K79M did not. Consistent with the HIRA KO, the H3.3G34V mutant also decreased the occupancy of various elongation factors and of ZMYND11 on the IgH variable and downstream switching regions. Our results reveal an unrecognized role of HIRA and the H3.3K36me3 modification in SHM and extend our knowledge of how transcription-associated chromatin structure and accessibility contribute to SHM in human B cells.

A Bayesian model based computational analysis of the relationship between bisulfite accessible single-stranded DNA in chromatin and somatic hypermutation of immunoglobulin genes

The B cells in our body generate protective antibodies by introducing somatic hypermutations (SHM) into the variable region of immunoglobulin genes (IgVs). The mutations are generated by activation induced deaminase (AID) that converts cytosine to uracil in single stranded DNA (ssDNA) generated during transcription. Attempts have been made to correlate SHM with ssDNA using bisulfite to chemically convert cytosines that are accessible in the intact chromatin of mutating B cells. These studies have been complicated by using different definitions of "bisulfite accessible regions" (BARs). Recently, deep-sequencing has provided much larger datasets of such regions but computational methods are needed to enable this analysis. Here we leveraged the deep-sequencing approach with unique molecular identifiers and developed a novel Hidden Markov Model based Bayesian Segmentation algorithm to characterize the ssDNA regions in the IGHV4-34 gene of the human Ramos B cell line. Combining hierarchical clustering and our new Bayesian model, we identified recurrent BARs in certain subregions of both top and bottom strands of this gene. Using this new system, the average size of BARs is about 15 bp. We also identified potential G-quadruplex DNA structures in this gene and found that the BARs co-locate with G-quadruplex structures in the opposite strand. Using various correlation analyses, there is not a direct site-to-site relationship between the bisulfite accessible ssDNA and all sites of SHM but most of the highly AID mutated sites are within 15 bp of a BAR. In summary, we developed a novel platform to study single stranded DNA in chromatin at a base pair resolution that reveals potential relationships among BARs, SHM and G-quadruplexes. This platform could be applied to genome wide studies in the future.

Role of Dot1L and H3K79 methylation in regulating somatic hypermutation of immunoglobulin genes

Somatic hypermutation (SHM) and class-switch recombination (CSR) of the immunoglobulin (Ig) genes allow B cells to make antibodies that protect us against a wide variety of pathogens. SHM is mediated by activation-induced deaminase (AID), occurs at a million times higher frequency than other mutations in the mammalian genome, and is largely restricted to the variable (V) and switch (S) regions of Ig genes. Using the Ramos human Burkitt's lymphoma cell line, we find that H3K79me2/3 and its methyltransferase Dot1L are more abundant on the V region than on the constant (C) region, which does not undergo mutation. In primary naïve mouse B cells examined ex vivo, the H3K79me2/3 modification appears constitutively in the donor Sμ and is inducible in the recipient Sγ1 upon CSR stimulation. Knockout and inhibition of Dot1L in Ramos cells significantly reduces V region mutation and the abundance of H3K79me2/3 on the V region and is associated with a decrease of polymerase II (Pol II) and its S2 phosphorylated form at the IgH locus. Knockout of Dot1L also decreases the abundance of BRD4 and CDK9 (a subunit of the P-TEFb complex) on the V region, and this is accompanied by decreased nascent transcripts throughout the IgH gene. Treatment with JQ1 (inhibitor of BRD4) or DRB (inhibitor of CDK9) decreases SHM and the abundance of Pol II S2P at the IgH locus. Since all these factors play a role in transcription elongation, our studies reinforce the idea that the chromatin context and dynamics of transcription are critical for SHM.

Topologically Associated Domains Delineate Susceptibility to Somatic Hypermutation

Somatic hypermutation (SHM) introduces point mutations into immunoglobulin (Ig) genes but also causes mutations in other parts of the genome. We have used lentiviral SHM reporter vectors to identify regions of the genome that are susceptible (“hot”) and resistant (“cold”) to SHM, revealing that SHM susceptibility and resistance are often properties of entire topologically associated domains (TADs). Comparison of hot and cold TADs reveals that while levels of transcription are equivalent, hot TADs are enriched for the cohesin loader NIPBL, super-enhancers, markers of paused/stalled RNA polymerase 2, and multiple important B cell transcription factors. We demonstrate that at least some hot TADs contain enhancers that possess SHM targeting activity and that insertion of a strong Ig SHM-targeting element into a cold TAD renders it hot. Our findings lead to a model for SHM susceptibility involving the cooperative action of cis-acting SHM targeting elements and the dynamic and architectural properties of TADs.

Senigl et al. show that genome susceptibility to somatic hypermutation (SHM) is confined within topologically associated domains (TADs) and is linked to markers of strong enhancers and stalled transcription and high levels of the cohesin loader NIPBL. Insertion of an ectopic SHM targeting element renders an entire TAD susceptible to SHM.

Complex sex-biased antibody responses: estrogen receptors bind estrogen response elements centered within immunoglobulin heavy chain gene enhancers

Nuclear hormone receptors including the estrogen receptor (ERα) and the retinoic acid receptor regulate a plethora of biological functions including reproduction, circulation and immunity. To understand how estrogen and other nuclear hormones influence antibody production, we characterized total serum antibody isotypes in female and male mice of C57BL/6J, BALB/cJ and C3H/HeJ mouse strains. Antibody levels were higher in females compared to males in all strains and there was a female preference for IgG2b production. Sex-biased patterns were influenced by vitamin levels, and by antigen specificity toward influenza virus or pneumococcus antigens. To help explain sex biases, we examined the direct effects of estrogen on immunoglobulin heavy chain sterile transcript production among purified, lipopolysaccharide-stimulated B cells. Supplemental estrogen in B-cell cultures significantly increased immunoglobulin heavy chain sterile transcripts. Chromatin immunoprecipitation analyses of activated B cells identified significant ERα binding to estrogen response elements (EREs) centered within enhancer elements of the immunoglobulin heavy chain locus, including the Eµ enhancer and hypersensitive site 1,2 (HS1,2) in the 3' regulatory region. The ERE in HS1,2 was conserved across animal species, and in humans marked a site of polymorphism associated with the estrogen-augmented autoimmune disease, lupus. Taken together, the results highlight: (i) the important targets of ERα in regulatory regions of the immunoglobulin heavy chain locus that influence antibody production, and (ii) the complexity of mechanisms by which estrogen instructs sex-biased antibody production profiles.

Transcription factor binding at Ig enhancers is linked to somatic hypermutation targeting

Secondary diversification of the Ig repertoire occurs through somatic hypermutation (SHM), gene conversion (GCV), and class switch recombination (CSR)-three processes that are initiated by activation-induced cytidine deaminase (AID). AID targets Ig genes at orders of magnitude higher than the rest of the genome, but the basis for this specificity is poorly understood. We have previously demonstrated that enhancers and enhancer-like sequences from Ig genes are capable of stimulating SHM of neighboring genes in a capacity distinct from their roles in increasing transcription. Here, we use an in vitro proteomics approach to identify E-box, MEF2, Ets, and Ikaros transcription factor family members as potential binders of these enhancers. ChIP assays in the hypermutating Ramos B cell line confirmed that many of these factors bound the endogenous Igλ enhancer and/or the IgH intronic enhancer (Eμ) in vivo. Further investigation using SHM reporter assays identified binding sites for E2A and MEF2B in Eμ and demonstrated an association between loss of factor binding and decreases in the SHM stimulating activity of Eμ mutants. Our results provide novel insights into trans-acting factors that dictate SHM targeting and link their activity to specific DNA binding sites within Ig enhancers.

Regulation of chromatin remodeling through RNA polymerase II stalling in the immune system

RNA polymerase II (Pol II) binds to promoter-proximal regions of inducible target genes that are controlled and not transcribed by several negative elongation factors, which is known as Pol II stalling. The occurrence of stalling is due to particular modification signatures and structural conformations of chromatin that affect Pol II elongation. The existence and physiological importance of Pol II stalling implies that there is a dynamic balance in chromatin regulation prior to endogenous or exogenous stimulation. In this review, we discuss the effects of ATP-dependent chromatin remodeling complexes and histone modification via transcriptional machinery Pol II C-terminal domain phosphorylated at serine 5 (S5P RNAPII) initiation and S2P RNAPII elongation on the expression or silence of specific genes after the production of activated or differentiated signals or cytokines. The response occurs immediately during immune cell development and function, and it also includes the generation of immunological memories. This summary suggests that the host immune response genes involve a novel mechanism of selectively regulatory chromatin remodeling, a fundamental and crucial aspect of epigenetic regulation.

RNA exosome-regulated long non-coding RNA transcription controls super-enhancer activity

We have ablated the cellular RNA degradation machinery in differentiated B cells and pluripotent embryonic stem cells (ESCs) by conditional mutagenesis of core (Exosc3) and nuclear RNase (Exosc10) components of RNA exosome and identified a vast number of long non-coding RNAs (lncRNAs) and enhancer RNAs (eRNAs) with emergent functionality. Unexpectedly, eRNA-expressing regions accumulate R-loop structures upon RNA exosome ablation, thus demonstrating the role of RNA exosome in resolving deleterious DNA/RNA hybrids arising from active enhancers. We have uncovered a distal divergent eRNA-expressing element (lncRNA-CSR) engaged in long-range DNA interactions and regulating IgH 3' regulatory region super-enhancer function. CRISPR-Cas9-mediated ablation of lncRNA-CSR transcription decreases its chromosomal looping-mediated association with the IgH 3' regulatory region super-enhancer and leads to decreased class switch recombination efficiency. We propose that the RNA exosome protects divergently transcribed lncRNA expressing enhancers by resolving deleterious transcription-coupled secondary DNA structures, while also regulating long-range super-enhancer chromosomal interactions important for cellular function.

RNA Exosome Regulates AID DNA Mutator Activity in the B Cell Genome

The immunoglobulin diversification processes of somatic hypermutation and class switch recombination critically rely on transcription-coupled targeting of activation-induced cytidine deaminase (AID) to Ig loci in activated B lymphocytes. AID catalyzes deamination of cytidine deoxynucleotides on exposed single-stranded DNA. In addition to driving immunoglobulin diversity, promiscuous targeting of AID mutagenic activity poses a deleterious threat to genomic stability. Recent genome-wide studies have uncovered pervasive AID activity throughout the B cell genome. It is increasingly apparent that AID activity is frequently targeted to genomic loci undergoing early transcription termination where RNA exosome promotes the resolution of stalled transcription complexes via cotranscriptional RNA degradation mechanisms. Here, we review aspects and consequences of eukaryotic transcription that lead to early termination, RNA exosome recruitment, and ultimately targeting of AID mutagenic activity.

Overlapping hotspots in CDRs are critical sites for V region diversification

Activation-induced deaminase (AID) mediates the somatic hypermutation (SHM) of Ig variable (V) regions that is required for the affinity maturation of the antibody response. An intensive analysis of a published database of somatic hypermutations that arose in the IGHV3-23*01 human V region expressed in vivo by human memory B cells revealed that the focus of mutations in complementary determining region (CDR)1 and CDR2 coincided with a combination of overlapping AGCT hotspots, the absence of AID cold spots, and an abundance of polymerase eta hotspots. If the overlapping hotspots in the CDR1 or CDR2 did not undergo mutation, the frequency of mutations throughout the V region was reduced. To model this result, we examined the mutation of the human IGHV3-23*01 biochemically and in the endogenous heavy chain locus of Ramos B cells. Deep sequencing revealed that IGHV3-23*01 in Ramos cells accumulates AID-induced mutations primarily in the AGCT in CDR2, which was also the most frequent site of mutation in vivo. Replacing the overlapping hotspots in CDR1 and CDR2 with neutral or cold motifs resulted in a reduction in mutations within the modified motifs and, to some degree, throughout the V region. In addition, some of the overlapping hotspots in the CDRs were at sites in which replacement mutations could change the structure of the CDR loops. Our analysis suggests that the local sequence environment of the V region, and especially of the CDR1 and CDR2, is highly evolved to recruit mutations to key residues in the CDRs of the IgV region.

Ubiquitination events that regulate recombination of immunoglobulin Loci gene segments

Programed DNA mutagenesis events in the immunoglobulin (Ig) loci of developing B cells utilize the common and conserved mechanism of protein ubiquitination for subsequent proteasomal degradation to generate the required antigen-receptor diversity. Recombinase proteins RAG1 and RAG2, necessary for V(D)J recombination, and activation-induced cytidine deaminase, an essential mutator protein for catalyzing class switch recombination and somatic hypermutation, are regulated by various ubiquitination events that affect protein stability and activity. Programed DNA breaks in the Ig loci can be identified by various components of DNA repair pathways, also regulated by protein ubiquitination. Errors in the ubiquitination pathways for any of the DNA double-strand break repair proteins can lead to inefficient recombination and repair events, resulting in a compromised adaptive immune system or development of cancer.

E3-ubiquitin ligase Nedd4 determines the fate of AID-associated RNA polymerase II in B cells

Programmed mutagenesis of the immunoglobulin locus of B lymphocytes during class switch recombination (CSR) and somatic hypermutation requires RNA polymerase II (polII) transcription complex-dependent targeting of the DNA mutator activation-induced cytidine deaminase (AID). AID deaminates cytidine residues on substrate sequences in the immunoglobulin (Ig) locus via a transcription-dependent mechanism, and this activity is stimulated by the RNA polII stalling cofactor Spt5 and the 11-subunit cellular noncoding RNA 3'-5' exonucleolytic processing complex RNA exosome. The mechanism by which the RNA exosome recognizes immunoglobulin locus RNA substrates to stimulate AID DNA deamination activity on its in vivo substrate sequences is an important question. Here we report that E3-ubiquitin ligase Nedd4 destabilizes AID-associated RNA polII by a ubiquitination event, leading to generation of 3' end free RNA exosome RNA substrates at the Ig locus and other AID target sequences genome-wide. We found that lack of Nedd4 activity in B cells leads to accumulation of RNA exosome substrates at AID target genes and defective CSR. Taken together, our study links noncoding RNA processing following RNA polII pausing with regulation of the mutator AID protein. Our study also identifies Nedd4 as a regulator of noncoding RNAs that are generated by stalled RNA polII genome-wide.