BACKGROUND
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  <cit> .
under an rna world scenario, the molecules in question were autocatalytically replicating polymers of nucleotides, such as rna or a chemical equivalent  <cit> .
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  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  <cit> .
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  <cit> .
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  <cit> , it soon was argued that many viral populations evolved with the same dynamic.
the first empirical demonstration of quasispecies behavior was produced in  <dig> with the qβ phage  <cit> .
since then many viral populations have been described as evolving as quasispecies collectives under mutation-selection balance  <cit> , although this interpretation has been debated in the case of viruses  <cit> .
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  <cit> .
such experiments confirm some of the key predictions of quasispecies dynamics as outlined by eigen  <cit> : 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 , 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  in vitro evolution  <cit>  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  <cit>  of small populations when the mutation rate is limited by the intrinsic error rate of viral polymerases  <cit> .
now, with the error rate greatly enhanced with the addition of a chemical mutagen  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.
RESULTS
ce experiments
four clonal 100-molecule populations of the rna class i ligase ribozyme b16- <dig>  were each evolved independently using the ce system.
the ce protocol is a means to induce the rapid evolution of ligase ribozymes using moloney murine leukaemia virus reverse transcriptase  and t <dig> rna polymerase to sustain rna populations through sequential serial transfers  <cit> .
each serial transfer involves roughly three cycles of amplification that produces a rapid proliferation of rna molecules, and hence is termed "burst" in this paper.
the experimental conditions used were the same as in earlier experiments  <cit> , with the exception that mncl <dig> was added to the reaction vessel at a final concentration of  <dig> μm to increase the mutational rate.
the in vitro per nucleotide error rate of mmlv-rt has been estimated at about 1/ <dig>   <cit> , but mn2+ ions lower the substrate specificity of rna- and dna-dependent dna polymerases both in vivo and in vitro, resulting in a 6-30-fold higher error rate  <cit> .
all four lineages  were carried out for  <dig> bursts without a sign of population decay via muller's ratchet, and thus a mutational meltdown was never observed .
these results were unlike those of previous data obtained in the absence of mncl <dig>  in which a meltdown was observed at an average of  <dig>  bursts  <cit> .
in fact, the data collected here were directly compared to a more extended analysis of one lineage from the previous study  that actually did survive to burst  <dig>  <cit> .
genotypic characterization
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  two lineages  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  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  <cit>  without added mgcl <dig>  were genotyped for comparison.
network analysis
alignment of the nucleotide sequencing data for lineages 6h and 6l showed a trend in the population dynamics in which a majority of the clones have the same genotype, while a minority have slightly different ones.
this observation was the first clue that quasispecies behavior may be present in these ligase populations.
phylogenetic networks were drawn  to find the genetic relationships among the mutants in both lineages, and hence the structures of each putative quasispecies  <cit> .
the structures of these networks show a dynamic characterized by a dominant and centered sequence that is present at a high frequency, around which are located other, less frequent mutant sequences.
this population structure is indeed characteristic of a quasispecies, in which the dominant sequence is called the "master sequence" and the surrounding mutants form the mutant cloud  <cit> .
these quasispecies are characterized by a relatively close connection among the mutants as indicated by the hamming distances calculated within bursts  based on the network diagrams .
the close connectivity between these mutants is suggestive of a mutational resiliency in the population that is a plausible cause of their observed extended persistence times; they did not go extinct prior to burst  <dig> unlike all the 100-molecule lineages examined previously  <cit> .
key: lineage 3d is from a previous study  <cit> , while lineages 6h and 6l are from this study.
the master sequence  names are given by order of appearance, simply  <dig>   <dig>  etc., except in the case of the wildtype , b16- <dig>  the frequency of the ms, the total number of sequences sampled, and the number of unique genotypes was calculated by the network software, and the shannon diversity index values were calculated using the formula given in the text.
for comparison,  <dig> and  <dig> clones, respectively, from bursts  <dig> and  <dig> of the smallest lineage  that survived in the absence of added mncl <dig> were genotyped.
unlike lineages 6h and 6l, the network diagrams of these bursts failed to exhibit a dominant master sequence .
instead, there appeared to be a greater number of less-common genotypes competing for existence, as would be expected under more typical conditions of directional selection and/or random genetic drift.
in addition to the network diagrams, one manner in which lineage 3d could be compared to the quasispecies in lineages 6h and 6l is to enumerate the frequency of the most common genotype.
in lineage 3d this value did not exceed 42%, while in lineages 6h and 6l this value dropped below 50% only once in seven cases: in burst  <dig> of lineage 6h , where the network diagrams indicated that a transition between one master sequence and another was in progress .
also the average hamming distance within the 3d and 3l lineages was  <dig> , significantly higher than in 6h or 6l .
evolutionary dynamics
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 .
these results indicate that each lineage can develop a different dynamic over the course of  <dig> 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  shows that a quasispecies-like cloud is already formed , in which the master sequence is still the "wildtype" b16- <dig>  with a frequency of 76%.
this cloud was constructed from a sample of  <dig> clones, which contained  <dig> different genotypes .
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- <dig> clones per burst were genotyped , and yet the harmonic mean population size of these bursts was ~ <dig> individuals and the frequency distributions of genotypes were highly skewed , 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 .
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  <dig> 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,  <dig> clones from burst  <dig> 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  <dig> , as measured by shannon diversity index  where pi is the frequency of the ith clone).
deeper examination of the bursts in this lineage  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 .
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  <dig> 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  <dig>  b16- <dig> ; burst  <dig>  ms <dig> ; burst  <dig>  ms <dig> .
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  <cit> .
a network drawn by combining the sequences of all bursts inspected in lineage 6h  shows the fluctuation of the equilibrium dynamic of the quasispecies in its evolutionary history.
in lineage 6l, four bursts were evaluated .
these specific bursts were chosen because the rflp analysis detected genotypic transitions near these time points .
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  reveal that this fluctuation can occur in a relatively short time .
interestingly, some of the master sequences that appear during this lineage  are the same as were observed in lineage 6h .
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 .
these dynamics of staggered dominant genotypes that fluctuate as the population evolves  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 , 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  <cit> .
some specific mutations, such as the u62a "insurance" mutation  <cit> , 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 , which conferred resistance to mutational meltdown in small population sizes in the absence of mncl <dig>  was not seen in the current study, although short insertions, deletions, and rearrangements were common in the last 10- <dig> nt at the 3' end of the ligase , and specific mutations in this region appeared in more than one quasispecies.
recurrence of mutations is actually a predicted characteristic of quasispecies  <cit> , thus providing additional support for this interpretation of these data.
in general, from all the networks drawn, it was observed that the number of mutations that appeared in more than one burst and/or lineage constituted 61% of the total number of sequences explored, and the great majority of these mutants became part of master sequences during the lineages' evolutionary history .
the fact that most of the mutants that have evolved in these lineages were able to persist in time could be a consequence of the mutational resiliency that evolved owing to the short hamming distance values .
most of these recurrent mutations belonged within the master sequences .
in lineage 6h, the initial change in master sequences from b16- <dig> to ms <dig> implies ten changes, and the further change from ms <dig> to ms <dig> implies one change .
similarly, in lineage 6l, the initial change from b16- <dig> to ms <dig> implies sixteen changes, but further changes in the master sequence transitions are seventeen and one .
the hamming distance values are generally low , eighteen being the maximum value.
the distance between the two most recurrent master sequences is actually only one, and these master sequences  are the most representative in time and sequence space .
the hamming distance values, in addition to the high frequency of recurrent mutations, support the idea of a form of mutational robustness  <cit>  evolving in the system through a quasispecies behavior in ligase populations evolved in vitro, although a formal test of this will require a comparison of fitness values.
discussion
"what is a quasispecies?"
is the exact title of at least two papers  <cit> .
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- <dig> ligase ribozymes were evolved using the continuous in vitro evolution  method  <cit>  and the relatively error-prone mmlv-rt.
in particular, mncl <dig> 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 mncl <dig>  <cit> .
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  <cit> .
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 v <dig> of the qβ virus performed by orgel and co-workers  <cit> .
in this case, the mutagen ethidium bromide  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  <cit> .
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  <cit> , and the group i ribozymes  <cit> , but quasispecies-like behavior was not adopted by any of those populations.
the current study employed the ribozyme b16- <dig>  which is highly proficient in catalyzing the ligation reaction necessary for the ce to occur  <cit> , and posses an efficient folding into an active conformation  <cit> .
these characteristics locate this genotype on the top of a high fitness peak  <cit> .
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  <cit> .
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  <cit> .
the shift in the equilibrium distribution is possible because of two main reasons, first stated and then explained.
catalytic rna sequences such as ligase ribozymes posses the property of buffering mutations through epistatic interactions between secondary structure arrangements  <cit> .
these arrangements strongly stabilize the structure and thus a broader range of mutations will have a neutrally selective effect, hence relaxing the error threshold  <cit> .
the fitness of each genotype in the population is normalized with the total number of genotypes in the system .
thus, the proportional contribution of each genotype to the total fitness decreases as the number of genotypes increases  <cit> .
the high degeneracy observed in ribozyme genotype-to-phenotype maps insures that the majority of point mutations are neutral  <cit> .
in spite of this, the mutational buffering of secondary structure epistatic interactions can only favor the fitness of the lower class mutants .
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  <cit>  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 ; 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  be regenerated from one another in the eventual case of being removed from the population by random drift.
in contrast, poorly connected genotypes  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  <cit> .
therefore, the quasispecies cloud itself is being the target of selection  <cit>  and not the individual sequences.
this process is analogous to the manner in which kin selection operates in animal societies  <cit> .
ribozyme populations therefore can -- by means of indirect reproduction effects -- evolve a mutational robustness  <cit> , a behavior that empowers selection with an advantage relative to other evolutionary forces .
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  is an advantageous feature to ensure survival of populations during the rna world  <cit> , 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  <cit> .
additionally, quasispecies have an organization structure with the properties proposed by kaufmann  <cit>  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  <cit> , 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 .
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.
CONCLUSIONS
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  <cit> .

 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. under an rna world scenario, the molecules in question were autocatalytically replicating polymers of nucleotides, such as rna or a chemical equivalent. 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 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, it soon was argued that many viral populations evolved with the same dynamic. the first empirical demonstration of quasispecies behavior was produced in with the qβ phage. since then many viral populations have been described as evolving as quasispecies collectives under mutation-selection balance, although this interpretation has been debated in the case of viruses. 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. such experiments confirm some of the key predictions of quasispecies dynamics as outlined by eigen: 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, 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 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 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.  ce experiments four clonal 100-molecule populations of the rna class i ligase ribozyme b16- were each evolved independently using the ce system. the ce protocol is a means to induce the rapid evolution of ligase ribozymes using moloney murine leukaemia virus reverse transcriptase and t rna polymerase to sustain rna populations through sequential serial transfers. each serial transfer involves roughly three cycles of amplification that produces a rapid proliferation of rna molecules, and hence is termed "burst" in this paper. the experimental conditions used were the same as in earlier experiments, with the exception that mncl was added to the reaction vessel at a final concentration of μm to increase the mutational rate. the in vitro per nucleotide error rate of mmlv-rt has been estimated at about 1/, but mn2+ ions lower the substrate specificity of rna- and dna-dependent dna polymerases both in vivo and in vitro, resulting in a 6-30-fold higher error rate. all four lineages were carried out for bursts without a sign of population decay via muller's ratchet, and thus a mutational meltdown was never observed. these results were unlike those of previous data obtained in the absence of mncl in which a meltdown was observed at an average of bursts. in fact, the data collected here were directly compared to a more extended analysis of one lineage from the previous study that actually did survive to burst. genotypic characterization 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 two lineages 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 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 mgcl were genotyped for comparison. network analysis alignment of the nucleotide sequencing data for lineages 6h and 6l showed a trend in the population dynamics in which a majority of the clones have the same genotype, while a minority have slightly different ones. this observation was the first clue that quasispecies behavior may be present in these ligase populations. phylogenetic networks were drawn to find the genetic relationships among the mutants in both lineages, and hence the structures of each putative quasispecies. the structures of these networks show a dynamic characterized by a dominant and centered sequence that is present at a high frequency, around which are located other, less frequent mutant sequences. this population structure is indeed characteristic of a quasispecies, in which the dominant sequence is called the "master sequence" and the surrounding mutants form the mutant cloud. these quasispecies are characterized by a relatively close connection among the mutants as indicated by the hamming distances calculated within bursts based on the network diagrams. the close connectivity between these mutants is suggestive of a mutational resiliency in the population that is a plausible cause of their observed extended persistence times; they did not go extinct prior to burst unlike all the 100-molecule lineages examined previously. key: lineage 3d is from a previous study, while lineages 6h and 6l are from this study. the master sequence names are given by order of appearance, simply etc., except in the case of the wildtype, b16- the frequency of the ms, the total number of sequences sampled, and the number of unique genotypes was calculated by the network software, and the shannon diversity index values were calculated using the formula given in the text. for comparison, and clones, respectively, from bursts and of the smallest lineage that survived in the absence of added mncl were genotyped. unlike lineages 6h and 6l, the network diagrams of these bursts failed to exhibit a dominant master sequence. instead, there appeared to be a greater number of less-common genotypes competing for existence, as would be expected under more typical conditions of directional selection and/or random genetic drift. in addition to the network diagrams, one manner in which lineage 3d could be compared to the quasispecies in lineages 6h and 6l is to enumerate the frequency of the most common genotype. in lineage 3d this value did not exceed 42%, while in lineages 6h and 6l this value dropped below 50% only once in seven cases: in burst of lineage 6h, where the network diagrams indicated that a transition between one master sequence and another was in progress. also the average hamming distance within the 3d and 3l lineages was, significantly higher than in 6h or 6l. evolutionary dynamics 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. these results indicate that each lineage can develop a different dynamic over the course of 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 shows that a quasispecies-like cloud is already formed, in which the master sequence is still the "wildtype" b16- with a frequency of 76%. this cloud was constructed from a sample of clones, which contained different genotypes. 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- clones per burst were genotyped, and yet the harmonic mean population size of these bursts was ~ individuals and the frequency distributions of genotypes were highly skewed, 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. 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 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, clones from burst 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, as measured by shannon diversity index where pi is the frequency of the ith clone). deeper examination of the bursts in this lineage 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. 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 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 b16-; burst ms; burst ms. 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. a network drawn by combining the sequences of all bursts inspected in lineage 6h shows the fluctuation of the equilibrium dynamic of the quasispecies in its evolutionary history. in lineage 6l, four bursts were evaluated. these specific bursts were chosen because the rflp analysis detected genotypic transitions near these time points. 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 reveal that this fluctuation can occur in a relatively short time. interestingly, some of the master sequences that appear during this lineage are the same as were observed in lineage 6h. 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. these dynamics of staggered dominant genotypes that fluctuate as the population evolves 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, 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, which conferred resistance to mutational meltdown in small population sizes in the absence of mncl was not seen in the current study, although short insertions, deletions, and rearrangements were common in the last 10- nt at the 3' end of the ligase, 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. in general, from all the networks drawn, it was observed that the number of mutations that appeared in more than one burst and/or lineage constituted 61% of the total number of sequences explored, and the great majority of these mutants became part of master sequences during the lineages' evolutionary history. the fact that most of the mutants that have evolved in these lineages were able to persist in time could be a consequence of the mutational resiliency that evolved owing to the short hamming distance values. most of these recurrent mutations belonged within the master sequences. in lineage 6h, the initial change in master sequences from b16- to ms implies ten changes, and the further change from ms to ms implies one change. similarly, in lineage 6l, the initial change from b16- to ms implies sixteen changes, but further changes in the master sequence transitions are seventeen and one. the hamming distance values are generally low, eighteen being the maximum value. the distance between the two most recurrent master sequences is actually only one, and these master sequences are the most representative in time and sequence space. the hamming distance values, in addition to the high frequency of recurrent mutations, support the idea of a form of mutational robustness evolving in the system through a quasispecies behavior in ligase populations evolved in vitro, although a formal test of this will require a comparison of fitness values. discussion "what is a quasispecies?" is the exact title of at least two papers. 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- ligase ribozymes were evolved using the continuous in vitro evolution method and the relatively error-prone mmlv-rt. in particular, mncl 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 mncl. 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. 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 v of the qβ virus performed by orgel and co-workers. in this case, the mutagen ethidium bromide 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. 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- which is highly proficient in catalyzing the ligation reaction necessary for the ce to occur, and posses an efficient folding into an active conformation. these characteristics locate this genotype on the top of a high fitness peak. 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. the shift in the equilibrium distribution is possible because of two main reasons, first stated and then explained. catalytic rna sequences such as ligase ribozymes posses the property of buffering mutations through epistatic interactions between secondary structure arrangements. these arrangements strongly stabilize the structure and thus a broader range of mutations will have a neutrally selective effect, hence relaxing the error threshold. the fitness of each genotype in the population is normalized with the total number of genotypes in the system. 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. in spite of this, the mutational buffering of secondary structure epistatic interactions can only favor the fitness of the lower class mutants. 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; 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 be regenerated from one another in the eventual case of being removed from the population by random drift. in contrast, poorly connected genotypes 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. 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. 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. 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 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. 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. 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. 
using continuous in vitro evolution experimentation with an increased mutation rate imposed by mncl it was found that clonal 100-molecule populations of ribozymes clearly exhibit certain characteristics of a quasispecies. 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. extensive genotypic sampling from two replicate lineages was gathered and phylogenetic networks were constructed to elucidate the structure of the evolving rna populations. during the rna world, molecular populations were probably very small and highly susceptible to the force of strong random drift. the quasispecies nature of mutant rna clouds facilitates the recovery of genotypes under pressure of being removed from the population by random drift. most of the mutants in these clouds were connected by short hamming distance values, indicating their close relatedness. 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 the first time this has been seen with a catalytic rna. in conjunction with muller's ratchet, this would have imposed difficulties for the preservation of the genetic information and the survival of the populations.
 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 mncl 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. 