Quasispecies-like behavior observed in catalytic RNA populations evolving in a test tube
© Arenas and Lehman; licensee BioMed Central Ltd. 2010
Received: 24 August 2009
Accepted: 23 March 2010
Published: 23 March 2010
During the RNA World, molecular populations were probably very small and highly susceptible to the force of strong random drift. In conjunction with Muller's Ratchet, this would have imposed difficulties for the preservation of the genetic information and the survival of the populations. Mechanisms that allowed these nascent populations to overcome this problem must have been advantageous.
Using continuous in vitro evolution experimentation with an increased mutation rate imposed by MnCl2, it was found that clonal 100-molecule populations of ribozymes clearly exhibit certain characteristics of a quasispecies. This is the first time this has been seen with a catalytic RNA. Extensive genotypic sampling from two replicate lineages was gathered and phylogenetic networks were constructed to elucidate the structure of the evolving RNA populations. A common distribution was found in which a mutant sequence was present at high frequency, surrounded by a cloud of mutant with lower frequencies. This is a typical distribution of quasispecies. Most of the mutants in these clouds were connected by short Hamming distance values, indicating their close relatedness.
The quasispecies nature of mutant RNA clouds facilitates the recovery of genotypes under pressure of being removed from the population by random drift. The empirical populations therefore evolved a genotypic resiliency despite a high mutation rate by adopting the characteristics of quasispecies, implying that primordial RNA pools could have used this strategy to avoid extinction.
During the origins and early evolution of life on the Earth, the contemporary notion of well-delineated species was not yet realized. Instead, the adaptive fate of populations of naked molecules would have been more accurately described by the quasispecies model of mutant distributions, which was introduced in the 1970's by Eigen [1–3]. Under an RNA World scenario, the molecules in question were autocatalytically replicating polymers of nucleotides, such as RNA or a chemical equivalent [4, 5].
A quasispecies is basically a steady-state dynamic of mutant molecules distributed around a parental genotype, the so-called master sequence, which occupies a central position in the genotypic network space. This dynamic occurs after a sufficient amount of time at high mutational rates, such that the progeny of an individual genotype (the mutant cloud) can be rapidly produced. Importantly, the target of selection is the genotypic distribution as a whole, not single genotypes. Eigen realized that simple genetic entities would gain an evolutionary advantage by a primitive form of group selection in which replication rate deficiencies could be offset by the heightened production of new mutant types . Different genotypes can form mutant clouds of various sizes that can compete for survivorship during the evolution of the population. This can generate a fluctuating equilibrium dynamic as clouds of mutants are replaced by other ones at the interplay of selection and random drift .
Although the quasispecies concept was originally intended as a description of molecular replicators cooperating and competing for survival prior to their historical encapsulation in membranes that inexorably linked genotype with phenotype [5, 7], it soon was argued that many viral populations evolved with the same dynamic. The first empirical demonstration of quasispecies behavior was produced in 1978 with the Qβ phage . Since then many viral populations have been described as evolving as quasispecies collectives under mutation-selection balance [9–15], although this interpretation has been debated in the case of viruses [7, 16]. Moreover, many experiments, particularly by Biebricher and co-workers, have clearly highlighted the power of the quasispecies concept to describe the evolutionary progress of naked RNA molecules in the test tube as they are replicated by error-prone protein polymerases [17–19]. Such experiments confirm some of the key predictions of quasispecies dynamics as outlined by Eigen [5, 7]: a cloud of neutral mutations exists that surround a central genotype, selection operates on this cloud as a whole, in small and error-prone populations many specific mutants occur reproducibly (recur), and the greatest amount of genotypic diversity exists just below an error threshold above which information decays into chaos.
To date however, the use of a catalytic RNA in a system capable of forming quasispecies has not been reported. Currently available data are restricted to genomic or genome-fragment RNAs that can carry genetic information but not perform a catalytic function. The aim of the current work is to extend the quasispecies concept back at least as far as to a system in which the catalytic function of an RNA molecule is integral to its own replication process, even if the actual polymerization of nucleotides is carried out by an exogenous polymerase. Here, the continuous evolution (CE) in vitro evolution  of class I ligase ribozymes is used to create an evolving RNA population that can be tracked over hundreds of generations in a very short time frame. This system has previously been shown to exhibit mutational meltdown  of small populations when the mutation rate is limited by the intrinsic error rate of viral polymerases . Now, with the error rate greatly enhanced with the addition of a chemical mutagen (manganese (II) ion) such that prebiotic conditions are better simulated, it is shown that ribozyme populations do indeed display quasispecies-like dynamics while evolving in a test tube.
To investigate the cause of the extended time to extinction, a preliminary inspection of the populational genetic variability using RFLP was performed. In general, the cDNA in a population at any burst can be amplified via PCR and then genotyped by either RFLP or nucleotide sequence analysis. The fixation, or nearly so, of mutant forms was evidenced in all the lineages at selected bursts. Based on RFLP assays [Additional file 1: Supplemental Figure S2] two lineages (6H and 6L) were selected for a more extensive characterization of the genotypes present in the populations. Genotypes from three or four bursts, respectively, in lineages 6H and 6L (Figure 1B) were cloned and sequenced exhaustively enough to gather a representative sample of the population diversity. In addition, two bursts from the smallest surviving lineage from the previous study  without added MgCl2 (the 600-molecule lineage 3D) were genotyped for comparison.
Summarized characteristics of the quasispecies in CE lineages.
mutagen (MnCl2) added?
master sequence identity
frequency of most common sequence (%)
sample size (sequences)
# of unique genotypes
Shannon Diversity Index
mean Hamming distance
The observed clouds also showed that there existed a fluctuating dynamic in these populations, as the shape of the clouds -- and gross amount of mutant sequences -- changed from one burst to another (Figures 2D and 3E). These results indicate that each lineage can develop a different dynamic over the course of 50 bursts. It is likely that the high mutation rate used, coupled with the relative stability of a putative quasispecies structure, allows the populations to explore multiple viable alternatives of sequence space during the course of their evolutionary history. In lineage 6H, an early time point (burst 5 out of 50) shows that a quasispecies-like cloud is already formed (Figure 2A), in which the master sequence is still the "wildtype" B16-19, with a frequency of 76%. This cloud was constructed from a sample of 102 clones, which contained 16 different genotypes (Table 1).
The numbers of clones in each burst that needed to be genotyped by nucleotide sequence analysis in order to sample the bulk of unique genotypes present in the population was evaluated by constructing rarefaction plots. Because between 41-108 clones per burst were genotyped (Table 1), and yet the harmonic mean population size of these bursts was ~100 individuals and the frequency distributions of genotypes were highly skewed (see above), these rarefaction analyses indicated that this sampling was exhaustive enough to capture all but the rarest of genotypes. The required number of clones to be inspected was similar when comparing the bursts within the lineages, but different when comparing the two lineages [Additional file 1: Supplemental Table S1]. This result is not surprising, because a quasispecies would be a fairly stable equilibrium dynamic, and the environmental conditions are nearly constant during the CE experiments. It should be noted that, although the population size was ostensibly kept constant at 100 molecules throughout each lineage, the cloning procedure did involve PCR amplification. Therefore, the sampling of genotypes from the population is effectively sampling with replacement of the total diversity. For example, 102 clones from burst 5 of lineage 6H were obtained, but this did not mean that all individuals were sampled. Nevertheless, the observed sample diversity in this burst was estimated at 1.49, as measured by Shannon Diversity Index (H = -Σp i (lnp i ) where p i is the frequency of the ith clone). Deeper examination of the bursts in this lineage (Figures 2B and 2C) revealed a change in the master sequence identity, frequency, number of genotypes, and Shannon diversity values. In general, the quasispecies formed in each burst presented different identities from each other, but their characteristics are fairly similar (Table 1). This similarity is perhaps the result of the fact that the sequence space available for exploration by the populations is bounded by a unique starting point; they are all genotypically identical at the beginning of the evolution experiment. Therefore, the area of sequence space that can be explored in 50 serial transfers would be relatively short, and the lineages may be not very different in their diversity values.
Genotypes that are present at a higher frequency in one burst can become a master sequence at a later burst. Conversely, a master sequence that, having once been displaced, was never observed to come back to high frequency in the population. For example, in lineage 6H, the following transition can be observed in the master sequence identity and frequency: burst 5, B16-19 (86%); burst 35, MS1 (91%); burst 50, MS2 (48%). This dynamic of master sequences being displaced by one another resembles that of clonal interference, in which advantageous mutants have to compete for resources and some get displaced [31–33]. A network drawn by combining the sequences of all bursts inspected in lineage 6H (Figure 2D) shows the fluctuation of the equilibrium dynamic of the quasispecies in its evolutionary history.
In lineage 6L, four bursts were evaluated (Figure 3). These specific bursts were chosen because the RFLP analysis detected genotypic transitions near these time points [Additional file 1: Supplemental Figure S2]. The dynamic of the lineage is similar to that of lineage 6H, in that different quasispecies clouds emerge and evolve through time. Two successive bursts (#22 and #23; Figure 1B) reveal that this fluctuation can occur in a relatively short time (Figures 3A and 3B). Interestingly, some of the master sequences that appear during this lineage (Figures 3B, C, and 3D) are the same as were observed in lineage 6H (Figures 2B and 2C). These results suggest that some quasispecies may develop a stronger mutant coupling than others that enables them to recurrently out-compete other quasispecies present during the lineages evolution. Similar to the pattern in lineage 6H, the characteristics of the quasispecies change over the course of the evolutionary history as indicated by the master sequence identity, frequency and number of genotypes, and diversity values calculated (Table 1). These dynamics of staggered dominant genotypes that fluctuate as the population evolves (Figure 3E) may be a reflection of the interplay of Darwinian selection and random genetic drift acting on the quasispecies.
Other than the genotype used to seed these experiments (B16-19 = the "wildtype"), the genotypes that appeared in lineages 6H and 6L were not identical to the dominant genotypes that appeared in lineages evolved without added mutational pressure . Some specific mutations, such as the U62A "insurance" mutation , did appear recurrently in the quasispecies in lineages 6H and 6L, but the composite genotypes seen in this study were not those observed previously. In particular, the "immunity" mutation (a change from CUGAACCUUA to AAUCG at positions 123-132), which conferred resistance to mutational meltdown in small population sizes in the absence of MnCl2, was not seen in the current study, although short insertions, deletions, and rearrangements were common in the last 10-12 nt at the 3' end of the ligase (Figure 1), and specific mutations in this region appeared in more than one quasispecies. Recurrence of mutations is actually a predicted characteristic of quasispecies , thus providing additional support for this interpretation of these data.
"What is a quasispecies?" is the exact title of at least two papers [37, 38]. At times this phenomenon has clearly been difficult to detect and apply to real populations with absolute certainty. However the defining characteristics of a quasispecies, as described by Eigen, unequivocally include a spectrum of closely related genotypes and a population that is struggling to survive under a relatively high error rate. These features seem to apply to primordial collections of RNA near the origins of life, and thus prompted the current study. Here, clonal 100-molecule populations of B16-19 ligase ribozymes were evolved using the continuous in vitro evolution (CE) method  and the relatively error-prone MMLV-RT. In particular, MnCl2 was added to the reaction vessel to increase the error rate of protein enzymes. Populations evolved under these conditions did not show a shortened extinction time, as was observed previously when mutational meltdown conditions prevailed under only a weak mutational pressure of no added MnCl2 . The data suggest that this is a consequence of the advent of quasispecies in these ribozyme populations.
Quasispecies behavior has never before been demonstrated during in vitro evolution experiments with catalytic RNA. Other in vitro experiments with ribozymes have shown either convergence on a phenotype or recurrence of a genotype or motif, but not the type of dynamic of quasispecies that we are documenting here. For example, Yingfu Li and colleagues studied how the composition of a population of RNA-cleaving DNAzymes changed over time in response to selective pressures acting on the phenotype [39, 40]. Similar to the findings reported here, they found a dynamic fluctuation in the structure of the population. Many sequence classes peaked in frequency at different rounds of selection, but in this case, one class appeared to consistently maintain a high frequency. It will be interesting to explore the population structure that these DNAzymes would adopt if the mutation rate of the replication were increased. Perhaps mutational coupling may arise in these molecular populations as well. Another classic study of interest was the evolution of the RNA variant V2 of the Qβ virus performed by Orgel and co-workers . In this case, the mutagen ethidium bromide (EthBr) was added to the reaction vessel during the serial transfers. RNAs resistant to EthBr evolved and adapted to increasing drug concentrations. However in this case, in contrast to the CE/Mn2+ mutagenesis experiments presented here, the mutagen has a direct effect on the RNA structure, and therefore selection favored variants that mutated away the EthBr-binding sites, and a single "winner" emerged. Examples of experiments in which there is a convergent to a relatively more efficient solution for the population have been reviewed [23, 42–45]. Other in vitro evolution studies have been initiated from pools of limited diversity but in which a recurrent solution would appear upon independent trials. This was the evolution dynamic that gave rise to the class I ligase , and the group I ribozymes , but quasispecies-like behavior was not adopted by any of those populations.
The current study employed the ribozyme B16-19, which is highly proficient in catalyzing the ligation reaction necessary for the CE to occur [47, 48], and posses an efficient folding into an active conformation [23, 25, 26]. These characteristics locate this genotype on the top of a high fitness peak [23, 26]. Thus, the mutation accumulation that occurred during each serial transfer can cause structural changes in the individuals, which in turn can further cause a decrease in the mean fitness of the population via Muller's Ratchet and random drift. The population can then fall into a fitness valley and become extinct. The addition of Mn2+ to the reaction vessels increases the mutation rate of the replication process and consequently it can alter the equilibrium distribution of the population . Quasispecies behavior sustains this shift in equilibrium distribution by a mutational coupling, allowing the population to stay extant in spite of the strong random drift and the ostensible lack of recombination in the population, as was observed recently in empirical populations of viroids .
Catalytic RNA sequences such as ligase ribozymes posses the property of buffering mutations through epistatic interactions between secondary structure arrangements [35, 49]. These arrangements strongly stabilize the structure and thus a broader range of mutations will have a neutrally selective effect, hence relaxing the error threshold [49, 50].
The fitness of each genotype in the population is normalized with the total number of genotypes in the system (assuming single locus theory applies). Thus, the proportional contribution of each genotype to the total fitness decreases as the number of genotypes increases .
The high degeneracy observed in ribozyme genotype-to-phenotype maps insures that the majority of point mutations are neutral [34, 35]. In spite of this, the mutational buffering of secondary structure epistatic interactions can only favor the fitness of the lower class mutants (e.g., low Hamming distance values). The genetic load generated in higher-class mutants will likely disrupt secondary interactions and the stability of the individual ligases. Oddly enough, an increase in the mutational rate does not cause a proportional increase in the genetic load, and therefore the population does not become extinct at a faster pace. What could be happening in this case is that mutants of lower class emerge quickly, generating a wide low-mutant class in the early evolutionary pathway of the population. These mutants have a short Hamming distances and thus probably similar fitness values. Wilke  observed that the first couple of replication cycles mostly determine fixation or extinction for an invading sequence, or perhaps for a group of close connected sequences, such as the quasispecies mutant cloud. The major contribution to the fixation probability comes from the connectivity matrix of the local genetic neighborhood of the invading sequence (or mutant class); sequences farther away on the neutral network that are poorly connected become relatively unimportant, and may be drawn out of the population by genetic drift and mutation selection balance.
The level of connection in the matrix is determined by the Hamming distance values of the mutants in the network. Genotypes that are closely connected by short Hamming distances are closely related in the sense that they can rapidly (e.g., in a few generations) be regenerated from one another in the eventual case of being removed from the population by random drift. In contrast, poorly connected genotypes (e.g., only distant relatives) will have a slow recovery into the population, if at all. In this scenario, because mutants with short Hamming distance may have close fitness values, individual sequences are not essential for the survival of the population, rather the group of close-connected individuals with mutational robustness [13, 52, 53]. Therefore, the quasispecies cloud itself is being the target of selection  and not the individual sequences. This process is analogous to the manner in which kin selection operates in animal societies [54, 55]. Ribozyme populations therefore can -- by means of indirect reproduction effects -- evolve a mutational robustness , a behavior that empowers selection with an advantage relative to other evolutionary forces (e.g., the strength that random drift has in populations of small effective sizes). This genotypic malleability allows the population to avoid a mutational meltdown, and stay extant.
These results -- of ligase RNA molecules capable of forming population structures in which cooperation is more beneficial than competition -- suggest that altruistic behavior (e.g., cooperation) is an advantageous feature to ensure survival of populations during the RNA world , when the population size were small, when the mutational rate was high, and when random genetic drift had strong effect, conditions that certainly prevailed on the prebiotic Earth [52, 58]. Additionally, quasispecies have an organization structure with the properties proposed by Kaufmann  to be necessary for the origin and preservation of genetic information. In this structure, the closely connected cloud of the quasispecies can serve as an information-preserving core, and the distantly surrounding genotypes as ideal targets for random genetic drift because they are less frequent and the information loss through them does not negatively impact the survival of population. According to Kauffman , organized systems may have arisen as a consequence of the property of some elements to establish different levels of connectivity among each other. The highly interconnected elements can create organizational cores able to preserve the information relevant to survival of the system (e.g., autocatalytic function). In contrast, less interconnected elements can serve as a reservoir of mutations without a detrimental effect on this information. Thus, during the ancient acellular times at the biogenesis on the Earth, the assemblage of information cores, perhaps in the form of quasispecies clouds, may have provided the necessary route to increase population sizes, and allow enough time for information to mature into more sophisticated functions necessary for a cellular type of life.
The quasispecies is a population structure typically formed at high mutation rates that allow the mutants to stay closely connected and thus be easily regenerated from one another even if lost from the population through random genetic drift. This behavior empowers selection relative to other evolutionary forces. Consequently, information relevant to the survival of the population can be stored in a close-knit network of mutants and not in the individuals. It is likely that such a population structure would have greatly benefited primordial pools of nascent RNA molecules on the early Earth. Instead of relying on the fortuitous advent of specific self-replicating genotypes, the RNA World would have the luxury of swarms of quasispecies evolving over time, buffered against extinction through informational decay, as theorized by Eigen .
B16-19 ligase ribozymes were freshly prepared by transcription of PCR DNA of B16-19 clones obtained in a previous in vitro evolution experiment . The RNA transcripts were purified by PAGE. The concentrations of the RNAs obtained after gel purification were measured by UV spectroscopy at 260 nm. A dilution series was then performed to obtain the desired concentration of 100 molecules in the 8.20 μL aliquot used to seed the evolution experiments.
Continuous in vitroevolution
Ligase ribozyme populations were evolved using the continuous in vitro evolution methodology [20, 25, 26, 61]. To summarize, 2.03 × 10-8nM B16-19 ligase (100 molecules) were incubated with 64 pmol of the substrate oligo S-163 (5'-CTTGACGTCAGCCTGGACTAATACGACTCAC UAUA-3' = a DNA/RNA chimera, with ribonucleotides in boldface letters, and the T7 promoter in italics), 50 pmol of TAS 1.23 primer for reverse transcription (5'-GCTGAGCCTGCGATTGG-3'), 250 units of MMLV reverse transcriptase (United States Biochemicals, Cleveland), 50 units T7 RNA polymerase (Ambion, Austin, TX), 5 nmol each dNTP, 50 nmol each rNTP, 25 mM MgCl2, and 40 μM MnCl2, in a reaction buffer with 50 mM KCl, 30 mM 4-(2-hydroxyethyl)piperazine-1-propanesulfonic acid (EPPS), pH 8.3. The 25 μL reaction mix was incubated at 37°C for exactly 22 minutes, at the completion of which the reaction was stopped by the addition of 981 μL of water. An 8.2 μL aliquot was taken from the dilution tube and used to seed the next reaction cycle of the continuous evolution process, thereby preserving a constant harmonic mean of the population size .
The survival of the evolved populations was surveyed through PCR amplifications of all the bursts used to seed a reaction cycle. The PCR products were electrophoresed through 2% agarose gels containing ethidium bromide. Visualization of the gels by trans-illumination allowed the identification of a correctly band size when the population remained alive.
Preliminary genetic variability was evaluated by RFLP using the restriction enzymes employed in the previous publication . In the populations where genotypic variability was detected, a more extensive genotypic characterization was done, as described in the following section.
Bursts were chosen for in-depth sequence analysis based on whether they displayed variability by RFLP and/or on a general desire to examine early, middle, and later bursts of 50-burst lineages. In one case (lineage 6L), two adjacent bursts were analyzed because the RFLP analysis suggested a rapid genotypic shift at that time. Specific bursts of the evolving populations with genotypic diversity were cloned using the CloneJet™ PCR Cloning Kit (Fermentas, Maryland) and E. coli competent cells (Invitrogen, San Diego). Between 60-120 colonies per burst were chosen completely at random for genotyping. Colony PCR was used to isolate the insert from single clones and further sequencing was done with BDT v3.1 chemistry. The sequences were aligned with ClustalX 2.0.11 software; the alignments were edited with BioEdit sequence alignment editor v22.214.171.124 (Tom Hall, Ibis BioSciences, Carlsbad) and the chromatogram viewer FinchTV v1.4 (Geospiza Inc., Washington).
To estimate if the number of clones sampled in each selected burst contained all the unique genotypes present in the population, rarefaction plots were constructed. Here, the pool of known genotypes is entered into the computer and then drawn in a random order by a computer (using random numbers  imported into Microsoft Excel). Then a plot is made of total new genotypes found as a function of total number of genotypes sampled, following the method used in reference . From averages of these plots, non-linear curve fitting was performed using Origin Pro v8.0 software (OriginLab Corp, Massachusetts) to give the expected asymptotes, which are estimates of the theoretical total number of genotypes present in the population .
Phylogenetic network mapping
For each lineage, all the genetic variants that were detected were aligned using DNA alignment v126.96.36.199 (Fluxus Technology Ltd.) and plotted together using the median-joining method  implemented in NETWORK v188.8.131.52 software .
We would like to thank Aaron Burton, Eric Hayden, Nilesh Vaidya, Ken Stedman, Susan Masta, and Suzanne Estes for advice during this project. This work was supported by a grant from the National Science Foundation (DEB-0315286 to NL). The funders have no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
- Eigen M: Selforganization of matter and evolution of biological macromolecules. Naturwissenschaften. 1971, 58: 465-523. 10.1007/BF00623322.View ArticlePubMedGoogle Scholar
- Eigen M, Schuster P: The hypercycle. A principle of natural self-organization. Part A: Emergence of the hypercycle. Naturwissenschaften. 1977, 64: 541-565. 10.1007/BF00450633.View ArticlePubMedGoogle Scholar
- Eigen M, Schuster P: Stages of emerging life - five principles of early organization. J Mol Evol. 1982, 19: 47-61. 10.1007/BF02100223.View ArticlePubMedGoogle Scholar
- Gilbert W: The RNA world. Nature. 1986, 319: 618-10.1038/319618a0.View ArticleGoogle Scholar
- Eigen M, McCaskill J, Schuster P: Molecular quasi-species. J Phys Chem. 1988, 92: 6881-6891. 10.1021/j100335a010.View ArticleGoogle Scholar
- Bull JJ, Meyers LA, Lachmann M: Quasispecies made simple. PLoS Comput Biol. 2005, 1: e61-10.1371/journal.pcbi.0010061.PubMed CentralView ArticlePubMedGoogle Scholar
- Eigen M: On the nature of viral quasispecies. Trends Microbiol. 1996, 4: 216-218. 10.1016/0966-842X(96)20011-3.View ArticlePubMedGoogle Scholar
- Domingo E, Sabo D, Taniguchi T, Weissman C: Nucleotide sequence heterogeneity of an RNA phage population. Cell. 1978, 13: 735-744. 10.1016/0092-8674(78)90223-4.View ArticlePubMedGoogle Scholar
- Holland J, Spindler K, Horodyski F, Grabau E, Nichol S, VandePol S: Rapid evolution of RNA genomes. Science. 1982, 215: 1577-1585. 10.1126/science.7041255.View ArticlePubMedGoogle Scholar
- Steinhauer DA, de la Torre JC, Meier E, Holland JJ: Extreme heterogeneity in populations of vesicular stomatitis virus. J Virol. 1989, 63: 2072-2080.PubMed CentralPubMedGoogle Scholar
- Eigen M: Viral quasispecies. Sci Amer. 1993, 269: 42-49. 10.1038/scientificamerican0793-42.View ArticlePubMedGoogle Scholar
- Burch CL, Chao L: Evolvability of an RNA virus is determined by its mutational neighborhood. Nature. 2000, 406: 625-628. 10.1038/35020564.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 Pathogens. 2006, 2: e136-10.1371/journal.ppat.0020136.PubMed CentralView ArticlePubMedGoogle Scholar
- Domingo E, Martín V, Perales C, Grande-Pérez A, García-Arriaza J, Arias A: Viruses as quasispecies: biological implications. Curr Topics Microbiol Immunol. 2006, 299: 51-82. full_text.Google Scholar
- Fernandez G, Bonaventura C, Martinez MA: Fitness landscape of human immunodeficiency virus type 1 protease quasispecies. J Virol. 2007, 81: 2485-2496. 10.1128/JVI.01594-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Holmes EC, Moya A: Is the quasispecies concept relevant to RNA viruses?. J Virol. 2002, 76: 463-465. 10.1128/JVI.76.1.463-465.2002.View ArticleGoogle Scholar
- Biebricher CK: Replication and evolution of short-chained RNA species replicated by Q beta replicase. Cold Spring Harbor Symp Quant Biol. 1987, 52: 299-306.View ArticlePubMedGoogle Scholar
- Biebricher CK: Quantitative analysis of mutation and selection in self-replicating RNA. Adv Space Res. 1992, 12: 191-197. 10.1016/0273-1177(92)90172-T.View ArticlePubMedGoogle Scholar
- Biebricher CK, Luce R: Sequence analysis of RNA species synthesized by Qβ replicase without template. Biochemistry. 1993, 32: 4848-4854. 10.1021/bi00069a021.View ArticlePubMedGoogle Scholar
- Wright M, Joyce GF: Continuous in vitro evolution of catalytic function. Science. 1997, 276: 614-617. 10.1126/science.276.5312.614.View ArticlePubMedGoogle Scholar
- Lynch M, Burger R, Butcher D, Gabriel W: The mutational meltdown in asexual populations. J Hered. 1993, 84: 339-344.PubMedGoogle Scholar
- Soll S, Díaz Arenas C, Lehman N: Accumulation of deleterious mutations in small abiotic populations of RNA. Genetics. 2007, 175: 267-275. 10.1534/genetics.106.066142.PubMed CentralView ArticlePubMedGoogle Scholar
- Joyce GF: Directed evolution of nucleic acid enzymes. Annu Rev Biochem. 2004, 73: 791-836. 10.1146/annurev.biochem.73.011303.073717.View ArticlePubMedGoogle Scholar
- Voytek SB, Joyce GF: Emergence of a fast-reacting ribozymes that is capable of undergoing continuous evolution. Proc Natl Acad Sci USA. 2007, 104: 15288-93. 10.1073/pnas.0707490104.PubMed CentralView ArticlePubMedGoogle Scholar
- Schmitt T, Lehman N: Non-unity heritability demonstrated by continuous evolution in vitro. Chem Biol. 1999, 6: 857-869. 10.1016/S1074-5521(00)80005-8.View ArticlePubMedGoogle Scholar
- Lehman N: Assessing the likelihood of recurrence during RNA evolution in vitro. Artificial Life. 2004, 10: 1-22. 10.1162/106454604322875887.View ArticlePubMedGoogle Scholar
- Ji J, Loeb LA: Fidelity of HIV-1 reverse transcriptase copying RNA in vitro. Biochemistry. 1992, 31: 954-958. 10.1021/bi00119a002.View ArticlePubMedGoogle Scholar
- El-Deiry W, Downey K, So A: Molecular mechanisms of manganese mutagenesis. Proc Natl Acad Sci USA. 1984, 81: 7378-7382. 10.1073/pnas.81.23.7378.PubMed CentralView ArticlePubMedGoogle Scholar
- Lazcano A, Valverde V, Hernandez G, Gariglio P, Fox G, Oró J: On the early emergence of reverse transcription: theoretical basis and experimental evidence. J Mol Evol. 1992, 35: 524-536. 10.1007/BF00160213.View ArticlePubMedGoogle Scholar
- Vartanian JP, Sala M, Henry M, Wain-Hobson S, Meyerhans A: Manganese cations increase the mutational rate of human immunodeficiency virus type 1 ex vivo. J Gen Virol. 1999, 80: 1983-1986.View ArticlePubMedGoogle Scholar
- Fisher RA: The Genetical Theory of Natural Selection. 1930, Oxford: Oxford University PressView ArticleGoogle Scholar
- Muller HJ: Some genetic aspects of sex. Amer Nat. 1932, 66: 118-138. 10.1086/280418.View ArticleGoogle Scholar
- Hill WG, Robertson A: The effect of linkage on limits to artificial selection. Genet Res. 1966, 8: 269-294. 10.1017/S0016672300010156.View ArticlePubMedGoogle Scholar
- Nimwegen E, Crutchfild J, 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
- Lehman N, Delle Donne M, West M, Dewey G: The genotypic landscape during in vitro evolution of a catalytic RNA: implication for genotypic buffering. J Mol Evol. 2000, 50: 481-490.PubMedGoogle Scholar
- Sanjuán R, Cuevas JM, Furió V, Holmes EC, Moya A: Selection for robustness in mutagenized RNA viruses. PLoS Genet. 2007, 3 (6): e93-10.1371/journal.pgen.0030093.PubMed CentralView ArticlePubMedGoogle Scholar
- Nowak MA: What is a quasispecies?. Trends Ecol Evol. 1992, 7: 118-121. 10.1016/0169-5347(92)90145-2.View ArticlePubMedGoogle Scholar
- Biebricher CK, Eigen M: What is a quasispecies?. Curr Top Microbiol Immunol. 2006, 299: 1-31. full_text.PubMedGoogle Scholar
- Schlosser K, Li Y: Diverse evolutionary trajectories characterize a community of RNA-cleaving deoxyribozymes: a case study into the population dynamics of in vitro selection. J Mol Evol. 2005, 61: 192-206. 10.1007/s00239-004-0346-7.View ArticlePubMedGoogle Scholar
- Schlosser K, Lam J, Li Y: A genotype-to-phenotype map of in vitro selected RNA-cleaving DNAzymes: implications for accessing the target phenotype. Nucleic Acids Res. 2009, 37: 3545-3557. 10.1093/nar/gkp222.PubMed CentralView ArticlePubMedGoogle Scholar
- Orgel L: Selection in vitro. Proc Roy Soc Lond. 1979, 205: 435-442. 10.1098/rspb.1979.0077.View ArticleGoogle Scholar
- Tuerk C, Gold L: Systematic evolution of ligands by exponential enrichment: RNA ligands to bacteriophage T4 DNA polymerase. Science. 1990, 249: 505-510. 10.1126/science.2200121.View ArticlePubMedGoogle Scholar
- Ellington AD, Szostak JW: In vitro selection of RNA molecules that bind specific ligands. Nature. 1990, 346: 818-822. 10.1038/346818a0.View ArticlePubMedGoogle Scholar
- Lehman N, Joyce G: Evolution in vitro: analysis of a lineage of ribozymes. Curr Biol. 1993, 3: 723-734. 10.1016/0960-9822(93)90019-K.View ArticlePubMedGoogle Scholar
- Carothers J, Szostak JW: In vitro selection of functional oligonucleotides and the origins of biochemical activity. The Aptamer Handbook: Functional Oligonucleotides and their Applications. Edited by: Klussmann S. 2006, Weinheim: Wiley-VHC Publisher, 3-28.Google Scholar
- Hanczyc MM, Dorit RL: Replicability and recurrence in the experimental evolution of a group I ribozyme. Mol Biol Evol. 2000, 17: 1050-1060.View ArticlePubMedGoogle Scholar
- Bergman N, Lau N, Lehnert V, Westhof E, Bartel D: The three-dimensional architecture of the class I ligase ribozyme. RNA. 2004, 10: 176-184. 10.1261/rna.5177504.PubMed CentralView ArticlePubMedGoogle Scholar
- Bagby SC, Bergman N, Shechner DM, Yen C, Bartel D: A class I ligase ribozyme with reduced Mg2+ dependence: Selection, sequence analysis, and identification of functional tertiary interactions. RNA. 2009, 15: 2129-2146. 10.1261/rna.1912509.PubMed CentralView ArticlePubMedGoogle Scholar
- Kun Á, Santos M, Szathmáry E: Real ribozymes suggest a relaxed error threshold. Nature Genet. 2005, 37: 1008-10.1038/ng1621.View ArticlePubMedGoogle Scholar
- Holmes E: On being the right size. Nature Genet. 2005, 37: 923-10.1038/ng0905-923.View ArticlePubMedGoogle Scholar
- Wilke CO: Probability of fixation of an advantageous mutant in a viral quasispecies. Genetics. 2003, 163: 467-474.PubMed CentralPubMedGoogle Scholar
- Kimura M: The neutral theory of molecular evolution. 1983, Great Britain: Cambridge University PressView ArticleGoogle Scholar
- Wilke CO, Wang JL, Ofria C, Lenski RE, Adami C: Evolution of digital organism at high mutational rate leads to survival of the flattest. Nature. 2001, 412: 331-333. 10.1038/35085569.View ArticlePubMedGoogle Scholar
- Maynard Smith J: Group selection and kin selection. Nature. 1964, 201: 1145-1147. 10.1038/2011145a0.View ArticleGoogle Scholar
- Maynard Smith J, Szathmáry E: The origins of life. From the Birth of Life to the Origins of Language. 1999, UK: Oxford University PressGoogle Scholar
- Hayden E, Lehman N: Self-assembly of a group I intro from inactive oligonucleotide fragments. Chem Biol. 2006, 13: 909-918. 10.1016/j.chembiol.2006.06.014.View ArticlePubMedGoogle 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
- Santos M, Zintzaras E, Szathmáry E: Recombination in primeval genomes: A step forward but still a long leap from maintaining a sizable genome. J Mol Evol. 2004, 59: 507-519. 10.1007/s00239-004-2642-7.View ArticlePubMedGoogle Scholar
- Kauffman S: The origins of order: Self-organization and selection in evolution. 1993, USA: Oxford University Press, IncGoogle Scholar
- Kauffman S: Anticaos and adaptation. Scientific American. 1991, 265: 78-84. 10.1038/scientificamerican0891-78.View ArticlePubMedGoogle Scholar
- Díaz Arenas C, Lehman N: The continuous evolution in vitro technique. Current Protocols in Nucleic Acid Chemistry. Edited by: Herdewijn P, Hoboken NJ. 2010, John Wiley & Sons, Chapter 9: Unit 9.7.1-9.7.17.Google Scholar
- Random numbers. [http://www.random.org]
- Lehman N, Wayne RK: Analysis of coyote mitochondrial DNA genotype frequencies: Estimation of the effective number of alleles. Genetics. 1991, 128: 405-416.PubMed CentralPubMedGoogle Scholar
- Gotelli N, Colwell R: Quantifying biodiversity: procedures and pitfalls in measurement and comparison of species richness. Ecology Letters. 2001, 4: 379-391. 10.1046/j.1461-0248.2001.00230.x.View ArticleGoogle Scholar
- Bandelt H-J, Forster P, Röhl A: Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999, 16: 37-48.View ArticlePubMedGoogle Scholar
- Fluxus Technology Ltd: [http://www.fluxus-engineering.com]