Levels of Blood Organophosphorus Flame Retardants and

Jul 19, 2016 - Participants came to a mobile center for a physical examination and to provide blood and urine samples. All participants signed .... Pr...
1 downloads 10 Views 1MB Size
Subscriber access provided by CORNELL UNIVERSITY LIBRARY

Article

Levels of Blood Organophosphorus Flame Retardants and Association with Changes in Human Sphingolipid Homeostasis Fanrong Zhao, Yi Wan, Haoqi Zhao, Wenxin Hu, Di Mu, Thomas F Webster, and Jianying Hu Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.6b02474 • Publication Date (Web): 19 Jul 2016 Downloaded from http://pubs.acs.org on July 24, 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 30

Environmental Science & Technology

1 2 3

Levels of Blood Organophosphorus Flame Retardants and Association with Changes in Human Sphingolipid Homeostasis

4

Fanrong ZHAO1, Yi WAN1, Haoqi ZHAO1, Wenxin HU1, Di MU1, Thomas F. Webster2, and

5

Jianying HU1*

6 7

1

8

Peking University, Beijing 100871, China

9

2

10

Laboratory for Earth Surface Processes, College of Urban and Environmental Sciences,

Department of Environmental Health, Boston University School of Public Health, Boston,

Massachusetts, USA

11 12 13 14 15 16 17 18 19

Address for Correspondence

20

Dr. Jianying HU, College of Urban and Environmental Sciences, Peking University

21

Beijing 100871, China.

22

TEL & FAX: 86-10-62765520;

23

Email: [email protected].

1

ACS Paragon Plus Environment

Environmental Science & Technology

24

ABSTRACT: While a recent toxicological study has shown that organophosphorus flame

25

retardants (OPFRs) may disrupt sphingolipid homeostasis, epidemiologic evidence is

26

currently lacking. In this study, a total of 257 participants were recruited from Shenzhen,

27

China. Eleven OPFRs were for the first time simultaneously determined in the human blood

28

samples by ultraperformance liquid chromatography tandem mass spectrometry. Six OPFRs,

29

tributyl phosphate (TNBP), 2-ethylhexyl diphenyl phosphate (EHDPP), tris(2-chloroisopropyl)

30

phosphate (TCIPP), tris(2-butoxyethyl) phosphate (TBOEP), triethyl phosphate (TEP), and

31

TPHP were detectable in at least 90% of participants, with the median concentrations of 37.8,

32

1.22, 0.71, 0.54, 0.49, and 0.43 ng/mL, respectively. Sphingomyelin (SM) levels in the

33

highest quartile of EHDPP, TPHP, TNBP, TBOEP, TEP and TCIPP were 45.3% [95%

34

confidence interval: 38.1%, 53.0%]; 51.9% (45.5%, 58.6%); 153.6% (145.1%, 162.3%); 20.6%

35

(14.5%, 27.0%); 59.0% (52.1%, 66.2%) and 62.8% (55.2%, 70.6%) higher than those in the

36

lowest quartile, respectively, after adjusting for covariates. Sphingosine-1-phosphate (S1P)

37

levels in the highest quartile of EHDPP, TPHP and TNBP were 36% (-39%, -33%), 16%

38

(-19%, -14%) and 36% (-38%, -33%) lower than those in the lowest quartile, respectively. A

39

similar pattern emerged when exposures were modeled continuously. We for the first time

40

found the associations between OPFRs and changes in human sphingolipid homeostasis.

41

Keywords: Sphingomyelin, Sphingosine 1-phosphate; 2-Ethylhexyl diphenyl phosphate;

42

Triphenyl phosphate.

43

2

ACS Paragon Plus Environment

Page 2 of 30

Page 3 of 30

44

Environmental Science & Technology

INTRODUCTION

45

Organophosphorus flame retardants (OPFRs) are used as plasticizers, anti-foaming

46

agents and additives in floor polishes, glue, lubricants, food packaging, and hydraulic fluids.1,

47

2

48

significantly with the production phase-out and regulation of some brominated flame

49

retardants, with annual global production currently reaching approximately 200,000 tons.3

50

OPFRs have been detected in various environmental media including sediment, sewage water,

51

drinking water, dust, indoor air and biological samples.1, 4-9 Human may be exposed to OPFRs

52

through various exposure pathways including drinking water, food, indoor air, and indoor

53

dust.1,

54

neurotoxicity, and hemolytic, reproductive and cardiac effects, have been observed in animals

55

exposed

56

tris(1,3-dichloro-2-propyl) phosphate (TDCIPP) in house dust has been associated with

57

hormone levels and semen quality in men.15

The usage of OPFRs as alternative and replacement flame retardants has increased

2, 10

to

Adverse health issues, including skin irritation, carcinogenicity, dermatitis,

OPFRs,2,

11-15

and

exposure

to

triphenyl

phosphate

(TPHP)

and

58

The cardiotoxicity of OPFRs is of concern. Increased left ventricular wall thickening,

59

which is suggestive of poor cardiovascular performance, has been observed in male rats

60

exposed to Firemaster 550 (FM 550) containing TPHP as a major chemical.14 The impacts of

61

some OPFRs on cardiac looping and function have also been observed in zebrafish during

62

embryogenesis.16 In cardiovascular functions, sphingolipid homeostasis is of vital importance,

63

since sphingolipids and biosynthetic intermediates, including sphingomyelin (SM), ceramide

64

(Cer), sphingosine (Sph), and sphingosine-1-phosphate (S1P), play essential roles as both

65

structural components of cell membranes and signaling molecules that regulate cardiac 3

ACS Paragon Plus Environment

Environmental Science & Technology

66

development and barrier function of the vasculature.17 Sphigolipid levels are tightly controlled

67

by metabolic interplay of the de novo and recycling pathways (Figure S1). Especially, the

68

recycling metabolic pathway of sphingolipids including SM, Cer, Sph, and S1P, is the

69

dominant pathway to regulate the homeostasis of sphingolipid in most tissues,18, 19 and any

70

imbalance can cause stress to the cell and lead to cardiovascular disease.19 Increasing animal

71

and epidemiological evidences have shown that elevated SM plasma level is associated with

72

atherosclerosis, and sphingolipid homeostasis disorders are implicated in the pathogenesis of

73

atherosclerosis19-21 A recent study assessed sphingolipid homeostasis in mice exposed to

74

TPHP, in which sphingolipid homeostasis was disrupted as characterized by significant

75

increases in the levels of SM and decreases in its metabolite, Cer, and precursor,

76

sphinganine.22 However, to our best knowledge, there are no studies on the association of

77

OPFRs with the disruption of sphingolipid homeostasis in human.

78

While various OPFRs have been detected in milk, adipose tissue and seminal fluid and

79

their metabolites have been detected in human urine,1, 9, 23-26 only one paper made an attempt

80

to determine 9 OPFRs in human plasma and only TPHP was detected when the authors

81

developed a gas chromatography-nitrogen phosphorous detector (GC-NPD) method.9 In this

82

paper, 14 OPFRs were analyzed in 257 human blood samples using a sensitive

83

ultraperformance liquid chromatography tandem mass spectrometry (UPLC-ESI-MS/MS)

84

method. The associations between OPFR exposure levels and changes in sphingolipid levels

85

(SM, Cer, Sph, and S1P) were assessed to explore the potential impacts of OPFRs on human

86

sphingolipid homeostasis.

87 4

ACS Paragon Plus Environment

Page 4 of 30

Page 5 of 30

Environmental Science & Technology

88

MATERIALS AND METHODS

89

Study Population and Blood Collection. A total of 327 residents between 20 and 50 years of

90

age were recruited from the general population in Shenzhen, China, in November 2012. The

91

study

92

(IRB00001052-12058). Participants came to a mobile center for a physical examination and

93

to provide blood and urine samples. All participants signed informed consent when enrolled.

94

Volunteer participants were asked to fill out interview questionnaires by trained interviewers.

95

All subjects were asked to fast for 10 h overnight before collection of the fasting blood

96

sample. Blood samples were collected into heparinized brown glass bottles (CNW

97

Technologies GmbH, Germany) and stored at -80°C prior to extraction.

98

was

approved

by

the

Human

Ethics

Committee

of

Peking

University

We collected data on potential confounding variables through questionnaires. Our models

99

included a number of covariates that are important predictors of cardiovascular disease (CVD):

100

personal characteristics [age, gender and body mass index (BMI: calculated as weight in

101

kilograms divided by height in meters squared)], socioeconomic status (SES: Household

102

income was mainly considered) and life-style habits (dietary structure, alcohol intake and

103

tobacco use). Of the 327 participants, we finally got 296 blood samples with matching

104

questionnaires (Figure S2).

105

Materials. Purities of all analytical standards used in this study were ≥95%. The

106

standards of tris(2-chloroethyl) phosphate (TCEP), tripropyl phosphate (TPrP), triisopropyl

107

phosphate (TiPP), tris(2-chloroisopropyl)phosphate (TCIPP), tributyl phosphate (TNBP), and

108

tricresyl phosphate (TMPP) were purchased from Sigma-Aldrich (St. Louis, MO, USA). The

109

triethyl phosphate (TEP), tris(1,3-dichloro-2-propyl) phosphate (TDCIPP), tris(2-butoxyethyl) 5

ACS Paragon Plus Environment

Environmental Science & Technology

Page 6 of 30

110

phosphate (TBOEP), 2-ethylhexyl diphenyl phosphate (EHDPP), and tris(2-ethylhexyl)

111

phosphate (TEHP) were purchased from TCI Corp. (Tokyo, Japan). Tris(2,3-dibromopropyl)

112

phosphate (TDBPP) was supplied by AccuStandard Inc. (New Haven, CT, USA). Bisphenol A

113

bis(diphenyl phosphate) (BPA-BDPP) was purchased from Toronto Research Chemicals Inc.

114

(Toronto, Ontario, Canada). The SM, Cer, Sph and S1P were purchased from Avanti Polar

115

Lipids (Alabastar, AL, USA). The stable isotope labeled standards, including triethyl-d15

116

phosphate (TEP-d15), tri-n-propyl-d21 phosphate (TPrP-d21) and triphenyl-d15 phosphate

117

(TPHP-d15), were supplied by C/D/N Isotopes Inc. (Pointe-Claire, Quebec, Canada).

118

Tri-n-butyl-d27 phosphate (TNBP-d27) was purchased from Cambridge Isotope Laboratories

119

Inc.

120

tris(1,3-dichloro-2-propyl)-d15 phosphate (TDCIPP-d15) were obtained from Toronto Research

121

Chemicals Inc. (Toronto, Ontario, Canada). Tri-p-cresyl-d21 phosphate (TMPP-d21) and

122

tris(2-ethylhexyl)-d51 phosphate (TEHP-d51) were supplied by Hayashi Pure Chemical Ind.,

123

Ltd. (Osaka, Japan). Internal standards SM-d31, Cer-d31, Sph-d7 and S1P-d7 were purchased

124

from Avanti Polar Lipids (Alabastar, AL, USA). Solvents, including n-hexane, ethyl acetate

125

and methanol (MeOH), were of pesticide residue grade and obtained from Fisher Chemicals

126

(New Jersey, USA). Acetone and dichloromethane (pesticide residue grade) were purchased

127

from Mallinckrodt Baker Inc. (Phillipsburg, NJ, USA). Formic acid (HPLC grade) was from

128

Dikma Technologies Inc. (California, USA). Sep-Pak® Silica (3 cm3, 200 mg) and Sep-Pak®

129

C18 (3 cm3, 200 mg) solid phase extraction (SPE) cartridges were purchased from Waters

130

(Milford, MA, USA). Lithium-heparin blood collection tubes were from Corning

131

Incorporated (Tewksbury, MA, USA). Ultrapure water was prepared using a Milli-Q

(Tewksbury,

MA,

USA).

Tris(2-chloroethyl)-d12

phosphate

6

ACS Paragon Plus Environment

(TCEP-d12)

and

Page 7 of 30

Environmental Science & Technology

132

Synthesis water purification system (Millipore, Bedford, MA, USA). A full list of the 14

133

target compounds, along with their full chemical names and structures, can be found in Figure

134

S3.

135

OPFRs analysis. Blood samples were thawed on ice, and prepared immediately for

136

analysis. Each whole blood sample (0.4 mL) was transferred into an 8-mL glass centrifuge

137

tube, and 20 µL of internal standard solution (10 ng/mL for each surrogate) was added. Ethyl

138

acetate (2 ml) was added to the sample, which was shaken for 20 min on an orbital shaker and

139

centrifuged at 4000 rpm for 10 min, with the ethyl acetate layer then transferred to a clean

140

glass bottle. Extraction from the residue was repeated twice and the organic layers were

141

combined, concentrated to near dryness under a gentle stream of nitrogen and redissolved in

142

500 µL of n-hexane. The concentrated extract was loaded on a silica cartridge preconditioned

143

with 6 mL of hexane/acetone (1:1, v/v), 3 mL of hexane/dichloromethane (7:3, v/v) and 3 mL

144

of hexane. After the cartridge was rinsed with 3 mL of hexane and 3 mL of

145

hexane/dichloromethane (7:3, v/v), the target analytes was eluted by 3 mL of

146

n-hexane/acetone (1:1, v/v). The extracts were then dried and redissolved in 2 mL of methanol

147

to pass through a C18-SPE cartridge preconditioned by 6 mL of methanol. The filtrate was

148

collected and evaporated to dryness under a gentle stream of nitrogen and reconstituted with

149

100 µL of methanol for UPLC-MS/MS analysis.

150

Analysis of OPFRs was performed using a Waters Acquity UPLCTM system (Waters,

151

Milford, MA, USA). All OPFRs were separated using a Waters Acquity UPLC BEH C8

152

column (2.1 mm × 100 mm × 1.8 µm) preceded by a Waters Acquity UPLC BEH C18 guard

153

column (2.1 mm × 50 mm × 1.7 µm). The column was maintained at 40°C, with a flow rate of 7

ACS Paragon Plus Environment

Environmental Science & Technology

154

0.2 mL/min and injection volume of 5 µL. Methanol (A) and ultrapure water containing 0.1%

155

(v/v) formic acid (B) were used as the mobile phases. Detailed information on the UPLC

156

gradient conditions is shown in the Supporting Information.

157

Mass spectrometry was performed using a Waters Micromass Quattro Premier XE triple

158

quadruple instrument detector equipped with an electrospray ionization source (Micromass,

159

Manchester, UK) in positive ion mode. The optimized parameters were: source temperature,

160

110°C; desolvation temperature, 350°C; capillary voltage, 3.50 kV; desolvation gas flow, 800

161

L/h; cone gas flow, 50 L/h; and multiplier, 650 V. Finally, MS/MS data acquisition was

162

performed in multiple-reaction monitoring (MRM) mode, and time-segmented scanning in

163

four functions was used based on the chromatographic separation of target compounds to

164

maximize detection sensitivity. The MS/MS parameters for the analytes, including their

165

precursors and product ions, cone voltage and collision energy, are summarized in Table S1.

166

All analytical procedures were checked for accuracy, precision, reproducibility, linearity,

167

blank contamination, matrix spikes, method limits of detection (LODs) and limits of

168

quantification (LOQs). Matrix-spiked recoveries of individual OPFRs through the analytical

169

procedure were estimated by spiking the target analytes at 5, 50 and 100 ng/mL for TNBP, 1,

170

10 and 20 ng/mL for EHDPP, TCEP and TCIPP, 0.25, 1 and 10 ng/mL for TEP, and 0.1, 1 and

171

10 ng/mL for the others into sample matrices (n=6). The low concentrations of spiked

172

analytes were similar to those in samples. The method recoveries of all target analytes were

173

calculated by subtracting background concentrations in nonspiked samples from spiked

174

samples, and the recoveries ranged from 74% to 98% (with a RSD of 3−15%), 73% to 101%

175

(RSD: 2−7%) and 73% to 99% (RSD: 2−11%) for the low, medium and high concentrations, 8

ACS Paragon Plus Environment

Page 8 of 30

Page 9 of 30

Environmental Science & Technology

176

respectively. Relative standard deviation (RSD) was used to evaluate precision. The inter-day

177

precisions were calculated based on the means of six spiked samples at three different levels

178

during five days and the inter-day precisions of all substances were within 15%. To prevent

179

possible specimen contamination, only pretreated glassware (500°C, 6 h) were used

180

throughout the study, and aluminized paper was used in all plastic seals to minimize possible

181

contamination of the samples during sampling, storage, transport and extraction. The SPE

182

cartridges, silica and C18 were pre-rinsed with n-hexane/acetone (1:1, v/v) and MeOH prior

183

to use, respectively, to minimize the contamination of the SPE procedure. For each batch of

184

20 samples analyzed, 2 procedural blanks were processed. Procedural blanks were prepared

185

by substitution of 0.4 mL of Milli-Q water for blood, followed by passage through the entire

186

analytical procedure. TEP, TCEP, TCIPP, TPHP, TDCIPP, TNBP, TMPP, TBOEP, EHDPP,

187

BPA-BDPP, and TEHP were detected in the procedural blanks at the concentrations of

188

0.076±0.010, 0.106±0.031, 0.53±0.062, 0.084±0.007, 0.052±0.007, 0.64±0.17, 0.082±0.006,

189

0.012±0.005, 0.024±0.004, 0.018±0.002 and 0.016±0.002 ng/mL, respectively.

190

(n = 6) were performed by transferring 8 mL of Milli-Q water (with the same vacuum tubes

191

used for collecting blood samples) into sampling containers, storing and processing them as

192

samples. The concentrations of target compounds in field blanks were almost equal to those in

193

procedural blanks. Calibration curves of standards of target analytes were calculated with a

194

concentration series of 0.01, 0.05, 0.1, 0.5, 1, 5, 10, 50 and 100 ng/mL, except for TNBP and

195

TCEP (0.1, 0.5, 1, 5, 10, 50, 100, 500 and 1000 ng/mL). All calibration curves provided

196

adequate linearity (r2>0.995), and the signal-to-noise ratios for the lowest concentration on

197

calibration curves of OPFRs were 13−98. Concentrations of TCEP and TNBP in two samples 9

ACS Paragon Plus Environment

Field blanks

Environmental Science & Technology

198

were reanalyzed after diluting the extracts by 4-10 times due to the extremely high

199

concentrations.

200

Identification of the target analytes was accomplished by comparing the retention time

201

(within 2%) and ratio (within 20%) of the two selected precursor ion-produced ion transitions

202

with those of the standards. Quantification was accomplished using the MRM transitions. To

203

automatically correct for the losses of analytes during extraction or sample preparation and to

204

compensate for variations in instrumental response from injection to injection, quantification

205

of the analytes was achieved using an internal standard method with calibration against

206

standard solutions. TEP-d15 was used as the surrogate standard for TEP; TCEP-d12 for TCEP

207

and TCIPP; TDCIPP-d15 for TDCIPP and TDBPP; TPrP-d21 for TiPP and TPrP; TPHP-d15 for

208

TPHP and EHDPP; TNBP-d27 for TNBP; TMPP-d21 for TMPP; and TEHP-d51 for TBOEP,

209

BPA-BDPP and TEHP. For chemicals with detectable blank contamination, the LODs and

210

LOQs were calculated as three and ten times the standard deviations of procedural blanks,

211

respectively, and final concentrations were calculated by the initial concentrations subtracting

212

the blank values. For TiPP, TPrP and TDBPP, which were not detected in the blanks, LODs

213

and LOQs were calculated based on signal-to-noise ratios of 3 and 10 in matrix-spiked

214

samples, respectively. The LODs and LOQs were 0.004–0.52 and 0.02–1.73 ng/mL,

215

respectively (Table S2), which are significantly lower than the LODs of OPFRs in human

216

plasma (0.2–1.8 ng/mL) in a previous paper,27 showing high sensitivity using developed

217

method.

218

Sphingolipids Analysis. Sphingolipids were analyzed by UPLC-MS/MS using the

219

method previously reported with minor modifications.16,17 Briefly, after loading blood 10

ACS Paragon Plus Environment

Page 10 of 30

Page 11 of 30

Environmental Science & Technology

220

samples (100 µL for each test) into a 1.5 ml centrifuge tube, MeOH (880 µL) and 20 µl

221

internal standard solution (2000 ng/mL for each) were added to the same tube, which was then

222

shaken vigorously for 1 min. The concentrated extract (1 mL) was stored at -20°C for 24 h,

223

and most of the precipitated or suspended lipids were easily removed by filtration. After

224

centrifugation at 12,000 rpm for 15 min, the supernatant was collected and analysis of

225

shingolipid was performed by a Waters Acquity UPLCTM system. A Waters Acquity UPLC

226

BEH phenyl column (2.1 mm × 100 mm × 1.8 µm) was used for chromatographic separation.

227

Methanol (A) and ultrapure water containing 0.5% formic acid (B) were used as the mobile

228

phases. The gradient started at 10% A and then increased linearly to 60% in 6 min, 100% at 6

229

min and kept for 2 min, followed by a decrease to the initial conditions of 10% A and held for

230

2 min to allow for equilibration. The column was maintained at 40°C, with a flow rate of 0.3

231

mL/min and injection volume of 5 µL. Quantification of sphingolipids was accomplished

232

using the MRM transitions. The optimized MS/MS parameters for the analyses, including

233

precursor and product ions, cone voltage and collision energy, are shown in Table S3.

234

Mean recoveries of SM-d31, Cer-d3, Sph-d7 and S1P-d7 were 88% ± 9%, 91% ± 8%, 106%

235

± 4% and 103% ± 6%, respectively, by spiking internal standards into blood samples at low,

236

median and high levels (100, 500, and 1000 ng/mL for SM-d31; 1000, 5000 and 10000 ng/mL

237

for Cer-d31; 200, 1000 and 2000 ng/mL for Sph-d7 and 2000, 10000 and 20000 ng/mL for

238

S1P-d7) (n = 6 for each level). The LOQs for internal standards of SM, Cer, Sph, and S1P in

239

matrix-spiked samples were estimated to be 7.3, 40, 12 and 25 ng/mL, respectively based on

240

the peak-to-peak noise of the baseline and on a minimum signal-to-noise value of 10,

241

respectively. Calibration curves of standards of target analytes were calculated with a 11

ACS Paragon Plus Environment

Environmental Science & Technology

242

concentration series of 50, 100, 500, 1000, 2000, 5000, 10000, and 50000 ng/mL. All

243

calibration curves showed strong linearity (correlation coefficients >0.99) with good precision

244

(RSD ≤ 5%).

245

Statistical Analysis. Basic descriptive statistics were derived for population

246

characteristics, blood sphingolipid levels and blood OPFR levels. Whole blood samples with

247

non-detectable OPFR concentrations were assigned a value as the LOQ divided by the square

248

root of 2. Spearman correlation coefficients were calculated to assess bivariate relationships

249

between OPFR concentrations and sphingolipid levels. Sphingolipid concentrations were

250

transformed to the natural log for statistical analyses as they were log-normally distributed.

251

For OPFRs detected at over 90%, we used multivariate linear regression to assess associations

252

between OPFR concentrations and sphingolipid levels, while adjusting for all the relevant

253

covariates (age, gender, BMI, SES, dietary structure, alcohol intake and tobacco use). Quartile

254

variables were constructed for individual OPFRs in multivariate linear regression models. We

255

presented effect estimates for each quartile compared with the first quartile and their

256

corresponding 95% confidence intervals (CIs). Tests for trends in quartile analyses were

257

performed by treating the OPFR category as a linear predictor in the models. To further

258

confirm the associations between sphingolipid and OPFR concentrations, OPFRs were also

259

modeled as continuous predictors.

260

All analyses excluded those with a history of cardiovascular events (coronary event or

261

stroke) or heart failure. In addition, subjects who reported current use of medications or who

262

were missing this variable were also excluded. The details of people in the exclusion group

263

were shown in Figure S2. 12

ACS Paragon Plus Environment

Page 12 of 30

Page 13 of 30

Environmental Science & Technology

264

Data analysis was performed using SPSS, Version 22.0 (IBM, Corp, 2013). Models were

265

adjusted for relevant covariates. For all tests, an alpha level of 0.05 was chosen; p-values of
90% of participants; TEHP, TMPP, TCEP, TDCIPP and

283

BPA-BDPP were detected in 71.6%, 66.1%, 63.0%, 47.1% and 41.2% of participants,

284

respectively (Table 3). Among the 11 detected OPFRs, TNBP was the most abundant with a

285

median value of 37.8 ng/mL (