Journal of Oceanology and Limnology   2022, Vol. 40 issue(6): 2120-2145     PDF       
http://dx.doi.org/10.1007/s00343-022-1322-z
Institute of Oceanology, Chinese Academy of Sciences
0

Article Information

XIE Hang, ZOU Jian, ZHENG Chengzhi, QU Yuchen, HUANG Kaixuan, LÜ Songhui
Biodiversity and distribution of benthic dinoflagellates in tropical Zhongsha Islands, South China Sea
Journal of Oceanology and Limnology, 40(6): 2120-2145
http://dx.doi.org/10.1007/s00343-022-1322-z

Article History

Received Oct. 1, 2021
accepted in principle Dec. 8, 2021
accepted for publication Feb. 14, 2022
Biodiversity and distribution of benthic dinoflagellates in tropical Zhongsha Islands, South China Sea
Hang XIE1,3#, Jian ZOU1,3#, Chengzhi ZHENG1,3, Yuchen QU1,3, Kaixuan HUANG1,3, Songhui LÜ1,2,3     
1 Research Center of Harmful Algae and Marine Biology, College of Life Science and Technology, Jinan University, Guangzhou 510362, China;
2 Southern Marine Science and Engineering Guangdong Laboratory, Zhuhai 519000, China;
3 Key Laboratory of Eutrophication and Red Tide Prevention of Guangdong Higher Education Institutes, Jinan University, Guangzhou 510632, China
Abstract: Benthic dinoflagellates have attracted increasing attention in recent years because of their toxicity and ability to form extensive harmful algal blooms. Ostreopsis producing palytoxin and its analogs, Gambierdiscus producing ciguatera toxins, and Prorocentrum producing okadaic acid and dinophysis toxins, have been concerned as serious human poisonings. We explored the benthic dinoflagellate biodiversity and distribution characteristics of a series of tropical reefs in 20–40-m water depth in wet season in the Zhongsha Islands in South China Sea using morphological, phylogenetic, and cell counting methods. Results show that benthic dinoflagellates in the islands are rich in biodiversity and 15 species from genera Amphidinium, Coolia, Ostreopsis, and Prorocentrum were identified: Amphidinium carterae, A. magnum, A. massartii, A. operculatum, Coolia canariensis, C. malayensis, C. palmyrensis, C. tropicalis, Ostreopsis cf. ovata, Prorocentrum concavum, P. cf. sculptile, P. emarginatum, P. hoffmannianum, P. lima, and P. rhathymum. Among them, A. magnum is reported for the first time in Chinese waters. The abundance of benthic dinoflagellates was relatively low at 88–4 345 cells/100 cm2 on sediment and 10–91 cells/g on macroalgae. Prorocentrum and Amphidinium were the dominant and subdominant genera, respectively. It is speculated that the low abundance of benthic dinoflagellates is closely related to the scarcity of macroalgae and stronger water motion at the depth >15 m in Zhongsha Islands. This study expanded the study in biodiversity of benthic dinoflagellates in Chinese waters, and revealed the distribution characteristics of harmful benthic microalgae in reef habitats.
Keywords: benthic dinoflagellates    Zhongsha Islands    South China Sea    biodiversity    distribution    
1 INTRODUCTION

Benthic dinoflagellates, e.g., species of Ostreopsis, Gambierdiscus, Fukuyoa, Prorocentrum, Coolia, and Amphidinium, are group of microalgae that attach to substrates such as macroalgae, dead coral, sand, and rocks (Fukuyo, 1981; Boisnoir et al., 2019a). Harmful benthic algal blooms have aroused widespread concern due to benthic-dinoflagellate-produced toxins and bloom-induced anoxia, which endangers human health and marine ecosystem (Yasumoto et al., 1987; Litaker et al., 2010; Mangialajo et al., 2011; Kudela et al., 2017). Palytoxin and its analogs, the metabolites of some Ostreopsis species, are the most toxic non-protein algal toxins (Amzil et al., 2012; Ben Gharbia et al., 2016). The ciguatoxins and maitotoxins generated by species of Gambierdiscus and Fukuyoa, which are responsible for ciguatera fish poisoning (CFP) in humans, are the most serious non-bacterial illness associated with seafood consumption (Litaker et al., 2010). Some species of Prorocentrum produce okadaic acid and its analogs, which are the culprit of diarrheic shellfish poisoning (DSP) in humans (Tripuraneni et al., 1997; Foden et al., 2005). Furthermore, some Coolia and Amphidinium species can threaten marine benthic animals (Pagliara and Caroppo, 2012; Karafas et al., 2017).

In 1979, with the discovery of Gambierdiscus toxicus — a toxic species isolated from a ciguatera-endemic area — the study of harmful benthic dinoflagellates was initiated (Yasumoto et al., 1977; Adachi and Fukuyo, 1979). An increase in number of CFP incidences and Ostreopsis blooms reported from around the world has led to the increasing studies highlighting their ecological impacts (Litaker et al., 2010; Rhodes, 2011). In recent decades, scholars have paid more attention to the taxonomy, toxicity, and ecology of benthic dinoflagellates, as well as the geographical distribution of toxic and harmful species (Kudela et al., 2017). CFP producers, species of Gambierdiscus and Fukuyoa, are mainly reported in tropical and subtropical waters (Bienfang et al., 2008; Chinain et al., 2020). A review summarized that the Caribbean Sea and the Western Pacific Ocean are the areas where these two genera are most requently reported (Chinain et al., 2021). This distribution pattern has also expanded to higher latitude temperate waters, including Korea (Jang et al., 2018) and the northeast Mediterranean Sea in the past decade (Aligizaki and Nikolaidis, 2008). Ostreopsis species are distribution globally. Blooms of Ostreopsis have frequently recorded in waters of Mediterranean Sea (Mangialajo et al., 2011; Rhodes, 2011; Kang et al., 2013; Acaf et al., 2020; Nascimento et al., 2020). Besides Ostreopsis, other genera Coolia, Amphidinium, Prorocentrum, and Sinophysis is with higher biodiversity of benthic dinoflagellates (Rhodes, 2011; Ben Gharbia et al., 2016, 2019; Reñé et al., 2020). Benthic dinoflagellates are abundantly distributed on coral reefs and island ecosystems. Fukuyo (1981) isolated eleven benthic dinoflagellate species belonging to five genera from coral reefs and systematically described the whole benthic dinoflagellate community. Some islands in the South Pacific, such as French Polynesia, the Cook Islands, and the Republic of Kiribati, are considered to be "hot spots" of Gambierdiscus biodiversity (Chinain et al., 2016). So far studies of benthic dinoflagellates have been mostly focused on a certain genus or toxic species. There have been relatively few studies on entire benthic dinoflagellate communities (Nishimura et al., 2013; Boisnoir et al., 2018; Lim and Jeong, 2021). Those studies of community composition and substrate preferences have only been carried out in a few areas (Yong et al., 2018; Mustapa et al., 2019; Accoroni et al., 2020).

To date, three species of Gambierdiscus (Zhang et al., 2016), one species of Fukuyoa (Leung et al., 2018), two species of Ostreopsis (Zhang et al., 2018), four species of Coolia (Leung et al., 2017), eleven species of Amphidinium (Luo et al., 2022), and more than ten Prorocentrum species, in which some of the species were okadaic acid and/or DSP toxins producers (Luo et al., 2017; Zou et al., 2021), have been reported in the South China Sea. Recently, a harmful Prorocentrum concavum bloom was reported at Hainan Island, which is the first benthic dinoflagellate bloom in China at present (Zou et al., 2020).

To date, most studies on benthic dinoflagellates have been concentrated in intertidal areas, bays, lagoons, and islands, especially in the Caribbean Sea and Mediterranean Sea (Mangialajo et al., 2011; Durán-Riveroll et al., 2019; Chinain et al., 2021). A study of how water movement, depth, and different habitats influenced on the abundance and distribution of benthic dinoflagellate communities was carried out at Johnston Atoll, Pacific Ocean (Richlen and Lobel, 2011). It was found that moderate hydrodynamic disturbance promoted the distribution of benthic dinoflagellate community. The sampling depth of this study was concentrated below 5 m and the maximum depth was only 13 m. Lee et al. (2020) studied the effects of substratum and depth on benthic harmful dinoflagellate assemblages at Perhentian Islands and Terengganu, Malaysia. Though the deepest sampling depth was 25 m, sampling sites were mostly on the coasts of the islands. In short, there are few studies on reefs of more than 20-m depth, especially the areas far away from the mainland. Compared with the coasts of continents or large islands, reef ecosystems with water depths of more than 20 m are subject to stronger hydrodynamic disturbances. The Zhongsha Islands are the typical tropical coral reef ecosystems composed of reef-building corals and the marine organisms that grow in and around the reefs. The main body of the Zhongsha Islands, Zhongsha Great Atoll, is fully submerged beneath water depth of 13–26 m. The seabed surface of Zhongsha Great Atoll is coral sand and shell debris. In this study, we investigated the biodiversity and distribution of benthic dinoflagellates in the Zhongsha Islands based on morphology, phylogeny, and cell counting. Samples of sediments and macroalgae were collected to analyze the epibenthic dinoflagellate community composition and abundance at different shoals.

2 MATERIAL AND METHOD 2.1 Study area and sampling strategy

The Zhongsha Islands (13°57′N–19°33′N, 113°02′E–118°45′E), are composed of 26 named shoals on Zhongsha Great Atoll and Huangyan Island, and four scattered shoals: the Yitong, Xianfa, Shenhu, and Zhongnan Shoals. The Zhongsha Islands belong to a tropical monsoon climate zone, with the nature characteristics of strong sunshine exposure, low temperature variation, and high temperature and salinity. In addition, the Zhongsha Islands have distinct dry and wet seasons, and abundant precipitation that increases from north to south.

Sediments and macroalgal samples were collected by scuba diving between June 24 and July 5, 2020, at the Zhongsha Great Atoll (six stations), Yitong Shoal, and Shenhu Shoal, China (Fig. 1 & Table 1). Sediment samples were collected in triplicate using a 10-cm×10-cm×3-cm mud sampler and were transferred with seawater into a 1-L sample bag (Fig. 2). Sediment samples were placed for 10 min for sediment setting, and then seawater was slowly collected and the volume was recorded. Later, the seawater was filtered with a 120-μm sieve to collect the epibenthiccells. Benthic dinoflagellates were retained by a sieve with 20-μm mesh, then immediately fixed in 1.5% Lugo's solution. Water samples were collected using a water sampler. Macroalgae were sampled carefully with the surrounding seawater and transferred into sampling bags. Macroalgal samples were shaken for 30 s to detach attached dinoflagellate cells and then rinsed twice with 0.22-μm filtered seawater. The seawater containing epibenthic algal cells were filtered using a 120-μm sieve and then epibenthic algal cells were retained by a 20-μm mesh, and immediately fixed in 1.5% Lugol's solution. The separated macroalgae was dried using absorbent paper before weighing (Boisnoir et al., 2019b).

Fig.1 The location of South China Sea (a) and sampling sites in Zhongsha Islands (b), Yitong Shoal, and Shenhu Shoal (c) Map review No. GS(2016)1569.
Fig.2 The method of sediment sampling
Table 1 Locations of sampling sites, sampling depth, and sample types
2.2 Dinoflagellate cell abundance

Cells were counted by 1-mL Sedgewick Rafter Counting Cell using the light microscope (Olympus BX 61, Tokyo, Japan). The cell abundance of the sediment was described as cells per 100 cm2. Abundance on macroalgae was measured as cells per gram of fresh macroalgal weight.

2.3 Isolation and culture

Live cells during the exponential period were isolated with micropipette under the light microscope. Isolated benthic dinoflagellates cells were cultured in a 96-well plate with L1 medium (Guillard and Hargraves, 1993). The culture plate was under a dark꞉light cycle of 12 h꞉12 h with 100-μmol/(m2·s) irradiance, 25 ℃, and salinity of 30. Strains that divided a number of cells were ported to a 25-mL culture flask with a same culture conditions.

2.4 Morphological observation

Cell during the exponential period was observed under a microscope. Morphological data (cell shape, length, and width) were obtained using Image-Pro Plus 6.0 Image.

The morphology and distribution of chloroplasts and nuclei were observed using a fluorescence microscope, then image acquisition was conducted using the QImaging Retiga 400R digital camera (QImaging, Surrey, BC, Canada) at 400 magnifications.

Sample preparation for the scanning electron microscope (SEM) was performed according to Lim et al. (2019). Cells grown in the exponential phase were fixed in glutaraldehyde (final concentration: 2.5%) for 12 h at 4 ℃. One milliliter of algal cells was pipetted onto a round glass slide and soaked in 10% polylysine for more than 2 h. The slide was observed under a microscope to ensure that a large number of cells were adsorbed on the surface. The cells were soaked using a gradient of 100%, 80%, 60%, 40%, and 20% of the seawater and Milli-Q water for 5 min. Dehydration was carried out through a series of ethanol gradients (10%, 30%, 50%, 70%, and 90% and three times at 100%, for 5 min each). The dehydrated cells were dried by a CO2 critical point dryer (Leica Microsystems, Mannheim, Germany) and putter-coated with gold (Leica, Germany). Finally, the morphology and structure of cell was presented and photographed under a scanning electron microscope (Carl Zeiss Inc., Germany).

2.5 DNA isolation and gene amplification

The clonal cells of all benthic dinoflagellate strains were harvested at 12 000×g for 1 min at 4 ℃. Total genomic DNA of exponential cells were extracted by MiniBEST Universal Genomic

DNA Extraction Kit (TaKaRa, Japan) according to the instructions. A 50-μL PCR reaction system includes: 2-μL DNA, 25-μL Accurate TaqMaster Mix (2X) (Zoman Biotechnology, China), each 1 μL of the forward primer and reverse primer, and 21-μL ultrapure water. The D1–D3 regions of large subunit (LSU) was amplified using general primers (D1R: 5ʹ-ACCCGCTGAATTTAAGCATA-3ʹ, D3B: 5ʹ-TCGGAGGGAACCAGCTACTA-3ʹ) described in Scholin et al. (1994). The PCR cycle procedure: pre-denaturation at 94 ℃, 4 min and followed by 36 cycles, each cycle included denaturation at 94 ℃ for 20 s, annealing at 56 ℃ for 30 s, and extension at 72 ℃ for 1.5 min, and a final extension step at 72 ℃ for 10 min (Zou et al., 2020). All sequences were uploaded to GenBank and the accession numbers were listed in Supplementary Table S4.

2.6 Phylogenetic analysis

34 sequences of the LSU (D1–D3) regions obtained by sequencing was aligned and edited by ContigExpress software. The obtained sequences from our study and related species sequences from GenBank, which consist of 41 Amphidinium sequences, 43 Coolia sequences, 61 Ostreopsis sequences, and 64 Prorocentrum sequences, were analyzed for multiple alignments at the website (https://www.ebi.ac.uk/Tools/msa/muscle/). Sequence group of different genera were aligned separately. Then the aligned results were break off both ends using BioEdit (Version 7.2.9) and then were analyzed using the maximum likelihood (ML) analysis in an online program (https://www.phylo.org/) with RAxML-HPC2 and XSEDE V. 8.2.10. The best models were determined by the program MrModeltest2.3 under the Akaike Information Criterion. In addition, Bayesian inference (BI) was conducted at Mr. Bayes 3.2.7 with four Markov Chain Monte Carlo (MCMC) chains and each chain was run for 1 000 000 cycles and sampling every 100 cycles using the best fit model (GTR). The uncorrected P-distances among all species and its closely related species were used by the software Mega X and the results are displayed in the Supplementary Table S5.

3 RESULT 3.1 Biodiversity of benthic dinoflagellates species

Fifteen species of benthic dinoflagellates belonging to four genera were identified in morphology: Amphidinium carterae, A. magnum, A. massartii, A. operculatum, Coolia canariensis, C. malayensis, C. palmyrensis, C. tropicalis, Ostreopsis cf. ovata, Prorocentrum concavum, P. cf. sculptile, P. emarginatum, P. hocmffcmmannianum, P. lima, and P. rhathymum. We focus the morphology of a species first discovered in the waters of China (A. magnum), and potentially widespread toxic species (O. cf. ovata, P. lima, and P. hocmffcmmannianum). The morphological micrographs of P. cf. sculptile and A. carterae were not obtained because of cultivation failure.

3.1.1 Morphology of benthic dinoflagellates 3.1.1.1 Amphidinium magnum (Fig. 3)
Fig.3 Light microscopy (LM) and scanning electron microscopy (SEM) photographs of Amphidinium magnum a. LM, showing the shape of cell; b–c. fluorescence LM, showing the chloroplasts (b) and the position of nucleus (c); d–e. SEM, showing the micrographs of dorsal view (d) and ventral view (e). Scale bars: a–c:  10 μm; d–e: 2 μm.

Cells are 26.94–38.10 μm long (mean 34.05±2.31 μm, n=30) and 20.04–29.70 μm wide (mean 24.04±1.93 μm, n=30). The length/width (L/W) ratio is from 1.03 to 1.67 (mean 1.42±0.10, n=30). The cells are generally ovoid, and the epicone extends to the right with an approximate triangular structure. The antapex of the hypocone is rounder, while its right is almost vertical (Fig. 3a). No pyrenoid is observed but numerous chloroplasts fill the whole cell with a block distribution (Fig. 3b). The nucleus is posterior (Fig. 3c). The longitudinal flagellum is located approximately a third of the way up the cell (Fig. 3e). Irregular scales are scattered across the hypocone surface. The ventral ridge is long and narrow, which is about half the length of the cell, connecting the two flagella insertion points (Fig. 3e).

3.1.1.2 Amphidinium massartii and Amphidinium operculatum (Fig. 4; Supplementary Fig.S1 and Table S1)
Fig.4 Light microscopy (LM) photographs of Amphidinium massartii (a–c) and Amphidinium operculatum (d–f) a–b. LM, showing the shape of cell; b, c, e, and f. fluorescence LM, showing the chloroplasts (b and e) and the position of nucleus (c and f). Scale bars: 10 μm.

Amphidinium massartii (Fig. 4ac; Supplementary Fig.S1a–b): Cells are 16.44–21.92 μm long (mean 19.52±1.57 μm, n=30) and 10.31–18.08 μm wide (mean 15.34±1.72 μm, n=30). The ratio of L/W is from 1.21 to 1.59 (mean 1.29±0.20, n=30). Cell shape varies greatly from round to elliptical and the epicone extends to the left with a crescent structure (Fig. 4a). No pyrenoid is observed. The nucleus is oval and locates at posterior (Fig. 4c). The antapex of hypocone is ovoid generally while its right is almost vertical (Fig. 4a; Supplementary Fig.S1a–b).

Amphidinium operculatum (Fig. 4df; Supplementary Fig.S1c–d): Cells are 27.57–37.95 μm long (mean 31.89±2.66, n=30) and 20.04–29.70 μm wide (mean 18.73–26.39 μm, n=30). The ratio of L/W is from 1.24 to 1.63 (mean 1.43±0.09, n=30). Cell shape varies greatly from round to elliptical and the epicone extends to the right with an approximate triangular structure and the hypocone is relatively round and broad (Fig. 4d; Supplementary Fig.S1c–d). Numerous chloroplasts are fill the whole cell with a block distribution (Fig. 4e). The nucleus is oval and located at posterior (Fig. 4f).

3.1.1.3 Coolia canariensis, Coolia malayensis, Coolia palmyrensis, and Coolia tropicalis (Fig. 5; Supplementary Fig.S2 and Table S2)
Fig.5 Light microscopy (LM) photographs of Coolia canariensis, Coolia malayensis, Coolia palmyrensis, and Coolia tropicalis a, d, g, and j. LM, showing the shape of cells; b, c, e, f, h, i, k, and l. fluorescence LM, showing the chloroplasts (b, e, h, and k) and the position of nucleus (c, f, i, and l). Scale bars: 10 μm.

Coolia canariensis (Fig. 5ac; Supplementary Fig. S2a): Cells are round or oval in ventral view (Fig. 5a; Supplementary Fig.S2a). They are 25.28–36.34 μm long (dorsal to ventral, mean 30.93±2.46 μm, n=30), 23.53–33.53 μm wide (transdiameter, mean 25.51±2.93 μm, n=30). The ratio of L/W is from 1.02 to 1.24 (mean 1.11±0.06; n=30). The chloroplast is granular and fills the entire cell (Fig. 5b). The elongated nucleus locates at the posterior (Fig. 5c).

Coolia malayensis (Fig. 5df; Supplementary Fig. S2b): Cells are round or oval in ventral view (Fig. 5d; Supplementary Fig.S2b). They are 20.21–33.79 μm long (dorsal to ventral, mean 25.56±2.87 μm, n=30), 17.47–26.90 μm wide (transdiameter, mean 22.38±2.36 μm, n=30). The ratio of L/W is from 1.01 to 1.30 (mean 1.14±0.08; n=30). The chloroplast is granular and fills the entire cell (Fig. 5e). A bean-shaped nucleus locates at the right posterior end of cell. (Fig. 5f).

Coolia palmyrensis (Fig. 5gi; Supplementary Fig.S2c): Shape of cells is nearly spherical (Fig. 5g; Supplementary Fig.S2c). They are 21.01–25.97 μm long (dorsal to ventral, mean 23.89±1.41 μm, n=30), 17.43–24.04 μm wide (transdiameter, mean 21.46±1.67 μm, n=30). The ratio of L/W is from 1.01 to 1.22 (mean 1.12±0.07; n=30). The granular chloroplast was filled the entire cell (Fig. 5h). A bean-shaped nucleus locates at the right posterior end of cell (Fig. 5i).

Coolia tropicalis (Fig. 5jl; Supplementary Fig. S2d): Shape of cells is nearly spherical (Fig. 5j; Supplementary Fig.S2d). They are 28.74–40.84 μm long (dorsal to ventral, mean 34.82±3.08 μm, n=30), 25.34–35.71 μm wide (transdiameter, mean 32.07±3.33 μm, n=30). The ratio of L/W is from 1.02 to 1.19 (mean 1.09±0.05; n=30). The granular chloroplast are fill the entire cell (Fig. 5k). A U-shaped nucleus locates at the middle posterior end of cell (Fig. 5l).

3.1.1.4 Ostreopsis cf. ovata(Fig. 6)
Fig.6 Light microscopy (LM) and scanning electron microscopy (SEM) photographs of Ostreopsis cf. ovata a. LM, showing the shape of cell; b–c. fluorescence LM, showing the chloroplasts (b) and the position of nucleus (c); d–f. SEM, showing the micrographs of apical view (d), antapical view (e), and the apical pore plate (f). Po: pore plate. Scale bars: a–c:  10 μm; d–e: 2 μm; f: 1 μm.

The cells of this species are compressed and tear-shaped in general (Fig. 6a). In the ventral view, they are 32.55–42.87 μm long (dorsal to ventral, mean 36.04±2.72 μm, n=30), 22.30–32.95 μm wide (transdiameter, mean 25.51±2.93 μm, n=30). The L/W ratio ranges from 1.24 to 1.60 (1.42±0.09; n=30). An oval nucleus located at the posterior end of the cell (Fig. 6c). Thecal plate pattern: Po, 3ʹ, 7ʹʹ, pyrenoid (arrowhead), 5ʹʹʹ, 2ʹʹʹʹ. The plate surface is smooth with scattered pores (Fig. 6df). The pore plate (Po) is narrow and long in shape with a few pores and generally looks like a slit (Fig. 6f). In the epitheca, it is hard to observe apical plate 2ʹ, which is very close to Po in the ventral view. Compared with apical plate 2ʹ, the other two apical plates are much larger. Apical plate 1ʹ is a broad hexagon. The first apical plate (1ʹ) is much longer and broader than the other two apical plates (Fig. 6d). Additionally, apical plate 1ʹ occupies the middle in the ventral view. The hypotheca is made up of eight plates. Antapical plates 2ʹʹʹʹ are hexagonal and smaller than postcingular plates 2ʹʹʹ, 3ʹʹʹ, and 4ʹʹʹ (Fig. 6e). All postcingular plates are about the same size, except postcingular plate 1ʹʹʹ. The ventral pore is not observed.

3.1.1.5 Prorocentrum hoffmannianum (Fig. 7)
Fig.7 Light microscopy (LM) and scanning electron microscopy (SEM) photographs of Prorocentrum hoffmannianum a. LM, showing the pyrenoid (arrowhead) and shape of cell; b, c. fluorescence LM, showing the chloroplasts (b) and the position of nucleus (c); d–f. SEM, showing the micrographs of left thecal plate (d) and right thecal plate (e) and V-shaped periflagellar area (f). Scale bars: a–e:  10 μm; f: 1 μm.

Cells of P. hocmffcmmannianum are oval and symmetrical (Fig. 7a). They are 46.81–53.91 μm long (mean 50.52±1.60 μm, n=30) and 40.05–51.60 μm wide (mean 47.32±2.66 μm, n=30). The L/W ratio varies from 0.98 to 1.27 (mean 1.13±0.05, n=30). The pyrenoid is locates at the cell center, with abundant chloroplasts scattered all around the cell (Fig. 7b). The nucleus is located at the posterior (Fig. 7c). Kidney-shaped pores are scattered across the thecal surface, with no pore in the central area (Fig. 7de). The number of thecal pores varies from 139 to 159 and 0.52–0.96 μm length, 0.23–0.42 μm. Numerous small and shallow depressions are observed in the thecal plate and the intercalary band is surrounded by marginal pores varies from 82 to 93 (Fig. 7d). Additionally, there is a big depression in the plate center. Small round to ovoid pores were found within some depressions. The left thecal margin exhibited a flared and flattened curved apical collar that bordered the periflagellar area (Fig. 7e). An apical collar is just above the flagellar region. The periflagellar area is wide and V-shaped, with a large flagellar pore (fp) and an accessory pore (ap) (Fig. 7f). The plate patten is 1, 2, 3, 4, 5, 6, 7, and 8.

3.1.1.6 Prorocentrum lima (Fig. 8)
Fig.8 Light microscopy (LM) and scanning electron microscopy (SEM) photographs of Prorocentrum lima a. LM, showing the pyrenoid (arrowhead) and shape of cell; b–c. fluorescence LM, showing the chloroplasts (b) and the position of nucleus (c); d–f. SEM, showing the micrographs of left thecal plate (d) and right thecal plate (e) and V-shaped periflagellar area (f). Scale bars: a–c:  10 μm; d–e: 2 μm; f: 1 μm.

Cells of P. lima are oval, slightly thin, and symmetrical (Fig. 8a). They are 34.39–41.17 μm long (mean 38.27±1.68 μm, n=50) and 21.36–29.65 μm wide (mean 30.29±1.69 μm, n=50) with a L/W ratio varies from 1.34 to 1.67 (mean 1.45±0.07, n=50). The pyrenoid is observed in cell center with a starch ring (Fig. 8a). Numerous chloroplasts fill the whole cell (Fig. 8b). The thecal plate is smooth, and round or kidney-shaped pores are dispersed across the surface (Fig. 8df). Approximately 56–71 pores are observed, but no pore is observed in the center of the plate (Fig. 8d). The nucleus is oval and locates at the posterior (Fig. 8c). The intercalary bands are surrounded by a series of 44–56 neatly arranged pores (Fig. 8de). The periflagellar area is wide V-shaped. Eight platelets and two large pores are observed (Fig. 8f).

3.1.1.7 Prorocentrum concavum, Prorocentrum emarginatum, and Prorocentrum rhathymum (Fig. 9; Supplementary Fig.S3 and Table S3)
Fig.9 Light microscopy (LM) photographs of Prorocentrum concavum, Prorocentrum emarginatum, and Prorocentrum rhathymum a, d, and g. LM, showing the shape of cell; b, c, e, f, h, and i. fluorescence LM, showing the chloroplasts (b, e, and h) and the position of nucleus (c, f, and i). Scale bars: 10 μm.

Prorocentrum concavum (Fig. 9ac; Supplementary Fig.S3a): cells of Prorocentrum concavum are oval and symmetric (Fig. 9a; Supplementary Fig.S3a). They are 43.77–49.29 μm long (mean 46.7±1.67 μm, n=30) and 35.80–44.70 μm wide (mean 41.52±2.17 μm, n=30). The ratio of L/W is from 1.04 to 1.27 (mean 1.13±0.05, n=30). A pyrenoid could be brightly observed in the center of cell, which be rounded by a starch ring (Fig. 9a). Numerous chloroplasts fill the whole cell (Fig. 9b). Nucleus locates at the posterior of the cell (Fig. 9c).

Prorocentrum emarginatum (Fig. 9df; Supplementary Fig.S3b): cells of Prorocentrum emarginatum are broad, oval and asymmetric (Fig. 9d; Supplementary Fig.S3b). They are 28.56–40.56 μm long (mean 35.36±2.63 μm, n=30) and 20.67– 33.77 μm wide (mean 27.63±3.22 μm, n=30). The ratio of L/W is from 1.07 to 1.73 (mean 1.29±0.16, n=30). No pyrenoid is observed, while numerous chloroplasts fill the whole cell (Fig. 9e). The nucleus is oval and locates at the posterior of the cell (Fig. 9f).

Prorocentrum rhathymum (Fig. 9gi; Supplementary Fig.S3c): cells of this species are oval and asymmetric (Fig. 9g; Supplementary Fig.S3c). They are 27.79–36.86 μm long (mean 31.73±2.23 μm, n=30) and 19.22–26.43 μm wide (mean 22.55±1.98 μm, n=30). The ratio of L/W is from 1.25 to1.73 (mean 1.41±0.11, n=30). A pyrenoid could be observed in the center of cell, which be rounded by a starch ring (Fig. 9g). Numerous chloroplasts fill the whole cell (Fig. 9h). The nucleus is oval and located at the posterior of the cell (Fig. 9i).

3.1.2 Phylogenetic analysis 3.1.2.1 Amphidinium

In the present study seven strains, comprising four species, were used for phylogenetic analysis of the Amphidinium genus with related species sequences (Fig. 10). Three A. massartii strains isolated from this study were grouped together and they formed a sister clade to the Brazilian strain Am_Cub_1. The A. carterae strain ZS130 differed from the Chinese strain 3WZD8, Australian strain CS740, and New Zealand strain CAWD57. Amphidinium magnum was clustered with three other Bahamian strains with maximal support. Two A. operculatum strains were used for phylogenetic analysis. Strain ZS701 differed from strain WZD225, all of which were isolated from the South China Sea. Strain ZS716 was clustered with the Brazilian strain Ao-Rec-1 and Australian strain K_0663.

Fig.10 Phylogenetic tree of partial LSU rDNA (D1–D3) sequences of Amphidinium based on the maximum likelihood (ML) and Bayesian inference (BI) Heterocapsa sp. (CCMP424) is used as outgroup. Only bootstrap values > 50% and posterior probabilities > 0.8 are shown. Sequences obtained in the present study are in bold.
3.1.2.2 Coolia

In the present study, a total of eight strains, consisting of four species, were used for phylogenetic analysis of the Coolia genus with related species sequences (Fig. 11). The C. malayensis strains ZS719 and ZS733 were grouped with two Hong Kong strains SKLMP_W059 and SKLMP_W072, and the USA strain CCMP1345 with a 53 Bootstrap value. In addition, C. palmyensis strain ZS601 clustered together with the northeast Atlantic strain Cp1412_1 and two Hong Kong strains, SKLMP_S017 and SKLMP_W085. Three C. canariensis strains acquired in the present study clustered together and shared identical sequences with the Hainan Island strain D1C2. Two C. tropical strains, ZS602 and ZS619, grouped together with strain XS554 isolated from the Xisha Islands, South China Sea.

Fig.11 Phylogenetic tree of partial LSU rDNA (D1–D3) sequences of Coolia based on the maximum likelihood (ML) and Bayesian inference (BI) Ostreopsis ovata (OvPR04) is used as outgroup. Only bootstrap values > 50% and posterior probabilities > 0.8 are shown. Sequences obtained in the present study are in bold.
3.1.2.3 Ostreopsis

Five Ostreopsis sequences belonging to this study were aligned with related species sequences from GenBank. The phylogenetic relationship is shown in Fig. 12. The O. cf. ovata clade contained three groups, one consisted of some strains isolated from the western Pacific (South China Sea and Gulf of Thailand) and Indonesia. A Brazil and two Portugal (East Atlantic) strains were formed the Atlantic group. The remaining strains isolated from the Mediterranean Sea (Spain, Greece, Italy, Tunisia, and Manoca) and West Atlantic area (Brazil) clustered together as the Mediterranean/West Atlantic group. All O. cf. ovata strains acquired in the present study grouped together with other strains from the western Pacific.

Fig.12 Phylogenetic tree of partial LSU rDNA (D1–D3) sequences of Ostreopsis based on the maximum likelihood (ML) and Bayesian inference (BI) Coolia canariensis (CMJJ1) is used as outgroup. Only bootstrap values > 50% and posterior probabilities > 0.8 are shown. Sequences obtained in the present study are in bold.
3.1.2.4 Prorocentrum

In this study, 14 strains belonging to four Prorocentrum species, were used for phylogenetic analysis of the Prorocentrum genus with related species sequences (Fig. 13). There was one point for which the D1–D3 region sequence of P. concavum was not obtained. Based on the phylogenetic tree, the P. lima complex separated into four subclades, which roughly showed different geographical distributions. Nine strains belonging to P. lima were isolated from this study. Prorocentrum lima strain ZS102 was grouped together with other Chinese strains, and P. lima strain ZS516 was grouped together with strain TIO115a and strain S4, isolated from Guangxi, China, and Australia, respectively. Another seven P. lima strains were clustered with strain NMN07, and they were a sister clade of a Martinique strain. Two P. hocmffcmmannianum strains isolated from this study were the sister clade of the La Réunion strain PBMA_01. Prorocentrum rhathymum strain ZS405 clustered with the Xisha Islands strain XS501. The P. cf. sculptile strain ZS615 formed a clade and clustered with two Martinique strains. P. cf. emarginatum strain ZS251 was the sister clade of La Réunion strain PERN06. Besides, strain ZS251 grouped with three other strains, two Xisha Islands strains, and one was isolated from the USA.

Fig.13 Phylogenetic tree of partial LSU rDNA (D1–D3) sequences of Prorocentrum based on the maximum likelihood (ML) and Bayesian inference (BI) Adenoides eludens (ADE15) is used as outgroup. Only bootstrap values > 50% and posterior probabilities > 0.8 are shown. Sequences obtained in the present study are in bold.
3.2 Distribution of benthic dinoflagellates 3.2.1 Spatial distribution of benthic dinoflagellates

At Oliver Shoal (15 m), Hand Shoal (30 m), and Shenhu Shoal (17 m), no macroalgal sample was collected, so all species at those sites were isolated from coral sand (Table 1). Hand Shoal was located at the middle of Zhongsha Great Atoll, which had a deepest sampled depth of 30 m and only 3 species were isolated from coral sand at this Shoal. Additionally, 8 species were only identified from Padina australis at Magpie Shoal (20 m) (Table 2). At south Zhongsha Great Atoll, 10 benthic dinoflagellates species were identified from all types of samples at Addington Patch, which had a medium sampled depth of 20 m. Also, 7 species were identified from coral sand or Padina australis at Cawston Shoal with a middle sampled depth of 23 m (Table 2). Amphidinium carterae and A. magnum were only isolated from macroalgae Padina australis at Magpie Shoal and coral sand at Oliver Shoal, respectively, and A. massartii appeared at half of the Shoals, whatever sediment or all macroalgae (Table 2). Except Oliver Shoal and Hand Shoal, A. operculatum was found at all shoals from coral sand, Padina australis, and Halimeda sp. In Shenhu Shoal, C. canariensis and C. malayensis were present. C. palmyrensis and C. tropicalis only isolated from coral sand and Padina australis at Addington Patch (Table 2). Additionally, except Oliver Shoal, O. cf. ovata appeared at all Shoal. Prorocentrum cf. sculptile appeared at half of the shoals, including Oliver Shoal, east of Addington Patch, Magpie Shoal, and Yitong Shoal. This species was found in all types of samples (sediment or macroalgae). In addition, P. concavum was only identified from coral sand at Hand Shoal with a 30-m sampled depth. At all shoals, Prorocentrum lima was found in coral sand, Amphiroa fragiroa, Padina australis, and Halimeda sp. (Table 2). However, P. hocmffcmmannianum only appeared in coral sand and Padina australis at Oliver Shoal (Table 2). Only identified from coral sand, Padina australis, and Halimeda sp., P. rhathymum were found at Magpie Shoal, Cawston Shoal, and Yitong Shoal (Table 2). Furthermore, although classification data for Gambierdiscus were not obtained in this study, Gambierdiscus spp. appeared at all shoals, no matter coral sand or other macroalgae except Oliver and Cawston Shoal (Table 2). In short, Benthic dinoflagellates are most diversely represented on south of Zhongsha Great Atoll and Yitong Shoal, as well as macroalgae. Additionally, more species were identified from macroalgae samples Padina australis.

Table 2 Distribution of benthic dinoflagellate species in different shoals
3.2.2 Quantitative distribution of benthic dinoflagellates

The abundance of benthic dinoflagellates in the sediment was between 88 and 4 345 cells/100 cm2 (Fig. 14). Except for Shenhu Shoal and west of Addington Patch, the genus Prorocentrum was dominant at all sites. The lowest cell number averaged 88 cells/100 cm2 was at Oliver Shoal, which has the lowest sampled depth of 15 m, while the maximum cell abundance averaged 88 cells/100 cm2 was found at Cawston Shoal with a 23-m sampled depth (Fig. 14). The two sites were located at the northernmost and southernmost part of the Zhongsha Great Atoll (Fig. 14a). Only Prorocentrum and Ostreopsis species were found at Oliver Shoal, with an average of 82 cells/100 cm2 and 5 cells/100 cm2, respectively. At Hand Shoal with the deepest sampled depth, a cell abundance of averaged 334 cells/100 cm2 was collected and Coolia species were absent in samples collected (Fig. 14a). At southern Zhongsha Great Atoll, Amphidinium was dominant, with cell abundance averaging 54 cells/100 cm2 at west of Addington Patch. However, at east of Addington Patch, cell abundance was as high as 1 815 cells/100 cm2, which was more than 10 times that west of Addington Patch (142 cells/100 cm2), although their sampling depths are very close (20 m and 19 m, respectively). Yitong Shoal (15 m) and Shenhu Shoal (17 m) had a lower sampled depth and there were similar levels of cell abundance of benthic dinoflagellates were recorded, respectively averaging 642 cells/100 cm2 and 534 cells/100 cm2 (Fig. 14b). Compared with that at other shoals, the maximum cell abundance of Coolia was found at Shenhu Shoal, with an average of 408 cells/100 cm2, representing approximately 76.4% of all benthic dinoflagellates (Fig. 14b). Gambierdiscus species were not observed in samples from Cawston Shoal, or Shenhu Shoal (Fig. 14ab). As the dominant genus, Prorocentrum occupied more than 50% of the total abundance of benthic dinoflagellates at all shoals, except west of Addington Patch, where they comprised only 6.8%. Different distribution patterns in the abundance of benthic dinoflagellates were observed on the east and west sides of Addington Patch. The benthic dinoflagellates abundance near the outside of Zhongsha Great Atoll was higher, while that of the inside was much lower. At the same time, the abundance in the south of Zhongsha Great Atoll was higher than that in the northern (Fig. 14a). In addition, the relationship between depth and abundance was not very close, but the abundance at a medium depth of 20 m was much higher (averaged 1 815 cells/100 cm2 at East of Addington Patch and averaged 4 345 cells/100 cm2 at Cawston Shoal).

Fig.14 The cell abundances (cells/100 cm2) of five benthic dinoflagellates on sediment in Zhongsha Great Atoll (a), Yitong Shoal, and Shenhu Shoal (b)

Since only three kinds of macroalgae (Padina australis, Halimeda sp., and Amphiroa fragiroa) were collected from all sites, macroalgal preference is not discussed in this study. The abundance of benthic dinoflagellates on macroalgae was between 10 and 91 cells/g and Prorocentrum was dominant on macroalgae at all sites (Fig. 15). The northernmost site of Zhongsha Great Atoll, Magpie Shoal, had a total abundance of benthic dinoflagellates averaged 57 cells/g on Padina australis with a sampled depth of 20 m. However, the abundance of Gambierdiscus and Ostreopsis was as low as 1 cell/g (Fig. 15a). At west of Addington Patch, which had a sampled depth of 19 m, benthic dinoflagellates cell abundance averaged 49 cells/g on Padina australis, the difference was that Coolia comprised 33.1% of total abundance. At east of Addington Patch, which had the same depth as Magpie Shoal, the cell abundance of benthic dinoflagellates was lowest on Amphiroa fragiroa, averaging 10 cells/g, which was significantly different from that at other sampling sites (Fig. 15a). Additionally, no Gambierdiscus spp. or Coolia spp. were observed on Amphiroa fragiroa. At Yitong Shoal, with the lowest sampled depth of 15 m, the highest abundance of benthic dinoflagellates was found on Halimeda sp., with an average of 91 cells/g (Fig. 15a). In addition, all five benthic dinoflagellate genera were recorded from macroalgal samples. The benthic dinoflagellates abundance of Yitong Shoal was higher than that of Zhongsha Great Atoll because of its low water depth (Fig. 15ab). At north Zhongsha Great Atoll, Addington Patch had the lowest abundance. On macroalgae samples, the highest abundance was observed at Halimeda sp. (91 cells/g at Yitong Shoal), while the lowest abundance appeared at Amphiroa fragiroa (10 cells/g at east of Addington Patch).

Fig.15 Cell abundances (cells/g) of five benthic dinoflagellates on macroalgae in Zhongsha Great Atoll (a), Yitong Shoal, and Shenhu Shoal (b)
4 DISCUSSION 4.1 Biodiversity of benthic dinoflagellates in Zhongsha Islands 4.1.1 Amphidinum

Amphidinium magnum identified in Zhongsha Islands is the first record in Chinese waters. It is also the first report in addition to type species locality, Grand Bahama Beach, Bahamas (Karafas et al., 2017). Its morphological characteristics were consistent with those originally described. A long and narrow ventral ridge was clearly observable, but the pusule was absent in our strain. Scales were visible on the hypocone surface, which were absent in the strain from the type locality.

The phylogenetic tree based on the LSU region showed that the species isolated from Zhongsha Islands can be divided into two clades, one clade consisting of A. massartii and A. carterae, both of which are smaller than A. operculatum, which formed a clade on its own, with 100% bootstrap support in the phylogeny. Amphidinium magnum grouped together with the three type locality strains, with maximum support. Additionally, the position of A. magnum was close to that of A. cf. thermaeum and A. theodori, which is consistent with the results of Karafas et al. (2017).

4.1.2 Coolia

Morphological characteristics (size, plate arrangement, and shape) of all four Coolia species, including C. canariensis, C. malayensis, C. palmyrensis, and C. tropicalis, identified from Zhongsha Islands, were generally similar to that of the type species descriptions. The four species were reported previously from Hong Kong and Hainan Island in China (Leung et al., 2017; Zhang et al., 2021).

The phylogenetic tree of Coolia showed that the C. malayensis closely related to the Hong Kong strains and New Zealand strain, which was consistent with that conducted by Zhang et al. (2021) in Hainan Island. The results reveal that all C. malayensis strains isolated from Chinese waters have a common ancestor. Strains of C. palmyensis from Zhongsha Islands and from NE Atlantic were a sister clade of Hong Kong strain, which suggested that there were few differences between Zhongsha Islands strain and Hong Kong strain in phylogeny (Leung et al., 2017). Two strains of C. tropical isolated from Zhongsha Islands grouped together with Xisha Islands, which also showed that there are differences among strains even they are geographically adjacent.

4.1.3 Ostreopsis

Morphologically, it is difficult to distinguish the various species of Ostreopsis, which are similar to the type species O. siamensis. The cell size of O. cf. ovata strain from the Zhongsha Islands was 32.55–42.87 μm long (dorsal to ventral, mean 36.04±2.72 μm), 22.30–32.95 μm wide (transdiameter, mean 25.51±2.93 μm), which is significantly smaller than that of O. siamensis (dorso-ventral 60–100 μm; transdiameter 45–90 μm). In addition, the morphological characteristics of the Zhongsha Islands strain (ZS718) in this article were generally consistent with those of strains isolated by Fukuyo (1981), which were generally tear-shaped with a L/W ratio of 1.5–2.0 (Fig. 6). In original description, a large antapical plate of O. cf. ovata was designated as postcingular intercalary plate (1p) (Fukuyo, 1981), now it is antapical plate (2ʹʹʹʹ). Compared with the original description, our plate pattern (Po, 3ʹ, 7ʹʹ, ?C, ?S, 5ʹʹʹ, 2ʹʹʹʹ) has been reported in recent literatures (Chomérat et al., 2020; Nascimento et al., 2020).

Penna et al. (2010) revealed that O. cf. ovata could phylogenetically divide into three clades, including Atlantic/Mediterranean/Pacific clade, Indian/ Pacific clade, and Atlantic/Indian/Pacific clade. All O. cf. ovata strains from this study were clustered together to form an Indian/Pacific clade. The Indian/ Pacific clade included all strains collected from the South China Sea (including the newly discovered Zhongsha Islands strain, Vietnamese strain, and Malaysian strain) and adjacent waters (Thai strain and Indonesian strain). The consistency in genetic characteristics and geographic distribution may be because of the similarity in ecosystems and the influence of ocean circulation. The islands and land surrounding the South China Sea make it a marginal sea, and the unique genetic characteristics of O. cf. ovata suggest that it is genetically separated from the Mediterranean/Indian Ocean O. cf. ovata strain. Thus, geographical isolation has had a strong influence on phylogenetic divergence.

4.1.4 Prorocentrum

A number of studies have confirmed that P. lima shows genetic diversity and phenotypic plasticity (Hoppenrath et al., 2013; Zhang et al., 2015). Similar to O. cf. ovata, P. lima is also a cosmopolitan distributed species. Based on morphology and phylogeny constructed with LSU and internal transcribed spacer (ITS) regions, of strains of P. lima from Hainan Island, China were classified as five morphotypes (Zhang et al., 2015). Then Nascimento et al. (2017) defined morphotype 4 of P. lima as a new species, P. caipirignum. A study also discovered two morphological types (morphotypes 1 and 2) of P. lima in the Caribbean (Chomérat et al., 2019). The results showed that these two types had obvious morphological differences. Morphotype 2 was oblong oval, the cell length and width were 37.1–38.2 μm and 27.8–32.3 μm, respectively (L/W ratio was 1.18– 1.36). The morphology and cell measurements of a P. lima strain (ZS102) from Zhongsha Islands was in accordance with the morphotype 2 in the Caribbean. The cell of P. hocmffmannianum strain from Zhongsha Islands was broader than that of the type species, but it was similar to that of P. belizeanum (Faust, 1990). Herrera-Sepúlveda et al. (2015) summarized that P. hoffmannianum and P. belizeanum have some micromorphological similarities (cell size, surface ornamentation, and periflagellar area architecture) and molecular data lack sufficient evidence to distinguish these two species. Therefore, they are considered to be conspecifics.

On the base of phylogenetic tree of LSU rDNA (D1–D3 region), we also identified our strains (ZS102 and ZS516) from Zhongsha Islands belonged to P. lima morphotype 2. Furthermore, the P. lima complex could be clearly distinguished from the similar species P. caipirignum and had high genetic diversity, based on phylogenetic analysis of D1–D3 region. This same result is observed in previous studies (Zhang et al., 2015; Chomérat et al., 2019). Normally, the same geographical origin have similar phylogenies. However, the P. lima complex in this study is divided into three groups. The first group has a neighboring geographical origin, the second is distributed in the northern and southern hemispheres, and the third strains are geographically separated from each other. This result suggests that P. lima in the same area also has different genetic characteristics.

Hoppenrath et al. (2013) summarized that P. belizeanum, P. hocmffcmmannianum, and P. maculosum were highly similar in morphology with small genetic differences. But their toxin-producing characteristics were different, which blurred the boundaries between these species. However, with the application of more precise molecular techniques, these species are all attributed to the P. hocmffcmmannianum species complex. In this study, the Zhongsha Islands strains of P. hocmffcmmannianum were a sister clade of a La Réunion strain, which was once defined as P. belizeanum. Therefore, we further demonstrate that our P. hocmffcmmannianum strains are closer to P. belizeanum, which is consistent with the morphological results. However, these were our subclades of P. hocmffcmmannianum and these were directly related to geographic locations, indicating that geographical isolation has had a considerable influence on phylogenetic divergence. Therefore, another marker is needed to identify P. hocmffcmmannianum species more accurately.

4.2 Distribution of benthic dinoflagellates from Zhongsha Islands 4.2.1 Geographical distribution of benthic dinoflagellates

There are approximately 200 species of Amphidinium, both freshwater and seawater (Guiry and Guiry, 2022). Many species were heterotrophic (Taylor, 1971; Calado et al., 1998). Most of Amphidinium species identified in Zhongsha Islands were widespread and commonly reported in tropical and subtropical waters, as well as in temperate higher latitude regions (Murray and Patterson, 2002; Lee et al., 2013; Jiang, 2019; Luo et al., 2022).

Leaw et al. (2016) provided a new scenario for the phylogeographic history of Coolia and it was speculated that the divergence of the two species, C. malayensis and C. monotis, may be due to geographic segmentation caused by vicariant events. Four Coolia species, C. canariensis, C. malayensis, C. palmyrensis, and C. tropicalis identified from the Zhongsha Islands, were widely distributed in the western Pacific (C. tropicalis and C. canariensis in Vietnam, C. canariensis and C. malayensis in the Korean Sea, C. malayensis, C. palmyrensis, and C. tropicalis in Australia) (Jeong et al., 2012; Ho and Nguyen, 2014; Larsson et al., 2019). Aside from the western Pacific, C. malayensis was also found in Florida, Belize, and the Cook Islands (Leaw et al., 2016). Coolia palmyrensis has been found in the Caribbean Sea (Karafas et al., 2015) and western Atlantic Ocean (De Azevedo Tibiriçá et al., 2020). Coolia tropicalis was first isolated from an Atlantic barrier reef ecosystem (Faust, 1995) and has then been expanded to Brazil (Nascimento et al., 2019) and many areas of the western Pacific (Larsson et al., 2019; Zhang et al., 2021).

The globally distributed species O. cf. ovata has been reported in the Caribbean Sea, Mediterranean Sea, Atlantic, Indian, and Pacific Oceans (Monti et al., 2007; Mangialajo et al., 2008; Rhodes, 2011; Kang et al., 2013; Zhang et al., 2018). Parsons et al. (2012) summarized the global distribution of Ostreopsis and found that Ostreospsis species were also widely present in the western and southern Pacific, as well as western and northern Atlantic, except for the Mediterranean and Caribbean regions where the species was widespread. The genetic analyses indicated that O. cf. ovata were divided into three geographical areas: Atlantic/Mediterranean/ Pacific, Atlantic/Indian/Pacific, and Pacific (Penna et al., 2010; Mangialajo et al., 2011; Sato et al., 2011). Ostreopsis species were commonly observed in Caribbean and Western Atlantic region such as Florida Keys (Accoroni et al., 2020), Brazil (Nascimento et al., 2020), and Bezile (Faust, 1995). Zhang et al. (2018) reported O. cf. ovata and O. lenticularis from Hainan Island, South China Sea. Zheng et al. (2017) found that the O. cf. ovata isolated in Beibu Bay, southern China, were contain higher concentrations of ovatoxins-c (OVTX-c) and OVTX-d/e, as well as small amounts of OVTX-g and OVTX-f. So far there were only two Ostreopsis species have been reported in tropical waters of the southern China.

Prorocentrum as the dominant genus is a consistent feature of tropical waters (Mohammad-Noor et al., 2005; Parsons and Preskitt, 2007). Prorocentrum lima is a widely distributed Prorocentrum species, which is often observed in water column, on macroalgal surfaces and sediments (Fukuyo, 1981; Faust, 1990; Aligizaki et al., 2008). Prorocentrum hocmffcmmannianum was first isolated in Belize (Faust, 1990). Subsequently, P. hocmffcmmannianum, P. belizeanum, and P. maculosum were revised as a P. hocmffcmmannianum species complex (Herrera-Sepúlveda et al., 2015; Rodríguez et al., 2018). P. hocmffcmmannianum species complex were mainly tropical and subtropical distribution, such as the Mediterranean (mez, 2003), Australia, New Zealand (Rhodes and Smith, 2019), and North America (Torres-Ariño et al., 2019). Prorocentrum concavum was widely distributed in the water column and reef ecosystems from tropical to high latitude temperate zones (Faust, 1994; Luo et al., 2017).

4.2.2 Quantitative distribution of benthic dinoflagellates

Although the Zhongsha Islands are called an archipelago, its main body, Zhongsha Great Atoll, is actually a series of reefs located below sea surface, far from human activities in water depths of approximately 20–40 m. The habitat is mainly coral reef ecosystems, which are greatly affected by ocean water dynamics such as turbulence. In this study, the environment in which benthic dinoflagellates grow is characterized by weak light intensity, low temperature relative to the intertidal zone, and low nutrient concentrations.

Benthic dinoflagellates are generally epiphyte on the seafloor matrix (sand, silt, or debris shell) or on the surface of macrophytes, and are less abundant in water column unless bloom occurs (Fukuyo, 1981; Boisnoir et al., 2019a). Prorocentrum was dominant on coral sand and three macroalgaes (Padinaaustralis, Halimeda sp., and Amphiroa fragiroa) except for those to the west of Addington Patch and Shenhu Shoal. A similar phenomenon that Prorocentrum was the first dominant genus has also observed in previous researches (Richlen and Lobel, 2011; Widiarti and Anggrain, 2012; Accoroni et al., 2020). Interestingly, our results indicate that Ostreopsis was the least abundant group in the Zhongsha Islands. The maximum density of Ostreopsis was appeared at Yitong Shoal with averaged 64 cells/100 cm2 on coral sand. Lee et al. (2020) studied the effect of substratum and depth on benthic dinoflagellates assemblages and found that Ostreopsis were the second most dominant genera of benthic dinoflagellate assemblages, both in inorganic matrix and macrophytes. The maximum abundance of Ostreopsis were observed in hard coral and turf algae. A similar phenomenon was observed at Johnston Atoll in the Pacific, which has a similar geographical characteristics and environmental conditions with Zhongsha Islands (Richlen and Lobel, 2011). Ostreopsis spp. from Johnston Atoll accounted for the largest proportion at several of hard coral reef sites, an opposite result to our study. Boisnoir et al. (2019a) indicated that Ostreopsis preferred to attach to blade-ramified macrophytes, particularly Florideophyceae. However, the minimum abundance was presented on Padina australis with averaged 8 cells/g. We speculate that other factors may have contributed to this phenomenon, such as the sampling depth. The maximum sampling depth at Johnston Atoll was 13 m, which was lower than the least depth in our study. Although Ostreopsis had generally low abundance on the three macroalgae and coral sand, they were significantly more abundant on blades with larger surface areas, e.g. Halimeda sp. (8 cells/g). Whereas the abundance on branched Amphiroa fragiroa and foliated Padina australis were as low as 1 cell/g.

Like Prorocentrum, Amphidinium was widely distributed in various substrates with a high abundance. Amphidinium was the genus with the second highest abundance in the benthic dinoflagellate community of Zhongsha Islands, both on coral sands and macroalgae. This was inconsistent with those of other studies that the abundance of Amphidinium was generally the lowest (Richlen and Lobel, 2011; Lee et al., 2020). Coolia is present in low abundance. The same results are consistent with previous studies in the Gulf of Mexico (Okolodkov et al., 2007; Hachani et al., 2018; Ben Gharbia et al., 2019).

A study on the depth distribution of benthic dinoflagellates in the Caribbean demonstrated that their distribution was significantly correlated with light intensity. The cell density decreased significantly with increasing depth (Boisnoir et al., 2018). The benthic dinoflagellate density was < 100 cells/g in both wet and dry seasons at a depth of 20 m. It was found that no Ostreopsis cells were observed when the water depth exceeds 10 m. Interestingly, in this study, Ostreopsis spp. were observed at all depths on coral sands and macroalgae except that on Padina australis at Cawston Shoal. A different trend was observed by Lee et al. (2020) that at depths greater than 12 m, only Prorocentrum and Coolia were presented. In our study, on coral sand samples, Gambierdiscus were only observed at 15 m (Yitong Shoal) and 19 m (West of Addington Patch). Although Gambierdiscus was also present in the samples of the macroalgae Padina australis at Cawston Shoal (20 m) and Magpie Shoal (23 m), Their abundance was as low as 1 cell/g. On macroalgae samples, the maximum abundance of Prorocentrum were observed at the minimum depth of 15 m (Yitong Shoal). At a maximum depth of 30 m, a second minimum Prorocentrum abundance of 334 cells/100 cm2 on coral sand were presented. However, no Coolia spp. was observed. Macroalgae generally grow in abundance in the intertidal zone with water depth less than 3 m. The Zhongsha Islands have a water depth of more than 15 m, and the bottom is a coral reef ecosystem. Increased depth means less macroalgae. Compared to the coasts of continents or large islands, reef ecosystems with water depths of more than 20 m are subject to stronger hydrodynamic disturbances. Benthic dinoflagellates were significantly more abundant in calm lagoons than on reefs heavily affected by waves (Boisnoir et al., 2018). Hydrodynamic conditions may be the most important factor affecting the growth of benthic-epiphytic dinoflagellates. The abundance in sheltered areas is significantly higher than that in unsheltered areas (Shears and Ross, 2009; Mangialajo et al., 2011). It is speculated that the low abundance of benthic dinoflagellates is closely related to the scarcity of macroalgae and stronger water motion at the depth > 15 m in Zhongsha Islands.

5 CONCLUSION

The biodiversity and distribution of benthic dinoflagellates in the Zhongsha Islands have been revealed. In general, 15 species of benthic dinoflagellates belonging to four genera were identified: A. carterae, A. magnum, A. massartii, A. operculatum, C. canariensis, C. malayensis, C. palmyrensis, C. tropicalis, O. cf. ovata, P. cf. sculptile, P. concavum, P. emarginatum, P. hocmffcmmannianum, P. lima, and P. rhathymum. Benthic dinoflagellates were distributed broadly on coral sands and three species of macroalgae (Padina australis, Halimeda sp., and Amphiroa fragiroa). However, the total abundance of benthic dinoflagellates in the Zhongsha Islands, which was 88–4 345 cells/100 cm2 on sediment and 10– 91 cells/g on macroalgae, is much lower than that in other areas. Prorocentrum and Amphidinium were the dominant and subdominant genera, respectively. It is speculated that the low abundance of benthic dinoflagellates is closely related to the scarcity of macroalgae and stronger water motion at the depth > 15 m in Zhongsha Islands. This study has provided a deeper understanding of benthic dinoflagellate assemblages in the Zhongsha Islands.

6 DATA AVAILABILITY STATEMENT

The datasets generated and/or analyzed during the current study are available from the corresponding author on reasonable request.

Electronic supplementary material

Supplementary material (Supplementary Tables S1–S5 and Figs.S1–S6) is available in the online version of this article at https://doi.org/10.1007/s00343-022-1322-z .

References
Acaf L, Saab M A A, Khoury-Hanna M, et al. 2020. Bloom dynamics of the newly described toxic benthic dinoflagellate Ostreopsis fattorussoi along the Lebanese coast (Eastern Mediterranean). Regional Studies in Marine Science, 38: 101338. DOI:10.1016/j.rsma.2020.101338
Accoroni S, Totti C, Romagnoli T, et al. 2020. Distribution and potential toxicity of benthic harmful dinoflagellates in waters of Florida Bay and the Florida Keys. Marine Environmental Research, 155: 104891. DOI:10.1016/j.marenvres.2020.104891
Adachi R, Fukuyo Y. 1979. The thecal structure of a marine toxic dinoflagellate Gambierdiscus toxicus gen. et sp. nov. collected in a ciguatera-endemic area.. Nippon Suisan Gakkaishi, 45(1): 67-71. DOI:10.2331/SUISAN.45.67
Aligizaki K, Nikolaidis G F. 2008. Morphological identification of two tropical dinoflagellates of the genera Gambierdiscus and Sinophysis in the Mediterranean Sea. Journal of Biological Research-Thessaloniki, 9: 75-82.
Amzil Z, Sibat M, Chomerat N, et al. 2012. Ovatoxin-a and palytoxin accumulation in seafood in relation to Ostreopsis cf.ovata blooms on the French Mediterranean coast. Marine Drugs, 10(2): 477-496. DOI:10.3390/md10020477
Ben Gharbia H, Laabir M, Ben Mhamed A, et al. 2019. Occurrence of epibenthic dinoflagellates in relation to biotic substrates and to environmental factors in Southern Mediterranean (Bizerte Bay and Lagoon, Tunisia): an emphasis on the harmful Ostreopsis spp, Prorocentrum lima and Coolia monotis.. Harmful Algae, 90: 101704. DOI:10.1016/j.hal.2019.101704
Ben-Gharbia H, Yahia O K, Amzil Z, et al. 2016. Toxicity and growth assessments of three thermophilic benthic dinoflagellates (Ostreopsis cf. ovata, Prorocentrum lima and Coolia monotis) developing in the southernMediterranean basin.. Toxins, 8(10): 297. DOI:10.3390/toxins8100297
Bienfang P, Oben B, DeFelice S, et al. 2008. Ciguatera: the detection of neurotoxins in carnivorous reef fish from the coast of Cameroon, West Africa. African Journal of Marine Science, 30(3): 533-540. DOI:10.2989/ajms.2008.30.3.8.642
Boisnoir A, Pascal P Y, Cordonnier S, et al. 2018. Depth distribution of benthic dinoflagellates in the Caribbean Sea. Journal of Sea Research, 135: 74-83. DOI:10.1016/j.seares.2018.02.001
Boisnoir A, Pascal P Y, Cordonnier S, et al. 2019a. Spatio-temporal dynamics and biotic substrate preferences of benthic dinoflagellates in the Lesser Antilles, Caribbean sea. Harmful Algae, 81: 18-29. DOI:10.1016/j.hal.2018.11.012
Boisnoir A, Pascal P Y, Marro S, et al. 2019b. First spatial distribution of potentially toxic benthic dinoflagellates in the Lesser Antilles (Guadeloupe and Martinique), Caribbean Sea. Botanica Marina, 62(4): 309-322. DOI:10.1515/bot-2018-0002
Calado A J, Craveiro S C, Moestrup Ø. 1998. Taxonomy and ultrastructure of a freshwater, heterotrophic Amphidinium (Dinophyceae) that feeds on unicellular protists. Journal of Phycology, 34(3): 536-554. DOI:10.1046/j.1529-8817.1998.340536.x
Chinain M, Darius H T, Gatti C M, et al. 2016. Update on ciguatera research in French Polynesia. SPC Fisheries Newsletter, 150: 42-51.
Chinain M, Gatti C M I, Darius H T, et al. 2021. Ciguatera poisonings: a global review of occurrences and trends. Harmful Algae, 102: 101873. DOI:10.1016/j.hal.2020.101873
Chinain M, Gatti C M, Roué M et al. 2020. Ciguatera-causing dinoflagellates in the genera Gambierdiscus and Fukuyoa: distribution, ecophysiology and toxicology. In: Rao D V S ed. Dinoflagellates. Nova Science Publishers, Inc., p. 405-457.
Chomérat N, Bilien G, Viallon J, et al. 2020. Taxonomy and toxicity of a bloom-forming Ostreopsis species (Dinophyceae, Gonyaulacales) in Tahiti island (South Pacific Ocean): one step further towards resolving the identity of O. siamensis. Harmful Algae, 98: 101888. DOI:10.1016/j.hal.2020.101888
Chomérat N, Bilien G, Zentz F. 2019. A taxonomical study of benthic Prorocentrum species (Prorocentrales, Dinophyceae) from Anse Dufour (Martinique Island, eastern Caribbean Sea). Marine Biodiversity, 49(3): 1299-1319. DOI:10.1007/s12526-018-0913-6
De Azevedo Tibiriçá C E J, Sibat M, Fernandes L F, et al. 2020. Diversity and toxicity of the genus Coolia Meunier in Brazil, and detection of 44-methyl gambierone in Coolia tropicalis. Toxins, 12(5): 327. DOI:10.3390/toxins12050327
Durán-Riveroll L M, Cembella A D, Okolodkov Y B. 2019. A Review on the biodiversity and biogeography of toxigenic benthic marine dinoflagellates of the coasts of Latin America. Frontiers in Marine Science, 6: 148. DOI:10.3389/fmars.2019.00148
Faust M A. 1990. Morphologic details of six benthic species of Prorocentrum (Pyrrophyta) from a Mangrove Island, Twin Cays, Belize, including two new species. Journal of Phycology, 26(3): 548-558. DOI:10.1111/j.0022-3646.1990.00548.x
Faust M A. 1994. Three new benthic species of Prorocentrum (Dinophyceae) from carrie bow Cay, Belize: P. sabulosum sp. nov., P. sculptile sp. nov., and P. arenarium sp. nov. Journal of Phycology, 30(4): 755-763. DOI:10.1111/j.0022-3646.1994.00755.x
Faust M A. 1995. Observation of sand-dwelling toxic dinoflagellates (Dinophyceae) from widely differing sites, including two new species. Journal of Phycology, 31(6): 996-1003. DOI:10.1111/j.0022-3646.1995.00996.x
Foden J, Purdie D A, Morris S, et al. 2005. Epiphytic abundance and toxicity of Prorocentrum lima populations in the Fleet Lagoon, UK. Harmful Algae, 4(6): 1063-1074. DOI:10.1016/j.hal.2005.03.004
Fukuyo Y. 1981. Taxonomical study on benthic dinoflagellates collected in coral reefs. Nippon Suisan Gakkaishi, 47(8): 967-978. DOI:10.2331/SUISAN.47.967
Gómez F. 2003. Checklist of mediterranean free-living dinoflagellates. Botanica Marina, 46(3): 215-242.
Guillard R R L, Hargraves P E. 1993. Stichochrysis immobilis is a diatom, not a chrysophyte. Phycologia, 32(3): 234-236. DOI:10.2216/i0031-8884-32-3-234.1
Guiry M D, Guiry G M. 2022. AlgaeBase. World-Wide Electronic Publication, National University of Ireland, Galway.
Hachani M A, Dhib A, Fathalli A, et al. 2018. Harmful epiphytic dinoflagellate assemblages on macrophytes in the Gulf of Tunis. Harmful Algae, 77: 29-42. DOI:10.1016/j.hal.2018.06.006
Herrera-Sepúlveda A, Medlin L K, Murugan G, et al. 2015. Are Prorocentrum homannianum and Prorocentrum belizeanum (Dinophyceae, Prorocentrales), the samespecies? An integration of morphological and molecular data. Journal of Phycology, 51(1): 173-188. DOI:10.1111/jpy.12265
Ho T V, Nguyen L N. 2014. Morphology and distribution of the three epiphytic dinoflagellate species Coolia monotis, C. tropicalis, and C. canariensis (Ostreopsidaceae, Gonyaulacales, Dinophyceae) from Vietnamese coastal waters. Ocean Science Journal, 49(3): 211-221. DOI:10.1007/s12601-014-0021-x
Hoppenrath M, Chomérat N, Horiguchi T, et al. 2013. Taxonomy and phylogeny of the benthic Prorocentrum species (Dinophyceae)—a proposal and review. Harmful Algae, 27: 1-28. DOI:10.1016/j.hal.2013.03.006
Jang S H, Jeong H J, Yoo Y D. 2018. Gambierdiscus jejuensis sp. nov., an epiphytic dinoflagellate from the waters of Jeju Island, Korea, effect of temperature on the growth, and its global distribution.. Harmful Algae, 80: 149-157. DOI:10.1016/j.hal.2018.11.007
Jeong H J, Yih W, Kang N S, et al. 2012. First report of the epiphytic benthic dinoflagellates Coolia canariensis and Coolia malayensis in the waters off Jeju Island, Korea: morphology and rDNA sequences. . J.Eukaryot. Microbiol., 59(2): 114-133. DOI:10.1111/j.1550-7408.2012.00610.x
Jiang L C. 2019. Taxonomic Research of the Benthic Dinoflagellate Gambierdiscus and Amphidinum from the South China Sea and Studies on the Toxin Constituents of Gambierdiscus. Jinan University. p. 1-88. (in Chinese)
Kang N S, Jeong H J, Lee S Y, et al. 2013. Morphology and molecular characterization of the epiphytic benthic dinoflagellate Ostreopsis cf. ovata in the temperate waters off Jeju Island, Korea. Harmful Algae, 27: 98-112. DOI:10.1016/j.hal.2013.05.006
Karafas S, Teng S T, Leaw C P, et al. 2017. An evaluation of the genus Amphidinium (Dinophyceae) combining evidence from morphology, phylogenetics, and toxin production, with the introduction of six novel species. Harmful Algae, 68: 128-151. DOI:10.1016/j.hal.2017.08.001
Karafas S, York R, Tomas C. 2015. Morphological and genetic analysis of the Coolia monotis species complex with the introduction of two new species, Coolia santacroce sp. nov. and Coolia palmyrensis sp. nov. (Dinophyceae). Harmful Algae, 46: 18-33. DOI:10.1016/j.hal.2015.05.002
Kudela R M, Berdalet E, Enevoldsen H, et al. 2017. GEOHAB-the global ecology and oceanography of harmful algal blooms program: motivation, goals, and legacy. Oceanography, 30(1): 12-21. DOI:10.5670/oceanog.2017.106
Larsson M E, Smith K F, Doblin M A. 2019. First description of the environmental niche of the epibenthic dinoflagellate species Coolia palmyrensis, C.malayensis, and C. tropicalis (Dinophyceae) from Eastern Australia. Journal of Phycology, 55(3): 565-577. DOI:10.1111/jpy.12833
Leaw C P, Tan T H, Lim H C, et al. 2016. New scenario for speciation in the benthic dinoflagellate genus Coolia (Dinophyceae). Harmful Algae, 55: 137-149. DOI:10.1016/j.hal.2016.02.010
Lee K H, Jeong H J, Park K, et al. 2013. Morphology and molecular characterization of the epiphytic dinoflagellate Amphidinium massartii, isolated from the temperatewaters off Jeju Island, Korea. Algae, 28(3): 213-231. DOI:10.4490/algae.2013.28.3.213
Lee L K, Lim Z F, Gu H, et al. 2020. Effects of substratum and depth on benthic harmful dinoflagellate assemblages. Science Reports, 10(1): 11251. DOI:10.1038/s41598-020-68136-6
Leung P T Y, Yan M, Lam V T T, et al. 2018. Phylogeny, morphology and toxicity of benthic dinoflagellates of the genus Fukuyoa (Goniodomataceae, Dinophyceae) from a subtropical reef ecosystem in the South China Sea. Harmful Algae, 74: 78-97. DOI:10.1016/j.hal.2018.03.003
Leung P T Y, Yan M, Yiu S K F, et al. 2017. Molecular phylogeny and toxicity of harmful benthic dinoflagellates Coolia (Ostreopsidaceae, Dinophyceae) in a sub-tropicalmarine ecosystem: the first record from Hong Kong. Marine Pollution Bulletin, 124(2): 878-889. DOI:10.1016/j.marpolbul.2017.01.017
Lim A S, Jeong H J. 2021. Benthic dinoflagellates in Korean waters. Algae, 36(2): 91-109. DOI:10.4490/algae.2021.36.5.31
Lim Z F, Luo Z H, Lee L K, et al. 2019. Taxonomy and toxicity of Prorocentrum from Perhentian Islands (Malaysia), with a description of a non-toxigenic species Prorocentrum malayense sp. nov. (Dinophyceae). Harmful Algae, 83: 95-108. DOI:10.1016/j.hal.2019.01.007
Litaker R W, Vandersea M W, Faust M A, et al. 2010. Global distribution of ciguatera causing dinoflagellates in the genus Gambierdiscus. Toxicon, 56(5): 711-730. DOI:10.1016/j.toxicon.2010.05.017
Luo Z H, Zhang H, Krock B, et al. 2017. Morphology, molecular phylogeny and okadaic acid production of epibenthic Prorocentrum (Dinophyceae) species from the northernSouth China Sea. Algal Research, 22: 14-30. DOI:10.1016/j.algal.2016.11.020
Luo Z H, Zhang H, Li Q et al. 2022. Characterization of Amphidinium (Amphidiniales, Dinophyceae) speciesfrom the China Sea based on morphological, molecular, and pigment data. Journal of Oceanology and Limnology, https://doi.org/10.1007/s00343-021-1049-2.
Mangialajo L, Ganzin N, Accoroni S, et al. 2011. Trends in Ostreopsis proliferation along the Northern Mediterraneancoasts. Toxicon, 57(3): 408-420. DOI:10.1016/j.toxicon.2010.11.019
Mangialajo L, Bertolotto R, Cattaneo-Vietti R, et al. 2008. The toxic benthic dinoflagellate Ostreopsis ovata: quantification of proliferation along the coastline of Genoa, Italy. Marine Pollution Bullton, 56(6): 1209-1214. DOI:10.1016/j.marpolbul.2008.02.028
Mohammad-Noor N, Daugbjerg N, Moestrup Ø, et al. 2005. Marine epibenthic dinoflagellates from Malaysia-a study of live cultures and preserved samples based on light and scanning electron microscopy. Nordic Journal of Botany, 24(6): 629-690. DOI:10.1111/j.1756-1051.2004.tb01938.x
Monti M, Minocci M, Beran A, et al. 2007. First record of Ostreopsis cfr. ovata on macroalgae in the NorthernAdriatic Sea. Marine Pollution Bulletin, 54(5): 598-601. DOI:10.1016/j.marpolbul.2007.01.013
Murray S, Patterson D J. 2002. The benthic dinoflagellate genus Amphidinium in south-eastern Australian waters, including three new species. European Journal of Phycology, 37(2): 279-298. DOI:10.1017/s0967026202003591
Mustapa N I, Yong H L, Lee L K, et al. 2019. Growth and epiphytic behavior of three Gambierdiscus species (Dinophyceae) associated with various macroalgal substrates. Harmful Algae, 89: 101671. DOI:10.1016/j.hal.2019.101671
Nascimento S M, Da Silva R A F, Oliveira F, et al. 2019. Morphology and molecular phylogeny of Coolia tropicalis, Coolia malayensis and a new lineage of the Coolia canariensis species complex (Dinophyceae)isolated from Brazil. European Journal of Phycology, 54(3): 484-496. DOI:10.1080/09670262.2019.1599449
Nascimento S M, Mendes M C Q, Menezes M, et al. 2017. Morphology and phylogeny of Prorocentrum caipirignum sp. nov. (Dinophyceae), a new tropical toxic benthic dinoflagellate. Harmful Algae, 70: 73-89. DOI:10.1016/j.hal.2017.11.001
Nascimento S M, Neves R A F, De'Carli G A L, et al. 2020. Ostreopsis cf. ovata (Dinophyceae) molecular phylogeny, morphology, and detection of ovatoxins in strains and field samples from Brazil. Toxins, 12(2): 70. DOI:10.3390/toxins12020070
Nishimura T, Sato S, Tawong W, et al. 2013. Genetic diversity and distribution of the ciguatera-causing dinoflagellate Gambierdiscus spp. (Dinophyceae) in coastal areas ofJapan. PLoS One, 8(4): ee60882. DOI:10.1371/journal.pone.0060882
Okolodkov Y B, Campos-Bautista G, Gárate-Lizárraga I, et al. 2007. Seasonal changes of benthic and epiphytic dinoflagellates in the Veracruz reef zone, Gulf of Mexico. Aquatic Microbial Ecology, 47: 223-237. DOI:10.3354/ame047223
Pagliara P, Caroppo C. 2012. Toxicity assessment of Amphidinium carterae, Coolia cfr. monotis and Ostreopsis cfr. ovata (Dinophyta) isolated from the northern Ionian Sea (Mediterranean Sea). Toxicon, 60(6): 1203-1214. DOI:10.1016/j.toxicon.2012.08.005
Parsons M L, Aligizaki K, Bottein M-Y D, et al. 2012. Gambierdiscus and Ostreopsis: reassessment of the state ofknowledge of their taxonomy, geography, ecophysiology, and toxicology. Harmful Algae, 14: 107-129.
Parsons M L, Preskitt L B J H. 2007. A survey of epiphytic dinoflagellates from the coastal waters of the island of Hawai'i. Harmful Algae, 6(5): 658-669. DOI:10.1016/j.hal.2007.01.001
Penna A, Fraga S, Battocchi C, et al. 2010. A phylogeographical study of the toxic benthic dinoflagellate genus Ostreopsis Schmidt. Journal of Biogeography, 37(5): 830-841. DOI:10.1111/j.1365-2699.2009.02265.x
Reñé A, Satta C T, López-García P, et al. 2020. Re-evaluation of Amphidiniopsis (Dinophyceae) morphogroups based on phylogenetic relationships, and description of three new sand-dwelling species from the NW Mediterranean. Journal of Phycology, 56(1): 68-84. DOI:10.1111/jpy.12938
Rhodes L, Smith K F. 2019. A checklist of the benthic and epiphytic marine dinoflagellates of New Zealand, including Rangitāhua/Kermadec Islands. New Zealand Journal of Marine and Freshwater Research, 53(2): 258-277. DOI:10.1080/00288330.2018.1492425
Rhodes L. 2011. World-wide occurrence of the toxic dinoflagellate genus Ostreopsis Schmidt. Toxicon, 57(3): 400-407. DOI:10.1016/j.toxicon.2010.05.010
Richlen M L, Lobel P S. 2011. Effects of depth, habitat, and water motion on the abundance and distribution of ciguatera dinoflagellates at Johnston Atoll, Pacific Ocean. Marine Ecology Progress Series, 421: 51-66. DOI:10.3354/meps08854
Rodríguez F, Riobó P, Crespín G D, et al. 2018. The toxic benthic dinoflagellate Prorocentrum maculosum Faust is a synonym of Prorocentrum homannianum Faust. Harmful Algae, 78: 1-8. DOI:10.1016/j.hal.2018.06.009
Sato S, Nishimura T, Uehara K, et al. 2011. Phylogeography of Ostreopsis along west Pacific coast, with special referenceto a novel clade from Japan. PLoS One, 6(12): e27983. DOI:10.1371/journal.pone.0027983
Scholin C A, Herzog M, Sogin M, et al. 1994. Identification of group-and strain-specific genetic markers for globally distributed Alexandrium (Dinophyceae). Ⅱ. Sequence analysis of a fragment of the LSU rRNA gene. Journal of Phycology, 30(6): 999-1011. DOI:10.1111/j.0022-3646.1994.00999.x
Shears N T, Ross P M. 2009. Blooms of benthic dinoflagellates of the genus Ostreopsis; an increasing and ecologically important phenomenon on temperate reefs in New Zealand and worldwide. Harmful Algae, 8(6): 916-925. DOI:10.1016/j.hal.2009.05.003
Taylor D L. 1971. On the symbiosis between Amphidinium klebsii[Dinophyceae] and Amphiscolops langerhansi [Turbellaria: Acoela].. Journal of the Marine Biological Association of the United Kingdom, 51(2): 301-313. DOI:10.1017/S0025315400031799
Torres-Ariño A, Okolodkov Y B, Herrera-Herrera N V. 2019. A checklist of phytoplankton and microphytobenthos of the southeastern Mexican Pacific. Cymbella, 5(1): 1-97.
Tripuraneni J, Koutsouris A, Pestic L, et al. 1997. The toxin of diarrheic shellfish poisoning, okadaic acid, increases intestinal epithelial paracellular permeability. Gastroenterology, 112(1): 100-108. DOI:10.1016/s0016-5085(97)70224-5
Widiarti R, Anggrain F. 2012. The distribution of toxic dinoflagellates on sea grass Enhalus acoroides at Pari Island, Seribu Islands. Jurnal Ilmu dan Teknologi Kelautan Tropis, 4(2): 259-266. DOI:10.29244/jitkt.v4i2.7788
Yasumoto T, Nakajima I, Bagnis R, et al. 1977. Finding of a dinoflagellate as a likely culprit of ciguatera. Nippon Suisan Gakkaishi, 43(8): 1021-1026. DOI:10.2331/suisan.43.1021
Yasumoto T, Seino N, Murakami Y, et al. 1987. Toxins produced by benthic dinoflagellates. Biological Bulletin, 172(1): 128-131. DOI:10.2307/1541612
Yong H L, Mustapa N I, Lee L K, et al. 2018. Habitat complexity affects benthic harmful dinoflagellate assemblages in the fringing reef of Rawa Island, Malaysia. Harmful Algae, 78: 56-68. DOI:10.1016/j.hal.2018.07.009
Zhang H, Li Y, Cen J Y, et al. 2015. Morphotypes of Prorocentrum lima (Dinophyceae) from HainanIsland, South China Sea: morphological and molecular characterization. Phycologia, 54(5): 503-516. DOI:10.2216/15-8.1
Zhang H, Lü S H, Cen J Y, et al. 2021. Morphology and molecular phylogeny of three species of Coolia (Dinophyceae) from Hainan Island, South China Sea. Journal of Oceanology and Limnology, 39(3): 1020-1032. DOI:10.1007/s00343-020-9326-z
Zhang H, Lu S H, Li Y, et al. 2018. Morphology and molecular phylogeny of Ostreopsis cf.ovata and O. lenticularis (Dinophyceae) from Hainan Island, South China Sea. Phycological Research, 66(1): 3-14. DOI:10.1111/pre.12192
Zhang H, Wu Z, Cen J Y, et al. 2016. First report of three benthic dinoflagellates, Gambierdiscus pacificus, G.australes and G. Caribaeus (Dinophyceae), from HainanIsland, South China Sea. Phycological Research, 64(4): 259-273. DOI:10.1111/pre.12137
Zheng J P, Li Q, Lv S H, et al. 2017. Production of Palytoxins by two strains of Ostreopsis sp.(Dinophyceae) from Weizhou Island, Guangxi, China.. Ecology and Environment Sciences, 26(4): 663-670. DOI:10.16258/j.cnki.1674-5906.2017.04.017.(inChinesewithEnglishabstract)
Zou J, Li Q, Liu H et al. 2021. Taxonomy and toxin profile of harmful benthic Prorocentrum (Dinophyceae) species from the Xisha Islands, South China Sea. Journal of Oceanology and Limnology, https://doi.org/10.1007/s00343-021-1045-6.
Zou J, Li Q, Lu S H, et al. 2020. The first benthic harmful dinoflagellate bloom in China: morphology and toxicology of Prorocentrum concavum. Marine Pollution Bulletin, 158: 111313. DOI:10.1016/j.marpolbul.2020.111313