A Quantitative Clinicopathological Signature for Predicting Recurrence Risk of Pancreatic Ductal Adenocarcinoma After Radical Resection

Recurrence and distant metastases were main reasons of unfavorable outcomes for patients with pancreatic ductal adenocarcinoma (PDAC) after surgery. The aim of this study was to describe the patterns, timing, and predictors of recurrence or metastasis in PDAC patients after curative surgery. Patients with PDAC who underwent radical pancreatectomy were included. Associations between clinicopathological and radiological characteristics and specific pattern of progression were investigated. Least absolute shrinkage and selection operator (LASSO) and Cox regression were applied to assess the prognostic factors for overall survival (OS) and progression-free survival (PFS). A total of 302 patients were included into present study, and 173 patients were documented as recurrence after a median survival of 24.7 months. More than half of patients recurred after 12 months after surgery, and the liver was the most common metastatic site. Decreased time interval to progression, elevated carbohydrate antigen 19-9 (CA19-9) level, and lymph node (LN)16 metastasis were independent predictors for reduced OS. Independent prognostic factors for PFS included elevated carcinoembryonic antigen (CEA) level, local progression, liver or lung-only metastasis, local + distant progression, multiple metastases, LN16 metastasis, imaging tumor size, chemotherapy, and tumor–node–metastasis (TNM) stage. The predictive system showed valuable prediction performance with values of concordance indexes (C-indexes) and the area under the receiver operating characteristic curve (AUC) over 0.80. Different survival curves and predictive factors for specific patterns of disease progression suggested the biological heterogeneity, providing new versions into personal management of recurrence in PDAC patients after surgery.

Introduction

Pancreatic ductal adenocarcinoma (PDAC) is a lethal disease and is predicted to become the second leading cause of cancer-specific death by 2030 (1). Surgery followed by adjuvant chemotherapy has been widely established as the best mean to obtain longer survival. However, this combination therapy can only be applied to 20% of patients, whereas most patients suffered from locally advanced or metastatic diseases, owing to the lack of early clinical symptoms and effective screening methods. Moreover, even after curative resection, up to 80% of patients suffered from recurrence soon after surgery (2–4), and the 5-year survival rate was <6% (5).

Progression had a truly negative effect on prognoses of patients with PDAC. However, the variations of biological behaviors and clinicopathological factors of tumors would contribute to different patterns and timing of progression even when diseases were classified as the same stages. Although multiple studies illustrated the risk factors of progression, such as resection margin status and lymph node (LN) metastasis (6, 7), the relationship between the prognosis and progression was rarely evaluated for patients with PDAC. The prognosis might be changing among patients with different patterns and timing of progression, whereas significant heterogeneity existed among the current reports regarding patterns and timing of recurrence owing to the small sample sizes and limited period of follow-up (8, 9). Understanding both the risk factors and the patterns of progression of PDAC patients can provide an insight into optimization of the treatment, as well as the surveillance strategies. Although recurrence was associated with decreased survival, whether the sites and timing of recurrence had different influences on survival remained controversial. Thus, the aim of this study was to evaluate the risk factors for different patterns of recurrences and compare the survival differences in PDAC patients with varied patterns or timing of disease progression.

Materials and Methods

Patients

This study included consecutive patients with PDAC who underwent surgical resection at Sun Yat-sen University Cancer Center (SYSUCC) between 2008 and 2018. Excluded patients were those with metastatic diseases detected at diagnosis by radiological examination, such as computed tomography (CT) and magnetic resonance imaging (MRI). Positron emission tomography/CT (PET/CT) and diagnostic laparoscopy were also selectively performed to detect metastases on the basis of the recommendation of the pancreatic multidisciplinary team. The resection margin for radical margin was defined as 1.5–2 mm, which was the same as that of previous studies (10, 11). Excluded were also patients with microscopic or macroscopic incomplete resection, a history of secondary tumors, period of follow-up <1 year, and missing follow-up records.

Data Collection

Resectability was judged on the basis of CT or MRI, and staging was determined by the pathological factors in accordance with the 8th edition of American Joint Committee on Cancer staging system (12, 13). A team specialized in pancreatic surgery performed all radical pancreatic resection. An experienced pancreatic pathologist assessed all the surgical specimens, made the diagnosis of PDAC, and described the pathological variables, including tumor site, tumor size, tumor differentiation, T-stage, LN status (N-stage), LN total number, positive LN number, macrovascular invasion, microvascular invasion, lymph vessel invasion, perineural invasion, adjacent organ invasion, and satellite foci. LN ratio (LNR) was defined as the number of LNs with metastases divided by the total number of excised LNs. Several radiological variables, including imaging tumor size, LN metastasis, vascular invasion, and LN size, were analyzed. Inflammation-based indexes, such as neutrophil-to-lymphocyte ratio (NLR), platelet-to-lymphocyte ratio (PLR), modified Glasgow Prognostic Score (mGPS), prognostic nutritional index (PNI), prognostic index (PI), and systemic immune-inflammation index (SII), were included in this study and calculated according to previous studies (14, 15). Clinical data were also analyzed in this study, including age, gender, white blood cell (WBC) count, C-reactive protein (CRP), albumin (ALB), serum levels of carbohydrate antigen 19-9 (CA19-9), and carcinoembryonic antigen (CEA).

Follow-Up and Recurrence

The follow-up of patients occurred at the outpatient clinic of our hospital. In general, follow-up strategies consisted of regular chest CT, abdominal CT, and CA19-9 test at least every 2 months during the first year after surgical resection and every 3 months thereafter. Occasional additional imaging modalities, such as MRI and PET/CT, were selectively performed to determine patterns of recurrence. Follow-up data were retrieved at the end of May 2019. The categories of regression patterns in the study conducted by Groot et al. (4) were adopted in this study. When considering patterns of recurrence, only the first location of recurrence was documented, and local recurrence and distant recurrence were registered, separately. In addition, distant recurrences were judged as “liver-only” and “lung-only” recurrences for the isolated hepatic and pulmonic recurrence, respectively, and “others” for isolated recurrence in other less common locations. If both local recurrence and isolated distant metastasis occurred or multiple distant metastases were detected at the same time, recurrences were defined as “local + distant” and “multiple” recurrences, respectively.

Survival Outcomes and Statistical Analysis

Progression-free survival (PFS) and overall survival (OS) were defined as the duration from the date of surgery until the date when tumor progression was diagnosed and death, respectively, or last follow-up. Post-progression survival (PPS) was defined as the time from the first recurrence to either death or last follow-up. Survival time was estimated using the Kaplan–Meier method, and the subgroup differences were compared with log-rank test. Univariate analyses were performed to describe the association between clinical, pathological, and radiological factors and specific patterns of recurrence. For PFS and OS prediction, multivariate logistic regression was conducted on the basis of clinical characteristics and pathological or radiological variables selected by least absolute shrinkage and selection operator (LASSO) logistic regression model. The prediction algorithms were further validated using receiver operating characteristic (ROC) curves. Area under the ROC curve (AUC) and concordance index (C-index) of the multi-marker algorithms were calculated and compared. Two-tailed P < 0.05 were considered statistically significant. All statistical analyses were conducted using R software version 3.2.5 (R Development Core Team; http://www.r-project.org).

Results

Patient Characteristic

From 2008 to 2018, a total of 355 patients underwent radical pancreaticoduodenectomy (PD) or distal pancreatectomy for histologically confirmed PDAC. Excluded from this cohort were 10 patients with microscopic or macroscopic incomplete resection, 12 patients with second primary tumors, and 31 patients with incomplete follow-up information. Consequently, 302 patients were included into this study. All patients were followed up at least 1 year. At the end of follow-up, 195 patients (64.6%) were alive after a median follow-up of 24.7 months (95% confidence interval [CI] 20.3–29.1) from surgery. Recurrence was documented in 173 patients (57.3%), whereas 129 patients (42.7%) had no signs of recurrence. The median follow-up time for patients with and without tumor progression was 13.8 and 40.6 months, respectively.

Timing of Recurrence

Among 173 patients who had recurrence, 18 patients had done so within 6 months, 26 within 6–12 months, 57 within 12–24 months, and 72 beyond 24 months after surgery. There were no significant differences in ages and sexes among patients in different recurrent time groups. Primary tumors in early recurrence groups were larger, more likely to be poorly differentiated, and diagnosed at more advanced local stages. Patients with early recurrence had more often T4 tumors, more metastatic LNs, and more often para-aortic LNs (LN16) metastasis than had those in late recurrence groups (Table 1). Median PFS was 11.8 months (95% CI 10.2–15.3) for the whole cohort and 7.0 months (95% CI 6.2–8.4) for those who developed recurrences. For patients who developed recurrences, the comparisons of PPS and OS stratified by different time intervals of recurrences are shown in Figure 1. It was shown that median OS and PPS for patients who developed recurrences beyond 24 months over surgery (OS, 45.1 months, 95% CI 40.2–52.6; PPS, 17.1 months, 95% CI 11.1–17.5) were significantly longer than for those who had recurrence within 24 months since surgery. Also, patients had similar OS and PPS when their recurrences developed within 6, 6 to 12, or 12 to 24 months since surgery.

TABLE 1

Table 1. Clinicopathological characteristics of patients with PDAC stratified by time of metastases.

FIGURE 1

Figure 1. Post progression survival (A) and overall survival (B) stratified by time period to tumor progression diagnosis counted from the date of surgery.

Patterns of Recurrence

Overall, there were six different patterns of recurrence for all radiological or pathological evidence of progression. Most of patients first recurred at the liver (n = 69, 39.9%), followed by local progression (n = 55, 31.8%), and lung metastases (n = 17, 9.8%). There were 20 (11.6%) patients who had both local and distant progression, and multiple recurrences were observed in 12 (6.9%) patients as the first progression. Liver and lung metastases were the most common distant metastases, compared with the local recurrence, and also contributed to most of the multiple progressions. The proportions of recurrence locations differed significantly at progressive time points. Distribution of these recurrent patterns is shown in Figure 2. Liver-only progressions occupied the majority of all progressions within 6 months, whereas they were responsible for just 12.5% of all recurrences after 24 months since surgery (P < 0.001). Also, liver-only progression diminished over time, and recurrences of other sites became more and more common 1 year later since surgery.

FIGURE 2

Figure 2. Distribution of tumor progression pattern at different time points.

Patients with different progression patterns had significantly different cumulative recurrence rates in different time periods after surgery (Supplement Figure 1). It was shown that cumulative rates of liver metastasis were significantly higher than those of local and other sites of progression, whereas the cumulative rates of liver, lung, and local plus distant and multiple metastases were comparable. The pairwise comparisons of OS (Figure 3), PPS (Supplement Figure 2), and PFS (Figure 4) for patients with different recurrence patterns were conducted. Median OS for patients with local recurrence (29.4 months, 95% CI 24.5–39.6) was significantly longer than that of patients with multiple progressions (17.5 months, 95% CI 11.2–19.5), whereas patients with other recurrence patterns of progression had similar OS rates. Similar results of survival comparisons were observed for PPS. In terms of PFS, patients with local (9.0 months, 95% CI 6.4–10.6) and other sites of progressions (12.7 months, 95% CI 9.1–28.7) had similar median survival, whereas they were both higher than those with other patterns of progressions.

FIGURE 3

Figure 3. Pairwise comparison of overall survival for different tumor progression patterns. (A) Stratification of patients using different progression patterns of progression free, local, liver only, lung only, others, local and distant and multiple progressions. (B–P) Stratification of patients by comparing the following patterns of progression: local vs. liver only, local vs. lung only, local vs. others, local vs. local + distant, local vs. multiple, liver only vs. lung only, liver only vs. others, liver only vs. local + distant, liver only vs. multiple, lung only vs. others, lung only vs. local + distant, lung only vs. multiple, others vs. local + distant, others vs. multiple and local + distant vs. multiple.

FIGURE 4

Figure 4. Pairwise comparison of progression-free survival for different tumor progression patterns. (A) Stratification of patients using different progression patterns of progression free, local, liver only, lung only, others, local and distant and multiple progressions. (B–P) Stratification of patients by comparing the following patterns of progression: local vs. liver only, local vs. lung only, local vs. others, local vs. local + distant, local vs. multiple, liver only vs. lung only, liver only vs. others, liver only vs. local + distant, liver only vs. multiple, lung only vs. others, lung only vs. local + distant, lung only vs. multiple, others vs. local + distant, others vs. multiple and local + distant vs. multiple.

Risk Factors for Different Patterns of Recurrence

Results of univariate and multivariate logistic regression models for local recurrence and liver-only metastasis are shown in Tables 2, 3, respectively. Also, risk factors of lung only, other sites of metastasis, local + distant, and multiple metastases are shown in Supplement Tables 1–4, respectively. Age older than 60 years was a strong predictor for both liver-only metastasis (hazard ratio [HR] = 1.35, 95% CI 1.21–1.73, P = 0.031) and multiple metastases (HR = 9.82, 95% CI 1.20–80.66, P = 0.033). Specific stations of LN metastases were significantly associated with different patterns of progressions, including LN15 metastasis as a predictor for liver-only (HR = 6.39, 95% CI 1.29–31.52, P = 0.023) and local + distant metastases (HR = 8.51, 95% CI 1.27–59.11, P = 0.030), LN18 metastasis as a predictor for local progression (HR = 8.97, 95% CI 1.48–54.23, P = 0.017), LN10 metastasis as a predictor for lung-only metastasis (HR = 15.96, 95% CI 1.89–134.86, P = 0.011), and LN14 metastasis as a predictor for multiple metastases (HR = 7.38, 95% CI 1.61–33.74, P = 0.010). Patients receiving adjuvant chemotherapy had a decreased likelihood of local progression (HR = 0.18, 95% CI 0.08–0.42, P < 0.001) and lung-only metastasis (HR = 0.14, 95% CI 0.02–0.83, P = 0.031) than are those who did not receive adjuvant chemotherapy. Also, PLR was the only independent predictor for other sites of metastases (HR = 0.13, 95% CI 0.02–0.87, P = 0.036), and enlarged imaging LN size was found to increase the likelihood of local + distant metastases (HR = 4.57, 95% CI 1.34–15.60, P = 0.015).

TABLE 2

Table 2. Risk factors for local recurrence in PDAC patients after surgery.

TABLE 3

Table 3. Risk factors for liver metastases in PDAC patients after surgery.

Risk Factors for Progression-Free Survival and Overall Survival

For all included patients, 1-, 2-, and 3-year OS and PFS were 81.3, 58.4, and 47.0% and 49.7, 36.0, and 29.7%, respectively. In order to investigate the prognostic factors of survival, a total of 48 high-dimensional radiological and pathological data were incorporated in the LASSO regression (Figure 5). Three best predictors for OS, including LN16 metastasis, tumor differentiation, and imaging tumor size, and another eight predictors for PFS, including the eighth edition tumor–node–metastasis (TNM) stage, liver-only metastasis, lung-only metastasis, local progression, multiple metastases, LN16 metastasis, imaging tumor size, and LNR, were identified. The predictors selected by LASSO regression, along with the associated clinical factors identified by a univariate analysis, were incorporated to the multivariable analysis. Subsequent analyses illustrated that decreased time interval to progression (HR = 4.30, 95% CI 2.57–7.20, P < 0.001), elevated CA19-9 level (HR = 1.92, 95% CI 1.03–3.58, P = 0.039), and LN16 metastasis (HR = 3.63, 95% CI 1.68–7.82, P = 0.001) were independent predictors for reduced OS (Table 4). Independent prognostic factors for PFS included elevated CEA level (HR = 1.78, 95% CI 1.25–2.53, P = 0.002), local progression (HR = 8.84, 95% CI 5.25–14.87, P < 0.001), liver-only metastasis (HR = 14.74, 95% CI 9.12–23.84, P < 0.001), lung-only metastasis (HR = 9.41, 95% CI 4.45–19.91, P < 0.001), local + distant progression (HR = 11.69, 95% CI 5.79–23.58, P < 0.001), multiple metastases (HR = 19.51, 95% CI 8.78–43.38, P < 0.001), LN16 metastasis (HR = 3.04, 95% CI 1.58–5.99, P < 0.001), imaging tumor size (HR = 1.76, 95% CI 1.16–2.67, P = 0.008), chemotherapy (HR = 0.60, 95% CI 0.42–0.86, P = 0.005), and TNM stage (HR = 2.40, 95% CI 1.17–4.92, P = 0.017) (Table 5).

FIGURE 5

Figure 5. Feature selection using the least absolute shrinkage and selection operator (LASSO) Cox regression model. LASSO coefficient profiles of 48 variables against the log (Lambda) sequence for PFS (A) and tuning parameter (Lambda) selection in the LASSO model used 10-fold cross-validation via minimum criteria for PFS (B). LASSO coefficient profiles of 48 variables against the log (Lambda) sequence for OS (C) and tuning parameter (Lambda) selection in the LASSO model used 10-fold cross-validation via minimum criteria for OS (D). PFS, progression-free survival; OS, overall survival.

TABLE 4

Table 4. Independent prognostic factors for OS.

TABLE 5

Table 5. Independent prognostic factors for PFS.

Performance of Prediction for Overall Survival and Progression-Free Survival

The comparisons of ROC curves of the predictive systems on the basis of the risk factors and TNM stage system are shown in Figure 6. The values of AUC for 1-, 2-, and 3-year OS and PFS prediction were 0.823, 0.844, and 0.858 and 0.789, 0.829, and 0.863, respectively, which were significant higher than those of the TNM stage system (OS, 1 year, 0.614; 2 years, 0.592; 3 years, 0.599; PFS, 1 year, 0.669; 2 years, 0.647; 3 years, 0.630). The predictive system also demonstrated significantly more valuable prediction performance with the C-indexes of 0.829 (95% CI 0.760–0.898) for OS and 0.797 (95% CI 0.723–0. 871) for PFS, respectively, than did the TNM stage system (C-index, OS, 0.588 [95% CI 0.465–0.711]; PFS, 0.619 [95% CI 0.524–0.713]).

FIGURE 6

Figure 6. Comparisons of receiver operating characteristic (ROC) curves of both the predictive system and TNM stage system for predicting 1-, 2-, and 3-year OS (A–C) and PFS (D–F) for LAPC patients after surgery, respectively. TNM, tumor–node–metastasis; PFS, progression-free survival; OS, overall survival; LAPC, locally advanced pancreatic cancer.

Discussion

Pancreatic cancer has an extremely poor prognosis even after surgical resection. Recurrence was observed in more than 60% of all PDAC patients after surgery (4, 16) and remained the main reason of poor prognosis in these patients. In this study, recurrence was observed in 57.3% of patients. In addition, 68.2% of recurrences occurred at a distant site, illustrating that there were systemic diseases in these patients at the time of surgery. Also, 41.6% of recurrences occurred 2 years after surgery. Maybe recurrence-free survival for 2 years did not mean cure, and regular follow-up was also needed for these patients. Furthermore, it was shown that different time intervals or patterns of recurrence would both have different survival. These results suggested that maybe recurrence time interval and patterns were important aspects of recurrence, and the evaluation of factors associated with time intervals and patterns of recurrence opened the door to the exploration of unique biological behaviors of PDAC.

Although the prognostic value of recurrence in survival had been illustrated by previous studies, the current published results differed considerably. For instance, the recurrence rates of PDAC patients after surgery ranged from 38 to 88% in previous studies (17–21). The discrepancies might differ greatly owing to the variations of neoadjuvant treatment regimen and differences of time periods of follow-up. Additionally, the patterns and timing of recurrence of PDAC patients were also not clearly illustrated owing to small population size and limited period of follow-up (22, 23). Moreover, only “local,” “distant,” and “local + distant” groups were analyzed in most of these studies (8, 9, 24–26), and specific recurrence sites were seldom illustrated. In the present study, our detailed recurrence data allowed for further stratification of recurrence patterns in six separated groups: local, liver-only, lung-only, other, local + distant, and multiple metastases. Similar with those in other studies (27, 28), liver-only metastasis and local recurrence contributed to most of the disease progressions. Considering the time period to tumor progress, our study further illustrated that liver-only metastasis occurred mainly in early phase after surgery and diminished over time. Oppositely, other patterns of progressions, including local recurrence and lung metastasis, were more and more common along with time. Following the variations of progression patterns over time, patients might benefit from changes of therapy focus during the period of follow-up.

Progression patterns and time period were two important natural aspects of progression. Apart from the changes of progression patterns over time, it was shown that survival differences were significant when they were stratified by different sites of first recurrence and time periods to tumor progression. In the current study, liver-only metastasis led to the shortest median PFS of only 5.1 months, which was comparable with that of local + distant progression or multiple metastases. Owing to the high rates of occurrence, liver-only metastasis contributed to most of the local + distant progression and multiple metastases. This may partly explain the similar PFS among these three patterns of progression. Similar results were also observed in a study conducted by Suenaga et al. (3), which reported that the median PFS of PDAC patients after surgery was 6.0 months. Apart from liver-only and lung-only metastases, other sites of sole metastasis contributed the longest median PFS (median 12.7 months) among all patterns of progression, followed by local recurrence with a median PFS of 9.0 months. A similar result was also achieved in Vincent's study (4). Additionally, survival differences of OS and PPS were also explored in the present study. Compared with patients with liver-only metastasis, although patients with other sites of distant metastases had slightly short median PPS, they finally achieved longer median OS owing to the significantly extended PFS. Moreover, compared with other patterns of progression, local recurrence contributed to better OS, followed by other sites of sole metastasis, and better PPS, followed by lung-only metastasis. A complete understanding of why local recurrence and lung-only metastasis were associated with relatively favorable PPS remains elusive. A hypothesis assumed that the large capacity of tumor bed and lung allowed patients to endure a greater tumor burden, leading to extended survival (29). Considering the slow growth pattern and apparently less aggressive tumor biology of local progression and lung-only metastasis, maybe locally advanced pancreatic cancer (LAPC) patients can benefit from additional treatment of the subsequent lung and local recurrence after surgery. Additionally, the inherent nature of organ-specific metastasis might be explained by the distinct genetic signatures of both primary PDAC and metastatic lesions. The analysis of biological mechanisms would potentially provide personal therapeutic approaches.

The exploration of risk factors for organ-specific recurrences and predictive factors for survival formed another important finding of this study. Several characteristics were risk factors for liver-only metastasis, such as age older than 60 years and LN15 metastasis. Presence of specific stations of LN metastases could be interpreted as signs of increasing probabilities of progressions. LN18, LN15, and LN14 metastases were identified as predictors for local progression, local + distant metastases, and multiple metastases, respectively. Additionally, as an effective adjuvant therapy to increase survival, the effects of chemotherapy on patterns of progression were poorly understood. Similar with previous study (4), the current study showed that chemotherapy significantly reduced the likelihood of recurrence, especially for local recurrence and lung-only metastasis. Additionally, the prognostic factors were also explored. Apart from the conventional recurrence patterns, elevated levels of CEA, enlarged imaging tumor size, poor differentiation, and advanced TNM stages were all predictive factors of decreased PFS. The exact relation of poor differentiation and poor PFS remained unclear. Maybe this could be partly explained by the ability of PDAC to develop distant metastases, which could be enhanced by the molecules released by the poorly differentiated tumor, including epidermal growth factor, E-cadherin (24). On the other hand, an increasing time prior to tumor progression was also a predictive factor of improved OS, indicating more favorable tumor behavior in patients with late progression. After other risk factors were controlled, the multivariate analysis also illustrated that elevated level of CA19-9 and LN16 metastasis were significantly associated with decreased OS, suggesting that patients with these unfavorable characteristics needed to receive adjuvant therapy after surgery to earn prolonged survival. Similar with study conducted by Groot et al. (30), our results showed that chemotherapy was associated with less local progression and lung-only metastasis and was an independent predictor for PFS. However, the significant associations between chemotherapy and other patterns of recurrences were not observed, and chemotherapy failed to act as a predictor of OS in this study. Owing to the heterogeneity in the length and regimen of the chemotherapy, data on the adjuvant or neoadjuvant chemotherapy in the current literatures were often limited and contradictory. A previous study based on 1,375 patients did not show survival benefit from adjuvant chemotherapy (31), whereas in another study, the additional survival benefit from adjuvant chemotherapy was reported in PDAC patients (32). The selection bias partly contributed to this discrepancy in retrospective study, and maybe more insights concerning the survival benefit of chemotherapy were available from prospective studies.

It is important to note that the precise prediction of progression is essential for the individual treatment. An important advantage of this study was the use of a relatively large cohort to determine the risk factors for different patterns of recurrences and survival. Several independent prognostic factors were selected by evaluating high-dimensional radiological and clinicopathological variables in the current study. In addition, analyses of ROC curves and comparisons of the associated values of AUC and C-indexes of the predictive system and TNM stage system showed a strong predictive strength of the predictive system on the basis of risk factors for OS and PFS. The inclusion of additional clinicopathological variables guaranteed that the established predictive system was better in predicting OS and PFS than did the eighth edition of the TNM stage system. On the other hand, the different clinicopathological features of progression patterns and timing suggested that there might be unique biological features in different progressions. Currently, the molecular feature, SMAD4, was shown to have a close relationship with progression patterns. Tumors with SMAD4 up-regulated tended to be localized, whereas the down-regulation or silence of this gene was likely to promote metastasis (33). Moreover, different regulation of specific genes was associated closely with different patterns of progressions in an animal model (34, 35). Therefore, maybe the combination of clinicopathological characteristics and genetic features would have more meaningful implications in predicting progressions. Clinicians could perform evaluation of recurrence risks and survival on the basis of individual risk factors of patients and specialize the adjuvant therapy, which fitted the current trend to personalized medicine.

This study has several limitations. First, the specific adjuvant therapies after surgery and the associated response to adjuvant therapy were unavailable. More detailed information of length and regimen of chemotherapy would further illustrate the association between therapy and progression. Second, this study only focused on the first recurrence, and subsequent progressions were not taken into accounted. Third, it was well-known that more progressions would be observed over time. In this study, the period of follow-up for all included patients was longer than 1 year, but this time period was not relatively long enough. Although patients were followed up with a median time of 2 years in this study, the whole view of progression in patients could be changed if patients were followed up even longer. A prospective study with an even longer period of follow-up is also needed to validate results of this study. Last, sometimes diagnoses of progression on the basis of imaging were challenging, and it was possible to overestimate the probabilities of progression in PDAC patients after surgery.

In conclusion, for PDAC patients after radical operation, the different patterns and timing of recurrence were accurately described in the present study. This study further identified the risk factors of different recurrence patterns, which could help to predict the occurrence of first tumor progression. Furthermore, individual predictors of OS and PFS were also identified and validated for these patients. These findings further suggested the linkages between different progression patterns and biological heterogeneity, and the exploration might provide new versions into the prediction of tumor progression, prognosis stratification, and a more personalized management for PDAC patients after surgery.

Data Availability Statement

The authenticity of this article has been validated by uploading the key raw data onto the Research Data Deposit public platform (http://www.researchdata.org.cn), with the Approval Number as RDDA2019001267.

Ethics Statement

This study was approved by the Institutional Review Board of SYSUCC. All procedures performed in present study involving human participants were in accordance with the ethical standards of institutional and/or national research committees and the 1964 Declaration of Helsinki and its later amendments or similar ethical standards. Written informed consent for inclusion in this study was obtained from patients prior to treatment.

Author Contributions

SL was responsible for conception, design, and quality control of this study. CH, XH, and YZ performed the study selection, data extraction, statistical analyses and were major contributors in writing the manuscript. CH and XH participated in study selection and statistical analyses. CH, XH, YZ, ZC, and XL contributed in classification criteria discussion. CH, XH, and YZ contributed to the writing of manuscript. SL reviewed and edited the manuscript. All authors have read and approved the final version of the manuscript.

Funding

This work was supported by grants from the National Natural Science Funds (no. 81672390) and the National Key Research and Development Plan (no. 2017YFC0910002).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2019.01197/full#supplementary-material

References

1. Rahib L, Smith BD, Aizenberg R, Rosenzweig AB, Fleshman JM, Matrisian LM. Projecting cancer incidence and deaths to 2030: the unexpected burden of thyroid, liver, and pancreas cancers in the United States. Cancer Res. (2014) 74:2913–21. doi: 10.1158/0008-5472.CAN-14-0155

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Parikh AA, Maiga A, Bentrem D, Squires MH III, Kooby DA, Maithel SK, et al. Adjuvant therapy in pancreas cancer: does it influence patterns of recurrence? J Am Coll Surg. (2016) 222:448–56. doi: 10.1016/j.jamcollsurg.2015.12.031

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Suenaga M, Fujii T, Kanda M, Takami H, Okumura N, Inokawa Y, et al. Pattern of first recurrent lesions in pancreatic cancer: hepatic relapse is associated with dismal prognosis and portal vein invasion. Hepato Gastroenterol. (2014) 61:1756–61. doi: 10.5754/hge13637

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Groot VP, Rezaee N, Wu W, Cameron JL, Fishman EK, Hruban RH, et al. Patterns, timing, and predictors of recurrence following pancreatectomy for pancreatic ductal adenocarcinoma. Ann Surg. (2018) 267:936–45. doi: 10.1097/SLA.0000000000002234

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Ellison LF, Wilkins K. An update on cancer survival. Health Rep. (2010) 21:55–60.

PubMed Abstract | Google Scholar

6. Tummers WS, Groen JV, Sibinga Mulder BG, Farina-Sarasqueta A, Morreau J, Putter H, et al. Impact of resection margin status on recurrence and survival in pancreatic cancer surgery. Br J Surg. (2019) 106:1055–65. doi: 10.1002/bjs.11115

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Kim NH, Kim HJ. Preoperative risk factors for early recurrence in patients with resectable pancreatic ductal adenocarcinoma after curative intent surgical resection. Hepatobiliary Pancreat Dis Int. (2018) 17:450–5. doi: 10.1016/j.hbpd.2018.09.003

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Van den Broeck A, Sergeant G, Ectors N, Van Steenbergen W, Aerts R, Topal B. Patterns of recurrence after curative resection of pancreatic ductal adenocarcinoma. Eur J Surg Oncol. (2009) 35:600–4. doi: 10.1016/j.ejso.2008.12.006

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Sperti C, Pasquali C, Piccoli A, Pedrazzoli S. Recurrence after resection for ductal adenocarcinoma of the pancreas. World J Surg. (1997) 21:195–200. doi: 10.1007/s002689900215

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Chang DK, Johns AL, Merrett ND, Gill AJ, Colvin EK, Scarlett CJ, et al. Margin clearance and outcome in resected pancreatic cancer. J Clin Oncol. (2009) 27:2855–62. doi: 10.1200/JCO.2008.20.5104

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Gebauer F, Tachezy M, Vashist YK, Marx AH, Yekebas E, Izbicki JR, et al. Resection margin clearance in pancreatic cancer after implementation of the Leeds Pathology Protocol (LEEPP): clinically relevant or just academic? World J Surg. (2015) 39:493–9. doi: 10.1007/s00268-014-2808-4

CrossRef Full Text | Google Scholar

12. Tempero MA, Malafa MP, Al-Hawary M, Asbun H, Bain A, Behrman SW, et al. Pancreatic adenocarcinoma, version 2.2017, NCCN clinical practice guidelines in oncology. J Natl Compr Cancer Netw. (2017) 15:1028–61. doi: 10.6004/jnccn.2017.0131

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Amin MB ES, Greene F. AJCC Cancer Staging Manual. 8th ed. Chicago, IL: Springer (2017).

PubMed Abstract

14. He CB, Lin XJ. Inflammation scores predict the survival of patients with hepatocellular carcinoma who were treated with transarterial chemoembolization and recombinant human type-5 adenovirus H101. PLoS ONE. (2017) 12:e0174769. doi: 10.1371/journal.pone.0174769

PubMed Abstract | CrossRef Full Text | Google Scholar

15. De Giorgi U, Procopio G, Giannarelli D, Sabbatini R, Bearz A, Buti S, et al. Association of systemic inflammation index and body mass index with survival in patients with renal cell cancer treated with nivolumab. Clin Cancer Res. (2019) 25:3839–46. doi: 10.1158/1078-0432.CCR-18-3661

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Kim YI, Song KB, Lee YJ, Park KM, Hwang DW, Lee JH, et al. Management of isolated recurrence after surgery for pancreatic adenocarcinoma. Br J Surg. (2019) 106:898–909. doi: 10.1002/bjs.11144

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Katz MH, Pisters PW, Evans DB, Sun CC, Lee JE, Fleming JB, et al. Borderline resectable pancreatic cancer: the importance of this emerging stage of disease. J Am Coll Surg. (2008) 206:833–46; discussion: 46–8. doi: 10.1016/j.jamcollsurg.2007.12.020

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Kadera BE, Sunjaya DB, Isacoff WH, Li L, Hines OJ, Tomlinson JS, et al. Locally advanced pancreatic cancer: association between prolonged preoperative treatment and lymph-node negativity and overall survival. JAMA Surg. (2014) 149:145–53. doi: 10.1001/jamasurg.2013.2690

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Ferrone CR, Marchegiani G, Hong TS, Ryan DP, Deshpande V, McDonnell EI, et al. Radiological and surgical implications of neoadjuvant treatment with FOLFIRINOX for locally advanced and borderline resectable pancreatic cancer. Ann Surg. (2015) 261:12–7. doi: 10.1097/SLA.0000000000000867

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Kim HS, Jang JY, Han Y, Lee KB, Joo I, Lee DH, et al. Survival outcome and prognostic factors of neoadjuvant treatment followed by resection for borderline resectable pancreatic cancer. Ann Surg Treat Res. (2017) 93:186–94. doi: 10.4174/astr.2017.93.4.186

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Jang JY, Han Y, Lee H, Kim SW, Kwon W, Lee KH, et al. Oncological benefits of neoadjuvant chemoradiation with gemcitabine versus upfront surgery in patients with borderline resectable pancreatic cancer: a prospective, randomized, open-label, multicenter phase 2/3 trial. Ann Surg. (2018) 268:215–22. doi: 10.1097/SLA.0000000000002705

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Zhao Y, Wang C. Clinicopathological features, recurrence patterns, and prognosis of pancreatic adenocarcinoma with normal serum CA19–9. A consecutive series of 154 cases from a single institute. J Gastrointest Surg. (2019). doi: 10.1007/s11605-019-04209-w. [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Li S, Xu H, Wang W, Gao H, Li H, Zhang S, et al. The systemic inflammation response index predicts survival and recurrence in patients with resectable pancreatic ductal adenocarcinoma. Cancer Manage Res. (2019) 11:3327–37. doi: 10.2147/CMAR.S197911

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Shibata K, Matsumoto T, Yada K, Sasaki A, Ohta M, Kitano S. Factors predicting recurrence after resection of pancreatic ductal carcinoma. Pancreas. (2005) 31:69–73. doi: 10.1097/01.mpa.0000166998.04266.88

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Zhang Y, Frampton AE, Kyriakides C, Bong JJ, Habib N, Ahmad R, et al. Loco-recurrence after resection for ductal adenocarcinoma of the pancreas: predictors and implications for adjuvant chemoradiotherapy. J Cancer Res Clin Oncol. (2012) 138:1063–71. doi: 10.1007/s00432-012-1165-7

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Asiyanbola B, Gleisner A, Herman JM, Choti MA, Wolfgang CL, Swartz M, et al. Determining pattern of recurrence following pancreaticoduodenectomy and adjuvant 5-flurouracil-based chemoradiation therapy: effect of number of metastatic lymph nodes and lymph node ratio. J Gastrointest Surg. (2009) 13:752–9. doi: 10.1007/s11605-008-0762-x

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Tani M, Kawai M, Miyazawa M, Hirono S, Ina S, Nishioka R, et al. Liver metastasis as an initial recurrence has no impact on the survival of patients with resectable pancreatic adenocarcinoma. Langenbeck's Arch Surg. (2009) 394:249–53. doi: 10.1007/s00423-008-0296-4

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Hernandez JM, Morton CA, Al-Saadi S, Villadolid D, Cooper J, Bowers C, et al. The natural history of resected pancreatic cancer without adjuvant chemotherapy. Am Surg. (2010) 76:480–5.

PubMed Abstract | Google Scholar

29. Zheng B, Ohuchida K, Yan Z, Okumura T, Ohtsuka T, Nakamura M. Primary recurrence in the lung is related to favorable prognosis in patients with pancreatic cancer and postoperative recurrence. World J Surg. (2017) 41:2858–66. doi: 10.1007/s00268-017-4068-6

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Groot VP, Blair AB, Gemenetzis G, Ding D, Burkhart RA, Yu J, et al. Recurrence after neoadjuvant therapy and resection of borderline resectable and locally advanced pancreatic cancer. Eur J Surg Oncol. (2019) 45:1674–83. doi: 10.1016/j.ejso.2019.04.007

PubMed Abstract | CrossRef Full Text | Google Scholar

31. de Geus SWL, Kasumova GG, Eskander MF, Ng SC, Kent TS, James Moser A, et al. Is neoadjuvant therapy sufficient in resected pancreatic cancer patients? A national study. J Gastrointest Surg. (2018) 22:214–25. doi: 10.1007/s11605-017-3541-8

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Roland CL, Katz MH, Tzeng CW, Lin H, Varadhachary GR, Shroff R, et al. The addition of postoperative chemotherapy is associated with improved survival in patients with pancreatic cancer treated with preoperative therapy. Ann Surg Oncol. (2015) 22(Suppl 3):S1221–8. doi: 10.1245/s10434-015-4854-z

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Iacobuzio-Donahue CA, Fu B, Yachida S, Luo M, Abe H, Henderson CM, et al. DPC4 gene status of the primary carcinoma correlates with patterns of failure in patients with pancreatic cancer. J Clin Oncol. (2009) 27:1806–13. doi: 10.1200/JCO.2008.17.7188

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Morales M, Arenas EJ, Urosevic J, Guiu M, Fernández E, Planet E, et al. RARRES3 suppresses breast cancer lung metastasis by regulating adhesion and differentiation. EMBO Mol Med. (2014) 6:865–81. doi: 10.15252/emmm.201303675

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Valiente M, Obenauf AC, Jin X, Chen Q, Zhang XH, Lee DJ, et al. Serpins promote cancer cell survival and vascular co-option in brain metastasis. Cell. (2014) 156:1002–16. doi: 10.1016/j.cell.2014.01.040

PubMed Abstract | CrossRef Full Text | Google Scholar