3.1 Toxin profiles
Table 2 presents the concentration data set from all the PSTs quantified in both species across the six sites. The samples collected from the Fujian coast in China were found to contain a range of PSTs analogues, including GTX1, GTX2, GTX3, GTX4, GTX5, dcSTX, STX, neoSTX, dcGTX2, and dcGTX3. In terms of toxicity proportion, the profiles were dominated by GTX5 (detection rate 4.46%), followed by neoSTX (3.48%) and dcGTX2 (2.93%). The remaining analogues were present at relatively low levels, with a detection rate of approximately 1.10–2.08%. The toxicity was higher than the European Food Safety Authority (EFSA) legal limit in 4 out of the 2355 detected samples (0.17%) (EFSA, 2009). The 2137.10 µg STXeq/kg was the strongest content, with a mean concentration from 9.59 µg/kg to 16.24 µg/kg, while the maximum content was 45.90–4080 µg/kg.
Table 2
Concentration data set of PSTs of all samples (N = 2355).
| Maximum concentration (ug/kg) | Mean concentration (ug/kg) |
STX | 239.0 | 10.01 |
dcSTX | 440.0 | 9.81 |
neoSTX | 130.0 | 9.66 |
dcGTX2 | 2220.0 | 11.72 |
dcGTX3 | 342.0 | 9.97 |
GTX1 | 162.0 | 9.76 |
GTX2 | 315.0 | 9.89 |
GTX3 | 103.0 | 9.70 |
GTX4 | 45.9 | 9.59 |
GTX5 | 4080.0 | 16.24 |
The shellfish species with a high detection rate were A. granosa, C. farreri, P. viridis, and C. gigas at 14.63%, 13.54%, 11.75%, and 10.79%, respectively. The highest concentrations were quantified in P. viridis and C. gigas from QZ and ZZ, with values over 800 µg STXeq/kg in both sites, with a maximum concentration of 2137.10 µg STXeq/kg in the P. viridis collected in June 2017 from ZZ. Significantly lower total PSTs concentrations were obtained in the C. gigas collected from QZ and ZZ, with the highest concentrations of 1469.46 and 1983.89 µg STXeq/kg, respectively. At the A. granosa and C. farreri sites, lower (< 800 µg STXeq/kg) levels than those of the P. viridis and C. gigas of PSTs were found between May 2017 and December 2021.
3.2 Temporal variation
Figure 3 shows the highest levels of PSTs in 2017, followed by 2018 and 2019, while almost no PSTs were detected in 2020 and 2021. The dominance of GTX5 was observed in 2017 (21.77%), followed by the detection rate of dcSTX (8.06%). The highest toxin level in 2018 was also the GTX5 (10.21%), followed by dcGTX2 (4.08%).The detection rate between 2019 and 2020 was less than 5%, and PSTs were not detected in 2021.
The distribution of toxins for each month in 2018 and 2019 is shown in Fig. 4. PSTs occurs mainly in the periods from May to July and November to December. GTX5 was present in almost every month, dcGTX2 was mainly present from May to August, while neoSTX showed high levels in November, December and January. The recent work of Wang et al. (2017) showed evidence of PSTs in the area of Shenzhen through the radial basis function difference model and concluded that the accumulation time of PSTs was from January to June (Wang et al., 2017). Yao et al. (2019) studied seven elements in China and found that it was rarely detected in spring (April to May) and detected in summer (June to August). This may be related to temperature; warm climate is conducive to the reproduction and secretion of toxins by poisonous phytoplankton (Onofrio et al., 2021).
3.3 Interspecific variation
PSTs undergo transformation from one form to another through different processes. Such processes include reduction, epimerization, oxidation, and desulfation, all of which may potentially result in changes to the overall shellfish toxicity (Andres et al., 2019). Six kinds of bivalve shellfish, namely, Perna viridis, S. constricta, R. philippinarum, Crassostyea gigas, C. farreri, Arca granosa, were collected from six sampling locations along the coastline of Fujian, China and analyzed and pretreated under optimized analytical conditions. The toxins, STX, neoSTX, GTX4, dominated the A. granosa samples throughout the observation period, with GTX5 observed every month in the C. gigas samples. The highest detection rate of 14.63% was found in the A. granosa samples, followed by the C. gigas, C. farreri, and P. viridis samples, with detection rates of 10.79%, 13.54%, and 11.75%, respectively, whereas the R. philippinarum samples exhibited only 2.26%. Moreover, differences among species were observed in the toxins content in this research. The PSTs content exceeding 800 µg STXeq/kg was recorded in 2 out of the 6 marine organism species from 6 sites, including the C. gigas and P. viridis mollusks samples.
The mean concentration of all the detected PSTs in bivalve whole soft tissues collected in five years in the study area was obtained. Figure 5 illustrates that mean concentration of GTX5 in S. constricta (12.39 µg/kg) and neoSTX, dcSTX, dcGTX2&3, and GTX1-4 have the same mean concentration (10 µg/kg). The detected samples in other shellfish A. granosa, C. gigas, C. farreri, and P. viridis were higher than that of R. philippinarum. Two toxin analogues (STX, GTX2) were rich in A. granosa, whereas the profile observed in six cockle samples appeared to have a lower mean concentration of GTX5 than that found on average in the C. gigas data set. The C. gigas samples exhibited the dominance of the GTX5 toxin (22.25 µg/kg), followed by dcGTX2 (mean concentration 14.01%). GTX5 was the highest mean concentration (14.04 µg/kg), followed by dcGTX2&3 and GTX1 (mean concentration 10.54, 10.51, and 10.83 µg/kg, respectively), however, the remainder were lower levels (9 ~ 10 µg/kg) in C. farreri. Also, mean concentration of GTX5 (13.04 µg/kg) present in the Perna viridis was dominant toxin, followed by dcSTX, dcGTX2&3. S. constricta could not be compared because of the values below the detection limit. The PST components detected in shellfish toxins are mainly carbamates, and the toxic components of shellfish are related to the species that feed on poisonous algae, which are usually similar to the PST components contained in shellfish and the toxic components of poisonous algae. However, due to the metabolism and transformation of the toxins in shellfish, the proportion and components of the toxins in shellfish will change.
3.4 Spatial variation
Significant variation in PSTs was quantified between different geographical regions within Fujian province. Figures 6 and 7 illustrate the total PSTs quantified in the sampled areas, namely, FZ, ND, PT, QZ, XM, and ZZ, each progressively exhibiting increased variation. More than half of the study samples were harvested from the most ND region of Fujian Province (21.61%). From these, the highest toxicity was 403.26 µg STXeq/kg, with no samples above the RL (800 µg STXeq/kg). The highest toxicities were observed in the sampled areas of QZ and ZZ, with the latter of the three containing the samples showing the highest toxicities out of all the samples studied (toxic equivalent = 2137.10, 1835.25, and 1983.89 µg STXeq/kg respectively), with 0.17% of the total samples above the RL. Although shellfish products in every city were contaminated with at least some of paralytic marine toxins, the toxins concentration and composition varied widely among the various cities.
The four cities (FZ, QZ, XM, and ZZ) showed a similar trend of paralytic marine biotoxin contamination, and GTX5 is the most ubiquitous toxin. Most of the paralytic marine biotoxins, except for GTX4, were present in many specimens from two cities, namely, ND and XM. The detection rate of FZ was 4.23%, and the maximum concentration was 344 µg/kg, but STX, neoSTX, dcSTX, and GTX3 were not found in the area. neoSTX was the most detected in ND, accounting for 71.74% of the total toxins, and the maximum concentration was 330 µg/kg. Furthermore, dcGTX3 was not discovered in PT, STX, neoSTX, and GTX4 were not recognized in QZ. GTX5 was the most prevalent toxins, accounting for 41.33% of the total toxins in this area. The result suggests that the detection rate of ND and QZ is high, and the difference in detection rate is insignificant, while the detection rate of FZ is much lower than that in the other areas.
Figure 8 summarizes the proportion of toxin analogues quantified in terms of STX equivalents in the bivalve mollusks harvested across each of the six regional areas along the Southeastern China coast. The figure shows that most of the toxin analogues were detected in May to August, occasionally in November to December, and never quantified in February to March. PSTs were detected in XM and PT mainly in May to June, accounting for 33–67% of the total toxicity, followed by QZ and FZ, accounting for 10–20% of the total toxicity. The mollusks from ND exhibited low relative proportions in January to November, displaying the highest relative proportion (42%).