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

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.     

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.     

Occurrence of multiple, independent gene fusion events for the fifth and sixth enzymes of pyrimidine biosynthesis in different eukaryotic groups

The genes encoding orotate phosphoribosyltransferase (OPRT) and orotidine-5'-monophosphate decarboxylase (OMPDC), the fifth and sixth enzymes in the de novo pyrimidine biosynthetic pathway, are fused as OPRT-OMPDC in most eukaryotic groups. On the other hand, the inversely linked OMPDC-OPRT fusion is present in trypanosomatids, belonging to kinetoplastids together with bodonids in a supergroup, Euglenozoa. Here, we show the presence of OMPDC-OPRT in the bodonid, Bodo caudatus, while OPRT-OMPDC in Euglena gracilis, another euglenozoan species belonging to euglenoids. These results suggest that the OMPDC-OPRT fusion event occurred in a common ancestor of kinetoplastids. Genome sequence database searches further revealed the presence of OMPDC-OPRT in stramenopiles and cyanobacteria. Phylogenetic reconstruction of OPRT and OMPDC rejected statistically the monophyly of the OPRT domains of stramenopile and kinetoplastid OMPDC-OPRT, demonstrating that these gene fusions do not share a common evolutionary origin, despite the identical gene order. Thus, the OMPDC-OPRT fusion is likely to have emerged independently in these eukaryotic groups. Phylogenetic analyses also suggested that cyanobacterial OMPDC-OPRT arose via lateral transfer. We conclude that gene fusion events occur more frequently than previously thought and that lateral gene transfer has made a marked contribution to establishment of the rearranged structure of OPRT and OMPDC genes in eukaryotes.     

Phylogenetic affinities of Diplonema within the Euglenozoa as inferred from the SSU rRNA gene and partial COI protein sequences

In order to shed light on the phylogenetic position of diplonemids within the phylum Euglenozoa, we have sequenced small subunit rRNA (SSU rRNA) genes from Diplonema (syn. Isonema) papillatum and Diplonema sp. We have also analyzed a partial sequence of the mitochondrial gene for cytochrome c oxidase subunit I from D. papillatum. With both markers, the maximum likelihood method favored a closer grouping of diplonemids with kinetoplastids, while the parsimony and distance suggested a closer relationship of diplonemids with euglenoids. In each case, the differences between the best tree and the alternative trees were small. The frequency of codon usage in the partial D. papillatum COI was different from both related groups; however, as is the case in kinetoplastids but not in Euglena, both the non-canonical UGA codon and the canonical UGG codon were used to encode tryptophan in Diplonema.     

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.     

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.     

The Origin of Dihydroorotate Dehydrogenase Genes of Kinetoplastids, with Special Reference to Their Biological Significance and Adaptation to Anaerobic, Parasitic Conditions

Trypanosoma cruzi dihydroorotate dehydrogenase (DHOD), the fourth enzyme of the de novo pyrimidine biosynthetic pathway, is localized in the cytosol and utilizes fumarate as electron acceptor (fumarate reductase activity), while the enzyme from other various eukaryotes is mitochondrial membrane-linked. Here we report that DHOD-knockout T. cruzi did not express the enzyme protein and could not survive even in the presence of pyrimidine nucleosides, substrates for the potentially active salvage pathway, suggesting a vital role of fumarate reductase activity in the regulation of cellular redox balance. Cloning and phylogenetic analysis of euglenozoan DHOD genes showed that the euglenoid Euglena gracilis had a mitochondrial DHOD and that biflagellated bodonids, a sister group of trypanosomatids within kinetoplastids, harbor the cytosolic DHOD. Further, Bodo saliens, a bodonid, had an ACT/DHOD gene fusion encoding aspartate carbamoyltransferase (ACT), the second enzyme of the de novo pyrimidine pathway, and DHOD. This is the first report of this novel gene structure. These results are consistent with suggestions that an ancient common ancestor of Euglenozoa had a mitochondrial DHOD whose descendant exists in E. gracilis and that a common ancestor of kinetoplastids (bodonids and trypanosomatids) subsequently acquired a cytosolic DHOD by horizontal gene transfer. The cytosolic DHOD gene thus acquired may have contributed to adaptation to anaerobiosis in the kinetoplastid lineage and further contributed to the subsequent establishment of parasitism in a trypanosomatid ancestor. Different molecular strategies for anaerobic adaptation in pyrimidine biosynthesis, used by kinetoplastids and by euglenoids, are discussed. Evolutionary implications of the ACT/DHOD gene fusion are also discussed.Sequence availability: The nucleotide sequence data reported here appear in the GenBank, EMBL, and DDBJ databases with the accession numbers AB120414, AB159227, and AB159228 for Euglena gracilis dihydroorotate dehydrogenase (DHOD), Bodo saliens aspartate carbamoyltransferase/dihydroorotate dehydrogenase (ACT/DHOD), and B. caudatus DHOD, respectively.Reviewing Editor: Dr. Patrick Keeling     

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.     

PHYLOGENETIC RELATIONSHIPS OF SELECTED EUGLENOID GENERA BASED ON MORPHOLOGICAL AND MOLECULAR DATA1

The small subunit rRNA (SSU rRNA) coding regions sequenced from the euglenoids Petalomonas cantuscygni, Peranema trichophorum, and Khawkinea quartana were used to assess the phylogenetic relationships of these genera within the Euglenozoa. Phylogenies derived from distance, parsimony, and maximum likelihood methods infer that the euglenoids and kinetoplastids form sister clades within a monophyletic assemblage. Distances representative of closely related lineages separate the genera within the Kinetoplastida, whereas larger distance values separate genera within the euglenoid assemblage. The results of the morphological and molecular studies suggest that phagotrophy arose early in the euglenozoan lineage with the subsequent acquisition of phototrophy, osmotrophy, and parasitism. Phagotrophic euglenoids with a pellicle composed of longitudinal strips appear to have diverged prior to genera with helically arranged strips. This study suggests that the hypothetical ancestor to the Euglenozoa was a phagotroph with two flagella, both containing paraxonemal rods. Furthermore, its basal bodies contained proximal cartwheels, were connected by a prominent fiber, and were anchored with three asymmetrically arranged flagellar roots.     

Phylogeny of phagotrophic euglenids (Euglenozoa) as inferred from hsp90 gene sequences

Molecular phylogenies of euglenids are usually based on ribosomal RNA genes that do not resolve the branching order among the deeper lineages. We addressed deep euglenid phylogeny using the cytosolic form of the heat-shock protein 90 gene (hsp90), which has already been employed with some success in other groups of euglenozoans and eukaryotes in general. Hsp90 sequences were generated from three taxa of euglenids representing different degrees of ultrastructural complexity, namely Petalomonas cantuscygni and wild isolates of Entosiphon sulcatum, and Peranema trichophorum. The hsp90 gene sequence of P. trichophorum contained three short introns (ranging from 27 to 31 bp), two of which had non-canonical borders GG-GG and GG-TG and two 10-bp inverted repeats, suggesting a structure similar to that of the non-canonical introns described in Euglena gracilis. Phylogenetic analyses confirmed a closer relationship between kinetoplastids and diplonemids than to euglenids, and supported previous views regarding the branching order among primarily bacteriovorous, primarily eukaryovorous, and photosynthetic euglenids. The position of P. cantuscygni within Euglenozoa, as well as the relative support for the nodes including it were strongly dependent on outgroup selection. The results were most consistent when the jakobid Reclinomonas americana was used as the outgroup. The most robust phylogenies place P. cantuscygni as the most basal branch within the euglenid clade. However, the presence of a kinetoplast-like mitochondrial inclusion in P. cantuscygni deviates from the currently accepted apomorphy-based definition of the kinetoplastid clade and highlights the necessity of detailed studies addressing the molecular nature of the euglenid and diplonemid mitochondrial genome.     

Phylogenetic analyses of various euglenoid taxa (euglenozoa) based on 18s rdna sequence data

18S rRNA genes (SSU rDNA) of five newly sequenced species were used as molecular markers to infer phylogenetic relationships within the euglenoids. Two members of the order Euglenales ( Lepocinclis ovata Playfair , Phacus similis Christen), two of the order Eutreptiales ( Distigma proteus Ehrenberg, , D. curvata Pringsheim) and Gyropaigne lefévrei Bourelly et Georges of the order Rhabdomonadales were used in parsimony, maximum likelihood, and distance analyses. All trees derived from SSU rRNA data strongly supported the monophyletic origin of the Euglenozoa, with kinetoplastids as sister clade to the euglenoids and Petalomonas cantuscygni Cann et Pennick diverging at the base of the monophyletic euglenoid lineage. The data also supported the theory that phagotrophic euglenoids arose prior to osmotrophs and phototrophs. A lineage of Peranema trichophorum Ehrenberg and all sequenced Euglenales formed a sister clade to the osmotrophs. This suggests that the evolution of phototrophy within the euglenoids radiated from a single event.     

Intrageneric Variability Between the Chloroplast Genomes of Trachelomonas grandis and Trachelomonas volvocina and Phylogenomic Analysis of Phototrophic Euglenoids

The latest studies of chloroplast genomes of phototrophic euglenoids yielded different results according to intrageneric variability such as cluster arrangement or diversity of introns. Although the genera Euglena and Monomorphina in those studies show high syntenic arrangements at the intrageneric level, the two investigated Eutreptiella species comprise low synteny. Furthermore Trachelomonas volvocina show low synteny to the chloroplast genomes of the sister genera Monomorphina aenigmatica, M. parapyrum, Cryptoglena skujae, Euglenaria anabaena, Strombomonas acuminata, all of which were highly syntenic. Consequently, this study aims at the analysis of the cpGenome of Trachelomonas grandis and a comparative examination of T. volvocina to investigate whether the cpGenomes are of such resemblance as could be expected for a genus within the Euglenaceae. Although these analyses resulted in almost identical gene content to other Euglenaceae, the chloroplast genome showed significant novelties: In the rRNA operon, we detected group II introns, not yet found in any other cpGenome of Euglenaceae and a substantially heterogeneous cluster arrangement in the genus Trachelomonas. The phylogenomic analysis with 84 genes of 19 phototrophic euglenoids and 18 cpGenome sequences from Chlorophyta and Streptophyta resulted in a well‐supported cpGenome phylogeny, which is in accordance to former phylogenetic analyses.     

Ultrastructure and molecular phylogeny of Calkinsia aureus: cellular identity of a novel clade of deep-sea euglenozoans with epibiotic bacteria

Background  

The Euglenozoa is a large group of eukaryotic flagellates with diverse modes of nutrition. The group consists of three main subclades - euglenids, kinetoplastids and diplonemids - that have been confirmed with both molecular phylogenetic analyses and a combination of shared ultrastructural characteristics. Several poorly understood lineages of putative euglenozoans live in anoxic environments, such as Calkinsia aureus, and have yet to be characterized at the molecular and ultrastructural levels. Improved understanding of these lineages is expected to shed considerable light onto the ultrastructure of prokaryote-eukaryote symbioses and the associated cellular innovations found within the Euglenozoa and beyond.     

The gene map of the Norway rat (Rattus norvegicus) and comparative mapping with mouse and man

The current status of the rat gene map is presented. Mapping information is now available for a total of 214 loci and the number of mapped genes is increasing steadily. The corresponding number of loci quoted at HGM10 was 128. Genes have been assigned to 20 of the 22 chromosomes in the rat. Some aspects of comparative mapping with mouse and man are also discussed. It was found that there is a good correlation between the morphological homologies detectable in rat and mouse chromosomes, on the one hand, and homology at the gene level on the other. For 10 rat synteny groups all the genes so far mapped are syntenic also in the mouse. For the remaining rat synteny groups it appears that the majority of the genes will be syntenic on specific (homologous) mouse chromosomes, with only a few genes dispersed to other members of the mouse karyotype. Furthermore, the data indicate that mouse chromosome 1 genetically corresponds to two rat chromosomes, viz., 9 and 13, equalizing the difference in chromosome number between the two species. Further mappings will show whether the genetic homology will prove to be as extensive as these preliminary results indicate. As might be expected from evolutionary considerations, rat synteny groups are much more dispersed in the human genome. It is clear, however, that many groups of genes have remained syntenic during the period since man and rat shared a common ancestor. One further point was noted. In two cases groups of genes were syntenic in the mouse but dispersed to two chromosomes in rat and man, whereas in a third case a group of genes was syntenic in the rat but dispersed to two chromosomes in mouse and man. This finding argues in favor of the notion that the original gene groups were on separate ancestral chromosomes, which have fused in one rodent species but remained separate in the other and in man.     

Phylogenomics of T4 cyanophages: lateral gene transfer in the 'core' and origins of host genes

The last two decades have revealed that phages (viruses that infect bacteria) are abundant and play fundamental roles in the Earth System, with the T4-like myoviruses (herein T4-like phages) emerging as a dominant 'signal' in wild populations. Here we examine 27 T4-like phage genomes, with a focus on 17 that infect ocean picocyanobacteria (cyanophages), to evaluate lateral gene transfer (LGT) in this group. First, we establish a reference tree by evaluating concatenated core gene supertrees and whole genome gene content trees. Next, we evaluate what fraction of these 'core genes' shared by all 17 cyanophages appear prone to LGT. Most (47 out of 57 core genes) were vertically transferred as inferred from tree tests and genomic synteny. Of those 10 core genes that failed the tree tests, the bulk (8 of 10) remain syntenic in the genomes with only a few (3 of the 10) having identifiable signatures of mobile elements. Notably, only one of these 10 is shared not only by the 17 cyanophages, but also by all 27 T4-like phages (thymidylate synthase); its evolutionary history suggests cyanophages may be the origin of these genes to Prochlorococcus. Next, we examined intragenic recombination among the core genes and found that it did occur, even among these core genes, but that the rate was significantly higher between closely related phages, perhaps reducing any detectable LGT signal and leading to taxon cohesion. Finally, among 18 auxiliary metabolic genes (AMGs, a.k.a. 'host' genes), we found that half originated from their immediate hosts, in some cases multiple times (e.g. psbA, psbD, pstS), while the remaining have less clear evolutionary origins ranging from cyanobacteria (4 genes) or microbes (5 genes), with particular diversity among viral TalC and Hsp20 sequences. Together, these findings highlight the patterns and limits of vertical evolution, as well as the ecological and evolutionary roles of LGT in shaping T4-like phage genomes.     

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.     

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.