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

Homologous Recombination Shapes the Architecture and Evolution of Bacterial Genomes

Homologous recombination is a key evolutionary force that varies considerably across bacterial species. However, how the landscape of homologous recombination varies across genes and within individual genomes has only been studied in a few species. Here, we used Approximate Bayesian Computation to estimate the recombination rate along the genomes of 145 bacterial species. Our results show that homologous recombination varies greatly along bacterial genomes and shapes many aspects of genome architecture and evolution. The genomic landscape of recombination presents several key signatures: rates are highest near the origin of replication in most species, patterns of recombination generally appear symmetrical in both replichores (i.e. replicational halves of circular chromosomes) and most species have genomic hotpots of recombination. Furthermore, many closely related species share conserved landscapes of recombination across orthologs indicating that recombination landscapes are conserved over significant evolutionary distances. We show evidence that recombination drives the evolution of GC-content through increasing the effectiveness of selection and not through biased gene conversion, thereby contributing to an ongoing debate. Finally, we demonstrate that the rate of recombination varies across gene function and that many hotspots of recombination are associated with adaptive and mobile regions often encoding genes involved in pathogenicity.

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.

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'.

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.

Selection in a growing colony biases results of mutation accumulation experiments

Mutations provide the raw material for natural selection to act. Therefore, understanding the variety and relative frequency of different type of mutations is critical to understanding the nature of genetic diversity in a population. Mutation accumulation (MA) experiments have been used in this context to estimate parameters defining mutation rates, distribution of fitness effects (DFE), and spectrum of mutations. MA experiments can be performed with different effective population sizes. In MA experiments with bacteria, a single founder is grown to a size of a colony (~ 10⁸). It is assumed that natural selection plays a minimal role in dictating the dynamics of colony growth. In this work, we simulate colony growth via a mathematical model, and use our model to mimic an MA experiment. We demonstrate that selection ensures that, in an MA experiment, fraction of all mutations that are beneficial is over-represented by a factor of almost two, and that the distribution of fitness effects of beneficial and deleterious mutations are inaccurately captured in an MA experiment. Given this, the estimate of mutation rates from MA experiments is non-trivial. We then perform an MA experiment with 160 lines of E. coli, and show that due to the effect of selection in a growing colony, the size and sector of a colony from which the experiment is propagated impacts the results. Overall, we demonstrate that the results of MA experiments need to be revisited taking into account the action of selection in a growing colony.

Speed variations of bacterial replisomes

Replisomes are multi-protein complexes that replicate genomes with remarkable speed and accuracy. Despite their importance, their dynamics is poorly characterized, especially in vivo. In this paper, we present an approach to infer the replisome dynamics from the DNA abundance distribution measured in a growing bacterial population. Our method is sensitive enough to detect subtle variations of the replisome speed along the genome. As an application, we experimentally measured the DNA abundance distribution in Escherichia coli populations growing at different temperatures using deep sequencing. We find that the average replisome speed increases nearly five-fold between 17°C and 37°C. Further, we observe wave-like variations of the replisome speed along the genome. These variations correlate with previously observed variations of the mutation rate, suggesting a common dynamical origin. Our approach has the potential to elucidate replication dynamics in E. coli mutants and in other bacterial species.

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.

Temperature dependence of spontaneous mutation rates

Mutation is the source of genetic variation and the fundament of evolution. Temperature has long been suggested to have a direct impact on realized spontaneous mutation rates. If mutation rates vary in response to environmental conditions, such as the variation of the ambient temperature through space and time, they should no longer be described as species-specific constants. By combining mutation accumulation with whole-genome sequencing in a multicellular organism, we provide empirical support to reject the null hypothesis of a constant, temperature-independent mutation rate. Instead, mutation rates depended on temperature in a U-shaped manner with increasing rates toward both temperature extremes. This relation has important implications for mutation-dependent processes in molecular evolution, processes shaping the evolution of mutation rates, and even the evolution of biodiversity as such.

Divergent Evolution of Mutation Rates and Biases in the Long-Term Evolution Experiment with Escherichia coli

All organisms encode enzymes that replicate, maintain, pack, recombine, and repair their genetic material. For this reason, mutation rates and biases also evolve by mutation, variation, and natural selection. By examining metagenomic time series of the Lenski long-term evolution experiment (LTEE) with Escherichia coli (Good BH, McDonald MJ, Barrick JE, Lenski RE, Desai MM. 2017. The dynamics of molecular evolution over 60,000 generations. Nature 551(7678):45–50.), we find that local mutation rate variation has evolved during the LTEE. Each LTEE population has evolved idiosyncratic differences in their rates of point mutations, indels, and mobile element insertions, due to the fixation of various hypermutator and antimutator alleles. One LTEE population, called Ara+3, shows a strong, symmetric wave pattern in its density of point mutations, radiating from the origin of replication. This pattern is largely missing from the other LTEE populations, most of which evolved missense, indel, or structural mutations in topA, fis, and dusB—loci that all affect DNA topology. The distribution of mutations in those genes over time suggests epistasis and historical contingency in the evolution of DNA topology, which may have in turn affected local mutation rates. Overall, the replicate populations of the LTEE have largely diverged in their mutation rates and biases, even though they have adapted to identical abiotic conditions.

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.

Bacterial phenotypic heterogeneity in DNA repair and mutagenesis

Genetically identical cells frequently exhibit striking heterogeneity in various phenotypic traits such as their morphology, growth rate, or gene expression. Such non-genetic diversity can help clonal bacterial populations overcome transient environmental challenges without compromising genome stability, while genetic change is required for long-term heritable adaptation. At the heart of the balance between genome stability and plasticity are the DNA repair pathways that shield DNA from lesions and reverse errors arising from the imperfect DNA replication machinery. In principle, phenotypic heterogeneity in the expression and activity of DNA repair pathways can modulate mutation rates in single cells and thus be a source of heritable genetic diversity, effectively reversing the genotype-to-phenotype dogma. Long-standing evidence for mutation rate heterogeneity comes from genetics experiments on cell populations, which are now complemented by direct measurements on individual living cells. These measurements are increasingly performed using fluorescence microscopy with a temporal and spatial resolution that enables localising, tracking, and counting proteins with single-molecule sensitivity. In this review, we discuss which molecular processes lead to phenotypic heterogeneity in DNA repair and consider the potential consequences on genome stability and dynamics in bacteria. We further inspect these concepts in the context of DNA damage and mutation induced by antibiotics.

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.

The Location of Substitutions and Bacterial Genome Arrangements

Increasing evidence supports the notion that different regions of a genome have unique rates of molecular change. This variation is particularly evident in bacterial genomes where previous studies have reported gene expression and essentiality tend to decrease, whereas substitution rates usually increase with increasing distance from the origin of replication. Genomic reorganization such as rearrangements occur frequently in bacteria and allow for the introduction and restructuring of genetic content, creating gradients of molecular traits along genomes. Here, we explore the interplay of these phenomena by mapping substitutions to the genomes of Escherichia coli, Bacillus subtilis, Streptomyces, and Sinorhizobium meliloti, quantifying how many substitutions have occurred at each position in the genome. Preceding work indicates that substitution rate significantly increases with distance from the origin. Using a larger sample size and accounting for genome rearrangements through ancestral reconstruction, our analysis demonstrates that the correlation between the number of substitutions and the distance from the origin of replication is significant but small and inconsistent in direction. Some replicons had a significantly decreasing trend (E. coli and the chromosome of S. meliloti), whereas others showed the opposite significant trend (B. subtilis, Streptomyces, pSymA and pSymB in S. meliloti). dN, dS, and ω were examined across all genes and there was no significant correlation between those values and distance from the origin. This study highlights the impact that genomic rearrangements and location have on molecular trends in some bacteria, illustrating the importance of considering spatial trends in molecular evolutionary analysis. Assuming that molecular trends are exclusively in one direction can be problematic.

The roles of replication-transcription conflict in mutagenesis and evolution of genome organization

Replication-transcription conflicts promote mutagenesis and give rise to evolutionary signatures, with fundamental importance to genome stability ranging from bacteria to metastatic cancer cells. This review focuses on the interplay between replication-transcription conflicts and the evolution of gene directionality. In most bacteria, the majority of genes are encoded on the leading strand of replication such that their transcription is co-directional with the direction of DNA replication fork movement. This gene strand bias arises primarily due to negative selection against deleterious consequences of head-on replication-transcription conflict. However, many genes remain head-on. Can head-on orientation provide some benefit? We combine insights from both mechanistic and evolutionary studies, review published work, and analyze gene expression data to evaluate an emerging model that head-on genes are temporal targets for adaptive mutagenesis during stress. We highlight the alternative explanation that genes in the head-on orientation may simply be the result of genomic inversions and relaxed selection acting on nonessential genes. We seek to clarify how the mechanisms of replication-transcription conflict, in concert with other mutagenic mechanisms, balanced by natural selection, have shaped bacterial genome evolution.

Experimental Evolution In Vivo To Identify Selective Pressures during Pneumococcal Colonization

Experimental evolution is a powerful technique to understand how populations evolve from selective pressures imparted by the surrounding environment. With the advancement of whole-population genomic sequencing, it is possible to identify and track multiple contending genotypes associated with adaptations to specific selective pressures. This approach has been used repeatedly with model species in vitro, but only rarely in vivo Herein we report results of replicate experimentally evolved populations of Streptococcus pneumoniae propagated by repeated murine nasal colonization with the aim of identifying gene products under strong selection as well as the population genetic dynamics of infection cycles. Frameshift mutations in one gene, dltB, responsible for incorporation of d-alanine into teichoic acids on the bacterial surface, evolved repeatedly and swept to high frequency. Targeted deletions of dltB produced a fitness advantage during initial nasal colonization coupled with a corresponding fitness disadvantage in the lungs during pulmonary infection. The underlying mechanism behind the fitness trade-off between these two niches was found to be enhanced adherence to respiratory cells balanced by increased sensitivity to host-derived antimicrobial peptides, a finding recapitulated in the murine model. Additional mutations that are predicted to affect trace metal transport, central metabolism, and regulation of biofilm production and competence were also selected. These data indicate that experimental evolution can be applied to murine models of pathogenesis to gain insight into organism-specific tissue tropisms.IMPORTANCE Evolution is a powerful force that can be experimentally harnessed to gain insight into how populations evolve in response to selective pressures. Herein we tested the applicability of experimental evolutionary approaches to gain insight into how the major human pathogen Streptococcus pneumoniae responds to repeated colonization events using a murine model. These studies revealed the population dynamics of repeated colonization events and demonstrated that in vivo experimental evolution resulted in highly reproducible trajectories that reflect the environmental niche encountered during nasal colonization. Mutations impacting the surface charge of the bacteria were repeatedly selected during colonization and provided a fitness benefit in this niche that was counterbalanced by a corresponding fitness defect during lung infection. These data indicate that experimental evolution can be applied to models of pathogenesis to gain insight into organism-specific tissue tropisms.

Variability in the durability of CRISPR-Cas immunity

The durability of host resistance is challenged by the ability of pathogens to escape the defence of their hosts. Understanding the variability in the durability of host resistance is of paramount importance for designing more effective control strategies against infectious diseases. Here, we study the durability of various clustered regularly interspaced short palindromic repeats-Cas (CRISPR-Cas) alleles of the bacteria Streptococcus thermophilus against lytic phages. We found substantial variability in durability among different resistant bacteria. Since the escape of the phage is driven by a mutation in the phage sequence targeted by CRISPR-Cas, we explored the fitness costs associated with these escape mutations. We found that, on average, escape mutations decrease the fitness of the phage. Yet, the magnitude of this fitness cost does not predict the durability of CRISPR-Cas immunity. We contend that this variability in the durability of resistance may be because of variations in phage mutation rate or in the proportion of lethal mutations across the phage genome. These results have important implications on the coevolutionary dynamics between bacteria and phages and for the optimal deployment of resistance strategies against pathogens and pests. Understanding the durability of CRISPR-Cas immunity may also help develop more effective gene-drive strategies based on CRISPR-Cas9 technology. This article is part of a discussion meeting issue 'The ecology and evolution of prokaryotic CRISPR-Cas adaptive immune systems'.

Low mutational load and high mutation rate variation in gut commensal bacteria

Bacteria generally live in species-rich communities, such as the gut microbiota. Yet little is known about bacterial evolution in natural ecosystems. Here, we followed the long-term evolution of commensal Escherichia coli in the mouse gut. We observe the emergence of mutation rate polymorphism, ranging from wild-type levels to 1,000-fold higher. By combining experiments, whole-genome sequencing, and in silico simulations, we identify the molecular causes and explore the evolutionary conditions allowing these hypermutators to emerge and coexist within the microbiota. The hypermutator phenotype is caused by mutations in DNA polymerase III proofreading and catalytic subunits, which increase mutation rate by approximately 1,000-fold and stabilise hypermutator fitness, respectively. Strong mutation rate variation persists for >1,000 generations, with coexistence between lineages carrying 4 to >600 mutations. The in vivo molecular evolution pattern is consistent with fitness effects of deleterious mutations s_d ≤ 10⁻⁴/generation, assuming a constant effect or exponentially distributed effects with a constant mean. Such effects are lower than typical in vitro estimates, leading to a low mutational load, an inference that is observed in in vivo and in vitro competitions. Despite large numbers of deleterious mutations, we identify multiple beneficial mutations that do not reach fixation over long periods of time. This indicates that the dynamics of beneficial mutations are not shaped by constant positive Darwinian selection but could be explained by other evolutionary mechanisms that maintain genetic diversity. Thus, microbial evolution in the gut is likely characterised by partial sweeps of beneficial mutations combined with hitchhiking of slightly deleterious mutations, which take a long time to be purged because they impose a low mutational load. The combination of these two processes could allow for the long-term maintenance of intraspecies genetic diversity, including mutation rate polymorphism. These results are consistent with the pattern of genetic polymorphism that is emerging from metagenomics studies of the human gut microbiota, suggesting that we have identified key evolutionary processes shaping the genetic composition of this community.

Weak-effect deleterious mutations and negative frequency–dependent selection, acting on beneficial mutations, shape the dynamics of molecular evolution within the mouse gut microbiota.

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.

Integration Host Factor IHF facilitates homologous recombination and mutagenic processes in Pseudomonas putida

Nucleoid-associated proteins (NAPs) such as IHF, HU, Fis, and H-NS alter the topology of bound DNA and may thereby affect accessibility of DNA to repair and recombination processes. To examine this possibility, we investigated the effect of IHF on the frequency of homologous recombination (HR) and point mutations in soil bacterium Pseudomonas putida by using plasmidial and chromosomal assays. We observed positive effect of IHF on the frequency of HR, whereas this effect varied depending both on the chromosomal location of the HR target and the type of plasmid used in the assay. The occurrence of point mutations in plasmid was also facilitated by IHF, whereas in the chromosome the positive effect of IHF appeared only at certain DNA sequences and/or chromosomal positions. We did not observe any significant effects of IHF on the spectrum of mutations. However, despite of the presence or absence of IHF, different mutational hot spots appeared both in plasmid and in chromosome. Additionally, the frequency of frameshift mutations in the chromosome was also strongly affected by the location of the mutational target sequence. Taking together, our results indicate that IHF facilitates the occurrence of genetic changes in P. putida, whereas the location of the target sequence affects both the IHF-dependent and IHF-independent mechanisms.

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.

What is mutation? A chapter in the series: How microbes "jeopardize" the modern synthesis

Mutations drive evolution and were assumed to occur by chance: constantly, gradually, roughly uniformly in genomes, and without regard to environmental inputs, but this view is being revised by discoveries of molecular mechanisms of mutation in bacteria, now translated across the tree of life. These mechanisms reveal a picture of highly regulated mutagenesis, up-regulated temporally by stress responses and activated when cells/organisms are maladapted to their environments-when stressed-potentially accelerating adaptation. Mutation is also nonrandom in genomic space, with multiple simultaneous mutations falling in local clusters, which may allow concerted evolution-the multiple changes needed to adapt protein functions and protein machines encoded by linked genes. Molecular mechanisms of stress-inducible mutation change ideas about evolution and suggest different ways to model and address cancer development, infectious disease, and evolution generally.