Abstract

Method validation is an essential technique for ensuring the reliability and accuracy of an analytical method. This study aimed to optimize and validate a fast, reliable, and accurate method for quantitatively determining pesticide residues of diverse chemical classes in the tomato matrix. Various method performance characteristics were tested and compared with predefined criteria. Twenty-six different pesticides of diverse chemical classes were selected based on their use in tomato cultivation and the availability of reference materials. The pesticide residues in tomato samples were extracted with the QuEChERS technique with some modifications, followed by injection into an LC-MS/MS system operating in an optimized method. The validated method demonstrated reasonable specificity, as there were no interferences from matrix components at the retention times of pesticides. The calibration curves for all pesticides exhibited excellent linearities, with correlation coefficients exceeding 0.99. No significant matrix effect was observed for all pesticides in tomatoes, as the values fell within the range of ±20%. All pesticides were quantified successfully at a concentration of 5 μg/kg except for carbaryl, with an average recovery of more than 70% and a relative standard deviation of less than 20%. Similarly, measurement uncertainties were also estimated based on the validation data, and the values were found below the default limit of 50%. Subsequently, the validated method was applied to analyze 52 locally collected tomato samples. Study findings revealed that only four of the studied pesticides were detected in these samples, and their concentrations were below the maximum residue limits (500 µg/kg each for carbendazim, imidacloprid, and metalaxyl) established for tomatoes by the Government of Nepal and the Codex Alimentarius Commission.

1. Introduction

Pesticides are widely used in agricultural production, including in the cultivation of tomatoes. Several groups of pesticides are used to manage weeds, insects, and other pests. However, some pesticides persist as residues in agricultural produce, and upon consumption of such contaminated produce, these pesticide residues enter the human body, posing potential adverse health impacts [1, 2]. Due to the hazardous nature of these pesticides, their presence is closely monitored in food to inhibit the potential adverse effects on human health. Consequently, many countries have established the maximum residue limits (MRLs) for the residues in various foods. Thus, monitoring the level of pesticide residues in food is desirable to ensure food safety.

Pesticide residue extraction from a food matrix is challenging due to its low concentration and potential interferences from the complex sample matrix [3]. Several sample extraction protocols are employed for this purpose [4]. However, the QuEChERS (Quick, Easy, Cheap, Effective, Rugged, and Safe) extraction technique is the obvious choice for multiresidue analysis due to its ease of use and fast extraction time [57]. Additionally, pesticides with a broad range of chemical nature can be extracted simultaneously. It involves two steps: sample extraction and sample cleanup. In the first step, residues are extracted from the homogeneous sample by adding acetonitrile and a blend of salts. The salt mixture allows the normally miscible organic solvent to separate from the water in the sample. In the second step, an aliquot of the organic phase from step 1 is taken for the cleanup through the use of dSPE (dispersive solid-phase extraction), where PSA (primary secondary amine) and anhydrous MgSO4 are used to remove residual water and many potential interfering substances such as sugar and organic acids from the extract. Based on the nature of the sample, other reagents, such as C18 and graphitized carbon black (GCB), can be added in the second step to facilitate the effective cleanup [811].

Several analytical techniques, such as gas chromatography-mass spectrometry, gas chromatography coupled to an electron capture detector, high-performance liquid chromatography, and liquid chromatography-tandem mass spectrometry, are available to monitor the level of pesticide residue in a food matrix [7, 12]. Among them, LC-MS/MS is one of the most widely used analytical techniques for pesticide residue analysis in food due to its ability to perform multiresidue analysis quickly with remarkable sensitivity. This study optimized an LC-MS/MS-based method, recognizing the necessity for a fast and reliable method for analyzing multi residue. Various method parameters were optimized to get an optimum response for each of the 26 selected pesticides. These pesticides belong to various chemical classes, including but not limited to carbamates, organophosphates, benzimidazoles, and neonicotinoids. The selection of pesticides in the study was based on the probable use in tomato cultivation and the availability of reference standards. Similarly, tomatoes were chosen as a matrix for this study due to their widespread consumption and availability of already established MRLs for pesticide residues. The Government of Nepal has set MRLs for 75 pesticides in tomatoes, the maximum number of pesticides regulated for any food category in the country. After optimization, the method was validated to assess its applicability. Method validation is an important technique that assesses an analytical method’s suitability for its intended purpose [13]. In validation, various performance criteria are examined to determine the fitness of purpose, and the results are compared against the predefined criteria. If all the tested parameters met the predefined acceptance criteria, then the method is considered fit for purpose. It is essential for ensuring the reliability of the test results. Although several method validation protocols are available depending on the nature of the analytical method, the widely used protocol for pesticide residue analysis is the SANTE guideline [14].

2. Materials and Methods

2.1. Chemicals and Reagents

Individual 16 pesticide standards were obtained from Sigma Aldrich, and a mixed standard of 10 carbamate pesticides was purchased from Restek. HPLC grade of acetonitrile, methanol, formic acid, ammonium formate, and acetic acid and analytical grade of anhydrous magnesium sulfate and primary secondary amine (PSA) were purchased from local suppliers. The individual stock solution other than the mixed standard solution was prepared at a concentration of 1000 mg/L by dissolving in an appropriate amount of methanol. All standard solutions were stored in a refrigerator at 4°C before use.

2.2. Instrumentation

The study was performed in an Agilent 1290 Infinity LC system connected to an Agilent 6460 triple quadrupole mass spectrometer equipped with Agilent Jet Stream electrospray ionization (AJS-ESI). The chromatographic and mass spectrometric data were acquired and analyzed by using MassHunter software. The chromatographic separation of pesticides was carried out in an Agilent Poroshell 120 EC-C18 analytical column having dimensions of 3.0 × 50 mm and 2.7 μm particle size.

2.3. Method Parameters

The method parameters were optimized before starting actual validation. Separation was carried out in a gradient mode with mobile phase A consisting of 0.1% formic acid and 5 mM ammonium formate in water. In contrast, mobile phase B comprised the same composition of formic acid and ammonium formate in methanol. The gradient started with 5% of mobile phase B, remaining constant until 0.5 minutes. Subsequently, it increased linearly to 65% at 5 minutes and rose to 95% at 6.5 minutes, maintaining the same composition until 9.0 minutes. At 9.1 minutes, the value of mobile phase B decreased sharply to 5% and remained the same until 12 minutes. The optimum flow rate, column temperature, and injection volume were 0.5 mL/min, 40°C, and 3 μL, respectively. Each analysis was completed in 12 minutes.

Pesticides were ionized in positive electrospray ionization (ESI) mode and acquired in dynamic multiple reaction monitoring (dMRM) mode. The optimized flow rates for drying gas and sheath gas were 10 L/min and 11 L/min, respectively, with their temperatures maintained at 250°C and 350°C, respectively. The nebulizer gas pressure was constant at 40 psi, while capillary and nozzle voltages were set at 4000 V and 300 V, respectively.

2.4. Sample Preparation

Pesticide residues incurred in tomato samples were extracted using the QuEChERS AOAC 2007.01 protocol [15] with some modifications. Homogenized samples underwent extraction using a solution comprising 1% acetic acid in acetonitrile. The phase separation of acetonitrile and water layers was achieved using a mixture of anhydrous magnesium sulfate and sodium acetate. The acetonitrile layer containing the extracted pesticides was cleaned using a primary secondary amine and anhydrous magnesium sulfate mixture. Generally, the final extract would undergo evaporation followed by reconstitution in an appropriate solvent that matches the mobile phase composition. However, in this study, the final extract was diluted with water in a 1 : 3 ratio, bypassing the time-consuming evaporation step.

2.5. Method Validation

Method validation is a process of demonstrating the fitness of a method for its intended purpose by examining and providing objective evidence [13]. The method was validated as per analytical quality control and method validation procedures for pesticides in food and feed [14]. Critical method performance parameters such as specificity, linearity, limit of quantification, trueness, and precision were rigorously assessed experimentally in a blank tomato matrix. The obtained data were compared against the predefined criteria outlined in the protocol.

2.6. Measurement Uncertainty

Every measurement is associated with a dispersion known as measurement uncertainty (MU). It gives a range of values within which the true value of the measured quantity is expected to lie. It reflects the inherent variability in a measurement process. In pesticide residue analysis, measurement uncertainty is critical during compliance statements against a standard. Measurement uncertainty values were estimated using a top-down approach based on the validation data [16, 17]. This approach uses trueness and precision data generated in the method validation experiment to estimate the MU value.

3. Results and Discussion

3.1. Optimization of Analyte-Dependent Mass Parameters

The mix working standard solution of all the pesticides was prepared at 1000 ng/mL and used to optimize their precursor ions, product ions, fragmentor, and collision voltages. The optimization work was performed by injecting pesticide standard solution without a column, and data were acquired using MassHunter Optimizer software. All pesticides showed good response in positive ESI mode, with the majority forming protonated ions [M + H]+, except for aldicarb, aldicarb sulfone, and propargite, which formed ammonium adducts [M + NH4]+ as a precursor ion. The optimized fragmentor voltages were found in the range of 50–145 V for all precursor ions, while collision energies ranged from 0 to 48 V for all product ions. Two product ions per precursor ion, their ion ratios, and retention time were utilized for pesticide confirmation [14]. The mass spectra and the ion ratio of carbaryl are shown in Figure 1. Initially, pesticides were monitored in multiple reaction monitoring (MRM) mode to obtain the retention times of each pesticide. Subsequently, it was upgraded to dynamic multiple reaction monitoring (dMRM) mode for maximizing sensitivity. In dMRM mode, specific MS transitions are monitored in a narrow time window corresponding to their expected elution from LC rather than throughout the full analysis time. The transition from MRM to dMRM drastically enhances the sensitivity of pesticides, as the MS duty cycle is not wasted by monitoring them when they are not expected to elute from LC [18]. The optimum acquisition parameters of all the target pesticides are summarized in Table 1.

3.2. Optimization of Chromatographic Parameters

Various mobile phase gradient programs were studied in search of optimal resolution and sensitivity. However, due to the diverse chemical nature of multiresidues, some peaks were coeluted even after numerous experiments. The extracted ion chromatogram of mixed standard pesticides in an optimum method (presented in Section 2.3) is shown in Figure 2.

3.3. Modification of Sample Preparation

This study followed the QuEChERS AOAC 2007.01 protocol [15] without deviation until the cleanup step. However, the solvent evaporation of the final extract was skipped to shorten the sample preparation time. The evaporation processes involving a nitrogen turbo evaporator are time consuming and involve the risk of degradation of certain pesticides if temperature is not maintained carefully. An alternative approach was employed to avoid the probable degradation of pesticides and reduce the lengthy evaporation time. In this approach, the final acetonitrile extract was diluted with water in a ratio of 1 : 3 instead of the final evaporation step. This dilution considerably shortened the sample preparation time while ensuring the solvent matching with the mobile phase, resulting in satisfactory peak shapes, especially for early eluting pesticides.

3.4. Method Validation
3.4.1. Specificity

The specificity of an analytical method is examined to verify the absence of potential interfering compounds at the retention time of the target analytes. The tomato blank matrix and matrix-matched standard solutions were analyzed simultaneously to assess the specificity of the method. Comparison of total ion chromatograms (TICs) of both samples, as shown in Figure 3, revealed the absence of any significant interfering peaks at the retention times of all target pesticides, thus indicating the specificity of the method.

3.4.2. Linearity

Linearity was evaluated from the calibration curve constructed from a series of eight duplicate concentrations ranging from 0.5 ng/mL to 100 ng/mL in a solvent with a final composition of 1 : 3 acetonitrile and water. The calibration curves were best fitted to a linear curve with weight 1/x. Good linear relationships were observed with regression coefficients (R2) of 0.99 or higher across the examined concentration range for all the pesticides. Calibration curves of some selected pesticides are presented in Figure 4, while regression coefficients (R2) and slopes of all the pesticides are summarized in Table 2.

3.4.3. Matrix Effect

The matrix effect is frequently encountered in LC-MS/MS analysis due to ion suppression or enhancement effect. It arises from coextracted compounds from the sample matrix, influencing analyte concentration measurement [19]. If not adequately compensated for, the matrix effect can significantly affect the trueness of analytical results. The matrix effect was evaluated by comparing the slopes of the solvent standard (SS) calibration curve with the slopes of the matrix-matched (MM) calibration curve, constructed using standards prepared in a blank sample matrix [20]. The percentage matrix effect was determined using the following formula [21]:

Matrix effect values between 0 and 20% are considered low matrix effects, and the correction for the matrix effect is not necessary. However, if more than a 20% matrix effect is observed, it should be addressed to ensure more accurate results [14]. The negative value represents ion suppression, whereas the positive value represents ion enhancement. The matrix effects for all the analyzed pesticides were found within ±20%, as shown in Table 2, indicating the absence of significant matrix effects caused by the tomato matrix. In such instances, solvent standards could be used for quantifying pesticide residues incurred in unknown samples. This is particularly useful in routine analysis when no blank tomato matrix is available for preparing matrix-matched calibration standards.

3.4.4. Limit of Quantification (LOQ)

The limit of quantification (LOQ) is a critical performance characteristic for pesticide residues. The value should be at or below the maximum residue level (MRL) to draw meaningful conclusions about pesticide residues incurred in a sample. LOQ can be determined using various methods [22, 23]. LOQ was determined by spiking a series of low concentrations in a blank tomato matrix. All pesticides, except carbaryl, were successfully analyzed at the 5 µg/kg level, meeting all identification, recovery, and precision criteria. For carbaryl, it was 10 µg/kg. These LOQs were much lower than the MRLs established for tomatoes by the Government of Nepal [24] and the Codex Alimentarius Commission [25]. The obtained values of LOQ are summarized in Table 2.

3.4.5. Trueness and Precision

Various approaches, such as recovery experiments and analysis of certified reference material, have been employed to assess an analytical method’s trueness [13, 26]. In this study, trueness was evaluated by recovery experiments, in which blank tomato samples were spiked at three different concentrations: low (5 µg/kg), medium (10 µg/kg), and high (40 µg/kg), each level in 6 replicates. The recoveries in the 70–120% range were obtained, exhibiting good trueness of the method. Similarly, precision was assessed as repeatability and within-laboratory reproducibility. Both were below 20% relative standard deviation (RSD) for all pesticides. The obtained values of trueness and precision are summarized in Table 3.

3.4.6. Measurement Uncertainty (MU)

Measurement uncertainty (MU) is inherently associated with any measurement. There are various approaches available for the estimation of measurement uncertainty. The most widely used approach is the ISO GUM approach [22, 27, 28]. However, estimating using the bottom-up approach for multiresidue is impractical. Thus, the present study estimated MUs using a top-down approach. Intralaboratory validation data were used to estimate the standard uncertainty and then expressed at a 95% confidence level. Two significant sources, precision and trueness (as bias), were taken as the main contributors to measurement uncertainty. The standard uncertainty of each pesticide was estimated using the following formula:

Uncertainty due to bias was estimated from the recovery experiment, while the percentage relative standard deviation of within-laboratory reproducibility was used to estimate the uncertainty arising from precision. The expanded uncertainty was calculated by multiplying the standard uncertainty with a coverage factor (k) of 2, which approximately gives a confidence level of 95%:

The obtained values for all the target pesticides, as presented in Figure 5, were lower than the 50% default value employed by many regulatory authorities for enforcement decisions.

3.5. Application of the Validated Method to the Real Sample

The validated method was further applied to 52 tomato samples collected from various local vegetable markets in Kathmandu, Nepal, for the simultaneous analysis of pesticide residues. The findings of the analyzed samples are shown in Figure 6. None of the target pesticides were detected in 10 (19.2%) samples. Pesticides were detected in 80.8% of the samples analyzed, ranging from one to four pesticides (carbendazim, imidacloprid, metalaxyl, and thiabendazole). Carbendazim, imidacloprid, metalaxyl, and thiabendazole were detected in 34, 28, 7, and 1 samples, respectively, either individually or in combination. Of 42 positive samples, 24 (57.1%) samples were found to contain at least one pesticide residue at or above the LOQ level, while in the remaining 18 (42.9%) samples, pesticides were found below the LOQ level for each detected pesticide. The concentrations of pesticide residues in analyzed tomato samples were found to be in the range of 8–268 µg/kg, with an average of 63 µg/kg for carbendazim, 5–125 µg/kg, with an average of 34 µg/kg for imidacloprid, and 8–67 µg/kg, with an average 32 µg/kg for metalaxyl. These pesticide residue concentrations were significantly lower than the maximum residue limit (MRL) set by the Codex Alimentarius Commission [25] and the Nepal Government [24] for tomatoes (500 µg/kg for each of the pesticides mentioned above). Similar results were observed in a previous study conducted by Bhandari et al. [29] in selected tomato samples (N = 32) from Nepal.

4. Conclusions

Pesticide residues in food have always been a public concern, with many people aware of harmful health effects caused by consuming food containing these residues. Analyzing the food samples for pesticide residues is essential to safeguard public health. Routine analysis for residues requires a quick and accurate method with high sensitivity. Thus, we optimized and validated an LC-MS/MS method to analyze 26 pesticide residues in tomato samples. All method validation performance characteristics were satisfactory, indicating the method’s reliability.

Furthermore, MUs were also estimated, which would be helpful while deciding the compliance of tomato samples against the established MRL. The validated method was successfully applied to analyze actual tomato samples, and some positive samples were detected and subsequently quantified. Even though 80.8% of the analyzed samples (N = 52) had one to four target pesticide residues, these amounts were significantly lower than the MRLs set for tomatoes. The validated method offers several advantages, including simple extraction, cleanup without an additional evaporation step of the final extract, high sensitivity, and good accuracy. Additionally, more pesticides can be included in the validated method to cover a broader range of pesticides. It offers an efficient solution for a quality control laboratory that analyzes multiresidues in tomato samples.

Data Availability

The experimental data used to support the findings of the study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Suraj Shrestha conceptualized and validated study, collected the sample, wrote the original draft, and finalized the manuscript. Bandana Lamichhane validated the study and collected the sample. Nibedita Chaudhary reviewed the draft. All the authors have approved the final manuscript.

Acknowledgments

The authors would like to acknowledge the Director General and Deputy Director Generals of the Department of Food Technology and Quality Control, Nepal, for constant support throughout the research.