Skip to main content

Distributions and determinants of urinary biomarkers of organophosphate pesticide exposure in a prospective Spanish birth cohort study



Prenatal exposure to organophosphate pesticides (OPs) has been associated with impaired child development. Pesticide exposure determinants need to be studied in order to identify sources and pathways of pesticide exposure. The aim of this paper is to describe prenatal exposure to OPs and evaluate the associated factors in pregnant women.


The study population consisted of pregnant women (n = 573) who participated in the INMA birth cohort study in Valencia (Spain, 2003–2006). OP metabolites were analyzed in maternal urine at the 32nd week of gestation using a liquid chromatography-high resolution mass spectrometry method. The analysis included non-specific (diethyl phosphate [DEP], diethyl thiophosphate [DETP], dimethyl thiophosphate [DMTP], dimethyl dithiophosphate [DMDTP]) and specific metabolites (2-diethylamino-6-methyl-4-pyrimidinol [DEAMPY], 2-isopropyl-4-methyl-6-hydroxypyrimidine [IMPY], para-nitrophenol [PNP], and 3,5,6-trichloro-2-pyridinol [TCPY]). Information about the sociodemographic, environmental, and dietary characteristics was obtained by questionnaire. The association between log-transformed OPs and covariates was analyzed using multivariable interval censored regression.


The detection frequencies were low, DMTP and TCPY being the most frequently detected metabolites (53.8% and 39.1%, respectively). All the OP metabolites were positively associated with maternal intake of fruits and vegetables. Other maternal characteristics related to the OPs were body mass index (BMI) before pregnancy and smoking habit during pregnancy. Women with lower BMI and those who did not smoke presented higher OP concentrations. Moreover, mothers who had a yard or garden with plants at home or who lived in an urban area were also more exposed to OPs.


The OP detection frequencies and the concentrations observed in our study population were low, compared with most of the previously published studies. Given the high vulnerability of the fetus to neurotoxicant exposure, further research on the determinants of the body burden of OPs during pregnancy would be necessary. The knowledge gained from such studies would enhance the effectiveness of public health control and future recommendations in order to reduce the risk to both the health of pregnant women and the health and development of their children.

Peer Review reports


Pesticides are widely used in the agricultural areas of Spain. The sales of these products over the last decade have exceeded 600 million euros annually, thus accounting for 10% of the total sales in Europe [1]. The Valencia Region is the second largest agricultural area in Spain and one of the largest pesticide users in this country. This region was responsible for more than 12% of the total national pesticide consumption in 2009 [2]. The organophosphate pesticides (OPs) are the most widely used active substance in insecticides, followed by pyrethroids and carbamates [3]. As a result of their massive use, especially chlorpyrifos, OPs spread through the environment and contaminate water, soil, and atmosphere, resulting in a potential risk to humans and the environment [4].

Pesticides are also used in domestic settings, pyrethroids being the most common active ingredient used in residential insecticides [5]. However, some OPs, such as chlorpyrifos, were commonly found in the composition of domestic insecticides in Spain [6]. As of 2008, the use of chlorpyrifos as a domestic pesticide was phased out in the EU (Directive 98/8/CE).

Prenatal exposure to OPs was first described in three birth cohort studies conducted in USA. Chlorpyrifos was detected in cord blood samples from newborns participating in the Columbia Center for Children’s Environmental Health in New York, showing a high correlation with maternal blood, which may indicate that this pesticide crosses the placenta barrier [7]. Urinary levels of some OP metabolites have also been measured among pregnant women from the Children’s Environmental Health Study in New York [8] and from the Health Assessment of Mothers and Children of Salinas (CHAMACOS study), an agricultural area in California [9]. These latter studies found detectable levels of different OP metabolites in nearly all urine samples.

There is increasing public concern about the effects associated to the exposure to pesticides during early development [10]. Fetuses and children are especially more vulnerable to exposure to environmental pollutants in comparison to adults, since their organs and systems are still developing and their detoxification mechanisms are not yet fully mature [11]. There is increasing evidence of a relationship between prenatal exposure to OPs and child neurodevelopment [12]. In addition, OPs have also been linked in epidemiology studies to shorter time of gestation [13, 14], low birth weight [15, 16], increased child blood pressure [17], respiratory outcomes [18], obesity and diabetes [19, 20]. To date, very few longitudinal studies have investigated factors associated with pesticide exposure in pregnant women [8, 21,22,23,24].

The aim of this study is to describe prenatal exposure to OPs and evaluate the associated factors in pregnant women participating in a birth cohort study located in a Spanish area with intense agricultural activity. Studies on pesticide exposure determinants are needed to identify sources and pathways of pesticide exposure and to contribute to policies aimed at reducing exposure.


Study population

This study is framed within the INfancia y Medio Ambiente (INMA) Project (Environment and Childhood), the aim of which is to investigate the effects of environmental exposure, diet, and genetics on fetal and child development in a cohort of pregnant women and their offspring in different regions of Spain ( The study protocol has been reported elsewhere [25]. Briefly, pregnant women were recruited at the beginning of their pregnancy (10–13 weeks of gestation) in the region of Valencia (2003–2005, n = 855). These women were followed up until the third trimester of pregnancy (n = 794). The final study population consisted of 573 pregnant women with complete data on OP exposure. The main reason for the decrease in the study population was the unavailability of urine samples for OP measurement (n = 122) and the limited resources for the analysis of OPs in all available samples (n = 99).

Analysis of OP metabolites in urine

Urine samples were collected at the third trimester of pregnancy (mean ± sd: 32.2 ± 1.4 weeks of gestation). Samples were separated into aliquots of 10 mL and then frozen at −20 °C until analysis in the Public Health Laboratory of Valencia. One aliquot was used for the analytical determination of metabolites in urine by using an ultra-high pressure liquid chromatography coupled with high resolution mass spectrometry method (UPLC-HRMS) [26]. Briefly, metabolites from hydrolyzed urine were extracted by liquid-liquid extraction with 10 ml of acetonitrile followed by a QuEChERS extraction [27]. The acetonitrile layer obtained was evaporated to dryness, dissolved in 200 μL of methanol: water (10/90, v/v) containing 0.1% of acetic acid, ultra-centrifuged (11,000 rpm, 3 min and 10 °C) and transferred into an injection vial for analysis on an Accela liquid chromatography UHPLC system from ThermoFisher Scientific (Bremen, Germany). Separation was performed by using a Hypersil Gold column (100 × 116 2.1 mm, 1.9 μ) with a flow rate of 400 μL min−1 and an injection volume of 10 μL. Mobile phase components were 0.1% acetic acid aqueous solution (A) and methanol containing 0.1% acetic acid (B). Targeted mass analysis was performed on an Orbitrap mass spectrometer Exactive™ analyzer (Thermo Scientific, Bremen, Germany) operating in both positive and negative ESI modes. The ion source parameters were 3.5 kV (positive mode) and 2.5 kV (negative mode); sheath gas flow-rate: 55; auxiliary gas flow-rate: 10; skimmer voltage: 23 V; heater temperature: 300 °C; capillary temperature: 150 °C; capillary voltage: 45 V, and tube lens voltage: 120 V. The system operated at a resolving power of 50,000 (250 ms). Each metabolite was identified and confirmed following the criteria of relative retention time (RRT), mass tolerance value of the molecular ion <5 ppm, isotopic pattern, and fragment ions.

The analysis included four non-specific metabolites of exposure to OPs (Table 1): diethyl phosphate (DEP), diethyl thiophosphate (DETP), dimethyl thiophosphate (DMTP), dimethyl dithiophosphate (DMDTP), and the four specific metabolites: 2-diethylamino-6-methyl-4-pyrimidinol (DEAMPY), 2-isopropyl-4-methyl-6-hydroxypyrimidine (IMPY), para-nitrophenol (PNP), and 3,5,6-trichloro-2-pyridinol (TCPY). Sum (molar) variables were also calculated: sumDEP: DEP + DETP, sumDMP: DMTP + DMDTP, sumDAP: sumDEP + sumDMP, and sumOPs: DEAMPY + IMPY + PNP + TCPY.

Table 1 OP metabolites analyzed and possible precursor compounds

The analytical suitability of the method was checked each working day through the analysis of the quality controls previously described by Roca et al. [26]. In accordance with this quality procedure, in each analytical batch, various quality control samples (QC) were prepared by spiking blank urine samples at LoQ, QCinter, and QChigh levels to calculate the extraction efficiency and ensure a good quantification of real samples. The QCs were subjected to the same extraction and analysis procedures as real samples and calibration curve points.

OP concentrations were expressed in μg/L and in μg/g of creatinine in order to correct for urinary dilution. Creatinine was determined by the Jaffé method (kinetic with target measurement, compensated method) with Roche reagents in a Hitachi 911 auto-analyzer.


The women filled in two questionnaires during pregnancy, one at the first trimester (mean ± sd: 12.1 ± 1.5 weeks of gestation) and the other at the third trimester (mean ± sd: 31.8 ± 1.7 weeks of gestation). The questionnaires were administered by trained interviewers and focused on sociodemographic, environmental, and lifestyle information during pregnancy. The maternal covariates obtained were maternal age (years), country of birth (Spain, other), educational level (up to primary, secondary, university), body mass index before pregnancy (BMI, kg/m2), parity (0, 1, >1), employment status during pregnancy (employed, not employed), smoking habit during pregnancy (smoker, non-smoker), season of sampling (winter, spring, summer, fall), area of residence (urban, metropolitan, semi-urban, rural), use of indoor pesticides (no, yes), yard with plants at home (no, yes), outdoor pesticides application (no, yes), residence near fields or greenhouses (no, residence near fields, residence near fields sprayed with pesticides). We also obtained information about the paternal or maternal occupation and identified those related to the use of pesticides (agriculture, gardening).

We defined social class based on the maternal or paternal occupation during pregnancy with the highest social class, according to a widely used Spanish adaptation of the coding system of the International Standard Classification of Occupations approved in 1988 (ISCO88) [28]. Class I + II included managerial jobs, senior technical staff, and commercial managers; class III included skilled non-manual workers; and class IV + V, manual workers.

The Valencia Region was divided into four sub-areas according to the population density and land uses: urban, metropolitan, semi-urban, and rural.

Dietary information was collected by using a validated semi-quantitative food frequency questionnaire (FFQ) of 101 food items [29]. The FFQ was administered during the third trimester of pregnancy (the same time point as the general questionnaire) and covered average intakes from the previous 3 months. Information about the intake of different fruits and vegetables was obtained and grouped as follows (expressed as daily servings): citrus fruits (orange and orange juice), stone fruits (prune, peach, nectarine, and apricot), kiwis, apples and pears, bananas, watermelon, fruit vegetables (courgette, aubergine, cucumber, tomato, and pepper), green vegetables, (lettuce, spinach, chard, and cabbage), other vegetables (green beans, onion, carrot, and pumpkin). Total fruit and vegetable intakes were also calculated.

Statistical analysis

Maximum likelihood estimation of censored linear regression models was used to estimate the geometric mean (GM) and 95% confidence intervals (95%CI) for OP metabolites [30, 31]. The compounds were log-transformed so as to approach normality. This scheme also applies to molar summed variables. The OP metabolite concentrations were expressed as μg/g of creatinine, as μg/L, and as molar (nm/g of creatinine and nm/L). All individual compounds were included in the definition of sum variables, although, when studied separately, we only considered those with a detection frequency > 30%.

The association between OP metabolite concentrations and the sociodemographic, environmental, and dietary characteristics was evaluated by multivariable interval censored regression models in order to deal with values below the LOD. These models were built using a backward elimination procedure: first, all the covariates related to the OP metabolite concentrations at the univariate level (p < 0.1) were included in the initial models. Then, those variables not associated with the OP metabolite concentrations in the multivariable model using the Likelihood Ratio Test (p > 0.1) were sequentially excluded. For comparability purposes, final models were adjusted for all the covariates that were retained for at least one of the compounds. For this analysis the unadjusted OP concentrations were used including creatinine as a covariate.

In a further analysis, total fruit and vegetable intake variables were replaced by the subgroup dietary variables (citrus fruits, stone fruits, kiwis, apples and pears, bananas, watermelon, fruit vegetables, green vegetables, other vegetables) in order to evaluate their influence on OP metabolite concentrations.


Sociodemographic, environmental, and dietary characteristics of the study population are shown in Table 2. Nearly all the participants were Spaniards, 25% of them had finished university studies, around 80% worked during pregnancy, and nearly half of them belonged to the lowest social class. More than 60% of the women used domestic pesticides at home during pregnancy. Only 16 of the fathers and none of the mothers had an occupation related with the use of pesticides (agriculture or gardening).

Table 2 Sociodemographic, environmental, and dietary characteristics of the study population, INMA-Valencia, Spain, 2003–2006

Differences between the study population (n = 573) and the excluded women (those who arrived at the third trimester visit but for whom OP metabolite measurements are unavailable, n = 221) were assessed. We observed statistically significant differences according to the parity, social class, and season of sampling. Among the excluded women there was a higher proportion of parity > 1, lower social class, and urine samples taken during winter.

The frequencies of detection of OP metabolites were low, in fact only DMTP and TCPY were detected in more than 30% of the samples. The GM of the concentrations ranged from 0.02 μg/g of creatinine for IMPY and DEAMPY to 2.76 μg/g of creatinine for DEP. Regarding the sum variables, the GM of the concentrations ranged from 3.07 μg/g of creatinine for sumDEP to 16.29 μg/g of creatinine for sumDAP (Table 3).

Table 3 Urinary OP concentrations (creatinine adjusted and unadjusted and molar) in pregnant women from INMA-Valencia cohort (Spain, 2003–2006, n = 573)

The characteristics of the study population related to the OP metabolite concentrations can be found in Table 4. Smoking habit was significantly statistically associated with DMTP concentrations, women who smoked during pregnancy being the ones who presented lower levels. The pattern with the other compounds was similar but not significant. We observed that women who smoked during pregnancy consumed less fruit (mean = 2.24; standard deviation [sd] = 1.29 daily servings) than non-smoker women (mean = 2.58; sd = 1.62 daily servings, p-value Mann Whitney test =0.001). Smokers also consumed fewer vegetables than non-smokers (mean = 2.22; sd = 1.42 vs. mean = 2.24; sd = 1.29 daily servings), but the differences were not statistically significant (Additional file 1: Table S1).

Table 4 Sociodemographic, environmental, and dietary characteristics associated with the OP concentrations. INMA-Valencia cohort (Spain, 2003–2006)a

Women with university studies had the lowest concentrations of sumDAP. Season of urine sampling was associated with nearly all the compounds, samples taken during summer or fall being those which had the highest concentrations. Women who lived in an urban zone had higher OP metabolites, as did women who lived near fields or greenhouses sprayed with pesticides. Women who had a yard or garden with plants at home had the highest concentrations of DMTP, TCPy, and sumDMP, and women who applied outdoor pesticides had the highest concentrations of sumDEP. The intake of fruit and vegetables during the third trimester of pregnancy was positively associated with the OP metabolite concentrations, and the BMI before pregnancy was inversely associated with nearly all the compounds. We analyzed the differences in the intake of fruits and vegetables according to the maternal BMI and observed that women with BMI > 25 Kg/m2 consumed more vegetables (2.35 weekly servings for women with BMI = 25–30 Kg/m2 and 2.30 weekly servings for women with BMI ≥ 30 Kg/m2) than women with BMI < 25 Kg/m2 (2.19 weekly servings), but differences were not statistically significant (p-value Kruskal Wallis test = 0.416) (Additional file 1: Table S1).

Women whose partners had an occupation related with the use of pesticides had higher levels of TCPY (mean [sd] = 7.2 [15.3] μg/g of creatinine) than those whose partners did not work with pesticides (mean [sd] = 2.8 [8.1] μg/g of creatinine, p-value Mann Whitney test = 0.046).

The relationship between the intake of specific fruits and vegetables with the OP metabolite concentrations has been also evaluated (Fig. 1). The intake of green vegetables during pregnancy was significantly and positively associated with all the compounds, except for TCPy. Other groups of fruits and vegetables also associated with a considerable number of compounds (>3) were stone fruits, kiwis, apples and pears, and fruit vegetables. TCPy was the OP metabolite least associated with the intake of fruits and vegetables, and sumDAP was the OPs metabolite group associated with the intake of the highest number of different fruits and vegetables.

Fig. 1

Association between the intake of different types of fruits and vegetables during pregnancy and the OP concentrations. All models were adjusted by smoking habit, educational level, season of sampling, zone of residence, yard with plants at home, application of outdoor pesticides, residence near fields or greenhouses, fruit intake, vegetable intake, body mass index before pregnancy, and creatinine.


In this Spanish birth cohort study we assessed maternal exposure to OPs and evaluated the associated factors. The detection frequencies were low, DMTP and TCPY being the most frequently detected metabolites. The concentrations were positively associated with maternal intake of fruits and vegetables during pregnancy, especially the intake of green and fruit vegetables, stone fruits, kiwis, and apples and pears. Other maternal characteristics related to the OP metabolite concentrations were BMI before pregnancy and smoking habit during pregnancy. Women with lower BMI and those who did not smoke presented higher OP metabolite concentrations. Moreover, mothers who had a yard or garden with plants at home or lived in the urban area were also more exposed to OPs.

Other epidemiological studies have assessed the maternal exposure to OPs during pregnancy (Table 5). The first birth cohort studies where OP metabolites were analyzed in maternal urine were conducted in USA, one in New York City [8], and the other in an agricultural community in Salinas Valley, California [32, 33]. In both studies most of the metabolite concentrations were higher than in our study; DEP was the only one with concentrations a little higher in our study than in the USA cohorts. However, the detection frequency of this metabolite in our cohort was lower, probably due to a higher limit of detection. In our cohort the percentage of women using indoor pesticides (63%) was higher than in the New York cohort (46%), but urinary OP metabolite concentrations did not associate with this indoor use obtained by questionnaire. The possible reason could be related with differences in the time window of exposure reflected by both urinary metabolites and questionnaire. Whilst urinary biomarkers reflect acute or short-term exposure better, questionnaires may assess exposure during a period of time greater than the half-life of these substances. Another possible reason is that pyrethroids, but not OPs, were the main active substance found in the domestic insecticides [5], although some OPs, such as chlorpyrifos, diazinon, and fenitrothion were still allowed as domestic insecticides when urine samples were taken.

Table 5 OP concentrations and associated factors among pregnant women observed in previous studies

Ye et al. [34] compared urinary OP metabolites among pregnant women from Norway, Netherlands, and USA. The OP metabolite concentrations were found to be higher in women from Rotterdam (Netherlands) than women from Oslo (Norway) and USA (NHANES) [34]. Regarding the comparison with our study, the metabolite concentrations observed in our study population were lower than those in the women from Norway, Netherlands, and USA. TCPy, DMTP, DEP, and DETP concentrations in our cohort were much lower than those reported for pregnant women from Shanghai (China) [14], Cincinnati (USA) [21], Mexico DF [35], and Canada [22, 24]. Similar concentrations of DEP and TCPy to those of our study were found in the Caribbean Islands [36] and Puerto Rico [23] studies, respectively, although the detection frequencies in our study were lower.

Maternal educational level was associated with TCPy concentrations in women from New York [8], women with lower/middle studies being the ones who had lower concentrations in comparison to women with a higher educational level. A similar pattern was observed in women from Cincinnati [21] and from Canada (Mirec Study) [24]: those with a higher education presented higher urinary concentrations of sumDEP and sumDMP. Other factors associated with exposure to OPs in the Cincinnati study were race, marital status, employment, and vegetable and fruit intake. Lower sumDEP, sumDMP, and sumDAP concentrations were observed among black, unmarried, and unemployed women. Similarly to our study, the intakes of vegetables and fruits were positively associated with OP metabolite concentrations. OP exposure in pregnant women from Puerto Rico was studied in relationship with their sociodemographic characteristics, consumption of specific vegetables and fruits 48 h prior to urine collection, and home pest-related issues [23]. Overall, OP metabolite concentrations were inversely associated with maternal age and positively associated with some vegetable and fruit items, such as raisins, collards, peanuts, grape juice, and cherries.

We observed that maternal BMI before pregnancy was inversely and consistently associated with nearly all the metabolites. This association has also been described in the Mirec Study [24], where BMI before pregnancy was inversely associated with both sumDMP and sumDEP. The authors commented that a possible explanation could be the awareness on the importance of good nutrition among women with BMI < 25 Kg/m2. In our study, we did not find statistically significant differences in the intake of vegetables and fruits according to the pre-pregnancy BMI. These results suggest that there could be another factor, maybe metabolic, that might be promoting the association between BMI and OP metabolites.

Another maternal characteristic related to the OP exposure in our study population was the place of residence. Women who lived in the urban area were more exposed to OPs. This finding was expected since Valencia is a city with a long agricultural tradition and the urban area is widely surrounded by fields with an intensive use of agriculture. In addition, women who lived near fields sprayed with pesticides presented higher sumOP concentrations. In fact, this variable has been considered a useful alternative to biomonitoring and a good proxy of OP concentrations [37].

This is one of the largest birth cohorts to evaluate prenatal exposure to OPs and to study the associated factors, and, as far as we know, the first one conducted in Spain. The longitudinal design of our study will allow evaluation of postnatal exposure and the possible associated effects on child development.

A limitation in this study could be that a single spot measurement may not adequately characterize exposure over the whole pregnancy, given that OPs have a fast clearance from the body. In fact, results from the CHAMACOS cohort showed impairment in children’s cognitive performance at 7 years old [38] associated to an average of two DAP metabolite measures during pregnancy, but not when the one-point-in-time measures were evaluated. We collected urine samples from all the INMA participants over two different periods during pregnancy, the first and third trimesters, but due to limited resources to measure OP metabolites in all the urine samples we only have one measurement. The average of two OP metabolite measurements would have probably been a more adequate proxy of prenatal exposure. Additionally, there was a long time lag between collection of samples and analysis. The low concentrations observed in this population could be related with this limitation; in fact, Hoppin et al. [39] observed that urinary TCPY concentrations decreased with the time elapsed since collection at room temperature; however, our samples were frozen at −20 °C immediately after collection and this drop in OP metabolite concentrations is expected not be so important. Our FFQ covered average intakes during the 3 months before the sampling, when the half-life of these compounds is shorter. In any case, we observed several statistical associations with the intake of some fruits and vegetables.

Another important consideration is the limitation of measuring OP metabolites instead of parent compounds. The metabolism of parent compounds results in the production of active oxon forms, responsible for the developmental neurotoxicity [40]. These oxon forms are metabolically converted to less toxic compounds (DAPs) via hydrolysis. However, the hydrolysis of parent compounds present in food and the environment also generates DAP metabolites that are not toxic when consumed [41,42,43]. Therefore urinary non-specific metabolites may represent exposure to parent pesticides and to preformed derivatives from food and the environment, and the measurement of DAP metabolites could overestimate the real exposure to parent compounds. In fact, Yolton et al. [21] observed that higher urinary concentrations of DE metabolites were associated with improved attention in infants. This positive association seems to be related to the DAP metabolites present in vegetables and fruits that could be acting as a proxy of the nutrients and antioxidants present in diet. Despite this limitation, DAP metabolites have been extensively used in epidemiological studies because these metabolites are common to the majority of OP pesticides, and laboratory methods are not available for pesticide-specific metabolites. Consequently, the non-specific DAPs provide valuable information about cumulative exposure to the OP class [44].

Finally, we were able to provide an estimate of the geometric mean for all metabolites using censored regression. However, this estimation could be imprecise due to the distributional assumptions of these models and the high proportion of undetected values.


In conclusion, the OP metabolite detection frequencies and the concentrations observed in our study population were low, compared with previously published studies. For TCPy, the metabolite of chlorpyrifos, the concentrations were low despite the fact that it was still allowed as a domestic pesticide in Spain when the urine sampling was performed. Some maternal characteristics were associated with prenatal OP exposure. More in-depth knowledge of the determinants of the body burden of OPs during pregnancy would enhance the effectiveness of public health control and future recommendations for decreasing the levels of these compounds, thus reducing the risk for the health of pregnant women and the health and development of their children.



Body mass index before pregnancy




Diethyl phosphate


Diethyl thiophosphate


Dimethyl dithiophosphate


Dimethyl thiophosphate


Food frequency questionnaire


Geometric mean




INfancia y Medio Ambiente


Limit of determination


Organophosphate pesticides


Phenoxybenzoic acid






sumDEP + sumDMP










High resolution mass spectrometry


  1. 1.

    European Crop Protection. Industry Statistics 2001–2010. 2010. Accessed 12 May 2016.

  2. 2.

    Asociación de empresarios para la producción de las plantas (AEPLA). Report 2009. 2009. Accessed 15 Apr 2016.

  3. 3.

    Spanish Ministry of Agriculture F and E. Annual statistics about the pesticide use in the agriculture. 2013. Accessed 20 May 2016.

  4. 4.

    Ccanccapa A, Masia A, Andreu V, Pico Y. Spatio-temporal patterns of pesticide residues in the Turia and Jucar rivers (Spain). Sci Total Environ. 2016;540:200–10.

    CAS  Article  Google Scholar 

  5. 5.

    Moreno MJ, Melia LA, Oltra Moscardo MT, Jimenez PR. Current situation in Spain of aerosol insect sprays registered for household use by the environmental health authorities. Rev EspSalud Publica. 2003;77:383–91.

    Google Scholar 

  6. 6.

    Sociedad Española de Sanidad Ambiental. Study of the use of organophosphate insecticide Chlorpyrifos in urban spaces. Rev Salud Ambient. 2009;9(Supl 1):1–20.

    Google Scholar 

  7. 7.

    Whyatt RM, Barr DB, Camann DE, Kinney PL, Barr JR, Andrews HF, et al. Contemporary-use pesticides in personal air samples during pregnancy and blood samples at delivery among urban minority mothers and newborns. Environ Health Perspect. 2003;111:749–56.

    CAS  Article  Google Scholar 

  8. 8.

    Berkowitz GS, Obel J, Deych E, Lapinski R, Godbold J, Liu Z, et al. Exposure to indoor pesticides during pregnancy in a multiethnic, urban cohort. Environ Health Perspect. 2003;111:79–84.

    Article  Google Scholar 

  9. 9.

    Bradman A, Eskenazi B, Barr DB, Bravo R, Castorina R, Chevrier J, et al. Organophosphate urinary metabolite levels during pregnancy and after delivery in women living in an agricultural community. Environ Health Perspect. 2005;113:1802–7.

    CAS  Article  Google Scholar 

  10. 10.

    Roberts JR, Karr CJ, Council On Environmental Health. Pesticide exposure in children. Pediatrics. 2012;130:e1765–88.

    Article  Google Scholar 

  11. 11.

    Selevan SG, Kimmel CA, Mendola P. IDentifying critical windows of exposure for children’s health. EnvironHealth Perspect. 2000;108 Suppl 3:451–455.

  12. 12.

    Gonzalez-Alzaga B, Lacasana M, Guilar-Garduno C, Rodriguez-Barranco M, Ballester F, Rebagliato M, et al. A systematic review of neurodevelopmental effects of prenatal and postnatal organophosphate pesticide exposure. Toxicol Lett. 2014;230:104–21.

    CAS  Article  Google Scholar 

  13. 13.

    Eskenazi B, Harley K, Bradman A, Weltzien E, Jewell NP, Barr DB, et al. Association of in utero organophosphate pesticide exposure and fetal growth and length of gestation in an agricultural population. Environ Health Perspect. 2004;112:1116–24.

    CAS  Article  Google Scholar 

  14. 14.

    Wang P, Tian Y, Wang XJ, Gao Y, Shi R, Wang GQ, et al. Organophosphate pesticide exposure and perinatal outcomes in shanghai, China. Environ Int. 2012;42:100–4.

    CAS  Article  Google Scholar 

  15. 15.

    Rauch SA, Braun JM, Barr DB, Calafat AM, Khoury J, Montesano AM, et al. Associations of prenatal exposure to organophosphate pesticide metabolites with gestational age and birth weight. Environ Health Perspect. 2012;120:1055–60.

    CAS  Article  Google Scholar 

  16. 16.

    Whyatt RM, Rauh V, Barr DB, Camann DE, Andrews HF, Garfinkel R, et al. Prenatal insecticide exposures and birth weight and length among an urban minority cohort. Environ Health Perspect. 2004;112:1125–32.

    CAS  Article  Google Scholar 

  17. 17.

    Harari R, Julvez J, Murata K, Barr D, Bellinger DC, Debes F, et al. Neurobehavioral deficits and increased blood pressure in school-age children prenatally exposed to pesticides. Environ Health Perspect. 2010;118:890–6.

    CAS  Article  Google Scholar 

  18. 18.

    Raanan R, Balmes JR, Harley KG, Gunier RB, Magzamen S, Bradman A, et al. Decreased lung function in 7-year-old children with early-life organophosphate exposure. Thorax. 2016;71:148–53.

    Article  Google Scholar 

  19. 19.

    Bost-Legrand A, Warembourg C, Massart C, Chevrier C, Bonvallot N, Monfort C, et al. Prenatal exposure to persistent organic pollutants and organophosphate pesticides, and markers of glucose metabolism at birth. Environ Res. 2016;146:207–17.

    Article  Google Scholar 

  20. 20.

    Slotkin TA. Does early-life exposure to organophosphate insecticides lead to prediabetes and obesity? ReprodToxicol. 2011;31:297–301.

    CAS  Google Scholar 

  21. 21.

    Yolton K, Xu Y, Sucharew H, Succop P, Altaye M, Popelar A, et al. Impact of low-level gestational exposure to organophosphate pesticides on neurobehavior in early infancy: a prospective study. Environ Health. 2013;12:79.

    Article  Google Scholar 

  22. 22.

    Colapinto CK, Arbuckle TE, Dubois L, Fraser W. Tea consumption in pregnancy as a predictor of pesticide exposure and adverse birth outcomes: the MIREC study. Environ Res. 2015;142:77–83.

    CAS  Article  Google Scholar 

  23. 23.

    Lewis RC, Cantonwine DE, Del Toro LV, Calafat AM, Valentin-Blasini L, Davis MD, et al. Distribution and determinants of urinary biomarkers of exposure to organophosphate insecticides in Puerto Rican pregnant women. Sci Total Environ. 2015;512–513:337–44.

    Article  Google Scholar 

  24. 24.

    Sokoloff K, Fraser W, Arbuckle TE, Fisher M, Gaudreau E, LeBlanc A, et al. Determinants of urinary concentrations of dialkyl phosphates among pregnant women in Canada - results from the MIREC study. Environ Int. 2016;94:133–40.

    CAS  Article  Google Scholar 

  25. 25.

    Guxens M, Ballester F, Espada M, Fernandez MF, Grimalt JO, Ibarluzea J, et al. Cohort profile: the INMA--INfancia y Medio Ambiente--(environment and childhood) project. Int J Epidemiol. 2012;41:930–40.

    Article  Google Scholar 

  26. 26.

    Roca M, Leon N, Pastor A, Yusa V. Comprehensive analytical strategy for biomonitoring of pesticides in urine by liquid chromatography-orbitrap high resolution masss pectrometry. J Chromatogr A. 2014;1374:66–76.

    CAS  Article  Google Scholar 

  27. 27.

    Anastassiades M, Lehotay SJ, Stajnbaher D, Schenck FJ. Fast and easy multiresidue method employing acetonitrile extraction/partitioning and “dispersive solid-phase extraction” for the determination of pesticide residues in produce. J AOAC Int. 2003;86:412–31.

    CAS  Google Scholar 

  28. 28.

    Domingo-Salvany A, Regidor E, Alonso J, Varez-Dardet C. Proposal for a social class measure. Working Group of the Spanish Society of epidemiology and the Spanish Society of Family and Community Medicine. Aten Primaria. 2000;25:350–63.

    CAS  Article  Google Scholar 

  29. 29.

    Vioque J, Navarrete-Munoz EM, Gimenez-Monzo D, Garcia-de-la-Hera M, Granado F, Young IS, et al. Reproducibility and validity of a food frequency questionnaire among pregnant women in a Mediterranean area. Nutr J. 2013;12:26.

    CAS  Article  Google Scholar 

  30. 30.

    Newman MC, Dixon PM, Looney BB, Pinder JE. Estimating mean and variance for environmental samples with below detection limit observations. Water Resources Butlletin. 1989;25:905–16.

    CAS  Article  Google Scholar 

  31. 31.

    Lubin JH, Colt JS, Camann D, Davis S, Cerhan JR, Severson RK, et al. Epidemiologic evaluation of measurement data in the presence of detection limits. Environ Health Perspect. 2004;112:1691–6.

    CAS  Article  Google Scholar 

  32. 32.

    Castorina R, Bradman A, Fenster L, Barr DB, Bravo R, Vedar MG, et al. Comparison of current-use pesticide and other toxicant urinary metabolite levels among pregnant women in the CHAMACOS cohort and NHANES. Environ Health Perspect. 2010;118:856–63.

    CAS  Article  Google Scholar 

  33. 33.

    Castorina R, Bradman A, McKone TE, Barr DB, Harnly ME, Eskenazi B. Cumulative organophosphate pesticide exposure and risk assessment among pregnant women living in an agricultural community: a case study from the CHAMACOS cohort. Environ Health Perspect. 2003;111:1640–8.

    Article  Google Scholar 

  34. 34.

    Ye X, Pierik FH, Angerer J, Meltzer HM, Jaddoe VW, Tiemeier H, et al. Levels of metabolites of organophosphate pesticides, phthalates, and bisphenol a in pooled urine specimens from pregnant women participating in the Norwegian mother and child cohort study (MoBa). Int J Hyg Environ Health. 2009;212:481–91.

    CAS  Article  Google Scholar 

  35. 35.

    Fortenberry GZ, Meeker JD, Sanchez BN, Barr DB, Panuwet P, Bellinger D, et al. Urinary 3,5,6-trichloro-2-pyridinol (TCPY) in pregnant women from Mexico City: distribution, temporal variability, and relationship with child attention and hyperactivity. Int J Hyg Environ Health. 2014;217:405–12.

    CAS  Article  Google Scholar 

  36. 36.

    Forde MS, Robertson L, Laouan Sidi EA, Cote S, Gaudreau E, Drescher O, et al. Evaluation of exposure to organophosphate, carbamate, phenoxy acid, and chlorophenol pesticides in pregnant women from 10 Caribbean countries. Environ Sci Process Impacts. 2015;17:1661–71.

    CAS  Article  Google Scholar 

  37. 37.

    Gonzalez-Alzaga B, Hernandez AF, Rodriguez-Barranco M, Gomez I, Guilar-Garduno C, Lopez-Flores I, et al. Pre- and postnatal exposures to pesticides and neurodevelopmental effects in children living in agricultural communities from south-eastern Spain. Environ Int. 2015;85:229–37.

    CAS  Article  Google Scholar 

  38. 38.

    Bouchard MF, Chevrier J, Harley KG, Kogut K, Vedar M, Calderon N, et al. Prenatal exposure to organophosphate pesticides and IQ in 7-year-old children. Environ Health Perspect. 2011;119:1189–95.

    CAS  Article  Google Scholar 

  39. 39.

    Hoppin JA, Ulmer R, Calafat AM, Barr DB, Baker SV, Meltzer HM, et al. Impact of urine preservation methods and duration of storage on measured levels of environmental contaminants. J Expo Sci Environ Epidemiol. 2006;16:39–48.

    CAS  Article  Google Scholar 

  40. 40.

    Flaskos J. The developmental neurotoxicity of organophosphorus insecticides: a direct role for the oxon metabolites. Toxicol Lett. 2012;209:86–93.

    CAS  Article  Google Scholar 

  41. 41.

    Chen L, Zhao T, Pan C, Ross JH, Krieger RI. Preformed biomarkers including dialkylphosphates (DAPs) in produce may confound biomonitoring in pesticide exposure and risk assessment. J Agric Food Chem. 2012;60:9342–51.

    CAS  Article  Google Scholar 

  42. 42.

    Quirós-Alcalá L, Bradman A, Smith K, Weerasekera G, Odetokun M, Barr DB, et al. Organophosphorous pesticide breakdown products in house dust and children’s urine. J Expo Sci Environ Epidemiol. 2012;22:559–68.

    Article  Google Scholar 

  43. 43.

    Lu C, Bravo R, Caltabiano LM, Irish RM, Weerasekera G, Barr DB. The presence of dialkylphosphates in fresh fruit juices: implication for organophosphorus pesticide exposure and risk assessments. J Toxicol Environ Health Part A. 2005;68:209–27.

    CAS  Article  Google Scholar 

  44. 44.

    Bravo R, Caltabiano LM, Weerasekera G, Whitehead RD, Fernandez C, Needham LL, et al. Measurement of dialkyl phosphate metabolites of organophosphorus pesticides in human urine using lyophilization with gas chromatography-tandem mass spectrometry and isotope dilution quantification. J Expo Anal Environ Epidemiol. 2004;14:249–59.

    CAS  Article  Google Scholar 

Download references


The authors would particularly like to thank all the participants for their generous collaboration.


This study was supported by grants from the Instituto de Salud Carlos III [FIS-FEDER 13/1944, 13/2032, 14/0891, 14/1687, 16/1288 and Miguel Servet-FEDER CP15/0025]; Conselleria de Sanitat, Generalitat Valenciana (FISABIO UGP 15–230) and EU (FP7-ENV-2011 DENAMIC cod 282,957).

Availability of data and materials

The datasets generated and/or analyzed during the current study are not publicly available due to ethical concerns but are available from the corresponding author on reasonable request. All INMA questionnaires and protocols are available at

Authors’ contributions

Authors contributed to the article as follows: SL conceived the study, supervised the data collection and data analysis, and prepared the manuscript. MM conducted data analysis of the association of interest. MR and LLG participated in data collection. All authors provided critical revision of the manuscript and helped with data interpretation and manuscript preparation. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Informed consent was obtained from all participants and the study was approved by the Ethic Committees of the La Fe hospital (Valencia).

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information



Corresponding author

Correspondence to Sabrina Llop.

Additional file

Additional file 1: Table S1.

Intake of fruit and vegetables according to the maternal sociodemographic and life style characteristics. (DOC 49 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, 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. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Llop, S., Murcia, M., Iñiguez, C. et al. Distributions and determinants of urinary biomarkers of organophosphate pesticide exposure in a prospective Spanish birth cohort study. Environ Health 16, 46 (2017).

Download citation


  • Pesticides
  • Diet
  • Vegetables
  • Fruit
  • Agriculture
  • Development