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

Transcriptome response of a marine copepod in response to environmentally-relevant concentrations of saxitoxin

Paralytic shellfish toxins (PSTs) can pose a serious threat to human health. Among them, saxitoxin (STX) is one of the most potent natural neurotoxins. Here, the copepod Tigriopus japonicus, was exposed to environmentally relevant concentrations (2.5 and 25 μg/L) STX for 48 h. Although no lethal effects were observed at both concentrations, the transcriptome was significantly altered, and displayed a concentration-dependent response. STX exposure decreased the copepod's metabolism and compromised immune defense and detoxification. Additionally, STX disturbed signal transduction, which might affect other cellular processes. STX exposure could inhibit the copepod's chitin metabolism, disrupting its molting process. Also, the processes related to damage repair and protection were up-regulated to fight against high concentration exposure. Collectively, this study has provided an early warning of PSTs for coastal ecosystem not only because of their potent toxicity effect but also their bioaccumulation that can transfer up the food chain after ingestion by copepods.

Transoceanic ships as a source of alien dinoflagellate invasions of inland freshwater ecosystems

Ships' ballast water and sediments have long been linked to the global transport and expansion of invasive species and thus have become a hot research topic and administrative challenge in the past decades. The relevant concerns, however, have been mainly about the ocean-to-ocean invasion and sampling practices have been almost exclusively conducted onboard. We examined and compared the dinoflagellate cysts assemblages in 49 sediment samples collected from ballast tanks of international and domestic routes ships, washing basins associated with a ship-repair yard, Jiangyin Port (PS), and the nearby area of Yangtze River (YR) during 2017-2018. A total of 43 dinoflagellates were fully identified to species level by metabarcoding, single-cyst PCR-based sequencing, cyst germination and phylogenetic analyses, including 12 species never reported from waters of China, 14 HABs-causing, 9 toxic, and 10 not strictly marine species. Our metabarcoding and single-cyst sequencing also detected many OTUs and cysts of dinoflagellates that could not be fully identified, indicating ballast tank sediments being a risky repository of currently unrecognizable invasive species. Particularly important, 10 brackish and fresh water species of dinoflagellate cysts (such as Tyrannodinium edax) were detected from the transoceanic ships, indicating these species may function as alien species potentially invading the inland rivers and adjacent lakes if these ships conduct deballast and other practices in fresh waterbodies. Significantly higher numbers of reads and OTUs of dinoflagellates in the ballast tanks and washing basins than that in PS and YR indicate a risk of releasing cysts by ships and the associated ship-repair yards to the surrounding waters. Phylogenetic analyses revealed high intra-species genetic diversity for multiple cyst species from different ballast tanks. Our work provides novel insights into the risk of bio-invasion to fresh waters conveyed in ship's ballast tank sediments and washing basins of shipyards.

Uncovering nutrient regeneration, transformation pattern, and its contribution to harmful algal blooms in mariculture waters

The prevalence of harmful algal blooms (HABs), especially in mariculture waters, has become a concern for environmental and human health worldwide. Notably, the frequent occurrence of HABs relies upon a substantial supply of available nutrients, which are influenced by nutrient recycling. However, nutrient regeneration, transformation pattern, and their contribution to HABs in mariculture waters remain largely unknown. In this study, by combining field investigation and incubation experiments from June to September 2020, the temporal variations in nutrients and algal composition were revealed. In addition, the nutrient regeneration and assimilation rates in the water column during two continuous algal blooms were measured. The results indicated that organic nutrients, which were the dominant components, strongly stimulated nutrient regeneration. High regeneration rates were observed, with dissolved inorganic nitrogen (DIN) and phosphorous (DIP) regeneration rates ranging from 0.25 to 2.64 μmol/L·h and 0.01 to 0.09 μmol/L·h, respectively. Compared to the direct uptake of organic nutrients, the rapid regeneration of inorganic nutrients played a vital role in sustaining continuous algal blooms, as regenerated DIN contributed 100 % while regenerated DIP contributed 72-100 % of the algal assimilation demand. Furthermore, the redundancy analysis and inverse solution equations indicated that different N transformation patterns and utilization strategies occurred during Heterosigma and Nannochloris blooms. The shorter N recycling pathway and faster NH4+ supply rates provided favorable conditions for the dominance of Nannochloris over Heterosigma, which had a preference for the uptake of NO3-. In conclusion, we propose that nutrient regeneration is a key maintenance mechanism underlying the maintenance of continuous algal blooms, and different N transformation patterns and utilization strategies regulate algal communities in mariculture waters.

Diatom-dinoflagellate succession in the Bohai Sea: The role of N/P ratios and dissolved organic nitrogen components

Dissolved organic nitrogen (DON) is a pivotal component of total dissolved nitrogen pools, serving as a crucial nitrogen source for phytoplankton. This study investigated the impact of nitrogen-to-phosphorus (N/P) ratios and different DON components (hydrophilic vs hydrophobic DON) on diatom-dinoflagellate succession through field culture experiments. Results showed that dinoflagellates have a competitive advantage under high N/P ratios and phosphorus limitation, regardless of DON or DIN treatments. Hydrophilic DON exhibits greater bioavailability than hydrophobic DON (40.6% vs. 21.7 %), resulting in increased algal biomass and diatoms dominance in the community. Additionally, DON was categorized into labile and refractory components (LDON and RDON) based on bioavailability. LDON primarily consists of protein-like components that can be readily consumed by algae, whereas RDON is primarily composed of humic-like components that are less accessible to algae. Diatoms and dinoflagellates exhibited differential responses to LDON and RDON, with diatoms thriving in high LDON environments, while dinoflagellates gained a competitive advantage when RDON was the predominant nitrogen source. Furthermore, a significant negative correlation was observed between bioavailable nitrogen concentration (BAN: DIN + LDON) and the ratio of dinoflagellates to diatoms (p<0.05). In conclusion, our study highlights the role of LDON in promoting diatom dominance, whereas environments dominated by RDON foster dinoflagellate success. These findings enhance our comprehension of diatom-dinoflagellate succession dynamics.

Transcriptome analysis of the harmful alga Heterosigma akashiwo under a 24-hour light-dark cycle

The photoperiod, which is defined as the period of time within a 24-hour time frame that light is available, is an important environmental regulator of several physiological processes in phytoplankton, including harmful bloom-forming phytoplankton. The ichthyotoxic raphidophyte Heterosigma akashiwo is a globally distributed bloom-forming phytoplankton. Despite extensive studies on the ecological impact of H. akashiwo, the influence of the photoperiod on crucial biological processes of this species remains unclear. In this study, gene expression in H. akashiwo was analyzed over a 24-hour light-dark (14:10) treatment period. Approximately 36 % of unigenes in H. akashiwo were differentially expressed during this 24-hour treatment period, which is indicative of their involvement in the response to light-dark variation. Notably, the number of differentially expressed genes exhibited an initial increase followed by a subsequent decrease as the sampling time progressed (T0 vs. other time points). Unigenes associated with photosynthesis and photoprotection reached their peak expression levels after 2-4 h of illumination (T12-T14). In contrast, the expression of unigenes associated with DNA replication peaked at the starting point of the dark period (T0). Furthermore, although several unigenes annotated to photoreceptors displayed potential diel periodicity, genes from various photoreceptor families (such as phytochrome and cryptochrome) showed unique expression patterns. Collectively, our findings offer novel perspectives on the response of H. akashiwo to the light-dark cycle, serving as a valuable resource for investigating the physiology and ecology of this species.

Harmful algal blooms in China: History, recent expansion, current status, and future prospects

The impacts of harmful algal blooms (HABs) on economies, public health, ecosystems, and aquaculture across the globe have all increased in recent decades, and this has been acutely the case in China. Here, we review the history of HABs and HABs research in China, as well as recent trends in HABs and future prospects of HAB science in China. The most updated analyses demonstrated that the number of HAB events, the number of HAB species, the aerial coverage of HABs, and the impacts of HABs in Chinese waters during the 21st century were all higher than that during the last two decades of the 20th century. The increase in the number of HABs in China has been significantly correlated with the increased discharge of ammonium and total phosphorus into coastal waters (p < 0.01 for both). Notable newly recognized events this century have included chronic HABs caused by Prorocentrum donghaiense and Karenia mikimotoi, a paralytic shellfish poisoning event caused by Gymnodinium catenatum that sickened 80 people, brown tides caused by Aureococcus anophagefferens, green tides caused by Ulva prolifera, golden tides caused by Sargassum horneri, and the disruption of a nuclear power plant caused by a bloom of Phaeocystis globosa. A series of key discoveries regarding HABs has been made this century including documentation of nearly all known HAB toxins in Chinese waters, discovery of novel cyst-formation and/or life stages of multiple HABs-causing species, identification of the chemical and physical oceanographic drivers of multiple HABs including those formed by P. donghaiense, K. mikimotoi, and U. prolifera, and the successful mitigation of HABs via the use of modified clay approaches. Future research prospects highlighted include the use of macroalgae as a means to prevent, mitigate, and control (PCM) HABs and the process by which multi-disciplinary studies involving molecular approaches (omics), remote in situ detection, artificial intelligence, and mega-data analyses might be used to develop refined and realistic HAB forecasting platforms. Collectively, this review demonstrates the significant evolution of HAB science since the 20th century in China and demonstrates that while HABs in China are complex and widespread, recent and on-going discoveries make the development of detailed understanding and effective measures to mitigate the negative effects of HABs a hopeful outcome in the coming years.

The molecular mechanisms and environmental effects of modified clay control algal blooms in aquacultural water

Modified clay (MC) technology is an effective method for controlling harmful algal blooms (HABs). Based on field experience, a bloom does not continue after treatment with MC, even though the residual HAB biomass accounts for 20-30% of the initial biomass. Laboratory studies using unialgal cultures have found that MC could inhibit the growth of the residual algal cells to prevent HABs. Nevertheless, the phytoplankton in field waters is diverse. Therefore, unclassified complex mechanisms may exist. To illustrate the molecular mechanisms through which MC controls HABs in the field and verify the previous laboratory findings, a series of experiments and bioinformatics analyses were conducted using bloom waters from aquacultural ponds. The results showed that a 72.29% removal efficiency of algal biomass could effectively control blooms. The metatranscriptomic results revealed that the number of downregulated genes (131,546) was greater than that of upregulated genes (24,318) at 3 h after MC addition. Among these genes, several genes related to DNA replication were downregulated; however, genes involved in DNA repair were upregulated. Metabolism-related pathways were the most significantly upregulated (q < 0.05), including photosynthesis and oxidative phosphorylation. The results also showed that MC reduced most of the biomass of the dominant phytoplankton species, likely by removing apical dominance, which increased the diversity and stability of the phytoplankton community. In addition to reducing the pathogenic bacterial density, MC reduced the concentrations of PO₄^3- (96.22%) and SiO₃^2- (66.77%), thus improving the aquaculture water quality, altering the phytoplankton community structure (the proportion of Diatomea decreased, and that of Chlorophyta increased), and inhibiting phytoplankton growth. These effects hindered the rapid development of large phytoplankton biomasses and allowed the community structure to remain stable, reducing HAB threats. This study illustrates the molecular mechanisms through which MC controls HABs in the field and provides a scientific method for removing HABs in aquacultural waters.

Multi-trait analysis reveals large interspecific differences for phytoplankton in response to thermal change

Understanding the responses of multiple traits in phytoplankton, and identifying interspecific variabilities to thermal changes is crucial for predicting the impacts of ocean warming on phytoplankton distributions and community structures in future scenarios. Here, we applied a trait-based approach by examining the patterns in multi-traits variations (eight traits) and interspecific variabilities in five phytoplankton species (two diatoms, three dinoflagellates) in response to a wide range of ecologically relevant temperatures (14-30 °C). Our results show large inter-traits and interspecific variabilities of thermal reaction norms in all of the tested traits. We also found that the interspecific variability exceeded the variations induced by thermal changes. Constrained variations and trade-offs between traits both revealed substantial interspecific differences and shifted as the temperature changed. Our study helps to understand the species-specific response patterns of multiple traits to ocean warming and to investigate the implications of these responses in the context of global change.

Occurrence and distribution of lipophilic marine algal toxins in the coastal seawater of Southeast China and the South China Sea

The composition, levels, and spatial distribution of dissolved lipophilic marine algal toxins (LMATs) including cyclic imines (CIs), yessotoxins (YTXs), okadaic acid (OA) and its derivatives, pectenotoxins (PTXs), azaspiracids (AZAs), and brevetoxins (BTXs) in the coastal waters of Southeast China (Xiamen) and the South China Sea (Hainan Island and Beibu Gulf) were investigated and compared for the first time. Dissolved AZA3 was firstly detected in the coastal seawater of China. OA and PTX2 were widely distributed in the three areas studied. Gymnodimine (GYM), 13-desmethyl spirolide C (SPX1), YTX, and homo-yessotoxins (h-YTX) were found mainly in the South China Sea. The average ∑LMAT concentrations in the coastal waters of Xiamen, Hainan Island, and Beibu Gulf were 10.02 ng/L, 4.21 ng/L, and 44.27 ng/L, respectively. More groups and much higher concentrations of LMATs occurred in the South China Sea than that in the other sea areas of China.

Effect of ferrous-activated calcium peroxide oxidation on forward osmosis treatment of algae-laden water: Membrane fouling mitigation and mechanism

Forward osmosis (FO) is a high-efficiency and low-energy consumption way for algae-laden water treatment, whereas membrane fouling is still an unavoidable problem in its practical application. In this work, a strategy of ferrous-activated calcium peroxide (Fe(II)/CaO₂) was proposed to control FO membrane fouling in the purification of algae-laden water. With the treatment of Fe(II)/CaO₂, the aggregation of algal contaminants was promoted, the cell viability and integrity were well preserved, and the fluorescent organics were efficiently removed. With respect to the fouling of FO membrane, the flux decline was generally alleviated, and the flux recovery was promoted to varying degrees under different process conditions. It could be revealed through the extended Derjaguin-Landau-Verwey-Overbeek theory that the adhesion of contaminants and membrane surfaces was reduced by Fe(II)/CaO₂ treatment. The interface morphologies and functional groups of membrane verified that Fe(II)/CaO₂ could mitigate the fouling by reducing the amount of algal contaminants adhering to the FO membrane. The co-coagulation of in-situ Fe(III) together with Ca(OH)₂, as well as the oxidation of ^•OH were the main mechanisms for fouling mitigation. In sum, the Fe(II)/CaO₂ process could effectively improve the efficiency of FO for algae-laden water treatment, and has broad application prospects.

Impacts of Climate Change on the Biogeography of Three Amnesic Shellfish Toxin Producing Diatom Species

Harmful algal blooms (HABs) are considered one of the main risks for marine ecosystems and human health worldwide. Climate change is projected to induce significant changes in species geographic distribution, and, in this sense, it is paramount to accurately predict how it will affect toxin-producing microalgae. In this context, the present study was intended to project the potential biogeographical changes in habitat suitability and occurrence distribution of three key amnesic shellfish toxin (AST)-producing diatom species (i.e., Pseudo-nitzschia australis, P. seriata, and P. fraudulenta) under four different climate change scenarios (i.e., RCP-2.6, 4.5, 6.0, and 8.5) up to 2050 and 2100. For this purpose, we applied species distribution models (SDMs) using four abiotic predictors (i.e., sea surface temperature, salinity, current velocity, and bathymetry) in a MaxEnt framework. Overall, considerable contraction and potential extirpation were projected for all species at lower latitudes together with projected poleward expansions into higher latitudes, mainly in the northern hemisphere. The present study aims to contribute to the knowledge on the impacts of climate change on the biogeography of toxin-producing microalgae species while at the same time advising the correct environmental management of coastal habitats and ecosystems.

Projecting Future Climate Change-Mediated Impacts in Three Paralytic Shellfish Toxins-Producing Dinoflagellate Species

Simple Summary

Harmful algal blooms present a particular risk for marine ecosystems and human health alike. In this sense, it is important to accurately predict how toxin-producing microalgae could be affected by future climate change. The present study applied species distribution models (SDMs) to project the potential changes in the habitat suitability and distribution of three key paralytic shellfish toxin (PST)-producing dinoflagellate species (i.e., Alexandrium catenella, A. minutum, and Gymnodinium catenatum), up to 2040/50 and 2090/2100, across four different greenhouse gas emission scenarios, and using four abiotic predictors (i.e., sea surface temperature, salinity, current velocity, and bathymetry). In general, considerable contractions were observed for all three species in the lower latitudes of their distribution, together with projected expansions into higher latitudes, particularly in the Northern Hemisphere. This study aims to entice further research on the future biogeographical impacts of climate change in toxin-producing microalgae species while, at the same time, helping to advise the correct environmental management of coastal habitats and ecosystems.

Abstract

Toxin-producing microalgae present a significant environmental risk for ecosystems and human societies when they reach concentrations that affect other aquatic organisms or human health. Harmful algal blooms (HAB) have been linked to mass wildlife die-offs and human food poisoning episodes, and climate change has the potential to alter the frequency, magnitude, and geographical extent of such events. Thus, a framework of species distribution models (SDMs), employing MaxEnt modeling, was used to project changes in habitat suitability and distribution of three key paralytic shellfish toxin (PST)-producing dinoflagellate species (i.e., Alexandrium catenella, A. minutum, and Gymnodinium catenatum), up to 2050 and 2100, across four representative concentration pathway scenarios (RCP-2.6, 4.5, 6.0, and 8.5; CMIP5). Despite slightly different responses at the regional level, the global habitat suitability has decreased for all the species, leading to an overall contraction in their tropical and sub-tropical ranges, while considerable expansions are projected in higher latitudes, particularly in the Northern Hemisphere, suggesting poleward distributional shifts. Such trends were exacerbated with increasing RCP severity. Yet, further research is required, with a greater assemblage of environmental predictors and improved occurrence datasets, to gain a more holistic understanding of the potential impacts of climate change on PST-producing species.

Spatiotemporal changes of bacterial communities during a cyanobacterial bloom in a subtropical water source reservoir ecosystem in China

Cyanobacterial blooms, which not only threaten the health and stability of aquatic ecosystems but also influence the microbial community within, emerges as one of the most concerning problems in China. However, how cyanobacterial blooms affect the spatiotemporal variation of aquatic microbial communities remains relatively unclear. In this study, we used high-throughput sequencing to investigate how the cyanobacterial and bacterial community spatiotemporally vary along with main cyanobacterial bloom phases in upstream rivers of a eutrophicated water source reservoir. Both cyanobacterial and bacterial diversities in each river were significantly lower (P < 0.05) during the bloom outbreak phase, showing the apparent influence of cyanobacterial bloom. Dominant cyanobacterial taxa included Cyanobacteriales and Synechococcales, and dominant bacterial taxa comprised Acinetobacter, CL500-29, hgcI clade, Limnohabitans, Flavobacterium, Rhodoluna, Porphyrobacter, Rhodobacter, Pseudomonas, and Rhizobiales, whose changes of relative abundance along with the bloom indicated distinct community composition. Non-metric multidimensional scaling analysis proved that community composition had significant difference amongst bloom phases. Linear discriminant analysis (LDA) with LDA effect size analysis (LEfSe) identified unique dominant cyanobacterial and bacterial OTUs at different phases in each river, indicating spatiotemporal variations of communities. Canonical correlation analysis or redundancy analysis revealed that at different bloom phases communities of each river had distinct correlation patterns with the environmental parameters (temperature, ammonium, nitrate, and total phosphorus etc.), implying the spatial variations of microbial communities. Overall, these results expand current understanding on the spatiotemporal variations of microbial communities due to cyanobacterial blooms. Microbial interactions during the bloom may shed light on controlling cyanobacterial blooms in the similar aquatic ecosystems.

Expression Plasticity of Peroxisomal Acyl-Coenzyme A Oxidase Genes Implies Their Involvement in Redox Regulation in Scallops Exposed to PST-Producing Alexandrium

Filter-feeding bivalves can accumulate paralytic shellfish toxins (PST) produced by toxic microalgae, which may induce oxidative stress and lipid peroxidation. Peroxisomal acyl-coenzyme A oxidases (ACOXs) are key enzymes functioning in maintaining redox and lipid homeostasis, but their roles in PST response in bivalves are less understood. Herein, a total of six and six ACOXs were identified in the Chlamys farreri and Patinopecten yessoensis genome, respectively, and the expansion of ACOX1s was observed. Gene expression analysis revealed an organ/tissue-specific expression pattern in both scallops, with all ACOXs being predominantly expressed in the two most toxic organs, digestive glands and kidneys. The regulation patterns of scallop ACOXs after exposure to different PST-producing algaes Alexandrium catenella (ACDH) and A. minutum (AM-1) were revealed. After ACDH exposure, more differentially expressed genes (DEGs) were identified in C. farreri digestive glands (three) and kidneys (five) than that in P. yessoensis (two), but the up-regulated DEGs showed similar expression patterns in both species. In C. farreri, three DEGs were found in both digestive glands and kidneys after AM-1 exposure, with two same CfACOX1s being acutely and chronically induced, respectively. Notably, these two CfACOX1s also showed different expression patterns in kidneys between ACDH (acute response) and AM-1 (chronic response) exposure. Moreover, inductive expression of CfACOXs after AM-1 exposure was observed in gills and mantles, and all DEGs in both tissues were up-regulated and their common DEGs exhibited both acute and chronic induction. These results indicate the involvement of scallop ACOXs in PST response, and their plasticity expression patterns between scallop species, among tissues, and between the exposure of different PST analogs.

Co-occurrence of chromophytic phytoplankton and the Vibrio community during Phaeocystis globosa blooms in the Beibu Gulf

Accumulating research evidence has revealed that harmful algal blooms (HABs) can substantially affect the community structures of phytoplankton and heterotrophic bacteria in marine ecosystems. However, little is known about their species-specific interactions between phytoplankton and heterotrophic bacteria during the HABs period and about their interaction shifts in response to blooms. From this perspective, we investigated the co-occurrence of chromophytic phytoplankton and Vibrio during Phaeocystis globosa blooms in the Beibu Gulf. The results showed that Vibrio communities were distinct during the blooms, and P. globosa blooms resulted in a decline in phytoplankton alpha diversity, revealing that the blooms could affect their community compositions. The regression lines between the Shannon indices and Bray-Curtis distances of phytoplankton and Vibrio showed positive correlations with each other (p < 0.001), suggesting that they may have intrageneric symbiotic interactions overall. In addition, network analysis further demonstrated that relationships between phytoplankton and Vibrio were dominated by positive correlations, and more interaction modules were observed during the blooms, revealing that the blooms intensified synergistic association and mutual symbiotic interactions between them. Environmental factors (SiO₃^2-, NH₄⁺, NO₃^- and TN,) and P. globosa density more deeply affected network interactions between phytoplankton and Vibrio during the periods of P. globosa blooms than those before the blooms and after the blooms. This study provided new insight to elucidate community structure and interaction relationships between phytoplankton and Vibrio in response to P. globosa blooms and their ecological effects in marine ecosystems.