Heritable Methylation Research Articles

O03 Genetic impacts on the skin methylome

Abstract Introduction and aims DNA methylation at CpG dinucleotides is a dynamic epigenetic mechanism that has been linked to numerous skin phenotypes and diseases. Previous studies in blood show that a substantial proportion of the human DNA methylome exhibits evidence for genetic influence, and that 15–17% of narrow-sense methylation heritability is explained by genetic variants, or methylation quantitative trait loci (meQTLs). However, no genome-wide meQTL studies have previously been performed in skin. The primary purpose of this work was to identify the genetic basis of DNA methylation levels in whole skin. Methods Utilizing a twin-based heritability model, we estimated DNA methylation heritability using 414 bulk skin tissue samples from healthy older female twins from the TwinsUK cohort, profiled using the Illumina Infinium HumanMethylation450 BeadChip (Illumina, San Diego, CA, USA). We then integrated genomic data for 394 of these twins to perform genome-wide meQTL analyses. Finally, we integrated these meQTL Results with skin eQTL Results from an overlapping sample (N = 664). Results Over 8.2% of measured CpG sites showed evidence for strong heritability (h2 > 0.4), with a mean additive heritability estimate of 11.07% across all tested sites. We identified 76 442 CpG sites (18.8% of measured CpG sites) associated with a meQTL in cis (< 1 MB variant-CpG distance) and 1438 (0.35%) in trans. Both CpG sites under meQTL effects, and genetic variants underlying these meQTL effects in skin, are enriched in regions previously identified from genome-wide and epigenome-wide association studies for skin-related phenotypes, including psoriasis, melanoma and vitiligo. Additionally, some meQTLs influence gene expression via mediation of expression quantitative trait loci effects in skin, including in genes linked to melanoma and psoriasis. Conclusions Our study is the first genome-wide DNA methylation heritability and meQTL analysis of whole skin, giving insights into the regulatory landscape of epigenomic variation in skin.

Mixed Patterns of Intergenerational DNA Methylation Inheritance in Acropora.

For sessile organisms at high risk from climate change, phenotypic plasticity can be critical to rapid acclimation. Epigenetic markers like DNA methylation are hypothesized as mediators of plasticity; methylation is associated with the regulation of gene expression, can change in response to ecological cues, and is a proposed basis for the inheritance of acquired traits. Within reef-building corals, gene-body methylation (gbM) can change in response to ecological stressors. If coral DNA methylation is transmissible across generations, this could potentially facilitate rapid acclimation to environmental change. We investigated methylation heritability in Acropora, a stony reef-building coral. Two Acropora millepora and two Acropora selago adults were crossed, producing eight offspring crosses (four hybrid, two of each species). We used whole-genome bisulfite sequencing to identify methylated loci and allele-specific alignments to quantify per-locus inheritance. If methylation is heritable, differential methylation (DM) between the parents should equal DM between paired offspring alleles at a given locus. We found a mixture of heritable and nonheritable loci, with heritable portions ranging from 44% to 90% among crosses. gBM was more heritable than intergenic methylation, and most loci had a consistent degree of heritability between crosses (i.e. the deviation between parental and offspring DM were of similar magnitude and direction). Our results provide evidence that coral methylation can be inherited but that heritability is heterogenous throughout the genome. Future investigations into this heterogeneity and its phenotypic implications will be important to understanding the potential capability of intergenerational environmental acclimation in reef building corals.

18S rRNA methyltransferases DIMT1 and BUD23 drive intergenerational hormesis

Heritable non-genetic information can regulate a variety of complex phenotypes. However, what specific non-genetic cues are transmitted from parents to their descendants are poorly understood. Here, we perform metabolic methyl-labeling experiments to track the heritable transmission of methylation from ancestors to their descendants in the nematode Caenorhabditis elegans (C.elegans). We find heritable methylation in DNA, RNA, proteins, and lipids. We find that parental starvation elicits reduced fertility, increased heat stress resistance, and extended longevity in fed, naïve progeny. This intergenerational hormesis is accompanied by a heritable increase in N6'-dimethyl adenosine (m6,2A) on the 18S ribosomal RNA at adenosines 1735 and 1736. We identified DIMT-1/DIMT1 as the m6,2A and BUD-23/BUD23 as the m7G methyltransferases in C.elegans that are both required for intergenerational hormesis, while other rRNA methyltransferases are dispensable. This study labels and tracks heritable non-genetic material across generations and demonstrates the importance of rRNA methylation for regulating epigenetic inheritance.

Revisiting genetic artifacts on DNA methylation microarrays exposes novel biological implications

BackgroundIllumina DNA methylation microarrays enable epigenome-wide analysis vastly used for the discovery of novel DNA methylation variation in health and disease. However, the microarrays’ probe design cannot fully consider the vast human genetic diversity, leading to genetic artifacts. Distinguishing genuine from artifactual genetic influence is of particular relevance in the study of DNA methylation heritability and methylation quantitative trait loci. But despite its importance, current strategies to account for genetic artifacts are lagging due to a limited mechanistic understanding on how such artifacts operate.ResultsTo address this, we develop and benchmark UMtools, an R-package containing novel methods for the quantification and qualification of genetic artifacts based on fluorescence intensity signals. With our approach, we model and validate known SNPs/indels on a genetically controlled dataset of monozygotic twins, and we estimate minor allele frequency from DNA methylation data and empirically detect variants not included in dbSNP. Moreover, we identify examples where genetic artifacts interact with each other or with imprinting, X-inactivation, or tissue-specific regulation. Finally, we propose a novel strategy based on co-methylation that can discern between genetic artifacts and genuine genomic influence.ConclusionsWe provide an atlas to navigate through the huge diversity of genetic artifacts encountered on DNA methylation microarrays. Overall, our study sets the ground for a paradigm shift in the study of the genetic component of epigenetic variation in DNA methylation microarrays.

Ectopic targeting of CG DNA methylation in Arabidopsis with the bacterial SssI methyltransferase

The ability to target epigenetic marks like DNA methylation to specific loci is important in both basic research and in crop plant engineering. However, heritability of targeted DNA methylation, how it impacts gene expression, and which epigenetic features are required for proper establishment are mostly unknown. Here, we show that targeting the CG-specific methyltransferase M.SssI with an artificial zinc finger protein can establish heritable CG methylation and silencing of a targeted locus in Arabidopsis. In addition, we observe highly heritable widespread ectopic CG methylation mainly over euchromatic regions. This hypermethylation shows little effect on transcription while it triggers a mild but significant reduction in the accumulation of H2A.Z and H3K27me3. Moreover, ectopic methylation occurs preferentially at less open chromatin that lacks positive histone marks. These results outline general principles of the heritability and interaction of CG methylation with other epigenomic features that should help guide future efforts to engineer epigenomes.

VTRNA2-1: Genetic Variation, Heritable Methylation and Disease Association.

VTRNA2-1 is a metastable epiallele with accumulating evidence that methylation at this region is heritable, modifiable and associated with disease including risk and progression of cancer. This study investigated the influence of genetic variation and other factors such as age and adult lifestyle on blood DNA methylation in this region. We first sequenced the VTRNA2-1 gene region in multiple-case breast cancer families in which VTRNA2-1 methylation was identified as heritable and associated with breast cancer risk. Methylation quantitative trait loci (mQTL) were investigated using a prospective cohort study (4500 participants with genotyping and methylation data). The cis-mQTL analysis (334 variants ± 50 kb of the most heritable CpG site) identified 43 variants associated with VTRNA2-1 methylation (p < 1.5 × 10−4); however, these explained little of the methylation variation (R2 < 0.5% for each of these variants). No genetic variants elsewhere in the genome were found to strongly influence VTRNA2-1 methylation. SNP-based heritability estimates were consistent with the mQTL findings (h2 = 0, 95%CI: −0.14 to 0.14). We found no evidence that age, sex, country of birth, smoking, body mass index, alcohol consumption or diet influenced blood DNA methylation at VTRNA2-1. Genetic factors and adult lifestyle play a minimal role in explaining methylation variability at the heritable VTRNA2-1 cluster.

Tissue culture-induced DNA methylation in crop plants: a review.

Plant tissue culture techniques have been extensively employed in commercial micropropagation to provide year-round production. Tissue culture regenerants are not always genotypically and phenotypically similar. Due to the changes in the tissue culture microenvironment, plant cells are exposed to additional stress which induces genetic and epigenetic instabilities in the regenerants. These changes lead to tissue culture-induced variations (TCIV) which are also known as somaclonal variations to categorically specify the inducing environment. TCIV includes molecular and phenotypic changes persuaded in the in vitro culture due to continuous sub-culturing and tissue culture-derived stress. Epigenetic variations such as altered DNA methylation pattern are induced due to the above-mentioned factors. Reportedly, alteration in DNA methylation pattern is much more frequent in the plant genome during the tissue culture process. DNA methylation plays an important role in gene expression and regulation of plant development. Variants originated in tissue culture process due to heritable methylation changes, can contribute to intra-species phenotypic variation. Several molecular techniques are available to detect DNA methylation at different stages of in vitro culture. Here, we review the aspects of TCIV with respect to DNA methylation and its effect on crop improvement programs. It is anticipated that a precise and comprehensive knowledge of molecular basis of in vitro-derived DNA methylation will help to design strategies to overcome the bottlenecks of micropropagation system and maintain the clonal fidelity of the regenerants.

DNA methylomic homogeneity and heterogeneity in muscles and testes throughout pig adulthood.

DNA methylome pattern is significantly different among tissues, ages, breeds, and genders. We assessed 20 methylome and transcriptome data in longissimus dorsi (LD) or testicles from Bamaxiang (BMX) and Large White pigs (LW) by deep sequencing technology. We identified ~55.7M CpGs and 5.30M, 0.20M, 1.20M, and 0.16M differential CpGs (P<0.01) between tissues, ages, breeds, and genders, respectively. Interestingly, 7.54% of differentially methylated regions (DMRs) are co-localized with promoters, which potentially regulate gene expression. RNA-seq analysis revealed that 23.42% CpGs are significantly correlated with gene expression (mean |r|=0.58, P<0.01), most of which are enriched in tissue-specific functions. Specially, we also found that the methylation levels in promoters of 655 genes were strongly associated with their expression levels (mean |r|=0.66, P<0.01). In addition, differentially methylated CpGs (DMCpGs) between breeds in HOXC gene cluster imply important regulatory roles in myocytes hypertrophy and intermuscular fat (IMF) deposition. Dramatically, higher similarity of methylation pattern was observed within pedigree than across pedigrees, which indicates the existence of heritable methylation regions. In summary, a part of CpGs in promoter can change its methylation pattern and play a marked regulatory function in different physiological or natural environments.

Segregation of an MSH1 RNAi transgene produces heritable non-genetic memory in association with methylome reprogramming

MSH1 is a plant-specific protein. RNAi suppression of MSH1 results in phenotype variability for developmental and stress response pathways. Segregation of the RNAi transgene produces non-genetic msh1 ‘memory’ with multi-generational inheritance. First-generation memory versus non-memory comparison, and six-generation inheritance studies, identifies gene-associated, heritable methylation repatterning. Genome-wide methylome analysis integrated with RNAseq and network-based enrichment studies identifies altered circadian clock networks, and phytohormone and stress response pathways that intersect with circadian control. A total of 373 differentially methylated loci comprising these networks are sufficient to discriminate memory from nonmemory full sibs. Methylation inhibitor 5-azacytidine diminishes the differences between memory and wild type for growth, gene expression and methylation patterning. The msh1 reprogramming is dependent on functional HISTONE DEACETYLASE 6 and methyltransferase MET1, and transition to memory requires the RNA-directed DNA methylation pathway. This system of phenotypic plasticity may serve as a potent model for defining accelerated plant adaptation during environmental change.

DNA methylation QTL analysis identifies new regulators of human longevity.

Human longevity is a complex trait influenced by both genetic and environmental factors, whose interaction is mediated by epigenetic mechanisms like DNA methylation. Here, we generated genome-wide whole-blood methylome data from 267 individuals, of which 71 were long-lived (90–104 years), by applying reduced representation bisulfite sequencing. We followed a stringent two-stage analysis procedure using discovery and replication samples to detect differentially methylated sites (DMSs) between young and long-lived study participants. Additionally, we performed a DNA methylation quantitative trait loci analysis to identify DMSs that underlie the longevity phenotype. We combined the DMSs results with gene expression data as an indicator of functional relevance. This approach yielded 21 new candidate genes, the majority of which are involved in neurophysiological processes or cancer. Notably, two candidates (PVRL2, ERCC1) are located on chromosome 19q, in close proximity to the well-known longevity- and Alzheimer’s disease-associated loci APOE and TOMM40. We propose this region as a longevity hub, operating on both a genetic (APOE, TOMM40) and an epigenetic (PVRL2, ERCC1) level. We hypothesize that the heritable methylation and associated gene expression changes reported here are overall advantageous for the LLI and may prevent/postpone age-related diseases and facilitate survival into very old age.

Convergence of DNA Methylation Profiles of the Reef Coral Porites astreoides in a Novel Environment

Phenotypic acclimatization is an organismal response to environmental change that may be rooted in epigenetic mechanisms. In reef building corals, organisms that are severely threatened by environmental change, some evidence suggests that DNA methylation is an environmentally responsive mediator of acclimatization. We investigated changes in DNA methylation of the reef coral Porites astreoides in response to simulated environmental change. Coral colonies were sampled from a variety of habitats on the Belize Barrier Reef and transplanted to a common garden for one year. We used restriction site associated DNA sequencing, including a methylation-sensitive variant, to subsample the genome and assess changes in DNA methylation levels after a year in the common garden. Methylation changes among the 629 CpG loci we recovered were subtle, yet coral methylomes were more similar to each other after a year in the common garden together, indicating convergence of methylation profiles in the common environment. Differentially methylated loci showed matches with both coding and non-coding RNA sequences with putative roles in intracellular signaling, apoptosis, gene regulation, and epigenetic crosstalk. There was a positive and significant relationship between genetic and epigenetic variation, providing evidence of methylation heritability. Altogether, our results suggest that DNA methylation in P. astreoides is at least somewhat responsive to environmental change, reflective of the environment, and heritable, characteristics necessary for methylation to be implicated as part of potential transgenerational acclimatization responses.

Characterization of the contribution of shared environmental and genetic factors to metabolic syndrome methylation heritability and familial correlations

Characterization of the contribution of shared environmental and genetic factors to metabolic syndrome methylation heritability and familial correlations

A Bayesian mixed modeling approach for estimating heritability

BackgroundA Bayesian mixed model approach using integrated nested Laplace approximations (INLA) allows us to construct flexible models that can account for pedigree structure. Using these models, we estimate genome-wide patterns of DNA methylation heritability (h2), which are currently not well understood, as well as h2 of blood lipid measurements.MethodsWe included individuals from the Genetics of Lipid Lowering Drugs and Diet Network (GOLDN) study with Infinium 450 K cytosine-phosphate-guanine (CpG) methylation and blood lipid data pre- and posttreatment with fenofibrate in families with up to three-generation pedigrees. For genome-wide patterns, we constructed 1 model per CpG with methylation as the response variable, with a random effect to model kinship, and age and gender as fixed effects.ResultsIn total, 425,791 CpG sites pre-, but only 199,027 CpG sites posttreatment were found to have nonzero heritability. Across these CpG sites, the distributions of h2 estimates are similar in pre- and posttreatment (pre: median = 0.31, interquartile range [IQR] = 0.16; post: median = 0.34, IQR = 0.20). Blood lipid h2 estimates were similar pre- and posttreatment with overlapping 95% credibility intervals. Heritability was nonzero for treatment effect, that is, the difference between pre- and posttreatment blood lipids. Estimates for triglycerides h2 are 0.48 (pre), 0.42 (post), and 0.21 (difference); likewise for high-density lipoprotein cholesterol h2 the estimates are 0.61, 0.68, and 0.10.ConclusionsWe show that with INLA, a fully Bayesian approach to estimate DNA methylation h2 is possible on a genome-wide scale. This provides uncertainty assessment of the estimates, and allows us to perform model selection via deviance information criterion (DIC) to identify CpGs with strong evidence for nonzero heritability.

Characterization of the contribution of shared environmental and genetic factors to metabolic syndrome methylation heritability and familial correlations

BackgroundTransgenerational epigenetic inheritance has been posited as a possible contributor to the observed heritability of metabolic syndrome (MetS). Yet the extent to which estimates of epigenetic inheritance for DNA methylation sites are inflated by environmental and genetic covariance within families is still unclear. We applied current methods to quantify the environmental and genetic contributors to the observed heritability and familial correlations of four previously associated MetS methylation sites at three genes (CPT1A, SOCS3 and ABCG1) using real data made available through the GAW20.ResultsOur findings support the role of both shared environment and genetic variation in explaining the heritability of MetS and the four MetS cytosine-phosphate-guanine (CpG) sites, although the resulting heritability estimates were indistinguishable from one another. Familial correlations by type of relative pair generally followed our expectation based on relatedness, but in the case of sister and parent pairs we observed nonsignificant trends toward greater correlation than expected, as would be consistent with the role of shared environmental factors in the inflation of our estimated correlations.ConclusionsOur work provides an interesting and flexible statistical framework for testing models of epigenetic inheritance in the context of human family studies. Future work should endeavor to replicate our findings and advance these methods to more robustly describe epigenetic inheritance patterns in human populations.

Abstract 3314: Heritable methylation marks associated with prostate and breast cancer risk

Abstract While most epigenetic marks are reprogrammed during early embryogenesis, some studies have reported Mendelian-like inheritance of germline DNA methylation in particular in cancer susceptibility genes. For instance, individuals with MLH1 silenced throughout the soma fit the clinical criteria for hereditary nonpolyposis colorectal cancer that is indistinguishable from the syndrome resulting from germline mutations in MLH1. Research using multiple-case breast cancer families has shown that LINE-1 and Sat2 DNA methylation levels are lower in individuals with a strong family history. Family clustering of cancer could therefore be due to epigenetic as well as genetic and shared environmental factors. We have recently identified heritable methylation marks associated with breast and/or prostate cancer susceptibility by conducting a study involving 45 Australian multi-generational families with multiple cases of breast or prostate cancer who are not known to carry genetic mutations in cancer susceptibility genes. We developed and applied a new statistical method to identify heritable methylation marks based on complex segregation analysis and identified 24 and 41 methylation marks significantly associated with breast and prostate cancer risk respectively. Several marks across VTRNA2-1, a gene located in a differentially methylated region that is involved in imprinting and shows allele-specific methylation, were associated with heritable risk of both cancer types. A proportion of all identified marks were found to be associated with cancer risk in independent nested case-control studies (ie outside of the multiple-case family setting). We are expanding these successful studies to include additional families and to estimate HRs and age-specific cumulative risks of cancer associated with these marks to enable the incorporation of this information into clinical tools for risk prediction. Citation Format: Melissa C. Southey, Jihoon E. Joo, James G. Dowty, Roger L. Milne, EE Ming Wong, Pierre-Antoine Dugué, Dallas English, John L. Hopper, David E. Goldgar, Graham G. Giles. Heritable methylation marks associated with prostate and breast cancer risk [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 3314.

Heritable methylation marks associated with breast and prostate cancer risk.

DNA methylation can mimic the effects of germline mutations in cancer predisposition genes. Recently, we identified twenty-four heritable methylation marks associated with breast cancer risk. As breast and prostate cancer share genetic risk factors, including rare, high-risk mutations (eg, in BRCA2), we hypothesized that some of these heritable methylation marks might also be associated with the risk of prostate cancer. We studied 869 incident prostate cancers (430 aggressive and 439 non-aggressive) and 869 matched controls nested within a prospective cohort study. DNA methylation was measured in pre-diagnostic blood samples using the Illumina Infinium HM450K BeadChip. Conditional logistic regression models, adjusted for prostate cancer risk factors and blood cell composition, were used to estimate odds ratios and 95% confidence intervals for the association between the 24 methylation marks and the risk of prostate cancer. Five methylation marks within the VTRNA2-1 promoter region (cg06536614, cg00124993, cg26328633, cg25340688, and cg26896946), and one in the body of CLGN (cg22901919) were associated with the risk of prostate cancer. In stratified analyses, the five VTRNA2-1 marks were associated with the risk of aggressive prostate cancer. This work highlights a potentially important new area of investigation for prostate cancer susceptibility and adds to our knowledge about shared risk factors for breast and prostate cancer.

Heritable DNA methylation marks associated with susceptibility to breast cancer

Mendelian-like inheritance of germline DNA methylation in cancer susceptibility genes has been previously reported. We aimed to scan the genome for heritable methylation marks associated with breast cancer susceptibility by studying 25 Australian multiple-case breast cancer families. Here we report genome-wide DNA methylation measured in 210 peripheral blood DNA samples provided by family members using the Infinium HumanMethylation450. We develop and apply a new statistical method to identify heritable methylation marks based on complex segregation analysis. We estimate carrier probabilities for the 1000 most heritable methylation marks based on family structure, and we use Cox proportional hazards survival analysis to identify 24 methylation marks with corresponding carrier probabilities significantly associated with breast cancer. We replicate an association with breast cancer risk for four of the 24 marks using an independent nested case–control study. Here, we report a novel approach for identifying heritable DNA methylation marks associated with breast cancer risk.

Testing Two Evolutionary Theories of Human Aging with DNA Methylation Data

The evolutionary theories of mutation accumulation (MA) and disposable soma (DS) provide possible explanations for the existence of human aging. To better understand the relative importance of these theories, we devised a test to identify MA- and DS-consistent sites across the genome using familial DNA methylation data. Two key characteristics of DNA methylation allowed us to do so. First, DNA methylation exhibits distinct and widespread changes with age, with numerous age-differentially-methylated sites observed across the genome. Second, many sites show heritable DNA methylation patterns within families. We extended heritability predictions of MA and DS to DNA methylation, predicting that MA-consistent age-differentially-methylated sites will show increasing heritability with age, while DS-consistent sites will show the opposite. Variance components models were used to test for changing heritability of methylation with age at 48,601 age-differentially-methylated sites across the genome in 610 individuals from 176 families. Of these, 102 sites showed significant MA-consistent increases in heritability with age, while 2266 showed significant DS-consistent decreases in heritability. These results suggest that both MA and DS play a role in explaining aging and aging-related changes, and that while the majority of DNA methylation changes observed in aging are consistent with epigenetic drift, targeted changes exist and may mediate effects of aging-related genes.

Abstract 3357: Heritable methylation marks associated with breast cancer risk

Abstract While most epigenetic marks are reprogrammed during early embryogenesis, some studies have reported Mendelian-like inheritance of germline DNA methylation in particular cancer susceptibility genes. We aimed to identify heritable methylation marks associated with breast cancer susceptibility by studying 25 Australian multi-generational families with multiple cases of breast cancer who were not known to carry mutations in breast cancer susceptibility genes. Peripheral blood DNA methylation was measured at approximately 480,000 genetic loci (methylation sites) for 210 members of these families using the Infinium HumanMethylation450 BeadArray. We developed and applied a new statistical method to identify heritable methylation marks based on complex segregation analysis. Carrier probabilities were estimated for the 1000 most heritable methylation marks, based on family structure but not disease status, and Cox proportional hazards survival analysis was used to test for associations between these carrier probabilities and breast cancer. After correcting for multiple testing, we identified 11 methylation marks whose corresponding carriers probabilities were significantly associated with breast cancer. Three marks clustered within 200 base pairs of the noncoding RNA miR886. We found evidence for these marks to be associated with aggressive (but not non-aggressive) prostate cancer risk in a population-based prostate cancer study. The association with breast cancer risk was replicated for three of the eleven marks in an independent population-based breast cancer study. Our work demonstrates the efficacy of our novel approach for identifying heritable methylation marks associated with breast cancer. This methodology may provide insights into other disease settings where known lifestyle exposures and genetic risk variants do not fully explain the familial relative risk. Citation Format: Melissa C. Southey, Jihoon E. Joo, James G. Dowty, Roger L. Milne, EE Ming Wong, Pierre-Antoine Dugué, Dallas English, John L. Hopper, Daivd E. Goldgar, Graham G. Giles. Heritable methylation marks associated with breast cancer risk [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 3357. doi:10.1158/1538-7445.AM2017-3357

DNA methylation changes at infertility genes in newborn twins conceived by in vitro fertilisation

BackgroundThe association of in vitro fertilisation (IVF) and DNA methylation has been studied predominantly at regulatory regions of imprinted genes and at just thousands of the ~28 million CpG sites in the human genome.MethodsWe investigated the links between IVF and DNA methylation patterns in whole cord blood cells (n = 98) and cord blood mononuclear cells (n = 82) from newborn twins using genome-wide methylated DNA immunoprecipitation coupled with deep sequencing.ResultsAt a false discovery rate (FDR) of 5%, we identified one significant whole blood DNA methylation change linked to conception via IVF, which was located ~3 kb upstream of TNP1, a gene previously linked to male infertility. The 46 most strongly associated signals (FDR of 25%) included a second region in a gene also previously linked to infertility, C9orf3, suggesting that our findings may in part capture the effect of parental subfertility. Using twin modelling, we observed that individual-specific environmental factors appear to be the main overall contributors of methylation variability at the FDR 25% IVF-associated differentially methylated regions, although evidence for methylation heritability was also obtained at several of these regions. We replicated previous findings of differential methylation associated with IVF at the H19/IGF2 region in cord blood mononuclear cells, and we validated the signal at C9orf3 in monozygotic twins. We also explored the impact of intracytoplasmic sperm injection on the FDR 25% signals for potential effects specific to male or female infertility factors.ConclusionsTo our knowledge, this is the most comprehensive study of DNA methylation profiles at birth and IVF conception to date, and our results show evidence for epigenetic modifications that may in part reflect parental subfertility.