Tempo and mode of genome evolution in a 50,000-generation experiment

Single-cell genomics reveals population structures from in vitro evolutionary studies of Salmonella

Single-cell DNA sequencing has the potential to reveal detailed hierarchical structures in evolving populations of cells. Single cell approaches are increasingly used to study clonal evolution in human ageing and cancer but have not yet been deployed to study evolving clonal microbial populations. Here, we present an approach for single bacterial genomic analysis for in vitro evolution experiments using FACS isolation of individual bacteria followed by whole-genome amplification and sequencing. We apply this to the experimental evolution of a hypermutator strain of Salmonella in response to antibiotic stress (ciprofloxacin). By analysing sequence polymorphisms in individual cells from populations we identified the presence and prevalence of sub-populations which have acquired polymorphisms in genes previously demonstrated to be associated with ciprofloxacin susceptibility. We were also able to identify that the population exposed to antibiotic stress was able to develop resistance whilst maintaining diversity. This population structure could not be resolved from bulk sequence data, and our results show how high-throughput single-cell sequencing can enhance experimental studies of bacterial evolution.

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

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

Population genetics of clonally transmissible cancers

Populations of cancer cells are subject to the same core evolutionary processes as asexually reproducing, unicellular organisms. Transmissible cancers are particularly striking examples of these processes. These unusual cancers are clonal lineages that can spread through populations via physical transfer of living cancer cells from one host individual to another, and they have achieved long-term success in the colonization of at least eight different host species. Population genetic theory provides a useful framework for understanding the shift from a multicellular sexual animal into a unicellular asexual clone and its long-term effects on the genomes of these cancers. In this Review, we consider recent findings from transmissible cancer research with the goals of developing an evolutionarily informed perspective on transmissible cancers, examining possible implications for their long-term fate and identifying areas for future research on these exceptional lineages.

Quantitative contribution of the spacer length in the supercoiling-sensitivity of bacterial promoters

DNA supercoiling acts as a global transcriptional regulator in bacteria, but the promoter sequence or structural determinants controlling its effect remain unclear. It was previously proposed to modulate the torsional angle between the -10 and -35 hexamers, and thereby regulate the formation of the closed-complex depending on the length of the 'spacer' between them. Here, we develop a thermodynamic model of this notion based on DNA elasticity, providing quantitative and parameter-free predictions of the relative activation of promoters containing a short versus long spacer when the DNA supercoiling level is varied. The model is tested through an analysis of in vitro and in vivo expression assays of mutant promoters with variable spacer lengths, confirming its accuracy for spacers ranging from 15 to 19 nucleotides, except those of 16 nucleotides where other regulatory mechanisms likely overcome the effect of this specific step. An analysis at the whole-genome scale in Escherichia coli then demonstrates a significant effect of the spacer length on the genomic expression after transient or inheritable superhelical variations, validating the model's predictions. Altogether, this study shows an example of mechanical constraints associated to promoter binding by RNA Polymerase underpinning a basal and global regulatory mechanism.

Theragnostic application of nanoparticle and CRISPR against food-borne multi-drug resistant pathogens

Foodborne infection is one of the leading sources of infections spreading across the world. Foodborne pathogens are recognized as multidrug-resistant (MDR) pathogens posing a significant problem in the food industry and healthy consumers resulting in enhanced economic burden, and nosocomial infections. The continued search for enhanced microbial detection tools has piqued the interest of the CRISPR-Cas system and Nanoparticles. CRISPR-Cas system is present in the bacterial genome of some prokaryotes and is repurposed as a theragnostic tool against MDR pathogens. Nanoparticles and composites have also emerged as an efficient tool in theragnostic applications against MDR pathogens. The diagnostic limitations of the CRISPR-Cas system are believed to be overcome by a synergistic combination of the nanoparticles system and CRISPR-Cas using nanoparticles as vehicles. In this review, we have discussed the diagnostic application of CRISPR-Cas technologies along with their potential usage in applications like phage resistance, phage vaccination, strain typing, genome editing, and antimicrobial. we have also elucidated the antimicrobial and detection role of nanoparticles against foodborne MDR pathogens. Moreover, the novel combinatorial approach of CRISPR-Cas and nanoparticles for their synergistic effects in pathogen clearance and drug delivery vehicles has also been discussed.

Long-term experimental evolution decouples size and production costs in Escherichia coli

Body size covaries with population dynamics across life’s domains. Metabolism may impose fundamental constraints on the coevolution of size and demography, but experimental tests of the causal links remain elusive. We leverage a 60,000-generation experiment in which Escherichia coli populations evolved larger cells to examine intraspecific metabolic scaling and correlations with demographic parameters. Over the course of their evolution, the cells have roughly doubled in size relative to their ancestors. These larger cells have metabolic rates that are absolutely higher, but relative to their size, they are lower. Metabolic theory successfully predicted the relations between size, metabolism, and maximum population density, including support for Damuth’s law of energy equivalence, such that populations of larger cells achieved lower maximum densities but higher maximum biomasses than populations of smaller cells. The scaling of metabolism with cell size thus predicted the scaling of size with maximum population density. In stark contrast to standard theory, however, populations of larger cells grew faster than those of smaller cells, contradicting the fundamental and intuitive assumption that the costs of building new individuals should scale directly with their size. The finding that the costs of production can be decoupled from size necessitates a reevaluation of the evolutionary drivers and ecological consequences of biological size more generally.

Complex Ecotype Dynamics Evolve in Response to Fluctuating Resources

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

The role of sigma factor competition in bacterial adaptation under prolonged starvation

The study of adaptive microbial evolution in the laboratory can illuminate the genetic mechanisms of gaining fitness under a pre-defined set of selection factors. Laboratory evolution of bacteria under long-term starvation has gained importance in recent years because of its ability to uncover adaptive strategies that overcome prolonged nutrient limitation, a condition often encountered by natural microbes. In this evolutionary paradigm, bacteria are maintained in an energy-restricted environment in a growth phase called long-term stationary phase (LTSP). This phase is characterized by a stable, viable population size and highly dynamic genetic changes. Multiple independent iterations of LTSP evolution experiments have given rise to mutants that are slow-growing compared to the ancestor. Although the antagonistic regulation between rapid growth and the stress response is well-known in bacteria (especially Escherichia coli), the growth deficit of many LTSP-adapted mutants has not been explored in detail. In this review, I pinpoint the trade-off between growth and stress response as a dominant driver of evolutionary strategies under prolonged starvation. Focusing on mainly E. coli-based research, I discuss the various affectors and regulators of the competition between sigma factors to occupy their targets on the genome, and assess its effect on growth advantage in stationary phase (GASP). Finally, I comment on some crucial issues that hinder the progress of the field, including identification of novel metabolites in nutrient-depleted media, and the importance of using multidisciplinary research to resolve them.

Evolutionary transition from a single RNA replicator to a multiple replicator network

In prebiotic evolution, self-replicating molecules are believed to have evolved into complex living systems by expanding their information and functions open-endedly. Theoretically, such evolutionary complexification could occur through successive appearance of novel replicators that interact with one another to form replication networks. Here we perform long-term evolution experiments of RNA that replicates using a self-encoded RNA replicase. The RNA diversifies into multiple coexisting host and parasite lineages, whose frequencies in the population initially fluctuate and gradually stabilize. The final population, comprising five RNA lineages, forms a replicator network with diverse interactions, including cooperation to help the replication of all other members. These results support the capability of molecular replicators to spontaneously develop complexity through Darwinian evolution, a critical step for the emergence of life.

Predicting Parallelism and Quantifying Divergence in Microbial Evolution Experiments

The degree to which independent populations subjected to identical environmental conditions evolve in similar ways is a fundamental question in evolution. To address this question, microbial populations are often experimentally passaged in a given environment and sequenced to examine the tendency for similar mutations to repeatedly arise. However, there remains the need to develop an appropriate statistical framework to identify genes that acquired more mutations in one environment than in another (i.e., divergent evolution), genes that serve as genetic candidates of adaptation. Here, we develop a mathematical model to evaluate evolutionary outcomes among replicate populations in the same environment (i.e., parallel evolution), which can then be used to identify genes that contribute to divergent evolution. Applying this approach to data sets from evolve-and-resequence experiments, we found that the distribution of mutation counts among genes can be predicted as an ensemble of independent Poisson random variables with zero free parameters. Building on this result, we propose that the degree of divergent evolution at a given gene between populations from two different environments can be modeled as the difference between two Poisson random variables, known as the Skellam distribution. We then propose and apply a statistical test to identify specific genes that contribute to divergent evolution. By focusing on predicting patterns among replicate populations in a given environment, we are able to identify an appropriate test for divergence between environments that is grounded in first principles. IMPORTANCE There is currently no universally accepted framework for identifying genes that contribute to molecular divergence between microbial populations in different environments. To address this absence, we developed a null model to describe the distribution of mutation counts among genes. We find that divergent evolution within a given gene can be modeled as the absolute difference in the total number of mutations observed between two environments. This quantity is effectively captured by a probability distribution known as the Skellam distribution, providing an appropriate statistical test for researchers seeking to identify the set of genes that contribute to divergent evolution in microbial evolution experiments.

Mutators Enhance Adaptive Micro-Evolution in Pathogenic Microbes

Adaptation to the changing environmental conditions experienced within a host requires genetic diversity within a microbial population. Genetic diversity arises from mutations which occur due to DNA damage from exposure to exogenous environmental stresses or generated endogenously through respiration or DNA replication errors. As mutations can be deleterious, a delicate balance must be obtained between generating enough mutations for micro-evolution to occur while maintaining fitness and genomic integrity. Pathogenic microorganisms can actively modify their mutation rate to enhance adaptive micro-evolution by increasing expression of error-prone DNA polymerases or by mutating or decreasing expression of genes required for DNA repair. Strains which exhibit an elevated mutation rate are termed mutators. Mutators are found in varying prevalence in clinical populations where large-effect beneficial mutations enhance survival and are predominately caused by defects in the DNA mismatch repair (MMR) pathway. Mutators can facilitate the emergence of antibiotic resistance, allow phenotypic modifications to prevent recognition and destruction by the host immune system and enable switching to metabolic and cellular morphologies better able to survive in the given environment. This review will focus on recent advances in understanding the phenotypic and genotypic changes occurring in MMR mutators in both prokaryotic and eukaryotic pathogens.

Anatomy promotes neutral coexistence of strains in the human skin microbiome

What enables strains of the same species to coexist in a microbiome? Here, we investigate whether host anatomy can explain strain co-residence of Cutibacterium acnes, the most abundant species on human skin. We reconstruct on-person evolution and migration using whole-genome sequencing of C. acnes colonies acquired from healthy subjects, including from individual skin pores, and find considerable spatial structure at the level of pores. Although lineages (sets of colonies separated by <100 mutations) with in vitro fitness differences coexist within centimeter-scale regions, each pore is dominated by a single lineage. Moreover, colonies from a pore typically have identical genomes. An absence of adaptive signatures suggests a genotype-independent source of low within-pore diversity. We therefore propose that pore anatomy imposes random single-cell bottlenecks; the resulting population fragmentation reduces competition and promotes coexistence. Our findings suggest that therapeutic interventions involving pore-dwelling species might focus on removing resident populations over optimizing probiotic fitness.

Comparative Population Genetics in the Human Gut Microbiome

Genetic variation in the human gut microbiome is responsible for conferring a number of crucial phenotypes like the ability to digest food and metabolize drugs. Yet, our understanding of how this variation arises and is maintained remains relatively poor. Thus, the microbiome remains a largely untapped resource, as the large number of coexisting species in the microbiome presents a unique opportunity to compare and contrast evolutionary processes across species to identify universal trends and deviations. Here we outline features of the human gut microbiome that, while not unique in isolation, as an assemblage make it a system with unparalleled potential for comparative population genomics studies. We consciously take a broad view of comparative population genetics, emphasizing how sampling a large number of species allows researchers to identify universal evolutionary dynamics in addition to new genes, which can then be leveraged to identify exceptional species that deviate from general patterns. To highlight the potential power of comparative population genetics in the microbiome, we reanalyze patterns of purifying selection across ∼40 prevalent species in the human gut microbiome to identify intriguing trends which highlight functional categories in the microbiome that may be under more or less constraint.

Evolutionary Dynamics of Asexual Hypermutators Adapting to a Novel Environment

How microbes adapt to a novel environment is a central question in evolutionary biology. Although adaptive evolution must be fueled by beneficial mutations, whether higher mutation rates facilitate the rate of adaptive evolution remains unclear. To address this question, we cultured Escherichia coli hypermutating populations, in which a defective methyl-directed mismatch repair pathway causes a 140-fold increase in single-nucleotide mutation rates. In parallel with wild-type E. coli, populations were cultured in tubes containing Luria-Bertani broth, a complex medium known to promote the evolution of subpopulation structure. After 900 days of evolution, in three transfer schemes with different population-size bottlenecks, hypermutators always exhibited similar levels of improved fitness as controls. Fluctuation tests revealed that the mutation rates of hypermutator lines converged evolutionarily on those of wild-type populations, which may have contributed to the absence of fitness differences. Further genome-sequence analysis revealed that, although hypermutator populations have higher rates of genomic evolution, this largely reflects strong genetic linkage. Despite these linkage effects, the evolved population exhibits parallelism in fixed mutations, including those potentially related to biofilm formation, transcription regulation, and mutation-rate evolution. Together, these results are generally inconsistent with a hypothesized positive relationship between the mutation rate and the adaptive speed of evolution, and provide insight into how clonal adaptation occurs in novel environments.

Emerging Adaptive Strategies Under Temperature Fluctuations in a Laboratory Evolution Experiment of Escherichia Coli

Generalists and specialists are types of strategies individuals can employ that can evolve in fluctuating environments depending on the extremity and periodicity of the fluctuation. To evaluate whether the evolution of specialists or generalists occurs under environmental fluctuation regimes with different levels of periodicity, 24 populations of Escherichia coli underwent laboratory evolution with temperatures alternating between 15 and 43°C in three fluctuation regimes: two periodic regimes dependent on culture's cell density and one random (non-periodic) regime with no such dependency, serving as a control. To investigate contingencies on the genetic background, we seeded our experiment with two different strains. After the experiment, growth rate measurements at the two temperatures showed that the evolution of specialists was favored in the random regime, while generalists were favored in the periodic regimes. Whole genome sequencing demonstrated that several gene mutations were selected in parallel in the evolving populations with some dependency on the starting genetic background. Given the genes mutated, we hypothesized that the driving force behind the observed adaptations is the restoration of the internal physiology of the starting strains' unstressed states at 37°C, which may be a means of improving fitness in the new environments. Phenotypic array measurements supported our hypothesis by demonstrating a tendency of the phenotypic response of the evolved strains to move closer to the starting strains' response at the optimum of 37°C, especially for strains classified as generalists.

Adaptation and genomic erosion in fragmented Pseudomonas aeruginosa populations in the sinuses of people with cystic fibrosis

Pseudomonas aeruginosa notoriously adapts to the airways of people with cystic fibrosis (CF), yet how infection-site biogeography and associated evolutionary processes vary as lifelong infections progress remains unclear. Here we test the hypothesis that early adaptations promoting aggregation influence evolutionary-genetic trajectories by examining longitudinal P. aeruginosa from the sinuses of six adults with CF. Highly host-adapted lineages harbored mutator genotypes displaying signatures of early genome degradation associated with recent host restriction. Using an advanced imaging technique (MiPACT-HCR [microbial identification after passive clarity technique]), we find population structure tracks with genome degradation, with the most host-adapted, genome-degraded P. aeruginosa (the mutators) residing in small, sparse aggregates. We propose that following initial adaptive evolution in larger populations under strong selection for aggregation, P. aeruginosa persists in small, fragmented populations that experience stronger effects of genetic drift. These conditions enrich for mutators and promote degenerative genome evolution. Our findings underscore the importance of infection-site biogeography to pathogen evolution.

The Transporter-Mediated Cellular Uptake and Efflux of Pharmaceutical Drugs and Biotechnology Products: How and Why Phospholipid Bilayer Transport Is Negligible in Real Biomembranes

Over the years, my colleagues and I have come to realise that the likelihood of pharmaceutical drugs being able to diffuse through whatever unhindered phospholipid bilayer may exist in intact biological membranes in vivo is vanishingly low. This is because (i) most real biomembranes are mostly protein, not lipid, (ii) unlike purely lipid bilayers that can form transient aqueous channels, the high concentrations of proteins serve to stop such activity, (iii) natural evolution long ago selected against transport methods that just let any undesirable products enter a cell, (iv) transporters have now been identified for all kinds of molecules (even water) that were once thought not to require them, (v) many experiments show a massive variation in the uptake of drugs between different cells, tissues, and organisms, that cannot be explained if lipid bilayer transport is significant or if efflux were the only differentiator, and (vi) many experiments that manipulate the expression level of individual transporters as an independent variable demonstrate their role in drug and nutrient uptake (including in cytotoxicity or adverse drug reactions). This makes such transporters valuable both as a means of targeting drugs (not least anti-infectives) to selected cells or tissues and also as drug targets. The same considerations apply to the exploitation of substrate uptake and product efflux transporters in biotechnology. We are also beginning to recognise that transporters are more promiscuous, and antiporter activity is much more widespread, than had been realised, and that such processes are adaptive (i.e., were selected by natural evolution). The purpose of the present review is to summarise the above, and to rehearse and update readers on recent developments. These developments lead us to retain and indeed to strengthen our contention that for transmembrane pharmaceutical drug transport "phospholipid bilayer transport is negligible".

A modified fluctuation assay reveals a natural mutator phenotype that drives mutation spectrum variation within Saccharomyces cerevisiae

Although studies of Saccharomyces cerevisiae have provided many insights into mutagenesis and DNA repair, most of this work has focused on a few laboratory strains. Much less is known about the phenotypic effects of natural variation within S. cerevisiae’s DNA repair pathways. Here, we use natural polymorphisms to detect historical mutation spectrum differences among several wild and domesticated S. cerevisiae strains. To determine whether these differences are likely caused by genetic mutation rate modifiers, we use a modified fluctuation assay with a CAN1 reporter to measure de novo mutation rates and spectra in 16 of the analyzed strains. We measure a 10-fold range of mutation rates and identify two strains with distinctive mutation spectra. These strains, known as AEQ and AAR, come from the panel’s ‘Mosaic beer’ clade and share an enrichment for C > A mutations that is also observed in rare variation segregating throughout the genomes of several Mosaic beer and Mixed origin strains. Both AEQ and AAR are haploid derivatives of the diploid natural isolate CBS 1782, whose rare polymorphisms are enriched for C > A as well, suggesting that the underlying mutator allele is likely active in nature. We use a plasmid complementation test to show that AAR and AEQ share a mutator allele in the DNA repair gene OGG1, which excises 8-oxoguanine lesions that can cause C > A mutations if left unrepaired.

High-throughput characterization of mutations in genes that drive clonal evolution using multiplex adaptome capture sequencing

Understanding how cells are likely to evolve can guide medical interventions and bioengineering efforts that must contend with unwanted mutations. The adaptome of a cell-the neighborhood of genetic changes that are most likely to drive adaptation in a given environment-can be mapped by tracking rare beneficial variants during the early stages of clonal evolution. We used multiplex adaptome capture sequencing (mAdCap-seq), a procedure that combines unique molecular identifiers and hybridization-based enrichment, to characterize mutations in eight Escherichia coli genes known to be under selection in a laboratory environment. We tracked 301 mutations at frequencies as low as 0.01% and inferred the fitness effects of 240 of these mutations. There were distinct molecular signatures of selection on protein structure and function for the three genes with the most beneficial mutations. Our results demonstrate how mAdCap-seq can be used to deeply profile a targeted portion of a cell's adaptome.

The cost of evolved constitutive lac gene expression is usually, but not always, maintained during evolution of generalist populations

Beneficial mutations can become costly following an environmental change. Compensatory mutations can relieve these costs, while not affecting the selected function, so that the benefits are retained if the environment shifts back to be similar to the one in which the beneficial mutation was originally selected. Compensatory mutations have been extensively studied in the context of antibiotic resistance, responses to specific genetic perturbations, and in the determination of interacting gene network components. Few studies have focused on the role of compensatory mutations during more general adaptation, especially as the result of selection in fluctuating environments where adaptations to different environment components may often involve trade-offs. We examine whether costs of a mutation in lacI, which deregulated the expression of the lac operon in evolving populations of Escherichia coli bacteria, were compensated. This mutation occurred in multiple replicate populations selected in environments that fluctuated between growth on lactose, where the mutation was beneficial, and on glucose, where it was deleterious. We found that compensation for the cost of the lacI mutation was rare, but, when it did occur, it did not negatively affect the selected benefit. Compensation was not more likely to occur in a particular evolution environment. Compensation has the potential to remove pleiotropic costs of adaptation, but its rarity indicates that the circumstances to bring about the phenomenon may be peculiar to each individual or impeded by other selected mutations.

Role of the Discriminator Sequence in the Supercoiling Sensitivity of Bacterial Promoters

DNA supercoiling acts as a global transcriptional regulator that contributes to the rapid transcriptional response of bacteria to many environmental changes. Although a large fraction of promoters from phylogenetically distant species respond to superhelical variations, the sequence or structural determinants of this behavior remain elusive. Here, we focus on the sequence of the "discriminator" element that was shown to modulate this response in several promoters. We develop a quantitative thermodynamic model of this regulatory effect, focusing on open complex formation during transcription initiation independently from promoter-specific regulatory proteins. We analyze previous and new expression data and show that the model predictions quantitatively match the in vitro and in vivo supercoiling response of selected promoters with mutated discriminator sequences. We then test the universality of this mechanism by a statistical analysis of promoter sequences from transcriptomes of phylogenetically distant bacteria under conditions of supercoiling variations (i) by gyrase inhibitors, (ii) by environmental stresses, or (iii) inherited in the longest-running evolution experiment. In all cases, we identify a robust and significant sequence signature in the discriminator region, suggesting that supercoiling-modulated promoter opening underpins a ubiquitous regulatory mechanism in the prokaryotic kingdom based on the fundamental mechanical properties of DNA and its basal interaction with RNA polymerase. IMPORTANCE In this study, we highlight the role of the discriminator as a global sensor of supercoiling variations and propose the first quantitative regulatory model of this principle, based on the specific step of promoter opening during transcription initiation. It defines the predictive rule by which DNA supercoiling quantitatively modulates the expression rate of bacterial promoters, depending on the G/C content of their discriminator and independently from promoter-specific regulatory proteins. This basal mechanism affects a wide range of species, which is tested by an extensive analysis of global high-throughput expression data. Altogether, ours results confirm and provide a quantitative framework for the long-proposed notion that the discriminator sequence is a significant determinant of promoter supercoiling sensitivity, underpinning the ubiquitous regulatory action of DNA supercoiling on the core transcriptional machinery, in particular in response to quick environmental changes.

Changes in the distribution of fitness effects and adaptive mutational spectra following a single first step towards adaptation

Historical contingency and diminishing returns epistasis have been typically studied for relatively divergent genotypes and/or over long evolutionary timescales. Here, we use Saccharomyces cerevisiae to study the extent of diminishing returns and the changes in the adaptive mutational spectra following a single first adaptive mutational step. We further evolve three clones that arose under identical conditions from a common ancestor. We follow their evolutionary dynamics by lineage tracking and determine adaptive outcomes using fitness assays and whole genome sequencing. We find that diminishing returns manifests as smaller fitness gains during the 2^nd step of adaptation compared to the 1^st step, mainly due to a compressed distribution of fitness effects. We also find that the beneficial mutational spectra for the 2^nd adaptive step are contingent on the 1^st step, as we see both shared and diverging adaptive strategies. Finally, we find that adaptive loss-of-function mutations, such as nonsense and frameshift mutations, are less common in the second step of adaptation than in the first step.

Analyses of both natural and experimental evolution suggest that adaptation depends on the evolutionary past and adaptive potential decreases over time. Here, by tracking yeast adaptation with DNA barcoding, the authors show that such evolutionary phenomena can be observed even after a single adaptive step.

Forecasting of phenotypic and genetic outcomes of experimental evolution in Pseudomonas protegens

Experimental evolution with microbes is often highly repeatable under identical conditions, suggesting the possibility to predict short-term evolution. However, it is not clear to what degree evolutionary forecasts can be extended to related species in non-identical environments, which would allow testing of general predictive models and fundamental biological assumptions. To develop an extended model system for evolutionary forecasting, we used previous data and models of the genotype-to-phenotype map from the wrinkly spreader system in Pseudomonas fluorescens SBW25 to make predictions of evolutionary outcomes on different biological levels for Pseudomonas protegens Pf-5. In addition to sequence divergence (78% amino acid and 81% nucleotide identity) for the genes targeted by mutations, these species also differ in the inability of Pf-5 to make cellulose, which is the main structural basis for the adaptive phenotype in SBW25. The experimental conditions were changed compared to the SBW25 system to test if forecasts were extendable to a non-identical environment. Forty-three mutants with increased ability to colonize the air-liquid interface were isolated, and the majority had reduced motility and was partly dependent on the Pel exopolysaccharide as a structural component. Most (38/43) mutations are expected to disrupt negative regulation of the same three diguanylate cyclases as in SBW25, with a smaller number of mutations in promoter regions, including an uncharacterized polysaccharide synthase operon. A mathematical model developed for SBW25 predicted the order of the three main pathways and the genes targeted by mutations, but differences in fitness between mutants and mutational biases also appear to influence outcomes. Mutated regions in proteins could be predicted in most cases (16/22), but parallelism at the nucleotide level was low and mutational hot spot sites were not conserved. This study demonstrates the potential of short-term evolutionary forecasting in experimental populations and provides testable predictions for evolutionary outcomes in other Pseudomonas species.

Reversible gene silencing through frameshift indels and frameshift scars provide adaptive plasticity for Mycobacterium tuberculosis

Mycobacterium tuberculosis can adapt to changing environments by non-heritable mechanisms. Frame-shifting insertions and deletions (indels) may also participate in adaptation through gene disruption, which could be reversed by secondary introduction of a frame-restoring indel. We present ScarTrek, a program that scans genomic data for indels, including those that together disrupt and restore a gene's reading frame, producing "frame-shift scars" suggestive of reversible gene inactivation. We use ScarTrek to analyze 5977 clinical M. tuberculosis isolates. We show that indel frequency inversely correlates with genomic linguistic complexity and varies with gene-position and gene-essentiality. Using ScarTrek, we detect 74 unique frame-shift scars in 48 genes, with a 3.74% population-level incidence of unique scar events. We find multiple scars in the ESX-1 gene cluster. Six scars show evidence of convergent evolution while the rest shared a common ancestor. Our results suggest that sequential indels are a mechanism for reversible gene silencing and adaptation in M. tuberculosis.

Synergistic epistasis enhances the co-operativity of mutualistic interspecies interactions

Early evolution of mutualism is characterized by big and predictable adaptive changes, including the specialization of interacting partners, such as through deleterious mutations in genes not required for metabolic cross-feeding. We sought to investigate whether these early mutations improve cooperativity by manifesting in synergistic epistasis between genomes of the mutually interacting species. Specifically, we have characterized evolutionary trajectories of syntrophic interactions of Desulfovibrio vulgaris (Dv) with Methanococcus maripaludis (Mm) by longitudinally monitoring mutations accumulated over 1000 generations of nine independently evolved communities with analysis of the genotypic structure of one community down to the single-cell level. We discovered extensive parallelism across communities despite considerable variance in their evolutionary trajectories and the perseverance within many evolution lines of a rare lineage of Dv that retained sulfate-respiration (SR+) capability, which is not required for metabolic cross-feeding. An in-depth investigation revealed that synergistic epistasis across pairings of Dv and Mm genotypes had enhanced cooperativity within SR- and SR+ assemblages, enabling their coexistence within the same community. Thus, our findings demonstrate that cooperativity of a mutualism can improve through synergistic epistasis between genomes of the interacting species, enabling the coexistence of mutualistic assemblages of generalists and their specialized variants.

Dynamics of extended-spectrum beta-lactamase-producing Enterobacterales colonization in long-term carriers following travel abroad

Travel to tropical regions is associated with high risk of acquiring extended-spectrum beta-lactamase-producing Enterobacterales (ESBL-E) that are typically cleared in less than 3 months following return. The conditions leading to persistent carriage that exceeds 3 months in some travellers require investigation. Whole-genome sequencing (Illumina MiSeq) was performed on the 82 ESBL-E isolates detected upon return and 1, 2, 3, 6 and 12 months later from the stools of 11 long-term (>3 months) ESBL-E carriers following travel abroad. One to five different ESBL Escherichia coli strains were detected per traveller upon return, and this diminished to one after 3 months. Long-term carriage was due to the presence of the same ESBL E. coli strain, for more than 3 months, in 9 out of 11 travellers, belonging to epidemic sequence type complexes (STc 10, 14, 38, 69, 131 and 648). The mean carriage duration of strains belonging to phylogroups B2/D/F, associated with extra-intestinal virulence, was higher than that for commensal-associated A/B1/E phylogroups (3.5 vs 0.5 months, P=0.021). Genes encoding iron capture systems (fyuA, irp), toxins (senB, sat), adhesins (flu, daaF, afa/nfaE, pap, ecpA) and colicin (cjrA) were more often present in persistent strains than in transient ones. Single-nucleotide polymorphism (SNP) analysis in persistent strains showed a maximum divergence of eight SNPs over 12 months without signs of adaptation. Genomic plasticity was observed during the follow-up with the loss or gain of mobile genetic elements such as plasmids, integrons and/or transposons that may contain resistance genes at different points in the follow-up. Long-term colonization of ESBL-E following travel is primarily due to the acquisition of E. coli strains belonging to epidemic clones and harbouring ‘virulence genes’, allowing good adaptation to the intestinal microbiota.

Choice of β-Lactam Resistance Pathway Depends Critically on Initial Antibiotic Concentration

Antibiotic resistance trajectories with different final resistance may critically depend on the first mutation, due to epistatic interactions. Here, we study the effect of mutation bias and the concentration-dependent effects on fitness of two clinically important mutations in TEM-1 β-lactamase in initiating alternative trajectories to cefotaxime resistance. We show that at low cefotaxime concentrations, the R164S mutation (a mutation of arginine to serine at position 164), which confers relatively low resistance, is competitively superior to the G238S mutation, conferring higher resistance, thus highlighting a critical influence of antibiotic concentration on long-term resistance evolution.

Population and comparative genetics of thermotolerance divergence between yeast species

Many familiar traits in the natural world-from lions' manes to the longevity of bristlecone pine trees-arose in the distant past, and have long since fixed in their respective species. A key challenge in evolutionary genetics is to figure out how and why species-defining traits have come to be. We used the thermotolerance growth advantage of the yeast Saccharomyces cerevisiae over its sister species Saccharomyces paradoxus as a model for addressing these questions. Analyzing loci at which the S. cerevisiae allele promotes thermotolerance, we detected robust evidence for positive selection, including amino acid divergence between the species and conservation within S. cerevisiae populations. Because such signatures were particularly strong at the chromosome segregation gene ESP1, we used this locus as a case study for focused mechanistic follow-up. Experiments revealed that, in culture at high temperature, the S. paradoxus ESP1 allele conferred a qualitative defect in biomass accumulation and cell division relative to the S. cerevisiae allele. Only genetic divergence in the ESP1 coding region mattered phenotypically, with no functional impact detectable from the promoter. Our data support a model in which an ancient ancestor of S. cerevisiae, under selection to boost viability at high temperature, acquired amino acid variants at ESP1 and many other loci, which have been constrained since then. Complex adaptations of this type hold promise as a paradigm for interspecies genetics, especially in deeply diverged traits that may have taken millions of years to evolve.

Divergent Evolution of Mutation Rates and Biases in the Long-Term Evolution Experiment with Escherichia coli

All organisms encode enzymes that replicate, maintain, pack, recombine, and repair their genetic material. For this reason, mutation rates and biases also evolve by mutation, variation, and natural selection. By examining metagenomic time series of the Lenski long-term evolution experiment (LTEE) with Escherichia coli (Good BH, McDonald MJ, Barrick JE, Lenski RE, Desai MM. 2017. The dynamics of molecular evolution over 60,000 generations. Nature 551(7678):45–50.), we find that local mutation rate variation has evolved during the LTEE. Each LTEE population has evolved idiosyncratic differences in their rates of point mutations, indels, and mobile element insertions, due to the fixation of various hypermutator and antimutator alleles. One LTEE population, called Ara+3, shows a strong, symmetric wave pattern in its density of point mutations, radiating from the origin of replication. This pattern is largely missing from the other LTEE populations, most of which evolved missense, indel, or structural mutations in topA, fis, and dusB—loci that all affect DNA topology. The distribution of mutations in those genes over time suggests epistasis and historical contingency in the evolution of DNA topology, which may have in turn affected local mutation rates. Overall, the replicate populations of the LTEE have largely diverged in their mutation rates and biases, even though they have adapted to identical abiotic conditions.

Genomic evolution of antibiotic resistance is contingent on genetic background following a long-term experiment with Escherichia coli

Antibiotic resistance is a growing health concern. Efforts to control resistance would benefit from an improved ability to forecast when and how it will evolve. Epistatic interactions between mutations can promote divergent evolutionary trajectories, which complicates our ability to predict evolution. We recently showed that differences between genetic backgrounds can lead to idiosyncratic responses in the evolvability of phenotypic resistance, even among closely related Escherichia coli strains. In this study, we examined whether a strain's genetic background also influences the genotypic evolution of resistance. Do lineages founded by different genotypes take parallel or divergent mutational paths to achieve their evolved resistance states? We addressed this question by sequencing the complete genomes of antibiotic-resistant clones that evolved from several different genetic starting points during our earlier experiments. We first validated our statistical approach by quantifying the specificity of genomic evolution with respect to antibiotic treatment. As expected, mutations in particular genes were strongly associated with each drug. Then, we determined that replicate lines evolved from the same founding genotypes had more parallel mutations at the gene level than lines evolved from different founding genotypes, although these effects were more subtle than those showing antibiotic specificity. Taken together with our previous work, we conclude that historical contingency can alter both genotypic and phenotypic pathways to antibiotic resistance.

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.

Experimental evolution of extremophile resistance to ionizing radiation

A growing number of known species possess a remarkable characteristic - extreme resistance to the effects of ionizing radiation (IR). This review examines our current understanding of how organisms can adapt to and survive exposure to IR, one of the most toxic stressors known. The study of natural extremophiles such as Deinococcus radiodurans has revealed much. However, the evolution of Deinococcus was not driven by IR. Another approach, pioneered by Evelyn Witkin in 1946, is to utilize experimental evolution. Contributions to the IR-resistance phenotype affect multiple aspects of cell physiology, including DNA repair, removal of reactive oxygen species, the structure and packaging of DNA and the cell itself, and repair of iron-sulfur centers. Based on progress to date, we overview the diversity of mechanisms that can contribute to biological IR resistance arising as a result of either natural or experimental evolution.

Compensatory evolution of Pseudomonas aeruginosa's slow growth phenotype suggests mechanisms of adaptation in cystic fibrosis

Long-term infection of the airways of cystic fibrosis patients with Pseudomonas aeruginosa is often accompanied by a reduction in bacterial growth rate. This reduction has been hypothesised to increase within-patient fitness and overall persistence of the pathogen. Here, we apply adaptive laboratory evolution to revert the slow growth phenotype of P. aeruginosa clinical strains back to a high growth rate. We identify several evolutionary trajectories and mechanisms leading to fast growth caused by transcriptional and mutational changes, which depend on the stage of adaptation of the strain. Return to high growth rate increases antibiotic susceptibility, which is only partially dependent on reversion of mutations or changes in the transcriptional profile of genes known to be linked to antibiotic resistance. We propose that similar mechanisms and evolutionary trajectories, in reverse direction, may be involved in pathogen adaptation and the establishment of chronic infections in the antibiotic-treated airways of cystic fibrosis patients.

Maintenance of Metabolic Plasticity despite Relaxed Selection in a Long-Term Evolution Experiment with Escherichia coli

AbstractTraits that are unused in a given environment are subject to processes that tend to erode them, leading to reduced fitness in other environments. Although this general tendency is clear, we know much less about why some traits are lost while others are retained and about the roles of mutation and selection in generating different responses. We addressed these issues by examining populations of a facultative anaerobe, Escherichia coli, that have evolved for >30 years in the presence of oxygen, with relaxed selection for anaerobic growth and the associated metabolic plasticity. We asked whether evolution led to the loss, improvement, or maintenance of anaerobic growth, and we analyzed gene expression and mutational data sets to understand the outcomes. We identified genomic signatures of both positive and purifying selection on aerobic-specific genes, while anaerobic-specific genes showed clear evidence of relaxed selection. We also found parallel evolution at two interacting loci that regulate anaerobic growth. We competed the ancestor and evolved clones from each population in an anoxic environment, and we found that anaerobic fitness had not decayed, despite relaxed selection. In summary, relaxed selection does not necessarily reduce an organism's fitness in other environments. Instead, the genetic architecture of the traits under relaxed selection and their correlations with traits under positive and purifying selection may sometimes determine evolutionary outcomes.

Mutations as Levy flights

Data from a long time evolution experiment with Escherichia Coli and from a large study on copy number variations in subjects with European ancestry are analyzed in order to argue that mutations can be described as Levy flights in the mutation space. These Levy flights have at least two components: random single-base substitutions and large DNA rearrangements. From the data, we get estimations for the time rates of both events and the size distribution function of large rearrangements.

High-throughput functional variant screens via in vivo production of single-stranded DNA

Creating and characterizing individual genetic variants remains limited in scale, compared to the tremendous variation both existing in nature and envisioned by genome engineers. Here we introduce retron library recombineering (RLR), a methodology for high-throughput functional screens that surpasses the scale and specificity of CRISPR-Cas methods. We use the targeted reverse-transcription activity of retrons to produce single-stranded DNA (ssDNA) in vivo, incorporating edits at >90% efficiency and enabling multiplexed applications. RLR simultaneously introduces many genomic variants, producing pooled and barcoded variant libraries addressable by targeted deep sequencing. We use RLR for pooled phenotyping of synthesized antibiotic resistance alleles, demonstrating quantitative measurement of relative growth rates. We also perform RLR using the sheared genomic DNA of an evolved bacterium, experimentally querying millions of sequences for causal variants, demonstrating that RLR is uniquely suited to utilize large pools of natural variation. Using ssDNA produced in vivo for pooled experiments presents avenues for exploring variation across the genome.

Quantifying the Mutational Robustness of Protein-Coding Genes

We use large-scale mutagenesis data and computer simulations to quantify the mutational robustness of protein-coding genes by taking into account constraints arising from protein function and the genetic code. Analyses of the distribution of amino acid substitutions from 18 mutagenesis studies revealed an average of 45% of neutral variants; while mutagenesis data of 12 proteins artificially designed under no other constraints but stability, reach an average of 60%. Simulations using a lattice protein model allow us to contrast these estimates to the expected mutational robustness of protein families by generating unbiased samples of foldable sequences, which we find to have 30% of neutral variants. In agreement with mutagenesis data of designed proteins, the model shows that maximally robust protein families might access up to twice the amount of neutral variants observed in the unbiased samples (i.e. 60%). A biophysical model of protein-ligand binding suggests that constraints associated to molecular function have only a moderate impact on robustness of approximately 5 to 10% of neutral variants; and that the direction of this effect depends on the relation between functional performance and thermodynamic stability. Although the genetic code constraints the access of a gene's nucleotide sequence to only 30% of the full distribution of amino acid mutations, it provides an extra 15 to 20% of neutral variants to the estimations above, such that the expected, observed, and maximal robustness of protein-coding genes are approximately 50, 65, and 75%, respectively. We discuss our results in the light of three main hypothesis put forward to explain the existence of mutationally robust genes.

Dynamics of bacterial adaptation

Determining pattern in the dynamics of population evolution is a long-standing focus of evolutionary biology. Complementing the study of natural populations, microbial laboratory evolution experiments have become an important tool for addressing these dynamics because they allow detailed and replicated analysis of evolution in response to controlled environmental and genetic conditions. Key findings include a tendency for smoothly declining rates of adaptation during selection in constant environments, at least in part a reflection of antagonism between accumulating beneficial mutations, and a large number of beneficial mutations available to replicate populations leading to significant, but relatively low genetic parallelism, even as phenotypic characteristics show high similarity. Together, there is a picture of adaptation as a process with a varied and largely unpredictable genetic basis leading to much more similar phenotypic outcomes. Increasing sophistication of sequencing and genetic tools will allow insight into mechanisms behind these and other patterns.

Adaptive differentiation and rapid evolution of a soil bacterium along a climate gradient

Increasing evidence suggests that evolutionary processes frequently shape ecological patterns; however, most microbiome studies thus far have focused on only the ecological responses of these communities. By using parallel field experiments and focusing in on a model soil bacterium, we showed that bacterial “species” are differentially adapted to local climates, leading to changes in their composition. Furthermore, we detected strain-level evolution, providing direct evidence that both ecological and evolutionary processes operate on annual timescales. The consideration of eco-evolutionary dynamics may therefore be important to understand the response of soil microbiomes to future environmental change.

Microbial community responses to environmental change are largely associated with ecological processes; however, the potential for microbes to rapidly evolve and adapt remains relatively unexplored in natural environments. To assess how ecological and evolutionary processes simultaneously alter the genetic diversity of a microbiome, we conducted two concurrent experiments in the leaf litter layer of soil over 18 mo across a climate gradient in Southern California. In the first experiment, we reciprocally transplanted microbial communities from five sites to test whether ecological shifts in ecotypes of the abundant bacterium, Curtobacterium, corresponded to past adaptive differentiation. In the transplanted communities, ecotypes converged toward that of the native communities growing on a common litter substrate. Moreover, these shifts were correlated with community-weighted mean trait values of the Curtobacterium ecotypes, indicating that some of the trait variation among ecotypes could be explained by local adaptation to climate conditions. In the second experiment, we transplanted an isogenic Curtobacterium strain and tracked genomic mutations associated with the sites across the same climate gradient. Using a combination of genomic and metagenomic approaches, we identified a variety of nonrandom, parallel mutations associated with transplantation, including mutations in genes related to nutrient acquisition, stress response, and exopolysaccharide production. Together, the field experiments demonstrate how both demographic shifts of previously adapted ecotypes and contemporary evolution can alter the diversity of a soil microbiome on the same timescale.

Changes in Cell Size and Shape during 50,000 Generations of Experimental Evolution with Escherichia coli

Bacteria adopt a wide variety of sizes and shapes, with many species exhibiting stereotypical morphologies. How morphology changes, and over what timescales, is less clear. Previous work examining cell morphology in an experiment with Escherichia coli showed that populations evolved larger cells and, in some cases, cells that were less rod-like. That experiment has now run for over two more decades. Meanwhile, genome sequence data are available for these populations, and new computational methods enable high-throughput microscopic analyses. In this study, we measured stationary-phase cell volumes for the ancestor and 12 populations at 2,000, 10,000, and 50,000 generations, including measurements during exponential growth at the last time point. We measured the distribution of cell volumes for each sample using a Coulter counter and microscopy, the latter of which also provided data on cell shape. Our data confirm the trend toward larger cells while also revealing substantial variation in size and shape across replicate populations. Most populations first evolved wider cells but later reverted to the ancestral length-to-width ratio. All but one population evolved mutations in rod shape maintenance genes. We also observed many ghost-like cells in the only population that evolved the novel ability to grow on citrate, supporting the hypothesis that this lineage struggles with maintaining balanced growth. Lastly, we show that cell size and fitness remain correlated across 50,000 generations. Our results suggest that larger cells are beneficial in the experimental environment, while the reversion toward ancestral length-to-width ratios suggests partial compensation for the less favorable surface area-to-volume ratios of the evolved cells.IMPORTANCE Bacteria exhibit great morphological diversity, yet we have only a limited understanding of how their cell sizes and shapes evolve and of how these features affect organismal fitness. This knowledge gap reflects, in part, the paucity of the fossil record for bacteria. In this study, we revived and analyzed samples extending over 50,000 generations from 12 populations of experimentally evolving Escherichia coli to investigate the relation between cell size, shape, and fitness. Using this "frozen fossil record," we show that all 12 populations evolved larger cells concomitant with increased fitness, with substantial heterogeneity in cell size and shape across the replicate lines. Our work demonstrates that cell morphology can readily evolve and diversify, even among populations living in identical environments.

Selection Maintains Protein Interactome Resilience in the Long-Term Evolution Experiment with Escherichia coli

Most cellular functions are carried out by a dynamic network of interacting proteins. An open question is whether the network properties of protein interactomes represent phenotypes under natural selection. One proposal is that protein interactomes have evolved to be resilient, such that they tend to maintain connectivity when proteins are removed from the network. This hypothesis predicts that interactome resilience should be maintained by natural selection during long-term experimental evolution. I tested this prediction by modeling the evolution of protein-protein interaction (PPI) networks in Lenski's long-term evolution experiment with Escherichia coli (LTEE). In this test, I removed proteins affected by nonsense, insertion, deletion, and transposon mutations in evolved LTEE strains, and measured the resilience of the resulting networks. I compared the rate of change of network resilience in each LTEE population to the rate of change of network resilience for corresponding randomized networks. The evolved PPI networks are significantly more resilient than networks in which random proteins have been deleted. Moreover, the evolved networks are generally more resilient than networks in which the random deletion of proteins was restricted to those disrupted in LTEE. These results suggest that evolution in the LTEE has favored PPI networks that are, on average, more resilient than expected from the genetic variation across the evolved strains. My findings therefore support the hypothesis that selection maintains protein interactome resilience over evolutionary time.

MPL resolves genetic linkage in fitness inference from complex evolutionary histories

Genetic linkage causes the fate of new mutations in a population to be contingent on the genetic background on which they appear. This makes it challenging to identify how individual mutations affect fitness. To overcome this challenge, we developed marginal path likelihood (MPL), a method to infer selection from evolutionary histories that resolves genetic linkage. Validation on real and simulated data sets shows that MPL is fast and accurate, outperforming existing inference approaches. We found that resolving linkage is crucial for accurately quantifying selection in complex evolving populations, which we demonstrate through a quantitative analysis of intrahost HIV-1 evolution using multiple patient data sets. Linkage effects generated by variants that sweep rapidly through the population are particularly strong, extending far across the genome. Taken together, our results argue for the importance of resolving linkage in studies of natural selection.

Global epistasis emerges from a generic model of a complex trait

Epistasis between mutations can make adaptation contingent on evolutionary history. Yet despite widespread ‘microscopic’ epistasis between the mutations involved, microbial evolution experiments show consistent patterns of fitness increase between replicate lines. Recent work shows that this consistency is driven in part by global patterns of diminishing-returns and increasing-costs epistasis, which make mutations systematically less beneficial (or more deleterious) on fitter genetic backgrounds. However, the origin of this ‘global’ epistasis remains unknown. Here, we show that diminishing-returns and increasing-costs epistasis emerge generically as a consequence of pervasive microscopic epistasis. Our model predicts a specific quantitative relationship between the magnitude of global epistasis and the stochastic effects of microscopic epistasis, which we confirm by reanalyzing existing data. We further show that the distribution of fitness effects takes on a universal form when epistasis is widespread and introduce a novel fitness landscape model to show how phenotypic evolution can be repeatable despite sequence-level stochasticity.

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.

Idiosyncratic variation in the fitness costs of tetracycline-resistance mutations in Escherichia coli

A bacterium's fitness relative to its competitors, both in the presence and absence of antibiotics, plays a key role in its ecological success and clinical impact. In this study, we examine whether tetracycline-resistant mutants are less fit in the absence of the drug than their sensitive parents, and whether the fitness cost of resistance is constant or variable across independently derived lines. Tetracycline-resistant lines suffered, on average, a reduction in fitness of almost 8%. There was substantial among-line variation in the fitness cost. This variation was not associated with the level of resistance conferred by the mutations, nor did it vary significantly across several genetic backgrounds. The two resistant lines with the most extreme fitness costs involved functionally unrelated mutations on different genetic backgrounds. However, there was also significant variation in the fitness costs for mutations affecting the same pathway and even different alleles of the same gene. Our findings demonstrate that the fitness costs of antibiotic resistance do not always correlate with the phenotypic level of resistance or the underlying genetic changes. Instead, these costs reflect the idiosyncratic effects of particular resistance mutations and the genetic backgrounds in which they occur.

A conserved Zn2Cys6 transcription factor, identified in a spontaneous mutant from in vitro passaging, is involved in pathogenicity of the blackleg fungus Leptosphaeria maculans

Continuous passaging in vitro can lead to the accumulation of changes in DNA sequence that potentially affect the properties of microbes, making them different from the original isolates. The identification of such genetic alterations is rare in fungi. A set of insertional mutants in the plant pathogenic fungus Leptosphaeria maculans, all derived from the same transformation experiment, had independent Agrobacterium T-DNA insertions and reduced pathogenicity on canola (Brassica napus). None of the insertions co-segregated in progeny from crosses with the reduction in pathogenicity. Genome sequences of three strains were analysed, and a mutation identified in a gene (ptf1, for pathogenicity-associated transcription factor 1) encoding a putative Zn₂(II)Cys₆ transcription factor. Homologs are found in other ascomycetes, and are required for pathogenicity by Fusarium graminearum, Fusarium oxysporum and Magnaporthe oryzae. The mutation in the L. maculans ptf1 gene co-segregates in progeny from crosses with the reduction in pathogenicity, a strain with an independent mutant allele isolated using CRISPR-Cas9 editing has reduced pathogenicity, and addition of wild type copies of the gene restores pathogenicity. Thus, this work defines a base pair substitution that occurred during in vitro passaging of a fungus that contributed to an attenuation of pathogenicity.

Beyond Standardization: Improving External Validity and Reproducibility in Experimental Evolution

Discussions of reproducibility are casting doubts on the credibility of experimental outcomes in the life sciences. Although experimental evolution is not typically included in these discussions, this field is also subject to low reproducibility, partly because of the inherent contingencies affecting the evolutionary process. A received view in experimental studies more generally is that standardization (i.e., rigorous homogenization of experimental conditions) is a solution to some issues of significance and internal validity. However, this solution hides several difficulties, including a reduction of external validity and reproducibility. After explaining the meaning of these two notions in the context of experimental evolution, we import from the fields of animal research and ecology and suggests that systematic heterogenization of experimental factors could prove a promising alternative. We also incorporate into our analysis some philosophical reflections on the nature and diversity of research objectives in experimental evolution.