Next Article in Journal
Connections among Land Use, Water Quality, Biodiversity of Aquatic Invertebrates, and Fish Behavior in Amazon Rivers
Previous Article in Journal
Head-to-Head Study of Developmental Neurotoxicity and Resultant Phenotype in Rats: α-Hexabromocyclododecane versus Valproic Acid, a Recognized Model of Reference for Autism Spectrum Disorders
Previous Article in Special Issue
Residential Green and Blue Spaces and Type 2 Diabetes Mellitus: A Population-Based Health Study in China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Diurnal Variation in Biomarkers of Exposure to Endocrine-Disrupting Chemicals and Their Association with Oxidative Damage in Norwegian Adults: The EuroMix Study

1
Cyprus International Institute for Environmental and Public Health, Cyprus University of Technology, Limassol 3041, Cyprus
2
Norwegian Institute of Public Health (NIPH), P.O. Box 222 Skøyen, N-0213 Oslo, Norway
3
Center for Environmental and Respiratory Health Research (CERH), Faculty of Medicine, University of Oulu, Biocenter, FI-90014 Oulu, Finland
*
Author to whom correspondence should be addressed.
Toxics 2022, 10(4), 181; https://doi.org/10.3390/toxics10040181
Submission received: 6 March 2022 / Revised: 27 March 2022 / Accepted: 2 April 2022 / Published: 7 April 2022

Abstract

:
Much evidence on the adverse health effects of endocrine-disrupting chemicals (EDCs) has accumulated during recent decades. EDCs are commonly found in various foods and personal care products (PCP). Data documenting a diurnally varying EDC metabolism in humans is scarce. This study examined (i) the time-of-day effect on the diurnal magnitude and variance of urinary biomarkers of exposure to EDCs, and (ii) the association between EDC exposures and oxidative damage in a Norwegian adult subpopulation. This was a cross-sectional panel study using biobanked samples from the EuroMix project. During a typical weekday, participants were asked to collect all day’s urine voids and record dietary and PCP habitual uses in a diary. Collected time stamps of urine voids were classified into three distinct periods in the day (morning 6 a.m.–12 p.m., mid-day 12 p.m.–6 p.m., evening 6 p.m.–6 a.m.). Questionnaires regarding demographic characteristics, personal care product usage, and dietary habits were completed. Urinary levels of EDCs (phthalates, parabens, and bisphenols) were measured using mass spectrometry and adjusted for urinary volume using specific gravity. Urinary 4-hydroxynonenal (4HNE), a lipid peroxidation marker, was measured using an immunoassay kit. Linear mixed-effect models identified EDCs under the influence of a diurnal variation effect that was adjusted for dietary habits and PCP use and examined associations between EDC and 4HNE. p-values were FDR-adjusted. Most phthalates appeared to be diurnally varying with higher urinary levels towards the evening (q < 0.001) than those measured during mid-day; this strong diurnal variation effect was not present for parabens and bisphenols. Significant (q < 0.001) positive associations were observed between all phthalates, parabens, and bisphenols (except bisphenol S) and 4HNE. This study’s findings highlighted the diurnal variation of excretion for certain EDC, but not for others, in real-life conditions. The degree of EDC chronotoxicity in distinct diurnal windows of the day warrants further investigation with longitudinal human studies.

1. Introduction

Much scientific interest and discussion about the impact of endocrine-disrupting chemicals (EDCs) on humans and the environment have emerged during recent decades [1]. According to the Organization for Economic Co-operation and Development (OECD), “An endocrine disruptor is an exogenous substance or mixture that alters function(s) of the endocrine system and consequently causes adverse health effects in an intact organism, or its progeny, or (sub)populations” [2]. The economic costs of diseases and dysfunction associated with exposure to EDCs in the European Union (EU) have been estimated to be hundreds of billions of euros per year [3].
EDCs are ubiquitous in most environmental media, while humans may come in contact with them via all routes and various sources of exposure. Prominent examples of EDC classes are phthalates, bisphenols, and parabens, which hold specific characteristics of non-persistent chemicals, such as limited lipophilicity and relatively short half-lives of elimination [4,5]; indeed, the biological half-life of metabolites of phthalate/DINCH in urine was <12 h for all measured compounds [6,7,8]. BPA has a half-life of <6–7 h [9,10], while the half-lives of most parabens range between 13 h and 29 h [7,11].
Phthalates have been widely used as plasticizers, providing flexibility to rigid plastics. They are used in food-contact materials, toys, as emulsifying agents and solvents in cosmetics, and excipients in pharmaceuticals [12,13]. Phthalates have been associated with a series of non-communicable diseases, ranging from reproductive toxicity in both adults and adolescents across genders, impaired infant reproductive health, insulin resistance and type II diabetes, hormonal cancer, etc. [4,14,15,16]. Parabens are widely used as preservatives in cosmetics and pharmaceutical products [17]. Bisphenols are used primarily in the production of polycarbonate plastics and epoxy resins [18]. Parabens and bisphenols have been associated with respiratory diseases, atopic dermatitis, allergic rhinitis, reproductive dysfunction, or cancer [4,19,20,21,22].
The widespread usage of phthalates, parabens, and bisphenols leads to ubiquitous, constant, and potentially inevitable exposure in humans [4,23]. Due to the short half-life of the aforementioned EDCs, variations in EDCs’ body burden may depend on the source of EDC exposure and the time of exposure [5,8,24]. Therefore, the time of the sampling is important for detecting within- and between-day variations in human exposure. Interestingly, temporal within- and between-day variations in concentrations have been shown for metabolites of phthalates, bisphenols, and parabens [5,8,9,11,23,24,25,26]. One of these studies examined the impact of time-of-day expression of enzymatic/metabolic activity (CYP2E1) on the observed diurnal (within-day) variability of exposure to disinfection byproducts [24]. Non-persistent EDCs may be characterized by low intraclass correlation coefficient (ICC) estimates (<0.5), where single-spot urine data would likely lead to a higher probability of exposure misclassification [5,8,23]. Time-resolved, repeated measures in human studies would be needed to reliably estimate the magnitude and variance of exposure to non-persistent EDC [5,23].
Synchronization between the central circadian clock and the peripheral liver clock’s enzymatic activity may be key in controlling the hepatic enzyme-driven EDC metabolism patterns [24,27]. The interplay between EDC exposure patterns and their Phase I/II metabolism rates as influenced by the endogenous clock-dictated CYP450 enzymatic activity patterns remains much unexplored. These EDCs have been shown to induce oxidative stress in human studies, initiating lipid peroxidation in cellular membranes [28]. Lipid peroxidation may generate reactive-oxygen-species-induced tissue injury via the formation of lipid peroxidation products, such as 4HNE [29]. 4HNE has been involved in various pathologies, such as metabolic diseases, neurodegenerative diseases [30], cancers [31], alcoholic liver disease [32], and chronic obstructive pulmonary disease [33].
The biological hypothesis that the liver-clock-driven metabolism of EDC diurnally differs by chemical class-specific enzyme systems was tested in a repeated-measures adult study in Norway. The lipid peroxidation potential of these EDC exposures associated with dietary and personal-care-product use (PCP) was also considered. Thus, the main objectives of this study were to (i) investigate the influence of the time-of-day effects on the magnitude and variation of biomarkers of exposure to EDCs in Norwegian adults and (ii) evaluate the association between the biomarkers of exposure to EDCs with the biomarkers of lipid peroxidation (4-hydroxynonenal, 4HNE) measurements.

2. Materials and Methods

2.1. Data Collection

A detailed description of the EuroMix study, including the study’s inclusion/exclusion criteria, can be found in Husøy et al. [34]. In short, a total of 144 participants (44 men and 100 women) were recruited amongst employees of governmental institutes and authorities and universities in the counties of Oslo and Akershus in Norway between September 2016 and November 2017. The recruited participants might not be a representative sample of the overall Norwegian population, but this was not the aim of the biomonitoring study. A sample size of approximately 150 was considered sufficient as previous studies conducted by our team using comparable sample collection and study design had similar sample sizes [34,35,36]. Participants collected every urination during a 24-h period in separate containers and marked them with the time and date. The day after collection, the samples were measured and stored as pools of urine combining all the voids into three time periods/pools (morning time group 1 = 06:00–12:00, mid-day time group 2 = 12:00–18:00, evening time group 3 = 18:00–06:00 the next day) (Figure 1). Of the 144 participants, six did not submit urine for one of the time periods; these were excluded from the analyses.
A total of 144 participants (44 males aged 25–72 years and 100 females aged 24–72 years) completed a diary of their weighed food record, a cosmetic diary, a food frequency questionnaire (FFQ) (for the 24 h period), and a questionnaire on socio-demographic and lifestyle characteristics (Figure 1). The study was approved by the Regional Committees for Medical and Health Research Ethics (REK ID no 2015/1868) on the 19th of November, 2015, and all the participants provided their written informed consent.

2.2. Biomarkers

Data on biomarkers of exposure to EDCs based on the 3 urine time group pools were included in this analysis. Twelve phthalate metabolites (MEP, MiBP, MnBP, MBzP, MEHP, MEHHP, MEOHP, MECPP, oh-MiNP, oxo-MiNP, cx-MiNP, oh-MPHP) and two metabolites of DINCH (oh-MINCH and oxo-MINCH]) were determined in the three urinary time pools of the day, using liquid chromatography coupled to mass spectrometry (LC-MS-MS) as described in Sabaredzovic et al. [37]. These metabolites were from eight parent compounds, as shown in Table S1. The concentrations of DEHP, DiNP, and DINCH in urine were given as an estimated sum of their metabolites (total ug), after adjusting for molecular weights, as sumDEHP (MEHP, MEHHP, MEOHP, MECPP, MMCHP), sumDiNP (oh-MiNP, oxo-MiNP, cx-MiNP), and sumDINCH (oh-MINCH, oxo-MINCH). Four parabens (MEPA, ETPA, PRPA, and BUPA) and three environmental phenols (BPA, BPS, and BPF) were also included (units ug/L) and were determined using the LC-MS-MS method [38]. Quality control and quality assurance details of the measured biomarkers can be found in the original study by Husøy et al. [34]. Urinary 4-hydroxynonenal (4-HNE, ug/L) was determined using immunoassay following the manufacturer’s instructions (Wuhan Fine Biotech Co., Ltd., China). The inter- and intra-coefficients of variability for the 4-HNE assay were 3–4% and 5–7%, respectively.

2.3. Data Handling

2.3.1. Data Cleaning

The following datasets were obtained, including their timestamps in the day: Demographics, 24 h food diaries, 24 h PCP diaries, EDC, and 4HNE measurements. Although data on thermal paper handling was collected, it was not included in this analysis as the number of participants who reported handling thermal paper was small (24%). In order to carry out the analyses, the datasets of the food and PCP diaries were modified. The times recording each meal and PCP use were grouped into a three-level categorical variable to match the time groups for urine collection: Time group 1 (morning, 06:00–12:00), time group 2 (mid-day 12:00–18:00), and time group 3 (evening 18:00–06:00). Time group 2 (mid-day) was set as the reference level.
The 44 different food categories recorded in the study were reclassified into 15 group categories, as follows: Beverages, sweets, bread, fish, grain, dairy, egg, cakes, fruit, vegetables, butter and oil, potatoes, meat, cheese, and others. Each food group consumed was transformed into a binary categorical variable with 0 = ‘did not eat’ and 1 = ‘eat’. All PCP products were categorized into the following 26 new classes: Shower gel, shampoo, conditioner, deodorant, facial cleanser, facial moisturizer, body lotion, anti-wrinkle cream, sunscreen, mouth wash, toothpaste, perfume, lip gloss/stick/balm, foundation, intimate soap, hand cream, foot cream, hair styling product, hair treatment, eye make-up product, rouge and powder product, make-up remover product, shaving product, antibacterial product, ‘oils’ product, and hand soap. Certain products that could not be classified into any of the categories due to the ambiguity of the product description in the dataset were not used in our analysis; these products were used by less than 5% of the participants: Vaseline, wax, serum, lenses, condom, nail polish, paint (sens) on hand/arms, contact lenses, waxing, concealer pencil, nasal spray, physiological saltwater, carnival make-up, and floss. Each PCP product was defined as a binary variable with 0 = ‘did not use it’ and 1 = ‘used it’.

2.3.2. Dealing with Missing Values

We removed variables with > 80% of their values as missing values. As for the variables that had ≥20% but less than 80% of their values as missing values, we proceeded with imputation using Monte Carlo imputation [39].

2.3.3. Data below Limit of Detection

Outcome measurements below the limit of detection (LOD) were substituted with LOD/2. Outcomes with >60% but <90% values below detection were used as binary variables (above LOD vs. below LOD). Outcomes with values >90% below LOD were not used in the final analysis.

2.3.4. Data Adjustment

To correct for urinary dilution, specific gravity (SG) was measured in all urine samples. SG-adjusted concentrations were used in the statistical analyses. Unlike creatinine-adjusted values, SG is less affected by factors such as diet, season, age, muscle mass, and gender [40]. In essence, all measurement values for the biomarkers of exposure to EDCs and the biomarker of effect (4HNE) were SG-adjusted using the Formula (1):
CHEMSG_Adj (µg/L) = CHEMi ((SGm − 1)/(SGi − 1))
where CHEMSG_Adj is the SG-standardized analyte concentration, CHEMi is the observed analyte concentration, SGi is the specific gravity of the urine sample, and SGm is the median specific gravity for the study population.

2.4. Statistical Analysis

Descriptive statistics (percentiles, min, and max) were calculated for the biomarkers of exposure to EDCs and 4HNE. Boxplots were constructed to visualize the time-of-day effect on EDCs and 4HNE.
Linear mixed-effect (LME) models were constructed for each biomarker of exposure to EDCs (log-transformed and SG-adjusted) accounting for the repeated measures within subjects. Fixed-effect variables were the time groups (mid-day time group 2 was the reference) along with the covariates of sex, age, and BMI. Select food groups and PCP types (transformed and imputed as binary variables) earlier found to be associated with the EDC levels in the original project analysis [34] were also included in the regressions: (i) For phthalates, meat, bread, beverages, butter, oil and PCPs of shower gel, hand cream, toothpaste, antiwrinkle cream, shaving products; and (ii) for both parabens and bisphenols, meat, bread, beverages, butter and oil, and lip products (BPA only). Participant id was included in the LME as the random-effect variable. Smoking status, along with educational level, was not included in the model, as the sample population was considered highly educated, with nearly 80% of the participants having at least a university degree and >84% of the participants not being active smokers. A correlation coefficient matrix was constructed to observe correlations between the studied EDCs, including 4HNE.
To explore the association between the biomarkers of exposure to EDCs and 4HNE, the second set of regression models was constructed using 4HNE (log-transformed and SG-adjusted) as the response variable. Each biomarker of exposure to EDCs was used as a fixed-effect variable, including covariates, such as sex, age, and BMI. The participants were considered random effects in the model.
In both sets of models, multiple testing correction was performed using the Benjamini–Hochberg (false discovery rate, FDR) method considering all regression models. A q-value < 0.05 indicates the statistical significance of an association after controlling the FDR at 5%. All analyses were performed in R (v.4.0.5) with RStudio (1.4.1717).
A sensitivity analysis was performed for the first set of regression models, using the frequency of consumption (in each time group of the day) for select food groups as a continuous function; associations and the regression coefficients remained essentially the same.

3. Results

3.1. Population Characteristics

The EuroMix study population (n = 144) had a mean age of 43 years with normal BMI, being mostly females, highly educated, with nearly 80% of the study participants having a university degree. There were no regular smokers, with 87% either having never smoked before or having quit smoking and the remaining 13% being active but only occasional smokers. The mean weights of the males and the females were 82 kg (standard deviation, SD = 8.5) and 65 kg (SD = 8.9), respectively (Table 1). Detailed descriptive information on diet from weighed food records and FFQ, as well as PCP use by the study population, was previously reported [34].

3.2. Biomarkers of Exposure/Effect

3.2.1. Urinary Biomarkers of Exposure to Phenols, Phthalates, and Parabens

Only 4% of the urine samples had BPA levels <LOD, while 71% and 96% of the samples had levels < LOD for BPS and BPF, respectively. There were varying detection levels for the parabens; all urine samples had MEPA concentrations above the LOD, while 1%, 35%, and 50% of the samples were below LOD for ETPA, PRPA, and BUPA, respectively (Table S2). The non-quantified rates (samples below LOD) for the metabolites of phthalates/DINCH ranged between 0% and 24%; all values were above LOD for 10 phthalates—MeP, MiBP, MnBP, MBzP, MEHHP, MEOHP, MECPP, oh_MiNP, oxo_MiNP, and cx_MiNP. The <LOD percentages for MEHP, MMCHP, oh-MINCH, oxo-MINCH, and oh-MPHP were 12%, 3%, 11%, 4%, and 24%, respectively (Table S1).
Of all the EDC biomarkers, the highest concentrations were observed in MEPA (morning, 18519 ug/L), ETPA (night, 2038 ug/L), and MEP (mid-day, 1990 ug). More details on the percentiles of the urinary EDCs by the time of day and overall are presented in Table 2. The concentrations of all the metabolites of phthalates/DINCH in the urine samples followed a diurnal variation pattern with the lowest concentrations observed during mid-day, with levels being ordered by evening > morning > mid-day (Figure S2). Urinary concentrations were higher in males than females for all phthalates except for MEP (Figure S3). Urinary BPA levels were higher in the morning, with a median value of 1.39 ug/L (IQR:0.62, 2.44), than at night (1.17 ug/L; IQR:0.68, 1.97) (Table 2, Figure S2). The diurnal concentration pattern for BPS could not be clarified as the lower concentration quartiles were below the LOD (Table 2). Concentrations of BPA in males were slightly higher than in females (Figure S3). The diurnal profile for MEPA followed the order of evening (7.32 ug/L; IQR: 3.2, 27.7) > morning (6.31 ug/L; IQR: 3.4, 37.2) > mid-day (5.85 ug/L; IQR: 2.95, 24.3). The diurnal profile of urinary PRPA and ETPA followed the order of morning (0.25 ug/L; IQR: 0.02, 2.8) > night (0.17 ug/L; IQR: 0.03, 3.40) > mid-day (0.10 ug/L; IQR: 0.02, 1.16), and morning (1.24 ug/L; IQR: 0.52, 3.43) > mid-day (1.08 ug/L; IQR: 0.40,3.19) > night (0.84 ug/L; IQR: 0.34, 3.26), respectively. The concentrations of all parabens were higher in females than in males (Figure S3).
The magnitude of urinary 4HNE levels by time group in the 24 h period did not diurnally vary (p > 0.05), because of its anticipated slower rate of formation at the expense of the faster EDC metabolism. However, higher levels of 4HNE were observed among males than in females (Figure S3).

3.2.2. Association between Time of the Day and the Biomarkers of EDC Exposure and Effect

The LME models for most phthalates/DINCH metabolites showed significant associations between urinary EDC concentrations and the time of the day (Figure 2, Table S4i–xiii). All models were adjusted for age, sex, BMI, dietary habits, and PCP use. Statistically significant positive associations were found between time (morning) and MiBP (β = 0.32, 95% CI: 0.15, 0.49, q < 0.001), MnBP (β = 0.38, 95% CI: 0.21, 0.55, q < 0.001), MBzP (β = 0.32, 95% CI: 0.14, 0.50, q = 0.006), MEP (β = 0.39, 95% CI: 0.16, 0.61, q = 0.006), and SumDEHP (β = 0.30, 95% CI: 0.11, 0.48, q = 0.011), but not with SumDiNP, SumDINCH, and oh-MPHP. Statistically significant associations (q < 0.001) were also found between time (night) and all phthalate metabolites (Table S4). There were statistically significant differences in the concentration of some phthalates (MiBP, MnBP, MBzP, SumDEHP, SumDiNP, and SumDINCH) by sex (Table S4). However, only MiBP, MnBP, SumDEHP, and SumDINCH remained significantly different (q < 0.05) after applying the Benjamini–Hochberg correction. The effect of time on the concentration of phthalates was unaffected by BMI (Table S4).
In the case of parabens, a statistically significant positive association was only found for PRPA and time (morning) (β = 0.57, 95% CI: 0.27, 0.86, q < 0.001) as well as time (night) (β = 0.50, 95% CI: 0.21, 0.86, q = 0.006). For the rest of the parabens, including bisphenols, there were no statistically significant associations between the urinary concentration levels of MEPA, ETPA, BUPA, and BPA and time of urination groups of the day; BPS was significantly associated with time (morning) (β = 0.41, 95% CI: 0.21, 0.81, q = 0.047) but not with time (night) (Figure 3, Table S4). Statistically significant differences by sex were observed in the concentrations of PRPA and MEPA (q < 0.001 and q = 0.043, respectively). There was little differentiation in the 4HNE levels among the three time groups of the day (Table S4).

3.2.3. Association between Biomarkers of Exposure to EDCs and the Biomarker of Oxidative Damage

An LME model was employed to test the hypothesis that the diurnally varying biomarkers of exposure would be associated with the biomarker of enzymatic activity/oxidative damage (4HNE). Statistically significant associations were found between 4HNE and all metabolites of phthalates/DINCH, phenols, and parabens, except for BPS. Significant associations remained for all biomarkers of exposure to EDCs after applying the Benjamini–Hochberg correction (q value < 0.001) (Table 3). The regression models on the associations between biomarkers of exposure to EDCs and the biomarker of effect (4HNE) can be found in Table S5.
Sex significantly impacted the association between 4HNE and the biomarkers of exposure to phthalates, phenols, and parabens; the magnitude of 4HNE associations was consistently higher in males for all EDCs tested (q-value < 0.001) (Table S5). Age and BMI had no impact on the association between the biomarkers of exposure and the biomarker of effect.

4. Discussion

We found that all urinary biomarkers of exposure to phthalates/DINCH followed a diurnal pattern, peaking in the evening and continuing to be high in the morning, and then gradually declining during the day with the lowest concentrations measured at mid-day. This consistent pattern was not so obvious for parabens or bisphenols. Earlier research by Gängler et al. [23] looked at the effect of the time of day on the variability of urinary biomarkers of exposure to disinfectant by-products. They observed that the levels of urinary trihalomethanes (THMs) during the day followed a diurnal pattern; higher levels of THM were measured in urine samples collected during the evening and lower levels in those collected during mid-day before and after engaging in household cleaning activities, respectively [23].
It was also speculated that the chrono-varying diurnal trend for phthalates could be associated with the variance in the day liver-clock-driven enzymatic expression of CYP450. Cytochromes P450 (CYPs) are a superfamily of enzymes that play a leading role in drug metabolism and detoxification, contributing to the metabolism of up to 75% of drugs. Biological pathways in which CYPs are involved, including drug and xenobiotic metabolism, have been shown to display circadian variation [41,42]. Several previous studies indicated that these CYPs are not directly regulated by the CLOCK/BMAL1 heterodimer, but indirectly through clock-regulated genes, such as PARbZip [43] and nuclear receptors [44]. While animal studies have shown circadian-based hepatotoxic effects of CYP450 activity as a result of xenobiotics treatment, only three human studies [23,45,46] have previously shown the influence of circadian/diurnal variation on the metabolism and possible chronotoxicity of environmental chemicals, including EDC. The diurnal variation in the magnitude of the EDC metabolites was evident in our study, but we are unsure whether this would be exclusively ascribed to the diurnal time-dependent expression of such liver CYPs enzyme systems. The possibility that food intake and its timestamp in the day may not have been sufficiently adjusted for remains, and it could be partially held responsible for some of the EDCs’ diurnal differences observed. It is also possible that the lack of considering the previous day’s food intake habits and quantities would have influenced the observed diurnal trends on the biomarkers of exposure to EDCs. Further, we did not take into account the EDC concentrations in the studied food items, which might lead to residual confounding. The sensitivity-based regression models using the number of food items consumed in each time group of the day did not alter the direction or strength of the diurnal time effect on the biomarkers of exposure to EDC after adjusting for the possible confounding effect of dietary intake. Select food items consumed in each time group of the day were chosen for inclusion in the regression models in two ways: (i) Added as a yes/no variable and (ii) added as a continuous variable in grams consumed per time group. In both cases, the food intake did not influence the effect size of the diurnal effects on urinary EDC.
Interestingly, this diurnal variation in the concentration of EDC metabolites was not consistently observed for parabens and bisphenols. Except for paraben PRPA, the rest of the parabens (MEPA, ETPA, BUPA) and bisphenols (BPA and BPS) did not follow the same diurnal variation as seen in phthalates. A notable difference between the phthalates and the parabens is that, on average, the urinary half-lives of excretion for parabens are longer (ranging between 13 and 29 h) [5,9] when compared with those of phthalates (<12 h) [4,5,6]. This difference may be important for the monitoring of diurnal (24-h) variation effects on the metabolism of short-lived chemicals. Studies also suggest that parabens are metabolized by carboxylesterases (CES) [47,48,49], with carboxylesterase1 (CES1) being the predominant CES enzyme in the liver, responsible for the metabolism of short-chain parabens (e.g., ETPA, PRPA, BUPA, and MEPA), while carboxylesterase2 (CES2) is more prevalent in skin tissue [49]. As such, the route of exposure plays an important role in the metabolism of parabens, as the level of absorption following oral intake has been shown to be higher than the levels following dermal uptake [50]. Findings on the exposure pattern of the participants to parabens based on the exposure to food or PCPs that were significantly correlated with the observed concentrations have been previously reported [34]. Hepatic CES enzymes have been shown to exhibit circadian variation in animal studies with higher levels in the dark phase than in the light phase [42,51,52]. However, Zhang et al. [43] clearly illustrated that although variations followed this pattern of diurnal expression, fluctuations in the expression of CES transcripts were low, with levels of hepatic CES isoenzymes during the dark phase being only slightly higher than levels during the light phase with additional minor troughs during the dark phase. In addition, the inhibitory influence of parabens on CYPs enzyme systems is also well known. There is some evidence that parabens may have an inhibitory effect on CYP450 enzymatic activity [53], thereby altering the pharmacokinetics of co-occurring EDC, such as BPA [54]. However, this panel study was not designed to study the biological mechanism of the observed diurnal variation patterns in the EDC metabolism as a result of the corresponding hepatic enzyme systems for each EDC class. There is a need for further studies to establish the potential impact of co-exposure on the metabolic/pharmacokinetic fate of individual EDCs.
To the best of our knowledge, this seems to be the first human study to report significant associations between 4HNE and EDCs (bisphenols, phthalates, and parabens). Most of the EDCs in this study showed associations with the marker of lipid oxidative damage (4HNE). 4HNE is a well-known biomarker of oxidative stress-induced lipid peroxidation both in animal models [55,56] and in human studies [57], and it has been associated with the pathogenesis and onset of cancer [58], arthritis [59], cardiovascular diseases, and diabetes [60]. The association between EDCs and oxidative damage observed in this study is supported by findings in earlier studies [61]. Urinary BPA was positively associated with 8-hydroxy-2′-deoxyguanosine (8-OHdG), 8-isoprostane, and 4-hydroxy-2-nonenal-mercapturic acid (HNE-MA), while BPF was associated with HNE-MA and 8-isoprostane in a Chinese cross-sectional study [62,63]. Similar associations between bisphenols and oxidative damage have been found in pregnant women and fetuses/newborns [64,65], or children and adolescents [66,67]. DEHP exposure in adolescents affected lipid metabolism by regulating the expression of lipid metabolism-related genes [68]. Urinary parabens have also been associated with oxidative stress markers in human studies [28].
In addition to chemical exposures, the oxidative stress imbalance exceeding the antioxidant defense mechanisms may also be associated with a suite of lifestyle and behavioral risk factors, such as alcohol or smoking [69]. A positive association was found earlier between smoking/tobacco exposure [70,71,72] and alcohol [73,74] consumption with increased levels of oxidative damage (4HNE). Moreover, a significant association was found between BMI and oxidative damage, with studies reporting a positive effect of obesity on increased 4HNE [73,75]. However, the majority of our study population was composed of non-active smokers and non-obese, thus their lifestyle characteristics were anticipated to have little influence on the variation of the observed 4HNE levels.
4-HNE-associated oxidative stress may also be accompanied by changes in enzyme activity patterns. For example, elevated formation of 4HNE was associated with the expression of GSTA4–4, CYP2E1, and Hsp70 proteins in the mitochondria of PC12 cells [76]. Based on a study with CYP2E1-deficient mice, 4HNE was elevated in mice exposed to environmental chemicals when compared with the CYP2E1-knock out mice [77]. Disinfection by-product chemical toxicity was mediated by its free-radical metabolism by CYP2E1, an isoform of cytochrome P450, forming 4HNE adducts [77,78]. Increased CYP2E1 activity led to significantly higher 4-HNE levels across cell types in the liver [79]. In biological systems where the 4HNE formation rate seems to be much slower than the EDC excretion rates (e.g., short-lived EDC with half-lives of elimination in the order of a few hours), 4HNE may not be used as a chrono-impacted marker in such a context, as also suggested in our study.
This study had some strengths and limitations. The use of repeatedly collected urine samples during the day, instead of a single spot urine void in the day, in analyzing exposure patterns and concentrations of the sample population to environmental compounds was a strength, as it allowed for observing the chrono-varying diurnal EDC biomarker variation and reduced the probability of exposure misclassification. Another strength was the simultaneous analysis of a suite of different EDC classes and various chemicals within each EDC class. However, there were also a few limitations. The lack of controlling for the time stamp of each food or PCP usage activity combined with varying half-lives of elimination for the different EDCs monitored made it difficult to match food/PCP exposures with the urinary EDC metabolites. An earlier analysis of the same study demonstrated that the dietary estimated exposure correlated better with the evening urinary phthalate metabolite levels [80]. The study population was not representative of EDC exposures of the Norwegian population based on demographics and dietary habits. Moreover, previous-day exposures, including dietary and PCP usage habits, were not recorded. The urination time stamps were matched with the group time stamps of eating and the use of PCP without having designed the study to consider the time needed for absorption and metabolism of the compounds before urinary excretion.
More studies are needed to effectively control for the time of food intake in relation to the monitored diurnal variation in the biomarkers of exposure to non-persistent EDCs (with relatively short half-lives of elimination) and their downstream biological perturbations. It becomes obvious that the totality of environmental agents, such as chemicals and diet, including behaviors and lifestyles throughout one’s lifetime or in critical life stages, including the endogenous response, shall be characterized in detail if one were to further advance our understanding of the interactions occurring at the interface of the external and internal exposome domains [27]. This may be achieved with the methodological framework and tools of the human exposome concept that, by definition, incorporates the global dimensions of chronos and location [27,81,82].

5. Conclusions

Evidence obtained in this study suggests that the concentrations of some urinary EDCs metabolites followed a diurnal pattern (phthalates), while this was less consistent for the studied parabens and bisphenols. Diurnal (24-h) variation in urinary biomarkers of exposure may be particularly important for short-lived chemicals. The biological mechanism responsible for this interesting phenomenon is not yet clear; the influence of metabolism and associated time-dependent enzymatic activity patterns is possible, but more data are needed. The impact of exposure to EDCs, even at low concentrations, was shown to have a significant impact on lipid-based oxidative damage. The timing of human exposure to certain EDCs and the associated (chrono)toxicity window may be worth studying in more detail. The human exposome concept and its critical windows of susceptibility are key methodological concepts that would help in better designing such human chronotoxicity studies. To this extent, exposome studies with a longitudinal design are needed to further elucidate the human chronotoxicity potential of EDCs and the associated biological mechanisms and pathways.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/toxics10040181/s1, Figure S1: Correlation matrix plot of the SG-adjusted, logarithm-transformed biomarkers of exposure to EDCs and the biomarker of effect or lipid peroxidation, 4HNE. From light blue to deep blue means positive correlation; from white to dark red means no or negative correlation; Figure S2: Box-plots for specific gravity-adjusted, logarithm-transformed values of EDCs and 4HNE per time group measured in urine samples; Figure S3: Visualization of the association between SG-adjusted, log-transformed EDCs and 4HNE with time group based on sex; Table S1: Phthalate metabolites determined in 24-hour urine of study day; Table S2: Environmental phenol and paraben substances determined in 24-hour urine of study day; Table S3: Median concentrations of the specific gravity-adjusted biomarkers of exposure to EDC and 4HNE by the time of urine sampling; Table S4 (i)–(xv): Linear mixed-effect models of specific gravity-adjusted biomarker of exposure & effect as a function of time of urine sampling in the day adjusted for dietary habits and personal care product use (reference: noon); Table S5 (i)–(xiv): Linear mixed-effect models for the association between the specific gravity-adjusted log(4HNE) and logarithm-transformed biomarker of exposure.

Author Contributions

Conceptualization, K.C.M.; Data curation, S.O.A., C.X. and K.C.M.; Formal analysis, S.O.A. and C.X.; Funding acquisition, A.K.S., T.H. and H.D.; Investigation, S.O.A., C.X., A.K.S., B.H., T.H., H.D. and K.C.M.; Methodology, K.C.M., A.K.S. and B.H.; Project administration, K.C.M., H.D., T.H. and B.H.; Resources, K.C.M., A.K.S., T.H. and H.D.; Supervision, K.C.M.; Validation, A.K.S., B.H., T.H. and K.C.M.; Writing—original draft, S.O.A. and K.C.M.; Writing—review & editing, S.O.A., C.X., A.K.S., B.H., T.H., H.D. and K.C.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the EXPOSOWORK project OPPORTUNITY/0916/MSCA/0077, co-financed by the European Regional Development Fund and the Republic of Cyprus through the Research and Innovation Foundation; EXPOSOWORK proposal received the Seal of Excellence by the European Commission as submitted in the Marie Skłodowska-Curie action H2020-MSCA-IF-2018.

Institutional Review Board Statement

The study was approved by the Regional Committees for Medical and Health Research Ethics (REK ID no 2015/1868) on the 19 November 2015.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

All data are available upon request to the corresponding author via email ([email protected]).

Acknowledgments

We thank all members of the Water and Health Laboratory at CII, Pantelis Charisiadis, Thibaut Delplancke, Corina Konstantinou, and Shamah Agboola for their critical contributions to the study execution. We also thank Eleni Zoumpoulia Papadopoulou for her critical insights into the study design of this work.

Conflicts of Interest

Authors declare no conflict of interest.

Abbreviations

4HNE, 4-Hydroxynonenal; BBzP, butyl benzyl phthalate; BMI, body mass index; BPA, bisphenol A; BPF, bisphenol F; BPs, bisphenols; BPS, bisphenol S; BUPA, butyl paraben; CES, carboxylesterases; CYPs, cytochromes P450; cx-MiNP, mono-4-methyl-7- carboxyoctyl phthalate; DEHP, diethyl phthalate; DEP, diethyl phthalate; DiBP, di-iso-butyl phthalate; DINCH, 1,2-cyclohexane dicarboxylic acid diisononyl ester; DiNP, di-iso-nonyl phthalate; DnBP, di-n-butyl phthalate; DPHP, di(2-propyl heptyl) phthalate; EDC, endocrine-disrupting chemical; ETPA, ethyl paraben; FDR, false discovery rate; FFQ, food frequency questionnaire; HNE-MA, 4-hydroxy-2-nonenal-mercapturic acid; ICC, intraclass correlation coefficient; LME, linear mixed-effect; LOD, limit of detection; LOQ, limit of quantification; MBzP, mono benzyl phthalate; MECPP, Mono-2-ethyl 5-carboxypentyl phthalate; MEHHP, mono-2-ethyl-5-hydroxyhexyl phthalate; MEHP, mono-2-ethylhexyl phthalate; MEOHP, mono-2-ethyl-5-oxohexyl phthalate; MEP, monoethyl phthalate; MiBP, mono-isobutyl phthalate; MMCHP, mono-2-carboxymethyl hexyl phthalate; MnBP, mono-n-butyl phthalate; MEP, methylparaben; OHdG, 8-hydroxy-2′-deoxyguanosine, 8-; OH-EtP, ethyl-protocatechuic acid; OH-MeP, methyl-protocatechuic acid; oh-MINCH, (2-(((Hydroxy-4-methyloctyl)oxy)carbonyl)-cyclohexane carboxylic acid; oh-MiNP, mono-4-methyl-7-hydroxyoctyl phthalate; oh-MPHP, 6-hydroxy monopropylheptyl phthalate; oxo-MINCH, 2-(((4-Methyl-7-oxyooctyl)oxy)carbonyl)-cyclohexanecarboxylic acid; oxo-MiNP, mono-4-methyl-7-oxooctyl phthalate; PBs, parabens; PCP, personal care product; PVC, Polyvinyl chloride; SCN, suprachiasmatic nucleus; SD, standard deviation; SG, specific gravity; THM, trihalomethanes.

References

  1. Gore, A.C.; Chappell, V.A.; Fenton, S.E.; Flaws, J.A.; Nadal, A.; Prins, G.S.; Toppari, J.; Zoeller, R.T. EDC-2: The Endocrine Society’s Second Scientific Statement on Endocrine-Disrupting Chemicals. Endocr. Rev. 2015, 36, E1–E150. [Google Scholar] [CrossRef] [PubMed]
  2. OECD. Revised Guidance Document 150 on Standardised Test Guidelines for Evaluating Chemicals for Endocrine Disruption; OECD Series on Testing and Assessment, No. 150; OECD Publishing: Paris, France, 2018. [Google Scholar] [CrossRef]
  3. Trasande, L.; Zoeller, R.T.; Hass, U.; Kortenkamp, A.; Grandjean, P.; Myers, J.P.; DiGangi, J.; Bellanger, M.; Hauser, R.; Legler, J.; et al. Estimating Burden and Disease Costs of Exposure to Endocrine-Disrupting Chemicals in the European Union. J. Clin. Endocrinol. Metab. 2015, 100, 1245–1255. [Google Scholar] [CrossRef] [PubMed]
  4. Diamanti-Kandarakis, E.; Bourguignon, J.-P.; Giudice, L.C.; Hauser, R.; Prins, G.S.; Soto, A.M.; Zoeller, R.T.; Gore, A.C. Endocrine-Disrupting Chemicals: An Endocrine Society Scientific Statement. Endocr. Rev. 2009, 30, 293–342. [Google Scholar] [CrossRef] [PubMed]
  5. LaKind, J.; Idri, F.; Naiman, D.Q.; Verner, M.-A. Biomonitoring and Nonpersistent Chemicals—Understanding and Addressing Variability and Exposure Misclassification. Curr. Environ. Health Rep. 2019, 6, 16–21. [Google Scholar] [CrossRef] [PubMed]
  6. Frederiksen, H.; Søgaard Nielsen, J.K.; Aarøe Mørck, T.; Hansen, P.W.; Jensen, J.F.; Nielsen, O.; Andersson, A.-M.; Knudsen, L.E. Urinary excretion of phthalate metabolites, phenols and parabens in rural and urban Danish mother–child pairs. Int. J. Hyg. Environ. Health 2013, 216, 772–783. [Google Scholar] [CrossRef] [PubMed]
  7. Larsson, K.; Björklund, K.L.; Palm, B.; Wennberg, M.; Kaj, L.; Lindh, C.H.; Jönsson, B.A.; Berglund, M. Exposure determinants of phthalates, parabens, bisphenol A and triclosan in Swedish mothers and their children. Environ. Int. 2014, 73, 323–333. [Google Scholar] [CrossRef] [Green Version]
  8. Koch, H.M.; Bolt, H.M.; Preuss, R.; Angerer, J. New metabolites of di(2-ethylhexyl)phthalate (DEHP) in human urine and serum after single oral doses of deuterium-labelled DEHP. Arch. Toxicol. 2005, 79, 367–376. [Google Scholar] [CrossRef]
  9. Morgan, M.K.; Nash, M.; Barr, D.B.; Starr, J.M.; Clifton, M.S.; Sobus, J. Distribution, variability, and predictors of urinary bisphenol A levels in 50 North Carolina adults over a six-week monitoring period. Environ. Int. 2017, 112, 85–99. [Google Scholar] [CrossRef]
  10. Völkel, W.; Colnot, T.; Csanády, G.A.; Filser, J.G.; Dekant, W. Metabolism and Kinetics of Bisphenol A in Humans at Low Doses Following Oral Administration. Chem. Res. Toxicol. 2002, 15, 1281–1287. [Google Scholar] [CrossRef]
  11. Ye, X.; Wong, L.-Y.; Bishop, A.M.; Calafat, A.M. Variability of Urinary Concentrations of Bisphenol A in Spot Samples, First Morning Voids, and 24-Hour Collections. Environ. Health Perspect. 2011, 119, 983–988. [Google Scholar] [CrossRef]
  12. Kwak, E.S.; Just, A.; Whyatt, R.; Miller, R.L. Phthalates, pesticides, and bisphenol-A exposure and the development of nonoccupational asthma and allergies: How valid are the links? Open Allergy J. 2009, 2, 45–50. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Crinnion, W.J. The CDC fourth national report on human exposure to environmental chemicals: What it tells us about our toxic burden and how it assists environmental medicine physicians. Altern. Med. Rev. 2010, 15, 101–109. [Google Scholar] [PubMed]
  14. Benjamin, S.; Masai, E.; Kamimura, N.; Takahashi, K.; Anderson, R.C.; Faisal, P.A. Phthalates impact human health: Epidemiological evidences and plausible mechanism of action. J. Hazard. Mater. 2017, 340, 360–383. [Google Scholar] [CrossRef] [PubMed]
  15. Lottrup, G.; Andersson, A.-M.; Leffers, H.; Mortensen, G.K.; Toppari, J.; Skakkebaek, N.E.; Main, K.M. Possible impact of phthalates on infant reproductive health. Int. J. Androl. 2006, 29, 172–180. [Google Scholar] [CrossRef]
  16. Ramadan, M.; Cooper, B.; Posnack, N.G. Bisphenols and phthalates: Plastic chemical exposures can contribute to adverse cardiovascular health outcomes. Birth Defects Res. 2020, 112, 1362–1385. [Google Scholar] [CrossRef]
  17. Leppert, B.; Strunz, S.; Seiwert, B.; Schlittenbauer, L.; Schlichting, R.; Pfeiffer, C.; Röder, S.; Bauer, M.; Borte, M.; Stangl, G.I.; et al. Maternal paraben exposure triggers childhood overweight development. Nat. Commun. 2020, 11, 561. [Google Scholar] [CrossRef]
  18. Ben-Jonathan, N.; Hugo, E.R. Bisphenols Come in Different Flavors: Is “S” Better Than “A”? Endocrinology 2016, 157, 1321–1323. [Google Scholar] [CrossRef]
  19. Anderson, G.; Beischlag, T.V.; Vinciguerra, M.; Mazzoccoli, G. The circadian clock circuitry and the AHR signaling pathway in physiology and pathology. Biochem. Pharmacol. 2013, 85, 1405–1416. [Google Scholar] [CrossRef]
  20. Vernet, C.; Pin, I.; Giorgis-Allemand, L.; Philippat, C.; Benmerad, M.; Quentin, J.; Calafat, A.M.; Ye, X.; Annesi-Maesano, I.; Siroux, V.; et al. In Utero Exposure to Select Phenols and Phthalates and Respiratory Health in Five-Year-Old Boys: A Prospective Study. Environ. Health Perspect. 2017, 125, 097006. [Google Scholar] [CrossRef] [Green Version]
  21. Hwang, M.; Choi, K.; Park, C. Urinary levels of phthalate, bisphenol, and paraben and allergic outcomes in children: Korean National Environmental Health Survey 2015–2017. Sci. Total Environ. 2021, 818, 151703. [Google Scholar] [CrossRef]
  22. Pugazhendhi, D.; Pope, G.S.; Darbre, P.D. Oestrogenic activity of p-hydroxybenzoic acid (common metabolite of paraben esters) and methylparaben in human breast cancer cell lines. J. Appl. Toxicol. 2005, 25, 301–309. [Google Scholar] [CrossRef] [PubMed]
  23. Fisher, M.; E Arbuckle, T.; Mallick, R.; Leblanc, A.; Hauser, R.; Feeley, M.; Koniecki, D.; Ramsay, T.; Provencher, G.; Bérubé, R.; et al. Bisphenol A and phthalate metabolite urinary concentrations: Daily and across pregnancy variability. J. Expo. Sci. Environ. Epidemiol. 2014, 25, 231–239. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Gängler, S.; Charisiadis, P.; Seth, R.; Chatterjee, S.; Makris, K.C. Time of the day dictates the variability of biomarkers of exposure to disinfection byproducts. Environ. Int. 2018, 112, 33–40. [Google Scholar] [CrossRef] [PubMed]
  25. Kim, J.H.; Kang, D.R.; Kwak, J.M.; Lee, J.K. Concentration and Variability of Urinary Phthalate Metabolites, Bisphenol A, Triclosan, and Parabens in Korean Mother–Infant Pairs. Sustainability 2020, 12, 8516. [Google Scholar] [CrossRef]
  26. Vernet, C.; Philippat, C.; Calafat, A.M.; Ye, X.; Lyon-Caen, S.; Siroux, V.; Schisterman, E.; Slama, R. Within-Day, Between-Day, and Between-Week Variability of Urinary Concentrations of Phenol Biomarkers in Pregnant Women. Environ. Health Perspect. 2018, 126, 037005. [Google Scholar] [CrossRef] [PubMed]
  27. Makris, K.C. Desynchronized circadian clock and exposures to xenobiotics are associated with differentiated disease phenotypes. BioEssays 2021, 43, 2100159. [Google Scholar] [CrossRef]
  28. Kang, S.; Kim, S.; Park, J.; Kim, H.-J.; Lee, J.; Choi, G.; Choi, S.; Kim, S.; Kim, S.Y.; Moon, H.-B.; et al. Urinary paraben concentrations among pregnant women and their matching newborn infants of Korea, and the association with oxidative stress biomarkers. Sci. Total Environ. 2013, 461–462, 214–221. [Google Scholar] [CrossRef]
  29. Zheng, R.; Dragomir, A.-C.; Mishin, V.; Richardson, J.R.; Heck, D.E.; Laskin, D.L.; Laskin, J.D. Differential metabolism of 4-hydroxynonenal in liver, lung and brain of mice and rats. Toxicol. Appl. Pharmacol. 2014, 279, 43–52. [Google Scholar] [CrossRef] [Green Version]
  30. Zarkovic, K. 4-Hydroxynonenal and neurodegenerative diseases. Mol. Asp. Med. 2003, 24, 293–303. [Google Scholar] [CrossRef]
  31. Dalleau, S.; Baradat, M.; Guéraud, F.; Huc, L. Cell death and diseases related to oxidative stress: 4-hydroxynonenal (HNE) in the balance. Cell Death Differ. 2013, 20, 1615–1630. [Google Scholar] [CrossRef] [Green Version]
  32. Paradis, V.; Kollinger, M.; Fabre, M.; Holstege, A.; Poynard, T.; Bedossa, P. In situ detection of lipid peroxidation by-products in chronic liver diseases. Hepatology 1997, 26, 135–142. [Google Scholar] [CrossRef] [PubMed]
  33. Rahman, I.; Van Schadewijk, A.A.M.; Crowther, A.J.L.; Hiemstra, P.S.; Stolk, J.; MacNee, W.; de Boer, W. 4-Hydroxy-2-Nonenal, a Specific Lipid Peroxidation Product, Is Elevated in Lungs of Patients with Chronic Obstructive Pulmonary Disease. Am. J. Respir. Crit. Care Med. 2002, 166, 490–495. [Google Scholar] [CrossRef] [PubMed]
  34. Husøy, T.; Andreassen, M.; Hjertholm, H.; Carlsen, M.; Norberg, N.; Sprong, C.; Papadopoulou, E.; Sakhi, A.; Sabaredzovic, A.; Dirven, H. The Norwegian biomonitoring study from the EU project EuroMix: Levels of phenols and phthalates in 24-hour urine samples and exposure sources from food and personal care products. Environ. Int. 2019, 132, 105103. [Google Scholar] [CrossRef] [PubMed]
  35. Papadopoulou, E.; Sabaredzovic, A.; Namork, E.; Nygaard, U.C.; Granum, B.; Haug, L.S. Exposure of Norwegian toddlers to perfluoroalkyl substances (PFAS): The association with breastfeeding and maternal PFAS concentrations. Environ. Int. 2016, 94, 687–694. [Google Scholar] [CrossRef] [PubMed]
  36. Sakhi, A.K.; Sabaredzovic, A.; Cequier, E.; Thomsen, C. Phthalate metabolites in Norwegian mothers and children: Levels, diurnal variation and use of personal care products. Sci. Total Environ. 2017, 599, 1984–1992. [Google Scholar] [CrossRef]
  37. Sabaredzovic, A.; Sakhi, A.K.; Brantsæter, A.L.; Thomsen, C. Determination of 12 urinary phthalate metabolites in Norwegian pregnant women by core-shell high performance liquid chromatography with on-line solid-phase extraction, column switching and tandem mass spectrometry. J. Chromatogr. B 2015, 1002, 343–352. [Google Scholar] [CrossRef]
  38. Sakhi, A.K.; Sabaredzovic, A.; Papadopoulou, E.; Cequier, E.; Thomsen, C. Levels, variability and determinants of environmental phenols in pairs of Norwegian mothers and children. Environ. Int. 2018, 114, 242–251. [Google Scholar] [CrossRef]
  39. Karangwa, I. Using the Markov Chain Monte Carlo Method to Make Inferences on Items of Data Contaminated by Missing Values. Am. J. Theor. Appl. Stat. 2013, 2, 48–53. [Google Scholar] [CrossRef]
  40. Moon, J.-Y.; Jung, H.-J.; Moon, M.H.; Chung, B.C.; Choi, M.H. Heat-map visualization of gas chromatography-mass spectrometry based quantitative signatures on steroid metabolism. J. Am. Soc. Mass Spectrom. 2009, 20, 1626–1637. [Google Scholar] [CrossRef] [Green Version]
  41. Sifakis, S.; Androutsopoulos, V.P.; Tsatsakis, A.M.; Spandidos, D.A. Human exposure to endocrine disrupting chemicals: Effects on the male and female reproductive systems. Environ. Toxicol. Pharmacol. 2017, 51, 56–70. [Google Scholar] [CrossRef]
  42. Bray, M.S.; Shaw, C.A.; Moore, M.W.S.; Garcia, R.A.P.; Zanquetta, M.M.; Durgan, D.J.; Jeong, W.J.; Tsai, J.-Y.; Bugger, H.; Zhang, D.; et al. Disruption of the circadian clock within the cardiomyocyte influences myocardial contractile function, metabolism, and gene expression. Am. J. Physiol. Circ. Physiol. 2008, 294, H1036–H1047. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Schindler, B.K.; Esteban, M.; Koch, H.M.; Castano, A.; Koslitz, S.; Cañas, A.; Casteleyn, L.; Kolossa-Gehring, M.; Schwedler, G.; Schoeters, G.; et al. The European COPHES/DEMOCOPHES project: Towards transnational comparability and reliability of human biomonitoring results. Int. J. Hyg. Environ. Health 2014, 217, 653–661. [Google Scholar] [CrossRef] [PubMed]
  44. Chan, E.C.Y.; Pasikanti, K.K.; Nicholson, J. Global urinary metabolic profiling procedures using gas chromatography–mass spectrometry. Nat. Protoc. 2011, 6, 1483–1499. [Google Scholar] [CrossRef] [PubMed]
  45. Curtin, P.; Austin, C.; Curtin, A.; Gennings, C.; Arora, M.; Tammimies, K.; Willfors, C.; Berggren, S.; Siper, P.; Rai, D.; et al. Dynamical features in fetal and postnatal zinc-copper metabolic cycles predict the emergence of autism spectrum disorder. Sci. Adv. 2018, 4, eaat1293. [Google Scholar] [CrossRef] [Green Version]
  46. Carroll, R.; Metcalfe, C.; Gunnell, D.; Mohamed, F.; Eddleston, M. Diurnal variation in probability of death following self-poisoning in Sri Lanka--evidence for chronotoxicity in humans. Int. J. Epidemiology 2012, 41, 1821–1828. [Google Scholar] [CrossRef] [Green Version]
  47. Imai, T. Human Carboxylesterase Isozymes: Catalytic Properties and Rational Drug Design. Drug Metab. Pharmacokinet. 2006, 21, 173–185. [Google Scholar] [CrossRef]
  48. Jewell, C.; Prusakiewicz, J.J.; Ackermann, C.; Payne, N.A.; Fate, G.; Voorman, R.; Williams, F.M. Hydrolysis of a series of parabens by skin microsomes and cytosol from human and minipigs and in whole skin in short-term culture. Toxicol. Appl. Pharmacol. 2007, 225, 221–228. [Google Scholar] [CrossRef]
  49. Obringer, C.; Wu, S.; Troutman, J.; Karb, M.; Lester, C. Effect of chain length and branching on the in vitro metabolism of a series of parabens in human liver S9, human skin S9, and human plasma. Regul. Toxicol. Pharmacol. 2021, 122, 104918. [Google Scholar] [CrossRef]
  50. Aubert, N.; Ameller, T.; Legrand, J.-J. Systemic exposure to parabens: Pharmacokinetics, tissue distribution, excretion balance and plasma metabolites of [14C]-methyl-, propyl- and butylparaben in rats after oral, topical or subcutaneous administration. Food Chem. Toxicol. 2012, 50, 445–454. [Google Scholar] [CrossRef]
  51. Zhang, Y.-K.J.; Yeager, R.L.; Klaassen, C.D. Circadian Expression Profiles of Drug-Processing Genes and Transcription Factors in Mouse Liver. Drug Metab. Dispos. 2008, 37, 106–115. [Google Scholar] [CrossRef] [Green Version]
  52. Zhao, M.; Zhang, T.; Yu, F.; Guo, L.; Wu, B. E4bp4 regulates carboxylesterase 2 enzymes through repression of the nuclear receptor Rev-erbα in mice. Biochem. Pharmacol. 2018, 152, 293–301. [Google Scholar] [CrossRef] [PubMed]
  53. Ozaki, H.; Sugihara, K.; Watanabe, Y.; Ohta, S.; Kitamura, S. Cytochrome P450-inhibitory activity of parabens and phthalates used in consumer products. J. Toxicol. Sci. 2016, 41, 551–560. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Pollock, T.; Weaver, R.E.; Ghasemi, R.; Decatanzaro, D. Butyl paraben and propyl paraben modulate bisphenol A and estradiol concentrations in female and male mice. Toxicol. Appl. Pharmacol. 2017, 325, 18–24. [Google Scholar] [CrossRef] [PubMed]
  55. Lu, Y.; Cederbaum, A.I. CYP2E1 and oxidative liver injury by alcohol. Free Radic. Biol. Med. 2008, 44, 723–738. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Seth, R.K.; Kumar, A.; Das, S.; Kadiiska, M.B.; Michelotti, G.; Diehl, A.M.; Chatterjee, S. Environmental Toxin–Linked Nonalcoholic Steatohepatitis and Hepatic Metabolic Reprogramming in Obese Mice. Toxicol. Sci. 2013, 134, 291–303. [Google Scholar] [CrossRef] [Green Version]
  57. Weber, D.; Milkovic, L.; Bennett, S.J.; Griffiths, H.R.; Zarkovic, N.; Grune, T. Measurement of HNE-protein adducts in human plasma and serum by ELISA—Comparison of two primary antibodies. Redox Biol. 2013, 1, 226–233. [Google Scholar] [CrossRef] [Green Version]
  58. Zhong, H.; Yin, H. Role of lipid peroxidation derived 4-hydroxynonenal (4-HNE) in cancer: Focusing on mitochondria. Redox Biol. 2014, 4, 193–199. [Google Scholar] [CrossRef] [Green Version]
  59. Łuczaj, W.; Gindzienska-Sieskiewicz, E.; Jarocka-Karpowicz, I.; Andrisic, L.; Sierakowski, S.; Zarkovic, N.; Waeg, G.; Skrzydlewska, E. The onset of lipid peroxidation in rheumatoid arthritis: Consequences and monitoring. Free Radic. Res. 2016, 50, 304–313. [Google Scholar] [CrossRef]
  60. Castro, J.P.; Jung, T.; Grune, T.; Siems, W. 4-Hydroxynonenal (HNE) modified proteins in metabolic diseases. Free Radic. Biol. Med. 2017, 111, 309–315. [Google Scholar] [CrossRef] [Green Version]
  61. Steffensen, I.-L.; Dirven, H.; Couderq, S.; David, A.; D’Cruz, S.C.; Fernández, M.F.; Mustieles, V.; Rodríguez-Carrillo, A.; Hofer, T. Bisphenols and Oxidative Stress Biomarkers—Associations Found in Human Studies, Evaluation of Methods Used, and Strengths and Weaknesses of the Biomarkers. Int. J. Environ. Res. Public Health 2020, 17, 3609. [Google Scholar] [CrossRef]
  62. Wang, Y.-X.; Liu, C.; Shen, Y.; Wang, Q.; Pan, A.; Yang, P.; Chen, Y.-J.; Deng, Y.-L.; Lu, Q.; Cheng, L.-M.; et al. Urinary levels of bisphenol A, F and S and markers of oxidative stress among healthy adult men: Variability and association analysis. Environ. Int. 2018, 123, 301–309. [Google Scholar] [CrossRef] [PubMed]
  63. Asimakopoulos, A.G.; Xue, J.; De Carvalho, B.P.; Iyer, A.; Abualnaja, K.O.; Yaghmoor, S.S.; Kumosani, T.A.; Kannan, K. Urinary biomarkers of exposure to 57 xenobiotics and its association with oxidative stress in a population in Jeddah, Saudi Arabia. Environ. Res. 2016, 150, 573–581. [Google Scholar] [CrossRef] [PubMed]
  64. Ferguson, K.; Cantonwine, D.E.; McElrath, T.F.; Mukherjee, B.; Meeker, J.D. Repeated measures analysis of associations between urinary bisphenol-A concentrations and biomarkers of inflammation and oxidative stress in pregnancy. Reprod. Toxicol. 2016, 66, 93–98. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Watkins, D.J.; Ferguson, K.; Del Toro, L.V.A.; Alshawabkeh, A.N.; Cordero, J.F.; Meeker, J.D. Associations between urinary phenol and paraben concentrations and markers of oxidative stress and inflammation among pregnant women in Puerto Rico. Int. J. Hyg. Environ. Health 2014, 218, 212–219. [Google Scholar] [CrossRef] [Green Version]
  66. Rocha, B.A.; Asimakopoulos, A.G.; Honda, M.; da Costa, N.L.; Barbosa, R.M.; Barbosa, F.; Kannan, K. Advanced data mining approaches in the assessment of urinary concentrations of bisphenols, chlorophenols, parabens and benzophenones in Brazilian children and their association to DNA damage. Environ. Int. 2018, 116, 269–277. [Google Scholar] [CrossRef]
  67. Lv, Y.; Rui, C.; Dai, Y.; Pang, Q.; Li, Y.; Fan, R.; Lu, S. Exposure of children to BPA through dust and the association of urinary BPA and triclosan with oxidative stress in Guangzhou, China. Environ. Sci. Process. Impacts 2016, 18, 1492–1499. [Google Scholar] [CrossRef]
  68. Ding, S.; Qi, W.; Xu, Q.; Zhao, T.; Li, X.; Yin, J.; Zhang, R.; Huo, C.; Zhou, L.; Ye, L. Relationships between di-(2-ethylhexyl) phthalate exposure and lipid metabolism in adolescents: Human data and experimental rat model analyses. Environ. Pollut. 2021, 286, 117570. [Google Scholar] [CrossRef]
  69. Aseervatham, G.S.B.; Sivasudha, T.; Jeyadevi, R.; Ananth, D.A. Environmental factors and unhealthy lifestyle influence oxidative stress in humans—An overview. Environ. Sci. Pollut. Res. 2013, 20, 4356–4369. [Google Scholar] [CrossRef]
  70. Rytilä, P.; Rehn, T.; Ilumets, H.; Rouhos, A.; Sovijärvi, A.; Myllärniemi, M.; Kinnula, V.L. Increased oxidative stress in asymptomatic current chronic smokers and GOLD stage 0 COPD. Respir. Res. 2006, 7, 69. [Google Scholar] [CrossRef] [Green Version]
  71. Kode, A.; Yang, S.-R.; Rahman, I. Differential effects of cigarette smoke on oxidative stress and proinflammatory cytokine release in primary human airway epithelial cells and in a variety of transformed alveolar epithelial cells. Respir. Res. 2006, 7, 132. [Google Scholar] [CrossRef] [Green Version]
  72. Sbrana, E.; Suter, M.A.; Abramovici, A.R.; Hawkins, H.K.; Moss, J.E.; Patterson, L.; Shope, C.; Aagaard-Tillery, K. Maternal tobacco use is associated with increased markers of oxidative stress in the placenta. Am. J. Obstet. Gynecol. 2011, 205, 246.e1–246.e7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Minato, T.; Tsutsumi, M.; Tsuchishima, M.; Hayashi, N.; Saito, T.; Matsue, Y.; Toshikuni, N.; Arisawa, T.; George, J. Binge Alcohol Consumption Aggravates Oxidative Stress and Promotes Pathogenesis of NASH from Obesity-Induced Simple Steatosis. Mol. Med. 2014, 20, 490–502. [Google Scholar] [CrossRef] [PubMed]
  74. Watt, J.; Schuon, J.; Davis, J.; Ferguson, T.F.; Welsh, D.A.; Molina, P.E.; Ronis, M.J.J. Reduced Serum Osteocalcin in High-Risk Alcohol Using People Living With HIV Does Not Correlate With Systemic Oxidative Stress or Inflammation: Data From the New Orleans Alcohol Use in HIV Study. Alcohol. Clin. Exp. Res. 2019, 43, 2374–2383. [Google Scholar] [CrossRef] [PubMed]
  75. Wegner, S.A.; Pollard, K.A.; Kharazia, V.; Darevsky, D.; Perez, L.; Roychowdhury, S.; Xu, A.; Ron, R.; Nagy, L.E.; Hopf, F.W. Limited Excessive Voluntary Alcohol Drinking Leads to Liver Dysfunction in Mice. Alcohol. Clin. Exp. Res. 2017, 41, 345–358. [Google Scholar] [CrossRef] [Green Version]
  76. Raza, H.; John, A. 4-Hydroxynonenal induces mitochondrial oxidative stress, apoptosis and expression of glutathione S-transferase A4-4 and cytochrome P450 2E1 in PC12 cells. Toxicol. Appl. Pharmacol. 2006, 216, 309–318. [Google Scholar] [CrossRef]
  77. Das, S.; Kumar, A.; Seth, R.; Tokar, E.J.; Kadiiska, M.B.; Waalkes, M.P.; Mason, R.P.; Chatterjee, S. Proinflammatory adipokine leptin mediates disinfection byproduct bromodichloromethane-induced early steatohepatitic injury in obesity. Toxicol. Appl. Pharmacol. 2013, 269, 297–306. [Google Scholar] [CrossRef] [Green Version]
  78. Tomasi, A.; Albano, E.; Biasi, F.; Slater, T.F.; Vannini, V.; Dianzani, M.U. Activation of chloroform and related trihalomethanes to free radical intermediates in isolated hepatocytes and in the rat in vivo as detected by the ESR-spin trapping technique. Chem. Interact. 1985, 55, 303–316. [Google Scholar] [CrossRef]
  79. Seth, R.K.; Das, S.; Dattaroy, D.; Chandrashekaran, V.; Alhasson, F.; Michelotti, G.; Nagarkatti, M.; Nagarkatti, P.; Diehl, A.M.; Bell, P.D.; et al. TRPV4 activation of endothelial nitric oxide synthase resists nonalcoholic fatty liver disease by blocking CYP2E1-mediated redox toxicity. Free Radic. Biol. Med. 2016, 102, 260–273. [Google Scholar] [CrossRef]
  80. Gkrillas, A.; Dirven, H.; Papadopoulou, E.; Andreassen, M.; Hjertholm, H.; Husøy, T. Exposure estimates of phthalates and DINCH from foods and personal care products in comparison with biomonitoring data in 24-hour urine from the Norwegian EuroMix biomonitoring study. Environ. Int. 2021, 155, 106598. [Google Scholar] [CrossRef]
  81. Haddad, N.; Andrianou, X.D.; Makris, K.C. A Scoping Review on the Characteristics of Human Exposome Studies. Curr. Pollut. Rep. 2019, 5, 378–393. [Google Scholar] [CrossRef] [Green Version]
  82. Vermeulen, R.; Schymanski, E.L.; Barabási, A.-L.; Miller, G.W. The exposome and health: Where chemistry meets biology. Science 2020, 367, 392–396. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Graphical presentation of the human biomonitoring study design.
Figure 1. Graphical presentation of the human biomonitoring study design.
Toxics 10 00181 g001
Figure 2. Forest plot visualization of the linear mixed-effect models of SG-adjusted metabolites of phthalates/DINCH as a function of time of urine sampling in the day adjusted for dietary habits and personal-care-product use (reference: mid-day = zero point).
Figure 2. Forest plot visualization of the linear mixed-effect models of SG-adjusted metabolites of phthalates/DINCH as a function of time of urine sampling in the day adjusted for dietary habits and personal-care-product use (reference: mid-day = zero point).
Toxics 10 00181 g002
Figure 3. Forest plot visualization of the linear mixed-effect models of SG-adjusted parabens and bisphenols as a function of time of urine sampling in the day adjusted for dietary habits and personal-care-product use (reference: mid-day = zero point on x-axis).
Figure 3. Forest plot visualization of the linear mixed-effect models of SG-adjusted parabens and bisphenols as a function of time of urine sampling in the day adjusted for dietary habits and personal-care-product use (reference: mid-day = zero point on x-axis).
Toxics 10 00181 g003
Table 1. Demographic characteristics of the study participants.
Table 1. Demographic characteristics of the study participants.
Basic CharacteristicsOverall (n = 144)Males (n = 44)Females (n = 100)
Age (years, mean ± SD)42.56± 12.0543.36 ± 11.6542.20 ± 12.26
Weight (kg, mean ± SD)70.45 ± 11.7181.98 ± 8.5465.22 ± 8.88
Height (m, mean ± SD)1.72 ± 0.081.81 ± 0.061.68 ± 0.06
BMI (kg/m2, mean ± SD)23.47 ± 3.5224.94 ± 2.3322.82 ± 3.77
Smoking status (n, %)
Non-smokers90 (62.5)26 (59.1)64 (64)
Ex-smokers35 (24.3)11 (25.0)24 (24)
Occasional smokers19 (13.2)7 (15.9)12 (12)
Education (n, %)
University/college
up to 4 years30 (20.8)8 (18.2)22 (22)
>4 years114 (79.2)36 (81.8)78 (78)
Table 2. Percentiles of the specific gravity-adjusted biomarkers of exposure to EDCs and the biomarker of effect/lipid peroxidation (4HNE) by the time of urine sampling.
Table 2. Percentiles of the specific gravity-adjusted biomarkers of exposure to EDCs and the biomarker of effect/lipid peroxidation (4HNE) by the time of urine sampling.
BiomarkersTime GroupMin5th25th50th75th90th95thMax
MeP (ug)6:00–12:000.4271.0753.0135.93712.58223.52144.215418
12:00–18:000.4060.8151.7513.8728.52418.27734.0461990
18:00–6:000.6071.9234.0998.21918.51734.43860.88364
Overall 0.406 2.8015.95512.715 1990
MiBP (ug)6:00–12:000.5041.0842.3944.3416.96911.78917.681106
12:00–18:000.110.8221.542.5695.0057.90610.94723.5
18:00–6:000.9172.2143.9075.7389.37117.1524.39978.4
Overall 0.110 2.2894.3707.340 106.4
MnBP (ug)6:00–12:000.8091.5283.1195.1749.14717.03321.59142.1
12:00–18:000.4790.9792.1293.4315.3789.06111.9720.6
18:00–6:001.4762.5375.0368.03614.62121.36224.60330.5
Overall 0.480 2.9955.1929.531 42.1
MBzP (ug)6:00–12:000.0640.2210.4180.6861.1672.0173.24740.9
12:00–18:000.0430.1330.3080.4780.8181.5052.1748.7
18:00–6:000.1660.3630.7101.0681.9112.5453.53240.6
Overall 0.043 0.4250.7331.248 40.9
SumDEHP (ug)6:00–12:002.3793.528.46312.52622.37233.67242.4364.7
12:00–18:001.4553.0796.0419.53714.74126.00735.354176.8
18:00–6:002.6656.97112.97919.20130.55944.15658.501127.7
Overall 1.455 8.46713.29523.315 176.8
SumDiNP (ug)6:00–12:000.5681.3682.9165.44810.46620.74338.159526
12:00–18:000.6611.1872.4354.8057.99814.91529.127156
18:00–6:000.6442.6455.1968.27115.65337.22169.2391276
Overall 0.568 3.1856.21911.136 1276
SumDINCH (ug)6:00–12:000.0660.1580.3780.6911.2392.7045.52612.5
12:00–18:000.0250.120.2950.4830.9571.6372.3457.8
18:00–6:000.1240.2220.5550.9411.6542.7624.38993.3
Overall 0.025 0.3990.7101.216 93.3
oh-MPHP (ug)6:00–12:000.0050.020.0780.1960.3810.5820.7712.7
12:00–18:000.0070.0170.070.160.2990.5070.6841.18
18:00–6:000.0180.0380.1350.3090.5380.8511.3155.95
Overall 0.005 0.0840.2080.410 5.95
MEPA (ug/L)6:00–12:000.8161.5123.3946.30537.19299.038158.27518519
12:00–18:000.3361.112.9545.8524.27381.163134.8533960
18:00–6:000.7561.4563.2777.32027.65387.744123.8856728
Overall 0.336 3.2466.46725.362 18519
ETPA (ug/L)6:00–12:000.0240.1640.5201.2353.43214.11826.195761.7
12:00–18:000.0260.1440.3901.0833.19412.17736.465271.1
18:00–6:000.0290.1450.3480.8413.26919.47140.1272038
Overall 0.024 0.4221.0833.264 2038
PRPA (ug/L)6:00–12:000.0080.0120.0240.2512.82720.07646.206384.9
12:00–18:000.0080.0120.0240.0951.15910.66019.738225.7
18:00–6:000.0080.0120.0270.1693.40416.41632.363175.7
Overall 0.008 0.0240.1732.366 384.9
BUPA (ug/L)6:00–12:000.0140.0180.0320.0650.1790.5281.0839.82
12:00–18:000.0160.0190.0410.0740.1770.3470.5855.45
18:00–6:000.0140.0190.0320.0510.1620.3640.68147.1
Overall 0.014 0.0350.0650.176 47.1
BPA (ug/L)6:00–12:000.0140.2190.6211.3932.4393.9655.64621.1
12:00–18:000.0110.0600.6021.0592.1754.1764.98324.6
18:00–6:000.0140.1840.6781.1701.9653.7745.31919.6
Overall 0.011 0.6211.1632.194 24.6
BPS (ug/L)6:00–12:00<LODLOD<LOD<LOD<LOD0.5011.0832.8
12:00–18:00<LOD<LOD<LOD<LOD0.1820.4550.6937.04
18:00–6:00<LOD<LOD<LOD<LOD0.1650.6911.08255.7
Overall <LOD<LOD<LOD<LOD0.138 55.7
BPF (ug/L)6:00–12:00<LOD<LOD<LOD<LOD<LOD<LOD<LOD6.0
12:00–18:00<LOD<LOD<LOD<LOD<LOD<LOD<LOD7.7
18:00–6:00<LOD<LOD<LOD<LOD<LOD<LOD0.58012.8
Overall <LOD <LOD<LOD<LOD 12.8
4HNE (ug/L)6:00–12:003.0283.2623.7634.0424.364.5834.7845.6
12:00–18:002.7663.2893.8134.0834.4584.7715.0265.8
18:00–6:002.9053.3363.7624.0444.2944.6184.9486.5
Overall 2.766 3.7754.0484.375 6.5
Table 3. Mixed-effect models for the association between the SG-adjusted and logarithm transformed biomarkers of exposure to EDCs and the biomarker of effect (4HNE).
Table 3. Mixed-effect models for the association between the SG-adjusted and logarithm transformed biomarkers of exposure to EDCs and the biomarker of effect (4HNE).
Log(4HNE)
Parent EDCsBiomarkerEstimate95% CIq Value
PhthalatesMeP0.090.05–0.13<0.001
MiBP0.140.09–0.20<0.001
MnBP0.1430.088–0.199<0.001
MBzP0.1110.059–0.163<0.001
SumDEHP0.1830.125–0.240<0.001
SumDiNP0.1030.061–0.146<0.001
SumDINCH0.1130.065–0.161<0.001
oh-MPHP0.1250.083–0.167<0.001
ParabensMEPA0.130.10–0.17<0.001
ETPA0.0960.064–0.128<0.001
PRPA0.0390.017–0.061<0.001
BUPA0.1230.064–0.181<0.001
BisphenolsBPA0.1280.091–0.165<0.001
BPS−0.015−0.114–0.0840.767
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Abimbola, S.O.; Xeni, C.; Sakhi, A.K.; Heibati, B.; Husøy, T.; Dirven, H.; Makris, K.C. Diurnal Variation in Biomarkers of Exposure to Endocrine-Disrupting Chemicals and Their Association with Oxidative Damage in Norwegian Adults: The EuroMix Study. Toxics 2022, 10, 181. https://doi.org/10.3390/toxics10040181

AMA Style

Abimbola SO, Xeni C, Sakhi AK, Heibati B, Husøy T, Dirven H, Makris KC. Diurnal Variation in Biomarkers of Exposure to Endocrine-Disrupting Chemicals and Their Association with Oxidative Damage in Norwegian Adults: The EuroMix Study. Toxics. 2022; 10(4):181. https://doi.org/10.3390/toxics10040181

Chicago/Turabian Style

Abimbola, Samuel Olushola, Christina Xeni, Amrit Kaur Sakhi, Behzad Heibati, Trine Husøy, Hubert Dirven, and Konstantinos C. Makris. 2022. "Diurnal Variation in Biomarkers of Exposure to Endocrine-Disrupting Chemicals and Their Association with Oxidative Damage in Norwegian Adults: The EuroMix Study" Toxics 10, no. 4: 181. https://doi.org/10.3390/toxics10040181

APA Style

Abimbola, S. O., Xeni, C., Sakhi, A. K., Heibati, B., Husøy, T., Dirven, H., & Makris, K. C. (2022). Diurnal Variation in Biomarkers of Exposure to Endocrine-Disrupting Chemicals and Their Association with Oxidative Damage in Norwegian Adults: The EuroMix Study. Toxics, 10(4), 181. https://doi.org/10.3390/toxics10040181

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop