Skip to main content
Download PDF

Novel DNA methylation signatures of tobacco smoking with trans-ethnic effects

Clinical Epigenetics volume 13, Article number: 36 (2021) Cite this article

Abstract

Background

Smoking remains one of the leading preventable causes of death. Smoking leaves a strong signature on the blood methylome as shown in multiple studies using the Infinium HumanMethylation450 BeadChip. Here, we explore novel blood methylation smoking signals on the Illumina MethylationEPIC BeadChip (EPIC) array, which also targets novel CpG-sites in enhancers.

Method

A smoking-methylation meta-analysis was carried out using EPIC DNA methylation profiles in 1407 blood samples from four UK population-based cohorts, including the MRC National Survey for Health and Development (NSHD) or 1946 British birth cohort, the National Child Development Study (NCDS) or 1958 birth cohort, the 1970 British Cohort Study (BCS70), and the TwinsUK cohort (TwinsUK). The overall discovery sample included 269 current, 497 former, and 643 never smokers. Replication was pursued in 3425 trans-ethnic samples, including 2325 American Indian individuals participating in the Strong Heart Study (SHS) in 1989–1991 and 1100 African-American participants in the Genetic Epidemiology Network of Arteriopathy Study (GENOA).

Results

Altogether 952 CpG-sites in 500 genes were differentially methylated between smokers and never smokers after Bonferroni correction. There were 526 novel smoking-associated CpG-sites only profiled by the EPIC array, of which 486 (92%) replicated in a meta-analysis of the American Indian and African-American samples. Novel CpG sites mapped both to genes containing previously identified smoking-methylation signals and to 80 novel genes not previously linked to smoking, with the strongest novel signal in SLAMF7. Comparison of former versus never smokers identified that 37 of these sites were persistently differentially methylated after cessation, where 16 represented novel signals only profiled by the EPIC array. We observed a depletion of smoking-associated signals in CpG islands and an enrichment in enhancer regions, consistent with previous results.

Conclusion

This study identified novel smoking-associated signals as possible biomarkers of exposure to smoking and may help improve our understanding of smoking-related disease risk.

Background

Tobacco smoking is one of the leading preventable causes of death [46] and a leading risk factor for disease burden [16]. Smoking damages the airways and induces lung disease, such as chronic obstructive pulmonary disease (COPD), lung cancer and increased risk of multiple long-term conditions including heart disease and stroke [43]. Further wide-ranging health risks associated with smoking relate to several other cancers (including oral, gastro-intestinal, urinary tract and reproductive cancers), bone health, gum disease, macular degeneration, type 2 diabetes, rheumatoid arthritis and altered immune function [15, 16].

Smoking exposure has been robustly associated with changes in DNA methylation. Many smoking-associated CpG-sites, or differentially methylated positions (DMPs), have been identified and replicated to date [18, 25]. The majority of recent studies linking blood DNA methylation levels with smoking exposure have used the Infinium HumanMethylation450 BeadChip (450k array) to generate genome-wide DNA methylation profiles. In the largest smoking-methylation study of adults to date [25], Illumina 450k methylation profiles for 15,907 blood samples from 16 cohorts were used to identify 2623 smoking-DMPs in or near to 1405 genes. Many of these DMPs were already linked to smoking by previous studies [7, 49]. The most consistently associated CpG-sites across independent studies are located in or near the AHRR, RARA, PRSS23, F2RL3, GPR15 and GNG12 genes and in chromosomal regions 2q37.1 and 6p21.33. These epigenetic signals have also been used to develop DNA methylation-based biomarkers of smoking status such as EpiSmokEr, currently the most robust smoking classifier [5].

Recently human genome-wide DNA methylation studies have started to explore the additional genome coverage afforded by the new Infinium MethylationEPIC BeadChip (EPIC array), which assays 850,000 CpG-sites at nearly double the coverage of the Illumina 450k array. Recent findings based on DNA methylation profiles in lung tissue [37] and saliva [2] identified smoking-associated methylation signals profiled only by the EPIC array, which suggests that novel smoking EPIC methylation signals may also be identified in blood. Here, we assessed the association between smoking and DNA methylation levels profiled using the Illumina EPIC array in 1407 whole blood samples from four UK population cohorts and pursued replication in independent cohorts. Our analyses compared DNA methylation levels between current and never smokers in a cross-sectional analysis to identify novel smoking-associated differentially methylated positions (smoking-DMPs). The observed signals were subsequently compared between former and never smokers to assess whether these alterations persist after smoking cessation. The results have potential to improve existing biomarkers of smoking that could be used to infer smoking exposure where it is not known. Furthermore, the findings may identify novel genes which may explain specific disease risk mechanisms in smokers.

Methods

Discovery phase participants

The participants in the discovery phase included 1407 individuals from four UK population cohorts—TwinsUK and three birth cohorts, including the 1970 British Cohort Study (BCS70), the National Child Development Study (NCDS) or 1958 British birth cohort, and the MRC National Survey of Health and Development (NSHD) or 1946 British birth cohort (Additional file 1: Note). The sample included 235 individuals from BCS70, 529 individuals from NCDS, 236 individuals from the NSHD and 407 individuals from TwinsUK (Table 1). The majority of participants included for DNA methylation profiling were not selected for a specific phenotype distribution or environmental exposure, although sample selection strategies included minimizing exposure or phenotype data missingness across samples (see Additional file 1: Note). An exception was a subset of 294 participants in the NCDS cohort selected for extremes of child and adulthood adversity [6, 41] (see Additional file 1: Note). Smoking information was obtained from questionnaire data collected at the time of DNA methylation profiling. All study participants provided informed consent, and ethical approval was granted by local research ethics committees (see Additional file 1: Note). The overall discovery sample included 269 current, 495 former and 643 never smokers from the UK population (see Additional file 1: Note).

Table 1 Sample characteristics
Full size table

DNA methylation profiles

DNA samples extracted from whole blood were profiled using the Illumina MethylationEPIC BeadChip (Illumina EPIC) array, which targets over 850,000 CpG-sites including more than 90% of the probes on the 450k and additional CpG-sites predominantly in enhancers [36]. DNA methylation levels were determined using methylation beta-values, defined as the ratio of the methylated bead signal to the sum of the unmethylated bead signal plus the methylated bead signal plus 100 [10]. Methylation beta-values range between 0 at unmethylated CpG-sites and 1 at fully methylated CpG-sites. DNA methylation quality control measures included multiple checks and normalizations. Altogether, 72,471 probes were identified as cross-reactive or polymorphic and excluded from the analysis. Cross-reactive probes were defined as those mapping to multiple locations of the in silico bisulfite converted human genome allowing for two mismatches. Probes were also excluded if they targeted polymorphic CpG sites with minor allele frequency (MAF) > 5% in the UK10K haplotype reference panel. Probes with greater than 5% missingness were also excluded from the analysis, spanning 1348 probes. DNA methylation levels were then normalized using ENmix [48]. The number of CpG-sites included in downstream analyses was 710,658.

Peripheral blood cell proportions

A significant difference in blood cell-type proportions has previously been reported between smokers and never smokers [32]. This finding, along with the observation that DNA methylation profiles also vary by blood cell type, highlights the need to take into account blood cell proportion differences in the analyses. Blood cell-type proportions were estimated for monocytes, granulocytes, immune cells (Natural Killer (NK) cells, CD8 and CD4) and plasmablasts using the approach proposed by Houseman et al. [22]. Correlations between estimated blood cell types were investigated to inform a set of covariates for inclusion in epigenetic linear models of association. As a result, blood cell subtype covariates used in downstream analyses included monocytes, granulocytes, NK cells and CD8-naive cells.

Epigenetic association analyses

Epigenome-wide association scans (EWAS) were carried out within each cohort dataset, followed by meta-analyses. DNA methylation values at each CpG-site were normalized to N(0, 1) prior to fitting linear models. EWAS focused on the comparison of DNA methylation profiles between smokers and never smokers at 710,658 CpG-sites genome-wide. Linear models of association, run separately for each cohort dataset, compared DNA methylation values as the response variable to smoking status as a predictor. For the three birth cohorts, linear models were fitted (lm function in R) where normalized methylation level at each CpG-site was the response variable and predictors included smoking status, sex, blood cell-type proportion, BMI and methylation chip and position of the sample on the chip. EWAS were run separately in the two NCDS cohort subsets, due to different sample selection strategies and proportion of smokers. In the TwinsUK sample linear mixed-effects models were fitted using lme4 [3] and lmerTest [27] in R, using the same fixed-effects covariates as for the birth cohort models excluding sex because all twins were female, but including age as well as random effects variables for family and zygosity.

Following individual cohort sample EWAS, a fixed-effects inverse variance weighted meta-analysis was applied to combine results across cohorts. Meta-analysis, performed using GWAMA [34], was carried out across the 5 datasets of 1407 subjects in total. To minimize effects attributed to heterogeneity across cohorts we only considered meta-analysis results that did not exhibit strong evidence for heterogeneity (Qp < 10% and an I2 > 50%, [21]). Multiple-testing adjustment of the meta-analysis results was performed using both the Benjamini and Hochberg false discovery rate (FDR) threshold (FDR = 1%) and a Bonferroni-corrected threshold (P value = 6.25 × 10−8 ≈ 0.05/800,000), with both sets of results reported. Methylation effect sizes were calculated using the same linear models, but without normalizing DNA methylation levels to N(0, 1) prior to data analysis. Evidence for genomic inflation was assessed using λ, which is the ratio of the median of the empirically observed distribution of the test statistic to the expected median, thus quantifying the extent of the bulk inflation and the excess false-positive rate. Results from the largest smoking-methylation study to date, Joehanes et al. [25] (at FDR < 0.05), were used to assess whether smoking-DMPs are novel.

Once smoking differentially methylated positions (smoking-DMPs) were identified from the meta-analysis, follow-up analysis included a comparison of DNA methylation levels between former and never smokers at the 952 smoking-DMPs. In these follow-up analyses a Bonferroni-corrected threshold was applied for multiple testing (P value = 5.25 × 10−5 ≈ 0.05/952).

Replication in trans-ethnic samples

Many of the previously established smoke exposure effects on DNA methylation have been robustly replicated across different ethnicities. With this in mind and to assess whether the signals that we detected are conserved across different ethnic populations, we pursued replication in two samples of non-European ancestry. Replication of the novel smoking-DMPs was carried out firstly in American Indian participants from the Strong Heart Study [28] and secondly in African-American participants from the Genetic Epidemiology Network of Arteriopathy study [8] (see Additional file 1: Note).

The first replication sample consisted of 2325 American Indians aged 45–74 (893 smokers and 684 never smokers). DNA methylation was measured using the EPIC array in whole blood samples collected in 1989–1991. Pre-processing was conducted according to Illumina’s recommendations, and snoob and Regression on Correlated Probes (RCP) normalizations were applied [13, 35]. Batch effects by sample plate, sample row and DNA isolation time were corrected using comBat (sva R package). Peripheral blood cell proportions were estimated as in the other cohorts. Replication was pursued at 525 of novel smoking-DMPs profiled by the EPIC array only that also passed quality control assessment in the replication dataset. The same covariates and modeling approach as for the birth year cohorts was used in the replication sample.

The second replication sample consisted of 1100 African-Americans with a mean age of 56 (179 smokers and 666 never smokers). A total of 1106 samples at GENOA Phase I and 304 samples at GENOA Phase II were assessed using the Illumina HumanMethylationEPIC BeadChip. Raw IDAT files were imported using Minfi R package [1]. The shinyMethyl R package [13] was used to visualize the raw intensity data and identify sex mismatches and outliers were removed. Individual probes with detection P value < 10e−16 were considered to be detected successfully [29], and samples and probes with detection rate < 10% were removed. Samples with incomplete bisulfite conversion identified using the QCinfo() function in the ENmix R package were removed [48]. Sample identity was checked using the 59 SNP probes implemented in the EPIC chip and mismatched samples removed. Next, Noob was used for individual background and dye-bias normalization [14]. Since two types of probes are present on the EPIC BeadChip (Infinium I and Infinium II), we used the RCP method to adjust for probe-type bias [35]. After exclusions, a total of 857,121 probes in 1100 samples at Phase I and 294 samples at Phase II were available for analysis. Peripheral blood cell proportions were estimated as in the other cohorts. Replication was pursued at 526 of novel smoking-DMPs profiled by the EPIC array only. Linear mixed-effects models were fitted using the same methodology as for the TwinsUK cohort (lme4 and lmerTest), and the same fixed-effects covariates as for the birth cohort models but including age as well as a random effect variable for sibship. Meta-analysis across the two replication samples was carried out using the same methodology as for the main analysis. A Bonferroni-corrected threshold was applied (P value = 9.5 × 10−5 ≈ 0.05/525), and only results showing the same direction of association as in the discovery sample were considered.

SLAMF7 gene expression analysis

Gene expression analysis for SLAMF7 was carried out in a sample from the TwinsUK cohort with available whole blood gene expression data. RNA-seq data generation and pre-processing have been previously described in detail [17]. In summary, STAR software v2.4.0.1 [12] was used to align reads to the hg19 reference genome. Samples with fewer than 10 million aligned reads were excluded. Following this process, there were 383 whole blood samples remaining including 227 never smokers, 30 current smokers and 126 former smokers, where 162 individuals overlapped with the TwinsUK sample used in the main analysis. Gene counts were transformed into trimmed mean of M-values (TMM)-adjusted counts per million (CPMs) and inverse-normalized prior to all downstream analyses.

A mixed effect linear model was fitted (lme4) with smoking as the predictor and gene expression as the response. Covariates included fixed effects: insert-size median, mean GC content and random effects: primer index, date of sequencing, zygosity, family and RNA extraction batch. A P value was determined using lmerTest with a significance threshold of 0.05.

Genomic annotation and pathway analysis

Genomic annotation of smoking-DMPs was carried out initially using the EPIC Illumina manifest, for the purpose of identifying novel CpG-sites that are specific to the EPIC-array, for allocating CpG-sites to genes, and relative to CpG-density, including CpG island (CGI), CGI-shore, CGI-shelf, and open sea. Further genomic annotations took into account data from the ENCODE project [20]. We explored if smoking-DMPs mapped within ChromHMM [11] categories to assess enrichment or depletion in smoking-DMPs relative to different functional genomic domains, including insulators, enhancers and specific transcription factors binding sites. In the enrichment analysis, we considered all smoking-DMP probes mapping to a specific annotation category compared to the total number of probes tested that mapped to that category. For each genome annotation category, the results show the log fold change for smoking-DMP probes compared to total probes tested, and the significance of the difference is based on a Fisher’s exact test. Pathway analysis was carried out for genes annotated to smoking-DMPs using Ingenuity Pathway Analysis (IPA; QIAGEN Inc. https://www.qiagenbioinformatics.com/products/ingenuitypathway-analysis).

Prediction of smoker status

We tested several models for classifying smoking status based on subsets of the newly identified smoking-DMPs. The sensitivity and specificity of each classification was assessed using receiver operative curve (ROC), implemented using the pROC package in R [38]. The analyses were carried out in the combined dataset of 1113 subjects who were not selected for phenotypic extremes, excluding the NCDS subset of 294 individuals. Training datasets were created taking 60% of the combined dataset at random. A test dataset was created with the remaining 40%. Altogether, 20 random samples were taken creating 20 random training and test set combinations.

For each of the training datasets, a generalized linear model was fitted based on the predictors, including un-adjusted DNA methylation levels at the candidate CpG-site(s) and covariates (age, sex, blood cell-type proportions, BMI, methylation chip and position of the sample on the chip) using the R glm function. The test dataset was then loaded into the derived model with outcomes predicted using the R predict function. The average AUC was determined for each methylation value combination.

Methylation combinations for distinguishing between current smokers and never smokers included cg05575921 (AHRR) and cg00045592 (SLAMF7) on their own, and then in combination.

Methylation combinations for distinguishing any smoke exposure (that is, either current or former smoker) and never smokers were explored in three models. The first model included the 5 ex-smoking DMPs with the largest effect size in the 450k array (cg21566642, cg05575921, cg01940273, cg25189904, cg12803068), the second model included the 5 ex-smoking DMPs with largest effect size sites from the EPIC array (cg14391737, cg21566642, cg05575921, cg25189904, cg05533761), and the third model included the 5 ex-smoking DMPs with the largest effect size overall (cg14391737, cg21566642, cg05575921, cg25189904, cg05533761).

Results

DNA methylation profiles and smoking were explored in 1407 total whole blood samples from 4 UK population cohorts in the discovery stage (269 current smokers, 495 ex-smokers and 643 never smokers), and in 3425 whole blood samples from American Indians and African-Americans in the replication stage (Table 1). The primary analyses focused on identification of smoking differentially methylated positions (smoking-DMPs) between current smokers and never smokers, and follow-ups explored smoking-DMPs genomic distribution, pathway analysis, stability upon smoking cessation, and predictive value.

Epigenome-wide association analysis: current versus never smokers

Meta-analysis comparing DNA methylation profiles in current smokers (N = 269) and never smokers (N = 643) identified 952 CpG sites or smoking-DMPs in 500 genes that were statistically differentially methylated at a Bonferroni-adjusted threshold (P value = 6.25 × 10−8). At a more relaxed threshold (FDR 1%), there were 3348 CpG sites in 1632 genes (Additional file 2: Table S1). There was evidence for genomic inflation (see Methods) with an overall λ of 1.28 (Fig. 1a), which is consistent with other meta-analyses of smoke exposure [25]. Smoking-DMPs are spread across the genome, consistent with previous observations (Fig. 1b).

Fig. 1
figure 1

Methylation association results in current versus never smokers. a Manhattan plot of genome-wide results for methylation association with smoking. Smoking-DMPs are indicated above the Bonferroni-adjusted threshold (red line). b Quantile–quantile (QQ) plot for CpG-site association in current versus never smoker

Full size image

Of the 952 smoking-associated CpG-sites, 422 have previously been identified [25] and are represented on both the 450 k array and the EPIC array (Additional file 2: Table S2). The strongest association was observed at cg21566642 in the 2q37.1 region (P value = 1.6 × 10−121), which has been observed previously in multiple studies. Furthermore, cg05575921 in AHRR, which is the most frequently observed association in previous studies, was the most significantly associated CpG site in our study that was also annotated to a gene (P = 2.9 × 10−99). The effect size at this site, based on unadjusted DNA methylation beta values, was also the largest where current smokers exhibited on average 25% lower DNA methylation values compared to never smokers, which is also consistent with previous studies [25, 49].

Of the 952 significant CpG-sites, 526 smoking-DMPs represent novel CpG sites only profiled by the EPIC array (Additional file 2: Table S2). The 526 smoking-DMPs were annotated to 277 genes, of which 80 genes have not previously been linked to differential methylation with smoking. The most significant novel CpG-site, cg00045592, was annotated to gene SLAMF7, (P value = 1.39 × 10−59), which has not been previously linked to methylation changes with smoking (Fig. 2). The effect size, or the unadjusted mean methylation in non-smokers compared to smokers, was 10% change in mean methylation, representing one of the larger effect sizes (13th largest). Furthermore, a second DMP at this gene, site cg04009575, was also identified as a smoking-DMP (P value = 3.02 × 10−10). In addition to novel DMPs mapping to the 80 novel genes, many other novel DMPs mapped to genes linked to smoking by CpG sites measured by the Illumina 450k array. Examples include cg1431737 in PRSS23 (P value = 6.5 × 10–92) and cg17739917 in RARA (P value = 1.97 × 10–92). Altogether, 68% of novel DMPs are upstream of or within genes, with 20% in new genes and 48% in genes where methylation changes have previously been associated with smoking.

Fig. 2
figure 2

Novel smoking-associated DNA methylation signals in SLAMF7. a coMET plot [33] describing the genomic region of epigenome-wide association between smoking and SLAMF7 methylation (top panel) showing the two smoking-DMPs cg00045592 and cg04009575, along with functional annotation of the region (middle panel) where broad ChromHMM regions are displayed using UCSC genome browser color schemes, and pattern of co-methylation at the 12 CpG sites in the EPIC array annotated to SLAMF7 (bottom panel). b Boxplot showing a comparison of DNA methylation levels at cg00045592 between smokers, former smokers and never smokers in the combined TwinsUK, NCDS, NSHD, BCS70 data set. c Boxplot showing a comparison of DNA methylation levels at cg04009575 between smokers, former smokers and never smokers in the Strong Heart Study

Full size image

The majority of smoking-DMP effects were hypomethylated in current smokers (74% of Bonferroni-adjusted sites, and 90% of the 100 sites with the lowest associated P values). The average effect sizes, measured as the difference in mean unadjusted methylation levels between never smokers and current smokers, were broadly similar between the hypomethylated (average effect size of 3.3% for hypomethylated sites) and hypermethylated (average effect size of 2.2% for hypermethylated) signals. However, the largest effect sizes detected overall (up to 25% mean difference) were observed at hypomethylated sites.

Replication of novel smoking signals

Novel smoking-DMPs that were profiled only by the EPIC array were evaluated for replication. Many smoking signals have previously been replicated across ancestries, and therefore we tested whether our results were robust across different ancestries. Replication was pursued in two independent samples, the first included 2325 American Indian participants from the Strong Heart Study [28], and the second included 1100 African-American participants in the Genetic Epidemiology Network of Arteriopathy (GENOA) study. A recent study in the SHS sample by Domingo-Relloso et al. [9] analyzed EPIC array profiles of whole blood samples for smoke exposure in connection with cadmium levels in urine and also identified novel smoking-DMPs. In addition to the distinct ancestry, this cohort also has a higher proportion of smokers, and furthermore one of the leading causes of death for American Indians is CVD, for which smoking is a risk factor. Of the 526 novel sites, 389 (74%) replicated in the SHS sample alone at a Bonferroni-corrected threshold (P value = 9.5 × 10–5) with the same direction of association. This included cg00045592, annotated to gene SLAMF7, which replicated with a P value of 1.45 × 10–23. At nominal significance (P value = 0.05), there were 500 (95%) novel sites with the same direction of association effect. The GENOA sample also has distinct ancestry, and African-Americans have the highest incidence of hypertension, for which smoking is a risk factor. Recent analysis of the dataset validated smoking-DMPs identified in saliva [2]. Of the 526 novel sites, 418 (79%) replicated at a Bonferroni-corrected threshold (P value = 9.5 × 10–5) with the same direction of association in the GENOA sample alone. This also included cg00045592, annotated to gene SLAMF7, which replicated with a P value of 7.57 × 10–31. At nominal significance (P value = 0.05), there were 508 (97%) novel sites with the same direction of association effect.

In a meta-analysis of the two replication cohorts in altogether 3425 individuals, 486 (92%) of the 526 novel smoking-DMPs replicated at a Bonferroni-corrected threshold (P value = 9.5 × 10–5) with the same direction of association. This also included cg00045592, annotated to gene SLAMF7, which replicated with a P value of 3.54 × 10–51. At nominal significance (P value = 0.05), there were 522 (99%) novel sites with the same direction of association effect (Additional file 2: Table S4).

Genome annotation analysis

We assessed the genome distribution for all 952 smoking-DMPs identified in our study, relative to all probes assayed by the EPIC array. We explored enrichment and depletion of smoking-DMPs across different genomic categories. The enrichment analysis considered 18 genomic annotation categories (Fig. 3a). The strongest effect was a clear enrichment of smoking-DMPs in enhancer regions as predicted by ChromHMM (36% of smoking-DMPs relative to 12% of probes tested), consistent with previous studies [25] and with the hypothesis that smoking exposure impacts regulatory genomic features. Furthermore, also consistent with previous work, there was an enrichment of smoking-DMPs in gene bodies (50% of smoking-DMPs relative to 42% of probes tested). Finally, there was a depletion of smoking-DMPs in CpG islands (7% relative to 18% of tested probes), which is consistent with previous observations that CpG islands are less dynamic in response to exposure [50].

Fig. 3
figure 3

Smoking-DMPs annotation and pathway analysis. a Log fold changes relating to the proportion of probes annotated to particular genomic locations and the proportion of smoking-DMPs. Only annotation categories with statistically significant results based on Fisher’s exact test are shown b. IPA canonical pathway analysis, 20 lowest P values shown

Full size image

Pathway analysis

The 952 smoking-DMPs were annotated to 500 genes. We explored evidence for enrichment of these genes within different biological processes, focusing on canonical pathways using IPA (see Methods). Altogether, the 500 genes identified enrichment for 101 molecular pathways (P value = 0.01, Additional file 2: Table S5). Among the 20 most-enriched pathways, 8 (40%) relate to cancer, with the remainder relating to cell signaling and growth pathways, neuronal health, cardiovascular health and insulin receptor activity (Fig. 3b).

SLAMF7 gene expression follow-up analysis

To assess potential functional impacts of the novel smoking-DMPs identified in SLAMF7, we explored SLAMF7 gene expression levels in smokers and never smokers. The analysis was carried out in 383 blood samples from the TwinsUK cohort with available blood RNAseq levels [17]. A nominally significant difference in SLAMF7 gene expression (P = 0.02) was detected where current smokers had reduced levels of expression of SLAMF7.

Epigenome-wide association analysis: former versus never smokers

To assess how smoking-DMPs behave after smoking cessation, we carried out a cross-sectional analysis comparing DNA methylation levels in 497 former and 643 never smokers at the 952 smoking DMPs. If DNA methylation levels at smoking DMPs persist after smoking cessation, we would expect significant differences in DNA methylation at all 952 signals, as we did in the comparison of current and never smokers. However, we observed that at the majority of smoking-DMPs there was no significant difference between former and never smokers, suggesting reversal of smoking-associated DNA methylation levels upon smoking cessation. Altogether, there were 37 differentially methylated sites at a Bonferroni-adjusted threshold (P value = 1.5 × 10–5), at which our results were consistent with persistent effects after smoking cessation (41 signals at FDR 1%, Additional file 2: Table S3). Of the 37 Bonferroni signals, 16 represented novel signals only profiled on the EPIC array, 20 replicated previously identified persistent smoking-methylation signals observed by Joehanes et al. [25] present on both arrays, and one signal profiled on both arrays was novel. The 16 novel EPIC-specific signals annotated to 11 genes, and the most significant signal in a novel gene was in SLAMF7 (cg00045592, P value = 8.8 × 10–7), where former smokers have on average around 2.4% lower unadjusted methylation levels than never smokers. The majority (92%) of the 37 signals are hypomethylated in former smokers, and the effect sizes were smaller than observed between the smokers and never smokers. On average, the mean effect size across all 37 signals was a hypomethylation of 2% in former smokers, and the largest effect size obtained was 6% (cg14391737).

Methylation-based detection of smoking status

We assessed the performance of the peak smoking-DMPs as classifiers of smoking status. The DNA methylation level at our peak novel smoking-DMP, cg00045592 (SLAMF7), is a moderately good classifier for distinguishing between current smokers and never smokers with an AUC of 0.87. The overall peak smoking-DMP and widely replicated cg05575921 (AHRR) is a much stronger classifier between smokers and never smokers with an AUC of 0.95. Adding SLAMF7 to AHRR marginally improves the AUC to 0.96. We next explored the predictive value of different DMPs to distinguish current or former exposure to smoking (‘ever smoked’), by comparing current and former smokers to never smokers. Neither cg00045592 (SLAMF7) or cg05575921 (AHRR) are very good predictors of current or former smoke exposure compared to never smokers, with AUC values of 0.68 and 0.76, respectively. Using the 5 450 k-only array signals with the largest effect sizes from the comparison of former and never smokers (cg21566642, cg05575921, cg01940273, cg25189904, cg12803068) results in a slightly better classification of ever smoked exposure with an average AUC of 0.79 (Additional file 3: Figure S1a). The classifier is marginally improved by including novel EPIC-specific sites, where the 5 overall signals with the largest effect sizes from the comparison of former and never smokers (cg14391737, cg21566642, cg05575921, cg25189904, cg05533761) result in an average AUC of 0.80 (Additional file 3: Figure S1b).

Discussion

The current study explored smoking status and DNA methylation levels profiled on the EPIC array in 1407 whole blood samples from individuals across 4 independent UK population-based cohorts. We identified 952 CpG smoking-DMPs between smokers and never smokers at Bonferroni adjustment. Of these, 422 replicated previous findings and 526 signals were novel, involving CpG sites profiled only by the EPIC array. Of the 526 novel signals, 486 replicated in the American Indian and African-American replication sample. The different ethnicity of the replication samples, along with the high proportion of smokers and their propensity for CVD in the American Indian population and prevalence of hypertension in the African-American population, indicates that our smoking-DMPs signals are robust and trans-ethnic. While many of the novel signals were in genes previously found to be associated with smoke exposure, 80 novel genes were also identified where the top novel signal was in the SLAMF7 gene. The genes that harbored smoking-related methylation signals were enriched to fall in biological pathways related to cancer, cell signaling and growth pathways, neuronal health, cardiovascular health and insulin receptor activity. Although the majority of signals showed at least some evidence for reversal in DNA methylation levels with smoking cessation, we also identified smoking-DMPs which persisted after smoking cessation. Overall, this is consistent with previous studies and observations that upon smoking cessation risk of smoking-related diseases reverts to non-smoker levels over time [45], but for some outcomes risk is not fully reversible.

The most strongly associated novel smoking-DMP (cg00045592) is in a gene not previously linked to smoking, SLAMF7, which also contains a secondary smoking-DMP signal (cg04009575, ranked 489th of the 952 smoking-DMPs). Both sites are located in an enhancer and therefore likely to influence the expression of genes including SLAMF7. A follow-up analysis showed that smokers had significantly lower levels of expression of SLAMF7. Both signals replicated in both replication cohorts. DNA methylation at cg00045592 alone can distinguish current smokers and never smokers with an average AUC of 0.87. SLAMF7 is a protein coding gene and a member of the Signalling Lymphocyte Activation Molecule Family (SLAMF), a family of receptors with a role in both innate and adaptive immunity. SLAMF7 is expressed on immune cells, recognizes and binds to itself, in turn leading to the downstream activation of natural killer cells [19]. Smokers have been shown to have suppressed NK activation [23], which in turn has been shown to lead to a reduction in NK tumor surveillance in smokers, and consequently an increased propensity for tumors [30]. SLAMF7 has been linked to several diseases, for example, it is associated with systemic lupus erythematosus [26], multiple myeloma [24] and immunoregulatory interactions [31]. Methylation of SLAMF7 has also been identified as a regulator in atherosclerosis [47]. SLAMF7 has also been suggested as a possible therapeutic target for rheumatoid arthritis, a disease more common in smokers [44]. Our analysis showed that the differential methylation at cg0045592 persisted after smoking cessation, retaining on average at 2.4% lower DNA methylation levels compared to never smokers, which given its potential functional role may have implications for altered immune function in current and former smokers.

Genome annotation analysis showed a clear enrichment in smoking-DMPs in enhancer regions. Although this result has previously been observed on the 450k array in blood [25], a key strength of the EPIC array is the more detailed genome coverage of enhancers. This enrichment provides evidence that smoke exposure leads to methylation changes in regulatory regions, which in turn likely result in functional changes. This is consistent with previous studies which have shown gene expression changes as a result of smoke exposure [25, 42]. In addition, the rest of the annotation analysis showed results consistent with previous work [25, 50].

The results from the current study can be used to not only understand possible implications of smoke exposure on gene function, but also as potentially useful biomarkers of smoke exposure. Smoking is a confounder in many epidemiological studies given its wide-ranging impact on human health [43]. Not only is a molecular biomarker of smoking useful to confirm self-reported smoke exposure data, but it can be also valuable as a predictor of smoking in samples with missing smoking data, errors in self-reported status or in forensic settings. Previous work has developed predictors of current smoking based on 450k array smoking-DMPs, observing good predictive value (AUC of 0.95) of AHRR in particular as a smoking predictor. Our novel signals do not significantly improve smoking classification compared to using AHRR alone. However, a subset of smoking effects persists after cessation. Although we are unable to assess exactly how long these effects last for, they can be used as measures of current or previous smoke exposure. Therefore, our smoking-DMPs may be useful toward developing a classifier of present or past smoke exposure, that is, distinguishing current and former smokers from never smokers. In an attempt to address this, we observed an AUC of 0.80 for distinguishing present or past smoke exposure using the top 5 signals from both arrays. Including novel sites present only on the EPIC array affected prediction only marginally (0.79 vs 0.8). The results could also be used for further studies toward prediction of future smoking-related health outcomes.

Our findings may provide insights into molecular changes underlying tobacco smoke exposure risk effects in disease. Pathway analysis of the smoking-DMPs indicated that the genes annotated to smoking-DMPs are strongly linked to cancer, with 40% of the top 20 significant pathways relating to cancer. These included a number of different cancer types, not solely those relating to lung cancer, which concurs with epidemiological studies which show that smoking is a risk factor for several types of cancer. The methylome varies by tissue type, but we see a relationship between smoking-DMPs found in blood and those found in tissues directly relevant to some of these cancers, for example, in lung tissue. Specifically, there is an overlap between the smoking-DMPs identified here and those found by Ringh et al. [37] in lung tissue, where 32 of the 952 sites overlap particularly at sites annotated to AHRR and in genes previously found to be associated with cancer (KCNMA1, CDH23, LRP5) [[4, 39, 40]]. In addition to implications for cancer, the pathway analysis results also revealed strong links between smoking-DMPs and neuronal health. The novel gene findings are a starting point for analyzing whether and how DNA methylation alterations could lead to smoking-related disease.

There are some limitations to the current study. While the birth cohorts are not in themselves selective and should therefore lead to a broad sample representative of the population, not all individuals will continue to participate in cohort sweeps. In this study the majority of subjects were selected based on minimizing missing data across a range of variables, with low level of oversampling of specific subgroups. This approach could lead to some bias in sampling more engaged cohort participants. We observed genomic inflation in the smoking epigenome-wide analyses, which has also been observed in previous smoking-methylation meta-analysis reports [25]. Although this observation raises potential concerns about false positives, our results are in line with previous studies, replicating many previously reported signals, as well as 92% of novel signals in the trans-ethnic samples. We took a conservative approach by removing meta-analysis results that may be due to heterogeneity. This approach is likely to have removed some results that may well be true positives in some sample subsets. Whole blood samples are a mixture of different cell types, and we used methylation estimates of cell proportions to address this in the analysis. However, these are likely to be over generalizations and ideally analyses in specific cell subpopulations should be carried out. Another limitation is that DNA methylation was only studied in blood and not in other tissues (e.g., lung, adipose), which have previously been explored in the context of smoking on the 450k array. Other tissues such as lung tissue are likely to produce more insight into disease mechanisms, although with the exception of saliva, they lend themselves less well to biomarker detection due to the relative difficulty of sampling. Follow-up analyses, which our study did not explore, include a full gene-expression analysis to assess functional impacts at the novel smoking-DMPs, longitudinal analyses to characterize the stability of smoking-DMPs upon cessation over long timescales, and studies of these signals in samples of other ethnicities.

In conclusion, our study identified hundreds of novel smoking-methylation signals, including those annotated to genes not previously associated with smoke exposure. Some of the novel signals persist after cessation of smoking. The findings have potential to act as biomarkers of exposure to smoking and may improve our understanding of smoking-related disease risk.

Availability of data and materials

The majority of DNA methylation datasets in the current study, except those from the Strong Heart Study replication cohort, are available in the public domain. TwinsUK methylation data are uploaded on the ReShare UK Data Service, under Data collection id 853,526. NSHD methylation data access is through https://doi.org/10.5522/NSHD/S202. NCDS and BCS70 methylation data access is through https://doi.org/10.5255/UKDA-SN-5594-2. Access to further individual-level data including phenotype data can be applied for through each cohort data access committee. For information on access and how to apply, see https://twinsuk.ac.uk/resources-for-researchers/access-our-data/ (TwinsUK), http://www.nshd.mrc.ac.uk/data (NSHD), https://beta.ukdataservice.ac.uk/datacatalogue/series/series?id=2000032 (NCDS) and https://beta.ukdataservice.ac.uk/datacatalogue/series/series?id=200001 (BCS70). The Strong Heart Study DNA methylation data can be available to external investigators by following the procedures established by the Strong Heart Study Steering Committee in agreement with the study’s Tribal partners. These procedures are available at https://strongheartstudy.org/. Phenotype data for GENOA participants are available from the Database of Genotypes and Phenotypes (dbGaP): phs001401.v2.p1. Methylation data are from the Gene Expression Omnibus (GEO): GSE157131. Due to IRB restriction, mapping of the sample IDs between genotype data (dbGaP) and methylation data (GEO) cannot be provided publicly, but are available upon written request to JAS and SLRK.

References

  1. Aryee M, Jaffe A, Corrada-Bravo H, Ladd-Acosta C, Feinberg A, Hansen K, Irizarry R. Minfi: a flexible and comprehensive bioconductor package for the analysis of infinium DNA methylation microarrays. Bioinformatics. 2014;30:1363–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  2. Barcelona V, Huang Y, Brown K, Liu J, Zhao W, Yu M, Kardia S, Smith J, Taylor J, Sun Y. Novel DNA methylation sites associated with cigarette smoking among African Americans. Epigenetics. 2019;14(4):383–91.

    PubMed  PubMed Central  Article  Google Scholar 

  3. Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67(1):1–48.

    Article  Google Scholar 

  4. Bloch M, Ousingsawat J, Simon R, et al. KCNMA1 gene amplification promotes tumor cell proliferation in human prostate cancer. Oncogene. 2007;26:2525–34.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  5. Bollepalli S, Korhonen T, Kaprio J, Anders S, Ollikainen M. EpiSmokEr: a robust classifier to determine smoking status from DNA methylation data. Epigenomics. 2019;11(13):1469–86.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  6. Borghol M, Suderman M, McArdle W, Racine A, Hallett M, Pembrey M, Hertzman C, Power C, Szyf M. Associations with early-life socio-economic position in adult DNA methylation. Int J Epidemiol. 2012;41(1):62–74.

    PubMed  Article  PubMed Central  Google Scholar 

  7. Breitling L, Yang R, Korn B, Burwinkel B, Brenner H. Tobacco-smoking-related differential DNA methylation: 27K discovery and replication. Am J Hum Genet. 2011;88(4):450–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  8. Daniels P, Kardia S, Hanis C, Brown C, Hutchinson R, Boerwinkle E, Turner S, Network GE, of Arteriopathy study, . Familial aggregation of hypertension treatment and control in the genetic epidemiology network of arteriopathy (GENOA) study. Am J Med. 2004;116:676–81.

    PubMed  Article  Google Scholar 

  9. Domingo-Relloso A, Riffo-Campos A, Haack K, Rentero-Garrodo P, Ladd-Acosta C, Fallin D, Tang W, Herreros-Martinez M, Gonzalez J, Bozack A, Cole S, Navas-Acien A, Tellez-Plaza M. Cadmium, smoking, and human blood DNA methylation profiles in adults from the strong heart study. Environ Health Perspect. 2020;128(6):67005.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  10. Du P, Zhang X, Huang C-C, Jafari N, Kibbe W, Hou L, Lin S. Comparison of Beta-value and M-value methods for quantifying methylation levels by microarray analysis. BMC Bioinform. 2010;11:587.

    CAS  Article  Google Scholar 

  11. Ernst J, Kellis M. ChromHMM: automating chromatin-state discovery and characterization. Nat Methods. 2012;9:215–6.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  12. Fairfax B, Humburg P, Makino S, Naranbhai V, Wong D, Lau E, Jostins L, Plant K, Andrews R, McGee C, Knight J. Innate immune activity conditions the effect of regulatory variants upon monocyte gene expression. Science. 2014;343:1246949–1246949.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  13. Fortin J, Fertig E, Hansen K. shinyMethyl: interactive quality control of illumina 450k DNA methylation arrays in R. F1000Res. 2014;3:175.

    PubMed  PubMed Central  Article  Google Scholar 

  14. Fortin J, Triche T, Hansen K. Preprocessing, normalization and integration of the Illumina HumanMethylationEPIC array with minfi. Bioinformatics (Oxford, England). 2017;33(4):558–60.

    CAS  Google Scholar 

  15. GBD 2015 Tobacco Collaborators. Smoking prevalence and attributable disease burden in 195 countries and territories, 1990–2015: a systematic analysis from the Global Burden of Disease Study 2015. Lancet. 2017;389(10082):1885–906.

    Article  Google Scholar 

  16. GBD 2017 Disease and Injury Incidence and Prevalence Collaborators. Global, regional, and national incidence, prevalence, and years lived with disability for 354 diseases and injuries for 195 countries and territories, 1990–2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet Glob Health Metr. 2018;392(10159):1789–858.

    Google Scholar 

  17. Glastonbury A, et al. Cell-type heterogeneity in adipose tissue is associated with complex traits and reveals disease-relevant cell-specific eQTLs. Am J Hum Genet. 2019;104(6):1013–24.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  18. Guida F, Sandanger T, Castane R, Campanella G, Polidoro S, Palli D, Krogh V, Tumino R, Scerdote C, Panico S, Severi G, Kyrtopoulos S, Georgiadis P, Vermeulen R, Lund E, Vineis P, Chadeau-Hyam M. Dynamics of smoking-induced genome-wide methylation changes with time since smoking cessation. Hum Mol Genet. 2015;24(8):2349–59.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  19. Guo H, Cruz-Munoz M-E, Wu N, Robbins M, Veillette A. Immune cell inhibition by SLAMF7 is mediated by a mechanism requiring Src Kinases, CD45, and SHIP-1 that is defective in multiple myeloma cells. Mol Cell Biol. 2014;35(1):41–51.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  20. Harrow J, Frankish A, Gonzalez J, Tapanari E, Diekhans M, Kokocinski F, Aken B, Barrell Zadissa A, Searle S, Barnes I, Bignell A, Boychenko V, Hunt T, Kay M, Mukherjee G, Rajan J, Despacio-Reyes G, Saunders G, Steward C, Harte R, Lin M, Howald C, Tanzer A, Derrien A, Chrast J, Walters N, Balasubramanian S, Pei B, Tress M, Rodriguez J, Ezkurdia L, van Baren J, Brent M, Haussler D, Kellis M, Valencia A, Reymond A, Gerstein M, Guigo R, Hubbard T. GENCODE: the reference human genome annotation for The ENCODE Project. Cold Spring Harb Lab Press. 2012;22:1760–74.

    CAS  Google Scholar 

  21. Higgins J, Thompson S, Deeks J, Altman D. Measuring inconsistency in meta-analyses. BMJ. 2003;327(7414):557–60.

    PubMed  PubMed Central  Article  Google Scholar 

  22. Houseman E, Accomando W, Keostler D, Christensen B, Marsit C, Nelson H, Wiencke J, Kelsey K. DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinform. 2012;13:86.

    Article  Google Scholar 

  23. Inoue C, Takeshita T, Kondo H, et al. Cigarette smoking is associated with the reduction of lymphokine-activated killer cell and natural killer cell activities. Environ Health Prev Med. 1996;1:14.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  24. Ishibashi M, Soeda S, Sasaki M, Handa H, Imai Y, Tanaka N, Tanosaki S, Ito S, Odajima T, Sugimori H, Asayama T, Sunakawa M, Kaito Y, Kinoshita R, Kuribayashi Y, Onodera A, Moriya K, Tanaka J, Tsukune Y, Komatsu N. Clinical impact of serum soluble SLAMF7 in multiple myeloma. Oncotarget. 2018;9(78):34784–93.

    PubMed  PubMed Central  Article  Google Scholar 

  25. Joehanes R, Just AC, Marioni RE, Pilling LC, Reynolds LM, Mandaviya PR, Guan W, Xu T, Elks CE, Aslibekyan S, Moreno-Macias H, Smith JA, Brody JA, Dhingra R, Yousefi P, Pankow JS, Kunze S, Shah SH, McRae AF, Lohman K, Sha J, Absher DM, Ferrucci L, Zhao W, Demerath EW, Bressler J, Grove ML, Huan T, Liu C, Mendelson MM, Yao C, Kiel DP, Peters A, Wang-Sattler R, Visscher PM, Wray NR, Starr JM, Ding J, Rodriguez CJ, Wareham NJ, Irvin MR, Zhi D, Barrdahl M, Vineis P, Ambatipudi S, Uitterlinden AG, Hofman A, Schwartz J, Colicino E, Hou L, Vokonas PS, Hernandez DG, Singleton AB, Bandinelli S, Turner ST, Ware EB, Smith AK, Klengel T, Binder EB, Psaty BM, Taylor KD, Gharib SA, Swenson BR, Liang L, DeMeo DL, O’Connor GT, Herceg Z, Ressler KJ, Conneely KN, Sotoodehnia N, Kardia SL, Melzer D, Baccarelli AA, van Meurs JB, Romieu I, Arnett DK, Ong KK, Liu Y, Waldenberger M, Deary IJ, Fornage M, Levy D, London SJ. Epigenetic signatures of cigarette smoking. Cardiovasc Genet. 2016;9(5):436–47.

    CAS  Article  Google Scholar 

  26. Karampetsou M, Comte D, Kis-Toth K, Kyttaris V, Tsokos G. Expression patterns of signaling lymphocytic activation molecule family members in peripheral blood mononuclear cell subsets in patients with systemic lupus erythematosus. PLoS ONE. 2017;12(10):e01860.

    Article  CAS  Google Scholar 

  27. Kuznetsova A, Brockhoff PB, Christensen RHB. lmerTest package: tests in linear mixed effects models. J Stat Softw. 2017;82(13):1–26.

    Article  Google Scholar 

  28. Lee E, Welty T, Fabsitz R, Cowan L, Le N, Oopik A, Cucchiara A, Savage P, Howard B. The strong heart study: a study of cardiovascular disease in American Indians: design and methods. Am J Epidemiol. 1990;132(6):1141–55.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  29. Lehne B, Drong A, Loh M, Zhang W, Scott W, Tan S, Afzal U, Scott J, Jarvelin M, Elliott P, McCarthy M, Kooner J, Chambers JC. A coherent approach for analysis of the Illumina HumanMethylation450 BeadChip improves data quality and performance in epigenome-wide association studies. Genome Biol. 2015;16:37.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  30. Lu L, Zavitz C, Chen B, Sussan K. Cigarette smoke impairs NK cell-dependent tumour immune surveillance. J Immunol. 2007;178:936–43.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  31. Malaer JD, Mathew PA. CS1 (SLAMF7, CD319) is an effective immunotherapeutic target for multiple myeloma. Am J Cancer Res. 2017;7(8):1637–41.

    CAS  PubMed  PubMed Central  Google Scholar 

  32. Malenica M, Pmjavorac B, Bego T, Dujic T, Semiz S, Skrbo S, Gusic A, Hadzic A, Causevic A. Effect of cigarette smoking on haematological parameters in healthy population. Med Arch. 2017;71(2):132–6.

    PubMed  PubMed Central  Article  Google Scholar 

  33. Martin TC, Yet I, Tsai PC, Bell JT. coMET: visualisation of regional epigenome-wide association scan results and DNA co-methylation patterns. BMC Bioinform. 2015;16:131.

    CAS  Article  Google Scholar 

  34. Mägi R, Morris AP. GWAMA: software for genome-wide association meta-analysis. BMC Bioinform. 2016;11:288.

    Article  Google Scholar 

  35. Niu L, Xu Z, Taylor J. RCP: a novel probe design bias correction method for Illumina Methylation BeadChip. Bioinformatics (Oxford, England). 2016;32(17):2659–63.

    CAS  Article  Google Scholar 

  36. Pidsley R, Zotenko E, Peters T, Lawrence M, Risbridger G, Molloy P, Van Djik S, Muhlhausler B, Stirzaker C, Clark S. Critical evaluation of the Illumina MethylationEPIC BeadChip microarray for whole-genome DNA methylation profiling. Genome Biol. 2016;17:208.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  37. Ringh M, Hagemann-Jensen M, Needhamsen M, Kular L, Breeze C, Sjöholm L, Slavec L, Kullberg S, Wahlström J, Grunewald J, Brynedal B, Liu Y, Almgren M, Jagodic M, Öckinger J, Ekström T. Tobacco smoking induces changes in true DNA methylation, hydroxymethylation and gene expression in bronchoalveolar lavage cells. EBioMedicine. 2019;46:290–304.

    PubMed  PubMed Central  Article  Google Scholar 

  38. Robin X, Turck N, Hainard A, Tiberti N, Lisacek F, Sanchez J-C, Müller M. pROC: an open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinform. 2011;12:77.

    Article  Google Scholar 

  39. Roslan Z, Muhamad M, Selvaratnam L, Ab-Rahim S. The roles of low-density lipoprotein receptor-related proteins 5, 6, and 8 in cancer: a review. J Oncol. 2019;2019:4536302.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  40. Sannigrahi MK, Srinivas CS, Deokate N, Rakshit S. The strong propensity of Cadherin-23 for aggregation inhibits cell migration. Mol Oncol. 2019;13(5):1092–109. https://doi.org/10.1002/1878-0261.12469.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  41. Suderman M, Borghol N, Pappas J et al. Childhood abuse is associated with methylation of multiple loci in adult DNA. BMC Med Genomics. 2014;7(13)

  42. Tsai P, Glastonbury C, Eliot M, Bollepalli S, Yet I, Castillo-Fernandez J, Carnero-Montero E, Hardiman H, Martin T, Vickers A, Mangino M, Ward K, Pietilainen K, Deloukas P, Spector TD, Vinuela A, Loucks E, Ollikainen M, Kelsey K, Small K, Bell JT. Smoking induces coordinated DNA methylation and gene expression changes in adipose tissue with consequences for metabolic health. Clin Epigenet. 2018;10:26.

    Article  CAS  Google Scholar 

  43. U.S. Department of Health and Human Services. The health consequences of smoking—50 years of progress: a report of the surgeon general. 2014. https://www.cdc.gov/tobacco/data_statistics/sgr/50th-anniversary/index.htm. Accessed May 2019.

  44. Woo J, Vierboom M, Kwon H, Chao D, Ye S, Li J, Lin K, Tang I, Belmar N, Hartman T, Breedveld E, Vexler V, Hart B, Law D, Starling G. PDL241, a novel humanized monoclonal antibody, reveals CD319 as a therapeutic target for rheumatoid arthritis. Arithritis Res Ther. 2013;15(6):R207.

    Article  CAS  Google Scholar 

  45. Woodward M, Lam TH, Barzi F. Smoking, quitting, and the risk of cardiovascular disease among women and men in the Asia-Pacific region. Int J Epidemiol. 2015;34:1036–45.

    Google Scholar 

  46. World Health Organization. WHO global report on trends in prevalence of tobacco smoking. 2015. https://apps.who.int/iris/handle/10665/156262. Accessed April 2019.

  47. Xia X, Gu M, Jia X, Wang X, Wu C, Guo J, Zhang L, Du Y, Wang J. Integrated DNA methylation and gene expression analysis identifies SLAMF7 as a key regulator of atherosclerosis. Aging (Albany NY). 2018;10(6):1324–37.

    CAS  Article  Google Scholar 

  48. Xu Z, Niu L, Li L, Taylor J. ENmix: a novel background correction method for Illumina HumanMethylation450 BeadChip. Nucleic Acids Res. 2016;44(3):e20.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  49. Zeilinger S, Kuhnel B, Klopp N, Baurecht H, Kleinschmidt A, Gieger C. Tobacco smoking leads to extensive genome-wide changes in DNA methylation. PLoS ONE. 2013;8(5):e63812.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  50. Ziller MJ, Gu H, Müller F, Donaghey J, Tsai LT, Kohlbacher O, De Jager PL, Rosen ED, Bennett DA, Bernstein BE, Gnirke A, Meissner A. Charting a dynamic DNA methylation landscape of the human genome. Nature. 2013;500:477–81.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

Download references

Acknowledgements

The authors would like to thank the research volunteers who participated in the study. Cohort-specific acknowledgments are provided in Additional file 1: Note.

Funding

This study was supported by the Economic and Social Research Council (ESRC, Grant Number ES/N000404/1 to J.T.B). Data governance for access to UK birth cohort data was provided by the METADAC data access committee, funded by ESRC, Wellcome and MRC (2015–2018: Grant Number MR/N01104X/1 2018–2020: Grant Number ES/S008349/1). This work made use of data and samples generated by the 1958 Birth Cohort (NCDS), which is managed by the Centre for Longitudinal Studies at the UCL Institute of Education, funded by the ESRC (Grant Number ES/M001660/1); blood sample collection was funded by MRC (Grant G0000934). Access to these resources was enabled via the Wellcome Trust and MRC: 58FORWARDS Grant (108439/Z/15/Z). Before 2015 biomedical resources were maintained under the Wellcome Trust and Medical Research Council 58READIE Project (Grant Numbers WT095219MA and G1001799). Work undertaken at Great Ormond Street Hospital/University College London, Institute of Child Health, is in part supported by the National Institute for Health Research Biomedical Research Centre at Great Ormond Street Hospital for Children NHS Foundation Trust and University College London. The UK Medical Research Council also provided core funding for the MRC National Survey of Health and Development (Grant MC_UU_00019/1). The TwinsUK study is supported by the Wellcome Trust; European Community’s Seventh Framework Programme (FP7/2007–2013); the National Institute for Health Research (NIHR) Biomedical Research Centre based at Guy’s and St Thomas’ NHS Foundation Trust and King’s College London and/or the NIHR Clinical Research Facility. The SHS study received support from the National Heart, Lung, and Blood Institute (NHLBI) of the USA (under contract numbers 75N92019D00027, 75N92019D00028, 75N92019D00029 and 75N92019D00030) and previous Grants (R01HL090863, R01HL109315, R01HL109301, R01HL109284, R01HL109282 and R01HL109319 and Cooperative Agreements: U01HL41642, U01HL41652, U01HL41654, U01HL65520 and U01HL65521), from the National Institutes of Health Sciences of the USA (R01ES021367, R01ES025216, P42ES010349, P30ES009089) and from “la Caixa” Foundation (ID 100010434) (fellowship code “LCF/BQ/DR19/11740016”). The GENOA study received funding from U01HL054457, RC1HL100185, R01HL119443, R01HL133221, R01HL141292.

Author information

Affiliations

  1. Department of Twin Research and Genetic Epidemiology, King’s College London, London, UK

    C. Christiansen, J. E. Castillo-Fernandez, J. S. El-Sayed Moustafa, P.-C. Tsai, K. S. Small & J. T. Bell

  2. Department of Environmental Health Sciences, Columbia University Mailman School of Public Health, New York, USA

    A. Domingo-Relloso & A. Navas-Acien

  3. Department of Chronic Diseases Epidemiology, National Center for Epidemiology, Carlos III Health Institute, Madrid, Spain

    A. Domingo-Relloso & M. Tellez-Plaza

  4. Department of Statistics and Operative Research, University of Valencia, Valencia, Spain

    A. Domingo-Relloso

  5. Department of Epidemiology, School of Public Health, University of Michigan, Ann Arbor, USA

    W. Zhao, S. L. R. Kardia & J. A. Smith

  6. MRC Unit for Lifelong Health and Ageing, Institute of Cardiovascular Science, University College London, London, UK

    J. Maddock, A. Wong, D. Kuh & R. Hardy

  7. Population Health Program, Texas Biomedical Research Institute, San Antonio, USA

    K. Haack & S. A. Cole

  8. School of Population Health and Environmental Sciences, King’s College London, London, UK

    M. Molokhia

  9. MRC Integrative Epidemiology Unit, University of Bristol, Bristol, UK

    M. Suderman & C. Relton

  10. Population, Policy and Practice Research and Teaching Department, UCL Great Ormond Street Institute of Child Health, London, UK

    C. Power & C. Relton

  11. Centre for Longitudinal Studies, UCL Social Research Institute, University College London, London, UK

    A. Goodman & G. B. Ploubidis

  12. Department of Biomedical Sciences, Chang Gung University, Taoyuan, Taiwan

    P.-C. Tsai

  13. Genomic Medicine Research Core Laboratory, Chang Gung Memorial Hospital, Linkou, Taiwan

    P.-C. Tsai

Authors
  1. C. Christiansen
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. J. E. Castillo-Fernandez
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. A. Domingo-Relloso
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. W. Zhao
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. J. S. El-Sayed Moustafa
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. P.-C. Tsai
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. J. Maddock
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. K. Haack
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. S. A. Cole
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. S. L. R. Kardia
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. M. Molokhia
    View author publications

    You can also search for this author in PubMed Google Scholar

  12. M. Suderman
    View author publications

    You can also search for this author in PubMed Google Scholar

  13. C. Power
    View author publications

    You can also search for this author in PubMed Google Scholar

  14. C. Relton
    View author publications

    You can also search for this author in PubMed Google Scholar

  15. A. Wong
    View author publications

    You can also search for this author in PubMed Google Scholar

  16. D. Kuh
    View author publications

    You can also search for this author in PubMed Google Scholar

  17. A. Goodman
    View author publications

    You can also search for this author in PubMed Google Scholar

  18. K. S. Small
    View author publications

    You can also search for this author in PubMed Google Scholar

  19. J. A. Smith
    View author publications

    You can also search for this author in PubMed Google Scholar

  20. M. Tellez-Plaza
    View author publications

    You can also search for this author in PubMed Google Scholar

  21. A. Navas-Acien
    View author publications

    You can also search for this author in PubMed Google Scholar

  22. G. B. Ploubidis
    View author publications

    You can also search for this author in PubMed Google Scholar

  23. R. Hardy
    View author publications

    You can also search for this author in PubMed Google Scholar

  24. J. T. Bell
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

CC and JTB designed the study and outlined the main conceptual ideas. JTB, KH, SAC, KS, MS, CP, CR, AW, DK, AG, SLRK, JAS, MT-P, A-N-A, GBP, RH generated the primary datasets or contributed toward data generation. CC led the data analysis. JEC-F, AD-R, P-CT, JM, MM, WZ and JSE-SM contributed to data analysis and results interpretation. CC and JTB wrote the article, and all authors provided critical feedback and helped shape the research, analysis and manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to J. T. Bell.

Ethics declarations

Ethics approval and consent to participate

Ethical approval for TwinsUK was granted by the National Research Ethics Service London-Westminster, the St Thomas’ Hospital Research Ethics Committee (EC04/015 and 07/H0802/84). Ethical approval for NSHD was granted by the Central Manchester Research Ethics Committee (07/H1008/168 and 07/H1008/245) and the Scotland A Research Ethics Committee (08/MRE00/12). Ethical approval for NCDS was granted by the London Central REC (14/LO/0097, 12/LO/2010 and 08/H0718/29) and by the South East MREC (01/1/44). Ethical approval for the BCS70 was granted by the London Central REC (14/LO/1371, 11/LO/1560) and the Southampton and South West Hampshire REC (08/H0504/144). All research participants have signed informed consent prior to taking part in any research activities. The Strong Heart Study is approved by Columbia University institutional research board (IRB), as well as other IRBs at multiple institutions, local tribal IRBs and Indian Health Service IRBs. Institutional review board approval was received for the GENOA study protocol through the University of Mississippi Medical Center and the University of Michigan IRBs.

Consent for publication

All authors have approved the manuscript for publication.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Note.

Additional file 2: Table S1.

Statistically significant CpGs in relation to current versus never smoking status at false discovery rate (FDR) < 0.01. Table S2. Statistically significant CpGs in relation to current versus never smoking status at a Bonferroni-adjusted threshold (P value = 6.25 × 10–8). Table S3. Statistically significant CpGs in relation to former versus never smoking status at false discovery rate (FDR) < 0.01. Table S4. Statistically significant CpGs in replication meta-analysis relation to current versus non-smoking status at a Bonferroni-adjusted threshold (P value = 9.5 × 10−5). Table S5. Pathway Analysis (IPA) results.

Additional file 3: Figure S1.

Receiver operating characteristic (ROC) curves predicting current smokers and smoke exposure. a Top 5 CpG from 450 k only predicting smoke exposure and b Top 5 CpG from EPIC predicting smoke exposure.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Christiansen, C., Castillo-Fernandez, J.E., Domingo-Relloso, A. et al. Novel DNA methylation signatures of tobacco smoking with trans-ethnic effects. Clin Epigenet 13, 36 (2021). https://doi.org/10.1186/s13148-021-01018-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13148-021-01018-4

Keywords

  • Smoking
  • Epigenetics
  • DNA methylation
  • Environment
  • Lifestyle
  • SLAMF7

Clinical Epigenetics

ISSN: 1868-7083

Contact us