Mutation rate, spectrum, topology, and context-dependency in the DNA mismatch repair-deficient Pseudomonas fluorescens ATCC948

Mutational features of chromids and chromosomes in Pseudoalteromonas provide new insights into the evolution of secondary replicons

The genomes of multi-replicon bacteria are composed of a primary replicon (the chromosome) and secondary replicons (chromids). Currently, there is a lack of understanding of the mutation features and evolutionary patterns of these different replicons. Specifically, in the genus Pseudoalteromonas, the chromids of multi-replicon species exhibit both unidirectional and bidirectional replication. Here, we investigated the similarities and differences between chromosomes and chromids in sequence composition and gene synteny of Pseudoalteromonas species by comparative genomic analysis, as well as the spontaneous mutation features of different replicons by mutation accumulation (MA) experiments combined with whole-genome sequencing strategy (MA-WGS). MA-WGS analysis revealed that there was no significant difference between chromids and chromosomes in the mutation rate or mutation spectrum of P. sp. LC0214 (where the chromid is unidirectional in replication) and P. sp. JCM12884T (where the chromid is bidirectional in replication). In addition, the context-dependence and variation pattern of the base-pair substitutions (BPSs) rates of the entire replicons exhibited differences that may be caused by the different replication directions of the chromids. The results of this study provide a new theoretical foundation for an in-depth understanding of the origin and evolution of chromids in multi-replicon bacterial species and facilitate further exploration of the complex mechanisms of bacterial diversity.IMPORTANCEDe novo mutations are a critical driving force in species evolution. Currently, there is a lack of sufficient research on the influence of replicon types on the occurrence of genomic mutations in bacteria. Moreover, the scarcity in systematic analysis and comparison of spontaneous mutation features between different replicons results in the limited information on the evolutionary dynamics of multi-replicon species. The diversity of replication direction in the multi-replicon species of the genus Pseudoalteromonas provides a unique opportunity for studying the impact of replication direction on the patterns of mutation. In addition to the composition characteristics between chromosomes and chromids, the spontaneous mutation rates in the context-dependence and variation pattern of the base-pair substitutions (BPSs) across different replicons within Pseudoalteromonas species revealed in this study provide valuable insights into the evolutionary dynamics of bacterial secondary replicons.

Elevated rates and biased spectra of mutations in anaerobically cultured lactic acid bacteria

The rate, spectrum, and biases of mutations represent a fundamental force shaping biological evolution. Convention often attributes oxidative DNA damage as a major driver of spontaneous mutations. Yet, despite the contribution of oxygen to mutagenesis and the ecological, industrial, and biomedical importance of anaerobic organisms, relatively little is known about the mutation rates and spectra of anaerobic species. Here, we present the rates and spectra of spontaneous mutations assessed anaerobically over 1000 generations for three fermentative lactic acid bacteria species with varying levels of aerotolerance: Lactobacillus acidophilus, Lactobacillus crispatus, and Lactococcus lactis. Our findings reveal highly elevated mutation rates compared to the average rates observed in aerobically respiring bacteria with mutations strongly biased towards transitions, emphasizing the prevalence of spontaneous deamination in these anaerobic species and highlighting the inherent fragility of purines even under conditions that minimize oxidative stress. Beyond these overarching patterns, we identify several novel mutation dynamics: positional mutation bias around the origin of replication in Lb. acidophilus, a significant disparity between observed and equilibrium GC content in Lc. lactis, and repeated independent deletions of spacer sequences from within the CRISPR locus in Lb. crispatus providing mechanistic insights into the evolution of bacterial adaptive immunity. Overall, our study provides new insights into the mutational landscape of anaerobes, revealing how non-oxygenic factors shape mutation rates and influence genome evolution.

Mutation bias and adaptation in bacteria

Genetic mutation, which provides the raw material for evolutionary adaptation, is largely a stochastic force. However, there is ample evidence showing that mutations can also exhibit strong biases, with some mutation types and certain genomic positions mutating more often than others. It is becoming increasingly clear that mutational bias can play a role in determining adaptive outcomes in bacteria in both the laboratory and the clinic. As such, understanding the causes and consequences of mutation bias can help microbiologists to anticipate and predict adaptive outcomes. In this review, we provide an overview of the mechanisms and features of the bacterial genome that cause mutational biases to occur. We then describe the environmental triggers that drive these mechanisms to be more potent and outline the adaptive scenarios where mutation bias can synergize with natural selection to define evolutionary outcomes. We conclude by describing how understanding mutagenic genomic features can help microbiologists predict areas sensitive to mutational bias, and finish by outlining future work that will help us achieve more accurate evolutionary forecasts.

The divergence of mutation rates and spectra across the Tree of Life

Owing to advances in genome sequencing, genome stability has become one of the most scrutinized cellular traits across the Tree of Life. Despite its centrality to all things biological, the mutation rate (per nucleotide site per generation) ranges over three orders of magnitude among species and several-fold within individual phylogenetic lineages. Within all major organismal groups, mutation rates scale negatively with the effective population size of a species and with the amount of functional DNA in the genome. This relationship is most parsimoniously explained by the drift-barrier hypothesis, which postulates that natural selection typically operates to reduce mutation rates until further improvement is thwarted by the power of random genetic drift. Despite this constraint, the molecular mechanisms underlying DNA replication fidelity and repair are free to wander, provided the performance of the entire system is maintained at the prevailing level. The evolutionary flexibility of the mutation rate bears on the resolution of several prior conundrums in phylogenetic and population-genetic analysis and raises challenges for future applications in these areas.

Generalised interrelations among mutation rates drive the genomic compliance of Chargaff's second parity rule

Chargaff's second parity rule (PR-2), where the complementary base and k-mer contents are matching within the same strand of a double stranded DNA (dsDNA), is a phenomenon that invited many explanations. The strict compliance of nearly all nuclear dsDNA to PR-2 implies that the explanation should also be similarly adamant. In this work, we revisited the possibility of mutation rates driving PR-2 compliance. Starting from the assumption-free approach, we constructed kinetic equations for unconstrained simulations. The results were analysed for their PR-2 compliance by employing symbolic regression and machine learning techniques. We arrived to a generalised set of mutation rate interrelations in place in most species that allow for their full PR-2 compliance. Importantly, our constraints explain PR-2 in genomes out of the scope of the prior explanations based on the equilibration under mutation rates with simpler no-strand-bias constraints. We thus reinstate the role of mutation rates in PR-2 through its molecular core, now shown, under our formulation, to be tolerant to previously noted strand biases and incomplete compositional equilibration. We further investigate the time for any genome to reach PR-2, showing that it is generally earlier than the compositional equilibrium, and well within the age of life on Earth.

Mutation bias and the predictability of evolution

Predicting evolutionary outcomes is an important research goal in a diversity of contexts. The focus of evolutionary forecasting is usually on adaptive processes, and efforts to improve prediction typically focus on selection. However, adaptive processes often rely on new mutations, which can be strongly influenced by predictable biases in mutation. Here, we provide an overview of existing theory and evidence for such mutation-biased adaptation and consider the implications of these results for the problem of prediction, in regard to topics such as the evolution of infectious diseases, resistance to biochemical agents, as well as cancer and other kinds of somatic evolution. We argue that empirical knowledge of mutational biases is likely to improve in the near future, and that this knowledge is readily applicable to the challenges of short-term prediction. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

Silver nanoparticles elevate mutagenesis of eukaryotic genomes

Metal nanoparticles, especially silver, have been used in various medical scenarios, due to their excellent antimicrobial effects. Recent studies have shown that AgNPs do not exert mutagenic effects on target bacteria, but the degree to which they compromise eukaryotic genomes remains unclear. To study this, we evaluated the mutagenic effects of AgNPs on the fission yeast Schizosaccharomyces pombe ATCC-16979, of which ∼23% genes are homologous to human ones, at single-nucleotide resolution, and whole-genome scale by running 283 mutation accumulation lines for ∼260,000 cell divisions in total. We also explored the action and mutagenesis mechanisms using differential gene-expression analysis based on RNAseq. Upon AgNPs treatment, the genomic base-substitution mutation rate of S. pombe at four-fold degenerate sites increased by 3.46×, and small indels were prone to occur in genomic regions that are not simple sequence repeats. The G:C → T:A transversion rate was also significantly increased, likely mostly from oxidative damage. Thus, in addition to their antimicrobial potency, AgNPs might pose slight genotoxicity threats to eukaryotic and possibly human genomes, though at a low magnitude.

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.

Hypermutator strains of Pseudomonas aeruginosa reveal novel pathways of resistance to combinations of cephalosporin antibiotics and beta-lactamase inhibitors

Hypermutation due to DNA mismatch repair (MMR) deficiencies can accelerate the development of antibiotic resistance in Pseudomonas aeruginosa. Whether hypermutators generate resistance through predominantly similar molecular mechanisms to wild-type (WT) strains is not fully understood. Here, we show that MMR-deficient P. aeruginosa can evolve resistance to important broad-spectrum cephalosporin/beta-lactamase inhibitor combination antibiotics through novel mechanisms not commonly observed in WT lineages. Using whole-genome sequencing (WGS) and transcriptional profiling of isolates that underwent in vitro adaptation to ceftazidime/avibactam (CZA), we characterized the detailed sequence of mutational and transcriptional changes underlying the development of resistance. Surprisingly, MMR-deficient lineages rapidly developed high-level resistance (>256 μg/mL) largely without corresponding fixed mutations or transcriptional changes in well-established resistance genes. Further investigation revealed that these isolates had paradoxically generated an early inactivating mutation in the mexB gene of the MexAB-OprM efflux pump, a primary mediator of CZA resistance in P. aeruginosa, potentially driving an evolutionary search for alternative resistance mechanisms. In addition to alterations in a number of genes not known to be associated with resistance, 2 mutations were observed in the operon encoding the RND efflux pump MexVW. These mutations resulted in a 4- to 6-fold increase in resistance to ceftazidime, CZA, cefepime, and ceftolozane-tazobactam when engineered into a WT strain, demonstrating a potentially important and previously unappreciated mechanism of resistance to these antibiotics in P. aeruginosa. Our results suggest that MMR-deficient isolates may rapidly evolve novel resistance mechanisms, sometimes with complex dynamics that reflect gene inactivation that occurs with hypermutation. The apparent ease with which hypermutators may switch to alternative resistance mechanisms for which antibiotics have not been developed may carry important clinical implications.

The Impact of RNA-DNA Hybrids on Genome Integrity in Bacteria

During the essential processes of DNA replication and transcription, RNA-DNA hybrid intermediates are formed that pose significant risks to genome integrity when left unresolved. To manage RNA-DNA hybrids, all cells rely on RNase H family enzymes that specifically cleave the RNA portion of the many different types of hybrids that form in vivo. Recent experimental advances have provided new insight into how RNA-DNA hybrids form and the consequences to genome integrity that ensue when persistent hybrids remain unresolved. Here we review the types of RNA-DNA hybrids, including R-loops, RNA primers, and ribonucleotide misincorporations, that form during DNA replication and transcription and discuss how each type of hybrid can contribute to genome instability in bacteria. Further, we discuss how bacterial RNase HI, HII, and HIII and bacterial FEN enzymes contribute to genome maintenance through the resolution of hybrids.

A near-deterministic mutational hotspot in Pseudomonas fluorescens is constructed by multiple interacting genomic features

Mutation - whilst stochastic - is frequently biased toward certain loci. When combined with selection this results in highly repeatable and predictable evolutionary outcomes. Immotile variants of the bacterium Pseudomonas fluorescens (SBW25) possess a 'mutational hotspot' that facilitates repeated occurrences of an identical de novo single nucleotide polymorphism when re-evolving motility, where ≥95% independent lines fix the mutation ntrB A289C. Identifying hotspots of similar potency in other genes and genomic backgrounds would prove valuable for predictive evolutionary models, but to do so we must understand the genomic features that enable such a hotspot to form. Here we reveal that genomic location, local nucleotide sequence, gene strandedness and presence of mismatch repair proteins operate in combination to facilitate the formation of this mutational hotspot. Our study therefore provides a framework for utilising genomic features to predict and identify hotspot positions capable of enforcing near-deterministic evolution.

Rates of Mutations and Transcript Errors in the Foodborne Pathogen Salmonella enterica subsp. enterica

Because errors at the DNA level power pathogen evolution, a systematic understanding of the rate and molecular spectra of mutations could guide the avoidance and treatment of infectious diseases. We thus accumulated tens of thousands of spontaneous mutations in 768 repeatedly bottlenecked lineages of 18 strains from various geographical sites, temporal spread, and genetic backgrounds. Entailing over ∼1.36 million generations, the resultant data yield an average mutation rate of ∼0.0005 per genome per generation, with a significant within-species variation. This is one of the lowest bacterial mutation rates reported, giving direct support for a high genome stability in this pathogen resulting from high DNA-mismatch-repair efficiency and replication-machinery fidelity. Pathogenicity genes do not exhibit an accelerated mutation rate, and thus, elevated mutation rates may not be the major determinant for the diversification of toxin and secretion systems. Intriguingly, a low error rate at the transcript level is not observed, suggesting distinct fidelity of the replication and transcription machinery. This study urges more attention on the most basic evolutionary processes of even the best-known human pathogens and deepens the understanding of their genome evolution.

Context-Dependent Substitution Dynamics in Plastid DNA Across a Wide Range of Taxonomic Groups

The influence of neighboring base composition, or context, on substitution bias at fourfold degenerate coding sites and in intergenic regions in plastid DNA is compared across the angiosperms, gymnosperms, ferns, liverworts, chlorophytes, stramenopiles and rhodophytes. An influence of flanking base G + C content on the relative rates of transitions and transversions is observed in all lineages and extends up to four nucleotides from the site of substitution in some. Despite finding context effects in all lineages, significant differences were observed between lineages. Overall, the data suggest that context is a general factor affecting mutation bias in plastid DNA but that the dynamics of the influence have evolved over time. It is also shown that, although there are similar effects of context on substitution bias at fourfold degenerate coding sites and at sites within intergenic regions, there are also small but significant differences, suggesting that there could be some selection on some of these sites and that there could be some difference in the mutation and/or repair process between coding and noncoding DNA.

Mutation rate and spectrum in obligately outcrossing Caenorhabditis elegans mutation accumulation lines subjected to RNAi-induced knockdown of the mismatch repair gene msh-2

DNA mismatch repair (MMR), an evolutionarily conserved repair pathway shared by prokaryotic and eukaryotic species alike, influences molecular evolution by detecting and correcting mismatches, thereby protecting genetic fidelity, reducing the mutational load, and preventing lethality. Herein we conduct the first genome-wide evaluation of the alterations to the mutation rate and spectrum under impaired activity of the MutSα homolog, msh-2, in Caenorhabditis elegans male–female fog-2(lf) lines. We performed mutation accumulation (MA) under RNAi-induced knockdown of msh-2 for up to 50 generations, followed by next-generation sequencing of 19 MA lines and the ancestral control. msh-2 impairment in the male–female background substantially increased the frequency of nuclear base substitutions (∼23×) and small indels (∼328×) relative to wildtype hermaphrodites. However, we observed no increase in the mutation rates of mtDNA, and copy-number changes of single-copy genes. There was a marked increase in copy-number variation of rDNA genes under MMR impairment. In C. elegans, msh-2 repairs transitions more efficiently than transversions and increases the AT mutational bias relative to wildtype. The local sequence context, including sequence complexity, G + C-content, and flanking bases influenced the mutation rate. The X chromosome exhibited lower substitution and higher indel rates than autosomes, which can either result from sex-specific mutation rates or a nonrandom distribution of mutable sites between chromosomes. Provided the observed difference in mutational pattern is mostly due to MMR impairment, our results indicate that the specificity of MMR varies between taxa, and is more efficient in detecting and repairing small indels in eukaryotes relative to prokaryotes.

DNA Replication-Transcription Conflicts Do Not Significantly Contribute to Spontaneous Mutations Due to Replication Errors in Escherichia coli

Encounters between DNA replication and transcription can cause genomic disruption, particularly when the two meet head-on. Whether these conflicts produce point mutations is debated. This paper presents detailed analyses of a large collection of mutations generated during mutation accumulation experiments with mismatch repair (MMR)-defective Escherichia coli. With MMR absent, mutations are primarily due to DNA replication errors. Overall, there were no differences in the frequencies of base pair substitutions or small indels (i.e., insertion and deletions of ≤4 bp) in the coding sequences or promoters of genes oriented codirectionally versus head-on to replication. Among a subset of highly expressed genes, there was a 2- to 3-fold bias for indels in genes oriented head-on to replication, but this difference was almost entirely due to the asymmetrical genomic locations of tRNA genes containing mononucleotide runs, which are hot spots for indels. No additional orientation bias in mutation frequencies occurred when MMR⁻ strains were also defective for transcription-coupled repair (TCR). However, in contrast to other reports, loss of TCR slightly increased the overall mutation rate, meaning that TCR is antimutagenic. There was no orientation bias in mutation frequencies among the stress response genes that are regulated by RpoS or induced by DNA damage. Thus, biases in the locations of mutational targets can account for most, if not all, apparent biases in mutation frequencies between genes oriented head-on versus codirectional to replication. In addition, the data revealed a strong correlation of the frequency of base pair substitutions with gene length but no correlation with gene expression levels.

A mutational hotspot that determines highly repeatable evolution can be built and broken by silent genetic changes

Mutational hotspots can determine evolutionary outcomes and make evolution repeatable. Hotspots are products of multiple evolutionary forces including mutation rate heterogeneity, but this variable is often hard to identify. In this work, we reveal that a near-deterministic genetic hotspot can be built and broken by a handful of silent mutations. We observe this when studying homologous immotile variants of the bacteria Pseudomonas fluorescens, AR2 and Pf0-2x. AR2 resurrects motility through highly repeatable de novo mutation of the same nucleotide in >95% lines in minimal media (ntrB A289C). Pf0-2x, however, evolves via a number of mutations meaning the two strains diverge significantly during adaptation. We determine that this evolutionary disparity is owed to just 6 synonymous variations within the ntrB locus, which we demonstrate by swapping the sites and observing that we are able to both break (>95% to 0%) and build (0% to 80%) a deterministic mutational hotspot. Our work reveals a key role for silent genetic variation in determining adaptive outcomes.

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.

Spatiotemporal Manipulation of the Mismatch Repair System of Pseudomonas putida Accelerates Phenotype Emergence

The development of complex phenotypes in industrially relevant bacteria is a major goal of metabolic engineering, which encompasses the implementation of both rational and random approaches. In the latter case, several tools have been developed toward increasing mutation frequencies, yet the precise control of mutagenesis processes in cell factories continues to represent a significant technical challenge. Pseudomonas species are endowed with one of the most efficient DNA mismatch repair (MMR) systems found in the bacterial domain. Here, we investigated if the endogenous MMR system could be manipulated as a general strategy to artificially alter mutation rates in Pseudomonas species. To bestow a conditional mutator phenotype in the platform bacterium Pseudomonas putida, we constructed inducible mutator devices to modulate the expression of the dominant-negative mutL^E36K allele. Regulatable overexpression of mutL^E36K in a broad-host-range, easy-to-cure plasmid format resulted in a transitory inhibition of the MMR machinery, leading to a significant increase (up to 438-fold) in DNA mutation frequencies and a heritable fixation of mutations in the genome. Following such an accelerated mutagenesis-followed by selection approach, three phenotypes were successfully evolved: resistance to antibiotics streptomycin and rifampicin (either individually or combined) and reversion of a synthetic uracil auxotrophy. Thus, these mutator devices could be applied to accelerate the evolution of metabolic pathways in long-term evolutionary experiments, alternating cycles of (inducible) mutagenesis coupled to selection schemes toward the desired phenotype(s).

Biodiversity-based development and evolution: the emerging research systems in model and non-model organisms

Evolutionary developmental biology, or Evo-Devo for short, has become an established field that, broadly speaking, seeks to understand how changes in development drive major transitions and innovation in organismal evolution. It does so via integrating the principles and methods of many subdisciplines of biology. Although we have gained unprecedented knowledge from the studies on model organisms in the past decades, many fundamental and crucially essential processes remain a mystery. Considering the tremendous biodiversity of our planet, the current model organisms seem insufficient for us to understand the evolutionary and physiological processes of life and its adaptation to exterior environments. The currently increasing genomic data and the recently available gene-editing tools make it possible to extend our studies to non-model organisms. In this review, we review the recent work on the regulatory signaling of developmental and regeneration processes, environmental adaptation, and evolutionary mechanisms using both the existing model animals such as zebrafish and Drosophila, and the emerging nonstandard model organisms including amphioxus, ascidian, ciliates, single-celled phytoplankton, and marine nematode. In addition, the challenging questions and new directions in these systems are outlined as well.

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.

Evolutionary Trajectory of the Replication Mode of Bacterial Replicons

Chromosome replication is an essential process for cell division. The mode of chromosome replication has important impacts on the structure of the chromosome and replication speed.

As typical bacterial replicons, circular chromosomes replicate bidirectionally and circular plasmids replicate either bidirectionally or unidirectionally. Whereas the finding of chromids (plasmid-derived chromosomes) in multiple bacterial lineages provides circumstantial evidence that chromosomes likely evolved from plasmids, all experimentally assayed chromids were shown to use bidirectional replication. Here, we employed a model system, the marine bacterial genus Pseudoalteromonas, members of which consistently carry a chromosome and a chromid. We provide experimental and bioinformatic evidence that while chromids in a few strains replicate bidirectionally, most replicate unidirectionally. This is the first experimental demonstration of the unidirectional replication mode in bacterial chromids. Phylogenomic and comparative genomic analyses showed that the bidirectional replication evolved only once from a unidirectional ancestor and that this transition was associated with insertions of exogenous DNA and relocation of the replication terminus region (ter2) from near the origin site (ori2) to a position roughly opposite it. This process enables a plasmid-derived chromosome to increase its size and expand the bacterium’s metabolic versatility while keeping its replication synchronized with that of the main chromosome. A major implication of our study is that the uni- and bidirectionally replicating chromids may represent two stages on the evolutionary trajectory from unidirectionally replicating plasmids to bidirectionally replicating chromosomes in bacteria. Further bioinformatic analyses predicted unidirectionally replicating chromids in several unrelated bacterial phyla, suggesting that evolution from unidirectionally to bidirectionally replicating replicons occurred multiple times in bacteria.

Pseudouridines of tRNA Anticodon Stem-Loop Have Unexpected Role in Mutagenesis in Pseudomonas sp

Pseudouridines are known to be important for optimal translation. In this study we demonstrate an unexpected link between pseudouridylation of tRNA and mutation frequency in Pseudomonas species. We observed that the lack of pseudouridylation activity of pseudouridine synthases TruA or RluA elevates the mutation frequency in Pseudomonas putida 3 to 5-fold. The absence of TruA but not RluA elevates mutation frequency also in Pseudomonas aeruginosa. Based on the results of genetic studies and analysis of proteome data, the mutagenic effect of the pseudouridylation deficiency cannot be ascribed to the involvement of error-prone DNA polymerases or malfunctioning of DNA repair pathways. In addition, although the deficiency in TruA-dependent pseudouridylation made P. putida cells more sensitive to antimicrobial compounds that may cause oxidative stress and DNA damage, cultivation of bacteria in the presence of reactive oxygen species (ROS)-scavenging compounds did not eliminate the mutator phenotype. Thus, the elevated mutation frequency in the absence of tRNA pseudouridylation could be the result of a more specific response or, alternatively, of a cumulative effect of several small effects disturbing distinct cellular functions, which remain undetected when studied independently. This work suggests that pseudouridines link the translation machinery to mutation frequency.

Genomic Landscape of Mutational Biases in the Pacific Oyster Crassostrea gigas

Mutation is a driving force of evolution that has been shaped by natural selection and is universally biased. Previous studies determined genome-wide mutational patterns for several species and investigated the heterogeneity of mutational patterns at fine-scale levels. However, little evidence of the heterogeneity of mutation rates over large genomic regions was shown. Hence, the mutational patterns of different large-scale genomic regions and their association with selective pressures still need to be explored. As the second most species-rich animal phylum, little is known about the mutational patterns in Mollusca, especially oysters. In this study, the mutational bias patterns are characterized by using whole-genome resequencing data in the Crassostrea gigas genome. I studied the genome-wide relative rates of the pair mutations and found that the predominant mutation is GC -> AT, irrespective of the genomic regions. This analysis reveals that mutational biases were associated with gene expression levels across the C. gigas genome. Genes with higher expression levels and breadth expression patterns, longer coding length, and more exon numbers had relatively higher GC -> AT rates. I also found that genes with larger dN/dS values had relatively higher GC -> AT rates. This work represents the first comprehensive research on the mutational biases in Mollusca species. Here, I comprehensively investigated the relationships between mutational biases with some intrinsic genetic factors and evolutionary indicators and proposed that selective pressures are important forces shaping the mutational biases across the C. gigas genome.

Estimation of the Genome-Wide Mutation Rate and Spectrum in the Archaeal Species Haloferax volcanii

Organisms adapted to life in extreme habitats (extremophiles) can further our understanding of the mechanisms of genetic stability, particularly replication and repair. Despite the harsh environmental conditions they endure, these extremophiles represent a great deal of the Earth's biodiversity. Here, for the first time in a member of the archaeal domain, we report a genome-wide assay of spontaneous mutations in the halophilic species Haloferax volcanii using a direct and unbiased method: mutation accumulation experiments combined with deep whole-genome sequencing. H. volcanii is a key model organism not only for the study of halophilicity, but also for archaeal biology in general. Our methods measure the genome-wide rate, spectrum, and spatial distribution of spontaneous mutations. The estimated base substitution rate of 3.15 × 10^-10 per site per generation, or 0.0012 per genome per generation, is similar to the value found in mesophilic prokaryotes (optimal growth at ∼20-45°). This study contributes to a comprehensive phylogenetic view of how evolutionary forces and molecular mechanisms shape the rate and molecular spectrum of mutations across the tree of life.

Universally high transcript error rates in bacteria

Errors can occur at any level during the replication and transcription of genetic information. Genetic mutations derived mainly from replication errors have been extensively studied. However, fundamental details of transcript errors, such as their rate, molecular spectrum, and functional effects, remain largely unknown. To globally identify transcript errors, we applied an adapted rolling-circle sequencing approach to Escherichia coli, Bacillus subtilis, Agrobacterium tumefaciens, and Mesoplasma florum, revealing transcript-error rates 3 to 4 orders of magnitude higher than the corresponding genetic mutation rates. The majority of detected errors would result in amino-acid changes, if translated. With errors identified from 9929 loci, the molecular spectrum and distribution of errors were uncovered in great detail. A G→A substitution bias was observed in M. florum, which apparently has an error-prone RNA polymerase. Surprisingly, an increased frequency of nonsense errors towards the 3' end of mRNAs was observed, suggesting a Nonsense-Mediated Decay-like quality-control mechanism in prokaryotes.

Parallel Evolution of Tobramycin Resistance across Species and Environments

The rise of antimicrobial resistance is a leading medical threat, motivating efforts to forecast both its evolutionary dynamics and its genetic causes. Aminoglycosides are a major class of antibiotics that disrupt translation, but resistance may occur by a number of mechanisms. Here, we show the repeated evolution of resistance to the aminoglycoside tobramycin in both P. aeruginosa and A. baumannii via mutations in fusA1, encoding elongation factor G, and ptsP, encoding the nitrogen-specific phosphotransferase system. Laboratory evolution and whole-population genome sequencing were used to identify these targets, but mutations at identical amino acid positions were also found in published genomes of diverse bacterial species and clinical isolates. We also identified other resistance mechanisms associated with growth in biofilms that likely interfere with drug binding or uptake. Characterizing the evolution of multiple species in the presence of antibiotics can identify new, repeatable causes of resistance that may be predicted and counteracted by alternative treatment.

Different species exposed to a common stress may adapt by mutations in shared pathways or in unique systems, depending on how past environments have molded their genomes. Understanding how diverse bacterial pathogens evolve in response to an antimicrobial treatment is a pressing example of this problem, where discovery of molecular parallelism could lead to clinically useful predictions. Evolution experiments with pathogens in environments containing antibiotics, combined with periodic whole-population genome sequencing, can be used to identify many contending routes to antimicrobial resistance. We separately propagated two clinically relevant Gram-negative pathogens, Pseudomonas aeruginosa and Acinetobacter baumannii, in increasing concentrations of tobramycin in two different environments each: planktonic and biofilm. Independently of the pathogen, the populations adapted to tobramycin selection by parallel evolution of mutations in fusA1, encoding elongation factor G, and ptsP, encoding phosphoenolpyruvate phosphotransferase. As neither gene is a direct target of this aminoglycoside, mutations to either are unexpected and underreported causes of resistance. Additionally, both species acquired antibiotic resistance-associated mutations that were more prevalent in the biofilm lifestyle than in the planktonic lifestyle; these mutations were in electron transport chain components in A. baumannii and lipopolysaccharide biosynthesis enzymes in P. aeruginosa populations. Using existing databases, we discovered site-specific parallelism of fusA1 mutations that extends across bacterial phyla and clinical isolates. This study suggests that strong selective pressures, such as antibiotic treatment, may result in high levels of predictability in molecular targets of evolution, despite differences between organisms’ genetic backgrounds and environments.

Ancestral niche separation and evolutionary rate differentiation between sister marine flavobacteria lineages

Marine flavobacteria are specialists for polysaccharide degradation. They dominate in habitats enriched with polysaccharides, but are also prevalent in pelagic environments where polysaccharides are less available. These niches are likely occupied by distinct lineages, but evolutionary processes underlying their niche differentiation remain elusive. Here, genomic analyses and physiological assays indicate that the sister flavobacteria lineages Leeuwenhoekiella and Nonlabens likely explore polysaccharide-rich macroalgae and polysaccharide-poor pelagic niches respectively. Phylogenomic analyses inferred that the niche separation likely occurred anciently and coincided with increased sequence evolutionary rate in Nonlabens compared with Leeuwenhoekiella. Further analyses ruled out the known mechanisms likely driving evolutionary rate acceleration, including reduced selection efficiency, decreased generation time and increased mutation rate. In particular, the mutation rates were determined using an unbiased experimental method, which measures the present-day populations and may not reflect ancestral populations. These data collectively lead to a new hypothesis that an ancestral and transient mutation rate increase resulted in evolutionary rate increase in Nonlabens. This hypothesis was supported by inferring that gains and losses of genes involved in SOS response, a mechanism known to drive transiently increased mutation rate, coincided with evolutionary rate acceleration. Our analyses highlight the evolutionary mechanisms underlying niche differentiation of flavobacteria lineages.

Mutation and Recombination Rates Vary Across Bacterial Chromosome

Bacteria evolve as a result of mutations and acquisition of foreign DNA by recombination processes. A growing body of evidence suggests that mutation and recombination rates are not constant across the bacterial chromosome. Bacterial chromosomal DNA is organized into a compact nucleoid structure which is established by binding of the nucleoid-associated proteins (NAPs) and other proteins. This review gives an overview of recent findings indicating that the mutagenic and recombination processes in bacteria vary at different chromosomal positions. Involvement of NAPs and other possible mechanisms in these regional differences are discussed. Variations in mutation and recombination rates across the bacterial chromosome may have implications in the evolution of bacteria.

The Symmetrical Wave Pattern of Base-Pair Substitution Rates across the Escherichia coli Chromosome Has Multiple Causes

Mutation accumulation experiments followed by whole-genome sequencing have revealed that, for several bacterial species, the rate of base-pair substitutions (BPSs) is not constant across the chromosome but varies in a wave-like pattern that is symmetrical about the origin of replication. The experiments reported here demonstrated that, in Escherichia coli, several interacting factors determine the wave. The origin is a major driver of BPS rates. When it is relocated, the BPS rates in a 1,000-kb region surrounding the new origin reproduce the pattern that surrounds the normal origin. However, the pattern across distant regions of the chromosome is unaltered and thus must be determined by other factors. Increasing the deoxynucleoside triphosphate (dNTP) concentration shifts the wave pattern away from the origin, supporting the hypothesis that fluctuations in dNTP pools coincident with replication firing contribute to the variations in the mutation rate. The nucleoid binding proteins (HU and Fis) and the terminus organizing protein (MatP) are also major factors. These proteins alter the three-dimensional structure of the DNA, and results suggest that mutation rates increase when highly structured DNA is replicated. Biases in error correction by proofreading and mismatch repair, both of which may be responsive to dNTP concentrations and DNA structure, also are major determinants of the wave pattern. These factors should apply to most bacterial and, possibly, eukaryotic genomes and suggest that different areas of the genome evolve at different rates.IMPORTANCE It has been found in several species of bacteria that the rate at which single base pairs are mutated is not constant across the genome but varies in a wave-like pattern that is symmetrical about the origin of replication. Using Escherichia coli as our model system, we show that this pattern is the result of several interconnected factors. First, the timing and progression of replication are important in determining the wave pattern. Second, the three-dimensional structure of the DNA is also a factor, and the results suggest that mutation rates increase when highly structured DNA is replicated. Finally, biases in error correction, which may be responsive both to the progression of DNA synthesis and to DNA structure, are major determinants of the wave pattern. These factors should apply to most bacterial and, possibly, eukaryotic genomes and suggest that different areas of the genome evolve at different rates.

Specificity of the DNA Mismatch Repair System (MMR) and Mutagenesis Bias in Bacteria

The mutation rate of an organism is influenced by the interaction of evolutionary forces such as natural selection and genetic drift. However, the mutation spectrum (i.e., the frequency distribution of different types of mutations) can be heavily influenced by DNA repair. Using mutation-accumulation lines of the extremophile bacterium Deinococcus radiodurans ΔmutS1 and the model soil bacterium Pseudomonas fluorescens wild-type and MMR- (Methyl-dependent Mismatch Repair-deficient) strains, we report the mutational features of these two important bacteria. We find that P. fluorescens has one of the highest MMR repair efficiencies among tested bacteria. We also discover that MMR of D. radiodurans preferentially repairs deletions, contrary to all other bacteria examined. We then, for the first time, quantify genome-wide efficiency and specificity of MMR in repairing different genomic regions and mutation types, by evaluating the P. fluorescens and D. radiodurans mutation data sets, along with previously reported ones of Bacillus subtilis subsp. subtilis, Escherichia coli, Vibrio cholerae, and V. fischeri. MMR in all six bacteria shares two general features: 1) repair efficiency is influenced by the neighboring base composition for both transitions and transversions, not limited to transversions as previously reported; and 2) MMR only recognizes indels <4 bp in length. This study demonstrates the power of mutation accumulation lines in quantifying DNA repair and mutagenesis patterns.

Influence of Electron-Holes on DNA Sequence-Specific Mutation Rates

Biases in mutation rate can influence molecular evolution, yielding rates of evolution that vary widely in different parts of the genome and even among neighboring nucleotides. Here, we explore one possible mechanism of influence on sequence-specific mutation rates, the electron–hole, which can localize and potentially trigger a replication mismatch. A hole is a mobile site of positive charge created during one-electron oxidation by, for example, radiation, contact with a mutagenic agent, or oxidative stress. Its quantum wavelike properties cause it to localize at various sites with probabilities that vary widely, by orders of magnitude, and depend strongly on the local sequence. We find significant correlations between hole probabilities and mutation rates within base triplets, observed in published mutation accumulation experiments on four species of bacteria. We have also computed hole probability spectra for hypervariable segment I of the human mtDNA control region, which contains several mutational hotspots, and for heptanucleotides in noncoding regions of the human genome, whose polymorphism levels have recently been reported. We observe significant correlations between hole probabilities, and context-specific mutation and substitution rates. The correlation with hole probability cannot be explained entirely by CpG methylation in the heptanucleotide data. Peaks in hole probability tend to coincide with mutational hotspots, even in mtDNA where CpG methylation is rare. Our results suggest that hole-enhanced mutational mechanisms, such as oxidation-stabilized tautomerization and base deamination, contribute to molecular evolution.

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.

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.

Bacterial EndoMS/NucS acts as a clamp-mediated mismatch endonuclease to prevent asymmetric accumulation of replication errors

Mismatch repair (MMR) systems based on MutS eliminate mismatches originating from replication errors. Despite extensive conservation of mutS homologues throughout the three domains of life, Actinobacteria and some archaea do not have genes homologous to mutS. Here, we report that EndoMS/NucS of Corynebacterium glutamicum is the mismatch-specific endonuclease that functions cooperatively with a sliding clamp. EndoMS/NucS function in MMR was fully dependent on physical interaction between EndoMS/NucS and sliding clamp. A combination of endoMS/nucS gene disruption and a mutation in dnaE, which reduced the fidelity of DNA polymerase, increased the mutation rate synergistically and confirmed the participation of EndoMS in replication error correction. EndoMS specifically cleaved G/T, G/G and T/T mismatches in vitro, and such substrate specificity was consistent with the mutation spectrum observed in genome-wide analyses. The observed substrate specificity of EndoMS, together with the effects of endoMS gene disruption, led us to speculate that the MMR system, regardless of the types of proteins in the system, evolved to address asymmetrically occurring replication errors in which G/T mismatches occur much more frequently than C/A mismatches.

Real-time dynamics of mutagenesis reveal the chronology of DNA repair and damage tolerance responses in single cells

Evolutionary processes are driven by diverse molecular mechanisms that act in the creation and prevention of mutations. It remains unclear how these mechanisms are regulated because limitations of existing mutation assays have precluded measuring how mutation rates vary over time in single cells. Toward this goal, I detected nascent DNA mismatches as a proxy for mutagenesis and simultaneously followed gene expression dynamics in single Escherichia coli cells using microfluidics. This general microscopy-based approach revealed the real-time dynamics of mutagenesis in response to DNA alkylation damage and antibiotic treatments. It also enabled relating the creation of DNA mismatches to the chronology of the underlying molecular processes. By avoiding population averaging, I discovered cell-to-cell variation in mutagenesis that correlated with heterogeneity in the expression of alternative responses to DNA damage. Pulses of mutagenesis are shown to arise from transient DNA repair deficiency. Constitutive expression of DNA repair pathways and induction of damage tolerance by the SOS response compensate for delays in the activation of inducible DNA repair mechanisms, together providing robustness against the toxic and mutagenic effects of DNA alkylation damage.

Similar Mutation Rates but Highly Diverse Mutation Spectra in Ascomycete and Basidiomycete Yeasts

Yeast species are extremely diverse and not monophyletic. Because the majority of yeast research focuses on ascomycetes, the mutational determinants of genetic diversity across yeast species are not well understood. By combining mutation-accumulation techniques with whole-genome sequencing, we resolved the genomic mutation rate and spectrum of the oleaginous (oil-producing) ‘red yeast’ Rhodotorula toruloides, the first such study in the fungal phylum Basidiomycota. We find that the mutation spectrum is quite different from what has been observed in all other studied unicellular eukaryotes, but similar to that in most bacteria—a predominance of transitions relative to transversions. Rhodotorula toruloides has a significantly higher A:T→G:C transition rate—possibly elevated by the abundant flanking G/C nucleotides in the GC-rich genome, as well as a much lower G:C→T:A transversion rate. In spite of these striking differences, there are substantial consistencies between R. toruloides and the ascomycete model yeasts: a spontaneous base-substitution mutation rate of 1.90 × 10 ⁻¹⁰ per site per cell division as well as an elevated mutation rate at non-methylated 5'CpG3' sites. These results imply the evolution of variable mutation spectra in the face of similar mutation rates in yeasts.

Sources of spontaneous mutagenesis in bacteria

Mutations in an organism's genome can arise spontaneously, that is, in the absence of exogenous stress and prior to selection. Mutations are often neutral or deleterious to individual fitness but can also provide genetic diversity driving evolution. Mutagenesis in bacteria contributes to the already serious and growing problem of antibiotic resistance. However, the negative impacts of spontaneous mutagenesis on human health are not limited to bacterial antibiotic resistance. Spontaneous mutations also underlie tumorigenesis and evolution of drug resistance. To better understand the causes of genetic change and how they may be manipulated in order to curb antibiotic resistance or the development of cancer, we must acquire a mechanistic understanding of the major sources of mutagenesis. Bacterial systems are particularly well-suited to studying mutagenesis because of their fast growth rate and the panoply of available experimental tools, but efforts to understand mutagenic mechanisms can be complicated by the experimental system employed. Here, we review our current understanding of mutagenic mechanisms in bacteria and describe the methods used to study mutagenesis in bacterial systems.

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.

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.

Anaerobically Grown Escherichia coli Has an Enhanced Mutation Rate and Distinct Mutational Spectra

Oxidative stress is a major cause of mutation but little is known about how growth in the absence of oxygen impacts the rate and spectrum of mutations. We employed long-term mutation accumulation experiments to directly measure the rates and spectra of spontaneous mutation events in Escherichia coli populations propagated under aerobic and anaerobic conditions. To detect mutations, whole genome sequencing was coupled with methods of analysis sufficient to identify a broad range of mutational classes, including structural variants (SVs) generated by movement of repetitive elements. The anaerobically grown populations displayed a mutation rate nearly twice that of the aerobic populations, showed distinct asymmetric mutational strand biases, and greater insertion element activity. Consistent with mutation rate and spectra observations, genes for transposition and recombination repair associated with SVs were up-regulated during anaerobic growth. Together, these results define differences in mutational spectra affecting the evolution of facultative anaerobes.

Population Genomics of Daphnia pulex

Using data from 83 isolates from a single population, the population genomics of the microcrustacean Daphnia pulex are described and compared to current knowledge for the only other well-studied invertebrate, Drosophila melanogaster These two species are quite similar with respect to effective population sizes and mutation rates, although some features of recombination appear to be different, with linkage disequilibrium being elevated at short ([Formula: see text] bp) distances in D. melanogaster and at long distances in D. pulex The study population adheres closely to the expectations under Hardy-Weinberg equilibrium, and reflects a past population history of no more than a twofold range of variation in effective population size. Fourfold redundant silent sites and a restricted region of intronic sites appear to evolve in a nearly neutral fashion, providing a powerful tool for population genetic analyses. Amino acid replacement sites are predominantly under strong purifying selection, as are a large fraction of sites in UTRs and intergenic regions, but the majority of SNPs at such sites that rise to frequencies [Formula: see text] appear to evolve in a nearly neutral fashion. All forms of genomic sites (including replacement sites within codons, and intergenic and UTR regions) appear to be experiencing an [Formula: see text] higher level of selection scaled to the power of drift in D. melanogaster, but this may in part be a consequence of recent demographic changes. These results establish D. pulex as an excellent system for future work on the evolutionary genomics of natural populations.

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.

Loss of mutL homolog-1 (MLH1) expression promotes acquisition of oncogenic and inhibitor-resistant point mutations in tyrosine kinases

Genomic instability drives cancer progression by promoting genetic abnormalities that allow for the multi-step clonal selection of cells with growth advantages. We previously reported that the IL-9-dependent TS1 cell line sequentially acquired activating substitutions in JAK1 and JAK3 upon successive selections for growth factor independent and JAK inhibitor-resistant cells, suggestive of a defect in mutation avoidance mechanisms. In the first part of this paper, we discovered that the gene encoding mutL homolog-1 (MLH1), a key component of the DNA mismatch repair system, is silenced by promoter methylation in TS1 cells. By means of stable ectopic expression and RNA interference methods, we showed that the high frequencies of growth factor-independent and inhibitor-resistant cells with activating JAK mutations can be attributed to the absence of MLH1 expression. In the second part of this paper, we confirm the clinical relevance of our findings by showing that chronic myeloid leukemia relapses upon ABL-targeted therapy correlated with a lower expression of MLH1 messenger RNA. Interestingly, the mutational profile observed in our TS1 model, characterized by a strong predominance of T:A>C:G transitions, was identical to the one described in the literature for primitive cells derived from chronic myeloid leukemia patients. Taken together, our observations demonstrate for the first time a causal relationship between MLH1-deficiency and incidence of oncogenic point mutations in tyrosine kinases driving cell transformation and acquired resistance to kinase-targeted cancer therapies.

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.

The Fitness Effects of Spontaneous Mutations Nearly Unseen by Selection in a Bacterium with Multiple Chromosomes

Mutation accumulation (MA) experiments employ the strategy of minimizing the population size of evolving lineages to greatly reduce effects of selection on newly arising mutations. Thus, most mutations fix within MA lines independently of their fitness effects. This approach, more recently combined with genome sequencing, has detailed the rates, spectra, and biases of different mutational processes. However, a quantitative understanding of the fitness effects of mutations virtually unseen by selection has remained an untapped opportunity. Here, we analyzed the fitness of 43 sequenced MA lines of the multi-chromosome bacterium Burkholderia cenocepacia that had each undergone 5554 generations of MA and accumulated an average of 6.73 spontaneous mutations. Most lineages exhibited either neutral or deleterious fitness in three different environments in comparison with their common ancestor. The only mutational class that was significantly overrepresented in lineages with reduced fitness was the loss of the plasmid, though nonsense mutations, missense mutations, and coding insertion-deletions were also overrepresented in MA lineages whose fitness had significantly declined. Although the overall distribution of fitness effects was similar between the three environments, the magnitude and even the sign of the fitness of a number of lineages changed with the environment, demonstrating that the fitness of some genotypes was environmentally dependent. These results present an unprecedented picture of the fitness effects of spontaneous mutations in a bacterium with multiple chromosomes and provide greater quantitative support for the theory that the vast majority of spontaneous mutations are neutral or deleterious.

The repeatability of genome-wide mutation rate and spectrum estimates

Over the last decade, mutation studies have grown in popularity due to the affordability and accessibility of whole genome sequencing. As the number of species in which spontaneous mutation has been directly estimated approaches 20 across two domains of life, questions arise over the repeatability of results in such experiments. Five species were identified in which duplicate mutation studies have been performed. Across these studies the difference in estimated spontaneous mutation rate is at most, weakly significant (p < 0.01). However, a highly significant (p < 10(-5)), threefold difference in the rate of insertions/deletions (indels) exists between two recent studies in Schizosaccharomyces pombe. Upon investigation of the ancestral genome sequence for both studies, a possible anti-mutator allele was identified. The observed variation in indel rate may imply that the use of indel markers, such as microsatellites, for the investigation of genetic diversity within and among populations may be inappropriate because of the assumption of uniform mutation rate within a species.

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⁻¹⁰ 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.

Synonymous Genetic Variation in Natural Isolates of Escherichia coli Does Not Predict Where Synonymous Substitutions Occur in a Long-Term Experiment

Synonymous genetic differences vary by more than 20-fold among genes in natural isolates of Escherichia coli. One hypothesis to explain this heterogeneity is that genes with high levels of synonymous variation mutate at higher rates than genes with low synonymous variation. If so, then one would expect to observe similar mutational patterns in evolution experiments. In fact, however, the pattern of synonymous substitutions in a long-term evolution experiment with E. coli does not support this hypothesis. In particular, the extent of synonymous variation across genes in that experiment does not reflect the variation observed in natural isolates of E. coli. Instead, gene length alone predicts with high accuracy the prevalence of synonymous changes in the experimental populations. We hypothesize that patterns of synonymous variation in natural E. coli populations are instead caused by differences across genomic regions in their effective population size that, in turn, reflect different histories of recombination, horizontal gene transfer, selection, and population structure.

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.

Dissecting genetic and environmental mutation signatures with model organisms

Deep sequencing has impacted on cancer research by enabling routine sequencing of genomes and exomes to identify genetic changes associated with carcinogenesis. Researchers can now use the frequency, type, and context of all mutations in tumor genomes to extract mutation signatures that reflect the driving mutational processes. Identifying mutation signatures, however, may not immediately suggest a mechanism. Consequently, several recent studies have employed deep sequencing of model organisms exposed to discrete genetic or environmental perturbations. These studies exploit the simpler genomes and availability of powerful genetic tools in model organisms to analyze mutation signatures under controlled conditions, forging mechanistic links between mutational processes and signatures. We discuss the power of this approach and suggest that many such studies may be on the horizon.