Phylogenomic analyses support the monophyly of Excavata and resolve relationships among eukaryotic "supergroups"

The phage shock protein (PSP) envelope stress response: discovery of novel partners and evolutionary history

Bacterial phage shock protein (PSP) systems stabilize the bacterial cell membrane and protect against envelope stress. These systems have been associated with virulence, but despite their critical roles, PSP components are not well characterized outside proteobacteria. Using comparative genomics and protein sequence-structure-function analyses, we systematically identified and analyzed PSP homologs, phyletic patterns, domain architectures, and gene neighborhoods. This approach underscored the evolutionary significance of the system, revealing that its core protein PspA (Snf7 in ESCRT outside bacteria) was present in the last universal common ancestor and that this ancestral functionality has since diversified into multiple novel, distinct PSP systems across life. Several novel partners of the PSP system were identified: (i) the Toastrack domain, likely facilitating assembly of sub-membrane stress-sensing and signaling complexes, (ii) the newly defined HTH-associated α-helical signaling domain-PadR-like transcriptional regulator pair system, and (iii) multiple independent associations with ATPase, CesT/Tir-like chaperone, and Band-7 domains in proteins thought to mediate sub-membrane dynamics. Our work also uncovered links between the PSP components and other domains, such as novel variants of SHOCT-like domains, suggesting roles in assembling membrane-associated complexes of proteins with disparate biochemical functions. Results are available at our interactive web app, https://jravilab.org/psp.IMPORTANCEPhage shock proteins (PSP) are virulence-associated, cell membrane stress-protective systems. They have mostly been characterized in Proteobacteria and Firmicutes. We now show that a minimal PSP system was present in the last universal common ancestor that evolved and diversified into newly identified functional contexts. Recognizing the conservation and evolution of PSP systems across bacterial phyla contributes to our understanding of stress response mechanisms in prokaryotes. Moreover, the newly discovered PSP modularity will likely prompt new studies of lineage-specific cell envelope structures, lifestyles, and adaptation mechanisms. Finally, our results validate the use of domain architecture and genetic context for discovery in comparative genomics.

Biogenesis of extracellular vesicles from the pathogen perspective: Transkingdom strategies for delivering messages

EVs are nanoparticles enclosing proteins, nucleic acids and lipids released by cells and are essential for their metabolism and useful for intercellular communication. The importance of EVs has been highlighted by their use as biomarkers or as vaccine antigens. The release of vesicles is exploited by a wide range of organisms: from unicellular bacteria or protozoa to multicellular prokaryotes like fungi, helminths and arthropods. The mechanisms elucidated to date in each biological group are presented, as well as a discussion of interesting directions for future EV studies.

Ancient eukaryotic protein interactions illuminate modern genetic traits and disorders

All eukaryotes share a common ancestor from roughly 1.5 - 1.8 billion years ago, a single-celled, swimming microbe known as LECA, the Last Eukaryotic Common Ancestor. Nearly half of the genes in modern eukaryotes were present in LECA, and many current genetic diseases and traits stem from these ancient molecular systems. To better understand these systems, we compared genes across modern organisms and identified a core set of 10,092 shared protein-coding gene families likely present in LECA, a quarter of which are uncharacterized. We then integrated >26,000 mass spectrometry proteomics analyses from 31 species to infer how these proteins interact in higher-order complexes. The resulting interactome describes the biochemical organization of LECA, revealing both known and new assemblies. We analyzed these ancient protein interactions to find new human gene-disease relationships for bone density and congenital birth defects, demonstrating the value of ancestral protein interactions for guiding functional genetics today.

Lighting lantern above Psalteriomonadidae: Unveiling novel diversity within the genus Psalteriomonas (Discoba: Heterolobosea)

Psalteriomonadidae are a small family of anaerobic free-living protists belonging to Heterolobosea, Discoba. We cultured 74 new strains of mostly amoeboid Psalteriomonadidae obtained from mainly freshwater habitats and sequenced their 18S rRNA gene. Based on the phylogenetic analysis and genetic distances, we report multiple novel species, four of which we formally describe based on the light-microscopic morphology (Psalteriomonas minuta, P. australis, P. fimbriata, and P. parva). We also examined the ultrastructure of two Psalteriomonas species using transmission electron microscopy. We transfer Sawyeria marylandensis into the genus Psalteriomonas and synonymize Sawyeria with Psalteriomonas. In addition, we studied the flagellate stage of P. marylandensis comb. nov. for the first time, using light and scanning electron microscopy.

Euglena's atypical respiratory chain adapts to the discoidal cristae and flexible metabolism

Euglena gracilis, a model organism of the eukaryotic supergroup Discoba harbouring also clinically important parasitic species, possesses diverse metabolic strategies and an atypical electron transport chain. While structures of the electron transport chain complexes and supercomplexes of most other eukaryotic clades have been reported, no similar structure is currently available for Discoba, limiting the understandings of its core metabolism and leaving a gap in the evolutionary tree of eukaryotic bioenergetics. Here, we report high-resolution cryo-EM structures of Euglena's respirasome I + III2 + IV and supercomplex III2 + IV2. A previously unreported fatty acid synthesis domain locates on the tip of complex I's peripheral arm, providing a clear picture of its atypical subunit composition identified previously. Individual complexes are re-arranged in the respirasome to adapt to the non-uniform membrane curvature of the discoidal cristae. Furthermore, Euglena's conformationally rigid complex I is deactivated by restricting ubiquinone's access to its substrate tunnel. Our findings provide structural insights for therapeutic developments against euglenozoan parasite infections.

16S and 18S rDNA Amplicon Sequencing Analysis of Aesthetically Problematic Microbial Mats on the Walls of the Petralona Cave: The Use of Essential Oils as a Cleaning Method

The presence of microbial communities on cave walls and speleothems is an issue that requires attention. Traditional cleaning methods using water, brushes, and steam can spread the infection and cause damage to the cave structures, while chemical agents can lead to the formation of toxic compounds and damage the cave walls. Essential oils (EOs) have shown promising results in disrupting the cell membrane of bacteria and affecting their membrane permeability. In this study, we identified the microorganisms forming unwanted microbial communities on the walls and speleothems of Petralona Cave using 16S and 18S rDNA amplicon sequencing approaches and evaluated the efficacy of EOs in reducing the ATP levels of these ecosystems. The samples exhibited a variety of both prokaryotic and eukaryotic microorganisms, including Proteobacteria, Actinobacteria, Bacteroidetes, Chloroflexi, Firmicutes, the SAR supergroup, Opisthokonta, Excavata, Archaeplastida, and Amoebozoa. These phyla are often found in various habitats, including caves, and contribute to the ecological intricacy of cave ecosystems. In terms of the order and genus taxonomy, the identified biota showed abundances that varied significantly among the samples. Functional predictions were also conducted to estimate the differences in expressed genes among the samples. Oregano EO was found to reduce ATP levels by 87% and 46% for black and green spots, respectively. Consecutive spraying with cinnamon EO further reduced ATP levels, with reductions of 89% for black and 88% for green spots. The application of a mixture solution caused a significant reduction up to 96% in ATP levels of both areas. Our results indicate that EOs could be a promising solution for the treatment of microbial communities on cave walls and speleothems.

The free-living flagellate Paratrimastix pyriformis uses a distinct mitochondrial carrier to balance adenine nucleotide pools

Paratrimastix pyriformis is a free-living flagellate thriving in low-oxygen freshwater sediments. It belongs to the group Metamonada along with human parasites, such as Giardia and Trichomonas. Like other metamonads, P. pyriformis has a mitochondrion-related organelle (MRO) which in this protist is primarily involved in one-carbon folate metabolism. The MRO contains four members of the solute carrier family 25 (SLC25) responsible for the exchange of metabolites across the mitochondrial inner membrane. Here, we characterise the function of the adenine nucleotide carrier PpMC1 by thermostability shift and transport assays. We show that it transports ATP, ADP and, to a lesser extent, AMP, but not phosphate. The carrier is distinct in function and origin from both ADP/ATP carriers and ATP-Mg/phosphate carriers, and it most likely represents a distinct class of adenine nucleotide carriers.

Evolutionary diversification of the autophagy-related ubiquitin-like conjugation systems

Two autophagy-related (ATG) ubiquitin-like conjugation systems, the ATG12 and ATG8 systems, play important roles in macroautophagy. While multiple duplications and losses of the ATG conjugation system proteins are found in different lineages, the extent to which the underlying systems diversified across eukaryotes is not fully understood. Here, in order to understand the evolution of the ATG conjugation systems, we constructed a transcriptome database consisting of 94 eukaryotic species covering major eukaryotic clades and systematically identified ATG conjugation system components. Both ATG10 and the C-terminal glycine of ATG12 are essential for the canonical ubiquitin-like conjugation of ATG12 and ATG5. However, loss of ATG10 or the C-terminal glycine of ATG12 occurred at least 16 times in a wide range of lineages, suggesting that possible covalent-to-non-covalent transition is not limited to the species that we previously reported such as Alveolata and some yeast species. Some species have only the ATG8 system (with conjugation enzymes) or only ATG8 (without conjugation enzymes). More than 10 species have ATG8 homologs without the conserved C-terminal glycine, and Tetrahymena has an ATG8 homolog with a predicted transmembrane domain, which may be able to anchor to the membrane independent of the ATG conjugation systems. We discuss the possibility that the ancestor of the ATG12 and ATG8 systems is more similar to ATG8. Overall, our study offers a whole picture of the evolution and diversity of the ATG conjugation systems among eukaryotes, and provides evidence that functional diversifications of the systems are more common than previously thought.Abbreviations: APEAR: ATG8-PE association region; ATG: autophagy-related; LIR: LC3-interacting region; NEDD8: neural precursor cell expressed, developmentally down-regulated gene 8; PE: phosphatidylethanolamine; SAMP: small archaeal modifier protein; SAR: Stramenopiles, Alveolata, and Rhizaria; SMC: structural maintenance of chromosomes; SUMO: small ubiquitin like modifier; TACK: Thaumarchaeota, Aigarchaeota, Crenarchaeota, and Korarchaeota; UBA: ubiquitin like modifier activating enzyme; UFM: ubiquitin fold modifier; URM: ubiquitin related modifier.

Biogeography and Population Divergence of Microeukaryotes Associated with Fluids and Chimneys in the Hydrothermal Vents of the Southwest Indian Ocean

Deep-sea hydrothermal vents have been proposed as oases for microbes, but microeukaryotes as key components of the microbial loop have not been well studied. Based on high-throughput sequencing and network analysis of the 18S rRNA gene, distinct biogeographical distribution patterns and impacting factors were revealed from samples in the three hydrothermal fields of the southwest Indian Ocean, where higher gene abundance of microeukaryotes appeared in chimneys. The microeukaryotes in the fluids might be explained by hydrogeochemical heterogeneity, especially that of the nitrate and silicate concentrations, while the microeukaryotes in the chimneys coated with either Fe oxides or Fe-Si oxyhydroxides might be explained by potentially different associated prokaryotic groups. Population divergence of microeukaryotes, especially clades of parasitic Syndiniales, was observed among different hydrothermal fluids and chimneys and deserves further exploration to gain a deeper understanding of the trophic relationships and potential ecological function of microeukaryotes in the deep-sea extreme ecosystems, especially in the complex deep-sea chemoautotrophic habitats. IMPORTANCE Deep-sea hydrothermal vents have been proposed as oases for microbes, but microeukaryotes as key components of the microbial loop have not been well studied. Based on high-throughput sequencing and network analysis of the 18S rRNA gene, population divergence of microeukaryotes, especially clades of parasitic Syndiniales, was observed among different hydrothermal fields. This might be attributed to the hydrogeochemical heterogeneity of fluids and to the potentially different associated prokaryotic groups in chimneys.

Free-Living Amoebas in Extreme Environments: The True Survival in our Planet

Free-living amoebas (FLAs) are microorganisms, unicellular protozoa widely distributed in nature and present in different environments, such as water or soil; they are maintained in ecosystems and play a fundamental role in the biological control of bacteria, other protozoa, and mushrooms. In particular circumstances, some can reach humans or animals, promoting several health complications. Notably, FLAs are characterized by a robust capacity to survive in extreme environments. However, currently, there is no updated information on the existence and distribution of this protozoan in inhospitable places. Undoubtedly, the cellular physiology of these protozoan microorganisms is very particular. They can resist and live in extreme environments due to their encysting capacity and tolerance to different osmolarities, temperatures, and other environmental factors, which give them excellent adaptative resistance. In this review, we summarized the most relevant evidence related to FLAs and the possible mechanism, which could explain their adaptative capacity to several extreme environments.

Sex in protists: A new perspective on the reproduction mechanisms of trypanosomatids

The Protist kingdom individuals are the most ancestral representatives of eukaryotes. They have inhabited Earth since ancient times and are currently found in the most diverse environments presenting a great heterogeneity of life forms. The unicellular and multicellular algae, photosynthetic and heterotrophic organisms, as well as free-living and pathogenic protozoa represents the protist group. The evolution of sex is directly associated with the origin of eukaryotes being protists the earliest protagonists of sexual reproduction on earth. In eukaryotes, the recombination through genetic exchange is a ubiquitous mechanism that can be stimulated by DNA damage. Scientific evidences support the hypothesis that reactive oxygen species (ROS) induced DNA damage can promote sexual recombination in eukaryotes which might have been a decisive factor for the origin of sex. The fact that some recombination enzymes also participate in meiotic sex in modern eukaryotes reinforces the idea that sexual reproduction emerged as consequence of specific mechanisms to cope with mutations and alterations in genetic material. In this review we will discuss about origin of sex and different strategies of evolve sexual reproduction in some protists such that cause human diseases like malaria, toxoplasmosis, sleeping sickness, Chagas disease, and leishmaniasis.

The Porifera microeukaryome: Addressing the neglected associations between sponges and protists

While bacterial and archaeal communities of sponges are intensively studied, given their importance to the animal's physiology as well as sources of several new bioactive molecules, the potential and roles of associated protists remain poorly known. Historically, culture-dependent approaches dominated the investigations of sponge-protist interactions. With the advances in omics techniques, these associations could be visualized at other equally important scales. Of the few existing studies, there is a strong tendency to focus on interactions with photosynthesizing taxa such as dinoflagellates and diatoms, with fewer works dissecting the interactions with other less common groups. In addition, there are bottlenecks and inherent biases in using primer pairs and bioinformatics approaches in the most commonly used metabarcoding studies. Thus, this review addresses the issues underlying this association, using the term "microeukaryome" to refer exclusively to protists associated with an animal host. We aim to highlight the diversity and community composition of protists associated with sponges and place them on the same level as other microorganisms already well studied in this context. Among other shortcomings, it could be observed that the biotechnological potential of the microeukaryome is still largely unexplored, possibly being a valuable source of new pharmacological compounds, enzymes and metabolic processes.

Combined nanometric and phylogenetic analysis of unique endocytic compartments in Giardia lamblia sheds light on the evolution of endocytosis in Metamonada

Background

Giardia lamblia, a parasitic protist of the Metamonada supergroup, has evolved one of the most diverged endocytic compartment systems investigated so far. Peripheral endocytic compartments, currently known as peripheral vesicles or vacuoles (PVs), perform bulk uptake of fluid phase material which is then digested and sorted either to the cell cytosol or back to the extracellular space.

Results

Here, we present a quantitative morphological characterization of these organelles using volumetric electron microscopy and super-resolution microscopy (SRM). We defined a morphological classification for the heterogenous population of PVs and performed a comparative analysis of PVs and endosome-like organelles in representatives of phylogenetically related taxa, Spironucleus spp. and Tritrichomonas foetus. To investigate the as-yet insufficiently understood connection between PVs and clathrin assemblies in G. lamblia, we further performed an in-depth search for two key elements of the endocytic machinery, clathrin heavy chain (CHC) and clathrin light chain (CLC), across different lineages in Metamonada. Our data point to the loss of a bona fide CLC in the last Fornicata common ancestor (LFCA) with the emergence of a protein analogous to CLC (GlACLC) in the Giardia genus. Finally, the location of clathrin in the various compartments was quantified.

Conclusions

Taken together, this provides the first comprehensive nanometric view of Giardia’s endocytic system architecture and sheds light on the evolution of GlACLC analogues in the Fornicata supergroup and, specific to Giardia, as a possible adaptation to the formation and maintenance of stable clathrin assemblies at PVs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01402-3.

Phylogenomic Analyses of 2,786 Genes in 158 Lineages Support a Root of the Eukaryotic Tree of Life between Opisthokonts and All Other Lineages

Advances in phylogenomics and high-throughput sequencing have allowed the reconstruction of deep phylogenetic relationships in the evolution of eukaryotes. Yet, the root of the eukaryotic tree of life remains elusive. The most popular hypothesis in textbooks and reviews is a root between Unikonta (Opisthokonta + Amoebozoa) and Bikonta (all other eukaryotes), which emerged from analyses of a single-gene fusion. Subsequent, highly cited studies based on concatenation of genes supported this hypothesis with some variations or proposed a root within Excavata. However, concatenation of genes does not consider phylogenetically-informative events like gene duplications and losses. A recent study using gene tree parsimony (GTP) suggested the root lies between Opisthokonta and all other eukaryotes, but only including 59 taxa and 20 genes. Here we use GTP with a duplication-loss model in a gene-rich and taxon-rich dataset (i.e., 2,786 gene families from two sets of 155 and 158 diverse eukaryotic lineages) to assess the root, and we iterate each analysis 100 times to quantify tree space uncertainty. We also contrasted our results and discarded alternative hypotheses from the literature using GTP and the likelihood-based method SpeciesRax. Our estimates suggest a root between Fungi or Opisthokonta and all other eukaryotes; but based on further analysis of genome size, we propose that the root between Opisthokonta and all other eukaryotes is the most likely.

A natural fusion of flavodiiron, rubredoxin, and rubredoxin oxidoreductase domains is a self-sufficient water-forming oxidase of Trichomonas vaginalis

Microaerophilic pathogens such as Giardia lamblia, Entamoeba histolytica, and Trichomonas vaginalis have robust oxygen consumption systems to detoxify oxygen and maintain intracellular redox balance. This oxygen consumption results from H2O-forming NADH oxidase (NOX) activity of two distinct flavin-containing systems: H2O-forming NOXes and multicomponent flavodiiron proteins (FDPs). Neither system is membrane bound, and both recycle NADH into oxidized NAD+ while simultaneously removing O2 from the local environment. However, little is known about the specific contributions of these systems in T. vaginalis. In this study, we use bioinformatics and biochemical analyses to show that T. vaginalis lacks a NOX-like enzyme and instead harbors three paralogous genes (FDPF1-3), each encoding a natural fusion product between the N-terminal FDP, central rubredoxin (Rb), and C-terminal NADH:Rb oxidoreductase domains. Unlike a "stand-alone" FDP that lacks Rb and oxidoreductase domains, this natural fusion protein with fully populated flavin redox centers directly accepts reducing equivalents of NADH to catalyze the four-electron reduction of oxygen to water within a single polypeptide with an extremely high turnover. Furthermore, using single-particle cryo-EM, we present structural insights into the spatial organization of the FDP core within this multidomain fusion protein. Together, these results contribute to our understanding of systems that allow protozoan parasites to maintain optimal redox balance and survive transient exposure to oxic conditions.

Ciliary transition zone evolution and the root of the eukaryote tree: implications for opisthokont origin and classification of kingdoms Protozoa, Plantae, and Fungi

I thoroughly discuss ciliary transition zone (TZ) evolution, highlighting many overlooked evolutionarily significant ultrastructural details. I establish fundamental principles of TZ ultrastructure and evolution throughout eukaryotes, inferring unrecognised ancestral TZ patterns for Fungi, opisthokonts, and Corticata (i.e., kingdoms Plantae and Chromista). Typical TZs have a dense transitional plate (TP), with a previously overlooked complex lattice as skeleton. I show most eukaryotes have centriole/TZ junction acorn-V filaments (whose ancestral function was arguably supporting central pair microtubule-nucleating sites; I discuss their role in centriole growth). Uniquely simple malawimonad TZs (without TP, simpler acorn) pinpoint the eukaryote tree's root between them and TP-bearers, highlighting novel superclades. I integrate TZ/ciliary evolution with the best multiprotein trees, naming newly recognised major eukaryote clades and revise megaclassification of basal kingdom Protozoa. Recent discovery of non-photosynthetic phagotrophic flagellates with genome-free plastids (Rhodelphis), the sister group to phylum Rhodophyta (red algae), illuminates plant and chromist early evolution. I show previously overlooked marked similarities in cell ultrastructure between Rhodelphis and Picomonas, formerly considered an early diverging chromist. In both a nonagonal tube lies between their TP and an annular septum surrounding their 9+2 ciliary axoneme. Mitochondrial dense condensations and mitochondrion-linked smooth endomembrane cytoplasmic partitioning cisternae further support grouping Picomonadea and Rhodelphea as new plant phylum Pararhoda. As Pararhoda/Rhodophyta form a robust clade on site-heterogeneous multiprotein trees, I group Pararhoda and Rhodophyta as new infrakingdom Rhodaria of Plantae within subkingdom Biliphyta, which also includes Glaucophyta with fundamentally similar TZ, uniquely in eukaryotes. I explain how biliphyte TZs generated viridiplant stellate-structures.

The Giardia ventrolateral flange is a lamellar membrane protrusion that supports attachment

Attachment to the intestinal epithelium is critical to the lifestyle of the ubiquitous parasite Giardia lamblia. The ventrolateral flange is a sheet-like membrane protrusion at the interface between parasites and attached surfaces. This structure has been implicated in attachment, but its role has been poorly defined. Here, we identified a novel actin associated protein with putative WH2-like actin binding domains we named Flangin. Flangin complexes with Giardia actin (GlActin) and is enriched in the ventrolateral flange making it a valuable marker for studying the flanges' role in Giardia biology. Live imaging revealed that the flange grows to around 1 μm in width after cytokinesis, then remains uniform in size during interphase, grows in mitosis, and is resorbed during cytokinesis. A flangin truncation mutant stabilizes the flange and blocks cytokinesis, indicating that flange disassembly is necessary for rapid myosin-independent cytokinesis in Giardia. Rho family GTPases are important regulators of membrane protrusions and GlRac, the sole Rho family GTPase in Giardia, was localized to the flange. Knockdown of Flangin, GlActin, and GlRac result in flange formation defects. This indicates a conserved role for GlRac and GlActin in forming membrane protrusions, despite the absence of canonical actin binding proteins that link Rho GTPase signaling to lamellipodia formation. Flangin-depleted parasites had reduced surface contact and when challenged with fluid shear force in flow chambers they had a reduced ability to remain attached, confirming a role for the flange in attachment. This secondary attachment mechanism complements the microtubule based adhesive ventral disc, a feature that may be particularly important during mitosis when the parental ventral disc disassembles in preparation for cytokinesis. This work supports the emerging view that Giardia's unconventional actin cytoskeleton has an important role in supporting parasite attachment.

SpeciesRax: A tool for maximum likelihood species tree inference from gene family trees under duplication, transfer, and loss

Species tree inference from gene family trees is becoming increasingly popular because it can account for discordance between the species tree and the corresponding gene family trees. In particular, methods that can account for multiple-copy gene families exhibit potential to leverage paralogy as informative signal. At present, there does not exist any widely adopted inference method for this purpose. Here, we present SpeciesRax, the first maximum likelihood method that can infer a rooted species tree from a set of gene family trees and can account for gene duplication, loss, and transfer events. By explicitly modeling events by which gene trees can depart from the species tree, SpeciesRax leverages the phylogenetic rooting signal in gene trees. SpeciesRax infers species tree branch lengths in units of expected substitutions per site and branch support values via paralogy-aware quartets extracted from the gene family trees. Using both empirical and simulated data sets we show that SpeciesRax is at least as accurate as the best competing methods while being one order of magnitude faster on large data sets at the same time. We used SpeciesRax to infer a biologically plausible rooted phylogeny of the vertebrates comprising 188 species from 31,612 gene families in 1 h using 40 cores. SpeciesRax is available under GNU GPL at https://github.com/BenoitMorel/GeneRax and on BioConda.

Anaeramoebae are a divergent lineage of eukaryotes that shed light on the transition from anaerobic mitochondria to hydrogenosomes

Discoveries of diverse microbial eukaryotes and their inclusion in comprehensive phylogenomic analyses have crucially re-shaped the eukaryotic tree of life in the 21st century.¹ At the deepest level, eukaryotic diversity comprises 9-10 "supergroups." One of these supergroups, the Metamonada, is particularly important to our understanding of the evolutionary dynamics of eukaryotic cells, including the remodeling of mitochondrial function. All metamonads thrive in low-oxygen environments and lack classical aerobic mitochondria, instead possessing mitochondrion-related organelles (MROs) with metabolisms that are adapted to low-oxygen conditions. These MROs lack an organellar genome, do not participate in the Krebs cycle and oxidative phosphorylation,² and often synthesize ATP by substrate-level phosphorylation coupled to hydrogen production.^3,4 The events that occurred during the transition from an oxygen-respiring mitochondrion to a functionally streamlined MRO early in metamonad evolution remain largely unknown. Here, we report transcriptomes of two recently described, enigmatic, anaerobic protists from the genus Anaeramoeba.⁵ Using phylogenomic analysis, we show that these species represent a divergent, phylum-level lineage in the tree of metamonads, emerging as a sister group of the Parabasalia and reordering the deep branching order of the metamonad tree. Metabolic reconstructions of the Anaeramoeba MROs reveal many "classical" mitochondrial features previously not seen in metamonads, including a disulfide relay import system, propionate production, and amino acid metabolism. Our findings suggest that the cenancestor of Metamonada likely had MROs with more classical mitochondrial features than previously anticipated and demonstrate how discoveries of novel lineages of high taxonomic rank continue to transform our understanding of early eukaryote evolution.

Genomic analysis finds no evidence of canonical eukaryotic DNA processing complexes in a free-living protist

Cells replicate and segregate their DNA with precision. Previous studies showed that these regulated cell-cycle processes were present in the last eukaryotic common ancestor and that their core molecular parts are conserved across eukaryotes. However, some metamonad parasites have secondarily lost components of the DNA processing and segregation apparatuses. To clarify the evolutionary history of these systems in these unusual eukaryotes, we generated a genome assembly for the free-living metamonad Carpediemonas membranifera and carried out a comparative genomics analysis. Here, we show that parasitic and free-living metamonads harbor an incomplete set of proteins for processing and segregating DNA. Unexpectedly, Carpediemonas species are further streamlined, lacking the origin recognition complex, Cdc6 and most structural kinetochore subunits. Carpediemonas species are thus the first known eukaryotes that appear to lack this suite of conserved complexes, suggesting that they likely rely on yet-to-be-discovered or alternative mechanisms to carry out these fundamental processes.

An emerging paradigm for the origin and evolution of shelled amoebae, integrating advances from molecular phylogenetics, morphology and paleontology

The phylogenetic paradigm of eukaryotic evolution has changed dramatically over the past two decades, with profound reflections on the understanding of life on earth. Arcellinida testate (shelled) amoebae lineages represent some of the oldest fossils of eukaryotes, and the elucidation of their phylogenetic relationships opened a window to the distant past, with important implications for understanding the evolution of life on earth. This four-part essay summarises advances made in the past 20 years regarding: (i) the phylogenetic relationships among amoebae with shells evolving in concert with the advances made in the phylogeny of eukaryotes; (ii) paleobiological studies unraveling the biological affinities of Neoproterozoic vase-shaped microfossils (VSMs); (iii) the interwoven interpretation of these different sets of data concluding that the Neoproterozoic contains a surprising diversity of organisms, in turn demanding a reinterpretation of the most profound events we know in the history of eukaryotes, and; (iv) a synthesis of the current knowledge about the evolution of Arcellinida, together with the possibilities and pitfalls of their interpretation.

Mitochondrial Genomic Landscape: A Portrait of the Mitochondrial Genome 40 Years after the First Complete Sequence

Notwithstanding the initial claims of general conservation, mitochondrial genomes are a largely heterogeneous set of organellar chromosomes which displays a bewildering diversity in terms of structure, architecture, gene content, and functionality. The mitochondrial genome is typically described as a single chromosome, yet many examples of multipartite genomes have been found (for example, among sponges and diplonemeans); the mitochondrial genome is typically depicted as circular, yet many linear genomes are known (for example, among jellyfish, alveolates, and apicomplexans); the chromosome is normally said to be “small”, yet there is a huge variation between the smallest and the largest known genomes (found, for example, in ctenophores and vascular plants, respectively); even the gene content is highly unconserved, ranging from the 13 oxidative phosphorylation-related enzymatic subunits encoded by animal mitochondria to the wider set of mitochondrial genes found in jakobids. In the present paper, we compile and describe a large database of 27,873 mitochondrial genomes currently available in GenBank, encompassing the whole eukaryotic domain. We discuss the major features of mitochondrial molecular diversity, with special reference to nucleotide composition and compositional biases; moreover, the database is made publicly available for future analyses on the MoZoo Lab GitHub page.

Prokaryotic and eukaryotic diversity in hydrothermal continental systems

The term extremophile was suggested more than 30 years ago and represents microorganisms that are capable of developing and living under extreme conditions, these conditions being particularly hostile to other types of microorganisms and to humankind. In terrestrial hydrothermal sites, like hot springs, "mud pools", solfataras, and geysers, the dominant extreme conditions are high temperature, low or high pH, and high levels of salinity. The diversity of microorganisms inhabiting these sites is determined by the conditions of the environment. Organisms belonging to the domains Archaea and Bacteria are more represented than the one belonging to Eukarya. Eukarya members tend to be less present because of their lower tolerance to higher temperatures, however, they perform important ecosystem processes when present. Both prokaryotes and eukaryotes have morphological and physical adaptations that allow them to colonize extreme environments. Microbial mats are complex associations of microorganisms that help the colonization of more extreme systems. In this review, a characterization of prokaryotic and eukaryotic organisms that populate terrestrial hydrothermal systems are made.

Phylogenomic fingerprinting of tempo and functions of horizontal gene transfer within ochrophytes

Horizontal gene transfer (HGT) is an important source of novelty in eukaryotic genomes. This is particularly true for the ochrophytes, a diverse and important group of algae. Previous studies have shown that ochrophytes possess a mosaic of genes derived from bacteria and eukaryotic algae, acquired through chloroplast endosymbiosis and from HGTs, although understanding of the time points and mechanisms underpinning these transfers has been restricted by the depth of taxonomic sampling possible. We harness an expanded set of ochrophyte sequence libraries, alongside automated and manual phylogenetic annotation, in silico modeling, and experimental techniques, to assess the frequency and functions of HGT across this lineage. Through manual annotation of thousands of single-gene trees, we identify continuous bacterial HGT as the predominant source of recently arrived genes in the model diatom Phaeodactylum tricornutum Using a large-scale automated dataset, a multigene ochrophyte reference tree, and mathematical reconciliation of gene trees, we note a probable elevation of bacterial HGTs at foundational points in diatom evolution, following their divergence from other ochrophytes. Finally, we demonstrate that throughout ochrophyte evolutionary history, bacterial HGTs have been enriched in genes encoding secreted proteins. Our study provides insights into the sources and frequency of HGTs, and functional contributions that HGT has made to algal evolution.

Sequential production of gametes during meiosis in trypanosomes

Meiosis is a core feature of eukaryotes that occurs in all major groups, including the early diverging excavates. In this group, meiosis and production of haploid gametes have been described in the pathogenic protist, Trypanosoma brucei, and mating occurs in the salivary glands of the insect vector, the tsetse fly. Here, we searched for intermediate meiotic stages among trypanosomes from tsetse salivary glands. Many different cell types were recovered, including trypanosomes in Meiosis I and gametes. Significantly, we found trypanosomes containing three nuclei with a 1:2:1 ratio of DNA contents. Some of these cells were undergoing cytokinesis, yielding a mononucleate gamete and a binucleate cell with a nuclear DNA content ratio of 1:2. This cell subsequently produced three more gametes in two further rounds of division. Expression of the cell fusion protein HAP2 (GCS1) was not confined to gametes, but also extended to meiotic intermediates. We propose a model whereby the two nuclei resulting from Meiosis I undergo asynchronous Meiosis II divisions with sequential production of haploid gametes.

Gene Duplications Trace Mitochondria to the Onset of Eukaryote Complexity

The last eukaryote common ancestor (LECA) possessed mitochondria and all key traits that make eukaryotic cells more complex than their prokaryotic ancestors, yet the timing of mitochondrial acquisition and the role of mitochondria in the origin of eukaryote complexity remain debated. Here, we report evidence from gene duplications in LECA indicating an early origin of mitochondria. Among 163,545 duplications in 24,571 gene trees spanning 150 sequenced eukaryotic genomes, we identify 713 gene duplication events that occurred in LECA. LECA's bacterial-derived genes include numerous mitochondrial functions and were duplicated significantly more often than archaeal-derived and eukaryote-specific genes. The surplus of bacterial-derived duplications in LECA most likely reflects the serial copying of genes from the mitochondrial endosymbiont to the archaeal host's chromosomes. Clustering, phylogenies and likelihood ratio tests for 22.4 million genes from 5,655 prokaryotic and 150 eukaryotic genomes reveal no evidence for lineage-specific gene acquisitions in eukaryotes, except from the plastid in the plant lineage. That finding, and the functions of bacterial genes duplicated in LECA, suggests that the bacterial genes in eukaryotes are acquisitions from the mitochondrion, followed by vertical gene evolution and differential loss across eukaryotic lineages, flanked by concomitant lateral gene transfer among prokaryotes. Overall, the data indicate that recurrent gene transfer via the copying of genes from a resident mitochondrial endosymbiont to archaeal host chromosomes preceded the onset of eukaryotic cellular complexity, favoring mitochondria-early over mitochondria-late hypotheses for eukaryote origin.

The sexual side of parasitic protists

Much of the vast evolutionary landscape occupied by Eukaryotes is dominated by protists. Though parasitism has arisen in many lineages, there are three main groups of parasitic protists of relevance to human and livestock health: the Apicomplexa, including the malaria parasite Plasmodium and coccidian pathogens of livestock such as Eimeria; the excavate flagellates, encompassing a diverse range of protist pathogens including trypanosomes, Leishmania, Giardia and Trichomonas; and the Amoebozoa, including pathogenic amoebae such as Entamoeba. These three groups represent separate, deep branches of the eukaryote tree, underlining their divergent evolutionary histories. Here, I explore what is known about sex in these three main groups of parasitic protists.

HERMES: An improved method to test mitochondrial genome molecular synapomorphies among clades

Mitochondrial chromosomes have diversified among eukaryotes and many different architectures and features are now acknowledged for this genome. Here we present the improved HERMES index, which can measure and quantify the amount of molecular change experienced by mitochondrial genomes. We test the improved approach with ten molecular phylogenetic studies based on complete mitochondrial genomes, representing six bilaterian Phyla. In most cases, HERMES analysis spotted out clades or single species with peculiar molecular synapomorphies, allowing to identify phylogenetic and ecological patterns. The software presented herein handles linear, circular, and multi-chromosome genomes, thus widening the HERMES scope to the complete eukaryotic domain.

Multigene phylogenetics of euglenids based on single-cell transcriptomics of diverse phagotrophs

Euglenids are a well-known group of single-celled eukaryotes, with phototrophic, osmotrophic and phagotrophic members. Phagotrophs represent most of the phylogenetic diversity of euglenids, and gave rise to the phototrophs and osmotrophs, but their evolutionary relationships are poorly understood. Symbiontids, in contrast, are anaerobes that are alternatively inferred to be derived euglenids, or a separate euglenozoan group. Most phylogenetic studies of euglenids have examined the SSU rDNA only, which is often highly divergent. Also, many phagotrophic euglenids (and symbiontids) are uncultured, restricting collection of other molecular data. We generated transcriptome data for 28 taxa, mostly using a single-cell approach, and conducted the first multigene phylogenetic analyses of euglenids to include phagotrophs and symbiontids. Euglenids are recovered as monophyletic, with symbiontids forming an independent branch within Euglenozoa. Spirocuta, the clade of flexible euglenids that contains both the phototrophs (Euglenophyceae) and osmotrophs (Aphagea), is robustly resolved, with the ploeotid Olkasia as its sister group, forming the new taxon Olkaspira. Ploeotids are paraphyletic, although Ploeotiidae (represented by Ploeotia spp.), Lentomonas, and Keelungia form a robust clade (new taxon Alistosa). Petalomonadida branches robustly as sister to other euglenids in outgroup-rooted analyses. Within Spirocuta, Euglenophyceae is a robust clade that includes Rapaza, and Anisonemia is a well-supported monophyletic group containing Anisonemidae (Anisonema and Dinema spp.), 'Heteronema II' (represented by H. vittatum), and a clade of Neometanema plus Aphagea. Among 'peranemid' phagotrophs, Chasmostoma branches with included Urceolus, and Peranema with the undescribed 'Jenningsia II', while other relationships are weakly supported and consequently the closest sister group to Euglenophyceae remains unresolved. Our results are inconsistent with recent inferences that Entosiphon is the evolutionarily pivotal sister either to other euglenids, or to Spirocuta. At least three transitions between posterior and anterior flagellar gliding occurred in euglenids, with the phylogenetic positions and directions of those transitions remaining ambiguous.

Nucleosome landscape reflects phenotypic differences in Trypanosoma cruzi life forms

Trypanosoma cruzi alternates between replicative and nonreplicative life forms, accompanied by a shift in global transcription levels and by changes in the nuclear architecture, the chromatin proteome and histone posttranslational modifications. To gain further insights into the epigenetic regulation that accompanies life form changes, we performed genome-wide high-resolution nucleosome mapping using two T. cruzi life forms (epimastigotes and cellular trypomastigotes). By combining a powerful pipeline that allowed us to faithfully compare nucleosome positioning and occupancy, more than 125 thousand nucleosomes were mapped, and approximately 20% of them differed between replicative and nonreplicative forms. The nonreplicative forms have less dynamic nucleosomes, possibly reflecting their lower global transcription levels and DNA replication arrest. However, dynamic nucleosomes are enriched at nonreplicative regulatory transcription initiation regions and at multigenic family members, which are associated with infective-stage and virulence factors. Strikingly, dynamic nucleosome regions are associated with GO terms related to nuclear division, translation, gene regulation and metabolism and, notably, associated with transcripts with different expression levels among life forms. Finally, the nucleosome landscape reflects the steady-state transcription expression: more abundant genes have a more deeply nucleosome-depleted region at putative 5' splice sites, likely associated with trans-splicing efficiency. Taken together, our results indicate that chromatin architecture, defined primarily by nucleosome positioning and occupancy, reflects the phenotypic differences found among T. cruzi life forms despite the lack of a canonical transcriptional control context.

Identification of Evolutionarily Conserved Nuclear Matrix Proteins and Their Prokaryotic Origins

Compared to prokaryotic cells, a typical eukaryotic cell is much more complex along with its endomembrane system and membrane-bound organelles. Although the endosymbiosis theories convincingly explain the evolution of membrane-bound organelles such as mitochondria and chloroplasts, very little is understood about the evolutionary origins of the nucleus, the defining feature of eukaryotes. Most studies on nuclear evolution have not been able to take into consideration the underlying structural framework of the nucleus, attributed to the nuclear matrix (NuMat), a ribonucleoproteinaceous structure. This can largely be attributed to the lack of annotation of its core components. Since NuMat has been shown to provide a structural platform for facilitating a variety of nuclear functions such as replication, transcription, and splicing, it is important to identify its protein components to better understand these processes. In this study, we address this issue using the developing embryos of Drosophila melanogaster and Danio rerio and identify 362 core NuMat proteins that are conserved between the two organisms. We further compare our results with publicly available Mus musculus NuMat dataset and Homo sapiens cellular localization dataset to define the core homologous NuMat proteins consisting of 252 proteins. We find that of them, 86 protein groups have originated from pre-existing proteins in prokaryotes. While 36 were conserved across all eukaryotic supergroups, 14 new proteins evolved before the evolution of the last eukaryotic common ancestor and together, these 50 proteins out of the 252 core conserved NuMat proteins are conserved across all eukaryotes, indicating their indispensable nature for nuclear function for over 1.5 billion years of eukaryotic history. Our analysis paves the way to understand the evolution of the complex internal nuclear architecture and its functions.

Description of Imasa heleensis, gen. nov., sp. nov. (Imasidae, fam. nov.), a Deep-Branching Marine Malawimonad and Possible Key Taxon in Understanding Early Eukaryotic Evolution

Malawimonadida is a deep-level (arguably "kingdom-scale") lineage of eukaryotes whose phylogenetic affinities are uncertain but of great evolutionary interest, as the group is suspected to branch close to the root of the tree of eukaryotes. Part of the difficulty in placing Malawimonadida phylogenetically is its tiny circumscription: at present, it comprises only two described and one cultured but undescribed species, all of them are freshwater suspension-feeding nanoflagellates. In this study, we cultivated and characterised Imasa heleensis gen. nov., sp. nov. (Imasidae fam. nov.), the first marine malawimonad to be described. Light and electron microscopy observations show that Imasa is largely similar to other malawimonads, but more frequently adheres to the substrate, often by means of a pliable posterior extension. Phylogenetic analyses based on two ribosomal RNA genes and four translated protein-coding genes using three different taxon sets place Imasa as sister to the three freshwater malawimonad strains with strong support. Imasa's mitochondrial genome is circular-mapping and shows a similar gene complement to other known malawimonads. We conclude that Imasa represents an important expansion of the range of taxa available for future evolutionary study.

Evolutionary History of Mitochondrial Genomes in Discoba, Including the Extreme Halophile Pleurostomum flabellatum (Heterolobosea)

Data from Discoba (Heterolobosea, Euglenozoa, Tsukubamonadida, and Jakobida) are essential to understand the evolution of mitochondrial genomes (mitogenomes), because this clade includes the most primitive-looking mitogenomes known, as well some extremely divergent genome information systems. Heterolobosea encompasses more than 150 described species, many of them from extreme habitats, but only six heterolobosean mitogenomes have been fully sequenced to date. Here we complete the mitogenome of the heterolobosean Pleurostomum flabellatum, which is extremely halophilic and reportedly also lacks classical mitochondrial cristae, hinting at reduction or loss of respiratory function. The mitogenome of P. flabellatum maps as a 57,829-bp-long circular molecule, including 40 coding sequences (19 tRNA, two rRNA, and 19 orfs). The gene content and gene arrangement are similar to Naegleria gruberi and Naegleria fowleri, the closest relatives with sequenced mitogenomes. The P. flabellatum mitogenome contains genes that encode components of the electron transport chain similar to those of Naegleria mitogenomes. Homology searches against a draft nuclear genome showed that P. flabellatum has two homologs of the highly conserved Mic60 subunit of the MICOS complex, and likely lost Mic19 and Mic10. However, electron microscopy showed no cristae structures. We infer that P. flabellatum, which originates from high salinity (313‰) water where the dissolved oxygen concentration is low, possesses a mitochondrion capable of aerobic respiration, but with reduced development of cristae structure reflecting limited use of this aerobic capacity (e.g., microaerophily).

Phylogenies from mitochondrial genomes of 120 species of ticks: Insights into the evolution of the families of ticks and of the genus Amblyomma

The evolution and phylogenetic relationships of the ticks at both the family and genus levels are contested. The genus Amblyomma and its subgenera are in a state of flux; moreover, the relationships among the three tick families are controversial due to conflicting phylogenetic support for different arrangements of the three families of living ticks. With 18 newly sequenced mitochondrial (mt) genomes of ticks included, we executed the largest mt genome phylogenetic study of ticks so far. Phylogenetic trees were inferred from one sea spider mt genome, one horseshoe crab, five mite mt genomes and 146 tick mt genomes from 120 species: 153 mt genomes in total. Sixteen phylogenetic trees were inferred from 10 datasets using both maximum likelihood and Bayesian inference methods. We describe the first novel mt gene-arrangement for the metastriate Ixodidae in Amblyomma (Africaniella) transversale. Also, three unusual partial 16S rRNA gene inserts were found in the mt genome of Haemaphysalis (Alloceraea) kitaokai: we consider the possible role of past genome translocation events in the formation of these inserts. Our phylogenies revealed evidence that: (i) the genus Amblyomma is polyphyletic with respect to Amblyomma (Africaniella) transversale; (ii) the subgenus Aponomma is apparently embedded in the genus Amblyomma; (iii) Haemaphysalis (Segalia) parva and Haemaphysalis (Alloceraea) kitaokai form a clade to the exclusion of other Haemaphysalis species; and (iv) the phylogenetic position of the family Nuttalliellidae is unstable among phylogenies from different datasets.

The Molecular Diversity of Phagotrophic Euglenids Examined Using Single-cell Methods

Euglenids are a diverse group of euglenozoan flagellates that includes phototrophs, osmotrophs, and phagotrophs. Despite making up most of the phylogenetic diversity of euglenids, phagotrophs remain understudied, and recent work has focused on 'deep-branching' groups. Spirocuta is the large clade encompassing all flexible euglenids including the phototroph and primary osmotroph clades, plus various phagotrophs. Understanding the phylogenetic diversity of phagotrophic spirocutes is crucial for tracing euglenid evolution, including how phototrophs arose. We used single-cell approaches to greatly increase sampling of SSU rDNA for phagotrophic euglenids, particularly spirocutes, including the first sequences from Urceolus, Jenningsia, Chasmostoma, and Sphenomonas, and expanded coverage for Dinema and Heteronema sensu lato, amongst others. Urceolus monophyly is unconfirmed. Organisms referred to Jenningsia form two distinct clades. Heteronema vittatum and similar cells branch separately from Heteronema (c.f.) globuliferum and Teloprocta/Heteronema scaphurum, while Dinema appears as 2-3 clades. Sphenomonas is monophyletic and the deepest branch within Petalomonadida. The census of genera markedly underestimates the phylogenetic diversity of phagotrophs, but taxonomic restraint is necessary when sequences are not available from type species or reasonable surrogates. SSU rDNA phylogenies do not resolve most deep relationships within Spirocuta, but identify units of diversity to sample in future multigene analyses.

Progress towards the Tree of Eukaryotes

Developing a detailed understanding of how all known forms of life are related to one another in the tree of life has been a major preoccupation of biology since the idea of tree-like evolution first took hold. Since most life is microbial, our intuitive use of morphological comparisons to infer relatedness only goes so far, and molecular sequence data, most recently from genomes and transcriptomes, has been the primary means to infer these relationships. For prokaryotes this presented new challenges, since the degree of horizontal gene transfer led some to question the tree-like depiction of evolution altogether. Most eukaryotes are also microbial, but in contrast to prokaryotic life, the application of large-scale molecular data to the tree of eukaryotes has largely been a constructive process, leading to a small number of very diverse lineages, or 'supergroups'. The tree is not completely resolved, and contentious problems remain, but many well-established supergroups now encompass much more diversity than the traditional kingdoms. Some of the most exciting recent developments come from the discovery of branches in the tree that we previously had no inkling even existed, many of which are of great ecological or evolutionary interest. These new branches highlight the need for more exploration, by high-throughput molecular surveys, but also more traditional means of observations and cultivation.

Identification and evolutionary analysis of the nucleolar proteome of Giardia lamblia

Background

The nucleoli, including their proteomes, of higher eukaryotes have been extensively studied, while few studies about the nucleoli of the lower eukaryotes – protists were reported. Giardia lamblia, a protist with the controversy of whether it is an extreme primitive eukaryote or just a highly evolved parasite, might be an interesting object for carrying out the nucleolar proteome study of protists and for further examining the controversy.

Results

Using bioinformatics methods, we reconstructed G. lamblia nucleolar proteome (GiNuP) and the common nucleolar proteome of the three representative higher eukaryotes (human, Arabidopsis, yeast) (HEBNuP). Comparisons of the two proteomes revealed that: 1) GiNuP is much smaller than HEBNuP, but 78.4% of its proteins have orthologs in the latter; 2) More than 68% of the GiNuP proteins are involved in the “Ribosome related” function, and the others participate in the other functions, and these two groups of proteins are much larger and much smaller than those in HEBNuP, respectively; 3) Both GiNuP and HEBNuP have their own specific proteins, but HEBNuP has a much higher proportion of such proteins to participate in more categories of nucleolar functions.

Conclusion

For the first time the nucleolar proteome of a protist - Giardia was reconstructed. The results of comparison of it with the common proteome of three representative higher eukaryotes -- HEBNuP indicated that the simplicity of GiNuP is most probably a reflection of primitiveness but not just parasitic reduction of Giardia, and simultaneously revealed some interesting evolutionary phenomena about the nucleolus and even the eukaryotic cell, compositionally and functionally.

Evolution of metabolic capabilities and molecular features of diplonemids, kinetoplastids, and euglenids

Background

The Euglenozoa are a protist group with an especially rich history of evolutionary diversity. They include diplonemids, representing arguably the most species-rich clade of marine planktonic eukaryotes; trypanosomatids, which are notorious parasites of medical and veterinary importance; and free-living euglenids. These different lifestyles, and particularly the transition from free-living to parasitic, likely require different metabolic capabilities. We carried out a comparative genomic analysis across euglenozoan diversity to see how changing repertoires of enzymes and structural features correspond to major changes in lifestyles.

Results

We find a gradual loss of genes encoding enzymes in the evolution of kinetoplastids, rather than a sudden decrease in metabolic capabilities corresponding to the origin of parasitism, while diplonemids and euglenids maintain more metabolic versatility. Distinctive characteristics of molecular machines such as kinetochores and the pre-replication complex that were previously considered specific to parasitic kinetoplastids were also identified in their free-living relatives. Therefore, we argue that they represent an ancestral rather than a derived state, as thought until the present. We also found evidence of ancient redundancy in systems such as NADPH-dependent thiol-redox. Only the genus Euglena possesses the combination of trypanothione-, glutathione-, and thioredoxin-based systems supposedly present in the euglenozoan common ancestor, while other representatives of the phylum have lost one or two of these systems. Lastly, we identified convergent losses of specific metabolic capabilities between free-living kinetoplastids and ciliates. Although this observation requires further examination, it suggests that certain eukaryotic lineages are predisposed to such convergent losses of key enzymes or whole pathways.

Conclusions

The loss of metabolic capabilities might not be associated with the switch to parasitic lifestyle in kinetoplastids, and the presence of a highly divergent (or unconventional) kinetochore machinery might not be restricted to this protist group. The data derived from the transcriptomes of free-living early branching prokinetoplastids suggests that the pre-replication complex of Trypanosomatidae is a highly divergent version of the conventional machinery. Our findings shed light on trends in the evolution of metabolism in protists in general and open multiple avenues for future research.

The draft nuclear genome sequence and predicted mitochondrial proteome of Andalucia godoyi, a protist with the most gene-rich and bacteria-like mitochondrial genome

Background

Comparative analyses have indicated that the mitochondrion of the last eukaryotic common ancestor likely possessed all the key core structures and functions that are widely conserved throughout the domain Eucarya. To date, such studies have largely focused on animals, fungi, and land plants (primarily multicellular eukaryotes); relatively few mitochondrial proteomes from protists (primarily unicellular eukaryotic microbes) have been examined. To gauge the full extent of mitochondrial structural and functional complexity and to identify potential evolutionary trends in mitochondrial proteomes, more comprehensive explorations of phylogenetically diverse mitochondrial proteomes are required. In this regard, a key group is the jakobids, a clade of protists belonging to the eukaryotic supergroup Discoba, distinguished by having the most gene-rich and most bacteria-like mitochondrial genomes discovered to date.

Results

In this study, we assembled the draft nuclear genome sequence for the jakobid Andalucia godoyi and used a comprehensive in silico approach to infer the nucleus-encoded portion of the mitochondrial proteome of this protist, identifying 864 candidate mitochondrial proteins. The A. godoyi mitochondrial proteome has a complexity that parallels that of other eukaryotes, while exhibiting an unusually large number of ancestral features that have been lost particularly in opisthokont (animal and fungal) mitochondria. Notably, we find no evidence that the A. godoyi nuclear genome has or had a gene encoding a single-subunit, T3/T7 bacteriophage-like RNA polymerase, which functions as the mitochondrial transcriptase in all eukaryotes except the jakobids.

Conclusions

As genome and mitochondrial proteome data have become more widely available, a strikingly punctuate phylogenetic distribution of different mitochondrial components has been revealed, emphasizing that the pathways of mitochondrial proteome evolution are likely complex and lineage-specific. Unraveling this complexity will require comprehensive comparative analyses of mitochondrial proteomes from a phylogenetically broad range of eukaryotes, especially protists. The systematic in silico approach described here offers a valuable adjunct to direct proteomic analysis (e.g., via mass spectrometry), particularly in cases where the latter approach is constrained by sample limitation or other practical considerations.

Dur3 and nrt2 genes in the bloom-forming dinoflagellate Prorocentrum minimum: Transcriptional responses to available nitrogen sources

The increasing inflow of nitrogen (N) substrates into marine nearshore ecosystems induces proliferation of harmful algal blooms (HABs) of dinoflagellates, such as potentially toxic invasive species Prorocentrum minimum. In this study, we estimated the influence of NO3-, NH4+ and urea on transcription levels and urea transporter dur3 and nitrate transporter nrt2 genes expression in these dinoflagellates. We identified dur3 and nrt2 genes sequences in unannotated transcriptomes of P. minimum and other dinoflagellates presented in MMETSP database. Phylogenetic analysis showed that these genes of dinoflagellates clustered to the distinct clade demonstrating evolutionary relationship with the other known dur3 and nrt2 genes of microalgae. The evaluation of expression levels of dur3 and nrt2 genes by RT-qPCR revealed their sensitivity to input of the studied N sources. Dur3 expression levels were downregulated after the supplementation of additional N sources and were 1.7-2.6-fold lower than in the nitrate-grown culture. Nrt2 expression levels decreased 1.9-fold in the presence of NH4+. We estimated total RNA and DNA synthesis rates by the analysis of incorporation of 3H-thymidine and 3H-uridine in batch and continuous cultures. Addition of N compounds did not affect the DNA synthesis rates. Transcription levels increased up to 12.5-fold after the N supplementation in urea-limited treatments. Investigation of various nitrogen sources as biomarkers of dinoflagellate proliferation due to their differentiated impact on expression of dur3 and nrt2 genes and transcription rates in P. minimum cells allowed concluding about high potential of the studied parameters for future modeling of HABs under global N pollution.

EukRef-excavates: seven curated SSU ribosomal RNA gene databases

The small subunit ribosomal RNA (SSU rRNA) gene is a widely used molecular marker to study the diversity of life. Sequencing of SSU rRNA gene amplicons has become a standard approach for the investigation of the ecology and diversity of microbes. However, a well-curated database is necessary for correct classification of these data. While available for many groups of Bacteria and Archaea, such reference databases are absent for most eukaryotes. The primary goal of the EukRef project (eukref.org) is to close this gap and generate well-curated reference databases for major groups of eukaryotes, especially protists. Here we present a set of EukRef-curated databases for the excavate protists—a large assemblage that includes numerous taxa with divergent SSU rRNA gene sequences, which are prone to misclassification. We identified 6121 sequences, 625 of which were obtained from cultures, 3053 from cell isolations or enrichments and 2419 from environmental samples. We have corrected the classification for the majority of these curated sequences. The resulting publicly available databases will provide phylogenetically based standards for the improved identification of excavates in ecological and microbiome studies, as well as resources to classify new discoveries in excavate diversity.

Cryo electron tomography with volta phase plate reveals novel structural foundations of the 96-nm axonemal repeat in the pathogen Trypanosoma brucei

The 96-nm axonemal repeat includes dynein motors and accessory structures as the foundation for motility of eukaryotic flagella and cilia. However, high-resolution 3D axoneme structures are unavailable for organisms among the Excavates, which include pathogens of medical and economic importance. Here we report cryo electron tomography structures of the 96-nm repeat from Trypanosoma brucei, a protozoan parasite in the Excavate lineage that causes African trypanosomiasis. We examined bloodstream and procyclic life cycle stages, and a knockdown lacking DRC11/CMF22 of the nexin dynein regulatory complex (NDRC). Sub-tomogram averaging yields a resolution of 21.8 Å for the 96-nm repeat. We discovered several lineage-specific structures, including novel inter-doublet linkages and microtubule inner proteins (MIPs). We establish that DRC11/CMF22 is required for the NDRC proximal lobe that binds the adjacent doublet microtubule. We propose that lineage-specific elaboration of axoneme structure in T. brucei reflects adaptations to support unique motility needs in diverse host environments.

Monitoring of the toxic dinoflagellate Alexandrium catenella in Osaka Bay, Japan using a massively parallel sequencing (MPS)-based technique

Since 2002, blooms of Alexandrium catenella sensu Fraga et al. (2015) and paralytic shellfish toxicity events have occurred almost yearly in Osaka Bay, Japan. To better understand the triggers for reoccurring A. catenella blooms in Osaka Bay, phytoplankton community was monitored during the spring seasons of 2012-2015. Monitoring was performed using massively parallel sequencing (MPS)-based technique on amplicon sequences of the 18S rRNA gene. Dense blooms of A. catenella occurred every year except in 2012, however, there was no significant correlation with the environmental parameters investigated. Plankton community diversity decreased before and middle of the A. catenella blooms, suggesting that the decline in diversity could be an indicator for the bloom occurrence. The yearly abundance pattern of A. catenella cells obtained by morphology-based counting coincided with the relative sequence abundances, which supports the effectiveness of MPS-based phytoplankton monitoring.

The New Tree of Eukaryotes

For 15 years, the eukaryote Tree of Life (eToL) has been divided into five to eight major groupings, known as 'supergroups'. However, the tree has been profoundly rearranged during this time. The new eToL results from the widespread application of phylogenomics and numerous discoveries of major lineages of eukaryotes, mostly free-living heterotrophic protists. The evidence that supports the tree has transitioned from a synthesis of molecular phylogenetics and biological characters to purely molecular phylogenetics. Most current supergroups lack defining morphological or cell-biological characteristics, making the supergroup label even more arbitrary than before. Going forward, the combination of traditional culturing with maturing culture-free approaches and phylogenomics should accelerate the process of completing and resolving the eToL at its deepest levels.

Fe-S Cluster Assembly in Oxymonads and Related Protists

The oxymonad Monocercomonoides exilis was recently reported to be the first eukaryote that has completely lost the mitochondrial compartment. It was proposed that an important prerequisite for such a radical evolutionary step was the acquisition of the SUF Fe–S cluster assembly pathway from prokaryotes, making the mitochondrial ISC pathway dispensable. We have investigated genomic and transcriptomic data from six oxymonad species and their relatives, composing the group Preaxostyla (Metamonada, Excavata), for the presence and absence of enzymes involved in Fe–S cluster biosynthesis. None possesses enzymes of mitochondrial ISC pathway and all apparently possess the SUF pathway, composed of SufB, C, D, S, and U proteins, altogether suggesting that the transition from ISC to SUF preceded their last common ancestor. Interestingly, we observed that SufDSU were fused in all three oxymonad genomes, and in the genome of Paratrimastix pyriformis. The donor of the SUF genes is not clear from phylogenetic analyses, but the enzyme composition of the pathway and the presence of SufDSU fusion suggests Firmicutes, Thermotogae, Spirochaetes, Proteobacteria, or Chloroflexi as donors. The inventory of the downstream CIA pathway enzymes is consistent with that of closely related species that retain ISC, indicating that the switch from ISC to SUF did not markedly affect the downstream process of maturation of cytosolic and nuclear Fe–S proteins.

Comparative molecular cell biology of phototrophic euglenids and parasitic trypanosomatids sheds light on the ancestor of Euglenozoa

Parasitic trypanosomatids and phototrophic euglenids are among the most extensively studied euglenozoans. The phototrophic euglenid lineage arose relatively recently through secondary endosymbiosis between a phagotrophic euglenid and a prasinophyte green alga that evolved into the euglenid secondary chloroplast. The parasitic trypanosomatids (i.e. Trypanosoma spp. and Leishmania spp.) and the freshwater phototrophic euglenids (i.e. Euglena gracilis) are the most evolutionary distant lineages in the Euglenozoa phylogenetic tree. The molecular and cell biological traits they share can thus be considered as ancestral traits originating in the common euglenozoan ancestor. These euglenozoan ancestral traits include common mitochondrial presequence motifs, respiratory chain complexes containing various unique subunits, a unique ATP synthase structure, the absence of mitochondria-encoded transfer RNAs (tRNAs), a nucleus with a centrally positioned nucleolus, closed mitosis without dissolution of the nuclear membrane and nucleoli, a nuclear genome containing the unusual 'J' base (β-D-glucosyl-hydroxymethyluracil), processing of nucleus-encoded precursor messenger RNAs (pre-mRNAs) via spliced-leader RNA (SL-RNA) trans-splicing, post-transcriptional gene silencing by the RNA interference (RNAi) pathway and the absence of transcriptional regulation of nuclear gene expression. Mitochondrial uridine insertion/deletion RNA editing directed by guide RNAs (gRNAs) evolved in the ancestor of the kinetoplastid lineage. The evolutionary origin of other molecular features known to be present only in either kinetoplastids (i.e. polycistronic transcripts, compaction of nuclear genomes) or euglenids (i.e. monocistronic transcripts, huge genomes, many nuclear cis-spliced introns, polyproteins) is unclear.

Comparative genomic analysis of the 'pseudofungus' Hyphochytrium catenoides

Eukaryotic microbes have three primary mechanisms for obtaining nutrients and energy: phagotrophy, photosynthesis and osmotrophy. Traits associated with the latter two functions arose independently multiple times in the eukaryotes. The Fungi successfully coupled osmotrophy with filamentous growth, and similar traits are also manifested in the Pseudofungi (oomycetes and hyphochytriomycetes). Both the Fungi and the Pseudofungi encompass a diversity of plant and animal parasites. Genome-sequencing efforts have focused on host-associated microbes (mutualistic symbionts or parasites), providing limited comparisons with free-living relatives. Here we report the first draft genome sequence of a hyphochytriomycete ‘pseudofungus’; Hyphochytrium catenoides. Using phylogenomic approaches, we identify genes of recent viral ancestry, with related viral derived genes also present on the genomes of oomycetes, suggesting a complex history of viral coevolution and integration across the Pseudofungi. H. catenoides has a complex life cycle involving diverse filamentous structures and a flagellated zoospore with a single anterior tinselate flagellum. We use genome comparisons, drug sensitivity analysis and high-throughput culture arrays to investigate the ancestry of oomycete/pseudofungal characteristics, demonstrating that many of the genetic features associated with parasitic traits evolved specifically within the oomycete radiation. Comparative genomics also identified differences in the repertoire of genes associated with filamentous growth between the Fungi and the Pseudofungi, including differences in vesicle trafficking systems, cell-wall synthesis pathways and motor protein repertoire, demonstrating that unique cellular systems underpinned the convergent evolution of filamentous osmotrophic growth in these two eukaryotic groups.

Targeting of tail-anchored proteins to Trichomonas vaginalis hydrogenosomes

Tail-anchored (TA) proteins are membrane proteins that are found in all domains of life. They consist of an N-terminal domain that performs various functions and a single transmembrane domain (TMD) near the C-terminus. In eukaryotes, TA proteins are targeted to the membranes of mitochondria, the endoplasmic reticulum (ER), peroxisomes and in plants, chloroplasts. The targeting of these proteins to their specific destinations correlates with the properties of the C-terminal domain, mainly the TMD hydrophobicity and the net charge of the flanking regions. Trichomonas vaginalis is a human parasite that has adapted to oxygen-poor environment. This adaptation is reflected by the presence of highly modified mitochondria (hydrogenosomes) and the absence of peroxisomes. The proteome of hydrogenosomes is considerably reduced; however, our bioinformatic analysis predicted 120 putative hydrogenosomal TA proteins. Seven proteins were selected to prove their localization. The elimination of the net positive charge in the C-tail of the hydrogenosomal TA4 protein resulted in its dual localization to hydrogenosomes and the ER, causing changes in ER morphology. Domain mutation and swap experiments with hydrogenosomal (TA4) and ER (TAPDI) proteins indicated that the general principles for specific targeting are conserved across eukaryotic lineages, including T. vaginalis; however, there are also significant lineage-specific differences.