PipeIT2: A tumor-only somatic variant calling workflow for molecular diagnostic Ion Torrent sequencing data

Precision oncology relies on the accurate identification of somatic mutations in cancer patients. While the sequencing of the tumoral tissue is frequently part of routine clinical care, the healthy counterparts are rarely sequenced. We previously published PipeIT, a somatic variant calling workflow specific for Ion Torrent sequencing data enclosed in a Singularity container. PipeIT combines user-friendly execution, reproducibility and reliable mutation identification, but relies on matched germline sequencing data to exclude germline variants. Expanding on the original PipeIT, here we describe PipeIT2 to address the clinical need to define somatic mutations in the absence of germline control. We show that PipeIT2 achieves a > 95% recall for variants with variant allele fraction >10%, reliably detects driver and actionable mutations and filters out most of the germline mutations and sequencing artifacts. With its performance, reproducibility, and ease of execution, PipeIT2 is a valuable addition to molecular diagnostics laboratories.     

Mutation Rate Variability across Human Y-Chromosome Haplogroups

A common assumption in dating patrilineal events using Y-chromosome sequencing data is that the Y-chromosome mutation rate is invariant across haplogroups. Previous studies revealed interhaplogroup heterogeneity in phylogenetic branch length. Whether this heterogeneity is caused by interhaplogroup mutation rate variation or nongenetic confounders remains unknown. Here, we analyzed whole-genome sequences from cultured cells derived from >1,700 males. We confirmed the presence of branch length heterogeneity. We demonstrate that sex-chromosome mutations that appear within cell lines, which likely occurred somatically or in vitro (and are thus not influenced by nongenetic confounders) are informative for germline mutational processes. Using within-cell-line mutations, we computed a relative Y-chromosome somatic mutation rate, and uncovered substantial variation (up to 83.3%) in this proxy for germline mutation rate among haplogroups. This rate positively correlates with phylogenetic branch length, indicating that interhaplogroup mutation rate variation is a likely cause of branch length heterogeneity.     

High-Throughput Genome Editing and Phenotyping Facilitated by High Resolution Melting Curve Analysis

With the goal to generate and characterize the phenotypes of null alleles in all genes within an organism and the recent advances in custom nucleases, genome editing limitations have moved from mutation generation to mutation detection. We previously demonstrated that High Resolution Melting (HRM) analysis is a rapid and efficient means of genotyping known zebrafish mutants. Here we establish optimized conditions for HRM based detection of novel mutant alleles. Using these conditions, we demonstrate that HRM is highly efficient at mutation detection across multiple genome editing platforms (ZFNs, TALENs, and CRISPRs); we observed nuclease generated HRM positive targeting in 1 of 6 (16%) open pool derived ZFNs, 14 of 23 (60%) TALENs, and 58 of 77 (75%) CRISPR nucleases. Successful targeting, based on HRM of G0 embryos correlates well with successful germline transmission (46 of 47 nucleases); yet, surprisingly mutations in the somatic tail DNA weakly correlate with mutations in the germline F1 progeny DNA. This suggests that analysis of G0 tail DNA is a good indicator of the efficiency of the nuclease, but not necessarily a good indicator of germline alleles that will be present in the F1s. However, we demonstrate that small amplicon HRM curve profiles of F1 progeny DNA can be used to differentiate between specific mutant alleles, facilitating rare allele identification and isolation; and that HRM is a powerful technique for screening possible off-target mutations that may be generated by the nucleases. Our data suggest that micro-homology based alternative NHEJ repair is primarily utilized in the generation of CRISPR mutant alleles and allows us to predict likelihood of generating a null allele. Lastly, we demonstrate that HRM can be used to quickly distinguish genotype-phenotype correlations within F1 embryos derived from G0 intercrosses. Together these data indicate that custom nucleases, in conjunction with the ease and speed of HRM, will facilitate future high-throughput mutation generation and analysis needed to establish mutants in all genes of an organism.     

A method to reduce ancestry related germline false positives in tumor only somatic variant calling

Background

Significant clinical and research applications are driving large scale adoption of individualized tumor sequencing in cancer in order to identify tumors-specific mutations. When a matched germline sample is available, somatic mutations may be identified using comparative callers. However, matched germline samples are frequently not available such as with archival tissues, which makes it difficult to distinguish somatic from germline variants. While population databases may be used to filter out known germline variants, recent studies have shown private germline variants result in an inflated false positive rate in unmatched tumor samples, and the number germline false positives in an individual may be related to ancestry.

Methods

First, we examined the relationship between the germline false positives and ancestry. Then we developed and implemented a tumor only caller (LumosVar) that leverages differences in allelic frequency between somatic and germline variants in impure tumors. We used simulated data to systematically examine how copy number alterations, tumor purity, and sequencing depth should affect the sensitivity of our caller. Finally, we evaluated the caller on real data.

Results

We find the germline false-positive rate is significantly higher for individuals of non-European Ancestry largely due to the limited diversity in public polymorphism databases and due to population-specific characteristics such as admixture or recent expansions. Our Bayesian tumor only caller (LumosVar) is able to greatly reduce false positives from private germline variants, and our sensitivity is similar to predictions based on simulated data.

Conclusions

Taken together, our results suggest that studies of individuals of non-European ancestry would most benefit from our approach. However, high sensitivity requires sufficiently impure tumors and adequate sequencing depth. Even in impure tumors, there are copy number alterations that result in germline and somatic variants having similar allele frequencies, limiting the sensitivity of the approach. We believe our approach could greatly improve the analysis of archival samples in a research setting where the normal is not available.

     

multiSNV: a probabilistic approach for improving detection of somatic point mutations from multiple related tumour samples

Somatic variant analysis of a tumour sample and its matched normal has been widely used in cancer research to distinguish germline polymorphisms from somatic mutations. However, due to the extensive intratumour heterogeneity of cancer, sequencing data from a single tumour sample may greatly underestimate the overall mutational landscape. In recent studies, multiple spatially or temporally separated tumour samples from the same patient were sequenced to identify the regional distribution of somatic mutations and study intratumour heterogeneity. There are a number of tools to perform somatic variant calling from matched tumour-normal next-generation sequencing (NGS) data; however none of these allow joint analysis of multiple same-patient samples. We discuss the benefits and challenges of multisample somatic variant calling and present multiSNV, a software package for calling single nucleotide variants (SNVs) using NGS data from multiple same-patient samples. Instead of performing multiple pairwise analyses of a single tumour sample and a matched normal, multiSNV jointly considers all available samples under a Bayesian framework to increase sensitivity of calling shared SNVs. By leveraging information from all available samples, multiSNV is able to detect rare mutations with variant allele frequencies down to 3% from whole-exome sequencing experiments.     

Mammalian Male Mutation Bias: Impacts of Generation Time and Regional Variation in Substitution Rates

In mammals, males undergo a greater number of germline cell divisions compared with females. Thus, the male germline accumulates more DNA replication errors, which result in male mutation bias—a higher mutation rate for males than for females. The phenomenon of male mutation bias has been investigated mostly for rodents and primates, however, it has not been studied in detail for other mammalian orders. Here we sequenced and analyzed five introns of three genes (DBX/DBY, UTX/UTY, and ZFX/ZFY) homologous between X and Y chromosomes in several species of perissodactyls (horses and rhinos) and of primates. Male mutation bias was evident: substitution rate was higher for a Y chromosome intron than for its X chromosome homologue for all five intron pairs studied. Substitution rates varied regionally among introns sequenced on the same chromosome and this variation influenced male mutation bias inferred from each intron pair. Interestingly, we observed a positive correlation in substitution rates between homologous X and homologous Y introns as well as between orthologous primate and perissodactyl introns. The male-to-female mutation rate ratio estimated from concatenated sequences of five perissodactyl introns was 3.88 (95% CI = 2.90–6.07). Using the data generated here and estimates available in the literature, we compared male mutation bias among several mammalian orders. We conclude that male mutation bias is significantly higher for organisms with long generation times (primates, perissodactyls, and felids) than for organisms with short generation times (e.g., rodents) since the former undergo a greater number of male germline cell divisions. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Dr. Deborah Charlesworth]     

UV-exposure,endogenous DNA damage,and DNA replication errors shape the spectra of genome changes in human skin

Human skin is continuously exposed to environmental DNA damage leading to the accumulation of somatic mutations over the lifetime of an individual. Mutagenesis in human skin cells can be also caused by endogenous DNA damage and by DNA replication errors. The contributions of these processes to the somatic mutation load in the skin of healthy humans has so far not been accurately assessed because the low numbers of mutations from current sequencing methodologies preclude the distinction between sequencing errors and true somatic genome changes. In this work, we sequenced genomes of single cell-derived clonal lineages obtained from primary skin cells of a large cohort of healthy individuals across a wide range of ages. We report here the range of mutation load and a comprehensive view of the various somatic genome changes that accumulate in skin cells. We demonstrate that UV-induced base substitutions, insertions and deletions are prominent even in sun-shielded skin. In addition, we detect accumulation of mutations due to spontaneous deamination of methylated cytosines as well as insertions and deletions characteristic of DNA replication errors in these cells. The endogenously induced somatic mutations and indels also demonstrate a linear increase with age, while UV-induced mutation load is age-independent. Finally, we show that DNA replication stalling at common fragile sites are potent sources of gross chromosomal rearrangements in human cells. Thus, somatic mutations in skin of healthy individuals reflect the interplay of environmental and endogenous factors in facilitating genome instability and carcinogenesis.     

Stem cells are units of natural selection in a colonial ascidian

Stem cells are highly conserved biological units of development and regeneration. Here we formally demonstrate that stem cell lineages are also legitimate units of natural selection. In a colonial ascidian, Botryllus schlosseri, vascular fusion between genetically distinct individuals results in cellular parasitism of somatic tissues, gametes, or both. We show that genetic hierarchies of somatic and gametic parasitism following fusion can be replicated by transplanting cells between colonies. We prospectively isolate a population of multipotent, self-renewing stem cells that retain their competitive phenotype upon transplantation. Their single-cell contribution to either somatic or germline fates, but not to both, is consistent with separate lineages of somatic and germline stem cells or pluripotent stem cells that differentiate according to the niche in which they land. Since fusion is restricted to individuals that share a fusion/histocompatibility allele, these data suggest that histocompatibility genes in Botryllus evolved to protect the body from parasitic stem cells usurping asexual or sexual inheritance.     

Germline and somatic mutation profile in Cancer patients revealed by a medium-sized pan-Cancer panel

Gene mutation detection and the resulted precision-medicine therapy is transforming clinical practice. Here, we report the use of a custom-developed, medium-sized, pan-cancer probe panel for the detection of somatic and germline mutations. We used a hybridization capture-based NGS assay for targeted deep sequencing of all exons and selected introns of 181 key cancer driver genes, covering both inherited risks and somatic mutations. We performed paired-variant calling on tumor samples and their matched normal samples. We processed clinical patient samples of formalin-fixed, paraffin embedded tumors (FFPE samples) and cell-free peripheral blood (cfDNA samples). We found germline mutations of inherited cancer risk at 9%; and discovered a novel germline mutation in BRCA1. Somatic mutation rate in driver genes is at 73.1%, much higher than previously reported. On recommending precision-medicine therapeutics, we achieved 91.6% for patients with FFPE samples.     

Mosaike bei monogenen Erkrankungen

Mosaicism is defined as the simultaneous presence of cells with different genotypes that originate from a common zygote. Mutations can either be present in germline or somatic cells. Monogenic disorders apparently caused by a de novo mutation may show a recurrence risk due to germline mosaicism in a parent. Duchenne muscular dystrophy is a well investigated example with a high frequency of germline mosaicism and the estimation for the risk of recurrence is based on theoretical models and empirical data. Recently, somatic mutations have been uncovered in various syndromic disorders, such as Proteus syndrome or hemimegalencephaly and respective mutations often show gain-of-function properties. Genetic testing is mainly based on next generation sequencing technologies but still remains challenging; however, detection of somatic mosaicism is expected to be of increasing relevance in the diagnosis of monogenic disorders. Somatic mosaicism may also play a hitherto underestimated role in common disorders.     

Inhibition of human cytochrome P450 1B1, 1A1 and 1A2 by antigenotoxic compounds,purpurin and alizarin

Expanded simple tandem repeat (ESTR) loci include some of the most unstable DNA in the mouse genome and have been extensively used in pedigree studies of germline mutation. We now show that repeat DNA instability at the mouse ESTR locus Ms6-hm can also be monitored by single molecule PCR analysis of genomic DNA. Unlike unstable human minisatellites which mutate almost exclusively in the germline by a meiotic recombination-based process, mouse Ms6-hm shows repeat instability both in germinal (sperm) DNA and in somatic (spleen, brain) DNA. There is no significant variation in mutation frequency between mice of the same inbred strain. However, significant variation occurs between tissues, with mice showing the highest mutation frequency in sperm. The size spectra of somatic and sperm mutants are indistinguishable and heavily biased towards gains and losses of only a few repeat units, suggesting repeat turnover by a mitotic replication slippage process operating both in the soma and in the germline. Analysis of male mice following acute pre-meiotic exposure to X-rays showed a significant increase in sperm but not somatic mutation frequency, though no change in the size spectrum of mutants. The level of radiation-induced mutation at Ms6-hm was indistinguishable from that established by conventional pedigree analysis following paternal irradiation. This confirms that mouse ESTR loci are very sensitive to ionizing radiation and establishes that induced germline mutation results from radiation-induced mutant alleles being present in sperm, rather than from unrepaired sperm DNA lesions that subsequently lead to the appearance of mutants in the early embryo. This single molecule monitoring system has the potential to substantially reduce the number of mice needed for germline mutation monitoring, and can be used to study not only germline mutation but also somatic mutation in vivo and in cell culture.     

Evaluating somatic tumor mutation detection without matched normal samples

Background

Observations of recurrent somatic mutations in tumors have led to identification and definition of signaling and other pathways that are important for cancer progression and therapeutic targeting. As tumor cells contain both an individual’s inherited genetic variants and somatic mutations, challenges arise in distinguishing these events in massively parallel sequencing datasets. Typically, both a tumor sample and a “normal” sample from the same individual are sequenced and compared; variants observed only in the tumor are considered to be somatic mutations. However, this approach requires two samples for each individual.

Results

We evaluate a method of detecting somatic mutations in tumor samples for which only a subset of normal samples are available. We describe tuning of the method for detection of mutations in tumors, filtering to remove inherited variants, and comparison of detected mutations to several matched tumor/normal analysis methods. Filtering steps include the use of population variation datasets to remove inherited variants as well a subset of normal samples to remove technical artifacts. We then directly compare mutation detection with tumor-only and tumor-normal approaches using the same sets of samples. Comparisons are performed using an internal targeted gene sequencing dataset (n = 3380) as well as whole exome sequencing data from The Cancer Genome Atlas project (n = 250). Tumor-only mutation detection shows similar recall (43–60%) but lesser precision (20–21%) to current matched tumor/normal approaches (recall 43–73%, precision 30–82%) when compared to a “gold-standard” tumor/normal approach. The inclusion of a small pool of normal samples improves precision, although many variants are still uniquely detected in the tumor-only analysis.

Conclusions

A detailed method for somatic mutation detection without matched normal samples enables study of larger numbers of tumor samples, as well as tumor samples for which a matched normal is not available. As sensitivity/recall is similar to tumor/normal mutation detection but precision is lower, tumor-only detection is more appropriate for classification of samples based on known mutations. Although matched tumor-normal analysis is preferred due to higher precision, we demonstrate that mutation detection without matched normal samples is possible for certain applications.

     

Quantification of the paternal allele bias for new germline mutations in the retinoblastoma gene

New germline mutations in the human retinoblastoma gene are known to arise preferentially on paternally derived chromosomes, but the magnitude of that bias has not been measured. We evaluated 49 cases with a new germline mutation and found that in 40 cases (82%) the mutation arose on the paternally derived allele. We also evaluated 48 cases likely to have a somatic initial mutation; in this group the initial mutation arose on paternal or maternal chromosomes with approximately equal frequency. There was no statistically significant difference in the average age of fathers of children with new paternal germline mutations from the average age of fathers of children with new maternal germline mutations or somatic initial mutations. Combining the data with that from previous reports from other groups, the proportion of new germline mutations arising on a paternally derived allele is 85% (based on 72 cases; 95% confidence interval = 76–93%). This number can be useful in the genetic counseling of some families with retinoblastoma. Received: 18 December 1996 / Accepted: 30 April 1997     

Analysis of somatic mutations at human minisatellite loci in tumors and cell lines

Hypervariable human minisatellite loci show a substantial level of germline instability, and spontaneous mutation rates to new length alleles have been measured directly by pedigree analysis. We now show that mutation events altering the number of minisatellite repeat units are not restricted to the germline, but also arise in other tissues. Mutant alleles can be detected at a very low frequency in lymphoblastoid cell lines and at much higher frequencies in clonal tumor cell populations, most particularly in gastrointestinal adenocarcinomas. Mutant alleles in these tumors are usually present at a dosage equal to or greater than that of the progenitor allele, indicating that most or all of the tumor cells carry the same clonally derived mutant allele. As with germline mutation, the incidence of somatic mutations in tumors varies from locus to locus, with the same locus showing the highest level of germline and somatic instability. Most length changes, as those in the germline, are of only a few repeat units; however, very large changes are also observed, implying that such mutations can occur in the absence of meiosis.     

Intervention of somatic mutational events in vivo by a germline defect at the adenine phosphoribosyltransferase locus

Both germline and somatic mutations are known to affect phenotypes of human cells in vivo. In previous studies, we cloned mutant peripheral blood T cells from germline heterozygous humans for adenine phosphoribosyltransferase (APRT) (EC 2.4.2.7) deficiency and found that approximately 1.3 × 10^–4 peripheral T cells had undergone in vivo somatic mutations. Loss of heterozygosity (LOH) was the major cause of the mutations at the APRT locus since approximately 80% of the mutant T cell clones exhibited loss of normal alleles. In the present study, we identified three heterozygous individuals for APRT deficiency (representing two separate families), in whom none of the somatic mutant cells exhibited LOH at the APRT locus. The germline mutant APRT alleles of these heterozygotes from two unrelated families had the same gross DNA abnormalities detectable by Southern blotting. None of the germline mutant APRT alleles so far reported had such gross DNA abnormalities. The data suggest that the germline mutation unique to these heterozygous individuals is associated with the abrogation of LOH in somatic cells. The absence of LOH at a different locus has already been reported in vitro in an established cell line but the present study describes the first such event in vivo in human individuals. Received: 10 May 1996     

Condition‐dependent mutation rates and sexual selection

‘Good genes’ models of sexual selection show that females can gain indirect benefits for their offspring if male ornaments are condition‐dependent signals of genetic quality. Recurrent deleterious mutation is viewed as a major contributor to variance in genetic quality, and previous theoretical treatments of ‘good genes’ processes have assumed that the influx of new mutations is constant. I propose that this assumption is too simplistic, and that mutation rates vary in ways that are important for sexual selection. Recent data have shown that individuals in poor condition can have higher mutation rates, and I argue that if both male sexual ornaments and mutation rates are condition‐dependent, then females can use male ornamentation to evaluate their mate’s mutation rate. As most mutations are deleterious, females benefit from choosing well‐ornamented mates, as they are less likely to contribute germline‐derived mutations to offspring. I discuss some of the evolutionary ramifications of condition‐dependent mutation rates and sexual selection.     

Somatic versus germline mutation processes at minisatellite CEB1 (D2S90) in humans and transgenic mice

The most variable human minisatellites show extreme germline instability dominated by complex intra-allelic rearrangements plus a lower frequency of inter-allelic transfers of repeat units. In contrast, little is known about somatic instability at such loci. We have therefore used single-molecule PCR to analyze mutation at minisatellite CEB1 (D2S90) in human blood DNA. Somatic mutants were rare and involved only relatively simple intra-allelic events, with no bias toward expansions, in sharp contrast to the complex gain-biased rearrangements seen in sperm. Somatic and germline mutation processes were further analyzed in mice transgenic for a cosmid insert containing CEB1. Mutant molecules in transgenic sperm and blood were detected but only at the low frequencies seen in human blood and arose mainly by simple duplications and deletions as seen for somatic mutations in human. These data suggest distinct pathways for germline and somatic CEB1 mutations with germline instability involving recombination-based repair of meiotic double-strand breaks and somatic mutation arising by replication slippage or mitotic recombination. The problem of transferring germline-specific features of minisatellite instability from human to mouse suggests, with other recent observations, that long-range chromatin conformation may be required for the recombination-based mode of germline instability at human minisatellites.     

Germline inactivation of PTEN and dysregulation of the phosphoinositol-3-kinase/Akt pathway cause human Lhermitte-Duclos disease in adults

Lhermitte-Duclos disease (LDD), or dysplastic gangliocytoma of the cerebellum, is an unusual hamartomatous overgrowth disorder. LDD can be familial or, more commonly, sporadic. It has been only recently recognized that LDD may be associated with Cowden syndrome (CS). Over 80% of patients with CS carry germline mutations in PTEN. It remains unclear whether all cases of LDD, even without features of CS, are caused by germline PTEN mutation and whether somatic PTEN mutation occurs in sporadic LDD. We obtained paraffin-embedded LDD lesions from 18 unselected, unrelated patients and performed mutational analysis of PTEN. Overall, 15 (83%) of 18 samples were found to carry a PTEN mutation. All individuals with mutations were adult-onset patients, but the three without mutations were diagnosed at the ages of 1, 3, and 11 years. Germline DNA was available from six adult-onset cases, and all had germline PTEN mutations. Of these six, two had CS features, one did not have CS features, and three were of unknown CS status. Immunohistochemistry revealed that 75% of the LDD samples had complete or partial loss of PTEN expression accompanied by elevated phosphorylated Akt, specifically in the dysplastic gangliocytoma cells. These data suggest that the loss of PTEN function is sufficient to cause LDD. The high frequency and spectrum of germline PTEN mutations in patients ascertaining by LDD alone confirm that LDD is an important defining feature of CS. Individuals with LDD, even without apparent CS features, should be counseled as in CS.     

Germline and somatic mutations in hMSH6 and hMSH3 in gastrointestinal cancers of the microsatellite mutator phenotype

Hereditary and sporadic gastrointestinal cancer of the microsatellite mutator phenotype (MMP) is characterized by a remarkable genomic instability at simple repeated sequences. The genomic instability is often caused by germline and somatic mutations in DNA mismatch repair (MMR) genes hMSH2 and hMLH1. The MMP can be also caused by epigenetic inactivation of hMLH1. The MMP generates many somatic frameshift mutations in genes containing mononucleotide repeats. We previously reported that in MMP tumors the hMSH6 and hMSH3 MMR genes often carry frameshift mutations in their (C)(8) and (A)(8) tracks, respectively. We proposed that these 'secondary mutator mutations' contribute to a gradual manifestation of the MMP. Here we report the detection of other frameshift, nonsense, and missense mutations in these genes in colon and gastric cancers of the MMP. A germline frameshift mutation was found in hMSH6 in a colon tumor harboring another somatic frameshift mutation. Several germline sequence variants and somatic missense mutations at conserved residues were detected in hMSH6 and only one was detected in hMSH3. Of the three hMSH6 germline variants in conserved residues, one coexisted with a somatic mutation at the (C)(8) track and another had a somatic missense mutation. We suggest that some of these germline and somatic missense variants are pathogenic. While biallelic hMSH6 and hMSH3 frameshift mutations were found in some tumors, many tumors seemed to contain only monoallelic mutations. In some tumors, these somatic monoallelic frameshift mutations at the (C)(8) and (A)(8) tracks were found to coexist with other somatic mutations in the other allele, supporting their functionality during tumorigenesis. However, the low incidence of these additional somatic mutations in hMSH6 and hMSH3 leaves many tumors with only monoallelic mutations. The impact of the frameshift mutations in gene expression was studied by comparative analysis of RNA and protein expression in different tumor cell clones with different genotypes. The results show that the hMSH6 (C)(8) frameshift mutation abolishes protein expression, ruling out a dominant negative effect by a truncated protein. We suggest the functionality of these secondary monoallelic mutator mutations in the context of an accumulative haploinsufficiency model.     

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.