CDMAP/CDVIS: context-dependent mutation analysis package and visualization software.

The Context-dependent Mutation Analysis Package and Visualization Software (CDMAP/CDVIS) is an automated, modular toolkit used for the analysis and visualization of context-dependent mutation patterns (site-specific variation in mutation rate from neighboring-nucleotide effects). The CDMAP computes context-dependent mutation rates using a Variant Call File (VCF), Genbank file, and reference genome and can generate high-resolution figures to analyze variation in mutation rate across spatiotemporal scales. This algorithm has been benchmarked against mutation accumulation data but can also be used to calculate context-dependent mutation rates for polymorphism or closely related species as long as the input requirements are met. Output from CDMAP can be integrated into CDVIS, an interactive database for visualizing mutation patterns across multiple taxa simultaneously.

Precise measurement of the fitness effects of spontaneous mutations by droplet digital PCR in Burkholderia cenocepacia.

Understanding how mutations affect survivability is a key component to knowing how organisms and complex traits evolve. However, most mutations have a minor effect on fitness and these effects are difficult to resolve using traditional molecular techniques. Therefore, there is a dire need for more accurate and precise fitness measurements methods. Here, we measured the fitness effects in Burkholderia cenocepacia HI2424 mutation accumulation (MA) lines using droplet-digital polymerase chain reaction (ddPCR). Overall, the fitness measurements from ddPCR-MA are correlated positively with fitness measurements derived from traditional phenotypic marker assays (r = 0.297, P = 0.05), but showed some differences. First, ddPCR had significantly lower measurement variance in fitness (F = 3.78, P < 2.6 × 10-13) in control experiments. Second, the mean fitness from ddPCR-MA measurements were significantly lower than phenotypic marker assays (-0.0041 vs -0.0071, P = 0.006). Consistent with phenotypic marker assays, ddPCR-MA measurements observed multiple (27/43) lineages that significantly deviated from mean fitness, suggesting that a majority of the mutations are neutral or slightly deleterious and intermixed with a few mutations that have extremely large effects. Of these mutations, we found a significant excess of mutations within DNA excinuclease and Lys R transcriptional regulators that have extreme deleterious and beneficial effects, indicating that modifications to transcription and replication may have a strong effect on organismal fitness. This study demonstrates the power of ddPCR as a ubiquitous method for high-throughput fitness measurements in both DNA- and RNA-based organisms regardless of cell type or physiology.

Unexpected Discovery of Hypermutator Phenotype Sounds the Alarm for Quality Control Strains.

Microbial strains with high genomic stability are particularly sought after for testing the quality of commercial microbiological products, such as biological media and antibiotics. Yet, using mutation-accumulation experiments and de novo assembled complete genomes based on Nanopore long-read sequencing, we find that the widely used quality-control strain Shewanella putrefaciens ATCC-8071, also a facultative pathogen, is a hypermutator, with a base-pair substitution mutation rate of 2.42 × 10-8 per nucleotide site per cell division, ∼146-fold greater than that of the wild-type strain CGMCC-1.6515. Using complementation experiments, we confirm that mutL dysfunction, which was a recent evolutionary event, is the cause for the high mutation rate of ATCC-8071. Further analyses also give insight into possible relationships between mutation and genome evolution in this important bacterium. This discovery of a well-known strain being a hypermutator necessitates screening the mutation rate of bacterial strains before any quality control or experiments.

The insect-killing bacterium Photorhabdus luminescens has the lowest mutation rate among bacteria.

Mutation is a primary source of genetic variation that is used to power evolution. Many studies, however, have shown that most mutations are deleterious and, as a result, extremely low mutation rates might be beneficial for survival. Using a mutation accumulation experiment, an unbiased method for mutation study, we found an extremely low base-substitution mutation rate of 5.94 × 10-11 per nucleotide site per cell division (95% Poisson confidence intervals: 4.65 × 10-11, 7.48 × 10-11) and indel mutation rate of 8.25 × 10-12 per site per cell division (95% confidence intervals: 3.96 × 10-12, 1.52 × 10-11) in the bacterium Photorhabdus luminescens ATCC29999. The mutations are strongly A/T-biased with a mutation bias of 10.28 in the A/T direction. It has been hypothesized that the ability for selection to lower mutation rates is inversely proportional to the effective population size (drift-barrier hypothesis) and we found that the effective population size of this bacterium is significantly greater than most other bacteria. This finding further decreases the lower-bounds of bacterial mutation rates and provides evidence that extreme levels of replication fidelity can evolve within organisms that maintain large effective population sizes.

Low Base-Substitution Mutation Rate but High Rate of Slippage Mutations in the Sequence Repeat-Rich Genome of Dictyostelium discoideum.

We describe the rate and spectrum of spontaneous mutations for the social amoeba Dictyostelium discoideum, a key model organism in molecular, cellular, evolutionary and developmental biology. Whole-genome sequencing of 37 mutation accumulation lines of D. discoideum after an average of 1,500 cell divisions yields a base-substitution mutation rate of 2.47 × 10-11 per site per generation, substantially lower than that of most eukaryotic and prokaryotic organisms, and of the same order of magnitude as in the ciliates Paramecium tetraurelia and Tetrahymena thermophila Known for its high genomic AT content and abundance of simple sequence repeats, we observe that base-substitution mutations in D. discoideum are highly A/T biased. This bias likely contributes both to the high genomic AT content and to the formation of simple sequence repeats in the AT-rich genome of Dictyostelium discoideum In contrast to the situation in other surveyed unicellular eukaryotes, indel rates far exceed the base-substitution mutation rate in this organism with a high proportion of 3n indels, particularly in regions without simple sequence repeats. Like ciliates, D. discoideum has a large effective population size, reducing the power of random genetic drift, magnifying the effect of selection on replication fidelity, in principle allowing D. discoideum to evolve an extremely low base-substitution mutation rate.

Periodic Variation of Mutation Rates in Bacterial Genomes Associated with Replication Timing.

The causes and consequences of spatiotemporal variation in mutation rates remain to be explored in nearly all organisms. Here we examine relationships between local mutation rates and replication timing in three bacterial species whose genomes have multiple chromosomes: Vibrio fischeri, Vibrio cholerae, and Burkholderia cenocepacia Following five mutation accumulation experiments with these bacteria conducted in the near absence of natural selection, the genomes of clones from each lineage were sequenced and analyzed to identify variation in mutation rates and spectra. In lineages lacking mismatch repair, base substitution mutation rates vary in a mirrored wave-like pattern on opposing replichores of the large chromosomes of V. fischeri and V. cholerae, where concurrently replicated regions experience similar base substitution mutation rates. The base substitution mutation rates on the small chromosome are less variable in both species but occur at similar rates to those in the concurrently replicated regions of the large chromosome. Neither nucleotide composition nor frequency of nucleotide motifs differed among regions experiencing high and low base substitution rates, which along with the inferred ~800-kb wave period suggests that the source of the periodicity is not sequence specific but rather a systematic process related to the cell cycle. These results support the notion that base substitution mutation rates are likely to vary systematically across many bacterial genomes, which exposes certain genes to elevated deleterious mutational load.IMPORTANCE That mutation rates vary within bacterial genomes is well known, but the detailed study of these biases has been made possible only recently with contemporary sequencing methods. We applied these methods to understand how bacterial genomes with multiple chromosomes, like those of Vibrio and Burkholderia, might experience heterogeneous mutation rates because of their unusual replication and the greater genetic diversity found on smaller chromosomes. This study captured thousands of mutations and revealed wave-like rate variation that is synchronized with replication timing and not explained by sequence context. The scale of this rate variation over hundreds of kilobases of DNA strongly suggests that a temporally regulated cellular process may generate wave-like variation in mutation risk. These findings add to our understanding of how mutation risk is distributed across bacterial and likely also eukaryotic genomes, owing to their highly conserved replication and repair machinery.

An Unbiased Genome-Wide View of the Mutation Rate and Spectrum of the Endosymbiotic Bacterium Teredinibacter turnerae.

Mutations contribute to genetic variation in all living systems. Thus, precise estimates of mutation rates and spectra across a diversity of organisms are required for a full comprehension of evolution. Here, a mutation-accumulation (MA) assay was carried out on the endosymbiotic bacterium Teredinibacter turnerae. After ∼3,025 generations, base-pair substitutions (BPSs) and insertion-deletion (indel) events were characterized by whole-genome sequencing analysis of 47 independent MA lines, yielding a BPS rate of 1.14 × 10-9 per site per generation and indel rate of 1.55 × 10-10 events per site per generation, which are among the highest within free-living and facultative intracellular bacteria. As in other endosymbionts, a significant bias of BPSs toward A/T and an excess of deletion mutations over insertion mutations are observed for these MA lines. However, even with a deletion bias, the genome remains relatively large (∼5.2 Mb) for an endosymbiotic bacterium. The estimate of the effective population size (Ne) in T. turnerae is quite high and comparable to free-living bacteria (∼4.5 × 107), suggesting that the heavy bottlenecking associated with many endosymbiotic relationships is not prevalent during the life of this endosymbiont. The efficiency of selection scales with increasing Ne and such strong selection may have been operating against the deletion bias, preventing genome erosion. The observed mutation rate in this endosymbiont is of the same order of magnitude of those with similar Ne, consistent with the idea that population size is a primary determinant of mutation-rate evolution within endosymbionts, and that not all endosymbionts have low Ne.

Evolutionary determinants of genome-wide nucleotide composition.

One of the long-standing mysteries of evolutionary genomics is the source of the wide phylogenetic diversity in genome nucleotide composition (G + C versus A + T), which must be a consequence of interspecific differences in mutation bias, the efficiency of selection for different nucleotides or a combination of the two. We demonstrate that although genomic G + C composition is strongly driven by mutation bias, it is also substantially modified by direct selection and/or as a by-product of biased gene conversion. Moreover, G + C composition at fourfold redundant sites is consistently elevated above the neutral expectation-more so than for any other class of sites.

Single-sperm sequencing reveals the accelerated mitochondrial mutation rate in male Daphnia pulex (Crustacea, Cladocera).

Mutation rate in the nuclear genome differs between sexes, with males contributing more mutations than females to their offspring. The male-biased mutation rates in the nuclear genome is most likely to be driven by a higher number of cell divisions in spermatogenesis than in oogenesis, generating more opportunities for DNA replication errors. However, it remains unknown whether male-biased mutation rates are present in mitochondrial DNA (mtDNA). Although mtDNA is maternally inherited and male mtDNA mutation typically does not contribute to genetic variation in offspring, male mtDNA mutations are critical for male reproductive health. In this study, we measured male mtDNA mutation rate using publicly available whole-genome sequences of single sperm of the freshwater microcrustacean Daphnia pulex Using a stringent mutation detection pipeline, we found that the male mtDNA mutation rate is 3.32 × 10-6 per site per generation. All the detected mutations are heteroplasmic base substitutions, with 57% of mutations converting G/C to A/T nucleotides. Consistent with the male-biased mutation in the nuclear genome, the male mtDNA mutation rate in D. pulex is approximately 20 times higher than the female rate per generation. We propose that the elevated mutation rate per generation in male mtDNA is consistent with an increased number of cell divisions during male gametogenesis.

Spontaneous mutations of a model heterotrophic marine bacterium.

Heterotrophic marine bacterioplankton populations display substantive genomic diversity that is commonly explained to be the result of selective forces imposed by resource limitation or interactions with phage and predators. Here we use a mutation-accumulation experiment followed by whole-genome sequencing of mutation lines to determine an unbiased rate and molecular spectrum of spontaneous mutations for a model heterotrophic marine bacterium in the globally important Roseobacter clade, Ruegeria pomeroyi DSS-3. We find evidence for mutational bias towards deletions over insertions, and this process alone could account for a sizable portion of genome size diversity among roseobacters and also implies that lateral gene transfer and/or selection must also play a role in maintaining roseobacters with large genome sizes. We also find evidence for a mutational bias in favor of changes from A/T to G/C nucleobases, which explains widespread occurrences of G/C-enriched Roseobacter genomes. Using the calculated mutation rate of 1.39 × 10-10 per base per generation, we implement a 'mutation-rate clock' approach to date the evolution of roseobacters by assuming a constant mutation rate along their evolutionary history. This approach gives an estimated date of Roseobacter genome expansion in good agreement with an earlier fossil-based estimate of ~250 million years ago and is consistent with a hypothesis of a correlated evolutionary history between roseobacters and marine eukaryotic phytoplankton groups.

Genome-Wide Biases in the Rate and Molecular Spectrum of Spontaneous Mutations in Vibrio cholerae and Vibrio fischeri.

The vast diversity in nucleotide composition and architecture among bacterial genomes may be partly explained by inherent biases in the rates and spectra of spontaneous mutations. Bacterial genomes with multiple chromosomes are relatively unusual but some are relevant to human health, none more so than the causative agent of cholera, Vibrio cholerae Here, we present the genome-wide mutation spectra in wild-type and mismatch repair (MMR) defective backgrounds of two Vibrio species, the low-%GC squid symbiont V. fischeri and the pathogen V. cholerae, collected under conditions that greatly minimize the efficiency of natural selection. In apparent contrast to their high diversity in nature, both wild-type V. fischeri and V. cholerae have among the lowest rates for base-substitution mutations (bpsms) and insertion-deletion mutations (indels) that have been measured, below 10-3/genome/generation. Vibrio fischeri and V. cholerae have distinct mutation spectra, but both are AT-biased and produce a surprising number of multi-nucleotide indels. Furthermore, the loss of a functional MMR system caused the mutation spectra of these species to converge, implying that the MMR system itself contributes to species-specific mutation patterns. Bpsm and indel rates varied among genome regions, but do not explain the more rapid evolutionary rates of genes on chromosome 2, which likely result from weaker purifying selection. More generally, the very low mutation rates of Vibrio species correlate inversely with their immense population sizes and suggest that selection may not only have maximized replication fidelity but also optimized other polygenic traits relative to the constraints of genetic drift.

Genetic drift, selection and the evolution of the mutation rate.

As one of the few cellular traits that can be quantified across the tree of life, DNA-replication fidelity provides an excellent platform for understanding fundamental evolutionary processes. Furthermore, because mutation is the ultimate source of all genetic variation, clarifying why mutation rates vary is crucial for understanding all areas of biology. A potentially revealing hypothesis for mutation-rate evolution is that natural selection primarily operates to improve replication fidelity, with the ultimate limits to what can be achieved set by the power of random genetic drift. This drift-barrier hypothesis is consistent with comparative measures of mutation rates, provides a simple explanation for the existence of error-prone polymerases and yields a formal counter-argument to the view that selection fine-tunes gene-specific mutation rates.

Low Base-Substitution Mutation Rate in the Germline Genome of the Ciliate Tetrahymena thermophil.

Mutation is the ultimate source of all genetic variation and is, therefore, central to evolutionary change. Previous work on Paramecium tetraurelia found an unusually low germline base-substitution mutation rate in this ciliate. Here, we tested the generality of this result among ciliates using Tetrahymena thermophila. We sequenced the genomes of 10 lines of T. thermophila that had each undergone approximately 1,000 generations of mutation accumulation (MA). We applied an existing mutation-calling pipeline and developed a new probabilistic mutation detection approach that directly models the design of an MA experiment and accommodates the noise introduced by mismapped reads. Our probabilistic mutation-calling method provides a straightforward way of estimating the number of sites at which a mutation could have been called if one was present, providing the denominator for our mutation rate calculations. From these methods, we find that T. thermophila has a germline base-substitution mutation rate of 7.61 × 10 - 12 per-site, per cell division, which is consistent with the low base-substitution mutation rate in P. tetraurelia. Over the course of the evolution experiment, genomic exclusion lines derived from the MA lines experienced a fitness decline that cannot be accounted for by germline base-substitution mutations alone, suggesting that other genetic or epigenetic factors must be involved. Because selection can only operate to reduce mutation rates based upon the "visible" mutational load, asexual reproduction with a transcriptionally silent germline may allow ciliates to evolve extremely low germline mutation rates.

Evolution of the Insertion-Deletion Mutation Rate Across the Tree of Life.

Mutations are the ultimate source of variation used for evolutionary adaptation, while also being predominantly deleterious and a source of genetic disorders. Understanding the rate of insertion-deletion mutations (indels) is essential to understanding evolutionary processes, especially in coding regions, where such mutations can disrupt production of essential proteins. Using direct estimates of indel rates from 14 phylogenetically diverse eukaryotic and bacterial species, along with measures of standing variation in such species, we obtain results that imply an inverse relationship of mutation rate and effective population size. These results, which corroborate earlier observations on the base-substitution mutation rate, appear most compatible with the hypothesis that natural selection reduces mutation rates per effective genome to the point at which the power of random genetic drift (approximated by the inverse of effective population size) becomes overwhelming. Given the substantial differences in DNA metabolism pathways that give rise to these two types of mutations, this consistency of results raises the possibility that refinement of other molecular and cellular traits may be inversely related to species-specific levels of random genetic drift.

The Rate and Spectrum of Spontaneous Mutations in Mycobacterium smegmatis, a Bacterium Naturally Devoid of the Postreplicative Mismatch Repair Pathway.

Mycobacterium smegmatis is a bacterium that is naturally devoid of known postreplicative DNA mismatch repair (MMR) homologs, mutS and mutL, providing an opportunity to investigate how the mutation rate and spectrum has evolved in the absence of a highly conserved primary repair pathway. Mutation accumulation experiments of M. smegmatis yielded a base-substitution mutation rate of 5.27 × 10(-10) per site per generation, or 0.0036 per genome per generation, which is surprisingly similar to the mutation rate in MMR-functional unicellular organisms. Transitions were found more frequently than transversions, with the A:T→G:C transition rate significantly higher than the G:C→A:T transition rate, opposite to what is observed in most studied bacteria. We also found that the transition-mutation rate of M. smegmatis is significantly lower than that of other naturally MMR-devoid or MMR-knockout organisms. Two possible candidates that could be responsible for maintaining high DNA fidelity in this MMR-deficient organism are the ancestral-like DNA polymerase DnaE1, which contains a highly efficient DNA proofreading histidinol phosphatase (PHP) domain, and/or the existence of a uracil-DNA glycosylase B (UdgB) homolog that might protect the GC-rich M. smegmatis genome against DNA damage arising from oxidation or deamination. Our results suggest that M. smegmatis has a noncanonical Dam (DNA adenine methylase) methylation system, with target motifs differing from those previously reported. The mutation features of M. smegmatis provide further evidence that genomes harbor alternative routes for improving replication fidelity, even in the absence of major repair pathways.

High mutational rates of large-scale duplication and deletion in Daphnia pulex.

Knowledge of the genome-wide rate and spectrum of mutations is necessary to understand the origin of disease and the genetic variation driving all evolutionary processes. Here, we provide a genome-wide analysis of the rate and spectrum of mutations obtained in two Daphnia pulex genotypes via separate mutation-accumulation (MA) experiments. Unlike most MA studies that utilize haploid, homozygous, or self-fertilizing lines, D. pulex can be propagated ameiotically while maintaining a naturally heterozygous, diploid genome, allowing the capture of the full spectrum of genomic changes that arise in a heterozygous state. While base-substitution mutation rates are similar to those in other multicellular eukaryotes (about 4 × 10(-9) per site per generation), we find that the rates of large-scale (>100 kb) de novo copy-number variants (CNVs) are significantly elevated relative to those seen in previous MA studies. The heterozygosity maintained in this experiment allowed for estimates of gene-conversion processes. While most of the conversion tract lengths we report are similar to those generated by meiotic processes, we also find larger tract lengths that are indicative of mitotic processes. Comparison of MA lines to natural isolates reveals that a majority of large-scale CNVs in natural populations are removed by purifying selection. The mutations observed here share similarities with disease-causing, complex, large-scale CNVs, thereby demonstrating that MA studies in D. pulex serve as a system for studying the processes leading to such alterations.

The Spontaneous Mutation Rate in the Fission Yeast Schizosaccharomyces pombe.

The rate at which new mutations arise in the genome is a key factor in the evolution and adaptation of species. Here we describe the rate and spectrum of spontaneous mutations for the fission yeast Schizosaccharomyces pombe, a key model organism with many similarities to higher eukaryotes. We undertook an ∼1700-generation mutation accumulation (MA) experiment with a haploid S. pombe, generating 422 single-base substitutions and 119 insertion-deletion mutations (indels) across the 96 replicates. This equates to a base-substitution mutation rate of 2.00 × 10(-10) mutations per site per generation, similar to that reported for the distantly related budding yeast Saccharomyces cerevisiae. However, these two yeast species differ dramatically in their spectrum of base substitutions, the types of indels (S. pombe is more prone to insertions), and the pattern of selection required to counteract a strong AT-biased mutation rate. Overall, our results indicate that GC-biased gene conversion does not play a major role in shaping the nucleotide composition of the S. pombe genome and suggest that the mechanisms of DNA maintenance may have diverged significantly between fission and budding yeasts. Unexpectedly, CpG sites appear to be excessively liable to mutation in both species despite the likely absence of DNA methylation.

The Rate and Molecular Spectrum of Spontaneous Mutations in the GC-Rich Multichromosome Genome of Burkholderia cenocepacia.

Spontaneous mutations are ultimately essential for evolutionary change and are also the root cause of many diseases. However, until recently, both biological and technical barriers have prevented detailed analyses of mutation profiles, constraining our understanding of the mutation process to a few model organisms and leaving major gaps in our understanding of the role of genome content and structure on mutation. Here, we present a genome-wide view of the molecular mutation spectrum in Burkholderia cenocepacia, a clinically relevant pathogen with high %GC content and multiple chromosomes. We find that B. cenocepacia has low genome-wide mutation rates with insertion-deletion mutations biased toward deletions, consistent with the idea that deletion pressure reduces prokaryotic genome sizes. Unlike prior studies of other organisms, mutations in B. cenocepacia are not AT biased, which suggests that at least some genomes with high %GC content experience unusual base-substitution mutation pressure. Importantly, we also observe variation in both the rates and spectra of mutations among chromosomes and elevated G:C > T:A transversions in late-replicating regions. Thus, although some patterns of mutation appear to be highly conserved across cellular life, others vary between species and even between chromosomes of the same species, potentially influencing the evolution of nucleotide composition and genome architecture.

Background Mutational Features of the Radiation-Resistant Bacterium Deinococcus radiodurans.

Deinococcus bacteria are extremely resistant to radiation, oxidation, and desiccation. Resilience to these factors has been suggested to be due to enhanced damage prevention and repair mechanisms, as well as highly efficient antioxidant protection systems. Here, using mutation-accumulation experiments, we find that the GC-rich Deinococcus radiodurans has an overall background genomic mutation rate similar to that of E. coli, but differs in mutation spectrum, with the A/T to G/C mutation rate (based on a total count of 88 A:T → G:C transitions and 82 A:T → C:G transversions) per site per generation higher than that in the other direction (based on a total count of 157 G:C → A:T transitions and 33 G:C → T:A transversions). We propose that this unique spectrum is shaped mainly by the abundant uracil DNA glycosylases reducing G:C → A:T transitions, adenine methylation elevating A:T → C:G transversions, and absence of cytosine methylation decreasing G:C → A:T transitions. As opposed to the greater than 100× elevation of the mutation rate in MMR(-) (DNA Mismatch Repair deficient) strains of most other organisms, MMR(-) D. radiodurans only exhibits a 4-fold elevation, raising the possibility that other DNA repair mechanisms compensate for a relatively low-efficiency DNA MMR pathway. As D. radiodurans has plentiful insertion sequence (IS) elements in the genome and the activities of IS elements are rarely directly explored, we also estimated the insertion (transposition) rate of the IS elements to be 2.50 × 10(-3) per genome per generation in the wild-type strain; knocking out MMR did not elevate the IS element insertion rate in this organism.

Asymmetric Context-Dependent Mutation Patterns Revealed through Mutation-Accumulation Experiments.

Despite the general assumption that site-specific mutation rates are independent of the local sequence context, a growing body of evidence suggests otherwise. To further examine context-dependent patterns of mutation, we amassed 5,645 spontaneous mutations in wild- type (WT) and mismatch-repair deficient (MMR(-)) mutation-accumulation (MA) lines of the gram-positive model organism Bacillus subtilis. We then analyzed>7,500 spontaneous base-substitution mutations across B. subtilis, Escherichia coli, and Mesoplasma florum WT and MMR(-) MA lines, finding a context-dependent mutation pattern that is asymmetric around the origin of replication. Different neighboring nucleotides can alter site-specific mutation rates by as much as 75-fold, with sites neighboring G:C base pairs or dimers involving alternating pyrimidine-purine and purine-pyrimidine nucleotides having significantly elevated mutation rates. The influence of context-dependent mutation on genome architecture is strongest in M. florum, consistent with the reduced efficiency of selection in organisms with low effective population size. If not properly accounted for, the disparities arising from patterns of context-dependent mutation can significantly influence interpretations of positive and purifying selection.