Next Article in Journal
Enhanced Antimicrobial Activity of AgCu Nanoparticles: The Role of Particle Size and Alloy Composition
Previous Article in Journal
Inhibition of ACE2–S Protein Interaction by a Short Functional Peptide with a Boomerang Structure
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Support Vector Machine-Assisted Metabolomics Approach for Non-Targeted Screening of Multi-Class Pesticides and Veterinary Drugs in Maize

Technology Centre of Dalian Customs, Dalian 116000, China
*
Authors to whom correspondence should be addressed.
Molecules 2024, 29(13), 3026; https://doi.org/10.3390/molecules29133026
Submission received: 8 May 2024 / Revised: 15 June 2024 / Accepted: 19 June 2024 / Published: 26 June 2024

Abstract

:
The contamination risks of plant-derived foods due to the co-existence of pesticides and veterinary drugs (P&VDs) have not been fully understood. With an increasing number of unexpected P&VDs illegally added to foods, it is essential to develop a non-targeted screening method for P&VDs for their comprehensive risk assessment. In this study, a modified support vector machine (SVM)-assisted metabolomics approach by screening eligible variables to represent marker compounds of 124 multi-class P&VDs in maize was developed based on the results of high-performance liquid chromatography–tandem mass spectrometry. Principal component analysis and orthogonal partial least squares discriminant analysis indicate the existence of variables with obvious inter-group differences, which were further investigated by S-plot plots, permutation tests, and variable importance in projection to obtain eligible variables. Meanwhile, SVM recursive feature elimination under the radial basis function was employed to obtain the weight-squared values of all the variables ranging from large to small for the screening of eligible variables as well. Pairwise t-tests and fold changes of concentration were further employed to confirm these eligible variables to represent marker compounds. The results indicate that 120 out of 124 P&VDs can be identified by the SVM-assisted metabolomics method, while only 109 P&VDs can be found by the metabolomics method alone, implying that SVM can promote the screening accuracy of the metabolomics method. In addition, the method’s practicability was validated by the real contaminated maize samples, which provide a bright application prospect in non-targeted screening of contaminants. The limits of detection for 120 P&VDs in maize samples were calculated to be 0.3~1.5 µg/kg.

1. Introduction

It has always been a headache to resolve pesticide contamination in plant-derived foods, which poses a considerable threat to global food safety. Some national governments (e.g., China [1], the United States [2], and Japan [3]) and international authorities (e.g., the European Union [4]) have issued a series of formal regulatory documents on the maximum residue limits (MRLs) of pesticides in plant-derived foods. With the deepening of research, more and more evidence has shown that plant-derived foods are also suffering from serious contamination by veterinary drugs [5,6,7,8,9,10,11,12,13]. Previous studies [5,6,7,8,9,10,11,12,13] have shown that veterinary drugs excreted by livestock and poultry are commonly used as fertilizers on farmland, and some plant-derived foods (e.g., maize, potatoes, cucumbers, and lettuce) are easy to absorb from the soil, resulting in a variety of veterinary drugs (e.g., tetracycline, quinolones, and sulfonamides) accumulating in the foods, with a total concentration up to several mg/kg levels [8,14,15,16,17]. Due to the lack of risk evaluation standards for veterinary drugs in plant-derived foods at home and abroad, it is difficult to directly judge whether the concentration of veterinary drugs will cause adverse effects. Referring to the regulatory documents of MRLs of veterinary drugs in animal-derived foods [18,19], the concentration of 10 μg/kg was deemed the safety threshold for most veterinary drugs, and it is inferred that the reported concentration level of veterinary drugs in some plant-derived foods is likely to cause serious food safety incidents. For example, some organophosphorus pesticides can damage the nervous system, while sulfonamide antibiotics can cause allergic reactions after entering the human body. To sum up, resolving the food safety issue caused by the co-existence of pesticides and veterinary drugs (P&VDs) is urgent, and developing an effective screening method is a feasible way to prevent the occurrence of potential contamination risks.
Traditional methods to screen contaminants in foods are mostly based on the establishment of a set of databases containing a variety of P&VDs [20,21,22]. These methods have proved to be successful in discerning contaminants in the database, but they are helpless for contaminant identification outside the database, which casts a shadow on food safety. For example, some infamous food safety incidents took place in the 21st century, including melamine milk powder and fipronil eggs, which arose from some enterprises or individuals illegally adding unexpected contaminants not in the routine test scope of foods. Even if the concentrations of these contaminants were highly above the red line of food safety, these defective products were also considered to be qualified. As a result, they have gradually developed into extremely serious food safety incidents. It can be seen that the constant establishment of new databases of contaminants to screen unexpected P&VDs in plant-derived foods is a passive way that consumes a lot of human and material resources and is also unable to meet the increasing demands for the detection of unexpected P&VDs in foods. Therefore, it has become a popular trend to develop highly efficient non-targeted screening methods for contaminants, as the European Commission NORMAN network proposed [23,24].
Metabolomics has shown great promise in screening unexpected contaminants in the food safety field due to its advantages in handling massive data with complicated characteristics, such as small sample sizes, mass interferents, and high noise [25,26,27,28,29]. Marker compounds on behalf of unexpected contaminants in foods were screened and identified by in-house or network databases (e.g., Massbank, SciFinder, ChemSpider, PubChem, and Metlin) during metabolomics analysis in order to achieve non-targeted screening of contaminants [28,29].
As we know, metabolomics data are composed of fewer samples but massive features (i.e., variables). How to extract meaningful information and find the eligible variables on behalf of marker compounds is a primary issue for researchers to resolve. Generally speaking, the distinguishable factors for metabolomics data between any two groups (e.g., experiment vs. control) are a combination of variables rather than a single one. Therefore, several chemometric approaches, such as principal component analysis (PCA) [30,31,32], orthogonal partial least squares discriminant analysis (OPLS-DA) [30,31,32], random forest (RF) [33,34], and support vector machine (SVM) [33,35,36], have widely been used to deal with these metabolomics data. Among all these modeling methods, SVM is gaining popularity in a wide variety of metabolomics studies due to its prediction performance. SVM is known to have excellent generalization ability and can be applied to nonlinear cases with the assistance of kernels relative to PCA and OPLS-DA, with only an assumption of linearity. SVM recursive feature elimination (SVMRFE) is a wrapper approach that adopts the rule of weight (w) to rank the variables and works well when the number of samples is small but the number of variables is large. Due to this feature of SVMRFE, which is similar to that of metabolomics data, SVMRFE is often used in metabolomics studies [37,38,39,40]. In addition, SVMRFE has shown better predictability than OPLS-DA. To date, SVMRFE-involved metabolomics studies are mainly applied to seek eligible biomarkers for disease diagnosis, but they remain insufficient to complete the non-targeted screening of contaminants in foods. In this work, we try to apply SVM to metabolomics analysis to better screen and identify marker compounds on behalf of 124 multi-class P&VDs in maize so as to realize non-targeted screening of residual contaminants in plant-derived foods.

2. Results and Discussion

2.1. Data Preprocessing

The total ion chromatograms containing 124 P&VDs in three concentration groups are shown in Figure 1. The matrix complexity of maize may cause the recoveries of each contaminant in different samples to be significantly discrepant. As a result, no evident correlation between peak intensities and concentration levels was observed, as shown in Figure 1. To eliminate the peak intensity errors of variables arising from different recoveries of P&VDs during the pretreatment process, enrofloxacin-d5 (parent ion m/z 365.21092; fragment ions m/z 348.19692, 322.21804, and 246.11012; retention time 7.01 min) and atrazine-d5 (parent ion m/z 221.14017; fragment ions m/z 179.08602, 137.06393, and 101.08703; retention time 5.77 min) were jointly spiked for recovery calibration. Blank maize extract solutions were employed to prepare the standard curves (5, 10, 25, 50, and 100 ng/mL) of enrofloxacin-d5 and atrazine-d5 to calculate the recoveries of these two deuterated compounds in maize samples, with the recoveries to be 75.6%~93.1%, 72.9%~93.4%, and 74.9%~93.6% for enrofloxacin-d5 and 73.6%~95.2%, 81.2%~94.7%, and 72.9%~96.4% for atrazine-d5 in the 20, 50, and 100 ng/mL groups, respectively (Table S1, Supplementary Materials).
For each sample, we adopted the formula of 2 × 100%/((100%/recovery of enrofloxacin-d5) + (100%/recovery of atrazine-d5)) to calculate the final recovery of internal standards (Table S1), whose average value from the same concentration group needs to multiply a calibration coefficient to obtain 100% recovery. Peak intensities for internal standards and all other variables were also calibrated after multiplying the coefficient. After this, eligible variables with a relative standard deviation of peak intensity less than 30% in QC and three concentration groups were obtained to establish a 1447 × 39 data matrix for further analysis [41].

2.2. Multivariate Analysis

2.2.1. PCA Results

As indicated in Figure 2, all samples showed their credibility within a 95% confidence level. QC group samples have a close gathering, implying that the high quality of the data was available and deserved further analysis [42]. Nine samples from the same concentration group clustered closely, indicating the good intra-group similarity of these samples. Inter-group samples separated obviously on the first principal component axis, implying distinguishable differences in variables existed among these groups, which provided the possibility to hunt for marker compounds among concentration groups.

2.2.2. Cluster Analysis Results

Cluster analysis can directly describe the similarity level of samples, meaning that samples with the highest similarity gather together preferentially, level by level, until all the samples finish their aggregation. As observed in Figure S1 (Supplementary Materials), the samples from the same concentration group were obviously separated from others, showing their high intra-group congeniality. The evident inter-group separation can be attributed to the differences in peak intensities of variables. The findings of cluster analysis were consistent with those of PCA, and both supported the existence of variables with significant differences.

2.2.3. OPLS-DA Results

As we know, OPLS-DA is a two-class classification model; the size imbalance of two classes may introduce a bias in the computation of the decision rule, which can result in the unreliability of the OPLS-DA model [43]. To resolve this issue, we adopted the synthetic minority over-sampling technique (SMOTE) [44] to increase the sample size of the class with fewer samples. Herein, when the 20 ng/mL concentration group was designed as Class 1 (only 9 samples), and the 50 and 100 ng/mL concentration groups as a whole (18 samples in total) were designated Class 2, it needed to increase the size of Class 1 from 9 to 18 samples (Figure 3a). Similarly, the 100 ng/mL concentration group in Class 1 also needed to increase the class size to match that of the 20 and 50 ng/mL concentration groups in Class 2 (Figure 3b). As indicated in Figure 3, the obvious separation between Class 1 and Class 2 on the first principal component axis in OPLS-DA score plots was observed, implying that the differences in variables were significant between the two classes. R2Y and Q2 are two key parameters in the OPLS-DA model to evaluate the interpretability along the Y-axis and the predictability of the model, respectively [45,46]. When R2Y and Q2 values are close to 1, it means that the OPLS-DA model has high reliability and predictability. In Figure 3, the R2Y and Q2 values were both equal to 1, which favored the robustness of these OPLS-DA models.
A series of data points to represent the variables were shown in the S-plot plots (Figure 4), in which the points near the two tips of the ‘S’ plots signified the larger contribution and higher confidence levels of the corresponding variables in differentiating two classes of the OPLS-DA model. These variables were qualified to be marker compound candidates. Herein, marker compounds in the significantly low (20 ng/mL) and high (100 ng/mL) concentration groups should be sought at the right end of Figure 4a and the left end of Figure 4b, respectively.
As we know, the Y-axis in the permutation test plot is on behalf of R2Y and Q2Y for the model, and the X-axis represents the correlation coefficient between original and permuted response data. It is desirable to summarize the results of the permutation test in a quantitative manner. One way to do this is to conduct conventional regression analysis on the two sets of points, i.e., one regression line is fitted among the R2Y points (green circles) and another line among the Q2Y points (blue squares), as shown in Figure 5. The intercepts of the resulting regression lines are interpretable as measures of ‘background’ R2Y and Q2Y obtained by fit to random data. Experience shows that the R2Y-intercept should not exceed 0.3~0.4 and that the O2Y-intercept should not exceed 0.05. Intercepts below these limits indicate valid models. A 200-iteration permutation test was a common method to evaluate whether the OPLS-DA model underwent over-fitting by the significance test (p-value) of Q2Y metrics [47]. If p < 0.05, it meant that the OPLS-DA models were free of over-fitting, as reflected by the Q2Y-intercept being less than 0.05 [47]. Figure 5 clearly points out that all OPLS-DA models were immune to over-fitting under the principles mentioned above.
VIP, as a measure of the importance of the variable in the OPLS-DA model, was computed for each extracted variable. Usually, VIP > 1 was considered to be the crucial threshold to pick eligible variables on behalf of marker compounds. In this study, a total of 228 variables with VIP > 1 from 20 concentration groups vs. 50 and 100 concentration groups (20 vs. 50 and 100 in abbreviation applied below) were chosen to be the marker compound candidates, but only 188 variables were eligible from 100 concentration groups vs. 20 and 50 concentration groups (100 vs. 20 and 50 in abbreviation applied below). One hundred and thirty-nine variables were observed to be overlapped between two VIP lists.

2.2.4. SVM Results

After running SVM codes, all the classification accuracy of the 10-fold cross-validation was computed to be 100%, in favor of the models with good generalization ability. The best penalty factor (c) and best kernel function parameter (g) from 3-fold cross-validation at 100% accuracy rate were computed to be 0.0039 and 0.0039 for 20 vs. 50 and 100, and 0.0037 and 0.0037 for 100 vs. 20 and 50. Two hundred and twenty-eight variables with VIP > 1 obtained from 20 vs. 50 and 100 were included in the weight table of two hundred and fifty-three variables. Similarly, 188 variables with VIP > 1 corresponded to those in the weight table, which contained an additional 26 variables with VIP < 1 from 100 vs. 20 and 50. From the above, we can see that the overlap ratio of variables obtained between VIP and weight ranking methods was 90.1% (i.e., 228 × 100%/253) from 20 vs. 50 and 100, while this ratio fell to 87.8% (i.e., 188 × 100%/214) from 100 vs. 20 and 50, indicating that the variables sought by the two methods were of high consistency, which laid a solid foundation to accurately seek eligible variables to represent marker compounds. One hundred and fifty-four overlapped variables from the two weight lists mentioned above were obtained for univariate analysis to further confirm the validity of those variables.

2.3. Univariate Analysis

Pairwise t-tests [47,48,49] and fold changes of concentration [27,50] have proved to be appropriate for univariate analysis of the metabolomics data. A pairwise t-test focuses on the significant differences of variables in peak intensity between two concentration groups, while a fold change of concentration considers the peak intensity ratio of variables between the high and low concentration groups. In general, the variables with a significance level (p) of pairwise t-tests below 0.05 and fold change (FC) of concentration above 2 are deemed eligible. In this study, p values of 154 variables between 20 and 50, 20 and 100, as well as 50 and 100 ng/mL concentration groups, i.e., p20vs.50, p20vs.100, and p50vs.100, were calculated to be lower than 0.05 in a pairwise t-test (Table S2), indicating that the significant inter-group differences of 154 variables indeed existed. In addition, the fold changes of concentration between 50 and 20 (FC50vs.20), as well as 100 and 20 ng/mL groups (FC100vs.20), were also calculated to be 2.333~2.671 and 4.555~5.332 (Table S3), which were above 2 to further support the eligibility of 154 variables as marker compound candidates.
As shown in Table 1, 120 out of 154 variables were confirmed as marker compounds (124 in total) by the weight ranking method, with a screening rate of over 96% (i.e., 120 × 100%/124), while the VIP method can only find 109 marker compounds with a screening rate of about 88% (i.e., 109 × 100%/124), showing that the weight ranking method expressed better performance in selecting eligible variables to represent marker compounds than the VIP method. Fleroxacin, lincomycin, clindamycin, and mebendazole, as four veterinary drugs, have no variables equal to their identities as marker compounds in weight and VIP lists. To be different with these four contaminants, three pesticides (i.e., pyridaben, hexazinone, and phosmet) and another eight veterinary drugs, including sulfaphenazole, sulfamethoxazole, sulfamethoxypyridazine, lomefloxacin, fenthion, danofloxacin, sparfloxacin, and tilmicosin, played the role of unique marker compounds, which corresponded to the variables with VIP < 1 from 20 vs. 50 and 100, but were also included in the weight lists (Table 1). Thirty-four in one hundred and fifty-four variables not identified as marker compounds in weight lists showed more complicated and erratic LC–MS/MS and metabolomics information, which needed more in-depth studies to discern their origins from a matrix or foreign impurity.
A blank maize sample was also investigated by the SVM-assisted metabolomics method with no corresponding variables of 124 P&VDs found, eliminating the inherent (rather than spiked) interference of these contaminants in maize to screen marker compounds and the existence of false positives.

2.4. Limits of Detection

According to the method proposed by the US Environmental Protection Agency to calculate the limits of detection (LODs) of 120 P&VDs [51], we first prepared a 2.0 g blank maize sample for the same pretreatment as introduced above to obtain a 1 mL extract solution. In this extract, we spiked a mixed solution (1 µg/mL, 20 µL) containing 120 P&VDs, resulting in a final concentration of 20 ng/mL. This process was repeated to gain seven replicates, which were also subjected to metabolomics analysis to obtain peak intensities of 120 P&VDs. For each P&VD, the concentration level of 20 ng/mL corresponds to the peak intensity average value of seven replicates. Therefore, the concentration (ng/mL) of each P&VD can be calculated by the formula for the respective peak intensity × 20/average peak intensity. As shown in Table 1, the LOD results for 120 P&VDs varied between 0.3 and 1.5 µg/kg.
The metabolomics method has successfully been applied to non-targeted contaminant screening in plant-derived foods such as tea [27], lettuce [29], maize [29], and orange juice [33]. In contrast to these studies, which only focused on pesticides or veterinary drugs, our study expands the domain of target compounds to simultaneously screen multi-class P&VDs. This is the first time to prove the introduction of SVM as an effective tool for promoting the screening rate of the metabolomics approach. In addition, method LODs were calculated to be significantly below 10 µg/kg for 120 P&VDs, which even have the upper hand over some targeted multi-residue methods [20,21,22].

2.5. Practicability Test

For plant-derived foods, the contamination pathways of pesticides are almost explicit; however, the contamination sources of veterinary drugs have not been fully understood. As reported in previous studies [52,53,54], the application of fertilizer and irrigation water contaminated with veterinary drugs onto farmland is a primary reason to cause the contamination of plant-derived foods; we hypothesized that some plant-derived foods from rural areas exposed to the abuse of veterinary drugs in livestock farming may suffer higher contamination risks. To investigate the practicability of our proposed SVM-assisted metabolomics screening method, we adopted maize samples from Qili and Bali villages as a case study. These two villages in Jinpu New Area, affiliated with Dalian City, are not only crucial maize suppliers but also pivotal providers of fresh livestock and poultry meat, which makes the maize samples highly susceptible to contamination by veterinary drugs.
Maize samples were collected after the maturity period in late September 2023. After metabolomics analysis, four contaminants, including norfloxacin, enrofloxacin, imidacloprid, and carbendazim, were all detected in two villages, with their concentrations being 10.7~17.7 µg/kg (Table S4). The abundance of norfloxacin and enrofloxacin in maize samples can be attributed to their extensive use as medicines and feed additives in animal husbandry [55], while the residual imidacloprid and carbendazim in maize samples with relatively high concentrations can be due to their widespread application to combat aphid [56] and bacterial wilt [56], respectively. As indicated in Table S4, the contamination of pesticides remained predominant over that of veterinary drugs, and imidacloprid presented the highest concentration among the four contaminants in both villages. In addition, Qili village suffered more serious overall contamination from four contaminants than Bali village. The test results supported the practicability of our proposed SVM-assisted metabolomics method in non-targeted screening of contaminants. It is worth noting that all four contaminants identified by the method presented a concentration above 10 µg/kg, which was understandable to show the feasibility of the method. But for those P&VDs with concentrations below 10 µg/kg, it is still unclear how to prove the applicability of the method, which probably led to those contaminants being unable to be detected. That is to say, the actual contamination of maize by P&VDs in the Jinpu New Area can be even worse than what we expected, and it deserves more attention to defuse the underlying risks. Food safety has no borders. In view of some food issues that initially occurred in local areas on a small scale, it is still essential for any country and international organization to come to notice and to strengthen cooperation in the prevention and control of the issues, such as through information sharing, experience exchange, and even the publication of standardized documents to ensure food safety from the source.

3. Materials and Methods

3.1. Chemicals and Materials

Maize samples were purchased from a local farmer’s market in Dalian City, China. Methanol and acetonitrile (HPLC grade) were purchased from Merck Corporation (Darmstadt, Germany). Formic acid and acetic acid (HPLC grade) were provided by China Shanghai ANPEL Laboratory Technologies Inc. (Shanghai, China). Ammonium formate, anhydrous sodium sulfate (Na2SO4), anhydrous sodium acetate (NaAC), and primary secondary amine (PSA) sorbent were obtained from China Sinopharm Chemical Reagent Co., Ltd. (Shanghai, China). Filter membrane (0.22 µm) was purchased from Agilent Technologies (Santa Clara, CA, USA). Ultrapure water was produced by the Milli-Q ultrapure water system from Merck Corporation (Darmstadt, Germany). Enrofloxacin-d5 and atrazine-d5 (100 μg/mL in methanol) for recovery calibration were provided by First Standard (Worcester, MA, USA). Sixty-nine pesticides and fifty-five veterinary drugs were purchased from the corporations First Standard (Worcester, MA, USA), Sigma (St. Louis, MO, USA), Dr. Ehrenstorfer (Wesel, Germany), and TRC (Montréal, QC, Canada), with a purity greater than 98%. The detailed information on these P&VDs is shown in Table 2.

3.2. Solution Preparation

A total of 124 P&VDs was separately prepared in methanol with a concentration of 100 µg/mL as the stock solution, 1 mL of which was taken, mixed together, and diluted with methanol to prepare a 1 µg/mL working solution. A 100 ng/mL mixed working solution of enrofloxacin-d5 and atrazine-d5 was prepared by diluting their 100 µg/mL stock solution with methanol.

3.3. Sample Preparation and Pretreatment Process

(a) The maize sample was ground into powder by a grinder; (b) 2.0, 5.0, and 10.0 g of maize powder, together with corresponding 20, 50, and 100 µL of 124 P&VDs mixed solutions (1 µg/mL), were poured into 50 mL polypropylene centrifuge tubes, respectively. To calibrate the recovery during the sample pretreatment process, a mixed solution (0.5 mL, 100 ng/mL) of enrofloxacin-d5 and atrazine-d5 was further added as recovery internal standards; (c) a total of 20 mL of acetonitrile/water (80/20, v/v) with 1% acetic acid was dumped into the tube by vortex and ultrasonic extraction for 1 and 30 min, respectively. A total of 6.0 g anhydrous Na2SO4 and 1.5 g NaAC were added, shocking violently for 1 min and centrifuging at 6000 r/min for 5 min; (d) all the supernatant solutions were taken into a new tube containing 2.0 g anhydrous Na2SO4 and 0.3 g PSA, centrifuging at 6000 r/min for 2 min; (e) all the solutions were extracted, dried with nitrogen gas flow, and redissolved with 1 mL 40% (v/v) methanol–4 mmol/L ammonium formate buffer solution, vortexing for 1 min; (f) filtered with a 0.22 µm filter membrane, the sample solutions of 124 P&VDs at the theoretical concentrations of 20, 50, and 100 ng/mL were prepared. Each concentration experiment was completed in nonuplicate to obtain replicate samples.

3.4. Sample Grouping and Naming

Twenty-seven samples named 20 ng/mL−1~20 ng/mL−9, 50 ng/mL−1~50 ng/mL−9, and 100 ng/mL−1~100 ng/mL−9 were obtained from three concentration groups. A quality control (QC) sample was prepared by mixing 30 µL of each sample from the above-mentioned concentration groups [27,57,58] and underwent three repeated injections into the LC–MS/MS system before and after the injection of each concentration group. Twelve injections of QC samples were named after QC-1, QC-2, …, and QC-12 to evaluate the stability of LC–MS/MS.

3.5. Analytical Method

A quadrupole/electrostatic field orbitrap LC–MS/MS system (Q Exactive Plus, Thermo Fisher Scientific Inc., Waltham, MA, USA) equipped with a heatable electrospray ion (ESI) source was used to analyze all pesticides, veterinary drugs, enrofloxacin-d5, and atrazine-d5. The ESI-positive mode was adopted to obtain M+H adducts for metabolomics analysis. An Accucore RP-MS column (100 mm × 2.1 mm, 2.6 µm particle diameter) purchased from Thermo Fisher Corporation (Waltham, MA, USA) was employed to separate components after an injection of 10 µL sample solution. The mobile phases were composed of 0.1% (v/v) formic acid and ammonium formate at 4 mmol/L in water (eluent A) and 0.1% (v/v) formic acid and ammonium formate at 4 mmol/L in methanol (eluent B). The flow rate was kept at 0.3 mL/min. The elution program was as follows: 100% A (0 min), 100% A (2 min), 0% A (8 min), 0% A (13 min), and 100% A (14 min) until the end of the run. The oven temperature was kept at 40 °C. The Q Exactive Plus parameter settings include (a) heating and capillary temperature at 320 °C; (b) sheath and auxiliary gas from N2, with flow rates of 40 and 10 arb, respectively; (c) lens and spray voltage at 50 and 3200 V, respectively; (d) scan mode to be full-scan/data-dependent two-stage scanning; (e) MS parameters to be full-scan resolution at 70,000, AGC target at 1 × 106, maximum dwell time at 100 ms, and scan range of m/z between 100 and 1000; (f) MS/MS parameters to be resolution at 17,500, AGC target at 2 × 105, and maximum dwell time at 50 ms.
Trace Finder software (Version 3.3) installed on the LC–MS/MS system was employed for residue analysis, with specific screening conditions as follows: (a) for primary parent ions, the response intensity threshold at 10,000, the ratio of signal to noise at 5.0, and mass error at 5 ppm; (b) for secondary fragment ions, the minimum number of matching ions at 1, the response intensity threshold at 10,000, and mass error at 5 ppm. The quantification of enrofloxacin-d5 and atrazine-d5 was finished by standard curves on the basis of the peak areas of primary parent ions for recovery calibration.

3.6. Metabolomics Data Processing

LC–MS/MS output files cannot be directly employed for metabolomics analysis and need to be converted from .RAW format to .mzXML format [59]. The new-formatted files were uploaded onto the Workflow4Metabolomics (W4M) platform (https://workflow4metabolomics.usegalaxy.fr/, accessed on 19 October 2022) for further analysis [60]. Chromatographic peak detection, alignment, and retention time calibration were performed on the platform, with further operations including normalization, centralizing, scaling, and data transformation of peak intensity to acquire the data matrix, in which variable and sample names were designated abscissa and ordinate, respectively [60,61]. PCA [62,63,64] and OPLS-DA [65,66] embedded in SIMCA software (Version 14.1) [67], together with cluster analysis in Heml software (Version 1.0.3.7), were performed for multivariate analysis of the data matrix. Over-fitting of the OPLS-DA model was evaluated by permutation tests [67,68]. Variable importance in projection (VIP) > 1 is the threshold to screen variables as marker compound candidates [67,68,69]. Under this principle, eligible variables from three concentration groups were included in the VIP lists.
SVM can perform nonlinear classification by mapping the input data into a high-dimensional feature space, the process of which is known to be the kernel trick [70]. Many kernel functions are available, but the most commonly used one is the radial basis function (RBF). Therefore, we chose SVMRFE under RBF to run the program in MATLAB software (Version R2019b). A 10-fold cross-validation was performed by randomly dividing the dataset into 10 subsets and computing the mean value of 10 accuracy values. This process was repeated 10 times, and the mean accuracy was obtained as the final classification accuracy. Kernel parameters were chosen through a 3-fold cross-validation approach, which divided the training dataset randomly into three subsets. The classifier was trained on either of the two subsets and tested on the third one. A set of parameters that provided the best cross-validation accuracy was employed for further analysis using the LIBSVM interface [71]. The kernel parameters with the best cross-validation accuracy were chosen to perform the actual classification. Under the conditions described above, weight values (w) and weight squared values (w2) of variables were obtained, with the latter ones ranging from large to small to form a variable weight table for further analysis. The MATLAB codes used to generate the results were provided in the Supplementary Materials.
All VIP > 1 variables filled their corresponding vacancies in the variable weight table. For those variables, their sequence in the weight table was sequential, implying that two methods, including VIP and weight ranking, to seek eligible variables have completely consistent findings. But for those variables, their VIP > 1 sequence in the weight table was not consecutive, implying some variables with VIP < 1 were also listed in the weight table. As mentioned above, SVMRFE has higher predictability than OPLS-DA in general, i.e., the weight ranking method is more reliable than the VIP method in picking eligible variables. Therefore, we still kept those variables with VIP < 1 in the weight table to prevent the possible loss of some valid variables during VIP computation for further analysis. Overlapped variables in 20 and 100 ng/mL groups, on behalf of their significantly low and high concentrations, underwent further verification by pairwise t-test [47,48,49] performed in SPSS Statistics software (Version 17.0) and fold change calculation of concentration during univariate analysis. Finally, eligible variables were confirmed to represent the marker compounds by comparing the retention time and precise molecular weight (absolute value of error less than 5 ppm) of 124 P&VDs (Table 2).

4. Conclusions

This study developed an SVM-assisted metabolomics method to screen the marker compounds of 124 P&VDs in maize samples. One hundred and twenty out of one hundred and twenty-four P&VDs can be identified by our proposed method, while only one hundred and nine P&VDs can be found by the metabolomics method alone, implying that SVM can promote the screening accuracy of the metabolomics method. It is the first time to apply the SVM-assisted metabolomics method to non-targeted screening of P&VDs in plant-derived foods. Method LODs were calculated for 120 P&VDs, with values significantly below 10 µg/kg for all of them, which were favorably comparable with a targeted multi-residue method. Our approach, developed on a simple and self-designed case, has been validated on a more complicated and realistic condition. This study promoted insight into the development of truly non-targeted screening approaches on the basis of metabolomics (particularly LC–MS/MS and chemometrics) for food safety assessment, which may be in need of non-targeted methods as a complement to targeted methods in view of screening potentially contaminated food products in the near future. Since our proposed method was rather generic, it was applicable with only a few modifications to any other LC–MS/MS dataset or even to other fields like authenticity or origin issues to complement existing methods. In conclusion, there are reasons to believe that these findings will encourage more breakthroughs in analytical issues, both in the methods and the tools.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules29133026/s1, MATLAB codes to run SVM; Table S1: Recovery of internal standards in maize sample groups; Table S2: Pairwise t-test results of 154 variables; Table S3: Fold change results of 154 variables; Table S4: Concentration of P&VDs in maize samples from Jinpu New Area. Figure S1. Cluster analysis plot of spiked maize sample groups.

Author Contributions

Conceptualization, W.X.; data curation, W.X.; investigation, W.X., F.L. and X.L.; methodology, W.X.; writing—original draft, W.X.; writing—review and editing, W.X. and Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Natural Science Foundation of Liaoning Province of China (No. 2023-MS-349) and the Science and Technology Project of the General Administration of Customs of the PRC (No. 2023HK124).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in this article and Supplementary Materials.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. GB 2763-2021; National Food Safety Standard-Maximum Residue Limits for Pesticides in Food. Ministry of Agriculture and Rural Affairs of the People’s Republic of China (MARAPRC): Beijing, China, 2021.
  2. National Archives and Records Administration (NARA). Available online: https://www.ecfr.gov/current/title-40/chapter-I/subchapter-E/part-180 (accessed on 15 October 2022).
  3. The Japan Food Chemical Research Foundation (JFCRF). Available online: http://db.ffcr.or.jp/front/ (accessed on 15 October 2022).
  4. European Union (EU). Commission Regulation (EU) No 396/2005 of 23 February 2005. 2005. Available online: https://eur-lex.europa.eu/eli/reg/2005/396/oj (accessed on 15 October 2022).
  5. Ahmed, M.B.M.; Rajapaksha, A.U.; Lim, J.E.; Vu, N.T.; Kim, S.; Kang, H.M.; Lee, S.S.; Ok, Y.S. Distribution and accumulative pattern of tetracyclines and sulfonamides in edible vegetables of cucumber, tomato and lettuce. J. Agric. Food Chem. 2015, 63, 398–405. [Google Scholar] [CrossRef] [PubMed]
  6. Boxall, A.B.A.; Johnson, P.; Smith, E.J.; Sinclair, C.J.; Stutt, E.; Levy, L.S. Uptake of veterinary medicines from soils into plants. J. Agric. Food Chem. 2006, 54, 2288–2297. [Google Scholar] [CrossRef] [PubMed]
  7. Hawker, D.W.; Cropp, R.; Boonsaner, M. Uptake of zwitterionic antibiotics by rice (Oryza sativa L.) in contaminated soil. J. Hazard. Mater. 2013, 263, 458–466. [Google Scholar] [CrossRef] [PubMed]
  8. Hu, X.; Zhou, Q.; Luo, Y. Occurrence and source analysis of typical veterinary antibiotics in manure, soil, vegetables and groundwater from organic vegetable bases, northern China. Environ. Pollut. 2010, 158, 2992–2998. [Google Scholar] [CrossRef] [PubMed]
  9. Jones-Lepp, T.L.; Sanchez, C.A.; Moy, T.; Kazemi, R. Method development and application to determine potential plant uptake of antibiotics and other drugs in irrigated crop production systems. J. Agric. Food Chem. 2010, 58, 11568–11573. [Google Scholar] [CrossRef] [PubMed]
  10. Tanoue, R.; Sato, Y.; Motoyama, M.; Nakagawa, S.; Shinohara, R.; Nomiyama, K. Plant uptake of pharmaceutical chemicals detected in recycled organic manure and reclaimed wastewater. J. Agric. Food Chem. 2012, 60, 10203–10211. [Google Scholar] [CrossRef] [PubMed]
  11. Tasho, R.P.; Cho, J.Y. Veterinary antibiotics in animal waste, its distribution in soil and uptake by plants: A review. Sci. Total Environ. 2016, 563–564, 366–376. [Google Scholar] [CrossRef] [PubMed]
  12. Wang, C.; Luo, Y.; Mao, D. Sources, fate, ecological risks and mitigation strategies of antibiotics in the soil environment. Environ. Chem. 2014, 33, 19–29. [Google Scholar]
  13. Zhang, M.; Xu, Q. Review on pollution and behavior of veterinary antibiotics in agricultural system. Acta Agric. Zhejiang 2013, 25, 416–424. (In Chinese) [Google Scholar]
  14. Bao, Y.; Li, Y.; Mo, C.; Yao, Y.; Tai, Y.; Wu, X.; Zhang, Y. Determination of six sulfonamide antibiotics in vegetables by solid phase extraction and high performance liquid chromatography. Environ. Chem. 2010, 29, 513–518. (In Chinese) [Google Scholar]
  15. Wu, X.; Xiang, L.; Mo, C.; Li, Y.; Jiang, Y.; Yan, Q.; Lv, X.; Huang, X. Determination of quinolones in vegetables using ultra performance liquid chromatography-electrospray ionization tandem mass spectrometry. Chin. J. Anal. Chem. 2013, 41, 876–881. (In Chinese) [Google Scholar] [CrossRef]
  16. Wu, X.; Bao, Y.; Xiang, L.; Yan, Q.; Jiang, Y.; Li, Y.; Huang, X.; Mo, C. Simultaneous sulfonamide antibiotics in vegetables using solid phase extraction and high performance liquid chromatography coupled with mass spectrometry. Environ. Chem. 2013, 32, 1038–1044. (In Chinese) [Google Scholar]
  17. Yao, Y.; Mo, C.; Li, Y.; Bao, Y.; Tai, Y.; Wu, X. Determination of tetracyclines in vegetables using solid phase extraction and HPLC with fluorescence detection. Environ. Chem. 2010, 29, 536–541. (In Chinese) [Google Scholar]
  18. European Union (EU). Commission Regulation (EU) No 37/2010 of 22 December 2009. 2010. Available online: https://eur-lex.europa.eu/legal-content/EN/TXT/?uri=CELEX%3A32010R0037 (accessed on 15 October 2022).
  19. GB 31650-2019; National Food Safety Standard-Maximum Residue Limits for Veterinary Drugs in Foods. Ministry of Agriculture and Rural Affairs of the People’s Republic of China (MARAPRC): Beijing, China, 2019.
  20. Jadhav, M.R.; Pudale, A.; Raut, P.; Utture, S.; Shabeer, T.P.A.; Banerjee, K. A unified approach for high-throughput quantitative analysis of the residues of multi-class veterinary drugs and pesticides in bovine milk using LC-MS/MS and GC-MS/MS. Food Chem. 2019, 272, 292–305. [Google Scholar] [CrossRef] [PubMed]
  21. Jia, W.; Chu, X.; Ling, Y.; Huang, J.; Chang, J. High-throughput screening of pesticide and veterinary drug residues in baby food by liquid chromatography coupled to quadrupole Orbitrap mass spectrometry. J. Chromatogr. A 2014, 1347, 122–128. [Google Scholar] [CrossRef] [PubMed]
  22. Gao, F.; Zhao, Y.; Shao, B.; Zhang, J. Determination of residues of pesticides and veterinary drugs in milk by ultra performance liquid chromatography coupled with quadrupole-time of flight mass spectrometry. Chin. J. Chromatogr. 2012, 30, 560–567. (In Chinese) [Google Scholar] [CrossRef] [PubMed]
  23. Brack, W.; Dulio, V.; Slobodnik, J. The NORMAN network and its activities on emerging environmental substances with a focus on effect-directed analysis of complex environmental contamination. Environ. Sci. Eur. 2012, 24, 29–33. [Google Scholar] [CrossRef]
  24. Fu, Y.; Zhou, Z.; Kong, H.; Lu, X.; Zhao, X.; Chen, Y.; Chen, J.; Wu, Z.; Xu, Z.; Zhao, C.; et al. Nontargeted screening method for illegal additives based on ultrahigh-performance liquid chromatography-high-resolution mass spectrometry. Anal. Chem. 2016, 88, 8870–8877. [Google Scholar] [CrossRef] [PubMed]
  25. Kunzelmann, M.; Winter, M.; Aberg, M.; Hellenäs, K.E.; Rosén, J. Non-targeted analysis of unexpected food contaminants using LC-HRMS. Anal. Bioanal. Chem. 2018, 410, 5593–5602. [Google Scholar] [CrossRef]
  26. Tengstrand, E.; Rosén, J.; Hellenäs, K.E.; Aberg, K.M. A concept study on non-targeted screening for chemical contaminants in food using liquid chromatography-mass spectrometry in combination with a metabolomics approach. Anal. Bioanal. Chem. 2013, 405, 1237–1243. [Google Scholar] [CrossRef]
  27. Delaporte, G.; Cladière, M.; Jouan-Rimbaud Bouveresse, D.; Camel, V. Untargeted food contaminant detection using UHPLC-HRMS combined with multivariate analysis: Feasibility study on tea. Food Chem. 2019, 277, 54–62. [Google Scholar] [CrossRef] [PubMed]
  28. Xue, W.; Zhang, H.; Wang, M.; Liu, Y.; Liu, M.; Shen, B. Metabolomics-based non-targeted screening analysis of 34 PPCPs in bovine and piscine muscles. Anal. Methods 2022, 14, 233–240. [Google Scholar] [CrossRef] [PubMed]
  29. Xue, W.; Yang, C.; Liu, M.; Lin, X.; Wang, M.; Wang, X. Metabolomics approach on non-targeted screening of 50 PPCPs in lettuce and maize. Molecules 2022, 27, 4711. [Google Scholar] [CrossRef] [PubMed]
  30. Rivera-Pérez, A.; Romero-González, R.; Frenich, A.G. A metabolomics approach based on 1H NMR fingerprinting and chemometrics for quality control and geographical discrimination of black pepper. J. Food Compos. Anal. 2022, 105, 104235–104244. [Google Scholar] [CrossRef]
  31. Azizan, K.A.; Zamani, A.I.; Muhammad, N.A.N.; Khairudin, K.; Yusoff, N.; Nawawi, M.F. Dose-dependent effect of wedelia trilobata essential oil (EO) on lettuce (Lactuca sativa L.) with multivariate analysis. Chem. Biodivers. 2022, 19, e202100833. [Google Scholar] [CrossRef] [PubMed]
  32. Cui, Y.; Ge, L.; Lu, W.; Wang, S.; Li, Y.; Wang, H.; Huang, M.; Xie, H.; Liao, J.; Tao, Y. Real-time profiling and distinction of lipids from different mammalian milks using rapid evaporative ionization mass spectrometry combined with chemometric analysis. J. Agric. Food Chem. 2022, 70, 7786–7795. [Google Scholar] [CrossRef] [PubMed]
  33. Sun, R.; Xing, R.; Zhang, J.; Yu, N.; Ge, Y.; Zhang, W.; Chen, Y. UPLC-QTOF-MS coupled with machine learning to discriminate between NFC and FC orange juice. Food Control. 2023, 145, 109487. [Google Scholar] [CrossRef]
  34. Yang, S.; Yang, H.; Huang, Y.; Chen, G.; Shen, C.; Zheng, S. Serum metabolomic signatures of hirschsprung’s disease based on GC–MS and LC–MS. J. Proteome Res. 2023, 22, 2307–2318. [Google Scholar] [CrossRef] [PubMed]
  35. Zheng, J.; Zheng, Y.; Li, W.; Zhi, J.; Huang, X.; Zhu, W.; Liu, Z.; Gong, L. Combined metabolomics with transcriptomics reveals potential plasma biomarkers correlated with non-small-cell lung cancer proliferation through the Akt pathway. Clin. Chim. Acta 2022, 530, 66–73. [Google Scholar] [CrossRef]
  36. Lu, Y.; Yao, G.; Wang, X.; Zhang, Y.; Zhao, J.; Yu, Y.J.; Wang, H. Chemometric discrimination of the geographical origin of licorice in China by untargeted metabolomics. Food Chem. 2022, 380, 132235. [Google Scholar] [CrossRef]
  37. Lin, X.; Yang, F.; Zhou, L.; Yin, P.; Kong, H.; Xing, W.; Lu, X.; Jia, L.; Wang, Q.; Xu, G. A support vector machine-recursive feature elimination feature selection method based on artificial contrast variables and mutual information. J. Chromatogr. B 2012, 910, 149–155. [Google Scholar] [CrossRef]
  38. Huang, X.; Lin, X.; Zeng, J.; Wang, L.; Yin, P.; Zhou, L.; Hu, C.; Yao, W. A computational method of defining potential biomarkers based on differential sub-networks. Sci. Rep. 2017, 7, 14339. [Google Scholar] [CrossRef] [PubMed]
  39. Smolinska, A.; Blanchet, L.; Coulier, L.; Ampt, K.A.M.; Luider, T.; Hintzen, R.Q.; Wijmenga, S.S.; Buydens, L.M.C. Interpretation and visualization of non-linear data fusion in kernel space: Study on metabolomic characterization of progression of multiple sclerosis. PLoS ONE 2012, 7, e38163. [Google Scholar] [CrossRef] [PubMed]
  40. Li, N.; Xiong, Z.; Deng, T.; Ren, K. Automated brain tumor segmentation from multimodal MRI data based on Tamura texture feature and an ensemble SVM classifier. Int. J. Intel. Comput. Cyber. 2019, 12, 466–480. [Google Scholar]
  41. Chen, Q.; Dai, W.; Lin, Z.; Xie, D.; Lv, M.; Lin, Z. Effects of shading on main quality components in tea (Camellia Sinensis (L) O. Kuntze) leaves based on metabolomics analysis. Sci. Agric. Sin. 2019, 52, 1066–1077. (In Chinese) [Google Scholar]
  42. Pasikanti, K.K.; Esuvaranathan, K.; Ho, P.C.; Mahendran, R.; Kamaraj, R.; Wu, Q.H.; Chiong, E.; Chan, E.C.Y. Noninvasive urinary metabonomic diagnosis of human bladder cancer. J. Proteome Res. 2010, 9, 2988–2995. [Google Scholar] [CrossRef]
  43. Brereton, R.G.; Lloyd, G.R. Partial least squares discriminant analysis: Taking the magic away. J. Chemom. 2014, 28, 213–225. [Google Scholar] [CrossRef]
  44. Chawla, N.V.; Bowyer, K.W.; Hall, L.O.; Kegelmeyer, W.P. ‘SMOTE: Synthetic minority over-sampling technique’. J. Artif. Intell. Res. 2002, 16, 321–357. [Google Scholar] [CrossRef]
  45. Broadhurst, D.I.; Kella, D.B. Statistical strategies for avoiding false discoveries in metabolomics and related experiments. Metabolomics 2006, 2, 171–196. [Google Scholar] [CrossRef]
  46. Slupsky, C.M.; Steed, H.; Wells, T.H.; Dabbs, K.; Schepansky, A.; Capstick, V.; Faught, W.; Sawyer, M.B. Urine metabolite analysis offers potential early diagnosis of ovarian and breast cancers. Clin. Cancer Res. 2010, 16, 5835–5841. [Google Scholar] [CrossRef]
  47. Thévenot, E.A.; Roux, A.; Xu, Y.; Ezan, E.; Junot, C. Analysis of the human adult urinary metabolome variations with age, body mass index, and gender by implementing a comprehensive workflow for univariate and OPLS statistical analyses. J. Proteome Res. 2015, 14, 3322–3335. [Google Scholar] [CrossRef]
  48. Gika, H.G.; Theodoridis, G.A.; Plumb, R.S.; Wilson, D. Current practice of liquid chromatography-mass spectrometry in metabolomics and metabonomics. J. Pharm. Biomed. Anal. 2014, 87, 12–25. [Google Scholar] [CrossRef]
  49. Rubert, J.; Righetti, L.; Stranska-Zachariasova, M.; Dzuman, Z.; Chrpova, J.; Dall’Asta, C.; Hajslova, J. Untargeted metabolomics based on ultra-high-performance liquid chromatography-high-resolution mass spectrometry merged with chemometrics: A new predictable tool for an early detection of mycotoxins. Food Chem. 2017, 224, 423–431. [Google Scholar] [CrossRef] [PubMed]
  50. Ortmayr, K.; Charwat, V.; Kasper, C.; Hann, S.; Koellensperger, G. Uncertainty budgeting in fold change determination and implications for non-targeted metabolomics studies in model systems. Analyst 2017, 142, 80–90. [Google Scholar] [CrossRef]
  51. Environmental Protection Agency (EPA). Method 8061A Phthalate Esters by Gas Chromatography with Electron Capture Detection (GC/ECD); Environmental Protection Agency: Washington, DC, USA, 1996. [Google Scholar]
  52. Tai, Y.; Mo, C.; Li, Y.; Wu, X.; Wang, J.; Su, Q. Concentration and distribution of tetracycline antibiotics in soils from vegetable fields of Dongguan City. China Environ. Sci. 2011, 31, 90–95. (In Chinese) [Google Scholar]
  53. Tang, X.; Lou, C.; Wang, S.; Lu, Y.; Liu, M.; Hashmi, M.Z.; Liang, X.; Li, Z.; Liao, Y.; Qin, W.; et al. Effects of long-term manure applications on the occurrence of antibiotics and antibiotic resistance genes (ARGs) in paddy soils: Evidence from four field experiments in south of China. Soil Biol. Biochem. 2015, 90, 179–187. [Google Scholar] [CrossRef]
  54. Zhu, X.; Zeng, Q.; Xie, Q.; Ding, D.; Ru, S. Residues and risk assessment of tetracycline antibiotics in vegetable-growing soils from suburban areas of northern Guangzhou. J. Agro-Environ. Sci. 2017, 36, 2257–2266. (In Chinese) [Google Scholar]
  55. Hang, L.; Zhao, Y.; Liu, C.; Yu, Y.; He, Y.; Xu, J.; Lu, Z. Determine multiple classes of veterinary antibiotics in soil: Comparing dispersive and solid-phase extraction for sample cleanup. Chromatographia 2021, 84, 833–844. [Google Scholar] [CrossRef]
  56. Lian, Y.L. Control of maize pests and diseases and measures of rational use of pesticides. World Trop. Agric. Info. 2023, 10, 33–34. (In Chinese) [Google Scholar]
  57. Want, E.J.; Wilson, I.D.; Gika, H.; Theodoridis, G.; Plumb, R.S.; Shockcor, J.; Holmes, E.; Nicholson, J.K. Global metabolic profiling procedures for urine using UPLC-MS. Nat. Protoc. 2010, 5, 1005–1018. [Google Scholar] [CrossRef]
  58. Dunn, W.B.; Broadhurst, D.; Begley, P.; Zelena, E.; Francis-McIntyre, S.; Anderson, N.; Brown, M.; Knowles, J.D.; Halsall, A.; Haselden, J.N.; et al. Procedures for large-scale metabolic profiling of serum and plasma using gas chromatography and liquid chromatography coupled to mass spectrometry. Nat. Protoc. 2011, 6, 1060–1083. [Google Scholar] [CrossRef] [PubMed]
  59. Chambers, M.C.; MacLean, B.; Burke, R.; Amodei, D.; Ruderman, D.L.; Neumann, S.; Gatto, L.; Fischer, B.; Pratt, B.; Egertson, J.; et al. A cross-platform toolkit for mass spectrometry and proteomics. Nat. Biotechnol. 2012, 30, 918–920. [Google Scholar] [CrossRef] [PubMed]
  60. Giacomoni, F.; Le Corguille, G.; Monsoor, M.; Landi, M.; Pericard, P.; Petera, M.; Duperier, C.; Tremblay-Franco, M.; Martin, J.F.; Jacob, D.; et al. Workflow4Metabolomics: A collaborative research infrastructure for computational metabolomics. Bioinformatics 2015, 31, 1493–1495. [Google Scholar] [CrossRef] [PubMed]
  61. Smith, C.A.; Want, E.J.; O’Maille, G.; Abagyan, R.; Siuzdak, G. XCMS: Processing mass spectrometry data for metabolite profiling using nonlinear peak alignment, matching, and identification. Anal. Chem. 2006, 78, 779–787. [Google Scholar] [CrossRef] [PubMed]
  62. Gika, H.G.; Theodoridis, G.A.; Wilson, I.D. Liquid chromatography and ultra-performance liquid chromatography-mass spectrometry fingerprinting of human urine: Sample stability under different handling and storage conditions for metabonomics studies. J. Chromatogr. A 2008, 1189, 314–322. [Google Scholar] [CrossRef]
  63. Winnike, J.H.; Busby, M.G.; Watkins, P.B.; O’Connell, T.M. Effects of a prolonged standardized diet on normalizing the human metabolome. Am. J. Clin. Nutr. 2009, 90, 1496–1501. [Google Scholar] [CrossRef] [PubMed]
  64. Yin, P.; Peter, A.; Franken, H.; Zhao, X.; Neukamm, S.S.; Rosenbaum, L.; Lucio, M.; Zell, A.; Haring, H.U.; Xu, G.; et al. Preanalytical aspects and sample quality assessment in metabolomics studies of human blood. Clin. Chem. 2013, 59, 833–845. [Google Scholar] [CrossRef] [PubMed]
  65. Schievano, E.; Morelato, E.; Facchin, C.; Mammi, S. Characterization of markers of botanical origin and other compounds extracted from unifloral honeys. J. Agric. Food Chem. 2013, 61, 1747–1755. [Google Scholar] [CrossRef] [PubMed]
  66. Pontes, T.A.; Barbosa, A.D.; Silva, R.D.; Melo-Junior, M.R.; Silva, R.O. Osteopenia-osteoporosis discrimination in postmenopausal women by 1H NMR-based metabonomics. PLoS ONE 2019, 14, e0217348. [Google Scholar] [CrossRef]
  67. Zhao, W.; Wang, G.; Xun, W.; Yu, Y.; Ge, C.; Liao, G. Selection of water-soluble compounds by characteristic flavor in Chahua Chicken muscles based on metabolomics. Sci. Agric. Sin. 2020, 53, 1627–1642. (In Chinese) [Google Scholar]
  68. Dai, Y.; Lü, C.; He, L.; Yi, C.; Liu, X.; Huang, W.; Chen, J. Metabolic changes in the processing of Yunkang 10 sun-dried green tea based on metabolomics. Sci. Agric. Sin. 2020, 53, 357–370. (In Chinese) [Google Scholar]
  69. Yang, J.; Chen, T.; Sun, L.; Zhao, Z.; Wan, C. Potential metabolite markers of schizophrenia. Mol. Psychiatry 2013, 18, 67–78. [Google Scholar] [CrossRef] [PubMed]
  70. Mendez, K.M.; Reinke, S.N.; Broadhurst, D.I. A comparative evaluation of the generalised predictive ability of eight machine learning algorithms across ten clinical metabolomics data sets for binary classification. Metabolomics 2019, 15, 150. [Google Scholar] [CrossRef] [PubMed]
  71. Mahadevan, S.; Shah, S.L.; Marrie, T.J.; Slupsky, C.M. Analysis of metabolomic data using support vector machines. Anal. Chem. 2008, 80, 7562–7570. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Total ion chromatograms of spiked maize sample groups on the W4M platform.
Figure 1. Total ion chromatograms of spiked maize sample groups on the W4M platform.
Molecules 29 03026 g001
Figure 2. PCA score plot of spiked maize sample groups.
Figure 2. PCA score plot of spiked maize sample groups.
Molecules 29 03026 g002
Figure 3. OPLS-DA score plots of spiked maize sample groups ((a): 20 ng/mL vs. 50 & 100 ng/mL; (b): 100 ng/mL vs. 20 & 50 ng/mL).
Figure 3. OPLS-DA score plots of spiked maize sample groups ((a): 20 ng/mL vs. 50 & 100 ng/mL; (b): 100 ng/mL vs. 20 & 50 ng/mL).
Molecules 29 03026 g003
Figure 4. S-plot plots of spiked maize sample groups ((a): 20 ng/mL vs. 50 & 100 ng/mL; (b): 100 ng/mL vs. 20 & 50 ng/mL).
Figure 4. S-plot plots of spiked maize sample groups ((a): 20 ng/mL vs. 50 & 100 ng/mL; (b): 100 ng/mL vs. 20 & 50 ng/mL).
Molecules 29 03026 g004
Figure 5. Permutation test plots of spiked maize sample groups ((a): 20 ng/mL vs. 50 & 100 ng/mL; (b): 100 ng/mL vs. 20 & 50 ng/mL).
Figure 5. Permutation test plots of spiked maize sample groups ((a): 20 ng/mL vs. 50 & 100 ng/mL; (b): 100 ng/mL vs. 20 & 50 ng/mL).
Molecules 29 03026 g005
Table 1. Marker compounds were screened in maize sample groups.
Table 1. Marker compounds were screened in maize sample groups.
No.Var ID
(Primary)
Marker
Compounds
VIP
Pred a
(Weight Squared
Value × 105) b
m/z ct
(min) d
Mass Error
(ppm) e
LOD
(µg/kg)
1M406T575Difenoconazole7.378/12.94256,490/23,814406.071309.58−1.6540.5
2M326T524Benalaxyl5.798/5.9759548/20,867326.175168.730.2791.1
3M336T541Zoxamide4.658/5.73211,861/11,744336.032839.022.6570.7
4M218T216Pymetrozine7.828/5.64925,040/13,579218.103713.600.3370.7
5M372T574Profenofos6.189/5.53126,626/4524372.942689.570.6990.9
6M353T444Tebufenozide4.262/4.7613052/4239353.223837.404.1781.0
7M365T664Pyridaben0.514/4.49273/3656365.1440811.07−2.2180.9
8M256T210Trichlorfon2.470/4.263689/3233256.929413.50−1.7530.5
9M218T300Propanil2.321/4.243599/2611218.013155.00−1.1601.3
10M220T226Dichlorvos4.161/4.2252493/2742220.953673.772.2200.9
11M304T410Fenpropimorph3.241/4.1911248/2663304.263616.830.3861.5
12M223T241Acetamiprid4.561/4.1703386/2965223.074894.021.7610.5
13M349T621Clorpyrifos3.678/4.1411708/2579349.9339310.351.0490.5
14M192T245Carbendazim2.651/4.094746/2573192.076374.08−1.9571.1
15M732T421Spinosad4.227/4.0592608/2462732.468557.020.5881.0
16M294T344Paclobutrazol2.783/4.041815/2289294.136895.730.4200.9
17M345T650 Oxadiazon3.066/3.9971010/2326345.0761310.83−1.7241.0
18M279T211Oxadixyl3.629/3.9801621/2292279.133633.52−1.0630.9
19M276T326Dimethenamid3.496/3.9681471/2338276.081425.43−1.9040.8
20M282T660Pendimethalin4.213/3.8732396/2265282.1448011.00−0.0920.7
21M305T541Diazinon4.289/3.8442528/2250305.108309.02−0.1140.4
22M302T228Flutriafol4.602/3.8133265/2173302.110463.801.6950.4
23M330T429Epoxiconazol2.132/3.739394/2123330.080777.151.1410.4
24M324T499Flutolanil5.143/3.5694610/1575324.119928.32−2.0680.3
25M292T384Cyproconazol3.236/3.4881016/1409292.121316.400.6450.5
26M253T239Hexazinone0.434/3.44644/1392253.165203.98−2.7740.8
27M319T370Pyriftalid2.171/3.416347/1241319.074436.17−0.8850.9
28M321T561Chlorpyrifos-methyl3.301/3.3491028/1237321.902819.351.7021.1
29M343T308Thiophanate-methyl1.896/3.345249/1184343.052045.13−2.5991.1
30M306T540Pirimiphos-methyl7.863/3.34324,423/1100306.103519.00−0.2170.6
31M888T502Emamectin benzoate1.008/3.257106/1067888.546038.37−0.8301.1
32M337T436Fenbuconazole1.937/3.235245/1010337.121037.27−1.2600.9
33M289T354Myclobutanil2.280/3.197347/895289.121435.90−0.0820.7
34M224T228Monocrotophos4.047/3.1451844/885224.067603.80−2.8480.8
35M222T196Carbofuran2.358/3.126362/856222.112303.27−0.7730.5
36M304T200Fenamiphos1.974/3.103237/790304.113163.330.2670.9
37M299T446Quinalphos4.042/3.0291806/729299.061837.431.5141.2
38M338T555Bitertanol2.887/2.982588/752338.186719.251.2230.9
39M300T306Phosphamidon1.897/2.956199/789300.076685.101.5360.8
40M256T277Phosfolan1.855/2.946188/728256.022834.621.0781.1
41M271T540Cadusafos3.111/2.930729/640271.094089.00−3.3561.1
42M376T359Prochloraz1.725/2.915158/619376.038415.980.8500.5
43M226T450Cyprodinil3.535/2.7861080/608226.133157.50−3.1960.5
44M243T411Mocap1.673/2.746132/591243.063676.85−0.0800.9
45M215T432Metribuzin2.828/2.720507/396215.096617.202.3120.9
46M436T523Fipronil4.356/2.6982316/545436.945898.72−0.1990.9
47M250T234Clothianidin1.491/2.63095/531250.016933.903.7080.9
48M253T261Thiacloprid1.262/2.62168/460253.030374.35−2.1580.5
49M292T214Thiamethoxam1.650/2.596116/582292.026553.57−0.0660.4
50M294T360Triadimefon3.591/2.5721107/393294.100896.001.7460.5
51M302T266Methidathion3.877/2.5381467/386302.969134.43−0.0490.9
52M368T509Anilofos1.299/2.49165/377368.030638.480.2810.9
53M299T346Phoxim1.463/2.43079/367299.061425.770.1491.0
54M318T521Phosmet0.495/2.37431/326318.001808.68−0.0500.6
55M293T309Etrimfos2.431/2.362279/281293.071305.15−2.2000.8
56M284T503Metolachlor1.320/2.28957/312284.141468.380.9930.7
57M330T488Iprodione1.610/2.27988/299330.040778.130.2930.6
58M216T346Atrazine2.591/2.276334/284216.101925.774.0240.5
59M307T499Sulfotep2.147/2.267186/321307.052318.32−1.7960.6
60M208T297Fenobucarb1.281/2.24751/261208.133204.95−0.0591.0
61M256T248Imidacloprid3.372/2.230845/260256.059434.13−0.5951.0
62M203T178Dinotefuran2.636/2.201345/241203.113812.97−0.3051.1
63M249T290Linuron2.276/2.189521/222249.019044.83−0.6890.8
64M242T356Prometryn1.801/2.173108/198242.144145.933.0790.5
65M230T309Terbutylazine2.165/2.136178/220230.117035.151.4430.8
66M329T323Pencycuron1.550/2.12471/221329.142035.381.5350.7
67M318T276Azinphos-methyl1.632/1.67560/94318.014304.603.9160.3
68M279T490Fenthion0.372/1.4388/54279.027308.17−0.0060.9
69M214T317Simetryn3.354/2.111815/190214.112645.282.5810.5
70M313T320Praziquantel1.624/2.04576/183313.191305.330.7931.0
71M479T347Chlortetracycline3.939/2.0201530/182479.122165.781.2341.1
72M463T377Tetracycline2.347/2.005221/176463.171836.281.5600.3
73M445T384Doxycycline3.592/1.9901061/167445.160716.400.3900.5
74M277T489Sulfabenzamide2.596/1.981310/162277.064278.150.4620.4
75M275T541Ormetoprim1.289/1.98041/159275.150839.022.0930.6
76M281T212Sulfamonomethoxine2.221/1.979181/147281.070143.53−0.5190.9
77M315T378Sulfaphenazole0.518/1.90828/144315.091296.300.8700.9
78M215T306Sulfacetamide1.919/1.893112/125215.049155.103.0560.9
79M281T361Sulfameter2.102/1.822142/112281.070376.020.2740.8
80M279T225Sulfamethazine2.008/1.818123/111279.091613.752.1040.9
81M254T427Sulfamethoxazole1.564/1.81558/118254.059897.121.9791.0
82M265T308Sulfamerazine1.827/1.81391/106265.075935.132.1020.8
83M254T429Sulfamethoxazole0.616/1.81116/105254.059377.15−0.0630.8
84M311T666Sulfamethazine1.671/1.80071/103311.0805511.10−0.9620.8
85M301T301Sulfaquinoxaline1.493/1.79851/102301.075895.021.7391.0
86M285T412sulfachloropyridazine1.001/1.75819/95285.021136.871.3160.8
87M215T293Sulfaguanidine1.926/1.734103/97215.059634.88−0.4030.9
88M251T288Sulfadiazine2.213/1.692164/96251.059424.80−1.1770.4
89M256T275Sulfathiazole2.800/1.687396/99256.021804.583.5310.4
90M279T377Sulfisomidine2.312/1.671192/77279.091466.281.5840.4
91M268T208Sulfafurazole1.659/1.66864/87268.075773.472.7111.0
92M291T365Trimethoprim1.707/1.65270/77291.145926.082.5750.9
93M281T380Sulfamethoxypyridazine0.735/1.63626/89281.070316.330.0660.8
94M360T414Enrofloxacin3.164/1.616637/75360.171986.900.4910.8
95M320T382Norfloxacin1.198/1.60124/75320.141946.374.5070.7
96M334T383Pefloxacin2.511/1.588260/73334.155816.38−1.0230.9
97M332T404Ciprofloxacin1.473/1.57337/70332.141046.731.6210.7
98M362T385Ofloxacin1.266/1.53527/68362.151516.421.2520.9
99M386T456Sarafloxacin1.148/1.50747/67386.132037.602.4920.4
100M352T427Lomefloxacin0.728/1.4517/62352.146037.12−1.9730.3
101M233T608Nalidixic acid1.036/1.44915/56233.0924310.131.5280.3
102M262T619Flumequine2.101/1.441128/55262.0876010.320.7700.5
103M358T411Danofloxacin0.814/1.4157/52358.157166.852.8230.8
104M400T442Difloxacin1.830/1.41477/50400.146837.370.2820.8
105M396T433Orbifloxacin1.500/1.39538/48396.153717.221.9271.0
106M393T494Sparfloxacin0.671/1.39310/48393.172688.23−1.4980.9
107M734T616Erythromycin1.106/1.36516/46734.4680810.27−0.5990.9
108M837T659Roxithromycin1.236/1.35423/38837.5316710.98−0.2141.0
109M828T630Josamycin1.349/1.35168/36828.4746110.500.7330.8
110M916T611Tylosin1.039/1.34813/36916.5268910.180.5050.9
111M702T622Kitasamycin1.758/1.34565/36702.4059310.370.0100.8
112M869T566Tilmicosin0.491/1.3437/35869.573749.430.4841.0
113M128T2542-methyl-5-nitroimidazole2.095/1.336269/35128.045514.230.4440.9
114M114T2654-nitroimidazole1.036/1.33314/32114.029894.420.8211.0
115M164T2465-nitrobenzimidazole1.122/1.31916/32164.046134.104.1151.3
116M142T207Dimetridazole1.492/1.31136/32142.061633.453.7550.5
117M172T218Metronidazole2.251/1.281164/31172.071803.630.7791.1
118M201T240Ronidazole1.356/1.28026/31201.062304.002.3150.4
119M158T329Hydroxy dimetridazole2.126/1.262131/31158.056465.482.7970.4
120M170T362Ipronidazole1.305/1.25823/31170.092776.032.1560.9
Note: a two VIP values from 20 vs. 50 and 100 and 100 vs. 20 and 50, respectively; b two weight squared values (×105) from 20 vs. 50 and 100 and 100 vs. 20 and 50, respectively; c m/z represents extracted molecular weight from W4M platform; d t represents retention time from W4M platform; e Mass error (ppm) = (extracted molecular weight from W4M platform − extracted molecular weight from LC–MS/MS) × 106/extracted molecular weight from LC–MS/MS.
Table 2. Basic information on 69 pesticides (No. 1~69) and 55 veterinary drugs (No. 70~124).
Table 2. Basic information on 69 pesticides (No. 1~69) and 55 veterinary drugs (No. 70~124).
No.CompoundsCategoryCAS No.Extracted Molecular WeightRetention Time (min)
1Dinotefuraninsecticide165252-70-0203.113873.06
2Carbofuraninsecticide1563-66-2222.112473.21
3Fenamiphosinsecticide22224-92-6304.113083.44
4Trichlorfoninsecticide52-68-6256.929863.47
5Thiamethoxaminsecticide153719-23-4292.026573.54
6Pymetrozineinsecticide123312-89-0218.103643.61
7Dichlorvosinsecticide62-73-7220.953183.78
8Monocrotophosinsecticide2157-98-4224.068243.81
9Clothianidininsecticide210880-92-5250.016003.95
10Acetamipridinsecticide135410-20-7223.074504.03
11Imidaclopridinsecticide105827-78-9256.059584.09
12Thiaclopridinsecticide111988-49-9253.030924.33
13Methidathioninsecticide950-37-8302.969144.37
14Phosfolaninsecticide947-02-4256.022554.54
15Azinphos-methylinsecticide86-50-0318.013054.59
16Fenobucarbinsecticide3766-81-2208.133214.96
17Phosphamidoninsecticide13171-21-6300.076225.05
18Etrimfosinsecticide38260-54-7293.071945.16
19Phoximinsecticide14816-18-3299.061385.71
20Mocapinsecticide13194-48-4243.063696.89
21Spinosadinsecticide131929-60-7732.468127.15
22Tebufenozideinsecticide112410-23-8353.222357.35
23Quinalphosinsecticide13593-03-8299.061387.41
24Fenthioninsecticide55-38-9279.027308.07
25Sulfotepinsecticide3689-24-5307.052868.26
26Emamectin benzoateinsecticide155569-91-8888.546778.54
27Phosmetinsecticide732-11-6318.001828.57
28Fipronilinsecticide120068-37-3436.945988.73
29Cadusafosinsecticide95465-99-9271.094998.95
30Diazinoninsecticide333-41-5305.108338.96
31Pirimiphos-methylinsecticide29232-93-7306.103588.96
32Chlorpyrifos-methylinsecticide5598-13-0321.902269.32
33Profenofosinsecticide41198-08-7372.942429.49
34Clorpyrifosinsecticide2921-88-2349.9335610.26
35Pyridabeninsecticide96489-71-3365.1448911.06
36Oxadixylbactericide77732-09-3279.133933.49
37Flutriafolbactericide76674-21-0302.109953.79
38Carbendazimbactericide10605-21-7192.076754.11
39Thiophanate-methylbactericide23564-05-8343.052935.14
40Pencycuronbactericide66063-05-6329.141525.47
41Prochlorazbactericide67747-09-5376.038095.89
42Myclobutanilbactericide88671-89-0289.121456.06
43Triadimefonbactericide43121-43-3294.100386.11
44Cyproconazolbactericide113096-99-4292.121126.36
45Fenpropimorphbactericide67306-03-0304.263496.72
46Epoxiconazolbactericide106325-08-0330.080397.11
47Fenbuconazolebactericide114369-43-6337.121457.33
48Cyprodinilbactericide121552-61-2226.133877.64
49Iprodionebactericide36734-19-7330.040678.13
50Flutolanilbactericide66332-96-5324.120598.23
51Benalaxylbactericide71626-11-4326.175078.75
52Zoxamidebactericide156052-68-5336.031948.99
53Bitertanolbactericide55179-31-2338.186309.17
54Difenoconazolebactericide119446-68-3406.071979.64
55Hexazinoneherbicide51235-04-2253.165903.94
56Linuronherbicide330-55-2249.019214.92
57Propanilherbicide709-98-8218.013404.98
58Terbutylazineherbicide5915-41-3230.116705.07
59Simetrynherbicide1014-70-6214.112095.32
60Dimethenamidherbicide87674-68-8276.081955.41
61Atrazineherbicide102029-43-6216.101055.77
62Prometrynherbicide7287-19-6242.143395.97
63Pyriftalidherbicide135186-78-6319.074716.02
64Metribuzinherbicide21087-64-9215.096117.14
65Metolachlorherbicide51218-45-2284.141188.34
66Anilofosherbicide64249-01-0368.030538.49
67 Oxadiazonherbicide19666-30-9345.0767210.77
68Pendimethalinherbicide40487-42-1282.1448310.92
69Paclobutrazolgrowth regulator76738-62-0294.136775.77
70Sulfafurazolesulfamido127-69-5268.075043.38
71Sulfamonomethoxinesulfamido1220-83-3281.070293.54
72Sulfamethazinesulfamido57-68-1279.091023.76
73Sulfathiazolesulfamido72-14-0256.020904.54
74Sulfadiazinesulfamido68-35-9251.059724.73
75Sulfaguanidinesulfamido57-67-0215.059724.85
76Sulfaquinoxalinesulfamido59-40-5301.075375.05
77Sulfacetamidesulfamido144-80-9215.048495.06
78Sulfamerazinesulfamido127-79-7265.075375.11
79Sulfametersulfamido651-06-9281.070295.95
80Trimethoprimsulfamido738-70-5291.145176.01
81Sulfaphenazolesulfamido526-08-9315.091026.21
82Sulfisomidinesulfamido515-64-0279.091026.24
83Sulfamethoxypyridazinesulfamido80-35-3281.070296.51
84Sulfachloropyridazinesulfamido80-32-0285.020756.87
85Sulfamethoxazolesulfamido723-46-6254.059397.09
86Sulfamethoxazolesulfamido723-46-6254.059397.09
87Sulfabenzamidesulfamido127-71-9277.064148.05
88Ormetoprimsulfamido6981-18-6275.150258.94
89Sulfamethazinesulfamido57-68-1311.0808511.05
90Dimetridazolenitroimidazoles551-92-8142.061103.46
91Metronidazolenitroimidazoles443-48-1172.071673.58
92Ronidazolenitroimidazoles7681-76-7201.061833.96
935-nitrobenzimidazolenitroimidazoles94-52-0164.045454.05
942-methyl-5-nitroimidazolenitroimidazoles88054-22-2128.045454.17
954-nitroimidazolenitroimidazoles3034-38-6114.029804.41
96Hydroxy dimetridazolenitroimidazoles936-05-0158.056025.55
97Ipronidazolenitroimidazoles14885-29-1170.092406.05
98Mebendazolenitroimidazoles31431-39-7296.102977.74
99Fleroxacinquinolones79660-72-3370.137305.99
100Ofloxacinquinolones82419-36-1362.151066.39
101Pefloxacinquinolones70458-92-3334.156156.41
102Norfloxacinquinolones70458-96-7320.140506.51
103Ciprofloxacinquinolones85721-33-1332.140506.73
104Enrofloxacinquinolones93106-60-6360.171807.01
105Danofloxacinquinolones112398-08-0358.156157.02
106Lomefloxacinquinolones98079-51-7352.146727.12
107Orbifloxacinquinolones113617-63-3396.152957.27
108Difloxacinquinolones98106-17-3400.146727.42
109Sarafloxacinquinolones98105-99-8386.131077.64
110Sparfloxacinquinolones110871-86-8393.173278.34
111Nalidixic acidquinolones389-08-2233.0920710.01
112Flumequinequinolones42835-25-6262.0874010.15
113Praziquantelvermifuge55268-74-1313.191055.31
114Chlortetracyclinetetracyclines57-62-5479.121575.68
115Tetracyclinetetracyclines60-54-8463.171116.21
116Doxycyclinetetracyclines564-25-0445.160546.34
117Lincomycinmacrolides154-21-2407.221049.41
118Tilmicosinmacrolides108050-54-0869.573329.49
119Clindamycinmacrolides18323-44-9425.187159.68
120Tylosinmacrolides1401-69-0916.5264310.23
121Erythromycinmacrolides114-07-8734.4685210.25
122Kitasamycinmacrolides1392-21-8702.4059210.48
123Josamycinmacrolides16846-24-5828.4740010.65
124Roxithromycinmacrolides80214-83-1837.5318510.84
125Atrazine-d5 163165-75-1221.140175.77
126Enrofloxacin-d5 1173021-92-5365.210927.01
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xue, W.; Li, F.; Li, X.; Liu, Y. A Support Vector Machine-Assisted Metabolomics Approach for Non-Targeted Screening of Multi-Class Pesticides and Veterinary Drugs in Maize. Molecules 2024, 29, 3026. https://doi.org/10.3390/molecules29133026

AMA Style

Xue W, Li F, Li X, Liu Y. A Support Vector Machine-Assisted Metabolomics Approach for Non-Targeted Screening of Multi-Class Pesticides and Veterinary Drugs in Maize. Molecules. 2024; 29(13):3026. https://doi.org/10.3390/molecules29133026

Chicago/Turabian Style

Xue, Weifeng, Fang Li, Xuemei Li, and Ying Liu. 2024. "A Support Vector Machine-Assisted Metabolomics Approach for Non-Targeted Screening of Multi-Class Pesticides and Veterinary Drugs in Maize" Molecules 29, no. 13: 3026. https://doi.org/10.3390/molecules29133026

APA Style

Xue, W., Li, F., Li, X., & Liu, Y. (2024). A Support Vector Machine-Assisted Metabolomics Approach for Non-Targeted Screening of Multi-Class Pesticides and Veterinary Drugs in Maize. Molecules, 29(13), 3026. https://doi.org/10.3390/molecules29133026

Article Metrics

Back to TopTop