Compartmentalization of a glycolytic enzyme in Diplonema, a non-kinetoplastid euglenozoan

Glycosomes are peroxisome-related organelles containing glycolytic enzymes that have been found only in kinetoplastids. We show here that a glycolytic enzyme is compartmentalized in diplonemids, the sister group of kinetoplastids. We found that, similar to kinetoplastid aldolases, the fructose 1,6-bisphosphate aldolase of Diplonema papillatum possesses a type 2-peroxisomal targeting signal. Western blotting showed that this aldolase was present predominantly in the membrane/organellar fraction. Immunofluorescence analysis showed that this aldolase had a scattered distribution in the cytosol, suggesting its compartmentalization. In contrast, orotidine-5'-monophosphate decarboxylase, a non-glycolytic glycosomal enzyme in kinetoplastids, was shown to be a cytosolic enzyme in D. papillatum. Since euglenoids, the earliest diverging branch of Euglenozoa, do not possess glycolytic compartments, these findings suggest that the routing of glycolytic enzymes into peroxisomes may have occurred in a common ancestor of diplonemids and kinetoplastids, followed by diversification of these newly established organelles in each of these euglenozoan lineages.     

Evolutionary analysis of synteny and gene fusion for pyrimidine biosynthetic enzymes in Euglenozoa: an extraordinary gap between kinetoplastids and diplonemids

A unique feature of the genome architecture in the parasitic trypanosomatid protists is large-scale synteny. We addressed the evolutionary trait of synteny in the eukaryotic group, Euglenozoa, which consists of euglenoids (earliest branching), diplonemids, and kinetoplastids (trypanosomatids and bodonids). Synteny of the pyrimidine biosynthetic (pyr) gene cluster, which constitutes part of a large syntenic cluster in trypanosomatids and includes four separate genes (pyr1-pyr4) and one fused gene (pyr6/pyr5 fusion), was conserved in the bodonid, Parabodo caudatus. In the diplonemid, Diplonema papillatum, we identified pyr4 and pyr6 genes. Phylogenetic analyses of pyr4 and pyr6 showed the separate origin of each in kinetoplastids and euglenoids/diplonemids and suggested that kinetoplastids have acquired these genes via lateral gene transfer (LGT). Because replacement of genes by non-orthologs within the syntenic cluster is highly unlikely, we concluded that, after separation of the line leading to diplonemids, the syntenic pyr gene cluster was established in the common ancestor of kinetoplastids, preceded by their acquisition via LGT. Notably, we found that diplonemid pyr6 is a stand-alone gene, inconsistent with both euglenoid pyr5/pyr6 and kinetoplastid pyr6/pyr5 fusions. Our findings provide insights into the evolutionary gaps within Euglenozoa and the evolutionary trait of rearrangement of gene fusion in this lineage.     

Unique mitochondrial genome structure in diplonemids, the sister group of kinetoplastids

Kinetoplastid flagellates are characterized by uniquely massed mitochondrial DNAs (mtDNAs), the kinetoplasts. Kinetoplastids of the trypanosomatid group possess two types of mtDNA molecules: maxicircles bearing protein and mitoribosomal genes and minicircles specifying guide RNAs, which mediate uridine insertion/deletion RNA editing. These circles are interlocked with one another to form dense networks. Whether these peculiar mtDNA features are restricted to kinetoplastids or prevail throughout Euglenozoa (euglenids, diplonemids, and kinetoplastids) is unknown. Here, we describe the mitochondrial genome and the mitochondrial ultrastructure of Diplonema papillatum, a member of the diplonemid flagellates, the sister group of kinetoplastids. Fluorescence and electron microscopy show a single mitochondrion per cell with an ultrastructure atypical for Euglenozoa. In addition, DNA is evenly distributed throughout the organelle rather than compacted. Molecular and electron microscopy studies distinguish numerous 6- and 7-kbp-sized mitochondrial chromosomes of monomeric circular topology and relaxed conformation in vivo. Remarkably, the cox1 gene (and probably other mitochondrial genes) is fragmented, with separate gene pieces encoded on different chromosomes. Generation of the contiguous cox1 mRNA requires trans-splicing, the precise mechanism of which remains to be determined. Taken together, the mitochondrial gene/genome structure of Diplonema is not only different from that of kinetoplastids but unique among eukaryotes as a whole.     

Protein phylogenies robustly resolve the deep-level relationships within Euglenozoa

The deepest-level relationships amongst Euglenozoa remain poorly resolved, despite a rich history of morphological examination and numerous molecular phylogenetic studies of small subunit ribosomal RNA (SSU rRNA) data. We address this question using two nuclear-encoded proteins, the cytosolic isoforms of heat shock protein 90 (hsp90) and heat shock protein 70 (hsp70). For both proteins we examined sequences from the three primary groups within Euglenozoa (euglenids, diplonemids, and kinetoplastids), and from their close relatives, Heterolobosea. Maximum likelihood (ML) and ML distance analyses of these proteins support a close relationship between diplonemids and kinetoplastids to the exclusion of the euglenid Euglena gracilis. In hsp90 and combined protein analyses bootstrap support is very strong and alternative topologies are generally rejected by 'approximately unbiased' (AU) tests. This result is consistent with recent molecular biological and morphological data, but contradicts early structural accounts and many SSU rRNA analyses that favour a closer relationship between diplonemids and euglenids. However, a re-examination of an important SSU rRNA data set highlights the instability of the inferences from this marker. The protein analyses also suggest that bodonids are paraphyletic, with trypanosomatids grouping with 'clade 2' and 'clade 3' bodonids to the exclusion of 'clade 1' bodonids.     

Ribosomal RNA phylogeny of bodonid and diplonemid flagellates and the evolution of euglenozoa

Euglenozoa is a major phylum of excavate protozoa (comprising euglenoids, kinetoplastids, and diplonemids) with highly unusual nuclear, mitochondrial, and chloroplast genomes. To improve understanding of euglenozoan evolution, we sequenced nuclear small-subunit rRNA genes from 34 bodonids (Bodo, Neobodo, Parabodo, Dimastigella-like, Rhynchobodo, Rhynchomonas, and unidentified strains), nine diplonemids (Diplonema, Rhynchopus), and a euglenoid (Entosiphon). Phylogenetic analysis reveals that diplonemids and bodonids are more diverse than previously recognised, but does not clearly establish the branching order of kinetoplastids, euglenoids, and diplonemids. Rhynchopus is holophyletic; parasitic species arose from within free-living species. Kinetoplastea (bodonids and trypanosomatids) are robustly holophyletic and comprise a major clade including all trypanosomatids and most bodonids ('core bodonids') and a very divergent minor one including Ichthyobodo. The root of the major kinetoplastid clade is probably between trypanosomatids and core bodonids. Core bodonids have three distinct subclades. Clade 1 has two distinct Rhynchobodo-like lineages; a lineage comprising Dimastigella and Rhynchomonas; and another including Cruzella and Neobodo. Clade 2 comprises Cryptobia/ Trypanoplasma, Procryptobia, and Parabodo. Clade 3 is an extensive Bodo saltans species complex. Neobodo designis is a vast genetically divergent species complex with mutually exclusive marine and freshwater subclades. Our analysis supports three phagotrophic euglenoid orders: Petalomonadida (holophyletic), Ploeotiida (probably holophyletic), Peranemida (paraphyletic).     

Phylogenetic position of Rhynchopus sp. and Diplonema ambulator as indicated by analyses of euglenozoan small subunit ribosomal DNA

The taxa Rhynchopus Skuja and Diplonema Griessmann were first described as remarkable protists with euglenid affinities. Later on, the placement of Diplonema within the Euglenozoa was confirmed by molecular data. For this study two new sequences were added to the euglenozoan data set. The uncertainly placed Rhynchopus can be identified as a close relative to Diplonema by small subunit ribosomal DNA (SSU rDNA) analysis. The new sequence of Diplonema ambulator is in close relationship to two other Diplonema species. Our molecular analyses clearly support the monophyly of the diplonemids comprising Rhynchopus and Diplonema. Yet the topology at the base of the euglenozoan tree remains unresolved, and especially the monophyly of the euglenids is arguable. SSU rDNA sequence analyses suggest that significantly different GC contents, high mutational saturation in the euglenids, and different evolutionary rates in the euglenozoan clades make it difficult to identify any sister group to the diplonemids.     

Phylogenomic analysis of kinetoplastids supports that trypanosomatids arose from within bodonids

Kinetoplastids are a large group of free-living and parasitic eukaryotic flagellates, including the medically important trypanosomatids (e.g., Trypanosoma and Leishmania) and the widespread free-living and parasitic bodonids. Small subunit rRNA- and conserved protein-based phylogenies support the division of kinetoplastids into five orders (Prokinetoplastida, Neobodonida, Parabodonida, Eubodonida, and Trypanosomatida), but they produce incongruent results regarding their relative branching order, in particular for the position of the Trypanosomatida. In general, small subunit rRNA tends to support their early emergence, whereas protein phylogenies most often support a more recent origin from within bodonids. In order to resolve this question through a phylogenomic approach, we carried out massive parallel sequencing of cDNA from representatives of three bodonid orders (Bodo saltans -Eubodonida-, Procryptobia sorokini -Parabodonida-, and Rhynchomonas nasuta -Neobodonida-). We identified 64 well-conserved proteins shared by these species, four trypanosomatids, and two closely related outgroup species (Euglena gracilis and Diplonema papillatum). Phylogenetic analysis of a concatenated data set yielded a strongly supported tree showing the late emergence of trypanosomatids as a sister group of the Eubodonida. In addition, we identified homologues of proteins involved in trypanosomatid mitochondrial mRNA editing in the three bodonid species, suggesting that editing may be widespread in kinetoplastids. Comparison of expressed sequences from mitochondrial genes showed variability at U positions, in agreement with the existence of editing activity in the three bodonid orders most closely related to trypanosomatids (Neobodonida, Parabodonida, and Eubodonida). Mitochondrial mRNA editing appears to be an ancient phenomenon in kinetoplastids.     

Targeted integration by homologous recombination enables in situ tagging and replacement of genes in the marine microeukaryote Diplonema papillatum

Diplonemids are a group of highly diverse and abundant marine microeukaryotes that belong to the phylum Euglenozoa and form a sister clade to the well-studied, mostly parasitic kinetoplastids. Very little is known about the biology of diplonemids, as few species have been formally described and just one, Diplonema papillatum, has been studied to a decent extent at the molecular level. Following up on our previous results showing stable but random integration of delivered extraneous DNA, we demonstrate here homologous recombination in D. papillatum. Targeting various constructs to the intended position in the nuclear genome was successful when 5′ and 3′ homologous regions longer than 1 kbp were used, achieving N-terminal tagging with mCherry and gene replacement of α- and β-tubulins. For more convenient genetic manipulation, we designed a modular plasmid, pDP002, which bears a protein-A tag and used it to generate and express a C-terminally tagged mitoribosomal protein. Lastly, we developed an improved transformation protocol for broader applicability across laboratories. Our robust methodology allows the replacement, integration as well as endogenous tagging of D. papillatum genes, thus opening the door to functional studies in this species and establishing a basic toolkit for reverse genetics of diplonemids in general.     

Environmental determinants of the distribution of planktonic diplonemids and kinetoplastids in the oceans

We analysed a widely used barcode, the V9 region of the 18S rRNA gene, to study the effect of environmental conditions on the distribution of two related heterotrophic protistan lineages in marine plankton, kinetoplastids and diplonemids. We relied on a major published dataset (Tara Oceans) where samples from the mesopelagic zone were available from just 32 of 123 locations, and both groups are most abundant in this zone. To close sampling gaps and obtain more information from the deeper ocean, we collected 57 new samples targeting especially the mesopelagic zone. We sampled in three geographic regions: the Arctic, two depth transects in the Adriatic Sea, and the anoxic Cariaco Basin. In agreement with previous studies, both protist groups are most abundant and diverse in the mesopelagic zone. In addition to that, we found that their abundance, richness, and community structure also depend on geography, oxygen concentration, salinity, temperature, and other environmental variables reflecting the abundance of algae and nutrients. Both groups studied here demonstrated similar patterns, although some differences were also observed. Kinetoplastids and diplonemids prefer tropical regions and nutrient-rich conditions and avoid high oxygen concentration, high salinity, and high density of algae.     

Early evolution within kinetoplastids (euglenozoa), and the late emergence of trypanosomatids

Many important relationships amongst kinetoplastids, including the position of trypanosomatids, remain uncertain, with limited taxon sampling of markers other than small subunit ribosomal RNA (SSUrRNA). We report gene sequences for cytosolic heat shock proteins 90 and/or 70 (HSP90, HSP70) from the potentially early-diverging kinetoplastids Ichthyobodo necator and Rhynchobodo sp., and from bodonid clades ‘2’ (Parabodonidae) and ‘3’ (Eubodonidae). Some of the new cytosolic HSP70 sequences represent a distinct paralog family (HSP70-B), which is related to yet another paralog known from trypanosomatids (HSP70-C). The (HSP70-B, HSP70-C) clade seemingly diverged before the separation between kinetoplastids and diplonemids. Protein phylogenies support the basal placement of Ichthyobodo within kinetoplastids. Unexpectedly, Rhynchobodo usually forms the next most basal group, separated from the clade ‘1’ bodonids with which it has been allied. Bootstrap support is often weak, but the possibility that Rhynchobodo represents a separate early-diverging lineage within core kinetoplastids deserves further testing. Trypanosomatids always fall remote from the root of kinetoplastids, forming a specific relationship with bodonid clades 2 (and 3), generally with strong bootstrap support. These protein trees with improved taxon sampling provide the best evidence to date for a ‘late’ emergence of trypanosomatids, contradicting recent SSUrRNA-based proposals for a relatively early divergence of this group.     

A freshwater radiation of diplonemids

Diplonemids are considered marine protists and have been reported among the most abundant and diverse eukaryotes in the world oceans. Recently we detected the presence of freshwater diplonemids in Japanese deep freshwater lakes. However, their distribution and abundances in freshwater ecosystems remain unknown. We assessed abundance and diversity of diplonemids from several geographically distant deep freshwater lakes of the world by amplicon-sequencing, shotgun metagenomics and catalysed reporter deposition-fluorescent in situ hybridization (CARD-FISH). We found diplonemids in all the studied lakes, albeit with low abundances and diversity. We assembled long 18S rRNA sequences from freshwater diplonemids and showed that they form a new lineage distinct from the diverse marine clades. Freshwater diplonemids are a sister-group to a marine clade, which are mainly isolates from coastal and bay areas, suggesting a recent habitat transition from marine to freshwater habitats. Images of CARD-FISH targeted freshwater diplonemids suggest they feed on bacteria. Our analyses of 18S rRNA sequences retrieved from single-cell genomes of marine diplonemids show they encode multiple rRNA copies that may be very divergent from each other, suggesting that marine diplonemid abundance and diversity both have been overestimated. These results have wider implications on assessing eukaryotic abundances in natural habitats by using amplicon-sequencing alone.     

Phylogeny and Reclassification of Hemistasia phaeocysticola (Scherffel) Elbrächter & Schnepf, 1996

Hemistasia phaeocysticola is a marine flagellate that preys on diatoms and dinoflagellates among others. Although its morphology and ultrastructure were previously observed and characterized, its phylogenetic position has not been analyzed using molecular sequence data. This flagellate was classified as a kinetoplastid on the basis of the presence of a kinetoplast in the mitochondrion. However, several morphological characteristics similar to those of diplonemids, a sister group of kinetoplastids, have also been noted. Herein, we report that H. phaeocysticola branches within the diplonemid clade in the phylogenetic tree reconstructed by analyzing 18S rRNA gene sequences. Its systematic placement based on this finding is also discussed.     

Distribution and phylogenetic relationships of freshwater Thaumatomonads with a description of the new species Thaumatomonas coloniensis n. sp

The order Thaumatomonadida includes biflagellated heterotrophic flagellates that form filopodia and typically possess siliceous surface scales. We found thaumatomonads to contribute on average about 5%-10% to flagellate abundance in different benthic habitats. A new species of thaumatomonads, Thaumatomonas coloniensis n. sp., is described on the basis of morphological and molecular biological features. This new species was isolated both from groundwater at Appeldorn near Rees (Germany) and from the Rhine River at Cologne (Germany). We have sequenced the small subunit rRNA (ssu rRNA) gene and a fragment of the large subunit rRNA (lsu rRNA) gene (D3-D5 region) from the isolates of the new species, including the first sequence of a representative of the thaumatomonad genus Gyromitus. In agreement with previous studies, the differences in ribosomal genes of different thaumatomonad species are very small. For understanding the phylogenetic relationships of Thaumatomonadida and to explore their sister group relationships, we have created three sequence data sets (ssu rRNA, partial lsu rRNA, concatenated alignment of both) with the same composition of isolates (from Thaumatomonadida, Euglyphida, Cercomonadidae, and Heteromitidae). According to a Kishino-Hasegawa test, Thaumatomonadida evolved within the Cercozoa as a sister taxon to the Heteromitidae. A possibly close relationship to the Euglyphida, recently grouped together with the Thaumatomonadida in the class Imbricatea/Silicofilosea based on the rRNA data sets was not supported by our analyses.     

Diplonemid glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and prokaryote-to-eukaryote lateral gene transfer.

Lateral gene transfer refers to the movement of genetic information from one genome to another, and the integration of that foreign DNA into its new genetic environment. There are currently only a few well-supported cases of prokaryote-to-eukaryote transfer known that do not involve mitochondria or plastids, but it is not clear whether this reflects a lack of such transfer events, or poor sampling of diverse eukaryotes. One gene where this process is apparently active is glyceraldehyde-3-phosphate dehydrogenase (GAPDH), where lateral transfer has been implicated in the origin of euglenoid and kinetoplastid genes. We have characterised GAPDH genes from diplonemids, heterotrophic flagellates that are closely related to kinetoplastids and euglenoids. Two distinct classes of diplonemid GAPDH genes were found in diplonemids, however, neither class is closely related to any other euglenozoan GAPDH. One diplonemid GAPDH is related to the cytosolic gapC of eukaryotes, although not to either euglenoids or kinetoplastids, and the second is related to cyanobacterial and proteobacterial gap3. The bacterial gap3 gene in diplonemids provides one of the most well-supported examples of lateral gene transfer from a bacterium to a eukaryote characterised to date, and may indicate that diplonemids have acquired a novel biochemical capacity through lateral transfer.     

Structure of the Large Subunit rDNA from a Diatom, and Comparison Between Small and Large Subunit Ribosomal RNA for Studying Stramenopile Evolution

The aim of this study was to compare the usefulness of complete small and large subunit rRNA, and a combination of both molecules, for reconstructing stramenopile evolution. To this end, phylogenies from species of which both sequences are known Acre constructed with the neighbor-joining, maximum parsimony, and maximum likelihood methods. Also the use of structural features of the rRNAs was evaluated. The large subunit rRNA from the diatom Skeletonema pseudocostatum was sequenced in order to have a more complete taxon sampling, and a group I intron was identified. Our results indicated that heterokont algae are monophyletic, with diatoms diverging first. However, as the analysis was restricted to a particular data set containing merely six taxa, the outcome has limited value for elucidating stramenopile relationships. On the other hand, this approach permits comparison of the performance of both rRNA molecules without interference from other factors, such as a different species selection for each molecule. For the taxa used, the large subunit rRNA clearly contained more phylogenetic information than the small subunit rRNA. Although this result can definitely not be generalized and depends on the phvlogeny to be studied, in some cases determining complete large subunit rRNA sequences certainly seems worthwhile.