Low-Level Environmental Phthalate Exposure ... - ACS Publications

May 3, 2016 - The study cohort consisted of the volunteers enrolled in NJMU Fertility Study, which was designed to investigate the association between...
0 downloads 0 Views 1MB Size
Subscriber access provided by University of Sussex Library

Article

Low-level environmental phthalate exposure associates with urine metabolome alteration in a Chinese male cohort Jie Zhang, Liangpo Liu, Xiaofei Wang, Qingyu Huang, Meiping Tian, and Heqing Shen Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.6b00034 • Publication Date (Web): 03 May 2016 Downloaded from http://pubs.acs.org on May 8, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Environmental Science & Technology is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 31

Environmental Science & Technology

1

Low-level environmental phthalate exposure associates with urine metabolome

2

alteration in a Chinese male cohort

3

Jie Zhang*, Liangpo Liu, Xiaofei Wang, Qingyu Huang, Meiping Tian, Heqing Shen*

4

Key Lab of Urban Environment and Health, Institute of Urban Environment, Chinese Academy of

5

Sciences, China

6

*To whom correspondence may be addressed:

7

Jie Zhang, Institute of Urban Environment, Chinese Academy of Sciences, 1799 Jimei Road,

8

Xiamen, 361021, China; Telephone/Fax: (86)-592-6190523; E-mail: [email protected]

9

Heqing Shen, Institute of Urban Environment, Chinese Academy of Sciences, 1799 Jimei Road,

10

Xiamen, 361021, China; Telephone/Fax: (86)-592-6190771; E-mail: [email protected]

11 12

ACS Paragon Plus Environment

Environmental Science & Technology

Page 2 of 31

13

Abstract

14

The general population is exposed to phthalates through various sources and routes. Integration of

15

omics data and epidemiological data is a key step towards directly linking phthalate

16

bio-monitoring data with biological response. Urine metabolomics is a powerful tool to identify

17

exposure biomarkers and delineate the modes of action of environmental stressors. The objectives

18

of this study are to investigate the association between low-level environmental phthalate

19

exposure and urine metabolome alteration in male population, and to unveil the metabolic

20

pathways involved in the mechanisms of phthalate toxicity. In this retrospective cross-sectional

21

study, we studied the urine metabolomic profiles of 364 male subjects exposed to low-level

22

environmental phthalates. Di(2-ethylhexyl) phthalate (DEHP) and dibutyl phthalate (DBP) are the

23

most widely used phthalates. ΣDEHP and MBP (the major metabolite of DBP) were associated

24

with significant alteration of global urine metabolome in the male population. We observed

25

significant increase in the levels of acetylneuraminic acid, carnitine C8:1, carnitine C18:0, cystine,

26

phenylglycine, phenylpyruvic acid and glutamylphenylalanine; and meanwhile, decrease in the

27

levels

28

methylglutaconic acid, hydroxyl-PEG2 and keto-PGE2 in high exposure group. The observations

29

indicated that low-level environmental phthalate exposure associated with increased oxidative

30

stress and fatty acid oxidation and decreased prostaglandin metabolism. Urea cycle, tryptophan

31

and phenylalanine metabolism disruption was also observed. The urine metabolome disruption

32

effects associated with ΣDEHP and MEP were similar, but not identical. The multi-biomarker

33

models presented AUC values of 0.845 and 0.834 for ΣDEHP and MEP, respectively. The

34

predictive accuracy rates of established models were 81% for ΣDEHP and 73% for MEP. Our

of

carnitine

C16:2,

diacetylspermine,

alanine,

taurine,

ACS Paragon Plus Environment

tryptophan,

ornithine,

Page 3 of 31

Environmental Science & Technology

35

results suggest that low-level environmental phthalate exposure associates with urine metabolome

36

disruption in male population, providing new insight into the early molecular events of phthalate

37

exposure.

38

ACS Paragon Plus Environment

Environmental Science & Technology

39

Introduction

40

The ubiquitousness of phthalates results in widespread exposure to these chemicals among the

41

general population. Phthalates has been tightly linked to adverse health outcomes, such as diabetes,

42

allergy symptoms and reproductive abnormalities.1-4 Phthalates are rapidly converted into

43

mono-phthalates in human body,

44

bio-monitoring indicators. However, mono-phthalates are not enough to characterize adverse

45

effects due to their irrelevance with molecular response to the exposure. Mechanisms of response

46

and biomarkers of susceptibility to phthalates are interrelated, thus the development of new

47

strategies- especially the identification of early effect biomarkers using omics approaches (e.g.

48

metabolomics) is of particular importance.6 Furthermore, the integration of omics data with

49

mechanistic and epidemiological data is a key step towards linking exposure biomarkers and

50

susceptibility to disease mechanisms and outcomes.7

51

Metabolomics studies from our group

52

disrupted many important metabolic pathways in rodents and other experimental models. Animal

53

models used in toxicological studies have been criticized due to the lack of similarities to human

54

exposure conditions, in particular through the use of high doses and short term exposure. The best

55

strategy for identification of metabolic disturbances in environmental health is to directly target on

56

human population6. Metabolomics offers three distinguished advantages in studying human–

57

environment interactions at molecular level: (1) metabolic fingerprint represents actual functional

58

status of the organism, which can be mechanistically related to organism phenotype; (2)

59

metabolomics focuses more on questions than hypotheses, thus it can discover unexpected

60

associations between metabolite responses and environmental exposure 8; (3) metabolic changes

5

and thus urinary mono-phthalates are commonly used as

9

and others

10-14

have demonstrated phthalate exposure

ACS Paragon Plus Environment

Page 4 of 31

Page 5 of 31

Environmental Science & Technology

61

can be observed in biological fluids (e.g. urine, blood), making it possible to directly identify

62

potential biomarkers that may indicate the early effects of environmental exposure in humans.

63

Recently, researchers began to apply environmental metabolomics to epidemiological studies.

64

Differential urine metabolomes have been characterized for male welders vs. office workers

65

volunteers living near a closed zinc smelter of cadmium pollution16, pregnant women exposed to

66

low-dose/complex pesticides17, and diabetic/non-diabetic population exposed to inorganic arsenic

67

through drinking water18. These studies investigated the metabolome alteration of the populations

68

with occupational exposure

69

known about low-level environmental exposure in general population.

70

In this study, we aimed at answering the following questions: What is the influence of low-level

71

environmental phthalate exposure on the urine metabolome in general male population? Which

72

metabolic pathways are involved in the mechanisms of phthalate toxicity? We performed liquid

73

chromatography - mass spectrometry based urine metabolomic analyses of 364 Chinese men in

74

order to identify discriminant metabolites, and we further investigated disrupted metabolic

75

pathways associated with the early effects of phthalates exposure.

15

or the populations living in highly polluted areas,

16-18

15

,

but little is

76 77

Methods

78

Subject demographics and sample collection

79

The retrospective, cross-sectional study was approved by the institutional ethics committee and

80

conducted under the Helsinki Declaration. Three hundred and sixty four Chinese adult men

81

(ethnically Han) were enrolled from affiliated hospitals of Nanjing Medical University (NJMU).

82

They were the volunteers of NJMU Fertility Study. Informed consent was obtained from all

ACS Paragon Plus Environment

Environmental Science & Technology

83

participants. The data including age, weight, height, smoking and alcohol drinking status,

84

education level, annual income, and occupation were collected by questionnaire. The subjects with

85

diabetes mellitus were excluded to avoid any ambiguity. Spot morning urine samples were

86

collected for mono-phthalate determination and metabolic profile acquisition. After collection, the

87

samples were transported in ice to the lab in 24 hours, and then stored in -80°C prior to analysis.

88

Urinary phthalate metabolite analysis

89

Seven major phthalate metabolites were determined using liquid chromatography - mass

90

spectrometry, including monomethyl phthalate (MMP), monoethyl phthalate (MEP), monobutyl

91

phthalate (MBP), monobenzyl phthalate (MBzP), mono-2-ethylhexyl phthalate (MEHP),

92

mono-2-ethyl-5-oxohexyl phthalate (MEOHP) and mono-(2-ethyl-5-hydroxylhexyl) phthalate

93

(MEHHP). The details of sample preparation, standard solutions and instrumental analysis have

94

been described previously

95

β-glucuronididase, purified by solid-phase extraction, and analyzed by isotope dilution

96

high-performance liquid chromatography - tandem mass spectrometry. The limits of detection

97

(LODs) of MMP, MEP, MBP, MBzP, MEHP, MEOHP and MEHHP were 0.08, 0.10, 0.10, 0.14,

98

0.07, 0.06 and 0.10 µg/L, respectively. The relative standard deviation (RSD) values of peak area

99

for these phthalate metabolites ranged from 2.2 to 13.9%. Urinary concentrations of phthalate

100

metabolites were expressed as creatinine-adjusted concentrations. Di(2-ethylhexyl) phthalate

101

(DEHP) and dibutyl phthalate (DBP) are the most widely used phthalates. MBP is the major

102

metabolite of DBP. ΣDEHP was calculated based on mass sum of major metabolite concentrations

103

(MEHP, MEHHP and MEOHP) 20.

104

Urine metabolome analysis

19

. Briefly, urinary phthalate metabolites were deconjugated using

ACS Paragon Plus Environment

Page 6 of 31

Page 7 of 31

Environmental Science & Technology

105

Details of sample preparation, metabolic profiling acquisition, data processing, biomarker

106

screening and identification, quality control procedures were described in Supporting Information.

107

Briefly, the urine metabolic profiles were acquired using a liquid chromatography/time-of-flight

108

mass spectrometer, and then processed with Profile Analysis 2.0 to obtain a feature table. The

109

table was Pareto-scaled and introduced to SIMCA-P v11.5 software (Umetrics, Sweden) for

110

multivariate statistical analysis. Principal component analysis (PCA) was performed to cluster the

111

samples, and identified the outliers. The outliners were the samples far away from the cluster

112

center in PCA score plot. The samples were possibly contaminated during the storage and

113

preparation processes. After careful inspection of raw chromatograms, the outliers (n=2) were

114

removed from the dataset. ΣDEHP, MBP and MEP concentration was classified into tertiles.

115

Accordingly, the dataset were separately categorized into three groups (1st tertile, 2nd tertile and 3rd

116

tertile), and defined as low, medium and high exposure groups, respectively. MBzP was detected

117

in 52.7% of the subjects, thus the dataset was divided into two groups (detected vs. undetected).

118

Phthalate exposure-oriented PLS-DA (partial least squares discriminant analysis) models were

119

established, in which internal levels of phthalates were used as classifiers. The biomarkers were

120

screened from the PLS-DA models according to the following criteria: (1) VIP scores > 2; (2)

121

jack-knifing confidence interval > 0; (3) p-value between low and high exposure groups < 0.05;

122

and (4) the features significantly (p < 0.05) correlated with phthalate concentrations after adjusted

123

by age, BMI, smoke and alcohol drinking status; (5) the features significantly (p < 0.05) correlated

124

with semen volume, concentration, motility, progression, and motion parameters were excluded.

125

Statistical analysis

126

All the statistical analysis was performed using SPSS 18 (SPSS Inc.). The nonparametric

ACS Paragon Plus Environment

Environmental Science & Technology

127

Mann-Whitney test was used to evaluate the significant difference of each potential biomarker

128

between the groups. Spearman correlation was performed to investigate the association between

129

demographic factors, phthalates and biomarkers. Partial correlation analysis was performed to

130

investigate the association between the biomarkers and phthalate exposure after adjustment of age,

131

BMI, smoking and alcohol drinking status. P-values of < 0.05 were considered statistically

132

significant.

133

Receiver operating characteristic (ROC) curves

134

Classical univariate ROC analyses of potential biomarkers were performed using SPSS 18 (SPSS

135

Inc.). Multivariate analyses of combinational biomarker patterns were performed using online

136

ROCCET (ROC Curve Explorer & Tester) software (http://www.roccet.ca/ROCCET/).

137 138

Results

139

Population characteristics

140

The demographic characteristics of study population are listed in Table 1. The participants are

141

male residents of Nanjing city aging from 19 to 44 years. BMI exhibited a wide range of 16.5–

142

32.9 kg/m2 (median 23.7 kg/m2). Half of the participants never consume tobacco and alcohol.

143

Urinary phthalate levels

144

Urinary phthalate levels are listed in Table 2. MEP, MBP, MEHP, MEHHP, and MEOHP were

145

detected in all the samples. MBzP were detected in 52.7% of the subjects. MMP were only

146

detected in 3.6% of the subjects, thus it was not included in metabolomics analysis. We detected

147

much higher concentrations of MBP, MEP and MMP than other phthalate metabolites. Statistical

148

analysis showed ΣDEHP were highly correlated with its major metabolites (r= 0.90, 0.99 and 0.98

ACS Paragon Plus Environment

Page 8 of 31

Page 9 of 31

Environmental Science & Technology

149

for MEHP, MEHHP and MEOHP, respectively) (see SI Table S1). In addition, we also observed a

150

moderate correlation between ΣDEHP and MBP (r=0.55) (see SI Table S1). Correlation analysis

151

showed a slight positive correlation (r ≤ 0.22) between phthalate metabolites and BMI (except

152

MBzP and MEHP), but no significant correlation was observed between phthalates and age. The

153

tertiles of ΣDEHP and MBP (see SI Table S2) also showed no significant difference in education

154

level, annual income, occupation, smoking and alcohol drinking status. Moreover, no significant

155

difference of age, BMI and phthalate levels was observed between smokers and non-smokers, and

156

neither was observed between alcohol drinkers and non-alcohol drinkers (see SI Table S3 and S4).

157

Urine metabolome

158

PCA and PLS-DA models were established to investigate the metabolic disruptions (see SI Table

159

S5). In PCA scoring plots, low exposure group was mostly separated from high exposure group,

160

but medium exposure group could not be separated from low and high exposure groups (see

161

Figure 1 and SI Figure S1). We also investigated the influence of smoking and alcohol drinking on

162

urine metabolome. The developed PCA models could not provide differentiation for non-smokers

163

vs. smokers and non-alcohol drinkers vs. alcohol drinkers (see SI Figure S1). In order to further

164

differentiate the groups and screen potential biomarkers, PLS-DA analysis was separately

165

performed for phthalates (low exposure group vs. high exposure group) and smoking/alcohol

166

drinking status. The strict criteria of permutation tests (see SI Figure S2) were used to validate the

167

models: (1) the R2 and Q2 values in the permutation test were lower than the original ones; and (2)

168

the Y-axis intercept for Q2 is below zero. Finally, robust PLS-DA models were only established

169

for MBP and ΣDEHP. Good separations of metabolic profiles were observed between low and

170

high exposure groups in MBP and ΣDEHP models (Figure 1). No satisfactory PLS-DA models

ACS Paragon Plus Environment

Environmental Science & Technology

171

were generated for MEP, MBzP and smoking/alcohol drinking status (see SI Table S5), thus these

172

factors were not examined further in analysis.

173

Biomarker screening was conducted according to the protocol described in experimental section.

174

The study cohort consisted of the volunteers enrolled in NJMU Fertility Study, which was

175

designed to investigate the association between environmental exposure and semen quality. In this

176

study, we aimed to investigate the association between metabolome alteration and phthalate

177

exposure, and identify the disrupted metabolic pathways associated with early effects of phthalates

178

exposure. Since male infertility status may alter the metabolome

179

correlated with semen parameters were removed from the biomarker candidate list (n=2 for MBP,

180

n=1 for ΣDEHP). These removed metabolites have little influence on the overall findings of this

181

study. Finally, twelve biomarkers were selected from MBP model, and five biomarkers were

182

selected from ΣDEHP model (Table 3). These tentatively identified biomarkers were key nodes of

183

the early effects associated with phthalate exposure in male population. The biomarkers were

184

closely involved in amino acid metabolism, mitochondrial beta-oxidation of fatty acids, and

185

prostaglandin metabolism. Except for acetylneuraminic acid and phenylglycine, the biomarkers

186

showed similar trend change in both MBP and ΣDEHP models. As shown in Table 3, in high

187

exposure group, five amino acids and their related metabolites (i.e. alanine, taurine, tryptophan,

188

ornithine, and methylglutaconic acid) significantly decreased, while other amino acid derivatives

189

(i.e. cystine, phenylglycine, phenylpyruvic acid and glutamylphenylalanine) significantly

190

increased. Carnitines play a key role in mitochondrial beta-oxidation of fatty acids. Carnitine C8:1

191

and carnitine C18:0 increased after phthalate exposure. However, another long chain carnitine -

192

carnitine C16:2 exhibited significant decrease in high MEP and ΣDEHP group. Two

21

, the metabolites significantly

ACS Paragon Plus Environment

Page 10 of 31

Page 11 of 31

Environmental Science & Technology

193

prostaglandins – hydroxyl-PGE2 and keto-PGE2 dramatically decreased in high MBP or ΣDEHP

194

group. The change of acetylneuraminic acid was observed only in high MBP group, while the

195

decrease of diacetylspermine was noted in both groups. Most biomarkers showed significant

196

dose-dependent changes in three groups (see SI Figure S3).

197

ROC analysis

198

The closer the AUC value approaches to 1, the better the model provides diagnostic performance.

199

For MBP model, 12 biomarkers gave AUC values between 0.7 and 0.9 (Table 3), indicating their

200

moderate discrimination ability. For ΣDEHP model, 10 biomarkers gave AUC values between 0.7

201

and 0.8 (Table 3). Multi-biomarker models were further established for MBP and ΣDEHP.

202

Combinations of top ten biomarkers turned out to be the best indicators for both MBP and ΣDEHP

203

models (Figure 2). Confusion matrix showed the predictive accuracy of the model as the

204

percentage of correctly classified samples in a given class. The predicative accuracy was about 72%

205

for MBP model and 80% for ΣDEHP model (Figure 2).

206 207

Discussion

208

Humans are daily exposed to phthalates through various sources and routes (diet, dust, air, oral

209

and respiratory). Together with the use of population characteristics data, understanding human

210

metabolome responses to phthalate exposure will provide new fingerprinting tools in assessing

211

health risk and understanding their toxicity.

212

The noticeable phthalates are DEHP with more than 75% of its total usage in plastic products,

213

followed by customer and personal care related chemicals i.e., DBP and DEP 22, 23. In this study,

214

we detected higher concentrations of DBP and DEHP metabolites than others. Furthermore,

ACS Paragon Plus Environment

Environmental Science & Technology

Page 12 of 31

215

correlation analysis revealed strong associations between these two phthalates (R2=0.546),

216

indicating they may have common exposure source (mainly through diet) for the investigated

217

population. The levels of urinary phthalate metabolite in the study were significantly lower than

218

occupational exposure

219

Chinese report22, but were lower than those observed in the National Health and Nutrition

220

Examination Survey 4.

221

Metabolome analysis

222

After phthalate determination, we followed a non-targeted metabolic profiling approach to

223

characterize urine metabolome, to further offer insights into metabolic disruption in male

224

population. The metabolome represents the fingerprint of all small, low molecular weight

225

metabolites in a biological cell, tissue, organ or organism, which are the end products of cellular

226

processes. In this study, an environmental exposure-oriented metabolome strategy was used,

227

which was reported in our previous paper regarding arsenic exposure

228

method, the low exposure group was well discriminated from the high exposure group, suggesting

229

that the exposure to DBP and DEHP associated with substantial global metabolic variations in

230

male population. The individual’s metabolism differs from person to person, resulting in

231

metabolome variation. The phthalate exposure-oriented metabolomic results were derived from

232

statistical analysis of hundreds of participants, thus the uncertainty induced by individual’s

233

metabolism variation were reduced to a certain extent. Furthermore, phthalate exposure-oriented

234

PLS-DA analysis and sequent statistical analysis could exclude the metabolic features associated

235

with the variation of individual’s metabolism.

236

Recently, a few study reported cigarette smoking 26, 27 and alcohol consumption 28 induced serum

24

. The data were comparable to those of our previous report19 and other

ACS Paragon Plus Environment

25

. By using PLS-DA

Page 13 of 31

Environmental Science & Technology

237

or blood metabolomic differences in humans. In this study, we did not observe significant urine

238

metabolome alteration. The main reason is the composition of urine is significantly different from

239

serum and blood. Study population difference our study and previous reports 26-28 may be another

240

reason. No significant difference of MBP/ΣDEHP biomarkers was observed between alcohol

241

drinkers and non-alcohol drinkers; neither was observed between non-smokers and smokers. This

242

indicated the biomarkers were not associated with smoking/alcohol drinking status.

243

Totally 16 differential metabolites were putatively identified in MBP/ΣDEHP statistical models.

244

Several important metabolic pathways were associated with phthalate exposure, for example, an

245

increase in fatty acid beta-oxidation (carnitine C8:1 and carnitine C18:0), a decrease in

246

prostaglandin metabolism (hydroxyl-PGE2 and keto-PGE2) and disturbed amino acid metabolism

247

(alanine, taurine, tryptophan, ornithine, cystine and derivatives). The identified metabolites also

248

showed similar trend in both MBP and ΣDEHP models, indicating these two phthalates may exert

249

similar disruption effect on human metabolism network. This finding is consistent with previous

250

rodent-based metabolomics study, where the authors observed similar change of 39 lipid and

251

amino acid metabolites for the DEHP and DBP exposure in male rats. 12

252

Phthalate exposure associates with increased fatty acid oxidation and oxidative stress

253

Carnitines have been used as potential metabolic modulators to modulate mitochondrial

254

metabolism and alleviate oxidative stress

255

C8:1) were significantly increased in the subjects exposed to high level of phthalates, pointing to

256

possible changes in fatty acid β-oxidation associated with phthalate exposure. Altered

257

acylcarnitine levels have also been observed in adipocytes, brain and reproductive organs of

258

phthalate treated mice or rats. 13, 30 Recently, DEHP treatment was reported to elevate carnitine

29

. Main acylcarnitines (carnitine C18:0 and carnitine

ACS Paragon Plus Environment

Environmental Science & Technology

Page 14 of 31

259

acetyltransferase/palmitoyltransferase activities and up-regulate the expression of peroxisomal and

260

mitochondrial β-oxidation enzymes in rats

261

(PPARs) are closely involved in the regulation of carnitine homeostasis and mitochondrial

262

function. Therefore, we propose that PPARs dis-regulation is a potential mechanism by which

263

phthalates can exert their effects on lipid metabolism in male population.

264

Dis-regulated carnitine metabolism observed in our study also connects phthalate exposure with

265

increased oxidative stress. Elevated acylcarnitine expression is a positive response to oxidative

266

stress induced by phthalate exposure. Glutathione is a key antioxidant, preventing reactive oxygen

267

species to damage cellular components. We observed increased levels of cysteine and

268

glutamylphenylalanine which are important precursors and intermediates in glutathione synthesis.

269

These findings support the hypothesis that environmental phthalate exposure causes oxidative

270

stress-related effects in male population.

271

Phthalate exposure associates with decreased prostaglandin metabolism

272

Our data revealed the association between phthalate exposure and prostaglandin pathway, which is

273

tightly involved in sexual development, inflammatory responses and hormone regulation.

274

Prostaglandin E2 (PGE2) is the most biologically active prostaglandin and widely distributed in

275

various tissues. In our study, the PGE2 metabolites hydroxy- and keto- PGE2 were observed to

276

decrease significantly, indicating that phthalate exposure associates with decreased prostaglandin

277

metabolism. This finding is consistent with previous in-vivo and in-vitro studies

278

that phthalates are endocrine-disrupting compounds. Phthalates inhibit the prostaglandin pathway

279

through releasing cytokine, activating PGD and PGE receptors, activating PPAR pathways, and

280

inhibiting the COX-mediated conversion of arachidonic acid to prostaglandins 34, 35.

31,32

. Peroxisome proliferator-activated receptors

ACS Paragon Plus Environment

33, 34

supporting

Page 15 of 31

Environmental Science & Technology

281

Phthalate exposure associates with amino acid metabolism disruption

282

Serum uric acid has been recently reported to have a significant positive correlation with serum

283

DEHP level in humans.

284

ornithine cycle). Combining our observations, we concluded that phthalate exposure associated

285

with disrupted urea cycle. Alanine and ornithine metabolism-related enzymes (e.g. alanine

286

aminotransferase ALT, ornithine decarboxylase ODC) may be important toxic target. Phthalate

287

exposure significantly elevated the ALT level

288

disrupt alanine synthesis and polyamine pathway. The disregulation of polyamine pathway was

289

further supported by decreased diacetylspermine level observed in our study.

290

Tryptophan is another important essential amino acid. Previous reports revealed DBP and DEHP

291

may increase the conversion ratio of tryptophan to niacin through the inhibition of

292

α-amino-β-carboxymuconate-esemialdehyde decarboxylase38. Although niacin level was not

293

measured, the decreased level of tryptophan observed in our study partially supported previous

294

findings. Phenylalanine metabolism was also probably disturbed by phthalates, indicated by

295

increased

296

glutamylphenylalanine. The increased phenylalanine was observed in male rats treated with 3000

297

ppm DEHP for 4 weeks 12, which was consistent with our hypothesis.

298

ROC analysis

299

The identified biomarkers were the indicators of phthalate exposure because they represented the

300

disrupted pathways. The biomarkers distributed in several important pathways, thus the

301

combinational biomarker pattern represented the overall disrupted metabolic network. As expected,

302

combined biomarker pattern provided higher discrimination power than single biomarkers. The

levels

of

36

its

Alanine and ornithine are involved in urea cycle (also known as

two

important

36

and ODC activity

intermediates

-

ACS Paragon Plus Environment

37

, which were capable to

phenylpyruvic

acid

and

Environmental Science & Technology

303

multi-biomarker models (top ten biomarker combination) presented AUC values of 0.834 and

304

0.845 for MBP and ΣDEHP, respectively. Usually, a biomarker with AUC >0.85 is acceptable for

305

most clinical applications. Therefore, the combinational biomarker pattern has great potential in

306

assessing health risk associated with phthalate exposure in the actual environment.

307 308

Our study suggests a tight relationship between phthalate exposure and urine metabolome

309

disruption in male population, and highlights potential mechanisms that deserve further study.

310

However, it has several limitations: (1) Control of cohort quality remains a big challenge for

311

epidemiological metabolome studies. Due to the difficulty in performing full clinic diagnosis for

312

all participants, a small number of the subjects with diseases were possibly included in the study.

313

Although the developed PLS-DA model in combination with correlation analysis could greatly

314

minimize the influence of these factors on metabolome analysis, excluding these participants will

315

enhance the credibility of this study. (2) Individual metabolic differences may impact both urine

316

metabolome and phthalate metabolite concentrations. Although the biomarkers were screened with

317

strict procedure, we could not completely eliminate the uncertainty induced by individual

318

metabolic differences. (3) Phthalates have different modes of toxic action, thus they may have

319

different influence on urine metabolome. So far, more than 20 phthalates have been used in a large

320

variety of products. Major phthalate metabolites have been investigated in this study, but further

321

study is needed to elucidate the impact of unmeasured phthalates. (4) Humans are daily

322

simultaneously exposed to a variety of harmful chemicals, and some of them possibly have similar

323

influences on metabolic pathway. In this study, the identified phthalate biomarkers were

324

completely different from those obtained from arsenic study 25, and the combinational biomarker

ACS Paragon Plus Environment

Page 16 of 31

Page 17 of 31

Environmental Science & Technology

325

pattern provided accurate assessment for phthalate exposure. However, uncertainty might be

326

introduced into the results due to the incapability of characterizing disruption effect of single

327

exposure for all harmful chemicals. Nevertheless, our data still provided direct evidence that

328

phthalate exposure was associated with urinary metabolome disruption in male adult population,

329

which could be useful in population based risk assessments of phthalate exposure.

330 331

In summary, this study demonstrated the association between low-level environmental phthalate

332

exposure and global urine metabolome disruption in male population. Our results also suggested

333

that phthalate exposure associated with increased oxidative stress and fatty acid oxidation,

334

decreased prostaglandin metabolism, and disrupted amino acid metabolism. The multi-biomarker

335

ROC models showed great potential in assessing health risk associated with phthalate exposure in

336

the actual environment.

337 338

Supporting Information Available

339

Spearman correlations among creatinine-adjusted phthalate concentrations and demographic

340

factors (Table S1); distribution of smoking/alcohol drinking status, education level, annual

341

income, profession for tertiles of phthalates (Table S2); comparison of subject characteristics

342

between smokers and non-smokers, and between alcohol drinkers and non- alcohol drinkers

343

(Table S3); comparison of urinary phthalate levels between smokers and non-smokers, and

344

between alcohol drinkers and non- alcohol drinkers (Table S4); established PCA and PLS-DA

345

models by using phthalates or smoking and alcohol drinking history as the classifiers (Table S5); a

346

PCA scores plot of human urine samples and QCs obtained by HPLC/QTOF-MS (Figure S1);

347

score plots of MEP, smoking status, alcohol drinking status based PCA and PLS-DA models

ACS Paragon Plus Environment

Environmental Science & Technology

348

(Figure S2); the 999-time permutation tests conducted for each PLS-DA model (Figure S3);

349

representative dot plots for fold changes of identified differential metabolites from ΣDEHP and

350

MBP models (Figure S4). This material is available free of charge via the Internet at

351

http://pubs.acs.org.

352 353

Acknowledgements

354

This work was financially supported by the National Natural Science Foundation of China

355

(21407143, 21407144 and 21307126), Youth Innovation Promotion Association of CAS

356

(2015246), Natural Science Foundation of Fujian Province (2013J01063) and the Ningbo Science

357

and Technology Fund (2014A610284).

358 359

Running title: Urine metabolome of phthalate exposure

360 361

Competing financial interests: The authors declare no competing interests.

362 363

References

364

(1) Chen, X.; Xu, S.; Tan, T.; Lee, S.T.; Cheng, S.H.; Lee, F.W.; Xu, L.S.J.; Ho K.C. Toxicity and

365

estrogenic endocrine disrupting activity of phthalates and their mixtures. Int J Environ Res Public

366

Health 2014, 11(3), 3156-3168.

367

(2) Latini, G.; Del Vecchio, A.; Massaro, M.; Verrotti, A.; De Felice, C. Phthalate exposure and

368

male infertility. Toxicology 2006, 226(2-3), 90-98.

369

(3) Sun, Q.; Cornelis, M.C.; Townsend, M.K.; Tobias, D.K.; Eliassen, A.H.; Franke, A.A.;

ACS Paragon Plus Environment

Page 18 of 31

Page 19 of 31

Environmental Science & Technology

370

Hauser, R.; Hu F.B. Association of urinary concentrations of bisphenol A and phthalate

371

metabolites with risk of type 2 diabetes: a prospective investigation in the Nurses' Health Study

372

(NHS) and NHSII Cohorts. Environ Health Perspect. 2014, 122(6), 616-623.

373

(4) Hoppin, J.A.; Jaramillo, R.; London, S.J.; Bertelsen, R.J.; Salo, P.M.; Sandler, D.P.; Zeldin,

374

D.C. Phthalate exposure and allergy in the U.S. population: results from NHANES 2005-2006.

375

Environ Health Perspect 2013,121(10),1129-1134.

376

(5) Lorber, M.; Angerer, J.; Koch, H. A simple pharmacokinetic model to characterize exposure of

377

Americans to Di-2-ethylhexyl phthalate. J Expo Sci Env Epid 2009, 20(1), 38-53.

378

(6) Bonvallot, N.; Tremblay-Franco, M.; Chevrier, C.; Canlet, C.; Debrauwer, L.; Cravedi, J.P.;

379

Cordier, S. Potential input from metabolomics for exploring and understanding the links between

380

environment and health. J Toxicol Environ Health B Crit Rev 2014, 17(1), 21-44.

381

(7) Carlin, D.J.; Naujokas, M.F.; Bradham, K.D.; Cowden, J.; Heacock, M.; Henry, H.F.; Lee, J.S.;

382

Thomas, D.J.; Thompson, C.; Tokar, E.J.; Waalkes, M.P.; Birnbaum, L.S.; Suk, W.A. Arsenic and

383

environmental health: state of the science and future research opportunities. Environ Health

384

Perspect 2015. doi:10.1289/ehp.1510209 [Online 20 November 2015].

385

(8) Bundy, J.G.; Davey, P.M.; Viant, M.R. Environmental metabolomics: a critical review and

386

future perspectives. Metabolomics 2009, 5, 3-21.

387

(9) Zhang, J.; Yan, L.J.; Tian, M.P.; Huang, Q.Y.; Peng, S.Y.; Dong, S.J.; Shen, H.Q. The

388

metabonomics of combined dietary exposure to phthalates and polychlorinated biphenyls in mice.

389

J Pharm Biomed Anal 2012, 66, 287-297.

390

(10) Eveillard, A.; Lasserre, F.; de Tayrac, M.; Polizzi, A.; Claus, S.; Canlet, C.; Mselli-Lakhal, L.;

391

Gotardi, G.; Paris, A.; Guillou, H.; Martin, P.G.; Pineau, T. Identification of potential mechanisms

ACS Paragon Plus Environment

Environmental Science & Technology

392

of toxicity after di-(2-ethylhexyl)-phthalate (DEHP) adult exposure in the liver using a systems

393

biology approach. Toxicol Appl Pharmacol 2009, 236, 282-292.

394

(11) Sumner, S.; Snyder, R.; Burgess, J.; Myers, C.; Tyl, R.; Sloan, C.; Fennell, T. Metabolomics

395

in the assessment of chemical-induced reproductive and developmental outcomes using

396

non-invasive biological fluids: application to the study of butylbenzyl phthalate. J Appl Toxicol

397

2009, 29(8), 703-714.

398

(12) van Ravenzwaay, B.; Coelho-Palermo Cunha, G.; Strauss, V.; Wiemer, J.; Leibold, E.; Kamp,

399

H.; Walk, T.; Mellert, W.; Looser, R.; Prokoudine, A.; Fabian, E.; Krennrich, G.; Herold, M. The

400

individual and combined metabolite profiles (metabolomics) of dibutylphthalate and

401

di(2-ethylhexyl)phthalate following a 28-day dietary exposure in rats. Toxicol Lett 2010, 198(2),

402

159-170.

403

(13) Ellero-Simatos, S.; Claus, S.P.; Benelli, C.; Forest, C.; Letourneur, F.; Cagnard, N.; Beaune,

404

P.H.; de Waziers, I. Combined transcriptomic–1H NMR metabonomic study reveals that

405

monoethylhexyl phthalate stimulates adipogenesis and glyceroneogenesis in human adipocytes. J

406

Proteome Res 2011, 10(12), 5493-5502.

407

(14) Xia, H.; Chi, Y.; Qi, X.; Su, M.; Cao, Y.; Song, P.; Li, X.; Chen, T.; Zhao, A.; Zhang Y.; Cao,

408

Y.; Ma, X.; Jia, W. Metabolomic evaluation of di-n-butyl phthalate-induced teratogenesis in mice.

409

Metabolomics 2011, 7, 559-571.

410

(15) Wang, K.C.; Kuo, C.H.; Tian, T.F.; Tsai, M.H.; Chiung, Y.M.; Hsiech, C.M.; Tsai, S.J.; Wang,

411

S.Y.; Tsai, D.M.; Huang, C.C.; Tseng, Y.J. Metabolomic characterization of laborers exposed to

412

welding fumes. Chem Res Toxicol 2012, 25(3), 676-686.

413

(16) Ellis, J.K.; Athersuch, T.J.; Thomas, L.D.; Teichert, F.; Pérez-Trujillo, M.; Svendsen, C.;

ACS Paragon Plus Environment

Page 20 of 31

Page 21 of 31

Environmental Science & Technology

414

Spurgeon, D.J.; Singh, R.; Järup, L.; Bundy, J.G.; Keun, H.C. Metabolic profiling detects early

415

effects of environmental and lifestyle exposure to cadmium in a human population. BMC

416

Medicine 2012, 10, 61.

417

(17) Bonvallot, N.; Tremblay-Franco, M.; Chevrier, C.; Canlet, C.; Warembourg, C.; Cravedi, J.P.;

418

Cordier, S. Metabolomics tools for describing complex pesticide exposure in pregnant women in

419

Brittany (France). PLoS One 2013, 8(5), e64433.

420

(18) Martin, E.; González-Horta, C.; Rager, J.; Bailey, K.A.; Sánchez-Ramírez, B.;

421

Ballinas-Casarrubias, L.; Ishida, M.C.; Gutiérrez-Torres, D.S.; Hernández Cerón, R.; Viniegra

422

Morales, D.; Baeza Terrazas, F.A.; Saunders, R.J.; Drobná, Z.; Mendez, M.A.; Buse, J.B.; Loomis,

423

D.; Jia, W.; García-Vargas, G.G.; Del Razo, L.M.; Stýblo, M.; Fry, R. Metabolomic characteristics

424

of arsenic-associated diabetes in a prospective cohort in Chihuahua, Mexico. Toxicol Sci. 2015,

425

144(2), 338-346.

426

(19) Liu, L.P.; Bao, H.Q.; Liu, F.; Zhang, J.; Shen, H.Q. Phthalates exposure of Chinese

427

reproductive age couples and its effect on male semen quality, a primary study. Environ Int. 2012,

428

42, 78–83.

429

(20) Joensen, U.N.; Frederiksen, H.; Jensen, M.B.; Lauritsen, M.P.; Olesen, I.A.; Lassen, T.H.;

430

Andersson, A.M.; Jørgensen, N. Phthalate excretion pattern and testicular function: a study of 881

431

healthy Danish men. Environ Health Perspect 2012, 120(10), 1397-1403.

432

(21) Kovac, J.R.; Pastuszak, A.W.; Lamb, D.J. The use of genomics, proteomics, and

433

metabolomics in identifying biomarkers of male infertility. Fertil Steril 2013, 99(4), 998-1007.

434

(22) Guo, Y.; Wu, Q.; Kannan, K. Phthalate metabolites in urine from China, and implications for

435

human exposures. Environ Int 2011, 37(5), 893-898.

ACS Paragon Plus Environment

Environmental Science & Technology

436

(23) Ji, K.; Kho, Y.L.; Park, Y.; Choi, K. Influence of a five-day vegetarian diet on urinary levels

437

of antibiotics and phthalate metabolites: a pilot study with “Temple Stay” participants. Environ

438

Res 2010, 110(4), 375-382.

439

(24) Pan, G.; Hanaoka, T.; Yoshimura, M.; Zhang, S.; Wang, P.; Tsukino, H.; Inoue, K.; Nakazawa,

440

H.; Tsugane, S.; Takahashi, K. Decreased serum free testosterone in workers exposed to high

441

levels of di-n-butyl phthalate (DBP) and di-2-ethylhexyl phthalate (DEHP): a cross-sectional

442

study in China. Environ Health Perspect 2006, 114(11):1643-1648.

443

(25) Zhang, J.; Shen, H.Q.; Xu, W.P.; Xia, Y.K.; Barr, D.B.; Mu, X.L.; Wang, X.; Liu, L.; Huang,

444

Q.Y.; Tian, M.P. Urinary metabolomics revealed arsenic internal dose-related metabolic alterations:

445

a proof-of-concept study in a Chinese male cohort. Environ Sci Technol 2014, 48(20),

446

12265-12274.

447

(26) Xu, T.; Holzapfel, C.; Dong, X.; Bader, E.; Yu, Z.; Prehn, C.; Perstorfer, K.; Jaremek, M.;

448

Roemisch-Marg, W.; Rathmann, W.; Li, Y.; Wichmann, H.E.; Wallaschofski, H.; Ladwig, K.H.;

449

Theis, F.; Suhre, K.; Adamski, J.; Illig, T.; Peters, A.; Wang-Sattler, R. Effects of smoking and

450

smoking cessation on human serum metabolite profile: results from the KORA cohort study. BMC

451

Med 2013, 11, 60.

452

(27) Hsu, P.C.; Zhou, B.; Zhao, Y.; Ressom, H.W.; Cheema, A.K.; Pickworth, W.; Shields, P.G.

453

Feasibility of identifying the tobacco-related global metabolome in blood by UPLC-QTOF-MS. J

454

Proteome Res 2013, 12(2), 679-691.

455

(28) Jaremek, M.; Yu, Z.; Mangino, M.; Mittelstrass, K.; Prehn, C.; Singmann, P.; Xu, T.; Dahmen,

456

N.; Weinberger, K.M.; Suhre, K.; Peters, A.; Döring, A.; Hauner, H.; Adamski, J.; Illig, T.; Spector,

457

T.D.; Wang-Sattler, R. Alcohol-induced metabolomic differences in humans. Transl Psychiatry

ACS Paragon Plus Environment

Page 22 of 31

Page 23 of 31

Environmental Science & Technology

458

2013, 3:e276.

459

(29) Marcovina, S.M.; Sirtori, C.; Peracinom, A.; Gheorghiade, M.; Borum, P.; Remuzzi, G.;

460

Ardehali, H. Translating the basic knowledge of mitochondrial functions to metabolic therapy:

461

role of L-carnitine. Transl Res 2013, 161(2), 73-84.

462

(30) Banerjee, R.; Pathmasiri, W.; Snyder, R.; McRitchie, S.; Sumner, S. Metabolomics of brain

463

and reproductive organs: characterizing the impact of gestational exposure to butylbenzyl

464

phthalate on dams and resultant offspring. Metabolomics 2012, 8, 1012-1025.

465

(31) Reubsaet, F.A.; Veerkamp, J.H.; Dirven, H.A.; Brückwilder, M.L.; Hashimoto, T.; Trijbels,

466

J.M.; Monnens, L.A. The effect of di(ethylhexyl)phthalate on fatty acid oxidation and carnitine

467

palmitoyltransferase in various rat tissues. Biochim Biophys Acta 1990, 1047(3), 264-270.

468

(32) Reddy, J.K.; Moody, D.E.; Azarnoff, D.L.; Rao, M.S. Di-(2-ethylhexyl)phthalate: an

469

industrial plasticizer induces hypolipidemia and enhances hepatic catalase and carnitine

470

acetyltransferase activities in rat and mice. Life Sci 1976, 18(9), 941-945.

471

(33) Xu, Y.; Agrawal, S.; Cook, T.J.; Knipp, G.T. Maternal di-(2-ethylhexyl)-phthalate exposure

472

influences essential fatty acid homeostasis in rat placenta. Placenta 2008, 29, 962–969.

473

(34) Kristensen, D.M.; Skalkam, M.L.; Audouze, K.; Lesné, L.; Desdoits-Lethimonier, C.;

474

Frederiksen, H.; Brunak, S.; Skakkebæk, N.E.; Jégou, B.; Hansen, J.B.; Junker, S.; Leffers, H.

475

Many putative endocrine disruptors inhibit prostaglandin synthesis. Environ Health Perspect 2011,

476

119(4), 534-541.

477

(35) Meeker, J.D.; Hu, H.; Cantonwine, D.E.; Lamadrid-Figueroa, H.; Calafat, A.M.; Ettinger,

478

A.S.; Hernandez-Avila, M.; Loch-Caruso, R.; Téllez-Rojo M. Urinary phthalate metabolites in

479

relation to preterm birth in Mexico City. Environ Health Perspect. 2009,117(10), 1587–1592.

ACS Paragon Plus Environment

Environmental Science & Technology

480

(36) Kim, S.H.; Park, M.J. Serum levels of di(2-ethylhexyl) phthalate are increased in obese

481

versus nonobese children. Annual meeting of The Endocrine Society, 23-26 June 2012, Houston,

482

TX.

483

(37) Dhalluin, S.; Gate, L.; Vasseur, P.; Tapiero, H.; Nguyen-Ba, G. Dysregulation of ornithine

484

decarboxylase activity, apoptosis and Bcl-2 oncoprotein in Syrian hamster embryo cells

485

stage-exposed to di(2-ethylhexyl)phthalate and tetradecanoylphorbol acetate. Carcinogenesis 1997,

486

18(11), 2217-2223.

487

(38) Fukuwatari, T.; Ohsaki, S.; Suzuki, Y.; Fukuoka, S.; Sasaki, R.; Shibata, K. 2003 The effects

488

of phthalate esters on the tryptophan-niacin metabolism. In: Developments in tryptophan and

489

serotonin metabolism. New York: Springer, 659-664.

490

ACS Paragon Plus Environment

Page 24 of 31

Page 25 of 31

Environmental Science & Technology

491

Figure Legends

492

Figure 1. Scoring plots of the developed PCA and PLS-DA models of MBP (A) and ΣDEHP (B).

493

▲ low exposure group; ◆ medium exposure group; □ high exposure group.

494

Figure 2. ROC curves, probability views and confusion matrix of the combined biomarker patterns

495

for MBP (A) and ΣDEHP (B). ROC curves are generated by Monte-Carlo cross validation using

496

balanced subsampling. The predicted class probabilities were calculated for each sample using the

497

developed ROC models.

498 499

ACS Paragon Plus Environment

Environmental Science & Technology

Page 26 of 31

Table 1. Characteristics of the study population (n = 364) Characteristic

Mean ± SD

Median (5th, 95th)

Age (years) Height (cm) Weight (kg) BMI (kg/m2) Smoking status Never Past Current Alcohol drinking status Never Past Current

29.0 ± 4.7 173.0 ± 5.5 71.1 ± 10.9 23.7 ± 3.2

28.5 (22.6, 37.5) 173.0 (165.0, 181.8) 70.0 (55.0, 90.0) 23.7 (19.0, 29.3)

n (%)

194 (53.3%) 23 (6.3%) 147 (40.4%) 172 (47.3%) 23 (6.3%) 169 (46.4%)

Education level < College > College Missing data Annual income < 20,000 RMB > 20,000 RMB Missing data

86(23.7%) 117(32.1%) 161(44.2%)

Occupation White collar Blue collar Missing data

119(32.7%) 146(40.1%) 99(27.2%)

171(47.0%) 136(37.4%) 57(15.6%)

ACS Paragon Plus Environment

Page 27 of 31

Environmental Science & Technology

Table 2. Urinary concentrations (µg/g creatinine) of phthalates (n = 364). Metabolite

Percent>LOD

Median (5th, 95th)

Low tertile

Medium tertile

High tertile

MEP MBP MBzP MEHP MEHHP MEOHP MMP

100 100 52.7 100 100 100 3.6

28.2 (3.2, 339.5) 47.1 (6.9, 448.0) 0.8 (0.1, 7.6) 8.0 (1.1, 38.5) 22.5 (4.1, 90.7) 10.7 (2.0, 46.8) 24.4 (0.3, 112.0) a

0.8-15.3 1.5-32.1 0-0.4 0.2-5.3 0.9-14.9 0.5-7.7 0

15.3-54.6 32.1-80.5 0.4-1.2 5.3-12.8 14.9-31.8 7.7-17.2 0-8.5

54.6-10322.3 80.5-1350.4 1.2-623.9 12.8-545.8 31.8-1195.4 17.2-570.9 8.5-112

42.0 (8.5, 166.8)

1.6-28.2

28.2-65.7

65.7-2312.1

ΣDEHP a a

ΣDEHP was the mass sum of individual metabolite concentrations (MEHP, MEHHP and MEOHP).

ACS Paragon Plus Environment

Environmental Science & Technology

Page 28 of 31

Table 3. Differentially expressed urinary biomarkers associated with phthalate exposure for male population. ΣDEHP

MBP No.

Biomarkers FC (M/L)

a

b

FC (H/L)

AUC(95% CI)

c

FC (M/L)

a

FC (H/L)b

AUC(95% CI)c

1

Acetylneuraminic acid d

1.14

1.34**

0.66(0.59-0.73)**

1.05**

1.08

0.53(0.46,0.60)

2

Carnitine C8:1 d

1.09**

1.16**

0.69(0.62-0.75) **

1.07*

1.10**

0.63(0.56-0.70) **

3

Carnitine C18:0 e

1.59**

1.74**

0.76(0.70-0.82) **

1.51**

2.06**

0.79(0.73-0.84) **

4

Carnitine C16:2

d

**

**

**

**

**

**

5

Diacetylspermine d

0.66

0.51

0.81(0.76-0.87)

0.80**

0.75**

0.70(0.64-0.77) **

0.72

**

0.64

**

0.73(0.66-0.79)

**

0.74

**

0.59

**

0.72(0.66-0.79)

**

0.69

0.90

0.74(0.67-0.80)

Amino sugar metabolism

Carnitine metabolism, fatty acid beta oxidation, oxidative stress

0.84**

0.62(0.55-0.69) **

0.70

**

0.70(0.63-0.76)

**

Alanine metabolism, glycine and serine metabolism, urea cycle

0.71

**

0.65(0.58-0.71)

**

Bile acid biosynthesis, taurine and hypotaurine metabolism

Polyamine metabolism

Alanine

f

7

Taurine

d

8

Tryptophan d

0.67**

0.6**

0.77(0.71-0.83) **

0.75**

0.63**

0.73(0.66-0.79) **

9

Ornithine e

0.68**

0.54**

0.83(0.77-0.88) **

0.71**

0.61**

0.77(0.71-0.83) **

10

Cystine e

1.21**

1.3**

0.69(0.62-0.76) **

1.27**

1.37**

0.75(0.69-0.81) **

Glutathione metabolism, oxidative stress

11

Phenylglycine d

1.11

1.53*

0.58(0.50-0.65)*

1.06

1.19

0.52(0.44-0.59)

Amino acid metabolism

6

d

**

**

Phenylpyruvic acid

1.34

1.51

0.76**

0.70**

0.77(0.71-0.83) **

0.89**

0.78**

0.70(0.64-0.77) **

14

e

1.58

**

1.72

**

0.74(0.67-0.80)

**

**

1.96

**

0.79(0.73-0.85)

**

0.72

**

0.63

**

0.72(0.66-0.78)

**

0.66

0.54

**

0.73(0.67-0.80)

**

0.81**

0.71**

15

Hydroxyl-PGE2 Keto-PGE2 d

0.75**

0.62**

0.77(0.71-0.83) **

1.36

**

0.70(0.63-0.76)

**

Methylglutaconic acid d

16

0.75(0.69-0.81)

**

0.79

*

13

d

1.7

**

0.87

12

Glutamylphenylalanine

1.53

**

**

0.62

Pathway

0.69(0.62-0.76) **

a

the ratio of the means of medium exposure group and low exposure group. the ratio of the means of high exposure group and low exposure group. c area under curve (95% confidence interval, lower−upper bound) . d screened from MBP-based PLS-DA model. e screened from ΣDEHP-based PLS-DA model. f screened from both MBP and ΣDEHP-based PLS-DA models. * indicates p < 0.05, ** indicates p