Hypermutation at A-T base pairs: the A nucleotide replacement spectrum is affected by adjacent nucleotides and there is no reverse complementarity of sequences flanking mutated A and T nucleotides

Hypermutation is thought to be a two-phase process. The first phase is via the action of activation-induced cytidine deaminase (AID), which deaminates C nucleotides in WRC motifs. This results in the RGYW/WRCY hot spot motifs for mutation from G and C observed in vivo. The resemblance between the hot spot for C mutations and the reverse complement of that for G mutations implies a process acting equally on both strands of DNA. The second phase of hypermutation generates mutations from A and T and exhibits strand bias, with more mutations from A than T. Although this does not concur with the idea of one mechanism acting equally on both strands, it has been suggested that the AT mutator also has a reversible motif; WA/TW. We show here that the motifs surrounding the different substitutions from A vary significantly; there is no single targeting motif for all A mutations. Sequence preferences associated with mutations from A more likely reflect an influence of adjacent nucleotides over what the A mutates "to." This influence tends toward "like" replacements: Purines (A or G) in the 5' position bias toward replacement by another purine (G), whereas replacement with pyrimidines (C or T) is more likely if the preceding base is also a pyrimidine. There is no reverse complementarity in these observations, in that similar influences of nucleotides adjacent to T are not seen. Hence, WA and TW should not be considered as reverse complement hot spot motifs for A and T mutations.     

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.     

Amino acid changes coded by bacteriophage T4 DNA polymerase mutator mutants. Relating structure to function

Previous studies on the selection of bacteriophage T4 mutator mutants have been extended and a method to regulate the mutator activity of DNA polymerase mutator strains has been developed. The nucleotide changes of 17 bacteriophage T4 DNA polymerase mutations that confer a mutator phenotype and the nucleotide substitutions of several other T4 DNA polymerase mutations have been determined. The most striking observation is that the distribution of DNA polymerase mutator mutations is not random; almost all mutator mutations are located in the N-terminal half of the DNA polymerase. It has been shown that the T4 DNA polymerase shares several regions of homology at the protein sequence level with DNA polymerases of herpes, adeno and pox viruses. From studies of bacteriophage T4 and herpes DNA polymerase mutants, and from analyses of similar protein sequences from several organisms, we conclude that DNA polymerase synthetic activities are located in the C-terminal half of the DNA polymerase and that exonucleolytic activity is located nearer the N terminus.     

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.      

CpG/CpNpG motifs in the coding region are preferred sites for mutagenesis in the breast cancer susceptibility genes

The range of BRCA1/BRCA2 gene mutations is diverse and the mechanism accounting for this heterogeneity is obscure. To gain insight into the endogenous mutational mechanisms involved, we evaluated the association of specific sequences (i.e. CpG/CpNpG motifs, homonucleotides, short repeats) and mutations within the genes. We classified 1337 published mutations in BRCA1 (1765 BRCA2 mutations) for each specific sequence, and employed computer simulation combined with mathematical calculations to estimate the true underlying tendency of mutation occurrence. Interestingly, we found no mutational bias to homonucleotides and repeats in deletions/insertions and substitutions but striking bias to CpG/CpNpG in substitutions in both genes. This suggests that methylation-dependent DNA alterations would be a major mechanism for mutagenesis.     

Simple Sequence Repeats Together with Mismatch Repair Deficiency Can Bias Mutagenic Pathways in Pseudomonas aeruginosa during Chronic Lung Infection

Pseudomonas aeruginosa is an opportunistic pathogen that chronically infects the airways of cystic fibrosis (CF) patients and undergoes a process of genetic adaptation based on mutagenesis. We evaluated the role of mononucleotide G:C and A:T simple sequence repeats (SSRs) in this adaptive process. An in silico survey of the genome sequences of 7 P. aeruginosa strains showed that mononucleotide G:C SSRs but not A:T SSRs were greatly under-represented in coding regions, suggesting a strong counterselection process for G:C SSRs with lengths >5 bp but not for A:T SSRs. A meta-analysis of published whole genome sequence data for a P. aeruginosa strain from a CF patient with chronic airway infection showed that G:C SSRs but not A:T SSRs were frequently mutated during the infection process through the insertion or deletion of one or more SSR subunits. The mutation tendency of G:C SSRs was length-dependent and increased exponentially as a function of SSR length. When this strain naturally became a stable Mismatch Repair System (MRS)-deficient mutator, the degree of increase of G:C SSRs mutations (5-fold) was much higher than that of other types of mutation (2.2-fold or less). Sequence analysis of several mutated genes reported for two different collections, both containing mutator and non-mutator strains of P. aeruginosa from CF chronic infections, showed that the proportion of G:C SSR mutations was significantly higher in mutators than in non-mutators, whereas no such difference was observed for A:T SSR mutations. Our findings, taken together, provide genome-scale evidences that under a MRS-deficient background, long G:C SSRs are able to stochastically bias mutagenic pathways by making the genes in which they are harbored more prone to mutation. The combination of MRS deficiency and virulence-related genes that contain long G:C SSRs is therefore a matter of concern in P. aeruginosa CF chronic infection.     

DNA polymerase eta contributes to strand bias of mutations of A versus T in immunoglobulin genes

DNA polymerase (pol) eta participates in hypermutation of A:T bases in Ig genes because humans deficient for the polymerase have fewer substitutions of these bases. To determine whether polymerase eta is also responsible for the well-known preference for mutations of A vs T on the nontranscribed strand, we sequenced variable regions from three patients with xeroderma pigmentosum variant (XP-V) disease, who lack polymerase eta. The frequency of mutations in the intronic region downstream of rearranged J(H)4 gene segments was similar between XP-V and control clones; however, there were fewer mutations of A:T bases and correspondingly more substitutions of C:G bases in the XP-V clones (p < 10(-7)). There was significantly less of a bias for mutations of A compared with T nucleotides in the XP-V clones compared with control clones, whereas the frequencies for mutations of C and G were identical in both groups. An analysis of mutations in the WA sequence motif suggests that polymerase eta generates more mutations of A than T on the nontranscribed strand. This in vivo data from polymerase eta-deficient B cells correlates well with the in vitro specificity of the enzyme. Because polymerase eta inserts more mutations opposite template T than template A, it would generate more substitutions of A on the newly synthesized strand.     

DNA sequence context affects UV-induced mutagenesis in Escherichia coli

We investigated the effect of altering the DNA sequence surrounding a mutable target site on the production of ultraviolet light (UV) induced mutations. Site-directed base substitutions were incorporated on both sides of a TAA sequence encoding a UAA nonsense defect in the tyrA14 allele of Escherichia coli. This allele is readily revertable by UV and a total of eight different base substitution mutations can be recovered. Five different strains harboring DNA sequences allowing the formation of 5'-TT, 5'-CT and 5'-TA* photoproducts were constructed and exposed to UV. DNA sequence analysis was used to determine the spectrum of the revertants that were recovered. The results showed that changes at the 3'-base of a TT site were predominantly T to C transitions and T to A transversions. However, unlike the TT site, a 5'-CT site produced a relatively high frequency of T to G transversions. In addition, T to A transversions that could not have been targeted by a cyclobutane-type or [6-4]-type pyrimidine dimer were produced; this result suggested that these mutations may be targeted by a TA* photoproduct. Also, a distinct strand bias was noted for two mechanistically identical base substitutions in a strain having a palindromic target sequence; this result may reflect an unequal damage distribution or processing of photoproducts as a consequence of asymmetric DNA replication. Finally, our results show that DNA sequences expected to allow the greatest density of UV-induced DNA damage produce the highest mutation frequencies. Overall, these findings provide new insights regarding the role of DNA photoproducts in UV mutagenesis.     

Evaluation of the flanking nucleotide sequences of sarcomeric hypertrophic cardiomyopathy substitution mutations

Hypertrophic cardiomyopathy (HCM) is a familial myocardial disease with a prevalence of 1 in 500. More than 400 causative mutations have been identified in 13 sarcomeric and myofilament related genes, 350 of these are substitution mutations within eight sarcomeric genes. Within a population, examples of recurring identical disease causing mutations that appear to have arisen independently have been noted as well as those that appear to have been inherited from a common ancestor. The large number of novel HCM mutations could suggest a mechanism of increased mutability within the sarcomeric genes. The objective of this study was to evaluate the most commonly reported HCM genes, beta myosin heavy chain (MYH7), myosin binding protein C, troponin I, troponin T, cardiac regulatory myosin light chain, cardiac essential myosin light chain, alpha tropomyosin and cardiac alpha-actin for sequence patterns surrounding the substitution mutations that may suggest a mechanism of increased mutability. The mutations as well as the 10 flanking nucleotides were evaluated for frequency of di-, tri- and tetranucleotides containing the mutation as well as for the presence of certain tri- and tetranculeotide motifs. The most common substitutions were guanine (G) to adenine (A) and cytosine (C) to thymidine (T). The CG dinucleotide had a significantly higher relative mutability than any other dinucleotide (p<0.05). The relative mutability of each possible trinucleotide and tetranucleotide sequence containing the mutation was calculated; none were at a statistically higher frequency than the others. The large number of G to A and C to T mutations as well as the relative mutability of CG may suggest that deamination of methylated CpG is an important mechanism for mutation development in at least some of these cardiac genes.     

Mutational specificity of the dnaE173 mutator associated with a defect in the catalytic subunit of DNA polymerase III of Escherichia coli.

We developed a system to examine forward mutations that occurred in the rpsL gene of Escherichia coli placed on a multicopy plasmid. Using this system we determined the mutational specificity for a dnaE173 mutator strain in which the editing function of DNA polymerase III is impeded. The frequency of rpsL- mutations increased 32,000-fold, due to the dnaE173 mutator, and 87 independent rpsL- mutations in the mutator strain were analyzed by DNA sequencing, together with 100 mutants recovered from dnaE+ strain, as the control. While half the number of mutations that occurred in the wild-type strain were caused by insertion elements, no such mutations were recovered from the mutator strain. A novel class of mutation, named "sequence substitution" was present in mutants raised in the dnaE173 strain; seven sequence substitutions induced in the mutator strain occurred at six sites, and all were located in quasipalindromic sequences, carrying the GTG or CAC sequence at one or both endpoints. While other types of mutation were found in both strains, single-base frameshifts were the most frequent events in the mutator strain. Thus, the mutator effect on this class of mutation was 175,000-fold. A total of 95% of the single-base frameshifts in the mutator strain were additions, most of which occurred at runs of A or C bases so as to increase the number of identical residues. Base substitutions, the frequency of which was enhanced 25,000-fold by the mutator effect, occurred primarily at several hotspots in the mutator strain, whereas those induced in the wild-type strain were more randomly distributed throughout the rpsL sequence. The dnaE173 mutator also increased the frequency of duplications 28,000-fold. Of the three duplications recovered from the mutator strain, one was a simple duplication, the region of which was flanked by direct repeats. The other duplications were complex, one half part of which was in the inverted orientation of a region containing two sets of inverted repeats. The same duplications were also recovered from the wild-type strain. The present data suggest that dnaE173 is a novel class of mutator that sharply induces sequence-directed mutagenesis, yielding high frequencies of single base frameshifts, duplications with inversions, sequence substitutions and base substitutions at hotspots.     

Base-pair substitutions alter the site-specific mutagenicity of UV and MNNG in the SUP4-o gene of yeast

Yeast strains carrying SUP4-o genes that have base-pair substitutions at hotspots for UV or MNNG mutagenesis were treated with these agents. In both cases, the induced mutation frequencies were substantially reduced. Furthermore, specific substitutions at positions in SUP4-o that had not been mutated by MNNG resulted in the recovery of MNNG-induced mutations at these sites. These results demonstrate that base-pair identity is an important factor determining the site-specific mutagenicity of UV and MNNG in yeast. For UV, our findings suggest that the type of lesion that is induced, but not flanking DNA sequences, plays a role in specifying mutability at the sites examined. In contrast, DNA sequence context seems to be an important factor for MNNG mutagenesis.     

Position of the rearranged V kappa and its 5' flanking sequences determines the location of somatic mutations in the J kappa locus

Somatic hypermutation is known to occur in the VJ kappa exon and its flanking sequences, yet little is known about the hypermutation mechanism or its exact target within the rearranged locus. Mutations may occur at the same frequency, spanning a region from the leader intron to 3' of J kappa 5, regardless of which J is chosen for VJ rearrangement. Another possibility is that mutations may be limited to the rearranged VJ kappa and its immediate flanking sequences. To distinguish between these possibilities, the JC introns of 21 alleles with V kappa rearranged to J kappa 1 were sequenced, and mutations were located. The frequency of mutations was determined for different sections of the intron and compared with the frequencies of mutations found in the JC intron of a set of VJ kappa 5 alleles. The results showed that mutations were concentrated in and around the rearranged VJ, regardless of whether J kappa 1 or J kappa 5 was used. These data imply that the hypermutational mechanism focuses on rearranged V genes.     

Analysis of alkaptonuria (AKU) mutations and polymorphisms reveals that the CCC sequence motif is a mutational hot spot in the homogentisate 1,2 dioxygenase gene (HGO)

We recently showed that alkaptonuria (AKU) is caused by loss-of-function mutations in the homogentisate 1,2 dioxygenase gene (HGO). Herein we describe haplotype and mutational analyses of HGO in seven new AKU pedigrees. These analyses identified two novel single-nucleotide polymorphisms (INV4+31A-->G and INV11+18A-->G) and six novel AKU mutations (INV1-1G-->A, W60G, Y62C, A122D, P230T, and D291E), which further illustrates the remarkable allelic heterogeneity found in AKU. Reexamination of all 29 mutations and polymorphisms thus far described in HGO shows that these nucleotide changes are not randomly distributed; the CCC sequence motif and its inverted complement, GGG, are preferentially mutated. These analyses also demonstrated that the nucleotide substitutions in HGO do not involve CpG dinucleotides, which illustrates important differences between HGO and other genes for the occurrence of mutation at specific short-sequence motifs. Because the CCC sequence motifs comprise a significant proportion (34.5%) of all mutated bases that have been observed in HGO, we conclude that the CCC triplet is a mutational hot spot in HGO.     

The bacteriophage T4 insertion/substitution vector system. A method for introducing site-specific mutations into the virus chromosome

A bacteriophage T4 insertion/substitution vector system has been developed as a means of introducing in vitro generated mutations into the T4 chromosome. The insertion/substitution vector is a 2638-base pair plasmid containing the pBR322 origin of replication and ampicillin resistance determinant, a T4 gene 23 promoter/synthetic supF tRNA gene fusion, and a polylinker with eight unique restriction enzyme recognition sites. A T4 chromosomal "target" DNA sequence is cloned into this vector and mutated by standard recombinant DNA techniques. Escherichia coli cells containing this plasmid are then infected with T4 bacteriophage that carry amber mutations in two essential genes. The plasmid integrates into the T4 chromosome by recombination between the plasmid-borne T4 target sequence and its homologous chromosomal counterpart. The resulting phage, termed "integrants," are selectable by the supF-mediated suppression of their two amber mutations. Thus, although the integrants comprise 1-3% or less of the total phage progeny, growth on a nonsuppressing host permits their direct selection. The pure integrant phage can be either analyzed directly for a possible mutant phenotype or transferred to nonselective growth conditions. In the latter case, plasmid-free phage segregants rapidly accumulate due to homologous recombination between the duplicated target sequences surrounding the supF sequence in each integrant chromosome. A major fraction of these segregants will retain the in vitro generated mutation within their otherwise unchanged chromosomes and are isolated as stable mutant bacteriophage. The insertion/substitution vector system thereby allows any in vitro mutated gene to be readily substituted for its wild-type counterpart in the bacteriophage T4 genome.     

Flanking nucleotide specificity for DNA mismatch repair-deficient frameshifts within activin receptor 2 (ACVR2)

We previously demonstrated that exonic selectivity for frameshift mutation (exon 10 over exon 3) of ACVR2 in mismatch repair (MMR)-deficient cells is partially determined by 6 nucleotides flanking 5' and 3' of each microsatellite. Substitution of flanking nucleotides surrounding the exon 10 microsatellite with those surrounding the exon 3 microsatellite greatly diminished heteroduplex (A(7)/T(8)) and full (A(7)/T(7)) mutation, while substitution of flanking nucleotides from exon 3 with those from exon 10 enhanced frameshift mutation. We hypothesized that specific individual nucleotide(s) within these flanking sequences control ACVR2 frameshift mutation rates. Only the 3rd nucleotide 5' of the microsatellite, and 3rd, 4th, and 5th nucleotides 3' of the microsatellite were altered from the native flanking sequences and these locations were individually altered (sites A, B, C, and D, respectively). Constructs were cloned +1bp out-of-frame of EGFP, allowing a -1bp frameshift to express EGFP. Plasmids were stably transfected into MMR-deficient cells. Non-fluorescent cells were sorted, cultured for 35 days, and harvested for flow cytometry and DNA-sequencing. Site A (C to T) and B (G to C) in ACVR2 exon 10 decreased both heteroduplex and full mutant as much as the construct containing all 4 alterations. For ACVR2 exon 3, site A (T to C), C (A to G), and D (G to C) are responsible for increased heteroduplex formation, whereas site D is responsible for full mutant formation by ACVR2 exon 10 flanking sequences. Exonic selectivity for frameshift mutation within ACVR2's sequence context is specifically controlled by individual nucleotides flanking each microsatellite.     

Escherichia coli DNA polymerase IV mutator activity: genetic requirements and mutational specificity

The dinB gene of Escherichia coli is known to be involved in the untargeted mutagenesis of lambda phage. Recently, we have demonstrated that this damage-inducible and SOS-controlled gene encodes a novel DNA polymerase, DNA Pol IV, which is able to dramatically increase the untargeted mutagenesis of F' plasmid. At the amino acid level, DNA Pol IV shares sequence homologies with E. coli UmuC (DNA Pol V), Rev1p, and Rad30p (DNA polymerase eta) of Saccharomyces cerevisiae and human Rad30A (XPV) proteins, all of which are involved in translesion DNA synthesis. To better characterize the Pol IV-dependent untargeted mutagenesis, i.e., the DNA Pol IV mutator activity, we analyzed the genetic requirements of this activity and determined the forward mutation spectrum generated by this protein within the cII gene of lambda phage. The results indicated that the DNA Pol IV mutator activity is independent of polA, polB, recA, umuDC, uvrA, and mutS functions. The analysis of more than 300 independent mutations obtained in the wild-type or mutS background revealed that the mutator activity clearly promotes single-nucleotide substitutions as well as one-base deletions in the ratio of about 1:2. The base changes were strikingly biased for substitutions toward G:C base pairs, and about 70% of them occurred in 5'-GX-3' sequences, where X represents the base (T, A, or C) that is mutated to G. These results are discussed with respect to the recently described biochemical characteristics of DNA Pol IV.     

The nucleotide-replacement spectrum under somatic hypermutation exhibits microsequence dependence that is strand-symmetric and distinct from that under germline mutation

Somatic mutation is a fundamental component of acquired immunity. Although its molecular basis remains undetermined, the sequence specificity with which mutations are introduced has provided clues to the mechanism. We have analyzed data representing over 1700 unselected mutations in V gene introns and nonproductively rearranged V genes to identify the sequence specificity of the mutation spectrum-the distribution of resultant nucleotides. In other words, we sought to determine what effects the neighboring bases have on what a given base mutates "to." We find that both neighboring bases have a significant effect on the mutation spectrum. Their influences are complicated, but much of the effect can be characterized as enhancing homogeneity of the mutated DNA sequence. In contrast to what has been reported for the sequence specificity of the "targeting" mechanism, that of the spectrum is notably symmetric under complementation, indicating little if any strand bias. We compared the spectrum to that found previously for germline mutations as revealed by analyzing pseudogene sequences. We find that the influences of nearest neighbors are quite different in the two datasets. Altogether, our findings suggest that the mechanism of somatic hypermutation is complex, involving two or more stages: introduction of mis-pairs and their subsequent resolution, each with distinct sequence specificity and strand bias.     

Spontaneous Mutation at the Mtr Locus in Neurospora: The Molecular Spectrum in Wild-Type and a Mutator Strain

Sequence analysis of 34 mtr mutations has yielded the first molecular spectrum of spontaneous mutants in Neurospora crassa. The great majority of the mutations are base substitutions (48%) or deletions (35%). In addition, sequence analysis of the entire mtr region, including the 1472-base pair open reading frame and 1205 base pairs of flanking DNA, was performed in both the Oak Ridge and Mauriceville strains of Neurospora, which are known to be divergent at the DNA level. Sixteen sequence differences between these two strains have been found in the mtr region, with 13 of these in DNA flanking the open reading frame. The differences consisted of base substitutions and small frameshifts at monotonic runs. This set of sequence differences has allowed a comparison of mutations in unselected DNA to those mutations that produce a phenotypic signal. We have isolated a mutator strain (mut-1) of Neurospora in which the spontaneous mutation rate at various loci is as much as 80-fold higher than in the non-mutator (wild type). Twenty-one mtr mutations in the mutator background have been sequenced and compared to the non-mutator spectrum, revealing a striking increase in -1 frameshift mutations. These frameshifts occur exclusively within or adjacent to monotonic runs and can be explained by small slippage events during DNA replication. This argues for a role of the mut-1 gene in this process.