NommPred: Prediction of Mitochondrial and Mitochondrion-Related Organelle Proteins of Nonmodel Organisms

Distribution and Functional Analysis of Isocitrate Dehydrogenases across Kinetoplastids

Isocitrate dehydrogenase is an enzyme converting isocitrate to α-ketoglutarate in the canonical tricarboxylic acid (TCA) cycle. There are three different types of isocitrate dehydrogenase documented in eukaryotes. Our study points out the complex evolutionary history of isocitrate dehydrogenases across kinetoplastids, where the common ancestor of Trypanosomatidae and Bodonidae was equipped with two isoforms of the isocitrate dehydrogenase enzyme: the NADP+-dependent isocitrate dehydrogenase 1 with possibly dual localization in the cytosol and mitochondrion and NADP+-dependent mitochondrial isocitrate dehydrogenase 2. In the extant trypanosomatids, isocitrate dehydrogenase 1 is present only in a few species suggesting that it was lost upon separation of Trypanosoma spp. and replaced by the mainly NADP+-dependent cytosolic isocitrate dehydrogenase 3 of bacterial origin in all the derived lineages. In this study, we experimentally demonstrate that the omnipresent isocitrate dehydrogenase 2 has a dual localization in both mitochondrion and cytosol in at least four species that possess only this isoform. The apparent lack of the NAD+-dependent isocitrate dehydrogenase activity in trypanosomatid mitochondrion provides further support to the existence of the noncanonical TCA cycle across trypanosomatids and the bidirectional activity of isocitrate dehydrogenase 3 when operating with NADP+ cofactor instead of NAD+. This observation can be extended to all 17 species analyzed in this study, except for Leishmania mexicana, which showed only low isocitrate dehydrogenase activity in the cytosol. The variability in isocitrate oxidation capacity among species may reflect the distinct metabolic strategies and needs for reduced cofactors in particular environments.

Encyclopedia of Family A DNA Polymerases Localized in Organelles: Evolutionary Contribution of Bacteria Including the Proto-Mitochondrion

DNA polymerases synthesize DNA from deoxyribonucleotides in a semiconservative manner and serve as the core of DNA replication and repair machinery. In eukaryotic cells, there are 2 genome-containing organelles, mitochondria, and plastids, which were derived from an alphaproteobacterium and a cyanobacterium, respectively. Except for rare cases of genome-lacking mitochondria and plastids, both organelles must be served by nucleus-encoded DNA polymerases that localize and work in them to maintain their genomes. The evolution of organellar DNA polymerases has yet to be fully understood because of 2 unsettled issues. First, the diversity of organellar DNA polymerases has not been elucidated in the full spectrum of eukaryotes. Second, it is unclear when the DNA polymerases that were used originally in the endosymbiotic bacteria giving rise to mitochondria and plastids were discarded, as the organellar DNA polymerases known to date show no phylogenetic affinity to those of the extant alphaproteobacteria or cyanobacteria. In this study, we identified from diverse eukaryotes 134 family A DNA polymerase sequences, which were classified into 10 novel types, and explored their evolutionary origins. The subcellular localizations of selected DNA polymerases were further examined experimentally. The results presented here suggest that the diversity of organellar DNA polymerases has been shaped by multiple transfers of the PolI gene from phylogenetically broad bacteria, and their occurrence in eukaryotes was additionally impacted by secondary plastid endosymbioses. Finally, we propose that the last eukaryotic common ancestor may have possessed 2 mitochondrial DNA polymerases, POP, and a candidate of the direct descendant of the proto-mitochondrial DNA polymerase I, rdxPolA, identified in this study.

Gleaning Euglenozoa-specific DNA polymerases in public single-cell transcriptome data

Multiple genes encoding family A DNA polymerases (famA DNAPs), which are evolutionary relatives of DNA polymerase I (PolI) in bacteria and phages, have been found in eukaryotic genomes, and many of these proteins are used mainly in organelles. Among members of the phylum Euglenozoa, distinct types of famA DNAP, PolIA, PolIBCD+, POP, and eugPolA, have been found. It is intriguing how the suite of famA DNAPs had been established during the evolution of Euglenozoa, but the DNAP data have not been sampled from the taxa that sufficiently represent the diversity of this phylum. In particular, little sequence data were available for basal branching species in Euglenozoa until recently. Thanks to the single-cell transcriptome data from symbiontids and phagotrophic euglenids, we have an opportunity to cover the "hole" in the repertory of famA DNAPs in the deep branches in Euglenozoa. The current study identified 16 new famA DNAP sequences in the transcriptome data from 33 phagotrophic euglenids and two symbiontids, respectively. Based on the new famA DNAP sequences, the updated diversity and evolution of famA DNAPs in Euglenozoa are discussed.

Comparative analysis of mitochondrion-related organelles in anaerobic amoebozoans

Archamoebae comprises free-living or endobiotic amoebiform protists that inhabit anaerobic or microaerophilic environments and possess mitochondrion-related organelles (MROs) adapted to function anaerobically. We compared in silico reconstructed MRO proteomes of eight species (six genera) and found that the common ancestor of Archamoebae possessed very few typical components of the protein translocation machinery, electron transport chain and tricarboxylic acid cycle. On the other hand, it contained a sulphate activation pathway and bacterial iron-sulphur (Fe-S) assembly system of MIS-type. The metabolic capacity of the MROs, however, varies markedly within this clade. The glycine cleavage system is widely conserved among Archamoebae, except in Entamoeba, probably owing to its role in catabolic function or one-carbon metabolism. MRO-based pyruvate metabolism was dispensed within subgroups Entamoebidae and Rhizomastixidae, whereas sulphate activation could have been lost in isolated cases of Rhizomastix libera, Mastigamoeba abducta and Endolimax sp. The MIS (Fe-S) assembly system was duplicated in the common ancestor of Mastigamoebidae and Pelomyxidae, and one of the copies took over Fe-S assembly in their MRO. In Entamoebidae and Rhizomastixidae, we hypothesize that Fe-S cluster assembly in both compartments may be facilitated by dual localization of the single system. We could not find evidence for changes in metabolic functions of the MRO in response to changes in habitat; it appears that such environmental drivers do not strongly affect MRO reduction in this group of eukaryotes.

Evolutionary analysis of cellular reduction and anaerobicity in the hyper-prevalent gut microbe Blastocystis

Blastocystis is the most prevalent microbial eukaryote in the human and animal gut, yet its role as commensal or parasite is still under debate. Blastocystis has clearly undergone evolutionary adaptation to the gut environment and possesses minimal cellular compartmentalization, reduced anaerobic mitochondria, no flagella, and no reported peroxisomes. To address this poorly understood evolutionary transition, we have taken a multi-disciplinary approach to characterize Proteromonas lacertae, the closest canonical stramenopile relative of Blastocystis. Genomic data reveal an abundance of unique genes in P. lacertae but also reductive evolution of the genomic complement in Blastocystis. Comparative genomic analysis sheds light on flagellar evolution, including 37 new candidate components implicated with mastigonemes, the stramenopile morphological hallmark. The P. lacertae membrane-trafficking system (MTS) complement is only slightly more canonical than that of Blastocystis, but notably, we identified that both organisms encode the complete enigmatic endocytic TSET complex, a first for the entire stramenopile lineage. Investigation also details the modulation of mitochondrial composition and metabolism in both P. lacertae and Blastocystis. Unexpectedly, we identify in P. lacertae the most reduced peroxisome-derived organelle reported to date, which leads us to speculate on a mechanism of constraint guiding the dynamics of peroxisome-mitochondrion reductive evolution on the path to anaerobiosis. Overall, these analyses provide a launching point to investigate organellar evolution and reveal in detail the evolutionary path that Blastocystis has taken from a canonical flagellated protist to the hyper-divergent and hyper-prevalent animal and human gut microbe.

Reduced mitochondria provide an essential function for the cytosolic methionine cycle

The loss of mitochondria in oxymonad protists has been associated with the redirection of the essential Fe-S cluster assembly to the cytosol. Yet as our knowledge of diverse free-living protists broadens, the list of functions of their mitochondrial-related organelles (MROs) expands. We revealed another such function in the closest oxymonad relative, Paratrimastix pyriformis, after we solved the proteome of its MRO with high accuracy, using localization of organelle proteins by isotope tagging (LOPIT). The newly assigned enzymes connect to the glycine cleavage system (GCS) and produce folate derivatives with one-carbon units and formate. These are likely to be used by the cytosolic methionine cycle involved in S-adenosyl methionine recycling. The data provide consistency with the presence of the GCS in MROs of free-living species and its absence in most endobionts, which typically lose the methionine cycle and, in the case of oxymonads, the mitochondria.

Evidence for an Independent Hydrogenosome-to-Mitosome Transition in the CL3 Lineage of Fornicates

Fornicata, a lineage of a broader and ancient anaerobic eukaryotic clade Metamonada, contains diverse taxa that are ideally suited for evolutionary studies addressing various fundamental biological questions, such as the evolutionary trajectory of mitochondrion-related organelles (MROs), the transition between free-living and endobiotic lifestyles, and the derivation of alternative genetic codes. To this end, we conducted detailed microscopic and transcriptome analyses in a poorly documented strain of an anaerobic free-living marine flagellate, PCS, in the so-called CL3 fornicate lineage. Fortuitously, we discovered that the original culture contained two morphologically similar and closely related CL3 representatives, which doubles the taxon representation within this lineage. We obtained a monoeukaryotic culture of one of them and formally describe it as a new member of the family Caviomonadidae, Euthynema mutabile gen. et sp. nov. In contrast to previously studied caviomonads, the endobiotic Caviomonas mobilis and Iotanema spirale, E. mutabile possesses an ultrastructurally discernible MRO. We sequenced and assembled the transcriptome of E. mutabile, and by sequence subtraction, obtained transcriptome data from the other CL3 clade representative present in the original PCS culture, denoted PCS-ghost. Transcriptome analyses showed that the reassignment of only one of the UAR stop codons to encode Gln previously reported from I. spirale does not extend to its free-living relatives and is likely due to a unique amino acid substitution in I. spirale’s eRF1 protein domain responsible for termination codon recognition. The backbone fornicate phylogeny was robustly resolved in a phylogenomic analysis, with the CL3 clade amongst the earliest branching lineages. Metabolic and MRO functional reconstructions of CL3 clade members revealed that all three, including I. spirale, encode homologs of key components of the mitochondrial protein import apparatus and the ISC pathway, indicating the presence of a MRO in all of them. In silico evidence indicates that the organelles of E. mutabile and PCS-ghost host ATP and H₂ production, unlike the cryptic MRO of I. spirale. These data suggest that the CL3 clade has experienced a hydrogenosome-to-mitosome transition independent from that previously documented for the lineage leading to Giardia.

The Pet127 protein is a mitochondrial 5'-to-3' exoribonuclease from the PD-(D/E)XK superfamily involved in RNA maturation and intron degradation in yeasts

Pet127 is a mitochondrial protein found in multiple eukaryotic lineages, but absent from several taxa, including plants and animals. Distant homology suggests that it belongs to the divergent PD-(D/E)XK superfamily which includes various nucleases and related proteins. Earlier yeast genetics experiments suggest that it plays a nonessential role in RNA degradation and 5' end processing. Our phylogenetic analysis suggests that it is a primordial eukaryotic invention that was retained in diverse groups, and independently lost several times in the evolution of other organisms. We demonstrate for the first time that the fungal Pet127 protein in vitro is a processive 5'-to-3' exoribonuclease capable of digesting various substrates in a sequence nonspecific manner. Mutations in conserved residues essential in the PD-(D/E)XK superfamily active site abolish the activity of Pet127. Deletion of the PET127 gene in the pathogenic yeast Candida albicans results in a moderate increase in the steady-state levels of several transcripts and in accumulation of unspliced precursors and intronic sequences of three introns. Mutations in the active site residues result in a phenotype identical to that of the deletant, confirming that the exoribonuclease activity is related to the physiological role of the Pet127 protein. Pet127 activity is, however, not essential for maintaining the mitochondrial respiratory activity in C. albicans.

Genome evolution of a nonparasitic secondary heterotroph, the diatom Nitzschia putrida

Secondary loss of photosynthesis is observed across almost all plastid-bearing branches of the eukaryotic tree of life. However, genome-based insights into the transition from a phototroph into a secondary heterotroph have so far only been revealed for parasitic species. Free-living organisms can yield unique insights into the evolutionary consequence of the loss of photosynthesis, as the parasitic lifestyle requires specific adaptations to host environments. Here, we report on the diploid genome of the free-living diatom Nitzschia putrida (35 Mbp), a nonphotosynthetic osmotroph whose photosynthetic relatives contribute ca. 40% of net oceanic primary production. Comparative analyses with photosynthetic diatoms and heterotrophic algae with parasitic lifestyle revealed that a combination of gene loss, the accumulation of genes involved in organic carbon degradation, a unique secretome, and the rapid divergence of conserved gene families involved in cell wall and extracellular metabolism appear to have facilitated the lifestyle of a free-living secondary heterotroph.

Anaerobic derivates of mitochondria and peroxisomes in the free-living amoeba Pelomyxa schiedti revealed by single-cell genomics

Background

Mitochondria and peroxisomes are the two organelles that are most affected during adaptation to microoxic or anoxic environments. Mitochondria are known to transform into anaerobic mitochondria, hydrogenosomes, mitosomes, and various transition stages in between, collectively called mitochondrion-related organelles (MROs), which vary in enzymatic capacity. Anaerobic peroxisomes were identified only recently, and their putatively most conserved function seems to be the metabolism of inositol. The group Archamoebae includes anaerobes bearing both anaerobic peroxisomes and MROs, specifically hydrogenosomes in free-living Mastigamoeba balamuthi and mitosomes in the human pathogen Entamoeba histolytica, while the organelles within the third lineage represented by Pelomyxa remain uncharacterized.

Results

We generated high-quality genome and transcriptome drafts from Pelomyxa schiedti using single-cell omics. These data provided clear evidence for anaerobic derivates of mitochondria and peroxisomes in this species, and corresponding vesicles were tentatively identified in electron micrographs. In silico reconstructed MRO metabolism harbors respiratory complex II, electron-transferring flavoprotein, a partial TCA cycle running presumably in the reductive direction, pyruvate:ferredoxin oxidoreductase, [FeFe]-hydrogenases, a glycine cleavage system, a sulfate activation pathway, and an expanded set of NIF enzymes for iron-sulfur cluster assembly. When expressed in the heterologous system of yeast, some of these candidates localized into mitochondria, supporting their involvement in the MRO metabolism. The putative functions of P. schiedti peroxisomes could be pyridoxal 5′-phosphate biosynthesis, amino acid and carbohydrate metabolism, and hydrolase activities. Unexpectedly, out of 67 predicted peroxisomal enzymes, only four were also reported in M. balamuthi, namely peroxisomal processing peptidase, nudix hydrolase, inositol 2-dehydrogenase, and d-lactate dehydrogenase. Localizations in yeast corroborated peroxisomal functions of the latter two.

Conclusions

This study revealed the presence and partially annotated the function of anaerobic derivates of mitochondria and peroxisomes in P. schiedti using single-cell genomics, localizations in yeast heterologous systems, and transmission electron microscopy. The MRO metabolism resembles that of M. balamuthi and most likely reflects the state in the common ancestor of Archamoebae. The peroxisomal metabolism is strikingly richer in P. schiedti. The presence of myo-inositol 2-dehydrogenase in the predicted peroxisomal proteome corroborates the situation in other Archamoebae, but future experimental evidence is needed to verify additional functions of this organelle.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01247-w.

Highly flexible metabolism of the marine euglenozoan protist Diplonema papillatum

BACKGROUND

The phylum Euglenozoa is a group of flagellated protists comprising the diplonemids, euglenids, symbiontids, and kinetoplastids. The diplonemids are highly abundant and speciose, and recent tools have rendered the best studied representative, Diplonema papillatum, genetically tractable. However, despite the high diversity of diplonemids, their lifestyles, ecological functions, and even primary energy source are mostly unknown.

RESULTS

We designed a metabolic map of D. papillatum cellular bioenergetic pathways based on the alterations of transcriptomic, proteomic, and metabolomic profiles obtained from cells grown under different conditions. Comparative analysis in the nutrient-rich and nutrient-poor media, as well as the absence and presence of oxygen, revealed its capacity for extensive metabolic reprogramming that occurs predominantly on the proteomic rather than the transcriptomic level. D. papillatum is equipped with fundamental metabolic routes such as glycolysis, gluconeogenesis, TCA cycle, pentose phosphate pathway, respiratory complexes, β-oxidation, and synthesis of fatty acids. Gluconeogenesis is uniquely dominant over glycolysis under all surveyed conditions, while the TCA cycle represents an eclectic combination of standard and unusual enzymes.

CONCLUSIONS

The identification of conventional anaerobic enzymes reflects the ability of this protist to survive in low-oxygen environments. Furthermore, its metabolism quickly reacts to restricted carbon availability, suggesting a high metabolic flexibility of diplonemids, which is further reflected in cell morphology and motility, correlating well with their extreme ecological valence.

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.

Barthelonids represent a deep-branching metamonad clade with mitochondrion-related organelles predicted to generate no ATP

We here report the phylogenetic position of barthelonids, small anaerobic flagellates previously examined using light microscopy alone. Barthelona spp. were isolated from geographically distinct regions and we established five laboratory strains. Transcriptomic data generated from one Barthelona strain (PAP020) were used for large-scale, multi-gene phylogenetic (phylogenomic) analyses. Our analyses robustly placed strain PAP020 at the base of the Fornicata clade, indicating that barthelonids represent a deep-branching metamonad clade. Considering the anaerobic/microaerophilic nature of barthelonids and preliminary electron microscopy observations on strain PAP020, we suspected that barthelonids possess functionally and structurally reduced mitochondria (i.e. mitochondrion-related organelles or MROs). The metabolic pathways localized in the MRO of strain PAP020 were predicted based on its transcriptomic data and compared with those in the MROs of fornicates. We here propose that strain PAP020 is incapable of generating ATP in the MRO, as no mitochondrial/MRO enzymes involved in substrate-level phosphorylation were detected. Instead, we detected a putative cytosolic ATP-generating enzyme (acetyl-CoA synthetase), suggesting that strain PAP020 depends on ATP generated in the cytosol. We propose two separate losses of substrate-level phosphorylation from the MRO in the clade containing barthelonids and (other) fornicates.

Catalase and Ascorbate Peroxidase in Euglenozoan Protists

In this work, we studied the biochemical properties and evolutionary histories of catalase (CAT) and ascorbate peroxidase (APX), two central enzymes of reactive oxygen species detoxification, across the highly diverse clade Eugenozoa. This clade encompasses free-living phototrophic and heterotrophic flagellates, as well as obligate parasites of insects, vertebrates, and plants. We present evidence of several independent acquisitions of CAT by horizontal gene transfers and evolutionary novelties associated with the APX presence. We posit that Euglenozoa recruit these detoxifying enzymes for specific molecular tasks, such as photosynthesis in euglenids and membrane-bound peroxidase activity in kinetoplastids and some diplonemids.

Inventory and Evolution of Mitochondrion-localized Family A DNA Polymerases in Euglenozoa

The order Trypanosomatida has been well studied due to its pathogenicity and the unique biology of the mitochondrion. In Trypanosoma brucei, four DNA polymerases, namely PolIA, PolIB, PolIC, and PolID, related to bacterial DNA polymerase I (PolI), were shown to be localized in mitochondria experimentally. These mitochondrion-localized DNA polymerases are phylogenetically distinct from other family A DNA polymerases, such as bacterial PolI, DNA polymerase gamma (Polγ) in human and yeasts, “plant and protist organellar DNA polymerase (POP)” in diverse eukaryotes. However, the diversity of mitochondrion-localized DNA polymerases in Euglenozoa other than Trypanosomatida is poorly understood. In this study, we discovered putative mitochondrion-localized DNA polymerases in broad members of three major classes of Euglenozoa—Kinetoplastea, Diplonemea, and Euglenida—to explore the origin and evolution of trypanosomatid PolIA-D. We unveiled distinct inventories of mitochondrion-localized DNA polymerases in the three classes: (1) PolIA is ubiquitous across the three euglenozoan classes, (2) PolIB, C, and D are restricted in kinetoplastids, (3) new types of mitochondrion-localized DNA polymerases were identified in a prokinetoplastid and diplonemids, and (4) evolutionarily distinct types of POP were found in euglenids. We finally propose scenarios to explain the inventories of mitochondrion-localized DNA polymerases in Kinetoplastea, Diplonemea, and Euglenida.