Why do B cells mutate their immunoglobulin receptors?

Diversification of the VH3-53 immunoglobulin gene segment by somatic hypermutation results in neutralization of SARS-CoV-2 virus variants

COVID-19 induces re-circulating long-lived memory B cells (MBC) that, upon re-encounter with the pathogen, are induced to mount immunoglobulin responses. During convalescence, antibodies are subjected to affinity maturation, which enhances the antibody binding strength and generates new specificities that neutralize virus variants. Here, we performed a single-cell RNA sequencing analysis of spike-specific B cells from a SARS-CoV-2 convalescent subject. After COVID-19 vaccination, matured infection-induced MBC underwent recall and differentiated into plasmablasts. Furthermore, the transcriptomic profiles of newly activated B cells transiently shifted toward the ones of atypical and CXCR3+ B cells and several B-cell clonotypes massively expanded. We expressed monoclonal antibodies (mAbs) from all B-cell clones from the largest clonotype that used the VH3-53 gene segment. The in vitro analysis revealed that some somatic hypermutations enhanced the neutralization breadth of mAbs in a putatively stochastic manner. Thus, somatic hypermutation of B-cell clonotypes generates an anticipatory memory that can neutralize new virus variants.

Somatic hypermutation introduces bystander mutations that prepare SARS-CoV-2 antibodies for emerging variants

Somatic hypermutation (SHM) drives affinity maturation and continues over months in SARS-CoV-2-neutralizing antibodies (nAbs). However, several potent SARS-CoV-2 antibodies carry no or only a few mutations, leaving the question of how ongoing SHM affects neutralization unclear. Here, we reverted variable region mutations of 92 antibodies and tested their impact on SARS-CoV-2 binding and neutralization. Reverting higher numbers of mutations correlated with decreasing antibody functionality. However, for some antibodies, including antibodies of the public clonotype VH1-58, neutralization of Wu01 remained unaffected. Although mutations were dispensable for Wu01-induced VH1-58 antibodies to neutralize Alpha, Beta, and Delta variants, they were critical for Omicron BA.1/BA.2 neutralization. We exploited this knowledge to convert the clinical antibody tixagevimab into a BA.1/BA.2 neutralizer. These findings broaden our understanding of SHM as a mechanism that not only improves antibody responses during affinity maturation but also contributes to antibody diversification, thus increasing the chances of neutralizing viral escape variants.

450 million years in the making: mapping the evolutionary foundations of germinal centers

Germinal centers (GCs) are distinct microanatomical structures that form in the secondary lymphoid organs of endothermic vertebrates (i.e., mammals and some birds). Within GCs, B cells undergo a Darwinian selection process to identify clones which can respond to pathogen insult as well as affinity mature the B cell repertoire. The GC response ultimately generates memory B cells and bone marrow plasma cells which facilitate humoral immunological memory, the basis for successful vaccination programs. GCs have not been observed in the secondary lymphoid organs of ectothermic jawed vertebrates (i.e., fishes, reptiles, and amphibians). However, abundant research over the past decades has indicated these organisms can produce antigen specific B cell responses and some degree of affinity maturation. This review examines data demonstrating that the fundamentals of B cell selection may be more conserved across vertebrate phylogeny than previously anticipated. Further, research in both conventional mammalian model systems and comparative models raises the question of what evolutionary benefit GCs provide endotherms if they are seemingly unnecessary for generating the basic functional components of jawed vertebrate humoral adaptive immune responses.

Organized B cell sites in cartilaginous fishes reveal the evolutionary foundation of germinal centers

The absence of germinal centers (GCs) in cartilaginous fishes lies at odds with data showing that nurse sharks can produce robust antigen-specific responses and affinity mature their B cell repertoires. To investigate this apparent incongruity, we performed RNA sequencing on single nuclei, allowing us to characterize the cell types present in the nurse shark spleen, and RNAscope to provide in situ cellular resolution of key marker gene expression following immunization with R-phycoerythrin (PE). We tracked PE to the splenic follicles where it co-localizes with CXCR5high centrocyte-like B cells and a population of putative T follicular helper (Tfh) cells, surrounded by a peripheral ring of Ki67+ AID+ CXCR4+ centroblast-like B cells. Further, we reveal selection of mutations in B cell clones dissected from these follicles. We propose that the B cell sites identified here represent the evolutionary foundation of GCs, dating back to the jawed vertebrate ancestor.

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.

Translatability and transferability of in silico models: context of use switching to predict the effects of environmental chemicals on the immune system

Immunotoxicity hazard identification of chemicals aims to evaluate the potential for unintended effects of chemical exposure on the immune system. Perfluorinated alkylate substances (PFAS), such as perfluorooctane sulfonate (PFOS) and perfluorooctanoic acid (PFOA), are persistent, globally disseminated environmental contaminants known to be immunotoxic. Elevated PFAS exposure is associated with lower antibody responses to vaccinations in children and in adults. In addition, some studies have reported a correlation between PFAS levels in the body and lower resistance to disease, in other words an increased risk of infections or cancers.

In this context, modelling and simulation platforms could be used to simulate the human immune system with the aim to evaluate the adverse effects that immunotoxicants may have.

Here, we show the conditions under which a mathematical model developed for one purpose and application (e.g., in the pharmaceutical domain) can be successfully translated and transferred to another (e.g., in the chemicals domain) without undergoing significant adaptation. In particular, we demonstrate that the Universal Immune System Simulator was able to simulate the effects of PFAS on the immune system, introducing entities and new interactions that are biologically involved in the phenomenon. This also revealed a potentially exploitable pathway for assessing immunotoxicity through a computational model.

Large deletions in immunoglobulin genes are associated with a sustained absence of DNA Polymerase η

Somatic hypermutation of immunoglobulin genes is a highly mutagenic process that is B cell-specific and occurs during antigen-driven responses leading to antigen specificity and antibody affinity maturation. Mutations at the Ig locus are initiated by Activation-Induced cytidine Deaminase and are equally distributed at G/C and A/T bases. This requires the establishment of error-prone repair pathways involving the activity of several low fidelity DNA polymerases. In the physiological context, the G/C base pair mutations involve multiple error-prone DNA polymerases, while the generation of mutations at A/T base pairs depends exclusively on the activity of DNA polymerase η. Using two large cohorts of individuals with xeroderma pigmentosum variant (XP-V), we report that the pattern of mutations at Ig genes becomes highly enriched with large deletions. This observation is more striking for patients older than 50 years. We propose that the absence of Pol η allows the recruitment of other DNA polymerases that profoundly affect the Ig genomic landscape.

A Sensitivity Analysis Comparison of Three Models for the Dynamics of Germinal Centers

Germinal centers (GCs) are transient anatomical microenvironments where antibody affinity maturation and memory B cells generation takes place. In the past, models of Germinal Center (GC) dynamics have focused on understanding antibody affinity maturation rather than on the main mechanism(s) driving their rise-and-fall dynamics. Here, based on a population dynamics model core, we compare three mechanisms potentially responsible for this GC biphasic behavior dependent on follicular dendritic cell (FDC) maturation, follicular T helper (Tfh) cell maturation, and antigen depletion. Analyzing the kinetics of B and T cells, as well as its parameter sensitivities, we found that only the FDC-maturation-based model could describe realistic GC dynamics, whereas the simple Tfh-maturation and antigen-depletion mechanisms, as implemented here, could not. We also found that in all models the processes directly related to Tfh cell kinetics have the highest impact on GC dynamics. This suggests the existence of some still unknown mechanism(s) tuning GC dynamics by affecting Tfh cell response to proliferation-inducing stimuli.

Diagnosis of Xeroderma pigmentosum variant in a young patient with two novel mutations in the POLH gene

We describe the characterization of Xeroderma Pigmentosum variant (XPV) in a young Caucasian patient with phototype I, who exhibited a high sensitivity to sunburn and multiple cutaneous tumors at the age of 15 years. Two novel mutations in the POLH gene, which encodes the translesion DNA polymerase η, with loss of function due to two independent exon skippings, are reported to be associated as a compound heterozygous state in the patient. Western blot analysis performed on proteins from dermal fibroblasts derived from the patient and analysis of the mutation spectrum on immunoglobulin genes produced during the somatic hypermutation process in his memory B cells, show the total absence of translesion polymerase η activity in the patient. The total lack of Polη activity, necessary to bypass in an error-free manner UVR-induced pyrimidine dimers following sun exposure, explains the early unusual clinical appearance of this patient.

Integrating B cell lineage information into statistical tests for detecting selection in Ig sequences

Detecting selection in B cell Ig sequences is critical to understanding affinity maturation and can provide insights into Ag-driven selection in normal and pathologic immune responses. The most common sequence-based methods for detecting selection analyze the ratio of replacement and silent mutations using a binomial statistical analysis. However, these approaches have been criticized for low sensitivity. An alternative method is based on the analysis of lineage trees constructed from sets of clonally related Ig sequences. Several tree shape measures have been proposed as indicators of selection that can be statistically compared across cohorts. However, we show that tree shape analysis is confounded by underlying experimental factors that are difficult to control for in practice, including the sequencing depth and number of generations in each clone. Thus, although lineage tree shapes may reflect selection, their analysis alone is an unreliable measure of in vivo selection. To usefully capture the information provided by lineage trees, we propose a new method that applies the binomial statistical framework to mutations identified based on lineage tree structure. This hybrid method is able to detect selection with increased sensitivity in both simulated and experimental data sets. We anticipate that this approach will be especially useful in the analysis of large-scale Ig sequencing data sets generated by high-throughput sequencing technologies.

Quantifying selection in high-throughput Immunoglobulin sequencing data sets

High-throughput immunoglobulin sequencing promises new insights into the somatic hypermutation and antigen-driven selection processes that underlie B-cell affinity maturation and adaptive immunity. The ability to estimate positive and negative selection from these sequence data has broad applications not only for understanding the immune response to pathogens, but is also critical to determining the role of somatic hypermutation in autoimmunity and B-cell cancers. Here, we develop a statistical framework for Bayesian estimation of Antigen-driven SELectIoN (BASELINe) based on the analysis of somatic mutation patterns. Our approach represents a fundamental advance over previous methods by shifting the problem from one of simply detecting selection to one of quantifying selection. Along with providing a more intuitive means to assess and visualize selection, our approach allows, for the first time, comparative analysis between groups of sequences derived from different germline V(D)J segments. Application of this approach to next-generation sequencing data demonstrates different selection pressures for memory cells of different isotypes. This framework can easily be adapted to analyze other types of DNA mutation patterns resulting from a mutator that displays hot/cold-spots, substitution preference or other intrinsic biases.

Factors important in evolutionary shaping of immunoglobulin gene loci

Background

The extraordinary diversity characterizing the antibody repertoire is generated by both evolution and lymphocyte development. Much of this diversity is due to the existence of immunoglobulin (Ig) variable region gene segment libraries, which were diversified during evolution and, in higher vertebrates, are used in generating the combinatorial diversity of antibody genes. The aim of the present study was to address the following questions: What evolutionary parameters affect the size and structure of gene libraries? Are the number of genes in libraries of contemporary species, and the corresponding gene locus structure, a random result of evolutionary history, or have these properties been optimized with respect to individual or population fitness? If a larger number of genes or different genome structures do not increase the fitness, then the current structure is probably optimized.

Results

We used a simulation of variable region gene library evolution. We measured the effect of different parameters on gene library size and diversity, and the corresponding fitness. We found compensating relationships between parameters, which optimized Ig library size and diversity.

Conclusions

We conclude that contemporary species' Ig libraries have been optimized by evolution in terms of Ig sequence lengths, the number and diversity of Ig genes, and antibody-antigen affinities.

Error-prone DNA repair activity during somatic hypermutation in shark B lymphocytes

Sharks are representatives of the earliest vertebrates that possess an immune system utilizing V(D)J recombination to generate Ag receptors. Their Ab repertoire diversity is based in part on a somatic hypermutation process that introduces adjacent nucleotide substitutions of 2-5 bp. We have isolated mutant nonfunctional Ig rearrangements and intronic flank sequences to characterize the nonselected, intrinsic properties of this phenomenon; changes unique to shark were observed. Duplications and deletions were associated with N additions, suggesting participation of a DNA polymerase with some degree of template independence during the repair of DNA breaks initiated by activation-induced cytidine deaminase. Other mutations were consistent with some in vitro activities of mammalian translesion DNA polymerase η: tandem base substitutions, strand slippage, and small insertions/deletions. The nature of substitution patterns shows that DNA lesions at shark Ig genes recruit DNA repair factors with a species-specific repertoire of activities. We speculate that the tandem mutations are introduced by direct sequential misinsertions and that, in shark B cells, the mispairs tend to be extended rather than proofread. Despite extensive changes undergone by some mutants, the physical range of mutational activity remained restricted to VDJ and within the first 2-kb portion of the 6.8-kb J-C intron, perhaps a self-regulating aspect of activation-induced cytidine deaminase action that is conserved in evolution.

Toll-like receptors and B-cell receptors synergize to induce immunoglobulin class-switch DNA recombination: relevance to microbial antibody responses

Differentiation of naïve B cells, including immunoglobulin class-switch DNA recombination, is critical for the immune response and depends on the extensive integration of signals from the B-cell receptor (BCR), tumor necrosis factor (TNF) family members, Toll-like receptors (TLRs), and cytokine receptors. TLRs and BCR synergize to induce class-switch DNA recombination in T cell-dependent and T cell-independent antibody responses to microbial pathogens. BCR triggering together with simultaneous endosomal TLR engagement leads to enhanced B-cell differentiation and antibody responses. Te requirement of both BCR and TLR engagement would ensure appropriate antigen-specific activation in an infection. Co-stimulation of TLRs and BCR likely plays a significant role in anti-microbial antibody responses to contain pathogen loads until the T cell-dependent antibody responses peak. Furthermore, the temporal sequence of different signals is also critical for optimal B cell responses, as exemplified by the activation of B cells by initial TLR engagement, leading to the up-regulation of co-stimulatory CD80 and MCH-II receptors, which result in more efficient interactions with T cells, thereby enhancing the germinal center reaction and antibody affinity maturation. Overall, BCR and TLR stimulation and the integration with signals from the pathogen or immune cells and their products determine the ensuing B-cell antibody response.

Analysis of somatic hypermutation in X-linked hyper-IgM syndrome shows specific deficiencies in mutational targeting

Subjects with X-linked hyper-IgM syndrome (X-HIgM) have a markedly reduced frequency of CD27(+) memory B cells, and their Ig genes have a low level of somatic hypermutation (SHM). To analyze the nature of SHM in X-HIgM, we sequenced 209 nonproductive and 926 productive Ig heavy chain genes. In nonproductive rearrangements that were not subjected to selection, as well as productive rearrangements, most of the mutations were within targeted RGYW, WRCY, WA, or TW motifs (R = purine, Y = pyrimidine, and W = A or T). However, there was significantly decreased targeting of the hypermutable G in RGYW motifs. Moreover, the ratio of transitions to transversions was markedly increased compared with normal. Microarray analysis documented that specific genes involved in SHM, including activation-induced cytidine deaminase (AICDA) and uracil-DNA glycosylase (UNG2), were up-regulated in normal germinal center (GC) B cells, but not induced by CD40 ligation. Similar results were obtained from light chain rearrangements. These results indicate that in the absence of CD40-CD154 interactions, there is a marked reduction in SHM and, specifically, mutations of AICDA-targeted G residues in RGYW motifs along with a decrease in transversions normally related to UNG2 activity.

Antigen-driven selection in germinal centers as reflected by the shape characteristics of immunoglobulin gene lineage trees: a large-scale simulation study

During the immune response, the generation of memory B lymphocytes in germinal centers involves affinity maturation of the cells' antigen receptors, based on somatic hypermutation of receptor genes and antigen-driven selection of the resulting mutants. Affinity maturation is vital for immune protection, and is the basis of humoral immune learning and memory. Lineage trees of somatically hypermutated immunoglobulin genes often serve to qualitatively illustrate claims concerning the dynamics of affinity maturation in germinal centers. Here, we derive the quantitative relationships between parameters characterizing affinity maturation dynamics (proliferation, differentiation and mutation rates, initial affinity of the Ig to the antigen, and selection thresholds) and the mathematical properties of lineage trees, using a computer simulation which combines mathematical models for all mature B cell populations, stochastic models of hypermutation and selection, lineage tree generation and measurement of graphical tree characteristics. We identified seven key lineage tree properties, and found correlations of these with initial clone affinity and with the selection threshold. These two parameters were found to be the main factors affecting lineage tree shapes in both primary and secondary response trees. The results also confirm that recycling from centrocytes back to centroblasts is highly likely.

Stereotyped patterns of somatic hypermutation in subsets of patients with chronic lymphocytic leukemia: implications for the role of antigen selection in leukemogenesis

Somatic hypermutation (SHM) features in a series of 1967 immunoglobulin heavy chain gene (IGH) rearrangements obtained from patients with chronic lymphocytic leukemia (CLL) were examined and compared with IGH sequences from non-CLL B cells available in public databases. SHM analysis was performed for all 1290 CLL sequences in this cohort with less than 100% identity to germ line. At the cohort level, SHM patterns were typical of a canonical SHM process. However, important differences emerged from the analysis of certain subgroups of CLL sequences defined by: (1) IGHV gene usage, (2) presence of stereotyped heavy chain complementarity-determining region 3 (HCDR3) sequences, and (3) mutational load. Recurrent, “stereotyped” amino acid changes occurred across the entire IGHV region in CLL subsets carrying stereotyped HCDR3 sequences, especially those expressing the IGHV3-21 and IGHV4-34 genes. These mutations are underrepresented among non-CLL sequences and thus can be considered as CLL-biased. Furthermore, it was shown that even a low level of mutations may be functionally relevant, given that stereotyped amino acid changes can be found in subsets of minimally mutated cases. The precise targeting and distinctive features of somatic hypermutation (SHM) in selected subgroups of CLL patients provide further evidence for selection by specific antigenic element(s).

Impact of spiramycin treatment and gestational age on maturation of Toxoplasma gondii immunoglobulin G avidity in pregnant women

The objective of the present study was to investigate the maturation of immunoglobulin G (IgG) avidity after Toxoplasma gondii seroconversion during pregnancy and the factors that affect IgG avidity over time. The study used 309 serum samples from 117 women and a multiple linear mixed regression analysis to show the patterns of variation of IgG avidity throughout gestation. The IgG avidity ratios and the patterns of their evolution with time were quite diverse among the women and were statistically heterogeneous (P = 0.011); however, the trend was toward a statistically significant increase (P < 0.0001). On average, a 1.0167-fold increase was observed for each additional gestational week after the putative date of infection. At 12 weeks after putative infection (the expected IgG avidity maturation time), the mean avidity ratio was 16.6% (95% confidence interval, 15.4 to 17.9%). At all times, the avidity ratio remained significantly heterogeneous among the women (P < 0.05); for 95% of them, that ratio ranged from 7.8 to 35.3% at 12 weeks after putative infection. Maternal age at the putative time of infection did not influence the maturation of IgG avidity. However, on average, a 1.009-fold decrease (P = 0.03) in that avidity was observed for each additional week of gestational age before infection and a 1.03-fold increase (P = 0.0003) was observed for each additional week of delay to the onset of spiramycin treatment. The rate of increase in the avidity ratio was lower if infection occurred late in pregnancy and higher if the delay to treatment was long. This information cannot allow accurate determination of the delay since the time of infection. The present results provide support for interpretation of the assay and caution against overinterpretation.

AID-initiated purposeful mutations in immunoglobulin genes

Exposure brings risk to all living organisms. Using a remarkably effective strategy, higher vertebrates mitigate risk by mounting a complex and sophisticated immune response to counter the potentially toxic invasion by a virtually limitless army of chemical and biological antagonists. Mutations are almost always deleterious, but in the case of antibody diversification there are mutations occurring at hugely elevated rates within the variable (V) and switch regions (SR) of the immunoglobulin (Ig) genes that are responsible for binding to and neutralizing foreign antigens throughout the body. These mutations are truly purposeful. This chapter is centered on activation-induced cytidine deaminase (AID). AID is required for initiating somatic hypermutation (SHM) in the V regions and class switch recombination (CSR) in the SR portions of Ig genes. By converting C --> U, while transcription takes place, AID instigates a cascade of mutational events involving error-prone DNA polymerases, base excision and mismatch repair enzymes, and recombination pathways. Together, these processes culminate in highly mutated antibody genes and the B cells expressing antibodies that have achieved optimal antigenic binding undergo positive selection in germinal centers. We will discuss the biological role of AID in this complex process, primarily in terms of its biochemical properties in relation to SHM in vivo. The chapter also discusses recent advances in experimental methods to characterize antibody dynamics as a function of SHM to help elucidate the role that the AID-induced mutations play in tailoring molecular recognition. The emerging experimental techniques help to address long-standing conundrums concerning evolution-imposed constraints on antibody structure and function.

DNA polymerase eta is the sole contributor of A/T modifications during immunoglobulin gene hypermutation in the mouse

Mutations at A/T bases within immunoglobulin genes have been shown to be generated by a repair pathway involving the DNA-binding moiety of the mismatch repair complex constituted by the MSH2–MSH6 proteins, together with DNA polymerase η (pol η). However, residual A/T mutagenesis is still observed upon inactivation in the mouse of each of these factors, suggesting that the panel of activities involved might be more complex. We reported previously (Delbos, F., A. De Smet, A. Faili, S. Aoufouchi, J.-C. Weill, and C.-A. Reynaud. 2005. J. Exp. Med. 201:1191–1196) that residual A/T mutagenesis in pol η–deficient mice was likely contributed by another enzyme not normally involved in hypermutation, DNA polymerase κ, which is mobilized in the absence of the normal polymerase partner. We report the complete absence of A/T mutations in MSH2–pol η double-deficient mice, thus indicating that the residual A/T mutagenesis in MSH2-deficient mice is contributed by pol η, now recruited by uracil N-glycosylase, the second DNA repair pathway involved in hypermutation. We propose that this particular recruitment of pol η corresponds to a profound modification of the function of uracil glycosylase in the absence of the mismatch repair complex, suggesting that MSH2–MSH6 actively prevent uracil glycosylase from error-free repair during hypermutation. pol η thus appears to be the sole contributor of A/T mutations in the normal physiological context.