Abstract

Antibiotics and pesticides are widespread in most rivers and lakes due to the overuse of antibiotics and pesticides, but there are few methods for simultaneous analysis of antibiotics and pesticides in aquatic environments. To address this knowledge gap, a concise and sensitive analytical method is proposed in which three classes of human and veterinary drugs (sulfonamides, macrolides, and hormones) and two classes of pesticides (organophosphorus and neonicotinoids) are simultaneously extracted and determined in surface water. The solid-phase extraction column with Cleanert PEP-2 was preconditioned sequentially with 6 mL of methanol, ultrapure water, and citric acid buffer (pH 3.0) each for simultaneous extraction and further purification. The forty-seven target analytes were analysed by LC-MS/MS in positive and negative ion modes. The LC separation was performed using a Sigma-Aldrich C18 column with 0.1% formic acid in water and acetonitrile as a gradient eluting mobile phase in positive ion mode. The internal standard method was used to overcome the inevitable matrix effects in LC-MS/MS analysis. The matrix effects of most target analytes were in the range of 27–151%. The recoveries of forty analytes in the three concentrations (10, 50, and 100 ng L−1) of surface water spiked samples ranged from 41 to 127%. The method quantitative limits of the analytes were in the range of 0.40–5.49 ng L−1. Application of the method to analyze samples in the eight runoff outlets of the Pearl River Delta showed that some antibiotics and pesticides were detected, and the concentration of parathion was as high as 154 ng L−1. A powerful tool for quickly and efficiently screening for contaminants in surface water has been presented.

1. Introduction

Human and veterinary drugs (HVDs), considered emerging pollutants, are attracting increased attention worldwide [1, 2]. They have been increasingly administered to prevent and treat microbial infections in humans and animals [35] and are commonly added to animal feeds as supplements to promote the growth of animals in livestock and aquaculture [68]. Upon administration, fractions of pharmaceuticals can enter into different environmental compartments in their parent forms or as metabolites through excretion of urine and feces [9]. These chemical residues can be further discharged into aquatic environments through the direct disposal of human and animal wastes [10, 11] or effluents from municipal sewage treatment plants [12, 13]. HVDs are not easy to combine with solid materials. Parent drugs and their metabolites are often hydrophilic and easily reach the aquatic environment, posing a threat to surface and ground water. Consequently, HVDs have been detected in surface water [14] and groundwater [15] in bioactive and bioavailable forms, leading to undesirable effects on nontarget organisms and drinking water supplies, as well as contamination of food [16]. For example, HVDs have been detected in a number of river systems in China [1726], especially sulfonamides (SAs) and macrolides (MLs) which have been detected frequently, indicating that bacterial resistance and potential ecological risks caused by HVDs in water should not be overlooked.

Organophosphates (OPs) and neonicotinoids (NNIs) have been used to protect crops against pests as they exhibit a broad spectrum of insecticidal activities and low toxicity to mammals [27]. The use of organophosphorus pesticides and neonicotinoid pesticides accounted for one-third of the global pesticide market share. Although some OPs have been banned and replaced by NNIs, pollution of the aquatic environment has become a serious problem due to the toxicity and abuse of OPs, causing a risk to human health [28]. NNIs may have potential ecotoxic effects on pollinators such as bees and invertebrates [29]. The European Union banned the use of NNIs on flowering crops that attract bees in 2013 [30]. However, NNIs remain widely used around the world due to their low price and significant effects. These pesticides may remain in agricultural products consumed by humans and enter into the digestive tract. They can also be washed off plant surfaces by rain and discharged into river systems. Numerous methods have been developed to detect OPs and NNIs in fruits [31, 32], vegetables [33, 34], and honey [35, 36].

Despite the availability of numerous methods for separate detection of HVDs and pesticide residues in a variety of sample matrices [31, 32, 3639], relatively fewer methods are available for simultaneously purifying and measuring SAs, MLs, hormones (HMs), OPs, and NNIs in environmental samples. Previous studies on pretreatment and instrumental analysis of surface water samples have targeted one group or a narrow range of compounds only [39, 40]. For combined extraction, clean-up, and analytical separation, there seems to be a considerable challenge due to the amphoteric nature, the wide range of chemical families, and the wide polarity range of pollutants present.

To address this issue, the present study is set out to develop a rapid, sensitive, and validated method for simultaneously extracting, purifying, and detecting 47 analytes extracted from surface water. Target analytes include 47 human and veterinary drugs and pesticide residues (Supplemental Data (SD); Table S1), which are commonly used chemicals in aquaculture and agriculture. As a worldwide population and agricultural country, the dramatic increase in the use of HVDs and pesticides has led to the frequent detection of these chemicals in most rivers in China [1727]. Therefore, it is very urgent and necessary to provide an analytical method for rapid screening of multiple pollutants. Considering that simultaneous detection of multiple analytes may cause mutual interference, the instrument methods were first optimized, and the analytes were effectively separated via liquid chromatography. For sample pretreatment, a solid-phase extraction method was selected and optimized. It was verified by a matrix-spiked experiment and successfully applied to environmental water samples. The results showed that this pretreatment method has acceptable and stable recoveries.

2. Materials and Methods

2.1. Materials and Chemicals

Standards of SAs and OPs were obtained from O2si Smart Solutions (Charleston, SC, USA). Other MLs and HMs were purchased from ChemTek (Worcester, MA, USA). The ten NNIs compounds were obtained from Manhage Bio-Tech (Beijing, China) (Supplemental Data; Table S1). Dinotefuran, imidaclothiz, and isotope-labeled standards (sulfachloropyridazine-d3 and clothianidin-d3) were purchased from Dr. Ehrenstorfer (Augsburg, Germany). Isotope-labeled standards imidaclothiz-d4 and acetamiprid-d3 were obtained from C/D/N isotopes (Quebec, Canada), while estradiol-d2 and diethylstilbrstrol-d8 were obtained from Manhage Bio-Tech (Beijing, China). Another three isotope-labeled standards tilmicosin-d3, sulfamethoxzole-d4, and sulfamethazine-d4 were obtained from Toronto Research Chemicals (North York, ON, Canada). All standards of the target analytes are solutions (concentration of 100 mg/L), and the isotope-labeled standards are solid. The solid standards were weighed with an electronic balance (Ohaus Instruments, Changzhou, China) with an accuracy of 0.00001 g and dissolved in methanol (concentration of 100 mg/L) prior to use.

Formic acid (LCMS grade) and disodium ethylenediamine tetraacetate (Na2EDTA, LC grade) were purchased from Sigma-Aldrich (St. Louis, MO, USA). Citric acid monohydrate and trisodium citrate dihydrate were purchased from Chemical Reagent Factory (Guangzhou, China). Ultrapure water was provided by a MilliQ water purification system (Millipore, Billerica, Germany). Methanol and acetonitrile of HPLC and LCMS grades were purchased from Oceanpak (Gothenburg, Sweden) and Thermo Fisher Scientific (Waltham, MA, USA), respectively. The Cleanert PEP-2 cartridges (6 mL, 200 mg) were supplied by Agela Technologies (Tianjin, China). Glass fiber filters (pore size 0.7 μm) were purchased from Whatman International (Maidstone, England) and heated at 450°C for 4 h in order to pyrolyze any organic material present.

2.2. Optimization of Sample Pretreatment

Solid-phase extraction cartridges were selected for extraction and purification. The pH value of water samples and the amount of Na2EDTA used were significant to the recovery of antibiotics [39]. These two parameters were optimized as part of method development efforts. The separation efficacy of a liquid chromatography column is directly related to its diameter and length, the particle size of the packing materials, additives, mobile phase gradients, and flow rates. These parameters were thus rigorously examined with the analytes selected.

Water samples were filtered through 0.7 μm glass fiber filters. One liter of filtered water was acidified to pH 3.0 with 1.0 mol L−1 H2SO4, followed by the addition of 0.2 g Na2EDTA. Surrogate standards (sulfamethazine-d4, sulfamethoxypyridazine-d3, acetamiprid-d3, clothianidin-d3, and diethylstilbrstrol-d8 at 1 μg mL-1; 20 μL) were added to each sample prior to extraction. A Cleanert PEP-2 cartridge (6 mL, 200 mg) was preconditioned sequentially with 6 mL of methanol, ultrapure water, and citric acid buffer (pH 3.0) each. Thereafter, the samples were passed through the cartridge at a flow rate of 3.0 mL min−1. The cartridge was then rinsed with 6 mL of ultrapure water and dried under N2 gas for 1 h. Each dried cartridge was eluted with 6 mL of methanol. Target analytes were collected in a 10 mL glass vial, volume-reduced under N2 purge to approximately 10 μL, and dissolved in a mixture of methanol : water (1 : 9 in volume) containing 0.1% formic acid to a final volume of 1 mL. Suspended particles were removed by 0.22 μm membrane filtration. Internal standards (sulfamethoxazole-d4, tilmicosin-d3, imidacloprid-d4, and estradiol-d2 at 1 μg mL−1; 20 μL) were added to each extract which was stored at −18°C until LC-MS/MS analysis. The retention time of the internal standard is close to that of the analyte and evenly distributed.

Matrix-spiked samples were prepared from water (collected from the park of Qi XingGang, Guangdong province) spiked with all analytes at three concentrations (10, 50, and 100 ng L−1), and each recovery test was conducted in triplicate. These three concentrations were comparable to those in the environment and also provided an appropriate concentration range for examining the robustness of the method.

2.3. Instrumental Analysis

All analytes were measured with a triple quadrupole 5500 electrospray ionization-tandem mass spectrometer (AB SCIEX; Redwood City, CA) equipped with an LC-20A high-performance liquid chromatograph (Shimadzu, Japan) in the multiple-reaction monitoring mode. Ion source gas, curtain gas, collision gas, temperature, and ionspray voltage were optimized by flow injection analysis tuning. The detailed parameters are shown in Table 1. The mass spectrometer parameters were optimized for each analyte by direct infusion of each standard solution in methanol at a concentration of 100 ng mL−1. The precursor and product ions, declustering potential, collision energy, entrance potential, and collision cell exit potential were optimized and confirmed to have the best instrument response through the liquid chromatography-mass spectrometry system (Tables 2 and 3). Two ion fragments were selected for qualitative and quantitative determination of each analyte to obtain high selectivity and sensitivity. The choice of positive and negative ion modes was based mainly on the physicochemical properties and instrumental responses of the analytes. The majority were analysed in positive ionisation mode, with the exception of the four steroids listed (diethylstilbestrol, estradiol, estrone, and estriol), which were analysed in negative ionisation mode. All analytes (except for diethylstilbestrol, estradiol, estrone, and estriol) were detected in the positive ion mode. The selected precursor ions for all the analytes in the positive mode were [M+H]+ while those in the negative mode were [M−H].

2.4. Site Description and Sample Collection

As one of the most economically developed regions in China, the Pearl River Delta (PRD) has experienced explosive economic growth over the past 30 years, leading to rapid industrialization and urbanization and resulting in a large population increase [41]. The living and consumption levels of residents have been greatly improved, which has brought about non-negligible impacts on the ecological environment of the region. The Pearl River Basin is a complex river basin composed of numerous main streams and tributaries, which are greatly affected by the surrounding human activities [42]. The eight major runoff outlets are imported into the South China Sea. Water samples were collected from the eight major runoff outlets at the PRD, located in Guangdong Province, Southern China, on August 11, 2018. The collected samples are odorless and contain a small amount of particles. The position of the sampling sites has been given by the previous study [41] (Supplemental Data; Table S2). The collected water samples were pretreated with a filter membrane within 48 hours. The water samples used for matrix spikes come from the park of Qi Xing Gang, which is located in Guangdong Province and has a target-free environmental pollution in the landscape protection zone. The water samples are odorless and free from interference.

2.5. Quantification and Method Validation

Solvent blanks, process blanks, and independent test standards (50 ng mL−1 standard solution) were run in sequence to check for residue, background contamination, and system performance. The reported quantitative values of each analyte in the samples were required to have the same retention time as its calibration standard and the same ion ratio. An independent test standard was inserted approximately every 10 injections to verify the stability of the instrument. If the concentration computed and the retention time of the instrument exceeded the standard value, the instrument was rebalanced and the sample quantified again to ensure the accuracy of the data results. All data processing was performed using Analyst Software Version 1.6.2. Limits of detection and limits of quantification of all analytes were calculated with signal/noise ratios of 3 and 10, respectively. The signal/noise ratios were obtained by averaging the results of three samples with the lowest spiked concentration [40]. The estimated limits of quantification values were used as the sample addition concentrations for seven environmental matrix spiked extracts to obtain the method detection limits (MDLs) and method quantitative limits (MQLs) [38]. The MDL of a target analyte was calculated according to the following equation:where and SD are, respectively, the t-distribution with a degree of freedom of n-1 and a confidence level of 99% and the standard deviation of the results of measurements.

In addition, matrix effects were estimated by the analyte concentration differences in environmental substrate sample extracts and solvents both spiked with 10 ng mL−1 of all analytes [37], i.e.,where and are chromatographic peak areas of matrix blank and standard solution, respectively. The matrix effects values of less or greater than 100% represent signal suppression or enhancement, respectively (Supplemental Data; Table S3).

3. Results and Discussion

3.1. Optimization of LC-MS/MS Conditions

The separation of 47 target analytes was optimized with an Agilent Eclipse Plus C18 column (100 mm × 2.1 mm, 1.8 μm), a Sigma-Aldrich C18 column (100 mm × 4.6 mm, 3 μm), and a Thermo Betasil C18 column (100 mm × 2.1 mm, 5 μm) with the mobile phase of MilliQ water containing 0.1% formic (mobile phase A) and acetonitrile (mobile phase B). When the flow rates of the 1.8 μm, 3 μm, and 5 μm columns were 0.25, 0.5, and 0.5 mL min−1, respectively, the analytes were well resolved and had sharp peak shapes. In terms of separation efficiency, the separation of 1.8 μm and 3 μm columns was better than that of 5 μm column. The three isomers (sulfameter, sulfamethoxypyridazine, and sulfamonomethoxine) were not well separated under 5 μm column conditions, and some of the analytes have a large half-width. Considering the chromatographic resolution and analysis time as well as reducing the high-pressure injection protection instrument, a comprehensive comparison of the final selection was a 3 μm column.

The addition of additives to the mobile phase could improve the peak shapes of the analytes and their responses on the instrument [39]. Some studies added formic acid and ammonium acetate to the mobile phase to increase the response of the antibiotics on the instrument, and the organic phase was selected from both methanol and acetonitrile [38, 39]. We compared the difference between methanol and acetonitrile in the mobile phase. In general, the system pressure of methanol was higher than that of acetonitrile. Although acetonitrile is not recommended as a toxic solvent, the higher system pressure was not conducive to the long-term stable operation of the column and instrument. It is reasonable to choose acetonitrile as the mobile phase. We also compared the effects of the proportion of additive formic acid (0.1%, 0.2%, and 0.5% in volume) on the analytes in the mobile phase A. The addition of formic acid has improved the peak shape of the analytes, but the increase in the ratio of formic acid does not significantly enhance the response of the analytes on the instrument. In contrast, we found that by adding 0.1% ammonia to mobile phase A in negative ion mode, the response of the analytes on the instrument significantly improved. In order to reduce the switching time of the positive and negative ion modes and the use of the mobile phase, a 5 μm column as the separation column for the negative ion mode was chosen. The optimized total elution time of the analytes is 17 and 8 min in positive and negative ion modes, respectively (Figures 1 and 2).

3.2. Optimization of Pretreatment Method Conditions

Studies have indicated that the addition of Na2EDTA to aqueous samples has increased the recoveries of MLs [43]. These studies [38, 39] also compared the effects of the quality of the addition of Na2EDTA on antibiotics. Na2EDTA (0.2 g) was added to promote the recovery of antibiotics. Extraction of all analytes in water samples was optimized by testing with the spiked samples at pH (3.0, 5.0, and 7.0). It showed that the recoveries of tilmicosin and erythromycin at pH 3.0 were higher than at other pH conditions, but lincomycin, clindamycin, testosterone, and methytestosterone at pH 7.0 were higher than at pH 3.0 (Figure 3). Shao et al. also verified that the extraction efficiency of solid-phase extraction cartridges with different pH values (3.0, 5.0, 7.0, and 9.0) indicated that the extraction efficiency of sulfonamides was not affected by pH [38]. But the extraction efficiency of tilmicosin and erythromycin was significantly better than other conditions at pH 3.0, which was similar to the optimization results in this study. The pH had slight effects on the spiked recoveries of OPs and NNIs. Although the spiked recoveries of some of the analytes were inhibited, in order to consider the low spiked recoveries of MLs and ensure acceptable recoveries of other analytes, a pH of 3.0 was finally selected. Therefore, conditions (0.2 g Na2EDTA and pH 3.0) were finally selected as they provided acceptable recoveries from the spiked samples.

3.3. Method Performance Validation

Method performance was assessed in terms of calibration linearity and range, recovery efficacy, matrix interferences, and method detection limits. Accurate quantification was performed using a seven-point calibration curve and showed a good linear relationship for analytes in the range of 0.1 to 100 ng mL−1. The linearity of all calibration curves was S > 0.99. The right linear range was selected for accurate quantification of environmental samples. If the quantitative concentration of a sample exceeded the linear range, a diluted sample was reinjected for accurate quantitation.

The recoveries of analytes were between 41% and 127%, indicating that the separation and purification of the samples met the requirements. The loss of analytes in the solid-phase extraction process is mainly due to the key factors such as the activation of the filler, the amount of adsorption, and the control of the loading flow rates. Quantitative use of the internal standard method minimizes the negative effects of matrix effects. Experimental analysis of water samples in the eight major runoff outlets of the PRD showed that some of the target analytes were detected, and the recoveries (average ± standard deviation) of sulfamethazine-d4, sulfamethoxypyridazine-d3, acetamiprid-d3, clothianidin-d3, and diethylstilbrstrol-d8 were 68 ± 7%, 79 ± 11%, 73 ± 5%, 95 ± 11%, and 71 ± 10%. The relative standard deviation was equal to or less than 20% except for individual analytes.

The matrix effects in LC-MS/MS are caused by the ionisation efficiency of the electrospray interface affected by the coeluting component of the analytes, which are manifested by ion enhancement or inhibition. It is usually corrected by the internal standard method or matrix standard solution to compensate for matrix effects. The matrix effects of most target analytes were in the range of 27–151%. The matrix effects were compensated by adding four suitable internal standards to the matrix. The matrix effects are effectively compensated by matrix-matched calibration curves and isotope dilution techniques. Due to the complexity and variability of the samples, it is not appropriate to have a matrix-matched calibration curve that takes into account economic issues.

MDLs are closely related to instrument sensitivity, concentration factor, and injection volume. This method has a concentration factor of 1000 and an injection volume of 3 μL with a Triple Quadrupole 5500 electrospray ionization-tandem mass spectrometer. We have calculated the MQLs of all analytes by (1). MQLs of the SAs, MLs, HMs, OPs, and NNIs were in the range of 0.53–1.49 ng L−1, 0.40–5.49 ng L−1, 0.40–3.16 ng L−1, 0.54–1.98 ng L−1, and 0.50–2.07 ng L−1.

3.4. Application of Environmental Samples

This method was used to analyze the environmental water samples collected at eight major runoff outlets of the PRD, South China, on August 11, 2018. The concentrations of target analytes in the environmental samples are presented in Table 4. The five sulfonamides detected in the eight major runoff outlets are also found in most rivers in China. And erythromycin is also present in many large rivers in China, but lincomycin and clindamycin, which have lower detection rates in most rivers, were detected in water samples. The concentration of parathion was as high as 154 ng L−1. The use of parathion was banned in 2003, but high levels of parathion were still detected in the water samples, indicating a high background in the environment and slow release through surface runoff and rain washout. The three neonicotinoids (acetamiprid, clothianidin, and imidacloprid), which are used in large amounts, were also presented in the water sample. It can be seen that the concentration of NNIs was generally higher than that of OPs, which may restrict the use of OPs and the replacement of NNIs as OPs in agricultural production. The average recoveries of the method were 68–95%, indicating that the method is simple and effective for actual environmental waters. This optimization method has been successfully applied to the analysis of the waters at eight major runoff outlets in the PRD. The frequent detection of various human veterinary drugs and pesticide residues indicated that the PRD has received pollution, and there may be certain ecological risks that should be given sufficient attention. A simple, economical, and powerful tool has been provided for other environmental waters.

Data Availability

The data used to support the findings of this study are included within the supplementary information file.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The present study was financially supported by the Youth Talent Innovation Foundation of Hubei Provincial Academy of Eco-Environmental Sciences (2021-15).

Supplementary Materials

Table S1: analyte details. Table S2: geographic information of sampling sites and sampling time. Table S3: recoveries (n = 3), method detection limits, and method quantitation limits of human and veterinary drugs and pesticide residues from surface water. (Supplementary Materials)