Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Urinary Phthalate Metabolites Are Associated with Body Mass Index and Waist Circumference in Chinese School Children

  • Hexing Wang,

    Affiliation Key Laboratory of Public Health Safety of Ministry of Education, School of Public Health, Fudan University, Shanghai, China

  • Ying Zhou ,

    yingchou@fudan.edu.cn

    Affiliation Key Laboratory of Public Health Safety of Ministry of Education, School of Public Health, Fudan University, Shanghai, China

  • Chuanxi Tang,

    Affiliation Centers for Disease Control and Prevention of Changning District, Shanghai, China

  • Yanhong He,

    Affiliation Department of Laboratory Medicine, Anting Hospital, Shanghai, China

  • Jingui Wu,

    Affiliation Centers for Disease Control and Prevention of Changning District, Shanghai, China

  • Yue Chen,

    Affiliation Department of Epidemiology and Community Medicine, Faculty of Medicine, University of Ottawa, Ottawa, Ontario, Canada

  • Qingwu Jiang

    Affiliation Key Laboratory of Public Health Safety of Ministry of Education, School of Public Health, Fudan University, Shanghai, China

Abstract

Background

Lab studies have suggested that ubiquitous phthalate exposures are related to obesity, but relevant epidemiological studies are scarce, especially for children.

Objective

To investigate the association of phthalate exposures with body mass index (BMI) and waist circumference (WC) in Chinese school children.

Methods

A cross-sectional study was conducted in three primary and three middle schools randomly selected from Changning District of Shanghai City of China in 2011–2012. According to the physical examination data in October, 2011, 124 normal weight, 53 overweight, and 82 obese students 8–15 years of age were randomly chosen from these schools on the basis of BMI-based age- and sex-specific criterion. First morning urine was collected in January, 2012, and fourteen urine phthalate metabolites (free plus conjugated) were determined by ultra-performance liquid chromatography coupled to tandem mass spectrometry. Multiple linear regression was used to explore the associations between naturally log-transformed urine phthalate metabolites and BMI or WC.

Results

The urine specific gravity-corrected concentrations of nine urine phthalate metabolites and five molar sums were positively associated with BMI or WC in Chinese school children after adjustment for age and sex. However, when other urine phthalate metabolites were included in the models together with age and sex as covariables, most of these significant associations disappeared except for mono (2-ethylhexyl) phthalate (MEHP) and monoethyl phthalate (MEP). Additionally, some associations showed sex- or age-specific differences.

Conclusions

Some phthalate exposures were associated with BMI or WC in Chinese school children. Given the cross-sectional nature of this study and lack of some important obesity-related covariables, further studies are needed to confirm the associations.

Introduction

Phthalates, the diesters of phthalic acid, are extensively used as plasticizers and solvents in industrial, medical, and consumer products [1]. Applications of phthalates depend on their chemical structures, by which they are classified into two groups of high and low molecular weight [2]. High molecular weight phthalates, such as di(2-ethylhexyl) phthalate (DEHP), butyl benzyl phthalate (BBzP), di-iso-nonyl phthalate (DiNP), di-iso-decyl phthalate (DiDP) and dioctyl phthalate (DOP), are primarily used in the production of polyvinylchloride plastics (PVC) [3], [4]. Low molecular weight phthalates, such as dimethyl phthalate (DMP), diethyl phthalate (DEP), di-n-butyl phthalate (DBP), and di-iso-butyl phthalate (DiBP), are often used in making personal care products (e.g. perfumes, lotions and cosmetics), varnishes, and coatings [5], [6]. Humans can be exposed to phthalates through ingestion, inhalation, and dermal contact as well as during intrauterine development through placenta [7]. As a result of ubiquitous use of these phthalate-related products, human exposures to phthalates have been considered to be extensive [8], [9]. Experimental studies have demonstrated that some phthalates are endocrine disrupting chemicals (EDCs), and have reproductive or developmental toxicity [10].

There is increasing evidences suggesting that some EDCs can act as obesogens and interfere with body’s natural weight control mechanisms by disrupting adipose tissue biology, endocrine hormone systems, or central hypothalamic-pituitary-adrenal axis [11], [12]. For this, the National Toxicology Program (NTP) organized a special workshop to discuss the phthalates as potential obesogen and recommended subsequent research strategy in 2011 [13]. Lab and population studies have suggested that phthalate exposures can promote adipogenesis by diverse biologic mechanisms, such as interfering with the functions of steroid hormones or thyroid hormones, and the inappropriate activation of peroxisome proliferator-activated receptors (PPARs) [14][17].

Sex hormones are critical to the regulation of body composition, fuel homeostasis, and reproduction in humans [18]. Some lab studies have observed that DEHP and (mono(2-ethylhexyl) phthalate (MEHP) decrease the biosynthesis of steroid hormones, such as estradiol and testosterone, by suppressing aromatase transcript levels [19], [20]. Some phthalates, such as DiBP, DBP, dicyclohexyl phthalate (DCHP) and DEHP, were found to activate human estrogen receptors, human androgen receptors, and glucocorticoid receptors [21]. Thyroid function is critical to the maintenance of basal metabolism. Shen et al. found that DBP and MBP altered expression of selected thyroid hormone response genes by Xenopus laevis [22]. Epidemiological studies also observed an inverse association between urinary phthalate metabolites, such as monoethyl phthalate (MEP), mono-n-butyl phthalate (MBP), monobenzyl phthalate (MBzP), and MEHP, and thyroid hormone serum levels in children and in men [23], [24]. Kitahara et al. observed that levels of serum TSH and free T3 were positively associated with body mass index (BMI) and waist circumference (WC) among euthyroid adults [25]. PPARs play a major role in adipocyte differentiation and energy storage. Bility et al. observed that several phthalate metabolites, such as MBzP, monoisononyl phthalate (MiNP), monooctyl phthalate (MOP) and MEHP, could activate mouse and human PPAR-α, PPAR-β, and PPAR-γ although the ability of these metabolites to activate these receptors considerably varied [26]. Feige et al. reported that MEHP was a selective PPAR-γ modulator that promotes adipogenesis [27].

However, epidemiological studies that investigate the associations of postnatal phthalate exposures with obesity or body size measurements are scarce. To our knowledge, there were four epidemiological studies for this, including two cross-sectional studies and one longitudinal study involved in children, another cross-sectional one in men [24], [28][30]. Among these studies, phthalate metabolites, such as MEP, MiBP, MBP, MBzP, and some metabolites of DEHP, were found to be significantly associated with BMI or WC.

Children and adolescents are in rapid growth and development and are more susceptible to external disturbance than adults. Moreover, it has been reported that children are extensively exposed to phthalates and the exposure levels are higher compared with adults [31], [32]. Body mass index is a reasonable measure of fatness in children and adolescents and used as an indicator of overall obesity, while WC could reflect the accumulation of abdominal fat tissue and is used as an indicator of central obesity [33], [34]. In this study we examined the association of urine phthalate metabolite concentrations with BMI and WC in Chinese school children.

Materials and Methods

Ethics Statement

The study was approved by the Ethical Review Board of the School of Public Health of Fudan University. Guardians provided written informed consent before the study.

Study Population

The study design and study population has been described in detail in a previous report of the association of urine bisphenol a (BPA) with BMI [35]. Briefly, three primary schools and three middle schools were randomly selected from 26 primary and 30 middle candidate schools in Changning District of Shanghai City in China. From each selected school, 20 obese, 10 overweight, and 30 normal weight students aged 8–15 years were randomly chosen in the basis of the most recent yearly regular physical examination conducted in October, 2011. A total of 360 eligible students were selected and they were all Han people. Of these 360 subjects, 53 did not provide urine samples and 48 provided urine samples but not first morning urine. Therefore, 259 participants were available for the current analysis. The distributions of participants by age, sex, and BMI subgroups (normal weight, overweight, and obesity) were listed in Table 1.

thumbnail
Table 1. Arithmetic mean (standard deviation) of waist circumference (WC) and body mass index (BMI) in 259 school children of Shanghai, China in 2011–2012.

https://doi.org/10.1371/journal.pone.0056800.t001

Anthropometric Measurements

Body weight (kg) and height (cm) were measured using the same type of apparatus and followed the standard procedures recommended by Cameron [36]. Each technician was required to pass a training course for the anthropometric measurements. Subjects were required to wear light clothes and stand straight and barefoot when being measured. Waist circumference (cm) was measured midway between the lowest rib and the superior border of the iliac crest with an inelastic measuring tape at the end of normal expiration to the nearest 0.1 cm [34]. Body mass index was calculated as weight in kilograms divided by height in meters squared. Normal weight, overweight, and obese participants in this study were identified according to BMI-based age- and sex-specific criterion proposed by the Working Group on Obesity in China (WGOC). This criterion was based on a large representative Chinese population 7–18 years of age and the specific BMI cut-off values were listed in Table S1 in Supporting Information File S1 [37].

Measurement of Urine Phthalate Metabolites

For a better comparability of urine phthalate metabolites among individuals, first morning urine was collected and participants themselves carried out the urine collection. In order to minimize the contamination of phthalates during urine sampling and storing, urine samples were collected in glass centrifuge tubes, which was rinsed by acetone and baked at 350°C for 2 hours, and the participants were told to avoid a contact of urine with plastic products in the process of urine collection. The collected urine samples were transported to lab as soon as possible and stored in dark at −20°C until analysis. All the urine samplings were collected between January 3 and 6, 2012.

Fourteen phthalate metabolites studied here included MEHP, mono(2-ethyl-5-oxohexyl) phthalate (MEOHP), mono(2-ethyl-5-carboxypentyl) phthalate (MECPP), mono(2-ethyl-5-hydroxyhexyl) phthalate (MEHHP), MBP, monoisobutyl phthalate (MiBP), MEP, mono[(2-carboxymethyl) hexyl] phthalate (MCMHP), mono(4-hydroxybutyl) phthalate (MHBP), monomethyl phthalate (MMP), monocyclohexyl phthalate (MCHP), MBzP, monoisononyl phthalate (MiNP), and monooctyl phthalate (MOP). MEHP is the primary metabolite of DEHP, and MEHHP, MECPP, MCMHP, and MEOHP are secondary metabolites of DEHP. MBP and MHBP are primary and secondary metabolites of DBP, respectively. MiBP is derived from DiBP, MMP from DMP, MEP from DEP, MCHP from DCHP, MBzP from BBzP, MiNP from DiNP, and MOP from DOP.

Fourteen urine phthalate metabolites (free and conjugated) were analyzed using a method modified from Guo et al. [1]. Briefly, after an aliquot (1.0 mL) of urine sample was enzymatically hydrolyzed and purified by solid-phase extraction, the phthalate metabolites in urine extract were resolved by reversed-phase ultra-performance liquid chromatography, detected by electrospray ionization tandem mass spectrometry, and quantified by an isotope internal standard curve method.

All urine analyses were carried out in random and information on BMI subgroup was blinded to analytical technicians. To minimize a potential systematic drift caused by analytical process, all the analyses of urine samples were performed in a short time span (all the analyses of urine samples were completed in May 2012) and by the same analytical team. The recoveries of phthalate metabolites of spiked urine samples used for quality control in the process of whole urine analyses varied between 70.3% and 124%, and their standard deviations varied between 9.2% and 19.1%. The detailed analytical method was supplied in Supporting Information File S1 (Tables S2 and S3 in Supporting Information File S1).

Statistical Analyses

For measurements below the limit of detection (LOD, see supplementary information File S1), a default value of LOD divided by the square root of 2 was assigned, which is a method of handling nondetectable values that produces reasonably nonbiased means and standard deviations [38]. In addition to twelve individual phthalate metabolites, we combined some phthalate metabolites into molar sums based on their molecular-weight or parent phthalates. Sum of low molecular-weight phthalate metabolites (MWP) comprised of MBP, MHBP, MiBP, MMP, and MEP, and sum of high MWP comprised of five metabolites of DEHP (MEHP, MECPP, MEHHP, MEOHP, MCMHP), MCHP, and MBzP. By respective molar sum of each metabolites of constructed sums, sum of metabolites of DBP (MBP and MHBP) and sum of low MWP were converted to MBP (molecular weight 222) mass concentrations (ng/mL), and sum of metabolites of DEHP (MEHP, MEHHP, MECPP, MCMHP, and MEOHP), sum of high MWP, and sum of all phthalate metabolites were converted to MEHP (molecular weight 278) mass concentrations (ng/mL). These sums have the same unit as phthalate metabolites.

Due to the glucuronidation of phthalate metabolites in the liver and its elimination by active tubular secretion, creatinine correction may not be appropriate for urine phthalate metabolite concentration [8], [39]. Additionally, BMI could predict the urine creatinine concentration, and therefore, the associations between phthalate metabolites and BMI may be altered by creatinine adjustment [40]. For this reason, we used specific gravity to correct for urinary dilution, as recommended by Hauser et al. [39]. Specific gravity was measured using a handheld refractometer (Atago PAL 10-S, Tokyo, Japan). The correction formula was Pc = P × [(1.024–1)/(SG –1)], where Pc is the specific gravity-corrected phthalate metabolite concentration (ng/mL), P is the experimental phthalate metabolite concentration (ng/mL), and SG is the specific gravity of the urine sample [41].

For descriptive statistics, arithmetic mean and standard deviation (SD) of BMI and WC, and geometric mean (GM), 95% confidence interval (CI), and median of SG-uncorrected urine phthalate metabolite concentrations were calculated for all subjects or by sex, age, and BMI subgroups (normal weight, overweight, and obesity) (Table 1). Because the urinary concentration values of phthalate metabolite were not normally distributed, a natural log-transformation was applied to normalize the data. Analysis of variance was used to investigate the association of SG-uncorrected urine phthalate metabolite concentrations with age (8–11 and 12–15 age groups), sex, and BMI subgroups (normal weight, overweight, and obesity) in all subjects (Tables 2 and 3). When heterogeneity of variance was encountered, rank test was used instead of analysis of variance.

thumbnail
Table 2. Geometric mean (median), ng/mL) of urine phthalate metabolite concentrations before urine specific gravity correction associated with age and sex in 259 school children of Shanghai, China in 2011–2012.

https://doi.org/10.1371/journal.pone.0056800.t002

thumbnail
Table 3. Geometric mean (median), ng/mL) of urine phthalate metabolite concentrations before specific gravity correction associated with excess body weight after adjustment for sex and age in 259 school children of Shanghai, China in 2011–2012.

https://doi.org/10.1371/journal.pone.0056800.t003

Multiple linear regres­sion analysis was also performed to study the association of SG-corrected urine phthalate metabolite concentrations with BMI and WC values (Tables 4 and 5). Linear regression analysis included only phthalate metabolites detectable in more than 50% of children. Because a fixed proportion of BMI subgroups (30 normal weight, 10 overweight, and 20 obese children) was used to select participants from each school with different selection probabilities for the BMI subgroups, it might generate a bias when estimating the associations in combined population of three BMI subgroups. We used sampling weights to calculate weighted parameter estimates in multiple linear regres­sion analysis [42], [43]. Because urinary phthalate metabolite concentrations were correlated with each other, especially between metabolites derived from one parent phthalate (see Table S4 in Supporting Information File S1), urine phthalate metabolites were first examined separately in crude analysis and adjusted model 1 with age and sex included as covariables, and then, in adjusted model 2, urine phthalate metabolites were also included as covariables besides age and sex. In addition to the regression models including all the subjects, stratified modeling by age (8–11 and 12–15 age groups) and sex was also conducted. The regression coefficient (β) represented the change of naturally log-transformed BMI or WC value per one unit of naturally log-transformed urine phthalate metabolite concentration. Analysis of variance was used to test for the interactions of urine phthalate metabolites with age subgroups (8–11 and 12–15 years of age) or sex by modeling them as main effects and interaction terms of naturally log-transformed BMI.

thumbnail
Table 4. Log-transformed urinary phthalate metabolite concentrations corrected by specific gravity (SG) associated with body mass index (BMI) in 259 school children of Shanghai, China in 2011–2012.

https://doi.org/10.1371/journal.pone.0056800.t004

thumbnail
Table 5. Log-transformed urinary phthalate metabolite concentrations corrected by specific gravity (SG) associated with waist circumference (WC) in 259 school children of Shanghai, China in 2011–2012.

https://doi.org/10.1371/journal.pone.0056800.t005

In adjusted model 2, in order to avoid the collinearity derived from strong correlations among daughter metabolites of one parent phthalates in linear regression model, only one daughter metabolite was included, whereas all others not. In this model, sum of DEHP was used to represent for five metabolites of DEHP, and sum of DBP was used to represent for MBP and MHBP. When one of the phthalate metabolites from DEHP or DBP was included in this model, sum of DBP or sum of DEHP was not added as covariable into the model. In addition, sum of high MWP, sum of low MWP, and sum of all were not used as independent covariables in all these models. The detail information about covariable for specific metabolites modelling was provided in Tables 4 and 5. Data analysis was performed using SPSS (version 17; SPSS, Inc., Chicago, IL, USA); p<0.05 was considered significant, and all statistical tests were two-sided.

Results

Table 1 shows demographic characteristics and anthropometric measurements by age (8–9, 10–11, 12–13, and 14–15 years), sex and BMI subgroups (normal weight, overweight, and obesity). As expected, WC and/or BMI increased with age and BMI subgroups (p<0.001). There were no significant differences between males and females.

As shown in Table 2, phthalate exposure was prevalent. Except for MBzP, detection frequency of other eleven metabolites was 94.6% (MCHP) or above. Analysis of variance showed that four (MEHP, MCHP, MMP, and MEP) of twelve urine phthalate metabolites were significantly higher in children 12–15 years of age than those 8–11 years of age, and three (MCMHP, MBP, and MHBP) were marginally higher (p<0.10). No significant sex-related difference was observed for all urine phthalate metabolites. MOP and MiNP were not detected in all urine samples.

Table 3 shows the association of urinary metabolites with BMI grouped as normal weight, overweight, and obesity. Analysis of variance and Kruskai-Wallis rank test indicated that seven of eleven phthalate metabolites and five sums were significantly different among BMI groups. The average levels of MEHP, MiBP, MEP, sum of DBP, sum of low MW, and sum of all metabolites showed an increasing trend with BMI.

In multiple linear regression analysis, similar results were obtained before and after sampling weights were accounted for, and weighted estimates were presented. The detail sampling weights of each selected school was provided in Table S5 in Supporting Information File S1. Tables 4 and 5 showed that the regression results for the association of urine phthalate metabolites with BMI and WC, respectively. As shown in crude analysis, adjusted model 1, and adjusted model 2, the associations of phthalate metabolites with BMI or WC became weaker, and some association directions were even reversed after adjustment for more covariables.

Table 4 indicated that seven of eleven metabolites and five sums were positively associated with increasing BMI after adjustment for age and sex in adjusted model 1, but after additional adjustment for urine phthalate metabolites in adjusted model 2, only MEHP and MEP remained significant, and sum of low MWP showed a similar trend (p<0.1). As shown in Table 5, the significant associations of phthalate metabolites with WC decreased in the basis of those with BMI in crude analysis and adjusted model 1.and MEHP and MEP were still positively associated with WC in adjusted model 2, but compared to BMI, the association of MiBP with WC showed a positive trend (p<0.1) and MMP was inversely associated with WC.

As shown in Figure 1, and Tables S6 and S7 in Supporting Information File S1, the associations between urine metabolites and BMI showed some sex- or age-specific differences in adjusted model 1. Most of these associations tended to be stronger in younger children than in older ones and in males than in females although that interaction test suggested that the associations of phthalate metabolites with BMI did not change by age or sex with an exception of MCMHP by age (p = 0.061). Interestingly, both MEHHP and MEOHP were found to be significantly associated with BMI in the 8–11 year age group (Table S6 in Supporting Information File S1), but these associations were not significant in all subjects (Table 4). After adjusted for additional urine phthalate metabolites in adjusted model 2, these significant associations disappeared (data not shown).

thumbnail
Figure 1. Comparison of the associations of urine phthalate metabolites corrected by specific gravity with body mass index between age subgroups (A) or sex (B) by multiple linear regression in 259 school children of Shanghai, China in 2012 (A: adjustment for age in years and sex; B: adjustment for age; regression coefficient: log-transformed BMI/log-transformed metabolite concentration.; Error bar: 95% confidence interval).

https://doi.org/10.1371/journal.pone.0056800.g001

Discussion

On the whole, urine phthalate metabolite concentrations in this study were comparable to those from previous studies conducted in other countries (Table S8 in Supporting Information File S1), although the profiles of urinary phthalate metabolites demonstrated some differences in terms of quantity and species of metabolites. For example, MBzP was detected only in 38.6% of subjects in this study, but detected in almost 100% in the studies of U.S., German, and Spanish children [32], [44], [45]. These differences may primarily reflect regional differences in use of phthalates. Additionally, different race and urine type (first morning urine in this study) may contribute to these differences as well.

As a result of an extensive use of phthalate-related products, human can be exposed to phthalates through almost all possible pathways, such as ingestion, inhalation, dermal contact, and placenta transmission [7]. Generally speaking, ingestion of food is the major exposure route of phthalates, in particular for high molecular weight phthalates, such as DEHP and BBzP [46]. For low molecular weight phthalates, such as DBP, DiBP, and DMP, inhalation of indoor air or dermal contact may be relevant exposure pathways next to food ingestion [47]. For DEP, one of low molecular weight phthalates, the major exposure route may largely be through dermal contact (especially by use of personal-care products) [48]. It was expected that exposure pathway-related variables, such as caloric intake and food composition related to food ingestion, and daily indoor time related to inhalation of indoor air, were positively associated with corresponding phthalate exposures, and then urine phthalate metabolites. Sum of DEHP and MCHP may be proxies for food ingestion-related variables, such as caloric intake and food composition, MMP may be for daily indoor time, and MEP may be for use of personal-care products, and sum of DBP or MiBP may be related to both food ingestion and daily indoor time.

In the present study we found nine urine phthalate metabolites and five sums were positively associated with BMI or WC in Chinese school children when adjusted for individual SG, age, and sex (in adjusted model 1), but most of these significant associations disappeared and some even reversed when further adjusted for other urine phthalate metabolites (in adjusted model 2). Only MEHP and MEP remained significant. Close correlations among phthalate metabolite measurements may be one reason for these association changes (Table S4 in Supporting Information File S1), and these changes may mean other possibilities. One possibility was that some phthalate metabolites might act as confounders. When these metabolites were included as covariables, they corrected the associations of metabolites of interest with BMI or WC. Another possibility was that some phthalates might have modifying effects on BMI or WC [49], [50]. A simple adjustment might diminish true associations of metabolites of interest with BMI or WC. So it was expected that a still significant association might be more reliable after adjustment for other urine phthalate metabolites.

Given the cross-sectional nature of this study, the possibility of reverse causation cannot be ruled out. It is possible that increased overall caloric intake, different diet composition, and longer daily indoor time for obese individuals might increase the risk of phthalate exposures, compared to normal weight ones [51]. These variables related to obesity were also confounders for the relationships between phthalate exposures and obesity. Lack of information on these obesity-related variables was an important limitation for this study. Fortunately, since urine phthalate metabolites could be proxies for these confounders as discussed above, the inclusion of urine phthalate metabolites as covariables in multiple linear regression model (adjusted model 2) corrected potential effects of these confounders unavailable in this study to some extent. Therefore, it should be more valid for the presence of significant associations of MEHP and MEP with BMI or WC after adjustment for age, sex as well as other urine phthalate metabolites.

An underlying assumption in this study and other similar studies is that one single measurement of spot urine could represent a long-term exposure of corresponding participant to no-persistent environmental chemicals [52]. Given short biologic half-lives of phthalate metabolites (12 hours or less), multiple sources and exposure routes of phthalates, varying diet composition and behavior of participants, urine collection time in a day and urine dilution, one single measurement of spot urine is not likely to perfectly represent a long-term exposure [53]. Some researches assessed the prediction of single spot urine for a long-term exposure and reported to have a moderate predictabilityapp:ds:power for long-term exposure and this ability differed by phthalate metabolites [53]. For example, Hauser et al. reported that single spot urine was moderately predictive of each subject’s exposure tertile categorization over three months with the sensitivities ranging from 0.56 to 0.74, and most predictive for MEP and least predictive for MEHP among five metabolites in 11 adult men [39]. Teitelbaum et al. assessed the variability of eight creatinine-corrected phthalate metabolites over six months and found moderate reproducibilities with intraclass correlation coefficients (ICCs) ranging from 0.21 (mono(3-carboxypropyl) phthalate, MCPP) to 0.62 (MBzP) in 35 children 6–10 years of age [54]. Additionally, misclassification due to one measurement of spot urine might result in an attenuated estimate of association of phthalate metabolite exposures with BMI or WC, which meant that significant associations based on spot urine were conservative [51]. In this study, we used first morning urine instead of spot one. Because the urine collection time is similar in a day for first morning urine, it decreases the variations of urine phthalate metabolite concentration. This strengthened the predictive ability of one measurement for long-term exposure than pure spot urine [53].

In our studies, we observed some age- or sex-related differences in the associations of phthalate metabolites with BMI or WC. MEHHP and MEOHP were found to be significantly associated with BMI only in the 8–11 year age group, not in all the subjects. These age- or sex-related differences were also reported in some previous studies [29], [30]. These differences might be related to different status of sex hormone endocrine of participants between sex or age subgroups (8–11 or 12–15 year age group). Puberty usually begins between 10 and 13 years of age, and serum levels of endogenous sex hormone in puberty significantly rise, especially for estradiol in girls and testosterone in boys. The differences in sex hormone function may lead to different susceptibility of participants to sex hormone perturbation caused by phthalate exposure [35]. More accurate measures of sex hormone levels and puberty status are needed to further explore the biologic mechanisms for the interrelationship among sex hormone, phthalate exposures and BMI.

There were some epidemiologic studies supporting the associations of postnatal phthalate exposures and body size measures in children. Using data from the National Health and Nutrition Examination Survey (NHANES) 1999–2002, Hatch et al. found that BMI increased with urinary MEP quartiles and inversely related to MEHP quartiles in girls 12–19 years of age, but no metabolite quartiles were significantly associated with BMI in the 6–11 year age group [29]. Boas et al. observed that some metabolites, such as MEP, MEOHP, MECPP, and sum of DEHP metabo­lites, negatively correlated with BMI in 845 Danish children 4–9 years of age although this study was primarily designed to investigate the association of urine phthalate metabolites and thyroid function, insulin-like factor 1, and growth [24]. Teitelbaum et al. reported no significant associations of urinary nine phthalate metabolites at baseline with BMI or WC measured one year later in 387 Hispanic and Black children 6–8 years of age in New York City; however, after divided into normal weight and overweight groups, the data showed significant associations of MEP and sum of low molecular-weight (MEP, MBP, and MiBP) with BMI and WC in overweight children [30]. Our study results were comparable to those from Teitelbaum et al. and Hatch et al., except for an inverse association of MEHP with BMI reported by Hatch et al., but distinctly different from those by Boas et al.

There are some plausible explanations for these discrepancies. Firstly, creatinine correction for urine phthalate metabolite concentrations can weaken or even reverse the association due to the positive associations of urine phthalate metabolite concentrations and BMI with urinary creatinine [30], [40]. In this circumference, SG correction may be a better choice [55]. The studies conducted by Teitelbaum et al. and Hatch et al. used creatinine to correct the urine dilution. Secondly, some studies including the study by Boas et al. reported that phthalates, such as DEHP and DiNP, might suppress thyroid hormone levels and insulin-like growth factor I (IGF-I) in serum [23], [24]. This suppression would decrease the growth of children, and therefore BMI values. When this trend exceeded the adipogenesis effect caused by phthalates, phthalate exposures would be negatively associated with BMI. Because younger children experience more rapid growth, their growth is more likely to be affected by external disturbance. This might explain the negative associations of phthalate exposures with BMI in the study of younger study population by Boas et al. [24]. There are some mice studies that have shown negative associations between prenatal phthalate exposures, such as DBP, BBP, and DEHP, and offspring weight [56], [57]. If these explanations are reasonable, DEP, DEHP, DBP, and DiBP are main phthalates associated with BMI or WC.

Conclusions

In this study we found nine urinary phthalate metabolites and three sums were positively associated with BMI or WC in Chinese school children, and it was more evident for MEHP and MEP. Some associations showed sex and/or age related difference. Given the cross-sectional nature of this study and lack of some important obesity-related variables, the associations discovered here need to be further investigated.

Supporting Information

Supporting Information File S1.

Table S1 BMI-based criteria by age and sex for screening of overweight and obese school children proposed by Working Group on Obesity in China (WGOC) in 2005. Table S2 Settings for multiple reaction monitoring (MRM) mode. Table S3 Results of quality assurance/quality control. Table S4 Pearson correlation coefficients among phthalate metabolites uncorrected by specific gravity. Table S5 Calculation of sampling weights of each selected school (sampling weights = number of source population/number of participants). Table S6 Association of urinary phthalate metabolite concentrations corrected by specific gravity with body mass index by age groups after adjustment for age in years and sex in 259 school children of Shanghai, China in 2011–2012. Table S7 Association of urinary phthalate metabolite concentrations corrected by specific gravity with body mass index by sex after adjustment for age in 259 school children of Shanghai, China in 2011–2012. Table S8 Comparison of urinary phthalate metabolite concentrations (geometric mean, ng/mL) (detection frequency, %) between this study and other several populations.

https://doi.org/10.1371/journal.pone.0056800.s001

(DOC)

Acknowledgments

We thank Dr. H. Kan and B. Chen, School of Public Health, Fudan University for the comments on all or part of our manuscript; Dr J. Luo, School of Public Health, Fudan University for the help on statistics.

Author Contributions

Critical revision of the manuscript for important intellectual content: YC QJ. Conceived and designed the experiments: YZ YC QJ. Performed the experiments: CT JW YH HW YZ. Analyzed the data: HW YZ YC QJ. Contributed reagents/materials/analysis tools: CT JW YH. Wrote the paper: HW YZ YC.

References

  1. 1. Guo Y, Alomirah H, Cho HS, Minh TB, Mohd MA, et al. (2011) Occurrence of phthalate metabolites in human urine from several Asian countries. Environ Sci Technol 45: 3138–3144.
  2. 2. Wittassek M, Koch HM, Angerer J, Bruning T (2011) Assessing exposure to phthalates - the human biomonitoring approach. Mol Nutr Food Res 55: 7–31.
  3. 3. Zhou Y, Wang H, Chen Y, Jiang Q (2011) Environmental and food contamination with plasticisers in China. Lancet 378: E4.
  4. 4. Latini G (2005) Monitoring phthalate exposure in humans. Clin Chim Acta 361: 20–29.
  5. 5. Yen T, Lin-Tan D, Lin J (2011) Food safety involving ingestion of foods and beverages prepared with phthalate-plasticizer-containing clouding agents. J Formos Med Assoc 110: 671–684.
  6. 6. Guo Y, Zhang Z, Liu L, Li Y, Ren N, et al.. (2012) Occurrence and profiles of phthalates in foodstuffs from China and their Implications for human exposure. J Agric Food Chem.
  7. 7. Swan SH (2008) Environmental phthalate exposure in relation to reproductive outcomes and other health endpoints in humans. Environ Res 108: 177–184.
  8. 8. Hines EP, Calafat AM, Silva MJ, Mendola P, Fenton SE (2009) Concentrations of phthalate metabolites in milk, urine, saliva, and Serum of lactating North Carolina women. Environ Health Perspect 117: 86–92.
  9. 9. Guo Y, Wu Q, Kannan K (2011) Phthalate metabolites in urine from China, and implications for human exposures. Environ Int 37: 893–898.
  10. 10. Meeker JD, Sathyanarayana S, Swan SH (2009) Phthalates and other additives in plastics: human exposure and associated health outcomes. Philos Trans R Soc Lond B Biol Sci 364: 2097–2113.
  11. 11. Gruen F, Blumberg B (2009) Endocrine disrupters as obesogens. Mol Cell Endocrinol 304: 19–29.
  12. 12. Holtcamp W (2012) Obesogens: an environmental link to obesity. Environ Health Perspect 120: a62–a68.
  13. 13. Thayer KA, Heindel JJ, Bucher JR, Gallo MA (2012) Role of environmental chemicals in diabetes and obesity: a national toxicology program workshop review. Environ Health Persp 120: 779–789.
  14. 14. Hatch EE, Nelson JW, Stahlhut RW, Webster TF (2010) Association of endocrine disruptors and obesity: perspectives from epidemiological studies. Int J Androl 33: 324–331.
  15. 15. Sargis RM, Johnson DN, Choudhury RA, Brady MJ (2010) Environmental endocrine disruptors promote adipogenesis in the 3T3-L1 cell line through glucocorticoid receptor activation. Obesity 18: 1283–1288.
  16. 16. Casals-Casas C, Desvergne B (2011) Endocrine disruptors: from endocrine to metabolic disruption. Annu Rev Physiol 73: 135–162.
  17. 17. Desvergne B, Feige JN, Casals-Casas C (2009) PPAR-mediated activity of phthalates: A link to the obesity epidemic? Mol Cell Endocrinol 304: 43–48.
  18. 18. Pasquali R (2006) Obesity and androgens: facts and perspectives. Fertil Steril 85: 1319–1340.
  19. 19. Lovekamp TN, Davis BJ (2001) Mono-(2-ethylhexyl) phthalate suppresses aromatase transcript levels and estradiol production in cultured rat granulosa cells. Toxicol Appl Pharmacol 172: 217–224.
  20. 20. Desdoits-Lethimonier C, Albert O, Le Bizec B, Perdu E, Zalko D, et al. (2012) Human testis steroidogenesis is inhibited by phthalates. Hum Reprod 27: 1451–1459.
  21. 21. Takeuchi S, Iida M, Kobayashi S, Jin K, Matsuda T, et al. (2005) Differential effects of phthalate esters on transcriptional activities via human estrogen receptors alpha and beta, and androgen receptor. Toxicology 210: 223–233.
  22. 22. Shen O, Wu W, Du G, Liu R, Yu L, et al. (2011) Thyroid disruption by Di-n-butyl phthalate (DBP) and mono-n-butyl phthalate (MBP) in Xenopus laevis. PLoS One 6: e19159.
  23. 23. Meeker JD, Calafat AM, Hauser R (2007) Di(2-ethylhexyl) phthalate metabolites may alter thyroid hormone levels in men. Environ Health Perspect 115: 1029–1034.
  24. 24. Boas M, Frederiksen H, Feldt-Rasmussen U, Skakkebaek NE, Hegedus L, et al. (2010) Childhood exposure to phthalates: associations with thyroid function, insulin-like growth factor I, and growth. Environ Health Perspect 118: 1458–1464.
  25. 25. Kitahara CM, Platz EA, Ladenson PW, Mondul AM, Menke A, et al. (2012) Body fatness and markers of thyroid function among U.S. men and women. PLoS One 7: e34979.
  26. 26. Bility MT, Thompson JT, Mckee RH, David RM, Butala JH, et al. (2004) Activation of mouse and human peroxisome proliferator-activated receptors (PPARs) by phthalate monoesters. Toxicol Sci 82: 170–182.
  27. 27. Feige JN, Gelman L, Rossi D, Zoete V, Metivier R, et al. (2007) The endocrine disruptor monoethyl-hexyl-phthalate is a selective peroxisome proliferator-activated receptor gamma modulator that promotes adipogenesis. J Biol Chem 282: 19152–19166.
  28. 28. Stahlhut RW, van Wijngaarden E, Dye TD, Cook S, Swan SH (2007) Concentrations of urinary phthalate metabolites are associated with increased waist circumference and insulin resistance in adult U.S. males. Environ Health Perspect 115: 876–882.
  29. 29. Hatch EE, Nelson JW, Qureshi MM, Weinberg J, Moore LL, et al. (2008) Association of urinary phthalate metabolite concentrations with body mass index and waist circumference: a cross-sectional study of NHANES data, 1999–2002. Environ Health 7: 27.
  30. 30. Teitelbaum SL, Mervish N, Moshier EL, Vangeepuram N, Galvez MP, et al. (2012) Associations between phthalate metabolite urinary concentrations and body size measures in New York City children. Environ Res 112: 186–193.
  31. 31. Silva MJ, Barr DB, Reidy JA, Malek NA, Hodge CC, et al. (2004) Urinary levels of seven phthalate metabolites in the U.S. population from the National Health and Nutrition Examination Survey (NHANES) 1999–2000. Environ Health Perspect 112: 331–338.
  32. 32. Becker K, Goen T, Seiwert M, Conrad A, Pick-Fuss H, et al. (2009) GerES IV: phthalate metabolites and bisphenol A in urine of German children. Int J Hyg Environ Health 212: 685–692.
  33. 33. Pietrobelli A, Faith MS, Allison DB, Gallagher D, Chiumello G, et al. (1998) Body mass index as a measure of adiposity among children and adolescents: A validation study. J Pediatr 132: 204–210.
  34. 34. Ji CY, Yt SR, Ma GS, Ma J, He ZH, et al. (2010) Waist circumference distribution of Chinese school-age children and adolescents. Biomed Environ Sci 23: 12–20.
  35. 35. Wang HX, Zhou Y, Tang CX, Wu JG, Chen Y, et al. (2012) Association between bisphenol A exposure and body mass index in Chinese school children: a cross-sectional study. Environ Health 11: 79.
  36. 36. Cameron N (1978) The methods of axiological anthropometry. In: Falkner F, Tanner JM, editors. Human Growth, Vol 2, 1st ed. New York: Plenum. 35–90.
  37. 37. Ji CY (2005) Report on childhood obesity in China (1) - Body mass index reference for screening overweight and obesity in Chinese school-age children. Biomed Environ Sci 18: 390–400.
  38. 38. Hornung RW, Reed LD (1990) Estimation of average concentration in the presence of nondetectable Values. Applied Occupational and Environmental Hygiene 5: 46–51.
  39. 39. Hauser R, Meeker JD, Park S, Silva MJ, Calafat AM (2004) Temporal variability of urinary phthalate metabolite levels in men of reproductive age. Environ Health Perspect 112: 1734–1740.
  40. 40. Barr DB, Wilder LC, Caudill SP, Gonzalez AJ, Needham LL, et al. (2005) Urinary creatinine concentrations in the US population: Implications for urinary biologic monitoring measurements. Environ Health Persp 113: 192–200.
  41. 41. Boeniger MF, Lowry LK, Rosenberg J (1993) Interpretation of urine results used to assess chemical exposure with emphasis on creatinine adjustments: a review. Am Ind Hyg Assoc J 54: 615–627.
  42. 42. Winship C, Radbill L (1994) Sampling weights and regression analysis. Sociol Methods Res 23: 230–257.
  43. 43. Korn EL, Graubard BI (1991) Epidemiologic studies utilizing surveys-accounting for the sampling design. Am J Public Health 81: 1166–1173.
  44. 44. Casas L, Fernandez MF, Llop S, Guxens M, Ballester F, et al. (2011) Urinary concentrations of phthalates and phenols in a population of Spanish pregnant women and children. Environ Int 37: 858–866.
  45. 45. CDC (Centers for Disease Control and Prevention) (2012) Fourth national report on human exposure to environmental chemicals. Available: http://www.cdc.gov/exposurereport/pdf/FourthReport_UpdatedTables_Feb2012.pdf. Accessed 2012 july 1.
  46. 46. Clark K, Cousins I, Mac Kay D (2003) Assessment of critical exposure pathways.In: Staples C, ed. The Handbook of Environmental Chemistry, Vol 3Q: Phthalate Esters. New York: Springer. 227–262.
  47. 47. Wormuth M, Scheringer M, Vollenweider M, Hungerbuhler K (2006) What are the sources of exposure to eight frequently used phthalic acid esters in Europeans? Risk Anal 26: 803–824.
  48. 48. Duty SM, Ackerman RM, Calafat AM, Hauser R (2005) Personal care product use predicts urinary concentrations of some phthalate monoesters. Environ Health Perspect 113: 1530–1535.
  49. 49. Rajapakse N, Silva E, Kortenkamp A (2002) Combining xenoestrogens at levels below individual No-observed-effect concentrations dramatically enhances steroid hormone action. Environ Health Persp 110: 917–921.
  50. 50. Rider CV, Furr J, Wilson VS, Jr. Gray LE (2008) A mixture of seven antiandrogens induces reproductive malformations in rats. Int J Androl 31: 249–262.
  51. 51. Trasande L, Attina TM, Blustein J (2012) Association between urinary bisphenol A concentration and obesity prevalence in children and adolescents. JAMA 308: 1113–1121.
  52. 52. Aylward LL, Kirman CR, Adgate JL, Mckenzie LM, Hays SM (2012) Interpreting variability in population biomonitoring data: role of elimination kinetics. J Expo Sci Environ Epidemiol 22: 398–408.
  53. 53. Preau JJ, Wong LY, Silva MJ, Needham LL, Calafat AM (2010) Variability over 1 week in the urinary concentrations of metabolites of diethyl phthalate and di(2-ethylhexyl) phthalate among eight adults: an observational study. Environ Health Perspect 118: 1748–1754.
  54. 54. Teitelbaum SL, Britton JA, Calafat AM, Ye X, Silva MJ, et al. (2008) Temporal variability in urinary concentrations of phthalate metabolites, phytoestrogens and phenols among minority children in the United States. Environ Res 106: 257–269.
  55. 55. Braun JM, Smith KW, Williams PL, Calafat AM, Berry K, et al. (2012) Variability of urinary phthalate metabolite and bisphenol A concentrations before and during pregnancy. Environ Health Perspect 120: 739–745.
  56. 56. Tanaka T (2005) Reproductive and neurobehavioural effects of bis(2-ethylhexyl) phthalate (DEHP) in a cross-mating toxicity study of mice. Food Chem Toxicol 43: 581–589.
  57. 57. Tyl RW, Myers CB, Marr MC, Fail PA, Seely JC, et al. (2004) Reproductive toxicity evaluation of dietary butyl benzyl phthalate (BBP) in rats. Reprod Toxicol 18: 241–264.