De Novo Transcriptome Meta-Assembly of the Mixotrophic Freshwater Microalga Euglena gracilis.

Euglena gracilis is a well-known photosynthetic microeukaryote considered as the product of a secondary endosymbiosis between a green alga and a phagotrophic unicellular belonging to the same eukaryotic phylum as the parasitic trypanosomatids. As its nuclear genome has proven difficult to sequence, reliable transcriptomes are important for functional studies. In this work, we assembled a new consensus transcriptome by combining sequencing reads from five independent studies. Based on a detailed comparison with two previously released transcriptomes, our consensus transcriptome appears to be the most complete so far. Remapping the reads on it allowed us to compare the expression of the transcripts across multiple culture conditions at once and to infer a functionally annotated network of co-expressed genes. Although the emergence of meaningful gene clusters indicates that some biological signal lies in gene expression levels, our analyses confirm that gene regulation in euglenozoans is not primarily controlled at the transcriptional level. Regarding the origin of E. gracilis, we observe a heavily mixed gene ancestry, as previously reported, and rule out sequence contamination as a possible explanation for these observations. Instead, they indicate that this complex alga has evolved through a convoluted process involving much more than two partners.

Orientation of FtsH protease homologs in Trypanosoma brucei inner mitochondrial membrane and its evolutionary implications.

Trypanosoma brucei is an important human pathogen. In this study, we have focused on the characterization of FtsH protease, ATP-dependent membrane-bound mitochondrial enzyme important for regulation of protein abundance. We have determined localization and orientation of all six putative T.brucei FtsH homologs in the inner mitochondrial membrane by in silico analyses, by immunofluorescence, and with protease assay. The evolutionary origin of these homologs has been tested by comparative phylogenetic analysis. Surprisingly, some kinetoplastid FtsH proteins display inverted orientation in the mitochondrial membrane compared to related proteins of other examined eukaryotes. Moreover, our data strongly suggest that during evolution the orientation of FtsH protease in T. brucei varied due to both loss and acquisition of the transmembrane domain.

Shape-based separation of microalga Euglena gracilis using inertial microfluidics.

Euglena gracilis (E. gracilis) has been proposed as one of the most attractive microalgae species for biodiesel and biomass production, which exhibits a number of shapes, such as spherical, spindle-shaped, and elongated. Shape is an important biomarker for E. gracilis, serving as an indicator of biological clock status, photosynthetic and respiratory capacity, cell-cycle phase, and environmental condition. The ability to prepare E. gracilis of uniform shape at high purities has significant implications for various applications in biological research and industrial processes. Here, we adopt a label-free, high-throughput, and continuous technique utilizing inertial microfluidics to separate E. gracilis by a key shape parameter-cell aspect ratio (AR). The microfluidic device consists of a straight rectangular microchannel, a gradually expanding region, and five outlets with fluidic resistors, allowing for inertial focusing and ordering, enhancement of the differences in cell lateral positions, and accurate separation, respectively. By making use of the shape-activated differences in lateral inertial focusing dynamic equilibrium positions, E. gracilis with different ARs ranging from 1 to 7 are directed to different outlets.

Sequence evidence for the presence of two tetrapyrrole pathways in Euglena gracilis.

Genes encoding enzymes of the tetrapyrrole biosynthetic pathway were searched within Euglena gracilis EST databases and 454 genome reads and their 5' end regions were sequenced when not available. Phylogenetic analyses and protein localization predictions support the hypothesis concerning the presence of two separated tetrapyrrole pathways in E. gracilis. One of these pathways resembles the heme synthesis in primarily heterotrophic eukaryotes and was presumably present in the host cell prior to secondary endosymbiosis with a green alga. The second pathway is similar to the plastid-localized tetrapyrrole syntheses in plants and photosynthetic algae. It appears to be localized to the secondary plastid, presumably derived from an algal endosymbiont and probably serves only for the production of plastidial heme and chlorophyll. Thus, E. gracilis represents an evolutionary intermediate in a metabolic transformation of a primary heterotroph to a photoautotroph through secondary endosymbiosis. We propose here that the tetrapyrrole pathway serves as a highly informative marker for the evolution of plastids and plays a crucial role in the loss of plastids.

Euglena gracilis ascorbate peroxidase forms an intramolecular dimeric structure: its unique molecular characterization.

Euglena gracilis lacks a catalase and contains a single APX (ascorbate peroxidase) and enzymes related to the redox cycle of ascorbate in the cytosol. In the present study, a full-length cDNA clone encoding the Euglena APX was isolated and found to contain an open reading frame encoding a protein of 649 amino acids with a calculated molecular mass of 70.5 kDa. Interestingly, the enzyme consisted of two entirely homologous catalytic domains, designated APX-N and APX-C, and an 102 amino acid extension in the N-terminal region, which had a typical class II signal proposed for plastid targeting in Euglena. A computer-assisted analysis indicated a novel protein structure with an intramolecular dimeric structure. The analysis of cell fractionation showed that the APX protein is distributed in the cytosol, but not the plastids, suggesting that Euglena APX becomes mature in the cytosol after processing of the precursor. The kinetics of the recombinant mature FL (full-length)-APX and the APX-N and APX-C domains with ascorbate and H2O2 were almost the same as that of the native enzyme. However, the substrate specificity of the mature FL-APX and the native enzyme was different from that of APX-N and APX-C. The mature FL-APX, but not the truncated forms, could reduce alkyl hydroperoxides, suggesting that the dimeric structure is correlated with substrate recognition. In Euglena cells transfected with double-stranded RNA, the silencing of APX expression resulted in a significant increase in the cellular level of H2O2, indicating the physiological importance of APX to the metabolism of H2O2.

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.

Phylogenetic position and inter-relationships of the osmotrophic euglenids based on SSU rDNA data, with emphasis on the Rhabdomonadales (Euglenozoa).

In order to reconstruct the evolution of euglenid flagellates, euglenozoan SSU rDNA data have been used to investigate phylogenetic relationships with a focus on osmotrophic taxa and especially on the Rhabdomonadales. The dataset consisting of the SSU rDNAs of osmotrophic, phagotrophic and phototrophic taxa was used in parsimony, maximum-likelihood and distance analyses. Five genera make up the Rhabdomonadales, all of them osmotrophic: Gyropaigne, Menoidium, Parmidium, Rhabdomonas and Rhabdospira. According to our analyses they form a strongly supported monophyletic assemblage which is characterized by a low sequence divergence compared to the euglenids in general. Closest relatives are the members of the osmotrophic genus Distigma. All primary osmotrophic species constitute a larger monophyletic group with the phototrophic euglenids and the phagotroph Peranema trichophorum. The combination of three rhabdomonadalian species Rhabdomonas gibba, Rhabdomonas spiralis and Rhabdospira spiralis with nearly identical SSU rDNA sequences is strongly recommended. The phagotroph Petalomonas cantuscygni branches at the bottom of the euglenid subtree with significantly weaker support. The inter-relationship of the three distinct euglenozoan taxa (euglenids, kinetoplastids and diplonemids) could not be convincingly resolved by this study.

Evolutionary appearance of genes encoding proteins associated with box H/ACA snoRNAs: cbf5p in Euglena gracilis, an early diverging eukaryote, and candidate Gar1p and Nop10p homologs in archaebacteria.

A reverse transcription-polymerase chain reaction (RT-PCR) approach was used to clone a cDNA encoding the Euglena gracilis homolog of yeast Cbf5p, a protein component of the box H/ACA class of snoRNPs that mediate pseudouridine formation in eukaryotic rRNA. Cbf5p is a putative pseudouridine synthase, and the Euglena homolog is the first full-length Cbf5p sequence to be reported for an early diverging unicellular eukaryote (protist). Phylogenetic analysis of putative pseudouridine synthase sequences confirms that archaebacterial and eukaryotic (including Euglena) Cbf5p proteins are specifically related and are distinct from the TruB/Pus4p clade that is responsible for formation of pseudouridine at position 55 in eubacterial (TruB) and eukaryotic (Pus4p) tRNAs. Using a bioinformatics approach, we also identified archaebacterial genes encoding candidate homologs of yeast Gar1p and Nop10p, two additional proteins known to be associated with eukaryotic box H/ACA snoRNPs. These observations raise the possibility that pseudouridine formation in archaebacterial rRNA may be dependent on analogs of the eukaryotic box H/ACA snoRNPs, whose evolutionary origin may therefore predate the split between Archaea (archaebacteria) and Eucarya (eukaryotes). Database searches further revealed, in archaebacterial and some eukaryotic genomes, two previously unrecognized groups of genes (here designated 'PsuX' and 'PsuY') distantly related to the Cbf5p/TruB gene family.

The nucleotide sequence of Euglena cytoplasmic phenylalanine transfer RNA. Evidence for possible classifications of Euglena among the animal rather than the plant kingdom.

The nucleotide sequence of cytoplasmic phenylalanine tRNA from Euglena gracilis has been elucidated using procedures described previously for the corresponding chloroplastic tRNA [Cell, 9, 717 (1976)]. The sequence is: pG-C-C-G-A-C-U-U-A-m(2)G-C-U-Cm-A-G-D-D-G-G-G-A-G-A-G-C-m(2)2G-psi-psi-A-G-A-Cm -U-Gm-A-A-Y-A-psi-C-U-A-A-A-G-m(7)G-U-C-*C-C-U-G-G-T-psi-C-G-m(1)A-U-C-C-C-G-G- G-A-G-psi-C-G-G-C-A-C-C-A. Like other tRNA Phes thus far sequenced, this tRNA has a chain length of 76 nucleotides. The sequence of E. gracilis cytoplasmic tRNA Phe is quite different (27 nucleotides out of 76 different) from that of the corresponding chloroplastic tRNA but is surprisingly similar (72 out of 76 nucleotides identical) to that of tRNA Phe from mammalian cytoplasm. This extent of sequence homology even exceeds that found between E. gracilis and wheat germ cytoplasmic tRNA Phe. These findings raise interesting questions on the evolution of tRNAs and the taxonomy of Euglena.

Tryptophan synthetase in Euglena gracilis strain G.

The five enzyme activities in the synthesis of l-tryptophan have been obtained in extracts of Euglena gracilis. One of these, tryptophan synthetase, has been studied in detail. The general catalytic properties of tryptophan synthetase, including the range of reactions catalyzed and its substrate and cofactor affinities, are similar to those reported for other organisms. The Euglena enzyme has two properties never previously observed for tryptophan synthetase. First, the rate of catalysis of the conversion of indole-glycerol phosphate to l-tryptophan remained at its maximal value and was unaffected by the ionic environment up to 0.3 m KCl. In contrast, the conversion of indole to tryptophan showed a sharp maximum at 0.08 m KCl. Second, the enzyme is a component of a complex that includes every enzyme in the pathway committed to tryptophan biosynthesis with the exception of anthranilate synthetase, the regulatory enzyme.