- Technical advance
- Open Access
An improved sequencing-based strategy to estimate locus-specific DNA methylation
BMC Cancer volume 15, Article number: 639 (2015)
DNA methylation is an important epigenetic mechanism of transcriptional control that plays an essential role in several cellular functions. Aberrant DNA methylation in cancer has been frequently associated with downregulation of microRNAs and protein coding genes, such as miR-200c/miR-141 cluster and E-cadherin. Current strategies to assess DNA methylation, including bisulfite treatment-based assays, tend to be time-consuming and may be quite expensive when a precise appraisal is required. The Sanger-sequencing of the amplified bisulfite-treated DNA (BSP) might represent a practical option to measure DNA methylation at single CpG resolution. However, this strategy often produces noisy data, which affects accurate quantification. Here we propose an improved, reliable and cost-effective BSP-based protocol that allows proper DNA methylation assessment.
Our strategy, named normalized-BSP (NBSP), takes advantage of tailed C-balanced primers and a normalization procedure based on C/T ratio to overcome BSP-associated noise problems and nucleotide signal unbalance. NBSP was applied to estimate miR-200c/miR-141 locus methylation in serial dilution experiments and was compared to conventional methods. Besides, it was applied in the analysis of FFPE breast cancer samples and further validated in the context of the E-cadherin promoter.
NBSP strategy outperformed conventional BSP in the estimate of the fraction of methylated cytosine in serial dilution experiments, providing data in agreement with the widely used but cumbersome cloning-based protocol. This held true for both miR-200c/miR-141 locus and E-cadherin promoter analyses. Moreover, the miR-200c/miR-141 locus methylation reflected the decrease in miRNA expression both in breast cancer cell lines and in the FFPE samples.
NBSP is a rapid and economical method to estimate the extent of methylation at each CpG of a given locus. Notably, NBSP works efficiently on FFPE samples, thus disclosing the perspective of its application also in the diagnostic setting.
DNA methylation, one of the best-characterized epigenetic modifications, consists in the addition of a methyl group to cytosines included in CpG dinucleotides. The methylation of CpG islands (CGI), which are common in promoter regions, correlates with gene transcriptional repression [1, 2]. Aberrant DNA methylation is typically observed in tumors where it occurs at both protein coding gene and microRNA (miRNA) loci [3–5].
Several technologies have been developed to profile the methylation at CGI. These include comprehensive but expensive next-generation sequencing-based approaches (i.e.: WGBS [6, 7], RRBS , MethylCap-seq  or MBD-seq  as well as array- and PCR-based methods, more affordable and still used [11, 12]. Most techniques rely on the bisulfite conversion of unmethylated cytosine to uracil, and thus to thymine after PCR, leaving unaltered the methylated cytosine .
Rapid and simple methods to detect the ratio between C and T include the Sanger sequencing of PCR products of bisulfite-treated DNA (BSP). However, this approach fails to provide a quantitative measure of methylation because of high background noise and overscaled cytosine signals due to the DNA sequencing software that artificially adjusts signal strengths of underrepresented bases . On the other hand, the cloning and subsequent Sanger sequencing of the PCR clones (cloning-based sequencing method) , although more accurate, is time-consuming and expensive, as it needs the sequencing of a significant number of clones for statistical accuracy .
Here we report an enhanced Sanger sequencing-based protocol for quantifying CGI promoter methylation based on the use of 5’-end tailed PCR primers that allow for the improvement of both signal-to-noise and C/T ratio. The method was successfully applied to assess the methylation status of both a miRNA locus (miR-200c/miR-141) and the promoter of E-cadherin and was also suitable for the analysis of FFPE tumor samples.
MiR-200 is a tumor suppressor miRNA family that includes five members clustered and expressed as two separate polycistronic pri-miRNAs: the miR-200a/miR-200b/miR-429 cluster, mapping at 1p36; and the miR-200c/miR-141 cluster, at 12p13 [17–19]. Promoter hypermethylation has been reported to play a crucial role in the downregulation of miR-200 [20–22] that has been associated with malignancy, increased chemo- and radio-resistance, invasiveness and transition of carcinomas from epithelial towards a mesenchymal phenotype (EMT) [23–29]. A hallmark of EMT is the downregulation of the cell-cell adhesion protein E-cadherin (E-cad) , whose low expression, as a result of promoter hypermethylation, has been described in diverse carcinoma subtypes and is associated with poor prognosis [31–33].
MDA-MB-231, MDA-MB-157 and MCF7 were obtained from ATCC (American Type Culture Collection) and cultured as previously described .
Patients and samples
Formalin-fixed paraffin-embedded (FFPE) specimens from 3 breast cancers were collected at the CRO Aviano National Cancer Institute Biobank under patients’ informed consent. The use of tumor samples for this study was approved by the Institutional Review Board. Two 20 μm-slides with tumor cellularity greater than 70 %, as evaluated by a breast cancer pathologist (TP), were used per each case. Total RNA and DNA were isolated using the Recover All Total Nucleic Acid Isolation Kit (Life Technologies) according to the manufacturer’s instructions.
RNA extraction and qRT-PCR
Total RNA from cell lines was isolated using TRIzol (Life Technologies). MiRNA was reverse-transcribed and qRT-PCR performed using the TaqMan MicroRNA Assay kits specific for miR-200c and RNU48 (Life Technologies) and TaqMan Universal Master Mix (Life Technologies) according to the manufacturer's guidelines. miRNA levels were normalized with RNU48 and relative levels were calculated using the ΔΔCt method. Three independent experiments were performed in triplicate.
DNA extraction and bisulfite conversion
Genomic DNA was extracted from cell lines using the EZ1 DNA Tissue Kit (Qiagen). Bisulfite conversion of DNA (500 ng - 1 μg), obtained from cell lines and tissues, was carried out with the EpiTect Bisulfite kit (Qiagen), according to the manufacturer’s instructions.
Bisulfite PCR amplification
The region of the miR-200c/miR-141 locus, spanning from position −353 to −108 relative to the pre-miRNA-200c first nucleotide (chromosome12:7,072,510:7,072,755; Fig. 1a) and the promoter region from −115 to +54 nucleotide relative to the transcriptional start site of E-cad (CDH1 gene; chromosome16:68,771,079: 68,771,249; Fig. 5a) were amplified with primers specifically designed by MethPrimer (Additional file 1) .
5’-end tailed primers were obtained by adding at the 5’-end of the 200c-BSP-F and 200c-BSP-R a tail derived from the M13 (Tail1) or from the Decipher Project barcode library (Tail2-6; http://www.decipherproject.net). Tail1, Tail3 and Tail5 were added to the forward oligo and Tail2, Tail4 and Tail6 to the reverse oligo (Additional file 1). Tails 2–6 were randomly chosen among barcodes devoid of C or G at the 5’-end and in which each base is roughly equally represented (22-28 %). The tails, by contributing with C and T (G and A in the reverse primer) allow for compensation in the elaboration process. Primers with Tail1, Tail3 and Tail5 were used in combination with primers with Tail2, Tail4 and Tail6, respectively. All the three couples of primers well amplified miR-200c/miR-141 locus (Additional file 2). We selected Tail1- and Tail2-primers for this work. Tail1 and Tail2 were also added to E-cad-BSP forward and reverse oligo, respectively.
0.7-1 μl of bisulfite-treated DNA were amplified by using GoTaq® Polymerase (Promega) if not otherwise specified. The PCR amplification was performed in 20 μl reaction volume containing GoTaq® Green Master Mix 1X, 250 nM forward and reverse primers and with the following protocol: 95 °C for 4 min, 40 X [95 °C for 45 s, 60 °C (E-cad) or 62 °C (miR-200c) for 1 min and 30 s, 72 °C for 1 min 30 s], 72 °C for 4 min. Phusion U Hot Start DNA polymerase (Thermoscientific) was tested for the amplification of miR-200c/miR-141 locus of genomic DNA standards. The PCR amplification of 20 μl mixture containing 0.7-1 μl of bisulfite-treated DNA, 0.4 U Phusion U Hot Start DNA Polymerase, 400 nM forward and reverse primers, Phusion HF Buffer 1X and 200 μM dNTPs was performed with the following protocol: 98 °C for 1 min, 37 X [98 °C for 10 s, 64 °C for 15 s, 72 °C for 30 s] 72 °C for 5 min. 10 μl of PCR products were size-checked on a 2 % agarose gel and 5 μl were purified with 2 μl of ExoSap-IT (Affymetrix).
PCR cloning and assessment of methylation
Bisulfite-treated DNA was amplified by PCR with untailed primers (Additional file 1) and 1 μl of the PCR was directly cloned into the pCR2.1 vector using TA Cloning Kit (Life Technologies), according to the manufacturer’s protocol. Plasmids DNA from at least 20 colonies were isolated using the QiaPrep Spin Plasmid Miniprep kit (Qiagen) and sequenced. The methylation level for each CpG was deducted by dividing the number of C at each CpG site for the total number of clones sequenced.
Generation of DNA standards
We generated plasmid and genomic DNA standards to mimic different methylation levels of miR-200c/miR-141 locus. To obtain the plasmid DNA standards, miR-200c/miR-141 locus was amplified from bisulfite-converted DNA of MCF7 and MDA-MB-157 (unmethylated and 97 % methylated, respectively, as determined by the cloning method) and cloned into the pCR2.1 vector (TA Cloning Kit, Life Technologies) according to the manufacturer’s protocol. Two of these clones derived from completely methylated and unmethylated (for all CpG sites) template, respectively, were mixed to mimic different DNA methylation percentages: 0, 12.5, 25, 55, 75, 87.5 and 100 %. The C/T ratio, calculated as described below, was confirmed by plasmid direct sequencing (Additional file 3).
Moreover, a set of the genomic DNA standards was generated by mixing the bisulfite-treated DNAs of the aforementioned cell lines in order to obtain the following methylation levels: 0, 12.1, 24.2, 48.4, 72.6 and 97 %.
Sequencing reactions (10 μl) were performed using 1 μl of ExoSap-IT-purified PCR amplicons or 500 ng of plasmids, 2 μl of BigDye Terminator v.3.1 kit (Life Technologies), 300 nM sequencing primer, corresponding to 200c-BSP-F, Tail1 or Tail2 (Additional file 1), and the following protocol: 95 °C for 5 min, 25 X [95 °C 30 s, 50 °C for 30 s and 60 °C for 1 min and 30 s]. The sequencing reactions were then purified using the BigDye XTerminator Purification kit and ran on an ABI prism 3130 Genetic Analyzer (Applied Biosystems). SeqScape® Software v2.5 with the KB™ basecaller software or Chromas Lite Version 2.1.1 were used for data analysis.
Quantification of methylation by BSP
DNA standards and bisulfite-treated DNA were amplified by PCR with tailed primers (Additional file 1) and sequenced as described. The percentage of methylation at each CpG site was calculated as 100 ∗ C/(C + T), i.e. 100 times the ratio between the peak height of C on the sequencing chromatograms and the sum of peak height for C and T .
Quantification of methylation by NBSP
DNA standards and bisulfite-treated DNA were amplified by PCR with 5’-end tailed primers (Additional file 1) and sequenced as above. To adjust the overscaled C signals in the sequencing chromatograms we introduced a normalization factor (NF), based on the ratio of the signals for the C and T encoded by the tails of primers. Specifically, NF was calculated as the ratio between the mean of the peak height of the C and T read in sense direction on the sequence of Tail2 (corresponding to G and A in Tail2 reverse primer sequence; Additional file 1, Figs. 1c and 5c).
Then, the peak height of each C (Ci) included in the target sequence was corrected for this NF as follow: Cnorm = Ci/NF. Finally, the normalized C signals were used to determine the methylation percentage as described above, i.e. 100 ∗ Cnorm/(Cnorm + T).
The concordance between observed and expected values was analyzed by using the approach recommended by Bland and Altman [37, 38]. For all Bland–Altman plots, the mean percentage difference between the observed and expected results (mean bias) with associated 97.5 % confidence intervals and limits of agreement (±1.96 SD) were calculated (GraphPad Prism software).
For the analyses of miR-200c/miR-141 promoter methylation we focused on the region referred to as relevant for transcription (−353 to −108, relative to the pre-miRNA-200c first nucleotide) and that comprised 14 CpG sites (Fig. 1a) [17, 18].
We first performed the analysis of this region in a set of plasmid DNA standards obtained by mixing defined amount of clones corresponding to methylated and unmethylated DNA (see Methods). The direct sequencing of the PCR products of these standards displayed overscaled C signals and a high background noise that prevented the actual estimate of miR-200c/miR-141 promoter methylation (Fig. 1b and Additional file 4A).
In order to improve the quality of the sequencing traces, we amplified the aforementioned standards with 5’-end tailed primers (Fig. 1c and Additional files 1 and 4B) characterized by at least 4 C in the tails. The sequencing of these PCR products (BSP) provided chromatograms without any or only minimal background (Fig. 2a and b). Still, the measure of methylation extent was unsatisfactory. In fact, especially in the presence of low-intermediate levels of methylation, the C signals (i.e. non-converted, methylated cytosines) were overscaled, which resulted in an overestimate of DNA methylation. In fact, the mean bias (i.e. average percent difference between the observed and expected methylation levels) was 7.93 (limits of agreement from −13.66 to 29.52; Fig. 2c left panel). It is worth to note that the clone from unmethylated DNA displayed a G > A transition (at position 7,072,604 in the miR-200c/miR-141 locus). The ratio between G and A of each standard reflected the expected methylation levels suggesting the goodness of the standards (data not shown).
To overcome the C overestimation, we introduced a normalization strategy (referred in text as Normalized BSP, NBSP) that took into account the elaboration of overall nucleotide signals by the DNA sequencing software. Based on the assumption that, for any given sequence and in the absence of altering factors, the relation between mean of the peak height of two nucleotides, namely C and T, should be relatively constant, we calculated the ratio between C and T within the tail of the primers and used this ratio to normalize the overscaled C signals of the sequence (see Methods). The introduction of this normalization step significantly improved the estimate of methylation rate reducing the mean bias to −1.02 (Fig. 2c right panel; limits of agreement from −2.71 to 0.66).
Next we validated our strategy on the genomic DNA standards. Uracil present in the bisulfite-converted DNA may impair the DNA polymerase activity of Taq polymerase. Thus we compared the results obtained with Taq polymerase and with an uracil tolerant enzyme (Phusion U Hot start DNA polymerase). The two DNA polymerases showed similar results (Additional file 4C-F) and, importantly, NBSP displayed an improvement in the assessment of the methylation rate of genomic DNA standards compared to BSP in both analyses (Fig. 2d, Additional file 4G-H).
To further validate our signal normalization approach, we compared the performance of BSP and NBSP to the cloning-based sequencing method. According to the standard BSP procedure, the percentages of methylation at each CpG of the miR-200c/miR-141 locus in the MDA-MB-231 breast cancer cell line ranged from ~ 30 to 80 %. These were globally greater than those gauged by the cloning-based method, particularly for low and intermediate CpG methylation (Fig. 3a-c). NBSP outperformed the BSP, providing estimate close to those of the cloning procedure for the majority of CpG sites. Forward and reverse tailed primers provided similar results, both in terms of percentages of methylation and extent of the normalization factors (Additional file 5).
The partial methylation of the miR-200c/miR-141 locus in MDA-MB-231 corresponded to a limited expression of miR-200c compared to the unmethylated MCF7 and the fully methylated MDA-MB-157 (Additional file 6).
A similar inverse association between miR-200c expression and locus methylation was observed also when NBSP was applied to FFPE breast tumor samples, particularly for the CpG from −223 to −135 (Fig. 4a-d). A normal breast tissue sample showed only one methylated CpG and, as expected, expressed high levels of miR-200c.
Finally, we investigated the methylation pattern of E-cad, a typical gene silenced by DNA-hypermethylation. Our study focused on a well-defined CGI spanning between −115 and +54 nucleotides from transcription start site of the E-cad promoter (Fig. 5a). Again, NBSP outperformed BSP in the measure of E-cad promoter methylation in MDA-MB-231 and provided data similar to those obtained with standard cloning-based method (Fig. 5b-d).
Epigenetic inactivation of tumor suppressor genes is a frequent event that drives tumorigenic initiation and progression [39–41]. The increasing interest in the evaluation of miR-200c/miR-141 locus methylation as a measure of cancer progression [42, 43], prompted us to set up a reliable, fast and affordable method for the assessment of DNA methylation.
The NBSP method here proposed relies on the use of 5’-end tailed primers that reintroduce ‘true’ C, improve the quality of sequencing traces and allow C/T signal normalization. We implemented the normalization procedure because of the overscaled C signals engendered by the sequencing software which, during raw data elaboration, tends to artificially enhance the signal of underrepresented C resulted from the bisulfite treatment. Overestimation of C may also be caused by preferential amplification of methylated alleles, though it occurs more rarely than the PCR bias favoring the unmethylated ones [44, 45]. Furthermore, it has been reported that tailed primers could introduce amplification bias depending on the template . However, we can exclude these biases since the amplification of plasmid DNA standards harboring a G > A variant produced the expected G/A ratio. Nevertheless, we cannot rule out that the chosen tails, which work well with the two genes we analyzed, unevenly perform with other genes.
A number of studies have proposed alternative solutions for analyzing the Sanger-sequencing data, but their algorithms are often overwhelming [14, 44]. Our approach can be easily used and, importantly, it yields an estimate of methylation at each CpG site in agreement with data obtained with the conventional but cumbersome cloning-based method. Moreover, locus methylation as assessed by NBSP well reflected the miRNA expression in FFPE breast cancer samples. Importantly, NBSP allowed an accurate detection of methylation rate close to 10 %, a level below which methylation has negligible effects on miR-200c/miR-141 expression . Finally, NBSP can be applied to other genes, such here shown for E-cadherin.
We have presented here a reliable and cost-effective method to detect the methylation level of several CpGs. Our approach well performed in the analysis of the miR-200c/miR-141 locus and of the E-cad promoter, genes downregulated by methylation in a number of carcinoma. Besides, NBSP also works with FFPE tissues and thus may provide a viable and affordable tool to detect DNA methylation both for research and for diagnostic purposes.
Direct sequencing of PCR from bisulfite-treated DNA
Epithelial to mesenchymal transition
Formalin-Fixed Paraffin-Embedded (tissue)
Methylated DNA Binding Domain sequencing
Methylated-DNA capture sequencing
Reduced representation bisulfite sequencing
Whole-genome bisulfite sequencing
Baylin SB, Jones PA. A decade of exploring the cancer epigenome — biological and translational implications. Nat Rev Cancer. 2011;11:726–34.
You JS, Jones PA. Cancer genetics and epigenetics: two sides of the same coin? Cancer Cell. 2012;22:9–20.
Li J, Jin H, Wang X. Epigenetic Biomarkers: Potential Applications in Gastrointestinal Cancers. ISRN Gastroenterol. 2014;2014:464015 .
Suzuki H, Maruyama R, Yamamoto E, Kai M. Epigenetic alteration and microRNA dysregulation in cancer. Front Genet. 2013;4:258.
Klajic J, Fleischer T, Dejeux E, Edvardsen H, Warnberg F, Bukholm I, et al. Quantitative DNA methylation analyses reveal stage dependent DNA methylation and association to clinico-pathological factors in breast tumors. BMC Cancer. 2013;13:456.
Cokus SJ, Feng S, Zhang X, Chen Z, Merriman B, Haudenschild CD, et al. Shotgun bisulphite sequencing of the Arabidopsis genome reveals DNA methylation patterning. Nature. 2008;452:215–9.
Lister R, O’Malley RC, Tonti-Filippini J, Gregory BD, Berry CC, Millar AH, et al. Highly integrated single-base resolution maps of the epigenome in arabidopsis. Cell. 2008;133:523–36.
Meissner A, Gnirke A, Bell GW, Ramsahoye B, Lander ES, Jaenisch R. Reduced representation bisulfite sequencing for comparative high-resolution DNA methylation analysis. Nucleic Acids Res. 2005;33:5868–77.
Brinkman AB, Simmer F, Ma K, Kaan A, Zhu J, Stunnenberg HG. Whole-genome DNA methylation profiling using MethylCap-seq. Methods. 2010;52:232–6. DNA Methylation Analysis.
Serre D, Lee BH, Ting AH. MBD-isolated Genome Sequencing provides a high-throughput and comprehensive survey of DNA methylation in the human genome. Nucleic Acids Res. 2010;38:391–9.
Mensaert K, Denil S, Trooskens G, Van Criekinge W, Thas O, De Meyer T. Next-generation technologies and data analytical approaches for epigenomics. Environ Mol Mutagen. 2014;55:155–70.
Kim Y-J, Park S-W, Kim T-H, Park J-S, Cheong HS, Shin HD, et al. Genome-wide methylation profiling of the bronchial mucosa of asthmatics: relationship to atopy. BMC Med Genet. 2013;14:39.
Frommer M, McDonald LE, Millar DS, Collis CM, Watt F, Grigg GW, et al. A genomic sequencing protocol that yields a positive display of 5-methylcytosine residues in individual DNA strands. Proc Natl Acad Sci U S A. 1992;89:1827–31.
Lewin J, Schmitt AO, Adorján P, Hildmann T, Piepenbrock C. Quantitative DNA methylation analysis based on four-dye trace data from direct sequencing of PCR amplificates. Bioinforma Oxf Engl. 2004;20:3005–12.
Clark SJ, Statham A, Stirzaker C, Molloy PL, Frommer M. DNA methylation: Bisulphite modification and analysis. Nat Protoc. 2006;1:2353–64.
Myöhänen S, Wahlfors J, Jänne J. Automated fluorescent genomic sequencing as applied to the methylation analysis of the human ornithine decarboxylase gene. DNA Seq J DNA Seq Mapp. 1994;5:1–8.
Bracken CP, Gregory PA, Kolesnikoff N, Bert AG, Wang J, Shannon MF, et al. A double-negative feedback loop between ZEB1-SIP1 and the microRNA-200 family regulates epithelial-mesenchymal transition. Cancer Res. 2008;68:7846–54.
Burk U, Schubert J, Wellner U, Schmalhofer O, Vincan E, Spaderna S, et al. A reciprocal repression between ZEB1 and members of the miR-200 family promotes EMT and invasion in cancer cells. EMBO Rep. 2008;9:582–9.
Saini HK, Griffiths-Jones S, Enright AJ. Genomic analysis of human microRNA transcripts. Proc Natl Acad Sci. 2007;104:17719–24.
Castilla MÁ, Díaz-Martín J, Sarrió D, Romero-Pérez L, López-García MÁ, Vieites B, et al. MicroRNA-200 family modulation in distinct breast cancer phenotypes. PloS One. 2012;7, e47709.
Ceppi P, Mudduluru G, Kumarswamy R, Rapa I, Scagliotti GV, Papotti M, et al. Loss of miR-200c expression induces an aggressive, invasive, and chemoresistant phenotype in non-small cell lung cancer. Mol Cancer Res MCR. 2010;8:1207–16.
Neves R, Scheel C, Weinhold S, Honisch E, Iwaniuk KM, Trompeter H-I, et al. Role of DNA methylation in miR-200c/141 cluster silencing in invasive breast cancer cells. BMC Res Notes. 2010;3:219.
Bojmar L, Karlsson E, Ellegård S, Olsson H, Björnsson B, Hallböök O, et al. The role of microRNA-200 in progression of human colorectal and breast cancer. PloS One. 2013;8, e84815.
Chen J, Tian W, Cai H, He H, Deng Y. Down-regulation of microRNA-200c is associated with drug resistance in human breast cancer. Med Oncol Northwood Lond Engl. 2012;29:2527–34.
Tang H, Deng M, Tang Y, Xie X, Guo J, Kong Y, et al. miR-200b and miR-200c as prognostic factors and mediators of gastric cancer cell progression. Clin Cancer Res Off J Am Assoc Cancer Res. 2013;19:5602–12.
Liu Y-N, Yin JJ, Abou-Kheir W, Hynes PG, Casey OM, Fang L, et al. MiR-1 and miR-200 inhibit EMT via Slug-dependent and tumorigenesis via Slug-independent mechanisms. Oncogene. 2013;32:296–306.
Samavarchi-Tehrani P, Golipour A, David L, Sung H, Beyer TA, Datti A, et al. Functional genomics reveals a BMP-Driven mesenchymal-to-epithelial transition in the initiation of somatic cell reprogramming. Cell Stem Cell. 2010;7:64–77.
Shimono Y, Zabala M, Cho RW, Lobo N, Dalerba P, Qian D, et al. Downregulation of miRNA-200c links breast cancer stem cells with normal stem cells. Cell. 2009;138:592–603.
Gregory PA, Bert AG, Paterson EL, Barry SC, Tsykin A, Farshid G, et al. The miR-200 family and miR-205 regulate epithelial to mesenchymal transition by targeting ZEB1 and SIP1. Nat Cell Biol. 2008;10:593–601.
Gheldof A, Berx G. Cadherins and Epithelial-to-Mesenchymal Transition. Prog Mol Biol and Transl Sci. 2013;116:317–336.
Berx G, Staes K, van Hengel J, Molemans F, Bussemakers MJG, van Bokhoven A, et al. Cloning and characterization of the human invasion suppressor gene E-cadherin (CDH1). Genomics. 1995;26:281–9.
Dhawan D, Hamdy FC, Rehman I, Patterson J, Cross SS, Feeley KM, et al. Evidence for the early onset of aberrant promoter methylation in urothelial carcinoma. J Pathol. 2006;209:336–43.
Galván JA, Astudillo A, Vallina A, Crespo G, Folgueras MV, González MV. Prognostic and diagnostic value of epithelial to mesenchymal transition markers in pulmonary neuroendocrine tumors. BMC Cancer. 2014;14:855.
Borgna S, Armellin M, di Gennaro A, Maestro R, Santarosa M. Mesenchymal traits are selected along with stem features in breast cancer cells grown as mammospheres. Cell Cycle Georget Tex. 2012;11:4242–51.
Li L-C, Dahiya R. MethPrimer: designing primers for methylation PCRs. Bioinforma Oxf Engl. 2002;18:1427–31.
Jiang M, Zhang Y, Fei J, Chang X, Fan W, Qian X, et al. Rapid quantification of DNA methylation by measuring relative peak heights in direct bisulfite-PCR sequencing traces. Lab Investig J Tech Methods Pathol. 2010;90:282–90.
Kwiecien R, Kopp-Schneider A, Blettner M. Concordance analysis. Dtsch Ärztebl Int. 2011;108:515–21.
Bland JM, Altman DG. Statistical methods for assessing agreement between two methods of clinical measurement. Lancet. 1986;1:307–10.
Bethge N, Lothe RA, Honne H, Andresen K, Trøen G, Eknæs M, et al. Colorectal cancer DNA methylation marker panel validated with high performance in Non-Hodgkin lymphoma. Epigenetics Off J DNA Methylation Soc. 2014;9:428–36.
Ying J, Li H, Murray P, Gao Z, Chen Y-W, Wang Y, et al. Tumor-specific methylation of the 8p22 tumor suppressor gene DLC1 is an epigenetic biomarker for Hodgkin, nasal NK/T-cell and other types of lymphomas. Epigenetics. 2007;2:15–21.
Kanemoto M, Shirahata M, Nakauma A, Nakanishi K, Taniguchi K, Kukita Y, et al. Prognostic prediction of glioblastoma by quantitative assessment of the methylation status of the entire MGMT promoter region. BMC Cancer. 2014;14:641.
Feng X, Wang Z, Fillmore R, Xi Y. MiR-200, a new star miRNA in human cancer. Cancer Lett. 2014;344:166–73.
Davalos V, Moutinho C, Villanueva A, Boque R, Silva P, Carneiro F, et al. Dynamic epigenetic regulation of the microRNA-200 family mediates epithelial and mesenchymal transitions in human tumorigenesis. Oncogene. 2012;31:2062–74.
Carr IM, Valleley EMA, Cordery SF, Markham AF, Bonthron DT. Sequence analysis and editing for bisulphite genomic sequencing projects. Nucleic Acids Res. 2007;35, e79.
Wojdacz TK, Hansen LL, Dobrovic A. A new approach to primer design for the control of PCR bias in methylation studies. BMC Res Notes. 2008;1:54.
Berry D, Ben Mahfoudh K, Wagner M, Loy A. Barcoded primers used in multiplex amplicon pyrosequencing bias amplification. Appl Environ Microbiol. 2011;77:7846–9.
Vrba L, Jensen TJ, Garbe JC, Heimark RL, Cress AE, Dickinson S, et al. Role for DNA methylation in the regulation of miR-200c and miR-141 expression in normal and cancer cells. PLoS ONE. 2010;5, e8697.
This study was supported by grants from Italian Ministry of Health (J31J11000480001), Associazione Via di Natale and Associazione Italiana per la Ricerca sul Cancro (AIRC, MCO-10016).
The authors declare that they have no competing interest.
GB participated in the study design, carried out the methylation assays, performed data analysis and drafted the paper. AdG and VD performed expression analyses on cell lines and human tissues. MA collaborated with the methylation analyses and cell cultures. TP provided the breast cancer and normal tissues. RM participated in data interpretation and drafting of the paper. MS participated in the study design and data interpretation, supervised the study and the drafting of the paper. All authors read and approved the final amanuscript.
Roberta Maestro and Manuela Santarosa contributed equally to this work.
Additional file 1:
Primers used in the methylation analysis. (PDF 101 kb)
Additional file 2:
Amplification of miR-200/miR-141 locus performed with three couples of 5’-end tailed primers. (PDF 96 kb)
Additional file 3:
Validation of the plasmid DNA standards. (PDF 277 kb)
Additional file 4:
Amplification of miR-200c/miR-141 locus and methylation analysis. (PDF 3603 kb)
Additional file 5:
miR-200c/miR-141 locus methylation of MDA-MB-231 breast cancer cell line determined by BSP and NBSP performed with forward and reverse primers. (PDF 455 kb)
Additional file 6:
Inverse association between miR-200c/miR-141 locus methylation and miR-200c expression. (PDF 1483 kb)
Rights and permissions
Open Access This 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.
About this article
Cite this article
Brisotto, G., di Gennaro, A., Damiano, V. et al. An improved sequencing-based strategy to estimate locus-specific DNA methylation. BMC Cancer 15, 639 (2015). https://doi.org/10.1186/s12885-015-1646-6
- miR-200c/miR-141 locus
- Bisulfite treatment