De novo assembly and comparative transcriptome analysis of Euglena gracilis in response to anaerobic conditions

Advances in euglenoid genomics: unravelling the fascinating biology of a complex clade

Euglenids have long been studied due to their unique physiology and versatile metabolism, providing underpinnings for much of our understanding of photosynthesis and biochemistry, and a growing opportunity in biotechnology. Until recently there has been a lack of genetic studies due to their large and complex genomes, but recently new technologies have begun to unveil their genetic capabilities. Whilst much research has focused on the model organism Euglena gracilis, other members of the euglenids have now started to receive due attention. Currently only poor nuclear genome assemblies of E. gracilis and Rhabdomonas costata are available, but there are many more plastid genome sequences and an increasing number of transcriptomes. As more assemblies become available, there are great opportunities to understand the fundamental biology of these organisms and to exploit them for biotechnology.

Genome editing-based mutagenesis stably modifies composition of wax esters synthesized by Euglena gracilis under anaerobic conditions

Microalgal oil production represents a promising renewable biofuel source. Metabolic engineering can enhance its utility, transforming it into an improved biofuel and expanding its applications as a feedstock for commodity chemicals, thereby increasing their value in biorefineries. This study focused on anaerobic wax ester production by the microalga Euglena gracilis, aiming to develop stable mutant strains with altered wax ester profiles through genome editing. Two enzymes in the fatty acid beta-oxidation pathway involved in wax ester production were targeted-3-ketoacyl-CoA thiolase and acyl-CoA dehydrogenase-using clustered regularly interspaced short palindromic repeats/Cas9. The results revealed one genetic mutation that lengthened and three that shortened the distribution of wax ester compositions compared to the wild-type (WT). The triple-KO mutant, combining mutations that shorten wax ester chains, produced wax esters with acyl chains two carbons shorter than WT. This study established a methodology to stably modify wax ester composition in E. gracilis.

Biochemical and Functional Profiling of Thioredoxin-Dependent Cytosolic GPX-like Proteins in Euglena gracilis

Unlike plants and animals, the phytoflagellate Euglena gracilis lacks catalase and contains a non-selenocysteine glutathione peroxidase-like protein (EgGPXL), two peroxiredoxins (EgPrx1 and EgPrx4), and one ascorbate peroxidase in the cytosol to maintain reactive oxygen species (ROS) homeostasis. In the present study, the full-length cDNA of three cytosolic EgGPXLs was obtained and further characterized biochemically and functionally. These EgGPXLs used thioredoxin instead of glutathione as an electron donor to reduce the levels of H2O2 and t-BOOH. The specific peroxidase activities of these enzymes for H2O2 and t-BOOH were 1.3 to 4.9 and 0.79 to 3.5 µmol/min/mg protein, respectively. Cytosolic EgGPXLs and EgPrx1/EgPrx4 were silenced simultaneously to investigate the synergistic effects of these genes on the physiological function of E. gracilis. The suppression of cytosolic EgGPXL genes was unable to induce any critical phenomena in Euglena under normal (100 μmol photons m-2 s-1) and high-light conditions (350 μmol photons m-2 s-1) at both autotrophic and heterotrophic states. Unexpectedly, the suppression of EgGPXL genes was able to rescue the EgPrx1/EgPrx4-silenced cell line from a critical situation. This study explored the potential resilience of Euglena to ROS, even with restriction of the cytosolic antioxidant system, indicating the involvement of some compensatory mechanisms.

Ras superfamily GTPases and signal transduction in Euglena gracilis

Biological complexity is challenging to define, but can be considered through one or more features, including overall genome size, number of genes, morphological features, multicellularity, number of life cycle stages and the ability to adapt to different environments. Euglena gracilis meets several of these criteria, with a large genome of ∼38,000 protein coding genes and a considerable ability to survive under many different conditions, some of which can be described as challenging or harsh. Potential molecular exemplars of complexity tying these aspects together are signalling pathways, including GTPases, kinases and ubiquitylation, which increase the functionality of the gene-encoded proteome manyfold. Each of these examples can modulate both protein activity and gene expression. To address the connection between genome size and complexity I have undertaken a brief, and somewhat qualitative, survey of the small ras-like GTPase superfamily of E. gracilis. Unexpectedly, apart from Rab-GTPases which control intracellular transport and organelle identify, the size of the GTPase cohort is modest, and, for example, has not scaled with gene number when compared to the close relatives, trypanosomatids. I suggest that understanding the functions of this protein family will be vital to uncovering the complexity of E. gracilis biology.

Transcriptome analysis of the bloom-forming dinoflagellate Prorocentrum donghaiense exposed to Ginkgo biloba leaf extract, with an emphasis on photosynthesis

Ginkgo biloba leaf extract (GBE) can effectively treat bloom-forming freshwater algae. However, there is limited information about the underlying suppression mechanism of the marine bloom-forming Prorocentrum donghaiense-the most dominant algal bloom species in the East China Sea. We investigated the effect of GBE on P. donghaiense in terms of its response to photosynthesis at the molecular/omic level. In total, 93,743 unigenes were annotated using six functional databases. Furthermore, 67,203 differentially expressed genes (DEGs) were identified in algae treated with 1.8 g∙L-1 GBE. Among these DEGs, we identified the genes involved in photosynthesis. PsbA, PsbB and PsbD in photosystem II, PsaA in photosystem I, and PetB and PetD in the cytochrome b6/f complex were downregulated. Other related genes, such as PsaC, PsaE, and PsaF in photosystem I; PetA in the cytochrome b6/f complex; and atpA, atpD, atpH, atpG, and atpE in the F-type H+-ATPase were upregulated. These results suggest that the structure and activity of the complexes were destroyed by GBE, thereby inhibiting the electron flow between the primary and secondary quinone electron acceptors, primary quinone electron acceptor, and oxygen-evolving complex in the PSII complex, and interrupting the electron flow between PSII and PSI, ultimately leading to a decline in algal cell photosynthesis. These findings provide a basis for understanding the molecular mechanisms underlying P. donghaiense exposure to GBE and a theoretical basis for the prevention and control of harmful algal blooms.

Understanding wax ester synthesis in Euglena gracilis: Insights into mitochondrial anaerobic respiration

Euglena gracilis, photosynthetic protist, has a unique ability to generate wax esters in the absence of oxygen, employing a distinctive fatty acid synthesis mechanism. Through comprehensive inhibitor assays and gene-silencing techniques, our research clearly emphasized the indispensable role of the mitochondrial anaerobic respiratory chain in this biosynthesis. We identified acyl-CoA dehydrogenase, electron transfer flavoprotein (ETF), and rhodoquinone (RQ) as central molecular components in the pathway. These findings strongly indicated a potential reversal of beta-oxidation occurring within mitochondria for fatty acid production in anaerobic conditions. Furthermore, our analysis revealed the pivotal function of nicotinamide nucleotide transhydrogenase (NNT) in efficiently managing the NADPH/NAD+ conversion essential for sustaining anaerobic metabolism. This review outlines our key findings and provides a comprehensive understanding of the molecular mechanisms that enable E. gracilis to produce wax ester anaerobically.

Transcriptomic and genomic identification of spliceosomal genes from Euglena gracilis

Diverse splicing types in nuclear and chloroplast genes of protist Euglena gracilis have been recognized for decades. However, the splicing machinery responsible for processing nuclear precursor messenger RNA introns, including trans-splicing of the 5' terminal outron and spliced leader (SL) RNA, remains elusive. Here, we identify 166 spliceosomal protein genes and two snRNA genes from E. gracilis by performing bioinformatics analysis from a combination of next-generation and full-length transcriptomic RNA sequencing (RNAseq) data as well as draft genomic data. With the spliceosomal proteins we identified in hand, the insensitivity of E. gracilis to some splicing modulators is revealed at the sequence level. The prevalence of SL RNA-mediated trans-splicing is estimated to be more than 70% from our full-length RNAseq data. Finally, the splicing proteomes between E. gracilis and its three evolutionary cousins within the same Excavata group are compared. In conclusion, our study characterizes the spliceosomal components in E. gracilis and provides the molecular basis for further exploration of underlying splicing mechanisms.

Biofuel production from Euglena: Current status and techno-economic perspectives

Sustainable aviation fuels (SAFs) can contribute reduce greenhouse gas emissions compared to conventional fuel. With the increasing SAFs demand, various generations of resources have been shifted from the 1st generation (oil crops), the 2nd generation (agricultural waste), to the 3rd generation (microalgae). Microalgae are the most suitable feedstock for jet biofuel production than other resources because of their productivity and capability to capture carbon dioxide. However, microalgae-based biofuel has a limitation of high freezing point. Recently, a jet biofuel derived from Euglena wax ester has been paying attention due to its low freezing point. Challenges still remain to enhance production yields in both upstream and downstream processes. Studies on downstream processes as well as techno-economic analysis on biofuel production using Euglena are highly limited to date. Economic aspects for the biofuel production will be ensured via valorization of industrial byproducts such as food wastes.

Metabolic network reconstruction of Euglena gracilis: Current state, challenges, and applications

A metabolic model, representing all biochemical reactions in a cell, is a prerequisite for several approaches in systems biology used to explore the metabolic phenotype of an organism. Despite the use of Euglena in diverse industrial applications and as a biological model, there is limited understanding of its metabolic network capacity. The unavailability of the completed genome data and the highly complex evolution of Euglena are significant obstacles to the reconstruction and analysis of its genome-scale metabolic model. In this mini-review, we discuss the current state and challenges of metabolic network reconstruction in Euglena gracilis. We have collated and present the available relevant data for the metabolic network reconstruction of E. gracilis, which could be used to improve the quality of the metabolic model of E. gracilis. Furthermore, we deliver the potential applications of the model in metabolic engineering. Altogether, it is supposed that this mini-review would facilitate the investigation of metabolic networks in Euglena and further lay out a direction for model-assisted metabolic engineering.

Functional and structural characterization of an endo-β-1,3-glucanase from Euglena gracilis

Endo-β-1,3-glucanases from several organisms have attracted much attention in recent years because of their capability for in vitro degrading β-1,3-glucan as a critical step for both biofuels production and short-chain oligosaccharides synthesis. In this study, we biochemically characterized a putative endo-β-1,3-glucanase (EgrGH64) belonging to the family GH64 from the single-cell protist Euglena gracilis. The gene coding for the enzyme was heterologously expressed in a prokaryotic expression system supplemented with 3% (v/v) ethanol to optimize the recombinant protein right folding. Thus, the produced enzyme was highly purified by immobilized-metal affinity and gel filtration chromatography. The enzymatic study demonstrated that EgrGH64 could hydrolyze laminarin (K_M 23.5 mg ml^-1,k_cat 1.20 s^-1) and also, but with less enzymatic efficiency, paramylon (K_M 20.2 mg ml^-1,k_cat 0.23 ml mg^-1 s^-1). The major product of the hydrolysis of both substrates was laminaripentaose. The enzyme could also use ramified β-glucan from the baker's yeast cell wall as a substrate (K_M 2.10 mg ml^-1, k_cat 0.88 ml mg^-1 s^-1). This latter result, combined with interfacial kinetic analysis evidenced a protein's greater efficiency for the yeast polysaccharide, and a higher number of hydrolysis sites in the β-1,3/β-1,6-glucan. Concurrently, the enzyme efficiently inhibited the fungal growth when used at 1.0 mg/mL (15.4 μM). This study contributes to assigning a correct function and determining the enzymatic specificity of EgrGH64, which emerges as a relevant biotechnological tool for processing β-glucans.

Fatty Acid Accumulations and Transcriptome Analyses Under Different Treatments in a Model Microalga Euglena gracilis

With the continuous growth of the world’s population and the increasing development of industrialization, the demand for energy by human beings has been expanding, resulting in an increasingly severe energy crisis. Microalgae are considered the most potential alternatives to traditional fossil fuels due to their many advantages, like fast growth rate, strong carbon sequestration capacity, and low growth environment requirements. Euglena can use carbon sources such as glucose, ethanol, and others for heterotrophic growth. Moreover, Euglena is highly adaptable to the environment and has a high tolerance to various environmental stresses, such as salinity, heavy metals, antibiotics, etc. Different treatments of Euglena cells could affect their growth and the accumulation of bioactive substances, especially fatty acids. To expand the industrial application of Euglena as a potential biodiesel candidate, we determine the physiological responses of Euglena against environmental stresses (antibiotics, heavy metals, salinity) or carbon resources (glucose and ethanol), and evaluate the potential for higher quality and yield of fatty acid with a high growth rate. Adding glucose into the culture media increases cell biomass and fatty acid production with high-quality biodiesel characters. The transcriptome analysis helped explore the possible regulation and biosynthesis of fatty acids under different treatments and exploited in the improvement of biodiesel production. This study provides insights for further improvement and various culture treatments for Euglena-based biodiesel and jet fuels.

Suppression of the Lycopene Cyclase Gene Causes Downregulation of Ascorbate Peroxidase Activity and Decreased Glutathione Pool Size, Leading to H2O2 Accumulation in Euglena gracilis

Carotenoids are photosynthetic pigments and hydrophobic antioxidants that are necessary for the survival of photosynthetic organisms, including the microalga Euglena gracilis. In the present study, we identified an uncharacterized gene encoding the E. gracilis β-carotene synthetic enzyme lycopene cyclase (EgLCY) and discovered a relationship between EgLCY-mediated carotenoid synthesis and the reactive oxygen species (ROS) scavenging system ascorbate-glutathione cycle. The EgLCY cDNA sequence was obtained via homology searching E. gracilis transcriptome data. An enzyme assay using Escherichia coli demonstrated that EgLCY converts lycopene to β-carotene. E. gracilis treated with EgLCY double-stranded RNA (dsRNA) produced colorless cells with hypertrophic appearance, inhibited growth, and marked decrease in carotenoid and chlorophyll content, suggesting that EgLCY is essential for the synthesis of β-carotene and downstream carotenoids, which are abundant and physiologically functional. In EgLCY dsRNA-treated cells, the ascorbate-glutathione cycle, composed of ascorbate peroxidase (APX), dehydroascorbate reductase (DHAR), monodehydroascorbate reductase (MDAR), and glutathione reductase (GR), was unusually modulated; APX and GR activities significantly decreased, whereas DHAR and MDAR activities increased. Ascorbate content was significantly increased and glutathione content significantly decreased in EgLCY dsRNA-treated cells and was correlated with their recycling enzyme activities. Fluorescent imaging demonstrated that EgLCY dsRNA-treated cells accumulated higher levels of H2O2 compared to wild-type cells. Taken together, this study revealed that EgLCY-mediated synthesis of β-carotene and downstream carotenoid species upregulates APX activity and increases glutathione pool size for H2O2 scavenging. Our study suggests a possible relationship between carotenoid synthesis and the ascorbate-glutathione cycle for ROS scavenging in E. gracilis.

NADPH-to-NADH conversion by mitochondrial transhydrogenase is indispensable for sustaining anaerobic metabolism in Euglena gracilis

Euglena gracilis produces ATP in the anaerobic mitochondria with concomitant wax ester formation, and NADH is essential for ATP formation and fatty acid synthesis in the mitochondria. This study demonstrated that mitochondrial cofactor conversion by nicotinamide nucleotide transhydrogenase (NNT), converting NADPH/NAD⁺ to NADP⁺ /NADH, is indispensable for sustaining anaerobic metabolism. Silencing of NNT genes significantly decreased wax ester production and cellular viability during anaerobiosis but had no such marked effects under aerobic conditions. An analogous phenotype was observed in the silencing of the gene encoding a mitochondrial NADP⁺ -dependent malic enzyme. These results suggest that the reducing equivalents produced in glycolysis are shuttled to the mitochondria as malate, where cytosolic NAD⁺ regeneration is coupled with mitochondrial NADPH generation.

Heterotrophic euglenid Rhabdomonas costata resembles its phototrophic relatives in many aspects of molecular and cell biology

Euglenids represent a group of protists with diverse modes of feeding. To date, only a partial genomic sequence of Euglena gracilis and transcriptomes of several phototrophic and secondarily osmotrophic species are available, while primarily heterotrophic euglenids are seriously undersampled. In this work, we begin to fill this gap by presenting genomic and transcriptomic drafts of a primary osmotroph, Rhabdomonas costata. The current genomic assembly length of 100 Mbp is 14× smaller than that of E. gracilis. Despite being too fragmented for comprehensive gene prediction it provided fragments of the mitochondrial genome and comparison of the transcriptomic and genomic data revealed features of its introns, including several candidates for nonconventional types. A set of 39,456 putative R. costata proteins was predicted from the transcriptome. Annotation of the mitochondrial core metabolism provides the first data on the facultatively anaerobic mitochondrion of R. costata, which in most respects resembles the mitochondrion of E. gracilis with a certain level of streamlining. R. costata can synthetise thiamine by enzymes of heterogenous provenances and haem by a mitochondrial-cytoplasmic C4 pathway with enzymes orthologous to those found in E. gracilis. The low percentage of green algae-affiliated genes supports the ancestrally osmotrophic status of this species.

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.

Elucidating carbohydrate metabolism in Euglena gracilis: Reverse genetics-based evaluation of genes coding for enzymes linked to paramylon accumulation

Euglena gracilis is a eukaryotic single-celled and photosynthetic organism grouped under the kingdom Protista. This phytoflagellate can accumulate the carbon photoassimilate as a linear β-1,3-glucan chain called paramylon. This storage polysaccharide can undergo degradation to provide glucose units to obtain ATP and reducing power both in aerobic and anaerobic growth conditions. Our group has recently characterized an essential enzyme for accumulating the polysaccharide, the UDP-glucose pyrophosphorylase (Biochimie vol 154, 2018, 176-186), which catalyzes the synthesis of UDP-glucose (the substrate for paramylon synthase). Additionally, the identification of nucleotide sequences coding for putative UDP-sugar pyrophosphorylases suggests the occurrence of an alternative source of UDP-glucose. In this study, we demonstrate the active involvement of both pyrophosphorylases in paramylon accumulation. Using techniques of single and combined knockdown of transcripts coding for these proteins, we evidenced a substantial decrease in the polysaccharide synthesis from 39 ± 7 μg/10⁶ cells determined in the control at day 21st of growth. Thus, the paramylon accumulation in Euglena gracilis cells decreased by 60% and 30% after a single knockdown of the expression of genes coding for UDP-glucose pyrophosphorylase and UDP-sugar pyrophosphorylase, respectively. Besides, the combined knockdown of both genes resulted in a ca. 65% reduction in the level of the storage polysaccharide. Our findings indicate the existence of a physiological dependence between paramylon accumulation and the partitioning of sugar nucleotides into other metabolic routes, including the Leloir pathway's functionality in Euglena gracilis.

The Lipid Composition of Euglena gracilis Middle Plastid Membrane Resembles That of Primary Plastid Envelopes

Euglena gracilis is a photosynthetic flagellate possessing chlorophyte-derived secondary plastids that are enclosed by only three enveloping membranes, unlike most secondary plastids, which are surrounded by four membranes. It has generally been assumed that the two innermost E. gracilis plastid envelopes originated from the primary plastid, while the outermost is of eukaryotic origin. It was suggested that nucleus-encoded plastid proteins pass through the middle and innermost plastid envelopes of E. gracilis by machinery homologous to the translocons of outer and inner chloroplast membranes, respectively. Although recent genomic, transcriptomic, and proteomic data proved the presence of a reduced form of the translocon of inner membrane, they failed to identify any outer-membrane translocon homologs, which raised the question of the origin of E. gracilis's middle plastid envelope. Here, we compared the lipid composition of whole cells of the pigmented E. gracilis strain Z and two bleached mutants that lack detectable plastid structures, W10BSmL and WgmZOflL We determined the lipid composition of E. gracilis strain Z mitochondria and plastids, and of plastid subfractions (thylakoids and envelopes), using HPLC high-resolution tandem mass spectrometry, thin-layer chromatography, and gas chromatography-flame ionization detection analytical techniques. Phosphoglycerolipids are the main structural lipids in mitochondria, while glycosyldiacylglycerols are the major structural lipids of plastids and also predominate in extracts of whole mixotrophic cells. Glycosyldiacylglycerols were detected in both bleached mutants, indicating that mutant cells retain some plastid remnants. Additionally, we discuss the origin of the E. gracilis middle plastid envelope based on the lipid composition of envelope fraction.

The Cryptic Plastid of Euglena longa Defines a New Type of Nonphotosynthetic Plastid Organelle

Most secondary nonphotosynthetic eukaryotes have retained residual plastids whose physiological role is often still unknown. One such example is Euglena longa, a close nonphotosynthetic relative of Euglena gracilis harboring a plastid organelle of enigmatic function. By mining transcriptome data from E. longa, we finally provide an overview of metabolic processes localized to its elusive plastid. The organelle plays no role in the biosynthesis of isoprenoid precursors and fatty acids and has a very limited repertoire of pathways concerning nitrogen-containing metabolites. In contrast, the synthesis of phospholipids and glycolipids has been preserved, curiously with the last step of sulfoquinovosyldiacylglycerol synthesis being catalyzed by the SqdX form of an enzyme so far known only from bacteria. Notably, we show that the E. longa plastid synthesizes tocopherols and a phylloquinone derivative, the first such report for nonphotosynthetic plastids studied so far. The most striking attribute of the organelle could be the presence of a linearized Calvin-Benson (CB) pathway, including RuBisCO yet lacking the gluconeogenetic part of the standard cycle, together with ferredoxin-NADP⁺ reductase (FNR) and the ferredoxin/thioredoxin system. We hypothesize that the ferredoxin/thioredoxin system activates the linear CB pathway in response to the redox status of the E. longa cell and speculate on the role of the pathway in keeping the redox balance of the cell. Altogether, the E. longa plastid defines a new class of relic plastids that is drastically different from the best-studied organelle of this category, the apicoplast.IMPORTANCE Colorless plastids incapable of photosynthesis evolved in many plant and algal groups, but what functions they perform is still unknown in many cases. Here, we study the elusive plastid of Euglena longa, a nonphotosynthetic cousin of the familiar green flagellate Euglena gracilis We document an unprecedented combination of metabolic functions that the E. longa plastid exhibits in comparison with previously characterized nonphotosynthetic plastids. For example, and truly surprisingly, it has retained the synthesis of tocopherols (vitamin E) and a phylloquinone (vitamin K) derivative. In addition, we offer a possible solution of the long-standing conundrum of the presence of the CO₂-fixing enzyme RuBisCO in E. longa Our work provides a detailed account on a unique variant of relic plastids, the first among nonphotosynthetic plastids that evolved by secondary endosymbiosis from a green algal ancestor, and suggests that it has persisted for reasons not previously considered in relation to nonphotosynthetic plastids.

Combined transcriptome sequencing and prokaryotic expression to investigate the key enzyme in the 2-C-methylerythritol-4-phosphate pathway of Osmanthus fragrans

Terpenoids are one of the main components of plant aromas. In the present study we investigated these compounds in Osmanthus fragrans Lour., which is a fragrant plant widely used for the production of essential oils. Quantitative real-time PCR (qRT-PCR) results of enzymes associated with the 2-C-methylerythritol-4-phosphate pathway confirmed that the TPS is a key enzyme for terpenoid synthesis in O. fragrans. In a series of experiments, we identified the TPS candidate genes in O. fragrans and revealed the underlying catalytic activity and subcellular localisation of the encoded proteins. Because there is no available O. fragrans reference genome, we sequenced and analysed its transcriptome and identified two putative TPS genes, OfTPS1 and OfTPS2. According to qRT-PCR analysis, both genes were most highly expressed at the full-bloom stage, suggesting that OfTPS1 and OfTPS2 contribute to O. fragrans terpenoid synthesis. To verify this hypothesis, we constructed prokaryotic expression vectors to obtain protein. In order to study the function of OfTPS1 and OfTPS2 in the synthesis of monoterpenes, the obtained proteins were reacted with geranyl pyrophosphate. As a result, two kinds of monoterpenes, (E)-β-ocimene and linalool, were detected from reaction products, respectively. In conclusion, OfTPS1 and OfTPS2 are both monoterpene synthases.

A Uniquely Complex Mitochondrial Proteome from Euglena gracilis

Euglena gracilis is a metabolically flexible, photosynthetic, and adaptable free-living protist of considerable environmental importance and biotechnological value. By label-free liquid chromatography tandem mass spectrometry, a total of 1,786 proteins were identified from the E. gracilis purified mitochondria, representing one of the largest mitochondrial proteomes so far described. Despite this apparent complexity, protein machinery responsible for the extensive RNA editing, splicing, and processing in the sister clades diplonemids and kinetoplastids is absent. This strongly suggests that the complex mechanisms of mitochondrial gene expression in diplonemids and kinetoplastids occurred late in euglenozoan evolution, arising independently. By contrast, the alternative oxidase pathway and numerous ribosomal subunits presumed to be specific for parasitic trypanosomes are present in E. gracilis. We investigated the evolution of unexplored protein families, including import complexes, cristae formation proteins, and translation termination factors, as well as canonical and unique metabolic pathways. We additionally compare this mitoproteome with the transcriptome of Eutreptiella gymnastica, illuminating conserved features of Euglenida mitochondria as well as those exclusive to E. gracilis. This is the first mitochondrial proteome of a free-living protist from the Excavata and one of few available for protists as a whole. This study alters our views of the evolution of the mitochondrion and indicates early emergence of complexity within euglenozoan mitochondria, independent of parasitism.

Light dependent accumulation of β-carotene enhances photo-acclimation of Euglena gracilis

Carotenoids are essential components of photosynthetic organisms including land plants, algae, cyanobacteria, and photosynthetic bacteria. Although the light-mediated regulation of carotenoid biosynthesis, including the light/dark cycle as well as the dependence of carotenoid biosynthesis-related gene translation on light wavelength, has been investigated in land plants, these aspects have not been studied in microalgae. Here, we investigated carotenoid biosynthesis in Euglena gracilis and found that zeaxanthin accumulates in the dark. The major carotenoid species in E. gracilis, namely β-carotene, neoxanthin, diadinoxanthin and diatoxanthin, accumulated corresponding to the duration of light irradiation under the light/dark cycle, although the translation of carotenoid biosynthesis genes hardly changed. Irradiation with either blue or red-light (3 μmol photons m^-2 s^-1) caused a 1.3-fold increase in β-carotene content compared with the dark control. Blue-light irradiation (300 μmol photons m^-2 s^-1) caused an increase in the cellular content of both zeaxanthin and all trans-diatoxanthin, and this increase was proportional to blue-light intensity. In addition, pre-irradiation with blue-light of 3 or 30 μmol photons m^-2 s^-1 enhanced the photosynthetic activity and tolerance to high-light stress. These findings suggest that the accumulation of β-carotene is regulated by the intensity of light, which may contribute to the acclimation of E. gracilis to the light environment in day night conditions.

A Global Analysis of Enzyme Compartmentalization to Glycosomes

In kinetoplastids, the first seven steps of glycolysis are compartmentalized into a glycosome along with parts of other metabolic pathways. This organelle shares a common ancestor with the better-understood eukaryotic peroxisome. Much of our understanding of the emergence, evolution, and maintenance of glycosomes is limited to explorations of the dixenous parasites, including the enzymatic contents of the organelle. Our objective was to determine the extent that we could leverage existing studies in model kinetoplastids to determine the composition of glycosomes in species lacking evidence of experimental localization. These include diverse monoxenous species and dixenous species with very different hosts. For many of these, genome or transcriptome sequences are available. Our approach initiated with a meta-analysis of existing studies to generate a subset of enzymes with highest evidence of glycosome localization. From this dataset we extracted the best possible glycosome signal peptide identification scheme for in silico identification of glycosomal proteins from any kinetoplastid species. Validation suggested that a high glycosome localization score from our algorithm would be indicative of a glycosomal protein. We found that while metabolic pathways were consistently represented across kinetoplastids, individual proteins within those pathways may not universally exhibit evidence of glycosome localization.

Evolution of metabolic capabilities and molecular features of diplonemids, kinetoplastids, and euglenids

Background

The Euglenozoa are a protist group with an especially rich history of evolutionary diversity. They include diplonemids, representing arguably the most species-rich clade of marine planktonic eukaryotes; trypanosomatids, which are notorious parasites of medical and veterinary importance; and free-living euglenids. These different lifestyles, and particularly the transition from free-living to parasitic, likely require different metabolic capabilities. We carried out a comparative genomic analysis across euglenozoan diversity to see how changing repertoires of enzymes and structural features correspond to major changes in lifestyles.

Results

We find a gradual loss of genes encoding enzymes in the evolution of kinetoplastids, rather than a sudden decrease in metabolic capabilities corresponding to the origin of parasitism, while diplonemids and euglenids maintain more metabolic versatility. Distinctive characteristics of molecular machines such as kinetochores and the pre-replication complex that were previously considered specific to parasitic kinetoplastids were also identified in their free-living relatives. Therefore, we argue that they represent an ancestral rather than a derived state, as thought until the present. We also found evidence of ancient redundancy in systems such as NADPH-dependent thiol-redox. Only the genus Euglena possesses the combination of trypanothione-, glutathione-, and thioredoxin-based systems supposedly present in the euglenozoan common ancestor, while other representatives of the phylum have lost one or two of these systems. Lastly, we identified convergent losses of specific metabolic capabilities between free-living kinetoplastids and ciliates. Although this observation requires further examination, it suggests that certain eukaryotic lineages are predisposed to such convergent losses of key enzymes or whole pathways.

Conclusions

The loss of metabolic capabilities might not be associated with the switch to parasitic lifestyle in kinetoplastids, and the presence of a highly divergent (or unconventional) kinetochore machinery might not be restricted to this protist group. The data derived from the transcriptomes of free-living early branching prokinetoplastids suggests that the pre-replication complex of Trypanosomatidae is a highly divergent version of the conventional machinery. Our findings shed light on trends in the evolution of metabolism in protists in general and open multiple avenues for future research.

Oxygenic Phototrophs Need ζ-Carotene Isomerase (Z-ISO) for Carotene Synthesis: Functional Analysis in Arthrospira and Euglena

For carotenogenesis, two biosynthetic pathways from phytoene to lycopene are known. Most bacteria and fungi require only phytoene desaturase (PDS, CrtI), whereas land plants require four enzymes: PDS (CrtP), ζ-carotene desaturase (ZDS, CrtQ), ζ-carotene isomerase (Z-ISO) and cis-carotene isomerase (CrtISO, CrtH). The gene encoding Z-ISO has been functionally identified in only two species, Arabidopsis thaliana and Zea mays, and has been little studied in other organisms. In this study, we found that the deduced amino acid sequences of Arthrospira Z-ISO and Euglena Z-ISO have 58% and 62% identity, respectively, with functional Z-ISO from Arabidopsis. We studied the function of Z-ISO genes from the cyanobacterium Arthrospira platensis and eukaryotic microalga Euglena gracilis. The Z-ISO genes of Arthrospira and Euglena were transformed into Escherichia coli strains that produced mainly 9,15,9'-tri-cis-ζ-carotene in darkness. In the resulting E. coli transformants cultured under darkness, 9,9'-di-cis-ζ-carotene was accumulated predominantly as Z-ISO in Arabidopsis. This indicates that the Z-ISO genes were involved in the isomerization of 9,15,9'-tri-cis-ζ-carotene to 9,9'-di-cis-ζ-carotene in darkness. This is the first functional analysis of Z-ISO as a ζ-carotene isomerase in cyanobacteria and eukaryotic microalgae. Green sulfur bacteria and Chloracidobacterium also use CrtP, CrtQ and CrtH for lycopene synthesis as cyanobacteria, but their genomes did not comprise Z-ISO genes. Consequently, Z-ISO is needed in oxygenic phototrophs, whereas it is not found in anoxygenic species.

Metabolic quirks and the colourful history of the Euglena gracilis secondary plastid

Euglena spp. are phototrophic flagellates with considerable ecological presence and impact. Euglena gracilis harbours secondary green plastids, but an incompletely characterised proteome precludes accurate understanding of both plastid function and evolutionary history. Using subcellular fractionation, an improved sequence database and MS we determined the composition, evolutionary relationships and hence predicted functions of the E. gracilis plastid proteome. We confidently identified 1345 distinct plastid protein groups and found that at least 100 proteins represent horizontal acquisitions from organisms other than green algae or prokaryotes. Metabolic reconstruction confirmed previously studied/predicted enzymes/pathways and provided evidence for multiple unusual features, including uncoupling of carotenoid and phytol metabolism, a limited role in amino acid metabolism, and dual sets of the SUF pathway for FeS cluster assembly, one of which was acquired by lateral gene transfer from Chlamydiae. Plastid paralogues of trafficking-associated proteins potentially mediating fusion of transport vesicles with the outermost plastid membrane were identified, together with derlin-related proteins, potential translocases across the middle membrane, and an extremely simplified TIC complex. The Euglena plastid, as the product of many genomes, combines novel and conserved features of metabolism and transport.

Comparative proteomic analysis of mitochondria isolated from Euglena gracilis under aerobic and hypoxic conditions

The unicellular microalga Euglena gracilis produces wax esters for ATP acquisition under low-oxygen conditions. The regulatory mechanism of wax ester production is not yet understood. Indeed, our previous transcriptomic analysis showed that transcript levels of genes involved in the wax ester synthesis hardly changed under hypoxic conditions, suggesting contribution of post-transcriptional regulation. In this study, we conducted a proteome analysis of E. gracilis mitochondria, as this organelle employs the fatty-acid synthesis pathway under hypoxic conditions. Mitochondria were isolated from E. gracilis SM-ZK strain treated with both aerobic and hypoxic conditions and used for shotgun proteomic analysis. Three independent proteomic analyses succeeded in identifying a total of 714 non-redundant proteins. Of these, 229 were detected in common to all experiments, and 116 were significantly recognized as differentially expressed proteins. GO enrichment analysis suggested dynamic changes in mitochondrial metabolic pathways and redox reactions under aerobic and hypoxic conditions. Protein levels of bifunctional enzymes isocitrate lyase and malate synthase in glyoxylate cycle were 1.35-fold higher under hypoxic conditions. Abundances of the propionyl-CoA synthetic enzymes, succinyl-CoA synthetase and propionyl-CoA carboxylase, were also 1.35- and 1.47-fold higher, respectively, under hypoxic conditions. Protein levels of pyruvate:NADP⁺ oxidoreductase, a key enzyme for anaerobic synthesis of acetyl-CoA, which serves as a C2 donor for fatty acids, showed a 1.68-fold increase under hypoxic conditions, whereas those of pyruvate dehydrogenase subunits showed a 0.77–0.81-fold decrease. Protein levels of the fatty-acid synthesis enzymes, 3-ketoacyl-CoA thiolase isoforms (KAT1 and KAT2), 3-hydroxyacyl-CoA dehydrogenases, and acyl-CoA dehydrogenase were up-regulated by 1.20- to 1.42-fold in response to hypoxic treatment. Overall, our proteomic analysis revealed that wax ester synthesis-related enzymes are up-regulated at the protein level post-transcriptionally to promote wax ester production in E. gracilis under low-oxygen conditions.

Changes of Gene Expression in Euglena gracilis Obtained During the 29th DLR Parabolic Flight Campaign

Parabolic flight maneuvers of Novespace's Airbus A310 ZERO-G produce subsequent phases of hypergravity (about 20 s), microgravity (about 22 s) and another 20 s hypergravity on experiments located in the experiment area of the aircraft. The 29th DLR parabolic flight campaign consisted of four consecutive flight days with thirty-one parabolas each day. Euglena gracilis cells were fixed with TRIzol during different acceleration conditions at the first and the last parabola of each flight. Samples were collected and analyzed with microarrays for one-color gene expression analysis. The data indicate significant changes in gene expression in E. gracilis within short time. Hierarchical clustering shows that changes induced by the different accelerations yield reproducible effects at independent flight days. Transcription differed between the first and last parabolas indicating adaptation effects in the course of the flight. Different gene groups were found to be affected in different phases of the parabolic flight, among others, genes involved in signal transduction, calcium signaling, transport mechanisms, metabolic pathways, and stress-response as well as membrane and cytoskeletal proteins. In addition, transcripts of other areas, e.g., DNA and protein modification, were altered. The study contributes to the understanding of short-term effects of microgravity and different accelerations on cells at a molecular level.

Euglena Central Metabolic Pathways and Their Subcellular Locations

Euglenids are a group of algae of great interest for biotechnology, with a large and complex metabolic capability. To study the metabolic network, it is necessary to know where the component enzymes are in the cell, but despite a long history of research into Euglena, the subcellular locations of many major pathways are only poorly defined. Euglena is phylogenetically distant from other commonly studied algae, they have secondary plastids bounded by three membranes, and they can survive after destruction of their plastids. These unusual features make it difficult to assume that the subcellular organization of the metabolic network will be equivalent to that of other photosynthetic organisms. We analysed bioinformatic, biochemical, and proteomic information from a variety of sources to assess the subcellular location of the enzymes of the central metabolic pathways, and we use these assignments to propose a model of the metabolic network of Euglena. Other than photosynthesis, all major pathways present in the chloroplast are also present elsewhere in the cell. Our model demonstrates how Euglena can synthesise all the metabolites required for growth from simple carbon inputs, and can survive in the absence of chloroplasts.

Comparative molecular cell biology of phototrophic euglenids and parasitic trypanosomatids sheds light on the ancestor of Euglenozoa

Parasitic trypanosomatids and phototrophic euglenids are among the most extensively studied euglenozoans. The phototrophic euglenid lineage arose relatively recently through secondary endosymbiosis between a phagotrophic euglenid and a prasinophyte green alga that evolved into the euglenid secondary chloroplast. The parasitic trypanosomatids (i.e. Trypanosoma spp. and Leishmania spp.) and the freshwater phototrophic euglenids (i.e. Euglena gracilis) are the most evolutionary distant lineages in the Euglenozoa phylogenetic tree. The molecular and cell biological traits they share can thus be considered as ancestral traits originating in the common euglenozoan ancestor. These euglenozoan ancestral traits include common mitochondrial presequence motifs, respiratory chain complexes containing various unique subunits, a unique ATP synthase structure, the absence of mitochondria-encoded transfer RNAs (tRNAs), a nucleus with a centrally positioned nucleolus, closed mitosis without dissolution of the nuclear membrane and nucleoli, a nuclear genome containing the unusual 'J' base (β-D-glucosyl-hydroxymethyluracil), processing of nucleus-encoded precursor messenger RNAs (pre-mRNAs) via spliced-leader RNA (SL-RNA) trans-splicing, post-transcriptional gene silencing by the RNA interference (RNAi) pathway and the absence of transcriptional regulation of nuclear gene expression. Mitochondrial uridine insertion/deletion RNA editing directed by guide RNAs (gRNAs) evolved in the ancestor of the kinetoplastid lineage. The evolutionary origin of other molecular features known to be present only in either kinetoplastids (i.e. polycistronic transcripts, compaction of nuclear genomes) or euglenids (i.e. monocistronic transcripts, huge genomes, many nuclear cis-spliced introns, polyproteins) is unclear.

Bioproducts From Euglena gracilis: Synthesis and Applications

In recent years, the versatile phototrophic protist Euglena gracilis has emerged as an interesting candidate for application-driven research and commercialisation, as it is an excellent source of dietary protein, pro(vitamins), lipids, and the β-1,3-glucan paramylon only found in euglenoids. From these, paramylon is already marketed as an immunostimulatory agent in nutraceuticals. Bioproducts from E. gracilis can be produced under various cultivation conditions discussed in this review, and their yields are relatively high when compared with those achieved in microalgal systems. Future challenges include achieving the economy of large-scale cultivation. Recent insights into the complex metabolism of E. gracilis have highlighted unique metabolic pathways, which could provide new leads for product enhancement by genetic modification of the organism. Also, development of molecular tools for strain improvement are emerging rapidly, making E. gracilis a noteworthy challenger for microalgae such as Chlorella spp. and their products currently on the market.

De Novo Assembly and Comparative Transcriptome Profiling of Anser anser and Anser cygnoides Geese Species' Embryonic Skin Feather Follicles

Geese feather production and the quality of downy feathers are additional economically important traits in the geese industry. However, little information is available about the molecular mechanisms fundamental to feather formation and the quality of feathers in geese. This study conducted de novo transcriptome sequencing analysis of two related geese species using the Illumina 4000 platform to determine the genes involved in embryonic skin feather follicle development. A total of 165,564,278 for Anser anser and 144,595,262 for Anser cygnoides clean reads were generated, which were further assembled into 77,134 unigenes with an average length of 906 base pairs in Anser anser and 66,041 unigenes with an average length of 922 base pairs in Anser cygnoides. To recognize the potential regulatory roles of differentially expressed genes (DEGs) during geese embryonic skin feather follicle development, the obtained unigenes were annotated to Gene Ontology (GO), Eukaryotic Orthologous Groups (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) for functional analysis. In both species, GO and KOG had shown similar distribution patterns during functional annotation except for KEGG, which showed significant variation in signaling enrichment. Anser asnser was significantly enriched in the calcium signaling pathway, whereas Anser cygnoides was significantly enriched with glycerolipid metabolism. Further analysis indicated that 14,227 gene families were conserved between the species, among which a total of 20,715 specific gene families were identified. Comparative RNA-Seq data analysis may reveal inclusive knowledge to assist in the identification of genetic regulators at a molecular level to improve feather quality production in geese and other poultry species.

Transcriptome, proteome and draft genome of Euglena gracilis

BACKGROUND

Photosynthetic euglenids are major contributors to fresh water ecosystems. Euglena gracilis in particular has noted metabolic flexibility, reflected by an ability to thrive in a range of harsh environments. E. gracilis has been a popular model organism and of considerable biotechnological interest, but the absence of a gene catalogue has hampered both basic research and translational efforts.

RESULTS

We report a detailed transcriptome and partial genome for E. gracilis Z1. The nuclear genome is estimated to be around 500 Mb in size, and the transcriptome encodes over 36,000 proteins and the genome possesses less than 1% coding sequence. Annotation of coding sequences indicates a highly sophisticated endomembrane system, RNA processing mechanisms and nuclear genome contributions from several photosynthetic lineages. Multiple gene families, including likely signal transduction components, have been massively expanded. Alterations in protein abundance are controlled post-transcriptionally between light and dark conditions, surprisingly similar to trypanosomatids.

CONCLUSIONS

Our data provide evidence that a range of photosynthetic eukaryotes contributed to the Euglena nuclear genome, evidence in support of the 'shopping bag' hypothesis for plastid acquisition. We also suggest that euglenids possess unique regulatory mechanisms for achieving extreme adaptability, through mechanisms of paralog expansion and gene acquisition.

Low Temperature Stress Alters the Expression of Phytoene Desaturase Genes (crtP1 and crtP2) and the ζ-Carotene Desaturase Gene (crtQ) Together with the Cellular Carotenoid Content of Euglena gracilis

Carotenoids participate in photosynthesis and photoprotection in oxygenic phototrophs. Euglena gracilis, a eukaryotic phytoflagellate, synthesizes several carotenoids: β-carotene, neoxanthin, diadinoxanthin and diatoxanthin. Temperature is one of the most striking external stimuli altering carotenoid production. In the present study, to elucidate the regulation of carotenoid synthesis of E. gracilis in response to environmental stimuli, we functionally identified phytoene desaturase genes (crtP1 and crtP2) and the ζ-carotene desaturase gene (crtQ) of this alga and analyzed expression of those genes and the composition of major carotenoids in cells grown under cold (20�C) and high-intensity light (HL; 240 �mol photon m-2 s-1) conditions. 20�C-HL treatment increased the transcriptional level of the phytoene synthase gene (crtB), and crtP1 and crtP2, whose products catalyze the early steps of carotenoid biosynthesis in this alga. Cultivation at 20�C under illumination at 55 �mol photon m-2 s-1 (low-intensity light; LL) decreased the cell concentration, Chl and total major carotenoid content by 61, 75 and 50%, respectively, relative to control (25�C-LL) cells. When grown at 20�C-HL, the cells showed a greater decrease in cell concentration and photosynthetic pigment contents than those in 20�C-LL. β-Carotene, neoxanthin and diadinoxanthin contents were decreased by more than half in 20�C-LL and 20�C-HL treatments. On the other hand, when subjected to 20�C-LL and 20�C-HL, the cells retained a diatoxanthin content comparable with control cells. Our findings suggested that diatoxanthin plays crucial roles in the acclimation to cold and intense light condition. To the best of our knowledge, this is the first report on a photosynthetic organism possessing dual crtP genes.

Characterization of sulfur-compound metabolism underlying wax-ester fermentation in Euglena gracilis

Euglena gracilis is a microalga, which has been used as a model organism for decades. Recent technological advances have enabled mass cultivation of this species for industrial applications such as feedstock in nutritional foods and cosmetics. E. gracilis degrades its storage polysaccharide (paramylon) under hypoxic conditions for energy acquisition by an oxygen-independent process and accumulates high amount of wax-ester as a by-product. Using this sequence of reactions referred to as wax-ester fermentation, E. gracilis is studied for its application in biofuel production. Although the wax-ester production pathway is well characterized, little is known regarding the biochemical reactions underlying the main metabolic route, especially, the existence of an unknown sulfur-compound metabolism implied by the nasty odor generation accompanying the wax-ester fermentation. In this study, we show sulfur-metabolomics of E. gracilis in aerobic and hypoxic conditions, to reveal the biochemical reactions that occur during wax-ester synthesis. Our results helped us in identifying hydrogen sulfide (H2S) as the nasty odor-producing component in wax-ester fermentation. In addition, the results indicate that glutathione and protein degrades during hypoxia, whereas cysteine, methionine, and their metabolites increase in the cells. This indicates that this shift of abundance in sulfur compounds is the cause of H2S synthesis.

A major isoform of mitochondrial trans-2-enoyl-CoA reductase is dispensable for wax ester production in Euglena gracilis under anaerobic conditions

Under anaerobic conditions, Euglena gracilis produces a large amount of wax ester through mitochondrial fatty acid synthesis from storage polysaccharides termed paramylon, to generate ATP. Trans-2-enoyl-CoA reductases (TERs) in mitochondria have been considered to play a key role in this process, because the enzymes catalyze the reduction of short chain length CoA-substrates (such as crotonyl-CoA). A TER enzyme (EgTER1) has been previously identified and enzymologically characterized; however, its physiological significance remained to be evaluated by genetic analysis. We herein generated EgTER1-knockdown Euglena cells, in which total crotonyl-CoA reductase activity was decreased to 10% of control value. Notably, the knockdown cells showed a severe bleaching phenotype with deficiencies in chlorophylls and glycolipids, but grew normally under heterotrophic conditions (with glucose supplementation). Moreover, the knockdown cells accumulated much greater quantities of wax ester than control cells before and after transfer to anaerobic conditions, which was accompanied by a large metabolomic change. Furthermore, we failed to find any contribution of other potential TER genes in wax ester production. Our findings propose a novel role of EgTER1 in the greening process and demonstrate that this enzyme is dispensable for wax ester production under anaerobic conditions.

Physiological role of β-carotene monohydroxylase (CYP97H1) in carotenoid biosynthesis in Euglena gracilis

Some carotenoids are found in the Euglena gracilis, including β-carotene, diadinoxanthin, diatoxanthins, and neoxanthin as the major species; however, the molecular mechanism underlying carotenoid biosynthesis in E. gracilis is not well understood. To clarify the pathway and regulation of carotenoid biosynthesis in this alga, we functionally characterized the cytochrome P450 (CYP)-type carotene hydroxylase gene EgCYP97H1. Heterologous in vivo enzyme assay in E. coli indicated that EgCYP97H1 hydroxylated β-carotene to β-cryptoxanthin. E. gracilis cells suppressing EgCYP97H1 resulted in marked growth inhibition and reductions in total carotenoid and chlorophyll contents. Analysis of carotenoid composition revealed that suppression of EgCYP97H1 resulted in higher level of β-carotene, suggesting that EgCYP97H1 is physiologically essential for carotenoid biosynthesis and thus normal cell growth. To our knowledge, this is the first time EgCYP97H1 has been suggested to be β-carotene monohydroxylase, but not β-carotene dihydroxylase. Moreover, during light adaptation of dark-grown E. gracilis, transcript levels of the carotenoid biosynthetic genes (EgCYP97H1, geranylgeranyl pyrophosphate synthase EgcrtE, and phytoene synthase EgcrtB) remained virtually unchanged. In contrast, carotenoid accumulation in E. gracilis grown under the same conditions was inhibited by treatment with a translational inhibitor but not a transcriptional inhibitor, indicating that photo-responsive carotenoid biosynthesis is regulated post-transcriptionally in this alga.

Anaerobic respiration coupled with mitochondrial fatty acid synthesis in wax ester fermentation by Euglena gracilis

In Euglena gracilis, wax ester fermentation produces ATP during anaerobiosis. Here, we report that anaerobic wax ester production is suppressed when the mitochondrial electron transport chain complex I is inhibited by rotenone, whereas it is increased by the uncoupler carbonyl cyanide m-chlorophenylhydrazone (CCCP). The ADP/ATP ratio in anaerobic cells is elevated by treatment with either rotenone or CCCP. Gene silencing experiments indicate that acyl-CoA dehydrogenase, electron transfer flavoprotein (ETF), and rhodoquinone (RQ) participate in wax ester production. These results suggest that fatty acids are synthesized in mitochondria by the reversal of β-oxidation, where trans-2-enoyl-CoA is reduced mainly by acyl-CoA dehydrogenase using the electrons provided by NADH via the electron transport chain complex I, RQ, and ETF, and that ATP production is highly supported by anaerobic respiration utilizing trans-2-enoyl-CoA as a terminal electron acceptor.

Peculiar features of the plastids of the colourless alga Euglena longa and photosynthetic euglenophytes unveiled by transcriptome analyses

Euglenophytes are a familiar algal group with green alga-derived secondary plastids, but the knowledge of euglenophyte plastid function and evolution is still highly incomplete. With this in mind we sequenced and analysed the transcriptome of the non-photosynthetic species Euglena longa. The transcriptomic data confirmed the absence of genes for the photosynthetic machinery, but provided candidate plastid-localised proteins bearing N-terminal bipartite topogenic signals (BTSs) of the characteristic euglenophyte type. Further comparative analyses including transcriptome assemblies available for photosynthetic euglenophytes enabled us to unveil salient aspects of the basic euglenophyte plastid infrastructure, such as plastidial targeting of several proteins as C-terminal translational fusions with other BTS-bearing proteins or replacement of the conventional eubacteria-derived plastidial ribosomal protein L24 by homologs of archaeo-eukaryotic origin. Strikingly, no homologs of any key component of the TOC/TIC system and the plastid division apparatus are discernible in euglenophytes, and the machinery for intraplastidial protein targeting has been simplified by the loss of the cpSRP/cpFtsY system and the SEC2 translocon. Lastly, euglenophytes proved to encode a plastid-targeted homolog of the termination factor Rho horizontally acquired from a Lambdaproteobacteria-related donor. Our study thus further documents a substantial remodelling of the euglenophyte plastid compared to its green algal progenitor.

RNA-Seq employing a novel rRNA depletion strategy reveals a rich repertoire of snoRNAs in Euglena gracilis including box C/D and Ψ-guide RNAs targeting the modification of rRNA extremities

Previous mRNA transcriptome studies of Euglena gracilis have shown that this organism possesses a large and diverse complement of protein coding genes; however, the study of non-coding RNA classes has been limited. The natural extensive fragmentation of the E. gracilis large subunit ribosomal RNA presents additional barriers to the identification of non-coding RNAs as size-selected small RNA libraries will be dominated by rRNA sequences. In this study we have developed a strategy to significantly reduce rRNA amplification prior to RNA-Seq analysis thereby producing a ncRNA library allowing for the identification of many new E. gracilis small RNAs. Library analysis reveals 113 unique new small nucleolar (sno) RNAs and a large collection of snoRNA isoforms, as well as the first significant collection of nuclear tRNAs in this organism. A 3' end AGAUGN consensus motif and conserved structural features can now be defined for E. gracilis pseudouridine guide RNAs. snoRNAs of both classes were identified that target modification of the 3' extremities of rRNAs utilizing predicted base-pairing interactions with internally transcribed spacers (ITS), providing insight into the timing of steps in rRNA maturation. Cumulatively, this represents the most comprehensive analysis of small ncRNAs in Euglena gracilis to date.

Integrated analysis of transcription factors and targets co-expression profiles reveals reduced correlation between transcription factors and target genes in cancer

Transcription factors are recognized as the key regulators of gene expression. However, the changes in the correlation of transcription factors and their target genes between normal and tumor tissues are usually ignored. In this research, we used mRNA expression profile data from The Cancer Genome Atlas which included 5726 samples across 11 major human cancers to perform co-expression analysis by the Pearson correlation coefficients. Then, integrating 81,357 pairs of transcription factors and target genes from transcription factors databases to find out the changes in the co-expression correlation of these gene pairs from normal to tumor tissues. Based on the changes in the number of co-expressed TF-TG pairs and changes in the level of co-expression, we found the generally reduced correlation between transcription factors and their target genes in cancer. Additionally, we screened out universal and specific transcription factors-target genes pairs which may significant influence particular cancer. Then, we obtained 423 cancer cell line expression profiles from Broad Institute Cancer Cell Line Encyclopedia to verify our results. Some of these pairs like XRCC5-XRCC6 have been reported to involve in multiple cancers, while pairs like IRF1-PSMB9 without any previous articles related to tumor but involve in the biological processes of cancer, which are of great potential to be therapeutic targets. Our research may provide insights to better understand the tumor development mechanisms and find potential therapeutic targets.

Comparative assessment of the Euglena gracilis var. saccharophila variant strain as a producer of the β-1,3-glucan paramylon under varying light conditions

Euglena gracilis Z and a "sugar loving" variant strain E. gracilis var. saccharophila were investigated as producers of paramylon, a β-1,3-glucan polysaccharide with potential medicinal and industrial applications. The strains were grown under diurnal or dark growth conditions on a glucose-yeast extract medium supporting high-level paramylon production. Both strains produced the highest paramylon yields (7.4-8 g · L^-1 , respectively) while grown in the dark, but the maximum yield was achieved faster by E. gracilis var. saccharophila (48 h vs. 72 h). The glucose-to-paramylon yield coefficient Y_par/glu  = 0.46 ± 0.03 in the E. gracilis var. saccharophila cultivation, obtained in this study, is the highest reported to date. Proteomic analysis of the metabolic pathways provided molecular clues for the strain behavior observed during cultivation. For example, overexpression of enzymes in the gluconeogenesis/glycolysis pathways including fructokinase-1 and chloroplastic fructose-1,6-bisphosphatase (FBP) may have contributed to the faster rate of paramylon accumulation in E. gracilis var. saccharophila. Differentially expressed proteins in the early steps of chloroplastogenesis pathway including plastid uroporphyrinogen decarboxylases, photoreceptors, and a highly abundant (68-fold increase) plastid transketolase may have provided the E. gracilis var. saccharophila strain an advantage in paramylon production during diurnal cultivations. In conclusion, the variant strain E. gracilis var. saccharophila seems to be well suited for producing large amounts of paramylon. This work has also resulted in the identification of molecular targets for future improvement of paramylon production in E. gracilis, including the FBP and phosophofructokinase 1, the latter being a key regulator of glycolysis.

Alterations of Membrane Lipid Content Correlated With Chloroplast and Mitochondria Development in Euglena gracilis

Euglenoids are unique protists that can grow photoautotrophically, photomixotrophically, and heterotrophically. Here we grew Euglena gracilis under these different growth conditions and determined cellular contents of seven membrane lipids and one storage lipid (triacylglycerol), which account for more than 94 mol% of total membrane lipids. We also describe the relationship among chloroplast and mitochondria developments with lipid contents, protein contents, and oxygen evolution/consumption rates. In photoautotrophic growth conditions, E. gracilis cells accumulated chlorophyll, photosynthetic proteins, and glycolipids typical to thylakoid membranes. The same occurred for the cells grown under photomixotrophic conditions with higher respiration rates. In heterotrophic conditions, E. gracilis cells had higher respiration rates compared to cells grown in other conditions with the accumulation of pyruvate: NADP+ oxidoreductase, a mitochondrial protein and phospholipid common in mitochondria. Cells were also observed using a confocal laser scanning microscope and found to show more chlorophyll autofluorescence when grown photoautotrophically and photomixotrophycally, and fluorescence of MitoTracker when grown photomixotrophically and heterotrophically. These results suggest that under illumination, E. gracilis develops functional thylakoid membranes with membrane lipids and proteins for photosynthesis. In the medium with glucose, the cells develop mitochondria with phospholipids and proteins for respiration. Possible application based on lipid analysis for the enhancement of wax ester or alkene synthesis is discussed.

Identification of Euglena gracilis β-1,3-glucan phosphorylase and establishment of a new glycoside hydrolase (GH) family GH149

Glycoside phosphorylases (EC 2.4.x.x) carry out the reversible phosphorolysis of glucan polymers, producing the corresponding sugar 1-phosphate and a shortened glycan chain. β-1,3-Glucan phosphorylase activities have been reported in the photosynthetic euglenozoan Euglena gracilis, but the cognate protein sequences have not been identified to date. Continuing our efforts to understand the glycobiology of E. gracilis, we identified a candidate phosphorylase sequence, designated EgP1, by proteomic analysis of an enriched cellular protein lysate. We expressed recombinant EgP1 in Escherichia coli and characterized it in vitro as a β-1,3-glucan phosphorylase. BLASTP identified several hundred EgP1 orthologs, most of which were from Gram-negative bacteria and had 37-91% sequence identity to EgP1. We heterologously expressed a bacterial metagenomic sequence, Pro_7066 in E. coli and confirmed it as a β-1,3-glucan phosphorylase, albeit with kinetics parameters distinct from those of EgP1. EgP1, Pro_7066, and their orthologs are classified as a new glycoside hydrolase (GH) family, designated GH149. Comparisons between GH94, EgP1, and Pro_7066 sequences revealed conservation of key amino acids required for the phosphorylase activity, suggesting a phosphorylase mechanism that is conserved between GH94 and GH149. We found bacterial GH149 genes in gene clusters containing sugar transporter and several other GH family genes, suggesting that bacterial GH149 proteins have roles in the degradation of complex carbohydrates. The Bacteroidetes GH149 genes located to previously identified polysaccharide utilization loci, implicated in the degradation of complex carbohydrates. In summary, we have identified a eukaryotic and a bacterial β-1,3-glucan phosphorylase and uncovered a new family of phosphorylases that we name GH149.

The Mitochondrion of Euglena gracilis

In the presence of oxygen, Euglena gracilis mitochondria function much like mammalian mitochondria. Under anaerobiosis, E. gracilis mitochondria perform a malonyl-CoA independent synthesis of fatty acids leading to accumulation of wax esters, which serve as the sink for electrons stemming from glycolytic ATP synthesis and pyruvate oxidation. Some components (enzymes and cofactors) of Euglena's anaerobic energy metabolism are found among the anaerobic mitochondria of invertebrates, others are found among hydrogenosomes, the H₂-producing anaerobic mitochondria of protists.