High-frequency aberrantly methylated targets in pancreatic adenocarcinoma identified via global DNA methylation analysis using methylCap-seq
Clinical Epigeneticsvolume 6, Article number: 18 (2014)
Extensive reprogramming and dysregulation of DNA methylation is an important characteristic of pancreatic cancer (PC). Our study aimed to characterize the genomic methylation patterns in various genomic contexts of PC. The methyl capture sequencing (methylCap-seq) method was used to map differently methylated regions (DMRs) in pooled samples from ten PC tissues and ten adjacent non-tumor (PN) tissues. A selection of DMRs was validated in an independent set of PC and PN samples using methylation-specific PCR (MSP), bisulfite sequencing PCR (BSP), and methylation sensitive restriction enzyme-based qPCR (MSRE-qPCR). The mRNA and expressed sequence tag (EST) expression of the corresponding genes was investigated using RT-qPCR.
A total of 1,131 PC-specific and 727 PN-specific hypermethylated DMRs were identified in association with CpG islands (CGIs), including gene-associated CGIs and orphan CGIs; 2,955 PC-specific and 2,386 PN-specific hypermethylated DMRs were associated with gene promoters, including promoters containing or lacking CGIs. Moreover, 1,744 PC-specific and 1,488 PN-specific hypermethylated DMRs were found to be associated with CGIs or CGI shores. These results suggested that aberrant hypermethylation in PC typically occurs in regions surrounding the transcription start site (TSS). The BSP, MSP, MSRE-qPCR, and RT-qPCR data indicated that the aberrant DNA methylation in PC tissue and in PC cell lines was associated with gene (or corresponding EST) expression.
Our study characterized the genome-wide DNA methylation patterns in PC and identified DMRs that were distributed among various genomic contexts that might influence the expression of corresponding genes or transcripts to promote PC. These DMRs might serve as diagnostic biomarkers or therapeutic targets for PC.
Pancreatic cancer (PC), a highly malignant tumor of the digestive system, is a type of solid tumor that currently has one of the worst prognoses, with a postoperative 5-year survival rate of less than 25%. Nearly 100,000 people die from PC every year in the USA and Europe; PC ranks fourth and fifth for cancer mortality in those countries, respectively . The incidence of PC in China is also displaying an increasing trend .
From a biological perspective, genetics is predominantly responsible for the stable transfer of hereditary information between generations, whereas the distinct somatic phenotypes in different tissues and cells are influenced by epigenetics. Because most tumors that develop display particular acquired biological phenotypes, epigenetic changes must surely play important roles during tumor development . DNA methylation, a well-studied epigenetic phenomenon, has already been extensively studied in PC. At the gene particular level, genes such as p14ARF and p16INK4a[4, 5] were found to display aberrant promoter methylation in PC, leading to abnormalities in gene transcription. At the genomic level, by combining the techniques of methylated CGI amplification with Agilent 244 K Human Promoter ChIP-on-chip microarrays, the genome-wide methylation abnormalities in PC have been identified [6, 7].
However, methylation array technology platform-based studies are typically focused on CpG islands (CGIs) and presumably provide less coverage of the entire genome than studies using next-generation sequencing technology . Therefore, the details of the genome-wide methylation profile of PC reported by these studies should be supplemented further, especially in regions such as CGI shores (2 kb regions flanking a CGI), non-CGI promoter regions, and non-gene-associated CGIs (orphan CGIs), as it has already been suggested that methylation changes in these particular regions are associated with certain tumor phenotypes or with tissue specificity [9, 10].
In light of these next-generation sequencing technologies, the methylated portion of the genome identified using methyl capture sequencing (methylCap-seq)  has been profiled in greater detail than using an array-based platform, revealing many novel regions that are differently methylated in a biological sample. Here, we report a comparison of global DNA methylation patterns between pooled PC tissue and pooled adjacent (PN) tissue samples to identify the critical epigenetic effectors responsible for the malignant phenotype of PC. Our study characterized the genome-wide methylation profile of PC and identified the genomic regions displaying a high frequency of aberrant methylation, including regions of gene-associated CGIs, orphan CGIs, CpG shores, and gene promoters lacking CGIs. The aberrant DNA methylation that occurred in these regions was separated into two categories: aberrant DNA methylations that downregulated gene expression, and those that did not affect gene expression. The former category might be related to a tumorigenesis mechanism; thus, these methylations should be studied biofunctionally and might represent targets for tumor treatment, and the latter methylations may be considered potential biomarkers for PC diagnosis.
Wide-spread aberrant hypermethylation in PC and PN revealed via genomic methylation profiling
The clinical characteristics of the patients enrolled in this study are listed in Table 1. The whole-genome methylation profiles of the PC and PN samples were successfully identified using the methylCap-seq method. Exogenous fully methylated and unmethylated spike DNA fragments were used as controls to confirm the capture accuracy of hypermethylated DNA fragments (Additional file 1: Figure S1). We acquired 33,784,358 raw reads in the PC group and 30,868,151 raw reads in the PN group. Based on their alignment with the human genome (hg19) sequence, 16,267,025 (48.15%) raw reads in the PC group and 15,033,135 (48.70%) raw reads in the PN group were uniquely positioned. The reads mapping to 28,691 CGIs, which were defined using the University of California, Santa Cruz (UCSC) Genome Browser, were investigated: 3.57% of the reads in the PC group and 4.25% of the reads in the PN group were positioned at CGIs, resulting in a CGI coverage rate of 64.31% in the PC group and 64.36% in the PN group. These data indicated that our experiment provided considerable information regarding genomic CGIs (Figure 1A).Accumulation of the mapped reads formed peaks. In total, 276,442 and 255,743 peaks were found in the PC and PN samples, respectively, displaying distinct distributions of hypermethylated regions between PC and PN throughout the chromosomes (Figure 1B). An analysis of the hypermethylated peaks approximately 5 kb from the transcription start site (TSS) revealed that methylation peaks accumulated near TSSs, and that more of these peaks were detected in the PC group than in the PN group (Figure 1C). After removing the peaks common to both PC and PN samples (approximately 209,000), 66,807 PC-specific and 46,815 PN-specific hypermethylated differently methylated regions (DMRs) were identified (Figure 1D). Of these hypermethylated DMRs, 36,959 PC-specific and 25,605 PN-specific DMRs were located within genes (Figure 1E), and 1,131 hypermethylated DMRs in PC tissue and 727 hypermethylated DMRs in PN tissue were associated with CGIs. The hypermethylated DMRs were separated into three categories: TSS, intragenic or intergenic (Figure 1E). Subsequently, the location of the TSS DMRs and intragenic DMRs were further determined using structural annotations of the human genome, such as downstream, enhancer, exon, intron, miRNA, promoter, and 5′ UTR (Figure 1F).
This analysis yielded 1,859 DMRs that mapped to CGIs in both the PC and PN samples, and these DMRs represented the vast majority of the hypermethylated CGIs in the refGene category for PC (88%) and PN (87%); orphan CGIs, which are not associated with any known refGene, accounted only for 12% (133 orphan CGIs of 1,131 affected CGIs; Additional file 2: Table S1) and 13% (96 orphan CGIs of 727 affected CGIs; Additional file 3: Table S2) of the CGI-related DMRs in PC and PN tissue, respectively (Figure 1G). There was no difference in the frequency of orphan CGIs among the DMRs between the PC and PN samples. When the gene-associated DMRs were considered, 5,341 such DMRs were identified; we found a higher frequency of CGI-containing promoters in PC tissue (21%, 609/2,955; Additional file 4: Table S3) than in PN tissue (13%, 312/2,386; Additional file 5: Table S4) (Figure 1G), suggesting the importance of CGIs in aberrant DNA methylation during PC tumorigenesis.
It is generally accepted that methylation of a CGI within a promoter is responsible for gene expression silencing. However, two recent studies discovered that certain types of tumor- and tissue-specific methylation occur in CGI shores that play important biological roles. In this study, aside from examining the methylation status of CGIs, we investigated the methylation status of CGI shores, particularly those associated with gene promoters. Regarding the hypermethylated genes and their sites of hypermethylation (in the CGI only, the CGI shore only, or both) in PC tissue, 527 genes (Additional file 6: Table S5) were hypermethylated in a CGI and a CGI shore (502 DMRs); 111 genes (Additional file 7: Table S6) were hypermethylated in a CGI alone (108 DMRs), and 1,278 genes (Additional file 8: Table S7) were hypermethylated in a CGI shore alone (1.242 DMRs). In contrast, based on analysis of the hypermethylated genes and sites of hypermethylation in PN tissue, 333 genes (Additional file 9: Table S8) were hypermethylated in both a CGI and a CGI shore (305 DMRs), 51 genes (Additional file 10: Table S9) were hypermethylated in a CGI alone (47 DMRs), and 1,298 genes (Additional file 11: Table S10) were hypermethylated in a CGI shore only (1,183 DMRs) (Figure 1H). In general, the number of genes aberrantly methylated in a CGI was much fewer, accounting for 20% to 25% of the total genes affected. More frequently, the genes were methylated in a CGI shore, which is probably due to the extended genomic regions defined by the CGI shore. We detected more abnormally hypermethylated genes in a CGI in PC tissue than in PN tissue (P = 0.0002, Chi-square test), suggesting that CGIs are more likely to contain an aberrant DNA methylation target during PC development.
Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of the aberrantly methylated genes in PC
It is well accepted that methylation abnormalities within promoters can influence the expression of the corresponding genes. Therefore, we conducted GO analysis of the genes that displayed promoter hypermethylation in PC and PN tissue. A significance of P < 0.05 indicated gene enrichment in several GO categories (Table 2). We determined that 668 hypermethylated genes in PC tissue (Additional file 12: Table S11) were enriched in ‘sequence-specific DNA binding’ (GO:0043565), ‘neuron differentiation’ (GO:0030182), ‘regulation of transcription, DNA-dependent’ (GO:0006355), or ‘cell morphogenesis involved in differentiation’ (GO:0000904) and that 383 hypermethylated genes in PN tissue (Additional file 13: Table S12) were enriched in ‘plasma membrane part’ (GO:0044459), ‘channel regulator activity’ (GO:0016247), ‘positive regulation of bone morphogenetic protein (BMP) signaling pathway’ (GO:0030513), ‘protein homo oligomerization’ (GO:0051260), or ‘neuron differentiation’ (GO:0030182). Furthermore, we identified 111 genes containing hypermethylated promoters in PC that were enriched in ‘regulation of transcription term’; among these genes, the methylation status of DLX4, ELAVL2, IRX1, PITX2, SIM2, TBX5, and TFAP2C was subsequently validated in the tissue samples by methylation-specific PCR (MSP) (Figure 2B).
A KEGG pathway analysis of the aforementioned genes revealed that hypermethylated genes in both PC and PN tissue were enriched in ‘neuroactive ligand-receptor interaction’ (hsa04080) (Table 2). Furthermore, several miRNAs were aberrantly methylated in PC tissue. The PC-related hypermethylated miRNAs included mir-9-3, mir-9-1, mir-124-3, mir-10b, mir-124-2, mir-718, and mir-203; the PN-related hypermethylated miRNAs included mir-210, mir-1469, mir-130b, mir-149, mir-1224, and mir-564 (Table 3).
We analyzed the hypermethylated genes that were enriched in PC and PN tissue, focusing on their functional involvement in tumorigenesis, and found that the aberrantly methylated genes could be categorized as either tumor promoters or tumor suppressors. The corresponding genes included 20 tumor promoters and 10 tumor suppressors among the PC-related hypermethylated genes and 10 tumor promoters and 5 tumor suppressors among the PN-related hypermethylated genes (Additional file 14: Table S13). These results suggest that aberrant DNA methylation plays an important role in tumor development via the important biological pathways related to the regulation of tumorigenesis.
Verification of PC-specific DMRs identified in methylCap-seq
The accuracy and precision of the DMR profiles were validated in two sample sets, which consisted of the pooled samples used to generate the methylCap-seq library and another independent sample set. In one validation vignette, the DNA methylation status of the ten most significant DMRs that were located in a promoter region (P <10-15) were evaluated using bisulfite sequencing PCR (BSP) in the same set of PC and PN samples used to generate the methylCap-seq libraries. Six candidate DMR genes, C5orf38, EMX1, NPR3, VSTM2B, ELAVL2, and TFAP2C, were validated to be significantly hypermethylated in PC tissue compared with PN tissue using the BSP technique (representative results are shown in Figure 2A).
Preliminary detection of the DMRs identified by genome methylation profiling in limited clinical PC samples
In another validation vignette, the methylation status of 20 gene-associated DMRs scattered throughout various genetic elements, such as promoters, miRNAs, introns, exons, and CGI shores, were analyzed using MSP in paired samples of PC and PN tissue. The results revealed that seven gene-associated DMRs displayed clear differences in methylation between the PC and PN samples (representative results are shown in Figure 2B). The DMRs in the promoter regions of TRADD, AGAP2, and FAM115A displayed a loss of methylation in PC tissue (The MSP results are presented in Additional file 15: Table S14).
MSRE-qPCR and RT-qPCR validation of the methylation of orphan CGIs and the expression of corresponding ESTs in PC cell lines treated with 5-aza-2′-deoxycytidine (5-aza-dc)
Three PC cell lines were treated with 5-aza-2′-deoxycytidine (5-aza-dc). The methylation status of ten hypermethylated DMRs (for gene locus information, see Figure 3A) in promoter CGIs and orphan CGIs was quantitatively analyzed via methylation sensitive restriction enzyme-based qPCR (MSRE-qPCR) (Figure 3B) in these three PC cell lines before and after treatment with 5-aza-dc. The expression levels of gene-associated CGI-containing genes and orphan CGI-containing expressed sequence tags (ESTs) were analyzed using RT-qPCR to ascertain the correlation between the aberrant DMRs and the corresponding mRNA expression levels (Figure 3C). The results indicated that the methylation levels of four orphan CGIs and one promoter CGI were decreased and that the mRNA expression of the corresponding genes or ESTs increased upon 5-aza-dc treatment, suggesting that the expression of these genes or ESTs might be regulated by DNA methylation. Quantitative analysis of the methylation status of these particular DMRs in an independent set of samples (testing group: eight PC samples, five PN samples, and three PC cell lines) confirmed the differences in methylation at these four DMRs in the clinical samples (Figure 3D).
Pooling strategies have been utilized in many previous genomic studies to investigate the phenotypic similarities between specific models, owing to the advantage of conserving samples [12, 13]. In this study, genome-wide methylation profiles of PC and PN tissues were established using corresponding pooled samples.
We identified 5,280 and 3,488 hypermethylated DMRs in PC and PN tissue, respectively, that were closely associated with genes and CGIs. Gene ontology analysis of the genes associated with these DMRs revealed that the aberrantly hypermethylated genes primarily belonged to categories related to nucleic acid binding, DNA binding, and activation of transcription factors, suggesting that the methylation of these affected genes coupled with the downregulation of RNA expression resulted in the decreased expression of other genes. Studies of non-small-cell lung carcinoma by Helman et al. and Zhao et al.  demonstrated that methylation-enriched genes displayed aberrant methylation and RNA expression in multiple tumor types; these genes were referred to as tumor suppressor genes. The methylation-enriched genes in PC associated with ‘cell morphogenesis involved in differentiation’ (GO: 0000904) may participate in the morphological changes and late-stage differentiation in PC tissue. In contrast, the hypermethylated genes in PN tissue were predominantly enriched in ‘plasma membrane part’ (GO:0044459) and ‘channel regulator activity’ (GO:0016247), and it has been confirmed that these genes are enriched in other tumor types in genomic methylation studies, suggesting that abnormal formation of the plasma membrane might be a common characteristic of tumor differentiation and maintenance [16, 17].
We also conducted KEGG pathway analysis of the genes containing DMRs in their promoter; surprisingly, the hypermethylated genes in both the PC and PN samples were enriched in the same pathway, ‘neuroactive ligand-receptor interaction’, which primarily participates in the endocrine and exocrine functions of cells. Functional abnormalities in these genes have been demonstrated in studies of meningioma and PC [18, 19]. Our study revealed the identical direction of methylation changes in this pathway in PC and PN tissue, but the detailed mechanism by which this pathway participates in PC development requires further investigation (Table 2).
As epigenetic factors, miRNAs play an important role in the regulation of cellular biophysical functions, and have been shown to be associated with the apoptosis, invasion, metastasis, recurrence, and drug resistance of tumor cells. The regulation of miRNAs by DNA methylation has been extensively studied. In this study, we identified particular miRNAs that were aberrantly methylated in PC tissue. Among these miRNAs, hsa-mir-124-3 has been shown to be hypermethylated in PC and is involved in the progression, metastasis, and recurrence of PC [20, 21]. Alternatively, mir-130b and mir-210 are two hypermethylated miRNAs in PN tissue that were discovered in this study, both of which were found to be more strongly expressed in PC and have been associated with the proliferation and invasion of PC cells [22, 23]. Certain miRNAs that we identified as aberrantly methylated in PC tissue have already been studied in other tumor types and are regulated by DNA methylation. For example, the hypermethylated miRNAs identified in PC tissue, such as miR-9-3, mir-9-1, miR-124, and miR-203, are also hypermethylated in non-small-cell lung carcinoma, breast cancer, cervical cancer, and hematological cancer, respectively; hypermethylation of these sites decreases miRNA expression, which promotes tumor development and tumor cell proliferation [24–27]. In addition, mir-10b, which was found to be hypermethylated in PC tissue in this study, was thought to be an inhibitor of tumor metastasis in animal models  and was found to be more strongly expressed in PC . Therefore, further investigation of the remaining miRNAs that were aberrantly methylated in PC will shed light on the mechanisms underlying pancreatic carcinogenesis, as the related research is currently very limited.
The top 40 genes, based on their P value, that contained methylation changes in their promoter region were selected and examined using MSP in the test group samples. Of these genes, 18 displayed significant differences in methylation between PC and PN tissues (or PC cell lines). Interestingly, among these genes, seven (DLX4, ELAVL2, IRX1, PITX2, SIM2, TBX5, and TFAP2C) were enriched for the annotation of regulation of transcription (GO: 0006355), which corresponds to the results obtained in our previous GO analysis. All of the genes discussed above have been investigated in PC and other tumor types, and their involvement in carcinogenesis has been confirmed.
At both the genetic and genomic levels, many hypermethylated genes previously reported in PC studies were identified as hypermethylated DMRs in the present study, including LHX1, FOXE1, PAX6, BNIP3, ALPP, CEBPA, CACNA1G, CCND2, BAI1, NRN1, PENK, FAM84A, and ZNF415. In addition, our study identified other genes that are frequently hypermethylated in different types of cancer, such as RASSF1a, CDKN2A, hHML1, and CDH1[34, 35]. Thus, we have established a relatively extensive database of abnormally methylated sites in PC. We also compared our data with those reported by Omura et al. , who analyzed nine pairs of PC and PN samples using human CGI microarray 244 k chips, obtaining (after data filtering using the appropriate thresholds) 1,658 differently methylated known loci. This comparison revealed the following. (1) Regarding the ability to capture aberrantly methylated gene targets, methylCap-seq identified more hypermethylated genes than the array method in PC tissue (1983 versus 1206) and PN tissue (1692 versus 379), indicating that methylated DNA fragment enrichment followed by deep sequencing identifies additional aberrant gene loci, although it is more labor-intensive and time-consuming. (2) A total of 737 genes (Additional file 16: Table S15) was identified by both the methylCap-seq and array methods, accounting for 46.7% of the total genes recovered. This high recovery rate between the two methods reflects the reliability of these methods for this purpose. However, the unique genes that were identified suggest that these two methods each have their own particular advantages. (3) The high percentage of commonality among the PC-specific hypermethylated DMR genes (46.2%) and low percentage of commonality among the PN-specific hypermethylated DMR genes (30 genes, 2.9% of all the PN-related hypermethylated DMR genes) (Additional file 17: Table S16) between the study of Omura et al.  and this study suggests that during the entire process of PC development, hypermethylation is a relatively defined and destined process, whereas hypermethylation at the initiation of tumorigenesis is relatively random or perhaps stochastic.
We also compared our aberrant methylated gene targets determined in this study with the results of previous studies. Among all 3,911 differently methylated genes (DMGs) identified in this study, 728 DMGs were reported by Nourse et al., 339 DMGs were reported by Vincent et al., and 55 DMGs were reported by Tan et al.  (Additional file 18: Table S17). This discrepancy in the number of DMGs obtained between the four groups might result from the different technological platforms adopted by each study group, as well as the different ethnic backgrounds of the enrolled patients. Furthermore, this discrepancy emphasizes that the array-based and sequencing-based DNA methylation assay methods must be applied alternately to complement one another, to elucidate DNA methylation at the genomic level.
Deaton and Bird separated CGIs into three categories, TSS, intragenic, and intergenic, with the latter two categories defined as orphan CGIs . Despite poor understanding of the functions of orphan CGIs until recently, studies have shown that orphan CGIs are involved in the regulation of gene transcription, genomic imprinting, and non-coding RNA transcription and that orphan CGIs might display tissue-specific methylation profiles . In this study, we investigated the methylation status of particular orphan CGIs in PC. Hypermethylated orphan CGIs have been found in PC. The methylation status of orphan CGIs was closely associated with the transcription levels of nearby non-annotated ESTs. Further studies should be conducted to clarify whether the methylation-regulated ESTs containing orphan CGIs are potential genes or gene elements. It is well known that merely 6.8% of CpGs are located in CGIs. The methylation status and biological functions of the other 93.2% of CpGs have yet to be adequately studied. The study by Yu et al. indicated that the methylation of CpGs in CGI shores is involved in regulating gene transcription or establishing tissue-specific methylation patterns, and changes in the methylation status of CpGs in CGI shores might occur at an earlier stage in carcinogenesis than the changes that occur in gene-associated CGIs. Our study suggests that the methylation changes in all of these CpG regions in PC are indispensable components of the genomic methylation profile of PC and may influence the transcription of PC-related genes and non-coding RNAs, potentially affecting tissue-specific cell differentiation and ultimately leading to carcinogenesis.
It is generally accepted that abnormal DNA hypermethylation can either downregulate gene expression (gene silencing) or exert no influence on gene expression. Gene silencing-related DNA methylation might be involved in PC development, and, therefore, these methylation sites may be examined in tumor development studies and considered as treatment targets. Alternatively, the DNA methylations that are not associated with gene expression might serve as biomarkers of the specific state of PC, and, therefore, these methylation sites could be used for clinical diagnosis. We anticipate that our comprehensive analysis of PC gene methylation will facilitate the further investigation of PC biomarkers for diagnostic, prognostic, and therapeutic applications, by: (1) improving the understanding of the epigenetic importance of DNA methylation in PC tumorigenesis, which might be located outside or within CGIs, including both orphan and gene-associated CGIs; and (2) providing additional candidate targets for PC diagnosis, prognosis and treatment (not limited to the previously reported CGIs and promoters). Moreover, the enormous number of targets (many thousands) obtained suggests a vast heterogeneity among PC patients, thus requiring a large patient cohort for chip-based analyses to validate, optimize, and establish potential targets for ultimate clinical application.
In the present study, the genome-wide methylation profiles of PC and PN tissues were established using methylCap-seq, revealing globally reprogrammed and deregulated DNA methylation in PC. Compared with PN tissue, there were many PC-specific aberrations in the hypermethylation of CpGs in TSS CGIs, orphan CGIs, CGI shores, and promoter regions lacking CGIs. These findings will be helpful in elucidating the mechanisms underlying pancreatic carcinogenesis related to the DNA methylation-regulated expression of genes and non-coding RNAs. Furthermore, the aberrantly methylated genes in PC identified in this study might serve as potential biomarkers for the early diagnosis and treatment of this deadly disease.
Pancreatic cancer tissue samples were collected from 18 patients who had undergone surgical treatment without receiving preoperative chemotherapy or radiotherapy from May 2009 to March 2011 in Renji Hospital, School of Medicine, Shanghai Jiao Tong University. A diagnosis of PC was confirmed by histological examination. Resected tumor tissues and matched normal tissues at least 2 cm away from the tumor tissues were collected during the operation, labeled, and stored at -80°C. The sixth edition of the Tumor, Node, and Metastasis (TNM) Staging System proposed by the International Union against Cancer  was utilized to stage the tumor tissue samples (Table 1). The study was approved by the medical ethics committee at Renji Hospital, School of Medicine, Shanghai Jiao Tong University. All the patients signed an informed consent form. DNA was isolated from frozen tissues or cell lines using a conventional proteinase K and organic extraction method, as previously described .
Genome-wide methylation profiling by methylCap-seq
Genomic DNA was extracted from ten PC tissues and ten matched normal tissues. Equal amounts of DNA were mixed to form the PC and PN groups. Pooled DNA (1.2 μg) from each group was used to generate the library for methylCap-seq as previously described .
Mapping the sequence reads and DMR identification and annotation
We used Burrows-Wheeler alignment tools  with the default settings to map the 36 bp unpaired reads to the hg19 human genome reference assembly . After removing PCR duplicates using Picard, the aligned data were converted, sorted, and indexed using Samtools  and Picard .
Methylation peaks (hypermethylated regions) were identified using model-based analysis of ChIP-Seq in the PC and PN samples, as previously described . The DMRs between PC and PN were identified using two methods, model-based analysis of ChIP-Seq  and a bi-asymmetric-Laplace model (BALM) , to increase the detective power of methylCap-seq. To decrease the false positive detection of DMR using BALM, the dual-threshold strategy was applied. A high-confidence threshold (0.975) was utilized in the PC hypermethylated region screening, and a low-confidence threshold (0.950) was utilized in the PN hypermethylated region screening. Cancer-specific methylation peaks were defined as hypermethylated regions. Similarly, normal tissue-specific methylation peaks were defined as hypomethylated regions with the reverse settings. Whole-genome methylation (methylation of each CpG) was inferred using BALM, which was processed for a Pearson correlation analysis among all the samples in the R environment. The refSeq genes (UCSC genes) and corresponding CGIs were downloaded from the table browser of the UCSC database . The browser extensible data (BED) file operations were performed using BEDTools  and other Perl scripts. All the scripts are available upon email request. The generated genomic methylation profile was uploaded to a public database (Gene Expression Omnibus: GSE54854). Gene ontology analysis was performed using DAVID Bioinformatics Resources 6.7 .
In this study, BSP was utilized to determine the methylation status at single CpG resolution of DMRs identified by genomic methylation profiling; MSP was performed for qualitative methylation screening in a small set of PC samples. Using MSRE-qPCR, the DNA methylation status in orphan CGIs was quantitatively analyzed in PC cell lines before and after 5-aza-2′-deoxycytidine treatment and in small samples of clinical PC tissues, as described previously . Approximately 1.0 μg of genomic DNA extracted from PC or PN samples or PC cell lines was bisulfate-treated using EpiTect Kit (Qiagen, Hilden, Germany). Primers for MSP and BSP were designed using MethPrimer, an online primer design tool . The MSRE-qPCR primers were designed using Primer3 . The sequences of the primers utilized in this study are listed in Additional files 19 and 20: Tables S18 and S19. Jumpstart Taq (Sigma-Aldrich, St. Louis, MO, USA) was used in BSP and MSP with a 20 μl reaction volume per tube. The BSP and MSP reaction conditions were as follows: 94°C for 3 min; 35 cycles of 94°C for 20 s, annealing for 20 s, and 72°C for 20 s; and 72°C for 5 min. The PCR products were analyzed by electrophoresis in 1.5% agarose gels. The PCR products were TA cloned and verified by sequencing. At least five clones were sequenced for each BSP reaction.
Cell culture and 5-aza-2′-deoxycytidine treatment
Three pancreatic adenocarcinoma cell lines were used: BxPC-3 ATCC, CRL-1687), PANC-1 (ATCC, CRL-1469), and CFPAC-1 (ATCC, CRL-1918). All the cell lines were cultured in RPMI1640 supplemented with 10% FBS, 100 U/ml penicillin, and 100 U/ml streptomycin. All the cell lines were maintained at 37°C in a humidified atmosphere with 5% CO2.
The restoration of gene expression by demethylation was evaluated in the BxPC-3, CFPAC-1, and CFPAC-1 cell lines. For the CpG demethylation analysis, exponentially growing cells were seeded at a density of 1.5 × 106 cells/100 mm dish and allowed to attach overnight. The cells were then treated with freshly prepared 5-aza-dC (5.0 μM; Sigma-Aldrich, St. Louis, MO, USA) for 3 days.
RNA isolation and real-time PCR
Total RNA was prepared from cultured cells using Trizol reagent according to the manufacturer’s instructions (Invitrogen, USA) and then reverse transcribed using an oligo (dT) primer and SuperscriptII RNase H-Reverse Transcriptase (Invitrogen, USA). Real-time PCR was performed with primer pairs for the EST expression assay, and GAPDH was used as the internal control. Real-time PCR was performed as follows: 94°C for 3 min followed by 40 cycles of 94°C for 10 s, 62°C for 10 s, and 72°C for 15 s. Real-time qPCR was performed to detect GAPDH expression with an SYBR Green PCR Kit (Applied Biosystems, Foster city, CA, USA) on a ROTOR-GENE 6000 Real-Time PCR System (ROTOR-GENE).
Statistical calculations were performed using the SPSS statistical software package (Version 13.0; SPSS, Inc. Chicago, IL). The measurement data were analyzed using one-way ANOVA. Statistical significance was considered for P < 0.05.
analysis of variance
browser extensible data
bone morphogenetic protein
bisulfite sequencing PCR
differently methylated gene
differently methylated region
expressed sequence tag
fetal bovine syndrome
Kyoto Encyclopedia of Genes and Genomes
methylation sensitive restriction enzyme-based qPCR
polymerase chain reaction
non-tumor tissue adjacent to pancreatic cancer
quantitative polymerase chain reaction
reverse transcription quantitative polymerase chain reaction
tumor, node, and metastasis
transcription start site, UCSC, University of California, Santa Cruz
Jemal A, Siegel R, Ward E, Murray T, Xu J, Thun MJ: Cancer statistics, 2007. CA Cancer J Clin. 2007, 57 (1): 43-66. 10.3322/canjclin.57.1.43.
Ma C, Jiang YX, Liu SZ, Quan PL, Sun XB, Zheng RS, Zhang SW, Chen WQ: Trend and prediction on the incidence of pancreatic cancer in China. Zhonghua Liu Xing Bing Xue Za Zhi. 2013, 34 (2): 160-163.
Kanwal R, Gupta S: Epigenetic modifications in cancer. Clin Genet. 2012, 81: 303-311. 10.1111/j.1399-0004.2011.01809.x.
Klump B, Hsieh CJ, Nehls O, Dette S, Holzmann K, Kiesslich R, Jung M, Sinn U, Ortner M, Porschen R, Gregor M: Methylation status of p14ARF and p16INK4a as detected in pancreatic secretions. Br J Cancer. 2003, 88: 217-222. 10.1038/sj.bjc.6600734.
Dammann R, Schagdarsurengin U, Liu L, Otto N, Gimm O, Dralle H, Boehm BO, Pfeifer GP, Hoang-Vu C: Frequent RASSF1A promoter hypermethylation and K-ras mutations in pancreatic carcinoma. Oncogene. 2003, 22: 3806-3812. 10.1038/sj.onc.1206582.
Omura N, Li CP, Li A, Hong SM, Walter K, Jimeno A, Hidalgo M, Goggins M: Genome-wide profiling of methylated promoters in pancreatic adenocarcinoma. Cancer Biol Ther. 2008, 7: 1146-1156. 10.4161/cbt.7.7.6208.
Vincent A, Omura N, Hong SM, Jaffe A, Eshleman J, Goggins M: Genome-wide analysis of promoter methylation associated with gene expression profile in pancreatic adenocarcinoma. Clin Cancer Res. 2011, 17 (13): 4341-4354. 10.1158/1078-0432.CCR-10-3431.
Bock C, Tomazou EM, Brinkman AB, Müller F, Simmer F, Gu H, Jäger N, Gnirke A, Stunnenberg HG, Meissner A: Quantitative comparison of genome-wide DNA methylation mapping technologies. Nat Biotechnol. 2010, 28: 1106-1114. 10.1038/nbt.1681.
Irizarry RA, Ladd-Acosta C, Wen B, Wu Z, Montano C, Onyango P, Cui H, Gabo K, Rongione M, Webster M, Ji H, Potash JB, Sabunciyan S, Feinberg AP: The human colon cancer methylome shows similar hypo- and hypermethylation at conserved tissue-specific CpG island shores. Nat Genet. 2009, 41: 178-186. 10.1038/ng.298.
Doi A, Park IH, Wen B, Murakami P, Aryee MJ, Irizarry R, Herb B, Ladd-Acosta C, Rho J, Loewer S, Miller J, Schlaeger T, Daley GQ, Feinberg AP: Differential methylation of tissue and cancer-specific CpG island shores distinguishes human induced pluripotent stem cells, embryonic stem cells and fibroblasts. Nat Genet. 2009, 41: 1350-1353. 10.1038/ng.471.
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-399. 10.1093/nar/gkp992.
Li LM, Hu ZB, Zhou ZX, Chen X, Liu FY, Zhang JF, Shen HB, Zhang CY, Zen K: Serum microRNA profiles serve as novel biomarkers for HBV infection and diagnosis of HBV-positive hepatocarcinoma. Cancer Res. 2010, 70 (23): 9798-9807. 10.1158/0008-5472.CAN-10-1001.
Deaton AM, Webb S, Kerr AR, Illingworth RS, Guy J, Andrews R, Bird A: Cell type-specific DNA methylation at intragenic CpG islands in the immune system. Genome Res. 2011, 21 (7): 1074-1086. 10.1101/gr.118703.110.
Helman E, Naxerova K, Kohane IS: DNA hypermethylation in lung cancer is targeted at differentiation-associated genes. Oncogene. 2012, 31 (9): 1181-1188. 10.1038/onc.2011.307.
Zhao M, Sun J, Zhao Z: TSGene: a web resource for tumor suppressor genes. Nucleic Acids Res. 2013, 41 (Database issue): D970-D976.
Cicek MS, Koestler DC, Fridley BL, Kalli KR, Armasu SM, Larson MC, Wang C, Winham SJ, Vierkant RA, Rider DN, Block MS, Klotzle B, Konecny G, Winterhoff BJ, Hamidi H, Shridhar V, Fan JB, Visscher DW, Olson JE, Hartmann LC, Bibikova M, Chien J, Cunningham JM, Goode EL: Epigenome-wide ovarian cancer analysis identifies a methylation profile differentiating clear-cell histology with epigenetic silencing of the HERG K+ channel. Hum Mol Genet. 2013, 22 (15): 3038-3047. 10.1093/hmg/ddt160.
Sriraksa R, Zeller C, Dai W, Siddiq A, Walley AJ, Limpaiboon T, Brown R: Aberrant DNA methylation at genes associated with a stem cell-like phenotype in cholangiocarcinoma tumors. Cancer Prev Res (Phila). 2013, 6 (12): 1348-1355. 10.1158/1940-6207.CAPR-13-0104.
Wei P, Tang H, Li D: Insights into pancreatic cancer etiology from pathway analysis of genome-wide association study data. PLoS One. 2012, 7 (10): e46887-10.1371/journal.pone.0046887.
Keller A, Ludwig N, Backes C, Romeike BF, Comtesse N, Henn W, Steudel WI, Mawrin C, Lenhof HP, Meese E: Genome wide expression profiling identifies specific deregulated pathways in meningioma. Int J Cancer. 2009, 124 (2): 346-351. 10.1002/ijc.23942.
Wang P, Chen L, Zhang J, Chen H, Fan J, Wang K, Luo J, Chen Z, Meng Z, Liu L: Methylation-mediated silencing of the miR-124 genes facilitates pancreatic cancer progression and metastasis by targeting Rac1. Oncogene. 2014, 33 (4): 514-10.1038/onc.2012.598.
Gebauer K, Peters I, Dubrowinskaja N, Hennenlotter J, Abbas M, Scherer R, Tezval H, Merseburger AS, Stenzl A, Kuczyk MA, Serth J: Hsa-mir-124-3 CpG island methylation is associated with advanced tumours and disease recurrence of patients with clear cell renal cell carcinoma. Br J Cancer. 2013, 108 (1): 131-138. 10.1038/bjc.2012.537.
Zhao G, Zhang JG, Shi Y, Qin Q, Liu Y, Wang B, Tian K, Deng SC, Li X, Zhu S, Gong Q, Niu Y, Wang CY: MiR-130b is a prognostic marker and inhibits cell proliferation and invasion in pancreatic cancer through targeting STAT3. PLoS One. 2013, 8: e73803-10.1371/journal.pone.0073803.
Takikawa T, Masamune A, Hamada S, Nakano E, Yoshida N, Shimosegawa T: miR-210 regulates the interaction between pancreatic cancer cells and stellate cells. Biochem Biophys Res Commun. 2013, 437 (3): 433-439. 10.1016/j.bbrc.2013.06.097.
Heller G, Weinzierl M, Noll C, Babinsky V, Ziegler B, Altenberger C, Minichsdorfer C, Lang G, Döme B, End-Pfützenreuter A, Arns BM, Grin Y, Klepetko W, Zielinski CC, Zöchbauer-Müller S: Genome-wide miRNA expression profiling identifies miR-9-3 and miR-193a as targets for DNA methylation in non-small cell lung cancers. Clin Cancer Res. 2012, 18 (6): 1619-1629. 10.1158/1078-0432.CCR-11-2450.
Lehmann U, Hasemeier B, Christgen M, Müller M, Römermann D, Länger F, Kreipe H: Epigenetic inactivation of microRNA gene hsa-mir-9-1 in human breast cancer. J Pathol. 2008, 214 (1): 17-24. 10.1002/path.2251.
Wilting SM, van Boerdonk RA, Henken FE, Meijer CJ, Diosdado B, Meijer GA, le Sage C, Agami R, Snijders PJ, Steenbergen RD: Methylation-mediated silencing and tumour suppressive function of hsa-miR-124 in cervical cancer. Mol Cancer. 2010, 9: 167-
Chim CS, Wong KY, Leung CY, Chung LP, Hui PK, Chan SY, Yu L: Epigenetic inactivation of the hsa-miR-203 in haematological malignancies. J Cell Mol Med. 2011, 15 (12): 2760-2767. 10.1111/j.1582-4934.2011.01274.x.
Ma L, Reinhardt F, Pan E, Soutschek J, Bhat B, Marcusson EG, Teruya-Feldstein J, Bell GW, Weinberg RA: Therapeutic silencing of miR-10b inhibits metastasis in a mouse mammary tumor model. Nat Biotechnol. 2010, 28: 341-347. 10.1038/nbt.1618.
Nakata K, Ohuchida K, Mizumoto K, Kayashima T, Ikenaga N, Sakai H, Lin C, Fujita H, Otsuka T, Aishima S, Nagai E, Oda Y, Tanaka M: MicroRNA-10b is overexpressed in pancreatic cancer, promotes its invasiveness, and correlates with a poor prognosis. Surgery. 2011, 150 (5): 916-922. 10.1016/j.surg.2011.06.017.
Sato N, Fukushima N, Maitra A, Matsubayashi H, Yeo CJ, Cameron JL, Hruban RH, Goggins M: Discovery of novel targets for aberrant methylation in pancreatic carcinoma using high-throughput microarrays. Cancer Res. 2003, 63: 3735-3742.
Kumagai T, Akagi T, Desmond JC, Kawamata N, Gery S, Imai Y, Song JH, Gui D, Said J, Koeffler HP: Epigenetic regulation and molecular characterization of C/EBPα in pancreatic cancer cells. Int J Cancer. 2009, 124: 827-833. 10.1002/ijc.23994.
Ueki T, Toyota M, Sohn T, Yeo CJ, Issa JP, Hruban RH, Goggins M: Hypermethylation of multiple genes in pancreatic adenocarcinoma. Cancer Res. 2000, 60: 1835-1839.
Matsubayashi H, Sato N, Fukushima N, Yeo CJ, Walter KM, Brune K, Sahin F, Hruban RH, Goggins M: Methylation of cyclin D2 is observed frequently in pancreatic cancer but is also an age-related phenomenon in gastrointestinal tissues. Clin Cancer Res. 2003, 9: 1446-1452.
Omura N, Goggins M: Epigenetics and epigenetic alterations in pancreatic cancer. Int J Clin Exp Pathol. 2009, 2: 310-326.
Schutte M, Hruban RH, Geradts J, Maynard R, Hilgers W, Rabindran SK, Moskaluk CA, Hahn SA, Schwarte-Waldhoff I, Schmiegel W, Baylin SB, Kern SE, Herman JG: Abrogation of the Rb/p16 tumor-suppressive path way in virtually all pancreatic carcinomas. Cancer Res. 1997, 57: 3126-3130.
Nourse C, Kazakoff S, Taylor D, Leonard C, Chang DK, Jones MD, Thomas M, Watson C, Pinese M, Cowley M, Rooman I, Pajic M, Butturini G, Malpaga A, Corbo V, Crippa S, Falconi M, Zamboni G, Castelli P, Lawlor RT, Gill AJ, Scarpa A, Pearson JV, Biankin AV, Grimmond SM, APGI: Genome-wide DNA methylation patterns in pancreatic ductal adenocarcinoma reveal epigenetic deregulation of SLIT-ROBO, ITGA2 and MET signaling. Int J Cancer. 2014, 135: 1110-1118. 10.1002/ijc.28765.
Tan AC, Jimeno A, Lin SH, Wheelhouse J, Chan F, Solomon A, Rajeshkumar NV, Rubio-Viqueira B, Hidalgo M: Characterizing DNA methylation patterns in pancreatic cancer genome. Mol Oncol. 2009, 3: 425-438. 10.1016/j.molonc.2009.03.004.
Deaton AM, Bird A: CpG islands and the regulation of transcription. Genes & Dev. 2011, 25: 1010-1022. 10.1101/gad.2037511.
Illingworth RS, Gruenewald-Schneider U, Webb S, Kerr AR, James KD, Turner DJ, Smith C, Harrison DJ, Andrews R, Bird AP: Orphan CpG islands identify numerous conserved promoters in mammalian genome. PLoS Genet. 2010, 6: e1001134-10.1371/journal.pgen.1001134.
Yu W, Jin C, Lou X, Han X, Li L, He Y, Zhang H, Ma K, Zhu J, Cheng L, Lin B: Global analysis of DNA methylation by methyl-capture sequencing reveals epigenetic control of cisplatin resistance in ovarian cancer cell. PLoS One. 2011, 6 (12): e29450-10.1371/journal.pone.0029450.
Greene FL, Page DL, Fleming ID, Fritz AG, Balch CM, Haller DG, Morrow M: AJCC cancer staging handbook: TNM classification of malignant tumors. 2002, New York: Springer-Verlag, 171-177. 6
Zhao Y, Guo S, Sun J, Huang Z, Zhu T, Zhang H, Gu J, He Y, Wang W, Ma K, Wang J, Yu J: Methylcap-Seq reveals novel DNA methylation markers for the diagnosis and recurrence prediction of bladder cancer in a Chinese population. PLoS One. 2012, 7: e35175-10.1371/journal.pone.0035175.
Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25 (14): 1754-1760. 10.1093/bioinformatics/btp324.
Kuhn RM, Haussler D, Kent WJ: The UCSC genome browser and associated tools. Brief Bioinform. 2012, 14 (2): 144-161.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R: The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009, 25 (16): 2078-2079. 10.1093/bioinformatics/btp352.
Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, Nusbaum C, Myers RM, Brown M, Li W, Liu XS: Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008, 9: R137-10.1186/gb-2008-9-9-r137.
Lan X, Adams C, Landers M, Dudas M, Krissinger D, Marnellos G, Bonneville R, Xu M, Wang J, Huang TH, Meredith G, Jin VX: High resolution detection and analysis of CpG dinucleotides methylation using MBD-Seq technology. PLoS One. 2011, 6 (7): e22226-10.1371/journal.pone.0022226.
Quinlan AR, Hall IM: BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010, 26 (6): 841-842. 10.1093/bioinformatics/btq033.
David Bioinformatics Resources 6.7. [http://david.abcc.ncifcrf.gov/]
Huang ZH, Hu Y, Hua D, Wu YY, Song MX, Cheng ZH: Quantitative analysis of multiple methylated genes in plasma for the diagnosis and prognosis of hepatocellular carcinoma. Exp Mol Pathol. 2011, 91 (3): 702-707. 10.1016/j.yexmp.2011.08.004.
We are grateful to all the staff of Biliary and Pancreatic Surgery Department in Shanghai Renji Hospital for providing valuable samples and clinical information. The authors would like to acknowledge the following funding source: National Natural Science Foundation of China (81372368, JY); Shanghai Municipal Commission of Health and Family Planning Foundation (20134055, YXZ); the State Key Laboratory of Oncogenes and Related Genes Foundation (No.91-11-01, NO.91-13-10, JY).
The authors declare that they have no competing interests.
The experiments were conceived and designed by JY. The experiments were performed by YZ, JS, JG, WW, and XZ. Data were analyzed by SG, HZ, and NT. The paper was written by JY, YZ, and JS. All authors read and approved the final manuscript.
Yangxing Zhao, Jinfeng Sun, Hongyu Zhang, Shicheng Guo contributed equally to this work.