Comparative genomic and transcriptomic characterization of the toxigenic marine dinoflagellate Alexandrium ostenfeldii

Unveiling the functional nature of retrogenes in dinoflagellates

Retroposition is a gene duplication mechanism that uses RNA molecules as intermediaries to generate new gene copies. Dinoflagellates are proposed as an ideal model for exploring this process due to the tagging of retrogenes with DNA-encoded remnants of the dinoflagellate-specific splice-leader motif at their 5' end. We conducted a comprehensive search for retrogenes in dinoflagellate transcriptomes to uncover their functional nature and the processes underlying their redundancy. We obtained a high-confidence set of hypothetical functional retrogenes widespread through the dinoflagellate lineage. Through annotations and gene ontology enrichment analysis, we found that the functional diversity of retrogenes reflects the most prevalent and active processes during stress periods, particularly those involving post-translational modifications and cell signalling pathways. Additionally, the significant presence of retrogenes linked to specific biological processes involved in symbiosis and toxin production underscores the role of retrogenes in adaptation. The expression profile and codon composition similar to protein-coding genes confirm the operational status of retrogenes and strengthen the idea that retrogenes recapitulate parental gene expression and function. This study provides new evidence supporting widespread gene retroposition across dinoflagellates and highlights the functional link of retrogenes with the core activity of the cell.

A decade of dinoflagellate genomics illuminating an enigmatic eukaryote cell

Dinoflagellates are a remarkable group of protists, not only for their association with harmful algal blooms and coral reefs but also for their numerous characteristics deviating from the rules of eukaryotic biology. Genome research on dinoflagellates has lagged due to their immense genome sizes in most species (~ 1-250 Gbp). Nevertheless, the last decade marked a fruitful era of dinoflagellate genomics, with 27 genomes sequenced and many insights attained. This review aims to synthesize information from these genomes, along with other omic data, to reflect on where we are now in understanding dinoflagellates and where we are heading in the future. The most notable insights from the decade-long genomics work include: (1) dinoflagellate genomes have been expanded in multiple times independently, probably by a combination of rampant retroposition, accumulation of repetitive DNA, and genome duplication; (2) Symbiodiniacean genomes are highly divergent, but share about 3,445 core unigenes concentrated in 219 KEGG pathways; (3) Most dinoflagellate genes are encoded unidirectionally and are not intron-poor; (4) The dinoflagellate nucleus has undergone extreme evolutionary changes, including complete or nearly complete loss of nucleosome and histone H1, and acquisition of dinoflagellate viral nuclear protein (DVNP); (5) Major basic nuclear protein (MBNP), histone-like protein (HLP), and bacterial HU-like protein (HCc) belong to the same protein family, and MBNP can be the unifying name; (6) Dinoflagellate gene expression is regulated by poorly understood mechanisms, but microRNA and other epigenetic mechanisms are likely important; (7) Over 50% of dinoflagellate genes are "dark" and their functions remain to be deciphered using functional genetics; (8) Initial insights into the genomic basis of parasitism and mutualism have emerged. The review then highlights functionally unique and interesting genes. Future research needs to obtain a finished genome, tackle large genomes, characterize the unknown genes, and develop a quantitative molecular ecological model for addressing ecological questions.

Polyphyletic origin of saxitoxin biosynthesis genes in the marine dinoflagellate Alexandrium revealed by comparative transcriptomics

The marine dinoflagellate Alexandrium is known to form harmful algal blooms, and at least 14 species within the genus can produce saxitoxins (STXs). STX biosynthesis genes (sxt) are individually revealed in toxic dinoflagellates; however, the evolutionary history remains controversial. Herein, we determined the transcriptome sequences of toxic Alexandrium (A. catenella and A. pacificum) and non-toxic Alexandrium (A. fraterculus and A. fragae) and characterized their sxt by focusing on evolutionary events and STX production. Comparative transcriptome analysis revealed higher homology of the sxt in toxic Alexandrium than in non-toxic species. Notably, non-toxic Alexandrium spp. were found to have lost two sxt core genes, namely sxtA4 and sxtG. Expression levels of 28 transcripts related to eight sxt core genes showed that sxtA, sxtG, and sxtI were relatively high (>1.5) in the toxic group compared to the non-toxic group. In contrast, the non-toxic group showed high expression levels in sxtU (1.9) and sxtD (1.7). Phylogenetic tree comparisons revealed distinct evolutionary patterns between 28S rDNA and sxtA, sxtB, sxtI, sxtD, and sxtU. However, similar topology was observed between 28S rDNA, sxtS, and sxtH/T. In the sxtB and sxtI phylogeny trees, toxic Alexandrium and cyanobacteria were clustered together, separating from non-toxic species. These suggest that Alexandrium may acquire sxt genes independently via horizontal gene transfer from toxic cyanobacteria and other multiple sources, demonstrating monocistronic transcripts of sxt in dinoflagellates.

Nuclear genome of dinoflagellates: Size variation and insights into evolutionary mechanisms

Recent progress in high-throughput sequencing technologies has dramatically increased availability of genome data for prokaryotes and eukaryotes. Dinoflagellates have distinct chromosomes and a huge genome size, which make their genomic analysis complicated. Here, we reviewed the nuclear genomes of core dinoflagellates, focusing on the genome and cell size. Till now, the genome sizes of several dinoflagellates (more than 25) have been measured by certain methods (e.g., flow cytometry), showing a range of 3-250 pg of genomic DNA per cell. In contrast to their relatively small cell size, their genomes are huge (about 1-80 times the human haploid genome). In the present study, we collected the genome and cell size data of dinoflagellates and compared their relationships. We found that dinoflagellate genome size exhibits a positive correlation with cell size. On the other hand, we recognized that the genome size is not correlated with phylogenetic relatedness. These may be caused by genome duplication, increased gene copy number, repetitive non-coding DNA, transposon expansion, horizontal gene transfer, organelle-to-nucleus gene transfer, and/or mRNA reintegration into the genome. Ultimate verification of these factors as potential causative mechanisms would require sequencing of more dinoflagellate genomes in the future.

Exchange or Eliminate: The Secrets of Algal-Bacterial Relationships

Algae and bacteria have co-occurred and coevolved in common habitats for hundreds of millions of years, fostering specific associations and interactions such as mutualism or antagonism. These interactions are shaped through exchanges of primary and secondary metabolites provided by one of the partners. Metabolites, such as N-sources or vitamins, can be beneficial to the partner and they may be assimilated through chemotaxis towards the partner producing these metabolites. Other metabolites, especially many natural products synthesized by bacteria, can act as toxins and damage or kill the partner. For instance, the green microalga Chlamydomonas reinhardtii establishes a mutualistic partnership with a Methylobacterium, in stark contrast to its antagonistic relationship with the toxin producing Pseudomonas protegens. In other cases, as with a coccolithophore haptophyte alga and a Phaeobacter bacterium, the same alga and bacterium can even be subject to both processes, depending on the secreted bacterial and algal metabolites. Some bacteria also influence algal morphology by producing specific metabolites and micronutrients, as is observed in some macroalgae. This review focuses on algal-bacterial interactions with micro- and macroalgal models from marine, freshwater, and terrestrial environments and summarizes the advances in the field. It also highlights the effects of temperature on these interactions as it is presently known.

Alexandrium spp.: From Toxicity to Potential Biotechnological Benefits

Many dinoflagellates of the genus Alexandrium are well known for being responsible for harmful algal blooms (HABs), producing potent toxins that cause damages to other marine organisms, aquaculture, fishery, tourism, as well as induce human intoxications and even death after consumption of contaminated shellfish or fish. In this review, we summarize potential bioprospecting associated to the genus Alexandrium, including which Alexandrium spp. produce metabolites with anticancer, antimicrobial, antiviral, as well as anti-Alzheimer applications. When available, we report their mechanisms of action and targets. We also discuss recent progress on the identification of secondary metabolites with biological properties favorable to human health and aquaculture. Altogether, this information highlights the importance of studying which culturing conditions induce the activation of enzymatic pathways responsible for the synthesis of bioactive metabolites. It also suggests considering and comparing clones collected in different locations for toxin monitoring and marine bioprospecting. This review can be of interest not only for the scientific community, but also for the entire population and industries.

Contribution of the satellitome to the exceptionally large genome of dinoflagellates: The case of the harmful alga Alexandrium minutum

Dinoflagellates are known to possess an exceptionally large genome organized in permanently condensed chromosomes. Focusing on the contribution of satellite DNA (satDNA) to the whole DNA content of genomes and its potential role in the architecture of the chromosomes, we present the characterization of the satellitome of Alexandriun minutum strain VGO577. To achieve this, we analyzed Illumina reads using graph-based clustering and performed complementary bioinformatic analyses. In this way, we discovered 180 satDNAs occupying 17.38 % of the genome. The 12 most abundant satDNAs represent the half of the satellitome but no satDNA is overrepresented, with the most abundant contributing ∼1.56 % of the genome. The largest repeat unit is 517 bp long but more than the half of the satDNAs (101) have repeat units shorter than 20 bp. We used FISH to map a selected set of 26 satDNAs. Although some satDNAs generate discrete hybridization signals at specific chromosomal locations (hybridization sites, HS), our cytological analysis showed that most satDNAs are dispersed throughout the genome, probably forming short arrays. Two satDNAs co-localize with the 45S rDNA. With the exception of telomeric DNA, no other satDNA yields HS on all chromosomes. In addition, we analyzed nine satDNAs yielding HS in VGO577 in four other A. minutum strains. Polymorphism at the intraspecific level was found for the presence/absence and/or abundance of some satDNAs, suggesting the amplification/deletion of these satDNAs following geographic separation or during culture maintenance of the strains. We also discuss how these results contribute to the understanding of chromosome architecture and evolution of dinoflagellate genomes.

Bioluminescent Dinoflagellates as a Bioassay for Toxicity Assessment

Dinoflagellates bioluminescence mechanism depends upon a luciferin–luciferase reaction that promotes blue light emission (480 nm) in specialized luminogenic organelles called scintillons. The scintillons contain luciferin, luciferase and, in some cases, a luciferin-binding protein (LBP), which prevents luciferin from non-enzymatic oxidation in vivo. Even though dinoflagellate bioluminescence has been studied since the 1950s, there is still a lack of mechanistic understanding on whether the light emission process involves a peroxidic intermediate or not. Still, bioassays employing luminous dinoflagellates, usually from Gonyaulax or Pyrocystis genus, can be used to assess the toxicity of metals or organic compounds. In these dinoflagellates, the response to toxicity is observed as a change in luminescence, which is linked to cellular respiration. As a result, these changes can be used to calculate a percentage of light inhibition that correlates directly with toxicity. This current approach, which lies in between fast bacterial assays and more complex toxicity tests involving vertebrates and invertebrates, can provide a valuable tool for detecting certain pollutants, e.g., metals, in marine sediment and seawater. Thus, the present review focuses on how the dinoflagellates bioluminescence can be applied to evaluate the risks caused by contaminants in the marine environment.

Metagenomic Analysis of the Species Composition and Seasonal Distribution of Marine Dinoflagellate Communities in Four Korean Coastal Regions

Biomonitoring of dinoflagellate communities in marine ecosystems is essential for efficient water quality management and limiting ecosystem disturbances. Current identification and monitoring of toxic dinoflagellates, which cause harmful algal blooms, primarily involves light or scanning electron microscopy; however, these techniques are limited in their ability to monitor dinoflagellates and plankton, leaving an incomplete analysis. In this study, we analyzed the species composition and seasonal distribution of the dinoflagellate communities in four Korean coastal regions using 18S rRNA amplicon sequencing. The results showed significantly high diversity in the dinoflagellate communities in all regions and seasons. Furthermore, we found seasonally dominant species and causative species of harmful algal blooms (Cochlodinium sp., Alexandrium sp., Dinophysis sp., and Gymnodinium sp.). Moreover, dominant species were classified by region and season according to the difference in geographical and environmental parameters. The molecular analysis of the dinoflagellate community based on metagenomics revealed more diverse species compositions that could not be identified by microscopy and revealed potentially harmful or recently introduced dinoflagellate species. In conclusion, metagenomic analysis of dinoflagellate communities was more precise and obtained results faster than microscopic analysis, and could improve the existing monitoring techniques for community analysis.

Hypothesis: Trans-splicing Generates Evolutionary Novelty in the Photosynthetic Amoeba Paulinella

Plastid primary endosymbiosis has occurred twice, once in the Archaeplastida ancestor and once in the Paulinella (Rhizaria) lineage. Both events precipitated massive evolutionary changes, including the recruitment and activation of genes that are horizontally acquired (HGT) and the redeployment of existing genes and pathways in novel contexts. Here we address the latter aspect in Paulinella micropora KR01 (hereafter, KR01) that has independently evolved spliced leader (SL) trans-splicing (SLTS) of nuclear-derived transcripts. We investigated the role of this process in gene regulation, novel gene origination, and endosymbiont integration. Our analysis shows that 20% of KR01 genes give rise to transcripts with at least one (but in some cases, multiple) sites of SL addition. This process, which often occurs at canonical cis-splicing acceptor sites (internal introns), results in shorter transcripts that may produce 5'-truncated proteins with novel functions. SL-truncated transcripts fall into four categories that may show: (i) altered protein localization, (ii) altered protein function, structure, or regulation, (iii) loss of valid alternative start codons, preventing translation, or (iv) multiple SL addition sites at the 5'-terminus. The SL RNA genes required for SLTS are putatively absent in the heterotrophic sister lineage of photosynthetic Paulinella species. Moreover, a high proportion of transcripts derived from genes of endosymbiotic gene transfer (EGT) and HGT origin contain SL sequences. We hypothesize that truncation of transcripts by SL addition may facilitate the generation and expression of novel gene variants and that SLTS may have enhanced the activation and fixation of foreign genes in the host genome of the photosynthetic lineages, playing a key role in primary endosymbiont integration.

First record of the spatial organization of the nucleosome-less chromatin of dinoflagellates: The nonrandom distribution of microsatellites and bipolar arrangement of telomeres in the nucleus of Gambierdiscus australes (Dinophyceae)

Dinoflagellates are a group of protists whose exceptionally large genome is organized in permanently condensed nucleosome-less chromosomes. In this study, we examined the potential role of repetitive DNAs in both the structure of dinoflagellate chromosomes and the architecture of the dinoflagellate nucleus. Non-denaturing fluorescent in situ hybridization (ND-FSH) was used to determine the abundance and physical distribution of telomeric DNA and 16 microsatellites (1- to 4-bp repeats) in the nucleus of Gambierdiscus australes. The results showed an increased relative abundance of the different microsatellite motifs with increasing GC content. Two ND-FISH probes, (A)₂₀ and (AAT)₅ , did not yield signals whereas the remainder revealed a dispersed but nonrandom distribution of the microsatellites, mostly in clusters. The bean-shaped interphase nucleus of G. australes contained a region with a high density of trinucleotides. This nuclear compartment was located between the nucleolar organizer region (NOR), located on the concave side of the nucleus, and the convex side. Telomeric DNA was grouped in multiple foci and distributed in two polarized compartments: one associated with the NOR and the other peripherally located along the convex side of the nucleus. Changes in the position of the telomeres during cell division evidenced their dynamic distribution and thus that of the chromosomes during dinomitosis. These insights into the spatial organization of microsatellites and telomeres and thus into the nuclear architecture of G. australes will open up new lines of research into the structure and function of the nucleosome-less chromatin of dinoflagellates.

Environmental estrogens in surface water and their interaction with microalgae: A review

Environmental estrogens (EEs) have received extensive attention because they interfere with biological endocrine and reproduction systems by mimicking, antagonizing, or otherwise affecting the actions of endogenous hormones. Additionally, harmful algal blooms have become a global problem in surface water. Microalgae, as an essential primary producer, is especially important for aquatic life and the entire ecosystem. The presence of EEs in surface water may be a potential promoting factor for algal blooms, and microalgae may have effects on the degradation of EEs. This review focuses on the distribution and pollution characteristics of EEs in global surface waters, effects of single and mixed EEs on microalgae regarding growth and toxin production, mechanisms of EEs on microalgae at the cellular and molecular level. The impacts of microalgae on EEs were also discussed. This review provides a risk assessment of EEs and identifies essential clues that will aid in formulating and revising the relevant standards of surface water regarding EEs, which is significant for ecosystems and human health.

Seasonal genotype dynamics of a marine dinoflagellate: Pelagic populations are homogeneous and as diverse as benthic seed banks

Genetic diversity is the basis for evolutionary adaptation and selection under changing environmental conditions. Phytoplankton populations are genotypically diverse, can become genetically differentiated within small spatiotemporal scales and many species form resting stages. Resting stage accumulations in sediments (seed banks) are expected to serve as reservoirs for genetic information, but so far their role in maintaining phytoplankton diversity and in evolution has remained unclear. In this study we used the toxic dinoflagellate Alexandrium ostenfeldii (Dinophyceae) as a model organism to investigate if (i) the benthic seed bank is more diverse than the pelagic population and (ii) the pelagic population is seasonally differentiated. Resting stages (benthic) and plankton (pelagic) samples were collected at a coastal bloom site in the Baltic Sea, followed by cell isolation and genotyping using microsatellite markers (MS) and restriction site associated DNA sequencing (RAD). High clonal diversity (98%-100%) combined with intermediate to low gene diversity (0.58-0.03, depending on the marker) was found. Surprisingly, the benthic and pelagic fractions of the population were equally diverse, and the pelagic fraction was temporally homogeneous, despite seasonal fluctuation of environmental selection pressures. The results of this study suggest that continuous benthic-pelagic coupling, combined with frequent sexual reproduction, as indicated by persistent linkage equilibrium, prevent the dominance of single clonal lineages in a dynamic environment. Both processes harmonize the pelagic with the benthic population and thus prevent seasonal population differentiation. At the same time, frequent sexual reproduction and benthic-pelagic coupling maintain high clonal diversity in both habitats.

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.

Comparison of 15 dinoflagellate genomes reveals extensive sequence and structural divergence in family Symbiodiniaceae and genus Symbiodinium

Background

Dinoflagellates in the family Symbiodiniaceae are important photosynthetic symbionts in cnidarians (such as corals) and other coral reef organisms. Breakdown of the coral-dinoflagellate symbiosis due to environmental stress (i.e. coral bleaching) can lead to coral death and the potential collapse of reef ecosystems. However, evolution of Symbiodiniaceae genomes, and its implications for the coral, is little understood. Genome sequences of Symbiodiniaceae remain scarce due in part to their large genome sizes (1–5 Gbp) and idiosyncratic genome features.

Results

Here, we present de novo genome assemblies of seven members of the genus Symbiodinium, of which two are free-living, one is an opportunistic symbiont, and the remainder are mutualistic symbionts. Integrating other available data, we compare 15 dinoflagellate genomes revealing high sequence and structural divergence. Divergence among some Symbiodinium isolates is comparable to that among distinct genera of Symbiodiniaceae. We also recovered hundreds of gene families specific to each lineage, many of which encode unknown functions. An in-depth comparison between the genomes of the symbiotic Symbiodinium tridacnidorum (isolated from a coral) and the free-living Symbiodinium natans reveals a greater prevalence of transposable elements, genetic duplication, structural rearrangements, and pseudogenisation in the symbiotic species.

Conclusions

Our results underscore the potential impact of lifestyle on lineage-specific gene-function innovation, genome divergence, and the diversification of Symbiodinium and Symbiodiniaceae. The divergent features we report, and their putative causes, may also apply to other microbial eukaryotes that have undergone symbiotic phases in their evolutionary history.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-00994-6.

From the sxtA4 Gene to Saxitoxin Production: What Controls the Variability Among Alexandrium minutum and Alexandrium pacificum Strains?

Paralytic shellfish poisoning (PSP) is a human foodborne syndrome caused by the consumption of shellfish that accumulate paralytic shellfish toxins (PSTs, saxitoxin group). In PST-producing dinoflagellates such as Alexandrium spp., toxin synthesis is encoded in the nuclear genome via a gene cluster (sxt). Toxin production is supposedly associated with the presence of a 4th domain in the sxtA gene (sxtA4), one of the core genes of the PST gene cluster. It is postulated that gene expression in dinoflagellates is partially constitutive, with both transcriptional and post-transcriptional processes potentially co-occurring. Therefore, gene structure and expression mode are two important features to explore in order to fully understand toxin production processes in dinoflagellates. In this study, we determined the intracellular toxin contents of twenty European Alexandrium minutum and Alexandrium pacificum strains that we compared with their genome size and sxtA4 gene copy numbers. We observed a significant correlation between the sxtA4 gene copy number and toxin content, as well as a moderate positive correlation between the sxtA4 gene copy number and genome size. The 18 toxic strains had several sxtA4 gene copies (9-187), whereas only one copy was found in the two observed non-toxin producing strains. Exploration of allelic frequencies and expression of sxtA4 mRNA in 11 A. minutum strains showed both a differential expression and specific allelic forms in the non-toxic strains compared with the toxic ones. Also, the toxic strains exhibited a polymorphic sxtA4 mRNA sequence between strains and between gene copies within strains. Finally, our study supported the hypothesis of a genetic determinism of toxin synthesis (i.e., the existence of several genetic isoforms of the sxtA4 gene and their copy numbers), and was also consistent with the hypothesis that constitutive gene expression and moderation by transcriptional and post-transcriptional regulation mechanisms are the cause of the observed variability in the production of toxins by A. minutum.

The 5S rRNA genes in Alexandrium: their use as a FISH chromosomal marker in studies of the diversity, cell cycle and sexuality of dinoflagellates

Chromosomal markers of the diversity and evolution of dinoflagellates are scarce because the genomes of these organisms are unique among eukaryotes in terms of their base composition and chromosomal structure. Similarly, a lack of appropriate tools has hindered studies of the chromosomal localization of 5S ribosomal DNA (rDNA) in the nucleosome-less chromosomes of dinoflagellates. In this study, we isolated and cloned 5S rDNA sequences from various toxin-producing species of the genus Alexandrium and developed a fluorescence in situ hybridization (FISH) probe that allows their chromosomal localization. Our results can be summarized as follows: 1) The 5S rDNA unit is composed of a highly conserved 122-bp coding region and an intergenic spacer (IGS), the length and sequence of which are variable even within strains. 2) Three different IGS types, one containing the U6 small nuclear RNA (snRNA) gene, were found among four of the studied species (A. minutum, A. tamarense, A. catenella and A. pacificum). 3) In all strains investigated by FISH (A. minutum, A. tamarense, A. pacificum, A. catenella, A. andersonii and A. ostenfeldii), 5S rDNA gene arrays were separate from the nucleolar organizer region, which contains the genes for the large 45S pre-ribosomal RNA. 4) One to three 5S rDNA sites per haploid genome were detected, depending on the strains/species. Intraspecific variability in the number of 5S rDNA sites was determined among strains of A. minutum and A. pacificum. 5) 5S rDNA is a useful chromosomal marker of mitosis progression and can be employed to differentiate vegetative (haploid) vs. planozygotes (diploid) cells. Thus, the FISH probe (oligo-Dino5Smix5) developed in this study facilitates analyses of the diversity, cell cycle and life stages of the genus Alexandrium.

Genomes of the dinoflagellate Polarella glacialis encode tandemly repeated single-exon genes with adaptive functions

BACKGROUND

Dinoflagellates are taxonomically diverse and ecologically important phytoplankton that are ubiquitously present in marine and freshwater environments. Mostly photosynthetic, dinoflagellates provide the basis of aquatic primary production; most taxa are free-living, while some can form symbiotic and parasitic associations with other organisms. However, knowledge of the molecular mechanisms that underpin the adaptation of these organisms to diverse ecological niches is limited by the scarce availability of genomic data, partly due to their large genome sizes estimated up to 250 Gbp. Currently available dinoflagellate genome data are restricted to Symbiodiniaceae (particularly symbionts of reef-building corals) and parasitic lineages, from taxa that have smaller genome size ranges, while genomic information from more diverse free-living species is still lacking.

RESULTS

Here, we present two draft diploid genome assemblies of the free-living dinoflagellate Polarella glacialis, isolated from the Arctic and Antarctica. We found that about 68% of the genomes are composed of repetitive sequence, with long terminal repeats likely contributing to intra-species structural divergence and distinct genome sizes (3.0 and 2.7 Gbp). For each genome, guided using full-length transcriptome data, we predicted > 50,000 high-quality protein-coding genes, of which ~40% are in unidirectional gene clusters and ~25% comprise single exons. Multi-genome comparison unveiled genes specific to P. glacialis and a common, putatively bacterial origin of ice-binding domains in cold-adapted dinoflagellates.

CONCLUSIONS

Our results elucidate how selection acts within the context of a complex genome structure to facilitate local adaptation. Because most dinoflagellate genes are constitutively expressed, Polarella glacialis has enhanced transcriptional responses via unidirectional, tandem duplication of single-exon genes that encode functions critical to survival in cold, low-light polar environments. These genomes provide a foundational reference for future research on dinoflagellate evolution.

Transcriptomic analysis of polyketide synthases in a highly ciguatoxic dinoflagellate, Gambierdiscus polynesiensis and low toxicity Gambierdiscus pacificus, from French Polynesia

Marine dinoflagellates produce a diversity of polyketide toxins that are accumulated in marine food webs and are responsible for a variety of seafood poisonings. Reef-associated dinoflagellates of the genus Gambierdiscus produce toxins responsible for ciguatera poisoning (CP), which causes over 50,000 cases of illness annually worldwide. The biosynthetic machinery for dinoflagellate polyketides remains poorly understood. Recent transcriptomic and genomic sequencing projects have revealed the presence of Type I modular polyketide synthases in dinoflagellates, as well as a plethora of single domain transcripts with Type I sequence homology. The current transcriptome analysis compares polyketide synthase (PKS) gene transcripts expressed in two species of Gambierdiscus from French Polynesia: a highly toxic ciguatoxin producer, G. polynesiensis, versus a non-ciguatoxic species G. pacificus, each assembled from approximately 180 million Illumina 125 nt reads using Trinity, and compares their PKS content with previously published data from other Gambierdiscus species and more distantly related dinoflagellates. Both modular and single-domain PKS transcripts were present. Single domain β-ketoacyl synthase (KS) transcripts were highly amplified in both species (98 in G. polynesiensis, 99 in G. pacificus), with smaller numbers of standalone acyl transferase (AT), ketoacyl reductase (KR), dehydratase (DH), enoyl reductase (ER), and thioesterase (TE) domains. G. polynesiensis expressed both a larger number of multidomain PKSs, and larger numbers of modules per transcript, than the non-ciguatoxic G. pacificus. The largest PKS transcript in G. polynesiensis encoded a 10,516 aa, 7 module protein, predicted to synthesize part of the polyether backbone. Transcripts and gene models representing portions of this PKS are present in other species, suggesting that its function may be performed in those species by multiple interacting proteins. This study contributes to the building consensus that dinoflagellates utilize a combination of Type I modular and single domain PKS proteins, in an as yet undefined manner, to synthesize polyketides.

New Perspectives Related to the Bioluminescent System in Dinoflagellates: Pyrocystis lunula, a Case Study

Pyrocystis lunula is considered a model organism due to its bioluminescence capacity linked to circadian rhythms. The mechanisms underlying the bioluminescent phenomenon have been well characterized in dinoflagellates; however, there are still some aspects that remain an enigma. Such is the case of the presence and diversity of the luciferin-binding protein (LBP), as well as the synthesis process of luciferin. Here we carry out a review of the literature in relation to the molecular players responsible for bioluminescence in dinoflagellates, with particular interest in P. lunula. We also carried out a phylogenetic analysis of the conservation of protein sequence, structure and evolutionary pattern of these key players. The basic structure of the luciferase (LCF) is quite conserved among the sequences reported to date for dinoflagellate species, but not in the case of the LBP, which has proven to be more variable in terms of sequence and structure. In the case of luciferin, its synthesis has been shown to be complex process with more than one metabolic pathway involved. The glutathione S-transferase (GST) and the P630 or blue compound, seem to be involved in this process. In the same way, various hypotheses regarding the role of bioluminescence in dinoflagellates are exposed.

Identification of Prostaglandin Pathway in Dinoflagellates by Transcriptome Data Mining

Dinoflagellates, a major class of marine eukaryote microalgae composing the phytoplankton, are widely recognised as producers of a large variety of toxic molecules, particularly neurotoxins, which can also act as potent bioactive pharmacological mediators. In addition, similarly to other microalgae, they are also good producers of polyunsaturated fatty acids (PUFAs), important precursors of key molecules involved in cell physiology. Among PUFA derivatives are the prostaglandins (Pgs), important physiological mediators in several physiological and pathological processes in humans, also used as “biological” drugs. Their synthesis is very expensive because of the elevated number of reaction steps required, thus the search for new Pgs production methods is of great relevance. One possibility is their extraction from microorganisms (e.g., diatoms), which have been proved to produce the same Pgs as humans. In the present study, we took advantage of the available transcriptomes for dinoflagellates in the iMicrobe database to search for the Pgs biosynthetic pathway using a bioinformatic approach. Here we show that dinoflagellates express nine Pg-metabolism related enzymes involved in both Pgs synthesis and reduction. Not all of the enzymes were expressed simultaneously in all the species analysed and their expression was influenced by culturing conditions, especially salinity of the growth medium. These results confirm the existence of a biosynthetic pathway for these important molecules in unicellular microalgae other than diatoms, suggesting a broad diffusion and conservation of the Pgs pathway, which further strengthen their importance in living organisms.

Biosynthesis of Saxitoxin in Marine Dinoflagellates: An Omics Perspective

Saxitoxin is an alkaloid neurotoxin originally isolated from the clam Saxidomus giganteus in 1957. This group of neurotoxins is produced by several species of freshwater cyanobacteria and marine dinoflagellates. The saxitoxin biosynthesis pathway was described for the first time in the 1980s and, since then, it was studied in more than seven cyanobacterial genera, comprising 26 genes that form a cluster ranging from 25.7 kb to 35 kb in sequence length. Due to the complexity of the genomic landscape, saxitoxin biosynthesis in dinoflagellates remains unknown. In order to reveal and understand the dynamics of the activity in such impressive unicellular organisms with a complex genome, a strategy that can carefully engage them in a systems view is necessary. Advances in omics technology (the collective tools of biological sciences) facilitated high-throughput studies of the genome, transcriptome, proteome, and metabolome of dinoflagellates. The omics approach was utilized to address saxitoxin-producing dinoflagellates in response to environmental stresses to improve understanding of dinoflagellates gene–environment interactions. Therefore, in this review, the progress in understanding dinoflagellate saxitoxin biosynthesis using an omics approach is emphasized. Further potential applications of metabolomics and genomics to unravel novel insights into saxitoxin biosynthesis in dinoflagellates are also reviewed.

Nuclear Gene Transformation in the Dinoflagellate Oxyrrhis marina

The lack of a robust gene transformation tool that allows proper expression of foreign genes and functional testing for the vast number of nuclear genes in dinoflagellates has greatly hampered our understanding of the fundamental biology in this ecologically important and evolutionarily unique lineage of microeukaryotes. Here, we report the development of a dinoflagellate expression vector containing various DNA elements from phylogenetically separate dinoflagellate lineages, an electroporation protocol, and successful expression of introduced genes in an early branching dinoflagellate, Oxyrrhis marina. This protocol, involving the use of Lonza's Nucleofector and a codon-optimized antibiotic resistance gene, has been successfully used to produce consistent results in several independent experiments for O. marina. It is anticipated that this protocol will be adaptable for other dinoflagellates and will allow characterization of many novel dinoflagellate genes.

Progress and promise of omics for predicting the impacts of climate change on harmful algal blooms

Climate change is predicted to increase the severity and prevalence of harmful algal blooms (HABs). In the past twenty years, omics techniques such as genomics, transcriptomics, proteomics and metabolomics have transformed that data landscape of many fields including the study of HABs. Advances in technology have facilitated the creation of many publicly available omics datasets that are complementary and shed new light on the mechanisms of HAB formation and toxin production. Genomics have been used to reveal differences in toxicity and nutritional requirements, while transcriptomics and proteomics have been used to explore HAB species responses to environmental stressors, and metabolomics can reveal mechanisms of allelopathy and toxicity. In this review, we explore how omics data may be leveraged to improve predictions of how climate change will impact HAB dynamics. We also highlight important gaps in our knowledge of HAB prediction, which include swimming behaviors, microbial interactions and evolution that can be addressed by future studies with omics tools. Lastly, we discuss approaches to incorporate current omics datasets into predictive numerical models that may enhance HAB prediction in a changing world. With the ever-increasing omics databases, leveraging these data for understanding climate-driven HAB dynamics will be increasingly powerful.

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.

The Genetic Basis of Toxin Biosynthesis in Dinoflagellates

In marine ecosystems, dinoflagellates can become highly abundant and even dominant at times, despite their comparatively slow growth rates. One factor that may play a role in their ecological success is the production of complex secondary metabolite compounds that can have anti-predator, allelopathic, or other toxic effects on marine organisms, and also cause seafood poisoning in humans. Our knowledge about the genes involved in toxin biosynthesis in dinoflagellates is currently limited due to the complex genomic features of these organisms. Most recently, the sequencing of dinoflagellate transcriptomes has provided us with valuable insights into the biosynthesis of polyketide and alkaloid-based toxin molecules in dinoflagellate species. This review synthesizes the recent progress that has been made in understanding the evolution, biosynthetic pathways, and gene regulation in dinoflagellates with the aid of transcriptomic and other molecular genetic tools, and provides a pathway for future studies of dinoflagellates in this exciting omics era.

Fugacium Spliced Leader Genes Identified from Stranded RNA-Seq Datasets

Trans-splicing mechanisms have been documented in many lineages that are widely distributed phylogenetically, including dinoflagellates. The spliced leader (SL) sequence itself is conserved in dinoflagellates, although its gene sequences and arrangements have diversified within or across different species. In this study, we present 18 Fugacium kawagutii SL genes identified from stranded RNA-seq reads. These genes typically have a single SL but can contain several partial SLs with lengths ranging from 103 to 292 bp. Unexpectedly, we find the SL gene transcripts contain sequences upstream of the canonical SL, suggesting that generation of mature transcripts will require additional modifications following trans-splicing. We have also identified 13 SL-like genes whose expression levels and length are comparable to Dino-SL genes. Lastly, introns in these genes were identified and a new site for Sm-protein binding was proposed. Overall, this study provides a strategy for fast identification of SL genes and identifies new sequences of F. kawagutii SL genes to supplement our understanding of trans-splicing.

Chromosomal markers in the genus Karenia: Towards an understanding of the evolution of the chromosomes, life cycle patterns and phylogenetic relationships in dinoflagellates

Dinoflagellates are a group of protists whose genome is unique among eukaryotes in terms of base composition, chromosomal structure and gene expression. Even after decades of research, the structure and behavior of their amazing chromosomes-which without nucleosomes exist in a liquid crystalline state-are still poorly understood. We used flow cytometry and fluorescence in situ hybridization (FISH) to analyze the genome size of three species of the toxic dinoflagellate genus Karenia as well the organization and behavior of the chromosomes in different cell-cycle stages. FISH was also used to study the distribution patterns of ribosomal DNA (45S rDNA), telomeric and microsatellites repeats in order to develop chromosomal markers. The results revealed several novel and important features regarding dinoflagellate chromosomes during mitosis, including their telocentric behavior and radial arrangement along the nuclear envelope. Additionally, using the (AG)10 probe we identified an unusual chromosome in K. selliformis and especially in K. mikimotoi that is characterized by AG repeats along its entire length. This feature was employed to easily differentiate morphologically indistinguishable life-cycle stages. The evolutionary relationship between Karenia species is discussed with respect to differences in both DNA content and the chromosomal distribution patterns of the DNA sequences analyzed.

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.

Dinoflagellates, a Unique Lineage for Retrogene Research

The birth and evolution of retrogenes have played crucial roles in genome evolution. Dinoflagellates represent a unique lineage for retrogene research because the retrogenes can be reliably identified by the presence of a 22 nucleotide splice leader called DinoSL, which is post-transcriptionally added to the 5' terminus of all mRNAs. Compared to studies of retrogenes conducted in other model genomes, dinoflagellate retrogenes can potentially be more comprehensively characterized because intron-containing retrogenes have already been detected. Unfortunately, dinoflagellate retrogene research has long been neglected. Here, we review the work on dinoflagellate retrogenes and show their distinct character. Like the dinoflagellate genome itself, dinoflagellate retrogenes are also characterized by many unusual features, including a high survival rate and large numbers in the genome. These data are critical complements to what we know about retrogenes, and will further frame our understanding of retroposition and its roles in genome evolution, as well as providing new insights into retrogene studies in other genomes.

Chlorophyll a fluorescence and transcriptome reveal the toxicological effects of bisphenol A on an invasive cyanobacterium, Cylindrospermopsis raciborskii

Bisphenol A has attracted worldwide attention due to its harmful effects on humans, animals and plants. In this study, the toxicological effects of BPA on Cylindrospermopsis raciborskii were assessed based on chlorophyll a fluorescence and transcriptome analyses. The results showed that the growth of C. raciborskii was significantly inhibited when BPA exceeded 0.1 mg L^-1. A marked rise of phase J was observed at a concentration greater than 0.1 mg L^-1, while a K phase appeared at 20 mg L^-1. The chlorophyll a fluorescence parameters of RC/CS₀, F₀, φ_P0, φ_E0, and ψ₀, underwent a significant decline under all treatments of BPA, whereas a significant increase in both V_J and M₀ occurred under all concentrations of BPA. Additionally, ABS/RC and DIo/RC markedly increased at 10 mg L^-1 and 20 mg L^-1. The transcriptome analysis revealed that the genes of photosynthesis, including psbA, psbB, psbC, psbD, apcA, apcB, cpcA, and cpcB, as well as those of chlorophyll and carotenoid biosynthesis, namely hemN, acsF, chlL, chlN, chlP, crtB, pds, were all down-regulated. Moreover, BPA also inhibited the oxidative phosphorylation, glycolysis/gluconeogenesis, citrate cycle (TCA cycle), and fatty acid metabolism in C. raciborskii. Taken together, these results suggest BPA can negatively affect the expression of multiple genes and the vital energy metabolism process to arrest the growth and photosynthesis of C. raciborskii.

A transposable element annotation pipeline and expression analysis reveal potentially active elements in the microalga Tisochrysis lutea

BACKGROUND

Transposable elements (TEs) are mobile DNA sequences known as drivers of genome evolution. Their impacts have been widely studied in animals, plants and insects, but little is known about them in microalgae. In a previous study, we compared the genetic polymorphisms between strains of the haptophyte microalga Tisochrysis lutea and suggested the involvement of active autonomous TEs in their genome evolution.

RESULTS

To identify potentially autonomous TEs, we designed a pipeline named PiRATE (Pipeline to Retrieve and Annotate Transposable Elements, download: https://doi.org/10.17882/51795 ), and conducted an accurate TE annotation on a new genome assembly of T. lutea. PiRATE is composed of detection, classification and annotation steps. Its detection step combines multiple, existing analysis packages representing all major approaches for TE detection and its classification step was optimized for microalgal genomes. The efficiency of the detection and classification steps was evaluated with data on the model species Arabidopsis thaliana. PiRATE detected 81% of the TE families of A. thaliana and correctly classified 75% of them. We applied PiRATE to T. lutea genomic data and established that its genome contains 15.89% Class I and 4.95% Class II TEs. In these, 3.79 and 17.05% correspond to potentially autonomous and non-autonomous TEs, respectively. Annotation data was combined with transcriptomic and proteomic data to identify potentially active autonomous TEs. We identified 17 expressed TE families and, among these, a TIR/Mariner and a TIR/hAT family were able to synthesize their transposase. Both these TE families were among the three highest expressed genes in a previous transcriptomic study and are composed of highly similar copies throughout the genome of T. lutea. This sum of evidence reveals that both these TE families could be capable of transposing or triggering the transposition of potential related MITE elements.

CONCLUSION

This manuscript provides an example of a de novo transposable element annotation of a non-model organism characterized by a fragmented genome assembly and belonging to a poorly studied phylum at genomic level. Integration of multi-omics data enabled the discovery of potential mobile TEs and opens the way for new discoveries on the role of these repeated elements in genomic evolution of microalgae.

A transposable element annotation pipeline and expression analysis reveal potentially active elements in the microalga Tisochrysis lutea

Background

Transposable elements (TEs) are mobile DNA sequences known as drivers of genome evolution. Their impacts have been widely studied in animals, plants and insects, but little is known about them in microalgae. In a previous study, we compared the genetic polymorphisms between strains of the haptophyte microalga Tisochrysis lutea and suggested the involvement of active autonomous TEs in their genome evolution.

Results

To identify potentially autonomous TEs, we designed a pipeline named PiRATE (Pipeline to Retrieve and Annotate Transposable Elements, download: 10.17882/51795), and conducted an accurate TE annotation on a new genome assembly of T. lutea. PiRATE is composed of detection, classification and annotation steps. Its detection step combines multiple, existing analysis packages representing all major approaches for TE detection and its classification step was optimized for microalgal genomes. The efficiency of the detection and classification steps was evaluated with data on the model species Arabidopsis thaliana. PiRATE detected 81% of the TE families of A. thaliana and correctly classified 75% of them. We applied PiRATE to T. lutea genomic data and established that its genome contains 15.89% Class I and 4.95% Class II TEs. In these, 3.79 and 17.05% correspond to potentially autonomous and non-autonomous TEs, respectively. Annotation data was combined with transcriptomic and proteomic data to identify potentially active autonomous TEs. We identified 17 expressed TE families and, among these, a TIR/Mariner and a TIR/hAT family were able to synthesize their transposase. Both these TE families were among the three highest expressed genes in a previous transcriptomic study and are composed of highly similar copies throughout the genome of T. lutea. This sum of evidence reveals that both these TE families could be capable of transposing or triggering the transposition of potential related MITE elements.

Conclusion

This manuscript provides an example of a de novo transposable element annotation of a non-model organism characterized by a fragmented genome assembly and belonging to a poorly studied phylum at genomic level. Integration of multi-omics data enabled the discovery of potential mobile TEs and opens the way for new discoveries on the role of these repeated elements in genomic evolution of microalgae.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4763-1) contains supplementary material, which is available to authorized users.

Comparative Genomics Reveals Two Major Bouts of Gene Retroposition Coinciding with Crucial Periods of Symbiodinium Evolution

Gene retroposition is an important mechanism of genome evolution but the role it plays in dinoflagellates, a critical player in marine ecosystems, is not known. Until recently, when the genomes of two coral-symbiotic dinoflagellate genomes, Symbiodinium kawagutii and S. minutum, were released, it has not been possible to systematically study these retrogenes. Here we examine the abundant retrogenes (∼23% of the total genes) in these species. The hallmark of retrogenes in the genome is the presence of DCCGTAGCCATTTTGGCTCAAG, a spliced leader (DinoSL) constitutively trans-spliced to the 5′-end of all nucleus-encoded mRNAs. Although the retrogenes have often lost part of the 22-nt DinoSL, the putative promoter motif from the DinoSL, TTT(T/G), is consistently retained in the upstream region of these genes, providing an explanation for the high survival rate of retrogenes in dinoflagellates. Our analysis of DinoSL sequence divergence revealed two major bursts of retroposition in the evolutionary history of Symbiodinium, occurring at ∼60 and ∼6 Ma. Reconstruction of the evolutionary trajectory of the Symbiodinium genomes mapped these 2 times to the origin and rapid radiation of this dinoflagellate lineage, respectively. GO analysis revealed differential functional enrichment of the retrogenes between the two episodes, with a broad impact on transport in the first bout and more localized influence on symbiosis-related processes such as cell adhesion in the second bout. This study provides the first evidence of large-scale retroposition as a major mechanism of genome evolution for any organism and sheds light on evolution of coral symbiosis.

Analysis of the genomic basis of functional diversity in dinoflagellates using a transcriptome-based sequence similarity network

Dinoflagellates are one of the most abundant and functionally diverse groups of eukaryotes. Despite an overall scarcity of genomic information for dinoflagellates, constantly emerging high-throughput sequencing resources can be used to characterize and compare these organisms. We assembled de novo and processed 46 dinoflagellate transcriptomes and used a sequence similarity network (SSN) to compare the underlying genomic basis of functional features within the group. This approach constitutes the most comprehensive picture to date of the genomic potential of dinoflagellates. A core-predicted proteome composed of 252 connected components (CCs) of putative conserved protein domains (pCDs) was identified. Of these, 206 were novel and 16 lacked any functional annotation in public databases. Integration of functional information in our network analyses allowed investigation of pCDs specifically associated with functional traits. With respect to toxicity, sequences homologous to those of proteins found in species with toxicity potential (e.g., sxtA4 and sxtG) were not specific to known toxin-producing species. Although not fully specific to symbiosis, the most represented functions associated with proteins involved in the symbiotic trait were related to membrane processes and ion transport. Overall, our SSN approach led to identification of 45,207 and 90,794 specific and constitutive pCDs of, respectively, the toxic and symbiotic species represented in our analyses. Of these, 56% and 57%, respectively (i.e., 25,393 and 52,193 pCDs), completely lacked annotation in public databases. This stresses the extent of our lack of knowledge, while emphasizing the potential of SSNs to identify candidate pCDs for further functional genomic characterization.

Biotechnological and Pharmacological Applications of Biotoxins and Other Bioactive Molecules from Dinoflagellates

The long-lasting interest in bioactive molecules (namely toxins) produced by (microalga) dinoflagellates has risen in recent years. Exhibiting wide diversity and complexity, said compounds are well-recognized for their biological features, with great potential for use as pharmaceutical therapies and biological research probes. Unfortunately, provision of those compounds is still far from sufficient, especially in view of an increasing demand for preclinical testing. Despite the difficulties to establish dinoflagellate cultures and obtain reasonable productivities of such compounds, intensive research has permitted a number of advances in the field. This paper accordingly reviews the characteristics of some of the most important biotoxins (and other bioactive substances) produced by dinoflagellates. It also presents and discusses (to some length) the main advances pertaining to dinoflagellate production, from bench to large scale-with an emphasis on material published since the latest review available on the subject. Such advances encompass improvements in nutrient formulation and light supply as major operational conditions; they have permitted adaptation of classical designs, and aided the development of novel configurations for dinoflagellate growth-even though shearing-related issues remain a major challenge.

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.

Both modular and single-domain Type I polyketide synthases are expressed in the brevetoxin-producing dinoflagellate, Karenia brevis (Dinophyceae)

Dinoflagellates are prolific producers of polyketide compounds, many of which are potent toxins with adverse impacts on human and marine animal health. To identify polyketide synthase (PKS) genes in the brevetoxin-producing dinoflagellate, Karenia brevis, we assembled a transcriptome from 595 million Illumina reads, sampled under different growth conditions. The assembly included 125,687 transcripts greater than 300 nt in length, with over half having >100× coverage. We found 121 transcripts encoding Type I ketosynthase (KS) domains, of which 99 encoded single KS domains, while 22 contained multiple KS domains arranged in 1-3 protein modules. Phylogenetic analysis placed all single domain and a majority of multidomain KSs within a monophyletic clade of protist PKSs. In contrast with the highly amplified single-domain KSs, only eight single-domain ketoreductase transcripts were found in the assembly, suggesting that they are more evolutionarily conserved. The multidomain PKSs were dominated by trans-acyltransferase architectures, which were recently shown to be prevalent in other algal protists. Karenia brevis also expressed several hybrid nonribosomal peptide synthetase (NRPS)/PKS sequences, including a burA-like sequence previously reported in a wide variety of dinoflagellates. This contrasts with a similarly deep transcriptome of Gambierdiscus polynesiensis, which lacked NRPS/PKS other than the burA-like transcript, and may reflect the presence of amide-containing polyketides in K. brevis and their absence from G. polynesiensis. In concert with other recent transcriptome analyses, this study provides evidence for both single domain and multidomain PKSs in the synthesis of polyketide compounds in dinoflagellates.

Building a locally diploid genome and transcriptome of the diatom Fragilariopsis cylindrus

The genome of the cold-adapted diatom Fragilariopsis cylindrus is characterized by highly diverged haplotypes that intersperse its homozygous genome. Here, we describe how a combination of PacBio DNA and Illumina RNA sequencing can be used to resolve this complex genomic landscape locally into the highly diverged haplotypes, and how to map various environmentally controlled transcripts onto individual haplotypes. We assembled PacBio sequence data with the FALCON assembler and created a haplotype resolved annotation of the assembly using annotations of a Sanger sequenced F. cylindrus genome. RNA-seq datasets from six different growth conditions were used to resolve allele-specifc gene expression in F. cylindrus. This approach enables to study differential expression of alleles in a complex genomic landscape and provides a useful tool to study how diverged haplotypes in diploid organisms are used for adaptation and evolution to highly variable environments.

Role of Modular Polyketide Synthases in the Production of Polyether Ladder Compounds in Ciguatoxin-Producing Gambierdiscus polynesiensis and G. excentricus (Dinophyceae)

Gambierdiscus, a benthic dinoflagellate, produces ciguatoxins that cause the human illness Ciguatera. Ciguatoxins are polyether ladder compounds that have a polyketide origin, indicating that polyketide synthases (PKS) are involved in their production. We sequenced transcriptomes of Gambierdiscus excentricus and Gambierdiscus polynesiensis and found 264 contigs encoding single domain ketoacyl synthases (KS; G. excentricus: 106, G. polynesiensis: 143) and ketoreductases (KR; G. excentricus: 7, G. polynesiensis: 8) with sequence similarity to type I PKSs, as reported in other dinoflagellates. In addition, 24 contigs (G. excentricus: 3, G. polynesiensis: 21) encoding multiple PKS domains (forming typical type I PKSs modules) were found. The proposed structure produced by one of these megasynthases resembles a partial carbon backbone of a polyether ladder compound. Seventeen contigs encoding single domain KS, KR, s-malonyltransacylase, dehydratase and enoyl reductase with sequence similarity to type II fatty acid synthases (FAS) in plants were found. Type I PKS and type II FAS genes were distinguished based on the arrangement of domains on the contigs and their sequence similarity and phylogenetic clustering with known PKS/FAS genes in other organisms. This differentiation of PKS and FAS pathways in Gambierdiscus is important, as it will facilitate approaches to investigating toxin biosynthesis pathways in dinoflagellates.

Condition-specific RNA editing in the coral symbiont Symbiodinium microadriaticum

RNA editing is a rare post-transcriptional event that provides cells with an additional level of gene expression regulation. It has been implicated in various processes including adaptation, viral defence and RNA interference; however, its potential role as a mechanism in acclimatization has just recently been recognised. Here, we show that RNA editing occurs in 1.6% of all nuclear-encoded genes of Symbiodinium microadriaticum, a dinoflagellate symbiont of reef-building corals. All base-substitution edit types were present, and statistically significant motifs were associated with three edit types. Strikingly, a subset of genes exhibited condition-specific editing patterns in response to different stressors that resulted in significant increases of non-synonymous changes. We posit that this previously unrecognised mechanism extends this organism's capability to respond to stress beyond what is encoded by the genome. This in turn may provide further acclimatization capacity to these organisms, and by extension, their coral hosts.

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.

Transcriptome Analysis of Scrippsiella trochoidea CCMP 3099 Reveals Physiological Changes Related to Nitrate Depletion

Dinoflagellates are a major component of marine phytoplankton and many species are recognized for their ability to produce harmful algal blooms (HABs). Scrippsiella trochoidea is a non-toxic, marine dinoflagellate that can be found in both cold and tropic waters where it is known to produce “red tide” events. Little is known about the genomic makeup of S. trochoidea and a transcriptome study was conducted to shed light on the biochemical and physiological adaptations related to nutrient depletion. Cultures were grown under N and P limiting conditions and transcriptomes were generated via RNAseq technology. De novo assembly reconstructed 107,415 putative transcripts of which only 41% could be annotated. No significant transcriptomic response was observed in response to initial P depletion, however, a strong transcriptional response to N depletion was detected. Among the down-regulated pathways were those for glutamine/glutamate metabolism as well as urea and nitrate/nitrite transporters. Transcripts for ammonia transporters displayed both up- and down-regulation, perhaps related to a shift to higher affinity transporters. Genes for the utilization of DON compounds were up-regulated. These included transcripts for amino acids transporters, polyamine oxidase, and extracellular proteinase and peptidases. N depletion also triggered down regulation of transcripts related to the production of Photosystems I & II and related proteins. These data are consistent with a metabolic strategy that conserves N while maximizing sustained metabolism by emphasizing the relative contribution of organic N sources. Surprisingly, the transcriptome also contained transcripts potentially related to secondary metabolite production, including a homolog to the Short Isoform Saxitoxin gene (sxtA) from Alexandrium fundyense, which was significantly up-regulated under N-depletion. A total of 113 unique hits to Sxt genes, covering 17 of the 34 genes found in C. raciborskii were detected, indicating that S. trochoidea has previously unrecognized potential for the production of secondary metabolites with potential toxicity.

Uniquely designed nuclear structures of lower eukaryotes

The nuclear structures of lower eukaryotes, specifically protists, often vary from those of yeasts and metazoans. Several studies have demonstrated the unique and fascinating features of these nuclear structures, such as a histone-independent condensed chromatin in dinoflagellates and two structurally distinct nuclear pore complexes in ciliates. Despite their unique molecular/structural features, functions required for formation of their cognate molecules/structures are highly conserved. This provides important information about the structure-function relationship of the nuclear structures. In this review, we highlight characteristic nuclear structures found in lower eukaryotes, and discuss their attractiveness as potential biological systems for studying nuclear structures.

Divinyl chlorophyll a in the marine eukaryotic protist Alexandrium ostenfeldii (Dinophyceae)

Here it is reported the first detection of DV-chl a together with the usual chl a in the marine dinoflagellate Alexandrium ostenfeldii from the Baltic Sea. Growth response and photosynthetic parameters were examined at two irradiances (80 and 240 μmol photons m(-2) s(-1)) and temperatures (15 °C and 19 °C) in a divinylic strain (AOTV-OS20) versus a monovinylic one (AOTV-OS16), using in vivo chl a fluorescence kinetics of PSII to characterize photosynthetic parameters by pulse amplitude modulated fluorescence, (14)C assimilation rates and toxin analyses. The divinylic isolate exhibited slower growth and stronger sensitivity to high irradiance than normal chl a strain. DV-chl a : chl a ratios decreased along time (from 11.3 to < 0.5 after 10 months) and to restore them sub-cloning and selection of strains with highest DV-chl a content was required. A mutation and/or epigenetic changes in the expression of divinyl reductase gene/s in A. ostenfeldii may explain this altered pigment composition. Despite quite severe limitations (reduced fitness and gradual loss of DV-chl a content), the DV-chl a-containing line in A. ostenfeldii could provide a model organism in photosynthetic studies related with chl biosynthesis and evolution.

Multifunctional polyketide synthase genes identified by genomic survey of the symbiotic dinoflagellate, Symbiodinium minutum

Background

Dinoflagellates are unicellular marine and freshwater eukaryotes. They possess large nuclear genomes (1.5–245 gigabases) and produce structurally unique and biologically active polyketide secondary metabolites. Although polyketide biosynthesis is well studied in terrestrial and freshwater organisms, only recently have dinoflagellate polyketides been investigated. Transcriptomic analyses have characterized dinoflagellate polyketide synthase genes having single domains. The Genus Symbiodinium, with a comparatively small genome, is a group of major coral symbionts, and the S. minutum nuclear genome has been decoded.

Results

The present survey investigated the assembled S. minutum genome and identified 25 candidate polyketide synthase (PKS) genes that encode proteins with mono- and multifunctional domains. Predicted proteins retain functionally important amino acids in the catalytic ketosynthase (KS) domain. Molecular phylogenetic analyses of KS domains form a clade in which S. minutum domains cluster within the protist Type I PKS clade with those of other dinoflagellates and other eukaryotes. Single-domain PKS genes are likely expanded in dinoflagellate lineage. Two PKS genes of bacterial origin are found in the S. minutum genome. Interestingly, the largest enzyme is likely expressed as a hybrid non-ribosomal peptide synthetase-polyketide synthase (NRPS-PKS) assembly of 10,601 amino acids, containing NRPS and PKS modules and a thioesterase (TE) domain. We also found intron-rich genes with the minimal set of catalytic domains needed to produce polyketides. Ketosynthase (KS), acyltransferase (AT), and acyl carrier protein (ACP) along with other optional domains are present. Mapping of transcripts to the genome with the dinoflagellate-specific spliced leader sequence, supports expression of multifunctional PKS genes. Metabolite profiling of cultured S. minutum confirmed production of zooxanthellamide D, a polyhydroxy amide polyketide and other unknown polyketide secondary metabolites.

Conclusion

This genomic survey demonstrates that S. minutum contains genes with the minimal set of catalytic domains needed to produce polyketides and provides evidence of the modular nature of Type I PKS, unlike monofunctional Type I PKS from other dinoflagellates. In addition, our study suggests that diversification of dinoflagellate PKS genes comprises dinoflagellate-specific PKS genes with single domains, multifunctional PKS genes with KS domains orthologous to those of other protists, and PKS genes of bacterial origin.

Electronic supplementary material

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

The diversity of algal phospholipase D homologs revealed by biocomputational analysis

Phospholipase D (PLD) participates in the formation of phosphatidic acid, a precursor in glycerolipid biosynthesis and a second messenger. PLDs are part of a superfamily of proteins that hydrolyze phosphodiesters and share a catalytic motif, HxKxxxxD, and hence a mechanism of action. Although HKD-PLDs have been thoroughly characterized in plants, animals and bacteria, very little is known about these enzymes in algae. To fill this gap in knowledge, we performed a biocomputational analysis by means of HMMER iterative profiling, using most eukaryotic algae genomes available. Phylogenetic analysis revealed that algae exhibit very few eukaryotic-type PLDs but possess, instead, many bacteria-like PLDs. Among algae eukaryotic-type PLDs, we identified C2-PLDs and PXPH-like PLDs. In addition, the dinoflagellate Alexandrium tamarense features several proteins phylogenetically related to oomycete PLDs. Our phylogenetic analysis also showed that algae bacteria-like PLDs (proteins with putative PLD activity) fall into five clades, three of which are novel lineages in eukaryotes, composed almost entirely of algae. Specifically, Clade II is almost exclusive to diatoms, whereas Clade I and IV are mainly represented by proteins from prasinophytes. The other two clades are composed of mitochondrial PLDs (Clade V or Mito-PLDs), previously found in mammals, and a subfamily of potentially secreted proteins (Clade III or SP-PLDs), which includes a homolog formerly characterized in rice. In addition, our phylogenetic analysis shows that algae have non-PLD members within the bacteria-like HKD superfamily with putative cardiolipin synthase and phosphatidylserine/phosphatidylglycerophosphate synthase activities. Altogether, our results show that eukaryotic algae possess a moderate number of PLDs that belong to very diverse phylogenetic groups.

Paralytic shellfish toxin biosynthesis in cyanobacteria and dinoflagellates: A molecular overview

Paralytic shellfish toxins (PSTs) are a group of water soluble neurotoxic alkaloids produced by two different kingdoms of life, prokaryotic cyanobacteria and eukaryotic dinoflagellates. Owing to the wide distribution of these organisms, these toxic secondary metabolites account for paralytic shellfish poisonings around the world. On the other hand, their specific binding to voltage-gated sodium channels makes these toxins potentially useful in pharmacological and toxicological applications. Much effort has been devoted to the biosynthetic mechanism of PSTs, and gene clusters encoding 26 proteins involved in PST biosynthesis have been unveiled in several cyanobacterial species. Functional analysis of toxin genes indicates that PST biosynthesis in cyanobacteria is a complex process including biosynthesis, regulation, modification and export. However, less is known about the toxin biosynthesis in dinoflagellates owing to our poor understanding of the massive genome and unique chromosomal characteristics [1]. So far, few genes involved in PST biosynthesis have been identified from dinoflagellates. Moreover, the proteins involved in PST production are far from being totally explored. Thus, the origin and evolution of PST biosynthesis in these two kingdoms are still controversial. In this review, we summarize the recent progress on the characterization of genes and proteins involved in PST biosynthesis in cyanobacteria and dinoflagellates, and discuss the standing evolutionary hypotheses concerning the origin of toxin biosynthesis as well as future perspectives in PST biosynthesis.

SCIENTIFIC QUESTION

Paralytic shellfish toxins (PSTs) are a group of potent neurotoxins which specifically block voltage-gated sodium channels in excitable cells and result in paralytic shellfish poisonings (PSPs) around the world. Two different kingdoms of life, cyanobacteria and dinoflagellates are able to produce PSTs. However, in contrast with cyanobacteria, our understanding of PST biosynthesis in dinoflagellates is extremely limited owing to their unique features. The origin and evolution of PST biosynthesis in these two kingdoms are still controversial.

TECHNICAL SIGNIFICANCE

High-throughput omics technologies, such as genomics, transcriptomics and proteomics provide powerful tools for the study of PST biosynthesis in cyanobacteria and dinoflagellates, and have shown their powerful potential with regard to revealing genes and proteins involved in PST biosynthesis in two kingdoms.

SCIENTIFIC SIGNIFICANCE

This review summarizes the recent progress in PST biosynthesis in cyanobacteria and dinoflagellates with focusing on the novel insights from omics technologies, and discusses the evolutionary relationship of toxin biosynthesis genes between these two kingdoms.