Skip to main content

Integrated omics unveil the secondary metabolic landscape of a basal dinoflagellate



Some dinoflagellates cause harmful algal blooms, releasing toxic secondary metabolites, to the detriment of marine ecosystems and human health. Our understanding of dinoflagellate toxin biosynthesis has been hampered by their unusually large genomes. To overcome this challenge, for the first time, we sequenced the genome, microRNAs, and mRNA isoforms of a basal dinoflagellate, Amphidinium gibbosum, and employed an integrated omics approach to understand its secondary metabolite biosynthesis.


We assembled the ~ 6.4-Gb A. gibbosum genome, and by probing decoded dinoflagellate genomes and transcriptomes, we identified the non-ribosomal peptide synthetase adenylation domain as essential for generation of specialized metabolites. Upon starving the cells of phosphate and nitrogen, we observed pronounced shifts in metabolite biosynthesis, suggestive of post-transcriptional regulation by microRNAs. Using Iso-Seq and RNA-seq data, we found that alternative splicing and polycistronic expression generate different transcripts for secondary metabolism.


Our genomic findings suggest intricate integration of various metabolic enzymes that function iteratively to synthesize metabolites, providing mechanistic insights into how dinoflagellates synthesize secondary metabolites, depending upon nutrient availability. This study provides insights into toxin production associated with dinoflagellate blooms. The genome of this basal dinoflagellate provides important clues about dinoflagellate evolution and overcomes the large genome size, which has been a challenge previously.


Phytoplankton communities are essential components of marine ecosystems, and dinoflagellates are of special interest because they exhibit morphological diversity, high species richness, and the capacity to survive in different ecological niches [1]. They are also infamous contributors to harmful algal blooms (HABs), often producing toxins that are deadly to aquatic organisms and humans [2]. Dinoflagellates exhibit many genetic and cellular features that are highly unusual for eukaryotes. The persistent condensed state of dinoflagellate chromosomes and their liquid crystalline organization, loss of nucleosomal chromatin packaging, use of 5-hydroxymethyluracil in nuclear genomic DNA, and huge genomes of some dinoflagellates (≥ 100 Gbp) are anomalous for eukaryotes [3,4,5]. Recently, the critical role of tandem-duplicated, unidirectional, single-exon genes to survive in cold, low-light environments was reported in two draft genomes (~ 2.8 Gb and ~ 3.0 Gb) of the free-living dinoflagellate, Polarella glacialis [6]. Even with ongoing genomic efforts, understanding of dinoflagellate toxin biosynthesis remains elusive due to their unusually large genomes and limited biosynthetic surveys [4,5,6,7,8,9,10].

Toxic compounds associated with HABs have a polyketide backbone, are synthesized by polyketide synthases (PKSs), and can be linked to non-ribosomal peptide synthases (NRPSs), resulting in hybrid molecules [11]. Several evolutionary events have enabled production of novel polyketides and non-ribosomal peptides [12]. To explore molecular mechanisms involved in secondary metabolite biosynthesis, we sequenced the genome of a basal dinoflagellate, Amphidinium gibbosum, belonging to a genus associated with HABs [3, 13,14,15,16]. Amphidinium species (Gymnodiniales: Gymnodiniaceae) possess intricate secondary metabolic pathways that synthesize unique macrolides with unusual, odd-numbered lactone rings, but their biosynthesis has remained unresolved [17,18,19]. Changes in environmental levels of nitrogen and phosphorus heavily influence the production of toxic metabolites during HABs [20,21,22], and an understanding of nutrient dynamics is critical to any attempt to understand molecular mechanisms associated with toxin production.

Biosynthesis of secondary metabolites having diverse structures and biological activities depends on environmental stresses and is sometimes restricted to specialized structures. Regulation of toxin biosynthesis tends to be coordinated principally at the transcriptional level [23]. Transcriptome analysis of toxic dinoflagellates has been performed [24], but the regulatory mechanisms involved in secondary metabolism during nutrient stress have not been fully explored. While individual omics datasets offer overviews of static states of dinoflagellate systems, integrating several kinds of datasets can strengthen inferences and preclude false assumptions. By sequencing the A. gibbosum genome, transcriptome, and microRNAome, we investigated genomic features and post-transcriptional regulation during nutrient stress, to globally comprehend its secondary metabolism. We identified several miRNAs from the assembled genome and their targets in the transcriptome under phosphate and nitrate starvation. Our integrated omics approach reveals the contributions of repetitive elements and introns in this dinoflagellate genome. It also illustrates the effects of alternative splicing and polycistronic expression and suggests possible implications of miRNA-mediated post-transcriptional regulation of secondary metabolism.

Results and discussion

What accounts for the large genome size and genomic features of the basal dinoflagellate, A. gibbosum?

We estimated that the 6.4-Gb A. gibbosum genome (~ 6.4 Gb by flow cytometry and ~ 6.3 Gb by k-mer analysis) encodes 85,139 genes, of which ~ 48% had matches in available databases (Fig. 1a, b; Table 1; and Additional file 1: Supplementary Fig. 1a-e, Additional file 2: Supplementary Table 1). The size difference between the estimated and assembled genomes may be due to the liquid crystalline structures of dinoflagellate chromosomes [3,4,5]. Genomic data showed the utilization of GC and GA (5′ donor splice sites) in addition to GT and clustering of unidirectional genes, consistent with other dinoflagellate genomes [4, 5, 25] (Fig. 1c, d). This genome included ~ 30% repetitive elements composed of simple repeats (1.97%), low complexity repeats (0.39%), satellite repeats (0.02%), LINEs (0.02%), LTR elements (0.03%), DNA elements (0.1%), and unclassified repeats (27.4%) (Additional file 2: Supplementary Tables 3 and 4). The abundance of repetitive elements may drive genome evolution in dinoflagellates, as reported in Symbiodiniaceae and Polarella glacialis genomes (16–68%) [6, 7]. Comparative analysis of intron and exon features of A. gibbosum provides additional insights into expansion of dinoflagellate genomes (Table 1). Intronic length in A. gibbosum genome is ~ 1.7 Gb, so the intronic region accounts for ~ 27% of the genome, whereas in the Symbiodiniaceae and Polarella glacialis genomes, the average total intronic lengths are 411.5 kb and 737.1 kb, respectively. Despite average exon lengths ranging from 99 to 185 bp, A. gibbosum has the lowest dinoflagellate exon density, with 8.1 exons per gene, compared with 11.3–19.6 exons per gene for other species (Table 1). Large introns have several biological implications, including high energy requirements during transcription, delays in protein production, and greater potential for errors in intron splicing [26, 27]. It follows that some advantage must compensate for such long introns.

Fig. 1

Genomic features of the basal dinoflagellate, Amphidinium gibbosum. a Phylogenetic analysis of dinoflagellates using partial LSU rDNA sequences by maximum likelihood, with red dots at nodes indicating bootstrap support ≥ 80%. b Transmission electron microscopy of A. gibbosum with a lower insert showing a detailed region of condensed chromosomes (lower left: ~ 11 chromosomes in nuclei; lower right: a chromosome). c Non-canonical splice sites show the use of GC and GA, in addition to GT, at the 5′ donor splice site in A. gibbosum, a unique feature of dinoflagellates. d Gene orientation changes using a 9-gene sliding window and 9-gene steps confirm the unidirectional alignment of genes in dinoflagellates. e KEGG pathways recovered from A. gibbosum in comparison with other eukaryotes show biosynthesis of secondary metabolites among top 10 hits. Numbers in brackets indicate the number of enzymes recovered from each pathway category

Table 1 Statistics of the A. gibbosum genome assembly and those of some available dinoflagellate genome assemblies

To understand whether A. gibbosum gene models are conserved at the pathway level, predicted genes were mapped to KEGG reference pathways and compared with those of other dinoflagellates and eukaryotes. This resulted in the recovery of 388 KEGG pathways, indicating that the A. gibbosum genome has most of the pathways present in other eukaryotes (Fig. 1e). Pfam analysis showed Leucine-rich repeat (LRR), Ankyrin, Tetratricopeptide (TPR), and Pentatricopeptide repeat (PPR) domains as the most abundant domains in A. gibbosum (Additional file 2: Supplementary Table 2). Compared with eukaryotes, these repeat domain families, which often contribute to duplication events and to protein-protein interactions, are more abundant in dinoflagellates [9, 28].

Diversified roles of NRPS adenylation domains in dinoflagellates

In order to understand evolution and functions of secondary metabolite genes in A. gibbosum, we conducted molecular phylogenetic analyses of the PKS and NRPS gene families. This confirmed the extensive diversification of these enzyme genes, as previously reported (Fig. 2 and Additional file 1: Supplementary Fig. 2) [10]. Detailed analysis of the adenylation (A) domain of NRPS revealed how specialized metabolites arise in dinoflagellates. The NRPS adenylation domain is the first enzyme in the NRPS complex that selectively incorporates amino acids into NRPSs for biosynthesis of peptide-based natural products, as well as hybrid PKS/NRPS metabolites [11]. The adenylation (A) domain can function as a freestanding protein (Additional file 1: Supplementary Fig. 3), a clear deviation from the usual assembly-line enzymology, known in bacterial genomes [29]. We found that freestanding A domains in A. gibbosum utilize cysteine, valine, and phenylalanine as substrates (Fig. 2), instead of glycine, tryptophan, and phenylalanine, the main substrates utilized by the Symbiodiniaceae [10].

Fig. 2

Affinities of adenylation domains from dinoflagellates show the importance of glycine as a substrate for biosynthesis of specialized toxin secondary metabolites. A molecular phylogenetic tree of adenylation domains indicates protein diversification in Symbiodiniaceae and A. gibbosum. Green- and orange-shaded regions indicate adenylation-domain affinities in Symbiodiniaceae and A. gibbosum, respectively. The Symbiodiniaceae can incorporate glycine (green box) during specialized toxin secondary metabolite biosynthesis such zooxanthellatoxin B (ZT-B) and zooxanthellamide D (ZAD-D), whereas A. gibbosum does not utilize glycine, yielding the simple nitrogen-lacking polyketides, amphidinin A and amphidinolide P. A. gibbosum adenylation sequences are denoted in blue. Red dots indicate a posterior probability ≥ 0.75 using Bayesian inference

Glycine is incorporated into complex metabolites in the Symbiodiniaceae by bridging and forming hybrid molecules, such as zooxanthellatoxin B (ZT-B) and zooxanthellamide D (ZAD-D) [30, 31]; however, none of the amphidinolides and related polyketides [17, 32] isolated from A. gibbosum (amphidinin A and amphidinolide P) contains glycine, resulting in smaller, simpler molecules. Marine dinoflagellates synthesize polyketides that are usually polyol in nature [33]. The carbon skeleton of these polyketides is commonly assembled from acetate, with the rare addition of glycine to form hybrid polyketides [34]. Glycine remains the only amino acid substrate reported in metabolites isolated from dinoflagellates [35, 36], and our analysis suggests that the unique substrate affinities of the NRPS adenylation domain contribute to metabolite complexity in dinoflagellates.

Secondary metabolite biosynthesis responses depend on nutrient starvation regimes

Several studies have demonstrated that nitrogen and phosphorus sources and their availabilities impact both biomass and secondary metabolite production in marine organisms [20,21,22]. It remains unclear which nutrient combinations or limitations drive toxin formation, and this motivated us to investigate whether nutrient starvation affects secondary metabolism in A. gibbosum. We performed deep transcriptome sequencing, recovering 422 pathways, with “metabolic pathways” and “biosynthesis of secondary metabolites” accounting for 1187 proteins (Additional file 1: Supplementary Fig. 1f and Additional file 2: Supplementary Table 5). Under nitrogen starvation, only 16 secondary metabolism genes (PKS and NRPS) were differentially expressed (|log2(FC)| > 2, q < 0.05) (Fig. 3a, b). Gene ontology (GO) enrichment showed that nitrogen starvation has significant effects on nitrogen transport and metabolism (AMT, NRT, NIA, and NRT genes were upregulated) and on anion export (Band 3 gene was downregulated) (|log2(FC)| > 1, p < 0.001) (Fig. 3a and Additional file 1: Supplementary Fig. 4a, b). KEGG pathway enrichment confirmed nitrogen metabolism as the most enriched pathway among upregulated genes, while pathways related to bicarbonate release were the most downregulated genes (p < 0.001) (Additional file 2: Supplementary Table 6). Our analysis revealed novel details about gene expression changes under nitrogen starvation [37]. A. gibbosum apparently tunes its carbon level and nitrogen intake during starvation by downregulating the bicarbonate export system (Band 3 gene) (Fig. 3a). Overall, our data indicate that A. gibbosum modulates incorporation and utilization of several forms of dissolved organic and inorganic nitrogen to respond to nitrogen availability.

Fig. 3

Differentially expressed genes (mRNAs and microRNAs) during nitrogen and phosphate starvation in Amphidinium gibbosum. a Schematic cellular overview of the main differentially expressed genes during nitrogen and phosphate starvation. Orange and blue coloring indicate up- and downregulation, respectively. Green ovals represent plastids, and red boxes indicate mitochondria. A detailed description of proteins is given in Additional file 2: Supplementary Table 9. b Expression profile of PKS and NRPS genes (q < 0.05 and |log2(FC)| > 2) under nitrogen and phosphate starvation. Values show fold changes while N1, N2, and N3; P1, P2, and P3; and NC1, NC2, and NC3 denote triplicate nitrogen, phosphate, and control samples, respectively. Details of the genes are provided in Additional file 2: Supplementary Table 10. NRPS and PKS genes are denoted in red and black, respectively, along the y-axis. c The presence of Dicer (DCL), HEN1, and AGO proteins indicates functional RNAi machinery in A. gibbosum, supported by genomic and transcriptomic data. Whether mature miRNAs in A. gibbosum are methylated is unknown (shaded gray). d Enrichment of miRNA targets during nitrogen starvation shows lactate metabolism as an enriched target process. e The miRNA, agi-miR7721-5p, targets pyruvate metabolism under nitrogen starvation, affecting secondary metabolite biosynthesis. Orange coloring indicates upregulation

Under phosphate starvation, however, 108 PKS and NRPS unigenes were differentially expressed at |log2(FC)| > 2 and q < 0.05 (Fig. 3b and Additional file 1: Supplementary Fig. 5a). Gene ontology (GO) enrichment showed that phosphate starvation upregulates small molecule biosynthesis and downregulates anion release (|log2(FC)| > 2, p value < 0.001) (Fig. 3a and Additional file 2: Supplementary Table 6b). KEGG pathway enrichment confirmed that ribosome, metabolic pathways, and biosynthesis of secondary metabolite pathways are the most enriched pathways among upregulated genes (p < 0.001) (Additional file 2: Supplementary Table 6). During phosphate starvation, membrane transporters (STP, ZIP, AMT, NRT, and AAT) involved in uptake of amino acids, ammonium, dissolved organic phosphate (DOP), metal ions, and nitrate were significantly upregulated. Insufficient dissolved inorganic phosphate can be overcome by utilizing DOPs, which are hydrolysed to release phosphate [38]. This suggests that A. gibbosum can utilize various sources of phosphorus while downregulating genes involved in bicarbonate export, similar to the response observed during nitrogen starvation. Key components of the ATP-consuming glycolytic pathway (e.g., glucokinase, glyceraldehyde-3-phosphate dehydrogenase, and pyruvate kinase) and several ribosomal proteins were significantly upregulated since they are involved in ATP-driven protein synthesis to meet cellular demand for metabolism and phosphate uptake. In both starvation treatments, hierarchical clustering of NRPS and PKS gene expression values revealed two main clusters (Fig. 3b), indicative of a set of co-expressed genes needed for secondary metabolite biosynthesis.

Dinoflagellate carbon-fixing potential increased during phosphate starvation, with several key plastid components (Fig. 3a) being upregulated, including phosphate transporters. This increase may be necessary to fuel augmented cellular processes, as observed in the alga, Prymnesium parvum [39]. Dinoflagellate toxin production changes when environmental parameters such as light, temperature, salinity, and nutrient levels shift [40]. The present analysis shows that the PKS and NRPS genes are upregulated when dinoflagellates are subjected to phosphorus starvation (Fig. 3b) and this can be explained evolutionarily, where microalgal growth slows under nutrient limitation, as cells divert carbon resources for defense [41] (Fig. 3a). Consistent with this theory, increased photosynthetic activity observed during phosphorus starvation in A. gibbosum would be a coordinated physiological response to provide energy necessary for secondary metabolite biosynthesis.

Possible regulation of toxin biosynthesis by microRNAs during nutrient starvation

Based on the low expression of PKS and NRPS unigenes under nitrogen starvation (Fig. 3b), we questioned whether post-transcriptional regulation by microRNAs could be involved. We found expected components of RNAi machinery in A. gibbosum consistent with previous reports [7, 42,43,44,45] (Fig. 3c and Additional file 1: Supplementary Fig. 6). Using the sequenced genome and expressed small RNA data, under phosphate starvation, we found that two miRNAs (agi-miR-6874-5p-2 and a new miRNA denoted, aginovel-mir-0021) were differentially expressed (q value < 0.05, log2(FC) > 2). Upregulation of the two miRNAs was > 18× compared to the control, suggesting that they could have significant effects during phosphate starvation. Indeed, under phosphate starvation, the two upregulated miRNAs targeted pathways involved in fructose-mannose metabolism, proteoglycan synthesis and N-glycan biosynthesis (enrichment > 4×, p < 0.01, Fisher’s exact test) (Additional file 2: Supplementary Table 7). Under nitrogen starvation, we found one miRNA (agi-miR7721-5p) that was differentially expressed (q value < 0.05, log2(FC) > 2). Amphidinium gibbosum had 303 potential target genes, and KEGG pathway target enrichment identified pyruvate-lactate metabolism as a major target (38.4× enrichment, p < 0.001, Fisher’s exact test) (Fig. 3d, e, Additional file 1: Supplementary Fig. 7, and Additional file 2: Supplementary Table 7). This would directly affect production of acetyl-CoA, which is synthesized from pyruvate, a key substrate for polyketide biosynthesis [46], thereby regulating secondary metabolism. No significant PKS and NRPS gene upregulation was observed under nitrogen starvation, in which miRNA-mediated post-transcriptional regulation might affect secondary metabolism by targeting pyruvate biosynthesis. miRNA effects on secondary metabolite biosynthesis have been reported in plants [47, 48].

Transcriptome sequencing reveals diversity of PKS transcripts

Alternative splicing (AS) is an important post-transcriptional regulatory mechanism, whereby a single gene can generate multiple mRNAs, increasing their diversity and complexity [49]. We surveyed five major AS types using rMATS [50] and identified 6970 AS events across 5417 genes, with skipped exons (SE) being the most common AS event (77.2%) (Fig. 4a), followed by alternative 3′splice sites (A3SS) and alternative 5′splice sites (A5SS) (6.8% and 11.3%, respectively). In order to determine biological processes of genes associated with alternative splicing, identified by rMATS [50], GO enrichment was performed. This revealed that ion transport, nucleic acid metabolism, and RNA metabolic process are the most enriched terms (Fig. 4b). Subsequently, we assessed whether AS events were associated with PKS genes. AS landscape analysis at the genome-wide level revealed one PKS gene (g70808) that underwent two AS events, A3SS and SE (Fig. 4c, Additional file 1: Supplementary Fig. 8a). With differential exon usage (DEU) analysis, we found 1 exon (E026) that was differentially expressed (q value < 0.05) during nitrogen starvation (Additional file 1: Supplementary Fig. 8b). AS events function in plant growth and stress responses [51]. Proteins resulting from differently spliced isoforms of the same gene can have different subcellular localization and can inhibit formation of alternative homo- and hetero-dimers [52, 53].

Fig. 4

Alternatively spliced isoforms and polycistronic PKS gene expression in Amphidinium gibbosum. a AS events and their frequencies. SE “skipped exon,” RI “retained intron,” MXE “mutually exclusive exon,” and A3SS and A5SS “alternative 3′ and 5′ splice events”. Black boxes indicate constitutively spliced exons while blue boxes represent alternatively spliced exons. b Gene ontology (GO) biological processes showing significant enrichment of all genes undergoing alternative splicing. c Alternative 3′ splice sites (i) and skipped exons (ii) were identified on a ketosynthase gene (g70808) on scaffold 13486. Phosphate and nitrate experiments are shown in red while controls are in orange. Expression is plotted on the y-axis, genomic coordinates on the x-axis, and isoforms are at bottom in black, with exons depicted in black boxes. Read coverage is represented with numbers. d Sashimi plot showing three uni-directionally aligned PKS genes on scaffold1342 (colored in blue) with multiple polycistronic transcripts (red lines) spanning these genes. PKS module organization within genes is based on PFAM annotation. Iso-Seq read coverage is represented by red vertical blocks, and splicing junction support is shown with numbers. Exons are shown in blue blocks, and lines between blue blocks represent introns. KS "ketosynthase," DH "dehydratase," ER "enoylreductase," KR "ketoreductase"  

To understand how splice junctions contribute to multifunctional polyketide synthase (PKS) isoforms, we conducted Pacbio Isoform sequencing and recovered several transcripts that contained all PKS domains except the acyltransferase (AT) domain, suggesting the trans-acting nature of these enzymes (Additional file 1: Supplementary Fig. 8c). AT genes were indeed trans-acting and belong mainly to the family of malonyl-CoA ACP transferase, contributing malonyl-CoA for chain elongation (Additional file 1: Supplementary Fig. 2b). By mapping these isoforms on the Amphidinium genome, we identified PKS polycistronic transcripts span multiple genes (Fig. 4d). Based on the presence of multiple PKS genes in the genome and their predicted signal peptides (Additional file 1: Supplementary Fig. 2), we asked whether these proteins are localized within the cell. Immunolocalization of ketosynthase and ketoreductase proteins showed that they are localized in mitochondria, chloroplasts, and secretory bodies, as previously reported (Additional file 1: Supplementary Fig. 9) [54]. Additionally, we detected PKS proteins in membrane vesicles, suggesting possible new functions, as demonstrated by their facilitation of nucleation in otolith mineralization [55]. Further functional studies of these proteins will be revealing. By combining different sequencing technologies, we detected polycistronic PKS transcripts, as well as AS events in PKS genes, deepening our understanding of dinoflagellate secondary metabolism. Based on long Iso-Seq reads, we investigated whether secondary metabolite biosynthetic genes contain spliced leader (SL) sequences at their 5′ ends. In dinoflagellates, mRNA maturation is thought to require trans-splicing of the SL sequence [56]. We recovered 548 sequences containing the SL and the relict SL signature, but no PKS transcripts contained it. This could be due to transcript degradation or to a lack of SL sequences at 5′ ends of these transcripts.

Iterative secondary metabolite biosynthesis in dinoflagellates

Polyketide biosynthesis resembles that of fatty acids. The chain is initiated with acetyl-CoA, extended in a series of Claisen ester condensation reactions with malonyl-CoA, and terminated when the required length is reached [10]. While amphidinolides are unique in structure and bioactivity, some similarities exist among them [17], suggesting a common biogenic origin. Complete biosynthesis of an amphidinolide would require all genes present in a cluster, representing up to 500 kb of genomic DNA [11, 18]. Our genomic survey of A. gibbosum confirmed that such long clusters of PKS genes are not present. Each ketosynthase enzyme contributes two carbons to a growing polyketide chain, so a 26-membered polyketide would require at least twelve rounds of carbon addition, implying that such a long cluster is not present in A. gibbosum. Thus, secondary metabolite biosynthesis in dinoflagellates can occur in two ways: (1) monofunctional, separate PKS proteins form an enzyme complex and iteratively catalyze addition of substrate, or (2) multifunctional small PKS proteins utilize substrate in many cycles, to yield a product stabilized by repeat domains that assist such protein-protein interactions (Fig. 5) [57,58,59]. Both these strategies resemble the iterative mono- and multifunctional PKSs of bacterial and fungal systems [60, 61], acquired by horizontal gene transfer [10]. Cross talk between these two co-occurring strategies in dinoflagellates could be mediated by the trans-acting acyltransferase (AT) and NRPS domains, considering that sets of secondary metabolic genes tend to be co-expressed during metabolite biosynthesis (Fig. 3b).

Fig. 5

Strategies for secondary metabolism in dinoflagellates based on a genomic survey. Acetyltransferase acts in trans to provide activated substrates to acyl carrier protein (ACP) with extensions and modifications by optional domains, terminating with hydrolysis by thioesterase. The adenylation domain activates the amino acyl substrate and bridges intermediate products, acting as a mediator. KS, ketosynthase; KR, ketoreductase, AT, acetyltransferase; DH, dehydratase, ER, enoylreductase; TE, thioesterase; A, adenylation. ACP are omitted for clarity


In this study, we applied an integrated omics approach to understand dinoflagellate secondary metabolite biosynthesis. To this end, we sequenced the genome of A. gibbosum and identified key features that regulate secondary metabolite levels and structural diversity. We hypothesize that miRNA-mediated, post-transcriptional regulation in A. gibbosum, which targets primary pyruvate metabolism, subsequently affects secondary metabolism. This study represents a first step to illuminate key molecular events involved in dinoflagellate secondary metabolism, and it should facilitate studies of HAB formation and associated toxin production. Ongoing high-throughput sequencing of dinoflagellate genomes promises to be informative, not only for understanding toxin secondary metabolism genes, but also for better insights into their genome organization. The availability of this first basal dinoflagellate genome provides important clues about dinoflagellate evolution and extends the genome size limit that has been a challenge for several years.


Biological sample

Amphidinium gibbosum was isolated from inner cells of a marine acoelomorph, Amphiscolops sp., collected near Ishigaki Island, Japan. The culture was maintained in artificial seawater (ASW) containing 1X Guillard’s (F/2) marine-water enrichment solution and an antibiotic-antimycotic mix in a 25 °C incubator under a 12:12 light and dark cycle. Subculture was performed with fresh medium approximately every 4 weeks and was handled aseptically. For transmission electron microscopy (TEM), cells were fixed in 2.5% glutaraldehyde for 1 h, washed 3× with 0.1 M cacodylate buffer, and incubated in 1% osmium tetroxide for 30 min. Cells were then washed and dehydrated in an ethanol series (70%, 80%,90%, 95%, 100%, 100%, 100%), at 5-min intervals. Samples were infiltrated with ethanol-Epon resin for 30 min and steeped in 100% resin overnight. The resin was polymerized at 60 °C for 2 days. Sections were cut using a diamond knife and viewed under a JEM-1230R JEOL microscope. The phylogenetic position of A. gibbosum was confirmed by aligning and trimming partial LSU rDNA sequences of several dinoflagellates and performing maximum likelihood analysis using RaxML [62]. Phylogenetic assignment was consistent with the taxonomic description [63].

Genome size estimation

For A. gibbosum genome size estimation, nuclear DNA from three replicates was measured using fluorescence-activated cell sorting (FACS) with Xenopus laevis (n = 3) as an internal control of known genome size. Nuclear extraction and staining were performed using a Partec CyStainPI absolute T kit (Partec #05-5023), following the manufacturer’s protocol, and fluorescence signals were measured with a BD Accuri C6 cell analyzer (BD Bioscience). The reported measurement for A. gibbosum reflects the 1C genome content, as Amphidinium is reportedly haploid in culture. K-mer analysis was performed using Jellyfish (v2.1.3) [64], and resulting histograms were visualized using GenomeScope [65] to survey the genome size and repeat content.

DNA sample preparation and sequencing

Cells were centrifuged at 3000g for 10 min and washed using TEN buffer (100 mM Tris-Cl pH 8, 100 mM EDTA pH 8, 1.5 M NaCl, 0.5 mg/mL proteinase K, and 7% SDS) for 2 h at 65 °C so as to lyse bacterial contaminants. DNA was extracted using a modified protocol [66] of gentle rotation for 1 h after addition of chloroform-isoamyl alcohol (24:1) before ethanol precipitation [4]. Isolated DNA was further cleaned using ethanol precipitation. DNA was fragmented and paired-end libraries with an insert size of 620–820 bp were prepared. Libraries were quantified by qPCR and sequenced using an Illumina Miseq, according to the manufacturer’s protocols. This generated ~ 10 Gb of 2 × 300 bp paired-end data. The same library was further sequenced using a Hiseq 2500, generating ~ 586 Gb of 2 × 125 bp of data. Reads were merged and trimmed using Trimmomatic (v0.35) [67] and were quality-checked using FastQC (v0.11.4) [68]. Additionally, 12 mate-pair libraries were constructed using Nextera technology with 2–18-kb inserts selected using the Bluepippin and SageELF systems. Mate-pair libraries were sequenced with a Hiseq 4000, generating ~ 200 Gb of data. Raw mate-paired reads were filtered using NextClip (v1.31) [69]. Genome assembly employed Platanus (v2.1.4) [70], and the assembled genome was subjected to two rounds of scaffolding with SSPACE (V3.0) [71]. Gaps in scaffolds were filled using GapCloser (v1.12) [72] (Additional file 1: Supplementary Fig. 10A).

Evaluation of genome assembly completeness and removal of contaminating sequences

The scaffolded Amphidinium genome was checked for genome completeness using BUSCO 303 highly conserved eukaryotic genes (CEGs) [73]. Additionally, the BLAST suite was used to recover 458 CEGs from CEGMA [74] against the Amphidinium genome to identify potential homologs at a cutoff value of 1e−5. To identify bacterial and viral contaminants, we conducted a BLASTN search against several databases that we built by retrieving draft and complete bacterial genomes and viral genomes from NCBI and PhanToME. A combination of cutoffs (total bit score > 1000, E ≤ 10−20) was used to identify scaffolds with similarities to bacterial and viral sequences.

cDNA construction, Iso-Seq sequencing, and data processing

RNA was extracted from cells growing under standard conditions (12:12 light and dark cycle), and a cDNA library was constructed using a TruSeq Stranded RNA Sample Prep Kit (Illumina). Libraries were quantified and validated by qPCR and with a 2100 Agilent Bioanalyzer, respectively. The validated library was subsequently sequenced using two lanes of Hiseq 2500 (Illumina). Reads were trimmed using Trimmomatic (v0.35) [67], quality-checked using FastQC (v0.11.4) [68], and assembled de novo using Trinity (v2.3.2) [75]. For Iso-Seq sequencing, RNA was extracted from several culture treatments and pooled. High-quality RNAs (RIN > 7.0) were used for cDNA synthesis using a Clontech SMARTer PCR cDNA kit. Size fractionation (0.7–2.5, 2.5–7, and > 7 kb) was conducted using the SageELF system (Sage Science, Beverly, MA, USA). Libraries were sequenced with the Pacific Biosciences RS II platform (P6-P4 chemistry) and a 360-min movie length. In total, 16 SMRT cells were sequenced. Raw sequencing data were processed using the RS_Iso-Seq protocol. HQ and LQ reads were error-corrected by employing proovread (v2.14) [76] using Illumina RNA-seq data. Reads were then merged, and “cd-hit-est” from CD-HIT (v4.6) [77] was used to remove redundancy with parameters: -c 0.99 -G 0 -aL 0.00 -aS 0.99 -AS 30 -M 0 -d 0 -p 1 -T 24. Non-redundant transcripts were further processed with Cogent ( Polished Iso-Seq sequences were surveyed for the dinoflagellate spliced leader (CCGTAGCCATTTTGGCTCAAG) and the relict dinoSL (CCGTAGCCATTTTGGCTCAAGCCATTTTGGCTCAAG) [78] sequences using BLAST with no gaps and up to 1 mismatch permitted.

Repetitive element annotation and gene model prediction

In order to confirm splice sites, the assembled transcriptome was mapped to the assembled genome using GMAP [79]. For annotating transposable elements (TEs), de novo repeats within the genome were identified using an l-mer size of 17 bp with RepeatScout [80]. A combined library was made, consisting of de novo repeats and eukaryotic TEs from RepBase. This library was then used to locate and annotate repetitive elements in the assembled genome using RepeatMasker [81]. RNA-seq reads were mapped to a soft-masked genome using STAR [82] and the BRAKER2 pipeline [83]. UTR and gene model prediction were performed with Augustus (v3.2.3) [84]. To improve gene prediction accuracy, intron and exon hints were generated as additional evidence of gene structure and location by mapping Illumina and Iso-Seq transcripts to the genome with GMAP [79] and STAR [82]. Hints were then used to perform final gene prediction using a modified version of Augustus (v3.2.3) [84], in which the source code was changed in consideration of non-canonical exon-intron boundaries. The final set of predicted proteins was annotated against UniProt [85] and PFAM [86]. Briefly, BLASTP searches for all protein models were performed with the SwissProt and TrEMBL databases (October 2018 release). Amino acid sequences were subjected to PFAM [86] domain searches using HMMER (v3.1b2) [87], and hits larger than 1e−5 were discarded. For KEGG pathway analysis, the online service on the KEGG Automatic Server (KAAS) was used to assign predicted genes to KEGG orthologs (bi-directional best hit method) and mapped orthologs to KEGG pathways.

Phylogenetic analysis of PKS and NRPS proteins and prediction of substrate specificities

The dataset used previously [10] was repopulated with ketosynthase, acyltransferase, adenylation, and condensation protein sequences from the A. gibbosum genome. Briefly, four datasets were created, consisting of 244 KS sequences (225 aa), 104 AT sequences (208 aa), 121 A-sequences (272 aa), and 111 C-sequences (253 aa). Mono- and multifunctional domain-containing sequences were aligned using MUSCLE [88], and domain areas with best alignment were retained while regions with ambiguity were removed. Two methods for phylogenetic reconstruction were used, maximum likelihood using RaxML [62] (1000 bootstraps and LG + G model) and Bayesian inference (run to a maximum of 6 million generations plus 4 chains, or until probability approached 0.01), using MrBayes (v3.2) [89]. Substrate specificity of A. gibbosum AT sequences was generated using I-TASSER [90]. In order to determine the A-domain specificity and C-domain types, the LSI-based A-domain predictor and NaPDos were used, respectively [91, 92]. The phylogenetic analysis of A-domain and a part of its substrate specificity are depicted in Fig. 2. Sequence alignment of the A-domain is provided as Additional File 3. PKS protein subcellular localization was detected using ChloroP 1.1 and TargetP 1.1 and was further confirmed with DeepLoc [93,94,95].

Nutrient starvation experiment

For a nitrate-starved culture, the culture medium was prepared by supplementing artificial seawater (ASW) with F/2 medium containing a reduced nitrate concentration (150 μM). For a phosphate-starved culture, the phosphate level was 22 μM. A phosphate and nitrate-replete treatment was set up as the control, in which nitrate and phosphate concentrations were 880 and 36 μM, respectively. Both starvation (depleted) and control treatments were conducted in triplicate (n = 3). First, measurements were started after 24 h of stabilization, and this was counted as day 1. Nitrate and phosphate levels were monitored using the Griess and phosphomolybdenum blue spectrophotometric methods, respectively [96, 97], until their concentrations were undetectable. Other physiological parameters, such as cell concentration, chlorophyll a, and photochemical efficiency (Fv/Fm ratio), were also monitored (Additional File 1: Supplementary Fig. 10B). Cell counts were obtained by fixing cells in formalin and using a hemocytometer for visualization. 1-mL samples were centrifuged, and cell pellets were immersed in N,N-dimethylformamide (DMF) and kept at − 20 °C for at least 12 h in order to extract chlorophyll a, which was measured using a Turner Trilogy (Turner Designs fluorometer, USA) and then averaged to content per cell. Photochemical efficiency was monitored with a Xe-PAM (Walz, Germany).

Gene expression analysis during nutrient starvation

When dissolved nitrate and phosphate reached an undetectable level, ~ 107 cells were collected, snap frozen in liquid nitrogen, and ground using a cryopress. RNA was extracted from 3 control, 3 nitrate-starved, and 3 phosphate-starved samples using PureLink reagent. Four micrograms of RNA was used for cDNA library construction with a TruSeq Stranded RNA Sample Prep Kit (Illumina). Libraries were quantified and validated by qPCR and with a 2100 Agilent Bioanalyzer, respectively, and sequenced in two lanes of a Hiseq 4000 (Illumina). Reads were trimmed using Trimmomatic (v0.35) [67], quality-checked using FastQC (v0.11.4) [68], and assembled using Trinity (v2.3.2) [75]. The assembly was processed with CD-HIT-EST (v4.6.7) [77] using a clustering threshold of 0.95. Functional annotation of non-redundant contigs was performed using BLAST with several databases: UniProt, GeneBank non-redundant (nr), Kyoto Encyclopedia of Genes and Genomes (KEGG), and eggNOG (E value cutoff of 10−5) [85, 98]. Transcriptomic gene completeness was evaluated using BUSCO (v3.0.2) [73]. For identification of differentially expressed transcripts, expression abundance was quantified using RSEM [99]. The R package, EdgeR [100], was used to identify differentially expressed genes with adjusted p values (q value) determined with the Benjamini, Krieger, and Yekutieli correction of the PRISM package. Figure 3a, b depicts the results of this analysis. Gene ontology term functional enrichment was performed using Fisher’s exact test in topGo with the parent-child analysis to categorize whether differentially expressed genes were enriched in molecular function, cellular components, and biological processes [101]. KEGG pathway enrichment was performed using DAVID [102] by applying Fisher’s exact test.

Small RNA sequencing for the nutrient starvation experiment

Small RNAs were isolated from the same RNA pellet (n = 3) collected from the depleted-replete experiments using the NEXTflexTM Small RNA-seq Kit V3 (Bioo Scientific). Single-end reads (1 × 50 bp) were generated on a Hiseq 2500 platform. Reads were cleaned by removing adapter and polyA/N sequences using Cutadapt-1.4.1 [103], and reads within the range of 17–25 were retained. Reads were further collapsed using the script of the MiRDeep2 package [104]. Sequences having hits to various non-coding RNAs (rRNAs, tRNAs, snRNAs, snoRNAs, and scRNAs) of the RNAcentral database (The RNAcentral Consortium, 2015) were discarded. Bowtie (v1.1.12) [105] was used to map clean, small RNA reads to the Amphidinium gibbosum genome with no mismatches and 1 alignment. Mapped reads were further queried against known miRNAs in miRBase 22.0 ( miRNAs were annotated using the miRdeep2 package. Previous miRNA criteria [42] were applied to the list of annotated miRNAs. miRNA expression level profiling was conducted and normalized using the script of the miRdeep2 package where processed reads were mapped to identified miRNA precursors. EdgeR [100] was then used to identify differentially expressed miRNAs at FDR < 0.05 (adjusted p value), as determined by Benjamini, Krieger, and Yekutieli of the PRISM package and |log2(FC)| > 1. Only miRNAs present in at least 2 replicates were considered further. For predicting mRNA targets of the miRNAs, 3′UTR sequences of unigenes were used by miRanda [106] under strict criteria. GO and KEGG pathway enrichment was performed for predicted target unigenes of differentially expressed miRNAs using topGO and DAVID, respectively [101, 102]. Figure 3c–e depicts the results of this analysis.

Identification of key proteins in microRNA biogenesis pathways

In order to confirm the presence of a miRNA biogenesis pathway, sequences of three core protein families involved in RNA interference (i.e., Argonaute, Dicer, and HEN1) were retrieved for model organisms (H. sapiens, C. elegans, S. pombe, D. melanogaster, and A. thaliana) from UniProtKB [85]. Sequences were then queried against predicted proteins from the A. gibbosum transcriptome using BLASTP (E value cutoff of 10−10). Hits were then searched for specific domains (a PAZ domain and a pair of RNase III domains for Dicer, Piwi and Dicer domains for Argonaute, and a methyltransferase domain for HEN1) needed for functional activity using InterProScan [107]. Alignment of homologs against retrieved RNAi proteins from model organisms was conducted using Clustal Omega [108] and visualized using Jalview [109].

Alternative splicing (AS) and enrichment analyses

In order to identify alternative splicing events (Skipped exon [SE], alternative 5′ splice site [A5SS], alternative 3′ splice site [A3SS], mutually exclusive exons [MXE], retained intron [RI]), rMATS [50] was used. Briefly, processed RNA-seq reads from nutrient stress experiments were mapped to the genome using STAR [82] and MISO [110] was employed to verify AS events. Iso-Seq reads were also mapped to the genome using STAR [82] to confirm the presence of exons. To evaluate differential exon usage, DEXSeq (version 1.28.3) [111] was used. Exon expression counts for each replicate in nutrient stress experiments were quantified using the Amphidinium genome annotation and BAM files generated from STAR [81] mapping. Default normalization of libraries was performed, and p values were corrected using FDR with a p-adjust cutoff of < 0.05. Gene ontology term functional enrichment of all genes showing alternative splicing was performed using the GOstats R package [112] and visualized using REVIGO [113]. Figure 4 depicts the results of these analyses.

PKS protein immunolocalization

Cells grown in normal ASW were first fixed in 2% paraformaldehyde in seawater, washed three times with PBS, and incubated in 50% methanol:PBS (5 min). Cells were then deposited on poly-l-lysine-coated coverslips, blocked with 5% normal goat serum for 1 h, and incubated with primary anti-PKS antibodies (KS and KR) at 1:100 dilution overnight at 4 °C. Cells were then incubated with Alexa Fluor-488-conjugated secondary antibodies for 1 h at room temperature. Coverslips were then mounted with Vectashield on glass slides and observed under a Zeiss Axio-Observer Z1 LSM 780 microscope. Data were collected using ZEN software (version For negative controls, cells were treated with PBS instead of primary antibodies. Stacks were analyzed using ImageJ [114].

Availability of data and materials

Sequence data from this study are available in the NCBI Short Read Archive (SRA) Bioproject ID PRJNA551917 [115]. Assembled genome, transcriptome, predicted gene models, and proteins are available at: [116].


  1. 1.

    Smayda TJ, Reynolds CS. Strategies of marine dinoflagellate survival and some rules of assembly. J Sea Res. 2013;49:95–106.

    Article  Google Scholar 

  2. 2.

    Wang D-Z. Neurotoxins from marine dinoflagellates: a brief review. Mar Drugs. 2008;6:349–71.

    PubMed  PubMed Central  Article  Google Scholar 

  3. 3.

    Wisecaver JH, Hackett JD. Dinoflagellate genome evolution. Annu Rev Microbiol. 2011;65:369–87.

    CAS  PubMed  Article  Google Scholar 

  4. 4.

    Shoguchi E, Shinzato C, Kawashima T, Gyoja F, Mungpakdee S, Koyanagi R, et al. Draft assembly of the Symbiodinium minutum nuclear genome reveals Dinoflagellate gene structure. Curr Biol. 2013;23:1399–408.

    CAS  PubMed  Article  Google Scholar 

  5. 5.

    Aranda M, Li Y, Liew YJ, Baumgarten S, Simakov O, Wilson MC, et al. Genomes of coral dinoflagellate symbionts highlight evolutionary adaptations conducive to a symbiotic lifestyle. Sci Rep. 2016;6:39734.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  6. 6.

    Stephens TG, González-Pech RA, Cheng Y, Mohamed AR, Burt DW, Bhattacharya D, et al. Genomes of the dinoflagellate Polarella glacialis encode tandemly repeated single-exon genes with adaptive functions. BMC Biol. 2020;18:56.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  7. 7.

    Lin S, Cheng S, Song B, Zhong X, Lin X, Li W, et al. The Symbiodinium kawagutii genome illuminates dinoflagellate gene expression and coral symbiosis. Science. 2015;350:691–4.

    CAS  PubMed  Article  Google Scholar 

  8. 8.

    Liu H, Stephens TG, Gonzalez-Pech RA, Beltran VH, Lapeyre B, Bongaerts P, et al. Symbiodinium genomes reveal adaptive evolution of functions related to coral-dinoflagellate symbiosis. Commun Biol. 2018;1:95.

    PubMed  PubMed Central  Article  Google Scholar 

  9. 9.

    Shoguchi E, Beedessee G, Tada I, Hisata K, Kawashima T, Takeuchi T, et al. Two divergent Symbiodinium genomes reveal conservation of a gene cluster for sunscreen biosynthesis and recently lost genes. BMC Genomics. 2018;19:458.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  10. 10.

    Beedessee G, Hisata K, Roy MC, Van Dolah FM, Satoh N, Shoguchi E. Diversified secondary metabolite biosynthesis gene repertoire revealed in symbiotic dinoflagellates. Sci Rep. 2019;9:1204.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  11. 11.

    Kellmann R, Stüken A, Orr RJS, Svendsen HM, Jakobsen KS. Biosynthesis and molecular genetics of Polyketides in marine dinoflagellates. Mar Drugs. 2010;8:1011–48.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  12. 12.

    Fischbach M, Walsh CT, Clardy J. The evolution of gene collectives: how natural selection drives chemical innovation. Proc Natl Acad Sci U S A. 2008;105:4601–8.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  13. 13.

    Lee JJ, Olea R, Cevasco M, Pochon X, Correia M, Shpigel M, et al. A marine dinoflagellate, Amphidinium eilatiensis n. sp., from the benthos of a Mariculture sedimentation pond in Eilat, Israel. J Eukaryot Microbiol. 2003;50:439–48.

    PubMed  Article  Google Scholar 

  14. 14.

    Baig HS, Saifullah SM, Dar A. Occurrence and toxicity of Amphidinium carterae Hulburt in the north Arabian Sea. Harmful Algae. 2006;5:133–40.

    CAS  Article  Google Scholar 

  15. 15.

    Gárate-Lizárraga I. Proliferation of Amphidinium carterae (Gymnodiniales: Gymnodiniaceae) in Bahía de La Paz, Gulf of California. CICIMAR Oceánides. 2012;27:37–49.

    Google Scholar 

  16. 16.

    Murray SA, Kohli GS, Farrell H, Spiers ZB, Place AR, Doranres-Aranda JJ, et al. A fish kill associated with a bloom of Amphidinium carterae in a coastal lagoon in Sydney, Australia. Harmful Algae. 2015;49:19–28.

    PubMed  PubMed Central  Article  Google Scholar 

  17. 17.

    Kobayashi J, Kubota T. Bioactive macrolides and polyketides from marine dinoflagellates of the genus Amphidinium. J Nat Prod. 2007;70:451–60.

    CAS  PubMed  Article  Google Scholar 

  18. 18.

    Kubota T, Iinuma Y, Kobayashi J. Cloning of polyketide synthase genes from Amphidinolide-producing dinoflagellate Amphidinium sp. Biol Pharm Bull. 2006;29:1314–8.

    CAS  PubMed  Article  Google Scholar 

  19. 19.

    Murray SA, Garby T, Hoppenrath M, Neilan BA. Genetic diversity, morphological uniformity and Polyketide production in Dinoflagellates (Amphidinium, Dinoflagellata). PLoS One. 2012;7:e38253.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  20. 20.

    Wang D, Ho AYT, Hsieh DPH. Production of C2 toxin by Alexandrium tamarense CI01 using different culture methods. J Appl Phycol. 2002;14:461–8.

    CAS  Article  Google Scholar 

  21. 21.

    Erdner DL, Anderson DM. Global transcriptional profiling of the toxic dinoflagellate Alexandrium fundyense using massively parallel signature sequencing. BMC Genomics. 2006;7:88.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  22. 22.

    Falkowski PG, Barber RT, Smetacek V. Production biogeochemical controls and feedbacks on ocean primary biogeochemical controls and feedbacks on ocean primary production. Science. 1998;281:200–7.

    CAS  PubMed  Article  Google Scholar 

  23. 23.

    Colinas M, Goossens A. Combinatorial transcriptional control of plant specialized metabolism. Trends Plant Sci. 2018;23:324–36.

    CAS  PubMed  Article  Google Scholar 

  24. 24.

    Moustafa A, Evans AN, Kulis DM, Hackett JD, Erdner DL, Anderson DM, Bhattacharya D. Transcriptome profiling of a toxic dinoflagellate reveals a gene-rich protist and a potential impact on gene expression due to bacterial presence. PLoS One. 2010;5:e9688.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  25. 25.

    Bachvaroff TR, Place AR. From stop to start: tandem gene arrangement, copy number and trans-splicing sites in the dinoflagellate Amphidinium carterae. PLoS One. 2008;3:e2929.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  26. 26.

    Fedorova L, Fedorov A. Puzzles of the human genome: why do we need our introns? Current Genomics. 2005;6:589–95.

    CAS  Article  Google Scholar 

  27. 27.

    Sun H, Chasin LA. Multiple splicing defects in an intronic false exon. Mol Cell Biol. 2000;20:6414–25.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  28. 28.

    Schaper E, Anisimova M. The evolution and function of protein tandem repeats in plants. New Phytol. 2015;206:397–410.

    CAS  PubMed  Article  Google Scholar 

  29. 29.

    Lin S, Lanen SGV, Shen B. A free-standing condensation enzyme catalyzing ester bond formation in C-1027 biosynthesis. Proc Natl Acad Sci U S A. 2009;106:4183–8.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  30. 30.

    Nakamura H, Asari T, Fujimaki K, Maruyama K, Murai A, Ohizumi Y, Kan Y. Zooxanthellatoxin-B, vasoconstrictive congener of zooxanthellatoxin-a from a symbiotic dinoflagellate Symbiodinium sp. Tetrahedron Lett. 1995;36:7255–8.

    CAS  Article  Google Scholar 

  31. 31.

    Fukatsu T, Onodera K, Ohta Y, Oba Y, Nakamura H, Shintani T, et al. Zooxanthellamide D, a polyhydroxy polyene amide from a marine dinoflagellate, and chemotaxonomic perspective of the symbiodinium polyols. J Nat Prod. 2007;70:407–11.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  32. 32.

    Kubota T, Sato H, Iwai T, Kobayashi J. Biosynthetic study of Amphidinin a and Amphidinolide P. Chem Pharm Bull. 2016;64:979–81.

    CAS  Article  Google Scholar 

  33. 33.

    Van Wagoner RM, Satake M, Wright JL. Polyketide biosynthesis in dinoflagellates: what makes it different? Nat Prod Rep. 2014;31:1101–37.

    PubMed  Article  PubMed Central  Google Scholar 

  34. 34.

    Walsh CT, O'Brien RV, Khosla C. Nonproteinogenic amino acid building blocks for nonribosomal peptide and hybrid Polyketide scaffolds. Angew Chem Int Ed. 2013;52:7098–124.

    CAS  Article  Google Scholar 

  35. 35.

    Jones AC, Monroe EA, Eisman EB, Gerwick L, Sherman DH, Gerwick WH. The unique mechanistic transformations involved in the biosynthesis of modular natural products from marine cyanobacteria. Nat Prod Rep. 2010;27:1048–65.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  36. 36.

    Wenzel SC, Muller R. Myxobacterial natural product assembly lines: fascinating examples of curious biochemistry. Nat Prod Rep. 2007;24:1211–24.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  37. 37.

    Lauritano C, De Luca D, Ferrarini A, Avanzato C, Minio A, Esposito F, et al. De novo transcriptome of the cosmopolitan dinoflagellate Amphidinium carterae to identify enzymes with biotechnological potential. Sci Rep. 2017;7:11701.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  38. 38.

    Lin S, Litaker RW, Sunda WG, Wood M. Phosphorus physiological ecology and molecular mechanisms in marine phytoplankton. J Phycol. 2016;52:10–36.

    CAS  PubMed  Article  Google Scholar 

  39. 39.

    Liu Z, Koid AE, Terrado R, Campbell V, Caron DA, Heidelberg KB. Changes in gene expression of Prymnesium parvum induced by nitrogen and phosphorus limitation. Front Microbiol. 2015;6:631.

    PubMed  PubMed Central  Google Scholar 

  40. 40.

    Han K, Lee H, Anderson DM, Kim B. Paralytic shellfish toxin production by the dinoflagellate Alexandrium pacificum (Chinhae Bay, Korea) in axenic, nutrient-limited chemostat cultures and nutrient-enriched batch cultures. Mar Pollut Bull. 2016;104:34–43.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  41. 41.

    Ianora A, Boersma M, Cassoti R, Fontana A, Harder J, Hoffmann F, et al. New trends in marine chemical ecology. Estuaries Coast. 2006;29:531–51.

    CAS  Article  Google Scholar 

  42. 42.

    Baumgarten S, Bayer T, Aranda M, Liew YJ, Carr A, Micklem G, et al. Integrating microRNA and mRNA expression profiling in Symbiodinium microadriaticum, a dinoflagellate symbiont of reef-building corals. BMC Genomics. 2013;14:704.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  43. 43.

    Gao D, Qiu L, Hou Z, Zhang Q, Wu J, Gao Q, Song L. Computational identification of microRNAs from the expressed sequence tags of toxic dinoflagellate Alexandrium Tamarense. Evol Bioinforma. 2013;9:479–85.

    CAS  Article  Google Scholar 

  44. 44.

    Geng H, Sui Z, Zhang S, Du Q, Ren Y, Liu Y, et al. Identification of microRNAs in the toxigenic dinoflagellate Alexandrium catenella by high-throughput Illumina sequencing and bioinformatic analysis. PLoS One. 2015;10:e0138709.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  45. 45.

    Dagenais-Bellefeuille S, Beauchemin, Morse, D miRNAs do not regulate circadian protein synthesis in the dinoflagellate Lingulodinium polyedrum PLoS ONE 2017; 12: e0168817.

  46. 46.

    Hopwood DA. Cracking the Polyketide code. PLoS Biol. 2004;2:e35.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  47. 47.

    Biswas S, Hazra S, Chattopadhyay S. Identification of conserved miRNAs and their putative target genes in Podophyllum hexandrum (Himalayan Mayapple). Plant Gene. 2016;6:82–9.

    CAS  Article  Google Scholar 

  48. 48.

    Liu J, Yuan Y, Wang Y, Jiang C, Chen T, Zhu F, et al. Regulation of fatty acid and flavonoid biosynthesis by miRNAs in Lonicera japonica. RSC Adv. 2017;7:35426–37.

    CAS  Article  Google Scholar 

  49. 49.

    Kalsotra A, Cooper TA. Functional consequences of developmentally regulated alternative splicing. Nature Rev Genet. 2011;12:715–29.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  50. 50.

    Shen S, Park JW, Lu ZX, Lin L, Henry MD, Wu YN, et al. rMATS: robust and flexible detection of differential alternative splicing from replicate RNA-Seq data. Proc Natl Acad Sci U S A. 2014;111:E5593–601.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  51. 51.

    Staiger D, Brown JW. Alternative splicing at the intersection of biological, development, and stress responses. Plant Cell. 2013;25:3640–56.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  52. 52.

    Zhu J, Wang X, Guo L, Xu Q, Zhao S, Li F, et al. Characterization and alternative splicing profiles of lipoxygenase gene family in tea plant (Camellia sinensis). Plant Cell Physiol. 2018;59:1765–81.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  53. 53.

    Seo PJ, Hong S-Y, Ryu JY, Jeong E-Y, Kim S-G, Baldwin IT, et al. Targeted inactivation of transcription factors by overexpression of their truncated forms in plants. Plant J. 2012;72:162–72.

    CAS  PubMed  Article  Google Scholar 

  54. 54.

    Monroe EA, Johnson JG, Wang Z, Pierce RK, Van Dolah FM. Characterization and expression of nuclear-encoded polyketide synthases in the brevetoxin-producing dinoflagellate Karenia brevis. J Phycol. 2010;46:541–52.

    CAS  Article  Google Scholar 

  55. 55.

    Hojo M, Omi A, Hamanaka G, Shindo K, Shimada A, Kondo M, et al. Unexpected link between polyketide synthase and calcium carbonate biomineralization. Zoological Lett. 2015;1:3.

    PubMed  PubMed Central  Article  Google Scholar 

  56. 56.

    Zhang H, Hou Y, Miranda L, Campbell DA, Sturm NR, Gaasterland T, et al. Spliced leader RNA trans-splicing in dinoflagellates. Proc Natl Acad Sci U S A. 2007;104:4618–23.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  57. 57.

    Blatch GL, Lassle M. The tetratricopeptide repeat: a structural motif mediating protein-protein interactions. Bioessays. 1999;21:932–9.

    CAS  PubMed  Article  Google Scholar 

  58. 58.

    Kobe B, Kajaba AV. The leucine-rich repeat as a protein recognition motif. Curr Opin Struct Biol. 2001;11:725–32.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  59. 59.

    Mosavi LK, Cammett TJ, Desrosiers DC, Peng ZY. The ankyrin repeat as molecular architechture for protein recognition. Protein Sci. 2004;13:1435–48.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  60. 60.

    Bretschneider T, Zocher G, Unger M, Scherlach K, Stehle T, Hertweck C. A ketosynthase homolog uses malonyl units to form esters in cervimycin biosynthesis. Nat Chem Biol. 2011;8:154–61.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  61. 61.

    Weissman KJ. Peering into the black box of fungal polyketide biosynthesis. ChemBioChem. 2010;11:485–8.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  62. 62.

    Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 30:1312–3.

  63. 63.

    Horiguchi T. Diversity and phylogeny of marine parasitic dinoflagellates. In: Ohtsuka S, Suzaki T, Horiguchi T, Suzuki N, Not F, editors. Marine protists: diversity and dynamics. Tokyo: Springer Japan; 2015. p. 397–419.

    Chapter  Google Scholar 

  64. 64.

    Marcais G, Kingsford C. A fast, lock-free approach for efficient parallel counting of occurrences of k-mers. Bioinformatics. 2011;27:764–70.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  65. 65.

    Vurture GW, Sedlazeck FJ, Nattestad M, Underwood CJ, Fang H, Gurtowski J, Schatz MC. GenomeScope: fast reference-free genome profiling from short reads. Bioinformatics. 2017;33:2202–4.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  66. 66.

    Doyle JJ, Doyle JL. A rapid DNA isolation procedure forsmall quantities of fresh leaf tissue. Phytochem Bull. 1987;19:11–5.

    Google Scholar 

  67. 67.

    Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  68. 68.

    Andrews S. FastQC: a quality control tool for high throughput sequence data. 2010; Available online at

  69. 69.

    Leggett RM, Clavijo BJ, Clissold L, Clark MD, Caccamo M. NextClip: an analysis and read preparation tool for Nextera Long mate pair libraries. Bioinformatics. 2014;30:566–8.

    CAS  PubMed  Article  Google Scholar 

  70. 70.

    Kajitani R, Toshimoto K, Noguchi H, Toyoda A, Ogura Y, Okuno M, et al. Efficient de novo assembly of highly heterozygous genomes from whole-genome shotgun short reads. Genome Res. 2014;24:1384–95.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  71. 71.

    Boetzer M, Henkel CV, Jansen HJ, Butler D, Pirovano W. Scaffolding pre-assembled contigs using SSPACE. Bioinformatics. 2011;27:578–9.

    CAS  PubMed  Article  Google Scholar 

  72. 72.

    Luo R, Liu B, Xie Y, Li Z, Huang W, Yuan J, et al. SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler. GigaScience. 2012;1:18.

    PubMed  PubMed Central  Article  Google Scholar 

  73. 73.

    Simão FA, Waterhouse RM, Ioannidis P, Kriventseva EV, Zdobnov EM. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics. 2015;31:3210–2.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  74. 74.

    Parra G, Bradnam K, Korf I. CEGMA: a pipeline to accurately annotate core genes in eukaryotic genomes. Bioinformatics. 2007;23:1061–7.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  75. 75.

    Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, et al. De novo transcript sequence reconstruction from RNA-seq using the trinity platform for reference generation and analysis. Nat Protoc. 2013;8:1494–512.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  76. 76.

    Hackl T, Hedrich R, Schultz J, Foerster F. Proovread: large-scale high accuracy PacBio correction through iterative short read consensus. Bioinformatics. 2014;30:3004–11.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  77. 77.

    Li W, Godzik A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006;22:1658–9.

    CAS  Article  Google Scholar 

  78. 78.

    Slamovits CH, Keeling PJ. Widespread recycling of processed cDNAs in dinoflagellates. Curr Biol. 2008;18:R550–2.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  79. 79.

    Wu TD, Watanabe CK. GMAP: a genomic mapping and alignment program for mRNA and EST sequences. Bioinformatics. 2005;21:1859–75.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  80. 80.

    Price AL, Jones NC, Pevzner PA. De novo identification of repeat families in large genomes. Bioinformatics. 2005;21:i351–8.

    CAS  PubMed  Article  Google Scholar 

  81. 81.

    Smit AFA, Hubley R, Green P. (1996–2010) RepeatMasker Open-3.0. (

  82. 82.

    Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29:15–21.

    CAS  PubMed  Article  Google Scholar 

  83. 83.

    Hoff KJ, Lange S, Lomsadze A, Borodovsky M, Stanke M. BRAKER1: unsupervised RNA-Seq-based genome annotation with GeneMark-ET and AUGUSTUS. Bioinformatics. 2016;32:767–9.

    CAS  PubMed  Article  Google Scholar 

  84. 84.

    Stanke M, Diekhans M, Baertsch R, Haussler D. Using native and syntenically mapped cDNA alignments to improve de novo gene finding. Bioinformatics. 2008;24:637–44.

    CAS  PubMed  Article  Google Scholar 

  85. 85.

    Magrane M, C. UniProt. UniProt Knowledgebase: a hub of integrated protein data. Database (Oxford), 2011; bar009.

  86. 86.

    Punta M, Coggill PC, Eberhardt RY, Mistry J, Tate J, Boursnell C, et al. The Pfam protein families database. Nucleic Acids Res. 2012;40:D290–301.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  87. 87.

    Finn RD, Clements J, Eddy SR. HMMER web server: interactive sequence similarity searching. Nucleic Acids Res. 2011;39:W29–37.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  88. 88.

    Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  89. 89.

    Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Höhna S. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61:539–42.

    PubMed  PubMed Central  Article  Google Scholar 

  90. 90.

    Zhang Y. I-TASSER server for protein 3D structure prediction. BMC Bioinformatics. 2008;9:40.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  91. 91.

    Baranašić D, Zucko J, Diminic J, Gacesa R, Long PF, Cullum J, et al. Predicting substrate specificity of adenylation domains of nonribosomal peptide synthetases and other protein properties by latent semantic indexing. J Ind Microbiol Biotechnol. 2014;41:461–7.

    PubMed  Article  CAS  Google Scholar 

  92. 92.

    Ziemert N, Podell S, Penn K, Badger JH, Allen E, Jensen PR. The natural product domain seeker NaPDoS: a phylogeny based bioinformatic tool to classify secondary metabolite gene diversity. PLoS One. 2012;7:e34064.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  93. 93.

    Emanuelsson O, Nielsen H, von Heijne G. ChloroP, a neural network-based method for predicting chloroplast transit peptides and their cleavage sites. Protein Sci. 1999;8:978–84.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  94. 94.

    Emanuelsson O, Brunak S, von Heijne G, Nielsen H. Locating proteins in the cell using TargetP, SignalP and related tools. Nat Protoc. 2007;2:953–71.

    CAS  PubMed  Article  Google Scholar 

  95. 95.

    Armenteros JJA, Sønderby CK, Sønderby SK, Nielsen H, Winther O. DeepLoc: prediction of protein subcellular localization using deep learning. Bioinformatics. 2017;33:3387–95.

    CAS  Article  Google Scholar 

  96. 96.

    Miranda KM, Espey MG, Wink DA. A rapid, simple spectrophotometric method for simultaneous detection of nitrate and nitrite. Nitric Oxide. 2001;5:62–71.

    CAS  PubMed  Article  Google Scholar 

  97. 97.

    Parsons TR. A manual of chemical & biological methods for seawater analysis. New York: Pergamon Press; 1984.

    Google Scholar 

  98. 98.

    Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2012;40:D109–14.

    CAS  PubMed  Article  Google Scholar 

  99. 99.

    Li B, Dewey CN. RSEM:accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12:323.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  100. 100.

    Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  101. 101.

    Alexa A, Rahnenfuhrer J. topGO: enrichment analysis for Gene Ontology. 2010; R package version 2.22.0.

  102. 102.

    Huang D, Sherman BT, Tan Q, Collins JR, Alvord WG, Roayaei J, et al. The DAVID gene functional classification tool: a novel biological module-centric algorithm to functionally analyze large gene lists. Genome Biol. 2007;8:R183.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  103. 103.

    Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17:10–2.

    Article  Google Scholar 

  104. 104.

    Friedländer MR, Mackowiak SD, Li N, Chen W, Rajewsky N. miRDeep2 accurately identifies known and hundreds of novel microRNA genes in seven animal clades. Nucleic Acids Res. 2012;40:37–52.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  105. 105.

    Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  106. 106.

    Enright AJ, John B, Gaul U, Tuschl T, Sander C, Marks DS. MicroRNA targets in drosophila. Genome Biol. 2003;5:R1.

    PubMed  PubMed Central  Article  Google Scholar 

  107. 107.

    Jones P, Binns D, Chang HY, Fraser M, Li W, McAnulla C, et al. InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014;30:1236–40.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  108. 108.

    Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li W, et al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal omega. Mol Sys Biol. 2014;7:539.

    Article  Google Scholar 

  109. 109.

    Clamp M, Cuff J, Searle SM, Barton GJ. The Jalview Java alignment editor. Bioinformatics. 2004;20:426–7.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  110. 110.

    Katz Y, Wang ET, Airoldi EM, Burge CB. Analysis and design of RNA sequencing experiments for identifying isoform regulation. Nat Methods. 2010;7:1009–15.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  111. 111.

    Anders S, Reyes A, Huber W. Detecting differential usage of exons from RNA-seq data. Genome Res. 2012;22:2008.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  112. 112.

    Falcon S, Gentleman R. Using GOstats to test gene lists for GO term association. Bioinformatics. 2007;23:257–8.

    CAS  PubMed  Article  Google Scholar 

  113. 113.

    Supek F, Bošnjak M, Škunca N, Šmuc T. REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS One. 2011;6:e21800.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  114. 114.

    Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T, et al. Fiji: an open-source platform for biological-image analysis. Nat Methods. 2012;9:676–82.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  115. 115.

    Beedessee G, Kubota T, Arimoto A, Nishitsuji K, Waller RF, Hisata K, et al. Integrated omics unveil the secondary metabolic landscape of a basal dinoflagellate. NCBI accession number PRJNA551917. 2020.

  116. 116.

    Beedessee G, Kubota T, Arimoto A, Nishitsuji K, Waller RF, Hisata K, et al. Integrated omics unveil the secondary metabolic landscape of a basal dinoflagellate. Amphidinium data repository. 2020.

Download references


The authors thank Ms. Haruhi Narisoko (OIST) for the assistance in culturing the alga, Dr. Toshio Sasaki and Dr. Koji Koizumi (IMG Section, OIST) for supporting microscopy imaging, Dr. Miyuki Kanda (SQC, OIST) for library preparation, and Dr. Frances van Dolah (College of Charleston, SC, USA) for providing PKS antibodies. We are also thankful to members of the Scientific Computing and Data Analysis of OIST for their support. We are grateful to anonymous reviewers for their valuable comments and to Dr. Steven D. Aird for editing the manuscript.


GB was supported by a Japanese Society for the Promotion of Science (JSPS) Research Fellowship for Young Scientists and a JSPS Grant-in-Aid for Fellows (17J00597). This work was supported by generous funding from Okinawa Institute of Science and Technology (OIST) Graduate University to the Marine Genomics Unit.

Author information




GB, ES, and NS conceptualized the study. GB, AA, KN, RFW, and ES analyzed the data and interpreted the results. GB and SY prepared the sequencing libraries. GB and ES prepared the figures and tables. KH, JK, and TK contributed reagents/analytic tools. GB and ES wrote the paper with input from all authors. The authors read and approved the final manuscript.

Corresponding author

Correspondence to Girish Beedessee.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Fig. S1.

Genome and transcriptome features of A. gibbosum. Fig. S2. Phylogenetic analysis of ketosynthase [KS], acyltransferase [AT] and condensation domains [C] using Bayesian inference. Fig. S3. Phylogenetic organization of adenylation domains from dinoflagellates. Fig. S4. Global expression profiles and enrichment of differentially expressed genes under nitrogen starvation (q-value < 0.001 and |log2(FC)| > 1). Fig. S5. Global expression profile and enrichment of differentially expressed genes under phosphate starvation (q-value < 0.001 and |log2(FC)| > 2). Fig. S6. Alignment of functional domains of the A. gibbosum homolog. Fig. S7. Length, distribution, and enrichment analysis of microRNAs detected from A. gibbosum. Fig. S8. Mapping of Illumina and Isoseq reads to g70808 and the presence of exons. Fig. S9. Immunofluorescent staining of Amphidinium with anti-KS and anti-KR antibodies. Fig. S10. Genome and transcriptome assembly workflows for Amphidinium gibbosum.

Additional file 2: Supplementary Table 1.

(a) Details of genome assembly based on statistics of scaffolds (b). Annotation statistics for gene models. Supplementary Table 2. The 30 most abundant domains in Amphidinium gibbosum. Supplementary Table 3. Amphidinium gibbosum repeat content. Supplementary Table 4. Comparison of major repeat content in Symbiodiniaceae and A. gibbosum. Supplementary Table 5. Top 10 KEGG pathways in A. gibbosum transcriptome. Supplementary Table 6. Significantly enriched KEGG pathways upregulated or downregulated under N and P starvation. Supplementary Table 7. miRNA KEGG pathway target enrichment under nitrogen and phosphate starvation. Supplementary Table 8. Details of miRNAs predicted from the A. gibbosum genome. Supplementary Table 9. Main differentially expressed genes during nutrient starvation in A. gibbosum, as shown in Fig. 3a. Supplementary Table 10. Annotation of PKS and NRPS genes under nitrogen and phosphate starvation, as shown in Fig. 3b.

Additional file 3.

Sequence alignment of the A-domain.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Beedessee, G., Kubota, T., Arimoto, A. et al. Integrated omics unveil the secondary metabolic landscape of a basal dinoflagellate. BMC Biol 18, 139 (2020).

Download citation


  • Polyketide synthases
  • Harmful algal blooms
  • Dinoflagellates
  • Iso-Seq
  • Duplication
  • Amphidinium