Transcriptome de novo assembly sequencing and analysis of the toxic dinoflagellate Alexandrium catenella using the Illumina platform

Proteomic and toxicological analysis of the response of dinoflagellate Alexandrium catenella to changes in NaNO3 concentration

Dinoflagellates of the genus Alexandrium cause Harmful Algal Blooms (HABs) in coastal waters worldwide, damaging marine environments, aquaculture, and human health. They synthesize potent neurotoxic alkaloids known as PSTs (i.e., Paralytic Shellfish Toxins), the etiological agents of PSP (i.e., Paralytic Shellfish Poisoning). In recent decades, the eutrophication of coastal waters with inorganic nitrogen (e.g., nitrate, nitrite, and ammonia) has increased the frequency and scale of HABs. PSTs concentrations within Alexandrium cells can increase by up to 76% after a nitrogen enrichment event; however, the mechanisms that underlie their biosynthesis in dinoflagellates remains unclear. This study combines mass spectrometry, bioinformatics, and toxicology and investigates the expression profiles of PSTs in Alexandrium catenella grown in 0.4, 0.9 and 1.3 mM NaNO₃. Pathway analysis of protein expression revealed that tRNA amino acylation, glycolysis, TCA cycle and pigment biosynthesis were upregulated in 0.4 mM and downregulated in 1.3 mM NaNO₃ compared to those grown in 0.9 mM NaNO₃. Conversely, ATP synthesis, photosynthesis and arginine biosynthesis were downregulated in 0.4 mM and upregulated in 1.3 mM NaNO₃. Additionally, the expression of proteins involved in PST biosynthesis (sxtA, sxtG, sxtV, sxtW and sxtZ) and overall PST production like STX, NEO, C1, C2, GTX1-6 and dcGTX2 was higher at lower nitrate concentrations. Therefore, increased nitrogen concentrations increase protein synthesis, photosynthesis, and energy metabolism and decrease enzyme expression in PST biosynthesis and production. This research provides new clues about how the changes in the nitrate concentration can modulate different metabolic pathways and the expression of PST biosynthesis in toxigenic dinoflagellates.

Establishment and preliminary study of electrophysiological techniques in a typical red tide species

Electrophysiology studies the electrical properties of cells and tissues including bioelectrical signals and membrane ion channel activities. As an important means to reveal ion channel related physiological functions and the underlying mechanisms, electrophysiological techniques have been widely used in studies of animals, higher plants and algae that are closely related to higher plants. However, few electrophysiological studies have been carried out in red tide organisms, especially in dinoflagellates, which is mainly due to the complex surface structure of dinoflagellate amphiesma. In this study, the surface amphiesma of Alexandrium pacificum, a typical red tide species, was removed by centrifugation, low-temperature treatment and enzymatic treatment. In all three treatments, low-temperature treatment with 4 °C for 2 h had high ecdysis rate and high fixation rate, and the treated cells were easy to puncture, so low-temperature treatment was used as a preprocessing treatment for subsequent current recording. Acquired protoplasts of A. pacificum were identified by calcofluor fluorescence and immobilized by poly-lysine. A modified "puncture" single-electrode voltage-clamp recording was first applied to dinoflagellates, and voltage-gated currents, which had the characteristics of outward K⁺ current and inward Cl^- current, were recorded and confirmed by ion replacement, indicating the voltage-gated currents were mixed. This method can be used as a technical basis for the electrophysiological study of dinoflagellates and provides a new perspective for the study of stress tolerance, red tide succession, and the regulation of physiological function of dinoflagellates.

Isolation from a fish kill and transcriptomic characterization of Gyrodinium jinhaense off Long Island Sound

First found in Korean coastal water, the dinoflagellate Gyrodinium jinhaense is a recently established species with unclear global distribution and unexplored genomic characteristics. From a laboratory fish mortality event off Long Island Sound, USA, we isolated a dinoflagellate, and by microscopic and molecular (18S rRNA gene; >99% identical) analyses found that it resembles G. jinhaense, hence named G. jinhaense strain AP17. Towards developing a genetic database for this dinoflagellate, a transcriptome of this species was sequenced using RNA-seq, producing 6 Gbp of data that was assembled into over 70,000 unigenes. The assembled transcriptome GC content was approximately 56% and the total Benchmarking Universal Single-Copy Orthologs for Eukaryota and Alveolata databases were approximately 50% and 57%, respectfully. Genes involved in grazing, energy generation, genome architecture, and protein synthesis, processing, and degradation were highly represented in the transcriptome. Moreover, fragments of polyketide synthase and saxitoxin genes were found but saxitoxins were not detected in high performance liquid chromatography measurements. With the first reported transcriptome for the Gyrodinium genus, this study will serve as a baseline for future Gyrodinium genomics and toxicological studies.

Finding Species-Specific Extracellular Surface-Facing Proteomes in Toxic Dinoflagellates

As a sequel to our previous report of the existence of species-specific protein/peptide expression profiles (PEPs) acquired by mass spectrometry in some dinoflagellates, we established, with the help of a plasma-membrane-impermeable labeling agent, a surface amphiesmal protein extraction method (SAPE) to label and capture species-specific surface proteins (SSSPs) as well as saxitoxins-producing-species-specific surface proteins (Stx-SSPs) that face the extracellular space (i.e., SSSPsEf and Stx-SSPsEf). Five selected toxic dinoflagellates, Alexandrium minutum, A. lusitanicum, A. tamarense, Gymnodinium catenatum, and Karenia mikimotoi, were used in this study. Transcriptomic databases of these five species were also constructed. With the aid of liquid chromatography linked-tandem mass spectrometry (LC-MS/MS) and the transcriptomic databases of these species, extracellularly facing membrane proteomes of the five different species were identified. Within these proteomes, 16 extracellular-facing and functionally significant transport proteins were found. Furthermore, 10 SSSPs and 6 Stx-SSPs were identified as amphiesmal proteins but not facing outward to the extracellular environment. We also found SSSPsEf and Stx-SSPsEf in the proteomes. The potential functional correlation of these proteins towards the production of saxitoxins in dinoflagellates and the degree of species specificity were discussed accordingly.

Transcriptome and metabolome analyses of cold and darkness-induced pellicle cysts of Scrippsiella trochoidea

BACKGROUND

Dinoflagellates are a group of unicellular organisms that are a major component of aquatic eukaryotes and important contributors to marine primary production. Nevertheless, many dinoflagellates are considered harmful algal bloom (HAB) species due to their detrimental environmental and human health impacts. Cyst formation is widely perceived as an adaptive strategy of cyst-forming dinoflagellates in response to adverse environmental conditions. Dinoflagellate cysts play critical roles in bloom dynamics. However, our insight into the underlying molecular basis of encystment is still limited. To investigate the molecular processes regulating encystment in dinoflagellates, transcriptome and metabolome investigations were performed on cold and darkness-induced pellicle cysts of Scrippsiella trochoidea.

RESULTS

No significant transcriptional response was observed at 2 h; however, massive transcriptome and metabolome reprogramming occurred at 5 h and in pellicle cysts. The gene-to-metabolite network demonstrated that the initial transformation from vegetative cells into pellicle cysts was highly energy demanding through the activation of catabolism, including glycolysis, β-oxidation, TCA cycle and oxidative phosphorylation, to cope with cold-darkness-induced stress. However, after transformation into pellicle cysts, the metabolism was greatly reduced, and various sugars, polyunsaturated fatty acids and amino acids accumulated to prolong survival. The identification of 56 differentially expressed genes (DEGs) related to signal transduction indicated that S. trochoidea received a cold-darkness signal that activated multiple signal transduction pathways, leading to encystment. The elevated expression of genes encoding enzymes involved in ROS stress suggested that pellicle cysts respond to increased oxidative stress. Several cell cycle-related genes were repressed. Intriguingly, 11 DEGs associated with sexual reproduction suggested that pellicle cysts (or some portion thereof) may be a product of sexual reproduction.

CONCLUSIONS

This study provides the first transcriptome and metabolome analyses conducted during the encystment of S. trochoidea, an event that requires complex regulatory mechanisms and impacts on population dynamics. The results reveal comprehensive molecular regulatory processes underlying life cycle regulation in dinoflagellates involving signal transduction, gene expression and metabolite profile, which will improve our ability to understand and monitor dinoflagellate blooms.

Transcriptome profiling reveals versatile dissolved organic nitrogen utilization, mixotrophy, and N conservation in the dinoflagellate Prorocentrum shikokuense under N deficiency

Harmful algal blooms formed by certain dinoflagellate species often occur when environmental nitrogen nutrients (N) are limited. However, the molecular mechanism by which dinoflagellates adapt to low N environments is poorly understood. In this study, we characterized the transcriptomic responses of Prorocentrum shikokuense to N deficiency, along with its physiological impact. Under N deficiency, P. shikokuense cultures exhibited growth inhibition, a reduction in cell size, and decreases in cellular chlorophyll a and nitrogen contents but an increase in carbon content. Accordingly, gene expression profiles indicated that carbon fixation and catabolism and fatty acid metabolism were enhanced. Transporter genes of nitrate/nitrite, ammonium, urea, and amino acids were significantly upregulated, indicating that P. shikokuense cells invest to enhance the uptake of available dissolved N. Notably, upregulated genes included those involved in endocytosis and phagosomes, evidence that P. shikokuense is a mixotrophic organism that activates phagotrophy to overcome N deficiency. Additionally, vacuolar amino acid transporters, the urea cycle, and urea hydrolysis genes were upregulated, indicating N recycling within the cells under N deficiency. Our study indicates that P. shikokuense copes with N deficiency by economizing nitrogen use and adopting multiple strategies to maximize N acquisition and reuse while maintaining carbon fixation. The remarkable low N adaptability may confer competitive advantages to P. shikokuense for forming harmful blooms in DIN-limited environments.

Recent Progress in the Study of Peroxiredoxin in the Harmful Algal Bloom Species Chattonella marina

Peroxiredoxin (Prx) is a relatively recently discovered antioxidant enzyme family that scavenges peroxides and is known to be present in organisms from biological taxa ranging from bacteria to multicellular eukaryotes, including photosynthetic organisms. Although there have been many studies of the Prx family in higher plants, green algae, and cyanobacteria, few studies have concerned raphidophytes and dinoflagellates, which are among the eukaryotic algae that cause harmful algal blooms (HABs). In our proteomic study using 2-D electrophoresis, we found a highly expressed 2-Cys peroxiredoxin (2-CysPrx) in the raphidophyte Chattonella marina var. antiqua, a species that induces mass mortality of aquacultured fish. The abundance of the C. marina 2-CysPrx enzyme was highest in the exponential growth phase, during which photosynthetic activity was high, and it then decreased by about a factor of two during the late stationary growth phase. This pattern suggested that 2-CysPrx is a key enzyme involved in the maintenance of high photosynthesis activity. In addition, the fact that the depression of photosynthesis by excessively high irradiance was more severe in the 2-CysPrx low-expression strain (wild type) than in the normal-expression strain (wild type) of C. marina suggested that 2-CysPrx played a critical role in protecting the cell from oxidative stress caused by exposure to excessively high irradiance. In the field of HAB research, estimates of growth potential have been desired to predict the population dynamics of HABs for mitigating damage to fisheries. Therefore, omics approaches have recently begun to be applied to elucidate the physiology of the growth of HAB species. In this review, we describe the progress we have made using a molecular physiological approach to identify the roles of 2-CysPrx and other antioxidant enzymes in mitigating environmental stress associated with strong light and high temperatures and resultant oxidative stress. We also describe results of a survey of expressed Prx genes and their growth-phase-dependent behavior in C. marina using RNA-seq analysis. Finally, we speculate about the function of these genes and the ecological significance of 2-CysPrx, such as its involvement in circadian rhythms and the toxicity of C. marina to fish.

Differentiating Two Closely Related Alexandrium Species Using Comparative Quantitative Proteomics

Alexandrium minutum and Alexandrium tamutum are two closely related harmful algal bloom (HAB)-causing species with different toxicity. Using isobaric tags for relative and absolute quantitation (iTRAQ)-based quantitative proteomics and two-dimensional differential gel electrophoresis (2D-DIGE), a comprehensive characterization of the proteomes of A. minutum and A. tamutum was performed to identify the cellular and molecular underpinnings for the dissimilarity between these two species. A total of 1436 proteins and 420 protein spots were identified using iTRAQ-based proteomics and 2D-DIGE, respectively. Both methods revealed little difference (10-12%) between the proteomes of A. minutum and A. tamutum, highlighting that these organisms follow similar cellular and biological processes at the exponential stage. Toxin biosynthetic enzymes were present in both organisms. However, the gonyautoxin-producing A. minutum showed higher levels of osmotic growth proteins, Zn-dependent alcohol dehydrogenase and type-I polyketide synthase compared to the non-toxic A. tamutum. Further, A. tamutum had increased S-adenosylmethionine transferase that may potentially have a negative feedback mechanism to toxin biosynthesis. The complementary proteomics approach provided insights into the biochemistry of these two closely related HAB-causing organisms. The identified proteins are potential biomarkers for organismal toxicity and could be explored for environmental monitoring.

De novo transcriptome sequencing of Paecilomyces tenuipes revealed genes involved in adenosine biosynthesis

The use of Paecilomyces tenuipes (P. tenuipes), a Chinese medicinal fungus in scientific research, is limited due to its low adenosine content. To improve adenosine production, the present study investigated the gene network of adenosine biosynthesis in P. tenuipes via transcriptome analysis. Mycelia of P. tenuipes cultured for 24 h (PT24), 102 h (PT102) and 196 h (PT192) were subjected to RNA sequencing. In total, 13,353 unigenes were obtained. Based on sequence similarity, 8,099 unigenes were annotated with known proteins. Of these 8,099 unigenes, 5,123 had functions assigned based on Gene Ontology terms while 4,158 were annotated based on the Eukaryotic Orthologous Groups database. Moreover, 1,272 unigenes were mapped to 281 Kyoto Encyclopedia of Genes and Genomes pathways. In addition, the differential gene expression of the three libraries was also performed. A total of 601, 1,658 and 628 differentially expressed genes (DEGs) were detected in PT24 vs. PT102, PT24 vs. PT192 and PT102 vs. PT192 groups, respectively. Reverse transcription‑quantitative PCR was performed to analyze the expression levels of 14 DEGs putatively associated with adenosine biosynthesis in P. tenuipes. The results showed that two DEGs were closely associated with adenosine accumulation of P. tenuipes. The present study not only provides an improved understanding of the genetic information of P. tenuipes but also the findings can be used to aid research into P. tenuipes.

Symbiotic lifestyle triggers drastic changes in the gene expression of the algal endosymbiont Breviolum minutum (Symbiodiniaceae)

Coral-dinoflagellate symbiosis underpins the evolutionary success of corals reefs. Successful exchange of molecules between the cnidarian host and the Symbiodiniaceae algae enables the mutualistic partnership. The algae translocate photosynthate to their host in exchange for nutrients and shelter. The photosynthate must traverse multiple membranes, most likely facilitated by transporters. Here, we compared gene expression profiles of cultured, free-living Breviolum minutum with those of the homologous symbionts freshly isolated from the sea anemone Exaiptasia diaphana, a widely used model for coral hosts. Additionally, we assessed expression levels of a list of candidate host transporters of interest in anemones with and without symbionts. Our transcriptome analyses highlight the distinctive nature of the two algal life stages, with many gene expression level changes correlating to the different morphologies, cell cycles, and metabolisms adopted in hospite versus free-living. Morphogenesis-related genes that likely underpin the metamorphosis process observed when symbionts enter a host cell were up-regulated. Conversely, many down-regulated genes appear to be indicative of the protective and confined nature of the symbiosome. Our results emphasize the significance of transmembrane transport to the symbiosis, and in particular of ammonium and sugar transport. Further, we pinpoint and characterize candidate transporters-predicted to be localized variously to the algal plasma membrane, the host plasma membrane, and the symbiosome membrane-that likely serve pivotal roles in the interchange of material during symbiosis. Our study provides new insights that expand our understanding of the molecular exchanges that underpin the cnidarian-algal symbiotic relationship.

Transcriptomic Profile and Sexual Reproduction-Relevant Genes of Alexandrium minutum in Response to Nutritional Deficiency

Alexandrium minutum is a typical marine toxic dinoflagellate responsible for producing paralytic shellfish poisoning (PSP) toxins. Until now, we know little about the genomic information of A. minutum, so a transcriptome study was conducted to clarify the physiological adaptations related to nutritional deficiency. Here, we performed RNA-Seq analysis to assess the gene expression patterns of A. minutum under N and P deficient conditions for 0 (control), 6, and 72 h. Main differences between the control and experimental groups were observed in hydrolase activity and fatty acid, lipid, protein, and P metabolism. Activities of photosystem I (PSI) and PSII were significantly down-regulated, and the endocytosis pathway (clathrin-dependent endocytosis) was significantly enriched under N and P stress compared with the control, indicating that A. minutum shifts its trophy pattern under N and P stress. We also identified several unigenes related to the process of sexual reproduction, including sex determination, sperm-egg recognition, sex differentiation, mating, and fertilization. Approximately 50% of the successfully annotated unigenes were differentially expressed between the short-term stimulated sample (6 h) and control (R). However, the expression level of most unigenes returned to normal levels after 72 h, indicating that N and P stress plays a limited role in the induction of sexual reproduction. Furthermore, the quantitative real-time PCR (qRT-PCR) results of the five representative sex-related unigenes were consistent with sequencing data, which confirmed the authenticity of transcriptomic analysis. Also, qRT-PCR analysis showed that the long and short form transcripts of the saxitoxin biosynthesis gene (sxtA) were down-regulated under the nutrient deficient condition compared with the control, indicating that N and P stress regulates sxtA expression. Overall, transcriptome analysis of A. minutum revealed that N and P deficiency induced responses associated with stress response, photosynthetic efficiency, toxin biosynthesis, and sexual reproduction. Our data indicate that algae change their trophic modes (to facultative mixotrophy) and related physiological reactions under stress conditions; this possibly represents an ecological adaption strategy in the algal life cycle.

Transcriptome sequencing of a toxic dinoflagellate, Karenia mikimotoi subjected to stress from solar ultraviolet radiation

Solar ultraviolet radiation (UVR) is a stress factor in aquatic environments and may act directly or indirectly on orgnisms in the upper layers of the water column. However, UVR effects are usually species-specific and difficult to extrapolate. Here we use the HAB-forming, toxic dinoflagellate Karenia mikimotoi (which was found to be relatively resistant in previous studies) to investigate its transcriptional responses to a one-week UVR exposure. For this, batch cultures of K. mikimotoi were grown with and without UVR, and their transcriptomes (generated via RNAseq technology) were compared. RNA-seq generated 45.31 million reads, which were further assembled to 202600 unigenes (>300bp). Among these, ca. 61% were annotated with NCBI, NR, GO, KOG, PFAM, Swiss-Prot, and KEGG database. Transcriptomic analysis revealed 722 differentially expressed unigenes (DEGs, defined as being within a |log2 fold change| ≥ 2 and padj < 0.05) responding to solar UVR, which were only 0.36% of all unigenes. 716 unigenes were down-regulated, and only 6 unigenes were up-regulated in the UVR compared to non-UVR treatment. KEGG pathway further analysis revealed DEGs were involved in the different pathway; genes involved in the ribosome, endocytosis and steroid biosynthesis pathways were highly down-regulated, but this was not the case for those involved in the energy metabolisms (including photosynthesis, oxidative phosphorylation) which may contribute to the sustainable growth observed in UVR treatment. The up-regulated expression of both zinc-finger proteins (ZFPs) and ribosomal protein L11 (RPL11) may be one of the acclimated mechanisms against UVR. In addition, this work identified down-regulated genes involved in fatty acid degradation and the hydrophobic branched chain amino acids (e.g., Valine, leucine, and isoleucine), which act as structural components of cell membranes modulating lipid homeostasis or turnover. In conclusion, the present study suggests that the toxic dinoflagellate K. mikimotoi has limited transcriptomic regulation but confirms that it appears as a tolerant species in response to solar UVR. These findings expand current knowledge of gene expression in HAB-forming species in response to natural environment factors such as solar radiation.

Omics Analysis for Dinoflagellates Biology Research

Dinoflagellates are important primary producers for marine ecosystems and are also responsible for certain essential components in human foods. However, they are also notorious for their ability to form harmful algal blooms, and cause shellfish poisoning. Although much work has been devoted to dinoflagellates in recent decades, our understanding of them at a molecular level is still limited owing to some of their challenging biological properties, such as large genome size, permanently condensed liquid-crystalline chromosomes, and the 10-fold lower ratio of protein to DNA than other eukaryotic species. In recent years, omics technologies, such as genomics, transcriptomics, proteomics, and metabolomics, have been applied to the study of marine dinoflagellates and have uncovered many new physiological and metabolic characteristics of dinoflagellates. In this article, we review recent application of omics technologies in revealing some of the unusual features of dinoflagellate genomes and molecular mechanisms relevant to their biology, including the mechanism of harmful algal bloom formations, toxin biosynthesis, symbiosis, lipid biosynthesis, as well as species identification and evolution. We also discuss the challenges and provide prospective further study directions and applications of dinoflagellates.

RNA-Seq Analysis Reveals Genes Related to Photoreception, Nutrient Uptake, and Toxicity in a Noxious Red-Tide Raphidophyte Chattonella antiqua

Aquaculture industries are under threat from noxious red tides, but harm can be mitigated by precautions such as early harvesting and restricting fish feeding to just before the outbreak of a red tide. Therefore, accurate techniques for forecasting red-tide outbreaks are strongly needed. Omics analyses have the potential to expand our understanding of the eco-physiology of these organisms at the molecular level, and to facilitate identification of molecular markers for forecasting their population dynamics and occurrence of damages to fisheries. Red tides of marine raphidophytes, especially Chattonella species, often extensively harm aquaculture industries in regions with a temperate climate around the world. A red tide of Chattonella tends to develop just after an input of nutrients along the coast. Chattonella displays diurnal vertical migration regulated by a weak blue light, so it photosynthesizes in the surface layer during the daytime and takes up nutrients in the bottom layer during the nighttime. Superoxide produced by Chattonella cells is a strong candidate for the cause of its toxicity to bacteria and fishes. Here we conducted mRNA-seq of Chattonella antiqua to identify genes with functions closely related to the dynamics of the noxious red tide, such as photosynthesis, photoreception, nutrient uptake, and superoxide production. The genes related to photosynthetic pigment biosynthesis and nutrient uptake had high similarity with those of model organisms of plants and algae and other red-tide microalgae. We identified orthologous genes of photoreceptors such as aureochrome (newly five genes), the cryptochrome/photolyase (CRY/PHR) family (6-4PHR, plant CRY or cyclobutane pyrimidine dimer [CPD] Class III, CPD Class II, and CRY-DASH), and phytochrome (four genes), which regulate various physiological processes such as flagellar motion and cell cycle in model organisms. Six orthologous genes of NADPH oxidase, which produces superoxide on the cell membrane, were found and divided into two types: one with 5–6 transmembrane domains and another with 11 transmembrane domains. The present study should open the way for analyzing the eco-physiological features of marine raphidophytes at the molecular level.

Detection, characterization and expression dynamics of histone proteins in the dinoflagellate Alexandrium pacificum during growth regulation

Histones are the most abundant proteins associated with eukaryotic nuclear DNA. The exception is dinoflagellates, which have histone protein expression that is mostly reported to be below detectable levels. In this study, we investigated the presence of histone proteins and their functions in the dinoflagellate, Alexandrium pacificum. Histone protein sequences were analyzed, focusing on phylogenetic analysis and histone code. Histone expression was analyzed during the cell cycle and under nutritionally enhanced conditions using quantitative-PCR and western blots. Acid-soluble proteins were subjected to mass spectrometry analysis. To our knowledge, this is the first report of immunological detection of histone proteins (H2B and H4) in any dinoflagellate species. Absolute quantification of histone transcript in activily dividing cells revealed significant transcription in cells. The stable expression of histones during the cell cycle suggested that the histone genes in A. pacificum belonged to a replication-independent class and appeared to have a limited role in DNA packaging. The conservation of numerous post-translationally modified residues of multiple histone variants and differential expression of histones under nutritionally enhanced conditions suggested their functional significance in dinoflagellates. However, we detected histone H2B protein only via mass spectrometry. Histone-like protein was identified as most abundant acid-soluble protein of the cells.

A brief history of circannual time

Innate circannual timing is an ancestral trait that first evolved in free-living eukaryotic cells some 2000 million years ago, with marine algae of the genus Allexandrium providing a living unicellular model. This species shows the primitive trait of 'alternation of generations', where the organism alternates between fast replicating vegetative cells in the summer and a dormant cystic cell over the winter. The resistant cysts sink into the cold ocean sediments. Remarkably, excystment in spring is governed by an endogenous circannual timing mechanism. Thus, a tiny, short-lived unicell can utilise a circannual clock as part of the life-history programme of the species. Innate timing allows for major adjustments in physiology and behaviour in anticipation of the seasons, and provides an internalised sense of seasonal time for the many species where standard environmental cues are weak or ambiguous. This is a highly adaptive strategy irrespective of the size and longevity of an organism. Circannual rhythms are expressed by a diverse range of organisms, from flowering plants to mammals, interwoven into the life-history programme of each species, being a consequence of forever living in a periodic world. In complex vertebrates, the early division of the zygote potentially carries circannual timer genes into all progeny cells and tissues. This supports the concept of a 'clock-shop' where cell-autonomous long-term rhythms are generated in each tissue, orchestrated by a central circannual pacemaker system. This is analogous to the organisation of the circadian timing system. For the circannual time-scale, specialised thyrotroph cells located in the pars tuberalis of the pituitary gland and adjacent tanycyte cells located in the ependymal wall of the third cerebral ventricle of the brain act as putative central circannual pacemakers. At a molecular level, epigenetically regulated, cyclical remodelling of chromatin, which determines whether specific circannual timer genes are transcriptionally active, or not, is considered to drive the oscillation between the summer and winter phenotypes.

Distinctive Nuclear Features of Dinoflagellates with A Particular Focus on Histone and Histone-Replacement Proteins

Dinoflagellates are important eukaryotic microorganisms that play critical roles as producers and grazers, and cause harmful algal blooms. The unusual nuclei of dinoflagellates "dinokaryon" have led researchers to investigate their enigmatic nuclear features. Their nuclei are unusual in terms of their permanently condensed nucleosome-less chromatin, immense genome, low protein to DNA ratio, guanine-cytosine rich methylated DNA, and unique mitosis process. Furthermore, dinoflagellates are the only known group of eukaryotes that apparently lack histone proteins. Over the course of evolution, dinoflagellates have recruited other proteins, e.g., histone-like proteins (HLPs), from bacteria and dinoflagellates/viral nucleoproteins (DVNPs) from viruses as histone substitutes. Expression diversity of these nucleoproteins has greatly influenced the chromatin structure and gene expression regulation in dinoflagellates. Histone replacement proteins (HLPs and DVNPs) are hypothesized to perform a few similar roles as histone proteins do in other eukaryotes, i.e., gene expression regulation and repairing DNA. However, their role in bulk packaging of DNA is not significant as low amounts of proteins are associated with the gigantic genome. This review intends to summarize the discoveries encompassing unique nuclear features of dinoflagellates, particularly focusing on histone and histone replacement proteins. In addition, a comprehensive view of the evolution of dinoflagellate nuclei is presented.

Transcriptional responses of Acropora hyacinthus embryo under the benzo(a)pyrene stress by deep sequencing

Coral embryos are a critical and sensitive period for the early growth and development of coral. Benzo(a)pyrene (BaP) is widely distributed in the ocean and has strong toxicity, but there is little information on the toxic effects to coral embryos exposed to this widespread environmental contaminant. Thus, in this study, we utilized the Illumina Hiseq™ 4000 platform to explore the gene response of Acropora hyacinthus embryos under the BaP stress. A total of 130,042 Unigenes were obtained and analyzed, and approximately 37.67% of those matched with sequences from four different species. In total, 2606 Unigenes were up-regulated, and 3872 Unigenes were down-regulated. After Gene Ontology (GO) annotation, the results show that the "cellular process" and "metabolic process" were leading in the category of biological processes, which the "binding" and "catalytic activity" were the most abundant subcategories in molecular function. Based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, the most differentially expressed genes (DEGs) were enriched, as well as down-regulated in the pathways of oxidative phosphorylation, metabolism of xenobiotics, immune-related genes, apoptosis and human disease genes. At the same time, 388,197 of Single-nucleotide Polymorphisms (SNPs) and 6164 of Simple Sequence Repeats (SSRs) were obtained, which can be served as the richer and more valuable SSRs molecular markers in the future. The results of this study can help to better understand the toxicological mechanism of coral embryo exposed to BaP, and it is also essential for the protection and restoration of coral reef ecosystem in the future.

6.0 K microarray reveals differential transcriptomic responses in the dinoflagellate Prorocentrum minimum exposed to polychlorinated biphenyl (PCB)

Endocrine disrupting chemicals (EDCs) have toxic effects on algae; however, their molecular genomic responses have not been sufficiently elucidated. Here, we evaluated genome-scaled responses of the dinoflagellate alga Prorocentrum minimum exposed to an EDC, polychlorinated biphenyl (PCB), using a 6.0 K microarray. Based on two-fold change cut-off, we identified that 609 genes (∼10.2%) responded to the PCB treatment. KEGG pathway analysis showed that differentially expressed genes (DEGs) were related to ribosomes, biosynthesis of amino acids, spliceosomes, and cellular processes. Many DEGs were involved in cell cycle progression, apoptosis, signal transduction, ion binding, and cellular transportation. In contrast, only a few genes related to photosynthesis and oxidative stress were expressed in response to PCB exposure. This was supported by that fact that there were no obvious changes in the photosynthetic efficiency and reactive oxygen species (ROS) production. These results suggest that PCB might not cause chloroplast and oxidative damage, but could lead to cell cycle arrest and apoptosis. In addition, various signal transduction and transport pathways might be disrupted in the cells, which could further contribute to cell death. These results expand the genomic understanding of the effects of EDCs on this dinoflagellate protist.

Transcriptomic Analyses of Scrippsiella trochoidea Reveals Processes Regulating Encystment and Dormancy in the Life Cycle of a Dinoflagellate, with a Particular Attention to the Role of Abscisic Acid

Due to the vital importance of resting cysts in the biology and ecology of many dinoflagellates, a transcriptomic investigation on Scrippsiella trochoidea was conducted with the aim to reveal the molecular processes and relevant functional genes regulating encystment and dormancy in dinoflagellates. We identified via RNA-seq 3,874 (out of 166,575) differentially expressed genes (DEGs) between resting cysts and vegetative cells; a pause of photosynthesis (confirmed via direct measurement of photosynthetic efficiency); an active catabolism including β-oxidation, glycolysis, glyoxylate pathway, and TCA in resting cysts (tested via measurements of respiration rate); 12 DEGs encoding meiotic recombination proteins and members of MEI2-like family potentially involved in sexual reproduction and encystment; elevated expressions in genes encoding enzymes responding to pathogens (chitin deacetylase) and ROS stress in cysts; and 134 unigenes specifically expressed in cysts. We paid particular attention to genes pertaining to phytohormone signaling and identified 4 key genes regulating abscisic acid (ABA) biosynthesis and catabolism, with further characterization based on their full-length cDNA obtained via RACE-PCR. The qPCR results demonstrated elevated biosynthesis and repressed catabolism of ABA during the courses of encystment and cyst dormancy, which was significantly enhanced by lower temperature (4 ± 1°C) and darkness. Direct measurements of ABA using UHPLC-MS/MS and ELISA in vegetative cells and cysts both fully supported qPCR results. These results collectively suggest a vital role of ABA in regulating encystment and maintenance of dormancy, akin to its function in seed dormancy of higher plants. Our results provided a critical advancement in understanding molecular processes in resting cysts of dinoflagellates.

Whole Transcriptomic Analysis Provides Insights into Molecular Mechanisms for Toxin Biosynthesis in a Toxic Dinoflagellate Alexandrium catenella (ACHK-T)

Paralytic shellfish toxins (PSTs), a group of neurotoxic alkaloids, are the most potent biotoxins for aquatic ecosystems and human health. Marine dinoflagellates and freshwater cyanobacteria are two producers of PSTs. The biosynthesis mechanism of PSTs has been well elucidated in cyanobacteria; however, it remains ambiguous in dinoflagellates. Here, we compared the transcriptome profiles of a toxin-producing dinoflagellate Alexandrium catenella (ACHK-T) at different toxin biosynthesis stages within the cell cycle using RNA-seq. The intracellular toxin content increased gradually in the middle G1 phase and rapidly in the late G1 phase, and then remained relatively stable in other phases. Samples from four toxin biosynthesis stages were selected for sequencing, and finally yielded 110,370 unigenes, of which 66,141 were successfully annotated in the known databases. An analysis of differentially expressed genes revealed that 2866 genes altered significantly and 297 were co-expressed throughout the four stages. These genes participated mainly in protein metabolism, carbohydrate metabolism, and the oxidation-reduction process. A total of 138 homologues of toxin genes were identified, but they altered insignificantly among different stages, indicating that toxin biosynthesis might be regulated translationally or post-translationally. Our results will serve as an important transcriptomic resource to characterize key molecular processes underlying dinoflagellate toxin biosynthesis.

Comprehensive transcriptome analysis provides new insights into nutritional strategies and phylogenetic relationships of chrysophytes

Background

Chrysophytes are protist model species in ecology and ecophysiology and important grazers of bacteria-sized microorganisms and primary producers. However, they have not yet been investigated in detail at the molecular level, and no genomic and only little transcriptomic information is available. Chrysophytes exhibit different trophic modes: while phototrophic chrysophytes perform only photosynthesis, mixotrophs can gain carbon from bacterial food as well as from photosynthesis, and heterotrophs solely feed on bacteria-sized microorganisms. Recent phylogenies and megasystematics demonstrate an immense complexity of eukaryotic diversity with numerous transitions between phototrophic and heterotrophic organisms. The question we aim to answer is how the diverse nutritional strategies, accompanied or brought about by a reduction of the plasmid and size reduction in heterotrophic strains, affect physiology and molecular processes.

Results

We sequenced the mRNA of 18 chrysophyte strains on the Illumina HiSeq platform and analysed the transcriptomes to determine relations between the trophic mode (mixotrophic vs. heterotrophic) and gene expression. We observed an enrichment of genes for photosynthesis, porphyrin and chlorophyll metabolism for phototrophic and mixotrophic strains that can perform photosynthesis. Genes involved in nutrient absorption, environmental information processing and various transporters (e.g., monosaccharide, peptide, lipid transporters) were present or highly expressed only in heterotrophic strains that have to sense, digest and absorb bacterial food. We furthermore present a transcriptome-based alignment-free phylogeny construction approach using transcripts assembled from short reads to determine the evolutionary relationships between the strains and the possible influence of nutritional strategies on the reconstructed phylogeny. We discuss the resulting phylogenies in comparison to those from established approaches based on ribosomal RNA and orthologous genes. Finally, we make functionally annotated reference transcriptomes of each strain available to the community, significantly enhancing publicly available data on Chrysophyceae.

Conclusions

Our study is the first comprehensive transcriptomic characterisation of a diverse set of Chrysophyceaen strains. In addition, we showcase the possibility of inferring phylogenies from assembled transcriptomes using an alignment-free approach. The raw and functionally annotated data we provide will prove beneficial for further examination of the diversity within this taxon. Our molecular characterisation of different trophic modes presents a first such example.

The Mechanism of Diarrhetic Shellfish Poisoning Toxin Production in Prorocentrum spp.: Physiological and Molecular Perspectives

Diarrhetic shellfish poisoning (DSP) is a gastrointestinal disorder caused by the consumption of seafood contaminated with okadaic acid (OA) and dinophysistoxins (DTXs). OA and DTXs are potent inhibitors of protein phosphatases 2A, 1B, and 2B, which may promote cancer in the human digestive system. Their expression in dinoflagellates is strongly affected by nutritional and environmental factors. Studies have indicated that the level of these biotoxins is inversely associated with the growth of dinoflagellates at low concentrations of nitrogen or phosphorus, or at extreme temperature. However, the presence of leucine or glycerophosphate enhances both growth and cellular toxin level. Moreover, the presence of ammonia and incubation in continuous darkness do not favor the toxin production. Currently, studies on the mechanism of this biotoxin production are scant. Full genome sequencing of dinoflagellates is challenging because of the massive genomic size; however, current advanced molecular and omics technologies may provide valuable insight into the biotoxin production mechanism and novel research perspectives on microalgae. This review presents a comprehensive analysis on the effects of various nutritional and physical factors on the OA and DTX production in the DSP toxin-producing Prorocentrum spp. Moreover, the applications of the current molecular technologies in the study on the mechanism of DSP toxin production are discussed.

Genome-wide analysis of transcription and photosynthesis inhibition in the harmful dinoflagellate Prorocentrum minimum in response to the biocide copper sulfate

Copper is an essential trace metal for organisms; however, excess copper may damage cellular processes. Their efficiency and physiological effects of biocides have been well documented; however, molecular transcriptome responses to biocides are insufficiently studied. In the present study, a 6.0K oligonucleotide chip was developed to investigate the molecular responses of the harmful dinoflagellate Prorocentrum minimum to copper sulfate (CuSO₄) treatment. The results revealed that 515 genes (approximately 8.6%) responded to CuSO₄, defined as being within a 2-fold change. Further, KEGG pathway analysis showed that differentially expressed genes (DEGs) were involved in ribosomal function, RNA transport, carbon metabolism, biosynthesis of amino acids, photosystem maintenance, and other cellular processes. Among the DEGs, 49 genes were related to chloroplasts and mitochondria. Furthermore, the genes involved in the RAS signaling pathway, MAPK signaling pathway, and transport pathways were identified. An additional experiment showed that the photosynthesis efficiency decreased considerably, and reactive oxygen species (ROS) production increased in P. minimum after CuSO₄ exposure. These results suggest that CuSO₄ caused cellular oxidative stress in P. minimum, affecting the ribosome and mitochondria, and severely damaged the photosystem. These effects may potentially lead to cell death, although the dinoflagellate has developed a complex signal transduction process to combat copper toxicity.

Transcriptomic profiles reveal the genome-wide responses of the harmful dinoflagellate Cochlodinium polykrikoides when exposed to the algicide copper sulfate

Background

Harmful algal blooms (HABs) caused by the dinoflagellate Cochlodinium polykrikoides lead to severe environmental impacts in oceans worldwide followed by huge economic losses. Algicide agent copper sulfate (CuSO₄) is regard as an economical and effective agent for HABs mitigation; its biochemical and physiological effects were revealed in C. polykrikoides. However, molecular mechanisms of CuSO₄ effect on the C. polykrikoides, even other HAB species, have not been investigated. The present study investigated the transcriptional response of C. polykrikoides against CuSO₄ treatments, with the aim of providing certain molecular mechanism of CuSO₄ effect on the C. polykrikoides blooms.

Results

RNA-seq generated 173 million reads, which were further assembled to 191,212 contigs. 43.3 %, 33.9 %, and 15.6 % of contigs were annotated with NCBI NR, GO, and KEGG database, respectively. Transcriptomic analysis revealed 20.6 % differential expressed contigs, which grouped into 8 clusters according to K-means clustering analysis, responding to CuSO₄; 848 contigs were up-regulated and 746 contigs were down-regulated more than 2-fold changes from 12 h to 48 h exposure. KEGG pathway analysis of eukaryotic homologous genes revealed the differentially expressed genes (DEGs) were involved in diverse pathway; amongst, the genes involved in the translation, spliceosome, and/or signal transduction genes were highly regulated. Most of photosystem related genes were down-regulated and most of mitochondria related genes were up-regulated. In addition, the genes involved in the copper ion binding or transporting and antioxidant systems were identified. Measurement of chlorophyll fluorescence showed that photosynthesis was significantly inhibited by CuSO₄ exposure.

Conclusions

This study reported the first transcriptome of the C. polykrikoides. The widely differential expressed photosystem genes suggested photosynthetic machinery were severely affected, and may further contribute to the cell death. Furthermore, gene translation and transcription processes may be disrupted, inhibiting cell growth and proliferation, and possibly accelerating cell death. However, antioxidant systems resistant to CuSO₄ caused stress; mitochondrion may compensate for photosynthesis efficiency decreasing caused energy deficiency. In addition, various signal transduction pathways may be involved in the CuSO₄ induced regulation network in the C. polykrikoides. These data provide the potential transcriptomic mechanism to explain the algicide CuSO₄ effect on the harmful dinoflagellate C. polykrikoides.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2341-3) contains supplementary material, which is available to authorized users.

Diversity and Divergence of Dinoflagellate Histone Proteins

Histone proteins and the nucleosomal organization of chromatin are near-universal eukaroytic features, with the exception of dinoflagellates. Previous studies have suggested that histones do not play a major role in the packaging of dinoflagellate genomes, although several genomic and transcriptomic surveys have detected a full set of core histone genes. Here, transcriptomic and genomic sequence data from multiple dinoflagellate lineages are analyzed, and the diversity of histone proteins and their variants characterized, with particular focus on their potential post-translational modifications and the conservation of the histone code. In addition, the set of putative epigenetic mark readers and writers, chromatin remodelers and histone chaperones are examined. Dinoflagellates clearly express the most derived set of histones among all autonomous eukaryote nuclei, consistent with a combination of relaxation of sequence constraints imposed by the histone code and the presence of numerous specialized histone variants. The histone code itself appears to have diverged significantly in some of its components, yet others are conserved, implying conservation of the associated biochemical processes. Specifically, and with major implications for the function of histones in dinoflagellates, the results presented here strongly suggest that transcription through nucleosomal arrays happens in dinoflagellates. Finally, the plausible roles of histones in dinoflagellate nuclei are discussed.

Identification of microRNAs in the Toxigenic Dinoflagellate Alexandrium catenella by High-Throughput Illumina Sequencing and Bioinformatic Analysis

Micro-ribonucleic acids (miRNAs) are a large group of endogenous, tiny, non-coding RNAs consisting of 19–25 nucleotides that regulate gene expression at either the transcriptional or post-transcriptional level by mediating gene silencing in eukaryotes. They are considered to be important regulators that affect growth, development, and response to various stresses in plants. Alexandrium catenella is an important marine toxic phytoplankton species that can cause harmful algal blooms (HABs). To date, identification and function analysis of miRNAs in A. catenella remain largely unexamined. In this study, high-throughput sequencing was performed on A. catenella to identify and quantitatively profile the repertoire of small RNAs from two different growth phases. A total of 38,092,056 and 32,969,156 raw reads were obtained from the two small RNA libraries, respectively. In total, 88 mature miRNAs belonging to 32 miRNA families were identified. Significant differences were found in the member number, expression level of various families, and expression abundance of each member within a family. A total of 15 potentially novel miRNAs were identified. Comparative profiling showed that 12 known miRNAs exhibited differential expression between the lag phase and the logarithmic phase. Real-time quantitative RT-PCR (qPCR) was performed to confirm the expression of two differentially expressed miRNAs that were one up-regulated novel miRNA (aca-miR-3p-456915), and one down-regulated conserved miRNA (tae-miR159a). The expression trend of the qPCR assay was generally consistent with the deep sequencing result. Target predictions of the 12 differentially expressed miRNAs resulted in 1813target genes. Gene ontology (GO) analysis and the Kyoto Encyclopedia of Genes and Genomes pathway database (KEGG) annotations revealed that some miRNAs were associated with growth and developmental processes of the alga. These results provide insights into the roles that miRNAs play in the growth of A. catenella, and they provide the basis for further studies of the molecular mechanisms that underlie bloom growth in red tides species.

De novo transcriptome of the European brittle star Amphiura filiformis pluteus larvae

BACKGROUND

In non-classical model species, Next Generation Sequencing increases the ability to analyze the expression of transcripts/genes. In this study, paired-end Illumina HiSeq sequencing technology has been employed to describe a larval transcriptome generated from 64 h post-fertilization pluteus larvae of the brittle star Amphiura filiformis. We focused our analysis on the detection of actors involved in the opsin based light perception, respectively the opsins and the phototransduction actors.

METHODS & RESULTS

In this research, about 47 million high quality reads were generated and 86,572 total unigenes were predicted after de novo assembly. Of all the larval unigenes, 18% show significant matches with reference online databases. 46% of annotated larval unigenes were significantly similar to transcripts from the purple sea urchin. COG, GO and KEGG analyses were performed on predicted unigenes. Regarding the opsin-based photoreception process, even if possible actors of ciliary and rhabdomeric phototransduction cascades were detected, no ciliary or rhabdomeric opsin was identified in these larvae. Additionally, partial non-visual RGR (retinal G protein coupled receptor) opsin mRNAs were identified,possibly indicating the presence of visual cycle reaction in early pluteus larvae. The eye morphogene Pax 6 was also identified in the pluteus transcriptome.

CONCLUSIONS

Contrary to sea-urchin larvae, brittle star larvae appear to be characterized by an absence of visual-like opsins. These RNA-seq data also provide a useful resource for the echinoderm research community and researchers with an interest in larval biology.

Identification of the Full 46 Cytochrome P450 (CYP) Complement and Modulation of CYP Expression in Response to Water-Accommodated Fractions of Crude Oil in the Cyclopoid Copepod Paracyclopina nana

The 46 cytochrome P450 (CYP) gene superfamily was identified in the marine copepod Paracyclopina nana after searching an RNA-seq database and comparing it with other copepod CYP gene families. To annotate the 46 Pn-CYP genes, a phylogenetic analysis of CYP genes was performed using a Bayesian method. Pn-CYP genes were separated into five different clans: CYP2, CYP3, CYP20, CYP26, and mitochondrial. Among these, the principal Pn-CYP genes involved in detoxification were identified by comparing them with those of the copepod Tigriopus japonicus and were examined with respect to their responses to exposure to a water-accommodated fraction (WAF) of crude oil and to the alkylated forms of two polycyclic aromatic hydrocarbons (PAHs; phenanthrene and fluorene). The expression of two Pn-CYP3027 genes (CYP3027F1 and CYP3027F2) was increased in response to WAF exposure and also was upregulated in response to the two alkylated PAHs. In particular, Pn-CYP3027F2 showed the most notable increase in response to 80% WAF exposure. These two responsive CYP genes (Pn-CYP3027F1 and CYP3027F2) were also phylogenetically clustered into the same clade of the WAF- and alkylated PAH-specific CYP genes of the copepod T. japonicus, suggesting that these CYP genes would be those chiefly involved in detoxification in response to WAF exposure in copepods. In this paper, we provide information on the copepod P. nana CYP gene superfamily and also speculate on its potential role in the detoxification of PAHs in marine copepods. Despite the nonlethality of WAF, Pn-CYP3027F2 was rapidly and significantly upregulated in response to WAF that may serve as a useful biomarker of 40% or higher concentration of WAF exposure. This paper will be helpful to better understand the molecular mechanistic events underlying the metabolism of environmental toxicants in copepods.

Unfolding the secrets of coral-algal symbiosis

Dinoflagellates from the genus Symbiodinium form a mutualistic symbiotic relationship with reef-building corals. Here we applied massively parallel Illumina sequencing to assess genetic similarity and diversity among four phylogenetically diverse dinoflagellate clades (A, B, C and D) that are commonly associated with corals. We obtained more than 30,000 predicted genes for each Symbiodinium clade, with a majority of the aligned transcripts corresponding to sequence data sets of symbiotic dinoflagellates and <2% of sequences having bacterial or other foreign origin. We report 1053 genes, orthologous among four Symbiodinium clades, that share a high level of sequence identity to known proteins from the SwissProt (SP) database. Approximately 80% of the transcripts aligning to the 1053 SP genes were unique to Symbiodinium species and did not align to other dinoflagellates and unrelated eukaryotic transcriptomes/genomes. Six pathways were common to all four Symbiodinium clades including the phosphatidylinositol signaling system and inositol phosphate metabolism pathways. The list of Symbiodinium transcripts common to all four clades included conserved genes such as heat shock proteins (Hsp70 and Hsp90), calmodulin, actin and tubulin, several ribosomal, photosynthetic and cytochrome genes and chloroplast-based heme-containing cytochrome P450, involved in the biosynthesis of xanthophylls. Antioxidant genes, which are important in stress responses, were also preserved, as were a number of calcium-dependent and calcium/calmodulin-dependent protein kinases that may play a role in the establishment of symbiosis. Our findings disclose new knowledge about the genetic uniqueness of symbiotic dinoflagellates and provide a list of homologous genes important for the foundation of coral-algal symbiosis.

The alveolate translation initiation factor 4E family reveals a custom toolkit for translational control in core dinoflagellates

Background

Dinoflagellates are eukaryotes with unusual cell biology and appear to rely on translational rather than transcriptional control of gene expression. The eukaryotic translation initiation factor 4E (eIF4E) plays an important role in regulating gene expression because eIF4E binding to the mRNA cap is a control point for translation. eIF4E is part of an extended, eukaryote-specific family with different members having specific functions, based on studies of model organisms. Dinoflagellate eIF4E diversity could provide a mechanism for dinoflagellates to regulate gene expression in a post-transcriptional manner. Accordingly, eIF4E family members from eleven core dinoflagellate transcriptomes were surveyed to determine the diversity and phylogeny of the eIF4E family in dinoflagellates and related lineages including apicomplexans, ciliates and heterokonts.

Results

The survey uncovered eight to fifteen (on average eleven) different eIF4E family members in each core dinoflagellate species. The eIF4E family members from heterokonts and dinoflagellates segregated into three clades, suggesting at least three eIF4E cognates were present in their common ancestor. However, these three clades are distinct from the three previously described eIF4E classes, reflecting diverse approaches to a central eukaryotic function. Heterokonts contain four clades, ciliates two and apicomplexans only a single recognizable eIF4E clade. In the core dinoflagellates, the three clades were further divided into nine sub-clades based on the phylogenetic analysis and species representation. Six of the sub-clades included at least one member from all eleven core dinoflagellate species, suggesting duplication in their shared ancestor. Conservation within sub-clades varied, suggesting different selection pressures.

Conclusions

Phylogenetic analysis of eIF4E in core dinoflagellates revealed complex layering of duplication and conservation when compared to other eukaryotes. Our results suggest that the diverse eIF4E family in core dinoflagellates may provide a toolkit to enable selective translation as a strategy for controlling gene expression in these enigmatic eukaryotes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0301-9) contains supplementary material, which is available to authorized users.

De novo transcriptome analysis of a medicinal fungi Phellinus linteus and identification of SSR markers

The aim of this study was to facilitate gene discovery for functional genome studies and to identify simple sequence repeat (SSR) markers for molecular-assisted selection in Phellinus linteus. The transcriptome of Phellinus linteus was sequenced using а high-throughput RNA sequencing system - the Illumina Hiseq 2000. A total of 16,383,818 clean sequencing reads, 35,532 contigs and 25,811 unigenes were postulated. Based on similarity searches with known proteins, 19,350 genes (74.97% of the unigenes) were annotated. In the present research, 19,266, 10,978 and 7831 unigenes were mapped in Nr, Swiss-Prot and clusters of orthologous groups (COG) classifications, respectively. Of all unigenes, 6845 were categorized into three functional groups, namely biological process, cellular components and molecular function and 11,088 were annotated to 108 pathways by searching the Kyoto Encyclopedia of Genes and Genomes pathway database. A total of 1129 SSRs were identified in these unigenes. In addition, 23 candidate genes, potentially involved in sterol biosynthesis, were identified and were worthy of further investigation.

Comparative transcriptome analysis of a toxin-producing dinoflagellate Alexandrium catenella and its non-toxic mutant

The dinoflagellates and cyanobacteria are two major kingdoms of life producing paralytic shellfish toxins (PSTs), a large group of neurotoxic alkaloids causing paralytic shellfish poisonings around the world. In contrast to the well elucidated PST biosynthetic genes in cyanobacteria, little is known about the dinoflagellates. This study compared transcriptome profiles of a toxin-producing dinoflagellate, Alexandrium catenella (ACHK-T), and its non-toxic mutant form (ACHK-NT) using RNA-seq. All clean reads were assembled de novo into a total of 113,674 unigenes, and 66,812 unigenes were annotated in the known databases. Out of them, 35 genes were found to express differentially between the two strains. The up-regulated genes in ACHK-NT were involved in photosynthesis, carbon fixation and amino acid metabolism processes, indicating that more carbon and energy were utilized for cell growth. Among the down-regulated genes, expression of a unigene assigned to the long isoform of sxtA, the initiator of toxin biosynthesis in cyanobacteria, was significantly depressed, suggesting that this long transcript of sxtA might be directly involved in toxin biosynthesis and its depression resulted in the loss of the ability to synthesize PSTs in ACHK-NT. In addition, 101 putative homologs of 12 cyanobacterial sxt genes were identified, and the sxtO and sxtZ genes were identified in dinoflagellates for the first time. The findings of this study should shed light on the biosynthesis of PSTs in the dinoflagellates.