Next Article in Journal
Temporal Variation and Industry-Specific Differences of the Use of Volatile Organic Compounds from 2018 to 2023 and Their Health Risks in a Typical Industrially Concentrated Area in South China
Previous Article in Journal
The Association between Metals and Thyroid Cancer in Puerto Rico—A National Health and Nutrition Examination Survey Analysis and Ecological Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Robust Method for Simultaneous Determination and Risk Assessment of Multiresidual Pesticides in Fishery Products

1
Department of Food Engineering, Daegu University, Gyeongsan 38453, Republic of Korea
2
Analysis Technology and Tomorrow, Daegu 42703, Republic of Korea
3
Food Safety and Processing Research Division, National Institute Fisheries Science, Busan 46083, Republic of Korea
4
Department of Horticulture and Life Science, Yeungnam University, Gyeongsan 38541, Republic of Korea
*
Author to whom correspondence should be addressed.
Toxics 2024, 12(9), 633; https://doi.org/10.3390/toxics12090633
Submission received: 5 August 2024 / Revised: 26 August 2024 / Accepted: 26 August 2024 / Published: 28 August 2024
(This article belongs to the Section Agrochemicals and Food Toxicology)

Abstract

:
In this study, we developed and validated a multiresidue analytical method for the simultaneous detection of 24 pesticides in fishery products. Using the EN15662 extraction method and C18 as the adsorbent for purification, the validation results complied with Codex guidelines, achieving recovery rates between 70% and 120% and relative standard deviation values (%RSD) within 20%, indicating excellent performance. The limit of detection ranged from 0.25 to 0.8 ng/kg, and the limit of quantification was between 3 and 10 ng/g, providing sufficient sensitivity to comply with future regulatory standards. The calibration curves for all 24 pesticides exhibited great linearity (R2 > 0.98), also satisfying the Codex requirements. The matrix effect was less than 30% for some pesticides—within ±20%—indicating minimal interference from impurities. An analysis of 300 fishery samples from nine regions across South Korea detected lufenuron at 10 ng/g in eels; however, the risk assessment was below 0.19%, posing no significant hazard to public health. This newly developed analytical method proved effective for the multi-analysis of pesticide residues in fishery products, offering rapid and reliable monitoring of the import and export safety of fishery products.

1. Introduction

The need for alternative protein sources has increased owing to the limited productivity of livestock and the resulting shortage of protein supply. With the increase in the production of cultured fish as an alternative protein source, fishery products have become a significant component of the global food industry. Fish consumption is recommended for its high content of omega-3 fatty acids, proteins, polyunsaturated fatty acids, and vitamins, which together help prevent various diseases, including cardiovascular disorders [1,2,3,4]. According to the Food and Agriculture Organization (FAO) of the United Nations, South Korea ranks within the top 3% of 239 countries worldwide for fishery product consumption [1,5]. Fishery production in South Korea is classified into four categories: marine aquaculture, coastal fisheries, deep-sea fisheries, and inland fisheries [6]. To produce high-quality fishery products, drugs are used in fish farming for promoting growth and preventing or treating diseases [7,8]. However, pesticides used in nearby agricultural areas or golf courses can contaminate surrounding rivers, reservoirs, and lakes, leading to unintended fish contamination [9,10]. Since most pesticides are toxic synthetic organic compounds, the consumption of food contaminated due to inadequate pesticide residue management can harm human health [11,12,13,14].
Organochlorine pesticides have been globally restricted since the early 1970s owing to their potential to cause congenital malformations, reproductive and endocrine disorders, immune system dysfunction, and cancer [15]. These compounds are structurally stable, resistant to degradation, exhibit long half-lives, and display high octanol–water partition coefficients, which together lead to bioaccumulation in the muscles and fat tissues of aquatic organisms, resulting in chronic toxicity in humans through long-term exposure [12,16,17]. Consequently, systematic management of aquatic environments and fish farms, which requires effective analytical methods for monitoring the presence of pesticide residues in these ecosystems, is essential to prevent pesticide residue contamination [18].
Worldwide, South Korea, along with the United States, Australia, the European Union, and Japan, have established pesticide maximum residue levels (MRLs) for domestic and imported fishery products to ensure food safety. Additionally, the Codex Alimentarius Commission (Codex) has set international standards to ensure fairness in global food trade.
Regarding pesticide management in food, the positive list system (PLS) was implemented for all agricultural products beginning 1 January 2019, applying a limit of detection (LOD) of 0.01 mg/kg for food products for which no pesticide MRLs have been established. Since 1 January 2024, the PLS of veterinary medicines has been applied to livestock products, including cattle, pigs, chickens, milk, and eggs. Following this trend, the PLS of pesticides can also be applied to fishery products.
Various analytical methods, such as solid–liquid extraction, Soxhlet extraction, sonication-assisted extraction, solid-phase extraction, and accelerated solvent extraction related to gel permeation chromatography, have been introduced [19,20,21]. However, these methods are often expensive, labor intensive, and time consuming; furthermore, they show low sample throughput and require advanced analytical expertise and expensive specialized equipment [22]. To address these issues, the quick, easy, cheap, effective, rugged, and safe (QuEChERS) method was developed, and it has been improved to suit various compounds and matrices [23,24,25]. Recent research includes a study on pesticide residue analysis in tilapia using the QuEChERS method [26], showing an improved multiresidue analysis method for 32 pesticides in aquaculture fish using acetonitrile with 1% formic acid and acetonitrile/methanol (5:5 v/v), with 1% formic acid as the extraction solvent [27]. Another study compared three extraction and purification methods for shrimp and eel [28]. However, most studies have focused on a single fish species, with limited research on methods that can simultaneously analyze multiple types of fishery products.
In this study, we aimed to develop an effective analytical method for the simultaneous analysis of 24 pesticides in fishery products, including triazines, triazoles, amides, carbamates, neonicotinoids, and anilines, by referencing and modifying the QuEChERS method. The new method was successfully applied to various types of fishery products, followed by comprehensive monitoring and risk assessment.

2. Materials and Methods

2.1. Reagents

The 24 pesticides analyzed in this study (acetamiprid, azimsulfuron, azinphos-methyl, azoxystrobin, carbendazim, carbofuran, clothianidin, daimuron, dichlorvos, difenoconazole, dinotefuran, diuron, fenobucarb, flubendiamide, hexaconazole, indoxacarb, lufenuron, pyrimisulfan, tebufenozide, terbuthylazine, thiacloprid, thiamethoxam, tricyclazole, and trifloxystrobin) were purchased from Kemidas (Gunpo, Republic of Korea). Acetonitrile (HPLC-grade) was obtained from J. T. Baker (Center Valley, PA, USA), and acetic acid (99%) was purchased from Duksan Chemicals (Ansan, Republic of Korea). Magnesium sulfate anhydrous (MgSO4; 99.5%) and sodium chloride (NaCl; 99.5%) were acquired from Junsei (Tokyo, Japan), and sodium citrate dihydrate (Na3Cit∙1.5H2O; ≥99.0) and sodium acetate (NaOAc; ≥99.0) were purchased from Sigma-Aldrich (St. Louis, MO, USA). Primary secondary amine (PSA) and C18 were acquired from Agilent Technology (Santa Clara, CA, USA).

2.2. Sample Preparation

Fishery product samples were purchased from fish farms in nine regions across the country: Chungcheongbuk-do, Chungcheongnam-do, Jeju-do, Jeollabuk-do, Jeollanam-do, Gangwon-do, Gyeonggi-do, Gyeongsangbuk-do, and Gyeongsangnam-do. Following the Korean Food Code, 300 samples were collected, including 20 samples each of 13 fish species (carp, Chinese muddy loach, crucian carp, eel, far eastern catfish, flathead mullet, Korean rockfish, mirror carp, olive flounder, rainbow trout, red sea bream, sea bass, starry flounder), one crustacean species (whiteleg shrimp), and one shellfish species (abalone) [29]. Crustaceans and shellfish were shelled and eviscerated, and fish samples were prepared by removing fins, bones, and heads, followed by homogenization of the edible parts, including the skin, with dry ice using a mixer (Grinmic gold-DA10000G, Daesung Artlon, Paju, Republic of Korea). All samples were stored at −20 °C until analysis.

2.3. Optimization of Analytical Methods

We used eel, chosen for its high fat content, as a control sample, treating it with 10 ng/g of a mixed pesticide standard. The liquid–liquid partitioning efficiencies of the EN15662 and AOAC methods were compared to establish the extraction method [30]. The EN15662 method used 20 mL of acetonitrile (+0.1% acetic acid), 4 g of MgSO4, 1 g of NaCl, 1.5 g of Na3Cit∙1.5H2O, and 5 mL of water. The AOAC method used 20 mL of acetonitrile (+0.1% acetic acid), 6 g of MgSO4, and 1.5 g of NaOAc. For purification, the efficacy of 50 mg of PSA, combined with 50 mg of C18, was compared to that of 50 mg of C18 alone. The optimal analytical method was determined based on the results of the recovery tests.

2.4. LC–MS/MS Analytical Conditions

Pesticide residues in fishery products were analyzed using an LC–MS/MS system (LCMS8060NX, Shimadzu, Kyoto, Japan), and data were processed using Labsolution SP2 software (version 5.75) (LCMS8060NX, Shimadzu, Kyoto, Japan). LC was performed using a C18 column (2.1 mm × 150 mm, 3.0 μm, Shimadzu Shim-pack GIST-HP) (Kyoto, Japan). Two mobile phases, A (aqueous, 0.1% formic acid, and 5 mM ammonium formate) and B (methanol, 0.1% formic acid, and 5 mM ammonium formate), were used at a flow rate of 0.2 mL/min. The gradient started at 5% B for 1 min, increased to 60% B over 2 min, then to 100% B within 10 min, and was maintained for 3 min. Finally, it was returned to 5% B over 0.1 min and held for 3.9 min, for a total of 20 min. MS/MS was conducted in the positive electron spray ionization mode with an interface voltage of 4.0 kV (Table 1). Multiple reaction monitoring conditions were established for each of the 24 pesticides by identifying the precursor ions and scanning for product ions. Based on the sensitivity and selectivity, at least two product ions for each pesticide were selected for quantitative and qualitative analyses, and their optimal collision energies were determined.

2.5. Method Validation

The newly developed analytical method presented in this study was validated according the Codex (2023) guidelines [31], assessing linearity, LOD, limit of quantitation (LOQ), accuracy, and precision. Olive flounder, eel, abalone, and whiteleg shrimp with no pesticide residues were used as control samples. Linearity was determined using matrix-matched calibration curves, with calculation of the coefficient of determination (R2). LOD and LOQ were determined by analyzing pesticide-free samples and calculating the signal-to-noise (S/N) ratio of the chromatographic peaks. LOD was defined as the concentration with an S/N ratio of 3 or higher, and LOQ was defined as that with an S/N ratio of 10 or higher. Recovery tests were conducted to assess accuracy and precision by spiking the control samples at the LOQ, 10 times the LOQ, and 50 times the LOQ with a mixed standard solution, followed by five replicate analyses.

2.6. Risk Assessment

The risk assessment was conducted based on the method reported by Kim et al. [32]. Estimated daily intake (EDI) for the detected pesticides and fishery products consumption data were calculated using nine different scenarios (Equation (1)), followed by exposure assessment.
EDI (ng/person/day) = DFI a (g/person/day) (1, 2, 3) × DPC b (ng/g) (A, B, C)
where a—daily food intake, b—detected pesticide concentration
1—average intake by fish species,
2—extreme (99th percentile) intake by detected fish,
3—extreme (99th percentile) intake by fish species,
A—(detected pesticide + LOQ value for non-detection)/number of test (20),
B—detected pesticide/number of detections, and
C—maximum detected pesticide.
Consumption data were based on the Korea National Health and Nutrition Examination Survey (KNHANES) for 2017–2021 [33], considering the average and extreme (99th percentile) intake of fish and shellfish (Table 2). For species with no available consumption data, intake values were estimated as half the minimum mean and extreme (99th percentile) intake of similar species. The average body weight of Koreans was assumed to be 60 kg. The health risk index (%ADI) was calculated by comparing the EDI with the acceptable daily intake (ADI) and expressing it as a percentage of ADI (Equation (2)).
% ADI = EDI   ng / person / day ADI   ng / person / day × 100

2.7. Statistical Analysis

Statistical analyses were performed using Excel 2013 (Microsoft Corporation, Seattle, WA, USA). Data are presented as the mean and standard deviation, resulting from the analysis of five replicates.

3. Results and Discussion

3.1. Specificity

The total ion chromatograms (TIC) of the 24 pesticides analyzed are shown in Figure 1. The qualitative and quantitative ions of these pesticides were sufficiently separated, allowing for distinct identification based on their retention times. Pesticides with similar retention times were differentiated by their m/z values for precursor and product ions.

3.2. Establishment of Sample Extraction and Purification Methods

The efficiency of pesticide extraction from fishery products was evaluated by comparing the EN15662 [34] and AOAC [30] methods. Additionally, to assess the effectiveness of removing interfering substances, such as fats, during purification, two purification approaches were compared: the use of C18 (50 mg) alone and a combination of PSA (50 mg) and C18 (50 mg). Table 3 summarizes the extraction and purification effectiveness results.
The EN15662 method demonstrated recoveries above 70% for all 24 pesticides, whereas the AOAC method resulted in recoveries below 60% for 8 pesticides, including carbendazimhad. Kwak et al. [35] and Lehotay et al. [36] reported that the acetate buffer in the AOAC method, which lowers sample pH, hampered pesticide transfer to acetonitrile, reducing recovery rates. Additionally, the EN15662 method showed a lower matrix effect than the AOAC method. Consequently, the EN15662 method, which demonstrated superior recoveries, was selected for the extraction. Since the amine groups of the PSA adsorbent interacted with the polar parts of azimsulfuron and thiamethoxam through ionic interactions or hydrogen bonding, PSA/C18 led to the adsorption of azimsulfuron and thiamethoxam during the purification, as well as lower recovery [37]. Purification with C18 alone achieved recoveries above 70% for all pesticides and was expected to effectively remove impurities, such as sugars, fatty acids, and organic acids, from the samples [38].
The EN15662 extraction method and C18 purification were selected for optimal performance, as they met the Codex guidelines for recovery rates.

3.3. Final Extraction and Purification

The optimized sample preparation procedure is detailed in Figure 2. Briefly, 5 g of the homogenized sample was added to a 50 mL Teflon tube, followed by 5 mL of water and 20 mL of acetonitrile containing 0.1 acetic acid. After shaking for 30 min, 4 g of MgSO4, 1 g of NaCl, 1 g of Na3Cit, and 1.5 g of Na3Cit·1.5H2O were added, and the mixture was shaken for an additional 30 min. The mixture was centrifuged at 4000 rpm for 5 min at 4 °C using a Megafuge 1.0 centrifuge (Thermo Fisher Scientific, Waltham, MA, USA). The acetonitrile and water layers were separated, and 1 mL of the supernatant (acetonitrile) was transferred to a dispersive solid-phase extraction (d-SPE) tube containing 50 mg of C18. After vortexing for 1 min and centrifugation at 10,000 rpm for 5 min, the supernatant was filtered through a 0.2 μm syringe filter to prepare the test solution for LC–MS/MS analysis.

3.4. Recovery, Linearity, LOD, LOQ, Matrix Effect

Method validation was conducted according to Codex guidelines [31] (Table 4). Recovery experiments were performed on representative fish species by spiking 24 pesticides at three different concentrations (LOQ, 10 × LOQ, 50 × LOQ) using the established method. Recovery rates for all the fortified samples were above 80%, meeting the Codex acceptance range (70–120%) for multiresidue analysis recovery.
Linearity was determined using standard calibration curves produced by LC–MS/MS. Matrix-matched calibration curves at seven different concentrations, including LOQ levels, were analyzed using LC–MS/MS. The concentration ranges for the calibration curves were 0.25, 0.75, 1.5, 3.75, 7.5, 12, 15 ng/g for 3 ng/g LOQ; 0.3, 1, 1.6, 2, 5, 10, 16, 20 ng/g for 4 ng/g LOQ; 0.4, 1.25, 2, 2.5, 6.25, 12.5, 20, 25 ng/g for 5 ng/g LOQ; 0.6, 1.75, 2.8, 3.5, 8.75, 17.5, 28, 35 ng/g for 7 ng/g LOQ; and 0.8, 2.5, 4, 5, 12.5, 25, 40, 50 ng/g for 10 ng/g LOQ. The coefficients of determination (R2) values for the 24 pesticides ranged from 0.9921 to 1.000, demonstrating excellent linearity in accordance with Codex guidelines (R2 > 0.98).
The LOD (S/N ≥ 3) ranged from 0.25 to 0.8 ng/g, and the LOQ (S/N ≥ 10) from 3 to 10 ng/g, depending on the pesticides and samples. Some pesticides, such as carbendazim, dichlorvos, and thiamethoxam, showed a higher LOQ (7–10 ng/g) owing to lower sensitivity of the LC–MS/MS. These LOQ results are sufficient for future applications in fishery product pesticide residue monitoring under PLS regulations.
To evaluate the matrix effect (ME) on pesticide ionization, matrix-matched calibration curves were compared to solvent-based curves. ME values between −20% and +20% indicated minimal interference from co-eluting substances, suggesting a weak matrix effect. ME values outside this range (>+50% or <−50%) indicated significant interference. Overall, most pesticides displayed low ME values, with only a few exceptions exceeding the ±20% range in specific samples, such as azinphos-methyl and carbofuran in eel, and tricyclazole in olive flounder and whiteleg shrimp.

3.5. Accuracy and Precision

Accuracy was evaluated based on recovery rates, whereas precision was assessed using the relative standard deviation (% RSD). Standard solutions of 24 pesticides at concentrations of LOQ, 10 × LOQ, and 50 × LOQ were spiked into untreated samples, and each sample was analyzed five times. Recoveries ranged from 71.7% to 106.3% at all concentration levels, meeting the Codex accuracy guidelines (70–120%). Precision (%RSD) ranged from 1.0% to 14.2% at LOQ, 0.6% to 9.0% at 10 × LOQ, and 0.4% to 6.9% at 50 × LOQ, all within the Codex requirement of ≤20%.

3.6. Monitoring Results for Pesticides

Based on the monitoring results, lufenuron, a benzoylurea insecticide, was detected at 10 ng/g in one out of the 300 eel samples analyzed. This pesticide is known for its non-systemic properties and inhibition of chitin synthesis [39], which effectively controls pests, such as diamondback moths. In South Korea, lufenuron is used to control pests in various agricultural products including sweet potatoes, peppers, blueberries, cherries, oranges, and grapes [40]. Pesticides with lufenuron as the active ingredient might have entered the aquatic environment and accumulated in eels. This finding underscores the necessity for the continuous monitoring of pesticide residues in fishery products to ensure food safety.

3.7. Risk Assessment for Fishery Product Consumption

Risk assessment was based on pesticide residue detection and consumption data. The EDI and health risk index (%ADI) for various scenarios are shown in Figure 3. With an ADI of 0.015 mg/kg bw/day for lufenuron, the %ADI for average consumers ranged from 0.0017% to 0.0072 (scenarios 1–3). For extreme consumers of eel in which lufenuron was detected, the %ADI values (scenarios 4–6) ranged from 0.0141% to 0.0269%. The %ADI values for extreme consumers of different fishery product species (scenarios 7–9) ranged from 0.0337% to 0.1917%, with the highest value (scenario 9) still indicating a very low risk. Overall, the risk assessment across all the scenarios indicated that the detected pesticide residues posed a very low risk. According to the FAO and WHO guidelines, a %ADI below 10% signifies minimal risk [41]. Kim et al. [32] also reported low risk levels (0.00% to 1.07%) for farmed fish in South Korea. This study confirms that even a high consumption of contaminated fishery products poses a very low risk to human health. Nonetheless, the detection of pesticides that are not intended for fish highlights the need for continued monitoring and systematic management of pesticide residues.

4. Conclusions

We developed and validated an efficient pretreatment and analytical method for the detection of pesticide residues in fishery products. The EN15662 method achieved higher recovery rates than did the AOAC method, attaining a greater than 70% recovery rate for all the 24 pesticides analyzed. Although the combination of PSA and C18 resulted in lower recovery rates for some pesticides, the use of C18 alone consistently resulted in higher recovery rates for all pesticides. The newly developed analytical method complies with the Codex guidelines, demonstrating high accuracy, precision, and minimal matrix effects, which ensures minimal ionization interference from the sample matrix. This demonstrates that our method is suitable for the analysis of pesticide residues in various fishery products.
Risk assessment analysis indicated that all scenarios fell well below the 10% safety threshold set by the FAO and WHO, suggesting a very low potential risk to human health from the consumption of fishery products. However, the detection of pesticides not intended for use with fish underscores the importance of improving the safety management and monitoring of fishery products.
The method developed in this study provides essential baseline data for the safety assessment and management of fishery products and contributes to the safety of these products through continuous monitoring. Additionally, our findings support maintenance of the safety of fishery products in international trade and can serve as a valuable reference for establishing related regulations and standards.

Author Contributions

Conceptualization: M.K., Y.L., M.-R.J. and M.-H.I.; formal analysis: M.K., T.-h.K. and J.-W.P.; funding acquisition: Y.L. and M.-R.J.; investigation: M.K. and Y.-S.M.; methodology: M.K., T.-h.K. and J.-W.P.; project administration: M.-H.I.; supervision: M.-H.I.; validation: M.K., T.-h.K. and J.-W.P.; writing—original draft: M.K.; writing—review and editing: Y.-S.M. and M.-H.I. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Institute of Fisheries Science, Ministry of Oceans and Fisheries, Korea (grant number R2024057).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available upon request from the corresponding author owing to legal restrictions.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. FAO. The State of World Fisheries and Aquaculture 2022. In Towards Blue Transformation; FAO: Rome, Italy, 2022. [Google Scholar] [CrossRef]
  2. Magalhães, K.M.; Carreira, R.S.; Rosa Filho, J.S.; Rocha, P.P.; Santana, F.M.; Yogui, G.T. Polycyclic aromatic hydrocarbons (PAHs) in fishery resources affected by the 2019 oil spill in Brazil: Short-term environmental health and seafood safety. Mar. Pollut. Bull. 2022, 175, 113334. [Google Scholar] [CrossRef]
  3. Udofia, U.S.; Ameh, C.; Miller, E.; Ekpenyong, M.S. Investigating the origin and tissue concentration of polycyclic aromatic hydrocarbons in seafood and health risk in Niger Delta, Nigeria. Environ. Sci. Process. Impacts 2021, 23, 1803–1814. [Google Scholar] [CrossRef] [PubMed]
  4. Shariatifar, N.; Sharifiarab, G.; Kargarghomsheh, P.; Moazzen, M.; Arabameri, M.; Seddighi, M.; Pirhadi, M. Polycyclic aromatic hydrocarbons (PAHs) in potato and related products in Tehran: A health risk assessment study. Int. J. Environ. Anal. Chem. 2022, 1–14. Available online: https://www.tandfonline.com/doi/citedby/10.1080/03067319.2022.2100258?scroll=top&needAccess=true (accessed on 30 April 2022). [CrossRef]
  5. MOF (Ministry of Oceans and Fisheries). Seafood Consumption (Per Person Per Year) and Self-Sufficiency. Available online: https://www.index.go.kr/potal/main/EachDtlPageDetail.do?idx_cd=1317 (accessed on 28 January 2022).
  6. MOF (Ministry of Oceans and Fisheries). Fishery Production and Aquaculture. Available online: https://www.index.go.kr/unity/potal/main/EachDtlPageDetail.do?idx_cd=2748 (accessed on 28 April 2022).
  7. MFDS (Ministry of Food and Drug Safety). The Veterinary Drugs Monitoring and Their Analytical Methods Development in Fishery Products from Domestic Market. 2017. Available online: https://scienceon.kisti.re.kr/srch/selectPORSrchReport.do?cn=TRKO201800035693 (accessed on 4 August 2024). [CrossRef]
  8. Oh, J.H.; Kwon, C.H.; Jeon, J.S.; Choi, D.M. Management of veterinary drug residues in food. Korean J. Environ. Agric. 2009, 28, 310–325. [Google Scholar] [CrossRef]
  9. KHIDI (Korea Health Industry Development Institute). Intake by Food Group. 2022. Available online: https://www.khidi.or.kr/eps (accessed on 30 April 2022).
  10. Kim, D.; Yoon, J.; Yoo, J.; Kim, S.J.; Yang, J.E. Status and management strategy of pesticide use in golf courses in Korea. J. Appl. Biol. Chem. 2014, 57, 267–277. [Google Scholar] [CrossRef]
  11. Aktar, W.; Sengupta, D.; Chowdhury, A. Impact of pesticides use in agriculture: Their benefits and hazards. Interdiscip. Toxicol. 2009, 2, 1–12. [Google Scholar] [CrossRef]
  12. Afful, S.; Anim, A.K.; Serfor-Armah, Y. Spectrum of organochlorine pesticide residues in fish samples from the Densu Basin. Res. J. Environ. Earth Sci. 2010, 2, 133–138. [Google Scholar]
  13. Hwang, I.S.; Oh, Y.J.; Kwon, H.Y.; Ro, J.H.; Kim, D.B.; Moon, B.C.; Lee, H.S. Monitoring of pesticide residues concerned in stream water. Korean J. Environ. Agric. 2019, 38, 173–184. [Google Scholar] [CrossRef]
  14. Ccanccapa, A.; Masiá, A.; Navarro-Ortega, A.; Picó, Y.; Barceló, D. Pesticides in the Ebro River basin: Occurrence and risk assessment. Environ. Pollut. 2016, 211, 414–424. [Google Scholar] [CrossRef]
  15. Garabrant, D.H.; Held, J.; Langholz, B.; Peters, J.M.; Mack, T.M. DDT and related compounds and risk of pancreatic cancer. JNCI J. Natl. Cancer Inst. 1992, 84, 764–771. [Google Scholar] [CrossRef]
  16. Magnusson, K.; Ekelund, R.; Grabic, R.; Bergqvist, P.A. Bioaccumulation of PCB congeners in marine benthic infauna. Mar. Environ. Res. 2006, 61, 379–395. [Google Scholar] [CrossRef] [PubMed]
  17. Adami, H.O.; Lipworth, L.; Titus-Ernstoff, L.; Hsieh, C.C.; Hanberg, A.; Ahlborg, U.; Trichopoulos, D. Organochlorine compounds and estrogen-related cancers in women. Cancer Causes Control 1995, 6, 551–566. [Google Scholar] [CrossRef] [PubMed]
  18. Lazartigues, A.; Fratta, C.; Baudot, R.; Wiest, L.; Feidt, C.; Thomas, M.; Cren-Olivé, C. Multiresidue method for the determination of 13 pesticides in three environmental matrices: Water, sediments and fish muscle. Talanta 2011, 85, 1500–1507. [Google Scholar] [CrossRef] [PubMed]
  19. Forsberg, N.D.; Wilson, G.R.; Anderson, K.A. Determination of parent and substituted polycyclic aromatic hydrocarbons in high-fat salmon using a modified QuEChERS extraction, dispersive SPE and GC–MS. J. Agric. Food Chem. 2011, 59, 8108–8116. [Google Scholar] [CrossRef]
  20. Berntssen, M.H.; Måge, A.; Julshamn, K.; Oeye, B.E.; Lundebye, A.K. Carry-over of dietary organochlorine pesticides, PCDD/Fs, PCBs, and brominated flame retardants to Atlantic salmon (Salmo salar L.) fillets. Chemosphere 2011, 83, 95–103. [Google Scholar] [CrossRef]
  21. Schenck, F.J.; Calderon, L.; Podhorniak, L.V. Determination of organochlorine pesticide and polychlorinated biphenyl residues in fatty fish by tandem solid-phase extraction cleanup. J. AOAC Int. 1996, 79, 1209–1214. [Google Scholar] [CrossRef]
  22. Holmes, B.; Dunkin, A.; Schoen, R.; Wiseman, C. Single-laboratory ruggedness testing and validation of a modified QuEChERS approach to quantify 185 pesticide residues in salmon by liquid chromatography–and gas chromatography–tandem mass spectrometry. J. Agric. Food Chem. 2015, 63, 5100–5106. [Google Scholar] [CrossRef]
  23. Lazartigues, A.; Wiest, L.; Baudot, R.; Thomas, M.; Feidt, C.; Cren-Olivé, C. Multiresidue method to quantify pesticides in fish muscle by QuEChERS-based extraction and LC-MS/MS. Anal. Bioanal. Chem. 2011, 400, 2185–2193. [Google Scholar] [CrossRef]
  24. Rawn, D.F.; Judge, J.; Roscoe, V. Application of the QuEChERS method for the analysis of pyrethrins and pyrethroids in fish tissues. Anal. Bioanal. Chem. 2010, 397, 2525–2531. [Google Scholar] [CrossRef]
  25. Chamkasem, N.; Ollis, L.W.; Harmon, T.; Lee, S.; Mercer, G. Analysis of 136 pesticides in avocado using a modified QuEChERS method with LC-MS/MS and GC-MS/MS. J. Agric. Food Chem. 2013, 61, 2315–2329. [Google Scholar] [CrossRef]
  26. Abbassy, M.A.; Khalifa, M.A.; Nassar, A.M.; El-Deen, E.E.N.; Salim, Y.M. Analysis of organochlorine pesticides residues in fish from Edko Lake (North of Egypt) using eco-friendly method and their health implications for humans. Toxicol. Res. 2021, 37, 495–503. [Google Scholar] [CrossRef]
  27. Kim, K.; Choi, Y.; Mok, S.; Moon, H.B.; Jeon, J. Optimization of the QuEChERS method for multi-residue analysis of pharmaceuticals and pesticides in aquaculture products. Food Chem. 2023, 399, 133958. [Google Scholar] [CrossRef]
  28. Cho, H.R.; Park, J.S.; Kim, J.; Han, S.B.; Choi, Y.S. Multiresidue method for the quantitation of 20 pesticides in aquatic products. Anal. Bioanal. Chem. 2015, 407, 9043–9052. [Google Scholar] [CrossRef] [PubMed]
  29. Ministry of Food and Drug Safety (MFDS). Food Code (N. 2021-54). Available online: https://www.mfds.go.kr/eng/brd/m_15/view.do?seq=72437&srchFr=&srchTo=&srchWord=&srchTp=&itm_seq_1=0&itm_seq_2=0&multi_itm_seq=0&company_cd=&company_nm=&page=1 (accessed on 7 June 2022).
  30. Lehotay, S.J. Determination of pesticide residues in foods by acetonitrile extraction and partitioning with magnesium sulfate: Collaborative study. J. AOAC Int. 2007, 90, 485–520. [Google Scholar] [CrossRef] [PubMed]
  31. FAO; WHO. Codex Alimentarius Commission Procedural Manual, 28th ed.; Revised; FAO: Rome, Italy; WHO: Geneva, Switzerland, 2023. [Google Scholar] [CrossRef]
  32. Kim, M.; Cho, M.; Kim, S.H.; Lee, Y.; Jo, M.R.; Moon, Y.S.; Im, M.H. Monitoring and risk assessment of pesticide residues in fishery products using GC–MS/MS in South Korea. Toxics 2024, 12, 299. [Google Scholar] [CrossRef]
  33. KNHANES (Korea National Health and Nutrition Examination Survey). Available online: https://knhanes.kdca.go.kr/knhanes/main.do (accessed on 29 July 2021).
  34. BS EN 15662; Foods of Plant Origin-Determination of Pesticide Residues Using GC-MS and. or LC-MS/MS Following Acetonitrile Extraction/Partitioning and Clean-up by Dispersive SPE-QuEChERS-Method. 2008. Available online: http://esearch.cen.euesearch (accessed on 29 July 2021).
  35. Kwak, S.Y.; Lee, S.H.; Jeong, H.R.; Nam, A.J.; Sarker, A.; Kim, H.Y.; Kim, J.E. Variation of pesticide residues in strawberries by washing and boiling processes. Korean J. Environ. Agric. 2019, 38, 281–290. [Google Scholar] [CrossRef]
  36. Lehotay, S.J.; Son, K.A.; Kwon, H.; Koesukwiwat, U.; Fu, W.; Mastovska, K.; Leepipatpiboon, N. Comparison of QuEChERS sample preparation methods for the analysis of pesticide residues in fruits and vegetables. J. Chromatogr. A 2010, 1217, 2548–2560. [Google Scholar] [CrossRef]
  37. Anastassiades, M.; Lehotay, S.J.; Štajnbaher, D.; Schenck, F.J. Fast and easy multiresidue method employing acetonitrile extraction/partitioning and “dispersive solid-phase extraction” for the determination of pesticide residues in produce. J. AOAC Int. 2003, 86, 412–431. [Google Scholar] [CrossRef]
  38. Herrmann, S.S.; Poulsen, M.E. Clean-up of cereal extracts for gas chromatography and mass spectrometry. J. Chromatogr. A 2015, 1423, 47–53. [Google Scholar] [CrossRef]
  39. Tuner, J.A. The Pesticide Manual: A World Compendium, 17th ed.; BCPC: Alton, UK, 2015. [Google Scholar]
  40. Korea Crop Protection Association. Crop Protection Agent Guidelines. 2022. Available online: http://www.koreacpa.org/korea/bbs/board.php?bo_table=3_3 (accessed on 30 April 2022).
  41. Joint FAO/WHO Expert Committee on Food Additives (JECFA). Procedures for Recommending Maximum Residue Limits: Residues of Veterinary Drugs in Food; JECFA: Geneva, Switzerland, 2000. [Google Scholar]
Figure 1. Total ion chromatogram (TIC) of a 100 ng/g mixture of 24 pesticides in solvent. The chromatogram displays the separation of 24 pesticide compounds in a solvent using LC–MS/MS, demonstrating their distinct retention times and ionization profiles. The numbers of the peaks correspond to the sequence listed in Table 1.
Figure 1. Total ion chromatogram (TIC) of a 100 ng/g mixture of 24 pesticides in solvent. The chromatogram displays the separation of 24 pesticide compounds in a solvent using LC–MS/MS, demonstrating their distinct retention times and ionization profiles. The numbers of the peaks correspond to the sequence listed in Table 1.
Toxics 12 00633 g001
Figure 2. Analytical method for the extraction and purification of pesticide residues in fishery products. The figure illustrates the procedures and steps involved in the developed method for isolating and purifying pesticide residues from fishery products.
Figure 2. Analytical method for the extraction and purification of pesticide residues in fishery products. The figure illustrates the procedures and steps involved in the developed method for isolating and purifying pesticide residues from fishery products.
Toxics 12 00633 g002
Figure 3. Assessment of dietary exposure to lufenuron in fishery products. (A) Estimated daily intake (EDI) and (B) health risk index (HI).
Figure 3. Assessment of dietary exposure to lufenuron in fishery products. (A) Estimated daily intake (EDI) and (B) health risk index (HI).
Toxics 12 00633 g003
Table 1. Experimental conditions for LC–MS/MS analysis in multiple reaction monitoring mode.
Table 1. Experimental conditions for LC–MS/MS analysis in multiple reaction monitoring mode.
NoPesticideRetention Time
(min)
Precursor Ion
(m/z)
Product Ion
(m/z)
Collision Energy
(eV)
1Acetamiprid5.448223.1126.1−11
55.9−16
2Azinphos-methyl8.218318.0132.0−16
160.1−7
3Azoxystrobin8.04404.0328.9−26
343.9−20
4Carbendazim5.269192.1132.1−25
105.1−37
5Carbofuran6.449222.1165.0−6
123.1−11
6Clothianidin5.274250.0132.0−16
169.1−13
7Dichlorvos6.487238.0109.1−21
220.9−11
8Dinotefuran4.603203.1114.1−13
87.0−16
9Diuron (DCMU)7.699233.046.1−17
−8.039
10Flubendiamide9.605683.0408.0−11
274.0−29
11Hexaconazole10.210314.170.0−22
159.0−30
12Pyrimisulfan7.421420.1370.1−19
388.1−13
13Terbuthylazine8.652230.1174.0−7
68.1−38
14Thiacloprid5.741253.0126.0−11
90.1−39
15Thiamethoxam4.943292.0181.1−22
132.0−22
16Tricyclazole6.006190.1136.0−24
109.0−36
17Azimsulfuron7.42424.9182.1−17
156.0−34
18Trifroxystrobin11.222409.4186.0−20
144.9−45
19Difenoconazole10.816405.9250.9−25
337.0−19
20Fenobucarb8.04208.095.1−14
151.9−10
21Tebufenozide9.454353.0133.0−11
297.1−9
22Lufenuron12.291510.8158.1−21
141.1−42
23Indoxacarb10.823528.1203.0−38
150.1−22
24Daimuron8.717269.2151.2−12
119.2−21
Table 2. Korean dietary consumption of 15 fish species.
Table 2. Korean dietary consumption of 15 fish species.
Fish SpeciesFood Consumption (g/Person/Day) in KNHANES a
(2017–2021)
MeanExtreme (99th Percentile)
Abalone0.600013.8000
Carp0.480012.1200
Chinese muddy loach0.960043.6800
Crucian carp0.480012.1200
Eel1.320024.2400
Far eastern catfish0.480012.1200
Flathead mullet0.480012.1200
Korean rockfish1.200031.2000
Mirror carp0.480012.1200
Olive flounder1.350048.4500
Rainbow trout0.480012.1200
Red seabream0.480012.1200
Sea bass0.480012.1200
Starry flounder0.480012.1200
Whiteleg shrimp1.800050.4000
a Korea National Health and Nutrition Examination Survey.
Table 3. Comparison of pesticide recovery rates for different extraction and clean-up methods.
Table 3. Comparison of pesticide recovery rates for different extraction and clean-up methods.
PesticideRecoveries (Mean ± Standard Deviation) (%)
Extraction MethodClean-Up Materials
EN 15662 aAOAC bC18 50 mgPSA 50 mg +
C18 50 mg
Acetamiprid87.9 ± 4.677.8 ± 6.775.9 ± 0.675.4 ± 1.3
Azinphos-methyl84.0 ± 1.868.3 ± 7.475.4 ± 4.874.2 ± 3.4
Azoxystrobin80.8 ± 1.690.9 ± 4.288.1 ± 2.580.5 ± 2.2
Carbendazim89.7 ± 3.149.9 ± 2.686.2 ± 2.875.0 ± 2.7
Carbofuran69.1 ± 4.270.5 ± 7.781.8 ± 6.580.5 ± 1.9
Clothianidin72.1 ± 3.374.7 ± 6.880.8 ± 1.584.2 ± 6.0
Dichlorvos63.5 ± 6.454.4 ± 7.077.2 ± 3.069.7 ± 2.0
Dinotefuran72.7 ± 1.964.3 ± 6.077.9 ± 3.892.8 ± 3.4
Diuron70.3 ± 4.368.0 ± 5.680.2 ± 0.375.4 ± 2.7
Flubendiamide89.3 ± 4.855.2 ± 1.680.2 ± 2.573.4 ± 4.9
Hexaconazole70.4 ± 5.165.2 ± 5.181.9 ± 1.777.6 ± 2.8
Pyrimisulfan86.5 ± 4.077.9 ± 4.483.6 ± 4.6265.7 ± 26.1
Terbuthylazine76.6 ± 3.557.6 ± 3.281.4 ± 4.076.6 ± 1.9
Thiacloprid81.9 ± 5.975.4 ± 7.280.3 ± 3.776.8 ± 2.9
Thiamethoxam78.7 ± 6.530.0 ± 3.186.2 ± 2.854.8 ± 3.4
Tricyclazole78.4 ± 3.260.8 ± 4.785.4 ± 3.777.5 ± 2.4
Azimsulfuron77.6 ± 0.360.8 ± 3.183.3 ± 0.5−3.5 ± 4.9
Trifroxystrobin74.6 ± 4.259.5 ± 1.582.2 ± 1.880.6 ± 1.9
Difenoconazole76.9 ± 3.848.4 ± 5.073.6 ± 4.579.5 ± 2.0
Fenobucarb76.5 ± 5.974.6 ± 2.886.4 ± 3.180.7 ± 3.4
Tebufenozide76.7 ± 3.790.4 ± 3.476.0 ± 2.776.6 ± 0.5
Lufenuron96.7 ± 5.074.5 ± 3.977.5 ± 5.680.0 ± 3.6
Indoxacarb79.4 ± 2.867.8 ± 3.679.1 ± 3.398.4 ± 3.2
Daimuron80.5 ± 1.736.5 ± 2.287.8 ± 3.480.7 ± 1.3
a 20 mL of acetonitrile containing 0.1% of acetic acid, 4 g of MgSO4, 1 g of NaCl, and 1.5 g of Na3Cit were used in the EN 15662 method for the extraction of pesticides. b 20 mL of acetonitrile containing 0.1% of acetic acid, 6 g of MgSO4, and 1.5 g of sodium acetate were used in the AOAC method for the extraction of pesticides.
Table 4. Linearity, matrix effect, limit of detection (LOD), limit of quantification (LOQ), recovery, and precision of multiclass pesticides.
Table 4. Linearity, matrix effect, limit of detection (LOD), limit of quantification (LOQ), recovery, and precision of multiclass pesticides.
PesticidesMatrixLinearity
(R2)
Matrix
Effect
LOD
(ng/g Wet Weight)
LOQ
(ng/g Wet Weight)
Relative Recovery (%)RSD a (%)
LOQ10 × LOQ50 × LOQLOQ10 × LOQ50 × LOQ
AcetamipridEel0.999916.800.4588.890.192.72.21.42.8
Flatfish0.9993−2.730.4584.478.993.21.59.01.1
Avalone1.0000−5.600.25388.990.187.92.81.12.9
Shrimp0.9997−9.550.6790.892.893.71.51.42.1
Azinphos-methylEel0.999524.870.25389.789.890.712.11.92.5
Flatfish0.99905.850.25384.086.693.16.13.14.6
Avalone0.9989−1.870.3485.191.790.57.32.55.0
Shrimp0.99980.770.4595.490.092.74.12.31.0
AzoxystrobinEel0.999819.440.3489.490.592.12.21.14.1
Flatfish0.9994−0.0020.25383.987.994.65.71.42.1
Avalone0.9998 −4.590.25390.491.491.47.13.32.9
Shrimp0.9995−5.990.3485.893.894.05.71.53.3
CarbendazimEel0.999515.060.3486.785.888.64.53.44.9
Flatfish0.9989−8.800.6780.279.087.74.13.52.1
Avalone0.9999 −11.900.4586.788.385.73.91.94.7
Shrimp0.9992−15.110.25387.691.089.18.43.35.2
CarbofuranEel0.999624.360.25391.888.892.74.71.23.8
Flatfish0.99940.830.25387.287.993.74.21.41.4
Avalone0.9998 −0.310.3493.489.588.45.00.72.9
Shrimp0.9995−4.040.4590.092.394.15.22.12.2
ClothianidinEel0.999423.750.25395.491.093.95.14.24.0
Flatfish0.9992−2.450.3475.283.491.011.21.92.9
Avalone0.9997 −6.900.3485.888.685.910.53.72.4
Shrimp0.9983−6.530.4584.292.591.52.82.72.0
DichlorvosEel0.999615.870.4584.088.192.77.12.11.7
Flatfish0.9988−4.160.4582.785.193.78.51.53.2
Avalone0.9999 −6.480.81088.990.891.45.62.42.7
Shrimp0.9993−9.030.81071.789.792.69.91.62.3
DinotefuranEel0.99985.450.3487.188.991.45.31.61.6
Flatfish0.9990−1.800.3482.184.890.23.92.51.6
Avalone0.9999 −9.930.4589.288.787.16.51.33.0
Shrimp0.999625.450.4588.189.591.91.62.72.0
DiuronEel0.999627.790.4591.589.092.04.31.61.6
Flatfish0.9991−1.960.25384.883.592.05.71.80.5
Avalone0.9999 −3.490.3493.289.389.63.61.92.1
Shrimp0.9991−3.490.3493.289.389.63.61.92.1
FlubendiamideEel0.999121.290.25388.387.388.56.22.75.9
Flatfish0.99767.120.4492.888.895.15.63.13.8
Avalone0.9995 2.040.4486.293.490.78.11.54.3
Shrimp0.99650.890.25395.392.893.710.74.14.9
HexaconazoleEel0.9921−1.770.25388.393.997.14.91.23.8
Flatfish0.9988−2.580.25383.784.989.35.83.22.9
Avalone0.9996 −8.250.25377.690.587.913.22.33.1
Shrimp0.9990−7.470.25382.991.492.95.82.95.2
PyrimisulfanEel0.999813.620.4590.389.092.21.91.23.8
Flatfish0.99906.350.3487.686.693.42.01.21.1
Avalone0.9996 1.080.25390.290.289.22.71.12.5
Shrimp0.99973.160.3488.690.992.81.21.82.1
TerbuthylazineEel0.999528.460.4588.886.790.12.60.63.7
Flatfish0.9994−1.350.25387.085.790.95.11.81.5
Avalone0.9999 −5.930.3490.191.689.74.41.32.9
Shrimp0.9994−4.810.3482.492.191.56.23.53.8
ThiaclopridEel0.999920.480.4584.890.391.92.81.24.2
Flatfish0.99950.820.4585.086.092.52.31.91.1
Avalone0.9939 −2.120.3488.489.986.72.81.70.4
Shrimp0.9994−6.370.6786.193.093.02.41.31.4
ThiamethoxamEel0.999716.730.25393.490.293.97.71.12.7
Flatfish0.9995−11.660.6781.282.591.92.82.61.8
Avalone1.0000 −14.100.3490.689.087.55.11.12.3
Shrimp0.9991−6.820.4585.092.194.24.52.32.4
TricyclazoleEel0.9997−0.080.25390.190.7101.33.30.82.4
Flatfish0.9993−20.330.4586.387.3106.31.02.01.9
Avalone0.9998 −19.230.4585.688.993.95.10.82.2
Shrimp0.9979−26.810.3490.695.4101.03.61.11.7
AzimsulfuronEel0.999815.920.3487.588.291.13.11.52.2
Flatfish0.99936.150.3483.083.490.53.61.71.9
Avalone0.9997 1.800.3485.890.089.23.72.12.1
Shrimp0.99974.490.25385.188.290.93.81.72.8
TrifroxystrobinEel0.999027.810.4596.890.992.74.53.03.8
Flatfish0.9984−3.560.25390.490.595.82.72.31.8
Avalone0.9998 −3.160.25387.890.588.74.32.12.1
Shrimp0.9986−10.760.25387.194.6100.57.82.11.5
DifenoconazoleEel0.999316.220.3486.388.792.51.81.42.2
Flatfish0.99901.250.3483.086.191.93.40.71.0
Avalone1.0000 −5.080.25387.691.890.55.80.82.5
Shrimp0.9996−4.780.25386.392.391.74.81.52.4
FenobucarbEel0.999718.710.3487.688.992.23.91.82.2
Flatfish0.9990−0.450.3481.786.091.77.41.22.3
Avalone0.9997 −4.370.3488.792.289.15.31.82.9
Shrimp0.9994−2.130.25390.193.793.13.11.72.0
TebufenozideEel0.999822.510.25392.190.693.92.81.83.9
Flatfish0.99930.710.4583.186.894.43.61.71.6
Avalone0.9998 −3.820.3484.291.589.04.61.71.9
Shrimp0.9989−3.610.4590.591.693.83.72.14.3
LufenuronEel0.999320.250.4587.686.687.25.32.46.7
Flatfish0.9997−10.190.4578.891.991.612.07.03.1
Avalone0.9997 −17.850.3480.491.984.114.22.26.9
Shrimp0.9992−3.160.25393.794.095.012.92.54.9
IndoxacarbEel0.999715.400.4590.889.292.62.22.45.1
Flatfish0.99765.460.3481.688.793.712.52.62.7
Avalone0.9997 −1.230.4594.190.990.57.72.01.2
Shrimp0.9988−2.550.3489.290.892.87.22.36.6
DaimuronEel0.999816.930.3487.288.991.41.31.53.6
Flatfish0.99931.610.3487.286.793.11.70.80.9
Avalone0.9999 −2.360.25388.791.490.41.60.92.3
Shrimp0.9997−3.470.25389.191.693.33.01.62.7
a Relative standard deviation of repeatability (intra-day tests).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kim, M.; Kim, T.-h.; Park, J.-W.; Lee, Y.; Jo, M.-R.; Moon, Y.-S.; Im, M.-H. A Robust Method for Simultaneous Determination and Risk Assessment of Multiresidual Pesticides in Fishery Products. Toxics 2024, 12, 633. https://doi.org/10.3390/toxics12090633

AMA Style

Kim M, Kim T-h, Park J-W, Lee Y, Jo M-R, Moon Y-S, Im M-H. A Robust Method for Simultaneous Determination and Risk Assessment of Multiresidual Pesticides in Fishery Products. Toxics. 2024; 12(9):633. https://doi.org/10.3390/toxics12090633

Chicago/Turabian Style

Kim, Myungheon, Tae-hwa Kim, Jong-Woo Park, Yoonmi Lee, Mi-Ra Jo, Yong-Sun Moon, and Moo-Hyeog Im. 2024. "A Robust Method for Simultaneous Determination and Risk Assessment of Multiresidual Pesticides in Fishery Products" Toxics 12, no. 9: 633. https://doi.org/10.3390/toxics12090633

APA Style

Kim, M., Kim, T. -h., Park, J. -W., Lee, Y., Jo, M. -R., Moon, Y. -S., & Im, M. -H. (2024). A Robust Method for Simultaneous Determination and Risk Assessment of Multiresidual Pesticides in Fishery Products. Toxics, 12(9), 633. https://doi.org/10.3390/toxics12090633

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