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.     

Somatic hypermutation of the JC intron is markedly reduced in unrearranged kappa and H alleles and is unevenly distributed in rearranged alleles

Somatic hypermutation of the Ig genes occurs in rearranged V(D)J and its flanking sequences after Ag stimulation. Even though C regions and unrearranged V segments have been found to lack mutations, it is not known whether the mutational mechanism can be active in unrearranged J segments and their flanking regions. By polymerase chain reaction and direct sequencing of the 500 bp at the 5' end of the JC intron of hybridoma DNA derived from splenic B cells, we show that the frequency of mutations in unrearranged J regions of kappa and H chain genes is 0/7849 bp (upper 95% confidence interval, less than 0.00038) and 1/3209 bp (upper 95% confidence interval, less than 0.0015), respectively. The frequency (f) for the same region of rearranged kappa and H chain genes was 29/9380 bp (95% confidence, 0.0021 less than f less than 0.0044) and 16/2750 bp (95% confidence, 0.0033 less than f less than 0.0094), respectively. The significantly higher frequency of mutations in the rearranged alleles indicates that rearrangement is needed to effect full activation of the mutational mechanism. The data also show that mutations occur predominantly in the 5'-most 250 bp of the JC kappa intron. Statistical analysis of the distribution of mutations within the 5'-most 521 bp of the JC kappa intron reveals significant deviation (p = 0.000085) from a theoretically determined uniform distribution, indicating that mutations are not evenly distributed within this region.     

Characterization of higher avidity monoclonal antibodies produced by murine B-cell hybridoma variants selected for increased antigen binding of membrane Ig

Somatic mutation in the Ig genes plays a major role in the increase of antibody affinity observed in secondary immunologic responses. It has been shown that the mechanism responsible for the high rate of somatic mutation in the Ig genes was active not only in normal B lymphocytes but also in B-cell hybridomas secreting mAb. Also, it has been reported that B-cell hybridomas were positive for membrane Ig of the same specificity as the secreted mAb. The presence of membrane Ig suggested that somatic variants secreting mAb of higher affinity could be selected by the increased capacity of these hybridoma cells to bind immobilized Ag. This hypothesis was tested with hybridoma cells secreting an IgM mAb reacting with the A Ag of the ABO blood group system. In two selection experiments, we have isolated several variant cell lines secreting mAb of increased avidity for the A Ag under similar IgM concentrations. Biochemical characterization of one of the variant mAb indicated that the mutation responsible for the increased avidity has occurred in the heavy chain gene. The method developed may have profound implications for the diagnostic and therapeutic use of mAb and will permit the study, in an in vitro system, of the role of somatic mutations in antibody diversity.     

Fitness effects of derived deleterious mutations in four closely related wild tomato species with spatial structure

A key issue in evolutionary biology is an improved understanding of the genetic mechanisms by which species adapt to various environments. Using DNA sequence data, it is possible to quantify the number of adaptive and deleterious mutations, and the distribution of fitness effects of new mutations (its mean and variance) by simultaneously taking into account the demography of a given species. We investigated how selection functions at eight housekeeping genes of four closely related, outcrossing species of wild tomatoes that are native to diverse environments in western South America (Solanum arcanum, S. chilense, S. habrochaites and S. peruvianum). We found little evidence for adaptive mutations but pervasive evidence for strong purifying selection in coding regions of the four species. In contrast, the strength of purifying selection seems to vary among the four species in non-coding (NC) regions (introns). Using F(ST)-based measures of fixation in subdivided populations, we suggest that weak purifying selection has affected the NC regions of S. habrochaites, S. chilense and S. peruvianum. In contrast, NC regions in S. arcanum show a distribution of fitness effects with mutations being either nearly neutral or very strongly deleterious. These results suggest that closely related species with similar genetic backgrounds but experiencing contrasting environments differ in the variance of deleterious fitness effects.     

The log multinomial regression model for nominal outcomes with more than two attributes

An estimate of the risk or prevalence ratio, adjusted for confounders, can be obtained from a log binomial model (binomial errors, log link) fitted to binary outcome data. We propose a modification of the log binomial model to obtain relative risk estimates for nominal outcomes with more than two attributes (the "log multinomial model"). Extensive data simulations were undertaken to compare the performance of the log multinomial model with that of an expanded data multinomial logistic regression method based on the approach proposed by Schouten et al. (1993) for binary data, and with that of separate fits of a Poisson regression model based on the approach proposed by Zou (2004) and Carter, Lipsitz and Tilley (2005) for binary data. Log multinomial regression resulted in "inadmissable" solutions (out-of-bounds probabilities) exceeding 50% in some data settings. Coefficient estimates by the alternative methods produced out-of-bounds probabilities for the log multinomial model in up to 27% of samples to which a log multinomial model had been successfully fitted. The log multinomial coefficient estimates generally had lesser relative bias and mean squared error than the alternative methods. The practical utility of the log multinomial regression model was demonstrated with a real data example. The log multinomial model offers a practical solution to the problem of obtaining adjusted estimates of the risk ratio in the multinomial setting, but must be used with some care and attention to detail.     

Evolution of Ig DNA sequence to target specific base positions within codons for somatic hypermutation

Ig variable (V) region genes are subjected to a somatic hypermutation process as B lymphocytes participate in immune reactions to protein Ags. Although little is known regarding the mechanism of mutagenesis, a consistent hierarchy of trinucleotide target preferences is evident. Analysis of trinucleotide regional distributions predicted and we now empirically confirm the surprising finding that the framework 2 region of kappa V region genes is highly mutable despite its importance to the structural integrity and function of the Ab molecule. Interestingly, much of this mutability appears to be focused on the third codon position where synonymous substitutions are most likely to occur. We also observed a trend for high predicted mutability for codon positions 1 and 2 in complementarity-determining regions. Consequently, amino acid replacements should occur at a higher rate in complementarity-determining regions than in framework regions due to the distribution and subsequent targeting of microsequences by the mutation mechanism. Our results reveal a subtle tier of V region gene evolution in which DNA sequence has been molded to direct mutations to specific base positions within codons in a manner that minimizes damage and maximizes the benefits of the somatic hypermutation process.     

Background selection in single genes may explain patterns of codon bias

Background selection involves the reduction in effective population size caused by the removal of recurrent deleterious mutations from a population. Previous work has examined this process for large genomic regions. Here we focus on the level of a single gene or small group of genes and investigate how the effects of background selection caused by nonsynonymous mutations are influenced by the lengths of coding sequences, the number and length of introns, intergenic distances, neighboring genes, mutation rate, and recombination rate. We generate our predictions from estimates of the distribution of the fitness effects of nonsynonymous mutations, obtained from DNA sequence diversity data in Drosophila. Results for genes in regions with typical frequencies of crossing over in Drosophila melanogaster suggest that background selection may influence the effective population sizes of different regions of the same gene, consistent with observed differences in codon usage bias along genes. It may also help to cause the observed effects of gene length and introns on codon usage. Gene conversion plays a crucial role in determining the sizes of these effects. The model overpredicts the effects of background selection with large groups of nonrecombining genes, because it ignores Hill-Robertson interference among the mutations involved.     

Quasi-equilibrium theory for the distribution of rare alleles in a subdivided population: justification and implications

This paper examines a quasi-equilibrium theory of rare alleles for subdivided populations that follow an island-model version of the Wright-Fisher model of evolution. All mutations are assumed to create new alleles. We present four results: (1) conditions for the theory to apply are formally established using properties of the moments of the binomial distribution; (2) approximations currently in the literature can be replaced with exact results that are in better agreement with our simulations; (3) a modified maximum likelihood estimator of migration rate exhibits the same good performance on island-model data or on data simulated from the multinomial mixed with the Dirichlet distribution, and (4) a connection between the rare-allele method and the Ewens Sampling Formula for the infinite-allele mutation model is made. This introduces a new and simpler proof for the expected number of alleles implied by the Ewens Sampling Formula.     

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.     

Characterization of Ig gene somatic hypermutation in the absence of activation-induced cytidine deaminase

Somatic hypermutation (SHM) of Ig genes depends upon the deamination of C nucleotides in WRCY (W = A/T, R = A/G, Y = C/T) motifs by activation-induced cytidine deaminase (AICDA). Despite this, a large number of mutations occur in WA motifs that can be accounted for by the activity of polymerase eta (POL eta). To determine whether there are AICDA-independent mutations and to characterize the relationship between AICDA- and POL eta-mediated mutations, 1470 H chain and 1313 kappa- and lambda-chain rearrangements from three AICDA(-/-) patients were analyzed. The Ig mutation frequency of all V(H) genes from AICDA(-/-) patients was 40-fold less than that of normal donors, whereas the mutation frequency of mutated V(H) sequences from AICDA(-/-) patients was 6.8-fold less than that of normal donors. AICDA(-/-) B cells lack mutations in WRCY/RGYW motifs as well as replacement mutations and mutational targeting in complementarity-determining regions. A significantly reduced mutation frequency in WA motifs compared with normal donors and an increased percentage of transitions, which may relate to reduced uracil DNA-glycosylase activity, suggest a role for AICDA in regulating POL eta and uracil DNA-glycosylase activity. Similar results were observed in V(L) rearrangements. The residual mutations were predominantly G:C substitutions, indicating that AICDA-independent cytidine deamination was a likely, yet inefficient, mechanism for mutating Ig genes.     

Adaptive evolution of variable region genes encoding an unusual type of immunoglobulin in camelids

A typical immunoglobulin (Ig) molecule is composed of four polypeptide chains: two identical heavy (H) chains and two identical light (L) chains. This tetrameric structure is conserved in almost all jawed vertebrate species. However, it has been discovered that camels and llamas (family: Camelidae) possess a type of dimeric Ig that consists of two H chains only. These H chains do not associate with L chains, and they do not have the first constant region (CH1), which is present in the conventional Ig. In spite of these changes, the dimeric Ig maintains the normal immune function. To understand the evolution of the dimeric Ig, we studied the phylogenetic relationships of the variable region (V(H)H) genes of the dimeric Ig from Camelidae and those (V(H)) of the conventional Ig from mammals. The results showed that the V(H)H genes form a monophyletic cluster within one of the mammalian V(H) groups, group C. We examined the type of selective force in complementarity-determining regions (CDRs) and framework regions (FRs) by comparing the rate of synonymous (dS) and nonsynonymous (dN) substitutions. We found that the results obtained from V(H)H genes were similar to those from V(H) genes in that CDRs showed an excess of dN over dS (indicating positive selection), whereas the reverse was true for FRs (purifying selection). However, when the extent of positive selection or purifying selection was investigated at each codon site, three major differences between V(H)H and V(H) genes were found. That is, very different types of selective force were observed between V(H)H and V(H) genes (1) at the sites that contact the L chain in the conventional Ig, (2) at the sites that interact with the CH1 region in the conventional Ig, and (3) in the H1 loop. Our findings suggest that adaptive evolution has occurred in the functionally important sites of the V(H)H genes to maintain the normal immune function in the dimeric Ig.     

The Question of the Total Gene Number in DROSOPHILA MELANOGASTER

A statistical analysis has been carried out on the distribution and allelism of nearly 500 sex-linked, X-ray-induced, cytologically normal and rearranged lethal mutations in Drosophila melanogaster that were obtained by G. Lefevre. The mutations were induced in four different regions of the X chromosome: (1) 1A1-3E8, (2) 6D1-8A5, (3) 9E1-11A7 and (4) 19A1-20F4, which together comprise more than one-third of the entire chromosome.--The analysis shows that the number of alleles found at different loci does not fit a Poisson distribution, even when the proper procedures are taken to accommodate the truncated nature of the data. However, the allele distribution fits a truncated negative binomial distribution quite well, with cytologically normal mutations fitting better than rearrangement mutations. This indicates that genes are not equimutable, as required for the data to fit a Poisson distribution.--Using the negative binomial parameters to estimate the number of genes that did not produce a detectable lethal mutation in our experiment (n0) gave a larger number than that derived from the use of the Poisson parameter. Unfortunately, we cannot estimate the total numbers of nonvital loci, loci with undetectable phenotypes and loci having extremely low mutabilities. In any event, our estimate of the total vital gene number was far short of the total number of bands in the analyzed regions; yet, in several short intervals, we have found more vital genes than bands; in other intervals, fewer. We conclude that the one-band, one-gene hypothesis, in its literal sense, is not true; furthermore, it is difficult to support, even approximately.--The question of the total gene number in Drosophila will, not doubt, eventually be solved by molecular analyses, not by statistical analysis of mutation data or saturation studies.     

A human systemic lupus erythematosus-related anti-cardiolipin/single-stranded DNA autoantibody is encoded by a somatically mutated variant of the developmentally restricted 51P1 VH gene.

We report the Ig H and L chain V region sequences from the cDNAs encoding a monoclonal human IgG anti-cardiolipin/ssDNA autoantibody (R149) derived from a patient with active SLE. Comparison with the germ-line V-gene repertoire of this patient revealed that R149 likely arose as a consequence of an Ag-driven selection process. The Ag-binding portions of the V regions were characterized by a high number of arginine residues, a property that has been associated with anti-dsDNA autoantibodies from lupus-prone mice and patients with SLE. The VH gene encoding autoantibody R149 was a somatically mutated variant of the 51P1 gene segment, which is frequently associated with the restricted fetal B cell repertoire, malignant CD5 B cells, and natural autoantibodies. These data suggest that in SLE patients a common antigenic stimulus may evoke anti-DNA and anti-cardiolipin autoantibodies and provide further evidence that a small set of developmentally restricted VH genes can give rise to disease-associated autoantibodies through Ag-selected somatic mutations.     

A simple binomial test for estimating sequencing errors in public repository 16S rRNA sequences

Sequences in public databases may contain a number of sequencing errors. A double binomial model describing the distribution of indel-excluded similarity coefficients (S) among repeatedly sequenced 16S rRNA was previously developed and it produced a confidence interval of S useful for testing sequence identity among sequences of 400-bp length. We characterized patterns in sequencing errors found in nearly complete 16S rRNA sequences of Vibrionaceae as highly variable in reported sequence length and containing a small number of indels. To accommodate these characteristics, a simple binomial model for distribution of the similarity coefficient (H) that included indels was derived from the double binomial model for S. The model showed good fit to empirical data. By using either a pre-determined or bootstrapping estimated standard probability of base matching, we were able to use the exact binomial test to determine the relative level of sequencing error for a given pair of duplicated sequences. A limitation of the method is the requirement that duplicated sequences for the same template sequence be paired, but this can be overcome by using only conserved regions of 16S rRNA sequences and pairing a given sequence with its highest scoring BLAST search hit from the nr database of GenBank.     

Inference of site frequency spectra from high-throughput sequence data: quantification of selection on nonsynonymous and synonymous sites in humans

Sequencing errors and random sampling of nucleotide types among sequencing reads at heterozygous sites present challenges for accurate, unbiased inference of single-nucleotide polymorphism genotypes from high-throughput sequence data. Here, we develop a maximum-likelihood approach to estimate the frequency distribution of the number of alleles in a sample of individuals (the site frequency spectrum), using high-throughput sequence data. Our method assumes binomial sampling of nucleotide types in heterozygotes and random sequencing error. By simulations, we show that close to unbiased estimates of the site frequency spectrum can be obtained if the error rate per base read does not exceed the population nucleotide diversity. We also show that these estimates are reasonably robust if errors are nonrandom. We then apply the method to infer site frequency spectra for zerofold degenerate, fourfold degenerate, and intronic sites of protein-coding genes using the low coverage human sequence data produced by the 1000 Genomes Project phase-one pilot. By fitting a model to the inferred site frequency spectra that estimates parameters of the distribution of fitness effects of new mutations, we find evidence for significant natural selection operating on fourfold sites. We also find that a model with variable effects of mutations at synonymous sites fits the data significantly better than a model with equal mutational effects. Under the variable effects model, we infer that 11% of synonymous mutations are subject to strong purifying selection.     

Activation-induced cytidine deaminase induces reproducible DNA breaks at many non-Ig Loci in activated B cells

After immunization or infection, activation-induced cytidine deaminase (AID) initiates diversification of immunoglobulin (Ig) genes in B cells, introducing mutations within the antigen-binding V regions (somatic hypermutation, SHM) and double-strand DNA breaks (DSBs) into switch (S) regions, leading to antibody class switch recombination (CSR). We asked if, during B cell activation, AID also induces DNA breaks at genes other than IgH genes. Using a nonbiased genome-wide approach, we have identified hundreds of reproducible, AID-dependent DSBs in mouse splenic B cells shortly after induction of CSR in culture. Most interestingly, AID induces DSBs at sites syntenic with sites of translocations, deletions, and amplifications found in human B cell lymphomas, including within the oncogene B cell lymphoma11a (bcl11a)/evi9. Unlike AID-induced DSBs in Ig genes, genome-wide AID-dependent DSBs are not restricted to transcribed regions and frequently occur within repeated sequence elements, including CA repeats, non-CA tandem repeats, and SINEs.     

The contested role of uracil DNA glycosylase in immunoglobulin gene diversification

Class switch recombination (CSR) and somatic hypermutation (SHM) of immunoglobulin (Ig) genes are initiated by the activation-induced cytosine deaminase AID. The resulting uracils in Ig genes were believed to be removed by the uracil glycosylase (UNG) and the resulting abasic sites treated in an error-prone fashion, creating breaks in the Ig switch regions and mutations in the variable regions. A recent report suggests that UNG does not act as a glycosylase in CSR and SHM but rather has unknown activity subsequent to DNA breaks that were created by other mechanisms.     

A composite-likelihood approach for detecting directional selection from DNA sequence data

We present a novel composite-likelihood-ratio test (CLRT) for detecting genes and genomic regions that are subject to recurrent natural selection (either positive or negative). The method uses the likelihood functions of Hartl et al. (1994) for inference in a Wright-Fisher genic selection model and corrects for nonindependence among sites by application of coalescent simulations with recombination. Here, we (1) characterize the distribution of the CLRT statistic (Lambda) as a function of the population recombination rate (R=4Ner); (2) explore the effects of bias in estimation of R on the size (type I error) of the CLRT; (3) explore the robustness of the model to population growth, bottlenecks, and migration; (4) explore the power of the CLRT under varying levels of mutation, selection, and recombination; (5) explore the discriminatory power of the test in distinguishing negative selection from population growth; and (6) evaluate the performance of maximum composite-likelihood estimation (MCLE) of the selection coefficient. We find that the test has excellent power to detect weak negative selection and moderate power to detect positive selection. Moreover, the test is quite robust to bias in the estimate of local recombination rate, but not to certain demographic scenarios such as population growth or a recent bottleneck. Last, we demonstrate that the MCLE of the selection parameter has little bias for weak negative selection and has downward bias for positively selected mutations.     

SAIVGeM: spreadsheet analysis of immunoglobulin VH gene mutations

The analysis of mutations in immunoglobulin heavy chain variable (IGHV) region genes is a tedious process when performed by hand on multiple sequences. This report describes a set of linked Microsoft Excel files that perform several common analyses on large numbers of IGHV sequences. The spreadsheet analysis of immunoglobulin VH gene mutations (SAIVGeM) package determines the distribution of mutations among each nucleotide, the nature of the mutation at both the nucleotide and amino acid level, the frequency of mutation in the A/G G C/T A/T (RGYW) hotspot motifs of both strand polarity, and the distribution of replacement and silent mutations among the complementarity determining regions (CDRs) and the framework regions (FRs) of the immunoglobulin gene as defined by either the Kabat or IMGT conventions. These parameters are summarized and graphically presented where appropriate. In addition, the SAIVGeM package analyzes those mutations that occur in third positions of redundant codons. Because any nucleotide change in these positions is inherently silent, these positions can be used to study the mutational spectra without biases from the selection of protein structure.