Detecting rare structural variation in evolving microbial populations from new sequence junctions using breseq

Experimental evolution for the recovery of growth loss due to genome reduction

As the genome encodes the information crucial for cell growth, a sizeable genomic deficiency often causes a significant decrease in growth fitness. Whether and how the decreased growth fitness caused by genome reduction could be compensated by evolution was investigated here. Experimental evolution with an Escherichia coli strain carrying a reduced genome was conducted in multiple lineages for approximately 1000 generations. The growth rate, which largely declined due to genome reduction, was considerably recovered, associated with the improved carrying capacity. Genome mutations accumulated during evolution were significantly varied across the evolutionary lineages and were randomly localized on the reduced genome. Transcriptome reorganization showed a common evolutionary direction and conserved the chromosomal periodicity, regardless of highly diversified gene categories, regulons, and pathways enriched in the differentially expressed genes. Genome mutations and transcriptome reorganization caused by evolution, which were found to be dissimilar to those caused by genome reduction, must have followed divergent mechanisms in individual evolutionary lineages. Gene network reconstruction successfully identified three gene modules functionally differentiated, which were responsible for the evolutionary changes of the reduced genome in growth fitness, genome mutation, and gene expression, respectively. The diversity in evolutionary approaches improved the growth fitness associated with the homeostatic transcriptome architecture as if the evolutionary compensation for genome reduction was like all roads leading to Rome.

Large-scale duplication events underpin population-level flexibility in tRNA gene copy number in Pseudomonas fluorescens SBW25

The complement of tRNA genes within a genome is typically considered to be a (relatively) stable characteristic of an organism. Here, we demonstrate that bacterial tRNA gene set composition can be more flexible than previously appreciated, particularly regarding tRNA gene copy number. We report the high-rate occurrence of spontaneous, large-scale, tandem duplication events in laboratory populations of the bacterium Pseudomonas fluorescens SBW25. The identified duplications are up to ∼1 Mb in size (∼15% of the wildtype genome) and are predicted to change the copy number of up to 917 genes, including several tRNA genes. The observed duplications are inherently unstable: they occur, and are subsequently lost, at extremely high rates. We propose that this unusually plastic type of mutation provides a mechanism by which tRNA gene set diversity can be rapidly generated, while simultaneously preserving the underlying tRNA gene set in the absence of continued selection. That is, if a tRNA set variant provides no fitness advantage, then high-rate segregation of the duplication ensures the maintenance of the original tRNA gene set. However, if a tRNA gene set variant is beneficial, the underlying duplication fragment(s) may persist for longer and provide raw material for further, more stable, evolutionary change.

Flagella-related gene mutations in Vibrio cholerae during extended cultivation in nutrient-limited media impair cell motility and prolong culturability

IMPORTANCE

Vibrio cholerae undergoes a transition to a viable but non-culturable (VNC) state when subjected to various environmental stresses. We showed here that flagellar motility was involved in the development of the VNC state of V. cholerae. In this study, motility-defective isolates with mutations in various flagella-related genes, but not motile isolates, were predominantly obtained under the stress of long-term batch culture. Other genomic regions were highly conserved, suggesting that the mutations were selective. During the stationary phase of long-term culture, V. cholerae isolates with mutations in the acetate kinase and flagella-related genes were predominant. This study suggests that genes involved in specific functions in V. cholerae undergo mutations under certain environmental conditions.

Within-host adaptation alters priority effects within the tomato phyllosphere microbiome

To predict the composition and function of ecological communities over time, it is essential to understand how in situ evolution alters priority effects between resident and invading species. Phyllosphere microbial communities are a useful model system to explore priority effects because the system is clearly spatially delineated and can be manipulated experimentally. We conducted an experimental evolution study with tomato plants and the early-colonizing bacterium species Pantoea dispersa, exploring priority effects when P. dispersa was introduced before, simultaneously with or after competitor species. P. dispersa rapidly evolved to invade a new niche within the plant tissue and altered its ecological interactions with other members of the plant microbiome and its effect on the host. Prevailing models have assumed that adaptation primarily improves the efficiency of resident species within their existing niches, yet in our study system, the resident species expanded its niche instead. This finding suggests potential limitations to the application of existing ecological theory to microbial communities.

The Hypersaline Archaeal Histones HpyA and HstA Are DNA Binding Proteins That Defy Categorization According to Commonly Used Functional Criteria

Histone proteins are found across diverse lineages of Archaea, many of which package DNA and form chromatin. However, previous research has led to the hypothesis that the histone-like proteins of high-salt-adapted archaea, or halophiles, function differently. The sole histone protein encoded by the model halophilic species Halobacterium salinarum, HpyA, is nonessential and expressed at levels too low to enable genome-wide DNA packaging. Instead, HpyA mediates the transcriptional response to salt stress. Here we compare the features of genome-wide binding of HpyA to those of HstA, the sole histone of another model halophile, Haloferax volcanii. hstA, like hpyA, is a nonessential gene. To better understand HpyA and HstA functions, protein-DNA binding data (chromatin immunoprecipitation sequencing [ChIP-seq]) of these halophilic histones are compared to publicly available ChIP-seq data from DNA binding proteins across all domains of life, including transcription factors (TFs), nucleoid-associated proteins (NAPs), and histones. These analyses demonstrate that HpyA and HstA bind the genome infrequently in discrete regions, which is similar to TFs but unlike NAPs, which bind a much larger genomic fraction. However, unlike TFs that typically bind in intergenic regions, HpyA and HstA binding sites are located in both coding and intergenic regions. The genome-wide dinucleotide periodicity known to facilitate histone binding was undetectable in the genomes of both species. Instead, TF-like and histone-like binding sequence preferences were detected for HstA and HpyA, respectively. Taken together, these data suggest that halophilic archaeal histones are unlikely to facilitate genome-wide chromatin formation and that their function defies categorization as a TF, NAP, or histone. IMPORTANCE Most cells in eukaryotic species-from yeast to humans-possess histone proteins that pack and unpack DNA in response to environmental cues. These essential proteins regulate genes necessary for important cellular processes, including development and stress protection. Although the histone fold domain originated in the domain of life Archaea, the function of archaeal histone-like proteins is not well understood relative to those of eukaryotes. We recently discovered that, unlike histones of eukaryotes, histones in hypersaline-adapted archaeal species do not package DNA and can act as transcription factors (TFs) to regulate stress response gene expression. However, the function of histones across species of hypersaline-adapted archaea still remains unclear. Here, we compare hypersaline histone function to a variety of DNA binding proteins across the tree of life, revealing histone-like behavior in some respects and specific transcriptional regulatory function in others.

Two modes of evolution shape bacterial strain diversity in the mammalian gut for thousands of generations

How and at what pace bacteria evolve when colonizing healthy hosts remains unclear. Here, by monitoring evolution for more than six thousand generations in the mouse gut, we show that the successful colonization of an invader Escherichia coli depends on the diversity of the existing microbiota and the presence of a closely related strain. Following colonization, two modes of evolution were observed: one in which diversifying selection leads to long-term coexistence of ecotypes and a second in which directional selection propels selective sweeps. These modes can be quantitatively distinguished by the statistics of mutation trajectories. In our experiments, diversifying selection was marked by the emergence of metabolic mutations, and directional selection by acquisition of prophages, which bring their own benefits and costs. In both modes, we observed parallel evolution, with mutation accumulation rates comparable to those typically observed in vitro on similar time scales. Our results show how rapid ecotype formation and phage domestication can be in the mammalian gut.

Here, the authors show that a colonizing bacterial strain evolves in the gut by either generating ecotypes or continuously fixing beneficial mutations. They associate the first mode to metabolic mutations and the second to domestication of bacteriophages that are incorporated into the bacterial genome.

Discovery of positive and purifying selection in metagenomic time series of hypermutator microbial populations

A general method to infer both positive and purifying selection during the real-time evolution of hypermutator pathogens would be broadly useful. To this end, we introduce a Simple Test to Infer Mode of Selection (STIMS) from metagenomic time series of evolving microbial populations. We test STIMS on metagenomic data generated by simulations of bacterial evolution, and on metagenomic data spanning 62,750 generations of Lenski’s long-term evolution experiment with Escherichia coli (LTEE). This benchmarking shows that STIMS detects positive selection in both nonmutator and hypermutator populations, and purifying selection in hypermutator populations. Using STIMS, we find strong evidence of ongoing positive selection on key regulators of the E. coli gene regulatory network, even in some hypermutator populations. STIMS also detects positive selection on regulatory genes in hypermutator populations of Pseudomonas aeruginosa that adapted to subinhibitory concentrations of colistin–an antibiotic of last resort–for just twenty-six days of laboratory evolution. Our results show that the fine-tuning of gene regulatory networks is a general mechanism for rapid and ongoing adaptation. The simplicity of STIMS, together with its intuitive visual interpretation, make it a useful test for positive and purifying selection in metagenomic data sets that track microbial evolution in real-time.

Organisms often evolve elevated mutation rates as they adapt to new environments. That said, the genomic basis of adaptation in asexual hypermutator populations has been difficult to discern, because the few mutations driving adaptation are often associated with large numbers of non-adaptive “hitchhiker” mutations elsewhere in the genome. To address this research gap, we present a Simple Test to Infer Mode of Selection (STIMS) which aggregates the overall signature of selection seen over hundreds of genes. Using STIMS, we find strong positive selection on key regulators of the E. coli gene regulatory network across nonmutator and hypermutator populations of Lenski’s long-term evolution experiment with Escherichia coli. We also find positive selection on regulatory genes in populations of Pseudomonas aeruginosa that evolved hypermutability as they evolved to resist colistin, an antibiotic of last resort. Our work shows that positive and purifying selection can be resolved in hypermutator populations, given sufficient genomic data, and show that regulatory mutations are a major driver of adaptation in experimental populations of microbes.

Host translesion polymerases are required for viral genome integrity

Human cells encode up to 15 DNA polymerases with specialized functions in chromosomal DNA synthesis and damage repair. In contrast, complex DNA viruses, such as those of the herpesviridae family, encode a single B-family DNA polymerase. This disparity raises the possibility that DNA viruses may rely on host polymerases for synthesis through complex DNA geometries. We tested the importance of error-prone Y-family polymerases involved in translesion synthesis (TLS) to human cytomegalovirus (HCMV) infection. We find most Y-family polymerases involved in the nucleotide insertion and bypass of lesions restrict HCMV genome synthesis and replication. In contrast, other TLS polymerases, such as the polymerase ζ complex, which extends past lesions, was required for optimal genome synthesis and replication. Depletion of either the polζ complex or the suite of insertion polymerases demonstrate that TLS polymerases suppress the frequency of viral genome rearrangements, particularly at GC-rich sites and repeat sequences. Moreover, while distinct from HCMV, replication of the related herpes simplex virus type 1 is impacted by host TLS polymerases, suggesting a broader requirement for host polymerases for DNA virus replication. These findings reveal an unexpected role for host DNA polymerases in ensuring viral genome stability.

Adaptation to novel spatially-structured environments is driven by the capsule and alters virulence-associated traits

The extracellular capsule is a major virulence factor, but its ubiquity in free-living bacteria with large environmental breadths suggests that it shapes adaptation to novel niches. Yet, how it does so, remains unexplored. Here, we evolve three Klebsiella strains and their capsule mutants in parallel. Their comparison reveals different phenotypic and genotypic evolutionary changes that alter virulence-associated traits. Non-capsulated populations accumulate mutations that reduce exopolysaccharide production and increase biofilm formation and yield, whereas most capsulated populations become hypermucoviscous, a signature of hypervirulence. Hence, adaptation to novel environments primarily occurs by fine-tuning expression of the capsular locus. The same evolutionary conditions selecting for mutations in the capsular gene wzc leading to hypermucoviscosity also result in increased susceptibility to antibiotics by mutations in the ramA regulon. This implies that general adaptive processes outside the host can affect capsule evolution and its role in virulence and infection outcomes may be a by-product of such adaptation.

Phenotypic and genotypic evolution in worrisome Klebsiella spp. is influenced by the capsule. Here the authors show that adaptation outside the host can impact virulence-associated traits, including de novo emergence of hypermucoviscosity.

Bacteriophages evolve enhanced persistence to a mucosal surface

The majority of viruses within the gut are obligate bacterial viruses known as bacteriophages (phages). Their bacteriotropism underscores the study of phage ecology in the gut, where they modulate and coevolve with gut bacterial communities. Traditionally, these ecological and evolutionary questions were investigated empirically via in vitro experimental evolution and, more recently, in vivo models were adopted to account for physiologically relevant conditions of the gut. Here, we probed beyond conventional phage-bacteria coevolution to investigate potential tripartite evolutionary interactions between phages, their bacterial hosts, and the mammalian gut mucosa. To capture the role of the mammalian gut, we recapitulated a life-like gut mucosal layer using in vitro lab-on-a-chip devices (to wit, the gut-on-a-chip) and showed that the mucosal environment supports stable phage-bacteria coexistence. Next, we experimentally coevolved lytic phage populations within the gut-on-a-chip devices alongside their bacterial hosts. We found that while phages adapt to the mucosal environment via de novo mutations, genetic recombination was the key evolutionary force in driving mutational fitness. A single mutation in the phage capsid protein Hoc-known to facilitate phage adherence to mucus-caused altered phage binding to fucosylated mucin glycans. We demonstrated that the altered glycan-binding phenotype provided the evolved mutant phage a competitive fitness advantage over its ancestral wild-type phage in the gut-on-a-chip mucosal environment. Collectively, our findings revealed that phages-in addition to their evolutionary relationship with bacteria-are able to evolve in response to a mammalian-derived mucosal environment.

Bistable auto-aggregation phenotype in Lactiplantibacillus plantarum emerges after cultivation in in vitro colonic microbiota

Background

Auto-aggregation is a desired property for probiotic strains because it is suggested to promote colonization of the human intestine, to prevent pathogen infections and to modulate the colonic mucosa. We recently reported the generation of adapted mutants of Lactiplantibacillus plantarum NZ3400, a derivative of the model strain WCFS1, for colonization under adult colonic conditions of PolyFermS continuous intestinal fermentation models. Here we describe and characterize the emerge of an auto-aggregating phenotype in L. plantarum NZ3400 derivatives recovered from the modelled gut microbiota.

Results

L. plantarum isolates were recovered from reactor effluent of four different adult microbiota and from spontaneously formed reactor biofilms. Auto-aggregation was observed in L. plantarum recovered from all microbiota and at higher percentage when recovered from biofilm than from effluent. Further, auto-aggregation percentage increased over time of cultivation in the microbiota. Starvation of the gut microbiota by interrupting the inflow of nutritive medium enhanced auto-aggregation, suggesting a link to nutrient availability. Auto-aggregation was lost under standard cultivation conditions for lactobacilli in MRS medium. However, it was reestablished during growth on sucrose and maltose and in a medium that simulates the abiotic gut environment. Remarkably, none of these conditions resulted in an auto-aggregation phenotype in the wild type strain NZ3400 nor other non-aggregating L. plantarum, indicating that auto-aggregation depends on the strain history. Whole genome sequencing analysis did not reveal any mutation responsible for the auto-aggregation phenotype. Transcriptome analysis showed highly significant upregulation of LP_RS05225 (msa) at 4.1–4.4 log₂-fold-change and LP_RS05230 (marR) at 4.5–5.4 log₂-fold-change in all auto-aggregating strains compared to non-aggregating. These co-expressed genes encode a mannose-specific adhesin protein and transcriptional regulator, respectively. Mapping of the RNA-sequence reads to the promoter region of the msa-marR operon reveled a DNA inversion in this region that is predominant in auto-aggregating but not in non-aggregating strains. This strongly suggests a role of this inversion in the auto-aggregation phenotype.

Conclusions

L. plantarum NZ3400 adapts to the in vitro colonic environment by developing an auto-aggregation phenotype. Similar aggregation phenotypes may promote gut colonization and efficacy of other probiotics and should be further investigated by using validated continuous models of gut fermentation such as PolyFermS.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-021-02331-x.

Coping with darkness: The adaptive response of marine picocyanobacteria to repeated light energy deprivation

The picocyanobacteria Prochlorococcus and Synechococcus are found throughout the ocean's euphotic zone, where the daily light:dark cycle drives their physiology. Periodic deep mixing events can, however, move cells below this region, depriving them of light for extended periods of time. Here, we demonstrate that members of these genera can adapt to tolerate repeated periods of light energy deprivation. Strains kept in the dark for 3 d and then returned to the light initially required 18-26 d to resume growth, but after multiple rounds of dark exposure they began to regrow after only 1-2 d. This dark-tolerant phenotype was stable and heritable; some cultures retained the trait for over 132 generations even when grown in a standard 13:11 light:dark cycle. We found no genetic differences between the dark-tolerant and parental strains of Prochlorococcus NATL2A, indicating that an epigenetic change is likely responsible for the adaptation. To begin to explore this possibility, we asked whether DNA methylation-one potential mechanism mediating epigenetic inheritance in bacteria-occurs in Prochlorococcus. LC-MS/MS analysis showed that while DNA methylations, including 6 mA and 5 mC, are found in some other Prochlorococcus strains, there were no methylations detected in either the parental or dark-tolerant NATL2A strains. These findings suggest that Prochlorococcus utilizes a yet-to-be-determined epigenetic mechanism to adapt to the stress of extended light energy deprivation, and highlights phenotypic heterogeneity as an additional dimension of Prochlorococcus diversity.

O-Specific Antigen-Dependent Surface Hydrophobicity Mediates Aggregate Assembly Type in Pseudomonas aeruginosa

Bacteria live in spatially organized aggregates during chronic infections, where they adapt to the host environment, evade immune responses, and resist therapeutic interventions. Although it is known that environmental factors such as polymers influence bacterial aggregation, it is not clear how bacterial adaptation during chronic infection impacts the formation and spatial organization of aggregates in the presence of polymers. Here, we show that in an in vitro model of cystic fibrosis (CF) containing the polymers extracellular DNA (eDNA) and mucin, O-specific antigen is a major factor determining the formation of two distinct aggregate assembly types of Pseudomonas aeruginosa due to alterations in cell surface hydrophobicity. Our findings suggest that during chronic infection, the interplay between cell surface properties and polymers in the environment may influence the formation and structure of bacterial aggregates, which would shed new light on the fitness costs and benefits of O-antigen production in environments such as CF lungs.

A Comprehensive Map of Mycobacterium tuberculosis Complex Regions of Difference

Mycobacterium tuberculosis complex (MTBC) species are classic examples of genetically monomorphic microorganisms due to their low genetic variability. Whole-genome sequencing made it possible to describe both the main species within the complex and M. tuberculosis lineages and sublineages. This differentiation is based on single nucleotide polymorphisms (SNPs) and large sequence polymorphisms in the so-called regions of difference (RDs). Although a number of studies have been performed to elucidate RD localizations, their distribution among MTBC species, and their role in the bacterial life cycle, there are some inconsistencies and ambiguities in the localization of RDs in different members of the complex. To address this issue, we conducted a thorough search for all possible deletions in the WGS data collection comprising 721 samples representing the full MTBC diversity. Discovered deletions were compared with a list of all previously described RDs. As with the SNP-based analysis, we confirmed the specificities of 79 regions at the species, lineage, or sublineage level, 17 of which are described for the first time. We also present RDscan (https://github.com/dbespiatykh/RDscan), an open-source workflow, which detects deletions from short-read sequencing data and correlates the results with high-specificity RDs, curated in this study. Testing of the workflow on a collection comprising ∼7,000 samples showed a high specificity of the found RDs. This study provides novel details that can contribute to a better understanding of the species differentiation within the MTBC and can help to determine how individual clusters evolve within various MTBC species.

IMPORTANCE Reductive genome evolution is one of the most important and intriguing adaptation strategies of different living organisms to their environment. Mycobacterium offers several notorious examples of either naturally reduced (Mycobacterium leprae) or laboratory-reduced (Mycobacterium bovis BCG) genomes. Mycobacterium tuberculosis complex has its phylogeny unambiguously framed by large sequence polymorphisms that present unidirectional unique event changes. In the present study, we curated all known regions of difference and analyzed both Mycobacterium tuberculosis and animal-adapted MTBC species. For 79 loci, we have shown a relationship with phylogenetic units, which can serve as a marker for diagnosing or studying biological effects. Moreover, intersections were found for some loci, which may indicate the nonrandomness of these processes and the involvement of these regions in the adaptation of bacteria to external conditions.

Experimental evolution in morbidostat reveals converging genomic trajectories on the path to triclosan resistance

Understanding the dynamics and mechanisms of acquired drug resistance across major classes of antibiotics and bacterial pathogens is of critical importance for the optimization of current anti-infective therapies and the development of novel ones. To systematically address this challenge, we developed a workflow combining experimental evolution in a morbidostat continuous culturing device with deep genomic sequencing of population samples collected in time series. This approach was applied to the experimental evolution of six populations of Escherichia coli BW25113 towards acquiring resistance to triclosan (TCS), an antibacterial agent in various consumer products. This study revealed the rapid emergence and expansion (up to 100% in each culture within 4 days) of missense mutations in the fabI gene, encoding enoyl-acyl carrier protein reductase, the known TCS molecular target. A follow-up analysis of isolated clones showed that distinct amino acid substitutions increased the drug IC₉₀ in a 3-16-fold range, reflecting their proximity to the TCS-binding site. In contrast to other antibiotics, efflux-upregulating mutations occurred only rarely and with low abundance. Mutations in several other genes were detected at an earlier stage of evolution. Most notably, three distinct amino acid substitutions were mapped in the C-terminal periplasmic domain of CadC protein, an acid stress-responsive transcriptional regulator. While these mutations do not confer robust TCS resistance, they appear to play a certain, yet unknown, role in adaptation to relatively low drug pressure. Overall, the observed evolutionary trajectories suggest that the FabI enzyme is the sole target of TCS (at least up to the ~50 µm level), and amino acid substitutions in the TCS-binding site represent the main mechanism of robust TCS resistance in E. coli. This model study illustrates the potential utility of the established morbidostat-based approach for uncovering resistance mechanisms and target identification for novel drug candidates with yet unknown mechanisms of action.

Nutrient conditions are primary drivers of bacterial capsule maintenance in Klebsiella

The fitness cost associated with the production of bacterial capsules is considered to be offset by the protection provided by these extracellular structures against biotic aggressions or abiotic stress. However, it is unknown if the capsule contributes to fitness in the absence of these. Here, we explored conditions favouring the maintenance of the capsule in Klebsiella pneumoniae, where the capsule is known to be a major virulence factor. Using short-term experimental evolution on different Klebsiella strains, we showed that small environmental variations have a strong impact on the maintenance of the capsule. Capsule inactivation is frequent in nutrient-rich, but scarce in nutrient-poor media. Competitions between wild-type and capsule mutants in nine different strains confirmed that the capsule is costly in nutrient-rich media. Surprisingly, these results also showed that the presence of a capsule provides a clear fitness advantage in nutrient-poor conditions by increasing both growth rates and population yields. The comparative analyses of the wild-type and capsule mutants reveal complex interactions between the environment, genetic background and serotype even in relation to traits known to be relevant during pathogenesis. In conclusion, our data suggest there are novel roles for bacterial capsules yet to be discovered and further supports the notion that the capsule's role in virulence may be a by-product of its contribution to bacterial adaptation outside the host.

Short and long-read ultra-deep sequencing profiles emerging heterogeneity across five platform Escherichia coli strains

Reprogramming organisms for large-scale bioproduction counters their evolutionary objectives of fast growth and often leads to mutational collapse of the engineered production pathways during cultivation. Yet, the mutational susceptibility of academic and industrial Escherichia coli bioproduction host strains are poorly understood. In this study, we apply 2^nd and 3^rd generation deep sequencing to profile simultaneous modes of genetic heterogeneity that decimate engineered biosynthetic production in five popular E. coli hosts BL21(DE3), TOP10, MG1655, W, and W3110 producing 2,3-butanediol and mevalonic acid. Combining short-read and long-read sequencing, we detect strain and sequence-specific mutational modes including single nucleotide polymorphism, inversion, and mobile element transposition, as well as complex structural variations that disrupt the integrity of the engineered biosynthetic pathway. Our analysis suggests that organism engineers should avoid chassis strains hosting active insertion sequence (IS) subfamilies such as IS1 and IS10 present in popular E. coli TOP10. We also recommend monitoring for increased mutagenicity in the pathway transcription initiation regions and recombinogenic repeats. Together, short and long sequencing reads identified latent low-frequency mutation events such as a short detrimental inversion within a pathway gene, driven by 8-bp short inverted repeats. This demonstrates the power of combining ultra-deep DNA sequencing technologies to profile genetic heterogeneities of engineered constructs and explore the markedly different mutational landscapes of common E. coli host strains. The observed multitude of evolving variants underlines the usefulness of early mutational profiling for new synthetic pathways designed to sustain in organisms over long cultivation scales.

The birth of a bacterial tRNA gene by large-scale, tandem duplication events

Organisms differ in the types and numbers of tRNA genes that they carry. While the evolutionary mechanisms behind tRNA gene set evolution have been investigated theoretically and computationally, direct observations of tRNA gene set evolution remain rare. Here, we report the evolution of a tRNA gene set in laboratory populations of the bacterium Pseudomonas fluorescens SBW25. The growth defect caused by deleting the single-copy tRNA gene, serCGA, is rapidly compensated by large-scale (45–290 kb) duplications in the chromosome. Each duplication encompasses a second, compensatory tRNA gene (serTGA) and is associated with a rise in tRNA-Ser(UGA) in the mature tRNA pool. We postulate that tRNA-Ser(CGA) elimination increases the translational demand for tRNA-Ser(UGA), a pressure relieved by increasing serTGA copy number. This work demonstrates that tRNA gene sets can evolve through duplication of existing tRNA genes, a phenomenon that may contribute to the presence of multiple, identical tRNA gene copies within genomes.

Network Rewiring: Physiological Consequences of Reciprocally Exchanging the Physical Locations and Growth-Phase-Dependent Expression Patterns of the Salmonella fis and dps Genes

We assessed the impact on Salmonella physiology of reciprocally translocating the genes encoding the Fis and Dps nucleoid-associated proteins (NAPs) and of inverting their growth-phase production patterns such that Fis was produced in stationary phase (like Dps) and Dps was produced in exponential phase (like Fis). Changes to peak binding of Fis were detected by ChIP-seq on the chromosome, as were widespread impacts on the transcriptome, especially when Fis production mimicked Dps production. Virulence gene expression and the expression of a virulence phenotype were altered. Overall, these radical changes to NAP gene expression were well tolerated, revealing the robust and well-buffered nature of global gene regulation networks in the bacterium.

The Fis nucleoid-associated protein controls the expression of a large and diverse regulon of genes in Gram-negative bacteria. Fis production is normally maximal in bacteria during the early exponential phase of batch culture growth, becoming almost undetectable by the onset of stationary phase. We tested the effect on the Fis regulatory network in Salmonella of moving the complete fis gene from its usual location near the origin of chromosomal replication to the position normally occupied by the dps gene in the right macrodomain of the chromosome, and vice versa, creating the gene exchange (GX) strain. In a parallel experiment, we tested the effect of rewiring the Fis regulatory network by placing the fis open reading frame under the control of the stationary-phase-activated dps promoter at the dps genetic location within the right macrodomain, and vice versa, creating the open reading frame exchange (OX) strain. Chromatin immunoprecipitation sequencing (ChIP-seq) was used to measure global Fis protein binding levels and to determine gene expression patterns. Strain GX showed few changes compared with the wild type, although we did detect increased Fis binding at Ter, accompanied by reduced binding at Ori. Strain OX displayed a more pronounced version of this distorted Fis protein-binding pattern together with numerous alterations in the expression of genes in the Fis regulon. OX, but not GX, had a reduced ability to infect cultured mammalian cells. These findings illustrate the inherent robustness of the Fis regulatory network with respect to the effects of rewiring based on gene repositioning alone and emphasize the importance of fis expression signals in phenotypic determination.

A Comparison of Whole Genome Sequencing of SARS-CoV-2 Using Amplicon-Based Sequencing, Random Hexamers, and Bait Capture

Genome sequencing of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is increasingly important to monitor the transmission and adaptive evolution of the virus. The accessibility of high-throughput methods and polymerase chain reaction (PCR) has facilitated a growing ecosystem of protocols. Two differing protocols are tiling multiplex PCR and bait capture enrichment. Each method has advantages and disadvantages but a direct comparison with different viral RNA concentrations has not been performed to assess the performance of these approaches. Here we compare Liverpool amplification, ARTIC amplification, and bait capture using clinical diagnostics samples. All libraries were sequenced using an Illumina MiniSeq with data analyzed using a standardized bioinformatics workflow (SARS-CoV-2 Illumina GeNome Assembly Line; SIGNAL). One sample showed poor SARS-CoV-2 genome coverage and consensus, reflective of low viral RNA concentration. In contrast, the second sample had a higher viral RNA concentration, which yielded good genome coverage and consensus. ARTIC amplification showed the highest depth of coverage results for both samples, suggesting this protocol is effective for low concentrations. Liverpool amplification provided a more even read coverage of the SARS-CoV-2 genome, but at a lower depth of coverage. Bait capture enrichment of SARS-CoV-2 cDNA provided results on par with amplification. While only two clinical samples were examined in this comparative analysis, both the Liverpool and ARTIC amplification methods showed differing efficacy for high and low concentration samples. In addition, amplification-free bait capture enriched sequencing of cDNA is a viable method for generating a SARS-CoV-2 genome sequence and for identification of amplification artifacts.

Mutation accumulation and horizontal gene transfer in Escherichia coli colonizing the gut of old mice

The ecology and environment of the microbes that inhabit the mammalian intestine undergoes several changes as the host ages. Here, we ask if the selection pressure experienced by a new strain colonizing the aging gut differs from that in the gut of young adults. Using experimental evolution in mice after a short antibiotic treatment, as a model for a common clinical situation, we show that a new colonizing E. coli strain rapidly adapts to the aging gut via both mutation accumulation and bacteriophage-mediated horizontal gene transfer (HGT). The pattern of evolution of E. coli in aging mice is characterized by a larger number of transposable element insertions and intergenic mutations compared to that in young mice, which is consistent with the gut of aging hosts harboring a stressful and iron limiting environment.

Acid Experimental Evolution of the Haloarchaeon Halobacterium sp. NRC-1 Selects Mutations Affecting Arginine Transport and Catabolism

Halobacterium sp. NRC-1 (NRC-1) is an extremely halophilic archaeon that is adapted to multiple stressors such as UV, ionizing radiation and arsenic exposure; it is considered a model organism for the feasibility of microbial life in iron-rich brine on Mars. We conducted experimental evolution of NRC-1 under acid and iron stress. NRC-1 was serially cultured in CM+ medium modified by four conditions: optimal pH (pH 7.5), acid stress (pH 6.3), iron amendment (600 μM ferrous sulfate, pH 7.5), and acid plus iron (pH 6.3, with 600 μM ferrous sulfate). For each condition, four independent lineages of evolving populations were propagated. After 500 generations, 16 clones were isolated for phenotypic characterization and genomic sequencing. Genome sequences of all 16 clones revealed 378 mutations, of which 90% were haloarchaeal insertion sequences (ISH) and ISH-mediated large deletions. This proportion of ISH events in NRC-1 was five-fold greater than that reported for comparable evolution of Escherichia coli. One acid-evolved clone had increased fitness compared to the ancestral strain when cultured at low pH. Seven of eight acid-evolved clones had a mutation within or upstream of arcD, which encodes an arginine-ornithine antiporter; no non-acid adapted strains had arcD mutations. Mutations also affected the arcR regulator of arginine catabolism, which protects bacteria from acid stress by release of ammonia. Two acid-adapted strains shared a common mutation in bop, which encodes bacterio-opsin, apoprotein for the bacteriorhodopsin light-driven proton pump. Thus, in the haloarchaeon NRC-1, as in bacteria, pH adaptation was associated with genes involved in arginine catabolism and proton transport. Our study is among the first to report experimental evolution with multiple resequenced genomes of an archaeon. Haloarchaea are polyextremophiles capable of growth under environmental conditions such as concentrated NaCl and desiccation, but little is known about pH stress. Interesting parallels appear between the molecular basis of pH adaptation in NRC-1 and in bacteria, particularly the acid-responsive arginine-ornithine system found in oral streptococci.

Evolution of satellite plasmids can prolong the maintenance of newly acquired accessory genes in bacteria

Transmissible plasmids spread genes encoding antibiotic resistance and other traits to new bacterial species. Here we report that laboratory populations of Escherichia coli with a newly acquired IncQ plasmid often evolve 'satellite plasmids' with deletions of accessory genes and genes required for plasmid replication. Satellite plasmids are molecular parasites: their presence reduces the copy number of the full-length plasmid on which they rely for their continued replication. Cells with satellite plasmids gain an immediate fitness advantage from reducing burdensome expression of accessory genes. Yet, they maintain copies of these genes and the complete plasmid, which potentially enables them to benefit from and transmit the traits they encode in the future. Evolution of satellite plasmids is transient. Cells that entirely lose accessory gene function or plasmid mobility dominate in the long run. Satellite plasmids also evolve in Snodgrassella alvi colonizing the honey bee gut, suggesting that this mechanism may broadly contribute to the importance of IncQ plasmids as agents of bacterial gene transfer in nature.

Adaptation to mutational inactivation of an essential gene converges to an accessible suboptimal fitness peak

The mechanisms of adaptation to inactivation of essential genes remain unknown. Here we inactivate E. coli dihydrofolate reductase (DHFR) by introducing D27G,N,F chromosomal mutations in a key catalytic residue with subsequent adaptation by an automated serial transfer protocol. The partial reversal G27- > C occurred in three evolutionary trajectories. Conversely, in one trajectory for D27G and in all trajectories for D27F,N strains adapted to grow at very low metabolic supplement (folAmix) concentrations but did not escape entirely from supplement auxotrophy. Major global shifts in metabolome and proteome occurred upon DHFR inactivation, which were partially reversed in adapted strains. Loss-of-function mutations in two genes, thyA and deoB, ensured adaptation to low folAmix by rerouting the 2-Deoxy-D-ribose-phosphate metabolism from glycolysis towards synthesis of dTMP. Multiple evolutionary pathways of adaptation converged to a suboptimal solution due to the high accessibility to loss-of-function mutations that block the path to the highest, yet least accessible, fitness peak.

Predicting how viruses and bacteria evolve remains a challenge. The ability to anticipate when and how bacteria might develop drug resistance would make treating life-threatening diseases easier and could potentially help prevent drug resistance altogether. Studying bacterial evolution under different conditions and cataloguing all possible DNA mutations that allow these bacteria to survive are crucial steps in predicting the appearance of drug resistance.

Studies have revealed that bacteria can adapt to sources of stress, such as antibiotics, in different ways – each involving mutations in distinct genes. However, not all the mutations provide the same benefits to the organism, and the factors that influence how bacteria will adapt are unclear.

Now, Rodrigues and Shakhnovich have used a new approach to push the adaptation ability of the bacterium Escherichia coli to the limit. First, they genetically engineered different E. coli strains by introducing distinct mutations to an enzyme the bacterium needs to make DNA. These mutations make the resulting strains dependent on external molecules to synthesize new DNA. Next, the cells were grown in conditions where the supply of these DNA precursors was progressively decreased, forcing them to adapt. The obvious way for bacteria to adapt to these conditions would be to ‘revert’ the mutations that Rodrigues and Shakhnovich introduced in the first place. By using this approach, Rodrigues and Shakhnovich were able to test how often the obvious evolutionary solution happens compared with presumably less-preferred alternative routes.

In rare cases, a specific mutation did restore the activity of the enzyme that enabled DNA synthesis. However, in most cases the bacteria found a different evolutionary solution whereby they all adapt to the decrease in external DNA precursors in the same way, but not by reverting the original mutation. Instead, further mutations disrupt the activity of two metabolic genes, allowing the cells to use the external DNA precursors more efficiently, and keep building DNA. These cells are therefore able to survive even when the levels of the external DNA components are very low, but they will die in the complete absence of these precursor molecules.

This evolutionary solution leads to a non-optimal effect: mutations that restore the activity of the original enzyme, which are the best solution when the two metabolic genes are intact, are no longer as effective. This finding represents a clear example of interactions between genes determining evolutionary outcomes. Rodrigues and Shakhnovich showed that, since it is more likely for a random mutation to disrupt a gene than to revert a previous mutation, adaptations that are less-than-optimal but still work might predominate simply because they happen faster.

Understanding why certain evolutionary adaptations prevail is an important step in predicting evolution and may lead to breakthroughs in many areas. For example, if scientists can identify mutations likely to make bacteria resistant to drugs, it may be possible to act proactively against the bacterial strains that carry those mutations. Eventually, if the factors that lead to specific adaptations are known, it may be possible to exploit this knowledge to create weaknesses in the bacteria’s own defences.

Experimental Evolution of Escherichia coli K-12 in the Presence of Proton Motive Force (PMF) Uncoupler Carbonyl Cyanide m-Chlorophenylhydrazone Selects for Mutations Affecting PMF-Driven Drug Efflux Pumps

Experimental evolution of Escherichia coli K-12 with benzoate, a partial uncoupler of the proton motive force (PMF), selects for mutations that decrease antibiotic resistance. We conducted experimental evolution in the presence of carbonyl cyanide m-chlorophenylhydrazone (CCCP), a strong uncoupler. Cultures were serially diluted daily 1:100 in LBK medium containing 20 to 150 µM CCCP buffered at pH 6.5 or at pH 8.0. After 1,000 generations, the populations tolerated up to 150 µM CCCP. Sequenced isolates had mutations in mprA (emrR), which downregulates the EmrAB-TolC pump that exports CCCP. A mprA::kanR deletion conferred growth at 60 μM CCCP, though not at the higher levels resisted by evolved strains (150 µM). Some mprA mutant strains also had point mutations affecting emrA, but deletion of emrA abolished the CCCP resistance. Thus, CCCP-evolved isolates contained additional adaptations. One isolate lacked emrA or mprA mutations but had mutations in cecR (ybiH), whose product upregulates drug pumps YbhG and YbhFSR, and in gadE, which upregulates the multidrug pump MdtEF. A cecR::kanR deletion conferred partial resistance to CCCP. Other multidrug efflux genes that had mutations included ybhR and acrAB The acrB isolate was sensitive to the AcrAB substrates chloramphenicol and tetracycline. Other mutant genes in CCCP-evolved strains include rng (RNase G) and cyaA (adenylate cyclase). Overall, experimental evolution revealed a CCCP-dependent fitness advantage for mutations increasing CCCP efflux via EmrA and for mutations that may deactivate proton-driven pumps for drugs not present (cecR, gadE, acrAB, and ybhR). These results are consistent with our previous report of drug sensitivity associated with evolved benzoate tolerance.IMPORTANCE The genetic responses of bacteria to depletion of proton motive force (PMF), and their effects on drug resistance, are poorly understood. PMF drives export of many antibiotics, but the energy cost may decrease fitness when antibiotics are absent. Our evolution experiment reveals genetic mechanisms of adaptation to the PMF uncoupler CCCP, including selection for increased CCCP efflux but also against the expression of PMF-driven pumps for drugs not present. The results have implications for our understanding of the gut microbiome, which experiences high levels of organic acids that decrease PMF.

Continuous culture of Escherichia coli, under selective pressure by a novel antimicrobial complex, does not result in development of resistance

We attempted to generate de novo resistance to a newly described biocidal complex, ITC (iodo-thiocyanate complex), and to levofloxacin (LVX) in Escherichia coli ATCC 25922, by means of selective chemostat culture. We measured resistance by determining the minimum inhibitory concentrations (MICs) for these agents. E. coli underwent 20-day parallel adaptive evolution routes under no antimicrobial selection, and gradually increasing ITC and LVX selection pressure. Long-term exposure of E. coli to ITC did not induce resistance to ITC, or cross-resistance to LVX. No distinct mutational pattern was evidenced from whole-genome sequence (WGS)-based comparisons of ITC-challenged and unchallenged bacterial populations. Moreover, the exposed E. coli population could not survive a 2 × MIC challenge of ITC. By contrast, resistance to LVX was rapidly induced (on day 1 the MIC had increased 16-fold), selected for (by day 14 the MIC had increased 64-fold) and enriched with a highly characteristic genome mutational pattern. WGS of this evolving population revealed that the majority of mutations appeared in the genes of LVX target proteins (GyrA, ParC, ParE) and drug influx (OmpF). This study suggests that the usage of ITC may not trigger the emergence of facile resistance or cross-resistance, in contrast to common antibiotics.

Genomic correlates of extraintestinal infection are linked with changes in cell morphology in Campylobacter jejuni

Campylobacter jejuni is the most common cause of bacterial diarrheal disease in the world. Clinical outcomes of infection can range from asymptomatic infection to life-threatening extraintestinal infections. This variability in outcomes for infected patients has raised questions as to whether genetic differences between C. jejuni isolates contribute to their likelihood of causing severe disease. In this study, we compare the genomes of ten C. jejuni isolates that were implicated in extraintestinal infections with reference gastrointestinal isolates, in order to identify unusual patterns of sequence variation associated with infection outcome. We identified a collection of genes that display a higher burden of uncommon mutations in invasive isolates compared with gastrointestinal close relatives, including some that have been previously linked to virulence and invasiveness in C. jejuni. Among the top genes identified were mreB and pgp1, which are both involved in determining cell shape. Electron microscopy confirmed morphological differences in isolates carrying unusual sequence variants of these genes, indicating a possible relationship between extraintestinal infection and changes in cell morphology.

Experimental Evolution of Escherichia coli K-12 at High pH and with RpoS Induction

Experimental evolution of Escherichia coli K-12 W3110 by serial dilutions for 2,200 generations at high pH extended the range of sustained growth from pH 9.0 to pH 9.3. pH 9.3-adapted isolates showed mutations in DNA-binding regulators and envelope proteins. One population showed an IS1 knockout of phoB (encoding the positive regulator of the phosphate regulon). A phoB::kanR knockout increased growth at high pH. phoB mutants are known to increase production of fermentation acids, which could enhance fitness at high pH. Mutations in pcnB [poly(A) polymerase] also increased growth at high pH. Three out of four populations showed deletions of torI, an inhibitor of TorR, which activates expression of torCAD (trimethylamine N-oxide respiration) at high pH. All populations showed point mutations affecting the stationary-phase sigma factor RpoS, either in the coding gene or in genes for regulators of RpoS expression. RpoS is required for survival at extremely high pH. In our microplate assay, rpoS deletion slightly decreased growth at pH 9.1. RpoS protein accumulated faster at pH 9 than at pH 7. The RpoS accumulation at high pH required the presence of one or more antiadaptors that block degradation (IraM, IraD, and IraP). Other genes with mutations after high-pH evolution encode regulators, such as those encoded by yobG (mgrB) (PhoPQ regulator), rpoN (nitrogen starvation sigma factor), malI, and purR, as well as envelope proteins, such as those encoded by ompT and yahO Overall, E. coli evolution at high pH selects for mutations in key transcriptional regulators, including phoB and the stationary-phase sigma factor RpoS.IMPORTANCEEscherichia coli in its native habitat encounters high-pH stress such as that of pancreatic secretions. Experimental evolution over 2,000 generations showed selection for mutations in regulatory factors, such as deletion of the phosphate regulator PhoB and mutations that alter the function of the global stress regulator RpoS. RpoS is induced at high pH via multiple mechanisms.

Analysis of bacterial genomes from an evolution experiment with horizontal gene transfer shows that recombination can sometimes overwhelm selection

Few experimental studies have examined the role that sexual recombination plays in bacterial evolution, including the effects of horizontal gene transfer on genome structure. To address this limitation, we analyzed genomes from an experiment in which Escherichia coli K-12 Hfr (high frequency recombination) donors were periodically introduced into 12 evolving populations of E. coli B and allowed to conjugate repeatedly over the course of 1000 generations. Previous analyses of the evolved strains from this experiment showed that recombination did not accelerate adaptation, despite increasing genetic variation relative to asexual controls. However, the resolution in that previous work was limited to only a few genetic markers. We sought to clarify and understand these puzzling results by sequencing complete genomes from each population. The effects of recombination were highly variable: one lineage was mostly derived from the donors, while another acquired almost no donor DNA. In most lineages, some regions showed repeated introgression and others almost none. Regions with high introgression tended to be near the donors' origin of transfer sites. To determine whether introgressed alleles imposed a genetic load, we extended the experiment for 200 generations without recombination and sequenced whole-population samples. Beneficial alleles in the recipient populations were occasionally driven extinct by maladaptive donor-derived alleles. On balance, our analyses indicate that the plasmid-mediated recombination was sufficiently frequent to drive donor alleles to fixation without providing much, if any, selective advantage.

History of antibiotic adaptation influences microbial evolutionary dynamics during subsequent treatment

Antibiotic regimens often include the sequential changing of drugs to limit the development and evolution of resistance of bacterial pathogens. It remains unclear how history of adaptation to one antibiotic can influence the resistance profiles when bacteria subsequently adapt to a different antibiotic. Here, we experimentally evolved Pseudomonas aeruginosa to six 2-drug sequences. We observed drug order–specific effects, whereby adaptation to the first drug can limit the rate of subsequent adaptation to the second drug, adaptation to the second drug can restore susceptibility to the first drug, or final resistance levels depend on the order of the 2-drug sequence. These findings demonstrate how resistance not only depends on the current drug regimen but also the history of past regimens. These order-specific effects may allow for rational forecasting of the evolutionary dynamics of bacteria given knowledge of past adaptations and provide support for the need to consider the history of past drug exposure when designing strategies to mitigate resistance and combat bacterial infections.

Bacteria readily adapt to their environments and can develop ways to survive and grow in the presence of antibiotics. While many studies have investigated how bacteria evolve to become resistant to single drugs, it is unclear how adaptation to other drugs and environments in the past affect the way bacteria adapt to new drugs and environments. In this study, we allowed bacteria in a laboratory setting to adapt to three different antibiotics. We first exposed wild-type susceptible bacteria to high concentrations of the three antibiotics individually and then exposed these populations to each of the other drugs. By tracking the levels of resistance to all three drugs in all of the treatments, we identified cases in which past adaptation to one treatment influenced subsequent evolutionary dynamics with regard to both phenotypes (levels of resistance) and genotypes (genes that became mutated). Additionally, by allowing bacterial isolates originating from human patients to adapt to the three drugs, we recapitulated a subset of the adaptation history-dependent evolutionary dynamics. Overall, this study sheds light on how adaptation history in bacteria can potentiate or constrain the rise of multidrug resistance, depending on the particular order of drugs used in a sequential therapy.

The role of recombination in evolutionary adaptation of Escherichia coli to a novel nutrient

The benefits and detriments of recombination for adaptive evolution have been studied both theoretically and experimentally, with conflicting predictions and observations. Most pertinent experiments examine recombination's effects in an unchanging environment and do not study its genomewide effects. Here, we evolved six replicate populations of either highly recombining R⁺ or lowly recombining R^- E. coli strains in a changing environment, by introducing the novel nutrients L-arabinose or indole into the environment. The experiment's ancestral strains are not viable on these nutrients, but 130 generations of adaptive evolution were sufficient to render them viable. Recombination conferred a more pronounced advantage to populations adapting to indole. To study the genomic changes associated with this advantage, we sequenced the genomes of 384 clones isolated from selected replicates at the end of the experiment. These genomes harbour complex changes that range from point mutations to large-scale DNA amplifications. Among several candidate adaptive mutations, those in the tryptophanase regulator tnaC stand out, because the tna operon in which it resides has a known role in indole metabolism. One of the highly recombining populations also shows a significant excess of large-scale segmental DNA amplifications that include the tna operon. This lineage also shows a unique and potentially adaptive combination of point mutations and DNA amplifications that may have originated independently from one another, to be joined later by recombination. Our data illustrate that the advantages of recombination for adaptive evolution strongly depend on the environment and that they can be associated with complex genomic changes.

Acid Evolution of Escherichia coli K-12 Eliminates Amino Acid Decarboxylases and Reregulates Catabolism

Acid-adapted strains of Escherichia coli K-12 W3110 were obtained by serial culture in medium buffered at pH 4.6 (M. M. Harden, A. He, K. Creamer, M. W. Clark, I. Hamdallah, K. A. Martinez, R. L. Kresslein, S. P. Bush, and J. L. Slonczewski, Appl Environ Microbiol 81:1932-1941, 2015, https://doi.org/10.1128/AEM.03494-14). Revised genomic analysis of these strains revealed insertion sequence (IS)-driven insertions and deletions that knocked out regulators CadC (acid induction of lysine decarboxylase), GadX (acid induction of glutamate decarboxylase), and FNR (anaerobic regulator). Each acid-evolved strain showed loss of one or more amino acid decarboxylase systems, which normally help neutralize external acid (pH 5 to 6) and increase survival in extreme acid (pH 2). Strains from populations B11, H9, and F11 had an IS5 insertion or IS-mediated deletion in cadC, while population B11 had a point mutation affecting the arginine activator adiY The cadC and adiY mutants failed to neutralize acid in the presence of exogenous lysine or arginine. In strain B11-1, reversion of an rpoC (RNA polymerase) mutation partly restored arginine-dependent neutralization. All eight strains showed deletion or downregulation of the Gad acid fitness island. Strains with the Gad deletion lost the ability to produce GABA (gamma-aminobutyric acid) and failed to survive extreme acid. Transcriptome sequencing (RNA-seq) of strain B11-1 showed upregulated genes for catabolism of diverse substrates but downregulated acid stress genes (the biofilm regulator ariR, yhiM, and Gad). Other strains showed downregulation of H₂ consumption mediated by hydrogenases (hya and hyb) which release acid. Strains F9-2 and F9-3 had a deletion of fnr and showed downregulation of FNR-dependent genes (dmsABC, frdABCD, hybABO, nikABCDE, and nrfAC). Overall, strains that had evolved in buffered acid showed loss or downregulation of systems that neutralize unbuffered acid and showed altered regulation of catabolism.IMPORTANCE Experimental evolution of an enteric bacterium under a narrow buffered range of acid pH leads to loss of genes that enhance fitness above or below the buffered pH range, including loss of enzymes that may raise external pH in the absence of buffer. Prominent modes of evolutionary change involve IS-mediated insertions and deletions that knock out key regulators. Over generations of acid stress, catabolism undergoes reregulation in ways that differ for each evolving strain.

Benzoate- and Salicylate-Tolerant Strains of Escherichia coli K-12 Lose Antibiotic Resistance during Laboratory Evolution

Escherichia coli K-12 W3110 grows in the presence of membrane-permeant organic acids that can depress cytoplasmic pH and accumulate in the cytoplasm. We conducted experimental evolution by daily diluting cultures in increasing concentrations of benzoic acid (up to 20 mM) buffered at external pH 6.5, a pH at which permeant acids concentrate in the cytoplasm. By 2,000 generations, clones isolated from evolving populations showed increasing tolerance to benzoate but were sensitive to chloramphenicol and tetracycline. Sixteen clones grew to stationary phase in 20 mM benzoate, whereas the ancestral strain W3110 peaked and declined. Similar growth occurred in 10 mM salicylate. Benzoate-evolved strains grew like W3110 in the absence of benzoate, in media buffered at pH 4.8, pH 7.0, or pH 9.0, or in 20 mM acetate or sorbate at pH 6.5. Genomes of 16 strains revealed over 100 mutations, including single-nucleotide polymorphisms (SNPs), large deletions, and insertion knockouts. Most strains acquired deletions in the benzoate-induced multiple antibiotic resistance (Mar) regulon or in associated regulators such as rob and cpxA, as well as the multidrug resistance (MDR) efflux pumps emrA, emrY, and mdtA Strains also lost or downregulated the Gad acid fitness regulon. In 5 mM benzoate or in 2 mM salicylate (2-hydroxybenzoate), most strains showed increased sensitivity to the antibiotics chloramphenicol and tetracycline; some strains were more sensitive than a marA knockout strain. Thus, our benzoate-evolved strains may reveal additional unknown drug resistance components. Benzoate or salicylate selection pressure may cause general loss of MDR genes and regulators.

IMPORTANCE

Benzoate is a common food preservative, and salicylate is the primary active metabolite of aspirin. In the gut microbiome, genetic adaptation to salicylate may involve loss or downregulation of inducible multidrug resistance systems. This discovery implies that aspirin therapy may modulate the human gut microbiome to favor salicylate tolerance at the expense of drug resistance. Similar aspirin-associated loss of drug resistance might occur in bacterial pathogens found in arterial plaques.

Protein Homeostasis Imposes a Barrier on Functional Integration of Horizontally Transferred Genes in Bacteria

Horizontal gene transfer (HGT) plays a central role in bacterial evolution, yet the molecular and cellular constraints on functional integration of the foreign genes are poorly understood. Here we performed inter-species replacement of the chromosomal folA gene, encoding an essential metabolic enzyme dihydrofolate reductase (DHFR), with orthologs from 35 other mesophilic bacteria. The orthologous inter-species replacements caused a marked drop (in the range 10–90%) in bacterial growth rate despite the fact that most orthologous DHFRs are as stable as E.coli DHFR at 37°C and are more catalytically active than E. coli DHFR. Although phylogenetic distance between E. coli and orthologous DHFRs as well as their individual molecular properties correlate poorly with growth rates, the product of the intracellular DHFR abundance and catalytic activity (k_cat/K_M), correlates strongly with growth rates, indicating that the drop in DHFR abundance constitutes the major fitness barrier to HGT. Serial propagation of the orthologous strains for ~600 generations dramatically improved growth rates by largely alleviating the fitness barriers. Whole genome sequencing and global proteome quantification revealed that the evolved strains with the largest fitness improvements have accumulated mutations that inactivated the ATP-dependent Lon protease, causing an increase in the intracellular DHFR abundance. In one case DHFR abundance increased further due to mutations accumulated in folA promoter, but only after the lon inactivating mutations were fixed in the population. Thus, by apparently distinguishing between self and non-self proteins, protein homeostasis imposes an immediate and global barrier to the functional integration of foreign genes by decreasing the intracellular abundance of their products. Once this barrier is alleviated, more fine-tuned evolution occurs to adjust the function/expression of the transferred proteins to the constraints imposed by the intracellular environment of the host organism.

Horizontal gene transfer (HGT) is central to bacterial evolution. The outcome of an HGT event (fixation in a population, elimination, or separation as a subdominant clone) depends not only on the availability of a new gene but crucially on the fitness cost or benefit of the genomic incorporation of the foreign gene and its expression in recipient bacteria. Here we studied the fitness landscape for inter-species chromosomal replacement of an essential protein, dihydrofolate reductase (DHFR) encoded by the folA gene, by its orthologs from other mesophilic bacteria. We purified and biochemically characterized 33 out of 35 orthologous DHFRs and found that most of them are stable and more catalytically active than E. coli DHFR. However, the inter-species replacement of DHFR caused significant fitness loss for most transgenic strains due to low abundance of orthologous DHFRs in E. coli cytoplasm. Laboratory evolution resulted in an increase in orthologous DHFR abundance leading to a dramatic fitness improvement. Genomic and proteomic analyses of “naive” and evolved strains suggest a new function of protein homeostasis to discriminate between “self” and “non-self” proteins, thus creating fitness barriers to HGT.

Genomic investigations of evolutionary dynamics and epistasis in microbial evolution experiments

Microbial evolution experiments enable us to watch adaptation in real time, and to quantify the repeatability and predictability of evolution by comparing identical replicate populations. Further, we can resurrect ancestral types to examine changes over evolutionary time. Until recently, experimental evolution has been limited to measuring phenotypic changes, or to tracking a few genetic markers over time. However, recent advances in sequencing technology now make it possible to extensively sequence clones or whole-population samples from microbial evolution experiments. Here, we review recent work exploiting these techniques to understand the genomic basis of evolutionary change in experimental systems. We first focus on studies that analyze the dynamics of genome evolution in microbial systems. We then survey work that uses observations of sequence evolution to infer aspects of the underlying fitness landscape, concentrating on the epistatic interactions between mutations and the constraints these interactions impose on adaptation.

Overcoming Challenges in Engineering the Genetic Code

Withstanding 3.5 billion years of genetic drift, the canonical genetic code remains such a fundamental foundation for the complexity of life that it is highly conserved across all three phylogenetic domains. Genome engineering technologies are now making it possible to rationally change the genetic code, offering resistance to viruses, genetic isolation from horizontal gene transfer, and prevention of environmental escape by genetically modified organisms. We discuss the biochemical, genetic, and technological challenges that must be overcome in order to engineer the genetic code.

Adaptation, Clonal Interference, and Frequency-Dependent Interactions in a Long-Term Evolution Experiment with Escherichia coli

Twelve replicate populations of Escherichia coli have been evolving in the laboratory for >25 years and 60,000 generations. We analyzed bacteria from whole-population samples frozen every 500 generations through 20,000 generations for one well-studied population, called Ara-1. By tracking 42 known mutations in these samples, we reconstructed the history of this population's genotypic evolution over this period. The evolutionary dynamics of Ara-1 show strong evidence of selective sweeps as well as clonal interference between competing lineages bearing different beneficial mutations. In some cases, sets of several mutations approached fixation simultaneously, often conveying no information about their order of origination; we present several possible explanations for the existence of these mutational cohorts. Against a backdrop of rapid selective sweeps both earlier and later, two genetically diverged clades coexisted for >6000 generations before one went extinct. In that time, many additional mutations arose in the clade that eventually prevailed. We show that the clades evolved a frequency-dependent interaction, which prevented the immediate competitive exclusion of either clade, but which collapsed as beneficial mutations accumulated in the clade that prevailed. Clonal interference and frequency dependence can occur even in the simplest microbial populations. Furthermore, frequency dependence may generate dynamics that extend the period of coexistence that would otherwise be sustained by clonal interference alone.