Ig Enhancers Increase RNA Polymerase II Stalling at Somatic Hypermutation Target Sequences

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.

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.

RNA processing mechanisms contribute to genome organization and stability in B cells

RNA processing includes post-transcriptional mechanisms controlling RNA quality and quantity to ensure cellular homeostasis. Noncoding (nc) RNAs that are regulated by these dynamic processes may themselves fulfill effector and/or regulatory functions, and recent studies demonstrated the critical role of RNAs in organizing both chromatin and genome architectures. Furthermore, RNAs can threaten genome integrity when accumulating as DNA:RNA hybrids, but could also facilitate DNA repair depending on the molecular context. Therefore, by qualitatively and quantitatively fine-tuning RNAs, RNA processing contributes directly or indirectly to chromatin states, genome organization, and genome stability. B lymphocytes represent a unique model to study these interconnected mechanisms as they express ncRNAs transcribed from key specific sequences before undergoing physiological genetic remodeling processes, including V(D)J recombination, somatic hypermutation, and class switch recombination. RNA processing actors ensure the regulation and degradation of these ncRNAs for efficient DNA repair and immunoglobulin gene remodeling while failure leads to B cell development alterations, aberrant DNA repair, and pathological translocations. This review highlights how RNA processing mechanisms contribute to genome architecture and stability, with emphasis on their critical roles during B cell development, enabling physiological DNA remodeling while preventing lymphomagenesis.

The SV40 virus enhancer functions as a somatic hypermutation-targeting element with potential oncogenic activity

Simian virus 40 (SV40) is a monkey virus associated with several types of human cancers. SV40 is most frequently detected in mesotheliomas, brain and bone tumors and lymphomas, but the mechanism for SV40 tumorigenesis in humans is not clear. SV40 relative Merkel cell polyomavirus (MCPyV) causes Merkel cell carcinoma (MCC) in humans by expressing truncated large tumor antigen (LT) caused by APOBEC cytidine deaminase family enzymes induced mutations. AID (activation-induced cytidine deaminase), a member of the APOBEC family, is the initiator of the antibody diversification process known as somatic hypermutation (SHM) and its aberrant expression and targeting is a frequent source of lymphomagenesis. In this study, we investigated whether AID-induced mutations could cause truncation of SV40 LT. We demonstrate that the SV40 enhancer has strong SHM targeting activity in several cell types and that AID-induced mutations accumulate to SV40 LT in B cells and kidney cells and cause truncated LT expression in B cells. Our results argue that the ability of the SV40 enhancer to target SHM to LT is a potential source of LT truncation events in various cell types that could contribute to carcinogenesis.

Evolutionary analysis of a complete chicken genome

Microchromosomes are prevalent in nonmammalian vertebrates [P. D. Waters et al., Proc. Natl. Acad. Sci. U.S.A. 118 (2021)], but a few of them are missing in bird genome assemblies. Here, we present a new chicken reference genome containing all autosomes, a Z and a W chromosome, with all gaps closed except for the W. We identified ten small microchromosomes (termed dot chromosomes) with distinct sequence and epigenetic features, among which six were newly assembled. Those dot chromosomes exhibit extremely high GC content and a high level of DNA methylation and are enriched for housekeeping genes. The pericentromeric heterochromatin of dot chromosomes is disproportionately large and continues to expand with the proliferation of satellite DNA and testis-expressed genes. Our analyses revealed that the 41-bp CNM repeat frequently forms higher-order repeats (HORs) at the centromeres of acrocentric chromosomes. The centromere core regions where the kinetochore attaches often encompass telomeric sequence (TTAGGG)n, and in a one of the dot chromosomes, the centromere core recruits an endogenous retrovirus (ERV). We further demonstrate that the W chromosome shares some common features with dot chromosomes, having large arrays of hypermethylated tandem repeats. Finally, using the complete chicken chromosome models, we reconstructed a fine picture of chordate karyotype evolution, revealing frequent chromosomal fusions before and after vertebrate whole-genome duplications. Our sequence and epigenetic characterization of chicken chromosomes shed insights into the understanding of vertebrate genome evolution and chromosome biology.

The IgH Eµ-MAR regions promote UNG-dependent error-prone repair to optimize somatic hypermutation

Intoduction

Two scaffold/matrix attachment regions (5'- and 3'-MARs_Eµ ) flank the intronic core enhancer (cEµ) within the immunoglobulin heavy chain locus (IgH). Besides their conservation in mice and humans, the physiological role of MARs_Eµ is still unclear and their involvement in somatic hypermutation (SHM) has never been deeply evaluated.

Methods

Our study analyzed SHM and its transcriptional control in a mouse model devoid of MARs_Eµ , further combined to relevant models deficient for base excision repair and mismatch repair.

Results

We observed an inverted substitution pattern in of MARs_Eµ -deficient animals: SHM being decreased upstream from cEµ and increased downstream of it. Strikingly, the SHM defect induced by MARs_Eµ -deletion was accompanied by an increase of sense transcription of the IgH V region, excluding a direct transcription-coupled effect. Interestingly, by breeding to DNA repair-deficient backgrounds, we showed that the SHM defect, observed upstream from cEµ in this model, was not due to a decrease in AID deamination but rather the consequence of a defect in base excision repair-associated unfaithful repair process.

Discussion

Our study pointed out an unexpected "fence" function of MARs_Eµ regions in limiting the error-prone repair machinery to the variable region of Ig gene loci.

Temporal regulation of head-on transcription at replication initiation sites

Head-on (HO) collisions between the DNA replication machinery and RNA polymerase over R-loop forming sequences (RLFS) are genotoxic, leading to replication fork blockage and DNA breaks. Current models suggest that HO collisions are avoided through replication initiation site (RIS) positioning upstream of active genes, ensuring co-orientation of replication fork movement and genic transcription. However, this model does not account for pervasive transcription, or intragenic RIS. Moreover, pervasive transcription initiation and CG-rich DNA is a feature of RIS, suggesting that HO transcription units (HO TUs) capable of forming R-loops might occur. Through mining phased GRO-seq data, and developing an informatics strategy to stringently identify RIS, we demonstrate that HO TUs containing RLFS occur at RIS in MCF-7 cells, and are downregulated at the G1/S phase boundary. Our analysis reveals a novel spatiotemporal relationship between transcription and replication, and supports the idea that HO collisions are avoided through transcriptional regulatory mechanisms.

Oncogenic Merkel Cell Polyomavirus T Antigen Truncating Mutations are Mediated by APOBEC3 Activity in Merkel Cell Carcinoma

Merkel cell carcinoma (MCC) is an aggressive skin cancer, which is frequently caused by Merkel cell polyomavirus (MCPyV). Mutations of MCPyV tumor (T) antigens are major pathologic events of virus-positive (MCPyV+) MCCs, but their source is unclear. Activation-induced cytidine deaminase (AID)/APOBEC family cytidine deaminases contribute to antiviral immunity by mutating viral genomes and are potential carcinogenic mutators. We studied the contribution of AID/APOBEC cytidine deaminases to MCPyV large T (LT) truncation events. The MCPyV LT area in MCCs was enriched with cytosine-targeting mutations, and a strong APOBEC3 mutation signature was observed in MCC sequences. AICDA and APOBEC3 expression were detected in the Finnish MCC sample cohort, and LT expression correlated with APOBEC3H and APOBEC3G. Marginal but statistically significant somatic hypermutation targeting activity was detected in the MCPyV regulatory region. Our results suggest that APOBEC3 cytidine deaminases are a plausible cause of the LT truncating mutations in MCPyV+ MCC, while the role of AID in MCC carcinogenesis is unlikely.

Significance:

We uncover APOBEC3 mutation signature in MCPyV LT that reveals the likely cause of mutations underlying MCPyV+ MCC. We further reveal an expression pattern of APOBECs in a large Finnish MCC sample cohort. Thus, the findings presented here suggest a molecular mechanism underlying an aggressive carcinoma with poor prognosis.

OUP accepted manuscript

Activation-induced deaminase (AID) is a DNA-cytosine deaminase that mediates maturation of antibodies through somatic hypermutation and class-switch recombination. While it causes mutations in immunoglobulin heavy and light chain genes and strand breaks in the switch regions of the immunoglobulin heavy chain gene, it largely avoids causing such damage in the rest of the genome. To help understand targeting by human AID, we expressed it in repair-deficient Escherichia coli and mapped the created uracils in the genomic DNA using uracil pull-down and sequencing, UPD-seq. We found that both AID and the human APOBEC3A preferentially target tRNA genes and transcription start sites, but do not show preference for highly transcribed genes. Unlike A3A, AID did not show a strong replicative strand bias or a preference for hairpin loops. Overlapping uracilation peaks between these enzymes contained binding sites for a protein, FIS, that helps create topological domains in the E. coli genome. To confirm whether these findings were relevant to B cells, we examined mutations from lymphoma and leukemia genomes within AID-preferred sequences. These mutations also lacked replicative strand bias or a hairpin loop preference. We propose here a model for how AID avoids causing mutations in the single-stranded DNA found within replication forks.