Biomarker implication of kallikrein-related peptidases as prognostic tissue substrates of poor survival in colorectal cancer

Background Recent studies have demonstrated that the kallikrein and kallikrein-related peptidases (KLKs) exhibit aberrant expression in patients with colorectal cancer (CRC) and might be considered as potential prognostic biomarkers of CRC. However, inconsistent findings have been reported, which promote us to summarize the global prognostic roles of KLKs for survival in CRC patients. Methods Eligible published studies were identified by searching electronic databases with several search strategies. The patients’ baseline characteristics and survival results were extracted from enrolled studies and pooled as combined hazard ratio (HR) with 95% confidence interval (95% CI) to estimate the effect size. Results A total of 25 and 22 eligible studies were included in the meta-analysis to evaluate the prognostic roles of KLKs on overall survival (OS) and disease-free survival (DFS), respectively. KLKs overexpression was significantly associated with worse OS (pooled HR = 1.43, 95% CI 1.27–1.60, P < 0.001) and short DFS (pooled HR = 1.35, 95% CI 1.21–1.51, P < 0.001). Importantly, subgroup and meta-regression analyses revealed the survival differences among different races and detection methods of KLKs. Furthermore, several specific members of KLKs were identified to be more significantly related to worse OS and DFS compared with other members. Conclusion The present study demonstrated that KLKs may have the potential to serve as promising biomarkers to monitor CRC prognosis and progression. The promising results concerning the utility of KLKs in clinical practice encourage the further investigation of their clinical utility applicability as tumor markers of CRC.

hydrolyze peptide bonds, playing important and vital role in cell growth regulation, invasion, and angiogenesis [4]. Among all serine proteases within the human genome, the kallikrein and kallikrein-related peptidases (KLKs) clusters are the largest [5]. The family of tissue KLKs consists of 15 secreted serine proteases encoded by conserved genes (KLK1-KLK15), which are localized in tandem on chromosomal region 19q13. 4. Accumulating evidence has indicated that many members of the human tissue KLKs are differentially expressed in many pathological conditions and a number of malignancies and may have clinical utility as cancer diagnostic/prognostic biomarkers [6]. Since then, the prognostic significance of several KLKs in CRC has been intensively investigated. Protein expression analysis of a panel of KLKs in cytosolic extracts from CRC tissues revealed that KLKs are aberrantly expressed in colorectal tumors, compared with their noncancerous counterparts [7]. Given their critical involvement in the vital biological processes and unique biomarker features mentioned above, KLKs could be considered as good candidates in increasing the accuracy of prediction of patients' survival beyond the traditional clinical information.
Although some of the studies evaluated the prognostic value of KLKs in CRC patients, the relationships between KLKs and CRC remain controversial as individual studies were not comprehensive for they involved only small study populations. For instance, Christodoulou et al. [8] showed that increased level of KLK6 was associated with worse overall survival of CRC patients; nevertheless, this association has not been detected in the study illustrated by Vakrakou et al. [9]. Therefore, in the present study, we conducted this comprehensive meta-analysis to compare the survival outcome between CRC patients with high levels of KLKs and those with low levels of KLKs. We aimed to overcome the limitation of the single study and to obtain a better understanding of the prognostic value of KLKs in CRC.

Publication search
The present study was conducted and reported under the guidelines formulated in Preferred Reporting Items for Systematic Reviews and Meta-analyses (PRISMA). A comprehensive literature search was carried out based on the electronic databases including PubMed, EMBASE, Cochrane Library, and Web of Science databases (up to March 12, 2019) by using the following keywords: ("Kallikrein"OR "KLK" OR "Kallikrein-related peptidase"), ("rectal" OR "rectum" OR "colon" OR "colorectal" OR "CRC"), and ("cancer" OR "tumor" OR "neoplasm" OR "carcinoma"). All potentially eligible studies were identified and their bibliographies were carefully examined to retrieve other eligible studies.
In case of omission, additional relevant studies were identified by scanning the references cited in the original studies. Two reviewers (Peng and Shen) independently performed the publication search and the following steps.

Inclusion criteria
The studies qualified to be included had to meet the following criteria: (1) they investigated the relationships between KLKs expression and CRC prognosis; (2) they reported survival data for OS or DFS; (3) they directly provided HRs with 95% CIs or providing adequate statistics to conjecture HRs with a corresponding 95% CI.

Exclusion criteria
The studies were excluded if (1) they were not pertinent to KLKs; (2) they published as reviews, letters, case reports, editorials, or expert opinions; (3) they were non-English publications; or (4) they lacked sufficient data for further quantification.

Data extraction
According to the inclusion and exclusion criteria above, data were collected carefully and independently by two researchers (Peng and Shen) from all eligible publications based on standardized forms. Any disagreement between the researchers was resolved by consulting with a third investigator (Zhao) through independently extracting data from the enrolled publication and then reaching a consensus by discussions. The following characteristics from each study were extracted: first author, publication year, study population, patient characteristics (age, gender, cancer type, etc.), number of patients, methods of KLKs detection, prognostic results including follow-up time and HRs estimates with 95% CIs for DFS or OS. If the survival data (HRs and 95% CIs) were not directly reported in the original study, they were obtained from Kaplan-Meier survival curves with the Engauge Digitizer V4.1 and estimated using the method previously introduced by Tierney et al. [10]. We also contacted the authors of eligible studies by email for additional information and the essential data required for the meta-analytic calculations.

Quality assessment
The methodological quality assessment of each eligible article was conducted following the guidelines of the Newcastle-Ottawa Scale (NOS), which assessed studies with 9 items including the selection of the patient population, study comparability, outcome of interest, followup et al. [11]. Studies with an NOS score ranging from 6 to 8 were considered of high quality.

Statistical analysis
In order to assess the associations between KLKs expression and the survival of CRC patients, pooled HRs with the corresponding 95% CIs were used to assess the strength of the associations between KLKs and the clinical prognosis of CRC patients. Heterogeneity across studies was checked using Cochran's Q test (significant at P < 0.05) and Higgins's I 2 statistic (ranging from 0 to 100%) [12]. For the presence of heterogeneity (P < 0.05, I 2 > 50%), a random-effect model was employed to calculate the pooled HRs and 95% CIs; otherwise, a fixed effect model was selected (P > 0.05, I 2 < 50%). Potential sources of heterogeneity were explored by performing meta-regression, subgroup, and sensitivity analyses [13]. At last, Begg's funnel plots and Egger's test were utilized to assess the included studies for the possible publication bias [14]. The statistical analyses were carried out using STATA (version 14.0) statistical software. Values of P < 0.05 were considered as statistical significance, except those for heterogeneity.

Demographic characteristics
Using different combinations of key terms, the initial search from the selected literature databases (PubMed, EMBASE, Cochrane Library, and Web of Science databases) and other sources retrieved a total of 297 records. As shown in the flow diagram for the literature (Fig. 1), after careful exclusion of inappropriate ones in each step, 14 articles including 25 studies for OS and 22 studies for DFS that met the inclusion norm were finally enrolled for the evidence synthesis, which evaluated the relevance between KLKs expression and CRC prognosis [8,9,[15][16][17][18][19][20][21][22][23][24][25]. The main features of all the eligible studies are summarized in Tables 1 and 2. Among all cohorts, Caucasian (20 studies both for OS and DFS) became the major race of literatures, followed by Asian (5 for OS and 2 for DFS). Three methods were applied to measure the expression including quantitative real-time polymerase chain reaction (qRT-PCR), immunohistochemistry (IHC) and enzyme-linked immunosorbent assay (ELISA). Several members of KLKs were evaluated by more studies    The Newcastle-Ottawa Scale (NOS) score of each study included for the evaluation of OS and DFS ranged from 7 to 9, which indicated that the quality of the included studies was moderate to high.

Impact of KLKs expression on DFS
A random-effects model was applied to estimate the pooled HR and corresponding 95% CI as the heterogeneity test reported the P value of 0.001 and I 2 values of 53.8%. The pooled result revealed that elevated KLKs expression was significantly associated with worse DFS of patients with CRC (pooled HR = 1.35, 95% CI = 1.21-1.51, P < 0.001) (Fig. 2).
Meta-regression analysis was also performed to explore the potential factors of the heterogeneity. We considered 4 covariates (ethnicity, sample size, KLKs member and detection method) may contribute to the heterogeneity. The results revealed that neither sample size, nor KLKs member was the source of heterogeneity, but the ethnicity populations and detection methods have influence on the pooled results (P < 0.05).

Impact of KLKs expression on OS
For studies evaluating OS, moderate heterogeneity across studies was also observed (I 2 = 51.9%, P = 0.001). Therefore, we also calculated the combined HR and the corresponding 95% CI based on a random model. According to the final pooled results, a significant correlation between KLKs overexpression and a worse OS was shown in patients with CRC with the pooled HR of 1.43 (95% CI 1.27-1.60, P < 0.001).
Then, a meta-regression analysis was carried out for the sources of the heterogeneity. It was indicated from the results that ethnicity and detection methods may have contributed to the heterogeneity.

Publication bias and sensitivity analysis
Begg's funnel plot and Egger's test were utilized to assess the presence of publication bias among the included literatures (Fig. 3). The funnel plots pointed out symmetry for all enrolled studies and Deeks' test revealed potential heterogeneity in the present analysis about DFS and OS (P < 0.001). Then, we applied the trim and fill method to further investigate the publication bias and make the pooled HR more reliable both with the P value less than 0.01. Sensitivity analyses were further carried out to evaluate the stability of the pooled results for DFS and OS (Fig. 4). However, no specific study dominated the evidence synthesis as the removal of any individual study had no significant influence on the pooled results of DFS and OS, which indicates that the conclusions from our meta-analysis were relatively reliable.

Discussion
Accumulating evidence has indicated that elevated KLKs promote cancer progression and predict poor prognosis of CRC patients. A series of quantitative analyses were conducted to investigate the prognostic value of KLKs overexpression in CRC. However, the sample sizes in most studies are small. Besides, it is inconclusive about the association between KLKs expression and progression of CRC. Thus, we performed this comprehensive and up-to-date research to draw a complete overview of To the best of our knowledge, the current study was the first systematic evaluation of the literatures to investigate the prognostic role of KLKs expression in CRC patients. We evaluated survival data from 25 different studies assessing OS and from 22 different studies estimating DFS. Our results suggest that the elevated level of KLKs is indeed a poor prognostic biomarker for CRC in OS (pooled HR = 1.43, 95% CI 1.27-1.60, P < 0.001) and  In particular, the predictive roles for DFS and OS were more significant in Asians than in Caucasians. Additionally, when data was stratified according to detection methods, the results showed that the prognostic value of KLKs over-expression was significant in mRNA level and the protein level detected by qRT-PCR, IHC and ELISA, respectively. What's more, several members of KLKs evaluated by more studies showed great promise for survival prediction including KLK6, KLK7 and KLK10. Of course, they along with the other members of KLKs worth further study. Moderate heterogeneity remained in the meta-analyses of the data for DFS and OS, which can be explained by ethnicity and detection methods in the meta-regression. Sensitivity analysis failed to identify any deviated study, indicating that there was high robustness in our meta-analysis.
Recently, accumulating researches have indicated KLKs to be potential predictors for CRC prognosis. Our data also revealed that KLKs are very promising for survival outcome prediction. However, several factors should be considered to be the necessary future directions of the application of KLKs in clinical practice. Firstly, an appropriate definition should be made about the standard cut-off value of KLKs level for increased survival risk. To a large extent, the methodological inconsistency contributed to the divergence of contemporary findings about the prognostic value of KLKs. Most researchers prefer median or mean value in their study as the approach to setting cut-off value of KLKs expression varied among different studies. Determination of standard patterns of KLKs expression will significantly prompt achievement of final consensus about the prognostic value of KLKs. Secondly, which should be used, the protein level or mRNA level? IHC, qRT-PCR and ELISA are all widely selected, detecting different levels of KLKs, respectively. Our results in this study also indicated that KLKs measured by these three methods all could be used to predict the survival outcome of CRC. Thirdly, which is better for clinical application, a single KLK or a panel of KLKs? CRC is a complex disease; single or limited biomarker is unlikely to reveal the complicated evolutionary process at the systemic level and lack of specificity and sensitivity. Our preliminary results show that combination biomarkers may be more reliable with greater power as they help to explain the internal mechanisms of CRC as well as the external factors influencing it [26,27]. As a future perspective, combination of KLKs may be considered for the further improvement of the prognostic role and large-scale prospective studies are still necessary for further validation.
The associations between KLKs expression and cancer prognosis may be partly caused by the biological function of the KLKs. KLKs are secreted serine proteases with distinct expression patterns and physiological functions in several systems, especially in the digestive system. KLKs have been demonstrated to take part in numerous physiological processes such as cell growth regulation, angiogenesis, invasion, and metastasis. Aberrant expression of KLK family members is highly associated with various clinic-pathological parameters of patients suffered with colorectal, gastric, pancreatic, hepatic, and esophageal cancer. Accumulating evidence has revealed that KLKs facilitate CRC progression due to their ability to degrade extracellular matrix proteins, thereby promoting tumor invasion as well as metastasis [28]. A large number of studies have also demonstrated that KLKs play vital roles in regulating proteinase-activated receptors (PARs), which is a protein family containing four G-protein-coupled receptor members including PAR1, PAR2, PAR3, and PAR4 [29,30]. Recent findings have indicated that KLKs has been inextricably linked to the cleavage of PARs, which may bring about coupling of the receptors to heterotrimeric G proteins, thereby generate signal transduction and thus enhance tumor cell proliferation [31]. It is also important to note that KLKs may play a significant part in uncoupling of the receptor from the signal transduction pathway with a strong implication on the initiation and progression of CRC as they are able to separate downstream of the PAR activation site or within an extracellular loop [32]. According to the previous evidence from studies on molecular mechanism, it is not difficult to understand why KLKs may become promising biomarkers for CRC prognosis prediction.
There are several important strengths from results of the current study. First, KLKs were confirmed to act as reliable prognostic biomarkers for CRC. Our data indicated that CRC patients with elevated expression levels of KLKs may suffer from an increased risk of poor survival, which was 1.43-fold higher for OS and 1.35-fold higher for DFS when compared with CRC patients with low KLKs expression. Second, we demonstrated that high KLKs expression correlated with poor OS and DFS both in Asian and Caucasian patients and this finding may be extended to other ethnic groups. Moreover, since some members of KLKs have been identified to have a biomarker role in the prognosis of CRC, future studies should be studied and uncovered the roles of other members of KLKs in CRC. In addition, as the easiest detection method, qRT-PCR may be recommended as the first choice since the prognostic value of KLKs by qRT-PCR was most significant in our study compared with ELISA and IHC. Finally, it underlines the potential to develop KLKs as a promising therapeutic target and prognostic biomarker for CRC. In the future, the KLKs may be accepted for clinical application like the common used tumor markers such as CEA, CA125, and CA199.
Apart from the inspiring outcomes, our conclusions should be interpreted with caution since there were also several limitations in the present work. To begin with, most of the included studies were designed as retrospective studies, which inherently contained greater potential for confounding than do randomised controlled trials. Next, TNM stage of the enrolled studies is a potent selection bias factor. Without patient-level data, it would be very difficult to tell whether KLKs are actually useful as a prognostic factor. A patient-level meta-analysis can be used to analyze all of the data in a consistent manner and includes data from unpublished studies. Moreover, only Asians and Caucasians were in the meta-analysis, no African population included in the analysis, which may cause potential heterogeneity from ethnicity. In addition, the numbers of studied were inconsistent among the KLKs including KLK1 (n = 5), KLK2 (n = 5), KLK4 (n = 1) and KLK5 (n = 1). Accordingly, subgroup analysis by specific member of KLKs could not be performed for the limited individual sample size. We ignored the variety of KLKs in the study. The meta-analysis was performed with studies of the selected KLKs, regardless of their types. Furthermore, the potential publication bias may exist in our analysis. At last, the current study is not registered and there may be a small offset, but we still strictly follow the steps of the systematic review.

Conclusion
Taken together, in this study, it is concluded that tissue KLKs may be effectively predictive biomarkers for CRC prognosis. In future, more clinical studies are warranted to confirm the prognostic role of KLKs before its practical implementation in management of CRC.
Abbreviations CRC : Colorectal cancer; HR: Hazard ratio; CI: Confidence interval; PRISMA: Preferred Reporting Items for Systematic Reviews and Meta-analyses; NOS: Newcastle-Ottawa Scale.