The dynamics of smoking-related disturbed methylation: a two time-point study of methylation change in smokers, non-smokers and former smokers.

BACKGROUND: The evidence for epigenome-wide associations between smoking and DNA methylation continues to grow through cross-sectional studies. However, few large-scale investigations have explored the associations using observations for individuals at multiple time-points. Here, through the use of the Illumina 450K BeadChip and data collected at two time-points separated by approximately 7 years, we investigate changes in methylation over time associated with quitting smoking or remaining a former smoker, and those associated with continued smoking. RESULTS: Our results indicate that after quitting smoking the most rapid reversion of altered methylation occurs within the first two decades, with reversion rates related to the initial differences in methylation. For 52 CpG sites, the change in methylation from baseline to follow-up is significantly different for former smokers relative to the change for never smokers (lowest p-value 3.61 x 10-39 for cg26703534, gene AHRR). Most of these sites' respective regions have been previously implicated in smoking-associated diseases. Despite the early rapid change, dynamism of methylation appears greater in former smokers vs never smokers even four decades after cessation. Furthermore, our study reveals the heterogeneous effect of continued smoking: the methylation levels of some loci further diverge between smokers and non-smokers, while others re-approach. Though intensity of smoking habit appears more significant than duration, results remain inconclusive. CONCLUSIONS: This study improves the understanding of the dynamic link between cigarette smoking and methylation, revealing the continued fluctuation of methylation levels decades after smoking cessation and demonstrating that continuing smoking can have an array of effects. The results can facilitate insights into the molecular mechanisms behind smoking-induced disturbed methylation, improving the possibility for development of biomarkers of past smoking behavior and increasing the understanding of the molecular path from exposure to disease.

Selected single-nucleotide polymorphisms in FOXE1, SERPINA5, FTO, EVPL, TICAM1 and SCARB1 are associated with papillary and follicular thyroid cancer risk: replication study in a German population.

Several single-nucleotide polymorphisms (SNPs) have been associated with papillary and follicular thyroid cancer (PTC and FTC, respectively) risk, but few have replicated. After analyzing 17525 tag SNPs in 1129 candidate genes, we found associations with PTC risk in SERPINA5, FTO, HEMGN (near FOXE1) and other genes. Here, we report results from a replication effort in a large independent PTC/FTC case-control study conducted in Germany. We evaluated the best tagging SNPs from our previous PTC study and additionally included SNPs in or near FOXE1 and NKX2-1 genes, known susceptibility loci for thyroid cancer. We genotyped 422 PTC and 130 FTC cases and 752 controls recruited from three German clinical centers. We used polytomous logistic regression to simultaneously estimate PTC and FTC associations for 79 SNPs based on log-additive models. We assessed effect modification by body mass index (BMI), gender and age for all SNPs, and selected SNP by SNP interactions. We confirmed associations with PTC and SNPs in FOXE1/HEMGN, SERPINA5 (rs2069974), FTO (rs8047395), EVPL (rs2071194), TICAM1 (rs8120) and SCARB1 (rs11057820) genes. We found associations with SNPs in FOXE1, SERPINA5, FTO, TICAM1 and HSPA6 and FTC. We found two significant interactions between FTO (rs8047395) and BMI (P = 0.0321) and between TICAM1 (rs8120) and FOXE1 (rs10984377) (P = 0.0006). Besides the known associations with FOXE1 SNPs, we confirmed additional PTC SNP associations reported previously. We also found several new associations with FTC risk and noteworthy interactions. We conclude that multiple variants and host factors might interact in complex ways to increase risk of PTC and FTC.

Genetic variation influencing DNA methylation provides insights into molecular mechanisms regulating genomic function.

We determined the relationships between DNA sequence variation and DNA methylation using blood samples from 3,799 Europeans and 3,195 South Asians. We identify 11,165,559 SNP-CpG associations (methylation quantitative trait loci (meQTL), P < 10-14), including 467,915 meQTL that operate in trans. The meQTL are enriched for functionally relevant characteristics, including shared chromatin state, High-throuhgput chromosome conformation interaction, and association with gene expression, metabolic variation and clinical traits. We use molecular interaction and colocalization analyses to identify multiple nuclear regulatory pathways linking meQTL loci to phenotypic variation, including UBASH3B (body mass index), NFKBIE (rheumatoid arthritis), MGA (blood pressure) and COMMD7 (white cell counts). For rs6511961 , chromatin immunoprecipitation followed by sequencing (ChIP-seq) validates zinc finger protein (ZNF)333 as the likely trans acting effector protein. Finally, we used interaction analyses to identify population- and lineage-specific meQTL, including rs174548 in FADS1, with the strongest effect in CD8+ T cells, thus linking fatty acid metabolism with immune dysregulation and asthma. Our study advances understanding of the potential pathways linking genetic variation to human phenotype.

DNA methylation and lipid metabolism: an EWAS of 226 metabolic measures.

BACKGROUND: The discovery of robust and trans-ethnically replicated DNA methylation markers of metabolic phenotypes, has hinted at a potential role of epigenetic mechanisms in lipid metabolism. However, DNA methylation and the lipid compositions and lipid concentrations of lipoprotein sizes have been scarcely studied. Here, we present an epigenome-wide association study (EWAS) (N = 5414 total) of mostly lipid-related metabolic measures, including a fine profiling of lipoproteins. As lipoproteins are the main players in the different stages of lipid metabolism, examination of epigenetic markers of detailed lipoprotein features might improve the diagnosis, prognosis, and treatment of metabolic disturbances. RESULTS: We conducted an EWAS of leukocyte DNA methylation and 226 metabolic measurements determined by nuclear magnetic resonance spectroscopy in the population-based KORA F4 study (N = 1662) and replicated the results in the LOLIPOP, NFBC1966, and YFS cohorts (N = 3752). Follow-up analyses in the discovery cohort included investigations into gene transcripts, metabolic-measure ratios for pathway analysis, and disease endpoints. We identified 161 associations (p value < 4.7 × 10-10), covering 16 CpG sites at 11 loci and 57 metabolic measures. Identified metabolic measures were primarily medium and small lipoproteins, and fatty acids. For apolipoprotein B-containing lipoproteins, the associations mainly involved triglyceride composition and concentrations of cholesterol esters, triglycerides, free cholesterol, and phospholipids. All associations for HDL lipoproteins involved triglyceride measures only. Associated metabolic measure ratios, proxies of enzymatic activity, highlight amino acid, glucose, and lipid pathways as being potentially epigenetically implicated. Five CpG sites in four genes were associated with differential expression of transcripts in blood or adipose tissue. CpG sites in ABCG1 and PHGDH showed associations with metabolic measures, gene transcription, and metabolic measure ratios and were additionally linked to obesity or previous myocardial infarction, extending previously reported observations. CONCLUSION: Our study provides evidence of a link between DNA methylation and the lipid compositions and lipid concentrations of different lipoprotein size subclasses, thus offering in-depth insights into well-known associations of DNA methylation with total serum lipids. The results support detailed profiling of lipid metabolism to improve the molecular understanding of dyslipidemia and related disease mechanisms.

Methylation of the FKBP5 gene in association with FKBP5 genotypes, childhood maltreatment and depression.

DNA methylation of the FKBP5 gene is assumed to alter FKBP5 expression and hence the synthesis of the FK506 binding protein 51, a central element of a genomic negative feedback loop for glucocorticoid receptor signaling. The present study aimed to replicate and extend previously reported influences of FKBP5 genotypes, childhood maltreatment and depression on methylation levels of five CpG sites in intron 7 of the FKBP5 gene in a large population-based sample. Besides the single nucleotide polymorphism (SNP) rs1360780, associations of the FKBP5 methylation with 22 other, unlinked FKBP5 SNPs as well as associations between FKBP5 methylation levels and transcription levels were investigated. Using whole-blood methylation of 3965 subjects of the Study of Health in Pomerania (SHIP) reduced methylation levels in TT allele carriers of rs1360780 (OR = 0.975, p = .005) and currently depressed subjects (OR = 0.995, p = 0.005) were found. Further, an impact of two yet undescribed SNPs (rs6910300, rs7771727) on methylation levels was observed. However, main and interactive effects for childhood maltreatment and lifetime major depressive disorder observed in previous studies could not be replicated. Finally, FKBP5 methylation levels were not related to FKBP5 transcription levels in whole blood. Thus, the present study verified the associations of FKBP5 genotypes and state depression on the FKBP5 methylation levels of five CpG sites in intron 7. However, FKBP5 methylation of these five CpG sites could not be validated as a valuable clinical biomarker for biological long-term effects of childhood maltreatment or lifetime depression.

Analysis of repeated leukocyte DNA methylation assessments reveals persistent epigenetic alterations after an incident myocardial infarction.

BACKGROUND: Most research into myocardial infarctions (MIs) have focused on preventative efforts. For survivors, the occurrence of an MI represents a major clinical event that can have long-lasting consequences. There has been little to no research into the molecular changes that can occur as a result of an incident MI. Here, we use three cohorts to identify epigenetic changes that are indicative of an incident MI and their association with gene expression and metabolomics. RESULTS: Using paired samples from the KORA cohort, we screened for DNA methylation loci (CpGs) whose change in methylation is potentially indicative of the occurrence of an incident MI between the baseline and follow-up exams. We used paired samples from the NAS cohort to identify 11 CpGs which were predictive in an independent cohort. After removing two CpGs associated with medication usage, we were left with an "epigenetic fingerprint" of MI composed of nine CpGs. We tested this fingerprint in the InCHIANTI cohort where it moderately discriminated incident MI occurrence (AUC = 0.61, P = 6.5 × 10-3). Returning to KORA, we associated the epigenetic fingerprint loci with cis-gene expression and integrated it into a gene expression-metabolomic network, which revealed links between the epigenetic fingerprint CpGs and branched chain amino acid (BCAA) metabolism. CONCLUSIONS: There are significant changes in DNA methylation after an incident MI. Nine of these CpGs show consistent changes in multiple cohorts, significantly discriminate MI in independent cohorts, and were independent of medication usage. Integration with gene expression and metabolomics data indicates a link between MI-associated epigenetic changes and BCAA metabolism.

Anxiety Associated Increased CpG Methylation in the Promoter of Asb1: A Translational Approach Evidenced by Epidemiological and Clinical Studies and a Murine Model.

Epigenetic regulation in anxiety is suggested, but evidence from large studies is needed. We conducted an epigenome-wide association study (EWAS) on anxiety in a population-based cohort and validated our finding in a clinical cohort as well as a murine model. In the KORA cohort, participants (n=1522, age 32-72 years) were administered the Generalized Anxiety Disorder (GAD-7) instrument, whole blood DNA methylation was measured (Illumina 450K BeadChip), and circulating levels of hs-CRP and IL-18 were assessed in the association between anxiety and methylation. DNA methylation was measured using the same instrument in a study of patients with anxiety disorders recruited at the Max Planck Institute of Psychiatry (MPIP, 131 non-medicated cases and 169 controls). To expand our mechanistic understanding, these findings were reverse translated in a mouse model of acute social defeat stress. In the KORA study, participants were classified according to mild, moderate, or severe levels of anxiety (29.4%/6.0%/1.5%, respectively). Severe anxiety was associated with 48.5% increased methylation at a single CpG site (cg12701571) located in the promoter of the gene encoding Asb1 (ß-coefficient=0.56 standard error (SE)=0.10, p (Bonferroni)=0.005), a protein hypothetically involved in regulation of cytokine signaling. An interaction between IL-18 and severe anxiety with methylation of this CpG cite showed a tendency towards significance in the total population (p=0.083) and a significant interaction among women (p=0.014). Methylation of the same CpG was positively associated with Panic and Agoraphobia scale (PAS) scores (ß=0.005, SE=0.002, p=0.021, n=131) among cases in the MPIP study. In a murine model of acute social defeat stress, Asb1 gene expression was significantly upregulated in a tissue-specific manner (p=0.006), which correlated with upregulation of the neuroimmunomodulating cytokine interleukin 1 beta. Our findings suggest epigenetic regulation of the stress-responsive Asb1 gene in anxiety-related phenotypes. Further studies are necessary to elucidate the causal direction of this association and the potential role of Asb1-mediated immune dysregulation in anxiety disorders.

Evidence for Stress-like Alterations in the HPA-Axis in Women Taking Oral Contraceptives.

Using oral contraceptives has been implicated in the aetiology of stress-related disorders like depression. Here, we followed the hypothesis that oral contraceptives deregulate the HPA-axis by elevating circulating cortisol levels. We report for a sample of 233 pre-menopausal women increased circulating cortisol levels in those using oral contraceptives. For women taking oral contraceptives, we observed alterations in circulating phospholipid levels and elevated triglycerides and found evidence for increased glucocorticoid signalling as the transcript levels of the glucocorticoid-regulated genes DDIT4 and FKBP5 were increased in whole blood. The effects were statistically mediated by cortisol. The associations of oral contraceptives with higher FKBP5 mRNA and altered phospholipid levels were modified by rs1360780, a genetic variance implicated in psychiatric diseases. Accordingly, the methylation pattern of FKBP5 intron 7 was altered in women taking oral contraceptives depending on the rs1360780 genotype. Moreover, oral contraceptives modified the association of circulating cortisol with depressive symptoms, potentially explaining conflicting results in the literature. Finally, women taking oral contraceptives displayed smaller hippocampal volumes than non-using women. In conclusion, the integrative analyses of different types of physiological data provided converging evidence indicating that oral contraceptives may cause effects analogous to chronic psychological stressors regarding the regulation of the HPA axis.

DNA methylation of lipid-related genes affects blood lipid levels.

BACKGROUND: Epigenetic mechanisms might be involved in the regulation of interindividual lipid level variability and thus may contribute to the cardiovascular risk profile. The aim of this study was to investigate the association between genome-wide DNA methylation and blood lipid levels high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, triglycerides, and total cholesterol. Observed DNA methylation changes were also further analyzed to examine their relationship with previous hospitalized myocardial infarction. METHODS AND RESULTS: Genome-wide DNA methylation patterns were determined in whole blood samples of 1776 subjects of the Cooperative Health Research in the Region of Augsburg F4 cohort using the Infinium HumanMethylation450 BeadChip (Illumina). Ten novel lipid-related CpG sites annotated to various genes including ABCG1, MIR33B/SREBF1, and TNIP1 were identified. CpG cg06500161, located in ABCG1, was associated in opposite directions with both high-density lipoprotein cholesterol (ß coefficient=-0.049; P=8.26E-17) and triglyceride levels (ß=0.070; P=1.21E-27). Eight associations were confirmed by replication in the Cooperative Health Research in the Region of Augsburg F3 study (n=499) and in the Invecchiare in Chianti, Aging in the Chianti Area study (n=472). Associations between triglyceride levels and SREBF1 and ABCG1 were also found in adipose tissue of the Multiple Tissue Human Expression Resource cohort (n=634). Expression analysis revealed an association between ABCG1 methylation and lipid levels that might be partly mediated by ABCG1 expression. DNA methylation of ABCG1 might also play a role in previous hospitalized myocardial infarction (odds ratio, 1.15; 95% confidence interval=1.06-1.25). CONCLUSIONS: Epigenetic modifications of the newly identified loci might regulate disturbed blood lipid levels and thus contribute to the development of complex lipid-related diseases.

Genistein ameliorates hyperglycemia in a mouse model of nongenetic type 2 diabetes.

While peripheral insulin resistance is common during obesity and aging in mice and people, the progression to type 2 diabetes (T2D) is largely due to loss of ß-cell mass and function through apoptosis. We recently reported that genistein, a soy derived isoflavone, can improve glycemic control and ß-cell function in insulin-deficient diabetic mice. However, whether it can prevent ß-cell loss and diabetes in T2D mice is unknown. Our current study aimed to investigate the effect of dietary supplemented genistein in a nongenetic T2D mouse model. Nongenetic, middle-aged obese diabetic mice were generated by high fat diet and a low dose of streptozotocin injection. The effect of dietary supplementation of genistein on glycemic control and ß-cell mass and function was determined. Dietary intake of genistein (250 mg·kg(-1) diet) improved hyperglycemia, glucose tolerance, and blood insulin level in obese diabetic mice, whereas it did not affect body weight gain, food intake, fat deposit, plasma lipid profile, and peripheral insulin sensitivity. Genistein increased the number of insulin-positive ß-cell in islets, promoted islet ß-cell survival, and preserved islet mass. In conclusion, dietary intake of genistein could prevent T2D via a direct protective action on ß-cells without alteration of periphery insulin sensitivity.