The effect of genetic robustness on evolvability in digital organisms
© Elena and Sanjuán; licensee BioMed Central Ltd. 2008
Received: 15 April 2008
Accepted: 14 October 2008
Published: 14 October 2008
Recent work has revealed that many biological systems keep functioning in the face of mutations and therefore can be considered genetically robust. However, several issues related to robustness remain poorly understood, such as its implications for evolvability (the ability to produce adaptive evolutionary innovations).
Here, we use the Avida digital evolution platform to explore the effects of genetic robustness on evolvability. First, we obtained digital organisms with varying levels of robustness by evolving them under combinations of mutation rates and population sizes previously shown to select for different levels of robustness. Then, we assessed the ability of these organisms to adapt to novel environments in a variety of experimental conditions. The data consistently support that, for simple environments, genetic robustness fosters long-term evolvability, whereas, in the short-term, robustness is not beneficial for evolvability but may even be a counterproductive trait. For more complex environments, however, results are less conclusive.
The finding that the effect of robustness on evolvability is time-dependent is compatible with previous results obtained using RNA folding algorithms and transcriptional regulation models. A likely scenario is that, in the short-term, genetic robustness hampers evolvability because it reduces the intensity of selection, but that, in the long-term, relaxed selection facilitates the accumulation of genetic diversity and thus, promotes evolutionary innovation.
Evidence has accumulated during recent years suggesting that organisms can maintain their performance in the face of a broad range of perturbations [1, 2]. This includes the tolerance of proteins to amino acid replacements , the ability of genetic networks to withstand alterations to its parts , the stability of cellular processes to the stochastic variations of gene expression levels , or the resilience of embryonic development to environmental or genetic changes . In general, the term 'robustness' is used to describe this kind of behavior and, when mutations are the cause of perturbation, the term 'genetic robustness' or 'mutational robustness' is used. Many issues related to genetic robustness remain unresolved. For example, asserting that elevated robustness is a fundamental property of living organisms is problematic because we often ignore what normal robustness should be . Still, we can try to identify the genetic and ecological factors associated with differences in robustness between species or genotypes [8, 9]. Also, it remains unclear whether the evolutionary transition to a robust state occurs as a direct product of selection [10–13] or merely as a by-product of selection acting on correlated traits [14–16].
A system is said to be evolvable if it can be modified through genetic change in a way that enhances survival and reproduction. For natural selection to act, the system must show heritable phenotypic variation. Yet, genetic robustness implies that the system produces little phenotypic variation in response to genetic variation. Therefore, robustness might limit evolutionary optimization and innovation . In this vein, theoretical work has postulated that buffering mechanisms can lead to maladaptation compared to what would be achieved in the absence of these mechanisms . Also, the analysis of gene expression noise in yeast suggests that noise control may indirectly increase mutational robustness, which might in turn hamper evolvability at the level of gene expression . On the other hand, genetic robustness facilitates the accumulation of neutral or nearly neutral variation by relaxing the intensity of natural selection. This accumulated diversity can become visible to selection upon changes in the environment or the genetic background, and thus be the source of evolutionary innovation. Computer simulations on simple population genetics models predict that genetic robustness can sometimes facilitate access to new adaptive peaks provided that occasional failures of robustness mechanisms occur . The view that robustness can foster evolvability has also been supported by lattice protein models and PCR-based mutagenesis experiments showing that protein variants with increased thermodynamic stability have increased genetic robustness and are more likely to evolve new catalytic capabilities .
Here, we address the relationship between robustness and evolvability by directly comparing the ability of genotypes with different degrees of robustness to adapt to novel environments. At present, it is difficult to envisage such experiments with natural organisms. Current work in this area is mainly focused on the characterization of robustness in different species or genotypes from the same species [22–24]. Instead, we perform our experiments using digital organisms. Digital organisms are self-replicating computer programs that inhabit a virtual world where they reproduce, mutate, compete for resources, and evolve . Since digital and natural organisms evolve under the same fundamental rules, the evolution of digital organisms should be informative about the natural world. It must be noted, however, that the particular physicochemical properties of natural systems are not always captured in silico, despite the potentially important role of these properties in evolution. On the other hand, the use of digital organisms offers several advantages: first, experiments can be performed on scales that are beyond reach with any biological entity; second, the study of evolution with digital organisms allows one to collect extremely precise information about the evolutionary process; third, one can easily perform certain genetic manipulations that would be exceedingly laborious on natural organisms.
Using the Avida platform , we first obtained digital organisms with varying levels of genetic robustness, as previously described . We then evaluated the ability of these organisms to adapt to new environments. The results suggest that, in a simple environment, genetic robustness retards the first adaptive steps but clearly fosters long-term evolvability. Consequently, at any given time point, whether robustness promotes evolvability depends on the rate of adaptation and thus, indirectly, on the mutation rate, but also on the complexity of the environment. We argue that, in the short-term, genetic robustness can slow down adaptive evolution by relaxing the intensity natural selection. In the long-term, however, it can favor the accumulation of genetic variation and thus provide access to more of the fitness landscape.
Preliminary experiment: evolution of genotypes with different levels of genetic robustness
We evolved the Avida default organism under different regimes previously shown to favor the divergence of robustness . Nine regimes, resulting from the combination of three genomic mutation rates (0.02, 0.2 and 2) and three population sizes (102, 103 and 5 × 104) were tested, and 20 independent lineages were seeded for each regime, yielding 180 lineages in total. Organisms were maintained for 105 updates in an environment in which four logic tasks (the digital equivalent to metabolic resources) were rewarded (NOT, AND, OR, and NOR). At the end of the evolution experiment, the most frequent organism was isolated from each lineage and its adaptation to the new environment was judged from its ability to perform these four tasks, a requirement which was met by 62 of the 180 evolved organisms. Organisms not performing the four tasks were discarded from subsequent analyses. The use of organisms homogeneous with respect to their ability to perform tasks ensured that none of them would be a priori better predisposed than others to learn new tasks.
The robustness of these organisms was computed and an analysis of variance showed that, consistent with previous work [27, 28], high mutation rates tended to favor the evolution of robustness (P < 0.001), whereas population size had no direct effect on robustness (P = 0.818), but and indirect one conditional on mutation rate (interaction term, P < 0.001). We also observed a weak negative correlation between robustness and log fitness (Spearman's ρ = -0.270, 60 d.f., P = 0.034), indicating that robustness paid some fitness cost, an observation which is also in line with theoretical considerations  and previous experiments [22, 24, 28]. As expected, robustness positively correlated with the fraction of selectively neutral mutations (ρ = 0.899, 60 d.f., P < 0.001) and negatively correlated with the magnitude of deleterious fitness effects (ρ = -0.376, 60 d.f., P = 0.003). Interestingly, the fraction of beneficial mutations also increased with robustness (ρ = 0.626, 60 d.f., P < 0.001). This later association was not a trivial consequence of the fact that more robust genotypes were less fit, because the correlation held after controlling for log fitness (partial correlation: r = 0.564, 59 d.f., P < 0.001). Finally, the magnitude of beneficial mutations did not depend on robustness after controlling for fitness (partial correlation: r = -0.140, 59 d.f., P = 0.282).
Analysis of the genotypes showing the largest and smallest robustness values
The two organisms with the most extreme genetic robustness values obtained from the preliminary experiment
5 × 104
5 × 104
We placed R and F in a novel environment in which four new tasks (NAND, ORN, ANDN, and XOR) were rewarded in addition to the four original ones (for a total of 8 tasks rewarded). Differences in mutational effects held in the new environment (F: 80.84% deleterious, 18.84% neutral, 0.32% beneficial; R: 52.84% deleterious, 40.52% neutral, and 6.64% beneficial; χ2 = 20.913, 2 d.f., P < 0.001). We seeded 50 independent lineages of each R and F and let them evolve for 1000 updates at a population size of 104 (this size was used in all subsequent experiments). Given that F and R had been previously evolved at different mutation rates (Table 1), we explored various genomic mutation rates (U), ranging 0.03 – 3. Also, to evaluate the stability of robustness, we recalculated it at the end of the run for each of the 50 populations evolved at U = 0.3. Even though genotypes evolved from R lost 24.38% of their robustness and those evolved from F increased in robustness by 20.09%, differences between the two groups were still highly significant (Wilcoxon's signed ranks test: P < 0.001).
General correlations between robustness and evolvability
In the 8-task environment, the results are in line with the experiments performed using the (R, F) pair. For U = 0.03, the correlation between robustness and evolvability remained negative during approximately the first 2500 updates and became positive from this point onwards. After update 11,500 the correlation reached a steady-state median value of 0.172, with fluctuations in the range 0.117 – 0.220 (Fig. 3). For U = 0.3 the correlation was null or slightly negative before update 2500, and then became positive, reaching its maximum value (0.263) around update 104. From this time point on, the correlation declined until reaching a plateau value around update 15,000, and fluctuated afterwards around a median value of 0.109 (range 0.075 – 0.154). Finally, for U = 3, the advantage of robustness was evident right from the beginning, as supported by the positive correlation between robustness and evolvability observed along the entire experiment. In this case, the correlation fluctuated in the range 0.011 – 0.216, with a median value of 0.104. Altogether, these data are compatible with our above conclusion that robustness fosters evolvability in the long-term. The distinction between long- and short-term logically depends on the rate of adaptation which, in turn depends on the rate of mutation (see above). This is why the higher the mutation rate, the earlier the payoffs of robustness. For the three mutation rates explored, the correlation between robustness and evolvability converged to a positive value as time increased, although the exact values of these plateaus seemed to decrease as mutation rate increases.
In the more complex 77-task environment, the results were more difficult to interpret. For U = 0.03, the correlation between robustness and evolvability was always negative (median – 0.079; range – 0.159 to – 0.004) throughout the 5 × 104 updates of evolution, similar to the results obtained using the (R, F) pair. It is possible that the time required for robustness to confer a fitness advantage was longer than 5 × 104 updates, but after such long time the initial differences in robustness may have vanished or, alternatively, populations might have already been confined to a stable fitness value. For U = 0.3, robustness fostered evolvability from approximately update 5500 on, similar to what was observed in the 8-task environment. However, the benefit of robustness declined with evolutionary time, and the correlation tended to converge to zero at the end of the experiment. For U = 3, the correlation was always negative (median value of – 0.121; range – 0.017 to – 0.176), a result which was not consistent with the pattern observed in the 8-task environment.
Using the Avida digital evolution platform , we have investigated the effects of genetic robustness on evolvability. We first obtained organisms with varying levels of robustness by evolving a common ancestor at different combinations of mutation rates and population sizes. As previously shown [28, 29], high mutation rates directly favored the evolution of robustness, whereas the effect of population size was only evident through its interplay with mutation rate. Differences in robustness evolved from populations that, apart from mutation rate, were subjected to exactly the same conditions, indicating that high mutation rates can directly select for robustness without the contribution of additional factors.
Our analyses indicate that the effect of robustness on evolvability is time-dependent, and probably also dependent on the environment and the mutation rate. In the simple environment, we observed that robustness tended to hinder evolvability in the short-term. The distinction between short- and long-term does not refer to generations but instead to stages of adaptation. As such, the effect of robustness indirectly depends on the mutation rate, since the latter affects the rate of adaptation. For low mutation rates, adaptation was slow and hence, the short-term cost of robustness endured for a relatively large number of generations. In contrast, for high mutation rates, the cost was brief or absent. We also observed that robustness consistently fostered long-term adaptation in the simple environment, a result which was weakly dependent on the mutation rate. However, for the complex environment, such a long-term benefit was less clear. The analysis of the (R, F) pair suggested that the results in the complex environment paralleled those in the simpler one, although long-term benefits occurred after larger numbers of generations. However, the analysis of the 62 genotypes obtained from the preliminary experiment indicated that the correlation between robustness and evolvability could be either positive or negative within the timeframe studied. It is possible that longer evolutionary times were required to observe a sustained positive association between robustness and evolvability. However, the genetic inheritance of robustness might be lost after such long times, due to genetic drift or to competition with fitter, less robust, genotypes.
It is noteworthy that the fraction of beneficial mutations increased with robustness, at least in the simple environment, whereas the magnitude of beneficial mutations did not significantly change. Looking at these data, one may conclude that robustness should invariantly promote evolvability, simply because beneficial mutations are more readily found. Furthermore, if beneficial mutations are rare, i.e. for short evolutionary times or for low mutation rates, the rate of adaptation ought to be mainly determined by the time required to find beneficial mutations , but this contradicts our observation that robust genotypes were less evolvable in the short-term, or at low mutation rates. To explain this apparent contradiction, deleterious mutations have to be incorporated into the picture. Recall that in the robust genotypes the intensity of deleterious mutational effects decreased. Given that the magnitude of fitness effects is higher on average for deleterious mutations than for beneficial ones , and given that the digital organisms used in our experiment reproduced asexually, the chance that a beneficial mutation was successful critically depended on whether it first appeared on a genetic background free of deleterious mutations. Owing to the weaker intensity of selection, deleterious mutations are removed less efficiently from a population of genetically robust genotypes  and therefore, the rate at which beneficial mutations are fixed can be slowed down in a robust genotypic background, even if such mutations appear more frequently.
Recent work has attempted to resolve the apparently paradoxical relationship between robustness and evolvability . To do that, a distinction between genotypic and phenotypic robustness has been proposed. The former simply refers to the standard definition of genetic robustness and is the one used in this article, whereas the latter refers to the average genetic robustness of all genotypes encoding a phenotype. Hence, whereas genotypic robustness applies to individuals, phenotypic robustness concerns sets of neutral genotypes. If evolvability is defined as the probability of finding novel functions in the immediate mutational neighborhood (i.e. innovations involving a single mutation), genotypic robustness necessarily hampers evolvability . However, phenotypic robustness would tend to promote evolvability because this kind of robustness should be associated with larger neutral network sizes, which would allow populations to explore vaster regions of the fitness landscape [33, 34]. The distinction between genotypic and phenotypic robustness parallels with the distinction between short-term and long-term evolution. A founder genotype with elevated genetic robustness will perform its first adaptive steps at a slower path than a less robust one. However, as generations go on, it will degenerate on a population of genetically more heterogeneous individuals. Hence, despite decreasing phenotypic variability, robustness might foster long-term evolvability.
Two conditions have been identified as required for a system to be robust and evolvable : first, genotypes need to be integrated in a highly connected neutral network, thereby implying robustness. It has been suggested that the connectedness of networks might be a very general property of evolving systems, since it has been observed independently in different processes as transcriptional regulation and RNA folding . Second, neutral networks need to span a large fraction of the total genotypic space. Without this second condition, the system would be poorly evolvable, because only a small region of the fitness landscape could be explored throughout neutral mutations. In other words, the 'memory' of past robustness would be rapidly lost along evolutionary paths. This seems to be the case of RNA secondary structures, for which the effects of robustness on evolvability can only be predicted in the very short-term . However, simulations suggest that neutral networks could be much more pervasive in the case of transcriptional regulation . In Avida, robustness seems to be relatively stable, as suggested by the fact that initial differences in robustness were maintained after thousands of generations. In a simple environment, this stability was sufficient for the benefits of increased evolvability to be paid off before non-robust types could invade. However, we could not conclude whether the same would be true in a more complex environment. Future work will elucidate whether this condition is attainable for other evolving systems and, therefore, whether robustness can generally promote long-term evolvability.
The effect of robustness on evolvability in a simple environment is time-dependent. Robustness tends to retard the first adaptive steps, whereas it fosters long-term adaptation. Furthermore, this effect also depends on mutation rate: at low mutation rates, relatively large numbers of generations are required to overcome the short-term cost of robustness. In contrast, for high mutation rates, this cost can be brief or absent. A possible explanation for these results is that deleterious mutations are removed less efficiently from a population of genetically robust genotypes, therefore interfering with the fixation of beneficial mutations. However, in the long-term, robust populations would accumulate higher levels of genetic variation and hence would be able to explore vaster regions of the fitness landscape. The complexity of the environment also plays a role in the relationship between robustness and evolvability: in a more complex environment, the correlation between robustness and evolvability could be either positive or negative within the timeframe studied.
The Avida platform
The Avida platform , was used for all experiments. In short, the genome of the starting digital organism consists of 100 instructions, with 26 possible values for each instruction . Genomes are copied line by line by repeatedly executing the "copy" instruction. A subsequent execution of the "divide" instruction produces two independent genomes. Each new offspring replaces a random organism from the population, and organisms die when other organisms replace them or if 2000 instructions are executed before completing replication.
Each organism's phenotype depends upon the complex rules governing how the instructions encoded in its genome are executed. CPU time needed to execute genomic instructions is the basic resource organisms compete for. Access to CPU time depends on the metabolic rate, which was initially set to the arbitrary value 100 for each individual, but varies as organisms mutate. The investigator defines an environment by choosing a series of rewarded logical operations (e.g. AND, OR, NAND) that represent available resources. Organisms that spontaneously acquire the ability to perform these operations increase their own metabolic rate (i.e. can use more CPU time) and thus tend to increase their own replication rate. In this study, all logical resources were equally rewarded. It is important to clarify that replication rates are not user-defined but instead result from the expression of genotypes in the environment. Finally, the realized fitness of each individual depends not only on its own replication rate, but also on that of the other members of the population . As in biological systems, the ability to perform tasks evolves by spontaneous mutation, selection, and drift.
Each run was initiated with a population of identical individuals. The population size was chosen by the investigator and remained constant throughout the evolution experiment (hereafter, a run). Updates are the normal time unit in Avida and represent the average time required to execute 30 instructions . Here, an update corresponded to roughly 0.1 generations. Replication is subject to a fraction of point mutations, with a probability of μ per copied instruction or U per genome (U = 100 μ). The mutation rate was user-defined and constant throughout each run. In all cases, genome lengths were kept constant at 100 instructions and reproduction was strictly asexual. In any given experimental treatment, each replicate run had identical initial conditions except for a random-number seed. This initial seed causes runs to differ at all subsequent points where stochastic events occur, including mutations and the physical location where each new offspring is placed.
Experiments were performed using version 2.7.0 of Avida. This software is freely available from devolab.cse.msu.edu/software/avida. Specific details about Avida can be found elsewhere  and at the above website. Runs were performed using the Linux operating system on a cluster of 160 AMD Athlon 1600+ processors.
Measures of genetic robustness and evolvability
For any given population, we first identified the most abundant genotype and then systematically constructed every possible single mutation on this genotype and tested its fitness effects, as previously described [27, 28]. The signed selection coefficient of mutation i was expressed as s i = W i /W0 - 1, where W i is the fitness of the ith mutant and W0 the fitness of the non-mutated organism. If a genotype has increased robustness, then mutations should, on average, have a little or no impact on fitness and thus, the average selection coefficient should be close to zero. By contrast, the less robust a genotype is, the more negative the average effect across all possible mutations. Genetic robustness thus can be measured as the signed average selection coefficient, excluding any possible beneficial mutations, . Finally, notice that our estimates of the average selection coefficient were obtained without error, since the fitness consequence of every possible single mutation was evaluated.
To measure evolvability we placed genotypes in a novel environment, defined by a combination of rewarded tasks (metabolic resources). The genotype was then let to evolve at a population size of 104. This population size was chosen because it is large enough to allow selection to prevail over genetic drift while not compromising computational efficiency. Evolvability was measured as the observed increase in log fitness relative to the ancestor.
We thank Richard E. Lenski and Jeff Clune for helpful comments and Javier Forment for technical assistance. This work was supported by grants BFU2006-14819-C02-01/BMC and BFU2008-03978/BMC, and the Ramón y Cajal program from the Spanish MICINN.
- de Visser JA, Hermisson J, Wagner GP, Ancel ML, Bagheri-Chaichian H, Blanchard JL, Chao L, Cheverud JM, Elena SF, Fontana W, et al: Perspective: Evolution and detection of genetic robustness. Evolution. 2003, 57: 1959-1972.PubMedGoogle Scholar
- Wagner A: Robustness and evolvability in living systems. 2005, Princeton University Press, New JerseyGoogle Scholar
- Sinha N, Nussinov R: Point mutations and sequence variability in proteins: redistributions of preexisting populations. Proc Natl Acad Sci USA. 2001, 98: 3139-3144. 10.1073/pnas.051399098.PubMed CentralView ArticlePubMedGoogle Scholar
- Aldana M, Balleza E, Kauffman S, Resendiz O: Robustness and evolvability in genetic regulatory networks. J Theor Biol. 2007, 245: 433-448. 10.1016/j.jtbi.2006.10.027.View ArticlePubMedGoogle Scholar
- Batada NN, Hurst LD: Evolution of chromosome organization driven by selection for reduced gene expression noise. Nat Genet. 2007, 39: 945-949. 10.1038/ng2071.View ArticlePubMedGoogle Scholar
- von Dassow G, Meir E, Munro EM, Odell GM: The segment polarity network is a robust developmental module. Nature. 2000, 406: 188-192. 10.1038/35018085.View ArticlePubMedGoogle Scholar
- Ciliberti S, Martin OC, Wagner A: Robustness can evolve gradually in complex regulatory gene networks with varying topology. PLoS Comput Biol. 2007, 3: e15-10.1371/journal.pcbi.0030015.PubMed CentralView ArticlePubMedGoogle Scholar
- Krakauer DC, Plotkin JB: Redundancy, antiredundancy, and the robustness of genomes. Proc Natl Acad Sci USA. 2002, 99: 1405-1409. 10.1073/pnas.032668599.PubMed CentralView ArticlePubMedGoogle Scholar
- Sanjuán R, Elena SF: Epistasis correlates to genomic complexity. Proc Natl Acad Sci USA. 2006, 103: 14402-14405. 10.1073/pnas.0604543103.PubMed CentralView ArticlePubMedGoogle Scholar
- Proulx SR: The opportunity for canalization and the evolution of genetic networks. Am Nat. 2005, 165: 147-162. 10.1086/426873.View ArticlePubMedGoogle Scholar
- Proulx SR, Nuzhdin S, Promislow DE: Direct selection on genetic robustness revealed in the yeast transcriptome. PLoS ONE. 2007, 2: e911-10.1371/journal.pone.0000911.PubMed CentralView ArticlePubMedGoogle Scholar
- van Nimwegen E, Crutchfield JP, Huynen M: Neutral evolution of mutational robustness. Proc Natl Acad Sci USA. 1999, 96: 9716-9720. 10.1073/pnas.96.17.9716.PubMed CentralView ArticlePubMedGoogle Scholar
- Wagner GP, Booth G, Bagheri-Chaichian H: A population genetic theory of canalization. Evolution. 1997, 51: 329-347. 10.2307/2411105.View ArticleGoogle Scholar
- Ancel LW, Fontana W: Plasticity, evolvability, and modularity in RNA. J Exp Zool. 2000, 288: 242-283. 10.1002/1097-010X(20001015)288:3<242::AID-JEZ5>3.0.CO;2-O.View ArticlePubMedGoogle Scholar
- Stearns SC, Kaiser M, Kawecki TJ: The differentical genetic and environmental canalization of fitness components in Drosophila melanogaster. J Evol Biol. 1995, 8: 539-557. 10.1046/j.1420-9101.1995.8050539.x.View ArticleGoogle Scholar
- Stearns SC: Progress on canalization. Proc Natl Acad Sci USA. 2002, 99: 10229-10230. 10.1073/pnas.172388999.PubMed CentralView ArticlePubMedGoogle 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
- Frank SA: Maladaptation and the paradox of robustness in evolution. PLoS ONE. 2007, 2: e1021-10.1371/journal.pone.0001021.PubMed CentralView ArticlePubMedGoogle Scholar
- Lehner B: Selection to minimise noise in living systems and its implications for the evolution of gene expression. Mol Syst Biol. 2008, 4: 170-10.1038/msb.2008.11.PubMed CentralView ArticlePubMedGoogle Scholar
- Kim Y: Rate of adaptive peak shifts with partial genetic robustness. Evolution. 2007, 61: 1847-1856. 10.1111/j.1558-5646.2007.00166.x.View ArticlePubMedGoogle Scholar
- Bloom JD, Labthavikul ST, Otey CR, Arnold FH: Protein stability promotes evolvability. Proc Natl Acad Sci USA. 2006, 103: 5869-5874. 10.1073/pnas.0510098103.PubMed CentralView ArticlePubMedGoogle Scholar
- Codoñer FM, Daros JA, Sole RV, Elena SF: The fittest versus the flattest: experimental confirmation of the quasispecies effect with subviral pathogens. PLoS Pathog. 2006, 2: e136-10.1371/journal.ppat.0020136.PubMed CentralView ArticlePubMedGoogle Scholar
- Montville R, Froissart R, Remold SK, Tenaillon O, Turner PE: Evolution of mutational robustness in an RNA virus. PLoS Biol. 2005, 3: e381-10.1371/journal.pbio.0030381.PubMed CentralView ArticlePubMedGoogle Scholar
- Sanjuán R, Cuevas JM, Furio V, Holmes EC, Moya A: Selection for robustness in mutagenized RNA viruses. PLoS Genet. 2007, 3: e93-10.1371/journal.pgen.0030093.PubMed CentralView ArticlePubMedGoogle Scholar
- Ray TS: An approach to the synthesis of life. Artificial Life II. Edited by: Langton CG, Taylor C, Farmer JD, Rasmussen S. 1991, Redwood City, CA: Addison-WesleyGoogle Scholar
- Ofria C, Wilke CO: Avida: a software platform for research in computational evolutionary biology. Artif Life. 2004, 10: 191-229. 10.1162/106454604773563612.View ArticlePubMedGoogle 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
- 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
- Elena SF, Sanjuán R: Virus Evolution: Insights From an Experimental Approach. Annu Rev Ecol Sys Evol. 2007, 38: 27-52. 10.1146/annurev.ecolsys.38.091206.095637.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 ArticlePubMedGoogle Scholar
- Orr HA: The genetic theory of adaptation: a brief history. Nat Rev Genet. 2005, 6: 119-127. 10.1038/nrg1523.View ArticlePubMedGoogle Scholar
- Orr HA: The rate of adaptation in asexuals. Genetics. 2000, 155: 961-968.PubMed CentralPubMedGoogle Scholar
- Wagner A: Robustness and evolvability: a paradox resolved. Proc Biol Sci. 2008, 275: 91-100. 10.1098/rspb.2007.1137.PubMed CentralView ArticlePubMedGoogle Scholar
- Ciliberti S, Martin OC, Wagner A: Innovation and robustness in complex regulatory gene networks. Proc Natl Acad Sci USA. 2007, 104: 13591-13596. 10.1073/pnas.0705396104.PubMed CentralView 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
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.