- Research article
- Open Access
- Open Peer Review
Comparative effectiveness of primary tumor resection in patients with stage III pancreatic adenocarcinoma
BMC Cancervolume 19, Article number: 761 (2019)
Previous studies comparing primary tumor resection (PTR) to palliative treatment for advanced-stage pancreatic ductal adenocarcinoma (PDA) were limited by strong selection bias. We used multiple methods to control for confounding and selection bias to estimate the effect of PTR on survival for late-stage PDA.
Surveillance, Epidemiology, and End Results (SEER) 18 registry database for 2004 through 2014 was retrieved for the present study. A total of 4322 patients with stage III (AJCC, 6th) PDA were included in this study. Propensity score matching (PSM) was performed to eliminate possible bias. In addition, instrumental variable (IV) analysis was utilized to adjust for both measured and unmeasured confounders.
A total of 4322 patients with stage III PDA including 552 (12.8%) who underwent PTR, 3770 (87.2%) without PTR, were identified. In the multivariable cohort, a clear prognostic advantage of PTR was observed in overall survival (OS) (P < 0.001) and disease-specific survival (DSS) (P < 0.001) compared to patients after non-surgery therapy. In the PSM cohort, patients in PTR group showed a better OS and DSS (both P values < 0.001) compared to patients in non-surgery group. The survival benefit of PTR for stage III PDA was not observed in the two-stage residual inclusion (2SRI) model. Estimates based on this instrument indicated that patients treated with PTR had similar OS (P = 0.448) and DSS (P = 0.719). In IV analyses stratified by chemotherapy and tumor location, patients undergoing PTR had similar OS and DSS compared to patients in non-surgery group across all subgroups.
Survival with PTR did not differ significantly from palliative treatment in marginal patients with stage III pancreatic adenocarcinoma. High-quality randomized trials are needed to validate these results.
The incidence of pancreatic ductal adenocarcinoma (PDA) continues to increase. In 2017, an estimated 53670 new cases (female: 25700; male: 27970) of pancreatic cancer were diagnosed within the US and 43090 individuals (female: 20790; male: 22300) were expected to die of the tumor . Primary tumor resection (PTR) is the only curative modality, while more than 80% of tumors were unresectable when present . The 5-year survival of patients with PDA after surgery is approximately 20% (the median survival is 15–23 months) [3,4,5,6,7,8].
According to American Joint Committee on Cancer (AJCC) classification version 6, patients with stage III PDA (tumors involved celiac axis and/or superior mesenteric artery) can be divided into borderline resectable and unresectable, depending on the extent of the tumor encasement of major vessels [2, 5, 6]. Previously, chemoradiotherapy has been carried out to reduce the risk of a positive surgical margin and distant metastasis [2, 4, 9]. Since the FOLFIRINOX regimen (irinotecan, oxaliplatin, leucovorin, and fluorouracil) was introduced in 2011 by a prospective randomized controlled trial , it has been reported to result in objective response rates that were 2–3 fold higher than other regimens in PDA . Several studies have confirmed that a large number of cases even with locally advanced and unresectable PDA can be converted to be resectable by FOLFIRINOX [12,13,14].
Previous publications have reported a discrepant overall survival of PDA patients with advanced disease undergoing PTR and vascular reconstruction (the median survival ranged from 12 to 35 months) [15,16,17,18]. Owing to the varying outcomes regarding to the long-term survival of underline resectable PDA receiving PTR, we designed a population-based cohort study to explore the independent role of PTR in patients with stage III PDA (S-III PDA). We utilized an instrumental variable (IV) analysis to determine variation in outcomes across geographical areas that were different in PTR rates. The IV analysis is aimed to control for potential unknown confounding factors in decision making for surgeries [19, 20]. In the present study, PTR rates in various Health Service Areas (HSA) was employed as our instrument. The treatment option (PTR or non-PTR) for marginal patients (those with a borderline or uncertain need for PTR) may be affected by preferences, beliefs, or surgical skills of surgeons in their HSAs. Patients with S-III PDA would be performed PTR in a high-use HSA, while not in a low-use HSA [21, 22]. The coefficient in the IV analysis represents the adjusted treatment effect for the marginal population rather than the average treatment effect .
Surveillance, Epidemiology, and End Results (SEER) 18 database for 2004 through 2014 was retrieved for this study (seer.cancer.gov/about/overview.html). The SEER population-based cancer registries covers approximately 28% of the US population, which collects data of tumor incidence, demographics, tumor characteristics and patient survival. Firstly, 107544 patients with PDA was identified based on the pathological diagnosis. The ICD-O-3 (International Classification of Diseases for Oncology, 3rd Edition) site code is C25 and histologic type codes are 8140, 8500, 8010, 8000, 8480, 8481, 8490, 8255, 8021, 8020, 8521, 8141, 8022, 8144. Tumor, node, and metastasis stage of PDA in SEER was based on AJCC stage version 6. The flow diagram of patient selection is shown in Fig. 1. Finally, a total of 4322 cases with stage III PDA were included in the analysis. The following codes related to PDA treatment were selected: PTR: 30 (partial pancreatectomy), 35–37 (Whipple), 40 (whole pancreatectomy), 60 (whole pancreatectomy with subtotal gastrectomy/duodenectomy), 70 (an extended pancreatoduodenectomy) and 80 (pancreatectomy, NOS); none surgical treatment: 0. This study has been approved by the Institutional Review Board at the West China Hospital. All patient data from SEER database is public available and anonymous.
Continuous data are shown as mean ± SD and categorical variables are presented as number (%). The continuous variable was examined by t-test or Kruskal-Wallis H test the categorical data was tested by Chi-square test or Fisher’s exact test. Overall survival (OS) and disease-specific survival (DSS) were the primary endpoints. The former was defined as the time from the date of treatment to the date of death with any cause and the latter was defined as interval until death caused by PDA. The Kaplan-Meier method was used to analyze survival data (compared by the log-rank test). Multivariable analyses were performed by the Cox proportional hazards regressions. Based on previous study , the interaction tests were also carried out to identify the interactive factors influencing the relationship between treatment methods and survival.
To further adjust for potential baseline confounding factors, the propensity score matching (PSM) was carried out according to the following parameters: sex, age, race, year of diagnosis, tumor differentiation, tumor size, and chemotherapy. Cases with PTR were matched to those without PTR with a matching ratio of 1:3. The nearest-neighbor PSM was performed by logistic regression.
PTR rates in HSAs were utilized as an IV. In this study, the PTR rate in HSAs is a qualified instrument because it is highly associated with the probability of a patient’s exposure to primary tumor resection (F statistic > 10) and is also not related to patient survival. In addition, covariate balance across quintiles was also examined. We did not include patients living in HSAs with fewer than 20 patients, given the difficulty to confirm the PTR rates in those HSAs . The two-stage residual inclusion (2SRI) method was utilized for instrumental variable analyses .
A total of 4322 S-III PDA patients including 552 (12.8%) who underwent PTR, 3770 (87.2%) without PTR, were identified. Table 1 displayed the general demographics of the final cohort of 4322 S-III PDA patients with available variables. The mean age of patients undergoing PTR and none was 64.6 and 66.7 years, respectively. Compared with the non-surgery group, the PTR group had more pancreatic head tumor (79.3% vs. 70.9%) and smaller tumor size (38.5 ± 16.3 mm vs. 41.0 ± 15.4 mm). The other clinicopathologic characteristics including sex, race and number of patients receiving chemotherapy showed no statistically significance between PTR and non-surgery groups (all P > 0.05).
In multivariable analyses, we included a total of 4322 patients with known prognostic data. The mean overall survival time (in the total cohort) for patients who underwent PTR and patients undergoing non-surgery were 23.9 months and 14.5 months, respectively. The mean DSS time for cases after PTR and none were 23.9 and 15.4 months, respectively. Patients with PTR had longer OS (P < 0.001) and DSS (P < 0.001) compared to patients with non-surgery treatment (Fig. 2a and c).
In the cohort for multivariable analyses (OS: n = 4322; DSS: n = 4055), after adjusting for potential confounding factors, a clear prognostic advantage of PTR was observed in OS (HR, 0.59; 95% CI, 0.53 to 0.66; P < 0.001) and DSS (HR, 0.57; 95% CI, 0.51 to 0.65; P < 0.001) compared to patients after non-surgery therapy (Table 2).
Instrumental variable analysis
To be valid, an instrumental variable must meet two conditions: 1) the variable must be highly associated with the treatment of interest (in this study receipt of PTR); and 2) the instrumental variable cannot be related to the outcomes (in this study survival) except through its effect on the treatment received .
The average PTR rate in HSAs fluctuated from a low of 16% (quintile 1) to a high of 27% (quintile 5). The F-statistics was 3325.8 (P < 0.001), indicating that the instrument was strongly related to the treatment. In addition, in a standard COX regression, no significant correlation was observed between the IV and OS (HR 1.18, 95% CI 0.27–5.05, P = 0.827). We divided patients into quintiles according to the proportion of cases within each HSA who underwent PTR (Additional file 1: Table S1). Most of the clinicopathologic features were balanced across quintiles. Consequently, these observations suggest that HSA PTR rate meets the two requirements for a valid instrument.
For patients with S-III PDA, the salutary benefit of PTR in survival was not observed in the 2SRI model. In IV analysis, results indicated that patients who underwent PTR had similar OS (HR 0.74, 95% CI 0.34–1.61, P = 0.448) and DSS (HR 0.86, 95% CI 0.38–1.94, P = 0.719) after adjusting confounding factors (Table 3).
In IV analyses stratified by chemotherapy, we found that the similar effects of PTR vs. none on patient survival were consistent across both subgroups (Table 4). In IV analyses, patients in the PTR group receiving chemotherapy had similar OS (HR 0.43, 95% CI 0.09–2.15, P = 0.304) and DSS (HR 0.56, 95% CI 0.10–3.14, P = 0.508) compared to patients in the non-PTR group receiving chemotherapy. In IV analyses stratified by tumor location, we found that the similar effects of PTR vs. non-PTR on survival (both OS and DSS) were consistent across all subgroups with different tumor location (Table 5).
In IV analyses stratified by the other clinicopathologic characteristics (sex, race, age, year of diagnosis and tumor size.), we found that the treatment effect of PTR (OS and DSS) was consistent across all the subgroups (data not shown).
Propensity score matched analyses
In the propensity-matched population, all the potential prognosis-relevant characteristics were well-balanced for most of the baseline features (Table 1). In the PSM cohort, results from the univariate analysis indicated that cases with PTR had better OS (HR 0.59, 95% CI 0.50–0.69 P < 0.001) and DSS (HR 0.57, 95% CI 0.47–0.67, P < 0.001) compared to patients with non-surgery treatment (Table 2). In the PSM-adjusted population, patients in PTR group still showed a better OS and DSS (both P values < 0.001) compared to patients in non-surgery group by Kaplan-Meier method (Fig. 2b and d).
The HRs (PTR vs. none) adjusted by propensity score showed both longer OS (continuous: HR 0.61, 95% CI 0.55–0.67, P < 0.001; quintile: HR 0.62, 95% CI 0.56–0.69, P < 0.001) and DSS (continuous: HR 0.61, 95% CI 0.54–0.68, P < 0.001; quintile: HR 0.60, 95% CI 0.53–0.67, P < 0.001) associated with PTR (Table 2).
Pancreatic resection is associated with better outcomes for early-stage PDA.  However, for patients with underline resectable PDA, though chemotherapy regimen such as FOLFIRINOX increased the tumor resectability, the long-term survival (OS and DSS) in this study was comparable to those receiving non-surgery treatment in IV analyses. This conclusion is inconsistent with previous studies demonstrating that S-III PDA patients had a better survival after PTR compared to those without surgery [15, 18, 28,29,30].
Among previous studies related to PTR versus non-surgical management in cases with advanced PDA, there were two randomized controlled trials (RCTs) [28, 29] comparing PTR versus non-surgical treatments in cases with underline resectable PDA. Both RCTs included patients with locally advanced PDA invading the serosa anteriorly or retroperitoneum posteriorly or involving the major vascular structures. One RCT enrolled cases with PDA in different location of the pancreas and another included only cases with tumor in the pancreatic head or neck. Both studies demonstrated that patients receiving PTR and vascular resection and reconstruction had longer survival compared to patients only undergoing non-surgery treatment such as chemoradiotherapy. However, both researches were at high risk of bias and only a small number of patients (47 and 51 patients in two studies, respectively) were included.
In this study, utilizing IV analyses, we concluded that patients with S-III PDA receiving PTR had a similar long-term prognosis compared to those without PTR (only receive non-surgery treatments). We have applied both traditional regression analyses and propensity score methods to explore relations between surgical methods and long-term survival. However, these analytic methods cannot adjust for unknown confounding factors [15, 26]. In contrast, results from IV analyses (2SRI model) were observed to be more close to outcomes from RCTs . Given the lack of high-quality RCTs associated with PTR vs. none for PDA patients, results of IV analyses may represent the best evidence available to guide treatment decision-making. It should be noted that IV analyses estimate the treatment effect on the marginal population rather than the average treatment effect of PTR [19, 21]. The marginal population represents the population that would receive PTR in a high-use HSA but not in a low-use HSA. The IV analysis does not rely on defining the specific clinical parameters of these populations. Instead, it is based on the hypothesis that patients reside randomly around hospitals and some patients are treated differently in different centers.
There are several limitations to this study. First, although we have acquired the data related to the chemotherapy from SEER database, the detailed regimens and the timing of chemotherapy were yet inaccessible. Patients without chemotherapy or with unknown data of chemotherapy accounted for nearly 80% of patients with pancreatic cancer. In addition, we cannot divide S-III PDA into borderline resectable and unresectable groups based on the extension of tumor invasion, thus we cannot assess whether patients with surgical resection was well chosen. The proportion of patients with R0 tumor resection could not be confirmed in this study. Second, patient performance status and presence of comorbidities are risk factors for patient prognosis. However, the SEER database does not provide these data, thus we cannot adjust these factors by multivariable analyses. Third, details on postoperative morbidities were extremely limited, thus we could not evaluate the influence of treatment methods to the short-term prognosis. Fourth, the observations of this study should be interpreted cautiously, due to some patients with unknown covariates in the SEER database were excluded from the analyses. Finally, even though treatment rate (PTR rate) is a useful practical IV, there remains potential for instrument-outcome confounding, such as receipt of other treatments also associated with our instrument and the outcome. In addition, IV analysis only estimates the effect on marginal population, while the marginal population excludes patients who would always or never receive PTR, focusing on PDA patients whose indications for PTR are more uncertain (Fig. 3). [20, 26, 31]
In conclusion, by integrating results from multivariate COX regression, PSM models and IV analysis, our study demonstrated that PTR provided similar overall and tumor-specific survival benefits in cases with stage III PDA compared to patients with palliative treatments. Further high-quality prospective randomized trials are needed to validate this conclusion and further investigations are required to identify late-stage patients suitable for PTR.
Availability of data and materials
All primary data is available by sending email to: firstname.lastname@example.org or downloading from SEER database.
American Joint Committee on Cancer
National Comprehensive Cancer Network
Propensity score matching
Primary tumor resection
Surveillance, Epidemiology, and End Results
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2017. CA Cancer J Clin. 2017;67(1):7–30.
Tempero MA, Malafa MP, Al-Hawary M, Asbun H, Bain A, Behrman SW, Benson AB 3rd, Binder E, Cardin DB, Cha C, et al. Pancreatic adenocarcinoma, version 2.2017, NCCN clinical practice guidelines in oncology. J Natl Compr Canc Netw. 2017;15(8):1028–61.
Yeo CJ, Abrams RA, Grochow LB, Sohn TA, Ord SE, Hruban RH, Zahurak ML, Dooley WC, Coleman J, Sauter PK, et al. Pancreaticoduodenectomy for pancreatic adenocarcinoma: postoperative adjuvant chemoradiation improves survival. A prospective, single-institution experience. Ann Surg. 1997;225(5):621–33 discussion 633-626.
Neoptolemos JP, Stocken DD, Bassi C, Ghaneh P, Cunningham D, Goldstein D, Padbury R, Moore MJ, Gallinger S, Mariette C, et al. Adjuvant chemotherapy with fluorouracil plus folinic acid vs gemcitabine following pancreatic cancer resection: a randomized controlled trial. Jama. 2010;304(10):1073–81.
Sabater L, Munoz E, Rosello S, Dorcaratto D, Garces-Albir M, Huerta M, Roda D, Gomez-Mateo MC, Ferrandez-Izquierdo A, Darder A, et al. Borderline resectable pancreatic cancer. Challenges and controversies. Cancer Treat Rev. 2018;68:124–35.
Lopez NE, Prendergast C, Lowy AM. Borderline resectable pancreatic cancer: definitions and management. World J Gastroenterol. 2014;20(31):10740–51.
Varadhachary GR, Tamm EP, Abbruzzese JL, Xiong HQ, Crane CH, Wang H, Lee JE, Pisters PW, Evans DB, Wolff RA. Borderline resectable pancreatic cancer: definitions, management, and role of preoperative therapy. Ann Surg Oncol. 2006;13(8):1035–46.
Springett GM, Hoffe SE. Borderline resectable pancreatic cancer: on the edge of survival. Cancer control : journal of the Moffitt Cancer Center. 2008;15(4):295–307.
Chin V, Nagrial A, Sjoquist K, O'Connor CA, Chantrill L, Biankin AV, Scholten RJ, Yip D. Chemotherapy and radiotherapy for advanced pancreatic cancer. Cochrane Database Syst Rev. 2018;3:Cd011044.
Vaccaro V, Sperduti I, Milella M. FOLFIRINOX versus gemcitabine for metastatic pancreatic cancer. N Engl J Med. 2011;365(8):768–9 author reply 769.
Sadot E, Doussot A, O'Reilly EM, Lowery MA, Goodman KA, Do RK, Tang LH, Gonen M, D'Angelica MI, DeMatteo RP, et al. FOLFIRINOX induction therapy for stage 3 pancreatic adenocarcinoma. Ann Surg Oncol. 2015;22(11):3512–21.
Luu AM, Herzog T, Hoehn P, Reinacher-Schick A, Munding J, Uhl W, Braumann C. FOLFIRINOX treatment leading to pathologic complete response of a locally advanced pancreatic cancer. J Gastrointest Oncol. 2018;9(2):E9–e12.
Nitsche U, Wenzel P, Siveke JT, Braren R, Holzapfel K, Schlitter AM, Stoss C, Kong B, Esposito I, Erkan M, et al. Resectability after first-line FOLFIRINOX in initially Unresectable locally advanced pancreatic Cancer: a single-center experience. Ann Surg Oncol. 2015;22(Suppl 3):S1212–20.
Ferrone CR, Marchegiani G, Hong TS, Ryan DP, Deshpande V, McDonnell EI, Sabbatino F, Santos DD, Allen JN, Blaszkowsky LS, et al. Radiological and surgical implications of neoadjuvant treatment with FOLFIRINOX for locally advanced and borderline resectable pancreatic cancer. Ann Surg. 2015;261(1):12–7.
Wang L, Cheng CS, Chen L, Chen Z. Benefit from the inclusion of surgery in the treatment of patients with stage III pancreatic cancer: a propensity-adjusted, population-based SEER analysis. Cancer Manag Res. 2018;10:1907–18.
Beane JD, House MG, Pitt SC, Kilbane EM, Hall BL, Parmar AD, Riall TS, Pitt HA. Distal pancreatectomy with celiac axis resection: what are the added risks? HPB (Oxford). 2015;17(9):777–84.
Denecke T, Andreou A, Podrabsky P, Grieser C, Warnick P, Bahra M, Klein F, Hamm B, Neuhaus P, Glanemann M. Distal pancreatectomy with en bloc resection of the celiac trunk for extended pancreatic tumor disease: an interdisciplinary approach. Cardiovasc Intervent Radiol. 2011;34(5):1058–64.
Klompmaker S, de Rooij T, Korteweg JJ, van Dieren S, van Lienden KP, van Gulik TM, Busch OR, Besselink MG. Systematic review of outcomes after distal pancreatectomy with coeliac axis resection for locally advanced pancreatic cancer. Br J Surg. 2016;103(8):941–9.
McDowell BD, Chapman CG, Smith BJ, Button AM, Chrischilles EA, Mezhir JJ. Pancreatectomy predicts improved survival for pancreatic adenocarcinoma: results of an instrumental variable analysis. Ann Surg. 2015;261(4):740–5.
Baiocchi M, Cheng J, Small DS. Instrumental variable methods for causal inference. Stat Med. 2014;33(13):2297–340.
Valley TS, Sjoding MW, Ryan AM, Iwashyna TJ, Cooke CR. Association of Intensive Care Unit Admission with Mortality among Older Patients with Pneumonia. Jama. 2015;314(12):1272–9.
Tan HJ, Norton EC, Ye Z, Hafez KS, Gore JL, Miller DC. Long-term survival following partial vs radical nephrectomy among older patients with early-stage kidney cancer. Jama. 2012;307(15):1629–35.
Zhao J, Mao J, Li W. Association of Tumor Grade with Long-Term Survival in patients with hepatocellular carcinoma after liver transplantation. Transplant Proc. 2019;51(3):813–9.
Xu H, Xia Z, Jia X, Chen K, Li D, Dai Y, Tao M, Mao Y. Primary tumor resection is associated with improved survival in stage IV colorectal Cancer: an instrumental variable analysis. Sci Rep. 2015;5:16516.
Gore JL, Litwin MS, Lai J, Yano EM, Madison R, Setodji C, Adams JL, Saigal CS. Use of radical cystectomy for patients with invasive bladder cancer. J Natl Cancer Inst. 2010;102(11):802–11.
Terza JV, Basu A, Rathouz PJ. Two-stage residual inclusion estimation: addressing endogeneity in health econometric modeling. J Health Econ. 2008;27(3):531–43.
Mohammed S, Van Buren G 2nd, Fisher WE. Pancreatic cancer: advances in treatment. World J Gastroenterol. 2014;20(28):9354–60.
Doi R, Imamura M, Hosotani R, Imaizumi T, Hatori T, Takasaki K, Funakoshi A, Wakasugi H, Asano T, Hishinuma S, et al. Surgery versus radiochemotherapy for resectable locally invasive pancreatic cancer: final results of a randomized multi-institutional trial. Surg Today. 2008;38(11):1021–8.
Lygidakis NJ, Singh G, Bardaxoglou E, Dedemadi G, Sgourakis G, Nestoridis J, Malliotakis A, Pedonomou M, Solomou EK, Safioleas M, et al. Mono-bloc total spleno-pancreaticoduodenectomy for pancreatic head carcinoma with portal-mesenteric venous invasion. A prospective randomized study. Hepato-gastroenterology. 2004;51(56):427–33.
Scemama U, Birnbaum DJ, Ouaissi M, Turrini O, Moutardier V, Soussan J. Portal vein stent placement in five patients with chronic portal vein thrombosis prior to pancreatic surgery. J Vasc Interv Radiol. 2016;27(6):889–94.
Wan F, Small D. Bias in estimating the causal hazard ratio when using two-stage instrumental variable methods. Stat Med. 2015;34(14):2235–65.
We would like to thank the following: Dr. Chen Chi for providing guidance on data analyses.
Ethics approval and consent to participate
This study is in accordance with the Declaration of Helsinki and has been approved by the Institutional Review Board at the West China Hospital. The data was retrieved after our application was approved by the SEER database.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Characteristics of patients by quintile of Health Services Area PTR rates. (DOCX 18 kb)