Emergence of linkage between cooperative RNA replicators encoding replication and metabolic enzymes through experimental evolution.

The integration of individually replicating genes into a primitive chromosome is a key evolutionary transition in the development of life, allowing the simultaneous inheritance of genes. However, how this transition occurred is unclear because the extended size of primitive chromosomes replicate slower than unlinked genes. Theoretical studies have suggested that a primitive chromosome can evolve in the presence of cell-like compartments, as the physical linkage prevents the stochastic loss of essential genes upon division, but experimental support for this is lacking. Here, we demonstrate the evolution of a chromosome-like RNA from two cooperative RNA replicators encoding replication and metabolic enzymes. Through their long-term replication in cell-like compartments, linked RNAs emerged with the two cooperative RNAs connected end-to-end. The linked RNAs had different mutation patterns than the two unlinked RNAs, suggesting that they were maintained as partially distinct lineages in the population. Our results provide experimental evidence supporting the plausibility of the evolution of a primitive chromosome from unlinked gene fragments, an important step in the emergence of complex biological systems.

Experimental evidence for the correlation between RNA structural fluctuations and the frequency of beneficial mutations.

RNA has been used as a model molecule to understand the adaptive evolution process owing to the simple relationship between the structure (i.e., phenotype) and sequence (i.e., genotype). RNA usually forms multiple substructures with similar thermodynamic stabilities, called structural fluctuations. Ancel and Fontana theoretically proposed that structural fluctuation is directly related to the ease of change in structures by mutations and thus works as a source of adaptive evolution; however, experimental verification is limited. Here, we analyzed 76 RNA genotypes that appeared in our previous in vitro evolution to examine whether (i) RNA fluctuation decreases as adaptive evolution proceeds and (ii) RNAs that have larger fluctuations tend to have higher frequencies of beneficial mutations. We first computationally estimated the structural fluctuations of all RNAs and observed that they tended to decrease as their fitness increased. We next measured the frequency of beneficial mutations for 10 RNA genotypes and observed that the total number of beneficial mutations was correlated with the size of the structural fluctuations. These results consistently support the idea that the structural fluctuation of RNA, at least those evaluated in our study, works as a source of adaptive evolution.

Plausible pathway for a host-parasite molecular replication network to increase its complexity through Darwinian evolution.

How the complexity of primitive self-replication molecules develops through Darwinian evolution remains a mystery with regards to the origin of life. Theoretical studies have proposed that coevolution with parasitic replicators increases network complexity by inducing inter-dependent replication. Particularly, Takeuchi and Hogeweg proposed a complexification process of replicator networks by successive appearance of a parasitic replicator followed by the addition of a new host replicator that is resistant to the parasitic replicator. However, the feasibility of such complexification with biologically relevant molecules is still unknown owing to the lack of an experimental model. Here, we investigated the plausible complexification pathway of host-parasite replicators using both an experimental host-parasite RNA replication system and a theoretical model based on the experimental system. We first analyzed the parameter space that allows for sustainable replication in various replication networks ranging from a single molecule to three-member networks using computer simulation. The analysis shows that the most plausible complexification pathway from a single host replicator is the addition of a parasitic replicator, followed by the addition of a new host replicator that is resistant to the parasite, consistent with the previous study by Takeuchi and Hogeweg. We also provide evidence that the pathway actually occurred in our previous evolutionary experiment. These results provide experimental evidence that a population of a single replicator spontaneously evolves into multi-replicator networks through coevolution with parasitic replicators.

Structural transition of replicable RNAs during in vitro evolution with Qß replicase.

Single-stranded RNAs (ssRNAs) are utilized as genomes in some viruses and also in experimental models of ancient life-forms, owing to their simplicity. One of the largest problems for ssRNA replication is the formation of double-stranded RNA (dsRNA), a dead-end product for ssRNA replication. A possible strategy to avoid dsRNA formation is to create strong intramolecular secondary structures of ssRNA. To design ssRNAs that efficiently replicate by Qß replicase with minimum dsRNA formation, we previously proposed the "fewer unpaired GC rule." According to this rule, ssRNAs that have fewer unpaired G and C bases in the secondary structure should efficiently replicate with less dsRNA formation. However, the validity of this rule still needs to be examined, especially for longer ssRNAs. Here, we analyze nine long ssRNAs that successively appeared during an in vitro evolution of replicable ssRNA by Qß replicase and examine whether this rule can explain the structural transitions of the RNAs. We found that these ssRNAs improved their template abilities step-by-step with decreasing dsRNA formation as mutations accumulated. We then examine the secondary structures of all the RNAs by a chemical modification method. The analysis of the structures revealed that the probabilities of unpaired G and C bases tended to decrease gradually in the course of evolution. The decreases were caused by the local structural changes around the mutation sites in most of the cases. These results support the validity of the "fewer unpaired GC rule" to efficiently design replicable ssRNAs by Qß replicase, useful for more complex ssRNA replication systems.

Exchange of Proteins in Liposomes through Streptolysin O Pores.

Liposomes, which are vesicles surrounded by lipid membranes, can be used as biochemical reactors by encapsulating various reactions. Accordingly, they are useful for studying cellular functions under controlled conditions that mimic the environment within a cell. However, one of the shortcomings of liposomes as biochemical reactors is the difficulty of introducing or removing proteins due to the impermeability of the membrane. In this study, we established a method for exchanging proteins in liposomes by forming reversible pores in the membrane. We used the toxic protein streptolysin O (SLO); this forms pores in membranes made of phospholipids containing cholesterol that can be closed by the addition of calcium ions. After optimizing the experimental procedure and lipid composition, we observed the exchange of fluorescent proteins (transferrin Alexa Fluor 488 and 647) in 9.9 % of liposomes. We also introduced T7 RNA polymerase, a 98-kDa enzyme, and observed RNA synthesis in ∼8 % of liposomes. Our findings establish a new method for controlling the internal protein composition of liposomes, thereby increasing their utility as bioreactors.

A Fusion Method to Develop an Expanded Artificial Genomic RNA Replicable by Qß Replicase.

RNA-based genomes are used to synthesize artificial cells that harbor genome replication systems. Previously, continuous replication of an artificial genomic RNA that encoded an RNA replicase was demonstrated. The next important challenge is to expand such genomes by increasing the number of encoded genes. However, technical difficulties are encountered during such expansions because the introduction of new genes disrupts the secondary structure of RNA and makes RNA nonreplicable through replicase. Herein, a fusion method that enables the construction of a longer RNA from two replicable RNAs, while retaining replication capability, is proposed. Two replicable RNAs that encode different genes at various positions are fused, and a new parameter, the unreplicable index, which negatively correlates with the replication ability of the fused RNAs better than that of the previous parameter, is determined. The unreplicable index represents the expected value of the number of G or C nucleotides that are unpaired in both the template and complementary strands. It is also observed that some of the constructed fused RNAs replicate efficiently by using the internally translated replicase. The method proposed herein could contribute to the development of an expanded RNA genome that can be used for the purpose of artificial cell synthesis.

Host-parasite oscillation dynamics and evolution in a compartmentalized RNA replication system.

To date, various cellular functions have been reconstituted in vitro such as self-replication systems using DNA, RNA, and proteins. The next important challenges include the reconstitution of the interactive networks of self-replicating species and investigating how such interactions generate complex ecological behaviors observed in nature. Here, we synthesized a simple replication system composed of two self-replicating host and parasitic RNA species. We found that the parasitic RNA eradicates the host RNA under bulk conditions; however, when the system is compartmentalized, a continuous oscillation pattern in the population dynamics of the two RNAs emerges. The oscillation pattern changed as replication proceeded mainly owing to the evolution of the host RNA. These results demonstrate that a cell-like compartment plays an important role in host-parasite ecological dynamics and suggest that the origin of the host-parasite coevolution might date back to the very early stages of the evolution of life.

A design principle for a single-stranded RNA genome that replicates with less double-strand formation.

Single-stranded RNA (ssRNA) is the simplest form of genetic molecule and constitutes the genome in some viruses and presumably in primitive life-forms. However, an innate and unsolved problem regarding the ssRNA genome is formation of inactive double-stranded RNA (dsRNA) during replication. Here, we addressed this problem by focusing on the secondary structure. We systematically designed RNAs with various structures and observed dsRNA formation during replication using an RNA replicase (Qß replicase). From the results, we extracted a simple rule regarding ssRNA genome replication with less dsRNA formation (less GC number in loops) and then designed an artificial RNA that encodes a domain of the ß-galactosidase gene based on this rule. We also obtained evidence that this rule governs the natural genomes of all bacterial and most fungal viruses presently known. This study revealed one of the structural design principles of an ssRNA genome that replicates continuously with less dsRNA formation.

Periodic Pattern of Genetic and Fitness Diversity during Evolution of an Artificial Cell-Like System.

Genetic and phenotypic diversity are the basis of evolution. Despite their importance, however, little is known about how they change over the course of evolution. In this study, we analyzed the dynamics of the adaptive evolution of a simple evolvable artificial cell-like system using single-molecule real-time sequencing technology that reads an entire single artificial genome. We found that the genomic RNA population increases in fitness intermittently, correlating with a periodic pattern of genetic and fitness diversity produced by repeated diversification and domination. In the diversification phase, a genomic RNA population spreads within a genetic space by accumulating mutations until mutants with higher fitness are generated, resulting in an increase in fitness diversity. In the domination phase, the mutants with higher fitness dominate, decreasing both the fitness and genetic diversity. This study reveals the dynamic nature of genetic and fitness diversity during adaptive evolution and demonstrates the utility of a simplified artificial cell-like system to study evolution at an unprecedented resolution.

Adaptation and Diversification of an RNA Replication System under Initiation- or Termination-Impaired Translational Conditions.

Adaptation to various environments is a remarkable characteristic of life. Is this limited to extant complex living organisms, or is it also possible for a simpler self-replication system to adapt? In this study, we addressed this question by using a translation-coupled RNA replication system that comprised a reconstituted translation system and an RNA "genome" that encoded a replicase gene. We performed RNA replication reactions under four conditions, under which different components of translation were partly inhibited. We found that replication efficiency increased with the number of rounds of replication under all the tested conditions. The types of dominant mutations differed depending on the condition, thus indicating that this simple system adapted to different environments in different ways. This suggests that even a primitive self-replication system composed of a small number of genes on the early earth could have had the ability to adapt to various environments.

Effect of Liposome Size on Internal RNA Replication Coupled with Replicase Translation.

Cell membranes inhibit the diffusion of intracellular materials, and compartment size can strongly affect the intracellular biochemical reactions. To assess the effect of the size of microcompartments on intracellular reactions, we constructed a primitive cell model consisting of giant liposomes and a translation-coupled RNA replication (TcRR) system. The RNA was replicated by Qß replicase, which was translated from the RNA in giant liposomes encapsulating the cell-free translation system. A reporter RNA encoding the antisense strand of ß-glucuronidase was introduced into the system to yield a TcRR read-out (green fluorescence). We demonstrate that TcRR was hardly detectable in larger liposomes (230 fL) but was more effective in smaller (7.7 fL) liposomes. Our experimental and theoretical results show that smaller microcompartments considerably enhance TcRR because the synthesized molecules, such as RNA and replicases, are more concentrated in smaller liposomes.

Inference of fitness values and putative appearance time points for evolvable self-replicating molecules from time series of occurrence frequencies in an evolution reactor.

We have established a translation-coupled RNA replication system within a cell-like compartment, and conducted an experimental evolution of the RNA molecules in the system. Then, we obtained a time series of occurrence frequencies of 91 individual genotypes through random sampling and next-generation sequencing. The time series showed a complex clonal interference and a polymorphic population called the "quasispecies". By fitting a deterministic kinetic model of evolvable simple self-replicators to the time series, we estimated the fitness value and "putative appearance time point" for each of the 91 major genotypes identified, where the putative appearance time point is defined as a certain time point at which a certain mutant genotype is supposed to appear in the deterministic kinetic model. As a result, the kinetic model was well fitted and additionally we confirmed that the estimated fitness values for 11 genotypes were considerably close to the experimentally measured ones (Ichihashi et al., 2015). In this sequel paper, with the theoretical basis of the deterministic kinetic model, we present the details of inference of the fitness values and putative appearance time points for the 91 genotypes. It may be possible to apply this methodology to other self-replicating molecules, viruses and bacteria.

Replication of partial double-stranded RNAs by Qß replicase.

Qß replicase, an RNA-dependent RNA polymerase of bacteriophage Qß, uses single-stranded RNA as a template to synthesize the complementary strand. A single-stranded RNA template may contain rigid secondary structures, such as long stems, intermolecular double-stranded RNA regions. Presently, the effect of the size of such double-stranded regions on the replication of RNA by Qß replicase is unknown. In this study, we prepared RNA templates hybridized with complementary RNA or DNA strands of various sizes and analyzed their replication by Qß replicase. We found that Qß replicase synthesizes the complementary strand as long as the template RNA is hybridized with no more than 200 nt fragments, although the replication amounts were decreased. This is important information to evaluate processivity of Qß replicase.

The evolutionary enhancement of genotype-phenotype linkages in the presence of multiple copies of genetic material.

Genetic evolutionary mechanisms employed by protolife developed without accompanying regulatory mechanisms for the amounts of genetic material in protocells. When many copies of genetic material are present, inactive copies generated by mutations are not effectively excluded through phenotypic selection. We demonstrate a model of gene evolution initiated with different amounts of DNA inside artificial protocells. We adopted transcription- and translation-coupled RNA replication and liposome-based in vitro compartmentalization. Despite the fact that the average number of DNA copies in each liposome was 6.4, DNA encoding active genes was maintained until the 16th selection round. Our experimental and theoretical results indicated that gene evolution can occur in the presence of multiple DNA copies. Most genetic material became junk code through gene mutations, and consequently the linkage between genotype and phenotype was enhanced through the associated decreases in active genetic material.

Kinetic analysis of aptazyme-regulated gene expression in a cell-free translation system: modeling of ligand-dependent and -independent expression.

Aptazymes are useful as RNA-based switches of gene expression responsive to several types of compounds. One of the most important properties of the switching ability is the signal/noise (S/N) ratio, i.e., the ratio of gene expression in the presence of ligand to that in the absence of ligand. The present study was performed to gain a quantitative understanding of how the aptazyme S/N ratio is determined by factors involved in gene expression, such as transcription, RNA self-cleavage, RNA degradation, protein translation, and their ligand dependencies. We performed switching of gene expression using two on-switch aptazymes with different properties in a cell-free translation system, and constructed a kinetic model that quantitatively describes the dynamics of RNA and protein species involved in switching. Both theoretical and experimental analyses consistently demonstrated that factors determining both the absolute value and the dynamics of the S/N ratio are highly dependent on the routes of translation in the absence of ligand: translation from the ligand-independently cleaved RNA or leaky translation from the noncleaved RNA. The model obtained here is useful to assess the factors that restrict the S/N ratio and to improve aptazymes more efficiently.

Mutagenic effect of ultraviolet radiation and trimethyl psoralen in Mycoplasma towards minimal genome.

Mycoplasma, solely culturable bacteria with the smallest genome, is an important organism to understand the minimal form of life. Mutagenesis using mutagens is a useful methodology for understanding the essential regions of genomic information. Ultraviolet light and trimethyl psoralen are mutagens known to induce various mutations; the latter is reported to specifically induce deletions in nematodes. However, their mutagenic effects on mycoplasma are not known. Here, we exposed Metamycoplasma salivarium to ultraviolet (UV) light or trimethyl psoralen and UV as mutagens, and analyzed the mutational pattern after several rounds of serial cultivation ranging from 34 to 56 for different lineages. Our results showed that more deletions, but fewer point mutations, were induced with TMP and UV-A than with UV alone, indicating the usefulness of TMP in inducing deletions. In addition, we compared our results with mutational data from other studies, which suggested that both TMP-UVA and UV exposure induced point mutations that were highly biased toward C to T and G to A transitions. These data provide useful basic knowledge for mutational studies on M. salivarium.

Investigation of Compatibility between DNA Replication, Transcription, and Translation for in Vitro Central Dogma.

Recent advances in in vitro synthetic biology have made it possible to reconstitute various cellular functions in a test tube. However, the integration of these functions remains a major challenge. This study aimed to identify a suitable condition to achieve all three reactions that constitute the central dogma: transcription, translation, and DNA replication. Specifically, we investigated the effect of the concentrations of 11 nonprotein factors required for in vitro transcription, translation, and DNA replication on each of these reactions. Our results indicate that certain factors have opposing effects on the three reactions. For example, while dNTP is necessary for DNA replication, it inhibited translation, and both rNTP and tRNA, which are essential for transcription and translation, inhibited DNA replication with several DNA polymerases. We also found that these opposing effects were partially alleviated by optimizing the magnesium concentration. Using this knowledge, we successfully demonstrated transcription/translation-coupled DNA replication with higher levels of transcription and translation while maintaining a certain level of DNA replication. These findings not only provide useful insights for the development of a complex artificial system with the central dogma but also raise the question of how natural cells overcome the incompatibility between the three reactions.

Minimal RNA self-reproduction discovered from a random pool of oligomers.

The emergence of RNA self-reproduction from prebiotic components would have been crucial in developing a genetic system during the origins of life. However, all known self-reproducing RNA molecules are complex ribozymes, and how they could have arisen from abiotic materials remains unclear. Therefore, it has been proposed that the first self-reproducing RNA may have been short oligomers that assemble their components as templates. Here, we sought such minimal RNA self-reproduction in prebiotically accessible short random RNA pools that undergo spontaneous ligation and recombination. By examining enriched RNA families with common motifs, we identified a 20-nucleotide (nt) RNA variant that self-reproduces via template-directed ligation of two 10 nt oligonucleotides. The RNA oligomer contains a 2'-5' phosphodiester bond, which typically forms during prebiotically plausible RNA synthesis. This non-canonical linkage helps prevent the formation of inactive complexes between self-complementary oligomers while decreasing the ligation efficiency. The system appears to possess an autocatalytic property consistent with exponential self-reproduction despite the limitation of forming a ternary complex of the template and two substrates, similar to the behavior of a much larger ligase ribozyme. Such a minimal, ribozyme-independent RNA self-reproduction may represent the first step in the emergence of an RNA-based genetic system from primordial components. Simultaneously, our examination of random RNA pools highlights the likelihood that complex species interactions were necessary to initiate RNA reproduction.

In Vitro Transcription/Translation-Coupled DNA Replication through Partial Regeneration of 20 Aminoacyl-tRNA Synthetases.

The in vitro reconstruction of life-like self-reproducing systems is a major challenge in in vitro synthetic biology. Self-reproduction requires regeneration of all molecules involved in DNA replication, transcription, and translation. This study demonstrated the continuous DNA replication and partial regeneration of major translation factors, 20 aminoacyl-tRNA synthetases (aaRS), in a reconstituted transcription/translation system (PURE system) for the first time. First, we replicated each DNA that encodes one of the 20 aaRSs through aaRS expression from the DNA by serial transfer experiments. Thereafter, we successively increased the number of aaRS genes and achieved simultaneous, continuous replication of DNA that encodes all 20 aaRSs, which comprised approximately half the number of protein factors in the PURE system, except for ribosomes, by employing dialyzed reaction and sequence optimization. This study provides a step-by-step methodology for continuous DNA replication with an increasing number of self-regenerative genes toward self-reproducing artificial systems.