Phytoplankton community and HAB species in the South China Sea detected by morphological and metabarcoding approaches

Early warning on the potential harmful algal bloom species in Beibu Gulf of South China Sea under the background of climate change and human activity

Human activity and global climate change increasingly affect marine environments, leading to increases in harmful algal blooms (HABs) caused by phytoplankton. These blooms pose significant threats to public health, tourism, fisheries, and ecosystems. As an important fishing ground and tourist destination, the Beibu Gulf faces growing environmental pressure. This study sought to assess the phytoplankton community structure and status of HABs, with a focus on potential HAB species. Using environmental DNA (eDNA) metabarcoding, summer and winter surveys at both coastal and offshore waters revealed 66 potential HAB species, 23 of which were newly recorded in the Beibu Gulf. The potential HAB species exhibited greater richness and relative abundance in summer than in winter. Offshore areas showed greater diversity, whereas coastal areas showed greater relative abundance. Temperature emerged as the most influential factor shaping phytoplankton composition, and pH was found to play an important role in coastal areas. Nutrients such as silicate and ammonium are critical for the distribution of potential HAB species. Among the potential HAB species, Cyclotella cryptica predominated in coastal areas during winter, whereas Chaetoceros tenuissimus predominated in summer. Some species that caused severe HAB events in other oceanic regions were first detected in this study, including Margalefidinium polykrikoides, Karlodinium veneficum, and Prorocentrum concavum. This study revealed the diversity and complexity of the phytoplankton community in the Beibu Gulf, emphasizing the critical importance of monitoring and early warning of potential HAB species, particularly those driven by human activities and climate change.

The combination of microscopic and metabarcoding analyses revealed a high diversity of algal bloom species in the Taiwan Strait

The Taiwan Strait is located between Fujian and Taiwan Provinces, China, and is susceptible to various currents and human activities, which result in frequent algal blooms. Surface water samples were collected in the Taiwan Strait in the spring of 2022, and phytoplankton community structure was assessed by both microscopic and metabarcoding analyses. Diatoms dominated in phytoplankton community based on microscopy, while dinoflagellates dominated in the eukaryotic phytoplankton based on metabarcoding. Water temperature was the main environmental factor affecting the phytoplankton community structure. A total of 133 algal bloom (AB) species, including 32 harmful algal bloom (HAB) species, were detected via the combination of both methods. Scrippsiella acuminata and Noctiluca scintillans were the dominant AB species. N. scintillans occurs abundantly in shallow coastal waters with low temperatures, low salinities and high nutrient levels. The distribution of AB species in our study is consistent with historical algal blooms in the Taiwan Strait.

Spatiotemporal distribution, composition, and influencing factors of lipophilic marine algal toxins in the Beibu Gulf, China: Implications for necessity of offshore seawater monitoring

Lipophilic marine algal toxins (LMATs) are persistent, bioaccumulative metabolites that pose significant threats to marine biodiversity and seafood safety. To comprehensively understand LMATs pollution in offshore natural fishing grounds, we investigated the spatiotemporal distribution and composition of 13 LMATs in surface seawater in Beibu Gulf's offshore area, along with their key influencing factors and latitudinal variations across coastal regions in China. We detected 20-methyl spirolide G (SPX20G) for the first time in coastal waters in China, alongside other frequently detected LMATs, including okadaic acid (OA), dinophysistoxin-1 (DTX1), pectenotoxin-2 (PTX2), azaspiracid-1 (AZA1), azaspiracid-2 (AZA2), gymnodimine (GYM), and 13-desmethyl spirolide C (SPX1). Significant seasonal variations were observed, with summer exhibiting eight LMAT types and higher total LMAT concentrations (ΣLMATs) (16.92 ± 20.45 ng/L) compared to winter's five LMAT types and ΣLMATs concentrations at 5.63 ± 3.30 ng/L. Predominant PTX2 and OA concentrations showed a decreasing trend from northern to southern Beibu Gulf, while other LMATs were distributed in patches, particularly in summer. Notably, cyclic imines (GYM and SPX1) were detected with high frequencies and concentrations, underscoring the need for their monitoring. Spearman correlation and redundancy analyses identified water depth, chlorophyll a (Chla), temperature, salinity, nutrients (DIN:DIP and SiO₃²⁻), and dissolved oxygen (DO) as key factors influencing LMATs distribution in the Beibu Gulf. Latitudinal analysis of LMAT contents across Chinese waters revealed DTX1 as more prevalent at higher latitudes, and GYM, SPX1, and AZA1 more common at lower latitudes, while ΣLMATs concentrations from this study ranked at moderate to high levels compared to other regions. Therefore, special attention is required for LMATs pollution in offshore fishing areas.

Kareniaceae in focus: A molecular survey of harmful algal dinoflagellates in the South China Sea

In recent years, species of the family Kareniaceae have been forming algal blooms worldwide, leading to significant fish mortality. The biodiversity of this group of dinoflagellates has been underestimated, primarily due to their highly similar morphologies. In this study, using high-throughput sequencing (HTS) targeted the large subunit ribosomal DNA (LSU rDNA) of dinoflagellates, 21 Kareniaceae species were identified from the South China Sea (SCS), including nine species newly recorded in Chinese coastal waters. Additionally, species-specific real-time fluorescent quantitative PCR (qPCR) assays targeting the internal transcribed spacer (ITS) for three Karenia species (Karenia longicanalis, K. papilionacea, and K. selliformis) isolated from the SCS were established. Karenia longicanalis and K. papilionacea were identified as the most prevalent Karenia species in summer in the Zhujiang Estuary via qPCR. Further qPCR analysis during a Karenia spp. bloom in the Beibu Gulf on August 3, 2023, revealed K. selliformis as the most abundant species, followed by K. longicanalis, K. mikimotoi, and K. papilionacea. Both HTS and qPCR methods successfully detected the targeted Karenia species in field samples; however, qPCR assays identified them at more stations than HTS. The integration of qPCR assays with HTS provides a comprehensive and precise assessment of the diversity of Kareniaceae species and harmful dinoflagellates. This multi-method approach enhances our understanding and management of harmful algal blooms, thereby contributing to the improved protection of marine ecosystems.

Performance of DNA Metabarcoding vs. Morphological Methods for Assessing Intertidal Turf and Foliose Algae Diversity

Large biogeographical shifts in marine communities are taking place in response to climate change and biological invasions yet we still lack a full understanding of their diversity and distribution. An important example of this is turf and foliose algae that are key coastal primary producers in several regions and are expanding into new environments. Traditionally, monitoring turf and foliose algae communities involves species identification based on morphological traits, which is challenging due to their reduced dimensions and highly variable morphology. Molecular methods promise to revolutionise this field, but their effectiveness in detecting turf and foliose algae has yet to be tested. Here, we evaluate the performance of DNA metabarcoding (COI and rbcL markers) and morphological identification (in situ and photoquadrat) to describe intertidal turf and foliose algae communities along the Portuguese coast. Both molecular markers detected more taxa than the morphological methods and showed greater discrimination of turf and foliose algae communities between regions, matching our knowledge of the geographical and climatic patterns for the region. In sum, our multi-marker metabarcoding approach was more efficient than morphology-based methods in characterising turf and foliose algae communities along the Portuguese coast, differentiating morphologically similar species, and detecting unicellular organisms. However, certain taxa that were identified by in situ and photoquadrat approaches were not detected through metabarcoding, partly due to lack of reference barcodes or taxonomic resolution. Metabarcoding emerges as a valuable tool for monitoring these communities, particularly in long-term programmes requiring accuracy, speed, and reproducibility.

Morphological and Molecular Diversity of Phytoplankton in Beibu Gulf, Northern South China Sea

The Beibu Gulf, a vital region for marine biodiversity and aquaculture, is increasingly affected by nutrient-driven ecological shifts in the phytoplankton community. This study combined morphology and eDNA metabarcoding (18S rDNA V4) to investigate phytoplankton diversity and environmental drivers during summer and winter in the Beibu Gulf. Metabarcoding detected 3.5 times more phytoplankton species, contributing to higher species diversity and richness than morphology. Metabarcoding identified 200 phytoplankton genera from eight phyla, while morphology only identified 49 genera from six phyla. Both methods revealed different dominant phytoplankton communities. Bacillariophyta and Haptophyta dominated the phytoplankton community based on morphology, in summer and winter, respectively; meanwhile, Dinophyta dominated in both seasons under metabarcoding due to their high 18S rDNA copy number. Altogether, 83 HAB and/or toxic species were identified, among which 10 were dominant, suggesting a high risk of HAB outbreaks in the Beibu Gulf. Phytoplankton abundance increased from south to north and west to east in both seasons, following the high input of dissolved inorganic nitrogen (DIN) and silicate. Excess ammonium input can promote the dominance of Scrippsiella trochoidea and Heterocapsa circularisquama, positioning them as emerging HAB species, while excess DIN caused extreme phosphorus limitation and favored the dominance of Phaeocystis globosa in the Beibu Gulf. This study provided a comprehensive description of the influence of environmental drivers on the phytoplankton community in the Beibu Gulf.

Spatial-temporal distribution of phytoplankton HAB species and contamination status of oyster toxins under intensive oyster farming in Jiangmen coasts, the South China Sea

The diversity and spatial-temporal distribution of phytoplankton HAB species, contamination status of oyster toxin, and their sources were investigated in Jiangmen oyster farming area based on morphological observation and liquid chromatography-tandem mass spectrometry analysis. The results revealed there were 28 HAB species, including 19 harmless HAB species, two toxic species (Akashiwo sanguinea and Karenia brevis), and seven toxin-producing microalgae (Alexandrium pacificum, Dinophysis caudata, D. miles, D. fortii, Gonyaulax spinifera, Gymnodinium catenatum, and K. mikimotoi). The mean abundance of total HAB species generally showed a trend of increasing from winter to autumn. The total average abundances of toxic HAB species were 269, 265, 321 and 2.6 × 103 cells L-1 in winter, spring, summer and autumn, respectively. Redundancy analysis showed temperature, dissolved oxygen, silicate and phosphate were the key factors related with variations of HAB species. Only spring oyster samples were detected paralytic shellfish toxins (PSTs) (1/15), and the composition included gonyautoxins (GTX1&2) and decarbamoyl gonyautoxin 2 (dcGTX2), with a total toxicity level of 9.96 µg STXeq kg-1. N-sulfocarbamyl (C1) and decarbamoyl gonyautoxin 3 (dcGTX3) were observed in the net-concentrated phytoplankton samples. It is inferred that these five types of PSTs derived from A. pacificum. The cultured oysters were observed five types of low-concentration lipophilic marine toxins (LMTs) including okadaic acid (OA), dinophysis-1 (DTX1), pectenotoxin-2 (PTX2), gymnodimine (GYM), and homo-yessotoxin (homo-YTX). Though the levels of PSTs and LMTs in the cultured oysters were low, the presence highlights a potential threat to the safety of oyster products from HAB species.

Spatial-temporal characteristics and the importance of environmental factors in relation to algal blooms in coastal seas

In the prediction of algal blooms with the artificial intelligence (AI) technology, the proper selection of the indices is crucial to the quality of the results. Although eutrophication and climate change have been identified as the main contributors to algal blooms, limited knowledge of the relationships among environmental factors hinders the accurate and correct AI prediction. To explore the importance of environmental factors in the occurrence of algal blooms in coastal seas and facilitate the selection of inputs for the prediction model, this study examines observed water quality data and analyzes their relationship using machine learning methods. It was found that nutrients, water temperature, and salinity play important roles in controlling algal blooms in coastal seas, and their influence varies depending on locations. The high-risk cases of algal bloom mainly occur in environments with medium salinity (22-28 psu), particularly in open estuaries with high nutrient levels. According to the feature importance analysis, water temperature presented the most significant impact on algal bloom risk compared to other environmental factors. The variations of water temperature and Chl-a concentration are synchronous when the water temperature is below 32 °C. Additionally, the relationships between inorganic nutrients (NO3--NO2- and Ortho-P) and Chl-a exhibit asynchronous across multi-temporal scales. The uptake of inorganic nutrients occurs in large quantities within the 24 h prior to the boom in Chl-a concentration. This demonstrates that variations in concentration of inorganic nutrients can be considered as crucial factors in algal bloom predictions, and their significant decrease can be served as an important signal of an impending algal bloom. This study enhances our understanding of the mechanisms behind algal blooms and facilitates their management and prediction strategies in coastal seas.

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.

Molecular diversity and seasonal dynamics of Ostreococcus (Mamiellophyceae, Chlorophyta) in typical mariculture bays based on metabarcoding analysis

Ostreococcus (Mamiellophyceae, Chlorophyta) is a cosmopolitan genus of marine pico-phytoplankton and the smallest free-living photosynthetic eukaryotes with cell size of 1-2 μm. To understand the diversity and spatio-temporal distribution of Ostreococcus in the Rongcheng coastal regions in northern China, metabarcoding analysis based on the 18S rDNA V4 molecular marker was applied to study the molecular diversity and seasonal dynamics of Ostreococcus in three typical mariculture bays (Rongcheng Bay, Lidao Bay and Sanggou Bay). A total of 103 amplicon sequence variants (ASVs) annotated as Ostreococcus were detected in these three typical mariculture bays throughout the year. The top five ASVs in terms of abundance were ASV4, ASV9, ASV14, ASV28 and ASV109, totally occupying 99.1% of Ostreococcus reads. Phylogenetic analysis showed that these five dominant ASVs represented two Ostreococcus ecotypes (OI and OII) and were grouped into four Ostreococcus clades including Ostreococcus lucimarinus (ASV9) and Ostreococcus tauri (ASV28 and ASV109) in OI, and Ostreococcus sp. RC1 (ASV4) and Ostreococcus sp. RC2 (ASV14) in OII, which provided direct evidence to support the co-existence of two ecotypes in the Rongcheng coastal regions. Five dominant ASVs in OI and OII exhibited two distinct seasonal distribution patterns. Three dominant ASVs (ASV9, ASV28 and ASV109) in OI could be detected in all four seasons of the year, exhibiting native distribution properties, while two ASVs (ASV4 and ASV14) in OII decreased sharply in winter and could not be detected in spring, exhibiting characteristics of alien inputs. The composition, succession and association of Ostreococcus community were mainly driven by water temperature in these mariculture bays. This study helps us systematically understand the molecular diversity and distribution patterns of Ostreococcus in typical mariculture bays in northern China, laying the foundation for understanding and revealing the ecological functions of pico-phytoplankton.

Composition and temporal dynamics of the phytoplankton community in Laizhou Bay revealed by microscopic observation and rbcL gene sequencing

Laizhou Bay, a major breeding ground for economic marine organisms in the northern waters of China, is facing rapid environmental degradation. In this study, field surveys in this area were conducted in the spring, summer, and autumn of 2020. Microscopic observation and RuBisCO large subunit (rbcL) gene analysis were employed to understand the community structure and temporal dynamics of phytoplankton. The phytoplankton community structures detected by the two methods showed significant differences. Microscopic observation revealed the dominance of dinoflagellates in spring that shifted to the dominance of diatoms in summer and autumn. However, rbcL gene sequencing consistently identified diatoms as dominant throughout all three seasons, with their relative abundance showing an increasing trend. Conversely, the relative abundance of the second- and third-most abundant taxa, namely, haptophytes and ochrophytes, decreased as the seasons transitioned. rbcL gene sequencing annotated more species than microscopy. It could detect haptophytes and cryptophytes, which were overlooked by microscopy. In addition, rbcL gene sequencing detected a remarkable amount of Thalassiosira profunda, which was previously unidentified in this sea area. However, it appeared to underestimate the contribution of dinoflagellates considerably, with most taxa being only identified through microscopic identification. The two methods jointly identified 28 harmful algal bloom taxa with similar detection quantities but substantial differences in species composition. Phytoplankton communities were influenced by temperature, salinity, and nutrients. The results of this work suggest that a combination of multiple techniques is necessary for a comprehensive understanding of phytoplankton.

Metabarcoding and co-occurrence network reveal significant effects of mariculture on benthic eukaryotic microalgal community: A case study in Daya Bay of the South China Sea

Benthic eukaryotic microalgae were analyzed by metabarcoding the partial 18S rRNA gene in Daya Bay bi-monthly in 2021. Altogether 941 eukaryotic microalgal OTUs were detected, belonging to 27 classes of 8 phyla. Dinophyta and Chlorophyta were the dominant phyla. Microalgal community in the mariculture zone differed significantly from those in non-mariculture zone, reflected by low alpha diversity indexes and increasing abundance and richness of chlorophytes and correspondingly decreasing of dinoflagellates. The abundant occurrences of the pico- and nano-sized taxa such as the chlorophyte Picochlorum in the mariculture zone suggested that nutrient enrichment might result in the miniaturization of the benthic eukaryotic microalgae. The co-occurrence network suggested more negative interactions between taxa in the mariculture zone. A total of 41 algal bloom and/or harmful algal bloom (HAB) species were detected in this study, suggesting a high potential risk of HABs in Daya Bay, especially for the recurrent bloom species Scrippsiella acuminata.

Harmful Microalgae Exhibit Broad Environmental Adaptability in High-Salinity Area Across the Dafengjiang River Estuary

Harmful algal blooms (HABs) often occur in estuaries due to their unique environmental heterogeneity, posing significant environmental and human health risks. However, there is limited understanding of the community composition and community-level change points (thresholds) of harmful microalgae in subtropical estuaries. This study explored harmful microalgae community structure and thresholds in the Dafengjiang River estuary using a metabarcoding approach. The results revealed 63 harmful microalgae species, and major species included Guinardia flaccida, Prorocentrum cordatum, Thalassiosira punctigera, Pseudo-nitzschia galaxiae and T. gravida. Nonparametric change-point analysis and threshold indicator taxa analysis (TITAN) showed threshold responses of harmful microalgae community structure to ammonium (57.5-60 μg·L-1), total phosphorus (27.8-28.5 μg·L-1) and dissolved inorganic phosphorus (14.5-28 μg·L-1) along the salinity gradient. Wider environmental thresholds were also found in hypersaline areas. Additionally, Pyrodinium bahamense, Pfiesteria piscicida, Skeletonema tropicum and T. punctigera were sensitive to environmental changes and thus could be used as bioindicators. Overall, our study unveiled diverse abrupt transitions of harmful microalgal communities, providing a risk assessment for human health and ecological safety in subtropical estuary ecosystems.

Diversity and seasonal occurrence of haptophyta in northern South China Sea through size-fractionated metabarcoding

Haptophyta plays a key role in marine pico-nanoeukaryote communities but information on their diversity and ecology is extremely limited. A total of 103 water samples were collected in northern South China Sea to assess the diversity of haptophyta through metabarcoding targeting 18S V4 rDNA. Furthermore, we investigated the potential genetic differentiation among seasonal occurring Phaeocystis globosa using the high resolution molecular marker pgcp1. 18S V4 rDNA metabarcoding dataset revealed 41 species of haptophytes, with 16 of them as the first record in this region. Notably, six harmful species were detected, including Chrysochromulina leadbeateri, Phaeocystis globosa, and Prymnesium parvum. The pgcp1 marker revealed two clades of Phaeocystis globosa and both of them were present around the year. Clade I was found to predominate in warm season, while Clade III tended to bloom in cold waters. Our results highlight the risk potential of harmful haptophytes in the northern South China Sea.

Distribution, community structure and assembly patterns of phytoplankton in the northern South China Sea

A cruise was conducted in the summer of 2023 from the Pearl River Estuary (PRE) to the adjacent waters of the Xisha Islands in the northern South China Sea (NSCS) to investigate the distribution, community structure, and assembly patterns of eukaryotic and prokaryotic phytoplankton using high-throughput sequencing (HTS) and microscopic observation. Dinophyta were the most abundant phylum in the eukaryotic phytoplankton community based on HTS, accounting for 92.17% of the total amplicon sequence variants (ASVs). Syndiniales was the most abundant order among eukaryotic phytoplankton, whereas Prochlorococcus was the most abundant genus within cyanobacteria. The alpha diversity showed the lowest values in the PRE area and decreased gradually with depth, while cyanobacteria exhibited higher alpha diversity indices in the PRE and at depths ranging from 75 m to 750 m. The morphological results were different from the data based on HTS. Diatoms (37 species) dominated the phytoplankton community, with an average abundance of 3.01 × 104 cells L-1, but only six species of dinoflagellate were observed. Spearman correlation analysis and redundancy analysis (RDA) showed that the distribution and community structure of phytoplankton were largely influenced by geographical location and environmental parameters in the NSCS. The neutral community model (NCM) and null model indicated that deterministic processes played a significant role in the assembly of eukaryotic phytoplankton, with heterogeneous selection and homogeneous selection accounting for 47.27 and 29.95%, respectively. However, stochastic processes (over 60%) dominated the assembly of cyanobacteria and undominated processes accounted for 63.44%. In summary, the formation of eukaryotic phytoplankton was mainly influenced by environmental factors and geographic location, but the assembly of cyanobacteria was shaped by both stochastic processes, which accounted for over 60%, and environmental selection in the NSCS.

Environmental DNA metabarcoding reveals diverse phytoplankton assemblages and potentially harmful algal distribution along the urban coasts of Türkiye

Marine phytoplankton are widely used to monitor the state of the water column due to their rapid changes in response to environmental conditions. In this study, we aimed to investigate the coastal phytoplankton assemblages, including bloom-forming species using high-throughput sequencing of 18S rRNA genes targeting the V4 region and their relationship with environmental variables along the Istanbul coasts of the Sea of Marmara. A total of 118 genera belonging to six phyla were detected. Among them, Dinoflagellata (36) and Bacillariophyta (26) were represented with the highest number of genera. According to the relative abundance of DNA reads, the most abundant taxa were Dinoflagellata_phylum (18.1%), Emiliania (8.4%), Biecheleria (8.4), and Noctiluca (8.1%). The ANOSIM test showed that there was a significant temporal difference in the assemblages, while the driving environmental factors were pH, water temperature, and salinity. According to the TRIX index, the trophic state of the coasts was highly mesotrophic and eutrophic. In addition, 45 bloom-forming and HAB taxa were detected and two species of Noctiluca and Emiliania, which frequently cause blooms in the area, were recorded in high abundance. Our results provide insight into the phytoplankton assemblages along the urbanized coastlines by analysing the V4 region of 18S rRNA. This data can support future studies that use both traditional methods and metabarcoding, employing various primers and targeting different genes and regions.

Metabarcoding reveals a high diversity and complex eukaryotic microalgal community in coastal waters of the northern Beibu Gulf, China

Beibu Gulf is an important semi-enclosed bay located in the northwestern South China Sea, and is famous for its high bio-productivity and rich bio-diversity. The fast development along the Beibu Gulf Economical Rim has brought pressure to the environment, and algal blooms occurred frequently in the gulf. In this study, surface water samples and micro-plankton samples (20-200 μm) were collected in the northern Beibu Gulf coast. Diversity and distribution of eukaryotic planktonic microalgae were analyzed by both metabarcoding and microscopic analyses. Metabarcoding revealed much higher diversity and species richness of microalgae than morphological observation, especially for dinoflagellates. Metabarcoding detected 144 microalgal genera in 8 phyla, while microscopy only detected 40 genera in 2 phyla. The two methods revealed different microalgal community structures. Dinoflagellates dominated in microalgal community based on metabarcoding due to their high copies of 18 s rRNA gene, and diatoms dominated under microscopy. Altogether 48 algal bloom and/or toxic species were detected in this study, 34 species by metabarcoding and 19 species by microscopy. Our result suggested a high potential risk of HABs in the Beibu Gulf. Microalgal community in the surface water samples demonstrated significantly higher OTU/species richness, alpha diversity, and abundance than those in the micro-plankton samples, although more HAB taxa were detected by microscopic observations in the micro-plankton samples. Furthermore, nano-sized taxa, such as those in chlorophytes, haptophytes, and chrysophyceans, occurred more abundantly in the surface water samples. This study provided a comprehensive morphological and molecular description of microalgal community in the northern Beibu Gulf.

A comparison of two gene regions for assessing community composition of eukaryotic marine microalgae from coastal ecosystems

Two gene regions commonly used to characterise the diversity of eukaryotic communities using metabarcoding are the 18S ribosomal DNA V4 and V9 gene regions. We assessed the effectiveness of these two regions for characterising diverisity of coastal eukaryotic microalgae communities (EMCs) from tropical and temperate sites. We binned amplicon sequence variants (ASVs) into the high level taxonomic groups: dinoflagellates, pennate diatoms, radial centric diatoms, polar centric diatoms, chlorophytes, haptophytes and 'other microalgae'. When V4 and V9 generated ASV abundances were compared, the V9 region generated a higher number of raw reads, captured more diversity from all high level taxonomic groups and was more closely aligned with the community composition determined using light microscopy. The V4 region did resolve more ASVs to a deeper taxonomic resolution within the dinoflagellates, but did not effectively resolve other major taxonomic divisions. When characterising these communities via metabarcoding, the use of multiple gene regions is recommended, but the V9 gene region can be used in isolation to provide high-level community biodiversity to reflect relative abundances within groups. This approach reduces the cost of sequencing multiple gene regions whilst still providing important baseline ecosystem function information.

Enantioselective transformation of phytoplankton-derived dihydroxypropanesulfonate by marine bacteria

Chirality, a fundamental property of matter, is often overlooked in the studies of marine organic matter cycles. Dihydroxypropanesulfonate (DHPS), a globally abundant organosulfur compound, serves as an ecologically important currency for nutrient and energy transfer from phytoplankton to bacteria in the ocean. However, the chirality of DHPS in nature and its transformation remain unclear. Here, we developed a novel approach using chiral phosphorus-reagent labeling to separate DHPS enantiomers. Our findings demonstrated that at least one enantiomer of DHPS is present in marine diatoms and coccolithophores, and that both enantiomers are widespread in marine environments. A novel chiral-selective DHPS catabolic pathway was identified in marine Roseobacteraceae strains, where HpsO and HpsP dehydrogenases at the gateway to DHPS catabolism act specifically on R-DHPS and S-DHPS, respectively. R-DHPS is also a substrate for the dehydrogenase HpsN. All three dehydrogenases generate stable hydrogen bonds between the chirality-center hydroxyls of DHPS and highly conserved residues, and HpsP also form coordinate-covalent bonds between the chirality-center hydroxyls and Zn2+, which determines the mechanistic basis of strict stereoselectivity. We further illustrated the role of enzymatic promiscuity in the evolution of DHPS metabolism in Roseobacteraceae and SAR11. This study provides the first evidence of chirality's involvement in phytoplankton-bacteria metabolic currencies, opening a new avenue for understanding the ocean organosulfur cycle.

Diversity and distribution of benthic dinoflagellates in Tonga include the potentially harmful genera Gambierdiscus and Fukuyoa

Benthic dinoflagellates that can cause illness, such as ciguatera poisoning (CP), are prevalent around the Pacific but are poorly described in many locations. This study represents the first ecological assessment of benthic harmful algae species in the Kingdom of Tonga, a country where CP occurs regularly. Surveys were conducted in June 2016 in the Tongatapu island group, and in June 2017 across three island groups: Ha'apai, Vava'u, and Tongatapu. Shallow subtidal coastal habitats were investigated by measuring water quality parameters and conducting quadrat surveys. Microalgae samples were collected using either macrophyte collection or the artificial substrate method. Benthic dinoflagellates (Gambierdiscus and/or Fukuyoa, Ostreopsis, and Prorocentrum) were counted using light microscopy, followed by molecular analyses (real-time PCR in 2016 and high throughput sequencing (metabarcoding) in 2017) to identify Gambierdiscus and Fukuyoa to species level. Six species were detected from the Tongatapu island group in 2016 (G. australes, G. carpenteri, G. honu, G. pacificus, F. paulensis, and F. ruetzleri) using real-time PCR. Using the metabarcoding approach in 2017, a total of eight species (G. australes, G. carpenteri, G. honu, G. pacificus, G. cheloniae, G. lewisii, G. polynesiensis, and F. yasumotoi) were detected. Species were detected in mixed assemblages of up to six species, with G. pacificus and G. carpenteri being the most frequently observed. Ha'apai had the highest diversity with eight species detected, which identifies this area as a Gambierdiscus diversity 'hotspot'. Vava'u and Tongatapu had three and six species found respectively. Gambierdiscus polynesiensis, a described ciguatoxin producer and proposed causative agent of CP was found only in Ha'apai and Vava'u in 2017, but not in Tongatapu in either year. Ostreopsis spp. and Prorocentrum spp. were also frequently observed, with Prorocentrum most abundant at the majority of sites. In 2016, the highest number of Gambierdiscus and/or Fukuyoa cells were observed on seagrass (Halodule uninervis) from Sopu, Tongatapu. In 2017, the highest numbers of Gambierdiscus and/or Fukuyoa from artificial substrate samples were recorded in the Halimeda dominant habitat at Neiafu Tahi, Vava'u, a low energy site. This raised the question of the effect of wave motion or currents on abundance measurements from artificial substrates. Differences in detection were noticed between macrophytes and artificial substrates, with higher numbers of species found on artificial substrates. This study provides a baseline of benthic dinoflagellate distributions and diversity for Tonga that may be used for future studies and the development of monitoring programmes.

Mapping harmful microalgal species by eDNA monitoring: A large-scale survey across the southwestern South China Sea

A large-scale sampling was undertaken during a research cruise across the South China Sea in August 2016, covering an area of about 100,000 km2 to investigate the molecular diversity and distributions of micro-eukaryotic protists, with a focus on the potentially harmful microalgal (HAB) species along the east coast of Peninsular Malaysia. Environmental DNAs from 30 stations were extracted and DNA metabarcoding targeting the V4 and V9 markers in the 18S rDNA was performed. Many protistan molecular units, including previously unreported HAB taxa, were discovered for the first time in the water. Our findings also revealed interesting spatial distribution patterns, with a marked signal of compositional turnover between latitudinal regimes of water masses, where dinophytes and diatom compositions were among the most strongly enhanced at the fronts, leading to distinct niches. Our results further confirmed the widespread distribution of HAB species, such as the toxigenic Alexandrium tamiyavaichii and Pseudo-nitzschia species, and the fish-killing Margalefidinium polykrikoides and Karlodinium veneficum. The molecular information obtained from this study provides an updated HAB species inventory and a toolset that could facilitate existing HAB monitoring schemes in the region to better inform management decisions.

Phytoplankton community structure and environmental factors during the outbreak of Crown-of-Thorns Starfish in Xisha Islands, South China Sea

The "larval starvation hypothesis" proposed that the growing frequency of Crown-of-Thorns Starfish (CoTS) outbreaks could be attributed to increased availability of phytoplankton. However, comprehensive field investigation on the living environment of CoTS larvae and the availability of phytoplankton are still lacking. A cruise was conducted in June 2022 in Xisha Islands, South China Sea, to study the interaction between environmental conditions and phytoplankton communities during CoTS outbreak period. The average concentrations of dissolved inorganic phosphorus (0.05 ± 0.01 μmol L-1), dissolved inorganic nitrogen (0.66 ± 0.8 μmol L-1) and chlorophyll a (0.05 ± 0.05 μg L-1) suggested that phytoplankton may be limited for CoTS larvae in Xisha Islands. Microscopic observation and high-throughput sequencing were used to study the composition and structure of the phytoplankton communities. Bacillariophyta predominated in phytoplankton communities with the highest abundance and species richness. 29 dominant species, including 4 species with size-range preferred by CoTS larvae, were identified in Xisha Islands. The diversity index of all stations indicated a species-rich and structure-stable phytoplankton community in Xisha Islands during the period of CoTS outbreak, which may contribute to CoTS outbreak. These findings revealed the structure of phytoplankton community and environmental factors in the study area during CoTS outbreak, providing the groundwork for future research into the causes and processes of CoTS outbreak.

Phytoplankton community composition in relation to environmental variability in the Urdaibai estuary (SE Bay of Biscay): Microscopy and eDNA metabarcoding

Phytoplankton monitoring is essential for the global understanding of aquatic ecosystems. The present research studies the phytoplankton community of the Urdaibai estuary, combining microscopy and eDNA metabarcoding for the first time in the area. The main aims were to describe the phytoplankton community composition in relation to the environmental conditions of the estuary, and to compare the two methods used. Diatoms Minutocellus polymorphus and Chaetoceros tenuissimus dominated the outer estuary, being replaced by Teleaulax acuta (cryptophyte), Kryptoperidinium foliaceum (dinoflagellate) and Cyclotella spp. (diatom) towards the inner area. This change was mainly prompted by salinity and nutrients. Metabarcoding revealed the presence of 223 species that were not observed by microscopy in previous studies in the estuary. However, several characteristic species (e.g., K. foliaceum) were only detected with microscopy. Additionally, microscopy covered the limitations of eDNA metabarcoding concerning quantification. Thus, to give a full insight, a combination of techniques is recommended.

Annual variation in domoic acid in phytoplankton and shellfish samples from Daya Bay of the South China Sea

Domoic acid (DA) is a well-known phycotoxin that causes amnesic shellfish poisoning (ASP) and is mainly produced by diatom species belonging to the genus Pseudo-nitzschia. An annual survey was conducted monthly over the period of September 2020 to August 2021 in Daya Bay of the South China Sea to investigate the dynamics of particulate and shellfish DA and their relationships with the abundance of Pseudo-nitzschia spp. and environmental parameters. Pseudo-nitzschia spp. was one of the most dominant phytoplankton taxa, and a Pseudo-nitzschia bloom occurred during the survey with the highest abundance of 1.91 × 106 cells L-1. DA was detected in almost all plankton samples with the highest value of 120.7 ng L-1, and high DA concentrations coincided with the abundant presence of Pseudo-nitzschia. DA is prevalent in Daya Bay throughout the year, with detection rates of 98.3%, 82.6%, and 82.6% in plankton samples, in-situ and purchased shellfish, respectively. Higher DA concentrations were detected in the scallop (Chamys nobilis), with the highest concentration of 5.34 µg g-1. High water temperature and low DSi:DIN ratio promoted the growth of Pseudo-nitzschia and DA production. The results suggest that the increasing nitrogen loading and silicate limitation during Pseudo-nitzschia blooms together with the increase in water temperature may increase the risk of DA contamination in Daya Bay.

Metabarcoding analysis identifies high diversity of harmful algal bloom species in the coastal waters of the Beibu Gulf

Harmful algal blooms (HABs) have occurred more frequently in recent years. In this study, to investigate their potential impact in the Beibu Gulf, short-read and long-read metabarcoding analyses were combined for annual marine phytoplankton community and HAB species identification. Short-read metabarcoding showed a high level of phytoplankton biodiversity in this area, with Dinophyceae dominating, especially Gymnodiniales. Multiple small phytoplankton, including Prymnesiophyceae and Prasinophyceae, were also identified, which complements the previous lack of identifying small phytoplankton and those unstable after fixation. Of the top 20 phytoplankton genera identified, 15 were HAB-forming genera, which accounted for 47.3%-71.5% of the relative abundance of phytoplankton. Based on long-read metabarcoding, a total of 147 OTUs (PID > 97%) belonging to phytoplankton were identified at the species level, including 118 species. Among them, 37 species belonged to HAB-forming species, and 98 species were reported for the first time in the Beibu Gulf. Contrasting the two metabarcoding approaches at the class level, they both showed a predominance of Dinophyceae, and both included high abundances of Bacillariophyceae, Prasinophyceae, and Prymnesiophyceae, but the relative contents of the classes varied. Notably, the results of the two metabarcoding approaches were quite different below the genus level. The high abundance and diversity of HAB species were probably due to their special life history and multiple nutritional modes. Annual HAB species variation revealed in this study provided a basis for evaluating their potential impact on aquaculture and even nuclear power plant safety in the Beibu Gulf.

The Effects of Intraguild Predation on Phytoplankton Assemblage Composition and Diversity: A Mesocosm Experiment

Intraguild predation (IGP) can have a significant impact on phytoplankton biomass, but its effects on their diversity and assemblage composition are not well understood. In this study, we constructed an IGP model based on the common three-trophic food chain of "fish (or shrimp)-large branchiopods (Daphnia)-phytoplankton", and investigated the effects of IGP on phytoplankton assemblage composition and diversity in outdoor mesocosms using environmental DNA high-throughput sequencing. Our results indicated that the alpha diversities (number of amplicon sequence variants and Faith's phylogenetic diversity) of phytoplankton and the relative abundance of Chlorophyceae increased with the addition of Pelteobagrus fulvidraco, while similar trends were found in alpha diversities but with a decrease in the relative abundance of Chlorophyceae in the Exopalaemon modestus treatment. When both predators were added to the community, the strength of collective cascading effects on phytoplankton alpha diversities and assemblage composition were weaker than the sum of the individual predator effects. Network analysis further showed that this IGP effect also decreased the strength of collective cascading effects in reducing the complexity and stability of the phytoplankton assemblages. These findings contribute to a better understanding of the mechanisms underlying the impacts of IGP on lake biodiversity, and provide further knowledge relevant to lake management and conservation.

Long-Term Toxicity of 50-nm and 1-μm Surface-Charged Polystyrene Microbeads in the Brine Shrimp Artemia parthenogenetica and Role of Food Availability

Micro and nanoplastics (MNPs) as emerging contaminants have become a global environmental issue due to their small size and high bioavailability. However, very little information is available regarding their impact on zooplankton, especially when food availability is a limiting factor. Therefore, the present study aims at evaluating the long-term effects of two different sizes (50 nm and 1 μm) of amnio-modified polystyrene (PS-NH₂) particles on brine shrimp, Artemia parthenogenetica, by providing different levels of food (microalgae) supply. Larvae were exposed to three environmentally relevant concentrations (5.5, 55, and 550 μg/L) of MNPs over a 14-days of exposure with two food levels, high (3 × 10⁵~1 × 10⁷ cells/mL), and low (1 × 10⁵ cells/mL) food conditions. When exposed to high food levels, the survival, growth, and development of A. parthenogenetica were not negatively affected at the studied exposure concentrations. By comparison, when exposed to a low food level, a U shape trend was observed for the three measured effects (survival rate, body length, and instar). Significant interactions between food level and exposure concentration were found for all three measured effects (three-way ANOVA, p < 0.05). The activities of additives extracted from 50 nm PS-NH₂ suspensions were below toxic levels, while those from 1-μm PS-NH₂ showed an impact on artemia growth and development. Our results demonstrate the long-term risks posed by MNPs when zooplankton have low levels of food intake.

Ostreopsis Schmidt and Coolia Meunier (Dinophyceae, Gonyaulacales) from Cook Islands and Niue (South Pacific Ocean), including description of Ostreopsis tairoto sp. nov

It is important to decipher the diversity and distribution of benthic dinoflagellates, as there are many morphologically indistinct taxa that differ from one another in production of potent toxins. To date, the genus Ostreopsis comprises twelve described species, of which seven are potentially toxic and produce compounds presenting a threat to human and environmental health. In this study, isolates previously identified as "Ostreopsis sp. 3" were sampled from the area where it was first reported, Rarotonga, Cook Islands, and have been taxonomically and phylogenetically characterised as Ostreopsis tairoto sp. nov. Phylogenetically, the species is closely related to "Ostreopsis sp. 8", O. mascarenensis, "O. sp. 4", O. fattorussoi, O. rhodesiae and O. cf. siamensis. Previously, it was considered a part of the O. cf. ovata complex but can be distinguished from O. cf. ovata based on the small pores identified on this study, and from O. fattorussoi and O. rhodesiae based on relative lengths of the 2' plates. No known palytoxin -like compounds were detected in strains investigated in this study. Strains of O. lenticularis, Coolia malayensis and C. tropicalis were also identified and described. This study advances our knowledge of biogeography, distribution, and toxins of Ostreopsis and Coolia species.

The biogeography and co-occurrence network patterns of bacteria and microeukaryotes in the estuarine and coastal waters

Community and diversity shifts of bacteria and microeukaryotes with strong environmental and spatial variations have been unveiled in the Pearl River Estuary (PRE) and northern coastal part of South China Sea (SCS). However, it is not clear what the determining factors shape the microbial community and how the biotic interactions respond to the estuarine and oceanic environment. Here, we established the multiple regression models (MRM) and co-occurrence networks on microbial communities in PRE and SCS habitats. The results showed that there were significant differences of the abiotic factors affecting the bacterial and microeukaryotic communities between PRE and SCS habitats. Salinity explained the largest variations to the microbial community dissimilarities in PRE. Whereas spatial and environmental factors determined the microbial community dissimilarities in SCS. Positive relations between parasitic lineages (e.g. Perkinsea and Cercozoa) and algal taxa (Dinophyceae, Cryptophyta, Chlorophyta and Ochrophyta) dominated in the PRE network. While parasites Syndiniales positively correlated with other Syndiniales and protists in SCS. Strong positive associations among autotrophic and heterotrophic groups were revealed in both niches. Therefore, the biotic interactions are also important and may be responsible for the unexplained variations of the abiotic factors from MRM models. Microbial network in the PRE estuarine water had weakened resistance to environmental disturbances, while the SCS network had greater capacity to maintain network stability. This study shed light on the different mechanisms of abiotic and biotic factors in shaping the compositions of bacteria and microeukaryotes between PRE and SCS niches, and highlights the weakening effect of environmental disturbances on the microbial network stability.

Comparison of photosynthetic responses between haptophyte Phaeocystis globosa and diatom Skeletonema costatum under phosphorus limitation

The diatom Skeletonema costatum and the haptophyte Phaeocystis globosa often form blooms in the coastal waters of the South China Sea. Skeletonema costatum commonly dominates in nutrient enrichment coastal waters, whereas P. globosa starts flourishing after the diatom blooms when phosphorus (P) is limited. Therefore, P limitation was proposed to be a critical factor affecting diatom-haptophyte transition. To elucidate the tolerance to P limitation in P. globosa compared with S. costatum, the effect of P limitation on their photosystem II (PSII) performance was investigated and their photosynthesis acclimation strategies in response to P limitation were evaluated. P limitation did not affect the growth of P. globosa over 7 days but decreased it for S. costatum. Correspondingly, the PSII activity of S. costatum was significantly inhibited by P limitation. The decline in PSII activity in S. costatum under P limitation was associated with the impairment of the oxygen-evolving complex (the donor side of PSII), the hindrance of electron transport from Q_A ^- to Q_B (the acceptor side of PSII), and the inhibition of electron transport to photosystem I (PSI). The 100% decrease in D1 protein level of S. costatum after P limitation for 6 days and PsbO protein level after 2 days of P limitation were attributed to its enhanced photoinhibition. In contrast, P. globosa maintained its photosynthetic activity with minor impairment of the function of PSII. With accelerated PSII repair and highly increased non-photochemical quenching (NPQ), P. globosa can avoid serious PSII damage under P limitation. On the contrary, S. costatum decreased its D1 restoration under P limitation, and the maximum NPQ value in S. costatum was only one-sixth of that in P. globosa. The present work provides extensive evidence that a close interaction exists between the tolerance to P limitation and photosynthetic responses of S. costatum and P. globosa.

Karenia selliformis bloom dynamics and growth rate estimation in the Sfax harbour (Tunisia), by using automated flow cytometry equipped with image in flow, during autumn 2019

A Karenia selliformis bloom event in the Gulf of Gabès (Mediterranean Sea), was monitored over 9 days at high frequency during fall 2019, by using an automated flow cytometer (Cytosense, Cytobuoy b.v.) with an image-in-flow attachment. The instrument recorded the shape of the optical signals that lead to the resolution of six cell groups of pico-, nano- and microphytoplankton, during the Harmful Algal Bloom (HAB). K. selliformis cell dimensions derived from the hourly records, enabled to estimate the daily division rate over the bloom period. Results revealed that K. selliformis was the only bloom-forming species and it reached its highest mean abundance the fourth day of the survey. A shift in the nutrient composition occurred with a potential P limitation during the bloom growth and N limitation during the bloom collapse. The co-inertia analysis revealed opposite patterns for K. selliformis and heterotrophic prokaryotes suggesting trophic interactions and possible mixotrophic behaviour of K. selliformis at the end of the bloom. K. selliformis exhibited low growth rates generally < 1 division day^-1, which could not explain the observed high abundance. The tide played a crucial role in the dynamics of K. selliformis at a semi-diurnal scale and at spring-neap tide scale and was probably enhancing K. selliformis accumulation.