- Open Access
Cell-free DNA promoter hypermethylation in plasma as a diagnostic marker for pancreatic adenocarcinoma
Clinical Epigeneticsvolume 8, Article number: 117 (2016)
Pancreatic cancer has a 5-year survival rate of only 5–7%. Difficulties in detecting pancreatic cancer at early stages results in the high mortality and substantiates the need for additional diagnostic tools. Surgery is the only curative treatment and unfortunately only possible in localized tumours. A diagnostic biomarker for pancreatic cancer will have a major impact on patient survival by facilitating early detection and the possibility for curative treatment. DNA promoter hypermethylation is a mechanism of early carcinogenesis, which can cause inactivation of tumour suppressor genes. The aim of this study was to examine promoter hypermethylation in a panel of selected genes from cell-free DNA, as a diagnostic marker for pancreatic adenocarcinoma.
Patients with suspected or biopsy-verified pancreatic cancer were included prospectively and consecutively. Patients with chronic/acute pancreatitis were included as additional benign control groups. Based on an optimized accelerated bisulfite treatment protocol, methylation-specific PCR of a 28 gene panel was performed on plasma samples. A diagnostic prediction model was developed by multivariable logistic regression analysis using backward stepwise elimination.
Patients with pancreatic adenocarcinoma (n = 95), chronic pancreatitis (n = 97) and acute pancreatitis (n = 59) and patients screened, but negative for pancreatic adenocarcinoma (n = 27), were included. The difference in mean number of methylated genes in the cancer group (8.41 (95% CI 7.62–9.20)) vs the total control group (4.74 (95% CI 4.40–5.08)) was highly significant (p < 0.001). A diagnostic prediction model (age >65, BMP3, RASSF1A, BNC1, MESTv2, TFPI2, APC, SFRP1 and SFRP2) had an area under the curve of 0.86 (sensitivity 76%, specificity 83%). The model performance was independent of cancer stage.
Cell-free DNA promoter hypermethylation has the potential to be a diagnostic marker for pancreatic adenocarcinoma and differentiate between malignant and benign pancreatic disease. This study brings us closer to a clinical useful diagnostic marker for pancreatic cancer, which is urgently needed. External validation is, however, required before the test can be applied in the clinic.
Pancreatic cancer is the fourth leading cause of cancer death in the world , with a 5-year survival rate of approximately 5–7% [1, 2]. The only curative treatment is complete tumour resection. Unfortunately, only 10–20% of patients receive treatment with the intend to cure. Despite surgery, 50% of patients experience recurrence . Difficulties in detecting the disease at an early stage results in high mortality. This is mainly due to lacking or non-specific symptoms, which are also related to chronic pancreatitis, an essential differential diagnosis and a known risk factor for pancreatic cancer [3, 4]. Often, several complex or invasive techniques such as PET scan (positron emission tomography), CT scan (computed tomography), endoscopic or laparoscopic ultrasound and ERCP (endoscopic retrograde cholangiopancreatography) are needed for the diagnosis and many patients also need a histological evaluation. However, the differentiation between malignant and benign pancreatic disease can be difficult, and even surgery may be needed to establish a definite diagnosis. The only useful biomarker is CA-19-9, which is unspecific as patients with chronic pancreatitis and particularly benign biliary obstruction tend to express high levels of CA-19-9. Moreover, 10% of the population lack the ability to produce CA-19-9, making its utility less apparent [5–7]. It would be a major advance for the patients if a blood-based diagnostic marker was available.
During the development of pancreatic cancer, genetic and epigenetic changes take place. Epigenetic modifications occur at a genomic level, which does not change the DNA sequence. Epigenetic modifications change the DNA conformation and therefore the gene expression. DNA hypermethylation is an epigenetic phenomenon, where a methyl (CH3) residue is added to cytosines preceding guanosines (CpGs) [8–11]. Hypermethylation in the promoter region results in gene silencing, which may be associated with cancer formation [8, 9, 12, 13].
Cancer cells may release cell-free DNA into the blood [14, 15]. DNA hypermethylation can be detected in cell-free DNA in plasma and serum and is potentially tumour specific and useable as blood-based diagnostic markers for pancreatic cancer [14–16].
Thus far, only a few studies with small numbers of patients have evaluated cell-free DNA hypermethylation as a blood-based marker for pancreatic cancer, testing the methylation status of only a single gene or small gene panel . These data have shown a significant difference in DNA hypermethylation between patients with pancreatic cancer and healthy controls [4, 17]. However, the studies had difficulties in differentiating between malignant and benign pancreatic disease . None of the previously examined genes have the potential to serve as an individual diagnostic marker . When developing and testing a biomarker for pancreatic cancer, inclusion of relevant control groups with benign pancreatic disease is very important to enable differentiation of pancreatic cancer-specific hypermethylation and hypermethylation related to pancreatic disease in general .
The aim of this study was to test (by methylation-specific polymerase chain reaction (PCR)) cell-free DNA promoter hypermethylation of a panel of 28 genes as a blood-based diagnostic marker for pancreatic adenocarcinoma, including clinical relevant control groups of patients with benign pancreatic disease.
This study was conducted as a prospective observational cohort study of patients with suspected or biopsy-verified pancreatic cancer admitted to the Department of Gastrointestinal Surgery, Aalborg University Hospital between February 2008 and February 2011 . Additional benign control groups were patients with chronic pancreatitis treated at the hospital or at the outpatient clinic at Aalborg University Hospital between August 2013 and August 2014 and patients admitted with acute pancreatitis at the Department of Gastrointestinal Surgery, Aalborg University Hospital or the Department of General Surgery, Hospital of Vendsyssel between November 2013 and May 2015.
The study was approved by the Research Ethics Committee for the North Denmark Region (N-2013037) and registered in ClinicalTrials.gov (NCT02079363). All participants gave written informed consent.
Consecutive patients with suspected or biopsy-verified upper gastrointestinal cancer were included prospectively in a study on gastrointestinal cancer and venous thromboembolism . Patients had blood drawn on admission before diagnostic work-up and before any treatment. Patients were divided into the following groups (Fig. 1a). Only patients with pancreatic adenocarcinoma (cancer group) and patients screened, but negative for upper gastrointestinal cancer (control group 1), were included in this study. Patients with chronic pancreatitis (control group 2) had blood drawn during hospitalization or at a scheduled visit in the outpatient clinic. Patients diagnosed with acute pancreatitis (control group 3) were enrolled during the first three days of hospitalization. Patients with chronic pancreatitis and acute pancreatitis were excluded if they had previous cancer history or ongoing anticoagulant therapy.
Blood sampling and analytical method
Blood samples were obtained by skilled technicians, using peripheral venipuncture from an antebrachial vein, according to the guidelines recommended by the European Concerted Action on Thrombosis . EDTA plasma samples for methylation analysis were centrifuged for 20 min at 4000 rpm at 4 °C and stored within 2 h after sampling in a biobank at −80 °C until methylation analysis.
All methylation analyses were performed by a single expert laboratory scientist.
Extraction and deamination
Plasma nucleic acids were extracted using the EasyMag platform (Biomerieux) according to manufacturer’s instruction. Approximately 500-μl EDTA plasma was used for the extraction, and purified nucleic acids were eluted in 35-μl elution buffer (Biomerieux). Five microliters was used for quantitation of extracted DNA, and the remainder was deaminated as previously described by our group . In brief, 30-μl DNA extract was mixed with 60-μl deamination solution and deaminated for 10 min at 90 °C, followed by purification using EasyMag and elution in 25 μl 10 mM KOH .
First round of PCR amplification
In order to amplify the amount of deaminated DNA of interest, a first round of PCR was conducted using a mix of outer methylation-specific primers (Additional file 1) for all promoter regions tested. The reaction buffer for each sample consisted of 25 μl containing PCR stock, 13 μM MgCl2, 0.6 mM dNTP, 250 nM of each outer primer (Additional file 1), 1.5 U Taq polymerase (Bioline) and 0.3 U UNG (Invitrogen). The first round reaction mix was distributed to individual 200-μl PCR tubes and incubated for 5 min at 37 °C (UNG activity), followed by incubation at 95 °C for 5 min and cooling to room temperature. Twenty-five microliters of purified deamination product was added to each tube containing the first round reaction mix. PCR was performed for 20 cycles at 92 °C for 15 s, 55 °C for 30 s and 72 °C for 30 s.
Second round of PCR
Ten microliters of mix containing 0.4 μM inner methylation-specific primers and methylation-specific probes (Additional file 1) was distributed in 30 individual wells in a 96-well PCR plate. Ten microliters of first round PCR product was added to 710 μl of reaction mix containing PCR stock, 250 μM dNTP, 10 μM MgCl2, and 15 U Taq polymerase (Bioline). Twenty microliters of the reaction mix was added to each of the 30 wells containing primers and probes. Real-time PCR was carried out for 45 cycles at 94 °C for 15 s, 55 °C for 30 s (annealing and detection) and 72 °C for 30 s.
Twenty-eight genes were selected for methylation analysis (Additional file 2). The genes were primarily selected based on a literature review performed by our group prior to this study . The genes selected for the panel had previously been detected as hypermethylated in either cell-free DNA in plasma or serum, pancreatic juice or tumour tissue from patients with pancreatic cancer and in addition unmethylated in samples from healthy individuals. Few additional genes were chosen based on a pilot study on cell-free DNA hypermethylation in colorectal adenocarcinoma (unpublished data).
Primer and probe design
All primers and probes were designed using Beacon Designer® (PREMIER Biosoft International, Palo Alto, CA) software and evaluated to be hypermethylation specific by MethPrimer® (The Li Lab, Peking, China) . Primers were designed to be rich in CpGs and to be located up-stream of exon one, which was interpreted as the promoter regions of the genes. The aim was to design PCR products with a length less than 140–150 base pairs, because the cell-free DNA fragments most likely have a length of 160 base pairs consistent with nucleosomal DNA size . The primers and probes were designed and optimized for the present study; however, effort was made to design primers for previously tested promoter sequences (Additional file 1).
Hemi-methylated MEST transcript variant 1 was used as reference gene in both the first and second round PCR.
The primary aim was to establish a prediction model for pancreatic adenocarcinoma, enabling differentiation of pancreatic adenocarcinoma patients and a clinical relevant control group of patients screened, but negative for upper gastrointestinal cancer and patients with chronic pancreatitis.
Each gene in the gene panel was analysed as a binary variable (hypermethylated or non-methylated).
Validation of dichotomous data
We calculated the differences between the threshold cycle (Ct) values of the hemimethylated reference gene MEST transcript variant 1 and the Ct values of each gene for which Ct > 0 (∆ Ct). To assess the amount of information lost in the dichotomization, histograms of ∆ Ct for the cancer group and control group 1 combined with control group 2 were produced.
Level of cell-free DNA
We calculated the median level (ng/ml) of cell-free DNA for each group. Nonparametric Wilcoxon rank sum test was used for comparison of the cancer group and the benign control groups.
The total number of hypermethylated genes was calculated for each patient. The Kendall’s rank test was used for correlation analysis of total number of hypermethylated genes and level of cell-free DNA.
The methylation frequency of each gene and the (exact) 95% confidence interval (CI) were calculated for each group. The mean number of hypermethylated genes in each group and the 95% CI was calculated. The means were compared as numerical data with the nonparametric Wilcoxon rank sum test. p values less than 0.05 were considered statistically significant.
Prediction model development
Screening of each individual variable as a diagnostic marker for pancreatic adenocarcinoma: Logistic regression was performed separately for each gene in the gene panel and for smoking status, gender and patient age >65. The p value and the area under the receiver operating characteristic curve (AUC) were calculated.
The selection of variables: Variables having a p value less than 0.2 were selected for further analysis.
Model selection: Stepwise backward elimination in logistic regression models was performed to select the relevant variables using 0.05 as the significance level for removal from the model. In the backward elimination algorithm, variables were eliminated one by one to identify the optimal combination of variables representing the highest predictive power. The least significant variable in the variable combination was eliminated in the stepwise procedure. For each intermediate model, the AUC value was calculated.
Determination of the best model: The decision was based on the model complexity combined with the model performance according to the AUC.
Interactions between the variables: The significance of interactions between all pairs of variables was assessed in the final model. Interactions with a p value less than 0.01 were considered statistically significant.
Validation: To account for optimism in the internal validation of discriminative model performance (measured by the AUC), “leave pair out cross validation” was used . For the calibration performance, Hosmer–Lemeshow test was performed.
Probability score: For each patient, a probability score was calculated.
All data were analysed using Stata 14.0 software [StataCorp LP, TX].
All authors had full access to the study data and had reviewed and approved the final manuscript.
Ninety-five patients with confirmed pancreatic adenocarcinoma were included in the study (Fig. 1a). After diagnostic work-up (gastroscopy, endoscopic ultrasound, magnetic resonance (MR) or CT scan), 35 patients without evidence of malignancy were categorized as patients screened but negative for pancreatic adenocarcinoma (control group 1). Eight patients were subsequently excluded from this group (Fig. 1a). Two additional groups of control patients with benign pancreatic disease were included. Overall, 103 patients with chronic pancreatitis (control group 2) and 62 patients with acute pancreatitis (control group 3) were included. Subsequently, six patients from control group 2 and three patients from control group 3 were excluded (Fig. 1b, c). Descriptive data of the four groups are shown in Table 1.
Validation of the dichotomous data
There was no clear difference in ∆ Ct between the cancer group and control group 1 combined with control group 2, which indicated that no significant amount of information was lost by dichotomizing the genes as hypermethylated or non-methylated genes regardless of the observed Ct value. Additional file 3a, b lists the distribution of Ct values (0, 0–25, 25–30 and >30) for each gene in within patient group and illustrates a slightly difference in Ct values between the groups, with a tendency towards Ct values in the cancer group being lower compared to the benign control groups. However, due to limited power, the effect of this difference could not be evaluated in the multivariable logistic regression model; consequently, we treated hypermethylation as a dichotomised variable.
Level of cell-free DNA
Patients with pancreatic adenocarcinoma had a significant higher median level of cell-free DNA (11.60 ng/ml (range 0.60–957.17)) compared to control group 1 with 6.17 ng/ml (range 1.06–48.43), control group 2 with 2.18 ng/ml (0.11–115.44) and 4.09 ng/ml (range 0.65–62.42) for control group 3 (Additional file 4). In addition, the correlation between level of cell-free DNA and number of hypermethylated genes was statistically significant with a Kendall’s τ of 0.34 (Additional file 5).
The hypermethylation profile for each patient is illustrated on the heat map plot in Additional file 6. The methylation frequency of each gene is presented in Table 2. The mean number of methylated genes of the whole gene panel (28 genes) was 8.41 (95% CI 7.62–9.20) for the cancer group compared to 4.34 (95% CI 3.85–4.83) for patients with chronic pancreatitis (control group 2), 4.89 (95% CI 4.07–5.71) for patients screened, but negative for pancreatic cancer (control group 1) and 5.34 (95% CI 4.76–5.91) for patients with acute pancreatitis (control group 3). The difference between the cancer group and the three benign control groups was highly statistically significant (Table 3).
Prediction model development
In the following analyses, we chose to combine control group 1 and 2, as the combined group has symptoms resembling those of pancreatic cancer, which makes a biomarker to distinguish these from pancreatic cancer of utmost clinical relevance. For the remainder of the analysis, patients with acute pancreatitis were excluded because a clinical picture of acute inflammation is rarely seen in pancreatic cancer.
There was a highly significant difference (p < 0.001) between the cancer group and control groups 1 + 2 with regard to hypermethylation frequency of ten genes (APC, BMP3, BNC1, MESTv2, NPTX2, RASSF1A, SFRP1, SST, TFPI2 and TAC1) (Table 4) and significant difference (p < 0.05) in seven other genes (ALX4, ESR1, HIC1, RARB, SFRP2, SEPT9v2 and WNT5A) (Table 4). VIM and PENK could not be evaluated by logistic regression, as none of the patients in the control group had hypermethylation of these two genes; however, chi-square test found significant difference between the cancer group and the control groups 1 + 2. Despite that, VIM and PENK were excluded from the following analysis because only very few cancer patients had VIM or PENK hypermethylation. (Tables 2 and 4). There was no significant difference in gender; consequently, this variable was excluded from the subsequent analysis. Smoking, however, was a preventive factor for cancer when comparing patients with pancreatic cancer and patients with chronic pancreatitis. Smoking was therefore excluded from the model because it is a known risk factor for cancer. By stratifying the patients into groups according to age (>65 years old, ≤65 years old), a statistically significant difference was found between the cancer group and control groups 1 + 2. Consequently, patient age >65 years old was included as a covariate in the multivariable logistic regression analysis.
All genes with an individual p value below 0.20 (20 genes out of 28 examined genes) and patient age >65 were included in the multivariable logistic regression model. Backward stepwise selection was performed. Figure 2 illustrates the stepwise elimination of variables from the model and the corresponding AUC. The initial model (model 1) with 20 genes had an AUC of 0.87 (Fig. 2). Removing the 12 least significant genes from the model and leaving eight genes (model 13; age >65, BMP3, RASSF1A, BNC1, MESTv2, TFPI2, APC, SFRP1 and SFRP2) resulted in an AUC of 0.86 (95% CI 0.81–0.91) (Figs. 2 and 3a). The mean probability for having pancreatic adenocarcinoma was 0.67 (0.61–0.72) for cancer patients and 0.26 (0.22–0.29) for the control groups 1 + 3. Model 13 was determined as the model with the best performance (probability cut point of 0.50; sensitivity 76% and specificity 83%). There were no statistically significant interactions between variables in model 13. The model was well calibrated (p = 0.40) and had an estimated optimism in AUC of 0.03.
Forty patients had stage I or II tumours. Model 13 had an apparent AUC of 0.86 (95% CI 0.79–0.92) for stage I/II tumours (probability cut point of 0.50; sensitivity 73% and specificity 83%) (Fig. 3b) with an optimism in AUC of 0.06.
We examined cell-free DNA promoter hypermethylation of 28 genes in the plasma of a large cohort of patients with pancreatic adenocarcinoma and compared it to different clinical relevant control groups. We designed the gene-panel primary based on our literature review addressing genes aberrantly methylated in pancreatic adenocarcinoma . This approach was used, to evaluate the overall diagnostic performance of genes which previously had been examined separately as diagnostic markers for pancreatic cancer. The panel was composed of genes previously detected as hypermethylated in plasma/serum and tumour tissue in relation to pancreatic adenocarcinoma (BNC1 [17, 24, 25], NPTX2 [4, 24–26], PENK [4, 14, 25], CDKN2A [4, 26, 27], RASSF1A [4, 24, 27], SFRP1 (SARP2) [4, 25], APC [24, 27], BRCA1 [28, 29], CDKN2B [28, 30], ESR1 [25, 28], MGMT [24, 28], MLH1 [28, 31] and RARB [28, 32]), genes earlier found to be hypermethylated in pancreatic juice or tumour tissue from patients with pancreatic adenocarcinoma (BMP3 , EYA2 , GSTP1 , HIC1 [25, 33], SFRP2 , TFPI2 , VIM , NEUROG1 [24, 25], TAC1 [24, 25], CHFR  and WNT5a [24, 25]) and genes found based on a pilot study on cell-free DNA hypermethylation in colorectal adenocarcinoma (ALX4, MESTv2, SEPT9v2 and SST). To our knowledge, this is the first study to examine cell-free DNA hypermethylation in a wide selection of genes by methylation-specific PCR in a large group of patients with either malignant or benign pancreatic disease.
A statistically significant difference in the hypermethylation status in 19 out of the 28 genes was found when comparing pancreatic adenocarcinoma patients and a control group containing patients screened, but negative for pancreatic cancer, as well as in patients with chronic pancreatitis. Cell-free DNA hypermethylation of BMP3, MESTv2, SST, TFPI2, TAC1, ALX4, HIC1, SFRP2, SEPT9v2 and WNT5A has not previously been described in the literature in relation to pancreatic cancer. Yi et al. described BNC1 hypermethylation to have a sensitivity of 79% and a specificity of 89% when comparing pancreatic cancer and healthy individuals . We found BNC1 to be hypermethylated in only 36% of pancreatic cancer patients with a specificity of 94%. Park et al. examined hypermethylation of a small gene panel (NPTX2, RASSF1A, SFRP1, UCHL1, PENK and CDKN2A) by methylation-specific PCR . The gene panel could differentiate pancreatic cancer from healthy controls; however, it was not able to discriminate benign and malignant pancreatic disease.
Our study shows that cell-free DNA hypermethylation is detectable in both malignant and benign pancreatic disease. However, patients with pancreatic adenocarcinoma have a higher level of hypermethylated genes in plasma-derived cell-free DNA. Consistent with previous studies, our gene panel did not demonstrate a single gene, which could be used as an individual diagnostic marker for pancreatic cancer. This result suggests that a larger gene panel is needed to achieve sufficient accuracy . We developed a diagnostic prediction model (age >65, BMP3, RASSF1A, BNC1, MESTv2, TFPI2, APC, SFRP1 and SFRP2), which was able to differentiate between pancreatic adenocarcinoma and a large control group of great clinical relevance. The control group included patients with chronic pancreatitis or patients referred to the hospital with symptoms of pancreatic cancer. The AUC was high, and the predictive value of our model is superior to the predictive value of CA-19-9, which currently is the only blood-based biomarker for pancreatic cancer. Particularly keeping in mind that CA-19-9 is highly dependent on the Lewis blood group status of the patients. Only Lea+b− or Lea-b+ individuals are able to express CA-19-9 but not Lea−b− individuals, which represent 5–10% of the Caucasian population . In a recent study, CA-19-9 could differentiate patients with stage I–II pancreatic cancer from patients with chronic pancreatitis with an AUC of 0.77 (sensitivity of 53% and a specificity of 92%) and pancreatic cancer patients from patients with benign biliary obstruction with an AUC of only 0.65 . Our study included patients with stage I–IV pancreatic adenocarcinoma. It is most important to diagnose patients with stage I and II disease as early detection at this stage of the disease has the potential to improve the outcome of surgery. We tested our model on stage I and II disease and found an AUC of 0.86. This finding shows that the performance of the prediction model is independent of the cancer stage. DNA hypermethylation is detectable in plasma even in an early stage of the disease and thereby potentially usable as an early blood-based diagnostic marker.
In order to further differentiate DNA hypermethylation related to malignant and benign pancreatic disease, patients with acute pancreatitis were included. The aim was to achieve a more basic understanding of hypermethylated DNA during the course of an acute pancreatic inflammatory reaction, which has not been described earlier in literature. Our study shows that DNA hypermethylation takes place during pathological conditions in the pancreas including acute inflammation. However, the changes are more pronounced in patients with pancreatic adenocarcinoma.
Our study has some limitations. The study was exploratory, showing training data only, which is known to produce an overestimation of the test performance due to overfitting. Validation of the results in an independent cohort is needed to substantiate the results.
Patients were not matched according to age, which one should be aware of because epigenetic changes can be a part of ageing . To address this problem, we incorporated age as a covariate in our prediction model.
In addition, comparison of the performance of our prediction model to CA-19-9 would have been relevant. Unfortunately, this was impossible, as CA-19-9 was not available on two thirds of the patients as this test was first implemented in 2010 at our department.
The difference in sensitivity of the genes analysed in our study and the sensitivity of genes examined in previous studies by others might be due to the use of different primer sequences. Several methods are described for methylation analysis which furthermore makes inter-study comparison difficult .
We performed methylation-specific PCR, which is a quantitative method using hemimethylated MEST transcript variant 1 as a reference gene . However, our study lacked sufficient power to conduct a quantitative analysis. Therefore, we analysed hypermethylation as a binary variable, which unfortunately results in loss of the quantitative information.
At the end of the analyses, we discovered that the use of UNG (Invitrogen) had a tendency to lower the sensitivity compared to the use of COD UNG (ArcticZymes). All our samples are analysed using UNG (Invitrogen) because it was not possible to repeat all analyses with COD UNG (ArcticZymes) due to the lack of sample material.
However, the study also has several strengths. We tested cell-free DNA hypermethylation of a broad gene panel in the plasma from a large group of patients with pancreatic adenocarcinoma, all included prospectively and consecutively, before the diagnostic workup and before any treatment.
In addition, we included a large utmost relevant group of control patients with either benign pancreatic disease or with symptoms mimicking pancreatic cancer, which all are patients clinically hard to differentiate from patients with pancreatic cancer.
We performed methylation-specific PCR based on an optimized bisulfite treatment protocol . This method has several advantages, due to a high recovery from samples with minute amounts of DNA (<0.01 ng/ml) and a rapid deamination of DNA (less than 2 h) .
We developed a diagnostic prediction model for pancreatic adenocarcinoma with a high performance, independent of cancer stage. In addition, the diagnostic prediction model only had a modest optimism in performance by intern validation.
Diagnostic biomarkers for pancreatic cancer are lacking. We developed a diagnostic test, which has the great advantage of being blood based and thereby minimally invasive. In general, blood-based markers are of great benefit to the patients compared to tissue-based markers, as the latter entail a risk of complications. Furthermore, due to the deep location of the pancreas in the upper abdomen, biopsies may be difficult to obtain, why blood-based markers are of utmost importance regarding pancreatic disease.
Our study demonstrates statistically significant differences in cell-free DNA hypermethylation of several genes between malignant and benign pancreatic diseases. Patients with pancreatic adenocarcinoma have a highly significant number of hypermethylated genes compared to patients with benign pancreatic diseases. A panel of hypermethylated genes (BMP3, RASSF1A, BNC1, MESTv2, TFPI2, APC, SFRP1 and SFRP2) are able to differentiate between patients with pancreatic adenocarcinoma and a most relevant control group. Based on our study, alterations in cell-free DNA hypermethylation have the potential of serving as blood-based biomarkers for the diagnosis of pancreatic adenocarcinoma. External validation is however required before the biomarker can be applied in daily clinical practice.
American Cancer Society: Pancreatic cancer survival rates, by stage [Online material]. http://www.cancer.org/cancer/pancreaticcancer/detailedguide/pancreatic-cancer-survival-rates Accessed Nov. 2016.
Jemal A, Bray F, Ferlay J, Al E. Global Cancer Statistics. CA Cancer J CLIN. 2011;61:69–90.
Michl P, Gress TM. Current concepts and novel targets in advanced pancreatic cancer. Gut. 2012;62(figure 1):317–26.
Park JW, Baek IH, Kim YT. Preliminary study analyzing the methylated genes in the plasma of patients with pancreatic cancer. Scand J Surg. 2012;101:38–44.
Haab BB, Huang Y, Balasenthil S, Partyka K, Tang H, Anderson M, Allen P, Sasson A, Zeh H, Kaul K, Kletter D, Ge S, Bern M, Kwon R, Blasutig I, Srivastava S, Frazier ML, Sen S, Hollingsworth MA, Rinaudo JA, Killary AM, Brand RE. Definitive characterization of CA 19-9 in resectable pancreatic cancer using a reference set of serum and plasma specimens. PLoS One. 2015;10:e0139049.
Hartwig W, Strobel O, Hinz U, Fritz S, Hackert T, Roth C, Büchler MW, Werner J. CA19-9 in potentially resectable pancreatic cancer: perspective to adjust surgical and perioperative therapy. Ann Surg Oncol. 2013;7:2188–96.
Kim J-E, Lee KT, Lee JK, Paik SW, Rhee JC, Choi KW. Clinical usefulness of carbohydrate antigen 19-9 as a screening test for pancreatic cancer in an asymptomatic population. J Gastroenterol Hepatol. 2004;19:182–6.
Delpu Y, Hanoun N, Lulka H, Sicard F, Selves J, Buscail L, Torrisani J, Cordelier P. Genetic and epigenetic alterations in pancreatic carcinogenesis. Curr Genomics. 2011;12:15–24.
Lomberk GA. Epigenetic silencing of tumor suppressor genes in pancreatic cancer. J Gastrointest Cancer. 2011;42:93–9.
Lomberk G, Mathison AJ, Grzenda A, Urrutia R. The sunset of somatic genetics and the dawn of epigenetics: a new frontier in pancreatic cancer research. Curr Opin Gastroenterol. 2008;24:597–602.
Mulero-Navarro S, Esteller M. Epigenetic biomarkers for human cancer: the time is now. Crit Rev Oncol Hematol. 2008;68:1–11.
Costa FF. Epigenomics in cancer management. Cancer Manag Res. 2010;2:255–65.
Sebova K, Fridrichova I. Epigenetic tools in potential anticancer therapy. Anticancer Drugs. 2010;21:565–77.
Jiao L, Zhu J, Hassan MM, Evans DB, Abbruzzese JL, Li D. K-ras mutation and p16 and preproenkephalin promoter hypermethylation in plasma DNA of pancreatic cancer patients: in relation to cigarette smoking. Pancreas. 2007;34:55–62.
Esteller M, Rosell R, Sidransky D, Baylin SB. Detection of aberrant promoter hypermethylation of tumor suppressor genes in serum DNA from non-small cell lung cancer patients. Cancer Res. 1999;59:67–70.
Henriksen SD, Madsen PH, Krarup H, Thorlacius-Ussing O. DNA hypermethylation as a blood-based marker for pancreatic cancer: a literature review. Pancreas. 2015;44:1036–45.
Yi JM, Guzzetta A a, Bailey VJ, Downing SR, Van Neste L, Chiappinelli KB, Keeley BP, Stark A, Herrera A, Wolfgang C, Pappou EP, Iacobuzio-Donahue C a, Goggins MG, Herman JG, Wang T-H, Baylin SB, Ahuja N. Novel methylation biomarker panel for the early detection of pancreatic cancer. Clin Cancer Res. 2013;19:6544–55.
Larsen AC, Dabrowski T, Frøkjær JB, Fisker RV, Iyer VV, Møller BK, Kristensen SR, Thorlacius-Ussing O. Prevalence of venous thromboembolism at diagnosis of upper gastrointestinal cancer. Br J Surg. 2014;101:246–53.
Jespersen J, Bertina R, Haverkate F. Laboratory techniques in thrombosis—a manual: second revised edition of the ECAT assay procedures. The Netherlands: Kluwer Academic Publishers, Dordrecht; 1999.
Pedersen IS, Krarup HB, Thorlacius-Ussing O, Madsen PH. High recovery of cell-free methylated DNA based on a rapid bisulfite-treatment protocol. BMC Mol Biol. 2012;13:12.
Li L-C, Dahiya R. MethPrimer: designing primers for methylation PCRs. Bioinformatics. 2002;18:1427–31.
Mouliere F, Rosenfeld N. Circulating tumor-derived DNA is shorter than somatic DNA in plasma. Proc Natl Acad Sci. 2015;112:3178–9.
Smith GCS, Seaman SR, Wood AM, Royston P, White IR. Correcting for optimistic prediction in small data sets. Am J Epidemiol. 2014;180:318–24.
Zhao Y, Sun J, Zhang H, Guo S, Gu J, Wang W, Tang N, Zhou X, Yu J. High-frequency aberrantly methylated targets in pancreatic adenocarcinoma identified via global DNA methylation analysis using methylCap-seq. Clin Epigenetics. 2014;6:18.
Vincent A, Omura N, Hong SM-M, 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:4341–54.
Park JK, Ryu JK, Yoon WJ, Lee SH, Lee GY, Jeong KS-S, Kim YT-T, Yoon YB. The role of quantitative NPTX2 hypermethylation as a novel serum diagnostic marker in pancreatic cancer. Pancreas. 2012;41:95–101.
Kawasaki H, Igawa E, Kohosozawa R, Kobayashi M, Nishiko R, Abe H. Detection of aberrant methylation of tumor suppressor genes in plasma from cancer patients. Pers Med Universe. 2013;2:20–4.
Liggett T, Melnikov A, Yi Q-LL, Replogle C, Brand R, Kaul K, Talamonti M, Abrams RA, Levenson V. Differential methylation of cell-free circulating DNA among patients with pancreatic cancer versus chronic pancreatitis. Cancer. 2010;116:1674–80.
Peng DF-F, Kanai Y, Sawada M, Ushijima S, Hiraoka N, Kitazawa S, Hirohashi S. DNA methylation of multiple tumor-related genes in association with overexpression of DNA methyltransferase 1 (DNMT1) during multistage carcinogenesis of the pancreas. Carcinogenesis. 2006;27:1160–8.
Li G, Ji Y, Liu C, Li J, Zhou Y. Reduced levels of p15INK4b, p16INK4a, p21cip1 and p27kip1 in pancreatic carcinoma. Mol Med Rep. 2012;5:1106–10.
Li M, Zhao ZW. Clinical implications of mismatched repair gene promoter methylation in pancreatic cancer. Med Oncol. 2012;29:970–6.
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–9.
Zhao G, Qin Q, Zhang J, Liu Y, Deng S, Liu L, Wang B, Tian K, Wang C. Hypermethylation of HIC1 promoter and aberrant expression of HIC1/SIRT1 might contribute to the carcinogenesis of pancreatic cancer. Ann Surg Oncol. 2013;20:S301-S311.
Ballehaninna UK, Chamberlain RS. Serum CA 19-9 as a biomarker for pancreatic cancer—a comprehensive review. Indian J Surg Oncol. 2011;2:88–100.
Sinsheimer JS, Bocklandt S, Lin W, Sehl ME, Sa FJ, Vilain E. Epigenetic predictor of age. PLoS One. 2011;6:1–6.
Kurdyukov S, Bullock M. DNA methylation analysis: choosing the right method. Biology (Basel). 2016;5:3.
We are grateful to all the patients who have participated in this study. We thank the private foundations that supported the study: Medical Advancement of Science Foundation, Speciallæge Heinrick Koops Foundation, Aase and Ejnar Danielsens Foundation, Marie Pedersen and Jensine Heibergs Foundation, Beckett Foundation, Blegdalens Foundation and Resident Foundation at Aalborg University Hospital. We also thank J. Lundtoft, a laboratory technician at the Research Unit of the Department of Gastrointestinal Surgery, Aalborg University Hospital for obtaining the blood samples. We thank the staff at the Mech-Sense, Department of Gastroenterology, Aalborg University Hospital for assistance with the inclusion of patients with chronic pancreatitis and the Department of General Surgery, Hospital of Vendsyssel for assistance with the inclusion of patients with acute pancreatitis.
Funding was provided by the Department of Gastrointestinal Surgery, Aalborg University Hospital; Department of Clinical Medicine, Aalborg University; and Section of Molecular Diagnostics, Clinical Biochemistry, Aalborg University Hospital. Private foundations are as follows: Medical Advancement of Science Foundation, Speciallæge Heinrick Koops Foundation, Aase and Ejnar Danielsens Foundation, Marie Pedersen and Jensine Heibergs Foundation, Beckett Foundation, Blegdalens Foundation and Resident Foundation at Aalborg University Hospital.
The foundations had no influence on the study design, data analysis, data interpretation or manuscript writing. The corresponding author confirms that she had full access to all the data in the study and had final responsibility for the decision to submit for publication.
Availability of data and materials
The datasets analysed during the current study are available from the corresponding author on reasonable request.
OTU proposed the idea of the project. SDH, OTU, HK and PHM performed the literature review and designed the study. ACL included all patients with suspected or biopsy-verified pancreatic cancer. SDH included patients with acute pancreatitis. SDH and AMD included patients with chronic pancreatitis. ACL acquired sample material from patients with pancreatic cancer and was responsible for the database management of this patient group. SDH was responsible the acquisition of sample material and database management of patient with benign pancreatic disease. ISP, HK and PHM developed the optimized method of bisulfite treatment. PHM performed the methylation analyses. MBJ and SDH performed the statistical analysis and data interpretation. SDH drafted the first version of the paper. All authors revised the draft and approved the final version of the paper.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
The study was approved by the Research Ethics Committee for the North Denmark Region (N-2013037) in June 2013. All participants gave written informed consent.
DNA sequences for probes and primers. (DOCX 24 kb)
Characteristics of genes used in the gene panel. (DOCX 23 kb)
a: Distribution of Ct values in the cancer group and control group 1. b: Distribution of Ct values in the control group 2 and control group 3. (ZIP 26 kb)
Level of cell-free DNA. (DOCX 15 kb)
Correlation between level of cell-free DNA (ng/ml) and total number of hypermethylated genes. (DOCX 15 kb)
Heat map of hypermethylated genes in each patient. (XLSX 38 kb)