Genome-wide DNA methylation profiling of non-small cell lung carcinomas
- Rejane Hughes Carvalho1, 2,
- Vanja Haberle3,
- Jun Hou1, 4,
- Teus van Gent1, 4,
- Supat Thongjuea3,
- Wilfred van IJcken5,
- Christel Kockx5,
- Rutger Brouwer5,
- Erikjan Rijkers5,
- Anieta Sieuwerts6,
- John Foekens6,
- Mirjam van Vroonhoven5, 7,
- Joachim Aerts8,
- Frank Grosveld1, 2, 4, 9,
- Boris Lenhard3Email author and
- Sjaak Philipsen1, 2, 4Email author
© Hughes Carvalho et al.; licensee BioMed Central Ltd. 2012
Received: 3 April 2012
Accepted: 22 June 2012
Published: 22 June 2012
Non-small cell lung carcinoma (NSCLC) is a complex malignancy that owing to its heterogeneity and poor prognosis poses many challenges to diagnosis, prognosis and patient treatment. DNA methylation is an important mechanism of epigenetic regulation involved in normal development and cancer. It is a very stable and specific modification and therefore in principle a very suitable marker for epigenetic phenotyping of tumors. Here we present a genome-wide DNA methylation analysis of NSCLC samples and paired lung tissues, where we combine MethylCap and next generation sequencing (MethylCap-seq) to provide comprehensive DNA methylation maps of the tumor and paired lung samples. The MethylCap-seq data were validated by bisulfite sequencing and methyl-specific polymerase chain reaction of selected regions.
Analysis of the MethylCap-seq data revealed a strong positive correlation between replicate experiments and between paired tumor/lung samples. We identified 57 differentially methylated regions (DMRs) present in all NSCLC tumors analyzed by MethylCap-seq. While hypomethylated DMRs did not correlate to any particular functional category of genes, the hypermethylated DMRs were strongly associated with genes encoding transcriptional regulators. Furthermore, subtelomeric regions and satellite repeats were hypomethylated in the NSCLC samples. We also identified DMRs that were specific to two of the major subtypes of NSCLC, adenocarcinomas and squamous cell carcinomas.
Collectively, we provide a resource containing genome-wide DNA methylation maps of NSCLC and their paired lung tissues, and comprehensive lists of known and novel DMRs and associated genes in NSCLC.
KeywordsDNA Methylation Epigenetics MethylCap Next generation sequencing Non-small cell lung Cancer
Non-small cell lung carcinoma (NSCLC) is a common malignancy characterized by a worldwide high incidence and low survival rate . NSCLC is a heterogenic disease which is broadly classified into three major histopathological subtypes: adenocarcinoma (ADC), squamous cell carcinoma (SCC) and large cell carcinoma (LCC). This heterogeneity poses challenges for diagnosis and treatment, since each subtype presents with a distinctive prognosis  and the choice of therapeutic regimen is predominantly based on tumor subtype and staging parameters . The development of personalized diagnostics and therapy is leading the way to a new era that may see us overcome some of the difficulties in treating complex diseases such as NSCLC.
In the past decade, comparative gene expression profiles of tumors have been extensively studied [4–6], yielding useful insights into the molecular hallmarks of carcinogenesis [7, 8]. With the advent of next generation sequencing, genome-wide screening has become an attractive tool for profiling tumors versus lung tissues [7, 9]. DNA methylation is a very stable epigenetic mark and next generation sequencing studies have recently shown that many genes are aberrantly methylated in cancer [10, 11]. Tissue specific DNA methylation patterns are stabilized during embryonic development, and faithfully maintained through cell divisions [12–14]. Nevertheless, established methylation patterns can be reprogrammed, with tumor cells undergoing DNA demethylation and de novo methylation through mechanisms not yet completely understood. CDKN2A and RASSF1 are examples of genes found to be aberrantly methylated in a wide variety of tumors [15–18], and epigenetic silencing of these genes has also been reported in NSCLC [19–22]. The great majority of DNA methylation studies are concentrated on the analysis of CpG islands located in the promoter areas of pre-selected genes. However, differentially methylated areas may be located within genes and at large distances from the nearest neighboring genes [23, 24]. Although data on methylated genes in NSCLC are rapidly accumulating, unbiased data concerning specificity of the genome-wide distribution of methylated loci are still scarce.
In this study, we used Methyl-DNA Capture (MethylCap) and high-throughput sequencing (MethylCap-seq, ) to perform a genome-wide DNA methylation screening of NSCLC tumors and paired adjacent lung tissues. With this approach, we sought to identify genome-wide aberrant methylation patterns of NSCLC. Specific differentially methylated regions would be promising candidate molecular markers for non-invasive diagnostics using circulating tumor DNA, and increase the number of possible targets for epigenetic therapy.
Methylation profiles in NSCLC-study outline
Data from patients used for MethylCap-seq and bisulfite sequencing validation
Overall survival (month)
Age at diagnosis
Global analysis of genome-wide methylation patterns of NSCLC
The DMRs were analyzed according to general hallmarks of their genomic localization (Figure 2E). Compared with the general distribution of randomly sampled regions across the genome (label “Random regions”), we found that promoter areas were particularly enriched among the DMRs, both hypo- (P=5.7 × 10-3) and hypermethylated (p=5.3 × 10-4). Although we removed most signals from repetitive regions by keeping only uniquely mapped reads, the signals from their uniquely mappable reads and flanking regions were still above background and allowed us to test for differential methylation. This revealed that some specific repeat classes such as LINE and satellites were also differentially represented in the DMRs, with satellite repeats being particularly hypomethylated in tumors when compared with lung tissues (P=1.3 × 10-5). By contrast, LINE repeats are relatively underrepresented in the DMRs (hypo and hypermethylation), indicating that the methylation status of the bulk of the LINE elements was similar between the tumors and paired lung counterparts.
Differentially methylated regions in NSCLC
We found 57 DMRs present in all tumor samples, of which 50 were hypermethylated and 7 were hypomethylated. An example of two DMRs in all samples, including the replicate experiment can be found in Additional file 10. We also analyzed DMRs in relation to NSCLC subtypes. We found 287 DMRs unique for SCC and 26 DMRs unique for ADC (Additioonal File 11). These DMRs were classified by hyper- and hypomethylation and their genomic locations relative to genes (Figure 3B). We found that the ADC-specific DMRs had a genomic distribution similar to that observed for the DMRs present in all seven tumors, while for SCCs, DMRs within genes were concentrated in the gene body. SCC also presented considerably more subgroup-specific DMRs than ADC and ‘All-sample’; this could be explained by the notion that SCC is a relatively homogeneous tumor type.
For the DMRs that fall in gene areas, we analyzed which types of genes were present in each of the three groups (Figure 3C and Additional file 11). We found that for the DMRs present in all tumor samples, more than 75% of the associated genes belonged to the class of transcriptional regulators. For SCC-specific DMRs, we observed a considerable heterogeneity in the functional categories of associated genes. Nevertheless, there was an overrepresentation of genes involved in transcriptional regulation, organization of the cytoskeleton and cell cycle regulation. We note that while some of the DMRs were associated with genes that had been previously reported to be differentially methylated in lung tumors, such as APC CDH13 CDKN2A DAPK hMLH1 HOX genes, OTX1 HOX2 and many others [22, 29–32], we also found DMRs which, to the best of our knowledge, have not been previously reported to be associated with methylation in NSCLC or other types of cancer.
Validation of differentially methylated regions in NSCLC by bisulfite sequencing
To assure sufficient depth of coverage for quantitative analysis, we sequenced 36 individual clones from each fragment and each bisulfite-converted DNA sample. Figure 4B shows an example of the methylation status for one fragment of a paired sample. The methylation status of the five fragments in all seven paired NSCLC/lung tissue samples is summarized in Figure 4C; full data for all samples is available in Additional files 12 and 13. We also examined potential differences in methylation according to tumor histology. To perform this analysis we grouped the samples in four categories: all clones from the lung tissues of patients (ADC_N and SCC_N) and clones from tumors (ADC_T and SCC_T), (Figure 4D and Additional file 14). We observed a statistically significant difference (P≤ 0.0001), not only between tumors and lung tissues, but also between ADC and SCC for all fragments analyzed by bisulfite sequencing. This identifies these fragments as candidate markers for differentiating tumors versus paired lung tissues, and also between tumor subtypes. Moreover this finding validates the MethylCap-seq data, since the methylation status determined by bisulfite sequencing correlated quantitatively with the MethylCap-seq data (Figure 4E).
Screening by methylation-specific PCR verifies methylation differences between tumor and paired lung tissues, and between tumor subtypes
Frequency of methylation by methylation-specific polymerase chain reaction in the 96 samples (48 tumor and matching lung tissue) in five fragments
Squamous cell carcinoma
MSP is a more qualitative approach than bisulfite sequencing. In MSP, selective amplification depends on the methylation status of only the CpGs present in the designed primers, which were two or three CpGs depending on the fragment. For each methylated (M) primer set, amplification of the fragment is dependent on methylation of the CpGs that are present in the sequence targeted by the “M” primer; while each unmethylated (U) primer set requires unmethylated CpGs in the referred sequence for amplification. In the 48 paired samples analyzed, CDKN2A showed frequent amplification with the M primer set in lung samples (24 out of 48), although the amplified bands were often of low intensity. Tumors had a much larger number displaying amplification with both M and U primer sets. Of the 47 tumors, 36 displayed amplification with the M primer set, and two of the tumors displayed amplification only with the M primer set. For RASSF1 and ZIC4, 41 paired samples were available. Both fragments differentiated between lung and tumor samples, whereas no methylation was detected in lung samples. ZIC4 showed methylation in more than 70% of the tumors, which is much higher than the 41% found for RASSF1. The findings for RASSF1 are in accordance with other studies . Also, Frag_01 presented no amplification with the M primer set in the lung samples, in the Frag_02 M primer set only 5 out of 48 were amplified. In tumors, we observed fragment amplification in 29out of 48 samples for Frag_01 and 36out of 48 samples for Frag_02 with the M primer sets. Differences in methylation status were significant in tumors versus lung samples for all 5 fragments, with a p-value < 0.02 for CDKN2A and P < 0.0001 for the other four fragments (Figure 5B). Frag_01, Frag_02 and ZIC4 present a much lower P-value, indicating they may be better tumor markers than CDKN2A and RASSF1. We also observed that significant differences between SCC and ADC were only found for Frag_01 (Figure 5C). In conclusion, the MSP results validate the MethylCap-seq data, since the results obtained for all five fragments were consistent with the significant DMRs in all tumors, and Frag_01 appeared as part of the SCC-specific DMRs in the MethylCap-seq data.
To date, genome-wide methylation studies in NSCLC have concentrated on specific areas of the genome , promoter regions [34–36] or sets of pre-selected candidate genes [19, 37]. Recently Kwon et al.  published a report on genome-wide analysis of DNA methylation in NSCLC. However, they also concentrate the analysis on a selection of candidate genes. In our study, we performed an unbiased genome-wide DNA methylation analysis of seven individual patients with NSCLC and their adjacent lung tissues. We combined MethylCap with high-throughput sequencing (MethylCap-seq) to draw detailed methylome maps of NSCLC tumors and paired lung tissues. The replicate experiments showed that the depth of coverage generated by MethylCap-seq was sufficient to capture the methylation status of the entire genome. Moreover, there was a very high degree of correlation between the replicate experiments, demonstrating the robustness of the MethylCap-seq approach. As controls for normalization and validation purposes we made use of fully methylated and fully unmethylated DNA samples. The MethylCap-seq data of these controls provide a useful resource for benchmarking of the fully methylated versus fully unmethylated status of any particular fragment in the genome. Using the MEDIPS computational tool , we generated methylation profiles of the seven paired NSCLC/lung cases which can be viewed in the genome browsers. This broad view of DNA methylation in NSCLC will provide new opportunities for the identification of specific epigenetic markers of NSCLC that could, for instance, be used for early detection of the disease.
The occurrence of methylated CpGs in the fragments selected for bisulfite sequencing showed a high quantitative correlation with the reads found by MethylCap-seq in the same area. In conjunction with the MethylCap-seq data on the fully methylated and fully unmethylated controls, this demonstrates that MethylCap-seq yields reliable quantitative information on the methylation status of the fragments in a particular DNA sample. This is important considering that the tumor cell content of samples to be analyzed will be variable. Provided that sufficient sequencing depth is achieved, it should be possible to use MethylCap-seq for the identification of tumor-specific hypermethylated DNA regions in samples containing only a minority of tumor cells, and possibly in circulating tumor DNA isolated from serum samples. The results obtained with MSP are qualitative in nature, but this PCR-based approach is very sensitive and may therefore be used to detect tumor-specific hypermethylated DNA fragments in serum DNA. As a first step in this direction, we used MSP to investigate the presence of five DMRs in a larger set of paired samples. The results confirmed the MethylCap-seq data and demonstrated that these DMRs are frequently hypermethylated in NSCLC. In particular, the Frag_01 fragment was found to be hypermethylated in 13 out of 16 SCC cases and in none of the ADC, suggesting that this might be a useful marker for this type of NSCLC, and ZIC4 showed 100% specificity for tumors and higher sensitivity (>70%) than RASSF1. In addition, CDKN2A amplification with the M primer was often observed in lung tissue samples, although at much lower amplification levels when compared with tumor samples. This could reflect an early event in the methylation status of CDKN2A in tumorigenesis.
The mapped reads point to differentially methylated areas across the genome, with some marked preferences for hypermethylation or hypomethylation in certain genomic regions. There is a strong correlation between CpG islands and hypermethylation in tumors, with the most significantly hypermethylated regions associated with promoter regions of genes. As previously reported by many groups [40–43], gene promoters are a target of methylation as an epigenetic regulatory mechanism. However, our data also reveal a high degree of hypermethylation outside promoter areas at intergenic regions, and at long distances from genes. This indicates that patterns of DNA methylation may play a role not only by silencing promoters of key tumor suppressor genes, but also by regulating gene expression in a more complex manner through distant regulatory elements such as insulators and enhancers.
We found that satellite regions were particularly hypomethylated in the NSCLC tumors when compared with lung tissues. Recently, Ting et al.  have shown that overexpression of satellite repeats is associated with different types of cancer. The observed hypomethylation of these regions in NSCLC would be compatible with their results. We also found that hypomethylated regions are often present in subtelomeric regions, in agreement with the observation that hypomethylation is more frequently located near chromosome ends in lung cancer . While CpG islands and genes are present in the subtelomeric regions of chromosomes, the impact of subtelomeric hypomethylation on the regulation of these genes is currently unclear.
In an alternative approach, we selected the most significant DMRs in the individual NSCLC cases (14,742) and performed unsupervised cluster analysis of the hyper- and hypomethylated regions in the tumor samples. We observed that the samples clustered according to their histological classification. This indicates that there are DMRs that could not only be used as NSCLC markers but also as markers for histological classification of NSCLC tumors. Subsequently we searched for the specific areas that were differentially methylated in all tumors when compared with the paired lung tissues, and specific for the subtypes. We observed that SCCs displayed more DMRs than did ADCs did. This observation could be explained by the heterogeneity of ADCs, which would result in fewer DMRs shared by all three ADC samples. The contrary is seen in SCC, which is a more homogeneous type of NSCLC. This characteristic has been observed previously in gene expression profiling studies [4, 46].
Remarkably, the majority (>75%) of the hypermethylated genes present in all seven samples encoded for transcription factors, while the shared hypomethylated genes were spread among many different functional categories. Transcription factors play a central role in maintaining or modifying cell fate, either in normal processes such as development and cell differentiation, or in cancer. Mechanisms of epigenetic regulation via DNA methylation are not yet completely understood; several factors such as DNA methyltransferases, chromatin remodelling proteins and DNA-binding transcription factors are involved [47, 48]. Therefore, any perturbation that leads to a decrease in the expression status of those factors may disrupt important pathways for maintenance of the differentiated cellular state. Aberrant methylation can change chromatin structure, rendering DNA binding sites accessible or inaccessible to transcription factors leading to activation or silencing of genes important in cell differentiation and maintenance. It has been recently shown that hypermethylated genes found in lung tumors are associated with cellular morphogenetic differentiation . These cellular mechanisms are orchestrated by transcriptional regulators; therefore changes in the methylation status of their binding sites may disrupt those processes thus contributing to oncogenesis.
List of most differentially methylated genes in non-small cell carcinoma
Gene class/type/ biological function
Common in all tumors
BARX homeobox 1
paired box 9
orthodenticle homeobox 1
G-protein coupled receptor
natriuretic peptide receptor C/guanylate cyclase C (atrionatriuretic peptide receptor C)
fibroblast growth factor 12
one cut homeobox 2
PR domain containing 14
retina and anterior neural fold homeobox
short stature homeobox 2
DMRT-like family A2
fer-1-like 4 (Caenorhabditis. elegans) pseudogene
SIX homeobox 6
GATA binding protein 3
SKI familytranscriptional corepressor 1
sal-like 1 (Drosophila)
iroquois homeobox 2
glutamate receptor, ionotropic, kainate 2
SATB homeobox 2
Meis homeobox 1
ventral anterior homeobox 1
cyclin N-terminal domain containing 2
zinc finger, MYND-type containing 10
v-myc myelocytomatosis viral oncogene homolog (avian)
neuron navigator 1
cytoplasmic polyadenylation element binding protein 3
Unique for squamous cell carcinomas
trafficking protein particle complex 9
abhydrolase domain containing 2
catenin (cadherin-associated protein), delta 1
histone cluster 1, H2bb
epithelial membrane protein 1
transducin (beta)-like 1 X-linked receptor 1
Zic family member 4
acyloxyacyl hydrolase (neutrophil)
actinin, alpha 4
chromosome 1 open reading frame 21
protein kinase C and casein kinase substrate in neurons 2
DOT1-like, histone H3 methyltransferase (Sacchormyces. cerevisiae)
e3 ubiquitin-protein ligase
WW domain containing E3 ubiquitin protein ligase 2
generaltranscription factor IIIC, polypeptide 1, alpha 220 kDa
MDN1, midasin homolog (yeast)
death inducer-obliterator 1
histone cluster 1, H3c
ankyrin repeat domain 13B
calcitonin-related polypeptide beta
protein tyrosine phosphatase, receptor type, A
signal transducer and activator of transcription 5A
LIM domain kinase 1
solute carrier family 23 (nucleobase transporters), member 2
BARX homeobox 1
Nance-Horan syndrome (congenital cataracts and dental anomalies)
forkhead box K1
protein-L-isoaspartate (D-aspartate) O-methyltransferase
SET domain containing 1A
centromere protein P
enhances neuronal dendrite outgrowth
SLIT and NTRK-like family, member 1
RAR-related orphan receptor B
postsynaptic scaffold in neuronal cells
discs, large (Drosophila) homolog-associated protein 1
chromosome 3 open reading frame 21
ST6 (alpha-N-acetyl-neuraminyl-2,3-beta-galactosyl-1,3)-N-acetylgalactosaminide alpha-2,6-sialyltransferase 1
Zic family member 3 (odd-paired homolog, Drosophila)
transmembrane protein 132D
microtubule-associated protein 1 light chain 3 beta 2
chibby homolog 3 (Drosophila)
Ras interacting protein 1
protein kinase, cGMP-dependent, type I
WD repeat domain 72
potassium voltage-gated channel, KQT-like subfamily, member 2
G-protein coupled receptor
BAI1-associated protein 3
mitogen-activated protein kinase kinase kinase 10
Unique for Adenocarcinomas
family with sequence similarity 78, member B
growth arrest-specific 1
protein tyrosine phosphatase, receptor type, N polypeptide 2
Wnt receptor signaling pathway
R-spondin 2 homolog (Xenopus laevis)
POU class 3 homeobox 3
transient receptor potential cation channel, subfamily A, member 1
solute carrier family 6 (neurotransmitter transporter, noradrenalin), member 2
LIM homeobox 1
small GTPase mediated signal transduction
Rap guanine nucleotide exchange factor (GEF) 5
growth differentiation factor 10
chromosome 3 open reading frame 45
slit homolog 2 (Drosophila)
zinc finger protein 423
actin filament organization
ras homolog gene family, member F (in filopodia)
Samples were obtained from the patients with NSCLC (n = 48) who had undergone surgical lung resection at Erasmus University Medical Center Rotterdam. Specimens were collected from the tumor and adjacent non-cancerous lung tissue and studied under an anonymous tissue protocol approved by the medical ethical committee of Erasmus University Medical Center Rotterdam. Tissues were snap-frozen within two hours after surgical resection in liquid nitrogen pre-cooled isopentane, and stored at 196 °C or 80 °C until DNA extraction.
Patient samples were independently reviewed by two pathologists. The cohort included 17 with ADC; 16 with SCC, seven with LCC , eight unclassified samples, and paired lung tissues for each tumor sample.
The MRC-5 lung fibroblast-like cell line was used as control. Cells were cultured under standard conditions using minimum essential medium supplemented with 10% heat-inactivated fetal bovine serum, 2 mM L-glutamine, 1% non-essential amino acids and penicillin/streptomycin. Cells were harvested when they reached 90% confluence.
Genomic DNA from patient tissues and cultured cells were extracted by overnight treatment incubation with lysis buffer and proteinase K, followed by phenol-chloroform extraction, ethanol precipitation and RNase digestion.
Artificial demethylation and methylation of genomic DNA
We used DNA extracted from the MRC-5 cell line and commercially available Universal unmethylated DNA (UUD; Millipore(Billerica, MA, USA). Fully unmethylated DNA was obtained by whole-genome amplification using the REPLI-g kit (QiagenQiagen – Germantown, MD, USA) according to manufacturer protocol, followed by phenol-chloroform extraction. Fully methylated DNA was prepared by treating MRC-5 DNA and UUD with M.SssI enzyme (New England Biolabs - Ipswich, MA, USA) according to the manufacturer’s protocol. In short, 10 μg of DNA was incubated for 2 hours at 37°C with 40 U of M.SssI and 640 μM of S-adenosylmethionine. DNA was then treated with phenol-chloroform and recovered by ethanol precipitation.
Enrichment of methylated DNA was carried out using the MethylCap technique. DNA samples were sheared by sonication to obtain fragments between 200 and 800 bp. Methylated DNA capture was carried out using an adapted protocol from the MethylMiner Methylated DNA Enrichment kit (Invitrogen – Carlsbad, CA, USA). Two adaptations were made. Firstly, 20 μL of beads and 14 μL (7 μg) of MBD-biotin protein were used for 5 μg of sheared genomic DNA. Secondly, a single elution fraction was obtained by resuspending the beads in 200 μL of 1X binding/washing buffer containing 2 μL of proteinase K (20 mg/ml). The samples were then incubated for 90 minutes at 57 °C with 800 rpm agitation. Remaining steps were performed following the manufacturer’s protocol scaled up to 5 μg of DNA. Before the samples were used for high-throughput sequencing, methyl-enrichment was tested by quantitative (q) PCR (data not shown). qPCR was carried out using SYBR r green, Phire Hot Start DNA polymerase (Finnzymes - Vantaa, Finland)), and 1 μL of the enriched samples. Primer sequences are available in Additional file 16. qPCR was performed in triplicate and enrichment levels were calculated as previously described .
High-throughput sequencing (HTS) was carried out on 18 samples: seven tumors (three ADC, three SCC and one LCC), seven paired lung tissues and four controls (two artificially fully methylated DNA and two fully unmethylated DNA). For each sample, DNA recovered from two independent Methyl Capture experiments was sequenced. The Illumina Genome Analyzer IIx platform (San Diego, CA, USA was used for both replicates, one sample per lane, according to the manufacturer’s protocol. In short: after MethylCap, fragments of the methyl-enriched fraction were end-repaired and ligated to single-read adaptors. Samples were size-selected to an average of 320 bp, PCR-amplified and 36 bp sequenced.
The Illumina Casava pipeline was used for base calling, alignment and quality control. Results were mapped against the reference Human_UCSChg18_AllChromosomes using eland_extended by Illumina pipeline 1.6.0.
Detailed information about the evaluation and selection of normalization parameters, and identification, characterization and selection of DMRs can be found in Additional file 3.
Bisulfite sequencing was carried out on the same 16 samples used for MethylCap HTS. Detailed information is provided in Additional file 3. In short: bisulfite conversion was carried out using the Epitect kit from Qiagen following the manufacturer’s protocol. Samples were amplified using primers specific for the fragments of interest ( Additional file 16). Amplified fragments were loaded on 2% agarose gels and extracted from the gel using a NucleoSpin Extract II kit (Machery-Nagel – Düren, Germany). Fragments were then ligated to pGEM-t easy vector (Promega - Madison, WI, USA) and cloned into DH12-β competent cells. Colony PCR was performed on 36 colonies and PCR fragments sequenced. Analysis of methylated and unmethylated CpGs was executed using two online programs: BISMA (http://biochem.jacobs-university.de/BDPC/BISMA) and QUMA (http://quma.cdb.riken.jp).
To determining the methylation status of the selected regions, we used MSP  in 96 samples (48 tumors and matching lung tissues). Samples were bisulfite converted as described for bisulfite sequencing, and amplified using two different set of primers designed for the methylated and unmethylated sequences ( Additional file 16). The fully methylated and fully unmethylated DNA samples were used as controls, and a water blank reaction used as control for contamination. After amplification, products were resolved on 2% agarose gels containing ethidium bromide and visualized under UV transillumination.
Statistical analysis for bisulfite sequencing was provided by the BISMA program. In short, Fisher’s exact test was used for the statistical significance between two bisulfite sequence groups at each CpG, and a Mann-Whitney U-test for the statistical significance between two groups of the entire set of CpG sites. Two-tailed P-value of Fisher’s exact test was calculated from the 2 × 2 tables at each CpG site. A Mann-Whitney U-test was used for the statistical significance of the entire set of CpG sites between the tumor and lung groups. For determination of significance in the MSP, we used chi-square distribution. Two-tailed P-values were determined by the counts of the amplified bands for each primer in the group of samples lung/tumor or SCC/ADC. Groups that showed a P-value < 0.05 were considered significantly different.
differentially methylated region
large cell carcinoma
long interspersed element
non-small cell lung carcinoma
polymerase chain reaction
squamous cell carcinoma
Universal unmethylated DNA
We thank Peter Riegman and Monique Oomen for providing patient samples stored in the Erasmus MC tissue bank, and Ilya Grigorev for the MRC-5 cell line. We thank Nynke Gillemans for technical assistance and Thamar van Dijk, Cíntia Bombardieri, Ricardo Leite, Vojka Gouka and Marco Gouka for helpful discussions during the preparation of this manuscript. We also thank Alex Pindyurin for this help with the R script. This work was supported by the Netherlands Genomics Initiative.
- Jemal A, Siegel R, Xu J, Ward E: Cancer statistics, 2010. CA Cancer J Clin. 2010, 60: 277-300.View ArticlePubMedGoogle Scholar
- Little AG: No nodes is good nodes. Ann Thorac Surg. 2006, 82: 4-5.View ArticlePubMedGoogle Scholar
- Molina JR, Yang P, Cassivi SD, Schild SE, Adjei AA: Non-small cell lung cancer: epidemiology, risk factors, treatment, and survivorship. Mayo Clin Proc. 2008, 83: 584-594.PubMed CentralView ArticlePubMedGoogle Scholar
- Hou J, Aerts J, den Hamer B, van Ijcken W, den Bakker M, Riegman P, van der Leest C, van der Spek P, Foekens JA, Hoogsteden HC, Grosveld F, Philipsen S: Gene expression-based classification of non-small cell lung carcinomas and survival prediction. PLoS One. 2010, 5: e10312.PubMed CentralView ArticlePubMedGoogle Scholar
- Larsen JE, Pavey SJ, Passmore LH, Bowman R, Clarke BE, Hayward NK, Fong KM: Expression profiling defines a recurrence signature in lung squamous cell carcinoma. Carcinogenesis. 2007, 28: 760-766.View ArticlePubMedGoogle Scholar
- Liang P, Pardee AB: Analysing differential gene expression in cancer. Nat Rev Cancer. 2003, 3: 869-876.View ArticlePubMedGoogle Scholar
- Kameyama K, Huang CL, Liu D, Masuya D, Nakashima T, Sumitomo S, Takami Y, Kinoshita M, Yokomise H: Reduced ING1b gene expression plays an important role in carcinogenesis of non-small cell lung cancer patients. Clin Cancer Res. 2003, 9: 4926-4934.PubMedGoogle Scholar
- Petty RD, Nicolson MC, Kerr KM, Collie-Duguid E, Murray GI: Gene expression profiling in non-small cell lung cancer: from molecular mechanisms to clinical application. Clin Cancer Res. 2004, 10: 3237-3248.View ArticlePubMedGoogle Scholar
- Ding L, Mardis ER, Li S, Larson DE, Chen K, Wallis JW, Harris CC, McLellan MD, Fulton RS, Fulton LL, Abbott RM, Hoog J, Dooling DJ, Koboldt DC, Schmidt H, Kalicki J, Zhang Q, Chen L, Lin L, Wendl MC, McMichael JF, Magrini VJ, Cook L, McGrath SD, Vickery TL, Appelbaum E, Deschryver K, Davies S, Guintoli T, Lin L: Genome remodelling in a basal-like breast cancer metastasis and xenograft. Nature. 2010, 464: 999-1005.PubMed CentralView ArticlePubMedGoogle Scholar
- Lister R, Pelizzola M, Dowen RH, Hawkins RD, Hon G, Tonti-Filippini J, Nery JR, Lee L, Ye Z, Ngo QM, Edsall L, Antosiewicz-Bourget J, Stewart R, Ruotti V, Millar AH, Thomson JA, Ren B, Ecker JR: Human DNA methylomes at base resolution show widespread epigenomic differences. Nature. 2009, 462: 315-322.PubMed CentralView ArticlePubMedGoogle Scholar
- Weber M, Davies JJ, Wittig D, Oakeley EJ, Haase M, Lam WL, Schubeler D: Chromosome-wide and promoter-specific analyses identify sites of differential DNA methylation in normal and transformed human cells. Nat Genet. 2005, 37: 853-862.View ArticlePubMedGoogle Scholar
- Bird A: DNA methylation patterns and epigenetic memory. Genes Dev. 2002, 16: 6-21.View ArticlePubMedGoogle Scholar
- Bird AP, Wolffe AP: Methylation-induced repression–belts, braces, and chromatin. Cell. 1999, 99: 451-454.View ArticlePubMedGoogle Scholar
- Leonhardt H, Page AW, Weier HU, Bestor TH: A targeting sequence directs DNA methyltransferase to sites of DNA replication in mammalian nuclei. Cell. 1992, 71: 865-873.View ArticlePubMedGoogle Scholar
- Agathanggelou A, Honorio S, Macartney DP, Martinez A, Dallol A, Rader J, Fullwood P, Chauhan A, Walker R, Shaw JA, Hosoe S, Lerman MI, Minna JD, Maher ER, Latif F: Methylation associated inactivation of RASSF1A from region 3p21.3 in lung, breast and ovarian tumours. Oncogene. 2001, 20: 1509-1518.View ArticlePubMedGoogle Scholar
- Kohonen-Corish MR, Sigglekow ND, Susanto J, Chapuis PH, Bokey EL, Dent OF, Chan C, Lin BP, Seng TJ, Laird PW, Young J, Leggett BA, Jass JR, Sutherland RL: Promoter methylation of the mutated in colorectal cancer gene is a frequent early event in colorectal cancer. Oncogene. 2007, 26: 4435-4441.View ArticlePubMedGoogle Scholar
- Liu L, Broaddus RR, Yao JC, Xie S, White JA, Wu TT, Hamilton SR, Rashid A: Epigenetic alterations in neuroendocrine tumors: methylation of RAS-association domain family 1, isoform A and p16 genes are associated with metastasis. Mod Pathol. 2005, 18: 1632-1640.PubMedGoogle Scholar
- Pfeifer GP, Yoon JH, Liu L, Tommasi S, Wilczynski SP, Dammann R: Methylation of the RASSF1A gene in human cancers. Biol Chem. 2002, 383: 907-914.View ArticlePubMedGoogle Scholar
- Feng Q, Hawes SE, Stern JE, Wiens L, Lu H, Dong ZM, Jordan CD, Kiviat NB, Vesselle H: DNA methylation in tumor and matched normal tissues from non-small cell lung cancer patients. Cancer Epidemiol Biomarkers Prev. 2008, 17: 645-654.PubMed CentralView ArticlePubMedGoogle Scholar
- Hsu HS, Chen TP, Hung CH, Wen CK, Lin RK, Lee HC, Wang YC: Characterization of a multiple epigenetic marker panel for lung cancer detection and risk assessment in plasma. Cancer. 2007, 110: 2019-2026.View ArticlePubMedGoogle Scholar
- Vaissiere T, Hung RJ, Zaridze D, Moukeria A, Cuenin C, Fasolo V, Ferro G, Paliwal A, Hainaut P, Brennan P, Tost J, Boffetta P, Herceg Z: Quantitative analysis of DNA methylation profiles in lung cancer identifies aberrant DNA methylation of specific genes and its association with gender and cancer risk factors. Cancer Res. 2009, 69: 243-252.PubMed CentralView ArticlePubMedGoogle Scholar
- Dammann R, Li C, Yoon JH, Chin PL, Bates S, Pfeifer GP: Epigenetic inactivation of a RAS association domain family protein from the lung tumour suppressor locus 3p21.3. Nat Genet. 2000, 25: 315-319.View ArticlePubMedGoogle Scholar
- Clark SJ: Action at a distance: epigenetic silencing of large chromosomal regions in carcinogenesis. Hum Mol Genet. 2007, 16 Spec No 1: R88-R95.View ArticlePubMedGoogle Scholar
- Mayor R, Casadome L, Azuara D, Moreno V, Clark SJ, Capella G, Peinado MA: Long-range epigenetic silencing at 2q14.2 affects most human colorectal cancers and may have application as a non-invasive biomarker of disease. Br J Cancer. 2009, 100: 1534-1539.PubMed CentralView ArticlePubMedGoogle Scholar
- Brinkman AB, Simmer F, Ma K, Kaan A, Zhu J, Stunnenberg HG: Whole-genome DNA methylation profiling using MethylCap-seq. Methods. 2010, 52: 232-236.View ArticlePubMedGoogle Scholar
- Teodoridis JM, Hardie C, Brown R: CpG island methylator phenotype (CIMP) in cancer: causes and implications. Cancer Lett. 2008, 268: 177-186.View ArticlePubMedGoogle Scholar
- Keshet I, Schlesinger Y, Farkash S, Rand E, Hecht M, Segal E, Pikarski E, Young RA, Niveleau A, Cedar H, Simon I: Evidence for an instructive mechanism of de novo methylation in cancer cells. Nat Genet. 2006, 38: 149-153.View ArticlePubMedGoogle Scholar
- Rauch TA, Pfeifer GP: DNA methylation profiling using the methylated-CpG island recovery assay (MIRA). Methods. 2010, 52: 213-217.PubMed CentralView ArticlePubMedGoogle Scholar
- Rauch T, Wang Z, Zhang X, Zhong X, Wu X, Lau SK, Kernstine KH, Riggs AD, Pfeifer GP: Homeobox gene methylation in lung cancer studied by genome-wide analysis with a microarray-based methylated CpG island recovery assay. Proc Natl Acad Sci U S A. 2007, 104: 5527-5532.PubMed CentralView ArticlePubMedGoogle Scholar
- Rauch TA, Zhong X, Wu X, Wang M, Kernstine KH, Wang Z, Riggs AD, Pfeifer GP: High-resolution mapping of DNA hypermethylation and hypomethylation in lung cancer. Proc Natl Acad Sci U S A. 2008, 105: 252-257.PubMed CentralView ArticlePubMedGoogle Scholar
- Schmidt B, Liebenberg V, Dietrich D, Schlegel T, Kneip C, Seegebarth A, Flemming N, Seemann S, Distler J, Lewin J, Tetzner R, Weickmann S, Wille U, Liloglou T, Raji O, Walshaw M, Fleischhacker M, Witt C, Field JK: SHOX2 DNA methylation is a biomarker for the diagnosis of lung cancer based on bronchial aspirates. BMC Cancer. 2010, 10: 600.PubMed CentralView ArticlePubMedGoogle Scholar
- Zochbauer-Muller S, Fong KM, Virmani AK, Geradts J, Gazdar AF, Minna JD: Aberrant promoter methylation of multiple genes in non-small cell lung cancers. Cancer Res. 2001, 61: 249-255.PubMedGoogle Scholar
- Herman JG, Graff JR, Myohanen S, Nelkin BD, Baylin SB: Methylation-specific PCR: a novel PCR assay for methylation status of CpG islands. Proc Natl Acad Sci U S A. 1996, 93: 9821-9826.PubMed CentralView ArticlePubMedGoogle Scholar
- Bibikova M, Lin Z, Zhou L, Chudin E, Garcia EW, Wu B, Doucet D, Thomas NJ, Wang Y, Vollmer E, Goldmann T, Seifart C, Jiang W, Barker DL, Chee MS, Floros J, Fan JB: High-throughput DNA methylation profiling using universal bead arrays. Genome Res. 2006, 16: 383-393.PubMed CentralView ArticlePubMedGoogle Scholar
- Sano A, Kage H, Sugimoto K, Kitagawa H, Aki N, Goto A, Fukayama M, Nakajima J, Takamoto S, Nagase T, Yatomi Y, Ohishi N, Takai D: A second-generation profiling system for quantitative methylation analysis of multiple gene promoters: application to lung cancer. Oncogene. 2007, 26: 6518-6525.View ArticlePubMedGoogle Scholar
- Shames DS, Girard L, Gao B, Sato M, Lewis CM, Shivapurkar N, Jiang A, Perou CM, Kim YH, Pollack JR, Fong KM, Lam CL, Wong M, Shyr Y, Nanda R, Olopade OI, Gerald W, Euhus DM, Shay JW, Gazdar AF, Minna JD: A genome-wide screen for promoter methylation in lung cancer identifies novel methylation markers for multiple malignancies. PLoS Med. 2006, 3: 486.View ArticleGoogle Scholar
- Zhang Y, Wang R, Song H, Huang G, Yi J, Zheng Y, Wang J, Chen L: Methylation of multiple genes as a candidate biomarker in non-small cell lung cancer. Cancer Letters. 2011, 303: 21-28.View ArticlePubMedGoogle Scholar
- Kwon YJ, Lee SJ, Koh JS, Kim SH, Lee HW, Kang MC, Bae JB, Kim YJ, Park JH: Genome-wide analysis of DNA methylation and the gene expression change in lung cancer. J Thorac Oncol. 2012, 7: 20-33.View ArticlePubMedGoogle Scholar
- Chavez L, Jozefczuk J, Grimm C, Dietrich J, Timmermann B, Lehrach H, Herwig R, Adjaye J: Computational analysis of genome-wide DNA methylation during the differentiation of human embryonic stem cells along the endodermal lineage. Genome Res. 2010, 20: 1441-1450.PubMed CentralView ArticlePubMedGoogle Scholar
- Chung JH, Lee HJ, Kim BH, Cho NY, Kang GH: DNA methylation profile during multistage progression of pulmonary adenocarcinomas. Virchows Arch. 2011, 459: 201-211.View ArticlePubMedGoogle Scholar
- Esteller M: Epigenetic gene silencing in cancer: the DNA hypermethylome. Hum Mol Genet. 2007, 16 Spec No 1: R50-R59.View ArticlePubMedGoogle Scholar
- Liu Y, Pang JC, Dong S, Mao B, Poon WS, Ng HK: Aberrant CpG island hypermethylation profile is associated with atypical and anaplastic meningiomas. Hum Pathol. 2005, 36: 416-425.View ArticlePubMedGoogle Scholar
- Marsit CJ, Houseman EA, Christensen BC, Eddy K, Bueno R, Sugarbaker DJ, Nelson HH, Karagas MR, Kelsey KT: Examination of a CpG island methylator phenotype and implications of methylation profiles in solid tumors. Cancer Res. 2006, 66: 10621-10629.View ArticlePubMedGoogle Scholar
- Ting DT, Lipson D, Paul S, Brannigan BW, Akhavanfard S, Coffman EJ, Contino G, Deshpande V, Iafrate AJ, Letovsky S, Rivera MN, Bardeesy N, Maheswaran S, Haber DA: Aberrant overexpression of satellite repeats in pancreatic and other epithelial cancers. Science. 2011, 331: 593-596.PubMed CentralView ArticlePubMedGoogle Scholar
- Pfeifer GP, Rauch TA: DNA methylation patterns in lung carcinomas. Semin Cancer Biol. 2009, 19: 181-187.PubMed CentralView ArticlePubMedGoogle Scholar
- Hou J, Lambers M, den Hamer B, den Bakker M, Hoogsteden MC, Grosveld F, Hegmans J, Aerts J, Philipsen S: Expression profiling-based subtyping identifies novel non-small cell lung cancer subgroups and implicates putative resistance to pemetrexed therapy. J Thorac Oncol. 2012, 7: 105-114.View ArticlePubMedGoogle Scholar
- Chatagnon A, Ballestar E, Esteller M, Dante R: A role for methyl-CpG binding domain protein 2 in the modulation of the estrogen response of pS2/TFF1 gene. PLoS One. 2010, 5: e9665.PubMed CentralView ArticlePubMedGoogle Scholar
- Clouaire T, de Las Heras JI, Merusi C, Stancheva I: Recruitment of MBD1 to target genes requires sequence-specific interaction of the MBD domain with methylated DNA. Nucleic Acids Res. 2010, 38: 4620-4634.PubMed CentralView ArticlePubMedGoogle Scholar
- Helman E, Naxerova K, Kohane IS: DNA hypermethylation in lung cancer is targeted at differentiation-associated genes. Oncogene. 2011, 31: 1181-1188.View ArticlePubMedGoogle Scholar
- Karlen Y, McNair A, Perseguers S, Mazza C, Mermod N: Statistical significance of quantitative PCR. BMC Bioinformatics. 2007, 8: 131.PubMed CentralView ArticlePubMedGoogle Scholar
- R Development Core Team: R: A language and environment for statistical computing. Reference index version: 2.11.1. 2010, R Foundation for Statistical Computing, Vienna, Austria,http://www.R-project.org.Google Scholar
- Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5: R80.PubMed CentralView ArticlePubMedGoogle Scholar
- Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, Zahler AM, Haussler D: The human genome browser at UCSC. Genome Res. 2002, 12: 996-1006.PubMed CentralView ArticlePubMedGoogle Scholar
- Down TA, Rakyan VK, Turner DJ, Flicek P, Li H, Kulesha E, Graf S, Johnson N, Herrero J, Tomazou EM, Thorne NP, Bäckdahl L, Herberth M, Howe KL, Jackson DK, Miretti MM, Marioni JC, Birney E, Hubbard TJ, Durbin R, Tavaré S, Beck S: A Bayesian deconvolution strategy for immunoprecipitation-based DNA methylome analysis. Nat Biotechnol. 2008, 26: 779-785.PubMed CentralView ArticlePubMedGoogle Scholar
- Pelizzola M, Koga Y, Urban AE, Krauthammer M, Weissman S, Halaban R, Molinaro AM: MEDME: an experimental and analytical methodology for the estimation of DNA methylation levels based on microarray derived MeDIP-enrichment. Genome Res. 2008, 18: 1652-1659.PubMed CentralView ArticlePubMedGoogle 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/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.