The effects of low-impact mutations in digital organisms
© Nelson and Sanford; licensee BioMed Central Ltd. 2011
Received: 27 January 2011
Accepted: 18 April 2011
Published: 18 April 2011
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.
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.
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.
The standard explanation for the origin of biological complexity is that it arises through the Darwinian process of mutation and natural selection. Beneficial mutations accumulate through positive selection, and deleterious mutations tend to be eliminated by purifying selection. However, developments in genomics suggest theoretical problems with this view, and many features of living systems cannot be explained without recourse to nonadaptive processes [1–4].
Because of the slow pace of evolutionary change, it has generally been difficult to empirically test long-term evolutionary scenarios. A computational approach known as digital genetics [5, 6] attempts to overcome this limitation by using digital organisms, short computer programs that replicate and compete in a virtual environment. Generations take only a few seconds, making it possible to observe the outcome of large numbers of mutation and replication events in relatively short periods of real time. Further, the user is able to alter parameters of interest (e.g., mutation rates) to observe their influence on important population factors (e.g., fitness).
Early versions of digital life culminated in the program Tierra , which demonstrated adaptive genome shrinkage, cooperation, and parasitism. Genomes were simulated as computer code, distinguishing the software from numerical simulation. Mutating digital organisms competed for computer processing time, undergoing adaptive change over many generations. Recognizing the importance of local interactions, the program Avida [7, 8] advanced the field by implementing a virtual world in which organisms were housed on a two-dimensional grid and underwent interactions with neighbors.
Researchers have claimed a high degree of biological relevance for Avida, comparing its digital organisms to organic viruses . Titles like "The biology of digital organisms" , "Evolution of biological complexity" , and "Testing Darwin"  evidence Avida's impact on biological theory. In addition to the evolution of biological complexity [11, 13], the software has been used to study the evolution of sex [9, 14, 15], the evolution of altruism , the dynamics of long-term adaptation [17–21], ecosystem dynamics [19, 22–24], and the effects of mutation on genetic architecture [14, 25–28], among other topics.
Avida is used in the present study to better understand the evolutionary consequences of low-impact mutations in digital organisms. Though many studies report the occurrence of neutral mutations, Eyre-Walker & Keightley  note that:
... it seems unlikely that any mutation is truly neutral in the sense that it has no effect on fitness. All mutations must have some effect, even if that effect is vanishingly small. However, there is a class of mutations that we can term effectively neutral... As such, the definition of neutrality is operational rather than functional; it depends on whether natural selection is effective on the mutation in the population or the genomic context in which it segregates, not solely on the effect of the mutation on fitness.
This point applies to viruses as well as more complex systems . The term selection threshold has been introduced [Gibson P, et al., in preparation] to describe the mutational fitness effect that marks the "tipping point" between natural selection and random genetic drift in an evolving system. Mutations with fitness effects below this critical value are primarily affected by random genetic drift. One of the first to allude to this phenomenon was Muller , who noted: "There comes a level of advantage... that is too small to be effectively seized upon by selection."
The selection threshold is elevated by any factor that influences replication rate in a manner independent of the genotype, decreasing the efficacy of selection as more mutations behave in a neutral fashion. Population size has typically been the primary focus of these factors , and its role is described in Kimura's  well-known expression, |s| < 1/(2N e ). This inequality states that random genetic drift will dominate a mutation's fate if its selection coefficient (s) is less than the reciprocal of twice the effective population size (N e ). However, numerous other factors also influence the selection threshold, including environmental noise and developmental canalization, and the efficacy of selection is highly dependent on the complexity of the system under study.
The present study takes an empirical approach to determining the selection threshold by measuring the mutational fitness effect at which selection successfully captures half of the beneficial mutations that arise. Previous experiments using Avida have studied the evolutionary emergence of complex features resulting from high-impact beneficial mutations . Avida's default settings provide mutational fitness effects of 1.0 - 31.0 for beneficial mutations that give rise to certain computational operations, where fitness effects are measured as w - 1, and w is the relative fitness of the organism expressing a given operation. For example, a mutation producing the NAND operation will multiply an organism's fitness by 2, corresponding to a fitness effect of 1.0. However, fitness effects this large are extremely rare in nature (see Discussion). In the present study, we approximate the selection threshold in Avida by performing experiments with more biologically common mutational fitness effects of 1.0 and below. The effects of low-impact mutations are explored and the biological relevance of digital life is discussed.
An experiment with Avida begins by seeding a two-dimensional grid with a short computer program (the ancestral organism) that has been designed to self-replicate. By default, a 60 × 60 grid is seeded with a single Avidian organism that consists of 100 computational instructions. This artificial geography allows the population to grow to a maximum of 3,600 organisms. Avidians replicate asexually for approximately 10,000 generations, incurring an average of 0.85 mutations per genome per generation. Mutations randomly substitute, insert, or delete single instructions in an Avidian genome, drawing upon 26 available instructions defined in the software. The ancestral genome devotes about 15 instructions to the essential replication code, while the remaining 85 positions are occupied by benign no-operation instructions, analogous to inert "junk DNA" that can be used as raw material for evolutionary tinkering.
Once an experiment begins, replication ensues, and multiple organisms arise and compete with one another. When an Avidian replicates, its offspring is randomly placed in one of eight positions surrounding the parent organism, effectively killing the previous resident. Speed of replication therefore defines fitness in Avida; the programs that replicate fastest replace their slower counterparts and increase in number.
Speed of replication is itself determined by two factors. The first and primary way that Avidians replicate faster is by earning additional computer resources. The allocation of computer time is based upon an organism's merit, a numerical value that reflects its ability to perform one or more simple computational tasks. Specifically, Avidians may evolve any of nine logic operations, for which they are rewarded with additional computer time to execute and replicate their genomes. Secondarily, speed of replication in Avida is influenced by genome size. Organisms with larger genomes naturally require more computer time and replicate at a slightly slower rate. However, under default settings, this factor is offset by artificially rewarding larger genomes with additional computer time, such that genome size is not under direct selection in most experiments. More detailed descriptions of the software are available elsewhere [33–35].
Default rewards for performing nine logic operations in Avida
Number of NAND operations needed (n)
Default multiplicative bonus (2 n )
Default fitness effect (w- 1)
~(A and B)
A and B
(A or ~B); (~A or B)
A or B
(A and ~B); (~A and B)
~A and ~B
(A and ~B) or (~A and B)
(A and B) or (~A and ~B)
Avida may be conceptualized as a computational Darwinian search designed to discover the EQU operation. The simplest operations in Avida are easy to evolve, i.e., NAND and NOT are performed by a single genomic instruction, provided instructions for correctly inputting and outputting numbers are present. Any logic operation can be built using different combinations of NAND and NOT. Therefore, EQU can itself be constructed using any of the eight simpler operations as precursors, providing a scalable fitness landscape for the evolution of complexity - beneficial changes are useful for constructing more complex beneficial features. When NAND or NOT arises, the software rewards the lucky organism by doubling its fitness. Fitness bonuses for the other operations increase exponentially with complexity (Table 1). The evolution of EQU may therefore proceed one advantageous step at a time, each step requiring relatively few mutations. Dembski and Marks  have suggested the term "stair step active information" to describe this type of reward scheme.
Some of the ways Avida has been implemented (e.g., its parameter settings) are distinctly "un-biological" . These factors include the distribution of mutational fitness effects, the fitness terrain, and the artificial rewards given to organisms with larger genomes. The present study pursues several lines of experimentation with altered mutational fitness effects to improve biological relevance and aid in the interpretation of Avida results. The first set of experiments removed merit bonuses to determine which logic operations arise by mutation alone, without selection. The second set of experiments examined Avida's default settings to quantify typical aspects of evolutionary change in this system. In order to test the hypothesis that mutation pressure prevents the fixation of beneficial operations in Avida, a third set of experiments examined logic operation frequencies at a reduced mutation rate. Finally, a fourth set of experiments implemented fitness effects falling in the normal biological range (0.01 - 1.0), rather than Avida's default range (1.0 - 31.0). The effects on evolutionary dynamics were observed.
Mutation and drift
Twenty experiments were performed in which no logic operations were rewarded. Across these experiments, an average of 6.4 (± 0.8) operations drifted into a population at least once over the course of 10,000 generations, indicating that they are easily produced by random mutation. Because of this, a distinction was made between those operations that arose by chance in Avida (those that arose) and those that selection was able to propagate (those that successfully evolved, i.e., rose to a frequency of 50% or greater, following the precedent of biological studies [37, 38]).
Dynamics of mutation and drift for nine logic operations in Avida
Proportion of experiments in which operation arose by mutation
Average maximum frequency in population
Average maximum number of organisms
Maximum frequency observed
Maximum number of organisms observed
0.027 (± 0.0062)
0.017 (± 0.0046)
0.0015 (± 0.00099)
0.0063 (± 0.0020)
0.00089 (± 0.00091)
0.0030 (± 0.0019)
0.00053 (± 0.00062)
0 (± 0)
0 (± 0)
Evolution under default settings
Mutation pressure and clonal interference
Interestingly, no operations reached fixation under default settings, despite their remarkably high fitness bonuses. The average end-of-experiment frequency for operations that successfully evolved was only 84.5% (± 13.5%). This contrasts with the rapid fixation of high-impact beneficial mutations observed in biological experiments. For example, in one study of E. coli, the Rbs- mutation increased fitness only by about 1.4%, yet reached fixation (97-100%) in only 2,000 generations.
The effects of mutation rate on phenotype frequencies
Frequency with default mutation rate
Frequency with reduced mutation rate
Evolutionary consequences of low-impact mutational fitness effects
To explore the evolutionary consequences of low-impact mutational fitness effects in Avida, experiments were performed with multiplicative fitness effects of 0, 0.01, 0.05, 0.075, 0.1, 0.25, 0.5, and 1.0, with 0 being neutral and 1.0 corresponding to a doubling of fitness (100% increase). This allowed an empirical estimation of Avida's selection threshold, the critical "tipping point" between random genetic drift and natural selection. Because most operations arise readily by chance in Avida, evolution of an individual operation was again considered successful only if its end-of-experiment frequency was 50% or greater. Two sets of 20 replicates were performed, one for beneficial mutations and one for deleterious mutations, with each replicate consisting of eight experiments (one experiment for each fitness effect). For beneficial mutations, experiments were simply initiated with uniform fitness effects of the specified value (e.g., for a fitness effect of 0.1, all nine operations multiplied fitness by 1.1). For deleterious mutations, experiments were performed first under Avida's default settings to allow the evolution of complexity, and then continued for an additional 10,000 generations with the alternative beneficial fitness effects. A range of fitness effects could also have been used, with rare operations incurring greater benefits; however, uniform fitness effects were ideal for the purpose of approximating the selection threshold in Avida, and using a range would not appreciably alter our results. Since mutation pressure is a significant force in Avida, it was expected that the existing operations would incur deactivating mutations, and that the fitness bonuses would determine selection's efficacy in maintaining those operations.
Although Avida has routinely been used to address biological questions, some aspects of the program are not amenable to direct biological comparison. For example, key terms such as nucleotide, gene, heritability, selection, and fertility lack a clear equivalent in the software. Because of this, several approximations were necessary in this study. Allele frequencies were measured as phenotype frequencies, ignoring the potential for chance performance. Mutation rates were measured as the rate of random substitution of single instructions, though these monomers can perform multiple computations and are not comparable to biological nucleotides. Generation times changed substantially over the course of a typical experiment, so the average end-of-experiment generation time was used to measure experiment length. Finally, genome size also fluctuated in these experiments, causing the genomic mutation rate to change. For simplicity, the mutation rates reported were those for the ancestral genome size (100).
In these experiments, all but two logic operations in Avida arose via mutation alone, despite conferring no fitness rewards (Table 2). Most operations are therefore very simple to produce in the Avida environment, with relatively short waiting times. The genomic monomers (instructions) themselves do most of the computational work that these operations require; this underlying information is included in the artificial physics of Avida and is not subject to mutational change. Interestingly, un-rewarded operations did not accumulate to produce the more complex operations XOR and EQU. This suggests difficulties for traditional models of evolution by gene duplication in which novel functions arise by neofunctionalization of unconstrained loci [41, 42]. Previous work has explored the evolution of EQU when other operations are made neutral [13, 39], and further Avida studies should explore the dynamics of neutral evolution in digital organisms.
Several studies have focused on the evolution of "robustness" in Avida under elevated mutation rates [25–28, 43]. These studies have shown that, when functional genomes experience high mutation rates, functionality is generally lost, with some operations evolving to utilize fewer genomic positions. This is consistent with the results reported here, which suggest that mutation pressure is a significant force preventing the fixation of beneficial genotypes in Avida (Table 3). Reduced mutation rates allowed advantageous phenotypes to reach higher frequencies; however, fewer operations evolved, evidencing a tradeoff between reducing genetic load and increasing the waiting time to beneficial mutation.
The decelerating rate of adaptive change in Avida (Figure 1) is somewhat reminiscent of biological evolution experiments, e.g., with bacteriophage  and E. coli[45, 46]. However, the explosive fitness increases observed in Avida are roughly seven orders of magnitude greater than those observed in biological experiments of similar duration. Because fitness is defined as relative replication rate in Avida, the program's results may be directly compared with those from biological studies. For example, in experiments with E. coli, growth rate increased by an average of ~37% after 2,000 generations , ~48% after 10,000 generations , and ~75% after 20,000 generations . These changes, resulting from numerous mutations, are negligible compared to those observed under Avida's default settings. Yet the fitness leaps observed in Avida are due primarily to the large multiplicative fitness effects of just nine simple innovations. For example, when fitness effects for all logic operations were set to 1.0, the average end-of-experiment fitness plummeted from almost 20,000,000 to just 180 (still an immense increase relative to biological organisms).
An analogy will help to elucidate the preceding point. Consider species A, a large mammal with a generation time of 30 years, and species B, a bacterial species with a generation time of 1 day. In terms of replication rate, species B is about 10,950 times fitter than species A. Yet this number pales in comparison to the increases observed in Avida. After only 10,000 generations, the fitness (replication rate) of digital organisms in Avida increased by 20 million. Such an increase would allow mammalian species A to evolve a generation time of just 1.6 minutes in this time. This phenomenon occurs because the bonuses readily available to digital organisms in Avida are large and multiplicative, producing exorbitant gains in fitness (i.e., the product of all possible bonuses is 22 × 42 × 82 × 162 × 32 = 33,554,432). Fitness bonuses this large are extremely rare in nature (but see references [48, 49]).
Mutations of smaller effect (i.e., fitness effects of ≤ 1.0) can occur in Avida when the generation time is altered by insertions or deletions within an organism's replication loop. However, the rewards gained by performing logic operations dominate fitness dynamics in Avida, and these are the only fitness effects that can be user-specified. Mutations disabling any of the evolved operations have similarly large (but not identical) deleterious effects. It is our view that the distribution of fitness effects used in Avida has severely limited its relevance to biological systems.
Though many details of the biological distribution of mutational fitness effects have yet to be understood , a general picture has emerged. There is a continuum of fitness effects and, with few exceptions [51, 52], advantageous mutations are exponentially distributed, being much more rare than deleterious mutations [29, 30, 53–56]. The distribution of deleterious mutations is likely multimodal, with a distinct class being lethal and another class having very small effects . In most systems studied, deleterious mutations of small effect are more abundant than those of large effect [29, 54], such that selection coefficients in the range of 0.01 to 0.1 are considered large . For example, over 90% of gene knockouts in E. coli are viable , decreasing fitness by an average of ≤ 3% . Similarly, a recent study of mutations in Salmonella typhimurium reported average deleterious selection coefficients of 0.0096 and 0.0131 for synonymous and nonsynonymous mutations, respectively. No significantly advantageous mutations were found, and no mutations caused a complete loss of function.
Viruses are somewhat peculiar because of their high mutational sensitivity, with approximately 20 to 41% of mutations being lethal, and many mutations being neutral or nearly neutral [30, 55, 56, 59]. However, viable deleterious mutations of small effect are still more common than those of large effect. A recent review  of several viruses reported a mean fitness effect of 0.10 to 0.13 (though some estimates have been considerably lower ). Lind et al. note that the high frequencies of neutral mutations reported in some studies may be a consequence of assays that lack sufficient sensitivity, and that nearly neutral mutations may be more common than previously thought. Some studies have reported the fixation of highly beneficial mutations in viruses [48, 49], but they have not consistently measured fitness effects as they are defined here. These reports have suggested that beneficial mutations in viruses may be described by a uniform distribution with an upper bound . It is clear that mutational fitness effects in biological organisms are substantially smaller than those used heretofore in digital life research.
The term selection threshold has been introduced to describe the critical mutational fitness effect for which natural selection and random genetic drift contribute equally to a mutation's fate [Gibson P, et al., in preparation]. Though many studies report the occurrence of neutral mutations, it is unlikely that any mutation is truly neutral , including those in viruses . It is increasingly being recognized that most mutations in multicellular organisms fall far below the selection threshold, having fitness effects so slight that they cannot be measured [29, 49]. It is also noteworthy that highly beneficial mutations have gone undetected in most evolution experiments with eukaryotes . Clearly, the results of evolution experiments with microorganisms cannot be extrapolated to eukaryotes with larger genomes, greater phenotypic complexity, and smaller population sizes. In light of this, some may ask whether the results of experiments with digital organisms have any relevance to living systems. We conclude that digital genetics is a valid platform for studying some biological questions, but that the applicability of results will depend critically upon the parameters used.
Population size has routinely been used as the sole predictor of selection efficacy. To our knowledge, the present study is the first that uses an empirical approach to estimate the selection threshold in an evolving system. This approach implicitly considers all factors affecting selection, including (but not limited to) population size, the probabilistic nature of selection, and environmental effects. We find that, given the sources of noise inherent in the Avida world, mutations with fitness effects below the 0.075 - 0.1 range are entirely invisible to selection, despite arising frequently. Fitness effects of approximately 0.2 are necessary for selection to successfully capture half of the beneficial mutations that arise, corresponding to the selection threshold. Though the value of this threshold is certain to differ among biological and digital systems, its existence has important theoretical and medical implications. Other Avida experiments  using single organisms and truncation selection have improved the program's performance, suggesting that local interactions and the probabilistic nature of selection are important sources of noise in the Avida world. Some readers might object that, due to this noise, fitness effects in Avida do not directly correspond to fitness effects in biological organisms. If there is truly no correspondence, then experiments using Avida are not capable of shedding light on biological questions. However, there are several reasons why these experiments are broadly relevant to biology. Importantly, there is also noise in biology. The information contained in the heritable material is processed through multiple levels, including transcription, mRNA processing, protein folding, physiological interactions, and more. Each level is subject to mechanisms of canalization and homeostasis that obscure the effects of mutations on fitness. Further, most noise in Avida may be attributed to the probabilistic nature of selection, yet probability selection is also operative in nature, and may be considerably weaker  than the scheme implemented in Avida. It is therefore probable that biological organisms experience more fitness noise than digital organisms.
The present study used uniform mutational fitness effects. A range of fitness effects could also have been used, with rare logic operations incurring greater benefits. However, uniform fitness effects are often employed, and were ideal for the purpose of approximating the selection threshold. Even the simplest operations did not evolve at fitness effects of ≤ 0.075, demonstrating that the threshold exists independent of a mutation's rarity or the length of an evolutionary experiment. Moreover, because each operation is rewarded only once per organism, the evolution of simpler operations should not prevent the subsequent evolution of additional complexity. Further research should attempt to better approximate the selection threshold in this and other systems using varying fitness effect distributions.
The distribution of mutational fitness effects has serious implications for genetic disease. Numerous analyses have confirmed that the accumulation of slightly deleterious mutations can cause gradual fitness loss leading to extinction in asexual species [31, 62–66], and similar processes are relevant to sexual species [67, 68], including humans [69–74]. The results reported here reveal a quantifiable selection threshold, below which random genetic drift dominates the behavior of advantageous and deleterious mutations alike. Biological studies elucidating the extent of sequence-dependent functionality in nonprotein-coding DNA regions will continue to inform estimates of the rates of mutations affecting fitness in various species, allowing a realistic evaluation of the severity of genetic decline that can be expected in coming generations.
We observed that, when fitness effects in Avida are small, all advantageous logic operations are lost. Though digital organisms are peculiar in that they can survive such a loss, these data confirm that the accumulation of slightly deleterious mutations can lead to decreasing biological functionality and potentially eventual extinction. Because deleterious mutations are much more common than advantageous mutations in most systems studied, reduction in the efficacy of selection imposes strong directionality on evolution by favoring the fixation of deleterious mutations . The conditions under which fitness recovery may be possible  should be studied more thoroughly using computational approaches. An understanding of these issues may be applicable to the lethal mutagenesis of pathogen populations [66, 76], and is relevant to human health . It is clear that mutation accumulation may affect human health at various levels, including the nuclear and mitochondrial genomes as well as immune cells. For example, mutation accumulation may play a key role in the deterioration of the immune system during HIV infection . It may also influence the longevity of pathogen populations. Various factors relevant to pathogen attenuation, including the consequences of periodic bottlenecking and elevated mutation rates, should be studied using computational models. Understanding the interplay of these co-evolutionary processes may allow for substantial advances in medicine, including novel treatments and an increased awareness of the role of mutation in disease.
Avida has previously been used as a powerful demonstration of adaptation resulting from high-impact beneficial mutations. However, there are several ways in which Avida's default settings produce results which conflict with observations from biological experiments. Precursors necessary for the most complex logic operation in the program, EQU, are frequently produced by random mutation, yet confer very large fitness rewards. Fitness effects of beneficial mutations under Avida's default settings range from 1.0 to 31.0, values that are extremely rare in the natural world. As a result, fitness increases by an average of 20 million in only 10,000 generations. This is roughly seven orders of magnitude greater than the changes observed in biological evolution experiments.
In contrast to Avida's default settings, most mutations in biological organisms are low-impact , and this class of mutations may dominate evolutionary change [1, 2]. When Avida is used with more realistic mutational fitness effects, it demonstrates a clear selection threshold. Mutations that influence fitness by approximately 20% or less come to be dominated by random genetic drift. Mutations that affect fitness by 7.5 - 10.0% or less are entirely invisible to selection in this system. These results provide evidence that low-impact mutations can present a substantial barrier to progressive evolution by natural selection. Understanding mutation is of primary importance, as selection depends on the mutational production of new genotypes. Numerous changes that would be beneficial may nevertheless fail to occur because mutation cannot produce them in the time available. Further, it is important for biologists to realistically appraise what selection can and cannot do under various circumstances. Selection may neither be necessary nor sufficient to explain numerous genomic or cellular features of complex organisms [2–4].
Future studies should explore the interaction of low-impact fitness effects with other evolutionary factors, such as alternative fitness terrains , to elucidate their synergistic effects on evolutionary dynamics. Additionally, researchers should attempt to further quantify the selection threshold for various systems, and determine the phenotypic consequences of accumulating low-impact mutations. The accumulation of slightly deleterious mutations may pose an important health risk for numerous species, including humans , and warrants further study using computational approaches. Reducing the rate of mutation in the human genome may be an important step in fighting genetic disease. Additionally, the connection between mutation accumulation and pathogen attenuation should be studied. Finally, we recommend that future experiments with digital organisms employ more biologically relevant mutational fitness effects.
This study used Avida version 2.8.1 . For all experiments, random number seeds were chosen randomly as an integer in the range 1 to 1000000000. These values are reported in the supplemental data. Settings were manipulated using the configuration files present in the Avida home folder.
In Avida, time is measured in arbitrary units called updates in which the population is allowed to execute about 30 instructions per organism. The ancestral organism executes 389 instructions per generation to execute and copy its genome. However, generation times change as organisms evolve. Average end-of-experiment generation time under default settings was 310, corresponding to approximately 9,679 generations over 100,000 updates. Thus 10,000 generations was used as an approximation of experiment length, with 1 generation corresponding to roughly 10 updates.
Fitness in Avida is measured as an organism's total merit divided by its generation time. Thus an increase in merit will increase fitness, while an increase in generation time will decrease it. However, this value has no intuitive meaning and the software does not consistently report it, e.g., at generation 0. For consistency and ease of biological comparison, fitness in this study was re-calculated and reported relative to the ancestral organism's fitness. Thus the average fitness of a generation was equivalent to the average merit divided by the product of average generation time and the ancestral organism's fitness (scaling the ancestral organism's fitness to 1.0). Fitness effects of mutations producing the logic operations were measured as w - 1, where w is the relative fitness of the organism carrying the mutation. Thus a mutation producing the EQU operation, multiplying fitness by 32, had a fitness effect of 31.0.
To study mutation and drift, 20 experiments were performed in which no logic operations were rewarded. Merit bonuses in the environment.cfg file were defined multiplicatively (type = mult) as 1.0 (value = 1.0), corresponding to fitness effects of 0. All other settings maintained their default values. The output file tasks.dat was examined to determine which operations arose in an experiment. As allele frequencies are not reported by Avida, phenotype frequencies were measured as the number of organisms performing a logic operation divided by the total number of organisms, the latter of which is reported in the count.dat file.
Logic operations readily arose independent of selection. In order to distinguish those mutations that selection promoted, a mutation was said to have simply arisen until its frequency reached 50%, at which point it was considered to have successfully evolved, a measure that has precedent [37, 38]. For analysis of end-of-experiment frequencies of logic operations, averages were taken using only nonzero values (operations that did not arise were not considered).
Thirty experiments were performed under default settings. The case study reported in Figure 1 occurred when using a random seed of 574423164. Thirty experiments were also performed with a reduced mutation rate. For these runs, the copy mutation rate was changed from the default of 0.0075 to 0.004 in the avida.cfg file (COPY_MUT_PROB 0.004). Probabilities of insertions and deletions at the time of replication each remained at 0.05. Thus the default mutation rate corresponded to an average rate of 0.85 mutations per genome per generation, and the lower rate to an average of 0.5 mutations per genome per generation. P-values reported in Table 3 were for two-tailed two-sample t-tests with equal variances (homoscedastic) comparing all nonzero end-of-experiment frequencies from the two environments. The case study shown in Figure 2 occurred when using default settings and a random seed of 13903545.
To study the consequences of alternative mutational fitness effects, merit bonuses were modified in the environment.cfg file. Fitness effects of 0, 0.01, 0.05, 0.075, 0.1, 0.25, 0.5, and 1.0 were used, with 0 being neutral and 1.0 corresponding to a doubling of fitness. Uniform effects were used, such that all advantageous operations conferred the same bonus. For example, for a fitness effect of 0.01, merit bonuses were defined multiplicatively (type = mult) as 1.01 (value = 1.01) for all nine operations. Two sets of 20 replicates were performed, one for beneficial mutations and one for deleterious mutations. Each replicate for beneficial mutations consisted of eight experiments, one for each of the fitness effects tested. Each replicate for deleterious mutations consisted of eight similar experiments in which evolution proceeded first under default settings, and then continued an additional 10,000 generations using the alternative fitness effects and a new random number seed. Because fitness bonuses were multiplicative in our experiments, fitness effects would in reality be slightly different for mutations creating and destroying the same operation. For example, a mutation creating an operation with a bonus of 1.25 would have a beneficial fitness effect of (1 - 1.25 / 1.0) = 0.25, but a mutation destroying the same operation would have a deleterious fitness effect of (1 - 1.0 / 1.25) = 0.2. For simplicity, the consequences of deleterious mutations were simply reported (Figure 3) relative to the beneficial fitness effect. Additionally, while deleterious fitness effects could have been indicated as negative values, the absolute values were used for simplicity.
This work was supported by Rainbow Technologies, Inc. The authors thank Winston Ewert for helpful comments regarding the Avida software.
- Kimura M: The Neutral Theory of Molecular Evolution. 1983, Cambridge: Cambridge University PressView ArticleGoogle Scholar
- Lynch M: The frailty of adaptive hypotheses for the origins of organismal complexity. Proc Natl Acad Sci USA. 2007, 104 (Suppl 1): 8597-8604. 10.1073/pnas.0702207104.PubMed CentralView ArticlePubMedGoogle Scholar
- Hughes AL: Looking for Darwin in all the wrong places: the misguided quest for positive selection at the nucleotide sequence level. Heredity. 2007, 99: 364-373. 10.1038/sj.hdy.6801031.View ArticlePubMedGoogle Scholar
- Koonin EV: Darwinian evolution in the light of genomics. Nucleic Acids Res. 2009, 37: 1011-1034. 10.1093/nar/gkp089.PubMed CentralView ArticlePubMedGoogle Scholar
- Ray TS: An approach to the synthesis of life. Artificial Life II, SFI Studies in the Sciences of Complexity, Vol. X: 1991. Edited by: Langton CG, Taylor C, Farmer JD, Rasmussen S. 1991, Addison-Wesley, 371-408.Google Scholar
- Adami C: Digital genetics: unravelling the genetic basis of evolution. Nat Rev Genet. 2006, 7: 109-118. 10.1038/nrg1771.View ArticlePubMedGoogle Scholar
- Adami C, Brown CT: Evolutionary learning in the 2D artificial life system "Avida.". Artificial Life IV: Proceedings of the Fourth International Workshop on the Synthesis and Simulation of Living Systems: 1994. Edited by: Brooks R, Maes P. 1994, Cambridge, MA: MIT Press, 377-381.Google Scholar
- Adami C: Introduction to Artificial Life. 1998, New York: Springer-VerlagView ArticleGoogle Scholar
- Misevic D, Lenski RE, Ofria C: Sexual reproduction and Muller's Ratchet in digital organisms. Artificial Life IX: Proceedings of the Ninth International Conference on the Simulation and Synthesis of Living Systems. Edited by: Pollack J, Bedau M, Husbands P, Ikegami T, Watson RA. 2004, Cambridge, MA: MIT Press, 340-345.Google Scholar
- Wilke CO, Adami C: The biology of digital organisms. Trends Ecol Evol. 2002, 17: 528-532. 10.1016/S0169-5347(02)02612-5.View ArticleGoogle Scholar
- Adami C, Ofria C, Collier TC: Evolution of biological complexity. Proc Natl Acad Sci USA. 2000, 97: 4463-4468. 10.1073/pnas.97.9.4463.PubMed CentralView ArticlePubMedGoogle Scholar
- Zimmer C: Testing Darwin. Discover. 2005, 26: 28-35.Google Scholar
- Lenski RE, Ofria C, Pennock RT, Adami C: The evolutionary origin of complex features. Nature. 2003, 423: 139-144. 10.1038/nature01568.View ArticlePubMedGoogle Scholar
- Misevic D, Ofria C, Lenski RE: Sexual reproduction reshapes the genetic architecture of digital organisms. Proc R Soc London B. 2006, 273: 457-464. 10.1098/rspb.2005.3338.View ArticleGoogle Scholar
- Misevic D, Ofria C, Lenski RE: Experiments with digital organisms on the origin and maintenance of sex in changing environments. J Hered. 2010, 101 (Suppl 1): S46-S54. 10.1093/jhered/esq017.View ArticlePubMedGoogle Scholar
- Goings S, Clune J, Ofria C, Pennock RT: Kin-selection: the rise and fall of kin-cheaters. Artificial Life IX. Edited by: Pollack, J., M. Bedau, P. Husbands, I. Ikegami, and R. A. Watson. 2004, Cambridge, MA: MIT Press, 303-308.Google Scholar
- Ofria C, Adami C, Collier TC: Selective pressures on genomes in molecular evolution. J Theor Biol. 2003, 222: 477-483.View ArticlePubMedGoogle Scholar
- Hagstrom GI, Hang DH, Ofria C, Torng E: Using Avida to test the effects of natural selection of phylogenetic reconstruction methods. Artificial Life. 2004, 10: 157-166. 10.1162/106454604773563586.View ArticlePubMedGoogle Scholar
- Ostrowski EA, Ofria C, Lenski RE: Ecological specialization and adaptive decay in digital organisms. Am Nat. 2007, 169: E1-E20. 10.1086/510211.View ArticlePubMedGoogle Scholar
- Yedid G, Ofria CA, Lenski RE: Historical and contingent factors affect re-evolution of a complex feature lost during mass extinction in communities of digital organisms. J Evol Biol. 2008, 21: 1335-1357. 10.1111/j.1420-9101.2008.01564.x.View ArticlePubMedGoogle Scholar
- Clune F, Misevic D, Ofria C, Lenski RE, Elena SF, Sanjuan R: Natural selection fails to optimize mutation rates for long-term adaptation on rugged fitness landscapes. PLoS Comput Biol. 2008, 4: e1000187-10.1371/journal.pcbi.1000187.PubMed CentralView ArticlePubMedGoogle Scholar
- Adami C: Ab initio modeling of ecosystems with artificial life. Nat Res Model. 2002, 15: 133-145. 10.1111/j.1939-7445.2002.tb00083.x.View ArticleGoogle Scholar
- Chow SS, Wilke CO, Ofria C, Lenski RE, Adami C: Adaptive radiation from resource competition in digital organisms. Science. 2004, 305: 84-86. 10.1126/science.1096307.View ArticlePubMedGoogle Scholar
- Cooper TF, Ofria C: Evolution of stable ecosystems in populations of digital organisms. Artificial Life VIII: Proceedings of the Eighth International Conference on Artificial Life, International Society for Artificial Life: 9-13 December 2002; Sydney, Australia. Edited by: Standish R, Bedau MA, Abbass HA. 2003, Cambridge, MA: MIT Press, 227-232.Google Scholar
- Elena SF, Wilke CO, Ofria C, Lenski RE: Effects of population size and mutation rate on the evolution of mutational robustness. Evolution. 2007, 61: 666-674. 10.1111/j.1558-5646.2007.00064.x.View ArticlePubMedGoogle Scholar
- Edlund JA, Adami C: Evolution of robustness in digital organisms. Artificial Life. 2004, 10: 167-179. 10.1162/106454604773563595.View ArticlePubMedGoogle Scholar
- Lenski RE, Ofria C, Collier TC, Adami C: Genome complexity, robustness and genetic interactions in digital organisms. Nature. 1999, 400: 661-664. 10.1038/23245.View ArticlePubMedGoogle Scholar
- Wilke CO, Wang JL, Ofria C, Lenski RE, Adami C: Evolution of digital organisms at high mutation rates leads to survival of the flattest. Nature. 2001, 412: 331-333. 10.1038/35085569.View ArticlePubMedGoogle Scholar
- Eyre-Walker A, Keightley PD: The distribution of fitness effects of new mutations. Nat Rev Genet. 2007, 8: 610-618. 10.1038/nrg2146.View ArticlePubMedGoogle Scholar
- Domingo-Calap P, Cuevas JM, Sanjuán R: The fitness effects of random mutations in single-stranded DNA and RNA bacteriophages. PLoS Genet. 2009, 5: e1000742-10.1371/journal.pgen.1000742.PubMed CentralView ArticlePubMedGoogle Scholar
- Muller HJ: The relation of recombination to mutational advance. Mutat Res. 1964, 1: 2-9.View ArticleGoogle Scholar
- Lynch M: The Origins of Genome Architecture. 2007, Sunderland, MA: Sinauer Associates, Inc. PublishersGoogle Scholar
- Ofria C, Brown CT, Adami C: The Avida user's manual. Introduction to Artificial Life. Edited by: Adami C. 1998, New York: Springer, 297-350.Google Scholar
- Ofria C, Wilke CO: Avida: a software platform for research in computational evolutionary biology. Artificial Life. 2004, 10: 191-229. 10.1162/106454604773563612.View ArticlePubMedGoogle Scholar
- ,: The Avida Digital Life Platform. accessed 27 January 2011, http://devolab.cse.msu.edu/software/avida/doc/http://devolab.cse.msu.edu/software/avida/doc/
- Dembski WA, Marks RJ: Conservation of information in search: measuring the cost of success. IEEE Trans Syst Man Cybern A Syst Hum. 2009, 5: 1051-1061. 10.1109/TSMCA.2009.2025027.View ArticleGoogle Scholar
- Cooper VS, Schneider D, Bolt M, Lenski RE: Mechanisms causing rapid and parallel losses of ribose catabolism in evolving populations of Escherichia coli B. J Bacteriol. 2001, 183: 2834-2841. 10.1128/JB.183.9.2834-2841.2001.PubMed CentralView ArticlePubMedGoogle Scholar
- Barrick JE, Yu DS, Yoon SH, Jeong H, Oh TK, Schneider D, Lenski RE, Kim JF: Genome evolution and adaptation in a long-term experiment with Escherichia coli. Science. 2009, 461: 1243-1247.Google Scholar
- Ewert W, Dembski WA, Marks RJ: Evolutionary synthesis of nand logic: dissecting a digital organism. Proceedings of the 2009 IEEE International Conference on Systems, Man, and Cybernetics: 11-14 October 2009. 2009, San Antonio, TX, 3047-3053.View ArticleGoogle Scholar
- Gerrish PJ, Lenski RE: The fate of competing beneficial mutations in an asexual population. Genetica. 1998, 102/103: 127-144. 10.1023/A:1017067816551.View ArticleGoogle Scholar
- Ohno S: Evolution By Gene Duplication. 1970, Springer-VerlagView ArticleGoogle Scholar
- Zhang J: Evolution by gene duplication: an update. Trends Ecol Evol. 2003, 18: 292-298. 10.1016/S0169-5347(03)00033-8.View ArticleGoogle Scholar
- Lenski RE, Barrick JE, Ofria C: Balancing robustness and evolvability. PLoS Biol. 2006, 4: e428-10.1371/journal.pbio.0040428.PubMed CentralView ArticlePubMedGoogle Scholar
- Mills DR, Peterson RL, Spiegelman S: An extracellular Darwinian experiment with a self-duplicating nucleic acid molecule. Biochemistry. 1967, 58: 217-224.Google Scholar
- Lenski RE, Travisano M: Dynamics of adaptation and diversification: A 10,000-generation experiment with bacterial populations. Proc Natl Acad Sci USA. 1994, 91: 6808-6814. 10.1073/pnas.91.15.6808.PubMed CentralView ArticlePubMedGoogle Scholar
- Cooper VS, Lenski RE: The population genetics of ecological specialization in evolving Escherichia coli populations. Nature. 2000, 407: 736-739. 10.1038/35037572.View ArticlePubMedGoogle Scholar
- Lenski RE, Rose MR, Simpson SC, Tadler SC: Long-term experimental evolution in Escherichia coli. I. Adaptation and divergence during 2,000 generations. Am Nat. 1991, 138: 1315-1341. 10.1086/285289.View ArticleGoogle Scholar
- Bull JJ, Badgett MR, Wichman HA: Big-benefit mutations in a bacteriophage inhibited with heat. Mol Biol Evol. 2000, 17: 942-950.View ArticlePubMedGoogle Scholar
- Rokyta DR, Beisel CJ, Joyce P, Ferris MT, Burch CL, Wichman HA: Beneficial effects are not exponential for two viruses. J Mol Evol. 2008, 67: 368-376. 10.1007/s00239-008-9153-x.PubMed CentralView ArticlePubMedGoogle Scholar
- Keightley PD, Eyre-Walker A: What can we learn about the distribution of fitness effects of new mutations from DNA sequence data?. Phil Trans R Soc B. 2010, 365: 1187-1193. 10.1098/rstb.2009.0266.PubMed CentralView ArticlePubMedGoogle Scholar
- Shaw RG, Byers DL, Darmo E: Spontaneous mutational effects on reproductive traits of Arabidopsis thaliana. Genetics. 2000, 155: 369-378.PubMed CentralPubMedGoogle Scholar
- Shaw FH, Geyer CJ, Shaw RG: A comprehensive model of mutations affecting fitness and inferences for Arabidopsis thaliana. Evolution. 2002, 56: 453-463.View ArticlePubMedGoogle Scholar
- Burke MK, Dunham JP, Shahrestani P, Thornton KR, Rose MR, Long AD: Genome-wide analysis of a long-term evolution experiment with Drosophila. Nature. 2010, 467: 587-590. 10.1038/nature09352.View ArticlePubMedGoogle Scholar
- Lind PA, Berg OG, Andersson DI: Mutational robustness of ribosomal proteins genes. Science. 2010, 330: 825-827. 10.1126/science.1194617.View ArticlePubMedGoogle Scholar
- Sanjuán R: Mutational fitness effects in RNA and single-stranded DNA viruses: common patterns revealed by site-directed mutagenesis studies. Phil Trans R Soc B. 2010, 365: 1975-1982.PubMed CentralView ArticlePubMedGoogle Scholar
- Peris JB, Davis P, Cuevas JM, Nebot MR, Sanjuán R: Distribution of fitness effects caused by single-nucleotide substitutions in bacteriophage f1. Genetics. 2010, 185: 603-609. 10.1534/genetics.110.115162.PubMed CentralView ArticlePubMedGoogle Scholar
- Baba T, Ara T, Hasegawa M, Takai Y, Okumura Y, Baba M, Datsenko KA, Tomita M, Wanner BL, Mori H: Construction of Escherichia coli K-12 in-frame, single-gene knockout mutants: the Keio collection. Mol Syst Biol. 2006, 2: 2006-10.1038/msb4100050.View ArticleGoogle Scholar
- Elena SF, Ekunwe L, Hajela N, Oden SA, Lenski RE: Distribution of fitness effects caused by random insertion mutations in Escherichia coli. Genetica. 1998, 102/103: 349-358. 10.1023/A:1017031008316.View ArticleGoogle Scholar
- Sanjuán R, Moya A, Elena SF: The distribution of fitness effects caused by single-nucleotide substitutions in an RNA virus. Proc Natl Acad Sci USA. 2004, 101: 8396-8401.PubMed CentralView ArticlePubMedGoogle Scholar
- Elena SF, Moya A: Rate of deleterious mutation and the distribution of its effects on fitness in vesicular stomatitis virus. J Evol Biol. 1999, 12: 1078-1088. 10.1046/j.1420-9101.1999.00110.x.View ArticleGoogle Scholar
- Kingsolver JG, Hoekstra HE, Hoekstra JM, Berrigan D, Vignieri SN, Hill CE, Hoang A, Gibert P, Beerli P: The strength of phenotypic selection in natural populations. Am Nat. 2001, 157: 245-261. 10.1086/319193.View ArticlePubMedGoogle Scholar
- Lynch M, Gabriel W: Mutation load and the survival of small populations. Evolution. 1990, 44: 1725-1737. 10.2307/2409502.View ArticleGoogle Scholar
- Gabriel W, Lynch M, Bürger R: Muller's ratchet and mutational meltdowns. Evolution. 1993, 47: 1744-1757. 10.2307/2410218.View ArticleGoogle Scholar
- Lynch M, Bürger R, Butcher D, Gabriel W: The mutational meltdown in asexual populations. J Hered. 1993, 84: 339-344.PubMedGoogle Scholar
- Butcher D: Muller's ratchet, epistasis and mutation effects. Genetics. 1995, 141: 431-437.PubMed CentralPubMedGoogle Scholar
- Bull JJ, Sanjuán R, Wilke CO: Theory of lethal mutagenesis for viruses. J Virol. 2007, 81: 2930-2939. 10.1128/JVI.01624-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Lande R: Risk of population extinction from fixation of new deleterious mutations. Evolution. 1994, 48: 1460-1469. 10.2307/2410240.View ArticleGoogle Scholar
- Lynch M, Conery J, Bürger R: Mutation accumulation and the extinction of small populations. Am Nat. 1995, 146: 489-518. 10.1086/285812.View ArticleGoogle Scholar
- Muller HJ: Our load of mutations. Am J Hum Genet. 1950, 2: 111-176.PubMed CentralPubMedGoogle Scholar
- Kondrashov AS: Contamination of the genome by very slightly deleterious mutations: why have we not died 100 times over?. J Theor Biol. 1995, 175: 583-594. 10.1006/jtbi.1995.0167.View ArticlePubMedGoogle Scholar
- Crow JF: The high spontaneous mutation rate: is it a health risk?. Proc Natl Acad Sci USA. 1997, 94: 8380-8386. 10.1073/pnas.94.16.8380.PubMed CentralView ArticlePubMedGoogle Scholar
- Eyre-Walker A, Keightley PD: High genomic deleterious mutation rates in Hominids. Nature. 1999, 397: 344-347. 10.1038/16915.View ArticlePubMedGoogle Scholar
- Eyre-Walker A, Woolfit M, Phelps T: The distribution of fitness effects of new deleterious amino acid mutations in humans. Genetics. 2006, 173: 891-900. 10.1534/genetics.106.057570.PubMed CentralView ArticlePubMedGoogle Scholar
- Lynch M: Rate, molecular spectrum, and consequences of human mutation. Proc Natl Acad Sci USA. 2010, 107: 961-968. 10.1073/pnas.0912629107.PubMed CentralView ArticlePubMedGoogle Scholar
- Estes S, Lynch M: Rapid fitness recovery in mutationally degraded lines of Caenorhabditis elegans. Evolution. 2003, 57: 1022-1030.View ArticlePubMedGoogle Scholar
- Martin G, Gandon S: Lethal mutagenesis and evolutionary epidemiology. Phil Trans R Soc B. 2010, 365: 1953-1963. 10.1098/rstb.2010.0058.PubMed CentralView ArticlePubMedGoogle Scholar
- Galvani AP: The role of mutation accumulation in HIV progression. Proc R Soc B. 2005, 272: 1851-1858. 10.1098/rspb.2005.3083.PubMed CentralView ArticlePubMedGoogle Scholar
- ,: Digital Evolution Laboratory. accessed 27 January 2011, http://devolab.msu.edu/http://devolab.msu.edu/
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.