A Hyper-IgM Syndrome Mutation in Activation-Induced Cytidine Deaminase Disrupts G-Quadruplex Binding and Genome-wide Chromatin Localization

B-Cell Epigenetic Modulation of IgA Response by 5-Azacytidine and IgA Nephropathy

Key Points

Dysregulated IgA production plays a key role in the pathogenesis of IgA nephropathy. Increased 5-methylcytosine modification, an epigenetic regulatory mechanism, exaggerated IgA nephropathy phenotype in mice. Conversely, inhibition of 5-methylcytosine modification ameliorated progression of IgA nephropathy–like kidney disease in mice.

Background

IgA nephropathy is an important global cause of kidney failure. Dysregulation of IgA production is believed to play a key role in IgA nephropathy pathogenesis; however, little is known about the epigenetic mechanisms, such as RNA 5-methylcytosine (5mC) modification, in regulating IgA synthesis.

Methods

To decipher the role of RNA 5mC in regulation of IgA class switch, the microRNA (miR)-23b−/− and Lactobacillus casei (Chinese Industrial Microbial Culture Collection Center) cell wall extract–induced Kawasaki disease mice were treated with 5-azacytidine. Trdmt1 −/− and double Trdmt1 −/−/miR-23b −/− mice and Aid −/− mice or Aid −/−/miR-23b −/− mice were also used.

Results

We showed that miR-23b downregulated expression of Transfer RNA Aspartic Acid Methyltransferase 1 and consequently reduced 5mC (m5C) RNA modification and IgA synthesis in B cells. Inhibition of m5C RNA modification normalized serum IgA levels and ameliorated progression of the IgA nephropathy–like kidney disease in miR-23b −/− and Kawasaki disease mice, while mesangial IgA and C3 deposition failed to develop in Trdmt1 −/−miR-23b −/− mice. By contrast, increased m5C RNA modification resulted in an exaggerated IgA nephropathy phenotype. miR-23b regulation of serum IgA levels and the development of an IgA nephropathy–like kidney disease in miR-23b −/− and Kawasaki disease mice is likely mediated through TRDMT1-driven 5mC RNA modification in B cells, resulting in impaired activation-induced cytidine deaminase activity and IgA class switch recombination.

Conclusions

This study revealed TRDMT1-induced RNA 5mC methylation regulated IgA class switch, and inhibition of RNA 5mC by 5-azacytidine ameliorated progression of IgA nephropathy.

Taming AID mutator activity in somatic hypermutation

Activation-induced cytidine deaminase (AID) initiates somatic hypermutation (SHM) by introducing base substitutions into antibody genes, a process enabling antibody affinity maturation in immune response. How a mutator is tamed to precisely and safely generate programmed DNA lesions in a physiological process remains unsettled, as its dysregulation drives lymphomagenesis. Recent research has revealed several hidden features of AID-initiated mutagenesis: preferential activity on flexible DNA substrates, restrained activity within chromatin loop domains, unique DNA repair factors to differentially decode AID-caused lesions, and diverse consequences of aberrant deamination. Here, we depict the multifaceted regulation of AID activity with a focus on emerging concepts/factors and discuss their implications for the design of base editors (BEs) that install somatic mutations to correct deleterious genomic variants.

Heterologous booster vaccination enhances antibody responses to SARS-CoV-2 by improving Tfh function and increasing B-cell clonotype SHM frequency

Heterologous prime-boost has broken the protective immune response bottleneck of the COVID-19 vaccines. however, the underlying mechanisms have not been fully elucidated. Here, we investigated antibody responses and explored the response of germinal center (GC) to priming with inactivated vaccines and boosting with heterologous adenoviral-vectored vaccines or homologous inactivated vaccines in mice. Antibody responses were dramatically enhanced by both boosting regimens. Heterologous immunization induced more robust GC activation, characterized by increased Tfh cell populations and enhanced helper function. Additionally, increased B-cell activation and antibody production were observed in a heterologous regimen. Libra-seq was used to compare the differences of S1-, S2- and NTD-specific B cells between homologous and heterologous vaccination, respectively. S2-specific CD19+ B cells presented increased somatic hypermutations (SHMs), which were mainly enriched in plasma cells. Moreover, a heterologous booster dose promoted the clonal expansion of B cells specific to S2 and NTD regions. In conclusion, the functional role of Tfh and B cells following SARS-CoV-2 heterologous vaccination may be important for modulating antibody responses. These findings provide new insights for the development of SARS-CoV-2 vaccines that induce more robust antibody response.

The immunological understanding on germinal center B cells in psoriasis

The development of psoriasis is mainly driven by the dysregulation of T cells within the skin, marking a primary involvement of these cells in the pathogenesis. Although B cells are integral components of the immune system, their role in the initiation and progression of psoriasis is not as pivotal as that of T cells. The paradox of B cell suggests that, while it is crucial for adaptive immunity, B cells may contribute to the exacerbation of psoriasis. Numerous ideas proposed that there are potential relationships between psoriasis and B cells especially within germinal centers (GCs). Recent research projected that B cells might be triggered by autoantigens which then induced molecular mimicry to alter B cells activity within GC and generate autoantibodies and pro-inflammatory cytokines, form ectopic GC, and dysregulate the proliferation of keratinocytes. Hence, in this review, we gathered potential evidence indicating the participation of B cells in psoriasis within the context of GC, aiming to enhance our comprehension and advance treatment strategies for psoriasis thus inviting many new researchers to investigate this issue.

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.

Alternative DNA structures in hematopoiesis and adaptive immunity

Besides the canonical B-form, DNA also adopts alternative non-B form conformations which are highly conserved in all domains of life. While extensive research over decades has centered on the genomic functions of B-form DNA, understanding how non-B-form conformations influence functional genomic states remains a fundamental and open question. Recent studies have ascribed alternative DNA conformations such as G-quadruplexes and R-loops as important functional features in eukaryotic genomes. This review delves into the biological importance of alternative DNA structures, with a specific focus on hematopoiesis and adaptive immunity. We discuss the emerging roles of G-quadruplex and R-loop structures, the two most well-studied alternative DNA conformations, in the hematopoietic compartment and present evidence for their functional roles in normal cellular physiology and associated pathologies.

Unique repetitive nucleic acid structures mirror switch regions in the human IgH locus

Immunoglobulin (Ig) genes carry the unique ability to be reshaped in peripheral B lymphocytes after these cells encounter a specific antigen. B cells can then further improve their affinity, acquire new functions as memory cells and eventually end up as antibody-secreting cells. Ig class switching is an important change that occurs in this context, thanks to local DNA lesions initiated by the enzyme activation-induced deaminase (AID). Several cis-acting elements of the Ig heavy (IgH) chain locus make it accessible to the AID-mediated lesions that promote class switch recombination (CSR). DNA repeats, with a non-template strand rich in G-quadruplexes (G4)-DNA, are prominent cis-targets of AID and define the so-called "switch" (S) regions specifically targeted for CSR. By analyzing the structure of the human IgH locus, we uncover that abundant DNA repeats, some with a putative G4-rich template strand, are additionally present in downstream portions of the IgH coding genes. These like-S (LS) regions stand as 3' mirror-images of S regions and also show analogies to some previously reported repeats associated with the IgH locus 3' super-enhancer. A regulatory role of LS repeats is strongly suggested by their specific localization close to exons encoding the membrane form of Ig molecules, and by their conservation during mammalian evolution.

Detection of alternative DNA structures and its implications for human disease

Around 3% of the genome consists of simple DNA repeats that are prone to forming alternative (non-B) DNA structures, such as hairpins, cruciforms, triplexes (H-DNA), four-stranded guanine quadruplexes (G4-DNA), and others, as well as composite RNA:DNA structures (e.g., R-loops, G-loops, and H-loops). These DNA structures are dynamic and favored by the unwinding of duplex DNA. For many years, the association of alternative DNA structures with genome function was limited by the lack of methods to detect them in vivo. Here, we review the recent advancements in the field and present state-of-the-art technologies and methods to study alternative DNA structures. We discuss the limitations of these methods as well as how they are beginning to provide insights into causal relationships between alternative DNA structures, genome function and stability, and human disease.

Distinct metabolic requirements regulate B cell activation and germinal center responses

Following infection or vaccination, activated B cells at extrafollicular sites or within germinal centers (GCs) undergo vigorous clonal proliferation. Proliferating lymphocytes have been shown to undertake lactate dehydrogenase A (LDHA)-dependent aerobic glycolysis; however, the specific role of this metabolic pathway in a B cell transitioning from a naïve to a highly proliferative, activated state remains poorly defined. Here, we deleted LDHA in a stage-specific and cell-specific manner. We find that ablation of LDHA in a naïve B cell did not profoundly affect its ability to undergo a bacterial lipopolysaccharide-induced extrafollicular B cell response. On the other hand, LDHA-deleted naïve B cells had a severe defect in their capacities to form GCs and mount GC-dependent antibody responses. In addition, loss of LDHA in T cells severely compromised B cell-dependent immune responses. Strikingly, when LDHA was deleted in activated, as opposed to naïve, B cells, there were only minimal effects on the GC reaction and in the generation of high-affinity antibodies. These findings strongly suggest that naïve and activated B cells have distinct metabolic requirements that are further regulated by niche and cellular interactions.

Conservation of the insert-2 motif confers Rev1 from different species with an ability to disrupt G-quadruplexes and stimulate translesion DNA synthesis

In some organisms, the replication of G-quadruplex (G4) structures is supported by the Rev1 DNA polymerase. We previously showed that residues in the insert-2 motif of human Rev1 (hRev1) increased the affinity of the enzyme for G4 DNA and mediated suppression of mutagenic replication near G4 motifs. We have now investigated the conservation of G4-selective properties in Rev1 from other species. We compared Rev1 from Danio rerio (zRev1), Saccharomyces cerevisiae (yRev1), and Leishmania donovani (lRev1) with hRev1, including an insert-2 mutant form of hRev1 (E466A/Y470A or EY). We found that zRev1 retained all of the G4-selective prowess of the human enzyme, but there was a marked attenuation of G4 binding affinity for the EY hRev1 mutant and the two Rev1 proteins lacking insert-2 (yRev1 and lRev1). Perhaps most strikingly, we found that insert-2 was important for disruption of the G4 structure and optimal stimulation of processive DNA synthesis across the guanine-rich motif by DNA polymerase kappa (pol κ). Our findings have implications for how Rev1 might contribute to G4 replication in different species spanning the evolutionary tree - signaling the importance of selection for enzymes with robust G4-selective properties in organisms where these non-B DNA structures may fulfill taxa-specific physiological functions.

Roles of G4-DNA and G4-RNA in Class Switch Recombination and Additional Regulations in B-Lymphocytes

Mature B cells notably diversify immunoglobulin (Ig) production through class switch recombination (CSR), allowing the junction of distant "switch" (S) regions. CSR is initiated by activation-induced deaminase (AID), which targets cytosines adequately exposed within single-stranded DNA of transcribed targeted S regions, with a specific affinity for WRCY motifs. In mammals, G-rich sequences are additionally present in S regions, forming canonical G-quadruplexes (G4s) DNA structures, which favor CSR. Small molecules interacting with G4-DNA (G4 ligands), proved able to regulate CSR in B lymphocytes, either positively (such as for nucleoside diphosphate kinase isoforms) or negatively (such as for RHPS4). G4-DNA is also implicated in the control of transcription, and due to their impact on both CSR and transcriptional regulation, G4-rich sequences likely play a role in the natural history of B cell malignancies. Since G4-DNA stands at multiple locations in the genome, notably within oncogene promoters, it remains to be clarified how it can more specifically promote legitimate CSR in physiology, rather than pathogenic translocation. The specific regulatory role of G4 structures in transcribed DNA and/or in corresponding transcripts and recombination hereby appears as a major issue for understanding immune responses and lymphomagenesis.

Heterologous expression of hAID in E. coli leads to the production of a splice isoform of AID: hAIDδC, a mystery to be explored

Activation-induced cytidine deaminase (AID) is a key player that initiates antibody diversification in activated B-cell. AID mediates somatic hypermutation (SHM) and class switch recombination (CSR) via the deamination of cytosine to uracil at the Ig locus, resulting in the production of high-affinity antibodies. AID is predominantly restricted to Ig genes, whereas off-targeting of AID leads to lymphocyte-related malignancies. Interestingly, apart from FL-AID other splice isoforms of AID are highly expressed in the lymphocyte malignancies. In our study, we found that the heterologous expression of hAID-FL in E. coli cells produced two induced bands of hAID as demonstrated by SDS-PAGE and western blotting. Remarkably, peptide mapping data predicted that one band is hAID-FL and the other is its splice isoform, hAIDδE4a. To get an insight into why E. coli cells expressed hAID-FL and hAID variant, we mutated the 5' and 3' splice site of a putative intron of hAID, but it failed to produce only hAID-FL. Incidentally, hAID expressed with fusion partners also displayed two bands, and peptide mapping data strongly suggest that besides hAID-FL, the lower band showed a significant number of amino acids missing towards the C-terminal domain (named as hAIDδC). Our results are the first report to show that expression of recombinant hAID alone or irrespective of solubilization tags in E. coli cells produced hAID-FL and hAIDδC. It will be fascinating to explore the potential mechanism underlying the expression of hAIDδC from recombinant hAID plasmid in E. coli cells.

AFF3, a susceptibility factor for autoimmune diseases, is a molecular facilitator of immunoglobulin class switch recombination

Immunoglobulin class switch recombination (CSR) plays critical roles in controlling infections and inflammatory tissue injuries. Here, we show that AFF3, a candidate gene for both rheumatoid arthritis and type 1 diabetes, is a molecular facilitator of CSR with an isotype preference. Aff3-deficient mice exhibit low serum levels of immunoglobulins, predominantly immunoglobulin G2c (IgG2c) followed by IgG1 and IgG3 but not IgM. Furthermore, Aff3-deficient mice show weak resistance to Plasmodium yoelii infection, confirming that Aff3 modulates immunity to this pathogen. Mechanistically, the AFF3 protein binds to the IgM and IgG1 switch regions via a C-terminal domain, and Aff3 deficiency reduces the binding of AID to the switch regions less efficiently. One AFF3 risk allele for rheumatoid arthritis is associated with high mRNA expression of AFF3, IGHG2, and IGHA2 in human B cells. These findings demonstrate that AFF3 directly regulates CSR by facilitating the recruitment of AID to the switch regions.

Triple-helix potential of the mouse genome

Certain DNA sequences, including mirror-symmetric polypyrimidine•polypurine runs, are capable of folding into a triple-helix–containing non–B-form DNA structure called H-DNA. Such H-DNA–forming sequences occur frequently in many eukaryotic genomes, including in mammals, and multiple lines of evidence indicate that these motifs are mutagenic and can impinge on DNA replication, transcription, and other aspects of genome function. In this study, we show that the triplex-forming potential of H-DNA motifs in the mouse genome can be evaluated using S1-sequencing (S1-seq), which uses the single-stranded DNA (ssDNA)–specific nuclease S1 to generate deep-sequencing libraries that report on the position of ssDNA throughout the genome. When S1-seq was applied to genomic DNA isolated from mouse testis cells and splenic B cells, we observed prominent clusters of S1-seq reads that appeared to be independent of endogenous double-strand breaks, that coincided with H-DNA motifs, and that correlated strongly with the triplex-forming potential of the motifs. Fine-scale patterns of S1-seq reads, including a pronounced strand asymmetry in favor of centrally positioned reads on the pyrimidine-containing strand, suggested that this S1-seq signal is specific for one of the four possible isomers of H-DNA (H-y5). By leveraging the abundance and complexity of naturally occurring H-DNA motifs across the mouse genome, we further defined how polypyrimidine repeat length and the presence of repeat-interrupting substitutions modify the structure of H-DNA. This study provides an approach for studying DNA secondary structure genome-wide at high spatial resolution.

The role of B cells in the development, progression, and treatment of lymphomas and solid tumors

B cells are integral components of the mammalian immune response as they have the ability to generate antibodies against an almost infinite array of antigens. Over the past several decades, significant scientific progress has been made in understanding that this enormous B cell diversity contributes to pathogen clearance. However, our understanding of the humoral response to solid tumors and to tumor-specific antigens is unclear. In this review, we first discuss how B cells interact with other cells in the tumor microenvironment and influence the development and progression of various solid tumors. The ability of B lymphocytes to generate antibodies against a diverse repertoire of antigens and subsequently tailor the humoral immune response to specific pathogens relies on their ability to undergo genomic alterations during their development and differentiation. We will discuss key transforming events that lead to the development of B cell lymphomas. Overall, this review provides a foundation for innovative therapeutic interventions for both lymphoma and solid tumor malignancies.

TET deficiency perturbs mature B cell homeostasis and promotes oncogenesis associated with accumulation of G-quadruplex and R-loop structures

Enzymes of the TET family are methylcytosine dioxygenases that undergo frequent mutational or functional inactivation in human cancers. Recurrent loss-of-function mutations in TET proteins are frequent in human diffuse large B cell lymphoma (DLBCL). Here, we investigate the role of TET proteins in B cell homeostasis and development of B cell lymphomas with features of DLBCL. We show that deletion of Tet2 and Tet3 genes in mature B cells in mice perturbs B cell homeostasis and results in spontaneous development of germinal center (GC)-derived B cell lymphomas with increased G-quadruplexes and R-loops. At a genome-wide level, G-quadruplexes and R-loops were associated with increased DNA double-strand breaks (DSBs) at immunoglobulin switch regions. Deletion of the DNA methyltransferase DNMT1 in TET-deficient B cells prevented expansion of GC B cells, diminished the accumulation of G-quadruplexes and R-loops and delayed B lymphoma development, consistent with the opposing functions of DNMT and TET enzymes in DNA methylation and demethylation. Clustered regularly interspaced short palindromic repeats (CRISPR)-mediated depletion of nucleases and helicases that regulate G-quadruplexes and R-loops decreased the viability of TET-deficient B cells. Our studies suggest a molecular mechanism by which TET loss of function might predispose to the development of B cell malignancies.

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.

The Relevance of G-Quadruplexes for DNA Repair

DNA molecules can adopt a variety of alternative structures. Among these structures are G-quadruplex DNA structures (G4s), which support cellular function by affecting transcription, translation, and telomere maintenance. These structures can also induce genome instability by stalling replication, increasing DNA damage, and recombination events. G-quadruplex-driven genome instability is connected to tumorigenesis and other genetic disorders. In recent years, the connection between genome stability, DNA repair and G4 formation was further underlined by the identification of multiple DNA repair proteins and ligands which bind and stabilize said G4 structures to block specific DNA repair pathways. The relevance of G4s for different DNA repair pathways is complex and depends on the repair pathway itself. G4 structures can induce DNA damage and block efficient DNA repair, but they can also support the activity and function of certain repair pathways. In this review, we highlight the roles and consequences of G4 DNA structures for DNA repair initiation, processing, and the efficiency of various DNA repair pathways.

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.

Characterization of DNA G-Quadruplex Structures in Human Immunoglobulin Heavy Variable (IGHV) Genes

Activation-induced deaminase (AID) is a key enzyme involved in antibody diversification by initiating somatic hypermutation (SHM) and class-switch recombination (CSR) of the Immunoglobulin (Ig) loci. AID preferentially targets WRC (W=A/T, R=A/G) hotspot motifs and avoids SYC (S=C/G, Y=C/T) coldspots. G-quadruplex (G4) structures are four-stranded DNA secondary structures with key functions in transcription, translation and replication. In vitro studies have shown G4s to form and bind AID in Ig switch (S) regions. Alterations in the gene encoding AID can further disrupt AID-G4 binding and reduce CSR in vivo. However, it is still unclear whether G4s form in the variable (V) region, or how they may affect SHM. To assess the possibility of G4 formation in human V regions, we analyzed germline human Ig heavy chain V (IGHV) sequences, using a pre-trained deep learning model that predicts G4 potential. This revealed that many genes from the IGHV3 and IGHV4 families are predicted to have high G4 potential in the top and bottom strand, respectively. Different IGHV alleles also showed variability in G4 potential. Using a high-resolution (G4-seq) dataset of biochemically confirmed potential G4s in IGHV genes, we validated our computational predictions. G4-seq also revealed variation between S and V regions in the distribution of potential G4s, with the V region having overall reduced G4 abundance compared to the S region. The density of AGCT motifs, where two AGC hotspots overlap on both strands, was roughly 2.6-fold greater in the V region than the Constant (C) region, which does not mutate despite having predicted G4s at similar levels. However, AGCT motifs in both V and C regions were less abundant than in S regions. In silico mutagenesis experiments showed that G4 potentials were generally robust to mutation, although large deviations from germline states were found, mostly in framework regions. G4 potential is also associated with higher mutability of certain WRC hotspots on the same strand. In addition, CCC coldspots opposite a predicted G4 were shown to be targeted significantly more for mutation. Our overall assessment reveals plausible evidence of functional G4s forming in the Ig V region.