Cognitive Performance and Exposure to Organophosphate Flame Retardants in Children: Evidence from a Cross-Sectional Analysis of Two European Mother–Child Cohorts
Authors: Rosolen Valentina , Giordani Elisa , Mariuz Marika , Parpinel Maria , Mustieles Vicente , Gilles Liese , Govarts Eva , Rodriguez Martin Laura , Baken Kirsten , Schoeters Greet , Sepai Ovnair , Sovcikova Eva , Fabelova Lucia , Kohoutek Jiři , Jensen Tina Kold , Covaci Adrian , Roggeman Maarten , Melymuk Lisa , Klánová Jana , Castano Argelia , Esteban López Marta , Barbone Fabio 
Published: October 2023

The knowledge of the effects of organophosphate flame retardants on children’s neurodevelopment is limited. The purpose of the present research is to evaluate the association between exposure to organophosphate flame retardants and children’s neurodevelopment in two European cohorts involved in the Human Biomonitoring Initiative Aligned Studies. The participants were school-aged children belonging to the Odense Child Cohort (Denmark) and the PCB cohort (Slovakia). In each cohort, the children’s neurodevelopment was assessed through the Full-Scale Intelligence Quotient score of the Wechsler Intelligence Scale for Children, using two different editions. The children’s urine samples, collected at one point in time, were analyzed for several metabolites of organophosphate flame retardants. The association between neurodevelopment and each organophosphate flame retardant metabolite was explored by applying separate multiple linear regressions based on the approach of MM-estimation in each cohort. In the Danish cohort, the mean ± standard deviation for the neurodevelopment score was 98 ± 12; the geometric mean (95% confidence interval (95% CI)) of bis(1,3-dichloro-2-propyl) phosphate (BDCIPP) standardized by creatinine (crt) was 0.52 µg/g crt (95% CI = 0.49; 0.60), while that of diphenyl phosphate (DPHP) standardized by crt was 1.44 µg/g crt (95% CI = 1.31; 1.58). The neurodevelopment score showed a small, negative, statistically imprecise trend with BDCIPP standardized by crt (β = −1.30; 95%CI = −2.72; 0.11; p-value = 0.07) and no clear association with DPHP standardized by crt (β = −0.98; 95%CI = −2.96; 0.99; p-value = 0.33). The neurodevelopment score showed a negative trend with BDCIPP (β = −1.42; 95% CI = −2.70; −0.06; p-value = 0.04) and no clear association with DPHP (β = −1.09; 95% CI = −2.87; 0.68; p-value = 0.23). In the Slovakian cohort, the mean ± standard deviation for the neurodevelopment score was 81 ± 15; the geometric mean of BDCIPP standardized by crt was 0.18 µg/g crt (95% CI = 0.16; 0.20), while that of DPHP standardized by crt was 2.24 µg/g crt (95% CI = 2.00; 3.52). The association of the neurodevelopment score with BDCIPP standardized by crt was −0.49 (95%CI = −1.85; 0.87; p-value = 0.48), and with DPHP standardized by crt it was −0.35 (95%CI = −1.90; 1.20; p-value = 0.66). No clear associations were observed between the neurodevelopment score and BDCIPP/DPHP concentrations that were not standardized by crt. No clear associations were observed with bis(1-chloro-2-propyl) phosphate (BCIPP) in either cohort, due to the low detection frequency of this compound. In conclusion, this study provides only limited evidence of an inverse association between neurodevelopment and exposure to BDCIPP and DPHP. The timing of exposure and effect modification of other organophosphate flame retardant metabolites and other substances should be the subject of further investigations that address this scientific hypothesis.