The properties of spontaneous mutations in the opportunistic pathogen Pseudomonas aeruginosa

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

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

Hypermutability bypasses genetic constraints in SCV phenotypic switching in Pseudomonas aeruginosa biofilms

Biofilms are critical in the persistence of Pseudomonas aeruginosa infections, particularly in cystic fibrosis patients. This study explores the adaptive mechanisms behind the phenotypic switching between Small Colony Variants (SCVs) and revertant states in P. aeruginosa biofilms, emphasizing hypermutability due to Mismatch Repair System (MRS) deficiencies. Through experimental evolution and whole-genome sequencing, we show that both wild-type and mutator strains undergo parallel evolution by accumulating compensatory mutations in factors regulating intracellular c-di-GMP levels, particularly in the Wsp and Yfi systems. While wild-type strains face genetic constraints, mutator strains bypass these by accessing alternative genetic pathways regulating c-di-GMP and biofilm formation. This increased genetic accessibility, driven by higher mutation rates and specific mutational biases, supports sustained cycles of SCV conversion and reversion. Our findings underscore the crucial role of hypermutability in P. aeruginosa adaptation, with significant implications for managing persistent infections in clinical settings.

The Effect of the Presence and Absence of DNA Repair Genes on the Rate and Pattern of Mutation in Bacteria

Bacteria lose and gain repair genes as they evolve. Here, we investigate the consequences of gain and loss of 11 DNA repair genes across a broad range of bacteria. Using synonymous polymorphisms from bacteria and a set of 50 phylogenetically independent contrasts, we find no evidence that the presence or absence of these 11 genes affects either the overall level of diversity or the pattern of mutation. Using phylogenetic generalized linear squares yields a similar conclusion. It seems likely that the lack of an effect is due to variation in the genetic background and the environment which obscures any effects that the presence or absence of individual genes might have.

Beyond antibiotics: Emerging antivirulence strategies to combat Pseudomonas aeruginosa in cystic fibrosis

Pseudomonas aeruginosa is an opportunistic pathogen that poses a significant threat to individuals suffering from cystic fibrosis (CF). The pathogen is highly prevalent in CF individuals and is responsible for chronic infection, resulting in severe tissue damage and poor patient outcome. Prolonged antibiotic administration has led to the emergence of multidrug resistance in P. aeruginosa. In this direction, antivirulence strategies achieving targeted inhibition of bacterial virulence pathways, including quorum sensing, efflux pumps, lectins, and iron chelators, have been explored against CF isolates of P. aeruginosa. Hence, this review article presents a bird's eye view on the pulmonary infections involving P. aeruginosa in CF patients by laying emphasis on factors contributing to bacterial colonization, persistence, and disease progression along with the current line of therapeutics against P. aeruginosa in CF. We further collate scientific literature and discusses various antivirulence strategies that have been tested against P. aeruginosa isolates from CF patients.

The evolution of antimicrobial peptide resistance in Pseudomonas aeruginosa is severely constrained by random peptide mixtures

The prevalence of antibiotic-resistant pathogens has become a major threat to public health, requiring swift initiatives for discovering new strategies to control bacterial infections. Hence, antibiotic stewardship and rapid diagnostics, but also the development, and prudent use, of novel effective antimicrobial agents are paramount. Ideally, these agents should be less likely to select for resistance in pathogens than currently available conventional antimicrobials. The usage of antimicrobial peptides (AMPs), key components of the innate immune response, and combination therapies, have been proposed as strategies to diminish the emergence of resistance. Herein, we investigated whether newly developed random antimicrobial peptide mixtures (RPMs) can significantly reduce the risk of resistance evolution in vitro to that of single sequence AMPs, using the ESKAPE pathogen Pseudomonas aeruginosa (P. aeruginosa) as a model gram-negative bacterium. Infections of this pathogen are difficult to treat due the inherent resistance to many drug classes, enhanced by the capacity to form biofilms. P. aeruginosa was experimentally evolved in the presence of AMPs or RPMs, subsequentially assessing the extent of resistance evolution and cross-resistance/collateral sensitivity between treatments. Furthermore, the fitness costs of resistance on bacterial growth were studied and whole-genome sequencing used to investigate which mutations could be candidates for causing resistant phenotypes. Lastly, changes in the pharmacodynamics of the evolved bacterial strains were examined. Our findings suggest that using RPMs bears a much lower risk of resistance evolution compared to AMPs and mostly prevents cross-resistance development to other treatments, while maintaining (or even improving) drug sensitivity. This strengthens the case for using random cocktails of AMPs in favour of single AMPs, against which resistance evolved in vitro, providing an alternative to classic antibiotics worth pursuing.

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.

A role for the stringent response in ciprofloxacin resistance in Pseudomonas aeruginosa

Pseudomonas aeruginosa is a major cause of nosocomial infections and the leading cause of chronic lung infections in cystic fibrosis and chronic obstructive pulmonary disease patients. Antibiotic treatment remains challenging because P. aeruginosa is resistant to high concentrations of antibiotics and has a remarkable ability to acquire mutations conferring resistance to multiple groups of antimicrobial agents. Here we report that when P. aeruginosa is plated on ciprofloxacin (cipro) plates, the majority of cipro-resistant (ciproR) colonies observed at and after 48 h of incubation carry mutations in genes related to the Stringent Response (SR). Mutations in one of the major SR components, spoT, were present in approximately 40% of the ciproR isolates. Compared to the wild-type strain, most of these isolates had decreased growth rate, longer lag phase and altered intracellular ppGpp content. Also, 75% of all sequenced mutations were insertions and deletions, with short deletions being the most frequently occurring mutation type. We present evidence that most of the observed mutations are induced on the selective plates in a subpopulation of cells that are not instantly killed by cipro. Our results suggests that the SR may be an important contributor to antibiotic resistance acquisition in P. aeruginosa.

Within-Host Evolution of Bacterial Pathogens in Acute and Chronic Infection

Bacterial pathogens undergo remarkable adaptive change in response to the selective forces they encounter during host colonization and infection. Studies performed over the past few decades have demonstrated that many general evolutionary processes can be discerned during the course of host adaptation, including genetic diversification of lineages, clonal succession events, convergent evolution, and balanced fitness trade-offs. In some cases, elevated mutation rates resulting from mismatch repair or proofreading deficiencies accelerate evolution, and active mobile genetic elements or phages may facilitate genome plasticity. The host immune response provides another critical component of the fitness landscapes guiding adaptation, and selection operating on pathogens at this level may lead to immune evasion and the establishment of chronic infection. This review summarizes recent advances in this field, with a special focus on different forms of bacterial genome plasticity in the context of infection, and considers clinical consequences of adaptive changes for the host.

Mutational spectra are associated with bacterial niche

As observed in cancers, individual mutagens and defects in DNA repair create distinctive mutational signatures that combine to form context-specific spectra within cells. We reasoned that similar processes must occur in bacterial lineages, potentially allowing decomposition analysis to detect both disruption of DNA repair processes and exposure to niche-specific mutagens. Here we reconstruct mutational spectra for 84 clades from 31 diverse bacterial species and find distinct mutational patterns. We extract signatures driven by specific DNA repair defects using hypermutator lineages, and further deconvolute the spectra into multiple signatures operating within different clades. We show that these signatures are explained by both bacterial phylogeny and replication niche. By comparing mutational spectra of clades from different environmental and biological locations, we identify niche-associated mutational signatures, and then employ these signatures to infer the predominant replication niches for several clades where this was previously obscure. Our results show that mutational spectra may be associated with sites of bacterial replication when mutagen exposures differ, and can be used in these cases to infer transmission routes for established and emergent human bacterial pathogens.

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.

Distribution of mutation rates challenges evolutionary predictability

Natural selection is commonly assumed to act on extensive standing genetic variation. Yet, accumulating evidence highlights the role of mutational processes creating this genetic variation: to become evolutionarily successful, adaptive mutants must not only reach fixation, but also emerge in the first place, i.e. have a high enough mutation rate. Here, we use numerical simulations to investigate how mutational biases impact our ability to observe rare mutational pathways in the laboratory and to predict outcomes in experimental evolution. We show that unevenness in the rates at which mutational pathways produce adaptive mutants means that most experimental studies lack power to directly observe the full range of adaptive mutations. Modelling mutation rates as a distribution, we show that a substantially larger target size ensures that a pathway mutates more commonly. Therefore, we predict that commonly mutated pathways are conserved between closely related species, but not rarely mutated pathways. This approach formalizes our proposal that most mutations have a lower mutation rate than the average mutation rate measured experimentally. We suggest that the extent of genetic variation is overestimated when based on the average mutation rate.

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

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

A complex interplay of evolutionary forces continues to shape ancient co-occurring symbiont genomes

Many insects depend on ancient associations with intracellular bacteria for essential nutrition. The genomes of these bacteria are often highly reduced. Although drift is a major driver of symbiont evolution, other evolutionary forces continue to influence them. To understand how ongoing molecular evolution and gene loss shape symbiont genomes, we sequenced two of the most ancient symbionts known, Sulcia and Nasuia, from 20 Hawaiian Nesophrosyne leafhoppers. We leveraged the parallel divergence of Nesophrosyne lineages throughout Hawaii as a natural experimental framework. Sulcia and Nasuia experience ongoing-but divergent-gene loss, often in a convergent fashion. Although some genes are under relaxed selection, purifying and positive selection are also important drivers of genome evolution, particularly in maintaining certain nutritional and cellular functions. Our results further demonstrate that symbionts experience dramatically different evolutionary environments, as evidenced by the finding that Sulcia and Nasuia have one of the slowest and fastest rates of molecular evolution known.

Metabolic basis for the evolution of a common pathogenic Pseudomonas aeruginosa variant

Microbes frequently evolve in reproducible ways. Here, we show that differences in specific metabolic regulation rather than inter-strain interactions explain the frequent presence of lasR loss-of-function (LOF) mutations in the bacterial pathogen Pseudomonas aeruginosa. While LasR contributes to virulence through its role in quorum sensing, lasR mutants have been associated with more severe disease. A model based on the intrinsic growth kinetics for a wild type strain and its LasR- derivative, in combination with an experimental evolution based genetic screen and further genetics analyses, indicated that differences in metabolism were sufficient to explain the rise of these common mutant types. The evolution of LasR- lineages in laboratory and clinical isolates depended on activity of the two-component system CbrAB, which modulates substrate prioritization through the catabolite repression control pathway. LasR- lineages frequently arise in cystic fibrosis lung infections and their detection correlates with disease severity. Our analysis of bronchoalveolar lavage fluid metabolomes identified compounds that negatively correlate with lung function, and we show that these compounds support enhanced growth of LasR- cells in a CbrB-controlled manner. We propose that in vivo metabolomes contribute to pathogen evolution, which may influence the progression of disease and its treatment.

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.

The Efficacy of Using Combination Therapy against Multi-Drug and Extensively Drug-Resistant Pseudomonas aeruginosa in Clinical Settings

Pseudomonas aeruginosa is a Gram-negative bacterium which is capable of developing a high level of antibiotic resistance. It has been placed on the WHO's critical priority pathogen list and it is commonly found in ventilator-associated pneumonia infections, blood stream infections and other largely hospital-acquired illnesses. These infections are difficult to effectively treat due to their increasing antibiotic resistance and as such patients are often treated with antibiotic combination regimens.

METHODS

We conducted a systematic search with screening criteria using the Ovid search engine and the Embase, Ovid Medline, and APA PsycInfo databases.

RESULTS

It was found that in many cases the combination therapies were able to match or outperform the monotherapies and none performed noticeably worse than the monotherapies. However, the clinical studies were mostly small, only a few were prospective randomized clinical trials and statistical significance was lacking.

CONCLUSIONS

It was concluded that combination therapies have a place in the treatment of these highly resistant bacteria and, in some cases, there is some evidence to suggest that they provide a more effective treatment than monotherapies.

The nutritional environment is sufficient to select coexisting biofilm and quorum-sensing mutants of Pseudomonas aeruginosa

The evolution of bacterial populations during infections can be influenced by various factors including available nutrients, the immune system, and competing microbes, rendering it difficult to identify the specific forces that select on evolved traits. The genomes of Pseudomonas aeruginosa isolated from the airway of patients with cystic fibrosis (CF), for example, have revealed commonly mutated genes, but which phenotypes led to their prevalence is often uncertain. Here, we focus on effects of nutritional components of the CF airway on genetic adaptations by P. aeruginosa grown in either well-mixed (planktonic) or biofilm-associated conditions. After only 80 generations of experimental evolution in a simple medium with glucose, lactate, and amino acids, all planktonic populations diversified into lineages with mutated genes common to CF infections: morA, encoding a regulator of biofilm formation, or lasR, encoding a quorum sensing regulator that modulates the expression of virulence factors. Although mutated quorum sensing is often thought to be selected in vivo due to altered virulence phenotypes or social cheating, isolates with lasR mutations demonstrated increased fitness when grown alone and outcompeted the ancestral PA14 strain. Nonsynonymous SNPs in morA increased fitness in a nutrient concentration-dependent manner during planktonic growth and surprisingly also increased biofilm production. Populations propagated in biofilm conditions also acquired mutations in loci associated with chronic infections, including lasR and cyclic-di-GMP regulators roeA and wspF. These findings demonstrate that nutrient conditions and biofilm selection are sufficient to select mutants with problematic clinical phenotypes including increased biofilm and altered quorum sensing. Importance Pseudomonas aeruginosa produces dangerous chronic infections that are known for their rapid diversification and recalcitrance to treatment. We performed evolution experiments to identify adaptations selected by two specific aspects of the CF respiratory environment: nutrient levels and surface attachment. Propagation of P. aeruginosa in nutrients present within the CF airway was sufficient to drive diversification into subpopulations with identical mutations in regulators of biofilm and quorum sensing to those arising during infection. Thus, the adaptation of opportunistic pathogens to nutrients found in the host may select mutants with phenotypes that complicate treatment and clearance of infection.

Low base-substitution mutation rate and predominance of insertion-deletion events in the acidophilic bacterium Acidobacterium capsulatum

Analyses of spontaneous mutation have shown that total genome-wide mutation rates are quantitatively similar for most prokaryotic organisms. However, this view is mainly based on organisms that grow best around neutral pH values (6.0-8.0). In particular, the whole-genome mutation rate has not been determined for an acidophilic organism. Here, we have determined the genome-wide rate of spontaneous mutation in the acidophilic Acidobacterium capsulatum using a direct and unbiased method: a mutation-accumulation experiment followed by whole-genome sequencing. Evaluation of 69 mutation accumulation lines of A. capsulatum after an average of ~2900 cell divisions yielded a base-substitution mutation rate of 1.22 × 10-10 per site per generation or 4 × 10-4 per genome per generation, which is significantly lower than the consensus value (2.5-4.6 × 10-3) of mesothermophilic (~15-40°C) and neutrophilic (pH 6-8) prokaryotic organisms. However, the insertion-deletion rate (0.43 × 10-10 per site per generation) is high relative to the base-substitution mutation rate. Organisms with a similar effective population size and a similar expected effect of genetic drift should have similar mutation rates. Because selection operates on the total mutation rate, it is suggested that the relatively high insertion-deletion rate may be balanced by a low base-substitution rate in A. capsulatum, with selection operating on the total mutation rate.

Evolutionary responses to codon usage of horizontally transferred genes in Pseudomonas aeruginosa: gene retention, amelioration and compensatory evolution

Prokaryote genome evolution is characterized by the frequent gain of genes through horizontal gene transfer (HGT). For a gene, being horizontally transferred can represent a strong change in its genomic and physiological context. If the codon usage of a transferred gene deviates from that of the receiving organism, the fitness benefits it provides can be reduced due to a mismatch with the expression machinery. Consequently, transferred genes with a deviating codon usage can be selected against or elicit evolutionary responses that enhance their integration, such as gene amelioration and compensatory evolution. Within bacterial species, the extent and relative importance of these different mechanisms has never been considered altogether. In this study, a phylogeny-based method was used to investigate the occurrence of these different evolutionary responses in Pseudomonas aeruginosa. Selection on codon usage of genes acquired through HGT was observed over evolutionary time, with the overall codon usage converging towards that of the core genome. Gene amelioration, through the accumulation of synonymous mutations after HGT, did not seem to systematically affect transferred genes. This pattern therefore seemed to be mainly driven by selective retention of transferred genes with an initial codon usage similar to that of the core genes. Additionally, variation in the copy number of tRNA genes was often associated with the acquisition of genes for which the observed variation could enhance their expression. This provides evidence that compensatory evolution might be an important mechanism for the integration of horizontally transferred genes.

The Building Blocks of Antimicrobial Resistance in Pseudomonas aeruginosa: Implications for Current Resistance-Breaking Therapies

P. aeruginosa is classified as a priority one pathogen by the World Health Organisation, and new drugs are urgently needed, due to the emergence of multidrug-resistant (MDR) strains. Antimicrobial-resistant nosocomial pathogens such as P. aeruginosa pose unwavering and increasing threats. Antimicrobial stewardship has been a challenge during the COVID-19 pandemic, with a majority of those hospitalized with SARS-CoV2 infection given antibiotics as a safeguard against secondary bacterial infection. This increased usage, along with increased handling of sanitizers and disinfectants globally, may further accelerate the development and spread of cross-resistance to antibiotics. In addition, P. aeruginosa is the primary causative agent of morbidity and mortality in people with the life-shortening genetic disease cystic fibrosis (CF). Prolonged periods of selective pressure, associated with extended antibiotic treatment and the actions of host immune effectors, results in widespread adaptive and acquired resistance in P. aeruginosa found colonizing the lungs of people with CF. This review discusses the arsenal of resistance mechanisms utilized by P. aeruginosa, how these operate under high-stress environments such as the CF lung and how their interconnectedness can result in resistance to multiple antibiotic classes. Intrinsic, adaptive and acquired resistance mechanisms will be described, with a focus on how each layer of resistance can serve as a building block, contributing to multi-tiered resistance to antimicrobial activity. Recent progress in the development of anti-resistance adjuvant therapies, targeting one or more of these building blocks, should lead to novel strategies for combatting multidrug resistant P. aeruginosa. Anti-resistance adjuvant therapy holds great promise, not least because resistance against such therapeutics is predicted to be rare. The non-bactericidal nature of anti-resistance adjuvants reduce the selective pressures that drive resistance. Anti-resistance adjuvant therapy may also be advantageous in facilitating efficacious use of traditional antimicrobials, through enhanced penetration of the antibiotic into the bacterial cell. Promising anti-resistance adjuvant therapeutics and targets will be described, and key remaining challenges highlighted. As antimicrobial stewardship becomes more challenging in an era of emerging and re-emerging infectious diseases and global conflict, innovation in antibiotic adjuvant therapy can play an important role in extending the shelf-life of our existing antimicrobial therapeutic agents.

cognac: rapid generation of concatenated gene alignments for phylogenetic inference from large, bacterial whole genome sequencing datasets

Background

The quantity of genomic data is expanding at an increasing rate. Tools for phylogenetic analysis which scale to the quantity of available data are required. To address this need, we present cognac, a user-friendly software package to rapidly generate concatenated gene alignments for phylogenetic analysis.

Results

We illustrate that cognac is able to rapidly identify phylogenetic marker genes using a data driven approach and efficiently generate concatenated gene alignments for very large genomic datasets. To benchmark our tool, we generated core gene alignments for eight unique genera of bacteria, including a dataset of over 11,000 genomes from the genus Escherichia producing an alignment with 1353 genes, which was constructed in less than 17 h.

Conclusions

We demonstrate that cognac presents an efficient method for generating concatenated gene alignments for phylogenetic analysis. We have released cognac as an R package (https://github.com/rdcrawford/cognac) with customizable parameters for adaptation to diverse applications.

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.

An interbacterial DNA deaminase toxin directly mutagenizes surviving target populations

When bacterial cells come in contact, antagonism mediated by the delivery of toxins frequently ensues. The potential for such encounters to have long-term beneficial consequences in recipient cells has not been investigated. Here, we examined the effects of intoxication by DddA, a cytosine deaminase delivered via the type VI secretion system (T6SS) of Burkholderia cenocepacia. Despite its killing potential, we observed that several bacterial species resist DddA and instead accumulate mutations. These mutations can lead to the acquisition of antibiotic resistance, indicating that even in the absence of killing, interbacterial antagonism can have profound consequences on target populations. Investigation of additional toxins from the deaminase superfamily revealed that mutagenic activity is a common feature of these proteins, including a representative we show targets single-stranded DNA and displays a markedly divergent structure. Our findings suggest that a surprising consequence of antagonistic interactions between bacteria could be the promotion of adaptation via the action of directly mutagenic toxins.

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.

Mutational Pressure Drives Differential Genome Conservation in Two Bacterial Endosymbionts of Sap-Feeding Insects

Compared with free-living bacteria, endosymbionts of sap-feeding insects have tiny and rapidly evolving genomes. Increased genetic drift, high mutation rates, and relaxed selection associated with host control of key cellular functions all likely contribute to genome decay. Phylogenetic comparisons have revealed massive variation in endosymbiont evolutionary rate, but such methods make it difficult to partition the effects of mutation versus selection. For example, the ancestor of Auchenorrhynchan insects contained two obligate endosymbionts, Sulcia and a betaproteobacterium (BetaSymb; called Nasuia in leafhoppers) that exhibit divergent rates of sequence evolution and different propensities for loss and replacement in the ensuing ∼300 Ma. Here, we use the auchenorrhynchan leafhopper Macrosteles sp. nr. severini, which retains both of the ancestral endosymbionts, to test the hypothesis that differences in evolutionary rate are driven by differential mutagenesis. We used a high-fidelity technique known as duplex sequencing to measure and compare low-frequency variants in each endosymbiont. Our direct detection of de novode novo mutations reveals that the rapidly evolving endosymbiont (Nasuia) has a much higher frequency of single-nucleotide variants than the more stable endosymbiont (Sulcia) and a mutation spectrum that is potentially even more AT-biased than implied by the 83.1% AT content of its genome. We show that indels are common in both endosymbionts but differ substantially in length and distribution around repetitive regions. Our results suggest that differences in long-term rates of sequence evolution in Sulcia versus BetaSymb, and perhaps the contrasting degrees of stability of their relationships with the host, are driven by differences in mutagenesis.

A hydrophobic ratchet entrenches molecular complexes

Most proteins assemble into multisubunit complexes1. The persistence of these complexes across evolutionary time is usually explained as the result of natural selection for functional properties that depend on multimerization, such as intersubunit allostery or the capacity to do mechanical work2. In many complexes, however, multimerization does not enable any known function3. An alternative explanation is that multimers could become entrenched if substitutions accumulate that are neutral in multimers but deleterious in monomers; purifying selection would then prevent reversion to the unassembled form, even if assembly per se does not enhance biological function3-7. Here we show that a hydrophobic mutational ratchet systematically entrenches molecular complexes. By applying ancestral protein reconstruction and biochemical assays to the evolution of steroid hormone receptors, we show that an ancient hydrophobic interface, conserved for hundreds of millions of years, is entrenched because exposure of this interface to solvent reduces protein stability and causes aggregation, even though the interface makes no detectable contribution to function. Using structural bioinformatics, we show that a universal mutational propensity drives sites that are buried in multimeric interfaces to accumulate hydrophobic substitutions to levels that are not tolerated in monomers. In a database of hundreds of families of multimers, most show signatures of long-term hydrophobic entrenchment. It is therefore likely that many protein complexes persist because a simple ratchet-like mechanism entrenches them across evolutionary time, even when they are functionally gratuitous.

A model of autowave self-organization as a hierarchy of active media in the biological evolution

Within the framework of the active media concept, we develop a biophysical model of autowave self-organization which is treated as a hierarchy of active media in the evolution of the biosphere. We also propose a mathematical model of the autowave process of speciation in a flow of mutations for the three main taxonometric groups (prokaryotes, unicellular and multicellular eukaryotes) with a naturally determined lower boundary of living matter (the appearance of prokaryotes) and an open upper boundary for the formation of new species. It is shown that the fluctuation-bifurcation description of the evolution for the formation of new taxonometric groups as a trajectory of transformation of small fluctuations into giant ones adequately reflects the process of self-organization during the formation of taxa. The major concepts of biological evolution, conditions of hierarchy formation as a fundamental manifestation of self-organization and complexity in the evolution of biological systems are considered.

Mutational analysis of the essential lipopolysaccharide-transport protein LptH of Pseudomonas aeruginosa to uncover critical oligomerization sites

Lipopolysaccharide (LPS) is a critical component of the outer membrane (OM) of many Gram-negative bacteria. LPS is translocated to the OM by the LPS transport (Lpt) system. In the human pathogen Pseudomonas aeruginosa, the periplasmic Lpt component, LptH, is essential for LPS transport, planktonic and biofilm growth, OM stability and infectivity. LptH has been proposed to oligomerize and form a protein bridge that accommodates LPS during transport. Based on the known LptH crystal structure, here we predicted by in silico modeling five different sites likely involved in LptH oligomerization. The relevance of these sites for LptH activity was verified through plasmid-mediated expression of site-specific mutant proteins in a P. aeruginosa lptH conditional mutant. Complementation and protein expression analyses provided evidence that all mutated sites are important for LptH activity in vivo. It was observed that the lptH conditional mutant overcomes the lethality of nonfunctional lptH variants through RecA-mediated homologous recombination between the wild-type lptH gene in the genome and mutated copies in the plasmid. Finally, biochemical assays on purified recombinant proteins showed that some LptH variants are indeed specifically impaired in oligomerization, while others appear to have defects in protein folding and/or stability.

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.

Parallel Evolution of Tobramycin Resistance across Species and Environments

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

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

Genes and Proteomes Associated With Increased Mutation Frequency and Multidrug Resistance of Naturally Occurring Mismatch Repair-Deficient Salmonella Hypermutators

The emergence of antibiotic-resistant Salmonella through mutations led to mismatch repair (MMR) deficiency that represents a potential hazard to public health. Here, four representative MMR-deficient Salmonella hypermutator strains and Salmonella Typhimurium LT2 were used to comprehensively reveal the influence of MMR deficiency on antibiotic resistance among Salmonella. Our results indicated that the mutation frequency ranged from 3.39 × 10^–4 to 5.46 × 10^–2 in the hypermutator. Mutation sites in MutS, MutL, MutT, and UvrD of the four hypermutators were all located in the essential and core functional regions. Mutation frequency of the hypermutator was most highly correlated with the extent of mutation in MutS. Mutations in MMR genes (mutS, mutT, mutL, and uvrD) were correlated with increased mutation in antibiotic resistance genes, and the extent of antibiotic resistance was significantly correlated with the number of mutation sites in MutL and in ParC. The number of mutation sites in MMR genes and antibiotic resistance genes exhibited a significant positive correlation with the number of antibiotics resisted and with expression levels of mutS, mutT, and mutL. Compared to Salmonella Typhimurium LT2, a total of 137 differentially expressed and 110 specifically expressed proteins were identified in the four hypermutators. Functional enrichment analysis indicated that the proteins significantly overexpressed in the hypermutators primarily associated with translation and stress response. Interaction network analysis revealed that the ribosome pathway might be a critical factor for high mutation frequency and multidrug resistance in MMR-deficient Salmonella hypermutators. These results help elucidate the mutational dynamics that lead to hypermutation, antibiotic resistance, and activation of stress response pathways in Salmonella.

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.

A Continuous-Flow Model for in vitro Cultivation of Mixed Microbial Populations Associated With Cystic Fibrosis Airway Infections

The airways of people with cystic fibrosis (CF) provide a nutrient-rich environment which favours colonisation by a variety of bacteria and fungi. Although the dominant pathogen associated with CF airway infections is Pseudomonas aeruginosa, it is becoming increasingly clear that inter-species interactions between P. aeruginosa and other colonists in the airways may have a large impact on microbial physiology and virulence. However, there are currently no suitable experimental models that permit long-term co-culture of P. aeruginosa with other CF-associated pathogens. Here, we redress this problem by describing a "3R's-compliant" continuous-flow in vitro culture model which enables long-term co-culture of three representative CF-associated microbes: P. aeruginosa, Staphylococcus aureus and Candida albicans. Although these species rapidly out-compete one another when grown together or in pairs in batch culture, we show that in a continuously-fed setup, they can be maintained in a very stable, steady-state community. We use our system to show that even numerically (0.1%) minor species can have a major impact on intercellular signalling by P. aeruginosa. Importantly, we also show that co-culturing does not appear to influence species mutation rates, further reinforcing the notion that the system favours stability rather than divergence. The model is experimentally tractable and offers an inexpensive yet robust means of investigating inter-species interactions between CF pathogens.

Dynamic Emergence of Mismatch Repair Deficiency Facilitates Rapid Evolution of Ceftazidime-Avibactam Resistance in Pseudomonas aeruginosa Acute Infection

Strains of Pseudomonas aeruginosa with deficiencies in DNA mismatch repair have been studied in the context of chronic infection, where elevated mutational rates ("hypermutation") may facilitate the acquisition of antimicrobial resistance. Whether P. aeruginosa hypermutation can also play an adaptive role in the more dynamic context of acute infection remains unclear. In this work, we demonstrate that evolved mismatch repair deficiencies may be exploited by P. aeruginosa to facilitate rapid acquisition of antimicrobial resistance in acute infection, and we directly document rapid clonal succession by such a hypermutating lineage in a patient. Whole-genome sequencing (WGS) was performed on nine serially cultured blood and respiratory isolates from a patient in whom ceftazidime-avibactam (CZA) resistance emerged in vivo over the course of days. The CZA-resistant clone was differentiated by 14 mutations, including a gain-of-function G183D substitution in the PDC-5 chromosomal AmpC cephalosporinase conferring CZA resistance. This lineage also contained a substitution (R656H) at a conserved position in the ATPase domain of the MutS mismatch repair (MMR) protein, and elevated mutational rates were confirmed by mutational accumulation experiments with WGS of evolved lineages in conjunction with rifampin resistance assays. To test whether MMR-deficient hypermutation could facilitate rapid acquisition of CZA resistance, in vitro adaptive evolution experiments were performed with a mutS-deficient strain. These experiments demonstrated rapid hypermutation-facilitated acquisition of CZA resistance compared with the isogenic wild-type strain. Our results suggest a possibly underappreciated role for evolved MMR deficiency in facilitating rapid adaptive evolution of P. aeruginosa in the context of acute infection.IMPORTANCE Antimicrobial resistance in bacteria represents one of the most consequential problems in modern medicine, and its emergence and spread threaten to compromise central advances in the treatment of infectious diseases. Ceftazidime-avibactam (CZA) belongs to a new class of broad-spectrum beta-lactam/beta-lactamase inhibitor combinations designed to treat infections caused by multidrug-resistant bacteria. Understanding the emergence of resistance to this important new drug class is of critical importance. In this work, we demonstrate that evolved mismatch repair deficiency in P. aeruginosa, an important pathogen responsible for significant morbidity and mortality among hospitalized patients, may facilitate rapid acquisition of resistance to CZA in the context of acute infection. These findings are relevant for both diagnosis and treatment of antimicrobial resistance emerging in acute infection in the hypermutator background and additionally have implications for the emergence of more virulent phenotypes.

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.

The Essential Role of Hypermutation in Rapid Adaptation to Antibiotic Stress

A common outcome of antibiotic exposure in patients and in vitro is the evolution of a hypermutator phenotype that enables rapid adaptation by pathogens. While hypermutation is a robust mechanism for rapid adaptation, it requires trade-offs between the adaptive mutations and the more common "hitchhiker" mutations that accumulate from the increased mutation rate. Using quantitative experimental evolution, we examined the role of hypermutation in driving the adaptation of Pseudomonas aeruginosa to colistin. Metagenomic deep sequencing revealed 2,657 mutations at ≥5% frequency in 1,197 genes and 761 mutations in 29 endpoint isolates. By combining genomic information, phylogenetic analyses, and statistical tests, we showed that evolutionary trajectories leading to resistance could be reliably discerned. In addition to known alleles such as pmrB, hypermutation allowed identification of additional adaptive alleles with epistatic relationships. Although hypermutation provided a short-term fitness benefit, it was detrimental to overall fitness. Alarmingly, a small fraction of the colistin-adapted population remained colistin susceptible and escaped hypermutation. In a clinical population, such cells could play a role in reestablishing infection upon withdrawal of colistin. We present here a framework for evaluating the complex evolutionary trajectories of hypermutators that applies to both current and emerging pathogen populations.

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.

Experimental Evolution of Innovation and Novelty

How does novelty, a new, genetically based function, evolve? A compelling answer has been elusive because there are few model systems where both the genetic mechanisms generating novel functions and the ecological conditions that govern their origin and spread can be studied in detail. This review article considers what we have learned about the evolution of novelty from microbial selection experiments. This work reveals that the genetic routes to novelty can be more highly variable than standard models have led us to believe and underscores the importance of considering both genetics and ecology in this process.

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.

Genetic Drift and Indel Mutation in the Evolution of Yeast Mitochondrial Genome Size

Mitochondrial genomes (mitogenomes) are remarkably diverse in genome size and organization, but the origins of dynamic mitogenome architectures are still poorly understood. For instance, the mutational burden hypothesis postulates that the drastic difference between large plant mitogenomes and streamlined animal mitogenomes can be driven by their different mutation rates. However, inconsistent trends between mitogenome sizes and mutation rates have been documented in several lineages. These conflicting results highlight the need of systematic and sophisticated investigations on the evolution and diversity of mitogenome architecture. This study took advantage of the strikingly variable mitogenome size among different yeast species and also among intraspecific strains, examined sequence dynamics of introns, GC-clusters, tandem repeats, mononucleotide repeats (homopolymers) and evaluated their contributions to genome size variation. The contributions of these sequence features to mitogenomic variation are dependent on the timescale, over which extant genomes evolved from their last common ancestor, perhaps due to a combination of different turnover rates of mobile sequences, variable insertion spaces, and functional constraints. We observed a positive correlation between mitogenome size and the level of genetic drift, suggesting that mitogenome expansion in yeast is likely driven by multiple types of sequence insertions in a primarily nonadaptive manner. Although these cannot be explained directly by the mutational burden hypothesis, our results support an important role of genetic drift in the evolution of yeast mitogenomes.

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.

Spatial Vulnerabilities of the Escherichia coli Genome to Spontaneous Mutations Revealed with Improved Duplex Sequencing

Investigation of spontaneous mutations by next-generation sequencing technology has attracted extensive attention lately due to the fundamental roles of spontaneous mutations in evolution and pathological processes. However, these studies only focused on the mutations accumulated through many generations during long-term (possibly be years of) culturing, but not the freshly generated mutations that occur at very low frequencies. In this study, we established a molecularly barcoded deep sequencing strategy to detect low abundant spontaneous mutations in genomes of bacteria cell cultures. Genome-wide spontaneous mutations in 15 Escherichia coli cell culture samples were defined with a high confidence (P < 0.01). We also developed a hotspot-calling approach based on the run-length encoding algorithm to find the genomic regions that are vulnerable to the spontaneous mutations. The hotspots for the mutations appeared to be highly conserved across the bacteria samples. Further biological annotation of these regions indicated that most of the spontaneous mutations were located at the repeat domains or nonfunctional domains of the genomes, suggesting the existence of mechanisms that could somehow prevent the occurrence of mutations in crucial genic areas. This study provides a more faithful picture of mutation occurrence and spectra in a single expansion process without long-term culturing.

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.

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

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

The distribution of bacterial doubling times in the wild

Generation time varies widely across organisms and is an important factor in the life cycle, life history and evolution of organisms. Although the doubling time (DT) has been estimated for many bacteria in the laboratory, it is nearly impossible to directly measure it in the natural environment. However, an estimate can be obtained by measuring the rate at which bacteria accumulate mutations per year in the wild and the rate at which they mutate per generation in the laboratory. If we assume the mutation rate per generation is the same in the wild and in the laboratory, and that all mutations in the wild are neutral, an assumption that we show is not very important, then an estimate of the DT can be obtained by dividing the latter by the former. We estimate the DT for five species of bacteria for which we have both an accumulation and a mutation rate estimate. We also infer the distribution of DTs across all bacteria from the distribution of the accumulation and mutation rates. Both analyses suggest that DTs for bacteria in the wild are substantially greater than those in the laboratory, that they vary by orders of magnitude between different species of bacteria and that a substantial fraction of bacteria double very slowly in the wild.

Heterogeneous Antimicrobial Susceptibility Characteristics in Pseudomonas aeruginosa Isolates from Cystic Fibrosis Patients

Patients with cystic fibrosis endure “chronic focal infections” with a variety of microorganisms. One microorganism, Pseudomonas aeruginosa, adapts to the host and develops resistance to a wide range of antimicrobials. Interestingly, as the infection progresses, multiple isogenic strains of P. aeruginosa emerge and coexist within the airways of these patients. Despite a common parental origin, the multiple strains of P. aeruginosa develop vastly different susceptibility patterns to actively used antimicrobial agents—a phenomenon we define as “heterogeneous MICs.” By sequencing pairs of P. aeruginosa isolates displaying heterogeneous MICs, we observed widespread isogenic gene lesions in drug transporters, DNA mismatch repair mechanisms, and many other structural or cellular functions. Coupled with the heterogeneous MICs, these genetic lesions demonstrated a symbiotic response to host selection and suggested evolution of a multicellular syntrophic bacterial lifestyle. Current laboratory standard interpretive criteria do not address the emergence of heterogeneous growth and susceptibilities in vitro with treatment implications.

Clinical isolates of Pseudomonas aeruginosa from patients with cystic fibrosis (CF) are known to differ from those associated with non-CF hosts by colony morphology, drug susceptibility patterns, and genomic hypermutability. Pseudomonas aeruginosa isolates from CF patients have long been recognized for their overall reduced rate of antimicrobial susceptibility, but their intraclonal MIC heterogeneity has long been overlooked. Using two distinct cohorts of clinical strains (n = 224 from 56 CF patients, n = 130 from 68 non-CF patients) isolated in 2013, we demonstrated profound Etest MIC heterogeneity in CF P. aeruginosa isolates in comparison to non-CF P. aeruginosa isolates. On the basis of whole-genome sequencing of 19 CF P. aeruginosa isolates from 9 patients with heterogeneous MICs, the core genome phylogenetic tree confirmed the within-patient CF P. aeruginosa clonal lineage along with considerable coding sequence variability. No extrachromosomal DNA elements or previously characterized antibiotic resistance mutations could account for the wide divergence in antimicrobial MICs between P. aeruginosa coisolates, though many heterogeneous mutations in efflux and porin genes and their regulators were present. A unique OprD sequence was conserved among the majority of isolates of CF P. aeruginosa analyzed, suggesting a pseudomonal response to selective pressure that is common to the isolates. Genomic sequence data also suggested that CF pseudomonal hypermutability was not entirely due to mutations in mutL, mutS, and uvr. We conclude that the net effect of hundreds of adaptive mutations, both shared between clonally related isolate pairs and unshared, accounts for their highly heterogeneous MIC variances. We hypothesize that this heterogeneity is indicative of the pseudomonal syntrophic-like lifestyle under conditions of being “locked” inside a host focal airway environment for prolonged periods.

IMPORTANCE Patients with cystic fibrosis endure “chronic focal infections” with a variety of microorganisms. One microorganism, Pseudomonas aeruginosa, adapts to the host and develops resistance to a wide range of antimicrobials. Interestingly, as the infection progresses, multiple isogenic strains of P. aeruginosa emerge and coexist within the airways of these patients. Despite a common parental origin, the multiple strains of P. aeruginosa develop vastly different susceptibility patterns to actively used antimicrobial agents—a phenomenon we define as “heterogeneous MICs.” By sequencing pairs of P. aeruginosa isolates displaying heterogeneous MICs, we observed widespread isogenic gene lesions in drug transporters, DNA mismatch repair machinery, and many other structural or cellular functions. Coupled with the heterogeneous MICs, these genetic lesions demonstrated a symbiotic response to host selection and suggested evolution of a multicellular syntrophic bacterial lifestyle. Current laboratory standard interpretive criteria do not address the emergence of heterogeneous growth and susceptibilities in vitro with treatment implications.

Mutational Biases Influence Parallel Adaptation

While mutational biases strongly influence neutral molecular evolution, the role of mutational biases in shaping the course of adaptation is less clear. Here we consider the frequency of transitions relative to transversions among adaptive substitutions. Because mutation rates for transitions are higher than those for transversions, if mutational biases influence the dynamics of adaptation, then transitions should be overrepresented among documented adaptive substitutions. To test this hypothesis, we assembled two sets of data on putatively adaptive amino acid replacements that have occurred in parallel during evolution, either in nature or in the laboratory. We find that the frequency of transitions in these data sets is much higher than would be predicted under a null model where mutation has no effect. Our results are qualitatively similar even if we restrict ourself to changes that have occurred, not merely twice, but three or more times. These results suggest that the course of adaptation is biased by mutation.