On the origin of DNA genomes: evolution of the division of labor between template and catalyst in model replicator systems

The division of labor between template and catalyst is a fundamental property of all living systems: DNA stores genetic information whereas proteins function as catalysts. The RNA world hypothesis, however, posits that, at the earlier stages of evolution, RNA acted as both template and catalyst. Why would such division of labor evolve in the RNA world? We investigated the evolution of DNA-like molecules, i.e. molecules that can function only as template, in minimal computational models of RNA replicator systems. In the models, RNA can function as both template-directed polymerase and template, whereas DNA can function only as template. Two classes of models were explored. In the surface models, replicators are attached to surfaces with finite diffusion. In the compartment models, replicators are compartmentalized by vesicle-like boundaries. Both models displayed the evolution of DNA and the ensuing division of labor between templates and catalysts. In the surface model, DNA provides the advantage of greater resistance against parasitic templates. However, this advantage is at least partially offset by the disadvantage of slower multiplication due to the increased complexity of the replication cycle. In the compartment model, DNA can significantly delay the intra-compartment evolution of RNA towards catalytic deterioration. These results are explained in terms of the trade-off between template and catalyst that is inherent in RNA-only replication cycles: DNA releases RNA from this trade-off by making it unnecessary for RNA to serve as template and so rendering the system more resistant against evolving parasitism. Our analysis of these simple models suggests that the lack of catalytic activity in DNA by itself can generate a sufficient selective advantage for RNA replicator systems to produce DNA. Given the widespread notion that DNA evolved owing to its superior chemical properties as a template, this study offers a novel insight into the evolutionary origin of DNA.     

Spatial Models of Prebiotic Evolution: Soup Before Pizza?

The problem of information integration andresistance to the invasion of parasitic mutants in prebiotic replicator systemsis a notorious issue of research on the origin of life.Almost all theoretical studies published so far havedemonstrated that some kind of spatial structure is indispensable forthe persistence and/or the parasite resistance of any feasible replicator system.Based on a detailed critical survey of spatial models on prebiotic informationintegration, we suggest a possible scenario for replicator system evolution leadingto the emergence of the first protocells capable of independent life.We show that even the spatial versions of the hypercycle model are vulnerable toselfish parasites in heterogeneous habitats. Contrary, the metabolic system remainspersistent and coexistent with its parasites both on heterogeneous surfaces andin chaotically mixing flowing media. Persistent metabolic parasites can beconverted to metabolic cooperators, or they can gradually obtain replicase activity.Our simulations show that, once replicase activity emerged, a gradual and simultaneousevolutionary improvement of replicase functionality (speed and fidelity) andtemplate efficiency is possible only on a surface that constrains the mobility ofmacromolecule replicators. Based on the results of the models reviewed, we suggestthat open chaotic flows (`soup') and surface dynamics (`pizza') both played keyroles in the sequence of evolutionary events ultimately concluding in theappearance of the first living cell on Earth.     

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.     

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.     

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.     

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 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.     

Selfishness versus functional cooperation in a stochastic protocell model

How to design an “evolvable” artificial system capable to increase in complexity? Although Darwin’s theory of evolution by natural selection obviously offers a firm foundation, little hope of success seems to be expected from the explanatory adequacy of modern evolutionary theory, which does a good job at explaining what has already happened but remains practically helpless at predicting what will occur. However, the study of the major transitions in evolution clearly suggests that increases in complexity have occurred on those occasions when the conflicting interests between competing individuals were partly subjugated. This immediately raises the issue about “levels of selection” in evolutionary biology, and the idea that multi-level selection scenarios are required for complexity to emerge. After analyzing the dynamical behaviour of competing replicators within compartments, we show here that a proliferation of differentiated catalysts and/or improvement of catalytic efficiency of ribozymes can potentially evolve in properly designed artificial cells where the strong internal competition between the different species of replicators is somewhat prevented (i.e., by choosing them with equal probability). Experimental evolution in these systems will likely stand as beautiful examples of artificial adaptive systems, and will provide new insights to understand possible evolutionary paths to the evolution of metabolic complexity.     

Analysis of chromosome III replicators reveals an unusual structure for the ARS318 silencer origin and a conserved WTW sequence within the origin recognition complex binding site

Saccharomyces cerevisiae chromosome III encodes 11 autonomously replicating sequence (ARS) elements that function as chromosomal replicators. The essential 11-bp ARS consensus sequence (ACS) that binds the origin recognition complex (ORC) has been experimentally defined for most of these replicators but not for ARS318 (HMR-I), which is one of the HMR silencers. In this study, we performed a comprehensive linker scan analysis of ARS318. Unexpectedly, this replicator depends on a 9/11-bp match to the ACS that positions the ORC binding site only 6 bp away from an Abf1p binding site. Although a largely inactive replicator on the chromosome, ARS318 becomes active if the nearby HMR-E silencer is deleted. We also performed a multiple sequence alignment of confirmed replicators on chromosomes III, VI, and VII. This analysis revealed a highly conserved WTW motif 17 to 19 bp from the ACS that is functionally important and is apparent in the 228 phylogenetically conserved ARS elements among the six sensu stricto Saccharomyces species.     

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     

Pattern selection in a cooperative biochemical in vitro amplification system: the role of parasites

Previously, numerical simulations have shown that evolving systems can be stabilized against emerging parasites by pattern formation in spatially extended flow reactors. Hence, it can be argued that pattern formation is a prerequisite for any experimental investigation of the biochemical evolution of cooperative function. Here, we study a model of an experimental biochemical system for the cooperative in vitro amplification of DNA strands and show that emerging parasites can induce a complex pattern formation even when no pattern formation occurs without parasites. In an adiabatic approximation where the cooperative amplification reaction is assumed to adapt fast to slowly emerging parasites, the parasite concentration itself acts as a Steuer parameter for the selection of various complex patterns. Without such an adiabatic approximation only transient patterns emerge. As any species can grow for very low concentrations, the parasite is able to infect the entire reactor and the system is finally diluted out. In the experimental biochemical system, however, the species are individual molecules and the growth of spatially separated, non-infected regions becomes feasible. Hence a cutoff threshold for the minimal concentration is applied. In these simulations the otherwise lethal infection by parasites induces the formation of spatiotemporal spirals, and this spatial structure help the host and parasitoid species to survive together. These theoretical results describe an inherent property of cooperative reactions and have an important impact on experimental investigations on the molecular evolution and complex function in spatially extended reactors. Since the formation of the complex pattern is restricted either to a rather large cutoff value or a special choice of the kinetic parameters, we, however, conclude that the persistence of evolving cooperative amplification is not possible in a simple reaction-diffusion reactor. Experimental set-ups with patchy environments, e.g. biomolecular amplification in coupled microstructured flow chambers or in microemulsion, are eligible candidates for the observation of such a self-organized pattern selection.     

Multilevel Selection in Models of Prebiotic Evolution: Compartments and Spatial Self-organization

In this paper we explore the impact of new levels ofselection in models of early evolution.We contrast two types of higher levels of selection. On the one hand we look at spatially explicitmodels of replicators in which, by a process of self-organization, new levels of selection ariseas large scale spatial patterns with a dynamics of their own. Alternatively externally imposedlevels of selection above the basic replicators are created by enclosing the replicators in vesicles.In this paper we first review some results on the impact of emerging higher levels of selection onthe evolutionary persistence of interacting co-evolvingreplicator systems. Moreover, we presenta vesicle model, which can potentially integrate emerging and imposed levels of selection. We use the models to examinethe classical problem information integration in early evolution.Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1025754907141     

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.     

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.     

Less can be more: RNA-adapters may enhance coding capacity of replicators

It is still not clear how prebiotic replicators evolved towards the complexity found in present day organisms. Within the most realistic scenario for prebiotic evolution, known as the RNA world hypothesis, such complexity has arisen from replicators consisting solely of RNA. Within contemporary life, remarkably many RNAs are involved in modifying other RNAs. In hindsight, such RNA-RNA modification might have helped in alleviating the limits of complexity posed by the information threshold for RNA-only replicators. Here we study the possible role of such self-modification in early evolution, by modeling the evolution of protocells as evolving replicators, which have the opportunity to incorporate these mechanisms as a molecular tool. Evolution is studied towards a set of 25 arbitrary 'functional' structures, while avoiding all other (misfolded) structures, which are considered to be toxic and increase the death-rate of a protocell. The modeled protocells contain a genotype of different RNA-sequences while their phenotype is the ensemble of secondary structures they can potentially produce from these RNA-sequences. One of the secondary structures explicitly codes for a simple sequence-modification tool. This 'RNA-adapter' can block certain positions on other RNA-sequences through antisense base-pairing. The altered sequence can produce an alternative secondary structure, which may or may not be functional. We show that the modifying potential of interacting RNA-sequences enables these protocells to evolve high fitness under high mutation rates. Moreover, our model shows that because of toxicity of misfolded molecules, redundant coding impedes the evolution of self-modification machinery, in effect restraining the evolvability of coding structures. Hence, high mutation rates can actually promote the evolution of complex coding structures by reducing redundant coding. Protocells can successfully use RNA-adapters to modify their genotype-phenotype mapping in order to enhance the coding capacity of their genome and fit more information on smaller sized genomes.     

Differential targeting of Tetrahymena ORC to ribosomal DNA and non‐rDNA replication origins

The Tetrahymena thermophila origin recognition complex (ORC) contains an integral RNA subunit, 26T RNA, which confers specificity to the amplified ribosomal DNA (rDNA) origin by base pairing with an essential cis‐acting replication determinant—the type I element. Using a plasmid maintenance assay, we identified a 6.7 kb non‐rDNA fragment containing two closely associated replicators, ARS1‐A (0.8 kb) and ARS1‐B (1.2 kb). Both replicators lack type I elements and hence complementarity to 26T RNA, suggesting that ORC is recruited to these sites by an RNA‐independent mechanism. Consistent with this prediction, although ORC associated exclusively with origin sequences in the 21 kb rDNA minichromosome, the interaction between ORC and the non‐rDNA ARS1 chromosome changed across the cell cycle. In G₂ phase, ORC bound to all tested sequences in a 60 kb interval spanning ARS1‐A/B. Remarkably, ORC and Mcm6 associated with just the ARS1‐A replicator in G₁ phase when pre‐replicative complexes assemble. We propose that ORC is stochastically deposited onto newly replicated non‐rDNA chromosomes and subsequently targeted to preferred initiation sites prior to the next S phase.     

Spatial Dynamics and the Evolution of Enzyme Production

We re-examine the problem of the evolution of protein synthesis or enzyme production using a stochastic cellular automaton model, where the replicators are fixed in the sites of a two-dimensional square lattice. In contrast with the classical chemical kinetics or mean-field predictions, we show that a small colony of mutant, protein-mediated (enzymatic) replicators has an appreciable probability to take over a resident population of simpler, direct-template replicators. In addition, we argue that the threshold phenomenon corresponding to the onset of invasion can be described quantitatively within the physics framework of nonequilibrium phase transitions. We study also the invasion of a resident population of enzymatic replicators by more efficient replicators of the same kind, and show that although slightly more efficient mutants cannot invade, invasion is a likely event if the productivity advantage of the mutants is large. In this sense, the establishment of a population of enzymatic replicators is not a `once-forever' evolutionary decision.     

Two neo-Darwinisms

There are two extant theories of evolution, each of which deserves the honourific "neo-Darwinism": Modern Synthesis Replicator theory and a theory I shall call Developmental Darwinism. The principal difference concerns the canonical unit of biological organization. Modern Synthesis replicator theory explains the process of evolution by appeal to the activities of genes or replicators. Developmental Darwinism explains the process of evolution by appeal to the capacities of organisms. In particular, it is the plasticity of organisms, manifested most distinctly during development, that causes adaptive evolution. Despite the fact that each, in its own way, traces its origin to the theory outlined by Darwin, they are radically different. The objectives of this essay are twofold: to underscore the differences between these theories, and to argue that Developmental Darwinism, though nascent, is a viable alternative to Modern Synthesis replicator theory.     

Selection without replicators: the origin of genes, and the replicator/interactor distinction in etiobiology

Genes are thought to have evolved from long-lived and multiply-interactive molecules in the early stages of the origins of life. However, at that stage there were no replicators, and the distinction between interactors and replicators did not yet apply. Nevertheless, the process of evolution that proceeded from initial autocatalytic hypercycles to full organisms was a Darwinian process of selection of favourable variants. We distinguish therefore between Neo-Darwinian evolution and the related Weismannian and Central Dogma divisions, on the one hand, and the more generic category of Darwinian evolution on the other. We argue that Hull’s and Dawkins’ replicator/interactor distinction of entities is a sufficient, but not necessary, condition for Darwinian evolution to take place. We conceive the origin of genes as a separation between different types of molecules in a thermodynamic state space, and employ a notion of reproducers.