The interaction effects of smoking and polycyclic aromatic hydrocarbons exposure on the prevalence of metabolic syndrome in coke oven workers

The interaction effects of smoking and polycyclic aromatic hydrocarbons exposure on the prevalence of metabolic syndrome in coke oven workers

Journal Pre-proof The interaction effects of smoking and polycyclic aromatic hydrocarbons exposure on the prevalence of metabolic syndrome in coke ove...

8MB Sizes 0 Downloads 71 Views

Journal Pre-proof The interaction effects of smoking and polycyclic aromatic hydrocarbons exposure on the prevalence of metabolic syndrome in coke oven workers Bin Zhang, Baolong Pan, Xinyu Zhao, Ye Fu, Xuejing Li, Aimin Yang, Qiang Li, Jun Dong, Jisheng Nie, Jin Yang PII:

S0045-6535(20)30072-2

DOI:

https://doi.org/10.1016/j.chemosphere.2020.125880

Reference:

CHEM 125880

To appear in:

ECSN

Received Date: 18 November 2019 Revised Date:

8 January 2020

Accepted Date: 8 January 2020

Please cite this article as: Zhang, B., Pan, B., Zhao, X., Fu, Y., Li, X., Yang, A., Li, Q., Dong, J., Nie, J., Yang, J., The interaction effects of smoking and polycyclic aromatic hydrocarbons exposure on the prevalence of metabolic syndrome in coke oven workers, Chemosphere (2020), doi: https:// doi.org/10.1016/j.chemosphere.2020.125880. This is a PDF file of an article that has undergone enhancements after acceptance, such as the addition of a cover page and metadata, and formatting for readability, but it is not yet the definitive version of record. This version will undergo additional copyediting, typesetting and review before it is published in its final form, but we are providing this version to give early visibility of the article. Please note that, during the production process, errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain. © 2020 Published by Elsevier Ltd.

1

The Interaction Effects of Smoking and Polycyclic Aromatic

2

Hydrocarbons Exposure on the Prevalence of Metabolic Syndrome in

3

Coke Oven Workers

4

Bin Zhang1, Baolong Pan1, 2, Xinyu Zhao1, Ye Fu1, Xuejing Li1, Aimin Yang3, Qiang

5

Li4, Jun Dong4, Jisheng Nie1, Jin Yang1*

6 7

1

8

University.

9

2

General Hospital of Taiyuan Iron & Steel (Group) Co., Ltd.

10

3

Hong Kong Institutes of Diabetes and Obesity, the Chinese University of Hong

11

Kong.

12

4

13

Ltd.

14

Bin Zhang and Baolong Pan contributed equally to this paper.

15

* Corresponding author. Department of Occupational Health, School of Public Health,

16

Shanxi Medical University, Taiyuan, China. Xinjiannan Road 56, Taiyuan 030001,

17

Shanxi, China.

18

Tel/fax: +86 (351) 4135 240

19

E-mail: [email protected]

Department of Occupational Health, School of Public Health, Shanxi Medical

Center of Occupational Disease Prevention, Xishan Coal Electricity (Group) Co.,

1

20 21

Abstract Introduction: Metabolic syndrome (MetS) is a cluster of interrelated risk factors,

22

which increase the risk of cardiovascular disease (CVD) and cancer. The prevalence

23

of MetS might be affected by environmental pollution and individual's poor lifestyles.

24

Methods: In this cross-sectional study, we aimed to evaluate the interactions effect of

25

PAHs exposure and poor lifestyles on MetS among coke oven workers. We measured

26

the concentrations of 11 urinary PAH metabolites among 682 coke oven workers by

27

HPLC-MS. China adult blood lipid abnormality prevention guide (2016) was

28

employed for diagnosing MetS. An interaction effect was tested in the modified

29

Poisson regression models. Results: The results showed that the urinary level of

30

1-NAP (P for trend = 0.043) and 2-FLU (P for trend = 0.037) had significant

31

dose-response relationships with increased PR of MetS. For 1-NAP, statistically

32

significant positive association was with low HDL among individuals (P for trend =

33

0.003). Besides, smoking was associated with a significantly increased risk of

34

prevalence ratio of MetS (PR = 1.07; 95% CI 1.00-1.13), high triglycerides (PR=1.13;

35

95% CI 1.05-1.20) and low HDL (PR=1.07; 95% CI 1.01-1.14). Smokers who with

36

high levels of 1-NAP and 2-FLU had higher prevalence ratio of MetS compared with

37

non-smokers who with low levels of 1-NAP [PR (95% CI): 1.17 (1.06-1.29); P for

38

trend = 0.002] and 2-FLU [PR (95% CI): 1.17 (1.06-1.29); P for trend = 0.004].

39

Conclusions: Our findings suggested PAHs exposure increased the prevalence ratio of

40

MetS and this effect can be modified by smoking status.

41

Keywords: Polycyclic aromatic hydrocarbons; Smoking; Metabolic Syndrome; 2

42

Interaction effects

43

Abbreviations

44

MetS: Metabolic syndrome

45

BMI: Body Mass Index

46

CVD: cardiovascular disease

47

PAHs: Polycyclic aromatic hydrocarbons

48

2-NAP: 2-hydroxynaphthalene

49

1-NAP: 1-hydroxynaphthalene

50

3-FLU: 3-hydroxyfluorene

51

2-FLU: 2-hydroxyfluorene

52

2-PHE: 2-hydroxyphenanthrene

53

9-PHE: 9-hydroxyphenanthrene

54

1-PHE: 1-hydroxyphenanthrene

55

1-PYR: 1-hydroxypyrene

56

3-CHR: 3-hydroxychrysene

57

6-CHR: 6-hydroxychrysene

58

9-BAP: 9-hydroxybenzpyrene

3

59 60

1. Introduction Metabolic syndrome (MetS) is a cluster of interrelated risk factors (high blood

61

pressure, dyslipidaemia, high glucose, and abdominal obesity), which increase the risk

62

of cardiovascular disease (CVD) and cancer (Alberti et al., 2009). While excessive

63

macronutrient intakes and physical inactivity have been identified as major

64

contributors, the high burden of MetS remains not fully explained (Grundy, 2016).

65

Polycyclic aromatic hydrocarbons (PAHs) are the principal pollutant of coke

66

oven emissions generated by incomplete combustion in coking production. PAHs

67

have strong lipophilic properties making them quickly absorbed by the fatty tissues

68

such as the kidney and liver within body and capable of being stored in fat cells and

69

tissues containing fat, and easily accumulated through repeated and long-term

70

exposures or bio-accumulate through the food chain (Simkhovich et al., 2008;

71

Scinicariello and Buser, 2014). Recent studies also studying the adverse health effects

72

of PAHs have focused on increasing risks of chronic diseases such as cardiovascular

73

diseases (CVD) and a variety of cancers (Burstyn et al., 2005; Coogan et al., 2012).

74

Previous studies also identified significant associations between PAH exposure and

75

MetS and increased risks of cardiovascular diseases (Chen et al., 2008; Kuang et al.,

76

2013; Liu et al., 2018). It is thus possible that PAHs exposure induce MetS by

77

affecting the metabolic disorder of the body.

78

Male workers account for large proportion in coke oven workers. Meanwhile,

79

coke oven workers were more likely to have poor lifestyles, like smoking, drinking

80

and fewer intakes of vegetables. And smoking was well-known risk factors for many 4

81

metabolic diseases, including cancer, chronic inflammation and endothelial

82

dysfunction (Lu et al., 2014; Lu et al., 2017). Cigarette smoking is also associated

83

with abdominal obesity and dyslipidaemia (Slagter et al., 2013; Keith et al., 2016).

84

Smoking and PAHs were both the common risk factors for coke oven workers.

85

However, there are few studies on the interaction between PAHs and poor lifestyles

86

on the prevalence of MetS yet. The present study aimed to evaluate the influence of

87

the interaction between PAHs and poor lifestyles on risk of MetS, using a modified

88

Poisson regression models from a population based cross sectional study in Chinese

89

coke oven workers.

90

2. Methods

91

2.1 Study population

92

The basic demographic data was collected from a coke oven plant in China using

93

a cross-sectional survey in 2017. 859 workers participated in the study. We restricted

94

our analyses to who had worked for more than 1 year. We excluded individuals who

95

were missing with sufficient blood samples and sufficient urine samples (n = 76), or

96

demographic characteristics (n = 83). We excluded individuals who were missing

97

waist circumference measurements (n = 2), fasting blood-glucose (n = 6), blood

98

pressure readings (n = 11), triglyceride measurements (n = 6), high-density

99

lipoprotein (n = 5). Thus, our final analytic sample was 682 participants.

100

Trained interviewers collected the information regarding sex, age, years of

101

working, cooking fumes, smoking and drinking status, body mass index and eating

102

habits (salt, vegetables and fruits) and occupational exposure history by a pre-tested 5

103

questionnaire. Smokers were defined as those who smoked at least 1 cigarette every

104

day and continuously more than six months, and drinkers were drank at least once a

105

week on average and continuously more than six months. Venous blood (5 ml) and

106

morning urine (20 ml) were provided by each participant. Every participant signed the

107

informed consent and the study was approved by the Medical Ethics Committee of the

108

Shanxi Medical University.

109

2.2 Urine PAH metabolites

110

According to the recommendation of the American Conference of Governmental

111

Industrial Hygienists (ACGIH), the morning urine samples of end-of-work-week were

112

obtained. All urine samples were freezed at -80°C until further processing. We used

113

high performance liquid chromatography mass spectrometry (HPLC-MS) to detect

114

urine biomarkers of PAHs exposure, including 2-hydroxynaphthalene (2-NAP),

115

1-hydroxynaphthalene (1-NAP), 3-hydroxyfluorene (3-FLU), 2-hydroxyfluorene

116

(2-FLU), 2-hydroxyphenanthrene (2-PHE), 9-hydroxyphenanthrene (9-PHE),

117

1-hydroxyphenanthrene (1-PHE), 1-hydroxypyrene (1-PYR), 3-hydroxychrysene

118

(3-CHR), 6-hydroxychrysene (6-CHR) and 9-hydroxybenzpyrene (9-BAP) levels.

119

The linearity (expressing as R2), limit of detection (LOD), reproducibility (expressing

120

as coefficient of variation (CV)) and mean recovery rate were 0.9989-1.000, 0.001 -

121

0.014 ng/ml, 0.55%-3.48% and 71.44 %- 121.20 %, respectively (Table S1). The

122

concentrations less than LOD were expressed with half a LOD value. Valid urine

123

concentrations of PAH metabolites were adjusted using urine gravity.

124

2.3 Metabolic syndrome and component conditions 6

125

China adult blood lipid abnormality prevention guide (2016) was employed for

126

diagnosing MetS. MetS was defined according to the harmonized definition as the

127

presence of at least 3 of following component conditions: 1) abdominal obesity (waist

128

circumference of ≥ 85cm for women or ≥ 90cm for men); 2) high fasting blood

129

glucose levels ( ≥ 6.10mmol/L or current use of medication to treat hyperglycaemia);

130

3) high blood pressure (systolic blood pressure ≥ 130 mmHg, diastolic blood pressure

131

≥ 85 mmHg, or current use of medication to treat high blood pressure); 4) high

132

triglyceride levels ( ≥ 1.7mmol/L, or current use of medication to treat elevated

133

triglycerides); or 5) low HDL cholesterol levels ( < 1.0mmol/L or current use of

134

medication to treat reduced HDL) .

135

Examinations and laboratory measures were conducted by trained nurse and

136

physicians in Center of Occupational Disease Prevention of Xishan Coal Electricity

137

(Group) Co., Ltd. Waist circumference measurements (cm) were taken at the level of

138

the umbilicus in mid-respiration while the participant was standing. After sitting for

139

10 min, systolic and diastolic blood pressures (mmHg) were measured twice within a

140

5-minute interval, and the two measurements were averaged. Venous blood samples

141

were obtained from each participant early in the morning, with at least 8 h of fasting.

142

Laboratory measures included glucose, total serum cholesterol, low-density

143

lipoprotein, high-density lipoprotein and triglycerides.

144

2.4 Statistical analysis

145 146

All statistical analyses were performed using SAS software version 9.4 (SAS Institute Inc., Cary, NC). We used the median (Med) and quartiles, frequency and 7

147

proportion to describe the basic characteristics of participants by MetS among 682

148

occupational workers. Crude data were compared by applying Chi-square test for

149

categorical variables, and the Kruskal-Wallis H test for numerical variable. Spearman

150

was used to assess the correlation of PAH metabolites.

151

We tested for linear trends across tertiles of PAH metabolites by including the

152

median of each tertile as a continuous variable in modified Poisson regression

153

models. Multivariable regression analyses included adjustment for sex, age (< 32,

154

32-46 and ≥ 46 years) and smoking status (yes or no), drinking status (yes or no)

155

(Model 1), and further adjustment for cooking fumes (few or much), salt (light,

156

medium and heavy), vegetables (< 2.5 and ≥ 2.5 kg/week), fruits (< 750 and ≥ 750

157

g/week) and other PAH metabolites. BMI (kg/m2) was entered as an additional

158

covariate, except for in models of abdominal obesity so as to avoid over adjustment

159

(Model 2). The PAHs metabolites were detectable in the majority of participants in

160

addition to 3-CHR, 6-CHR and 9-BAP. The percent of below LOD of 3-CHR and

161

6-CHR were more than 50% (58.4% and 74.0%). So, 3-CHR and 6-CHR were not

162

put in the regression analyses. And, 9-BAP was adjusted as a covariant in categorical

163

variable (< LOD, 0.010- 0.024, ≥0.024), because the detectable rate of 9-BAP close

164

to 50% (49.3%). In sensitivity analysis, we further performed additional 3-CHR and

165

6-CHR in the analysis as covariates to explore the associations between PAH

166

metabolites and the prevalence ratios of MetS. The dose-response relationship

167

between PAH metabolites with prevalence ratio of MetS were plotted using natural

168

cubic spline with three degree of freedom for PAH metabolites term in the fully 8

169

adjusted modified Poisson regression models. Linearity was tested by comparing the

170

model fit of the linear and the spline model using a log likelihood ratio Chi-square

171

test. Modified Poisson regression models were performed to estimate the prevalence

172

ratios of MetS for smoking status. To determine the interaction between smoking and

173

PAH metabolites levels, an interaction term was tested in the modified Poisson

174

regression models. Differences considered statistically significant for P values were

175

<0.05.

176

3. Results

177

3.1 Main characteristics and PAHs exposure of study subjects

178

Characteristics of the study population, stratified by MetS diagnosis, are depicted

179

in Table 1. Among the 682 people included in this study, 14.67% (100/682) had MetS,

180

and the mean age was 45.5 years old. Of these subjects, 61% and 38% were smokers

181

and drinkers, respectively. As expected, the median of age in the MetS group were

182

significantly older than the non-MetS group (P < 0.001), the group with MetS had

183

higher BMI (P < 0.001). There was more male (89/100) with MetS, while people who

184

have MetS were more likely to be exposed to cooking fumes (P < 0.05). And the

185

levels of 2-FLU were higher among the MetS group (P < 0.001).

186

Urine PAH metabolites were correlated with each other. The correlation ranged

187

from 0.74 to -0.06. The correlation between 9-PHE and 1-PHE was stronger (r = 0.74),

188

and there was a very slight negative correlation between 1-NAP and 9-BAP (r = -0.06)

189

(Figure S1).

190

3.2 Associations between PAHs exposure and MetS risk 9

191

We used the modified Poisson regression analyses to assess the association

192

between PAH metabolites concentration and the prevalence ratio of MetS (Table 2). In

193

model 2, after adjustment for demographic, lifestyles covariates (i.e. sex, age,

194

smoking, drinking, cooking fumes, eating habits and BMI) and other PAH metabolites,

195

the PRs (95% CI) of MetS for increasing tertiles of 2-FLU were 1.00 (reference), 1.07

196

(1.01-1.12) and 1.10 (1.02-1.18), respectively (P for trend = 0.037). 1-NAP was

197

positively related to the PR of MetS in a dose-dependent manner (3rd vs. 1st tertile:

198

PR=1.08, 95% CI: 1.00-1.16; P for trend = 0.043). No association was observed

199

between other PAHs metabolites and the prevalence of MetS. We also tested for linear

200

trends across tertiles of PAH metabolites by including the median of each tertile as a

201

continuous variable in logistic regression models (Table S2). The result was consistent

202

with modified Poisson regression models.

203

The associations between PAH metabolites and MetS components among the

204

study population are provided in Figure 1. For 1-NAP, the only statistically significant

205

positive association was with low HDL among individuals with the highest tertile

206

(PR=1.11, 95% CI: 1.04-1.19; P for trend = 0.003). Urinary 2-PHE in the 3rd tertile

207

were associated with an 8% (95% CI: 1.00-1.18; P for trend = 0.030) greater

208

prevalence ratio of abdominal obesity. For 1-PYR, the only statistically significant

209

positive association was with high glucose among individuals with the highest tertile

210

(PR=1.08, 95% CI: 1.01-1.16; P for trend = 0.008). 3-FLU concentrations were

211

negatively related to high triglycerides and abdominal obesity (P for trend < 0.05).

212

Significant linear associations of 9-PHE (3rd vs. 1st tertile: PR=0.92, 95% CI: 10

213

0.87-0.98; P for trend = 0.002) were observed for low HDL. And, significant linear

214

associations of 9-PHE (3rd vs. 1st tertile: PR=0.91, 95% CI: 0.85-0.98; P for trend =

215

0.021) was observed for high glucose.

216

Comparison of the linear and spline models suggested the exposure-response

217

relationship was essentially linear for the prevalence ratio of MetS (2-FLU: P = 0.877,

218

1-NAP: P = 0.947) and the prevalence ratio of low HDL (P = 0.988) (Figure S2).

219

Sensitivity analyses controlling for 3-CHR and 6-CHR did not appreciably alter

220

the results. 1-NAP (P for trend =0.044) and 2-FLU (P for trend =0.038) were still

221

positively related to MetS in a dose-dependent manner. 1-NAP concentrations were

222

associated with lower HDL (P for trend =0.001). The associations between PAH

223

metabolites and MetS components with additional adjustment for 3-CHR and 6-CHR

224

are provided in Supplemental Table 3.

225

3.3 Associations between lifestyles and MetS risk

226

Adjusted prevalence ratios for MetS component conditions by lifestyles were

227

assessed using the modified Poisson regression model in Figure 2. After adjusting for

228

multiple confounding factors, smoking was associated with a significantly increased

229

risk of prevalence ratio of MetS (PR=1.07; 95% CI 1.00-1.13), high triglycerides

230

(PR=1.13; 95% CI 1.05-1.20) and low HDL (PR=1.07; 95% CI 1.01-1.14),

231

respectively. Drinking (PR=1.22; 95% CI 1.02-1.46) and fruits intake (PR=1.20; 95%

232

CI 1.02-1.41) were associated with high blood pressure. However, no significant

233

association was observed between other lifestyles (salt and vegetables intake) and the

234

prevalence of MetS or components. 11

235 236

3.4 Effects of PAH metabolites and smoking on MetS components We tested the contribution rates of smoking on urine PAH metabolites

237

(Supplemental Table 4), and found smoking accounted for 0.5% of the urinary 1-NAP

238

variance, 0.1% of the urinary 2-FLU variance. There were the low contribution rates

239

of smoking on 1-NAP and 2-FLU, so we can choose them as biomarkers of PAH

240

exposure to explore the co-exposure effect of smoking and occupational PAH on the

241

risk of MetS and components.

242

The PRs for association of smoking and urine 1-NAP, 2-FLU co-exposure with

243

high prevalence of MetS were presented in Figure 3(A and B). After adjusting

244

covariates, we observed smokers who with high 1-NAP and 2-FLU levels had

245

significantly higher prevalence of MetS compared with non-smokers who with low

246

1-NAP [PR (95% CI): 1.17 (1.06-1.29); P for trend = 0.002] and 2-FLU levels [PR

247

(95% CI): 1.17 (1.06-1.29); P for trend = 0.004]. Smokers, no matter exposed to low

248

or high levels of urine 1-NAP and 2-FLU, had an increasing risk of high prevalence

249

ratio of MetS compared with non-smokers.

250

The PRs for association of smoking and urine 1-NAP co-exposure with low HDL

251

were presented in Figure 3C. After adjusting for multiple confounding factors,

252

smokers who with high 1-NAP levels had significantly low HDL compared with

253

non-smokers who with low 1-NAP [PR (95% CI): 1.22 (1.10-1.34); P for trend =

254

0.0001].

255

4. Discussion

256

In this study, we used a cross-sectional study of Chinese coke oven workers to 12

257

investigate the effects of PAHs exposure and smoking on the prevalence of MetS. The

258

present study showed that status of smoking can modify the association between

259

PAHs exposure and the prevalence of MetS. Smokers, no matter exposed to low or

260

high levels of urine 1-NAP and 2-FLU, had an increasing risk of high prevalence ratio

261

of MetS compared with non-smokers. We observed that smokers who with high

262

1-NAP or 2-FLU levels had significantly higher prevalence ratio of MetS compared

263

with non-smokers who with low 1-NAP or 2-FLU. And, smokers who with high

264

1-NAP levels had significantly low HDL compared with non-smokers who with low

265

1-NAP. These results provide evidence and potential explanations for the roles of

266

environmental pollution and poor lifestyles interactions on increased risks of MetS.

267

The etiology of MetS is very complex, and there may be obesity and adipose

268

tissue disease, insulin resistance, environmental and genetic factors, etc. Other factors

269

(such as age, pre-inflammatory state and hormone changes) also play a role. Recently,

270

some scholars believe that industrial toxicants may be one of the pathogenic factors of

271

MetS. The effects of PAHs exposure on lipids metabolism of mice have been

272

observed in the previously experimental studies. Lipids profiles were significantly

273

changed in mice serum after benzopyrene exposure (Li et al., 2019) and these results

274

indicated that the lipid metabolism response to PAHs exposure may contribute to the

275

MetS progression. In addition, other NHANES analyses have shown that increased

276

concentrations of 2-FLU were significantly associated with a higher prevalence of the

277

MetS (Hu et al., 2015) and fluorene metabolite also showed a marginally significant

278

linear trend with self-report CVD (Xu et al., 2010; Xu et al., 2013). It was consistent 13

279

with the results of our study. And the sensitivity analyses did not appreciably alter the

280

results, it indicates that the statistical analysis in this study has a reasonable control

281

over confounding factors and the analysis results are highly reliable. Meanwhile, the

282

associations of urinary PAHs with obesity and the expression of a number of

283

obesity-related cardiometabolic health risk factors have been reported in the literature

284

(Ranjbar et al., 2015), the positive dose-dependent association between obesity and

285

2-PHE was completely consistent with our results.

286

Exposure to cigarette smoke has negative effects on lipid metabolism and

287

oxidative stress status, even caused early glucose intolerance, the changes caused by

288

cigarette smoke exposure can trigger the earlier onset of metabolic disorders

289

associated with obesity, such as diabetes and metabolic syndrome (Damasceno et al.,

290

2013; Lietz et al., 2013). The evidence has been implicated the etiological role of

291

NADPH oxidase (NOX) in smoking-induced CVD, the dysregulations of reactive

292

oxygen species (ROS) generation and metabolism mainly contribute to the

293

development of diverse CVDs (Kim et al., 2014). Experimental studies have shown

294

that the obese rats exposed to tobacco cigarette smoke presented abnormal HDL-c

295

levels and higher DNA damage, triglycerides, total cholesterol, free fatty acids,

296

VLDL-c, and LDL-c levels compared to control and obese rats exposed to filtered air

297

(Damasceno et al., 2013). Our study found smoking status was associated with a

298

significantly increased risk of prevalence ratio of MetS, high triglycerides and low

299

HDL, which was consistent with the studies that cigarette smoking was associated

300

with MetS risk factors (Titz et al., 2016; Cheng et al., 2019). In summary, the previous 14

301

studies and ours have proved lipid metabolism and MetS were associated with PAHs

302

and smoking status. Exposure to PAHs and tobacco could stimulate a chain of internal

303

responses including xenobiotic metabolism, oxidative stress and DNA damages

304

(Shimada and Fujii-Kuriyama, 2004; Kuang et al., 2013; Kim et al., 2014). A link

305

between PAHs concentrations and MetS has also been reported by a NHANES study

306

(Hu et al., 2015), and consistent results were observed in the subgroup analysis among

307

non-smokers. These findings suggest that environmental exposure to PAHs

308

independent of cigarette smoking is associated with the increased prevalence of MetS.

309

Consequently, the impacts of smoking and PAHs on MetS were important and both

310

increase the prevalence. Meanwhile, there were more male in coke oven workers

311

whose lifestyles were correlated strongly with tobacco and alcohol. As expected, we

312

found there was an interaction effect of smoking and PAHs exposure on the

313

prevalence of MetS in this study and it’s consistent with the previously reported

314

studies.

315

The consistent inverse associations of 3-FLU with abdominal obesity and 9-PHE

316

with high glucose was unexpected. Meanwhile, 9-PHE was negative association with

317

lower level HDL. The underlying mechanisms for the association between PAHs

318

exposure and metabolic damage remain an open question. We have not found the

319

relationship between 3-FLU, 9-PHE and MetS, and as their concentration increased,

320

the change of related MetS components was marginally significant. The causes of

321

MetS are complex. Lipid metabolism is not only related to environmental pollution

322

and living habits, but also has a great relationship with genetic factors, 15

323

insulin-resistant and inflammatory factors (Maintinguer Norde et al., 2018). Our

324

research should then take these factors into account and make a further exploration.

325

Moreover, we monitored the environmental levels of PAH exposure and found that

326

the sum of PAH in the air of the worker places was dramatically lower than another

327

study (Kuang et al., 2013) in coke oven workers (0.38 mg/m3 vs. 1.13 mg/m3 for the

328

non-coke-oven workers, 1.45 mg/m3 for the coke-oven workers vs. 11.08 mg/m3 for

329

workers in the side and bottom of the coke oven, and 90.30 mg/m3 for workers at the

330

top of the coke oven). The concentration of this study is also lower than some of the

331

other coke-oven plant (Fu et al., 2018) and chimney sweeps workers (Alhamdow et al.,

332

2019). The difference in urinary PAH metabolites may be caused by lifestyle

333

behaviors, air pollution, laboratory methods and regional differences. In addition,

334

Table S1 showed that good quality control data for measurement of the urinary PAH

335

metabolites in present study.

336

Several limitations and strengths should be acknowledged in the current study.

337

First, a cross sectional study can show associations between PAHs exposure, smoking

338

status and the prevalence of MetS, but cannot assess the time order. Future studies

339

should consider repeatedly measuring PAHs exposure to better elucidate longitudinal

340

associations with MetS. Second, other pollutants, such as metals, nitrogen dioxide,

341

and ozone, are also present in the coke oven emissions and can increase the risk of

342

cardiometabolic health. Likewise, measurement of tobacco metabolites in urine as

343

better long-term markers of tobacco exposures, a more accurate picture of tobacco

344

exposure may be obtained. Then, we could further analyse the dose-response effects 16

345

of tobacco metabolites and PAH metabolites on MetS. And, we found no significant

346

effect of eating habits on MetS in this study, perhaps some workers suffer from

347

chronic metabolic diseases had changed their inherent eating habits. In addition, the

348

association between MetS and PAH metabolites was conducted in China occupational

349

workers, the results may not be generalized to populations with other ethnics and

350

different PAH exposure ranges. Despite these drawbacks, the current analysis has

351

several strengths. The modified Poisson regression models were used to analyse the

352

association between PAHs, smoking status and MetS (Thompson et al., 1998; Zou,

353

2004; Zou and Donner, 2013). If there are many covariate strata, the results for

354

estimating the PR by logistic regression may be quite cumbersome. Not only will this

355

conversion method provide invalid confidence limits, but it will also produce

356

inconsistent estimates for the relative risk. However, proportional hazards regression

357

to directly estimate the PR, but with wider or less precise interval estimates. Poisson

358

regression is usually regarded as an appropriate approach for analyzing rare events.

359

When Poisson regression is applied to binomial data, the error for the estimated

360

relative risk will be overestimated. However, this problem may be rectified by using a

361

robust error variance procedure known as sandwich estimation, thus leading to a

362

technique that I refer to as modified Poisson regression. Compared with application of

363

binomial regression, the modified Poisson regression procedure has no difficulty with

364

converging, and it provides results very similar to those obtained by using the

365

Mantel-Haenszel procedure when the covariate of interest is categorical. We also

366

tested the results in Logistic regression models, and the results were consistent with 17

367

modified Poisson regression models. However, the modified Poisson regression

368

models were more suitable for studying MetS. As we know, the etiology of MetS is

369

very complex. Meanwhile, there were 68.7% male workers in our study, 61.1% of

370

them were smokers, and the proportion of male in patients with MetS accounts for

371

88.8%. However, there was little study to evaluate environmental pollution and poor

372

lifestyles interactions on MetS among coke oven workers. In summary, this work

373

provides insights on PAHs exposures and tobacco exposure with adverse lipid

374

metabolic health in Chinese coke oven workers. Future prospective studies are needed

375

to confirm whether these findings represent causal associations.

376

5. Conclusions

377

Elevated exposure to PAHs can increase the prevalence ratio of MetS and this

378

effect can be modified by smoking status. These results may add potential evidence

379

for environment-lifestyle interactions on MetS.

380

Funding information

381

This work was supported by National Nature Science Foundation of China [No.

382

81273041 and 30901180] and Natural Science Foundation of Shanxi Province of

383

China [No. 201701D121146].

384

Conflicts of interest

385

The authors declare that they have no actual or potential competing financial

386

interest.

387

Acknowledgements

388

We are grateful for General Hospital of Taiyuan Iron & Steel (Group) Co., 18

389

Ltd.and Center of Occupational Disease Prevention of Xishan Coal Electricity (Group)

390

Co., Ltd. for their helps in collecting biological sample and interviewing the study

391

population participants.

392

19

393

Reference

394

Alberti, K.G., Eckel, R.H., Grundy, S.M., Zimmet, P.Z., Cleeman, J.I., Donato, K.A.,

395

Fruchart, J.C., James, W.P., Loria, C.M., Smith, S.C., Jr., 2009. Harmonizing the

396

metabolic syndrome: a joint interim statement of the International Diabetes

397

Federation Task Force on Epidemiology and Prevention; National Heart, Lung,

398

and Blood Institute; American Heart Association; World Heart Federation;

399

International Atherosclerosis Society; and International Association for the Study

400

of Obesity. Circulation 120, 1640-1645.

401

Alhamdow, A., Lindh, C., Albin, M., Gustavsson, P., Tinnerberg, H., Broberg, K.,

402

2019. Cardiovascular disease-related serum proteins in workers occupationally

403

exposed to polycyclic aromatic hydrocarbons. Toxicological sciences : an official

404

journal of the Society of Toxicology.

405

Burstyn, I., Kromhout, H., Partanen, T., Svane, O., Langard, S., Ahrens, W.,

406

Kauppinen, T., Stucker, I., Shaham, J., Heederik, D., Ferro, G., Heikkila, P.,

407

Hooiveld, M., Johansen, C., Randem, B.G., Boffetta, P., 2005. Polycyclic

408

aromatic hydrocarbons and fatal ischemic heart disease. Epidemiology 16,

409

744-750.

410

Chen, H., Goldberg, M.S., Villeneuve, P.J., 2008. A systematic review of the relation

411

between long-term exposure to ambient air pollution and chronic diseases.

412

Reviews on environmental health 23, 243-297.

413 414

Cheng, E., Burrows, R., Correa, P., Guichapani, C.G., Blanco, E., Gahagan, S., 2019. Light smoking is associated with metabolic syndrome risk factors in Chilean 20

415 416

young adults. Acta diabetologica 56, 473-479. Coogan, P.F., White, L.F., Jerrett, M., Brook, R.D., Su, J.G., Seto, E., Burnett, R.,

417

Palmer, J.R., Rosenberg, L., 2012. Air pollution and incidence of hypertension

418

and diabetes mellitus in black women living in Los Angeles. Circulation 125,

419

767-772.

420

Damasceno, D.C., Sinzato, Y.K., Bueno, A., Dallaqua, B., Lima, P.H., Calderon, I.M.,

421

Rudge, M.V., Campos, K.E., 2013. Metabolic profile and genotoxicity in obese

422

rats exposed to cigarette smoke. Obesity (Silver Spring, Md.) 21, 1596-1601.

423

Fu, W., Chen, Z., Bai, Y., Wu, X., Li, G., Chen, W., Wang, G., Wang, S., Li, X., He, M.,

424

Zhang, X., Wu, T., Guo, H., 2018. The interaction effects of polycyclic aromatic

425

hydrocarbons exposure and TERT- CLPTM1L variants on longitudinal telomere

426

length shortening: A prospective cohort study. Environmental pollution (Barking,

427

Essex : 1987) 242, 2100-2110.

428 429 430

Grundy, S.M., 2016. Metabolic syndrome update. Trends Cardiovasc Med 26, 364-373. Hu, H., Kan, H., Kearney, G.D., Xu, X., 2015. Associations between exposure to

431

polycyclic aromatic hydrocarbons and glucose homeostasis as well as metabolic

432

syndrome in nondiabetic adults. The Science of the total environment 505,

433

56-64.

434

Keith, R.J., Al Rifai, M., Carruba, C., De Jarnett, N., McEvoy, J.W., Bhatnagar, A.,

435

Blaha, M.J., Defilippis, A.P., 2016. Tobacco Use, Insulin Resistance, and Risk of

436

Type 2 Diabetes: Results from the Multi-Ethnic Study of Atherosclerosis. PloS 21

437 438 439 440

one 11, e0157592. Kim, M., Han, C.H., Lee, M.Y., 2014. NADPH oxidase and the cardiovascular toxicity associated with smoking. Toxicological research 30, 149-157. Kuang, D., Zhang, W., Deng, Q., Zhang, X., Huang, K., Guan, L., Hu, D., Wu, T.,

441

Guo, H., 2013. Dose-response relationships of polycyclic aromatic hydrocarbons

442

exposure and oxidative damage to DNA and lipid in coke oven workers.

443

Environmental science & technology 47, 7446-7456.

444

Li, F., Xiang, B., Jin, Y., Li, C., Li, J., Ren, S., Huang, H., Luo, Q., 2019.

445

Dysregulation of lipid metabolism induced by airway exposure to polycyclic

446

aromatic hydrocarbons in C57BL/6 mice. Environmental pollution (Barking,

447

Essex : 1987) 245, 986-993.

448

Lietz, M., Berges, A., Lebrun, S., Meurrens, K., Steffen, Y., Stolle, K., Schueller, J.,

449

Boue, S., Vuillaume, G., Vanscheeuwijck, P., Moehring, M., Schlage, W., De

450

Leon, H., Hoeng, J., Peitsch, M., 2013. Cigarette-smoke-induced atherogenic

451

lipid profiles in plasma and vascular tissue of apolipoprotein E-deficient mice are

452

attenuated by smoking cessation. Atherosclerosis 229, 86-93.

453

Liu, Y., Zhang, H., Zhang, H., Niu, Y., Fu, Y., Nie, J., Yang, A., Zhao, J., Yang, J.,

454

2018. Mediation effect of AhR expression between polycyclic aromatic

455

hydrocarbons exposure and oxidative DNA damage among Chinese occupational

456

workers. Environmental pollution (Barking, Essex : 1987) 243, 972-977.

457

Lu, L., Johnman, C., McGlynn, L., Mackay, D.F., Shiels, P.G., Pell, J.P., 2017.

458

Association between exposure to second-hand smoke and telomere length: 22

459

cross-sectional study of 1303 non-smokers. International journal of

460

epidemiology 46, 1978-1984.

461

Lu, L., Mackay, D.F., Pell, J.P., 2014. Meta-analysis of the association between

462

cigarette smoking and peripheral arterial disease. Heart (British Cardiac Society)

463

100, 414-423.

464

Maintinguer Norde, M., Oki, E., Ferreira Carioca, A.A., Teixeira Damasceno, N.R.,

465

Fisberg, R.M., Lobo Marchioni, D.M., Rogero, M.M., 2018. Influence of IL1B,

466

IL6 and IL10 gene variants and plasma fatty acid interaction on metabolic

467

syndrome risk in a cross-sectional population-based study. Clinical nutrition

468

(Edinburgh, Scotland) 37, 659-666.

469

Ranjbar, M., Rotondi, M.A., Ardern, C.I., Kuk, J.L., 2015. Urinary Biomarkers of

470

Polycyclic Aromatic Hydrocarbons Are Associated with Cardiometabolic Health

471

Risk. PloS one 10, e0137536.

472

Scinicariello, F., Buser, M.C., 2014. Urinary polycyclic aromatic hydrocarbons and

473

childhood obesity: NHANES (2001-2006). Environmental health perspectives

474

122, 299-303.

475

Shimada, T., Fujii-Kuriyama, Y., 2004. Metabolic activation of polycyclic aromatic

476

hydrocarbons to carcinogens by cytochromes P450 1A1 and 1B1. Cancer science

477

95, 1-6.

478

Simkhovich, B.Z., Kleinman, M.T., Kloner, R.A., 2008. Air pollution and

479

cardiovascular injury epidemiology, toxicology, and mechanisms. Journal of the

480

American College of Cardiology 52, 719-726. 23

481

Slagter, S.N., van Vliet-Ostaptchouk, J.V., Vonk, J.M., Boezen, H.M., Dullaart, R.P.,

482

Kobold, A.C., Feskens, E.J., van Beek, A.P., van der Klauw, M.M., Wolffenbuttel,

483

B.H., 2013. Associations between smoking, components of metabolic syndrome

484

and lipoprotein particle size. BMC medicine 11, 195.

485

Thompson, M.L., Myers, J.E., Kriebel, D., 1998. Prevalence odds ratio or prevalence

486

ratio in the analysis of cross sectional data: what is to be done? Occupational and

487

environmental medicine 55, 272-277.

488

Titz, B., Boue, S., Phillips, B., Talikka, M., Vihervaara, T., Schneider, T., Nury, C.,

489

Elamin, A., Guedj, E., Peck, M.J., Schlage, W.K., Cabanski, M., Leroy, P.,

490

Vuillaume, G., Martin, F., Ivanov, N.V., Veljkovic, E., Ekroos, K., Laaksonen, R.,

491

Vanscheeuwijck, P., Peitsch, M.C., Hoeng, J., 2016. Effects of Cigarette Smoke,

492

Cessation, and Switching to Two Heat-Not-Burn Tobacco Products on Lung

493

Lipid Metabolism in C57BL/6 and Apoe-/- Mice-An Integrative Systems

494

Toxicology Analysis. Toxicological sciences : an official journal of the Society of

495

Toxicology 149, 441-457.

496

Xu, X., Cook, R.L., Ilacqua, V.A., Kan, H., Talbott, E.O., Kearney, G., 2010. Studying

497

associations between urinary metabolites of polycyclic aromatic hydrocarbons

498

(PAHs) and cardiovascular diseases in the United States. The Science of the total

499

environment 408, 4943-4948.

500

Xu, X., Hu, H., Kearney, G.D., Kan, H., Sheps, D.S., 2013. Studying the effects of

501

polycyclic aromatic hydrocarbons on peripheral arterial disease in the United

502

States. The Science of the total environment 461-462, 341-347. 24

503 504

Zou, G., 2004. A modified poisson regression approach to prospective studies with binary data. American journal of epidemiology 159, 702-706.

505

Zou, G.Y., Donner, A., 2013. Extension of the modified Poisson regression model to

506

prospective studies with correlated binary data. Stat Methods Med Res 22,

507

661-670.

508

25

509

Figure Legends

510

Figure 1. The associations between PAH metabolites and MetS components

511

among the study population (n=682)

512

T = tertile

513

Adjusted prevalence ratios (95%) for metabolic syndrome and components conditions

514

by PAH metabolites

515

Adjusted for sex, age, smoking, drinking, cooking fumes, eating habits, BMI (except

516

for abdominal obesity), 2-NAP, 1-NAP, 3-FLU, 2-FLU, 2-PHE, 9-PHE, 1-PHE,

517

1-PYR and 9-BAP

518

Figure 2. The associations between lifestyles and MetS and components among

519

the study population (n=682)

520

Adjusted prevalence ratios (95%) for metabolic syndrome component conditions by

521

lifestyles (smoking, drinking and eating habits)

522

Adjusted for sex, age, smoking, drinking, cooking fumes, eating habits, BMI (except

523

for abdominal obesity), 2-NAP, 1-NAP, 3-FLU, 2-FLU, 2-PHE, 9-PHE, 1-PHE,

524

1-PYR and 9-BAP

525

Figure 3. Adjusted prevalence ratios of MetS and low HDL for different

526

combinations of smoking status and the levels of PAH metabolites (n=682)

527

T = tertile

528

Data were presented as prevalence ratio (PR)

529

Adjusted for sex, age, smoking, drinking, cooking fumes, eating habits, BMI, 2-NAP,

530

1-NAP, 3-FLU, 2-FLU, 2-PHE, 9-PHE, 1-PHE, 1-PYR and 9-BAP 26

531

(A) MetS & 1-NAP; (B) MetS & 2-FLU; (C) low HDL & 1-NAP

27

Table 1. Basic characteristics of participants by Metabolic Syndrome among 682 occupational workersa Variable

Metabolic Syndrome Total (n=682)

Yes (n=100)

No (n=582)

38 (31 - 48)

45.5 (33 - 50)

35 (30 - 47)

< 32

218 (23.5)

197 (25.3)

21 (14.3)

32 - 46

222 (31.0)

193 (31.9)

29 (26.4)

≥ 46

242 (45.4)

192 (42.8)

50 (59.2)

male

458 (68.7)

89 (88.8)

369 (64.9)

female

224 (31.3)

11 (11.2)

213 (35.1)

no

401 (57.7)

39 (40.3)

362 (61.0)

yes

281 (42.3)

61 (59.7)

220 (39.0)

no

494 (70.5)

62 (61.5)

432 (72.2)

yes

188 (29.5)

38 (38.5)

150 (27.8)

few

604 (87.7)

80 (78.9)

524 (89.3)

much

77 (12.2)

20 (21.1)

57 (10.5)

≤ 25

373 (54.2)

13 (12.4)

360 (62.1)

> 25

309 (45.8)

87 (87.6)

222 (37.9)

light

152 (22.5)

24 (24.8)

128 (22.0)

medium

316 (44.5)

40 (39.2)

276 (45.5)

heavy

214 (33.1)

36 (36.0)

178 (32.5)

< 2.5

223 (32.4)

31 (29.2)

192 (33)

≥ 2.5

459 (67.6)

69 (70.8)

390 (67)

Age (years)

P valueb

<0.001

Age (years)

0.003

Sex <0.001

Smoking <0.001

Drinking 0.012

Cooking fumes 0.003

BMI <0.001

Eating habits Salt

0.374

Vegetables (kg/week)

Fruits (g/week)

0.695

never

10 (1.7)

-

10 (2)

< 750

294 (43.6)

53 (51.8)

241 (42.1)

≥ 750

378 (54.6)

47 (48.2)

331 (55.9)

2-NAP

0.30 (0.13 - 0.69)

0.31 (0.15 - 1.04)

0.29 (0.13 - 0.68)

0.275

1-NAP

0.03 (0.01 - 0.06)

0.03 (0.01 - 0.07)

0.03 (0.01 - 0.06)

0.065

3-FLU

0.03 (0.01 - 0.06)

0.04 (0.02 - 0.06)

0.03 (0.01 - 0.05)

0.264

2-FLU

0.15 (0.07 - 0.29)

0.20 (0.12 - 0.34)

0.14 (0.07 - 0.28)

0.002

2-PHE

0.20 (0.12 - 0.32)

0.22 (0.12 - 0.42)

0.19 (0.12 - 0.32)

0.369

9-PHE

0.11 (0.07 - 0.18)

0.13 (0.07 - 0.18)

0.11 (0.07 - 0.18)

0.275

1-PHE

0.06 (0.03 - 0.11)

0.06 (0.03 - 0.12)

0.06 (0.03 - 0.11)

0.644

1-PYR

0.09 (0.06 - 0.16)

0.11 (0.06 - 0.22)

0.09 (0.06 - 0.16)

0.080

9-BAP

0.011 (0.003 - 0.029) 0.011 (0.004 - 0.029) 0.011 (0.003 - 0.029)

0.053

PAHs internal exposure biomarker (ng/ml)

0.727

Concentration unit of 2-NAP, 1-NAP 3-FLU, 2-FLU, 2-PHE, 9-PHE, 1-PHE, 1-PYE and 9-BAP is ng/ml. a : Data were presented as n (%) or Med (25th - 75th). b : P-values were calculated from Chi-square test for categorical variables and Kruskal-Wallis H test for numerical variables.

Table 2. The associations between PAH metabolites and Metabolic Syndrome among the study population (n=682). PAH metabolites (ng/ml) 2-NAP 1st tertile (< 0.17) 2nd tertile (0.17 - 0.50) 3rd tertile (≥ 0.50)

Metabolic Syndrome (n=100) PR* (95% CI) Unadjusted

Model1a

Model2b

1.00 ( Reference) 1.02 (0.96 - 1.08) 1.02 (0.96 - 1.08)

1.00 ( Reference) 1.00 (0.94 - 1.06) 0.96 (0.89 - 1.03)

1.00 ( Reference) 1.00 (0.94 - 1.05) 0.96 (0.90 - 1.03) 0.352

Ptrend 1-NAP 1st tertile (< 0.01) 2nd tertile (0.01 - 0.05) 3rd tertile (≥ 0.05)

1.00 ( Reference) 1.02 (0.97 - 1.08) 1.06 (1.00 - 1.12)

1.00 ( Reference) 1.02 (0.96 - 1.09) 1.08 (1.01 - 1.17)

Ptrend 3-FLU 1st tertile (< 0.02) 2nd tertile (0.02 - 0.04) 3rd tertile (≥ 0.04)

0.026 1.00 ( Reference) 1.03 (0.97 - 1.09) 1.02 (0.97 - 1.08)

1.00 ( Reference) 0.98 (0.92 - 1.04) 0.91 (0.85 - 0.99)

Ptrend 2-FLU 1st tertile (< 0.09) 2nd tertile (0.09 - 0.22) 3rd tertile (≥ 0.22)

1.00 ( Reference) 1.06 (1.01 - 1.12) 1.10 (1.04 - 1.16)

1.00 ( Reference) 1.06 (1.01 - 1.12) 1.12 (1.04 - 1.21)

1.00 ( Reference) 0.99 (0.94 - 1.05) 1.01 (0.96 - 1.07)

1.00 ( Reference) 0.97 (0.92 - 1.03) 0.97 (0.89 - 1.05)

Ptrend

1.00 ( Reference) 0.98 (0.92 - 1.04) 0.96 (0.89 - 1.04) 0.289

1.00 ( Reference) 0.98 (0.93 - 1.03) 1.03 (0.97 - 1.09)

1.00 ( Reference) 0.96 (0.90 - 1.02) 0.99 (0.91 - 1.06)

Ptrend 1-PHE 1st tertile (< 0.04) 2nd tertile (0.04 - 0.09) 3rd tertile (≥ 0.09)

1.00 ( Reference) 1.05 (1.00 - 1.11) 1.10 (1.02 - 1.18) 0.042

Ptrend 9-PHE 1st tertile (< 0.08) 2nd tertile (0.08 - 0.15) 3rd tertile (≥ 0.15)

1.00 ( Reference) 0.98 (0.93 - 1.04) 0.94 (0.87 - 1.01) 0.063

Ptrend 2-PHE 1st tertile (< 0.14) 2nd tertile (0.14 - 0.27) 3rd tertile (≥ 0.27)

1.00 ( Reference) 1.03 (0.97 - 1.09) 1.09 (1.01 - 1.17)

1.00 ( Reference) 0.96 (0.91 - 1.02) 0.99 (0.92 - 1.06) 0.859

1.00 ( Reference) 0.97 (0.92 - 1.03) 0.99 (0.94 - 1.05)

1.00 ( Reference) 0.98 (0.92 - 1.05) 0.99 (0.91 - 1.07)

1.00 ( Reference) 0.97 (0.91 - 1.03) 0.99 (0.92 - 1.07) 0.981

1-PYR 1st tertile (< 0.07) 2nd tertile (0.07 - 0.14) 3rd tertile (≥ 0.14)

1.00 ( Reference) 0.98 (0.93 - 1.03) 1.04 (0.98 - 1.10)

1.00 ( Reference) 0.97 (0.91 - 1.03) 0.99 (0.92 - 1.07)

1.00 ( Reference) 0.97 (0.92 - 1.03) 1.00 (0.93 - 1.07) 0.829

Ptrend

* A modified Poisson Regression was conducted. a : adjusted for sex, age, education. b : additional adjust for smoking, drinking, physical activities, eating habits, BMI, 2-NAP, 1-NAP, 3-FLU, 2-FLU, 2-PHE, 9-PHE, 1-PHE, 1-PYR, 9-BAP.

Highlights: 

Relationship of PAH exposure and MetS was analysed by modified Poisson regression.



PAHs exposure can increase the prevalence ratio of MetS in coke oven workers.



Smoking modifies the prevalence ratio of MetS induced by PAHs exposure.

Jin Yang: Funding acquisition, Conceptualization, Supervision, Project administration Bin Zhang: Writing- Original draft, Methodology and Formal analysis. Baolong Pan: Writing- Original draft, Methodology and Formal analysis. Xinyu Zhao: Investigation, Data Curation, Resources Ye Fu: Investigation, Data Curation, Validation Xuejing Li: Investigation, Resources Aimin Yang: Writing - Review & Editing, Software and Methodology Qiang Li: Investigation Jun Dong: Investigation Jisheng Nie: Project administration All authors have read and approved the final text.

Declaration of Interest Statement The authors declare that they have no actual or potential competing financial interest.