Early Eukaryotes Research Articles

Mechanisms of axoneme and centriole elimination in Naegleria gruberi.

The early branching eukaryote Naegleria gruberi can transform transiently from an amoeboid life form lacking centrioles and flagellato a flagellate life form where these elements are present, followed by reversion to the amoeboid state. The mechanisms imparting elimination of axonemes and centrioles during this reversion process are not known. Here, we uncover that flagella primarily fold onto the cell surface and fuse within milliseconds with the plasma membrane. Once internalized, axonemes are severed by Spastin into similarly-sized fragments that are then enclosed by membranes, before their contents are eliminated through the lysosomal pathway. Moreover, we discovered that centrioles undergo progressive K63 autophagy-linked poly-ubiquitination and K48 proteasome-promoting poly-ubiquitination, and that such ubiquitination occurs next to centriolar microtubules. Most centrioles are eliminated in either lysosomes or the cytoplasm in a lysosomal- and proteasome-dependent manner. Strikingly, we uncover in addition that centrioles can be shed in the extracellular milieu and taken up by other cells. Collectively, these findings reveal fundamental mechanisms governing the elimination of essential cellular constituents in Naegleria that may operate broadly in eukaryotic systems.

Genome-wide quantification of polycistronic transcription in Leishmania major.

Leishmania major is a human-pathogenic, obligate parasite and the etiological agent of the most prevalent, cutaneous form of leishmaniasis, which is an important neglected, tropical disease with ~1.2 million new infections per year. Leishmania, and the whole order Trypanosomatida, are early eukaryotes with highly diverged gene expression and regulation pathways, setting them apart from their mammalian hosts and from most other eukaryotes. Using precision run-on sequence analysis, we performed a genome-wide mapping and density analysis of RNA polymerases in isolated nuclei of the protozoan parasite Leishmania major. We map transcription initiation sites at divergent strand switch regions and head-tail regions within the chromosomes and correlate them with known sites of chromatin modifications. We confirm continuous, polycistronic RNA synthesis in all RNA polymerase II-dependent gene arrays but find small varying RNA polymerase activities in polycistronic transcription units (PTUs), excluding gene-specific transcription regulation, but not PTU-specific variations. Lastly, we find evidence for transcriptional pausing of all three RNA polymerase classes, hinting at a possible mechanism of transcriptional regulation.IMPORTANCELeishmania spp. are pathogens of humans and animals and cause one of the most important neglected tropical diseases. Regulation of gene expression in Leishmania but also in the related Trypanosoma is radically different from all eukaryotic model organisms, dispensing with regulated, gene-specific transcription, and relying instead on highly regulated translation. Our work sheds light on the initiation, elongation, and termination of transcription, maps unidirectional, polycistronic transcription units, provides evidence for transcriptional pausing at or near starting points of RNA synthesis, and quantifies the varying transcription rates of the polycistronic transcription units. Our results will further the understanding of these important pathogens and should provide a valuable resource for researchers in the field of eukaryotic microbiology.

NuSAP4 regulates chromosome segregation in Trypanosoma brucei by promoting bipolar spindle assembly

Faithful chromosome segregation in eukaryotes requires the assembly of a bipolar spindle and the faithful attachment of kinetochores to spindle microtubules, which are regulated by various spindle-associated proteins (SAPs) that play distinct functions in regulating spindle dynamics and microtubule-kinetochore attachment. The protozoan parasite Trypanosoma brucei employs evolutionarily conserved and kinetoplastid-specific proteins, including some kinetoplastid-specific nucleus- and spindle-associated proteins (NuSAPs), to regulate chromosome segregation. Here, we characterized NuSAP4 and its functional interplay with diverse SAPs in promoting chromosome segregation in T. brucei. NuSAP4 associates with the spindle during mitosis and concentrates at spindle poles where it interacts with SPB1 and MAP103. Knockdown of NuSAP4 impairs chromosome segregation by disrupting bipolar spindle assembly and spindle pole protein localization. These results uncover the mechanistic role of NuSAP4 in regulating chromosome segregation by promoting bipolar spindle assembly, and highlight the unusual features of mitotic regulation by spindle-associated proteins in this early divergent microbial eukaryote.

Proterozoic microfossils continue to provide new insights into the rise of complex eukaryotic life.

Eukaryotes have evolved to dominate the biosphere today, accounting for most documented living species and the vast majority of the Earth's biomass. Consequently, understanding how these biologically complex organisms initially diversified in the Proterozoic Eon over 539 million years ago is a foundational question in evolutionary biology. Over the last 70 years, palaeontologists have sought to document the rise of eukaryotes with fossil evidence. However, the delicate and microscopic nature of their sub-cellular features affords early eukaryotes diminished preservation potential. Chemical biomarker signatures of eukaryotes and the genetics of living eukaryotes have emerged as complementary tools for reconstructing eukaryote ancestry. In this review, we argue that exceptionally preserved Proterozoic microfossils are critical to interpreting these complementary tools, providing crucial calibrations to molecular clocks and testing hypotheses of palaeoecology. We highlight recent research on their preservation and biomolecular composition that offers new ways to enhance their utility.

Early Diversification of Membrane Intrinsic Proteins (MIPs) in Eukaryotes.

Membrane intrinsic proteins (MIPs), including aquaporins (AQPs) and aquaglyceroporins (GLPs), form an ancient family of transporters for water and small solutes across biological membranes. The evolutionary history and functions of MIPs have been extensively studied in vertebrates and land plants, but their widespread presence across the eukaryotic tree of life suggests both a more complex evolutionary history and a broader set of functions than previously thought. That said, the early evolution of MIPs remains obscure. The presence of one GLP and four AQP clades across both bacteria and archaea suggests that the first eukaryotes could have possessed up to five MIPs. Here, we report on a previously unknown richness in MIP diversity across all major eukaryotic lineages, including unicellular eukaryotes, which make up the bulk of eukaryotic diversity. Three MIP clades have likely deep evolutionary origins, dating back to the last eukaryotic common ancestor (LECA), and support the presence of a complex MIP repertoire in early eukaryotes. Overall, our findings highlight the growing complexity of the reconstructed LECA genome: the dynamic evolutionary history of MIPs was set in motion when eukaryotes were in their infancy followed by radiative bursts across all main eukaryotic lineages.

NADPH oxidase 5: Where are we now and which way to proceed?

Since the incorporation of mitochondria in early eukaryotes cells struggle to keep the deleterious effects of reactive oxygen species (ROS), mainly originating from the respiratory chain, at bay. Evolutionary adaptation to ROS burden went so far that by acting as messenger and effector molecules, ROS became important in maintaining homeostasis. The evolutionary success of this phenomenon is underscored by the arising of professional ROS-generating enzymes, namely the family of NADPH oxidases (NOXes). NOXes, by shaping ROS levels at different subcellular locations and in extracellular space, are involved in such fundamental functions as proliferation, differentiation, apoptosis, host defense, fertilization, and hormone biosynthesis. NOX5, being a calcium-regulated professional ROS source exerts its function at the crossroad of these two fundamental but potentially deleterious intracellular signaling pathways (i.e. Ca2+ and ROS). The expression of NOX5 in the adult human body under unchallenged conditions is restricted to very few sites, among which the two major tissue groups are genital organs (mainly testis) and immune tissues (mainly spleen). In cases of increased cellular proliferation and protein synthesis (e.g., diverse tumors, cultured primary cells, or sites of tissue damage) the expression and activity of NOX5 is often upregulated in various tissues. This and the evolutionary conserved nature of NOX5 would imply a very fundamental role for this enzyme, but intriguingly the genomes of rodents essentially lack the NOX5 gene. The latter fact had been a major obstacle in determining the physiological roles of NOX5 in normal tissues until the very recent generation of a NOX5-deficient rabbit model.

Evolutionary plasticity and functional repurposing of the essential metabolic enzyme MoeA.

MoeA, or gephyrin in higher eukaryotes, is crucial for molybdenum cofactor biosynthesis required in redox reactions. Gephyrin is a moonlighting protein also involved in postsynaptic receptor clustering, a feature thought to be a recent evolutionary trait. We showed previously that a repurposed copy of MoeA (Glp) is involved in bacterial cell division. To investigate how MoeA acquired multifunctionality, we used phylogenetic inference and protein structure analyses to understand the diversity and evolutionary history of MoeA. Glp-expressing Bacteria have at least two copies of the gene, and our analysis suggests that Glp has lost its enzymatic role. In Archaea we identified an ancestral duplication where one of the paralogs might bind tungsten instead of molybdenum. In eukaryotes, the acquisition of the moonlighting activity of gephyrin comprised three major events: first, MoeA was obtained from Bacteria by early eukaryotes, second, MogA was fused to the N-terminus of MoeA in the ancestor of opisthokonts, and finally, it acquired the function of anchoring GlyR receptors in neurons. Our results support the functional versatility and adaptive nature of the MoeA scaffold, which has been repurposed independently both in eukaryotes and bacteria to carry out analogous functions in network organization at the cell membrane.

Impact of steroid biosynthesis on the aerobic adaptation ofeukaryotes.

Steroids are indispensable components of the eukaryotic cellular membrane and the acquisition of steroid biosynthesis was a key factor that enabled the evolution of eukaryotes. The polycyclic carbon structures of steroids can be preserved in sedimentary rocks as chemical fossils for billions of years and thus provide invaluable clues to trace eukaryotic evolution from the distant past. Steroid biosynthesis consists of (1) the production of protosteroids and (2) the subsequent modifications toward "modern-type" steroids such as cholesterol and stigmasterol. While protosteroid biosynthesis requires only two genes for the cyclization of squalene, complete modification of protosteroids involves ~10 additional genes. Eukaryotes universally possess at least some of those additional genes and thus produce modern-type steroids as major final products. The geological biomarker records suggest a prolonged period of solely protosteroid production in the mid-Proterozoic before the advent of modern-type steroids in the Neoproterozoic. It has been proposed that mid-Proterozoic protosteroids were produced by hypothetical stem-group eukaryotes that presumably possessed genes only for protosteroid production, even though in modern environments protosteroid production as a final product is found exclusively in bacteria. The host identity of mid-Proterozoic steroid producers is crucial for understanding the early evolution of eukaryotes. In this perspective, we discuss how geological biomarker data and genetic data complement each other and potentially provide a more coherent scenario for the evolution of steroids and associated early eukaryotes. We further discuss the potential impacts that steroids had on the evolution of aerobic metabolism in eukaryotes, which may have been an important factor for the eventual ecological dominance of eukaryotes in many modern environments.

Foraging mechanisms in excavate flagellates shed light on the functional ecology of early eukaryotes

The phagotrophic flagellates described as "typical excavates" have been hypothesized to be morphologically similar to the Last Eukaryotic Common Ancestor and understanding the functional ecology of excavates may therefore help shed light on the ecology of these early eukaryotes. Typical excavates are characterized by a posterior flagellum equipped with a vane that beats in a ventral groove. Here, we combined flow visualization and observations of prey capture in representatives of the three clades of excavates with computational fluid dynamic modeling, to understand the functional significance of this cell architecture. We record substantial differences amongst species in the orientation of the vane and the beat plane of the posterior flagellum. Clearance rate magnitudes estimated from flow visualization and modeling are both like that of other similarly sized flagellates. The interaction between a vaned flagellum beating in a confinement is modeled to produce a very efficient feeding current at low energy costs, irrespective of the beat plane and vane orientation and of all other morphological variations. Given this predicted uniformity of function, we suggest that the foraging systems of typical excavates studied here may be good proxies to understand those potentially used by our distant ancestors more than 1 billion years ago.

Installation of LYRM proteins in early eukaryotes to regulate the metabolic capacity of the emerging mitochondrion.

Core mitochondrial processes such as the electron transport chain, protein translation and the formation of Fe-S clusters (ISC) are of prokaryotic origin and were present in the bacterial ancestor of mitochondria. In animal and fungal models, a family of small Leu-Tyr-Arg motif-containing proteins (LYRMs) uniformly regulates the function of mitochondrial complexes involved in these processes. The action of LYRMs is contingent upon their binding to the acylated form of acyl carrier protein (ACP). This study demonstrates that LYRMs are structurally and evolutionarily related proteins characterized by a core triplet of α-helices. Their widespread distribution across eukaryotes suggests that 12 specialized LYRMs were likely present in the last eukaryotic common ancestor to regulate the assembly and folding of the subunits that are conserved in bacteria but that lack LYRM homologues. The secondary reduction of mitochondria to anoxic environments has rendered the function of LYRMs and their interaction with acylated ACP dispensable. Consequently, these findings strongly suggest that early eukaryotes installed LYRMs in aerobic mitochondria as orchestrated switches, essential for regulating core metabolism and ATP production.

Dynamic localization of the chromosomal passenger complex in trypanosomes is controlled by the orphan kinesins KIN-A and KIN-B.

The chromosomal passenger complex (CPC) is an important regulator of cell division, which shows dynamic subcellular localization throughout mitosis, including kinetochores and the spindle midzone. In traditional model eukaryotes such as yeasts and humans, the CPC consists of the catalytic subunit Aurora B kinase, its activator INCENP, and the localization module proteins Borealin and Survivin. Intriguingly, Aurora B and INCENP as well as their localization pattern are conserved in kinetoplastids, an evolutionarily divergent group of eukaryotes that possess unique kinetochore proteins and lack homologs of Borealin or Survivin. It is not understood how the kinetoplastid CPC assembles nor how it is targeted to its subcellular destinations during the cell cycle. Here, we identify two orphan kinesins, KIN-A and KIN-B, as bona fide CPC proteins in Trypanosoma brucei, the kinetoplastid parasite that causes African sleeping sickness. KIN-A and KIN-B form a scaffold for the assembly of the remaining CPC subunits. We show that the C-terminal unstructured tail of KIN-A interacts with the KKT8 complex at kinetochores, while its N-terminal motor domain promotes CPC translocation to spindle microtubules. Thus, the KIN-A:KIN-B complex constitutes a unique 'two-in-one' CPC localization module, which directs the CPC to kinetochores from S phase until metaphase and to the central spindle in anaphase. Our findings highlight the evolutionary diversity of CPC proteins and raise the possibility that kinesins may have served as the original transport vehicles for Aurora kinases in early eukaryotes.

Dynamic localization of the chromosomal passenger complex in trypanosomes is controlled by the orphan kinesins KIN-A and KIN-B

The chromosomal passenger complex (CPC) is an important regulator of cell division, which shows dynamic subcellular localization throughout mitosis, including kinetochores and the spindle midzone. In traditional model eukaryotes such as yeasts and humans, the CPC consists of the catalytic subunit Aurora B kinase, its activator INCENP, and the localization module proteins Borealin and Survivin. Intriguingly, Aurora B and INCENP as well as their localization pattern are conserved in kinetoplastids, an evolutionarily divergent group of eukaryotes that possess unique kinetochore proteins and lack homologs of Borealin or Survivin. It is not understood how the kinetoplastid CPC assembles nor how it is targeted to its subcellular destinations during the cell cycle. Here, we identify two orphan kinesins, KIN-A and KIN-B, as bona fide CPC proteins in Trypanosoma brucei, the kinetoplastid parasite that causes African sleeping sickness. KIN-A and KIN-B form a scaffold for the assembly of the remaining CPC subunits. We show that the C-terminal unstructured tail of KIN-A interacts with the KKT8 complex at kinetochores, while its N-terminal motor domain promotes CPC translocation to spindle microtubules. Thus, the KIN-A:KIN-B complex constitutes a unique ‘two-in-one’ CPC localization module, which directs the CPC to kinetochores from S phase until metaphase and to the central spindle in anaphase. Our findings highlight the evolutionary diversity of CPC proteins and raise the possibility that kinesins may have served as the original transport vehicles for Aurora kinases in early eukaryotes.

A kinesin-13 family kinesin in Trypanosoma brucei regulates cytokinesis and cytoskeleton morphogenesis by promoting microtubule bundling.

The early branching eukaryote Trypanosoma brucei divides uni-directionally along the longitudinal cell axis from the cell anterior toward the cell posterior, and the cleavage furrow ingresses along the cell division plane between the new and the old flagella of a dividing bi-flagellated cell. Regulation of cytokinesis in T. brucei involves actomyosin-independent machineries and trypanosome-specific signaling pathways, but the molecular mechanisms underlying cell division plane positioning remain poorly understood. Here we report a kinesin-13 family protein, KIN13-5, that functions downstream of FPRC in the cytokinesis regulatory pathway and determines cell division plane placement. KIN13-5 localizes to multiple cytoskeletal structures, interacts with FPRC, and depends on FPRC for localization to the site of cytokinesis initiation. Knockdown of KIN13-5 causes loss of microtubule bundling at both ends of the cell division plane, leading to mis-placement of the cleavage furrow and unequal cytokinesis, and at the posterior cell tip, causing the formation of a blunt posterior. In vitro biochemical assays demonstrate that KIN13-5 bundles microtubules, providing mechanistic insights into the role of KIN13-5 in cytokinesis and posterior morphogenesis. Altogether, KIN13-5 promotes microtubule bundle formation to ensure cleavage furrow placement and to maintain posterior cytoskeleton morphology in T. brucei.

Chuaria and Tawuia fossils from ∼1.0 Ga rocks in North China: Implications for a polyphyletic origin of Chuaria and a potential biological link between these two widespread Proterozoic taxa

The Proterozoic witnessed the emergence and rise of early eukaryotes and Proterozoic carbonaceous compression macrofossils provide a valuable window for our understanding of this ecological expansion. Among them, the discoidal Chuaria and tomaculate Tawuia, which often co-occur in the same bedding layer, are perhaps the most common fossil forms in Proterozoic successions with extremely long temporal and broad geographical distributions. However, their phylogenetic interpretations and the relationship between them have been debated for decades, due to their simple morphologies lacking phylogenetically diagnostic features. In this study, we used an integrated application of reflected and transmitted light microscopy, scanning electron microscopy, biometric analysis, and Raman spectroscopy to characterize well-preserved carbonaceous compression macrofossils of Chuaria and Tawuia from the Tonian Liulaobei Formation in Huainan region, northern Anhui province, North China. The fossils are abundantly preserved in a ca. 1–2 cm thick homogeneous mudstone layer and can be considered to be closer to an ecological community than fossils from multiple different layers or localities. The comprehensive analysis shows that even Chuaria vesicles from the same layer may have a polyphyletic origin and three groups of Chuaria vesicles have been identified. Among them, small Chuaria vesicles with minimum diameter < 1.2 mm can be either related to the sphaeromorphic acritarch Leiosphaeridia or interpreted as multicellular aggregates. In contrast, large Chuaria vesicles with minimum diameter ≥ 1.2 mm are more comparable to the type material of Chuaria circularis from the Chuar Group and more closely related, either phylogenetically or ontogenetically, to Tawuia vesicles. This study takes a step on the path to solving the mystery of the most common and widely distributed Proterozoic macrofossil assemblage (i.e., Chuaria-Tawuia assemblage), and urges thorough and systematic evaluations to be taken on the possibility that Chuaria vesicles from the same Chuaria-Tawuia assemblages or sole Chuaria assemblages in other localities might also be polyphyletic.

Eukaryotic CD-NTase, STING, and viperin proteins evolved via domain shuffling, horizontal transfer, and ancient inheritance from prokaryotes.

Animals use a variety of cell-autonomous innate immune proteins to detect viral infections and prevent replication. Recent studies have discovered that a subset of mammalian antiviral proteins have homology to antiphage defense proteins in bacteria, implying that there are aspects of innate immunity that are shared across the Tree of Life. While the majority of these studies have focused on characterizing the diversity and biochemical functions of the bacterial proteins, the evolutionary relationships between animal and bacterial proteins are less clear. This ambiguity is partly due to the long evolutionary distances separating animal and bacterial proteins, which obscures their relationships. Here, we tackle this problem for 3 innate immune families (CD-NTases [including cGAS], STINGs, and viperins) by deeply sampling protein diversity across eukaryotes. We find that viperins and OAS family CD-NTases are ancient immune proteins, likely inherited since the earliest eukaryotes first arose. In contrast, we find other immune proteins that were acquired via at least 4 independent events of horizontal gene transfer (HGT) from bacteria. Two of these events allowed algae to acquire new bacterial viperins, while 2 more HGT events gave rise to distinct superfamilies of eukaryotic CD-NTases: the cGLR superfamily (containing cGAS) that has since diversified via a series of animal-specific duplications and a previously undefined eSMODS superfamily, which more closely resembles bacterial CD-NTases. Finally, we found that cGAS and STING proteins have substantially different histories, with STING protein domains undergoing convergent domain shuffling in bacteria and eukaryotes. Overall, our findings paint a picture of eukaryotic innate immunity as highly dynamic, where eukaryotes build upon their ancient antiviral repertoires through the reuse of protein domains and by repeatedly sampling a rich reservoir of bacterial antiphage genes.

RHO of plant signaling was established early in streptophyte evolution

The algal ancestors of land plants underwent a transition from a unicellular to a multicellular body plan.1 This transition likely took place early in streptophyte evolution, sometime after the divergence of the Chlorokybophyceae/Mesostigmatophyceae lineage, but before the divergence of the Klebsormidiophyceae lineage.2 How this transition was brought about is unknown; however, it was likely facilitated by the evolution of novel mechanisms to spatially regulate morphogenesis. In land plants, RHO of plant (ROP) signaling plays a conserved role in regulating polarized cell growth and cell division orientation to orchestrate morphogenesis.3,4,5,6,7,8 ROP constitutes a plant-specific subfamily of the RHO GTPases, which are more widely conserved throughout eukaryotes.9,10 Although the RHO family originated in early eukaryotes,11,12 how and when the ROP subfamily originated had remained elusive. Here, we demonstrate that ROP signaling was established early in the streptophyte lineage, sometime after the divergence of the Chlorokybophyceae/Mesostigmatophyceae lineage, but before the divergence of the Klebsormidiophyceae lineage. This period corresponds to when the unicellular-to-multicellular transition likely took place in the streptophytes. In addition to being critical for the complex morphogenesis of extant land plants, we speculate that ROP signaling contributed to morphological evolution in early streptophytes.

Frameworks for Interpreting the Early Fossil Record of Eukaryotes.

The origin of modern eukaryotes is one of the key transitions in life's history, and also one of the least understood. Although the fossil record provides the most direct view of this process, interpreting the fossils of early eukaryotes and eukaryote-grade organisms is not straightforward. We present two end-member models for the evolution of modern (i.e., crown) eukaryotes-one in which modern eukaryotes evolved early, and another in which they evolved late-and interpret key fossils within these frameworks, including where they might fit in eukaryote phylogeny and what they may tell us about the evolution of eukaryotic cell biology and ecology. Each model has different implications for understanding the rise of complex life on Earth, including different roles of Earth surface oxygenation, and makes different predictions that future paleontological studies can test.

The DUF3715 domain has a conserved role in RNA-directed transposon silencing.

RNA-directed transposon silencing operates in the mammalian soma and germline to safeguard genomic integrity. The piRNA pathway and the HUSH complex identify active transposons through recognition of their nascent transcripts, but mechanistic understanding of how these distinct pathways evolved is lacking. TASOR is an essential component of the HUSH complex. TASOR's DUF3715 domain adopts a pseudo-PARP structure and is required for transposon silencing in a manner independent of complex assembly. TEX15, an essential piRNA pathway factor, also contains the DUF3715 domain. Here, we show that TASOR's and TEX15's DUF3715 domain share extensive structural homology. We found that the DUF3715 domain arose in early eukaryotes and that in vertebrates it is restricted to TEX15, TASOR, and TASORB orthologs. While TASOR-like proteins are found throughout metazoa, TEX15 is vertebrate-specific. The branching of TEX15 and the TASOR-like DUF3715 domain likely occurred in early metazoan evolution. Remarkably, despite this vast evolutionary distance, the DUF3715 domain from divergent TEX15 sequences can functionally substitute the DUF3715 domain of TASOR and mediates transposon silencing. We have thus termed this domain of unknown function as the RNA-directed pseudo-PARP transposon silencing (RDTS) domain. In summary, we show an unexpected functional link between these critical transposon silencing pathways.

Titanium isotope evidence for the high topography of Nuna and Gondwana - Implications for Earth's redox and biological evolution

Titanium isotopes recorded in glacial diamictites with depositional ages between 2.9 and 0.3 Ga show that the upper continental crust became significantly more felsic relative to the present-day crust during the amalgamation of the Paleoproterozoic Nuna and the Neoproterozoic Gondwana supercontinents. This can be attributed to the continental collisions involved in the assembly of Nuna and Gondwana. The resulting high topographic relief of Nuna and Gondwana orogens must have resulted in an enhanced erosional supply from the continents to oceans. The step changes in the development of organismal complexity from prokaryotes to eukaryotes, and eventually metazoans, appear to be temporally correlated to instances where collisional mountain-building sustained an elevated nutrient supply from the continents to oceans. The nutrient surge associated with the rise of the Gondwana mountains likely provided the necessary impetus for the Neoproterozoic ecological expansion of eukaryotes and the eventual radiation of metazoans. A similar link between the enhanced nutrient supply from Nuna mountains and the radiation of early eukaryotes is plausible, although its mechanistic underpinnings remain unclear. The termination of Nuna orogeny and its transition to Rodinia without significant breakup and subsequent collisional orogenesis corresponds to the long lull in Earth's redox and biological evolution in its middle age.

Identification of hidden associations among eukaryotic genes through statistical analysis of coevolutionary transitions

Coevolution at the gene level, as reflected by correlated events of gene loss or gain, can be revealed by phylogenetic profile analysis. The optimal method and metric for comparing phylogenetic profiles, especially in eukaryotic genomes, are not yet established. Here, we describe a procedure suitable for large-scale analysis, which can reveal coevolution based on the assessment of the statistical significance of correlated presence/absence transitions between gene pairs. This metric can identify coevolution in profiles with low overall similarities and is not affected by similarities lacking coevolutionary information. We applied the procedure to a large collection of 60,912 orthologous gene groups (orthogroups) in 1,264 eukaryotic genomes extracted from OrthoDB. We found significant cotransition scores for 7,825 orthogroups associated in 2,401 coevolving modules linking known and unknown genes in protein complexes and biological pathways. To demonstrate the ability of the method to predict hidden gene associations, we validated through experiments the involvement of vertebrate malate synthase-like genes in the conversion of (S)-ureidoglycolate into glyoxylate and urea, the last step of purine catabolism. This identification explains the presence of glyoxylate cycle genes in metazoa and suggests an anaplerotic role of purine degradation in early eukaryotes.