Metatranscriptome analysis reveals environmental and diel regulation of a Heterosigma akashiwo (raphidophyceae) bloom

The allelopathic effects of Heterosigma akashiwo on Skeletonema costatum: Insights from gene expression and metabolomics analysis

The globally distributed harmful algal blooms (HAB) species, Heterosigma akashiwo, has been found to exhibit ichthyotoxicity. Previous studies have shown that H. akashiwo achieves a competitive edge during bloom occurrences by inhibiting the growth of a coexisting diatom, Skeletonema costatum, through allelopathy. However, the specific allelopathic mechanisms underlying the allelopathic effects of H. akashiwo on S. costatum remain unknown. To bridge this gap, our study utilized a combination of quantitative real-time PCR and metabolomics to examine the allelopathic processes of H. akashiwo on S. costatum. Our results demonstrate that the growth of S. costatum is hindered when co-cultured with H. akashiwo (initial cell concentration, 2 × 104 cell/mL). Gene expression investigation showed a substantial reduction in the mRNA levels of cytochrome b6, ribulose bisphosphate carboxylase large chain, and silicon transporter in S. costatum when grown in co-culture conditions. Furthermore, metabolic pathway analysis suggested that the allelopathic effects of H. akashiwo disrupted several vital metabolic pathways in S. costatum, including a reduction in purine and pyrimidine metabolism and an increase in fatty acid biosynthesis. Our investigation has revealed the intricate and substantial involvement of allelopathy in the formation of H. akashiwo blooms, demonstrating the complexity of the allelopathic interaction between H. akashiwo and S. costatum. These insights also contribute significantly to our understanding of the dynamics within HAB species.

Investigating the molecular mechanisms of Pseudalteromonas sp. LD-B1's algicidal effects on the harmful alga Heterosigma akashiwo

Heterosigma akashiwo is a harmful algal bloom species that causes significant detrimental effects on marine ecosystems worldwide. The algicidal bacterium Pseudalteromonas sp. LD-B1 has demonstrated potential effectiveness in mitigating these blooms. However, the molecular mechanisms underlying LD-B1's inhibitory effects on H. akashiwo remain poorly understood. In this study, we employed the comprehensive methodology, including morphological observation, assessment of photosynthetic efficiency (Fv/Fm), and transcriptomic analysis, to investigate the response of H. akashiwo to LD-B1. Exposure to LD-B1 resulted in a rapid decline of H. akashiwo's Fv/Fm ratio, with cells transitioning to a rounded shape within 2 hours, subsequently undergoing structural collapse and cytoplasmic leakage. Transcriptomic data revealed sustained downregulation of photosynthetic genes, indicating impaired functionality of the photosynthetic system. Additionally, genes related to the respiratory electron transfer chain and antioxidant defenses were consistently downregulated, suggesting prolonged oxidative stress beyond the cellular antioxidative capacity. Notably, upregulation of autophagy-related genes was observed, indicating autophagic responses in the algal cells. This study elucidates the molecular basis of LD-B1's algicidal effects on H. akashiwo, advancing our understanding of algicidal mechanisms and contributing to the development of effective strategies for controlling harmful algal blooms.

In situ community transcriptomics illuminates CO2-fixation potentials and supporting roles of phagotrophy and proton pump in plankton in a subtropical marginal sea

Lineage-wise physiological activities of plankton communities in the ocean are important but challenging to characterize. Here, we conducted whole-assemblage metatranscriptomic profiling at continental shelf and slope sites in the South China Sea to investigate carbon fixation potential in different lineages. RuBisCO expression, the proxy of Calvin carbon fixation (CCF) potential, was mainly contributed by Bacillariophyta, Chlorophyta, Cyanobacteria, and Haptophyta, which was differentially affected by environmental factors among lineages. CCF potential exhibited positive or negative correlations with phagotrophy gene expression, suggesting phagotrophy possibly enhances or complements CCF. Our data also reveal significant non-Calvin carbon fixation (NCF) potential, as indicated by the active expression of genes in all five currently recognized NCF pathways, mainly contributed by Flavobacteriales, Alteromonadales, and Oceanospirillales. Furthermore, in Flavobacteriales, Alteromonadales, Pelagibacterales, and Rhodobacterales, NCF potential was positively correlated with proton-pump rhodopsin (PPR) expression, suggesting that NCF might be energetically supported by PPR. The novel insights into the lineage-differential potential of carbon fixation, widespread mixotrophy, and PPR as an energy source for NCF lay a methodological and informational foundation for further research to understand carbon fixation and the trophic landscape in the ocean.IMPORTANCEMarine plankton plays an important role in global carbon cycling and climate regulation. Phytoplankton and cyanobacteria fix CO2 to produce organic compounds using solar energy and mainly by the Calvin cycle, whereas autotrophic bacteria and archaea may fix CO2 by non-Calvin cycle carbon fixation pathways. How active individual lineages are in carbon fixation and mixotrophy, and what energy source bacteria may employ in non-Calvin carbon fixation, in a natural plankton assemblage are poorly understood and underexplored. Using metatranscriptomics, we studied carbon fixation in marine plankton with lineage resolution in tropical marginal shelf and slope areas. Based on the sequencing results, we characterized the carbon fixation potential of different lineages and assessed Calvin- and non-Calvin- carbon fixation activities and energy sources. Data revealed a high number of unigenes (4.4 million), lineage-dependent differential potentials of Calvin carbon fixation and responses to environmental conditions, major contributors of non-Calvin carbon fixation, and their potential energy source.

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.

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.

Mechanisms of Phaeocystis globosa blooms in the Beibu Gulf revealed by metatranscriptome analysis

The haptophyceae Phaeocystis globosa is a species responsible for harmful algal blooms in the global ocean, forming blooms in the Beibu Gulf annually since 2011. This species can alternate between solitary free-living cells and colonies. Colonies are the dominant morphotype during blooms. To date, the underlying mechanism of P. globosa blooms in the Beibu Gulf is poorly understood. After combining results of ecological surveys, laboratory studies, and metatranscriptome and bioinformatics analyses, it was found that low temperatures, high nitrate, and low organic phosphorus induced P. globosa blooms in the Beibu Gulf. Additionally, the unique genetic and physiological characteristics that allow P. globosa to stand out as a dominant species in such an environment include (1) several genes encoding high-affinity nitrate transport proteins that could be highly expressed under sufficient nitrate conditions; (2) energy metabolism genes involved in photosynthesis and oxidative phosphorylation that were actively expressed at low temperatures to carry out carbon and energy reversion and produce sufficient ATP for various life activities, individually; (3) abundant glycan synthesis genes that were highly expressed at low temperatures, thus synthesizing large quantities of proteoglycans to construct the mucilaginous envelope forming the colony; (4) cells in colonies exhibited active gene expression in DNA replication contributing to a faster growth rate, which could help P. globosa occupy niches quickly; and (5) the energy and material expenditure was redistributed in colonial cells accompanied with chitin filaments and flagella degraded, more expenditure was used for the synthesis of the mucilaginous envelope and the rapid proliferation.

Establishment of a multiplex polymerase chain reaction detection assay for three common harmful microalgae in the East China Sea

It is urgent to develop techniques that can simultaneously detect multiple microalgae, due to the diversity of harmful algal blooms (HABs)-forming algal species. The target algae species in this study are Heterosigma akashiwo, Prorocentrum donghaiense and Karenia mikimotoi. These algae are the dominant species that cause HABs in the East China Sea, and the multiple detection technique focusing on these three algae is not common. Therefore, this study established a multiplex polymerase chain reaction(mPCR) to diagnose the three algae, which is simple and low cost. First, the corresponding specific primers were designed based on the D1-D2 region of the large subunit (LSU) ribosomal DNA sequence. Then, mPCR was established and the reaction conditions were optimized. And the specificity, sensitivity, and stability of mPCR were evaluated. The result of specificity test showed that the established mPCR had good specificity for the target microalgae and did not cross-react with eighteen non-target microalgae. The sensitivity of experiment was 3.3 × 10^-1 ng μL^-1, and the established mPCR was not affected by the interfering microalgae. Moreover, the practicability evaluation of mPCR by using the simulated natural water samples showed that the detection limit of target microalgae was 100 cells mL^-1, which could meet the demand for early warning of HABs. In summary, the established mPCR is characterized by strong specificity, good stability, and multiple analysis to detect H. akashiwo, P. donghaiense, and K. mikimotoi.

ROS-dependent cell death of Heterosigma akashiwo induced by algicidal bacterium Hahella sp. KA22

Marine algicidal bacteria and their metabolites are considered to be one of the most effective strategies to mitigate the harmful algal blooms (HABs). The bacterium Hahella sp. KA22 has previously been confirmed to have strong algicidal activity against the HABs causing microalgae, Heterosigma akashiwo. In this study, the molecular mechanism of microalgae cell death was detected. The results showed that the cell growth rate and photosynthetic efficiency were inhibited with addition of algicidal strain KA22, while the accumulation of reactive oxygen species (ROS) and oxidative damage in H. akashiwo cells increased. A total of 2056 unigenes were recognized to be differentially expressed in transcriptome sequences. In particular, the transcriptional levels of light-harvesting pigments and structural proteins in the oxygen-evolving-complex were continuously down-regulated, corresponding to the significant reduction of photosynthetic efficiency and the accumulation of ROS. Furthermore, glutamate dehydrogenase was significantly up-regulated in abundance. Meanwhile, calcium-dependent protein kinases were also detected with significant changes. Collectively, algicidal stress caused the suppressed electron transfer in chloroplast and impaired detoxification of intracellular oxidants by glutathione, which may subsequently result in multiple cell regulation and metabolic responses and ultimately lead to the ROS-dependent cell death of H. akashiwo.

Ocean warming and acidification affect the transitional C:N:P ratio and macromolecular accumulation in the harmful raphidophyte Heterosigma akashiwo

Despite an increase in ocean warming and acidification that is expected to increase the number of harmful algal species worldwide, the population of the raphidophyte Heterosigma akashiwo has been reported to be reduced. However, how this species physically and metabolically modifies transitional C:N:P ratio and macromolecule accumulation is unknown. Considering 1st, 10th, and 20th culture generations under present (low-temperature; low-carbon-dioxide [LTLC] 21 °C; pCO₂ 400 ppm) and future (high-temperature; high-carbon-dioxide [HTHC] 25 °C; pCO₂ 1000 ppm) ocean conditions, we examined transitional C:N:P ratio and macromolecule level changes and performed transcriptome sequencing. The results showed that compared to 1st generation cells, 20th generation cells under HTHC conditions showed a large decrease in carbon quota (Q_C: 34%), nitrogen quota (Q_N: 36%), and phosphorus quota (Q_P: 32%), which were reflected in an overall reduction in DNA and RNA quantity. Decreased activation of photosynthetic, carbon fixation and lipid metabolic pathways coincided with changes in photosynthetic efficiency, carbon concentration, and lipid accumulation after long-term (20th generation) exposure to HTHC conditions. We observed that these variations in internal metabolic pathways were caused by external changes in temperature, which activated the (Ca⁺) signaling pathway, and external changes in pCO₂, which altered proton exchange pathways. Our results suggest that H. akashiwo in a temperate environment will undergo profound changes in C:N:P ratio and macromolecular properties, leading to programmed cell death, in the future.

Full-Length Transcriptome Analysis of the Ichthyotoxic Harmful Alga Heterosigma akashiwo (Raphidophyceae) Using Single-Molecule Real-Time Sequencing

The raphidophyte Heterosigma akashiwo is a harmful algal species. The bloom of this organism has been associated with the massive mortality of fish in many coastal waters. To investigate the molecular mechanism of H. akashiwo blooms, having a reliable reference transcriptome of this species is essential. Therefore, in this study, a full-length transcriptome of H. akashiwo was obtained by single-molecule real-time sequencing. In total, 45.44 Gb subread bases were generated, and 16,668 unigenes were obtained after the sequencing data processing. A total of 8666 (52.00%) unigenes were successfully annotated using seven public databases. Among them, mostly phosphorus and nitrogen metabolism genes were detected. Moreover, there were 300 putative transcription factors, 4392 putative long non-coding RNAs, and 7851 simple sequence repeats predicted. This study provides a valuable reference transcriptome for understanding how H. akashiwo blooms at a molecular level.

Eukaryotic plankton community assembly and influencing factors between continental shelf and slope sites in the northern South China Sea

Eukaryotic plankton are pivotal members of marine ecosystems playing crucial roles in marine food webs and biogeochemical cycles. However, understanding the patterns and drivers of their community assembly remains a grand challenge. A study was conducted in the northern South China Sea (SCS) to address this issue. Here, 49 samples were collected and size-fractionated from discrete depths at continental shelf and continental slope in the northern SCS over a diel cycle. From high throughput sequencing of the 18S rDNA gene V4 region, 2463 operational taxonomic units (OTUs) were retrieved. Alveolata and Opisthokonta overwhelmingly dominated the assemblages in the abundance (44.76%, 31.08%) and species richness (59%, 12%). Biodiversity was higher in the slope than the shelf and increased with depth. Temperature and salinity appeared to be the most important deterministic drivers of taxon composition. Community structure was influenced by multiple factors in the importance order of: environmental factors (temperature + salinity) > spatial factor > water depth > sampling time. Furthermore, the neutral model explained more variations in the smaller-sized (0.22-3 μm) community (24%) than larger-sized (3-200 μm) community (16%) but generally explained less variations than did deterministic processes. Additionally, our data indicated that the larger plankton might be more environmentally filtered and less plastic whereas the smaller plankton had stronger dispersal ability. This study sheds light on the differential contributions of the deterministic process and stochastic process and complexities of assembly mechanisms in shaping the community assembly of micro-nano and pico-eukaryotic biospheres in a subtropical ocean.

Strong population genomic structure of the toxic dinoflagellate Alexandrium minutum inferred from meta-transcriptome samples

Despite theoretical expectations, marine microeukaryote population are often highly structured and the mechanisms behind such patterns remain to be elucidated. These organisms display huge census population sizes, yet genotyping usually requires clonal strains originating from single cells, hindering proper population sampling. Estimating allelic frequency directly from population wide samples, without any isolation step, offers an interesting alternative. Here, we validate the use of meta-transcriptome environmental samples to determine the population genetic structure of the dinoflagellate Alexandrium minutum. Strain and meta-transcriptome based results both indicated a strong genetic structure for A. minutum in Western Europe, to the level expected between cryptic species. The presence of numerous private alleles, and even fixed polymorphism, would indicate ancient divergence and absence of gene flow between populations. Single nucleotide polymorphisms (SNPs) displaying strong allele frequency differences were distributed throughout the genome, which might indicate pervasive selection from standing genetic variation (soft selective sweeps). However, a few genomic regions displayed extremely low diversity that could result from the fixation of adaptive de novo mutations (hard selective sweeps) within the populations.

Active reconfiguration of cytoplasmic lipid droplets governs migration of nutrient-limited phytoplankton

Nutrient availability, along with light and temperature, drives marine primary production. The ability to migrate vertically, a critical trait of motile phytoplankton, allows species to optimize nutrient uptake, storage, and growth. However, this traditional view discounts the possibility that migration in nutrient-limited waters may be actively modulated by the emergence of energy-storing organelles. Here, we report that bloom-forming raphidophytes harness energy-storing cytoplasmic lipid droplets (LDs) to biomechanically regulate vertical migration in nutrient-limited settings. LDs grow and translocate directionally within the cytoplasm, steering strain-specific shifts in the speed, trajectory, and stability of swimming cells. Nutrient reincorporation restores their swimming traits, mediated by an active reconfiguration of LD size and coordinates. A mathematical model of cell mechanics establishes the mechanistic coupling between intracellular changes and emergent migratory behavior. Amenable to the associated photophysiology, LD-governed behavioral shift highlights an exquisite microbial strategy toward niche expansion and resource optimization in nutrient-limited oceans.

Transcriptomic response of the harmful algae Heterosigma akashiwo to polyphosphate utilization and phosphate stress

Phosphorus (P) is one of the major macronutrients necessary for phytoplankton growth. In some parts of the ocean, however, P is frequently scarce, hence, there is limited phytoplankton growth. To cope with P deficiency, phytoplankton evolved a variety of strategies, including, utilization of different P sources. Polyphosphate (polyP) is ubiquitously present and serves an essential function in aquatic environments, but it is unclear if and how this polymer is utilized by phytoplankton. Here, we examined the physiological and molecular responses of the widely present harmful algal bloom (HAB) species, Heterosigma akashiwo in polyP utilization, and in coping with P-deficiency. Our results revealed that two forms of inorganic polyP, namely, sodium tripolyphosphate and sodium hexametaphosphate, support H. akashiwo growth as efficiently as orthophosphate. However, few genes involved in polyP utilization have been identified. Under P-deficient conditions, genes associated with P transport, dissolved organic P utilization, sulfolipid synthesis, and energy production, were markedly elevated. In summary, our results indicate that polyP is bioavailable to H. akashiwo, and this HAB species have evolved a comprehensive strategy to cope with P deficiency.

Isolation and characterization of an algicidal bacterium against the bloom-forming algae raphidophyte Heterosigma akashiwo

Harmful algae blooms (HABs) have increased in intensity and frequency worldwide, causing negative effects on public health and marine ecosystems. This study isolated and identified the bloom causing species and its associated algicidal bacterium during a phytoplankton bloom in coastal waters of Lianyungang, China. Morphological observations and DNA barcoding analysis indicate that the studied phytoplankton bloom was caused by the raphidophyte Heterosigma akashiwo, and the algicidal bacterium, strain LD-B1, was identified as a species belonging to the genus Pseudoalteromonas. Furthermore, the algicidal effects of strain LD-B1 against H. akashiwo were characterized; revealing strain LD-B1 show strong algicidal activity against H. akashiwo. After 48 h of bacterium culture addition, the algicidal rate reached up to 98.8% with a 1% final volume rate. Moreover, our findings indicate strain LD-B1's extracellular compounds involved in algicidal activity are likely not proteinaceous. These findings indicate that the isolated strain, LD-B1, is a promising algicidal bacterium to control H. akashiwo blooms.

Diversity and Evolution of Mamiellophyceae: Early-Diverging Phytoplanktonic Green Algae Containing Many Cosmopolitan Species

The genomic revolution has bridged a gap in our knowledge about the diversity, biology and evolution of unicellular photosynthetic eukaryotes, which bear very few discriminating morphological features among species from the same genus. The high-quality genome resources available in the class Mamiellophyceae (Chlorophyta) have been paramount to estimate species diversity and screen available metagenomic data to assess the biogeography and ecological niches of different species on a global scale. Here we review the current knowledge about the diversity, ecology and evolution of the Mamiellophyceae and the large double-stranded DNA prasinoviruses infecting them, brought by the combination of genomic and metagenomic analyses, including 26 metabarcoding environmental studies, as well as the pan-oceanic GOS and the Tara Oceans expeditions.

Emerging harmful algal bloom species over the last four decades in China

The first recorded micro-algae bloom in Chinese coastal waters dates back to 1933 and was caused by a mixture of Noctiluca scintillans and Skeletonema costatum sensu lato along the Zhejiang coast (the East China Sea). While well-documented harmful algal blooms (HABs) appeared to be extremely scarce from the 1950s to 1990, both the frequency and intensity have been reportedly increasing since 1990. Among them, the fish-killing HABs, mainly caused by Karenia mikimotoi, Karlodinium digitatum, Karlodinium veneficum, Margalefidinium polykrikoides, and Heterocapsa spp., have intensified. Karenia mikimotoi was responsible for at least two extremely serious events in the Pearl River Estuary in 1998 and the Taiwan Strait (in the East China Sea) in 2012, which appeared to be associated with abnormal climate conditions and excessive nutrients loading. Other major toxic algal blooms have been caused by the species responsible for paralytic shellfish poisoning (including Alexandrium catenella, Alexandrium pacificum, Gymnodinium catenatum) and diarrhetic shellfish poisoning (including Dinophysis spp., and a couple of benthic dinoflagellates). Consequent closures of shellfish farms have resulted in enormous economic losses, while consumption of contaminated shellfish has led to occasional human mortality in the Bohai Sea and the East China Sea. Expansions of these HABs species along the coastline of China have occurred over the last four decades and, due to the projected global changes in the climate and marine environments and other anthropological activities, there is potential for the emergence of new types of HABs in China in the future. This literature review aimed to present an updated overview of HABs species over the last four decades in China.

Toxicological effects of hypoxanthine on Heterosigmaakashiwo: Mechanism of growth inhibition and change in hemolytic toxin content

Heterosigmaakashiwo is an algal species that causes harmful algal blooms (HABs) with strong hemolytic toxicity on coastal aquatic organisms. This study investigated the mechanism of growth inhibition and changes in hemolytic toxin contents in algal culture after exposure to hypoxanthine, a compound secreted by algicidal bacterium Bacillus sp.strain B1. An algal inhibition rate of 86% was observed with 1.0 mM hypoxanthine treatment on day 15. The levels of superoxide dismutase and catalase in algal cell culture increased while that of glutathione decreased during the treatment. In addition,the level of hemolytic toxin contents increased on day 3 under hypoxanthine treatment, and significantly decreased on days 6, 9, 12, and 15. Twelve fatty acids in H.akashiwo were detected by GC-MS, and the changes in the contents of C16, C18, C18:4ω3, and C20:5ω3 in the treatment group were consistent with the change in hemolytic toxin content. The four fatty acids were tested for hemolysis and it was observed that the hemolysis rate of 25 μg/mL C18:4ω3 and 5 μg/mL C20:5ω3 reached more than 80%, but C16 and C18 exhibited no hemolytic capability.Therefore, our results showed that hypoxanthine inhibited the growth of H. akashiwo through the changes of levels of antioxidants and hemolytic toxin content in the cultures, and fatty acids C18:4ω3 and C20:5ω3 were contributors to hemolytic toxins. The results confirmed that hypoxanthine is a potential algal inhibitor to prevent HABs.

Growth and Nutrient Uptake Characteristics of Heterosigma akashiwo (Raphidophyceae) under Nitrogen and Phosphorus Concentrations in the East China Sea

Heterosigma akashiwo is classified as a harmful algal bloom (HAB) species that frequently occurs in eutrophic coastal waters and results in the contamination and mortality of fish and shellfish. The growth of H. akashiwo in four phosphate and nitrate concentration scenarios, representing the observed nutrient concentration ranges in the East China Sea (ECS), was evaluated to further understand the effect of nutrient concentrations on H. akashiwo blooms. The specific growth rate in the exponential growth phase (µ′) and the maximum cell density were lower (17–21% and 41%, respectively) under low phosphorus concentration scenarios, compared to the rates observed under high phosphorus concentration scenarios. The cellular nitrogen-to-phosphorus ratios of H. akashiwo were influenced by the initially supplied N:P ratio and the allocation strategy employed. Phosphorus concentration had a greater influence on the total growth of H. akashiwo than nitrate did, within the natural nutrient conditions of the ECS. These results could serve as a reference for coastal water management and marine ecological management and may be useful for further studies on the simulation and prediction of H. akashiwo blooms, particularly in the ECS.

Metabarcoding dissection of harmful algal bloom species in the East China Sea off Southern Zhejiang Province in late spring

The coastal region of the East China Sea (ECS) is a famous "hotspot" for harmful algal blooms (HABs) in China. We hypothesize that such frequent occurrences of diverse HABs in the ECS are determined by the presence of unique HAB species in this region. In this project, we identified 3966 amplicon sequence variants (ASVs) representing 35 classes in six protist phyla/divisions. Among the 237 annotated protist species, we identified 58 HAB species, of which 23 HAB species had never been previously reported in the ECS. Many HAB species also displayed unique spatial distribution patterns in the ECS. Notably, we identified three HAB species Prorocentrum donghaiense, Lebouridinium glaucum and Noctiluca scintillans in the site S05-1 with substantially elevated abundance, suggesting that this sampling site was experiencing a multiple-species HAB event. This study was the first attempt in applying ASV-based metabarcoding analysis in studying protist and HAB species in the ECS.

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.

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.

Transcriptomic and isotopic data reveal central role of ammonium in facilitating the growth of the mixotrophic dinoflagellate, Dinophysis acuminata

Dinophysis spp. are mixotrophs that are dependent on specific prey, but are also potentially reliant on dissolved nutrients. The extent to which Dinophysis relies on exogenous N and the specific biochemical pathways important for supporting its autotrophic and heterotrophic growth are unknown. Here, the nutritional ecology of Dinophysis was explored using two approaches: 1) ¹⁵N tracer experiments were conducted to quantify the concentration-dependent uptake rates and associated kinetics of various N compounds (nitrate, ammonium, urea) of Dinophysis cultures and 2) the transcriptomic responses of Dinophysis cultures grown with multiple combinations of prey and nutrients were assessed via dinoflagellate spliced leader-based transcriptome profiling. Of the N compounds examined, ammonium had the highest V_max and affinity coefficient, and lowest K_s for both pre-starved and pre-fed cultures, collectively demonstrating the preference of Dinophysis for this N source while little-to-no nitrate uptake was observed. During the transcriptome experiments, Dinophysis grown with nitrate and without prey had the largest number of genes with lower transcript abundances, did not increase abundance of transcripts associated with nitrate/nitrite uptake or reduction, and displayed no cellular growth, suggesting D. acuminata is not capable of growing on nitrate. When offered prey, the transcriptomic response of Dinophysis included the production of phagolysosomes, enzymes involved in protein and lipid catabolism, and N acquisition through amino acid degradation pathways. Compared with cultures only offered ammonium or prey, cultures offered both ammonium and prey had the largest number of genes with increased transcript abundances, the highest growth rate, and the unique activation of multiple pathways involved in cellular catabolism, further evidencing the ability of Dinophysis to grow optimally as a mixotroph. Collectively, this study evidences the key role ammonium plays in the mixotrophic growth of Dinophysis and reveals the precise biochemical pathways that facilitate its mixotrophic growth.

A metatranscriptomic analysis of changing dynamics in the plankton communities adjacent to aquaculture leases in southern Tasmania, Australia

Aquaculture releases nitrogen to the marine environment, potentially changing dynamics of local plankton populations and causing adverse impacts. Metatranscriptomics have been used to study planktonic nutrient cycles and community dynamics. We hypothesised that the metatranscriptome could be used to monitor changing phytoplankton physiology near leases. To test this hypothesis, opportunistic samples were collected from one oceanic location in winter and one estuarine location in spring and analysed via RNASeq. Transcriptomes from different locations were found to have little overlap, due to different community compositions in the oceanic and estuarine locations. Transcript function was similar at each location. Proximity to the salmon pen had little influence over the transcriptome at the estuarine location. In the oceanic environment, diatom-based activity decreased near pens and dinoflagellate-based activity increased as demonstrated through the abundance of carbon fixation and nitrogen-acquisition-related transcripts. Our initial results suggest that the use of the metatranscriptome in monitoring is promising.

Metatranscriptomics From a Small Aquatic System: Microeukaryotic Community Functions Through the Diurnal Cycle

Light is an important factor for the growth of planktonic organisms, and many of them depend on the diurnal light/dark cycle to regulate key metabolic processes. So far, most of the diel responses were only studied in single species or marine and large lake communities. Yet, we lack information on whether these processes are regulated similarly in small aquatic systems such as ponds. Here, we investigated the activity of a microeukaryotic community from a temperate, small freshwater pond in response to the diurnal cycle. For this, we took samples at midday and night during the Central European summer. We extracted pigments and RNA from samples and the sequencing of eukaryotic transcripts allowed us to obtain day and night metatranscriptomes. Differentially expressed transcripts primarily corresponded to photosynthesis-related and translational processes, and were found to be upregulated at midday with high light conditions compared to darkness. Unique gene ontology classes were found at each respective condition. During the day, ontology classes including photoreception for photosynthesis, defense, and stress mechanisms dominated, while motility, ribosomal assembly and other large, energy-consuming processes were restricted to the night. Euglenophyta and Chlorophyta dominated the active phototrophic community, as shown by the pigment composition analysis. Regarding the gene expression patterns, we could confirm that the pond community appears to follow similar diurnal dynamics as those described for larger aquatic ecosystems. Overall, combining pigment analyses, metatranscriptomics, and data on physicochemical factors yielded considerably more insight into the metabolic processes performed by the microeukaryotic community of a small freshwater ecosystem.

Species specific gene expression dynamics during harmful algal blooms

Harmful algal blooms are caused by specific members of microbial communities. Understanding the dynamics of these events requires comparing the strategies developed by the problematic species to cope with environmental fluctuations to the ones developed by the other members of the community. During three consecutive years, the meta-transcriptome of micro-eukaryote communities was sequenced during blooms of the toxic dinoflagellate Alexandrium minutum. The dataset was analyzed to investigate species specific gene expression dynamics. Major shifts in gene expression were explained by the succession of different species within the community. Although expression patterns were strongly correlated with fluctuation of the abiotic environment, and more specifically with nutrient concentration, transcripts specifically involved in nutrient uptake and metabolism did not display extensive changes in gene expression. Compared to the other members of the community, A. minutum displayed a very specific expression pattern, with lower expression of photosynthesis transcripts and central metabolism genes (TCA cycle, glucose metabolism, glycolysis…) and contrasting expression pattern of ion transporters across environmental conditions. These results suggest the importance of mixotrophy, cell motility and cell-to-cell interactions during A. minutum blooms.

Utilization of various forms of nitrogen and expression regulation of transporters in the harmful alga Heterosigma akashiwo (Raphidophyceae)

Nitrogen (N) is an essential nutrient for phytoplankton growth. There is ample evidence that N enrichment promotes harmful algae blooms (HABs) but molecular mechanisms regulating N-nutrient uptake and metabolism are not so clear, especially for the raphidophyte Heterosigma akashiwo, which forms ichthyotoxic HABs in many coastal waters. In this study, the utilization of three different chemical forms of N (nitrate, ammonium, and urea) by H. akashiwo CCMA 369 was investigated in batch culture conditions. Results showed that H. akashiwo grew well on all three N compounds, and the highest cell yield occurred in the NH₄⁺ culture group. Reverse transcription quantitative PCR analysis revealed that the expression of high-affinity NO₃^- transporter (NIT), NH₄⁺ transporter (AMT) and high-affinity urea active transporter (DUR3), were significantly up-regulated under N-limitation compared to the N-replete control. The mRNA levels of AMT and DUR3 also displayed a clear diel rhythm, with high levels at midnight. In addition, NH₄⁺ addition (5 μM) did not depress the transcript abundance of any of the three N transporters. Compared with the co-occurring immobile diatom Skeletonema costatum, the high expression of AMT in dark period in H. akashiwo is consistent with its diel vertical migration behavior, which may promote N-nutrient acquisition from deeper layers and give advantages for H. akashiwo to form blooms.

Identification and expression analysis of meiosis-related genes in the harmful alga Heterosigma akashiwo (Raphidophyceae)

Sexual reproduction (SR), hallmarked by meiosis, is widespread in eukaryotes. In phytoplankton, SR has been observed in many lineages, but molecular information on SR or meiosis of harmful algal bloom (HAB) species is scarce. The raphidophyte Heterosigma akashiwo is a globally distributed and devastating HAB species, but molecular evidence of its SR or meiosis is lacking. Here, to address the gap of knowledge, the presence of meiotic genes in H. akashiwo were examined. Interestingly, seven meiosis-specific or related genes (SPO11, MND, RAD21, RAD51, MSH2, MSH6 and MEI2) were identified from H. akashiwo transcriptomes. Furthermore, expression patterns of these genes except MSH6 (excluded due to primer failure) were investigated using quantitative reverse-transcription PCR. Results showed that the examined genes exhibited similar diel rhythms, typically, highest in early dark period and then gradually decreasing until mid-day. In addition, the expression of these six genes was not higher in the stationary growth stage than in the exponential stage, as would be expected if meiosis was to form cysts, and their elevated expression in response to colchicine treatment (arresting cells in the G2/M transition) indicated a potential role of these genes in cell division and population growth in H. akashiwo. Consistent with this, we also found a morning to afternoon shift in the expression of these genes during the bloom of H. akashiwo. This study documents a part of the typical SR gene repertoire and its potential role in regulating cell division in H. akashiwo, offering candidates for population growth markers for bloom monitoring although its linkage to meiosis and SR remain to be studied further in the future.

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.

Dynamic CO2 and pH levels in coastal, estuarine, and inland waters: Theoretical and observed effects on harmful algal blooms

Rising concentrations of atmospheric CO₂ results in higher equilibrium concentrations of dissolved CO₂ in natural waters, with corresponding increases in hydrogen ion and bicarbonate concentrations and decreases in hydroxyl ion and carbonate concentrations. Superimposed on these climate change effects is the dynamic nature of carbon cycling in coastal zones, which can lead to seasonal and diel changes in pH and CO₂ concentrations that can exceed changes expected for open ocean ecosystems by the end of the century. Among harmful algae, i.e. some species and/or strains of Cyanobacteria, Dinophyceae, Prymnesiophyceae, Bacillariophyceae, and Ulvophyceae, the occurrence of a CO₂ concentrating mechanisms (CCMs) is the most frequent mechanism of inorganic carbon acquisition in natural waters in equilibrium with the present atmosphere (400 μmol CO₂  mol^-1 total gas), with varying phenotypic modification of the CCM. No data on CCMs are available for Raphidophyceae or the brown tide Pelagophyceae. Several HAB species and/or strains respond to increased CO₂ concentrations with increases in growth rate and/or cellular toxin content, however, others are unaffected. Beyond the effects of altered C concentrations and speciation on HABs, changes in pH in natural waters are likely to have profound effects on algal physiology. This review outlines the implications of changes in inorganic cycling for HABs in coastal zones, and reviews the knowns and unknowns with regard to how HABs can be expected to ocean acidification. We further point to the large regions of uncertainty with regard to this evolving field.

Formation mechanism of the Microcystis aeruginosa bloom in the water with low dissolved phosphorus

The utilization of phosphorus by algae in the low-phosphorus state has drawn wide concerns due to the high risk of forming algal blooms. The cyanobacteria Microcystis aeruginosa (M. aeruginosa) grew well under low-phosphorus condition by hydrolyzing dissolved organic phosphorus (DOP) to dissolved inorganic phosphorus (DIP) through alkaline phosphatase (AP). There was a negative correlation between DIP concentration and AP activity of algae. AP activity significantly increased at 0-3 d (p < 0.05), and reached the peak values of 43.06 and 49.11 King unit/gprot on day 5 for DIP (0.1 mg/L) and DOP (4.0 mg/L), respectively. The relative expression of phosphate transporter gene increased with decreasing phosphorus concentrations. The catalase activity under low-phosphorus condition increased significantly (p < 0.05) after one week, and was generally higher than 0.15 U/mgprot on day 14. Understanding the utilization efficiency and mechanism of DIP and DOP in the low-phosphorus state would help to inhibit the formation of algal blooms.

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.

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.

Microbial Community Structure and Associations During a Marine Dinoflagellate Bloom

Interactions between microorganisms and algae during bloom events significantly impacts their physiology, alters ambient chemistry, and shapes ecosystem diversity. The potential role these interactions have in bloom development and decline are also of particular interest given the ecosystem impacts of algal blooms. We hypothesized that microbial community structure and succession is linked to specific bloom stages, and reflects complex interactions among taxa comprising the phycosphere environment. This investigation used pyrosequencing and correlation approaches to assess patterns and associations among bacteria, archaea, and microeukaryotes during a spring bloom of the dinoflagellate Alexandrium catenella. Within the bacterial community, Gammaproteobacteria and Bacteroidetes were predominant during the initial bloom stage, while Alphaproteobacteria, Cyanobacteria, and Actinobacteria were the most abundant taxa present during bloom onset and termination. In the archaea biosphere, methanogenic members were present during the early bloom period while the majority of species identified in the late bloom stage were ammonia-oxidizing archaea and Halobacteriales. Dinoflagellates were the major eukaryotic group present during most stages of the bloom, whereas a mixed assemblage comprising diatoms, green-algae, rotifera, and other microzooplankton were present during bloom termination. Temperature and salinity were key environmental factors associated with changes in bacterial and archaeal community structure, respectively, whereas inorganic nitrogen and inorganic phosphate were associated with eukaryotic variation. The relative contribution of environmental parameters measured during the bloom to variability among samples was 35.3%. Interaction analysis showed that Maxillopoda, Spirotrichea, Dinoflagellata, and Halobacteria were keystone taxa within the positive-correlation network, while Halobacteria, Dictyochophyceae, Mamiellophyceae, and Gammaproteobacteria were the main contributors to the negative-correlation network. The positive and negative relationships were the primary drivers of mutualist and competitive interactions that impacted algal bloom fate, respectively. Functional predictions showed that blooms enhance microbial carbohydrate and energy metabolism, and alter the sulfur cycle. Our results suggest that microbial community structure is strongly linked to bloom progression, although specific drivers of community interactions and responses are not well understood. The importance of considering biotic interactions (e.g., competition, symbiosis, and predation) when investigating the link between microbial ecological behavior and an algal bloom's trajectory is also highlighted.