Rate and molecular spectrum of spontaneous mutations in the bacterium Escherichia coli as determined by whole-genome sequencing

Model of neural development by differentiating human induced pluripotent stem cells into neural progenitor cells to study the neurodevelopmental toxicity of lead

Lead (Pb) exposure causes immeasurable damage to multiple human systems, particularly the central nervous system (CNS). In this study, human induced pluripotent stem cells (hiPSCs) were differentiated into neural progenitor cells (NPCs) to investigate the neurotoxic effects of Pb. The hiPSCs were treated with 0, 0.5, 1.0, 2.5, 5.0 and 10.0 μmol/L Pb for 7 days, whereas embryoid bodies (EBs) and NPCs were treated with 0, 0.1, 0.5, and 1.0 μmol/L Pb for 7 days. Pb exposure disrupted the cell cycle and caused apoptosis in hiPSCs, EBs, and NPCs. Besides, Pb inhibited the differentiation of NPCs and EBs. Whole exome sequencing revealed 2509, 2413, and 1984 single nucleotide variants (SNVs) caused by Pb in hiPSCs, EBs, and NPCs, respectively. The common mutation sites in the exon region were mostly nonsynonymous mutations. We identified 18, 19, and 18 common deleterious mutations in hiPSCs, EBs, and NPCs, respectively. Additionally, Online Mendelian Inheritance in Man database analysis revealed 30, 20, and 13 genes related to CNS disorders in hiPSCs, EBs, and NPCs, respectively. Our findings suggest that this in vitro model may supplement animal models and be applied to the study of neurodevelopmental toxicity in the future.

Phenotype switching of the mutation rate facilitates adaptive evolution

The mutation rate plays an important role in adaptive evolution. It can be modified by mutator and anti-mutator alleles. Recent empirical evidence hints that the mutation rate may vary among genetically identical individuals: evidence from bacteria suggests that the mutation rate can be affected by expression noise of a DNA repair protein and potentially also by translation errors in various proteins. Importantly, this non-genetic variation may be heritable via a transgenerational epigenetic mode of inheritance, giving rise to a mutator phenotype that is independent from mutator alleles. Here, we investigate mathematically how the rate of adaptive evolution is affected by the rate of mutation rate phenotype switching. We model an asexual population with two mutation rate phenotypes, non-mutator and mutator. An offspring may switch from its parental phenotype to the other phenotype. We find that switching rates that correspond to so-far empirically described non-genetic systems of inheritance of the mutation rate lead to higher rates of adaptation on both artificial and natural fitness landscapes. These switching rates can maintain within the same individuals both a mutator phenotype and intermediary mutations, a combination that facilitates adaptation. Moreover, non-genetic inheritance increases the proportion of mutators in the population, which in turn increases the probability of hitchhiking of the mutator phenotype with adaptive mutations. This in turns facilitates the acquisition of additional adaptive mutations. Our results rationalize recently observed noise in the expression of proteins that affect the mutation rate and suggest that non-genetic inheritance of this phenotype may facilitate evolutionary adaptive processes.

Retron-mediated multiplex genome editing and continuous evolution in Escherichia coli

While there are several genome editing techniques available, few are suitable for dynamic and simultaneous mutagenesis of arbitrary targeted sequences in prokaryotes. Here, to address these limitations, we present a versatile and multiplex retron-mediated genome editing system (REGES). First, through systematic optimization of REGES, we achieve efficiency of ∼100%, 85 ± 3%, 69 ± 14% and 25 ± 14% for single-, double-, triple- and quadruple-locus genome editing, respectively. In addition, we employ REGES to generate pooled and barcoded variant libraries with degenerate RBS sequences to fine-tune the expression level of endogenous and exogenous genes, such as transcriptional factors to improve ethanol tolerance and biotin biosynthesis. Finally, we demonstrate REGES-mediated continuous in vivo protein evolution, by combining retron, polymerase-mediated base editing and error-prone transcription. By these case studies, we demonstrate REGES as a powerful multiplex genome editing and continuous evolution tool with broad applications in synthetic biology and metabolic engineering.

Selective Pressure by Rifampicin Modulates Mutation Rates and Evolutionary Trajectories of Mycobacterial Genomes

Resistance to the frontline antibiotic rifampicin constitutes a challenge to the treatment and control of tuberculosis. Here, we analyzed the mutational landscape of Mycobacterium smegmatis during long-term evolution with increasing concentrations of rifampicin, using a mutation accumulation assay combined with whole-genome sequencing. Antibiotic treatment enhanced the acquisition of mutations, doubling the genome-wide mutation rate of the wild-type cells. While antibiotic exposure led to extinction of almost all wild-type lines, the hypermutable phenotype of the ΔnucS mutant strain (noncanonical mismatch repair deficient) provided an efficient response to the antibiotic, leading to high rates of survival. This adaptative advantage resulted in the emergence of higher levels of rifampicin resistance, an accelerated acquisition of drug resistance mutations in rpoB (β RNA polymerase), and a wider diversity of evolutionary pathways that led to drug resistance. Finally, this approach revealed a subset of adaptive genes under positive selection with rifampicin that could be associated with the development of antibiotic resistance. IMPORTANCE Rifampicin is the most important first-line antibiotic against mycobacterial infections, including tuberculosis, one of the top causes of death worldwide. Acquisition of rifampicin resistance constitutes a major global public health problem that makes the control of the disease challenging. Here, we performed an experimental evolution assay under antibiotic selection to analyze the response and adaptation of mycobacteria, leading to the acquisition of rifampicin resistance. This approach explored the total number of mutations that arose in the mycobacterial genomes under long-term rifampicin exposure, using whole-genome sequencing. Our results revealed the effect of rifampicin at a genomic level, identifying different mechanisms and multiple pathways leading to rifampicin resistance in mycobacteria. Moreover, this study detected that an increase in the rate of mutations led to enhanced levels of drug resistance and survival. In summary, all of these results could be useful to understand and prevent the emergence of drug-resistant isolates in mycobacterial infections.

Vibrio vulnificus mutation rate: an in vitro approach

Vibrio vulnificus is a multi-host pathogenic species currently subdivided into five phylogenetic lineages (L) plus one pathovar with the ability to infect fish due to a transmissible virulence plasmid. This plasmid (or a fragment of it) has been transmitted between lineages within the species, contributing to the evolution of V. vulnificus. This study aimed to provide an experimental approximation to the V. vulnificus mutation rate by determining spontaneous mutation rates from bacterial cultures of representants of the different lineages by whole-genome sequencing. To this purpose, synonymous SNP differences, i.e., spontaneous mutation not subjected to the evolutive forces, between initial and final culture after serial growth were evaluated and used for mutation rate calculation.

Hybridization breaks species barriers in long-term coevolution of a cyanobacterial population

Bacterial species often undergo rampant recombination yet maintain cohesive genomic identity. Ecological differences can generate recombination barriers between species and sustain genomic clusters in the short term. But can these forces prevent genomic mixing during long-term coevolution? Cyanobacteria in Yellowstone hot springs comprise several diverse species that have coevolved for hundreds of thousands of years, providing a rare natural experiment. By analyzing more than 300 single-cell genomes, we show that despite each species forming a distinct genomic cluster, much of the diversity within species is the result of hybridization driven by selection, which has mixed their ancestral genotypes. This widespread mixing is contrary to the prevailing view that ecological barriers can maintain cohesive bacterial species and highlights the importance of hybridization as a source of genomic diversity.

DNA Data Storage

The demand for data storage is growing at an unprecedented rate, and current methods are not sufficient to accommodate such rapid growth due to their cost, space requirements, and energy consumption. Therefore, there is a need for a new, long-lasting data storage medium with high capacity, high data density, and high durability against extreme conditions. DNA is one of the most promising next-generation data carriers, with a storage density of 10¹⁹ bits of data per cubic centimeter, and its three-dimensional structure makes it about eight orders of magnitude denser than other storage media. DNA amplification during PCR or replication during cell proliferation enables the quick and inexpensive copying of vast amounts of data. In addition, DNA can possibly endure millions of years if stored in optimal conditions and dehydrated, making it useful for data storage. Numerous space experiments on microorganisms have also proven their extraordinary durability in extreme conditions, which suggests that DNA could be a durable storage medium for data. Despite some remaining challenges, such as the need to refine methods for the fast and error-free synthesis of oligonucleotides, DNA is a promising candidate for future data storage.

Case Studies in the Assessment of Microbial Fitness: Seemingly Subtle Changes Can Have Major Effects on Phenotypic Outcomes

Following the completion of an adaptive evolution experiment, fitness evaluations are routinely conducted to assess the magnitude of adaptation. In doing so, proper consideration should be given when determining the appropriate methods as trade-offs may exist between accuracy and throughput. Here, we present three instances in which small changes in the framework or execution of fitness evaluations significantly impacted the outcomes. The first case illustrates that discrepancies in fitness conclusions can arise depending on the approach to evaluating fitness, the culture vessel used, and the sampling method. The second case reveals that variations in environmental conditions can occur associated with culture vessel material. Specifically, these subtle changes can greatly affect microbial physiology leading to changes in the culture pH and distorting fitness measurements. Finally, the last case reports that heterogeneity in CFU formation time can result in inaccurate fitness conclusions. Based on each case, considerations and recommendations are presented for future adaptive evolution experiments.

Mutation bias and the predictability of evolution

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

Adaptive mechanisms of Bacillus to near space extreme environments

Earth's near space is an extreme atmosphere environment with high levels of radiation, low atmospheric pressure and dramatic temperature fluctuations. The region is above the flight altitude of aircraft but below the orbit of satellites, which has special and Mars-like conditions for investigating the survival and evolution of life. Technical limitations including flight devices, payloads and technologies/methodologies hinder microbiological research in near space. In this study, we investigated microbial survival and adaptive strategies in near space using a scientific balloon fight mission and multi-omics analyses. Methods for sample preparation, storage, protector and vessel were optimized to prepare the exposed microbial samples. After 3 h 17 min of exposure at a float altitude of ~32 km, only Bacillus strains were alive with survival efficiencies of 0-10^-6. Diverse mutants with significantly altered metabolites were generated, firstly proving that Earth's near space could be used as a new powerful microbial breeding platform. Multi-omics analyses of mutants revealed cascade changes at the genome, transcriptome and proteome levels. In response to environmental stresses, two mutants had similar proteome changes caused by different genomic mutations and mRNA expression levels. Metabolic network analysis combined with proteins' expression levels revealed that metabolic fluxes of EMP, PPP and purine synthesis-related pathways were significantly altered to increase/decrease inosine production. Further analysis showed that proteins related to translation, molecular chaperones, cell wall/membrane, sporulation, DNA replication/repair and anti-oxidation were significantly upregulated, enabling cells to efficiently repair DNA/protein damages and improve viability against environmental stress. Overall, these results revealed genetic and metabolic responses of Bacillus to the harsh conditions in near space, providing a research basis for bacterial adaptive mechanisms in extreme environments.

A primer to directed evolution: current methodologies and future directions

Directed evolution is one of the most powerful tools for protein engineering and functions by harnessing natural evolution, but on a shorter timescale. It enables the rapid selection of variants of biomolecules with properties that make them more suitable for specific applications. Since the first in vitro evolution experiments performed by Sol Spiegelman in 1967, a wide range of techniques have been developed to tackle the main two steps of directed evolution: genetic diversification (library generation), and isolation of the variants of interest. This review covers the main modern methodologies, discussing the advantages and drawbacks of each, and hence the considerations for designing directed evolution experiments. Furthermore, the most recent developments are discussed, showing how advances in the handling of ever larger library sizes are enabling new research questions to be tackled.

Analysis of genome instability and implications for the consequent phenotype in Plasmodium falciparum containing mutated MSH2-1 (P513T)

Malarial parasites exhibit extensive genomic plasticity, which induces the antigen diversification and the development of antimalarial drug resistance. Only a few studies have examined the genome maintenance mechanisms of parasites. The study aimed at elucidating the impact of a mutation in a DNA mismatch repair gene on genome stability by maintaining the mutant and wild-type parasites through serial in vitro cultures for approximately 400 days and analysing the subsequent spontaneous mutations. A P513T mutant of the DNA mismatch repair protein PfMSH2-1 from Plasmodium falciparum 3D7 was created. The mutation did not influence the base substitution rate but significantly increased the insertion/deletion (indel) mutation rate in short tandem repeats (STRs) and minisatellite loci. STR mutability was affected by allele size, genomic category and certain repeat motifs. In the mutants, significant telomere healing and homologous recombination at chromosomal ends caused extensive gene loss and generation of chimeric genes, resulting in large-scale chromosomal alteration. Additionally, the mutant showed increased tolerance to N-methyl-N'-nitro-N-nitrosoguanidine, suggesting that PfMSH2-1 was involved in recognizing DNA methylation damage. This work provides valuable insights into the role of PfMSH2-1 in genome stability and demonstrates that the genomic destabilization caused by its dysfunction may lead to antigen diversification.

Detection of DNA replication errors and 8-oxo-dGTP-mediated mutations in E. coli by Duplex DNA Sequencing

Mutation is a phenomenon inescapable for all life-forms, including bacteria. While bacterial mutation rates are generally low due to the operation of error-avoidance systems, sometimes they are elevated by many orders of magnitude. Such a state, known as a hypermutable state, can result from exposure to stress or to harmful environments. Studies of bacterial mutation frequencies and analysis of the precise types of mutations can provide insights into the mechanisms by which mutations occur and the possible involvement of error-avoidance pathways. Several approaches have been used for this, like reporter assays involving non-essential genes or mutation accumulation over multiple generations. However, these approaches give an indirect estimation, and a more direct approach for determining mutations is desirable. With the recent development of a DNA sequencing technique known as Duplex Sequencing, it is possible to detect rare variants in a population at a frequency of 1 in 107 base pairs or less. Here, we have applied Duplex Sequencing to study spontaneous mutations in E. coli. We also investigated the production of replication errors by using a mismatch-repair defective (mutL) strain as well as oxidative-stress associated mutations using a mutT-defective strain. For DNA from a wild-type strain we obtained mutant frequencies in the range of 10-7 to 10-8 depending on the specific base-pair substitution, but we argue that these mutants merely represent a background of the system, rather than mutations that occurred in vivo. In contrast, bona-fide in vivo mutations were identified for DNA from both the mutL and mutT strains, as indicated by specific increases in base substitutions that are fully consistent with their established in vivo roles. Notably, the data reproduce the specific context effects of in vivo mutations as well as the leading vs. lagging strand bias among DNA replication errors.

Phenotypic and genetic alterations of Burkholderia pseudomallei in patients during relapse and persistent infections

The bacterium Burkholderia pseudomallei is the causative agent of melioidosis, a severe tropical disease associated with high mortality and relapse and persistent infections. Treatment of melioidosis requires prolonged antibiotic therapy; however, little is known about relapse and persistent infections, particularly the phenotypic and genetic alterations of B. pseudomallei in patients. In this study, we performed pulsed-field gel electrophoresis (PFGE) to compare the bacterial genotype between the initial isolate and the subsequent isolate from each of 23 suspected recurrent and persistent melioidosis patients in Northeast Thailand. We used whole-genome sequencing (WGS) to investigate multilocus sequence types and genetic alterations of within-host strain pairs. We also investigated the bacterial phenotypes associated with relapse and persistent infections, including multinucleated giant cell (MNGC) formation efficiency and intracellular multiplication. We first identified 13 (1.2%) relapse, 7 (0.7%) persistent, and 3 (0.3%) reinfection patients from 1,046 survivors. Each of the 20 within-host strain pairs from patients with relapse and persistent infections shared the same genotype, suggesting that the subsequent isolates arise from the infecting isolate. Logistic regression analysis of clinical data revealed regimen and duration of oral antibiotic therapies as risk factors associated with relapse and persistent infections. WGS analysis demonstrated 17 within-host genetic alteration events in 6 of 20 paired isolates, including a relatively large deletion and 16 single-nucleotide polymorphism (stocktickerSNP) mutations distributed across 12 genes. In 1 of 20 paired isolates, we observed significantly increased cell-to-cell fusion and intracellular replication in the second isolate compared with the initial isolate from a patient with persistent infection. WGS analysis suggested that a non-synonymous mutation in the tssB-5 gene, which encoded an essential component of the type VI secretion system, may be associated with the increased intracellular replication and MNGC formation efficiency of the second isolate of the patient. This information provides insights into genetic and phenotypic alterations in B. pseudomallei in human melioidosis, which may represent a bacterial strategy for persistent and relapse infections.

Rates and spectra of de novo structural mutations in Chlamydomonas reinhardtii

Genetic variation originates from several types of spontaneous mutation, including single-nucleotide substitutions, short insertions and deletions (indels), and larger structural changes. Structural mutations (SMs) drive genome evolution and are thought to play major roles in evolutionary adaptation, speciation, and genetic disease, including cancers. Sequencing of mutation accumulation (MA) lines has provided estimates of rates and spectra of single-nucleotide and indel mutations in many species, yet the rate of new SMs is largely unknown. Here, we use long-read sequencing to determine the full mutation spectrum in MA lines derived from two strains (CC-1952 and CC-2931) of the green alga Chlamydomonas reinhardtii The SM rate is highly variable between strains and between MA lines, and SMs represent a substantial proportion of all mutations in both strains (CC-1952 6%; CC-2931 12%). The SM spectra differ considerably between the two strains, with almost all inversions and translocations occurring in CC-2931 MA lines. This variation is associated with heterogeneity in the number and type of active transposable elements (TEs), which comprise major proportions of SMs in both strains (CC-1952 22%; CC-2931 38%). In CC-2931, a Crypton and a previously undescribed type of DNA element have caused 71% of chromosomal rearrangements, whereas in CC-1952, a Dualen LINE is associated with 87% of duplications. Other SMs, notably large duplications in CC-2931, are likely products of various double-strand break repair pathways. Our results show that diverse types of SMs occur at substantial rates, and support prominent roles for SMs and TEs in evolution.

Adaptive Laboratory Evolution of Microorganisms: Methodology and Application for Bioproduction

Adaptive laboratory evolution (ALE) is a useful experimental methodology for fundamental scientific research and industrial applications to create microbial cell factories. By using ALE, cells are adapted to the environment that researchers set based on their objectives through the serial transfer of cell populations in batch cultivations or continuous cultures and the fitness of the cells (i.e., cell growth) under such an environment increases. Then, omics analyses of the evolved mutants, including genome sequencing, transcriptome, proteome and metabolome analyses, are performed. It is expected that researchers can understand the evolutionary adaptation processes, and for industrial applications, researchers can create useful microorganisms that exhibit increased carbon source availability, stress tolerance, and production of target compounds based on omics analysis data. In this review article, the methodologies for ALE in microorganisms are introduced. Moreover, the application of ALE for the creation of useful microorganisms as cell factories has also been introduced.

Noncanonical Mismatch Repair Protein NucS Modulates the Emergence of Antibiotic Resistance in Mycobacterium abscessus

NucS/EndoMS-dependent noncanonical mismatch repair (MMR) ensures the stability of genomic DNA in mycobacteria and acts as a guardian of the genome by preventing the accumulation of point mutations. In order to address whether the inactivation of noncanonical MMR could increase the acquisition of drug resistance by mutation, a ΔnucS strain was constructed and explored in the emerging pathogen Mycobacterium abscessus. Deletion of nucS resulted in a mutator phenotype with increased acquisition of resistance to macrolides and aminoglycosides, the two main groups of antimycobacterial agents for M. abscessus treatment, and also to second-line drugs such as fluoroquinolones. Inactivation of the noncanonical MMR in M. abscessus led to increases of 10- to 22-fold in the appearance of spontaneous mutants resistant to the macrolide clarithromycin and the aminoglycosides amikacin, gentamicin, and apramycin, compared with the wild-type strain. Furthermore, emergence of fluoroquinolone (ciprofloxacin) resistance was detected in a nucS-deficient strain but not in a wild-type M. abscessus strain. Acquired drug resistance to macrolides and aminoglycosides was analyzed through sequencing of the 23S rRNA gene rrl and the 16S rRNA gene rrs from independent drug-resistant colonies of both strains. When the acquisition of clarithromycin resistance was examined, a different mutational profile was detected in the M. abscessus ΔnucS strain compared with the wild-type one. To summarize, M. abscessus requires the NucS-dependent noncanonical MMR pathway to prevent the emergence of drug-resistant isolates by mutation. To our knowledge, this is the first report that reveals the role of NucS in a human pathogen, and these findings have potential implications for the treatment of M. abscessus infections. IMPORTANCE Chronic infections caused by M. abscessus are an emerging challenge in public health, posing a substantial health and economic burden, especially in patients with cystic fibrosis. Treatment of M. abscessus infections with antibiotics is particularly challenging, as its complex drug resistance mechanisms, including constitutive resistance through DNA mutation, lead to high rates of treatment failure. To decipher the evolution of antibiotic resistance in M. abscessus, we studied NucS-dependent noncanonical MMR, a unique DNA repair pathway involved in genomic maintenance. Inactivation of NucS is linked to the increase of DNA mutations (hypermutation), which can confer drug resistance. Our analysis detected increased acquisition of mutations conferring resistance to first-line and second-line antibiotics. We believe that this study will improve the knowledge of how this pathogen could evolve into an untreatable infectious agent, and it uncovers a role for hypermutators in chronic infectious diseases under antibiotic pressure.

Genome stability during serial subculturing in hyperepidemic multidrug-resistant Klebsiella pneumoniae and Escherichia coli

OBJECTIVES

Core-genome single nucleotide variant (cgSNV) analysis represents a powerful tool for epidemiological investigations of multidrug-resistant (MDR) bacteria. However, cgSNV thresholds to confirm whether isolates are the same clone are not formally defined.

METHODS

We implemented hybrid whole-genome sequencing to study the genomic changes of four MDR isolates belonging to hyperepidemic sequence types (STs) during 20 propagation steps (T20) on MacConkey and CHROMID^(R) ESBL plates. The following strains were analyzed: Klebsiella pneumoniae AE-2247421 (OXA-48/NDM-1-producing, ST101), K. pneumoniae MCL-2017-2 (CTX-M-15-producing, ST307), Escherichia coli Ec-042 (OXA-181-producing, ST410), and E. coli Ec-050 (NDM-5-producing, ST167). The genome assembly at T5 and T20 was compared to that at time point zero (T0) and to two reference genomes.

RESULTS

At T20, AE-2247421 lost the IncL bla_OXA-48-carrying plasmid when grown on CHROMID^(R) ESBL plates, while a large fragment encompassing bla_NDM-1 was lost from its IncC plasmid when grown on both plates. In contrast, no structural changes were noted for the other three strains. Regarding the cgSNVs, the following results were obtained at T5 and T20 (ranges considering the different agar plates and reference genomes): AE-2247421 (1-8 and 2-12 cgSNVs), MCL-2017-2 (both 1-2 cgSNVs), Ec-042 (both 0 cgSNVs), and Ec-050 (0-6 and 0-9 cgSNVs).

CONCLUSION

We showed that structural changes and accumulation of cgSNVs can occur in few propagation steps under laboratory conditions. These changes might also arise in the clinical context in a short time, especially under antibiotics treatment. This phenomenon should be carefully considered because it might affect the final interpretation of epidemiological genomic analyses.

Assessing the relative importance of bacterial resistance, persistence and hyper-mutation for antibiotic treatment failure

To curb the rising threat of antimicrobial resistance, we need to understand the routes to antimicrobial treatment failure. Bacteria can survive treatment by using both genetic and phenotypic mechanisms to diminish the effect of antimicrobials. We assemble empirical data showing that, for example, Pseudomonas aeruginosa infections frequently contain persisters, transiently non-growing cells unaffected by antibiotics (AB) and hyper-mutators, mutants with elevated mutation rates, and thus higher probability of genetic resistance emergence. Resistance, persistence and hyper-mutation dynamics are difficult to disentangle experimentally. Hence, we use stochastic population modelling and deterministic fitness calculations to investigate the relative importance of genetic and phenotypic mechanisms for immediate treatment failure and establishment of prolonged, chronic infections. We find that persistence causes 'hidden' treatment failure with very low cell numbers if antimicrobial concentrations prevent growth of genetically resistant cells. Persister cells can regrow after treatment is discontinued and allow for resistance evolution in the absence of AB. This leads to different mutational routes during treatment and relapse of an infection. By contrast, hyper-mutation facilitates resistance evolution during treatment, but rarely contributes to treatment failure. Our findings highlight the time and concentration dependence of different bacterial mechanisms to escape AB killing, which should be considered when designing 'failure-proof' treatments.

Detecting bacterial adaptation within individual microbiomes

The human microbiome harbours a large capacity for within-person adaptive mutations. Commensal bacterial strains can stably colonize a person for decades, and billions of mutations are generated daily within each person's microbiome. Adaptive mutations emerging during health might be driven by selective forces that vary across individuals, vary within an individual, or are completely novel to the human population. Mutations emerging within individual microbiomes might impact the immune system, the metabolism of nutrients or drugs, and the stability of the community to perturbations. Despite this potential, relatively little attention has been paid to the possibility of adaptive evolution within complex human-associated microbiomes. This review discusses the promise of studying within-microbiome adaptation, the conceptual and technical limitations that may have contributed to an underappreciation of adaptive de novo mutations occurring within microbiomes to date, and methods for detecting recent adaptive evolution. This article is part of a discussion meeting issue 'Genomic population structures of microbial pathogens'.

Environmental complexity is more important than mutation in driving the evolution of latent novel traits in E. coli

Recent experiments show that adaptive Darwinian evolution in one environment can lead to the emergence of multiple new traits that provide no immediate benefit in this environment. Such latent non-adaptive traits, however, can become adaptive in future environments. We do not know whether mutation or environment-driven selection is more important for the emergence of such traits. To find out, we evolve multiple wild-type and mutator E. coli populations under two mutation rates in simple (single antibiotic) environments and in complex (multi-antibiotic) environments. We then assay the viability of evolved populations in dozens of new environments and show that all populations become viable in multiple new environments different from those they had evolved in. The number of these new environments increases with environmental complexity but not with the mutation rate. Genome sequencing demonstrates the reason: Different environments affect pleiotropic mutations differently. Our experiments show that the selection pressure provided by an environment can be more important for the evolution of novel traits than the mutational supply experienced by a wild-type and a mutator strain of E. coli.

Novel traits without immediate fitness benefit evolve frequently but we don’t know whether mutation or environment-driven selection drives this evolution. Here, using experimental evolution of E. coli populations, the authors demonstrate the importance of selection in the evolution of latent novel traits.

Mutagenic Characteristics of Six Heavy Metals in Escherichia coli: The Commonality and Specificity

The history of long-term environmental exposure to heavy metals can be recorded in the genome as sporadic and specific mutations. Variable environments introduce diverse and adaptive mutations to organisms. To reveal the information hidden in genomes about environmental exposure to heavy metals, we performed long-term mutation accumulation (MA) experiments with Escherichia coli, analyzed genomes from 36 populations across 1650 generations with 6 heavy metal exposure regimes (arsenic, cadmium, chromium, copper, nickel, and lead), and inferred metal-specific evolution modes at the genomic level. All heavy metals induced genetic mutations with a mean rate of 3.459 × 10-9 per nucleotide per generation. The mutational spectrum exhibited distinct signatures; however, heavy metals also shared common mutation signatures prominently associated with all cancer types. The mutated genes showed an average similarity of 54.4% within the same exposure regime, whereas only 38.8% between exposure regimes. In terms of biological insights, mutated genes were enriched to fundamental cellular processes such as metabolism, motility, and transport. Our study elucidates the mutagenic commonality and specificity of environmental heavy metals, which are highly specific at mutational features and locus, but conserved at gene and functional levels, and may play crucial roles in the convergence of adaptation to heavy metals.

Himalayan Microbiomes for Agro-environmental Sustainability: Current Perspectives and Future Challenges

The Himalayas are one of the most mystical, yet least studied terrains of the world. One of Earth's greatest multifaceted and diverse montane ecosystems is also one of the thirty-four global biodiversity hotspots of the world. These are supposed to have been uplifted about 60-70 million years ago and support, distinct environments, physiography, a variety of orogeny, and great biological diversity (plants, animals, and microbes). Microbes are the pioneer colonizer of the Himalayas that are involved in various bio-geological cycles and play various significant roles. The applications of Himalayan microbiomes inhabiting in lesser to greater Himalayas have been recognized. The researchers explored the applications of indigenous microbiomes in both agricultural and environmental sectors. In agriculture, microbiomes from Himalayan regions have been suggested as better biofertilizers and biopesticides for the crops growing at low temperature and mountainous areas as they help in the alleviation of cold stress and other biotic stresses. Along with alleviation of low temperature, Himalayan microbes also have the capability to enhance plant growth by availing the soluble form of nutrients like nitrogen, phosphorus, potassium, zinc, and iron. These microbes have been recognized for producing plant growth regulators (abscisic acid, auxin, cytokinin, ethylene, and gibberellins). These microbes have been reported for bioremediating the diverse pollutants (pesticides, heavy metals, and xenobiotics) for environmental sustainability. In the current perspectives, present review provides a detailed discussion on the ecology, biodiversity, and adaptive features of the native Himalayan microbiomes in view to achieve agro-environmental sustainability.

The effect of migration and variation on populations of Escherichia coli adapting to complex fluctuating environments

Migration, a critical evolutionary force, can have contrasting effects on adaptation. It can aid as well as impede adaptation. The effects of migration on microbial adaptation have been studied primarily in simple constant environments. Very little is known about the effects of migration on adaptation to complex, fluctuating environments. In our study, we subjected replicate populations of Escherichia coli, adapting to complex and unpredictably fluctuating environments to different proportions of clonal ancestral immigrants. Contrary to the results from simple/constant environments, the presence of clonal immigrants reduced all measured proxies of fitness. However, migration from a source population with a greater variance in fitness resulted in no change in fitness with respect to the no-migration control, except at the highest level of migration. Thus, the presence of variation in the immigrants could counter the adverse effects of migration in complex and unpredictably fluctuating environments. Our study demonstrates that the effects of migration are strongly dependent on the nature of the destination environment and the genetic makeup of immigrants. These results enhance our understanding of the influences of migrating populations, which could help better predict the consequences of migration.

Selection in a growing colony biases results of mutation accumulation experiments

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

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

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

Generation of an Escherichia coli strain growing on methanol via the ribulose monophosphate cycle

Methanol is a liquid with high energy storage capacity that holds promise as an alternative substrate to replace sugars in the biotechnology industry. It can be produced from CO₂ or methane and its use does not compete with food and animal feed production. However, there are currently only limited biotechnological options for the valorization of methanol, which hinders its widespread adoption. Here, we report the conversion of the industrial platform organism Escherichia coli into a synthetic methylotroph that assimilates methanol via the energy efficient ribulose monophosphate cycle. Methylotrophy is achieved after evolution of a methanol-dependent E. coli strain over 250 generations in continuous chemostat culture. We demonstrate growth on methanol and biomass formation exclusively from the one-carbon source by ¹³C isotopic tracer analysis. In line with computational modeling, the methylotrophic E. coli strain optimizes methanol oxidation by upregulation of an improved methanol dehydrogenase, increasing ribulose monophosphate cycle activity, channeling carbon flux through the Entner-Doudoroff pathway and downregulating tricarboxylic acid cycle enzymes. En route towards sustainable bioproduction processes, our work lays the foundation for the efficient utilization of methanol as the dominant carbon and energy resource.

Using one carbon compounds as feedstock is a promising approach in abating climate change. Here, the authors report the conversion of E. coli into a synthetic methylotroph that assimilates methanol via the ribulose monophosphate cycle and a set of distinctive mutations.

Biocides used as material preservatives modify rates of de novo mutation and horizontal gene transfer in bacteria

Antimicrobial resistance (AMR) is a global health problem with the environment being an important compartment for the evolution and transmission of AMR. Previous studies showed that de-novo mutagenesis and horizontal gene transfer (HGT) by conjugation or transformation - important processes underlying resistance evolution and spread - are affected by antibiotics, metals and pesticides. However, natural microbial communities are also frequently exposed to biocides used as material preservatives, but it is unknown if these substances induce mutagenesis and HGT. Here, we show that active substances used in material preservatives can increase rates of mutation and conjugation in a species- and substance-dependent manner, while rates of transformation are not increased. The bisbiguanide chlorhexidine digluconate, the quaternary ammonium compound didecyldimethylammonium chloride, the metal copper, the pyrethroid-insecticide permethrin, and the azole-fungicide propiconazole increase mutation rates in Escherichia coli, whereas no increases were identified for Bacillus subtilis and Acinetobacter baylyi. Benzalkonium chloride, chlorhexidine and permethrin increased conjugation in E. coli. Moreover, our results show a connection between the RpoS-mediated general stress and the RecA-linked SOS response with increased rates of mutation and conjugation, but not for all biocides. Taken together, our data show the importance of assessing the contribution of material preservatives on AMR evolution and spread.

Parallel Evolution towards Increased Motility in Long-Term Cultures of Escherichia coli, Even Though Motility was Not Required for Long-Term Survival

Escherichia coli can survive for long periods in batch culture in the laboratory, where they experience a stressful and heterogeneous environment. During this incubation, E. coli acquires mutations that are selected in response to this environment, ultimately leading to evolved populations that are better adapted to these complex conditions, which can lead to a better understanding of evolutionary mechanisms. Mutations in regulatory genes often play a role in adapting to heterogeneous environments. To identify such mutations, we examined transcriptional differences during log phase growth in unaged cells compared to those that had been aged for 10 days and regrown. We identified expression changes in genes involved in motility and chemotaxis after adaptation to long-term cultures. We hypothesized that aged populations would also have phenotypic changes in motility and that motility may play a role in survival and adaptation to long-term cultures. While aged populations did show an increase in motility, this increase was not essential for survival in long-term cultures. We identified mutations in the regulatory gene sspA and other genes that may contribute to the observed differences in motility. Taken together, these data provide an overall picture of the role of mutations in regulatory genes for adaptation while underscoring that all changes that occur during evolution in stressful environments are not necessarily adaptive. IMPORTANCE Understanding how bacteria adapt in long-term cultures aids in both better treatment options for bacterial infections and gives insight into the mechanisms involved in bacterial evolution. In the past, it has been difficult to study these organisms in their natural environments. By using experimental evolution in heterogeneous and stressful laboratory conditions, we can more closely mimic natural environments and examine evolutionary mechanisms. One way to observe these mechanisms is to look at transcriptomic and genomic data from cells adapted to these complex conditions. Here, we found that although aged cells increase motility, this increase is not essential for survival in these conditions. These data emphasize that not all changes that occur due to evolutionary processes are adaptive, but these observations could still lead to hypotheses about the causative mutations. The information gained here allow us to make inferences about general mechanisms underlying phenotypic changes due to evolution.

Determination of Mutation Rates with Two Symmetric and Asymmetric Mutation Types

We revisit our earlier paper, with two of the coauthors, in which we proposed an unbiased and consistent estimator μ^n for an unknown mutation rate μ of microorganisms. Previously, we proved that the associated sequence of estimators μ^n converges to μ almost surely pointwise on a nonextinct set Ω0. Here, we show that this sequence converges also in the mean square with respect to conditional probability measure P0·=P·∩Ω0/PΩ0 and that, with respect to P0, the estimator is asymptotically unbiased. We further assume that a microorganism can mutate or turn to a different variant of one of the two types. In particular, it can mean that bacteria under attack by a virus or chemical agent are either perishing or surviving, turning them to stronger variant. We propose estimators for their respective types and show that they are a.s. pointwise and L2-consistent and asymptotically unbiased with respect to measure P0.

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

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

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

The role of the Pseudomonas aeruginosa hypermutator phenotype on the shift from acute to chronic virulence during respiratory infection

Chronic respiratory infection (CRI) with Pseudomonas aeruginosa (Pa) presents many unique challenges that complicate treatment. One notable challenge is the hypermutator phenotype which is present in up to 60% of sampled CRI patient isolates. Hypermutation can be caused by deactivating mutations in DNA mismatch repair (MMR) genes including mutS, mutL, and uvrD. In vitro and in vivo studies have demonstrated hypermutator strains to be less virulent than wild-type Pa. However, patients colonized with hypermutators display poorer lung function and a higher incidence of treatment failure. Hypermutation and MMR-deficiency create increased genetic diversity and population heterogeneity due to elevated mutation rates. MMR-deficient strains demonstrate higher rates of mucoidy, a hallmark virulence determinant of Pa during CRI in cystic fibrosis patients. The mucoid phenotype results from simple sequence repeat mutations in the mucA gene made in the absence of functional MMR. Mutations in Pa are further increased in the absence of MMR, leading to microcolony biofilm formation, further lineage diversification, and population heterogeneity which enhance bacterial persistence and host immune evasion. Hypermutation facilitates the adaptation to the lung microenvironment, enabling survival among nutritional complexity and microaerobic or anaerobic conditions. Mutations in key acute-to-chronic virulence “switch” genes, such as retS, bfmS, and ampR, are also catalyzed by hypermutation. Consequently, strong positive selection for many loss-of-function pathoadaptive mutations is seen in hypermutators and enriched in genes such as lasR. This results in the characteristic loss of Pa acute infection virulence factors, including quorum sensing, flagellar motility, and type III secretion. Further study of the role of hypermutation on Pa chronic infection is needed to better inform treatment regimens against CRI with hypermutator strains.

Stop codon usage as a window into genome evolution: mutation, selection, biased gene conversion and the TAG paradox

Protein coding genes terminate with one of three stop codons (TAA, TGA, or TAG) that, like synonymous codons, are not employed equally. With TGA and TAG having identical nucleotide content, analysis of their differential usage provides an unusual window into the forces operating on what are ostensibly functionally identical residues. Across genomes and between isochores within the human genome, TGA usage increases with G + C content but, with a common G + C → A + T mutation bias, this cannot be explained by mutation bias-drift equilibrium. Increased usage of TGA in G + C-rich genomes or genomic regions is also unlikely to reflect selection for the optimal stop codon, as TAA appears to be universally optimal, probably because it has the lowest read-through rate. Despite TAA being favored by selection and mutation bias, as with codon usage bias G + C pressure is the prime determinant of between-species TGA usage trends. In species with strong G + C-biased gene conversion (gBGC), such as mammals and birds, the high usage and conservation of TGA is best explained by an A + T → G + C repair bias. How to explain TGA enrichment in other G + C-rich genomes is less clear. Enigmatically, across bacterial and archaeal species and between human isochores TAG usage is mostly unresponsive to G + C pressure. This unresponsiveness we dub the TAG paradox as currently no mutational, selective, or gBGC model provides a well-supported explanation. That TAG does increase with G + C usage across eukaryotes makes the usage elsewhere yet more enigmatic. We suggest resolution of the TAG paradox may provide insights into either an unknown but common selective preference (probably at the DNA/RNA level) or an unrecognized complexity to the action of gBGC.

A limit on the evolutionary rescue of an Antarctic bacterium from rising temperatures

Climate change is gradual, but it can also cause brief extreme heat waves that can exceed the upper thermal limit of any one organism. To study the evolutionary potential of upper thermal tolerance, we evolved the cold-adapted Antarctic bacterium Pseudoalteromonas haloplanktis to survive at 30°C, beyond its ancestral thermal limit. This high-temperature adaptation occurred rapidly and in multiple populations. It involved genomic changes that occurred in a highly parallel fashion and mitigated the effects of protein misfolding. However, it also confronted a physiological limit, because populations failed to grow beyond 30°C. Our experiments aimed to facilitate evolutionary rescue by using a small organism with large populations living at temperatures several degrees below their upper thermal limit. Larger organisms with smaller populations and living at temperatures closer to their upper thermal tolerances are even more likely to go extinct during extreme heat waves.

Comparative overview of emerging RNA viruses: Epidemiology, pathogenesis, diagnosis and current treatment

From many decades, emerging infections have threatened humanity. The pandemics caused by different CoVs have already claimed and will continue to claim millions of lives. The SARS, Ebola, MERS epidemics and the most recent emergence of COVID-19 pandemic have threatened populations across borders. Since a highly pathogenic CoV has been evolved into the human population in the twenty-first century known as SARS, scientific advancements and innovative methods to tackle these viruses have increased in order to improve response preparedness towards the unpredictable threat posed by these rapidly emerging pathogens. Recently published review articles on SARS-CoV-2 have mainly focused on its pathogenesis, epidemiology and available treatments. However, in this review, we have done a systematic comparison of all three CoVs i.e., SARS, MERS and SARS-CoV-2 along with Ebola and Zika in terms of their epidemiology, virology, clinical features and current treatment strategies. This review focuses on important emerging RNA viruses starting from Zika, Ebola and the CoVs which include SARS, MERS and SARS-CoV-2. Each of these viruses has been elaborated on the basis of their epidemiology, virulence, transmission and treatment. However, special attention has been given to SARS-CoV-2 and the disease caused by it i.e., COVID-19 due to current havoc caused worldwide. At the end, insights into the current understanding of the lessons learned from previous epidemics to combat emerging CoVs have been described. The travel-related viral spread, the unprecedented nosocomial outbreaks and the high case-fatality rates associated with these highly transmissible and pathogenic viruses highlight the need for new prophylactic and therapeutic actions which include but are not limited to clinical indicators, contact tracing, and laboratory investigations as important factors that need to be taken into account in order to arrive at the final conclusion.

•

Recently published review articles on SARS-CoV-2 have mainly focused on its pathogenesis, epidemiology and available treatments.

•

The pandemics caused by different CoVs have already claimed and will continue to claim millions of lives.

•

This review focuses on important emerging RNA viruses starting from Zika, Ebola and the CoVs which include SARS, MERS and SARS-CoV-2.

•

Globally, numerous studies and researchers have recently started fighting this virus.

Dynamics of Mismatch and Alternative Excision-Dependent Repair in Replicating Bacillus subtilis DNA Examined Under Conditions of Neutral Selection

Spontaneous DNA deamination is a potential source of transition mutations. In Bacillus subtilis, EndoV, a component of the alternative excision repair pathway (AER), counteracts the mutagenicity of base deamination-induced mispairs. Here, we report that the mismatch repair (MMR) system, MutSL, prevents the harmful effects of HNO2, a deaminating agent of Cytosine (C), Adenine (A), and Guanine (G). Using Maximum Depth Sequencing (MDS), which measures mutagenesis under conditions of neutral selection, in B. subtilis strains proficient or deficient in MutSL and/or EndoV, revealed asymmetric and heterogeneous patterns of mutations in both DNA template strands. While the lagging template strand showed a higher frequency of C → T substitutions; G → A mutations, occurred more frequently in the leading template strand in different genetic backgrounds. In summary, our results unveiled a role for MutSL in preventing the deleterious effects of base deamination and uncovered differential patterns of base deamination processing by the AER and MMR systems that are influenced by the sequence context and the replicating DNA strand.

A 16th century Escherichia coli draft genome associated with an opportunistic bile infection

Escherichia coli – one of the most characterized bacteria and a major public health concern – remains invisible across the temporal landscape. Here, we present the meticulous reconstruction of the first ancient E. coli genome from a 16^th century gallstone from an Italian mummy with chronic cholecystitis. We isolated ancient DNA and reconstructed the ancient E. coli genome. It consisted of one chromosome of 4446 genes and two putative plasmids with 52 genes. The E. coli strain belonged to the phylogroup A and an exceptionally rare sequence type 4995. The type VI secretion system component genes appears to be horizontally acquired from Klebsiella aerogenes, however we could not identify any pathovar specific genes nor any acquired antibiotic resistances. A sepsis mouse assay showed that a closely related contemporary E. coli strain was avirulent. Our reconstruction of this ancient E. coli helps paint a more complete picture of the burden of opportunistic infections of the past.

Ancient DNA from an Italian mummy’s gallstone provides insight into opportunistic E. coli infection.

Machine Learning Guided Batched Design of a Bacterial Ribosome Binding Site

Optimization of gene expression levels is an essential part of the organism design process. Fine control of this process can be achieved by engineering transcription and translation control elements, including the ribosome binding site (RBS). Unfortunately, the design of specific genetic parts remains challenging because of the lack of reliable design methods. To address this problem, we have created a machine learning guided Design-Build-Test-Learn (DBTL) cycle for the experimental design of bacterial RBSs to demonstrate how small genetic parts can be reliably designed using relatively small, high-quality data sets. We used Gaussian Process Regression for the Learn phase of the cycle and the Upper Confidence Bound multiarmed bandit algorithm for the Design of genetic variants to be tested in vivo. We have integrated these machine learning algorithms with laboratory automation and high-throughput processes for reliable data generation. Notably, by Testing a total of 450 RBS variants in four DBTL cycles, we have experimentally validated RBSs with high translation initiation rates equaling or exceeding our benchmark RBS by up to 34%. Overall, our results show that machine learning is a powerful tool for designing RBSs, and they pave the way toward more complicated genetic devices.

Evolutionary action of mutations reveals antimicrobial resistance genes in Escherichia coli

Since antibiotic development lags, we search for potential drug targets through directed evolution experiments. A challenge is that many resistance genes hide in a noisy mutational background as mutator clones emerge in the adaptive population. Here, to overcome this noise, we quantify the impact of mutations through evolutionary action (EA). After sequencing ciprofloxacin or colistin resistance strains grown under different mutational regimes, we find that an elevated sum of the evolutionary action of mutations in a gene identifies known resistance drivers. This EA integration approach also suggests new antibiotic resistance genes which are then shown to provide a fitness advantage in competition experiments. Moreover, EA integration analysis of clinical and environmental isolates of antibiotic resistant of E. coli identifies gene drivers of resistance where a standard approach fails. Together these results inform the genetic basis of de novo colistin resistance and support the robust discovery of phenotype-driving genes via the evolutionary action of genetic perturbations in fitness landscapes.

Complex Ecotype Dynamics Evolve in Response to Fluctuating Resources

Ecotypic diversification and its associated cooperative behaviors are frequently observed in natural microbial populations whose access to resources is often sporadic. However, the extent to which fluctuations in resource availability influence the emergence of cooperative ecotypes is not fully understood. To determine how exposure to repeated resource limitation affects the establishment and long-term maintenance of ecotypes in a structured environment, we followed 32 populations of Escherichia coli evolving to either 1-day or 10-day feast/famine cycles for 900 days. Population-level analysis revealed that compared to populations evolving to 1-day cycles, 10-day populations evolved increased biofilm density, higher parallelism in mutational targets, and increased mutation rates. As previous investigations of evolution in structured environments have identified biofilm formation as the earliest observable phenotype associated with diversification of ecotypes, we revived cultures midway through the evolutionary process and conducted additional genomic, transcriptional, and phenotypic analyses of clones isolated from these evolving populations. We found not only that 10-day feast/famine cycles support multiple ecotypes but also that these ecotypes exhibit cooperative behavior. Consistent with the black queen hypothesis, or evolution of cooperation by gene loss, transcriptomic evidence suggests the evolution of bidirectional cross-feeding behaviors based on essential resources. These results provide insight into how analogous cooperative relationships may emerge in natural microbial communities. IMPORTANCE Despite regular feast and famine conditions representing an environmental pressure that is commonly encountered by microbial communities, the evolutionary outcomes of repeated cycles of feast and famine have been less studied. By experimentally evolving initially isogenic Escherichia coli populations to 10-day feast/famine cycles, we observed rapid diversification into ecotypes with evidence of bidirectional cross-feeding on costly resources and frequency-dependent fitness. Although unidirectional cross-feeding has been repeatedly observed to evolve in laboratory culture, most investigations of bidirectional cooperative behaviors in microbial populations have been conducted in engineered communities. This work demonstrates the de novo evolution of black queen relationships in a microbial population originating from a single ancestor, providing a model for investigation of the eco-evolutionary processes leading to mutualistic cooperation.

The microbial ecology of Escherichia coli in the vertebrate gut

Escherichia coli has a rich history as biology's 'rock star', driving advances across many fields. In the wild, E. coli resides innocuously in the gut of humans and animals but is also a versatile pathogen commonly associated with intestinal and extraintestinal infections and antimicrobial resistance-including large foodborne outbreaks such as the one that swept across Europe in 2011, killing 54 individuals and causing approximately 4000 infections and 900 cases of haemolytic uraemic syndrome. Given that most E. coli are harmless gut colonizers, an important ecological question plaguing microbiologists is what makes E. coli an occasionally devastating pathogen? To address this question requires an enhanced understanding of the ecology of the organism as a commensal. Here, we review how our knowledge of the ecology and within-host diversity of this organism in the vertebrate gut has progressed in the 137 years since E. coli was first described. We also review current approaches to the study of within-host bacterial diversity. In closing, we discuss some of the outstanding questions yet to be addressed and prospects for future research.

Efficient and Iterative Retron-Mediated in vivo Recombineering in E. coli

Recombineering is an important tool in gene editing, enabling fast, precise and highly specific in vivo modification of microbial genomes. Oligonucleotide-mediated recombineering via the in vivo production of single-stranded DNA can overcome the limitations of traditional recombineering methods that rely on the exogenous delivery of editing templates. By modifying a previously reported plasmid-based system for fully in vivo single-stranded DNA recombineering, we demonstrate iterative editing of independent loci by utilizing a temperature-sensitive origin of replication for easy curing of the editing plasmid from recombinant cells. Optimization of the promoters driving the expression of the system’s functional components, combined with targeted counterselection against unedited cells with Cas9 nuclease, enabled editing efficiencies of 90–100%. The addition of a dominant-negative mutL allele to the system allowed single-nucleotide edits that were otherwise unachievable due to mismatch repair. Finally, we tested alternative recombinases and found that efficiency significantly increased for some targets. Requiring only a single cloning step for retargeting, our system provides an easy-to-use method for rapid, efficient construction of desired mutants.

Graphical Abstract

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.

Landscape of mobile genetic elements and their antibiotic resistance cargo in prokaryotic genomes

Prokaryotic Mobile Genetic Elements (MGEs) such as transposons, integrons, phages and plasmids, play important roles in prokaryotic evolution and in the dispersal of cargo functions like antibiotic resistance. However, each of these MGE types is usually annotated and analysed individually, hampering a global understanding of phylogenetic and environmental patterns of MGE dispersal. We thus developed a computational framework that captures diverse MGE types, their cargos and MGE-mediated horizontal transfer events, using recombinases as ubiquitous MGE marker genes and pangenome information for MGE boundary estimation. Applied to ∼84k genomes with habitat annotation, we mapped 2.8 million MGE-specific recombinases to six operational MGE types, which together contain on average 13% of all the genes in a genome. Transposable elements (TEs) dominated across all taxa (∼1.7 million occurrences), outnumbering phages and phage-like elements (<0.4 million). We recorded numerous MGE-mediated horizontal transfer events across diverse phyla and habitats involving all MGE types, disentangled and quantified the extent of hitchhiking of TEs (17%) and integrons (63%) with other MGE categories, and established TEs as dominant carriers of antibiotic resistance genes. We integrated all these findings into a resource (proMGE.embl.de), which should facilitate future studies on the large mobile part of genomes and its horizontal dispersal.

The selective breeding and mutagenesis mechanism of high-yielding surfactin Bacillus subtilis strains with atmospheric and room temperature plasma

BACKGROUND

Surfactin, a good biological surfactant, is derived from the metabolites of microorganisms. However, the ability of natural strains to produce surfactin is low, and so the presented study aimed to use a novel mutagenesis technology to increase their yields.

RESULTS

Atmospheric and room temperature plasma (ARTP) was used to conduct mutation breeding of Bacillus subtilis CICC 10721, and a mutant strain M45 with a higher surfactin yield of 34.2% and a stable subculture was screened out. From the fermentation kinetics study, it was found that the maximum cell dry weight, maximum growth rate and surfactin synthesis parameters of the mutant strain M45 were all greater than that of the original strain. Scanning electron microscope and laser scanning confocal microscope observations showed that the spore morphology changed after ARTP treating, and the intracellular Ca²⁺ concentration of the mutant increased. Genome resequencing analysis showed that 66 single nucleotide poymorphism non-synonymous mutation sites occurred in M45, and the identification results of the fermentation broth extract from M45 showed that it is composed of C₁₂ -C₁₆ surfactin.

CONCLUSION

ARTP mutagenesis was found to change the morphology of bacteria, membrane permeability and genes related to the synthesis and secretion of surfactin. The present study provides a basis for industrial production of surfactin and an understanding of the mutagenesis mechanism. © 2021 Society of Chemical Industry.

Evolutionary Stability Optimizer (ESO): A Novel Approach to Identify and Avoid Mutational Hotspots in DNA Sequences While Maintaining High Expression Levels

Modern synthetic biology procedures rely on the ability to generate stable genetic constructs that keep their functionality over long periods of time. However, maintenance of these constructs requires energy from the cell and thus reduces the host’s fitness. Natural selection results in loss-of-functionality mutations that negate the expression of the construct in the population. Current approaches for the prevention of this phenomenon focus on either small-scale, manual design of evolutionary stable constructs or the detection of mutational sites with unstable tendencies. We designed the Evolutionary Stability Optimizer (ESO), a software tool that enables the large-scale automatic design of evolutionarily stable constructs with respect to both mutational and epigenetic hotspots and allows users to define custom hotspots to avoid. Furthermore, our tool takes the expression of the input constructs into account by considering the guanine-cytosine (GC) content and codon usage of the host organism, balancing the trade-off between stability and gene expression, allowing to increase evolutionary stability while maintaining the high expression. In this study, we present the many features of the ESO and show that it accurately predicts the evolutionary stability of endogenous genes. The ESO was created as an easy-to-use, flexible platform based on the notion that directed genetic stability research will continue to evolve and revolutionize current applications of synthetic biology. The ESO is available at the following link: https://www.cs.tau.ac.il/~tamirtul/ESO/.

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

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

Mutation bias shapes the spectrum of adaptive substitutions

How do mutational biases influence the process of adaptation? A common assumption is that selection alone determines the course of adaptation from abundant preexisting variation. Yet, theoretical work shows broad conditions under which the mutation rate to a given type of variant strongly influences its probability of contributing to adaptation. Here we introduce a statistical approach to analyzing how mutation shapes protein sequence adaptation. Using large datasets from three different species, we show that the mutation spectrum has a proportional influence on the types of changes fixed in adaptation. We also show via computer simulations that a variety of factors can influence how closely the spectrum of adaptive substitutions reflects the spectrum of variants introduced by mutation.

Evolutionary adaptation often occurs by the fixation of beneficial mutations. This mode of adaptation can be characterized quantitatively by a spectrum of adaptive substitutions, i.e., a distribution for types of changes fixed in adaptation. Recent work establishes that the changes involved in adaptation reflect common types of mutations, raising the question of how strongly the mutation spectrum shapes the spectrum of adaptive substitutions. We address this question with a codon-based model for the spectrum of adaptive amino acid substitutions, applied to three large datasets covering thousands of amino acid changes identified in natural and experimental adaptation in Saccharomyces cerevisiae, Escherichia coli, and Mycobacterium tuberculosis. Using species-specific mutation spectra based on prior knowledge, we find that the mutation spectrum has a proportional influence on the spectrum of adaptive substitutions in all three species. Indeed, we find that by inferring the mutation rates that best explain the spectrum of adaptive substitutions, we can accurately recover the species-specific mutation spectra. However, we also find that the predictive power of the model differs substantially between the three species. To better understand these differences, we use population simulations to explore the factors that influence how closely the spectrum of adaptive substitutions mirrors the mutation spectrum. The results show that the influence of the mutation spectrum decreases with increasing mutational supply (Nμ) and that predictive power is strongly affected by the number and diversity of beneficial mutations.

Prochlorococcus have low global mutation rate and small effective population size

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

Mutation bias reflects natural selection in Arabidopsis thaliana

Since the first half of the twentieth century, evolutionary theory has been dominated by the idea that mutations occur randomly with respect to their consequences1. Here we test this assumption with large surveys of de novo mutations in the plant Arabidopsis thaliana. In contrast to expectations, we find that mutations occur less often in functionally constrained regions of the genome-mutation frequency is reduced by half inside gene bodies and by two-thirds in essential genes. With independent genomic mutation datasets, including from the largest Arabidopsis mutation accumulation experiment conducted to date, we demonstrate that epigenomic and physical features explain over 90% of variance in the genome-wide pattern of mutation bias surrounding genes. Observed mutation frequencies around genes in turn accurately predict patterns of genetic polymorphisms in natural Arabidopsis accessions (r = 0.96). That mutation bias is the primary force behind patterns of sequence evolution around genes in natural accessions is supported by analyses of allele frequencies. Finally, we find that genes subject to stronger purifying selection have a lower mutation rate. We conclude that epigenome-associated mutation bias2 reduces the occurrence of deleterious mutations in Arabidopsis, challenging the prevailing paradigm that mutation is a directionless force in evolution.