The effect of smoking on DNA methylation of peripheral blood mononuclear cells from African American women
© Dogan et al.; licensee BioMed Central Ltd. 2014
Received: 24 July 2013
Accepted: 17 February 2014
Published: 22 February 2014
Regular smoking is associated with a wide variety of syndromes with prominent inflammatory components such as cancer, obesity and type 2 diabetes. Heavy regular smoking is also associated with changes in the DNA methylation of peripheral mononuclear cells. However, in younger smokers, inflammatory epigenetic findings are largely absent which suggests the inflammatory response(s) to smoking may be dose dependent. To help understand whether peripheral mononuclear cells have a role in mediating these responses in older smokers with higher cumulative smoke exposure, we examined genome-wide DNA methylation in a group of well characterized adult African American subjects informative for smoking, as well as serum C-reactive protein (CRP) and interleukin-6 receptor (IL6R) levels. In addition, complementary bioinformatic analyses were conducted to delineate possible pathways affected by long-term smoking.
Genome-wide DNA methylation analysis with respect to smoking status yielded 910 significant loci after Benjamini-Hochberg correction. In particular, two loci from the AHRR gene (cg05575921 and cg23576855) and one locus from the GPR15 gene (cg19859270) were identified as highly significantly differentially methylated between smokers and non-smokers. The bioinformatic analyses showed that long-term chronic smoking is associated with altered promoter DNA methylation of genes coding for proteins mapping to critical sub-networks moderating inflammation, immune function, and coagulation.
We conclude that chronic regular smoking is associated with changes in peripheral mononuclear cell methylation signature which perturb inflammatory and immune function pathways and may contribute to increased vulnerability for complex illnesses with inflammatory components.
Smoking is the largest preventable cause of morbidity and mortality in the United States. It largely exerts these effects by increasing liability to complex disorders, such as cancer, chronic obstructive pulmonary disease (COPD), type 2 diabetes (T2DM) and obesity . Smoking driven chronic diseases contribute to early death, disabilities, and strain the health care system . Therefore, understanding the mechanism(s) through which smoking increases vulnerability to these disorders may establish new avenues for prevention or treatment of these complex disorders. Although some of the details remain unclear, one of the key mechanisms through which smoking may increase liability to these complex disorders is inflammation.
Although serological quantification of well characterized serum markers such as C-reactive protein (CRP) and interleukin 6 receptor (IL6R) may provide a partial understanding of inflammatory changes with respect to smoking, this approach provides limited comprehension of molecular perturbation at a genome-wide scale.
A surge in recent publications has suggested that smoking associated changes in DNA methylation may contribute to these perturbations. This surge began with sporadically published single gene studies that linked smoking to changes in MAOA promoter methylation as well as to an increased risk for coronary heart disease mediated through methylation changes at F2RL3[3, 4]. But in the past two years, these more limited examinations have been joined by several genome-wide investigations that have identified a growing number of loci whose methylation status is associated to smoking. The initial attempt at a more systematic approach was reported by Breitling and colleagues who used the Illumina HumanMethylation 27K BeadChip to probe DNA from peripheral mononuclear cell pellets and identified several candidate loci including F2RL3 (cg03636183), GPR15 (cg19859270) and ORAI2 (cg02564523) . The first truly genome-wide results using the then newly introduced Illumina HumanMethylation 450K BeadChip were first reported by Monick and colleagues who studied methylation in lymphoblast and lung macrophage DNA and found a large number of loci with particular emphasis on differential methylation at the Aryl Hydrocarbon Receptor Repressor (AHRR) . Six months later, using cord blood, Joubert and colleagues confirmed and extended these findings at AHRR and further nominated GF1 and CY1A1 as genes affected by smoking status . Finally, in a study just published, Zeilinger and colleagues identified a larger set of findings that confirmed the prior loci noted above and extended the gene list to include loci such as HIVEP3 (cg15542713), and CACNA1D (cg15417641 and cg21188533) .
These genome-wide finding present a potential portal for a better integrated understanding of pathways through which smoking potentially accelerates disease states. Previous studies have established several smoking associated disease pathways. One such is the cyclooxygenase-2 (COX-2) pathway where the expression of COX-2 induced by smoking leads to an increase in prostaglandin E2 (PGE2) that mediates tumor progression and an increase in thromboxane A2 (TxA2), contributing to tumor growth . The COX-2 pathway that has been shown to be involved in numerous smoking related cancers is only one in many similar pathways that could conceivably hold the key to designing new therapeutic interventions and in some instances aid in improving current medical care.
A more recently developed method of identifying pathways perturbed by smoking is by coalescing information from network biology and epigenetics, specifically DNA methylation. DNA methylation is a critical mediator between the genome and the environment. Although genome-wide DNA methylation analyses allow the identification of differentially methylated CpG sites and genes, there is growing evidence that genes do not function independently, but rather in networks. Hence, it is crucial to not only delineate individual loci affected by long-term smoking, but, in addition, use network theory to translate that single locus methylation information into a more holistic understanding of the effects of smoking on the function of the proteome and the cell in general.
This global understanding is necessary to help understand the reasons why some individuals are affected more severely than others. Even though the rate of adults smoking is roughly equivalent in Whites and African Americans, African Americans disproportionately experience smoking associated medical comorbidities . The reason for this disparity is not known but both genetic and sociological factors have been postulated. Epigenomic interrogation of this cohort would allow the disparity examination from a different perspective.
In prior studies, we have analyzed the relationship of smoking to peripheral mononuclear cell DNA methylation in young adult African American smokers . To determine whether smoking and peripheral mononuclear cells may have a role in the systemic inflammation in long-term smokers, we conducted a study to examine the relationship of smoking to genome-wide DNA methylation in an independent and well characterized African American cohort who are older and have greater cumulative smoke exposure. In addition, we use a recently developed network based analysis technique to understand the effect of these smoking associated changes in DNA methylation on protein interaction networks. Our data suggest that regular chronic smoking is associated with highly significant differential methylation at a large number of loci including AHRR (cg05575921 and cg23576855) and GPR15 (cg19859270). Furthermore, the bioinformatic analyses demonstrate enrichment of loci involved in regulation of inflammation, immune function and coagulation.
Clinical characteristics of the cohort
Clinical and demographic characteristics of female FACHS subjects participating in the genome-wide methylation studies
48.1 ± 7
48.7 ± 11
31.1 ± 8
35.3 ± 7
Current average cigarettes/day
5 ± 8
≥ ½ and < 10
≥ 10 and < 20
5.02 ± 7.6
2.89 ± 3.2
227 ± 58
209 ± 56
Genome-wide DNA methylation analyses
The thirty most significantly associated probes to smoking status after genome-wide correction
Average beta values
Two other CpG probes, cg05575921 and cg23576855 (ranked 2nd and 4th, respectively) were situated in the body of the aryl hydrocarbon receptor repressor (AHRR) gene. The protein encoded by this gene provides feedback inhibition of AHR activation of the xenobiotic pathway by several different mechanisms [12, 13]. Our analyses indicate that on average, at these two loci, smokers tended to be hypomethylated. Moreover, our analyses are consistent with prior results of young adult African American smokers. In young adult African Americans, cg05575921 was also highly significantly affiliated to smoking and smokers tended to be hypomethylated. However, although hypomethylation was apparent at cg23576855 in both young and older smokers, this site was only significant with respect to long-term smoking and therefore could be a potential cumulative smoking biomarker.
Validation of Illumina array findings
The top 30 differentially regulated gene ontology pathways from the GoMiner analysis
Phospholipid-translocating ATPase activity
Phospholipid transporter activity
Lipoprotein transporter activity
Cellular component organization
ATPase activity coupled
Retinoic acid receptor activity
Negative regulation of transcription DNA-dependent
Lipid transporter activity
Negative regulation of RNA metabolic process
Cellular component organization or biogenesis
Hydrolase activity acting on acid anhydrides
Auditory receptor cell differentiation
ATPase activity coupled to transmembrane movement of substances
ATPase activity coupled to movement of substances
Hydrolase activity acting on acid anhydrides catalyzing transmembrane movement of substances
Single-stranded RNA binding
Cellular component maintenance
Negative regulation of osteoblast differentiation
Hydrolase activity acting on acid anhydrides in phosphorus-containing anhydrides
Multicellular organism reproduction
Multicellular organismal reproductive process
Aldehyde reductase activity
Protein-protein interactions: network analyses using the miPALM algorithm
Even though conventional genome-wide DNA methylation analyses yielded consistent identification of differentially methylated CpG loci at the AHRR gene, this analysis did not capture the connectivity of co-regulated genes as a result of long-term smoking. Hence, we employed an alternative approach for understanding the possible effects of smoking associated DNA methylation changes on cellular function by reducing the dimensionality of data and taking advantage of our understanding of the protein interactome. We used the miPALM algorithm which employs a greedy search strategy on weighted networks to re-analyze the genome-wide methylation data utilizing only the subset of DNA methylation data mapped to gene promoter regions.
Top 10 pathways from BiNGO pathway analysis of protein sub-network depicted in Figure 3 (a)
Immune effector process
Activation of plasma proteins in acute inflammatory response
Humoral immune response
Peptide maturation by peptide bond cleavage
Acute inflammatory response
Activation of immune response
Positive regulation of immune response
Top 10 pathways from BiNGO pathway analysis of protein sub-network depicted in Figure 3 (b)
Response to chemical stimulus
Response to amino acid stimulus
Response to acid
Response to amine stimulus
Response to organic nitrogen
All pathways from BiNGO pathway analysis of protein sub-network depicted in Figure 3 (c)
JAK-STAT cascade in growth hormone signaling pathway
Growth hormone receptor signaling pathway
Cellular response to growth hormone stimulus
Response to growth hormone stimulus
Transmembrane receptor protein tyrosine kinase signaling pathway
Enzyme linked receptor protein signaling pathway
Top 10 pathways from BiNGO pathway analysis of protein sub-network depicted in Figure 3 (d)
Regulation of body fluid levels
Response to wounding
Regulation of biological quality
Response to stress
Response to stimulus
Acute inflammatory response
In summary, we confirm prior observations demonstrating that on average, the level of CRP is higher in smokers than non-smokers, consistent with increased inflammation levels in smokers. We also extend findings by ourselves and others on the effect of smoking on peripheral mononuclear cell DNA methylation to show that continued smoking is associated with altered methylation signatures at AHRR and in pathways associated with immune, coagulation and CNS function.
The strength of the genome-wide findings with respect to the effects of smoking on DNA methylation indicates a need for a better longitudinal understanding of the response of the genome to tobacco smoke. In two consecutive prior studies of younger African American smokers with considerably less total cumulative tobacco smoke exposure, we demonstrated that the principal changes of the DNA methylation signature in response to smoking were in the xenobiotic response pathway regulated by the AHRR gene . In contrast, in this group of African American female smokers, many of whom have been smoking for 25 or more years, the magnitude of the remodeling of the DNA methylation signature is significantly greater and extends to a number of other genes and pathways including those relevant to inflammation, immune, coagulation and CNS function. At the single gene level, the findings at GFI1 (probe rank 24), F2RL3 (probe rank 49) and GPR15 (probe rank 1) serve to directly confirm prior findings by other investigators [5, 7, 14]. At a more integrated level, the implication of alterations of pathways involved in immune/inflammatory activities and coagulation (see Figure 3) is consistent with the known effects of smoking on risk for inflammatory related diseases and stroke [1, 15].
Since there is a strong association between cigarette smoking and innate immunity, we hypothesized that networks enriching for immune system/inflammation would be generated. Two such networks are depicted in Figures 3(a) and (b) that highlight the potential role of peripheral mononuclear cells in moderating vulnerability to some but not all complex illnesses with inflammatory components. The BiNGO analysis of network 3(a) indicated that the proteins are involved in the response to acute inflammation and complement activation. Specifically, MASP2 with a differential methylation adjusted p-value of 3.78 × 10-6 is known to be a protease that activates the complement cascade, enabling the clearing of pathogens. It has been shown that MASP2 is also associated with the recurrence and survival rate of cancer (colorectal cancer in specific) . Therefore, studying the effects of smoking in a weighted network approach allows us to speculate that patients with colorectal cancer may increase their survival chances and prevent the recurrence of cancer by avoiding smoking. In clinical settings, knowing smoking history could assist physicians in making better judgments of the efficacy of treatments.
Based on network 3(a), it can also be seen that MASP2 interacts strongly with MASP1 and FCN1. MASP1 is known to support MASP2 in the complement activation cascade  and this network validates that a strong interaction is present between these two serine proteases. MASP2, MASP1, and FCN1 play key roles in innate immunity. Aberrant methylation of these genes as a result of smoking could possibly impair an individual’s innate immunity, leading to the development, progression and recurrence of cancer. The network shown in Figure 3(b) is enriched for chemotaxis. Here, it can be noted that CCL13 with a differential methylation adjusted p-value of 0.00012 forms relatively strong interactions with MMP3, CCL2 and CCBP2. Again, similar to network 3(a), these interactions corroborate the role of smoking in the progression of inflammation.
Our work also builds on the prior work by Breitling and colleagues who showed the role of smoking associated differential methylation at coagulation factor II (thrombin) receptor like-3 (F2RL3) in moderating coronary artery disease, which also has a strong inflammatory component and whose presence is positively associated with CRP [3, 18]. Furthermore, the network in Figure 3(d) that enriches for coagulation showed F2 as the central protein, forming interactions with other proteins including F11 and CRB2. The miPALM analyses which generated pathways with important implications demonstrate that the pathways involving IL6R and F2RL3 may not be the only important inflammatory genes/pathways operant in peripheral mononuclear cells and suggest the possibility that smoking associated DNA methylation may act through a variety of mechanisms to increase systemic inflammation and vulnerability to some illnesses with an inflammatory component.
The time frame in which these cells begin to experience smoking associated inflammatory changes in DNA methylation is not clear. In our prior examinations of young adult smokers, these changes were not evident on genome-wide basis. But these subjects were relatively young (19-22 years old) at the time of ascertainment and the cumulative smoking exposure of even the heaviest of the smokers was not as extensive as that the lightest of the current smokers . Understanding when that shift to inflammatory processes occurs may be important to the allocation of smoking cessation resources. Although it is well established that in general smoking cessation decreases risk for adverse outcomes such as myocardial infarctions and cancer, this may not be true for all patients who quit smoking. According to the “pulmonary overflow” hypothesis, many of the inflammatory changes associated with smoking are not due to the direct effects of smoke but are rather secondary to the “spillover” of inflammatory cytokines from activated macrophages in the lung parenchyma. Because this inflammatory process tends to be self-perpetuating once it has been initiated, if “overflow” hypothesis is correct, stopping smoking in individuals with severe diseases may have little effect on subsequent inflammatory changes in the peripheral mononuclear cell compartment. If so, focusing smoking cessation resources on those smokers who do not have substantial diseases as evidenced by either peripheral or direct measures, may lead to more effective use of prevention dollars.
In conclusion, we reiterate that smoking contributes to inflammation which is consistent with increase in the average serum inflammatory protein, CRP. We also demonstrate that methylation residues at the AHRR locus are sensitive to cumulative smoking and could be potential smoking biomarkers. However, the main finding in this well characterized cohort of long-term smokers is that differential methylation between smokers and non-smokers are key in delineating epigenetically contextual pathways and to identify highly connected networks that could be used for smoking related health diagnostics like cancer and coronary heart disease. We show the extensive effects of smoking on peripheral mononuclear cell DNA methylation on pathways involved in coagulation, CNS and immune function. Understanding aberrant methylation caused by smoking in a network setting provides an integrated stage for understanding complex diseases like cancer while designing potential therapies. Finally, smoking is an important confounder and should be included in future diagnostic models to accurately understand diseases.
The participants provided written consent to participate in this study. The written consent form, consent procedure, and protocols pertaining to this study were approved by the Institutional Review Board at the University of Iowa, the University of Georgia and Iowa State University.
The subjects described in the current study are adult African American females from the states of Iowa and Georgia who participated in data and biomaterial collection during the most recent waves of the Family and Community Health Study (FACHS). The FACHS study is a longitudinal examination of factors affecting the health and health related outcomes of rural African American families . These adults were identified because they had an 11 or 12-year-old child who resided in eligible recruitment area in either Iowa or Georgia. These geographical zones of eligibility varied with respect to socioeconomic status and were not selected for any other characteristic except relative enrichment in African American residents. Lists of these potential caretaker-offspring dyads living in these areas were then generated by community officials in Georgia and school officials in Iowa.
To recruit these individuals in the study, potential participant families were selected randomly from these lists and were then sent an introductory letter of invitation followed by a phone call. Families expressing interest in the study were then provided with a complete description of the study and if still willing, consented. The protocols and procedures pertaining to this study were approved by the Institutional Review Board at the University of Iowa, the University of Georgia and Iowa State University.
Each participant in the FACHS study is interviewed every two years in their homes with a battery of structured and semi-structured instruments that characterize the medical and psychological status of each subject and environment in which the subject is embedded. The data and biomaterials used in this study were collected during the most recent wave (Wave 10) of the project. As part of that examination, subjects were phlebotomized to provide biomaterials for the current study with a subset additionally being interviewed with the tobacco module from the Semi-Structured Assessment for the Genetics of Alcoholism-II (SSAGA-II) . Individuals who were actively smoking were characterized as smokers and those denying any use of tobacco products were characterized as non-smokers.
Genome-wide methylation profiling
For the current study, peripheral mononuclear cell DNA samples from FACHS subjects were assessed for genome-wide methylation status using our standard protocols [6, 11]. In brief, DNA was prepared from ficoll purified peripheral mononuclear cell DNA pellets using a Qiagen (Valencia, CA) DNA Mini Kit. Measurement of genome-wide peripheral mononuclear cell DNA methylation status of FACHS participants was conducted using the Illumina (San Diego, CA) HumanMethylation450 Beadchip by the University of Minnesota Genome Center (Minneapolis, MN) using the protocol specified by the manufacturer . This array contains 485,577 probes recognizing CpG positions of at least 20,216 known transcripts, potential transcripts or CpG islands. Subjects were randomly assigned to 12 sample “slides” with groups of 8 slides being bisulfite converted in a single batch. Four replicates of the same DNA were also included to monitor slide-to-slide and batch bisulfite conversion variability. The resulting data were inspected for complete bisulfite conversion and average beta values for each targeted CpG residue was determined using the Illumina Genome Studio Methylation Module, Version 3.2. The beta value at a CpG locus is the ratio between the intensity of the methylated probe to the sum of intensities of the methylated and unmethylated probes, also known as the total probe intensities. The resulting data was then processed using a custom PERL script to remove beta values with a detection p-value greater than 0.05.
Accounting for mixed cell population
In order to account for mixed cell population in peripheral mononuclear cell mixture, a regression calibration approach similar to that developed by Houseman and colleagues was performed . Since the method developed by Houseman and colleagues was based on the Illumina HumanMethylation 27K BeadChip, we instead utilized the Illumina HumanMethylation 450K BeadChip data of purified cells (CD4+ T cells, CD8+ T cells, CD14+ monocytes, CD19+ B cells and CD56+ Natural Killer cells) contributed by Reinius and colleagues (GEO database under accession number GSE35069) and performed a regression in MethLAB, Version 1.5 [22, 23]. The Benjamini-Hochberg method at a 0.05 significance level was used for genome-wide correction. This resulted in 815 sites being significantly differentially methylated with respect to cell types. The top 100 sites were then imported into JMP Version 10 (SAS Institute, Cary, USA) and principal components analysis (PCA) was performed. Based on the Scree plot of this PCA, the first 5 factors were chosen to be included in downstream genome-wide analyses to control for signal as a result of mixed cell population.
Analysis of genome-wide DNA methylation data
Genome-wide methylation data analyses were conducted using the R package, MethLAB, Version 1.5 as previously described . Briefly, the genome-wide differential methylation between smokers and non-smokers were evaluated while controlling for slide, plate and mixed cell population effects. In MethLAB, the phenotype dependent differential methylation p-value at each CpG site is determined using a linear model . Genome-wide false discovery rate was corrected using the method of Benjamini-Hochberg using the MethLAB default significance level of 0.05 .
Validation of differential methylation at cg05575921
An independent validation of the differential methylation at cg05575921 was performed on a subset of 62 individuals using quantitative-PCR primer probe set obtained from Behavioral Diagnostics Inc (Iowa City, IA). Quantitative-PCR was performed with an ABI7900 (Life Technologies, New York) using provided reagents, internal controls and standard manufacturer protocols. Interpolation of internal control methylation was used to deduce percent methylation of each sample. Pearson correlation was performed to determine the correlation between the methylation level obtained via the Illumina platform and that obtained from the quantitative-PCR. Also, the relationship between cg05575921 methylation from the quantitative-PCR and smoking status was assessed using an analysis of variance (ANOVA) test.
Protein network analyses
Protein sub-networks pertinent to smoking were identified using our previously described methods . In brief, the genome-wide methylation data were first divided into bins based on the position of the CpG site with respect to the gene. We then focused the analyses on CpG probes that mapped to the promoter region of the gene. This is due to the known influence of promoter methylation on transcription . This resulted in a total of 22,609 genes or potential gene transcripts with at least one CpG probe mapping to their putative promoter region. Using the data from all probes mapping to a given promoter region, two promoter methylation vectors were constructed for each gene; one for smokers and the other for non-smokers. Using the statistical programming language, R , for each gene, an unpaired t-test was performed between the two vectors resulting in p-values of differential methylation between smokers and non-smokers, that were adjusted for multiple comparisons using the Benjamini-Hochberg method .
Human protein-protein interaction data were downloaded from the iRefIndex database . This database provides a comprehensive and non-redundant version of protein-protein interactions (PPI) available in several other primary interaction databases including BIND, BioGRID, CORUM, DIP, HPRD, IntAct, MINT, MPact, MPPI, and OPHID. The 22,609 genes were then compared to the non-redundant PPI that were extracted. There were 29,912 interactions that mapped to the list of genes. In order to construct an edge-weighted network, the subset of interactions alongside the adjusted p-values was used to calculate the PPI edge weights.
Where w ij is the edge weight between genes i and j, p’ i is the adjusted p-value of the t-test between smokers and non-smokers for gene i, p’ j is the corresponding adjusted p-value for gene j, and p’ min is the minimum adjusted p-value among the 22,609 adjusted p-values calculated for the respective 22,609 genes.
Computed edge weights were assigned to their respective PPI and the resulting output was used to identify significant sub-networks using the miPALM algorithm implemented in JAVA. Briefly, this algorithm, which uses a greedy search method, identifies highly connected sub-networks using a novel graph theoretic measure, parametric local modularity [25, 29].
Statistically significant sub-networks at the 0.05 significance level were imported into Cytoscape, a bioinformatics platform that allows visualization and analysis of biological networks . The Cytoscape plugin BiNGO was used to determine Gene Ontology (GO) categories of statistically overrepresented genes . In addition, the adjusted p-values of individual genes were imported as a node attribute and visualized as a color gradient. Edge weights were used as an edge attribute to represent the strength of interactions illustrated by the width of the lines connecting nodes.
Enzyme linked immunosorbent assay (ELISA)
Sera were obtained using serum separator tubes. After centrifugation, sera were frozen at -80°C until use. ELISA assessments of CRP and IL6R serum levels were conducted using a DuoKit from R&D Systems (Minneapolis, USA) according to manufacturer’s recommendations.
The analyses of clinical and serological data were analyzed using the suite of general linear model algorithms (e.g. ANOVA, Bivariate, ordinal logistic regression) contained in JMP, Version 10 (SAS Institute, Cary, USA). Where indicated, the portion of the variance explained by the model tested is given by the Adjusted R square (Adj R2).
Availability of supporting data
The Illumina Human Methylation 450K BeadChip data used in this study has been deposited to the Gene Expression Omnibus (GEO) database. The GEO accession number is GSE53045.
The work in this study was supported by National Institutes of Health grants MH080898 and R21DA034457 to Dr. Philibert and 5R01HD030588-16A1 to Dr. Brody. Additional support for these studies was derived from the Center for Contextual Genetics and Prevention Science (Grant Number P30 DA027827, GB) funded by the National Institute on Drug Abuse. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
- Centers for Disease Control: Cigarette smoking among adults-United States 2006. Morbidity and Mortality Weekly Report. 2007, 56: 1157-1161.Google Scholar
- Centers for Disease C, Prevention: Smoking-attributable mortality, years of potential life lost, and productivity losses--United States, 2000-2004. MMWR Morb Mortal Wkly Rep. 2008, 57: 1226-1228.Google Scholar
- Breitling LP, Salzmann K, Rothenbacher D, Burwinkel B, Brenner H: Smoking, F2RL3 methylation, and prognosis in stable coronary heart disease. Eur Heart J. 2012, 33: 2841-2848. 10.1093/eurheartj/ehs091.PubMedView ArticleGoogle Scholar
- Philibert RA, Beach SR, Gunter TD, Brody GH, Madan A, Gerrard M: The effect of smoking on MAOA promoter methylation in DNA prepared from lymphoblasts and whole blood. Am J Med Genet B Neuropsychiatr Genet. 2010, 153B: 619-628.PubMedGoogle Scholar
- Breitling LP, Yang R, Korn B, Burwinkel B, Brenner H: Tobacco-smoking-related differential DNA methylation: 27K discovery and replication. American journal of human genetics. 2011, 88: 450-457. 10.1016/j.ajhg.2011.03.003.PubMed CentralPubMedView ArticleGoogle Scholar
- Monick MM, Beach SR, Plume J, Sears R, Gerrard M, Brody GH, Philibert RA: Coordinated changes in AHRR methylation in lymphoblasts and pulmonary macrophages from smokers. Am J Med Genet B Neuropsychiatr Genet. 2012, 159B: 141-151. 10.1002/ajmg.b.32021.PubMed CentralPubMedView ArticleGoogle Scholar
- Joubert BR, Håberg SE, Nilsen RM, Wang X, Vollset SE, Murphy SK, Huang Z, Hoyo C, Midttun O, Cupul-Uicab LA, et al: 450K epigenome-wide scan identifies differential DNA methylation in newborns related to maternal smoking during pregnancy. Environ Health Perspect. 2012, 120: 1425-1431. 10.1289/ehp.1205412.PubMed CentralPubMedView ArticleGoogle Scholar
- Zeilinger S, Kühnel B, Klopp N, Baurecht H, Kleinschmidt A, Gieger C, Weidinger S, Lattka E, Adamski J, Peters A, et al: Tobacco smoking leads to extensive genome-wide changes in DNA methylation. PLoS One. 2013, 8: e63812-10.1371/journal.pone.0063812.PubMed CentralPubMedView ArticleGoogle Scholar
- Huang RY, Chen GG: Cigarette smoking, cyclooxygenase-2 pathway and cancer. Biochim Biophys Acta. 2011, 1815: 158-169.PubMedGoogle Scholar
- Max W, Sung HY, Tucker LY, Stark B: The disproportionate cost of smoking for African Americans in California. Am J Public Health. 2010, 100: 152-158. 10.2105/AJPH.2008.149542.PubMed CentralPubMedView ArticleGoogle Scholar
- Philibert RA, Beach SR, Brody GH: Demethylation of the aryl hydrocarbon receptor repressor as a biomarker for nascent smokers. Epigenetics. 2012, 7: 1331-1338. 10.4161/epi.22520.PubMed CentralPubMedView ArticleGoogle Scholar
- Mimura J, Ema M, Sogawa K, Fujii-Kuriyama Y: Identification of a novel mechanism of regulation of Ah (dioxin) receptor function. Genes Dev. 1999, 13: 20-25. 10.1101/gad.13.1.20.PubMed CentralPubMedView ArticleGoogle Scholar
- Oshima M, Mimura J, Yamamoto M, Fujii-Kuriyama Y: Molecular mechanism of transcriptional repression of AhR repressor involving ANKRA2, HDAC4, and HDAC5. Biochem Biophys Res Commun. 2007, 364: 276-282. 10.1016/j.bbrc.2007.09.131.PubMedView ArticleGoogle Scholar
- Shenker NS, Polidoro S, van Veldhoven K, Sacerdote C, Ricceri F, Birrell MA, Belvisi MG, Brown R, Vineis P, Flanagan JM: Epigenome-wide association study in the European prospective investigation into cancer and nutrition (EPIC-Turin) identifies novel genetic loci associated with smoking. Hum Mol Genet. 2012, 22: 843-851.PubMedView ArticleGoogle Scholar
- Rose-John S, Waetzig GH, Scheller J, Grötzinger J, Seegert D: The IL-6/sIL-6R complex as a novel target for therapeutic approaches. Expert Opinion on Therapeutic Targets. 2007, 11: 613-624. 10.1517/14728184.108.40.2063.PubMedView ArticleGoogle Scholar
- Ytting H, Christensen IJ, Thiel S, Jensenius JC, Nielsen HJ: Serum mannan-binding lectin-associated serine protease 2 levels in colorectal cancer: relation to recurrence and mortality. Clin Cancer Res. 2005, 11: 1441-1446. 10.1158/1078-0432.CCR-04-1272.PubMedView ArticleGoogle Scholar
- Dobó J, Harmat V, Beinrohr L, Sebestyén E, Závodszky P, Gál P: MASP-1, a promiscuous complement protease: structure of its catalytic region reveals the basis of its broad specificity. J Immunol. 2009, 183: 1207-1214. 10.4049/jimmunol.0901141.PubMedView ArticleGoogle Scholar
- Pearson TA, Mensah GA, Alexander RW, Anderson JL, Cannon RO, Criqui M, Fadl YY, Fortmann SP, Hong Y, Myers GL, et al: Markers of inflammation and cardiovascular disease application to clinical and public health practice - a statement for healthcare professionals from the centers for disease control and prevention and the American Heart Association. Circulation. 2003, 107: 499-511. 10.1161/01.CIR.0000052939.59093.45.PubMedView ArticleGoogle Scholar
- Cutrona CE, Russell DW, Brown PA, Clark LA, Hessling RM, Gardner KA: Neighborhood context, personality, and stressful life events as predictors of depression among African American women. Journal of Abnormal Psychology. 2005, 114: 3-15.PubMed CentralPubMedView ArticleGoogle Scholar
- Bucholz KK, Cadoret R, Cloninger CR, Dinwiddie SH, Hesselbrock VM, Nurnberger JI, Reich T, Schmidt I, Schuckit MA: A new, semi-structured psychiatric interview for use in genetic linkage studies: a report on the reliability of the SSAGA. J Stud Alcohol. 1994, 55: 149-158.PubMedView ArticleGoogle Scholar
- Houseman EA, Accomando WP, Koestler DC, Christensen BC, Marsit CJ, Nelson HH, Wiencke JK, Kelsey KT: DNA methylation arrays as surrogate measures of cell mixture distribution. BMC Bioinformatics. 2012, 13: 86-10.1186/1471-2105-13-86.PubMed CentralPubMedView ArticleGoogle Scholar
- Reinius LE, Acevedo N, Joerink M, Pershagen G, Dahlén SE, Greco D, Söderhäll C, Scheynius A, Kere J: Differential DNA methylation in purified human blood cells: implications for cell lineage and studies on disease susceptibility. PLoS One. 2012, 7: e41361-10.1371/journal.pone.0041361.PubMed CentralPubMedView ArticleGoogle Scholar
- Kilaru V, Barfield RT, Schroeder JW, Smith AK, Conneely KN: MethLAB: a graphical user interface package for the analysis of array-based DNA methylation data. Epigenetics. 2012, 7: 225-229. 10.4161/epi.7.3.19284.PubMed CentralPubMedView ArticleGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate - a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society Series B-Methodological. 1995, 57: 289-300.Google Scholar
- Kim J, Gao L, Tan K: Multi-analyte network markers for tumor prognosis. PLoS One. 2012, 7: e52973-10.1371/journal.pone.0052973.PubMed CentralPubMedView ArticleGoogle Scholar
- Suzuki MM, Bird A: DNA methylation landscapes: provocative insights from epigenomics. Nat Rev Genet. 2008, 9: 465-476. 10.1038/nrg2341.PubMedView ArticleGoogle Scholar
- Team RC: R: A language and environment for statistical computing. 2012, Vienna, Austria: R Foundation for Statistical ComputingGoogle Scholar
- Razick S, Magklaras G, Donaldson IM: iRefIndex: a consolidated protein interaction database with provenance. Bmc Bioinformatics. 2008, 9: 405-10.1186/1471-2105-9-405.PubMed CentralPubMedView ArticleGoogle Scholar
- Kim J, Tan K: Discover protein complexes in protein-protein interaction networks using parametric local modularity. BMC Bioinformatics. 2010, 11: 521-10.1186/1471-2105-11-521.PubMed CentralPubMedView ArticleGoogle Scholar
- Smoot ME, Ono K, Ruscheinski J, Wang PL, Ideker T: Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics. 2011, 27: 431-432. 10.1093/bioinformatics/btq675.PubMed CentralPubMedView ArticleGoogle Scholar
- Maere S, Heymans K, Kuiper M: BiNGO: a cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics. 2005, 21: 3448-3449. 10.1093/bioinformatics/bti551.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.