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

Variation in the Spectrum of New Mutations among Inbred Strains of Mice

The mouse serves as a mammalian model for understanding the nature of variation from new mutations, a question that has both evolutionary and medical significance. Previous studies suggest that the rate of single-nucleotide mutations (SNMs) in mice is ∼50% of that in humans. However, information largely comes from studies involving the C57BL/6 strain, and there is little information from other mouse strains. Here, we study the mutations that accumulated in 59 mouse lines derived from four inbred strains that are commonly used in genetics and clinical research (BALB/cAnNRj, C57BL/6JRj, C3H/HeNRj, and FVB/NRj), maintained for eight to nine generations by brother-sister mating. By analyzing Illumina whole-genome sequencing data, we estimate that the average rate of new SNMs in mice is ∼μ = 6.7 × 10-9. However, there is substantial variation in the spectrum of SNMs among strains, so the burden from new mutations also varies among strains. For example, the FVB strain has a spectrum that is markedly skewed toward C→A transversions and is likely to experience a higher deleterious load than other strains, due to an increased frequency of nonsense mutations in glutamic acid codons. Finally, we observe substantial variation in the rate of new SNMs among DNA sequence contexts, CpG sites, and their adjacent nucleotides playing an important role.

Mutation Rate and Effective Population Size of the Model Cooperative Bacterium Myxococcus xanthus

Intrinsic rates of genetic mutation have diverged greatly across taxa and exhibit statistical associations with several other parameters and features. These include effective population size (Ne), genome size, and gametic multicellularity, with the latter being associated with both increased mutation rates and decreased effective population sizes. However, data sufficient to test for possible relationships between microbial multicellularity and mutation rate (µ) are lacking. Here, we report estimates of two key population-genetic parameters, Ne and µ, for Myxococcus xanthus, a bacterial model organism for the study of aggregative multicellular development, predation, and social swarming. To estimate µ, we conducted an ∼400-day mutation accumulation experiment with 46 lineages subjected to regular single colony bottlenecks prior to clonal regrowth. Upon conclusion, we sequenced one clonal-isolate genome per lineage. Given collective evolution for 85,323 generations across all lines, we calculate a per base-pair mutation rate of ∼5.5 × 10-10 per site per generation, one of the highest mutation rates among free-living eubacteria. Given our estimate of µ, we derived Ne at ∼107 from neutral diversity at four-fold degenerate sites across two dozen M. xanthus natural isolates. This estimate is below average for eubacteria and strengthens an already clear negative correlation between µ and Ne in prokaryotes. The higher and lower than average mutation rate and Ne for M. xanthus, respectively, amplify the question of whether any features of its multicellular life cycle-such as group-size reduction during fruiting-body development-or its highly structured spatial distribution have significantly influenced how these parameters have evolved.

Spontaneous mutations and mutational responses to penicillin treatment in the bacterial pathogen Streptococcus pneumoniae D39

Bacteria with functional DNA repair systems are expected to have low mutation rates due to strong natural selection for genomic stability. However, our study of the wild-type Streptococcus pneumoniae D39, a pathogen responsible for many common diseases, revealed a high spontaneous mutation rate of 0.02 per genome per cell division in mutation-accumulation (MA) lines. This rate is orders of magnitude higher than that of other non-mutator bacteria and is characterized by a high mutation bias in the A/T direction. The high mutation rate may have resulted from a reduction in the overall efficiency of selection, conferred by the tiny effective population size in nature. In line with this, S. pneumoniae D39 also exhibited the lowest DNA mismatch-repair (MMR) efficiency among bacteria. Treatment with the antibiotic penicillin did not elevate the mutation rate, as penicillin did not induce DNA damage and S. pneumoniae lacks a stress response pathway. Our findings suggested that the MA results are applicable to within-host scenarios and provide insights into pathogen evolution.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-024-00220-6.

Genome reduction occurred in early Prochlorococcus with an unusually low effective population size

In the oligotrophic sunlit ocean, the most abundant free-living planktonic bacterial lineages evolve convergently through genome reduction. The cyanobacterium Prochlorococcus responsible for 10% global oxygen production is a prominent example. The dominant theory known as "genome streamlining" posits that they have extremely large effective population sizes (Ne) such that selection for metabolic efficiency acts to drive genome reduction. Because genome reduction largely took place anciently, this theory builds on the assumption that their ancestors' Ne was similarly large. Constraining Ne for ancient ancestors is challenging because experimental measurements of extinct organisms are impossible and alternatively reconstructing ancestral Ne with phylogenetic models gives large uncertainties. Here, we develop a new strategy that leverages agent-based modeling to simulate the changes in the genome-wide ratio of radical to conservative nonsynonymous nucleotide substitution rate (dR/dC) in a possible range of Ne in ancestral populations. This proxy shows expected increases with decreases of Ne only when Ne falls to about 10 k - 100 k or lower, magnitudes characteristic of Ne of obligate endosymbiont species where drift drives genome reduction. Our simulations therefore strongly support a scenario where the primary force of Prochlorococcus genome reduction is drift rather than selection.

Challenges in estimating effective population sizes from metagenome-assembled genomes

Effective population size (Ne) plays a critical role in shaping the relative efficiency between natural selection and genetic drift, thereby serving as a cornerstone for understanding microbial ecological dynamics. Direct Ne estimation relies on neutral genetic diversity within closely related genomes, which is, however, often constrained by the culturing difficulties for the vast majority of prokaryotic lineages. Metagenome-assembled genomes (MAGs) offer a high-throughput alternative for genomic data acquisition, yet their accuracy in Ne estimation has not been fully verified. This study examines the Thermococcus genus, comprising 66 isolated strains and 29 MAGs, to evaluate the reliability of MAGs in Ne estimation. Despite the even distribution across the Thermococcus phylogeny and the comparable internal average nucleotide identity (ANI) between isolate populations and MAG populations, our results reveal consistently lower Ne estimates from MAG populations. This trend of underestimation is also observed in various MAG populations across three other bacterial genera. The underrepresentation of genetic variation in MAGs, including loss of allele frequency data and variable genomic segments, likely contributes to the underestimation of Ne. Our findings underscore the necessity for caution when employing MAGs for evolutionary studies, which often depend on high-quality genome assemblies and nucleotide-level diversity.

Molecular basis of the phenotypic variants arising from a Pseudoalteromonas lipolytica mutator

Bacterial deficiencies in the DNA repair system can produce mutator strains that promote adaptive microevolution. However, the role of mutator strains in marine Pseudoalteromonas, capable of generating various gain-of-function genetic variants within biofilms, remains largely unknown. In this study, inactivation of mutS in Pseudoalteromonas lipolytica conferred an approximately 100-fold increased resistance to various antibiotics, including ciprofloxacin, rifampicin and aminoglycoside. Furthermore, the mutator of P. lipolytica generated variants that displayed enhanced biofilm formation but reduced swimming motility, indicating a high phenotypic diversity within the ΔmutS population. Additionally, we observed a significant production rate of approximately 50 % for the translucent variants, which play important roles in biofilm formation, when the ΔmutS strain was cultured on agar plates or under shaking conditions. Using whole-genome deep-sequencing combined with genetic manipulation, we demonstrated that point mutations in AT00_17115 within the capsular biosynthesis cluster were responsible for the generation of translucent variants in the ΔmutS subpopulation, while mutations in flagellar genes fliI and flgP led to a decrease in swimming motility. Collectively, this study reveals a specific mutator-driven evolution in P. lipolytica, characterized by substantial genetic and phenotypic diversification, thereby offering a reservoir of genetic attributes associated with microbial fitness.

Evolution of a minimal cell

Possessing only essential genes, a minimal cell can reveal mechanisms and processes that are critical for the persistence and stability of life1,2. Here we report on how an engineered minimal cell3,4 contends with the forces of evolution compared with the Mycoplasma mycoides non-minimal cell from which it was synthetically derived. Mutation rates were the highest among all reported bacteria, but were not affected by genome minimization. Genome streamlining was costly, leading to a decrease in fitness of greater than 50%, but this deficit was regained during 2,000 generations of evolution. Despite selection acting on distinct genetic targets, increases in the maximum growth rate of the synthetic cells were comparable. Moreover, when performance was assessed by relative fitness, the minimal cell evolved 39% faster than the non-minimal cell. The only apparent constraint involved the evolution of cell size. The size of the non-minimal cell increased by 80%, whereas the minimal cell remained the same. This pattern reflected epistatic effects of mutations in ftsZ, which encodes a tubulin-homologue protein that regulates cell division and morphology5,6. Our findings demonstrate that natural selection can rapidly increase the fitness of one of the simplest autonomously growing organisms. Understanding how species with small genomes overcome evolutionary challenges provides critical insights into the persistence of host-associated endosymbionts, the stability of streamlined chassis for biotechnology and the targeted refinement of synthetically engineered cells2,7-9.

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.

Rapid evolution of mutation rate and spectrum in response to environmental and population-genetic challenges

Ecological and demographic factors can significantly shape the evolution of microbial populations both directly and indirectly, as when changes in the effective population size affect the efficiency of natural selection on the mutation rate. However, it remains unclear how rapidly the mutation-rate responds evolutionarily to the entanglement of ecological and population-genetic factors over time. Here, we directly assess the mutation rate and spectrum of Escherichia coli clones isolated from populations evolving in response to 1000 days of different transfer volumes and resource-replenishment intervals. The evolution of mutation rates proceeded rapidly in response to demographic and/or environmental changes, with substantial bidirectional shifts observed as early as 59 generations. These results highlight the remarkable rapidity by which mutation rates are shaped in asexual lineages in response to environmental and population-genetic forces, and are broadly consistent with the drift-barrier hypothesis for the evolution of mutation rates, while also highlighting situations in which mutator genotypes may be promoted by positive selection.

How rapidly the mutation rate responds evolutionarily to ecological and population-genetic factors over time is unclear. Here, the authors show that the evolution of mutation rates in E. coli proceeds rapidly in response to these factors with substantial bidirectional shifts.

Evolutionary jumps in bacterial GC content

Genomic GC (Guanine-Cytosine) content is a fundamental molecular trait linked with many key genomic features such as codon and amino acid use. Across bacteria, GC content is surprisingly diverse and has been studied for many decades; yet its evolution remains incompletely understood. Since it is difficult to observe GC content evolve on laboratory time scales, phylogenetic comparative approaches are instrumental; but this dimension is rarely studied systematically in the case of bacterial GC content. We applied phylogenetic comparative models to analyze GC content evolution in multiple bacterial groups across 2 major bacterial phyla. We find that GC content diversifies via a combination of gradual evolution and evolutionary "jumps." Surprisingly, unlike prior reports that solely focused on reductions in GC, we found a comparable number of jumps with both increased and decreased GC content. Overall, many of the identified jumps occur in lineages beyond the well-studied peculiar examples of endosymbiotic and AT-rich marine bacteria and do not support the predicted role of oxygen dependence. Our analysis of rapid and large shifts in GC content thus identifies new clades and novel contexts to further understand the ecological and evolutionary drivers of this important genomic trait.

Evolved resistance to a novel cationic peptide antibiotic requires high mutation supply

Background and Objectives

A key strategy for resolving the antibiotic resistance crisis is the development of new drugs with antimicrobial properties. The engineered cationic antimicrobial peptide WLBU2 (also known as PLG0206) is a promising broad-spectrum antimicrobial compound that has completed Phase I clinical studies. It has activity against Gram-negative and Gram-positive bacteria including infections associated with biofilm. No definitive mechanisms of resistance to WLBU2 have been identified.

Methodology

Here, we used experimental evolution under different levels of mutation supply and whole genome sequencing (WGS) to detect the genetic pathways and probable mechanisms of resistance to this peptide. We propagated populations of wild-type and hypermutator Pseudomonas aeruginosa in the presence of WLBU2 and performed WGS of evolved populations and clones.

Results

Populations that survived WLBU2 treatment acquired a minimum of two mutations, making the acquisition of resistance more difficult than for most antibiotics, which can be tolerated by mutation of a single target. Major targets of resistance to WLBU2 included the orfN and pmrB genes, previously described to confer resistance to other cationic peptides. More surprisingly, mutations that increase aggregation such as the wsp pathway were also selected despite the ability of WLBU2 to kill cells growing in a biofilm.

Conclusions and implications

The results show how experimental evolution and WGS can identify genetic targets and actions of new antimicrobial compounds and predict pathways to resistance of new antibiotics in clinical practice.

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.

Evidence of Genomic Diversification in a Natural Symbiotic Population Within Its Host

Planktonic cells of the luminous marine bacterium Vibrio fischeri establish themselves in the light-emitting organ of each generation of newly hatched Euprymna scolopes bobtail squid. A symbiont population is maintained within the 6 separated crypts of the organ for the ∼9-month life of the host. In the wild, the initial colonization step is typically accomplished by a handful of planktonic V. fischeri cells, leading to a species-specific, but often multi-strain, symbiont population. Within a few hours, the inoculating cells proliferate within the organ’s individual crypts, after which there is evidently no supernumerary colonization. Nevertheless, every day at dawn, the majority of the symbionts is expelled, and the regrowth of the remaining ∼5% of cells provides a daily opportunity for the population to evolve and diverge, thereby increasing its genomic diversity. To begin to understand the extent of this diversification, we characterized the light-organ population of an adult animal. First, we used 16S sequencing to determine that species in the V. fischeri clade were essentially the only ones detectable within a field-caught E. scolopes. Efforts to colonize the host with a minor species that appeared to be identified, V. litoralis, revealed that, although some cells could be imaged within the organ, they were <0.1% of the typical V. fischeri population, and did not persist. Next, we determined the genome sequences of seventy-two isolates from one side of the organ. While all these isolates were associated with one of three clusters of V. fischeri strains, there was considerable genomic diversity within this natural symbiotic population. Comparative analyses revealed a significant difference in both the number and the presence/absence of genes within each cluster; in contrast, there was little accumulation of single-nucleotide polymorphisms. These data suggest that, in nature, the light organ is colonized by a small number of V. fischeri strains that can undergo significant genetic diversification, including by horizontal-gene transfer, over the course of ∼1500 generations of growth in the organ. When the resulting population of symbionts is expelled into seawater, its genomic mix provides the genetic basis for selection during the subsequent environmental dispersal, and transmission to the next host.

Vibrio harveyi Exhibits the Growth Advantage in Stationary Phase Phenotype during Long-Term Incubation

The bioluminescent marine bacterium Vibrio harveyi can exist within a host, acting as a mutualist or a parasitic microbe, and as planktonic cells in open seawater. This study demonstrates the ability of V. harveyi populations to survive and adapt under nutrient stress conditions in the laboratory, starting in an initially rich medium. V. harveyi populations remain viable into long-term stationary phase, for at least 1 month, without the addition of nutrients. To determine whether these communities are dynamic, populations were sampled after 10, 20, and 30 days of incubation and examined for their competitive ability when cocultured with an unaged, parental population. While populations incubated for 10 or 20 days showed some fitness advantage over parental populations, only after 30 days of incubation did all populations examined outcompete parental populations in coculture, fully expressing the growth advantage in stationary phase (GASP) phenotype. The ability to express GASP, in the absence of additional nutrients after inoculation, verifies the dynamism of long-term stationary-phase V. harveyi populations, implies the ability to generate genetic diversity, and demonstrates the plasticity of the V. harveyi genome, allowing for rapid adaptation for survival in changing culture environments. Despite the dynamism, the adaptation to the changing culture environment occurs less rapidly than in Escherichia coli, possibly due to Vibrio harveyi's lower mutation frequency. IMPORTANCE Vibrio harveyi populations exist in many different niches within the ocean environment, as free-living cells, symbionts with particular squid and fish species, and parasites to other marine organisms. It is important to understand V. harveyi's ability to survive and evolve within each of these niches. This study focuses on V. harveyi's lifestyle outside the host environment, demonstrating this microbe's ability to survive long-term culturing after inoculation in an initially rich medium and revealing increased competitive fitness correlated with incubation time when aged V. harveyi populations are cocultured with unaged, parental cultures. Thus, this study highlights the development of the growth advantage in stationary phase (GASP) phenotype in V. harveyi populations suggesting a dynamic population with fluctuating genotype frequencies throughout long-term, host-independent incubation.

Prochlorococcus have low global mutation rate and small effective population size

Prochlorococcus are the most abundant free-living photosynthetic carbon-fixing organisms in the ocean. Prochlorococcus show small genome sizes, low genomic G+C content, reduced DNA repair gene pool and fast evolutionary rates, which are typical features of endosymbiotic bacteria. Nevertheless, their evolutionary mechanisms are believed to be different. Evolution of endosymbiotic bacteria is dominated by genetic drift owing to repeated population bottlenecks, whereas Prochlorococcus are postulated to have extremely large effective population sizes (N_e) and thus drift has rarely been considered. However, accurately extrapolating N_e requires measuring an unbiased global mutation rate through mutation accumulation, which is challenging for Prochlorococcus. Here, we managed this experiment over 1,065 days using Prochlorococcus marinus AS9601, sequenced genomes of 141 mutant lines and determined its mutation rate to be 3.50 × 10^-10 per site per generation. Extrapolating N_e additionally requires identifying population boundaries, which we defined using PopCOGenT and over 400 genomes related to AS9601. Accordingly, we calculated its N_e to be 1.68 × 10⁷, which is only reasonably greater than that of endosymbiotic bacteria but surprisingly smaller than that of many free-living bacteria extrapolated using the same approach. Our results therefore suggest that genetic drift is a key driver of Prochlorococcus evolution.

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.

Epigenetic modifications affect the rate of spontaneous mutations in a pathogenic fungus

Mutations are the source of genetic variation and the substrate for evolution. Genome-wide mutation rates appear to be affected by selection and are probably adaptive. Mutation rates are also known to vary along genomes, possibly in response to epigenetic modifications, but causality is only assumed. In this study we determine the direct impact of epigenetic modifications and temperature stress on mitotic mutation rates in a fungal pathogen using a mutation accumulation approach. Deletion mutants lacking epigenetic modifications confirm that histone mark H3K27me3 increases whereas H3K9me3 decreases the mutation rate. Furthermore, cytosine methylation in transposable elements (TE) increases the mutation rate 15-fold resulting in significantly less TE mobilization. Also accessory chromosomes have significantly higher mutation rates. Finally, we find that temperature stress substantially elevates the mutation rate. Taken together, we find that epigenetic modifications and environmental conditions modify the rate and the location of spontaneous mutations in the genome and alter its evolutionary trajectory.

De Novo Mutation Rate Variation and Its Determinants in Chlamydomonas

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

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.

A lasting symbiosis: how Vibrio fischeri finds a squid partner and persists within its natural host

As our understanding of the human microbiome progresses, so does the need for natural experimental animal models that promote a mechanistic understanding of beneficial microorganism-host interactions. Years of research into the exclusive symbiosis between the Hawaiian bobtail squid, Euprymna scolopes, and the bioluminescent bacterium Vibrio fischeri have permitted a detailed understanding of those bacterial genes underlying signal exchange and rhythmic activities that result in a persistent, beneficial association, as well as glimpses into the evolution of symbiotic competence. Migrating from the ambient seawater to regions deep inside the light-emitting organ of the squid, V. fischeri experiences, recognizes and adjusts to the changing environmental conditions. Here, we review key advances over the past 15 years that are deepening our understanding of these events.

Unexpectedly high mutation rate of a deep-sea hyperthermophilic anaerobic archaeon

Deep-sea hydrothermal vents resemble the early Earth, and thus the dominant Thermococcaceae inhabitants, which occupy an evolutionarily basal position of the archaeal tree and take an obligate anaerobic hyperthermophilic free-living lifestyle, are likely excellent models to study the evolution of early life. Here, we determined that unbiased mutation rate of a representative species, Thermococcus eurythermalis, exceeded that of all known free-living prokaryotes by 1-2 orders of magnitude, and thus rejected the long-standing hypothesis that low mutation rates were selectively favored in hyperthermophiles. We further sequenced multiple and diverse isolates of this species and calculated that T. eurythermalis has a lower effective population size than other free-living prokaryotes by 1-2 orders of magnitude. These data collectively indicate that the high mutation rate of this species is not selectively favored but instead driven by random genetic drift. The availability of these unusual data also helps explore mechanisms underlying microbial genome size evolution. We showed that genome size is negatively correlated with mutation rate and positively correlated with effective population size across 30 bacterial and archaeal lineages, suggesting that increased mutation rate and random genetic drift are likely two important mechanisms driving microbial genome reduction. Future determinations of the unbiased mutation rate of more representative lineages with highly reduced genomes such as Prochlorococcus and Pelagibacterales that dominate marine microbial communities are essential to test these hypotheses.

Control of Genome Stability by EndoMS/NucS-Mediated Non-Canonical Mismatch Repair

The DNA repair endonuclease EndoMS/NucS is highly conserved in Archaea and Actinobacteria. This enzyme is able to recognize and cleave dsDNA carrying a mismatched base pair, and its activity is enhanced by the interaction with the sliding clamp of the replisome. Today, EndoMS/NucS has been established as the key protein of a non-canonical mismatch repair (MMR) pathway, acting specifically in the repair of transitions and being essential for maintaining genome stability. Despite having some particularities, such as its lower activity on transversions and the inability to correct indels, EndoMS/NucS meets the main hallmarks of a MMR. Its absence leads to a hypermutator phenotype, a transition-biased mutational spectrum and an increase in homeologous recombination. Interestingly, polymorphic EndoMS/NucS variants with a possible effect in mutation rate have been detected in clinical isolates of the relevant actinobacterial pathogen Mycobacterium tuberculosis. Considering that MMR defects are often associated with the emergence of resistant bacteria, the existence of EndoMS/NucS-defective mutators could have an important role in the acquisition of antibiotic resistance in M. tuberculosis. Therefore, a further understanding of the EndoMS/NucS-mediated non-canonical MMR pathway may reveal new strategies to predict and fight drug resistance. This review is focused on the recent progress in NucS, with special emphasis on its effect on genome stability and evolvability in Actinobacteria.

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.

Comparative Evolutionary Patterns of Burkholderia cenocepacia and B. multivorans During Chronic Co-infection of a Cystic Fibrosis Patient Lung

During chronic respiratory infections of cystic fibrosis (CF) patients, bacteria adaptively evolve in response to the nutritional and immune environment as well as influence other infecting microbes. The present study was designed to gain insights into the genetic mechanisms underlying adaptation and diversification by the two most prevalent pathogenic species of the Burkholderia cepacia complex (Bcc), B. cenocepacia and B. multivorans. Herein, we study the evolution of both of these species during coinfection of a CF patient for 4.4 years using genome sequences of 9 B. multivorans and 11 B. cenocepacia. This co-infection spanned at least 3 years following initial infection by B. multivorans and ultimately ended in the patient's death by cepacia syndrome. Both species acquired several mutations with accumulation rates of 2.08 (B. cenocepacia) and 2.27 (B. multivorans) SNPs/year. Many of the mutated genes are associated with oxidative stress response, transition metal metabolism, defense mechanisms against antibiotics, and other metabolic alterations consistent with the idea that positive selection might be driven by the action of the host immune system, antibiotic therapy and low oxygen and iron concentrations. Two orthologous genes shared by B. cenocepacia and B. multivorans were found to be under strong selection and accumulated mutations associated with lineage diversification. One gene encodes a nucleotide sugar dehydratase involved in lipopolysaccharide O-antigen (OAg) biosynthesis (wbiI). The other gene encodes a putative two-component regulatory sensor kinase protein required to sense and adapt to oxidative- and heavy metal- inducing stresses. This study contributes to understanding of shared and species-specific evolutionary patterns of B. cenocepacia and B. multivorans evolving in the same CF lung environment.

Diversity within species: interpreting strains in microbiomes

Studying within-species variation has traditionally been limited to culturable bacterial isolates and low-resolution microbial community fingerprinting. Metagenomic sequencing and technical advances have enabled culture-free, high-resolution strain and subspecies analyses at high throughput and in complex environments. This holds great scientific promise but has also led to an overwhelming number of methods and terms to describe infraspecific variation. This Review aims to clarify these advances by focusing on the diversity within bacterial and archaeal species in the context of microbiomics. We cover foundational microevolutionary concepts relevant to population genetics and summarize how within-species variation can be studied and stratified directly within microbial communities with a focus on metagenomics. Finally, we describe how common applications of within-species variation can be achieved using metagenomic data. We aim to guide the selection of appropriate terms and analytical approaches to facilitate researchers in benefiting from the increasing availability of large, high-resolution microbiome genetic sequencing data.

A Hopeful Sea-Monster: A Very Large Homologous Recombination Event Impacting the Core Genome of the Marine Pathogen Vibrio anguillarum

Vibrio anguillarum is the causative agent of vibriosis in many species important to aquaculture. We generated whole genome sequence (WGS) data on a diverse collection of 64 V. anguillarum strains, which we supplemented with 41 publicly available genomes to produce a combined dataset of 105 strains. These WGS data resolved six major lineages (L1-L6), and the additional use of multilocus sequence analysis (MLSA) clarified the association of L1 with serotype O1 and Salmonidae hosts (salmon/trout), and L2 with serotypes O2a/O2b/O2c and Gadidae hosts (cod). Our analysis also revealed a large-scale homologous replacement of 526-kb of core genome in an L2 strain from a con-specific donor. Although the strains affected by this recombination event are exclusively associated with Gadidae, we find no clear genetic evidence that it has played a causal role in host specialism. Whilst it is established that Vibrio species freely recombine, to our knowledge this is the first report of a contiguous recombinational replacement of this magnitude in any Vibrio genome. We also note a smaller accessory region of high single nucleotide polymorphism (SNP) density and gene content variation that contains lipopolysaccharide biosynthesis genes which may play a role in determining serotype.

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.

Using Colonization Assays and Comparative Genomics To Discover Symbiosis Behaviors and Factors in Vibrio fischeri

The luminous marine Gram-negative bacterium Vibrio (Aliivibrio) fischeri is the natural light organ symbiont of several squid species, including the Hawaiian bobtail squid, Euprymna scolopes, and the Japanese bobtail squid, Euprymna morsei Work with E. scolopes has shown how the bacteria establish their niche in the light organ of the newly hatched host. Two types of V. fischeri strains have been distinguished based upon their behavior in cocolonization competition assays in juvenile E. scolopes, i.e., (i) niche-sharing or (ii) niche-dominant behavior. This study aimed to determine whether these behaviors are observed with other V. fischeri strains or whether they are specific to those isolated from E. scolopes light organs. Cocolonization competition assays between V. fischeri strains isolated from the congeneric squid E. morsei or from other marine animals revealed the same sharing or dominant behaviors. In addition, whole-genome sequencing of these strains showed that the dominant behavior is polyphyletic and not associated with the presence or absence of a single gene or genes. Comparative genomics of 44 squid light organ isolates from around the globe led to the identification of symbiosis-specific candidates in the genomes of these strains. Colonization assays using genetic derivatives with deletions of these candidates established the importance of two such genes in colonization. This study has allowed us to expand the concept of distinct colonization behaviors to strains isolated from a number of squid and fish hosts.IMPORTANCE There is an increasing recognition of the importance of strain differences in the ecology of a symbiotic bacterial species and, in particular, how these differences underlie crucial interactions with their host. Nevertheless, little is known about the genetic bases for these differences, how they manifest themselves in specific behaviors, and their distribution among symbionts of different host species. In this study, we sequenced the genomes of Vibrio fischeri isolated from the tissues of squids and fishes and applied comparative genomics approaches to look for patterns between symbiont lineages and host colonization behavior. In addition, we identified the only two genes that were exclusively present in all V. fischeri strains isolated from the light organs of sepiolid squid species. Mutational studies of these genes indicated that they both played a role in colonization of the squid light organ, emphasizing the value of applying a comparative genomics approach in the study of symbioses.

Specificity and mutagenesis bias of the mycobacterial alternative mismatch repair analyzed by mutation accumulation studies

The postreplicative mismatch repair (MMR) is an almost ubiquitous DNA repair essential for maintaining genome stability. It has been suggested that Mycobacteria have an alternative MMR in which NucS, an endonuclease with no structural homology to the canonical MMR proteins (MutS/MutL), is the key factor. Here, we analyze the spontaneous mutations accumulated in a neutral manner over thousands of generations by Mycobacterium smegmatis and its MMR-deficient derivative (ΔnucS). The base pair substitution rates per genome per generation are 0.004 and 0.165 for wild type and ΔnucS, respectively. By comparing the activity of different bacterial MMR pathways, we demonstrate that both MutS/L- and NucS-based systems display similar specificity and mutagenesis bias, revealing a functional evolutionary convergence. However, NucS is not able to repair indels in vivo. Our results provide an unparalleled view of how this mycobacterial system works in vivo to maintain genome stability and how it may affect Mycobacterium evolution.

First Estimation of the Spontaneous Mutation Rate in Diatoms

Mutations are the origin of genetic diversity, and the mutation rate is a fundamental parameter to understand all aspects of molecular evolution. The combination of mutation-accumulation experiments and high-throughput sequencing enabled the estimation of mutation rates in most model organisms, but several major eukaryotic lineages remain unexplored. Here, we report the first estimation of the spontaneous mutation rate in a model unicellular eukaryote from the Stramenopile kingdom, the diatom Phaeodactylum tricornutum (strain RCC2967). We sequenced 36 mutation accumulation lines for an average of 181 generations per line and identified 156 de novo mutations. The base substitution mutation rate per site per generation is μbs = 4.77 × 10-10 and the insertion-deletion mutation rate is μid = 1.58 × 10-11. The mutation rate varies as a function of the nucleotide context and is biased toward an excess of mutations from GC to AT, consistent with previous observations in other species. Interestingly, the mutation rates between the genomes of organelles and the nucleus differ, with a significantly higher mutation rate in the mitochondria. This confirms previous claims based on indirect estimations of the mutation rate in mitochondria of photosynthetic eukaryotes that acquired their plastid through a secondary endosymbiosis. This novel estimate enables us to infer the effective population size of P. tricornutum to be Ne∼8.72 × 106.

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.

Molecular evolution of genes encoding allergen proteins in the peanuts genus Arachis: Structural and functional implications

Food allergies are severe immune responses to plant and animal products mediated by immunoglobulin E (IgE). Peanuts (Arachis hypogaea L.) are among the top 15 crops that feed the world. However, peanuts is among the "big eight food allergens", and allergies induced by peanuts are a significant public health problem and a life-threatening concern. Targeted mutation studies in peanuts demonstrate that single residue alterations in these allergen proteins could result in substantial reduction in allergenicity. Knowledge of peanut allergen proteins is confined to the allotetraploid crop and its two progenitors. We explored frequencies and positions of natural mutations in the hyperallergenic homologues Ara h 2 and Ara h 6 in newly generated sequences for 24 Arachis wild species and the crop species, assessed potential mutational impact on allergenicity using immunoblots and structural modeling, and evaluated whether these mutations follow evolutionary trends. We uncovered a wealth of natural mutations, both substitutions and gaps, including the elimination of immunodominant epitopes in some species. These molecular alterations appear to be associated with substantial reductions in allergenicity. The study demonstrated that Ara h 2 and Ara h 6 follow contrasting modes of natural selection and opposing mutational patterns, particularly in epitope regions. Phylogenetic analysis revealed a progressive trend towards immunodominant epitope evolution in Ara h 2. The findings provide valuable insight into the interactions among mutations, protein structure and immune system response, thus presenting a valuable platform for future manipulation of allergens to minimize, treat or eliminate allergenicity. The study strongly encourages exploration of genepools of economically important plants in allergenicity research.

Evolutionary pathways to antibiotic resistance are dependent upon environmental structure and bacterial lifestyle

Bacterial populations vary in their stress tolerance and population structure depending upon whether growth occurs in well-mixed or structured environments. We hypothesized that evolution in biofilms would generate greater genetic diversity than well-mixed environments and lead to different pathways of antibiotic resistance. We used experimental evolution and whole genome sequencing to test how the biofilm lifestyle influenced the rate, genetic mechanisms, and pleiotropic effects of resistance to ciprofloxacin in Acinetobacter baumannii populations. Both evolutionary dynamics and the identities of mutations differed between lifestyle. Planktonic populations experienced selective sweeps of mutations including the primary topoisomerase drug targets, whereas biofilm-adapted populations acquired mutations in regulators of efflux pumps. An overall trade-off between fitness and resistance level emerged, wherein biofilm-adapted clones were less resistant than planktonic but more fit in the absence of drug. However, biofilm populations developed collateral sensitivity to cephalosporins, demonstrating the clinical relevance of lifestyle on the evolution of resistance.

Selection, periodicity and potential function for Highly Iterative Palindrome-1 (HIP1) in cyanobacterial genomes

Highly Iterated Palindrome 1 (HIP1, GCGATCGC) is hyper-abundant in most cyanobacterial genomes. In some cyanobacteria, average HIP1 abundance exceeds one motif per gene. Such high abundance suggests a significant role in cyanobacterial biology. However, 20 years of study have not revealed whether HIP1 has a function, much less what that function might be. We show that HIP1 is 15- to 300-fold over-represented in genomes analyzed. More importantly, HIP1 sites are conserved both within and between open reading frames, suggesting that their overabundance is maintained by selection rather than by continual replenishment by neutral processes, such as biased DNA repair. This evidence for selection suggests a functional role for HIP1. No evidence was found to support a functional role as a peptide or RNA motif or a role in the regulation of gene expression. Rather, we demonstrate that the distribution of HIP1 along cyanobacterial chromosomes is significantly periodic, with periods ranging from 10 to 90 kb, consistent in scale with periodicities reported for co-regulated, co-expressed and evolutionarily correlated genes. The periodicity we observe is also comparable in scale to chromosomal interaction domains previously described in other bacteria. In this context, our findings imply HIP1 functions associated with chromosome and nucleoid structure.

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.

Chromids Aid Genome Expansion and Functional Diversification in the Family Burkholderiaceae

Multipartite genomes, containing at least two large replicons, are found in diverse bacteria; however, the advantage of this genome structure remains incompletely understood. Here, we perform comparative genomics of hundreds of finished β-proteobacterial genomes to gain insights into the role and emergence of multipartite genomes. Almost all essential secondary replicons (chromids) of the β-proteobacteria are found in the family Burkholderiaceae. These replicons arose from just two plasmid acquisition events, and they were likely stabilized early in their evolution by the presence of core genes. On average, Burkholderiaceae genera with multipartite genomes had a larger total genome size, but smaller chromosome, than genera without secondary replicons. Pangenome-level functional enrichment analyses suggested that interreplicon functional biases are partially driven by the enrichment of secondary replicons in the accessory pangenome fraction. Nevertheless, the small overlap in orthologous groups present in each replicon's pangenome indicated a clear functional separation of the replicons. Chromids appeared biased to environmental adaptation, as the functional categories enriched on chromids were also overrepresented on the chromosomes of the environmental genera (Paraburkholderia and Cupriavidus) compared with the pathogenic genera (Burkholderia and Ralstonia). Using ancestral state reconstruction, it was predicted that the rate of accumulation of modern-day genes by chromids was more rapid than the rate of gene accumulation by the chromosomes. Overall, the data are consistent with a model where the primary advantage of secondary replicons is in facilitating increased rates of gene acquisition through horizontal gene transfer, consequently resulting in replicons enriched in genes associated with adaptation to novel environments.

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.

Transition bias influences the evolution of antibiotic resistance in Mycobacterium tuberculosis

Transition bias, an overabundance of transitions relative to transversions, has been widely reported among studies of the rates and spectra of spontaneous mutations. However, demonstrating the role of transition bias in adaptive evolution remains challenging. In particular, it is unclear whether such biases direct the evolution of bacterial pathogens adapting to treatment. We addressed this challenge by analyzing adaptive antibiotic-resistance mutations in the major human pathogen Mycobacterium tuberculosis (MTB). We found strong evidence for transition bias in two independently curated data sets comprising 152 and 208 antibiotic-resistance mutations. This was true at the level of mutational paths (distinct adaptive DNA sequence changes) and events (individual instances of the adaptive DNA sequence changes) and across different genes and gene promoters conferring resistance to a diversity of antibiotics. It was also true for mutations that do not code for amino acid changes (in gene promoters and the 16S ribosomal RNA gene rrs) and for mutations that are synonymous to each other and are therefore likely to have similar fitness effects, suggesting that transition bias can be caused by a bias in mutation supply. These results point to a central role for transition bias in determining which mutations drive adaptive antibiotic resistance evolution in a key pathogen.

Some types of mutations occur more frequently than expected. This study shows that such bias —an excess of transitions over transversions—influences the evolution of antibiotic resistance in a key global pathogen, Mycobacterium tuberculosis.

Vibrio cholerae genomic diversity within and between patients

Cholera is a severe, water-borne diarrhoeal disease caused by toxin-producing strains of the bacterium Vibrio cholerae. Comparative genomics has revealed 'waves' of cholera transmission and evolution, in which clones are successively replaced over decades and centuries. However, the extent of V. cholerae genetic diversity within an epidemic or even within an individual patient is poorly understood. Here, we characterized V. cholerae genomic diversity at a micro-epidemiological level within and between individual patients from Bangladesh and Haiti. To capture within-patient diversity, we isolated multiple (8 to 20) V. cholerae colonies from each of eight patients, sequenced their genomes and identified point mutations and gene gain/loss events. We found limited but detectable diversity at the level of point mutations within hosts (zero to three single nucleotide variants within each patient), and comparatively higher gene content variation within hosts (at least one gain/loss event per patient, and up to 103 events in one patient). Much of the gene content variation appeared to be due to gain and loss of phage and plasmids within the V. cholerae population, with occasional exchanges between V. cholerae and other members of the gut microbiota. We also show that certain intra-host variants have phenotypic consequences. For example, the acquisition of a Bacteroides plasmid and non-synonymous mutations in a sensor histidine kinase gene both reduced biofilm formation, an important trait for environmental survival. Together, our results show that V. cholerae is measurably evolving within patients, with possible implications for disease outcomes and transmission dynamics.

Old Trade, New Tricks: Insights into the Spontaneous Mutation Process from the Partnering of Classical Mutation Accumulation Experiments with High-Throughput Genomic Approaches

Mutations spawn genetic variation which, in turn, fuels evolution. Hence, experimental investigations into the rate and fitness effects of spontaneous mutations are central to the study of evolution. Mutation accumulation (MA) experiments have served as a cornerstone for furthering our understanding of spontaneous mutations for four decades. In the pregenomic era, phenotypic measurements of fitness-related traits in MA lines were used to indirectly estimate key mutational parameters, such as the genomic mutation rate, new mutational variance per generation, and the average fitness effect of mutations. Rapidly emerging next-generating sequencing technology has supplanted this phenotype-dependent approach, enabling direct empirical estimates of the mutation rate and a more nuanced understanding of the relative contributions of different classes of mutations to the standing genetic variation. Whole-genome sequencing of MA lines bears immense potential to provide a unified account of the evolutionary process at multiple levels-the genetic basis of variation, and the evolutionary dynamics of mutations under the forces of selection and drift. In this review, we have attempted to synthesize key insights into the spontaneous mutation process that are rapidly emerging from the partnering of classical MA experiments with high-throughput sequencing, with particular emphasis on the spontaneous rates and molecular properties of different mutational classes in nuclear and mitochondrial genomes of diverse taxa, the contribution of mutations to the evolution of gene expression, and the rate and stability of transgenerational epigenetic modifications. Future advances in sequencing technologies will enable greater species representation to further refine our understanding of mutational parameters and their functional consequences.

Patterns of Nucleotide Deletion and Insertion Inferred from Bacterial Pseudogenes

Pseudogenes are a paradigm of neutral evolution and their study has the potential to reveal intrinsic mutational biases. However, this potential is mitigated by the fact that pseudogenes are quickly purged from bacterial genomes. Here, we assembled a large set of pseudogenes from genomes experiencing reductive evolution as well as functional references for which we could establish reliable phylogenetic relationships. Using this unique dataset, we identified 857 independent insertion and deletion mutations and discover a pervasive bias towards deletions, but not insertions, with sizes multiples of 3 nt. We further show that selective constraints for the preservation of gene frame are unlikely to account for the observed mutational bias and propose that a mechanistic bias in alternative end-joining repair, a recombination-independent double strand break DNA repair mechanism, is responsible for the accumulation of 3n deletions.

Insertion sequence elements-mediated structural variations in bacterial genomes

Mobile genetic elements (MGEs) impact the evolution and stability of their host genomes. Insertion sequence (IS) elements are the most common MGEs in bacterial genomes and play a crucial role in mediating large-scale variations in bacterial genomes. It is understood that IS elements and MGEs in general coexist in a dynamical equilibrium with their respective hosts. Current studies indicate that the spontaneous movement of IS elements does not follow a constant rate in different bacterial genomes. However, due to the paucity and sparsity of the data, these observations are yet to be conclusive. In this paper, we conducted a comparative analysis of the IS-mediated genome structural variations in ten mutation accumulation (MA) experiments across eight strains of five bacterial species containing IS elements, including four strains of the E. coli. We used GRASPER algorithm, a denovo structural variation (SV) identification algorithm designed to detect SVs involving repetitive sequences in the genome. We observed highly diverse rates of IS insertions and IS-mediated recombinations across different bacterial species as well as across different strains of the same bacterial species. We also observed different rates of the elements from the same IS family in different bacterial genomes, suggesting that the distinction in rates might not be due to the different composition of IS elements across bacterial genomes.

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.

Determinants of Base-Pair Substitution Patterns Revealed by Whole-Genome Sequencing of DNA Mismatch Repair Defective Escherichia coli

Mismatch repair (MMR) is a major contributor to replication fidelity, but its impact varies with sequence context and the nature of the mismatch. Mutation accumulation experiments followed by whole-genome sequencing of MMR-defective Escherichia coli strains yielded ≈30,000 base-pair substitutions (BPSs), revealing mutational patterns across the entire chromosome. The BPS spectrum was dominated by A:T to G:C transitions, which occurred predominantly at the center base of 5'NAC3'+5'GTN3' triplets. Surprisingly, growth on minimal medium or at low temperature attenuated these mutations. Mononucleotide runs were also hotspots for BPSs, and the rate at which these occurred increased with run length. Comparison with ≈2000 BPSs accumulated in MMR-proficient strains revealed that both kinds of hotspots appeared in the wild-type spectrum and so are likely to be sites of frequent replication errors. In MMR-defective strains transitions were strand biased, occurring twice as often when A and C rather than T and G were on the lagging-strand template. Loss of nucleotide diphosphate kinase increases the cellular concentration of dCTP, which resulted in increased rates of mutations due to misinsertion of C opposite A and T. In an mmr ndk double mutant strain, these mutations were more frequent when the template A and T were on the leading strand, suggesting that lagging-strand synthesis was more error-prone, or less well corrected by proofreading, than was leading strand synthesis.

Mutation Landscape of Base Substitutions, Duplications, and Deletions in the Representative Current Cholera Pandemic Strain

Pandemic cholera is a major concern for public health because of its high mortality and morbidity. Mutation accumulation (MA) experiments were performed on a representative strain of the current cholera pandemic. Although the base-pair substitution mutation rates in Vibrio cholerae (1.24 × 10-10 per site per generation for wild-type lines and 3.29 × 10-8 for mismatch repair deficient lines) are lower than that previously reported in other bacteria using MA analysis, we discovered specific high rates (8.31 × 10-8 site/generation for wild-type lines and 1.82 × 10-6 for mismatch repair deficient lines) of base duplication or deletion driven by large-scale copy number variations (CNVs). These duplication-deletions are located in two pathogenic islands, IMEX and the large integron island. Each element of these islands has discrepant rate in rapid integration and excision, which provides clues to the pandemicity evolution of V. cholerae. These results also suggest that large-scale structural variants such as CNVs can accumulate rapidly during short-term evolution. Mismatch repair deficient lines exhibit a significantly increased mutation rate in the larger chromosome (Chr1) at specific regions, and this pattern is not observed in wild-type lines. We propose that the high frequency of GATC sites in Chr1 improves the efficiency of MMR, resulting in similar rates of mutation in the wild-type condition. In addition, different mutation rates and spectra were observed in the MA lines under distinct growth conditions, including minimal media, rich media and antibiotic treatments.

Experimental Evolution as a High-Throughput Screen for Genetic Adaptations

Experimental evolution is a method in which populations of organisms, often microbes, are founded by one or more ancestors of known genotype and then propagated under controlled conditions to study the evolutionary process. These evolving populations are influenced by all population genetic forces, including selection, mutation, drift, and recombination, and the relative contributions of these forces may be seen as mysterious. Here, I describe why the outcomes of experimental evolution should be viewed with greater certainty because the force of selection typically dominates. Importantly, any mutant rising rapidly to high frequency in large populations must have acquired adaptive traits in the selective environment. Sequencing the genomes of these mutants can identify genes or pathways that contribute to an adaptation. I review the logic and simple mathematics why this evolve-and-resequence approach is a powerful way to find the mutations or mutation combinations that best increase fitness in any new environment.

Carbon limitation drives GC content evolution of a marine bacterium in an individual-based genome-scale model

An important unanswered question in evolutionary genomics is the source of considerable variation of genomic base composition (GC content) even among organisms that share one habitat. Evolution toward GC-poor genomes has been considered a major adaptive pathway in the oligotrophic ocean, but GC-rich bacteria are also prevalent and highly successful in this environment. We quantify the contribution of multiple factors to the change of genomic GC content of Ruegeria pomeroyi DSS-3, a representative and GC-rich member in the globally abundant Roseobacter clade, using an agent-based model. The model simulates 2 × 10⁸ cells, which allows random genetic drift to act in a realistic manner. Each cell has a whole genome subject to base-substitution mutation and recombination, which affect the carbon and nitrogen requirements of DNA and protein pools. Nonsynonymous changes can be functionally deleterious. Together, these factors affect the growth and fitness. Simulations show that experimentally determined mutation bias toward GC is not sufficient to build the GC-rich genome of DSS-3. While nitrogen availability has been repeatedly hypothesized to drive the evolution of GC content in marine bacterioplankton, our model instead predicts that DSS-3 and its ancestors have been evolving in environments primarily limited by carbon.

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.

The Divided Bacterial Genome: Structure, Function, and Evolution

Approximately 10% of bacterial genomes are split between two or more large DNA fragments, a genome architecture referred to as a multipartite genome. This multipartite organization is found in many important organisms, including plant symbionts, such as the nitrogen-fixing rhizobia, and plant, animal, and human pathogens, including the genera Brucella, Vibrio, and Burkholderia. The availability of many complete bacterial genome sequences means that we can now examine on a broad scale the characteristics of the different types of DNA molecules in a genome. Recent work has begun to shed light on the unique properties of each class of replicon, the unique functional role of chromosomal and nonchromosomal DNA molecules, and how the exploitation of novel niches may have driven the evolution of the multipartite genome. The aims of this review are to (i) outline the literature regarding bacterial genomes that are divided into multiple fragments, (ii) provide a meta-analysis of completed bacterial genomes from 1,708 species as a way of reviewing the abundant information present in these genome sequences, and (iii) provide an encompassing model to explain the evolution and function of the multipartite genome structure. This review covers, among other topics, salient genome terminology; mechanisms of multipartite genome formation; the phylogenetic distribution of multipartite genomes; how each part of a genome differs with respect to genomic signatures, genetic variability, and gene functional annotation; how each DNA molecule may interact; as well as the costs and benefits of this genome structure.