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.     

Glyceraldehyde-3-phosphate dehydrogenase gene diversity in eubacteria and eukaryotes: evidence for intra- and inter-kingdom gene transfer.

Cyanobacteria contain up to three highly divergent glyceraldehyde-3-phosphate dehydrogenase (GAPDH) genes: gap1, gap2, and gap3. Genes gap1 and gap2 are closely related at the sequence level to the nuclear genes encoding cytosolic and chloroplast GAPDH of higher plants and have recently been shown to play distinct key roles in catabolic and anabolic carbon flow, respectively, of the unicellular cyanobacterium Synechocystis sp. PCC6803. In the present study, sequences of 10 GAPDH genes distributed across the cyanobacteria Prochloron didemni, Gloeobacter violaceus PCC7421, and Synechococcus PCC7942 and the alpha-proteobacterium Paracoccus denitrificans and the beta-proteobacterium Ralstonia solanacearum were determined. Prochloron didemni possesses homologs to the gap2 and gap3 genes from Anabaena, Gloeobacter harbors gap1 and gap2 homologs, and Synechococcus possesses gap1, gap2, and gap3. Paracoccus harbors two highly divergent gap genes that are related to gap3, and Ralstonia possesses a homolog of the gap1 gene. Phylogenetic analyses of these sequences in the context of other eubacterial and eukaryotic GAPDH genes reveal that divergence across eubacterial gap1, and gap2, and gap3 genes is greater than that between eubacterial gap1 and eukaroytic glycolytic GapC or between eubacterial gap2 and eukaryotic Calvin cycle GapAB. These data strongly support previous analyses which suggested that eukaryotes acquired their nuclear genes for GapC and GapAB via endosymbiotic gene transfer from the antecedents of mitochondria and chloroplasts, and extend the known range of sequence diversity of the antecedent eubacterial genes. Analyses of available GAPDH sequences from other eubacterial sources indicate that the glycosomal gap gene from trypanosomes (cytosolic in Euglena) and the gap gene from the spirochete Treponema pallidum are each other's closest relatives. This specific relationship can therefore not reflect organismal evolution but must be the result of an interkingdom gene transfer, the direction of which cannot be determined with certainty at present. Contrary to this, the origin of the cytosolic Gap gene from trypanosomes can now be clearly defined as gamma-proteobacterial, since the newly established Ralstonia sequence (beta-proteobacteria) branches basally to the gamma-proteobacterial/trypanosomal assemblage.     

Pan-oceanic distribution of new highly diverse clades of deep-sea diplonemids

Molecular rRNA gene surveys reveal a considerable diversity of microbial eukaryotes in different environments. Even within a single clade, the number of distinct phylotypes retrieved often goes beyond previous expectations. Here, we have used specific 18S rRNA PCR primers to investigate the diversity of diplonemids, a poorly known group of flagellates with only a few described species. We analysed surface and deep-sea plankton samples from different oceanic regions, including the water-column in the Marmara Sea. We retrieved a large diversity of diplonemid phylotypes, most of which formed two novel distinct clades without cultured representatives. Although most marine diplonemid phylotypes appeared to be cosmopolitan, they showed a marked stratified distribution through the water column, being very scarce or absent in surface waters. The small and specific diplonemid diversity found in surface samples and the fact that most sequences of uncultured diplonemids found in other studies came from deep-sea environments suggest that the two major uncultured diplonemid clades group species preferentially inhabit the deep ocean.     

Eukaryotic origin of glyceraldehyde-3-phosphate dehydrogenase genes in Clostridium thermocellum and Clostridium cellulolyticum genomes and putative fates of the exogenous gene in the subsequent genome evolution

Although lateral gene transfer (LGT) events have been frequently documented in the evolution of glyceraldehyde-3-phosphate dehydrogenase (GAPDH), no eukaryote-to-prokaryote transfer has been reported so far. Here we describe the first case of the GAPDH gene transfer from a eukaryote to a subset of Clostridium species (Bacteria, Firmicutes). A series of phylogenetic analyses of GAPDH homologues revealed that Clostridium thermocellum and Clostridium cellulolyticum homologues have the evolutionary affinity to the eukaryotic homologues, rather than to those of bacterial species closely related to the two Clostridium species in the organismal phylogeny. These results suggest that the GAPDH genes in the two Clostridium species are of eukaryotic origin, which is the first reported case of eukaryote-to-bacterium GAPDH gene transfer. Since a previously published 16S ribosomal DNA phylogeny and our GAPDH phylogeny commonly suggest an intimate evolutionary relationship between C. thermocellum and C. cellulolyticum, a common ancestor of the two species likely acquired the eukaryotic GAPDH gene. In the C. cellulolyticum genome, the exogenous GAPDH gene was physically separated from other glycolytic genes, suggesting that this gene organization was likely achieved by a random insertion of the laterally transferred gene. On the other hand, in the C. thermocellum genome, the laterally transferred GAPDH gene clusters with other bacterial glycolytic genes. We discuss possible scenarios for the evolutionarily chimeric glycolytic gene cluster in the C. thermocellum genome.     

Nucleus-encoded,plastid-targeted glyceraldehyde-3-phosphate dehydrogenase (GAPDH) indicates a single origin for chromalveolate plastids

Plastids (the photosynthetic organelles of plants and algae) originated through endosymbiosis between a cyanobacterium and a eukaryote and subsequently spread to other eukaryotes by secondary endosymbioses between two eukaryotes. Mounting evidence favors a single origin for plastids of apicomplexans, cryptophytes, dinoflagellates, haptophytes, and heterokonts (together with their nonphotosynthetic relatives, termed chromalveolates), but so far, no single molecular marker has been described that supports this common origin. One piece of evidence comes from plastid-targeted glyceraldehyde-3-phosphate dehydrogenase (GAPDH), which originated by a gene duplication of the cytosolic form. However, no plastid GAPDH has been characterized from haptophytes, leaving an important piece of the puzzle missing. We have sequenced genes encoding cytosolic, mitochondrion-targeted, and plastid-targeted GAPDH proteins from a number of haptophytes and heterokonts and found haptophyte homologs that branch within a strongly supported clade of chromalveolate plastid-targeted genes, being more closely related to an apicomplexan homolog than was expected. The evolution of plastid-targeted GAPDH supports red algal ancestry of apicomplexan plastids and raises a number of questions about the importance of plastid loss and the possibility of cryptic plastids in nonphotosynthetic lineages such as ciliates.     

Horizontal transfer of a eukaryotic plastid-targeted protein gene to cyanobacteria

Background  

Horizontal or lateral transfer of genetic material between distantly related prokaryotes has been shown to play a major role in the evolution of bacterial and archaeal genomes, but exchange of genes between prokaryotes and eukaryotes is not as well understood. In particular, gene flow from eukaryotes to prokaryotes is rarely documented with strong support, which is unusual since prokaryotic genomes appear to readily accept foreign genes.     

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.     

Horizontal gene transfer in eukaryotic evolution

Horizontal gene transfer (HGT; also known as lateral gene transfer) has had an important role in eukaryotic genome evolution, but its importance is often overshadowed by the greater prevalence and our more advanced understanding of gene transfer in prokaryotes. Recurrent endosymbioses and the generally poor sampling of most nuclear genes from diverse lineages have also complicated the search for transferred genes. Nevertheless, the number of well-supported cases of transfer from both prokaryotes and eukaryotes, many with significant functional implications, is now expanding rapidly. Major recent trends include the important role of HGT in adaptation to certain specialized niches and the highly variable impact of HGT in different lineages.     

Primary Structure and Phylogenetic Relationships of Glyceraldehyde-3-Phosphate Dehydrogenase Genes of Free-Living and Parasitic Diplomonad Flagellates

ABSTRACT. Complete nucleotide sequences have been established for two genes (gap1 and gap2) coding for glyceraldehyde-3-phosphate dehydrogenase (GAPDH, EC 1.2.1.12) homologs in the diplomonad Giardia lamblia. In addition, almost complete sequences of the GAPDH open reading frames were obtained from PCR products for two free-living diplomonad species, Trepomonas agilis and Hexamita inflata, and a parasite of Atlantic salmon, an as yet unnamed species with morphological affinities to Spironucleus. Giardia lamblia gap 1 and the genes from the three other diplomonad species show high similarity to each other and to other glycolytic GAPDH genes. All amino-acyl residues known to be highly conserved in this enzyme are also conserved in these sequences. Giardia lamblia gap2 gene is more divergent and its putative translation reveals the presence of a cysteine and serine-rich insertion resembling a metal binding finger. This motif has not yet been noted in other GAPDH molecules. All sequences contain an S-loop signature with characteristics close to those of eukaryotes. In phylogenetic reconstructions based on the derived amino acid sequences with neighborjoining, parsimony and maximum-likelihood methods the four typical GAPDH sequences of diplomonads cluster into a single clade. Within this clade, G. lamblia gap1 shares a common ancestor with the rest of the genes. The latter are more closely related to each other, indicating an early separation of the lineage leading to the genus Giardia from the lineage encompassing the morphologically less differentiated genera, Trepomonas, Hexamita and that of the unnamed species. This result is discordant with the orthogonal evolution of diplomonads suggested on the basis of comparative morphology. In neighbor-joining reconstructions G. lamblia gap2 occupies a variable position, due to its great divergence. In parsimony and maximum likelihood analysis however, it shares a most recent common ancestor with the typical G. lamblia gap1 gene, suggesting that it diverged after the separation of the Giardia lineage. The position of the diplomonad clade in broader phylogenetic reconstructions is firmly within the typical cytosolic glycolytic representatives of GAPDH of eukaryotes.     

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.     

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

Lateral transfer of eukaryotic ribosomal RNA genes: an emerging concern for molecular ecology of microbial eukaryotes

Ribosomal RNA (rRNA) genes are widely utilized in depicting organismal diversity and distribution in a wide range of environments. Although a few cases of lateral transfer of rRNA genes between closely related prokaryotes have been reported, it remains to be reported from eukaryotes. Here, we report the first case of lateral transfer of eukaryotic rRNA genes. Two distinct sequences of the 18S rRNA gene were detected from a clonal culture of the stramenopile, Ciliophrys infusionum. One was clearly derived from Ciliophrys, but the other gene originated from a perkinsid alveolate. Genome-walking analyses revealed that this alveolate-type rRNA gene is immediately adjacent to two protein-coding genes (ubc12 and usp39), and the origin of both genes was shown to be a stramenopile (that is, Ciliophrys) in our phylogenetic analyses. These findings indicate that the alveolate-type rRNA gene is encoded on the Ciliophrys genome and that eukaryotic rRNA genes can be transferred laterally.     

GAPDH gene diversity in spirochetes: a paradigm for genetic promiscuity.

In this study we have determined gap sequences from nine different spirochetes. Phylogenetic analyses of these sequences in the context of all other available eubacterial and a selection of eukaryotic Gap sequences demonstrated that the eubacterial glyceraldehyde-3-phosphate dehydrogenase (GAPDH) gene diversity encompasses at least five highly distinct gene families. Within these gene families, spirochetes show an extreme degree of sequence divergence that is probably the result of several lateral gene transfer events between spirochetes and other eubacterial phyla, and early gene duplications in the eubacterial ancestor. A Gap1 sequence from the syphilis spirochete Treponema pallidum has recently been shown to be closely related to GapC sequences from Euglenozoa. Here we demonstrate that several other spirochetal species are part of this cluster, supporting the conclusion that an interkingdom gene transfer from spirochetes to Euglenozoa must have occurred. Furthermore, we provide evidence that the GAPDH genes present in the protists Parabasalia may also be of spirochetal descent.     

Lateral gene transfer and the complex distribution of insertions in eukaryotic enolase

Insertions and deletions in protein-coding genes are relatively rare events compared with sequence substitutions because they are more likely to alter the tertiary structure of the protein. For this reason, insertions and deletions which are clearly homologous are considered to be stable characteristics of the proteins where they are found, and their presence and absence has been used extensively to infer large-scale evolutionary relationships and events. Recently, however, it has been shown that the pattern of highly conserved, clearly homologous insertions at positions with no other detectable homoplasy can be incongruent with the phylogeny of the genes or organisms in which they are found. One case where this has been reported is in the enolase genes of apicomplexan parasites and ciliates, which share homologous insertions in a highly conserved region of the gene with the apparently distantly related enolases of plants. Here we explore the distribution of this character in enolase genes from the third major alveolate group, the dinoflagellates, as well as two groups considered to be closely related to alveolates, haptophytes and heterokonts. With these data, all major groups of the chromalveolates are represented, and the distribution of these insertions is shown to be far more complicated than previously believed. The incongruence between this pattern, the known evolutionary relationships between the organisms, and enolase phylogeny itself cannot be explained by any single event or type of event. Instead, the distribution of enolase insertions is more likely the product of several forces that may have included lateral gene transfer, paralogy, and/or recombination. Of these, lateral gene transfer is the easiest to detect and some well-supported cases of eukaryote-to-eukaryote lateral transfer are evident from the phylogeny.     

64 Nucleus-encoded, plastid-targeted glyceraldehyde-3-phosphate dehydrogenase (GAPDH) indicates a single origin for chromalveolate plastids

Plastids (the photosynthetic organelles of plants and algae) ultimately originated through an endosymbiosis between a cyanobacterium and a eukaryote. Subsequently, plastids spread to other eukaryotes by secondary endosymbioses that took place between a eukaryotic alga and a second eukaryote. Recently, evidence has mounted in favour of a single origin for plastids of apicomplexans, cryptophytes, dinoflagellates, haptophytes, and heterokonts (together with their non-photosynthetic relatives, collectively termed chromalveolates). As of yet, however, no single molecular marker has been described which supports a common origin for all of these plastids. One piece of the evidence for a single origin of chromalveolate plastids came from plastid-targeted glyceraldehyde-3-phosphate dehydrogenase (GAPDH), which originated by a gene duplication of the cytosolic form. However, no plastid GAPDH has been characterized from haptophytes, leaving an important piece of the puzzle missing. We have sequenced genes encoding cytosolic, mitochondrial-targeted, and plastid-targeted GAPDH proteins from a number of haptophytes and heterokonts, and found the haptophyte homologues to branch within the strongly supported clade of chromalveolate plastid-targeted GAPDH genes. Interestingly, plastid-targeted GAPDH genes from the haptophytes were more closely related to apicomplexan genes than was expected. Overall, the evolution of plastid-targeted GAPDH reinforces other data for a red algal ancestry of apicomplexan plastids, and raises a number of questions about the importance of plastid loss and the possibility of cryptic plastids in non-photosynthetic lineages such as ciliates.     

The evolutionary origin of red algae as deduced from the nuclear genes encoding cytosolic and chloroplast glyceraldehyde-3-phosphate dehydrogenases from Chondrus crispus

Algae are a heterogeneous group of photosynthetic eukaryotes traditionally separated into three major subdivisions: rhodophytes, chlorophytes, and chromophytes. The evolutionary origin of rhodophytes or red algae and their links to other photosynthetic and nonphotosynthetic eukaryotes have been a matter of much controversy and speculation. Here we present the first cDNAs of nuclear protein genes from red algae: Those encoding cytosolic and chloroplast glyceraldehyde-3-phosphate dehydrogenases (GAPDH) from Chondrus crispus. A phylogenetic analysis including GAPDH gene sequences from a number of eukaryotic taxa, cyanobacteria, and purple bacteria suggests that chloroplasts and rhodoplasts together form a monophyletic group of cyanobacterial descent and that rhodophytes separated from chlorophytes at about the same time as animals and fungi. The composite GAPDH tree further demonstrates that chloroplast and cytosolic GAPDH genes are closely related to their homologs in cyanobacteria and purple bacteria, respectively, the presumptive ancestors of chloroplasts and mitochondria, thereby firmly establishing the endosymbiotic origin of these nuclear genes and their fixation in eukaryotic cells before the rhodophyte/chlorophyte separation. The present data are in conflict with phylogenetic inferences based on plastid-encoded rbcL sequences supporting a polyphyletic origin of rhodoplasts and chloroplasts. Comparison of rbcL to GAPDH phylogenies suggests that rbcL trees may be misleading because they are composed of branches representing ancient duplicated (paralogous) genes. Correspondence to: R. Cerff