Ecology and Evolution in the RNA World Dynamics and Stability of Prebiotic Replicator Systems

Group-selection via aggregative propagule-formation enables cooperative multicellularity in an individual based, spatial model

The emergence of multicellularity is one of the major transitions in evolution that happened multiple times independently. During aggregative multicellularity, genetically potentially unrelated lineages cooperate to form transient multicellular groups. Unlike clonal multicellularity, aggregative multicellular organisms do not rely on kin selection instead other mechanisms maintain cooperation against cheater phenotypes that benefit from cooperators but do not contribute to groups. Spatiality with limited diffusion can facilitate group selection, as interactions among individuals are restricted to local neighbourhoods only. Selection for larger size (e.g. avoiding predation) may facilitate the emergence of aggregation, though it is unknown, whether and how much role such selection played during the evolution of aggregative multicellularity. We have investigated the effect of spatiality and the necessity of predation on the stability of aggregative multicellularity via individual-based modelling on the ecological timescale. We have examined whether aggregation facilitates the survival of cooperators in a temporally heterogeneous environment against cheaters, where only a subset of the population is allowed to periodically colonize a new, resource-rich habitat. Cooperators constitutively produce adhesive molecules to promote aggregation and propagule-formation while cheaters spare this expense to grow faster but cannot aggregate on their own, hence depending on cooperators for long-term survival. We have compared different population-level reproduction modes with and without individual selection (predation) to evaluate the different hypotheses. In a temporally homogeneous environment without propagule-based colonization, cheaters always win. Predation can benefit cooperators, but it is not enough to maintain the necessary cooperator amount in successive dispersals, either randomly or by fragmentation. Aggregation-based propagation however can ensure the adequate ratio of cooperators-to-cheaters in the propagule and is sufficient to do so even without predation. Spatiality combined with temporal heterogeneity helps cooperators via group selection, thus facilitating aggregative multicellularity. External stress selecting for larger size (e.g. predation) may facilitate aggregation, however, according to our results, it is neither necessary nor sufficient for aggregative multicellularity to be maintained when there is effective group-selection.

Stochastic parabolic growth promotes coexistence and a relaxed error threshold in RNA-like replicator populations

The RNA world hypothesis proposes that during the early evolution of life, primordial genomes of the first self-propagating evolutionary units existed in the form of RNA-like polymers. Autonomous, non-enzymatic, and sustained replication of such information carriers presents a problem, because product formation and hybridization between template and copy strands reduces replication speed. Kinetics of growth is then parabolic with the benefit of entailing competitive coexistence, thereby maintaining diversity. Here, we test the information-maintaining ability of parabolic growth in stochastic multispecies population models under the constraints of constant total population size and chemostat conditions. We find that large population sizes and small differences in the replication rates favor the stable coexistence of the vast majority of replicator species ('genes'), while the error threshold problem is alleviated relative to exponential amplification. In addition, sequence properties (GC content) and the strength of resource competition mediated by the rate of resource inflow determine the number of coexisting variants, suggesting that fluctuations in building block availability favored repeated cycles of exploration and exploitation. Stochastic parabolic growth could thus have played a pivotal role in preserving viable sequences generated by random abiotic synthesis and providing diverse genetic raw material to the early evolution of functional ribozymes.

A Proposal of the Ur-RNAome

It is widely accepted that the earliest RNA molecules were folded into hairpins or mini-helixes. Herein, we depict the 2D and 3D conformations of those earliest RNA molecules with only RNY triplets, which Eigen proposed as the primeval genetic code. We selected 26 species (13 bacteria and 13 archaea). We found that the free energy of RNY hairpins was consistently lower than that of their corresponding shuffled controls. We found traces of the three ribosomal RNAs (16S, 23S, and 5S), tRNAs, 6S RNA, and the RNA moieties of RNase P and the signal recognition particle. Nevertheless, at this stage of evolution there was no genetic code (as seen in the absence of the peptidyl transferase centre and any vestiges of the anti-Shine-Dalgarno sequence). Interestingly, we detected the anticodons of both glycine (GCC) and threonine (GGU) in the hairpins of proto-tRNA.

The Evolution of Microbial Facilitation: Sociogenesis, Symbiogenesis, and Transition in Individuality

Metabolic cooperation is widespread, and it seems to be a ubiquitous and easily evolvable interaction in the microbial domain. Mutual metabolic cooperation, like syntrophy, is thought to have a crucial role in stabilizing interactions and communities, for example biofilms. Furthermore, cooperation is expected to feed back positively to the community under higher-level selection. In certain cases, cooperation can lead to a transition in individuality, when freely reproducing, unrelated entities (genes, microbes, etc.) irreversibly integrate to form a new evolutionary unit. The textbook example is endosymbiosis, prevalent among eukaryotes but virtually lacking among prokaryotes. Concerning the ubiquity of syntrophic microbial communities, it is intriguing why evolution has not lead to more transitions in individuality in the microbial domain. We set out to distinguish syntrophy-specific aspects of major transitions, to investigate why a transition in individuality within a syntrophic pair or community is so rare. We review the field of metabolic communities to identify potential evolutionary trajectories that may lead to a transition. Community properties, like joint metabolic capacity, functional profile, guild composition, assembly and interaction patterns are important concepts that may not only persist stably but according to thought-provoking theories, may provide the heritable information at a higher level of selection. We explore these ideas, relating to concepts of multilevel selection and of informational replication, to assess their relevance in the debate whether microbial communities may inherit community-level information or not.

RNA World Modeling: A Comparison of Two Complementary Approaches

Simple Summary

Despite years of dedicated research, scientists are still not sure what the first ”living” cell would have looked like. One of the most well-known hypotheses is the RNA world hypothesis, which assumes that, in the beginning, life relied on RNA molecules instead of DNA as information carriers and primitive enzymes. The population of such RNAs is made up of self-replicating molecules (replicases) that could make copies of themselves and parasite molecules that could only be copied by replicases. In this study, we further investigated the interplay between these hypothetical prebiotic RNA species, since it plays a crucial role in generating diversity and complexity in prebiotic molecular evolution. We compared two approaches that are commonly used to investigate such simple prebiotic systems, representing different modeling and observation scales—namely, microscopic and macroscopic. In both cases, we were able to obtain consistent results.

Abstract

The origin of life remains one of the major scientific questions in modern biology. Among many hypotheses aiming to explain how life on Earth started, RNA world is probably the most extensively studied. It assumes that, in the very beginning, RNA molecules served as both enzymes and as genetic information carriers. However, even if this is true, there are many questions that still need to be answered—for example, whether the population of such molecules could achieve stability and retain genetic information for many generations, which is necessary in order for evolution to start. In this paper, we try to answer this question based on the parasite–replicase model (RP model), which divides RNA molecules into enzymes (RNA replicases) capable of catalyzing replication and parasites that do not possess replicase activity but can be replicated by RNA replicases. We describe the aforementioned system using partial differential equations and, based on the analysis of the simulation, surmise general rules governing its evolution. We also compare this approach with one where the RP system is modeled and implemented using a multi-agent modeling technique. We show that approaching the description and analysis of the RP system from different perspectives (microscopic represented by MAS and macroscopic depicted by PDE) provides consistent results. Therefore, applying MAS does not lead to erroneous results and allows us to study more complex situations where many cases are concerned, which would not be possible through the PDE model.

Hitting Times of Some Critical Events in RNA Origins of Life

Can a replicase be found in the vast sequence space by random drift? We partially answer this question through a proof-of-concept study of the times of occurrence (hitting times) of some critical events in the origins of life for low-dimensional RNA sequences using a mathematical model and stochastic simulation studies from Python software. We parameterize fitness and similarity landscapes for polymerases and study a replicating population of sequences (randomly) participating in template-directed polymerization. Under the ansatz of localization where sequence proximity correlates with spatial proximity of sequences, we find that, for a replicating population of sequences, the hitting and establishment of a high-fidelity replicator depends critically on the polymerase fitness and sequence (spatial) similarity landscapes and on sequence dimension. Probability of hitting is dominated by landscape curvature, whereas hitting time is dominated by sequence dimension. Surface chemistries, compartmentalization, and decay increase hitting times. Compartmentalization by vesicles reveals a trade-off between vesicle formation rate and replicative mass, suggesting that compartmentalization is necessary to ensure sufficient concentration of precursors. Metabolism is thought to be necessary to replication by supplying precursors of nucleobase synthesis. We suggest that the dynamics of the search for a high-fidelity replicase evolved mostly during the final period and, upon hitting, would have been followed by genomic adaptation of genes and to compartmentalization and metabolism, effecting degree-of-freedom gains of replication channel control over domain and state to ensure the fidelity and safe operations of the primordial genetic communication system of life.

Novelty and imitation within the brain: a Darwinian neurodynamic approach to combinatorial problems

Efficient search in vast combinatorial spaces, such as those of possible action sequences, linguistic structures, or causal explanations, is an essential component of intelligence. Is there any computational domain that is flexible enough to provide solutions to such diverse problems and can be robustly implemented over neural substrates? Based on previous accounts, we propose that a Darwinian process, operating over sequential cycles of imperfect copying and selection of neural informational patterns, is a promising candidate. Here we implement imperfect information copying through one reservoir computing unit teaching another. Teacher and learner roles are assigned dynamically based on evaluation of the readout signal. We demonstrate that the emerging Darwinian population of readout activity patterns is capable of maintaining and continually improving upon existing solutions over rugged combinatorial reward landscapes. We also demonstrate the existence of a sharp error threshold, a neural noise level beyond which information accumulated by an evolutionary process cannot be maintained. We introduce a novel analysis method, neural phylogenies, that displays the unfolding of the neural-evolutionary process.

Catalytic promiscuity in the RNA World may have aided the evolution of prebiotic metabolism

The Metabolically Coupled Replicator System (MCRS) model of early chemical evolution offers a plausible and efficient mechanism for the self-assembly and the maintenance of prebiotic RNA replicator communities, the likely predecessors of all life forms on Earth. The MCRS can keep different replicator species together due to their mandatory metabolic cooperation and limited mobility on mineral surfaces, catalysing reaction steps of a coherent reaction network that produces their own monomers from externally supplied compounds. The complexity of the MCRS chemical engine can be increased by assuming that each replicator species may catalyse more than a single reaction of metabolism, with different catalytic activities of the same RNA sequence being in a trade-off relation: one catalytic activity of a promiscuous ribozyme can increase only at the expense of the others on the same RNA strand. Using extensive spatially explicit computer simulations we have studied the possibility and the conditions of evolving ribozyme promiscuity in an initial community of single-activity replicators attached to a 2D surface, assuming an additional trade-off between replicability and catalytic activity. We conclude that our promiscuous replicators evolve under weak catalytic trade-off, relatively strong activity/replicability trade-off and low surface mobility of the replicators and the metabolites they produce, whereas catalytic specialists benefit from very strong catalytic trade-off, weak activity/replicability trade-off and high mobility. We argue that the combination of conditions for evolving promiscuity are more probable to occur for surface-bound RNA replicators, suggesting that catalytic promiscuity may have been a significant factor in the diversification of prebiotic metabolic reaction networks.

Complex biochemical machineries responsible for maintaining the correct ratio of enzymes and genes were highly unlikely to exist at the wake of life. Individual genes must have been subject to competition for resources of replication leading to the competitive exclusion between them, and thus to the loss of genetic information. A feasible scenario that avoids competitive exclusion requires the assumption of mandatory cooperation between the enzymes.

A potentially dynamically important but mostly neglected feature of RNA enzymes (ribozymes) is their capacity to catalyse more than a single reaction. Here, we analyse the possibility that this “promiscous” nature of prebiotic ribozymes could have helped the maintenance of early replicator communities cooperating in running a simple metabolism. To do so, we have implemented a spatially explicit computer model simulating the dynamics of replicating entities on a mineral surface–an extension of the Metabolically Coupled Replicator System including the possibility of multiple catalytic activities within the same replicator. Our results suggest that under realistic assumptions of replicator and metabolite mobility and feasible trade-off relations between different catalytic activities of the same RNA replicator molecule, catalytic promiscuity may have indeed helped booting up life through supporting the assembly of minimal metabolisms.

Incorporating antagonistic pleiotropy into models for molecular replicators

In modern cells, chromosomal genes composed of DNA encode multi-subunit protein/RNA complexes that catalyze the replication of the chromosome and cell. One prevailing theory for the origin of life posits an early stage involving self-replicating macromolecules called replicators, which can be considered genes capable of self-replication. One prevailing theory for the genetics of aging in humans and other organisms is antagonistic pleiotropy, which posits that a gene can be beneficial in one context, and detrimental in another context. We previously reported that the conceptual simplicity of molecular replicators facilitates the generation of two simple models involving antagonistic pleiotropy. Here a third model is proposed, and each of the three models is presented with improved definition of the time variable. Computer simulations were used to calculate the proliferation of a hypothetical two-subunit replicator (AB), when one of the two subunits (B) exhibits antagonistic pleiotropy, leading to an advantage for B to be unstable. In model 1, instability of B yields free A subunits, which in turn stimulate the activity of other AB replicators. In model 2, B is lost and sometimes replaced by a more active mutant form, B'. In model 3, B becomes damaged and loses activity, and its instability allows it to be replaced by a new B. For each model, conditions were identified where instability of B was detrimental, and where instability of B was beneficial. The results are consistent with the hypothesis that antagonistic pleiotropy can promote molecular instability and system complexity, and provide further support for a model linking aging and evolution.

Evolution of linkage and genome expansion in protocells: The origin of chromosomes

Chromosomes are likely to have assembled from unlinked genes in early evolution. Genetic linkage reduces the assortment load and intragenomic conflict in reproducing protocell models to the extent that chromosomes can go to fixation even if chromosomes suffer from a replicative disadvantage, relative to unlinked genes, proportional to their length. Here we numerically show that chromosomes spread within protocells even if recurrent deleterious mutations affecting replicating genes (as ribozymes) are considered. Dosage effect selects for optimal genomic composition within protocells that carries over to the genic composition of emerging chromosomes. Lacking an accurate segregation mechanism, protocells continue to benefit from the stochastic corrector principle (group selection of early replicators), but now at the chromosome level. A remarkable feature of this process is the appearance of multigene families (in optimal genic proportions) on chromosomes. An added benefit of chromosome formation is an increase in the selectively maintainable genome size (number of different genes), primarily due to the marked reduction of the assortment load. The establishment of chromosomes is under strong positive selection in protocells harboring unlinked genes. The error threshold of replication is raised to higher genome size by linkage due to the fact that deleterious mutations affecting protocells metabolism (hence fitness) show antagonistic (diminishing return) epistasis. This result strengthens the established benefit conferred by chromosomes on protocells allowing for the fixation of highly specific and efficient enzymes.

The emergence of chromosomes harboring several genes is a crucial ingredient of the major evolutionary transition from naked replicators to cells. Linkage of replicating genes reduces conflict between them and alleviates the problem of chance loss of genes upon stochastic protocell fission. The emerging organization of protocells maintaining several segregating chromosomes with balanced gene composition also allows for an increase in the number of gene types despite recurrent deleterious mutations. We suggest that this interim genomic organization enabled protocells to evolve specific and efficient enzymes and paved the way toward an accurate mechanism for chromosome segregation later in evolution.

Multi-agent approach to sequence structure simulation in the RNA World hypothesis

The origins of life on Earth have been the subject of inquiry since the early days of philosophical thought and are still intensively investigated by the researchers around the world. One of the theories explaining the life emergence, that gained the most attention recently is the RNA World hypothesis, which assumes that life on Earth was sparked by replicating RNA chains. Since wet lab analysis is time-consuming, many mathematical and computational approaches have been proposed that try to explain the origins of life. Recently proposed one, based on the work by Takeuchi and Hogeweg, addresses the problem of interplay between RNA replicases and RNA parasitic species, which is crucial for understanding the first steps of prebiotic evolution. In this paper, the aforementioned model has been extended and modified by introducing RNA sequence (structure) information and mutation rate close to real one. It allowed to observe the simple evolution mechanisms, which could have led to the more complicated systems and eventually, to the formation of the first cells. The main goal of this study was to determine the conditions that allowed the spontaneous emergence and evolution of the prebiotic replicases equipped with simple functional domains within a large population. Here we show that polymerase ribozymes could have appeared randomly and then quickly started to copy themselves in order for the system to reach equilibrium. It has been shown that evolutionary selection works even in the simplest systems.

From self-replication to replicator systems en route to de novo life

The process by which chemistry can give rise to biology remains one of the biggest mysteries in contemporary science. The de novo synthesis and origin of life both require the functional integration of three key characteristics - replication, metabolism and compartmentalization - into a system that is maintained out of equilibrium and is capable of open-ended Darwinian evolution. This Review takes systems of self-replicating molecules as starting points and describes the steps necessary to integrate additional characteristics of life. We analyse how far experimental self-replicators have come in terms of Darwinian evolution. We also cover models of replicator communities that attempt to solve Eigen's paradox, whereby accurate replication needs complex machinery yet obtaining such complex self-replicators through evolution requires accurate replication. Successful models rely on a collective metabolism and a way of (transient) compartmentalization, suggesting that the invention and integration of these two characteristics is driven by evolution. Despite our growing knowledge, there remain numerous key challenges that may be addressed by a combined theoretical and experimental approach.

Dynamics and stability in prebiotic information integration: an RNA World model from first principles

The robust coevolution of catalytically active, metabolically cooperating prebiotic RNA replicators were investigated using an RNA World model of the origin of life based on physically and chemically plausible first principles. The Metabolically Coupled Replicator System assumes RNA replicators to supply metabolically essential catalytic activities indispensable to produce nucleotide monomers for their own template replication. Using external chemicals as the resource and the necessary ribozyme activities, Watson-Crick type replication produces complementary strands burdened by high-rate point mutations (insertions, deletions, substitutions). Metabolic ribozyme activities, replicabilities and decay rates are assigned to certain sequence and/or folding (thermodynamical) properties of single-stranded RNA molecules. Short and loosely folded sequences are given replication advantage, longer and tightly folded ones are better metabolic ribozymes and more resistant to hydrolytic decay. We show that the surface-bound MCRS evolves stable and metabolically functional communities of replicators of almost equal lengths, replicabilities and ribozyme activities. Being highly resistant to the invasion of parasitic (non-functional) replicators, it is also stable in the evolutionary sense. The template replication mechanism selects for catalytic “promiscuity”: the two (complementary) strands of the same evolved replicator will often carry more than a single catalytically active motif, thus maximizing functionality in a minimum of genetic information.

Multilevel selection as Bayesian inference, major transitions in individuality as structure learning

Complexity of life forms on the Earth has increased tremendously, primarily driven by subsequent evolutionary transitions in individuality, a mechanism in which units formerly being capable of independent replication combine to form higher-level evolutionary units. Although this process has been likened to the recursive combination of pre-adapted sub-solutions in the framework of learning theory, no general mathematical formalization of this analogy has been provided yet. Here we show, building on former results connecting replicator dynamics and Bayesian update, that (i) evolution of a hierarchical population under multilevel selection is equivalent to Bayesian inference in hierarchical Bayesian models and (ii) evolutionary transitions in individuality, driven by synergistic fitness interactions, is equivalent to learning the structure of hierarchical models via Bayesian model comparison. These correspondences support a learning theory-oriented narrative of evolutionary complexification: the complexity and depth of the hierarchical structure of individuality mirror the amount and complexity of data that have been integrated about the environment through the course of evolutionary history.

Low-Digit and High-Digit Polymers in the Origin of Life

Extant life uses two kinds of linear biopolymers that mutually control their own production, as well as the cellular metabolism and the production and homeostatic maintenance of other biopolymers. Nucleic acids are linear polymers composed of a relatively low structural variety of monomeric residues, and thus a low diversity per accessed volume. Proteins are more compact linear polymers that dispose of a huge compositional diversity even at the monomeric level, and thus bear a much higher catalytic potential. The fine-grained diversity of proteins makes an unambiguous information transfer from protein templates too error-prone, so they need to be resynthesized in every generation. But proteins can catalyse both their own reproduction as well as the efficient and faithful replication of nucleic acids, which resolves in a most straightforward way an issue termed “Eigen’s paradox”. Here the importance of the existence of both kinds of linear biopolymers is discussed in the context of the emergence of cellular life, be it for the historic orgin of life on Earth, on some other habitable planet, or in the test tube. An immediate consequence of this analysis is the necessity for translation to appear early during the evolution of life.

The Beginning of Systems Chemistry

Systems Chemistry has its roots in the research on the autocatalytic self-replication of biological macromolecules, first of all of synthetic deoxyribonucleic acids. A personal tour through the early works of the founder of Systems Chemistry, and of his first followers, recalls what's most important in this new era of chemistry: the growth and evolution of compartmented macromolecular populations, when provided with "food" and "fuel" and disposed of "waste".

Stable coevolutionary regimes for genetic parasites and their hosts: you must differ to coevolve

Background

Genetic parasites are ubiquitous satellites of cellular life forms most of which host a variety of mobile genetic elements including transposons, plasmids and viruses. Theoretical considerations and computer simulations suggest that emergence of genetic parasites is intrinsic to evolving replicator systems.

Results

Using methods of bifurcation analysis, we investigated the stability of simple models of replicator-parasite coevolution in a well-mixed environment. We first analyze what appears to be the simplest imaginable system of this type, one in which the parasite evolves during the replication of the host genome through a minimal mutation that renders the genome of the emerging parasite incapable of producing the replicase but able to recognize and recruit it for its own replication. This model has only trivial or “semi-trivial”, parasite-free equilibria: an inefficient parasite is outcompeted by the host and dies off, whereas an efficient one pushes the host out of existence, leading to the collapse of the entire system. We show that stable host-parasite coevolution (a non-trivial equilibrium) is possible in a modified model where the parasite is qualitatively distinct from the host replicator in that the replication of the parasite depends solely on the availability of the host but not on the carrying capacity of the environment.

Conclusions

We analytically determine the conditions for stable coevolution of genetic parasites and their hosts coevolution in simple mathematical models. It is shown that the evolutionary dynamics of a parasite that initially evolves from the host through the loss of the ability to replicate autonomously must substantially differ from that of the host, for a stable host-parasite coevolution regime to be established.

Electronic supplementary material

The online version of this article (10.1186/s13062-018-0230-9) contains supplementary material, which is available to authorized users.

Moderate sex between protocells can balance between a decrease in assortment load and an increase in parasite spread

Sexual reproduction is widespread in nature despite the different kinds of cost that it entails. We do not know exactly when the first sexual process took place and especially why it was beneficial at first. It is clearer why sex is advantageous for the prokaryotes and eukaryotes but the benefit of sex for protocells with individually replicating ribozymes is not yet fully understood. In this context sex is the simple horizontal gene transfer among two protocells that undergo transient fusion. Many authors argue that horizontal gene transfer (HGT) was very common in the early stage of evolution. However, HGT is a risky mechanism considering both the disruption of optimal compositions and the spread of parasites among protocells. In order to test the effects of HGT on the fitness of a protocell population, we explored by numerical simulations those conditions under which fusion might have been beneficial. We investigated multiple conceivable types of fusion in the stochastic corrector model framework and we considered the spread of parasites in every case. Protocells contain up to five species of unlinked, essential ribozymes; if a protocell has the same amount of each, it reaches maximum fitness. Fusion is dangerous not only due to the spread of parasites but also because it can ruin the cells with balanced ribozyme composition. We show that fusion can restore the ribozyme composition of the protocells under certain circumstances (high gene count, intermediate split size and low rate of fusion) and thus it can decrease the effect of the genetic load. Fusion could have been a useful early mechanism in contributing to the reliable coexistence of the different ribozymes before the spread of the chromosomes.

Mathematical modeling reveals spontaneous emergence of self-replication in chemical reaction systems

Explaining the origin of life requires us to elucidate how self-replication arises. To be specific, how can a self-replicating entity develop spontaneously from a chemical reaction system in which no reaction is self-replicating? Previously proposed mathematical models either supply an explicit framework for a minimal living system or consider only catalyzed reactions, and thus fail to provide a comprehensive theory. Here, we set up a general mathematical model for chemical reaction systems that properly accounts for energetics, kinetics, and the conservation law. We found that 1) some systems are collectively catalytic, a mode whereby reactants are transformed into end products with the assistance of intermediates (as in the citric acid cycle), whereas some others are self-replicating, that is, different parts replicate each other and the system self-replicates as a whole (as in the formose reaction, in which sugar is replicated from formaldehyde); 2) side reactions do not always inhibit such systems; 3) randomly chosen chemical universes (namely random artificial chemistries) often contain one or more such systems; 4) it is possible to construct a self-replicating system in which the entropy of some parts spontaneously decreases, in a manner similar to that discussed by Schrödinger; and 5) complex self-replicating molecules can emerge spontaneously and relatively easily from simple chemical reaction systems through a sequence of transitions. Together, these results start to explain the origins of prebiotic evolution.

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.

Sustainability of a Compartmentalized Host-Parasite Replicator System under Periodic Washout-Mixing Cycles

The emergence and dominance of parasitic replicators are among the major hurdles for the proliferation of primitive replicators. Compartmentalization of replicators is proposed to relieve the parasite dominance; however, it remains unclear under what conditions simple compartmentalization uncoupled with internal reaction secures the long-term survival of a population of primitive replicators against incessant parasite emergence. Here, we investigate the sustainability of a compartmentalized host-parasite replicator (CHPR) system undergoing periodic washout-mixing cycles, by constructing a mathematical model and performing extensive simulations. We describe sustainable landscapes of the CHPR system in the parameter space and elucidate the mechanism of phase transitions between sustainable and extinct regions. Our findings revealed that a large population size of compartments, a high mixing intensity, and a modest amount of nutrients are important factors for the robust survival of replicators. We also found two distinctive sustainable phases with different mixing intensities. These results suggest that a population of simple host-parasite replicators assumed before the origin of life can be sustained by a simple compartmentalization with periodic washout-mixing processes.