Pesticide traces in local crops of Sancti Spíritus, Cuba: risk assessment study

Vegetables and rice produced in Cuba may contain residue of pesticides that exceed the maximum residue limits (MRLs). Pesticide residues on crop samples from Sancti Spíritus province were analyzed. Based on these residue data, a risk assessment of consumer exposure was conducted. Samples collected during the 2017, 2018, and 2019 periods were analyzed using UPLC-MS/MS and GC/ECD chromatography techniques. A deterministic risk assessment was performed, using consumption data collected by the authors. Also, a probabilistic risk assessment was worked out with @RISK for the Monte Carlo simulation. A broad range of pesticides were detected in the samples. Most of the detected active ingredients were insecticides and fungicides. A considerable part of the detected pesticide residues exceeded the EU and Cuban MRLs. However, the deterministic risk assessment revealed that the consumer risk of pesticide exposure is low. Only in the maximum exposure scenario children and adults were exposed to methamidophos residues in tomatoes higher than the acceptable daily intake. When a more realistic consumption (50th percentile) and higher (95th percentile) consumption pattern for tomatoes was assumed, the deterministic risk assessment found no chronic risk for children and adults being exposed to methamidophos present in the tomato samples. Furthermore, the probabilistic risk assessment concerning methamidophos exposure in tomato for children and adult strata showed no risk. To perform in the future a more accurate risk assessment for Cuban consumers it is important to have more detailed consumption data and an expanded dataset of the residues in different crops.


Introduction
The use of synthetic pesticides worldwide has become a basic need for some crops to ensure quantity and quality (Kim et al. 2016;Souza Chaves et al. 2017). Pesticides have been a key factor in the fight against many plant diseases which has allowed access to high-quality food for a broader section of the population (Alonso- Hernández et al. 2015). It has been estimated that without synthetic pesticides a large amount of food production in the world can be lost (Räsänen et al. 2015). Despite the fact that pesticides are indispensable in current conventional agricultural practices, detectable levels of residues are present on the majority of foods in stores (JMPR 2015). Due to the adverse effects of certain active ingredients (AIs) on humans health and the environment (Wahid et al. 2017), the pesticide is therefore subject to the regulations (European Food Safety Authority 2014; Jeger et al. 2018;Markantonis et al. 2018).
In Cuba, to increase the productivity of agricultural systems, technological packages have been introduced whose main component is the use of synthetic pesticides (Fernandez et al. 2018). In the province of Sancti Spíritus, where agriculture is the main economic sector, its management is in line with the rest of the country. As described by Damalas and Koutroubas (2018) in their studies on agricultural development in developing countries, the need to increase yields of priority crops in the province to reduce imports led to the use of synthetic pesticides.
In a brief search was found that grain (rice), vegetables (mainly tomato), and tobacco are prioritized crops in the province agriculture, which receive a high amount of agrochemicals, i.e. synthetic pesticides (Fernandez et al. 2018;Perdomo Hernández et al. 2016). Coupled with the use of synthetic pesticides, mainly in fruits and vegetables, there is a constant concern in the local population regarding the risk to human health and the environment reflected in various journalistic studies. However, no recent study has been carried out assessing the content of synthetic pesticide in crops produced in Cuba farms, nor a health risk assessment has been conducted.
Other old studies are focused on the adverse effects of pesticide poisoning and accidents, such as the one presented by González Valiente and Conill Díaz (1999), who report acute poisonings caused by pesticides in Cuba. González Valiente and Conill Díaz (1999) reported a total of 629 deaths due to pesticide exposure for the period 1990-1994 in Cuba. These figures accentuated the discrete upward trend observed since 1987 (González Valiente and Conill Díaz 1999). Organophosphate insecticides were the main class of pesticides responsible for these accidents. The population groups with the highest risk of death from exposure to organophosphates during agricultural activities were boys and adults. Then in 2001, González Valiente et al. made a second Cuban report on these pesticide incidents. They reported the death of 576 people due to the use, handling, suicide, or accidental poisoning (confuse the container with a food container) for the period 1995-1997. In this second report, again adult males were the most affected population group (63%), and organophosphates the main class of pesticide involved.
Under the assumption that vegetables and rice produced in Sancti Spíritus province in Cuba contain amounts of pesticide residue exceeding the established norm, constituting a risk to the health of consumers. Pesticide residues on crop samples from Sancti Spíritus were analyzed. Based on the residue data, a risk assessment of consumer exposure to residues in Cuban vegetables and rice was conducted.

Materials and methods
Sampling Samples ( Table 1) were collected in August 2017 (season with low vegetable production), January 2018 (season with high vegetable production), and May 2019 (season with medium vegetable production). Sampling took place in three municipalities with high agricultural activities (La Sierpe, Sancti Spíritus [same name as the province] and Cabaiguán, Fig. 1). The vegetable samples from Cabaiguán were obtained from the farms during harvest and from the local market. Around Sancti Spíritus there are two sample sites, La Quinta and Banao. All farms that had the vegetables of interest planted were visited. Within each sampling field (about 1 ha), five individual plants (collecting two fruits per plant) were sampled using an "X" sampling procedure and then pooled to obtain an individual composite sample. In the case of the local market, five tomatoes per each sales point visited were collected.
La Sierpe is a broad municipality with large rice fields. The rice samples for this study were taken mainly from three active places; Mapo, Peralejo, and Las Nuevas. From each one, a combined sample of 1 k was collected, taking five random spikes (following an X principle) per plot (a quarter of a hectare) from the cultivated fields (10 ha). Due to scarce rainfall and thus accumulated water, no rice was planted in the sampling area in 2017 and 2019, therefore no samples could be collected in those years. However, the consideration of studying rice was maintained since it is one of the staple foods on the Cuban menu.

Sample preparation
Once the samples were taken, they were prepared to be stored and to be transported to the Laboratory of Crop Protection Chemistry at Ghent University, Belgium, for analysis. Vegetables were cut into small pieces (~1 cm 3 ) and stored in 50 ml centrifuge tubes and frozen (− 20°C). Rice (paddy) was also kept frozen. Once in the laboratory vegetable samples were crushed and homogenized with an Ultra-Turrax IKA T18 Basic (Metrohm N.V, Antwerpen, Belgium) apparatus. Rice samples were milled using a tabletop hammer mill (Glen Creston Stanmore 14-500S, London, England) and sieved by a 1.5 mm mesh.

Sample analysis Pesticide extraction
The QUEChERS method (Quick, Easy, Cheap, Effective, Rugged, and Safe), a highly beneficial analytical approach that simplifies the analysis of multiple pesticide residues in fruit, vegetables, and cereals (Masiá et al. 2015;Rajapakse et al. 2018) was used. QUEChERS offers high recoveries for a very wide range of pesticides belonging to different chemical classes (Donkor et al. 2015;Kiljanek et al. 2016;Rajapakse et al. 2018). Two grams of dry milled rice and 10 g of vegetables of a homogenous made sample were weighed in standard centrifuge tubes (50 ml). Milli-Q water was added until the dried samples of rice weighed 10 g, then 15 ml of ACN was added and this mixture was vigorously hand-shaken for 5 s. The following salts were added to each sample to remove coextracted contaminants: 1.5 g NaCl, 1.5 g C 6 H 5 Na 3 O 7 · 2H 2 O, 0.750 g C 6 H 6 Na 2 O 7 · 1.5H 2 O and 6.0 g MgSO 4 . Samples were shaken for 5 min by 300 rpm and centrifuged by Eppendorf 5804R Refrigerated Centrifuge (VWR International GmbH, Darmstadt, Germany) for 5 min at 10000 rpm. The solvent exchange is different for the UPLC-MS/ MS (polar) and GC-ECD (nonpolar) samples. For the UPLC-MS/MS samples, 1 ml of the upper layer was sampled and added to a volumetric flask of 10 ml. Milli-Q water (9 ml) was added to obtain a total volume of 10 ml. A subsample of ±1.5 ml was pipetted in a UPLC-MS/MS vial. For the GC-ECD samples, 5 ml of the upper layer was sampled to an evaporation bowl. The solvent (ACN) was evaporated in the rotary evaporator (Buchi SL 200, Essen, Germany) and 5 ml of n-hexane was added to the bowls to redissolve the analyte. A subsample of ±1.5 ml was pipetted in a GC-ECD vial.

Ultra-performance liquid chromatography operating conditions
A Waters ACQUITY UPLC™ (Zellik, Belgium), equipped with a quaternary pump and triple quadruple system with electrospray ionization (Waters Xevo® TQD) to perform sample analyses were used. The separation column, an ACQUITY UPLC BEH C18, 130 Å (1.7 μm × 2.1 mm × 50 mm) was kept at 40°C. A sample (10 μl) was automatically injected. The mobile phase components were (A) Milli-Q water with 0.1% formic acid and (B) ACN with 0.1% formic acid. A flow rate of 0.4 ml min − 1 of 98% mobile phase A for 0.25 min was used as a gradient set. From 0.25 min to 7 min, a linear gradient was used up to 98% mobile phase B, held for 1 min. Then a linear gradient was used to 98% mobile phase A and held for 1 min. The capillary needle was maintained at + 2 kV, curtain gas (N 2 ) at 7 bars, and temperature 500°C. The AIs were monitored and quantified using multiple reactions monitoring (MRM). Two different m/z transitions were selected for each analyte. The MS/MS-transitions, ionization mode, cone voltage, and collision energy are given in Supplementary Materials 1.

Gas chromatography with electron capture detection
An Agilent Technologies 6890 N gas chromatograph equipped with an Agilent Technologies 7683 Series autosampler injector, coupled to an electron capture detector (GC-ECD) was used. Separation was performed on a HP-5MS (5% phenyl methyl siloxane) capillary column (30 m 0.25 mm 0.25 μm). As operating conditions, the column was initially set at a temperature of 60°C and then the oven temperature was increased at a rate of 20°C min − 1 to 150°C. After that, it was increased at a rate of 15°C min − 1 to 250°C, held for 2 min at 250°C, followed by an increase at a rate of 30°C min − 1 to 270°C and held constant for 10 min at 270°C. Thereafter, it was increased at a rate of 30°C min − 1 to 280°C and finally, it was held at 280°C for 11 min. Injector and detector temperatures were maintained at 200°C and 250°C, respectively. Helium was used as a carrier gas at a flow rate of 1.1 ml min − 1 and the injections were made in the split mode with a split ratio of 52.7:1.
To test selectivity, individual standard solutions, and mixed standard solutions (5 for UPLC-MS/MS and 1 mixture for GC-ECD) were injected. The accuracy was evaluated using the spiked-placebo recovery method (Bernal 2014;DG SANTE 2018). Vegetables (tomato, onion, cucumber, and pepper) and rice were spiked eight times at a high (100 μl × 10 mg l − 1 ) and a low (100 μl × 1 mg l − 1 ) concentration to determine the recovery at the maximum and the minimum amount of the AIs expected in the samples. Linearity was determined using linear regression of the calibration curve from five concentration levels (0.004, 0.01, 0.02, 0.04 and 0.1 mg l − 1 ). The obtained recoveries were expressed in percentage. Values were corrected if the recoveries were below 70% or above 120%.
The limit of quantification LOQ and detection limit LOD were calculated by multiplying the standard deviation from the eight blanks replicates by 10 and 3, respectively, using the t 99sLLMV (lowest limit of method validation) method (Bernal 2014). LOQ was set up as 0.001 mg l − 1 for the AIs analyzed by the UPLC-MS/MS and 0.01 mg l − 1 for those analyzed by the GC-ECD, and LOD was 0.0003 mg l − 1 and 0.003 mg l − 1 respectively. The pesticide residues detected in the samples were compared with the MRL values. The Cuban MRLs were obtained from the Cuban National Bureau of Standards (2012) and the European MRLs come from the EU Pesticides database (European Commission 2016).

Dietary survey
To perform a risk assessment of consumer exposure to pesticide residues in vegetables and rice grown in Sancti Spíritus, it is necessary to know, in addition to the residue data on these crops, the intake levels of the crops by the local population. A descriptive consumption data of vegetables and rice were collected at random from 462 people in different neighborhoods of the province. The equation of Cochran (1977) was applied for the calculation of this sample size. Data were separated by strata (six age classes): toddlers (0-4 years old), children (5-14 years old), adolescents (15-24 years old), adults (25-64 years old), elderly (65-76 years old) and very elderly (77 years old and over) according to Cuban standards (ONEI 2019). Intake data were collected from a two-time 24-h recall interview conducted by trained interviewers through a face-to-face interview. This concerns a dietary recall method used to collect consumption data, which involves asking the consumers to recall the actual food intake for the past 24 h or the previous day. This recall was made during two non-consecutive days in 2 weeks. The average of both days of recall was used to estimate the daily vegetable and rice consumption. The consumption data results are shown in Supplementary Materials 2. In the calculation of the amount of food from the portion sizes and recipes declared by the consumer in the survey, the following assumptions were considered: Average weight (the weight of five crops divided by five): sweet pepper 2.5 ± 0.2 g, onion 10 ± 0.8 g, cucumber 350 ± 26.0 g, and tomato 108 ± 5.0 g. Onion recipe factor = 0.067 (average of bean soup consumed (200 ml) by one person, divided by the total volume of soup prepared (3000 ml)). In the case where the onion was consumed in another receipt, the value of its mass described above was taken into account. In all cases, it was multiplied by the total number of onions declared as consumed by the interviewee and in its variants. Sweet pepper recipe factor (soup ratio 0.067 multiplied by the average of sweet pepper declared as used as portion size per strata = 0.068 for toddlers, 0.13 for children, adolescents, and very elderly, and 0.27 for adults and elderly. Based on the cooked rice consumption portion size data collected, it was assumed 250 g of cooked rice as the maximum portion size for toddlers and very elderly people, 400 g was considered for children, and 500 g for adolescents, adults, and elderly. The average (rice and vegetable) daily consumption (kg food kg − 1 BW day − 1 ) was calculated from the frequency of daily intake (kg food day − 1 ), portion sizes, recipes or processing, and the kg body weight (BW) per day for each person.

Risk assessment
The consumption data and the residue concentration data were used to assess chronic consumer exposure to pesticides. The residue levels of some of the individual AIs in the vegetable samples were low, below the LOQ. Therefore, for the chronic exposure assessment a pessimistic approach was adopted (worst case scenario). This means that for the AI selected for the assessment, a LOQ value was used for the non-quantified residues and LOD for the non-detects. The estimated exposure was compared to ADI (acceptable daily intake) of the AI resulting in a ratio. To estimate the consumer's level of safety (Mekonen et al., 2015) the ratio was evaluated against the reference value of one. If the ratio was higher than one, it indicated that ADI was exceeded for that particular AI and that in turn there is a potential risk of exposure to the AI. A cumulative scenario was evaluated for AIs with the same mode of action, detected in the same crop or several of those evaluated, adding their exposure value per strata (toddler, children, adolescents, ...), per year.
With each selected AI, one exposure was assessed for each of the six age classes (scenario), based on three sub-scenarios: minimum (least realistic), medium (realistic), and maximum (worst case). The deterministic exposure was carried out based on a single point estimation, minimum exposure scenario (P5), median exposure scenario (P50), and maximum exposure scenario (P95) using Eq. 1. In this regard, the minimum, median, and maximum consumption values of vegetables and rice (kg kg − 1 BW day − 1 ) were respectively multiplied by the minimum, median, and maximum concentrations of the AI detected as in Mekonen et al. (2015).
The detected pesticide residue values were corrected according to the appropriate process factors (Table 2). For example, rice is cooked, onion and sweet pepper are cooked in soup or meat sauce and tomato and cucumber are washed (Holland et al. 1994;Keikotlhaile et al. 2010;Liang et al. 2013;Markantonis et al. 2018;Scholz et al. 2016;Scholz et al. 2018;Tomer and Sangha 2013). Based on the answers of the interviewees and the preference of the population, only raw tomatoes were considered, not their variants (soup, pasta, or pure ...). The values of processing factors are mainly according to the compilation of the Netherland National Institute for Public Health and the Environment (RIVM 2015). Others source also use are the studies from Holland et al. (1994), Liang et al. (2013) and Tomer and Sangha (2013) , In the case of obtaining the deterministic assessment a result that requires a more accurate analysis, a probabilistic assessment has been performed. The probabilistic analysis will follow the same methodology as the deterministic analysis, with each AI selected, independent of the crop, and exposure assessment will be made for each of the six age classes (scenarios), based on three subscenarios: minimum, medium, and maximum. In this way, individual risk can be assessed. Additionally, other analyses such as adding up the intake of various crops with residues of the same AI per age classes or adding up the exposure of AIs with the same MoA will also be performed. The probabilistic assessment of chronic risk is conducted with @RISK 7 (Palisade Corporation, USA), an add-in for Excel (Jeger et al. 2018;Mekonen et al., 2015), in which the consumption (kg food /kg BW day) and residue (mg A.I /kg food) distributions are combined into an exposure distribution (mg AI /kg BW day). In contrast to the deterministic approach, in the probabilistic assessments, suitable parametric distributions are fitted to the residues and the consumption data. Further information regarding @risk software and its operation is described in Supplementary Materials 3.

Statistical analysis
The absolute and relative frequencies expressed in percentages of each of the categories of the described variables were calculated. The Statistical Package for Social Sciences (SPSS) program (IBM Corp,v. 20. Armonk, NY, USA) was used. The Paired Samples T-test (p < 0.05) was used to evaluate differences in pesticide group between years, and the significant differences between the population strata (toddlers, children, adolescents, adults, elderly, and very elderly) concerning pesticide exposure.

Pesticide residue in crops
The residues range found using the LC-MS/MS and GC-ECD technique is given in Tables 3 and 4 respectively, full tables are in Supplementary Materials 4 and 5. The corresponding MRLs are shown at the bottom of the tables. In 18% of the analyzed samples (n = 106) no residue was detected. The majority of the samples (27%) contained three different AIs. At least one AI was detected in 14% of the samples, and up to six different AI residues were detected in 8% of the samples.
From the herbicides, only residues of one AI (alachlor) were detected. No significant differences were found between the frequency of the detection of fungicide and insecticide residues present in 2017 (p = 0.389) and 2018 (p = 0.680) vegetable samples. However, in 2019 there was a significant difference between the number of insecticides and fungicides residues detected (the difference does not include 0 (95% CI), therefore the p-value of the T-test is < 0.05). Secondly in the T-test analysis between the total concentrations of fungicides and insecticides per year, no significant differences were found (2017 p = 0.102, 2018 p = 0.764, and 2019 p = 0.979). The study with these results demonstrates that regardless of the year, total fungicide and insecticide residues are equally present in the vegetable samples. In 2017 the organophosphates methamidophos, acephate, and the organochlorine endosulfan, all have been banned for more than a decade in Europe (FAO/UNEP 2007, were detected. In 2018 the presence of methamidophos and endosulfan residues persisted. Lindane as well as alachlor, also banned in Europe, was detected that year. In 2019 residues of acephate were found again, and of two other AI banned in Europe (dimethoate and fenpropimorph). Also, cypermethrin, pyrimethanil, imidacloprid, tebuconazole, and azoxystrobin were detected in samples from all three sampling periods. Tomatoes contained most pesticide residues. In vegetables, 49% of the pesticides were detected in tomatoes (139) and 26% in onion (73) (whole data in Supplementary Materials 4 and 5).
It was also recognized by farmers that tomato and onion cultivation results in the greatest economic benefits. Also, farmers mentioned in the interview that motivated by the great economic benefit they overuse pesticides to guarantee the highest possible yield.   In certain cases, pesticide residues are detected but are below the limit of quantification (0.001 mg kg-1) and are therefore reported as "<LOQ". " a " Indicates lower limit of analytical determination (default value). This study finding high levels for some products, confirms that synthetic pesticides are still in use in vegetable production in Cuba. Of the samples, 31% contained at least one AI above the European MRL (most Cuban MRLs are missing, which difficulty a national analysis). This is despite previous Cuban studies which have warned of the adverse consequences of synthetic pesticides on human health and the environment (Figueroa González and Pérez-Consuegra 2012; Hernández Contreras et al. 2010; Hernández Núñez and Pérez-Consuegra 2012; Hernández et al. 2003aHernández et al. , 2003bPerdomo Hernández et al. 2016;Vega Bolaños et al. 1997).
In the onion samples, the following AIs were quantified: alachlor, pyrimethanil, chlorothalonil, azoxystrobin, cypermethrin, bifenthrin, imidacloprid, methamidophos, lindane, and endosulfan. Six of the ten AIs are insecticides; three are fungicides and one is the herbicide alachlor. Endosulfan and lindane have the same mode of action (MoA) and this also applies to bifenthrin and cypermethrin (IRAC 2019). They should be alternated with AIs with another MoA to prevent resistance development (IRAC 2019; Sparks and Nauen 2015).
In tomato samples, the following AIs were quantified: pyrimethanil, azoxystrobin, tebuconazole, imidacloprid, methamidophos, difenoconazole, triadimenol, bifenthrin, chlorothalonil, cypermethrin, and fenpropimorph. Seven out of the eleven AIs are fungicides, the others are insecticides. Except for methamidophos, the other AIs present in the samples were lower than the established MRLs. Due to the warm climate of the province of Sancti Spíritus (annual average 25.7°C, summer 30.1°C, and winter 21.0°C) with high relative humidity (annual average 78%), the pressure of fungal infection is high (Harel et al. 2014;Rawat 2015;Shamshiri et al. 2018). This is one of the reasons why farmers prefer to grow vegetables during winter (Ayyogari et al. 2014;Hernández Córdova and Soto Carreño 2012). It is also important to draw some attention to the fact that tebuconazole, difenoconazole, and triadimenol have the same MoA (FRAC 2018). Methamidophos is the only AI found in tomatoes that is no longer allowed to be used in Europe and the USA (FAO/UNEP 2015; Roberts and Routt Reigart 2013) and its use is also no longer allowed in the updated Cuban official list of authorized pesticides from 2016 (Registro Central de Plaguicidas 2016). Methamidophos was allowed to be used until 2016; the presence of methamidophos in the samples (from 2018) demonstrates unauthorized use. No previous reports of pesticide residues in Cuban tomatoes were found in the last decade. Before, Hernández et al. (2003aHernández et al. ( , 2003b have reported an average concentration of methamidophos of 67 μg kg − 1 in tomatoes, similar to the average value quantified in this study (73 μg kg − 1 ). Hernández et al. (2003) also have reported average values of dimethoate and carbaryl of 75 μg kg − 1 . In this study, dimethoate (<LOQ) was found in tomatoes and cucumbers. Furthermore, Hernández et al. (2003aHernández et al. ( , 2003b have reported average values of methamidophos in sweet pepper (270 μg kg − 1 ).  In certain cases, pesticide residues are detected but are below the limit of quantification (0.01 mg kg-1) and are therefore reported as "<LOQ". " a " Indicates lower limit of analytical determination (default value). In 2017 residues of two samples (18%) exceeded the MRL. This was observed for cypermethrin in onion (Cuban National Bureau of Standards 2012; JMPR 2015). For 2018 a higher percentage (65%) of exceedances was found. Then in 2019, the percentage of exceedances decreases again to only one AI, alachlor in cucumber. Alachlor is approved for use in Cuba (Registro Central de Plaguicidas 2016) on corn, beans, soya beans as well as on tobacco, but not in vegetables where it was found in this study. No MRL for alachlor is registered in the Cuban MRL database (Cuban National Bureau of Standards 2012). The residue of methamidophos and lindane, AIs which are banned for use in Europe, exceeded the European MRLs. The Cuban norm which registers the pesticide MRLs in food dates from 2012 (Cuban National Bureau of Standards 2012), and although both methamidophos and lindane were authorized for use at that time, (Lista Oficial de plaguicidas autorizados. Republica de Cuba 2008), no MRLs for methamidophos and lindane in vegetables were established (Cuban National Bureau of Standards 2012). As part of the national program to reduce the use of highly toxic AI, in the latest list of authorized pesticides in Cuba, dating from 2016 (Registro Central de Plaguicidas 2016), methamidophos and lindane do not appear anymore. In contrast to what was found in the samples, where it is evident that farmers maintain stocks of these products or due to its persistence in the soil some of these residues can still be detected today.
Endosulfan is still present on the list under restricted authorization use, mainly for coffee. It can also be used in onion, garlic, cucurbits, tomato, and sweet pepper. In the case of chlorothalonil, the residues quantified in onion samples exceeded the corresponding European MRL limit value (0.01 mg kg − 1 ). This was not the case for the Cuban MRL values (0.5 mg kg − 1 ) (Cuban National Bureau of Standards 2012). The same was observed for bifenthrin. Pyrimethanil was present in 2018 samples. In the region of Sancti Spíritus, pyrimethanil is authorized for use in banana crops to treat the Sigatoka spp. (Mycophaerella spp.), nevertheless, it is not authorized in vegetables.
In addition, few AIs (two) were found in the paddy rice and their concentrations were low (ten times lower than the corresponding MRL). The detected results were in correspondence with those reported by Teló et al. (2017) in Brazil. Also in previous studies from Barik et al. (2010) and Teló et al. (2015), no residues of imidacloprid and cyhalothrin were found in the rice grains cultivated under GAP.
The table is larger than the page length, the table is placed on page 34, as journal guidelines suggest. The table should appear in this location in the text file during production.
The table is larger than the page length, the table is placed on page 35, as journal guidelines suggest. The table should appear in this location in the text file during production.

Consumer risk assessment
As seen in Section 3.1, 25% of the quantified pesticide residues exceed the European and Cuban MRLs. The MRL is a value used to determine if a product can be placed in the market or not according to Good Agricultural Practices (it is related to food trade). It is not a toxicological safety limit such as ADI or ARfD, which indicates a health concern (Brancato et al. 2018;JMPR 2015). Therefore as long as exposure is estimated to be below the ADI or ARfD there is no perceived risk for the consumer. The vegetables and rice produced in Cuba are only for local consumption and do not go to foreign markets; therefore exceedances of the MRL in Cuba do not imply a real trade problem.
Vegetable consumption (Supplementary Materials 2) is low, only 42.4% of the surveyed population (average of both interview days) declared to have consumed at least one of the four vegetables evaluated (prevalent vegetables in the region). Mulet Pascual (2017) pointed out that the vegetable price increase is one of the main causes of the vegetable consumption decrease (1 k of tomatoes, cucumbers, onions, and sweet peppers respectively purchased in-state establishments represents 3.9% of the average salary "752 pesos" (Gaceta Oficial de la República de Cuba, Ministerio de Finanzas y Precios 2016)). An average household spends between 74 and 79% of their monthly income only on food (García Álvarez et al. 2014). Besides, the statistical yearbooks indicate a decrease in the total sales of vegetables from 97.9% (total sales in 2017 respect to 2016) to 60.1% in thousands of pesos for the period 2018-2019, in the same sampling period (January-June) (ONEI 2017(ONEI , 2019. Furthermore, the residue dataset is limited, due to not all possible IAs having been explored, for example, those used in fruits, other grains, or tobacco (Perdomo Hernández et al. 2016) or temephos in potable water (George et al. 2015). The accumulated risk of all the present AIs can imply a higher risk for consumers.

Deterministic
Residue data of 15 from the 18 AIs detected (Supplementary Materials 6) and consumption data of 462 people (Supplementary Materials 2) were used in the deterministic chronic dietary exposure assessment. The rest of the AIs (triadimenol, acephate, and dimethoate) were not considered, as there were more undetected (< LOD) than detected results. Further on, for several of the AI selected for the deterministic consumption exposure assessment (six different age classes: toddlers, children, adolescents, adults, elderly, and very elderly), only their maximum percentile (p95) have been considered in the evaluation since their median (p50) and minimum (p5) consumption values are zero and residual data concentrations are below the LOQ. These cases include azoxystrobin and tebuconazole in 2017, alachlor, endosulfan, chlorothalonil, tebuconazole, bifenthrin, imidacloprid, and cyhalothrin in 2018, and propiconazole, chlorothalonil, and cypermethrin in 2019. Figure 2 shows the results of the consumer deterministic chronic risk assessments (ADI exceedance), for the samples collected in 2017, 2018, and 2019. Since most of the risk values obtained are very low, it was decided to show in Fig. 2 only those that were equal or higher than 0.01 (only observed in maximum scenarios). Therefore, all the results of the minimum and median scenario evaluations, as well as many other maximum scenarios, from the six different age classes studied, were left out of Fig. 2. The detailed risk assessment results derived from the exposure analysis for each scenario with its ADI reference value, depending on the AI and the crop assessed, are provided in the tables of Supplementary material 7.
The T-test analysis performed showed that there were no significant differences (p = 0.555) between the strata (toddlers, children, adolescents, adults, elderly, and very elderly) concerning pesticide exposure via vegetable consumption, with an exception of methamidophos in tomatoes. The rest of the scenarios (other AIs and analyzed crops) showed that the evaluated strata are out of risk (exposure [minimum, medium and maximum] value below ADI indicated no health concern) (Fig. 2). Even if the risk indicators of the AIs with the same MoA are added in a combined diet (cumulative), such as is the case for endosulfan and lindane in onions and cucumbers, the pyrethroids cypermethrin, bifenthrin and λ-cyhalothrin in tomato, onion, and rice, chlorothalonil in tomato and onions, difenoconazole and propiconazole in sweet pepper and imidacloprid in sweet pepper and cucumber, there is a negligible human exposure risk. Figure 2 shows that only for the maximum exposure scenario of children, and adults a risk to methamidophos is seen due to the consumption of tomatoes from the 2018 period samples; children 2.3, and adults 1.2 times above the ADI. If it is considered that a maximum exposure scenario is the worst case (see Supplementary Materials 2 crop consumption), very few consumers (3 children and 17 adults) are exposed to higher levels than the ADI.
The hypothesis was that the average exposure scenario should be more realistic. However, due to the limitations of the low number of respondents and low consumption (more than 50% of the respondents did not consume vegetables) (Supplementary Materials 2), it was shown that it was not realistic enough.
To be able to assess realistic chronic exposure and risk for children and adults exposed to methamidophos residues in tomatoes, one assumption was made: increased tomato consumption. This is finally an additional exposure scenario; median exposure to hypothesized higher tomato consumption. To reach the median toxicological Fig. 2 ADI exceedance (chronic risk) of Cuban consumers (six different age classes: toddlers, children, adolescents, adults, elderly, and very elderly) for selected pesticide residues in crop samples collected in the period 2017-2019. Risk is the quotient between the value of the exposure and its reference (ADI). If the quotient is higher than 1, it means that the limit has been exceeded and indicates that a certain risk exists. ADI exceedance was calculated with deterministic modeling. "n" means the number of person threshold for children at an average residue concentration of 0.025 mg kg − 1 it is necessary to consume (mathematically) at least seven tomatoes daily. In the case of adults, thirteen tomatoes would be required. These are consumption data that are higher than real consumption. However, the minimal chronic risk is obtained when children consume daily 1.5 tomatoes and adults 2.6 tomatoes. Only one adult (0.4%), claimed having consumed that amount of tomatoes (three tomatoes). The deterministic assessment of dietary risk from rice consumption (Table 4) showed that there is a negligible risk for every stratum. The results were as expected, based on the results from other researchers, who mentioned that some pesticides can be removed through milling, cleaning, washing, and cooking processes, which leads to reduced exposure (Đorđević and Đurović-pejčev 2016;Shakoori et al. 2018;Teló et al. 2017). This result may depend on the level of the traces present and if GAP is respected (Teló et al. 2015).
A study developed by Lehmann et al. (2017) in Burkina Faso identified a risk for λ-cyhalothrin (five times higher than the concentration found in this study), as well as for other AIs, in acute and chronic exposure scenarios for the consumption of vegetables. In the same study, they mention that the processing factor largely influences the risk, suggesting that a simple washing of vegetables with water can considerably reduce the risk of hazardous exposure. Other authors also mention how household processing (washing, peeling, and cooking) removes pesticide residues (bifenthrin, diafenthiuron, imidacloprid, and emamectin benzoate) from vegetables (Keikotlhaile et al. 2010;Scholz et al. 2016;Sheikh et al. 2012;Sheikh et al. 2015;Skovgaard et al. 2017;Tomer and Sangha 2013). In this way, processing factors on pesticide residues for vegetables and cereals play a very important role in risk assessment (Boitshepo Miriam Keikotlhaile 2011).

Probabilistic
Based on the fact that only methamidophos residues in tomato exerted risk for children and adults population in the maximum exposure scenario of the deterministic approach. Furthermore, there was only a small fraction of the population (< 20%) that consumed tomatoes. Also, the only low residue of methamidophos was detected in the tomato samples. Because of this, it was decided to carry out a more accurate and realistic probabilistic assessment analysis of the aforementioned strata.
The residue data of methamidophos could not be fitted directly to distribution, because of the many LOD values (55%). Hence, the data were grouped into two groups, values above and below of LOD values. The first group was fitted to a distribution, following the rules of selection stable by the method and realistic according to the data. In this case, the most realistic distribution of the data was Invgauss, the option, Truncate minimum to 0, was also used, to avoid that the function will move to negative values. After which the IF function in @RISK was used to generate a random distribution for the entire residue data.
The IF functions for the methamidophos residue distribution was: IF(RAND() < fraction of LOD,LOD,methamidophos residue distribution).
The best-fitting distribution was also chosen for adult tomato consumption data. The data was divide into two, the percentage of zero consumption (84%) and the adult tomato consumption values, with a response to Weibull distribution. The IF function used to generate the random distribution for the entire consumption data corresponds to the following function.
In the case of children the IF function formulates is: IF(RAND() < fraction of 0,0,average consumption value).
Since only 5 children declared to have consumed tomatoes (low number of consumers) it is very limited to make a distribution. Instead of making a distribution, we opted for a fixed value (average) and then make an IF function.
From the probabilistic assessment of chronic exposure, shown in Supplementary Materials 8, it was found that unlike the results of the deterministic test, the daily intake of methamidophos (residue in tomatoes) was lower than the ADI in adults and children is very low. It only compromised 0.2% of that population, for average consumption. This implies that the dietary risk for the evaluated strata can be negligible. Considering that the probabilistic analysis is more precise than the deterministic one, it can be assumed that there was no risk at all from the consumption of tomatoes with methamidophos residues at the concentrations found and reported tomato consumption.
Out of these results, based on this case study in Sancti Spíritus, Cuba, and because the overall consumption of vegetables studied is low, it appears that the possible risk of pesticide residues in vegetables and rice is quite low. However, the work presents some weaknesses, which may also be responsible for the results obtained and which should be corrected in future research, to arrive at more realistic conclusions. The main weakness is the few variants and forms of consumption of the vegetables evaluated (e.g., tomato only raw), and the low number of foods studied.
To account for the study limitation, it could be proposed to evaluate the exposure of a typical Cuban daily recipe. Also consider different scenarios such as rural and urban, as well as seasonal crops. In this way, it would be possible to have a better idea of which crops, per season, in which stratum, and which active ingredients would require more attention, to keep consumers out of risks. Regardless of ensuring that good agricultural practices are followed, monitoring MRL values and advising farmers on the toxicity of the active ingredients they use and how to use them.

Conclusions
Most of the detected active ingredients (AIs) in the analyzed samples were insecticides and fungicides. Herbicides were almost absent. One-quarter of the quantified residues exceeded their corresponding maximum residue limit. In the deterministic assessment, only methamidophos in tomato samples for the maximum scenario in children and adult strata show a potential risk. No risk was found for the rest of the scenarios and strata. When more representative assumptions (median) and higher tomato consumption were assumed, the deterministic assessment found no chronic risk in children and adults due to methamidophos present in the tomato samples. The probabilistic risk assessment concerning methamidophos exposure in tomato for children and adult strata also showed negligible risk. Although the risk of exposure was not prominent, the presence of AIs at concentrations above the MRLs remains a problem. In the case of any government intervention where increased consumption of vegetables is considered, the use of pesticides in quantities exceeding those allowed could pose a problem. However, several weaknesses in the study limit a more accurate result. Such as few foods studied and consumption variants, like the type of food (raw or processed), of those studied. To improve a further accurate risk assessment for Cuban consumers it is important to have data on a more regular diet and food consumption varieties. In further research, this data could be collected so that the risk assessment is repeated. As seen in interviews with farmers, some pesticides that are frequently used in Cuba are not investigated in this work. In the future, exposure could be better determined when the dataset includes the consumption of fruits, grains, roots, and tubers and drinking water (which is treated frequently with temephos).
Additional file 1: S1. MS/MS-transitions, applied collision energy and dwell time for the active ingredients analyzed in the study. S2. Vegetable and rice consumption data for the evaluated strata obtained from the consumption survey and used in the deterministic and probabilistic risk assessment. S3. Technical explanation of @Risk. S4. Vegetable and rice samples with pesticide residues detected by LC-MS/MS in 2017, 2018, and 2019. In certain cases, pesticide residues are detected but lay under the limit of quantification (0.001 mg kg-1) and are therefore reported as "<LOQ". "*" Indicates lower limit of analytical determination (default value). "CoV" indicates the coefficient of variation. "e" Indicates exceedances of MRL values. S5. Vegetable and rice samples with pesticide residues detected by GC-ECD in 2017. In certain cases, pesticide residues are detected but lay under the limit of quantification (0.01 mg kg − 1 ) and are therefore reported as "<LOQ". "*" Indicates lower limit of analytical determination (default value). "CoV" Indicates the coefficient of variation. " e " Indicates exceedances of MRL values. S6. Residue data of the active ingredients selected for the deterministic consumption exposure and risk assessment. S7. Deterministic risk assessment for the chronic exposure of the consumer for samples collected in the period 2017-2019. S8. Probabilistic assessment of chronic risk for children and adults due to the consumption of tomato and its methamidophos residue. Exposure and ADI values are defined. The ADI of methamidophos (0.001 mg kg − 1 BW day − 1 ) is set on top of each graph. To the right of the ADI value, it can be found the percentage of the population exposed to risk.