The relationship between the error catastrophe, survival of the flattest, and natural selection
© Tejero et al; licensee BioMed Central Ltd. 2011
Received: 8 July 2010
Accepted: 4 January 2011
Published: 4 January 2011
The quasispecies model is a general model of evolution that is generally applicable to replication up to high mutation rates. It predicts that at a sufficiently high mutation rate, quasispecies with higher mutational robustness can displace quasispecies with higher replicative capacity, a phenomenon called "survival of the flattest". In some fitness landscapes it also predicts the existence of a maximum mutation rate, called the error threshold, beyond which the quasispecies enters into error catastrophe, losing its genetic information. The aim of this paper is to study the relationship between survival of the flattest and the transition to error catastrophe, as well as the connection between these concepts and natural selection.
By means of a very simplified model, we show that the transition to an error catastrophe corresponds to a value of zero for the selective coefficient of the mutant phenotype with respect to the master phenotype, indicating that transition to the error catastrophe is in this case similar to the selection of a more robust species. This correspondence has been confirmed by considering a single-peak landscape in which sequences are grouped with respect to their Hamming distant from the master sequence. When the robustness of a classe is changed by modification of its quality factor, the distribution of the population changes in accordance with the new value of the robustness, although an error catastrophe can be detected at the same values as in the general case. When two quasispecies of different robustness competes with one another, the entry of one of them into error catastrophe causes displacement of the other, because of the greater robustness of the former. Previous works are explicitly reinterpreted in the light of the results obtained in this paper.
The main conclusion of this paper is that the entry into error catastrophe is a specific case of survival of the flattest acting on phenotypes that differ in the trade-off between replicative ability and mutational robustness. In fact, entry into error catastrophe occurs when the mutant phenotype acquires a selective advantage over the master phenotype. As both entry into error catastrophe and survival of the flattest are caused by natural selection when mutation rate is increased, we propose differentiating between them by the level of selection at which natural selection acts. So we propose to consider the transition to error catastrophe as a phenomenon of intra-quasispecies selection, and survival of the flattest as a phenomenon of inter-quasispecies selection.
The quasispecies model was developed by Eigen in 1971 , and has been applied in many different fields, on account of its usefulness as a general evolutionary model for error-prone self-replicative systems. Prebiotic self-replicating molecules , RNA viruses , cancer cells , the immune system , etc. have all been modelled as quasispecies. One of the most important implications of the quasispecies theory, of particular relevance because of its originality and possibility for practical applications, is the concept of error catastrophe . In some fitness landscapes, an increase in the mutation rate beyond an error threshold causes the quasispecies to enter into error catastrophe. In the classical model developed by Eigen, several different phenomena can be observed when the population crosses the error threshold. First, the master sequence, i.e. the sequence with the highest replicative ability, is lost . Second, the quasispecies gets delocalized over the whole sequence space , in such a way that sequences become uniformly distributed. As this is materially impossible in finite populations, this delocalization should be interpreted as the drift of the population over the sequence space . Finally, several abrupt changes can be observed for different population traits, namely in the mean and the variance of the Hamming distance from the master sequence [9, 10], in some collective order parameters such as the consensus sequence  or in the so-called ancestral distribution [12, 13], etc. Although in the single-peak fitness landscape used by Eigen, all these phenomena take place at the same value of q, this is not necessarily the case for other, more complicated, fitness landscapes . As Hermisson and co-workers  have shown, these phenomena are different kinds of what they call "mutation thresholds". These thresholds were defined as "pronounced changes of the equilibrium distribution of some population trait or fitness values within a narrow change of mutation rates" . They identify four different thresholds: a, a "fitness threshold", which is the mean fitness of the population that suffers the pronounced change; b, a "wild-type threshold", which is defined as the loss of the fittest sequence; c, a "degradation threshold", beyond which fitness no longer changes with the mutation rate; and d, a "trait threshold", defined by a pronounced change in a population trait. Moreover, it has been shown that all these thresholds are completely dependent on the fitness and the mutation landscape [13, 15]. One problem derived from this fact is that the definitions of the error catastrophe, and of the error threshold, change from paper to paper [13, 14, 16]. In most papers the error catastrophe is equated to the loss of the master sequence [1, 17–19], in other papers it is taken to be the delocalization of the population over the sequence space [20–23], while in others it is defined as a transition to more robust regions of the sequence space . Actually, in some papers two of these definitions are used simultaneously. However, the usual mathematical criterion for defining the error threshold is the loss of master sequence in the absence of back mutation. The differences that are observed for the phenomenology and definition of error catastrophe, and error threshold, can also be observed in the interpretation of both concepts. Error threshold is said to be a critical mutation rate beyond which one of the following phenomena take place: mutation dominates over selection [13, 25, 26], natural selection ceases to operate [5, 13], there is no mutation-selection balance [18, 27], or evolutionary adaptation  or optimization  breaks down.
Furthermore, both the loss of the master sequence and the delocalization of the quasispecies over the sequence space have been related to an information "crisis" or "meltdown" [6, 16, 26, 29, 30]. In this sense, the entry into error catastrophe is supposed to establish a maximum limit on the information content that a self-replicative system can maintain at a given mutation rate [20, 31, 32]. This interpretation has two important practical consequences. In the field of the origin of life, it introduces what is known as "Eigen's paradox" [1, 29, 33], according to which the first self-replicating molecules would not be long enough, at prebiotic mutation rates, to encode the enzymes, or functions, required to copy the sequences more accurately. The second important consequence is in the field of RNA viruses [16, 34], since the possibility of pushing RNA virus replication into error catastrophe by means of mutagenic drugs was the origin of the first "lethal mutagenesis" experiments, i.e. virus extinction through increased mutagenesis, as well as a first explanation for the loss of viral infectivity [35–38].
Another implication of the quasispecies theory that we would like to emphasize is the "survival of the flattest" concept [39–42]. Basically, the "survival of the flattest" effect postulates that at high mutation rates a quasispecies with high reproductive capacity, but low mutational robustness can be displaced by another quasispecies with lower reproductive capacity but higher mutational robustness, that is to say, a greater insensitivity to the deleterious effect of mutations. This implies that mutational robustness may be optimised by natural selection and, consequently, that quasispecies with higher mutational robustness could be selected at high mutation rates. Apart from furthering our understanding of the natural history of populations of self-replicating species, this effect raises the possibility of a resistance mechanism to therapies based on "lethal mutagenesis", which would be distinct from and somehow complementary to drug resistance. However, this point is still a matter of contention [43–45].
Both transition to error catastrophe and survival of the flattest are related to the behaviour of quasispecies at high mutation rates [18, 24, 25, 46]. Some authors have explicitly related both concepts, pointing out that entry into error catastrophe "requires that some genotypes or phenotypes are more sensitive to mutation than others" , that it is "an evolutionary phenomenon in which (...) the population evolves to genotypes that are low in fitness but robust to the effects of mutations" , or that is a theory that shows that "the evolutionary potential of a phenotype depends on both its fitness relative to alternative phenotypes and its robustness to mutations" . However, to the best of our knowledge, neither an explicit study of this relationship, nor an explicit reinterpretation of its consequences, has as yet been taken.
This brief review of recent theoretical papers on this subject shows that the interpretation and meaning of error catastrophe, and its consequences, is a contentious subject. In this paper we will show that what has been called an entry into error catastrophe should actually be considered to be as a specific case of survival of the flattest. Accordingly, the transition into error catastrophe is caused by natural selection when the selection pressure is an increment in mutation rate. For this reason, an error threshold can only take place in fitness landscapes with flat enough regions. Thus, entry into error catastrophe is caused neither by the dominance of mutation over selection nor because natural selection ceases to operate.
In section 2.1 we show throughout a very simplified model that the entry into error catastrophe is formally equivalent to the survival of the flattest. In section 2.2 we study how changes the population distribution beyond the error threshold when the mutational robustness of some Hamming classes is changed by modifying its quality factor. In section 2.3 the possible relationship between the entry into error catastrophe and the survival of the flattest is considered. Finally, these and other previous results are interpreted in the discussion by considering the transition to error catastrophe as a specific case of intra-quasispecies survival of the flattest.
2.1. A simple model to show that entry into error catastrophe is caused by survival of the flattest
where is the mean population productivity.
The system of differential equations (1) can be expressed in terms of relative fitness and selection coefficients referred to the master phenotype, as usually done in population genetics. A coefficient of selection, s, indicates the selective advantage or disadvantage of a given phenotype with respect to another of reference. So, an s value greater than zero implies a selective advantage for a given phenotype with respect to that of the reference, a value less than zero implies a disadvantage, and a value of zero implies neutrality.
which is equivalent to dividing the equations of the system (1) by the effective fitness of the master phenotype, i.e. the product of its replicative ability and its mutational robustness: A m q ν .
Thus, the selection coefficient is the quotient between the effective fitness of the mutant phenotype over that of the master phenotypes, minus one. Actually, the effective fitness of the mutant phenotype should be considered as , where the tilde denotes that the quality factor is considered at a phenotypic level, with , as back mutation is neglected.
2.2. Beyond the error threshold the population evolves to more mutational robust regions of the sequence space
The matrix W is given by W = QA - D , where matrices A and D are diagonal matrices whose elements are the amplification and degradation factors of the different Hamming classes, respectively. Matrix Q is the mutational matrix whose elements Qij determine the probability of obtaining any class i sequence from a class j sequence. The largest eigenvalue, λ, is the mean fitness of the population at the steady state, and its associated right eigenvector v is the population distribution in this state. As this problem is not analytically tractable, even for small sequence lengths, it was solved numerically using MATLAB®.
Beyond the error threshold, the population is delocalized over the whole sequence space. This delocalization is due to the fact that all the sequences of the mutant phenotype have the same amplification factor and the same quality factor, so natural selection ceases to operate. However, this is an extreme restriction, and no general conclusion should be derived from it. When either the amplification factor or the robustness of some mutant Hamming classes is changed, the population distribution beyond the error threshold departs from the uniform distribution. However, the error threshold is not necessarily modified by these changes.
We will now study how the population distribution beyond the error threshold changes when the robustness of some Hamming classes is modified. As the number of sequences is an intrinsic property of the classification, and depends on the Hamming distance with respect to the master sequence, in order to change the robustness of each Hamming class, we assume that the quality factor per digit of each species, qi, depends on its sequence, as previously done in .
2.3. Two competing quasispecies with different robustness
In previous sections we showed that the entry intro error catastrophe implies the selection of a flatter phenotype when the population confronts a higher mutational pressure. In this section we will consider the possible relationship between terror catastrophe within a 12 quasispecies and the survival of the flattest between quasispecies. To study this relationship a simplified model of two quasispecies, A and B, similar to those presented in [49–51] to study the survival of the flattest, has been used. In this model, quasispecies A and B compete with each other but there is no mutational flow between them, i.e. a component of a quasispecies is very unlikely to produce a component of the other quasispecies by mutation. Each quasispecies is made up of a master phenotype (m), and a mutant phenotype (k). The four phenotypes have amplification factors AAm, AAk, ABm and ABk, in each of which the first subscript specifies the quasispecies and the second the phenotype. The degradation factor, D, is the same for all of them and equal to zero.
The per digit genotypic quality factor q is also assumed to be also the same for all of them, but each phenotype must have a different mutational robustness. The phenotypic quality factor of each phenotype i is determined by the expression , where ai is an arbitrary parameter inversely related to robustness [42, 52]. Consequently, the master and mutant phenotypes of quasispecies A have a robustness inversely related to aAm and aAk, respectively and, in a similar way, those of quasispecies B have a robustness inversely related to aBm and aBk.
3.1. The entry into error catastrophe is an specific case of survival of the flattest
The aim of this paper has been to show that the transition into error catastrophe is the result of natural selection acting on the differences in the trade-off between replicative ability and mutational robustness, also called effective fitness. Thus, entry into error catastrophe is a specific case of survival of the flattest. To show this more clearly, a minimal quasispecies model that displays an error threshold was reformulated in section 2.1 in terms of relative effective fitness and selection coefficients. Using this new formulation, we showed that the error threshold is the value of q for which the selection coefficient, s, equals zero (Figure 1). Therefore the selection coefficient is positive beyond the error threshold, which means that beyond that value of q the mutant phenotype has a selective advantage over the master phenotype.
Equation 6 shows that the selection coefficient is the quotient between the effective fitness of the mutant and the master phenotype. As back mutation has been neglected, the probability of going from any mutant sequence to any other sequence of the mutant phenotype is unity, so it does not appear in the numerator.
In section 2.2 we showed with the classical extended model in which sequences are grouped in Hamming classes that the transitions characteristic of an error threshold can be observed when the selection coefficient obtained in the previous section equals zero (Figure 2). This implies that the use of this formal representation for the selection coefficient is still valid for the extended model, when back mutation is taken into account.
Beyond the error threshold, a quasispecies is delocalized over the sequence space in such a way as to produce a uniform population distribution. This is hardly surprising, as the mutant phenotype comprises all but one of the sequences of the sequence space, so it is effectively flat both with respect to the amplification factor and the quality factor. However, the mutational robustness of some mutant sequences can be modified by considering a quality factor that depends on Hamming distance. In this case, a uniform distribution is not obtained beyond the error threshold, and the population evolves in such a way that the regions of the mutant phenotype with a higher quality factor are more populated. That is to say, a population evolves to regions with a greater mutational robustness, as has been shown previously for other similar situations [53, 54]. A remarkable result is the decoupling between the value of the error threshold and the delocalization in error catastrophe. Although the population distribution changes when the mutational robustness of the Hamming classes is modified, the error threshold is not altered (Figure 4b). This can be explained by the fact that changing the robustness distribution of the Hamming classes does not modify the effective fitness of the mutant phenotype as a whole, so the error threshold is not affected. As we will discuss below, this is not the case when lethality is introduced.
3.2. Revisiting previous results
Error catastrophe and error thresholds have been extensively studied since they were first postulated in 1971 . In this section, some previous results obtained in different papers on the error threshold and error catastrophe will be reinterpreted in terms of natural selection and survival of the flattest.
3.2.1. Revisiting the effect of lethality
The same paper studied the changes in the error catastrophe distribution caused by the introduction of lethality. Briefly, the introduction of lethality decreases the average Hamming distance at the error threshold, which subsequently increases linearly with the mutation rate. As a consequence, the uniform distribution is only obtained at q = 0.5. This result is analogous to the result obtained in section × of this paper, in which the effect of the quality factor that depends on the Hamming distance was considered. The introduction of the lethality scheme used in  implies that Hamming classes that are further away from the master sequence have more lethal sequences, and therefore a lesser effective fitness: as both robustness and replicative ability decrease. This means that Hamming classes with fewer lethal sequences have a selective advantage, so that the population distribution is displaced to lower average Hamming distances, once again reproducing previous results [53, 54]. However, unlike the case when the quality factor depends on the Hamming distance, the introduction of lethality decreases the error threshold, as the fitness of the mutant phenotype as a whole is changed.
3.2.2. Revisiting the effect of neutrality and canalization
Since error threshold is a selective transition to more robust phenotypes when the mutation rate is increased, we will study in this section how error threshold is affected by considering neutrality and/or canalization in the master and mutant phenotypes. For the sake of clarity, canalization is considered to be the existence of different genotypes that produce the same phenotype , and neutrality to be the existence of different phenotypes with the same fitness .
As far as we aware, the first account of neutrality in the quasispecies theory was given by Eigen and co-workers , who introduced a gene segment with no influence in the replicative ability of the self-replicating species. As mutations in this segment do not affect fitness, the segment can vary freely and it is not taken into account in the maximum length allowed by the error threshold. In our opinion, this points to the existing relation between natural selection and the error threshold.
The introduction of canalization in the master phenotype through neutral networks led the introduction of the distinction between genotypic and phenotypic error thresholds . When the master phenotype is composed of more than just one sequence, the genetic information of the master phenotype is lost for any value of mutation rate different from zero, beyond which the population wanders through the master neutral network [28, 59]. It is therefore necessary to define a phenotypic error threshold for which the master phenotype is lost, and the population begins a random walk through the phenotype space. However, this stricter definition does not change the considerations made in the previous section, namely, that a phenotypic error catastrophe has no influence at a higher selective-evolutionary level, as almost all those phenotypes have the same or almost the same fitness. In fact, the very idea of the genotypic error threshold reinforces this conclusion, as the loss of the genetic information of the master phenotypes has not critical consequences for their existence or evolutionary capacities.
Equation 16 shows that the effect of increasing the neutral network of the master phenotype is to increase its mutational robustness, and therefore the effective fitness of the master phenotype, at higher mutation rates. As a consequence the mutant phenotype outcompetes the master phenotype at a higher mutation rate than in the case with no neutrality, so error threshold decreases [31, 59, 60].
3.2.3. Revisiting multiple error thresholds
Several papers have described the existence of multiple error thresholds, also known as error cascades, in complex fitness landscapes [18, 61]. In point of fact, this phenomenon had probably been observed previously, although it was not identified as such, and was instead identified as either intermediate regions between high replication accuracy regions and error catastrophe  or as partial "delocalization transitions" to flatter regions of the sequence space .
From the point of view of the survival of the flattest and natural selection, the appearance of multiple error thresholds is just a consequence of the existence of multiple phenotypes with different trade-offs between replicative ability and mutational robustness . As the effective fitness of these different phenotypes can change differently with mutation rate, natural selection can induce several consecutive transitions between them. Actually, error cascades have been previously related with survival of the flattest [20, 61].
Finally, the existence of multiple error thresholds shows that delocalization over a given region of the sequence space is essentially dependent on the phenotype's degeneration. Therefore, delocalization can take place over the whole sequence space, or just over some limited regions corresponding to the neutral network of a given phenotype, as in the case of the so-called genotypic error catastrophe commented above.
3.3. Error catastrophe, survival of the flattest, and levels of selection
The main purpose of this paper has been to show that the so-called entry into error catastrophe is a specific form of survival of the flattest, that is to say: it is the consequence of natural selection acting in systems of self-replicating species at high mutation rates. As they are essentially equivalent phenomena, and thus have the same cause and features, we propose differentiating between them by considering that they refer to the two different levels of selection that appears in quasispecies models. The first selection level is that of individual self-replicating species, grouped in phenotypes, which, through selection and mutation, determine the population distribution of the quasispecies. The second is a higher level in which some quasispecies can compete with others, through their emerging biological fitness derived from the interaction of their components .
Therefore, we propose to use the term "survival of the flattest" to refer to situations in which two quasispecies compete, and there is no mutational coupling between them. This is either because the number of mutations between them implies that the possibility of obtaining one from the other in a reasonable time is negligible [39, 62], or because of structural or functional assumptions that are translated into "replicative isolation" in the model [40–42, 49]. On the other hand, entry into "error catastrophe" can be regarded as a selective transition resulting from the competition between two or more mutationally coupled phenotypes within the same quasispecies, each one with a different amplification factor and/or mutational robustness, in which a more robust phenotype displaces partially or totally another phenotype with more replicative ability but less robustness.
The study of the system of differential equations presented in section 2.4 shows that the changes produced by natural selection in the internal population structure of a quasispecies, i.e. entry into error catastrophe, can modify the result of natural selection on the immediately higher level of selection, that of competition between quasispecies. In fact, it can be said that entry into error catastrophe induces a survival of the flattest.
3.4. ¿Survival of the flattest vs. survival of the fittest?
At a metaphorical level, survival of the flattest and survival of the fittest are often used to denote two completely different concepts [40, 41, 49, 63]. We think that this is confusing, and can even give rise to serious misunderstandings. We do not believe that the term "survival of the flattest" should be discarded, but neither should it be used to signify the opposite of "survival of the fittest". We think it is more accurate to differentiate between "selection for replicative speed" and "selection for mutational robustness"  when comparing both processes. This differentiation stresses that natural selection always promotes selection of the fittest. The point is that when the mutation rate is small, the fittest entity is the entity with the greatest replicative ability, but when the mutation rate is high enough, the fittest entity could be one with the greatest mutational robustness. In the latter case, the fittest entity is the fittest because it is the flattest. Thus, mutational robustness should be regarded as another component of biological fitness, at the same level as the replicative ability, expressed in the amplification factor Ai. Accordingly, mutational robustness may play a key role in determining the effect of natural selection at high mutation rate conditions . This point can be seen more clearly in section 2.1. when, instead of the intrinsic replicative ability, Ai, the product of the replicative ability and mutational robustness (the effective fitness) are compared, and the transition to error catastrophe arises naturally as a result of natural selection acting on the difference of fitness.
3.5. Information crisis? What information crisis?
As briefly commented in the introduction, error catastrophe has been related either to an information crisis [1, 6, 16, 26, 29], or to a breakdown of evolutionary adaptation [2, 22, 28]. In light of the results and interpretations presented in this paper, we think that is quite inaccurate to speak of any kind of crisis beyond the error threshold.
The so-called "information crisis" beyond the error threshold sometimes is associated to the loss of the master sequence [1, 26, 29], whereas in other cases it is associated to the delocalization over the whole sequence space . When the quasispecies crosses the error threshold, information is said to melt, suffering a phase transition [6, 16, 26, 64]. As Hermisson et al  have shown, both phenomena coincide for the single-peak landscape, but not for other landscapes. In any case, neither the delocalization beyond the error threshold nor the loss of the master sequence can be associated to a critical loss of information.
In the first place, as Eigen has shown, natural selection implies displacements in the 'information space', which can take the form of phase transitions [62, 64]. A phase transition is a critical phenomenon, but it is critical because of some physico-mathematical characteristics of the transition and not necessarily because of the effects it produces. The error threshold shows a phase transition for collective order parameters such the consensus sequence [11, 16], but not, for example, for the average fitness of the population, whose value does not greatly vary at the error threshold, although its dependence on mutation rate does . Secondly, although application of the concept of "information" to biology can be very useful, it does have its limitations. Delocalization beyond the error threshold causes an "information melting", but only at the symbolic level (i.e. the genotype), the one to which the mathematical theory of communication can be applied . However, as all the mutant sequences can be grouped in the same phenotype, and with the same fitness, the symbolic delocalization has no critical consequences at the evolutionary, i.e. biological, level. As a matter of fact, taking the information metaphor further, it could be said that this "symbolic melting" has no consequences at higher informational levels: neither semantic (the meaning, i.e.: the phenotype) nor syntactic (relation between meanings or phenotypes).
Finally, the loss of the master sequence obviously implies the loss of the information contained in it. However, any natural selection process implies changes in the informational characteristics of the population. Obviously, if the result of natural selection is the disappearance of a given phenotype of the population, that implies an information loss, but, at the same time, it implies an increment in biological fitness. Again, no critical consequence can be expected.
A similar reasoning can be applied to the idea of a "breakdown of evolutionary adaptation" [22, 28]. As a result of natural selection, entry into error catastrophe does not imply the "breakdown of evolutionary adaptation" but, on the contrary, the adaptation to higher mutation rates by increasing the mutational robustness of the population. If there is no adaptation beyond that point it is because the population has reached the most adapted state allowed by the model.
In the specific case of the single-peak fitness landscape, entry into error catastrophe would be the most extreme case of selection for mutational robustness within a quasispecies. On the one hand, there is a master phenotype which, as it is made up of a single genotype, the master sequence, has null mutational robustness. On the other hand, there is a mutant phenotype which, as it is made up of the rest of the 2ν -1 possible sequences, has a virtually infinite mutational robustness (and thus, aiK = 0 in section 2.3). When confronting a mutational pressure, an infinite mutational robustness is probably one of the best adaptations that a system can reach. Although the mean fitness is insensitive to further increases in mutation rate beyond the error threshold, the condition of which Hermisson et al.  have termed a "degradation threshold", this insensitivity results from the fact that the populations has reached an infinite robust phenotype, i.e. the fittest phenotype, which means that calling this phenomenon a "complete mutational degradation" is, at the least, misleading.
3.6. Error threshold and RNA viruses
The error threshold concept acquired great importance, among other reasons, when the quasispecies concept was applied to RNA viruses , as it established a possible new antiviral mechanism which was radically different to those which had been studied before . In fact, the error threshold concept inspired a series of experimental papers which have demonstrated the possibility of extinguishing some RNA viruses using mutagens, in a process called lethal mutagenesis. However, several authors [18, 19] have recently called into question the idea that lethal mutagenesis is a result of entry into error catastrophe, correctly pointing out that error catastrophe is a genetic-evolutionary process while extinction is a demographic process . In this paper we have shown that error catastrophe can be regarded as a particular case of natural selection for mutational robustness within a quasispecies. Consequently, that conclusion is reinforced, i.e. entry into error catastrophe cannot explain viral extinction due to increased mutagenesis. The quasispecies model, however useful it has been for the study of RNA viruses [3, 67, 68], is still only an initial approach to complex intra- and extracellular viral dynamics.
The main conclusion of this paper is that the entry into error catastrophe is a specific case of survival of the flattest acting on phenotypes which differ in the trade-off between replicative ability and mutational robustness. In fact, the entry into error catastrophe takes place when the mutant phenotype acquires a selective advantage over the master phenotype. Moreover, beyond the error thresholds, changing the quality factor of some sequences modifies the population distribution at the error catastrophe, displacing it towards the flatter regions of the mutant phenotype. However, the value of the error threshold is not altered by these changes in the mutant phenotype as it depends on its effective fitness as a whole. Both neutrality and some lethality schemes increase the effective fitness of the master phenotype with respect to the mutant phenotype, so error threshold decreases. Taking this into account the notion of crisis information beyond error threshold does not make sense.
As both entry into error catastrophe and survival of the flattest are caused by natural selection when mutation rate is increased, we propose differentiating between them by the level of selection at which natural selection act. Thus, we propose to use the term "survival of the flattest" to refer to situations in which two quasispecies compete, and there is no mutational coupling between them; and the term "entry into error catastrophe" as the displacement of a phenotype with high replicative ability but less robustness by another flatter phenotype, when they are mutationally coupled and within the same quasispecies.
This paper has been supported in part by Grants no. BFU2009-12895-C02-02 from MEC (Spain). Hector Tejero is supported by AP2006-01044, from MEC (Spain); and Arturo Marín is supported by Consejería de Educación de la Comunidad de Madrid (Spain) and Fondo Social Europeo (FSE). We would like to thank Esteban Domingo for useful discussion and critical reading of this manuscript, and Athel Cornish-Bowden for the final improving of the written English.
- Eigen M: Self-organization of matter and the evolution of biological macromolecules. Naturwissenschaften. 1971, 58 (10): 465-523. 10.1007/BF00623322.View ArticlePubMedGoogle Scholar
- Schuster P, Stadler PF: Early replicons: Origin and evolution. Origin and Evolution of Viruses. Edited by: Domingo E, Parrish C, Holland JJ. 2008, Oxford: Elsevier, 1-42. full_text.View ArticleGoogle Scholar
- Domingo E, Martin V, Perales C, Grande-Perez A, Garcia-Arriaza J, Arias A: Viruses as quasispecies: biological implications. Curr Top Microbiol Immunol. 2006, 299: 51-82. full_text.PubMedGoogle Scholar
- Solé RV, Deisboeck TS: An error catastrophe in cancer?. J Theor Biol. 2004, 228 (1): 47-54.View ArticlePubMedGoogle Scholar
- Kamp C, Wilke CO, Adami C, Bornholdt S: Viral evolution under the pressure of an adaptive immune system: optimal mutation rates for viral escape. Complexity. 2002, 8 (2): 28-33. 10.1002/cplx.10067.View ArticleGoogle Scholar
- Biebricher C, Eigen M: The error threshold. Virus Res. 2005, 107 (2): 117-127. 10.1016/j.virusres.2004.11.002.View ArticlePubMedGoogle Scholar
- McCaskill JS: A localization threshold for macromolecular quasispecies from continuously distributed replication rates. J Chem Phys. 1984, 80 (10): 5194-5202. 10.1063/1.446590.View ArticleGoogle Scholar
- Swetina J, Schuster P: Self-replication with errors: A model for polynucleotide replication. Biophys Chem. 1982, 16 (4): 329-345. 10.1016/0301-4622(82)87037-3.View ArticlePubMedGoogle Scholar
- Swetina J: First and second moments and the mean Hamming distance in a stochastic replication-mutation model for biological macromolecules. J Math Biol. 1989, 27 (4): 463-483. 10.1007/BF00290640.View ArticlePubMedGoogle Scholar
- Wiehe T: Model dependency of error thresholds: the role of fitness functions and contrasts between the finite and infinite sites models. Genet Res. 1997, 69 (02): 127-136. 10.1017/S0016672397002619.View ArticleGoogle Scholar
- Bonhoeffer S, Stadler PF: Error Thresholds on Correlated Fitness Landscapes. J Theor Biol. 1993, 164 (3): 359-372. 10.1006/jtbi.1993.1160.View ArticleGoogle Scholar
- Tarazona P: Error thresholds for molecular quasispecies as phase transitions: From simple landscapes to spin-glass models. Phys Rev A. 1992, 45 (8): 6038-10.1103/PhysRevA.45.6038.View ArticlePubMedGoogle Scholar
- Hermisson J, Redner O, Wagner H, Baake E: Mutation-selection balance: ancestry, load, and maximum principle. Theor Popul Biol. 2002, 62 (1): 9-46. 10.1006/tpbi.2002.1582.View ArticlePubMedGoogle Scholar
- Baake E, Gabriel W: Biological evolution through mutation, selection, and drift: An introductory review. Ann Rev Comp Phys VII. Edited by: Stauffer D. 2000, Singapore: World Scientific, 203-264. full_text.Google Scholar
- Schuster P: Mathematical modeling of evolution. Solved and open problems. Theory Biosci. 2010Google Scholar
- Eigen M: Error catastrophe and antiviral strategy. Proc Natl Acad Sci USA. 2002, 99 (21): 13374-13376. 10.1073/pnas.212514799.View ArticlePubMedPubMed CentralGoogle Scholar
- Eigen M, Schuster P: Stages of emerging life --Five principles of early organization. J Mol Evol. 1982, 19 (1): 47-61. 10.1007/BF02100223.View ArticlePubMedGoogle Scholar
- Bull JJ, Meyers LA, Lachmann M: Quasispecies made simple. PLoS Comput Biol. 2005, 1 (6): e61-10.1371/journal.pcbi.0010061.View ArticlePubMedPubMed CentralGoogle Scholar
- Summers J, Litwin S: Examining The Theory of Error Catastrophe. J Virol. 2006, 80 (1): 20-26. 10.1128/JVI.80.1.20-26.2006.View ArticlePubMedPubMed CentralGoogle Scholar
- Wilke CO: Quasispecies theory in the context of population genetics. BMC Evol Biol. 2005, 5 (44):Google Scholar
- Campos PRA, Adami C, Wilke CO: Optimal adaptive performance and delocalization in NK fitness landscapes. Physica A. 2002, 304 (34): 495-506. 10.1016/S0378-4371(01)00572-6.View ArticleGoogle Scholar
- Takeuchi N, Hogeweg P: Error-threshold exists in fitness landscapes with lethal mutants. BMC Evol Biol. 2007, 7 (15):Google Scholar
- Nowak MA, Schuster P: Error thresholds of replication in finite populations mutation frequencies and the onset of muller's ratchet. J Theor Biol. 1989, 137 (4): 375-395. 10.1016/S0022-5193(89)80036-0.View ArticlePubMedGoogle Scholar
- Bull JJ, Sanjuan R, Wilke CO: Theory of lethal mutagenesis for viruses. J Virol. 2007, 81 (6): 2930-2939. 10.1128/JVI.01624-06.View ArticlePubMedPubMed CentralGoogle Scholar
- Cowperthwaite MC, Meyers LA: How Mutational Networks Shape Evolution: Lessons from RNA Models. Annu Rev Ecol Evol Syst. 2007, 38: 203-230. 10.1146/annurev.ecolsys.38.091206.095507.View ArticleGoogle Scholar
- Biebricher C, Eigen M: What is a quasispecies?. Curr Top Microbiol Immunol. 2006, 299: 131-Google Scholar
- Elena SF, Sanjuan R: Virus Evolution: Insights from an Experimental Approach. Annu Rev Ecol Evol Syst. 2007, 38: 27-52. 10.1146/annurev.ecolsys.38.091206.095637.View ArticleGoogle Scholar
- Huynen M, Stadler P, Fontana W: Smoothness within ruggedness: the role of neutrality in adaptation. Proc Natl Acad Sci USA. 1996, 93 (1): 397-401. 10.1073/pnas.93.1.397.View ArticlePubMedPubMed CentralGoogle Scholar
- Eigen M, Schuster P: The hypercycle: a principle of natural self-organization. 1979, Berlin: SpringerView ArticleGoogle Scholar
- Schuster P: The Origin and Processing of Biological Information. Frontiers of Life, The origins of life. Edited by: Baltimore D, Dulbecco R, Jacob F, Levi-Montalcini L. 2001, Academic Press, I:Google Scholar
- Takeuchi N, Poorthuis P, Hogeweg P: Phenotypic error threshold; additivity and epistasis in RNA evolution. BMC Evol Biol. 2005, 5 (1): 9-10.1186/1471-2148-5-9.View ArticlePubMedPubMed CentralGoogle Scholar
- Obermayer B, Frey E: Error thresholds for self- and cross-specific enzymatic replication. J Theor Biol. 2010Google Scholar
- Szathmáry E, Smith J: The major evolutionary transitions. Nature. 1995, 374 (6519): 227-232.View ArticlePubMedGoogle Scholar
- Eigen M: The fifth Paul Ehrlich lecture. Virus strains as models of molecular evolution. Med Res Rev. 1993, 13 (4): 385-398. 10.1002/med.2610130402.View ArticlePubMedGoogle Scholar
- Cameron CE, Castro C: The mechanism of action of ribavirin: lethal mutagenesis of RNA virus genomes mediated by the viral RNA-dependent RNA polymerase. Curr Opin Infect Dis. 2001, 14 (6): 757-764.View ArticlePubMedGoogle Scholar
- Holland JJ, Domingo E, de la Torre JC, Steinhauer DA: Mutation frequencies at defined single codon sites in vesicular stomatitis virus and poliovirus can be increased only slightly by chemical mutagenesis. J Virol. 1990, 64 (8): 3960-3962.PubMedPubMed CentralGoogle Scholar
- Loeb LA, Essigmann JM, Kazazi F, Zhang J, Rose KD, Mullins JI: Lethal mutagenesis of HIV with mutagenic nucleoside analogs. Proc Natl Acad Sci USA. 1999, 96 (4): 1492-1497. 10.1073/pnas.96.4.1492.View ArticlePubMedPubMed CentralGoogle Scholar
- Sierra S, Davila M, Lowenstein PR, Domingo E: Response of foot-and-mouth disease virus to increased mutagenesis: influence of viral load and fitness in loss of infectivity. J Virol. 2000, 74 (18): 8316-8323. 10.1128/JVI.74.18.8316-8323.2000.View ArticlePubMedPubMed CentralGoogle Scholar
- Schuster P, Swetina J: Stationary mutant distributions and evolutionary optimization. Bull Math Biol. 1988, 50 (6): 635-660.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 (6844): 331-333. 10.1038/35085569.View ArticlePubMedGoogle Scholar
- Wilke CO: Selection for fitness versus selection for robustness in RNA secondary structure folding. Evolution. 2001, 55 (12): 2412-2420.View ArticlePubMedGoogle Scholar
- Wilke CO, Adami C: Evolution of mutational robustness. Mutat Res. 2003, 522 (1-2): 3-11. 10.1016/S0027-5107(02)00307-X.View ArticlePubMedGoogle Scholar
- Sanjuán R, Cuevas JM, Furió V, Holmes EC, Moya A: Selection for Robustness in Mutagenized RNA Viruses. PLoS Genetics. 2007, 3 (6): e93-View ArticlePubMedPubMed CentralGoogle Scholar
- Martín V, Grande-Pérez A, Domingo E: No evidence of selection for mutational robustness during lethal mutagenesis of lymphocytic choriomeningitis virus. Virology. 2008, 378 (1): 185-192.View ArticlePubMedGoogle Scholar
- O'Dea EB, Keller TE, Wilke CO: Does mutational robustness inhibit extinction by lethal mutagenesis in viral populations?. PLoS Comput Biol. 2010, 6 (6): e1000811-View ArticlePubMedPubMed CentralGoogle Scholar
- Cowperthwaite MC, Economo EP, Harcombe WR, Miller EL, Meyers LA: The ascent of the abundant: how mutational networks constrain evolution. PLoS Comput Biol. 2008, 4 (7): e1000110-10.1371/journal.pcbi.1000110.View ArticlePubMedPubMed CentralGoogle Scholar
- Thompson CJ, McBride JL: On Eigen's theory of the self-organization of matter and the evolution of biological macromolecules. Math Biosci. 1974, 21 (12): 127-142.View ArticleGoogle Scholar
- Sasaki A, Nowak MA: Mutation landscapes. J Theor Biol. 2003, 224 (2): 241-247. 10.1016/S0022-5193(03)00161-9.View ArticlePubMedGoogle Scholar
- Sardanyés J, Elena SF, Solé RV: Simple quasispecies models for the survival-of-the-flattest effect: The role of space. J Theor Biol. 2008, 250 (3): 560-568.View ArticlePubMedGoogle Scholar
- O'Fallon BD, Adler FR, Proulx SR: Quasi-species evolution in subdivided populations favours maximally deleterious mutations. Proc R Soc Lond B Biol Sci. 2007, 274 (1629): 3159-3164.View ArticleGoogle Scholar
- Aguirre J, Manrubia SC: Out-of-equilibrium competitive dynamics of quasispecies. Europhys Lett. 2007, 77 (3): 38001-10.1209/0295-5075/77/38001.View ArticleGoogle Scholar
- Bloom JD, Arnold FH, Wilke CO: Breaking proteins with mutations: threads and thresholds in evolution. Mol Syst Biol. 2007, 3: 76-10.1038/msb4100119.View ArticlePubMedPubMed CentralGoogle 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Forster R, Adami C, Wilke CO: Selection for mutational robustness in finite populations. J Theor Biol. 2006, 243 (2): 181-190. 10.1016/j.jtbi.2006.06.020.View ArticlePubMedGoogle Scholar
- Tejero H, Marín A, Montero F: Effect of lethality on the extinction and on the error threshold of quasispecies. J Theor Biol. 2010, 262 (4): 733-741. 10.1016/j.jtbi.2009.10.011.View ArticlePubMedGoogle Scholar
- Waddington CH: Canalization of development and the inheritance of acquired characters. Nature. 1942, 150: 563-565. 10.1038/150563a0.View ArticleGoogle Scholar
- Krakauer D, Plotkin J: Principles and Parameters of Molecular Robustness. Robust Design: a repertoire of biological, ecological and engineering case studies. Edited by: Jen E. 2005, New York: Oxford University PressGoogle Scholar
- Eigen M, McCaskill JS, Schuster P: Molecular quasi-species. J Phys Chem. 1988, 92 (24): 6881-6891. 10.1021/j100335a010.View ArticleGoogle Scholar
- Schuster P, Fontana W: Chance and necessity in evolution: lessons from RNA. Physica D. 1999, 133 (14): 427-452. 10.1016/S0167-2789(99)00076-7.View ArticleGoogle Scholar
- Reidys C, Forst CV, Schuster P: Replication and mutation on neutral networks. Bull Math Biol. 2001, 63 (1): 57-94. 10.1006/bulm.2000.0206.View ArticlePubMedGoogle Scholar
- Tannenbaum E, Shakhnovich EI: Solution of the quasispecies model for an arbitrary gene network. Phys Rev E Stat Nonlin Soft Matter Phys. 2004, 70 (2 Pt 1): 021903-10.1103/PhysRevE.70.021903.View ArticlePubMedGoogle Scholar
- Wolff A, Krug J: Robustness and epistasis in mutation-selection models. Phys Biol. 2009, 6 (3): 036007-10.1088/1478-3975/6/3/036007.View ArticlePubMedGoogle Scholar
- Codoñer FM, Darós JA, Solé RV, Elena SF: The fittest versus the flattest: experimental confirmation of the quasispecies effect with subviral pathogens. PLoS Pathog. 2006, 2 (12): e136-View ArticlePubMedPubMed CentralGoogle Scholar
- Eigen M: Natural selection: a phase transition?. Biophys Chem. 2000, 85 (23): 101-123. 10.1016/S0301-4622(00)00122-8.View ArticlePubMedGoogle Scholar
- Shannon CE, Weaver W: The Mathematical Theory of Communication. 1949, Urbana: University of Illinois PressGoogle Scholar
- Domingo E, Sabo D, Taniguchi T, Weissmann C: Nucleotide sequence heterogeneity of an RNA phage population. Cell. 1978, 13 (4): 735-744. 10.1016/0092-8674(78)90223-4.View ArticlePubMedGoogle Scholar
- Domingo E, Biebricher CK, Eigen M, Holland JJ: Quasispecies and RNA virus Evolution: Principles and Consequences. 2001, Austin, TX: Landes BioscienceGoogle Scholar
- Domingo E: Quasispecies Theory in Virology. J Virol. 2002, 76 (1): 463-465. 10.1128/JVI.76.1.463-465.2002.View ArticlePubMed CentralGoogle 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.