Systems-level analysis of metabolic mechanism following nitrogen limitation in benthic dinoflagellate Prorocentrum lima

The characteristics of PtHSP40 gene family in Phaeodactylum tricornutum and its response to environmental stresses

Diatom has evolved response mechanisms to cope with multiple environmental stresses. Heat shock protein 40 (HSP40) plays a key role in these response mechanisms. HSP40 gene family in higher plants has been well-studied. However, the HSP40 gene family has not been systematically investigated in marine diatom. In this study, the bioinformatic characteristics, phylogenetic relationship, conserved motifs, gene structure, chromosome distribution and the transcriptional response of PtHSP40 to different environmental stresses were analyzed in the diatom Phaeodactylum tricornutum, and quantitative real-time PCR was conducted. Totally, 55 putative PtHSP40 genes are distributed to 21 chromosomes. All PtHSP40 proteins can be divided into four groups based on their evolutionary relationship, and 54 of them contain a conserved HPD (histidine-proline-aspartic acid tripeptide) motif. Additionally, six, eleven, ten and four PtHSP40 genes were significantly upregulated under the treatments of nitrogen starvation, phosphorus deprivation, 2,2',4,4'-tetrabrominated biphenyl ether (BDE-47) and ocean acidification, respectively. More interestingly, the expression level of 9 PtHSP40 genes was obviously upregulated in response to nickel stress, suggesting the sensitive to metal stress. The different expression models of PtHSP40 genes to environmental stresses imply the specificity of PtHSP40 proteins under different stresses. This study provides a systematic understanding of the PtHSP40 gene family in P. tricornutum and a comprehensive cognition in its functions and response mechanisms to environmental stresses.

Physiological and genetic responses of the benthic dinoflagellate Prorocentrum lima to polystyrene microplastics

Microplastics are well known as contaminants in marine environments. With the development of biofilms, most microplastics will eventually sink and deposit in benthic environment. However, little research has been done on benthic toxic dinoflagellates, and the effects of microplastics on benthic dinoflagellates are unknown. Prorocentrum lima is a cosmopolitan toxic benthic dinoflagellate, which can produce a range of polyether metabolites, such as diarrhetic shellfish poisoning (DSP) toxins. In order to explore the impact of microplastics on marine benthic dinoflagellates, in this paper, we studied the effects of polystyrene (PS) on the growth and toxin production of P. lima. The molecular response of P. lima to microplastic stress was analyzed by transcriptomics. We selected 100 nm, 10 μm and 100 μm PS, and set three concentrations of 1 mg L-1, 10 mg L-1 and 100 mg L-1. The results showed that PS exposure had limited effects on cell growth, but increased the OA and extracellular polysaccharide content at high concentrations. After exposure to PS MPs, genes associated with DSP toxins synthesis, carbohydrate synthesis and energy metabolism, such as glycolysis, TCA cycle and pyruvate metabolism, were significantly up-regulated. We speculated that after exposure to microplastics, P. lima may increase the synthesis of DSP toxins and extracellular polysaccharides, improve the level of energy metabolism and gene expression of ABC transporter, thereby protecting algal cells from damage. Our findings provide new insights into the effects of microplastics on toxic benthic dinoflagellates.

Toxin production and transcriptomic response to nitrate concentrations in the toxic dinoflagellate Alexandrium tamarense

The bloom-forming dinoflagellate Alexandrium tamarense is one of the most important producers of paralytic shellfish poisoning toxins. Annually recurrent blooms of this dinoflagellate species is associated with the incremental nitrogen influx, especially excessive nitrate input. However, limited studies have been conducted on the toxin production and underlying molecular regulation mechanisms of A. tamarense under various nitrate (N) conditions. Therefore, toxin production and transcriptomic responses of this species were investigated. The toxin profile of A. tamarense was consistently dominated by the C2-toxins, and the cellular toxicity increased with N concentrations peaking at 9.23 ± 0.03 fmol/cell in the 883 μM N-added group. Under lower N conditions, expressions of two STX-core genes, sxtA and sxtG, were significantly down-regulated, suggesting that N regulated sxt expression and triggered responses related to toxin biosynthesis. Results of this study provided valuable insights into the ecophysiology of A. tamarense, enhancing our understanding of the occurrence of toxification events in natural environments.

Molecular Phylogeny, Morphology, Growth and Toxicity of Three Benthic Dinoflagellates Ostreopsis sp. 9, Prorocentrum lima and Coolia monotis Developing in Strait of Gibraltar, Southwestern Mediterranean

Few works have been carried out on benthic harmful algal blooms (BHAB) species in the southern Mediterranean and no data are available for the highly dynamic Strait of Gibraltar (western Mediterranean waters). For the first time, Ostreopsis sp. 9, Prorocentrum lima and Coolia monotis were isolated in this key region in terms of exchanges between the Atlantic Ocean and the Mediterranean and subject to intense maritime traffic. Ribotyping confirmed the morphological identification of these three dinoflagellates species. Monoclonal cultures were established and the maximum growth rate and cell yield were measured at a temperature of 24 °C and an irradiance of 90 µmol photons m-2 s-1, for each species: 0.26 ± 0.02 d-1 (8.75 × 103 cell mL-1 after 28 days) for Ostreopsis sp. 9, 0.21 ± 0.01 d-1 (49 × 103 cell mL-1 after 145 days) for P. lima and 0.21 ± 0.01 d-1 (10.02 × 103 cell mL-1 after 28 days) for C. monotis. Only P. lima was toxic with concentrations of okadaic acid and dinophysistoxin-1 measured in optimal growth conditions ranging from 6.4 pg cell-1 to 26.97 pg cell-1 and from 5.19 to 25.27 pg cell-1, respectively. The toxin content of this species varied in function of the growth phase. Temperature influenced the growth and toxin content of P. lima. Results suggest that future warming of Mediterranean coastal waters may lead to higher growth rates and to increases in cellular toxin levels in P. lima. Nitrate and ammonia affected the toxin content of P. lima but no clear trend was noted. In further studies, we have to isolate other BHAB species and strains from Strait of Gibraltar waters to obtain more insight into their diversity and toxicity.

A study on the mechanism of the impact of phenthoate exposure on Prorocentrum lima

Extensive application of organophosphorus pesticides such as phenthoate results in its abundance in ecosystems, particularly in waterbodies, thereby providing the impetus to assess its role in aquatic organisms. However, the impact of phenthoate on marine algal physiological and proteomic response is yet to be explored despite its biological significance. In this study, we thus ought to investigate the impact of phenthoate in the marine dinoflagellate Prorocentrum lima, which is known for synthesizing okadaic acid (OA), the toxin responsible for diarrhetic shellfish poisoning (DSP). Our results showed that P. lima effectively absorbed phenthoate in seawater, with a reduction efficiency of 90.31% after 48 h. Surprisingly, the provision of phenthoate (100 and 1000 µg/L) substantially reduced the OA content of P. lima by 35.08% and 60.28% after 48 h, respectively. Meanwhile, phenthoate treatment significantly reduced the oxidative stress in P. lima. Proteomic analysis revealed that the expression level of seven crucial proteins involved in endocytosis was upregulated, suggesting that P. lima could absorb phenthoate via the endocytic signaling pathway. Importantly, phenthoate treatment resulted in the downregulation of proteins such as polyketide synthase (PKS)- 2, Cytochrome P450 (CYP450)- 1, and CYP450-2, involved in OA synthesis, thereby decreasing the OA biosynthesis by P. lima. Our results demonstrated the potential role of P. lima in the removal of phenthoate in water and exemplified the crucial proteins and their possible molecular mechanisms underpinning the phenthoate remediation by P. lima and also the regulatory role of phenthoate in restricting the OA metabolism. Collectively, these findings uncovered the synergistic mechanisms of phenthoate and P. lima in remediating phenthoate and reducing the toxic impact of P. lima.

Nutrient Deficiencies Impact on the Cellular and Metabolic Responses of Saxitoxin Producing Alexandrium minutum: A Transcriptomic Perspective

Dinoflagellate Alexandrium minutum Halim is commonly associated with harmful algal blooms (HABs) in tropical marine waters due to its saxitoxin production. However, limited information is available regarding the cellular and metabolic changes of A. minutum in nutrient-deficient environments. To fill this gap, our study aimed to investigate the transcriptomic responses of A. minutum under nitrogen and phosphorus deficiency. The induction of nitrogen and phosphorus deficiency resulted in the identification of 1049 and 763 differently expressed genes (DEGs), respectively. Further analysis using gene set enrichment analysis (GSEA) revealed 702 and 1251 enriched gene ontology (GO) terms associated with nitrogen and phosphorus deficiency, respectively. Our results indicate that in laboratory cultures, nitrogen deficiency primarily affects meiosis, carbohydrate catabolism, ammonium assimilation, ion homeostasis, and protein kinase activity. On the other hand, phosphorus deficiency primarily affects the carbon metabolic response, cellular ion transfer, actin-dependent cell movement, signalling pathways, and protein recycling. Our study provides valuable insights into biological processes and genes regulating A. minutum's response to nutrient deficiencies, furthering our understanding of the ecophysiological response of HABs to environmental change.

Transcriptomic Analysis of the Response of the Toxic Dinoflagellate Prorocentrum lima to Phosphorous Limitation

Some dinoflagellates cause harmful algal blooms, releasing toxic secondary metabolites, to the detriment of marine ecosystems and human health. Phosphorus (P) is a limiting macronutrient for dinoflagellate growth in the ocean. Previous studies have been focused on the physiological response of dinoflagellates to ambient P changes. However, the whole-genome's molecular mechanisms are poorly understood. In this study, RNA-Seq was utilized to compare the global gene expression patterns of a marine diarrheic shellfish poisoning (DSP) toxin-producing dinoflagellate, Prorocentrum lima, grown in inorganic P-replete and P-deficient conditions. A total of 148 unigenes were significantly up-regulated, and 30 unigenes were down-regulated under 1/4 P-limited conditions, while 2708 unigenes were significantly up-regulated, and 284 unigenes were down-regulated under 1/16 P-limited conditions. KEGG enrichment analysis of the differentially expressed genes shows that genes related to ribosomal proteins, glycolysis, fatty acid biosynthesis, phagosome formation, and ubiquitin-mediated proteolysis are found to be up-regulated, while most of the genes related to photosynthesis are down-regulated. Further analysis shows that genes encoding P transporters, organic P utilization, and endocytosis are significantly up-regulated in the P-limited cells, indicating a strong ability of P. lima to utilize dissolved inorganic P as well as intracellular organic P. These transcriptomic data are further corroborated by biochemical and physiological analyses, which reveals that under P deficiency, cellular contents of starch, lipid, and toxin increase, while photosynthetic efficiency declines. Our results indicate that has P. lima evolved diverse strategies to acclimatize to low P environments. The accumulation of carbon sources and DSP toxins could provide protection for P. lima to cope with adverse environmental conditions.

Phosphorus deficiency induces sexual reproduction in the dinoflagellate Prorocentrum cordatum

Nitrogen (N) and phosphorus (P) are essential elements whose availability promotes successful growth of phytoplankton and governs aquatic primary productivity. In this study, we investigated the effect of N and/or P deficiency on the sexual reproduction of Prorocentrum cordatum, the dinoflagellate with the haplontic life cycle which causes harmful algal blooms worldwide. In P. cordatum cultures, N and the combined N and P deficiency led to the arrest of the cell cycle in the G0/G1 phases and attenuation of cell culture growth. We observed, that P, but not N deficiency triggered the transition in the life cycle of P. cordatum from vegetative to the sexual stage. This resulted in a sharp increase in percentage of cells with relative nuclear DNA content 2C (zygotes) and the appearance of cells with relative nuclear DNA content 4C (dividing zygotes). Subsequent supplementation with phosphate stimulated meiosis and led to a noticeable increase in the 4C cell number (dividing zygotes). Additionally, we performed transcriptomic data analysis and identified putative phosphate transporters and enzymes involved in the phosphate uptake and regulation of its metabolism by P. cordatum. These include high- and low-affinity inorganic phosphate transporters, atypical alkaline phosphatase, purple acid phosphatases and SPX domain-containing proteins.

Adaptive responses of geographically distinct strains of the benthic dinoflagellate, Prorocentrum lima (Dinophyceae), to varying light intensity and photoperiod

The toxic Prorocentrum lima complex can potentially cause serious harm to the benthos and entire food chain. Studies have revealed physiological differences in strains from different regions related to local environment, while differences in the adaptive responses of P. lima complex should be urgently assessed. Hence, this study explored the adaptive responses to varying light intensities and photoperiods of two P. lima complex strains SHG101 and 3XS34, isolated from the Bohai Sea and the South China Sea, respectively. We found the highest cell density of 7.49 × 104 cells mL-1 recorded in the 3XS strain in the stationary phase with high light intensity exposure. No significant difference was observed in growth rate among SHG groups, however, significant differences were found among 3XS groups ranging from 0.176 to 0.311 d-1. Three key pigments Chl a, Peri, and Fuco accounted for up to 60% of the total pigments. Production and concentrations of pigments and Fv/Fm values exhibit a significant negative correlation with high light intensity and growth. Conversely, total diarrhetic shellfish toxin content and the proportion of diol esters increased to varying degrees after high intensity light exposure, with 3XS strain under high light intensity and a photoperiod of light and darkness (12L:12D) consistently exhibiting the highest levels, finally reaching a maximum (21.6 pg cell-1) at day 28. A shortened photoperiod of high light intensity (8L:16D) resulted in impaired recovery compared with 12L:12D. Furthermore, 3XS showed more delayed and intense adaptive responses, indicating a stronger tolerance compared to SHG. Collectively, these results directly characterized variation in the adaptive responses of geographically distinct strains of P. lima complex, highlighting the previously ignored potential risk diversity of this species.

Transcriptomic analysis of polyketide synthesis in dinoflagellate, Prorocentrum lima

The benthic dinoflagellate Prorocentrum lima is among the most common toxic morphospecies with a cosmopolitan distribution. P. lima can produce polyketide compounds, such as okadaic acid (OA), dinophysistoxin (DTX) and their analogues, which are responsible for diarrhetic shellfish poisoning (DSP). Studying the molecular mechanism of DSP toxin biosynthesis is crucial for understanding the environmental driver influencing toxin biosynthesis as well as for better monitoring of marine ecosystems. Commonly, polyketides are produced by polyketide synthases (PKS). However, no gene has been confirmatively assigned to DSP toxin production. Here, we assembled a transcriptome from 94,730,858 Illumina RNAseq reads using Trinity, resulting in 147,527 unigenes with average sequence length of 1035 nt. Using bioinformatics analysis methods, we found 210 unigenes encoding single-domain PKS with sequence similarity to type I PKSs, as reported in other dinoflagellates. In addition, 15 transcripts encoding multi-domain PKS (forming typical type I PKSs modules) and 5 transcripts encoding hybrid nonribosomal peptide synthetase (NRPS)/PKS were found. Using comparative transcriptome and differential expression analysis, a total of 16 PKS genes were identified to be up-regulated in phosphorus-limited cultures, which was related to the up regulation of toxin expression. In concert with other recent transcriptome analyses, this study contributes to the building consensus that dinoflagellates may utilize a combination of Type I multi-domain and single-domain PKS proteins, in an as yet undefined manner, to synthesize polyketides. Our study provides valuable genomic resource for future research in order to understand the complex mechanism of toxin production in this dinoflagellate.

The responses of harmful dinoflagellate Karenia mikimotoi to simulated ocean acidification at the transcriptional level

The HAB-forming, toxic dinoflagellate Karenia mikimotoi, previously found to benefit from ocean acidification (OA), was cultivated to investigate its transcriptional response to simulated OA for 30 generations. Batch cultures were grown under two CO₂ concentrations, 450 (control) and 1100 (simulated OA) μatm, and physiological parameters [growth, pigments, catalase (CAT), glutathione reductase (GR), and superoxide dismutase (SOD) activity], as well as transcriptomes (obtained via RNA-seq), were compared. Chlorophyll a (Chl a) and carotenoid (Caro) contents, as well as CAT and GR activities, were significantly increased under OA conditions. Transcriptomic analysis revealed 2,490 differentially expressed unigenes in response to OA, which comprised 1.54% of all unigenes. A total of 1,121 unigenes were upregulated, and 1,369 unigenes were downregulated in OA compared to control conditions. The downregulated expression of bicarbonate transporter and carbonic anhydrase genes was a landmark of OA acclimation. Key genes involved in energy metabolism, e.g., photosynthesis, tricarboxylic acid cycle, oxidative phosphorylation, and nitrogen metabolism, were highly upregulated under OA, contributing to increases in the Chl a (55.05%) and Caro (28.37%). The enhanced antioxidant enzyme activities (i.e. CAT, GR) and upregulated genes (i.e. glutathione peroxidase, ascorbate peroxidase, heat shock protein, 20S proteasome, aldehyde dehydrogenase, and apolipoprotein) benefit cells against the potential lower pH stress condition under OA. In addition, the downregulation of four genes associated with motility suggested that the preserved energy could further boost growth. In conclusion, the present study suggests that K. mikimotoi exhibits efficient gene expression regulation for the utilization of energy and resistance to OA-induced stress. Taken together, K. mikimotoi appeared as a tolerant species in response to OA. Thus, more extensive algal blooms that threaten marine organisms are likely in the future. These findings expand current knowledge on the gene expression of HAB-forming species in response to future OA.

Allelopathy of Alexandrium pacificum on Thalassiosira pseudonana in laboratory cultures

Alexandrium pacificum is a toxin-producing dinoflagellate with allelopathic effects. The elucidation of allelopathic mechanism of A. pacificum is of great significance for understanding A. pacificum blooms. To this end, using the model diatom Thalassiosira pseudonana as a target species, we observed changes in physiological, biochemical and gene transcription of T. pseudonana upon being co-cultured with A. pacificum. We found reciprocal effects between A. pacificum and T. pseudonana, and corroborated A. pacificum's allelopathy on T. pseudonana by observing inhibitory effects of filtrate from A. pacificum culture on the growth of T. pseudonana. We also found that co-culturing with A. pacificum, the expression of T. pseudonana genes related to photosynthesis, oxidative phosphorylation, antioxidant system, nutrient absorption and energy metabolism were drastically influenced. Coupled with the alterations in Fv/Fm (the variable/maximum fluorescence ratio), activity of superoxide dismutase, contents of malondialdehyde, neutral lipid and total protein in T. pseudonana co-cultured with A. pacificum, we propose that A. pacificum allelopathy could reduce the efficiency of photosynthesis and energy metabolism of T. pseudonana and caused the oxidative stress, while the nutrient absorption was also affected by allelopathic effects. The resultant data potentially uncovered the allelopathic molecular mechanism of A. pacificum to model alga T. pseudonana. The changes in nutrient uptake and even energy metabolism in T. pseudonana, as an adaptation to environmental conditions, may prevent it from stress-related injuries. Our finding might advance the understanding of allelopathic mechanism of A. pacificum.

Transcriptome responses of the dinoflagellate Karenia mikimotoi driven by nitrogen deficiency

The availability of ambient N nutrient is often correlated with the occurrences of harmful algal bloom formed by certain dinoflagellates, making it important to understand how these species might be responding to such conditions. Here, transcriptome sequencing of Karenia mikimotoi was conducted to understand the underlying molecular mechanisms by which this dinoflagellate copes with nitrogen (N) deficiency. Transcriptomic analysis revealed 8802 unigenes (3.56%) that were differentially expressed with ≥ 2-fold change. Under N-depleted conditions, genes involved in glycolysis, fatty acid metabolism, and the tricarboxylic acid (TCA) cycle as well as lipid accumulation were significantly upregulated. The elevated expression of enzymes used in protein degradation and turnover suggests possible metabolic reconfiguration towards accelerated N recycling. Moreover, a significant increase in urea transporter was observed, indicating increased assimilation of organic nitrogen resources as an alternative in N-depleted cultures of K. mikimotoi. The down-regulated glutamate synthase genes were also identified under N deficiency, suggesting suppression of primary amino acid synthesis to save N resource. Taken together, results of this study show enhanced multiple N resource acquisition and reuse of multiple N resources constitute a comprehensive strategy to cope with N deficiency in a dinoflagellate.

Classification of marine microalgae using low-resolution Mueller matrix images and convolutional neural network

In this paper, we used a convolutional neural network to study the classification of marine microalgae by using low-resolution Mueller matrix images. Mueller matrix images of 12 species of algae from 5 families were measured by a Mueller matrix microscopy with an LED light source at 514 nm wavelength. The data sets of seven resolution levels were generated by the bicubic interpolation algorithm. We conducted two groups of classification experiments; one group classified the algae into 12 classes according to species category, and the other group classified the algae into 5 classes according to family category. In each group of classification experiments, we compared the classification results of the Mueller matrix images with those of the first element (M11) images. The classification accuracy of Mueller matrix images declines gently with the decrease of image resolution, while the accuracy of M11 images declines sharply. The classification accuracy of Mueller matrix images is higher than that of M11 images at each resolution level. At the lowest resolution level, the accuracy of 12-class classification and 5-class classification of full Mueller matrix images is 29.89% and 35.83% higher than those of M11 images, respectively. In addition, we also found that the polarization information of different species had different contributions to the classification. These results show that the polarization information can greatly improve the classification accuracy of low-resolution microalgal images.

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

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

Omics Analysis for Dinoflagellates Biology Research

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