Next Article in Journal
Comment on Ultrasound Guidance for Botulinum Neurotoxin Chemodenervation Procedures. Toxins 2017, 10, 18—Quintessential Use of Ultrasound Guidance for Botulinum Toxin Injections—Muscle Innervation Zone Targeting Revisited
Next Article in Special Issue
Dinophysis acuta in Scottish Coastal Waters and Its Influence on Diarrhetic Shellfish Toxin Profiles
Previous Article in Journal
Transitioning from Unilateral to Bilateral Upper Limb Tremor Therapy for Parkinson’s Disease and Essential Tremor Using Botulinum Toxin: Case Series
Previous Article in Special Issue
Accumulation and Biotransformation of Dinophysis Toxins by the Surf Clam Mesodesma donacium
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Anatomical Distribution of Diarrhetic Shellfish Toxins (DSTs) in the Japanese Scallop Patinopecten yessoensis and Individual Variability in Scallops and Mytilus edulis Mussels: Statistical Considerations

1
National Research Institute of Fisheries Science, Japan Fisheries Research and Education Agency, Fukuura 2-12-4, Kanazawa-ku, Yokohama, Kanagawa 236-8648, Japan
2
Aomori Prefectural Industrial Technology Research Center, Fisheries Research Institute, Hiranai, Higashitsugarugun, Aomori 039-3381, Japan
3
Hokkaido Research Organization, Fisheries Research Department, Hakodate Fisheries Research Institute, Benten-cho 20-5, Hakodate, Hokkaido 040-0051, Japan
*
Author to whom correspondence should be addressed.
Submission received: 20 August 2018 / Revised: 21 September 2018 / Accepted: 21 September 2018 / Published: 27 September 2018
(This article belongs to the Special Issue Dinophysis Toxins: Distribution, Fate in Shellfish and Impacts)

Abstract

:
Diarrhetic shellfish toxins (DSTs) are a group of phycotoxins that include okadaic acid (OA)/dinophysistoxin (DTX) analogues. At present, detailed data on the distribution of DST is insufficient, and studies of the appropriate sample sizes are lacking. This study investigated the DST frequency distribution in scallops and mussels by liquid chromatography-tandem mass spectrometry (LC/MS/MS) and a resampling analysis of existing data was carried out. The DST population-interval and the necessary sample size were also estimated. DSTs are localized in the scallop digestive-gland, and the DST concentrations in scallops were water-depth-dependent. DST concentrations in scallops and mussels showed normal distributions, but mussels tended to contain more DSTs than scallops. In the statistical resampling analysis of the acquired data on scallops and mussels, especially that using the bootstrap method, sample size was difficult to estimate when the DST variation was large. Although the DST population-interval could be statistically estimated from the sample standard deviation of three samples, the sample size corresponded to the risk management level, and the use of 13 or more samples was preferable. The statistical methods used here to analyze individual contents and estimate population content-intervals could be applied in various situations and for shellfish toxins other than DSTs.
Key Contribution: This is the first detailed analysis of the distribution of individual concentrations of DSTs in shellfish samples, as well as the first report of a method for analyzing and evaluating the relationship between the individual concentrations and mean population concentrations based on statistical methods.

1. Introduction

Diarrhetic shellfish poisoning (DSP) is a severe gastrointestinal illness caused by the consumption of shellfish contaminated with diarrhetic shellfish toxins (DSTs) [1]. DSTs are a group of phycotoxins that include okadaic acid (OA) and dinophysistoxin (DTX) analogues [2,3]. OA, dinophysistoxin-1 (DTX1), and dinophysistoxin-2 (DTX2) are the most important DSTs because they cause severe diarrhea. These toxins have been shown to be potent protein-phosphatase inhibitors [4], a property that can cause inflammation of the intestinal tract and diarrhea [5], possibly leading to tumor promotion [6]. Okadaic acid analogues are metabolized to the esterified toxin in many bivalve species including Japanese scallops [7], and they are collectively called dinophysistoxin-3 (DTX3). In Japan, screening and quantification of DSTs are carried out on bivalves in accordance with the guidelines based on the official instrumental method [8] in production areas and markets. However, the Japanese guidelines do not provide detailed information on the distribution of DSTs between individual bivalves and have not established an appropriate sample size due to the lack of such data obtained by accurate analytical methods [9].
The Japanese scallop Patinopecten yessoensis (Mizuhopecten yessoensis), a major and important cultured species in Japan [10], has unique characteristics, including the metabolic transformation of lipophilic toxins [11]. In the present study, we analyzed the concentrations of DSTs in individuals of P. yessoensis and the mussel Mytilus edulis, and the validity of the size of sample were examined with statistical resampling analysis of the acquired DST data. Although some of our research has been presented in a previous work [12], more detailed data and novel results are provided in our present study. Furthermore, assuming conditions for investigating cultured scallops in the sea [10], estimation of the DST population-interval (interval of concentration of DST contained in population) were performed. Finally, based on our results, we consider and propose an adequate sample size.

2. Results

2.1. Concentrations and Distribution of DSTs

2.1.1. Anatomical Compartmentalization of DST in Scallops

The compartmentalization of DSTs in scallops collected at Nonai Station, Mutsu Bay, Aomori prefecture was investigated. From 14 to 20 individual scallops, the digestive gland, gonad, mantle, gill, and adductor muscle were separately dissected. The pieces were then grouped together by the body part (Table 1).
The dominant toxin in the scallops was DTX1, the highest concentration of which was found on 30 June, corresponding to about half of the regulation value of 0.16 mg/kg of whole meat (Figure 1).
The proportion of the DTX1 quantity corresponding to each tissue are shown in Figure 2.

2.1.2. DST Analysis of 30 Individual Scallops and Mussels

The concentrations of DTX1 in the digestive glands of 30 scallops or mussels collected at the Nonai Station were quantified for each individual (Figure 3 and Table 2).
The concentration of OA was not described due to the overall low concentrations found in the individual samples. The mean values of DTX1 differed even for the same date for both mussels and scallops, and the DTX1 values of the mussels were higher than those of the scallops (Figure 3). The distributions of the scallops and mussels were close to the normal distributions.

2.1.3. Analysis of DST Concentration in Scallop Samples from Different Water Depths

The DST concentrations of scallop digestive glands collected at different depths at Yakumo Station, in the western part of Funka Bay, Hokkaido prefecture, were investigated. Sampling information about the scallops is shown in Table 3.
The vertical gradient of DTX1 distribution (maximum at 5 m, minimum at 15 m) was reversed over the investigation period (Figure 4). Maximum levels of about half of the regulation value (corresponding to 0.16 mg/kg of whole meat) were found at 15 m on 11 August.
The distributions of scallops on 28 June and 11 August were close to the normal distribution (Figure 5).
Environmental conditions and vertical distribution of Dinophysis at Yakumo Station are represented in Table 4 and Figure 6, respectively.
D. acuminata was assumed to be main causative agent of the DST events. This species reached its maximal density on 28 June and declined on 11 August (Table 4). Other Dinophysis species were predominant from 10 m to the surface, and there was no clear relationship between the distribution of cells and reversal of the vertical distribution of DTX1 in scallops on 11 August (Figure 4 and Figure 6).

2.2. Statistical Analysis

2.2.1. Statistical Resampling Analysis of DSTs in Scallops and Mussels

The number of individuals necessary to correctly reflect the DST contents of collected samples was estimated by resampling analysis. For this purpose, we used a data set collected for 30 individual scallops and 30 individual mussels at Aomori on 2 June (Figure 3). Both values were highly variable. The means of 5–25 samples were calculated with random sampling and allowing or not allowing (bootstrap method) replacement. The sample means were converted to percentages as compared with those of 30 samples (Table 5).
In the resampling analysis without replacement, using ≥8 and ≥17 scallops fell within ±20% and 10% of the means of 30 individuals, respectively, with a probability of 98%. In the case of mussels, means of ≥13 and ≥19 individuals fell within ±30% and 20% of the means of 30 individuals, respectively, and with 98% probability (Table 5, underlined numbers in the left half). In the bootstrap method, only 11 scallops fell within ±20% of the mean value of 30 individuals with a probability of 98% (Table 5, double-underlined numbers in the right half), whereas the number of mussels was not obtained from 25 individuals.

2.2.2. Estimating the Mean Concentration of the Population (Cultured Scallops) When the Individual Concentration of a Sample Is Defined

In the former subsection, the ratio of the sample mean to the 30-individual population DST mean was analyzed using specific individual samples. Here, the mean DST concentration of the population (cultured scallops) was estimated using the sample mean concentration. Estimation of the standard normal distribution population mean with the confidence interval (CI, 95% = 100 (1 − α)%) is represented by the following equation [13]:
  X ¯ 1.96 × Var ( X ¯ ) μ ^ X ¯ + 1.96 × Var ( X ¯ )  
By the central limit theorem, variance Var ( X ¯ ) = σ 2 n [13].
  X ¯ 1.96 × σ 2 n μ ^ X ¯ + 1.96 × σ 2 n  
The unknown-population standard deviation σ can be replaced with the sample standard deviation s calculated from the sample data. Moreover, 1.96 = t0.05(∞) and hence is generalized with t α ( ν ) .
  X ¯ t α ( ν ) × s / n   μ ^ X ¯ + t α ( ν ) × s / n  
In Equation (3), t α ( ν ) × s / n   (CI of the population mean concentration) was estimated when s was 200 at [α = 0.05] using 10 samples. In accordance with n and α, t0.05(9) = 2.2622 from the Student’s t-distribution in Table 6 was assigned in Equation (3).
  t α ( ν ) × s / n =   2.2622   ×   200 / 10   =   143.1  
At α = 0.05 and 0.10, the interval estimations (±OA group ng/g) of μ with sample sizes of 3–20 and s of 100–650 were calculated and are presented in Table 7 and Table 8.
From Table 7 and Table 8, it is possible to estimate the confidence interval of the μ when s is 100–650 for 3–20 samples, when the individual concentration data of the sample are acquired.
From Equation (5), derived from Equation (3), the s value can be calculated to estimate the mean density of μ within interval ±160 ng/g (one-tenth of the digestive gland regulatory limit).
  t α ( ν ) × s / n   160  
  s 160 / t α ( ν ) × n .  
In Equation (5), the s was estimated using t0.05(9) = 2.2622 from Student’s t-distribution table.
  s 160 / t α ( ν ) × n =   160 / 2.2622   ×   3.1622   =   223.7  
At α = 0.05 and 0.10, the estimation of s with the mean density of μ within interval ±160 ng/g with sample sizes of 3–20 were calculated and are presented in Figure 7.
For an actual sample, when the calculated s is less than or equal to the graph value (striped zone of Figure 7), it is considered that the mean concentration of μ can be estimated with an interval of ±160 ng/g and 90% or 95% confidence.

2.2.3. Adequacy of Sample Size Based on the t-Value and Confidence Interval

The values in Table 7 and Table 8 and in Figure 7 are derived by equations with a t-value as a coefficient. n is a natural number, and the statistical degrees of freedom (ν) is a linear function. On the other hand, the t-value is defined by Equation (6) [14,15], and the upper cumulative probability (1/2 α) of the t distribution is derived from Equation (7).
  f ( t ) =   Γ ( ν + 1 2 ) ν π Γ ( ν 2 ) ( 1 + t 2 ν ) ν + 1 2   ( Γ   is   Gamma   function )  
  1 / 2 α = x f ( t ,   υ ) d t  
In Equation (7), although the t-value (x) is determined by the degrees of freedom ν and significance level α, it is inversely correlated with the ν and α, and the t-value (x) gradually approaches infinite degrees of freedom at each significance level. Figure 8 represents the t-value of α = 0.10 or 0.05 with ν = 2–30. The t-value increases as the confidence value becomes greater, without intersection.
Ultimately the line is nearly straight with zero slope passing through the t-value at each infinity ν (t0.10(∞) = 1.6449 or t0.05(∞) = 1.9600 equal to the standard normal distribution at each confidence level). When t = 0.05 or 0.10, points of n = 20, 25, 30 and approximate straight lines are drawn on a t-value graph (Figure 8, red double lines). Although the coefficient of approximate straight lines varies depending on the desired value and confidence, the risk of the obtained estimation value is greatly reduced as the sample size increases and the t-value approaches each linear approximation.
Both approximate straight lines start to diverge from the t-value around ν = 13 (n = 14), and the t-value increases exponentially as ν decreases. From the samples of about 14 or more, the relation between the risk of the estimated value and the size of the sample is assumed to describe a linear function. Thus, it is desirable to use a sample size of at least 14 when estimating the mean concentration of the population. When the sample size is 13 or less, the risk of the estimate value increases exponentially as the sample size decreases.
The graph of f(t, ν = 100) and the standard normal distribution f(t, ν = ∞) are almost identical (Figure 9c). The difference between t0.10(100) and t0.10(∞) or that between t0.05(100) and t0.05(∞) is 0.0153 or 0.024, respectively (Table 6). These levels are within a margin of error that does not matter practically. The ideal number of samples is 14 (ν = 13) or more, but considering the mathematical errors, 13 (ν = 12) or more samples is assumed to be a practical allowable range.

3. Discussion

Some Dinophysis species in Japan produces DTX1 as the dominant toxin [7,16]. In scallops, results showing that DSTs were detected exclusively in the digestive glands and that the adductor muscles were free of DSTs agree with the results of our previous feeding experiment study [11].
In our samples, the ratio of the digestive glands to the whole meat of the mussels was about 15% (data from three groups of 30 mussel samples). Some of the mussels collected on 2 June and 7 July exceeded the regulatory level, but none of the scallop samples exceeded the regulatory level. Even though the scallops and mussels were cultured in the same spot, the amounts of accumulated DSTs and their variabilities tended to be higher in mussels than in scallops. The mussels adhered to each other via the byssus, whereas the scallops were separated from each other in the net. Ecological factors and metabolism may be involved in the magnitude of DST variability in mussels and scallops. Concentration fluctuations and DST variation determined by water depth (Figure 4) seem to reflect the vertical distribution of Dinophysis in the water column and the individual metabolism of scallops. The peak of Dinophysis species cell density was assumed to have occurred in July, but unfortunately our data did not identify this trend.
Given a sufficient sample size, it is common practice to homogenize and prepare samples, and this makes sense from the viewpoint of equalizing the samples. Generally, according to the law of large numbers, when an appropriate sample size is collected, there is no problem in obtaining the mean value of the population even if the samples are combined. However, this theorem does not present the validity of the sample size, and in fact it is a problem that convergence of mean value requires plenty of samples. Hence, there is important question as to what the pooled sample size should be. It should also be noted that valuable information on the toxin concentrations in individual shellfish flesh is lost by homogenization and cannot be used to estimate risk.
In the current DST testing of scallops or mussels, it is not feasible to collect 30 individuals, so the number of individuals necessary to reflect the mean DST content of 30 samples collected were estimated in scallops and mussels by statistical resampling analysis using the actual values. Approximately 8 or more samples were considered adequate for scallops, as the variation in DTX1 levels was less than that of mussels. In the bootstrap analysis of more highly variable mussels, 25 samples were insufficient to fall within ±30% of mean value with a probability of 98%. The sample size required changes according to the desired degree of uncertainty level, and the calculation results are restricted to those specific already-known sample data.
As mentioned above, usually a sample size of 30 is difficult to obtain for scallops or mussels. Therefore, we used statistical methods to estimate the mean value of the general population under more practical conditions. Estimation of the population by statistical processing and estimation of the size of necessary samples are basic methods that are described in some textbooks [13,17]. However, no cases have been applied to the field of shellfish toxins. To estimate the mean value of the DST populations of shellfish samples, the samples must be distributed normally in order to apply the statistical parametric test equation. From the individual data group (Figure 3 and Figure 5), the DSTs of each sample basically have a one-peak distribution in which the histogram is almost symmetrical, and the mean and median values are almost matched. Thus, regardless of the concentration, when the variation was not extremely large, it was considered normally distributed. Moreover, it is necessary to take as random a sample as possible from the fact that the total of the depth distribution graph with 5–15 m takes a more normal distribution (Figure 5d,h).
According to the calculation results of Equations (3) and (5) (Table 7 and Figure 7), if there is individual information about the samples, the risk can be evaluated as a concrete figure for the mean concentration of the population. Even when the sample size is less than 5, the interval estimation against DST can be obtained from the value. This makes it possible to evaluate whether sampling is sufficient or not as well as the risk of estimating population concentration. Attention should be paid to the difference between σ and s in these statistical calculations. The statistical definition of σ is completely different from that of s, and s is required for estimating the value. Although in this report we did not deal directly with other shellfish toxins such as paralytic shellfish toxins (PST), the mean concentration of the population and the risk can be estimated and evaluated according to Equation (3) in the same manner.
As the sample size increases in t-value, the slope of the linear approximation decreases and the value becomes strict (Figure 8). Finally, it is a straight line with zero slope passing through the infinity t-value with the necessary reliability (ex, y = t0.10(∞) = 1.6449 or y = t0.05(∞) = 1.9600). The larger the sample size, the better, but in practice there are many cases where statistical ideals are not satisfied due to various restrictions. In the case of scallops or mussels, calculation of an approximate expression using an impractically large size of samples is irrational and not applicable to real-world conditions. Hence, the linear approximate expression at each reliability (α = 0.10 or 0.05) was calculated using samples of 20, 25, and 30 in this study (Figure 8). Samples of about 14 or more for scallops or mussels were derived as ideal sample sizes as the result of the estimation from Figure 8. Because the graph of f(t, ν = 100) and the standard normal distribution f(t, ν = ∞) are almost identical (Figure 9), 13 or more samples are considered to be a practical preferred range including mathematical errors.
On the other hand, in an investigation or analysis that can sample from 100 to 1000 individuals, another criterion corresponding to such sample size should be applied. A small sample size such as n = 3 and a high risk of α = 0.10 are presented in this report for research purposes. The desired minimum sample size is 13 or more, and it is necessary to carefully consider the risk corresponding to α = 0.10, meaning a rejection rate of 10%.
In conclusion, our study shows that DSTs in scallops and mussels are localized in the digestive gland and the DST concentrations have a normal distribution. Statistical analysis of the normal distribution data enables estimation of DST population-interval and shows that a sample size of 13 or more individuals is desirable. Simple evaluations and calculations using tables in this article can be applied in various situations. Although it is inevitable to combine samples at the time of the actual inspection of shellfish, in research, it is desirable to try to acquire as much important individual information as possible, to obtain more accurate values, and to evaluate the risks. Those results are used as an index for risk assessment and are expected to contribute to risk management in shellfish toxins.

4. Materials and Methods

4.1. Plankton Monitoring

One liter of seawater was sampled from May to August 2016 using Van Dorn bottles at 5 m depth intervals from Yakumo Station (42°16.208′ N, 140°20.568′ E) in Uchiura Bay (Funka Bay), Hokkaido, Japan. Each 1 L sample was concentrated and resuspended to 10 mL by filtration through a 20-μm mesh plankton net sieve and fixed with 1.25% glutaraldehyde. To estimate the densities of Dinophysis, 1 mL of each sample was stained with 0.01% of fluorescent dye (Whitex BB, Sumitomo Chemical, Chuo-ku, Tokyo, Japan) and observed with an inverted epifluorescence microscope (IX71, Olympus, Shinjuku-ku, Tokyo, Japan) under UV light excitation. Vertical profiles of temperature and salinity were obtained from CTD (RINKO-Profiler ASTD102, JFE Advantech, Nishinomiya, Hyogo, Japan) casts. Water transparency was recorded using the Secchi disk (30 cm in diameter, RIGO CO. LTD., Bunkyo-ku, Tokyo, Japan).

4.2. Scallops and Mussels

Scallops and mussels, several individuals grown on each lantern net [10] at the same point, were collected from Nonai Station (40°52′ N, 140°07′ E, Depth = 32 m) in Mutsu Bay, Aomori Prefecture, Japan, in 2014. The Aomori prefecture is located at the northern end of Honshu Island. Other scallops with ear-hanging [10] were harvested near Yakumo Station (42°16.558′ N, 140°20.000′ E) of Uchiura Bay, Hokkaido, Japan, in 2016. Hokkaido is the northernmost prefecture of Japan.

4.3. Extraction of DSTs and Hydrolysis of Esterified DSTs

Each dissected tissue was homogenized with 9 volumes of methanol-distilled water (9:1, v/v), and the homogenates were centrifuged at 3000 rpm for 5 min [18]. Alkaline hydrolysis of the OA group was carried out according to the EU harmonized standard operating procedure for lipophilic marine biotoxins in molluscs by LC-MS/MS, ver. 5 [19]. For hydrolysis, 125 μL of 2.5 M NaOH solution was added to a 1 mL aliquot of a methanolic extract of each sample. The mixture was kept at 80 °C for 30 min and neutralized with 125 μL of 2.5 M HCl. The hydrolyzed samples were analyzed by LC/MS/MS without further purification.

4.4. Standard Toxins

The National Metrology Institute of Japan certified reference material of okadaic acid (OA) and dinophysistoxin-1 (DTX1) [20] were dissolved in HPLC-grade methanol to prepare the calibration standards.

4.5. LC/MS/MS Analysis of DSTs

OA and DTX1 in sample extracts were analyzed and quantified by LC/MS/MS as reported previously [16,18]. Triplicate analyses were carried out for each sample extract. Multiple reaction monitoring (MRM) LC/MS/MS analysis for toxins was carried out using [M − H] as target parent ions in Q1 and particular fragment ions of each toxin in Q3, with a dwell time of 100 ms for each analogue as follows. OA: m/z 803.5 > 255.3; DTX1: m/z 817.5 > 255.3. LoD (limit of detection) of OA and DTX1 < 0.01 mg/kg. The proportion of the DTX1 quantity corresponding to each tissue was calculated by multiplying the concentrations by the total tissue weight.

4.6. Statistical Analyses

Statistical analysis program R [21] with boot (https://CRAN.R-project.org/package=boot) and MASS (https://CRAN.R-project.org/package=MASS) packages were employed for resampling analysis and for graphing the standard normal distribution and t-distribution. The values repeatedly calculated 10,000 times with 5–25 random samplings from 30 scallops or mussels were converted to percentages as compared with the mean value of 30 scallops or mussels. This computation is a kind of bootstrap method that has been modified so as not to allow replacement. This provides an estimate of the mean distribution, and how the mean varies depending on the size of samples is presented. The bootstrap method is also applied in the same manner except to allow replacement. Since the resampling analysis without replacement picks up different random samples and the bootstrap method may pick up random samples, including the same samples, more stringent results can be obtained.
Population is parent population.
Sample size is the number of individuals within a group. Number of samples is the number of groups.
  • N = number of individuals in the population.
  • n = number of individuals in the sample.
The population mean is represented by the Greek letter mu ( μ ), and μ ^ represents the estimator.
The x i is an individual sample value, and the sample mean is represented by x ¯ .
The population variance is denoted by σ 2 .
  σ 2 = 1 N i = 1 N ( x i μ ) 2  
The population standard deviation is denoted by σ .
  σ = σ 2  
The sample variance is denoted by s 2 .
  s 2 = 1 n 1 i = 1 n ( X i X ¯ ) 2  
The sample standard deviation is denoted by s .
  s = s 2 .  
  • Two-tailed significance level (α) = 0.05 or 0.10.
  • Confidence level = 1 − α.
  • Confidence interval (CI) = 100 (1 − α)%.
  • Statistical degrees of freedom = n − 1 = Greek letter nu (ν).
  • t is derived from Student’s t-distribution table using n − 1 and α (Table 6).
According to the central limit theorem, regardless of the distribution of the population, if the sample number is made sufficiently large, the error between the population mean and the sample mean follows a normal distribution.
According to the law of large numbers, if the sample size is made sufficiently large, the sample mean converges to the population mean.

Author Contributions

T.S., H.O., and T.A. interpreted data and supervised the research. I.O., Y.K., and M.K. collected and prepared the samples. R.M., H.U., and R.W. analyzed the toxins. R.M. and T.S. wrote the paper.

Funding

The series of studies was conducted and supported by the Ministry of Agriculture, Forestry, and Fisheries (Regulatory Research Projects for Food Safety).

Acknowledgments

We are grateful to Hidetsugu Yoshida of the Mariculture Fisheries Research Institute, Hokkaido Research Organization, and to Daisuke Achiya of the Yakumo Town Fisheries Cooperative for their helpful assistance in the field samplings and measurement. We would like to thank Lincoln Mackenzie for assistance with English corrections and his valuable comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yasumoto, T.; Oshima, Y.; Yamaguchi, M. Occurrence of a new type of shellfish poisoning in the Tohoku. Bull. Jpn. Soc. Sci. Fish. 1978, 44, 1249–1255. [Google Scholar] [CrossRef]
  2. Yasumoto, T.; Murata, M.; Oshima, Y.; Sano, M.; Matsumoto, G.K.; Clardy, J. Diarrhetic shellfish toxins. Tetrahedron 1985, 41, 1019–1025. [Google Scholar] [CrossRef]
  3. Yasumoto, T.; Murata, M. Marine toxins. Chem. Rev. 1993, 93, 1897–1909. [Google Scholar] [CrossRef]
  4. Bialojan, C.; Takagi, A. Inhibitory effect of a marine-sponge toxin, okadaic acid, on protein phosphatases. Specificity and kinetics. Biochem. J. 1988, 256, 283–290. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Terao, K.; Ito, E.; Yanagi, T.; Yasumoto, T. Histopathological studies on experimental marine toxin poisoning. I. Ultrastructural changes in the small intestine and liver of suckling mice induced by dinophysistoxin-1 and pectenotoxin-1. Toxicon 1986, 24, 1145–1151. [Google Scholar] [CrossRef]
  6. Fujiki, H.; Suganuma, M.; Suguri, H.; Yoshizawa, S.; Takagi, K.; Uda, N.; Wakamatsu, K.; Yamada, K.; Murata, M.; Yasumoto, T.; et al. Diarrhetic shellfish toxin, dinophysistoxin-1, is a potent tumor promoter on mouse skin. Jpn. J. Cancer Res. 1988, 79, 1089–1093. [Google Scholar] [CrossRef] [PubMed]
  7. Suzuki, T.; Ota, H.; Yamasaki, M. Direct evidence of transformation of dinophysistoxin-1 to 7-O-acyl-dinophysistoxin-1 (dinophysistoxin-3) in the scallop Patinopecten yessoensis. Toxicon 1999, 37, 187–198. [Google Scholar] [CrossRef]
  8. Guidelines for Risk Management of Shellfish Toxins in Bivalves. Available online: http://www.maff.go.jp/j/syouan/tikusui/gyokai/g_kenko/busitu/pdf/150306_kaidoku_guide.pdf (accessed on 25 September 2018).
  9. Standard for Liver and Raw Bivalve Molluscs (CODEX STAN 292-2008). Available online: http://www.fao.org/fao-who-codexalimentarius/sh-proxy/es/?lnk=1&url=https%253A%252F%252Fworkspace.fao.org%252Fsites%252Fcodex%252FStandards%252FCODEX%2BSTAN%2B292-2008%252FCXS_292e_2015.pdf (accessed on 25 September 2018).
  10. Sandra, E.S.; Parson, J.G. Scallops: Biology, Ecology, Aquaculture, and Fisheries, 3rd ed.; Elsevier Science: Amsterdam, The Netherlands, 2016; Volume 40, pp. 891–936. ISBN 978-0-444-62710-0. [Google Scholar]
  11. Matsushima, R.; Uchida, H.; Nagai, S.; Watanabe, R.; Kamio, M.; Nagai, H.; Kaneniwa, M.; Suzuki, T. Assimilation, accumulation and metabolism of dinophysistoxins (DTXs) and pectenotoxins (PTXs) in the Japanese scallop Patinopecten yessoensis. Toxins 2015, 7, 5141–5154. [Google Scholar] [CrossRef] [PubMed]
  12. Matsushima, R.; Uchida, H.; Watanabe, R.; Oikawa, H.; Kosaka, Y.; Tanabe, T.; Suzuki, T. Distribution of Diarrhetic Shellfish Toxins in Mussels, Scallops, and Ascidian. Food Saf. 2018, 6, 101–106. [Google Scholar] [CrossRef]
  13. Lemeshow, S.; Hosmer, D.W.; Klar, J.; Lwango, S.K. World Health Organization. Adequacy of Sample Size in Health Studies; John Wiley & Sons Ltd.: West Sussex, UK, 1990; pp. 36–62. ISBN 0471925179. [Google Scholar]
  14. Student. The probable error of a mean. Biometrika 1908, 6, 1–25. [Google Scholar] [CrossRef]
  15. Fisher, R.A. Applications of “tudent’s” distribution. Metron 1925, 5, 90–104. [Google Scholar]
  16. Suzuki, T.; Miyazono, A.; Baba, K.; Sugawara, R.; Kamiyama, T. LC–MS/MS analysis of okadaic acid analogues and other lipophilic toxins in single-cell isolates of several Dinophysis species collected in Hokkaido, Japan. Harmful Algae 2009, 8, 233–238. [Google Scholar] [CrossRef]
  17. Sampling Plans for Aflatoxin Analysis in Peanuts and Corn; Report for FAO Technical Consultation; FAO: Rome, Italy, May 1993.
  18. Suzuki, T.; Jin, T.; Shirota, Y.; Mitsuya, T.; Okumura, Y.; Kamiyama, T. Quantification of lipophilic toxins associated with diarrhetic shellfish poisoning in Japanese bivalves by liquid chromatography-mass spectrometry and comparison with mouse bioassay. Fish. Sci. 2005, 71, 1370–1378. [Google Scholar] [CrossRef]
  19. EU-Harmonised Standard Operating Procedure for Determination of Lipophilic Marine Biotoxins in Molluscs by LCMS/MS Version 5. Available online: http://www.aecosan.msssi.gob.es/AECOSAN/docs/documentos/laboratorios/LNRBM/ARCHIVO2EU-Harmonised-SOP-LIPO-LCMSMS_Version5.pdf (accessed on 25 September 2018).
  20. NMIJ CRM Catalog 2016–2017, National Institute of Advanced Industrial Science and Technology (AIST), National Metrology Institute of Japan (NMIJ). Available online: https://www.nmij.jp/english/service/C/CRM_Catalog_(JE)160901.pdf (accessed on 25 September 2018).
  21. R: A Language and Environment for Statistical Computing. Available online: https://www.R-project.org/ (accessed on 25 September 2018).
Figure 1. Concentrations of DTX1 and OA in the digestive glands of scallops. Fourteen to twenty individuals were combined into each sample set used for analysis. Black bars and white bars represent DTX1 and OA, respectively. The concentrations of toxins in the samples are shown on the vertical axis.
Figure 1. Concentrations of DTX1 and OA in the digestive glands of scallops. Fourteen to twenty individuals were combined into each sample set used for analysis. Black bars and white bars represent DTX1 and OA, respectively. The concentrations of toxins in the samples are shown on the vertical axis.
Toxins 10 00395 g001
Figure 2. Percentage (%) of DTX1 in each scallop tissue.
Figure 2. Percentage (%) of DTX1 in each scallop tissue.
Toxins 10 00395 g002
Figure 3. The 30 individual distributions, means ± σ of DTX1 in scallops and mussels. Panels (ad) show the results for scallops and (eg) show those for mussels.
Figure 3. The 30 individual distributions, means ± σ of DTX1 in scallops and mussels. Panels (ad) show the results for scallops and (eg) show those for mussels.
Toxins 10 00395 g003
Figure 4. Vertical distribution of DTX1 and σ in scallop digestive glands.
Figure 4. Vertical distribution of DTX1 and σ in scallop digestive glands.
Toxins 10 00395 g004
Figure 5. The distributions of scallop at each depth on 28 June and 11 August. The mean concentrations of DTX1 ± σ are shown for (ad) 28 June and (eh) 11 August.
Figure 5. The distributions of scallop at each depth on 28 June and 11 August. The mean concentrations of DTX1 ± σ are shown for (ad) 28 June and (eh) 11 August.
Toxins 10 00395 g005
Figure 6. Vertical distributions of Dinophysis species. (a) D. acuminata; (b) D. fortii and other Dinophysis species.
Figure 6. Vertical distributions of Dinophysis species. (a) D. acuminata; (b) D. fortii and other Dinophysis species.
Toxins 10 00395 g006
Figure 7. Estimation of s with the mean density of μ within interval ±160 ng/g at (a) 95% or (b) 90% confidence. The X axis n represents a sample size of 3–20.
Figure 7. Estimation of s with the mean density of μ within interval ±160 ng/g at (a) 95% or (b) 90% confidence. The X axis n represents a sample size of 3–20.
Toxins 10 00395 g007
Figure 8. The t-value of α = 0.05 or 0.10 with ν = 2–30. (a) Triangles are t0.05(2–30); (b) diamonds are t0.10(2–30). The dotted lines show (a) y = 1.9600 and (b) y = 1.6449. The red circles represent the t-values of n = (20, 25, 30) at α = 0.05 or 0.10, and red double lines approximate the straight line of each t0.05(20, 25, 30) and t0.10(20, 25, 30). The red equations on the graph represent a linear approximation line of t-values.
Figure 8. The t-value of α = 0.05 or 0.10 with ν = 2–30. (a) Triangles are t0.05(2–30); (b) diamonds are t0.10(2–30). The dotted lines show (a) y = 1.9600 and (b) y = 1.6449. The red circles represent the t-values of n = (20, 25, 30) at α = 0.05 or 0.10, and red double lines approximate the straight line of each t0.05(20, 25, 30) and t0.10(20, 25, 30). The red equations on the graph represent a linear approximation line of t-values.
Toxins 10 00395 g008
Figure 9. Overlay graphs of f(t, ν = 30, 50, and 100) red lines, and the standard normal distribution f(t, ν = ∞) black lines. (a) f(t, ν =30) and f(t, ν = ∞); (b) f(t, ν = 50) and f(t, ν = ∞); (c) f(t, ν = 100) and f(t, ν = ∞).
Figure 9. Overlay graphs of f(t, ν = 30, 50, and 100) red lines, and the standard normal distribution f(t, ν = ∞) black lines. (a) f(t, ν =30) and f(t, ν = ∞); (b) f(t, ν = 50) and f(t, ν = ∞); (c) f(t, ν = 100) and f(t, ν = ∞).
Toxins 10 00395 g009
Table 1. Sampling information and total weight (g) of each scallop tissue.
Table 1. Sampling information and total weight (g) of each scallop tissue.
201426 May2 June9 June30 June14 July22 July28 July
Number of Individuals16181718152014
Digestive gland72.9072.8772.9770.8060.8075.1658.56
Gonad42.3939.9940.7944.2248.0554.4639.83
Mantle155.01151.37163.48175.80163.66213.96169.30
Gill95.3193.43106.2892.2087.08116.4285.49
Adductor muscle301.69315.19311.42355.14318.80434.35340.26
The concentration of OA and DTX1 in each part was quantified by LC/MS/MS after hydrolysis.
Table 2. Sampling information and mean weight (g) of the digestive glands of 30 scallops or mussels. The mean values of 30 samples ± population standard deviation (σ).
Table 2. Sampling information and mean weight (g) of the digestive glands of 30 scallops or mussels. The mean values of 30 samples ± population standard deviation (σ).
201426 May2 June9 June16 June7 July8 August
Scallop
(Digestive gland/Whole meat %)
3.87 ± 1.07
(10.14%)
3.62 ± 0.81
(10.56%)
3.70 ± 0.98
(8.78%)
3.81 ± 0.72
(9.18%)
--
Mussel
(Digestive gland/Whole meat %)
-1.69 ± 0.45
(14.72%)
--1.40 ± 0.46
(14.89%)
1.51 ± 0.45
(12.90%)
Table 3. Sampling information and mean weight (g) of scallop digestive glands. The mean value of digestive glands ± σ at each depth.
Table 3. Sampling information and mean weight (g) of scallop digestive glands. The mean value of digestive glands ± σ at each depth.
201618 May28 June11 August
Number of Individuals101515
5 m6.45 ± 1.244.77 ± 0.821.40 ± 0.36
10 m6.16 ± 0.915.18 ± 1.331.38 ± 0.26
15 m5.30 ± 0.454.22 ± 0.751.46 ± 0.42
Table 4. Data on environmental conditions and densities of DST producing species at Yakumo Station.
Table 4. Data on environmental conditions and densities of DST producing species at Yakumo Station.
Transparency (m)Date (2016)Depth (m)Water Temperature (°C)Salinity (psu)D. fortii (Cells/L)D. acuminate (Cells/L)Other Dinophysis (Cells/L)
18 May011.230.840300
59.731.980600
5.0 109.232.0301000
158.132.320800
207.832.590300
257.432.670900
307.332.7001500
28 June016.629.58011200
514.831.153065060Dt60
4.0 1013.731.94180740190Dn150, Dt40
1513.432.093049080Dn60, Dr20
2013.032.184030010Dr10
2512.832.2207010Dn10
3012.332.3305010Dn10
11 August022.631.14000
520.731.670050Dt50
10.5 1016.932.3701030Dt30
1512.732.640030Dt20, Dr10
2010.832.951000
258.632.9550010Dt10
307.933.0290200
Dinophysis tripos (Dt); Dinophysis norvegica (Dn); Dinophysis rotundata (Dr).
Table 5. Resampling analysis of scallops and mussels without replacement and with the bootstrap method. The left half of the table is a resampling analysis without replacement, while the right half shows the data using the bootstrap method. The n columns represent 5–25 samples. 1 to 99 represent percentiles. Percentage: each percentile columns represents the mean value of a data set for each mean of the 30 individuals. >±30; <±30; <±20; <±10.
Table 5. Resampling analysis of scallops and mussels without replacement and with the bootstrap method. The left half of the table is a resampling analysis without replacement, while the right half shows the data using the bootstrap method. The n columns represent 5–25 samples. 1 to 99 represent percentiles. Percentage: each percentile columns represents the mean value of a data set for each mean of the 30 individuals. >±30; <±30; <±20; <±10.
ScallopMussel ScallopMussel
n159599159599n159599159599
575.682.1119.0126.349.560.4144.4163.7572.580.7120.7130.344.358.2149.8172.0
677.784.0117.0123.752.063.9139.4155.8675.682.3118.8127.248.861.3144.6164.2
780.085.4114.9120.456.266.9134.8148.8776.683.5116.8124.451.363.8139.8157.2
881.186.5113.9119.058.569.5132.6144.4878.484.7115.7122.954.966.5137.3154.3
982.987.6112.9117.661.971.8129.8140.3979.285.5114.3121.156.567.8133.9150.0
1084.488.5111.6116.064.173.7126.9138.61080.786.3114.0120.159.469.8132.6147.9
1185.589.3111.0114.767.075.9124.8135.41181.486.8113.3119.160.670.4131.4145.4
1286.189.8110.4114.169.477.4123.7132.11282.487.4112.9118.862.572.0130.7144.1
1387.490.7109.4112.970.678.5121.6129.71382.987.9112.2117.663.573.0128.8140.9
1487.891.2108.6112.172.579.7119.9127.71483.688.5112.1117.465.374.2128.7140.8
1588.791.7108.0111.174.281.2118.9126.31584.588.9111.6116.365.675.1126.8137.8
1689.792.5107.6110.375.082.2117.5124.21684.689.2111.0115.666.875.9126.0136.9
1790.092.8107.1109.777.683.7116.5122.41785.189.5110.7115.167.276.6124.9136.2
1890.893.2106.5109.178.584.8115.1120.41885.889.8110.6115.268.977.2124.3135.3
1991.493.6106.2108.480.085.6113.8119.21986.090.0110.2114.769.977.5123.9134.4
2092.094.0105.8107.881.186.7112.9117.32085.890.2109.9114.069.578.0123.2133.6
2192.694.5105.2107.282.287.4111.9116.02186.890.6109.6113.970.578.8122.6132.3
2293.194.9104.9106.683.788.4110.8114.52287.190.8109.4113.671.079.2122.0132.7
2393.595.2104.5106.384.789.2110.1113.52387.291.0109.3113.272.079.8121.6131.1
2494.195.7104.0105.586.190.2108.8111.72487.491.0109.2113.072.379.9121.3130.5
2594.696.1103.6105.087.391.1107.9110.32588.091.3109.1113.173.380.5121.1130.4
Table 6. Student’s t distribution.
Table 6. Student’s t distribution.
νTwo-Tailed Probability
0.100.05
22.92004.3027
32.35343.1824
42.13182.7764
52.01502.5706
61.94322.4469
71.89462.3646
81.85952.3060
91.83312.2622
101.81252.2281
111.79592.2010
121.78232.1788
131.77092.1604
141.76132.1448
151.75312.1314
161.74592.1199
171.73962.1098
181.73412.1009
191.72912.0930
201.72472.0860
211.72072.0796
221.71712.0739
231.71392.0687
241.71092.0639
251.70812.0595
261.70562.0555
271.70332.0518
281.70112.0484
291.69912.0452
301.69732.0423
501.67592.0086
1001.66021.9840
1.64491.9600
The two-tailed probability is 0.10 or 0.05.
Table 7. Interval estimation (±OA group ng/g digestive gland) of μ in α = 0.05. The n columns represent sample size 3–20 samples. Rows from 100 to 650 represent the s values.
Table 7. Interval estimation (±OA group ng/g digestive gland) of μ in α = 0.05. The n columns represent sample size 3–20 samples. Rows from 100 to 650 represent the s values.
nt (α = 0.05)Sample Standard Deviation (s)
100150200250300350400450500550600650
34.3027248.4372.6496.8621.0745.2869.5993.71117.91242.11366.31490.51614.7
43.1825159.1238.7318.3397.8477.4556.9636.5716.1795.6875.2954.81034.3
52.7764124.2186.2248.3310.4372.5434.6496.7558.7620.8682.9745.0807.1
62.5706104.9157.4209.9262.4314.8367.3419.8472.2524.7577.2629.7682.1
72.446992.5138.7185.0231.2277.5323.7369.9416.2462.4508.7554.9601.1
82.364683.6125.4167.2209.0250.8292.6334.4376.2418.0459.8501.6543.4
92.306076.9115.3153.7192.2230.6269.0307.5345.9384.3422.8461.2499.6
102.262271.5107.3143.1178.8214.6250.4286.1321.9357.7393.5429.2465.0
112.228167.2100.8134.4167.9201.5235.1268.7302.3335.9369.5403.1436.7
122.201063.595.3127.1158.8190.6222.4254.1285.9317.7349.5381.2413.0
132.178860.490.6120.9151.1181.3211.5241.7271.9302.1332.4362.6392.8
142.160457.786.6115.5144.3173.2202.1231.0259.8288.7317.6346.4375.3
152.144855.483.1110.8138.4166.1193.8221.5249.2276.9304.6332.3360.0
162.131553.379.9106.6133.2159.9186.5213.2239.8266.4293.1319.7346.4
172.119951.477.1102.8128.5154.2180.0205.7231.4257.1282.8308.5334.2
182.109849.774.699.5124.3149.2174.0198.9223.8248.6273.5298.4323.2
192.100948.272.396.4120.5144.6168.7192.8216.9241.0265.1289.2313.3
202.093046.870.293.6117.0140.4163.8187.2210.6234.0257.4280.8304.2
Table 8. Interval estimation (±OA group ng/g digestive gland) of μ in α = 0.10. The n columns represent sample size 3–20. Rows 100 to 650 represent the s values.
Table 8. Interval estimation (±OA group ng/g digestive gland) of μ in α = 0.10. The n columns represent sample size 3–20. Rows 100 to 650 represent the s values.
nt (α = 0.10)Sample Standard Deviation (s)
100150200250300350400450500550600650
32.9200168.6252.9337.2421.5505.8590.1674.3758.6842.9927.21011.51095.8
42.3534117.7176.5235.3294.2353.0411.8470.7529.5588.4647.2706.0764.9
52.131895.3143.0190.7238.3286.0333.7381.3429.0476.7524.4572.0619.7
62.015082.3123.4164.5205.7246.8287.9329.0370.2411.3452.4493.6534.7
71.943273.4110.2146.9183.6220.3257.1293.8330.5367.2404.0440.7477.4
81.894667.0100.5134.0167.5201.0234.4267.9301.4334.9368.4401.9435.4
91.859562.093.0124.0155.0186.0216.9247.9278.9309.9340.9371.9402.9
101.833158.087.0115.9144.9173.9202.9231.9260.9289.8318.8347.8376.8
111.812554.682.0109.3136.6163.9191.3218.6245.9273.2300.6327.9355.2
121.795951.877.8103.7129.6155.5181.5207.4233.3259.2285.1311.1337.0
131.782349.474.198.9123.6148.3173.0197.7222.4247.2271.9296.6321.3
141.770947.371.094.7118.3142.0165.7189.3213.0236.6260.3284.0307.6
151.761345.568.291.0113.7136.4159.2181.9204.6227.4250.1272.9295.6
161.753043.865.787.7109.6131.5153.4175.3197.2219.1241.0263.0284.9
171.745942.363.584.7105.9127.0148.2169.4190.5211.7232.9254.1275.2
181.739641.061.582.0102.5123.0143.5164.0184.5205.0225.5246.0266.5
191.734139.859.779.699.5119.3139.2159.1179.0198.9218.8238.7258.6
201.729138.758.077.396.7116.0135.3154.7174.0193.3212.7232.0251.3

Share and Cite

MDPI and ACS Style

Matsushima, R.; Uchida, H.; Watanabe, R.; Oikawa, H.; Oogida, I.; Kosaka, Y.; Kanamori, M.; Akamine, T.; Suzuki, T. Anatomical Distribution of Diarrhetic Shellfish Toxins (DSTs) in the Japanese Scallop Patinopecten yessoensis and Individual Variability in Scallops and Mytilus edulis Mussels: Statistical Considerations. Toxins 2018, 10, 395. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins10100395

AMA Style

Matsushima R, Uchida H, Watanabe R, Oikawa H, Oogida I, Kosaka Y, Kanamori M, Akamine T, Suzuki T. Anatomical Distribution of Diarrhetic Shellfish Toxins (DSTs) in the Japanese Scallop Patinopecten yessoensis and Individual Variability in Scallops and Mytilus edulis Mussels: Statistical Considerations. Toxins. 2018; 10(10):395. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins10100395

Chicago/Turabian Style

Matsushima, Ryoji, Hajime Uchida, Ryuichi Watanabe, Hiroshi Oikawa, Izumi Oogida, Yuki Kosaka, Makoto Kanamori, Tatsuro Akamine, and Toshiyuki Suzuki. 2018. "Anatomical Distribution of Diarrhetic Shellfish Toxins (DSTs) in the Japanese Scallop Patinopecten yessoensis and Individual Variability in Scallops and Mytilus edulis Mussels: Statistical Considerations" Toxins 10, no. 10: 395. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins10100395

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