Introduction

Ochratoxin A (OTA) is a secondary fungal metabolite produced mainly by Aspergillus and Penicillium species1,2. OTA exhibits potent nephrotoxicity, hepatotoxicity, cytotoxicity, and immunotoxicity and is a potential carcinogen3,4,5. OTA has been classified as a possible human carcinogen (group 2B) by the International Agency for Research on Cancer6. OTA may reach the human food chain in contaminated, mostly cereal-based feed ingested by animals that are fated to become food for human consumption7,8. Brera and colleagues9, Endou and colleagues10, and Krogh11 found that OTA contamination in foods correlates with endemic nephropathy in humans and pigs. OTA is found predominantly in grain, animal feed, fruit12, and processed foods of varying origin, with a wide range and low concentrations. The European Union has set a maximum concentration for OTA in plant-derived products, including cereal, coffee, and grape-based foods, as between 0.5 to 10 µg kg−1 13, but for animal-derived food no recommendation concerning a maximum level has been made. To better understand the extent to which agricultural products are contaminated with OTA and allow for risk assessment, a highly sensitive and reliable analytical method is needed to determine OTA levels in agricultural products, especially animal-derived products.

Extraction of target analyte from matrix and subsequent quantification are the two key steps for analysis of a contaminating compound. Most reports concerning OTA detection have focused mainly on contamination of cereals, coffee14, fruits15,16,17,18, and pork products19,20,21. These methods mainly have used a pretreatment step that relies on immunoaffinity column and analysis by HPLC-FLD. However, because the methods involves blending, quantitative transfers, wash the immunoaffinity column, elute the target compound, and evaporate and exchange the solvent, which are procedures that are complicated, time-consuming, and require large amounts of solvents. Also, each step is a potential source of systematic and random errors. In addition, since the average levels of OTA are lower in agricultural products, the detectability and sensitivity of analytical methods used should be very high, in the reported method, the limit of quantitation (LOQ) for all sample type were between 0.2 to 1 µg kg−1 14,15,20,21, it may higher than the OTA levels in real samples. Hence, development of a simple and sensitive pretreatment would be very useful.

QuEChERS (Quick, Easy, Cheap, Effective, Rugged, and Safe) methods usually use acetonitrile (MeCN) partitioning and a dispersive solid-phase extraction (d-SPE) clean-up22. Because of their simplicity, high efficiency, and cost-effectiveness, QuEChERS procedures have been applied to the analysis of mycotoxins in barley23, rice24, and animal feed25. Furthermore, QuEChERS methods are flexible, as the optimization of their pretreatment processes allows for the extraction of target compounds from complex matrices including foods of animal origin26. In addition, an ultra-performance liquid chromatography (UPLC) system coupled to a Xevo TQ-S triple quadrupole mass spectrometer serves as a powerful tool for the determination of trace levels of pesticides27 and mycotoxins28 from various matrices. The Xevo TQ-S29 is equipped with patented StepWave ion-guide technology, which improves ion sampling in the source and ion-transfer efficiency, resulting in excellent sensitivity. Therefore, we posited that the combined use of a QuEChERS pretreatment with a UPLC-Xevo TQ-S system should allow for quantification of trace levels of OTA in foods of plant and animal origins.

The aim of this study was thus to establish a rapid and sensitive analytical method to determine the extent to which OTA residue contaminates grapes, grape-related matrices (grape wine, grape juice, and raisins), and foods of animal origin (pork, porcine kidney, liver, and fat; chicken and eggs; cow milk). To maximize purification and recovery of OTA in the matrices, the ability of MeCN containing varying amounts of acetic acid (HOAc) to quantitatively extract OTA was assessed. To our knowledge, this is the first time a QuEChERS method and UPLC-tandem mass spectroscopy have been employed to quantify OTA levels in foods of animal origin. Furthermore, we collected 717 food samples from local markets, fields, and wholesale supply warehouses and analyzed them for their OTA content to make risk assessments concerning OTA contamination in foods available for human consumption.

Results

Linearity and LOQ

The optimized method was used to quantify OTA spiked into the various matrices. Calibration curves at concentration between 0.1 µg L−1 and 500 µg L−1, each using the results for six solutions spiked with the standard solution of OTA, were prepared for each matrix. The linear regression equations, and R2 values for the matrix-matched curves are given in Table 1. For each curve, an R2 value ≥ 0.9994 was obtained, and the LOQ was 0.1 µg kg−1.

Table 1 Calibration equations, R2 values, LOQ values, and matrix effects for Ochratoxin A detection by our optimized UPLC-MS/MS method.

Matrix effects

Matrix effects are a major challenge when attempting to quantify a contaminant by LC-MS/MS because co-eluting sample constituents may enhance or suppress the analyte response, which is also dependent on the instrument, the properties of the analyte and its concentration, the matrix, and the sample pretreatment procedure30. We assessed the effect of each matrix on MS/MS detection of OTA, the ratio of the slope of each of the curves for spiked OTA samples from each matrix to that of the standard solvent curve and the values of each matrix effect are shown in Table 1. Enhancement or suppression (from −6 to 83%) of the OTA signal was observed for the 11 matrices. An OTA signal enhancement effect (26% ≤matrix effect ≤83%) was observed for all matrices except for porcine liver, for which a signal suppression effect was evident (matrix effect = −6%). Therefore, matrix-matched calibration standards were necessary to accurately quantify OTA extracted from the different matrices.

Precision and accuracy

The precision and accuracy of the method were evaluated according to the recoveries obtained for blank samples analyzed in quintuplicate and spiked with OTA at 0.1, 1, 10, or and 500 µg kg−1. Mean recovery ranged from 70.3 to 114.7% with RSD values <19.2% for all samples (Table 2); the RSDr values ranged from 0.8 to 19.2%, and the RSDR values ranged from 2.0 to 16.9%. These results indicate that the method achieved satisfactory recovery, precision, and sensitivity values for the detection of OTA in all tested matrices.

Table 2 Accuracy and precision of the optimized method for the 11 matrices, each of which was spiked with one of four concentrations of Ochratoxin A.

Occurrence of OTA in real samples

The ability of the method to measure trace levels of OTA was evaluated by analyzing 717 samples obtained from fields, wholesale warehouses, and stores. OTA was found at concentrations of 0.35, 0.61, 1.11 µg kg−1 in three grape samples acquired from a field in Xinjiang province. For three raisin samples from a market in Tianjin, the concentration of OTA was 10.14 µg kg−1, which exceeds the maximum allowed concentration for OTA in raisins as established by the EU (10 μg kg−1)13, whereas the other two concentrations were 0.18 and 0.43 μg kg−1. For 6 of the 41 grape juice samples, the detected OTA concentration was between 0.1 and 0.2 μg kg−1, in other 35 grape juice samples, OTA could be detected but the concentration was below the LOQ, whereas for the remaining samples, OTA could not be detected. Information concerning the occurrence of OTA in real samples is shown in Table 3. The high occurrence of OTA in the grape juice and unacceptably high contamination level of one raisin sample indicate that supervision of the farm-to-store chain for grape juice and raisins is necessary, as is quantification of OTA in those foods.

Table 3 Occurrence of Ochratoxin A in samples obtained in the field.

Discussion

To select the multiple reaction-monitoring parameters, the standard solution of OTA (100 µg L−1) was directly injected into the UPLC-ESI-MS mass spectrum at different cone voltage. Excellent sensitivity was achieved in the negative ionization mode, the deprotonated molecule at m/z = 402.11 was selected as precursor MS ion, the final cone voltage was 34 V. The precursor ion was fragmented at different collision energies, and the selected reaction-monitoring parameters were optimized to obtain the highest degree of sensitivity. The two most intense fragment ions were observed in the product-ion spectra—one at m/z = 167.01 and the second at m/z = 211.02; these two ions were used for quantitative and qualitative analysis of the OTA content in the samples when the collision energy was 36 V or 26 V, respectively.

To optimize the peak shapes and retention behavior, different mobile-phase compositions MeCN and water, MeCN and 0.2% (v/v) formic acid(aq), methanol and water, and methanol and 0.2% (v/v) formic acid(aq) were assessed. The methanol/0.2% formic acid(aq) produced the best combination of satisfactory retention time (2.57 min) and peak shape for OTA with no interfering peaks (Fig. 1).

Figure 1
figure 1

Typical UPLC-MS/MS multiple reaction-monitoring chromatograms for Ochratoxin A (a) in a standard solution (1 μg L−1), (b) from a grape sample collected at a field in Xinjiang province, (c) from unspiked pork and grape samples, and (d) from pork and grape samples each spiked at 1 μg L−1.

The selection of an extraction solvent and sorbent is critical for the detection of a target compound in different matrices. We first considered MeCN as the extraction solvent because its use usually results in better recoveries with less co-extraction of matrix components31. D-SPE is a clean-up technique that involves mixing the sample extract with a mixture of anhydrous MgSO4 and sorbents to remove matrix components. Common sorbents include primary secondary amine (PSA), octadecylsilane (C18), graphitized carbon black (GCB), and Florisil, which remove different types of impurities based on chemical properties. PSA functions as a weak anion exchanger and has a strong sorption capacity for free fatty acids and other polar pigments, whereas C18 is a common reverse-phase sorbent and is suitable for extraction of nonpolar and partially polar compounds from polar samples. GCB is mainly used to remove pigments, and Florisil is used to separate and remove matrix secretions and grease. We assessed the ability of the following sorbents to remove contaminants from the grape and processed grape samples: 50 mg PSA; 50 mg PSA + 10 mg GCB; 50 mg C18; and 50 mg C18 + 10 mg GCB. We also assessed the following sorbents to remove contaminants from foods of animal origin: 50 mg PSA; 50 mg C18; and 50 mg Florisil. Pure MeCN resulted in satisfactory OTA recovery for the grape, grape wine, grape juice, and raisin samples when C18 was the sorbent (recoveries were between 84.56% and 112.05%; Fig. 2). However, extraction of OTA from foods of animal origin was not satisfactory when pure MeCN served as the extraction solvent (recoveries were <3.72%; Fig. 3). Then we assessed the ability of ethyl acetate to extract OTA from pork meat and porcine kidney, OTA recovery from porcine kidney was satisfactory when C18 was the sorbent (79.30%), but ethyl acetate could barely extract OTA from pork meat. Because organic solvents mixed with acidic solutions are often used to extract OTA from solid samples32, we assessed the ability of MeCN containing 4 or 5% (v/v) HOAc to extract OTA from pork meat and porcine kidney (Fig. 3). Recoveries increased as the concentration of HOAc increased when C18 was the sorbent. OTA recovery from pork was not satisfactory with use of MeCN containing 4% HOAc (recovery 55.03–69.04%); when HOAc was increased to 5%, however, mean recovery increased substantially (93.98% and 92.01% for the pork and kidney samples, respectively). The efficiency of this combination of extraction solution and sorbent was acceptable for all foods of animal origin. Given these results, we used MeCN as the extraction solvent for grapes and processed grape products, and MeCN containing 5% HOAc (v/v) was used for all foods of animal origin. As shown in Figs 2 and 3, Recovery and RSD values were always acceptable when 50 mg of C18 was used to purify OTA from a matrix (Figs 2 and 3). Therefore, C18 was used as the sorbent for all additional experiments.

Figure 2
figure 2

Effect of different types of sorbents on the detection of spiked Ochratoxin A (10 µg kg−1) in different processed grape products (n = 3 per matrix).

Figure 3
figure 3

Effect of different extraction solvents and sorbents on the detection of spiked Ochratoxin A (10 µg kg−1) in different foods of animal origin (n = 3 per matrix).

Conclusions

In the present study, we established and validated a simple and robust analytical method using UPLC-MS/MS in the negative ionization mode to detect OTA in grapes, processed grape products, and foods of animal origin. OTA eluted within 5 min with satisfactory recovery, linearity, accuracy, precision and a low LOQ. The detection of OTA in real samples using the method confirmed its reliability and efficacy for the routine monitoring of OTA in grapes, processed grape products, and foods of animal origin. Identification of OTA in samples from the food chain revealed that OTA may be present at high levels in raisins, mostly, at low levels in grapes and grape juice, although in one raisin sample the amount of OTA was above the maximum allowable limit, suggesting that supervision during grape product processing will be needed.

Materials and Methods

Chemicals and reagents

OTA (99.76%, analytical standard grade) was obtained from Fermentek Ltd. HPLC-grade methanol was purchased from Sigma-Aldrich (Steinheim, Germany). Anhydrous MgSO4 (analytical grade) NaCl, MeCN, and HOAc were purchased from Beihua Fine-Chemicals Co. (Beijing, China). Ultrapure water was prepared using a Milli-Q system (Bedford, MA, USA). PSA (40 µm), C18 (40 µm), Florisil, and GCB (40 µm) sorbents and 0.22-µm pore nylon syringe filters were provided by Agela Technologies Inc. (Tianjin, China).

A standard stock solution (100 mg L−1) of OTA was prepared in MeCN. Standard working OTA solutions at concentrations of 0.1, 1, 5, 10, 100 and 500 µg L−1 were prepared by diluting samples of the standard stock solution with MeCN. Accordingly, matrix-matched standard solutions of 0.1, 1, 5, 10, 100 and 500 µg L−1 OTA were obtained by adding an appropriate volume of a sample extract (grapes, grape wine, grape juice, raisins, pork, porcine kidney, liver, or fat, and chicken, eggs, or milk) to each serially diluted standard solution. All solutions were protected against light and stored in a refrigerator at −20 °C until used.

UPLC-MS/MS analysis

A Waters ACQUITY UPLC system (Milford, MA, USA) consisting of a UPLC binary solvent manager, a UPLC manager, and a column heater equipped with a UPLC HSS T3 column (2.1 × 100 mm, 1.8-μm particle size) was used. The mobile phase consisted of methanol (solvent A) and 0.2% (v/v) formic acid(aq) (solvent B). Elution was performed in the gradient mode (0 min, 10% A; 1.5 min, 90% A; 3.0 min, 90% A; 3.1 min, 10% A; 5.0 min, 10% A) for a total of 5 min. The flow rate was 300 µL/min, and the injection volume was 5 µL. Temperatures of the sample manager and column were 5 and 40 °C, respectively.

Qualitative and quantitative analyses of OTA were performed using a Waters XEVO TQ-S tandem quadrupole mass spectrometer equipped with an electrospray ionization source. Nitrogen (99.95%) and argon (99.999%) each at a pressure of 2 × 10–3 mbar served as the nebulizer and collision gas in the T-wave cell, respectively. Tandem mass spectral detection was performed by multiple reaction monitoring in the negative ionization mode. The monitoring conditions were optimized for OTA. All parameters for the multiple reaction-monitoring transitions, cone voltage, and collision energy were optimized to obtain the best sensitivity and resolution. Typical parameters were: capillary voltage, 3.1 kV; source temperature, 150 °C; desolvation temperature, 200 °C. Nitrogen was the cone and desolvation gas supplied at 150 and 800 L h−1, respectively. MassLynx software (version 4.1) was used to collect and analyze the data.

Field sample collection and preparation

A total of 717 samples were collected from the main production areas and/or the major provinces of China in 2016 (Table 4). The grape samples were collected from a field located in each county of the main production areas at harvest. Milk samples were collected from cows that had been pastured in the main production areas. Samples of processed grape products and foods of animal origin were collected from wholesaler supply warehouses or at large supermarkets in the major provinces and municipalities in China. For each collection, product was sampled throughout the location. Sample sizes were a minimum of 1 kg. Samples were sent to the laboratory immediately after acquisition where they were chopped, homogenized, and stored at −20 °C until used.

Table 4 Geographical locations and numbers of the samples acquired in different provinces of China.

Each thoroughly homogenized 10-g matrix sample was added into a 50-mL Teflon centrifuge tube. Next, 10 mL of MeCN was added into the tubes containing grape, grape wine, grape juice, or raisins, and 10 mL MeCN containing 5% (v/v) HOAc(aq) was added into the tubes containing pork, porcine kidney, liver, or fat, chicken, eggs, or milk. The tubes were shaken vigorously for 10 min. Next, 4 g of anhydrous MgSO4 and 1 g of NaCl were added, and the tubes were vortexed again for 5 min. Samples were then centrifuged at 2811 × g for 5 min. Next, the upper 1.5-mL volume of each sample was transferred into a 2.0-mL centrifuge tube containing 50 mg of C18 sorbent and 150 mg of anhydrous MgSO4. Then, the tubes were vortexed for 1 min and centrifuged for 5 min at 2400 × g. Each upper layer was filtered through a new 0.22-µm nylon syringe filter and then transferred into a brown sampler vial for UPLC-MS/MS.

Method validation

The method was validated in terms of its sensitivity, accuracy, precision, specificity, and linearity. The limit of quantification (LOQ), a measure of a method’s sensitivity, was established based on the lowest spiked concentration in each matrix that had a recovery value between 70 and 120% and a relative standard deviation (RSD) ≤20%33.

The assays used to determine recovery in terms of precision and accuracy were performed using spiked OTA samples at a concentration of 0. 1, 1, 10, or 500 µg kg−1. Five replicates of each spiked sample at a given OTA concentration were prepared by adding the appropriate volume of the standard solution to a matrix sample, allowing the samples to settle for 30 min, and then extracting and purifying the OTA in the samples according to the associated sample-preparation procedure. Accuracy is expressed as the percent recovery of the spiked samples and repeatability (precision) as the RSD. The precision was assessed as the intra-day repeatability (RSDr) and inter-day reproducibility (RSDR). All experiments were carried out in quintuplicate on three different days. The intra-day precision was measured by comparing the standard deviation of the OTA recovery percentages of spiked samples acquired on the same day. The inter-day precision was determined by comparing the standard deviation of the OTA recovery percentages of the spiked samples on three separate days.

Blank (unspiked) samples were extracted from the matrices and analyzed to determine the specificity of the method and to identify interfering peaks that co-eluted with OTA. Identification of co-extracted compounds was assessed by monitoring the typical ion-chromatograms at the retention time intervals expected for OTA.

The linearity of the method was determined by quantifying the MS/MS peak area results for the standard solutions and the matrix-matched standard solutions in triplicate at six concentrations: 0.1, 1, 5, 10, 100, 500 µg kg−1. The linear regression equation for the standard solution was plotted, and its slope, intercept, and the correlation coefficient (R2) calculated.

The stability of OTA in the stock solvent solutions and in the matrices containing spiked samples (10 µg kg−1) was assessed monthly. All aforementioned tested samples were stored at −20 °C.

When using electrospray ionization, the presence of coextractives may have positive or negative effects on the chromatographic response of the analyte. And the matrix effect could be calculated as follows: matrix effect (ME%) = (slope of calibration curves in matrix − slope of calibration curves in solvent)/slope of calibration curves in solvent (×100%).