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

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.     

The evolutionary history of kinetoplastids and their kinetoplasts

Despite extensive phylogenetic analysis of small subunit ribosomal RNA (SSUrRNA) genes, the deep-level relationships among kinetoplastids remain poorly understood, limiting our grasp of their evolutionary history, especially the origins of their bizarre mitochondrial genome organizations. In this study we examine the SSUrRNA data in the light of a new marker--cytoplasmic heat shock protein 90 (hsp90) sequences. Our phylogenetic analyses divide kinetoplastids into four main clades. Clades 1-3 include the various bodonid kinetoplastids. Trypanosomatids comprise the fourth clade. SSUrRNA analyses give vastly different and poorly supported positions for the root of the kinetoplastid tree, depending on the out-group and analysis method. This is probably due to the extraordinary length of the branch between kinetoplastids and any out-group. In contrast, almost all hsp90 analyses place the root between clade 1 (including Dimastigella, Rhynchomonas, several Bodo spp., and probably Rhynchobodo) and all other kinetoplastids. Maximum likelihood and maximum likelihood distance analyses of hsp90 protein and second codon-position nucleotides place trypanosomatids adjacent to Bodo saltans and Bodo cf. uncinatus (clade 3), as (weakly) do SSUrRNA analyses. Hsp90 first codon- plus second codon-position nucleotide analyses return a slightly different topology. We show that this may be an artifact caused, in part, by the different evolutionary behavior of first- and second-codon positions. This study provides the most robust evidence to date that trypanosomatids are descended from within bodonids and that B. saltans is a close relative of trypanosomatids. A total reevaluation of the high-level systematics within kinetoplastids is needed. We confirm that the interlocking network organization of kinetoplast DNA seen in trypanosomatids is a derived condition within kinetoplastids but suggest that open-conformation minicircles may have arisen early in kinetoplastid evolution. Further understanding of the evolution of kinetoplast structure and RNA editing is hampered by a paucity of data from basal (i.e., clade 1) bodonids.     

Phylogeny of the kinetoplastida: taxonomic problems and insights into the evolution of parasitism

To further investigate phylogeny of kinetoplastid protozoa, the sequences of small subunit (18S) ribosomal RNA of nine bodonid isolates and ten isolates of insect trypanosomatids have been determined. The root of the kinetoplastid tree was attached to the branch of Bodo designis and/or Cruzella marina. The suborder Trypanosomatina appeared as a monophyletic group, while the suborder Bodonina was paraphyletic. Among bodonid lineages, parasitic organisms were intermingled with free-living ones, implying multiple transitions to parasitism and supporting the 'vertebrate-first hypothesis'. The tree indicated that the genera Cryptobia and Bodo are artificial taxa. Separation of fish cryptobias and Trypanoplasma borreli as different genera was not supported. In trypanosomatids, the genera Leptomonas and Blastocrithidia were polyphyletic, similar to the genera Herpetomonas and Crithidia and in contrast to the monophyletic genera Trypanosoma and Phytomonas. This analysis has shown that the morphological classification of kinetoplastids does not in general reflect their genetic affinities and needs a revision.     

The evolution and diversity of kinetoplastid flagellates

Five years ago, little was known about kinetoplastid evolution. Recent improvements in the taxon sampling for nuclear rRNA genes and several protein markers have transformed this understanding. Parasitism evolved at least four times in kinetoplastids. Obligate parasitic trypanosomatids are a relatively 'derived' group within kinetoplastids; their closest relative is likely to be the free-living Bodo saltans, and the ancestral trypanosomatids were probably parasites of insects. Although subject to recent controversy, trypanosomes (genus Trypanosoma) probably constitute a monophyletic group. Several unusual features of trypanosomatid genomes (e.g. trans-splicing, mitochondrial RNA editing and intron poverty) are common in kinetoplastids and pre-date the adoption of parasitism. The framework of relationships is becoming robust enough for real comparative approaches to be used to understand kinetoplastid biology.     

基于核糖体小亚基的波豆类鞭毛虫分子系统发育研究

波豆类鞭毛虫(动基体目)是一类非常重要的原生动物,结构独特,分布广泛.但是,这一类群的系统发育关系尚存有很多争议.为了更好的了解该类群的系统发育关系,作者分离、纯化并培养了Bodo designis DH,测定了它的SSU rRNA序列.根据该序列和GenBank中的相关序列,用最大简约法和邻接法分别构建了基于全序列和保守区序列的系统树.结果如下:1)波豆属是多系发育的;2)波豆类鞭毛虫和锥虫类鞭毛虫的系统发育关系仍有待于深入研究;3)虽然所得的系统树都显示出一种基本的双歧式树型结构,但采用不同方法所构建的系统树之间有明显的差异.这些似乎表明,就波豆类鞭毛虫而言,SSU rRNA不是一个令人满意的系统发育标记.     

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.     

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.     

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.     

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     

Kinetoplastid phylogenetics, with special reference to the evolution of parasitic trypanosomes

To fully understand the evolutionary history of parasitic kinetoplastids and to understand the context within which the evolution of each parasite group has developed, an understanding not just of the parasites, but of all kinetoplastids is required. Accordingly, this paper provides an overview of kinetoplastid evolution and systematics, including coverage of the proposal by Moreira et al. (2004) to divide kinetoplasts into Prokinetoplastina (Ichthyobodo and Perkinsiella) and Metakinetoplastina (other bodonids and trypanosomatids). The implications of such a revision, with regard to correctly identifying outgroup taxa for studies of evolution within taxa of medical importance, are addressed, together with a more detailed review of the evolution and origins of the trypanosomes in the light of new phylogenies, new approaches and revisions in kinetoplastid systematics.     

RNA editing and mitochondrial genomic organization in the cryptobiid kinetoplastid protozoan Trypanoplasma borreli.

The bodonids and cryptobiids represent an early diverged sister group to the trypanosomatids among the kinetoplastid protozoa. The trypanosome type of uridine insertion-deletion RNA editing was found to occur in the cryptobiid fish parasite Trypanoplasma borreli. A pan-edited ribosomal protein, S12, and a novel 3'- and 5'-edited cytochrome b, in addition to an unedited cytochrome oxidase III gene and an apparently unedited 12S rRNA gene, were found in a 6-kb fragment of the 80- to 90-kb mitochondrial genome. The gene order differs from that in trypanosomatids, as does the organization of putative guide RNA genes; guide RNA-like molecules are transcribed from tandemly repeated 1-kb sequences organized in 200- and 170-kb molecules instead of minicircles. The presence of pan-editing in this lineage is consistent with an ancient evolutionary origin of this process.     

The evolutionary expansion of the trypanosomatid flagellates

The trypanosomatids combine a relatively uniform morphology with ability to parasitise a very diverse range of hosts including animals, plants and other protists. Along with their sister family, the biflagellate bodonids, they are set apart from other eukaryotes by distinctive organisational features, such as the kinetoplast-mitochondrion and RNA editing, isolation of glycolysis enzymes in the glycosome, use of the flagellar pocket for molecular traffic into and out of the cell, a unique method of generating cortical microtubules, and bizarre nuclear organisation. These features testify to the antiquity and isolation of the kinetoplast-bearing flagellates (Kinetoplastida). Molecular sequencing techniques (especially small subunit ribosomal RNA gene sequencing) are now radically reshaping previous ideas on the phylogeny of these organisms. The idea that the monogenetic (MG) trypanosomatids gave rise to the digenetic (DG) genera is losing ground to a view that, after the bodonids, the African trypanosomes (DG) represent the most ancient lineage, followed by Trypanosoma cruzi (DG), then Blastocrithidia (MG), Herpetomonas (MG) and Phytomonas (DG), with Leptomonas (MG), Crithidia (MG), Leishmania (DG) and Endotrypanum (DG) forming the crown of the evolutionary tree. Vast genetic distances (12% divergence) separate T. brucei and T. cruzi, while the Leishmania species are separated by very short distances (less than 1% divergence). These phylogenetic conclusions are supported by studies on RNA editing and on the nature of the parasite surface. The trypanosomatids seem to be able to adapt with ease their energy metabolism to the availability of substrates and oxygen, and this may give them the ability to institute new life cycles if host behaviour patterns allow. Sexual processes, though present in at least some trypanosomatids, may have played only a minor part in generating diversity during trypanosomatid evolution. On the other hand, the development of altruistic behaviour on the part of some life cycle stages may be a hitherto unconsidered way of maximising fitness in this group. It is concluded that, owing to organisational constraints, the trypanosomatids can undergo substantial molecular variation while registering very little in the way of morphological change.     

Mitochondrial minicircles in the free-living bodonid Bodo saltans contain two gRNA gene cassettes and are not found in large networks

In trypanosomatids, the majority of the guide (g) RNAs that provide the information for U-insertion/deletion RNA editing are encoded by minicircles that are catenated into large networks. In contrast, in the distantly related cryptobiid Trypanoplasma borreli, gRNA genes appear to reside in large 180-kb noncatenated DNA circles. To shed light on the evolutionary history and function of the minicircle network, we have analyzed minicircle organization in the free-living bodonid Bodo saltans, which is more closely related to trypanosomatids than T. borreli. We identified 1.4-kb circular DNAs as the B. saltans equivalent of minicircles via sequence analysis of 4 complete minicircles, 14 minicircle fragments, and 14 gRNAs. We show that each minicircle harbors two gRNA gene cassettes of opposite polarity residing in variable regions of about 200 nt in otherwise highly conserved molecules. In the conserved region, B. saltans minicircles contain a putative bent helix sequence and a degenerate dodecamer motif (CSB-3). Electron microscopy, sedimentation, and gel electrophoresis analyses showed no evidence for the existence of large minicircle networks in B. saltans, the large majority of the minicircles being present as circular and linear monomers (85-90%) with small amounts of catenated dimers and trimers. Our results provide the first example of a kinetoplastid species with noncatenated, gRNA gene-containing minicircles, which implies that the creation of minicircles and minicircle networks are separate evolutionary events.     

Thiol-based redox metabolism of protozoan parasites

The review considers redox enzymes of Plasmodium spp., Trypanosomatida, Trichomonas, Entamoeba and Giardia, with special emphasis on their potential use as targets for drug development. Thiol-based redox systems play pivotal roles in the success and survival of these parasitic protozoa. The synthesis of cysteine, the key molecule of any thiol metabolism, has been elucidated in trypanosomatids and anaerobes. In trypanosomatids, trypanothione replaces the more common glutathione system. The enzymes of trypanothione synthesis have recently been identified. The role of trypanothione in the detoxification of reactive oxygen species is reflected in the multiplicity of trypanothione-dependent peroxidases. In Plasmodium falciparum, the crystal structures of glutathione reductase and glutamate dehydrogenase are now available; another drug target, thioredoxin reductase, has been demonstrated to be essential for the malarial parasite.     

Mitochondrial and nuclear localization of topoisomerase II in the flagellate Bodo saltans (Kinetoplastida), a species with non-catenated kinetoplast DNA

We have studied topoisomerase II (topo II) in the cells of Bodo saltans, a free-living bodonid (Kinetoplastida). Phylogenetic analysis based on the sequence of the entire topo II gene, which is a single-copy gene, confirmed that B. saltans is a predecessor of parasitic trypanosomatids. Antibodies generated against either an overexpressed unique C-terminal region of topo II or a synthetic oligopeptide derived from the same region did not cross-react with cell lysates of related trypanosomatids, while they recognized a single specific band in the B. saltans lysate. Immunolocalization experiments using both antibodies showed that topo II is evenly dispersed throughout the kinetoplast. This is in striking difference from the localization of topo II in other flagellates, where it occurs in two antipodal centers flanking the kinetoplast disk. Moreover, the same topo II has a distinct localization in multiple loci at the periphery of the nucleus of B. saltans. With a minicircle probe derived from the conserved region we have shown that all relaxed non-catenated minicircles are confined to the globular kinetoplast DNA bundle. Therefore, in the mitochondrion of this primitive eukaryote topo II does not catenate relaxed DNA circles into a network in vivo, while a decatenating activity is present in partially purified cell lysates.     

Cophylogenetic relationships between penguins and their chewing lice

It is generally thought that the evolution of obligate parasites should be linked intimately to the evolution of their hosts and that speciation by the hosts should cause speciation of their parasites. The penguins and their chewing lice present a rare opportunity to examine codivergence between a complete host order and its parasitic lice. We estimated a phylogeny for all 15 species of lice parasitising all 17 species of penguins from the third domain of the mitochondrial 12S ribosomal rRNA gene, a portion of the mitochondrial cytochrome oxidase subunit 1 gene and 55 morphological characters. We found no evidence of extensive cospeciation between penguins and their chewing lice using TreeMap 2.02beta. Despite the paucity of cospeciation, there is support for significant congruence between the louse and penguin phylogenies due to possible failure to speciate events (parasites not speciating in response to their hosts speciating).