Digital genetics: unravelling the genetic basis of evolution

Reducing Epistasis and Pleiotropy Can Avoid the Survival of the Flattest Tragedy

This study investigates whether reducing epistasis and pleiotropy enhances mutational robustness in evolutionary adaptation, utilizing an indirect encoded model within the "survival of the flattest" (SoF) fitness landscape. By simulating genetic variations and their phenotypic consequences, we explore organisms' adaptive mechanisms to maintain positions on higher, narrower evolutionary peaks amidst environmental and genetic pressures. Our results reveal that organisms can indeed sustain their advantageous positions by minimizing the complexity of genetic interactions-specifically, by reducing the levels of epistasis and pleiotropy. This finding suggests a counterintuitive strategy for evolutionary stability: simpler genetic architectures, characterized by fewer gene interactions and multifunctional genes, confer a survival advantage by enhancing mutational robustness. This study contributes to our understanding of the genetic underpinnings of adaptability and robustness, challenging traditional views that equate complexity with fitness in dynamic environments.

Self-replicating artificial neural networks give rise to universal evolutionary dynamics

In evolutionary models, mutations are exogenously introduced by the modeler, rather than endogenously introduced by the replicator itself. We present a new deep-learning based computational model, the self-replicating artificial neural network (SeRANN). We train it to (i) copy its own genotype, like a biological organism, which introduces endogenous spontaneous mutations; and (ii) simultaneously perform a classification task that determines its fertility. Evolving 1,000 SeRANNs for 6,000 generations, we observed various evolutionary phenomena such as adaptation, clonal interference, epistasis, and evolution of both the mutation rate and the distribution of fitness effects of new mutations. Our results demonstrate that universal evolutionary phenomena can naturally emerge in a self-replicator model when both selection and mutation are implicit and endogenous. We therefore suggest that SeRANN can be applied to explore and test various evolutionary dynamics and hypotheses.

How Turing parasites expand the computational landscape of digital life

Why are living systems complex? Why does the biosphere contain living beings with complexity features beyond those of the simplest replicators? What kind of evolutionary pressures result in more complex life forms? These are key questions that pervade the problem of how complexity arises in evolution. One particular way of tackling this is grounded in an algorithmic description of life: living organisms can be seen as systems that extract and process information from their surroundings to reduce uncertainty. Here we take this computational approach using a simple bit string model of coevolving agents and their parasites. While agents try to predict their worlds, parasites do the same with their hosts. The result of this process is that, to escape their parasites, the host agents expand their computational complexity despite the cost of maintaining it. This, in turn, is followed by increasingly complex parasitic counterparts. Such arms races display several qualitative phases, from monotonous to punctuated evolution or even ecological collapse. Our minimal model illustrates the relevance of parasites in providing an active mechanism for expanding living complexity beyond simple replicators, suggesting that parasitic agents are likely to be a major evolutionary driver for biological complexity.

avidaR: an R library to perform complex queries on an ontology-based database of digital organisms

Digital evolution is a branch of artificial life in which self-replicating computer programs-digital organisms-mutate and evolve within a user-defined computational environment. In spite of its value in biology, we still lack an up-to-date and comprehensive database on digital organisms resulting from evolution experiments. Therefore, we have developed an ontology-based semantic database-avidaDB-and an R package-avidaR-that provides users of the R programming language with an easy-to-use tool for performing complex queries without specific knowledge of SPARQL or RDF. avidaR can be used to do research on robustness, evolvability, complexity, phenotypic plasticity, gene regulatory networks, and genomic architecture by retrieving the genomes, phenotypes, and transcriptomes of more than a million digital organisms available on avidaDB. avidaR is already accepted on CRAN (i.e., a comprehensive collection of R packages contributed by the R community) and will make biologists better equipped to embrace the field of digital evolution.

Ontology for the Avida digital evolution platform

The Ontology for Avida (OntoAvida) aims to develop an integrated vocabulary for the description of Avida, the most widely used computational approach for performing experimental evolution using digital organisms-self-replicating computer programs that evolve within a user-defined computational environment. The lack of a clearly defined vocabulary makes some biologists feel reluctant to embrace the field of digital evolution. This integrated framework empowers biologists by equipping them with the necessary tools to explore and analyze the field of digital evolution more effectively. By leveraging the vocabulary of Avida, researchers can gain deeper insights into the evolutionary processes and dynamics of digital organisms. In addition, OntoAvida allows researchers to make inference based on certain rules and constraints, facilitate the reproducibility of in silico evolution experiments and trace the provenance of the data stored in avidaDB-an RDF database containing the genomes, transcriptomes, and phenotypes of more than a million digital organisms. OntoAvida is part of the Open Biological and Biomedical Ontologies (OBO Foundry) and is available at http://www.obofoundry.org/ontology/ontoavida.html .

Emergence of functional information from multivariate correlations

The information content of symbolic sequences (such as nucleic or amino acid sequences, but also neuronal firings or strings of letters) can be calculated from an ensemble of such sequences, but because information cannot be assigned to single sequences, we cannot correlate information to other observables attached to the sequence. Here we show that an information score obtained from multivariate (multiple-variable) correlations within sequences of a 'training' ensemble can be used to predict observables of out-of-sample sequences with an accuracy that scales with the complexity of correlations, showing that functional information emerges from a hierarchy of multi-variable correlations. This article is part of the theme issue 'Emergent phenomena in complex physical and socio-technical systems: from cells to societies'.

A Study of the Coevolution of Digital Organisms with an Evolutionary Cellular Automaton

This paper presents an Evolutionary Cellular Automaton (ECA) that simulates the evolutionary dynamics of biological interactions by manipulating strategies of dispersion and associations between digital organisms. The parameterization of the different types of interaction and distribution strategies using configuration files generates easily interpretable results. In that respect, ECA is an effective instrument for measuring the effects of relative adaptive advantages and a good resource for studying natural selection. Although ECA works effectively in obtaining the expected results from most well-known biological interactions, some unexpected effects were observed. For example, organisms uniformly distributed in fragmented habitats do not favor eusociality, and mutualism evolved from parasitism simply by varying phenotypic flexibility. Finally, we have verified that natural selection represents a cost for the emergence of sex by destabilizing the stable evolutionary strategy of the 1:1 sex ratio after generating randomly different distributions in each generation.

The Surprising Creativity of Digital Evolution: A Collection of Anecdotes from the Evolutionary Computation and Artificial Life Research Communities

Evolution provides a creative fount of complex and subtle adaptations that often surprise the scientists who discover them. However, the creativity of evolution is not limited to the natural world: Artificial organisms evolving in computational environments have also elicited surprise and wonder from the researchers studying them. The process of evolution is an algorithmic process that transcends the substrate in which it occurs. Indeed, many researchers in the field of digital evolution can provide examples of how their evolving algorithms and organisms have creatively subverted their expectations or intentions, exposed unrecognized bugs in their code, produced unexpectedly adaptations, or engaged in behaviors and outcomes, uncannily convergent with ones found in nature. Such stories routinely reveal surprise and creativity by evolution in these digital worlds, but they rarely fit into the standard scientific narrative. Instead they are often treated as mere obstacles to be overcome, rather than results that warrant study in their own right. Bugs are fixed, experiments are refocused, and one-off surprises are collapsed into a single data point. The stories themselves are traded among researchers through oral tradition, but that mode of information transmission is inefficient and prone to error and outright loss. Moreover, the fact that these stories tend to be shared only among practitioners means that many natural scientists do not realize how interesting and lifelike digital organisms are and how natural their evolution can be. To our knowledge, no collection of such anecdotes has been published before. This article is the crowd-sourced product of researchers in the fields of artificial life and evolutionary computation who have provided first-hand accounts of such cases. It thus serves as a written, fact-checked collection of scientifically important and even entertaining stories. In doing so we also present here substantial evidence that the existence and importance of evolutionary surprises extends beyond the natural world, and may indeed be a universal property of all complex evolving systems.

Coevolution of Genome Architecture and Social Behavior

Although social behavior can have a strong genetic component, it can also result in selection on genome structure and function, thereby influencing the evolution of the genome itself. Here we explore the bidirectional links between social behavior and genome architecture by considering variation in social and/or mating behavior among populations (social polymorphisms) and across closely related species. We propose that social behavior can influence genome architecture via associated demographic changes due to social living. We establish guidelines to exploit emerging whole-genome sequences using analytical approaches that examine genome structure and function at different levels (regulatory vs structural variation) from the perspective of both molecular biology and population genetics in an ecological context.

Examining Community Stability in the Face of Mass Extinction in Communities of Digital Organisms

Digital evolution is a computer-based instantiation of Darwinian evolution in which short self-replicating computer programs compete, mutate, and evolve. It is an excellent platform for addressing topics in long-term evolution and paleobiology, such as mass extinction and recovery, with experimental evolutionary approaches. We evolved model communities with ecological interdependence among community members, which were subjected to two principal types of mass extinction: a pulse extinction that killed randomly, and a selective press extinction involving an alteration of the abiotic environment to which the communities had to adapt. These treatments were applied at two different strengths, along with unperturbed control experiments. We examined how stability in the digital communities was affected from the perspectives of division of labor, relative shift in rank abundance, and genealogical connectedness of the community's component ecotypes. Mass extinction that was due to a Strong Press treatment was most effective in producing reshaped communities that differed from the pre-treatment ones in all of the measured perspectives; weaker versions of the treatments did not generally produce significant departures from a Control treatment; and results for the Strong Pulse treatment generally fell between those extremes. The Strong Pulse treatment differed from others in that it produced a slight but detectable shift towards more generalized communities. Compared to Press treatments, Pulse treatments also showed a greater contribution from re-evolved ecological doppelgangers rather than new ecotypes. However, relatively few Control communities showed stability in any of these metrics over the whole course of the experiment, and most did not represent stable states (by some measure of stability) that were disrupted by the extinction treatments. Our results have interesting, broad qualitative parallels with findings from the paleontological record, and show the potential of digital evolution studies to illuminate many aspects of mass extinction and recovery by addressing them in a truly experimental manner.

Origin of life in a digital microcosm

While all organisms on Earth share a common descent, there is no consensus on whether the origin of the ancestral self-replicator was a one-off event or whether it only represented the final survivor of multiple origins. Here, we use the digital evolution system Avida to study the origin of self-replicating computer programs. By using a computational system, we avoid many of the uncertainties inherent in any biochemical system of self-replicators (while running the risk of ignoring a fundamental aspect of biochemistry). We generated the exhaustive set of minimal-genome self-replicators and analysed the network structure of this fitness landscape. We further examined the evolvability of these self-replicators and found that the evolvability of a self-replicator is dependent on its genomic architecture. We also studied the differential ability of replicators to take over the population when competed against each other, akin to a primordial-soup model of biogenesis, and found that the probability of a self-replicator outcompeting the others is not uniform. Instead, progenitor (most-recent common ancestor) genotypes are clustered in a small region of the replicator space. Our results demonstrate how computational systems can be used as test systems for hypotheses concerning the origin of life.This article is part of the themed issue 'Reconceptualizing the origins of life'.

Insights into resource consumption, cross-feeding, system collapse, stability and biodiversity from an artificial ecosystem

Community ecosystems at very different levels of biological organization often have similar properties. Coexistence of multiple species, cross-feeding, biodiversity and fluctuating population dynamics are just a few of the properties that arise in a range of ecological settings. Here we develop a bottom-up model of consumer-resource interactions, in the form of an artificial ecosystem 'number soup', which reflects basic properties of many bacterial and other community ecologies. We demonstrate four key properties of the number soup model: (i) communities self-organize so that all available resources are fully consumed; (ii) reciprocal cross-feeding is a common evolutionary outcome, which evolves in a number of stages, and many transitional species are involved; (iii) the evolved ecosystems are often 'robust yet fragile', with keystone species required to prevent the whole system from collapsing; (iv) non-equilibrium dynamics and chaotic patterns are general properties, readily generating rich biodiversity. These properties have been observed in empirical ecosystems, ranging from bacteria to rainforests. Establishing similar properties in an evolutionary model as simple as the number soup suggests that these four properties are ubiquitous features of all community ecosystems, and raises questions about how we interpret ecosystem structure in the context of natural selection.

The genotype-phenotype map of an evolving digital organism

To understand how evolving systems bring forth novel and useful phenotypes, it is essential to understand the relationship between genotypic and phenotypic change. Artificial evolving systems can help us understand whether the genotype-phenotype maps of natural evolving systems are highly unusual, and it may help create evolvable artificial systems. Here we characterize the genotype-phenotype map of digital organisms in Avida, a platform for digital evolution. We consider digital organisms from a vast space of 10¹⁴¹ genotypes (instruction sequences), which can form 512 different phenotypes. These phenotypes are distinguished by different Boolean logic functions they can compute, as well as by the complexity of these functions. We observe several properties with parallels in natural systems, such as connected genotype networks and asymmetric phenotypic transitions. The likely common cause is robustness to genotypic change. We describe an intriguing tension between phenotypic complexity and evolvability that may have implications for biological evolution. On the one hand, genotypic change is more likely to yield novel phenotypes in more complex organisms. On the other hand, the total number of novel phenotypes reachable through genotypic change is highest for organisms with simple phenotypes. Artificial evolving systems can help us study aspects of biological evolvability that are not accessible in vastly more complex natural systems. They can also help identify properties, such as robustness, that are required for both human-designed artificial systems and synthetic biological systems to be evolvable.

The phenotype of an organism comprises the set of morphological and functional traits encoded by its genome. In natural evolving systems, phenotypes are organized into mutationally connected networks of genotypes, which increase the likelihood for an evolving population to encounter novel adaptive phenotypes (i.e., its evolvability). We do not know whether artificial systems, such as self-replicating and evolving computer programs—digital organisms—are more or less evolvable than natural systems. By studying how genotypes map onto phenotypes in digital organisms, we characterize many commonalities between natural and artificial evolving systems. In addition, we show that phenotypic complexity can both facilitate and constrain evolution, which harbors lessons not only for designing evolvable artificial systems, but also for synthetic biology.

Different Evolutionary Paths to Complexity for Small and Large Populations of Digital Organisms

A major aim of evolutionary biology is to explain the respective roles of adaptive versus non-adaptive changes in the evolution of complexity. While selection is certainly responsible for the spread and maintenance of complex phenotypes, this does not automatically imply that strong selection enhances the chance for the emergence of novel traits, that is, the origination of complexity. Population size is one parameter that alters the relative importance of adaptive and non-adaptive processes: as population size decreases, selection weakens and genetic drift grows in importance. Because of this relationship, many theories invoke a role for population size in the evolution of complexity. Such theories are difficult to test empirically because of the time required for the evolution of complexity in biological populations. Here, we used digital experimental evolution to test whether large or small asexual populations tend to evolve greater complexity. We find that both small and large-but not intermediate-sized-populations are favored to evolve larger genomes, which provides the opportunity for subsequent increases in phenotypic complexity. However, small and large populations followed different evolutionary paths towards these novel traits. Small populations evolved larger genomes by fixing slightly deleterious insertions, while large populations fixed rare beneficial insertions that increased genome size. These results demonstrate that genetic drift can lead to the evolution of complexity in small populations and that purifying selection is not powerful enough to prevent the evolution of complexity in large populations.

Experimental Evolution as an Underutilized Tool for Studying Beneficial Animal-Microbe Interactions

Microorganisms play a significant role in the evolution and functioning of the eukaryotes with which they interact. Much of our understanding of beneficial host–microbe interactions stems from studying already established associations; we often infer the genotypic and environmental conditions that led to the existing host–microbe relationships. However, several outstanding questions remain, including understanding how host and microbial (internal) traits, and ecological and evolutionary (external) processes, influence the origin of beneficial host–microbe associations. Experimental evolution has helped address a range of evolutionary and ecological questions across different model systems; however, it has been greatly underutilized as a tool to study beneficial host–microbe associations. In this review, we suggest ways in which experimental evolution can further our understanding of the proximate and ultimate mechanisms shaping mutualistic interactions between eukaryotic hosts and microbes. By tracking beneficial interactions under defined conditions or evolving novel associations among hosts and microbes with little prior evolutionary interaction, we can link specific genotypes to phenotypes that can be directly measured. Moreover, this approach will help address existing puzzles in beneficial symbiosis research: how symbioses evolve, how symbioses are maintained, and how both host and microbe influence their partner’s evolutionary trajectories. By bridging theoretical predictions and empirical tests, experimental evolution provides us with another approach to test hypotheses regarding the evolution of beneficial host–microbe associations.

Evolution of Genome Size in Asexual Digital Organisms

Genome sizes have evolved to vary widely, from 250 bases in viroids to 670 billion bases in some amoebas. This remarkable variation in genome size is the outcome of complex interactions between various evolutionary factors such as mutation rate and population size. While comparative genomics has uncovered how some of these evolutionary factors influence genome size, we still do not understand what drives genome size evolution. Specifically, it is not clear how the primordial mutational processes of base substitutions, insertions, and deletions influence genome size evolution in asexual organisms. Here, we use digital evolution to investigate genome size evolution by tracking genome edits and their fitness effects in real time. In agreement with empirical data, we find that mutation rate is inversely correlated with genome size in asexual populations. We show that at low point mutation rate, insertions are significantly more beneficial than deletions, driving genome expansion and the acquisition of phenotypic complexity. Conversely, the high mutational load experienced at high mutation rates inhibits genome growth, forcing the genomes to compress their genetic information. Our analyses suggest that the inverse relationship between mutation rate and genome size is a result of the tradeoff between evolving phenotypic innovation and limiting the mutational load.

Information-Theoretic Considerations Concerning the Origin of Life

Research investigating the origins of life usually either focuses on exploring possible life-bearing chemistries in the pre-biotic Earth, or else on synthetic approaches. Comparatively little work has explored fundamental issues concerning the spontaneous emergence of life using only concepts (such as information and evolution) that are divorced from any particular chemistry. Here, I advocate studying the probability of spontaneous molecular self-replication as a function of the information contained in the replicator, and the environmental conditions that might enable this emergence. I show (under certain simplifying assumptions) that the probability to discover a self-replicator by chance depends exponentially on the relative rate of formation of the monomers. If the rate at which monomers are formed is somewhat similar to the rate at which they would occur in a self-replicating polymer, the likelihood to discover such a replicator by chance is increased by many orders of magnitude. I document such an increase in searches for a self-replicator within the digital life system avida.

A Model of Genome Size Evolution for Prokaryotes in Stable and Fluctuating Environments

Temporal variability in ecosystems significantly impacts species diversity and ecosystem productivity and therefore the evolution of organisms. Different levels of environmental perturbations such as seasonal fluctuations, natural disasters, and global change have different impacts on organisms and therefore their ability to acclimatize and adapt. Thus, to understand how organisms evolve under different perturbations is a key for predicting how environmental change will impact species diversity and ecosystem productivity. Here, we developed a computer simulation utilizing the individual-based model approach to investigate genome size evolution of a haploid, clonal and free-living prokaryotic population across different levels of environmental perturbations. Our results show that a greater variability of the environment resulted in genomes with a larger number of genes. Environmental perturbations were more effectively buffered by populations of individuals with relatively large genomes. Unpredictable changes of the environment led to a series of population bottlenecks followed by adaptive radiations. Our model shows that the evolution of genome size is indirectly driven by the temporal variability of the environment. This complements the effects of natural selection directly acting on genome optimization. Furthermore, species that have evolved in relatively stable environments may face the greatest risk of extinction under global change as genome streamlining genetically constrains their ability to acclimatize to the new environmental conditions, unless mechanisms of genetic diversification such as horizontal gene transfer will enrich their gene pool and therefore their potential to adapt.

Shape matters: lifecycle of cooperative patches promotes cooperation in bulky populations

Natural cooperative systems take many forms, ranging from one-dimensional cyanobacteria arrays to fractal-like biofilms. We use in silico experimental systems to study a previously overlooked factor in the evolution of cooperation, physical shape of the population. We compare the emergence and maintenance of cooperation in populations of digital organisms that inhabit bulky (100 × 100 cells) or slender (4 × 2500) toroidal grids. Although more isolated subpopulations of secretors in a slender population could be expected to favor cooperation, we find the opposite: secretion evolves to higher levels in bulky populations. We identify the mechanistic explanation for the shape effect by analyzing the lifecycle and dynamics of cooperator patches, from their emergence and growth, to invasion by noncooperators and extinction. Because they are constrained by the population shape, the cooperator patches expand less in slender than in bulky populations, leading to fewer cooperators, less public good secretion, and generally lower cooperation. The patch dynamics and mechanisms of shape effect are robust across several digital cooperation systems and independent of the underlying basis for cooperation (public good secretion or a cooperation game). Our results urge for a greater consideration of population shape in the study of the evolution of cooperation across experimental and modeling systems.

Resistant traits in digital organisms do not revert preselection status despite extended deselection: implications to microbial antibiotics resistance

Antibiotics resistance is a serious biomedical issue as formally susceptible organisms gain resistance under its selective pressure. There have been contradictory results regarding the prevalence of resistance following withdrawal and disuse of the specific antibiotics. Here, we use experimental evolution in “digital organisms” to examine the rate of gain and loss of resistance under the assumption that there is no fitness cost for maintaining resistance. Our results show that selective pressure is likely to result in maximum resistance with respect to the selective pressure. During deselection as a result of disuse of the specific antibiotics, a large initial loss and prolonged stabilization of resistance are observed, but resistance is not lost to the stage of preselection. This suggests that a pool of partial persists organisms persist long after withdrawal of selective pressure at a relatively constant proportion. Hence, contradictory results regarding the prevalence of resistance following withdrawal and disuse of the specific antibiotics may be a statistical variation about constant proportion. Our results also show that subsequent reintroduction of the same selective pressure results in rapid regain of maximal resistance. Thus, our simulation results suggest that complete elimination of specific antibiotics resistance is unlikely after the disuse of antibiotics once a resistant pool of microorganisms has been established.

Artificial evolution by viability rather than competition

Evolutionary algorithms are widespread heuristic methods inspired by natural evolution to solve difficult problems for which analytical approaches are not suitable. In many domains experimenters are not only interested in discovering optimal solutions, but also in finding the largest number of different solutions satisfying minimal requirements. However, the formulation of an effective performance measure describing these requirements, also known as fitness function, represents a major challenge. The difficulty of combining and weighting multiple problem objectives and constraints of possibly varying nature and scale into a single fitness function often leads to unsatisfactory solutions. Furthermore, selective reproduction of the fittest solutions, which is inspired by competition-based selection in nature, leads to loss of diversity within the evolving population and premature convergence of the algorithm, hindering the discovery of many different solutions.

Here we present an alternative abstraction of artificial evolution, which does not require the formulation of a composite fitness function. Inspired from viability theory in dynamical systems, natural evolution and ethology, the proposed method puts emphasis on the elimination of individuals that do not meet a set of changing criteria, which are defined on the problem objectives and constraints.

Experimental results show that the proposed method maintains higher diversity in the evolving population and generates more unique solutions when compared to classical competition-based evolutionary algorithms. Our findings suggest that incorporating viability principles into evolutionary algorithms can significantly improve the applicability and effectiveness of evolutionary methods to numerous complex problems of science and engineering, ranging from protein structure prediction to aircraft wing design.

Selection methods regulate evolution of cooperation in digital evolution

A key, yet often neglected, component of digital evolution and evolutionary models is the 'selection method' which assigns fitness (number of offspring) to individuals based on their performance scores (efficiency in performing tasks). Here, we study with formal analysis and numerical experiments the evolution of cooperation under the five most common selection methods (proportionate, rank, truncation-proportionate, truncation-uniform and tournament). We consider related individuals engaging in a Prisoner's Dilemma game where individuals can either cooperate or defect. A cooperator pays a cost, whereas its partner receives a benefit, which affect their performance scores. These performance scores are translated into fitness by one of the five selection methods. We show that cooperation is positively associated with the relatedness between individuals under all selection methods. By contrast, the change in the performance benefit of cooperation affects the populations' average level of cooperation only under the proportionate methods. We also demonstrate that the truncation and tournament methods may introduce negative frequency-dependence and lead to the evolution of polymorphic populations. Using the example of the evolution of cooperation, we show that the choice of selection method, though it is often marginalized, can considerably affect the evolutionary dynamics.

In silico experimental evolution: a tool to test evolutionary scenarios

Comparative genomics has revealed that some species have exceptional genomes, compared to their closest relatives. For instance, some species have undergone a strong reduction of their genome with a drastic reduction of their genic repertoire. Deciphering the causes of these atypical trajectories can be very difficult because of the many phenomena that are intertwined during their evolution (e.g. changes of population size, environment structure and dynamics, selection strength, mutation rates...). Here we propose a methodology based on synthetic experiments to test the individual effect of these phenomena on a population of simulated organisms. We developed an evolutionary model--aevol--in which evolutionary conditions can be changed one at a time to test their effects on genome size and organization (e.g. coding ratio). To illustrate the proposed approach, we used aevol to test the effects of a strong reduction in the selection strength on a population of (simulated) bacteria. Our results show that this reduction of selection strength leads to a genome reduction of ~35% with a slight loss of coding sequences (~15% of the genes are lost--mainly those for which the contribution to fitness is the lowest). More surprisingly, under a low selection strength, genomes undergo a strong reduction of the noncoding compartment (~55% of the noncoding sequences being lost). These results are consistent with what is observed in reduced Prochlorococcus strains (marine cyanobacteria) when compared to close relatives.

Microbial evolution in vivo and in silico: methods and applications

Microbial evolution has been extensively studied in the past fifty years, which has lead to seminal discoveries that have shaped our understanding of evolutionary forces and dynamics. It is only recently however, that transformative technologies and computational advances have enabled a larger in-scale and in-depth investigation of the genetic basis and mechanistic underpinnings of evolutionary adaptation. In this review we focus on the strengths and limitations of in vivo and in silico techniques for studying microbial evolution in the laboratory, and we discuss how these complementary approaches can be integrated in a unifying framework for elucidating microbial evolution.

Boosting forward-time population genetic simulators through genotype compression

BACKGROUND

Forward-time population genetic simulations play a central role in deriving and testing evolutionary hypotheses. Such simulations may be data-intensive, depending on the settings to the various parameters controlling them. In particular, for certain settings, the data footprint may quickly exceed the memory of a single compute node.

RESULTS

We develop a novel and general method for addressing the memory issue inherent in forward-time simulations by compressing and decompressing, in real-time, active and ancestral genotypes, while carefully accounting for the time overhead. We propose a general graph data structure for compressing the genotype space explored during a simulation run, along with efficient algorithms for constructing and updating compressed genotypes which support both mutation and recombination. We tested the performance of our method in very large-scale simulations. Results show that our method not only scales well, but that it also overcomes memory issues that would cripple existing tools.

CONCLUSIONS

As evolutionary analyses are being increasingly performed on genomes, pathways, and networks, particularly in the era of systems biology, scaling population genetic simulators to handle large-scale simulations is crucial. We believe our method offers a significant step in that direction. Further, the techniques we provide are generic and can be integrated with existing population genetic simulators to boost their performance in terms of memory usage.

Evolution of music by public choice

Music evolves as composers, performers, and consumers favor some musical variants over others. To investigate the role of consumer selection, we constructed a Darwinian music engine consisting of a population of short audio loops that sexually reproduce and mutate. This population evolved for 2,513 generations under the selective influence of 6,931 consumers who rated the loops' aesthetic qualities. We found that the loops quickly evolved into music attributable, in part, to the evolution of aesthetically pleasing chords and rhythms. Later, however, evolution slowed. Applying the Price equation, a general description of evolutionary processes, we found that this stasis was mostly attributable to a decrease in the fidelity of transmission. Our experiment shows how cultural dynamics can be explained in terms of competing evolutionary forces.

The use of information theory in evolutionary biology

Information is a key concept in evolutionary biology. Information stored in a biological organism's genome is used to generate the organism and to maintain and control it. Information is also that which evolves. When a population adapts to a local environment, information about this environment is fixed in a representative genome. However, when an environment changes, information can be lost. At the same time, information is processed by animal brains to survive in complex environments, and the capacity for information processing also evolves. Here, I review applications of information theory to the evolution of proteins and to the evolution of information processing in simulated agents that adapt to perform a complex task.

Historical contingency affects signaling strategies and competitive abilities in evolving populations of simulated robots

One of the key innovations during the evolution of life on earth has been the emergence of efficient communication systems, yet little is known about the causes and consequences of the great diversity within and between species. By conducting experimental evolution in 20 independently evolving populations of cooperatively foraging simulated robots, we found that historical contingency in the occurrence order of novel phenotypic traits resulted in the emergence of two distinct communication strategies. The more complex foraging strategy was less efficient than the simpler strategy. However, when the 20 populations were placed in competition with each other, the populations with the more complex strategy outperformed the populations with the less complex strategy. These results demonstrate a tradeoff between communication efficiency and robustness and suggest that stochastic events have important effects on signal evolution and the outcome of competition between distinct populations.

Robust monomer-distribution biosignatures in evolving digital biota

Because organisms synthesize component molecules at rates that reflect those molecules' adaptive utility, we expect a population of biota to leave a distinctive chemical signature on its environment that is anomalous given the local (abiotic) chemistry. We observe the same effect in the distribution of computer instructions used by an evolving population of digital organisms, and we characterize the robustness of the evolved signature with respect to a number of different changes in the system's physics. The observed instruction abundance anomaly has features that are consistent over a large number of evolutionary trials and alterations in system parameters, which makes it a candidate for a non-Earth-centric life diagnostic.

Systems-biology approaches for predicting genomic evolution

Is evolution predictable at the molecular level? The ambitious goal to answer this question requires an understanding of the mutational effects that govern the complex relationship between genotype and phenotype. In practice, it involves integrating systems-biology modelling, microbial laboratory evolution experiments and large-scale mutational analyses - a feat that is made possible by the recent availability of the necessary computational tools and experimental techniques. This Review investigates recent progresses in mapping evolutionary trajectories and discusses the degree to which these predictions are realistic.

Information content of colored motifs in complex networks

We study complex networks in which the nodes are tagged with different colors depending on their function (colored graphs), using information theory applied to the distribution of motifs in such networks. We find that colored motifs can be viewed as the building blocks of the networks (much more than the uncolored structural motifs can be) and that the relative frequency with which these motifs appear in the network can be used to define its information content. This information is defined in such a way that a network with random coloration (but keeping the relative number of nodes with different colors the same) has zero color information content. Thus, colored motif information captures the exceptionality of coloring in the motifs that is maintained via selection. We study the motif information content of the C. elegans brain as well as the evolution of colored motif information in networks that reflect the interaction between instructions in genomes of digital life organisms. While we find that colored motif information appears to capture essential functionality in the C. elegans brain (where the color assignment of nodes is straightforward), it is not obvious whether the colored motif information content always increases during evolution, as would be expected from a measure that captures network complexity. For a single choice of color assignment of instructions in the digital life form Avida, we find rather that colored motif information content increases or decreases during evolution, depending on how the genomes are organized, and therefore could be an interesting tool to dissect genomic rearrangements.

The effects of low-impact mutations in digital organisms

Background

Avida is a computer program that performs evolution experiments with digital organisms. Previous work has used the program to study the evolutionary origin of complex features, namely logic operations, but has consistently used extremely large mutational fitness effects. The present study uses Avida to better understand the role of low-impact mutations in evolution.

Results

When mutational fitness effects were approximately 0.075 or less, no new logic operations evolved, and those that had previously evolved were lost. When fitness effects were approximately 0.2, only half of the operations evolved, reflecting a threshold for selection breakdown. In contrast, when Avida's default fitness effects were used, all operations routinely evolved to high frequencies and fitness increased by an average of 20 million in only 10,000 generations.

Conclusions

Avidian organisms evolve new logic operations only when mutations producing them are assigned high-impact fitness effects. Furthermore, purifying selection cannot protect operations with low-impact benefits from mutational deterioration. These results suggest that selection breaks down for low-impact mutations below a certain fitness effect, the selection threshold. Experiments using biologically relevant parameter settings show the tendency for increasing genetic load to lead to loss of biological functionality. An understanding of such genetic deterioration is relevant to human disease, and may be applicable to the control of pathogens by use of lethal mutagenesis.

Monomer abundance distribution patterns as a universal biosignature: examples from terrestrial and digital life

Organisms leave a distinctive chemical signature in their environment because they synthesize those molecules that maximize their fitness. As a result, the relative concentrations of related chemical monomers in life-bearing environmental samples reflect, in part, those compounds' adaptive utility. In contrast, rates of molecular synthesis in a lifeless environment are dictated by reaction kinetics and thermodynamics, so concentrations of related monomers in abiotic samples tend to exhibit specific patterns dominated by small, easily formed, low-formation-energy molecules. We contend that this distinction can serve as a universal biosignature: the measurement of chemical concentration ratios that belie formation kinetics or equilibrium thermodynamics indicates the likely presence of life. We explore the features of this biosignature as observed in amino acids and carboxylic acids, using published data from numerous studies of terrestrial sediments, abiotic (spark, UV, and high-energy proton) synthesis experiments, and meteorite bodies. We then compare these data to the results of experimental studies of an evolving digital life system. We observe the robust and repeatable evolution of an analogous biosignature in a digital lifeform, suggesting that evolutionary selection necessarily constrains organism composition and that the monomer abundance biosignature phenomenon is universal to evolved biosystems.

Evolutionary dynamics, epistatic interactions, and biological information

We investigate a definition of biological information that connects population genetics with the tools of information theory by focusing on the distribution of genotypes found in a population. Previous research has treated loci as non-interacting by making specific approximations in the calculation of information-theoretic quantities. We expand earlier mathematical forms to include epistasis, or interactions between mutations at all pairs of loci. Application of our improved measure of biological information to evolution on two-locus, two-allele fitness landscapes demonstrates that mutual information between loci reflects epistatic interaction of mutations. Finally, we consider four-locus, two-allele fitness landscapes with modular structure. As modular interactions are inherently epistatic, we demonstrate that our refined approximation provides insight into the underlying structure of these non-trivial fitness landscapes.

Introduction to focus issue: genetic interactions

The perturbation of a gene in an organism's genome often causes changes in the organism's observable properties or phenotypes. It is not obvious a priori whether the simultaneous perturbation of two genes produces a phenotypic change that is easily predictable from the changes caused by individual perturbations. In fact, this is often not the case: the nonlinearity and interdependence between genetic variants in determining phenotypes, also known as epistasis, is a prevalent phenomenon in biological systems. This focus issue presents recent developments in the study of epistasis and genetic interactions, emphasizing the broad implications of this phenomenon in evolutionary biology, functional genomics, and human diseases.

Biological species is the only possible form of existence for higher organisms: the evolutionary meaning of sexual reproduction

Consistent holistic view of sexual species as the highest form of biological existence is presented. The Weismann's idea that sex and recombination provide the variation for the natural selection to act upon is dominated in most discussions of the biological meaning of the sexual reproduction. Here, the idea is substantiated that the main advantage of sex is the opposite: the ability to counteract not only extinction but further evolution as well. Living systems live long owing to their ability to reproduce themselves with a high fidelity. Simple organisms (like bacteria) reach the continued existence due to the high fidelity of individual genome replication. In organisms with a large genome and complex development, the achievable fidelity of DNA replication is not enough for the precise reproduction of the genome. Such species must be capable of surviving and must remain unchanged in spite of the continuous changes of their genes. This problem has no solution in the frame of asexual ("homeogenomic") lineages. They would rapidly degrade and become extinct or blurred out in the course of the reckless evolution. The core outcome of the transition to sexual reproduction was the creation of multiorganismic entity - biological species. Individual organisms forfeited their ability to reproduce autonomously. It implies that individual organisms forfeited their ability to substantive evolution. They evolve as a part of the biological species. In case of obligatory sexuality, there is no such a thing as synchronic multi-level selection. Natural selection cannot select anything that is not a unit of reproduction. Hierarchy in biology implies the functional predestination of the parts for the sake of the whole. A crucial feature of the sexual reproduction is the formation of genomes of individual organisms by random picking them over from the continuously shuffled gene pool instead of the direct replication of the ancestor's genome. A clear anti-evolutionary consequence of the sexuality is evident from the fact that the genotypes of the individuals with an enhanced competitiveness are not transmitted to the next generation. Instead, after mating with "ordinary" individuals, these genotypes scatter and rearrange in new gene combinations, thus preventing the winner from exploiting the success.

Quasispecies theory for finite populations

We present stochastic, finite-population formulations of the Crow-Kimura and Eigen models of quasispecies theory, for fitness functions that depend in an arbitrary way on the number of mutations from the wild type. We include back mutations in our description. We show that the fluctuation of the population numbers about the average values is exceedingly large in these physical models of evolution. We further show that horizontal gene transfer reduces by orders of magnitude the fluctuations in the population numbers and reduces the accumulation of deleterious mutations in the finite population due to Muller's ratchet. Indeed, the population sizes needed to converge to the infinite population limit are often larger than those found in nature for smooth fitness functions in the absence of horizontal gene transfer. These analytical results are derived for the steady state by means of a field-theoretic representation. Numerical results are presented that indicate horizontal gene transfer speeds up the dynamics of evolution as well.

The dynamics of adaptation on correlated fitness landscapes

Evolutionary theory predicts that a population in a new environment will accumulate adaptive substitutions, but precisely how they accumulate is poorly understood. The dynamics of adaptation depend on the underlying fitness landscape. Virtually nothing is known about fitness landscapes in nature, and few methods allow us to infer the landscape from empirical data. With a view toward this inference problem, we have developed a theory that, in the weak-mutation limit, predicts how a population's mean fitness and the number of accumulated substitutions are expected to increase over time, depending on the underlying fitness landscape. We find that fitness and substitution trajectories depend not on the full distribution of fitness effects of available mutations but rather on the expected fixation probability and the expected fitness increment of mutations. We introduce a scheme that classifies landscapes in terms of the qualitative evolutionary dynamics they produce. We show that linear substitution trajectories, long considered the hallmark of neutral evolution, can arise even when mutations are strongly selected. Our results provide a basis for understanding the dynamics of adaptation and for inferring properties of an organism's fitness landscape from temporal data. Applying these methods to data from a long-term experiment, we infer the sign and strength of epistasis among beneficial mutations in the Escherichia coli genome.

Gene divergence and pathway duplication in the metabolic network of yeast and digital organisms

We have studied the metabolic gene-function network in yeast and digital organisms evolved in the artificial life platform Avida. The gene-function network is a bipartite network in which a link exists between a gene and a function (pathway) if that function depends on that gene, and can also be viewed as a decomposition of the more traditional functional gene networks, where two genes are linked if they share any function. We show that the gene-function network exhibits two distinct degree distributions: the gene degree distribution is scale-free while the pathway distribution is exponential. This is true for both yeast and digital organisms, which suggests that this is a general property of evolving systems, and we propose that the scale-free gene degree distribution is due to pathway duplication, i.e. the development of a new pathway where the original function is still retained. Pathway duplication would serve as preferential attachment for the genes, and the experiments with Avida revealed precisely this; genes involved in many pathways are more likely to increase their connectivity. Measuring the overlap between different pathways, in terms of the genes that constitute them, showed that pathway duplication also is a likely mechanism in yeast evolution. This analysis sheds new light on the evolution of genes and functionality, and suggests that function duplication could be an important mechanism in evolution.