Interpretable deep learning reveals the role of an E-box motif in suppressing somatic hypermutation of AGCT motifs within human immunoglobulin variable regions.

Introduction: Somatic hypermutation (SHM) of immunoglobulin variable (V) regions by activation induced deaminase (AID) is essential for robust, long-term humoral immunity against pathogen and vaccine antigens. AID mutates cytosines preferentially within WRCH motifs (where W=A or T, R=A or G and H=A, C or T). However, it has been consistently observed that the mutability of WRCH motifs varies substantially, with large variations in mutation frequency even between multiple occurrences of the same motif within a single V region. This has led to the notion that the immediate sequence context of WRCH motifs contributes to mutability. Recent studies have highlighted the potential role of local DNA sequence features in promoting mutagenesis of AGCT, a commonly mutated WRCH motif. Intriguingly, AGCT motifs closer to 5' ends of V regions, within the framework 1 (FW1) sub-region1, mutate less frequently, suggesting an SHM-suppressing sequence context. Methods: Here, we systematically examined the basis of AGCT positional biases in human SHM datasets with DeepSHM, a machine-learning model designed to predict SHM patterns. This was combined with integrated gradients, an interpretability method, to interrogate the basis of DeepSHM predictions. Results: DeepSHM predicted the observed positional differences in mutation frequencies at AGCT motifs with high accuracy. For the conserved, lowly mutating AGCT motifs in FW1, integrated gradients predicted a large negative contribution of 5'C and 3'G flanking residues, suggesting that a CAGCTG context in this location was suppressive for SHM. CAGCTG is the recognition motif for E-box transcription factors, including E2A, which has been implicated in SHM. Indeed, we found a strong, inverse relationship between E-box motif fidelity and mutation frequency. Moreover, E2A was found to associate with the V region locale in two human B cell lines. Finally, analysis of human SHM datasets revealed that naturally occurring mutations in the 3'G flanking residues, which effectively ablate the E-box motif, were associated with a significantly increased rate of AGCT mutation. Discussion: Our results suggest an antagonistic relationship between mutation frequency and the binding of E-box factors like E2A at specific AGCT motif contexts and, therefore, highlight a new, suppressive mechanism regulating local SHM patterns in human V regions.

Fine-tuning spatial-temporal dynamics and surface receptor expression support plasma cell-intrinsic longevity.

Durable serological memory following vaccination is critically dependent on the production and survival of long-lived plasma cells (LLPCs). Yet, the factors that control LLPC specification and survival remain poorly resolved. Using intravital two-photon imaging, we find that in contrast to most plasma cells (PCs) in the bone marrow (BM), LLPCs are uniquely sessile and organized into clusters that are dependent on APRIL, an important survival factor. Using deep, bulk RNA sequencing, and surface protein flow-based phenotyping, we find that LLPCs express a unique transcriptome and phenotype compared to bulk PCs, fine-tuning expression of key cell surface molecules, CD93, CD81, CXCR4, CD326, CD44, and CD48, important for adhesion and homing. Conditional deletion of Cxcr4 in PCs following immunization leads to rapid mobilization from the BM, reduced survival of antigen-specific PCs, and ultimately accelerated decay of antibody titer. In naïve mice, the endogenous LLPCs BCR repertoire exhibits reduced diversity, reduced somatic mutations, and increased public clones and IgM isotypes, particularly in young mice, suggesting LLPC specification is non-random. As mice age, the BM PC compartment becomes enriched in LLPCs, which may outcompete and limit entry of new PCs into the LLPC niche and pool.

Evolvability of cancer-associated genes under APOBEC3A/B selection.

Evolvability is an emergent hallmark of cancer that depends on intra-tumor heterogeneity and genetic variation. Mutations generated by APOBEC3 contribute to genetic variation and tumor evolvability. However, the influence of APOBEC3 on the evolvability of the genome and its differential impact on cancer genes versus non-cancer genes remains unclear. Analyzing over 40,000 human protein-coding transcripts, we identified distinct distribution patterns of APOBEC3A/B TC motifs between cancer and non-cancer genes, suggesting unique associations with cancer. Studying a bat species with numerous APOBEC3 genes, we found distinct motif patterns in orthologs of cancer genes compared to non-cancer genes, as in humans, suggesting APOBEC3 evolution to reduce impacts on the genome rather than the converse. Simulations confirmed that APOBEC3-induced heterogeneity enhances cancer evolution through bimodal patterns of mutations in certain classes of genes. Our results suggest the bimodal distribution of APOBEC-induced mutations can significantly increase cancer heterogeneity.

Fine-tuning spatial-temporal dynamics and surface receptor expression support plasma cell-intrinsic longevity.

Durable serological memory following vaccination is critically dependent on the production and survival of long-lived plasma cells (LLPCs). Yet, the factors that control LLPC specification and survival remain poorly resolved. Using intra-vital two-photon imaging, we find that in contrast to most plasma cells in the bone marrow, LLPCs are uniquely sessile and organized into clusters that are dependent on April, an important survival factor. Using deep, bulk RNA sequencing, and surface protein flow-based phenotyping, we find that LLPCs express a unique transcriptome and proteome compared to bulk PCs, fine tuning expression of key cell surface molecules, CD93, CD81, CXCR4, CD326, CD44 and CD48, important for adhesion and homing, and phenotypically label LLPCs within mature PC pool. Conditional deletion of Cxcr4 in PCs following immunization leads to rapid mobilization from the BM, reduced survival of antigen-specific PCs, and ultimately accelerated decay of antibody titer. In naive mice, the endogenous LLPCs BCR repertoire exhibits reduced diversity, reduced somatic mutations, and increased public clones and IgM isotypes, particularly in young mice, suggesting LLPC specification is non-random. As mice age, the BM PC compartment becomes enriched in LLPCs, which may outcompete and limit entry of new PC into the LLPC niche and pool.

Modular cytosine base editing promotes epigenomic and genomic modifications.

Prokaryotic and eukaryotic adaptive immunity differ considerably. Yet, their fundamental mechanisms of gene editing via Cas9 and activation-induced deaminase (AID), respectively, can be conveniently complimentary. Cas9 is an RNA targeted dual nuclease expressed in several bacterial species. AID is a cytosine deaminase expressed in germinal centre B cells to mediate genomic antibody diversification. AID can also mediate epigenomic reprogramming via active DNA demethylation. It is known that sequence motifs, nucleic acid structures, and associated co-factors affect AID activity. But despite repeated attempts, deciphering AID's intrinsic catalytic activities and harnessing its targeted recruitment to DNA is still intractable. Even recent cytosine base editors are unable to fully recapitulate AID's genomic and epigenomic editing properties. Here, we describe the first instance of a modular AID-based editor that recapitulates the full spectrum of genomic and epigenomic editing activity. Our 'Swiss army knife' toolbox will help better understand AID biology per se as well as improve targeted genomic and epigenomic editing.

Evolvability of cancer-associated genes under APOBEC3A/B selection.

Evolvability is an emergent hallmark of cancer that depends on intra-tumor heterogeneity and, ultimately, genetic variation. Mutations generated by APOBEC3 cytidine deaminases can contribute to genetic variation and the consequences of APOBEC activation differ depending on the stage of cancer, with the most significant impact observed during the early stages. However, how APOBEC activity shapes evolutionary patterns of genes in the host genome and differential impacts on cancer-associated and non-cancer genes remain unclear. Analyzing over 40,000 human protein-coding transcripts, we identified distinct distribution patterns of APOBEC3A/B TC motifs between cancer-related genes and controls, suggesting unique associations with cancer. Studying a bat species with many more APOBEC3 genes, we found diverse motif patterns in orthologs of cancer genes compared to controls, similar to humans and suggesting APOBEC evolution to reduce impacts on the genome rather than the converse. Simulations confirmed that APOBEC-induced heterogeneity enhances cancer evolution, shaping clonal dynamics through bimodal introduction of mutations in certain classes of genes. Our results suggest that a major consequence of the bimodal distribution of APOBEC affects greater cancer heterogeneity.

Evolutionary potential of the monkeypox genome arising from interactions with human APOBEC3 enzymes.

APOBEC3, an enzyme subfamily that plays a role in virus restriction by generating mutations at particular DNA motifs or mutational 'hotspots', can drive viral mutagenesis with host-specific preferential hotspot mutations contributing to pathogen variation. While previous analysis of viral genomes from the 2022 Mpox (formerly Monkeypox) disease outbreak has shown a high frequency of C>T mutations at TC motifs, suggesting recent mutations are human APOBEC3-mediated, how emerging monkeypox virus (MPXV) strains will evolve as a consequence of APOBEC3-mediated mutations remains unknown. By measuring hotspot under-representation, depletion at synonymous sites, and a combination of the two, we analyzed APOBEC3-driven evolution in human poxvirus genomes, finding varying hotspot under-representation patterns. While the native poxvirus molluscum contagiosum exhibits a signature consistent with extensive coevolution with human APOBEC3, including depletion of TC hotspots, variola virus shows an intermediate effect consistent with ongoing evolution at the time of eradication. MPXV, likely the result of recent zoonosis, showed many genes with more TC hotspots than expected by chance (over-representation) and fewer GC hotspots than expected (under-representation). These results suggest the MPXV genome: (1) may have evolved in a host with a particular APOBEC GC hotspot preference, (2) has inverted terminal repeat (ITR) regions-which may be exposed to APOBEC3 for longer during viral replication-and longer genes likely to evolve faster, and therefore (3) has a heightened potential for future human APOBEC3-meditated evolution as the virus spreads in the human population. Our predictions of MPXV mutational potential can both help guide future vaccine development and identification of putative drug targets and add urgency to the task of containing human Mpox disease transmission and uncovering the ecology of the virus in its reservoir host.

Evolutionary potential of the monkeypox genome arising from interactions with human APOBEC3 enzymes.

APOBEC3, an enzyme subfamily that plays a role in virus restriction by generating mutations at particular DNA motifs or mutational "hotspots," can drive viral mutagenesis with host-specific preferential hotspot mutations contributing to pathogen variation. While previous analysis of viral genomes from the 2022 Mpox (formerly Monkeypox) disease outbreak has shown a high frequency of C>T mutations at T C motifs, suggesting recent mutations are human APOBEC3-mediated, how emerging monkeypox virus (MPXV) strains will evolve as a consequence of APOBEC3-mediated mutations remains unknown. By measuring hotspot under-representation, depletion at synonymous sites, and a combination of the two, we analyzed APOBEC3-driven evolution in human poxvirus genomes, finding varying hotspot under-representation patterns. While the native poxvirus molluscum contagiosum exhibits a signature consistent with extensive coevolution with human APOBEC3, including depletion of T C hotspots, variola virus shows an intermediate effect consistent with ongoing evolution at the time of eradication. MPXV, likely the result of recent zoonosis, showed many genes with more T C hotspots than expected by chance (over-representation) and fewer G C hotspots than expected (under-representation). These results suggest the MPXV genome: 1) may have evolved in a host with a particular APOBEC G C hotspot preference, 2) has inverted terminal repeat (ITR) regions -which may be exposed to APOBEC3 for longer during viral replication- and longer genes likely to evolve faster, and therefore 3) has a heightened potential for future human APOBEC3-meditated evolution as the virus spreads in the human population. Our predictions of MPXV mutational potential can both help guide future vaccine development and identification of putative drug targets and add urgency to the task of containing human Mpox disease transmission and uncovering the ecology of the virus in its reservoir host.

Beta human papillomavirus 8E6 promotes alternative end joining.

Double strand breaks (DSBs) are one of the most lethal DNA lesions in cells. The E6 protein of beta-human papillomavirus (HPV8 E6) impairs two critical DSB repair pathways: homologous recombination (HR) and non-homologous end joining (NHEJ). However, HPV8 E6 only delays DSB repair. How DSBs are repaired in cells with HPV8 E6 remains to be studied. We hypothesize that HPV8 E6 promotes a less commonly used DSB repair pathway, alternative end joining (Alt-EJ). Using CAS9-based Alt-EJ reporters, we show that HPV8 E6 promotes Alt-EJ. Further, using small molecule inhibitors, CRISPR/CAS9 gene knockout, and HPV8 E6 mutant, we find that HPV8 E6 promotes Alt-EJ by binding p300, an acetyltransferase that facilitates DSB repair by HR and NHEJ. At least some of this repair occurs through a subset of Alt-EJ known as polymerase theta dependent end joining. Finally, whole genome sequencing analysis showed HPV8 E6 caused an increased frequency of deletions bearing the microhomology signatures of Alt-EJ. This study fills the knowledge gap of how DSB is repaired in cells with HPV8 E6 and the mutagenic consequences of HPV8 E6 mediated p300 destabilization. Broadly, this study supports the hypothesis that beta-HPV promotes cancer formation by increasing genomic instability.

Role of EXO1 nuclease activity in genome maintenance, the immune response and tumor suppression in Exo1D173A mice.

DNA damage response pathways rely extensively on nuclease activity to process DNA intermediates. Exonuclease 1 (EXO1) is a pleiotropic evolutionary conserved DNA exonuclease involved in various DNA repair pathways, replication, antibody diversification, and meiosis. But, whether EXO1 facilitates these DNA metabolic processes through its enzymatic or scaffolding functions remains unclear. Here, we dissect the contribution of EXO1 enzymatic versus scaffolding activity by comparing Exo1DA/DA mice expressing a proven nuclease-dead mutant form of EXO1 to entirely EXO1-deficient Exo1-/- and EXO1 wild type Exo1+/+ mice. We show that Exo1DA/DA and Exo1-/- mice are compromised in canonical DNA repair processing, suggesting that the EXO1 enzymatic role is important for error-free DNA mismatch and double-strand break repair pathways. However, in non-canonical repair pathways, EXO1 appears to have a more nuanced function. Next-generation sequencing of heavy chain V region in B cells showed the mutation spectra of Exo1DA/DA mice to be intermediate between Exo1+/+ and Exo1-/- mice, suggesting that both catalytic and scaffolding roles of EXO1 are important for somatic hypermutation. Similarly, while overall class switch recombination in Exo1DA/DA and Exo1-/- mice was comparably defective, switch junction analysis suggests that EXO1 might fulfill an additional scaffolding function downstream of class switching. In contrast to Exo1-/- mice that are infertile, meiosis progressed normally in Exo1DA/DA and Exo1+/+ cohorts, indicating that a structural but not the nuclease function of EXO1 is critical for meiosis. However, both Exo1DA/DA and Exo1-/- mice displayed similar mortality and cancer predisposition profiles. Taken together, these data demonstrate that EXO1 has both scaffolding and enzymatic functions in distinct DNA repair processes and suggest a more composite and intricate role for EXO1 in DNA metabolic processes and disease.

Deep learning model of somatic hypermutation reveals importance of sequence context beyond hotspot targeting.

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 WWRCT 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 features that drive mutational targeting.

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.

Post-Transformation IGHV-IGHD-IGHJ Mutations in Chronic Lymphocytic Leukemia B Cells: Implications for Mutational Mechanisms and Impact on Clinical Course.

Analyses of IGHV gene mutations in chronic lymphocytic leukemia (CLL) have had a major impact on the prognostication and treatment of this disease. A hallmark of IGHV-mutation status is that it very rarely changes clonally over time. Nevertheless, targeted and deep DNA sequencing of IGHV-IGHD-IGHJ regions has revealed intraclonal heterogeneity. We used a DNA sequencing approach that achieves considerable depth and minimizes artefacts and amplification bias to identify IGHV-IGHD-IGHJ subclones in patients with prolonged temporal follow-up. Our findings extend previous studies, revealing intraclonal IGHV-IGHD-IGHJ diversification in almost all CLL clones. Also, they indicate that some subclones with additional IGHV-IGHD-IGHJ mutations can become a large fraction of the leukemic burden, reaching numerical criteria for monoclonal B-cell lymphocytosis. Notably, the occurrence and complexity of post-transformation IGHV-IGHD-IGHJ heterogeneity and the expansion of diversified subclones are similar among U-CLL and M-CLL patients. The molecular characteristics of the mutations present in the parental, clinically dominant CLL clone (CDC) differed from those developing post-transformation (post-CDC). Post-CDC mutations exhibit significantly lower fractions of mutations bearing signatures of activation induced deaminase (AID) and of error-prone repair by Polη, and most of the mutations were not ascribable to those enzymes. Additionally, post-CDC mutations displayed a lower percentage of nucleotide transitions compared with transversions that was also not like the action of AID. Finally, the post-CDC mutations led to significantly lower ratios of replacement to silent mutations in VH CDRs and higher ratios in VH FRs, distributions different from mutations found in normal B-cell subsets undergoing an AID-mediated process. Based on these findings, we propose that post-transformation mutations in CLL cells either reflect a dysfunctional standard somatic mutational process or point to the action of another mutational process not previously associated with IG V gene loci. If the former option is the case, post-CDC mutations could lead to a lesser dependence on antigen dependent BCR signaling and potentially a greater influence of off-target, non-IG genomic mutations. Alternatively, the latter activity could add a new stimulatory survival/growth advantage mediated by the BCR through structurally altered FRs, such as that occurring by superantigen binding and stimulation.

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.

Mutational pressure by host APOBEC3s more strongly affects genes expressed early in the lytic phase of herpes simplex virus-1 (HSV-1) and human polyomavirus (HPyV) infection.

Herpes-Simplex Virus 1 (HSV-1) infects most humans when they are young, sometimes with fatal consequences. Gene expression occurs in a temporal order upon lytic HSV-1 infection: immediate early (IE) genes are expressed, then early (E) genes, followed by late (L) genes. During this infection cycle, the HSV-1 genome has the potential for exposure to APOBEC3 (A3) proteins, a family of cytidine deaminases that cause C>U mutations on single-stranded DNA (ssDNA), often resulting in a C>T transition. We developed a computational model for the mutational pressure of A3 on the lytic cycle of HSV-1 to determine which viral kinetic gene class is most vulnerable to A3 mutations. Using in silico stochastic methods, we simulated the infectious cycle under varying intensities of A3 mutational pressure. We found that the IE and E genes are more vulnerable to A3 than L genes. We validated this model by analyzing the A3 evolutionary footprints in 25 HSV-1 isolates. We find that IE and E genes have evolved to underrepresent A3 hotspot motifs more so than L genes, consistent with greater selection pressure on IE and E genes. We extend this model to two-step infections, such as those of polyomavirus, and find that the same pattern holds for over 25 human Polyomavirus (HPyVs) genomes. Genes expressed earlier during infection are more vulnerable to mutations than those expressed later.

STAT3 imparts BRCAness by impairing homologous recombination repair in Epstein-Barr virus-transformed B lymphocytes.

Epstein-Barr virus (EBV) causes lymphomas and epithelial cell cancers. Though generally silent in B lymphocytes, this widely prevalent virus can cause endemic Burkitt lymphoma and post-transplant lymphoproliferative disorders/lymphomas in immunocompromised hosts. By learning how EBV breaches barriers to cell proliferation, we hope to undermine those strategies to treat EBV lymphomas and potentially other cancers. We had previously found that EBV, through activation of cellular STAT3 prevents phosphorylation of Chk1, and thereby, suppresses activation of the intra-S phase cell-cycle checkpoint, a potent barrier to oncogene-driven proliferation. This observation prompted us to examine the consequences on DNA repair since homologous recombination repair, the most error-free form, requires phosphoChk1. We now report that the defect in Chk1 phosphorylation also curtails RAD51 nucleation, and thereby, homologous recombination repair of DNA double strand breaks. The resulting reliance on error-prone microhomology-mediated end-joining (MMEJ) repair makes EBV-transformed cells susceptible to PARP inhibition and simultaneous accrual of genome-wide deletions and insertions resulting from synthesis-dependent MMEJ. Analysis of transcriptomic and drug susceptibility data from hundreds of cancer lines reveals a STAT3-dependent gene-set predictive of susceptibility of cancers to synthetic lethal PARP inhibition. These findings i) demonstrate how the tumor virus EBV re-shapes cellular DNA repair, ii) provide the first genome-wide evidence for insertions resulting from MMEJ in human cells, and iii) expand the range of cancers (EBV-related and -unrelated) that are likely to respond to synthetic lethal inhibitors given the high prevalence of cancers with constitutively active STAT3.

AID Overlapping and Polη Hotspots Are Key Features of Evolutionary Variation Within the Human Antibody Heavy Chain (IGHV) Genes.

Somatic hypermutation (SHM) of the immunoglobulin variable (IgV) loci is a key process in antibody affinity maturation. The enzyme activation-induced deaminase (AID), initiates SHM by creating C → U mismatches on single-stranded DNA (ssDNA). AID has preferential hotspot motif targets in the context of WRC/GYW (W = A/T, R = A/G, Y = C/T) and particularly at WGCW overlapping hotspots where hotspots appear opposite each other on both strands. Subsequent recruitment of the low-fidelity DNA repair enzyme, Polymerase eta (Polη), during mismatch repair, creates additional mutations at WA/TW sites. Although there are more than 50 functional immunoglobulin heavy chain variable (IGHV) segments in humans, the fundamental differences between these genes and their ability to respond to all possible foreign antigens is still poorly understood. To better understand this, we generated profiles of WGCW hotspots in each of the human IGHV genes and found the expected high frequency in complementarity determining regions (CDRs) that encode the antigen binding sites but also an unexpectedly high frequency of WGCW in certain framework (FW) sub-regions. Principal Components Analysis (PCA) of these overlapping AID hotspot profiles revealed that one major difference between IGHV families is the presence or absence of WGCW in a sub-region of FW3 sometimes referred to as "CDR4." Further differences between members of each family (e.g., IGHV1) are primarily determined by their WGCW densities in CDR1. We previously suggested that the co-localization of AID overlapping and Polη hotspots was associated with high mutability of certain IGHV sub-regions, such as the CDRs. To evaluate the importance of this feature, we extended the WGCW profiles, combining them with local densities of Polη (WA) hotspots, thus describing the co-localization of both types of hotspots across all IGHV genes. We also verified that co-localization is associated with higher mutability. PCA of the co-localization profiles showed CDR1 and CDR2 as being the main contributors to variance among IGHV genes, consistent with the importance of these sub-regions in antigen binding. Our results suggest that AID overlapping (WGCW) hotspots alone or in conjunction with Polη (WA/TW) hotspots are key features of evolutionary variation between IGHV genes.

Correlations in Somatic Hypermutation Between Sites in IGHV Genes Can Be Explained by Interactions Between AID and/or Polη Hotspots.

The somatic hypermutation (SHM) of Immunoglobulin (Ig) genes is a key process during antibody affinity maturation in B cells. The mutagenic enzyme activation induced deaminase (AID) is required for SHM and has a preference for WRC hotspots in DNA. Error-prone repair mechanisms acting downstream of AID introduce further mutations, including DNA polymerase eta (Polη), part of the non-canonical mismatch repair pathway (ncMMR), which preferentially generates mutations at WA hotspots. Previously proposed mechanistic models lead to a variety of predictions concerning interactions between hotspots, for example, how mutations in one hotspot will affect another hotspot. Using a large, high-quality, Ig repertoire sequencing dataset, we evaluated pairwise correlations between mutations site-by-site using an unbiased measure similar to mutual information which we termed "mutational association" (MA). Interactions are dominated by relatively strong correlations between nearby sites (short-range MAs), which can be almost entirely explained by interactions between overlapping hotspots for AID and/or Polη. We also found relatively weak dependencies between almost all sites throughout each gene (longer-range MAs), although these arise mostly as a statistical consequence of high pairwise mutation frequencies. The dominant short-range interactions are also highest within the most highly mutating IGHV sub-regions, such as the complementarity determining regions (CDRs), where there is a high hotspot density. Our results suggest that the hotspot preferences for AID and Polη have themselves evolved to allow for greater interactions between AID and/or Polη induced mutations.