HMCES protects immunoglobulin genes specifically from deletions during somatic hypermutation.

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

The evolutionary and functional significance of germline immunoglobulin gene variation.

The recombination between immunoglobulin (IG) gene segments determines an individual's naïve antibody repertoire and, consequently, (auto)antigen recognition. Emerging evidence suggests that mammalian IG germline variation impacts humoral immune responses associated with vaccination, infection, and autoimmunity - from the molecular level of epitope specificity, up to profound changes in the architecture of antibody repertoires. These links between IG germline variants and immunophenotype raise the question on the evolutionary causes and consequences of diversity within IG loci. We discuss why the extreme diversity in IG loci remains a mystery, why resolving this is important for the design of more effective vaccines and therapeutics, and how recent evidence from multiple lines of inquiry may help us do so.

Enhancers within the Ig V Gene Region Orchestrate Chromatin Topology and Regulate V Gene Rearrangement Frequency to Shape the B Cell Receptor Repertoire Specificities.

Effective Ab-mediated responses depend on a highly diverse Ab repertoire with the ability to bind a wide range of epitopes in disease-causing agents. The generation of this repertoire depends on the somatic recombination of the variable (V), diversity (D), and joining (J) genes in the Ig loci of developing B cells. It has been known for some time that individual V, D, and J gene segments rearrange at different frequencies, but the mechanisms behind this unequal V gene usage have not been well understood. However, recent work has revealed that newly described enhancers scattered throughout the V gene-containing portion of the Ig loci regulate the V gene recombination frequency in a regional manner. Deletion of three of these enhancers revealed that these elements exert many layers of control during V(D)J recombination, including long-range chromatin interactions, epigenetic milieu, chromatin accessibility, and compartmentalization.

Activation-Induced Cytidine Deaminase Expression in Human B Cell Precursors Is Essential for Central B Cell Tolerance.

Activation-induced cytidine deaminase (AID), the enzyme-mediating class-switch recombination (CSR) and somatic hypermutation (SHM) of immunoglobulin genes, is essential for the removal of developing autoreactive B cells. How AID mediates central B cell tolerance remains unknown. We report that AID enzymes were produced in a discrete population of immature B cells that expressed recombination-activating gene 2 (RAG2), suggesting that they undergo secondary recombination to edit autoreactive antibodies. However, most AID+ immature B cells lacked anti-apoptotic MCL-1 and were deleted by apoptosis. AID inhibition using lentiviral-encoded short hairpin (sh)RNA in B cells developing in humanized mice resulted in a failure to remove autoreactive clones. Hence, B cell intrinsic AID expression mediates central B cell tolerance potentially through its RAG-coupled genotoxic activity in self-reactive immature B cells.

Somatic Diversification of Rearranged Antibody Gene Segments by Intra- and Interchromosomal Templated Mutagenesis.

The ability of the humoral immune system to generate Abs capable of specifically binding a myriad of Ags is critically dependent on the somatic hypermutation program. This program induces both templated mutations (i.e., gene conversion) and untemplated mutations. In humans, somatic hypermutation is widely believed to result in untemplated point mutations. In this study, we demonstrate detection of large-scale templated events that occur in human memory B cells and circulating plasmablasts. We find that such mutations are templated intrachromosomally from IGHV genes and interchromosomally from IGHV pseudogenes as well as other homologous regions unrelated to IGHV genes. These same donor regions are used in multiple individuals, and they predominantly originate from chromosomes 14, 15, and 16. In addition, we find that exogenous sequences placed at the IgH locus, such as LAIR1, undergo templated mutagenesis and that homology appears to be the major determinant for donor choice. Furthermore, we find that donor tracts originate from areas in proximity with open chromatin, which are transcriptionally active, and are found in spatial proximity with the IgH locus during the germinal center reaction. These donor sequences are inserted into the Ig gene segment in association with overlapping activation-induced cytidine deaminase hotspots. Taken together, these studies suggest that diversity generated during the germinal center response is driven by untemplated point mutations as well as templated mutagenesis using local and distant regions of the genome.

Promoter Proximity Defines Mutation Window for VH and VΚ Genes Rearranged to Different J Genes.

Somatic hypermutation induced by activation-induced deaminase (AID) occurs at high densities between the Ig V gene promoter and intronic enhancer, which encompasses DNA encoding the rearranged V gene exon and J intron. It has been proposed that proximity between the promoter and enhancer defines the boundaries of mutation in V regions. However, depending on the J gene used, the distance between the promoter and enhancer is quite variable and may result in differential targeting around the V gene. To examine the effect of distance in mutation accumulation, we sequenced 320 clones containing different endogenous rearranged V genes in the IgH and Igκ loci from Peyer's patch B cells of mice. Clones were grouped by their use of different J genes. Distances between the V gene and enhancer ranged from ∼2.3 kb of intron DNA for rearrangements using J1, ∼2.0 kb for rearrangements using J2, ∼1.6 kb for rearrangements using J3 (H) or 4 (κ), and 1.1 kb for rearrangements using J4 (H) or 5 (κ). Strikingly, >90% of intron mutations occurred within 1 kb downstream of the J gene for both H and κ clones, regardless of which J gene was used. Thus, there is no evidence that the intron sequence or enhancer plays a role in determining the extent of mutation. The results indicate that V region intron mutations are targeted by their proximity to the promoter, suggesting they result from AID interactions with RNA polymerase II over a 1-kb region.

Genetic Characteristics of Non B Cell-Derived Immunoglobulin Genes.

It is widely acknowledged that immunoglobulins (Igs) are produced solely by B-lineage cells. The Ig gene is created by the rearrangement of a group of gene segments [variable (V), diversity (D), and joining (J) segments rearrangement, or V(D)J recombination], which results in the vast diversity of B cell-derived Ig responsible for recognising various antigens. Ig subsequently undergoes somatic hypermutation (SHM) and class switch recombination (CSR) after exposure to antigens, thus converting the low-affinity IgM to IgG, IgA, or IgE antibodies. IgM and IgD are primarily expressed in naïve B cells that have not been exposed to antigens, they do not undergo somatic hypermutation; hence, their variable region sequences remain the same as those in the germline. In contrast, IgG, IgA, and IgE are expressed in antigen-stimulated memory B cells or plasma cells, and thus, they often possess high-frequency mutations in their variable region sequences. Since the discovery that Ig can be produced by non-B cells, Qiu's group has investigated and compared the genetic characteristics of B cell-derived Ig and non-B cell-derived Ig. These findings demonstrated that non-B cell-derived Ig shares certain similarities with B cell-derived Ig in that the sequence of its constant region is identical to that of B cell-derived Ig, and its variable region is also strictly dependent on the rearrangement of V, D, and J gene segments. Moreover, akin to B cell-derived Ig, the V regions of IgM and IgD are rarely mutated, while IgG, IgA, and IgE produced by cancer cells are frequently mutated. However, the non-B cell-derived Ig V region sequence displays unique characteristics. (1) Unlike the vast diversity of B cell-derived Igs, non-B cell-derived Igs exhibit restricted diversity; cells from the same lineage always select the same V(D)J recombination patterns; (2) Both mRNA and proteins of RAG1/RAG2 recombinase have been detected in Ig positive cancer cell lines and normal tissues. But Ig recombination could also be found in RAG1-/- and RAG2-/- mice, suggesting that they are not necessary for the rearrangement of non-B cell-derived Igs. These features of non-B cell-derived Igs suggest a potentially undiscovered mechanism of V(D)J recombination, ligation, and SHM in non-B cells, which necessitates further investigation with advanced technology in molecular biology.

Deep-sequencing identification of the genomic targets of the cytidine deaminase AID and its cofactor RPA in B lymphocytes.

The cytidine deaminase AID hypermutates immunoglobulin genes but can also target oncogenes, leading to tumorigenesis. The extent of AID's promiscuity and its predilection for immunoglobulin genes are unknown. We report here that AID interacted broadly with promoter-proximal sequences associated with stalled polymerases and chromatin-activating marks. In contrast, genomic occupancy of replication protein A (RPA), an AID cofactor, was restricted to immunoglobulin genes. The recruitment of RPA to the immunoglobulin loci was facilitated by phosphorylation of AID at Ser38 and Thr140. We propose that stalled polymerases recruit AID, thereby resulting in low frequencies of hypermutation across the B cell genome. Efficient hypermutation and switch recombination required AID phosphorylation and correlated with recruitment of RPA. Our findings provide a rationale for the oncogenic role of AID in B cell malignancy.

An apoptosis-dependent checkpoint for autoimmunity in memory B and plasma cells.

B lymphocytes acquire self-reactivity as an unavoidable byproduct of antibody gene diversification in the bone marrow and in germinal centers (GCs). Autoreactive B cells emerging from the bone marrow are silenced in a series of well-defined checkpoints, but less is known about how self-reactivity that develops by somatic mutation in GCs is controlled. Here, we report the existence of an apoptosis-dependent tolerance checkpoint in post-GC B cells. Whereas defective GC B cell apoptosis has no measurable effect on autoantibody development, disruption of post-GC apoptosis results in accumulation of autoreactive memory B cells and plasma cells, antinuclear antibody production, and autoimmunity. The data presented shed light on mechanisms that regulate immune tolerance and the development of autoantibodies.

Role of STAT5 in controlling cell survival and immunoglobulin gene recombination during pro-B cell development.

STAT5 and interleukin 7 (IL-7) signaling are thought to control B lymphopoiesis by regulating the expression of key transcription factors and by activating variable (V(H)) gene segments at the immunoglobulin heavy-chain (Igh) locus. Using conditional mutagenesis to delete the gene encoding the transcription factor STAT5, we demonstrate that the development of pro-B cells was restored by transgenic expression of the prosurvival protein Bcl-2, which compensated for loss of the antiapoptotic protein Mcl-1. Expression of the genes encoding the B cell-specification factor EBF1 and the B cell-commitment protein Pax5 as well as V(H) gene recombination were normal in STAT5- or IL-7 receptor alpha-chain (IL-7Ralpha)-deficient pro-B cells rescued by Bcl-2. STAT5-expressing pro-B cells contained little or no active chromatin at most V(H) genes. In contrast, rearrangements of the immunoglobulin-kappa light-chain locus (Igk) were more abundant in STAT5- or IL-7Ralpha-deficient pro-B cells. Hence, STAT5 and IL-7 signaling control cell survival and the developmental ordering of immunoglobulin gene rearrangements by suppressing premature Igk recombination in pro-B cells.

B cell-specific and stimulation-responsive enhancers derepress Aicda by overcoming the effects of silencers.

Activation-induced cytidine deaminase (AID) is essential for the generation of antibody memory but also targets oncogenes, among other genes. We investigated the transcriptional regulation of Aicda (which encodes AID) in class switch-inducible CH12F3-2 cells and found that Aicda regulation involved derepression by several layers of positive regulatory elements in addition to the 5' promoter region. The 5' upstream region contained functional motifs for the response to signaling by cytokines, the ligand for the costimulatory molecule CD40 or stimuli that activated the transcription factor NF-kappaB. The first intron contained functional binding elements for the ubiquitous silencers c-Myb and E2f and for the B cell-specific activator Pax5 and E-box-binding proteins. Our results show that Aicda is regulated by the balance between B cell-specific and stimulation-responsive elements and ubiquitous silencers.

Lack of Evidence for a Substantial Rate of Templated Mutagenesis in B Cell Diversification.

BCR sequences diversify through mutations introduced by purpose-built cellular machinery. A recent paper has concluded that a "templated mutagenesis" process is a major contributor to somatic hypermutation and therefore Ig diversification in mice and humans. In this proposed process, mutations in the Ig locus are introduced by copying short segments from other Ig genes. If true, this would overturn decades of research on B cell diversification and would require a complete rewrite of computational methods to analyze B cell data for these species. In this paper, we re-evaluate the templated mutagenesis hypothesis. By applying the original inferential method using potential donor templates absent from B cell genomes, we obtain estimates of the methods' false positive rates. We find false positive rates of templated mutagenesis in murine and human Ig loci that are similar to or even higher than the original rate inferences, and by considering the bases used in substitution, we find evidence that if templated mutagenesis occurs, it is at a low rate. We also show that the statistically significant results in the original paper can easily result from a slight misspecification of the null model.

AID Phosphorylation Regulates Mismatch Repair-Dependent Class Switch Recombination and Affinity Maturation.

Activation-induced cytidine deaminase (AID) generates U:G mismatches in Ig genes that can be converted into untemplated mutations during somatic hypermutation or DNA double-strand breaks during class switch recombination (CSR). Null mutations in UNG and MSH2 demonstrate the complementary roles of the base excision repair (BER) and mismatch repair pathways, respectively, in CSR. Phosphorylation of AID at serine 38 was previously hypothesized to regulate BER during CSR, as the AID phosphorylation mutant, AID(S38A), cannot interact with APE1, a BER protein. Consistent with these findings, we observe a complete block in CSR in AIDS38A/S38AMSH2-/- mouse B cells that correlates with an impaired mutation frequency at 5'Sµ. Similarly, somatic hypermutation is almost negligible at the JH4 intron in AIDS38A/S38AMSH2-/- mouse B cells, and, consistent with this, NP-specific affinity maturation in AIDS38A/S38AMSH2-/- mice is not significantly elevated in response to NP-CGG immunization. Surprisingly, AIDS38A/S38AUNG-/- mouse B cells also cannot complete CSR or affinity maturation despite accumulating significant mutations in 5'Sµ as well as the JH4 intron. These data identify a novel role for phosphorylation of AID at serine 38 in mismatch repair-dependent CSR and affinity maturation.

Single-cell whole-genome sequencing reveals the functional landscape of somatic mutations in B lymphocytes across the human lifespan.

Accumulation of mutations in somatic cells has been implicated as a cause of aging since the 1950s. However, attempts to establish a causal relationship between somatic mutations and aging have been constrained by the lack of methods to directly identify mutational events in primary human tissues. Here we provide genome-wide mutation frequencies and spectra of human B lymphocytes from healthy individuals across the entire human lifespan using a highly accurate single-cell whole-genome sequencing method. We found that the number of somatic mutations increases from <500 per cell in newborns to >3,000 per cell in centenarians. We discovered mutational hotspot regions, some of which, as expected, were located at Ig genes associated with somatic hypermutation (SHM). B cell-specific mutation signatures associated with development, aging, or SHM were found. The SHM signature strongly correlated with the signature found in human B cell tumors, indicating that potential cancer-causing events are already present even in B cells of healthy individuals. We also identified multiple mutations in sequence features relevant to cellular function (i.e., transcribed genes and gene regulatory regions). Such mutations increased significantly during aging, but only at approximately one-half the rate of the genome average, indicating selection against mutations that impact B cell function. This full characterization of the landscape of somatic mutations in human B lymphocytes indicates that spontaneous somatic mutations accumulating with age can be deleterious and may contribute to both the increased risk for leukemia and the functional decline of B lymphocytes in the elderly.

Pronounced cohabitation of active immunoglobulin genes from three different chromosomes in transcription factories during maximal antibody synthesis.

To understand the relationships between nuclear organization and gene expression in a model system, we employed three-dimensional imaging and chromatin immunoprecipitation (ChIP)-chromosome conformation capture (3C) techniques to investigate the topographies of the immunoglobulin (Ig) genes and transcripts during B-cell development. Remarkably, in plasma cells, when antibody synthesis peaks, active Ig genes residing on three different chromosomes exhibit pronounced colocalizations in transcription factories, often near the nuclear periphery, and display trans-chromosomal enhancer interactions, and their transcripts frequently share interchromatin trafficking channels. Conceptually, these features of nuclear organization maximize coordinated transcriptional and transcript trafficking control for potentiating the optimal cytoplasmic assembly of the resulting translation products into protein multimers.

Immunoglobulins in teleosts.

Immunoglobulins are glycoproteins which are produced as membrane-bound receptors on B-cells or in a secreted form, known as antibodies. In teleosts, three immunoglobulin isotypes, IgM, IgT, and IgD, are present, each comprising two identical heavy and two identical light polypeptide chains. The basic mechanisms for generation of immunoglobulin diversity are similar in teleosts and higher vertebrates. The B-cell pre-immune repertoire is diversified by VDJ recombination, junctional flexibility, addition of nucleotides, and combinatorial association of light and heavy chains, while the post-immune repertoire undergoes somatic hypermutation during clonal expansion. Typically, the teleost immunoglobulin heavy chain gene complex has a modified translocon arrangement where the Dτ-Jτ-Cτ cluster of IgT is generally located between the variable heavy chain (VH) region and the Dµ/δ-Jµ/δ-Cµ-Cδ gene segments, or within the set of VH gene segments. However, multiple genome duplication and deletion events and loss of some individual genes through evolution has complicated the IgH gene organization. The IgH gene arrangement allows the expression of either IgT or IgM/IgD. Alternative splicing is responsible for the regulation of IgM/IgD expression and the secreted versus transmembrane forms of IgT, IgD, and IgM. The overall structure of IgM and IgT is usually conserved across species, whereas IgD has a large variety of structures. IgM is the main effector molecule in both systemic and mucosal immunity and shows a broad range of concentrations in different teleost species. Although IgM is usually present in higher concentrations under normal conditions, IgT is considered the main mucosal Ig.

Automated analysis of immunosequencing datasets reveals novel immunoglobulin D genes across diverse species.

Immunoglobulin genes are formed through V(D)J recombination, which joins the variable (V), diversity (D), and joining (J) germline genes. Since variations in germline genes have been linked to various diseases, personalized immunogenomics focuses on finding alleles of germline genes across various patients. Although reconstruction of V and J genes is a well-studied problem, the more challenging task of reconstructing D genes remained open until the IgScout algorithm was developed in 2019. In this work, we address limitations of IgScout by developing a probabilistic MINING-D algorithm for D gene reconstruction, apply it to hundreds of immunosequencing datasets from multiple species, and validate the newly inferred D genes by analyzing diverse whole genome sequencing datasets and haplotyping heterozygous V genes.

AncesTree: An interactive immunoglobulin lineage tree visualizer.

High-throughput sequencing of human immunoglobulin genes allows analysis of antibody repertoires and the reconstruction of clonal lineage evolution. The study of antibodies (Abs) affinity maturation is of specific interest to understand the generation of Abs with high affinity or broadly neutralizing activities. Moreover, phylogenic analysis enables the identification of the key somatic mutations required to achieve optimal antigen binding. The Immcantation framework provides a start-to-finish set of analytical methods for high-throughput adaptive immune receptor repertoire sequencing (AIRR-Seq; Rep-Seq) data. Furthermore, Immcantation's Change-O package has developed IgPhyML, an algorithm designed to build specifically immunoglobulin (Ig) phylogenic trees. Meanwhile Phylip, an algorithm that has been originally developed for applications in ecology and macroevolution, can also be used for the phylogenic reconstruction of antibodies maturation pathway. To complement Ig lineages made by IgPhyML or Dnaml (Phylip), we developed AncesTree, a graphic user interface (GUI) that aims to give researchers the opportunity to interactively explore antibodies clonal evolution. AncesTree displays interactive immunoglobulins phylogenic tree, Ig related mutations and sequence alignments using additional information coming from specialized antibody tools. The GUI is a Java standalone application allowing interaction with Ig tree that can run under Windows, Linux and Mac OS.

Identification of Two Nonrearranging IgSF Genes in Chicken Reveals a Novel Family of Putative Remnants of an Antigen Receptor Precursor.

In this study, we identified a pair of nonrearranging VJ-joined Ig superfamily genes, termed putative remnants of an Ag receptor precursor (PRARP) genes, in chicken. Both genes encode a single V-set Ig domain consisting of a canonical J-like segment and a potential immunoreceptor tyrosine-based inhibitory or switch motif in the cytoplasmic region. In vitro experiments showed that both genes were expressed at the cell surface as membrane proteins, and their recombinant products formed a monomer and a disulfide-linked homodimer or a heterodimer. These two genes were mainly expressed in B and T cells and were upregulated in response to stimulation with poly(I:C) in vitro and vaccination in vivo. Orthologs of PRARP have been identified in bony fish, amphibians, reptiles, and other birds, and a V-C1 structure similar to that of Ig or TCR chains was found in all these genes, with the exception of those in avian species, which appear to contain degenerated C1 domains or divergent Ig domains. Phylogenetic analyses suggested that the newly discovered genes do not belong to any known immune receptor family and appear to be a novel gene family. Further elucidation of the functions of PRARP and their origin might provide significant insights into the evolution of the immune system of jawed vertebrates.

Regulation of the Germinal Center Reaction and Somatic Hypermutation Dynamics by Homologous Recombination.

During somatic hypermutation (SHM) of Ig genes in germinal center B cells, lesions introduced by activation-induced cytidine deaminase are processed by multiple error-prone repair pathways. Although error-free repair by homologous recombination (HR) is crucial to prevent excessive DNA strand breakage at activation-induced cytidine deaminase off-target genes, its role at the hypermutating Ig locus in the germinal center is unexplored. Using B cell-specific inactivation of the critical HR factor Brca2, we detected decreased proliferation, survival, and thereby class switching of ex vivo-activated B cells. Intriguingly, an HR defect allowed for a germinal center reaction and affinity maturation in vivo, albeit at reduced amounts. Analysis of SHM revealed that a certain fraction of DNA lesions at C:G bp was indeed repaired in an error-free manner via Brca2 instead of being processed by error-prone translesion polymerases. By applying a novel pseudo-time in silico analysis of mutational processes, we found that the activity of A:T mutagenesis during SHM increased during a germinal center reaction, but this was in part defective in Brca2-deficient mice. These mutation pattern changes in Brca2-deficient B cells were mostly specific for the Ig V region, suggesting a local or time-dependent need for recombination repair to survive high rates of SHM and especially A:T mutagenesis.