2.8 MB Sample
2.8 MB Sample
2.8 MB Sample
DOI: 10.1289/ehp.1509834
Note to readers with disabilities: EHP strives to ensure that all journal content is accessible to all
readers. However, some figures and Supplemental Material published in EHP articles may not conform to
508 standards due to the complexity of the information being presented. If you need assistance accessing
journal content, please contact [email protected]. Our staff will work with you to assess and meet
your accessibility needs within 3 working days.
Supplemental Material
Xiaoyan Zhu, Jun Li, Siyun Deng, Kuai Yu, Xuezhen Liu, Qifei Deng, Huizhen Sun, Xiaomin
Zhang, Meian He, Huan Guo, Weihong Chen, Jing Yuan, Bing Zhang, Dan Kuang, Xiaosheng
He, Yansen Bai, Xu Han, Bing Liu, Xiaoliang Li, Liangle Yang, Haijing Jiang, Yizhi Zhang, Jie
Hu, Longxian Cheng, Xiaoting Luo, Wenhua Mei, Zhiming Zhou, Shunchang Sun, Liyun Zhang,
Chuanyao Liu, Yanjun Guo, Zhihong Zhang, Frank B. Hu, Liming Liang, and Tangchun Wu
Table of Contents
Table S1. The 318 identified CpGs associated with cigarette smoking in the genome-wide
meta-analysis (FDR < 0.05).
Table S2. The effect size (s.e.) of smoking on the 318 CpGs with inverse normal transformed
methylation values calculated in different groups.
Table S3. The median methylation values (inter quartile range) of the 318 smoking-related CpGs
in never, former and current smokers.
Table S4. Methylation-gene expression correlations of the smoking-related CpGs in SY (n =
144).
Table S5. Contribution rates of smoking on urinary PAHs metabolites calculated in male
individuals from the WHZH cohort (n = 206).
Table S6. Associations between the smoking-related CpGs and urinary 2-hydroxynaphthalene
levels in males from the WHZH Cohort and the Coke Oven Cohort.
Figure S1. The flowchart of the study.
References
2
Table S1. The 318 identified CpGs associated with cigarette smoking in the genome-wide meta-analysis (FDR < 0.05).
13
Table S2. The effect size (s.e.) of smoking on the 318 CpGs with inverse normal transformed methylation values calculated in
different groups.
24
cg03217253 ZNF677 0.213 (0.092) 0.023 0.172 (0.106) 0.107 0.253 (0.060) 3.10E-05 0.150 (0.105) 0.155
cg21733502 ZSCAN5B -0.268 (0.104) 0.012 -0.166 (0.103) 0.110 -0.155 (0.071) 0.030 -0.256 (0.111) 0.024
cg05007126 SDCBP2 -0.062 (0.083) 0.456 -0.100 (0.074) 0.184 -0.164 (0.047) 5.15E-04 -0.169 (0.079) 0.035
cg13724496 BMP2 0.048 (0.107) 0.655 0.248 (0.102) 0.017 0.240 (0.059) 6.04E-05 0.102 (0.103) 0.324
cg15207742 RIMS4 -0.055 (0.106) 0.602 0.310 (0.095) 0.002 0.275 (0.066) 4.17E-05 0.199 (0.092) 0.033
cg20344344 - -0.232 (0.090) 0.012 -0.024 (0.085) 0.783 -0.161 (0.051) 0.002 -0.216 (0.081) 0.009
cg12303084 ZMYND8 -0.103 (0.070) 0.144 -0.085 (0.056) 0.135 -0.110 (0.035) 0.002 -0.146 (0.072) 0.046
cg07339236 ATP9A -0.274 (0.086) 0.002 -0.270 (0.100) 0.009 -0.248 (0.060) 5.20E-05 -0.240 (0.101) 0.021
cg26736540 TFAP2C 0.012 (0.088) 0.894 0.136 (0.106) 0.203 0.187 (0.059) 0.002 0.367 (0.104) 6.65E-04
cg06135139 BHLHE23 0.152 (0.108) 0.162 -0.079 (0.122) 0.516 0.315 (0.068) 5.41E-06 0.305 (0.117) 0.011
cg15892280 ETS2 -0.117 (0.058) 0.049 -0.071 (0.042) 0.097 -0.142 (0.040) 4.38E-04 -0.068 (0.057) 0.236
cg23110422 ETS2 -0.357 (0.101) 7.04E-04 -0.164 (0.112) 0.147 -0.229 (0.071) 0.001 -0.282 (0.102) 0.007
cg10204884 PCBP3 -0.201 (0.095) 0.037 -0.223 (0.081) 0.007 -0.139 (0.058) 0.017 -0.122 (0.089) 0.176
cg14785479 SCARF2 0.159 (0.067) 0.020 0.094 (0.051) 0.071 0.121 (0.039) 0.002 0.135 (0.052) 0.012
cg08548559 PIK3IP1 -0.231 (0.083) 0.007 -0.197 (0.079) 0.014 -0.077 (0.052) 0.144 -0.234 (0.087) 0.009
cg02532700 NCF4 -0.251 (0.109) 0.024 -0.210 (0.112) 0.063 -0.182 (0.063) 0.004 -0.299 (0.114) 0.010
cg26364091 CHADL -0.180 (0.108) 0.099 -0.266 (0.113) 0.021 -0.180 (0.065) 0.006 -0.195 (0.108) 0.075
cg06008724 PHF21B 0.040 (0.105) 0.706 0.226 (0.100) 0.027 0.220 (0.056) 1.04E-04 0.115 (0.080) 0.156
cg04198308 FAM19A5 0.270 (0.098) 0.007 0.221 (0.112) 0.052 0.223 (0.064) 6.09E-04 0.018 (0.103) 0.862
25
Table S3. The median methylation values (inter quartile range) of the 318 smoking-related CpGs in never, former and current
smokers.
36
Table S4. Methylation-gene expression correlations of the smoking-related CpGs in SY (n = 144).
CpG information
Transcripts information
CpG-expression probe association
Chr Position CpG Gene Relation to Gene
Probe Gene
Effect s.e. p
1 11865352 cg05228408 MTHFR TSS1500
ILMN_1734830 MTHFR
7.733 2.148 4.41E-04
1 11865352 cg05228408 CLCN6 TSS1500
ILMN_1763036 CLCN6
7.137 2.156 0.001
1 16301562 cg04211179 ZBTB17 5'UTR
ILMN_1711048 ZBTB17
7.856 1.711 9.73E-06
1 21617442 cg26348226 ECE1 TSS1500
ILMN_1672174 ECE1
-14.337 1.926 9.11E-12
1 25254746 cg10951873 RUNX3 Body
ILMN_1787461 RUNX3
-26.344 4.769 1.57E-07
1 28521540 cg20460771 PTAFR TSS1500
ILMN_1746836 PTAFR
13.025 1.395 2.07E-16
1 90289611 cg20146909 LRRC8D 5'UTR
ILMN_1763409 LRRC8D
6.287 2.812 0.027
1 153746211 cg08129092 INTS3 3'UTR
ILMN_1756086 INTS3
-16.351 2.352 1.27E-10
1 154299179 cg06811467 ATP8B2 TSS1500
ILMN_1782057 ATP8B2
11.790 2.625 1.47E-05
1 154379696 cg09257526 IL6R Body
ILMN_1696394 IL6R
-16.010 2.494 1.98E-09
1 154379696 cg09257526 IL6R Body
ILMN_1754753 IL6R
-12.405 2.627 5.61E-06
2 231790037 cg16382047 GPR55 TSS200
ILMN_1723274 GPR55
7.689 1.974 1.52E-04
2 231790777 cg19827923 GPR55 TSS1500
ILMN_1723274 GPR55
-17.588 3.253 2.70E-07
3 49394622 cg18642234 GPX1 3'UTR
ILMN_1749662 GPX1
4.120 1.963 0.038
3 111260756 cg05655806 CD96 TSS200
ILMN_2415786 CD96
-7.176 1.715 5.05E-05
3 111260756 cg05655806 CD96 TSS200
ILMN_1711573 CD96
-9.435 1.646 5.89E-08
3 111260783 cg04039397 CD96 TSS200
ILMN_2415786 CD96
-6.825 1.525 1.57E-05
3 111260783 cg04039397 CD96 TSS200
ILMN_1711573 CD96
-9.026 1.449 5.36E-09
3 126242973 cg18165852 CHST13 TSS1500
ILMN_1734707 CHST13
-8.940 2.845 0.002
3 188425256 cg25853622 LPP Body
ILMN_1651254 LPP
-5.155 1.120 9.25E-06
3 194119885 cg13185177 GP5 5'UTR
ILMN_1780622 GP5
3.062 1.533 0.048
4 56813860 cg26542660 CEP135 TSS1500
ILMN_1693766 CEP135
-7.353 3.078 0.018
5 138861241 cg04232128 TMEM173 Body
ILMN_2145116 TMEM173
17.222 1.732 5.84E-18
6 35696870 cg25114611 FKBP5 TSS1500
ILMN_1778444 FKBP5
-6.479 2.793 0.022
37
6 36650733 cg01955533 CDKN1A 5'UTR
ILMN_1784602 CDKN1A
5.044 1.837 0.007
7 2847575 cg18446336 GNA12 Body
ILMN_2216157 GNA12
3.496 1.181 0.004
7 5457225 cg09022230 TNRC18 Body
ILMN_1678329 TNRC18
-9.621 1.664 4.61E-08
7 45002736 cg04180046 MYO1G Body
ILMN_1692295 MYO1G
5.040 1.758 0.005
7 51384528 cg09613161 COBL TSS200
ILMN_1711888 COBL
-2.977 0.990 0.003
7 51384609 cg04016086 COBL TSS200
ILMN_1711888 COBL
-3.274 1.548 0.036
7 56147257 cg19956914 SUMF2 Body
ILMN_1685371 SUMF2
2.725 1.084 0.013
7 100465833 cg22851200 TRIP6 Body
ILMN_1713990 TRIP6
-5.725 1.636 6.24E-04
7 116786606 cg22619824 ST7 Body
ILMN_1702175 ST7
7.387 1.314 9.94E-08
7 116786606 cg22619824 ST7 Body
ILMN_1707763 ST7
5.346 1.376 1.57E-04
8 31496644 cg18387156 NRG1 TSS1500
ILMN_1737252 NRG1
6.519 2.535 0.011
8 126446923 cg22644321 TRIB1 Body
ILMN_1803811 TRIB1
-18.246 2.262 2.93E-13
8 144576604 cg26361535 ZC3H3 Body
ILMN_1702389 ZC3H3
-3.217 1.597 0.046
8 145012748 cg25325005 PLEC1 Body
ILMN_1744268 PLEC1
-5.925 1.567 2.30E-04
8 145012748 cg25325005 PLEC1 Body
ILMN_1736792 PLEC1
7.896 1.510 6.06E-07
10 49892930 cg15164194 WDFY4 TSS1500
ILMN_3236551 WDFY4
-11.274 1.524 1.17E-11
10 73497514 cg10750182 CDH23 Body
ILMN_1779934 CDH23
-4.980 1.893 0.009
11 2722391 cg07123182 KCNQ1 Body
ILMN_2409520 KCNQ1
9.891 3.372 0.004
11 2722407 cg26963277 KCNQ1OT1 TSS1500
ILMN_3248707 KCNQ1OT1
5.587 2.563 0.031
11 86510915 cg11660018 PRSS23 TSS1500
ILMN_1797776 PRSS23
7.850 1.692 7.94E-06
11 86510998 cg23771366 PRSS23 TSS1500
ILMN_1797776 PRSS23
8.629 2.068 5.25E-05
11 86512100 cg23351584 PRSS23 5'UTR
ILMN_1797776 PRSS23
20.719 4.950 5.00E-05
11 121322456 cg10908953 SORL1 TSS1500
ILMN_1759818 SORL1
-18.779 2.766 2.95E-10
11 121322456 cg10908953 SORL1 TSS1500
ILMN_2060115 SORL1
-21.510 2.635 1.70E-13
12 4384888 cg07066369 CCND2 Body
ILMN_2067656 CCND2
28.294 6.613 3.47E-05
12 4384888 cg07066369 CCND2 Body
ILMN_1667081 CCND2
24.758 6.752 3.48E-04
12 7055657 cg23193870 PTPN6 TSS200
ILMN_1738675 PTPN6
-23.094 6.997 0.001
12 7055657 cg23193870 PTPN6 TSS200
ILMN_1716578 PTPN6
-21.788 7.095 0.003
38
12 11898284 cg07986378 ETV6 Body
ILMN_1789596 ETV6
-6.905 1.419 3.04E-06
12 124016861 cg21618017 RILPL1 Body
ILMN_1805643 RILPL1
-23.846 7.746 0.003
12 129337910 cg06419750 GLT1D1 TSS200
ILMN_1656327 GLT1D1
-4.899 1.281 1.96E-04
13 48987165 cg03646329 RB1 5'UTR
ILMN_1696591 RB1
3.697 1.372 0.008
13 114172890 cg00619505 TMCO3 Body
ILMN_2220739 TMCO3
11.309 2.788 8.24E-05
14 25079357 cg02150910 GZMH TSS1500
ILMN_1731233 GZMH
-14.102 1.462 3.38E-17
14 91818668 cg01055824 CCDC88C Body
ILMN_3248352 CCDC88C
-14.824 3.736 1.15E-04
14 91881497 cg20303561 CCDC88C Body
ILMN_3248352 CCDC88C
-10.949 1.963 1.21E-07
14 93552128 cg05284742 ITPK1 Body
ILMN_1715674 ITPK1
-9.874 2.247 2.18E-05
14 104190678 cg26242531 ZFYVE21 Body
ILMN_2197030 ZFYVE21
-5.361 1.477 3.97E-04
15 38857474 cg03603381 RASGRP1 TSS1500
ILMN_1768958 RASGRP1
5.513 2.317 0.019
15 55513621 cg24687805 RAB27A Body
ILMN_1665859 RAB27A
-20.736 4.910 4.31E-05
15 55513621 cg24687805 RAB27A Body
ILMN_1699878 RAB27A
-27.973 4.637 1.36E-08
15 55513621 cg24687805 RAB27A Body
ILMN_2329773 RAB27A
-25.988 4.564 7.01E-08
15 60959637 cg15451980 RORA Body
ILMN_2322499 RORA
-10.094 1.355 8.70E-12
15 60959637 cg15451980 RORA Body
ILMN_2322498 RORA
-10.437 1.349 1.81E-12
15 60959637 cg15451980 RORA Body
ILMN_1813134 RORA
-8.499 1.430 2.12E-08
15 70387217 cg26971042 TLE3 Body
ILMN_2234412 TLE3
-8.926 2.702 0.001
15 74862662 cg02384859 ARID3B Body
ILMN_1696956 ARID3B
9.939 1.669 2.01E-08
15 90357202 cg23161492 ANPEP 5'UTR
ILMN_1763837 ANPEP
-5.939 2.097 0.005
16 17562960 cg06321596 XYLT1 Body
ILMN_1830462 XYLT1
8.401 2.067 7.97E-05
16 84746995 cg08126789 USP10 Body
ILMN_1721116 USP10
6.098 2.780 0.030
16 89041793 cg04887172 CBFA2T3 Body
ILMN_1657627 CBFA2T3
-6.106 1.268 3.78E-06
16 89408248 cg01107178 ANKRD11 5'UTR
ILMN_2108709 ANKRD11
8.837 1.451 1.03E-08
17 8804279 cg05460226 PIK3R5 Body
ILMN_1681067 PIK3R5
-6.339 1.955 0.001
17 9921982 cg02018337 GAS7 Body
ILMN_1745994 GAS7
-3.668 1.795 0.043
17 27050723 cg18150958 RPL23A Body
ILMN_1788607 RPL23A
-7.159 1.797 1.09E-04
17 27401793 cg18960216 TIAF1 5'UTR
ILMN_2055634 TIAF1
-8.108 2.458 0.001
39
17 46622522 cg22807449 HOXB2 TSS200
ILMN_1810274 HOXB2
8.546 1.931 1.92E-05
17 56082867 cg08591265 SFRS1 3'UTR;Body
ILMN_1795341 SFRS1
-4.897 1.302 2.47E-04
17 58227267 cg05248618 CA4 TSS200
ILMN_1695157 CA4
-11.653 2.872 8.22E-05
17 75371764 cg07827420 SEPT9 5'UTR
ILMN_1769118 SEPT9
-8.857 3.586 0.015
17 75445905 cg07324245 SEPT9 Body
ILMN_1769118 SEPT9
-8.944 1.384 1.58E-09
19 51875451 cg10126923 NKG7 Body
ILMN_1682993 NKG7
9.290 1.462 2.75E-09
19 51876788 cg12916723 NKG7 TSS1500
ILMN_1682993 NKG7
-12.113 2.207 1.84E-07
20 1310884 cg05007126 SDCBP2 TSS1500
ILMN_1705107 SDCBP2
7.264 2.590 0.006
21 40180000 cg15892280 ETS2 5'UTR
ILMN_1720158 ETS2
-13.054 1.607 2.15E-13
22 31686097 cg08548559 PIK3IP1 Body
ILMN_1719986 PIK3IP1
-4.952 1.593 0.002
22 37257404 cg02532700 NCF4 Body
ILMN_2335704 NCF4
-8.136 2.809 0.004
22 37257404 cg02532700 NCF4 Body
ILMN_1757361 NCF4
-6.751 2.766 0.016
22 37257404 cg02532700 NCF4 Body
ILMN_1785005 NCF4
-8.782 2.835 0.002
Abbreviations: Chr, Chromosome; Body, gene body; TSS200, within 200 bps from transcription start site; TSS1500, within 1500 bps
from transcription start site; UTR, untranslated regions.
For each CpG, the significance threshold was defined as 0.05/number of expression probes of the corresponding gene.
40
Table S5. Contribution rates of smoking on urinary PAHs metabolites calculated in male individuals from the WHZH cohort (n =
206).
a
Model 1: Linear regression models with urinary metabolite as dependent variables, age, drinking, BMI, occupation, geographical
regions beadchip and operation date included as independent variables. Concentrations of urinary PAHs metabolites were
ln-transformed. b Model 2: additionally included smoking in Model 1.cThe variance of each urinary OH-PAH metabolite explained by
smoking, calculated as the variance explained by Model 2 minus the variance explained by Model 1.
41
Table S6. Associations between the smoking-related CpGs and urinary 2-hydroxynaphthalene levels in males from the WHZH Cohort
and the Coke Oven Cohorta.
53
Figure S1. The flowchart of the study.
54
References
Dogan MV, Shields B, Cutrona C, Gao L, Gibbons FX, Simons R, et al. 2014. The effect of
smoking on DNA methylation of peripheral blood mononuclear cells from african american
women. BMC Genomics 15:151; doi: 10.1186/1471-2164-15-151 [Online 3 February 2014].
Guida F, Sandanger TM, Castagne R, Campanella G, Polidoro S, Palli D, et al. 2015. Dynamics
of smoking-induced genome-wide methylation changes with time since smoking cessation.
Hum Mol Genet 24:2349-2359.
Philibert RA, Beach SR, Lei MK, Brody GH. 2013. Changes in DNA methylation at the aryl
hydrocarbon receptor repressor may be a new biomarker for smoking. Clin Epigenetics 5:19;
doi: 10.1186/1868-7083-5-19. [Online 10 February 2014].
Shenker NS, Polidoro S, van Veldhoven K, Sacerdote C, Ricceri F, Birrell MA, et al. 2013.
Epigenome-wide association study in the European Prospective Investigation into Cancer
and Nutrition (EPIC-Turin) identifies novel genetic loci associated with smoking. Hum Mol
Genet 22:843-851.
Sun YV, Smith AK, Conneely KN, Chang Q, Li W, Lazarus A, et al. 2013. Epigenomic
association analysis identifies smoking-related DNA methylation sites in African Americans.
Hum Genet 132:1027-1037.
Zeilinger S, Kuhnel B, Klopp N, Baurecht H, Kleinschmidt A, Gieger C, et al. 2013. Tobacco
smoking leads to extensive genome-wide changes in DNA methylation. PLoS One 8:e63812;
doi:10.1371/journal.pone.0063812 [Online 17 May 2013].
55