The mutation patterns in B-cell immunoglobulin receptors reflect the influence of selection acting at multiple time-scales

During the several-week course of an immune response, B cells undergo a process of clonal expansion, somatic hypermutation of the immunoglobulin (Ig) genes and affinity-dependent selection. Over a lifetime, each B cell may participate in multiple rounds of affinity maturation as part of different immune responses. These two time-scales for selection are apparent in the structure of B-cell lineage trees, which often contain a ‘trunk’ consisting of mutations that are shared across all members of a clone, and several branches that form a ‘canopy’ consisting of mutations that are shared by a subset of clone members. The influence of affinity maturation on the B-cell population can be inferred by analysing the pattern of somatic mutations in the Ig. While global analysis of mutation patterns has shown evidence of strong selection pressures shaping the B-cell population, the effect of different time-scales of selection and diversification has not yet been studied. Analysis of B cells from blood samples of three healthy individuals identifies a range of clone sizes with lineage trees that can contain long trunks and canopies indicating the significant diversity introduced by the affinity maturation process. We here show that observed mutation patterns in the framework regions (FWRs) are determined by an almost purely purifying selection on both short and long time-scales. By contrast, complementarity determining regions (CDRs) are affected by a combination of purifying and antigen-driven positive selection on the short term, which leads to a net positive selection in the long term. In both the FWRs and CDRs, long-term selection is strongly dependent on the heavy chain variable gene family.     

Inferring processes underlying B-cell repertoire diversity

We quantify the VDJ recombination and somatic hypermutation processes in human B cells using probabilistic inference methods on high-throughput DNA sequence repertoires of human B-cell receptor heavy chains. Our analysis captures the statistical properties of the naive repertoire, first after its initial generation via VDJ recombination and then after selection for functionality. We also infer statistical properties of the somatic hypermutation machinery (exclusive of subsequent effects of selection). Our main results are the following: the B-cell repertoire is substantially more diverse than T-cell repertoires, owing to longer junctional insertions; sequences that pass initial selection are distinguished by having a higher probability of being generated in a VDJ recombination event; somatic hypermutations have a non-uniform distribution along the V gene that is well explained by an independent site model for the sequence context around the hypermutation site.     

Antibody engineering for the analysis of affinity maturation of an anti-hapten response.

The influence of structural variation, previously observed in a panel of V186.2 VH/V lambda 1-expressing anti-NP antibodies from the secondary response, on the affinity of these antibodies was examined by site-specific mutagenesis and recombinant antibody construction. A tryptophan----leucine exchange at position 33 in the VH segment of all but one of the high-affinity antibodies is the most frequently observed somatic mutation and by itself leads to a 10-fold higher affinity; all other somatic exchanges are irrelevant for affinity selection. In the single case of a high-affinity antibody without this common exchange, high affinity is mediated by a combination of mutations (including a one-codon deletion) in VH and the particular D-JH rearrangement carried by this antibody. The data indicate that the pattern of somatic diversification through hypermutation is shaped by affinity selection, but that only a single point mutation is available in the VH and the VL gene of lambda 1 chain-bearing anti-NP antibodies which by itself leads to an increase of hapten-binding affinity. Based on the analysis of two secondary response antibodies from which somatic mutations in VH and VL have been eliminated, it is also concluded that the recruitment of B cell clones into the pathway of hypermutation involves a mechanism which is not based upon affinity differences towards the antigen.     

Mechanisms of B cell autoimmunity in SLE

Systemic lupus erythematosus (SLE) is a systemic autoimmune disease that is known to be associated with polyclonal B-cell hyperreactivity. The underlying causes of the diffuse B-cell over-reactivity are unclear, but potential candidates include (a) intrinsic hyper-reactivity leading to polyclonal B-cell activation with disturbed activation thresholds and ineffective negative selection; (b) lack of immunoregulatory functions; (c) secondary effects of an overactive inflammatory environment, such as overactive germinal center and ectopic follicular activity; and/or (d) disturbed cytokine production by non-B immune cells. These mechanisms are not mutually exclusive and may operate to varying extents and at varying times in SLE. Phenotypic and molecular studies as well as the results of recent clinical trials have begun to provide new insights to address these possibilities. Of importance, new information has made it possible to distinguish between the contribution played by abnormalities in central checkpoints that could lead to a pre-immune repertoire enriched in autoreactive B cells, on the one hand, and the possibility that autoimmunity arises in the periphery from somatic hypermutation and abnormal selection during T cell-dependent B-cell responses on the other. There is an intriguing possibility that apoptotic material bound to the surface of follicular dendritic cells positively selects autoreactive B cells that arise from non-autoreactive B-cell precursors as a result of somatic hypermutation and thereby promotes the peripheral emergence of autoimmunity.     

Lack of detectable somatic hypermutation in the V region of the Ig H chain gene of a human chronic B lymphocytic leukemia

Monoclonal human B cell tumors are a model system for the study of somatic hypermutation of the Ig genes of humans. It was previously shown that a number of B cell lymphomas exhibited striking V region point mutation, hypothesized to result from the somatic hypermutation mechanism. In this study we have extended the analysis to chronic lymphocytic leukemia. We have cloned and sequenced the productive Vh representing five different cells from a monoclonal chronic lymphocytic leukemia. All five Vh sequences were identical. Therefore, the Vh region in this leukemia was not the subject of detectable somatic mutation. These data suggest that chronic lymphocytic leukemia might lack the mechanism for somatic hypermutation and represent a stage of normal B lymphocyte differentiation in which the somatic hypermutation mechanism is not active.     

Estimating hypermutation rates from clonal tree data

To understand the mechanisms underlying the varying patterns of mutations that occur during immune and autoimmune responses, estimates of the somatic hypermutation rate are critical. However, despite its significance, precise estimates of the mutation rate do not currently exist. Microdissection studies of mutating B cell clones provide an opportunity to measure this rate more accurately than previously possible. Each microdissection provides a number of clonally related sequences that, through the analysis of shared mutations, can be genealogically related to each other. The shape of these clonal trees is influenced by many processes, including the hypermutation rate. We have developed two different methods to estimate the mutation rate based on these data. These methods are applied to two sets of experimental data, one from an autoimmune response and one from the antihapten response to (4-hydroxy-3-nitrophenyl)acetyl (NP). Comparable mutation rates are estimated for both responses, 0.7-0.9 x 10(-3) and 0.9-1.1 x 10(-3) bp(-1) division(-1) for the autoimmune and NP responses, respectively. In addition to comparing the results of the two procedures, we investigate the effect on our estimate of assumptions, such as the fraction of lethal mutations.     

The V<Subscript>H</Subscript> gene repertoire of splenic B cells and somatic hypermutation in systemic lupus erythematosus

In systemic lupus erythematosus (SLE) it has been hypothesized that self-reactive B cells arise from virgin B cells that express low-affinity, nonpathogenic germline V genes that are cross-reactive for self and microbial antigens, which convert to high-affinity autoantibodies via somatic hypermutation. The aim of the present study was to determine whether the V_H family repertoire and pattern of somatic hypermutation in germinal centre (GC) B cells deviates from normal in SLE. Rearranged immunoglobulin V_H genes were cloned and sequenced from GCs of a SLE patient's spleen. From these data the GC V gene repertoire and the pattern of somatic mutation during the proliferation of B-cell clones were determined. The results highlighted a bias in V_H5 gene family usage, previously unreported in SLE, and under-representation of the V_H1 family, which is expressed in 20–30% of IgM⁺ B cells of healthy adults and confirmed a defect in negative selection. This is the first study of the splenic GC response in human SLE.     

Maturation of the immune response in germinal centers.

Germinal centers develop in peripheral lymphatic tissue during the primary immune response and may play a crucial role in affinity maturation. We have compared the diversification of the antigen-specific repertoire of B cells, both from within and from outside the germinal centers, during the murine response to 2-phenyloxazolone (phOx). By sequencing V kappa Ox1 L-chains characteristic of phOx-specific antibodies, we show that somatic mutations accumulate in germinal center B cells and that a mutation conferring high affinity binding is found with increasing frequency. An analysis of V/D/J rearrangements suggests that this mutation occurred independently in many B cells, which were then preferentially expanded. We conclude that, although the hypermutation mechanism may be activated before germinal centers develop, affinity maturation by hypermutation and selection takes place in the germinal centers.     

Trends in antibody sequence changes during the somatic hypermutation process

Probable germline gene sequences from thousands of aligned mature Ab sequences are inferred using simple computational matching to known V(D)J genes. Comparison of the germline to mature sequences in a structural region-dependent fashion allows insights into the methods that nature uses to mature Abs during the somatic hypermutation process. Four factors determine the residue type mutation patterns: biases in the germline, accessibility from single base permutations, location of mutation hotspots, and functional pressures during selection. Germline repertoires at positions that commonly contact the Ag are biased with tyrosine, serine, and tryptophan. These residue types have a high tendency to be present in mutation hotspot motifs, and their abundance is decreased during maturation by a net conversion to other types. The heavy use of tyrosines on mature Ab interfaces is thus a reflection of the germline composition rather than being due to selection during maturation. Potentially stabilizing changes such as increased proline usage and a small number of double cysteine mutations capable of forming disulfide bonds are ascribed to somatic hypermutation. Histidine is the only residue type for which usage increases in each of the interface, core, and surface regions. The net overall effect is a conversion from residue types that could provide nonspecific initial binding into a diversity of types that improve affinity and stability. Average mutation probabilities are approximately 4% for core residues, approximately 5% for surface residues, and approximately 12% for residues in common Ag-contacting positions, excepting the those coded by the D gene.     

Quantifying evolutionary constraints on B-cell affinity maturation

The antibody repertoire of each individual is continuously updated by the evolutionary process of B-cell receptor (BCR) mutation and selection. It has recently become possible to gain detailed information concerning this process through high-throughput sequencing. Here, we develop modern statistical molecular evolution methods for the analysis of B-cell sequence data, and then apply them to a very deep short-read dataset of BCRs. We find that the substitution process is conserved across individuals but varies significantly across gene segments. We investigate selection on BCRs using a novel method that side-steps the difficulties encountered by previous work in differentiating between selection and motif-driven mutation; this is done through stochastic mapping and empirical Bayes estimators that compare the evolution of in-frame and out-of-frame rearrangements. We use this new method to derive a per-residue map of selection, which provides a more nuanced view of the constraints on framework and variable regions.     

Codon bias and plasticity in immunoglobulins

Immunoglobulin genes experience Darwinian evolution twice. In addition to the germline evolution all genes experience, immunoglobulins are subjected, upon exposure to antigen, to somatic hypermutation. This is accompanied by selection for high affinity to the eliciting antigen and frequently results in a significant increase in the specificity of the responding population. The hypermutation mechanism displays a strong sequence specificity. Thus arises the opportunity to manipulate codon bias in a site-specific manner so as to direct hypermutation to those parts of the gene that encode the antigen-binding portions of the molecule and away from those that encode the structurally conserved regions. This segregation of mutability would clearly be advantageous; it would enhance the generation of potentially useful variants while keeping mutational loss to acceptably low levels. But it is not clear that the advantage gained would be large enough to produce a measurable effect within the background stochasticity of the evolutionary process. I have performed a pair of statistical tests to determine whether site- specific codon bias in human immunoglobulin genes is correlated with the sequence specificity of the somatic mutation mechanism. The sequence specificity of the mutator was determined by analysis of a database of published immunoglobulin intron sequences that had experienced somatic mutation but not selection. The site-specific codon bias was determined by analysis of published sequences of human germline immunoglobulin V genes. Both tests strongly suggest that evolution has acted to enhance the plasticity of immunoglobulin genes under somatic hypermutation.      

The targeting of somatic hypermutation closely resembles that of meiotic mutation

We have compared the microsequence specificity of mutations introduced during somatic hypermutation (SH) and those introduced meiotically during neutral evolution. We have minimized the effects of selection by studying nonproductive (hence unselected) Ig V region genes for somatic mutations and processed pseudogenes for meiotic mutations. We find that the two sets of patterns are very similar: the mutabilities of nucleotide triplets are positively correlated between the somatic and meiotic sets. The major differences that do exist fall into three distinct categories: 1) The mutability is sharply higher at CG dinucleotides under meiotic but not somatic mutation. 2) The complementary triplets AGC and GCT are much more mutable under somatic than under meiotic mutation. 3) Triplets of the form WAN (W = T or A) are uniformly more mutable under somatic than under meiotic mutation. Nevertheless, the relative mutabilities both within this set and within the SAN (S = G or C) triplets are highly correlated with those under meiotic mutation. We also find that the somatic triplet specificity is strongly symmetric under strand exchange for A/T triplets as well as for G/C triplets in spite of the strong predominance of A over T mutations. Thus, we suggest that somatic mutation has at least two distinct components: one that specifically targets AGC/GCT triplets and another that acts as true catalysis of meiotic mutation.     

Method for generation of human hyperdiversified antibody fragment library

The selection of antibody fragments from libraries using in vitro screening technologies has proven to be a very good alternative to the classical hybridoma technology, and has overcome the laborious process of antibody humanization. However, the complexity of the library is critical in the probability of being able to directly isolate a high affinity antibody specific to a target. We report a method to make hyperdiversified antibody fragment libraries, based on human immunoglobulin variable genes mimicking the somatic hypermutation process. This mutagenesis technology, MutaGen, was used for the first time on the entire variable domain (frameworks and CDRs) of large repertoires of human variable antibody domains. Our MutaGen process uses low-fidelity human polymerases, known as mutases, suggested to be involved in the somatic hypermutation process of immunoglobulin genes. Depending on the mutases used, we generated complementary mutation patterns with randomly distributed mutations. The libraries were generated with an average of 1.8 mutations per 100 amino acids. The hyperdiversified antibody fragment libraries constructed with our process should enable the selection of antibody fragments specific to virtually any target.     

Complex Patterns of Chromosome 11 Aberrations in Myeloid Malignancies Target CBL,MLL, DDB1 and LMO2

Exome sequencing of primary tumors identifies complex somatic mutation patterns. Assignment of relevance of individual somatic mutations is difficult and poses the next challenge for interpretation of next generation sequencing data. Here we present an approach how exome sequencing in combination with SNP microarray data may identify targets of chromosomal aberrations in myeloid malignancies. The rationale of this approach is that hotspots of chromosomal aberrations might also harbor point mutations in the target genes of deletions, gains or uniparental disomies (UPDs). Chromosome 11 is a frequent target of lesions in myeloid malignancies. Therefore, we studied chromosome 11 in a total of 813 samples from 773 individual patients with different myeloid malignancies by SNP microarrays and complemented the data with exome sequencing in selected cases exhibiting chromosome 11 defects. We found gains, losses and UPDs of chromosome 11 in 52 of the 813 samples (6.4%). Chromosome 11q UPDs frequently associated with mutations of CBL. In one patient the 11qUPD amplified somatic mutations in both CBL and the DNA repair gene DDB1. A duplication within MLL exon 3 was detected in another patient with 11qUPD. We identified several common deleted regions (CDR) on chromosome 11. One of the CDRs associated with de novo acute myeloid leukemia (P=0.013). One patient with a deletion at the LMO2 locus harbored an additional point mutation on the other allele indicating that LMO2 might be a tumor suppressor frequently targeted by 11p deletions. Our chromosome-centered analysis indicates that chromosome 11 contains a number of tumor suppressor genes and that the role of this chromosome in myeloid malignancies is more complex than previously recognized.     

Mutation pattern of paired immunoglobulin heavy and light variable domains in chronic lymphocytic leukemia B cells

B-cell chronic lymphocytic leukemia (CLL) patients display leukemic clones bearing either germline or somatically mutated immunoglobulin heavy variable (IGHV ) genes. Most information on CLL immunoglobulins (Igs), such as the definition of stereotyped B-cell receptors (BCRs), was derived from germline unmutated Igs. In particular, detailed studies on the distribution and nature of mutations in paired heavy- and light-chain domains of CLL clones bearing mutated Igs are lacking. To address the somatic hyper-mutation dynamics of CLL Igs, we analyzed the mutation pattern of paired IGHV-diversity-joining (IGHV-D-J ) and immunoglobulin kappa/lambda variable-joining (IGK/LV-J ) rearrangements of 193 leukemic clones that displayed ≥ 2% mutations in at least one of the two immunoglobulin variable (IGV ) genes (IGHV and/or IGK/LV ). The relationship between the mutation frequency in IGHV and IGK/LV complementarity determining regions (CDRs) and framework regions (FRs) was evaluated by correlation analysis. Replacement (R) mutation frequency within IGK/LV chain CDRs correlated significantly with mutation frequency of paired IGHV CDRs in λ but not κ isotype CLL clones. CDRs of IGKV-J rearrangements displayed a lower percentage of R mutations than IGHVs. The frequency/pattern of mutations in kappa CLL Igs differed also from that in κ-expressing normal B cells described in the literature. Instead, the mutation frequency within the FRs of IGHV and either IGKV or IGLV was correlated. Notably, the amount of diversity introduced by replaced amino acids was comparable between IGHVs and IGKVs. The data indicate a different mutation pattern between κ and λ isotype CLL clones and suggest an antigenic selection that, in κ samples, operates against CDR variation.     

The challenges of modelling antibody repertoire dynamics in HIV infection

Antibody affinity maturation by somatic hypermutation of B-cell immunoglobulin variable region genes has been studied for decades in various model systems using well-defined antigens. While much is known about the molecular details of the process, our understanding of the selective forces that generate affinity maturation are less well developed, particularly in the case of a co-evolving pathogen such as HIV. Despite this gap in understanding, high-throughput antibody sequence data are increasingly being collected to investigate the evolutionary trajectories of antibody lineages in HIV-infected individuals. Here, we review what is known in controlled experimental systems about the mechanisms underlying antibody selection and compare this to the observed temporal patterns of antibody evolution in HIV infection. We describe how our current understanding of antibody selection mechanisms leaves questions about antibody dynamics in HIV infection unanswered. Without a mechanistic understanding of antibody selection in the context of a co-evolving viral population, modelling and analysis of antibody sequences in HIV-infected individuals will be limited in their interpretation and predictive ability.     

Trade-offs in antibody repertoires to complex antigens

Pathogens vary in their antigenic complexity. While some pathogens such as measles present a few relatively invariant targets to the immune system, others such as malaria display considerable antigenic diversity. How the immune response copes in the presence of multiple antigens, and whether a trade-off exists between the breadth and efficacy of antibody (Ab)-mediated immune responses, are unsolved problems. We present a theoretical model of affinity maturation of B-cell receptors (BCRs) during a primary infection and examine how variation in the number of accessible antigenic sites alters the Ab repertoire. Naive B cells with randomly generated receptor sequences initiate the germinal centre (GC) reaction. The binding affinity of a BCR to an antigen is quantified via a genotype–phenotype map, based on a random energy landscape, that combines local and distant interactions between residues. In the presence of numerous antigens or epitopes, B-cell clones with different specificities compete for stimulation during rounds of mutation within GCs. We find that the availability of many epitopes reduces the affinity and relative breadth of the Ab repertoire. Despite the stochasticity of somatic hypermutation, patterns of immunodominance are strongly shaped by chance selection of naive B cells with specificities for particular epitopes. Our model provides a mechanistic basis for the diversity of Ab repertoires and the evolutionary advantage of antigenically complex pathogens.     

Burkitt's lymphoma is a malignancy of mature B cells expressing somatically mutated V region genes.

BACKGROUND: The developmental stage from which stems the malignant B cell population in Burkitt's lymphoma (BL) is unclear. An approach to answering this question is provided by the sequence analysis of rear-ranged immunoglobulin (Ig) variable region (V) genes from BL for evidence of somatic mutations, together with a phenotypic characterization. As somatic hypermutation of Ig V region genes occurs in germinal center B cells, somatically mutated Ig genes are found in germinal center B cells and their descendents. MATERIALS AND METHODS: Rearranged V kappa region genes from 10 kappa-expressing sporadic and endemic BL-derived cell lines (9 IgM and 1 IgG positive) and three kappa-expressing endemic BL biopsy specimens were amplified by polymerase chain reaction and sequenced. In addition, VH region gene sequences from these cell lines were determined. RESULTS: All BL cell lines and the three biopsy specimens carried somatically mutated V region genes. The average mutation frequency of rearranged V kappa genes from eight BL cell lines established from sporadic BL was 1.8%. A higher frequency (6%) was found in five endemic cases (three biopsy specimens and two BL cell lines). CONCLUSIONS: The detection of somatic mutations in the rearranged V region genes suggests that both sporadic and endemic BL represent a B-cell malignancy originating from germinal center B cells or their descendants. Interestingly, the mutation frequency detected in sporadic BL is in a range similar to that characteristic for IgM-expressing B cells in the human peripheral blood and for mu chain-expressing germinal center B cells, whereas the mutation frequency found in endemic BL is significantly higher.