Unexpectedly Streamlined Mitochondrial Genome of the Euglenozoan Euglena gracilis

Dinoflagellate chloroplasts as a model for extreme genome reduction and fragmentation in organelles - The COCOA principle for gene retention

The genomes of peridinin-containing dinoflagellate chloroplasts have a very unusual organisation. These genomes are highly fragmented and greatly reduced, with most of the usual complement of chloroplast genes relocated to the nucleus. Dinoflagellate chloroplasts highlight evolutionary changes that are found to varying extents in a number of other organelle genomes. These include the chloroplast genome of the green alga Boodlea and other Cladophorales, and the mitochondrial genomes of blood-sucking and chewing lice, the parasitic plant Rhopalocnemis phalloides, the red alga Rhodosorus marinus and other members of the Stylonematophyceae, diplonemid flagellates, and some Cnidaria. Consideration of the coding content of the remnant chloroplast genomes indicates that organelles may preferentially retain genes for proteins important in initiating assembly of complexes, and the same is largely true for mitochondria. We propose a new principle, of CO-location for COntrol of Assembly (COCOA), indicating the importance of retaining these genes in the organelle. This adds to, but does not invalidate, the existing hypotheses of the multisubunit completion principle, CO-location for Redox Regulation (CORR) and Control by Epistasy of Synthesis (CES).

Unique features of conventional and nonconventional introns in Euglena gracilis

BACKGROUND

Nuclear introns in Euglenida have been understudied. This study aimed to investigate nuclear introns in Euglenida by identifying a large number of introns in Euglena gracilis (E. gracilis), including cis-spliced conventional and nonconventional introns, as well as trans-spliced outrons. We also examined the sequence characteristics of these introns.

RESULTS

A total of 28,337 introns and 11,921 outrons were identified. Conventional and nonconventional introns have distinct splice site features; the former harbour canonical GT/C-AG splice sites, whereas the latter are capable of forming structured motifs with their terminal sequences. We observed that short introns had a preference for canonical GT-AG introns. Notably, conventional introns and outrons in E. gracilis exhibited a distinct cytidine-rich polypyrimidine tract, in contrast to the thymidine-rich tracts observed in other organisms. Furthermore, the SL-RNAs in E. gracilis, as well as in other trans-splicing species, can form a recently discovered motif called the extended U6/5' ss duplex with the respective U6s. We also describe a novel type of alternative splicing pattern in E. gracilis. The tandem repeat sequences of introns in this protist were determined, and their contents were comparable to those in humans.

CONCLUSIONS

Our findings highlight the unique features of E. gracilis introns and provide insights into the splicing mechanism of these introns, as well as the genomics and evolution of Euglenida.

Circular extrachromosomal DNA in Euglena gracilis under normal and stress conditions

Extrachromosomal circular DNA (eccDNA) enhances genomic plasticity, augmenting its coding and regulatory potential. Advances in high-throughput sequencing have enabled the investigation of these structural variants. Although eccDNAs have been investigated in numerous taxa, they remained understudied in euglenids. Therefore, we examined eccDNAs predicted from Illumina sequencing data of Euglena gracilis Z SAG 1224-5/25, grown under optimal photoperiod and exposed to UV irradiation. We identified approximately 1000 unique eccDNA candidates, about 20% of which were shared across conditions. We also observed a significant enrichment of mitochondrially encoded eccDNA in the UV-irradiated sample. Furthermore, we found that the heterogeneity of eccDNA was reduced in UV-exposed samples compared to cells that were grown in optimal conditions. Hence, eccDNA appears to play a role in the response to oxidative stress in Euglena, as it does in other studied organisms. In addition to contributing to the understanding of Euglena genomes, our results contribute to the validation of bioinformatics pipelines on a large, non-model genome.

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.

Single-Cell Genomics Reveals the Divergent Mitochondrial Genomes of Retaria (Foraminifera and Radiolaria)

Mitochondria originated from an ancient bacterial endosymbiont that underwent reductive evolution by gene loss and endosymbiont gene transfer to the nuclear genome. The diversity of mitochondrial genomes published to date has revealed that gene loss and transfer processes are ongoing in many lineages. Most well-studied eukaryotic lineages are represented in mitochondrial genome databases, except for the superphylum Retaria-the lineage comprising Foraminifera and Radiolaria. Using single-cell approaches, we determined two complete mitochondrial genomes of Foraminifera and two nearly complete mitochondrial genomes of radiolarians. We report the complete coding content of an additional 14 foram species. We show that foraminiferan and radiolarian mitochondrial genomes contain a nearly fully overlapping but reduced mitochondrial gene complement compared to other sequenced rhizarians. In contrast to animals and fungi, many protists encode a diverse set of proteins on their mitochondrial genomes, including several ribosomal genes; however, some aerobic eukaryotic lineages (euglenids, myzozoans, and chlamydomonas-like algae) have reduced mitochondrial gene content and lack all ribosomal genes. Similar to these reduced outliers, we show that retarian mitochondrial genomes lack ribosomal protein and tRNA genes, contain truncated and divergent small and large rRNA genes, and contain only 14 or 15 protein-coding genes, including nad1, -3, -4, -4L, -5, and -7, cob, cox1, -2, and -3, and atp1, -6, and -9, with forams and radiolarians additionally carrying nad2 and nad6, respectively. In radiolarian mitogenomes, a noncanonical genetic code was identified in which all three stop codons encode amino acids. Collectively, these results add to our understanding of mitochondrial genome evolution and fill in one of the last major gaps in mitochondrial sequence databases. IMPORTANCE We present the reduced mitochondrial genomes of Retaria, the rhizarian lineage comprising the phyla Foraminifera and Radiolaria. By applying single-cell genomic approaches, we found that foraminiferan and radiolarian mitochondrial genomes contain an overlapping but reduced mitochondrial gene complement compared to other sequenced rhizarians. An alternative genetic code was identified in radiolarian mitogenomes in which all three stop codons encode amino acids. Collectively, these results shed light on the divergent nature of the mitochondrial genomes from an ecologically important group, warranting further questions into the biological underpinnings of gene content variability and genetic code variation between mitochondrial genomes.

Metabolic network reconstruction of Euglena gracilis: Current state, challenges, and applications

A metabolic model, representing all biochemical reactions in a cell, is a prerequisite for several approaches in systems biology used to explore the metabolic phenotype of an organism. Despite the use of Euglena in diverse industrial applications and as a biological model, there is limited understanding of its metabolic network capacity. The unavailability of the completed genome data and the highly complex evolution of Euglena are significant obstacles to the reconstruction and analysis of its genome-scale metabolic model. In this mini-review, we discuss the current state and challenges of metabolic network reconstruction in Euglena gracilis. We have collated and present the available relevant data for the metabolic network reconstruction of E. gracilis, which could be used to improve the quality of the metabolic model of E. gracilis. Furthermore, we deliver the potential applications of the model in metabolic engineering. Altogether, it is supposed that this mini-review would facilitate the investigation of metabolic networks in Euglena and further lay out a direction for model-assisted metabolic engineering.

Microbe Profile: Euglena gracilis: photogenic, flexible and hardy

Euglena gracilis is a unicellular photosynthetic eukaryotic flagellate of the Discoba supergroup, which also encompasses Kinetoplastida and Diplonema. Plastids have green algal origin and are secondarily acquired. The nuclear genome is extremely large and many genes suggest multiple endosymbiotic/gene transfer events, i.e. derivation from prokaryotes of various lineages. E. gracilis is remarkably robust and can proliferate in environments contaminated with heavy metals and acids. Extraordinary metabolic plasticity and a mixotrophic lifestyle confers an ability to thrive in a broad range of environments, as well as facilitating production of many novel metabolites, making Euglena of considerable biotechnological importance.

Discrimination of Euglena gracilis strains Z and bacillaris by MALDI-TOF MS

AIMS

Euglena gracilis is used as model organism for various microbiological, molecular biological and biotechnological studies. Its most studied wild type strains are Z and bacillaris, but their discrimination by standard molecular methods is difficult. Therefore, we decided to test the suitability of MALDI-TOF MS (matrix assisted laser desorption/ionization - time of flight mass spectrometry) for identification of E. gracilis and for discrimination of these two strains possessing functional chloroplasts. MALDI-TOF MS profiling was also tested for two white (non-photosynthetic) stable E. gracilis mutant strains W_gm ZOflL and W₁₀ BSmL.

METHODS AND RESULTS

We have successfully obtained main spectrum profiles (MSPs) of E. gracilis strains Z, SAG 1224-5/25 and bacillaris, SAG 1224-5/15 using protein extraction procedure. Subsequent MALDI-TOF MS profiling of a number of tested samples and the comparison of the obtained protein profiles with our in-house database including MSPs of both strains has revealed that these two strains can be easily distinguished by MALDI-TOF MS based on score values over two in most cases. This method has also confirmed the ancestry of white mutant strains W_gm ZOflL and W₁₀ BSmL, originally derived from strains Z and bacillaris, respectively.

CONCLUSIONS

MALDI-TOF MS is suitable, accurate and rapid method for discrimination of E. gracilis strains.

SIGNIFICANCE AND IMPACT OF STUDY

These results can have broad practical implications for laboratories cultivating various strains of euglenids, and they can be applied for their discrimination by MALDI-TOF MS.

Heterotrophic euglenid Rhabdomonas costata resembles its phototrophic relatives in many aspects of molecular and cell biology

Euglenids represent a group of protists with diverse modes of feeding. To date, only a partial genomic sequence of Euglena gracilis and transcriptomes of several phototrophic and secondarily osmotrophic species are available, while primarily heterotrophic euglenids are seriously undersampled. In this work, we begin to fill this gap by presenting genomic and transcriptomic drafts of a primary osmotroph, Rhabdomonas costata. The current genomic assembly length of 100 Mbp is 14× smaller than that of E. gracilis. Despite being too fragmented for comprehensive gene prediction it provided fragments of the mitochondrial genome and comparison of the transcriptomic and genomic data revealed features of its introns, including several candidates for nonconventional types. A set of 39,456 putative R. costata proteins was predicted from the transcriptome. Annotation of the mitochondrial core metabolism provides the first data on the facultatively anaerobic mitochondrion of R. costata, which in most respects resembles the mitochondrion of E. gracilis with a certain level of streamlining. R. costata can synthetise thiamine by enzymes of heterogenous provenances and haem by a mitochondrial-cytoplasmic C4 pathway with enzymes orthologous to those found in E. gracilis. The low percentage of green algae-affiliated genes supports the ancestrally osmotrophic status of this species.

De Novo Transcriptome Meta-Assembly of the Mixotrophic Freshwater Microalga Euglena gracilis

Euglena gracilis is a well-known photosynthetic microeukaryote considered as the product of a secondary endosymbiosis between a green alga and a phagotrophic unicellular belonging to the same eukaryotic phylum as the parasitic trypanosomatids. As its nuclear genome has proven difficult to sequence, reliable transcriptomes are important for functional studies. In this work, we assembled a new consensus transcriptome by combining sequencing reads from five independent studies. Based on a detailed comparison with two previously released transcriptomes, our consensus transcriptome appears to be the most complete so far. Remapping the reads on it allowed us to compare the expression of the transcripts across multiple culture conditions at once and to infer a functionally annotated network of co-expressed genes. Although the emergence of meaningful gene clusters indicates that some biological signal lies in gene expression levels, our analyses confirm that gene regulation in euglenozoans is not primarily controlled at the transcriptional level. Regarding the origin of E. gracilis, we observe a heavily mixed gene ancestry, as previously reported, and rule out sequence contamination as a possible explanation for these observations. Instead, they indicate that this complex alga has evolved through a convoluted process involving much more than two partners.

Single-cell genomics unveils a canonical origin of the diverse mitochondrial genomes of euglenozoans

Background

The supergroup Euglenozoa unites heterotrophic flagellates from three major clades, kinetoplastids, diplonemids, and euglenids, each of which exhibits extremely divergent mitochondrial characteristics. Mitochondrial genomes (mtDNAs) of euglenids comprise multiple linear chromosomes carrying single genes, whereas mitochondrial chromosomes are circular non-catenated in diplonemids, but circular and catenated in kinetoplastids. In diplonemids and kinetoplastids, mitochondrial mRNAs require extensive and diverse editing and/or trans-splicing to produce mature transcripts. All known euglenozoan mtDNAs exhibit extremely short mitochondrial small (rns) and large (rnl) subunit rRNA genes, and absence of tRNA genes. How these features evolved from an ancestral bacteria-like circular mitochondrial genome remains unanswered.

Results

We sequenced and assembled 20 euglenozoan single-cell amplified genomes (SAGs). In our phylogenetic and phylogenomic analyses, three SAGs were placed within kinetoplastids, 14 within diplonemids, one (EU2) within euglenids, and two SAGs with nearly identical small subunit rRNA gene (18S) sequences (EU17/18) branched as either a basal lineage of euglenids, or as a sister to all euglenozoans. Near-complete mitochondrial genomes were identified in EU2 and EU17/18. Surprisingly, both EU2 and EU17/18 mitochondrial contigs contained multiple genes and one tRNA gene. Furthermore, EU17/18 mtDNA possessed several features unique among euglenozoans including full-length rns and rnl genes, six mitoribosomal genes, and nad11, all likely on a single chromosome.

Conclusions

Our data strongly suggest that EU17/18 is an early-branching euglenozoan with numerous ancestral mitochondrial features. Collectively these data contribute to untangling the early evolution of euglenozoan mitochondria.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01035-y.

Euglenozoa: taxonomy, diversity and ecology, symbioses and viruses

Euglenozoa is a species-rich group of protists, which have extremely diverse lifestyles and a range of features that distinguish them from other eukaryotes. They are composed of free-living and parasitic kinetoplastids, mostly free-living diplonemids, heterotrophic and photosynthetic euglenids, as well as deep-sea symbiontids. Although they form a well-supported monophyletic group, these morphologically rather distinct groups are almost never treated together in a comparative manner, as attempted here. We present an updated taxonomy, complemented by photos of representative species, with notes on diversity, distribution and biology of euglenozoans. For kinetoplastids, we propose a significantly modified taxonomy that reflects the latest findings. Finally, we summarize what is known about viruses infecting euglenozoans, as well as their relationships with ecto- and endosymbiotic bacteria.

Phage Origin of Mitochondrion-Localized Family A DNA Polymerases in Kinetoplastids and Diplonemids

Mitochondria retain their own genomes as other bacterial endosymbiont-derived organelles. Nevertheless, no protein for DNA replication and repair is encoded in any mitochondrial genomes (mtDNAs) assessed to date, suggesting that the nucleus primarily governs the maintenance of mtDNA. As the proteins of diverse evolutionary origins occupy a large proportion of the current mitochondrial proteomes, we anticipate finding the same evolutionary trend in the nucleus-encoded machinery for mtDNA maintenance. Indeed, none of the DNA polymerases (DNAPs) in the mitochondrial endosymbiont, a putative α-proteobacterium, seemingly had been inherited by their descendants (mitochondria), as none of the known types of mitochondrion-localized DNAP showed a specific affinity to the α-proteobacterial DNAPs. Nevertheless, we currently have no concrete idea of how and when the known types of mitochondrion-localized DNAPs emerged. We here explored the origins of mitochondrion-localized DNAPs after the improvement of the samplings of DNAPs from bacteria and phages/viruses. Past studies have revealed that a set of mitochondrion-localized DNAPs in kinetoplastids and diplonemids, namely PolIB, PolIC, PolID, PolI-Perk1/2, and PolI-dipl (henceforth designated collectively as “PolIBCD+”) have emerged from a single DNAP. In this study, we recovered an intimate connection between PolIBCD+ and the DNAPs found in a particular group of phages. Thus, the common ancestor of kinetoplastids and diplonemids most likely converted a laterally acquired phage DNAP into a mitochondrion-localized DNAP that was ancestral to PolIBCD+. The phage origin of PolIBCD+ hints at a potentially large contribution of proteins acquired via nonvertical processes to the machinery for mtDNA maintenance in kinetoplastids and diplonemids.

The Lipid Composition of Euglena gracilis Middle Plastid Membrane Resembles That of Primary Plastid Envelopes

Euglena gracilis is a photosynthetic flagellate possessing chlorophyte-derived secondary plastids that are enclosed by only three enveloping membranes, unlike most secondary plastids, which are surrounded by four membranes. It has generally been assumed that the two innermost E. gracilis plastid envelopes originated from the primary plastid, while the outermost is of eukaryotic origin. It was suggested that nucleus-encoded plastid proteins pass through the middle and innermost plastid envelopes of E. gracilis by machinery homologous to the translocons of outer and inner chloroplast membranes, respectively. Although recent genomic, transcriptomic, and proteomic data proved the presence of a reduced form of the translocon of inner membrane, they failed to identify any outer-membrane translocon homologs, which raised the question of the origin of E. gracilis's middle plastid envelope. Here, we compared the lipid composition of whole cells of the pigmented E. gracilis strain Z and two bleached mutants that lack detectable plastid structures, W10BSmL and WgmZOflL We determined the lipid composition of E. gracilis strain Z mitochondria and plastids, and of plastid subfractions (thylakoids and envelopes), using HPLC high-resolution tandem mass spectrometry, thin-layer chromatography, and gas chromatography-flame ionization detection analytical techniques. Phosphoglycerolipids are the main structural lipids in mitochondria, while glycosyldiacylglycerols are the major structural lipids of plastids and also predominate in extracts of whole mixotrophic cells. Glycosyldiacylglycerols were detected in both bleached mutants, indicating that mutant cells retain some plastid remnants. Additionally, we discuss the origin of the E. gracilis middle plastid envelope based on the lipid composition of envelope fraction.

RNA Editing in Mitochondria and Plastids: Weird and Widespread

Though widespread, RNA editing is rare, except in endosymbiotic organelles. A combination of higher mutation rates, relaxation of energetic constraints, and high genetic drift is found within plastids and mitochondria and is conducive for evolution and expansion of editing processes, possibly starting as repair mechanisms. To illustrate this, we present an exhaustive phylogenetic overview of editing types.

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).

Reductionist Pathways for Parasitism in Euglenozoans? Expanded Datasets Provide New Insights

The unicellular trypanosomatids belong to the phylum Euglenozoa and all known species are obligate parasites. Distinct lineages infect plants, invertebrates, and vertebrates, including humans. Genome data for marine diplonemids, together with freshwater euglenids and free-living kinetoplastids, the closest known nonparasitic relatives to trypanosomatids, recently became available. Robust phylogenetic reconstructions across Euglenozoa are now possible and place the results of parasite-focused studies into an evolutionary context. Here we discuss recent advances in identifying the factors shaping the evolution of Euglenozoa, focusing on ancestral features generally considered parasite-specific. Remarkably, most of these predate the transition(s) to parasitism, suggesting that the presence of certain preconditions makes a significant lifestyle change more likely.

Molecular tools and applications of Euglena gracilis: From biorefineries to bioremediation

Euglena gracilis is a promising source of commercially important metabolites such as vitamins, wax esters, paramylon, and amino acids. However, the molecular tools available to create improved Euglena strains are limited compared to other microorganisms that are currently exploited in the biotechnology industry. The complex poly-endosymbiotic nature of the Euglena genome is a major bottleneck for obtaining a complete genome sequence and thus represents a notable shortcoming in gaining molecular information of this organism. Therefore, the studies and applications have been more focused on using the wild-type strain or its variants and optimizing the nutrient composition and cultivation conditions to enhance the production of biomass and valuable metabolites. In addition to producing metabolites, the E. gracilis biorefinery concept also provides means for the production of biofuels and biogas as well as residual biomass for the remediation of industrial and municipal wastewater. Using Euglena for bioremediation of environments contaminated with heavy metals is of special interest due to the strong ability of the organism to accumulate and sequester these compounds. The published draft genome and transcriptome will serve as a basis for further molecular studies of Euglena and provide a guide for the engineering of metabolic pathways of relevance for the already established as well as novel applications.

A Uniquely Complex Mitochondrial Proteome from Euglena gracilis

Euglena gracilis is a metabolically flexible, photosynthetic, and adaptable free-living protist of considerable environmental importance and biotechnological value. By label-free liquid chromatography tandem mass spectrometry, a total of 1,786 proteins were identified from the E. gracilis purified mitochondria, representing one of the largest mitochondrial proteomes so far described. Despite this apparent complexity, protein machinery responsible for the extensive RNA editing, splicing, and processing in the sister clades diplonemids and kinetoplastids is absent. This strongly suggests that the complex mechanisms of mitochondrial gene expression in diplonemids and kinetoplastids occurred late in euglenozoan evolution, arising independently. By contrast, the alternative oxidase pathway and numerous ribosomal subunits presumed to be specific for parasitic trypanosomes are present in E. gracilis. We investigated the evolution of unexplored protein families, including import complexes, cristae formation proteins, and translation termination factors, as well as canonical and unique metabolic pathways. We additionally compare this mitoproteome with the transcriptome of Eutreptiella gymnastica, illuminating conserved features of Euglenida mitochondria as well as those exclusive to E. gracilis. This is the first mitochondrial proteome of a free-living protist from the Excavata and one of few available for protists as a whole. This study alters our views of the evolution of the mitochondrion and indicates early emergence of complexity within euglenozoan mitochondria, independent of parasitism.

Genetic Engineering Strategies for Euglena gracilis and Its Industrial Contribution to Sustainable Development Goals: A Review

The sustainable development goals (SDGs) adopted at the 2015 United Nations Summit are globally applicable goals designed to help countries realize a sustainable future. To achieve these SDGs, it is necessary to utilize renewable biological resources. In recent years, bioeconomy has been an attractive concept for achieving the SDGs. Microalgae are one of the biological resources that show promise in realizing the "5F"s (food, fiber, feed, fertilizer, and fuel). Among the microalgae, Euglena gracilis has the potential for achieving the "5F"s strategy owing to its unique features, such as production of paramylon, that are lacking in other microalgae. E. gracilis has already been produced on an industrial scale for use as an ingredient in functional foods and cosmetics. In recent years, genetic engineering methods for breeding E. gracilis have been researched and developed to achieve higher yields. In this article, we summarize how microalgae contribute toward achieving the SDGs. We focus on the contribution of E. gracilis to the bioeconomy, including its advantages in industrial use as well as its unique characteristics. In addition, we review genetic engineering-related research trends centered on E. gracilis, including a complete nuclear genome determination project, genome editing technology using the CRISPR-Cas9 system, and the development of a screening method for selecting useful strains. In particular, genome editing in E. gracilis could be a breakthrough for molecular breeding of industrially useful strains because of its high efficiency.

The highly efficient NDH-dependent photosystem I cyclic electron flow pathway in the marine angiosperm Zostera marina

Zostera marina, a fully submerged marine angiosperm with a unique evolutionary history associated with its terrestrial origin, has distinct photochemical characteristics caused by its oxygen-evolving complex (OEC) being prone to deactivation in visible light. Based on the present phylogenetic analysis, the chloroplast NADPH dehydrogenase-like (NDH) complex was found to be completed in of Z. marina, unlike other marine plants, suggesting its crucial role. Thus, the responses of electron transport to irradiation were investigated through multiple chlorophyll fluorescence techniques and Western blot analysis. Moreover, the respective contribution of the two photosystem I cyclic electron flow (PSI-CEF) pathways to the generation of trans-thylakoid proton gradient (∆pH) was also examined using inhibitors. The contributions of the two PSI-CEF pathways to ∆pH were similar; furthermore, there was a trade-off between the two pathways under excess irradiation: the PGR5/L1-dependent PSI-CEF decreased gradually following its activation during the initial illumination, while NDH-dependent PSI-CEF was activated gradually with exposure duration. OEC inactivation was continuously under excess irradiation, which exhibits a positive linear correlation with the activation of NDH-dependent PSI-CEF. We suggest that PGR5/L1-dependent PSI-CEF was preferentially activated to handle the excess electron caused by the operation of OEC during the initial illumination. Subsequently, the increasing OEC inactivation with exposure duration resulted in a deficit of electrons. Limited electrons from PSI might preferentially synthesize NADPH, which could support the function of NDH-dependent PSI-CEF to generate ∆pH and ATP via reducing ferredoxin, thereby maintaining OEC stability.

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.

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.

Structure of a mitochondrial ATP synthase with bound native cardiolipin

The mitochondrial ATP synthase fuels eukaryotic cells with chemical energy. Here we report the cryo-EM structure of a divergent ATP synthase dimer from mitochondria of Euglena gracilis, a member of the phylum Euglenozoa that also includes human parasites. It features 29 different subunits, 8 of which are newly identified. The membrane region was determined to 2.8 Å resolution, enabling the identification of 37 associated lipids, including 25 cardiolipins, which provides insight into protein-lipid interactions and their functional roles. The rotor-stator interface comprises four membrane-embedded horizontal helices, including a distinct subunit a. The dimer interface is formed entirely by phylum-specific components, and a peripherally associated subcomplex contributes to the membrane curvature. The central and peripheral stalks directly interact with each other. Last, the ATPase inhibitory factor 1 (IF₁) binds in a mode that is different from human, but conserved in Trypanosomatids.

Every living thing uses the energy-rich molecule called adenosine triphosphate, or ATP, as fuel. It is the universal molecular currency for transferring energy. Cells trade it, mitochondria make it, and the energy extracted from it is used to drive chemical reactions, transport molecules across cell membranes, energize nerve impulses and contract muscles.

ATP synthase is the enzyme that makes ATP molecules. It is a multi-part complex that straddles the inner membrane of mitochondria, the energy factories in cells. The enzyme complex interacts with fatty molecules in the mitochondrial inner membrane, creating a curvature that is required to produce ATP more efficiently. The mitochondrial ATP synthase has been studied in many different organisms, including yeast, algae, plants, pigs, cows and humans. These studies show that most of these ATP synthases are similar to each other, but obtaining a high resolution structure has been a challenge.

Some single-cell organisms have unusual ATP synthases, which provide clues about how the enzyme evolved in pursuit of the most energy efficient arrangement. One such organism is the photosynthetic Euglena gracilis, which is closely related to the human parasites that cause sleeping sickness and Chagas disease.

Now, Mü̈hleip et al. have extracted ATP synthase from E. gracilis and reconstructed its structure using electron cryo-microscopy. The high resolution of this reconstruction allowed for the first time to examine the fatty molecules associated with ATP synthase, called cardiolipins. This is important, because cardiolipins are thought to modulate the rotating motor of the enzyme and affect how the complex sits in the membrane.

The analysis revealed that the ATP synthase in E. gracilis has 29 different protein subunits, 13 of which are only found in organisms of the same family. Some of the newly discovered subunits are glued together by fatty molecules and extend into the surrounding mitochondrial membrane. This distinctive structure suggests an adaptation which likely evolved independently in E. gracilis for efficiency.

These results represent an important advance in the field, and provide direct evidence for the functional roles of cardiolipin. This information will be used to reconstruct the evolution of this mighty molecule and to further study the roles of cardiolipin in energy conversion. Moreover, the analysis identified similarities between the ATP synthase in E. gracilis and human parasites, which could provide new therapeutic targets in disease-causing parasites.

A single-cell genome reveals diplonemid-like ancestry of kinetoplastid mitochondrial gene structure

Euglenozoa comprises euglenids, kinetoplastids, and diplonemids, with each group exhibiting different and highly unusual mitochondrial genome organizations. Although they are sister groups, kinetoplastids and diplonemids have very distinct mitochondrial genome architectures, requiring widespread insertion/deletion RNA editing and extensive trans-splicing, respectively, in order to generate functional transcripts. The evolutionary history by which these differing processes arose remains unclear. Using single-cell genomics, followed by small sub unit ribosomal DNA and multigene phylogenies, we identified an isolated marine cell that branches on phylogenetic trees as a sister to known kinetoplastids. Analysis of single-cell amplified genomic material identified multiple mitochondrial genome contigs. These revealed a gene architecture resembling that of diplonemid mitochondria, with small fragments of genes encoded out of order and or on different contigs, indicating that these genes require extensive trans-splicing. Conversely, no requirement for kinetoplastid-like insertion/deletion RNA-editing was detected. Additionally, while we identified some proteins so far only found in kinetoplastids, we could not unequivocally identify mitochondrial RNA editing proteins. These data invite the hypothesis that extensive genome fragmentation and trans-splicing were the ancestral states for the kinetoplastid-diplonemid clade but were lost during the kinetoplastid radiation. This study demonstrates that single-cell approaches can successfully retrieve lineages that represent important new branches on the tree of life, and thus can illuminate major evolutionary and functional transitions in eukaryotes.

This article is part of a discussion meeting issue ‘Single cell ecology’.

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.

Transcriptome, proteome and draft genome of Euglena gracilis

BACKGROUND

Photosynthetic euglenids are major contributors to fresh water ecosystems. Euglena gracilis in particular has noted metabolic flexibility, reflected by an ability to thrive in a range of harsh environments. E. gracilis has been a popular model organism and of considerable biotechnological interest, but the absence of a gene catalogue has hampered both basic research and translational efforts.

RESULTS

We report a detailed transcriptome and partial genome for E. gracilis Z1. The nuclear genome is estimated to be around 500 Mb in size, and the transcriptome encodes over 36,000 proteins and the genome possesses less than 1% coding sequence. Annotation of coding sequences indicates a highly sophisticated endomembrane system, RNA processing mechanisms and nuclear genome contributions from several photosynthetic lineages. Multiple gene families, including likely signal transduction components, have been massively expanded. Alterations in protein abundance are controlled post-transcriptionally between light and dark conditions, surprisingly similar to trypanosomatids.

CONCLUSIONS

Our data provide evidence that a range of photosynthetic eukaryotes contributed to the Euglena nuclear genome, evidence in support of the 'shopping bag' hypothesis for plastid acquisition. We also suggest that euglenids possess unique regulatory mechanisms for achieving extreme adaptability, through mechanisms of paralog expansion and gene acquisition.

RNA-Seq employing a novel rRNA depletion strategy reveals a rich repertoire of snoRNAs in Euglena gracilis including box C/D and Ψ-guide RNAs targeting the modification of rRNA extremities

Previous mRNA transcriptome studies of Euglena gracilis have shown that this organism possesses a large and diverse complement of protein coding genes; however, the study of non-coding RNA classes has been limited. The natural extensive fragmentation of the E. gracilis large subunit ribosomal RNA presents additional barriers to the identification of non-coding RNAs as size-selected small RNA libraries will be dominated by rRNA sequences. In this study we have developed a strategy to significantly reduce rRNA amplification prior to RNA-Seq analysis thereby producing a ncRNA library allowing for the identification of many new E. gracilis small RNAs. Library analysis reveals 113 unique new small nucleolar (sno) RNAs and a large collection of snoRNA isoforms, as well as the first significant collection of nuclear tRNAs in this organism. A 3' end AGAUGN consensus motif and conserved structural features can now be defined for E. gracilis pseudouridine guide RNAs. snoRNAs of both classes were identified that target modification of the 3' extremities of rRNAs utilizing predicted base-pairing interactions with internally transcribed spacers (ITS), providing insight into the timing of steps in rRNA maturation. Cumulatively, this represents the most comprehensive analysis of small ncRNAs in Euglena gracilis to date.

Massive mitochondrial DNA content in diplonemid and kinetoplastid protists

The mitochondrial DNA of diplonemid and kinetoplastid protists is known for its suite of bizarre features, including the presence of concatenated circular molecules, extensive trans‐splicing and various forms of RNA editing. Here we report on the existence of another remarkable characteristic: hyper‐inflated DNA content. We estimated the total amount of mitochondrial DNA in four kinetoplastid species (Trypanosoma brucei, Trypanoplasma borreli, Cryptobia helicis, and Perkinsela sp.) and the diplonemid Diplonema papillatum. Staining with 4′,6‐diamidino‐2‐phenylindole and RedDot1 followed by color deconvolution and quantification revealed massive inflation in the total amount of DNA in their organelles. This was further confirmed by electron microscopy. The most extreme case is the ∼260 Mbp of DNA in the mitochondrion of Diplonema, which greatly exceeds that in its nucleus; this is, to our knowledge, the largest amount of DNA described in any organelle. Perkinsela sp. has a total mitochondrial DNA content ~6.6× greater than its nuclear genome. This mass of DNA occupies most of the volume of the Perkinsela cell, despite the fact that it contains only six protein‐coding genes. Why so much DNA? We propose that these bloated mitochondrial DNAs accumulated by a ratchet‐like process. Despite their excessive nature, the synthesis and maintenance of these mtDNAs must incur a relatively low cost, considering that diplonemids are one of the most ubiquitous and speciose protist groups in the ocean. © 2018 The Authors. IUBMB Life published by Wiley Periodicals, Inc. on behalf of International Union of Biochemistry and Molecular Biology., 70(12):1267–1274, 2018

The Mitochondrion of Euglena gracilis

In the presence of oxygen, Euglena gracilis mitochondria function much like mammalian mitochondria. Under anaerobiosis, E. gracilis mitochondria perform a malonyl-CoA independent synthesis of fatty acids leading to accumulation of wax esters, which serve as the sink for electrons stemming from glycolytic ATP synthesis and pyruvate oxidation. Some components (enzymes and cofactors) of Euglena's anaerobic energy metabolism are found among the anaerobic mitochondria of invertebrates, others are found among hydrogenosomes, the H₂-producing anaerobic mitochondria of protists.

Mitochondrial Genome Evolution and a Novel RNA Editing System in Deep-Branching Heteroloboseids

Discoba (Excavata) is an evolutionarily important group of eukaryotes that includes Jakobida, with the most bacterial-like mitochondrial genomes known, and Euglenozoa, many of which have extensively fragmented mitochondrial genomes. However, little is known about the mitochondrial genomes of Heterolobosea, the third main group of Discoba. Here, we studied two heteroloboseids—an undescribed amoeba “BB2” and Pharyngomonas kirbyi. Phylogenomic analysis revealed that they form a clade that is a sister group to all other Heterolobosea. We characterized the mitochondrial genomes of BB2 and P. kirbyi, which encoded 44 and 48 putative protein-coding genes respectively. Their gene contents were similar to that of Naegleria. In BB2, mitochondrially encoded RNAs were heavily edited, with ∼500 mononucleotide insertion events, mostly guanosines. These insertions always have the same identity as an adjacent nucleotide. Editing occurs in all ribosomal RNAs and protein-coding transcripts except one, and half of the transfer RNAs. Analysis of Illumina deep-sequencing data suggested that this RNA editing is very accurate and efficient, and most likely co-transcriptional. The dissimilarity of this editing process to other RNA editing phenomena in discobids, as well as its apparent absence in P. kirbyi, suggest that this remarkably extensive system of insertional editing evolved independently in the BB2 lineage, after its divergence from the P. kirbyi lineage.

Biochemistry and Physiology of Heavy Metal Resistance and Accumulation in Euglena

Free-living microorganisms may become suitable models for removal of heavy metals from polluted water bodies, sediments, and soils by using and enhancing their metal accumulating abilities. The available research data indicate that protists of the genus Euglena are a highly promising group of microorganisms to be used in bio-remediation of heavy metal-polluted aerobic and anaerobic acidic aquatic environments. This chapter analyzes the variety of biochemical mechanisms evolved in E. gracilis to resist, accumulate and remove heavy metals from the environment, being the most relevant those involving (1) adsorption to the external cell pellicle; (2) intracellular binding by glutathione and glutathione polymers, and their further compartmentalization as heavy metal-complexes into chloroplasts and mitochondria; (3) polyphosphate biosynthesis; and (4) secretion of organic acids. The available data at the transcriptional, kinetic and metabolic levels on these metabolic/cellular processes are herein reviewed and analyzed to provide mechanistic basis for developing genetically engineered Euglena cells that may have a greater removal and accumulating capacity for bioremediation and recycling of heavy metals.

Euglena gracilis Genome and Transcriptome: Organelles, Nuclear Genome Assembly Strategies and Initial Features

Euglena gracilis is a major component of the aquatic ecosystem and together with closely related species, is ubiquitous worldwide. Euglenoids are an important group of protists, possessing a secondarily acquired plastid and are relatives to the Kinetoplastidae, which themselves have global impact as disease agents. To understand the biology of E. gracilis, as well as to provide further insight into the evolution and origins of the Kinetoplastidae, we embarked on sequencing the nuclear genome; the plastid and mitochondrial genomes are already in the public domain. Earlier studies suggested an extensive nuclear DNA content, with likely a high degree of repetitive sequence, together with significant extrachromosomal elements. To produce a list of coding sequences we have combined transcriptome data from both published and new sources, as well as embarked on de novo sequencing using a combination of 454, Illumina paired end libraries and long PacBio reads. Preliminary analysis suggests a surprisingly large genome approaching 2 Gbp, with a highly fragmented architecture and extensive repeat composition. Over 80% of the RNAseq reads from E. gracilis maps to the assembled genome sequence, which is comparable with the well assembled genomes of T. brucei and T. cruzi. In order to achieve this level of assembly we employed multiple informatics pipelines, which are discussed here. Finally, as a preliminary view of the genome architecture, we discuss the tubulin and calmodulin genes, which highlight potential novel splicing mechanisms.

Hyper-eccentric structural genes in the mitochondrial genome of the algal parasite Hemistasia phaeocysticola

Diplonemid mitochondria are considered to have very eccentric structural genes. Coding regions of individual diplonemid mitochondrial genes are fragmented into small pieces and found on different circular DNAs. Short RNAs transcribed from each DNA molecule mature through a unique RNA maturation process involving assembly and three types of RNA editing (i.e., U insertion and A-to-I and C-to-U substitutions), although the molecular mechanism(s) of RNA maturation and the evolutionary history of these eccentric structural genes still remain to be understood. Since the gene fragmentation pattern is generally conserved among the diplonemid species studied to date, it was considered that their structural complexity has plateaued and further gene fragmentation could not occur. Here, we show the mitochondrial gene structure of Hemistasia phaeocysticola, which was recently identified as a member of a novel lineage in diplonemids, by comparison of the mitochondrial DNA sequences with cDNA sequences synthesized from mature mRNA. The genes of H. phaeocysticola are fragmented much more finely than those of other diplonemids studied to date. Furthermore, in addition to all known types of RNA editing, it is suggested that a novel processing step (i.e., secondary RNA insertion) is involved in the RNA maturation in the mitochondria of H. phaeocysticola. Our findings demonstrate the tremendous plasticity of mitochondrial gene structures.

A modified single-cell electroporation method for molecule delivery into a motile protist, Euglena gracilis

Single-cell transfection is a powerful technique for delivering chemicals, drugs, or probes into arbitrary, specific single cells. This technique is especially important when the analysis of molecular function and cellular behavior in individual microscopic organisms such as protists requires the precise identification of the target cell, as fluorescence labeling of bulk populations makes tracking of individual motile protists virtually impossible. Herein, we have modified current single-cell electroporation techniques for delivering fluorescent markers into single Euglena gracilis, a motile photosynthetic microalga. Single-cell electroporation introduced molecules into individual living E. gracilis cells after a negative pressure was applied through a syringe connected to the micropipette to the target cell. The new method achieves high transfection efficiency and viability after electroporation. With the new technique, we successfully introduced a variety of molecules such as GFP, Alexa Fluor 488, and exciton-controlled hybridization-sensitive fluorescent oligonucleotide (ECHO) RNA probes into individual motile E. gracilis cells. We demonstrate imaging of endogenous mRNA in living E. gracilis without interfering with their physiological functions, such as swimming or division, over an extended period of time. Thus the modified single-cell electroporation technique is suitable for delivering versatile functional molecules into individual motile protists.

Genes in Hiding

Unrecognizable genes are an unsettling problem in genomics. Here, we survey the various types of cryptic genes and the corresponding deciphering strategies employed by cells. Encryption that renders genes substantially different from homologs in other species includes sequence substitution, insertion, deletion, fragmentation plus scrambling, and invasion by mobile genetic elements. Cells decode cryptic genes at the DNA, RNA or protein level. We will focus on a recently discovered case of unparalleled encryption involving massive gene fragmentation and nucleotide deletions and substitutions, occurring in the mitochondrial genome of a poorly understood protist group, the diplonemids. This example illustrates that comprehensive gene detection requires not only auxiliary sequence information - transcriptome and proteome data - but also knowledge about a cell's deciphering arsenal.

From simple to supercomplex: mitochondrial genomes of euglenozoan protists

Mitochondria are double membrane organelles of endosymbiotic origin, best known for constituting the centre of energetics of a eukaryotic cell. They contain their own mitochondrial genome, which as a consequence of gradual reduction during evolution typically contains less than two dozens of genes. In this review, we highlight the extremely diverse architecture of mitochondrial genomes and mechanisms of gene expression between the three sister groups constituting the phylum Euglenozoa - Euglenida, Diplonemea and Kinetoplastea. The earliest diverging euglenids possess a simplified mitochondrial genome and a conventional gene expression, whereas both are highly complex in the two other groups. The expression of their mitochondrial-encoded proteins requires extensive post-transcriptional modifications guided by complex protein machineries and multiple small RNA molecules. Moreover, the least studied diplonemids, which have been recently discovered as a highly abundant component of the world ocean plankton, possess one of the most complicated mitochondrial genome organisations known to date.

From simple to supercomplex: mitochondrial genomes of euglenozoan protists

Mitochondria are double membrane organelles of endosymbiotic origin, best known for constituting the centre of energetics of a eukaryotic cell. They contain their own mitochondrial genome, which as a consequence of gradual reduction during evolution typically contains less than two dozens of genes. In this review, we highlight the extremely diverse architecture of mitochondrial genomes and mechanisms of gene expression between the three sister groups constituting the phylum Euglenozoa - Euglenida, Diplonemea and Kinetoplastea. The earliest diverging euglenids possess a simplified mitochondrial genome and a conventional gene expression, whereas both are highly complex in the two other groups. The expression of their mitochondrial-encoded proteins requires extensive post-transcriptional modifications guided by complex protein machineries and multiple small RNA molecules. Moreover, the least studied diplonemids, which have been recently discovered as a highly abundant component of the world ocean plankton, possess one of the most complicated mitochondrial genome organisations known to date.

Aerobic mitochondria of parasitic protists: Diverse genomes and complex functions

In this review the main features of the mitochondria of aerobic parasitic protists are discussed. While the best characterized organelles are by far those of kinetoplastid flagellates and Plasmodium, we also consider amoebae Naegleria and Acanthamoeba, a ciliate Ichthyophthirius and related lineages. The simplistic view of the mitochondrion as just a power house of the cell has already been abandoned in multicellular organisms and available data indicate that this also does not apply for protists. We discuss in more details the following mitochondrial features: genomes, post-transcriptional processing, translation, biogenesis of iron-sulfur complexes, heme metabolism and the electron transport chain. Substantial differences in all these core mitochondrial features between lineages are compatible with the view that aerobic protists harbor organelles that are more complex and flexible than previously appreciated.