Next Article in Journal
Structure-Activity Relationships of Holothuroid’s Triterpene Glycosides and Some In Silico Insights Obtained by Molecular Dynamics Study on the Mechanisms of Their Membranolytic Action
Next Article in Special Issue
Quantitation Overcoming Matrix Effects of Lipophilic Toxins in Mytilus galloprovincialis by Liquid Chromatography-Full Scan High Resolution Mass Spectrometry Analysis (LC-HR-MS)
Previous Article in Journal
Bacillus spp. Inhibit Edwardsiella tarda Quorum-Sensing and Fish Infection
Previous Article in Special Issue
Toxicity Screening of a Gambierdiscus australes Strain from the Western Mediterranean Sea and Identification of a Novel Maitotoxin Analogue
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development of a Fast Liquid Chromatography Coupled to Mass Spectrometry Method (LC-MS/MS) to Determine Fourteen Lipophilic Shellfish Toxins Based on Fused–Core Technology: In-House Validation

Centro de Investigacións Mariñas (CIMA), Pedras de Corón s/n, 36620 Vilanova de Arousa, Spain
*
Author to whom correspondence should be addressed.
Mar. Drugs 2021, 19(11), 603; https://doi.org/10.3390/md19110603
Submission received: 22 September 2021 / Revised: 20 October 2021 / Accepted: 22 October 2021 / Published: 24 October 2021
(This article belongs to the Special Issue Novel Methods for Marine Toxins Detection and Quantification)

Abstract

:
Prevalence and incidence of the marine toxins (paralytic, amnesic, and lipophilic toxins) including the so-called emerging toxins (these are, gymnodimines, pinnatoxins, or spirolides among others) have increased in recent years all over the world. Climate change, which is affecting the distribution of their producing phytoplankton species, is probably one of the main causes. Early detection of the toxins present in a particular area, and linking the toxins to their causative phytoplankton species are key tools to minimize the risk they pose for human consumers. The development of both types of studies requires fast and highly sensitive analytical methods. In the present work, we have developed a highly sensitive liquid chromatography-mass spectrometry methodology (LC-MS/MS), using a column with fused-core particle technology, for the determination of fourteen lipophilic toxins in a single run of 3.6 min. The performance of the method was evaluated for specificity, linearity, precision (repeatability and reproducibility) and accuracy by analysing spiked and naturally contaminated samples. The in-house validation was successful, and the limit of detection (LOD) and quantification (LOQ) for all the toxins were far below their regulatory action limits. The method is suitable to be considered in monitoring systems of bivalves for food control.

1. Introduction

Lipophilic toxins are natural metabolites produced by dinoflagellates which can be extracted from bivalve tissues using organic solvents. Structurally, they belong to five different groups (Figure 1): okadaic acid, including okadaic acid (OA), and dinophysistoxins (DTXs); azaspiracids (AZAs); pectenotoxins (PTXs); yessotoxins (YTXs); and cyclic imines (CIs), which include spirolides (SPXs), pinnatoxins (PnTXs), pteriatoxins (PtTXs), and gymnodimines (GYMs).
Currently, some of these lipophilic toxins are regulated in the European Union. According to the current regulation [1,2], the live bivalve mollusks placed on the market for human consumption must not contain marine biotoxins in total quantities (measured in the whole body or any edible part separately) that exceed: OA, DTXs and PTXs together, 160 µg of okadaic acid equivalents per kilogram; YTXs, 3.75 mg of yessotoxin equivalent per kilogram; and AZAs, 160 micrograms of azaspiracid equivalents per kilogram. On the contrary, no limits are established for any compounds belonging to the CI group [3]. However, several clues indicate that it is more than likely that these “emerging” toxins will have to be considered in the not too distant future: Fundamentally, their increased occurrence in shellfish all over the world [4,5,6,7,8,9,10,11,12], probably due to climate change which is affecting the distribution of their producing phytoplankton species [3,13,14].
For effective monitoring programs, efficient, fast, and low-cost methods are crucial. The high sensitivity and selectivity of liquid chromatography-mass spectrometry make this technique the best choice for the simultaneous evaluation of different lipophilic toxins [16]. In the current European Union regulation [17], liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS) has been established as the reference methodology for monitoring the regulated lipophilic toxins. Because of this, the European Reference Laboratory in Marine Biotoxins (EURLMB) published an EU-Harmonized Standard Operating Procedure (SOP) [18] for the determination of lipophilic marine biotoxins in mollusks. This SOP only includes the regulated toxins and proposes six different chromatographic methods being the fast one 9 minutes long. In recent years, numerous LC-MS /MS methods for the detection of the regulated lipophilic toxins together with some cyclic imines have been published [19,20,21,22,23]. However, the run time of the required chromatograms is relatively long, and they do not include, at least, some especially important cyclic imines such as 13 desmethylspirolide C (13desmSPXC), gymnodimine A (GYMA), 13,19 didesmethylspirolide C (13,19didesmSPXC), 20 methylspirolide G (20MethylSPXG), or pinnatoxin G (PnTXG).
For the Galician monitoring program, Intecmar adapted and accredited (Accreditation N° 160/LE 394) an LC-MS/MS method considering the guidelines of the EURLMB SOP [18] and based on the methods described by Gerssen et al. [20] and Regueiro et al. [24]. The method uses a chromatographic column with sub-2 µm particles and allows for identifying and quantifying the regulated lipophilic toxins in a 9-minute run. The development of the core-shell or fused-core silica particles technology for the chromatographic columns allows for efficient chromatographic separations avoiding some of the drawbacks of the sub-2 µm technology [25], as high backpressure and clogging problems. These characteristics are especially interesting for monitoring systems because: (a) working with lower pressures reduces the stress of the chromatographic equipment making it them less prone to failures; and (b) using larger particles reduces the risk of aggregation of compounds in the column head which can interfere the analyses during long chromatographic batches.
This work describes the validation of a fast LC-MS/MS method, using a fused-core technology chromatographic column, for the determination, in a 3.6 minutes run, of fourteen lipophilic toxins (OA, DTX1, DTX2, PTX2, YTX, HomoYTX, AZA1, AZA2, AZA3, 13desmSPXC, PnTXG, GYMA, 13,19didesmSPXC and 20Methyl SPXG) in bivalve mollusks.

2. Results and Discussion

The LC-MS/MS method developed allows for, as far as we know for the first time, the determination of fifteen lipophilic toxins (Figure 2) (fourteen validated) in a single chromatographic run of 3.6 min. All toxins elute in the first 2.3 min of chromatogram. The peaks have good shape and resolution allowing for correct identification and quantification of the compounds. The 45-OH YTX was not included in the validation because no reference material is commercially available. Using 2.6 µm core-shell (present work) instead of sub-2 µm [22,23] particles allows for effectively reducing the analysis time because it is possible to increase the flow rate maintaining the backpressure, with only marginal losses in efficiency and resolution. Therefore, the analysis time per sample of the developed method is substantially reduced in comparison to previously proposed method [18,21,22,23,24,26,27].

2.1. LOQ/LOD

The limits of detection (LOD) ranged from 0.035 µg kg−1 to 24 µg kg−1 and the limits of quantification (LOQ) ranged from 0.12 µg kg−1 to 74 µg kg−1 (Table 1). The lowest and the highest values were for 13,19 didesmSPXC and HomoYTX, respectively.
Even injecting only 1 µL of sample, LOD and LOQ values are notably below the concentrations reported with other equipment/method, probably because a highly sensitive mass spectrometer was used. In this sense, there is a wide variety of LOQ and LOD data published for lipophilic toxin analyses using Agilent [22,23,26], Thermo [24,27] or Waters [20,21] mass spectrometers among others. Almost all of these papers reported LOQs and LODs values higher than those in the present work, except O‘Neill et al. [28] and Yang et al. [27] with similar and slightly lower values, respectively. However, none of them combines the 14 toxins in the same run and in such a short time
The LOQs values obtained were checked experimentally by fortifying the methanolic extracts of three blank samples (mussel, clam, and oyster) at concentrations close to the theoretical LOQ of each toxin (see selectivity/specificity section).

2.2. Selectivity/Specificity

The analyses to evaluate the recovery of fortified samples at levels around LOQ have been performed under repeatability conditions (n = 3, one day). Recoveries obtained for the three matrices were, in general, good (Table 2). On average for the analysed compounds, they ranged from 84.8 to 111.1%, except for 13,19 didesmSPXC (79.9 %), PnTXG (64.4%), and YTX (123.1 %). The recovery for YTX, was over 100 % for mussels (135.5%) and oysters (137.2%) but not for clams (96.7%). For PnTXG, it was below 100% for mussels (56.9 % of recovery) and oysters (57.7%). The high (for YTXs) and low (for PnTXG) recoveries for mussels and oysters, were the reason why the RSDr (species averaged) of these two toxins were close to 20% instead 15.3% as in the other studied toxins (Table 2). The high value obtained for 13desmSPXC in oysters (142.9%) was due to the presence of traces of this toxin in the blank sample, so it was not included in the computations of mean, SD and RSDr. The low recovery of GYMA in the oyster sample (42.4%) should attributed to its concentration after spiking, that was below its LOQ. This value has not been included in the mean, SD and RSDr calculations.
Since it was the extract that were spiked with toxins, the cause of the differences in recovery should be matrix effect. Matrix effect (ion suppression or ion enhancement) in LC-MS/MS analysis is well known, and can lead to under or overestimate toxin concentration. It is heavily dependent on the analytical system used [29] and, although different strategies have been used to minimize it [24,30] (e.g., SPE clean-up [31]), it is, still present in many published methods [22,23,26].

2.3. Linearity

The method has a good linearity. A straight line fitted well the calibration curve of each toxin, with r2 always above 0.989 (Table 3). Relative slope values were always between 80 and 120%. The random distribution of the residuals around the calibration line confirmed the linearity of these intervals.

2.4. Precision: Repeatability and Reproducibility

The obtained RSDr and RSDR values ranged from 3.1 to 21.9% and 3.0 to 17.4% respectively (Table 4) (except for Level 1 of GYMA and 13desmSPXC). These values comply with IUPAC, FDA and SANCO guidelines (15–20%) confirming the good precision of the method. The high standard relative deviations obtained in Level 1 for 20 MethylSPXG (17.8 and 16.5% for RSDr and RSDR respectively) and GYMA (36 and 30.9% for RSDr and RSDR) are not indicative since they were spiked to concentrations below their LOQ. RSDr and RSDR measured for 13desmSPXC in oysters at Level 1 were abnormally high due to the presence of traces of this toxin in the sample used as blank (it was not possible to obtain oyster completely free of this toxin), these data have not been included in the general computation. All other Level 1 data of this 13desmSPXC, as well as those of PnTXG, were also high (13desmSPXC: 21.2 and 30.6%, PnTXG: 21.9 and 14.8%, for RSDr and RSDR respectively). RSDr and RSDR obtained for OA, PTX2, AZA1, YTX were better than others previously published, for example, by Villar-González et al. [32].
Available naturally contaminated samples used for the precision study only contained OA and derivatives, so the analyses were carried out after performing alkaline hydrolysis (see material and methods section). The absence of natural samples contaminated with toxins other than OA, as it happens in most occasions in Galicia [33], made it impossible to carry out precision studies for components other than that compound. The precision, therefore, was evaluated using the estimations of total OA (measured as free OA after alkaline hydrolysis). For all species, the RSDr (repeatability conditions) were below 11.3%. The maximum RSDR was 14.5% in a mussel sample. Although all RSDR were acceptable, the values obtained in mussels were be slightly higher than in other species (Table 5). Additional information about precision in naturally contaminated samples is provided in Table S1.

2.5. Accuracy

Freeze-dried mussel (Mytilus edulis) tissue certified reference material (CRM -FDMT) has been analyzed seven times on two different days. Except for YTX, all the estimated values were between 92.4 and 132% of the certified concentration (recoveries) (Table 6). Large deviations have been observed for YTX with values of 154.1% the certified concentrations.
No commercially available certified reference materials have been found in matrices other than mussel. Therefore, to complete the studies, blank bivalve matrices (fortified at four or five concentration levels with certified solutions of lipophilic toxins were prepared (see material and methods section, point 3.6.4. for levels)). For each fortified level, recoveries have been calculated using the mean values obtained from three replicates analyzed in three different days in the three chosen blank matrices (n = 27) (Table 7). In the case of the regulated toxins, recoveries were good, always above 84% at the four evaluated levels. For the emerging toxins evaluated (13desmSPXC, 13,19didesmSPXC, 20MethylSPXG, GYMA and PnTXG) recoveries were also acceptable, ranging from 71.1 to 100.8%, except for PnTXG and 13,19didesmSPXC which were 62% and 67.1% respectively (Table 7). As already explained, 13desmSPXC data from oysters at Level 1 have not been considered, since the value was abnormally high, due to the presence of toxin traces in the blank sample. The relatively low recoveries obtained for PnTXG (62–77.1%) and 13,19didesmSPXC (67.1–73.7%) could be attributed to the different matrices used, being the recoveries obtained in clams higher than those in oyster and mussel.

3. Materials and Methods

3.1. Chemicals

Acetonitrile (LC-MS grade) and methanol (HPLC grade quality) were purchased from Scharlab (Spain) and VWR (Spain), respectively. Ultrapure water was obtained from a Milli-Q gradient system fed with an Elix Advantage-10 (Millipore Ibérica, Spain). Ammonium hydroxide (NH4OH, 25%) and sodium hydroxide (NaOH > 99%) were obtained from Merck (Barcelona, Spain), and hydrochloric acid (HCl, 37%) from Panreac (Barcelona, Spain).

3.2. Reference Materials

Certified reference standards (CRMs) for DTX1, YTX, HomoYTX, AZA1, AZA2 and AZA3, and quality control standards (QCSs) for 13desmSPXC, 13,19didesmSPXC and 20MethylSPXG were obtained from CIFGA, S.A. (Lugo, Spain). OA, DTX2, PTX2, GYMA and PnTXG certified solutions and the certified reference material CRM-FDMT1 (freeze-dried mussel tissue) containing OA, DTX2, DTX1, PTX2, 13desmSPXC, YTX, AZA1, AZA2 and AZA3, were acquired from the Institute for Marine Biosciences, National Research Council (NRC), Halifax, Nova Scotia, NS, Canada.

3.3. Shellfish Samples

Method validation has been performed with whole flesh of mussels (Mytilus galloprovincialis), clams (Ruditapes philippinarum) and oysters (Ostrea edulis). All of them were collected from shellfish harvesting areas of Galicia and fortified with certified reference solutions. Furthermore, naturally contaminated samples of mussels (M. galloprovincialis), clams (R. philippinarum), cockles (Cerastoderma edule) and razor clams (Ensis siliqua) from the Galician Rías were also analyzed.

3.4. Sample Extraction and Hydrolysis

Approximately 100 g of soft tissues of each tested bivalve (mussels, clams and oysters), not containing toxins (blank), were homogenized. For OA, DTXs, PTX2 and AZAs, aliquots of these blank samples were spiked with certified reference solutions to concentrations equivalent to LOQ, 0.5, 1, and 1.5 times the regulated levels (see point 3.6.4. for levels). For YTX and HomoYTX the concentrations were lower because their very high regulation limit makes it impossible to attain the required concentration without applying some concentration step to the currently available reference material (see point 3.6.4. for levels). Emerging toxins (13desmSPXC, GYMA, 13,19didesmSPXC, 20MethylSPXG and PnTXG) do not have regulated limits in the EU legislation so, fortification levels have been arbitrary established (see point 3.6.4. for levels). The freeze-dried reference material CRM-FDMT has been reconstituted following the procedure recommended by the manufacturer (NRC-CNRC). Briefly, CRM-FDMT (0.35 g) was reconstituted in a 50 mL centrifuge tube by adding 1.65 mL deionized water, followed by vortex mixing for 30 s and sonication for 1 min in an ultrasonic bath.
Raw, fortified samples and the CRM-FDMT material were extracted according to the standardized operating procedure of the EU-RL for the determination of marine lipophilic biotoxins in mollusks [18]. Briefly, for the analysis of free forms of the OA group (OA, DTXs) and most other lipophilic toxins (PTX2, 13desmSPXC, YTX, HomoYTX, AZA1, AZA2, AZA3, GYMA, 13,19didesmSPXC, 20MethylSPXG, PnTXG) 2-gram aliquots of homogenized tissues were extracted with 9 mL of MeOH 100% (twice) and centrifuged at 2000× g for 10 min. Both supernatants were combined and the final volume was adjusted to 20 mL. An aliquot was filtered through a 0.22 µm PVDF syringe filter, diluted 5/8 with methanol and analyzed by LC-MS/MS. Total OA concentration (main toxin plus its hydrolysable derivatives) was determined after alkaline hydrolysis. For it, 625 µL of 2.5M NaOH were added to a 5 mL aliquot of the methanolic extract, vortexed for 30 s and heated at 76 °C for 40 min; time elapsed, the hydrolysate was cooled to room temperature, weighed to check for solvent losses and neutralized by adding the same volume of 2.5M HCl. As for the analysis of the free forms, an aliquot was filtered through a 0.22 µm syringe filter and diluted with MeOH.

3.5. LC–MS/MS Procedure

The analyses have been carried out on an Exion LC AD™ System (SCIEX, Framingham, MA, USA) coupled to a Qtrap 6500+ mass spectrometer (SCIEX) through an IonDrive Turbo V interface in electrospray mode. The toxins were separated in a Phenomenex Kinetex EVO C18 “core-shell” column 50 mm (length) × 2.1 mm (id), 2.6 µm (particle size). Mobile Phase A was water and B MeCN 90%, both containing 6.7 mM NH4OH (pH 11) [20]. The gradient started with 22% B, was maintained for 0.1 min, followed by a linear increment to reach 95% B at minute 1.8, and maintaining this composition until minute 2.90. The composition was then returned linearly to the initial one in 0.20 min and maintained 0.5 min before the next injection. The flow rate was 1,000 µLmin−1, the injection volume was 1 µL and the column temperature was 40 °C.
The mass spectrometer parameters were optimized by direct infusion using toxin standards, when available, and were set to: ion source Gas 1, 75 (arbitrary units); ion source Gas 2, 75 (arbitrary units); ion spray voltage, 5000 (positive) and −4500 (negative); capillary temperature, 600 (°C); curtain gas, 30; collision Gas, medium. Specific MS/MS fragmentation conditions and collision energies for all the toxins validated are shown in Table 8. The transition with the highest intensity (qn) was used as for quantification, while the second transition (ql, the one with the lowest intensity) was used for confirmatory purposes.

3.6. Scope and Method Validation

The scope of the method includes lipophilic regulated toxins: OA, DTX2, DTX1, PTX2, YTX, HomoYTX, AZA1, AZA2, and AZA3 and five emerging unregulated toxins: 13desmSPXC, GYMA, 13,19diDesMethylSPXC, 20MethylSPXG and PnTXG.
The method has been validated according to the International Organization for Standardization, International Union of Pure and Applied Chemistry (IUPAC), and following the AOAC and Eurachem guidelines. The following parameters have been assessed: LOQ/LOD, specificity, linearity, precision (repeatability and reproducibility) and accuracy.

3.6.1. LOQ/LOD

For the theoretical estimation of LOQ and LOD, a standard curve of toxins has been prepared and expanded with dilutions (2, 4, 6 and 10×) of Level 1 (the one with the lowest concentration). Signal-to-noise ratio (S/N) has been determined directly in each chromatogram for the second transition (the one with the lowest intensity) with Analyst software (SCIEX). Theoretical LOQs for each toxin have been calculated as 10 times the S/N ratio. Theoretical LODs have been estimated by dividing the LOQ by 3.33 to attain an S/N of 3.

3.6.2. Selectivity/Specificity

The selectivity/specificity test was conducted by comparing chromatograms of three blank matrices (M. galloprovincialis, R. philippinarum and O. edulis) with their corresponding methanol extracts fortified with reference solutions at a level close to the theoretical LOQ of each of the toxins evaluated (Table 1). For emerging toxins (13,19didesmSPXC, 20MethylSPXG, GYMA and PnTXG) and to simplify the experimental design, the reference solution used for fortifying samples contained the same toxin concentrations for all of them (3 µg kg−1). For this reason, in the case of GYMA and 20MethylSPXG fortified values are slightly lower than their LOQ and for 13,19didesmSPXC and PnTXG are higher to it. Recoveries (%), mean, standard deviation (SD) and relative standard deviation of repeatability (RSDr %) of the three fortified species have been calculated for each toxin.

3.6.3. Linearity

In this study, linearity has been estimated using a calibration solution curves prepared by dilution CRMs in methanol, and analyzed in triplicate. The concentration levels used are those in Table 9. The calibration line of each toxin and its corresponding coefficient of determination (r2), were calculated by least squares regression weighed with 1/x. Linearity evaluation has been carried out by calculating how many points deviate from the calibration curve according to the next formula: relative slope = (Yi/Xi) × 100 / ∑ (Yi/Xi) / N, where: Xi: toxin concentration (ng mL−1); Yi: area of the peak, and N: number of pairs of Yi / Xi values. The acceptance range was defined as 80–120%.
Residual values have also been calculated, plotted and visually examined to detect suspicious deviations or curvature of the response.

3.6.4. Precision: Repeatability and Reproducibility

Precision has been evaluated using fortified methanolic extracts and naturally contaminated samples. Since precision frequently varies with concentration of analyte, it has been evaluated at different concentrations. The extracts of mussels, clams, and oysters from Galicia, were spiked at the four (or five for 13desmSPXC) concentration levels indicated in the Table 10. Three replicates of each fortified sample were analyzed under repeatability conditions (same day), and intra-laboratory reproducibility conditions (same equipment, same laboratory, three different days).
Contaminated natural samples consisted of mussels (M. galloprovincialis), cockles (C. edule), clams (R. philippinarum) and razor clams (E. siliqua) from Galician Rías, containing OA. They have been analyzed five times under repeatability conditions in one day and five different days for reproducibility. The concentration estimates have been compared to those obtained with other method and LC-MS/MS (Xevo TQ-S triple quadrupole mass, Waters) according to the accredited method (Accreditation N°160//LE 394) used for monitoring by Intecmar and based on Gerssen et al. [20] and Regueiro et al. [24]. After the analyses, mean concentration, SD, RSDr (for repeatability, 1 day), and RSDR (for reproducibility, 5 days) have been calculated.

3.6.5. Accuracy

A freeze-dried mussel tissue certified reference material for multiple marine toxins (CRM-FDMT1) has been used to evaluate the accuracy and traceability of the method. Eight replicates of a methanolic extract on two different days were analyzed in order to estimate the combined effects of the method and the particular laboratory bias. The recovery percentages and/or the bias of the estimates (in relation to the reference value) have been calculated.

4. Conclusions

The performance of the method developed in this work makes it suitable for accurate measurements of the lipophilic toxins regulated by EU legislation (OA, DTX2, DTX1, PTX2, YTX, HomoYTX, AZA1, AZA2, AZA3), and at least five emerging toxins (13desmSPXC, GYMA, 13,19didesmSPXC, 20MethylSPXG, PnTXG). The method allows for the detection and quantification of these fourteen toxins in 3.6 minutes with high specificity, robustness, and with LOQ and LOD low enough to ensure food safety. As far as we know, this method is the fastest one that allows the quantification of so many lipophilic toxins in a single run.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/md19110603/s1, Table S1: Precision results for the analysis of naturally contaminated samples.

Author Contributions

Conceptualization, A.E.R. and J.B.; formal analysis, A.E.R., H.M., C.M.; data curation, A.E.R. and J.B.; writing—original draft preparation, A.E.R.; writing—review and editing, A.E.R. and J.B.; funding acquisition, A.E.R. and J.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by TOXEMER project PGIDIT-CIMA 21/01.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are available on request in the Centro de Investigacións Mariñas.

Acknowledgments

The present work was supported by funds from the Conselleria do Mar (Xunta de Galicia) within the framework of the TOXEMER Project PGIDIT-CIMA 21/01. The authors greatly acknowledge to Intecmar, who supplied both blank and naturally contaminated samples from the official monitoring program of the Xunta de Galicia.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. European Comission. Regulation (EU) No 786/2013 of 16 August 2013 amending Annex III to Regulation (EC) No 853/2004 of the European Parliament and of the Council as regards the permitted limits of yessotoxins in live bivalve molluscs. Off. J. Eur. Union 2014, L 220, 14. [Google Scholar]
  2. European Comission. Regulation (EC) No 853/2004 of the European Parlamient and of the Council of 29 April 2004 laying down specific hygiene rules for on the hygiene of foodstuffs. Off. J. Eur. Union 2004, L 139, 55. [Google Scholar]
  3. EFSA Panel on Contaminants in the Food Chain. Scientific Opinion on marine biotoxins in shellfish—Cyclic imines (spirolides, gymnodimines, pinnatoxins and pteriatoxins). EFSA J. 2010, 8, 1–39. [Google Scholar]
  4. Aasen, J.; Mackinnon, S.L.; Leblanc, P.; Walter, J.A.; Hovgaard, P.; Aune, T.; Quilliam, M.A. Detection and identification of spirolides in Norwegian shellfish and plankton. Chem. Res. Toxicol. 2005, 18, 509–515. [Google Scholar] [CrossRef]
  5. Rundberget, T.; Aasen, J.A.B.; Selwood, A.I.; Miles, C.O. Pinnatoxins and spirolides in Norwegian blue mussels and seawater. Toxicon 2011, 58, 700–711. [Google Scholar] [CrossRef]
  6. Hess, P.; Abadie, E.; Herve, F.; Berteaux, T.; Sechet, V.; Araoz, R.; Molgo, J.; Zakarian, A.; Sibat, M.; Rundberget, T.; et al. Pinnatoxin G is responsible for atypical toxicity in mussels (Mytilus galloprovincialis) and clams (Venerupis decussata) from Ingril, a French Mediterranean lagoon. Toxicon 2013, 75, 16–26. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Rambla-Alegre, M.; Miles, C.O.; de la Iglesia, P.; Fernández-Tejedor, M.; Jacobs, S.; Sioen, I.; Verbeke, W.; Samdal, I.A.; Sandvik, M.; Barbosa, V.; et al. Occurrence of cyclic imines in European commercial seafood and consumers risk assessment. Environ. Res. 2018, 161, 392–398. [Google Scholar] [CrossRef] [Green Version]
  8. Ujevic, I.; Roje-Busatto, R.; Ezgeta-Balic, D. Comparison of amnesic, paralytic and lipophilic toxins profiles in cockle (Acanthocardia tuberculata) and smooth clam (Callista chione) from the central Adriatic Sea (Croatia). Toxicon 2019, 159, 32–37. [Google Scholar] [CrossRef]
  9. García-Altares, M.; Casanova, A.; Bane, V.; Diogene, J.; Furey, A.; de la Iglesia, P. Confirmation of pinnatoxins and spirolides in shellfish and passive samplers from Catalonia (Spain) by liquid chromatography coupled with triple quadrupole and high-resolution hybrid tandem mass spectrometry. Mar. Drugs 2014, 12, 3706–3732. [Google Scholar] [CrossRef]
  10. Villar González, A.; Rodríguez-Velasco, M.L.; Ben-Gigirey, B.; Botana, L.M. First evidence of spirolides in Spanish shellfish. Toxicon 2006, 48, 1068–1074. [Google Scholar] [CrossRef] [PubMed]
  11. Lamas, J.P.; Arévalo, F.; Moroño, A.; Correa, J.; Muñíz, S.; Blanco, J. Detection and spatio-temporal distribution of pinnatoxins in shellfish from the Atlantic and Cantabrian coasts of Spain. Toxins 2019, 11, 340. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Otero, P.; Miguens, N.; Rodríguez, I.; Botana, L.M. LC-MS/MS analysis of the emerging toxin pinnatoxin-G and high levels of esterified OA group toxins in Galician commercial mussels. Toxins 2019, 11, 394. [Google Scholar] [CrossRef] [Green Version]
  13. EFSA Panel on Contaminants in the Food Chain. Scientific Opinion of the Panel on Contaminants in the Food Chain on a request from the European Commission on marine biotoxins in shellfish—domoic acid. EFSA J. 2009, 1181, 1–61. [Google Scholar]
  14. EFSA Panel on Contaminants in the Food Chain. Scientific Opinion on marine biotoxins in shellfish—Emerging toxins: Brevetoxin group. EFSA J. 2010, 8, 1677. [Google Scholar] [CrossRef]
  15. Leyva-Valencia, I.; Hernández-Castro, J.E.; Band-Schmidt, C.J.; Turner, A.D.; O’Neill, A.; Núñez-Vázquez, E.J.; López-Cortés, D.J.; Bustillos-Guzmán, J.J.; Hernández-Sandoval, F.E. Lipophilic Toxins in Wild Bivalves from the Southern Gulf of California, Mexico. Mar. Drugs 2021, 19, 99. [Google Scholar] [CrossRef]
  16. Josić, D.; Rešetar, D.; Peršurić, Ž.; Martinović, T.; Kraljevic Pavelić, S. Detection of microbial toxins by -omics methods: A growing role of proteomics. In Proteomics for Microbial Toxin Identification; Elsevier: Amsterdam, The Netherlands, 2017; pp. 497–505. [Google Scholar]
  17. European Commission. Regulation (EU) 2019/627 of 15 March 2019 laying down uniform practical arrangements for the performance of official controls on products of animal origin intended for human consumption in accordance with Regulation (EU) 2017/625 of the European Parliament and of the Council and amending Commission Regulation (EC) No 2074/2005 as regards official controls. Off. J. Eur. Union 2019, L 131, 51. [Google Scholar]
  18. EURLMB. EU-Harmonised Standard Operating Procedure for Determination of Lipophilic Marine Biotoxins in Molluscs by LC-MS/MS.Version 5. EU-Harmonised Standard Operating Procedure for Determination of Lipophilic Marine Biotoxins in Molluscs by LC-MS/MS.Version 5. Available online: https://www.aesan.gob.es/CRLMB/docs/docs/metodos_analiticos_de_desarrollo/EU-Harmonised-SOP-LIPO-LCMSMS_Version5.pdf (accessed on 30 June 2021).
  19. McNabb, P.; Selwood, A.I.; Holland, P.T. Multiresidue method for determination of algal toxins in shellfish: Single-laboratory validation and interlaboratory study. J. AOAC Int. 2005, 88, 761–772. [Google Scholar] [CrossRef] [Green Version]
  20. Gerssen, A.; Mulder, P.P.J.; McElhinney, M.A.; de Boer, J. Liquid chromatography–tandem mass spectrometry method for the detection of marine lipophilic toxins under alkaline conditions. J. Chromatogr. A 2009, 1216, 1421–1430. [Google Scholar] [CrossRef] [PubMed]
  21. McCarron, P.; Wright, E.; Quilliam, M.A. Liquid chromatography/mass spectrometry of domoic acid and lipophilic shellfish toxins with selected reaction monitoring and optional confirmation by library searching of product ion spectra. J. AOAC Int. 2014, 97, 316–324. [Google Scholar] [CrossRef]
  22. Brana-Magdalena, A.; Leao-Martins, J.M.; Glauner, T.; Gago-Martinez, A. Intralaboratory validation of a fast and sensitive UHPLC/MS/MS method with fast polarity switching for the analysis of lipophilic shellfish toxins. J. AOAC Int. 2014, 97, 285–292. [Google Scholar] [CrossRef] [Green Version]
  23. Moreiras, G.; Leao, J.M.; Gago-Martinez, 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] [Green Version]
  24. Regueiro, J.; Rossignoli, A.E.; Álvarez, G.; Blanco, J. Automated on-line solid-phase extraction coupled to liquid chromatography tandem mass spectrometry for determination of lipophilic marine toxins in shellfish. Food Chem. 2011, 129, 533–540. [Google Scholar] [CrossRef] [PubMed]
  25. Guiochon, G.; Gritti, F. Shell particles, trials, tribulations and triumphs. J. Chromatogr. A 2011, 1218, 1915–1938. [Google Scholar] [CrossRef]
  26. Rodríguez, I.; Alfonso, A.; González-Jartín, J.M.; Vieytes, M.R.; Botana, L.M. A single run UPLC-MS/MS method for detection of all EU-regulated marine toxins. Talanta 2018, 189, 622–628. [Google Scholar] [CrossRef]
  27. Yang, L.; Singh, A.; Lankford, S.K.; Stuart, J.; Rice, D.; Wu, W.H.; Hungerford, J.M. A rapid method for the detection of diarrhetic shellfish toxins and azaspiracid shellfish toxins in Washington State shellfish by liquid chromatography tandem mass spectrometry. J. AOAC Int. 2020, 103, 792–799. [Google Scholar] [CrossRef]
  28. O’Neill, A.; Morrell, N.; Turner, A.D.; Maskrey, B.H. Method performance verification for the combined detection and quantitation of the marine neurotoxins cyclic imines and brevetoxin shellfish metabolites in mussels (Mytilus edulis) and oysters (Crassostrea gigas) by UHPLC-MS/MS. J. Chromatogr. B 2021, 1179. [Google Scholar] [CrossRef]
  29. Fux, E.; Rode, D.; Bire, R.; Hess, P. Approaches to the evaluation of matrix effects in the liquid chromatography-mass spectrometry (LC-MS) analysis of three regulated lipophilic toxin groups in mussel matrix (Mytilus edulis). Food Addit. Contam. 2008, 25, 1024–1032. [Google Scholar] [CrossRef]
  30. Kilcoyne, J.; Fux, E. Strategies for the elimination of matrix effects in the liquid chromatography tandem mass spectrometry analysis of the lipophilic toxins okadaic acid and azaspiracid-1 in molluscan shellfish. J. Chromatogr. A 2010, 1217, 7123–7130. [Google Scholar] [CrossRef]
  31. Gerssen, A.; McElhinney, M.; Mulder, P.; Bire, R.; Hess, P.; de Boer, J. Solid phase extraction for removal of matrix effects in lipophilic marine toxin analysis by liquid chromatography-tandem mass spectrometry. Anal. Bioanal. Chem. 2009, 394, 1213–1226. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Villar-González, A.; Rodríguez-Velasco, M.L.; Gago-Martínez, A. Determination of lipophilic toxins by LC/MS/MS: Single-laboratory validation. J. AOAC Int. 2011, 94, 909–922. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Blanco, J.; Arevalo, F.; Correa, J.; Morono, A. Lipophilic toxins in Galicia (NW Spain) between 2014 and 2017: Incidence on the main molluscan species and analysis of the monitoring efficiency. Toxins 2019, 11, 612. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Structure of the main lipophilic toxins. Source: Leyva-Valencia et al. [15].
Figure 1. Structure of the main lipophilic toxins. Source: Leyva-Valencia et al. [15].
Marinedrugs 19 00603 g001
Figure 2. Chromatogram with two MRM transitions (quantifier and qualifier) obtained from a mixture of reference solutions. Abbreviations: OA (okadaic acid), 45-OH YTX (45 hydroxyyessotoxin), DTX2 (dinophysistoxin 2), HomoYTX (Homoyessotoxin), YTX (yessotoxin), DTX1 (dinophysistoxin 1), AZA3 (azaspiracid 3), AZA1 (azaspiracid 1), AZA2 (azaspiracid 2), GYMA (gymnodimine A), 13,19diDesMetSPXC (13,19 didesmethylspirolide C), 13desmSPXC (13 desmethylspirolide C), PTX2 (pectenotoxin 2), 20MeSPXG (20 methylspirolide G), PnTXG (pinnatoxin G) (A); Zoom of second MRM transitions (qualifier) for OA, DTX2 and AZAs toxins (B).
Figure 2. Chromatogram with two MRM transitions (quantifier and qualifier) obtained from a mixture of reference solutions. Abbreviations: OA (okadaic acid), 45-OH YTX (45 hydroxyyessotoxin), DTX2 (dinophysistoxin 2), HomoYTX (Homoyessotoxin), YTX (yessotoxin), DTX1 (dinophysistoxin 1), AZA3 (azaspiracid 3), AZA1 (azaspiracid 1), AZA2 (azaspiracid 2), GYMA (gymnodimine A), 13,19diDesMetSPXC (13,19 didesmethylspirolide C), 13desmSPXC (13 desmethylspirolide C), PTX2 (pectenotoxin 2), 20MeSPXG (20 methylspirolide G), PnTXG (pinnatoxin G) (A); Zoom of second MRM transitions (qualifier) for OA, DTX2 and AZAs toxins (B).
Marinedrugs 19 00603 g002
Table 1. Limits of detection (LOD) and quantitation (LOQ) obtained for fourteen lipophilic toxins.
Table 1. Limits of detection (LOD) and quantitation (LOQ) obtained for fourteen lipophilic toxins.
ToxinLOD (µg kg−1)LOQ (µg kg−1)
OA310
DTX2929
DTX1929
AZA1929
AZA2929
AZA3929
YTX2272
HomoYTX2474
PTX27.325
13desmSPXC0.20.67
13,19didesmSPXC0.0350.12
20MethylSPXG1.13.8
GYMA1.23.9
PnTXG0.41.3
Table 2. Recovery values in fortified samples at levels around LOQ and performed under repeatability conditions.
Table 2. Recovery values in fortified samples at levels around LOQ and performed under repeatability conditions.
OADTX2DTX1AZA1AZA2AZA3PTX2YTXHomoYTX13desmSPXC13,19didesmSPXC20MethySPXGGYMAPnTXG
Fortified levels (µg kg−1)
30 *30303030303075750.73 *3 **3 **3 *
MatrixRecovery (%)
Mussel89.3100.3114.789.782.387.394.3135.5119.210084.3105.7103.356.9
Clam7180.796.775.791.779.381.396.79285.779.799.3101.778.6
Oyster94105.712210399.795100.3137.2111.1142.9 a75.791.742.4 a57.7
Mean (n = 3)84.895.6111.189.591.287.292.0123.1107.492.9 (n = 2)79.998.9102.5 (n = 2)64.4
SD (n = 3)12.213.213.013.78.77.99.722.914.010.1 (n = 2)4.307.01.1 (n = 2)12.3
RSDr, %14.313.811.715.39.6910.618.613.010.95.47.11.119.1
* Level of fortified slightly above its LOQ for simplicity of analysis. ** Level of fortified slightly below its LOQ for simplicity of analysis. a Data not considered for mean, SD and RSDr calculations.
Table 3. Coefficient of determination (r2) obtained for the fourteen lipophilic toxins.
Table 3. Coefficient of determination (r2) obtained for the fourteen lipophilic toxins.
Toxinr2
OA0.9922
DTX20.9917
DTX10.9918
AZA10.9972
AZA20.9972
AZA30.9966
YTX0.9891
HomoYTX0.9922
PTX20.9969
13desmSPXC0.9969
13,19didesmSPXC0.9921
20MethylSPXG0.9965
GYMA0.9970
PnTXG0.9958
Table 4. Average RSDr (repeatability) and RSDR (reproducibility) of the studied toxin in the three evaluated matrices.
Table 4. Average RSDr (repeatability) and RSDR (reproducibility) of the studied toxin in the three evaluated matrices.
Level 1Level 2Level 3Level 4Level 5
RSDr, %
(n = 9)
RSDR, %
(n = 27)
RSDr, %
(n = 9)
RSDR, %
(n = 27)
RSDr, %(n = 9)RSDR, %
(n = 27)
RSDr, %
(n = 9)
RSDR, %
(n = 27)
RSDr, %
(n = 9)
RSDR, %
(n = 27)
OA8.66.54.47.94.85.83.65.0
DTX28.57.94.36.44.76.75.76.7
DTX110.110.06.69.65.69.73.88.9
AZA111.310.73.14.13.73.23.63.4
AZA28.17.74.74.93.73.04.13.2
AZA39.38.54.54.83.02.74.03.2
YTX16.014.414.914.214.013.712.513.9
HomoYTX13.011.114.013.612.812.712.913.9
PTX28.18.59.37.27.55.15.85.0
13desmSPXC21.2 (n = 6) *30.6 (n = 18) *18.317.46.47.46.96.78.98.9
13,19didesmSPXC15.214.84.910.25.46.47.68.3
20MethylSPXG17.8 **16.5 **4.14.84.76.06.17.5
GYMA36.0 **30.9 **7.47.25.95.47.07.3
PnTXG21.914.84.95.16.37.07.58.7
* Data from oyster were not considered for RSD calculations. ** Level of fortified slightly below its LOQ for simplicity of analysis.
Table 5. Precision results for the analysis of naturally contaminated samples.
Table 5. Precision results for the analysis of naturally contaminated samples.
SampleMatrixToxinRSDr % (n = 5)RSDR % (n = 10)
508MusselTotal OA10.111.7
509MusselTotal OA8.614.5
510MusselTotal OA7.913.0
512MusselTotal OA9.710.3
513MusselTotal OA8.912.5
1532CockleTotal OA9.612.6 *
1534CockleTotal OA8.78.1 *
1540Razor clamTotal OA11.210.0 *
1541CockleTotal OA11.38.7 *
1547ClamTotal OA8.06.9 *
* n = 12
Table 6. Recovery and reproducibility for CRM-FDMT material.
Table 6. Recovery and reproducibility for CRM-FDMT material.
Type of ToxinMean (µg kg−1)SD (µg kg−1)RSDR, % (n = 8)Certified ValueRecovery %
OA306.722.87.4278.3110.2
DTX2795.164.68.1624.8129.1
DTX1155.711.67.5119.0132.0
AZA1809.170.58.7717.5112.3
AZA2194.312.96.6197.898.1
AZA3186.912.56.7168110.6
YTX670.0176.526.3435.8154.1
PTX2107.612.711.8115.592.4
13desmSPXC509.035.46.9472.5108.0
Table 7. Toxin recoveries for different fortification levels in the three matrices studied under reproducibility conditions.
Table 7. Toxin recoveries for different fortification levels in the three matrices studied under reproducibility conditions.
Recovery, % (Mean ± SD), n = 27
Level 1Level 2Level 3Level 4Level 5
OA84.8 ± 1.685.6 ± 5.491.7 ± 8.593.8 ± 14.9
DTX290.4 ± 2.192.5 ± 4.799.6 ± 10.6101.0 ± 21.6
DTX1100.7 ± 3.0105.7 ± 8.2110.7 ± 17.2112.9 ± 32.3
AZA189.9 ± 2.991.2 ± 3.096.0 ± 4.998.7 ± 10.7
AZA290.8 ± 2.191.8 ± 3.497.0 ± 4.699.3 ± 10.0
AZA388.2 ± 2.391.0 ± 3.597.4 ± 4.2099.8 ± 10.3
YTX107.9 ± 11.6113.0 ± 32.2121.1 ± 66.6120.5 ± 134.2
HomoYTX96.2 ± 8.0100.7 ± 27.3107.8 ± 54.7107.4 ± 119.0
PTX285.7 ± 2.184.9 ± 4.989.1 ± 7.291. 5 ± 14.6
13desmSPXC100.8 ± 0.2 (n = 18) *71.1 ± 1.076.5 ± 1.383.4 ± 2.587.5 ± 6.9
13,19didesmSPXC72.8 ± 0.371.8 ± 2.267.1 ± 3.273.7 ± 9.2
20MethylSPXG92.8 ** ± 0.598.9 ± 1.492.3 ± 4.1101.7 ± 11.5
GYMA83.5 ** ± 0.897.2 ± 2.190.6 ± 3.799.3 ± 10.8
PnTXG62.0 ± 0.375.8 ± 1.270.3 ± 3.777.1 ± 10.0
* Data from oyster have not been considered for recovery calculations. ** Level of fortified slightly below its LOQ for simplicity of analysis.
Table 8. MS/MS fragmentation conditions for lipophilic toxin determination. ESI = electrospray ionization mode, Q1 = m/z ratio in the first quadrupole, Q3 = m/z ratio in the third quadrupole, DEP(v) = declustering potential, EP(v) = entrance potential, CE(v) = collision energy, and CXP(v) = collision cell exit potential, qn = ion pair for quantitation and ql = ion pair for qualifier purpose.
Table 8. MS/MS fragmentation conditions for lipophilic toxin determination. ESI = electrospray ionization mode, Q1 = m/z ratio in the first quadrupole, Q3 = m/z ratio in the third quadrupole, DEP(v) = declustering potential, EP(v) = entrance potential, CE(v) = collision energy, and CXP(v) = collision cell exit potential, qn = ion pair for quantitation and ql = ion pair for qualifier purpose.
ToxinESIQ1Q3DEP (v)EP (v)CE (v)CXP(v)
OA_DTX2 (qn)NEG803.52255.15−80−15−62−11
OA_DTX2 (ql)NEG803.52563.40−80−15−60−11
YTX (qn)NEG570.43467.40−80−15−42−11
YTX (ql)NEG570.43396.40−80−15−42−11
HomoYTX (qn)NEG577.40474.40−80−15−42−11
HomoYTX (ql)NEG577.40403.40−80−15−42−11
DTX1 (qn)NEG817.50255.15−80−15−60−11
DTX1 (ql)NEG817.50563.45−80−15−52−11
AZA3 (qn)POS828.46810.580153010
AZA3 (ql)POS828.46658.480154310
AZA1 (qn)POS842.46824.580153010
AZA1 (ql)POS842.46672.480154310
AZA2 (qn)POS856.46838.580153010
AZA2 (ql)POS856.46672.480154310
GYMA (qn)POS508.33490.280155010
GYMA (ql)POS508.33136.0080155010
13,19didesmSPXC (qn)POS678.5164.0080156010
13,19didesmSPXC (ql)POS678.50430.380155010
13desmSPXC (qn)POS692.5164.380156010
13desmSPXC (ql)POS692.50444.380155510
PTX2 (qn)POS876.46823.580153510
PTX2 (ql)POS876.46213.180155010
20MethylSPXG (qn)POS706.5164.380156010
20MethylSPXG (ql)POS706.50348.380155510
PnTXG (qn)POS694.5164.380156010
PnTXG (ql)POS694.5440.380155010
Table 9. Toxin concentrations in ng mL−1 used to evaluate the linearity of the method.
Table 9. Toxin concentrations in ng mL−1 used to evaluate the linearity of the method.
OA, DTXs, AZAs, PTX2YTX, HomoYTX 13desmSPXC13,19didesmSPXC, 20MethylSPXG, GYMA, PnTXG
2.15.20.190.47
6.315.70.60.9
10.626.51.81.9
16.942.13.03.75
21.653.94.77.5
31.779.2615
45.7114.298.8
12.7
Table 10. Toxin concentrations in µg kg−1 used to fortify blank samples.
Table 10. Toxin concentrations in µg kg−1 used to fortify blank samples.
Level 1Level 2Level 3 Level 4Level 5
OA3080160320
DTXs3080160320
AZAs3080160320
PTX23080160320
YTX75200400800
HomoYTX75200400800
13desmSPXC0.78.322.244.489
GYMA33075150
PnTXG33075150
13,19didesmSPXC33075150
20MethylSPXG33075150
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rossignoli, A.E.; Mariño, C.; Martín, H.; Blanco, J. Development of a Fast Liquid Chromatography Coupled to Mass Spectrometry Method (LC-MS/MS) to Determine Fourteen Lipophilic Shellfish Toxins Based on Fused–Core Technology: In-House Validation. Mar. Drugs 2021, 19, 603. https://doi.org/10.3390/md19110603

AMA Style

Rossignoli AE, Mariño C, Martín H, Blanco J. Development of a Fast Liquid Chromatography Coupled to Mass Spectrometry Method (LC-MS/MS) to Determine Fourteen Lipophilic Shellfish Toxins Based on Fused–Core Technology: In-House Validation. Marine Drugs. 2021; 19(11):603. https://doi.org/10.3390/md19110603

Chicago/Turabian Style

Rossignoli, Araceli E., Carmen Mariño, Helena Martín, and Juan Blanco. 2021. "Development of a Fast Liquid Chromatography Coupled to Mass Spectrometry Method (LC-MS/MS) to Determine Fourteen Lipophilic Shellfish Toxins Based on Fused–Core Technology: In-House Validation" Marine Drugs 19, no. 11: 603. https://doi.org/10.3390/md19110603

APA Style

Rossignoli, A. E., Mariño, C., Martín, H., & Blanco, J. (2021). Development of a Fast Liquid Chromatography Coupled to Mass Spectrometry Method (LC-MS/MS) to Determine Fourteen Lipophilic Shellfish Toxins Based on Fused–Core Technology: In-House Validation. Marine Drugs, 19(11), 603. https://doi.org/10.3390/md19110603

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