Abstract

To establish diagnostic criteria for ductal adenocarcinomas of the pancreas (PCs), bacterial artificial chromosome (BAC) array-based methylated CpG island amplification was performed using 139 tissue samples. Twelve BAC clones, for which DNA methylation status was able to discriminate cancerous tissue (T) from noncancerous pancreatic tissue in the learning cohort with a specificity of 100%, were identified. Using criteria that combined the 12 BAC clones, T-samples were diagnosed as cancers with 100% sensitivity and specificity in both the learning and validation cohorts. DNA methylation status on 11 of the BAC clones, which was able to discriminate patients showing early relapse from those with no relapse in the learning cohort with 100% specificity, was correlated with the recurrence-free and overall survival rates in the validation cohort and was an independent prognostic factor by multivariate analysis. Genome-wide DNA methylation profiling may provide optimal diagnostic markers and prognostic indicators for patients with PCs.

1. Introduction

It is known that DNA hypomethylation results in chromosomal instability as a result of changes in chromatin structure and that DNA hypermethylation of CpG islands silences tumor-related genes in cooperation with histone modification in human cancers [15]. The incidence of DNA methylation alterations is generally high in cancers of various organs, and particular DNA methylation profiles are significantly associated with poorer tumor differentiation, tumor aggressiveness, and poor prognosis [68]. Moreover, unlike alterations of mRNA and protein expression, which can be easily affected by the microenvironment of cancer cells, DNA methylation alterations are stably preserved on DNA double strands by covalent bonds and can be detected using highly sensitive methodology. Therefore, alterations of DNA methylation can become optimal diagnostic markers of cancers and prognostic indicators for affected patients.

With regard to pancreatic carcinogenesis, we have reported that accumulation of DNA methylation of tumor-related genes [9] is associated with overexpression of DNA methyltransferase (DNMT) 1 [10], the major DNMT, even in peripheral pancreatic duct epithelia with an inflammatory background, in comparison with normal peripheral pancreatic duct epithelia. Ductal adenocarcinomas of the pancreas frequently develop after chronic damage due to pancreatitis, and at least a proportion of peripheral pancreatic duct epithelia with an inflammatory background are at the precancerous stage [11]. The average number of methylated tumor-related genes and the incidence of DNMT1 overexpression increase progressively with the progression of another precancerous lesion, pancreatic intraductal neoplasia [12], to well-differentiated ductal adenocarcinoma, and finally to poorly differentiated ductal adenocarcinoma, suggesting that DNA methylation alterations participate in multistage pancreatic carcinogenesis [9, 10]. However, even though we and other groups have examined the DNA methylation status of several specific tumor-related genes [9, 1317], only a few previous studies have employed recently developed array-based technology for analysis of DNA methylation in ductal adenocarcinomas of the pancreas [18, 19]. To our knowledge, no diagnostic criteria have yet been established for pancreatic cancers on the basis of such genome-wide DNA methylation profiling.

In the present study, in order to obtain diagnostic markers and prognostic indicators of ductal adenocarcinomas of the pancreas, we performed bacterial artificial chromosome (BAC) array-based methylated CpG island amplification (BAMCA) [2022], which is a technique suitable for overviewing the DNA methylation tendency of individual large regions among all chromosomes [23, 24], in samples of normal pancreatic tissue obtained from patients without ductal adenocarcinomas (C), noncancerous pancreatic tissue obtained from patients with ductal adenocarcinomas (N), and cancerous tissue (T).

2. Materials and Methods

2.1. Patients and Tissue Samples

Ninety-one T-samples were obtained from surgically resected specimens from patients with ductal adenocarcinomas who underwent pancreatectomy at the National Cancer Center Hospital, Tokyo, Japan, between 2003 and 2008. From 33 of the 91 patients, N-samples were also obtained from the same surgically resected specimens. Microscopic examination of the histological specimens taken from a region immediately adjoining that from which N-samples had been obtained revealed various degrees of chronic pancreatitis, but no contaminating cancer cells. Fifteen C-samples were obtained from patients without ductal adenocarcinomas who underwent pancreatectomy for metastasis of renal cell carcinoma (1 patient), adenocarcinoma of the gallbladder (3 patients), adenocarcinoma of the papilla of Vater (6 patients), serous cystadenoma (1 patient), mucinous cystadenoma (1 patient), solid-pseudopapillary neoplasm (1 patient), endocrine tumor (1 patient) of the pancreas, and lymphoplasmacytic pancreatitis (1 patient). The total samples were randomly divided into a learning cohort (8 C-, 17 N-, and 46 T-samples) and a validation cohort (7 C-, 16 N-, and 45 T-samples). In the learning cohort, patients from whom C-, N-, and T-samples were obtained comprised 5 men and 3 women with a mean age of 6 9 . 6 ± 8 . 1 (mean ± SD) years, 6 men and 11 women with a mean age of 6 7 . 6 ± 1 0 . 1 years, and 28 men and 18 women with a mean age of 6 4 . 2 ± 1 0 . 8 years, respectively. In the validation cohort, the patients from whom C-, N-, and T-samples were obtained comprised 3 men and 4 women with a mean age of 6 2 . 9 ± 1 8 . 2 years, 11 men and 5 women with a mean age of 6 5 . 0 ± 8 . 7 years, and 27 men and 18 women with a mean age of 6 4 . 6 ± 9 . 7 years, respectively. The clinicopathological parameters of patients who provided T-samples in both the learning and validation cohorts are summarized in Table 1. This study was approved by the Ethics Committee of the National Cancer Center, Tokyo, Japan, and was performed in accordance with the Declaration of Helsinki, 1995. All patients gave their informed consent prior to their inclusion in this study.

2.2. BAMCA

High-molecular-weight DNA from fresh frozen tissue samples was extracted using phenol-chloroform, followed by dialysis. DNA methylation status was analyzed by BAMCA using a custom-made array (MCG Whole Genome Array-4500) harboring 4361 BAC clones located throughout chromosomes 1 to 22, X and Y [25], as described previously [23, 26, 27]. Briefly, a mixture of normal pancreatic tissue DNA obtained from 8 C-samples in the learning cohort was used as a reference for all analyses of test DNA samples in both the learning and validation cohorts. Five-microgram aliquots of test or reference DNA were first digested with 100 units of methylation-sensitive restriction enzyme Sma I (NEB, Ipswich, MA) and subsequently with 20 units of methylation-insensitive Xma I (NEB). Adapters were ligated to Xma I-digested sticky ends, and polymerase chain reaction (PCR) was performed with an adapter primer set. Test and reference PCR products were labeled by random priming with Cy3- and Cy5-dCTP (GE Healthcare, Buckinghamshire, UK), respectively, and precipitated together with ethanol in the presence of Cot-I DNA (Invitrogen, Carlsbad, CA). The mixture was applied to array slides and incubated at 43°C for 63 h. Arrays were scanned with a GenePix Personal 4100A (Molecular Devices, Sunnyvale, CA) and analyzed using GenePix Pro 5.0 imaging software (Molecular Devices) and Acue 2 software (Mitsui Knowledge Industry, Tokyo, Japan). The signal ratios were normalized in each sample to make the mean signal ratios of all BAC clones 1.0. The reproducibility of BAMCA data was confirmed in representative samples by the duplicate study (data not shown).

2.3. Statistics

BAC clones whose signal ratios obtained by BAMCA differed significantly between the groups of samples were identified by Wilcoxon test. Survival curves of patient groups were calculated by the Kaplan-Meier method, and the differences were compared using the Log-rank test. The Cox proportional hazards multivariate model was used to examine the prognostic impact of DNA methylation status, surgical margin status (R0 versus R1 or R2) [28] and lymph node metastasis. Differences at 𝑃 < . 0 5 were considered significant.

3. Results

3.1. Genome-Wide DNA Methylation Alterations in Tissue Samples

Figure 1 shows representative examples of scanned array images and scattergrams of the signal ratios (test signal/ reference signal) for a C-sample, a N-sample, and the corresponding T-sample. In all C-samples, the signal ratios of 97% of the BAC clones were between 0.67 and 1.5 (red lines in Figure 1(b)). Therefore, in N- and T-samples, DNA methylation status corresponding to a signal ratio of less than 0.67 and more than 1.5 was defined as DNA hypo- and hypermethylation of each BAC clone relative to C-samples, respectively, as in our previous studies [23, 26, 27]. In N-samples, many BAC clones showed DNA hypo- or hypermethylation (Figure 1(b)). In T-samples, more BAC clones showed DNA hypo- or hypermethylation, and the degree of DNA hypo- or hypermethylation, that is, deviation of the signal ratio from 0.67 or 1.5, was increased in comparison with N-samples (Figure 1(b)).

3.2. Establishment of Criteria for Diagnosis of Ductal Adenocarcinomas of the Pancreas Based on DNA Methylation Profiles

Wilcoxon test ( 𝑃 < . 0 1 ) revealed that the average signal ratios of 331 BAC clones (Supplementary Table SI available at http://dx.doi.org/10.1155/2011/780836) in T-samples differed significantly from those in both C- and N-samples. Figure 2(a) shows scattergrams of the signal ratios for representative examples of the 331 BAC clones: RP11-88P10 and RP11-424K7 were able to discriminate T-samples from both C- and N-samples with 100% specificity (the ratio of the number of true negatives to the number of true negatives and false positives) using cutoff values of 0.86 and 1.29 (dotted lines in Figure 2(a)), respectively, (specificity was calculated as the ratio of the number of C- and N-samples showing signal ratios of 0.86 or more than 0.86 and 1.29 or less than 1.29 relative to the total number of C- and N-samples, resp.). The cutoff values of the signal ratios and sensitivities (the ratios of the number of true positives to the number of true positives and false negatives) of 12 BAC clones for which such discrimination was performed with 100% specificity are shown in Table 2. Genes located on the 12 BAC clones are summarized in Supplementary Table SII. A histogram showing the number of BAC clones satisfying the criteria listed in Table 2 in 8 C-samples, 17 N-samples, and 46 T-samples in the learning cohort is shown in Figure 2(b). Based on this histogram, we finally established the following criteria: when tissue samples satisfied the criteria in Table 2 for 1 or more BAC clones, they were judged to be ductal adenocarcinomas, and when tissue samples did not satisfy the criteria for any BAC clone, they were judged not to be ductal adenocarcinomas. Based on these criteria, both the sensitivity and specificity for diagnosis of T-samples in the learning cohort as ductal adenocarcinomas were 100% (sensitivity was calculated as the ratio of the number of T-samples satisfying the criteria in Table 2 for 1 or more BAC clones to the total number of T-samples, and specificity was calculated as the ratio of the number of C- and N-samples not satisfying the criteria in Table 2 for any BAC clone relative to the total number of C- and N-samples).

To confirm these criteria, 68 additional tissue samples in the validation cohort were analyzed. Forty-five samples satisfying the criteria listed in Table 2 for 1 or more BAC clones were all T-samples, and the other 23 samples not satisfying the criteria listed in Table 2 for any BAC clone were all C- or N-samples (Figure 2(c)). Our criteria enabled diagnosis of T-samples in the validation cohort as ductal adenocarcinomas with 100% sensitivity and specificity.

3.3. Establishment of Criteria for Prognostication of Patients with Ductal Adenocarcinomas of the Pancreas Based on DNA Methylation Profiles

To establish criteria for prognostication, 21 patients who had not undergone adjuvant chemotherapy with gemcitabine in the learning cohort were divided into two groups: 4 patients who had not suffered relapse for more than 4 years after pancreatectomy and 17 patients who had suffered relapse within 18 months after pancreatectomy were defined as the no-relapse group and early-relapse group, respectively. The period covered ranged from 215 to 1,846 days (mean, 823 days). Wilcoxon test ( 𝑃 < . 0 5 ) revealed that the average signal ratios of 64 BAC clones differed significantly between T-samples obtained from the no-relapse group and those from the early-relapse group.

Figure 3(a) shows scattergrams of the signal ratios for representative examples of the 64 BAC clones: RP11-165B11 and RP11-92A14 were able to discriminate T-samples from patients belonging to the early-relapse group from those belonging to the no-relapse group with 100% specificity (the ratio of the number of true negatives to the number of true negatives and false positives) using cutoff values of 1.29 and 1.03 (dotted lines in Figure 3(a)), respectively, (specificity was calculated as the ratio of the number of T-samples from patients belonging to the no-relapse group showing signal ratios of 1.29 or more than 1.29 and 1.03 or less than 1.03 relative to the total number of T-samples from patients belonging to the no-relapse group, resp.). The cutoff values of the signal ratios and sensitivities (the ratios of the number of true positives to the number of true positives and false negatives) of 11 BAC clones for which such discrimination was performed with 100% specificity are shown in Table 3. Genes located on the 11 BAC clones are summarized in Supplementary Table SII. A histogram showing the number of BAC clones satisfying the criteria listed in Table 3 in 4 T-samples from the no-relapse group and 17 T-samples from the early-relapse group in the learning cohort is shown in Figure 3(b). Based on these criteria (2 or more BAC clones versus less than 2 BAC clones listed in Table 3), both the sensitivity and specificity of discrimination of patients belonging to the early-relapse group from those belonging to the no-relapse group in the learning cohort were 100% (sensitivity was calculated as the ratio of the number of T-samples from patients belonging to the early-relapse group satisfying the criteria in Table 3 for 2 or more BAC clones relative to the total number of T-samples from patients belonging to the early-relapse group, and specificity was calculated as the ratio of the number of T-samples from patients belonging to the no-relapse group satisfying the criteria in Table 3 for less than 2 BAC clones relative to the total number of T-samples from patients belonging to the no-relapse group).

To confirm these criteria, 34 additional T-samples obtained from patients who had not undergone adjuvant chemotherapy with gemcitabine after surgery in the validation cohort were analyzed. The period covered ranged from 92 to 2,274 days (mean, 612 days). Both the recurrence-free and overall survival rates of 29 patients satisfying the criteria listed in Table 3 for 2 or more BAC clones were significantly lower than those of 5 patients satisfying the criteria listed in Table 3 for less than 2 BAC clones (Figure 3(c), 𝑃 = . 0 0 4 4 and 𝑃 = . 0 1 4 , resp.).

Moreover, multivariate analysis in all 91 patients with ductal adenocarcinomas revealed that satisfying the criteria listed in Table 3 for 2 or more BAC clones was a prognostic parameter for both recurrence-free and overall survival that was independent of surgical margin positivity (R1 or R2) [28] and lymph node metastasis at the time of surgery, which are known to have a prognostic impact [2933] (Table 4).

4. Discussion

Ductal adenocarcinoma of the pancreas, one of the most lethal of all human cancers, is now a common cause of cancer mortality in the United States and Japan [34]. Due to its aggressive growth behavior with early local spread into the surrounding tissues mostly along neural sheaths, peritoneal dissemination, and liver and lymph node metastasis, the prognosis remains poor. Surgical treatment still provides the only possibility of cure [35]. Although advances in preoperative diagnostic imaging have made it possible to detect tumors at an early stage when they are still resectable, diagnosis using pancreatic biopsy and/or specimens of pancreatic juice is indispensable before surgery. In general, pancreatic biopsy yields only a small amount of tissue, and in pancreatic juice specimens, the cellular morphology is not well preserved due to degeneration. Therefore, molecular diagnosis is advantageous for supporting the histological and/or cytological assessment of such specimens. DNA methylation profiles, which are stably preserved on DNA double strands by covalent bonds, even after degeneration of cellular morphology, may become diagnostic markers in pancreatic biopsy and/or pancreatic juice specimens.

We have previously established diagnostic criteria for cancers of the kidney [26], liver [27] and urinary tract [23] based on genome-wide DNA methylation profiles using the BAC array-based approach, BAMCA, which can assess DNA methylation status not only on promoter regions of specific tumor-related genes but also on genomic regions in which DNA hypomethylation affects chromosomal instability. Moreover, during human carcinogenesis, DNA methylation status is frequently altered in a coordinated manner, through processes such as long-range epigenetic silencing [36], in large chromosome regions. Since BAMCA is suitable for overviewing the DNA methylation tendency of individual large regions among all chromosomes [23, 24], we again employed this method to establish diagnostic criteria for ductal adenocarcinomas of the pancreas.

The results of BAMCA for C-samples reflected the DNA methylation profiles of normal peripheral pancreatic duct epithelia (the origin of ductal adenocarcinomas), acinar cells and islet cells. In N-samples, BAMCA revealed DNA hypo- or hypermethylation on many BAC clones in comparison to C-samples (Figure 1(b)). Microscopic observation of N-samples revealed lymphocytes and fibroblasts associated with various degrees of chronic pancreatitis, which is considered to be one of the precancerous conditions for ductal adencarcinomas [11]. Our previous studies using microdissection and immunohistochemistry revealed accumulation of DNA hypermethylation of tumor-related genes associated with DNMT1 overexpression, even in peripheral pancreatic duct epithelia at the precancerous stage [9, 10]. Therefore, the results of BAMCA for N-samples may reflect the DNA methylation profiles of peripheral pancreatic duct epithelia at the precancerous stage, lymphocytes, fibroblasts, acinar cells, and islet cells. In order to diagnose ductal adenocarcinomas in tissue samples, cancer-specific DNA methylation profiles should be discriminated from those of normal or precancerous peripheral pancreatic duct epithelia, lymphocytes, fibroblasts, acinar cells, and islet cells. Therefore, we identified 12 BAC clones whose DNA methylation status was able to discriminate T-samples from both C- and N-samples.

In both the learning and validation cohorts, the criteria combining the 12 BAC clones were able to diagnose T-samples as ductal adenocarcinomas with a sensitivity and specificity of 100%. Our criteria may be advantageous for supporting the histological diagnosis of tiny tissue samples obtained by pancreatic biopsy. Discrimination of cancer cells from exfoliated noncancerous epithelial cells and lymphocytes using the 12 BAC clones may be applicable for diagnosis using specimens of pancreatic juice. Development of methodology for assessing DNA methylation status on the 12 BAC clones in fewer cells may be more advantageous for clinical application, as we have already established a methodology for quantification of DNA methylation levels on specific CpG sites in a very small quantity of genomic DNA for estimation of carcinogenetic risk in patients with chronic liver diseases (unpublished data). Development of this methodology means that if DNA methylation alterations on the 12 BAC clones are not observed in circulating blood cells, our criteria may become applicable for noninvasive diagnosis of pancreatic cancers based on serum markers that differ from the widely used carbohydrate antigen 19-9, whose serum levels are also elevated in patients with chronic pancreatitis [37].

Even when resection with curative intent is performed for patients with pancreatic cancers, the rate of disease recurrence is high and the survival rate after surgery is poor. As surgical resection alone has limitations, development of nonsurgical treatments, including adjuvant therapy, is needed in order to improve the prognosis of patients with pancreatic cancers. Although previous studies have suggested the efficacy of adjuvant chemotherapy [38], it should be carried out carefully, paying close attention to adverse reactions [39]. In order to help decide the indications for such adjuvant chemotherapy after surgery, prognostic indicators should be explored. The criteria listed in Table 3 were able to discriminate the early-relapse group from the no-relapse group with 100% sensitivity and specificity in the learning cohort. Significant correlation between DNA methylation status on the 11 BAC clones and the recurrence-free and overall survival rates of patients with ductal adenocarcinomas in the validation cohort validated the criteria. Multivariate analysis revealed that our criteria were able to predict recurrence-free and overall patient outcome independently of parameters that had been reported to be significantly prognostic in many previous studies, such as surgical margin positivity (R1 or R2) [28] and lymph node metastasis. Therefore, prognostication based on our criteria may be promising for supportive use during followup after surgical resection in patients with ductal adenocarcinomas of the pancreas. Since histological heterogeneity is frequently observed even in a ductal adenocarcinoma of the pancreas from a single patient, the consistency of BAMCA data for multiple T-samples obtained from a single tumor should be carefully confirmed in a prospective validation study before clinical application of the prognostic criteria.

5. Conclusions

BAMCA revealed genome-wide DNA methylation alterations in ductal adenocarcinomas of the pancreas. Criteria combining the DNA methylation status on 12 BAC clones were able to discriminate T-samples from both C- and N-samples and to diagnose T-samples as ductal adenocarcinomas, with 100% sensitivity and specificity in both the learning and validation cohorts. Satisfying the criteria using 11 BAC clones was able to predict the recurrence-free and overall survival of patients with ductal adenocarcinomas independently of surgical margin positivity (R1 or R2) [28] and lymph node metastasis. Genome-wide DNA methylation profiling may provide optimal diagnostic markers for pancreatic cancers and prognostic indicators for affected patients.

Abbreviation

BAC:Bacterial artificial chromosome
BAMCA:BAC array-based methylated CpG island amplification
DNMT:DNA methyltransferase.

Acknowledgment

This study was supported by a Grant-in-Aid for the Third Term Comprehensive 10-Year Strategy for Cancer Control from the Ministry of Health, Labor and Welfare of Japan, a Grant-in-Aid for Cancer Research from the Ministry of Health, Labor and Welfare of Japan, a Grant from the New Energy and Industrial Technology Development Organization (NEDO), and the Program for Promotion of Fundamental Studies in Health Sciences of the National Institute of Biomedical Innovation (NiBio).

Conflict of Interests

There is no potential conflict of interests to disclose.

Supplementary Materials

Wilcoxon test (P < .01) revealed that the average signal ratios of the 331 BAC clones in samples of cancerous tissue (T) differed significantly from those in both samples of normal pancreatic tissue obtained from patients without ductal adenocarcinomas (C) and samples of noncancerous pancreatic tissue obtained from patients with ductal adenocarcinomas (N) in the learning cohort. Such 331 BAC clones are listed and their chromosomal loci based on the UCSC Genome Browser (http://genome.ucsc.edu/cgi-bin/hgGateway) are summarized in Supplementary Table SI.

In 12 BAC clones, using cutoff values of the signal ratios shown in Table 2, samples of cancerous tissue were discriminated from both samples of normal pancreatic tissue obtained from patients without ductal adenocarcinomas and samples of noncancerous pancreatic tissue obtained from patients with ductal adenocarcinomas in the learning cohort with 100% specificity. Genes located on the 12 BAC clones are summarized in Supplementary Table SII (a). In 11 BAC clones, using cutoff values of the signal ratios shown in Table 3, samples of cancerous tissue from patients belonging to the early-relapse group were discriminated from those belonging to the no-relapse group in the learning cohort with 100% specificity. Genes located on the 11 BAC clones are summarized in Supplementary Table SII (b).

  1. Supplementary Material