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.

Integration of Mutational Signature Analysis with 3D Chromatin Data Unveils Differential AID-Related Mutagenesis in Indolent Lymphomas

Activation-induced deaminase (AID) is required for somatic hypermutation in immunoglobulin genes, but also induces off-target mutations. Follicular lymphoma (FL) and chronic lymphocytic leukemia (CLL), the most frequent types of indolent B-cell tumors, are exposed to AID activity during lymphomagenesis. We designed a workflow integrating de novo mutational signatures extraction and fitting of COSMIC (Catalogue Of Somatic Mutations In Cancer) signatures, with tridimensional chromatin conformation data (Hi-C). We applied the workflow to exome sequencing data from lymphoma samples. In 33 FL and 30 CLL samples, 42% and 34% of the contextual mutations could be traced to a known AID motif. We demonstrate that both CLL and FL share mutational processes dominated by spontaneous deamination, failures in DNA repair, and AID activity. The processes had equiproportional distribution across active and nonactive chromatin compartments in CLL. In contrast, canonical AID activity and failures in DNA repair pathways in FL were significantly higher within the active chromatin compartment. Analysis of DNA repair genes revealed a higher prevalence of base excision repair gene mutations (p = 0.02) in FL than CLL. These data indicate that AID activity drives the genetic landscapes of FL and CLL. However, the final result of AID-induced mutagenesis differs between these lymphomas depending on chromatin compartmentalization and mutations in DNA repair pathways.

UnAIDed Class Switching in Activated B-Cells Reveals Intrinsic Features of a Self-Cleaving IgH Locus

Activation-induced deaminase (AID) is the major actor of immunoglobulin (Ig) gene diversification in germinal center B-cells. From its first description, it was considered as mandatory for class switch recombination (CSR), and this discovery initiated a long quest for all of the AID-interacting factors controlling its activity. The mechanisms focusing AID-mediated DNA lesions to given target sequences remain incompletely understood with regards the detailed characterization of optimal substrates in which cytidine deamination will lead to double strand breaks (DSBs) and chromosomal cleavage. In an effort to reconsider whether such CSR breaks absolutely require AID, we herein provide evidence, based on deep-sequencing approaches, showing that this dogma is not absolute in both human and mouse B lymphocytes. In activated B-cells from either AID-deficient mice or human AID-deficient patients, we report an intrinsic ability of the IgH locus to undergo “on-target” cleavage and subsequent synapsis of broken regions in conditions able to yield low-level CSR. DNA breaks occur in such conditions within the same repetitive S regions usually targeted by AID, but their repair follows a specific pathway with increased usage of microhomology-mediated repair. These data further demonstrate the role of AID machinery as not initiating de novo chromosomal cleavage but rather catalyzing a process which spontaneously initiates at low levels in an appropriately conformed IgH locus.

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.

Deep learning model of somatic hypermutation reveals importance of sequence context beyond targeting of AID and Polη hotspots

B-cells undergo somatic hypermutation (SHM) of the Immunoglobulin (Ig) variable region to generate high-affinity antibodies. SHM relies on the activity of activation-induced deaminase (AID), which mutates C>U preferentially targeting WRC (W=A/T, R=A/G) hotspots. Downstream mutations at WA Polymerase η hotspots contribute further mutations. Computational models of SHM can describe the probability of mutations essential for vaccine responses. Previous studies using short subsequences (k-mers) failed to explain divergent mutability for the same k-mer. We developed the DeepSHM (Deep learning on SHM) model using k-mers of size 5-21, improving accuracy over previous models. Interpretation of DeepSHM identified an extended DWRCT (D=A/G/T) motif with particularly high mutability. Increased mutability was further associated with lower surrounding G content. Our model also discovered a conserved AGYCTGGGGG (Y=C/T) motif within FW1 of IGHV3 family genes with unusually high T>G substitution rates. Thus, a wider sequence context increases predictive power and identifies novel features that drive mutational targeting.

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.

Auto-Antibody Production During Experimental Atherosclerosis in ApoE-/- Mice

Current models stipulate that B cells and antibodies function during atherosclerosis in two distinct ways based on antibody isotype, where IgM is protective and IgG is inflammatory. To examine this model, we generated ApoE-/- Aid-/- mice, which are unable to produce IgG antibodies due to the absence of activation-induced deaminase (AID) but maintain high plasma cholesterol due to the absence of apolipoprotein E (APOE). We saw a dramatic decrease in plaque formation in ApoE-/- Aid-/- mice compared to ApoE-/- mice. Rigorous analysis of serum antibodies revealed both ApoE-/- and ApoE-/- Aid-/- mice had substantially elevated titers of IgM antibodies compared to C57BL/6J controls, suggesting a more complex dynamic than previously described. Analysis of antigen specificity demonstrated that ApoE-/- Aid-/- mice had elevated titers of antibodies specific to malondialdehyde-oxidized low density lipoprotein (MDA-oxLDL), which has been shown to block macrophage recruitment into plaques. Conversely, ApoE-/- mice showed low levels of MDA-oxLDL specificity, but had antibodies specific to numerous self-proteins. We provide evidence for a hierarchical order of antibody specificity, where elevated levels of MDA-oxLDL specific IgM antibodies inhibit plaque formation. If the level of MDA-oxLDL specific IgM is insufficient, self-reactive IgM and IgG antibodies are generated against debris within the arterial plaque, resulting in increased inflammation and further plaque expansion.

Single cell clonal analysis identifies an AID-dependent pathway of plasma cell differentiation

Germinal centers (GC) are microstructures where B cells that have been activated by antigen can improve the affinity of their B cell receptors and differentiate into memory B cells (MBCs) or antibody secreting plasma cells. Activation Induced Deaminase (AID) initiates antibody diversification in GCs by somatic hypermutation and class switch recombination. Here we have addressed the role of AID in the terminal differentiation of GC B cells by combining single cell transcriptome and immunoglobulin clonal analysis in a mouse model that traces AID-experienced cells. We identified 8 transcriptional clusters that include dark zone and light zone GC subsets, plasmablasts/plasma cells (PB), 4 subsets of MBCs and a novel prePB subset, which shares the strongest clonal relationships with PBs. Mice lacking AID have various alterations in the size and expression profiles of these transcriptional clusters. We find that AID deficiency leads to a reduced proportion of prePB cells and severely impairs transitions between the prePB and the PB subsets. Thus, AID shapes the differentiation fate of GC B cells by enabling PB generation from a prePB state.

AID Contributes to Accelerated Disease Progression in the TCL1 Mouse Transplant Model for CLL

Adaptive somatic mutations conferring treatment resistance and accelerated disease progression is still a major problem in cancer therapy. Additionally in CLL, patients receiving novel, efficient drugs frequently become treatment refractory and eventually relapse. Activation-induced deaminase (AID) is a cytosine deaminase that catalyzes somatic hypermutation of genomic DNA at the immunoglobulin locus in activated B cells. As considerable off-target mutations by AID have been discerned in chronic lymphocytic leukemia, it is essential to investigate to which extent these mutations contribute to disease progression to estimate whether AID inhibition could counteract drug resistance mechanisms. In this study, we examined the TCL1 mouse model for CLL on an AID pro- and deficient background by comparing disease development and mutational landscapes. We provide evidence that AID contributes to the acquisition of somatic cancer-specific mutations also in the TCL1 model and accelerates CLL development particularly in the transplant setting. We conclude that AID is directly determining the fitness of the CLL clone, which prompts further studies to assess the effect of AID inhibition on the occurrence of drug resistance.