Many questions remain about the interplay between adaptive and neutral processes leading to genome expansion and the evolution of cellular complexity. Genome size appears to be tightly linked to the size of the regulatory repertoire of cells (van Nimwegen, 2006). In the context of gene regulation, we here study the interplay between adaptive and non-adaptive forces on genome and regulatory network in a computational model of cell-cycle adaptation to di_erent environments. Starting from the well-known Caulobacter crescentus network, we report on ten replicate in silico evolution experiments where cells evolve cell-cycle control by adapting to increasingly harsh spatial habitats. We find adaptive expansion of the regulatory repertoire of cells. Having a large genome is inherently costly, but also allows for improved cell-cycle behaviour. Replicates traverse difierent evolutionary trajectories leading to distinct eco-evolutionary strategies. In four replicates, cells evolve a generalist strategy to cope with a variety of nutrient levels; in two replicates, di_erent specialist cells evolve for speci_c nutrient levels; in the remaining four replicates, an intermediate strategy evolves. These diverse evolutionary outcomes reveal the role of contingency in a system under strong selective forces. This study shows that functionality of cells depends on the combination of regulatory network topology and genome organization. For example, the positions of dosage-sensitive genes are exploited to signal to the regulatory network when replication is completed, forming a de novo evolved cell-cycle checkpoint. Our results underline the importance of the integration of multiple organizational levels to understand complex gene regulation and the evolution thereof.

Eukaryotic genes are characterised by the presence of introns that are removed from the pre-mRNA by the spliceosome. This ribonucleoprotein complex is comprised of multiple RNA molecules and over a hundred proteins, which makes it one of the most complex molecular machines that originated during the prokaryote-to-eukaryote transition. Previous work has established that these introns and the spliceosomal core originated from self-splicing introns in prokaryotes. Yet it remains largely elusive how the spliceosomal core expanded by recruiting many additional proteins. In this study we use phylogenetic analyses to infer the evolutionary history of the 145 proteins that we could trace back to the spliceosome in the last eukaryotic common ancestor (LECA). We found that an overabundance of proteins derived from ribosome-related processes were added to the prokaryote-derived core. Extensive duplications of these proteins substantially increased the complexity of the emerging spliceosome. By comparing the intron positions between spliceosomal paralogs, we infer that most spliceosomal complexity postdates the spread of introns through the proto-eukaryotic genome. The reconstruction of early spliceosomal evolution provides insight into the driving forces behind the emergence of complexes with many proteins during eukaryogenesis.

Background: Convergent and parallel evolution provide unique insights into the mechanisms of natural selection. Some of the most striking convergent and parallel (collectively recurrent) amino acid substitutions in proteins are adaptive, but there are also many that are selectively neutral. Accordingly, genome-wide assessment has shown that recurrent sequence evolution in orthologs is chiefly explained by nearly neutral evolution. For paralogs, more frequent functional change is expected because additional copies are generally not retained if they do not acquire their own niche. Yet, it is unknown to what extent recurrent sequence differentiation is discernible after independent gene duplications in different eukaryotic taxa. Results: We develop a framework that detects patterns of recurrent sequence evolution in duplicated genes. This is used to analyze the genomes of 90 diverse eukaryotes. We find a remarkable number of families with a potentially predictable functional differentiation following gene duplication. In some protein families, more than ten independent duplications show a similar sequence-level differentiation between paralogs. Based on further analysis, the sequence divergence is found to be generally asymmetric. Moreover, about 6% of the recurrent sequence evolution between paralog pairs can be attributed to recurrent differentiation of subcellular localization. Finally, we reveal the specific recurrent patterns for the gene families Hint1/Hint2, Sco1/Sco2 and vma11/vma3. Conclusions: The presented methodology provides a means to study the biochemical underpinning of functional differentiation between paralogs. For instance, two abundantly repeated substitutions are identified between independently derived Sco1 and Sco2 paralogs. Such identified substitutions allow direct experimental testing of the biological role of these residues for the repeated functional differentiation. We also uncover a diverse set of families with recurrent sequence evolution and reveal trends in the functional and evolutionary trajectories of this hitherto understudied phenomenon.

The endosymbiosis of an alpha-proteobacterium that gave rise to mitochondria was one of the key events in eukaryogenesis. One striking outcome of eukaryogenesis was a much more complex cell with a large genome. Despite the existence of many alternative hypotheses for this and other patterns potentially related to endosymbiosis, a constructive evolutionary model in which these hypotheses can be studied is still lacking. Here, we present a theoretical approach in which we focus on the consequences rather than the causes of mitochondrial endosymbiosis. Using a constructive evolutionary model of cell-cycle regulation, we find that genome expansion and genome size asymmetry arise from emergent host–symbiont cell-cycle coordination. We also find that holobionts with large host and small symbiont genomes perform best on long timescales and mimic the outcome of eukaryogenesis. By designing and studying a constructive evolutionary model of obligate endosymbiosis, we uncovered some of the forces that may drive the patterns observed in nature. Our results provide a theoretical foundation for patterns related to mitochondrial endosymbiosis, such as genome size asymmetry, and reveal evolutionary outcomes that have not been considered so far, such as cell-cycle coordination without direct communication.

Molecules that replicate in trans are vulnerable to evolutionary extinction because they decrease the catalysis of replication to become more available as a template for replication. This problem can be alleviated with higher-level selection that clusters molecules of the same phenotype, favouring those groups that contain more catalysis. Here, we study a simple replicator model with implicit higher-level selection through space. We ask whether the functionality of such system can be enhanced when molecules reproduce through complementary replication, representing RNA-like replicators. For high diffusion, symmetry breaking between complementary strands occurs: one strand becomes a specialised catalyst and the other a specialised template. In ensemble, such replicators can modulate their catalytic activity depending on their environment, thereby mitigating the conflict between levels of selection. In addition, these replicators are more evolvable, facilitating survival in extreme conditions (i.e., for higher diffusion rates). Our model highlights that evolution with implicit higher-level selection—i.e., as a result of local interactions and spatial patterning—is very flexible. For different diffusion rates, different solutions to the selective conflict arise. Our results support an RNA World by showing that complementary replicators may have various ways to evolve more complexity.