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.     

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.     

Unusual mitochondrial genome structures throughout the Euglenozoa

Mitochondrial DNA of Kinetoplastea is composed of different chromosomes, the maxicircle (bearing 'regular' genes) and numerous minicircles (specifying guide RNAs involved in RNA editing). In trypanosomes [Kinetoplastea], DNA circles are compacted into a single dense body, the kinetoplast. This report addresses the question whether multi-chromosome mitochondrial genomes and compacted chromosome organization are restricted to Kinetoplastea or rather occur throughout Euglenozoa, i.e., Kinetoplastea, Euglenida and Diplonemea. To this end, we investigated the diplonemid Rhynchopus euleeides and the euglenids Petalomonas cantuscygni, Peranema trichophorum and Entosiphon sulcatum, using light and electron microscopy and molecular techniques. Our findings together with previously published data show that multi-chromosome mitochondrial genomes prevail across Euglenozoa, while kinetoplast-like mtDNA packaging is confined to trypanosomes.     

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.     

Rapid evolution of the compact and unusual mitochondrial genome in the ctenophore, Pleurobrachia bachei

Ctenophores are one of the most basally branching lineages of metazoans with the largest mitochondrial organelles in the animal kingdom. We sequenced the mitochondrial (mtDNA) genome from the Pacific cidipid ctenophore, Pleurobrachia bachei. The circular mitochondrial genome is 11,016 nts, with only 12 genes, and one of the smallest metazoan mtDNA genomes recorded. The protein coding genes are intronless cox1-3, cob, nad1, 3, 4, 4L and 5. The nad2 and 6 genes are represented as short fragments whereas the atp6 gene was found in the nuclear genome. Only the large ribosomal RNA subunit and two tRNAs were present with possibly the small subunit unidentifiable due to extensive fragmentation. The observed unique features of this mitochondrial genome suggest that nuclear and mitochondrial genomes have evolved at very different rates. This reduced mtDNA genome sharply contrasts with the very large sizes of mtDNA found in other basal metazoans including Porifera (sponges), and Placozoa (Trichoplax).     

Order within a mosaic distribution of mitochondrial c-type cytochrome biogenesis systems?

Mitochondrial cytochromes c and c(1) are present in all eukaryotes that use oxygen as the terminal electron acceptor in the respiratory chain. Maturation of c-type cytochromes requires covalent attachment of the heme cofactor to the protein, and there are at least five distinct biogenesis systems that catalyze this post-translational modification in different organisms and organelles. In this study, we use biochemical data, comparative genomic and structural bioinformatics investigations to provide a holistic view of mitochondrial c-type cytochrome biogenesis and its evolution. There are three pathways for mitochondrial c-type cytochrome maturation, only one of which is present in prokaryotes. We analyze the evolutionary distribution of these biogenesis systems, which include the Ccm system (System I) and the enzyme heme lyase (System III). We conclude that heme lyase evolved once and, in many lineages, replaced the multicomponent Ccm system (present in the proto-mitochondrial endosymbiont), probably as a consequence of lateral gene transfer. We find no evidence of a System III precursor in prokaryotes, and argue that System III is incompatible with multi-heme cytochromes common to bacteria, but absent from eukaryotes. The evolution of the eukaryotic-specific protein heme lyase is strikingly unusual, given that this protein provides a function (thioether bond formation) that is also ubiquitous in prokaryotes. The absence of any known c-type cytochrome biogenesis system from the sequenced genomes of various trypanosome species indicates the presence of a third distinct mitochondrial pathway. Interestingly, this system attaches heme to mitochondrial cytochromes c that contain only one cysteine residue, rather than the usual two, within the heme-binding motif. The isolation of single-cysteine-containing mitochondrial cytochromes c from free-living kinetoplastids, Euglena and the marine flagellate Diplonema papillatum suggests that this unique form of heme attachment is restricted to, but conserved throughout, the protist phylum Euglenozoa.     

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.     

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.     

The complete mitochondrial DNA sequence of Mesostigma viride identifies this green alga as the earliest green plant divergence and predicts a highly compact mitochondrial genome in the ancestor of all green plants.

To gain insights into the nature of the mitochondrial genome in the common ancestor of all green plants, we have completely sequenced the mitochondrial DNA (mtDNA) of Mesostigma viride. This green alga belongs to a morphologically heterogeneous class (Prasinophyceae) that includes descendants of the earliest diverging green plants. Recent phylogenetic analyses of ribosomal RNAs (rRNAs) and concatenated proteins encoded by the chloroplast genome identified Mesostigma as a basal branch relative to the Streptophyta and the Chlorophyta, the two phyla that were previously thought to contain all extant green plants. The circular mitochondrial genome of Mesostigma resembles the mtDNAs of green algae occupying a basal position within the Chlorophyta in displaying a small size (42,424 bp) and a high gene density (86.6% coding sequences). It contains 65 genes that are conserved in other mtDNAs. Although none of these genes represents a novel coding sequence among green plant mtDNAs, four of them (rps1, sdh3, sdh4, and trnL[caa]) have not been reported previously in chlorophyte mtDNAs, and two others (rpl14 and trnI[gau]) have not been identified in the streptophyte mtDNAs examined so far (land-plant mtDNAs). Phylogenetic analyses of 19 concatenated mtDNA-encoded proteins favor the hypothesis that Mesostigma represents the earliest branch of green plant evolution. Four group I introns (two in rnl and two in cox1) and three group II introns (two in nad3 and one in cox2), two of which are trans-spliced at the RNA level, reside in Mesostigma mtDNA. The insertion sites of the three group II introns are unique to this mtDNA, suggesting that trans-splicing arose independently in the Mesostigma lineage and in the Streptophyta. The few structural features that can be regarded as ancestral in Mesostigma mtDNA predict that the common ancestor of all green plants had a compact mtDNA containing a minimum of 75 genes and perhaps two group I introns. Considering that the mitochondrial genome is much larger in size in land plants than in Mesostigma, we infer that mtDNA size began to increase dramatically in the Streptophyta either during the evolution of charophyte green algae or during the transition from charophytes to land plants.     

The mitochondrial genome of the Basidiomycete fungus Pleurotus ostreatus (oyster mushroom)

In this study, the full mitochondrial genome of a basidiomycete fungus, Pleurotus ostreatus, was sequenced and analyzed. It is a circular DNA molecule of 73 242 bp and contains 44 known genes encoding 18 proteins and 26 RNA genes. The protein-coding genes include 14 common mitochondrial genes, one ribosomal small subunit protein 3 gene, one RNA polymerase gene and two DNA polymerase genes. In addition, one RNA and one DNA polymerase genes were identified in a mitochondrial plasmid. These two genes show relatively low similarities to their homologs in the mitochondrial genome but they are nearly identical to the known mitochondrial plasmid genes from another Pleurotus ostreatus strain. This suggests that the plasmid may mediate the horizontal gene transfer of the DNA and RNA polymerase genes into mitochondrial genome, and such a transfer may be an ancient event. Phylogenetic analysis based on the cox1 ORFs verified the traditional classification of Pleurotus ostreatus among fungi. However, the discordances were observed in the phylogenetic trees based on the six cox1 intronic ORFs of Pleurotus ostreatus and their homologs in other species, suggesting that these intronic ORFs are foreign DNA sequences obtained through HGT. In summary, this analysis provides valuable information towards the understanding of the evolution of fungal mtDNA.     

Complete mitochondrial genome sequence of D. littoralis (Diptera: Drosophilidae). Comparative analysis of mitochondrial genomes in the Drosophila virilis group

An analysis of the complete sequence of the mitochondrial genome (mt-genome) of D. littoralis is presented. Its basic characteristics, such as size, nucleotide composition, gene order, and the degree of purifying selection pressure on different parts of the genome are given. The details of the structure of proteinencoding genes and tRNA genes are discussed. The structure of nonencoding regions (control region and intergenic spacers) of mtDNA of the virilis group is given. Fragments of the mt-genes atp6 and cox3 were found in the nuclear genome of D. virilis. The evolutionary history of the mitochondrial and nuclear sequences of these genes indicates that the process of formation of mt-pseudogenes is currently taking place and is associated with the activity of retrotransposons.     

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

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.     

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.     

The complete mitochondrial genome of an ectoparasitic monopisthocotylean fluke Benedenia hoshinai (Monogenea: Platyhelminthes)

An exponential growth of mitochondrial genome information has brought significant progress in understanding the organismal phylogeny and mitochondrial genome evolution for many metazoans including platyhelminth groups. In this study, we determined the complete mitochondrial genome sequence for Benedenia hoshinai, an ectoparasitic monogenean species, and compared it with its congener Benedenia seriolae. The complete mitochondrial genome is 13,554 bp in length and contains 12 protein-coding genes (lacking the atp8 gene), 2 rRNA genes, and 22 tRNA genes, all encoded in the same direction as found in all other platyhelminth species sequenced to date. The gene arrangement of B. hoshinai mtDNA is almost identical to B. seriolae, differing only by the translocation of trnT between cox1 and rrnL. It is unclear whether the shared position of trnT between B. hoshinai and Gyrodactylus represents evidence for their phylogenetic affinity; testing this hypothesis requires further mitogenomic evidence.     

The complete mitochondrial genome of the verongid sponge Aplysina cauliformis: implications for DNA barcoding in demosponges

DNA “barcoding,” the determination of taxon-specific genetic variation typically within a fragment of the mitochondrial cytochrome oxidase 1 (cox1) gene, has emerged as a useful complement to morphological studies, and is routinely used by expert taxonomists to identify cryptic species and by non-experts to better identify samples collected during field surveys. The rate of molecular evolution in the mitochondrial genomes (mtDNA) of nonbilaterian animals (sponges, cnidarians, and placozoans) is much slower than in bilaterian animals for which DNA barcoding strategies were developed. If sequence divergence among nonbilaterian mtDNA and specifically cox1 is too slow to generate diagnostic variation, alternative genes for DNA barcoding and species-level phylogenies should be considered. Previous study across the Aplysinidae (Demospongiae, Verongida) family of sponges demonstrated no nucleotide substitutions in the traditional cox1 barcoding fragment among the Caribbean species of Aplysina. As the mitochondrial genome of Aplysina fulva has previously been sequenced, we are now able to make the first comparisons between complete mtDNA of congeneric demosponges to assess whether potentially informative variation exists in genes other than cox1. In this article, we present the complete mitochondrial genome of Aplysina cauliformis, a circular molecule 19620 bp in size. The mitochondrial genome of A. cauliformis is the same length as is A. fulva and shows six confirmed nucleotide differences and an additional 11 potential SNPs. Of the six confirmed SNPs, NADH dehydrogenase subunit 5 (nad5) and nad2 each contain two, and in nad2 both yield amino acid substitutions, suggesting balancing selection may act on this gene. Thus, while the low nucleotide diversity in Caribbean aplysinid cox1 extends to the entire mitochondrial genome, some genes do display variation. If these represent interspecific differences, then they may be useful alternative markers for studies in recently diverged sponge clades.