Long-term Evolution Experiment Research Articles

Template switching during DNA replication is a prevalent source of adaptive gene amplification.

Copy number variants (CNVs)-gains and losses of genomic sequences-are an important source of genetic variation underlying rapid adaptation and genome evolution. However, despite their central role in evolution little is known about the factors that contribute to the structure, size, formation rate, and fitness effects of adaptive CNVs. Local genomic sequences are likely to be an important determinant of these properties. Whereas it is known that point mutation rates vary with genomic location and local DNA sequence features, the role of genome architecture in the formation, selection, and the resulting evolutionary dynamics of CNVs is poorly understood. Previously, we have found that the GAP1 gene in Saccharomyces cerevisiae undergoes frequent and repeated amplification and selection under long-term experimental evolution in glutamine-limiting conditions. The GAP1 gene has a unique genomic architecture consisting of two flanking long terminal repeats (LTRs) and a proximate origin of DNA replication (autonomously replicating sequence, ARS), which are likely to promote rapid GAP1 CNV formation. To test the role of these genomic elements on CNV-mediated adaptive evolution, we performed experimental evolution in glutamine-limited chemostats using engineered strains lacking either the adjacent LTRs, ARS, or all elements. Using a CNV reporter system and neural network simulation-based inference (nnSBI) we quantified the formation rate and fitness effect of CNVs for each strain. We find that although GAP1 CNVs repeatedly form and sweep to high frequency in strains with modified genome architecture, removal of local DNA elements significantly impacts the rate and fitness effect of CNVs and the rate of adaptation. We performed genome sequence analysis to define the molecular mechanisms of CNV formation for 177 CNV lineages. We find that across all four strain backgrounds, between 26% and 80% of all GAP1 CNVs are mediated by Origin Dependent Inverted Repeat Amplification (ODIRA) which results from template switching between the leading and lagging strand during DNA synthesis. In the absence of the local ARS, a distal ARS can mediate CNV formation via ODIRA. In the absence of local LTRs, homologous recombination mechanisms still mediate gene amplification following de novo insertion of retrotransposon elements at the locus. Our study demonstrates the remarkable plasticity of the genome and reveals that template switching during DNA replication is a frequent source of adaptive CNVs.

Use of Viscous medium to study anthelmintic drug action in Caenorhabditis elegans

Caenorhabditis elegans is an appealing tool for experimental evolution and for working with antiparasitic drugs, from understanding the molecular mechanisms of drug action and resistance to uncover new drug targets. We present a new methodology for studying the impact of antiparasitic drugs in C. elegans. Viscous medium was initially designed for C. elegans maintenance during long-term evolution experiments. Viscous medium provides a less structured environment than the standard nematode growth media agar, yet the bacteria food source remains suspended. Further, the Viscous medium offers the worm population enough support to move freely, mate, and reproduce at a rate comparable to standard agar cultures. Here, the Viscous medium was adapted for use in antiparasitic research. We observed a similar sensitivity of C. elegans to anthelmintic drugs as in standard liquid media and statistical difference to the standard agar media through a larval development assay. Using Viscous medium in C. elegans studies will considerably improve antiparasitic resistance research, and this medium could be used in studies aimed at understanding long-term multigenerational drug activity.

Promoter recruitment drives the emergence of proto-genes in a long-term evolution experiment with Escherichia coli.

The phenomenon of de novo gene birth-the emergence of genes from non-genic sequences-has received considerable attention due to the widespread occurrence of genes that are unique to particular species or genomes. Most instances of de novo gene birth have been recognized through comparative analyses of genome sequences in eukaryotes, despite the abundance of novel, lineage-specific genes in bacteria and the relative ease with which bacteria can be studied in an experimental context. Here, we explore the genetic record of the Escherichia coli long-term evolution experiment (LTEE) for changes indicative of "proto-genic" phases of new gene birth in which non-genic sequences evolve stable transcription and/or translation. Over the time span of the LTEE, non-genic regions are frequently transcribed, translated and differentially expressed, with levels of transcription across low-expressed regions increasing in later generations of the experiment. Proto-genes formed downstream of new mutations result either from insertion element activity or chromosomal translocations that fused preexisting regulatory sequences to regions that were not expressed in the LTEE ancestor. Additionally, we identified instances of proto-gene emergence in which a previously unexpressed sequence was transcribed after formation of an upstream promoter, although such cases were rare compared to those caused by recruitment of preexisting promoters. Tracing the origin of the causative mutations, we discovered that most occurred early in the history of the LTEE, often within the first 20,000 generations, and became fixed soon after emergence. Our findings show that proto-genes emerge frequently within evolving populations, can persist stably, and can serve as potential substrates for new gene formation.

Accumulation Dynamics of Defective Genomes during Experimental Evolution of Two Betacoronaviruses.

Virus-encoded replicases often generate aberrant RNA genomes, known as defective viral genomes (DVGs). When co-infected with a helper virus providing necessary proteins, DVGs can multiply and spread. While DVGs depend on the helper virus for propagation, they can in some cases disrupt infectious virus replication, impact immune responses, and affect viral persistence or evolution. Understanding the dynamics of DVGs alongside standard viral genomes during infection remains unclear. To address this, we conducted a long-term experimental evolution of two betacoronaviruses, the human coronavirus OC43 (HCoV-OC43) and the murine hepatitis virus (MHV), in cell culture at both high and low multiplicities of infection (MOI). We then performed RNA-seq at regular time intervals, reconstructed DVGs, and analyzed their accumulation dynamics. Our findings indicate that DVGs evolved to exhibit greater diversity and abundance, with deletions and insertions being the most common types. Notably, some high MOI deletions showed very limited temporary existence, while others became prevalent over time. We observed differences in DVG abundance between high and low MOI conditions in HCoV-OC43 samples. The size distribution of HCoV-OC43 genomes with deletions differed between high and low MOI passages. In low MOI lineages, short and long DVGs were the most common, with an additional cluster in high MOI lineages which became more prevalent along evolutionary time. MHV also showed variations in DVG size distribution at different MOI conditions, though they were less pronounced compared to HCoV-OC43, suggesting a more random distribution of DVG sizes. We identified hotspot regions for deletions that evolved at a high MOI, primarily within cistrons encoding structural and accessory proteins. In conclusion, our study illustrates the widespread formation of DVGs during betacoronavirus evolution, influenced by MOI and cell- and virus-specific factors.

“R” U ready?: a case study using R to analyze changes in gene expression during evolution

As high-throughput methods become more common, training undergraduates to analyze data must include having them generate informative summaries of large datasets. This flexible case study provides an opportunity for undergraduate students to become familiar with the capabilities of R programming in the context of high-throughput evolutionary data collected using macroarrays. The story line introduces a recent graduate hired at a biotech firm and tasked with analysis and visualization of changes in gene expression from 20,000 generations of the Lenski Lab’s Long-Term Evolution Experiment (LTEE). Our main character is not familiar with R and is guided by a coworker to learn about this platform. Initially this involves a step-by-step analysis of the small Iris dataset built into R which includes sepal and petal length of three species of irises. Practice calculating summary statistics and correlations, and making histograms and scatter plots, prepares the protagonist to perform similar analyses with the LTEE dataset. In the LTEE module, students analyze gene expression data from the long-term evolutionary experiments, developing their skills in manipulating and interpreting large scientific datasets through visualizations and statistical analysis. Prerequisite knowledge is basic statistics, the Central Dogma, and basic evolutionary principles. The Iris module provides hands-on experience using R programming to explore and visualize a simple dataset; it can be used independently as an introduction to R for biological data or skipped if students already have some experience with R. Both modules emphasize understanding the utility of R, rather than creation of original code. Pilot testing showed the case study was well-received by students and faculty, who described it as a clear introduction to R and appreciated the value of R for visualizing and analyzing large datasets.

Correlated Allele Frequency Changes Reveal Clonal Structure and Selection in Temporal Genetic Data.

In evolving populations where the rate of beneficial mutations is large, subpopulations of individuals with competing beneficial mutations can be maintained over long times. Evolution with this kind of clonal structure is commonly observed in a wide range of microbial and viral populations. However, it can be difficult to completely resolve clonal dynamics in data. This is due to limited read lengths in high-throughput sequencing methods, which are often insufficient to directly measure linkage disequilibrium or determine clonal structure. Here, we develop a method to infer clonal structure using correlated allele frequency changes in time-series sequence data. Simulations show that our method recovers true, underlying clonal structures when they are known and accurately estimate linkage disequilibrium. This information can then be combined with other inference methods to improve estimates of the fitness effects of individual mutations. Applications to data suggest novel clonal structures in an E. coli long-term evolution experiment, and yield improved predictions of the effects of mutations on bacterial fitness and antibiotic resistance. Moreover, our method is computationally efficient, requiring orders of magnitude less run time for large data sets than existing methods. Overall, our method provides a powerful tool to infer clonal structures from data sets where only allele frequencies are available, which can also improve downstream analyses.

Emergence and maintenance of stable coexistence during a long-term multicellular evolution experiment.

The evolution of multicellular life spurred evolutionary radiations, fundamentally changing many of Earth's ecosystems. Yet little is known about how early steps in the evolution of multicellularity affect eco-evolutionary dynamics. Through long-term experimental evolution, we observed niche partitioning and the adaptive divergence of two specialized lineages from a single multicellular ancestor. Over 715 daily transfers, snowflake yeast were subjected to selection for rapid growth, followed by selection favouring larger group size. Small and large cluster-forming lineages evolved from a monomorphic ancestor, coexisting for over ~4,300 generations, specializing on divergent aspects of a trade-off between growth rate and survival. Through modelling and experimentation, we demonstrate that coexistence is maintained by a trade-off between organismal size and competitiveness for dissolved oxygen. Taken together, this work shows how the evolution of a new level of biological individuality can rapidly drive adaptive diversification and the expansion of a nascent multicellular niche, one of the most historically impactful emergent properties of this evolutionary transition.

Proteostatic tuning underpins the evolution of novel multicellular traits.

The evolution of multicellularity paved the way for the origin of complex life on Earth, but little is known about the mechanistic basis of early multicellular evolution. Here, we examine the molecular basis of multicellular adaptation in the multicellularity long-term evolution experiment (MuLTEE). We demonstrate that cellular elongation, a key adaptation underpinning increased biophysical toughness and organismal size, is convergently driven by down-regulation of the chaperone Hsp90. Mechanistically, Hsp90-mediated morphogenesis operates by destabilizing the cyclin-dependent kinase Cdc28, resulting in delayed mitosis and prolonged polarized growth. Reinstatement of Hsp90 or Cdc28 expression resulted in shortened cells that formed smaller groups with reduced multicellular fitness. Together, our results show how ancient protein folding systems can be tuned to drive rapid evolution at a new level of biological individuality by revealing novel developmental phenotypes.

Rediversification following ecotype isolation reveals hidden adaptive potential

Microbial communities play a critical role in ecological processes, and their diversity is key to their functioning. However, little is known about whether communities can regenerate ecological diversity following ecotype removal or extinction and how the rediversified communities would compare to the original ones. Here, we show that simple two-ecotype communities from the E.coli long-term evolution experiment (LTEE) consistently rediversified into two ecotypes following the isolation of one of the ecotypes, coexisting via negative frequency-dependent selection. Communities separated by more than 30,000 generations of evolutionary time rediversify in similar ways. The rediversified ecotype appears to share a number of growth traits with the ecotype it replaces. However, the rediversified community is also different from the original community in ways relevant to the mechanism of ecotype coexistence-for example, in stationary phase response and survival. We found substantial variation in the transcriptional states between the two original ecotypes, whereas the differences within the rediversified community were comparatively smaller, although the rediversified community showed unique patterns of differential expression. Our results suggest that evolution may leave room for alternative diversification processes even in a maximally reduced community of only two strains. We hypothesize that the presence of alternative evolutionary pathways may be even more pronounced in communities of many species where there are even more potential niches, highlighting an important role for perturbations, such as species removal, in evolving ecological communities.

The growth of microbial cultures complies with the laws of thermodynamics

The paper is the translation of the previously proposed growth model in a thermodynamic balance of the Gibbs free energy of the system (medium + microbes), based on a simple scheme of the cell duplication. In each duplication step, the cells garner a small extra Gibbs energy from the surrounding medium that loses also some energy through an exothermic effect. It turns out that the each duplication step implies an increase of the entropy of the system, but a decrease of the entropy of the involved cells. The overall number of duplication steps therefore determines the energy balance of the whole growth process. The growth model implies a relationship that links this number with the maximum specific growth rate and the no-growth latency that precedes the growth onset, namely, two parameters that reflect the biological efficiency of the cells. For this reason, the overall number of duplication steps, determined according to this model, seems the best proxy of the fitness of the microbial culture. In a Long Term Evolution Experiment (LTEE), the increasing fitness would therefore correspond to larger growth extent and specific rate, as well as to shorter pre-growth latency. This suggests that the gain of Gibbs free energy accumulated through the LTEE several-thousand generations leads to a faster attainment of the eventual steady state of the growth and a faster increase of the entropy of the system. If applied to a continuous LTEE carried out with a chemostat, this trend should reveal that the evolution of the culture (medium + cells) is an irreversible process.

Linking genotypic and phenotypic changes in the E. coli long-term evolution experiment using metabolomics

Changes in an organism’s environment, genome, or gene expression patterns can lead to changes in its metabolism. The metabolic phenotype can be under selection and contributes to adaptation. However, the networked and convoluted nature of an organism’s metabolism makes relating mutations, metabolic changes, and effects on fitness challenging. To overcome this challenge, we use the long-term evolution experiment (LTEE) with E. coli as a model to understand how mutations can eventually affect metabolism and perhaps fitness. We used mass spectrometry to broadly survey the metabolomes of the ancestral strains and all 12 evolved lines. We combined this metabolic data with mutation and expression data to suggest how mutations that alter specific reaction pathways, such as the biosynthesis of nicotinamide adenine dinucleotide, might increase fitness in the system. Our work provides a better understanding of how mutations might affect fitness through the metabolic changes in the LTEE and thus provides a major step in developing a complete genotype–phenotype map for this experimental system.

Linking genotypic and phenotypic changes in the E. coli long-term evolution experiment using metabolomics.

Changes in an organism's environment, genome, or gene expression patterns can lead to changes in its metabolism. The metabolic phenotype can be under selection and contributes to adaptation. However, the networked and convoluted nature of an organism's metabolism makes relating mutations, metabolic changes, and effects on fitness challenging. To overcome this challenge, we use the long-term evolution experiment (LTEE) with E. coli as a model to understand how mutations can eventually affect metabolism and perhaps fitness. We used mass spectrometry to broadly survey the metabolomes of the ancestral strains and all 12 evolved lines. We combined this metabolic data with mutation and expression data to suggest how mutations that alter specific reaction pathways, such as the biosynthesis of nicotinamide adenine dinucleotide, might increase fitness in the system. Our work provides a better understanding of how mutations might affect fitness through the metabolic changes in the LTEE and thus provides a major step in developing a complete genotype-phenotype map for this experimental system.

Heat-induced female biased sex ratio during development is not mitigated after prolonged thermal selection

BackgroundThe negative impacts of climate change on biodiversity are consistently increasing. Developmental stages are particularly sensitive in many ectotherms. Moreover, sex-specific differences in how organisms cope with thermal stress can produce biased sex ratios upon emergence, with potentially major impacts on population persistence. This is an issue that needs investigation, particularly testing whether thermal selection can alleviate sex ratio distortions in the long-term is a critical but neglected issue. Here, we report an experiment analyzing the sex ratio patterns at different developmental temperatures in Drosophila subobscura populations subjected to long-term experimental evolution (~ 30 generations) under a warming environment.ResultsWe show that exposure to high developmental temperatures consistently promotes sex ratio imbalance upon emergence, with a higher number of female than male offspring. Furthermore, we found that thermal selection resulting from evolution in a warming environment did not alleviate such sex ratio distortions generated by heat stress.ConclusionsWe demonstrate that heat stress during development can lead to clear sex ratio deviations upon emergence likely because of differential survival between sexes. In face of these findings, it is likely that sex ratio deviations of this sort occur in natural populations when facing environmental perturbation. The inability of many insects to avoid thermal shifts during their (more) sessile developmental stages makes this finding particularly troublesome for population subsistence in face of climate warming events.

Comparative genomics hints at dispensability of multiple essential genes in two Escherichia coli L-form strains.

Despite the critical role of bacterial cell walls in maintaining cell shapes, certain environmental stressors can induce the transition of many bacterial species into a wall-deficient state called L-form. Long-term induced Escherichia coli L-forms lose their rod shape and usually hold significant mutations that affect cell division and growth. Besides this, the genetic background of L-form bacteria is still poorly understood. In the present study, the genomes of two stable L-form strains of E. coli (NC-7 and LWF+) were sequenced and their gene mutation status was determined and compared with their parental strains. Comparative genomic analysis between two L-forms reveals both unique adaptions and common mutated genes, many of which belong to essential gene categories not involved in cell wall biosynthesis, indicating that L-form genetic adaptation impacts crucial metabolic pathways. Missense variants from L-forms and Lenski's long-term evolution experiment (LTEE) were analyzed in parallel using an optimized DeepSequence pipeline to investigate predicted mutation effects (α) on protein functions. We report that the two L-form strains analyzed display a frequency of 6-10% (0% for LTEE) in mutated essential genes where the missense variants have substantial impact on protein functions (α<0.5). This indicates the emergence of different survival strategies in L-forms through changes in essential genes during adaptions to cell wall deficiency. Collectively, our results shed light on the detailed genetic background of two E. coli L-forms and pave the way for further investigations of the gene functions in L-form bacterial models.

Open to Evolve Embodied Intelligence

For the goal of automatically evolving Embodied Intelligence (EI), we investigate an open software architecture inspired by the high surface area to volume ratio of animal lungs, which aims to avoid information theoretic limits on long term evolution experiments (LTEE) encountered with monolithic genetic programming trees. Instead individuals are teams composed of 1023 trees whose inputs and outputs are linked by a low entropy loss branching data (air) pathway. Most trees are shallow and software engineering’s failed disruption propagation (FDP) is observed in the small fraction of deep trees. After initial search, most improvements are at intermediate depths and performance is still rising even after 100 000 generations.Despite the use of double precision for the bifurcating data interconnect, some information loss is seen, particularly in early generations. The static optimisation benchmark, appears to encourage early convergence, which locks the population into possibly sub-optimal phenotypes. Later thousands of small improvements, sometimes in large bloated ensemble members, appear to compensate for early overfitting.Using tournament fitness selection and subtree crossover, we target pure nested side-effect free floating point functions, which are known to have low FDP, and high fidelity data paths, in the hope of generating code which is not too robust so as to prevent on going improvement. However, we again find genetic changes deep within trees are silent. For single precision, we find a maximum evolvability sweet spot with trees of depth 10 to 100. Accordingly, we suggest to evolve very large very complex programs needed for Embodied Intelligence, an open structure with a high surface area permitting most mutation sites to be within 10-100 levels of the organism’s environment, and many better placed test oracles to monitor the impact of mutated code, will be needed.

Robotic workflows for automated long-term adaptive laboratory evolution: improving ethanol utilization by Corynebacterium glutamicum

BackgroundAdaptive laboratory evolution (ALE) is known as a powerful tool for untargeted engineering of microbial strains and genomics research. It is particularly well suited for the adaptation of microorganisms to new environmental conditions, such as alternative substrate sources. Since the probability of generating beneficial mutations increases with the frequency of DNA replication, ALE experiments are ideally free of constraints on the required duration of cell proliferation.ResultsHere, we present an extended robotic workflow for performing long-term evolution experiments based on fully automated repetitive batch cultures (rbALE) in a well-controlled microbioreactor environment. Using a microtiter plate recycling approach, the number of batches and thus cell generations is technically unlimited. By applying the validated workflow in three parallel rbALE runs, ethanol utilization by Corynebacterium glutamicum ATCC 13032 (WT) was significantly improved. The evolved mutant strain WT_EtOH-Evo showed a specific ethanol uptake rate of 8.45 ± 0.12 mmolEtOH gCDW−1 h−1 and a growth rate of 0.15 ± 0.01 h−1 in lab-scale bioreactors. Genome sequencing of this strain revealed a striking single nucleotide variation (SNV) upstream of the ald gene (NCgl2698, cg3096) encoding acetaldehyde dehydrogenase (ALDH). The mutated basepair was previously predicted to be part of the binding site for the global transcriptional regulator GlxR, and re-engineering demonstrated that the identified SNV is key for enhanced ethanol assimilation. Decreased binding of GlxR leads to increased synthesis of the rate-limiting enzyme ALDH, which was confirmed by proteomics measurements.ConclusionsThe established rbALE technology is generally applicable to any microbial strain and selection pressure that fits the small-scale cultivation format. In addition, our specific results will enable improved production processes with C. glutamicum from ethanol, which is of particular interest for acetyl-CoA-derived products.

Evolution of a cross-feeding interaction following a key innovation in a long-term evolution experiment with Escherichia coli.

The evolution of a novel trait can profoundly change an organism's effects on its environment, which can in turn affect the further evolution of that organism and any coexisting organisms. We examine these effects and feedbacks following the evolution of a novel function in the Long-Term Evolution Experiment (LTEE) with Escherichia coli. A characteristic feature of E. coli is its inability to grow aerobically on citrate (Cit-). Nonetheless, a Cit+ variant with this capacity evolved in one LTEE population after 31 000 generations. The Cit+ clade then coexisted stably with another clade that retained the ancestral Cit- phenotype. This coexistence was shaped by the evolution of a cross-feeding relationship based on C4-dicarboxylic acids, particularly succinate, fumarate, and malate, that the Cit+ variants release into the medium. Both the Cit- and Cit+ cells evolved to grow on these excreted resources. The evolution of aerobic growth on citrate thus led to a transition from an ecosystem based on a single limiting resource, glucose, to one with at least five resources that were either shared or partitioned between the two coexisting clades. Our findings show that evolutionary novelties can change environmental conditions in ways that facilitate diversity by altering ecosystem structure and the evolutionary trajectories of coexisting lineages.

Daily Transfers, Archiving Populations, and Measuring Fitness in the Long-Term Evolution Experiment with Escherichia coli.

The Long-Term Evolution Experiment (LTEE) has followed twelve populations of Escherichia coli as they have adapted to a simple laboratory environment for more than 35 years and 77,000 bacterial generations. The setup and procedures used in the LTEE epitomize reliable and reproducible methods for studying microbial evolution. In this protocol, we first describe how the LTEE populations are transferred to fresh medium and cultured each day. Then, we describe how the LTEE populations are regularly checked for possible signs of contamination and archived to provide a permanent frozen "fossil record" for later study. Multiple safeguards included in these procedures are designed to prevent contamination, detect various problems when they occur, and recover from disruptions without appreciably setting back the progress of the experiment. One way that the overall tempo and character of evolutionary changes are monitored in the LTEE is by measuring the competitive fitness of populations and strains from the experiment. We describe how co-culture competition assays are conducted and provide both a spreadsheet and an R package (fitnessR) for calculating relative fitness from the results. Over the course of the LTEE, the behaviors of some populations have changed in interesting ways, and new technologies like whole-genome sequencing have provided additional avenues for investigating how the populations have evolved. We end by discussing how the original LTEE procedures have been updated to accommodate or take advantage of these changes. This protocol will be useful for researchers who use the LTEE as a model system for studying connections between evolution and genetics, molecular biology, systems biology, and ecology. More broadly, the LTEE provides a tried-and-true template for those who are beginning their own evolution experiments with new microbes, environments, and questions.

Dynamics of Podospora anserina Genome Evolution in a Long-Term Experiment.

The Podospora anserina long-term evolution experiment (PaLTEE) is the only running filamentous fungus study, which is still going on. The aim of our work is to trace the evolutionary dynamics of the accumulation of mutations in the genomes of eight haploid populations of P. anserina. The results of the genome-wide analysis of all of the lineages, performed 8 years after the start of the PaLTEE, are presented. Data analysis detected 312 single nucleotide polymorphisms (SNPs) and 39 short insertion-deletion mutations (indels) in total. There was a clear trend towards a linear increase in the number of SNPs depending on the experiment duration. Among 312 SNPs, 153 were fixed in the coding regions of P. anserina genome. Relatively few synonymous mutations were found, exactly 38; 42 were classified as nonsense mutations; 72 were assigned to missense mutations. In addition, 21 out of 39 indels identified were also localized in coding regions. Here, we also report the detection of parallel evolution at the paralog level in the P. anserina model system. Parallelism in evolution at the level of protein functions also occurs. The latter is especially true for various transcription factors, which may indicate selection leading to optimization of the wide range of cellular processes under experimental conditions.

Long-term evolution of proliferating cells using the eVOLVER platform

Experimental evolution using fast-growing unicellular organisms is a unique strategy for deciphering the principles and mechanisms underlying evolutionary processes as well as the architecture and wiring of basic biological functions. Over the past decade, this approach has benefited from the development of powerful systems for the continuous control of the growth of independently evolving cultures. While the first devices compatible with multiplexed experimental evolution remained challenging to implement and required constant user intervention, the recently developed eVOLVER framework represents a fully automated closed-loop system for laboratory evolution assays. However, it remained difficult to maintain and compare parallel evolving cultures in tightly controlled environments over long periods of time using eVOLVER. Furthermore, a number of tools were lacking to cope with the various issues that inevitably occur when conducting such long-term assays. Here we present a significant upgrade of the eVOLVER framework, providing major modifications of the experimental methodology, hardware and software as well as a new stand-alone protocol. Altogether, these adaptations and improvements make the eVOLVER a versatile and unparalleled set-up for long-term experimental evolution.