Genome-wide DNA methylation profiling with MeDIP-seq using archived dried blood spots
© The Author(s). 2016
Received: 4 May 2016
Accepted: 27 June 2016
Published: 26 July 2016
In utero and early-life experienced environmental exposures are suggested to play an important role in many multifactorial diseases potentially mediated through lasting effects on the epigenome. As the epigenome in addition remains modifiable throughout life, identifying specific disease-relevant biomarkers may prove challenging. This has led to an increased interest in epigenome-wide association studies using dried blood spots (DBS) routinely collected in perinatal screening programs. Such programs are in place in numerous countries around the world producing large and unique biobanks. However, availability of this biological material is highly limited as each DBS is made only from a few droplets of blood and storage conditions may be suboptimal for epigenetic studies. Furthermore, as relevant markers may reside outside gene bodies, epigenome-wide interrogation is needed.
Here we demonstrate, as a proof of principle, that genome-wide interrogation of the methylome based on methylated DNA immunoprecipitation coupled with next-generation sequencing (MeDIP-seq) is feasible using a single 3.2 mm DBS punch (60 ng DNA) from filter cards archived for up to 16 years. The enrichment profile, sequence quality and distribution of reads across genetic regions were comparable between samples archived 16 years, 4 years and a freshly prepared control sample.
In summary, we show that high-quality MeDIP-seq data is achievable from neonatal screening filter cards stored at room temperature, thereby providing information on annotated as well as on non-RefSeq genes and repetitive elements. Moreover, the quantity of DNA from one DBS punch proved sufficient allowing for multiple epigenome studies using one single DBS.
Epigenetic regulation has proven important in numerous cellular mechanisms including cell differentiation, gene expression and genome stability [1, 2]. With mounting evidence of epigenetic variations playing an important role in the etiology of common complex diseases such as cancer, diabetes and psychiatric disorders, the interest in epigenome-wide association studies (EWAS), specifically with regard to DNA methylation, has increased extensively in recent years [3–6]. However, identification of epigenetic risk variants with small effect sizes in a heterogenic population requires sample sets of considerable sizes. Additionally, given the modifiable nature of the epigenome distinguishing causality from consequence can prove challenging in adult samples, as marks dynamically change according to cellular environment and differentiation but also as a consequence of adverse exposures and disease. In fact, exposures experienced prenatally have shown to be able to inflict, even long-lasting, epigenetic changes. For example, maternal smoking or stress during pregnancy has shown to result in aberrant epigenetic changes in the offspring [7–9]. In addition, the inherent plasticity of the epigenome leads to changes not only over time per se but also in response to the use of certain pharmaceuticals and other potentially “epitoxic” substances. For instance, DNA methylation changes have been associated with the use of anti-psychotic drugs such as clozapine and haloperidol [10, 11]. A variety of common diseases including type II diabetes and multiple sclerosis have been suggested to include a component of in utero origin [12, 13]. The epigenome is especially vulnerable at this stage, laying the foundation for the “developmental origins of health and disease hypothesis”, suggesting that in utero experienced exposures affect the epigenome . Interrogation of the epigenetic profile is, therefore, preferably done before disease onset and ideally perinatally.
Screening for metabolic disorders in newborns is routinely undertaken in many countries, where blood from heel-pricks is collected on filter paper cards, also known as Guthrie cards. Subsequent storage of the cards essentially produces a population-wide biobank. Notable, since prenatal exposures are believed to have impact on multiple cell lineages, mirror effects may arise. Hence, methylome profiling in DNA extracted from dried blood spot (DBS) samples can potentially be informative both in the search for biomarkers and to delineate the etiology of diseases manifested in other tissues [15, 16]. Also, whereas the intra-individual variation in the epigenetic landscape is considerable, the inter-individual variation is relatively small, making not only longitudinal but also cross-individual studies feasible .
Nested case-control approaches taking advantage of large biobanks containing DBS filter cards can potentially satisfy both abovementioned needs of large and early-life sample sets provided that the epigenome remains stable under the storage conditions of the particular biobank and that sufficient biological material can be obtained. Previous studies have shown that DNA can be extracted from DBSs stored at −20°C for decades without considerable loss of quality paving the way for genome-wide association studies (GWAS) and whole genome sequencing (WGS) [18, 19]. Along the same line, both the use of methylation arrays and methylome-wide sequencing using DNA extracted from DBSs has previously shown to be feasible, however, requiring whole genome amplification of bisulfite converted DNA or large input amounts [20–24]. However, amplification of the bisulfitome (bisulfite-converted genomic DNA) has proven to introduce biases [25, 26]. Methylated DNA immunoprecipitation coupled with next-generation sequencing (MeDIP-seq) is a cost-effective methylome alternative, where methylated genomic DNA is immunoprecipitated with an antibody followed by sequencing of the enriched fragments. MeDIP-seq displays good positive correlation with array-based methods such as the HumanMethylation450 BeadChip [21, 27] and also provides information on non-RefSeq genes and repetitive elements, both of which may be of importance in disease development and not interrogated with current methylation arrays [28, 29]. The method of choice is purpose dependent. Arrays may be the appropriate choice for disease studies searching for common variants, whereas MeDIP-seq provides agnostic information on the methylation profile and can be used for a broader characterization of cell development and the impact of environmental exposure outside annotated genomic regions. The minimum DNA input requirements by the different methods is another limiting factor in the case of biobanked DBS filter cards (Additional file 1: Figure S1). With 50 ng of DNA sufficient for MeDIP-seq, this method is one of the least consuming whole genome approaches.
Storage conditions vary among the biobanks, thus, where the Danish National Screening Biobank (DNSB), the California Research-Ready Biospecimen Bank store the filter cards at around −20°C, the Swedish PKU-Biobank, the Michigan Neonatal Biobank and the Danish Twin Registry store the majority of their samples at room temperature.
In this study, we profiled the methylome of DNA stored on room temperature archived filter cards by the use MeDIP-seq using a single 3.2 mm punch. A freshly prepared filter card served as benchmark control. Despite prolonged (4–16 years) storage at room temperature the methylome on a global level appeared largely unchanged. These results underline that DNA extracted from only a small punch of an invaluable archival DBS on a filter card, can be used for genome-wide methylome analyses including methylome-wide association studies (MWAS).
MeDIP sequencing quality from DBS samples
hDBS (“homemade dried blood spot”)
1 month −20 °C
rDBS (“recent dried blood spot”)
4 years RT
oDBS (“old dried blood spot”)
16 years RT
Correlation between the three DBS samples revealed no global changes in methylation patterns
We looked in greater detail at the distribution of reads along genomic segments defined as 500 bp sliding windows with an overlap of 250 bp. The correlation was normalized only to the total number of reads and log transformed. Cumulating the number of reads from the lowest covered percentile of the genome to the highest should be largely identical between sets of the same type despite defined regional differences such as differently methylated regions (DMRs). The cumulative paths for all three sets appeared much alike, suggesting that no major technically or biologically introduced differences were present, such as large uncovered regions (Fig. 2b). The inconsistency at the lower end of the plot can likely be attributed to the low number of absolute counts of reads. Depicting the read density along the chromosomes in a domainogram again indicated that no systematic bias was evident comparing the three DBS samples. That is, the regional density of reads per segment indicated by a continuum from green (read deserts) to red is overall identical (Fig. 2c). As a check for biological validity, the X chromosome in rDBS (female) had a higher read density as expected, due to X chromosome inactivation, compared to hDBS and oDBS (both males). In addition, no coverage on the Y chromosome in the rDBS sample was observed, as expected.
Canonical DNA methylation patterns at defined genomic regions
About one-third of the genomic DNA methylation occurs in repetitive elements such as retrotransposons . Thus, analysis of the methylation level in such elements can be viewed as a proxy for the global genomic methylation status. Bisulfite pyrosequencing of three promoter CpG sites within the highly repetitive long interspersed element 1 (LINE-1) was performed on all three samples (Additional file 9: Figure S8). The LINE-1 CpG sites are normally heavily methylated to prevent retrotransposition. In all our DSB samples, all three sites were highly methylated and within the range of previously reported LINE-1 methylation levels in blood, approximately 80 % [32, 37].
No apparent systematic change in methylation observed between the DBS samples when clustering based on significant different segments
An unassisted hierarchical clustering of the significantly different segments of each genomic region did not show a clear systematic grouping. There was no general tendency of reduced read density in the archived DBSs compared to the control. Subdividing the clusters, based on a threshold for correlation set to 0.7, divided the cluster into five-six groups. In none of the groups did the segments aggregate at specific genomic locations.
Only few differentially methylated windows are common for the archived DBS samples
Retrospective methylome studies based on DBSs from biobanked filter cards hold great promise but are potentially hampered by poorer DNA quality and limited availability. In this study, we show that prolonged storage of DBSs, even at room temperature, has little overall influence on the methylome and that methylome-wide coverage can be obtained using DNA extracted from a single 3.2 mm punch, allowing for multiple studies on the same DBS. It has previously been shown that 50 ng of input DNA is sufficient for generating methylation libraries with adequate complexity using MeDIP-seq . A further decrease to 1 ng input is possible, but at the expense of an increased rate of duplicates and immunoprecipitation bias, which leads to loss of data and potentially false-positive findings [22, 39].
From single 3.2 mm DBS punches, we consistently extracted good quality DNA and with no quantitatively loss following prolonged storage, emphasizing the robustness of the extraction methodology applied. In the agnostic search for epigenetic markers or risk variants, a genome-wide approach such as MeDIP may prove valuable. Importantly, DNA extracted in this work proved adequate for MeDIP-seq. The genomic distribution of reads and coverage showed a high positive correlation between the three DBS samples. Moreover, the expected sharp and canonical drop in methylation level around the transcription start sites, as well as the gradual increase in methylation throughout the gene body was evident in all samples [40, 41].
A limitation of the study is that the compared filter cards were not from the same individual. The samples could also not be matched for sex or age [42, 43]. A follow-up study with more and better-matched samples is, therefore, warranted. Despite this, we were able to show that the overall methylation pattern is maintained in DNA isolated from old filter cards stored at room temperature and that this good quality of data can be obtained using MeDIP-seq on low amounts of DNA input. The analysis showed that only very few regions (<0.05 %) differed significantly between at least any two of the three DBS samples. Suggesting that the level of technical or storage-introduced noise is limited. Furthermore, although the three DBS samples originate from different individuals, and thus is not a longitudinal study, the comparison did not reveal any systematic change over the time stored.
Combining the archived DBS samples into one sample set and comparing this to the freshly prepared DBS defined only 149 significant 100 bp windows indicating that even without removal of confounding effects the difference appears limited. The majority of the significant windows were located in repetitive regions. This may be attributed not only to high biological variance at these sites, especially in pericentromeric satellites and subtelomeric elements, but also to alignment issues [44–46]. A previous inter-individual comparison based on MeDIP-seq and peripheral blood monocytes found differences to be highly enriched in repetitive elements, and especially in satellites . Pyrosequencing of three CpGs in the retrotransposable element LINE-1, a surrogate marker for global wide DNA methylation status, showed no apparent difference between the three DBS samples. This may again indicate that an overall loss of methylation has not taken place in the archived DBS samples.
It is important to note that defining true environment-specific or inter-individual DMRs requires that cell type composition and the genetic profile are taken into account. Means of cell composition correction in silico has been devised for array-generated data based on cell type specific methylation patterns [47, 48]. On the other hand, incorporation of genotyping data permits exploration of methylation quantitative trait loci (mQTLs), being a measure of the relationship between genetic polymorphisms and methylation . Notably, a considerable overlap of mQTLs across-tissue has been documented [15, 49]. Also, this allows for examination of gene and environment interactions and the effects on the epigenome .
In conclusion, we demonstrated that the methylome profile is highly comparable between DNA from archived DBS samples and DNA from a freshly prepared DBS sample. Moreover, we show, to our knowledge, for the first time that MeDIP-seq can be performed on neonatal screening cards even with small input amounts. Together, our findings add to the notion that archival DBS samples can be used for sequencing-based epigenome studies, thereby holding great potential in epidemiological research.
hDBS: Approximately 70 μL freshly drawn peripheral blood was spotted on Whatman 903 Protein Saver cards (GE Healthcare Life Sciences, Piscataway, NJ, USA). The cards were dried and stored at −20°C.
rDBS/oDBS: Randomly selected from anonymous DBS samples stored at room temperature for 4 years (rDBS) and 16 years (oDBS), respectively. For both, peripheral blood was spotted on the Whatman 903 Protein Saver cards by means of a finger prick.
Discs were punched in triplicates from the DBSs by the use of a puncher. Samples stored at RT were incubated in 180 μl PBS at 4 °C overnight in order to increase DNA extraction yield. Following the overnight incubation, PBS was removed and all samples were subjected to the DNA extraction protocol. DNA extraction was performed according to a previously published protocol , with introduction of additional centrifugation and ethanol washing steps. Concentration of the extracted DNA was measured using the Qubit instrument and the Qubit dsDNA HS Assay Kit (Thermo Fisher Scientific, Cleveland, OH, USA).
Extracted genomic DNA (gDNA) was sonicated in a Pico Bioruptor (Diagenode, Seraing, Belgium) at 10 ng/μL for 13 cycles of 30 s on 30 s off to a mean fragment size of 180 bp. Fragment length distribution was assessed by microelectrophoresis using the Qiaxcel instrument and a high-resolution gel cartridge (Qiagen, Hilden, Germany). Sonicated gDNA (in the range of 93–111 ng) was further used for end-repair and adaptor ligation employing the NEBNext Ultra DNA Library Prep Kit for Illumina (New England Biolabs, MA, USA). The reaction mix was purified using Ampure XP beads (Beckman-Coulter, CA, USA) after which methylated DNA immunoprecipitation (MeDIP) was set up on the SX-8G-IP-Star Compact robot (Diagenode) according to manufacturer’s instructions applying the Auto MeDIP kit (Diagenode) and including unmethylated and methylated spike-in controls. Effectively between 60 and 70 ng of adaptor-ligated DNA was used in the MeDIP process. Antibody incubation was performed at 4 °C for 15 h. Immunoprecipitated samples were magnetically purified on the robot using the Auto iPure v2 kit (Diagenode) according to the manufacturer’s instructions. Recovery and enrichment was evaluated by qPCR using primer sets specific for the spike-in controls. Minimum criteria were set to 10 % recovery and 25-fold enrichment. Based on the recovery rate, samples were PCR-amplified at 13 cycles (oDBS) or 14 cycles (hDBS and rDBS) using multiplex oligos (New England Biolabs, MA, USA). Samples were size-selected with Ampure XP beads on the SX-8G-IP-Star Compact robot (Diagenode) using a total of 90 μL beads per sample and eluted in 25 μL DNase-free H2O. Purity and fragment length distribution was evaluated on the Qiaxcel instrument. Post-amplification enrichment was verified by qPCR using primer pairs targeting the endogenous hypermethylated promoter region of testis specific histone 2B (TSH2B) (Diagenode, cat.nr. C17011041) or the hypomethylated TSS of glyceraldehyde 3-phosphate dehydrogenase (GAPDH) (Diagenode, cat.nr. C17011047).
Sequencing and bioinformatics
Samples were PE50 sequenced on a single lane on a HiSeq2000 instrument (Illumina, CA, USA), generating around 64–71 million clean reads per sample. The reads were cleaned using Soapnuke 1.5.0 tool (BGI, Shenzhen, China) applying the following filters: remove reads with adaptors, remove reads with >10 % unknown reads and remove reads with >50 % low quality bases (Q score <5). Clean reads only qualify if Q20 ≥ 85 %.
Using the Galaxy platform, clean reads were groomed and aligned to the human genome (build hg38) using BWA. SAM files were converted to BAM files and filtered to include only de-duplicated uniquely mapped reads resulting in ~54–58 million reads per sample. BAM files were imported to R and SeqMonk v0.32.1 (Barbraham Institute). The R package MEDIPS was used for quality control, genomic coverage estimation and differential coverage analysis . In SeqMonk segments were generated by quantitation of read counts in genetic windows corrected for total counts and with duplicates ignored. Filters include subdivision based on defined features (e.g. genes and CGIs) on any strand and calculating significantly different segments in pair-wise comparisons with a minimum intensity difference p value threshold of 0.05 with multiple testing correction (Bonferroni) applied. UCSC Genome Browser was used for visualization and interpretation of genomic regions.
Genomic DNA from two 3.2 mm DBS punches were pooled and bisulfite-treated using the EpiTect Bisulfite Kit (Qiagen) following the manufacturer’s protocol. A 146 bp fragment of the LINE-1 promoter was amplified by PCR from 100 ng bisulfite-treated DNA and preprocessed for sequencing using the PyroMark Q24 LINE-1 kit (Qiagen) following the manufacturer’s instructions. Prepared PCR products were sequenced on a PyroMark Q24 Advanced and analyzed on the appertaining software (Qiagen) according to the manufacturer’s instructions.
BWA, Burrows-Wheeler aligner; CGI, CpG Island; DBS, dried blood spots; DMR, differentially methylated region; DNSB, Danish Neonatal Screening Biobank; EWAS, epigenome-wide association study; GWAS, genome-wide association study; LINE-1, long interspersed nuclear element 1 ; MeDIP, methylated DNA immunoprecipitation; mQTLs, methylation quantitative trait loci; MWAS, methylome-wide association study; TSS, transcription start site; WGS, whole genome sequencing
We would like to thank Marit Nielsen and Tina Fuglsang Daugaard for the expert technical assistance at Translational Neuropsychiatric Unit, Aarhus University Hospital and Department of Biomedicine, University of Aarhus, respectively.
This work was funded by the Toyota Foundation and The Lundbeck Foundation, Denmark.
Availability of data and materials
The MeDIP-seq dataset supporting the conclusions of this article is available upon request.
NHS, AB and OM conceived and designed the experiments. LC collected and supplied samples. NHS and AS performed the experiments. NHS, ALN, MN and LC analyzed the data. NHS drafted the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
All procedures were performed in accordance with relevant regulation and approved by the Danish health research ethical committee (license no.: S-VF-19980072).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Putiri EL, Robertson KD. Epigenetic mechanisms and genome stability. Clin Epigenetics. 2011;2:299–314.View ArticlePubMedGoogle Scholar
- Meissner A. Epigenetic modifications in pluripotent and differentiated cells. Nat Biotechnol. 2010;28:1079–88.View ArticlePubMedGoogle Scholar
- Rakyan VK, Down TA, Balding DJ, Beck S. Epigenome-wide association studies for common human diseases. Nat Rev Genet. 2011;12:529–41.View ArticlePubMedPubMed CentralGoogle Scholar
- Laird PW. Principles and challenges of genomewide DNA methylation analysis. Nat Rev Genet. 2010;11:191–203.View ArticlePubMedGoogle Scholar
- Nishioka M, Bundo M, Kasai K, Iwamoto K. DNA methylation in schizophrenia: progress and challenges of epigenetic studies. Genome Med. 2012;4:96.View ArticlePubMedPubMed CentralGoogle Scholar
- Esteller M. Epigenetics in cancer. N Engl J Med. 2008;358:1148–59.View ArticlePubMedGoogle Scholar
- Guerrero-Preston R, Goldman LR, Brebi-Mieville P, Ili-Gangas C, Lebron C, Witter FR, Apelberg BJ, Hernandez-Roystacher M, Jaffe A, Halden RU, et al. Global DNA hypomethylation is associated with in utero exposure to cotinine and perfluorinated alkyl compounds. Epigenetics. 2010;5:539–46.View ArticlePubMedPubMed CentralGoogle Scholar
- Suter M, Ma J, Harris A, Patterson L, Brown KA, Shope C, Showalter L, Abramovici A, Aagaard-Tillery KM. Maternal tobacco use modestly alters correlated epigenome-wide placental DNA methylation and gene expression. Epigenetics. 2011;6:1284–94.View ArticlePubMedPubMed CentralGoogle Scholar
- Cao-Lei L, Massart R, Suderman MJ, Machnes Z, Elgbeili G, Laplante DP, Szyf M, King S. DNA methylation signatures triggered by prenatal maternal stress exposure to a natural disaster: Project Ice Storm. PLoS One. 2014;9:e107653.View ArticlePubMedPubMed CentralGoogle Scholar
- Melas PA, Rogdaki M, Osby U, Schalling M, Lavebratt C, Ekstrom TJ. Epigenetic aberrations in leukocytes of patients with schizophrenia: association of global DNA methylation with antipsychotic drug treatment and disease onset. FASEB J. 2012;26:2712–8.View ArticlePubMedGoogle Scholar
- Dong E, Nelson M, Grayson DR, Costa E, Guidotti A. Clozapine and sulpiride but not haloperidol or olanzapine activate brain DNA demethylation. Proc Natl Acad Sci U S A. 2008;105:13614–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Warner MJ, Ozanne SE. Mechanisms involved in the developmental programming of adulthood disease. Biochem J. 2010;427:333–47.View ArticlePubMedGoogle Scholar
- Tobi EW, Lumey LH, Talens RP, Kremer D, Putter H, Stein AD, Slagboom PE, Heijmans BT. DNA methylation differences after exposure to prenatal famine are common and timing- and sex-specific. Hum Mol Genet. 2009;18:4046–53.View ArticlePubMedPubMed CentralGoogle Scholar
- Meaney MJ, Ferguson-Smith AC. Epigenetic regulation of the neural transcriptome: the meaning of the marks. Nat Neurosci. 2010;13:1313–8.View ArticlePubMedGoogle Scholar
- Smith AK, Kilaru V, Kocak M, Almli LM, Mercer KB, Ressler KJ, Tylavsky FA, Conneely KN. Methylation quantitative trait loci (meQTLs) are consistently detected across ancestry, developmental stage, and tissue type. BMC Genomics. 2014;15:145.View ArticlePubMedPubMed CentralGoogle Scholar
- Aberg KA, Xie LY, McClay JL, Nerella S, Vunck S, Snider S, Beardsley PM, van den Oord EJ. Testing two models describing how methylome-wide studies in blood are informative for psychiatric conditions. Epigenomics. 2013;5:367–77.View ArticlePubMedPubMed CentralGoogle Scholar
- Bjornsson HT, Sigurdsson MI, Fallin MD, Irizarry RA, Aspelund T, Cui H, Yu W, Rongione MA, Ekstrom TJ, Harris TB, et al. Intra-individual change over time in DNA methylation with familial clustering. JAMA. 2008;299:2877–83.View ArticlePubMedPubMed CentralGoogle Scholar
- Hollegaard MV, Grauholm J, Nielsen R, Grove J, Mandrup S, Hougaard DM. Archived neonatal dried blood spot samples can be used for accurate whole genome and exome-targeted next-generation sequencing. Mol Genet Metab. 2013;110:65–72.View ArticlePubMedGoogle Scholar
- Bonnelykke K, Sleiman P, Nielsen K, Kreiner-Moller E, Mercader JM, Belgrave D, den Dekker HT, Husby A, Sevelsted A, Faura-Tellez G, et al. A genome-wide association study identifies CDHR3 as a susceptibility locus for early childhood asthma with severe exacerbations. Nat Genet. 2014;46:51–5.View ArticlePubMedGoogle Scholar
- Hollegaard MV, Grauholm J, Norgaard-Pedersen B, Hougaard DM. DNA methylome profiling using neonatal dried blood spot samples: a proof-of-principle study. Mol Genet Metab. 2013;108:225–31.View ArticlePubMedGoogle Scholar
- Beyan H, Down TA, Ramagopalan SV, Uvebrant K, Nilsson A, Holland ML, Gemma C, Giovannoni G, Boehm BO, Ebers GC, et al. Guthrie card methylomics identifies temporally stable epialleles that are present at birth in humans. Genome Res. 2012;22:2138–45.View ArticlePubMedPubMed CentralGoogle Scholar
- Aberg KA, Xie LY, Nerella S, Copeland WE, Costello EJ, van den Oord EJ. High quality methylome-wide investigations through next-generation sequencing of DNA from a single archived dry blood spot. Epigenetics. 2013;8:542–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Joo JE, Wong EM, Baglietto L, Jung CH, Tsimiklis H, Park DJ, Wong NC, English DR, Hopper JL, Severi G, et al. The use of DNA from archival dried blood spots with the Infinium HumanMethylation450 array. BMC Biotechnol. 2013;13:23.View ArticlePubMedPubMed CentralGoogle Scholar
- Cruickshank MN, Oshlack A, Theda C, Davis PG, Martino D, Sheehan P, Dai Y, Saffery R, Doyle LW, Craig JM. Analysis of epigenetic changes in survivors of preterm birth reveals the effect of gestational age and evidence for a long term legacy. Genome Med. 2013;5:96.View ArticlePubMedPubMed CentralGoogle Scholar
- Bundo M, Sunaga F, Ueda J, Kasai K, Kato T, Iwamoto K. A systematic evaluation of whole genome amplification of bisulfite-modified DNA. Clin Epigenetics. 2012;4:22.View ArticlePubMedPubMed CentralGoogle Scholar
- Ghantous A, Saffery R, Cros MP, Ponsonby AL, Hirschfeld S, Kasten C, Dwyer T, Herceg Z, Hernandez-Vargas H. Optimized DNA extraction from neonatal dried blood spots: application in methylome profiling. BMC Biotechnol. 2014;14:60.View ArticlePubMedPubMed CentralGoogle Scholar
- Clark C, Palta P, Joyce CJ, Scott C, Grundberg E, Deloukas P, Palotie A, Coffey AJ. A comparison of the whole genome approach of MeDIP-seq to the targeted approach of the Infinium HumanMethylation450 BeadChip((R)) for methylome profiling. PLoS One. 2012;7:e50233.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen G, Yu D, Chen J, Cao R, Yang J, Wang H, Ji X, Ning B, Shi T. Re-annotation of presumed noncoding disease/trait-associated genetic variants by integrative analyses. Sci Rep. 2015;5:9453.View ArticlePubMedPubMed CentralGoogle Scholar
- Feber A, Wilson GA, Zhang L, Presneau N, Idowu B, Down TA, Rakyan VK, Noon LA, Lloyd AC, Stupka E, et al. Comparative methylome analysis of benign and malignant peripheral nerve sheath tumors. Genome Res. 2011;21:515–24.View ArticlePubMedPubMed CentralGoogle Scholar
- Liu J, Morgan M, Hutchison K, Calhoun VD. A study of the influence of sex on genome wide methylation. PLoS One. 2010;5:e10028.View ArticlePubMedPubMed CentralGoogle Scholar
- Bock C, Tomazou EM, Brinkman AB, Muller F, Simmer F, Gu H, Jager N, Gnirke A, Stunnenberg HG, Meissner A. Quantitative comparison of genome-wide DNA methylation mapping technologies. Nat Biotechnol. 2010;28:1106–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Wu HC, Wang Q, Chung WK, Andrulis IL, Daly MB, John EM, Keegan TH, Knight J, Bradbury AR, Kappil MA, et al. Correlation of DNA methylation levels in blood and saliva DNA in young girls of the LEGACY Girls study. Epigenetics. 2014;9:929–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Shen H, Qiu C, Li J, Tian Q, Deng HW. Characterization of the DNA methylome and its interindividual variation in human peripheral blood monocytes. Epigenomics. 2013;5:255–69.View ArticlePubMedGoogle Scholar
- Lam LL, Emberly E, Fraser HB, Neumann SM, Chen E, Miller GE, Kobor MS. Factors underlying variable DNA methylation in a human community cohort. Proc Natl Acad Sci U S A. 2012;109 Suppl 2:17253–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Byun HM, Siegmund KD, Pan F, Weisenberger DJ, Kanel G, Laird PW, Yang AS. Epigenetic profiling of somatic tissues from human autopsy specimens identifies tissue- and individual-specific DNA methylation patterns. Hum Mol Genet. 2009;18:4808–17.View ArticlePubMedPubMed CentralGoogle Scholar
- Yang AS, Estecio MR, Doshi K, Kondo Y, Tajara EH, Issa JP. A simple method for estimating global DNA methylation using bisulfite PCR of repetitive DNA elements. Nucleic Acids Res. 2004;32:e38.View ArticlePubMedPubMed CentralGoogle Scholar
- Irahara N, Nosho K, Baba Y, Shima K, Lindeman NI, Hazra A, Schernhammer ES, Hunter DJ, Fuchs CS, Ogino S. Precision of pyrosequencing assay to measure LINE-1 methylation in colon cancer, normal colonic mucosa, and peripheral blood cells. J Mol Diagn. 2010;12:177–83.View ArticlePubMedPubMed CentralGoogle Scholar
- Taiwo O, Wilson GA, Morris T, Seisenberger S, Reik W, Pearce D, Beck S, Butcher LM. Methylome analysis using MeDIP-seq with low DNA concentrations. Nat Protoc. 2012;7:617–36.View ArticlePubMedGoogle Scholar
- Zhao MT, Whyte JJ, Hopkins GM, Kirk MD, Prather RS. Methylated DNA immunoprecipitation and high-throughput sequencing (MeDIP-seq) using low amounts of genomic DNA. Cell Reprogram. 2014;16:175–84.View ArticlePubMedGoogle Scholar
- Guo H, Zhu P, Guo F, Li X, Wu X, Fan X, Wen L, Tang F. Profiling DNA methylome landscapes of mammalian cells with single-cell reduced-representation bisulfite sequencing. Nat Protoc. 2015;10:645–59.View ArticlePubMedGoogle Scholar
- Guo H, Zhu P, Yan L, Li R, Hu B, Lian Y, Yan J, Ren X, Lin S, Li J, et al. The DNA methylation landscape of human early embryos. Nature. 2014;511:606–10.View ArticlePubMedGoogle Scholar
- McCarthy NS, Melton PE, Cadby G, Yazar S, Franchina M, Moses EK, Mackey DA, Hewitt AW. Meta-analysis of human methylation data for evidence of sex-specific autosomal patterns. BMC Genomics. 2014;15:981.View ArticlePubMedPubMed CentralGoogle Scholar
- Jones MJ, Goodman SJ, Kobor MS. DNA methylation and healthy human aging. Aging Cell. 2015;14:924–32.View ArticlePubMedPubMed CentralGoogle Scholar
- Flanagan JM, Popendikyte V, Pozdniakovaite N, Sobolev M, Assadzadeh A, Schumacher A, Zangeneh M, Lau L, Virtanen C, Wang SC, et al. Intra- and interindividual epigenetic variation in human germ cells. Am J Hum Genet. 2006;79:67–84.View ArticlePubMedPubMed CentralGoogle Scholar
- Sandovici I, Kassovska-Bratinova S, Loredo-Osti JC, Leppert M, Suarez A, Stewart R, Bautista FD, Schiraldi M, Sapienza C. Interindividual variability and parent of origin DNA methylation differences at specific human Alu elements. Hum Mol Genet. 2005;14:2135–43.View ArticlePubMedGoogle Scholar
- Salzberg SL, Phillippy AM, Zimin A, Puiu D, Magoc T, Koren S, Treangen TJ, Schatz MC, Delcher AL, Roberts M, et al. GAGE: a critical evaluation of genome assemblies and assembly algorithms. Genome Res. 2012;22:557–67.View ArticlePubMedPubMed CentralGoogle Scholar
- Jaffe AE, Irizarry RA. Accounting for cellular heterogeneity is critical in epigenome-wide association studies. Genome Biol. 2014;15:R31.View ArticlePubMedPubMed CentralGoogle Scholar
- Maksimovic J, Gagnon-Bartsch JA, Speed TP, Oshlack A. Removing unwanted variation in a differential methylation analysis of Illumina HumanMethylation450 array data. Nucleic Acids Res. 2015;43:e106.View ArticlePubMedPubMed CentralGoogle Scholar
- Hannon E, Lunnon K, Schalkwyk L, Mill J. Interindividual methylomic variation across blood, cortex, and cerebellum: implications for epigenetic studies of neurological and neuropsychiatric phenotypes. Epigenetics. 2015;10:1024–32.View ArticlePubMedPubMed CentralGoogle Scholar
- Klengel T, Mehta D, Anacker C, Rex-Haffner M, Pruessner JC, Pariante CM, Pace TW, Mercer KB, Mayberg HS, Bradley B, et al. Allele-specific FKBP5 DNA demethylation mediates gene-childhood trauma interactions. Nat Neurosci. 2013;16:33–41.View ArticlePubMedGoogle Scholar
- St Julien KR, Jelliffe-Pawlowski LL, Shaw GM, Stevenson DK, O’Brodovich HM, Krasnow MA, Stanford BPDSG. High quality genome-wide genotyping from archived dried blood spots without DNA amplification. PLoS One. 2013;8:e64710.View ArticlePubMedPubMed CentralGoogle Scholar
- Lienhard M, Grimm C, Morkel M, Herwig R, Chavez L. MEDIPS: genome-wide differential coverage analysis of sequencing data derived from DNA enrichment experiments. Bioinformatics. 2014;30:284–6.View ArticlePubMedGoogle Scholar