Why is the polymerase chain reaction resistant to in vitro evolution?

A variety of methods have been developed to amplify DNA and RNA. These methods vary in their susceptibility to evolve new molecular species differing from the starting template. PCR is exceptionally resistant to in vitro evolution, whereas methods such as Q replicase and 3SR are much less robust. This paper develops some simple mathematical models which suggest that PCR is resistant to in vitro evolution because the reaction controls replication in discrete cycles: fast replication is of little advantage during PCR because the reaction limits fast replicators as well as slow ones to a single copy per cycle. In contrast, continuous (isothermal) reactions, as in the Q replicase reaction, favor fast replicators. The advantage of fast replication is compounded in continuous reactions, because a fast replicator can complete many generations of replication during the time it takes a slow replicator to complete one generation. These models suggest that continuous amplication protocols will never achieve the robustness against in vitro evolution observed with PCR.Correspondence to: J.J. Bull     

Models of Replicator Proliferation Involving Differential Replicator Subunit Stability

Several models for the origin of life involve molecules that are capable of self-replication, such as self-replicating polymers composed of RNA or DNA or amino acids. Here we consider a hypothetical replicator (AB) composed of two subunits, A and B. Programs written in Python and C programming languages were used to model AB replicator abundance as a function of cycles of replication (iterations), under specified hypothetical conditions. Two non-exclusive models describe how a reduced stability for B relative to A can have an advantage for replicator activity and/or evolution by generating free A subunits. In model 1, free A subunits associate with AB replicators to create AAB replicators with greater activity. In simulations, reduced stability of B was beneficial when the replication activity of AAB was greater than two times the replication activity of AB. In model 2, the free A subunit is inactive for some number of iterations before it re-creates the B subunit. A re-creates the B subunit with an equal chance of creating B or B′, where B′ is a mutant that increases AB’ replicator activity relative to AB. In simulations, at moderate number of iterations (< 15), a shorter survival time for B is beneficial when the stability of B is greater than the inactive time of A. The results are consistent with the hypothesis that reduced stability for a replicator subunit can be advantageous under appropriate conditions.     

The Role of Complex Formation and Deleterious Mutations for the Stability of RNA-Like Replicator Systems

In the RNA world hypothesis, RNA(-like) self-replicators are suggested as the central player of prebiotic evolution. However, there is a serious problem in the evolution of complexity in such replicators, i.e., the problem of parasites. Parasites, which are replicated by catalytic replicators (catalysts), but do not replicate the others, can destroy a whole replicator system by exploitation. Recently, a theoretical study underlined complex formation between replicators--an often neglected but realistic process--as a stabilizing factor in a replicator system by demonstrating that complex formation can shift the viable range of diffusion intensity to higher values. In the current study, we extend the previous study of complex formation. Firstly, by investigating a well-mixed replicator system, we establish that complex formation gives parasites an implicit advantage over catalysts, which makes the system significantly more vulnerable to parasites. Secondly, by investigating a spatially extended replicator system, we show that the formation of traveling wave patterns plays a crucial role in the stability of the system against parasites, and that because of this the effect of complex formation is not straightforward; i.e., whether complex formation stabilizes or destabilizes the spatial system is a complex function of other parameters. We give a detailed analysis of the spatial system by considering the pattern dynamics of waves. Furthermore, we investigate the effect of deleterious mutations. Surprisingly, high mutation rates can weaken the exploitation of the catalyst by the parasite.     

Gause's Principle and the Effect of Resource Partitioning on the Dynamical Coexistence of Replicating Templates

Models of competitive template replication, although basic for replicator dynamics and primordial evolution, have not yet taken different sequences explicitly into account, neither have they analyzed the effect of resource partitioning (feeding on different resources) on coexistence. Here we show by analytical and numerical calculations that Gause''s principle of competitive exclusion holds for template replicators if resources (nucleotides) affect growth linearly and coexistence is at fixed point attractors. Cases of complementary or homologous pairing between building blocks with parallel or antiparallel strands show no deviation from the rule that the nucleotide compositions of stably coexisting species must be different and there cannot be more coexisting replicator species than nucleotide types. Besides this overlooked mechanism of template coexistence we show also that interesting sequence effects prevail as parts of sequences that are copied earlier affect coexistence more strongly due to the higher concentration of the corresponding replication intermediates. Template and copy always count as one species due their constraint of strict stoichiometric coupling. Stability of fixed-point coexistence tends to decrease with the length of sequences, although this effect is unlikely to be detrimental for sequences below 100 nucleotides. In sum, resource partitioning (niche differentiation) is the default form of competitive coexistence for replicating templates feeding on a cocktail of different nucleotides, as it may have been the case in the RNA world. Our analysis of different pairing and strand orientation schemes is relevant for artificial and potentially astrobiological genetics.     

Group selection of early replicators and the origin of life

A major problem of the origin of life has been that of information integration. As Eigen (1971) has shown, a mutant distribution of RNAs replicating without the aid of a replicase cannot integrate sufficient information for the functioning of a higher-level unit utilizing several types of encoded enzymes. He proposed the hypercycle model to bridge this gap in prebiology. It can be shown by a nonlinear game model, incorporating mutation of a hypercycle, that the selection properties of hypercycles make them inefficient information integrators as they cannot compete favourably with all kinds of less efficient information carriers or mutationally coupled hypercycles. The stochastic corrector model is presented as an alternative resolution of Eigen's paradox. It assumes that replicative templates are competing within replicative compartments, whose selective values depend on the internal template composition via a catalytic acid in replication and "metabolism". The dynamics of template replication are analyzed by numerical simulation of master equations. Due to the stochasticity in replication and compartment fission the best compartment types recur. An Eigen equation at the compartment level is set up and calculated. Even selfish template mutants cannot destroy the system though they make it less efficient. The genetic information of templates is evaluated at both levels, and the higher (compartment) level successfully constrains the lower (template) one. Compartmentation together with stochastic effects is sufficient to integrate information dispersed in competitive replicators. Compartment selection is considered to be group selection of replicators. Implications for the origin of life are discussed.     

The human beta-globin replication initiation region consists of two modular independent replicators

Previous studies have shown that mammalian cells contain replicator sequences, which can determine where DNA replication initiates. However, the specific sequences that confer replicator activity were not identified. Here we report a detailed analysis of replicator sequences that dictate initiation of DNA replication from the human beta-globin locus. This analysis suggests that the beta-globin replication initiation region contains two adjacent, redundant replicators. Each replicator was capable of initiating DNA replication independently at ectopic sites. Within each of these two replicators, we identified short, discrete, nonredundant sequences, which cooperatively determine replicator activity. Experiments with somatic cell hybrids further demonstrated that the requirements for initiation at ectopic sites were similar to the requirements for initiation within native human chromosomes. The replicator clustering and redundancy exemplified in the human beta-globin locus may account for the extreme difficulty in identifying replicator sequences in mammalian cells and suggest that mammalian replication initiation sites may be determined by cooperative sequence modules.     

DNA elements modulating the <Emphasis Type="Italic">KARS12</Emphasis> chromosomal replicator in <Emphasis Type="Italic">Kluyveromyces lactis</Emphasis>

Eukaryotic chromosomal DNA replication is initiated by a highly conserved set of proteins that interact with cis-acting elements on chromosomes called replicators. Despite the conservation of replication initiation proteins, replicator sequences show little similarity from species to species in the small number of organisms that have been examined. Examination of replicators in other species is likely to reveal common features of replicators. We have examined a Kluyeromyces lactis replicator, KARS12, that functions as origin of DNA replication on plasmids and in the chromosome. It contains a 50-bp region with similarity to two other K. lactis replicators, KARS101 and the pKD1 replication origin. Replacement of the 50-bp sequence with an EcoRI site completely abrogated the ability of KARS12 to support plasmid and chromosomal DNA replication origin activity, demonstrating this sequence is a common feature of K. lactis replicators and is essential for function, possibly as the initiator protein binding site. Additional sequences up to 1 kb in length are required for efficient KARS12 function. Within these sequences are a binding site for a global regulator, Abf1p, and a region of bent DNA, both of which contribute to the activity of KARS12. These elements may facilitate protein binding, protein/protein interaction and/or nucleosome positioning as has been proposed for other eukaryotic origins of DNA replication.     

Multilevel Selection in Models of Prebiotic Evolution II: A Direct Comparison of Compartmentalization and Spatial Self-Organization

Multilevel selection has been indicated as an essential factor for the evolution of complexity in interacting RNA-like replicator systems. There are two types of multilevel selection mechanisms: implicit and explicit. For implicit multilevel selection, spatial self-organization of replicator populations has been suggested, which leads to higher level selection among emergent mesoscopic spatial patterns (traveling waves). For explicit multilevel selection, compartmentalization of replicators by vesicles has been suggested, which leads to higher level evolutionary dynamics among explicitly imposed mesoscopic entities (protocells). Historically, these mechanisms have been given separate consideration for the interests on its own. Here, we make a direct comparison between spatial self-organization and compartmentalization in simulated RNA-like replicator systems. Firstly, we show that both mechanisms achieve the macroscopic stability of a replicator system through the evolutionary dynamics on mesoscopic entities that counteract that of microscopic entities. Secondly, we show that a striking difference exists between the two mechanisms regarding their possible influence on the long-term evolutionary dynamics, which happens under an emergent trade-off situation arising from the multilevel selection. The difference is explained in terms of the difference in the stability between self-organized mesoscopic entities and externally imposed mesoscopic entities. Thirdly, we show that a sharp transition happens in the long-term evolutionary dynamics of the compartmentalized system as a function of replicator mutation rate. Fourthly, the results imply that spatial self-organization can allow the evolution of stable folding in parasitic replicators without any specific functionality in the folding itself. Finally, the results are discussed in relation to the experimental synthesis of chemical Darwinian systems and to the multilevel selection theory of evolutionary biology in general. To conclude, novel evolutionary directions can emerge through interactions between the evolutionary dynamics on multiple levels of organization. Different multilevel selection mechanisms can produce a difference in the long-term evolutionary trend of identical microscopic entities.     

The replicon revisited: an old model learns new tricks in metazoan chromosomes

The origins of DNA replication were proposed in the replicon model to be specified genetically by replicator elements that coordinate the initiation of DNA synthesis with gene expression and cell growth. Recent studies have identified DNA sequences in mammalian cells that fulfil the genetic criteria for replicators and are beginning to uncover the sequence requirements for the initiation of DNA replication. Mammalian replicators are com- posed of non-redundant modules that cooperate to direct initiation to specific chromosomal sites. Conversely, replicators do not show strong sequence similarity, and their ability to initiate replication depends on the chromosomal context and epigenetic factors, as well as their primary sequence. Here, we review the properties of metazoan replicators, and discuss the genetic and epigenetic factors that determine where and when DNA replication is initiated.     

Self-organizing proto-replicators and the origin of life

Standard theories suggest that the first informational replicator involved RNA molecules (or a more primitive analogue) and that a preliminary step for the development of such replicating systems may have been the emergence of subunits capable of forming chains and interchain pairings. Following these hypotheses, this discussion describes various abstract simulations designed to investigate the structures resulting from such interactions for generalised subunits. Three classes of pairing strategy were considered for a range of subunit concentrations. The resulting dynamic self-organization of the systems produced high levels of structural complexity (some at low subunit concentrations and in the presence of disruptive subunits) and a significantly increased percentage of complementary base pairing (particularly in the more substantial structures). These properties of the systems, which did not require pre-existing replicators, templates or functional catalysts, were shown to be sensitive to the form of pairing strategy, subunit concentrations and various conditions that could theoretically be altered by products of the systems. Though no systems behaved as a replicator, some possessed collections of properties from which a replicating system might theoretically be constructed without requiring the introduction of additional classes of properties. The implications of such systems were considered with respect to the origin of life.     

Autocatalytic networks with translation

We consider the kinetics of an autocatalytic reaction network in which replication and catalytic actions are separated by a translation step. We find that the behaviour of such a system is closely related to second-order replicator equations, which describe the kinetics of autocatalytic reaction networks in which the replicators act also as catalysts. In fact, the qualitative dynamics seems to be described almost entirely be the second-order reaction rates of the replication step. For two species we recover the qualitative dynamics of the replicator equations. Larger networks show some deviations, however. A hypercyclic system consisting of three interacting species can converge toward a stable limit cycle in contrast to the replicator equation case. A singular perturbation analysis shows that the replication-translation system reduces to a second-order replicator equation if translation is fast. The influence of mutations on replication-translation networks is also very similar to the behavior of selection-mutation equations.     

The Paradox of Dual Roles in the RNA World: Resolving the Conflict Between Stable Folding and Templating Ability

The hypothesized dual roles of RNA as both information carrier and biocatalyst during the earliest stages of life require a combination of features: good templating ability (for replication) and stable folding (for ribozymes). However, this poses the following paradox: well-folded sequences are poor templates for copying, but poorly folded sequences are unlikely to be good ribozymes. Here, we describe a strategy to overcome this dilemma through G:U wobble pairing in RNA. Unlike Watson–Crick base pairs, wobble pairs contribute highly to the energetic stability of the folded structure of their sequence, but only slightly, if at all, to the stability of the folded reverse complement. Sequences in the RNA World might thereby combine stable folding of the ribozyme with an unstructured, reverse-complementary genome, resulting in a “division of labor” between the strands. We demonstrate this strategy using computational simulations of RNA folding and an experimental model of early replication, nonenzymatic template-directed RNA primer extension. Additional study is needed to solve other problems associated with a complete replication cycle, including separation of strands after copying. Interestingly, viroid RNA sequences, which have been suggested to be relics of an RNA World (Diener, Proc Natl Acad Sci USA 86:9370–9374, 1989), also show significant asymmetry in folding energy between the infectious (+) and template (?) strands due to G:U pairing, suggesting that this strategy may even be used by replicators in the present day.     

The origin of replicators and reproducers

Replicators are fundamental to the origin of life and evolvability. Their survival depends on the accuracy of replication and the efficiency of growth relative to spontaneous decay. Infrabiological systems are built of two coupled autocatalytic systems, in contrast to minimal living systems that must comprise at least a metabolic subsystem, a hereditary subsystem and a boundary, serving respective functions. Some scenarios prefer to unite all these functions into one primordial system, as illustrated in the lipid world scenario, which is considered as a didactic example in detail. Experimentally produced chemical replicators grow parabolically owing to product inhibition. A selection consequence is survival of everybody. The chromatographized replicator model predicts that such replicators spreading on surfaces can be selected for higher replication rate because double strands are washed away slower than single strands from the surface. Analysis of real ribozymes suggests that the error threshold of replication is less severe by about one order of magnitude than thought previously. Surface-bound dynamics is predicted to play a crucial role also for exponential replicators: unlinked genes belonging to the same genome do not displace each other by competition, and efficient and accurate replicases can spread. The most efficient form of such useful population structure is encapsulation by reproducing vesicles. The stochastic corrector model shows how such a bag of genes can survive, and what the role of chromosome formation and intragenic recombination could be. Prebiotic and early evolution cannot be understood without the models of dynamics.     

A role for a replicator dominance mechanism in silencing.

The role of the natural HMR-E silencer in modulating replication initiation and silencing by the origin recognition complex (ORC) was examined. When natural HMR-E was the only silencer controlling HMR, the silencer's ORC-binding site (ACS) was dispensable for replication initiation but essential for silencing, indicating that a non-silencer chromosomal replicator(s) existed in close proximity to the silencer. Further analysis revealed that regions flanking both sides of HMR-E contained replicators. In contrast to replication initiation by the intact silencer, initiation by the non-silencer replicator(s) was abolished in an orc2-1 mutant, indicating that these replicators were extremely sensitive to defects in ORC. Remarkably, the activity of one of the non-silencer replicators correlated with reduced silencing; inactivation of these replicators caused by either the orc2-1 mutation or the deletion of flanking sequences enhanced silencing. These data were consistent with a role for the ORC bound to the HMR-E silencer ACS in suppressing the function of neighboring ORC molecules capable of inhibiting silencing, and indicated that differences in ORC-binding sites within HMR itself had profound effects on ORC function. Moreover, replication initiation by natural HMR-E was inefficient, suggesting that closely spaced replicators within HMR contributed to an inhibition of replication initiation.     

Dividing the workload at a eukaryotic replication fork

Efficient and accurate replication of the eukaryotic nuclear genome requires DNA polymerases (Pols) alpha, delta and epsilon. In all current replication fork models, polymerase alpha initiates replication. However, several models have been proposed for the roles of Pol delta and Pol epsilon in subsequent chain elongation and the division of labor between these two polymerases is still unclear. Here, we revisit this issue, considering recent studies with diagnostic mutator polymerases that support a model wherein Pol epsilon is primarily responsible for copying the leading-strand template and Pol delta is primarily responsible for copying the lagging-strand template. We also review earlier studies in light of this model and then consider prospects for future investigations of possible variations on this simple division of labor.     

The evolution of enzyme specificity in the metabolic replicator model of prebiotic evolution

The chemical machinery of life must have been catalytic from the outset. Models of the chemical origins have attempted to explain the ecological mechanisms maintaining a minimum necessary diversity of prebiotic replicator enzymes, but little attention has been paid so far to the evolutionary initiation of that diversity. We propose a possible first step in this direction: based on our previous model of a surface-bound metabolic replicator system we try to explain how the adaptive specialization of enzymatic replicator populations might have led to more diverse and more efficient communities of cooperating replicators with two different enzyme activities. The key assumptions of the model are that mutations in the replicator population can lead towards a) both of the two different enzyme specificities in separate replicators: efficient "specialists" or b) a "generalist" replicator type with both enzyme specificities working at less efficiency, or c) a fast-replicating, non-enzymatic "parasite". We show that under realistic trade-off constraints on the phenotypic effects of these mutations the evolved replicator community will be usually composed of both types of specialists and of a limited abundance of parasites, provided that the replicators can slowly migrate on the mineral surface. It is only at very weak trade-offs that generalists take over in a phase-transition-like manner. The parasites do not seriously harm the system but can freely mutate, therefore they can be considered as pre-adaptations to later, useful functions that the metabolic system can adopt to increase its own fitness.     

A conserved motif of vertebrate Y RNAs essential for chromosomal DNA replication

Noncoding Y RNAs are required for the reconstitution of chromosomal DNA replication in late G1 phase template nuclei in a human cell-free system. Y RNA genes are present in all vertebrates and in some isolated nonvertebrates, but the conservation of Y RNA function and key determinants for its function are unknown. Here, we identify a determinant of Y RNA function in DNA replication, which is conserved throughout vertebrate evolution. Vertebrate Y RNAs are able to reconstitute chromosomal DNA replication in the human cell-free DNA replication system, but nonvertebrate Y RNAs are not. A conserved nucleotide sequence motif in the double-stranded stem of vertebrate Y RNAs correlates with Y RNA function. A functional screen of human Y1 RNA mutants identified this conserved motif as an essential determinant for reconstituting DNA replication in vitro. Double-stranded RNA oligonucleotides comprising this RNA motif are sufficient to reconstitute DNA replication, but corresponding DNA or random sequence RNA oligonucleotides are not. In intact cells, wild-type hY1 or the conserved RNA duplex can rescue an inhibition of DNA replication after RNA interference against hY3 RNA. Therefore, we have identified a new RNA motif that is conserved in vertebrate Y RNA evolution, and essential and sufficient for Y RNA function in human chromosomal DNA replication.