Selfishness versus functional cooperation in a stochastic protocell model

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.

Emergence of ribozyme and tRNA-like structures from mineral-rich muddy pools on prebiotic earth

The RNA world hypothesis, although a viable one regarding the origin of life on earth, has so far failed to provide a compelling explanation for the synthesis of RNA enzymes from free nucleotides via abiotic processes. To tackle this long-standing problem, we develop a realistic model for the onset of the RNA world, using experimentally determined rates for polymerization reactions. We start with minimal assumptions about the initial state that only requires the presence of short oligomers or just free nucleotides and consider the effects of environmental cycling by dividing a day into a dry, semi-wet and wet phases that are distinguished by the nature of reactions they support. Long polymers, with maximum lengths sometimes exceeding 100 nucleotides, spontaneously emerge due to a combination of non-enzymatic, non-templated polymer extension and template-directed primer extension processes. The former helps in increasing the lengths of RNA strands, whereas the later helps in producing complementary copies of the strands. Strands also undergo hydrolysis in a structure-dependent manner that favour breaking of bonds connecting unpaired nucleotides. We identify the most favourable conditions needed for the emergence of ribozyme and tRNA-like structures and double stranded RNA molecules, classify all RNA strands on the basis of their secondary structures and determine their abundance in the population. Our results indicate that under suitable environmental conditions, non-enzymatic processes would have been sufficient to lead to the emergence of a variety of ribozyme-like molecules with complex secondary structures and potential catalytic functions.

The social coevolution hypothesis for the origin of enzymatic cooperation

At the start of life, the origin of a primitive genome required individual replicators, or genes, to act like enzymes and cooperatively copy each other. The evolutionary stability of such enzymatic cooperation poses a problem, because it would have been susceptible to parasitic replicators that did not act like enzymes but could still benefit from the enzymatic behaviour of other replicators. Existing hypotheses to solve this problem require restrictive assumptions that may not be justified, such as the evolution of a cell membrane before the evolution of enzymatic cooperation. We show theoretically that, instead, selection itself can lead to replicators grouping themselves together in a way that favours cooperation. We show that the tendency to physically associate with others and cooperative enzymatic activity can coevolve, leading to the evolution of physically linked cooperative replicators. Our results shift the empirical problem from a search for special environmental conditions to questions about what types of phenotypes can be produced by simple replicators.

Survival of RNA Replicators is much Easier in Protocells than in Surface-Based, Spatial Systems

In RNA-World scenarios for the origin of life, replication is catalyzed by polymerase ribozymes. Replicating RNA systems are subject to invasion by non-functional parasitic strands. It is well-known that there are two ways to avoid the destruction of the system by parasites: spatial clustering in models with limited diffusion, or group selection in protocells. Here, we compare computational models of replication in spatial models and protocells as closely as possible in order to determine the relative importance of these mechanisms in the RNA World. For the survival of the polymerases, the replication rate must be greater than a minimum threshold value, k_min, and the mutation rate in replication must be less than a maximum value, M_max, which is known as the error threshold. For the protocell models, we find that k_min is substantially lower and M_max is substantially higher than for the equivalent spatial models; thus, the survival of polymerases is much easier in protocells than on surfaces. The results depend on the maximum number of strands permitted in one protocell or one lattice site in the spatial model, and on whether replication is limited by the supply of monomers or the population size of protocells. The substantial advantages that are seen in the protocell models relative to the spatial models are robust to changing these details. Thus, cooperative polymerases with limited accuracy would have found it much easier to operate inside lipid compartments, and this suggests that protocells may have been a very early step in the development of life. We consider cases where parasites have an equal replication rate to polymerases, and cases where parasites multiply twice as fast as polymerases. The advantage of protocell models over spatial models is increased when the parasites multiply faster.

Natural selection in compartmentalized environment with reshuffling

The emerging field of high-throughput compartmentalized in vitro evolution is a promising new approach to protein engineering. In these experiments, libraries of mutant genotypes are randomly distributed and expressed in microscopic compartments-droplets of an emulsion. The selection of desirable variants is performed according to the phenotype of each compartment. The random partitioning leads to a fraction of compartments receiving more than one genotype making the whole process a lab implementation of the group selection. From a practical point of view (where efficient selection is typically sought), it is important to know the impact of the increase in the mean occupancy of compartments on the selection efficiency. We carried out a theoretical investigation of this problem in the context of selection dynamics for an infinite non-mutating subdivided population that randomly colonizes an infinite number of patches (compartments) at each reproduction cycle. We derive here an update equation for any distribution of phenotypes and any value of the mean occupancy. Using this result, we demonstrate that, for the linear additive fitness, the best genotype is still selected regardless of the mean occupancy. Furthermore, the selection process is remarkably resilient to the presence of multiple genotypes per compartments, and slows down approximately inversely proportional to the mean occupancy at high values. We extend out results to more general expressions that cover nonadditive and non-linear fitnesses, as well non-Poissonian distribution among compartments. Our conclusions may also apply to natural genetic compartmentalized replicators, such as viruses or early trans-acting RNA replicators.

Sustainable replication and coevolution of cooperative RNAs in an artificial cell-like system

Cooperation among independently replicating molecules is a key phenomenon that allowed the development of complexity during the early evolution of life. Generally, this process is vulnerable to parasitic or selfish entities, which can easily appear and destroy such cooperation. It remains unclear how this fragile cooperation process appeared and has been sustained through evolution. Theoretical studies have indicated that spatial structures, such as compartments, allow sustainable replication and the evolution of cooperative replication, although this has yet to be confirmed experimentally. In this study, we constructed a molecular cooperative replication system, in which two types of RNA, encoding replication or metabolic enzymes, cooperate for their replication in compartments, and we performed long-term replication experiments to examine the sustainability and evolution of the RNAs. We demonstrate that the cooperative relationship of the two RNAs could be sustained at a certain range of RNA concentrations, even when parasitic RNA appeared in the system. We also found that more efficient cooperative RNA replication evolved during long-term replication through seemingly selfish evolution of each RNA. Our results provide experimental evidence supporting the sustainability and robustness of molecular cooperation on an evolutionary timescale.

The evolutionary dynamics of metabolic protocells

Protocell multilevel selection models have been proposed to study the evolutionary dynamics of vesicles encapsulating a set of replicating, competing and mutating sequences. The frequency of the different sequence types determines protocell survival through a fitness function. One of the defining features of these models is the genetic load generated when the protocell divides and its sequences are assorted between the offspring vesicles. However, these stochastic assortment effects disappear when the redundancy of each sequence type is sufficiently high. The fitness dependence of the vesicle with its sequence content is usually defined without considering a realistic account on how the lower level dynamics would specify the protocell fitness. Here, we present a protocell model with a fitness function determined by the output flux of a simple metabolic network with the aim of understanding how the evolution of both kinetic and topological features of metabolism would have been constrained by the particularities of the protocell evolutionary dynamics. In our model, the sequences inside the vesicle are both the carriers of information and Michaelis-Menten catalysts exhibiting saturation. We found that the saturation of the catalysts controlling the metabolic fluxes, achievable by modifying the kinetic or stoichiometric parameters, provides a mechanism to ameliorate the assortment load by increasing the redundancy of the catalytic sequences required to achieve the maximum flux. Regarding the network architecture, we conclude that combinations of parallel network motifs and bimolecular catalysts are a robust way to increase the complexity of the metabolism enclosed by the protocell.

The protocell hypothesis conjectures the existence of a vesicle containing catalytic and replicating sequences as the primordial cellular organization during the early stages of the evolution of life. Mathematical models of protocells traditionally consider RNA sequences being encapsulated and having both an informational and a catalytic role in the same molecule. Because of this dual function, the protocell sequences are evolutionarily constrained. Mathematical models have been extensively used to study the evolutionary dynamics of protocells with a focus on the processes, like mutation or stochastic sequence assortment upon division, that affect the protocell information capacity in terms of the coexistence of different sequence types. Here we introduce a simple model of metabolic networks whose output determines the survival of the protocell with the aim of studying the effect of modifying the kinetic and architectural properties of the network on sequence coexistence. We find that stochastic assortment and mutation limit the architectures able to be encapsulated by the protocell with a given fraction of the population harbouring all possible sequence types.

Co-operation between Polymerases and Nucleotide Synthetases in the RNA World

It is believed that life passed through an RNA World stage in which replication was sustained by catalytic RNAs (ribozymes). The two most obvious types of ribozymes are a polymerase, which uses a neighbouring strand as a template to make a complementary sequence to the template, and a nucleotide synthetase, which synthesizes monomers for use by the polymerase. When a chemical source of monomers is available, the polymerase can survive on its own. When the chemical supply of monomers is too low, nucleotide production by the synthetase is essential and the two ribozymes can only survive when they are together. Here we consider a computational model to investigate conditions under which coexistence and cooperation of these two types of ribozymes is possible. The model considers six types of strands: the two functional sequences, the complementary strands to these sequences (which are required as templates), and non-functional mutants of the two sequences (which act as parasites). Strands are distributed on a two-dimensional lattice. Polymerases replicate strands on neighbouring sites and synthetases produce monomers that diffuse in the local neighbourhood. We show that coexistence of unlinked polymerases and synthetases is possible in this spatial model under conditions in which neither sequence could survive alone; hence, there is a selective force for increasing complexity. Coexistence is dependent on the relative lengths of the two functional strands, the strand diffusion rate, the monomer diffusion rate, and the rate of deleterious mutations. The sensitivity of this two-ribozyme system suggests that evolution of a system of many types of ribozymes would be difficult in a purely spatial model with unlinked genes. We therefore speculate that linkage of genes onto mini-chromosomes and encapsulation of strands in protocells would have been important fairly early in the history of life as a means of enabling more complex systems to evolve.

Trans-acting polymerases are cooperative, because they copy neighbouring strands, and do not copy themselves directly. Inaccurate replication creates parasitic strands that act as templates but not ribozymes. It is known that in spatially distributed models with slow strand diffusion, clusters of cooperating polymerases arise that can survive in the presence of parasites provided that the error rate is less than a maximum limit (the error threshold). In the RNA World, we envisage multiple types of ribozymes working together. We would like to understand how a multi-ribozyme system could evolve from a system with a single type of polymerase ribozyme. As a first step in increasing complexity, we consider a two-ribozyme system in which there is one polymerase and one nucleotide synthetase that produces monomers for use by the polymerase. We are particularly interested to find conditions in which the chemical supply of monomers is too low for the polymerase to survive alone, but the additional monomers created by the synthetase allow the two-ribozyme system to survive where the single-ribozyme system could not. There is then a selective force for increasing the complexity of the system. Here we show that spatial clustering is sufficient to allow cooperation and survival of systems of unlinked ribozymes with different functions. Clusters form in which synthetases form fringes around the polymerases. Survival of the two-ribozyme system depends on several factors. The strand diffusion rate must be slow enough for cooperative clusters to emerge. The replication rate of the polymerase must be comparable to that of the synthetase. The diffusion rate of the monomers must be neither too slow nor too fast. The model considers the most difficult case for cooperation–unlinked genes with no compartments. The sensitivity of the two-ribozyme system that we study here suggests that evolution of a spatial system with multiple unlinked ribozymes would become increasingly more difficult as the number of components increased, and suggests that linkage and protocells would need to evolve relatively early in the history of life.

Maximal gene number maintainable by stochastic correction - The second error threshold

There is still no general solution to Eigen׳s Paradox, the chicken-or-egg problem of the origin of life: neither accurate copying, nor long genomes could have evolved without one another being established beforehand. But an array of small, individually replicating genes might offer a workaround, provided that multilevel selection assists the survival of the ensemble. There are two key difficulties that such a system has to overcome: the non-synchronous replication of genes, and their random assortment into daughter cells (the units of higher-level selection) upon fission. Here we find, using the Stochastic Corrector Model framework, that a large number (τ≥90) of genes can coexist. Furthermore, the system can tolerate about 10% replication rate asymmetry (competition) among the genes. On this basis, we put forward a plausible (and testable!) scenario for how novel genes could have been incorporated into early living systems: a route to complex metabolism.

The dynamics of the RNA world: insights and challenges

The RNA world hypothesis of the origin of life, in which RNA emerged as both enzyme and information carrier, is receiving solid experimental support. The prebiotic synthesis of biomolecules, the catalytic aid offered by mineral surfaces, and the vast enzymatic repertoire of ribozymes are only pieces of the origin of life puzzle; the full picture can only emerge if the pieces fit together by either following from one another or coexisting with each other. Here, we review the theory of the origin, maintenance, and enhancement of the RNA world as an evolving population of dynamical systems. The dynamical view of the origin of life allows us to pinpoint the missing and the not fitting pieces: (1) How can the first self-replicating ribozyme emerge in the absence of template-directed information replication? (2) How can nucleotide replicators avoid competitive exclusion despite utilizing the very same resources (nucleobases)? (3) How can the information catastrophe be avoided? (4) How can enough genes integrate into a cohesive system in order to transition to a cellular stage? (5) How can the way information is stored and metabolic complexity coevolve to pave to road leading out of the RNA world to the present protein-DNA world?

Computer simulations of cellular group selection reveal mechanism for sustaining cooperation

We present a computer simulation of group selection that is inspired by proto-cell division. Two types of replicating molecules, cooperators and defectors, reside inside membrane bound compartments. Cooperators pay a cost for other replicators in the cell to receive a benefit. Defectors pay no cost and distribute no benefits. The total population size fluctuates as a consequence of births and deaths of individual replicators. Replication requires activated substrates that are generated at a constant rate. Our model includes mutation between cooperators and defectors and selection on two levels: within proto-cells and between proto-cells. We find surprising similarities and differences between models with and without cell death. In both cases, a necessary requirement for group selection to favor some level of cooperation is the continuous formation of a minimum fraction of pure cooperator groups. Subsequently these groups become undermined by defectors, because of mutation and selection within cells. Cell division mechanisms which generate pure cooperator groups more efficiently are stronger promoters of cooperation. For example, division of a proto-cell into many daughter cells is more powerful in enhancing cooperation than division into two daughter cells. Our model differs from previous studies of group selection in that we explore a variety of different features and relax several restrictive assumptions that would be needed for analytic calculations.

Early evolution of efficient enzymes and genome organization

Background

Cellular life with complex metabolism probably evolved during the reign of RNA, when it served as both information carrier and enzyme. Jensen proposed that enzymes of primordial cells possessed broad specificities: they were generalist. When and under what conditions could primordial metabolism run by generalist enzymes evolve to contemporary-type metabolism run by specific enzymes?

Results

Here we show by numerical simulation of an enzyme-catalyzed reaction chain that specialist enzymes spread after the invention of the chromosome because protocells harbouring unlinked genes maintain largely non-specific enzymes to reduce their assortment load. When genes are linked on chromosomes, high enzyme specificity evolves because it increases biomass production, also by reducing taxation by side reactions.

Conclusion

The constitution of the genetic system has a profound influence on the limits of metabolic efficiency. The major evolutionary transition to chromosomes is thus proven to be a prerequisite for a complex metabolism. Furthermore, the appearance of specific enzymes opens the door for the evolution of their regulation.

Reviewers

This article was reviewed by Sándor Pongor, Gáspár Jékely, and Rob Knight.

Computer simulation on the cooperation of functional molecules during the early stages of evolution

It is very likely that life began with some RNA (or RNA-like) molecules, self-replicating by base-pairing and exhibiting enzyme-like functions that favored the self-replication. Different functional molecules may have emerged by favoring their own self-replication at different aspects. Then, a direct route towards complexity/efficiency may have been through the coexistence/cooperation of these molecules. However, the likelihood of this route remains quite unclear, especially because the molecules would be competing for limited common resources. By computer simulation using a Monte-Carlo model (with “micro-resolution” at the level of nucleotides and membrane components), we show that the coexistence/cooperation of these molecules can occur naturally, both in a naked form and in a protocell form. The results of the computer simulation also lead to quite a few deductions concerning the environment and history in the scenario. First, a naked stage (with functional molecules catalyzing template-replication and metabolism) may have occurred early in evolution but required high concentration and limited dispersal of the system (e.g., on some mineral surface); the emergence of protocells enabled a “habitat-shift” into bulk water. Second, the protocell stage started with a substage of “pseudo-protocells”, with functional molecules catalyzing template-replication and metabolism, but still missing the function involved in the synthesis of membrane components, the emergence of which would lead to a subsequent “true-protocell” substage. Third, the initial unstable membrane, composed of prebiotically available fatty acids, should have been superseded quite early by a more stable membrane (e.g., composed of phospholipids, like modern cells). Additionally, the membrane-takeover probably occurred at the transition of the two substages of the protocells. The scenario described in the present study should correspond to an episode in early evolution, after the emergence of single “genes”, but before the appearance of a “chromosome” with linked genes.

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.

At the core of all biological systems lies the division of labor between the storage of genetic information and its phenotypic implementation, in other words, the functional differentiation between templates (DNA) and catalysts (proteins). This fundamental property of life is believed to have been absent at the earliest stages of evolution. The RNA world hypothesis, the most realistic current scenario for the origin of life, posits that, in primordial replicating systems, RNA functioned both as template and as catalyst. How would such division of labor emerge through Darwinian evolution? We investigated the evolution of DNA-like molecules in minimal computational models of RNA replicator systems. Two models were considered: one where molecules are adsorbed on surfaces and another one where molecules are compartmentalized by dividing cellular boundaries. Both models exhibit the evolution of DNA and the ensuing division of labor, revealing the simple governing principle of these processes: DNA releases RNA from the trade-off between template and catalyst that is inevitable in the RNA world and thereby enhances the system's resistance against parasitic templates. Hence, this study offers a novel insight into the evolutionary origin of the division of labor between templates and catalysts in the RNA world.