Single Nucleotide Mutation Rate Research Articles

Global and local genomic features together modulate the spontaneous single nucleotide mutation rate

Spontaneous mutations are evolutionary engines as they generate variants for the evolutionary downstream processes that give rise to speciation and adaptation. Single nucleotide mutations (SNM) are the most abundant type of mutations among them. Here, we perform a meta-analysis to quantify the influence of selected global genomic parameters (genome size, genomic GC content, genomic repeat fraction, number of coding genes, gene count, and strand bias in prokaryotes) and local genomic features (local GC content, repeat content, CpG content and the number of SNM at CpG islands) on spontaneous SNM rates across the tree of life (prokaryotes, unicellular eukaryotes, multicellular eukaryotes) using wild-type sequence data in two different taxon classification systems. We find that the spontaneous SNM rates in our data are correlated with many genomic features in prokaryotes and unicellular eukaryotes irrespective of their sample sizes. On the other hand, only the number of coding genes was correlated with the spontaneous SNM rates in multicellular eukaryotes primarily contributed by vertebrates data. Considering local features, we notice that local GC content and CpG content significantly were correlated with the spontaneous SNM rates in the unicellular eukaryotes, while local repeat fraction is an important feature in prokaryotes and certain specific uni- and multi-cellular eukaryotes. Such predictive features of the spontaneous SNM rates often support non-linear models as the best fit compared to the linear model. We also observe that the strand asymmetry in prokaryotes plays an important role in determining the spontaneous SNM rates but the SNM spectrum does not.

Structural underpinnings of mutation rate variations in the human genome.

Single nucleotide mutation rates have critical implications for human evolution and genetic diseases. Importantly, the rates vary substantially across the genome and the principles underlying such variations remain poorly understood. A recent model explained much of this variation by considering higher-order nucleotide interactions in the 7-mer sequence context around mutated nucleotides. This model's success implicates a connection between DNA shape and mutation rates. DNA shape, i.e. structural properties like helical twist and tilt, is known to capture interactions between nucleotides within a local context. Thus, we hypothesized that changes in DNA shape features at and around mutated positions can explain mutation rate variations in the human genome. Indeed, DNA shape-based models of mutation rates showed similar or improved performance over current nucleotide sequence-based models. These models accurately characterized mutation hotspots in the human genome and revealed the shape features whose interactions underlie mutation rate variations. DNA shape also impacts mutation rates within putative functional regions like transcription factor binding sites where we find a strong association between DNA shape and position-specific mutation rates. This work demonstrates the structural underpinnings of nucleotide mutations in the human genome and lays the groundwork for future models of genetic variations to incorporate DNA shape.

Salt stress alters the spectrum of de novo mutation available to selection during experimental adaptation of Chlamydomonas reinhardtii.

The genetic basis of adaptation is driven by both selection and the spectrum of available mutations. Given that the rate of mutation is not uniformly distributed across the genome and varies depending on the environment, understanding the signatures of selection across the genome is aided by first establishing what the expectations of genetic change are from mutation. To determine the interaction between salt stress, selection, and mutation across the genome, we compared mutations observed in a selection experiment for salt tolerance in Chlamydomonas reinhardtii to those observed in mutation accumulation (MA) experiments with and without salt exposure. MA lines evolved under salt stress had a single-nucleotide mutation rate of , similar to that of MA lines under standard conditions ( ). However, we found that salt stress led to an increased rate of indel mutations, but that many of these mutations were removed under selection. Finally, lines adapted to salt also showed excess clustering of mutations in the genome and the co-expression network, suggesting a role for positive selection in retaining mutations in particular compartments of the genome during the evolution of salt tolerance. Our study shows that characterizing mutation rates and spectra expected under stress helps disentangle the effects of environment and selection duringadaptation.

Spontaneous rate of clonal single nucleotide mutations in Daphnia galeata.

Mutations are the ultimate source of heritable variation and therefore the fuel for evolution, but direct estimates of mutation rates exist only for few species. We estimated the spontaneous single nucleotide mutation rate among clonal generations in the waterflea Daphnia galeata with a short-term mutation accumulation approach. Individuals from eighteen mutation accumulation lines over five generations were deep sequenced to count de novo mutations that were not present in a pool of F1 individuals, representing the parental genotype. We identified 12 new nucleotide mutations in 90 clonal generational passages. This resulted in an estimated single nucleotide mutation rate of 0.745 x 10−9 (95% c.f. 0.39 x 10−9–1.26 x 10−9), which is slightly lower than recent estimates for other Daphnia species. We discuss the implications for the population genetics of Cladocerans.

Evolutionary Dynamics of Asexual Hypermutators Adapting to a Novel Environment.

How microbes adapt to a novel environment is a central question in evolutionary biology. Although adaptive evolution must be fueled by beneficial mutations, whether higher mutation rates facilitate the rate of adaptive evolution remains unclear. To address this question, we cultured Escherichia coli hypermutating populations, in which a defective methyl-directed mismatch repair pathway causes a 140-fold increase in single-nucleotide mutation rates. In parallel with wild-type E. coli, populations were cultured in tubes containing Luria-Bertani broth, a complex medium known to promote the evolution of subpopulation structure. After 900 days of evolution, in three transfer schemes with different population-size bottlenecks, hypermutators always exhibited similar levels of improved fitness as controls. Fluctuation tests revealed that the mutation rates of hypermutator lines converged evolutionarily on those of wild-type populations, which may have contributed to the absence of fitness differences. Further genome-sequence analysis revealed that, although hypermutator populations have higher rates of genomic evolution, this largely reflects strong genetic linkage. Despite these linkage effects, the evolved population exhibits parallelism in fixed mutations, including those potentially related to biofilm formation, transcription regulation, and mutation-rate evolution. Together, these results are generally inconsistent with a hypothesized positive relationship between the mutation rate and the adaptive speed of evolution, and provide insight into how clonal adaptation occurs in novel environments.

De Novo Mutation Rate Variation and Its Determinants in Chlamydomonas.

De novo mutations are central for evolution, since they provide the raw material for natural selection by regenerating genetic variation. However, studying de novo mutations is challenging and is generally restricted to model species, so we have a limited understanding of the evolution of the mutation rate and spectrum between closely related species. Here, we present a mutation accumulation (MA) experiment to study de novo mutation in the unicellular green alga Chlamydomonas incerta and perform comparative analyses with its closest known relative, Chlamydomonas reinhardtii. Using whole-genome sequencing data, we estimate that the median single nucleotide mutation (SNM) rate in C. incerta is μ = 7.6 × 10−10, and is highly variable between MA lines, ranging from μ = 0.35 × 10−10 to μ = 131.7 × 10−10. The SNM rate is strongly positively correlated with the mutation rate for insertions and deletions between lines (r > 0.97). We infer that the genomic factors associated with variation in the mutation rate are similar to those in C. reinhardtii, allowing for cross-prediction between species. Among these genomic factors, sequence context and complexity are more important than GC content. With the exception of a remarkably high C→T bias, the SNM spectrum differs markedly between the two Chlamydomonas species. Our results suggest that similar genomic and biological characteristics may result in a similar mutation rate in the two species, whereas the SNM spectrum has more freedom to diverge.

Loss of Heterozygosity and Base Mutation Rates Vary Among Saccharomyces cerevisiae Hybrid Strains.

A growing body of evidence suggests that mutation rates exhibit intra-species specific variation. We estimated genome-wide loss of heterozygosity (LOH), gross chromosomal changes, and single nucleotide mutation rates to determine intra-species specific differences in hybrid and homozygous strains of Saccharomyces cerevisiae. The mutation accumulation lines of the S. cerevisiae hybrid backgrounds - S288c/YJM789 (S/Y) and S288c/RM11-1a (S/R) were analyzed along with the homozygous diploids RM11, S288c, and YJM145. LOH was extensive in both S/Y and S/R hybrid backgrounds. The S/Y background also showed longer LOH tracts, gross chromosomal changes, and aneuploidy. Short copy number aberrations were observed in the S/R background. LOH data from the S/Y and S/R hybrids were used to construct a LOH map for S288c to identify hotspots. Further, we observe up to a sixfold difference in single nucleotide mutation rates among the S. cerevisiae S/Y and S/R genetic backgrounds. Our results demonstrate LOH is common during mitotic divisions in S. cerevisiae hybrids and also highlight genome-wide differences in LOH patterns and rates of single nucleotide mutations between commonly used S. cerevisiae hybrid genetic backgrounds.

Evolution of Listeria monocytogenes in a Food Processing Plant Involves Limited Single-Nucleotide Substitutions but Considerable Diversification by Gain and Loss of Prophages

Whole-genome sequencing (WGS) is becoming the standard method for subtyping Listeria monocytogenes Interpretation of WGS data for isolates from foods and associated environments is, however, challenging due to a lack of detailed data on Listeria evolution in processing facilities. Here, we used previously collected WGS data for 40 L. monocytogenes isolates obtained from a cold-smoked salmon processing facility between 1998 and 2015 to probe the L. monocytogenes molecular evolution in this facility, combined with phenotypic assessment of selected isolates. Isolates represented three clusters (1, 2, and 3); cluster 3 isolates (n = 32) were obtained over 18 years. The average mutation rate for cluster 3 was estimated as 1.15 × 10-7 changes per nucleotide per year (∼0.35 changes per genome per year); the most recent common ancestors (MRCAs) of subclusters 3a and 3b were estimated to have occurred around 1958 and 1974, respectively, within the age of the facility, suggesting long-term persistence in this facility. Extensive prophage diversity was observed within subclusters 3a and 3b, which have one shared and six unique prophage profiles for each subcluster (with 16 prophage profiles found among all 40 isolates). The plasmid-borne sanitizer tolerance operon bcrABC was found in all cluster 2 and 3 isolates, while the transposon-borne sanitizer tolerance gene qacH was found in one cluster 1 isolate; presence of these genes was correlated with the ability to survive increased concentrations of sanitizers. Selected isolates showed significant variation in the ability to attach to surfaces, with persistent isolates attaching better than transient isolates at 21°C.IMPORTANCE Knowledge about the genetic evolution of L. monocytogenes in food processing facilities over multiple years is generally lacking. This information is critical to interpret WGS findings involving food or food-associated isolates. This study suggests that L. monocytogenes that persists in processing facilities may evolve with a low single-nucleotide mutation rate mostly driven by negative (i.e., purifying) selection but with rapid diversification of prophages. Hence, isolation of L. monocytogenes with few single-nucleotide polymorphism (SNP) differences in different locations (e.g., supplier plants and receiving plants) is possible, highlighting the importance of epidemiological and detailed isolate metadata for interpreting WGS data in traceback investigation. Our study also shows how advanced WGS data analyses can be used to support root cause analysis efforts and may, for example, pinpoint the time when a persistence event started (which then potentially could be linked to facility changes, introduction of new equipment, etc.).

Prophage-Driven Genomic Structural Changes Promote Bartonella Vertical Evolution.

Bartonella is a genetically diverse group of vector-borne bacteria. Over 40 species have been characterized to date, mainly from mammalian reservoirs and arthropod vectors. Rodent reservoirs harbor one of the largest Bartonella diversity described to date, and novel species and genetic variants are continuously identified from these hosts. Yet, it is still unknown if this significant genetic diversity stems from adaptation to different niches or from intrinsic high mutation rates. Here, we explored the vertical occurrence of spontaneous genomic alterations in 18 lines derived from two rodent-associated Bartonella elizabethae-like strains, evolved in nonselective agar plates under conditions mimicking their vector- and mammalian-associated temperatures, and the transmission cycles between them (i.e., 26 °C, 37 °C, and alterations between the two), using mutation accumulation experiments. After ∼1,000 generations, evolved genomes revealed few point mutations (average of one-point mutation per line), evidencing conserved single-nucleotide mutation rates. Interestingly, three large structural genomic changes (two large deletions and an inversion) were identified over all lines, associated with prophages and surface adhesin genes. Particularly, a prophage, deleted during constant propagation at 37 °C, was associated with an increased autonomous replication at 26 °C (the flea-associated temperature). Complementary molecular analyses of wild strains, isolated from desert rodents and their fleas, further supported the occurrence of structural genomic variations and prophage-associated deletions in nature. Our findings suggest that structural genomic changes represent an effective intrinsic mechanism to generate diversity in slow-growing bacteria and emphasize the role of prophages as promoters of diversity in nature.

Small Insertions Are More Deleterious than Small Deletions in Human Genomes

Although lines of evidence suggest that small insertions and deletions differ in their mechanisms of formation, there remains the debate on whether natural selection acts differently on the two indel types. Currently available personal genomes and the 1000 Genomes Project permit population level and genome scale comparison of the selection regimes on the two indel types. We first developed a statistical model to evaluate the indel frequency biases of the 1000 Genomes Project phase 1 data. We then identified four independent lines of evidence demonstrating that human small (1-4bp) insertions are on average more deleterious than deletions. This genome-wide selection pattern is not affected by methodology, demography, and regional differences including indel density, introns versus exons, repeats versus nonrepeats, recombination rates, and the timing of DNA replication. This selection pattern has a profound effect on indel frequency spectra, deletional bias, and local single-nucleotide mutation rates. Finally, we observed that small insertions appear to be more actively implicated in shaping fast-evolving genomic sequences (or nonconserved regions).

Fine-Scale Signatures of Molecular Evolution Reconcile Models of Indel-Associated Mutation

Genomic structural alterations that vary within species, known as large copy number variants, represent an unanticipated and abundant source of genetic diversity that associates with variation in gene expression and susceptibility to disease. Even short insertions and deletions (indels) can exert important effects on genomes by locally increasing the mutation rate, with multiple mechanisms proposed to account for this pattern. To better understand how indels promote genome evolution, we demonstrate that the single nucleotide mutation rate is elevated in the vicinity of indels, with a resolution of tens of base pairs, for the two closely related nematode species Caenorhabditis remanei and C. sp. 23. In addition to indels being clustered with single nucleotide polymorphisms and fixed differences, we also show that transversion mutations are enriched in sequences that flank indels and that many indels associate with sequence repeats. These observations are compatible with a model that reconciles previously proposed mechanisms of indel-associated mutagenesis, implicating repeat sequences as a common driver of indel errors, which then recruit error-prone polymerases during DNA repair, resulting in a locally elevated single nucleotide mutation rate. The striking influence of indel variants on the molecular evolution of flanking sequences strengthens the emerging general view that mutations can induce further mutations.

Increased paternal age and the influence on burden of genomic copy number variation in the general population

Genomic copy number variations (CNVs) and increased parental age are both associated with the risk to develop a variety of clinical neuropsychiatric disorders such as autism, schizophrenia and bipolar disorder. At the same time, it has been shown that the rate of transmitted de novo single nucleotide mutations is increased with paternal age. To address whether paternal age also affects the burden of structural genomic deletions and duplications, we examined various types of CNV burden in a large population sample from the Netherlands. Healthy participants with parental age information (n=6,773) were collected at different University Medical Centers. CNVs were called with the PennCNV algorithm using Illumina genome-wide SNP array data. We observed no evidence in support of a paternal age effect on CNV load in the offspring. Our results were negative for global measures as well as several proxies for de novo CNV events in this unique sample. While recent studies suggest de novo single nucleotide mutation rate to be dominated by the age of the father at conception, our results strongly suggest that at the level of global CNV burden there is no influence of increased paternal age. While it remains possible that local genomic effects may exist for specific phenotypes, this study indicates that global CNV burden and increased father's age may be independent disease risk factors.

Longer time scale for human evolution

During the past few years, the best estimates of the human single nucleotide mutation rate have been cut in half. Until recently, estimates of mutation rate have relied on counting substitutions between primate species and assuming that fossil relatives of living species can accurately pin dates onto phylogenetic branches. This procedure allows very precise estimates, but introduces systematic bias toward higher substitution rates and longer branch lengths because a new lineage can leave a fossil record only after its origin, never beforehand (1). Now, widespread resequencing, initially of de novo Mendelian genetic disorders (2, 3) and later of whole genomes in parent–offspring trios (4), has allowed direct comparisons of parent and offspring genomes. The most commonly used, but now outdated, estimate of the mutation rate was 2.4 × 10−8 changes per nucleotide per generation (5). Current estimates of the same value based on resequencing data are much lower, approximately 1.1 to 1.28 × 10−8 (3, 4).

More single-nucleotide mutations surround small insertions than small deletions in primates

Early studies have shown that single-nucleotide mutation rates increase close to insertions and deletions, but it is not fully understood how natural selection shapes genome-wide patterns of indels and their nearby single-nucleotide mutations. In this study, we find that, in primates, more single-nucleotide mutations surround small insertions than small deletions. This pattern affects <150 base pair (bp) sequences close to indels and persists under different genomic properties, such as exon/intron/intergenic contexts, repeated/nonrepeated sequences, replication timing, recombination rates, indel density, and guanine-cytosine (GC) content. We propose two different, but not mutually exclusive, hypothetical mechanisms to explain the pattern. One mechanism is that the sequence context preferring insertion formation may also favor nucleotide substitutions. Another mechanism is related to a hypothesis in which indel heterozygosity tends to increase nearby nucleotide substitution rates. It means that if insertions spend more time in heterozygotes, insertions may accumulate more surrounding single-nucleotide changes. In conclusion, we characterize a special genome-wide evolutionary pattern for indels and nearby single-nucleotide changes. This pattern may be driven by natural selection and bias primates' genome evolution and phenotypic variations.

Analysis of the genome sequences of three Drosophila melanogaster spontaneous mutation accumulation lines

We inferred the rate and properties of new spontaneous mutations in Drosophila melanogaster by carrying out whole-genome shotgun sequencing-by-synthesis of three mutation accumulation (MA) lines that had been maintained by close inbreeding for an average of 262 generations. We tested for the presence of new mutations by generating alignments of each MA line to the D. melanogaster reference genome sequence and then compared these alignments base by base. We determined empirically that at least five reads at a site within each line are required for accurate single nucleotide mutation calling. We mapped a total of 174 single-nucleotide mutations, giving a single nucleotide mutation rate of 3.5 x 10(-9) per site per generation. There were no false positives in a random sample of 40 of these mutations checked by Sanger sequencing. Variation in the numbers of mutations among the MA lines was small and nonsignificant. Numbers of transition and transversion mutations were 86 and 88, respectively, implying that transition mutation rate is close to 2x the transversion rate. We observed 1.5x as many G or C --> A or T as A or T --> G or C mutations, implying that the G or C --> A or T mutation rate is close to 2x the A or T --> G or C mutation rate. The base composition of the genome is therefore not at an equilibrium determined solely by mutation. The predicted G + C content at mutational equilibrium (33%) is similar to that observed in transposable element remnants. Nearest-neighbor mutational context dependencies are nonsignificant, suggesting that this is a weak phenomenon in Drosophila. We also saw nonsignificant differences in the mutation rate between transcribed and untranscribed regions, implying that any transcription-coupled repair process is weak. Of seven short indel mutations confirmed, six were deletions, consistent with the deletion bias that is thought to exist in Drosophila.

Direct Estimation of the Mitochondrial DNA Mutation Rate in Drosophila melanogaster

Mitochondrial DNA (mtDNA) variants are widely used in evolutionary genetics as markers for population history and to estimate divergence times among taxa. Inferences of species history are generally based on phylogenetic comparisons, which assume that molecular evolution is clock-like. Between-species comparisons have also been used to estimate the mutation rate, using sites that are thought to evolve neutrally. We directly estimated the mtDNA mutation rate by scanning the mitochondrial genome of Drosophila melanogaster lines that had undergone approximately 200 generations of spontaneous mutation accumulation (MA). We detected a total of 28 point mutations and eight insertion-deletion (indel) mutations, yielding an estimate for the single-nucleotide mutation rate of 6.2 × 10−8 per site per fly generation. Most mutations were heteroplasmic within a line, and their frequency distribution suggests that the effective number of mitochondrial genomes transmitted per female per generation is about 30. We observed repeated occurrences of some indel mutations, suggesting that indel mutational hotspots are common. Among the point mutations, there is a large excess of G→A mutations on the major strand (the sense strand for the majority of mitochondrial genes). These mutations tend to occur at nonsynonymous sites of protein-coding genes, and they are expected to be deleterious, so do not become fixed between species. The overall mtDNA mutation rate per base pair per fly generation in Drosophila is estimated to be about 10× higher than the nuclear mutation rate, but the mitochondrial major strand G→A mutation rate is about 70× higher than the nuclear rate. Silent sites are substantially more strongly biased towards A and T than nonsynonymous sites, consistent with the extreme mutation bias towards A+T. Strand-asymmetric mutation bias, coupled with selection to maintain specific nonsynonymous bases, therefore provides an explanation for the extreme base composition of the mitochondrial genome of Drosophila.

Single-nucleotide mutation rate increases close to insertions/deletions in eukaryotes

Mutation hotspots are commonly observed in genomic sequences and certain human disease loci, but general mechanisms for their formation remain elusive. Here we investigate the distribution of single-nucleotide changes around insertions/deletions (indels) in six independent genome comparisons, including primates, rodents, fruitfly, rice and yeast. In each of these genomic comparisons, nucleotide divergence (D) is substantially elevated surrounding indels and decreases monotonically to near-background levels over several hundred bases. D is significantly correlated with both size and abundance of nearby indels. In comparisons of closely related species, derived nucleotide substitutions surrounding indels occur in significantly greater numbers in the lineage containing the indel than in the one containing the ancestral (non-indel) allele; the same holds within species for single-nucleotide mutations surrounding polymorphic indels. We propose that heterozygosity for an indel is mutagenic to surrounding sequences, and use yeast genome-wide polymorphism data to estimate the increase in mutation rate. The consistency of these patterns within and between species suggests that indel-associated substitution is a general mutational mechanism.