Next Article in Journal
Effects of a Single Session of OnabotulinumtoxinA Therapy on Sleep Quality and Psychological Measures: Preliminary Findings in a Population of Chronic Migraineurs
Next Article in Special Issue
Detection of Extremely Low Level Ciguatoxins through Monitoring of Lithium Adduct Ions by Liquid Chromatography-Triple Quadrupole Tandem Mass Spectrometry
Previous Article in Journal
Body Temperature Drop as a Humane Endpoint in Snake Venom-Lethality Neutralization Tests
Previous Article in Special Issue
Ciguatera Fish Poisoning in the Caribbean Sea and Atlantic Ocean: Reconciling the Multiplicity of Ciguatoxins and Analytical Chemistry Approach for Public Health Safety
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Measuring Marine Biotoxins in a Hypersaline Coastal Lagoon

Department of Chemical and Environmental Engineering, Technical University of Cartagena (UPCT), E-30203 Cartagena, Spain
*
Author to whom correspondence should be addressed.
Toxins 2023, 15(9), 526; https://doi.org/10.3390/toxins15090526
Submission received: 25 June 2023 / Revised: 14 August 2023 / Accepted: 20 August 2023 / Published: 26 August 2023
(This article belongs to the Special Issue Analytical Chemistry Techniques in Toxin Detection)

Abstract

:
Marine biotoxins have posed a persistent problem along various coasts for many years. Coastal lagoons are ecosystems prone to phytoplankton blooms when altered by eutrophication. The Mar Menor is the largest hypersaline coastal lagoon in Europe. Sixteen marine toxins, including lipophilic toxins, yessotoxins, and domoic acid (DA), in seawater samples from the Mar Menor coastal lagoon were measured in one year. Only DA was detected in the range of 44.9–173.8 ng L−1. Environmental stressors and mechanisms controlling the presence of DA in the lagoon are discussed. As an enrichment and clean-up method, we employed solid phase extraction to filter and acidify 75 mL of the sample, followed by pre-concentration through a C18 SPE cartridge. The analytes were recovered in aqueous solutions and directly injected into the liquid chromatography system (LC-MS), which was equipped with a C18 column. The system operated in gradient mode, and we used tandem mass spectrometry (MS/MS) with a triple quadrupole (QqQ) in the multiple reaction monitoring mode (MRM) for analysis. The absence of matrix effects was checked and the limits of detection for most toxins were low, ranging from 0.05 to 91.2 ng L−1, depending on the compound. To validate the measurements, we performed recovery studies, falling in the range of 74–122%, with an intraday precision below 14.9% RSD.
Key Contribution: The concentrations of marine biotoxins were evaluated in a hypersaline coastal lagoon; one toxin was detected at low concentration levels; phosphate limitation against silicates and dilution drive the amount of domoic acid toxin in the lagoon.

1. Introduction

Certain marine phytoplankton species produce phycotoxins that are causing serious health problems for both humans, by triggering epidemic outbreaks, and ecosystems, by altering their structure. These toxins can bio-accumulate in the food chain [1,2,3,4,5], and humans can ingest them when they accumulate in filter feeders, such as shellfish or other seafood [6,7]. To ensure food safety, many countries’ food safety regulators have established maximum tolerable levels for these toxins in seafood [8,9,10], and as a result, monitoring and control programs for harmful algae blooms (HABs) are typically established in aquaculture and shellfish areas. However, the impact of HABs is not limited to human health only; it also affects the economy when mass fish mortality occurs in aquaculture areas or precautionary closures, also in touristic areas, are implemented. These events carry an annual estimated cost of 3–4 billion USD worldwide [11], with millions of dollars in one single event [12,13]. Harmful algal blooms have a significant impact not only on the aquaculture and fisheries sectors but also on the quality of recreational waters. The tourism industry, in particular, is highly vulnerable to this problem. Beach closures resulting from HABs, such as those that occurred between 2005 and 2009 in popular tourist destinations along the Mediterranean coasts of France and Italy [14,15,16,17,18] due to Ostreopsis outbreaks, caused significant economic losses. While most attention is focused on monitoring the effects of biotoxins in seafood, less attention is given to recreational waters where the quality of the bathing water is typically determined by the presence of fecal coliforms or intestinal enterococci only [19,20].
Although the food safety agencies have regulated the maximum levels of biotoxins in seafood, with lipophilic marine toxins (LMTs) mainly studied in contaminated bivalves, there is no regulation for toxin concentrations in bathing waters, and the levels of LMTs in seawater remain largely unexplored. Special attention has been paid to domoic acid, a potent amnesic toxin [21,22], which has been reported in the literature in several members of the marine food chain, from small organisms like copepods and krill to top predators such as seabirds and sea lions [23,24,25,26].
The objective of this study is to measure these toxins in a hypersaline coastal lagoon that is experiencing severe eutrophication to better understand the processes keeping them in the water. There are a few records of in-water marine biotoxins in Mediterranean coastal lagoons [27,28,29,30,31]. The Mar Menor is the largest hypersaline coastal lagoon in Europe, and is a popular tourist area, with large fisheries, and although it is not a shellfish cultivation area, cockles and sea snails are frequently consumed by local residents. To perform this study, we used solid phase extraction (SPE) coupled to liquid chromatography (LC)–tandem mass spectrometry (MS/MS) with a triple quadrupole (QqQ) in multiple reaction monitoring (MRM) mode. Based on existing methods, we adjusted different parameters to optimize the measurement of sixteen marine toxins in seawater, including lipophilic marine toxins, yessotoxins, and domoic acid (DA) in the Mar Menor coastal lagoon for the first time over the course of one year. However, the only recorded toxin in the water was DA at trace (ng·L−1) level.
The use of LC–MS methodologies for the detection of phycotoxins in these types of waters is currently on the rise as an alternative strategy to determine the presence of toxins and their possible implications in the ecosystem [32,33,34,35,36]. However, the analysis of in-water toxins using LC-MS systems requires pre-treatment of samples, including desalting and clean-up steps, as well as preconcentration techniques such as solid phase extraction (SPE) to achieve detection at low levels [1,37,38,39,40,41,42].

2. Results and Discussion

2.1. Enhancements in the LC-MS Analytical Procedure

To separate the toxins, a C18 stationary phase was utilized with a mixture of 2 mM ammonium acetate and 0.1% FA in water and MeOH as the mobile phase, based on the literature [27,32,38]. Three LC methods were explored using Zorbax SB-C18 column for DA, YTXs, and the remaining lipophilic toxins.
The percentage of MeOH in the mobile phase was examined to elute DA through the column. After selecting 25% MeOH, the percentage was increased to 95% to clean the column after each injection, and then the initial conditions were restored.
The separation of the remaining toxins was optimized by changing the percentage of MeOH at different times to achieve good resolution of all compounds in a shorter method. Temperature also affected the peak shapes, and experiments conducted in the range of 20–50 °C indicated that peaks approached each other as the temperature increased. Therefore, the temperature was adjusted to 30 °C to maintain constant column conditions.
The retention times (tR) for all compounds, along with their multiple reaction monitoring (MRM) transitions, are listed in Table 1. The MS/MS parameters were adopted from previous works [27,32,38]. YTX and hYTX detection was included and showed maximum sensitivity in negative ion mode (ESI-) with the deprotonated molecular ion [M-H].

2.2. Enhancements in the SPE Procedure

A silica C18 stationary phase, as in the LC column, was used for the SPE cartridge, which allowed analytes to be retained from a large volume sample and preconcentrated in a small elution fraction. To optimize the procedure, a 2 ng/mL aqueous solution of toxins and a 200 mg, 10 mL Bond Elut C18 cartridge were employed. Modifications of the De la Iglesia et al. [38] protocol for DA were made to adjust it for lipophilic toxins to employ a same method for the 16 compounds. To study the changes, every solution that passed through the cartridge was analyzed using the LC-QqQ/MS system.
The SPE technique depends on multiple variables, which were optimized based on methods described in the literature. The SPE cartridges were conditioned with 20 mL of MeOH and 20 mL of 0.1% FA. The aqueous standard solution (50 mL) was acidified with 0.1% FA to keep DA with a neutral charge for retention. After loading the sample, a washed step was applied with 20 mL of 0.1% FA. Then, three elutions were performed using 3 mL of MeOH:water (1:9, v/v) 0.2 M ammonium acetate. The results showed that toxins were retained in the C18 phase, even after washing. Moreover, DA appeared only in the first elution, meaning that 3 mL was sufficient to recover it. However, lipophilic toxins appeared in all elutions, so their elution would require a higher percentage of MeOH or a higher concentration of ammonium acetate.
The percentage of MeOH was optimized from 10 to 100, with the recovery of DA decreasing (Figure 1A) and that of lipophilic toxins increasing (Figure 2A,B).
Then, two basic solutions, ammonium acetate and ammonia, were tested at several concentrations and combined with the organic phase. As DA and lipophilic toxins need to be eluted in a different percentage of the organic phase, two elution fractions (E1 and E2) were used. E1 contained 10% MeOH and 0.2 M ammonium acetate, while E2 contained 100% MeOH.
For E1 elution, three volumes (1, 2, and 5 mL) were tested, and the recoveries obtained for DA were 95%, 92%, and 78%, respectively. Therefore, elution E1 was optimized as 1 mL of MeOH:water (1:9, v/v) 0.2 M ammonium acetate.
Elution E2 was modified with 0.2 M and 1 M of ammonium acetate, and different amounts of NH3 (0.25, 0.50, and 1%). The compound signals decrease when the percentage of NH3 increases, and the optimal condition is with 2 mL of MeOH and 1 M of ammonium acetate.
Finally, the volume of seawater was studied in the 25–100 mL range. The sample was fortified with toxins at a concentration of 20 ng/mL. Peak areas were larger with larger seawater volumes for some toxins but not for others, so 75 mL was selected because linearity was maintained for all toxins.

2.3. Validation of the Measurements

To assure the validity of the measurements, the limits of quantification and detection were calculated and analyses for selectivity, linearity, recovery, and precision were performed [43,44].
Three calibration graphs for every analyte were created using the SPE procedure combined with LC-QqQ-MS/MS by least-squares linear regression analysis of the compound concentrations versus peak area, using five concentration levels. Table 2 presents the analytical characteristics of the optimized protocol, including the slopes obtained from calibration graphs using an aqueous standard solution, natural seawater sample, or synthetic seawater. ANOVA test results showed that there were no statistically significant differences in the slopes, indicating the absence of a matrix effect. Therefore, an aqueous standard addition is recommended for quantification.
The sensitivity of the procedure was evaluated by calculating the limits of quantification (LOQs) and detection (LODs) for a signal-to-noise ratio of 10 and 3, respectively. Table 2 shows that the LODs varied from 0.05 to 91.2 ng L−1 for 13,19-didesM and hYTX, respectively.
To determine the repeatability of the method, we calculated the RSD (average relative standard deviation) from 10 replicate analyses of a seawater sample fortified with the analytes at 100 ng L−1. The RSD intraday values ranged from 0.44% to 14.9% (Table 2).
The recovery of the SPE procedure was evaluated by fortifying three seawater samples with the analytes at two concentration levels in triplicate (Table 3).
This procedure was applied to the analysis of the toxin contents of 21 seawater samples. There were no interfering peaks at the retention times of the compounds. The analytes were identified by comparing the retention times, their transitions, and the transition ratios between the samples and those obtained from standard solutions. As mentioned above, only DA was found at levels above its LOD.

2.4. Domoic Acid and Environmental Variables

The evolution of the DA, the only toxin recorded, content in the water is shown in Figure 3 (yellow pale bars), ranging from 1.4 to 224.06 ng L−1. The density of all species of the genus Pseudo-nitzschia (blue diamonds) ranged from 11 to 9.6·106 cel/L. Finding a good correlation between producer cells and DA is not an easy task due to the many other processes involved as sources of variability. It has been reported that gene expression for DA production varies with environmental conditions [45,46]. For example, some species of Pseudo-nitzschia can produce higher levels of DA at salinities in the range found here [33], but also when lacking Fe [47,48], which is not the case as the Mar Menor is highly affected by historical mining waste [49,50], under nitrate from high groundwater discharge [51], under the limitation of Si and P [52,53,54,55], and with a varying pH [56].
The overlapping of the silicates to phosphate molar ratio (Si:P), which is an expression of phosphorous limitation against silicates used by diatoms to form its valves, suggests a positive correlation between Si:P and DA, although more data are required to confirm this hypothesis in our system. Furthermore, the life cycle for Pseudo-nitzschia can last from 2–3 days to several weeks [57,58], while the toxin can remain stable in the water for several weeks to a few months [59]. Moreover, not all Pseudo-nitzschia species produce the same amount of toxin, and they do not do so under the same exact environmental conditions.
The sediments in the Mar Menor are mostly clay and mud covered with extensive macroalgae meadows of Caulerpa prolifera [60,61]. These algae uptake large amounts of nitrogen that, once transformed into biomass, accumulate in the seabed when decomposition of the meadows occurs, mainly in the winter with the low temperatures. Moreover, these meadows act as sediment traps in calm conditions after storms, thus accumulating organic matter in the sediment. Once organic matter decomposes in the anoxic sediment, a large reservoir of ammonia and phosphate is retained in its interstitial waters. Therefore, resuspension is a key process that brings back ammonia and phosphates to the water column, thus fueling phytoplankton growth.
Several pieces of evidence also indicate that DA can accumulate in the sediment [62,63]. If this were the case in the Mar Menor, as it seems to be (Velazquez, pers.com.), resuspension by waves in the very shallow waters would bring this compound back to the water column, thus imposing a larger source for both spatial and temporal variability in the water column.
At the same time, a large spatial variability exists in the phytoplankton of the Mar Menor, which tends to grow forming patches of different densities according to the very local environmental conditions [64], namely nutrients, light field, and turbulence. Advective processes add even more variability to both spatial and temporal scales, being of vital importance in redistributing the phytoplankton patches over the lagoon. As a proxy for this advective process, we have used precipitation (vertical grey bars) and salinity (cyan dots and lines) in Figure 3.
The average annual precipitation in the Mar Menor is 300 mm/year, but much of it falls in torrential downpours over a few days each year. When a precipitation event occurs, the lagoon receives runoff from the water collected in its watershed, which spans 1244 km2. Although this freshwater tends to flow to the Mediterranean Sea by sliding over the denser water in the Mar Menor, the short-term effect of rainfall on salinity is significant. The overlap of precipitation with DA evolution over time suggests that precipitation dilutes and removes DA from the water column. Furthermore, the sampling station is located in the area of influence of the Las Encañizadas inlet, which, due to its shallowness and large cross-section, plays a key role in the exchange of water between the lagoon and the Mediterranean Sea [65]. Using salinity at the sampling station as a proxy for water exchange, one can see that several events of large amounts of clean Mediterranean water entering the lagoon cause the salinity to drop, and consequently, DA tends to decrease or even disappear completely, thus reinforcing the idea of the dilution effect as a regulating factor for this toxin in this environment.

3. Conclusions

LC-QqQ-MS/MS was applied for the first time to analyze lipophilic, yessotoxins, and domoic acid marine biotoxins in the Mar Menor hypersaline coastal lagoon. The use of this system in combination with the solid phase extraction (SPE) technique allowed for preconcentration at levels at ng L−1 and cleaning and desalting of the samples. The determination of the analytes was possible without interferences or matrix effects by employing the multiple reaction monitoring (MRM) mode and standards to determine the specific retention times. This optimized procedure is applicable for monitoring these toxins in hypersaline water.
Although no clear correlation between the concentrations of phytoplankton producers and DA exists due to the many different and large resources of both temporal and spatial variability, data suggests that phosphate limitation against silicates may be a driver for DA in the lagoon. On the other hand, the dilution effect, either by rainfall or by clean water entering the lagoon from the Mediterranean Sea, can play an important role in regulating DA concentrations in the Mar Menor system.

4. Material and Methods

4.1. Study Area

The Mar Menor is a microtidal lagoon located on the southeastern coast of Spain (Figure 4). It covers an area of 175 km2, with a maximum depth of 6.5 m (an average of 3.5 m), making it the largest hypersaline lagoon in Europe with a salinity range of 38.5 to 47.5 P.S.U. During the summer season, the highest mean temperature recorded is 32 °C, and the average annual evaporation exceeds precipitation.
Three inlets connect the lagoon with the Mediterranean Sea: Las Encañizadas, El Estacio, and Marchamalo. Las Encañizadas is a complex of channels that can be either dry or wet, depending on the sea water level. It has a cross-sectional area of 163 m2 (40% of the lagoon’s total cross-sectional area) and plays a crucial role in the exchange of water during storms. Due to its shallowness, its influence is greater during high water levels, mainly due to atmospheric pressure rather than tides, which are very small in the Mediterranean.
Since 2016, the Mar Menor lagoon has been experiencing severe intermittent eutrophication episodes [66,67] resulting in several anoxic events that have caused mass mortality of fishes in 2019 and 2021. Intensive agriculture in the watershed covering more than 1244 km2, urban development around the lagoon with deficiencies in the sewage system, and poor management practices in the groundwater level have resulted in high levels of nitrates and phosphates entering the lagoon. In 2018, a sampling station (Figure 4) in the northern part of the lagoon performed regular sampling throughout the year.

4.2. Environmental Variables

Precipitation data were obtained from the Spanish Met Agency (AEMET) station at San Javier Airport (Figure 4). Water temperature and salinity profiles were measured using a CTD (Castaway, Sontek - Xylem, WA, USA). Water samples were collected using a Lund’s tube sampler [68], which integrated the entire water column from surface to bottom. Three sub-samples were collected and filtered through GF/F filters of 25 mm diameter using a syringe, and the filtrate was stored at −24 °C until the analysis of nutrients. Another three subsamples were kept at 4 °C until the analysis for toxins. One more subsample was stored in 250 mL bottles and fixed with acid-lugol until further analysis.
The analysis of nutrients (nitrate, nitrite, ammonium, phosphate, and silicate) was carried out using a high-precision colorimeter on an AA3 segmented flow autoanalyzer (Seal Analytical, Germany) following the [69] method.
To estimate cell abundance and identify taxa, the fixed samples were sedimented in Utermöhl chambers [70] for a minimum of 24 h. The samples were then counted using an inverted microscope (Leica DM-IL) equipped with a 63× magnification objective.
A summary of the environmental variables recorded is given in Table 4.

4.3. LC-MS Analytical Method

Analytical methods of marine toxins include in vitro assays, the most widely used method of mouse bioassays and chromatographic methods [71,72], coupled with different detectors, such us fluorescence [37,73] and ultraviolet [39,74], and mass spectrometry (MS) [1,27,34,38,41,42,74].
To obtain an analytical method able to analyze toxins, which are expected to be in a concentration of the order of parts per trillion, it is necessary to apply a preconcentration stage to the samples, such as the SPE technique to eliminate interferences, the seawater matrix effect, and to simultaneously analyze a wide range of concentrations of several toxins.
Other extraction techniques, such as magnetic solid phase extraction (MSPE) [75], use of a molecularly imprinted polymer (MIP) [39,76], dispersive liquid–liquid microextraction (DLLME) [27], and dispersive micro solid phase extraction (DMSPE) [77], allow only the analysis of a small group of toxins.

4.4. Chemicals and Materials

High quality acetonitrile, methanol (MeOH) and water of LC-MS quality were obtained from Chem-Lab NV (Zedelgem, Belgium). Analytical standards of 13,19-didesmethyl spirolide C (13,19-didesM), 13-desmethyl spirolide C (13-desM), 20-methyl spirolide G (SPX20G), Azaspiracid 1 (AZA1), Azaspiracid 2 (AZA2), Azaspiracid 3 (AZA3), Azaspiracid 4 (AZA4), Azaspiracid 5 (AZA5), Dinophysistoxin 1 (DTX1), Dinophysistoxin 2 (DTX2), Okadoic acid (OA), Yessotoxin (YTX), and Hommo-Yessotoxin (hYTX) were obtained from Laboratorio Cifga, S.A. (Lugo, Spain). Pectenotoxin 2 (PTX2) and Gymnodimine (GYM) were purchased from National Research Council, Institute for Marine Biosciences (NRC CNRC, Halifax, NS, Canada) and DA form Sigma–Aldrich (St. Louis, MO, USA).
Other reagents that were required, such us formic acid (FA) (98% purity), were provided by Panreac (Barcelona, Spain). Water was deionized through a Milli-Q water purification system (Billerica, MA, USA).

4.5. Sampling and Pre-Treatment

Samples were collected in 500 mL capacity polystyrene container and stored at 4 °C. A 75 mL volume sample was filtered using a GF/F filter and acidified with 0.1% FA; the SPE procedure was then applied.
A VisiprepTM SPE Vacuum Manifold from Supelco (Sigma–Aldrich) was used to preconcentrate analytes in Bond Elut LRC C18 SPE cartridges, 200 mg, 10 mL from Agilent Technologies. Bond Elut C18 cartridge was conditioned with 20 mL of MeOH and 20 mL of 0.1% FA. The pretreated sample was passed through the cartridge at a flowrate of 1 mL min−1. The cartridge was washed with 20 mL of 0.1% FA and two solutions (E1 and E2) were needed to elute the analytes. Firstly, 1 mL of MeOH:water (1:9, v/v) and 0.2 M ammonium acetate (E1) preconcentrated the DA, and the lipophilic toxins were eluted with E2 (2 mL of methanol and 1 M ammonium acetate). Both fractions were filtered employing a 0.45 µm nylon filter and were transferred to two LC amber vials and, finally, injected into the LC system.
For recovery studies, three seawater samples were fortified at two concentration levels in triplicate at 5 and 10 ng L−1 of AZAs, 20 and 50 of SPX, YTXs, DTX1, and PTX2, 25 and 50 ng L−1 of DTX2, 50 and 100 ng L−1 of DA, 200 and 400 ng L−1 of OA, and 500 and 1000 ng L−1 of GYM.

4.6. Instrumentation

The analysis was performed using an HPLC Agilent 1260 Infinity system equipped with a quaternary pump (G1311B) and an autosampler (G1329B), with an 8 µL injection volume. The analytical column utilized for the reversed-phase technique was a Zorbax SB-C18 (2.1×75 mm, 3.5 µm) from Agilent Technologies, held in an oven at 30 °C with a flow rate of 0.3 mL/min, operating in gradient mode. For positive ionization mode (ESI+), the mobile phase consisted of a solvent mixture of 2 mM ammonium acetate and 0.1% FA in water (solvent A) and MeOH (solvent B). For negative mode (ESI-), the mobile phase comprised 2 mM ammonium acetate in MeOH (solvent C) and water (solvent D), as described in previous studies [27].
The SPE technique was used with two elution fractions (E1 and E2) to preconcentrate the toxins. The E1 elution was used to recover DA, while lipophilic toxins were extracted in the E2 elution. Since the E1 elution contains only one target compound, a different gradient program is required for each elution. To elute the DA toxin, the gradient program started at 25% B for 3 min, increased to 95% B in 4 min, and was maintained for 4 min, then returned to the initial conditions in 3 min. To separate lipophilic toxins, the gradient program started at 25% B for 2.5 min, increased to 60% B in half a min, and was maintained for 5 min, increased to 75% B and was maintained for 7 min, then the column was cleaned by establishing the percentage of B at 95 for 4 min, and initial conditions were recovered in 30 min.
The E2 elution also contains yessotoxins that are ionized in negative mode, requiring a different gradient program: the initial mobile phase consisted of 70% of C and 30% of D for 2 min, then the percentage of C was increased to 95 in half a minute and was maintained for 5 min, and initial conditions were recovered in 15 min.
The MS used was an Agilent G6410A triple quadrupole. The Agilent Mass Hunter Data Acquisition (Qualitative and Quantitative Analysis) was used for data acquisition and method development. Electrospray source ionization (ESI) was operated in negative mode for YTX and hYTX and in positive ion mode for the remaining compounds. Nitrogen was used as the nebulizer and collision gas at the following conditions: nebulizer pressure 40 psi, capillary voltage 4500 V, drying gas flow 8 L/min, and temperature 350 °C.
To optimize MRM transitions, each compound was injected into the LC-MS/MS-QqQ system at a concentration close to 1 µg/mL. Firstly, mass spectra were set in the 80–1000 amu m/z range, and the precursor ion was selected. Then, product ions were selected by applying different collision energies (CE). Finally, fragmentor voltages were studied from 70 up to 200 V. The selected conditions are shown in Table 1. The most sensitive transition was used for quantification purposes, and the rest of the MRM transitions were used for the identification of the compounds.

Author Contributions

A.O.-R. and J.G. conceived the study. A.O.-R. performed the chemical analysis and toxins data analysis. N.A.-O. performed phytoplankton analysis and processed environmental data. G.F. conducted sample collection and lab processing. A.O.-R., N.A.-O. and J.G. wrote the manuscript. All authors revised and edited the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

Partial financial support was received from Consejería de Agua, Agricultura, Ganadería, Pesca, Medio Ambiente y Emergencias of the Comunidad Autónoma de Murcia (Spain).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

On request to the corresponding author.

Acknowledgments

We thank members of the Ecosystem Research Group at the Technical University of Cartagena (UPCT) for helping with the field work and data access. We also thank an anonymous reviewer for helpful and constructive comments to a first manuscript.

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Orellana, G.; Van Meulebroek, L.; De Rijcke, M.; Janssen, C.R.; Vanhaecke, L. High resolution mass spectrometry-based screening reveals lipophilic toxins in multiple trophic levels from the North Sea. Harmful Algae 2017, 64, 30–41. [Google Scholar] [CrossRef] [PubMed]
  2. Turner, J.T.; Tester, P.A. Toxic marine phytoplankton, zooplankton grazers, and pelagic food webs. Limnol. Oceanogr. 1997, 42, 1203–1213. [Google Scholar] [CrossRef]
  3. Doucette, G.; Maneiro, I.; Riveiro, I.; Svensen, C. Phycotoxin Pathways in Aquatic Food Webs: Transfer, Accumulation and Degradation. In Ecological Studies; Springer: Berlin/Heidelberg, Germany, 2006; Volume 189, pp. 283–295. [Google Scholar]
  4. Li, J.; Ruan, Y.; Mak, Y.L.; Zhang, X.; Lam, J.C.W.; Leung, K.M.Y.; Lam, P.K.S. Occurrence and trophodynamics of marine lipophilic phycotoxins in a subtropical marine food web. Environ. Sci. Technol. 2021, 55, 8829–8838. [Google Scholar] [CrossRef]
  5. Lefebvre, K.A.; Fachon, E.; Bowers, E.K.; Kimmel, D.G.; Snyder, J.A.; Stimmelmayr, R.; Grebmeier, J.M.; Kibler, S.; Ransom Hardison, D.; Anderson, D.M.; et al. Paralytic shellfish toxins in Alaskan Arctic food webs during the anomalously warm ocean conditions of 2019 and estimated toxin doses to Pacific walruses and bowhead whales. Harmful Algae 2022, 114, 102205. [Google Scholar] [CrossRef]
  6. Fire, S.E.; Wang, Z.; Berman, M.; Langlois, G.W.; Morton, S.L.; Sekula-Wood, E.; Benitez-Nelson, C.R. Trophic transfer of the harmful algal toxin domoic acid as a cause of death in a minke whale (Balaenoptera acutorostrata) stranding in southern California. Aquat. Mamm. 2010, 36, 342–350. [Google Scholar] [CrossRef]
  7. Costa, P.R.; Garrido, S. Domoic acid accumulation in the sardine Sardina pilchardusand its relationship to Pseudo-nitzschia diatom ingestion. Mar. Ecol. Prog. Ser. 2004, 284, 261–268. [Google Scholar] [CrossRef]
  8. The European Parliament and the Council of the European Union. Regulation (EC) No 853/2004 of the European Parliament and of the Council of 29 April 2004 laying down specific hygiene rules for food of animal origin foodstuffs. Off. J. Eur. Union 2004, 139, 55–64. [Google Scholar]
  9. Vilariño, N.; Louzao, M.C.; Fraga, M.; Botana, L.M. From Science to Policy: Dynamic Adaptation of Legal Regulations on Aquatic Biotoxins. In Climate Change and Marine and Freshwater Toxins; Botana, L.M., Louzao, M.C., Vilarino, N., Eds.; De Gruyter Mouton: Berlin, Germany, 2020; pp. 607–654. ISBN 9783110625738. [Google Scholar]
  10. FAO. Marine Biotoxins. In Food and Nutrition Paper 80; FAO: Rome, Italy, 2004. [Google Scholar]
  11. Brown, A.R.; Lilley, M.; Shutler, J.; Lowe, C.; Artioli, Y.; Torres, R.; Berdalet, E.; Tyler, C.R. Assessing risks and mitigating impacts of harmful algal blooms on mariculture and marine fisheries. Rev. Aquac. 2020, 12, 1663–1688. [Google Scholar] [CrossRef]
  12. Sanseverino, I.; Conduto, D.; Pozzoli, L.; Dobricic, S.; Lettieri, T. Algal Bloom and Its Economic Impact; European Commission, Joint Research Centre Institute for Environment and Sustainability: Ispra, Italy, 2016. [Google Scholar]
  13. Díaz, P.A.; Álvarez, G.; Varela, D.; Pérez-Santos, I.; Díaz, M.; Molinet, C.; Seguel, M.; Aguilera-Belmonte, A.; Guzmán, L.; Uribe, E.; et al. Impacts of harmful algal blooms on the aquaculture industry: Chile as a case study. Perspect. Phycol. 2019, 6, 39–50. [Google Scholar] [CrossRef]
  14. Brescianini, C.; Grillo, C.; Melchiorre, N.; Bertolotto, R.; Ferrari, A.; Vivaldi, B.; Icardi, G.; Gramaccioni, L.; Funari, E.; Scardala, S. Ostreopsis ovata algal blooms affecting human health in Genova, Italy, 2005 and 2006. Eurosurveillance 2006, 11, 3040. [Google Scholar] [CrossRef]
  15. Durando, P.; Ansaldi, F.; Oreste, P.; Moscatelli, P.; Marensi, L.; Grillo, C.; Gasparini, R.; Icardi, G. Ostreopsis ovata and human health: Epidemiological and clinical features of respiratory syndrome outbreaks from a two-year syndromic surveillance, 2005–2006, in north-west Italy. Eurosurveillance 2007, 12, 2005. [Google Scholar] [CrossRef]
  16. Ciminiello, P.; Dell’Aversano, C.; Fattorusso, E.; Forino, M.; Tartaglione, L.; Grillo, C.; Melchiorre, N. Putative palytoxin and its new analogue, ovatoxin-a, in Ostreopsis ovata collected along the ligurian coasts during the 2006 toxic outbreak. J. Am. Soc. Mass Spectrom. 2008, 19, 111–120. [Google Scholar] [CrossRef]
  17. Kermarec, F.; Dor, F.; Armengaud, A.; Charlet, F.; Kantin, R.; Sauzade, D.; De Haro, L. Les risques sanitaires liés à la présence d’Ostreopsis ovata dans les eaux de baignade ou d’activités nautiques. Environ. Risques Et. Sante 2008, 7, 357–363. [Google Scholar]
  18. Tichadou, L.; Glaizal, M.; Armengaud, A.; Grossel, H.; Lemée, R.; Kantin, R.; Lasalle, J.L.; Drouet, G.; Rambaud, L.; Malfait, P.; et al. Health impact of unicellular algae of the Ostreopsis genus blooms in the Mediterranean Sea: Experience of the French Mediterranean coast surveillance network from 2006 to 2009. Clin. Toxicol. 2010, 48, 839–844. [Google Scholar] [CrossRef] [PubMed]
  19. The European Parliament and the Council of the European Union. European Parliament Directive 2006/7/EC of the European Parliament and of the Council of 15 February 2006 concerning the management of bathing water quality and repealing Directive 76/160/EEC. Off. J. Eur. Union 2006, 64, 37–51. [Google Scholar]
  20. USEPA. Water quality criteria. In Water Quality Standards Handbook; Office of Water: Washington, DC, USA, 2017; pp. 1–26. ISBN EPA-823-B-17-001. [Google Scholar]
  21. EFSA. Panel on contaminants in the food chain (CONTAM) Marine biotoxins in shellfish: Summary on regulated marine biotoxins. EFSA J. 2009, 1306, 1–23. [Google Scholar]
  22. EFSA. Panel on Contaminants in the Food Chain Marine biotoxins in shellfish—Domoic acid. EFSA J. 2009, 1181, 1–61. [Google Scholar]
  23. Work, T.M.; Barr, B.; Beale, A.M.; Fritz, L.; Michael, A.; Wright, J.L.C.; Url, S. Epidemiology of domoic acid poisoning in brown pelicans (Pelecanus occidentalis) and Brandt’s Cormorants (Phalacrocorax penicillatus) in California. J. Zoo Wildl. Med. 1993, 24, 54–62. [Google Scholar]
  24. Beltrán, A.S.; Palafox-Uribe, M.; Grajales-Montiel, J.; Cruz-Villacorta, A.; Ochoa, J.L. Sea bird mortality at Cabo San Lucas, Mexico: Evidence that toxic diatom blooms are spreading. Toxicon 1997, 35, 447–453. [Google Scholar] [CrossRef]
  25. Scholin, C.A.; Gulland, F.; Doucette, G.J.; Benson, S.; Busman, M.; Chavez, F.P.; Cordaro, J.; DeLong, R.; Lefebvre, K.A.; Lipscomb, T.; et al. Mortality of sea lions along the central California coast linked to a toxic diatom bloom. Lett. Nat. 2000, 403, 80–84. [Google Scholar] [CrossRef]
  26. Cook, K.B.; Lacaze, J.P.; MacHairopoulou, M.; Bresnan, E. Investigations into the relationship between domoic acid and copepods in Scottish waters. ICES J. Mar. Sci. 2022, 79, 963–973. [Google Scholar] [CrossRef]
  27. Oller-Ruiz, A.; Campillo, N.; Hernández-Córdoba, M.; Gilabert, J.; Viñas, P. Monitoring lipophilic toxins in seawater using dispersive liquid—Liquid microextraction and liquid chromatography with triple quadrupole mass spectrometry. Toxins 2021, 13, 57. [Google Scholar] [CrossRef] [PubMed]
  28. Sahraoui, I.; Grami, B.; Bates, S.S.; Bouchouicha, D.; Chikhaoui, M.A.; Mabrouk, H.H.; Hlaili, A.S. Response of potentially toxic Pseudo-nitzschia (Bacillariophyceae) populations and domoic acid to environmental conditions in a eutrophied, SW Mediterranean coastal lagoon (Tunisia). Estuar. Coast. Shelf Sci. 2012, 102–103, 95–104. [Google Scholar] [CrossRef]
  29. Barbaro, E.; Zangrando, R.; Rossi, S.; Cairns, W.R.L.; Piazza, R.; Corami, F.; Barbante, C.; Gambaro, A. Domoic acid at trace levels in lagoon waters: Assessment of a method using internal standard quantification. Anal. Bioanal. Chem. 2013, 405, 9113–9123. [Google Scholar] [CrossRef]
  30. Bazzoni, A.M.; Mudadu, A.G.; Esposito, G.; Urru, R.; Ortu, S.; Mara, L.; Uda, M.T.; Arras, I.; Lorenzoni, G.; Sanna, G.; et al. Bacterial and viral investigations combined with determination of phytoplankton and algal biotoxins in mussels and water from a Mediterranean coastal lagoon (Sardinia, Italy). J. Food Prot. 2019, 82, 1501–1511. [Google Scholar] [CrossRef]
  31. Marrouchi, R.; Abdennadher, M.; Feki-Sahnoun, W.; Marzougui, Z.; Hamza, I.S.; Hamza, A.; Kharrat, R. Marine biotoxin profile and Karenia selliformis and Alexandrium minitum occurrence in Boughrara lagoon during the last decade. Ecol. Chem. Eng. Sci. 2022, 29, 501–510. [Google Scholar] [CrossRef]
  32. Bosch-Orea, C.; Sanchís, J.; Farré, M.; Barceló, D. Analysis of lipophilic marine biotoxins by liquid chromatography coupled with high-resolution mass spectrometry in seawater from the Catalan Coast. Anal. Bioanal. Chem. 2017, 409, 5451–5462. [Google Scholar] [CrossRef]
  33. Doucette, G.J.; King, K.L.; Thessen, A.E.; Dortch, Q. The effect of salinity on domoic acid production by the diatom Pseudo-nitzschia multiseries. Nov. Hedwig. 2008, 133, 31–46. [Google Scholar]
  34. Estevez, P.; Gago-Martinez, A. Contribution of mass spectrometry to the advances in risk characterization of marine biotoxins: Towards the characterization of metabolites implied in human intoxications. Toxins 2023, 15, 103. [Google Scholar] [CrossRef] [PubMed]
  35. Moreiras, G.; Leão, J.M.; Gago-Martínez, A. Analysis of cyclic imines in mussels (Mytilus galloprovincialis) from Galicia (NW Spain) by LC-MS/MS. Int. J. Environ. Res. Public Health 2020, 17, 281. [Google Scholar] [CrossRef]
  36. Blanco, J.; Arévalo, F.; Moroño, Á.; Correa, J.; Muñíz, S.; Mariño, C.; Martín, H. Presence of azaspiracids in bivalve molluscs from Northern Spain. Toxicon 2017, 137, 135–143. [Google Scholar] [CrossRef]
  37. Chan, I.O.M.; Tsang, V.W.H.; Chu, K.K.; Leung, S.K.; Lam, M.H.W.; Lau, T.C.; Lam, P.K.S.; Wu, R.S.S. Solid-phase extraction-fluorimetric high performance liquid chromatographic determination of domoic acid in natural seawater mediated by an amorphous titania sorbent. Anal. Chim. Acta 2007, 583, 111–117. [Google Scholar] [CrossRef] [PubMed]
  38. De la Iglesia, P.; Giménez, G.; Diogène, J. Determination of dissolved domoic acid in seawater with reversed-phase extraction disks and rapid resolution liquid chromatography tandem mass spectrometry with head-column trapping. J. Chromatogr. A 2008, 1215, 116–124. [Google Scholar] [CrossRef] [PubMed]
  39. Ao, J.; Gu, J.; Yuan, T.; Li, D.; Ma, Y.; Shen, Z. Applying molecular modelling and experimental studies to develop molecularly imprinted polymer for domoic acid enrichment from both seawater and shellfish. Chemosphere 2018, 199, 98–106. [Google Scholar] [CrossRef] [PubMed]
  40. Liu, M.M.; Liu, S.; Shi, X.; Liu, M.M.; Zhang, Y.D.L.; Wang, Q. Quality evaluation of a herbal prescription through quantification of 40 components by HPLC-ESI-MS/MS. Phytochem. Anal. 2012, 23, 365–372. [Google Scholar] [CrossRef]
  41. D’amore, T.; Lo Magro, S.; Vita, V.; Di Taranto, A. Optimization and validation of a high throughput UHPLC-MS/MS method for determination of the EU regulated lipophilic marine toxins and occurrence in fresh and processed shellfish. Mar. Drugs 2022, 20, 173. [Google Scholar] [CrossRef]
  42. Wang, J.; Chen, J.; He, X.; Hao, S.; Wang, Y.; Zheng, X.; Wang, B. Simple determination of six groups of lipophilic marine algal toxins in seawater by automated on-line solid phase extraction coupled to liquid chromatography-tandem mass spectrometry. Chemosphere 2021, 262, 128374. [Google Scholar] [CrossRef] [PubMed]
  43. ISO 5725-1:1994(E); Accuracy (Trueness and Precision) of Measurements Methods and Results—Part 1: General Principles and Definitions. International Organization for Standarization (ISO): Geneva, Switzerland, 1994.
  44. European Commission. Commission Decision (2002/657/EC) of 12 August 2002 implementing Council Directive 96/23/EC concerning the performance of analytical methods and the interpretation of results. Off. J. Eur. Communities 2004, 180, 84–109. [Google Scholar] [CrossRef]
  45. Lelong, A.; Hégaret, H.; Soudant, P.; Bates, S.S. Pseudo-nitzschia (Bacillariophyceae) species, domoic acid and amnesic shellfish poisoning: Revisiting previous paradigms. Phycologia 2012, 51, 168–216. [Google Scholar] [CrossRef]
  46. Brunson, J.K.; McKinnie, S.M.K.; Chekan, J.R.; McCrow, J.P.; Miles, Z.D.; Bertrand, E.M.; Bielinski, V.A.; Luhavaya, H.; Oborník, M.; Smith, G.J.; et al. Biosynthesis of the neurotoxin domoic acid in a bloom-forming diatom. Science 2018, 361, 1356–1358. [Google Scholar] [CrossRef]
  47. Wells, M.L.; Trick, C.G.; Cochlan, W.P.; Hughes, M.P.; Trainer, V.L. Domoic acid: The synergy of iron, copper, and the toxicity of diatoms. Limnol. Oceanogr. 2005, 50, 1908–1917. [Google Scholar] [CrossRef]
  48. Silver, M.W.; Bargu, S.; Coale, S.L.; Benitez-Nelson, C.R.; Garcia, A.C.; Roberts, K.J.; Sekula-Wood, E.; Bruland, K.W.; Coale, K.H. Toxic diatoms and domoic acid in natural and iron enriched waters of the oceanic Pacific. Proc. Natl. Acad. Sci. USA 2010, 107, 20762–20767. [Google Scholar] [CrossRef] [PubMed]
  49. Alorda-Montiel, I.; Rodríguez-Puig, J.; Alorda-Kleinglass, A.; Diego-Feliu, M.; Rodellas, V.; Garcia-Orellana, J. Historical record and sources of contaminants in sediments from Mar Menor coastal lagoon. In Proceedings of the EGU General Assembly 2021, Online, 19–30 April 2021; p. EGU21-15883. [Google Scholar] [CrossRef]
  50. René, C.; Ruiz, G.; Rodellas, V.; Rodriguez-puig, J.; Santos-echeandía, J. Potentially toxic element input into Mar Menor coastal lagoon (Spain) through submarine groundwater discharge. In Proceedings of the EGU General Assembly 2023, Vienna, Austria, 24–28 April 2023; pp. 24–28. [Google Scholar]
  51. Liefer, J.D.; MacIntyre, H.L.; Novoveská, L.; Smith, W.L.; Dorsey, C.P. Temporal and spatial variability in Pseudo-nitzschia spp. in Alabama coastal waters: A “hot spot” linked to submarine groundwater discharge? Harmful Algae 2009, 8, 706–714. [Google Scholar] [CrossRef]
  52. Bates, S.S.; Garrison, D.L.; Horner, R.A. Bloom Dynamics and Physiology Producing Pseudo-nitzschia Species. In Physiological Ecology of Harmful Algal Blooms; Anderson, D.M., Cembella, A.D., Hallegraeff, G.M., Eds.; Springer: Berlin/Heidelberg, Germany, 1998; pp. 267–292. [Google Scholar]
  53. Kudela, R.; Roberts, A.; Armstrong, M. Laboratory Analyses of Nutrient Stress and Toxin Production in Pseudo-nitzschia spp. from Monterey Bay, California. In Harmful Algal Blooms 2002; Steidinger, K.A., Landsberg, J.H., Tomas, C.R., Vargo, G.A., Eds.; Florida and Wildlife Conservation Commission: Fort Myers, FL, USA; Institute of Oceanography: St. Petersburg, FL, USA; Intergovernmental Oceanographic Commission of UNESCO: Paris, France, 2003. [Google Scholar]
  54. Sun, J.; Hutchins, D.A.; Feng, Y.; Seubert, E.L.; Caron, D.A.; Fu, F.X. Effects of changing pCO2 and phosphate availability on domoic acid production and physiology of the marine harmful bloom diatom Pseudo-nitzschia multiseries. Limnol. Oceanogr. 2011, 56, 829–840. [Google Scholar] [CrossRef]
  55. Lema, K.A.; Latimier, M.; Nézan, É.; Fauchot, J.; Le Gac, M. Inter and intra-specific growth and domoic acid production in relation to nutrient ratios and concentrations in Pseudo-nitzschia: Phosphate an important factor. Harmful Algae 2017, 64, 11–19. [Google Scholar] [CrossRef]
  56. Lundholm, N.; Hansen, P.J.; Kotaki, Y. Effect of pH on growth and domoic acid production by potentially toxic diatoms of the genera Pseudo-nitzschia and Nitzschia. Mar. Ecol. Prog. Ser. 2004, 273, 1–15. [Google Scholar] [CrossRef]
  57. Hasle, G.R.; Syvertsen, E.E. Marine Diatoms. In The Diatoms: Biology and Morphology of the Genera; Round, F.E., Crawford, R.M., Eds.; University Press: Cambridge, UK, 1997; pp. 563–584. [Google Scholar]
  58. Bates, S.S.; Trainer, V.L. The ecology, evolution, and oceanography of harmful algal blooms. Limnol. Oceanogr. 2014, 59, 751–766. [Google Scholar]
  59. Turner, J.T.; Tester, P.A.; Gelsleichter, J. Persistence of domoic acid in seawater under laboratory and field conditions. Nat. Toxins 2000, 8, 271–279. [Google Scholar]
  60. Conesa, H.M.; Jiménez-Cárceles, F.J. The Mar Menor lagoon (SE Spain): A singular natural ecosystem threatened by human activities. Mar. Pollut. Bull. 2007, 54, 839–849. [Google Scholar] [CrossRef]
  61. Bernardeau-Esteller, J.; Sandoval-Gil, J.M.; Belando, M.D.; Ramos-Segura, A.; García-Muñoz, R.; Marín-Guirao, L.; Ruiz, J.M. The role of Cymodocea nodosa and Caulerpa prolifera meadows as nitrogen sinks in temperate coastal lagoons. Diversity 2023, 15, 172. [Google Scholar] [CrossRef]
  62. Burns, J.M.; Ferry, J.L. Adsorption of domoic acid to marine sediments and clays. J. Environ. Monit. 2007, 9, 1373–1377. [Google Scholar] [CrossRef]
  63. Smith, J.; Shultz, D.; Howard, M.D.A.; Robertson, G.; Phonsiri, V.; Renick, V.; Caron, D.A.; Kudela, R.M.; McLaughlin, K. Persistent domoic acid in marine sediments and benthic infauna along the coast of Southern California. Harmful Algae 2021, 108, 102103. [Google Scholar] [CrossRef]
  64. Gilabert, J. Seasonal plankton dynamics in a Mediterranean hypersaline coastal lagoon: The Mar Menor. J. Plankton Res. 2001, 23, 207–218. [Google Scholar] [CrossRef]
  65. López Castejón, F. Caracterización de la Hidrodinámica del Mar Menor y los Flujos de Intercambio con el Mediterráneo mediante Datos in situ y Modelado Numérico. Ph.D. Thesis, Universidad Politécnica de Cartagena, Cartagena, Spain, 2017; 208p. [Google Scholar] [CrossRef]
  66. Pérez-Ruzafa, A.; Campillo, S.; Fernández-Palacios, J.M.; García-Lacunza, A.; García-Oliva, M.; Ibañez, H.; Navarro-Martínez, P.C.; Pérez-Marcos, M.; Pérez-Ruzafa, I.M.; Quispe-Becerra, J.I.; et al. Long-term dynamic in nutrients, chlorophyll a, and water quality parameters in a coastal lagoon during a process of eutrophication for decades, a sudden break and a relatively rapid recovery. Front. Mar. Sci. 2019, 6, 26. [Google Scholar] [CrossRef]
  67. Mercado, J.M.; Cortés, D.; Gómez-Jakobsen, F.; García-Gómez, C.; Ouaissa, S.; Yebra, L.; Ferrera, I.; Valcárcel-Pérez, N.; López, M.; García-Muñoz, R.; et al. Role of small-sized phytoplankton in triggering an ecosystem disruptive algal bloom in a Mediterranean hypersaline coastal lagoon. Mar. Pollut. Bull. 2021, 164, 111989–111999. [Google Scholar] [CrossRef] [PubMed]
  68. Reguera, B.; Alonso, R.; Moreira, Á.; Méndez, S.; Bottein, M.-Y.D. Guide for Designing and Implementing a Plan to Monitor Toxin-Producing Microalgae. In IOC Manuals and Guides, 77; UNESCO: Paris, France; IAEA: Vienna, Austria, 2016; p. 66. [Google Scholar]
  69. Grasshoff, K.; Kremling, K.; Ehrhardt, M. Methods of Seawater Analysis; John Wiley & Sons: Hoboken, NJ, USA, 1999; ISBN 9783527613984. [Google Scholar]
  70. Utermöhl, H. Methods of collecting plankton for various purposes are discussed. SIL Commun. 1958, 9, 1–38. [Google Scholar] [CrossRef]
  71. McNabb, P. Chemistry, Metabolism, and Chemical Analysis of Okadaic Acid Group Toxins. In Seafood and Freshwater Toxins: Pharmacology, Physiology, and Detection; Botana, L.M., Ed.; CRC Press: Boca Raton, FL, USA; Taylor & Francis Group: Abingdon, UK, 2008; p. 209. [Google Scholar]
  72. Gerssen, A.; Mulder, P.P.J.; de Boer, J. Screening of lipophilic marine toxins in shellfish and algae: Development of a library using liquid chromatography coupled to orbitrap mass spectrometry. Anal. Chim. Acta 2011, 685, 176–185. [Google Scholar] [CrossRef]
  73. Pocklington, R.; Milley, J.E.; Bates, S.S.; Bird, C.J.; De Freitas, A.S.W.; Quilliam, M.A. Trace determination of domoic acid in seawater and phytoplankton by high-performance liquid chromatography of the fluorenylmethoxycarbonyl (FMOC) derivative. Int. J. Environ. Anal. Chem. 1990, 38, 351–368. [Google Scholar] [CrossRef]
  74. Regueiro, J.; Martín-Morales, E.; Álvarez, G.; Blanco, J. Sensitive determination of domoic acid in shellfish by on-line coupling of weak anion exchange solid-phase extraction and liquid chromatography-diode array detection-tandem mass spectrometry. Food Chem. 2011, 129, 672–678. [Google Scholar] [CrossRef]
  75. Zhang, W.; Lin, M.; Tong, P.; Lu, Q.; Zhang, L. Ferrite nanospheres-based magnetic solid-phase extraction for determination of domoic acid in seawater samples using high-performance liquid chromatography with tandem mass spectrometry. J. Chromatogr. A 2016, 1443, 54–61. [Google Scholar] [CrossRef] [PubMed]
  76. He, X.; Chen, J.; Wang, J.; Tan, L. Multipoint recognition of domoic acid from seawater by dummy template molecularly imprinted solid-phase extraction coupled with high-performance liquid chromatography. J. Chromatogr. A 2017, 1500, 61–68. [Google Scholar] [CrossRef] [PubMed]
  77. Zhang, Y.; Chen, D.; Hong, Z.; Zhou, S.; Zhao, Y. Polymeric ion exchange material based dispersive micro solid-phase extraction of lipophilic marine toxins in seawater followed by the Q Exactive mass spectrometer analysis using a scheduled high resolution parallel reaction monitoring. Microchem. J. 2018, 138, 526–532. [Google Scholar] [CrossRef]
Figure 1. Elution solvents tested in SPE extraction of domoic acid: (A) study of methanol percentage; (B) study of basic solutions.
Figure 1. Elution solvents tested in SPE extraction of domoic acid: (A) study of methanol percentage; (B) study of basic solutions.
Toxins 15 00526 g001
Figure 2. Elution solvents in SPE extraction of lipophilic toxins and yessotoxins tested: (A,B) study of percentage of methanol; (C,D) study of basic solutions.
Figure 2. Elution solvents in SPE extraction of lipophilic toxins and yessotoxins tested: (A,B) study of percentage of methanol; (C,D) study of basic solutions.
Toxins 15 00526 g002
Figure 3. Time series overlap of Pseudo-nitzschia density (blue diamonds), domoic acid concentration (pale green bars), silicates to phosphates (Si:P) ratio (red dots and line), precipitation (grey bars from above), and salinity (dark cyan dosts and line) in the Mar Menor lagoon. The water DA concentrations can be explained by a trade-off of bottom-up and top-down factors. Bottom-up factors are characterized by the Si:P ratio as a proxy for nutrients limitation; top-down factors are mainly dilution factor due to precipitation (as a proxy of rain over the lagoon plus run-off of the catchment area) and salinity as a proxy for the exchange of water with the Mediterranean Sea, mainly through the Encañizadas inlet. These factors, together with advection processes, explain the poor correlation between Pseudo-nitzschia density and DA concentration.
Figure 3. Time series overlap of Pseudo-nitzschia density (blue diamonds), domoic acid concentration (pale green bars), silicates to phosphates (Si:P) ratio (red dots and line), precipitation (grey bars from above), and salinity (dark cyan dosts and line) in the Mar Menor lagoon. The water DA concentrations can be explained by a trade-off of bottom-up and top-down factors. Bottom-up factors are characterized by the Si:P ratio as a proxy for nutrients limitation; top-down factors are mainly dilution factor due to precipitation (as a proxy of rain over the lagoon plus run-off of the catchment area) and salinity as a proxy for the exchange of water with the Mediterranean Sea, mainly through the Encañizadas inlet. These factors, together with advection processes, explain the poor correlation between Pseudo-nitzschia density and DA concentration.
Toxins 15 00526 g003
Figure 4. Locations of the Mar Menor coastal lagoon, sampling station, and meteorological station with detail of the three inlets connecting the Mar Menor coastal lagoon with the Mediterranean Sea. The cross-sectional area of the Las Encañizadas inlet is 40% of the total cross-sectional area of the lagoon, which has a heavy influence on the water exchange at the sampling point.
Figure 4. Locations of the Mar Menor coastal lagoon, sampling station, and meteorological station with detail of the three inlets connecting the Mar Menor coastal lagoon with the Mediterranean Sea. The cross-sectional area of the Las Encañizadas inlet is 40% of the total cross-sectional area of the lagoon, which has a heavy influence on the water exchange at the sampling point.
Toxins 15 00526 g004
Table 1. LC-MS/MS-QqQ parameters.
Table 1. LC-MS/MS-QqQ parameters.
CompoundMRM Transition (m/z)CE (V)Fragmentor Voltage (V)Retention Time (min)
DA312 → 266 Q101203.4
312 → 248 q10
312 → 161 q20
YTX1141.5 → 1061.3 Q351354.7
1141.5 → 925.5 q60
hYTX1155.4 → 1075.5 Q351354.7
GYM508.3 → 490.6 Q402006.9
13-desM692 → 164 Q502007.7
692 → 444 q20
13,19-didesM678 → 164 Q502007.3
678 → 430 q40
SPX20G706 → 688 Q301908.1
706 → 670 q35
OA827 → 723 Q5519013.8
827 → 809 q45
DTX2827 → 723 Q5519014.7
827 → 809 q45
AZA4844 → 826 Q3019014.7
844 → 808 q45
PTX2881.5 → 837.5 Q6023014.8
DTX1841 → 823 Q4519017.5
841 → 737 q55
AZA3828 → 810 Q3019018.6
828 → 792 q45
AZA5844 → 826 Q3018019.0
844 → 808 q40
AZA1842 → 824 Q3019019.4
842 → 806 q45
AZA2856 → 838 Q3019019.9
856 → 820 q45
Q = quantifier, q = qualifier; Described in [27].
Table 2. Analytical characteristics for SPE-LC-QqQ-MS/MS procedure.
Table 2. Analytical characteristics for SPE-LC-QqQ-MS/MS procedure.
CompoundAqueous Standard Slope
(L ng−1)
Natural Seawater Slope
(L ng−1)
Synthetic Seawater Slope
(L ng−1)
Lineal Range Studied
(ng L−1)
LOD
(ng L−1)
LOQ
(ng L−1)
RSD (%)
DA5.416.295.7445–90014.844.66.19
GYM25.420.222.816–10035.01052.31
13-desM46944447920–2000.060.1710.2
13,19-didesM56439755320–2000.050.1610.4
SPX21520025620–2000.150.448.09
OA0.640.690.75200–200080.62429.12
PTX24.375.134.1415–1002.296.8610.8
AZA45756148073–500.722.166.01
DTX269.566.871.426–1338.6325.96.34
AZA584.51581413–501.043.129.19
DTX14.916.185.2726–1338.9526.913.6
hYTX2.312.502.22100–110091.227412.5
YTX2.152.472.07100–110051.71556.34
AZA32352473533–501.715.1314.9
AZA14275957223–500.150.450.44
AZA23154084813–500.200.613.36
LOD: limit of detection; LOQ: limit of quantification; RSD: average relative standard deviation.
Table 3. SPE procedure recovery study.
Table 3. SPE procedure recovery study.
CompoundLevel Concentration
(ng L−1)
Recovery (%)
DA5083
10075
GYM500118
1000115
13-desM2080
5079
13,19-didesM2082
5074
SPX5091
20080
OA20090
400107
PTX250119
20090
AZA45112
1088
DTX225106
50107
AZA55110
1090
DTX150120
200122
hYTX5081
20077
YTX50108
20078
AZA3587
10105
AZA1599
10110
AZA25114
10120
Table 4. Environmental parameters at the sampling station in the Mar Menor.
Table 4. Environmental parameters at the sampling station in the Mar Menor.
MaxMinMeanStd
Temperature (°C)319.920.866.48
Salinity (PSU)45.6139.7543.061.37
Nitrate [NO3] (μM)3.830.0090.430.75
Nitrite [NO2] (μM)0.4770.010.10.12
Ammonia [NH4+] (μM)10.20.042.382.31
Phosphates [PO4] (μM)0.4340.0020.140.12
Silicates [Si(OH)4] (μM)33.60.0310.588.77
DIN11.0010.072.922.55
DIN:P1872.50.38992.72267.18
DIN:Si108.50.0055.8519.42
Si:P16900.084221.29326.99
Pseudo-nitzschia (cel/L)1.20 × 10701.30 × 1063.00 × 106
DA (ng/L)224.0614.8480.0956.78
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Oller-Ruiz, A.; Alcaraz-Oliver, N.; Férez, G.; Gilabert, J. Measuring Marine Biotoxins in a Hypersaline Coastal Lagoon. Toxins 2023, 15, 526. https://doi.org/10.3390/toxins15090526

AMA Style

Oller-Ruiz A, Alcaraz-Oliver N, Férez G, Gilabert J. Measuring Marine Biotoxins in a Hypersaline Coastal Lagoon. Toxins. 2023; 15(9):526. https://doi.org/10.3390/toxins15090526

Chicago/Turabian Style

Oller-Ruiz, Ainhoa, Nuria Alcaraz-Oliver, Gema Férez, and Javier Gilabert. 2023. "Measuring Marine Biotoxins in a Hypersaline Coastal Lagoon" Toxins 15, no. 9: 526. https://doi.org/10.3390/toxins15090526

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop