A Novel Group of Dynamin-Related Proteins Shared by Eukaryotes and Giant Viruses Is Able to Remodel Mitochondria From Within the Matrix.

The diverse GTPases of the dynamin superfamily play various roles in the cell, as exemplified by the dynamin-related proteins (DRPs) Mgm1 and Opa1, which remodel the mitochondrial inner membrane in fungi and metazoans, respectively. Via an exhaustive search of genomic and metagenomic databases, we found previously unknown DRP types occurring in diverse eukaryotes and giant viruses (phylum Nucleocytoviricota). One novel DRP clade, termed MidX, combined hitherto uncharacterized proteins from giant viruses and six distantly related eukaryote taxa (Stramenopiles, Telonemia, Picozoa, Amoebozoa, Apusomonadida, and Choanoflagellata). MidX stood out because it was not only predicted to be mitochondria-targeted but also to assume a tertiary structure not observed in other DRPs before. To understand how MidX affects mitochondria, we exogenously expressed MidX from Hyperionvirus in the kinetoplastid Trypanosoma brucei, which lacks Mgm1 or Opa1 orthologs. MidX massively affected mitochondrial morphology from inside the matrix, where it closely associates with the inner membrane. This unprecedented mode of action contrasts to those of Mgm1 and Opa1, which mediate inner membrane remodeling in the intermembrane space. We speculate that MidX was acquired in Nucleocytoviricota evolution by horizontal gene transfer from eukaryotes and is used by giant viruses to remodel host mitochondria during infection. MidX's unique structure may be an adaptation for reshaping mitochondria from the inside. Finally, Mgm1 forms a sister group to MidX and not Opa1 in our phylogenetic analysis, throwing into question the long-presumed homology of these DRPs with similar roles in sister lineages.

Single cell transcriptomics reveals UAR codon reassignment in Palmarella salina (Metopida, Armophorea) and confirms Armophorida belongs to APM clade.

Anaerobes have emerged in several major lineages of ciliates, but the number of independent transitions to anaerobiosis among ciliates is unknown. The APM clade (Armophorea, Muranotrichea, Parablepharismea) represents the largest clade of obligate anaerobes among ciliates and contains free-living marine and freshwater representatives as well as gut endobionts of animals. The evolution of APM group has only recently started getting attention, and our knowledge on its phylogeny and genetics is still limited to a fraction of taxa. While ciliates portray a wide array of alternatives to the standard genetic code across numerous classes, the APM ciliates were considered to be the largest group using exclusively standard nuclear genetic code. In this study, we present a pan-ciliate phylogenomic analysis with emphasis on the APM clade, bringing the first phylogenomic analysis of the family Tropidoatractidae (Armophorea) and confirming the position of Armophorida within Armophorea. We include five newly sequenced single cell transcriptomes from marine, freshwater, and endobiotic APM ciliates - Palmarella salina, Anteclevelandella constricta, Nyctotherus sp., Caenomorpha medusula, and Thigmothrix strigosa. We report the first discovery of an alternative nuclear genetic code among APM ciliates, used by Palmarella salina (Tropidoatractidae, Armophorea), but not by its close relative, Tropidoatractus sp., and provide a comparative analysis of stop codon identity and frequency indicating the precedency to the UAG codon loss/reassignment over the UAA codon reassignment in the specific ancestor of Palmarella. Comparative genomic and proteomic studies of this group may help explain the constraints that underlie UAR stop-to-sense reassignment, the most frequent type of alternative nuclear genetic code, not only in ciliates, but eukaryotes in general.

PhyloFisher: A phylogenomic package for resolving eukaryotic relationships.

Phylogenomic analyses of hundreds of protein-coding genes aimed at resolving phylogenetic relationships is now a common practice. However, no software currently exists that includes tools for dataset construction and subsequent analysis with diverse validation strategies to assess robustness. Furthermore, there are no publicly available high-quality curated databases designed to assess deep (>100 million years) relationships in the tree of eukaryotes. To address these issues, we developed an easy-to-use software package, PhyloFisher (https://github.com/TheBrownLab/PhyloFisher), written in Python 3. PhyloFisher includes a manually curated database of 240 protein-coding genes from 304 eukaryotic taxa covering known eukaryotic diversity, a novel tool for ortholog selection, and utilities that will perform diverse analyses required by state-of-the-art phylogenomic investigations. Through phylogenetic reconstructions of the tree of eukaryotes and of the Saccharomycetaceae clade of budding yeasts, we demonstrate the utility of the PhyloFisher workflow and the provided starting database to address phylogenetic questions across a large range of evolutionary time points for diverse groups of organisms. We also demonstrate that undetected paralogy can remain in phylogenomic "single-copy orthogroup" datasets constructed using widely accepted methods such as all vs. all BLAST searches followed by Markov Cluster Algorithm (MCL) clustering and application of automated tree pruning algorithms. Finally, we show how the PhyloFisher workflow helps detect inadvertent paralog inclusions, allowing the user to make more informed decisions regarding orthology assignments, leading to a more accurate final dataset.

An Enigmatic Stramenopile Sheds Light on Early Evolution in Ochrophyta Plastid Organellogenesis.

Ochrophyta is an algal group belonging to the Stramenopiles and comprises diverse lineages of algae which contribute significantly to the oceanic ecosystems as primary producers. However, early evolution of the plastid organelle in Ochrophyta is not fully understood. In this study, we provide a well-supported tree of the Stramenopiles inferred by the large-scale phylogenomic analysis that unveils the eukaryvorous (nonphotosynthetic) protist Actinophrys sol (Actinophryidae) is closely related to Ochrophyta. We used genomic and transcriptomic data generated from A. sol to detect molecular traits of its plastid and we found no evidence of plastid genome and plastid-mediated biosynthesis, consistent with previous ultrastructural studies that did not identify any plastids in Actinophryidae. Moreover, our phylogenetic analyses of particular biosynthetic pathways provide no evidence of a current and past plastid in A. sol. However, we found more than a dozen organellar aminoacyl-tRNA synthases (aaRSs) that are of algal origin. Close relationships between aaRS from A. sol and their ochrophyte homologs document gene transfer of algal genes that happened before the divergence of Actinophryidae and Ochrophyta lineages. We further showed experimentally that organellar aaRSs of A. sol are targeted exclusively to mitochondria, although organellar aaRSs in Ochrophyta are dually targeted to mitochondria and plastids. Together, our findings suggested that the last common ancestor of Actinophryidae and Ochrophyta had not yet completed the establishment of host-plastid partnership as seen in the current Ochrophyta species, but acquired at least certain nuclear-encoded genes for the plastid functions.

A new lineage of non-photosynthetic green algae with extreme organellar genomes.

BACKGROUND: The plastid genomes of the green algal order Chlamydomonadales tend to expand their non-coding regions, but this phenomenon is poorly understood. Here we shed new light on organellar genome evolution in Chlamydomonadales by studying a previously unknown non-photosynthetic lineage. We established cultures of two new Polytoma-like flagellates, defined their basic characteristics and phylogenetic position, and obtained complete organellar genome sequences and a transcriptome assembly for one of them. RESULTS: We discovered a novel deeply diverged chlamydomonadalean lineage that has no close photosynthetic relatives and represents an independent case of photosynthesis loss. To accommodate these organisms, we establish the new genus Leontynka, with two species (L. pallida and L. elongata) distinguishable through both their morphological and molecular characteristics. Notable features of the colourless plastid of L. pallida deduced from the plastid genome (plastome) sequence and transcriptome assembly include the retention of ATP synthase, thylakoid-associated proteins, the carotenoid biosynthesis pathway, and a plastoquinone-based electron transport chain, the latter two modules having an obvious functional link to the eyespot present in Leontynka. Most strikingly, the ~362 kbp plastome of L. pallida is by far the largest among the non-photosynthetic eukaryotes investigated to date due to an extreme proliferation of sequence repeats. These repeats are also present in coding sequences, with one repeat type found in the exons of 11 out of 34 protein-coding genes, with up to 36 copies per gene, thus affecting the encoded proteins. The mitochondrial genome of L. pallida is likewise exceptionally large, with its >104 kbp surpassed only by the mitogenome of Haematococcus lacustris among all members of Chlamydomonadales hitherto studied. It is also bloated with repeats, though entirely different from those in the L. pallida plastome, which contrasts with the situation in H. lacustris where both the organellar genomes have accumulated related repeats. Furthermore, the L. pallida mitogenome exhibits an extremely high GC content in both coding and non-coding regions and, strikingly, a high number of predicted G-quadruplexes. CONCLUSIONS: With its unprecedented combination of plastid and mitochondrial genome characteristics, Leontynka pushes the frontiers of organellar genome diversity and is an interesting model for studying organellar genome evolution.

Vestiges of the Bacterial Signal Recognition Particle-Based Protein Targeting in Mitochondria.

The main bacterial pathway for inserting proteins into the plasma membrane relies on the signal recognition particle (SRP), composed of the Ffh protein and an associated RNA component, and the SRP-docking protein FtsY. Eukaryotes use an equivalent system of archaeal origin to deliver proteins into the endoplasmic reticulum, whereas a bacteria-derived SRP and FtsY function in the plastid. Here we report on the presence of homologs of the bacterial Ffh and FtsY proteins in various unrelated plastid-lacking unicellular eukaryotes, namely Heterolobosea, Alveida, Goniomonas, and Hemimastigophora. The monophyly of novel eukaryotic Ffh and FtsY groups, predicted mitochondrial localization experimentally confirmed for Naegleria gruberi, and a strong alphaproteobacterial affinity of the Ffh group, collectively suggest that they constitute parts of an ancestral mitochondrial signal peptide-based protein-targeting system inherited from the last eukaryotic common ancestor, but lost from the majority of extant eukaryotes. The ability of putative signal peptides, predicted in a subset of mitochondrial-encoded N. gruberi proteins, to target a reporter fluorescent protein into the endoplasmic reticulum of Trypanosoma brucei, likely through their interaction with the cytosolic SRP, provided further support for this notion. We also illustrate that known mitochondrial ribosome-interacting proteins implicated in membrane protein targeting in opisthokonts (Mba1, Mdm38, and Mrx15) are broadly conserved in eukaryotes and nonredundant with the mitochondrial SRP system. Finally, we identified a novel mitochondrial protein (MAP67) present in diverse eukaryotes and related to the signal peptide-binding domain of Ffh, which may well be a hitherto unrecognized component of the mitochondrial membrane protein-targeting machinery.

The Oxymonad Genome Displays Canonical Eukaryotic Complexity in the Absence of a Mitochondrion.

The discovery that the protist Monocercomonoides exilis completely lacks mitochondria demonstrates that these organelles are not absolutely essential to eukaryotic cells. However, the degree to which the metabolism and cellular systems of this organism have adapted to the loss of mitochondria is unknown. Here, we report an extensive analysis of the M. exilis genome to address this question. Unexpectedly, we find that M. exilis genome structure and content is similar in complexity to other eukaryotes and less "reduced" than genomes of some other protists from the Metamonada group to which it belongs. Furthermore, the predicted cytoskeletal systems, the organization of endomembrane systems, and biosynthetic pathways also display canonical eukaryotic complexity. The only apparent preadaptation that permitted the loss of mitochondria was the acquisition of the SUF system for Fe-S cluster assembly and the loss of glycine cleavage system. Changes in other systems, including in amino acid metabolism and oxidative stress response, were coincident with the loss of mitochondria but are likely adaptations to the microaerophilic and endobiotic niche rather than the mitochondrial loss per se. Apart from the lack of mitochondria and peroxisomes, we show that M. exilis is a fully elaborated eukaryotic cell that is a promising model system in which eukaryotic cell biology can be investigated in the absence of mitochondria.

Between a Pod and a Hard Test: The Deep Evolution of Amoebae.

Amoebozoa is the eukaryotic supergroup sister to Obazoa, the lineage that contains the animals and Fungi, as well as their protistan relatives, and the breviate and apusomonad flagellates. Amoebozoa is extraordinarily diverse, encompassing important model organisms and significant pathogens. Although amoebozoans are integral to global nutrient cycles and present in nearly all environments, they remain vastly understudied. We present a robust phylogeny of Amoebozoa based on broad representative set of taxa in a phylogenomic framework (325 genes). By sampling 61 taxa using culture-based and single-cell transcriptomics, our analyses show two major clades of Amoebozoa, Discosea, and Tevosa. This phylogeny refutes previous studies in major respects. Our results support the hypothesis that the last common ancestor of Amoebozoa was sexual and flagellated, it also may have had the ability to disperse propagules from a sporocarp-type fruiting body. Overall, the main macroevolutionary patterns in Amoebozoa appear to result from the parallel losses of homologous characters of a multiphase life cycle that included flagella, sex, and sporocarps rather than independent acquisition of convergent features.

Nuclear genetic codes with a different meaning of the UAG and the UAA codon.

BACKGROUND: Departures from the standard genetic code in eukaryotic nuclear genomes are known for only a handful of lineages and only a few genetic code variants seem to exist outside the ciliates, the most creative group in this regard. Most frequent code modifications entail reassignment of the UAG and UAA codons, with evidence for at least 13 independent cases of a coordinated change in the meaning of both codons. However, no change affecting each of the two codons separately has been documented, suggesting the existence of underlying evolutionary or mechanistic constraints. RESULTS: Here, we present the discovery of two new variants of the nuclear genetic code, in which UAG is translated as an amino acid while UAA is kept as a termination codon (along with UGA). The first variant occurs in an organism noticed in a (meta)transcriptome from the heteropteran Lygus hesperus and demonstrated to be a novel insect-dwelling member of Rhizaria (specifically Sainouroidea). This first documented case of a rhizarian with a non-canonical genetic code employs UAG to encode leucine and represents an unprecedented change among nuclear codon reassignments. The second code variant was found in the recently described anaerobic flagellate Iotanema spirale (Metamonada: Fornicata). Analyses of transcriptomic data revealed that I. spirale uses UAG to encode glutamine, similarly to the most common variant of a non-canonical code known from several unrelated eukaryotic groups, including hexamitin diplomonads (also a lineage of fornicates). However, in these organisms, UAA also encodes glutamine, whereas it is the primary termination codon in I. spirale. Along with phylogenetic evidence for distant relationship of I. spirale and hexamitins, this indicates two independent genetic code changes in fornicates. CONCLUSIONS: Our study documents, for the first time, that evolutionary changes of the meaning of UAG and UAA codons in nuclear genomes can be decoupled and that the interpretation of the two codons by the cytoplasmic translation apparatus is mechanistically separable. The latter conclusion has interesting implications for possibilities of genetic code engineering in eukaryotes. We also present a newly developed generally applicable phylogeny-informed method for inferring the meaning of reassigned codons.

First multigene analysis of Archamoebae (Amoebozoa: Conosa) robustly reveals its phylogeny and shows that Entamoebidae represents a deep lineage of the group.

Archamoebae is an understudied group of anaerobic free-living or endobiotic protists that constitutes the major anaerobic lineage of the supergroup Amoebozoa. Hitherto, the phylogeny of Archamoebae was based solely on SSU rRNA and actin genes, which did not resolve relationships among the main lineages of the group. Because of this uncertainty, several different scenarios had been proposed for the phylogeny of the Archamoebae. In this study, we present the first multigene phylogenetic analysis that includes members of Pelomyxidae, and Rhizomastixidae. The analysis clearly shows that Mastigamoebidae, Pelomyxidae and Rhizomastixidae form a clade of mostly free-living, amoeboid flagellates, here called Pelobiontida. The predominantly endobiotic and aflagellated Entamoebidae represents a separate, deep-branching lineage, Entamoebida. Therefore, two unique evolutionary events, horizontal transfer of the nitrogen fixation system from bacteria and transfer of the sulfate activation pathway to mitochondrial derivatives, predate the radiation of recent lineages of Archamoebae. The endobiotic lifestyle has arisen at least three times independently during the evolution of the group. We also present new ultrastructural data that clarifies the primary divergence among the family Mastigamoebidae which had previously been inferred from phylogenetic analyses based on SSU rDNA.

Morphological Identities of Two Different Marine Stramenopile Environmental Sequence Clades: Bicosoeca kenaiensis (Hilliard, 1971) and Cantina marsupialis (Larsen and Patterson, 1990) gen. nov., comb. nov.

Although environmental DNA surveys improve our understanding of biodiversity, interpretation of unidentified lineages is limited by the absence of associated morphological traits and living cultures. Unidentified lineages of marine stramenopiles are called "MAST clades". Twenty-five MAST clades have been recognized: MAST-1 through MAST-25; seven of these have been subsequently discarded because the sequences representing those clades were found to either (1) be chimeric or (2) affiliate within previously described taxonomic groups. Eighteen MAST clades remain without a cellular identity. Moreover, the discarded "MAST-13" has been used in different studies to refer to two different environmental sequence clades. After establishing four cultures representing two different species of heterotrophic stramenopiles and then characterizing their morphology and molecular phylogenetic positions, we determined that the two different species represented the two different MAST-13 clades: (1) a lorica-bearing Bicosoeca kenaiensis and (2) a microaerophilic flagellate previously named "Cafeteria marsupialis". Both species were previously described with only light microscopy; no cultures, ultrastructural data or DNA sequences were available from these species prior to this study. The molecular phylogenetic position of three different "C. marsupialis" isolates was not closely related to the type species of Cafeteria; therefore, we established a new genus for these isolates, Cantina gen. nov.

Survey on diversity of marine/saline anaerobic Heterolobosea (Excavata: Discoba) with description of seven new species.

Diversity of the anaerobic Heterolobosea (Excavata: Discoba) is only poorly understood, especially in marine environments. We have isolated and cultured 16 strains of anaerobic heteroloboseid amoebae and flagellates from brackish, marine and saline anoxic habitats worldwide. Phylogenetic analyses of SSU rDNA sequences and light-microscopic observations showed that all the strains belong to the family Psalteriomonadidae, the main anaerobic lineage of Heterolobosea, and that they represent eight species from the genera Monopylocystis, Harpagon and Pseudoharpagon. Seven species are newly isolated and described here as Monopylocystis minor n. sp., Monopylocystis robusta n. sp., Monopylocystis elegans n. sp., Monopylocystis disparata n. sp., Harpagon salinus n. sp., Pseudoharpagon longus n. sp. and Pseudoharpagon tertius n. sp. Amoebae, cysts and the ultrastructure of the genus Pseudoharpagon are presented for the first time.

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.

Create, Analyze, and Visualize Phylogenomic Datasets Using PhyloFisher.

PhyloFisher is a software package written primarily in Python3 that can be used for the creation, analysis, and visualization of phylogenomic datasets that consist of protein sequences from eukaryotic organisms. Unlike many existing phylogenomic pipelines, PhyloFisher comes with a manually curated database of 240 protein-coding genes, a subset of a previous phylogenetic dataset sampled from 304 eukaryotic taxa. The software package can also utilize a user-created database of eukaryotic proteins, which may be more appropriate for shallow evolutionary questions. PhyloFisher is also equipped with a set of utilities to aid in running routine analyses, such as the prediction of alternative genetic codes, removal of genes and/or taxa based on occupancy/completeness of the dataset, testing for amino acid compositional heterogeneity among sequences, removal of heterotachious and/or fast-evolving sites, removal of fast-evolving taxa, supermatrix creation from randomly resampled genes, and supermatrix creation from nucleotide sequences. © 2024 Wiley Periodicals LLC. Basic Protocol 1: Constructing a phylogenomic dataset Basic Protocol 2: Performing phylogenomic analyses Support Protocol 1: Installing PhyloFisher Support Protocol 2: Creating a custom phylogenomic database.

The protist cultural renaissance.

Protists are key players in the biosphere. Here, we provide a perspective on integrating protist culturing with omics approaches, imaging, and high-throughput single-cell manipulation strategies, concluding with actions required for a successful return of the golden age of protist culturing.

Revisiting biocrystallization: purine crystalline inclusions are widespread in eukaryotes.

Despite the widespread occurrence of intracellular crystalline inclusions in unicellular eukaryotes, scant attention has been paid to their composition, functions, and evolutionary origins. Using Raman microscopy, we examined >200 species from all major eukaryotic supergroups. We detected cellular crystalline inclusions in 77% species out of which 80% is composed of purines, such as anhydrous guanine (62%), guanine monohydrate (2%), uric acid (12%) and xanthine (4%). Our findings shifts the paradigm assuming predominance of calcite and oxalates. Purine crystals emerge in microorganisms in all habitats, e.g., in freshwater algae, endosymbionts of reef-building corals, deadly parasites, anaerobes in termite guts, or slime molds. Hence, purine biocrystallization is a general and ancestral eukaryotic process likely present in the last eukaryotic common ancestor (LECA) and here we propose two proteins omnipresent in eukaryotes that are likely in charge of their metabolism: hypoxanthine-guanine phosphoribosyl transferase and equilibrative nucleoside transporter. Purine crystalline inclusions are multifunctional structures representing high-capacity and rapid-turnover reserves of nitrogen and optically active elements, e.g., used in light sensing. Thus, we anticipate our work to be a starting point for further studies spanning from cell biology to global ecology, with potential applications in biotechnologies, bio-optics, or in human medicine.

Large landslides cluster at the margin of a deglaciated mountain belt.

Landslides in deglaciated and deglaciating mountains represent a major hazard, but their distribution at the spatial scale of entire mountain belts has rarely been studied. Traditional models of landslide distribution assume that landslides are concentrated in the steepest, wettest, and most tectonically active parts of the orogens, where glaciers reached their greatest thickness. However, based on mapping large landslides (> 0.9 km2) over an unprecedentedly large area of Southern Patagonia (~ 305,000 km2), we show that the distribution of landslides can have the opposite trend. We show that the largest landslides within the limits of the former Patagonian Ice Sheet (PIS) cluster along its eastern margins occupying lower, tectonically less active, and arid part of the Patagonian Andes. In contrast to the heavily glaciated, highest elevations of the mountain range, the peripheral regions have been glaciated only episodically, leaving a larger volume of unstable sedimentary and volcanic rocks that are subject to ongoing slope instability.

Complex causes of landslides after ice sheet retreat: Post-LGM mass movements in the Northern Patagonian Icefield region.

Although the dynamics of individual rock-slope failures above recently shrinking glaciers have received increasing study, less is known about the spatial distribution of landslides in paraglacial settings. Here, we present a landslide inventory for large deglaciated area (~100,000 km2) situated within the Last Glacial Maximum (LGM) limits of the Northern Patagonian Icefield (NPI). Using satellite images and the TanDEM-X digital elevation model, we mapped a total of 15,543 landslides, among which 1006 are deep-seated landslides (DSLs) with area ≥0.01 km2. The distribution of DSLs is highly asymmetric in a W-E transect of the NPI region, with pronounced clustering along the semi-arid eastern front of the Patagonian Andes. The most strongly affected domain is volcanic tablelands overlying weak Miocene sedimentary rocks, but DSLs tend to also cluster along recently deglaciated (i.e. since the end of the 19th century) eastern margin of the NPI. Compared with other high mountain regions, alpine valleys of the Patagonian Andes are affected by DSLs only in <1% of their area, an order of magnitude lower than in other reported deglaciated mountains. The modest incidence of DSLs in the Patagonian Andes is due to dominance of hard granitoid rocks and relatively weak historical seismic activity. We conclude that 1) geological conditions control the distribution of DSLs and their types in the NPI region; 2) paraglacial effects play secondary (although locally important) roles in the origin of DSLs; 3) local clusters of large DSLs originate due to specifics of the post-LGM landscape evolution, involving drawdowns of glacial lakes and incision of rivers into the unconsolidated deposits; and 4) increased abundance of landslides above the recently shrinking margin of the NPI results from the repeated Holocene fluctuations of glacier snouts around the Little Ice Age (LIA) glacier limits and the spatial coincidence of glacial debuttressing effects with the presence of active faults.

Retortamonads from vertebrate hosts share features of anaerobic metabolism and pre-adaptations to parasitism with diplomonads.

Although the mitochondria of extant eukaryotes share a single origin, functionally these organelles diversified to a great extent, reflecting lifestyles of the organisms that host them. In anaerobic protists of the group Metamonada, mitochondria are present in reduced forms (also termed hydrogenosomes or mitosomes) and a complete loss of mitochondrion in Monocercomonoides exilis (Metamonada:Preaxostyla) has also been reported. Within metamonads, retortamonads from the gastrointestinal tract of vertebrates form a sister group to parasitic diplomonads (e.g. Giardia and Spironucleus) and have also been hypothesized to completely lack mitochondria. We obtained transcriptomic data from Retortamonas dobelli and R. caviae and searched for enzymes of the core metabolism as well as mitochondrion- and parasitism-related proteins. Our results indicate that retortamonads have a streamlined metabolism lacking pathways for metabolites they are probably capable of obtaining from prey bacteria or their environment, reminiscent of the biochemical arrangement in other metamonads. Retortamonads were surprisingly found do encode homologs of components of Giardia's remarkable ventral disk, as well as homologs of regulatory NEK kinases and secreted lytic enzymes known for involvement in host colonization by Giardia. These can be considered pre-adaptations of these intestinal microorganisms to parasitism. Furthermore, we found traces of the mitochondrial metabolism represented by iron­sulfur cluster assembly subunits, subunits of mitochondrial translocation and chaperone machinery and, importantly, [FeFe]­hydrogenases and hydrogenase maturases (HydE, HydF and HydG). Altogether, our results strongly suggest that a remnant mitochondrion is still present.

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.1 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,2 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.5 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.