Transcriptomic study reveals widespread spliced leader trans-splicing, short 5'-UTRs and potential complex carbon fixation mechanisms in the euglenoid Alga Eutreptiella sp

Monodopsis and Vischeria Genomes Shed New Light on the Biology of Eustigmatophyte Algae

Members of eustigmatophyte algae, especially Nannochloropsis and Microchloropsis, have been tapped for biofuel production owing to their exceptionally high lipid content. Although extensive genomic, transcriptomic, and synthetic biology toolkits have been made available for Nannochloropsis and Microchloropsis, very little is known about other eustigmatophytes. Here we present three near-chromosomal and gapless genome assemblies of Monodopsis strains C73 and C141 (60 Mb) and Vischeria strain C74 (106 Mb), which are the sister groups to Nannochloropsis and Microchloropsis in the order Eustigmatales. These genomes contain unusually high percentages of simple repeats, ranging from 12% to 21% of the total assembly size. Unlike Nannochloropsis and Microchloropsis, long interspersed nuclear element repeats are abundant in Monodopsis and Vischeria and might constitute the centromeric regions. We found that both mevalonate and nonmevalonate pathways for terpenoid biosynthesis are present in Monodopsis and Vischeria, which is different from Nannochloropsis and Microchloropsis that have only the latter. Our analysis further revealed extensive spliced leader trans-splicing in Monodopsis and Vischeria at 36-61% of genes. Altogether, the high-quality genomes of Monodopsis and Vischeria not only serve as the much-needed outgroups to advance Nannochloropsis and Microchloropsis research, but also shed new light on the biology and evolution of eustigmatophyte algae.

Abundant synthesis of long-chain polyunsaturated fatty acids in Eutreptiella sp. (Euglenozoa) revealed by chromatographic and transcriptomic analyses

Algal lipids are important molecules to store energy in algae and transfer energy in the marine food chain, and are potential materials for high value nutraceuticals (e.g., omega-3 fatty acids) or biofuel production. However, how lipid biosynthesis is regulated is not well understood in many species including Eutreptiella from the phylum of Euglenozoa. Here, we characterized the fatty acid (FA) profile of an Eutreptiella species isolated from Long Island Sound, USA, using gas chromatography-tandem mass spectrometry (GC/MS/MS) and investigated their biosynthesis pathways by transcriptome sequencing. We discovered 24 types of FAs including a relatively high proportion of long-chain unsaturated FAs. The abundances of C16, C18, and saturated FAs decreased when phosphate in the culture medium was depleted. Among the 24 FAs, docosahexaenoic acid (C22:6^{∆4,7,10,13,16,19} ) was most abundant, suggesting that Eutreptiella sp. preferentially invests in the synthesis of long-chain polyunsaturated fatty acids (LC-PFAs). Further transcriptomic analysis revealed that Eutreptiella sp. likely synthesizes LC-PFAs via ∆8 pathway and uses type I and II fatty acid synthases. Using RT-qPCR, we found that some of the lipid synthesis genes, such as β-ketoacyl-ACP reductase, fatty acid desaturase, acetyl-CoA carboxylase, acyl carrier protein, ∆8 desaturase, and Acyl-ACP thioesterase, were more actively expressed during light period, and two carbon fixation genes were up-regulated in the high-lipid illuminated cultures, suggesting a linkage between photosynthesis and lipid production. The lipid profile renders Eutreptiella sp. a nutritional prey and valuable source for nutraceuticals, and the biosynthesis pathway documented here will be useful for future research and applications.

PhycoCosm, a comparative algal genomics resource

Algae are a diverse, polyphyletic group of photosynthetic eukaryotes spanning nearly all eukaryotic lineages of life and collectively responsible for ∼50% of photosynthesis on Earth. Sequenced algal genomes, critical to understanding their complex biology, are growing in number and require efficient tools for analysis. PhycoCosm (https://phycocosm.jgi.doe.gov) is an algal multi-omics portal, developed by the US Department of Energy Joint Genome Institute to support analysis and distribution of algal genome sequences and other ‘omics’ data. PhycoCosm provides integration of genome sequence and annotation for >100 algal genomes with available multi-omics data and interactive web-based tools to enable algal research in bioenergy and the environment, encouraging community engagement and data exchange, and fostering new sequencing projects that will further these research goals.

A New Type of Circular RNA derived from Nonconventional Introns in Nuclear Genes of Euglenids

Nuclear protein-coding genes of euglenids (Discoba, Euglenozoa, Euglenida) contain conventional (spliceosomal) and nonconventional introns. The latter have been found only in euglenozoans. A unique feature of nonconventional introns is the ability to form a stable and slightly conserved RNA secondary structure bringing together intron ends and placing adjacent exons in proximity. To date, little is known about the mechanism of their excision (e.g. whether it involves the spliceosome or not). The tubA gene of Euglena gracilis harbors three conventional and three nonconventional introns. While the conventional introns are excised as lariats, nonconventional introns are present in the cell solely as circular RNAs with full-length ends. Based on this discovery as well as on previous observations indicating that nonconventional introns are observed frequently at unique positions of genes, we suggest that this new type of intronic circRNA might play a role in intron mobility.

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.

Differential Gene Expression Supports a Resource-Intensive, Defensive Role for Colony Production in the Bloom-Forming Haptophyte, Phaeocystis globosa

Phaeocystis globosa forms dense, monospecific blooms in temperate, northern waters. Blooms are usually dominated by the colonial morphotype—nonflagellated cells embedded in a secreted mucilaginous mass. Colonial Phaeocystis blooms significantly affect food‐web structure and function and negatively impact fisheries and aquaculture, but factors regulating colony formation remain enigmatic. Destructive P. globosa blooms have been reported in tropical and subtropical regions more recently and warm‐water blooms could become more common with continued climate change and coastal eutrophication. We therefore assessed genetic pathways associated with colony formation by investigating differential gene expression between colonial and solitary cells of a warm‐water P. globosa strain. Our results illustrate a transcriptional shift in colonial cells with most of the differentially expressed genes downregulated, supporting a reallocation of resources associated with forming and maintaining colonies. Dimethylsulfide and acrylate production and pathogen interaction pathways were upregulated in colonial cells, suggesting a defensive role for producing colonies. We identify several protein kinase signaling pathways that may influence the transition between morphotypes, providing targets for future research into factors affecting colony formation. This study provides novel insights into genetic mechanisms involved in Phaeocystis colony formation and provides new evidence supporting a defensive role for Phaeocystis colonies.

Landscape of the spliced leader trans-splicing mechanism in Schistosoma mansoni

Spliced leader dependent trans-splicing (SLTS) has been described as an important RNA regulatory process that occurs in different organisms, including the trematode Schistosoma mansoni. We identified more than seven thousand putative SLTS sites in the parasite, comprising genes with a wide spectrum of functional classes, which underlines the SLTS as a ubiquitous mechanism in the parasite. Also, SLTS gene expression levels span several orders of magnitude, showing that SLTS frequency is not determined by the expression level of the target gene, but by the presence of particular gene features facilitating or hindering the trans-splicing mechanism. Our in-depth investigation of SLTS events demonstrates widespread alternative trans-splicing (ATS) acceptor sites occurring in different regions along the entire gene body, highlighting another important role of SLTS generating alternative RNA isoforms in the parasite, besides the polycistron resolution. Particularly for introns where SLTS directly competes for the same acceptor substrate with cis-splicing, we identified for the first time additional and important features that might determine the type of splicing. Our study substantially extends the current knowledge of RNA processing by SLTS in S. mansoni, and provide basis for future studies on the trans-splicing mechanism in other eukaryotes.

Spliced leader-based analyses reveal the effects of polycyclic aromatic hydrocarbons on gene expression in the copepod Pseudodiaptomus poplesia

Polycyclic aromatic hydrocarbons (PAHs) are a group of toxic and carcinogenic pollutants that can adversely affect the development, growth and reproduction of marine organisms including copepods. However, knowledge on the molecular mechanisms regulating the response to PAH exposure in marine planktonic copepods is limited. In this study, we investigated the survival and gene expression of the calanoid copepod Pseudodiaptomus poplesia upon exposure to two PAHs, 1, 2-dimethylnaphthalene (1, 2-NAPH) and pyrene. Acute toxicity responses resulted in 96-h LC₅₀ of 788.98μgL^-1 and 54.68μgL^-1 for 1, 2-NAPH and pyrene, respectively. Using the recently discovered copepod spliced leader as a primer, we constructed full-length cDNA libraries from copepods exposed to sublethal concentrations and revealed 289 unique genes of diverse functions, including stress response genes and novel genes previously undocumented for this species. Eighty-three gene families were specifically expressed in PAH exposure libraries. We further analyzed the expression of seven target genes by reverse transcription-quantitative PCR in a time-course test with three sublethal concentrations. These target genes have primary roles in detoxification, oxidative defense, and signal transduction, and include different forms of glutathione S-transferase (GST), glutathione peroxidases (GPX), peroxiredoxin (PRDX), methylmalonate-semialdehyde dehydrogenase (MSDH) and ras-related C3 botulinum toxin substrate (RAC1). Expression stability of seven candidate reference genes were evaluated and the two most stable ones (RPL15 and RPS20 for 1, 2-NAPH exposure, RPL15 and EF1D for pyrene exposure) were used to normalize the expression levels of the target genes. Significant upregulation was detected in GST-T, GST-DE, GPX4, PRDX6 and RAC1 upon 1, 2-NAPH exposure, and GST-DE and MSDH upon pyrene exposure. These results indicated that the oxidative stress was induced and that signal transduction might be affected by PAH exposure in P. poplesia. However, gene upregulation was followed by a reduction in expression level towards 96h, indicating a threshold value of exposure time that leads to depressed gene expression. Prolonged exposure may cause dysfunction of detoxification and antioxidant machinery in P. poplesia. The transcriptional responses of GST-T, GPX2 and GPX4 upon pyrene exposure were minimal. Our results reveal the different sensitivity of P. poplesia to two PAHs at both the individual and transcriptional levels. As the first attempt, this study proved that copepod spliced leader is useful for obtaining full-length cDNA in P. poplesia exposed to PAHs and provided a valuable gene resource for this non-model species. This approach can be applied to other calanoid copepods exposed to various stressors, particularly under field conditions.

Protists and the Wild, Wild West of Gene Expression: New Frontiers, Lawlessness, and Misfits

The DNA double helix has been called one of life's most elegant structures, largely because of its universality, simplicity, and symmetry. The expression of information encoded within DNA, however, can be far from simple or symmetric and is sometimes surprisingly variable, convoluted, and wantonly inefficient. Although exceptions to the rules exist in certain model systems, the true extent to which life has stretched the limits of gene expression is made clear by nonmodel systems, particularly protists (microbial eukaryotes). The nuclear and organelle genomes of protists are subject to the most tangled forms of gene expression yet identified. The complicated and extravagant picture of the underlying genetics of eukaryotic microbial life changes how we think about the flow of genetic information and the evolutionary processes shaping it. Here, we discuss the origins, diversity, and growing interest in noncanonical protist gene expression and its relationship to genomic architecture.

Spliced leader RNA trans-splicing discovered in copepods

Copepods are one of the most abundant metazoans in the marine ecosystem, constituting a critical link in aquatic food webs and contributing significantly to the global carbon budget, yet molecular mechanisms of their gene expression are not well understood. Here we report the detection of spliced leader (SL) trans-splicing in calanoid copepods. We have examined nine species of wild-caught copepods from Jiaozhou Bay, China that represent the major families of the calanoids. All these species contained a common 46-nt SL (CopepodSL). We further determined the size of CopepodSL precursor RNA (slRNA; 108-158 nt) through genomic analysis and 3′-RACE technique, which was confirmed by RNA blot analysis. Structure modeling showed that the copepod slRNA folded into typical slRNA secondary structures. Using a CopepodSL-based primer set, we selectively enriched and sequenced copepod full-length cDNAs, which led to the characterization of copepod transcripts and the cataloging of the complete set of 79 eukaryotic cytoplasmic ribosomal proteins (cRPs) for a single copepod species. We uncovered the SL trans-splicing in copepod natural populations, and demonstrated that CopepodSL was a sensitive and specific tool for copepod transcriptomic studies at both the individual and population levels and that it would be useful for metatranscriptomic analysis of copepods.

Light-promoted rhodopsin expression and starvation survival in the marine dinoflagellate Oxyrrhis marina

The discovery of microbial rhodopsins in marine proteobacteria changed the dogma that photosynthesis is the only pathway to use the solar energy for biological utilization in the marine environment. Although homologs of these rhodopsins have been identified in dinoflagellates, the diversity of the encoding genes and their physiological roles remain unexplored. As an initial step toward addressing the gap, we conducted high-throughput transcriptome sequencing on Oxyrrhis marina to retrieve rhodopsin transcripts, rapid amplification of cDNA ends to isolate full-length cDNAs of dominant representatives, and quantitative reverse-transcription PCR to investigate their expression under varying conditions. Our phylogenetic analyses showed that O. marina contained both the proton-pumping type (PR) and sensory type (SR) rhodopsins, and the transcriptome data showed that the PR type dominated over the SR type. We compared rhodopsin gene expression for cultures kept under light: dark cycle and continuous darkness in a time course of 24 days without feeding. Although both types of rhodopsin were expressed under the two conditions, the expression levels of PR were much higher than SR, consistent with the transcriptomic data. Furthermore, relative to cultures kept in the dark, rhodopsin expression levels and cell survival rate were both higher in cultures grown in the light. This is the first report of light-dependent promotion of starvation survival and concomitant promotion of PR expression in a eukaryote. While direct evidence needs to come from functional test on rhodopsins in vitro or gene knockout/knockdown experiments, our results suggest that the proton-pumping rhodopsin might be responsible for the light-enhanced survival of O. marina, as previously demonstrated in bacteria.

De novo transcriptomes of a mixotrophic and a heterotrophic ciliate from marine plankton

Studying non-model organisms is crucial in the context of the current development of genomics and transcriptomics for both physiological experimentation and environmental characterization. We investigated the transcriptomes of two marine planktonic ciliates, the mixotrophic oligotrich Strombidium rassoulzadegani and the heterotrophic choreotrich Strombidinopsis sp., and their respective algal food using Illumina RNAseq. Our aim was to characterize the transcriptomes of these contrasting ciliates and to identify genes potentially involved in mixotrophy. We detected approximately 10,000 and 7,600 amino acid sequences for S. rassoulzadegani and Strombidinopsis sp., respectively. About half of these transcripts had significant BLASTP hits (E-value <10⁻⁶) against previously-characterized sequences, mostly from the model ciliate Oxytricha trifallax. Transcriptomes from both the mixotroph and the heterotroph species provided similar annotations for GO terms and KEGG pathways. Most of the identified genes were related to housekeeping activity and pathways such as the metabolism of carbohydrates, lipids, amino acids, nucleotides, and vitamins. Although S. rassoulzadegani can keep and use chloroplasts from its prey, we did not find genes clearly linked to chloroplast maintenance and functioning in the transcriptome of this ciliate. While chloroplasts are known sources of reactive oxygen species (ROS), we found the same complement of antioxidant pathways in both ciliates, except for one enzyme possibly linked to ascorbic acid recycling found exclusively in the mixotroph. Contrary to our expectations, we did not find qualitative differences in genes potentially related to mixotrophy. However, these transcriptomes will help to establish a basis for the evaluation of differential gene expression in oligotrichs and choreotrichs and experimental investigation of the costs and benefits of mixotrophy.

Signal recognition particle RNA in dinoflagellates and the Perkinsid Perkinsus marinus

In dinoflagellates and perkinsids, the molecular structure of the protein translocating machinery is unclear. Here, we identified several types of full-length signal recognition particle (SRP) RNA genes from Karenia brevis (dinoflagellate) and Perkinsus marinus (perkinsid). We also identified the four SRP S-domain proteins, but not the two Alu domain proteins, from P. marinus and several dinoflagellates. We mapped both ends of SRP RNA transcripts from K. brevis and P. marinus, and obtained the 3' end from four other dinoflagellates. The lengths of SRP RNA are predicted to be ∼260-300 nt in dinoflagellates and 280-285 nt in P. marinus. Although these SRP RNA sequences are substantially variable, the predicted structures are similar. The genomic organization of the SRP RNA gene differs among species. In K. brevis, this gene is located downstream of the spliced leader (SL) RNA, either as SL RNA-SRP RNA-tRNA gene tandem repeats, or within a SL RNA-SRP RNA-tRNA-U6-5S rRNA gene cluster. In other dinoflagellates, SRP RNA does not cluster with SL RNA or 5S rRNA genes. The majority of P. marinus SRP RNA genes array as tandem repeats without the above-mentioned small RNA genes. Our results capture a snapshot of a potentially complex evolutionary history of SRP RNA in alveolates.

Tandem repeats, high copy number and remarkable diel expression rhythm of form II RuBisCO in Prorocentrum donghaiense (Dinophyceae)

Gene structure and expression regulation of form II RuBisCO (rbcII) in dinoflagellates are still poorly understood. Here we isolated this gene (Pdrbc) and investigated its diel expression pattern in a harmful algal bloom forming dinoflagellate Prorocentrum donghaiense. We obtained cDNA sequences with triple tandem repeats of the coding unit (CU); the 5′ region has the sequence of a typical dinoflagellate plastid gene, encoding an N-terminus with two transmembrane regions separated by a plastid transit peptide. The CUs (1,455 bp except 1464 bp in last CU) are connected through a 63 bp spacer. Phylogenetic analysis showed that rbcII CUs within species formed monophyletic clusters, indicative of intraspecific gene duplication or purifying evolution. Using quantitative PCR (qPCR) we estimated 117±40 CUs of Pdrbc in the P. donghaiense genome. Although it is commonly believed that most dinoflagellate genes lack transcriptional regulation, our RT-qPCR analysis on synchronized cultures revealed remarkable diel rhythm of Pdrbc expression, showing significant correlations of transcript abundance with the timing of the dark-to-light transition and cell cycle G2M-phase. When the cultures were shifted to continuous light, Pdrbc expression remained significantly correlated with the G2M-phase. Under continuous darkness the cell cycle was arrested at the G1 phase, and the rhythm of Pdrbc transcription disappeared. Our results suggest that dinoflagellate rbcII 1) undergoes duplication or sequence purification within species, 2) is organized in tandem arrays in most species probably to facilitate efficient translation and import of the encoded enzyme, and 3) is regulated transcriptionally in a cell cycle-dependent fashion at least in some dinoflagellates.