The impact of methylation quantitative trait loci (mQTLs) on active smoking-related DNA methylation changes

Background Methylation quantitative trait loci (mQTLs) are the genetic variants that may affect the DNA methylation patterns of CpG sites. However, their roles in influencing the disturbances of smoking-related epigenetic changes have not been well established. This study was conducted to address whether mQTLs exist in the vicinity of smoking-related CpG sites (± 50 kb) and to examine their associations with smoking exposure and all-cause mortality in older adults. Results We obtained DNA methylation profiles in whole blood samples by Illumina Infinium Human Methylation 450 BeadChip array of two independent subsamples of the ESTHER study (discovery set, n = 581; validation set, n = 368) and their corresponding genotyping data using the Illumina Infinium OncoArray BeadChip. After correction for multiple testing (FDR), we successfully identified that 70 out of 151 previously reported smoking-related CpG sites were significantly associated with 192 SNPs within the 50 kb search window of each locus. The 192 mQTLs significantly influenced the active smoking-related DNA methylation changes, with percentage changes ranging from 0.01 to 18.96%, especially for the weakly/moderately smoking-related CpG sites. However, these identified mQTLs were not directly associated with active smoking exposure or all-cause mortality. Conclusions Our findings clearly demonstrated that if not dealt with properly, the mQTLs might impair the power of epigenetic-based models of smoking exposure to a certain extent. In addition, such genetic variants could be the key factor to distinguish between the heritable and smoking-induced impact on epigenome disparities. These mQTLs are of special importance when DNA methylation markers measured by Illumina Infinium assay are used for any comparative population studies related to smoking-related cancers and chronic diseases. Electronic supplementary material The online version of this article (doi:10.1186/s13148-017-0387-6) contains supplementary material, which is available to authorized users.


Background
Active smoking has been recognized as a critical lifestyle factor for cardiovascular, respiratory, and neoplastic diseases and contributes to the leading causes of preventable morbidity and mortality [1,2]. DNA methylation, one of the main forms of epigenetic modification, is involved in the pathways of smoking and smoking-induced diseases [3,4]. Previous epigenome-wide association studies (EWASs) based on whole blood samples have successfully discovered an increasing number of tobacco smoking-related CpG sites in various genes, such as AHRR and F2RL3 [5][6][7]. These DNA methylation patterns have been shown to be useful as quantitative biomarkers to reflect both current and lifetime smoking exposure and to enhance the prediction of smoking-related risks [8][9][10][11].
DNA methylation of particular genomic loci might be influenced by neighboring genetic sequence variants [12]. The single nucleotide polymorphisms (SNPs) that are associated with methylation levels of CpG sites are known as methylation quantitative trait loci (mQTLs) [13]. This genetic effect has been determined across different tissues [13][14][15][16] and has been highlighted in several diseases, including neurological disorders, arthritis, and cancer [17][18][19][20][21]. Recently, the mQTLs have been further reported to play a modifying role in the associations between DNA methylation levels at specific CpG sites and environmental exposures. For instance, Zhang et al. identified 238 mQTLs that were associated with 65 alcohol dependence-related CpG sites in African Americans and 305 mQTLs for 44 unique CpG sites in European Americans [22]. In 2016, Gonseth et al. found out that three of the strongest maternal smoking-related CpG sites in newborns were significantly associated with SNPs located in the vicinity of each gene [23]. Thus, these hereditary traits provide a possible mechanism by which methylation patterns could be different under environmental exposures, if the distribution of risk alleles differs between the exposed and the unexposed. In addition, the linkages of epigenetic signatures to genotypes might also further provide more mechanistic evidence on the genetic and environmental risk factors for various forms of diseases [24].
However, such genetic influences have not been well addressed or even overlooked by previous EWASs of active smoking exposure; to our knowledge, no study has so far investigated their contributions to the methylation intensities of active smoking-related CpG sites and smoking-related health outcomes in the general population. Therefore, we conducted a comprehensive analysis in a large population-based study of older adults in Germany with the aim of exploring the hitherto unknown association between active smoking-related DNA methylation and individual genetic variations. In particular, we aimed to identify the mQTLs within ± 50 kb from each of 151 previously reported active smoking-related CpG sites in whole blood samples [25] and to assess their relationships with active smoking exposure and allcause mortality.

Participant characteristics
Characteristics of the study population in the discovery and the validation panel with respect to smoking behaviors and lifestyle factors are summarized in Table 1. The average age of the participants of both subsets at the baseline was about 61 years. About half of the participants in each subset were ever smokers (current or former smokers), and around 18% still smoked at the time of recruitment. Female participants included a larger proportion of never smokers than males (discovery set, 67.9 vs. 28.6%; validation set, 63.3 vs. 21.4%). Average cumulative smoking exposure (pack-years) of current smokers was considerably higher than that of former smokers in both subsets (discovery set, 34.6 vs. 22.0; validation set, 33.1 vs. 19.4). Average time after smoking cessation (years) of former smokers in both subsets was also similar, approximately 17 years. The majority of participants in both subsets of the study Data missing for 5 and 7 participants, respectively, in discovery and validation panels population were overweight or obese, reported no or only low physical activity, and no or low amounts of alcohol drinking. During a median follow-up time of about 12 years (discovery set, 12.6 years; validation set, 12.2 years), 94 participants died in the discovery set (CVD = 30, cancer = 46, other diseases = 18) and 49 died in the validation set (CVD = 17, cancer = 21, other diseases = 11).

Identification of mQTLs for smoking-related CpG sites
For the 1396 SNP-CpG pairs consisting of 150 smokingrelated CpG sites and 909 corresponding SNPs (Fig. 1), 380 pairs were significant at a FDR < 0.05 in the discovery panel even after controlling for covariates (Additional file 1: Table S2; Additional file 2: Figure S1). These 380 pairs were then replicated in the validation panel by applying the fully adjusted mixed linear regression model. A subset of 246 pairs formed of 70 CpG sites and 192 SNPs reached the statistical significance level after FDR correction (FDR < 0.05; Table 2, Additional file 1: Table S3, Additional file 1: Table S4; Fig. 2). Eventually, 192 SNPs were designated as the mQTLs of 70 CpG sites. The pair cg23576855/rs75509302 showed the strongest inter-relationship (FDR-corrected p value = 8.86 e − 103). Among the 70 CpG sites, five were highly smoking-related loci (reported ≥ 6 times; Table 3), 14 were moderately smoking-related (4 or 5 times), and 51 were weakly smoking-related (2 or 3 times). These CpG sites with mQTLs were mainly located in the gene body (37/70), ten were located in transcription start sites (TSS1500) and 23 were in untranslated regions (UTR) or intergenetic regions (Additional file 1: Table S3). The largest number of mQTLs (n = 8) was found for locus cg06126421 within 6p21.33 (Table 3; Fig. 3). The coefficients of mQTLs ranged from −0.54 to 0.15.
The 192 mQTLs were mainly mapped on chromosomes 1 (16%), 6 (10%), and 7 (25%). Three SNPs were the most frequently identified mQTLs, rs75509392 (MAF = 0.144) for eight CpG sites within AHRR, rs79050605 (MAF = 0.202), and rs34835481 (MAF = 0.210) for five and six CpG sites located in GFI1, respectively (Table 4). We assessed the effects of the mQTLs on the DNA methylation changes as the absolute values of the coefficient changes of smoking status between the models without and with adjusting for corresponding mQTLs (carrier/non-carrier) among the 246 SNP-CpG pairs (Additional file 1: Table S5). Part of the mQTLs had opposite effects on different CpG sites (Additional file 1: Table S5, Additional file 3: Figure S2). For example, the minor allele of the SNP rs75509302 attenuated the association of smoking exposure with the methylation of cg11902777 by 5.2% (Additional file 1: Table S5, Additional file 3: Figure  S2). In contrast, this variant strengthened the demethylation of cg17287155 in response to different smoking behaviors by 2.44%.

Genetic contributions of mQTLs to the DNA methylation changes
As shown in Additional file 1: Table S5, the associations between smoking exposure and DNA methylation were changed by between 0.01 and 18.96% by the mQTLs and were categorized by the distances between genetic variants and CpG sites and the reported frequencies of CpG sites. We observed that the closest SNPs (distance < 10 kb) had a slightly lower impact on DNA methylation  levels than the mQTLs located ≥ 10 kb (Fig. 4a). Compared with the highly smoking-related CpG sites, the mQTLs affect the methylation levels of weakly smoking-related CpG sites the most, and the changes of moderately smoking-related CpG sites stayed in the intermediate position ( Fig. 4b, F = 4.91, p value = 0.008). Additionally, potential gene-environment interactions of the 192 mQTLs with smoking behaviors (current/never smoking) were assessed. Only for rs75509302, a significant interaction with smoking was observed regarding the methylation levels of cg23576855 (Tables 5, Additional file 1: Table S6).

Associations of mQTLs with active smoking exposure and all-cause mortality
Finally, we tested the associations of the 192 mQTLs (carrier/non-carrier) with different measurements of smoking exposure [ever smoking (current and former smoking) vs. never smoking, current smoking vs. never smoking, current smoking vs. former smoking, cumulative smoking (pack-years), durations of smoking (years), and the age of smoking initiation]. None of them was significantly associated with any smoking indicators after the correction of FDR (Additional file 1: Table S7). Similarly, the 192 SNPs were not significantly associated with all-cause mortality (death from CVD, cancer, and other chronic diseases) based on the results of the COX model (Additional file 1: Table S8).

Discussion
We conducted the first association study of 150 active smoking-related CpG sites and their corresponding SNPs located in the ± 50 kb region utilizing the genomic and epigenomic data of 949 participants from the ESTHER study. We found the DNA methylation levels of 70 CpG sites to be influenced by 192 proximal SNPs. These 192 mQTLs modified the DNA methylation changes in response to active smoking exposure, especially for the weakly/moderately active smoking-related CpG sites, but we did not observe any direct associations with active smoking exposure or all-cause mortality.   The mQTLs are presented in the vicinity of active smoking-related CpG sites. Locus cg06126421 (6p21.33) is one of the pronounced smoking-related CpG sites with eight mQTLs. Four of the SNPs could impair the hypomethylation of cg06126421, while others could accelerate this process. All of these mQTLs were located in a genomic region relating to inflammation and/or immune-related (malignant) diseases, such as allergies [26], multiple myeloma [27], diffuse large B cell lymphoma [28], and lung cancer [29]. We also identified a mQTL rs2227357 which could slightly promote the demethylation of cg03636183 (F2RL3), but no mQTLs were discovered near other well-established smokingrelated loci, like cg05575921 (AHRR) or cg19859270 (GPR15). However, eight CpG sites within the AHRR region were found to be modified by mQTLs. Among them, cg23576855 manifested the strongest connection with rs75509302 and was also the only hit which could be affected by the interaction between smoking status and corresponding mQTL. This phenomenon might be highly contributed by the genetic feature of this locus, which is also a CG → CA SNP annotated as rs6869832 [5,30], and shares a very high LD with rs75509302 in this study (R 2 = 0.94). Nevertheless, to our knowledge, the biological functions for most of the identified mQTLs in our study are not fully understood yet and need to be explored in further research.
The mQTLs may help to distinguish between the genetic and environmental effects on epigenome disparities. Researchers usually observed several outliers out of the predictive range of epigenetic signatures in EWASs [6,7,9,31]. One of the most plausible explanations is measurement bias that may result from recall bias or intentional underreporting [32]. Our finding provides another possibility that the deviations of DNA methylation levels might be caused by neighboring genetic variants. For highly smoking-related CpG sites, active smoking is still the main driver of DNA methylation changes. For instance, the SNP rs2227357 only contributed to about 0.01% of the changes of the methylation level of cg03636183 (F2RL3), and the SNP rs790051 altered only 0.37% of the methylation level of cg05951221 (2q37.1). However, the mQTLs affected the methylation levels of less robustly smoking-related loci much more. For instance, the SNPs rs78131 and rs2741302 explained nearly 19% of the changes of cg26963277 (KCNQ1) and a b Fig. 3 Locations of cg06126421 and eight mQTLs (carrier/non-carrier) (a) and distributions of methylation levels based on smoking status (b) in validation panel. Red line, FDR-corrected p value = 0.05; red dot, mQTLs; yellow triangle, cg06126421; blue bar, non-carriers of minor allele; red bar, carriers of minor allele cg27241845, respectively. While additional external validation studies certainly are needed, we speculate that this strong diversity of results might be a result of undetermined biological interactions between the SNP and CpG sites.
Smoking-related CpG sites have been recognized as informative signatures of smoking exposure and smokingrelated health outcomes [3,4]. Part of these 70 CpG sites with mQTLs have been reported to be highly associated with long-term smoking exposure [cg03636183 (F2RL3) and cg06126421 (6p21.33)] [9, 10], aging-related health outcomes, such as telomere length (cg21393163) [33] and the development of frailty [cg14753356, cg19589396, cg23667432 (ALPP) and cg25189904 (GNG12)] [34], and were even employed to construct a comprehensive index to predict smoking impact in buccal cells [35]. Together with previous studies [8,23], we suggest that future investigations utilizing smoking-related CpG sites might need to take the genotypes, especially the mQTLs of less robustly smoking-related loci, into consideration to account for their potential impact on DNA methylation levels.
Beyond the SNP-CpG associations, null associations of 192 mQTLs with active smoking and all-cause mortality additionally imply that these novel genetic variants might be independently associated with the DNA methylation changes and might not be involved in the pathophysiological development of smoking-related health outcomes. Therefore, these mQTLs might have the potential to be used in the causal inference tests between the CpG sites and smoking-related health outcomes as instrumental variables (Mendelian Randomization, MR) [36]. Recently, researchers have suggested a two-stage MR test to establish the causal role of epigenetic processes in pathways of diseases [37]. Larger population-based investigations with longitudinal design and repeated measurements of smoking exposure and epigenome data are warranted to evaluate these potential instrumental variables and obtain further insights into the plausibility of suggested causal effects of DNA methylation in the development of smoking-related diseases.
Major strengths of the present study include comprehensive information on a broad range of covariates in a population-based cohort and validation in an independent subgroup. Some limitations still have to be acknowledged in the interpretation of study results. First, smokingrelated shifts in leukocyte distribution might affect the The FDR-corrected p values of SNPs in fully adjusted mixed linear regression models, which controlled for age (years), sex, smoking status, random batch effect of methylation measurement, leukocyte distribution (Houseman algorithm), alcohol consumption (abstainer/low/intermediate/high), body mass index (BMI, underweight or normal weight/overweight/obese), physical activity (inactive/low/medium or high), prevalence of cardiovascular diseases (yes/no), prevalence of diabetes (yes/no), and prevalence of cancer (yes/no) associations of DNA methylation in whole blood samples with active smoking [38]. Hence, we adjusted for leukocyte distribution by the Houseman algorithm to restrict potential confounding from differential blood counts to the greatest possible extent [39]. Further studies are also needed to evaluate to what extent our results can be generalized to middle-aged individuals or non-Caucasians, as the ESTHER study was conducted in the older (aged 50-75 years), almost exclusively Caucasian population in southern Germany during a routine screening program. In addition, our study had limited power for detecting direct associations of mQTLs with smoking exposure and allcause mortality due to limited numbers of cases. Finally, we only measured a relatively small window of genetic regions (± 50 kb) in whole blood DNA due to the consideration of controlling for the pleiotropic effect or reverse causations from unknown genetic or epigenetic factors and the limited coverage of OncoArray [40], more mQTLs (cis-or trans-) for smoking-related CpG sites need to be established by expanding the search window. Since DNA methylation is highly tissue-specific, larger cohorts with various human tissues are also needed for more comprehensive evaluation of the whole landscape of genetic impact on the epigenome.

Conclusions
In conclusion, this study identified 192 mQTLs for 70 smoking-related CpG sites. These variants might theoretically reflect inherited differences in epigenetic states of people and their susceptibilities to smoking-related health outcomes. Incorporation of mQTLs might enhance the epigenetic-based assessments of smoking or smokingrelated health outcomes by accounting for potential confounding from genetic background. Our results need to be further validated and confirmed in additional studies with larger number of participants and more detailed assessment of genomic and epigenomic data, including the CpG sites that have not previously been replicated. Along with previous investigations on the epigenetic changes related to other environmental exposures or lifestyle factors, our study adds evidence for the complex interplays among genetic traits, epigenetic signatures, and environmental factors.

Study design and population
Study subjects were selected from the ESTHER study, an ongoing statewide population-based cohort study conducted in Saarland, a state located in southwestern Germany. Details of the study design have been reported previously [41]. Briefly, 9949 older adults (aged 50-75 years) were enrolled by their general practitioners during a routine health check-up between July 2000 and December 2002 and followed up thereafter. The current cross-sectional analysis is based on data and biospecimen collected at baseline. Two independent subgroups  The never smoking * genotype groups and current smoking * TT group were used as references; c The group of interaction between current smoking and listed genotype; d Never smoking was used as reference were selected as the discovery and the validation panel for DNA methylation analyses as previously described [33]. Briefly, the discovery panel included 581 participants recruited consecutively at the start of the ESTHER study between July and October 2000. The validation panel included 368 participants randomly selected from the participants recruited between October 2000 and March 2001. The study was approved by the ethics committees of the University of Heidelberg and the state medical board of Saarland, Germany. Written informed consent was obtained from all participants.

Data collection
Information on socio-demographic characteristics, lifestyle factors, and health status at baseline was obtained by standardized self-administered questionnaires. In particular, detailed information on lifetime smoking history was obtained, including current smoking status and intensity, age at initiation, and smoking intensities at various ages, as well as the age of quitting smoking for former smokers [42]. Additional information on body mass index (BMI) was extracted from a standardized form filled by the general practitioners during the health check-ups. Prevalent cardiovascular disease (CVD) at baseline was defined by either physician-reported coronary heart disease or a self-reported history of a major cardiovascular event, such as myocardial infarction, stroke, pulmonary embolism, or revascularization of coronary arteries. Prevalent diabetes was defined by physician diagnosis or the use of glucose-lowering drugs. Prevalent cancer [ICD-10 C00-C99 except non-melanoma skin cancer (C44)] was determined by self-report or record linkage with data from the Saarland Cancer Registry (http://www.krebsregister.saarland.de/ziele/ziel1.html; in German). Deaths during the follow-up (between 2000 and end of 2014) were identified by record linkage with population registries in Saarland. Participants migrated out of Saarland were censored at the date last known to be alive. Information about the major cause of death was obtained from death certificates provided by the local public health offices and was coded with ICD-10 codes.

DNA methylation data
Blood samples were taken during the health check-up and stored at −80°C until further processing. DNA from whole blood samples was collected using a salting out procedure [43]. DNA methylation profiles were extracted by the Illumina Human Methylation 450K BeadChip (Illumina, San Diego, CA, USA). As previously described [44], samples were analyzed following the manufacturer's instruction at the Genomics and Proteomics Core Facility of the German Cancer Research Center, Heidelberg, Germany. Illumina's GenomeStudio® (version 2011.1; Illumina, Inc.) was employed to extract DNA methylation signals from the scanned arrays (Module version 1.9.0; Illumina, Inc.). The methylation level of a specific CpG site was quantified as a β value ranging from 0 (no methylation) to 1 (full methylation). According to the manufacturer's protocol, no background correction was done and data were normalized to internal controls provided by the manufacturer. All controls were checked for inconsistencies in each measured plate. Probes with a detection p value > 0.05 were excluded from analysis. We utilized the Illumina normalization and preprocessing method implemented in Illumina's GenomeStudio®. We selected the profiles of 151 smoking-related loci which had been identified ≥ 2 times in previous smoking EWASs for the present analysis [25].

Genotyping data
Extracted DNA from blood cells was genotyped using the Illumina Infinium OncoArray BeadChip (Illumina, San Diego, CA, USA). General genotyping quality control assessment was as previously described [45]. Genotypes for common variants across the genome were imputed using data from 1000 Genomes Project (phase 3, Oct. 2014) with IMPUTE2 v2.3.2 after pre-phasing with SHAPEIT software v2.12. We set thresholds for imputation quality to retain both potential common and rare variants for validation. Specifically, poorly imputed SNPs defined by an information metric I < 0.70 were excluded. All genomic locations are given in NCBI Build 37/UCSC hg19 coordinates. All SNPs having a MAF < 1% were excluded. After imputation, the SNP set consisted of 9,198,808 genotyped and imputed SNPs. PLINK v1.90 was then used to extract SNPs for the required regions of interest [46]. As shown in Fig. 1, we first identified SNPs within 50 kb upstream and downstream from each of the 151 smoking-related CpG sites (470,582 SNP-CpG pairs), a window in which most SNPs with significant cis associations with CpG sites are located [13]. The locus cg01208318 was excluded without any corresponding SNPs in this restricted region. For each of the remaining 150 CpG sites, we excluded any SNPs with ≥ 1% missing values (n = 44,931), deviating from the Hardy-Weinberg equilibrium (HWE exact test's p value < 0.0001, n = 1), with a minor allele frequency ≤ 0.1 (n = 414,671) or with high linkage disequilibrium (LD, R 2 ≥ 0.5) (Additional file 1: Table S1). After the final quality control, 1396 SNP-CpG pairs with strongest SNPs remained for analysis, which were constituted of 150 CpG sites and 909 corresponding SNPs (Additional file 1: Table S2).

Statistical analyses
First, major socio-demographic characteristics, lifestyle factors, smoking behavior, and prevalence of major chronic diseases in both the discovery and the validation panel were summarized by descriptive statistics.
We then evaluated the associations between the methylation intensities of the 150 CpG sites and corresponding SNPs to identify mQTLs as follows. For all SNP-CpG pairs, we used a mixed linear regression model with methylation batch as a random effect in which the methylation level of CpG site was the outcome and each regional SNP was the predictor (categorical variable, coded into 0, 1, and 2 based on the numbers of the minor allele). The model was fully adjusted for the following covariates that have been shown to be associated with DNA methylation changes [47][48][49][50][51][52][53][54] [39], the prevalence of CVD (yes/ no), diabetes (yes/no), and cancer (yes/no) at the baseline. After correction for multiple testing by false discovery rate (FDR, Benjamini-Hochberg method [55]), SNP-CpG pairs with a FDR < 0.05 were selected and then analyzed in the validation panel. SNPs of the pairs with a FDR < 0.05 in the validation panel were eventually identified as the mQTL for the corresponding CpG site.
Furthermore, we tested the contributions of the identified mQTLs to the DNA methylation levels of corresponding CpG sites. Due to the limited number of individuals in the subgroup of minor homozygotes, we recoded the SNPs in order to use the dominant model, in which the heterozygote and minor homozygotes were combined as the carrier of minor allele and the major homozygotes were considered as non-carrier of the minor allele. We compared the coefficients of active smoking exposure (current vs. never smoking) in the fully adjusted model without the mQTLs (β 1 ) with the fully adjusted model including the mQTLs (β 2 ). The changes of coefficients were calculated as 100% * (β 1 -β 2 )/β 1 , and their absolute values were determined as the percentage change contributed by mQTLs. The percentage changes were categorized by the absolute distances (bp) between CpG sites and corresponding mQTLs and the reported frequencies of CpG sites. To explore whether the gene-environment interactions could modify the DNA methylation changes of smokingrelated CpG sites, we also tested whether interactions between the identified mQTLs and active smoking exposure (current vs. never smoking) could affect the impact of smoking on the methylation levels of corresponding CpG sites. The mQTLs, smoking status, and their interaction (mQTLs*smoking status) were added in the model as predictors, and the methylation levels of CpG sites were outcomes. After controlling for all the potential covariates, the interactions with a FDR < 0.05 were considered as methylation-related interactions for corresponding CpG sites.
Finally, we examined whether the identified mQTLs (carrier/non-carrier) were associated with six active smoking indicators: ever smoking (current and former smoking) vs. never smoking, current smoking vs. never smoking, current smoking vs. former smoking, cumulative smoking (packyears), durations of smoking (years), and the age of smoking initiation. The mixed linear models were fully adjusted for age (years), sex, smoking status, alcohol consumption, BMI, physical activity, the prevalence of CVD, diabetes, and cancer as described above. The mQTLs with a FDR < 0.05 were identified as smoking-related SNPs. We also assessed the associations of the significant mQTLs (carrier/non-carrier) with all-cause mortality in ESTHER study. Due to the limited number of deaths, we combined both subsets and performed the analysis using a multiple COX regression model. The model was adjusted for the above potential covariates and SNPs with a FDR < 0.05 were considered as all-cause mortality related variants.
Data cleaning and all aforementioned statistical analyses were performed by SAS version 9.4 (SAS Institute Inc., Cary, NC, USA). Manhattan plots for both panels were plotted by R package "ggplot2."

Additional files
Additional file 1: Table S1. Basic information of genomic data for 151 smoking-related CpG sites. Table S2 List of 1396 significant SNP-CpG pairs in discovery panel. Table S3 List of 246 significant SNP-CpG pairs in validation panel. Table S4 Results of regression analyses for 246 SNP-CpG pairs in validation panel. Table S5 Estimate changes of smoking status (never/ current smoking) between models with and without adjusting for mQTLs. Table S6 Impacts of mQTL-smoking interaction on the methylation levels of corresponding CpG sites in 949 participants (never/current smoking). Table S7 Associations of 192 mQTLs with active smoking indicators. Table S8