Complete sequence of Euglena gracilis chloroplast DNA

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.

Modes and mechanisms for the inheritance of mitochondria and plastids in pathogenic protists

Pathogenic protists are responsible for many diseases that significantly impact human and animal health across the globe. Almost all protists possess mitochondria or mitochondrion-related organelles, and many contain plastids. These endosymbiotic organelles are crucial to survival and provide well-validated and widely utilised drug targets in parasitic protists such as Plasmodium and Toxoplasma. However, mutations within the organellar genomes of mitochondria and plastids can lead to drug resistance. Such mutations ultimately challenge our ability to control and eradicate the diseases caused by these pathogenic protists. Therefore, it is important to understand how organellar genomes, and the resistance mutations encoded within them, are inherited during protist sexual reproduction and how this may impact the spread of drug resistance and future therapeutic approaches to target these organelles. In this review, we detail what is known about mitochondrial and plastid inheritance during sexual reproduction across different pathogenic protists, often turning to their better studied, nonpathogenic relatives for insight.

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.

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.

Hypothetical chloroplast reading frame 51 encodes a photosystem I assembly factor in cyanobacteria

Hypothetical chloroplast open reading frames (ycfs) are putative genes in the plastid genomes of photosynthetic eukaryotes. Many ycfs are also conserved in the genomes of cyanobacteria, the presumptive ancestors of present-day chloroplasts. The functions of many ycfs are still unknown. Here, we generated knock-out mutants for ycf51 (sll1702) in the cyanobacterium Synechocystis sp. PCC 6803. The mutants showed reduced photoautotrophic growth due to impaired electron transport between photosystem II (PSII) and PSI. This phenotype results from greatly reduced PSI content in the ycf51 mutant. The ycf51 disruption had little effect on the transcription of genes encoding photosynthetic complex components and the stabilization of the PSI complex. In vitro and in vivo analyses demonstrated that Ycf51 cooperates with PSI assembly factor Ycf3 to mediate PSI assembly. Furthermore, Ycf51 interacts with the PSI subunit PsaC. Together with its specific localization in the thylakoid membrane and the stromal exposure of its hydrophilic region, our data suggest that Ycf51 is involved in PSI complex assembly. Ycf51 is conserved in all sequenced cyanobacteria, including the earliest branching cyanobacteria of the Gloeobacter genus, and is also present in the plastid genomes of glaucophytes. However, Ycf51 has been lost from other photosynthetic eukaryotic lineages. Thus, Ycf51 is a PSI assembly factor that has been functionally replaced during the evolution of oxygenic photosynthetic eukaryotes.

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.

Comparison of Boraginales Plastomes: Insights into Codon Usage Bias, Adaptive Evolution, and Phylogenetic Relationships

The Boraginales (Boraginaceae a.l.) comprise more than 2450 species worldwide. However, little knowledge exists of the characteristics of the complete plastid genome. In this study, three new sequences representing the first pt genome of Heliotropiaceae and Cordiaceae were assembled and compared with other Boraginales species. The pt genome sizes of Cordia dichotoma, Heliotropium arborescens, and Tournefortia montana were 151,990 bp, 156,243 bp, and 155,891 bp, respectively. Multiple optimal codons were identified, which may provide meaningful information for enhancing the gene expression of Boraginales species. Furthermore, codon usage bias analyses revealed that natural selection and other factors may dominate codon usage patterns in the Boraginales species. The boundaries of the IR/LSC and IR/SSC regions were significantly different, and we also found a signal of obvious IR region expansion in the pt genome of Nonea vesicaria and Arnebia euchroma. Genes with high nucleic acid diversity (pi) values were also calculated, which may be used as potential DNA barcodes to investigate the phylogenetic relationships in Boraginales. psaI, rpl33, rpl36, and rps19 were found to be under positive selection, and these genes play an important role in our understanding of the adaptive evolution of the Boraginales species. Phylogenetic analyses implied that Boraginales can be divided into two groups. The existence of two subfamilies (Lithospermeae and Boragineae) in Boraginaceae is also strongly supported. Our study provides valuable information on pt genome evolution and phylogenetic relationships in the Boraginales species.

Challenging the Importance of Plastid Genome Structure Conservation: New Insights From Euglenophytes

Plastids, similar to mitochondria, are organelles of endosymbiotic origin, which retained their vestigial genomes (ptDNA). Their unique architecture, commonly referred to as the quadripartite (four-part) structure, is considered to be strictly conserved; however, the bulk of our knowledge on their variability and evolutionary transformations comes from studies of the primary plastids of green algae and land plants. To broaden our perspective, we obtained seven new ptDNA sequences from freshwater species of photosynthetic euglenids-a group that obtained secondary plastids, known to have dynamically evolving genome structure, via endosymbiosis with a green alga. Our analyses have demonstrated that the evolutionary history of euglenid plastid genome structure is exceptionally convoluted, with a patchy distribution of inverted ribosomal operon (rDNA) repeats, as well as several independent acquisitions of tandemly repeated rDNA copies. Moreover, we have shown that inverted repeats in euglenid ptDNA do not share their genome-stabilizing property documented in chlorophytes. We hypothesize that the degeneration of the quadripartite structure of euglenid plastid genomes is connected to the group II intron expansion. These findings challenge the current global paradigms of plastid genome architecture evolution and underscore the often-underestimated divergence between the functionality of shared traits in primary and complex plastid organelles.

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.

Comparative Plastid Genomics of Green-Colored Dinoflagellates Unveils Parallel Genome Compaction and RNA Editing

Dinoflagellates possess plastids that are diverse in both pigmentation and evolutionary background. One of the plastid types found in dinoflagellates is pigmented with chlorophylls a and b (Chl a + b) and originated from the endosymbionts belonging to a small group of green algae, Pedinophyceae. The Chl a + b-containing plastids have been found in three distantly related dinoflagellates Lepidodinium spp., strain MGD, and strain TGD, and were proposed to be derived from separate partnerships between a dinoflagellate (host) and a pedinophycean green alga (endosymbiont). Prior to this study, a plastid genome sequence was only available for L. chlorophorum, which was reported to bear the features that were not found in that of the pedinophycean green alga Pedinomonas minor, a putative close relative of the endosymbiont that gave rise to the current Chl a + b-containing plastid. In this study, we sequenced the plastid genomes of strains MGD and TGD to compare with those of L. chlorophorum as well as pedinophycean green algae. The mapping of the RNA-seq reads on the corresponding plastid genome identified RNA editing on plastid gene transcripts in the three dinoflagellates. Further, the comparative plastid genomics revealed that the plastid genomes of the three dinoflagellates achieved several features, which are not found in or much less obvious than the pedinophycean plastid genomes determined to date, in parallel.

Discrimination of Euglena gracilis strains Z and bacillaris by MALDI-TOF MS

AIMS

Euglena gracilis is used as model organism for various microbiological, molecular biological and biotechnological studies. Its most studied wild type strains are Z and bacillaris, but their discrimination by standard molecular methods is difficult. Therefore, we decided to test the suitability of MALDI-TOF MS (matrix assisted laser desorption/ionization - time of flight mass spectrometry) for identification of E. gracilis and for discrimination of these two strains possessing functional chloroplasts. MALDI-TOF MS profiling was also tested for two white (non-photosynthetic) stable E. gracilis mutant strains W_gm ZOflL and W₁₀ BSmL.

METHODS AND RESULTS

We have successfully obtained main spectrum profiles (MSPs) of E. gracilis strains Z, SAG 1224-5/25 and bacillaris, SAG 1224-5/15 using protein extraction procedure. Subsequent MALDI-TOF MS profiling of a number of tested samples and the comparison of the obtained protein profiles with our in-house database including MSPs of both strains has revealed that these two strains can be easily distinguished by MALDI-TOF MS based on score values over two in most cases. This method has also confirmed the ancestry of white mutant strains W_gm ZOflL and W₁₀ BSmL, originally derived from strains Z and bacillaris, respectively.

CONCLUSIONS

MALDI-TOF MS is suitable, accurate and rapid method for discrimination of E. gracilis strains.

SIGNIFICANCE AND IMPACT OF STUDY

These results can have broad practical implications for laboratories cultivating various strains of euglenids, and they can be applied for their discrimination by MALDI-TOF MS.

Proteomic Responses of Dark-Adapted Euglena gracilis and Bleached Mutant Against Light Stimuli

Euglena gracilis (E. gracilis) has secondary endosymbiotic chloroplasts derived from ancient green algae. Its chloroplasts are easily lost under numerous conditions to become permanently bleached mutants. Green cells adapted in the dark contain undeveloped proplastids and they will develop into mature chloroplasts after 3 days of light exposure. Thus, E. gracilis is an ideal model species for a chloroplast development study. Previous studies about chloroplast development in E. gracilis focused on morphology and physiology, whereas few studies have addressed the regulatory processes induced by light in the proteome. In this study, the whole-genome proteome of dark-adapted E. gracilis (WT) and permanently ofloxacin-bleached mutant (B2) was compared under the light exposure after 0, 12, and 72 h. The results showed that the photosynthesis-related proteins were up-regulated over time in both WT and B2. The B2 strain, with losing functional chloroplasts, seemed to possess a complete photosynthetic function system. Both WT and B2 exhibited significant light responses with similar alternation patterns, suggesting the sensitive responses to light in proteomic levels. The main metabolic activities for the utilization of carbon and energy in WT were up-regulated, while the proteins with calcium ion binding, cell cycle, and non-photosynthetic carbon fixation were down-regulated in B2. This study confirmed light-induced chloroplast development in WT from dark, and also for the first time investigates the light responses of a bleached mutant B2, providing more information about the unknown functions of residual plastids in Euglena bleached mutants.

Organellar Evolution: A Path from Benefit to Dependence

Eukaryotic organelles supposedly evolved from their bacterial ancestors because of their benefits to host cells. However, organelles are quite often retained, even when the beneficial metabolic pathway is lost, due to something other than the original beneficial function. The organellar function essential for cell survival is, in the end, the result of organellar evolution, particularly losses of redundant metabolic pathways present in both the host and endosymbiont, followed by a gradual distribution of metabolic functions between the organelle and host. Such biological division of metabolic labor leads to mutual dependence of the endosymbiont and host. Changing environmental conditions, such as the gradual shift of an organism from aerobic to anaerobic conditions or light to dark, can make the original benefit useless. Therefore, it can be challenging to deduce the original beneficial function, if there is any, underlying organellar acquisition. However, it is also possible that the organelle is retained because it simply resists being eliminated or digested untill it becomes indispensable.

The strange mitochondrial genomes of Metschnikowia yeasts

While sequencing and characterizing the mitochondrial genomes of 71 strains from the yeast genus Metschnikowia [1] (close cousin to the model species Candida albicans), we uncovered one of the most extreme examples of mitochondrial genome architectural diversity observed to date. These Metschnikowia mitochondrial DNAs (mtDNAs) capture nearly the entire known gene-size and intron-content range for cox1 and cob across all eukaryotic life and show remarkable differences in structure and noncoding content. This genomic variation can be seen both among species and between strains of the same species, raising the question: why are Metschnikowia mitogenomes so malleable?

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.

Occurrence, Evolution and Specificities of Iron-Sulfur Proteins and Maturation Factors in Chloroplasts from Algae

Iron-containing proteins, including iron-sulfur (Fe-S) proteins, are essential for numerous electron transfer and metabolic reactions. They are present in most subcellular compartments. In plastids, in addition to sustaining the linear and cyclic photosynthetic electron transfer chains, Fe-S proteins participate in carbon, nitrogen, and sulfur assimilation, tetrapyrrole and isoprenoid metabolism, and lipoic acid and thiamine synthesis. The synthesis of Fe-S clusters, their trafficking, and their insertion into chloroplastic proteins necessitate the so-called sulfur mobilization (SUF) protein machinery. In the first part, we describe the molecular mechanisms that allow Fe-S cluster synthesis and insertion into acceptor proteins by the SUF machinery and analyze the occurrence of the SUF components in microalgae, focusing in particular on the green alga Chlamydomonas reinhardtii. In the second part, we describe chloroplastic Fe-S protein-dependent pathways that are specific to Chlamydomonas or for which Chlamydomonas presents specificities compared to terrestrial plants, putting notable emphasis on the contribution of Fe-S proteins to chlorophyll synthesis in the dark and to the fermentative metabolism. The occurrence and evolutionary conservation of these enzymes and pathways have been analyzed in all supergroups of microalgae performing oxygenic photosynthesis.

Constructive Neutral Evolution 20 Years Later

Evolution has led to a great diversity that ranges from elegant simplicity to ornate complexity. Many complex features are often assumed to be more functional or adaptive than their simpler alternatives. However, in 1999, Arlin Stolzfus published a paper in the Journal of Molecular Evolution that outlined a framework in which complexity can arise through a series of non-adaptive steps. He called this framework Constructive Neutral Evolution (CNE). Despite its two-decade-old roots, many evolutionary biologists still appear to be unaware of this explanatory framework for the origins of complexity. In this perspective piece, we explain the theory of CNE and how it changes the order of events in narratives that describe the evolution of complexity. We also provide an extensive list of cellular features that may have become more complex through CNE. We end by discussing strategies to determine whether complexity arose through neutral or adaptive processes.

Molecular tools and applications of Euglena gracilis: From biorefineries to bioremediation

Euglena gracilis is a promising source of commercially important metabolites such as vitamins, wax esters, paramylon, and amino acids. However, the molecular tools available to create improved Euglena strains are limited compared to other microorganisms that are currently exploited in the biotechnology industry. The complex poly-endosymbiotic nature of the Euglena genome is a major bottleneck for obtaining a complete genome sequence and thus represents a notable shortcoming in gaining molecular information of this organism. Therefore, the studies and applications have been more focused on using the wild-type strain or its variants and optimizing the nutrient composition and cultivation conditions to enhance the production of biomass and valuable metabolites. In addition to producing metabolites, the E. gracilis biorefinery concept also provides means for the production of biofuels and biogas as well as residual biomass for the remediation of industrial and municipal wastewater. Using Euglena for bioremediation of environments contaminated with heavy metals is of special interest due to the strong ability of the organism to accumulate and sequester these compounds. The published draft genome and transcriptome will serve as a basis for further molecular studies of Euglena and provide a guide for the engineering of metabolic pathways of relevance for the already established as well as novel applications.

Genetic Engineering Strategies for Euglena gracilis and Its Industrial Contribution to Sustainable Development Goals: A Review

The sustainable development goals (SDGs) adopted at the 2015 United Nations Summit are globally applicable goals designed to help countries realize a sustainable future. To achieve these SDGs, it is necessary to utilize renewable biological resources. In recent years, bioeconomy has been an attractive concept for achieving the SDGs. Microalgae are one of the biological resources that show promise in realizing the "5F"s (food, fiber, feed, fertilizer, and fuel). Among the microalgae, Euglena gracilis has the potential for achieving the "5F"s strategy owing to its unique features, such as production of paramylon, that are lacking in other microalgae. E. gracilis has already been produced on an industrial scale for use as an ingredient in functional foods and cosmetics. In recent years, genetic engineering methods for breeding E. gracilis have been researched and developed to achieve higher yields. In this article, we summarize how microalgae contribute toward achieving the SDGs. We focus on the contribution of E. gracilis to the bioeconomy, including its advantages in industrial use as well as its unique characteristics. In addition, we review genetic engineering-related research trends centered on E. gracilis, including a complete nuclear genome determination project, genome editing technology using the CRISPR-Cas9 system, and the development of a screening method for selecting useful strains. In particular, genome editing in E. gracilis could be a breakthrough for molecular breeding of industrially useful strains because of its high efficiency.

The Functions of Chloroplast Glutamyl-tRNA in Translation and Tetrapyrrole Biosynthesis

The chloroplast glutamyl-tRNA (tRNA^Glu) is unique in that it has two entirely different functions. In addition to acting in translation, it serves as the substrate of glutamyl-tRNA reductase (GluTR), the enzyme catalyzing the committed step in the tetrapyrrole biosynthetic pathway. How the tRNA^Glu pool is distributed between the two pathways and whether tRNA^Glu allocation limits tetrapyrrole biosynthesis and/or protein biosynthesis remains poorly understood. We generated a series of transplastomic tobacco (Nicotiana tabacum) plants to alter tRNA^Glu expression levels and introduced a point mutation into the plastid trnE gene, which has been reported to uncouple protein biosynthesis from tetrapyrrole biosynthesis in chloroplasts of the protist Euglena gracilis We show that, rather than comparable uncoupling of the two pathways, the trnE mutation is lethal in tobacco because it inhibits tRNA processing, thus preventing translation of Glu codons. Ectopic expression of the mutated trnE gene uncovered an unexpected inhibition of glutamyl-tRNA reductase by immature tRNA^Glu We further demonstrate that whereas overexpression of tRNA^Glu does not affect tetrapyrrole biosynthesis, reduction of GluTR activity through inhibition by tRNA^Glu precursors causes tetrapyrrole synthesis to become limiting in early plant development when active photosystem biogenesis provokes a high demand for de novo chlorophyll biosynthesis. Taken together, our findings provide insight into the roles of tRNA^Glu at the intersection of protein biosynthesis and tetrapyrrole biosynthesis.

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 analysis of the complete chloroplast genome among Prunus mume, P. armeniaca, and P. salicina

Prunus mume Sieb. et Zucc., P. armeniaca L., and P. salicina L. are economically important fruit trees in temperate regions. These species are taxonomically perplexing because of shared interspecific morphological traits and variation, which are mainly attributed to hybridization. The chloroplast is cytoplasmically inherited and often used for evolutionary studies. We sequenced the complete chloroplast genomes of P. mume, P. armeniaca, and P. salicina using Illumina sequencing followed by de novo assembly. The three chloroplast genomes exhibit a typical quadripartite structure with conserved genome arrangement, structure, and moderate divergence. The lengths of the genomes are 157,815, 157,797, and 157,916 bp, respectively. The length of the large single-copy region (LSC) region is 86,113, 86,283, and 86,122 bp, and the length of the SSC region is 18,916, 18,734, and 19,028 bp; the IR region is 26,393, 26,390, and 26,383 bp, respectively. Each of the three chloroplast genomes encodes 133 genes, including 94 protein-coding, 31 tRNA, and eight rRNA genes. Differential gene analysis for the three species revealed that trnY-ATA is a unique gene in P. armeniaca; in contrast, the gene trnI-TAT is only present in P. mume and P. salicina, though the position of the gene in these chloroplast genomes differs. Further comparative analysis of the complete chloroplast genome sequences revealed that the ORF genes and the sequences of linked regions rps16 and atpA, atpH and atpI, trnc-GCA and psbD, ycf3 and atpB, and rpL32 and ndhD are significantly different and may be used as molecular markers in taxonomic studies. Phylogenetic evolution analysis of the three species suggests that P. mume has a closer genetic relationship to P. armeniaca than to P. salicina.

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.

Directed Repeats Co-occur with Few Short-Dispersed Repeats in Plastid Genome of a Spikemoss, Selaginella vardei (Selaginellaceae, Lycopodiopsida)

BACKGROUND

It is hypothesized that the highly conserved inverted repeats (IR) structure of land plant plastid genomes (plastomes) is beneficial for stabilizing plastome organization, whereas the mechanism of the occurrence and stability maintenance of the recently reported direct repeats (DR) structure is yet awaiting further exploration. Here we describe the DR structure of the Selaginella vardei (Selaginellaceae) plastome, to elucidate the mechanism of DR occurrence and stability maintenance.

RESULTS

The plastome of S. vardei is 121,254 bp in length and encodes 76 genes, of which 62 encode proteins, 10 encode tRNAs, and four encode rRNAs. Unexpectedly, the two identical rRNA gene regions (13,893 bp) are arranged in a direct orientation (DR), rather than inverted. Comparing to the IR organization in Isoetes flaccida (Isoetaceae, Lycopodiopsida) plastome, a ca. 50-kb trnN-trnF inversion that spans one DR copy was found in the plastome of S. vardei, which might cause the orientation change. In addition, we find extremely rare short dispersed repeats (SDRs) in the plastomes of S. vardei and its closely related species S. indica.

CONCLUSIONS

We suggest that the ca. 50-kb inversion resulted in the DR structure, and the reduction in SDRs plays a key role in maintaining the stability of plastomes with DR structure by avoiding potential secondary recombination. We further confirmed the presence of homologous recombination between DR regions, which are able to generate subgenomes and form diverse multimers. Our study deepens the understanding of Selaginella plastomes and provides new insights into the diverse plastome structures in land plants.

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.

Dynamic evolution of inverted repeats in Euglenophyta plastid genomes

Photosynthetic euglenids (Euglenophyta) are a monophyletic group of unicellular eukaryotes characterized by the presence of plastids, which arose as the result of the secondary endosymbiosis. Many Euglenophyta plastid (pt) genomes have been characterized recently, but they represented mainly one family - Euglenaceae. Here, we report a comparative analysis of plastid genomes from eight representatives of the family Phacaceae. Newly sequenced plastid genomes share a number of features including synteny and gene content, except for genes mat2 and mat5 encoding maturases. The observed diversity of intron number and presence/absence of maturases corroborated previously suggested correlation between the number of maturases in the pt genome and intron proliferation. Surprisingly, pt genomes of taxa belonging to Discoplastis and Lepocinclis encode two inverted repeat (IR) regions containing the rDNA operon, which are absent from the Euglenaceae. By mapping the presence/absence of IR region on the obtained phylogenomic tree, we reconstructed the most probable events in the evolution of IRs in the Euglenophyta. Our study highlights the dynamic nature of the Euglenophyta plastid genome, in particular with regards to the IR regions that underwent losses repeatedly.

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.

Unraveling chloroplast transcriptomes with ChloroSeq, an organelle RNA-Seq bioinformatics pipeline

Online sequence repositories are teeming with RNA sequencing (RNA-Seq) data from a wide range of eukaryotes. Although most of these data sets contain large numbers of organelle-derived reads, researchers tend to ignore these data, focusing instead on the nuclear-derived transcripts. Consequently, GenBank contains massive amounts of organelle RNA-Seq data that are just waiting to be downloaded and analyzed. Recently, a team of scientists designed an open-source bioinformatics program called ChloroSeq, which systemically analyzes an organelle transcriptome using RNA-Seq. The ChloroSeq pipeline uses RNA-Seq alignment data to deliver detailed analyses of organelle transcriptomes, which can be fed into statistical software for further analysis and for generating graphical representations of the data. In addition to providing data on expression levels via coverage statistics, ChloroSeq can examine splicing efficiency and RNA editing profiles. Ultimately, ChloroSeq provides a well-needed avenue for researchers of all stripes to start exploring organelle transcription and could be a key step toward a more thorough understanding of organelle gene expression.

Intrageneric Variability Between the Chloroplast Genomes of Trachelomonas grandis and Trachelomonas volvocina and Phylogenomic Analysis of Phototrophic Euglenoids

The latest studies of chloroplast genomes of phototrophic euglenoids yielded different results according to intrageneric variability such as cluster arrangement or diversity of introns. Although the genera Euglena and Monomorphina in those studies show high syntenic arrangements at the intrageneric level, the two investigated Eutreptiella species comprise low synteny. Furthermore Trachelomonas volvocina show low synteny to the chloroplast genomes of the sister genera Monomorphina aenigmatica, M. parapyrum, Cryptoglena skujae, Euglenaria anabaena, Strombomonas acuminata, all of which were highly syntenic. Consequently, this study aims at the analysis of the cpGenome of Trachelomonas grandis and a comparative examination of T. volvocina to investigate whether the cpGenomes are of such resemblance as could be expected for a genus within the Euglenaceae. Although these analyses resulted in almost identical gene content to other Euglenaceae, the chloroplast genome showed significant novelties: In the rRNA operon, we detected group II introns, not yet found in any other cpGenome of Euglenaceae and a substantially heterogeneous cluster arrangement in the genus Trachelomonas. The phylogenomic analysis with 84 genes of 19 phototrophic euglenoids and 18 cpGenome sequences from Chlorophyta and Streptophyta resulted in a well-supported cpGenome phylogeny, which is in accordance to former phylogenetic analyses.

Chloroplast genome expansion by intron multiplication in the basal psychrophilic euglenoid Eutreptiella pomquetensis

BACKGROUND

Over the last few years multiple studies have been published showing a great diversity in size of chloroplast genomes (cpGenomes), and in the arrangement of gene clusters, in the Euglenales. However, while these genomes provided important insights into the evolution of cpGenomes across the Euglenales and within their genera, only two genomes were analyzed in regard to genomic variability between and within Euglenales and Eutreptiales. To better understand the dynamics of chloroplast genome evolution in early evolving Eutreptiales, this study focused on the cpGenome of Eutreptiella pomquetensis, and the spread and peculiarities of introns.

METHODS

The Etl. pomquetensis cpGenome was sequenced, annotated and afterwards examined in structure, size, gene order and intron content. These features were compared with other euglenoid cpGenomes as well as those of prasinophyte green algae, including Pyramimonas parkeae.

RESULTS AND DISCUSSION

With about 130,561 bp the chloroplast genome of Etl. pomquetensis, a basal taxon in the phototrophic euglenoids, was considerably larger than the two other Eutreptiales cpGenomes sequenced so far. Although the detected quadripartite structure resembled most green algae and plant chloroplast genomes, the gene content of the single copy regions in Etl. pomquetensis was completely different from those observed in green algae and plants. The gene composition of Etl. pomquetensis was extensively changed and turned out to be almost identical to other Eutreptiales and Euglenales, and not to P. parkeae. Furthermore, the cpGenome of Etl. pomquetensis was unexpectedly permeated by a high number of introns, which led to a substantially larger genome. The 51 identified introns of Etl. pomquetensis showed two major unique features: (i) more than half of the introns displayed a high level of pairwise identities; (ii) no group III introns could be identified in the protein coding genes. These findings support the hypothesis that group III introns are degenerated group II introns and evolved later.

Mobile Group II Introns as Ancestral Eukaryotic Elements

The duality of group II introns, capable of carrying out both self-splicing and retromobility reactions, is hypothesized to have played a profound role in the evolution of eukaryotes. These introns likely provided the framework for the emergence of eukaryotic retroelements, spliceosomal introns and other key components of the spliceosome. Group II introns are found in all three domains of life and are therefore considered to be exceptionally successful mobile genetic elements. Initially identified in organellar genomes, group II introns are found in bacteria, chloroplasts, and mitochondria of plants and fungi, but not in nuclear genomes. Although there is no doubt that prokaryotic and organellar group II introns are evolutionary related, there are remarkable differences in survival strategies between them. Furthermore, an evolutionary relationship of group II introns to eukaryotic retroelements, including telomeres, and spliceosomes is unmistakable.

A rare case of plastid protein-coding gene duplication in the chloroplast genome of Euglena archaeoplastidiata (Euglenophyta)

Gene duplication is an important evolutionary process that allows duplicate functions to diverge, or, in some cases, allows for new functional gains. However, in contrast to the nuclear genome, gene duplications within the chloroplast are extremely rare. Here, we present the chloroplast genome of the photosynthetic protist Euglena archaeoplastidiata. Upon annotation, it was found that the chloroplast genome contained a novel tandem direct duplication that encoded a portion of RuBisCO large subunit (rbcL) followed by a complete copy of ribosomal protein L32 (rpl32), as well as the associated intergenic sequences. Analyses of the duplicated rpl32 were inconclusive regarding selective pressures, although it was found that substitutions in the duplicated region, all non-synonymous, likely had a neutral functional effect. The duplicated region did not exhibit patterns consistent with previously described mechanisms for tandem direct duplications, and demonstrated an unknown mechanism of duplication. In addition, a comparison of this chloroplast genome to other previously characterized chloroplast genomes from the same family revealed characteristics that indicated E. archaeoplastidiata was probably more closely related to taxa in the genera Monomorphina, Cryptoglena, and Euglenaria than it was to other Euglena taxa. Taken together, the chloroplast genome of E. archaeoplastidiata demonstrated multiple characteristics unique to the euglenoid world, and has justified the longstanding curiosity regarding this enigmatic taxon.

The first complete chloroplast genome sequences of Ulmus species by de novo sequencing: Genome comparative and taxonomic position analysis

Elm (Ulmus) has a long history of use as a high-quality heavy hardwood famous for its resistance to drought, cold, and salt. It grows in temperate, warm temperate, and subtropical regions. This is the first report of Ulmaceae chloroplast genomes by de novo sequencing. The Ulmus chloroplast genomes exhibited a typical quadripartite structure with two single-copy regions (long single copy [LSC] and short single copy [SSC] sections) separated by a pair of inverted repeats (IRs). The lengths of the chloroplast genomes from five Ulmus ranged from 158,953 to 159,453 bp, with the largest observed in Ulmus davidiana and the smallest in Ulmus laciniata. The genomes contained 137–145 protein-coding genes, of which Ulmus davidiana var. japonica and U. davidiana had the most and U. pumila had the fewest. The five Ulmus species exhibited different evolutionary routes, as some genes had been lost. In total, 18 genes contained introns, 13 of which (trnL-TAA⁺, trnL-TAA⁻, rpoC1^-, rpl2^-, ndhA^-, ycf1, rps12^-, rps12⁺, trnA-TGC⁺, trnA-TGC^-, trnV-TAC^-, trnI-GAT⁺, and trnI-GAT) were shared among all five species. The intron of ycf1 was the longest (5,675bp) while that of trnF-AAA was the smallest (53bp). All Ulmus species except U. davidiana exhibited the same degree of amplification in the IR region. To determine the phylogenetic positions of the Ulmus species, we performed phylogenetic analyses using common protein-coding genes in chloroplast sequences of 42 other species published in NCBI. The cluster results showed the closest plants to Ulmaceae were Moraceae and Cannabaceae, followed by Rosaceae. Ulmaceae and Moraceae both belonged to Urticales, and the chloroplast genome clustering results were consistent with their traditional taxonomy. The results strongly supported the position of Ulmaceae as a member of the order Urticales. In addition, we found a potential error in the traditional taxonomies of U. davidiana and U. davidiana var. japonica, which should be confirmed with a further analysis of their nuclear genomes. This study is the first report on Ulmus chloroplast genomes, which has significance for understanding photosynthesis, evolution, and chloroplast transgenic engineering.

Biochemistry and Physiology of Heavy Metal Resistance and Accumulation in Euglena

Free-living microorganisms may become suitable models for removal of heavy metals from polluted water bodies, sediments, and soils by using and enhancing their metal accumulating abilities. The available research data indicate that protists of the genus Euglena are a highly promising group of microorganisms to be used in bio-remediation of heavy metal-polluted aerobic and anaerobic acidic aquatic environments. This chapter analyzes the variety of biochemical mechanisms evolved in E. gracilis to resist, accumulate and remove heavy metals from the environment, being the most relevant those involving (1) adsorption to the external cell pellicle; (2) intracellular binding by glutathione and glutathione polymers, and their further compartmentalization as heavy metal-complexes into chloroplasts and mitochondria; (3) polyphosphate biosynthesis; and (4) secretion of organic acids. The available data at the transcriptional, kinetic and metabolic levels on these metabolic/cellular processes are herein reviewed and analyzed to provide mechanistic basis for developing genetically engineered Euglena cells that may have a greater removal and accumulating capacity for bioremediation and recycling of heavy metals.

Euglena gracilis Genome and Transcriptome: Organelles, Nuclear Genome Assembly Strategies and Initial Features

Euglena gracilis is a major component of the aquatic ecosystem and together with closely related species, is ubiquitous worldwide. Euglenoids are an important group of protists, possessing a secondarily acquired plastid and are relatives to the Kinetoplastidae, which themselves have global impact as disease agents. To understand the biology of E. gracilis, as well as to provide further insight into the evolution and origins of the Kinetoplastidae, we embarked on sequencing the nuclear genome; the plastid and mitochondrial genomes are already in the public domain. Earlier studies suggested an extensive nuclear DNA content, with likely a high degree of repetitive sequence, together with significant extrachromosomal elements. To produce a list of coding sequences we have combined transcriptome data from both published and new sources, as well as embarked on de novo sequencing using a combination of 454, Illumina paired end libraries and long PacBio reads. Preliminary analysis suggests a surprisingly large genome approaching 2 Gbp, with a highly fragmented architecture and extensive repeat composition. Over 80% of the RNAseq reads from E. gracilis maps to the assembled genome sequence, which is comparable with the well assembled genomes of T. brucei and T. cruzi. In order to achieve this level of assembly we employed multiple informatics pipelines, which are discussed here. Finally, as a preliminary view of the genome architecture, we discuss the tubulin and calmodulin genes, which highlight potential novel splicing mechanisms.

Proliferation of group II introns in the chloroplast genome of the green alga Oedocladium carolinianum (Chlorophyceae)

Background

The chloroplast genome sustained extensive changes in architecture during the evolution of the Chlorophyceae, a morphologically and ecologically diverse class of green algae belonging to the Chlorophyta; however, the forces driving these changes are poorly understood. The five orders recognized in the Chlorophyceae form two major clades: the CS clade consisting of the Chlamydomonadales and Sphaeropleales, and the OCC clade consisting of the Oedogoniales, Chaetophorales, and Chaetopeltidales. In the OCC clade, considerable variations in chloroplast DNA (cpDNA) structure, size, gene order, and intron content have been observed. The large inverted repeat (IR), an ancestral feature characteristic of most green plants, is present in Oedogonium cardiacum (Oedogoniales) but is lacking in the examined members of the Chaetophorales and Chaetopeltidales. Remarkably, the Oedogonium 35.5-kb IR houses genes that were putatively acquired through horizontal DNA transfer. To better understand the dynamics of chloroplast genome evolution in the Oedogoniales, we analyzed the cpDNA of a second representative of this order, Oedocladium carolinianum.

Methods

The Oedocladium cpDNA was sequenced and annotated. The evolutionary distances separating Oedocladium and Oedogonium cpDNAs and two other pairs of chlorophycean cpDNAs were estimated using a 61-gene data set. Phylogenetic analysis of an alignment of group IIA introns from members of the OCC clade was performed. Secondary structures and insertion sites of oedogonialean group IIA introns were analyzed.

Results

The 204,438-bp Oedocladium genome is 7.9 kb larger than the Oedogonium genome, but its repertoire of conserved genes is remarkably similar and gene order differs by only one reversal. Although the 23.7-kb IR is missing the putative foreign genes found in Oedogonium, it contains sequences coding for a putative phage or bacterial DNA primase and a hypothetical protein. Intergenic sequences are 1.5-fold longer and dispersed repeats are more abundant, but a smaller fraction of the Oedocladium genome is occupied by introns. Six additional group II introns are present, five of which lack ORFs and carry highly similar sequences to that of the ORF-less IIA intron shared with Oedogonium. Secondary structure analysis of the group IIA introns disclosed marked differences in the exon-binding sites; however, each intron showed perfect or nearly perfect base pairing interactions with its target site.

Discussion

Our results suggest that chloroplast genes rearrange more slowly in the Oedogoniales than in the Chaetophorales and raise questions as to what was the nature of the foreign coding sequences in the IR of the common ancestor of the Oedogoniales. They provide the first evidence for intragenomic proliferation of group IIA introns in the Viridiplantae, revealing that intron spread in the Oedocladium lineage likely occurred by retrohoming after sequence divergence of the exon-binding sites.

Convergent evolution of twintron-like configurations: One is never enough

Introns inserted within introns are commonly referred to as twintrons, however the original definition for twintron implied that splicing of the external member of the twintron could only proceed upon splicing of the internal member. This review examines the various types of twintron-like arrangements that have been reported and assigns them to either nested or twintron categories that are subdivided further into subtypes based on differences of their mode of splicing. Twintron-like arrangements evolved independently by fortuitous events among different types of introns but once formed they offer opportunities for the evolution of new regulatory strategies and/or novel genetic elements.

Mitochondrion-to-Chloroplast DNA Transfers and Intragenomic Proliferation of Chloroplast Group II Introns in Gloeotilopsis Green Algae (Ulotrichales, Ulvophyceae)

To probe organelle genome evolution in the Ulvales/Ulotrichales clade, the newly sequenced chloroplast and mitochondrial genomes of Gloeotilopsis planctonica and Gloeotilopsis sarcinoidea (Ulotrichales) were compared with those of Pseudendoclonium akinetum (Ulotrichales) and of the few other green algae previously sampled in the Ulvophyceae. At 105,236 bp, the G planctonica mitochondrial DNA (mtDNA) is the largest mitochondrial genome reported so far among chlorophytes, whereas the 221,431-bp G planctonica and 262,888-bp G sarcinoidea chloroplast DNAs (cpDNAs) are the largest chloroplast genomes analyzed among the Ulvophyceae. Gains of non-coding sequences largely account for the expansion of these genomes. Both Gloeotilopsis cpDNAs lack the inverted repeat (IR) typically found in green plants, indicating that two independent IR losses occurred in the Ulvales/Ulotrichales. Our comparison of the Pseudendoclonium and Gloeotilopsis cpDNAs offered clues regarding the mechanism of IR loss in the Ulotrichales, suggesting that internal sequences from the rDNA operon were differentially lost from the two original IR copies during this process. Our analyses also unveiled a number of genetic novelties. Short mtDNA fragments were discovered in two distinct regions of the G sarcinoidea cpDNA, providing the first evidence for intracellular inter-organelle gene migration in green algae. We identified for the first time in green algal organelles, group II introns with LAGLIDADG ORFs as well as group II introns inserted into untranslated gene regions. We discovered many group II introns occupying sites not previously documented for the chloroplast genome and demonstrated that a number of them arose by intragenomic proliferation, most likely through retrohoming.

Chloroplasts: state of research and practical applications of plastome sequencing

This review presents origins, structure and expression of chloroplast genomes. It also describes their sequencing, analysis and modification, focusing on potential practical uses and biggest challenges of chloroplast genome modification. During the evolution of eukaryotes, cyanobacteria are believed to have merged with host heterotrophic cell. Afterward, most of cyanobacterial genes from cyanobacteria were transferred to cell nucleus or lost in the process of endosymbiosis. As a result of these changes, a primary plastid was established. Nowadays, plastid genome (plastome) is almost always circular, has a size of 100-200 kbp (120-160 in land plants), and harbors 100-120 highly conserved unique genes. Plastids have their own gene expression system, which is similar to one of their cyanobacterial ancestors. Two different polymerases, plastid-derived PEP and nucleus-derived NEP, participate in transcription. Translation is similar to the one observed in cyanobacteria, but it also utilizes protein translation factors and positive regulatory mRNA elements absent from bacteria. Plastoms play an important role in genetic transformation. Transgenes are introduced into them either via gene gun (in undamaged tissues) or polyethylene glycol treatment (when protoplasts are targeted). Antibiotic resistance markers are the most common tool used for selection of transformed plants. In recent years, plastome transformation emerged as a promising alternative to nuclear transformation because of (1) high yield of target protein, (2) removing the risk of outcrossing with weeds, (3) lack of silencing mechanisms, and (4) ability to engineer the entire metabolic pathways rather than single gene traits. Currently, the main directions of such research regard: developing efficient enzyme, vaccine antigen, and biopharmaceutical protein production methods in plant cells and improving crops by increasing their resistance to a wide array of biotic and abiotic stresses. Because of that, the detailed knowledge of plastome structure and mechanism of functioning started to play a major role.

An intact plastid genome is essential for the survival of colorless Euglena longa but not Euglena gracilis

Euglena gracilis growth with antibacterial agents leads to bleaching, permanent plastid gene loss. Colorless Euglena (Astasia) longa resembles a bleached E. gracilis. To evaluate the role of bleaching in E. longa evolution, the effect of streptomycin, a plastid protein synthesis inhibitor, and ofloxacin, a plastid DNA gyrase inhibitor, on E. gracilis and E. longa growth and plastid DNA content were compared. E. gracilis growth was unaffected by streptomycin and ofloxacin. Quantitative PCR analyses revealed a time dependent loss of plastid genes in E. gracilis demonstrating that bleaching agents produce plastid gene deletions without affecting cell growth. Streptomycin and ofloxacin inhibited E. longa growth indicating that it requires plastid genes to survive. This suggests that evolutionary divergence of E. longa from E. gracilis was triggered by the loss of a cytoplasmic metabolic activity also occurring in the plastid. Plastid metabolism has become obligatory for E. longa cell growth. A process termed "intermittent bleaching", short term exposure to subsaturating concentrations of reversible bleaching agents followed by growth in the absence of a bleaching agent, is proposed as the molecular mechanism for E. longa plastid genome reduction. Various non-photosynthetic lineages could have independently arisen from their photosynthetic ancestors via a similar process.

The Chloroplast Genome of Euglena mutabilis-Cluster Arrangement, Intron Analysis, and Intrageneric Trends

A comparative analysis of the chloroplast genome of Euglena mutabilis underlined a high diversity in the evolution of plastids in euglenids. Gene clusters in more derived Euglenales increased in complexity with only a few, but remarkable changes in the genus Euglena. Euglena mutabilis differed from other Euglena species in a mirror-inverted arrangement of 12 from 15 identified clusters, making it very likely that the emergence at the base of the genus Euglena, which has been considered a long branch artifact, is truly a probable position. This was corroborated by many similarities in gene arrangement and orientation with Strombomonas and Monomorphina, rendering the genome organization of E. mutabilis in certain clusters as plesiomorphic feature. By RNA analysis exact exon-intron boundaries and the type of the 77 introns identified were mostly determined unambiguously. A detailed intron study of psbC pointed at two important issues: First, the number of introns varied even between species, and no trend from few to many introns could be observed. Second, mat1 was localized in Eutreptiales exclusively in intron 1, and mat2 was not identified. With the emergence of Euglenaceae in most species, a new intron containing mat2 inserted in front of the previous intron 1 and thereby became intron 2 with mat1.

RuBisCO in Non-Photosynthetic Alga Euglena longa: Divergent Features, Transcriptomic Analysis and Regulation of Complex Formation

Euglena longa, a close relative of the photosynthetic model alga Euglena gracilis, possesses an enigmatic non-photosynthetic plastid. Its genome has retained a gene for the large subunit of the enzyme RuBisCO (rbcL). Here we provide new data illuminating the putative role of RuBisCO in E. longa. We demonstrated that the E. longa RBCL protein sequence is extremely divergent compared to its homologs from the photosynthetic relatives, suggesting a possible functional shift upon the loss of photosynthesis. Similarly to E. gracilis, E. longa harbors a nuclear gene encoding the small subunit of RuBisCO (RBCS) as a precursor polyprotein comprising multiple RBCS repeats, but one of them is highly divergent. Both RBCL and the RBCS proteins are synthesized in E. longa, but their abundance is very low compared to E. gracilis. No RBCS monomers could be detected in E. longa, suggesting that processing of the precursor polyprotein is inefficient in this species. The abundance of RBCS is regulated post-transcriptionally. Indeed, blocking the cytoplasmic translation by cycloheximide has no immediate effect on the RBCS stability in photosynthetically grown E. gracilis, but in E. longa, the protein is rapidly degraded. Altogether, our results revealed signatures of evolutionary degradation (becoming defunct) of RuBisCO in E. longa and suggest that its biological role in this species may be rather unorthodox, if any.

Comparative Analysis of the Chloroplast Genomic Information of Cunninghamia lanceolata (Lamb.) Hook with Sibling Species from the Genera Cryptomeria D. Don, Taiwania Hayata, and Calocedrus Kurz

Chinese fir (Cunninghamia lanceolata (Lamb.) Hook) is an important coniferous tree species for timber production, which accounts for ~40% of log supply from plantations in southern China. Chloroplast genetic engineering is an exciting field to engineer several valuable tree traits. In this study, we revisited the published complete Chinese fir (NC_021437) and four other coniferous species chloroplast genome sequence in Taxodiaceae. Comparison of their chloroplast genomes revealed three unique inversions found in the downstream of the gene clusters and evolutionary divergence were found, although overall the chloroplast genomic structure of the Cupressaceae linage was conserved. We also investigated the phylogenetic position of Chinese fir among conifers by examining gene functions, selection forces, substitution rates, and the full chloroplast genome sequence. Consistent with previous molecular systematics analysis, the results provided a well-supported phylogeny framework for the Cupressaceae that strongly confirms the “basal” position of Cunninghamia lanceolata. The structure of the Cunninghamia lanceolata chloroplast genome showed a partial lack of one IR copy, rearrangements clearly occurred and slight evolutionary divergence appeared among the cp genome of C. lanceolata, Taiwania cryptomerioides, Taiwania flousiana, Calocedrus formosana and Cryptomeria japonica. The information from sequence divergence and length variation of genes could be further considered for bioengineering research.

The transcriptome of Euglena gracilis reveals unexpected metabolic capabilities for carbohydrate and natural product biochemistry

Euglena gracilis is a highly complex alga belonging to the green plant line that shows characteristics of both plants and animals, while in evolutionary terms it is most closely related to the protozoan parasites Trypanosoma and Leishmania. This well-studied organism has long been known as a rich source of vitamins A, C and E, as well as amino acids that are essential for the human diet. Here we present de novo transcriptome sequencing and preliminary analysis, providing a basis for the molecular and functional genomics studies that will be required to direct metabolic engineering efforts aimed at enhancing the quality and quantity of high value products from E. gracilis. The transcriptome contains over 30,000 protein-encoding genes, supporting metabolic pathways for lipids, amino acids, carbohydrates and vitamins, along with capabilities for polyketide and non-ribosomal peptide biosynthesis. The metabolic and environmental robustness of Euglena is supported by a substantial capacity for responding to biotic and abiotic stress: it has the capacity to deploy three separate pathways for vitamin C (ascorbate) production, as well as producing vitamin E (α-tocopherol) and, in addition to glutathione, the redox-active thiols nor-trypanothione and ovothiol.