Transcriptome, proteome and draft genome of Euglena gracilis

A lineage-specific protein network at the trypanosome nuclear envelope

The nuclear envelope (NE) separates translation and transcription and is the location of multiple functions, including chromatin organization and nucleocytoplasmic transport. The molecular basis for many of these functions have diverged between eukaryotic lineages. Trypanosoma brucei, a member of the early branching eukaryotic lineage Discoba, highlights many of these, including a distinct lamina and kinetochore composition. Here, we describe a cohort of proteins interacting with both the lamina and NPC, which we term lamina-associated proteins (LAPs). LAPs represent a diverse group of proteins, including two candidate NPC-anchoring pore membrane proteins (POMs) with architecture conserved with S. cerevisiae and H. sapiens, and additional peripheral components of the NPC. While many of the LAPs are Kinetoplastid specific, we also identified broadly conserved proteins, indicating an amalgam of divergence and conservation within the trypanosome NE proteome, highlighting the diversity of nuclear biology across the eukaryotes, increasing our understanding of eukaryotic and NPC evolution.

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

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.

The biomolecules of Euglena gracilis: Harnessing biology for natural solutions to future problems

Over the past decade, the autotrophic and heterotrophic protist Euglena gracilis (E. gracilis) has gained popularity across the studies of environmental science, biosynthesis experiments, and nutritional substitutes. The unique physiology and versatile metabolism of E. gracilis have been a recent topic of interest to many researchers who continue to understand the complexity and possibilities of using E. gracilis biomolecule production. In this review, we present a comprehensive representation of recent literature outlining the various uses of biomolecules derived from E. gracilis across the fields of natural product biosynthesis, as a nutritional substitute, and as bioremediation tools. In addition, we highlight effective strategies for altering metabolite production using abiotic stressors and growth conditions. To better understand metabolite biosynthesis and its role in E. gracilis, integrated studies involving genomics, metabolomics, and proteomics should be considered. Together, we show how the ongoing advancements in E. gracilis related research continue to broaden applications in the biosynthetic sector and highlight future works that would strengthen our understanding of overall Euglena metabolism.

Regulation of Microalgal Photosynthetic Electron Transfer

The global ecosystem relies on the metabolism of photosynthetic organisms, featuring the ability to harness light as an energy source. The most successful type of photosynthesis utilizes a virtually inexhaustible electron pool from water, but the driver of this oxidation, sunlight, varies on time and intensity scales of several orders of magnitude. Such rapid and steep changes in energy availability are potentially devastating for biological systems. To enable a safe and efficient light-harnessing process, photosynthetic organisms tune their light capturing, the redox connections between core complexes and auxiliary electron mediators, ion passages across the membrane, and functional coupling of energy transducing organelles. Here, microalgal species are the most diverse group, featuring both unique environmental adjustment strategies and ubiquitous protective mechanisms. In this review, we explore a selection of regulatory processes of the microalgal photosynthetic apparatus supporting smooth electron flow in variable environments.

A chromosome-level genome assembly for the paramylon-producing microalga Euglena gracilis

Euglena gracilis (E. gracilis), pivotal in the study of photosynthesis, endosymbiosis, and chloroplast development, is also an industrial microalga for paramylon production. Despite its importance, E. gracilis genome exploration faces challenges due to its intricate nature. In this study, we achieved a chromosome-level de novo assembly (2.37 Gb) using Illumina, PacBio, Bionano, and Hi-C data. The assembly exhibited a contig N50 of 619 Kb and scaffold N50 of 1.12 Mb, indicating superior continuity. Approximately 99.83% of the genome was anchored to 46 chromosomes, revealing structural insights. Repetitive elements constituted 58.84% of the sequences. Functional annotations were assigned to 39,362 proteins, enhancing interpretative power. BUSCO analysis confirmed assembly completeness at 80.39%. This first high-quality E. gracilis genome offers insights for genetics and genomics studies, overcoming previous limitations. The impact extends to academic and industrial research, providing a foundational resource.

Unique features of conventional and nonconventional introns in Euglena gracilis

BACKGROUND

Nuclear introns in Euglenida have been understudied. This study aimed to investigate nuclear introns in Euglenida by identifying a large number of introns in Euglena gracilis (E. gracilis), including cis-spliced conventional and nonconventional introns, as well as trans-spliced outrons. We also examined the sequence characteristics of these introns.

RESULTS

A total of 28,337 introns and 11,921 outrons were identified. Conventional and nonconventional introns have distinct splice site features; the former harbour canonical GT/C-AG splice sites, whereas the latter are capable of forming structured motifs with their terminal sequences. We observed that short introns had a preference for canonical GT-AG introns. Notably, conventional introns and outrons in E. gracilis exhibited a distinct cytidine-rich polypyrimidine tract, in contrast to the thymidine-rich tracts observed in other organisms. Furthermore, the SL-RNAs in E. gracilis, as well as in other trans-splicing species, can form a recently discovered motif called the extended U6/5' ss duplex with the respective U6s. We also describe a novel type of alternative splicing pattern in E. gracilis. The tandem repeat sequences of introns in this protist were determined, and their contents were comparable to those in humans.

CONCLUSIONS

Our findings highlight the unique features of E. gracilis introns and provide insights into the splicing mechanism of these introns, as well as the genomics and evolution of Euglenida.

Let's talk about sex: Why reproductive systems matter for understanding algae

Sex is a crucial process that has molecular, genetic, cellular, organismal, and population-level consequences for eukaryotic evolution. Eukaryotic life cycles are composed of alternating haploid and diploid phases but are constrained by the need to accommodate the phenotypes of these different phases. Critical gaps in our understanding of evolutionary drivers of the diversity in algae life cycles include how selection acts to stabilize and change features of the life cycle. Moreover, most eukaryotes are partially clonal, engaging in both sexual and asexual reproduction. Yet, our understanding of the variation in their reproductive systems is largely based on sexual reproduction in animals or angiosperms. The relative balance of sexual versus asexual reproduction not only controls but also is in turn controlled by standing genetic variability, thereby shaping evolutionary trajectories. Thus, we must quantitatively assess the consequences of the variation in life cycles on reproductive systems. Algae are a polyphyletic group spread across many of the major eukaryotic lineages, providing powerful models by which to resolve this knowledge gap. There is, however, an alarming lack of data about the population genetics of most algae and, therefore, the relative frequency of sexual versus asexual processes. For many algae, the occurrence of sexual reproduction is unknown, observations have been lost in overlooked papers, or data on population genetics do not yet exist. This greatly restricts our ability to forecast the consequences of climate change on algal populations inhabiting terrestrial, aquatic, and marine ecosystems. This perspective summarizes our extant knowledge and provides some future directions to pursue broadly across micro- and macroalgal species.

Circular extrachromosomal DNA in Euglena gracilis under normal and stress conditions

Extrachromosomal circular DNA (eccDNA) enhances genomic plasticity, augmenting its coding and regulatory potential. Advances in high-throughput sequencing have enabled the investigation of these structural variants. Although eccDNAs have been investigated in numerous taxa, they remained understudied in euglenids. Therefore, we examined eccDNAs predicted from Illumina sequencing data of Euglena gracilis Z SAG 1224-5/25, grown under optimal photoperiod and exposed to UV irradiation. We identified approximately 1000 unique eccDNA candidates, about 20% of which were shared across conditions. We also observed a significant enrichment of mitochondrially encoded eccDNA in the UV-irradiated sample. Furthermore, we found that the heterogeneity of eccDNA was reduced in UV-exposed samples compared to cells that were grown in optimal conditions. Hence, eccDNA appears to play a role in the response to oxidative stress in Euglena, as it does in other studied organisms. In addition to contributing to the understanding of Euglena genomes, our results contribute to the validation of bioinformatics pipelines on a large, non-model genome.

Metabolic responses of Euglena gracilis under photoheterotrophic and heterotrophic conditions

The protist Euglena gracilis has various trophic modes including heterotrophy and photoheterotrophy. To investigate how cultivation mode influences metabolic regulation, the chemical composition of cellular metabolites of Euglena gracilis grown under heterotrophic and photoheterotrophic conditions was monitored from the early exponential phase to the mid-stationary phase using two different techniques, i.e, nuclear magnetic resonance (NMR) spectroscopy and high-resolution mass spectrometry (HRMS). The combined metabolomics approach allowed an in-depth understanding of the mechanism of photoheterotrophic and heterotrophic growth for biomolecule production. Heterotrophic conditions promoted the production of polar amino and oxygenated compounds such as proteins and polyphenol compounds, especially at the end of the exponential phase while photoheterotrophic cells enhanced the production of organoheterocyclic compounds, carbohydrates, and alkaloids.

The impact of elevated sulfur and nitrogen levels on cadmium tolerance in Euglena species

Heavy metal (HM) pollution threatens human and ecosystem health. Current methods for remediating water contaminated with HMs are expensive and have limited effect. Therefore, bioremediation is being investigated as an environmentally and economically viable alternative. Freshwater protists Euglena gracilis and Euglena mutabilis were investigated for their tolerance to cadmium (Cd). A greater increase in cell numbers under Cd stress was noted for E. mutabilis but only E. gracilis showed an increase in Cd tolerance following pre-treatment with elevated concentrations of S or N. To gain insight regarding the nature of the increased tolerance RNA-sequencing was carried out on E. gracilis. This revealed transcript level changes among pretreated cells, and additional differences among cells exposed to CdCl2. Gene ontology (GO) enrichment analysis reflected changes in S and N metabolism, transmembrane transport, stress response, and physiological processes related to metal binding. Identifying these changes enhances our understanding of how these organisms adapt to HM polluted environments and allows us to target development of future pre-treatments to enhance the use of E. gracilis in bioremediation relating to heavy metals.

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.

New plastids, old proteins: repeated endosymbiotic acquisitions in kareniacean dinoflagellates

Dinoflagellates are a diverse group of ecologically significant micro-eukaryotes that can serve as a model system for plastid symbiogenesis due to their susceptibility to plastid loss and replacement via serial endosymbiosis. Kareniaceae harbor fucoxanthin-pigmented plastids instead of the ancestral peridinin-pigmented ones and support them with a diverse range of nucleus-encoded plastid-targeted proteins originating from the haptophyte endosymbiont, dinoflagellate host, and/or lateral gene transfers (LGT). Here, we present predicted plastid proteomes from seven distantly related kareniaceans in three genera (Karenia, Karlodinium, and Takayama) and analyze their evolutionary patterns using automated tree building and sorting. We project a relatively limited ( ~ 10%) haptophyte signal pointing towards a shared origin in the family Chrysochromulinaceae. Our data establish significant variations in the functional distributions of these signals, emphasizing the importance of micro-evolutionary processes in shaping the chimeric proteomes. Analysis of plastid genome sequences recontextualizes these results by a striking finding the extant kareniacean plastids are in fact not all of the same origin, as two of the studied species (Karlodinium armiger, Takayama helix) possess plastids from different haptophyte orders than the rest.

Euglena's atypical respiratory chain adapts to the discoidal cristae and flexible metabolism

Euglena gracilis, a model organism of the eukaryotic supergroup Discoba harbouring also clinically important parasitic species, possesses diverse metabolic strategies and an atypical electron transport chain. While structures of the electron transport chain complexes and supercomplexes of most other eukaryotic clades have been reported, no similar structure is currently available for Discoba, limiting the understandings of its core metabolism and leaving a gap in the evolutionary tree of eukaryotic bioenergetics. Here, we report high-resolution cryo-EM structures of Euglena's respirasome I + III2 + IV and supercomplex III2 + IV2. A previously unreported fatty acid synthesis domain locates on the tip of complex I's peripheral arm, providing a clear picture of its atypical subunit composition identified previously. Individual complexes are re-arranged in the respirasome to adapt to the non-uniform membrane curvature of the discoidal cristae. Furthermore, Euglena's conformationally rigid complex I is deactivated by restricting ubiquinone's access to its substrate tunnel. Our findings provide structural insights for therapeutic developments against euglenozoan parasite infections.

High-efficiency genome editing by Cas12a ribonucleoprotein complex in Euglena gracilis

Transgene-free genome editing based on clustered regularly interspaced short palindromic repeats (CRISPR) technology is key to achieving genetic engineering in microalgae for basic research and industrial applications. Euglena gracilis, a unicellular phytoflagellate microalga, is a promising biomaterial for foods, feeds, cosmetics and biofuels. However, methods for the genetic manipulation of E. gracilis are still limited. Here, we developed a high-efficiency, transgene-free genome editing method for E. gracilis using Lachnospiraceae bacterium CRISPR-associated protein 12a (LbCas12a) ribonucleoprotein (RNP) complex, which complements the previously established Cas9 RNP-based method. Through the direct delivery of LbCas12a-containing RNPs, our method reached mutagenesis rates of approximately 77.2-94.5% at two different E. gracilis target genes, Glucan synthase-like 2 (EgGSL2) and a phytoene synthase gene (EgcrtB). Moreover, in addition to targeted mutagenesis, we demonstrated efficient knock-in and base editing at the target site using LbCas12a-based RNPs with a single-stranded DNA donor template in E. gracilis. This study extends the genetic engineering capabilities of Euglena to accelerate its basic use for research and engineering for bioproduction.

Deviating from the norm: Nuclear organisation in trypanosomes

At first glance the nucleus is a highly conserved organelle. Overall nuclear morphology, the octagonal nuclear pore complex, the presence of peripheral heterochromatin and the nuclear envelope appear near constant features right down to the ultrastructural level. New work is revealing significant compositional divergence within these nuclear structures and their associated functions, likely reflecting adaptations and distinct mechanisms between eukaryotic lineages and especially the trypanosomatids. While many examples of mechanistic divergence currently lack obvious functional interpretations, these studies underscore the malleability of nuclear architecture. I will discuss some recent findings highlighting these facets within trypanosomes, together with the underlying evolutionary framework and make a call for the exploration of nuclear function in non-canonical experimental organisms.

Lessons from the deep: mechanisms behind diversification of eukaryotic protein complexes

Genetic variation is the major mechanism behind adaptation and evolutionary change. As most proteins operate through interactions with other proteins, changes in protein complex composition and subunit sequence provide potentially new functions. Comparative genomics can reveal expansions, losses and sequence divergence within protein-coding genes, but in silico analysis cannot detect subunit substitutions or replacements of entire protein complexes. Insights into these fundamental evolutionary processes require broad and extensive comparative analyses, from both in silico and experimental evidence. Here, we combine data from both approaches and consider the gamut of possible protein complex compositional changes that arise during evolution, citing examples of complete conservation to partial and total replacement by functional analogues. We focus in part on complexes in trypanosomes as they represent one of the better studied non-animal/non-fungal lineages, but extend insights across the eukaryotes by extensive comparative genomic analysis. We argue that gene loss plays an important role in diversification of protein complexes and hence enhancement of eukaryotic diversity.

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.

Comparative analysis of mitochondrion-related organelles in anaerobic amoebozoans

Archamoebae comprises free-living or endobiotic amoebiform protists that inhabit anaerobic or microaerophilic environments and possess mitochondrion-related organelles (MROs) adapted to function anaerobically. We compared in silico reconstructed MRO proteomes of eight species (six genera) and found that the common ancestor of Archamoebae possessed very few typical components of the protein translocation machinery, electron transport chain and tricarboxylic acid cycle. On the other hand, it contained a sulphate activation pathway and bacterial iron-sulphur (Fe-S) assembly system of MIS-type. The metabolic capacity of the MROs, however, varies markedly within this clade. The glycine cleavage system is widely conserved among Archamoebae, except in Entamoeba, probably owing to its role in catabolic function or one-carbon metabolism. MRO-based pyruvate metabolism was dispensed within subgroups Entamoebidae and Rhizomastixidae, whereas sulphate activation could have been lost in isolated cases of Rhizomastix libera, Mastigamoeba abducta and Endolimax sp. The MIS (Fe-S) assembly system was duplicated in the common ancestor of Mastigamoebidae and Pelomyxidae, and one of the copies took over Fe-S assembly in their MRO. In Entamoebidae and Rhizomastixidae, we hypothesize that Fe-S cluster assembly in both compartments may be facilitated by dual localization of the single system. We could not find evidence for changes in metabolic functions of the MRO in response to changes in habitat; it appears that such environmental drivers do not strongly affect MRO reduction in this group of eukaryotes.

Optimization of Heterotrophic Culture Conditions for the Microalgae Euglena gracilis to Produce Proteins

Euglena gracilis is one of the few permitted edible microalgae. Considering consumer acceptance, E. gracilis grown heterotrophically with yellow appearances have wider food industrial applications such as producing meat analogs than green cells. However, there is much room to improve the protein content of heterotrophic culture cells. In this study, the effects of nitrogen sources, temperature, initial pH, and C/N ratios on the protein production of E. gracilis were evaluated under heterotrophic cultivation. These results indicated that ammonium sulfate was the optimal nitrogen source for protein production. The protein content of E. gracilis cultured by ammonium sulfate increased by 113% and 44.7% compared with that cultured by yeast extract and monosodium glutamate, respectively. The manipulation of the low C/N ratio further improved E. gracilis protein content to 66.10% (w/w), which was 1.6-fold of that in the C/N = 25 group. Additionally, amino acid analysis revealed that the nitrogen-to-protein conversion factor (NTP) could be affected by nitrogen sources. A superior essential amino acid index (EAAI) of 1.62 and a balanced amino acid profile further confirmed the high nutritional value of E. gracilis protein fed by ammonium sulfate. This study highlighted the vast potency of heterotrophic cultured E. gracilis as an alternative dietary protein source.

A unique mRNA decapping complex in trypanosomes

Removal of the mRNA 5' cap primes transcripts for degradation and is central for regulating gene expression in eukaryotes. The canonical decapping enzyme Dcp2 is stringently controlled by assembly into a dynamic multi-protein complex together with the 5'-3'exoribonuclease Xrn1. Kinetoplastida lack Dcp2 orthologues but instead rely on the ApaH-like phosphatase ALPH1 for decapping. ALPH1 is composed of a catalytic domain flanked by C- and N-terminal extensions. We show that T. brucei ALPH1 is dimeric in vitro and functions within a complex composed of the trypanosome Xrn1 ortholog XRNA and four proteins unique to Kinetoplastida, including two RNA-binding proteins and a CMGC-family protein kinase. All ALPH1-associated proteins share a unique and dynamic localization to a structure at the posterior pole of the cell, anterior to the microtubule plus ends. XRNA affinity capture in T. cruzi recapitulates this interaction network. The ALPH1 N-terminus is not required for viability in culture, but essential for posterior pole localization. The C-terminus, in contrast, is required for localization to all RNA granule types, as well as for dimerization and interactions with XRNA and the CMGC kinase, suggesting possible regulatory mechanisms. Most significantly, the trypanosome decapping complex has a unique composition, differentiating the process from opisthokonts.

Euglena gracilis: Biochemical properties of a membrane bound ecto-phosphatase activity modulated by fluoroaluminate complexes and different trophic conditions

The ecto-phosphatases belong to a group of enzymes closely associated with the cell surface that has its catalytic site facing the extracellular environment, where different phosphorylated substrates can be hydrolyzed. In the present work, we biochemically characterized the ecto-phosphatase activity of the freshwater microalgae Euglena gracilis, a model microorganism, ubiquitously distributed and resistant to several environmental stressors. The ecto-phosphatase activity is acidic, stimulated by copper and presents the following apparent kinetic parameters: Km = 2.52 ± 0.12 mM p-NPP and Vmax = 3.62 ± 0.06 nmol p-NP × h-1 × 106 cells. We observed that zinc, orthovanadate, molybdate, fluoride, and inorganic phosphate inhibit the ecto-phosphatase activity with different magnitudes. Fluoroaluminate complexes are also inhibitors of this ecto-phosphatase activity. They can be formed in the enzyme reaction conditions and are likely to occur in a natural environment where E. gracilis can be found. The ecto-phosphatase activity is constant through the culture growth phases and is negatively modulated after continuous subculturing in the dark when a shift from phototrophic to the heterotrophic metabolism is likely. The analysis of those biochemical parameters may contribute to understanding the role of E. gracilis ecto-phosphatase activity in natural environments.

Subunit composition of mitochondrial dehydrogenase complexes in diplonemid flagellates

In eukaryotes, pyruvate, a key metabolite produced by glycolysis, is converted by a tripartite mitochondrial pyruvate dehydrogenase (PDH) complex to acetyl-coenzyme A, which is fed into the tricarboxylic acid cycle. Two additional enzyme complexes with analogous composition catalyze similar oxidative decarboxylation reactions albeit using different substrates, the branched-chain ketoacid dehydrogenase (BCKDH) complex and the 2-oxoglutarate dehydrogenase (OGDH) complex. Comparative transcriptome analyses of diplonemids, one of the most abundant and diverse groups of oceanic protists, indicate that the conventional E1, E2, and E3 subunits of the PDH complex are lacking. E1 was apparently replaced in the euglenozoan ancestor of diplonemids by an AceE protein of archaeal type, a substitution that we also document in dinoflagellates. Here we demonstrate that the mitochondrion of the model diplonemid Paradiplonema papillatum displays pyruvate and 2-oxoglutarate dehydrogenase activities. Protein mass spectrometry of mitochondria reveal that the AceE protein is as abundant as the E1 subunit of BCKDH. This corroborates the view that the AceE subunit is a functional component of the PDH complex. We hypothesize that by acquiring AceE, the diplonemid ancestor not only lost the eukaryotic-type E1, but also the E2 and E3 subunits of the PDH complex, which are present in other euglenozoans. We posit that the PDH activity in diplonemids seems to be carried out by a complex, in which the AceE protein partners with the E2 and E3 subunits from BCKDH and/or OGDH.

Deep Proteogenomics of a Photosynthetic Cyanobacterium

Cyanobacteria, the evolutionary ancestors of plant chloroplasts, contribute substantially to the Earth's biogeochemical cycles and are of great interest for a sustainable economy. Knowledge of protein expression is the key to understanding cyanobacterial metabolism; however, proteome studies in cyanobacteria are limited and cover only a fraction of the theoretical proteome. Here, we performed a comprehensive proteogenomic analysis of the model cyanobacterium Synechocystis sp. PCC 6803 to characterize the expressed (phospho)proteome, re-annotate known and discover novel open reading frames (ORFs). By mapping extensive shotgun mass spectrometry proteomics data onto a six-frame translation of the Synechocystis genome, we refined the genomic annotation of 64 ORFs, including eight completely novel ORFs. Our study presents the largest reported (phospho)proteome dataset for a unicellular cyanobacterium, covering the expression of about 80% of the theoretical proteome under various cultivation conditions, such as nitrogen or carbon limitation. We report 568 phosphorylated S/T/Y sites that are present on numerous regulatory proteins, including the transcriptional regulators cyAbrB1 and cyAbrB2. We also catalogue the proteins that have never been detected under laboratory conditions and found that a large portion of them is plasmid-encoded. This dataset will serve as a resource, providing dedicated information on growth condition-dependent protein expression and phosphorylation.

Functional differentiation of Sec13 paralogues in the euglenozoan protists

The β-propeller protein Sec13 plays roles in at least three distinct processes by virtue of being a component of the COPII endoplasmic reticulum export vesicle coat, the nuclear pore complex (NPC) and the Seh1-associated (SEA)/GATOR nutrient-sensing complex. This suggests that regulatory mechanisms coordinating these cellular activities may operate via Sec13. The NPC, COPII and SEA/GATOR are all ancient features of eukaryotic cells, and in the vast majority of eukaryotes, a single Sec13 gene is present. Here we report that the Euglenozoa, a lineage encompassing the diplonemid, kinetoplastid and euglenid protists, possess two Sec13 paralogues. Furthermore, based on protein interactions and localization studies we show that in diplonemids Sec13 functions are divided between the Sec13a and Sec13b paralogues. Specifically, Sec13a interacts with COPII and the NPC, while Sec13b interacts with Sec16 and components of the SEA/GATOR complex. We infer that euglenozoan Sec13a is responsible for NPC functions and canonical anterograde transport activities while Sec13b acts within nutrient and autophagy-related pathways, indicating a fundamentally distinct organization of coatomer complexes in euglenozoan flagellates.

Recent expansion of metabolic versatility in Diplonema papillatum, the model species of a highly speciose group of marine eukaryotes

BACKGROUND

Diplonemid flagellates are among the most abundant and species-rich of known marine microeukaryotes, colonizing all habitats, depths, and geographic regions of the world ocean. However, little is known about their genomes, biology, and ecological role.

RESULTS

We present the first nuclear genome sequence from a diplonemid, the type species Diplonema papillatum. The ~ 280-Mb genome assembly contains about 32,000 protein-coding genes, likely co-transcribed in groups of up to 100. Gene clusters are separated by long repetitive regions that include numerous transposable elements, which also reside within introns. Analysis of gene-family evolution reveals that the last common diplonemid ancestor underwent considerable metabolic expansion. D. papillatum-specific gains of carbohydrate-degradation capability were apparently acquired via horizontal gene transfer. The predicted breakdown of polysaccharides including pectin and xylan is at odds with reports of peptides being the predominant carbon source of this organism. Secretome analysis together with feeding experiments suggest that D. papillatum is predatory, able to degrade cell walls of live microeukaryotes, macroalgae, and water plants, not only for protoplast feeding but also for metabolizing cell-wall carbohydrates as an energy source. The analysis of environmental barcode samples shows that D. papillatum is confined to temperate coastal waters, presumably acting in bioremediation of eutrophication.

CONCLUSIONS

Nuclear genome information will allow systematic functional and cell-biology studies in D. papillatum. It will also serve as a reference for the highly diverse diplonemids and provide a point of comparison for studying gene complement evolution in the sister group of Kinetoplastida, including human-pathogenic taxa.

Versatile biotechnological applications of Euglena gracilis

Euglena gracilis is a freshwater protist possessing secondary chloroplasts of green algal origin. Various physical factors (e.g. UV) and chemical compounds (e.g. antibiotics) cause the bleaching of E. gracilis cells-the loss of plastid genes leading to the permanent inability to photosynthesize. Bleaching can be prevented by antimutagens (i.e. lignin, vitamin C and selenium). Besides screening the mutagenic and antimutagenic activity of chemicals, E. gracilis is also a suitable model for studying the biological effects of many organic pollutants. Due to its capability of heavy metal sequestration, it can be used for bioremediation. E. gracilis has been successfully transformed, offering the possibility of genetic modifications for synthesizing compounds of biotechnological interest. The novel design of the "next generation" transgenic expression cassettes with respect to the specificities of euglenid gene expression is proposed. Moreover, E. gracilis is a natural source of commercially relevant bioproducts such as (pro)vitamins, wax esters, polyunsaturated fatty acids and paramylon (β-1,3-glucan). One of the highest limitations of large-scale cultivation of E. gracilis is its disability to synthesize essential vitamins B₁ and B₁₂. This disadvantage can be overcome by co-cultivation of E. gracilis with other microorganisms, which can synthesize sufficient amounts of these vitamins. Such co-cultures can be used for the effective accumulation and harvesting of Euglena biomass by bioflocculation.

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.

Phytohormones enhance heavy metal responses in Euglena gracilis: Evidence from uptake of Ni, Pb and Cd and linkages to hormonomic and metabolomic dynamics

Over the last decade, significant effort has been made to understand phytohormonal functions (e.g., cytokinins (CKs) and abscisic acid (ABA)) in metal stress responses of higher plants and algae. Despite the potential for these phytohormones to improve industrial remediation by Euglena gracilis (Euglenophyceae), no such roles have been elucidated for this highly adaptive species and its response to heavy metals. This study demonstrates that toxic metals (nickel, lead, cadmium) modify hormonal activity profiles (i.e., CK forms and their concentrations) in E. gracilis. Furthermore, exogenous ABA or CK (tZ) enabled higher metal uptake efficiency (i.e., 9.35% in lead and 9.2% in cadmium uptake with CK) and alleviated metal toxicity through the regulation of endogenous CKs (i.e., total CK, isoprenoid CK) and gibberellin (GAs, GA₁ and GA₃) levels. These responses suggest that E. gracilis regulates multiple phytohormone signals during metal stress acclimation. A deeper approach, using untargeted metabolomic analyses, gave more detailed insight into phytohormone-controlled pathways and associated modified metabolites, which were frequently related to metal accumulation and the physiological acclimation to metal presence. Significant changes in the levels of cellular metabolites, especially those involved in acclimation to metal stress, were under the influence of phytohormones in algal cells. When grown under metal stress conditions, the presence of exogenous ABA or CKs, caused changes in cellular metabolites which included those from: lipid pathways, riboflavin metabolism, the biosynthesis of cofactors/vitamins, and carbohydrate metabolism. Also, bioactive secondary metabolites (e.g., terpenoids, alkaloids, flavonoids, carotenoids) were modified in algal cells treated with phytohormones. Thus, the study gives a detailed view on the regulatory functions of ABA and CKs in algal metal bioremediation strategies, which are attributed to enhanced metal uptake and in the fine-tuning of plant hormone levels during metal stress response. The results can guide efforts to develop efficient, low-cost and environmentally friendly methods for bioremediation.

Group-specific functional patterns of mitochondrion-related organelles shed light on their multiple transitions from mitochondria in ciliated protists

Adaptations of ciliates to hypoxic environments have arisen independently several times. Studies on mitochondrion-related organelle (MRO) metabolisms from distinct anaerobic ciliate groups provide evidence for understanding the transitions from mitochondria to MROs within eukaryotes. To deepen our knowledge about the evolutionary patterns of ciliate anaerobiosis, mass-culture and single-cell transcriptomes of two anaerobic species, Metopus laminarius (class Armophorea) and Plagiopyla cf. narasimhamurtii (class Plagiopylea), were sequenced and their MRO metabolic maps were compared. In addition, we carried out comparisons using publicly available predicted MRO proteomes from other ciliate classes (i.e., Armophorea, Litostomatea, Muranotrichea, Oligohymenophorea, Parablepharismea and Plagiopylea). We found that single-cell transcriptomes were similarly comparable to their mass-culture counterparts in predicting MRO metabolic pathways of ciliates. The patterns of the components of the MRO metabolic pathways might be divergent among anaerobic ciliates, even among closely related species. Notably, our findings indicate the existence of group-specific functional relics of electron transport chains (ETCs). Detailed group-specific ETC functional patterns are as follows: full oxidative phosphorylation in Oligohymenophorea and Muranotrichea; only electron-transfer machinery in Armophorea; either of these functional types in Parablepharismea; and ETC functional absence in Litostomatea and Plagiopylea. These findings suggest that adaptation of ciliates to anaerobic conditions is group-specific and has occurred multiple times. Our results also show the potential and the limitations of detecting ciliate MRO proteins using single-cell transcriptomes and improve the understanding of the multiple transitions from mitochondria to MROs within ciliates.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-022-00147-w.

Euglena International Network (EIN): Driving euglenoid biotechnology for the benefit of a challenged world

Euglenoids (Euglenida) are unicellular flagellates possessing exceptionally wide geographical and ecological distribution. Euglenoids combine a biotechnological potential with a unique position in the eukaryotic tree of life. In large part these microbes owe this success to diverse genetics including secondary endosymbiosis and likely additional sources of genes. Multiple euglenoid species have translational applications and show great promise in production of biofuels, nutraceuticals, bioremediation, cancer treatments and more exotically as robotics design simulators. An absence of reference genomes currently limits these applications, including development of efficient tools for identification of critical factors in regulation, growth or optimization of metabolic pathways. The Euglena International Network (EIN) seeks to provide a forum to overcome these challenges. EIN has agreed specific goals, mobilized scientists, established a clear roadmap (Grand Challenges), connected academic and industry stakeholders and is currently formulating policy and partnership principles to propel these efforts in a coordinated and efficient manner.

N-terminal modifications, the associated processing machinery, and their evolution in plastid-containing organisms

The N-terminus is a frequent site of protein modifications. Referring primarily to knowledge gained from land plants, here we review the modifications that change protein N-terminal residues and provide updated information about the associated machinery, including that in Archaeplastida. These N-terminal modifications include many proteolytic events as well as small group additions such as acylation or arginylation and oxidation. Compared with that of the mitochondrion, the plastid-dedicated N-terminal modification landscape is far more complex. In parallel, we extend this review to plastid-containing Chromalveolata including Stramenopiles, Apicomplexa, and Rhizaria. We report a well-conserved machinery, especially in the plastid. Consideration of the two most abundant proteins on Earth-Rubisco and actin-reveals the complexity of N-terminal modification processes. The progressive gene transfer from the plastid to the nuclear genome during evolution is exemplified by the N-terminus modification machinery, which appears to be one of the latest to have been transferred to the nuclear genome together with crucial major photosynthetic landmarks. This is evidenced by the greater number of plastid genes in Paulinellidae and red algae, the most recent and fossil recipients of primary endosymbiosis.

Sending the message: specialized RNA export mechanisms in trypanosomes

Export of RNA from the nucleus is essential for all eukaryotic cells and has emerged as a major step in the control of gene expression. mRNA molecules are required to complete a complex series of processing events and pass a quality control system to protect the cytoplasm from the translation of aberrant proteins. Many of these events are highly conserved across eukaryotes, reflecting their ancient origin, but significant deviation from a canonical pathway as described from animals and fungi has emerged in the trypanosomatids. With significant implications for the mechanisms that control gene expression and hence differentiation, responses to altered environments and fitness as a parasite, these deviations may also reveal additional, previously unsuspected, mRNA export pathways.

Microbe Profile: Euglena gracilis: photogenic, flexible and hardy

Euglena gracilis is a unicellular photosynthetic eukaryotic flagellate of the Discoba supergroup, which also encompasses Kinetoplastida and Diplonema. Plastids have green algal origin and are secondarily acquired. The nuclear genome is extremely large and many genes suggest multiple endosymbiotic/gene transfer events, i.e. derivation from prokaryotes of various lineages. E. gracilis is remarkably robust and can proliferate in environments contaminated with heavy metals and acids. Extraordinary metabolic plasticity and a mixotrophic lifestyle confers an ability to thrive in a broad range of environments, as well as facilitating production of many novel metabolites, making Euglena of considerable biotechnological importance.

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.

Euglenatides, Potent Antiproliferative Cyclic Peptides Isolated from the Freshwater Photosynthetic Microalga Euglena gracilis

By limiting the nitrogen source to glutamic acid, we isolated cyclic peptides from Euglena gracilis containing asparagine and non-proteinogenic amino acids. Structure elucidation was accomplished through spectroscopic methods, mass spectrometry and chemical degradation. The euglenatides potently inhibit pathogenic fungi and cancer cell lines e.g., euglenatide B exhibiting IC50 values of 4.3 μM in Aspergillus fumigatus and 0.29 μM in MCF-7 breast cancer cells. In an unprecedented convergence of non-ribosomal peptide synthetase and polyketide synthase assembly-line biosynthesis between unicellular species and the metazoan kingdom, euglenatides bear resemblance to nemamides from Caenorhabditis elegans and inhibited both producing organisms E. gracilis and C. elegans. By molecular network analysis, we detected over forty euglenatide-like metabolites in E. gracilis, E. sanguinea and E. mutabilis, suggesting an important biological role for these natural products.

A Strategy for Gene Knockdown in Dinoflagellates

Dinoflagellates are unicellular protists that display unusual nuclear features such as large genomes, condensed chromosomes and multiple gene copies organized as tandem gene arrays. Genetic regulation is believed to be controlled at the translational rather than transcriptional level. An important player in this process is initiation factor eIF4E which binds the 7-methylguanosine cap structure (m7G) at the 5′-end of mRNA. Transcriptome analysis of eleven dinoflagellate species has established that each species encodes between eight to fifteen eIF4E family members. Determining the role of eIF4E family members in gene expression requires a method of knocking down their expression. In other eukaryotes this can be accomplished using translational blocking morpholinos that bind to complementary strands of RNA, therefore inhibiting the mRNA processing. Previously, unmodified morpholinos lacked the ability to pass through cell membranes, however peptide-based reagents have been used to deliver substances into the cytosol of cells by an endocytosis-mediated process without damaging the cell membrane. We have successfully delivered fluorescently-tagged morpholinos to the cytosol of Amphidinium carterae by using a specific cell penetrating peptide with the goal to target an eIF4e-1a sequence to inhibit translation. Specific eIF4e knockdown success (up to 42%) has been characterized via microscopy and western blot analysis.

A Synthetic Biology Perspective on the Bioengineering Tools for an Industrial Microalga: Euglena gracilis

Euglena is a genus of single-celled eukaryotes that show both plant- and animal-like characteristics. Euglena gracilis, a model species, is of great academic interest for studying endosymbiosis and chloroplast development. As an industrial species, E. gracilis is also of primary biotechnological and economic importance as high value-added food, medicine, and cosmetic and high-quality feedstock for jet-fuel production because of its cells containing many high-value products, such as vitamins, amino acids, pigments, unsaturated fatty acids, and carbohydrate paramylon, as metabolites. For more than half a century, E. gracilis has been used as an industrial biotechnology platform for fundamental biology research, mainly exploring relevant physiological and biochemical method studies. Although many researchers focused on genetic engineering tools for E. gracilis in recent years, little progress has been achieved because of the lack of high-quality genome information and efficient techniques for genetic operation. This article reviewed the progress of the genetic transformation of E. gracilis, including methods for the delivery of exogenous materials and other advanced biotechnological tools for E. gracilis, such as CRISPR and RNA interference. We hope to provide a reference to improve the research in functional genomics and synthetic biology of Euglena.

Molecular Cross-Talk between Gravity- and Light-Sensing Mechanisms in Euglena gracilis

Euglena gracilis is a photosynthetic flagellate. To acquire a suitable position in its surrounding aquatic environment, it exploits light and gravity primarily as environmental cues. Several physiological studies have indicated a fine-tuned relationship between gravity sensing (gravitaxis) and light sensing in E. gracilis. However, the underlying molecular mechanism is largely unknown. The photoreceptor photoactivated adenylyl cyclase (PAC) has been studied for over a decade. Nevertheless, no direct/indirect interaction partner (upstream/downstream) has been reported for PAC. It has been shown that a specific protein, kinase A (PKA), showed to be involved in phototaxis and gravitaxis. The current study reports the localization of the specific PKA and its relationship with PAC.

Evolution and diversification of the nuclear envelope

Eukaryotic cells arose ~1.5 billion years ago, with the endomembrane system a central feature, facilitating evolution of intracellular compartments. Endomembranes include the nuclear envelope (NE) dividing the cytoplasm and nucleoplasm. The NE possesses universal features: a double lipid bilayer membrane, nuclear pore complexes (NPCs), and continuity with the endoplasmic reticulum, indicating common evolutionary origin. However, levels of specialization between lineages remains unclear, despite distinct mechanisms underpinning various nuclear activities. Several distinct modes of molecular evolution facilitate organellar diversification and  to understand which apply to the NE, we exploited proteomic datasets of purified nuclear envelopes from model systems for comparative analysis. We find enrichment of core nuclear functions amongst the widely conserved proteins to be less numerous than lineage-specific cohorts, but enriched in core nuclear functions. This, together with consideration of additional evidence, suggests that, despite a common origin, the NE has evolved as a highly diverse organelle with significant lineage-specific functionality.

Highly flexible metabolism of the marine euglenozoan protist Diplonema papillatum

BACKGROUND

The phylum Euglenozoa is a group of flagellated protists comprising the diplonemids, euglenids, symbiontids, and kinetoplastids. The diplonemids are highly abundant and speciose, and recent tools have rendered the best studied representative, Diplonema papillatum, genetically tractable. However, despite the high diversity of diplonemids, their lifestyles, ecological functions, and even primary energy source are mostly unknown.

RESULTS

We designed a metabolic map of D. papillatum cellular bioenergetic pathways based on the alterations of transcriptomic, proteomic, and metabolomic profiles obtained from cells grown under different conditions. Comparative analysis in the nutrient-rich and nutrient-poor media, as well as the absence and presence of oxygen, revealed its capacity for extensive metabolic reprogramming that occurs predominantly on the proteomic rather than the transcriptomic level. D. papillatum is equipped with fundamental metabolic routes such as glycolysis, gluconeogenesis, TCA cycle, pentose phosphate pathway, respiratory complexes, β-oxidation, and synthesis of fatty acids. Gluconeogenesis is uniquely dominant over glycolysis under all surveyed conditions, while the TCA cycle represents an eclectic combination of standard and unusual enzymes.

CONCLUSIONS

The identification of conventional anaerobic enzymes reflects the ability of this protist to survive in low-oxygen environments. Furthermore, its metabolism quickly reacts to restricted carbon availability, suggesting a high metabolic flexibility of diplonemids, which is further reflected in cell morphology and motility, correlating well with their extreme ecological valence.

The distinctive flagellar proteome of Euglena gracilis illuminates the complexities of protistan flagella adaptation

The eukaryotic flagellum/cilium is a prominent organelle with conserved structure and diverse functions. Euglena gracilis, a photosynthetic and highly adaptable protist, employs its flagella for both locomotion and environmental sensing. Using proteomics of isolated E. gracilis flagella we identify nearly 1700 protein groups, which challenges previous estimates of the protein complexity of motile eukaryotic flagella. We not only identified several unexpected similarities shared with mammalian flagella, including an entire glycolytic pathway and proteasome, but also document a vast array of flagella-based signal transduction components that coordinate gravitaxis and phototactic motility. By contrast, the pellicle was found to consist of > 900 protein groups, containing additional structural and signalling components. Our data identify significant adaptations within the E. gracilis flagellum, many of which are clearly linked to the highly flexible lifestyle.

Evolutionary Genomics of Sex-Related Chromosomes at the Base of the Green Lineage

Although sex is now accepted as a ubiquitous and ancestral feature of eukaryotes, direct observation of sex is still lacking in most unicellular eukaryotic lineages. Evidence of sex is frequently indirect and inferred from the identification of genes involved in meiosis from whole genome data and/or the detection of recombination signatures from genetic diversity in natural populations. In haploid unicellular eukaryotes, sex-related chromosomes are named mating-type (MTs) chromosomes and generally carry large genomic regions where recombination is suppressed. These regions have been characterized in Fungi and Chlorophyta and determine gamete compatibility and fusion. Two candidate MT+ and MT− alleles, spanning 450–650 kb, have recently been described in Ostreococcus tauri, a marine phytoplanktonic alga from the Mamiellophyceae class, an early diverging branch in the green lineage. Here, we investigate the architecture and evolution of these candidate MT+ and MT− alleles. We analyzed the phylogenetic profile and GC content of MT gene families in eight different genomes whose divergence has been previously estimated at up to 640 Myr, and found evidence that the divergence of the two MT alleles predates speciation in the Ostreococcus genus. Phylogenetic profiles of MT trans-specific polymorphisms in gametologs disclosed candidate MTs in two additional species, and possibly a third. These Mamiellales MT candidates are likely to be the oldest mating-type loci described to date, which makes them fascinating models to investigate the evolutionary mechanisms of haploid sex determination in eukaryotes.

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.

Agrobacterium tumefaciens-Mediated Nuclear Transformation of a Biotechnologically Important Microalga-Euglena gracilis

Euglena gracilis (E. gracilis) is an attractive organism due to its evolutionary history and substantial potential to produce biochemicals of commercial importance. This study describes the establishment of an optimized protocol for the genetic transformation of E. gracilis mediated by Agrobacterium (A. tumefaciens). E. gracilis was found to be highly sensitive to hygromycin and zeocin, thus offering a set of resistance marker genes for the selection of transformants. A. tumefaciens-mediated transformation (ATMT) yielded hygromycin-resistant cells. However, hygromycin-resistant cells hosting the gus gene (encoding β-glucuronidase (GUS)) were found to be GUS-negative, indicating that the gus gene had explicitly been silenced. To circumvent transgene silencing, GUS was expressed from the nuclear genome as transcriptional fusions with the hygromycin resistance gene (hptII) (encoding hygromycin phosphotransferase II) with the foot and mouth disease virus (FMDV)-derived 2A self-cleaving sequence placed between the coding sequences. ATMT of Euglena with the hptII-2A–gus gene yielded hygromycin-resistant, GUS-positive cells. The transformation was verified by PCR amplification of the T-DNA region genes, determination of GUS activity, and indirect immunofluorescence assays. Cocultivation factors optimization revealed that a higher number of transformants was obtained when A. tumefaciens LBA4404 (A₆₀₀ = 1.0) and E. gracilis (A₇₅₀ = 2.0) cultures were cocultured for 48 h at 19 °C in an organic medium (pH 6.5) containing 50 µM acetosyringone. Transformation efficiency of 8.26 ± 4.9% was achieved under the optimized cocultivation parameters. The molecular toolkits and method presented here can be used to bioengineer E. gracilis for producing high-value products and fundamental studies.

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.

Comparative Genomics of Peroxisome Biogenesis Proteins: Making Sense of the PEX Proteins

PEX genes encode proteins involved in peroxisome biogenesis and proliferation. Using a comparative genomics approach, we clarify the evolutionary relationships between the 37 known PEX proteins in a representative set of eukaryotes, including all common model organisms, pathogenic unicellular eukaryotes and human. A large number of previously unknown PEX orthologs were identified. We analyzed all PEX proteins, their conservation and domain architecture and defined the core set of PEX proteins that is required to make a peroxisome. The molecular processes in peroxisome biogenesis in different organisms were put into context, showing that peroxisomes are not static organelles in eukaryotic evolution. Organisms that lack peroxisomes still contain a few PEX proteins, which probably play a role in alternative processes. Finally, the relationships between PEX proteins of two large families, the Pex11 and Pex23 families, were analyzed, thereby contributing to the understanding of their complicated and sometimes incorrect nomenclature. We provide an exhaustive overview of this important eukaryotic organelle.