Pancreas Cancer Incidence and Pancreas Cancer-Associated Mortality Are Low in National Cohort of 7211 Pancreas Cyst Patients.

BACKGROUND AND AIMS: Pancreatic cancer incidence and mortality among patients with pancreas cysts are unclear. The aims of this study are to evaluate incidence of pancreatic cancer and cause-specific mortality among patients with pancreatic cysts using a large national cohort over a long follow-up period.
METHODS: We conducted a retrospective cohort study of US Veterans diagnosed with a pancreatic cyst 1999-2013, based on International Classification of Diseases, 9th edition (ICD9) coding within national Department of Veterans Affairs (VA) data. Pancreatic cancer incidence was ascertained using VA cancer registry data, ICD-9 codes, and the National Death Index, a national centralized database of death records, including cause-specific mortality.
RESULTS: Among 7211 Veterans with pancreatic cysts contributing 31,501 person-years of follow-up (median follow-up 4.4 years), 79 (1.1%) developed pancreatic cancer. A total of 1982 patients (27.5%) died during the study follow-up period. Sixty-three patients (3.2% of deaths; 0.9% of pancreas cyst cohort) died from pancreatic cancer, but the leading causes of death in the cohort were non-pancreatic cancer (n = 498, 25% of deaths) and cardiovascular disease (n = 398, 20% of deaths).
CONCLUSIONS: Pancreas cancer incidence and pancreatic cancer-associated mortality are very low in a large national cohort of VA pancreatic cyst patients with long-term follow-up. Most deaths were from non-pancreas cancers and cardiovascular causes, and only a minority (3.2%) were attributable to pancreas cancer. Given death from pancreas cancer is rare, future research should focus on identifying criteria for selecting individuals at high risk for death from pancreatic cancer for pancreatic cyst surveillance.

Introduction

Pancreatic cysts are commonly seen with an estimated prevalence of 2–15% [1-6]. Furthermore, with the increasing use of cross-sectional imaging, a greater number of incidental pancreatic cysts are being identified in patients without other signs or symptoms of pancreatic disease [7]. Once detected, these cysts can trigger significant anxiety about developing and dying from pancreatic cancer. Pancreatic cysts encompass benign cystic lesions such as pseudocysts and serous cystadenomas and potentially malignant cysts such as intraductal papillary mucinous neoplasms (IPMNs) and mucinous cystic neoplasms (MCNs). However, our ability to accurately determine cyst type, risk for developing cancer, and need for surgery is limited. Further, surgery for pancreatic cysts is associated with significant mortality and morbidity [8].

The incidence of pancreatic cancer among patients with incidentally found pancreatic cysts is postulated to be low but has not been well characterized. Most of the data regarding incidence of pancreatic cancer among patients with pancreatic cysts are based on surgical resection data. The major challenges to apply findings from surgical series include inherent problems with (a) selection bias and (b) an unclear understanding about the long-term outcomes in patients with non-malignant cysts at time of surgery (i.e., would some of these cysts eventually develop cancer?). Published non-surgical series are limited by small size, small numbers of patients developing pancreatic cancer and short-term follow-up resulting in imprecise estimates of pancreatic cancer incidence.

Even less is known about pancreas cancer-associated mortality among patients with pancreatic cysts due to a dearth of population-based studies evaluating the natural history of pancreatic cysts. Additional work is needed to clarify the pancreatic cancer-associated mortality after pancreatic cyst diagnosis and place risk of pancreatic cancer mortality in context of competing mortality risks from other causes.

Thus, the primary aim of this study is to calculate incidence of pancreatic cancer among patients with pancreatic cysts using a large national cohort over a long follow-up period. The secondary aim is to evaluate all-cause mortality and cause-specific mortality, including pancreas cancer-associated mortality among patients with pancreatic cysts.

Methods

Study Design

We conducted a retrospective cohort study of US Veterans with pancreatic cysts utilizing national Department of Veterans Affairs (VA) electronic health record data accessed through the VA Informatics and Computing Infrastructure (VINCI). VINCI allows secure access to the national VA corporate data warehouse (CDW) for research purposes. CDW data include discrete coded fields, such as prior diagnoses, medication history, clinical laboratories, and demographic characteristics, as well as free-text notes for all Veteran healthcare encounters, such as endoscopic ultrasound reports, clinical progress notes, radiology reports, and pathology reports. Furthermore, the complete electronic medical record of any Veteran can be reviewed via VINCI.

Pancreas Cyst Cohort

Our pancreas cyst cohort was created by identifying all adult (age ≥ 18 years) VA patients from January 1, 1999, to December 31, 2013, with International Classification of Diseases, Ninth Revision (ICD-9) code 577.2. Because the 577.2 code is assigned to individuals with pancreatic cyst or pseudocysts, individuals with a concurrent or prior ICD-9 codes for 577.0 (acute pancreatitis) and 577.1 (chronic pancreatitis) were excluded to minimize patients with pancreatic pseudocysts and to avoid patients with potentially symptomatic pancreatic cystic lesions causing pancreatitis. Patients diagnosed with pancreatic cancer within the first year of pancreas cyst diagnosis were excluded to minimize misclassification of patients who had prevalent cancer at time of cyst diagnosis. Patients with shorter than 1 year of follow-up were also excluded.

Pancreas Cancer Ascertainment

Pancreas cancer was ascertained based on meeting one or more of the following criteria: (1) presence of ICD-9 code 157.0, 157.1, 157.2, 157.3, 157.8, or 157.9 for pancreas cancer; (2) death from pancreas cancer based on the National Death Index (NDI); and (3) presence of an International Classification of Diseases for Oncology, Third Edition (ICD-O-3) code for pancreas cancer within the Oncology Domain. The NDI is a centralized database of death records which also contains information on cause-specific mortality and is considered the most accurate death record in the USA [9]. A major advantage of utilizing NDI data is to obtain information for Veterans who obtained care outside of the VA after cancer diagnosis. The Oncology Domain within VINCI contains data abstracted from local VA hospital cancer registrars. As part of the current study, manual chart review was performed for all patients with pancreatic cancer to confirm diagnoses. Patients who did not meet all inclusion criteria on chart review were excluded. Patients identified as having pathology other than adenocarcinoma (such as neuroendocrine tumors) were excluded from the analysis. For patients who were only identified utilizing NDI, the date of death from pancreas cancer was used for pancreas cancer diagnosis date. For three patients who were identified by ICD-9 code only, the diagnosis of pancreas cancer was based on high clinical suspicion by the treating physicians as a biopsy was never pursued/performed. These patients were included in the analysis because pancreas cancer is often found in an advanced state clinically, and biopsy is sometimes not pursued based on clinical considerations such as patient preferences.

All-Cause and Cause-Specific Mortality

All-cause and cause-specific mortality was ascertained using NDI. Cause of death was divided into the following categories: pancreas cancer, non-pancreas cancer, cardiovascular, respiratory, cerebrovascular, hepatic, diabetes, kidney, infection, and other.

Covariates

Covariates (at time of pancreas cyst diagnosis) including age, gender, race/ethnicity, body mass index (BMI), tobacco use, and Charlson Comorbidity Index Score (CCIS) were ascertained using data within the CDW. Diabetes was based on presence of at least one inpatient ICD-9 diagnosis of diabetes, at least 2 outpatient ICD-9 diagnoses of diabetes, at least 2 prescription fills of diabetes medications of 30 days or greater, or prescriptions for metformin and glucose test strips within a two-year timeframe [10].

Time to Cancer Analysis

The primary outcome was pancreas cancer diagnosis, defined by date of pancreas cancer diagnosis (based on ICD-9 code date, NDI pancreas cancer death date, or Oncology Domain ICD-O-3 date). For patients diagnosed as having pancreas cancer using multiple methods, the earliest date was utilized. Follow-up time began at time of pancreas cyst diagnosis, which was based on first ICD-9 code for pancreas cyst. For patients who did not develop pancreas cancer, they were censored in the following order: (1) date of pancreatic cyst surgery (based on Current Procedural Terminology codes 48,120, 48,140, 48,145, 48,150, 48,152, 48,153, 48,154); (2) date of death during follow-up; (3) date of last VA healthcare encounter; or (4) end of study period (12/31/2014). Date of pancreatic surgery was used to censor because the natural history of the underlying cyst is altered among pancreas cyst patients undergoing surgical resection. If pancreas cancer was present at time of surgery, patients were included in the analysis.

Pancreatic cancer incidence was computed in two ways: (a) using person-years follow-up approach; (b) 5-year, 10-year, and 15-year absolute incidence of pancreatic cancer. Continuous variables were compared using Wilcoxon rank sum test, and categorical variables were compared using Chi-square test between patients with pancreatic cancer and those without pancreatic cancer. Time to pancreatic cancer was estimated using Kaplan–Meier method. Univariable Cox proportional hazards regression was used to study the association between each baseline characteristic and time to pancreatic cancer. These variables were included as predictors in the multivariable Cox regression model if found moderately associated (p < 0.15) with the outcome in univariate analysis. Purposive variable selection was used to select main effects in the final model. All covariates significant at p < 0.10 were kept in the final model [11].

Time to Death Analysis

In the time to death analysis, the outcome of interest was date of death based on NDI. For patients who did not die during the study period, they were censored based on: (1) date of last VA healthcare encounter or (2) end of study period (12/31/2014). The survival probability among patients who develop pancreas cancer was compared to patients who did not develop cancer using the Kaplan–Meier method. Univariable Cox proportional hazards regression was used to study the association between each baseline characteristic and time to death. These variables were included as predictors in the multivariable Cox regression model if found moderately associated (p < 0.15) with the outcome in univariate analysis. Purposive variable selection was used to select main effects in the final model. All covariates significant at p < 0.10 were kept in the final model.

Results

A total of 7211 patients with pancreatic cysts were identified, contributing 31,501 person-years of follow-up (median follow-up 4.4 years). Seventy-nine of 7211 patients (1.1%) developed pancreas cancer. Thirty-nine of the 79 cancer patients had pathologically confirmed adenocarcinoma; 37 of the 79 patients with pancreatic cancer were identified based on NDI data alone; 3 of the 79 patients were included based on very high clinical suspicion for pancreas cancer, in whom biopsy was never pursued. Pancreas cyst patients who developed cancer were older at time of cyst diagnosis compared to patients who never developed cancer (74.1 years vs. 69.2 years; p < 0.01). Other demographic characteristics, including gender, race/ethnicity, diabetes, tobacco use, BMI, and CCIS, were similar between individuals with cysts who developed versus did not develop pancreatic cancer on follow-up (Table 1).

Table 1

Demographic and clinical characteristics of pancreas cyst cohort

	Pancreatic cyst patients with pancreatic cancer	Pancreatic cyst patients without pancreatic cancer	p value
Total N	79	7132
Mean age at cyst diagnosis (SD)	74.1 (10.3)	69.2 (11.8)	< 0.001
Gender, N (%)
Men	78 (98.7%)	6780 (95.1%)	0.19
Women	1 (1.3%)	352 (4.9%)
Race/ethnicity, N (%)
White	51 (64.6%)	5162 (72.4%)
African American	15 (19.0%)	1040 (14.6%)	0.27
Other	13 (16.5%)	930 (13.0%)
Diabetes mellitus, N (%)	25 (31.6%)	2228 (31.2%)	0.90
Tobacco use, N (%)
Current smoker	17 (21.5%)	1595 (22.4%)	0.91
Former smoker	24 (30.4%)	1920 (26.9%)
Never smoker	22 (27.8%)	2168 (30.4%)
Unknown	16 (20.3%)	1449 (20.3%)
Body mass index, N (%)
≤ 18.5 kg/m2	1 (1.3%)	137 (1.9%)	0.60
18.5–25 kg/m2	20 (25.3%)	1781 (25.0%)
25–30 kg/m2	25 (31.6%)	2397 (33.6%)
30–35 kg/m2	21 (26.6%)	1339 (18.8%)
35–40 kg/m2	5 (6.3%)	445 (6.2%)
> 40 kg/m2	0 (0.0%)	159 (2.2%)
Unknown	7 (8.9%)	874 (12.3%)
Charlson comorbidity Index score
0	18 (22.8%)	1715 (24.0%)	0.84
1	16 (20.3%)	1662 (23.3%)
2	14 (17.7%)	1351 (18.9%)
3	10 (12.7%)	893 (12.5%)
≥ 4	21 (26.6%)	1511 (21.2)

This table shows the demographic and clinical characteristics of pancreas cyst patients who developed versus did not develop pancreas cancer

Of the 79 patients who developed pancreas cancer on follow-up, the median time to cancer diagnosis after cyst diagnosis was 2.5 years (interquartile range 2.1 years). Sixty-eight patients (86.1% of cancer patients, 0.9% of pancreas cyst cohort) developed cancer 1–5 years after cyst diagnosis; 10 patients (12.7% of cancer patients, 0.14% of pancreas cyst cohort) developed cancer 5–10 years after cyst diagnosis; 1 patient (1.3% of cancer patients, 0.01% of pancreas cyst cohort) developed cancer more than 10 years after cyst diagnosis. Figure 1 illustrates the cumulative incidence of pancreatic cancer.

Fig. 1

Cumulative incidence of pancreatic cancer among individuals with pancreatic cysts (n = 7211). This figure shows the cumulative incidence of pancreas over time among patients with pancreas cysts utilizing the Kaplan–Meier method. The gray-shaded area represents the 95% confidence interval

Multivariable analysis of predictors of pancreas cancer identified older age (every 1 year; HR 1.05; 95% C.I. 1.03–1.08) and blacks compared to non-Hispanic whites (HR 1.92; 95% C.I. 1.07, 3.45) as the only statistically significant variables associated with development of pancreas cancer. Gender, diabetes, tobacco use, BMI, and CCIS were not associated with risk for incident pancreatic cancer on multivariable analysis.

Of the 7211 patients with pancreas cysts in our cohort, 1982 patients (27.5%) died during study follow-up. Sixty-three patients died from pancreas cancer (3.2% of patient deaths; 0.9% of pancreas cyst cohort). As expected, patients with pancreas cancer had shorter survival compared to patients who never developed pancreas cancer. Patients who died of any cause during follow-up were older at cyst diagnosis (73.7 years vs. 67.6 years; p < 0.01), more frequently men (97.8% vs. 94.1%, p < 0.01), more commonly had diabetes (33.8% vs. 30.3%, p < 0.01), had lower BMI (26.8 vs. 28.4, p < 0.01), and had higher mean CCIS at cyst diagnosis (2.8 vs. 1.9; p < 0.01) compared to patients who were alive during follow-up. Multivariable analysis identified older age (every 1 year; HR 1.05; 95% CI 1.04–1.05), male gender (HR 1.37, 95% CI 1.00–1.89), current smoking (HR 1.41, 95% CI 1.22–1.63, compared to never smoker), lower BMI (HR 0.96, 95% CI 0.95–0.97), and higher CCIS at time of cyst diagnosis as significant predictors of death among patients with pancreas cysts (Table 2).

Table 2

Multivariable model of factors associated with all-cause mortality among patients with pancreatic cysts

Covariate	HR (95% CI)	p value
Age (every 1 year)	1.05 (1.04, 1.05)	< 0.001
Men (Ref = Women)	1.37 (1.00, 1.89)	0.053
Tobacco use (Ref = Never Smoker)
Current smoker	1.41 (1.22, 1.63)	< 0.001
Former smoker	0.98 (0.86, 1.11)	0.75
Unknown	1.06 (0.92, 1.21)	0.43
Body mass index (every 1 kg/m2)	0.96 (0.95, 0.97)	< 0.001
Charlson comorbidity index score (Ref = 0)
CCIS 1	1.39 (1.17, 1.65)	< 0.001
CCIS 2	1.74 (1.46, 2.06)	< 0.001
CCIS 3	2.35 (1.96, 2.81)	< 0.001
CCIS ≥ 4	3.19 (2.73, 3.74)	< 0.001

Other variables including ethnicity and diabetes were not statistically significant in multivariable analysis

The most common cause of death among all individuals with pancreatic cysts (Fig. 2) was non-pancreas cancer (498 patients; 25.1% of deaths) followed by cardiovascular disease (398 patients; 20.1% of deaths). Among the 561 patients with death attributable to cancer, lung cancer was the most common cancer cause of death (137 patients) followed by pancreas (63 patients), prostate (51 patients), and colorectal cancer (44 patients).

Fig. 2

Causes of death among patients with pancreas cysts; N = 1982 This figure shows the most common causes of death among patients with pancreas cysts

Discussion

In this large national cohort of over 7000 US Veterans with 31,501 person-years of follow-up, we found incidence of pancreatic cancer was low, at just 1.1%, corresponding to 2.5 cancers per 1000 person-years follow-up. Further, among patients who develop pancreas cancer, the majority were diagnosed within the first five years of cyst diagnosis, with 14% developing cancer more than 5 years after pancreas cyst diagnosis. Importantly, among all patients with pancreatic cysts undergoing long-term follow-up, we found pancreas cancer-associated mortality was very low, at just 0.9%, with pancreatic cancer accounting for only 3.2% of all patient deaths occurring during follow-up. The vast majority of patients who died during follow-up died from non-pancreas cancers (25%) and cardiovascular causes (20%).

Our observation of low pancreatic cancer incidence among patients with pancreatic cysts confirms and extends findings from prior population-based studies [8, 12–14]. Scheiman et al. estimated the lifetime risk of developing pancreas cancer from a pancreas cyst as being less than 1% when utilizing national US population data and Surveillance, Epidemiology, and End Results (SEER) cancer statistics [8]. Wu et al. evaluated the incidence of pancreas cancer among a surveillance cohort of 1735 patient with pancreas cysts, of whom only 14 (0.8%; 4.1 cancers per 1000 person-years follow-up) developed pancreas cancer with average follow-up of 23.4 months (3384 person-years) [13]. The low incidence observed in our study and other population-based studies contrasts with higher pancreatic cancer incidence reported by non-population-based studies. Surgical series usually describe a higher rate of cancers, with a pooled proportion of cysts with invasive cancer closer to 15% [8]. Applicability of surgical series to clinical decision making is challenged by the limitations created by (a) selection bias, (b) unclear understanding about the long-term outcomes in patients with non-malignant cysts at time of surgery (i.e., would some of these cysts eventually develop cancer?), and (c) presence of substantial proportion of patients with symptomatic pancreas cysts. Other studies (single-center studies, endosonographic databases) also describe higher incidence rates of pancreas cancer (4–8%) [15, 16]. The higher pancreas cancer incidence rates in these studies may be related to referral bias as most of these studies are often from tertiary referral centers and may consist of patients with higher-risk cystic lesions. Other limitations of these studies are small size, small numbers of patients developing pancreatic cancer, and short-term follow-up resulting in imprecise estimates of pancreatic cancer incidence. In contrast to surgical series and other non-population-based studies, our study expands the number of population-based studies confirming that incidence of pancreatic cancer among individuals with pancreatic cysts is low.

Our study also provides new knowledge regarding timing of cancer development following pancreas cyst diagnosis. The American Gastroenterological Association guidelines from 2015 suggest stopping surveillance of pancreas cysts after 5 years if no change in cyst characteristics—a recommendation which was met with criticism [17]. The present study does not specifically address whether the 14% of patients who developed cancer after 5 years had change in cyst characteristics (such as cyst growth and development of mural nodule) but does demonstrate continued surveillance may be important for a substantial portion of patients with pancreas cysts.

Our study contributes new insights on the mortality after pancreatic cyst diagnosis, an area that has not been widely studied. Kwok et al. examined pancreas cancer-associated mortality among a cohort of patients with pancreas cysts [18]. Of 1800 pancreatic cyst patients, 402 patients (22.3%) died during a median 5.7-year follow-up. Pancreatic cancer-associated mortality was 2.4% (n = 43), and 11% of all deaths were due to pancreas cancer. Kwok et al. found a higher pancreas cancer-associated mortality compared to our study (2.4% vs. 0.9%), but this may be explained by differences in study design. In the Kwok study, pancreatic cancer was considered preexisting only if diagnosed within 3 months of cyst diagnosis, while our study excluded anyone with pancreas cancer within one year of pancreas cyst diagnosis due to challenges in cancer diagnosis among patient with cystic neoplasms. Median time to cancer diagnosis in the Kwok study for the 43 pancreas cancer deaths was 0.6 years (IQR 0.2, 2.1), raising concern some patients had cancer at time of cyst diagnosis. The inclusion of some prevalent cancers likely resulted in higher pancreas cancer-associated mortality compared to our study. Nonetheless, both Kwok et al. and our study support that mortality attributable to pancreatic cancer among patients with pancreatic cysts is low, and much lower that rates of deaths from other causes.

Several limitations may be considered in interpreting our findings. First, this is a retrospective study with inherent limitations. The ideal study design for evaluating the natural history of pancreas cysts would be a prospective cohort study. While such a study would be the gold standard, there are multiple issues with feasibility, cost, and length of follow-up. The outcome of interest, pancreatic cancer, is a rare event, and the size of the study and length of follow-up required would be very costly and inefficient. Second, the current study does not describe cyst-specific characteristics for the cohort of pancreas cyst patients. Due to the size of the cohort and need for manual chart review to extract cyst-specific characteristics (such as size, presence of nodule/solid component, pancreas duct dilation, and cyst growth) obtaining cyst-specific characteristics for the entire cohort is not feasible. Third, another limitation is the use of ICD-9 coding for creation of the pancreas cyst cohort. However, a review of over 400 charts revealed 79.1% of our cohort met all inclusion criteria. Four, the study is performed among US Veterans and may not be generalizable to other settings. Over 95% of the pancreas cyst cohort consists of male patients, and certain cystic neoplasms such as MCNs which occur almost exclusively in women are likely under-represented in this cohort. IPMNs, which are felt to be the most common pancreas cyst, likely account for the vast majority of cystic lesions in this study. Our current study has the potential to both underestimate and overestimate pancreas cancer incidence. Pancreas cancer incidence may be underestimated because approximately 20% of the pancreas cyst cohort did not meet inclusion criteria. However, even if the incidence rate we observed was doubled, the overall estimated incidence of pancreatic cancer would be low. Pancreas cancer incidence may also be overestimated in this study because ICD-9 codes were used to identify patients with pancreas cysts, but it is possible that patients with small cysts were not consistently coded, potentially leading to under ascertainment of true cyst diagnosis.

There are several strengths to this study. First, it represents the largest known cohort of pancreas cyst patients with long-term follow-up. Furthermore, this study allows the evaluation of natural history of pancreas cysts with long-term follow-up while minimizing referral bias. Second, compared to other administrative database studies, the use of VINCI allows manual chart review and verification of important diagnoses such as pancreas cancer. Third, there is a robust ascertainment strategy for pancreas cancer diagnosis with use of ICD-9 codes, NDI cause of death data, and cancer registry data. One advantage of using NDI cause of death data is obtaining information regarding patients who may have obtained care from outside of the VA. Furthermore, due to the poor prognosis of pancreas cancer (5-year survival < 9%), the vast majority of patients who develop pancreas cancer die within a short period of time and will likely be included within NDI cause of death data. This robust ascertainment strategy ensures high accuracy in identifying pancreas cancer cases.

In conclusion, the overall incidence of pancreas cancer among Veterans with pancreas cysts is low (1.1%). While most patients develop cancer within the first 5 years of cyst diagnosis, 14% of patients who developed cancer did so after five years. Pancreas cancer-associated mortality is very low (0.9%), and patients are six times more likely to die from a non-pancreas cancer than pancreas cancer. Given incidence of pancreas cancer and death from pancreas cancer are rare, future research should focus on identifying criteria for selecting individuals at high risk for developing and dying from pancreatic cancer.