Age-related changes in DNA methylation have been implicated in cellular senescence and longevity, yet the causes and functional consequences of these variants remain unclear. To elucidate the role of age-related epigenetic changes in healthy ageing and potential longevity, we tested for association between whole-blood DNA methylation patterns in 172 female twins aged 32 to 80 with age and age-related phenotypes. Twin-based DNA methylation levels at 26,690 CpG-sites showed evidence for mean genome-wide heritability of 18%, which was supported by the identification of 1,537 CpG-sites with methylation QTLs in cis at FDR 5%. We performed genome-wide analyses to discover differentially methylated regions (DMRs) for sixteen age-related phenotypes (ap-DMRs) and chronological age (a-DMRs). Epigenome-wide association scans (EWAS) identified age-related phenotype DMRs (ap-DMRs) associated with LDL (STAT5A), lung function (WT1), and maternal longevity (ARL4A, TBX20). In contrast, EWAS for chronological age identified hundreds of predominantly hyper-methylated age DMRs (490 a-DMRs at FDR 5%), of which only one (TBX20) was also associated with an age-related phenotype. Therefore, the majority of age-related changes in DNA methylation are not associated with phenotypic measures of healthy ageing in later life. We replicated a large proportion of a-DMRs in a sample of 44 younger adult MZ twins aged 20 to 61, suggesting that a-DMRs may initiate at an earlier age. We next explored potential genetic and environmental mechanisms underlying a-DMRs and ap-DMRs. Genome-wide overlap across cis-meQTLs, genotype-phenotype associations, and EWAS ap-DMRs identified CpG-sites that had cis-meQTLs with evidence for genotype–phenotype association, where the CpG-site was also an ap-DMR for the same phenotype. Monozygotic twin methylation difference analyses identified one potential environmentally-mediated ap-DMR associated with total cholesterol and LDL (CSMD1). Our results suggest that in a small set of genes DNA methylation may be a candidate mechanism of mediating not only environmental, but also genetic effects on age-related phenotypes.

While there have been studies exploring regulatory variation in one or more tissues, the complexity of tissue-specificity in multiple primary tissues is not yet well understood. We explore in depth the role of cis-regulatory variation in three human tissues: lymphoblastoid cell lines (LCL), skin, and fat. The samples (156 LCL, 160 skin, 166 fat) were derived simultaneously from a subset of well-phenotyped healthy female twins of the MuTHER resource. We discover an abundance of cis-eQTLs in each tissue similar to previous estimates (858 or 4.7% of genes). In addition, we apply factor analysis (FA) to remove effects of latent variables, thus more than doubling the number of our discoveries (1,822 eQTL genes). The unique study design (Matched Co-Twin Analysis—MCTA) permits immediate replication of eQTLs using co-twins (93%–98%) and validation of the considerable gain in eQTL discovery after FA correction. We highlight the challenges of comparing eQTLs between tissues. After verifying previous significance threshold-based estimates of tissue-specificity, we show their limitations given their dependency on statistical power. We propose that continuous estimates of the proportion of tissue-shared signals and direct comparison of the magnitude of effect on the fold change in expression are essential properties that jointly provide a biologically realistic view of tissue-specificity. Under this framework we demonstrate that 30% of eQTLs are shared among the three tissues studied, while another 29% appear exclusively tissue-specific. However, even among the shared eQTLs, a substantial proportion (10%–20%) have significant differences in the magnitude of fold change between genotypic classes across tissues. Our results underline the need to account for the complexity of eQTL tissue-specificity in an effort to assess consequences of such variants for complex traits.

Epigenetic modifications such as DNA methylation play a key role in gene regulation and disease susceptibility. However, little is known about the genome-wide frequency, localization, and function of methylation variation and how it is regulated by genetic and environmental factors. We utilized the Multiple Tissue Human Expression Resource (MuTHER) and generated Illumina 450K adipose methylome data from 648 twins. We found that individual CpGs had low variance and that variability was suppressed in promoters. We noted that DNA methylation variation was highly heritable (h(2)median = 0.34) and that shared environmental effects correlated with metabolic phenotype-associated CpGs. Analysis of methylation quantitative-trait loci (metQTL) revealed that 28% of CpGs were associated with nearby SNPs, and when overlapping them with adipose expression quantitative-trait loci (eQTL) from the same individuals, we found that 6% of the loci played a role in regulating both gene expression and DNA methylation. These associations were bidirectional, but there were pronounced negative associations for promoter CpGs. Integration of metQTL with adipose reference epigenomes and disease associations revealed significant enrichment of metQTL overlapping metabolic-trait or disease loci in enhancers (the strongest effects were for high-density lipoprotein cholesterol and body mass index [BMI]). We followed up with the BMI SNP rs713586, a cg01884057 metQTL that overlaps an enhancer upstream of ADCY3, and used bisulphite sequencing to refine this region. Our results showed widespread population invariability yet sequence dependence on adipose DNA methylation but that incorporating maps of regulatory elements aid in linking CpG variation to gene regulation and disease risk in a tissue-dependent manner.

Smoking is a major risk factor in many diseases. Genome wide association studies have linked genes for nicotine dependence and smoking behavior to increased risk of cardiovascular, pulmonary, and malignant diseases. We conducted an epigenome wide association study in peripheral-blood DNA in 464 individuals (22 current smokers and 263 ex-smokers), using the Human Methylation 450 K array. Upon replication in an independent sample of 356 twins (41 current and 104 ex-smokers), we identified 30 probes in 15 distinct loci, all of which reached genome-wide significance in the combined analysis P < 5 × 10(-8). All but one probe (cg17024919) remained significant after adjusting for blood cell counts. We replicated all 9 known loci and found an independent signal at CPOX near GPR15. In addition, we found 6 new loci at PRSS23, AVPR1B, PSEN2, LINC00299, RPS6KA2, and KIAA0087. Most of the lead probes (13 out of 15) associated with cigarette smoking, overlapped regions of open chromatin (FAIRE and DNaseI hypersensitive sites) or/and H3K27Ac peaks (ENCODE data set), which mark regulatory elements. The effect of smoking on DNA methylation was partially reversible upon smoking cessation for longer than 3 months. We report the first statistically significant interaction between a SNP (rs2697768) and cigarette smoking on DNA methylation (cg03329539). We provide evidence that the metSNP for cg03329539 regulates expression of the CHRND gene located circa 95 Kb downstream of the methylation site. Our findings suggest the existence of dynamic, reversible site-specific methylation changes in response to cigarette smoking , which may contribute to the extended health risks associated with cigarette smoking.

In the sentence beginning 'The meQTL are enriched for functionally relevant characteristics', please provide a de inition for Hi-C.2. In the author list, there are two authors with the name Panos Deloukas with different af iliations. Please con irm whether this is the same individual and whether the two authors can be merged. 3.Please check your article carefully, coordinate with any co-authors and enter all inal edits clearly in the eproof, remembering to save frequently. Once corrections are submitted, we cannot routinely make further changes to the article. 4. Note that the eproof should be amended in only one browser window at any one time; otherwise changes will be overwritten. 5. Author surnames have been highlighted. Please check these carefully and adjust if the irst name or surname is marked up incorrectly. Note that changes here will affect indexing of your article in public repositories such as PubMed. Also, carefully check the spelling and numbering of all author names and af iliations, and the corresponding email address(es).6. You cannot alter accepted Supplementary Information iles except for critical changes to scienti ic content. If you do resupply any iles, please also provide a brief (but complete) list of changes. If these are not considered scienti ic changes, any altered Supplementary iles will not be used, only the originally accepted version will be published. If applicable, please ensure that any accession codes and datasets whoseDOIs or other identi iers are mentioned in the paper are scheduled for public release as soon as possible, we recommend within a few days of submitting your proof, and update the database record with publication details from this article once available. 8. Your paper has been copy edited. Please review every sentence to ensure that it conveys your intended meaning; if changes are required, please provide further clari ication rather than reverting to the original text. Please note that formatting (including hyphenation, Latin words,and any reference citations that might be mistaken for exponents) has been made consistent with our house style. 9. Please con irm or correct the city name inserted in af iliations 5, 10, 12, 13, 21, 22, 39 and 44. 10. In the sentence beginning 'Our meQTL replicate in data generated by the Illumina', please provide a de inition for EPIC and con irm whether the de inition provided for MeDIP-seq is correct. 11. In the legend for Fig. 2g, please provide a de inition for SAT. 12. In the sentence beginning 'We used summary data-based Mendelian', please con irm whether it is correct to say '0.05 ÷ n tests'. 13. In the sentence beginning 'We then tested the 4,811 eQTL genes', please provide a de inition for coloc PP4. 14. Please note, we reserve ‘signi icant' and its derivatives for statistical signi icance. This has been reworded where this was not the intended meaning (for example to important, notable, ...