Perineural invasion as an independent predictor of biochemical recurrence in prostate cancer following radical prostatectomy or radiotherapy: a systematic review and meta-analysis.

BACKGROUND: Although numerous studies have shown that perineural invasion (PNI) is linked to prostate cancer (PCa) risk, the results have been inconsistent. This study aimed to explore the association between PNI and biochemical recurrence (BCR) in patients with PCa following radical prostatectomy (RP) or radiotherapy (RT).
METHODS: According to the PRISMA statement, we searched the PubMed, EMBASE, Chinese National Knowledge Infrastructure (CNKI) and Wan Fang databases from inception to May 2017. Hazard ratios (HRs) and 95% confidence intervals (95% CIs) were extracted from eligible studies. Fixed or random effects model were used to calculate pooled HRs and 95% CIs according to heterogeneity. Publication bias was calculated by Begg's test.
RESULTS: Ultimately, 19 cohort studies that met the eligibility criteria and that involved 13,412 patients (82-2,316 per study) were included in this meta-analysis. The results showed that PNI was associated with higher BCR rates in patients with PCa after RP (HR=1.23, 95% CI: 1.11, 1.36, p<0.001) or RT (HR=1.22, 95% CI: 1.12, 1.34, p<0.001). No potential publication bias was found among the included studies in the RP group (p-Begg = 0.124) or the RT group (p-Begg = 0.081).
CONCLUSIONS: This study suggests that the presence of PNI by histopathology is associated with higher risk of BCR in PCa following RP or RT, and could serve as an independent prognostic factor in patients with PCa.

Background

Prostate cancer (PCa) is the most common newly diagnosed cancer in males, with 1.6 million new cases per year; PCa is the third most common cause of cancer-related death in men [1]. Despite the use of radical prostatectomy (RP) and radiotherapy (RT) as initial therapies for localized PCa, approximately 18 % of patients eventually experience biochemical recurrence (BCR) [2]. Therefore, it is crucial to identify the patients who are at an increased risk of BCR after RP or RT. Preoperative prostate-specific antigen (p-PSA) levels [3], Gleason score (GS) [4] and pathological stage [5] are widely used as traditional risk factors for BCR. However, there is a growing interest in the identification of a new prognostic marker to improve the evaluation of the likelihood of BCR in PCa patients after local treatment.

Perineural invasion (PNI), which is considered a major mechanism for the extraprostatic spread of PCa [6], has been increasingly recognized as a novel prognostic marker [7]. However, whether PNI might be a prognostic factor for BCR is still under debate [8]. Some authors suggest that the presence of PNI is associated with adverse oncological outcomes and a higher risk of BCR, whereas others argue that PNI is not an independent predictor of BCR.

Therefore, to further clarify the relationship between PNI and the risk of BCR in PCa, we performed this systematic review and meta-analysis to evaluate whether the presence of PNI has a prognostic impact on BCR in patients following RP or RT.

Methods

Search strategy

According to the PRISMA guidelines [9], a systematic literature search of the PubMed, EMBASE, Chinese National Knowledge Infrastructure (CNKI) and Wan Fang databases was performed (up to May 2017). The search strategy used the following: (“prostate cancer” or “prostate AND neoplasms”) and (“radical prostatectomy” or “radiotherapy”) and (“perineural invasion”) and (“biochemical recurrence”). We also manually searched for potentially relevant studies from the references listed in the selected review articles. The language of the publications was limited to English and Chinese.

Inclusion and exclusion criteria

The inclusion criteria for the eligible studies were as follows: (1) all patients were diagnosed with histologically confirmed PCa, and PNI in RP specimens and biopsy specimens were assessed by pathologists; (2) all patients underwent RP or RT; (3) BCR after RP was defined as a detectable or rising PSA value after surgery that was ≥0.2 ng/ml with a second confirmation (American Urological Association [10]) ; BCR after RT was defined as a rise in PSA level of ≥2 ng/ml above the nadir (American Society for Radiation Oncology and Radiation Therapy Oncology Group in Phoenix [11]); (4) the risk of BCR was estimated as hazard ratios (HRs) with corresponding 95% confidence intervals (CIs) or the risk could be calculated from the original articles; (5) the study was of a prospective or retrospective cohort design; (6) the articles were published in English or Chinese. The exclusion criteria were as follows: (1) letters, reviews, case reports, editorials and author responses; (2) non-human studies; (3) studies that did not analyze the outcome after PNI and BCR; (4) duplicate articles; (5) articles contained elements that were inconsistent with the inclusion criteria.

Data extraction and Study Quality

The data of the eligible studies were extracted independently by two investigators (Zhen-lei Zha and Wei Qu). Any discrepancies were resolved by discussion with a third reviewer (Hu Zhao). The following data were extracted from the included studies: first author, year of publication, country, recruitment period, sample size, patient’s age, preoperative PSA level, Gleason score, pathological stage, positive percentage of PNI, definition of BCR, follow-up time and the HR(95%CI) of PNI for BCR. When the study provided the results of both the multivariate and univariate outcomes, we chose the multivariate model. The quality of the eligible studies was evaluated according to the Newcastle-Ottawa scale (NOS) [12], which include 3 domains with 8 items. Studies with scores of 7 or more stars were considered high- quality studies.

Statistical analyses

All statistical analyses in this meta-analysis were performed by Stata 12.0 software (Stat Corp, College Station, TX, USA).The association between PNI and BCR outcome was presented as summary relative risk estimates (SRREs) and 95% CIs. Heterogeneity between studies was assessed using Q and I2 statistics. P < 0.10 or I2 > 50% was considered significant heterogeneity. A random effects model was used when obvious heterogeneity was observed, but otherwise, a fixed-effects model was used. Sensitivity analysis was used to estimate the reliability of the pooled results via the sequential omission of each study. Subgroup analyses were performed to examine potential sources of heterogeneity according to the adjusted parameters. Publication bias was assessed by funnel plots and was statistically determined by Begg’s tests; a p value < 0.05 was considered statistically significant.

Results

Search results

The search and selection process for eligible studies is shown in Fig. 1. In all, 222 articles were initially identified. Among them, 73 duplicates were excluded. After the abstracts were screened, 93 articles were excluded for the following reasons: non-human studies, letters, case reports, reviews, and other obvious irrelevant studies. A further 37 articles were excluded after full-text review because the results did not reported PNI on the BCR. Finally, 19 retrospective cohort studies met our inclusion criteria and were included in the meta-analysis. Of these, 13 studies with 10,807patients were analyzed to investigate whether PNI acts as a predictive biomarker of BCR in PCa following RP. Six studies with 2,605 patients were evaluated to determine the relationship between PNI and the risk of BCR in PCa following RT.

Fig. 1

Flowchart of search and inclusion process for eligible studies

Characteristics of the included studies

The main characteristics of the included studies are shown in Tables 1 and 2. All articles included were published in English, except for one, which was published Chinese [13]. All studies were published between 2003 and 2017, and the median duration of follow-up varied from 18.4 to 108 months. Patients in these studies were all diagnosed with PCa and received RP (13 studies) or RT (6 studies). Of the 19studies, 9 originated in Asia, and 10 were conducted in the other regions (USA, Australia, Belgium, Canada). This meta-analysis was based on a total sample size of 13,412 patients, of which 4,197 patients were reported to have PNI. Regarding the GS, 8,306 patients presented with a GS ≥7. For quality evaluated by the NOS, all the studies were found to be of high quality.

Table 1

Characteristics of the included studies on PNI in patients with PCa following RP

Author	Year	Country	Recruitment period	sample size	Age (years)	p-PSA (ng/ml)	GS≥7 (n.%)	Pathological stage ≥ pT2 (n.%)	PNI (n.%)	Follow-up (months)	NOS (score)
Ito et al. [32]	2003	Japan	1989-1998	82	Mean 66.5	Mean 17.2	47 (57.3%)	50(61%)	64 (78%)	Mean 21.7	7
Jeon et al. [25]	2009	Korea	1995-2004	237	Mean 64.5	Mean 11.2	183 (77.2%)	237(100%)	100 (42.2%)	Median 21.6	8
Lee et al. [18]	2010	Korea	1999-2010	361	Mean 69	Mean 15.6	160 (44.3%)	361(100%)	188 (52.1%)	Mean 42.2	7
Loeb et al. [33]	2010	USA	2002-2007	1256	Mean 56.1	Mean 5.7	293 (23.3%)	297(23.6%)	188 (15%)	Mean 33.6	8
Muramaki et al. [34]	2010	Japan	2003-2007	174	Mean 67.5	Mean 9.5	150 (86.2%)	128(73.6%)	107 (83.6%)	Median 44.3	7
Jung et al. [35]	2011	Korea	2005-2009	407	Mean 63.2	Mean 10.3	247 (60.7%)	307(75.4%)	129 (31.7%)	Median 18.4	7
Lee et al. [36]	2015	Korea	2011-2012	752	Mean 66	Mean 9.3	605 (80.5%)	260(34.6%)	483 (64.2%)	Mean 28.3	7
Reeves et al. [27]	2015	Australia	2005-2012	1497	Median 62	NR	1306 (87.2%)	1173 (78%)	1173 (78.4%)	Median 14	7
Watkins et al. [37]	2015	USA	2003-2010	279	Median 62	Median 5.5	80 (28.7%)	279(100%)	193 (69%)	Median 52.7	8
Zhu et al. [13]	2016	China	2002-2014	502	Median 67	Median 12.9	248 (49.4%)	502(100%)	411 (81.9%)	Median 28.6	7
Kang et al. [24]	2016	Korea	2003-2014	2034	Median 67	Median 7.4	1726 (84.9%)	NR	252 (12.4%)	Median 48	7
Kim et al. [38]	2016	Korea	2005-2016	2316	Median 65	Median 7.4	1621 (69.9%)	2316(100%)	NR	Median 60	8
Aoun et al. [5]	2017	Belgium	1990-2013	910	Mean 64	Median 6.8	189 (20.8%)	910(100%)	305 (33.5%)	Median 108	7

GS: Gleason score; PNI: perineural invasion; NR: not reported;NOS: Newcastle–Ottawa Scale; p-PSA: preoperative PSA

Table 2

Characteristics of the included studies on PNI in patients with PCa following RT

Author	Year	Country	Recruitment period	sample size	Age (years)	p-PSA (ng/ml)	GS≥7 (n.%)	PNI (n.%)	Radiotherapy	Follow-up (months)	NOS (score)
Copp et al[39]	2005	USA	1997-2000	93	Median 69.1	Median 9.1	34(36.6%)	17(18.3%)	EBRT	Median 45	7
Yu et al[19]	2007	USA	1993-2002	657	Mean 68	NR	238(36.2%)	112 (19.1%)	EBRT	Median 68	7
Pina et al[40]	2010	Canada	1999-2008	339	Mean 61.9	Mean 5.7	45(13.3%)	89(26.3%)	BT	Median 32	8
Feng et al[41]	2011	USA	1998-2008	718	Median 69.4	Median 8.3	474(66%)	220(34%)	EBRT	Median 62.2	8
Ding et al[30]	2012	USA	1998-2003	185	Median 69	NR	47(25.4%)	10(5%)	BT/ EBRT	Median 80	7
Spratt et al[42]	2013	USA	1989-2011	613	Median 70	NR	613(100%)	156(25.4%)	EBRT	Median 72	7

GS: Gleason score; PNI: perineural invasion; EBRT: external beam radiotherapy; BT: brachytherapy; NR: not reported; NOS: Newcastle–Ottawa Scale; p-PSA: preoperative PSA

Relationship between PNI and BCR after RP

The forest plots of the meta-analyses are shown in Fig. 2. A random effects model was applied because the heterogeneity was evident among these studies (p= 0.025, I2= 48.6%). The pooled HR indicated that the presence of PNI was associated with a higher risk of BCR in patients with PCa after RP (HR=1.23, 95% CI: 1.11, 1.36, p<0.001). According to the sensitivity analysis, the SRRE ranged from 1.18 (95% CI: 1.09, 1.29) to 1.26 (95% CI: 1.14, 1.38) (Fig. 4a). In the subgroup analyses, we found that the heterogeneity was decreased significantly in some models, such as those of geographical region (“other regions”) and, age <65 years. It should also be noted that; the results of the subgroup analyses were consistent with the primary findings (Table 3). Furthermore, no statistical evidence of publication bias was found as assessed by Begg’s tests (p = 0.124) ( Fig. 5a).

Fig. 2

Forest plots for PNI and outcomes of BCR in PCa patients following RP

Fig. 4

a: Sensitivity analysis for PNI and BCR in PCa following RP; b: Sensitivity analysis for PNI and BCR in PCa following RT

Table 3

Summary and subgroup analyses of the eligible studies

Analysis specification	No. of studies	Study heterogeneity		Effects model	Pooled HR(95% CI)	P-Value
		I2(%)	P heterogeneity
Radical prostatectomy
Overall	13	48.6	0.025	Random	1.23 (1.11, 1.36)	<0.001
Geographical region
Asia	9	57.8	0.013	Random	1.26 (1.10, 1.45)	0.001
Other regions	4	6	0.363	Fixed	1.23 (1.11, 1.36)	0.027
Age(years)
≥65	7	67.2	0.006	Random	1.25 (1.04, 1.50)	0.015
<65	6	0	0.425	Fixed	1.20 (1.08, 1.32)	<0.001
p-PSA(ng/ml)
≥10	4	44.4	0.126	Fixed	1.21 (1.01, 1.46)	0.037
<10	8	57.8	0.027	Random	1.28 (1.10, 1.49)	0.001
Sample size(cases)
≥ 500	7	45.7	0.101	Fixed	1.23 (1.03, 1.47)	0.002
< 500	6	57.5	0.028	Random	1.24 (1.08, 1.41)	0.026
Median follow-up( months)
≥ 30	7	34.9	0.162	Fixed	1.34 (1.12, 1.61)	0.001
< 30	6	57.9	0.036	Random	1.16 (1.11, 1.36)	0.014
Radiotherapy
Overall	6	34.7	0.176	Fixed	1.22 (1.12, 1.34)	<0.001
Age(years)
≥65	4	44.8	0.143	Fixed	1.21 (1.11, 1.33)	<0.001
<65	2	38.8	0.201	Fixed	1.45 (0.93, 2.24)	0.097
Sample size(cases)
≥ 500	3	0	0.495	Fixed	1.20 (1.09, 1.31)	<0.001
< 500	3	31.0	0.235	Fixed	1.68 (1.18, 2.39)	0.004

Fig. 5

a: Funnel plots and Begg’s tests for the evaluation of publication bias in RP group; b: Funnel plots and Begg’s tests for the evaluation of publication bias in RT group

Association between PNI and BCR following RT

As shown in Fig. 3, the pooled HR from 6 studies indicated that PNI was an independent risk factor for BCA in PCa following RT (HR=1.22, 95% CI: 1.12, 1.34, p<0.001). Due to the lack of evidence of heterogeneity among the studies (p= 0.176, I2 = 34.7%), a fixed effects model was applied. The SRRE for BCR ranged from1.20 (95% CI: 1.07, 1.34) to 1.27(95% CI: 1.15, 1.42) according to the sensitivity analysis (Fig. 4b). The results of the subgroup analysis showed that PNI in a sample size ≥ 500 cases was significantly associated with BCR (n=3, HR=1.20, 95% CI: 1.09, 1.31, p <0.001) (Table 3). Begg’s tests (p = 0.081) provided no evidence of substantial publication bias in these studies (Fig. 5b).

Fig. 3

Forest plots for PNI and outcomes of BCR in PCa patients following RT

Discussion

With the wide spread use of the serum PSA level in PCa screening, most patients with newly diagnosed PCa present with clinically localized or locally advanced disease [14]. RP and RT have become the standard local treatments for localized PCa. However, approximately 1/3 of patients will develop BCR after RP, and 25~33% will experience BCR after RT [15]. Therefore, the characterization of the pathological features that may predict BCR is important for counseling patients and in directing the initial or adjuvant therapy. To date, seminal vesicle [16], surgical margins and extracapsular extension [17] are widely used as risk predictors of BCR in PCa. However, these are inferior and cannot accurately predict the risk of BCR in patients with PCa and complicated tumor backgrounds. Therefore, more sensitive and specific markers of risk are needed.

The phenomenon of peripheral nerve involvement has been overlooked for a long time, but it is now gaining recognition as a potential component of the cancer microenvironment [18]. PNI is defined as cancer cell tracking along or around a nerve within the perineural space [19]. It has been demonstrated that PNI may be a route of metastasis for many different cancers (pancreatic [20], bladder [21] and colorectal [22]). Patients with PNI have a higher risk of extraprostatic extension detected at the time of RP [23]; in other words, these patients have an increased risk of BCR.

However, the clinical significance of PNI in PCa following RP or RC is still controversial. Kang [24] and Jeon [25] showed that PNI is an adverse pathologic parameter and an independent predictor for BCR in PCa patients who undergo RP. Similarly, Yu [19] and Wong [26] considered that PNI is an independent risk factor associated with an increased risk of BCR in patients who undergo external beam radiotherapy. On the contrary, Reeves [27] and Freedland [28] reported that PNI is not correlated with extracapsular extension and BCR in PCa after RP. Nevertheless, Weight [29] suggested that the presence of PNI is not a significant predictor of BCR in patients undergoing brachytherapy (BT) for PCa, and Ding[30] demonstrated a significant independent association between PNI and an increased risk of biochemical failure in 185 PCa patients who received BT.

The conflicting results from different studies prompted us to perform this meta-analysis. For this study, we analyzed 19 studies that included 13,412 patients with PCa, of which 4,197 (31.2%) had PNI. Our results demonstrate that PNI was associated with a higher risk of BCR in PCa patients who underwent RP (p <0.001) or RT (p <0.001). These findings are similar to those reported by Kang [24] and Yu [19]. Notably, the overall findings in the present meta-analysis were consistently independent of geographical region, age, p-PSA, sample size, publication year, and follow-up in the RP group; the findings were also independent of age and sample size in the RT series. In the sensitivity analysis, we sequentially excluded each study, and the reliability and robustness of the results were confirmed. Moreover, no evidence of significant publication bias was found in this analysis according to Funnel plots and Begg’s tests. Although subgroup analyses were conducted in the RP and RT groups, significant heterogeneity still existed in some subgroups. The relevant results indicated that these factors were not potential sources of high heterogeneity.

Several potential limitations of our meta-analysis need to be addressed. First, all the included studies were retrospective cohort studies, and data extracted from those studies may have led to inherent potential bias as an unmeasured and uncontrolled confounder. Second, only papers in English and Chinese were selected, and thus some studies that were published in other languages may have otherwise been eligible. Therefore, both selection and publication bias are possible. Third, the pathological diagnosis of PCa and the detection method of PNI varied throughout the eligible studies. Collectively, these factors may promote significant heterogeneity. Fourth, substantial heterogeneity was identified across the studies, and although subgroup analyses were conducted to explore the source of the heterogeneity, heterogeneity still existed in certain categories of analysis. Thus, our results should be interpreted with caution.

Conclusions

To conclude, despite the limitations listed above, this meta-analysis suggested the prognostic and clinicopathological importance of PNI in PCa. The results demonstrated that the presence of PNI is associated with a high risk of BCR whether the patient undergoes RP or RT. Since BCR has been reported to lead to distant metastasis and cancer death [31], we suggested that some patients with PCa and PNI may benefit from adjuvant local or systemic therapy. However, more randomized controlled trials with standardized methods and long-term follow-up are needed to verify our results.