Environmentally induced DNA methylation is inherited across generations in an aquatic keystone species

Transgenerational inheritance of environmentally induced epigenetic marks can have significant impacts on eco-evolutionary dynamics, but the phenomenon remains controversial in ecological model systems. We used whole-genome bisulfite sequencing of individual water fleas (Daphnia magna) to assess whether environmentally induced DNA methylation is transgenerationally inherited. Genetically identical females were exposed to one of three natural stressors, or a de-methylating drug, and their offspring were propagated clonally for four generations under control conditions. We identified between 70 and 225 differentially methylated CpG positions (DMPs) in F1 individuals whose mothers were exposed to a natural stressor. Roughly half of these environmentally induced DMPs persisted until generation F4. In contrast, treatment with the drug demonstrated that pervasive hypomethylation upon exposure is reset almost completely after one generation. These results suggest that environmentally induced DNA methylation is non-random and stably inherited across generations in Daphnia, making epigenetic inheritance a putative factor in the eco-evolutionary dynamics of freshwater communities

The tissue-specific aspect of genome-wide DNA methylation in newborn and placental tissues: Implications for epigenetic epidemiologic studies

Epigenetic programming is essential for lineage differentiation, embryogenesis and placentation in early pregnancy. In epigenetic association studies, DNA methylation is often examined in DNA derived from white blood cells, although its validity to other tissues of interest remains questionable. Therefore, we investigated the tissue specificity of epigenome-wide DNA methylation in newborn and placental tissues. Umbilical cord white blood cells (UC-WBC, n = 25), umbilical cord blood mononuclear cells (UC-MNC, n = 10), human umbilical vein endothelial cells (HUVEC, n = 25) and placental tissue (n = 25) were obtained from 36 uncomplicated pregnancies. Genome-wide DNA methylation was measured by the Illumina HumanMethylation450K BeadChip. Using UC-WBC as a reference tissue, we identified 3595 HUVEC tissue-specific differentially methylated regions (tDMRs) and 11,938 placental tDMRs. Functional enrichment analysis showed that HUVEC and placental tDMRs were involved in embryogenesis, vascular development and regulation of gene expression. No tDMRs were identified in UC-MNC. In conclusion, the extensive amount of genome-wide HUVEC and placental tDMRs underlines the relevance of tissue-specific approaches in future epigenetic association studies, or the use of validated representative tissues for a certain disease of interest, if available. To this purpose, we herewith provide a relevant dataset of paired, tissue-specific, genome-wide methylation measurements in newborn tissues

DNA methylation differences at birth after conception through ART

STUDY QUESTION: Is there a relation between ART and DNA methylation (DNAm) patterns in cord blood, including any differences between IVF and ICSI? SUMMARY ANSWER: DNAm at 19 CpGs was associated with conception via ART, with no difference found between IVF and ICSI. WHAT IS KNOWN ALREADY: Prior studies on either IVF or ICSI show conflicting outcomes, as both widespread effects on DNAm and highly localized associations have been reported. No study on both IVF and ICSI and genome-wide neonatal DNAm has been performed. STUDY DESIGN, SIZE, DURATION: This was a cross-sectional study comprising 87 infants conceived with IVF or ICSI and 70 conceived following medically unassisted conception. The requirement for inclusion in the study was an understanding of the Swedish language and exclusion was the use of donor gametes. PARTICIPANTS/MATERIALS, SETTING, METHODS: Participants were from the UppstART study, which was recruited from fertility and reproductive health clinics, and the Born into Life cohort, which is recruited from the larger LifeGene study. We measured DNAm from DNA extracted from cord blood collected at birth using a micro-array (450k array). Group differences in DNAm at individual CpG dinucleotides (CpGs) were determined using robust linear models and post-hoc Tukey's tests. MAIN RESULTS AND THE ROLE OF CHANCE: We found no association of ART conception with global methylation levels, imprinted loci and meta-stable epialleles. In contrast, we identify 19 CpGs at which DNAm was associated with being conceived via ART (effect estimates: 0.5-4.9%, PFDR < 0.05), but no difference was found between IVF and ICSI. The associated CpGs map to genes related to brain function/development or genes connected to the plethora of conditions linked to subfertility, but funct

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

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

Blood lipids influence DNA methylation in circulating cells

Background: Cells can be primed by external stimuli to obtain a long-term epigenetic memory. We hypothesize that long-term exposure to elevated blood lipids can prime circulating immune cells through changes in DNA methylation, a process that may contribute to the development of atherosclerosis. To interrogate the causal relationship between triglyceride, low-density lipoprotein (LDL) cholesterol, and high-density lipoprotein (HDL) cholesterol levels and genome-wide DNA methylation while excluding confounding and pleiotropy, we perform a stepwise Mendelian randomization analysis in whole blood of 3296 individuals. Results: This analysis shows that differential methylation is the consequence of inter-individual variation in blood lipid levels and not vice versa. Specifically, we observe an effect of triglycerides on DNA methylation at three CpGs, of LDL cholesterol at one CpG, and of HDL cholesterol at two CpGs using multivariable Mendelian randomization. Using RNA-seq data available for a large subset of individuals (N = 2044), DNA methylation of these six CpGs is associated with the expression of CPT1A and SREBF1 (for triglycerides), DHCR24 (for LDL cholesterol) and

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

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

Autosomal genetic variation is associated with DNA methylation in regions variably escaping X-chromosome inactivation

X-chromosome inactivation (XCI), i.e., the inactivation of one of the female X chromosomes, restores equal expression of X-chromosomal genes between females and males. However, ~10% of genes show variable degrees of escape from XCI between females, although little is known about the causes of variable XCI. Using a discovery data-set of 1867 females and 1398 males and a replication sample of 3351 females, we show that genetic variation at three autosomal loci is associated with female-specific changes in X-chromosome methylation. Through cis-eQTL expression analysis, we map these loci to the genes SMCHD1/METTL4, TRIM6/HBG2, and ZSCAN9. Low-expression alleles of the loci are predominantly associated with mild hypomethylation of CpG islands near genes known to variably escape XCI, implicating the autosomal genes in variable XCI. Together, these results suggest a genetic basis for variable escape from XCI and highlight the potential of a population genomics approach to identify genes involved in XCI

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