Inter-individual methylation variability in differentially methylated regions between maternal whole blood and first trimester CVS
Molecular Cytogenetics volume 7, Article number: 73 (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.
This DMR validation study was conducted on a set of 50 CVS and 50 WBF samples using the MeDIP–qPCR methodology (Table 1), the efficiency of which was monitored using one hypermethylated (HYPER) and one hypomethylated (HYPO) control regions. The HYPER is a region that showed hypermethylation for both CVS and WBF, while the HYPO is a region that showed hypomethylation for the two tissues . Enrichment values for HYPO were low in WBF and CVS samples while the HYPER control region showed enrichment for CVS and WBF with mean enrichment values of 3.12 and 3.22 respectively, indicating that the MeDIP procedure was highly specific for the methylated regions. Moreover, the previously validated DMRs performed as previously described , exhibiting distinctively different enrichment between CVS and WBF.
All tested DMRs showed a significant enrichment (p < 0.01) in CVS compared to those of WBF (Table 1). We compared the performance of the 15 newly selected DMRs with the previously validated set and we were able to determine that 11 of 15 DMRs showed enrichment values higher than the lowest of the previously validated DMRs, ranging from 1.9 to 6.4. Additional comparison of the two DMR sets also illustrated that for 11 of these 15 regions the difference of means (mean enrichment CVS – mean enrichment WBF) was again higher than the respective values of the validated DMRs (ranging from 1.6 to 6.4) (Table 1). Our analysis also shows that the enrichment distributions for CVS and WBF have no overlap for these 11 DMRs (Figure 1).
To better investigate tissue specificity (CVS-WBF) in the 15 newly selected DMRs in relation to the previously validated DMRs, we also constructed a heat map and hierarchical clustering of the 50 CVS and 50 WBF samples based on the obtained enrichment values (Figure 2). This analysis shows a clear differentiation between the two tissue types based on the obtained enrichment values. Furthermore, DMR clustering analysis showed that there was no distinct clustering separation between the newly selected and the previously validated DMRs.
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 .
Candidate DMRs on chromosomes 18 and 21 were selected from a set of potential differentially methylated regions previously described  according to the following three criteria: a) the region included at least three consecutive microarray probes, b) array results showed consistent DNA hypermethylation in first and third trimester placentas and hypomethylation in WBF samples, c) the region did not include segmental duplications and copy number variable regions based on the Database for Genomic Variants (DGV) . The regions considered for this paper are shown in Table 2.
Real-Time quantitative PCR (qPCR)
Primer design, optimal primer concentration experiments, and efficiency (e) of each qPCR reaction were performed as previously described  with the following modifications. Each qPCR reaction was performed on 8 ng of template DNA using SYBR Green PCR mastermix (Eurogentec) in a final reaction volume of 10 μl, using a BIORAD CFX 384 Real time system (BIORAD, Hercules, California). Each MeDIP, or Input template DNA was used to prepare three replicate qPCR reactions that were used for calculating the average Ct value for each template. Primer3 software  was used to design the qPCR primer sets that were synthesized by Sigma-Aldrich (Munich, Germany). Primer sets utilized for this body of work are listed in Table 3.
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.
Bird A: DNA methylation patterns and epigenetic memory. Genes Dev 2002, 16: 6–21. 10.1101/gad.947102
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.006
Baylin SB, Jones PA: A decade of exploring the cancer epigenome―biological and translational implications. Nat Rev Cancer 2011, 11: 726–734. 10.1038/nrc3130
Laird PW: Principles and challenges of genomewide DNA methylation analysis. Nat Rev Cancer 2010, 11: 191–203. 10.1038/nrg2732
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.008
Jakovcevski M, Akbarian S: Epigenetic mechanisms in neurological disease. Nat Med 2012,18(8):1194–1204. 10.1038/nm.2828
Rawson JB, Bapat B: Epigenetic biomarkers in colorectal cancer diagnostics. Expert Rev Mol Diagn 2012, 12: 499–509. 10.1586/erm.12.39
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-0
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.090199
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-4
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.0015244
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.081038
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.2312
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.3947
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.108
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.003
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.0906183107
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.0016252
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.24384
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.13392
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.12226
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.1121249109
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/gkq126
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/gkn122
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/gkt958
Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol 2000, 132: 365–386.
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/1177730491
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.271
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).
GK and EAP are currently employed by, and own shares of NIPD Genetics LTD which is developing a new non-invasive prenatal diagnosis for trisomy 21 based on the MeDIP-qPCR methodology. MI, AK, MH and CS declare that they have no competing interests.
At the time that the study was carried out ET and CS were employed by NIPD Genetics LTD. Currently, they have no affiliations with NIPD Genetics LTD and they declare no competing interests.
EAP and PCP have filed a PCT patent application for the MeDIP real time qPCR based NIPD approach (PCT Patent Application No. PCT/1B2011/000217).
Electronic supplementary material
About this article
Cite this article
Ioannides, M., Papageorgiou, E.A., Keravnou, A. et al. Inter-individual methylation variability in differentially methylated regions between maternal whole blood and first trimester CVS. Mol Cytogenet 7, 73 (2014). https://doi.org/10.1186/s13039-014-0073-8
- Non-invasive prenatal diagnosis
- Inter-individual variability
- Differentially methylated regions