Next Article in Journal
Removal of Antibiotics and Nutrients by Vetiver Grass (Chrysopogon zizanioides) from a Plug Flow Reactor Based Constructed Wetland Model
Previous Article in Journal
Comparison of Simultaneous Quantitative Analysis of Methylmercury and Inorganic Mercury in Cord Blood Using LC-ICP-MS and LC-CVAFS: The Pilot Study of the Japan Environment and Children’s Study
Previous Article in Special Issue
A Toolbox for the Determination of Nitroaromatic Explosives in Marine Water, Sediment, and Biota Samples on Femtogram Levels by GC-MS/MS
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Toxics
Volume 9
Issue 4
10.3390/toxics9040083
toxics-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Easy Procedure to Quantify Anticoagulant Rodenticides and Pharmaceutical Active Compounds in Soils

by
Andrea Acosta-Dacal
1,†,
Cristian Rial-Berriel
1,†,
Ricardo Díaz-Díaz
2,
María del Mar Bernal-Suárez
2,
Manuel Zumbado
1,3,
Luis Alberto Henríquez-Hernández
1,3 and
Octavio P. Luzardo
1,3,*
1
Toxicology Unit, Research Institute of Biomedical and Health Sciences (IUIBS), Universidad de Las Palmas de Gran Canaria, Paseo Blas Cabrera s/n, 35016 Las Palmas de Gran Canaria, Spain
2
Department of Environmental Analysis, Technological Institute of the Canary Islands, C/Los Cactus No 68, Polígono Industrial de Arinaga, Agüimes, 35118 Las Palmas, Spain
3
Spanish Biomedical Research Center in Physiopathology of Obesity and Nutrition (CIBERObn), Instituto de Salud Carlos III, 28029 Madrid, Spain
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to the work, and therefore should be considered indistinctly as first authors.
Toxics 2021, 9(4), 83; https://0-doi-org.brum.beds.ac.uk/10.3390/toxics9040083
Submission received: 8 March 2021 / Revised: 2 April 2021 / Accepted: 7 April 2021 / Published: 10 April 2021
(This article belongs to the Special Issue New and Novel Analytical Approaches to Measurement of Emerging Environmental Contaminants)
Browse Figures
Versions Notes

Abstract

:
A modified QuEChERS (quick, easy, cheap, effective, rugged, and safe) extraction was validated for the extraction of seven coumarin anticoagulant rodenticides (ARs) and 36 pharmaceutical active compounds (PhACs) residues in soil samples using liquid chromatography tandem mass spectrometry (LC-MS/MS). The aim of this work was to develop a method for the monitoring of these compounds in agricultural lands as well as in forensic applications for the determination of ARs poisoning. As far as is known, this is the first time that a QuEChERS-based method is used for the extraction of ARs in soil, as well as on such a quantity of PhACs. A matrix effect study was carried out on samples of soil devoted to agriculture in the midland area of the Canary Islands (clay loam type). It was in house validated (accuracy, precision, and linearity) at seven spiked levels between 0.5 and 50 ng g−1. The limits of quantification (LOQ) ranged between 0.5 and 50.0 ng g−1 and the limits of detection (LOD) ranged from 0.024 to 6.25 ng g−1. The method was then successfully used for both the determination of the target analytes in the soils from the agricultural plots that had been irrigated with regenerated water, and in the soil collected from underneath wild bird carcasses (which had been the subject of forensic investigation).

1. Introduction

Soil is a very vulnerable pollution receptor environment from where pollutants can be emitted to other environmental compartments such as the atmosphere, ground or surface water, and biota [1]. Contaminated soil presents a large number of problems that may include the destruction of ecosystems, agricultural productivity losses, the contamination of groundwater, and danger to human and animal life due to the accidental ingestion of soil or consumption of food that has grown on contaminated soil. Agricultural soil contamination can be direct and intended, thorough the application of pesticides to crops, or indirect and nonintended, mainly due to irrigation using regenerated waters and the use of compost of wastewater treatment plants (WWTPs) sludge. Both water and sewage sludge may contain various contaminants of emerging concern (CECs), including anticoagulant rodenticides (ARs) and pharmaceutical active compounds (PhACs).
ARs are the most widely used rodent control agent in agriculture (usually in postharvest plants and warehouses), livestock farming, and in domestic and urban applications [2]. These substances act by preventing the formation of clotting factors by blocking the vitamin K cycle and, thus, causing spontaneous internal bleeding [3]. The first generation of ARs (FGARs) started with warfarin, soon followed by coumatetralyl, chlorophacinone, and diphacinone [4]. Following the widespread development of rodent resistance to FGARs, a second and more toxic generation of compounds (SGARs), including flocoumafen, difethialone, difenacoum, bromadiolone and brodifacoum, were developed in the 1970s [5]. On the other hand, some anticoagulants, such as warfarin or its analogues, are used in the medical treatment of blood hypercoagulability disorders [6]. Taking this into account, ARs can reach wastewater through runoff after agricultural and livestock farming applications, bait solubilization in urban infrastructures or through urinary excretion from medical treatments [7,8].
PhACs are essential for the prevention and treatment of disease in humans and animals [9,10]. After administration, these active substances are absorbed and can be either metabolized or converted back into active ingredients [11]. A part of them is excreted through urine or feces, and may enter WWTPs and end up in environmental compartments as both metabolites and parent compounds [12,13]. The presence of PhACs in the environment can favor the appearance and the development of various adverse effects, such as antibiotic-resistant microorganisms, even at low levels due to continuous exposure to antimicrobials [13].
Once in WWTPs, both ARs and PhACs can be degraded, accumulate in sewage sludge or remain in the water, even after tertiary treatments, which are released to effluent receiving waters or used to irrigate crops [14]. In addition, due to their hydrophobicity, ARs are more likely to accumulate in sewage sludge, which is commonly used to make compost [15] while PhACs can be found in both media depending on their physicochemical properties [16]. Similarly, ARs are relatively stable and have a moderate to high octanol/water partition coefficient (KOW) allowing them to associate with natural soil organic matter [17,18]. As for the fate of PhACs, depending on their interaction with soil matrix components, they may remain in or be transported to groundwater and surface water through surface runoff, infiltration/percolation [12,19]. As seen, the amounts of ARs and PhACs that may eventually accumulate in soils through irrigation with regenerated wastewater and/or through sludge and manure spreading or from direct bait application on fields, should be a matter of concern.
Both PhACs and ARs may pose a high risk to wildlife due to their potential for bioaccumulation and transfer through trophic webs. Although ARs are intended to control rodents, baits may be consumed by other nontarget species such as granivorous birds (primary exposure). Similarly, predators or scavengers may consume rodents and nontarget animals that have already been exposed to ARs (secondary exposure) [20]. Moreover, ARs baits are often intentionally used to poison and kill not rodents but other animals. Likewise, the exposure of wild animals to PhACs can occur in several ways: contaminated water [21], agricultural soils, plants and arthropods [22,23]. The carcasses of farm animals that have been previously medicated, as well as their slurry, can also be a relevant source of exposure to pollutants for wildlife [16,24]. In such cases, exposure is usually monitored by blood or liver tissue samples from dead animals [25,26]. However, it is quite usual to find only the remains of animal carcasses and not being able to collect tissue or fluid samples. Therefore, for these specific cases, an extraction method to test the soil sampled under the carcass may be an appropriate way to prove exposure to these compounds due to the capacity of the soil organic matter to retain them after being leached from the carcasses.
There are very few studies on extraction methods for anticoagulants in soils [27,28,29,30,31,32] and, to our knowledge, none of them use the QuEChERS methodology. Although there are more methods of extracting PhACs from soils, only a small number of them use QuEChERS procedures [33,34,35,36]. The QuEChERS method consists of an acetonitrile extraction/partitioning step followed by a solid-phase dispersive extraction as a clean-up step and was first developed for the extraction of pesticides from fruits and vegetables [37]. This user-friendly methodology provides high extraction yields using reduced amounts of samples and organic solvents, which has led to its use in different matrices and analytes [38]. Modified versions of QuEChERS have been successfully applied in soils for the extraction of organic contaminants, including pesticides [39,40,41], persistent organic pollutants [42,43], as well as in the extraction of ARs and PhACs from other samples such as food [44,45], blood [46,47], sediments and sludge [48], and the liver [49,50].
Along with a robust extraction method, adequate detection analysis is necessary. ARs and PhACs have been determined mainly by liquid chromatography (LC) with different detectors such as UV [32,51], diode array (DAD) [30,52], DAD with fluorescence detectors (FLD) [53,54] and single quadruple mass spectrometry (MS) detectors [28,55]. However, the preferred option when looking for sensitivity and selectivity at trace level is liquid chromatography-tandem mass spectrometry (LC-MS/MS) [7,35,56,57], which is the best option for low-level residues that can be expected to be found in soils.
The objective of this research was to evaluate and validate a QuEChERS-based and LC-MS/MS method for the quantification of ARs and PhACs in soils, as well as their verification in samples of agricultural origin and from wildlife carcasses. The work is proposed as a scope extension of the previously optimized method for the quantification of pesticide residues in soils [41]. The proposed method can be applied in agricultural, environmental and forensic monitoring.

2. Materials and Methods

2.1. Reagents and Chemicals

Individual certified standards of ARs and PhACs (purity 95.19% to 99.9%) were acquired from Dr. Ehrenstorfer (Augsburg, Germany), Sigma-Aldrich (Augsburg, Germany) and European Pharmacopoeia Reference Standards (Strasbourg, France). Atrazine-d5, Carbendazim-d3, Coumachlor, Cyromazine-d4, Linuron-d3 and Pirimicarb-d6 (Dr. Ehrenstorfer and Sigma-Aldrich, 99.3–99.9% purity) were used as procedural internal standards (P-IS) and were maintained as described in the reference method [41]. These internal standards were added to the samples at the beginning of the procedure to account for various sources of errors throughout all stages in the method [58]. The list of the selected analytes, their classification, and their legal status are in Table 1.
LC-MS grade methanol (MeOH), acetonitrile (ACN) and formic acid (FA, HCOOH) were obtained from Honeywell (Morristown, NJ, USA). Ammonium acetate (NH4CH3CO2) was purchased from Fisher Scientific (Loughborough, UK). AOAC method QuEChERS salts [59] (6 g of MgSO4 and 1.5 g of CH3COONa) were acquired in commercial premixes from Agilent Technologies (Palo Alto, CA, USA). The ultrapure water was produced in the laboratory using a Gradient A10 Milli-Q System (Millipore, Bedfore, MA, USA).

2.2. Standard Stock Solutions and Mixes

Individual standard stock solutions of all analytes and P-IS were prepared in ACN at a concentration of 1000 μg mL−1. Individual intermediate solutions of 1 μg mL−1 were prepared for spectrometry optimization. In addition, two mixed stock solutions were prepared, one containing target ARs and the other containing the selected PhACs, at 10 μg mL−1 followed by a mixed working solution diluted ten times more than 1 μg mL−1 of each. Finally, a P-IS solution was prepared at 1 μg mL−1.
Both solvent and matrix-matched calibration curves were prepared with the standard working mix solution in either soil extracted with the evaluated procedure or ACN 2.5% FA diluted with ultrapure water (1:1, v/v).
All standards, working mix solutions and matrix-matched calibrators were stored in glass amber vials at −20 °C and checked periodically for stability.

2.3. Sample Selection and Pretreatment

The extraction method was validated on a representative type of soil, which due to its physicochemical properties, can be classified as clay loam soil. The soil was chosen from two farms dedicated to organic production and sampled at various times throughout 2020. A composite sample of at least four subsamples collected at depths between 20 and 30 cm was prepared in each sampling plot. It was then thoroughly mixed and homogenized, air-dried at room temperature and sieved (2 mm mesh) before being considered suitable for analysis.
The physicochemical properties of this composite are as follows: pH 4.88, electrical conductivity 209 μS cm−1, oxidizable organic carbon 2.19% (approximately 3.9% organic matter), 6% moisture and particle size distribution: 29.5% clay, 28.3% fine silt, 11.3% thick slit, 11.5% coarse sand, and 19.4% fine sand. All these parameters were determined in the facilities of the Department of Department of Animal Biology, Edaphology and Geology of the University of La Laguna (Tenerife, Spain) [41].

2.4. Sample Preparation

The method used for the extraction of the target analytes was previously developed and validated in our laboratory for the extraction of pesticides in clay loam soil [41,60]. In summary, the samples were prepared by weighing 10 ± 0.05 g of dried and sieved soil into a 50 mL centrifuge. All samples were spiked with the appropriate volume of the P-IS solution to achieve a concentration of 5 ng g−1, and recoveries and Quality Controls (QCs) were added standard mix solutions and left to stand for 1 h prior to extraction. Then, 10 mL of ACN-2.5% FA were added and vigorously shaken for 1 min. In the same way, 6 g of MgSO4 and 1.5 g of CH3COONa were added, energetically shaken for another minute and sonicated for 15 min in an ultrasonic bath operating at 50/60 Hz and 120 W (VWR, Radnor, PA, USA). After that, samples were placed 25 min in a rotatory shaker (Ovan, Barcelona, Spain). They were then centrifuged for 10 min at 4200 rpm (3175.16× g) in a 5804 R Eppendorf centrifuge (Eppendorf, Hamburg, Germany). An aliquot of the supernatant extract was filtered through 0.20 µm Chromafil® PET filters (Macherey-Nagel, Düren, Germany) and diluted with H2O (1:1, v/v) before analysis in LC-MS/MS.

2.5. LC-MS/MS Analysis

LC-MS/MS equipment (1290 model UPLC tandem coupled to a 6460 model Triple Quadrupole mass spectrometer, Agilent Technologies, Palo Alto, CA, USA). The chromatographic column was an Agilent Poroshell 120 EC-C18 (2.1 × 100 mm, 2.7 µm), with a guard precolumn (2.1 × 5 mm, 1.8 µm) plus a and prefilter. The column oven was set at 50 °C. The mobile phases consisted of 2 mM ammonium acetate 0.1% formic acid in ultrapure water (A) and 2 mM ammonium acetate in MEOH (B). The binary gradient was set as follows: 5% B–0.5 min; 5% B–1 min; 40% B–2.5 min; 85% B–8 min; 100% B–10 to 14 min; 5% B–14.01 min. The flow rate was 0.4 mL min−1. The volume injected was 5 μL and the total run time was 18 min.
The parameter for the MS/MS analyzer were those described in the original method [41]. Data analysis was performed using Agilent software MassHunter Quantitative Analysis (for QQQ) version B.07.01 and MassHunter Qualitative Analysis version B.07.00.
The MS/MS conditions were optimized by injecting 5 μL of individual solutions of each compound at 1 μg mL−1 in ACN directly to the mass spectrometer replacing the chromatographic column with a stainless steel zero dead volume union. The mobile phases were set in isocratic mode (50:50, v/v), and these were 2 mM ammonium acetate and 0.1% FA in ultrapure water (phase A) and 2 mM ammonium acetate and 0.1% FA in methanol (B). The product ions were optimized in the MRM mode at different collision energies, from which that exhibiting the highest response was selected.

2.6. In House Validation Parameters

To validate the use of the proposed method and broaden it to the selected ARs and PhACs, we performed a single-laboratory validation. Since there is no a specific guidance for the analysis of ARs and PhACs residues in soil, the in house validation (hereafter validation) of the proposed method was performed following the guidelines for the analysis of pesticide residues in food and feed of the European Union [58,61].
The linear range of LC-MS/MS method was studied using calibration curves prepared in soil matrix-water (1:1, v/v), ranging from 0.024 to 50 ng g−1 using both correlation coefficient and Mandel test (95% confidence level) [62]. Accuracy was estimated by recovery assays (in quintuplicate) through the quantification of trueness and precision (% relative standard deviation) at 7 concentration levels: 0.5, 1, 2.5, 5, 10, 20 and 50 ng g−1. According to the SANTE and SANCO guides, recoveries between 70–120% are considered acceptable, when the RSDs are below 20%. Additionally, following the criteria specified in those guides the limit of quantification (LOQ) for each compound was the lowest level of the recovery experiments that met all the validation criteria: recovery between 80–120%, and RSD below 20%. On the other hand, the limits of detection (LOD) were calculated using calibration standards. For this purpose, triplicate matrix-matched calibration curves were prepared ranging from 0.024 to 100 ng g−1. Thus, the LOD was selected as the lowest point of the calibration curve that meets had a signal-to-noise ratio (S/N) > 3 (Peak-to-Peak algorithm) and an accuracy between 80–120%.
The identity of the compounds was based on the acquisition of two transitions in the MRM mode, one of them employed as the quantification transition (Q), and the other as the confirmation transition (q). According to the SANTE and SANCO guides an ion ratio tolerance of 30% was considered acceptable. The tolerance for the deviation of the retention time (tR) ± 0.1 with respect to that of the reference standard.

3. Results and Discussion

3.1. Optimization of LC-MS/MS Conditions

The seven target rodenticides, chloramphenicol, tolfenamic acid, and the P-IS coumachlor were determined in negative mode with precursor ions corresponding to [M − H]−. The remaining analytes and P-IS were analyzed in positive mode using precursor ions corresponding to [M + H]+ except for cefuroxime axetil, cloxacillin, eprinomectin, josamycin and penicillin V. It is known that, together with the protonated and deprotonated molecules already mentioned, the ionization process can produce a variety of adduct ions depending on the composition and concentration of mobile-phase additives and analytes, the pH, the mobile-phase flow rate, the temperature or even solvent impurities and glassware [63,64,65]. While cefuroxime axetil formed a sodium adduct [M + Na]+, cloxacillin, josamycin and penicillin V formed methanol adducts, which was the organic phase solvent. In the case of eprinomectin, the ion m/z 878 selected as parental for both transitions is its protonated derivative ion resulting from dehydrative aromatization [66]. In the case of brodifacoum and bromadiolone, which have Br in their chemical structures, an element with an abundant characteristic isotope distribution, transitions corresponding to 79Br and 81Br were selected. In the same way, chloramphenicol, which has two atoms of Cl in its structure, was optimized selecting the transitions corresponding to 35Cl and 37Cl. The transitions were chosen in terms of selectivity and sensitivity, selecting the most abundant ones for quantification proposals.
The parameters of the ion source were then optimized injecting standards at 100 ng mL−1 in ACN (in triplicate). The parameters were nebulizer and sheath gas flow and temperature, and capillary and nozzle voltages. The Agilent software (Mass Hunter Source Optimizer) was employed for this purpose.
Likewise, we also injected 20 ng mL−1 standards in triplicate, prepared in the same solvent, for the chromatographic condition optimization process. The final conditions selected for both the source and the chromatography were then tested with a matrix-matched standard at the same concentration. We decided to use Poroshell 120 EC-C18 column (2.1 × 100 mm, 2.7 µm; Agilent Technologies) due to the satisfactory performance for the variety of pesticide belonging to different chemical groups optimized in the previous work [41]. The mobile phase solvents and modifiers used for the compound optimization process were also set as the final phase composition to ensure the formation of the mentioned adducts. However, different concentrations and gradients of ammonium acetate (2 and 5 mM) and FA (none and 0.1%) were tested to achieve the best resolution and separation of the compounds. A higher amount of ammonium not only did not improve the sensitivity, but also increased the pressure in the column, so it was set as 2 mM. However, we eliminated FA from the organic phase as it worsens the sensitivity of ARs compounds, but we continued to use it in the aqueous phase because it was necessary for the analysis of some PhACs. This decision was made because ARs retention starts with warfarin at 7.87 min (with this combination of mobile phases and gradient) where the organic mobile phase is approximately 83% and continues to increase. This way, the acid composition in the column will decrease when the ARs are retained in the column.
The injection volume was set at 5 μL after testing from 1 to 20 μL because it produced good peak shapes while giving high responses. In addition, a higher volume would result in saturation of the MS/MS detector. The possible dilution of the extract with water was also evaluated at this point before injecting it in the LC-MS/MS equipment. The ratios tested were 1:1, 1:2, 1:5 y 1:10 (v/v) of 20 ng g−1 standards in ACN-2.5%FA-ultrapure water. The experiment was then repeated with standards in the soil matrix at the same dilutions. In order to achieve considerable sensitivity and selectivity, while improving the peak shape of some compounds such as chloramphenicol, dexamethasone and mebendazole, we chose to dilute the final extract 1:1, v/v with water.
All the chromatography and mass spectrometry conditions are shown in Table 1. In addition, a chromatogram of a blank soil sample spiked at 50 ng g−1 is shown in Figure 1.

3.2. Matrix Effect Study

Soil is an extremely complex matrix whose components can suppress or enhance the response of the target analytes in the mass spectrometer and can therefore condition the integrity of the analysis. Bearing this in mind, a study of the matrix effect was performed. The matrix effect (ME) was evaluated comparing the slopes of calibration curves in the solvent (SS) and in the matrix (Sm), which was extracted using the procedure described in Section 2.4, according to the equation:
ME (%) = (Sm/SS) × 100
Hence, the effect of the matrix components on the signal is qualified as the percentage of enhancement or suppression, whether the ME values are above or below 100%, respectively. No significant matrix effects were considered when ME was between 80–120% [61].
The calibration curves covered the range of 3.125 to 50 ng g−1 and were prepared either in soil matrix or ACN 2.5%FA, both diluted with ultrapure water (1:1, v/v). All curves were prepared in triplicate and adjusted to a linearity equation (y = ax + b). Figure 2 shows ME mean values and SD together with the tolerance interval in which no matrix effect is considered (in grey). As can be seen in the graph, the majority of the analytes did not present a significant matrix effect. However, moxidectin, sulfanilamide, and sulfapyridine showed strong, medium, and slightly significant signal suppression, respectively, and only tolfenamic acid showed strong signal enhancement. Although cloxacillin and imipenem are in the tolerance range of ME%, they presented a variable value with RSD% over 20%. Consequently, matrix-matched calibration was used for quantification in the subsequent experiments.

3.3. In House Validation Studies

The proposed extraction method was validated following the agreements stated in the “In house validation parameters” section.
Linearity response was satisfactory in the studied range for each compound, with R2 values higher than 0.99 and p-values over 0.05 after Mandel test application (Table 2). All ARs and the majority of the PhACs were satisfactorily extracted in terms of accuracy and precision from their LOQ to the highest concentration tested, with recoveries between 70–120% and RSD values below 20%, respectively. However, cloxacillin and tolfenamic acid showed recoveries above 120% for some of the concentration levels. Similarly, albendazole, dexamethasone, flunixin, oxfendazole, penicillin V, sulfacetamide, sulfametoxipiridacine, sulfamonomethoxine, sulfanilamide and sulfapyridine were recovered below 70% at some of these levels. However, all of these compounds were extracted with high reproducibility (RSD < 20%) in the range of 60–130%, which is contemplated in the SANTE guidelines for routine analysis [58]. On the other hand, other PhACs showed high variability (RSD > 20%) at 5 ng g−1 or below but were recovered between 70–120% which is also taken into account in these guidelines for concentration below 10 ng g−1 of residues in soils. The results of the recovery experiments are summarized in Table 3.
All analytes LOQs were equal to or below 20 ng g−1 except for imipenem with 50 ng g−1, which is the LOQ normally required for residues in soils according to the guidance for analytical methods for residues in soil [61]. As for the LOD, all analytes had a detection limit of 6.25 ng g−1 or even lower (Table 2). As mentioned above, there are very few studies that determine ARs in soil and none of them analyze difethialone residues. Moreover, the LOQ and LOD values obtained with our modified QuEChERS method are much lower than those obtained using other extraction methodologies [27,29,31] except for brodifacoum, coumatetralyl, difenacoum and flocoumafen, whose LOQs are in the same range than those obtained with a MeOH and ammonium format extraction [28]. However, it should be noted that the LOD and LOQ in those cases were established as three and ten times the S/N of the blank sample extract, respectively, while ours were calculated following the more restrictive criteria of the European Union SANTE guidance [58]. Furthermore, among the few QuEChERS-based extraction methods for PhACs in the literature for their determination in soils, our method is capable of extracting and analyzing more of these compounds. Of the compounds that coincide with those determined in the previously published methods, sulfanilamide, sulfadiazine, and sulfadimethoxine have lower LOQs in the method of Salvia et al. [35] but were determined as S/N ratio equal to 10. Similarly, Lee et al. also determined the LOQs as ten times the S/N ratio, and obtained a LOQ for sulfamethazine that was equal to ours (0.5 ng g−1), but the one for sulfamethoxazole was lower (0.5 vs. 2.5 ng g−1) [33]. Finally, Ferhi et al. determined the LOQs with the results of the validation data in a more similar way to ours for diclofenac, sulfamethoxazole, and sulfamethoxazole, our values being lower in the first two cases and slightly higher but in the same range in the latter (2.1 ng g−1 vs. 5 ng g−1) [34].
In summary, the proposed QuEChERS-based method without a clean-up step was expanded from 218 [41] to 261 analytes, and now proves to be accurate and reliable for the analysis of the ARs and selected PhACs in soil samples, in addition to agricultural pesticides.

3.4. Verifying the Method in Different Scenarios

Once the method was validated under the above conditions, it was applied for the analysis of ARs and PhACs in two different scenarios to verify its applicability: Agricultural soils samples irrigated with regenerated water and a soil sample from an environmental forensic investigation of a suspected wildlife poisoning episode.

3.4.1. Agricultural Soil Samples

We decided to analyze samples from mid-range farms on the island of Gran Canaria that use regenerated water from the WWTP of the nearby city (Las Palmas de Gran Canaria, 425,000 inhabitants). In the initial analysis of the soil samples, no PhACs residues were detected in any of them, although several pesticide residues from those included in the initial method were detected in the different samples tested (data not shown, see Acosta-Dacal et al., 2021). In addition, thanks to the extension of the analytical scope, brodifacoum was detected in one of the soil samples, and flocoumafen in another (Figure 3, panels B and C), at concentrations of 2.96 and 1.37 ng g−1, respectively.
Although these soil samples had been irrigated with regenerated wastewater, no PhACs residues were detected, so there was a possibility that the water was free of residues. In order to verify this, samples of the irrigation water from the WWTP were collected during four consecutive days. These water samples were directly analyzed in the LC-MS/MS where three authorized PhACs were detected: diclofenac, ketoprofen, and sulfamethoxazole (panel A of Figure 3). These PhACs are used for both humans and livestock. In contrast, no AR was found in the water samples. Therefore, the presence of brodifacoum and flocoumafen detected in the soil samples could have been due to a rodent-control application. It should be noted that the water analyzed, although it comes from the same treatment plant, is not the same as the water used in the irrigated plots, since it was subsequently collected. In other words, it is also possible that the water used in the irrigated plots at that time contained the rodenticides and not the PhACs. This possibility is also plausible taking into account that neither brodifacoum nor flocoumafen are EU nonauthorized pesticides as agricultural products (Table 1).
As mentioned, ARs tend to accumulate in organic matter due to their moderate to high KOW, so it is to be expected that most of them are eliminated in the active sludge treatment. This is consistent with what other authors have found, such as Gomez-Canela et al. who determined very few ARs above the detection limit in the effluents of the different treatments of the WWTPs analyzed. Nonetheless, they found flocoumafen in the effluent of the primary (18.1 ng L−1) and secondary (18.1 ng L−1) treatments from one of the WWTPs and both flocoumafen (29.3 ng L−1) and brodifacoum (38.4 ng L−1) in the effluent of the tertiary treatment of another one [8]. Therefore, we cannot discard the possibility that although we have not found them in the analyzed water, these compounds could have been present in previous days. On the other hand, the fact that the PhACs detected in water did not show up in the analyzed soil samples could have been due to the low concentrations found, which are very close to the method LOD. The mean concentrations of diclofenac, ketoprofen, and sulfamethoxazole were 1.03, 0.88, and 0.35 ng g−1, respectively. These were either below the limit of quantification we reached with our QuEChERS-based extraction method (2.5 ng g−1 for both diclofenac and sulfamethoxazole) or slightly higher in the case of ketoprofen (0.5 ng g−1). As already mentioned, wastewater content and values may vary on different days even if our samples were in the same range for the period studied.
Thus, we decided to investigate whether the soil could accumulate these compounds during irrigation. First, we prepared a pool with the water samples collected during four consecutive days and concentrated it ten times with the intention of achieving a concentration much higher that the LOQ of the three compounds detected. This was performed using a vacuum concentrator RVC 2-25 CD plus (Christ, Germany) at 35 °C and then reconstituting the sample in a volume 10 times lower of ultrapure water. After that, an irrigation simulation model was prepared with the soil used for the validation procedure. Approximately 500 g of soil was placed in a tray container and irrigated with 0.5 L of the concentrated water solution on alternate days for two weeks (Figure 4, inset). Then, the soil was left to air dry for two weeks, homogenized and extracted using the proposed method. Finally, both the concentrated water and the extract were analyzed. Diclofenac ketoprofen and sulfamethoxazole were extracted with an extraction efficiency of 93.4, 85.1, and 107.2%. Therefore, the soil is capable of accumulating these compounds in irrigation if the concentration is higher. A chromatogram of the irrigated soil together with an image of the experiment can be seen in Figure 4.

3.4.2. Environmental Forensic Investigation

Finally, we tested the validity of the method in another scenario within the context of a judicial forensic investigation. We had dealt with soil samples from a protected natural area on the island of Gran Canaria (Canary Islands, Spain) where the skeletonized remains of several protected birds (ancient corpses) were found. There were strong suspicions that a bait containing some type of poison had been maliciously placed. Therefore, several samples collected in the area were submitted, including a soil sample taken from underneath one of the corpses of a medium-sized bird of prey. This soil sample was analyzed with the method proposed in this article, and as shown in Figure 5, three ARs were identified: brodifacoum, bromadiolone and difenacoum at 223.62, 12.61 and 1.94 ng g−1, respectively.

4. Conclusions

The proposed methodology was successfully in house validated for the determination of 7 coumarin ARs and 36 PhACs in terms of linearity, trueness, and precision for all the analytes. The LODs were between 0.024 and 6.25 ng g−1 and LOQs in the range of 0.5–20.0 ng g−1, and only the imipenem presented the typically fixed LOQ for soil residues, 50 ng g−1. It was then applied to samples from agricultural plots irrigated with regenerated water, and brodifacoum and flocoumafen were found in two of them. It was also used in a soil sample collected from underneath wildlife carcasses in the context of an environmental forensic investigation where three ARs (brodifacoum, bromadiolone, and difenacoum) were identified.
The scope of this simple method has been extended and now allows the determination of 43 CECs in addition to the 218 pesticides in soil. Thus, it can be used for monitoring agricultural soils and for forensic purposes in soils found under wildlife carcasses. To our knowledge, this is the first time that a QuEChERS methodology is applied to ARs extraction in soil matrix.

Author Contributions

Conceptualization, A.A.-D. and O.P.L.; methodology, A.A.-D., C.R.-B. and O.P.L.; software, L.A.H.-H.; validation, A.A.-D., R.D.-D. and M.d.M.B.-S.; formal analysis, A.A.-D., C.R.-B., and M.Z.; investigation, A.A.-D., C.R.-B., and M.Z.; resources, O.P.L., and M.Z.; data curation, L.A.H.-H.; writing—original draft preparation, A.A.-D.; writing—review and editing, O.P.L., R.D.-D., and M.d.M.B.-S.; visualization, O.P.L.; supervision, O.P.L., and R.D.-D.; project administration, O.P.L.; funding acquisition, O.P.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Spanish Ministry of Education, Culture and Sports, grant number FPU16-01888, to A. Acosta-Dacal, and by the University of Las Palmas de Gran Canaria, grant number ULPGC-012-2016, to C. Rial-Berriel.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to some of the data have been obtained from legal proceedings subject to confidentiality.

Acknowledgments

The authors would like to thank Jesús Notario from the Department of Animal Biology, Edaphology and Geology of the University of La Laguna for his help in the identification of the soil samples as clay loam soil, and Ana Macías Montes for her technical assistance in the Toxicology Laboratory of the University of Las Palmas de Gran Canaria.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fernandez-Alvarez, M.; Lamas, J.P.; Garcia-Chao, M.; Garcia-Jares, C.; Llompart, M.; Lores, M.; Dagnac, T. Monitoring of pesticide residues in dairy cattle farms from NW Spain. J. Environ. Monit. 2010, 12, 1864. [Google Scholar] [CrossRef] [PubMed]
  2. Watt, B.E.; Proudfoot, A.T.; Bradberry, S.M.; Vale, J.A. Anticoagulant rodenticides. Toxicol. Rev. 2005, 24, 259–269. [Google Scholar] [CrossRef] [PubMed]
  3. Pelfrène, A.F. Rodenticides. In Hayes’ Handbook of Pesticide Toxicology; Elsevier: Amsterdam, The Netherlands, 2010; pp. 2153–2217. [Google Scholar]
  4. Murphy, M.J. Anticoagulant Rodenticides. In Veterinary Toxicology: Basic and Clinical Principles, 3rd ed.; Elsevier: Amsterdam, The Netherlands, 2018; ISBN 9780128114100. [Google Scholar]
  5. Jackson, W.B.; Brooks, J.E.; Bowerman, A.M.; Kaukeinen, D.E. Anticoagulant resistance in Norway rats as found in US cities. Pest Control 1975, 43, 12–16. [Google Scholar]
  6. Singer, D.E.; Chang, Y.; Fang, M.C.; Borowsky, L.H.; Pomernacki, N.K.; Udaltsova, N.; Go, A.S. The net clinical benefit of warfarin anticoagulation in atrial fibrillation. Ann. Intern. Med. 2009, 151, 297–305. [Google Scholar] [CrossRef] [PubMed]
  7. Gómez-Canela, C.; Lacorte, S. Comprehensive characterization of anticoagulant rodenticides in sludge by liquid chromatography-tandem mass spectrometry. Environ. Sci. Pollut. Res. 2016, 23, 15739–15748. [Google Scholar] [CrossRef] [PubMed]
  8. Gómez-Canela, C.; Barata, C.; Lacorte, S. Occurrence, elimination, and risk of anticoagulant rodenticides and drugs during wastewater treatment. Environ. Sci. Pollut. Res. 2014, 21, 7194–7203. [Google Scholar] [CrossRef]
  9. Jelić, A.; Petrović, M.; Barceló, D. Multi-residue method for trace level determination of pharmaceuticals in solid samples using pressurized liquid extraction followed by liquid chromatography/quadrupole-linear ion trap mass spectrometry. Talanta 2009, 80, 363–371. [Google Scholar] [CrossRef] [PubMed]
  10. Margalida, A.; Bogliani, G.; Bowden, C.G.R.; Donázar, J.A.; Genero, F.; Gilbert, M.; Karesh, W.B.; Kock, R.; Lubroth, J.; Manteca, X.; et al. One Health approach to use of veterinary pharmaceuticals. Science (80-) 2014, 346, 1296–1298. [Google Scholar] [CrossRef] [Green Version]
  11. Díaz-Cruz, M.S.; Barceló, D. LC-MS2 trace analysis of antimicrobials in water, sediment and soil. TrAC Trends Anal. Chem. 2005, 24, 645–657. [Google Scholar] [CrossRef]
  12. Kim, K.R.; Owens, G.; Kwon, S.I.; So, K.H.; Lee, D.B.; Ok, Y.S. Occurrence and environmental fate of veterinary antibiotics in the terrestrial environment. Water. Air. Soil Pollut. 2011, 214, 163–174. [Google Scholar] [CrossRef]
  13. Kaczala, F.; Blum, S.E. The Occurrence of Veterinary Pharmaceuticals in the Environment: A Review. Curr. Anal. Chem. 2015, 12, 169–182. [Google Scholar] [CrossRef] [Green Version]
  14. Gómez-Canela, C.; Vázquez-Chica, A.; Lacorte, S. Comprehensive characterization of rodenticides in wastewater by liquid chromatography-tandem mass spectrometry. Anal. Bioanal. Chem. 2014, 406, 345–358. [Google Scholar] [CrossRef] [PubMed]
  15. Gómez-Canela, C.; Barth, J.A.C.; Lacorte, S. Occurrence and fate of perfluorinated compounds in sewage sludge from Spain and Germany. Environ. Sci. Pollut. Res. 2012, 19, 4109–4119. [Google Scholar] [CrossRef]
  16. Cuthbert, R.J.; Taggart, M.A.; Prakash, V.; Chakraborty, S.S.; Deori, P.; Galligan, T.; Kulkarni, M.; Ranade, S.; Saini, M.; Sharma, A.K.; et al. Avian scavengers and the threat from veterinary pharmaceuticals. Philos. Trans. R. Soc. B Biol. Sci. 2014, 369. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Schaff, J.E.; Montgomery, M.A. An HPLC-HR-MS-MS Method for identification of anticoagulant rodenticides in blood *. J. Anal. Toxicol. 2013, 37, 321–325. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. André, C.; Guyon, C.; Thomassin, M.; Barbier, A.; Richert, L.; Guillaume, Y.C. Association mechanism between a series of rodenticide and humic acid: A frontal analysis to support the biological data. J. Chromatogr. B Anal. Technol. Biomed. Life Sci. 2005, 820, 9–14. [Google Scholar] [CrossRef] [PubMed]
  19. Kemper, N.; Färber, H.; Skutlarek, D.; Krieter, J. Analysis of antibiotic residues in liquid manure and leachate of dairy farms in Northern Germany. Agric. Water Manag. 2008, 95, 1288–1292. [Google Scholar] [CrossRef]
  20. Pitt, W.C.; Berentsen, A.R.; Shiels, A.B.; Volker, S.F.; Eisemann, J.D.; Wegmann, A.S.; Howald, G.R. Non-target species mortality and the measurement of brodifacoum rodenticide residues after a rat (Rattus rattus) eradication on Palmyra Atoll, tropical Pacific. Biol. Conserv. 2015, 185, 36–46. [Google Scholar] [CrossRef] [Green Version]
  21. Obimakinde, S.; Fatoki, O.; Opeolu, B.; Olatunji, O. Veterinary pharmaceuticals in aqueous systems and associated effects: An update. Environ. Sci. Pollut. Res. 2017, 24, 3274–3297. [Google Scholar] [CrossRef]
  22. Arnold, K.E.; Brown, A.R.; Brown, A.R.; Ankley, G.T.; Sumpter, J.P. Medicating the environment: Assessing risks of pharmaceuticals to wildlife and ecosystems. Philos. Trans. R. Soc. B Biol. Sci. 2014, 369. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Bártíková, H.; Skálová, L.; Stuchĺková, L.; Vokřál, I.; Vaněk, T.; Podlipná, R. Xenobiotic-metabolizing enzymes in plants and their role in uptake and biotransformation of veterinary drugs in the environment. Drug Metab. Rev. 2015, 47, 374–387. [Google Scholar] [PubMed]
  24. Blanco, G.; Junza, A.; Barrón, D. Occurrence of veterinary pharmaceuticals in golden eagle nestlings: Unnoticed scavenging on livestock carcasses and other potential exposure routes. Sci. Total Environ. 2017, 586, 355–361. [Google Scholar] [CrossRef]
  25. Ruiz-Suárez, N.; Melero, Y.; Giela, A.; Henríquez-Hernández, L.A.; Sharp, E.; Boada, L.D.; Taylor, M.J.; Camacho, M.; Lambin, X.; Luzardo, O.P.; et al. Rate of exposure of a sentinel species, invasive American mink (Neovison vison) in Scotland, to anticoagulant rodenticides. Sci. Total Environ. 2016, 569–570, 1013–1021. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Gao, X.; Li, H.; Li, H.; Dong, S.; Chu, J.; Guo, H.; Zhao, Q. Sensitive determination of nine anticoagulant rodenticides in blood by high resolution mass spectrometry with supported liquid extraction pretreatment. Forensic Sci. Int. 2018, 292, 39–44. [Google Scholar] [CrossRef] [PubMed]
  27. Barron, L.; Tobin, J.; Paull, B. Multi-residue determination of pharmaceuticals in sludge and sludge enriched soils using pressurized liquid extraction, solid phase extraction and liquid chromatography with tandem mass spectrometry. J. Environ. Monit. 2008, 10, 353–361. [Google Scholar] [CrossRef]
  28. Hernández, A.M.; Bernal, J.; Bernal, J.L.; Martín, M.T.; Caminero, C.; Nozal, M.J. Simultaneous determination of nine anticoagulant rodenticides in soil and water by LC-ESI-MS. J. Sep. Sci. 2013, 36, 2593–2601. [Google Scholar] [CrossRef] [PubMed]
  29. Ogilvie, S.C.; Pierce, R.J.; Wright, G.R.G.; Booth, L.H.; Eason, C.T. Brodifacoum residue analysis in water, soil, invertebrates, and birds after rat eradication on Lady Alice Island. N. Z. J. Ecol. 1997, 21, 195–197. [Google Scholar] [CrossRef]
  30. Medvedovici, A.; David, F.; Sandra, P. Determination of the rodenticides warfarin, diphenadione and chlorophacinone in soil samples by HPLC-DAD. Talanta 1997, 44, 1633–1640. [Google Scholar] [CrossRef]
  31. Chen, W.; Geng, M.; Zhang, L.; Xu, L.; Yuan, H.; Li, C.; Peng, C.; Wang, J.; Zhang, M. Determination of coumarin rodenticides in soils by high performance liquid chromatography. Chinese J. Chromatogr. (Se Pu) 2016, 34, 912–917. [Google Scholar] [CrossRef] [Green Version]
  32. Laktičová, K.; Hromada, R.; Ondrašovič, M.; Legáth, J.; Ďurečko, R.; Ondrašovičová, O.; Nowakovic-Dȩbek, B.; Saba, L. Occurrence of residues of warfarin after its application to cereal crops. Polish J. Environ. Stud. 2009, 18, 405–409. [Google Scholar]
  33. Lee, Y.J.; Choi, J.H.; Abd El-Aty, A.M.; Chung, H.S.; Lee, H.S.; Kim, S.W.; Rahman, M.M.; Park, B.J.; Kim, J.E.; Shin, H.C.; et al. Development of a single-run analytical method for the detection of ten multiclass emerging contaminants in agricultural soil using an acetate-buffered QuEChERS method coupled with LC–MS/MS. J. Sep. Sci. 2017, 40, 415–423. [Google Scholar] [CrossRef] [PubMed]
  34. Ferhi, S.; Bourdat-Deschamps, M.; Daudin, J.J.; Houot, S.; Nélieu, S. Factors influencing the extraction of pharmaceuticals from sewage sludge and soil: An experimental design approach. Anal. Bioanal. Chem. 2016, 408, 6153–6168. [Google Scholar] [CrossRef] [PubMed]
  35. Salvia, M.V.; Vulliet, E.; Wiest, L.; Baudot, R.; Cren-Olivé, C. Development of a multi-residue method using acetonitrile-based extraction followed by liquid chromatography-tandem mass spectrometry for the analysis of steroids and veterinary and human drugs at trace levels in soil. J. Chromatogr. A 2012, 1245, 122–133. [Google Scholar] [CrossRef] [PubMed]
  36. Bragança, I.; Plácido, A.; Paíga, P.; Domingues, V.F.; Delerue-Matos, C. QuEChERS: A new sample preparation approach for the determination of ibuprofen and its metabolites in soils. Sci. Total Environ. 2012, 433, 281–289. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  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] [PubMed] [Green Version]
  38. Perestrelo, R.; Silva, P.; Porto-Figueira, P.; Pereira, J.A.M.; Silva, C.; Medina, S.; Câmara, J.S. QuEChERS—Fundamentals, relevant improvements, applications and future trends. Anal. Chim. Acta 2019, 1070, 1–28. [Google Scholar] [CrossRef] [PubMed]
  39. Asensio-Ramos, M.; Hernández-Borges, J.; Ravelo-Pérez, L.M.; Rodríguez-Delgado, M.A. Evaluation of a modified QuEChERS method for the extraction of pesticides from agricultural, ornamental and forestal soils. Anal. Bioanal. Chem. 2010, 396, 2307–2319. [Google Scholar] [CrossRef] [PubMed]
  40. Vieira, N.M.; Paiva, M.C.G.; da Silva, A.A.; Costa, A.I.G.; de Freitas, F.C.L.; da Costa Marinho, M.I. Validation and application of a QuEChERS based method for estimation of half-life of imidazolinone herbicides in soils by LC-ESI-MS/MS. Ecotoxicol. Environ. Saf. 2018, 167, 212–217. [Google Scholar] [CrossRef]
  41. Acosta-Dacal, A.; Rial-Berriel, C.; Díaz-Díaz, R.; del Mar Bernal-Suárez, M.; Luzardo, O.P. Optimization and validation of a QuEChERS-based method for the simultaneous environmental monitoring of 218 pesticide residues in clay loam soil. Sci. Total Environ. 2021, 753, 142015. [Google Scholar] [CrossRef]
  42. He, Z.; Wang, L.; Peng, Y.; Luo, M.; Wang, W.; Liu, X. Determination of selected polychlorinated biphenyls in soil and earthworm (Eisenia fetida) using a QuEChERS-based method and gas chromatography with tandem MS. J. Sep. Sci. 2015, 38, 3766–3773. [Google Scholar] [CrossRef]
  43. Rouvire, F.; Buleté, A.; Cren-Olivé, C.; Arnaudguilhem, C. Multiresidue analysis of aromatic organochlorines in soil by gas chromatography-mass spectrometry and QuEChERS extraction based on water/dichloromethane partitioning. Comparison with accelerated solvent extraction. Talanta 2012, 93, 336–344. [Google Scholar] [CrossRef]
  44. Cao, X.; Yang, X.; Liu, Z.; Jiao, H.; Liu, S.; Liu, L.; Meng, Q. Rapid Simultaneous Screening and Detection of 12 Anticoagulant Rodenticides in Food by Ultra-performance Liquid Chromatography-Triple Quadrupole/Linear Ion Trap Tandem Mass Spectrometry. Food Anal. Methods 2017, 10, 3538–3547. [Google Scholar] [CrossRef]
  45. Zhu, F.; Liu, H.; Chen, B.; Rong, W.; Ma, Y. Rapid simultaneous screening and detection of ten anticoagulant rodenticides in foods by ultra-high performance liquid chromatography-tandem mass spectrometry. Chin. J. Chromatogr. 2013, 31, 473. [Google Scholar] [CrossRef] [PubMed]
  46. Gómez-Ramírez, P.; Martínez-López, E.; Navas, I.; María-Mojica, P.; García-Fernández, A.J. A modification of QuEChERS method to analyse anticoagulant rodenticides using small blood samples. Rev. Toxicol. 2012, 29, 10–14. [Google Scholar]
  47. Rial-Berriel, C.; Acosta-Dacal, A.; Zumbado, M.; Luzardo, O.P. Micro QuEChERS-based method for the simultaneous biomonitoring in whole blood of 360 toxicologically relevant pollutants for wildlife. Sci. Total Environ. 2020, 736. [Google Scholar] [CrossRef] [PubMed]
  48. Almeida, C.M.M. Overview of Sample Preparation and Chromatographic Methods to Analysis Pharmaceutical Active Compounds in Waters Matrices. Separations 2021, 8, 16. [Google Scholar] [CrossRef]
  49. Vudathala, D.; Cummings, M.; Murphy, L. Analysis of multiple anticoagulant rodenticides in animal blood and liver tissue using principles of QuEChERS method. J. Anal. Toxicol. 2010, 34, 273–279. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Orfanidis, A.; Gika, H.; Theodoridis, G.; Mastrogianni, O.; Raikos, N. Development of a UHPLC-MS/MS method for the determination of 84 pharmaceuticals and drugs of abuse in human liver. J. Chromatogr. B Anal. Technol. Biomed. Life Sci. 2020, 1151, 122192. [Google Scholar] [CrossRef] [PubMed]
  51. Martínez-Huelamo, M.; Jiménez-Gámez, E.; Hermo, M.P.; Barrón, D.; Barbosa, J. Determination of penicillins in milk using LC-UV, LC-MS and LC-MS/MS. J. Sep. Sci. 2009, 32, 2385–2393. [Google Scholar] [CrossRef]
  52. Hroboňová, K.; Sádecká, J.; Čižmárik, J. HPLC separation and determination of dicoumarol and other simple coumarins in sweet clover. Nov. Biotechnol. Chim. 2018, 17, 95–102. [Google Scholar] [CrossRef]
  53. Hernández, A.M.; Bernal, J.; Bernal, J.L.; Martín, M.T.; Caminero, C.; Nozal, M.J. Analysis of anticoagulant rodenticide residues in Microtus arvalis tissues by liquid chromatography with diode array, fluorescence and mass spectrometry detection. J. Chromatogr. B Anal. Technol. Biomed. Life Sci. 2013, 925, 76–85. [Google Scholar] [CrossRef]
  54. Gallocchio, F.; Basilicata, L.; Benetti, C.; Angeletti, R.; Binato, G. Multi-residue determination of eleven anticoagulant rodenticides by high-performance liquid chromatography with diode array/fluorimetric detection: Investigation of suspected animal poisoning in the period 2012–2013 in north-eastern Italy. Forensic Sci. Int. 2014, 244, 63–69. [Google Scholar] [CrossRef] [PubMed]
  55. Tuerk, J.; Reinders, M.; Dreyer, D.; Kiffmeyer, T.K.; Schmidt, K.G.; Kuss, H.M. Analysis of antibiotics in urine and wipe samples from environmental and biological monitoring—Comparison of HPLC with UV-, single MS- and tandem MS-detection. J. Chromatogr. B Anal. Technol. Biomed. Life Sci. 2006, 831, 72–80. [Google Scholar] [CrossRef] [PubMed]
  56. Taylor, M.J.; Giela, A.; Sharp, E.A.; Senior, C.C.; Vyas, D.S. A rapid multi-class, multi-residue UHPLC-MS/MS method for the simultaneous determination of anticoagulant rodenticides, pesticides and veterinary medicines in wild animals, pets and livestock. Anal. Methods 2019, 11, 1087–1101. [Google Scholar] [CrossRef]
  57. Bourdat-Deschamps, M.; Leang, S.; Bernet, N.; Daudin, J.J.; Nélieu, S. Multi-residue analysis of pharmaceuticals in aqueous environmental samples by online solid-phase extraction-ultra-high-performance liquid chromatography-tandem mass spectrometry: Optimisation and matrix effects reduction by quick, easy, cheap, effective, rugged and safe extraction. J. Chromatogr. A 2014, 1349, 11–23. [Google Scholar] [CrossRef] [PubMed]
  58. SANTE/12682 Guidance document on analytical quality control and validation procedures for pesticide residues analysis in food and feed. Eur. Comm. Health Consum. Prot. Dir. 2019, 2–44. Available online: https://www.eurl-pesticides.eu/docs/public/tmplt_article.asp?CntID=727 (accessed on 8 March 2021).
  59. Lehotay, S.J.; O’Neil, M.; Tully, J.; García, A.V.; Contreras, M.; Mol, H.; Heinke, V.; Anspach, T.; Lach, G.; Fussell, R.; et al. 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] [Green Version]
  60. Acosta-Dacal, A.; Rial-Berriel, C.; Díaz-Díaz, R.; del Mar Bernal-Suárez, M.; Zumbado, M.; Henríquez-Hernández, L.A.; Luzardo, O.P. Supporting dataset on the optimization and validation of a QuEChERS-based method for the determination of 218 pesticide residues in clay loam soil. Data Br. 2020, 33, 106393. [Google Scholar] [CrossRef] [PubMed]
  61. SANCO/825/00 Guidance Document on Pesticide Residue Analytical Methods. 2010. Available online: https://ec.europa.eu/food/sites/food/files/plant/docs/pesticides_ppp_app-proc_guide_res_post-reg-cont-monitor.pdf (accessed on 8 March 2021).
  62. Raposo, F. Evaluation of analytical calibration based on least-squares linear regression for instrumental techniques: A tutorial review. TrAC Trends Anal. Chem. 2016, 77, 167–185. [Google Scholar] [CrossRef]
  63. Kamel, A.M.; Brown, P.R.; Munson, B. Effects of mobile-phase additives, solution pH, ionization constant, and analyte concentration on the sensitivities and electrospray ionization mass spectra of nucleoside antiviral agents. Anal. Chem. 1999, 71, 5481–5492. [Google Scholar] [CrossRef]
  64. King, R.; Bonfiglio, R.; Fernandez-Metzler, C.; Miller-Stein, C.; Olah, T. Mechanistic investigation of ionization suppression in electrospray ionization. J. Am. Soc. Mass Spectrom. 2000, 11, 942–950. [Google Scholar] [CrossRef] [Green Version]
  65. Kruve, A.; Kaupmees, K. Adduct Formation in ESI/MS by Mobile Phase Additives. J. Am. Soc. Mass Spectrom. 2017, 28, 887–894. [Google Scholar] [CrossRef] [PubMed]
  66. Danaher, M.; O’Keeffe, M.; Glennon, J.D.; Howells, L. Development and optimisation of an improved derivatisation procedure for the determination of avermectins and milbemycins in bovine liver. Analyst 2001, 126, 576–580. [Google Scholar] [CrossRef] [PubMed]
Toxics 09 00083 g001 550
Figure 1. LC-MS/MS chromatograms of a blank soil sample spiked with the anticoagulant rodenticides (ARs) and the pharmaceutical active compounds (PhACs) at 50 ng g−1.
Figure 1. LC-MS/MS chromatograms of a blank soil sample spiked with the anticoagulant rodenticides (ARs) and the pharmaceutical active compounds (PhACs) at 50 ng g−1.
Toxics 09 00083 g001
Toxics 09 00083 g002 550
Figure 2. Matrix effect of the ARs and PhACs in the soil used for the validation process. Bars represent the mean recoveries and the SD of the 43 analytes. PhACs (orange) and ARs (green) are presented according to their number in Table 1. The suppression has been indicated with an oblique dashed pattern.
Figure 2. Matrix effect of the ARs and PhACs in the soil used for the validation process. Bars represent the mean recoveries and the SD of the 43 analytes. PhACs (orange) and ARs (green) are presented according to their number in Table 1. The suppression has been indicated with an oblique dashed pattern.
Toxics 09 00083 g002
Toxics 09 00083 g003 550
Figure 3. Chromatogram of the regenerated water (A) and the two agricultural soils with the identification of (B) brodifacoum and (C) flocoumafen.
Figure 3. Chromatogram of the regenerated water (A) and the two agricultural soils with the identification of (B) brodifacoum and (C) flocoumafen.
Toxics 09 00083 g003
Toxics 09 00083 g004 550
Figure 4. Chromatogram of the soil after the irrigation simulation experiment.
Figure 4. Chromatogram of the soil after the irrigation simulation experiment.
Toxics 09 00083 g004
Toxics 09 00083 g005 550
Figure 5. Chromatogram with the identification of brodifacoum, bromadiolone and difenacoum in soil underneath the remains of dead animal in the context of a forensic investigation.
Figure 5. Chromatogram with the identification of brodifacoum, bromadiolone and difenacoum in soil underneath the remains of dead animal in the context of a forensic investigation.
Toxics 09 00083 g005
Table
Table 1. Characterization of target compounds by therapeutical classes or biocides, legal status, and their chromatographic and mass spectrometric conditions.
Table 1. Characterization of target compounds by therapeutical classes or biocides, legal status, and their chromatographic and mass spectrometric conditions.
N°CompoundCategory aLegal Status in the EU btR (min)PolarityQuantificationConfirmationFragmentor Voltage (V)
MRM Transition (m/z)CE (eV)MRM transition (m/z)CE (eV)
1SulfanilamidePhACs, ABApproved1.06Positive173.0 → 93.024173.0 → 76.150126
2SulfacetamidePhACs, ABApproved2.28Positive215.3 → 92.020215.3 → 65.34590
3MetronidazolePhACs, ABApproved2.73Positive172.1 → 128.012172.1 → 82.12498
4SulfadiacinePhACs, ABApproved2.82Positive251.0 → 92.028251.0 → 156.012111
5SulfapyridinePhACs, ABApproved3.16Positive250.0 → 92.028250.0 → 156.012126
6SulfameracinePhACs, ABNot approved3.29Positive265.0 → 92.028265.0 → 156.012126
7SulfametizolePhACs, ABNot approved3.59Positive271.0 → 155.98271.0 → 92.028103
8SulfametacinePhACs, ABNot approved3.68Positive279.1 → 186.012279.1 → 92.032134
9SulfametoxipiridacinePhACs, ABNot approved3.76Positive281.0 → 155.912281.0 → 92.128121
10SulfachloropiridacinePhACs, ABNot approved3.89Positive285.0 → 156.012285.0 → 92.128101
11SulfametoxazolePhACs, ABApproved3.95Positive254.0 → 92.028254.0 → 156.012111
12SulfamonomethoxinePhACs, ABNot approved4.02Positive281.1 → 92.114281.1 → 156.032120
13SulfadoxinePhACs, ABApproved4.13Positive311.1 → 156.016311.1 → 92.032126
14SulfisoxazolePhACs, ABNot approved4.14Positive268.0 → 156.08268.0 → 92.124106
15ChloramphenicolPhACs, ABApproved4.76Negative321.0 → 152.14323.0 → 152.14113
16SulfadimetoxinePhACs, ABApproved4.83Positive311.0 → 156.016311.0 → 92.032139
17SulfaquinoxalinePhACs, ABApproved4.99Positive301.0 → 156.012301.0 → 92.132159
18Cefuroxime axetil (two isomers)PhACs, ABApproved5.40Positive533.0 → 447.015533.0 → 386.020160
19OxfendazolePhACs, AHApproved5.64Positive316.1 → 159.032316.1 → 191.116166
20Penicillin VPhACs, ABApproved6.48Positive383.2 → 159.910383.2 → 113.940130
21MebendazolePhACs, AHApproved6.73Positive296.1 → 264.120296.1 → 77.048151
22CloxacillinPhACs, ABApproved7.09Positive468.1 → 159.98468.1 → 177.820126
23DexamethasonePhACs, GCApproved7.18Positive393.2 → 373.22393.2 → 355.26103
24AlbendazolePhACs, AHApproved7.27Positive266.1 → 234.116266.1 → 191.032155
25KetoprofenPhACs, NSAIDApproved7.36Positive255.1 → 209.18255.1 → 77.148123
26JosamycinPhACs, ABApproved7.52Positive860.5 → 173.940860.5 → 108.940200
27NaproxenPhACs, NSAIDApproved7.60Positive231.1 → 185.010231.1 → 169.913120
28CortiscosteronePhACs, GCApproved7.91Positive389.1 → 329.013389.1 → 371.01380
29FenbendazolePhACs, AHApproved8.11Positive300.1 → 268.120300.1 → 159.036156
30FlunixinPhACs, NSAIDApproved8.11Positive297.1 → 279.124297.1 → 264.132141
31DiclofenacPhACs, NSAIDApproved8.75Positive296.0 → 215.116296.0 → 214.148103
32ImipenemPhACs, ABApproved8.75Positive300.0 → 125.915300.0 → 98.01550
33Mefenamic acidPhACs, NSAIDApproved9.52Positive242.1 → 209.128242.1 → 180.144108
34Tolfenamic acidPhACs, NSAIDApproved9.78Negative260.0 → 216.18260.0 → 35.120108
35EprinomectinPhACs, ABApproved10.83Positive878.5 → 186.015878.5 → 154.045160
36MoxidectinPhACs, AHApproved11.25Positive641.4 → 529.25641.4 → 499.25100
37WarfarinARsNot approved7.87Negative307.1 → 161.120307.1 → 250.120140
38CoumatetralylARsNot approved8.28Negative291.1 → 141.030291.1 → 247.020140
39BromadioloneARsApproved9.74Negative525.3 → 250.040527.3 → 250.040200
40DifenacoumARsNot approved10.27Negative443.2 → 135.040443.2 → 293.035200
41FlocoumafenARsNot approved10.36Negative541.3 → 382.025541.3 → 161.040230
42BrodifacoumARsNot approved10.67Negative521.3 → 79.050523.3 → 135.045220
43DifethialoneARsNot approved10.80Negative537.3 → 79.050537.3 → 151.045220
Cyromazine-d4P-IS–1.58Positive171.0 → 86.015171.0 → 129.015100
Carbendazim-d3P-IS–3.45Positive195.1 → 160.115195.1 → 131.930100
Pirimicarb-d6P-IS–5.12Positive245.2 → 78.25245.2 → 185.11570
Atrazine-d5P-IS–6.66Positive221.2 → 179.015221.2 → 101.03090
Linuron-d3P-IS–7.45Positive255.1 → 159.815255.1 → 185.015100
CoumachlorP-IS–8.52Negative341.1 → 161.015341.1 → 284.015120
CE: Collision Energy; tR: Retention time. a PhACs—pharmaceuticals active compound, ARs—anticoagulant rodenticides, AH—anthelminthic, AB—antibiotic, NSAID—nonsteroidal anti-inflamatory drug, GC—glucocorticoid, P-IS—Procedural Internal Standard b For human and veterinary drugs, the marketing status in Spain is specified, as shown in the Cima and Cimavet search engine of the Spanish agency for drugs and health products (https://cima.aemps.es/cima/publico/home.html accessed on 8 March 2021; https://cimavet.aemps.es/cimavet/publico/home.html accessed on 8 March 2021). For rodenticides, the legal status reflecting the EU Pesticide Database was considered (https://ec.europa.eu/food/plant/pesticides/eu-pesticidesdatabase/public/?event=activesubstance.selection&language=EN accessed on 8 March 2021), which is valid for the entire EU.—Not applicable.
Table
Table 2. Linear studies and threshold limits of the ARs and PhACs.
Table 2. Linear studies and threshold limits of the ARs and PhACs.
N°CompoundGroupLinear Range (ng g−1)R2p-Value (Mandel Test) aLOD (ng g−1)LOQ (ng g−1)
1SulfanilamidePhACs1.56–250.99890.16111.5602.5
2SulfacetamidePhACs0.195–500.99950.55110.1950.5
3MetronidazolePhACs0.195–250.99770.44880.1950.5
4SulfadiacinePhACs0.195–500.99960.41220.1950.5
5SulfapyridinePhACs0.39–500.9990.20550.3900.5
6SulfameracinePhACs0.39–250.99860.98100.3901.0
7SulfametizolePhACs1.56–250.99450.23241.56010.0
8SulfametacinePhACs0.39–250.99720.18260.3900.5
9SulfametoxipiridacinePhACs0.39–250.99740.92640.3900.5
10SulfachloropiridacinePhACs0.39–250.99310.41490.3902.5
11SulfametoxazolePhACs0.195–250.99890.17820.1952.5
12SulfamonomethoxinePhACs0.78–500.99560.06960.7801.0
13SulfadoxinePhACs0.195–250.99960.33570.1950.5
14SulfisoxazolePhACs0.78–250.99750.82310.7801.0
15ChloramphenicolPhACs6.25–500.99000.88666.25020.0
16SulfadimetoxinePhACs0.195–500.99980.19560.1950.5
17SulfaquinoxalinePhACs0.78–500.99870.05820.7801.0
18Cefuroxime axetil (two isomers)PhACs0.39–500.99290.37220.3901.0
19OxfendazolePhACs0.39–500.99980.39310.3900.5
20Penicillin VPhACs1.56–250.99870.85501.5602.5
21MebendazolePhACs0.195–500.99970.13910.1950.5
22CloxacillinPhACs1.56–500.99040.10941.5602.5
23DexamethasonePhACs1.56–500.99740.08261.5605.0
24AlbendazolePhACs0.097–250.99690.63870.0970.5
25KetoprofenPhACs0.39–500.99850.70220.3900.5
26JosamycinPhACs0.39–250.99940.46660.3901.0
27NaproxenPhACs3.125–500.99180.30953.12520.0
28CortiscosteronePhACs1.56–250.99740.67111.5605.0
29FenbendazolePhACs0.048–250.99790.55510.0480.5
30FlunixinPhACs0.097–250.99970.46750.0970.5
31ImipenemPhACs6.25–500.99960.24996.25050.0
32DiclofenacPhACs0.097–500.99380.73540.0972.5
33Mefenamic acidPhACs0.78–250.99670.14250.78010.0
34Tolfenamic acidPhACs6.25–500.99710.49076.25010.0
35EprinomectinPhACs0.78–250.99740.25000.7802.5
36MoxidectinPhACs1.56–250.99720.12691.56020.0
37WarfarinARs0.39–250.99410.14000.3901.0
38CoumatetralylARs0.195–250.99770.10550.1952.5
39BromadioloneARs0.195–250.99910.16910.1951.0
40DifenacoumARs0.097–500.99430.22210.0971.0
41FlocoumafenARs0.024–500.99830.21670.0241.0
42BrodifacoumARs0.097–500.99970.52000.1951.0
43DifethialoneARs3.125–500.99460.07973.12520.0
a According to Mandel’s test, the null hypothesis (H0) would indicate that the linear model is adequate to describe the calibration data. Since p-value > 0.05 in the range indicated in the table, we accept the H0 and would conclude that the linear model is adequate or reasonable to model the calibration data within this range.
Table
Table 3. Recoveries (Rec) and relative standard deviation (RSD).
Table 3. Recoveries (Rec) and relative standard deviation (RSD).
N°.CompoundGroup0.5 (ng g−1)1.0 (ng g−1)2.5 (ng g−1)5 (ng g−1)10.0 (ng g−1)20.0 (ng g−1)50.0 (ng g−1)
Rec (%)RSD (%)Rec (%)RSD (%)Rec (%)RSD (%)Rec (%)RSD (%)Rec (%)RSD (%)Rec (%)RSD (%)Rec (%)RSD (%)
1SulfanilamidePhACsN/AN/AN/AN/A112.17.990.724.291.317.979.99.063.312.8
2SulfacetamidePhACs83.315.163.210.359.88.178.016.558.65.462.85.975.53.7
3MetronidazolePhACs75.316.886.526.987.57.3102.513.194.48.097.75.592.71.6
4SulfadiacinePhACs88.113.876.08.274.45.083.813.775.15.178.16.977.62.1
5SulfapyridinePhACs111.018.978.216.172.25.279.221.875.29.768.77.769.73.2
6SulfameracinePhACsN/AN/A91.619.481.32.785.714.778.25.180.59.278.93.7
7SulfametizolePhACsN/AN/AN/AN/AN/AN/AN/AN/AN/A9.871.213.468.63.8
8SulfametacinePhACs106.89.984.017.878.511.483.918.878.93.283.97.785.41.5
9SulfametoxipiridacinePhACs84.110.572.021.662.57.882.717.975.56.174.25.875.02.6
10SulfachloropiridacinePhACsN/AN/AN/AN/A81.616.285.722.784.36.087.76.789.62.8
11SulfametoxazolePhACsN/AN/AN/AN/A89.313.196.416.793.37.2101.22.4103.23.0
12SulfamonomethoxinePhACsN/AN/A75.210.766.67.680.230.775.84.480.39.180.66.1
13SulfadoxinePhACs96.59.786.16.785.24.796.816.090.06.794.05.989.91.4
14SulfisoxazolePhACsN/AN/A89.813.977.415.782.818.971.310.680.17.080.83.5
15ChloramphenicolPhACsN/AN/AN/AN/AN/AN/AN/AN/AN/AN/A117.011.0117.89.0
16SulfadimetoxinePhACs85.710.386.77.484.84.091.915.686.26.090.95.788.21.3
17SulfaquinoxalinePhACsN/AN/A111.18.882.97.584.716.674.15.777.07.478.92.1
18Cefuroxime axetil (two isomers)PhACsN/AN/A83.418.296.16.8117.211.5109.56.7115.55.3104.30.7
19OxfendazolePhACs115.113.183.49.865.22.982.313.869.95.366.24.763.40.4
20Penicillin VPhACsN/AN/AN/AN/A66.517.669.726.763.018.265.912.873.29.9
21MebendazolePhACs92.415.479.813.475.14.696.018.080.75.680.76.479.31.5
22CloxacillinPhACsN/AN/AN/AN/A76.115.4122.318.698.79.2101.89.7126.011.7
23DexamethasonePhACsN/AN/AN/AN/AN/AN/A71.117.668.09.175.46.269.93.6
24AlbendazolePhACs76.39.076.38.269.53.899.516.582.96.376.25.975.81.1
25KetoprofenPhACs103.518.399.813.876.66.793.420.792.16.889.64.987.23.2
26JosamycinPhACsN/AN/A106.514.574.510.080.820.270.37.671.810.474.15.6
27NaproxenPhACsN/AN/AN/AN/AN/AN/AN/AN/AN/AN/A86.918.194.25.2
28CortiscosteronePhACsN/AN/AN/AN/AN/AN/A112.915.493.96.3101.310.899.75.4
29FenbendazolePhACs93.912.287.18.983.06.2107.014.692.25.390.66.285.31.4
30FlunixinPhACs81.011.370.217.368.64.281.917.071.84.269.84.265.02.7
31ImipenemPhACsN/AN/AN/AN/AN/AN/AN/AN/AN/AN/AN/AN/A91.717.2
32DiclofenacPhACsN/AN/AN/AN/A110.317.9107.128.884.417.793.98.185.12.1
33Mefenamic acidPhACsN/AN/AN/AN/AN/AN/AN/AN/A82.714.383.117.282.55.9
34Tolfenamic acidPhACsN/AN/AN/AN/AN/AN/AN/AN/A121.515.4121.27.8111.510.2
35EprinomectinPhACsN/AN/AN/AN/A102.112.7104.015.498.311.892.19.584.82.6
36MoxidectinPhACsN/AN/AN/AN/AN/AN/AN/AN/AN/AN/A95.014.996.39.6
37WarfarinARsN/AN/A81.29.2103.84.6106.014.7113.56.5105.54.985.85.0
38CoumatetralylARsN/AN/AN/AN/A104.613.5106.512.8107.313.099.33.387.65.9
39BromadioloneARsN/AN/A92.814.487.515.490.913.893.67.791.47.482.42.6
40DifenacoumARsN/AN/A110.212.990.66.5110.214.3107.85.4102.67.486.31.4
41FlocoumafenARsN/AN/A105.013.094.116.8108.111.6108.98.0100.85.085.11.9
42BrodifacoumARsN/AN/A106.511.587.66.598.414.898.08.693.44.284.34.0
43DifethialoneARsN/AN/AN/AN/AN/AN/AN/AN/AN/AN/A113.74.995.76.2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Acosta-Dacal, A.; Rial-Berriel, C.; Díaz-Díaz, R.; Bernal-Suárez, M.d.M.; Zumbado, M.; Henríquez-Hernández, L.A.; Luzardo, O.P. An Easy Procedure to Quantify Anticoagulant Rodenticides and Pharmaceutical Active Compounds in Soils. Toxics 2021, 9, 83. https://0-doi-org.brum.beds.ac.uk/10.3390/toxics9040083

AMA Style

Acosta-Dacal A, Rial-Berriel C, Díaz-Díaz R, Bernal-Suárez MdM, Zumbado M, Henríquez-Hernández LA, Luzardo OP. An Easy Procedure to Quantify Anticoagulant Rodenticides and Pharmaceutical Active Compounds in Soils. Toxics. 2021; 9(4):83. https://0-doi-org.brum.beds.ac.uk/10.3390/toxics9040083

Chicago/Turabian Style

Acosta-Dacal, Andrea, Cristian Rial-Berriel, Ricardo Díaz-Díaz, María del Mar Bernal-Suárez, Manuel Zumbado, Luis Alberto Henríquez-Hernández, and Octavio P. Luzardo. 2021. "An Easy Procedure to Quantify Anticoagulant Rodenticides and Pharmaceutical Active Compounds in Soils" Toxics 9, no. 4: 83. https://0-doi-org.brum.beds.ac.uk/10.3390/toxics9040083

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