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

Environmental Factors Modulate Saxitoxins (STXs) Production in Toxic Dinoflagellate Alexandrium: An Updated Review of STXs and Synthesis Gene Aspects

The marine dinoflagellate Alexandrium is known to form harmful algal blooms (HABs) and produces saxitoxin (STX) and its derivatives (STXs) that cause paralytic shellfish poisoning (PSP) in humans. Cell growth and cellular metabolism are affected by environmental conditions, including nutrients, temperature, light, and the salinity of aquatic systems. Abiotic factors not only engage in photosynthesis, but also modulate the production of toxic secondary metabolites, such as STXs, in dinoflagellates. STXs production is influenced by a variety of abiotic factors; however, the relationship between the regulation of these abiotic variables and STXs accumulation seems not to be consistent, and sometimes it is controversial. Few studies have suggested that abiotic factors may influence toxicity and STXs-biosynthesis gene (sxt) regulation in toxic Alexandrium, particularly in A. catenella, A. minutum, and A. pacificum. Hence, in this review, we focused on STXs production in toxic Alexandrium with respect to the major abiotic factors, such as temperature, salinity, nutrients, and light intensity. This review informs future research on more sxt genes involved in STXs production in relation to the abiotic factors in toxic dinoflagellates.

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.

Development of saxitoxin biosynthesis gene sxtB-targeted qPCR assay for the quantification of toxic dinoflagellates Alexandrium catenella (group I) and A. pacificum (group IV) occurring in the Korean coast

Toxic dinoflagellate Alexandrium can produce saxitoxins (STXs) and cause paralytic shellfish poisoning (PSP), and thus they are monitored for environmental safety management. Microscopic discrimination of dinoflagellates is difficult to distinguish between toxic and non-toxic species due to their similar morphology. Meanwhile, an alternative quantitative PCR (qPCR) assay is sensitive, rapid, and cost-effective for harmful species monitoring. Herein, we developed a novel qPCR assay to detect the STXs biosynthesis gene sxtB of Alexandrium catenella and A. pacificum, the leading cause of PSP outbreaks in Asian coasts and worldwide. The newly designed sxtB TaqMan probes target the species without any positive signal in other relative dinoflagellates. Deming regression analysis revealed that the sxtB copy number of A. catenella and A. pacificum was 3.6 and 4.1 copies per cell, respectively. During the blooming periods (April 13th-14th, 2020), only A. catenella cells were detected through the qPCR assay, ranging from 5.0 × 10 to 2.5 × 104 eq cells L-1. In addition, sxtB qPCR quantified more accurately compared to large subunit (LSU) rRNA targeting qPCR assay that overestimate cell density. Besides, the sensitivity of sxtB was higher compared to the microscope when the species were rarely present (5.0 × 102 cells L-1). These suggest that the sxtB qPCR assay can be applied to toxic Alexandrium monitoring in the Korean coast, even in the early stage of bloomings.

Uncovering the regulation effect of modified clay on toxin production in Alexandrium pacificum: From physiological insights

As a common dinoflagellate, Alexandrium pacificum can produce paralytic shellfish toxins (PSTs). It can be removed from water by Polyaluminium chloride modified clay (PAC-MC), but it is unclear whether PAC-MC can prevent PSTs content and toxicity from increasing and whether PAC-MC can stimulate PSTs biosynthesis by A. pacificum. Effect of PAC-MC on PSTs and the physiological mechanism were analysed here. The results showed total PSTs content and toxicity decreased respectively by 34.10 % and 48.59 % in 0.2 g/L PAC-MC group at 12 days compared with control group. And the restriction of total PSTs by PAC-MC was mainly achieved via inhibition of algal cell proliferation, by affecting A. pacificum physiological processes and changing phycosphere microbial community. Meanwhile, single-cell PSTs toxicity did not increase significantly throughout the experiment. Moreover, A. pacificum treated with PAC-MC tended to synthesize sulfated PSTs such as C1&2. Mechanistic analysis showed that PAC-MC induced upregulation of sulfotransferase sxtN (related to PSTs sulfation), and functional prediction of bacterial community also showed significant enrichment of "sulfur relay system" after PAC-MC treatment, which might also promote PSTs sulfation. The results will provide theoretical guidance for the application of PAC-MC to field control of toxic Alexandrium blooms.

Bacterial communities and toxin profiles of Ostreopsis (Dinophyceae) from the Pacific Island of Okinawa, Japan

Variations in toxicity of the benthic dinoflagellate Ostreopsis Schmidt 1901 have been attributed to specific molecular clades, biogeography of isolated strains, and the associated bacterial community. Here, we attempted to better understand the biodiversity and the basic biology influencing toxin production of Ostreopsis. Nine clonal cultures were established from Okinawa, Japan, and identified using phylogenetic analysis of the ITS-5.8S rRNA and 28S rRNA genes. Morphological analysis suggests that the apical pore complex L/W ratio could be a feature for differentiating Ostreopsis sp. 2 from the O. ovata species complex. We analyzed the toxicity and bacterial communities using liquid chromatography-mass spectrometry, and PCR-free metagenomic sequencing. Ovatoxin was detected in three of the seven strains of O. cf. ovata extracts, highlighting intraspecies variation in toxin production. Additionally, two new potential analogs of ovatoxin-a and ostreocin-A were identified. Commonly associated bacteria clades of Ostreopsis were identified from the established cultures. While some of these bacteria groups may be common to Ostreopsis (Rhodobacterales, Flavobacteria-Sphingobacteria, and Enterobacterales), it was not clear from our analysis if any one or more of these plays a role in toxin biosynthesis. Further examination of biosynthetic pathways in metagenomic data and additional experiments isolating specific bacteria from Ostreopsis would aid these efforts.

Genetic association of toxin production in the dinoflagellate Alexandrium minutum

Dinoflagellates of the genus Alexandrium are responsible for harmful algal blooms and produce paralytic shellfish toxins (PSTs). Their very large and complex genomes make it challenging to identify the genes responsible for toxin synthesis. A family-based genomic association study was developed to determine the inheritance of toxin production in Alexandrium minutum and identify genomic regions linked to this production. We show that the ability to produce toxins is inheritable in a Mendelian way, while the heritability of the toxin profile is more complex. We developed the first dinoflagellate genetic linkage map. Using this map, several major results were obtained: 1. A genomic region related to the ability to produce toxins was identified. 2. This region does not contain any polymorphic sxt genes, known to be involved in toxin production in cyanobacteria. 3. The sxt genes, known to be present in a single cluster in cyanobacteria, are scattered on different linkage groups in A. minutum. 4. The expression of two sxt genes not assigned to any linkage group, sxtI and sxtG, may be regulated by the genomic region related to the ability to produce toxins. Our results provide new insights into the organization of toxicity-related genes in A. minutum, suggesting a dissociated genetic mechanism for the production of the different analogues and the ability to produce toxins. However, most of the newly identified genes remain unannotated. This study therefore proposes new candidate genes to be further explored to understand how dinoflagellates synthesize their toxins.

Effects of modified clay on the formation of Phaeocystis globosa colony revealed by physiological and transcriptomic analyses

The harmful algal bloom (HAB) species Phaeocystis globosa is commonly observed in global temperate and tropical oceans, and colonies of P. globosa exhibit a dominant morphotype during blooms. The use of polyaluminium chloride modified clay (PAC-MC) is an effective mitigation strategy for P. globosa blooms. Although previous studies have found that PAC-MC can stimulate P. globosa colony formation at low concentrations and inhibit it at higher concentrations, the underlying mechanisms of these effects are poorly understood. Here, we comprehensively compared the physiochemical indices and transcriptomic response of residual P. globosa cells after treatment with two concentrations of PAC-MC. The results showed that PAC-MC induced oxidative stress, photosynthetic inhibition, and DNA damage in residual cells. Moreover, it could activate antioxidant responses and enhance the repair of photosynthetic structure and DNA damage in cells. The biosynthesis of polysaccharides was enhanced and genes associated with cell motility were down-regulated after treatment with PAC-MC, resulting in the accumulation of colonial matrixes. After treatment with a low concentration of PAC-MC (0.1 g/L), the residual cells were slightly stressed, including physical damage, oxidative stress and other damage, and polysaccharide synthesis was enhanced to promote colony formation to alleviate environmental stress. Moreover, the damage to residual cells was slight; thus, normal cell function provided abundant energy and matter for colony formation. After treatment with a high concentration of PAC-MC (0.5 g/L), the residual cells suffered severe damage, which disrupted normal physiological processes and inhibited cell proliferation and colony formation. The present study elucidated the concentration-dependent mechanism of PAC-MC affecting the formation of P. globosa colonies and provided a reference for the application of PAC-MC to control P. globosa blooms.

Cyanotoxin-encoding genes as powerful predictors of cyanotoxin production during harmful cyanobacterial blooms in an inland freshwater lake: Evaluating a novel early-warning system

Freshwater harmful cyanobacterial blooms (HCBs) potentially produce excessive cyanotoxins, mainly microcystins (MCs), significantly threatening aquatic ecosystems and public health. Accurately predicting HCBs is thus essential to developing effective HCB mitigation and prevention strategies. We previously developed a novel early-warning system that uses cyanotoxin-encoding genes to predict cyanotoxin production in Harsha Lake, Ohio, USA, in 2015. In this study, we evaluated the efficacy of the early-warning system in forecasting the 2016 HCB in the same lake. We also examined potential HCB drivers and cyanobacterial community composition. Our results revealed that the cyanobacterial community was stable at the phylum level but changed dynamically at the genus level over time. Microcystis and Planktothrix were the major MC-producing genera that thrived in June and July and produced high concentrations of MCs (peak level 10.22 μg·L^-1). The abundances of the MC-encoding gene cluster mcy and its transcript levels significantly correlated with total MC concentrations (before the MC concentrations peaked) and accurately predicted MC production as revealed by logistic equations. When the Microcystis-specific gene mcyG reached approximately 1.5 × 10³ copies·mL^-1 or when its transcript level reached approximately 2.4 copies·mL^-1, total MC level exceeded 0.3 μg L^-1 (a health advisory limit) approximately one week later (weekly sampling scheme). This study suggested that cyanotoxin-encoding genes are promising predictors of MC production in inland freshwater lakes, such as Harsha Lake. The evaluated early-warning system can be a useful tool to assist lake managers in predicting, mitigating, and/or preventing HCBs.

Metal stresses modify soluble proteomes and toxin profiles in two Mediterranean strains of the distributed dinoflagellate Alexandrium pacificum

HABs involving Alexandrium pacificum have been reported in metal-contaminated ecosystems, suggesting that this distributed species adapts to and/or can tolerate the effects of metals. Modifications in soluble proteomes and PST contents were characterized in two Mediterranean A. pacificum strains exposed to mono- or polymetallic stresses (zinc, lead, copper, cadmium). These strains were isolated from two anthropized locations: Santa Giusta Lagoon (Italy, SG C10-3) and the Tarragona seaport (Spain, TAR C5-4F). In both strains, metals primarily downregulated key photosynthesis proteins. Metals also upregulated other proteins involved in photosynthesis (PCP in both strains), the oxidative stress response (HSP 60, proteasome and SOD in SG C10-3; HSP 70 in TAR C5-4F), energy metabolism (AdK in TAR C5-4F), neoglucogenesis/glycolysis (GAPDH and PEP synthase in SG C10-3) and protein modification (PP in TAR C5-4F). These proteins, possibly involved in adaptive proteomic responses, may explain the development of these A. pacificum strains in metal-contaminated ecosystems. The two strains showed different proteomic responses to metals, with SG C10-3 upregulating more proteins, particularly PCP. Among the PSTs, regardless of the metal and the strain studied, C2 and GTX4 predominated, followed by GTX5. Under the polymetallic cocktail, (i) total PSTs, C2 and GTX4 reached the highest levels in SG C10-3 only, and (ii) total PSTs, C2, GTX5 and neoSTX were higher in SG C10-3 than in TAR C5-4F, whereas in SG C10-3 under copper stress, total PSTs, GTX5, GTX1 and C1 were higher than in the controls, revealing variability in PST biosynthesis between the two strains. Total PSTs, C2, GTX4 and GTX1 showed significant positive correlations with PCP, indicating that PST production may be positively related to photosynthesis. Our results showed that the A. pacificum strains adapt their proteomic and physiological responses to metals, which may contribute to their ecological success in highly anthropized areas.

Unveiling the genomic structures and evolutionary events of the saxitoxin biosynthetic gene sxtA in the marine toxic dinoflagellate Alexandrium

Marine dinoflagellates Alexandriumare known to produce saxitoxin (STX) and cause paralytic shellfish poisoning (PSP) which can result in mortality in human. SxtA is considered a core gene for the biosynthesis of STX. However, its gene coding structure and evolutionary history have yet to be fully elucidated. Here, we determined the full-length sequences of sxtA cDNA and genomic coding regions from two toxic dinoflagellates, Alexandrium catenella (LIMS-PS-2645 and LIMS-PS-2647) andA. pacificum (LMBE-C4), characterised their domain structures, and resolved evolutionary events. The sxtA gene was encoded on the genome without introns, and was identical in length (4002 bp) between two A. catenella strains, but their sequences differed from A. pacificum (5031 bp). SxtA consists of four domains, sxtA1, sxtA2, sxtA3, and sxtA4; however, A. pacificum has an extra domain TauD near sxtA1. Each domain had >64.4% GC content, with the highest being 71.6% in sxtA3. Molecular divergence was found to be significantly higher in sxtA4 than in the other domains. Phylogenetic trees of sxtA and separate domains showed that bacteria diverged earliest, followed by non-toxic, toxic cyanobacteria, toxic dinoflagellates. While sxtA domains in Alexandrium were similar to the PKS-like structure with the conserved sxtA1, sxtA2, and sxtA3. PKS_KS may be replaced by sxtA4 in toxic Alexandrium. These suggest that sxtA in Alexandrium may have evolved by acquiring specific domains, whose modification and complexity markedly affect toxin biosynthesis.

Temperature influences the content and biosynthesis gene expression of saxitoxins (STXs) in the toxigenic dinoflagellate Alexandrium pacificum

Temperature may affect the production of saxitoxin (STX) and its derivatives (STXs); however, this is still controversial. Further, STX-biosynthesis gene regulation and the relation of its toxicity with temperature are not clearly understood. In the present study, we evaluated the effects of different temperatures (12 °C, 16 °C, and 20 °C) on the growth, toxin profiles, and expression of two core STX-biosynthesis genes, sxtA and sxtG, in the toxic dinoflagellate Alexandrium pacificum Alex05, isolated from Korean coasts. We found that temperature significantly affected cell growth, with maximum growth recorded at 16 °C, followed by 20 °C and 12 °C. HPLC analysis revealed mostly 12 of STXs from the tested cultures. Interestingly, the contents of STXs increased in the cells cultured at 16 °C and exposed to cold stress, compared to the 20 °C culture and heat stress; however, toxin components were much more diverse under heat stress. These toxin profiles generally matched with the sxtA and sxtG expression levels. Incubation at lower temperatures (12 °C and 16 °C) and exposure to cold stress increased sxtA and sxtG expressions in the cells, whereas heat stress showed little change or downregulated the transcription of both genes. Principal component analysis (PCA) showed low correlation between STXs eq and expressional levels of sxtA and sxtG in heat-stressed cells. These results suggest that temperature might be a crucial factor affecting the level and biosynthesis of STXs in marine toxic dinoflagellates.

Biochemical Mapping of Pyrodinium bahamense Unveils Molecular Underpinnings behind Organismal Processes

Proteins, lipids, and carbohydrates from the harmful algal bloom (HAB)-causing organism Pyrodinium bahamense were characterized to obtain insights into the biochemical processes in this environmentally relevant dinoflagellate. Shotgun proteomics using label-free quantitation followed by proteome mapping using the P. bahamense transcriptome and translated protein databases of Marinovum algicola, Alexandrium sp., Cylindrospermopsis raciborskii, and Symbiodinium kawagutii for annotation enabled the characterization of the proteins in P. bahamense. The highest number of annotated hits were obtained from M. algicola and highlighted the contribution of microorganisms associated with P. bahamense. Proteins involved in dimethylsulfoniopropionate (DMSP) degradation such as propionyl CoA synthethase and acryloyl-CoA reductase were identified, suggesting the DMSP cleavage pathway as the preferred route in this dinoflagellate. Most of the annotated proteins were involved in amino acid biosynthesis and carbohydrate degradation and metabolism, indicating the active roles of these molecules in the vegetative stage of P. bahamense. This characterization provides baseline information on the cellular machinery and the molecular basis of the ecophysiology of P. bahamense.

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

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

Salinity Affects Saxitoxins (STXs) Toxicity in the Dinoflagellate Alexandrium pacificum, with Low Transcription of SXT-Biosynthesis Genes sxtA4 and sxtG

Salinity is an important factor for regulating metabolic processes in aquatic organisms; however, its effects on toxicity and STX biosynthesis gene responses in dinoflagellates require further elucidation. Herein, we evaluated the physiological responses, toxin production, and expression levels of two STX synthesis core genes, sxtA4 and sxtG, in the dinoflagellate Alexandrium pacificum Alex05 under different salinities (20, 25, 30, 35, and 40 psu). Optimal growth was observed at 30 psu (0.12 cell division/d), but cell growth significantly decreased at 20 psu and was irregular at 25 and 40 psu. The cell size increased at lower salinities, with the highest size of 31.5 µm at 20 psu. STXs eq was highest (35.8 fmol/cell) in the exponential phase at 30 psu. GTX4 and C2 were predominant at that time but were replaced by GTX1 and NeoSTX in the stationary phase. However, sxtA4 and sxtG mRNAs were induced, and their patterns were similar in all tested conditions. PCA showed that gene transcriptional levels were not correlated with toxin contents and salinity. These results suggest that A. pacificum may produce the highest amount of toxins at optimal salinity, but sxtA4 and sxtG may be only minimally affected by salinity, even under high salinity stress.

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

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

Insights into Alexandrium minutum Nutrient Acquisition, Metabolism and Saxitoxin Biosynthesis through Comprehensive Transcriptome Survey

Simple Summary

Alexandrium minutum is one of the causing organisms for the occurrence of harmful algae bloom (HABs) in marine ecosystems. This species produces saxitoxin, one of the deadliest neurotoxins which can cause human mortality. However, molecular information such as genes and proteins catalog on this species is still lacking. Therefore, this study has successfully characterized several new molecular mechanisms regarding A. minutum environmental adaptation and saxitoxin biosynthesis. Ultimately, this study provides a valuable resource for facilitating future dinoflagellates’ molecular response to environmental changes.

Abstract

The toxin-producing dinoflagellate Alexandrium minutum is responsible for the outbreaks of harmful algae bloom (HABs). It is a widely distributed species and is responsible for producing paralytic shellfish poisoning toxins. However, the information associated with the environmental adaptation pathway and toxin biosynthesis in this species is still lacking. Therefore, this study focuses on the functional characterization of A. minutum unigenes obtained from transcriptome sequencing using the Illumina Hiseq 4000 sequencing platform. A total of 58,802 (47.05%) unigenes were successfully annotated using public databases such as NCBI-Nr, UniprotKB, EggNOG, KEGG, InterPRO and Gene Ontology (GO). This study has successfully identified key features that enable A. minutum to adapt to the marine environment, including several carbon metabolic pathways, assimilation of various sources of nitrogen and phosphorus. A. minutum was found to encode homologues for several proteins involved in saxitoxin biosynthesis, including the first three proteins in the pathway of saxitoxin biosynthesis, namely sxtA, sxtG and sxtB. The comprehensive transcriptome analysis presented in this study represents a valuable resource for understanding the dinoflagellates molecular metabolic model regarding nutrient acquisition and biosynthesis of saxitoxin.

Omics study of harmful algal blooms in China: Current status, challenges, and future perspectives

In the past two decades, the frequency, scale, and scope of harmful algal blooms (HABs) have increased significantly in the coastal waters of China. HABs have become a major ecological and environmental problem in China that seriously threatens the structure and function of marine ecosystems, the sustainable development of mariculture, and the health of human beings. Much effort has been devoted to studying HABs in China, and great achievements have been made in understanding the oceanographic and ecological mechanisms of HABs as well as the biology and physiological ecology of HAB-causing species. Furthermore, state-of-the-art omics technologies, such as transcriptomics and proteomics, have been used to elucidate the physiological responses of HAB-causing species to environmental changes, the biosynthesis of paralytic shellfish toxin, and the mechanisms underlying the formation of HABs. This review summarizes omics studies of HABs in China over the past few years and discusses challenges and future perspectives of HAB research.

sxtA4+ and sxtA4- Genotypes Occur Together within Natural Pyrodinium bahamense Sub-Populations from the Western Atlantic

Saxitoxin (STX) is a secondary metabolite and potent neurotoxin produced by several genera of harmful algal bloom (HAB) marine dinoflagellates. The basis for variability in STX production within natural bloom populations is undefined as both toxic and non-toxic strains (of the same species) have been isolated from the same geographic locations. Pyrodinium bahamense is a STX-producing bioluminescent dinoflagellate that blooms along the east coast of Florida as well as the bioluminescent bays in Puerto Rico (PR), though no toxicity reports exist for PR populations. The core genes in the dinoflagellate STX biosynthetic pathway have been identified, and the sxtA4 gene is essential for toxin production. Using sxtA4 as a molecular proxy for the genetic capacity of STX production, we examined sxtA4+ and sxtA4- genotype frequency at the single cell level in P. bahamense populations from different locations in the Indian River Lagoon (IRL), FL, and Mosquito Bay (MB), a bioluminescent bay in PR. Multiplex PCR was performed on individual cells with Pyrodinium-specific primers targeting the 18S rRNA gene and sxtA4. The results reveal that within discrete natural populations of P. bahamense, both sxtA4+ and sxtA4- genotypes occur, and the sxtA4+ genotype dominates. In the IRL, the frequency of the sxtA4+ genotype ranged from ca. 80–100%. In MB, sxtA4+ genotype frequency ranged from ca 40–66%. To assess the extent of sxtA4 variation within individual cells, sxtA4 amplicons from single cells representative of the different sampling sites were cloned and sequenced. Overall, two variants were consistently obtained, one of which is likely a pseudogene based on alignment with cDNA sequences. These are the first data demonstrating the existence of both genotypes in natural P. bahamense sub-populations, as well as sxtA4 presence in P. bahamense from PR. These results provide insights on underlying genetic factors influencing the potential for toxin variability among natural sub-populations of HAB species and highlight the need to study the genetic diversity within HAB sub-populations at a fine level in order to identify the molecular mechanisms driving HAB evolution.

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.

SxtA localizes to chloroplasts and changes to its 3'UTR may reduce toxin biosynthesis in non-toxic Alexandrium catenella (Group I)✰

SxtA is the enzyme that catalyses the first step of saxitoxin biosynthesis. We developed an immunofluorescent method to detect SxtA using antibodies against SxtA peptides. Confocal microscopy revealed the presence of abundant, sub-cellularly localized signal in cells of toxic species and its absence in non-toxic species. Co-localization of SxtA with Rubisco II and ultra-structural observation by transmission electron microscopy strongly suggested the association of SxtA with chloroplasts. We also characterized a non-toxic sub-clone of Alexandrium catenella (Group I) to elucidate the mutation responsible for its loss of toxicity. Although sxtA4 gene copy number was indistinguishable in toxic and non-toxic sub-clones, mRNA and protein expression were significantly reduced in the non-toxic sub-clone and we uncovered sequence variation at the 3' untranslated region (3'UTR) of sxtA4 mRNA. We propose that differences in the sxtA4 mRNA 3'UTR lead to down-regulation of STX biosynthesis post-transcriptionally, thereby explaining the differences in toxicity amongst different A. catenella (Group I) sub-clones.

De novo Transcriptome of the Non-saxitoxin Producing Alexandrium tamutum Reveals New Insights on Harmful Dinoflagellates

Many dinoflagellates species, especially of the Alexandrium genus, produce a series of toxins with tremendous impacts on human and environmental health, and tourism economies. Alexandrium tamutum was discovered for the first time in the Gulf of Naples, and it is not known to produce saxitoxins. However, a clone of A. tamutum from the same Gulf showed copepod reproduction impairment and antiproliferative activity. In this study, the full transcriptome of the dinoflagellate A. tamutum is presented in both control and phosphate starvation conditions. RNA-seq approach was used for in silico identification of transcripts that can be involved in the synthesis of toxic compounds. Phosphate starvation was selected because it is known to induce toxin production for other Alexandrium spp. Results showed the presence of three transcripts related to saxitoxin synthesis (sxtA, sxtG and sxtU), and others potentially related to the synthesis of additional toxic compounds (e.g., 44 transcripts annotated as "polyketide synthase"). These data suggest that even if this A. tamutum clone does not produce saxitoxins, it has the potential to produce toxic metabolites, in line with the previously observed activity. These data give new insights into toxic microalgae, toxin production and their potential applications for the treatment of human pathologies.

Comparative transcriptomics reveals colony formation mechanism of a harmful algal bloom species Phaeocystis globosa

Phaeocystis globosa is a major causative agent of harmful algal blooms in the global ocean, featuring a complex polymorphic life cycle alternating between free-living solitary cells and colonial cells. Colony is the dominant morphotype during P. globosa bloom. However, the underlying mechanism of colony formation is poorly understood. Here, we comprehensively compared global transcriptomes of P. globosa cells at four distinctive colony formation stages: free-living solitary cells, two cell-, four cell- and multi-cell colonies, under low (20 °C) and high (32 °C) temperatures, and characterized the genes involved in colony formation. Glycosaminoglycan (GAG) synthesis was enhanced while its degradation was decreased during colony formation, resulting in the accumulation of GAGs that are an essential substrate of the colony matrix. Nitrogen metabolism and glutamine synthesis were remarkably increased in the colonial cells, which provided precursors for GAG synthesis. Furthermore, cell defense and motility were down-regulated in the colonial cells, thereby conserving energy for GAG synthesis. Notably, high temperature led to decreased synthesis and increased degradation of GAGs, resulting in insufficient substrates to form the colony. Our study indicates that GAGs accumulation is critical for colony formation of P. globosa, but high temperature inhibits GAGs' accumulation and colony formation.

Transcriptome survey and toxin measurements reveal evolutionary modification and loss of saxitoxin biosynthesis genes in the dinoflagellates Amphidinium carterae and Prorocentrum micans

In the present study, we characterized the potential toxin genes for polyketide synthase (PKS) and saxitoxin (STX) biosynthesis using the transcriptomes of two non-STX producing dinoflagellates Amphidinium carterae and Prorocentrum micans. RNA sequencing revealed 94 and 166 PKS contigs in A. carterae and P. micans, respectively. We first detected type III PKS, which was closely related to bacteria. In addition, dozens of homologs of 20 STX biosynthesis genes were identified. Interestingly, the core STX-synthesizing genes sxtA and sxtB were only found in P. micans, whereas sxtD was detected in A. carterae alone. Bioinformatic analysis showed that the first two core genes (sxtA and sxtG) had a low sequence similarity (37.0-67.6%) and different domain organization compared to those of other toxigenic dinoflagellates, such as Alexandrium pacificum. These might result in the breakdown of the initial reactions in STX production and ultimately the loss of the ability to synthesize the toxins in both dinoflagellates. Our findings suggest that toxin-related PKS and sxt genes are commonly found in non-STX producing dinoflagellates. In addition to their involvement in the synthesis of toxins, our result indicates that genes may also have other molecular metabolic functions.

Characterizing ciguatoxin (CTX)- and Non-CTX-producing strains of Gambierdiscus balechii using comparative transcriptomics

Gambierdiscus spp. can produce the polyketide compound, ciguatoxin (CTX), and are hence responsible for ciguatera fish poisoning (CFP). Studying the molecular mechanism that regulates CTX production is crucial for understanding the environmental trigger of CTX as well as for better informing fishery management. Commonly, polyketide synthases are important for polyketide synthesis; however, no gene has been confirmatively assigned to CTX production. Here, suppression subtractive hybridization (SSH) and transcriptome sequencing (RNA-Seq) were used to compare a CTX-producing strain with a non-CTX-producing strain. Using both methods, a total of 52 polyketide synthase (PKS) genes were identified to be up-regulated in the CTX-producing G. balechii, including transcripts encoding single-domain PKSs as well as transcripts encoding multi-domain PKSs. Using reverse transcription quantitative PCR, the expression of these genes in the CTX-producing strain and in nitrogen-limited cultures of the strain was further documented. These data suggest that PKSs are likely involved in polyketide synthesis and potentially in CTX synthesis in this dinoflagellate species. Our study provides the candidate biomarkers for the detection of CTXs or CFP in waters or any other organisms as well as a valuable genomic resource for the research on Gambierdiscus and other dinoflagellates.

Intraspecific variability in membrane proteome, cell growth, and morphometry of the invasive marine neurotoxic dinoflagellate Alexandrium pacificum grown in metal-contaminated conditions

Over the past decades, the occurrence, distribution and intensity of harmful algal blooms involving the dinoflagellate Alexandrium pacificum have increased in marine coastal areas disturbed by anthropogenic inputs. This invasive species produces saxitoxin, which causes the paralytic shellfish poisoning syndrome in humans upon consumption of contaminated seafood. Blooms of A. pacificum have been reported in metal-contaminated coastal ecosystems, suggesting some ability of these microorganisms to adapt to and/or resist in metal stress conditions. This study seeks to characterize the modifications in membrane proteomes (by 2-D electrophoresis coupled to LC-MS/MS), cell growth and morphometry (measured with an inverted microscope), in response to metal stress (addition of Zn²⁺, Pb²⁺, Cu²⁺ and Cd²⁺), in two Mediterranean A. pacificum strains: SG C10-3 and TAR C5-4F, respectively isolated from the Santa Giusta Lagoon (Sardinia, Italy) and from the Tarragona seaport (Spain), both metal-contaminated ecosystems. In the SG C10-3 cultures grown in a metal cocktail, cell growth was significantly delayed, and cell size increased (22% of 37.5 μm cells after 25 days of growth). Conversely, no substantial change was observed for cell growth or cell size in the TAR C5-4F cultures grown in a metal cocktail (P > 0.10), thus indicating intraspecific variability in the responses of A. pacificum strains to metal contamination. Regardless of the conditions tested, the total number of proteins constituting the membrane proteome was significantly higher for TAR C5-4F than for SG C10-3, which may help TAR C5-4F to thrive better in contaminated conditions. For both strains, the total number of proteins constituting the membrane proteomes was significantly lower in response to metal stress (29% decrease in the SG C10-3 proteome: 82 ± 12 proteins for controls, and 58 ± 12 in metal-contaminated cultures; 17% decrease in the TAR C5-4F proteome: 101 ± 8 proteins for controls, and 84 ± 5 in metal-contaminated cultures). Moreover, regardless of the strain, proteins with significantly modified expression in response to stress were mainly down-regulated (representing 45% of the proteome for SG C10-3 and 38% for TAR C5-4F), clearly showing the harmful effects of the metals. Protein down-regulation may affect cell transport (actin and phospholipid scramblase in SG C10-3), photosynthesis (RUBISCO in SG C10-3, light-harvesting protein in TAR C5-4F, and high-CO₂-inducing periplasmic protein in both strains), and finally energy metabolism (ATP synthase in both strains). However, other modifications in protein expression may confer to these A. pacificum strains a capacity for adaptation and/or resistance to metal stress conditions, for example by (i) limiting the metal entry through the plasma membrane of the SG C10-3 cells (via the down-regulation of scramblase) and/or (ii) reducing the oxidative stress generated by metals in SG C10-3 and TAR C5-4F cells (due to down-regulation of ATP-synthase).

Comparative transcriptomics of toxin synthesis genes between the non-toxin producing dinoflagellate Cochlodinium polykrikoides and toxigenic Alexandrium pacificum

In the present study, we extensively characterized potential toxin-related genes, including polyketide synthase (PKS), saxitoxin (STX) and fatty acid synthase (FAS) from the non-toxin producing marine dinoflagellate Cochlodinium polykrikoides, comparing to those of a toxigenic dinoflagellate Alexandrium pacificum. RNA sequencing revealed 50 and 271 PKS contigs from C. polykrikoides and A. pacificum, respectively. According to domain constitute and amino acid alteration, we further classified the dinoflagellate type I PKS genes into 4 sub-groups. Type III PKS was first identified in C. polykrikoides. Interestingly, we detected a large number (242 and 288) of homologs of 18 sxt genes from two studied dinoflagellates. Most of the eight key genes (sxtA, sxtB, sxtD, sxtG, sxtH/T, sxtI, sxtS and sxtU) for STX synthesis were detected in both dinoflatellates, whereas a core STX biosynthesis gene sxtG was not detected in C. polykrikoides. This may partially explain the absence of saxitoxin production in C. polykrikoides. In addition, we identified several type I and type II FAS genes, including FabD, FabF, FabG, FabH, FabI, and FabZ, whereas FabB was not found in C. polykrikoides. Overall, the numbers of the toxin-related genes in C. polykrikoides were less than that of A. pacificum. Phylogenetic analyses showed that type I PKS/FASs of dinoflagellates had close relationships with apicomplexans and bacteria. These suggest that the toxin-related PKS and sxt genes are commonly present in toxigenic and non-toxin producing dinoflagellates, and may be involved not only in the toxin synthesis, but also in other related molecular metabolic functions.

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.

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

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

Pyrene metabolites by bacterium enhancing cell division of green alga Selenastrum capricornutum

Collaborations between multiple microbial species are important for understanding natural clearance and ecological effects of toxic organic contaminants in the environment. However, the interactions between different species in the transformation and degradation of contaminants remain to address. In this study, the effects of pyrene and its bacterial metabolites on the algal growth (Selenastrum capricornutum) were examined. The specific growth rate of algal cells incubated with bacterial pyrene metabolites (1.18 d^-1) was highest among all treatment, followed by the controls (1.07 d^-1), treated with pyrene-free bacterial metabolites (1.04 d^-1) and those treated with pyrene (0.55 d^-1). G₁ phase is the key growth phase for the cells to synthesize biomolecules for subsequent cell division in the cell cycle. Approximately 76.9% of the cells treated with bacterial pyrene metabolites were at the G₁ phase and significantly lower than those with the controls (85.3%), pyrene-free bacterial metabolites (85.5%) and pyrene treatment (92.5%). Transcriptomic analysis of algae showed that the expression of 47 ribosomal unigenes was down-regulated by 5 mg L^-1 of pyrene, while 308 unigenes related to the preparation of cell division (DNA replication and protein synthesis) were up-regulated by bacterial pyrene metabolites. It indicated that basal metabolism associated with the growth and proliferation of algal cells could be significantly promoted by bacterial pyrene metabolites. Overall, this study suggests a close relationship between algae and bacteria in the transformation and ecological effects of toxic contaminants.

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.

In silico prediction of the secretome from the invasive neurotoxic marine dinoflagellate Alexandrium catenella

Alexandrium catenella, a marine dinoflagellate responsible for harmful algal blooms (HABs), proliferates with greater frequency, distribution and intensity, in disturbed marine coastal ecosystems. The proteins secreted into seawater may play a crucial role in maintaining this dinoflagellate in these ecosystems, but this possibility has never been investigated before. In this study, the A. catenella secretome was predicted from its transcriptome by combining several bioinformatics tools. Our results predict a secretome of 2 779 proteins, among which 79% contain less than 500 amino acids, suggesting that most secreted proteins are short in length. The predicted secretome includes 963 proteins (35%) with Pfam domains: 773 proteins with one Pfam domain and 190 proteins with two or more Pfam domains. Their functional annotation showed that they are mainly involved in (i) proteolysis, (ii) stress responses and (iii) primary metabolism. In addition, 47% of the secreted proteins appear to be enzymes, primarily peptidases, known to be biologically active in the extracellular medium during stress responses. Finally, this study provides a wealth of candidates of proteins secreted by A. catenella, which may interact with the marine environment and help this dinoflagellate develop in various environmental conditions.

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.

iTRAQ-Based Quantitative Proteomic Analysis of a Toxigenic Dinoflagellate Alexandrium catenella at Different Stages of Toxin Biosynthesis during the Cell Cycle

Paralytic shellfish toxins (PSTs) are a group of potent neurotoxic alkaloids that are produced mainly by marine dinoflagellates. PST biosynthesis in dinoflagellates is a discontinuous process that is coupled to the cell cycle. However, little is known about the molecular mechanism underlying this association. Here, we compared global protein expression profiles of a toxigenic dinoflagellate, Alexandrium catenella, collected at four different stages of toxin biosynthesis during the cell cycle, using an isobaric tags for relative and absolute quantification (iTRAQ)-based quantitative proteomic approach. The results showed that toxin biosynthesis occurred mainly in the G1 phase, especially the late G1 phase. In total, 7232 proteins were confidently identified, and 210 proteins exhibited differential expression among the four stages. Proteins involved in protein translation and photosynthetic pigment biosynthesis were significantly upregulated during toxin biosynthesis, indicating close associations among the three processes. Nine toxin-related proteins were detected, and two core toxin biosynthesis proteins, namely, sxtA and sxtI, were identified for the first time in dinoflagellates. Among these proteins, sxtI and ompR were significantly downregulated when toxin biosynthesis stopped, indicating that they played important roles in the regulation of PST biosynthesis. Our study provides new insights into toxin biosynthesis in marine dinoflagellates: nitrogen balance among different biological processes regulates toxin biosynthesis, and that glutamate might play a key modulatory role.

Venomix: a simple bioinformatic pipeline for identifying and characterizing toxin gene candidates from transcriptomic data

The advent of next-generation sequencing has resulted in transcriptome-based approaches to investigate functionally significant biological components in a variety of non-model organism. This has resulted in the area of “venomics”: a rapidly growing field using combined transcriptomic and proteomic datasets to characterize toxin diversity in a variety of venomous taxa. Ultimately, the transcriptomic portion of these analyses follows very similar pathways after transcriptome assembly often including candidate toxin identification using BLAST, expression level screening, protein sequence alignment, gene tree reconstruction, and characterization of potential toxin function. Here we describe the Python package Venomix, which streamlines these processes using common bioinformatic tools along with ToxProt, a publicly available annotated database comprised of characterized venom proteins. In this study, we use the Venomix pipeline to characterize candidate venom diversity in four phylogenetically distinct organisms, a cone snail (Conidae; Conus sponsalis), a snake (Viperidae; Echis coloratus), an ant (Formicidae; Tetramorium bicarinatum), and a scorpion (Scorpionidae; Urodacus yaschenkoi). Data on these organisms were sampled from public databases, with each original analysis using different approaches for transcriptome assembly, toxin identification, or gene expression quantification. Venomix recovered numerically more candidate toxin transcripts for three of the four transcriptomes than the original analyses and identified new toxin candidates. In summary, we show that the Venomix package is a useful tool to identify and characterize the diversity of toxin-like transcripts derived from transcriptomic datasets. Venomix is available at: https://bitbucket.org/JasonMacrander/Venomix/.

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.

Exploring dinoflagellate biology with high-throughput proteomics

Dinoflagellates are notorious for their ability to form the harmful algal blooms known as "red tides," yet the mechanisms underlying bloom formation remain poorly understood. Despite recent advances in nucleic acid sequencing, which have generated transcriptomes from a wide range of species exposed to a variety of different conditions, measuring changes in RNA levels have not generally produced great insight into dinoflagellate cell biology or environmental physiology, nor do we have a thorough grasp on the molecular events underpinning bloom formation. Not only is the transcriptomic response of dinoflagellates to environmental change generally muted, but there is a markedly low degree of congruency between mRNA expression and protein expression in dinoflagellates. Herein we discuss the application of high-throughput proteomics to the study of dinoflagellate biology. By profiling the cellular protein complement (the proteome) instead of mRNA (the transcriptome), the biomolecular events that underlie the changes of phenotypes can be more readily evaluated, as proteins directly determine the structure and the function of the cell. Recent advances in proteomics have seen this technique become a high-throughput method that is now able to provide a perspective different from the more commonly employed nucleic acid sequencing. We suggest that the time is ripe to exploit these new technologies in addressing the many mysteries of dinoflagellate biology, such as how the symbiotic dinoflagellate inhabiting reef corals acclimate to increases in temperature, as well as how harmful algal blooms are initiated at the sub-cellular level. Furthermore, as dinoflagellates are not the only eukaryotes that demonstrate muted transcriptional responses, the techniques addressed within this review are amenable to a wide array of organisms.

iTRAQ-Based Quantitative Proteomic Analysis of a Toxigenic Dinoflagellate Alexandrium catenella and Its Non-toxigenic Mutant Exposed to a Cell Cycle Inhibitor Colchicine

Paralytic shellfish toxins (PSTs) are a group of potent neurotoxic alkaloids mainly produced by marine dinoflagellates and their biosynthesis is associated with the cell cycle. Study shows that colchicine can cease cell division and inhibit PST production of dinoflagellates. However, the molecular mechanism behind this linkage is unknown. Here, we applied the iTRAQ-based proteomic approach to investigate protein expression profiles of a toxigenic dinoflagellate Alexandrium catenella (ACHK-T) and its non-toxigenic mutant (ACHK-NT) when treated with colchicine. The results showed that the cell cycles of both strains were arrested at the G1 phase by colchicine, and the toxin biosynthesis of ACHK-T was inhibited. Among 6,988 proteins identified, 113 and 253 proteins were differentially expressed in the colchicine-treated ACHK-T and ACHK-NT, respectively, compared with their non-colchicine treatments. Proteins involved in reactive oxygen species scavenging and protein degradation were upregulated in both strains while proteins participating in photosynthetic pigment biosynthesis and nitrogen metabolism presented different expressions. Nitrate reductase and glutamine synthetase were altered insignificantly in the colchicine-treated ACHK-T while both of them were remarkably downregulated in the colchicine-treated ACHK-NT, suggesting a feedback regulation between PST production and nitrogen metabolism in ACHK-T. Nitrogen originally for PST biosynthesis might be reallocated to photosynthetic pigment biosynthesis in the colchicine-treated ACHK-T. A total of 55 homologs of 7 toxin-related proteins were obtained; however, they altered insignificantly in both colchicine-treated strains, suggesting that toxin biosynthesis might be post-translationally regulated. Our study provided new insights into toxin biosynthesis in marine dinoflagellates.

Transcriptome sequencing and characterization of Symbiodinium muscatinei and Elliptochloris marina, symbionts found within the aggregating sea anemone Anthopleura elegantissima

There is a growing body of literature using transcriptomic data to study how tropical cnidarians and their photosynthetic endosymbionts respond to environmental stressors and participate in metabolic exchange. Despite these efforts, our understanding of how essential genes function to facilitate symbiosis establishment and maintenance remains limited. The inclusion of taxonomically and ecologically diverse endosymbionts will enhance our understanding of these interactions. Here we characterize the transcriptomes of two very different symbionts found within the temperate sea anemone Anthopleura elegantissima: the chlorophyte Elliptochloris marina and the dinoflagellate Symbiodinium muscatinei. We use a multi-level approach to assess the diversity of genes found across S. muscatinei and E. marina transcriptomes, and compare their overall protein domains with other dinoflagellates and chlorophytes. Our analysis identified several genes that are potentially involved in mitigating stress response (e.g., heat shock proteins pathways for mediating reactive oxygen species) and metabolic exchange (e.g., ion transporters). Finally, we show that S. muscatinei and other Symbiodinium strains are equipped with a high salt peridinin-chl-protein (HSPCP) gene previously identified only in free-living dinoflagellates. The addition of these transcriptomes to the cnidarian-symbiont molecular toolkit will aid in understanding how these vitally important symbiotic relationships are established and maintained across a variety of environmental conditions.

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

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

Genome complexity of harmful microalgae

During the past decade, next generation sequencing (NGS) technologies have provided new insights into the diversity, dynamics, and metabolic pathways of natural microbial communities. But, these new techniques face challenges related to the genome size and level of genome complexity of the species under investigation. Moreover, the coverage depth and the short-read length achieved by NGS based approaches also represent a major challenge for assembly. These factors could limit the use of these high-throughput sequencing methods for species lacking a reference genome and characterized by a high level of complexity. In the present work, the evolutionary history, mainly consisting of gene transfer events from bacteria and unicellular eukaryotes to microalgae, including harmful species, is discussed and reviewed as it relates to NGS application in microbial communities, with a particular focus on harmful algal bloom species and dinoflagellates. In the context of genetic population studies, genotyping-by-sequencing (GBS), an NGS based approach, could be used for the discovery and analysis of single nucleotide polymorphisms (SNPs). The NGS technologies are still relatively new and require further improvement. Specifically, there is a need to develop and standardize tools and approaches to handle large data sets, which have to be used for the majority of HAB species characterized by evolutionary highly dynamic genomes.

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.

Comparative Transcriptomic Analysis Reveals Novel Insights into the Adaptive Response of Skeletonema costatum to Changing Ambient Phosphorus

Phosphorus (P) is a limiting macronutrient for diatom growth and productivity in the ocean. Much effort has been devoted to the physiological response of marine diatoms to ambient P change, however, the whole-genome molecular mechanisms are poorly understood. Here, we utilized RNA-Seq to compare the global gene expression patterns of a marine diatom Skeletonema costatum grown in inorganic P-replete, P-deficient, and inorganic- and organic-P resupplied conditions. In total 34,942 unique genes were assembled and 20.8% of them altered significantly in abundance under different P conditions. Genes encoding key enzymes/proteins involved in P utilization, nucleotide metabolism, photosynthesis, glycolysis, and cell cycle regulation were significantly up-regulated in P-deficient cells. Genes participating in circadian rhythm regulation, such as circadian clock associated 1, were also up-regulated in P-deficient cells. The response of S. costatum to ambient P deficiency shows several similarities to the well-described responses of other marine diatom species, but also has its unique features. S. costatum has evolved the ability to re-program its circadian clock and intracellular biological processes in response to ambient P deficiency. This study provides new insights into the adaptive mechanisms to ambient P deficiency in marine diatoms.

Transcriptome of the Australian Mollusc Dicathais orbita Provides Insights into the Biosynthesis of Indoles and Choline Esters

Dicathais orbita is a mollusc of the Muricidae family and is well known for the production of the expensive dye Tyrian purple and its brominated precursors that have anticancer properties, in addition to choline esters with muscle-relaxing properties. However, the biosynthetic pathways that produce these secondary metabolites in D. orbita are not known. Illumina HiSeq 2000 transcriptome sequencing of hypobranchial glands, prostate glands, albumen glands, capsule glands, and mantle and foot tissues of D. orbita generated over 201 million high quality reads that were de novo assembled into 219,437 contigs. Annotation with reference to the Nr, Swiss-Prot and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases identified candidate-coding regions in 76,152 of these contigs, with transcripts for many enzymes in various metabolic pathways associated with secondary metabolite biosynthesis represented. This study revealed that D. orbita expresses a number of genes associated with indole, sulfur and histidine metabolism pathways that are relevant to Tyrian purple precursor biosynthesis, and many of which were not found in the fully annotated genomes of three other molluscs in the KEGG database. However, there were no matches to known bromoperoxidase enzymes within the D. orbita transcripts. These transcriptome data provide a significant molecular resource for gastropod research in general and Tyrian purple producing Muricidae in particular.

Physiological and transcriptional responses to inorganic nutrition in a tropical Pacific strain of Alexandrium minutum: Implications for the saxitoxin genes and toxin production

Saxitoxins (STXs) constitute a family of potent sodium channel blocking toxins, causative agents of paralytic shellfish poisoning (PSP), and are produced by several species of marine dinoflagellates and cyanobacteria. Two STX-core genes, sxtA and sxtG, have been well elucidated in Alexandrium but the expression of these genes under various nutritional modes in tropical species remains unclear. This study investigates the physiological responses of a tropical Pacific strain of Alexandrium minutum growing with nitrate or ammonium, and with various nitrogen to phosphorus (N:P) supply ratios. The transcriptional responses of the sxt genes were observed. Likewise, a putative sxtI encoding O-carbamoyltransferase (herein designated as AmsxtI) was recovered from the transcriptomic data, and its expression was investigated. The results revealed that the cellular toxin quota (Q_t) was higher in P-depleted, nitrate-grown cultures. With cultures at similar N:P (<16), cells grown with excess ammonium showed a higher Q_t than those grown with nitrate. sxtA1 was not expressed under any culture conditions, suggesting that this gene might not be involved in STX biosynthesis by this strain. Conversely, sxtA4 and sxtG showed positive correlations with Q_t, and were up-regulated in P-depleted, nitrate-grown cultures and with excess ambient ammonium. On the other hand, AmsxtI was expressed only when induced by P-depletion, suggesting that this gene may play an important role in P-recycling metabolism, while simultaneously enhancing toxin production.

Marine natural products

This review covers the literature published in 2014 for marine natural products (MNPs), with 1116 citations (753 for the period January to December 2014) referring to compounds isolated from marine microorganisms and phytoplankton, green, brown and red algae, sponges, cnidarians, bryozoans, molluscs, tunicates, echinoderms, mangroves and other intertidal plants and microorganisms. The emphasis is on new compounds (1378 in 456 papers for 2014), together with the relevant biological activities, source organisms and country of origin. Reviews, biosynthetic studies, first syntheses, and syntheses that lead to the revision of structures or stereochemistries, have been included.

Quantity of the dinoflagellate sxtA4 gene and cell density correlates with paralytic shellfish toxin production in Alexandrium ostenfeldii blooms

Many marine dinoflagellates, including several species of the genus Alexandrium, Gymnodinium catenatum, and Pyrodinium bahamense are known for their capability to produce paralytic shellfish toxins (PST), which can cause severe, most often food-related poisoning. The recent discovery of the first PST biosynthesis genes has laid the foundation for the development of molecular detection methods for monitoring and study of PST-producing dinoflagellates. In this study, a probe-based qPCR method for the detection and quantification of the sxtA4 gene present in Alexandrium spp. and Gymnodinium catenatum was designed. The focus was on Alexandrium ostenfeldii, a species which recurrently forms dense toxic blooms in areas within the Baltic Sea. A consistent, positive correlation between the presence of sxtA4 and PST biosynthesis was observed, and the species was found to maintain PST production with an average of 6 genomic copies of sxtA4. In August 2014, A. ostenfeldii populations were studied for cell densities, PST production, as well as sxtA4 and species-specific LSU copy numbers in Föglö, Åland, Finland, where an exceptionally dense bloom, consisting of 6.3×10⁶cellsL^-1, was observed. Cell concentrations, and copy numbers of both of the target genes were positively correlated with total STX, GTX2, and GTX3 concentrations in the environment, the cell density predicting toxin concentrations with the best accuracy (Spearman's ρ=0.93, p<0.01). The results indicated that all A. ostenfeldii cells in the blooms harbored the genetic capability of PST production, making the detection of sxtA4 a good indicator of toxicity.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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