Analysis of 6912 unselected somatic hypermutations in human VDJ rearrangements reveals lack of strand specificity and correlation between phase II substitution rates and distance to the nearest 3' activation-induced cytidine deaminase target

Ancestral reconstruction reveals catalytic inactivation of activation-induced cytidine deaminase concomitant with cold water adaption in the Gadiformes bony fish

BACKGROUND

Antibody affinity maturation in vertebrates requires the enzyme activation-induced cytidine deaminase (AID) which initiates secondary antibody diversification by mutating the immunoglobulin loci. AID-driven antibody diversification is conserved across jawed vertebrates since bony and cartilaginous fish. Two exceptions have recently been reported, the Pipefish and Anglerfish, in which the AID-encoding aicda gene has been lost. Both cases are associated with unusual reproductive behavior, including male pregnancy and sexual parasitism. Several cold water fish in the Atlantic cod (Gadinae) family carry an aicda gene that encodes for a full-length enzyme but lack affinity-matured antibodies and rely on antibodies of broad antigenic specificity. Hence, we examined the functionality of their AID.

RESULTS

By combining genomics, transcriptomics, immune responsiveness, and functional enzymology of AID from 36 extant species, we demonstrate that AID of that Atlantic cod and related fish have extremely lethargic or no catalytic activity. Through ancestral reconstruction and functional enzymology of 71 AID enzymes, we show that this enzymatic inactivation likely took place relatively recently at the emergence of the true cod family (Gadidae) from their ancestral Gadiformes order. We show that this AID inactivation is not only concordant with the previously shown loss of key adaptive immune genes and expansion of innate and cell-based immune genes in the Gadiformes but is further reflected in the genomes of these fish in the form of loss of AID-favored sequence motifs in their immunoglobulin variable region genes.

CONCLUSIONS

Recent demonstrations of the loss of the aicda gene in two fish species challenge the paradigm that AID-driven secondary antibody diversification is absolutely conserved in jawed vertebrates. These species have unusual reproductive behaviors forming an evolutionary pressure for a certain loss of immunity to avoid tissue rejection. We report here an instance of catalytic inactivation and functional loss of AID rather than gene loss in a conventionally reproducing vertebrate. Our data suggest that an expanded innate immunity, in addition to lower pathogenic pressures in a cold environment relieved the pressure to maintain robust secondary antibody diversification. We suggest that in this unique scenario, the AID-mediated collateral genome-wide damage would form an evolutionary pressure to lose AID function.

Immunoglobulin somatic hypermutation in a defined biochemical system recapitulates affinity maturation and permits antibody optimization

We describe a purified biochemical system to produce monoclonal antibodies (Abs) in vitro using activation-induced deoxycytidine deaminase (AID) and DNA polymerase η (Polη) to diversify immunoglobulin variable gene (IgV) libraries within a phage display format. AID and Polη function during B-cell affinity maturation by catalyzing somatic hypermutation (SHM) of immunoglobulin variable genes (IgV) to generate high-affinity Abs. The IgV mutational motif specificities observed in vivo are conserved in vitro. IgV mutations occurred in antibody complementary determining regions (CDRs) and less frequently in framework (FW) regions. A unique feature of our system is the use of AID and Polη to perform repetitive affinity maturation on libraries reconstructed from a preceding selection step. We have obtained scFv Abs against human glucagon-like peptide-1 receptor (GLP-1R), a target in the treatment of type 2 diabetes, and VHH nanobodies targeting Fatty Acid Amide Hydrolase (FAAH), involved in chronic pain, and artemin, a neurotropic factor that regulates cold pain. A round of in vitro affinity maturation typically resulted in a 2- to 4-fold enhancement in Ab-Ag binding, demonstrating the utility of the system. We tested one of the affinity matured nanobodies and found that it reduced injury-induced cold pain in a mouse model.

Recent Progress in the Discovery and Development of Monoclonal Antibodies against Viral Infections

Monoclonal antibodies (mAbs), the new revolutionary class of medications, are fast becoming tools against various diseases thanks to a unique structure and function that allow them to bind highly specific targets or receptors. These specialized proteins can be produced in large quantities via the hybridoma technique introduced in 1975 or by means of modern technologies. Additional methods have been developed to generate mAbs with new biological properties such as humanized, chimeric, or murine. The inclusion of mAbs in therapeutic regimens is a major medical advance and will hopefully lead to significant improvements in infectious disease management. Since the first therapeutic mAb, muromonab-CD3, was approved by the U.S. Food and Drug Administration (FDA) in 1986, the list of approved mAbs and their clinical indications and applications have been proliferating. New technologies have been developed to modify the structure of mAbs, thereby increasing efficacy and improving delivery routes. Gene delivery technologies, such as non-viral synthetic plasmid DNA and messenger RNA vectors (DMabs or mRNA-encoded mAbs), built to express tailored mAb genes, might help overcome some of the challenges of mAb therapy, including production restrictions, cold-chain storage, transportation requirements, and expensive manufacturing and distribution processes. This paper reviews some of the recent developments in mAb discovery against viral infections and illustrates how mAbs can help to combat viral diseases and outbreaks.

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.

Antibody heavy chain CDR3 length-dependent usage of human IGHJ4 and IGHJ6 germline genes

Therapeutic antibody discovery using synthetic diversity has been proved productive, especially for target proteins not suitable for traditional animal immunization-based antibody discovery approaches. Recently, many lines of evidences suggest that the quality of synthetic diversity design limits the development success of synthetic antibody hits. The aim of our study is to understand the quality limitation and to properly address the challenges with a better design. Using VH3–23 as a model framework, we observed and quantitatively mapped CDR-H3 loop length-dependent usage of human IGHJ4 and IGHJ6 germline genes in the natural human immune repertoire. Skewed usage of DH2-JH6 and DH3-JH6 rearrangements was quantitatively determined in a CDR-H3 length-dependent manner in natural human antibodies with long CDR-H3 loops. Structural modeling suggests choices of JH help to stabilize antibody CDR-H3 loop and JH only partially contributes to the paratope. Our observations shed light on the design of next-generation synthetic diversity with improved probability of success.

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.

AID-RNA polymerase II transcription-dependent deamination of IgV DNA

Activation-induced deoxycytidine deaminase (AID) initiates somatic hypermutation (SHM) in immunoglobulin variable (IgV) genes to produce high-affinity antibodies. SHM requires IgV transcription by RNA polymerase II (Pol II). A eukaryotic transcription system including AID has not been reported previously. Here, we reconstitute AID-catalyzed deamination during Pol II transcription elongation in conjunction with DSIF transcription factor. C→T mutations occur at similar frequencies on non-transcribed strand (NTS) and transcribed strand (TS) DNA. In contrast, bacteriophage T7 Pol generates NTS mutations predominantly. AID-Pol II mutations are strongly favored in WRC and WGCW overlapping hot motifs (W = A or T, R = A or G) on both DNA strands. Single mutations occur on 70% of transcribed DNA clones. Mutations are correlated over a 15 nt distance in multiply mutated clones, suggesting that deaminations are catalyzed processively within a stalled or backtracked transcription bubble. Site-by-site comparisons for biochemical and human memory B-cell mutational spectra in an IGHV3-23*01 target show strongly favored deaminations occurring in the antigen-binding complementarity determining regions (CDR) compared to the framework regions (FW). By exhibiting consistency with B-cell SHM, our in vitro data suggest that biochemically defined reconstituted Pol II transcription systems can be used to investigate how, when and where AID is targeted.

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.

ASAP-SML: An antibody sequence analysis pipeline using statistical testing and machine learning

Antibodies are capable of potently and specifically binding individual antigens and, in some cases, disrupting their functions. The key challenge in generating antibody-based inhibitors is the lack of fundamental information relating sequences of antibodies to their unique properties as inhibitors. We develop a pipeline, Antibody Sequence Analysis Pipeline using Statistical testing and Machine Learning (ASAP-SML), to identify features that distinguish one set of antibody sequences from antibody sequences in a reference set. The pipeline extracts feature fingerprints from sequences. The fingerprints represent germline, CDR canonical structure, isoelectric point and frequent positional motifs. Machine learning and statistical significance testing techniques are applied to antibody sequences and extracted feature fingerprints to identify distinguishing feature values and combinations thereof. To demonstrate how it works, we applied the pipeline on sets of antibody sequences known to bind or inhibit the activities of matrix metalloproteinases (MMPs), a family of zinc-dependent enzymes that promote cancer progression and undesired inflammation under pathological conditions, against reference datasets that do not bind or inhibit MMPs. ASAP-SML identifies features and combinations of feature values found in the MMP-targeting sets that are distinct from those in the reference sets.

Genomic pattern of intratumor heterogeneity predicts the risk of progression in early stage diffuse large B-cell lymphoma

Intratumor heterogeneity (ITH) plays an important role in tumor development, metastases, recurrence and has impact on clinical diagnosis, therapeutic responses. In this study, we used mutant-allele tumor heterogeneity (MATH) through whole-exome sequencing data to evaluate ITH in early stage diffuse large B-cell lymphoma (DLBCL). A discovery set of 22 patients and a validation set of 35 patients with primary DLBCL which staged I or II were chosen from The Cancer Genome Atlas database and The Gene Index Project database, respectively. Then patients were divided into low and high MATH score groups according to the median expression level. As a result, higher MATH score displayed an increasing risk of progression compared with lower MATH score both in the discovery set (P = 0.045) and the validation set (P = 0.025). Further, the genomic pattern according to MATH demonstrated that mutation rates of immunoglobulin lambda locus, B-cell translocation gene and membrane-associated guanylate kinase were the sites with the highest mutation rate. Moreover, Gene enrichment analysis showed that immunoglobulin lambda constant 2 (IGLC2) which belongs to immunoglobulin lambda locus was only enriched in high MATH group. Besides, BTG2 and Caspase recruitment domain-containing protein 11 were cancer driver genes which listed on the top three high rate of mutation. Our study revealed the prognostic value of MATH and relevant genomic pattern in early stage DLBCL. Potential mechanisms on genetic level of how genomic aberration influence the ITH of DLBCL are provided. MATH based on whole genome sequencing may be feasible for deciding clinical treatment of DLBCL patients in the future researches.

Beyond Hot Spots: Biases in Antibody Somatic Hypermutation and Implications for Vaccine Design

The evolution of antibodies in an individual during an immune response by somatic hypermutation (SHM) is essential for the ability of the immune system to recognize and remove the diverse spectrum of antigens that may be encountered. These mutations are not produced at random; nucleotide motifs that result in increased or decreased rates of mutation were first reported in 1992. Newer models that estimate the propensity for mutation for every possible 5- or 7-nucleotide motif have emphasized the complexity of SHM targeting and suggested possible new hot spot motifs. Even with these fine-grained approaches, however, non-local context matters, and the mutations observed at a specific nucleotide motif varies between species and even by locus, gene segment, and position along the gene segment within a single species. An alternative method has been provided to further abstract away the molecular mechanisms underpinning SHM, prompted by evidence that certain stereotypical amino acid substitutions are favored at each position of a particular V gene. These "substitution profiles," whether obtained from a single B cell lineage or an entire repertoire, offer a simplified approach to predict which substitutions will be well-tolerated and which will be disfavored, without the need to consider path-dependent effects from neighboring positions. However, this comes at the cost of merging the effects of two distinct biological processes, the generation of mutations, and the selection acting on those mutations. Since selection is contingent on the particular antigens an individual has been exposed to, this suggests that SHM may have evolved to prefer mutations that are most likely to be useful against pathogens that have co-evolved with us. Alternatively, the ability to select favorable mutations may be strongly limited by the biases of SHM targeting. In either scenario, the sequence space explored by SHM is significantly limited and this consequently has profound implications for the rational design of vaccine strategies.

Analysis of 5518 unique, productively rearranged human VH3-23*01 gene sequences reveals CDR-H3 length-dependent usage of the IGHD2 gene family

Clear and accurate understanding of diversity in antibody complementarity-determining regions (CDRs) is critical for antibody discovery and engineering. Previous observations of antibody CDR-H3 diversity were based on analyzing available antibody sequences in the public databases. The results may not accurately reflect that of natural antibody repertoire due to erroneous species annotation and the presence of man-made CDR loop diversity in public antibody sequence databases. In this study, in a precisely controlled germline context, we explored the relationship between amino acid composition and CDR-H3 length using 5518 unique productively rearranged human VH3-23*01 gene sequences. CDR-H3 length-dependent usage of the Cys-Xn-Cys motif is reported here.

The Number of Overlapping AID Hotspots in Germline IGHV Genes Is Inversely Correlated with Mutation Frequency in Chronic Lymphocytic Leukemia

The targeting of mutations by Activation-Induced Deaminase (AID) is a key step in generating antibody diversity at the Immunoglobulin (Ig) loci but is also implicated in B-cell malignancies such as chronic lymphocytic leukemia (CLL). AID has previously been shown to preferentially deaminate WRC (W = A/T, R = A/G) hotspots. WGCW sites, which contain an overlapping WRC hotspot on both DNA strands, mutate at much higher frequency than single hotspots. Human Ig heavy chain (IGHV) genes differ in terms of WGCW numbers, ranging from 4 for IGHV3-48*03 to as many as 12 in IGHV1-69*01. An absence of V-region mutations in CLL patients ("IGHV unmutated", or U-CLL) is associated with a poorer prognosis compared to "IGHV mutated" (M-CLL) patients. The reasons for this difference are still unclear, but it has been noted that particular IGHV genes associate with U-CLL vs M-CLL. For example, patients with IGHV1-69 clones tend to be U-CLL with a poor prognosis, whereas patients with IGHV3-30 tend to be M-CLL and have a better prognosis. Another distinctive feature of CLL is that ~30% of (mostly poor prognosis) patients can be classified into "stereotyped" subsets, each defined by HCDR3 similarity, suggesting selection, possibly for a self-antigen. We analyzed >1000 IGHV genes from CLL patients and found a highly significant statistical relationship between the number of WGCW hotspots in the germline V-region and the observed mutation frequency in patients. However, paradoxically, this correlation was inverse, with V-regions with more WGCW hotspots being less likely to be mutated, i.e., more likely to be U-CLL. The number of WGCW hotspots in particular, are more strongly correlated with mutation frequency than either non-overlapping (WRC) hotspots or more general models of mutability derived from somatic hypermutation data. Furthermore, this correlation is not observed in sequences from the B cell repertoires of normal individuals and those with autoimmune diseases.

Automated analysis of high-throughput B-cell sequencing data reveals a high frequency of novel immunoglobulin V gene segment alleles

Individual variation in germline and expressed B-cell immunoglobulin (Ig) repertoires has been associated with aging, disease susceptibility, and differential response to infection and vaccination. Repertoire properties can now be studied at large-scale through next-generation sequencing of rearranged Ig genes. Accurate analysis of these repertoire-sequencing (Rep-Seq) data requires identifying the germline variable (V), diversity (D), and joining (J) gene segments used by each Ig sequence. Current V(D)J assignment methods work by aligning sequences to a database of known germline V(D)J segment alleles. However, existing databases are likely to be incomplete and novel polymorphisms are hard to differentiate from the frequent occurrence of somatic hypermutations in Ig sequences. Here we develop a Tool for Ig Genotype Elucidation via Rep-Seq (TIgGER). TIgGER analyzes mutation patterns in Rep-Seq data to identify novel V segment alleles, and also constructs a personalized germline database containing the specific set of alleles carried by a subject. This information is then used to improve the initial V segment assignments from existing tools, like IMGT/HighV-QUEST. The application of TIgGER to Rep-Seq data from seven subjects identified 11 novel V segment alleles, including at least one in every subject examined. These novel alleles constituted 13% of the total number of unique alleles in these subjects, and impacted 3% of V(D)J segment assignments. These results reinforce the highly polymorphic nature of human Ig V genes, and suggest that many novel alleles remain to be discovered. The integration of TIgGER into Rep-Seq processing pipelines will increase the accuracy of V segment assignments, thus improving B-cell repertoire analyses.

Overlapping hotspots in CDRs are critical sites for V region diversification

Activation-induced deaminase (AID) mediates the somatic hypermutation (SHM) of Ig variable (V) regions that is required for the affinity maturation of the antibody response. An intensive analysis of a published database of somatic hypermutations that arose in the IGHV3-23*01 human V region expressed in vivo by human memory B cells revealed that the focus of mutations in complementary determining region (CDR)1 and CDR2 coincided with a combination of overlapping AGCT hotspots, the absence of AID cold spots, and an abundance of polymerase eta hotspots. If the overlapping hotspots in the CDR1 or CDR2 did not undergo mutation, the frequency of mutations throughout the V region was reduced. To model this result, we examined the mutation of the human IGHV3-23*01 biochemically and in the endogenous heavy chain locus of Ramos B cells. Deep sequencing revealed that IGHV3-23*01 in Ramos cells accumulates AID-induced mutations primarily in the AGCT in CDR2, which was also the most frequent site of mutation in vivo. Replacing the overlapping hotspots in CDR1 and CDR2 with neutral or cold motifs resulted in a reduction in mutations within the modified motifs and, to some degree, throughout the V region. In addition, some of the overlapping hotspots in the CDRs were at sites in which replacement mutations could change the structure of the CDR loops. Our analysis suggests that the local sequence environment of the V region, and especially of the CDR1 and CDR2, is highly evolved to recruit mutations to key residues in the CDRs of the IgV region.

Sterile DJH rearrangements reveal that distance between gene segments on the human Ig H chain locus influences their ability to rearrange

Rearrangement of the Ig locus occurs in two steps. First, a JH gene is rearranged to a D gene followed by a VH gene rearranging to the DJH rearrangement. By next generation sequencing, we analyzed 9969 unique DJH rearrangements and 5919 unique VHDJH rearrangements obtained from peripheral blood B cells from 110 healthy adult donors. We found that DJH rearrangements and nonproductive VHDJH rearrangements share many features but differ significantly in their use of D genes and propensity for somatic hypermutation. In D to JH gene rearrangements, the D genes proximal to the JH locus are used more frequently than JH locus distal D genes, whereas VH locus proximal D genes were observed more frequently in nonproductive VHDJH rearrangements. We further demonstrate that the distance between VH, D, and JH gene segments influence their ability to rearrange within the human Ig locus.

Affinity maturation of B cells involves not only a few but a whole spectrum of relevant mutations

Affinity maturation of B lymphocytes within germinal centers involves both diversification of their B-cell receptors (BCRs) by somatic hypermutation (SHM) and a crucial receptor-mediated selection step. However, in contrast to recent advances in revealing the molecular mechanism of SHM, the fundamentals of the selection process are still poorly understood, i.e. it is often not clear how and how many mutations contribute to improving a BCR during the response against a given antigen. A general drawback in assessing the mutations relevant to the selection process is the difficult task of rating the relative contributions of selection and intrinsic biases to the experimentally observed mutation patterns of BCRs. The approach proposed here is premised on statistical comparison of the frequency distributions of nucleotide substitutions as observed in datasets of hypermutated BCRs against their frequency distribution expected under the null hypothesis of no selection. Thereby, we show that the spectrum of mutations relevant to maturation of canonical anti-(4-hydroxy-3-nitrophenyl)acetyl BCRs is much broader than previously acknowledged, going beyond the scope of single key mutations. Moreover, our results suggest that maturation not only involves selection by means of affinity but likewise expression and stabilization of BCRs.

Coupling mammalian cell surface display with somatic hypermutation for the discovery and maturation of human antibodies

A novel approach has been developed for the isolation and maturation of human antibodies that replicates key features of the adaptive immune system by coupling in vitro somatic hypermutation (SHM) with mammalian cell display. SHM is dependent on the action of the B cell specific enzyme, activation-induced cytidine deaminase (AID), and can be replicated in non-B cells through expression of recombinant AID. A library of human antibodies, based on germline V-gene segments with recombined human regions was used to isolate low-affinity antibodies to human β nerve growth factor (hβNGF). These antibodies, initially naïve to SHM, were subjected to AID-directed SHM in vitro and selected using the same mammalian cell display system, as illustrated by the maturation of one of the antibodies to low pM K(D). This approach overcomes many of the previous limitations of mammalian cell display, enabling direct selection and maturation of antibodies as full-length, glycosylated IgGs.

Recombinase-mediated cassette exchange as a novel method to study somatic hypermutation in Ramos cells

Activation-induced cytidine deaminase (AID) mediates the somatic hypermutation (SHM) of immunoglobulin (Ig) variable (V) regions that is required for the generation of antibody diversity and for the affinity maturation of the antibody response against infectious agents and toxic substances. AID preferentially targets WRC (W = A/T, R = A/G) hot spot motifs, particularly WGCW motifs that create overlapping hot spots on both strands. In order to gain a better understanding of the generation of antibody diversity and to create a platform for the in vitro generation of affinity-matured antibodies, we have established a system involving recombinase-mediated cassette exchange (RMCE) to replace the V region and its flanking sequences. This makes it possible to easily manipulate the sequence of the Ig gene within the endogenous heavy chain of the Ramos human Burkitt’s lymphoma cell line. Here we show that the newly integrated wild-type (WT) VH regions introduced by RMCE undergo SHM similarly to non-RMCE-modified Ramos cells. Most importantly, we have shown that introducing a cluster of WGCW motifs into the complementary determining region 2 (CDR2) of the human heavy chain V region significantly raised the mutation frequency and number of mutations per sequence compared to WT controls. Thus, we have demonstrated a novel platform in Ramos cells whereby we can easily and quickly manipulate the endogenous human VH region to further explore the regulation and targeting of SHM. This platform will be useful for generating human antibodies with changes in affinity and specificity in vitro.

An effective immune response requires a highly diverse repertoire of affinity-matured antibodies. Activation-induced cytidine deaminase (AID) is required for somatic hypermutation (SHM) of immunoglobulin (Ig) genes. Although a great deal has been learned about the regulation of AID, it remains unclear how it is preferentially targeted to particular motifs, to certain locations within the Ig gene and not to other highly expressed genes in the germinal center B cell. This is an important question because AID is highly mutagenic and is sometimes mistargeted to other highly expressed genes, including proto-oncogenes, leading to B cell lymphomas. Here we describe how we utilize recombinase-mediated cassette exchange (RMCE) to modify the sequence of the endogenous heavy chain locus in the Ramos Burkitt’s lymphoma cell line. This platform can be used to explore the regulation and targeting of SHM and to generate human antibodies with changes in affinity and specificity in vitro.

The concerted action of Msh2 and UNG stimulates somatic hypermutation at A . T base pairs

Mismatch repair plays an essential role in reducing the cellular mutation load. Paradoxically, proteins in this pathway produce A . T mutations during the somatic hypermutation of immunoglobulin genes. Although recent evidence implicates the translesional DNA polymerase eta in producing these mutations, it is unknown how this or other translesional polymerases are recruited to immunoglobulin genes, since these enzymes are not normally utilized in conventional mismatch repair. In this report, we demonstrate that A . T mutations were closely associated with transversion mutations at a deoxycytidine. Furthermore, deficiency in uracil-N-glycolase (UNG) or mismatch repair reduced this association. These data reveal a previously unknown interaction between the base excision and mismatch repair pathways and indicate that an abasic site generated by UNG within the mismatch repair tract recruits an error-prone polymerase, which then introduces A . T mutations. Our analysis further indicates that repair tracts typically are approximately 200 nucleotides long and that polymerase eta makes approximately 1 error per 300 T nucleotides. The concerted action of Msh2 and UNG in stimulating A . T mutations also may have implications for mutagenesis at sites of spontaneous cytidine deamination.

Exploration of factors affecting the onset and maturation course of follicular lymphoma through simulations of the germinal center

Genetic mutations frequently observed in human follicular lymphoma (FL) B-cells result in aberrant expression of the anti-apoptotic protein bcl-2 and surface immunoglobulins (Igs) which display one or more novel variable (V) region N-glycosylation motifs. In the present study, we develop a simulation model of the germinal center (GC) to explore how these mutations might influence the emergence and clonal expansion of key mutants which provoke FL development. The simulations employ a stochastic method for calculating the cellular dynamics, which incorporates actual IgV region sequences and a simplified hypermutation scheme. We first bring our simulations into agreement with experimental data for well-characterized normal and bcl-2(+) anti-hapten GC responses in mice to provide a model for understanding how bcl-2 expression leads to permissive selection and memory cell differentiation of weakly competitive B-cells. However, as bcl-2 expression in the GC alone is thought to be insufficient for FL development, we next monitor simulated IgV region mutations to determine the emergence times of key mutants displaying aberrant N-glycosylation motifs recurrently observed in human FL IgV regions. Simulations of 26 germline V(H) gene segments indicate that particular IgV regions have a dynamical selective advantage by virtue of the speed with which one or more of their key sites can generate N-glycosylation motifs upon hypermutation. Separate calculations attribute the high occurrence frequency of such IgV regions in FL to an ability to produce key mutants at a fast enough rate to overcome stochastic processes in the GC that hinder clonal expansion. Altogether, these simulations characterize three pathways for FL maturation through positively selected N-glycosylations, namely, via one of two key sites within germline V(H) region gene segments, or via a site in the third heavy chain complementarity-determining region (CDR-H3) that is generated from VDJ recombination.

DNA polymerase beta is able to repair breaks in switch regions and plays an inhibitory role during immunoglobulin class switch recombination

Immunoglobulin (Ig) class switch recombination (CSR) is initiated by activation-induced cytidine deaminase (AID), which converts cytosines to uracils in switch (S) regions. Subsequent excision of dU by uracil DNA glycosylase (UNG) of the base excision repair (BER) pathway is required to obtain double-strand break (DSB) intermediates for CSR. Since UNG normally initiates faithful repair, it is unclear how the AID-instigated S region lesions are converted into DSBs rather than correctly repaired by BER. Normally, DNA polymerase beta (Polbeta) would replace the dC deaminated by AID, leading to correct repair of the single-strand break, thereby preventing CSR. We address the question of whether Polbeta might be specifically down-regulated during CSR or inhibited from accessing the AID-instigated lesions, or whether the numerous AID-initiated S region lesions might simply overwhelm the BER capacity. We find that nuclear Polbeta levels are induced upon activation of splenic B cells to undergo CSR. When Polbeta(-/-) B cells are activated to switch in culture, they switch slightly better to IgG2a, IgG2b, and IgG3 and have more S region DSBs and mutations than wild-type controls. We conclude that Polbeta attempts to faithfully repair S region lesions but fails to repair them all.