Epigenome-wide association study of dietary fatty acid intake.

BACKGROUND: Dietary intake of n-3 polyunsaturated fatty acids (PUFA) may have a protective effect on the development of cardiovascular diseases, diabetes, depression and cancer, while a high intake of n-6 PUFA was often reported to be associated with inflammation-related traits. The effect of PUFAs on health outcomes might be mediated by DNA methylation (DNAm). The aim of our study is to identify the impact of PUFA intake on DNAm in the Cooperative Health Research in the Region of Augsburg (KORA) FF4 cohort and the Leiden Longevity Study (LLS). RESULTS: DNA methylation levels were measured in whole blood from the population-based KORA FF4 study (N = 1354) and LLS (N = 448), using the Illumina MethylationEPIC BeadChip and Illumina HumanMethylation450 array, respectively. We assessed associations between DNAm and intake of eight and four PUFAs in KORA and LLS, respectively. Where possible, results were meta-analyzed. Below the Bonferroni correction threshold (p < 7.17 × 10-8), we identified two differentially methylated positions (DMPs) associated with PUFA intake in the KORA study. The DMP cg19937480, annotated to gene PRDX1, was positively associated with docosahexaenoic acid (DHA) in model 1 (beta: 2.00 × 10-5, 95%CI: 1.28 × 10-5-2.73 × 10-5, P value: 6.98 × 10-8), while cg05041783, annotated to gene MARK2, was positively associated with docosapentaenoic acid (DPA) in our fully adjusted model (beta: 9.80 × 10-5, 95%CI: 6.25 × 10-5-1.33 × 10-4, P value: 6.75 × 10-8). In the meta-analysis, we identified the CpG site (cg15951061), annotated to gene CDCA7L below Bonferroni correction (1.23 × 10-7) associated with eicosapentaenoic acid (EPA) intake in model 1 (beta: 2.00 × 10-5, 95% CI: 1.27 × 10-5-2.73 × 10-5, P value = 5.99 × 10-8) and we confirmed the association of cg19937480 with DHA in both models 1 and 2 (beta: 2.07 × 10-5, 95% CI: 1.31 × 10-5-2.83 × 10-5, P value = 1.00 × 10-7 and beta: 2.19 × 10-5, 95% CI: 1.41 × 10-5-2.97 × 10-5, P value = 5.91 × 10-8 respectively). CONCLUSIONS: Our study identified three CpG sites associated with PUFA intake. The mechanisms of these sites remain largely unexplored, highlighting the novelty of our findings. Further research is essential to understand the links between CpG site methylation and PUFA outcomes.

A robust and standardized method to isolate and expand mesenchymal stromal cells from human umbilical cord.

BACKGROUND AIMS: Human umbilical cord-derived mesenchymal stromal cells (hUC-MSCs) are increasingly used in research and therapy. To obtain hUC-MSCs, a diversity of isolation and expansion methods are applied. Here, we report on a robust and standardized method for hUC-MSC isolation and expansion. METHODS: Using 90 hUC donors, we compared and optimized critical variables during each phase of the multi-step procedure involving UC collection, processing, MSC isolation, expansion and characterization. Furthermore, we assessed the effect of donor-to-donor variability regarding UC morphology and donor attributes on hUC-MSC characteristics. RESULTS: We demonstrated robustness of our method across 90 UC donors at each step of the procedure. With our method, UCs can be collected up to 6 h after birth, and UC-processing can be initiated up to 48 h after collection without impacting on hUC-MSC characteristics. The removal of blood vessels before explant cultures improved hUC-MSC purity. Expansion in Minimum essential medium α supplemented with human platelet lysate increased reproducibility of the expansion rate and MSC characteristics as compared with Dulbecco's Modified Eagle's Medium supplemented with fetal bovine serum. The isolated hUC-MSCs showed a purity of ∼98.9%, a viability of >97% and a high proliferative capacity. Trilineage differentiation capacity of hUC-MSCs was reduced as compared with bone marrow-derived MSCs. Functional assays indicated that the hUC-MSCs were able to inhibit T-cell proliferation demonstrating their immune-modulatory capacity. CONCLUSIONS: We present a robust and standardized method to isolate and expand hUC-MSCs, minimizing technical variability and thereby lay a foundation to advance reliability and comparability of results obtained from different donors and different studies.

Mesenchymal stromal cells as a tool to unravel the developmental origins of disease.

The intrauterine environment can induce alterations of the epigenome that have a lasting impact on disease risk. Current human studies in the field focus on a single epigenetic mark, DNA methylation, measured in blood. For in-depth mechanistic insight into the developmental origins of disease, it will be crucial to consider innovative tissue types. Mesenchymal stromal cells (MSCs) may serve as a novel tool to investigate the full epigenome beyond DNA methylation, to explore other omics levels, and to perform functional assays. Moreover, MSCs can be differentiated into multiple cell types and thereby mimic otherwise inaccessible cell types. A first wave of studies supports the potential of MSCs and illustrates how the innovative use of this cell type may be incorporated in birth cohorts.

Respiratory distress syndrome and bronchopulmonary dysplasia after fetal growth restriction: Lessons from a natural experiment in identical twins.

BACKGROUND: Fetal growth restriction (FGR) is thought to negatively affect lung development resulting in increased respiratory morbidity. However, research performed in singletons is often limited by a certain level of bias caused by individual differences in genetic constitution, obstetrical and maternal factors. METHODS: Respiratory morbidity was compared between the smaller and the larger twin in monochorionic twins with selective fetal growth restriction (sFGR), defined as a birth weight discordance ≥ 20%, born in our center between 2010 and 2019 in this retrospective study. Respiratory distress syndrome (RDS) was diagnosed based on the clinical picture of a neonate with respiratory failure requiring mechanical ventilation and/or surfactant, confirmed by a chest X-ray. Bronchopulmonary dysplasia (BPD) was diagnosed when the neonate required treatment with >21% oxygen for at least 28 days. FINDINGS: Median gestational age at birth for the 94 included pregnancies was 32.4 (IQR 30.4-34.3) weeks. Within-pair analyses showed that the prevalence of RDS was lower in the smaller twin compared to the larger twin, 19.1% (18/94) vs 34.0% (32/94), respectively (p = 0.004). The odds of RDS for the larger twin was doubled (OR 2.1 (CI95% 1.3-3.5). In contrast, the rate of BPD in the smaller twin was higher as opposed to the larger twin, 16.7% (15/90) vs 6.7% (6/89), respectively (p = 0.008), with a more than doubled odds (OR 2.5 (CI95% 1.3-4.9)). INTERPRETATION: Despite being genetically identical, sFGR twins have different respiratory outcomes. Adverse growth condition in utero in the smaller twin is associated with a reduced odds of RDS at birth but a more than doubled odds of BPD, reflecting the pathophysiologic adverse effect of growth restriction on lung development. FUNDING: The Dutch Heart Foundation (2017T075).

Genome-wide identification of genes regulating DNA methylation using genetic anchors for causal inference.

BACKGROUND: DNA methylation is a key epigenetic modification in human development and disease, yet there is limited understanding of its highly coordinated regulation. Here, we identify 818 genes that affect DNA methylation patterns in blood using large-scale population genomics data. RESULTS: By employing genetic instruments as causal anchors, we establish directed associations between gene expression and distant DNA methylation levels, while ensuring specificity of the associations by correcting for linkage disequilibrium and pleiotropy among neighboring genes. The identified genes are enriched for transcription factors, of which many consistently increased or decreased DNA methylation levels at multiple CpG sites. In addition, we show that a substantial number of transcription factors affected DNA methylation at their experimentally determined binding sites. We also observe genes encoding proteins with heterogenous functions that have widespread effects on DNA methylation, e.g., NFKBIE, CDCA7(L), and NLRC5, and for several examples, we suggest plausible mechanisms underlying their effect on DNA methylation. CONCLUSION: We report hundreds of genes that affect DNA methylation and provide key insights in the principles underlying epigenetic regulation.

Genome-wide analysis of constitutional DNA methylation in familial melanoma.

BACKGROUND: Heritable epigenetic alterations have been proposed as an explanation for familial clustering of melanoma. Here we performed genome-wide DNA methylation analysis on affected family members not carrying pathogenic variants in established melanoma susceptibility genes, compared with healthy volunteers. RESULTS: All melanoma susceptibility genes showed the absence of epimutations in familial melanoma patients, and no loss of imprinting was detected. Unbiased genome-wide DNA methylation analysis revealed significantly different levels of methylation in single CpG sites. The methylation level differences were small and did not affect reported tumour predisposition genes. CONCLUSION: Our results provide no support for heritable epimutations as a cause of familial melanoma.

A characterization of cis- and trans-heritability of RNA-Seq-based gene expression.

Insights into individual differences in gene expression and its heritability (h2) can help in understanding pathways from DNA to phenotype. We estimated the heritability of gene expression of 52,844 genes measured in whole blood in the largest twin RNA-Seq sample to date (1497 individuals including 459 monozygotic twin pairs and 150 dizygotic twin pairs) from classical twin modeling and identity-by-state-based approaches. We estimated for each gene h2total, composed of cis-heritability (h2cis, the variance explained by single nucleotide polymorphisms in the cis-window of the gene), and trans-heritability (h2res, the residual variance explained by all other genome-wide variants). Mean h2total was 0.26, which was significantly higher than heritability estimates earlier found in a microarray-based study using largely overlapping (>60%) RNA samples (mean h2 = 0.14, p = 6.15 × 10-258). Mean h2cis was 0.06 and strongly correlated with beta of the top cis expression quantitative loci (eQTL, ρ = 0.76, p < 10-308) and with estimates from earlier RNA-Seq-based studies. Mean h2res was 0.20 and correlated with the beta of the corresponding trans-eQTL (ρ = 0.04, p < 1.89 × 10-3) and was significantly higher for genes involved in cytokine-cytokine interactions (p = 4.22 × 10-15), many other immune system pathways, and genes identified in genome-wide association studies for various traits including behavioral disorders and cancer. This study provides a thorough characterization of cis- and trans-h2 estimates of gene expression, which is of value for interpretation of GWAS and gene expression studies.

TwinLIFE: The Twin Longitudinal Investigation of FEtal Discordance.

Lifelong health is thought to be partially set during intrauterine life by persistent epigenetic changes induced by the prenatal environment. To evaluate this hypothesis, we initiated a prospective longitudinal study in monochorionic (MC) twins: the TwinLIFE study. MC twins are monozygotic, thus in origin genetically identical, and share a single placenta. Although MC twins have many environmental factors in common, in one-third of the MC twin pairs, one fetus has significantly less access to nutrients and resources during pregnancy than its co-twin often resulting in a significant discordance in prenatal growth. Hence, MC twins constitute a unique natural experiment to study the influence of the prenatal environment on health. In TwinLIFE, we will chart intrapair differences in DNA methylation focusing on mesenchymal stromal cells isolated from cord as an advanced proxy of epigenetic dysregulation relevant for long-term health consequences. Next, we will follow up the MC twins for growth, cardiovascular and neurodevelopmental outcomes during childhood and evaluate the impact of an epigenetic signature at birth on future health. The current target is to include 100 MC twin pairs, but we aim to continue enrollment after procuring additional funding. TwinLIFE will not only address an unmet clinical need in the high-risk group of MC twins, but may also advance early-life strategies to prevent adverse growth, cardiovascular and neurodevelopmental outcomes in the general population.

Association of dietary folate and vitamin B-12 intake with genome-wide DNA methylation in blood: a large-scale epigenome-wide association analysis in 5841 individuals.

BACKGROUND: Folate and vitamin B-12 are essential micronutrients involved in the donation of methyl groups in cellular metabolism. However, associations between intake of these nutrients and genome-wide DNA methylation levels have not been studied comprehensively in humans. OBJECTIVE: The aim of this study was to assess whether folate and/or vitamin B-12 intake are asssociated with genome-wide changes in DNA methylation in leukocytes. METHODS: A large-scale epigenome-wide association study of folate and vitamin B-12 intake was performed on DNA from 5841 participants from 10 cohorts using Illumina 450k arrays. Folate and vitamin B-12 intakes were calculated from food-frequency questionnaires (FFQs). Continuous and categorical (low compared with high intake) linear regression mixed models were applied per cohort, controlling for confounders. A meta-analysis was performed to identify significant differentially methylated positions (DMPs) and regions (DMRs), and a pathway analysis was performed on the DMR annotated genes. RESULTS: The categorical model resulted in 6 DMPs, which are all negatively associated with folate intake, annotated to FAM64A, WRAP73, FRMD8, CUX1, and LCN8 genes, which have a role in cellular processes including centrosome localization, cell proliferation, and tumorigenesis. Regional analysis showed 74 folate-associated DMRs, of which 73 were negatively associated with folate intake. The most significant folate-associated DMR was a 400-base pair (bp) spanning region annotated to the LGALS3BP gene. In the categorical model, vitamin B-12 intake was associated with 29 DMRs annotated to 48 genes, of which the most significant was a 1100-bp spanning region annotated to the calcium-binding tyrosine phosphorylation-regulated gene (CABYR). Vitamin B-12 intake was not associated with DMPs. CONCLUSIONS: We identified novel epigenetic loci that are associated with folate and vitamin B-12 intake. Interestingly, we found a negative association between folate and DNA methylation. Replication of these methylation loci is necessary in future studies.

Human monocyte-to-macrophage differentiation involves highly localized gain and loss of DNA methylation at transcription factor binding sites.

BACKGROUND: Macrophages and their precursors monocytes play a key role in inflammation and chronic inflammatory disorders. Monocyte-to-macrophage differentiation and activation programs are accompanied by significant epigenetic remodeling where DNA methylation associates with cell identity. Here we show that DNA methylation changes characteristic for monocyte-to-macrophage differentiation occur at transcription factor binding sites, and, in contrast to what was previously described, are generally highly localized and encompass both losses and gains of DNA methylation. RESULTS: We compared genome-wide DNA methylation across 440,292 CpG sites between human monocytes, naïve macrophages and macrophages further activated toward a pro-inflammatory state (using LPS/IFNγ), an anti-inflammatory state (IL-4) or foam cells (oxLDL and acLDL). Moreover, we integrated these data with public whole-genome sequencing data on monocytes and macrophages to demarcate differentially methylated regions. Our analysis showed that differential DNA methylation was most pronounced during monocyte-to-macrophage differentiation, was typically restricted to single CpGs or very short regions, and co-localized with lineage-specific enhancers irrespective of whether it concerns gain or loss of methylation. Furthermore, differentially methylated CpGs were located at sites characterized by increased binding of transcription factors known to be involved in monocyte-to-macrophage differentiation including C/EBP and ETS for gain and AP-1 for loss of methylation. CONCLUSION: Our study highlights the involvement of subtle, yet highly localized remodeling of DNA methylation at regulatory regions in cell differentiation.

Skewed X-inactivation is common in the general female population.

X-inactivation is a well-established dosage compensation mechanism ensuring that X-chromosomal genes are expressed at comparable levels in males and females. Skewed X-inactivation is often explained by negative selection of one of the alleles. We demonstrate that imbalanced expression of the paternal and maternal X-chromosomes is common in the general population and that the random nature of the X-inactivation mechanism can be sufficient to explain the imbalance. To this end, we analyzed blood-derived RNA and whole-genome sequencing data from 79 female children and their parents from the Genome of the Netherlands project. We calculated the median ratio of the paternal over total counts at all X-chromosomal heterozygous single-nucleotide variants with coverage ≥10. We identified two individuals where the same X-chromosome was inactivated in all cells. Imbalanced expression of the two X-chromosomes (ratios ≤0.35 or ≥0.65) was observed in nearly 50% of the population. The empirically observed skewing is explained by a theoretical model where X-inactivation takes place in an embryonic stage in which eight cells give rise to the hematopoietic compartment. Genes escaping X-inactivation are expressed from both alleles and therefore demonstrate less skewing than inactivated genes. Using this characteristic, we identified three novel escapee genes (SSR4, REPS2, and SEPT6), but did not find support for many previously reported escapee genes in blood. Our collective data suggest that skewed X-inactivation is common in the general population. This may contribute to manifestation of symptoms in carriers of recessive X-linked disorders. We recommend that X-inactivation results should not be used lightly in the interpretation of X-linked variants.

Association of maternal prenatal smoking GFI1-locus and cardio-metabolic phenotypes in 18,212 adults.

BACKGROUND: DNA methylation at the GFI1-locus has been repeatedly associated with exposure to smoking from the foetal period onwards. We explored whether DNA methylation may be a mechanism that links exposure to maternal prenatal smoking with offspring's adult cardio-metabolic health. METHODS: We meta-analysed the association between DNA methylation at GFI1-locus with maternal prenatal smoking, adult own smoking, and cardio-metabolic phenotypes in 22 population-based studies from Europe, Australia, and USA (n = 18,212). DNA methylation at the GFI1-locus was measured in whole-blood. Multivariable regression models were fitted to examine its association with exposure to prenatal and own adult smoking. DNA methylation levels were analysed in relation to body mass index (BMI), waist circumference (WC), fasting glucose (FG), high-density lipoprotein cholesterol (HDL-C), triglycerides (TG), diastolic, and systolic blood pressure (BP). FINDINGS: Lower DNA methylation at three out of eight GFI1-CpGs was associated with exposure to maternal prenatal smoking, whereas, all eight CpGs were associated with adult own smoking. Lower DNA methylation at cg14179389, the strongest maternal prenatal smoking locus, was associated with increased WC and BP when adjusted for sex, age, and adult smoking with Bonferroni-corrected P < 0·012. In contrast, lower DNA methylation at cg09935388, the strongest adult own smoking locus, was associated with decreased BMI, WC, and BP (adjusted 1 × 10-7 < P < 0.01). Similarly, lower DNA methylation at cg12876356, cg18316974, cg09662411, and cg18146737 was associated with decreased BMI and WC (5 × 10-8 < P < 0.001). Lower DNA methylation at all the CpGs was consistently associated with higher TG levels. INTERPRETATION: Epigenetic changes at the GFI1 were linked to smoking exposure in-utero/in-adulthood and robustly associated with cardio-metabolic risk factors. FUND: European Union's Horizon 2020 research and innovation programme under grant agreement no. 633595 DynaHEALTH.

Genome-wide identification of directed gene networks using large-scale population genomics data.

Identification of causal drivers behind regulatory gene networks is crucial in understanding gene function. Here, we develop a method for the large-scale inference of gene-gene interactions in observational population genomics data that are both directed (using local genetic instruments as causal anchors, akin to Mendelian Randomization) and specific (by controlling for linkage disequilibrium and pleiotropy). Analysis of genotype and whole-blood RNA-sequencing data from 3072 individuals identified 49 genes as drivers of downstream transcriptional changes (Wald P < 7 × 10-10), among which transcription factors were overrepresented (Fisher's P = 3.3 × 10-7). Our analysis suggests new gene functions and targets, including for SENP7 (zinc-finger genes involved in retroviral repression) and BCL2A1 (target genes possibly involved in auditory dysfunction). Our work highlights the utility of population genomics data in deriving directed gene expression networks. A resource of trans-effects for all 6600 genes with a genetic instrument can be explored individually using a web-based browser.

IgG glycosylation and DNA methylation are interconnected with smoking.

BACKGROUND: Glycosylation is one of the most common post-translation modifications with large influences on protein structure and function. The effector function of immunoglobulin G (IgG) alters between pro- and anti-inflammatory, based on its glycosylation. IgG glycan synthesis is highly complex and dynamic. METHODS: With the use of two different analytical methods for assessing IgG glycosylation, we aim to elucidate the link between DNA methylation and glycosylation of IgG by means of epigenome-wide association studies. In total, 3000 individuals from 4 cohorts were analyzed. RESULTS: The overlap of the results from the two glycan measurement panels yielded DNA methylation of 7 CpG-sites on 5 genomic locations to be associated with IgG glycosylation: cg25189904 (chr.1, GNG12); cg05951221, cg21566642 and cg01940273 (chr.2, ALPPL2); cg05575921 (chr.5, AHRR); cg06126421 (6p21.33); and cg03636183 (chr.19, F2RL3). Mediation analyses with respect to smoking revealed that the effect of smoking on IgG glycosylation may be at least partially mediated via DNA methylation levels at these 7 CpG-sites. CONCLUSION: Our results suggest the presence of an indirect link between DNA methylation and IgG glycosylation that may in part capture environmental exposures. GENERAL SIGNIFICANCE: An epigenome-wide analysis conducted in four population-based cohorts revealed an association between DNA methylation and IgG glycosylation patterns. Presumably, DNA methylation mediates the effect of smoking on IgG glycosylation.

DNA methylation signatures of educational attainment.

Educational attainment is a key behavioural measure in studies of cognitive and physical health, and socioeconomic status. We measured DNA methylation at 410,746 CpGs (N = 4152) and identified 58 CpGs associated with educational attainment at loci characterized by pleiotropic functions shared with neuronal, immune and developmental processes. Associations overlapped with those for smoking behaviour, but remained after accounting for smoking at many CpGs: Effect sizes were on average 28% smaller and genome-wide significant at 11 CpGs after adjusting for smoking and were 62% smaller in never smokers. We examined sources and biological implications of education-related methylation differences, demonstrating correlations with maternal prenatal folate, smoking and air pollution signatures, and associations with gene expression in cis, dynamic methylation in foetal brain, and correlations between blood and brain. Our findings show that the methylome of lower-educated people resembles that of smokers beyond effects of their own smoking behaviour and shows traces of various other exposures.

Tea and coffee consumption in relation to DNA methylation in four European cohorts.

Lifestyle factors, such as food choices and exposure to chemicals, can alter DNA methylation and lead to changes in gene activity. Two such exposures with pharmacologically active components are coffee and tea consumption. Both coffee and tea have been suggested to play an important role in modulating disease-risk in humans by suppressing tumour progression, decreasing inflammation and influencing estrogen metabolism. These mechanisms may be mediated by changes in DNA methylation. To investigate if DNA methylation in blood is associated with coffee and tea consumption, we performed a genome-wide DNA methylation study for coffee and tea consumption in four European cohorts (N = 3,096). DNA methylation was measured from whole blood at 421,695 CpG sites distributed throughout the genome and analysed in men and women both separately and together in each cohort. Meta-analyses of the results and additional regional-level analyses were performed. After adjusting for multiple testing, the meta-analysis revealed that two individual CpG-sites, mapping to DNAJC16 and TTC17, were differentially methylated in relation to tea consumption in women. No individual sites were associated with men or with the sex-combined analysis for tea or coffee. The regional analysis revealed that 28 regions were differentially methylated in relation to tea consumption in women. These regions contained genes known to interact with estradiol metabolism and cancer. No significant regions were found in the sex-combined and male-only analysis for either tea or coffee consumption.

Controlling bias and inflation in epigenome- and transcriptome-wide association studies using the empirical null distribution.

We show that epigenome- and transcriptome-wide association studies (EWAS and TWAS) are prone to significant inflation and bias of test statistics, an unrecognized phenomenon introducing spurious findings if left unaddressed. Neither GWAS-based methodology nor state-of-the-art confounder adjustment methods completely remove bias and inflation. We propose a Bayesian method to control bias and inflation in EWAS and TWAS based on estimation of the empirical null distribution. Using simulations and real data, we demonstrate that our method maximizes power while properly controlling the false positive rate. We illustrate the utility of our method in large-scale EWAS and TWAS meta-analyses of age and smoking.

Age-related accrual of methylomic variability is linked to fundamental ageing mechanisms.

BACKGROUND: Epigenetic change is a hallmark of ageing but its link to ageing mechanisms in humans remains poorly understood. While DNA methylation at many CpG sites closely tracks chronological age, DNA methylation changes relevant to biological age are expected to gradually dissociate from chronological age, mirroring the increased heterogeneity in health status at older ages. RESULTS: Here, we report on the large-scale identification of 6366 age-related variably methylated positions (aVMPs) identified in 3295 whole blood DNA methylation profiles, 2044 of which have a matching RNA-seq gene expression profile. aVMPs are enriched at polycomb repressed regions and, accordingly, methylation at those positions is associated with the expression of genes encoding components of polycomb repressive complex 2 (PRC2) in trans. Further analysis revealed trans-associations for 1816 aVMPs with an additional 854 genes. These trans-associated aVMPs are characterized by either an age-related gain of methylation at CpG islands marked by PRC2 or a loss of methylation at enhancers. This distinct pattern extends to other tissues and multiple cancer types. Finally, genes associated with aVMPs in trans whose expression is variably upregulated with age (733 genes) play a key role in DNA repair and apoptosis, whereas downregulated aVMP-associated genes (121 genes) are mapped to defined pathways in cellular metabolism. CONCLUSIONS: Our results link age-related changes in DNA methylation to fundamental mechanisms that are thought to drive human ageing.

Genetic and environmental influences interact with age and sex in shaping the human methylome.

The methylome is subject to genetic and environmental effects. Their impact may depend on sex and age, resulting in sex- and age-related physiological variation and disease susceptibility. Here we estimate the total heritability of DNA methylation levels in whole blood and estimate the variance explained by common single nucleotide polymorphisms at 411,169 sites in 2,603 individuals from twin families, to establish a catalogue of between-individual variation in DNA methylation. Heritability estimates vary across the genome (mean=19%) and interaction analyses reveal thousands of sites with sex-specific heritability as well as sites where the environmental variance increases with age. Integration with previously published data illustrates the impact of genome and environment across the lifespan at methylation sites associated with metabolic traits, smoking and ageing. These findings demonstrate that our catalogue holds valuable information on locations in the genome where methylation variation between people may reflect disease-relevant environmental exposures or genetic variation.

Prenatal parental tobacco smoking, gene specific DNA methylation, and newborns size: the Generation R study.

BACKGROUND: Deleterious effects of prenatal tobacco smoking on fetal growth and newborn weight are well-established. One of the proposed mechanisms underlying this relationship is alterations in epigenetic programming. We selected 506 newborns from a population-based prospective birth cohort in the Netherlands. Prenatal parental tobacco smoking was assessed using self-reporting questionnaires. Information on birth outcomes was obtained from medical records. The deoxyribonucleic acid (DNA) methylation of the growth genes IGF2DMR and H19 was measured in newborn umbilical cord white blood cells. Associations were assessed between parental tobacco smoking and DNA methylation using linear mixed models and adjusted for potential confounders. RESULTS: The DNA methylation levels of IGF2DMR and H19 in the non-smoking group were median (90 % range), 54.0 % (44.6-62.0), and 30.0 % (25.5-34.0), in the first trimester only smoking group 52.2 % (44.5-61.1) and 30.8 % (27.1-34.1), and in the continued smoking group 51.6 % (43.9-61.3) and 30.2 % (23.7-34.8), respectively. Continued prenatal maternal smoking was inversely associated with IGF2DMR methylation (ß = -1.03, 95 % CI -1.76; -0.30) in a dose-dependent manner (P-trend = 0.030). This association seemed to be slightly more profound among newborn girls (ß = -1.38, 95 % CI -2.63; -0.14) than boys (ß = -0.72, 95 % CI -1.68; 0.24). H19 methylation was also inversely associated continued smoking <5 cigarettes/day (ß = -0.96, 95 % CI -1.78; -0.14). Moreover, the association between maternal smoking and newborns small for gestational age seems to be partially explained by IGF2DMR methylation (ß = -0.095, 95 % CI -0.249; -0.018). Among non-smoking mothers, paternal tobacco smoking was not associated with IGF2DMR or H19 methylation. CONCLUSIONS: Maternal smoking is inversely associated with IGF2DMR methylation in newborns, which can be one of the underlying mechanisms through which smoking affects fetal growth.