Flexible B12 ecophysiology of Phaeocystis antarctica due to a fusion B12-independent methionine synthase with widespread homologues.

Coastal Antarctic marine ecosystems are significant in carbon cycling because of their intense seasonal phytoplankton blooms. Southern Ocean algae are primarily limited by light and iron (Fe) and can be co-limited by cobalamin (vitamin B12). Micronutrient limitation controls productivity and shapes the composition of blooms which are typically dominated by either diatoms or the haptophyte Phaeocystis antarctica. However, the vitamin requirements and ecophysiology of the keystone species P. antarctica remain poorly characterized. Using cultures, physiological analysis, and comparative omics, we examined the response of P. antarctica to a matrix of Fe-B12 conditions. We show that P. antarctica is not auxotrophic for B12, as previously suggested, and identify mechanisms underlying its B12 response in cultures of predominantly solitary and colonial cells. A combination of proteomics and proteogenomics reveals a B12-independent methionine synthase fusion protein (MetE-fusion) that is expressed under vitamin limitation and interreplaced with the B12-dependent isoform under replete conditions. Database searches return homologues of the MetE-fusion protein in multiple Phaeocystis species and in a wide range of marine microbes, including other photosynthetic eukaryotes with polymorphic life cycles as well as bacterioplankton. Furthermore, we find MetE-fusion homologues expressed in metaproteomic and metatranscriptomic field samples in polar and more geographically widespread regions. As climate change impacts micronutrient availability in the coastal Southern Ocean, our finding that P. antarctica has a flexible B12 metabolism has implications for its relative fitness compared to B12-auxotrophic diatoms and for the detection of B12-stress in a more diverse set of marine microbes.

Molecular forecasting of domoic acid during a pervasive toxic diatom bloom.

In 2015, the largest recorded harmful algal bloom (HAB) occurred in the Northeast Pacific, causing nearly 100 million dollars in damages to fisheries and killing many protected marine mammals. Dominated by the toxic diatom Pseudo-nitzschia australis, this bloom produced high levels of the neurotoxin domoic acid (DA). Through molecular and transcriptional characterization of 52 near-weekly phytoplankton net-tow samples collected at a bloom hotspot in Monterey Bay, California, we identified active transcription of known DA biosynthesis (dab) genes from the three identified toxigenic species, including P. australis as the primary origin of toxicity. Elevated expression of silicon transporters (sit1) during the bloom supports the previously hypothesized role of dissolved silica (Si) exhaustion in contributing to bloom physiology and toxicity. We find that coexpression of the dabA and sit1 genes serves as a robust predictor of DA one week in advance, potentially enabling the forecasting of DA-producing HABs. We additionally present evidence that low levels of iron could have colimited the diatom population along with low Si. Iron limitation represents an overlooked driver of both toxin production and ecological success of the low-iron-adapted Pseudo-nitzschia genus during the 2015 bloom, and increasing pervasiveness of iron limitation may fuel the escalating magnitude and frequency of toxic Pseudo-nitzschia blooms globally. Our results advance understanding of bloom physiology underlying toxin production, bloom prediction, and the impact of global change on toxic blooms.

Mlf mediates proteotoxic response via formation of cellular foci for protein folding and degradation in Giardia.

Myeloid leukemia factor 1 (Mlf1) was identified as a proto-oncoprotein that affects hematopoietic differentiation in humans. However, its cellular function remains elusive, spanning roles from cell cycle regulation to modulation of protein aggregate formation and participation in ciliogenesis. Given that structurally conserved homologs of Mlf1 can be found across the eukaryotic tree of life, we decided to characterize its cellular role underlying this phenotypic pleiotropy. Using a model of the unicellular eukaryote Giardia intestinalis, we demonstrate that its Mlf1 homolog (GiMlf) mainly localizes to two types of cytosolic foci: microtubular structures, where it interacts with Hsp40, and ubiquitin-rich, membraneless compartments, found adjacent to mitochondrion-related organelles known as mitosomes, containing the 26S proteasome regulatory subunit 4. Upon cellular stress, GiMlf either relocates to the affected compartment or disperses across the cytoplasm, subsequently accumulating into enlarged foci during the recovery phase. In vitro assays suggest that GiMlf can be recruited to membranes through its affinity for signaling phospholipids. Importantly, cytosolic foci diminish in the gimlf knockout strain, which exhibits extensive proteomic changes indicative of compromised proteostasis. Consistent with data from other cellular systems, we propose that Mlf acts in the response to proteotoxic stress by mediating the formation of function-specific foci for protein folding and degradation.

New plastids, old proteins: repeated endosymbiotic acquisitions in kareniacean dinoflagellates.

Dinoflagellates are a diverse group of ecologically significant micro-eukaryotes that can serve as a model system for plastid symbiogenesis due to their susceptibility to plastid loss and replacement via serial endosymbiosis. Kareniaceae harbor fucoxanthin-pigmented plastids instead of the ancestral peridinin-pigmented ones and support them with a diverse range of nucleus-encoded plastid-targeted proteins originating from the haptophyte endosymbiont, dinoflagellate host, and/or lateral gene transfers (LGT). Here, we present predicted plastid proteomes from seven distantly related kareniaceans in three genera (Karenia, Karlodinium, and Takayama) and analyze their evolutionary patterns using automated tree building and sorting. We project a relatively limited ( ~ 10%) haptophyte signal pointing towards a shared origin in the family Chrysochromulinaceae. Our data establish significant variations in the functional distributions of these signals, emphasizing the importance of micro-evolutionary processes in shaping the chimeric proteomes. Analysis of plastid genome sequences recontextualizes these results by a striking finding the extant kareniacean plastids are in fact not all of the same origin, as two of the studied species (Karlodinium armiger, Takayama helix) possess plastids from different haptophyte orders than the rest.

Fates of Sec, Tat, and YidC Translocases in Mitochondria and Other Eukaryotic Compartments.

Formation of mitochondria by the conversion of a bacterial endosymbiont was a key moment in the evolution of eukaryotes. It was made possible by outsourcing the endosymbiont's genetic control to the host nucleus, while developing the import machinery for proteins synthesized on cytosolic ribosomes. The original protein export machines of the nascent organelle remained to be repurposed or were completely abandoned. This review follows the evolutionary fates of three prokaryotic inner membrane translocases Sec, Tat, and YidC. Homologs of all three translocases can still be found in current mitochondria, but with different importance for mitochondrial function. Although the mitochondrial YidC homolog, Oxa1, became an omnipresent independent insertase, the other two remained only sporadically present in mitochondria. Only a single substrate is known for the mitochondrial Tat and no function has yet been assigned for the mitochondrial Sec. Finally, this review compares these ancestral mitochondrial proteins with their paralogs operating in the plastids and the endomembrane system.

Changes in Phenotypes and DNA Methylation of In Vitro Aging Sperm in Common Carp Cyprinus carpio.

The purpose of the current study was to analyze phenotypic and functional characteristics of common carp (Cyprinus carpio) spermatozoa during in vitro aging and to investigate whether global DNA methylation is affected by sperm aging. Milt was collected from five individual males, stored in vitro on ice in a refrigerator for up to 96 h post stripping (HPS) and used to fertilize eggs with intervals of 1, 24 and 96 h. Computer-assisted sperm analysis and a S3e Cell Sorter was employed to determine the spermatozoa phenotypic characteristics (motility, velocity, concentration and viability). In addition, pH and osmolality of the seminal fluid and the capacity of the spermatozoa to fertilize, hatching rate and health of the resulting embryos were examined at different aging times. Whole-genome bisulfite sequencing was used to compare the global and gene-specific DNA methylation in fresh and aged spermatozoa. The results demonstrated that spermatozoa aging in common carp significantly affects their performance and thus the success of artificial fertilization. The methylation level at the cytosine-phosphate-guanine (CpG) sites increased significantly with 24 HPS spermatozoa compared to the fresh group at 1 HPS and then decreased significantly at 96 HPS. A more detailed investigation of gene specific differences in the DNA methylation was hindered by incomplete annotation of the C. carpio genome in the public databases.

Metabolic quirks and the colourful history of the Euglena gracilis secondary plastid.

Euglena spp. are phototrophic flagellates with considerable ecological presence and impact. Euglena gracilis harbours secondary green plastids, but an incompletely characterised proteome precludes accurate understanding of both plastid function and evolutionary history. Using subcellular fractionation, an improved sequence database and MS we determined the composition, evolutionary relationships and hence predicted functions of the E. gracilis plastid proteome. We confidently identified 1345 distinct plastid protein groups and found that at least 100 proteins represent horizontal acquisitions from organisms other than green algae or prokaryotes. Metabolic reconstruction confirmed previously studied/predicted enzymes/pathways and provided evidence for multiple unusual features, including uncoupling of carotenoid and phytol metabolism, a limited role in amino acid metabolism, and dual sets of the SUF pathway for FeS cluster assembly, one of which was acquired by lateral gene transfer from Chlamydiae. Plastid paralogues of trafficking-associated proteins potentially mediating fusion of transport vesicles with the outermost plastid membrane were identified, together with derlin-related proteins, potential translocases across the middle membrane, and an extremely simplified TIC complex. The Euglena plastid, as the product of many genomes, combines novel and conserved features of metabolism and transport.

Isolation of plastids and mitochondria from Chromera velia.

MAIN CONCLUSION: We present an easy and effective procedure to purify plastids and mitochondria from Chromera velia. Our method enables downstream analyses of protein and metabolite content of the organelles. Chromerids are alveolate algae that are the closest known phototrophic relatives to apicomplexan parasites such as Plasmodium or Toxoplasma. While genomic and transcriptomic resources for chromerids are in place, tools and experimental conditions for proteomic studies have not been developed yet. Here we describe a rapid and efficient protocol for simultaneous isolation of plastids and mitochondria from the chromerid alga Chromera velia. This procedure involves enzymatic treatment and breakage of cells, followed by differential centrifugation. While plastids sediment in the first centrifugation step, mitochondria remain in the supernatant. Subsequently, plastids can be purified from the crude pellet by centrifugation on a discontinuous 60%/70% sucrose density gradient, while mitochondria can be obtained by centrifugation on a discontinuous 33%/80% Percoll density gradient. Isolated plastids are autofluorescent, and their multi-membrane structure was confirmed by transmission electron microscopy. Fluorescent optical microscopy was used to identify isolated mitochondria stained with MitoTrackerTM green, while their intactness and membrane potential were confirmed by staining with MitoTrackerTM orange CMTMRos. Total proteins were extracted from isolated organellar fractions, and the purity of isolated organelles was confirmed using immunoblotting. Antibodies against the beta subunit of the mitochondrial ATP synthase and the plastid protochlorophyllide oxidoreductase did not cross-react on immunoblots, suggesting that each organellar fraction is free of the residues of the other. The presented protocol represents an essential step for further proteomic, organellar, and cell biological studies of C. velia and can be employed, with minor optimizations, in other thick-walled unicellular algae.

CpSAT-1, a transcribed satellite sequence from the codling moth, Cydia pomonella.

Satellite DNA (satDNA) is a non-coding component of eukaryotic genomes, located mainly in heterochromatic regions. Relevance of satDNA began to emerge with accumulating evidence of its potential yet hardly comprehensible role that it can play in the genome of many organisms. We isolated the first satDNA of the codling moth (Cydia pomonella, Tortricidae, Lepidoptera), a species with holokinetic chromosomes and a single large heterochromatic element, the W chromosome in females. The satDNA, called CpSAT-1, is located on all chromosomes of the complement, although in different amounts. Surprisingly, the satellite is almost missing in the heterochromatic W chromosome. Additionally, we isolated mRNA from all developmental stages (1st-5th instar larva, pupa, adult), both sexes (adult male and female) and several tissues (Malpighian tubules, gut, heart, testes, and ovaries) of the codling moth and showed the CpSAT-1 sequence was transcribed in all tested samples. Using CpSAT-1 specific primers we amplified, cloned and sequenced 40 monomers from cDNA and gDNA, respectively. The sequence analysis revealed a high mutation rate and the presence of potentially functional motifs, mainly in non-conserved regions of the monomers. Both the chromosomal distribution and the sequence analysis suggest that CPSAT-1 has no function in the C. pomonella genome.

Complex Endosymbioses I: From Primary to Complex Plastids, Serial Endosymbiotic Events.

A considerable part of the diversity of eukaryotic phototrophs consists of algae with plastids that evolved from endosymbioses between two eukaryotes. These complex plastids are characterized by a high number of envelope membranes (more than two) and some of them contain a residual nucleus of the endosymbiotic alga called a nucleomorph. Complex plastid-bearing algae are thus chimeric cell assemblies, eukaryotic symbionts living in a eukaryotic host. In contrast, the primary plastids of the Archaeplastida (plants, green algae, red algae, and glaucophytes) possibly evolved from a single endosymbiosis with a cyanobacterium and are surrounded by two membranes. Complex plastids have been acquired several times by unrelated groups of eukaryotic heterotrophic hosts, suggesting that complex plastids are somewhat easier to obtain than primary plastids. Evidence suggests that complex plastids arose twice independently in the green lineage (euglenophytes and chlorarachniophytes) through secondary endosymbiosis, and four times in the red lineage, first through secondary endosymbiosis in cryptophytes, then by higher-order events in stramenopiles, alveolates, and haptophytes. Engulfment of primary and complex plastid-containing algae by eukaryotic hosts (secondary, tertiary, and higher-order endosymbioses) is also responsible for numerous plastid replacements in dinoflagellates. Plastid endosymbiosis is accompanied by massive gene transfer from the endosymbiont to the host nucleus and cell adaptation of both endosymbiotic partners, which is related to the trophic switch to phototrophy and loss of autonomy of the endosymbiont. Such a process is essential for the metabolic integration and division control of the endosymbiont in the host. Although photosynthesis is the main advantage of acquiring plastids, loss of photosynthesis often occurs in algae with complex plastids. This chapter summarizes the essential knowledge of the acquisition, evolution, and function of complex plastids.

Copper in parasitic protists - a hitherto neglected virulence factor.

Copper plays a fundamental role in aerobic metabolism, but its role is double-edged, given its toxicity. Our understanding of copper metabolism in parasites remains rudimentary, despite its significance in virulence. Here we discuss how parasitic protists control copper homeostasis and show the potential key players identified by our bioinformatic analysis.

Mixotrophic growth of a ubiquitous marine diatom.

Diatoms are major players in the global carbon cycle, and their metabolism is affected by ocean conditions. Understanding the impact of changing inorganic nutrients in the oceans on diatoms is crucial, given the changes in global carbon dioxide levels. Here, we present a genome-scale metabolic model (iMK1961) for Cylindrotheca closterium, an in silico resource to understand uncharacterized metabolic functions in this ubiquitous diatom. iMK1961 represents the largest diatom metabolic model to date, comprising 1961 open reading frames and 6718 reactions. With iMK1961, we identified the metabolic response signature to cope with drastic changes in growth conditions. Comparing model predictions with Tara Oceans transcriptomics data unraveled C. closterium's metabolism in situ. Unexpectedly, the diatom only grows photoautotrophically in 21% of the sunlit ocean samples, while the majority of the samples indicate a mixotrophic (71%) or, in some cases, even a heterotrophic (8%) lifestyle in the light. Our findings highlight C. closterium's metabolic flexibility and its potential role in global carbon cycling.

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.

Convergent evolution and horizontal gene transfer in Arctic Ocean microalgae.

Microbial communities in the world ocean are affected strongly by oceanic circulation, creating characteristic marine biomes. The high connectivity of most of the ocean makes it difficult to disentangle selective retention of colonizing genotypes (with traits suited to biome specific conditions) from evolutionary selection, which would act on founder genotypes over time. The Arctic Ocean is exceptional with limited exchange with other oceans and ice covered since the last ice age. To test whether Arctic microalgal lineages evolved apart from algae in the global ocean, we sequenced four lineages of microalgae isolated from Arctic waters and sea ice. Here we show convergent evolution and highlight geographically limited HGT as an ecological adaptive force in the form of PFAM complements and horizontal acquisition of key adaptive genes. Notably, ice-binding proteins were acquired and horizontally transferred among Arctic strains. A comparison with Tara Oceans metagenomes and metatranscriptomes confirmed mostly Arctic distributions of these IBPs. The phylogeny of Arctic-specific genes indicated that these events were independent of bacterial-sourced HGTs in Antarctic Southern Ocean microalgae.

Combinatorial analysis of lupulin gland transcription factors from R2R3Myb, bHLH and WDR families indicates a complex regulation of chs_H1 genes essential for prenylflavonoid biosynthesis in hop (Humulus Lupulus L.).

BACKGROUND: Lupulin glands of hop produce a specific metabolome including hop bitter acids valuable for the brewing process and prenylflavonoids with promising health-beneficial activities. The detailed analysis of the transcription factor (TF)-mediated regulation of the oligofamily of one of the key enzymes, i.e., chalcone synthase CHS_H1 that efficiently catalyzes the production of naringenin chalcone, a direct precursor of prenylflavonoids in hop, constitutes an important part of the dissection of the biosynthetic pathways leading to the accumulation of these compounds. RESULTS: Homologues of flavonoid-regulating TFs HlMyb2 (M2), HlbHLH2 (B2) and HlWDR1 (W1) from hop were cloned using a lupulin gland-specific cDNA library from the hop variety Osvald's 72. Using a "combinatorial" transient GUS expression system it was shown that these unique lupulin-gland-associated TFs significantly activated the promoter (P) of chs_H1 in ternary combinations of B2, W1 and either M2 or the previously characterized HlMyb3 (M3). The promoter activation was strongly dependent on the Myb-P binding box TCCTACC having a core sequence CCWACC positioned on its 5' end region and it seems that the complexity of the promoter plays an important role. M2B2W1-mediated activation significantly exceeded the strength of expression of native chs_H1 gene driven by the 35S promoter of CaMV, while M3B2W1 resulted in 30% of the 35S:chs_H1 expression level, as quantified by real-time PCR. Another newly cloned hop TF, HlMyb7, containing a transcriptional repressor-like motif pdLNLD/ELxiG/S (PDLNLELRIS), was identified as an efficient inhibitor of chs_H1-activating TFs. Comparative analyses of hop and A. thaliana TFs revealed a complex activation of Pchs_H1 and Pchs4 in combinatorial or independent manners. CONCLUSIONS: This study on the sequences and functions of various lupulin gland-specific transcription factors provides insight into the complex character of the regulation of the chs_H1 gene that depends on variable activation by combinations of R2R3Myb, bHLH and WDR TF homologues and inhibition by a Myb repressor.

Biological and molecular analysis of the pathogenic variant C3 of potato spindle tuber viroid (PSTVd) evolved during adaptation to chamomile (Matricaria chamomilla).

Viroid-caused pathogenesis is a specific process dependent on viroid and host genotype(s), and may involve viroid-specific small RNAs (vsRNAs). We describe a new PSTVd variant C3, evolved through sequence adaptation to the host chamomile (Matricaria chamomilla) after biolistic inoculation with PSTVd-KF440-2, which causes extraordinary strong ('lethal') symptoms. The deletion of a single adenine A in the oligoA stretch of the pathogenicity (P) domain appears characteristic of PSTVd-C3. The pathogenicity and the vsRNA pool of PSTVd-C3 were compared to those of lethal variant PSTVd-AS1, from which PSTVd-C3 differs by five mutations located in the P domain. Both lethal viroid variants showed higher stability and lower variation in analyzed vsRNA pools than the mild PSTVd-QFA. PSTVd-C3 and -AS1 caused similar symptoms on chamomile, tomato, and Nicotiana benthamiana, and exhibited similar but species-specific distributions of selected vsRNAs as quantified using TaqMan probes. Both lethal PSTVd variants block biosynthesis of lignin in roots of cultured chamomile and tomato. Four 'expression markers' (TCP3, CIPK, VSF-1, and VPE) were selected from a tomato EST library to quantify their expression upon viroid infection; these markers were strongly downregulated in tomato leaf blades infected by PSTVd-C3- and -AS1 but not by PSTVd-QFA.

The response of Naegleria gruberi to oxidative stress.

Aerobic organisms require oxygen for respiration but must simultaneously cope with oxidative damages inherently linked with this molecule. Unicellular amoeboflagellates of the genus Naegleria, containing both free-living species and opportunistic parasites, thrive in aerobic environments. However, they are also known to maintain typical features of anaerobic organisms. Here, we describe the mechanisms of oxidative damage mitigation in Naegleria gruberi and focus on the molecular characteristics of three noncanonical proteins interacting with oxygen and its derived reactive forms. We show that this protist expresses hemerythrin, protoglobin, and an aerobic-type rubrerythrin, with spectral properties characteristic of the cofactors they bind. We provide evidence that protoglobin and hemerythrin interact with oxygen in vitro and confirm the mitochondrial localization of rubrerythrin by immunolabeling. Our proteomic analysis and immunoblotting following heavy metal treatment revealed upregulation of hemerythrin, while rotenone treatment resulted in an increase in rubrerythrin protein levels together with a vast upregulation of alternative oxidase. Our study provided new insights into the mechanisms employed by N. gruberi to cope with different types of oxidative stress and allowed us to propose specific roles for three unique and understudied proteins: hemerythrin, protoglobin, and rubrerythrin.

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.

Evolutionary and Molecular Aspects of Plastid Endosymbioses.

Plastids are membrane-bound organelles that bestow phototrophic abilities to eukaryotes [...].

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.