Inter-individual methylation variability in differentially methylated regions between maternal whole blood and first trimester CVS
© Ioannides et al.; licensee BioMed Central Ltd. 2014
Received: 22 September 2014
Accepted: 12 October 2014
Published: 1 November 2014
DNA methylation is the most studied form of epigenetic regulation, a process by which chromatin composition and transcription factor binding is altered to influence tissue specific gene expression and differentiation. Such tissue specific methylation patterns are investigated as biomarkers for cancer and cell-free fetal DNA using various methodologies.
We have utilized methylation DNA immunoprecipitation (MeDIP) and real-time quantitative PCR to investigate the inter-individual methylation variability of differentially methylated regions (DMRs) on chromosomes 18 and 21. We have characterized 15 newly selected and seven previously validated DMRs in 50, 1st trimester Chorionic villus samplings (CVS) and 50 female non-pregnant peripheral blood (WBF) samples. qPCR results from MeDIP and genomic DNA (Input) assays were used to calculate fold enrichment values for each DMR. For all regions tested, enrichment was higher in CVS than in WBF samples with mean enrichments ranging from 0.22 to 6.4 and 0.017 to 1 respectively. Despite inter-individual variability, mean enrichment values for CVS were significantly different than those for WBF in all DMRs tested (p < 0.01). This observation is reinforced by the absence of overlap in CVS and WBF enrichment value distributions for 15 of 22 DMRs.
Our work provides an expansion in the biomarker panel available for non-invasive prenatal diagnosis (NIPD) using the MeDIP-qPCR methology for Down syndrome and can eventually provide the starting point towards the development for assays towards the detection of Edwards syndrome. Furthermore, our data indicate that inter-experimental and inter-individual variation in methylation is apparent, yet the difference in methylation status across tissues is large enough to allow for robust tissue specific methylation identification.
In vertebrates DNA methylation is a conserved epigenetic modification by which DNA methyltransferases add a methyl group to carbon 5 of cytosine residues present in CpG dinucleotides. This modification is the most studied form of epigenetic regulation and has been strongly associated with chromosomal stability and imprinting control . Furthermore, this epigenetic process also regulates chromatin composition and transcription factor binding to directly influence transcriptional activity ,.
DNA methylation occurs primarily in CpG islands (GGIs) and shores both in coding and non-coding regions of the genome, with gene regulatory regions such as promoters and first exons being a frequent methylation target . Due to this integral relationship with gene expression regulation, DNA methylation patterns are very closely associated with developmental processes and differentiation. Consequently, DNA methylation directly modulates phenotype, and distinct methylation patterns have been associated with tissue specificity and a variety of disease states ranging from cancer to neurological disorders ,. These tissue specific differentially methylated regions (tDMRs) are currently under investigation for their utility as biomarkers for disease progression and prognosis, particularly in the field of cancer research, disease detection and response to treatment .
The discovery of cell free fetal DNA (cffDNA) in the maternal circulation has greatly facilitated the development of non-invasive prenatal diagnosis (NIPD) . The direct correlation between phenotype and DNA methylation patterns has allowed the use of DMRs as possible biomarkers in prenatal diagnosis. Several groups have utilized the methylation differences between placenta-derived cffDNA and maternal DNA in order to identify highly specific fetal DMR biomarkers for non-invasive prenatal diagnosis of aneuploidies. Previous studies employed a variety of methods including sodium bisulfite conversion and methylation sensitive restriction digestion, but yielded a relatively small number of fetal specific DMRs including the SERPINB5, RASSF1A and U-PDE9A genes -.
In 2009, Papageorgiou et al.  applied methylation DNA immunoprecipitation (MeDIP) coupled with high resolution tiling oligonucleotide array (Chip) analysis to identify DMRs between Chorionic villus sampling (CVS) and female peripheral blood DNA (WBF). They were able to identify thousands of DMRs on chromosomes 13, 18, 21, X, and Y including methylation sensitive restriction sites, CGIs and promoter regions. This MeDIP-Chip approach was the trigger for investigating the utility of MeDIP followed by real-time qPCR (MeDIP-qPCR) for the non-invasive prenatal diagnosis of trisomy 21, yielding 100% sensitivity and specificity . This novel NIPD method was validated by a second study of 175 cases again yielding high sensitivity and specificity .
The current study utilizes the MeDIP-qPCR methodology to expand our range of fetal specific DMR biomarkers by selecting and screening 15 additional DMRs on chromosomes 21 and 18. Special emphasis is given on investigating the methylation variability in different samples from these newly selected and previously reported DMRs - by screening them in a set of 50, 1st trimester CVS and 50 WBF. Overall, this work confirms the distinctively different methylation status of these regions in CVS and WBF.
Using the above criteria we identified a set of 40 candidate DMRs between CVS and WBF from the microarray data . This set was subsequently screened in a cohort of six CVS and six WBF to calculate the enrichment values for each DMR (Additional file 1). Based on this initial screening we were able to select the 15 regions with the highest CVS enrichment for further validation/characterization, using seven previously validated DMRs by Papageorgiou et al.  and Tsaliki et al. , as a comparison standard.
Ranking of DMRs tested according to the difference between mean enrichment values for each DMRs
Coefficient of variation WBF
Coefficient of variation CVS
Our study aimed to validate and characterize a set of differentially methylated regions between CVS and WBF, obtained from MeDIP-Chip data . The methylation characteristics of the 15 candidate DMRs, located on chromosomes 18 and 21, were ascertained in 50 CVS and 50 WBF using the MeDIP-qPCR methodology. To our knowledge this is the first MeDIP based biomarker screening study utilizing such a large sample set. None of the selected DMRs were located on CGIs, but within intergenic or intragenic regions. Such DMR distribution in non-coding intergenic and intragenic sequences is in agreement with data from a large scale investigation of tissue specific methylation profiles by Rakyan et al. , who reported differential methylation of intergenic and CpG poor promoter regions in addition to CGIs.
Based on our results, all 15 DMRs showed differential enrichment between the two tissues and 11 out of these 15 were strongly and consistently hypermethylated in CVS samples. The ability of these DMRs to distinguish between CVS and WBF was equivalent to that of the seven previously validated DMRs used as performance standards. In fact, the tissue discriminating performance of the DMRs tested here, shows close similarities with the previously validated DMRs as it is illustrated by our heat map distribution and the unsupervised clustering patterns obtained.
DMR enrichment values showed variability among the different samples. This is likely caused by combination of both inter-experimental technical variability and inter-individual methylation variability. The presence of variability in MeDIP based assays has previously been described by Butcher et al. . In addition, the issue of inter-individual DNA methylation variability has been the focus of several studies -. This high inter-individual variability has been attributed to a variety of factors including environmental conditions, diet, age and psychosocial factors -. Furthermore, it has been documented that regions with low CpG density, as the DMRs under investigation, show higher inter-individual variability as compared to regions with high CpG density, such as CGIs .
Others have also shown that methylation variability can coincide with tissue specific DMRs without obscuring the tissue discriminating properties of those DMRs . It is therefore of no surprise that despite the DNA methylation variability in our study, the newly validated set of 11 DMRs clearly distinguishes between CVS and WBF tissues.
Our work here substantially increases the number of confirmed chromosome 21 fetal specific DMRs, and therefore provides a significant expansion in the biomarker panel available for MeDIP-qPCR-based NIPD of Down syndrome. Such an expansion is predicted to further improve the robustness of the methodology and bolster its diagnostic classification power. It is also very important to note that our current study is the first to validate chromosome 18 fetal specific DMRs in a relatively large sample set. This small panel of chromosome 18 DMRs can potentially provide a very valuable testing platform on which future NIPD assays for Edwards syndrome will be developed.
NIPD has gained a lot of interest the last few years. Utilizing the methylation differences between fetal and maternal DNA, several groups have managed to identify biomarkers using different approaches. This study aimed to characterize and validate fetal specific methylated regions using the MeDIP-qPCR methodology. We were able to show that the selected regions had distinct methylation patterns between fetal and maternal tissue, despite inter-individual and inter-experimental variability. In addition, we have expanded the panel of the existing DMRs on chromosome 21 and have characterized a new set of markers on chromosome 18 which can provide the starting point towards the development for assays towards the detection of Edwards syndrome.
Human Samples and DNA preparation
WBF samples were obtained anonymously from 50 normal non-pregnant females 20-40 years of age. Fifty, 1st trimester CVS were obtained from the Department of Cytogenetic and Genomics at the Cyprus Institute of Neurology and Genetics (Nicosia, Cyprus). Protocols used for collecting samples for our study were approved by the appropriate Bioethics Committees, and informed consent was obtained from all participants. WBF and CVS samples were used to extract DNA using the QIAamp DNA blood midi kit or the QIAmp DNA mini kit according to the manufacturer’s instructions (QIAGEN, Hilden, Germany). All CVS underwent karyotyping and Quantitative-Fluorescent PCR (QF-PCR) analysis in order to confirm their normal status.
Ligation-mediated PCR (LM-PCR) and MeDIP assay
LM-PCR and MeDIP assays were conducted as described previously . Briefly, 2.5 μg of genomic DNA were sonicated using the Bioruptor Twin sonication system (UCD-400, Diagenode, Liege, Belgium) into fragments, 300-1000 bp in size. Fragment size was verified using agarose gel electrophoresis. The fragments were blunt-ended using HPLC water, 1X NEB buffer 2 (New England BioLabs, Ipswich, UK), 10X bovine serum albumin (New England BioLabs) 100 mmol/L dNTP mix (GE Healthcare, Little Chalfont, UK) and T4 DNA polymerase (3 U/μl; New England BioLabs). Fragments were purified using the QiAquick PCR purification kit (Qiagen) and linkers were then ligated onto the blunt ends by overnight incubation at 16°C with T4 DNA ligase (New England) and T4 DNA ligase buffer (New England). Overhangs were subsequently filled in by incubating at 72°C for 10 minutes with 100 mmol/L dNTP mix (GE Healthcare), 1X PCR gold buffer (Roche, Mannheim, Germany), 1.5 mmol/L MgCl2 (Roche) HPLC water and AmpliTaq DNA polymerase (Applied Biosystems, Branchburg, New Jersey, USA). 50 ng of ligated DNA was removed and kept as input DNA. The remaining ligated DNA (800-1200 ng) was subjected to MeDIP using 3 μg mouse anti-5’methylCytosine (a-5mC) antibody (Eurogentec Saraing, Belgium). Hypermethylated DNA bound to a-mC antibodies was magnetically captured using Dynabeads® M-280 Sheep Anti-Mouse IgG magnetic beads (Life technologies, Carlsbad, California, USA) and subsequently released using proteinase K (Roche). LM-PCR used 12 ng of each input and MeDIP DNA as described earlier .
Characteristics of the regions tested
Real-Time quantitative PCR (qPCR)
Primer sequences on DMRs tested
MeDIP enrichment values of the CVS and WBF samples were calculated for each region using the following equation:
Enrichment = e ∆Ct where e corresponds to the efficiency obtained in each real-time PCR reaction e = 10(-1/slope of STD curve)) and ∆Ct indicates the cycle difference between input DNA and MeDIP DNA [Ct(IN) – Ct(IP) ].
The mean enrichment values of each DMR were compared between WBF and CVS samples using the Mann-Whitney U tests  and the corresponding p-values were used to decide whether there was significant evidence to claim that the mean enrichments of the two groups were different.
Hierarchical clustering of the DMRs was conducted using an iterative algorithm that joins similar clusters based on the set of dissimilarities of the 100 individuals (calculating the Euclidean distanced between clusters) and re-computing their distances at each stage by the Lance-Williams dissimilarity update formula .
MI assisted in the design of the study, collected samples, carried out all the experiments and drafted the manuscript. EAP assisted in the design of the study. AK participated in the Methylation DNA Immunoprecipitation and qPCR assays. ET assisted in the biomarker discovery. CS carried out the statistical analysis. CS assisted in the design of the study. GK and MH assisted in the design of the study and helped to edit the manuscript. PCP conceived and coordinated the study. All authors read and approved the final manuscript.
The authors would like to thank Andrea Chrysostomou for her assistance in MeDIP and qPCR assays. This research was partially funded by the 7th Framework Programme as part of the ANGELAB project (http://cordis.europa.eu/projects/rcn/105552_en.html, project number 317635).
- Bird A: DNA methylation patterns and epigenetic memory. Genes Dev 2002, 16: 6–21. 10.1101/gad.947102View ArticlePubMedGoogle Scholar
- Gopalakrishnan S, Van Emburgh BO, Robertson KD: DNA methylation in development and human disease. Mutat Res 2008, 647: 30–38. 10.1016/j.mrfmmm.2008.08.006PubMed CentralView ArticlePubMedGoogle Scholar
- Baylin SB, Jones PA: A decade of exploring the cancer epigenome―biological and translational implications. Nat Rev Cancer 2011, 11: 726–734. 10.1038/nrc3130PubMed CentralView ArticlePubMedGoogle Scholar
- Laird PW: Principles and challenges of genomewide DNA methylation analysis. Nat Rev Cancer 2010, 11: 191–203. 10.1038/nrg2732View ArticleGoogle Scholar
- You JS, Jones PA: Cancer genetics and epigenetics: two sides of the same coin? Cancer Cell 2012, 22: 9–20. 10.1016/j.ccr.2012.06.008PubMed CentralView ArticlePubMedGoogle Scholar
- Jakovcevski M, Akbarian S: Epigenetic mechanisms in neurological disease. Nat Med 2012,18(8):1194–1204. 10.1038/nm.2828PubMed CentralView ArticlePubMedGoogle Scholar
- Rawson JB, Bapat B: Epigenetic biomarkers in colorectal cancer diagnostics. Expert Rev Mol Diagn 2012, 12: 499–509. 10.1586/erm.12.39View ArticlePubMedGoogle Scholar
- Lo YM, Corbetta N, Chamberlain PF, Rai V, Sargent IL, Redman CW, Wainscoat JS: Presence of fetal DNA in maternal plasma and serum. Lancet 1997, 350: 485–487. 10.1016/S0140-6736(97)02174-0View ArticlePubMedGoogle Scholar
- Della Ragione F, Mastrovito P, Campanile C, Conti A, Papageorgiou EA, Hulten MA, Patsalis PC, Carter NP, D’Esposito M: Differential DNA methylation as a tool for noninvasive prenatal diagnosis (NIPD) of X chromosome aneuploidies. J Mol Diagn 2010, 12: 797–807. 10.2353/jmoldx.2010.090199PubMed CentralView ArticlePubMedGoogle Scholar
- Old RW, Crea F, Puszyk W, Hulten MA: Candidate epigenetic biomarkers for non-invasive prenatal diagnosis of Down syndrome. Reprod Biomed Online 2007, 15: 227–235. 10.1016/S1472-6483(10)60713-4View ArticlePubMedGoogle Scholar
- Tong YK, Chiu RW, Akolekar R, Leung TY, Lau TK, Nicolaides KH, Lo YM: Epigenetic-genetic chromosome dosage approach for fetal trisomy 21 detection using an autosomal genetic reference marker. PLoS One 2010, 5: e15244. 10.1371/journal.pone.0015244PubMed CentralView ArticlePubMedGoogle Scholar
- Papageorgiou EA, Fiegler H, Rakyan V, Beck S, Hulten M, Lamnissou K, Carter NP, Patsalis PC: Sites of differential DNA methylation between placenta and peripheral blood: molecular markers for noninvasive prenatal diagnosis of aneuploidies. Am J Pathol 2009, 174: 1609–1618. 10.2353/ajpath.2009.081038PubMed CentralView ArticlePubMedGoogle Scholar
- Papageorgiou EA, Karagrigoriou A, Tsaliki E, Velissariou V, Carter NP, Patsalis PC: Fetal-specific DNA methylation ratio permits noninvasive prenatal diagnosis of trisomy 21. Nat Med 2011, 17: 510–513. 10.1038/nm.2312PubMed CentralView ArticlePubMedGoogle Scholar
- Tsaliki E, Papageorgiou EA, Spyrou C, Koumbaris G, Kypri E, Kyriakou S, Sotiriou C, Touvana E, Keravnou A, Karagrigoriou A, Lamnissou K, Velissariou V, Patsalis PC: MeDIP real-time qPCR of maternal peripheral blood reliably identifies trisomy 21. Prenat Diagn 2012, 32: 996–1001. 10.1002/pd.3947View ArticlePubMedGoogle Scholar
- Rakyan VK, Down TA, Thorne NP, Flicek P, Kulesha E, Graf S, Tomazou EM, Backdahl L, Johnson N, Herberth M, Howe KL, Jackson DK, Miretti MM, Fiegler H, Marioni JC, Birney E, Hubbard TJ, Carter NP, Tavare S, Beck S: An integrated resource for genome-wide identification and analysis of human tissue-specific differentially methylated regions (tDMRs). Genome Res 2008,18(9):1518–1529. 10.1101/gr.077479.108PubMed CentralView ArticlePubMedGoogle Scholar
- Butcher LM, Beck S: AutoMeDIP-seq: a high-throughput, whole genome, DNA methylation assay. Methods 2010, 52: 223–231. 10.1016/j.ymeth.2010.04.003PubMed CentralView ArticlePubMedGoogle Scholar
- Feinberg AP, Irizarry RA: Evolution in health and medicine Sackler colloquium: Stochastic epigenetic variation as a driving force of development, evolutionary adaptation, and disease. Proc Natl Acad Sci U S A 2010,1(Suppl 107):1757–1764. 10.1073/pnas.0906183107View ArticleGoogle Scholar
- El-Maarri O, Walier M, Behne F, van Uum J, Singer H, Diaz-Lacava A, Nusgen N, Niemann B, Watzka M, Reinsberg J, van der Ven H, Wienker T, Stoffel-Wagner B, Schwaab R, Oldenburg J: Methylation at global LINE-1 repeats in human blood are affected by gender but not by age or natural hormone cycles. PLoS One 2011, 6: e16252. 10.1371/journal.pone.0016252PubMed CentralView ArticlePubMedGoogle Scholar
- Choi SH, Worswick S, Byun HM, Shear T, Soussa JC, Wolff EM, Douer D, Garcia-Manero G, Liang G, Yang AS: Changes in DNA methylation of tandem DNA repeats are different from interspersed repeats in cancer. Int J Cancer 2009, 125: 723–729. 10.1002/ijc.24384PubMed CentralView ArticlePubMedGoogle Scholar
- Fryer AA, Emes RD, Ismail KM, Haworth KE, Mein C, Carroll WD, Farrell WE: Quantitative, high-resolution epigenetic profiling of CpG loci identifies associations with cord blood plasma homocysteine and birth weight in humans. Epigenetics 2011, 6: 86–94. 10.4161/epi.6.1.13392PubMed CentralView ArticlePubMedGoogle Scholar
- Wong CC, Caspi A, Williams B, Craig IW, Houts R, Ambler A, Moffitt TE, Mill J: A longitudinal study of epigenetic variation in twins. Epigenetics 2010, 5: 516–526. 10.4161/epi.5.6.12226PubMed CentralView 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,2(Suppl 109):17253–17260. 10.1073/pnas.1121249109View ArticleGoogle Scholar
- Schneider E, Pliushch G, El Hajj N, Galetzka D, Puhl A, Schorsch M, Frauenknecht K, Riepert T, Tresch A, Muller AM, Coerdt W, Zechner U, Haaf T: Spatial, temporal and interindividual epigenetic variation of functionally important DNA methylation patterns. Nucleic Acids Res 2010, 38: 3880–3890. 10.1093/nar/gkq126PubMed CentralView ArticlePubMedGoogle Scholar
- Bock C, Walter J, Paulsen M, Lengauer T: Inter-individual variation of DNA methylation and its implications for large-scale epigenome mapping. Nucleic Acids Res 2008, 36: e55. 10.1093/nar/gkn122PubMed CentralView ArticlePubMedGoogle Scholar
- MacDonald JR, Ziman R, Yuen RK, Feuk L, Scherer SW: The Database of Genomic Variants: a curated collection of structural variation in the human genome. Nucleic Acids Re 2014, 42: D986-D992. 10.1093/nar/gkt958View ArticleGoogle Scholar
- Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol 2000, 132: 365–386.PubMedGoogle Scholar
- Mann HB, Whitney DR: On a test of whether one of two random variables is stochastically larger than the other. The annals of mathematical statistics 1947,18(1):50–60. 10.1214/aoms/1177730491View ArticleGoogle Scholar
- Lance GN, Williams WT: A general theory of classificatory sorting strategies II. Clustering systems. The computer journal 1967, 10: 271–277. 10.1093/comjnl/10.3.271View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.