Genome dynamics during experimental evolution

Reinventing the wheel? Reassessing the roles of gene flow, sorting and convergence in repeated evolution

Biologists have long been intrigued by apparently predictable and repetitive evolutionary trajectories inferred across a variety of lineages and systems. In recent years, high-throughput sequencing analyses have started to transform our understanding of such repetitive shifts. While researchers have traditionally categorized such shifts as either "convergent" or "parallel," based on relatedness of the lineages involved, emerging genomic insights provide an opportunity to better describe the actual evolutionary mechanisms at play. A synthesis of recent genomic analyses confirms that convergence is the predominant driver of repetitive evolution among species, whereas repeated sorting of standing variation is the major driver of repeated shifts within species. However, emerging data reveal numerous notable exceptions to these expectations, with recent examples of de novo mutations underpinning convergent shifts among even very closely related lineages, while repetitive sorting processes have occurred among even deeply divergent taxa, sometimes via introgression. A number of very recent analyses have found evidence for both processes occurring on different scales within taxa. We suggest that the relative importance of convergent versus sorting processes depends on the interplay between gene flow among populations, and phylogenetic relatedness of the lineages involved.

Microfluidic platform for monitoring Saccharomyces cerevisiae mutation accumulation

Mutations in DNA have large-ranging consequences, from evolution to disease. Many mechanisms contribute to mutational processes such as dysfunctions in DNA repair pathways and exogenous or endogenous mutagen exposures. Model organisms and mutation accumulation (MA) experiments are indispensable to study mutagenesis. Classical MA is, however, time consuming and laborious. To fill the need for more efficient approaches to characterize mutational profiles, we have developed an innovative microfluidic-based system that automatizes MA culturing over many generations in budding yeast. This unique experimental tool, coupled with high-throughput sequencing, reduces by one order of magnitude the time required for genome-wide measurements of mutational profiles, while also parallelizing and simplifying the cell culture. To validate our approach, we performed microfluidic MA experiments on two different genetic backgrounds, a wild-type strain and a base-excision DNA repair ung1 mutant characterized by a well-defined mutational profile. We show that the microfluidic device allows for mutation accumulation comparable to the traditional method on plate. Our approach thus paves the way to massively-parallel MA experiments with minimal human intervention that can be used to investigate mutational processes at the origin of human diseases and to identify mutagenic compounds relevant for medical and environmental research.

Universal Constraints on Protein Evolution in the Long-Term Evolution Experiment with Escherichia coli

Although it is well known that abundant proteins evolve slowly across the tree of life, there is little consensus for why this is true. Here, I report that abundant proteins evolve slowly in the hypermutator populations of Lenski's long-term evolution experiment with Escherichia coli (LTEE). Specifically, the density of all observed mutations per gene, as measured in metagenomic time series covering 60,000 generations of the LTEE, significantly anticorrelates with mRNA abundance, protein abundance, and degree of protein-protein interaction. The same pattern holds for nonsynonymous mutation density. However, synonymous mutation density, measured across the LTEE hypermutator populations, positively correlates with protein abundance. These results show that universal constraints on protein evolution are visible in data spanning three decades of experimental evolution. Therefore, it should be possible to design experiments to answer why abundant proteins evolve slowly.

Genetic Basis of Chromate Adaptation and the Role of the Pre-existing Genetic Divergence during an Experimental Evolution Study with Desulfovibrio vulgaris Populations

Hexavalent chromium [Cr(VI)] is a common environmental pollutant. However, little is known about the genetic basis of microbial evolution under Cr(VI) stress and the influence of the prior evolution histories on the subsequent evolution under Cr(VI) stress. In this study, Desulfovibrio vulgaris Hildenborough (DvH), a model sulfate-reducing bacterium, was experimentally evolved for 600 generations. By evolving the replicate populations of three genetically diverse DvH clones, including ancestor (AN, without prior experimental evolution history), non-stress-evolved EC3-10, and salt stress-evolved ES9-11, the contributions of adaptation, chance, and pre-existing genetic divergence to the evolution under Cr(VI) stress were able to be dissected. Significantly decreased lag phases under Cr(VI) stress were observed in most evolved populations, while increased Cr(VI) reduction rates were primarily observed in populations evolved from EC3-10 and ES9-11. The pre-existing genetic divergence in the starting clones showed strong influences on the changes in lag phases, growth rates, and Cr(VI) reduction rates. Additionally, the genomic mutation spectra in populations evolved from different starting clones were significantly different. A total of 14 newly mutated genes obtained mutations in at least two evolved populations, suggesting their importance in Cr(VI) adaptation. An in-frame deletion mutation of one of these genes, the chromate transporter gene DVU0426, demonstrated that it played an important role in Cr(VI) tolerance. Overall, our study identified potential key functional genes for Cr(VI) tolerance and demonstrated the important role of pre-existing genetic divergence in evolution under Cr(VI) stress conditions.

IMPORTANCE Chromium is one of the most common heavy metal pollutants of soil and groundwater. The potential of Desulfovibrio vulgaris Hildenborough in heavy metal bioremediation such as Cr(VI) reduction was reported previously; however, experimental evidence of key functional genes involved in Cr(VI) resistance are largely unknown. Given the genetic divergence of microbial populations in nature, knowledge on how this divergence affects the microbial adaptation to a new environment such as Cr(VI) stress is very limited. Taking advantage of our previous study, three groups of genetically diverse D. vulgaris Hildenborough populations with or without prior experimental evolution histories were propagated under Cr(VI) stress for 600 generations. Whole-population genome resequencing of the evolved populations revealed the genomic changes underlying the improved Cr(VI) tolerance. The strong influence of the pre-existing genetic divergence in the starting clones on evolution under Cr(VI) stress conditions was demonstrated at both phenotypic and genetic levels.

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

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

Single nucleotide polymorphisms (SNPs) in prostate cancer: its implications in diagnostics and therapeutics

Prostate cancer is one of the most frequently diagnosed malignancies in developed countries and approximately 248,530 new cases of prostate cancer are likely to be diagnosed in the United States in 2021. During the late 1990s and 2000s, the prostate cancer-related death rate has decreased by 4% per year on average because of advancements in prostate-specific antigen (PSA) testing. However, the non-specificity of PSA to distinguish between benign and malignant forms of cancer is a major concern in the management of prostate cancer. Despite other risk factors in the pathogenesis of prostate cancer, recent advancement in molecular genetics suggests that genetic heredity plays a crucial role in prostate carcinogenesis. Approximately, 60% of heritability and more than 100 well-recognized single-nucleotide-polymorphisms (SNPs) have been found to be associated with prostate cancer and constitute a major risk factor in the development of prostate cancer. Recent findings revealed that a low to moderate effect on the progression of prostate cancer of individual SNPs was observed compared to a strong progressive effect when SNPs were in combination. Here, in this review, we made an attempt to critically analyze the role of SNPs and associated genes in the development of prostate cancer and their implications in diagnostics and therapeutics. A better understanding of the role of SNPs in prostate cancer susceptibility may improve risk prediction, enhance fine-mapping, and furnish new insights into the underlying pathophysiology of prostate cancer.

Correlation between the spatial distribution and colony size was common for monogenetic bacteria in laboratory conditions

Background

Geographically separated population growth of microbes is a common phenomenon in microbial ecology. Colonies are representative of the morphological characteristics of this structured population growth. Pattern formation by single colonies has been intensively studied, whereas the spatial distribution of colonies is poorly investigated.

Results

The present study describes a first trial to address the questions of whether and how the spatial distribution of colonies determines the final colony size using the model microorganism Escherichia coli, colonies of which can be grown under well-controlled laboratory conditions. A computational tool for image processing was developed to evaluate colony density, colony size and size variation, and the Voronoi diagram was applied for spatial analysis of colonies with identical space resources. A positive correlation between the final colony size and the Voronoi area was commonly identified, independent of genomic and nutritional differences, which disturbed the colony size and size variation.

Conclusions

This novel finding of a universal correlation between the spatial distribution and colony size not only indicated the fair distribution of spatial resources for monogenetic colonies growing with identical space resources but also indicated that the initial localization of the microbial colonies decided by chance determined the fate of the subsequent population growth. This study provides a valuable example for quantitative analysis of the complex microbial ecosystems by means of experimental ecology.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-021-02180-8.

Dynamics of Adaptation During Three Years of Evolution Under Long-Term Stationary Phase

Many bacterial species that cannot sporulate, such as the model bacterium Escherichia coli, can nevertheless survive for years, following exhaustion of external resources, in a state termed long-term stationary phase (LTSP). Here we describe the dynamics of E. coli adaptation during the first three years spent under LTSP. We show that during this time, E. coli continuously adapts genetically through the accumulation of mutations. For nonmutator clones, the majority of mutations accumulated appear to be adaptive under LTSP, reflected in an extremely convergent pattern of mutation accumulation. Despite the rapid and convergent manner in which populations adapt under LTSP, they continue to harbor extensive genetic variation. The dynamics of evolution of mutation rates under LTSP are particularly interesting. The emergence of mutators affects overall mutation accumulation rates as well as the mutational spectra and the ultimate spectrum of adaptive alleles acquired under LTSP. With time, mutators can evolve even higher mutation rates through the acquisition of additional mutation rate-enhancing mutations. Different mutator and nonmutator clones within a single population and time point can display extreme variation in their mutation rates, resulting in differences in both the dynamics of adaptation and their associated deleterious burdens. Despite these differences, clones that vary greatly in their mutation rates tend to coexist within their populations for many years, under LTSP.

Daptomycin Resistance in Enterococcus faecium Can Be Delayed by Disruption of the LiaFSR Stress Response Pathway

LiaFSR signaling plays a major role in mediating daptomycin (DAP) resistance in enterococci, and the lack of a functional LiaFSR pathway leads to DAP hypersusceptibility. Using in vitro experimental evolution, we evaluated how Enterococcus faecium with a liaR response regulator gene deletion evolved DAP resistance. We found that knocking out LiaFSR signaling significantly delayed the onset of resistance, but resistance could emerge eventually through various alternate mechanisms that were influenced by the environment.

pH Adaptation Drives Diverse Phenotypes in a Beneficial Bacterium-Host Mutualism

Abiotic variation can influence the evolution of specific phenotypes that contribute to the diversity of bacterial strains observed in the natural environment. Environmentally transmitted symbiotic bacteria are particularly vulnerable to abiotic fluctuations, given that they must accommodate the transition between the free-living state and the host's internal environment. This type of life history strategy can strongly influence the success of a symbiont, and whether adapting to changes outside the host will allow a greater capacity to survive in symbiosis with the host partner. One example of how environmental breadth is advantageous to the symbiosis is the beneficial association between Vibrio fischeri and sepiolid squids (Cephalopoda: Sepiolidae). Since Vibrio bacteria are environmentally transmitted, they are subject to a wide variety of abiotic variables prior to infecting juvenile squids and must be poised to survive in the host light organ. In order to better understand how a changing abiotic factor (e.g., pH) influences the diversification of symbionts and their eventual symbiotic competence, we used an experimental evolution approach to ascertain how pH adaptation affects symbiont fitness. Results show that low pH adapted Vibrio strains have more efficient colonization rates compared to their ancestral strains. In addition, growth rates had significant differences compared to ancestral strains (pH 6.5–6.8, and 7.2). Bioluminescence production (a marker for symbiont competence) of pH evolved strains also improved at pH 6.5–7.2. Results imply that the evolution and diversification of Vibrio strains adapted to low pH outside the squid improves fitness inside the squid by allowing a higher success rate for host colonization and symbiotic competence.

Msh Pilus Mutations Increase the Ability of a Free-Living Bacterium to Colonize a Piscine Host

Symbioses between animals and bacteria are ubiquitous. To better understand these relationships, it is essential to unravel how bacteria evolve to colonize hosts. Previously, we serially passaged the free-living bacterium, Shewanella oneidensis, through the digestive tracts of germ-free larval zebrafish (Danio rerio) to uncover the evolutionary changes involved in the initiation of a novel symbiosis with a vertebrate host. After 20 passages, we discovered an adaptive missense mutation in the mshL gene of the msh pilus operon, which improved host colonization, increased swimming motility, and reduced surface adhesion. In the present study, we determined that this mutation was a loss-of-function mutation and found that it improved zebrafish colonization by augmenting S. oneidensis representation in the water column outside larvae through a reduced association with environmental surfaces. Additionally, we found that strains containing the mshL mutation were able to immigrate into host digestive tracts at higher rates per capita. However, mutant and evolved strains exhibited no evidence of a competitive advantage after colonizing hosts. Our results demonstrate that bacterial behaviors outside the host can play a dominant role in facilitating the onset of novel host associations.

Bacterial fitness landscapes stratify based on proteome allocation associated with discrete aero-types

The fitness landscape is a concept commonly used to describe evolution towards optimal phenotypes. It can be reduced to mechanistic detail using genome-scale models (GEMs) from systems biology. We use recently developed GEMs of Metabolism and protein Expression (ME-models) to study the distribution of Escherichia coli phenotypes on the rate-yield plane. We found that the measured phenotypes distribute non-uniformly to form a highly stratified fitness landscape. Systems analysis of the ME-model simulations suggest that this stratification results from discrete ATP generation strategies. Accordingly, we define "aero-types", a phenotypic trait that characterizes how a balanced proteome can achieve a given growth rate by modulating 1) the relative utilization of oxidative phosphorylation, glycolysis, and fermentation pathways; and 2) the differential employment of electron-transport-chain enzymes. This global, quantitative, and mechanistic systems biology interpretation of fitness landscape formed upon proteome allocation offers a fundamental understanding of bacterial physiology and evolution dynamics.

Transcriptional Regulatory Networks Involved in C3-C4 Alcohol Stress Response and Tolerance in Yeast

Alcohol toxicity significantly impacts the titer and productivity of industrially produced biofuels. To overcome this limitation, we must find and use strategies to improve stress tolerance in production strains. Previously, we developed a multiplex navigation of a global regulatory network (MINR) library that targeted 25 regulatory genes that are predicted to modify global regulation in yeast under different stress conditions. In this study, we expanded this concept to target the active sites of 47 transcriptional regulators using a saturation mutagenesis library. The 47 targeted regulators interact with more than half of all yeast genes. We then screened and selected for C3-C4 alcohol tolerance. We identified specific mutants that have resistance to isopropanol and isobutanol. Notably, the WAR1_K110N variant improved tolerance to both isopropanol and isobutanol. In addition, we investigated the mechanisms for improvement of isopropanol and isobutanol stress tolerance and found that genes related to glycolysis play a role in tolerance to isobutanol, while changes in ATP synthesis and mitochondrial respiration play a role in tolerance to both isobutanol and isopropanol. Overall, this work sheds light on basic mechanisms for isopropanol and isobutanol toxicity and demonstrates a promising strategy to improve tolerance to C3-C4 alcohols by perturbing the transcriptional regulatory network.

An interbacterial DNA deaminase toxin directly mutagenizes surviving target populations

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

Identification and Characterization of Base-Substitution Mutations in the Macronuclear Genome of the Ciliate Tetrahymena thermophila

Polyploidy can provide adaptive advantages and drive evolution. Amitotic division of the polyploid macronucleus (MAC) in ciliates acts as a nonsexual genetic mechanism to enhance adaptation to stress conditions and thus provides a unique model to investigate the evolutionary role of polyploidy. Mutation is the primary source of the variation responsible for evolution and adaptation; however, to date, de novo mutations that occur in ciliate MAC genomes during these processes have not been characterized and their biological impacts are undefined. Here, we carried out long-term evolution experiments to directly explore de novo MAC mutations and their molecular features in the model ciliate, Tetrahymena thermophila. A simple but effective method was established to detect base-substitution mutations in evolving populations whereas filtering out most of the false positive base-substitutions caused by repetitive sequences and the programmed genome rearrangements. The detected mutations were rigorously validated using the MassARRAY system. Validated mutations showed a strong G/C→A/T bias, consistent with observations in other species. Moreover, a progressive increase in growth rate of the evolving populations suggested that some of these mutations might be responsible for cell fitness. The established mutation identification and validation methods will be an invaluable resource to make ciliates an important model system to study the role of polyploidy in evolution.

Adaptive evolution of nontransitive fitness in yeast

A common misconception is that evolution is a linear ‘march of progress’, where each organism along a line of descent is more fit than all those that came before it. Rejecting this misconception implies that evolution is nontransitive: a series of adaptive events will, on occasion, produce organisms that are less fit compared to a distant ancestor. Here we identify a nontransitive evolutionary sequence in a 1000-generation yeast evolution experiment. We show that nontransitivity arises due to adaptation in the yeast nuclear genome combined with the stepwise deterioration of an intracellular virus, which provides an advantage over viral competitors within host cells. Extending our analysis, we find that nearly half of our ~140 populations experience multilevel selection, fixing adaptive mutations in both the nuclear and viral genomes. Our results provide a mechanistic case-study for the adaptive evolution of nontransitivity due to multilevel selection in a 1000-generation host/virus evolution experiment.

It is widely accepted in biology that all life on Earth gradually evolved over billions of years from a single ancestor. Yet, there is still much about this process that is not fully understood. Evolution is often thought of as progressing in a linear fashion, with each new generation being better adapted to its environment than the last. But it has been proposed that evolution is also nontransitive: this means even if each generation is ‘fitter’ than its immediate predecessor, these series of adaptive changes will occasionally result in organisms that are less fit than their distant ancestors.

Laboratory experiments of evolution are a good way to test evolutionary theories because they allow researchers to create scenarios that are impossible to observe in natural populations, such as an organism competing against its extinct ancestors. Buskirk et al. set up such an experiment using yeast to determine whether nontransitive effects can be observed in the direct descendants of an organism.

At the start of the experiment, the yeast cells were host to a non-infectious ‘killer’ virus that is common among yeast. Cells containing the virus produce a toxin that destroys other yeast that lack the virus. The populations of yeast were given a nutrient-rich broth in which to grow and subjected to a simple evolutionary pressure: to grow fast, which limits the amount of resources available.

As the yeast evolved, they gained beneficial genetic mutations that allowed them to outcompete their neighbors, and they passed these traits down to their descendants. Some of these mutations occurred not in the yeast genome, but in the genome of the killer virus, and this stopped the yeast infected with the virus from producing the killer toxin. Over time, other mutations resulted in the infected yeast no longer being immune to the toxin. Thus, when Buskirk et al. pitted these yeast against their distant ancestors, the new generation were destroyed by the toxins the older generation produced.

These findings provide the first experimental evidence for nontransitivity along a line of descent. The results have broad implications for our understanding of how evolution works, casting doubts over the idea that evolution always involves a direct progression towards new, improved traits.

Reflections on the Predictability of Evolution: Toward a Conceptual Framework

Evolution is generally considered to be unpredictable because genetic variations are known to occur randomly. However, remarkable patterns of repeated convergent evolution are observed, for instance, loss of pigments by organisms living in caves. Analogous phenotypes appear in similar environments, sometimes in response to similar constraints. Alongside randomness, a certain evolutionary determinism also exists, for instance, the selection of particular phenotypes subjected to particular environmental constraints in the “evolutionary funnel.” We pursue the idea that eco-evolutionary specialization is in some way determinist. The conceptual framework of phenotypic changes entailing specialization presented in this essay explains how evolution can be predicted. We also discuss how the predictability of evolution could be tested using the case of metabolic specialization through gene losses. We also put forward that microorganisms could be key models to test and possibly make headway evolutionary predictions and knowledge about evolution.

Evidence for the Accumulation of Nonsynonymous Mutations and Favorable Pleiotropic Alleles During Wheat Breeding

Plant breeding leads to the genetic improvement of target traits by selecting a small number of genotypes from among typically large numbers of candidate genotypes after careful evaluation. In this study, we first investigated how mutations at conserved nucleotide sites normally viewed as deleterious, such as nonsynonymous sites, accumulated in a wheat, Triticum aestivum, breeding lineage. By comparing a 150 year old ancestral and modern cultivar, we found recent nucleotide polymorphisms altered amino acids and occurred within conserved genes at frequencies expected in the absence of purifying selection. Mutations that are deleterious in other contexts likely had very small or no effects on target traits within the breeding lineage. Second, we investigated if breeders selected alleles with favorable effects on some traits and unfavorable effects on others and used different alleles to compensate for the latter. An analysis of a segregating population derived from the ancestral and modern parents provided one example of this phenomenon. The recent cultivar contains the Rht-B1b green revolution semi-dwarfing allele and compensatory alleles that reduce its negative effects. However, improvements in traits other than plant height were due to pleiotropic loci with favorable effects on traits and to favorable loci with no detectable pleiotropic effects. Wheat breeding appears to tolerate mutations at conserved nucleotide sites and to only select for alleles with both favorable and unfavorable effects on traits in exceptional situations.

Rampant loss of social traits during domestication of a Bacillus subtilis natural isolate

Most model bacteria have been domesticated in laboratory conditions. Yet, the tempo with which a natural isolate diverges from its ancestral phenotype under domestication to a novel laboratory environment is poorly understood. Such knowledge, however is essential to understanding the rate of evolution, the time scale over which a natural isolate can be propagated without loss of its natural adaptive traits, and the reliability of experimental results across labs. Using experimental evolution, phenotypic assays, and whole-genome sequencing, we show that within a week of propagation in a common laboratory environment, a natural isolate of Bacillus subtilis acquires mutations that cause changes in a multitude of traits. A single adaptive mutational step in the gene coding for the transcriptional regulator DegU impairs a DegU-dependent positive autoregulatory loop and leads to loss of robust biofilm architecture, impaired swarming motility, reduced secretion of exoproteases, and to changes in the dynamics of sporulation across environments. Importantly, domestication also resulted in improved survival when the bacteria face pressure from cells of the innate immune system. These results show that degU is a target for mutations during domestication and underscores the importance of performing careful and extremely short-term propagations of natural isolates to conserve the traits encoded in their original genomes.

Increased time sampling in an evolve-and-resequence experiment with outcrossing Saccharomyces cerevisiae reveals multiple paths of adaptive change

"Evolve and resequence" (E&R) studies combine experimental evolution and whole-genome sequencing to interrogate the genetics underlying adaptation. Due to ease of handling, E&R work with asexual organisms such as bacteria can employ optimized experimental design, with large experiments and many generations of selection. By contrast, E&R experiments with sexually reproducing organisms are more difficult to implement, and design parameters vary dramatically among studies. Thus, efforts have been made to assess how these differences, such as number of independent replicates, or size of experimental populations, impact inference. We add to this work by investigating the role of time sampling-the number of discrete time points sequence data are collected from evolving populations. Using data from an E&R experiment with outcrossing Saccharomyces cerevisiae in which populations were sequenced 17 times over ~540 generations, we address the following questions: (a) Do more time points improve the ability to identify candidate regions underlying selection? And (b) does high-resolution sampling provide unique insight into evolutionary processes driving adaptation? We find that while time sampling does not improve the ability to identify candidate regions, high-resolution sampling does provide valuable opportunities to characterize evolutionary dynamics. Increased time sampling reveals three distinct trajectories for adaptive alleles: one consistent with classic population genetic theory (i.e., models assuming constant selection coefficients), and two where trajectories suggest more context-dependent responses (i.e., models involving dynamic selection coefficients). We conclude that while time sampling has limited impact on candidate region identification, sampling eight or more time points has clear benefits for studying complex evolutionary dynamics.

Earlier onset of flowering and increased reproductive allocation of an annual invasive plant in the north of its novel range

BACKGROUND AND AIMS

It remains unclear whether invasive species can maintain both high biomass and reproductive output across their invaded range. Along latitudinal gradients, allocation theory predicts that faster flowering onset at high latitudes results in maturation at smaller size and thus reduced reproductive output. For annual invasive plants, more favourable environmental conditions at low latitudes probably result in stronger competition of co-occurring species, potentially driving selection for higher investment in vegetative biomass, while harsher climatic conditions and associated reproductive uncertainty at higher latitudes could reduce selection for vegetative biomass and increased selection for high reproductive investment (stress-gradient hypothesis). Combined, these drivers could result in increased or constant reproductive allocation with increasing latitude.

METHODS

We quantified life-history traits in the invasive annual plant Impatiens glandulifera along a latitudinal gradient in Europe. By growing two successive glasshouse generations, we assessed genetic differentiation in vegetative growth and reproductive output across six populations, and tested whether onset of flowering drives this divergence.

KEY RESULTS

Trait variation was mainly caused by genetic differentiation. As expected, flowering onset was progressively earlier in populations from higher latitudes. Plant height and vegetative biomass also decreased in populations from higher latitudes, as predicted by allocation theory, but their variation was independent of the variation in flowering onset. Reproductive output remained constant across latitudes, resulting in increased reproductive allocation towards higher latitudes, supporting the stress-gradient hypothesis. We also observed trait genetic differentiation among populations that was independent of latitude.

CONCLUSIONS

We show that an annual invasive plant evolved several life-history traits across its invaded range in ~150 years. The evolution of vegetative and reproductive traits seems unconstrained by evolution of flowering onset. This genetic decoupling between vegetative and reproductive traits possibly contributes to the invasion success of this species.

Molecular epidemiology of leprosy: An update

Molecular epidemiology investigations are notoriously challenging in the leprosy field mainly because the inherent characteristics of the disease as well as its yet uncultivated causative agents, Mycobacterium leprae and M. lepromatosis. Despite significant developments in understanding the biology of leprosy bacilli through genomic approaches, the exact mechanisms of transmission is still unclear and the factors underlying pathological variation of the disease in different patients remain as major gaps in our knowledge about leprosy. Despite these difficulties, the last two decades have seen the development of genotyping procedures based on PCR-sequencing of target loci as well as by the genome-wide analysis of an increasing number of geographically diverse isolates of leprosy bacilli. This has provided a foundation for molecular epidemiology studies that are bringing a better understanding of strain evolution associated with ancient human migrations, and phylogeographical insights about the spread of disease globally. This review discusses the advantages and drawbacks of the main tools available for molecular epidemiological investigations of leprosy and summarizes various methods ranging from PCR-based genotyping to genome-typing techniques. We also describe their main applications in analyzing the short-range and long-range transmission of the disease. Finally, we summarise the current gaps and challenges that remain in the field of molecular epidemiology of leprosy.

Experimental Evolution of Bacillus subtilis Reveals the Evolutionary Dynamics of Horizontal Gene Transfer and Suggests Adaptive and Neutral Effects

Tracing evolutionary processes that lead to fixation of genomic variation in wild bacterial populations is a prime challenge in molecular evolution. In particular, the relative contribution of horizontal gene transfer (HGT) vs.de novo mutations during adaptation to a new environment is poorly understood. To gain a better understanding of the dynamics of HGT and its effect on adaptation, we subjected several populations of competent Bacillus subtilis to a serial dilution evolution on a high-salt-containing medium, either with or without foreign DNA from diverse pre-adapted or naturally salt tolerant species. Following 504 generations of evolution, all populations improved growth yield on the medium. Sequencing of evolved populations revealed extensive acquisition of foreign DNA from close Bacillus donors but not from more remote donors. HGT occurred in bursts, whereby a single bacterial cell appears to have acquired dozens of fragments at once. In the largest burst, close to 2% of the genome has been replaced by HGT. Acquired segments tend to be clustered in integration hotspots. Other than HGT, genomes also acquired spontaneous mutations. Many of these mutations occurred within, and seem to alter, the sequence of flagellar proteins. Finally, we show that, while some HGT fragments could be neutral, others are adaptive and accelerate evolution.

Experimental evolution reveals nitrate tolerance mechanisms in Desulfovibrio vulgaris

Elevated nitrate in the environment inhibits sulfate reduction by important microorganisms of sulfate-reducing bacteria (SRB). Several SRB may respire nitrate to survive under elevated nitrate, but how SRB that lack nitrate reductase survive to elevated nitrate remains elusive. To understand nitrate adaptation mechanisms, we evolved 12 populations of a model SRB (i.e., Desulfovibrio vulgaris Hildenborough, DvH) under elevated NaNO₃ for 1000 generations, analyzed growth and acquired mutations, and linked their genotypes with phenotypes. Nitrate-evolved (EN) populations significantly (p < 0.05) increased nitrate tolerance, and whole-genome resequencing identified 119 new mutations in 44 genes of 12 EN populations, among which six functional gene groups were discovered with high mutation frequencies at the population level. We observed a high frequency of nonsense or frameshift mutations in nitrosative stress response genes (NSR: DVU2543, DVU2547, and DVU2548), nitrogen regulatory protein C family genes (NRC: DVU2394-2396, DVU2402, and DVU2405), and nitrate cluster (DVU0246-0249 and DVU0251). Mutagenesis analysis confirmed that loss-of-functions of NRC and NSR increased nitrate tolerance. Also, functional gene groups involved in fatty acid synthesis, iron regulation, and two-component system (LytR/LytS) known to be responsive to multiple stresses, had a high frequency of missense mutations. Mutations in those gene groups could increase nitrate tolerance through regulating energy metabolism, barring entry of nitrate into cells, altering cell membrane characteristics, or conferring growth advantages at the stationary phase. This study advances our understanding of nitrate tolerance mechanisms and has important implications for linking genotypes with phenotypes in DvH.

In vitro selection predicts malaria parasite resistance to dihydroorotate dehydrogenase inhibitors in a mouse infection model

Resistance has developed in Plasmodium malaria parasites to every antimalarial drug in clinical use, prompting the need to characterize the pathways mediating resistance. Here, we report a framework for assessing development of resistance of Plasmodium falciparum to new antimalarial therapeutics. We investigated development of resistance by P. falciparum to the dihydroorotate dehydrogenase (DHODH) inhibitors DSM265 and DSM267 in tissue culture and in a mouse model of P. falciparum infection. We found that resistance to these drugs arose rapidly both in vitro and in vivo. We identified 13 point mutations mediating resistance in the parasite DHODH in vitro that overlapped with the DHODH mutations that arose in the mouse infection model. Mutations in DHODH conferred increased resistance (ranging from 2- to ~400-fold) to DHODH inhibitors in P. falciparum in vitro and in vivo. We further demonstrated that the drug-resistant parasites carrying the C276Y mutation had mitochondrial energetics comparable to the wild-type parasite and also retained their fitness in competitive growth experiments. Our data suggest that in vitro selection of drug-resistant P. falciparum can predict development of resistance in a mouse model of malaria infection.

Host-parasite co-evolution and its genomic signature

Studies in diverse biological systems have indicated that host-parasite co-evolution is responsible for the extraordinary genetic diversity seen in some genomic regions, such as major histocompatibility (MHC) genes in jawed vertebrates and resistance genes in plants. This diversity is believed to evolve under balancing selection on hosts by parasites. However, the mechanisms that link the genomic signatures in these regions to the underlying co-evolutionary process are only slowly emerging. We still lack a clear picture of the co-evolutionary concepts and of the genetic basis of the co-evolving phenotypic traits in the interacting antagonists. Emerging genomic tools that provide new options for identifying underlying genes will contribute to a fuller understanding of the co-evolutionary process.

Engineering regulatory networks for complex phenotypes in E. coli

Regulatory networks describe the hierarchical relationship between transcription factors, associated proteins, and their target genes. Regulatory networks respond to environmental and genetic perturbations by reprogramming cellular metabolism. Here we design, construct, and map a comprehensive regulatory network library containing 110,120 specific mutations in 82 regulators expected to perturb metabolism. We screen the library for different targeted phenotypes, and identify mutants that confer strong resistance to various inhibitors, and/or enhanced production of target compounds. These improvements are identified in a single round of selection, showing that the regulatory network library is universally applicable and is convenient and effective for engineering targeted phenotypes. The facile construction and mapping of the regulatory network library provides a path for developing a more detailed understanding of global regulation in E. coli, with potential for adaptation and use in less-understood organisms, expanding toolkits for future strain engineering, synthetic biology, and broader efforts.

Adaptations Accumulated under Prolonged Resource Exhaustion Are Highly Transient

Many nonsporulating bacterial species can survive for years within exhausted growth media in a state termed long-term stationary phase (LTSP). We have been carrying out evolutionary experiments aimed at elucidating the dynamics of genetic adaptation under LTSP. We showed that Escherichia coli adapts to prolonged resource exhaustion through the highly convergent acquisition of mutations. In the most striking example of such convergent adaptation, we observed that across all independently evolving LTSP populations, over 90% of E. coli cells carry mutations to one of three specific sites of the RNA polymerase core enzyme (RNAPC). These LTSP adaptations reduce the ability of the cells carrying them to grow once fresh resources are again provided. Here, we examine how LTSP populations recover from costs associated with their adaptation once resources are again provided to them. We demonstrate that due to the ability of LTSP populations to maintain high levels of standing genetic variation during adaptation, costly adaptations are very rapidly purged from the population once they are provided with fresh resources. We further demonstrate that recovery from costs acquired during adaptation under LTSP occurs more rapidly than would be possible if LTSP adaptations had fixed during the time populations spent under resource exhaustion. Finally, we previously reported that under LTSP, some clones develop a mutator phenotype, greatly increasing their mutation accumulation rates. Here, we show that the mechanisms by which populations recover from costs associated with fixed adaptations may depend on mutator status.IMPORTANCE Many bacterial species can survive for decades under starvation, following the exhaustion of external growth resources. We have previously shown that bacteria genetically adapt under these conditions in a manner that reduces their ability to grow once resources again become available. Here, we study how populations that have been subject to very prolonged resource exhaustion recover from costs associated with their adaptation. We demonstrate that rapid adaptations acquired under prolonged starvation tend to be highly transient, rapidly reducing in frequency once bacteria are no longer starved. Our results shed light on the longer-term consequences of bacterial survival under prolonged starvation. More generally, these results may also be applicable to understanding longer-term consequences of rapid adaptation to additional conditions as well.

The omnistat: A flexible continuous-culture system for prolonged experimental evolution

Microbial evolution experiments provide a powerful tool to unravel the molecular basis of adaptive evolution but their outcomes can be difficult to interpret, unless the selective forces that are applied during the experiment are carefully controlled. In this respect, experimental evolution in continuous cultures provides advantages over commonly used sequential batch-culture protocols because continuous cultures allow for more accurate control over the induced selective environment. However, commercial continuous-culture systems are large and expensive, while available DIY continuous-culture systems are not versatile enough to allow for multiple sensors and rigorous stirring.We present a modular continuous-culture system that adopts the commonly used GL45 glass laboratory bottle as a bioreactor vessel. Our design offers three advantages: first, it is equipped with a large head plate, fitting two sensors and seven input/output ports, enabling the customization of the system for many running modes (chemostat, auxostat, etc.). Second, the bioreactor is small (25-250 ml), which makes it feasible to run many replicates in parallel. Third, bioreactor modules can be coupled by uni- or bi-directional flows to induce spatiotemporal variation in selection. These features result in a particularly flexible culturing platform that facilitates the investigation of a broad range of evolutionary and ecological questions.To illustrate the versatility of our culturing system, we outline two evolution experiments that impose a temporally or spatially variable regime of selection. The first experiment illustrates how controlled temporal variation in resource availability can be utilized to select for anticipatory switching. The second experiment illustrates a spatially structured morbidostat setup that is designed to probe epistatic interactions between adaptive mutations. Furthermore, we demonstrate how sensor data can be used to stabilize selection pressures or track evolutionary adaptation.Evolution experiments in which populations are exposed to controlled spatiotemporal variation, are essential to gain insight into the process of adaptation and the mechanisms that constrain evolution. Continuous-culture systems, like the one presented here, offer control over key environmental parameters and establish a well-defined regime of selection. As such, they create the opportunity to expose evolutionary constraints in the form of phenotypic trade-offs, contributing to a mechanistic understanding of adaptive evolution.

Progress in Neoantigen Targeted Cancer Immunotherapies

Immunotherapies that harness the immune system to kill cancer cells have showed significant therapeutic efficacy in many human malignancies. A growing number of studies have highlighted the relevance of neoantigens in recognizing cancer cells by intrinsic T cells. Cancer neoantigens are a direct consequence of somatic mutations presenting on the surface of individual cancer cells. Neoantigens are fully cancer-specific and exempt from central tolerance. In addition, neoantigens are important targets for checkpoint blockade therapy. Recently, technological innovations have made neoantigen discovery possible in a variety of malignancies, thus providing an impetus to develop novel immunotherapies that selectively enhance T cell reactivity for the destruction of cancer cells while leaving normal tissues unharmed. In this review, we aim to introduce the methods of the identification of neoantigens, the mutational patterns of human cancers, related clinical trials, neoantigen burden and sensitivity to immune checkpoint blockade. Moreover, we focus on relevant challenges of targeting neoantigens for cancer treatment.

Repeated horizontal gene transfers triggered parallel evolution of magnetotaxis in two evolutionary divergent lineages of magnetotactic bacteria

Under the same selection pressures, two genetically divergent populations may evolve in parallel toward the same adaptive solutions. Here, we hypothesized that magnetotaxis (i.e., magnetically guided chemotaxis) represents a key adaptation to micro-oxic habitats in aquatic sediments and that its parallel evolution homogenized the phenotypes of two evolutionary divergent clusters of freshwater spirilla. All magnetotactic bacteria affiliated to the Magnetospirillum genus (Alphaproteobacteria class) biomineralize the same magnetic particle chains and share highly similar physiological and ultrastructural features. We looked for the processes that could have contributed at shaping such an evolutionary pattern by reconciling species and gene trees using newly sequenced genomes of Magnetospirillum related bacteria. We showed that repeated horizontal gene transfers and homologous recombination of entire operons contributed to the parallel evolution of magnetotaxis. We propose that such processes could represent a more parsimonious and rapid solution for adaptation compared with independent and repeated de novo mutations, especially in the case of traits as complex as magnetotaxis involving tens of interacting proteins. Besides strengthening the idea about the importance of such a function in micro-oxic habitats, these results reinforce previous observations in experimental evolution suggesting that gene flow could alleviate clonal interference and speed up adaptation under some circumstances.

Myriocin-induced adaptive laboratory evolution of an industrial strain of Saccharomyces cerevisiae reveals its potential to remodel lipid composition and heat tolerance

The modification of lipid composition allows cells to adjust membrane biophysical properties in response to changes in environmental temperature. Here, we use adaptive laboratory evolution (ALE) in the presence of myriocin, a sphingolipid (SLs) biosynthesis inhibitor, to remodel the lipid profile of an industrial yeast strain (LH) of Saccharomyces cerevisiae. The approach enabled to obtain a heterogeneous population (LHev) of myriocin-tolerant evolved clones characterized by its growth capacity at high temperature. Myriocin exposure also caused tolerance to soraphen A, an inhibitor of the acetyl-CoA carboxylase Acc1, the rate-limiting enzyme in fatty acid de novo production, supporting a change in lipid metabolism during ALE. In line with this, characterization of two randomly selected clones, LH03 and LH09, showed the presence of lipids with increased saturation degree and reduced acyl length. In addition, the clone LH03, which displays the greater improvement in fitness at 40°C, exhibited higher SL content as compared with the parental strain. Analysis of the LH03 and LH09 genomes revealed a loss of chromosomes affecting genes that have a role in fatty acid synthesis and elongation. The link between ploidy level and growth at high temperature was further supported by the analysis of a fully isogenic set of yeast strains with ploidy between 1N and 4N which showed that the loss of genome content provides heat tolerance. Consistent with this, a thermotolerant evolved population (LH40°) generated from the parental LH strain by heat-driven ALE exhibited a reduction in the chromosome copy number. Thus, our results identify myriocin-driven evolution as a powerful approach to investigate the mechanisms of acquired thermotolerance and to generate improved strains.

Construction of a Robust Sphingomonas sp. Strain for Welan Gum Production via the Expression of Global Transcriptional Regulator IrrE

Welan gum is a widely used microbial polysaccharide produced by Sphingomonas sp. However, an important factor hindering the expansion of its production is the maladaptation of strain to fermentation conditions. In this work, the global transcriptional regulator gene irre was selected as a stress-resistant element. And it was integrated into the site of the genomic carotene synthesis key enzyme gene crtB to construct a robust carotenoid-free welan gum producing strain. Fermentation with the recombinant strain effectively reduced the ethanol consumption and pigment content in the product. The tolerance temperature increased by 10°C without the need for controlling the pH. Under this fermentation condition, welan gum concentration could still reach 20.26 ± 0.25 g/L, which was 187.38% higher than that of the wild-type strain (7.05 ± 0.15 g/L). Transcriptome analysis showed that with the control of IrrE, more than 1000 genes that are involved in multiple pathways, including two-component system, bacterial chemotaxis, flagellar assembly, and cell cycle, exhibited changes at the transcriptional level and jointly allowed the strain to protect against environmental stresses.

Microevolutionary Dynamics of Chicken Genomes under Divergent Selection for Adiposity

Decades of artificial selection have significantly improved performance and efficiency of animal production systems. However, little is known about the microevolution of genomes due to intensive breeding. Using whole-genome sequencing, we document dynamic changes of chicken genomes under divergent selection on adiposity over 19 generations. Directional selection reduced within-line but increased between-line genomic differences. We observed that artificial selection tended to result in recruitment of preexisting variations of genes related to adipose tissue growth. In addition, novel mutations contributed to divergence of phenotypes under selection but contributed significantly less than preexisting genomic variants. Integration of 15 generations genome sequencing, genome-wide association study, and multi-omics data further identified that genes involved in signaling pathways important to adipogenesis, such as autophagy and lysosome (URI1, MBL2), neural system (CHAT), and endocrine (PCSK1) pathways, were under strong selection. Our study provides insights into the microevolutionary dynamics of domestic animal genomes under artificial selection.

Multifunctional Supramolecular Assemblies with Aggregation-Induced Emission (AIE) for Cell Line Identification, Cell Contamination Evaluation, and Cancer Cell Discrimination

It is undoubted the important role of cells in biology and medicine, but worldwide misidentified and cross-contaminated cell lines have caused much trouble in related fields. Herein, three kinds of supramolecular AIE (aggregation-induced emission) nanoassemblies were constructed by the host-guest interaction between tetraphenylethene (TPE) derivatives and cucurbit[8]uril (CB[8]). Based on the recognized mechanism of AIE, the TPE derivatives could achieve stronger fluorescence emission and higher fluorescence quantum yield after assembling with CB[8]. Moreover, the constructed supramolecular AIE complexes obtained well-confirmed nanostructures and exhibited different sizes and shapes. Correspondingly, they generated characteristic biological properties and fluorescence enhancement of cells. Inspired by the concept of Big Data Analysis, these fluorescence signals were further transformed into a unique fingerprint of cells via linear discriminant analysis. Immediately, we realized the veracious identification between a normal cell line, two cancer cell lines, and two metastasized cancer cell lines in a qualitative analysis. More importantly, it was well used to monitor the evaluation of cross-contaminated cells and the discrimination of cancer cells. As a proper bioapplication of ideal supramolecular nanomaterials, this system was easy to learn and apply, and the whole procedure was kept to 20 min, without cell disruption, centrifugation, or washing steps.

How Do Cells Adapt? Stories Told in Landscapes

Cells adapt to changing environments. Perturb a cell and it returns to a point of homeostasis. Perturb a population and it evolves toward a fitness peak. We review quantitative models of the forces of adaptation and their visualizations on landscapes. While some adaptations result from single mutations or few-gene effects, others are more cooperative, more delocalized in the genome, and more universal and physical. For example, homeostasis and evolution depend on protein folding and aggregation, energy and protein production, protein diffusion, molecular motor speeds and efficiencies, and protein expression levels. Models provide a way to learn about the fitness of cells and cell populations by making and testing hypotheses.

Genomic and phenotypic evolution of Escherichia coli in a novel citrate-only resource environment

Evolutionary innovations allow populations to colonize new ecological niches. We previously reported that aerobic growth on citrate (Cit⁺) evolved in an Escherichia coli population during adaptation to a minimal glucose medium containing citrate (DM25). Cit⁺ variants can also grow in citrate-only medium (DM0), a novel environment for E. coli. To study adaptation to this niche, we founded two sets of Cit⁺ populations and evolved them for 2500 generations in DM0 or DM25. The evolved lineages acquired numerous parallel mutations, many mediated by transposable elements. Several also evolved amplifications of regions containing the maeA gene. Unexpectedly, some evolved populations and clones show apparent declines in fitness. We also found evidence of substantial cell death in Cit⁺ clones. Our results thus demonstrate rapid trait refinement and adaptation to the new citrate niche, while also suggesting a recalcitrant mismatch between E. coli physiology and growth on citrate.

Evolution of Microbial Growth Traits Under Serial Dilution

Selection of mutants in a microbial population depends on multiple cellular traits. In serial-dilution evolution experiments, three key traits are the lag time when transitioning from starvation to growth, the exponential growth rate, and the yield (number of cells per unit resource). Here, we investigate how these traits evolve in laboratory evolution experiments using a minimal model of population dynamics, where the only interaction between cells is competition for a single limiting resource. We find that the fixation probability of a beneficial mutation depends on a linear combination of its growth rate and lag time relative to its immediate ancestor, even under clonal interference. The relative selective pressure on growth rate and lag time is set by the dilution factor; a larger dilution factor favors the adaptation of growth rate over the adaptation of lag time. The model shows that yield, however, is under no direct selection. We also show how the adaptation speeds of growth and lag depend on experimental parameters and the underlying supply of mutations. Finally, we investigate the evolution of covariation between these traits across populations, which reveals that the population growth rate and lag time can evolve a nonzero correlation even if mutations have uncorrelated effects on the two traits. Altogether these results provide useful guidance to future experiments on microbial evolution.

How to better predict long-term benefits and risks in weed biocontrol: an evolutionary perspective

Classical biological control (also called importation biological control) of weeds has a remarkable track record for efficiency and safety, but further improvement is still needed, particularly to account for potential evolutionary changes after release. Here, we discuss the increasing yet limited evidence of post-introduction evolution and describe approaches to predict evolutionary change. Recent advances include using experimental evolution studies over several generations that combine -omics tools with behavioral bioassays. This novel approach in weed biocontrol is well suited to explore the potential for rapid evolutionary change in real-time and thus can be used to estimate more accurately potential benefits and risks of agents before their importation. We outline this approach with a chrysomelid beetle used to control invasive common ragweed.

Localized Hypermutation is the Major Driver of Meningococcal Genetic Variability during Persistent Asymptomatic Carriage

Host persistence of bacteria is facilitated by mutational and recombinatorial processes that counteract loss of genetic variation during transmission and selection from evolving host responses. Genetic variation was investigated during persistent asymptomatic carriage of Neisseria meningitidis Interrogation of whole-genome sequences for paired isolates from 25 carriers showed that de novo mutations were infrequent, while horizontal gene transfer occurred in 16% of carriers. Examination of multiple isolates per time point enabled separation of sporadic and transient allelic variation from directional variation. A comprehensive comparative analysis of directional allelic variation with hypermutation of simple sequence repeats and hyperrecombination of class 1 type IV pilus genes detected an average of seven events per carrier and 2:1 bias for changes due to localized hypermutation. Directional genetic variation was focused on the outer membrane with 69% of events occurring in genes encoding enzymatic modifiers of surface structures or outer membrane proteins. Multiple carriers exhibited directional and opposed switching of allelic variants of the surface-located Opa proteins that enables continuous expression of these adhesins alongside antigenic variation. A trend for switching from PilC1 to PilC2 expression was detected, indicating selection for specific alterations in the activities of the type IV pilus, whereas phase variation of restriction modification (RM) systems, as well as associated phasevarions, was infrequent. We conclude that asymptomatic meningococcal carriage on mucosal surfaces is facilitated by frequent localized hypermutation and horizontal gene transfer affecting genes encoding surface modifiers such that optimization of adhesive functions occurs alongside escape of immune responses by antigenic variation.IMPORTANCE Many bacterial pathogens coexist with host organisms, rarely causing disease while adapting to host responses. Neisseria meningitidis, a major cause of meningitis and septicemia, is a frequent persistent colonizer of asymptomatic teenagers/young adults. To assess how genetic variation contributes to host persistence, whole-genome sequencing and hypermutable sequence analyses were performed on multiple isolates obtained from students naturally colonized with meningococci. High frequencies of gene transfer were observed, occurring in 16% of carriers and affecting 51% of all nonhypermutable variable genes. Comparative analyses showed that hypermutable sequences were the major mechanism of variation, causing 2-fold more changes in gene function than other mechanisms. Genetic variation was focused on genes affecting the outer membrane, with directional changes in proteins responsible for bacterial adhesion to host surfaces. This comprehensive examination of genetic plasticity in individual hosts provides a significant new platform for rationale design of approaches to prevent the spread of this pathogen.

Dissection of the mutation accumulation process during bacterial range expansions

BACKGROUND

Recent experimental work has shown that the evolutionary dynamics of bacteria expanding across space can differ dramatically from what we expect under well-mixed conditions. During spatial expansion, deleterious mutations can accumulate due to inefficient selection on the expansion front, potentially interfering with and modifying adaptive evolutionary processes.

RESULTS

We used whole genome sequencing to follow the genomic evolution of 10 mutator Escherichia coli lines during 39 days ( ~ 1650 generations) of a spatial expansion, which allowed us to gain a temporal perspective on the interaction of adaptive and non-adaptive evolutionary processes during range expansions. We used elastic net regression to infer the positive or negative effects of mutations on colony growth. The colony size, measured after three day of growth, decreased at the end of the experiment in all 10 lines, and mutations accumulated at a nearly constant rate over the whole experiment. We find evidence that beneficial mutations accumulate primarily at an early stage of the experiment, leading to a non-linear change of colony size over time. Indeed, the rate of colony size expansion remains almost constant at the beginning of the experiment and then decreases after ~ 12 days of evolution. We also find that beneficial mutations are enriched in genes encoding transport proteins, and genes coding for the membrane structure, whereas deleterious mutations show no enrichment for any biological process.

CONCLUSIONS

Our experiment shows that beneficial mutations target specific biological functions mostly involved in inter or extra membrane processes, whereas deleterious mutations are randomly distributed over the whole genome. It thus appears that the interaction between genetic drift and the availability or depletion of beneficial mutations determines the change in fitness of bacterial populations during range expansion.

Understanding the evolution of interspecies interactions in microbial communities

Microbial communities are complex multi-species assemblages that are characterized by a multitude of interspecies interactions, which can range from mutualism to competition. The overall sign and strength of interspecies interactions have important consequences for emergent community-level properties such as productivity and stability. It is not well understood how interspecies interactions change over evolutionary timescales. Here, we review the empirical evidence that evolution is an important driver of microbial community properties and dynamics on timescales that have traditionally been regarded as purely ecological. Next, we briefly discuss different modelling approaches to study evolution of communities, emphasizing the similarities and differences between evolutionary and ecological perspectives. We then propose a simple conceptual model for the evolution of interspecies interactions in communities. Specifically, we propose that to understand the evolution of interspecies interactions, it is important to distinguish between direct and indirect fitness effects of a mutation. We predict that in well-mixed environments, traits will be selected exclusively for their direct fitness effects, while in spatially structured environments, traits may also be selected for their indirect fitness effects. Selection of indirectly beneficial traits should result in an increase in interaction strength over time, while selection of directly beneficial traits should not have such a systematic effect. We tested our intuitions using a simple quantitative model and found support for our hypotheses. The next step will be to test these hypotheses experimentally and provide input for a more refined version of the model in turn, thus closing the scientific cycle of models and experiments. This article is part of the theme issue 'Conceptual challenges in microbial community ecology'.

Transcription factor Hap5 induces gsh2 expression to enhance 2-phenylethanol tolerance and production in an industrial yeast Candida glycerinogenes

2-Phenylethanol (2-PE) is an important flavor compound but also impairs cell growth severely, which in turn blocks its bioproduction. However, the molecular mechanism of 2-PE tolerance is unclear. In this study, a superb 2-PE stress-tolerant and producing yeast, Candida glycerinogenes, was selected to uncover the underlying mechanism of 2-PE tolerance. We discovered that Hap5 is an essential regulator to 2-PE resistance, and its induction by 2-PE stress occurs at the post-transcriptional level, rather than at the transcriptional level. Under 2-PE stress, Hap5 is activated and imported into the nucleus rapidly. Then, the nuclear Hap5 binds to the glutathione synthetase (gsh2) promoter via CCAAT box, to induce the expression of gsh2 gene. The increased gsh2 expression contributes to enhanced cellular glutathione content, and consequently alleviates ROS accumulation, lipid peroxidation, and cell membrane damage caused by 2-PE toxicity. Specifically, increasing the expression of gsh2 is effective in improving not just 2-PE tolerance (33.7% higher biomass under 29 mM 2-PE), but also 2-PE production (16.2% higher). This study extends our knowledge of 2-PE tolerance mechanism and also provides a promising strategy to improve 2-PE production.

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

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

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

Chance and necessity in the pleiotropic consequences of adaptation for budding yeast

Mutations that a population accumulates during evolution in one 'home' environment may cause fitness gains or losses in other environments. Such pleiotropic fitness effects determine the evolutionary fate of the population in variable environments and can lead to ecological specialization. It is unclear how the pleiotropic outcomes of evolution are shaped by the intrinsic randomness of the evolutionary process and by the deterministic variation in selection pressures across environments. Here, to address this question, we evolved 20 replicate populations of the yeast Saccharomyces cerevisiae in 11 laboratory environments and measured their fitness across multiple conditions. We found that evolution led to diverse pleiotropic fitness gains and losses, driven by multiple types of mutations. Approximately 60% of this variation is explained by the home environment of a clone and the most common parallel genetic changes, whereas about 40% is attributed to the stochastic accumulation of mutations whose pleiotropic effects are unpredictable. Although populations are typically specialized to their home environment, generalists also evolved in almost all of the conditions. Our results suggest that the mutations that accumulate during evolution incur a variety of pleiotropic costs and benefits with different probabilities. Thus, whether a population evolves towards a specialist or a generalist phenotype is heavily influenced by chance.

Biological amplification of low frequency mutations unravels laboratory culture history of the bio-threat agent Francisella tularensis

Challenges of investigating a suspected bio attack include establishing if microorganisms have been cultured to produce attack material and to identify their source. Addressing both issues, we have investigated genetic variations that emerge during laboratory culturing of the bacterial pathogen Francisella tularensis. Key aims were to identify genetic variations that are characteristic of laboratory culturing and explore the possibility of using biological amplification to identify genetic variation present at exceedingly low frequencies in a source sample. We used parallel serial passage experiments and high-throughput sequencing of F. tularensis to explore the genetic variation. We found that during early laboratory culture passages of F. tularensis, gene duplications emerged in the pathogen genome followed by single-nucleotide polymorphisms in genes for bacterial capsule synthesis. Based on a biological enrichment scheme and the use of high-throughput sequencing, we identified genetic variation that likely pre-existed in a source sample. The results support that capsule synthesis gene mutations are common during laboratory culture, and that a biological amplification strategy is useful for linking a F. tularensis sample to a specific laboratory variant among many highly similar variants.

The fitness challenge of studying molecular adaptation

Advances in bioinformatics and high-throughput genetic analysis increasingly allow us to predict the genetic basis of adaptive traits. These predictions can be tested and confirmed, but the molecular-level changes - i.e. the molecular adaptation - that link genetic differences to organism fitness remain generally unknown. In recent years, a series of studies have started to unpick the mechanisms of adaptation at the molecular level. In particular, this work has examined how changes in protein function, activity, and regulation cause improved organismal fitness. Key to addressing molecular adaptations is identifying systems and designing experiments that integrate changes in the genome, protein chemistry (molecular phenotype), and fitness. Knowledge of the molecular changes underpinning adaptations allow new insight into the constraints on, and repeatability of adaptations, and of the basis of non-additive interactions between adaptive mutations. Here we critically discuss a series of studies that examine the molecular-level adaptations that connect genetic changes and fitness.

CRISPR-Cas, a highly effective tool for genome editing in Clostridium saccharoperbutylacetonicum N1-4(HMT)

The solventogenic clostridia have long been known for their ability to convert sugars from complex feedstocks into commercially important solvents. Although the acetone-butanol-ethanol process fell out of favour decades ago, renewed interest in sustainability and ‘green’ chemistry has re-established our appetite for reviving technologies such as these, albeit with 21st century improvements. As CRISPR-Cas genome editing tools are being developed and applied to the solventogenic clostridia, their industrial potential is growing. Through integration of new pathways, the beneficial traits and historical track record of clostridial fermentation can be exploited to generate a much wider range of industrially relevant products. Here we show the application of genome editing using the endogenous CRISPR-Cas mechanism of Clostridium saccharoperbutylacetonicum N1-4(HMT), to generate a deletion, SNP and to integrate new DNA into the genome. These technological advancements pave the way for application of clostridial species to the production of an array of products.

Accurate Tracking of the Mutational Landscape of Diploid Hybrid Genomes

Mutations, recombinations, and genome duplications may promote genetic diversity and trigger evolutionary processes. However, quantifying these events in diploid hybrid genomes is challenging. Here, we present an integrated experimental and computational workflow to accurately track the mutational landscape of yeast diploid hybrids (MuLoYDH) in terms of single-nucleotide variants, small insertions/deletions, copy-number variants, aneuploidies, and loss-of-heterozygosity. Pairs of haploid Saccharomyces parents were combined to generate ancestor hybrids with phased genomes and varying levels of heterozygosity. These diploids were evolved under different laboratory protocols, in particular mutation accumulation experiments. Variant simulations enabled the efficient integration of competitive and standard mapping of short reads, depending on local levels of heterozygosity. Experimental validations proved the high accuracy and resolution of our computational approach. Finally, applying MuLoYDH to four different diploids revealed striking genetic background effects. Homozygous Saccharomyces cerevisiae showed a ∼4-fold higher mutation rate compared with its closely related species S. paradoxus. Intraspecies hybrids unveiled that a substantial fraction of the genome (∼250 bp per generation) was shaped by loss-of-heterozygosity, a process strongly inhibited in interspecies hybrids by high levels of sequence divergence between homologous chromosomes. In contrast, interspecies hybrids exhibited higher single-nucleotide mutation rates compared with intraspecies hybrids. MuLoYDH provided an unprecedented quantitative insight into the evolutionary processes that mold diploid yeast genomes and can be generalized to other genetic systems.