Usage of the relationship between the application rates of the active ingredient of fungicides and their residue levels in mature apples to creating a coherent system of MRLs

Setting of maximum residue levels (MRLs; tolerances in the USA) in crops requires a big amount of data concerning residues from a number of supervised field trials for each pesticide/crop combination. This task is time-consuming, costly and fairly complicated. Therefore, we initiated a study on the utilization of the interpolation method in supporting and facilitating of assessing the value of the MRLs which are now being in use, and in creating their consistent system in the future. A mathematical formula for predicting the initial pesticide residue level in mature apples was developed by establishing the relationship between application rates and residue levels of the pesticide active ingredients in mature apples. This dependence was described by a linear equation R0 = 1.2593 × D, with coefficient of determination r2 = 0.984, where D is an application rate of a given substance. This relationship makes it possible to predict a residue level of other substances of fungicidal activity. Thus, residue levels (both predicted by the formula and obtained in the frame of Polish National Monitoring Program) of fungicides now used in apple orchards, in mature apples were typically lower than those of the statutory MRLs. Bearing in mind that the European Union policy is directed towards substantial reduction in pesticide use, the MRLs might be reconsidered in order to be set at lower values.


Introduction
Among the food safety hazards for the consumers, pesticide residues, some bacterial pathogens (e.g. Salmonella enterica) and foodborne viruses (e.g. norovirus) have been identified as three most important risk factors [4]. In order to safeguard consumers' health and to promote principles of good agricultural practice (GAP) in the use of pesticides, maximum residue levels (MRLs; tolerances in the USA) of pesticides have been set by all EU Member States applying the same evaluation procedures and authorization criteria in order to put a plant protection product (PPP) on the market. These MRLs represent the highest concentrations of pesticides (expressed in mg of active ingredient per kg) which are legally permitted in food commodities, and they are interpreted as the highest residue of a given substance which may be found if a pesticide, registered and authorized, is applied according to label and, therefore, also to principles of GAP [8]. So defined MRLs can be established only when the residues in a crop resulting from particular use pattern of the pesticide meet the public health risk assessment criteria [8,14,19]. Despite international discrepancies between food safety regulations, compliance with MRLs is still an essential prerequisite in trade of food and agricultural products [6,23,24].
The residue levels of pesticides occurring in crops at harvest time are dependent on many factors and are susceptible to influence of parameters such as a spray equipment used, a spray quality and conditions, crop management, crop varieties, practices of growing, growth stages at time of application and weather conditions [8]. But the most important determining factors are application rate of the active ingredient (AI) per hectare and an interval between the last application and the harvest (pre-harvest interval PHI). It appears fairly obvious that relatively high residues may be expected when pesticides are applied at high doses later in growing season when plant growth rate is slow [17].
For setting MRLs, residue data obtained from a number of supervised field trials are required for each pesticide/ crop combination. This task is time-consuming, costly and fairly complicated to perform [2,3,8,16]. However, development of a consistent approach for supporting and facilitating estimation of the MRLs based on effective application rate (dose) of the AI appears to be feasible. For this reason, we initiated the field trials to investigate whether the interpolation method may be efficient in predicting of the residue levels of fungicides now used in apple orchards against diseases that potentially may develop during fruit storage. Apples are particularly interesting study object, since these fruits are heavily exposed to pesticide contamination due to numerous pesticide treatments to whom they are subjected, and because they are the most highly consumed fruits (along with oranges) in the European Union countries and the USA.

Field trials
The field trials were carried out in a commercial orchard which was located near the Kraśnik town (Lublin Province, south-eastern Poland) in 2011-2012. The orchard is specialized in providing apples for baby food production and was certified in compliance with Integrated Pest Management (IPM) system by the Main Inspectorate of Plant Health and Seed Inspection (PIORiN-Państwowa Inspekcja Ochrony Roślin i Nasiennictwa).
Trial 1: To protect apples against fungal diseases that develop during the cold storage, apple trees of Gloster variety were sprayed with Switch 62.5 WG, a commercial product containing 375 g kg -1 of cyprodinil and 250 g kg -1 of fludioxonil in the form of water dispersible granules (WG), at a dose of 0.8 kg ha -1 . Switch 62.5 WG, which is a mixture of compounds belonging to anilinopyrimidine and phenylpyrrole chemical groups, exhibits a protective and systemic activity against infections of fruits and vegetables by Botrytis cinerea (grey mould).
Trial 2: One week before harvesting mature apples, and submitting them to the storage room, apple trees of Lobo variety were sprayed with Zato 50 WG, at a dose of 0.2 kg ha -1 . Trifloxystrobin, the AI of this plant protection product, belongs to the family of strobilurin fungicides and exhibits systemic mode of action. Trials 3 and 4: These two field trials were carried out on the apple trees of Gloster and Lobo varieties, which were sprayed with fungicides Merpan 80 WG and Captan 80 WG (in both cases, the active ingredient was captan) at a dose of 1.9 kg ha -1 . According to the label, Merpan 80 WG or 50 WP is recommended to control apple scab (Venturia inaequalis). Preparations containing captan (and folpet belonging to the same phthalimide group) are the fungicides commonly used to control the mentioned disease during the whole vegetative period.
Trial 5: Apple trees of Golden Delicious variety were sprayed with the Bellis 38 WG (25.5 % boscalid and 12 % pyraclostrobin) using 0.8 kg of the fungicide per ha. Bellis 38 WG is composed of two AIs: pyraclostrobin belonging to strobilurin group, and boscalid, preventive fungicide with curative action belonging to the anilide chemical group.
Trial 6: Similarly to the field trial 2, apple trees of Gloster variety were sprayed with Zato 50 WG, but at a lower dose of 0.15 kg ha -1 .

Sampling
Each time, sampling was made at the next day after pesticide application (approximately after 12 h). Each laboratory sample consisted of eight apples collected from randomly selected trees. Of each field trial, four laboratory samples were taken, each of them from one randomly selected row of the apple trees. Fruit samples were transported to the laboratory where they were homogenized, and then analytical portions of known weight (each of them weighing about a hundred grams) were used for the pesticide residue analysis.

Pesticide residue analysis
Reagents and standards: acetone and dichloromethane of HPLC grade were obtained from Sigma-Aldrich Sp. z o.o., Poznań, Poland. Petroleum ether (Chempur, Piekary Ś ląskie, Poland) and diethyl ether (POCH, Gliwice, Poland) were distilled from glass before use. Sodium sulphate(VI) (POCH, Gliwice, Poland) was heated at 550°C for 7 h. Florisil (Sigma-Aldrich Sp. z o.o., Poznań, Poland) was activated by heating at 130-135°C for 7 h and stored in a desiccator before use. Certified pesticide analytical standards were obtained from Ehrenstorfer (Augsburg, Germany). Stock solutions of approximately 1000 lg mL -1 were prepared in acetone and stored at -17°C of which intermediate concentration standards (10 lg mL -1 ) were prepared by dilution with acetone (stored at 4°C). Working standard mixtures were obtained by diluting the intermediate concentration solutions with appropriate volumes of petroleum ether and stored at 4°C.
Sample preparation: the sample preparation was carried out according to the analytical method which was previously described in detail [18,21]. Briefly, an analytical portion (about 100 g) was homogenized for 2 min with 150 mL of acetone. The extract was filtered, and then an aliquot of the filtrate (equivalent to 20 g of the analytical portion) was taken for further processing. After addition of 100 mL of 2.5 % aqueous solution of sodium sulphate(VI), it was partitioned sequentially liquid-liquid with 20, 10 and 10 mL dichloromethane. The combined extracts were evaporated to dryness and dissolved with 10 mL petroleum ether. An aliquot of 5 mL was transferred into a glass column which was previously packed with 1.1 g activated florisil and 4.5 g anhydrous sodium sulphate(VI). The residues were eluted with 70 mL diethyl ether/petroleum ether (3:7, v/v) and 70 mL acetone/petroleum ether (3:7, v/ v). The combined extracts were evaporated and diluted with petroleum ether. The concentration of matrix in the final extract was 1 g per mL.
Gas chromatographic determination: a gas chromatograph model 7890 (Agilent Technologies, Palo Alto, CA, USA) was used. The instrument was equipped with a micro-electron capture (lEC) and nitrogen-phosphorus (NP) detectors. The HP-5 MS Ultra Inert capillary column 30 m 9 0.32 mm 9 0.25 lm film thickness (Agilent Technologies, Palo Alto, CA, USA) was used. Sample extracts of 2 ll were injected splitless at 250°C. The column temperature was ramped from 100 to 180°C at 10°C min -1 (held for 4 min), then ramped to 220°C at 3°C min -1 (held for 15 min), then ramped to 260°C at 10°C min -1 (held for 11 min). Cyprodinil and fludioxonil were determined by using the nitrogen-phosphorus detector (NPD) which was operated at 300°C while the gas flows were as follows: hydrogen 3 mL min -1 , air 60 mL min -1 and makeup (nitrogen) 10 mL min -1 . Boscalid, captan, piraclostrobin and trifloxystrobin were determined by using the micro-electron capture detector (lECD) which was operated at 280°C, and the makeup gas flow (nitrogen) was 30 mL min -1 . Agilent Technologies ChemStation Rev.B04.03 software was used for instrument control, data acquisition and evaluation.

Pesticide residue estimation
Residues of the tested substances were determined by gas chromatographic analysis, and their values were expressed in mg of AI per kg of apples. Average residue levels (n = 4), standard deviations and coefficients of variations (CV)/relative standard deviations (RSDs) of the obtained values were calculated. Subsequently, the relationship between the two variables: average residue level (y) and application rate (x) of the tested substances was established in the form of linear equation of y = ax, based on the assumption that this relationship is linear, and the straight line goes through the point (0,0) on graph (y-intercept = 0).

Method validation study
To ensure quality of the analytical results, the method was validated before the use to determine the fungicide residues in the apple samples ( Table 1). The recovery of the fungicides was determined for samples of untreated apples which were spiked with the target compounds on two concentration levels (n = 5). The average recoveries were satisfactory because they were compliant with the European Union criteria (SANCO/12495/2011) which stipulate the average recoveries in the range 70-120 % with corresponding RSD less or equal 20 % (http://ec.europa.eu/food/plant/plant_ protection_products/guidance_documents/docs/qualcontrol_ en.pdf).
The limit of quantification (LOQ) was defined as the lowest spiking level (rounded to one significant number) for which the validation criteria were satisfied and it was equal 0.01 mg kg -1 for boscalid, captan, cyprodinil, piraclostrobin and trifloxystrobin, and 0.02 mg kg -1 for fludioxonil. Excellent linearity with the coefficients of determination (R 2 [ 0.99) was achieved for the studied pesticides when using matrix-matched standards.
Measurement uncertainty of the method was estimated based on the ''top-down'' approach by using the recovery and precision data obtained in the validation study [22]. The expanded uncertainty was calculated individually for each pesticide as twice the value of the uncertainty (k = 2, confidence level 95 %), and it ranged from 8 to 35 % (22 % on average), which was distinctively less than a maximum default value of ±50 % recommended by the European Union guidelines SANCO/12495/2011 (http://ec. europa.eu/food/plant/plant_protection_products/guidance_ documents/docs/qualcontrol_en.pdf).

Field trials
The post-harvest treatment of the apple fruits (e.g. dipping) against the storage diseases of fungal origin is not allowed in Poland. In order to keep fruit quality and to provide effective protection against storage diseases, 7-14 days before harvest (pre-harvest interval PHI), growers spray apples that are still on the trees with fungicides, to control pathogens responsible for developing these diseases. Because these treatments are performed in the orchard as the last ones before harvest, they generate the highest possible residue levels of the applied active ingredients (AIs).
To determine the overall relationship between application rate of the AI of plant protection product and its residue level in mature apples, six field trials were carried out on four varieties of dessert apples (Elise, Gloster, Golden Delicious and Lobo) using five plant protection products (Bellis 38 WG, Captan 80 WG, Merpan 80 WG, Switch 62.5 WG and Zato 50 WG) containing a total six different AIs ( Table 2). The studied AIs were: captan, pyraclostrobin, trifloxystrobin, boscalid, cyprodinil and fludioksonil, belonging to four different chemical groups.
Captan (trials 3 and 4) is a non-systemic phthalimide fungicide which is used to control many fungal diseases of fruit, mainly in apple production. Captan (Captan 80 WG, Merpan 80 WG) has a protective activity, and the main product of its degradation is tetrahydrofthalimid (THPI), which was also detected in some apples intended for manufacturing baby foods (unpublished data) and during apple processing [1]. Folpet (Folpan 80 WG) is a chemical analogue of captan, and it is also registered for use in apple orchards in Poland. Similarly to captan, it undergoes decomposition to a phthalimide degradation product [5].
In trials 2, 5 and 6, pyraclostrobin and trifloxystrobin were used. Both compounds belong to the strobilurin group.
They have systemic and protective activity [15]. In Poland, formulations containing these substances are often used. Compared with captan or folpet, these substances protect the fruit at lower doses not exceeding 100 g ha -1 . Another studied fungicides were: boscalid (an anilide fungicide of systemic action), cyprodinil (an anilinopyrimidine fungicide of systemic action) and fludioxonil (a phenylpyrrole fungicide) [7]. These three fungicides are effective at doses that can be classified as medium (200 g ha -1 ).
All the fungicides used in the field trials were applied in accordance with the label recommendations, and their application rates covered a wide range of doses (D), depending on the active ingredient. The lowest dose was used in the case of strobilurin fungicides

Application rates versus residue levels of fungicides
Data on field trials and residue levels of boscalid, captan, cyprodinil, fludioxonil, pyraclostrobin and trifloxystrobin which have been found in mature apples in the following days after their applications are summarized in Table 2. In estimation of the obtained results, we assumed that the relationship between the dose of any AI of plant protection product and its average residue level in mature apples may be well approximated by a linear function, and thus, based on obtained results, such correlation had been expressed by the following equation (see also Fig. 1): where R 0 is the average residue level (mg kg -1 ) of a given substance which may be found in mature apples just after its application at a D dose (kg ha -1 ), and the numeric value of 1.2593 represents residue level of any substance if it is applied at a dose of 1 kg ha -1 .
As can be seen in Table 3, this mathematical formula makes possible the calculation of the average residue level of any fungicide now in use in the apple orchards, as well for each substance possessing fungicidal properties for the future introduction to the plant protection programs. Moreover, knowing the initial residue level (R 0 ) of a given substance and constant rate (k) of its exponential dissipation, determined independently, we may predict its concentration at any t-time, according to the general formula: R t = R 0 9 e -kt [9] and then R t ¼ 1:2593 Â D Â e Àkt : The above formula, moreover, makes it possible to estimate residue level (R t = PHI) of a given substance just after a fixed pre-harvest interval (PHI) and to compare it with the monitoring data, to ensure the correctness of the pesticide application in crop protection.

Predicted residue levels of fungicides just after treatment versus their MRLs now in force
The reduction of yield losses caused by pests and pathogenic fungi is the major challenge for agricultural production [13]. Ideally, pesticides should maintain the pest population at lower level than that causing economically unacceptable damage or loss, but without leaving residues in the edible parts of crops [10]. This is possible very rarely, and for this reason, legally, maximum residue levels (MRLs; ''tolerances'' in the USA) must be established for particular combination of pesticide and commodity. These  Table 3 summarizes plant protection products recommended for control of diseases of fungal origin in apple orchards in Poland and their application rates along with MRLs currently in force in the European Union. As seen, the expected initial pesticide residue levels were all well below European MRLs, and therefore, they do not indicate any possibility of violation of regulatory limits. In general, our study indicates a significant discrepancy between MRLs and residue levels (R 0 ) immediately after treatments as calculated on the basis of the established relationship. We found that good consistency between the two values occurs in the case of kresoxim-methyl (MRL = 0.2 vs. R 0 = 0.13 mg kg -1 ), captan and folpet (MRL = 3 vs. R 0 \ 3 mg kg -1 ) as well as for some triazoles, namely for fluquinconazole (MRL = 0.1 vs. R 0 = 0.09 mg kg -1 in periods of the high infection pressure) and propiconazole (MRL = 0.15 vs., R 0 = 0.09 mg kg -1 ). However, even for those few substances, after passing PHI period their residues will be also significantly below the MRLs.
The MRLs of other substances belonging to the triazole group (tebuconazole, myclobutanil, difenoconazole, triadimenol, fenbuconazole, tetraconazole), strobilurines (pyraclostrobin and trifloxystrobin), anilinopyrimidines (cyprodinil and pirimethanil) as well as the dodine and dithianon seem to be overestimated from 2 to 13 times. Furthermore, the discrepancies between MRLs and the final residue levels would be even greater, when taking into account their actual declines during the PHIs.
The proposed model of estimation of the residue level of fungicides in apples may facilitate not only setting MRLs for those fungicides and creation of their coherent system  but also assessing the level of pesticide residues found in fruit and vegetables in the frame of National Monitoring Programs. Our results of the residue calculations in apples just after the treatment correspond well with the results of numerous data of pesticide residue monitoring. For example, in 2008, 141 of apple samples have been tested in Poland. Among the fungicides, the most frequently detected was captan (48.2 % samples), and the highest residue level was 1.9 mg kg -1 (the MRL value is 3 mg kg -1 ), while trifloxystrobin residue was found 0.02 mg kg -1 (the MRL value is 0.5 mg kg -1 ). The highest residue levels of cyprodinil, difenoconazole, dithiocarbamates, flusilazole, folpet and pyrimethanil were: 0.06, 0.08, 0.66, 0.09, 0.33 and 0.35 mg kg -1 , respectively, and all of them were well below the applicable MRLs [12]. These findings are also consistent with the results of other monitoring surveys conducted in Poland [11] and Czech Republic [20]. Pesticide MRLs provide an enforcement tool to ensure compliance with the registered pesticide labels and reassure that food products that move in commerce are safe for consumers. According to current registration practices, as described in Regulation ( , we proposed that certain MRLs could be reconsidered and set at lower values based on pesticide residue monitoring data from surveys on crops produced under typical commercial conditions. It must be highlighted that applying the proposed approach neither mean a significant reduction in the number of initial field trials nor eliminate the need for zonal evaluation in the pesticide registration process but could be very useful for reviewing the existing MRLs, setting temporary MRLs and setting MRLs for pesticide residues in minor crops.

Conclusions
We proposed a simplified mathematical model showing the relationship between the application rate (D dose) of a given plant protection product (PPP) and the level residue of active ingredient (AI) of the substance in mature apples. This model allows the assessment of MRL values for pesticides currently approved for the protection of apple orchards and prediction of the residue level of the substances whose application is currently under consideration, when their effective doses are known, as well as any other, which may be used in apple orchards in an undefined future. Using the proposed model, we could have assessed the current MRLs. Most of them, currently found in apple orchards protection programs in Poland, as it turned out, were fixed at too high levels. We have made a comparison between level of some pesticide residues found in the mature apples within the framework of National Monitoring Programs of residue surveys carried out in Poland and the Czech Republic, and the initial level of the residue of the same substances predicted by the model. This comparison clearly has shown that the MRLs currently in force in the EU countries require verification. In the creation of a coherent system MRL an equation R t = R 0 e -kt may be helpful, especially R t = 1.2593 9 D 9 e -kt , which takes into account the relationship between the application rate of the plant protection product and the residue level of an active ingredient in the mature apples.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://crea tivecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.