A study of the evolution of repeated DNA sequences in primates and the existence of a new class of repetitive sequences in primates.

A new class of lowly repetitive DNA sequences has been detected in the primate genome. The renaturation rate of this sequence class is practically indistinguishable from the renaturation rate of single-copy sequences. Consequently, this lowly repetitive sequence class has not been previously observed in DNA renaturation rate studies. This new sequence class is significant in that it might occupy a major fraction of the primate genome.Based on a study of the thermal stabilities of DNA heteroduplexes constructed from human DNA and either bonnet monkey or galago DNAs, we are able to compare the relative mutation rates of repetitive and single-copy sequences in the primate genome. We find that the mutation rate of short, interspersed repetitive sequences is either less than or approximately equal to the mutation rate of single-copy sequences. This implies that the base sequence of these repetitive sequences is important to their biological function.We also find that numerous mutations have accumulated in interspersed repeated sequences since the divergence of galago and human. These mutations are only recognizable because they occur at specific sites in the repeated sequence rather than at random sites in the sequence. Although interspersed repetitive sequences from human and galago can readily cross-hybridize, these site-specific mutations identify them as being two distinct classes. In contrast, far fewer site-specific mutations have occurred since the divergence of human and monkey.     

The evolution of mutation rates: separating causes from consequences

Natural selection can adjust the rate of mutation in a population by acting on allelic variation affecting processes of DNA replication and repair. Because mutation is the ultimate source of the genetic variation required for adaptation, it can be appealing to suppose that the genomic mutation rate is adjusted to a level that best promotes adaptation. Most mutations with phenotypic effects are harmful, however, and thus there is relentless selection within populations for lower genomic mutation rates. Selection on beneficial mutations can counter this effect by favoring alleles that raise the mutation rate, but the effect of beneficial mutations on the genomic mutation rate is extremely sensitive to recombination and is unlikely to be important in sexual populations. In contrast, high genomic mutation rates can evolve in asexual populations under the influence of beneficial mutations, but this phenomenon is probably of limited adaptive significance and represents, at best, a temporary reprieve from the continual selection pressure to reduce mutation. The physiological cost of reducing mutation below the low level observed in most populations may be the most important factor in setting the genomic mutation rate in sexual and asexual systems, regardless of the benefits of mutation in producing new adaptive variation. Maintenance of mutation rates higher than the minimum set by this "cost of fidelity" is likely only under special circumstances.     

Direct mutation analysis by high-throughput sequencing: from germline to low-abundant, somatic variants

DNA mutations are the source of genetic variation within populations. The majority of mutations with observable effects are deleterious. In humans mutations in the germ line can cause genetic disease. In somatic cells multiple rounds of mutations and selection lead to cancer. The study of genetic variation has progressed rapidly since the completion of the draft sequence of the human genome. Recent advances in sequencing technology, most importantly the introduction of massively parallel sequencing (MPS), have resulted in more than a hundred-fold reduction in the time and cost required for sequencing nucleic acids. These improvements have greatly expanded the use of sequencing as a practical tool for mutation analysis. While in the past the high cost of sequencing limited mutation analysis to selectable markers or small forward mutation targets assumed to be representative for the genome overall, current platforms allow whole genome sequencing for less than $5000. This has already given rise to direct estimates of germline mutation rates in multiple organisms including humans by comparing whole genome sequences between parents and offspring. Here we present a brief history of the field of mutation research, with a focus on classical tools for the measurement of mutation rates. We then review MPS, how it is currently applied and the new insight into human and animal mutation frequencies and spectra that has been obtained from whole genome sequencing. While great progress has been made, we note that the single most important limitation of current MPS approaches for mutation analysis is the inability to address low-abundance mutations that turn somatic tissues into mosaics of cells. Such mutations are at the basis of intra-tumor heterogeneity, with important implications for clinical diagnosis, and could also contribute to somatic diseases other than cancer, including aging. Some possible approaches to gain access to low-abundance mutations are discussed, with a brief overview of new sequencing platforms that are currently waiting in the wings to advance this exploding field even further.     

Rates of molecular evolution and the fraction of nucleotide positions free to vary

Summary Selective constraints on DNA sequence change were incorporated into a model of DNA divergence by restricting substitutions to a subset of nucleotide positions. A simple model showed that both mutation rate and the fraction of nucleotide positions free to vary are strong determinants of DNA divergence over time.When divergence between two species approaches the fraction of positions free to vary, standard methods that correct for multiple mutations yield severe underestimates of the number of substitutions per site. A modified method appropriate for use with DNA sequence, restriction site, or thermal renaturation data is derived taking this fraction into account. The model also showed that the ratio of divergence in two gene classes (e.g., nuclear and mitochondrial) may vary widely over time even if the ratio of mutation rates remains constant.DNA sequence divergence data are used increasingly to detect differences in rates of molecular evolution. Often, variation in divergence rate is assumed to represent variation in mutation rate. The present model suggests that differing divergence rates among comparisons (either among gene classes or taxa) should be interpreted cautiously. Differences in the fraction of nucleotide positions free to vary can serve as an important alternative hypothesis to explain differences in DNA divergence rates.     

Increased transmission of mutations by low-condition females: evidence for condition-dependent DNA repair

Evidence is mounting that mutation rates are sufficiently high for deleterious alleles to be a major evolutionary force affecting the evolution of sex, the maintenance of genetic variation, and many other evolutionary phenomena. Though point estimates of mutation rates are improving, we remain largely ignorant of the biological factors affecting these rates at the individual level. Of special importance is the possibility that mutation rates are condition-dependent with low-condition individuals experiencing more mutation. Theory predicts that such condition dependence would dramatically increase the rate at which populations adapt to new environments and the extent to which populations suffer from mutation load. Despite its importance, there has been little study of this phenomenon in multicellular organisms. Here, we examine whether DNA repair processes are condition-dependent in Drosophila melanogaster. In this species, damaged DNA in sperm can be repaired by maternal repair processes after fertilization. We exposed high- and low-condition females to sperm containing damaged DNA and then assessed the frequency of lethal mutations on paternally derived X chromosomes transmitted by these females. The rate of lethal mutations transmitted by low-condition females was 30% greater than that of high-condition females, indicating reduced repair capacity of low-condition females. A separate experiment provided no support for an alternative hypothesis based on sperm selection.     

Irregularities in genetic variation and mutation rates with environmental stresses

The appearance of new mutations is determined by the equilibrium between DNA error formation and repair. In bacteria like Escherichia coli, stresses are thought shift this balance towards increased mutagenesis. Recent findings, however, suggest a very uneven relationship between stress and mutations. Only a subset of stressful environments increase the net rate of mutation and different forms of nutritional stress (such as oxygen, carbon or phosphorus limitations) result in markedly different mutation rates after similar reductions in growth rate. Moreover, different stresses result in altered mutational spectra, with some increasing transposition and others increasing indel formation. Single-base substitution rates are lower with some stresses than in unstressed bacteria. Indeed, changes to the mix of mutations with stress are more widespread than a marked increase in net mutation rate. Much remains to be learned on how environments have unique mutational signatures and why some stresses are more mutagenic than others. Even beyond stress-induced genetic variation, the fundamental unresolved question in the stress–mutation relationship is the adaptive value of different types of mutations and mutation rates; is transposition, for example, more advantageous under anaerobic conditions? It remains to be investigated whether stress-specific genetic variation impacts on evolvability differentially in distinct environments.     

The construction and evolution of DNA sequences having the same doublet ratios as those of non-satellite, mammalian DNA.

The doublet or nearest-neighbour ratios of the nucleotides in various computer-generated sequences of DNA have been counted to find out which sequences would have the same ratios as those measured for guinea-pig DNA by Russell et al. (1976). Their data shows that the ratio patterns for all nuclear DNA fractions except satellite, ribosomal and tRNA coding DNA are similar irrespective of G+C content and are characterised by the amount of the doublet CpG being less than 30% of that expected on a random basis. To construct and analyse such theoretical sequences, methods have been developed which allow to be counted the doublet frequencies that random DNA and the DNA expected to code for any amino acid (AA) sequence would have were they analysed experimentally. These methods permit the C+G content to be altered and the frequency of the doublet CpG to be lowered without affecting the information stored in the DNA. The former is achieved by selecting codons for a given AA that are high or low in G and C while the latter requires selecting against triplets that either contain CpG or will cause a CpG to occur between codons.The results show that no DNA sequence that we have been able to construct using the unrestricted genetic code has doublet ratios similar to those observed. However, the DNA expected to code for a group of 27 vertebrate proteins (5237 AAs) of diverse functions has doublet ratios virtually identical to those measured experimentally for the 47% G+C fraction, provided that 77·5 % of the CpG is eliminated. The data for the 34–43% G+C fractions are matched well by the protein sequence provided that codons low in G and C are selected and, again, that 77·5% of CpG is eliminated. We have been unable to match the data for the satellite DNA. A perhaps surprising result was that DNA with a random sequence of nucleotides but subjected to the removal of 80% of CpG had doublet ratios that were similar to the experimental data but matched them less well than the doublets of protein-coding DNA. This result probably does no more than emphasise that a significant part of the match of the pattern of doublet ratios of guinea-pig DNA derives from the elimination of CpG.The effect of evolution (random mutation) on the doublet ratios of protein-coding DNA has been investigated by assuming a mutation rate of 3·10^?9 per base per generation (the mutation rate of haemoglobin) and seeing how a great many generations of such mutation affect the doublet ratios. The results show that it will take ~1·5·10⁶ generations for the number of termination codons to double and ~2·10⁷ generations for the doublet ratios to become indistinguishable from those of random DNA. This seems to imply either that selection acts over the whole of the DNA or that the mutation rate of haemoglobin DNA is unusually high. The results, as a whole, support the view that, whether or not all non-satellite DNA actually codes for protein, its sequences are similar to those that would code for proteins.     

Genomic mutation rates: what high‐throughput methods can tell us

High‐throughput DNA analyses are increasingly being used to detect rare mutations in moderately sized genomes. These methods have yielded genome mutation rates that are markedly higher than those obtained using pre‐genomic strategies. Recent work in a variety of organisms has shown that mutation rate is strongly affected by sequence context and genome position. These observations suggest that high‐throughput DNA analyses will ultimately allow researchers to identify trans‐acting factors and cis sequences that underlie mutation rate variation. Such work should provide insights on how mutation rate variability can impact genome organization and disease progression.     

Estimating the frequency of events that cause multiple-nucleotide changes

Existing mathematical models of DNA sequence evolution assume that all substitutions derive from point mutations. There is, however, increasing evidence that larger-scale events, involving two or more consecutive sites, may also be important. We describe a model, denoted SDT, that allows for single-nucleotide, doublet, and triplet mutations. Applied to protein-coding DNA, the SDT model allows doublet and triplet mutations to overlap codon boundaries but still permits data to be analyzed using the simplifying assumption of independence of sites. We have implemented the SDT model for maximum-likelihood phylogenetic inference and have applied it to an alignment of mammalian globin sequences and to 258 other protein-coding sequence alignments from the Pandit database. We find the SDT model's inclusion of doublet and triplet mutations to be overwhelmingly successful in giving statistically significant improvements in fit of model to data, indicating that larger-scale mutation events do occur. Distributions of inferred parameter values over all alignments analyzed suggest that these events are far more prevalent than previously thought. Detailed consideration of our results and the absence of any known mechanism causing three adjacent nucleotides to be substituted simultaneously, however, leads us to suggest that the actual evolutionary events occurring may include still-larger-scale events, such as gene conversion, inversion, or recombination, or a series of rapid compensatory changes.     

Evidence for widespread convergent evolution around human microsatellites

Microsatellites are a major component of the human genome, and their evolution has been much studied. However, the evolution of microsatellite flanking sequences has received less attention, with reports of both high and low mutation rates and of a tendency for microsatellites to cluster. From the human genome we generated a database of many thousands of (AC)_n flanking sequences within which we searched for common characteristics. Sequences flanking microsatellites of similar length show remarkable levels of convergent evolution, indicating shared mutational biases. These biases extend 25–50 bases either side of the microsatellite and may therefore affect more than 30% of the entire genome. To explore the extent and absolute strength of these effects, we quantified the observed convergence. We also compared homologous human and chimpanzee loci to look for evidence of changes in mutation rate around microsatellites. Most models of DNA sequence evolution assume that mutations are independent and occur randomly. Allowances may be made for sites mutating at different rates and for general mutation biases such as the faster rate of transitions over transversions. Our analysis suggests that these models may be inadequate, in that proximity to even very short microsatellites may alter the rate and distribution of mutations that occur. The elevated local mutation rate combined with sequence convergence, both of which we find evidence for, also provide a possible resolution for the apparently contradictory inferences of mutation rates in microsatellite flanking sequences.     

Volatility of Mutator Phenotypes at Single Cell Resolution

Mutator phenotypes accelerate the evolutionary process of neoplastic transformation. Historically, the measurement of mutation rates has relied on scoring the occurrence of rare mutations in target genes in large populations of cells. Averaging mutation rates over large cell populations assumes that new mutations arise at a constant rate during each cell division. If the mutation rate is not constant, an expanding mutator population may contain subclones with widely divergent rates of evolution. Here, we report mutation rate measurements of individual cell divisions of mutator yeast deficient in DNA polymerase ε proofreading and base-base mismatch repair. Our data are best fit by a model in which cells can assume one of two distinct mutator states, with mutation rates that differ by an order of magnitude. In error-prone cell divisions, mutations occurred on the same chromosome more frequently than expected by chance, often in DNA with similar predicted replication timing, consistent with a spatiotemporal dimension to the hypermutator state. Mapping of mutations onto predicted replicons revealed that mutations were enriched in the first half of the replicon as well as near termination zones. Taken together, our findings show that individual genome replication events exhibit an unexpected volatility that may deepen our understanding of the evolution of mutator-driven malignancies.     

How important is DNA replication for mutagenesis?

Rates of mutation and substitution in mammals are generally greater in the germ lines of males. This is usually explained as resulting from the larger number of germ cell divisions during spermatogenesis compared with oogenesis, with the assumption made that mutations occur primarily during DNA replication. However, the rate of cell division is not the only difference between male and female germ lines, and mechanisms are known that can give rise to mutations independently of DNA replication. We investigate the possibility that there are other causes of male-biased mutation. First, we show that patterns of variation at approximately 5,200 short tandem repeat (STR) loci indicate a higher mutation rate in males. We estimate a ratio of male-to-female mutation rates of approximately 1.9. This is significantly greater than 1 and supports a greater rate of mutation in males, affecting the evolution of these loci. Second, we show that there are chromosome-specific patterns of nucleotide and dinucleotide composition in mammals that have been shaped by mutation at CpG dinucleotides. Comparable patterns occur in birds. In mammals, male germ lines are more methylated than female germ lines, and these patterns indicate that differential methylation has played a role in male-biased vertebrate evolution. However, estimates of male mutation bias obtained from both classes of mutation are substantially lower than estimates of cell division bias from anatomical data. This discrepancy, along with published data indicating slipped-strand mispairing arising at STR loci in nonreplicating DNA, suggests that a substantial percentage of mutation may occur in nonreplicating DNA.     

Differential Relationship of DNA Replication Timing to Different Forms of Human Mutation and Variation

Human genetic variation is distributed nonrandomly across the genome, though the principles governing its distribution are only partially known. DNA replication creates opportunities for mutation, and the timing of DNA replication correlates with the density of SNPs across the human genome. To enable deeper investigation of how DNA replication timing relates to human mutation and variation, we generated a high-resolution map of the human genome’s replication timing program and analyzed its relationship to point mutations, copy number variations, and the meiotic recombination hotspots utilized by males and females. DNA replication timing associated with point mutations far more strongly than predicted from earlier analyses and showed a stronger relationship to transversion than transition mutations. Structural mutations arising from recombination-based mechanisms and recombination hotspots used more extensively by females were enriched in early-replicating parts of the genome, though these relationships appeared to relate more strongly to the genomic distribution of causative sequence features. These results indicate differential and sex-specific relationship of DNA replication timing to different forms of mutation and recombination.     

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.     

Bioenergetics in human evolution and disease: implications for the origins of biological complexity and the missing genetic variation of common diseases

Two major inconsistencies exist in the current neo-Darwinian evolutionary theory that random chromosomal mutations acted on by natural selection generate new species. First, natural selection does not require the evolution of ever increasing complexity, yet this is the hallmark of biology. Second, human chromosomal DNA sequence variation is predominantly either neutral or deleterious and is insufficient to provide the variation required for speciation or for predilection to common diseases. Complexity is explained by the continuous flow of energy through the biosphere that drives the accumulation of nucleic acids and information. Information then encodes complex forms. In animals, energy flow is primarily mediated by mitochondria whose maternally inherited mitochondrial DNA (mtDNA) codes for key genes for energy metabolism. In mammals, the mtDNA has a very high mutation rate, but the deleterious mutations are removed by an ovarian selection system. Hence, new mutations that subtly alter energy metabolism are continuously introduced into the species, permitting adaptation to regional differences in energy environments. Therefore, the most phenotypically significant gene variants arise in the mtDNA, are regional, and permit animals to occupy peripheral energy environments where rarer nuclear DNA (nDNA) variants can accumulate, leading to speciation. The neutralist–selectionist debate is then a consequence of mammals having two different evolutionary strategies: a fast mtDNA strategy for intra-specific radiation and a slow nDNA strategy for speciation. Furthermore, the missing genetic variation for common human diseases is primarily mtDNA variation plus regional nDNA variants, both of which have been missed by large, inter-population association studies.     

Environmental heterogeneity enhances clonal interference

Clonal interference (CI) is a phenomenon that may be important in several asexual microbes. It occurs when population sizes are large and mutation rates to new beneficial alleles are of significant magnitude. Here we explore the role of gene flow and spatial heterogeneity in selection strength in the adaptation of asexuals. We consider a subdivided population of individuals that are adapting, through new beneficial mutations, and that migrate between different patches. The fitness effect of each mutation depends on the patch and all mutations considered are assumed to be unconditionally beneficial. We find that spatial variation in selection pressure affects the rate of adaptive evolution and its qualitative effects depend on the level of gene flow. In particular, we find that both low migration and high levels of heterogeneity lead to enhanced CI. In contrast, for high levels of migration the rate of fixation of adaptive mutations is higher when environmental heterogeneity is present. In addition, we observe that the level of fitness variation is higher and simultaneous fixation of multiple mutations tends to occur in the regime of low migration rates and high heterogeneity.     

Clusters of nucleotide substitutions and insertion/deletion mutations are associated with repeat sequences

The genome-sequencing gold rush has facilitated the use of comparative genomics to uncover patterns of genome evolution, although their causal mechanisms remain elusive. One such trend, ubiquitous to prokarya and eukarya, is the association of insertion/deletion mutations (indels) with increases in the nucleotide substitution rate extending over hundreds of base pairs. The prevailing hypothesis is that indels are themselves mutagenic agents. Here, we employ population genomics data from Escherichia coli, Saccharomyces paradoxus, and Drosophila to provide evidence suggesting that it is not the indels per se but the sequence in which indels occur that causes the accumulation of nucleotide substitutions. We found that about two-thirds of indels are closely associated with repeat sequences and that repeat sequence abundance could be used to identify regions of elevated sequence diversity, independently of indels. Moreover, the mutational signature of indel-proximal nucleotide substitutions matches that of error-prone DNA polymerases. We propose that repeat sequences promote an increased probability of replication fork arrest, causing the persistent recruitment of error-prone DNA polymerases to specific sequence regions over evolutionary time scales. Experimental measures of the mutation rates of engineered DNA sequences and analyses of experimentally obtained collections of spontaneous mutations provide molecular evidence supporting our hypothesis. This study uncovers a new role for repeat sequences in genome evolution and provides an explanation of how fine-scale sequence contextual effects influence mutation rates and thereby evolution.     

Substitution Rate Heterogeneity and the Male Mutation Bias

Germline mutation rates have been found to be higher in males than in females in many organisms, a likely consequence of cell division being more frequent in spermatogenesis than in oogenesis. If the majority of mutations are due to DNA replication error, the male-to-female mutation rate ratio (α_m) is expected to be similar to the ratio of the number of germ line cell divisions in males and females (c), an assumption that can be tested with proper estimates of α_m and c. α_m is usually estimated by comparing substitution rates in putatively neutral sequences on the sex chromosomes. However, substantial regional variation in substitution rates across chromosomes may bias estimates of α_m based on the substitution rates of short sequences. To investigate regional substitution rate variation, we estimated sequence divergence in 16 gametologous introns located on the Z and W chromosomes of five bird species of the order Galliformes. Intron ends and potentially conserved blocks were excluded to reduce the effect of using sequences subject to negative selection. We found significant substitution rate variation within Z chromosome (G₁₅ = 37.6, p = 0.0010) as well as within W chromosome introns (G₁₅ = 44.0, p = 0.0001). This heterogeneity also affected the estimates of α_m, which varied significantly, from 1.53 to 3.51, among the introns (ANOVA: F_13,14 =2.68, p = 0.04). Our results suggest the importance of using extensive data sets from several genomic regions to avoid the effects of regional mutation rate variation and to ensure accurate estimates of α_m. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. [Reviewing Editor: Mr. Martin Kreitman] Nick G.C. Smith Deceased     

The role of background selection in shaping patterns of molecular evolution and variation: evidence from variability on the Drosophila X chromosome

In the putatively ancestral population of Drosophila melanogaster, the ratio of silent DNA sequence diversity for X-linked loci to that for autosomal loci is approximately one, instead of the expected "null" value of 3/4. One possible explanation is that background selection (the hitchhiking effect of deleterious mutations) is more effective on the autosomes than on the X chromosome, because of the lack of crossing over in male Drosophila. The expected effects of background selection on neutral variability at sites in the middle of an X chromosome or an autosomal arm were calculated for different models of chromosome organization and methods of approximation, using current estimates of the deleterious mutation rate and distributions of the fitness effects of deleterious mutations. The robustness of the results to different distributions of fitness effects, dominance coefficients, mutation rates, mapping functions, and chromosome size was investigated. The predicted ratio of X-linked to autosomal variability is relatively insensitive to these variables, except for the mutation rate and map length. Provided that the deleterious mutation rate per genome is sufficiently large, it seems likely that background selection can account for the observed X to autosome ratio of variability in the ancestral population of D. melanogaster. The fact that this ratio is much less than one in D. pseudoobscura is also consistent with the model's predictions, since this species has a high rate of crossing over. The results suggest that background selection may play a major role in shaping patterns of molecular evolution and variation.