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.     

The fittest versus the flattest: experimental confirmation of the quasispecies effect with subviral pathogens

The "survival of the fittest" is the paradigm of Darwinian evolution in which the best-adapted replicators are favored by natural selection. However, at high mutation rates, the fittest organisms are not necessarily the fastest replicators but rather are those that show the greatest robustness against deleterious mutational effects, even at the cost of a low replication rate. This scenario, dubbed the "survival of the flattest", has so far only been shown to operate in digital organisms. We show that "survival of the flattest" can also occur in biological entities by analyzing the outcome of competition between two viroid species coinfecting the same plant. Under optimal growth conditions, a viroid species characterized by fast population growth and genetic homogeneity outcompeted a viroid species with slow population growth and a high degree of variation. In contrast, the slow-growth species was able to outcompete the fast species when the mutation rate was increased. These experimental results were supported by an in silico model of competing viroid quasispecies.     

The evolution of replicators

Replicators of interest in chemistry, biology and culture are briefly surveyed from a conceptual point of view. Systems with limited heredity have only a limited evolutionary potential because the number of available types is too low. Chemical cycles, such as the formose reaction, are holistic replicators since replication is not based on the successive addition of modules. Replicator networks consisting of catalytic molecules (such as reflexively autocatalytic sets of proteins, or reproducing lipid vesicles) are hypothetical ensemble replicators, and their functioning rests on attractors of their dynamics. Ensemble replicators suffer from the paradox of specificity: while their abstract feasibility seems to require a high number of molecular types, the harmful effect of side reactions calls for a small system size. No satisfactory solution to this problem is known. Phenotypic replicators do not pass on their genotypes, only some aspects of the phenotype are transmitted. Phenotypic replicators with limited heredity include genetic membranes, prions and simple memetic systems. Memes in human culture are unlimited hereditary, phenotypic replicators, based on language. The typical path of evolution goes from limited to unlimited heredity, and from attractor-based to modular (digital) replicators.     

Simple quasispecies models for the survival-of-the-flattest effect: The role of space

The survival-of-the-flattest effect postulates that under high mutation rates natural selection does not necessarily favor the faster replicators. Under such conditions, genotypes which are robust against deleterious mutational effects may be favored instead, even at the cost of a slower replication. This tantalizing hypothesis has been recently proved using digital organisms, subviral RNA plant pathogens (viroids), and an animal RNA virus. In this work we study a simple theoretical system composed by two competing quasispecies which are located at two widely different fitness landscapes that represent, respectively, a fit and a flat quasispecies. The fit quasispecies is characterized by high replication rate and low mutational robustness, whereas the flat quasispecies is characterized by low replication rate but high mutational robustness. By using a mean field model, in silico simulations with digital replicons and a two-dimensional spatial model given by a stochastic cellular automata (CA), we predict the presence of an absorbing first-order phase transition with critical slowing down between selection for replication speed and selection for mutational robustness, where the surpassing of a critical mutation rate involves the outcompetition of the fit quasispecies by the flat one. Furthermore, it is shown that space, which involves a lower critical mutation rate, broadens the conditions at which the survival-of-the-flattest may occur.     

Evolution of replicators playing a strategic game

A mathematical model of replicator evolution is considered. Replicators are words of a formal language specifying a strategy for a parlour game. They replicate with mutations and are selected according to their pay-off against other replicators.     

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 natural selection of metabolism and mass selects lifeforms from viruses to multicellular animals

I show that the natural selection of metabolism and mass can select for the major life‐history and allometric transitions that define lifeforms from viruses, over prokaryotes and larger unicells, to multicellular animals. The proposed selection is driven by a mass‐specific metabolism that is selected as the pace of the resource handling that generates net energy for self‐replication. An initial selection of mass is given by a dependence of mass‐specific metabolism on mass in replicators that are close to a lower size limit. A sublinear maximum dependence selects for virus‐like replicators, with no intrinsic metabolism, no cell, and practically no mass. A superlinear dependence selects for prokaryote‐like self‐replicating cells, with asexual reproduction and incomplete metabolic pathways. These self‐replicators have selection for increased net energy, and this generates a gradual unfolding of population‐dynamic feed‐back selection from interactive competition. The incomplete feed‐back selects for larger unicells with more developed metabolic pathways, and the completely developed feed‐back for multicellular animals with sexual reproduction. This model unifies the natural selection of lifeforms from viruses to multicellular animals, and it provides a parsimonious explanation where allometries and major life histories evolve from the natural selection of metabolism and mass.     

The extended replicator

This paper evaluates and criticises the developmental systems conception of evolution and develops instead an extension of the gene's eye conception of evolution. We argue (i) Dawkin's attempt to segregate developmental and evolutionary issues about genes is unsatisfactory. On plausible views of development it is arbitrary to single out genes as the units of selection. (ii) The genotype does not carry information about the phenotype in any way that distinguishes the role of the genes in development from that other factors. (iii) There is no simple and general causal criterion which distinguishes the role of genes in development and evolution. (iv) There is, however, an important sense in which genes but not every other developmental factor represent the phenotype. (v) The idea that genes represent features of the phenotype forces us to recognise that genes are not the only, or almost the only, replicators. Many mechanisms of replication are involved in both development and evolution. (vi) A conception of evolutionary history which recognises both genetic and non-genetic replicators, lineages of replicators and interactors has advantages over both the radical rejection of the replicator/interactor distinction and the conservative restriction of replication to genetic replication.     

Bifurcations and phase transitions in spatially extended two-member hypercycles

Mounting theoretical and experimental evidence indicates that the success of molecular replicators is strongly tied to the local nature of their interactions. Local dispersal in a given spatial domain, particularly on surfaces, might strongly enhance the growth and selection of fit molecules and their resistance to parasites. In this work the spatial dynamics of a simple hypercycle model consisting of two molecular species is analysed. In order to characterize it, both mean field models and stochastic, spatially explicit approaches are considered. The mean field approach predicts the presence of a saddle-node bifurcation separating a phase involving stable hypercycles from extinction, consistently with spatially explicit models, where an absorbing first-order phase transition is shown to exist and diffusion is explicitly introduced. The saddle-node bifurcation is shown to leave a ghost in the phase plane. A metapopulation-based model is also developed in order to account for the observed phases when both diffusion and reaction are considered. The role of information and diffusion as well as the relevance of these phases and the underlying spatial structures are discussed, and their potential implications for the evolution of early replicators are outlined.     

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.     

Eco-evolutionary dynamics,coding structure and the information threshold

Background  

The amount of information that can be maintained in an evolutionary system of replicators is limited by genome length, the number of errors during replication (mutation rate) and various external factors that influence the selection pressure. To date, this phenomenon, known as the information threshold, has been studied (both genotypically and phenotypically) in a constant environment and with respect to maintenance (as opposed to accumulation) of information. Here we take a broader perspective on this problem by studying the accumulation of information in an ecosystem, given an evolvable coding structure. Moreover, our setup allows for individual based as well as ecosystem based solutions. That is, all functions can be performed by individual replicators, or complementing functions can be performed by different replicators. In this setup, where both the ecosystem and the individual genomes can evolve their structure, we study how populations cope with high mutation rates and accordingly how the information threshold might be alleviated.     

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.     

Linear derivatives of Saccharomyces cerevisiae chromosome III can be maintained in the absence of autonomously replicating sequence elements

Replication origins in Saccharomyces cerevisiae are spaced at intervals of approximately 40 kb. However, both measurements of replication fork rate and studies of hypomorphic alleles of genes encoding replication initiation proteins suggest the question of whether replication origins are more closely spaced than should be required. We approached this question by systematically deleting replicators from chromosome III. The first significant increase in loss rate detected for the 315-kb full-length chromosome occurred only after all five efficient chromosomal replicators in the left two-thirds of the chromosome (ARS305, ARS306, ARS307, ARS309, and ARS310) had been deleted. The removal of the inefficient replicator ARS308 from this originless region caused little or no additional increase in loss rate. Chromosome fragmentations that removed the normally inactive replicators on the left end of the chromosome or the replicators distal to ARS310 on the right arm showed that both groups of replicators contribute significantly to the maintenance of the originless chromosome. Surprisingly, a 142-kb derivative of chromosome III, lacking all sequences that function as autonomously replicating sequence elements in plasmids, replicated and segregated properly 97% of the time. Both the replication initiation protein ORC and telomeres or a linear topology were required for the maintenance of chromosome fragments lacking replicators.     

Replicator Selection and the Extended Phenotype3

Adaptations are often spoken of as ‘for the good of’ some entity, but what is that entity? Groups and species are now rightly unfashionable, so what are we left with? The prevailing answer is Darwin's: ‘the individual’. Individuals clearly do not maximise their own survival so the concept of fitness had to be invented. If fitness is correctly defined in Hamilton's way as ‘inclusive fitness’ it ceases to matter whether we speak of individuals maximising their inclusive fitness or of genes maximising their survival. The two formulations are mutually inter-translatable. Yet some serious mistranslations are quoted from the literature, which have led their authors into actual biological error. The present paper blames the prevailing concentration on the individual for these errors, and advocates a reversion to the replicator as the proper focus of evolutionary attention. A gene is an obvious replicator, but there are others, and the general properties of replicators are discussed. Defenders of the individual as the unit of selection often point to the unity and integration of the genome as expressed phenotypically. This paper ends by attacking even this assumption, not by a reductionist fragmentation of the phenotype, but, on the contrary, by extending it to include more than one individual. Replicators survive by virtue of their effects on the world, and these effects are not restricted to one individual body but constitute a wider ‘extended phenotype’.     

Protein-mediated Selective Enclosure of Early Replicators Inside of Membranous Vesicles: First Step Towards Cell Membranes

Containment in cell membranes is essential for all contemporary life, and apparently even the earliest life forms had to be somehow contained. It has been postulated that random enclosure of replicating molecules inside of spontaneously assembled vesicles would have formed the initial cellular ancestors. However, completely random re-formation or division of such primitive vesicles would have abolished the heritability of their contents, nullifying any selective advantage to them. We propose that the containment of the early replicators in membranous vesicles was adopted only after the invention of genetically encoded proteins, and that selective enclosure of target molecules was mediated by specific proteins. A similar containment process is still utilised by various RNA- and retroviruses to isolate their replication complexes from the host’s intracellular environment. Such selective encapsulation would have protected the replicators against competitor and parasitic sequences, and provided a strong positive selection within the replicator communities.     

Serins with intermediate brightness have a higher survival in the wild

More ornamented individuals are supposed to signal their superior quality and health, however the development of larger ornaments or brighter plumages can affect their survival probabilities. The survival of male serins Serinus serinus in relation to the color characteristics of carotenoid-derived breast coloration was analyzed in a six-year capture–recapture study in northeastern Spain. A total of 210 males were individually marked generating 340 captures, between 1997 and 2003. The model with the best fit indicated a quadratic relationship between brightness of carotenoid-derived coloration and survival, independent of time and age of the birds, so that individuals with intermediate plumage brightness survived longer. Hence, results stress that plumage color in male serins is under stabilizing natural selection, with reduced survival both in drabber and brighter individuals. While lower condition of drabber individuals can explain their lower survival, brighter birds can suffer more from predation risk or from a higher investment of resources into mate attraction and/or reproduction.     

Extinction, colonization, and species occupancy in tidepool fishes

Despite the increasing sophistication of ecological models with respect to the size and spatial arrangement of habitat, there is relatively little empirical documentation of how species dynamics change as a function of habitat size and the fraction of habitat occupied. In an assemblage of tidepool fishes, I used maximum-likelihood estimation to test whether models which included habitat size provided a better fit to empirical data on extinction and colonization probabilities than models that assumed constant probabilities over all habitats. I found species differences in how extinction and colonization probabilities scaled with habitat size (and hence local population size). However, there was little evidence for a relationship between extinction and colonization probabilities and the fraction of occupied tidepools, as assumed in simple metapopulation models. Instead, colonization and extinction were independent of the fraction of occupied tidepools, favoring a MacArthur-Wilson island-mainland model. When I incorporated declines in extinction probability with tidepool volume in a simple simulation model, I found that predicted occupancy could change greatly, especially when colonization was low. However, the predicted fraction of occupied patches in the simulation model changed little when I incorporated the range of values reported here for extinction and colonization and the rate at which they scale with habitat size. Quantifying extinction and colonization patterns of natural populations is fundamental to understanding how species are distributed spatially and whether metapopulation models of species occupancy provide explanatory power for field populations. Received: 14 March 1997 / Accepted: 21 September 1997     

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     

Mutation induced extinction in finite populations: lethal mutagenesis and lethal isolation

Reproduction is inherently risky, in part because genomic replication can introduce new mutations that are usually deleterious toward fitness. This risk is especially severe for organisms whose genomes replicate "semi-conservatively," e.g. viruses and bacteria, where no master copy of the genome is preserved. Lethal mutagenesis refers to extinction of populations due to an unbearably high mutation rate (U), and is important both theoretically and clinically, where drugs can extinguish pathogens by increasing their mutation rate. Previous theoretical models of lethal mutagenesis assume infinite population size (N). However, in addition to high U, small N can accelerate extinction by strengthening genetic drift and relaxing selection. Here, we examine how the time until extinction depends jointly on N and U. We first analytically compute the mean time until extinction (τ) in a simplistic model where all mutations are either lethal or neutral. The solution motivates the definition of two distinct regimes: a survival phase and an extinction phase, which differ dramatically in both how τ scales with N and in the coefficient of variation in time until extinction. Next, we perform stochastic population-genetics simulations on a realistic fitness landscape that both (i) features an epistatic distribution of fitness effects that agrees with experimental data on viruses and (ii) is based on the biophysics of protein folding. More specifically, we assume that mutations inflict fitness penalties proportional to the extent that they unfold proteins. We find that decreasing N can cause phase transition-like behavior from survival to extinction, which motivates the concept of "lethal isolation." Furthermore, we find that lethal mutagenesis and lethal isolation interact synergistically, which may have clinical implications for treating infections. Broadly, we conclude that stably folded proteins are only possible in ecological settings that support sufficiently large populations.     

Replicator dominance in a eukaryotic chromosome.

Replicators are genetic elements that control initiation at an origin of DNA replication (ori). They were first identified in the yeast Saccharomyces cerevisiae as autonomously replicating sequences (ARSs) that confer on a plasmid the ability to replicate in the S phase of the cell cycle. The DNA sequences required for ARS function on a plasmid have been defined, but because many sequences that participate in ARS activity are not components of chromosomal replicators, a mutational analysis of the ARS1 replicator located on chromosome IV of S. cerevisiae was performed. The results of this analysis indicate that four DNA elements (A, B1, B2 and B3) are either essential or important for ori activation in the chromosome. In a yeast strain containing two closely spaced and identical copies of the ARS1 replicator in the chromosome, only one is active. The mechanism of replicator repression requires the essential A element of the active replicator. This element is the binding site for the origin recognition complex (ORC), a putative initiator protein. The process that determines which replicator is used, however, depends entirely upon flanking DNA sequences.