Too Many Mutants with Multiple Mutations

ABSTRACT

It has recently become clear that the classical notion of the random nature of mutation does not hold for the distribution of mutations among genes: most collections of mutants contain more isolates with two or more mutations than predicted by the mutant frequency on the assumption of a random distribution of mutations. Excesses of multiples are seen in a wide range of organisms, including riboviruses, DNA viruses, prokaryotes, yeasts, and higher eukaryotic cell lines and tissues. In addition, such excesses are produced by DNA polymerases in vitro. These “multiples” appear to be generated by transient, localized hypermutation rather than by heritable mutator mutations. The components of multiples are sometimes scattered at random and sometimes display an excess of smaller distances between mutations. As yet, almost nothing is known about the mechanisms that generate multiples, but such mutations have the capacity to accelerate those evolutionary pathways that require multiple mutations where the individual mutations are neutral or deleterious. Examples that impinge on human health may include carcinogenesis and the adaptation of microbial pathogens as they move between individual hosts.     

Mutational clusters generated by non-processive polymerases: A case study using DNA polymerase β in vitro

Available DNA mutational spectra reveal that the number of mutants with multiple mutations (“multiples”) is usually greater than expected from a random distribution of mutations among mutants. These overloads imply the occurrence of non-random clusters of mutations, probably generated during episodes of low-fidelity DNA synthesis. Excess multiples have been reported not only for viruses, bacteria, and eukaryotic cells but also for the DNA polymerases of phages T4 and RB69 in vitro. In the simplest case of a purified polymerase, non-random clusters may be generated by a subfraction of phenotypic variants able to introduce more errors per cycle of DNA synthesis than the normal enzyme. According to this hypothesis, excess multiples are not expected with non-processive polymerases even if they harbor rare mutator variants. DNA polymerase β (Pol β) is a mammalian DNA-repair polymerase with very low processivity. Although several Pol β mutational spectra have been described, there is conflicting evidence on whether or not excess multiples occur, with spectra based on the HSV-tk system tending to show excess multiples. Excess multiples generated by Pol β or any of its mutants might imply that the excesses of multiples observed in numerous other systems, especially those with processive polymerases, could be artifactual. Here, the distributions of mutations generated by native and recombinant rat Pol β and by the Pol β^Y265C mutator were analyzed in the M13mp2 lacZα system. Our results present no evidence for a significant excess of multiples over the expected numbers with any of the Pol β enzymes tested in this system. The reported excess of Pol β-generated multiples in the HSV-tk system may reflect a reduced efficiency of detection of base substitutions that cause weak phenotypes, which in turn may artifactually increase the frequency of multiples.     

Multiple mutations in a shuttle vector modified by ultraviolet irradiation, (+/-)-7 beta,8 alpha-dihydroxy-9 alpha,10 alpha-epoxy-7,8,9,10-tetrahydrobenzo[a]pyrene, and aflatoxin B(1) have different properties than single mutations and may be generated during translesion synthesis

Shuttle vector-based systems are extensively employed to study the mutational properties of various mutagens in mammalian cells. Such vectors are designed for the detection of point mutations, that is small deletions and single base and tandem substitutions. However, mutant target genes carrying two or more point mutations, referred to as multiple mutations, can also be found in various proportions depending on the mutagen and the cells used. To evaluate the frequency and characteristics of multiple mutations, we used a system where the plasmid, pYZ289, was treated by ultraviolet irradiation, aflatoxin B(1) or (+/-)-7 beta,8 alpha-dihydroxy-9 alpha, 10 alpha-epoxy-7,8,9,10-tetrahydrobenzo[a]pyrene before transfection into mouse fibroblast cells. The kinds of mutations and the mutational spectra were different for single and multiple mutations. In addition, in at least 75% of the cases, mutations of multiples appeared to arise in the same strand. Furthermore, mutational spectra for multiple mutations were different for 5' and 3' members of multiple sets. These observations suggest that multiple mutations arise via a different mechanism than single mutations. Moreover, these findings suggest that multiples arise during translesion DNA synthesis and involve an error-prone polymerase able to introduce a base opposite misinstructive or noninstructional DNA lesions and subject to subsequent misincorporation errors.     

In Vitro Direct Repeats-mediated Deletion During PCR Amplification

While conducting our research on mutations in the human blood platelet glycoprotein Ib-alpha (GPIbalpha) gene, we detected an unusual deletion of 84 bp. This deletion took place in vitro, during PCR and between two direct repeats. It was observed that the deletion could be detected either by the direct sequencing of the PCR product or after the latter's cloning into a plasmid. After observing a series of four sequenced clones from the same individual, we noticed that while three had the same 84-bp deletion, the fourth exhibited a shorter one. We also noted that there were no cases wherein both deleted and undeleted amplicons coexisted and that several point mutations occurred in the sequence surrounding the deletion. Such Taq errors are statistically more frequent in the "deletion prone DNA" than usual. Interestingly, the deletion was observed only in a DNA, which we call here "deletion prone DNA", whose structure might have been particularly reorganized. Indeed, the mung bean nuclease pre-treatment of this DNA prior to PCR prevented the deletion, thus strengthening the hypothesis that an intra-strand hairpin structure was involved in the deletion process. Direct repeats-mediated deletion is well known in vivo but this is the first report of such "in vitro direct repeats deletion".     

Random oligonucleotide mutagenesis: application to a large protein coding sequence of a major histocompatibility complex class I gene, H-2DP.

We have used random oligonucleotide mutagenesis (or saturation mutagenesis) to create a library of point mutations in the alpha 1 protein domain of a Major Histocompatibility Complex (MHC) molecule. This protein domain is critical for T cell and B cell recognition. We altered the MHC class I H-2DP gene sequence such that synthetic mutant alpha 1 exons (270 bp of coding sequence), which contain mutations identified by sequence analysis, can replace the wild type alpha 1 exon. The synthetic exons were constructed from twelve overlapping oligonucleotides which contained an average of 1.3 random point mutations per intact exon. DNA sequence analysis of mutant alpha 1 exons has shown a point mutant distribution that fits a Poisson distribution, and thus emphasizes the utility of this mutagenesis technique to "scan" a large protein sequence for important mutations. We report our use of saturation mutagenesis to scan an entire exon of the H-2DP gene, a cassette strategy to replace the wild type alpha 1 exon with individual mutant alpha 1 exons, and analysis of mutant molecules expressed on the surface of transfected mouse L cells.     

The Intercellular Distribution of Mutations Induced in Oocytes of DROSOPHILA MELANOGASTER by Chemical and Physical Mutagens

When females of Drosophila melanogaster are treated with chemical or physical mutagens, not only in one but also in both of the two homologous X chromosomes of a given oocyte, a recessive sex-linked lethal mutation may be induced. A method is described that discriminates between such "single" and "double mutations". A theory is developed to show how a comparison between the expected and the observed frequency of double mutations yields an indication of the intercellular distribution (random or nonrandom) of recessive lethal mutations induced by mutagenic agents in oocytes and, consequently, of the distribution (homogeneous or nonhomogeneous) of those agents.--Three agents were tested: FUdR (12.5, 50.5 and 81.0 micrograms/ml), mitomycin C (130.0 micrograms/ml) and X rays (2000 R, 150 kV). After FUdR feeding, no increase in the mutation frequency usually observed in D. melanogaster without mutagenic treatment was obtained (u = 0.13%, namely three single mutations among 2332 chromosomes tested). After mitomycin C feeding, 104 single and three double mutations were obtained. All of the 50 mutations observed after X irradiation were single mutations. The results obtained in the mitomycin C and radiation experiments favor the assumption of a random intercellular distribution of recessive lethal mutations induced by these two agents in oocytes of D. melanogaster. Reasons are discussed why for other types of mutagenic agents nonrandom distributions may be observed with our technique.     

Codon usage bias as a function of generation time and life expectancy

It has recently been demonstrated that human natural codon usage bias is optimized towards a higher buffering capacity to mutations (measured as the tendency of single point mutations in a DNA sequence to yield the same or similar amino acids) compared to random sequences. In this work, we investigate this phenomenon further by analyzing the natural DNA of four different species (human, mouse, zebrafish and fruit fly) to determine whether such a tolerance to mutations is correlated with the life span and age of sexual maturation for the corresponding organisms. We also propose a new measure to quantify the buffering capacity of a DNA sequence to mutations that takes into account the observed mutation rates within every genome and the effect of the corresponding mutation.Our results suggest there is a propensity for tolerance to mutations that is positively correlated with the life expectancy of the considered organisms. Moreover, random sequences that are constrained to produce the same protein as the naturally occurring sequences are found to be more buffered than completely random sequences while being less buffered than the natural sequences. These results suggest that optimization toward protective mechanisms tolerant to mutations is correlated with both life expectancy and age to sexual maturity at both the levels of codon usage bias and the bias of the natural sequence of codons itself.     

Random movement pattern of the sea urchin Strongylocentrotus droebachiensis

We describe the fine-scale movement of the sea urchin Strongylocentrotus droebachiensis based on analyses of video recordings of undisturbed individuals in the two habitats which mainly differed in food availability, urchin barrens and grazing front. Urchin activity decreased as urchin density increased. Individuals alternated between moving and being stationary and their behaviour did not appear to be affected by either current velocity (within the range from 0 to 15 cm s^− 1) and temperature (2.3 to 6.0 °C). Movement of individuals at each location was compared to that predicted by a random walk model. Mean move length (linear distance between two stationary periods), turning angle and net squared displacement were calculated for each individual. The distribution of turning angles was uniform at each location and there was no evidence of a relationship between urchin density and either move length or urchin velocity. The random model predicted a higher dispersal rate at locations with low urchin densities, such as barrens habitats. However, the movement was sometimes greater or less than predicted by the model, suggesting the influence of local environmental factors. The deviation of individual paths from the model revealed that urchins can be stationary or adopt a local (displacement less than random), random or directional movement. The net daily distance displaced on the barrens, predicted by a random walk model, was similar to the observed movement recorded in our previous study of tagged urchins at one site, but less than that observed at a second site. We postulate that the random dispersal of urchins allows individuals on barrens to reach the kelp zone where food is more abundant although the time required to reach the kelp zone may be considerable (months to years). Urchins decrease their rate of dispersal once they reach the kelp zone so that they likely remain close to this abundant food sources for long periods.     

The mutational specificity of DNA polymerase-beta during in vitro DNA synthesis. Production of frameshift, base substitution, and deletion mutations

The frequency and specificity of mutations produced in vitro by eucaryotic DNA polymerase-beta have been determined in a forward mutation assay using a 250-base target sequence in M13mp2 DNA. Homogeneous DNA polymerase-beta, isolated from four different sources, produces mutations at a frequency of 4-6%/single round of gap-filling DNA synthesis. DNA sequence analyses of 460 independent mutants resulting from this error-prone DNA synthesis demonstrate a wide variety of mutational events. Frameshift and base substitutions are made at approximately equal frequency and together comprise about 90% of all mutations. Two mutational "hot spots" for frameshift and base substitution mutations were observed. The characteristics of the mutations at these sites suggest that certain base substitution errors result from dislocation of template bases rather than from direct mispair formation by DNA polymerase-beta. When considering the entire target sequence, single-base frameshift mutations occur primarily in runs of identical bases, usually pyrimidines. The loss of a single base occurs 20-80 times more frequently than single-base additions and much more frequently than the loss of two or more bases. Base substitutions occur at many sites throughout the target, representing a wide spectrum of mispair formations. Averaged over a large number of phenotypically detectable sites, the base substitution error frequency is greater than one mistake for every 5000 bases polymerized. Large deletion mutations are also observed, at a frequency more than 10-fold over background, indicating that purified DNA polymerases alone are capable of producing such deletions. These data are discussed in relation to the physical and kinetic properties of the purified enzymes and with respect to the proposed role for this DNA polymerase in vivo.     

Clonal expansions of mitochondrial genomes: implications for in vivo mutational spectra

It is often assumed mutant frequencies, as measured in a DNA sample, faithfully represent basic mutation rates associated with these mutations. This paradigm was extremely helpful for in vitro studies of the mechanisms of mutagenesis/repair and causes of mutations. However, in vivo, mutant fractions appear to vary dramatically and randomly from sample to sample. It's unlikely that basic mutational rates vary so much. Such variations are probably caused by clonal expansions of mutants within tissue. Whether a particular tissue sample includes an expansion or not, is a matter of chance, which explains the observed random fluctuations of mutant fractions. Well-known examples of clonal expansions involve pathological conditions such as cancer or mitochondrial disease. It is less appreciated that even in normal tissue, expansions of somatic mutants create local deviations from the "expected" mutant frequencies. The sizes of clonal expansions appear to span a wide range and thus, may affect samples of various sizes, from individual cells to individuals. In conclusion, human body appears to be a sort of a "gambling ground" for clonally expanding mutants. We speculate that expansion of early mutants rather than de novo mutation at old age may be the major source of at least some aging-specific mutants in our bodies.     

In vivo site-directed mutagenesis using oligonucleotides

Functional characterization of the genes of higher eukaryotes has been aided by their expression in model organisms and by analyzing site-specific changes in homologous genes in model systems such as the yeast Saccharomyces cerevisiae. Modifying sequences in yeast or other organisms such that no heterologous material is retained requires in vitro mutagenesis together with subcloning. PCR-based procedures that do not involve cloning are inefficient or require multistep reactions that increase the risk of additional mutations. An alternative approach, demonstrated in yeast, relies on transformation with an oligonucleotide, but the method is restricted to the generation of mutants with a selectable phenotype. Oligonucleotides, when combined with gap repair, have also been used to modify plasmids in yeast; however, this approach is limited by restriction-site availability. We have developed a mutagenesis approach in yeast based on transformation by unpurified oligonucleotides that allows the rapid creation of site-specific DNA mutations in vivo. A two-step, cloning-free process, referred to as delitto perfetto, generates products having only the desired mutation, such as a single or multiple base change, an insertion, a small or a large deletion, or even random mutations. The system provides for multiple rounds of mutation in a window up to 200 base pairs. The process is RAD52 dependent, is not constrained by the distribution of naturally occurring restriction sites, and requires minimal DNA sequencing. Because yeast is commonly used for random and selective cloning of genomic DNA from higher eukaryotes such as yeast artificial chromosomes, the delitto perfetto strategy also provides an efficient way to create precise changes in mammalian or other DNA sequences.     

Molecular evolution of minisatellites in hemiascomycetous yeasts

Minisatellites are DNA tandem repeats exhibiting size polymorphism among individuals of a population. This polymorphism is generated by two different mechanisms, both in human and yeast cells, "replication slippage" during S-phase DNA synthesis and "repair slippage" associated to meiotic gene conversion. The Saccharomyces cerevisiae genome contains numerous natural minisatellites. They are located on all chromosomes without any obvious distribution bias. Minisatellites found in protein-coding genes have longer repeat units and on the average more repeat units than minisatellites in noncoding regions. They show an excess of cytosines on the coding strand, as compared to guanines (negative GC skew). They are always multiples of three, encode serine- and threonine-rich amino acid repeats, and are found preferably within genes encoding cell wall proteins, suggesting that they are positively selected in this particular class of genes. Genome-wide, there is no statistically significant association between minisatellites and meiotic recombination hot spots. In addition, minisatellites that are located in the vicinity of a meiotic hot spot are not more polymorphic than minisatellites located far from any hot spot. This suggests that minisatellites, in S. cerevisiae, evolve probably by strand slippage during replication or mitotic recombination. Finally, evolution of minisatellites among hemiascomycetous yeasts shows that even though many minisatellite-containing genes are conserved, most of the time the minisatellite itself is not conserved. The diversity of minisatellite sequences found in orthologous genes of different species suggests that minisatellites are differentially acquired and lost during evolution of hemiascomycetous yeasts at a pace faster than the genes containing them.     

Special factors in biological strings

Biological macromolecules such as DNA, RNA, and proteins can be regarded as finite sequences of symbols (or words) over a finite alphabet. In this paper, we refer to DNA (RNA) sequences which are words on a four-letter alphabet. A comparison is made between some "genes", or fragments of them, with random sequences or random reshuffled sequences on the same alphabet and having the same length. Some combinatorial techniques of analysis of finite words are developed. A crucial role in the comparison is played by the so-called special factors of a given word. In all the analysed DNA (RNA) fragments the distribution on the length of the number of right (left) special factors differs, in a very typical way, from the corresponding distribution in a string on the same alphabet and having the same length generated by a random source or obtained by making a random alteration (=shuffling) of the original string. This kind of change is irrespective of the length in the range that we have considered <2650 bp and of the phylogenetic origin of the fragment.     

Frequency-dependent selection, metrical characters and molecular evolution

Computer models of selection acting on a quantitative character show that a combination of frequency-dependent and stabilizing selection can maintain many polymorphisms among the genes that determine the character. The models also show that the random order of mutations can give rise to selectively driven stochastic effects that are sometimes more important than random genetic drift. They suggest simple explanations for patterns of divergence between populations and species, and for apparent discrepancies between the rates of morphological and molecular evolution. They point towards a selective theory of 'molecular clocks'.     

Nonrandom chromosomal distribution of spontaneous breakpoints and satellite DNA clusters in two geographically distant populations of Chironomus riparius (Diptera: Chironomidae)

Two geographically distant populations of Chironomus riparius (syn. C. thummi) from two environmentally polluted sites (Santena, Italy and Varna, Bulgaria) show numerous somatic and inherited chromosomal aberrations (inversions, deletions and deficiencies). Fifty-five percent of the observed breakpoints occurred in at least two larvae from both populations. Breakpoints occurring twice or more were considered as common structural chromosomal breakpoints. We tested whether such common breakpoints in larvae of the two polluted populations had a random chromosomal distribution or occurred preferentially in specific heterochromatic regions. Distribution of common breakpoints was not random, and proximal regions of first and third chromosome had significantly more common breakpoints than distal ones. By FISH we identified and mapped 56 chromosomal sections containing clusters of two tandem-repetitive satellite DNA families called Hinf and Alu elements. Like the common breakpoints, these repetitive DNA clusters appeared to be significantly more abundant in regions of constitutive heterochromatin such as the pericentromeric regions, while in distal sections of chromosomal arms they were rare or absent. Twenty-four out of 45 common breakpoints (i.e., 53.3%) occurred in cytogenetic sections where Alu and Hinf satellite DNA probes hybridized. The frequency of co-localization between common breakpoints and repetitive DNA hybridization signals was significantly higher than expected by chance. We hypothesize that spontaneous or induced breaks occur more frequently in sections containing blocks of repetitive DNA.     

Intrachromosomal gene conversion, linkage, and the evolution of multigene families

The evolution of the probabilities of genetic identity within and between tandemly repeated loci of a multigene family is investigated analytically and numerically. Unbiased intrachromosomal gene conversion, equal crossing over, random genetic drift, and mutation to new alleles are incorporated. Generations are discrete and nonoverlapping; the diploid, monoecious population mates at random. Under the restriction that there is at most one crossover in the multigene family per individual per generation, the dependence on location of the probabilities of identity is treated exactly. In the "homogeneous" approximation to this "exact" model, end effects are disregarded; in the "exchangeable" approximation, to which all previous work was confined, all position dependence is neglected. Numerical results indicate that the exchangeable and homogeneous models are both qualitatively correct, the exchangeable model is sometimes too inaccurate for quantitative conclusions, and the homogeneous model is always more accurate than the exchangeable one and is always sufficiently accurate for quantitative conclusions.     

DNA associated with nucleosomes in plants.

50 to 55% of tobacco and barley nuclear DNA is accessible to micrococcal endonuclease digestion. The DNA fragments resulting from a mild endonuclease treatment are multiples of a basic unit of 194 +/- 6 base pairs in tobacco and 195 +/- 6 base pairs in barley. After extensive digestion, a DNA fragment of approximately 140 base pairs is predominant. Hence the "extra-core" or "linker"-DNA is 55 base pairs long. Other fragments having 158 and less than 140 base pairs are present as well. Treatment with DNase I results in multiples of 10 bases when analysed under denaturating conditions. These results show that the general organization of the DNA within the nucleosomes is about the same in higher plants as in other higher eukaryotes.     

Genetic variability of Hong Kong (H3N2) influenza viruses: spontaneous mutations and their location in the viral genome

The genetic heterogeneity of five influenza A (H3N2) strains isolated between 1968 and 1977 has been estimated by T1-oligonucleotide fingerprinting of 32P-labeled viral RNA. Assuming that the large T1-resistant oligonucleotides represent a random sample of the viral RNA, the genetic differences observed would affect 0.3 to 10.7% of the RNA positions of the genes studied, depending on the pair of viruses considered. A smaller degree of genetic heterogeneity was observed when six coetaneous viral samples were compared. The distribution of spontaneous mutations among the viral genes was studied by fingerprinting individual RNA segments isolated either by gel electrophoresis or hybridization with plasmids containing influenza-specific DNA sequences. No statistically significant differences were detected in the distribution of mutations among the viral genes studied. The mutation frequency at the hemagglutinin RNA region coding for the HA1 subunit was found to be two times higher than that at the region encoding that HA2 subunit. Our results suggest that the antigenic variability of influenza viruses may be a consequence of a general genetic variability which effects many of the viral genes.