Molecular insights into a dinoflagellate bloom

Comparative secretome and proteome analysis unveils the response mechanism in the phosphorus utilization of Alexandrium pacificum

Phosphorus (P) acts as a crucial limiting nutrient for the growth of marine phytoplankton cells and the formation of algal blooms. The dinoflagellate Alexandrium pacificum is known for causing frequent and intense blooms in specific estuarine and coastal regions. In this study, we investigated the growth and physiological transformations under conditions characterized by P-deficiency, NaH2PO4, and ATP. For the first time, an integrated comparative analysis of the secretome and proteome was performed to investigate the global protein expression profile of A. pacificum, with 355 and 2308 differentially expressed proteins (DEPs), respectively. The results demonstrated that P-deficiency led to a reduction in growth and notable decreases in metabolic processes in A. pacificum. In P-deficient and ATP groups, the expression of secretory protein alkaline phosphatase A (PhoA) was increased, while intracellular acid phosphatase (ACP) displayed significant upregulation in P-deficient group, indicating that A. pacificum has evolved multiple organic P utilization strategies to adapt to low-P environments. A. pacificum can utilize the intracellular carbohydrate storage pools via glycolysis and the TCA cycle to replenish Calvin cycle intermediates. However, the growth of the ATP and NaH2PO4 groups showed no significant alteration. These results suggest that A. pacificum possesses distinct adaptive strategies towards P-deficiency in the environment and employs specific mechanisms for utilizing organic P, which may be a crucial factor in the formation of blooms in low inorganic P environments.

Metaproteomics reveals metabolic activities potentially involved in bloom formation and succession during a mixed dinoflagellate bloom of Prorocentrum obtusidens and Karenia mikimotoi

Understanding metabolic activities involved in bloom formation during a single-species algal bloom has improved greatly. However, little is known about metabolic activities during a multi-species algal bloom. Here, we investigated protein expression profiles at different bloom stages of a mixed dinoflagellate bloom caused by Karenia mikimotoi and Prorocentrum obtusidens (syn. Prorocentrum donghaiense) using a metaproteomic approach. Our results indicated that both P. obtusidens and K. mikimotoi cells highly expressed proteins associated with essential cellular metabolisms such as cell growth and nutrient acquisition before their respective bloom occurrence. P. obtusidens preferentially enhanced uptake and utilization of ammonium, amino acid and organophosphorus-like phospholipid at the early bloom stage, and expressed highly abundant chloroplast peridinin-chlorophyll a-binding protein at the early and the P. obtusidens-dominated bloom stages, indicating their important roles in preferential occurrence and maintenance of P. obtusidens bloom. While absorption and utilization of nutrients, especially ammonium, urea, cyanate, phospholipid, and nucleotide, as well as endocytosis, in K. mikimotoi cells, were enhanced. Notably, both species increased photosynthesis, energy generation, cell proliferation, cell motility and cell defense before their respective blooms, which were beneficial to dealing with adverse external stresses, enabling them to be more competitive and advantageous in complex environments. Interestingly, diatom groups (Skeletonema, Pseudo-nitzschia, and Thalassiosira) decreased uptake and utilization of ambient nutrients and cell proliferation during the bloom period. This study demonstrates that niche differentiation and functional complementarity among phytoplankton species regulate bloom formation and succession during the mixed bloom.

Broad active metabolic pathways, autophagy, and antagonistic hormones regulate dinoflagellate cyst dormancy in marine sediments

This work aimed to reveal the molecular machinery regulating the dormancy of dinoflagellate resting cysts buried in marine sediments. Dinoflagellates play pivotal roles in marine ecosystems, particularly as major contributors of harmful algal blooms. Despite vital roles of cysts in blooming cycles and dinoflagellate ecology, the molecular processes controlling cyst dormancy have largely remained unexplored due to technological difficulties. Using DinoSL as a dinoflagellates-specific mRNA "hook" and SMRT sequencing, we analyzed metatranscriptomes of sediment-buried dinoflagellate cyst assemblages. The data show that most major metabolic and regulatory pathways, except photosynthesis, were transcriptionally active. This suggests the crucial importance of broad metabolic pathways in sustaining cyst viability and germination potential. Further expression analyses of 11 genes (relevant to autophagy and phytohormone gibberellin), lysosome/autolysosome staining, and germination experiments revealed vital roles of autophagy in energy generation, nutrient recycling, and of phytohormones abscisic acid/gibberellin in modulating dormancy/germination of resting cysts. Our findings lay a cornerstone for elucidating the molecular machinery regulating dinoflagellate cyst dormancy.

Assessment of dinoflagellate diversity using DNA metabarcoding reveals toxic dinoflagellate species in Australian coastal waters

Harmful algal blooms (HABs) present severe risks to marine ecosystems, wildlife, human health and economies globally. This study investigates the diversity and abundance of the primary HAB group-dinoflagellates-in Hervey Bay, Queensland, Australia, a region notable for its ecological uniqueness and high susceptibility to HABs. By employing DNA metabarcoding targeting the 18S rRNA V8-V9 regions, we aimed to provide a comprehensive overview of dinoflagellate communities across various substrates and locations, identifying both toxic and non-toxic species. Our findings revealed 66 distinct genera, including known toxin producers such as Alexandrium, Gambierdiscus, Karenia and Prorocentrum, with some toxic species detected in Australian waters for the first time. Additionally, we assessed the abundance of these toxic species and examined the influence of environmental factors on their occurrence. This study emphasises the importance of ongoing monitoring and ecological assessments to manage HAB impacts in vulnerable coastal areas such as Hervey Bay.

Divergent responses of an armored and an unarmored dinoflagellate to ocean acidification

Dinoflagellates, both armored and unarmored, with distinct cell wall difference, are being affected by elevated CO2-induced ocean acidification (OA). However, their specific responses to OA are not well understood. In this study, we investigated the physiological and molecular response of the armored species Prorocentrum obtusidens and the unarmored species Karenia mikimotoi to OA over a 28-day period. The results show that the two species responded differently to OA. Cell growth rate, particulate organic carbon (POC) content, and the activities of C4 pathway enzymes decreased in P. obtusidens under future acidified ocean condition (pH 7.8, 1000 μatm pCO2), but the activities of carbonic anhydrase (CA), ribulose-1,5-bisphosphate carboxylase/oxygenase (RubisCO), and superoxide dismutase (SOD) increased. Whereas cell growth rate, contents of Chl a and PON, and SOD activity altered insignificantly in K. mikimotoi, but contents of POC and total carbohydrate, and the activity of RubisCO increased while the activities of CA and C4 pathway enzymes decreased. Transcriptomic analysis indicates that genes associated with antioxidative response, heat shock protein, proteasome, signal transduction, ribosome, and pH regulation were up-regulated in P. obtusidens but down-regulated in K. mikimotoi. Notably, the synthesis of soluble organic matter (i.e., spermidine and trehalose) was enhanced in K. mikimotoi, thereby regulating intracellular pH and improving stress resistance. This study highlights the divergent response of the armored and unarmored dinoflagellates to OA, with the unarmored dinoflagellate exhibiting a higher ability to withstand this stressor. Therefore, caution should be exercised when predicting the behavior and the eventual fate of dinoflagellates in the future acidified ocean.

Ciguatera poisoning: A review of the ecology and detection methods for Gambierdiscus and Fukuyoa species

Ciguatera poisoning is the most prevalent non-bacterial seafood illness globally, with an estimated 10,000 to 50,000 human cases reported annually. While most symptoms are generally mild, some cases can result in severe and long-lasting neurological and psychological damage, and in some instances, even death. The known causative agents of ciguatera poisoning are benthic toxic dinoflagellate species belonging to the genera Gambierdiscus and Fukuyoa. These species produce highly potent ciguatoxins that bioaccumulate through the marine food chain, eventually reaching humans through seafood consumption. Although Gambierdiscus and Fukuyoa species are widespread in tropical waters worldwide, the full extent of their distribution remains uncertain. This review provides a detailed examination of the ecological dynamics of these dinoflagellates and explores the diverse range of detection methods used to monitor them. These include a focus on molecular techniques for detection, alongside morphological methods, emerging technologies, and a toxin detection overview. Additionally, we offer recommendations on how the field can advance, highlighting novel solutions and next steps for improving detection and monitoring practices. By assessing the strengths and limitations of current approaches and proposing directions for future research, this review aims to support efforts in better understanding and mitigating the risk of ciguatera poisoning.

The micro-ecological feature of colonies is a potential strategy for Phaeocystis globosa bloom formation

Phaeocystis globosa is among the dominant microalgae associated with harmful algal blooms. P. globosa has a polymorphic life cycle and its ecological success has been attributed to algal colony formation, however, few studies have assessed differences in microbial communities and their functional profiles between intra- and extra-colonies during P. globosa blooms. To address this, environmental and metagenomics tools were used to conduct a time-series analysis of the bacterial composition and metabolic characteristics of intra- and extra-colonies during a natural P. globosa bloom. The results show that bacterial composition, biodiversity, and network interactions differed significantly between intra- and extra-colonies. Dominant extra-colonial bacteria were Bacteroidia and Saccharimonadis, while dominant intra-colonial bacteria included Alphaproteobacteria and Gammaproteobacteria. Despite the lower richness and diversity observed in the intra-colonial bacterial community, relative to extra-colonies, the complexity and interconnectedness of the intra-colonial networks were higher. Regarding bacterial function, more functional genes were enriched in substance metabolism (polysaccharides, iron element and dimethylsulfoniopropionate) and signal communication (quorum sensing, indoleacetic acid-IAA) pathways in intra- than in extra-colonies. Conceptual model construction showed that microbial cooperative synthesis of ammonium, vitamin B12, IAA, and siderophores were strongly related to the P. globosa bloom, particularly in the intra-colonial environment. Overall, our data highlight the differences in bacterial structure and functions within and outside the colony during P. globosa blooms. These findings represent fundamental information indicating that phenotypic heterogeneity is a selective strategy that improves microbial population competitiveness and environmental adaptation, benefiting P. globosa bloom formation and persistence.

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.

Electrophysiological techniques in marine microalgae study: A new perspective for harmful algal bloom (HAB) research

Electrophysiological techniques, by measuring bioelectrical signals and ion channel activities in tissues and cells, are now widely utilized to study ion channel-related physiological functions and their underlying mechanisms. Electrophysiological techniques have been extensively employed in the investigation of animals, plants, and microorganisms; however, their application in marine algae lags behind that in other organisms. In this paper, we present an overview of current electrophysiological techniques applicable to algae while reviewing the historical usage of such techniques in this field. Furthermore, we explore the potential specific applications of electrophysiological technology in harmful algal bloom (HAB) research. The application prospects in the studies of stress tolerance, competitive advantage, nutrient absorption, toxin synthesis and secretion by HAB microalgae are discussed and anticipated herein with the aim of providing novel perspectives on HAB investigations.

Molecular physiology of Antarctic diatom natural assemblages and bloom event reveal insights into strategies contributing to their ecological success

The continental shelf of the Western Antarctic Peninsula (WAP) is a highly variable system characterized by strong cross-shelf gradients, rapid regional change, and large blooms of phytoplankton, notably diatoms. Rapid environmental changes coincide with shifts in plankton community composition and productivity, food web dynamics, and biogeochemistry. Despite the progress in identifying important environmental factors influencing plankton community composition in the WAP, the molecular basis for their survival in this oceanic region, as well as variations in species abundance, metabolism, and distribution, remains largely unresolved. Across a gradient of physicochemical parameters, we analyzed the metabolic profiles of phytoplankton as assessed through metatranscriptomic sequencing. Distinct phytoplankton communities and metabolisms closely mirrored the strong gradients in oceanographic parameters that existed from coastal to offshore regions. Diatoms were abundant in coastal, southern regions, where colder and fresher waters were conducive to a bloom of the centric diatom, Actinocyclus. Members of this genus invested heavily in growth and energy production; carbohydrate, amino acid, and nucleotide biosynthesis pathways; and coping with oxidative stress, resulting in uniquely expressed metabolic profiles compared to other diatoms. We observed strong molecular evidence for iron limitation in shelf and slope regions of the WAP, where diatoms in these regions employed iron-starvation induced proteins, a geranylgeranyl reductase, aquaporins, and urease, among other strategies, while limiting the use of iron-containing proteins. The metatranscriptomic survey performed here reveals functional differences in diatom communities and provides further insight into the environmental factors influencing the growth of diatoms and their predicted response to changes in ocean conditions.

IMPORTANCE

In the Southern Ocean, phytoplankton must cope with harsh environmental conditions such as low light and growth-limiting concentrations of the micronutrient iron. Using metratranscriptomics, we assessed the influence of oceanographic variables on the diversity of the phytoplankton community composition and on the metabolic strategies of diatoms along the Western Antarctic Peninsula, a region undergoing rapid climate change. We found that cross-shelf differences in oceanographic parameters such as temperature and variable nutrient concentrations account for most of the differences in phytoplankton community composition and metabolism. We opportunistically characterized the metabolic underpinnings of a large bloom of the centric diatom Actinocyclus in coastal waters of the WAP. Our results indicate that physicochemical differences from onshore to offshore are stronger than between southern and northern regions of the WAP; however, these trends could change in the future, resulting in poleward shifts in functional differences in diatom communities and phytoplankton blooms.

Eukaryotic genomes from a global metagenomic data set illuminate trophic modes and biogeography of ocean plankton

IMPORTANCE

Single-celled eukaryotes play ecologically significant roles in the marine environment, yet fundamental questions about their biodiversity, ecological function, and interactions remain. Environmental sequencing enables researchers to document naturally occurring protistan communities, without culturing bias, yet metagenomic and metatranscriptomic sequencing approaches cannot separate individual species from communities. To more completely capture the genomic content of mixed protistan populations, we can create bins of sequences that represent the same organism (metagenome-assembled genomes [MAGs]). We developed the EukHeist pipeline, which automates the binning of population-level eukaryotic and prokaryotic genomes from metagenomic reads. We show exciting insight into what protistan communities are present and their trophic roles in the ocean. Scalable computational tools, like EukHeist, may accelerate the identification of meaningful genetic signatures from large data sets and complement researchers' efforts to leverage MAG databases for addressing ecological questions, resolving evolutionary relationships, and discovering potentially novel biodiversity.

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.

Systematic exploration of transcriptional responses of interspecies interaction between Karenia mikimotoi and Prorocentrum shikokuense

Karenia mikimotoi and Prorocentrum shikokuense (also identified as P. donghaiense Lu and P. obtusidens Schiller) are two important harmful algal species which often form blooms in the coasts of China. Studies have shown that the allelopathy of K. mikimotoi and P. shikokuense plays an important role in inter-algal competition, though the underlying mechanisms remain largely unclear. Here, we observed reciprocal inhibitory effects between K. mikimotoi and P. shikokuense under co-cultures. Based on the reference sequences, we isolated RNA sequencing reads of K. mikimotoi and P. shikokuense from co-culture metatranscriptome, respectively. We found the genes involved in photosynthesis, carbon fixation, energy metabolism, nutrients absorption and assimilation were significantly up-regulated in K. mikimotoi after co-cultured with P. shikokuense. However, genes involved in DNA replication and cell cycle were significantly down-regulated. These results suggested that co-culture with P. shikokuense stimulated cell metabolism and nutrients competition activity of K. mikimotoi, and inhibited cell cycle. In contrast, genes involved in energy metabolism, cell cycle and nutrients uptake and assimilation were dramatically down-regulated in P. shikokuense under co-culture with K. mikimotoi, indicating that K. mikimotoi could highly affect the cellular activity of P. shikokuense. In addition, the expression of PLA2G12 (Group XII secretory phospholipase A2) that can catalyze the accumulation of linoleic acid or linolenic acid, and nitrate reductase that may be involved in nitric oxide production were significantly increased in K. mikimotoi, suggesting that PLA2G12 and nitrate reductase may play important roles in the allelopathy of K. mikimotoi. Our findings shed new light on the interspecies competition between K. mikimotoi and P. shikokuense, and provide a novel strategy for studying interspecific competition in complex systems.

In Situ Molecular Ecological Analyses Illuminate Distinct Factors Regulating Formation and Demise of a Harmful Dinoflagellate Bloom

The development and demise of a harmful algal bloom (HAB) are generally regulated by multiple processes; identifying specific critical drivers for a specific bloom is important yet challenging. Here, we conducted a whole-assemblage molecular ecological study on a dinoflagellate bloom to address the hypothesis that energy and nutrient acquisition, defense against grazing and microbial attacks, and sexual reproduction are critical to the rise and demise of the bloom. Microscopic and molecular analyses identified the bloom-causing species as Karenia longicanalis and showed that the ciliate Strombidinopsis sp. was dominant in a nonbloom plankton community, whereas the diatom Chaetoceros sp. dominated the after-bloom community, along with remarkable shifts in the community structure for both eukaryotes and prokaryotes. Metatranscriptomic analysis indicated that heightened energy and nutrient acquisition in K. longicanalis significantly contributed to bloom development. In contrast, active grazing by the ciliate Strombidinopsis sp. and attacks by algicidal bacteria (Rhodobacteracea, Cryomorphaceae, and Rhodobacteracea) and viruses prevented (at nonbloom stage) or collapsed the bloom (in after-bloom stage). Additionally, nutrition competition by the Chaetoceros diatoms plausibly contributed to bloom demise. The findings suggest the importance of energy and nutrients in promoting this K. longicanalis bloom and the failure of antimicrobial defense and competition of diatoms as the major bloom suppressor and terminator. This study provides novel insights into bloom-regulating mechanisms and the first transcriptomic data set of K. longicanalis, which will be a valuable resource and essential foundation for further elucidation of bloom regulators of this and related species of Kareniaceae in the future. IMPORTANCE HABs have increasingly occurred and impacted human health, aquatic ecosystems, and coastal economies. Despite great efforts, the factors that drive the development and termination of a bloom are poorly understood, largely due to inadequate in situ data about the physiology and metabolism of the causal species and the community. Using an integrative molecular ecological approach, we determined that heightened energy and nutrient acquisition promoted the bloom, while resource allocation in defense and failure to defend against grazing and microbial attacks likely prevented or terminated the bloom. Our findings reveal the differential roles of multiple abiotic and biotic environmental factors in driving the formation or demise of a toxic dinoflagellate bloom, suggesting the importance of a balanced biodiverse ecosystem in preventing a dinoflagellate bloom. The study also demonstrates the power of whole-assemblage metatranscriptomics coupled to DNA barcoding in illuminating plankton ecological processes and the underlying species and functional diversities.

Machine learning predicts the impact of antibiotic properties on the composition and functioning of bacterial community in aquatic habitats

In the past decades, hundreds of antibiotics have been isolated from microbial metabolites or have been artificially synthesized for protecting humans, animals and crops from microbial infections. Their everlasting usage results in impacts on the microbial community composition and causes well-known collateral damage to the functioning of microbial communities. Nevertheless, the impact of different antibiotic properties on aquatic microbial communities have so far only poorly been disentangled. Here we characterized the environmental risk of 50 main kinds of antibiotics from 9 classes at a concentration of 10 μg/L for aquatic bacterial communities via metadata analysis combined with machine learning. Metadata analysis showed that the alpha diversity of the bacterial community increased only after treatment with aminoglycoside and β-lactam antibiotics, while its structure was changed by almost all tested antibiotics. The antibiotic treatment also disturbed the functions of the bacterial community, especially with regard to metabolic pathways, including amino acids, cofactors, vitamins, xenobiotics and carbohydrate metabolism. The critical characteristics (atom stereocenter count, number of hydrogen atoms in the antibiotic, and the adipose water coefficient) of antibiotics affecting the composition of the bacterial community in aquatic habitats were screened by machine learning. The key characteristics of antibiotics affecting the function bacterial communities were the number of hydrogen atoms, molecular weight and complexity. In summary, by developing machine learning models and by performing metadata analysis, this study provides the relationship between the properties of antibiotics and their adverse impacts on aquatic microbial communities from a macro perspective. The study also provides guidance for the rational design of antibiotics.

Metaproteomics reveals unique metabolic niches of dominant bacterial groups in response to rapid regime shifts during a mixed dinoflagellate bloom

The dynamics of bacterial composition and metabolic activity during a distinct phytoplankton bloom have been reported. However, there is limited information on the bacterial community response to drastic environmental changes caused by species succession during a mixed-species bloom. This study investigated active bacterial groups and metabolic activity during a mixed bloom formed by dinoflagellates Prorocentrum obtusidens and Karenia mikimotoi using a metaproteomic approach. Bacterial community structure and dominant bacterial groups varied rapidly with the bloom regime shifts caused by species succession. Pseudoalteromonas and Vibrio dominated the bacterial community in the P. obtusidens-dominated regime, while Alteromonas, Cytophaga-Flavobacteria-Bacteroides (CFB) group, and marine Roseobacter clade (MRC) were the major contributors in other regimes, with the most abundant taxa being Alteromonas in the K. mikimotoi-dominated regime and the CFB group in the dissipation regime. Specific metabolic niches and unique substrate specificity of different bacterial groups enabled them to dominate and thrive in different bloom regimes. High metabolic plasticity in signal response, substrate utilization, motility, and adhesion are essential for bacteria to respond to drastic bloom regime shift, and the predominance of specific bacteria under unique bloom regimes may be the result of long-term coevolution between bacteria and bloom-forming phytoplankton species.

Comparative metagenomics of phytoplankton blooms after nutrient enrichment of oligotrophic marine waters

Increasing anthropogenic pressures on the coastal marine environments impact these ecosystems via a variety of mechanisms including nutrient loading, leading to eutrophication and increases in algal blooms. Here, we use a metagenomics approach to assess the taxonomic and functional changes of the microbial community throughout a nutrient enriched mesocosm phytoplankton bloom. We tested four different nutrient treatments consisting of either nitrate and phosphate or nitrate, phosphate and silicate, administered on the first day or continuously for the first two weeks of the experiment. Our results show a shift in the taxonomic composition of the community over time that is dependent on the nutrient addition regime. Significant differences in the functional potential of the communities were detected, with an interaction between bloom period (pre-bloom, bloom and post-bloom) and nutrient treatment (p = 0.004). A sharp drop in functional similarity was observed in the first week in all treatments and after 20 days had not returned to pre-bloom levels. Changes within energy metabolism pathways showed a remarkable enrichment of the dissimilatory nitrate reduction pathway in the post-bloom period. Eukaryotic oxidative phosphorylation and photosynthetic antenna proteins were more abundant during the bloom, especially in the continuous treatment with silicate. Our results suggest that continuous (i.e. chronic) nutrient enrichment has a larger effect on the functioning of marine systems compared to a single (i.e acute) addition. A deep understanding of the functional and taxonomic shifts in the community during blooms is essential to reverse or mitigate human impacts on coastal environments.

Distribution patterns and seasonal variations in phytoplankton communities of the hypersaline Pulicat lagoon, India

Phytoplankton structure and patterns are key components to forecast the net result of the gain and loss process that outline the resilience of the lagoon ecosystem. In order to understand the phytoplankton community structure and its relationship with the environmental variables in the shallow hypersaline Pulicat lagoon, east coast of India, observations were carried out during August 2018-January 2019 covering the three seasons: premonsoon (PrM), monsoon (M), and postmonsoon (PoM). The salinity of the lagoon varied with a minimum of 12.1 for the M and a maximum of 81.65 during the PoM. The clustering analysis performed on the phytoplankton abundance data separated the lagoon into three sectors: north sector (NS), central sector (CS), and south sector (SS). A total of 59 taxa/morphotypes from four taxonomic classes (Bacillariophyceae, Chlorophyceae, Cyanophyceae, and Dinophyceae) were recorded during the study period. The class Bacillariophyceae was dominant in the lagoon both spatially and temporally by 44.06% with Chaetoceros borealis as dominant species. Presence of characteristic species like Dunaliella sp. was observed in the higher salinity, whereas Pediastrum duplex and Scenedesmus sp. were dominant in the freshwater influx areas. The individual-based functional approach allowed grouping these taxa into 11 functional entities based on the derived 4 functional trait values (cell size, trophic regime, mobility, and coloniality). Formation of algal blooms of Protoperidinium sp. (3.3×10⁵ind L^-1) and Odentella sp. (2.8×10⁵ind L^-1) was observed in the SS during PoM as a result of reduced water exchange in the lagoon. During the same period, toxin-producing strains like Anabaena sp. and Nostoc sp. of Cyanophyceae were also recorded. Correlating the three sectors of the lake (NS, CS, and SS), it is observed that the physical, chemical, and biological properties of the lake varied continuously depending on the season and freshwater availability. Seasonal nutrient stoichiometry played a significant role in regulating the community structure and distribution pattern of phytoplankton communities of the lagoon.

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.

Microbial community transcriptional patterns vary in response to mesoscale forcing in the North Pacific Subtropical Gyre

The physical and biological dynamics that influence phytoplankton communities in the oligotrophic ocean are complex, changing across broad temporal and spatial scales. Eukaryotic phytoplankton (e.g., diatoms), despite their relatively low abundance in oligotrophic waters, are responsible for a large component of the organic matter flux to the ocean interior. Mesoscale eddies can impact both microbial community structure and function, enhancing primary production and carbon export, but the mechanisms that underpin these dynamics are still poorly understood. Here, mesoscale eddy influences on the taxonomic diversity and expressed functional profiles of surface communities of microeukaryotes and particle-associated heterotrophic bacteria from the North Pacific Subtropical Gyre were assessed over 2 years (spring 2016 and summer 2017). The taxonomic diversity of the microeukaryotes significantly differed by eddy polarity (cyclonic versus anticyclonic) and between sampling seasons/years and was significantly correlated with the taxonomic diversity of particle-associated heterotrophic bacteria. The expressed functional profile of these taxonomically distinct microeukaryotes varied consistently as a function of eddy polarity, with cyclones having a different expression pattern than anticyclones, and between sampling seasons/years. These data suggest that mesoscale forcing, and associated changes in biogeochemistry, could drive specific physiological responses in the resident microeukaryote community, independent of species composition.

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.

Metaproteomics reveals the molecular mechanism underlying bloom maintenance of a marine dinoflagellate under low ambient CO2 and inorganic nutrients

Dinoflagellate blooming periods are paradoxically characterized by high biomass growth rate and low ambient dissolved CO₂ and inorganic nutrients, however, the underlying mechanisms linking cell growth and nutrient acquisition are poorly understood. Here, we compared metaproteomes of non-bloom, mid-blooming and late-blooming cells of a marine dinoflagellate Prorocentrum donghaiense. Cell division, metabolism of carbon, nitrogen, phosphorus, lipid, porphyrin and chlorophyll were more active in blooming cells than in non-bloom cells. Up-regulation of carbonic anhydrase, ribulose-1,5-bisphosphate carboxylase/oxygenase II, and C₄-cycle proteins enhanced CO₂ assimilation of P. donghaiense. Proteins participating in external organic nutrient acquisition and conversion, such as transporters for fatty acids, peptides and amino acids, external- and internal-phosphomonoester hydrolase, and diverse peptidases and amino acid transaminases, exhibited higher expression in blooming cells relative to non-bloom cells. Interestingly, dissolved organic nitrogen (DON) such as urea and aspartate significantly down-regulated expression and activity of carbon assimilation proteins except for RuBisCO form II, suggesting that DON provided sufficient carbon source which reduced the need to concentrate internal CO₂. This study demonstrates that coupling of efficient CO₂ assimilation with DON utilization are essential for bloom maintenance of P. donghaiense, and future efforts should be devoted to dissolved organic nutrients for prevention and management of dinoflagelllate blooms.

Initiation of efficient C4 pathway in response to low ambient CO2 during the bloom period of a marine dinoflagellate

Dinoflagellates are important primary producers and major causative agents of harmful algal blooms in the global ocean. Despite the great ecological significance, the photosynthetic carbon acquisition by dinoflagellates is still poorly understood. The pathways of photosynthetic carbon assimilation in a marine dinoflagellate Prorocentrum donghaiense under both in situ and laboratory-simulated bloom conditions were investigated using a combination of metaproteomics, qPCR, stable carbon isotope and targeted metabolomics approaches. A rapid consumption of dissolved CO₂ to generate high biomass was observed as the bloom proceeded. The carbon assimilation genes and proteins including intracellular carbonic anhydrase 2, phosphoenolpyruvate carboxylase, phosphoenolpyruvate carboxykinase and RubisCO as well as their enzyme activities were all highly expressed at the low CO₂ level, indicating that C₄ photosynthetic pathway functioned in the blooming P. donghaiense cells. Furthermore, δ¹³ C values and content of C₄ compound (malate) significantly increased with the decreasing CO₂ concentration. The transition from C₃ to C₄ pathway minimizes the internal CO₂ leakage and guarantees efficient carbon fixation at the low CO₂ level. This study demonstrates the existence of C₄ photosynthetic pathway in a marine dinoflagellate and reveals its important complementary role to assist carbon assimilation for cell proliferation during the bloom period.

Autoactivation of Translation Causes the Bloom of Prorocentrum donghaiense in Harmful Algal Blooms

Harmful algal blooms (HABs) are symptomatic of ecosystem imbalance, leading to major worldwide marine natural disasters, and seriously threaten the human health. Some HAB algae's exceptional genome size prohibited the genomic investigations on molecular mechanisms, for example, Prorocentrum. This study performed translatome sequencing (RNC-seq) for Prorocentrum donghaiense to assemble the translatome reference sequences on appropriate cost to enable the global molecular study at translatome and proteome levels. By analyzing the translatome and proteome of P. donghaiense in phosphor-rich, phosphor-deficient, and phosphor-restored media, we found massive up-regulation of energy and material production pathways in phosphor-rich conditions that enables autoactivation of translation, which is the key to its exponential growth in HABs. To break down the autoactivation, we demonstrated that mild translation delay using very low concentrations of cycloheximide efficiently controls the blooming without harming other aquatic organisms and humans. Our result provides a novel hint for controlling HABs and demonstrated the RNC-seq as an economic strategy on investigating functions of organisms with large and unknown genomes.

Cobalamin and microbial plankton dynamics along a coastal to offshore transect in the Eastern North Atlantic Ocean

Cobalamin (B12) is an essential cofactor that is exclusively synthesized by some prokaryotes while many prokaryotes and eukaryotes require an external supply of B12. The spatial and temporal availability of B12 is poorly understood in marine ecosystems. Field measurements of B12 along with a large set of ancillary biotic and abiotic factors were obtained during three oceanographic cruises in the NW Iberian Peninsula, covering different spatial and temporal scales. B12 concentrations were remarkably low (<1.5 pM) in all samples, being significantly higher at the subsurface Eastern North Atlantic Central Water than at shallower depths, suggesting that B12 supply in this water mass is greater than demand. Multiple regression models excluded B12 concentration as predictive variable for phytoplankton biomass or production, regardless of the presence of B12-requiring algae. Prokaryote production was the best predictor for primary production, and eukaryote community composition was better correlated with prokaryote community composition than with nutritional resources, suggesting that biotic interactions play a significant role in regulating microbial communities. Interestingly, co-occurrence network analyses based on 16S and 18S rRNA sequences allowed the identification of significant associations between potential B12 producers and consumers (e.g. Thaumarchaeota and Dynophyceae, or Amylibacter and Ostreococcus respectively), which can now be investigated using model systems in the laboratory.

Picophytoplankton dynamics in a large temperate estuary and impacts of extreme storm events

Picophytoplankton (PicoP) are increasingly recognized as significant contributors to primary productivity and phytoplankton biomass in coastal and estuarine systems. Remarkably though, PicoP composition is unknown or not well-resolved in several large estuaries including the semi-lagoonal Neuse River Estuary (NRE), a tributary of the second largest estuary-system in the lower USA, the Pamlico-Albemarle Sound. The NRE is impacted by extreme weather events, including recent increases in precipitation and flooding associated with tropical cyclones. Here we examined the impacts of moderate to extreme (Hurricane Florence, September 2018) precipitation events on NRE PicoP abundances and composition using flow cytometry, over a 1.5 year period. Phycocyanin-rich Synechococcus-like cells were the most dominant PicoP, reaching ~ 106 cells mL-1, which highlights their importance as key primary producers in this relatively long residence-time estuary. Ephemeral "blooms" of picoeukaryotic phytoplankton (PEUK) during spring and after spikes in river flow were also detected, making PEUK periodically major contributors to PicoP biomass (up to ~ 80%). About half of the variation in PicoP abundance was explained by measured environmental variables. Temperature explained the most variation (24.5%). Change in total dissolved nitrogen concentration, an indication of increased river discharge, explained the second-most variation in PicoP abundance (15.9%). The short-term impacts of extreme river discharge from Hurricane Florence were particularly evident as PicoP biomass was reduced by ~ 100-fold for more than 2 weeks. We conclude that precipitation is a highly influential factor on estuarine PicoP biomass and composition, and show how 'wetter' future climate conditions will have ecosystem impacts down to the smallest of phytoplankton.

Phytoplankton composition in a eutrophic estuary: Comparison of multiple taxonomic approaches and influence of environmental factors

To assess the comparability between taxonomic identification methods for phytoplankton, multiple approaches were used to characterize phytoplankton community composition within the Neuse River Estuary (NRE), North Carolina. Small subunit 18S rRNA gene sequencing and accessory pigment analysis displayed similar trends, indicating chlorophytes were the dominant microalgal group during most of the year, whereas results from microscopic cell counts, biovolume analysis and metatranscriptomics suggested diatom and dinoflagellate-dominated communities. Spatial environmental gradients drove variation in taxonomic composition due to preferences for specific environmental conditions among different microalgal groups. Cryptophytes were a greater proportion of the phytoplankton community within high nutrient, fresher environments whereas diatoms and dinoflagellates dominated higher salinity sections of the estuary. This study provides a detailed examination of phytoplankton communities associated with environmental gradients present in the NRE. The high level of taxonomic resolution offered by DNA sequencing (i.e., species to sub-species level) provides a better understanding of population dynamics at the base of estuarine food webs.

Temporal heterogeneity of microbial communities and metabolic activities during a natural algal bloom

Elucidating the interactions between algae and associated microbial communities is critical for understanding the mechanisms that mediate the dynamic of harmful algal blooms (HABs) in marine environment. However, the microbial functional profiles and their biogeochemical potential in HABs process remains elusive, especially during a complete natural HAB cycle. Here, we used pyrosequencing and functional gene array (GeoChip) to investigate microbial community dynamics and metabolic potential during a natural dinoflagellate (Noctiluca scintillans) bloom. The results shown that bacterioplankton exhibited significant temporal heterogeneity over the course of the bloom stages. Microbial succession was co-driven by environmental parameters and biotic interactions. The functional analysis revealed significant variations in microbial metabolism during matter cycling. At bloom onset-stage, metabolic potential associated with iron oxidation and transport was elevated. Carbon fixation and degradation, denitrification, phosphorus acquisition, and sulfur transfer/oxidation were significantly enhanced at the plateau stage. During the decline and terminal stages, oxidative stress, lysis of compounds, and toxin degradation & protease synthesis increased. This work reveal phycosphere microorganisms can enhanced organic C decomposition capacity, altered N assimilation rate and S/P turnover efficiency, and balancing of the Fe budget during HAB process. The ecological linkage analysis has further shown that microbial composition and functional potential were significantly linked to algal blooms occurrence. It suggest that structural variability and functional plasticity of microbial communities influence HAB trajectory.

Functional profiles of phycospheric microorganisms during a marine dinoflagellate bloom

Harmful algal blooms (HABs) are an ecological concern but relatively few studies have investigated the functional potential of bacterioplankton over a complete algal bloom cycle, which is critical for determining their contribution to the fate of algal blooms. To address this point, we carried out a time-series metagenomic analysis of the functional features of microbial communities at three different Gymnodinium catenatum bloom stages (pre-, peak-, and post-bloom). Different microbial composition were observed during the blooming stages. The environmental parameters and correlation networks co-contribute to microbial variability, and the former explained 38.4% of total variations of the bacterioplankton community composition. Functionally, a range of pathways involved in carbon, nitrogen, phosphorus and sulfur cycling were significantly different during the various HAB stages. Genes associated with carbohydrate-active enzymes, denitrification, and iron oxidation were enriched at the pre-bloom stage; genes involved in reductive citrate cycle for carbon fixation, carbon degradation, nitrification and phosphate transport were enhanced at the peak stage; and relative gene abundance related to sulfur oxidation, vitamin synthesis, and iron transport and storage was increased at the post-bloom stage. The ecological linkage analysis has shown that microbial functional potential especially the C/P/Fe metabolism were significantly linked to the fate of the algal blooms. Taken together, our results demonstrated that microorganisms displayed successional patterns not only at the community level, but also in the metabolic potential on HAB's progression. This work contributes to a growing understanding of microbial structural elasticity and functional plasticity and shed light on the potential mechanisms of microbial-mediated HAB trajectory.

Insights into the transcriptional responses of a microbial community to silver nanoparticles in a freshwater microcosm

Silver nanoparticles (AgNPs) are widely used because of their excellent antibacterial properties. They are, however, easily discharged into the water environment, causing potential adverse environmental effects. Meta-transcriptomic analyses are helpful to study the transcriptional response of prokaryotic and eukaryotic aquatic microorganisms to AgNPs. In the present study, microcosms were used to investigate the toxicity of AgNPs to a natural aquatic microbial community. It was found that a 7-day exposure to 10 μg L^-1 silver nanoparticles (AgNPs) dramatically affected the structure of the microbial community. Aquatic micro eukaryota (including eukaryotic algae, fungi, and zooplankton) and bacteria (i.e., heterotrophic bacteria and cyanobacteria) responded differently to the AgNPs stress. Meta-transcriptomic analyses demonstrated that eukaryota could use multiple cellular strategies to cope with AgNPs stress, such as enhancing nitrogen and sulfur metabolism, over-expressing genes related to translation, amino acids biosynthesis, and promoting bacterial-eukaryotic algae interactions. By contrast, bacteria were negatively affected by AgNPs with less signs of detoxification than in case of eukaryota; various pathways related to energy metabolism, DNA replication and genetic repair were seriously inhibited by AgNPs. As a result, eukaryotic algae (mainly Chlorophyta) dominated over cyanobacteria in the AgNPs treated microcosms over the 7-d exposure. The present study helps to understand the effects of AgNPs on aquatic microorganisms and provides insights into the contrasting AgNPs toxicity in eukaryota and bacteria.

Dur3 and nrt2 genes in the bloom-forming dinoflagellate Prorocentrum minimum: Transcriptional responses to available nitrogen sources

The increasing inflow of nitrogen (N) substrates into marine nearshore ecosystems induces proliferation of harmful algal blooms (HABs) of dinoflagellates, such as potentially toxic invasive species Prorocentrum minimum. In this study, we estimated the influence of NO3-, NH4+ and urea on transcription levels and urea transporter dur3 and nitrate transporter nrt2 genes expression in these dinoflagellates. We identified dur3 and nrt2 genes sequences in unannotated transcriptomes of P. minimum and other dinoflagellates presented in MMETSP database. Phylogenetic analysis showed that these genes of dinoflagellates clustered to the distinct clade demonstrating evolutionary relationship with the other known dur3 and nrt2 genes of microalgae. The evaluation of expression levels of dur3 and nrt2 genes by RT-qPCR revealed their sensitivity to input of the studied N sources. Dur3 expression levels were downregulated after the supplementation of additional N sources and were 1.7-2.6-fold lower than in the nitrate-grown culture. Nrt2 expression levels decreased 1.9-fold in the presence of NH4+. We estimated total RNA and DNA synthesis rates by the analysis of incorporation of 3H-thymidine and 3H-uridine in batch and continuous cultures. Addition of N compounds did not affect the DNA synthesis rates. Transcription levels increased up to 12.5-fold after the N supplementation in urea-limited treatments. Investigation of various nitrogen sources as biomarkers of dinoflagellate proliferation due to their differentiated impact on expression of dur3 and nrt2 genes and transcription rates in P. minimum cells allowed concluding about high potential of the studied parameters for future modeling of HABs under global N pollution.

Nuclear Gene Transformation in the Dinoflagellate Oxyrrhis marina

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

Comparative metatranscriptomic profiling and microRNA sequencing to reveal active metabolic pathways associated with a dinoflagellate bloom

Harmful algal blooms (HABs) have increased as a result of global climate and environmental changes, exerting increasing impacts on the aquatic ecosystem, coastal economy, and human health. Despite great research efforts, our understanding on the drivers of HABs is still limited in part because HAB species' physiology is difficult to probe in situ. Here, we used molecular ecological analyses to characterize a dinoflagellate bloom at Xiamen Harbor, China. Prorocentrum donghaiense was identified as the culprit, which nutrient bioassays showed were not nutrient-limited. Metatranscriptome profiling revealed that P. donghaiense highly expressed genes related to N- and P-nutrient uptake, phagotrophy, energy metabolism (photosynthesis, oxidative phophorylation, and rhodopsin) and carbohydrate metabolism (glycolysis/gluconeogenesis, TCA cycle and pentose phosphate) during the bloom. Many genes in P. donghaiense were up-regulated at night, including phagotrophy and environmental communication genes, and showed active expression in mitosis. Eight microbial defense genes were up-regulated in the bloom compared with previously analyzed laboratory cultures. Furthermore, 76 P. donghaiense microRNA were identified from the bloom, and their target genes exhibited marked differences in amino acid metabolism between the bloom and cultures and the potential of up-regulated antibiotic and cell communication capabilities. These findings, consistent with and complementary to recent reports, reveal major metabolic processes in P. donghaiense potentially important for bloom formation and provide a gene repertoire for developing bloom markers in future research.

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

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

Transcriptomic response to changing ambient phosphorus in the marine dinoflagellate Prorocentrum donghaiense

Dinoflagellates represent major contributors to the harmful algal blooms in the oceans. Phosphorus (P) is an essential macronutrient that limits the growth and proliferation of dinoflagellates. However, the specific molecular mechanisms involved in the P acclimation of dinoflagellates remain poorly understood. Here, the transcriptomes of a dinoflagellate Prorocentrum donghaiense grown under inorganic P-replete, P-deficient, and inorganic- and organic P-resupplied conditions were compared. Genes encoding low- and high-affinity P transporters were significantly down-regulated in the P-deficient cells, while organic P utilization genes were significantly up-regulated, indicating strong ability of P. donghaiense to utilize organic P. Up-regulation of membrane phospholipid catabolism and endocytosis provided intracellular and extracellular organic P for the P-deficient cells. Physiological responses of P. donghaiense to dissolved inorganic P (DIP) or dissolved organic P (DOP) resupply exhibited insignificant differences. However, the corresponding transcriptomic responses significantly differed. Although the expression of multiple genes was significantly altered after DIP resupplementation, few biological processes varied. In contrast, various metabolic processes associated with cell growth, such as translation, transport, nucleotide, carbohydrate and lipid metabolisms, were significantly altered in the DOP-resupplied cells. Our results indicated that P. donghaiense evolved diverse DOP utilization strategies to adapt to low P environments, and that DOPs might play critical roles in the P. donghaiense bloom formation.

Proteomic analysis of the hepatotoxicity of Microcystis aeruginosa in adult zebrafish (Danio rerio) and its potential mechanisms

Microcystis aeruginosa is one of the main species of cyanobacteria that causes water blooms. M. aeruginosa can release into the water several types of microcystins (MCs), which are harmful to aquatic organisms and even humans. However, few studies have investigated the hepatotoxicity of M. aeruginosa itself in zebrafish in environments that simulate natural aquatic systems. The objective of this study was to evaluate the hepatotoxicity of M. aeruginosa in adult zebrafish (Danio rerio) after short-term (96 h) exposure and to elucidate the potential underlying mechanisms. Distinct histological changes in the liver, such as enlargement of the peripheral nuclei and sinusoids and the appearance of fibroblasts, were observed in zebrafish grown in M. aeruginosa culture. In addition, antioxidant enzyme activity was activated and protein phosphatase (PP) activity was significantly decreased with increasing microalgal density. A proteomic analysis revealed alterations in a number of protein pathways, including ribosome translation, immune response, energy metabolism and oxidative phosphorylation pathways. Western blot and real-time PCR analyses confirmed the results of the proteomic analysis. All results indicated that M. aeruginosa could disrupt hepatic functions in adult zebrafish, thus highlighting the necessity of ecotoxicity assessments for M. aeruginosa at environmentally relevant densities.

Functional Differences in the Blooming Phytoplankton Heterosigma akashiwo and Prorocentrum donghaiense Revealed by Comparative Metaproteomics

Phytoplankton blooms are natural phenomena in the ocean, which are the results of rapid cell growth of some phytoplankton species in a unique environment. However, little is known about the molecular events occurring during the bloom. Here, we compared metaproteomes of two phytoplankton Heterosigma akashiwo and Prorocentrum donghaiense in the coastal East China Sea. H. akashiwo and P. donghaiense accounted for 7.82% and 4.74% of the phytoplankton community protein abundances in the nonbloom sample, whereas they contributed to 60.13% and 78.09%, respectively, in their individual blooming samples. Compared with P. donghaiense, H. akashiwo possessed a significantly higher abundance of light-harvesting complex proteins, carbonic anhydrasem and RuBisCO. The blooming H. akashiwo cells expressed more proteins related to external nutrient acquisition, such as bicarbonate transporter SLC4, ammonium transporter, nitrite transporter, and alkaline phosphatase, while the blooming P. donghaiense cells highly expressed proteins related to extra- and intracellular organic nutrient utilization, such as amino acid transporter, 5'-nucleotidase, acid phosphatase, and tripeptidyl-peptidase. The strong capabilities of light harvesting, as well as acquisition and assimilation of inorganic carbon, nitrogen, and phosphorus, facilitated the formation of the H. akashiwo bloom under the high turbidity and inorganic nutrient-sufficient condition, whereas the competitive advantages in organic nutrient acquisition and reallocation guaranteed the occurrence of the P. donghaiense bloom under the inorganic nutrient-insufficient condition. This study highlights the power of metaproteomics for revealing the underlying molecular behaviors of different coexisting phytoplankton species and advances our knowledge on the formation of phytoplankton blooms.IMPORTANCE A deep understanding of the mechanisms driving bloom formation is a prerequisite for effective bloom management. Metaproteomics was applied in this study to reveal the adaptive and responsive strategies of two coexisting phytoplankton species, H. akashiwo and P. donghaiense, during their bloom periods. Metabolic features and niche divergence in light harvesting, as well as carbon, nitrogen, and phosphorus acquisition and assimilation likely promoted the bloom occurrence under different environments. The molecular behaviors of coexisting bloom-causing species will give clues for bloom monitoring and management in the oceans.

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.

Strategies among phytoplankton in response to alleviation of nutrient stress in a subtropical gyre

Despite generally low primary productivity and diatom abundances in oligotrophic subtropical gyres, the North Atlantic Subtropical Gyre (NASG) exhibits significant diatom-driven carbon export on an annual basis. Subsurface pulses of nutrients likely fuel brief episodes of diatom growth, but the exact mechanisms utilized by diatoms in response to these nutrient injections remain understudied within near-natural settings. Here we simulated delivery of subsurface nutrients and compare the response among eukaryotic phytoplankton using a combination of physiological techniques and metatranscriptomics. We show that eukaryotic phytoplankton groups exhibit differing levels of transcriptional responsiveness and expression of orthologous genes in response to release from nutrient limitation. In particular, strategies for use of newly delivered nutrients are distinct among phytoplankton groups. Diatoms channel new nitrate to growth-related strategies while physiological measurements and gene expression patterns of other groups suggest alternative strategies. The gene expression patterns displayed here provide insights into the cellular mechanisms that underlie diatom subsistence during chronic nitrogen-depleted conditions and growth upon nutrient delivery that can enhance carbon export from the surface ocean.

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.

Evaluation of the taxonomic and functional variation of freshwater plankton communities induced by trace amounts of the antibiotic ciprofloxacin

Ciprofloxacin (CIP), one of the most frequently detected antibiotics in water systems, has become an aquatic contaminant because of improper disposal and excretion by humans and animals. It is still unknown how trace amounts of CIP affect the aquatic microbial community diversity and function. We therefore investigated the effects of CIP on the structure and function of freshwater microbial communities via 16S/18S rRNA gene sequencing and metatranscriptomic analyses. CIP treatment (7 μg/L) did not significantly alter the physical and chemical condition of the water body as well as the composition of the main species in the community, but slightly increased the relative abundance of cyanobacteria and decreased the relative abundance of eukaryotes. Metatranscriptomic results showed that bacteria enhanced their phosphorus transport and photosynthesis after CIP exposure. The replication, transcription, translation and cell proliferation were all suppressed in eukaryotes, while the bacteria were not affected in any of these aspects. This interesting phenomenon was the exact opposite to both the antibacterial property of CIP and its safety for eukaryotes. We hypothesize that reciprocal and antagonistic interactions in the microcosm both contribute to this result: cyanobacteria may enhance their tolerance to CIP through benefiting from cross-feeding and some secreted substances that withstand bacterial CIP stress would also affect eukaryotic growth. The present study thus indicates that a detailed assessment of the aquatic ecotoxicity of CIP is essential, as the effects of CIP are much more complicated in microbial communities than in monocultures. CIP will continue to be an environmental contaminant due to its wide usage and production and more attention should be given to the negative effects of antibiotics as well as other bioactive pollutants on aquatic environments.

ABC Transporters in Prorocentrum lima and Their Expression Under Different Environmental Conditions Including Okadaic Acid Production

Prorocentrum lima is a typical benthic toxic dinoflagellate, which can produce phycotoxins such as okadaic acid (OA). In this study, we identified three ABC transporter genes (ABCB1, ABCC1 and ABCG2) and characterized their expression patterns, as well as OA production under different environmental conditions in P. lima. We found that the three ABC transporters all showed high identity with related ABC proteins from other species, and contained classical features of ABC transport proteins. Among them, ABCG2 was a half size transporter. The three ABC transporter genes displayed various expression profiles under different conditions. The high concentration of Cu²⁺ could up-regulate ABCB1, ABCC1 and ABCG2 transcripts in P. lima, suggesting the potential defensive role of ABC transporters against metal ions in surrounding waters. Cu²⁺, in some concentration, could induce OA production; meanwhile, tributyltin inhibited OA accumulation. The grazer Artemia salina could induce OA production, and P. lima displayed some toxicity to the grazer, indicating the possibility of OA as an anti-grazing chemical. Collectively, our results revealed intriguing data about OA production and the expression patterns of three ABC transporter genes. However, we could not find any significant correlation between OA production and expression pattern of the three ABC transporters in P. lima. Our results might provide new molecular insights on the defensive responses of P. lima to the surrounding environment.

Transcriptomic Responses to Thermal Stress and Varied Phosphorus Conditions in Fugacium kawagutii

Coral reef-associated Symbiodiniaceae live in tropical and oligotrophic environments and are prone to heat and nutrient stress. How their metabolic pathways respond to pulses of warming and phosphorus (P) depletion is underexplored. Here, we conducted RNA-seq analysis to investigate transcriptomic responses to thermal stress, phosphate deprivation, and organic phosphorus (OP) replacement in Fugacium kawagutii. Using dual-algorithm (edgeR and NOIseq) to remedy the problem of no replicates, we conservatively found 357 differentially expressed genes (DEGs) under heat stress, potentially regulating cell wall modulation and the transport of iron, oxygen, and major nutrients. About 396 DEGs were detected under P deprivation and 671 under OP utilization, both mostly up-regulated and potentially involved in photosystem and defensome, despite different KEGG pathway enrichments. Additionally, we identified 221 genes that showed relatively stable expression levels across all conditions (likely core genes), mostly catalytic and binding proteins. This study reveals a wide range of, and in many cases previously unrecognized, molecular mechanisms in F. kawagutii to cope with heat stress and phosphorus-deficiency stress. Their quantitative expression dynamics, however, requires further verification with triplicated experiments, and the data reported here only provide clues for generating testable hypotheses about molecular mechanisms underpinning responses and adaptation in F. kawagutii to temperature and nutrient stresses.

Unraveling the molecular mechanism of the response to changing ambient phosphorus in the dinoflagellate Alexandrium catenella with quantitative proteomics

Phosphorus (P) is a key macronutrient limiting cell growth and bloom formation of marine dinoflagellates. Physiological responses to changing ambient P have been investigated in dinoflagellates; however, the molecular mechanisms behind these responses remain limited. Here, we compared the protein expression profiles of a marine dinoflagellate Alexandrium catenella grown in inorganic P-replete, P-deficient, and inorganic- and organic-P resupplied conditions using an iTRAQ-based quantitative proteomic approach. P deficiency inhibited cell growth and enhanced alkaline phosphatase activity (APA) but had no effect on photosynthetic efficiency. After P resupply, the P-deficient cells recovered growth rapidly and APA decreased. Proteins involved in sphingolipid metabolism, organic P utilization, starch and sucrose metabolism, and photosynthesis were up-regulated in the P-deficient cells, while proteins associated with protein synthesis, nutrient assimilation and energy metabolism were down-regulated. The responses of the P-deficient A. catenella to the resupply of organic and inorganic P presented significant differences: more biological processes were enhanced in the organic P-resupplied cells than those in the inorganic P-resupplied cells; A. catenella might directly utilize G-6-P for nucleic acid synthesis through the pentose phosphate pathway. Our results indicate that A. catenella has evolved diverse adaptive strategies to ambient P deficiency and specific mechanisms to utilize dissolved organic P, which might be an important reason resulting in A. catenella bloom in the low inorganic P environment. BIOLOGICAL SIGNIFICANCE: The ability of marine dinoflagellates to utilize different phosphorus (P) species and adapt to ambient P deficiency determines their success in the ocean. In this study, we investigated the response mechanisms of a dinoflagellate Alexandrium catenella to ambient P deficiency, and resupply of inorganic- and organic-P at the proteome level. Our results indicated that A. catenella initiated multiple adaptive strategies to ambient P deficiency, e.g. utilizing nonphospholipids and glycosphingolipids instead of phospholipids, enhancing expression of acid phosphatase to utilize organic P, and reallocating intracellular energy. Proteome responses of the P-deficient A. catenella to resupply of inorganic- and organic-P differed significantly, indicating different utilization pathways of inorganic and organic P, A. catenella might directly utilize low molecular weight organic P, such as G-6-P as both P and carbon sources.

Metatranscriptomic Signatures Associated With Phytoplankton Regime Shift From Diatom Dominance to a Dinoflagellate Bloom

Diatoms and dinoflagellates dominate coastal marine phytoplankton communities as major players of marine biogeochemical cycles and their seasonal succession often leads to harmful algal blooms (HABs). What regulates their respective dominances and the development of the HABs remains elusive. Here we conducted time-sequential metatranscriptomic profiling on a natural assemblage that evolved from diatom dominance to a dinoflagellate bloom to interrogate the underlying major metabolic and ecological drivers. Data reveals similarity between diatoms and dinoflagellates in exhibiting high capacities of energy production, nutrient acquisition, and stress protection in their respective dominance stages. The diatom-to-dinoflagellate succession coincided with an increase in turbidity and sharp declines in silicate and phosphate availability, concomitant with the transcriptomic shift from expression of silicate uptake and urea utilization genes in diatoms to that of genes for light harvesting, diversified phosphorus acquisition and autophagy-based internal nutrient recycling in dinoflagellates. Furthermore, the diatom-dominant community featured strong potential to carbohydrate metabolism and a strikingly high expression of trypsin potentially promoting frustule building. In contrast, the dinoflagellate bloom featured elevated expression of xanthorhodopsin, and antimicrobial defensin genes, indicating potential importance of energy harnessing and microbial defense in bloom development. This study sheds light on mechanisms potentially governing diatom- and dinoflagellate-dominance and regulating bloom development in the natural environment and raises new questions to be addressed in future studies.

Transcriptional Shifts Highlight the Role of Nutrients in Harmful Brown Tide Dynamics

Harmful algal blooms (HABs) threaten ecosystems and human health worldwide. Controlling nitrogen inputs to coastal waters is a common HAB management strategy, as nutrient concentrations often suggest coastal blooms are nitrogen-limited. However, defining best nutrient management practices is a long-standing challenge: in part, because of difficulties in directly tracking the nutritional physiology of harmful species in mixed communities. Using metatranscriptome sequencing and incubation experiments, we addressed this challenge by assaying the in situ physiological ecology of the ecosystem destructive alga, Aureococcus anophagefferens. Here we show that gene markers of phosphorus deficiency were expressed in situ, and modulated by the enrichment of phosphorus, which was consistent with the observed growth rate responses. These data demonstrate the importance of phosphorus in controlling brown-tide dynamics, suggesting that phosphorus, in addition to nitrogen, should be evaluated in the management and mitigation of these blooms. Given that nutrient concentrations alone were suggestive of a nitrogen-limited ecosystem, this study underscores the value of directly assaying harmful algae in situ for the development of management strategies.

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

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

Genomic analysis of red-tide water bloomed with Heterosigma akashiwo in Geoje

Microbial community structures of harmful algal bloom (HAB) caused by Heterosigma akashiwo in Geoje were analyzed using the MiSeq platform. To analyze phytoplankton communities without cross-reactivity with predominant bacteria, a new phytoplankton-specific 23S universal primer set was designed by modifying two previously used ones. The new universal primer set turned out to be a useful tool for the analysis of the phytoplankton community; it showed a high specificity for phytoplankton without cross-reactivity to bacterial sequences as well as the wide taxon coverage presenting from prokaryotic cyanobacteria to eukaryotic algae. Next Generation Sequencing (NGS) data generated by two universal primer sets (16S and 23S) provided useful information about the H. akashiwo bloom. According to the 23S universal primer set, proportions of H. akashiwo increased by more than 200-fold as the bloom occurred and its numbers were high enough to detect in control sites. Its operational taxonomic units (OTUs) were detected in the bloom sites at low proportions suggesting that the 16S universal primer set may not be as effective for monitoring harmful algal blooming (HAB) as the 23S universal primer set. In addition, several abundant OTUs in Chlorophyta were not presented by the 16S universal primer set in this study. However, the 16S primer set was useful for detecting decreases in Foraminifera as HAB occurred suggesting that genomic analyses using two universal primer sets would provide more reliable data for understanding microbial community changes by various environmental or ecological events, including HAB. Genomic analyses using two universal primer sets was also useful for determining a correlation between microbial components as HAB occurred. Heterosigma akashiwo was positively correlated with other bloom species, including Karenia mikimotoi, Teleaulax amphioxeia, and bacteria in Verrucomicrobia.

Inorganic mercury (Hg2+) accumulation in autotrophic and mixotrophic planktonic protists: Implications for Hg trophodynamics in ultraoligotrophic Andean Patagonian lakes

Microbial assemblages are typical of deep ultraoligotrophic Andean Patagonian lakes and comprise picoplankton and protists (phytoflagellates and mixotrophic ciliates), having a central role in the C cycle, primary production and in the incorporation of dissolved inorganic mercury (Hg²⁺) into lake food webs. In this study we evaluated the mechanisms of Hg²⁺ incorporation in hetero- and autotrophic bacteria, in the autotrophic dinoflagellate (Gymnodinium paradoxum) and in two mixotrophic ciliates (Stentor araucanus and Ophrydium naumanni) dominating the planktonic microbial assemblage. The radioisotope ¹⁹⁷Hg was used to trace the Hg²⁺ incorporation in microbiota. Hg uptake was analyzed as a function of cell abundance (BCF: bioconcentration factor), cell surface (SCF: surface concentration factor) and cell volume (VCF: volume concentration factor). Overall, the results obtained showed that these organisms incorporate substantial amounts of dissolved Hg²⁺ passively (adsorption) and actively (bacteria consumption or attachment), displaying different Hg internalization and therefore, varying potential for Hg transfer. Surface area and quality, and surface:volume ratio (S:V) control the passive uptake in all the organisms. Active incorporation depends on bacteria consumption in the mixotrophic ciliates, or on bacteria association to surface in the autotrophic dinoflagellate. Hg bioaccumulated by pelagic protists can be transferred to higher trophic levels through plankton and fish feeding, regenerated to the dissolved phase by excretion, and/or transferred to the sediments by particle sinking. In ultraoligotrophic Andean Patagonian lakes, picoplankton and planktonic protists are key components of lake food webs, linking the pelagic and benthic Hg pathways, and thereby playing a central role in Hg trophodynamics.