Use of mutation spectra analysis software

Antibody Structure and Function: The Basis for Engineering Therapeutics

Antibodies and antibody-derived macromolecules have established themselves as the mainstay in protein-based therapeutic molecules (biologics). Our knowledge of the structure–function relationships of antibodies provides a platform for protein engineering that has been exploited to generate a wide range of biologics for a host of therapeutic indications. In this review, our basic understanding of the antibody structure is described along with how that knowledge has leveraged the engineering of antibody and antibody-related therapeutics having the appropriate antigen affinity, effector function, and biophysical properties. The platforms examined include the development of antibodies, antibody fragments, bispecific antibody, and antibody fusion products, whose efficacy and manufacturability can be improved via humanization, affinity modulation, and stability enhancement. We also review the design and selection of binding arms, and avidity modulation. Different strategies of preparing bispecific and multispecific molecules for an array of therapeutic applications are included.

UV-induced mutations in epidermal cells of mice defective in DNA polymerase η and/or ι

Xeroderma pigmentosum variant (XP-V) is a human rare inherited recessive disease, predisposed to sunlight-induced skin cancer, which is caused by deficiency in DNA polymerase η (Polη). Polη catalyzes accurate translesion synthesis (TLS) past pyrimidine dimers, the most prominent UV-induced lesions. DNA polymerase ι (Polι) is a paralog of Polη that has been suggested to participate in TLS past UV-induced lesions, but its function in vivo remains uncertain. We have previously reported that Polη-deficient and Polη/Polι double-deficient mice showed increased susceptibility to UV-induced carcinogenesis. Here, we investigated UV-induced mutation frequencies and spectra in the epidermal cells of Polη- and/or Polι-deficient mice. While Polη-deficient mice showed significantly higher UV-induced mutation frequencies than wild-type mice, Polι deficiency did not influence the frequencies in the presence of Polη. Interestingly, the frequencies in Polη/Polι double-deficient mice were statistically lower than those in Polη-deficient mice, although they were still higher than those of wild-type mice. Sequence analysis revealed that most of the UV-induced mutations in Polη-deficient and Polη/Polι double-deficient mice were base substitutions at dipyrimidine sites. An increase in UV-induced mutations at both G:C and A:T pairs associated with Polη deficiency suggests that Polη contributes to accurate TLS past both thymine- and cytosine-containing dimers in vivo. A significant decrease in G:C to A:T transition in Polη/Polι double-deficient mice when compared with Polη-deficient mice suggests that Polι is involved in error-prone TLS past cytosine-containing dimers when Polη is inactivated.

Comparative mutational analyses of influenza A viruses

The error-prone RNA-dependent RNA polymerase (RdRP) and external selective pressures are the driving forces for RNA viral diversity. When confounded by selective pressures, it is difficult to assess if influenza A viruses (IAV) that have a wide host range possess comparable or distinct spontaneous mutational frequency in their RdRPs. We used in-depth bioinformatics analyses to assess the spontaneous mutational frequencies of two RdRPs derived from human seasonal (A/Wuhan/359/95; Wuhan) and H5N1 (A/Vietnam/1203/04; VN1203) viruses using the mini-genome system with a common firefly luciferase reporter serving as the template. High-fidelity reverse transcriptase was applied to generate high-quality mutational spectra which allowed us to assess and compare the mutational frequencies and mutable motifs along a target sequence of the two RdRPs of two different subtypes. We observed correlated mutational spectra (τ correlation P < 0.0001), comparable mutational frequencies (H3N2:5.8 ± 0.9; H5N1:6.0 ± 0.5), and discovered a highly mutable motif "(A)AAG" for both Wuhan and VN1203 RdRPs. Results were then confirmed with two recombinant A/Puerto Rico/8/34 (PR8) viruses that possess RdRP derived from Wuhan or VN1203 (RG-PR8×Wuhan(PB2, PB1, PA, NP) and RG-PR8×VN1203(PB2, PB1, PA, NP)). Applying novel bioinformatics analysis on influenza mutational spectra, we provide a platform for a comprehensive analysis of the spontaneous mutation spectra for an RNA virus.

Molecular characterization of the piggyBac-like element, a candidate marker for phylogenetic research of Chilo suppressalis (Walker) in China

Background

Transposable elements (TEs, transposons) are mobile genetic DNA sequences. TEs can insert copies of themselves into new genomic locations and they have the capacity to multiply. Therefore, TEs have been crucial in the shaping of hosts’ current genomes. TEs can be utilized as genetic markers to study population genetic diversity. The rice stem borer Chilo suppressalis Walker is one of the most important insect pests of many subtropical and tropical paddy fields. This insect occurs in all the rice-growing areas in China. This research was carried out in order to find diversity between C. suppressalis field populations and detect the original settlement of C. suppressalis populations based on the piggyBac-like element (PLE). We also aim to provide insights into the evolution of PLEs in C. suppressalis and the phylogeography of C. suppressalis.

Results

Here we identify a new piggyBac-like element (PLE) in the rice stem borer Chilo suppressalis Walker, which is called CsuPLE1.1 (GenBank accession no. JX294476). CsuPLE1.1 is transcriptionally active. Additionally, the CsuPLE1.1 sequence varied slightly between field populations, with polymorphic indels (insertion/deletion) and hyper-variable regions including the identification of the 3′ region outside the open reading frame (ORF). CsuPLE1.1 insertion frequency varied between field populations. Sequences variation was found between CsuPLE1 copies and varied within and among field populations. Twenty-one different insertion sites for CsuPLE1 copies were identified with at least two insertion loci found in all populations.

Conclusions

Our results indicate that the initial invasion of CsuPLE1 into C. suppressalis occurred before C. suppressalis populations spread throughout China, and suggest that C. suppressalis populations have a common ancestor in China. Additionally, the lower reaches of the Yangtze River are probably the original settlement of C. suppressalis in China. Finally, the CsuPLE1 insertion site appears to be a candidate marker for phylogenetic research of C. suppressalis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12867-014-0028-y) contains supplementary material, which is available to authorized users.

Oncogenic potential is related to activating effect of cancer single and double somatic mutations in receptor tyrosine kinases

Aberrant activation of receptor tyrosine kinases (RTKs) is a common feature of many cancer cells. It was previously suggested that the mechanisms of kinase activation in cancer might be linked to transitions between active and inactive states. Here, we estimate the effects of single and double cancer mutations on the stability of active and inactive states of the kinase domains from different RTKs. We show that singleton cancer mutations destabilize active and inactive states; however, inactive states are destabilized more than the active ones, leading to kinase activation. We show that there exists a relationship between the estimate of oncogenic potential of cancer mutation and kinase activation. Namely, more frequent mutations have a higher activating effect, which might allow us to predict the activating effect of the mutations from the mutation spectra. Independent evolutionary analysis of mutation spectra complements this observation and finds the same frequency threshold defining mutation hotspots. We analyze double mutations and report a positive epistasis and additional advantage of doublets with respect to cancer cell fitness. The activation mechanisms of double mutations differ from those of single mutations and double mutation spectrum is found to be dissimilar to the mutation spectrum of singletons.

Altered spectrum of somatic hypermutation in common variable immunodeficiency disease characteristic of defective repair of mutations

Pathogenic common variable immunodeficiency diseases (CVID) are genetic, usually inherited diseases for which a limited number of genetic defects have been implicated. As CVID presents with a wide range of clinical characteristics, there are likely diverse and for the most part unidentified genetic causes. In some individuals, defects in somatic hypermutation (SHM) have been suggested as the underlying cause of CVID. To address the mechanisms of SHM defects in CVID, we conducted a comprehensive mutational analysis of immunoglobulin heavy chain sequences from CVID patients. We identified several remarkably specific alterations in the spectra of SHM in comparison to healthy individuals. We provide evidence that some CVID cases are associated with defective repair of AID-induced mutations by the DNA mismatch repair (MMR) machinery. Our findings together with reports of increased chromosomal radiosensitivity and associated lymphoproliferative disorders amongst CVID patients, suggest that altered DNA damage repair may be a cause of CVID.

Mutational spectrum of DMD mutations in dystrophinopathy patients: application of modern diagnostic techniques to a large cohort

Mutations in the DMD gene, encoding the dystrophin protein, are responsible for the dystrophinopathies Duchenne Muscular Dystrophy (DMD), Becker Muscular Dystrophy (BMD), and X-linked Dilated Cardiomyopathy (XLDC). Mutation analysis has traditionally been challenging, due to the large gene size (79 exons over 2.2 Mb of genomic DNA). We report a very large aggregate data set comprised of DMD mutations detected in samples from patients enrolled in the United Dystrophinopathy Project, a multicenter research consortium, and in referral samples submitted for mutation analysis with a diagnosis of dystrophinopathy. We report 1,111 mutations in the DMD gene, including 891 mutations with associated phenotypes. These results encompass 506 point mutations (including 294 nonsense mutations) and significantly expand the number of mutations associated with the dystrophinopathies, highlighting the utility of modern diagnostic techniques. Our data supports the uniform hypermutability of CGA>TGA mutations, establishes the frequency of polymorphic muscle (Dp427m) protein isoforms and reveals unique genomic haplotypes associated with "private" mutations. We note that 60% of these patients would be predicted to benefit from skipping of a single DMD exon using antisense oligonucleotide therapy, and 62% would be predicted to benefit from an inclusive multiexonskipping approach directed toward exons 45 through 55.

Mutation patterns in human α-galactosidase A

A way to study the mutation pattern is to convert a 20-letter protein sequence into a scalar protein sequence, because the 20-letter protein sequence is neither vector nor scalar while a promising way to study patterns is in numerical domain. In this study, we use the amino-acid pair predictability to convert α-galactosidase A with its 137 mutations into scalar sequences, and analyse which amino-acid pairs are more sensitive to mutation. Our results show that the unpredictable amino-acid pairs are more sensitive to mutation, and the mutation trend is to narrow the difference between predicted and actual frequency of amino-acid pairs.

Low-fidelity DNA synthesis by the L979F mutator derivative of Saccharomyces cerevisiae DNA polymerase zeta

To probe Pol ζ functions in vivo via its error signature, here we report the properties of Saccharomyces cerevisiae Pol ζ in which phenyalanine was substituted for the conserved Leu-979 in the catalytic (Rev3) subunit. We show that purified L979F Pol ζ is 30% as active as wild-type Pol ζ when replicating undamaged DNA. L979F Pol ζ shares with wild-type Pol ζ the ability to perform moderately processive DNA synthesis. When copying undamaged DNA, L979F Pol ζ is error-prone compared to wild-type Pol ζ, providing a biochemical rationale for the observed mutator phenotype of rev3-L979F yeast strains. Errors generated by L979F Pol ζ in vitro include single-base insertions, deletions and substitutions, with the highest error rates involving stable misincorporation of dAMP and dGMP. L979F Pol ζ also generates multiple errors in close proximity to each other. The frequency of these events far exceeds that expected for independent single changes, indicating that the first error increases the probability of additional errors within 10 nucleotides. Thus L979F Pol ζ, and perhaps wild-type Pol ζ, which also generates clustered mutations at a lower but significant rate, performs short patches of processive, error-prone DNA synthesis. This may explain the origin of some multiple clustered mutations observed in vivo.

Mutation analysis identifies GUCY2D as the major gene responsible for autosomal dominant progressive cone degeneration

PURPOSE

Heterozygous mutations in the GUCY2D gene, which encodes the membrane-bound retinal guanylyl cyclase-1 protein (RetGC-1), have been shown to cause autosomal dominant inherited cone degeneration and cone-rod degeneration (adCD, adCRD). The present study was a comprehensive screening of the GUCY2D gene in 27 adCD and adCRD unrelated families of these rare disorders.

METHODS

Mutation analysis was performed by direct sequencing as well as PCR and subsequent restriction length polymorphism analysis (PCR/RFLP). Haplotype analysis was performed in selected patients by using microsatellite markers.

RESULTS

GUCY2D gene mutations were identified in 11 (40%) of 27 patients, and all mutations clustered to codon 838, including two known and one novel missense mutation: p.R838C, p.R838H, and p.R838G. Haplotype analysis showed that among the studied patients only two of the six analyzed p.R838C mutation carriers shared a common haplotype and that none of the p.R838H mutation carriers did.

CONCLUSIONS

GUCY2D is a major gene responsible for progressive autosomal dominant cone degeneration. All identified mutations localize to codon 838. Haplotype analysis indicates that in most cases these mutations arise independently. Thus, codon 838 is likely to be a mutation hotspot in the GUCY2D gene.

Detection and characterization of DNA variants in the promoter regions of hundreds of human disease candidate genes

Understanding genetic variation might reveal the cause of individual susceptibility to a variety of complex diseases such as asthma, diabetes, and cancer. Current efforts to identify functional DNA variants have essentially been oriented toward single nucleotide polymorphisms (SNPs) found in coding regions of candidate genes since they have direct impact on the structure and function of the affected proteins. Abnormal expression of finely regulated genes could also lead to disequilibria in different metabolic pathways and/or biological processes. Thus investigation of SNPs in the promoter regions (pSNPs) of genes should improve our knowledge of the etiology of complex diseases. Unfortunately, little is known about the nature and the prevalence of pSNPs. We have analyzed 197 genes targeting the promoter region, arbitrarily defined as a 2-kb genomic segment upstream of the transcription initiation site, by screening by dHPLC for the presence of SNPs in a worldwide panel of 40 individuals. As a result 1838 pSNPs were detected, 75% of which modify (by either gain or loss) putative binding sites of known transcription factors. We also examined the distribution of these pSNPs among features such as conserved regions, repeats, and dinucleotides as well as Gene Ontology terms. This report supports the functional relevance of several of the pSNPs investigated and suggests a putative impact on disease susceptibility.

Diversity and function of adaptive immune receptors in a jawless vertebrate

Instead of the immunoglobulin-type antigen receptors of jawed vertebrates, jawless fish have variable lymphocyte receptors (VLRs), which consist of leucine-rich repeat (LRR) modules. Somatic diversification of the VLR gene is shown here to occur through a multistep assembly of LRR modules randomly selected from a large bank of flanking cassettes. The predicted concave surface of the VLR is lined with hypervariable positively selected residues, and computational analysis suggests a repertoire of about 10(14) unique receptors. Lamprey immunized with anthrax spores responded with the production of soluble antigen-specific VLRs. These findings reveal that two strikingly different modes of antigen recognition through rearranged lymphocyte receptors have evolved in the jawless and jawed vertebrates.

iMARS--mutation analysis reporting software: an analysis of spontaneous cII mutation spectra

The sensitivity of any mutational assay is determined by the level at which spontaneous mutations occur in the corresponding untreated controls. Establishing the type and frequency at which mutations occur naturally within a test system is essential if one is to draw scientifically sound conclusions regarding chemically induced mutations. Currently, mutation-spectra analysis is laborious and time-consuming. Thus, we have developed iMARS, a comprehensive mutation-spectrum analysis package that utilises routinely used methodologies and visualisation tools. To demonstrate the use and capabilities of iMARS, we have analysed the distribution, types and sequence context of spontaneous base substitutions derived from the cII gene mutation assay in transgenic animals. Analysis of spontaneous mutation spectra revealed variation both within and between the transgenic rodent test systems Big Blue Mouse, MutaMouse and Big Blue Rat. The most common spontaneous base substitutions were G:C-->A:T transitions and G:C-->T:A transversions. All Big Blue Mouse spectra were significantly different from each other by distribution and nearly all by mutation type, whereas the converse was true for the other test systems. Twenty-eight mutation hotspots were observed across all spectra generally occurring in CG, GA/TC, GG and GC dinucleotides. A mutation hotspot at nucleotide 212 occurred at a higher frequency in MutaMouse and Big Blue Rat. In addition, CG dinucleotides were the most mutable in all spectra except two Big Blue Mouse spectra. Thus, spontaneous base-substitution spectra showed more variation in distribution, type and sequence context in Big Blue Mouse relative to spectra derived from MutaMouse and Big Blue Rat. The results of our analysis provide a baseline reference for mutation studies utilising the cII gene in transgenic rodent models. The potential differences in spontaneous base-substitution spectra should be considered when making comparisons between these test systems. The ease at which iMARS has allowed us to carry out an exhaustive investigation to assess mutation distribution, mutation type, strand bias, target sequences and motifs, as well as predict mutation hotspots provides us with a valuable tool in helping to distinguish true chemically induced hotspots from background mutations and gives a true reflection of mutation frequency.

The cytidine deaminase AID exhibits similar functional properties in yeast and mammals

A recent work published in Molecular Immunology examined the editing activity of activation-induced deaminase (AID) in yeast (Krause, K., Marcu, K.B., Greeve, J., 2006. The cytidine deaminases AID and APOBEC-1 exhibit distinct functional properties in a novel yeast selectable system. Mol. Immunol.). It was proposed that expression of AID in yeast is not sufficient for the generation of point mutations in a highly transcribed gene due to the lack of cofactors for AID-induced somatic hypermutation, which are unique to B cells. It was suggested that, on its own, AID does not have an intrinsic specificity for its target sequences. However, it has been shown previously that expression of the human AID gene in yeast was moderately mutagenic in a wild-type strain and highly mutagenic in an ung1 uracil-DNA glycosylase-deficient strain (Mayorov, V.I., Rogozin, I.B., Adkison, L.R., Frahm, C.R., Kunkel T.A., Pavlov Y. I., 2005. Expression of human AID in yeast induces mutations in context similar to the context of somatic hypermutation at G-C pairs in immunoglobulin genes. BMC Immunol. 6, 10; Poltoratsky, V.P., Wilson, S.H., Kunkel, T.A., Pavlov, Y.I., 2004. Recombinogenic phenotype of human activation-induced cytosine deaminase. J. Immunol. 172, 4308-4313). The vast majority of mutations were at G-C pairs. Mutations showed a clear DNA sequence context specificity which resembled the specificity of somatic hypermutation at G-C pairs in immunoglobulin genes and AID mutation specificity in vitro. The inability to detect mutator effects of AID by Krause et al. is likely to be caused by the use of the wild-type yeast strain and a small sample of clones examined for the presence of mutations. In addition, we show that non-uniformity of the mutation hotspot distribution is a factor potentially decreasing the chances of detecting mutations.

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.

Expression of human AID in yeast induces mutations in context similar to the context of somatic hypermutation at G-C pairs in immunoglobulin genes

Background

Antibody genes are diversified by somatic hypermutation (SHM), gene conversion and class-switch recombination. All three processes are initiated by the activation-induced deaminase (AID). According to a DNA deamination model of SHM, AID converts cytosine to uracil in DNA sequences. The initial deamination of cytosine leads to mutation and recombination in pathways involving replication, DNA mismatch repair and possibly base excision repair. The DNA sequence context of mutation hotspots at G-C pairs during SHM is DGYW/WRCH (G-C is a hotspot position, R = A/G, Y = T/C, W = A/T, D = A/G/T).

Results

To investigate the mechanisms of AID-induced mutagenesis in a model system, we studied the genetic consequences of AID expression in yeast. We constructed a yeast vector with an artificially synthesized human AID gene insert using codons common to highly expressed yeast genes. We found that expression of the artificial hAIDSc gene was moderately mutagenic in a wild-type strain and highly mutagenic in an ung1 uracil-DNA glycosylase-deficient strain. A majority of mutations were at G-C pairs. In the ung1 strain, C-G to T-A transitions were found almost exclusively, while a mixture of transitions with 12% transversions was characteristic in the wild-type strain. In the ung1 strain mutations that could have originated from deamination of the transcribed stand were found more frequently. In the wild-type strain, the strand bias was reversed. DGYW/WRCH motifs were preferential sites of mutations.

Conclusion

The results are consistent with the hypothesis that AID-mediated deamination of DNA is a major cause of mutations at G-C base pairs in immunoglobulin genes during SHM. The sequence contexts of mutations in yeast induced by AID and those of somatic mutations at G-C pairs in immunoglobulin genes are significantly similar. This indicates that the intrinsic substrate specificity of AID itself is a primary determinant of mutational hotspots at G-C base pairs during SHM.

Mutagenesis by transient misalignment in the human mitochondrial DNA control region

To study spontaneous base substitutions in human mitochondrial DNA (mtDNA), we reconstructed the mutation spectra of the hypervariable segments I and II (HVS I and II) using published data on polymorphisms from various human populations. Classification analysis revealed numerous mutation hotspots in HVS I and II mutation spectra. Statistical analysis suggested that strand dislocation mutagenesis, operating in monotonous runs of nucleotides, plays an important role in generating base substitutions in the mtDNA control region. The frequency of mutations compatible with the primer strand dislocation in the HVS I region was almost twice as high as that for template strand dislocation. Frequencies of mutations compatible with the primer and template strand dislocation models are almost equal in the HVS II region. Further analysis of strand dislocation models suggested that an excess of pyrimidine transitions in mutation spectra, reconstructed on the basis of the L-strand sequence, is caused by an excess of both L-strand pyrimidine transitions and H-strand purine transitions. In general, no significant bias toward parent H-strand-specific dislocation mutagenesis was found in the HVS I and II regions.

Mutational-reporter transgenes rescued from mice lacking either Mgmt, or both Mgmt and Msh6 suggest that O6-alkylguanine-induced miscoding does not contribute to the spontaneous mutational spectrum

O6-methylguanine methyltransferase, Mgmt, constitutes the first line of defense against O6-alkylguanine, which can result in G : C to A : T transitions upon DNA replication. Mgmt has been found in organisms as diverse as archaebacteria and mammals. This evolutionary conservation suggests that all organisms may be exposed to either endogenous or environmental alkylating agents. We thus hypothesized that tissues of Mgmt-/- mice would exhibit elevated mutant frequencies. Employing the Big Blue trade mark transgenic system, we evaluated lacI mutants rescued from liver and small intestinal DNA of young Mgmt-/- mice. Interestingly, while there was a small difference between Mgmt-/- mice and controls with respect to lacI mutant frequency, no differences attributable to Mgmt deficiency were apparent in the mutational spectra. Although mutations stemming from O6-guanine alkylations would be predicted to be cumulative, we found no evidence of an Mgmt-dependent alteration in mutation spectrum in DNA samples from 12 month-old mice. To optimize our ability to detect mutations resulting from O6-alkylguanine-induced G : T mismatches, mice with combined deficiencies of Mgmt and the DNA mismatch repair molecule, Msh6, were analysed. In spite of this strategy, we observed no significant differences between Mgmt-/- Msh6-/- and Msh6-/- mouse lacI mutations, except for a trend towards a greater percentage (of total transitions) of G : C to A : T changes in Mgmt-/-Msh6-/- livers. Therefore, despite the striking evolutionary conservation of Mgmt, deficiency of this gene did not significantly impact the spontaneous lacI mutational spectrum in vivo.

Mutation hotspots in the p53 gene in tumors of different origin: correlation with evolutionary conservation and signs of positive selection

We present a classification analysis of the mutation spectra of the p53 gene and construct maps of hotspots for the germline (Li-Fraumein syndrome), different types of tumors and their derived cell lines. While spectra from solid tumors share common hotspots with the germline spectrum, they also contain unique sets of somatic hotspots that are not observed in the germline. All these hotspots correspond to amino acid replacements in the DNA-binding interface of p53. The mutation spectra of lymphomas and cell lines derived from lymphomas and lung cancers contained few hotspots compared to solid tumors. Thus, the distribution of hotspots in the p53 gene appears to depend on the tumor type and cell growth conditions; this specificity is missed by the bulk hotspot analysis. A negative correlation was detected between the amino acid replacement propensity in tumors and evolutionary variability: the hotspots are located in the positions that are highly conserved in p53 and its paralogs, p63 and p73. In all the mutation spectra, substitutions leading to amino acid replacements strongly dominate over silent substitutions, indicating that functional sites evolving under strong purifying selection are subject to intensive positive selection in p53-dependent tumors. These results are compatible with the gain-of-function concept of the role of p53 in tumorigenesis.

Context of deletions and insertions in human coding sequences

We studied the dependence of the rate of short deletions and insertions on their contexts using the data on mutations within coding exons at 19 human loci that cause mendelian diseases. We confirm that periodic sequences consisting of three to five or more nucleotides are mutagenic. Mutability of sequences with strongly biased nucleotide composition is also elevated, even when mutations within homonucleotide runs longer than three nucleotides are ignored. In contrast, no elevated mutation rates have been detected for imperfect direct or inverted repeats. Among known candidate contexts, the indel context GTAAGT and regions with purine-pyrimidine imbalance between the two DNA strands are mutagenic in our sample, and many others are not mutagenic. Data on mutation hot spots suggest two novel contexts that increase the deletion rate. Comprehensive analysis of mutability of all possible contexts of lengths four, six, and eight indicates a substantially elevated deletion rate within YYYTG and similar sequences, which is one of the two contexts revealed by the hot spots. Possible contexts that increase the insertion rate (AT(A/C)(A/C)GCC and TACCRC) and decrease deletion (TATCGC) or insertion (GCGG) rates have also been identified. Two-thirds of deletions remove a repeat, and over 80% of insertions create a repeat, i.e., they are duplications.

Cutting Edge: DGYW/WRCH Is a Better Predictor of Mutability at G:C Bases in Ig Hypermutation Than the Widely Accepted RGYW/WRCY Motif and Probably Reflects a Two-Step Activation-Induced Cytidine Deaminase-Triggered Process

A feature of Ig hypermutation is the presence of hypermutable DNA sequences that are preferentially found in the V regions of Ig genes. Among these, RGYW/WRCY is the most pronounced motif (G:C is a mutable position; R=A/G, Y=C/T, and W=A/T). However, a molecular basis for the high mutability of RGYW was not known until recently. The discovery that activation-induced cytidine deaminase targets the DNA encoding V regions, has enabled the analysis of its targeting properties when expressed outside of the context of hypermutation. We analyzed these data and found evidence that activation-induced cytidine deaminase is the major source of the RGYW mutable motif, but with a new twist: DGYW/WRCH (G:C is the mutable position; D=A/G/T, H=T/C/A) is a better descriptor of the Ig mutation hotspot than RGYW/WRCY. We also found evidence that a DNA repair enzyme may play a role in modifying the sequence of hypermutation hotspots.

Unique error signature of the four-subunit yeast DNA polymerase epsilon

We have purified wild type and exonuclease-deficient four-subunit DNA polymerase epsilon (Pol epsilon) complex from Saccharomyces cerevisiae and analyzed the fidelity of DNA synthesis by the two enzymes. Wild type Pol epsilon synthesizes DNA accurately, generating single-base substitutions and deletions at average error rates of </=2 x 10-5 and </=5 x 10-7, respectively. Pol epsilon lacking 3' --> 5' exonuclease activity is less accurate to a degree suggesting that wild type Pol epsilon proofreads at least 92% of base substitution errors and at least 99% of frameshift errors made by the polymerase. Surprisingly the base substitution fidelity of exonuclease-deficient Pol epsilon is severalfold lower than that of proofreading-deficient forms of other replicative polymerases. Moreover the spectrum of errors shows a feature not seen with other A, B, C, or X family polymerases: a high proportion of transversions resulting from T.dTTP, T.dCTP, and C.dTTP mispairs. This unique error specificity and amino acid sequence alignments suggest that the structure of the polymerase active site of Pol epsilon differs from those of other B family members. We observed both similarities and differences between the spectrum of substitutions generated by proofreading-deficient Pol epsilon in vitro and substitutions occurring in vivo in a yeast strain defective in Pol epsilon proofreading and DNA mismatch repair. We discuss the implications of these findings for the role of Pol epsilon polymerase activity in DNA replication.

The amino-acid mutational spectrum of human genetic disease

The human disease spectrum is compared to the spectra of mutual amino-acid mutation frequencies, non-disease polymorphisms in human genes, and substitutions fixed between species.

Background

Nonsynonymous mutations in the coding regions of human genes are responsible for phenotypic differences between humans and for susceptibility to genetic disease. Computational methods were recently used to predict deleterious effects of nonsynonymous human mutations and polymorphisms. Here we focus on understanding the amino-acid mutation spectrum of human genetic disease. We compare the disease spectrum to the spectra of mutual amino-acid mutation frequencies, non-disease polymorphisms in human genes, and substitutions fixed between species.

Results

We find that the disease spectrum correlates well with the amino-acid mutation frequencies based on the genetic code. Normalized by the mutation frequencies, the spectrum can be rationalized in terms of chemical similarities between amino acids. The disease spectrum is almost identical for membrane and non-membrane proteins. Mutations at arginine and glycine residues are together responsible for about 30% of genetic diseases, whereas random mutations at tryptophan and cysteine have the highest probability of causing disease.

Conclusions

The overall disease spectrum mainly reflects the mutability of the genetic code. We corroborate earlier results that the probability of a nonsynonymous mutation causing a genetic disease increases monotonically with an increase in the degree of evolutionary conservation of the mutation site and a decrease in the solvent-accessibility of the site; opposite trends are observed for non-disease polymorphisms. We estimate that the rate of nonsynonymous mutations with a negative impact on human health is less than one per diploid genome per generation.

Theoretical analysis of mutation hotspots and their DNA sequence context specificity

Mutation frequencies vary significantly along nucleotide sequences such that mutations often concentrate at certain positions called hotspots. Mutation hotspots in DNA reflect intrinsic properties of the mutation process, such as sequence specificity, that manifests itself at the level of interaction between mutagens, DNA, and the action of the repair and replication machineries. The hotspots might also reflect structural and functional features of the respective DNA sequences. When mutations in a gene are identified using a particular experimental system, resulting hotspots could reflect the properties of the gene product and the mutant selection scheme. Analysis of the nucleotide sequence context of hotspots can provide information on the molecular mechanisms of mutagenesis. However, the determinants of mutation frequency and specificity are complex, and there are many analytical methods for their study. Here we review computational approaches for analyzing mutation spectra (distribution of mutations along the target genes) that include many mutable (detectable) positions. The following methods are reviewed: derivation of a consensus sequence, application of regression approaches to correlate nucleotide sequence features with mutation frequency, mutation hotspot prediction, analysis of oligonucleotide composition of regions containing mutations, pairwise comparison of mutation spectra, analysis of multiple spectra, and analysis of "context-free" characteristics. The advantages and pitfalls of these methods are discussed and illustrated by examples from the literature. The most reliable analyses were obtained when several methods were combined and information from theoretical analysis and experimental observations was considered simultaneously. Simple, robust approaches should be used with small samples of mutations, whereas combinations of simple and complex approaches may be required for large samples. We discuss several well-documented studies where analysis of mutation spectra has substantially contributed to the current understanding of molecular mechanisms of mutagenesis. The nucleotide sequence context of mutational hotspots is a fingerprint of interactions between DNA and DNA repair, replication, and modification enzymes, and the analysis of hotspot context provides evidence of such interactions.

Neutral evolution of ten types of mariner transposons in the genomes of Caenorhabditis elegans and Caenorhabditis briggsae

Ten types of mariner transposable elements (232 individual sequences) are present in the completed genomic DNA sequence of Caenorhabditis elegans and the partial sequence of Caenorhabditis briggsae. We analyze these replicated instances of mariner evolution and find that elements of a type have evolved within their genomes under no selection on their transposase genes. Seven of the ten reconstructed ancestral mariners carry defective transposase genes. Selection has acted during the divergence of some ancestral elements. The neutrally-evolving mariners are used to analyze the pattern of molecular evolution in Caenorhabditis. There is a significant mutational bias against transversions and significant variation in rates of change across sites. Deletions accumulate at a rate of 0.034 events/bp per substitution/site, with an average size of 166 bp (173 gaps observed). Deletions appear to obliterate preexisting deletions over time, creating larger gaps. Insertions accumulate at a rate of 0.019 events/bp per substitution/site, with an average size of 151 bp (61 events). Although the rate of deletion is lower than most estimates in other species, the large size of deletions causes rapid elimination of neutral DNA: a mariner's "half-life" (the time by which half an element's sequence should have been deleted) is approximately 0.1 subsitutions/site. This high rate of DNA deletion may explain the compact nature of the nematode genome.

Correlation of somatic hypermutation specificity and A-T base pair substitution errors by DNA polymerase eta during copying of a mouse immunoglobulin kappa light chain transgene

To test the hypothesis that inaccurate DNA synthesis by mammalian DNA polymerase eta (pol eta) contributes to somatic hypermutation (SHM) of Ig genes, we measured the error specificity of mouse pol eta during synthesis of each strand of a mouse Ig kappa light chain transgene. We then compared the results to the base substitution specificity of SHM of this same gene in the mouse. The in vitro and in vivo base substitution spectra shared a number of common features. A highly significant correlation was observed for overall substitutions at A-T pairs but not for substitutions at G-C pairs. Sixteen mutational hotspots at A-T pairs observed in vivo were also found in spectra generated by mouse pol eta in vitro. The correlation was strongest for errors made by pol eta during synthesis of the non-transcribed strand, but it was also observed for synthesis of the transcribed strand. These facts, and the distribution of substitutions generated in vivo, support the hypothesis that pol eta contributes to SHM of Ig genes at A-T pairs via short patches of low fidelity DNA synthesis of both strands, but with a preference for the non-transcribed strand.

Characterization of mutant spectra generated by a forward mutational assay for gene A of Phi X174 from ENU-treated transgenic mouse embryonic cell line PX-2

The sensitivity of in vivo transgenic mutation assays benefits from the sequencing of mutations, although the large number of possible mutations hinders high throughput sequencing. A forward mutational assay exists for Phi X174 that requires an altered, functional Phi X174 protein and therefore should have fewer targets (sense, base-pair substitutions) than forward assays that inactivate a protein. We investigated this assay to determine the number of targets and their suitability for detecting a known mutagen, N-ethyl-N-nitrosourea (ENU). We identified 25 target sites and 33 different mutations in Phi X174 gene A after sequencing over 350 spontaneous and ENU-induced mutants, mostly from mouse embryonic cell line PX-2 isolated from mice transgenic for Phi X174 am3, cs70 (line 54). All six types of base-pair substitution were represented among both the spontaneous and ENU-treated mutant spectra. The mutant spectra from cells treated with 200 and 400 microg/ml ENU were both highly different from the spontaneous spectrum (P < 0.000001) but not from each other. The dose trend was significant (P < 0.0001) for a linear regression of mutant frequencies (R(2) = 0.79), with a ninefold increase in mutant frequency at the 400 microg/ml dose. The spontaneous mutant frequency was 1.9 x 10(-5) and the spontaneous spectrum occurred at 11 target base pairs with 15 different mutations. Thirteen mutations at 12 targets were identified only from ENU-treated cells. Seven mutations had highly significant increases with ENU treatment (P < 0.0001) and 15 showed significant increases. The results suggest that the Phi X174 forward assay might be developed into a sensitive, inexpensive in vivo mutagenicity assay.

Statistical inference of sequence-dependent mutation rates

Several lines of research are now converging towards an integrated understanding of mutational mechanisms and their evolutionary implications. Experimentally, crystal structures reveal the effect of sequence context on polymerase fidelity; large-scale sequencing projects generate vast amounts of sequence polymorphism data; and locus-specific databases are being constructed. Computationally, software and analytical tools have been developed to analyze mutational data, to identify mutational hot spots, and to compare the signatures of mutagenic agents.

Somatic mutation hotspots correlate with DNA polymerase eta error spectrum

Mutational spectra analysis of 15 immunoglobulin genes suggested that consensus motifs RGYW and WA were universal descriptors of somatic hypermutation. Highly mutable sites, "hotspots", that matched WA were preferentially found in one DNA strand and RGYW hotspots were found in both strands. Analysis of base-substitution hotspots in DNA polymerase error spectra showed that 33 of 36 hotspots in the human polymerase eta spectrum conformed to the WA consensus. This and four other characteristics of polymerase eta substitution specificity suggest that errors introduced by this enzyme during synthesis of the nontranscribed DNA strand in variable regions may contribute to strand-specific somatic hypermutagenesis of immunoglobulin genes at A-T base pairs.