Compensatory evolution of gene regulation in response to stress by Escherichia coli lacking RpoS

Adaptive laboratory evolution in a novel parallel shaken pH-auxostat

Adaptive laboratory evolution (ALE) is a widely used microbial strain development and optimization method. ALE experiments, to select for faster-growing strains, are commonly performed as serial batch cultivations in shake flasks, serum bottles, or microtiter plates or as continuous cultivations in bioreactors on a laboratory scale. To combine the advantages of higher throughput in parallel shaken cultures with continuous fermentations for conducting ALE experiments, a new Continuous parallel shaken pH-auxostat (CPA) was developed. The CPA consists of six autonomous parallel shaken cylindrical reactors, equipped with real-time pH control of the culture medium. The noninvasive pH measurement and control are realized by biocompatible pH sensor spots and a programmable pump module, to adjust the dilution rate of fresh medium for each reactor separately. Two different strains of the methylotrophic yeast Ogataea polymorpha were used as microbial model systems for parallel chemostat and pH-auxostat cultivations. During cultivation, the medium is acidified by the microbial activity of the yeast. For pH-auxostat cultivations, the growth-dependent acidification triggers the addition of fresh feed medium into the reactors, leading to a pH increase and thereby to the control of the pH to a predetermined set value. By controlling the pH to a predetermined set value, the dilution rate of the continuous cultivation is adjusted to values close to the washout point, in the range of the maximum specific growth rate of the yeast. The pH control was optimized by conducting a step-response experiment and obtaining tuned PI controller parameters by the Chien-Hrones-Reswick (CHR) PID tuning method. Two pH-auxostat cultivations were performed with two different O. polymorpha strains at high dilution rates for up to 18 days. As a result, up to 4.8-fold faster-growing strains were selected. The increased specific maximum growth rates of the selected strains were confirmed in subsequent batch cultivations.

RpoS and the bacterial general stress response

SUMMARYThe general stress response (GSR) is a widespread strategy developed by bacteria to adapt and respond to their changing environments. The GSR is induced by one or multiple simultaneous stresses, as well as during entry into stationary phase and leads to a global response that protects cells against multiple stresses. The alternative sigma factor RpoS is the central GSR regulator in E. coli and conserved in most γ-proteobacteria. In E. coli, RpoS is induced under conditions of nutrient deprivation and other stresses, primarily via the activation of RpoS translation and inhibition of RpoS proteolysis. This review includes recent advances in our understanding of how stresses lead to RpoS induction and a summary of the recent studies attempting to define RpoS-dependent genes and pathways.

Advance of tolerance engineering on microbes for industrial production

Industrial microbes have become the core of biological manufacturing, which utilized as the cell factory for production of plenty of chemicals, fuels and medicine. However, the challenge that the extreme stress conditions exist in production is unavoidable for cell factory. Consequently, to enhance robustness of the chassis cell lays the foundation for development of bio-manufacturing. Currently, the researches on cell tolerance covered various aspects, involving reshaping regulatory network, cell membrane modification and other stress response. In fact, the strategies employed to improve cell robustness could be summarized into two directions, irrational engineering and rational engineering. In this review, the metabolic engineering technologies on enhancement of microbe tolerance to industrial conditions are summarized. Meanwhile, the novel thoughts emerged with the development of biological instruments and synthetic biology are discussed.

The timing of transcription of RpoS-dependent genes varies across multiple stresses in Escherichia coli K-12

IMPORTANCE

Bacteria adapt to changing environments by altering the transcription of their genes. Specific proteins can regulate these changes. This study explored how a single protein called RpoS controls how many genes change expression during adaptation to three stresses. We found that: (i) RpoS is responsible for activating different genes in different stresses; (ii) that during a stress, the timing of gene activation depends on the what stress it is; and (iii) that how much RpoS a gene needs in order to be activated can predict when that gene will be activated during the stress of stationary phase.

Hypermutator emergence in experimental Escherichia coli populations is stress-type dependent

Genotypes exhibiting an increased mutation rate, called hypermutators, can propagate in microbial populations because they can have an advantage due to the higher supply of beneficial mutations needed for adaptation. Although this is a frequently observed phenomenon in natural and laboratory populations, little is known about the influence of parameters such as the degree of maladaptation, stress intensity, and the genetic architecture for adaptation on the emergence of hypermutators. To address this knowledge gap, we measured the emergence of hypermutators over ~1,000 generations in experimental Escherichia coli populations exposed to different levels of osmotic or antibiotic stress. Our stress types were chosen based on the assumption that the genetic architecture for adaptation differs between them. Indeed, we show that the size of the genetic basis for adaptation is larger for osmotic stress compared to antibiotic stress. During our experiment, we observed an increased emergence of hypermutators in populations exposed to osmotic stress but not in those exposed to antibiotic stress, indicating that hypermutator emergence rates are stress type dependent. These results support our hypothesis that hypermutator emergence is linked to the size of the genetic basis for adaptation. In addition, we identified other parameters that covaried with stress type (stress level and IS transposition rates) that might have contributed to an increased hypermutator provision and selection. Our results provide a first comparison of hypermutator emergence rates under varying stress conditions and point towards complex interactions of multiple stress-related factors on the evolution of mutation rates.

Application of Adaptive Laboratory Evolution in Lipid and Terpenoid Production in Yeast and Microalgae

Due to the complexity of metabolic and regulatory networks in microorganisms, it is difficult to obtain robust phenotypes through artificial rational design and genetic perturbation. Adaptive laboratory evolution (ALE) engineering plays an important role in the construction of stable microbial cell factories by simulating the natural evolution process and rapidly obtaining strains with stable traits through screening. This review summarizes the application of ALE technology in microbial breeding, describes the commonly used methods for ALE, and highlights the important applications of ALE technology in the production of lipids and terpenoids in yeast and microalgae. Overall, ALE technology provides a powerful tool for the construction of microbial cell factories, and it has been widely used in improving the level of target product synthesis, expanding the range of substrate utilization, and enhancing the tolerance of chassis cells. In addition, in order to improve the production of target compounds, ALE also employs environmental or nutritional stress strategies corresponding to the characteristics of different terpenoids, lipids, and strains.

Sequence-dependent model of genes with dual σ factor preference

Escherichia coli uses σ factors to quickly control large gene cohorts during stress conditions. While most of its genes respond to a single σ factor, approximately 5% of them have dual σ factor preference. The most common are those responsive to both σ⁷⁰, which controls housekeeping genes, and σ³⁸, which activates genes during stationary growth and stresses. Using RNA-seq and flow-cytometry measurements, we show that 'σ⁷⁰⁺³⁸ genes' are nearly as upregulated in stationary growth as 'σ³⁸ genes'. Moreover, we find a clear quantitative relationship between their promoter sequence and their response strength to changes in σ³⁸ levels. We then propose and validate a sequence dependent model of σ⁷⁰⁺³⁸ genes, with dual sensitivity to σ³⁸ and σ⁷⁰, that is applicable in the exponential and stationary growth phases, as well in the transient period in between. We further propose a general model, applicable to other stresses and σ factor combinations. Given this, promoters controlling σ⁷⁰⁺³⁸ genes (and variants) could become important building blocks of synthetic circuits with predictable, sequence-dependent sensitivity to transitions between the exponential and stationary growth phases.

Adaptive laboratory evolution triggers pathogen-dependent broad-spectrum antimicrobial potency in Streptomyces

BACKGROUND

In the present study, adaptive laboratory evolution was used to stimulate antibiotic production in a Streptomyces strain JB140 (wild-type) exhibiting very little antimicrobial activity against bacterial pathogens. The seven different competition experiments utilized three serial passages (3 cycles of adaptation-selection of 15 days each) in which Streptomyces strain (wild-type) was challenged repeatedly to one (bi-culture) or two (tri-culture) or three (quadri-culture) target pathogens. The study demonstrates a simple laboratory model to study the adaptive potential of evolved phenotypes and genotypes in Streptomyces to induce antibiotic production.

RESULTS

Competition experiments resulted in the evolution of the wild-type Streptomyces strain JB140 into the seven unique mutant phenotypes that acquired the ability to constitutively exhibit increased antimicrobial activity against three bacterial pathogens Salmonella Typhi (NCIM 2051), Staphylococcus aureus (NCIM 2079), and Proteus vulgaris (NCIM 2027). The mutant phenotypes not only effectively inhibited the growth of the tested pathogens but were also observed to exhibit improved antimicrobial responses against one clinical multidrug-resistant (MDR) uropathogenic Escherichia coli (UPEC 1021) isolate. In contrast to the adaptively evolved mutants, only a weak antimicrobial activity was detected in the wild-type parental strain. To get molecular evidence of evolution, RAPD profiles of the wild-type Streptomyces and its evolved mutants were compared which revealed significant polymorphism among them.

CONCLUSION

The competition-based adaptive laboratory evolution method can constitute a platform for evolutionary engineering to select improved phenotypes (mutants) with increased antibacterial profiles against targeted pathogens.

Mobile genetic elements mediate the mixotrophic evolution of novel Alicyclobacillus species for acid mine drainage adaptation

Alicyclobacillus species inhabit diverse environments and have adapted to broad ranges of pH and temperature. However, their adaptive evolutions remain elusive, especially regarding the role of mobile genetic elements (MGEs). Here, we characterized the distributions and functions of MGEs in Alicyclobacillus species across five environments, including acid mine drainage (AMD), beverages, hot springs, sediments, and soils. Nine Alicyclobacillus strains were isolated from AMD and possessed larger genome sizes and more genes than those from other environments. Four AMD strains evolved to be mixotrophic and fell into distinctive clusters in phylogenetic tree. Four types of MGEs including genomic island (GI), insertion sequence (IS), prophage, and integrative and conjugative element (ICE) were widely distributed in Alicyclobacillus species. Further, AMD strains did not possess CRISPR-Cas systems, but had more GI, IS, and ICE, as well as more MGE-associated genes involved in the oxidation of iron and sulfide and the resistance of heavy metal and low temperature. These findings highlight the differences in phenotypes and genotypes between strains isolated from AMD and other environments and the important role of MGEs in rapid environment niche expansions.

Synthetic biology-driven microbial production of folates: Advances and perspectives

With the development and application of synthetic biology, significant progress has been made in the production of folate by microbial fermentation using cell factories, especially for using generally regarded as safe (GRAS) microorganism as production host. In this review, the physiological functions and applications of folates were firstly discussed. Second, the current advances of folate-producing GRAS strains development were summarized. Third, the applications of synthetic biology-based metabolic regulatory tools in GRAS strains were introduced, and the progress in the application of these tools for folate production were summarized. Finally, the challenges to folates efficient production and corresponding emerging strategies to overcome them by synthetic biology were discussed, including the construction of biosensors using tetrahydrofolate riboswitches to regulate metabolic pathways, adaptive evolution to overcome the flux limitations of the folate pathway. The combination of new strategies and tools of synthetic biology is expected to further improve the efficiency of microbial folate synthesis.

Effect of dithiocyano-methane on hexose monophosphate pathway in the respiratory metabolism of Escherichia coli

This paper studied the inhibitory effects of dithiocyano-methane (DM) on the glucose decomposition pathway in the respiratory metabolism of Escherichia coli. We investigated the effects of DM on the activities of key enzymes (ATPase and glucose-6-phosphate dehydrogenase, G6PDH), the levels of key product (nicotinamide adenosine denucleotide hydro-phosphoric acid, NADPH), and gene expression in the hexose monophosphate pathway (HMP). The results showed that the minimum inhibitory concentration (MIC) and the minimum bactericide concentration (MBC) of DM against the tested strains were 5.86 mg/L and 11.72 mg/L, respectively. Bacteria exposed to DM at MIC demonstrated an increase in bacterial ATPase and G6PDH activities, NADPH levels, and gene expression in the HMP pathway compared to bacteria in the control group, which could be interpreted as a behavioral response to stress introduced by DM. However, DM at a lethal concentration of 10 × MIC affected glucose decomposition by inhibiting mainly the HMP pathway in E. coli.

A glimpse of potential transposable element impact on adaptation of the industrial yeast Saccharomyces cerevisiae

The adapted industrial yeast strain Saccharomyces cerevisiae NRRL Y-50049 is able to in situ detoxify major toxic aldehyde compounds derived from sugar conversion of lignocellulosic biomass while producing ethanol. Pathway-based studies on its mechanisms of tolerance have been reported previously, however, little is known about transposable element (TE) involvement in its adaptation to inhibitory compounds. This work presents a comparative dynamic transcription expression analysis in response to a toxic treatment between Y-50049 and its progenitor, an industrial type strain NRRL Y-12632, using a time-course study. At least 77 TEs from Y-50049 showed significantly increased expression compared with its progenitor, especially during the late lag phase. Sequence analysis revealed significant differences in TE sequences between the two strains. Y-50049 was also found to have a transposons of yeast 2 (Ty2) long terminal repeat-linked YAT1 gene showing significantly higher copy number changes than its progenitor. These results raise awareness of potential TE involvement in the adaptation of industrial yeast to the tolerance of toxic chemicals.

Evolutionary Repair Experiments as a Window to the Molecular Diversity of Life

Comparative genomics reveals an unexpected diversity in the molecular mechanisms underlying conserved cellular functions, such as DNA replication and cytokinesis. However, the genetic bases and evolutionary processes underlying this 'molecular diversity' remain to be explained. Here, we review a tool to generate alternative mechanisms for conserved cellular functions and test hypotheses concerning the generation of molecular diversity - evolutionary repair experiments, in which laboratory microbial populations adapt in response to a genetic perturbation. We summarize the insights gained from evolutionary repair experiments, the spectrum and dynamics of compensatory mutations, and the alternative molecular mechanisms used to repair perturbed cellular functions. We relate these experiments to the modifications of conserved functions that have occurred outside the laboratory. We end by proposing strategies to improve evolutionary repair experiments as a tool to explore the molecular diversity of life.

Rapid Accumulation of Motility-Activating Mutations in Resting Liquid Culture of Escherichia coli

Expression of motility genes is a potentially beneficial but costly process in bacteria. Interestingly, many isolate strains of Escherichia coli possess motility genes but have lost the ability to activate them under conditions in which motility is advantageous, raising the question of how they respond to these situations. Through transcriptome profiling of strains in the E. coli single-gene knockout Keio collection, we noticed drastic upregulation of motility genes in many of the deletion strains compared to levels in their weakly motile parent strain (BW25113). We show that this switch to a motile phenotype is not a direct consequence of the genes deleted but is instead due to a variety of secondary mutations that increase the expression of the major motility regulator, FlhDC. Importantly, we find that this switch can be reproduced by growing poorly motile E. coli strains in nonshaking liquid medium overnight but not in shaking liquid medium. Individual isolates after the nonshaking overnight incubations acquired distinct mutations upstream of the flhDC operon, including different insertion sequence (IS) elements and, to a lesser extent, point mutations. The rapidity with which genetic changes sweep through the populations grown without shaking shows that poorly motile strains can quickly adapt to a motile lifestyle by genetic rewiring.IMPORTANCE The ability to tune gene expression in times of need outside preordained regulatory networks is an essential evolutionary process that allows organisms to survive and compete. Here, we show that upon overnight incubation in liquid medium without shaking, populations of largely nonmotile Escherichia coli bacteria can rapidly accumulate mutants that have constitutive motility. This effect contributes to widespread secondary mutations in the single-gene knockout library, the Keio collection. As a result, 49/71 (69%) of the Keio strains tested exhibited various degrees of motility, whereas their parental strain is poorly motile. These observations highlight the plasticity of gene expression even in the absence of preexisting regulatory programs and should raise awareness of procedures for handling laboratory strains of E. coli.

Within-Host Adaptation Mediated by Intergenic Evolution in Pseudomonas aeruginosa

Bacterial pathogens evolve during the course of infection as they adapt to the selective pressures that confront them inside the host. Identification of adaptive mutations and their contributions to pathogen fitness remains a central challenge. Although mutations can either target intergenic or coding regions in the pathogen genome, studies of host adaptation have focused predominantly on molecular evolution within coding regions, whereas the role of intergenic mutations remains unclear. Here, we address this issue and investigate the extent to which intergenic mutations contribute to the evolutionary response of a clinically important bacterial pathogen, Pseudomonas aeruginosa, to the host environment, and whether intergenic mutations have distinct roles in host adaptation. We characterize intergenic evolution in 44 clonal lineages of P. aeruginosa and identify 77 intergenic regions in which parallel evolution occurs. At the genetic level, we find that mutations in regions under selection are located primarily within regulatory elements upstream of transcriptional start sites. At the functional level, we show that some of these mutations both increase or decrease transcription of genes and are directly responsible for evolution of important pathogenic phenotypes including antibiotic sensitivity. Importantly, we find that intergenic mutations facilitate essential genes to become targets of evolution. In summary, our results highlight the evolutionary significance of intergenic mutations in creating host-adapted strains, and that intergenic and coding regions have different qualitative contributions to this process.

Rapid decline of bacterial drug-resistance in an antibiotic-free environment through phenotypic reversion

Antibiotic resistance typically induces a fitness cost that shapes the fate of antibiotic-resistant bacterial populations. However, the cost of resistance can be mitigated by compensatory mutations elsewhere in the genome, and therefore the loss of resistance may proceed too slowly to be of practical importance. We present our study on the efficacy and phenotypic impact of compensatory evolution in Escherichia coli strains carrying multiple resistance mutations. We have demonstrated that drug-resistance frequently declines within 480 generations during exposure to an antibiotic-free environment. The extent of resistance loss was found to be generally antibiotic-specific, driven by mutations that reduce both resistance level and fitness costs of antibiotic-resistance mutations. We conclude that phenotypic reversion to the antibiotic-sensitive state can be mediated by the acquisition of additional mutations, while maintaining the original resistance mutations. Our study indicates that restricting antimicrobial usage could be a useful policy, but for certain antibiotics only.

The emergence of adaptive laboratory evolution as an efficient tool for biological discovery and industrial biotechnology

Harnessing the process of natural selection to obtain and understand new microbial phenotypes has become increasingly possible due to advances in culturing techniques, DNA sequencing, bioinformatics, and genetic engineering. Accordingly, Adaptive Laboratory Evolution (ALE) experiments represent a powerful approach both to investigate the evolutionary forces influencing strain phenotypes, performance, and stability, and to acquire production strains that contain beneficial mutations. In this review, we summarize and categorize the applications of ALE to various aspects of microbial physiology pertinent to industrial bioproduction by collecting case studies that highlight the multitude of ways in which evolution can facilitate the strain construction process. Further, we discuss principles that inform experimental design, complementary approaches such as computational modeling that help maximize utility, and the future of ALE as an efficient strain design and build tool driven by growing adoption and improvements in automation.

Plasticity of Promoter-Core Sequences Allows Bacteria to Compensate for the Loss of a Key Global Regulatory Gene

Transcription regulatory networks (TRNs) are of central importance for both short-term phenotypic adaptation in response to environmental fluctuations and long-term evolutionary adaptation, with global regulatory genes often being targets of natural selection in laboratory experiments. Here, we combined evolution experiments, whole-genome resequencing, and molecular genetics to investigate the driving forces, genetic constraints, and molecular mechanisms that dictate how bacteria can cope with a drastic perturbation of their TRNs. The crp gene, encoding a major global regulator in Escherichia coli, was deleted in four different genetic backgrounds, all derived from the Long-Term Evolution Experiment (LTEE) but with different TRN architectures. We confirmed that crp deletion had a more deleterious effect on growth rate in the LTEE-adapted genotypes; and we showed that the ptsG gene, which encodes the major glucose-PTS transporter, gained CRP (cyclic AMP receptor protein) dependence over time in the LTEE. We then further evolved the four crp-deleted genotypes in glucose minimal medium, and we found that they all quickly recovered from their growth defects by increasing glucose uptake. We showed that this recovery was specific to the selective environment and consistently relied on mutations in the cis-regulatory region of ptsG, regardless of the initial genotype. These mutations affected the interplay of transcription factors acting at the promoters, changed the intrinsic properties of the existing promoters, or produced new transcription initiation sites. Therefore, the plasticity of even a single promoter region can compensate by three different mechanisms for the loss of a key regulatory hub in the E. coli TRN.

Novel approach to quantitative spatial gene expression uncovers genetic stochasticity in the developing Drosophila eye

Robustness in development allows for the accumulation of genetically based variation in expression. However, this variation is usually examined in response to large perturbations, and examination of this variation has been limited to being spatial, or quantitative, but because of technical restrictions not both. Here we bridge these gaps by investigating replicated quantitative spatial gene expression using rigorous statistical models, in different genotypes, sexes, and species (Drosophila melanogaster and D. simulans). Using this type of quantitative approach with molecular developmental data allows for comparison among conditions, such as different genetic backgrounds. We apply this approach to the morphogenetic furrow, a wave of differentiation that patterns the developing eye disc. Within the morphogenetic furrow, we focus on four genes, hairy, atonal, hedgehog, and Delta. Hybridization chain reaction quantitatively measures spatial gene expression, co-staining for all four genes simultaneously. We find considerable variation in the spatial expression pattern of these genes in the eye between species, genotypes, and sexes. We also find that there has been evolution of the regulatory relationship between these genes, and that their spatial interrelationships have evolved between species. This variation has no phenotypic effect, and could be buffered by network thresholds or compensation from other genes. Both of these mechanisms could potentially be contributing to long term developmental systems drift.

Reduction of the fitness cost of antibiotic resistance caused by chromosomal mutations under poor nutrient conditions

The prevalence of antibiotic resistance in drinking water system is pressing public health risk. Antibiotic resistance conferred by chromosomal mutations often produces fitness cost, which may affect its spread and persistence. In this study, the rifampin-resistant strains were competed with their wild-type counterparts at different nutrient levels. It was observed that the ratio of the absolute number between resistant and wild-type cells quickly decreased under rich nutrient conditions, but it slowly reduced or remained stable in the poor nutrient medium. This finding suggested that poor nutrient conditions resulted in the reduction of fitness cost of antibiotic resistance, i.e. the resistant bacteria became more competitive. Implying mechanisms analysis found that the differences of metabolic activity between wild-type and rifampin-resistant strains was significant smaller (P < 0.05) at low nutrient levels. Additionally, distinguishable large colony size rifampin-resistant strains were observed during competition assay. DNA sequencing of RNA polymerase subunit genes further revealed that these colonies could be adaptive mutants from wild-type strain in rpoB gene. To our knowledge, this is the first study to reveal that the oligotrophic conditions facilitate the persistence of antibiotic resistance in drinking water by reducing the fitness cost of the resistant strains.

Experimental Design, Population Dynamics, and Diversity in Microbial Experimental Evolution

In experimental evolution, laboratory-controlled conditions select for the adaptation of species, which can be monitored in real time. Despite the current popularity of such experiments, nature's most pervasive biological force was long believed to be observable only on time scales that transcend a researcher's life-span, and studying evolution by natural selection was therefore carried out solely by comparative means. Eventually, microorganisms' propensity for fast evolutionary changes proved us wrong, displaying strong evolutionary adaptations over a limited time, nowadays massively exploited in laboratory evolution experiments. Here, we formulate a guide to experimental evolution with microorganisms, explaining experimental design and discussing evolutionary dynamics and outcomes and how it is used to assess ecoevolutionary theories, improve industrially important traits, and untangle complex phenotypes. Specifically, we give a comprehensive overview of the setups used in experimental evolution. Additionally, we address population dynamics and genetic or phenotypic diversity during evolution experiments and expand upon contributing factors, such as epistasis and the consequences of (a)sexual reproduction. Dynamics and outcomes of evolution are most profoundly affected by the spatiotemporal nature of the selective environment, where changing environments might lead to generalists and structured environments could foster diversity, aided by, for example, clonal interference and negative frequency-dependent selection. We conclude with future perspectives, with an emphasis on possibilities offered by fast-paced technological progress. This work is meant to serve as an introduction to those new to the field of experimental evolution, as a guide to the budding experimentalist, and as a reference work to the seasoned expert.

Random sequences rapidly evolve into de novo promoters

How new functions arise de novo is a fundamental question in evolution. We studied de novo evolution of promoters in Escherichia coli by replacing the lac promoter with various random sequences of the same size (~100 bp) and evolving the cells in the presence of lactose. We find that ~60% of random sequences can evolve expression comparable to the wild-type with only one mutation, and that ~10% of random sequences can serve as active promoters even without evolution. Such a short mutational distance between random sequences and active promoters may improve the evolvability, yet may also lead to accidental promoters inside genes that interfere with normal expression. Indeed, our bioinformatic analyses indicate that E. coli was under selection to reduce accidental promoters inside genes by avoiding promoter-like sequences. We suggest that a low threshold for functionality balanced by selection against undesired targets can increase the evolvability by making new beneficial features more accessible.

Bacterial promoters initiate gene transcription and have distinct sequence features. Here, the authors show that random sequences that contain no information are just on the verge of functioning as promoters in Escherichia coli.

Key Metabolites and Mechanistic Changes for Salt Tolerance in an Experimentally Evolved Sulfate-Reducing Bacterium, Desulfovibrio vulgaris

Rapid genetic and phenotypic adaptation of the sulfate-reducing bacterium Desulfovibrio vulgaris Hildenborough to salt stress was observed during experimental evolution. In order to identify key metabolites important for salt tolerance, a clone, ES10-5, which was isolated from population ES10 and allowed to experimentally evolve under salt stress for 5,000 generations, was analyzed and compared to clone ES9-11, which was isolated from population ES9 and had evolved under the same conditions for 1,200 generations. These two clones were chosen because they represented the best-adapted clones among six independently evolved populations. ES10-5 acquired new mutations in genes potentially involved in salt tolerance, in addition to the preexisting mutations and different mutations in the same genes as in ES9-11. Most basal abundance changes of metabolites and phospholipid fatty acids (PLFAs) were lower in ES10-5 than ES9-11, but an increase of glutamate and branched PLFA i17:1ω9c under high-salinity conditions was persistent. ES9-11 had decreased cell motility compared to the ancestor; in contrast, ES10-5 showed higher cell motility under both nonstress and high-salinity conditions. Both genotypes displayed better growth energy efficiencies than the ancestor under nonstress or high-salinity conditions. Consistently, ES10-5 did not display most of the basal transcriptional changes observed in ES9-11, but it showed increased expression of genes involved in glutamate biosynthesis, cation efflux, and energy metabolism under high salinity. These results demonstrated the role of glutamate as a key osmolyte and i17:1ω9c as the major PLFA for salt tolerance in D. vulgaris The mechanistic changes in evolved genotypes suggested that growth energy efficiency might be a key factor for selection.IMPORTANCE High salinity (e.g., elevated NaCl) is a stressor that affects many organisms. Salt tolerance, a complex trait involving multiple cellular pathways, is attractive for experimental evolutionary studies. Desulfovibrio vulgaris Hildenborough is a model sulfate-reducing bacterium (SRB) that is important in biogeochemical cycling of sulfur, carbon, and nitrogen, potentially for bio-corrosion, and for bioremediation of toxic heavy metals and radionuclides. The coexistence of SRB and high salinity in natural habitats and heavy metal-contaminated field sites laid the foundation for the study of salt adaptation of D. vulgaris Hildenborough with experimental evolution. Here, we analyzed a clone that evolved under salt stress for 5,000 generations and compared it to a clone evolved under the same condition for 1,200 generations. The results indicated the key roles of glutamate for osmoprotection and of i17:1ω9c for increasing membrane fluidity during salt adaptation. The findings provide valuable insights about the salt adaptation mechanism changes during long-term experimental evolution.

Subfunctionalization influences the expansion of bacterial multidrug antibiotic resistance

Background

Antibiotic resistance is a major problem for human health. Multidrug resistance efflux pumps, especially those of the Resistance-Nodulation-Cell Division (RND) family, are major contributors to high-level antibiotic resistance in Gram-negative bacteria. Most bacterial genomes contain several copies of the different classes of multidrug resistance efflux pumps. Gene duplication and gain of function by the duplicate copies of multidrug resistance efflux pump genes plays a key role in the expansion and diversification of drug-resistance mechanisms.

Results

We used two members of the Burkholderia RND superfamily as models to understand how duplication events affect the antibiotic resistance of these strains. First, we analyzed the conservation and distribution of these two RND systems and their regulators across the Burkholderia genus. Through genetic manipulations, we identified both the exact substrate range of these transporters and their eventual interchangeability. We also performed a directed evolution experiment, combined with next generation sequencing, to evaluate the role of antibiotics in the activation of the expression of these systems. Together, our results indicate that the first step to diversify the functions of these pumps arises from changes in their regulation (subfunctionalization) instead of functional mutations. Further, these pumps could rewire their regulation to respond to antibiotics, thus maintaining high genomic plasticity.

Conclusions

Studying the regulatory network that controls the expression of the RND pumps will help understand and eventually control the development and expansion of drug resistance.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4222-4) contains supplementary material, which is available to authorized users.

The adaptive landscape of wildtype and glycosylation-deficient populations of the industrial yeast Pichia pastoris

Background

The effects of long-term environmental adaptation and the implications of major cellular malfunctions are still poorly understood for non-model but biotechnologically relevant species. In this study we performed a large-scale laboratory evolution experiment with 48 populations of the yeast Pichia pastoris in order to establish a general adaptive landscape upon long-term selection in several glucose-based growth environments. As a model for a cellular malfunction the implications of OCH1 mannosyltransferase knockout-mediated glycosylation-deficiency were analyzed.

Results

In-depth growth profiling of evolved populations revealed several instances of genotype-dependent growth trade-off/cross-benefit correlations in non-evolutionary growth conditions. On the genome level a high degree of mutational convergence was observed among independent populations. Environment-dependent mutational hotspots were related to osmotic stress-, Rim - and cAMP signaling pathways. In agreement with the observed growth phenotypes, our data also suggest diverging compensatory mutations in glycosylation-deficient populations. High osmolarity glycerol (HOG) pathway loss-of-functions mutations, including genes such as SSK2 and SSK4, represented a major adaptive strategy during environmental adaptation. However, genotype-specific HOG-related mutations were predominantly observed in opposing environmental conditions. Surprisingly, such mutations emerged during salt stress adaptation in OCH1 knockout populations and led to growth trade-offs in non-adaptive conditions that were distinct from wildtype HOG-mutants. Further environment-dependent mutations were identified for a hitherto uncharacterized species-specific Gal4-like transcriptional regulator involved in environmental sensing.

Conclusion

We show that metabolic constraints such as glycosylation-deficiency can contribute to evolution on the molecular level, even in non-diverging growth environments. Our dataset suggests universal adaptive mechanisms involving cellular stress response and cAMP/PKA signaling but also the existence of highly species-specific strategies involving unique transcriptional regulators, improving our biological understanding of distinct Ascomycetes species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3952-7) contains supplementary material, which is available to authorized users.

Complex chromosomal neighborhood effects determine the adaptive potential of a gene under selection

How the organization of genes on a chromosome shapes adaptation is essential for understanding evolutionary paths. Here, we investigate how adaptation to rapidly increasing levels of antibiotic depends on the chromosomal neighborhood of a drug-resistance gene inserted at different positions of the Escherichia coli chromosome. Using a dual-fluorescence reporter that allows us to distinguish gene amplifications from other up-mutations, we track in real-time adaptive changes in expression of the drug-resistance gene. We find that the relative contribution of several mutation types differs systematically between loci due to properties of neighboring genes: essentiality, expression, orientation, termination, and presence of duplicates. These properties determine rate and fitness effects of gene amplification, deletions, and mutations compromising transcriptional termination. Thus, the adaptive potential of a gene under selection is a system-property with a complex genetic basis that is specific for each chromosomal locus, and it can be inferred from detailed functional and genomic data.

DOI:http://dx.doi.org/10.7554/eLife.25100.001

Transposable Elements Mediate Adaptive Debilitation of Flagella in Experimental Escherichia coli Populations

Although insertion sequence (IS) elements are generally considered genomic parasites, they can mediate adaptive genetic changes in bacterial genomes. We discovered that among 12 laboratory-evolved Escherichia coli populations, three had experienced at least six different IS1-mediated deletions of flagellar genes. These deletions all involved the master flagellar regulator flhDC, and as such completely incapacitate motility. Two lines of evidence strongly suggest that these deletions were adaptive in our evolution experiment: (1) parallel evolution in three independent populations is highly unlikely just by chance, and (2) one of these deletion mutations swept to fixation within ~1000 generations, which is over two million times faster than expected if this deletion was instead selectively neutral and thus evolving by genetic drift. Because flagella are energetically expensive to synthesize and operate, we suspect that debilitating their construction conferred a fitness advantage in our well-stirred evolution experiment. These findings underscore the important role that IS elements can play in mediating adaptive loss-of-function mutations in bacteria.

Inactivation of Transcriptional Regulators during Within-Household Evolution of Escherichia coli

We analyzed the within-household evolution of two household-associated Escherichia coli strains from pandemic clonal group ST131-H30, using isolates recovered from five individuals within two families, each of which had a distinct strain. Family 1's strain was represented by a urine isolate from the index patient (older sister) with recurrent cystitis and a blood isolate from her younger sister with fatal urosepsis. Family 2's strain was represented by a urine isolate from the index patient (father) with pyelonephritis and renal abscesses, blood and kidney drainage isolates from the daughter with emphysematous pyelonephritis, and urine and fecal isolates from the mother with cystitis. Collectively, the several variants of each family's strain had accumulated a total of 8 (family 1) and 39 (family 2) point mutations; no two isolates were identical. Of the 47 total mutations, 36 resulted in amino acid changes or truncation of coded proteins. Fourteen such mutations (39%) targeted genes encoding transcriptional regulators, and 9 (25%) involved DNA-binding transcription factors (TFs), which significantly exceeded the relative contribution of TF genes to the isolates' genomes (∼6%). At least one-half of the transcriptional regulator mutations were inactivating, based on phenotypic and/or transcriptional analysis. In particular, inactivating mutations in the global regulator LrhA (repressor of type 1 fimbriae and flagella) occurred in the blood isolates from both households and increased the virulence of E. coli strains in a murine sepsis model. The results indicate that E. coli undergoes adaptive evolution between and/or within hosts, generating subpopulations with distinctive phenotypes and virulence potential.IMPORTANCE The clonal evolution of bacterial strains associated with interhost transmission is poorly understood. We characterized the genome sequences of clonal descendants of two Escherichia coli strains, recovered at different time points from multiple individuals within two households who had different types of urinary tract infection. We found evidence that the E. coli strains underwent extensive mutational diversification between and within these individuals, driven disproportionately by inactivation of transcriptional regulators. In urosepsis isolates, the mutations observed in the global regulator LrhA increased bacterial virulence in a murine sepsis model. Our findings help in understanding the adaptive dynamics and strategies of E. coli during short-term natural evolution.

Environment-directed activation of the Escherichia coliflhDC operon by transposons

The flagellar system in Escherichia coli K12 is expressed under the control of the flhDC-encoded master regulator FlhDC. Transposition of insertion sequence (IS) elements to the upstream flhDC promoter region up-regulates transcription of this operon, resulting in a more rapid motility. Wang and Wood (ISME J 2011;5:1517-1525) provided evidence that insertion of IS5 into upstream activating sites occurs at higher rates in semi-solid agar media in which swarming behaviour is allowed as compared with liquid or solid media where swarming cannot occur. We confirm this conclusion and show that three IS elements, IS1, IS3 and IS5, transpose to multiple upstream sites within a 370 bp region of the flhDC operon control region. Hot spots for IS insertion correlate with positions of stress-induced DNA duplex destabilization (SIDD). We show that IS insertion occurs at maximal rates in 0.24 % agar, with rates decreasing dramatically with increasing or decreasing agar concentrations. In mixed cultures, we show that these mutations preferentially arise from the wild-type parent at frequencies of up to 3×10-3 cell-1 day-1 when the inoculated parental and co-existing IS-activated mutant cells are entering the stationary growth phase. We rigorously show that the apparent increased mutation frequencies cannot be accounted for by increased swimming or by increased growth under the selective conditions used. Thus, our data are consistent with the possibility that appropriate environmental conditions, namely those that permit but hinder flagellar rotation, result in the activation of a mutational pathway that involves IS element insertion upstream of the flhDC operon.

Regulation of solvent tolerance in Pseudomonas putida S12 mediated by mobile elements

Organic solvent‐tolerant bacteria are outstanding and versatile hosts for the bio‐based production of a broad range of generally toxic aromatic compounds. The energetically costly solvent tolerance mechanisms are subject to multiple levels of regulation, involving among other mobile genetic elements. The genome of the solvent‐tolerant Pseudomonas putida S12 contains many such mobile elements that play a major role in the regulation and adaptation to various stress conditions, including the regulation of expression of the solvent efflux pump SrpABC. We recently sequenced the genome of P. putida S12. Detailed annotation identified a threefold higher copy number of the mobile element ISS12 in contrast to earlier observations. In this study, we describe the mobile genetic elements and elaborate on the role of ISS12 in the establishment and maintenance of solvent tolerance in P. putida. We identified three different variants of ISS12 of which a single variant exhibits a high translocation rate. One copy of this variant caused a loss of solvent tolerance in the sequenced strain by disruption of srpA. Solvent tolerance could be restored by applying selective pressure, leading to a clean excision of the mobile element.

Genome-Wide Transcriptional Response to Varying RpoS Levels in Escherichia coli K-12

The alternative sigma factor RpoS is a central regulator of many stress responses in Escherichia coli The level of functional RpoS differs depending on the stress. The effect of these differing concentrations of RpoS on global transcriptional responses remains unclear. We investigated the effect of RpoS concentration on the transcriptome during stationary phase in rich media. We found that 23% of genes in the E. coli genome are regulated by RpoS, and we identified many RpoS-transcribed genes and promoters. We observed three distinct classes of response to RpoS by genes in the regulon: genes whose expression changes linearly with increasing RpoS level, genes whose expression changes dramatically with the production of only a little RpoS ("sensitive" genes), and genes whose expression changes very little with the production of a little RpoS ("insensitive"). We show that sequences outside the core promoter region determine whether an RpoS-regulated gene is sensitive or insensitive. Moreover, we show that sensitive and insensitive genes are enriched for specific functional classes and that the sensitivity of a gene to RpoS corresponds to the timing of induction as cells enter stationary phase. Thus, promoter sensitivity to RpoS is a mechanism to coordinate specific cellular processes with growth phase and may also contribute to the diversity of stress responses directed by RpoS.IMPORTANCE The sigma factor RpoS is a global regulator that controls the response to many stresses in Escherichia coli Different stresses result in different levels of RpoS production, but the consequences of this variation are unknown. We describe how changing the level of RpoS does not influence all RpoS-regulated genes equally. The cause of this variation is likely the action of transcription factors that bind the promoters of the genes. We show that the sensitivity of a gene to RpoS levels explains the timing of expression as cells enter stationary phase and that genes with different RpoS sensitivities are enriched for specific functional groups. Thus, promoter sensitivity to RpoS is a mechanism that coordinates specific cellular processes in response to stresses.

A novel point mutation in RpoB improves osmotolerance and succinic acid production in Escherichia coli

Background

Escherichia coli suffer from osmotic stress during succinic acid (SA) production, which reduces the performance of this microbial factory.

Results

Here, we report that a point mutation leading to a single amino acid change (D654Y) within the β-subunit of DNA-dependent RNA polymerase (RpoB) significantly improved the osmotolerance of E. coli. Importation of the D654Y mutation of RpoB into the parental strain, Suc-T110, increased cell growth and SA production by more than 40% compared to that of the control under high glucose osmolality. The transcriptome profile, determined by RNA-sequencing, showed two distinct stress responses elicited by the mutated RpoB that counterbalanced the osmotic stress. Under non-stressed conditions, genes involved in the synthesis and transport of compatible solutes such as glycine-betaine, glutamate or proline were upregulated even without osmotic stimulation, suggesting a “pre-defense” mechanism maybe formed in the rpoB mutant. Under osmotic stressed conditions, genes encoding diverse sugar transporters, which should be down-regulated in the presence of high osmotic pressure, were derepressed in the rpoB mutant. Additional genetic experiments showed that enhancing the expression of the mal regulon, especially for genes that encode the glycoporin LamB and maltose transporter, contributed to the osmotolerance phenotype.

Conclusions

The D654Y single amino acid substitution in RpoB rendered E. coli cells resistant to osmotic stress, probably due to improved cell growth and viability via enhanced sugar uptake under stressed conditions, and activated a potential “pre-defense” mechanism under non-stressed conditions. The findings of this work will be useful for bacterial host improvement to enhance its resistance to osmotic stress and facilitate bio-based organic acids production.

Electronic supplementary material

The online version of this article (doi:10.1186/s12896-017-0337-6) contains supplementary material, which is available to authorized users.

ATP-Dependent Persister Formation in Escherichia coli

Persisters are dormant variants that form a subpopulation of cells tolerant to antibiotics. Persisters are largely responsible for the recalcitrance of chronic infections to therapy. In Escherichia coli, one widely accepted model of persister formation holds that stochastic accumulation of ppGpp causes activation of the Lon protease that degrades antitoxins; active toxins then inhibit translation, resulting in dormant, drug-tolerant persisters. We found that various stresses induce toxin-antitoxin (TA) expression but that induction of TAs does not necessarily increase persisters. The 16S rRNA promoter rrnB P1 was proposed to be a persister reporter and an indicator of toxin activation regulated by ppGpp. Using fluorescence-activated cell sorting (FACS), we confirmed the enrichment for persisters in the fraction of rrnB P1-gfp dim cells; however, this is independent of toxin-antitoxins. rrnB P1 is coregulated by ppGpp and ATP. We show that rrnB P1 can report persisters in a relA/spoT deletion background, suggesting that rrnB P1 is a persister marker responding to ATP. Consistent with this finding, decreasing the level of ATP by arsenate treatment causes drug tolerance. Lowering ATP slows translation and prevents the formation of DNA double-strand breaks upon fluoroquinolone treatment. We conclude that variation in ATP levels leads to persister formation by decreasing the activity of antibiotic targets.

Persisters are a subpopulation of antibiotic-tolerant cells responsible for the recalcitrance of chronic infections. Our current understanding of persister formation is primarily based on studies of E. coli. The activation of toxin-antitoxin systems by ppGpp has become a widely accepted model for persister formation. In this study, we found that stress-induced activation of mRNA interferase-type toxins does not necessarily cause persister formation. We also found that the persister marker rrnB P1 reports persister cells because it detects a drop in cellular ATP levels. Consistent with this, lowering the ATP level decreases antibiotic target activity and, thus, leads to persister formation. We conclude that stochastic variation in ATP is the main mechanism of persister formation. A decrease in ATP provides a satisfactory explanation for the drug tolerance of persisters, since bactericidal antibiotics act by corrupting energy-dependent targets.

Genome-wide mapping of mutations at single-nucleotide resolution for protein, metabolic and genome engineering

Improvements in DNA synthesis and sequencing have underpinned comprehensive assessment of gene function in bacteria and eukaryotes. Genome-wide analyses require high-throughput methods to generate mutations and analyze their phenotypes, but approaches to date have been unable to efficiently link the effects of mutations in coding regions or promoter elements in a highly parallel fashion. We report that CRISPR-Cas9 gene editing in combination with massively parallel oligomer synthesis can enable trackable editing on a genome-wide scale. Our method, CRISPR-enabled trackable genome engineering (CREATE), links each guide RNA to homologous repair cassettes that both edit loci and function as barcodes to track genotype-phenotype relationships. We apply CREATE to site saturation mutagenesis for protein engineering, reconstruction of adaptive laboratory evolution experiments, and identification of stress tolerance and antibiotic resistance genes in bacteria. We provide preliminary evidence that CREATE will work in yeast. We also provide a webtool to design multiplex CREATE libraries.

Piecewise linear approximations to model the dynamics of adaptation to osmotic stress by food-borne pathogens

Addition of salt to food is one of the most ancient and most common methods of food preservation. However, little is known of how bacterial cells adapt to such conditions. We propose to use piecewise linear approximations to model the regulatory adaptation of Escherichiacoli to osmotic stress. We apply the method to eight selected genes representing the functions known to be at play during osmotic adaptation. The network is centred on the general stress response factor, sigma S, and also includes a module representing the catabolic repressor CRP-cAMP. Glutamate, potassium and supercoiling are combined to represent the intracellular regulatory signal during osmotic stress induced by salt. The output is a module where growth is represented by the concentration of stable RNAs and the transcription of the osmotic gene osmY. The time course of gene expression of transport of osmoprotectant represented by the symporter proP and of the osmY is successfully reproduced by the network. The behaviour of the rpoS mutant predicted by the model is in agreement with experimental data. We discuss the application of the model to food-borne pathogens such as Salmonella; although the genes considered have orthologs, it seems that supercoiling is not regulated in the same way. The model is limited to a few selected genes, but the regulatory interactions are numerous and span different time scales. In addition, they seem to be condition specific: the links that are important during the transition from exponential to stationary phase are not all needed during osmotic stress. This model is one of the first steps towards modelling adaptation to stress in food safety and has scope to be extended to other genes and pathways, other stresses relevant to the food industry, and food-borne pathogens. The method offers a good compromise between systems of ordinary differential equations, which would be unmanageable because of the size of the system and for which insufficient data are available, and the more abstract Boolean methods.

Memory and Combinatorial Logic Based on DNA Inversions: Dynamics and Evolutionary Stability

Genetic memory can be implemented using enzymes that catalyze DNA inversions, where each orientation corresponds to a "bit". Here, we use two DNA invertases (FimE and HbiF) that reorient DNA irreversibly between two states with opposite directionality. First, we construct memory that is set by FimE and reset by HbiF. Next, we build a NOT gate where the input promoter drives FimE and in the absence of signal the reverse state is maintained by the constitutive expression of HbiF. The gate requires ∼3 h to turn on and off. The evolutionary stabilities of these circuits are measured by passaging cells while cycling function. The memory switch is stable over 400 h (17 days, 14 state changes); however, the gate breaks after 54 h (>2 days) due to continuous invertase expression. Genome sequencing reveals that the circuit remains intact, but the host strain evolves to reduce invertase expression. This work highlights the need to evaluate the evolutionary robustness and failure modes of circuit designs, especially as more complex multigate circuits are implemented.

Re-engineering cellular physiology by rewiring high-level global regulatory genes

Knowledge of global regulatory networks has been exploited to rewire the gene control programmes of the model bacterium Salmonella enterica serovar Typhimurium. The product is an organism with competitive fitness that is superior to that of the wild type but tuneable under specific growth conditions. The paralogous hns and stpA global regulatory genes are located in distinct regions of the chromosome and control hundreds of target genes, many of which contribute to stress resistance. The locations of the hns and stpA open reading frames were exchanged reciprocally, each acquiring the transcription control signals of the other. The new strain had none of the compensatory mutations normally associated with alterations to hns expression in Salmonella; instead it displayed rescheduled expression of the stress and stationary phase sigma factor RpoS and its regulon. Thus the expression patterns of global regulators can be adjusted artificially to manipulate microbial physiology, creating a new and resilient organism.

Chloramphenicol Selection of IS10 Transposition in the cat Promoter Region of Widely Used Cloning Vectors

The widely used pSU8 family of cloning vectors is based on a p15A replicon and a chloramphenicol acetyltransferase (cat) gene conferring chloramphenicol resistance. We frequently observed an increase in the size of plasmids derived from these vectors. Analysis of the bigger molecular species shows that they have an IS10 copy inserted at a specific site between the promoter and the cat open reading frame. Promoter activity from both ends of IS10 has been reported, suggesting that the insertion events could lead to higher CAT production. Insertions were observed in certain constructions containing inserts that could lead to plasmid instability. To test the possibility that IS10 insertions were selected as a response to chloramphenicol selection, we have grown these constructs in the presence of different amounts of antibiotic and we observed that insertions arise promptly under higher chloramphenicol selective pressure. IS10 is present in many E. coli laboratory strains, so the possibility of insertion in constructions involving cat-containing vectors should be taken into account. Using lower chloramphenicol concentrations could solve this problem.

Lysine acetylation is a common post-translational modification of key metabolic pathway enzymes of the anaerobe Porphyromonas gingivalis

Porphyromonas gingivalis is a Gram-negative anaerobe considered to be a keystone pathogen in the development of the bacterial-associated inflammatory oral disease chronic periodontitis. Although post-translational modifications (PTMs) of proteins are commonly found to modify protein function in eukaryotes and prokaryotes, PTMs such as lysine acetylation have not been examined in P. gingivalis. Lysine acetylation is the addition of an acetyl group to a lysine which removes this amino acid's positive charge and can induce changes in a protein's secondary structure and reactivity. A proteomics based approach combining immune-affinity enrichment with high sensitivity Orbitrap mass spectrometry identified 130 lysine acetylated peptides from 92 P. gingivalis proteins. The majority of these peptides (71) were attributed to 45 proteins with predicted metabolic activity; these proteins could be mapped to several P. gingivalis metabolic pathways where enzymes catalysing sequential reactions within the same pathway were often found acetylated. In particular, the catabolic pathways of complex anaerobic fermentation of amino acids to produce energy had 12 enzymes lysine acetylated. The results suggest that lysine acetylation may be an important mechanism in metabolic regulation in P. gingivalis, which is vital for P. gingivalis survival and adaptation of its metabolism throughout infection. Statement of significance. Porphyromonas gingivalis is a keystone pathogen in the development of chronic periodontitis, an inflammatory disease of the supporting tissues of the teeth. The ability of the pathogen to induce dysbiosis and disease is related to an array of specific virulence factors and metabolic regulation that enables the bacterium to proliferate in an inflamed periodontal pocket. The mechanisms P. gingivalis uses to adapt to a changing and hostile environment are poorly understood and here we show, for the first time, that enzymes of critical metabolic pathways for energy production in this bacterium were acetylated on certain lysine residues. These enzymes were often found catalysing sequential reactions within the same catabolic pathway. The results suggest that lysine acetylation is an important mechanism of metabolic regulation in P. gingivalis vital for its adaptation and proliferation to produce disease.

Determination of single nucleotide variants in Escherichia coli DH5α by using short-read sequencing

Escherichia coli DH5α is a common laboratory strain that provides an important platform for routine use in cloning and synthetic biology applications. Many synthetic circuits have been constructed and successfully expressed in E. coli DH5α; however, its genome sequence has not been determined yet. Here, we determined E. coli DH5α genome sequence and identified genetic mutations that affect its phenotypic functions by using short-read sequencing. The sequencing results clearly described the genotypes of E. coli DH5α, which aid in further studies using the strain. Additionally, we observed 105 single nucleotide variants (SNVs), 83% of which were detected in protein-coding regions compared to the parental strain E. coli DH1. Interestingly, 23% of the protein-coding regions have mutations in their amino acid residues, whose biological functions were categorized into two-component systems, peptidoglycan biosynthesis and lipopolysaccharide biosynthesis. These results underscore the advantages of E. coli DH5α, which tolerates the components of transformation buffer and expresses foreign plasmids efficiently. Moreover, these SNVs were also observed in the commercially available strain. These data provide the genetic information of E. coli DH5α for its future application in metabolic engineering and synthetic biology.

Rapid selective sweep of pre-existing polymorphisms and slow fixation of new mutations in experimental evolution of Desulfovibrio vulgaris

To investigate the genetic basis of microbial evolutionary adaptation to salt (NaCl) stress, populations of Desulfovibrio vulgaris Hildenborough (DvH), a sulfate-reducing bacterium important for the biogeochemical cycling of sulfur, carbon and nitrogen, and potentially the bioremediation of toxic heavy metals and radionuclides, were propagated under salt stress or non-stress conditions for 1200 generations. Whole-genome sequencing revealed 11 mutations in salt stress-evolved clone ES9-11 and 14 mutations in non-stress-evolved clone EC3-10. Whole-population sequencing data suggested the rapid selective sweep of the pre-existing polymorphisms under salt stress within the first 100 generations and the slow fixation of new mutations. Population genotyping data demonstrated that the rapid selective sweep of pre-existing polymorphisms was common in salt stress-evolved populations. In contrast, the selection of pre-existing polymorphisms was largely random in EC populations. Consistently, at 100 generations, stress-evolved population ES9 showed improved salt tolerance, namely increased growth rate (2.0-fold), higher biomass yield (1.8-fold) and shorter lag phase (0.7-fold) under higher salinity conditions. The beneficial nature of several mutations was confirmed by site-directed mutagenesis. All four tested mutations contributed to the shortened lag phases under higher salinity condition. In particular, compared with the salt tolerance improvement in ES9-11, a mutation in a histidine kinase protein gene lytS contributed 27% of the growth rate increase and 23% of the biomass yield increase while a mutation in hypothetical gene DVU2472 contributed 24% of the biomass yield increase. Our results suggested that a few beneficial mutations could lead to dramatic improvements in salt tolerance.

Thermal and solvent stress cross-tolerance conferred to Corynebacterium glutamicum by adaptive laboratory evolution

Reinforcing microbial thermotolerance is a strategy to enable fermentation with flexible temperature settings and thereby to save cooling costs. Here, we report on adaptive laboratory evolution (ALE) of the amino acid-producing bacterium Corynebacterium glutamicum under thermal stress. After 65 days of serial passage of the transgenic strain GLY3, in which the glycolytic pathway is optimized for alanine production under oxygen deprivation, three strains adapted to supraoptimal temperatures were isolated, and all the mutations they acquired were identified by whole-genome resequencing. Of the 21 mutations common to the three strains, one large deletion and two missense mutations were found to promote growth of the parental strain under thermal stress. Additive effects on thermotolerance were observed among these mutations, and the combination of the deletion with the missense mutation on otsA, encoding a trehalose-6-phosphate synthase, allowed the parental strain to overcome the upper limit of growth temperature. Surprisingly, the three evolved strains acquired cross-tolerance for isobutanol, which turned out to be partly attributable to the genomic deletion associated with the enhanced thermotolerance. The deletion involved loss of two transgenes, pfk and pyk, encoding the glycolytic enzymes, in addition to six native genes, and elimination of the transgenes, but not the native genes, was shown to account for the positive effects on thermal and solvent stress tolerance, implying a link between energy-producing metabolism and bacterial stress tolerance. Overall, the present study provides evidence that ALE can be a powerful tool to refine the phenotype of C. glutamicum and to investigate the molecular bases of stress tolerance.

Evolution of Escherichia coli to 42 °C and subsequent genetic engineering reveals adaptive mechanisms and novel mutations

Adaptive laboratory evolution (ALE) has emerged as a valuable method by which to investigate microbial adaptation to a desired environment. Here, we performed ALE to 42 °C of ten parallel populations of Escherichia coli K-12 MG1655 grown in glucose minimal media. Tightly controlled experimental conditions allowed selection based on exponential-phase growth rate, yielding strains that uniformly converged toward a similar phenotype along distinct genetic paths. Adapted strains possessed as few as 6 and as many as 55 mutations, and of the 144 genes that mutated in total, 14 arose independently across two or more strains. This mutational recurrence pointed to the key genetic targets underlying the evolved fitness increase. Genome engineering was used to introduce the novel ALE-acquired alleles in random combinations into the ancestral strain, and competition between these engineered strains reaffirmed the impact of the key mutations on the growth rate at 42 °C. Interestingly, most of the identified key gene targets differed significantly from those found in similar temperature adaptation studies, highlighting the sensitivity of genetic evolution to experimental conditions and ancestral genotype. Additionally, transcriptomic analysis of the ancestral and evolved strains revealed a general trend for restoration of the global expression state back toward preheat stressed levels. This restorative effect was previously documented following evolution to metabolic perturbations, and thus may represent a general feature of ALE experiments. The widespread evolved expression shifts were enabled by a comparatively scant number of regulatory mutations, providing a net fitness benefit but causing suboptimal expression levels for certain genes, such as those governing flagellar formation, which then became targets for additional ameliorating mutations. Overall, the results of this study provide insight into the adaptation process and yield lessons important for the future implementation of ALE as a tool for scientific research and engineering.

Parallel genome-wide fixation of ancestral alleles in partially outcrossing experimental populations of Caenorhabditis elegans

Experimental evolution studies, coupled with new advances in DNA sequencing technology, have become a powerful tool for exploring how populations respond to selection at the genomic level. Recent experiments in microbes typically have found evidence for multiple novel mutations, which are usually fixed. In contrast, in animal model systems, evolutionary responses seem to involve more modest changes in the frequencies of pre-existing alleles, probably because these populations outcross and are usually initialized with greater levels of standing variation. In this experiment, I used whole-genome resequencing to estimate allele frequencies and look for novel substitutions in experimentally evolved populations of Caenorhabditis elegans. These populations were founded with a fixed pair of deleterious mutations introgressed into multiple wild genetic backgrounds and allowed to evolve for 50 generations with a mixed mating system. There is evidence for some recombination between ancestral haplotypes, but selective sweeps seem to have resulted in the fixation of large chromosomal segments throughout most of the genome. In addition, a few new mutations were detected. Simulations suggest that strong selection and low outcrossing rates are likely explanations for the observed outcomes, consistent with earlier work showing large fitness increases in these populations over 50 generations. These results also show clear parallels to population genetic patterns in C. elegans in nature: recent selective sweeps, high linkage disequilibrium, and low effective recombination rates. Thus, the genomic consequences of selection depend heavily on the biology of the organism in question, including its mating system and levels of genetic variation.

The dynamics of diverse segmental amplifications in populations of Saccharomyces cerevisiae adapting to strong selection

Population adaptation to strong selection can occur through the sequential or parallel accumulation of competing beneficial mutations. The dynamics, diversity, and rate of fixation of beneficial mutations within and between populations are still poorly understood. To study how the mutational landscape varies across populations during adaptation, we performed experimental evolution on seven parallel populations of Saccharomyces cerevisiae continuously cultured in limiting sulfate medium. By combining quantitative polymerase chain reaction, array comparative genomic hybridization, restriction digestion and contour-clamped homogeneous electric field gel electrophoresis, and whole-genome sequencing, we followed the trajectory of evolution to determine the identity and fate of beneficial mutations. During a period of 200 generations, the yeast populations displayed parallel evolutionary dynamics that were driven by the coexistence of independent beneficial mutations. Selective amplifications rapidly evolved under this selection pressure, in particular common inverted amplifications containing the sulfate transporter gene SUL1. Compared with single clones, detailed analysis of the populations uncovers a greater complexity whereby multiple subpopulations arise and compete despite a strong selection. The most common evolutionary adaptation to strong selection in these populations grown in sulfate limitation is determined by clonal interference, with adaptive variants both persisting and replacing one another.

Adaptive laboratory evolution -- principles and applications for biotechnology

Adaptive laboratory evolution is a frequent method in biological studies to gain insights into the basic mechanisms of molecular evolution and adaptive changes that accumulate in microbial populations during long term selection under specified growth conditions. Although regularly performed for more than 25 years, the advent of transcript and cheap next-generation sequencing technologies has resulted in many recent studies, which successfully applied this technique in order to engineer microbial cells for biotechnological applications. Adaptive laboratory evolution has some major benefits as compared with classical genetic engineering but also some inherent limitations. However, recent studies show how some of the limitations may be overcome in order to successfully incorporate adaptive laboratory evolution in microbial cell factory design. Over the last two decades important insights into nutrient and stress metabolism of relevant model species were acquired, whereas some other aspects such as niche-specific differences of non-conventional cell factories are not completely understood. Altogether the current status and its future perspectives highlight the importance and potential of adaptive laboratory evolution as approach in biotechnological engineering.

Effect of long-term laboratory propagation on Chlamydia trachomatis genome dynamics

It is assumed that bacterial strains maintained in the laboratory for long time shape their genome in a different fashion from the nature-circulating strains. Here, we analyzed the impact of long-term in vitro propagation on the genome of the obligate intracellular pathogen Chlamydia trachomatis. We fully-sequenced the genome of a historical prototype strain (L2/434/Bu) and a clinical isolate (E/CS88), before and after one-year of serial in vitro passaging (up to 3500 bacterial generations). We observed a slow adaptation of C. trachomatis to the in vitro environment, which was essentially governed by four mutations for L2/434/Bu and solely one mutation for E/CS88, corresponding to estimated mutation rates from 3.84 × 10(-10) to 1.10 × 10(-9) mutations per base pair per generation. In a speculative basis, the mutations likely conferred selective advantage as: (i) mathematical modeling showed that selective advantage is mandatory for frequency increase of a mutated clone; (ii) transversions and non-synonymous mutations were overrepresented; (iii) two non-synonymous mutations affected the genes CTL0084 and CTL0610, encoding a putative transferase and a protein likely implicated in transcription regulation respectively, which are families known to be highly prone to undergone laboratory-derived advantageous mutations in other bacteria; and (iv) the mutation for E/CS88 is located likely in the regulatory region of a virulence gene (CT115/incD) believed to play a role in subverting the host cell machinery. Nevertheless, we found no significant differences in the growth rate, plasmid load, and attachment/entry rate, between strains before and after their long-term laboratory propagation. Of note, from the mixture of clones in E/CS88 initial population, an inactivating mutation in the virulence gene CT135 evolved to 100% prevalence, unequivocally indicating that this gene is superfluous for C. trachomatis survival in vitro. Globally, C. trachomatis revealed a slow in vitro adaptation that only modestly modifies the in vivo-derived genomic evolutionary landscape.

Synchronous waves of failed soft sweeps in the laboratory: remarkably rampant clonal interference of alleles at a single locus

It has increasingly been recognized that adapting populations of microbes contain not one, but many lineages continually arising and competing at once. This process, termed "clonal interference," alters the rate and dynamics of adaptation and biases winning mutations toward those with the largest selective effect. Here we uncovered a dramatic example of clonal interference between multiple similar mutations occurring at the same locus within replicate populations of Methylobacterium extorquens AM1. Because these mutational events involved the transposition of an insertion sequence into a narrow window of a single gene, they were both readily detectable at low frequencies and could be distinguished due to differences in insertion sites. This allowed us to detect up to 17 beneficial alleles of this type coexisting in a single population. Despite conferring a large selective benefit, the majority of these alleles rose and then fell in frequency due to other lineages emerging that were more fit. By comparing allele-frequency dynamics to the trajectories of fitness gains by these populations, we estimated the fitness values of the genotypes that contained these mutations. Collectively across all populations, these alleles arose upon backgrounds with a wide range of fitness values. Within any single population, however, multiple alleles tended to rise and fall synchronously during a single wave of multiple genotypes with nearly identical fitness values. These results suggest that alleles of large benefit arose repeatedly in failed "soft sweeps" during narrow windows of adaptation due to the combined effects of epistasis and clonal interference.