KLK8, 10, and 11 mRNA expression in tumor tissues of triple-negative breast cancer (TNBC) patients and their correlation with clinical parameters
KLK8, 10, and 11 mRNA expression levels were quantified by established qPCR assays [25,33] in 123 cases of triple-negative breast cancer tissues. KLK8 mRNA levels ranged from 0.00 to 19.34 (median: 0.21), KLK10 mRNA expression ranged from 0.00 to 10.53 (median: 0.15), and KLK11 mRNA expression ranged from 0.00 to 2.70 (median: 0.01, a supplementary figure file shows this in more detail [see Suppl. file 1]).
mRNA expression of KLK8 and KLK11 was categorized by the 50th percentile (median) into a low-expressing versus a high-expressing group, and KLK10 by the 67th percentile in a low expression group (tertiles 1+2) versus a high expression group (tertile 3).
By applying Spearman correlation analysis, a strong positive correlation was observed between KLK8 and KLK10 (rs = 0.612, p < 0.001) as well as between KLK10 and KLK11 (rs = 0.722, p < 0.001), whereas the correlation between KLK8 and KLK11 mRNA was less pronounced (rs = 0.460, p < 0.001; an additional figure file shows this in more detail [see Suppl. file 2]).
The associations between low and high expression groups of the three KLK pairs were also validated by box plot analysis (Mann-Whitney test; p-value in all cases < 0.001, Fig. 1) indicating co-regulation of the three KLK genes in TNBC. Due to this observation, we further dichotomized the mRNA levels into either a KLK8+KLK10, KLK8+KLK11, KLK10+KLK11, or KLK8+KLK10+KLK11 low-expressing group (mRNA levels below the respective cutoff value) versus high-expressing groups (KLK8 and/or KLK10 and/or KLK11 mRNA levels above the respective cutoff value) for statistical analyses.
The associations between relative KLK mRNA expression levels and the established clinical variables, including age, lymph node status, and tumor size, are summarized in Table 1. Neither KLK8 nor KLK10 nor KLK11 mRNA expression levels were associated with clinicopathological parameters of TNBC. Similarly, no statistically significant associations were observed between KLK combinations and clinical parameters (see Suppl. Table 1).
Table 1 Association between KLKs mRNA expression levels and clinicopathological parameters in patients with triple-negative breast cancer
Clinicopathological parameters
|
KLK8a
|
KLK10a
|
KLK11a
|
low/high
|
low/high
|
low/high
|
Age
|
p=0.907
|
p=0.294
|
p=0.761
|
≤ 60 years
|
31/30
|
43/25
|
33/28
|
> 60 years
|
27/25
|
39/15
|
23/22
|
Lymph node status
|
p=0.756
|
p=0.990
|
p=0.473
|
N0
|
31/31
|
45/22
|
33/26
|
N1/N2/N3
|
27/24
|
37/18
|
23/24
|
Tumor Size
|
p=0.843
|
p=0.207
|
p=0.665
|
≤20 mm
|
16/14
|
25/8
|
15/15
|
>20 mm
|
42/40
|
56/32
|
41/34
|
a Chi-square test (cut-off point: KLK8=50th percentile, KLK10=67th percentile, KLK11=50th percentile).
Due to missing values, numbers do not always add up to n = 123.
Association of KLK8, 10, and 11 mRNA expression levels with disease-free survival (DFS) and overall survival (OS) in univariate Cox regression analysis
Univariate Cox regression analysis was performed to evaluate the prognostic relevance of KLKs, the combination thereof, and the clinicopathological parameters regarding the patient outcome in the TNBC cohort (Table 2). Among the clinical variables, advanced age and a positive lymph node status indicated a significantly shorter DFS and OS.
Table 2 Univariate Cox regression analysis of clinical outcome in triple-negative breast cancer for tumor biological factors.
Clinicopathological
parameters
|
DFS
DFS
|
OS
OS
|
Noa
|
HR (95% CI)b
|
p
|
Noa
|
HR (95% CI)b
|
p
|
Age
|
|
|
0.005
|
|
|
< 0.001
|
≤ 60 years
|
65
|
1
|
|
65
|
1
|
|
> 60 years
|
52
|
2.09 (1.25-3.52)
|
|
56
|
2.76 (1.62-4.70)
|
|
Lymph node status
|
|
|
0.023
|
|
|
0.018
|
N0
|
63
|
1
|
|
65
|
1
|
|
N1/N2/N3
|
54
|
1.82 (1.09-3.05)
|
|
56
|
1.86 (1.11-3.10)
|
|
Tumor Size
|
|
|
0.058
|
|
|
0.076
|
≤20 mm
|
32
|
1
|
|
33
|
1
|
|
>20 mm
|
84
|
1.89 (0.98-3.64)
|
|
87
|
1.85 (0.94-3.66)
|
|
KLK8 mRNAc
|
|
|
0.037
|
|
|
0.313
|
low
|
54
|
1
|
|
57
|
1
|
|
high
|
54
|
1.88 (1.04-3.41)
|
|
54
|
1.37 (0.74-2.52)
|
|
KLK10 mRNAc
|
|
|
0.045
|
|
|
0.076
|
low
|
78
|
1
|
|
81
|
1
|
|
high
|
37
|
1.78 (1.01-3.12)
|
|
38
|
1.68 (0.95-2.99)
|
|
KLK11 mRNAc
|
|
|
0.053
|
|
|
0.016
|
low
|
52
|
1
|
|
54
|
1
|
|
high
|
50
|
1.85 (0.99-3.44)
|
|
50
|
2.29 (1.17-4.49)
|
|
KLK8+KLK10d
|
|
|
0.013
|
|
|
0.164
|
low
|
47
|
1
|
|
50
|
1
|
|
high
|
60
|
2.28 (1.19-4.34)
|
|
60
|
1.57 (0.83-2.95)
|
|
KLK10+KLK11e
|
|
|
0.046
|
|
|
0.035
|
low
|
46
|
1
|
|
48
|
1
|
|
high
|
54
|
1.95 (1.01-3.78)
|
|
54
|
2.12 (1.05-4.28)
|
|
KLK8+KLK11f
|
|
|
0.042
|
|
|
0.147
|
low
|
31
|
1
|
|
33
|
1
|
|
high
|
66
|
2.18 (1.03-4.61)
|
|
66
|
1.76 (0.82-3.79)
|
|
KLK8+KLK10+KLK11g
|
|
|
0.019
|
|
|
0.106
|
low
|
27
|
1
|
|
29
|
1
|
|
high
|
69
|
2.85 (1.19-6.86)
|
|
69
|
1.99 (0.86-4.59)
|
|
Chi-square test, significant p-values (p < 0.05) are indicated in bold, trends towards significance (p < 0.08) in italics.
a Number of patients.
b HR: hazard ratio (CI: confidence interval) of univariate Cox regression analysis.
c Dichotomized into low and high levels by the 50th percentile for KLK8 and KLK11, by 67th percentile for KLK10.
d Dichotomized into low level by KLK8 low and KLK10 low, and high level by KLK8 high and/or KLK10 high.
e Dichotomized into low level by KLK10 low and KLK11 low, and high level by KLK10 high and/or KLK11 high.
f Dichotomized into low level by KLK8 low and KLK11 low, and high level by KLK8 high and/or KLK11 high.
g Dichotomized into low level by KLK8 low, KLK10 low as well as KLK11 low, and high level by KLK8 high, KLK10 high and/or KLK11 high.
Due to missing values, numbers do not always add up to n = 117 (DFS) and n = 121 (OS).
Elevated KLK8 mRNA expression levels were found to be a significant predictive factor for worse DFS (HR = 1.88, p = 0.037), but not for OS, indicating an about two-fold increased probability of tumor progression in the KLK8 high-expressing group. KLK10 mRNA levels also represented a significant predictive marker for shortened DFS (HR, 1.78; p = 0.045), but only showed a trend towards significance for OS (HR: 1.68; p = 0.076). Regarding KLK11, high expression levels showed a trend towards significance for DFS (HR: 1.85; p = 0.053) only, but were notably associated with worse OS (HR: 2.29; p = 0.016).
KLK8+KLK10 high expression, as well as KLK8+KLK11 high expression were also significantly associated with an elevated probability of disease recurrence (HR: 2.28, p = 0.013 and HR: 2.18, p = 0.042), while KLK10+KLK11 high expression was significantly not only correlated with shorter DFS (HR: 1.95, p = 0.046) but also with shorter OS (HR: 2.12, p = 0.035). Compared to the KLKs alone, the combined paired factors values, thus, could mildly increase the predictive power for DFS and/or OS. Finally, combination of all three biological factors (KLK8+KLK10+KLK11 low/low/low versus high and/or high and/or high) revealed that patients displaying tumor tissue-associated low expression of all three factors had a 2.85-fold reduced risk for disease recurrence compared to the group with high expression of at least one of these KLKs.
The impact of these factors on patient survival was also confirmed by Kaplan-Meier estimation. As shown by the respective survival curves, high KLK8 (Fig. 2A) and high KLK10 levels (Fig. 2B) were significantly associated with shortened DFS (p = 0.033, p = 0.042, respectively), whereas high KLK11 expression (Fig. 2C) was notably associated with both worse DFS (p = 0.049) and OS (p = 0.013). In case of the combined factors, all combinations represented prognostic factors for DFS (Fig. 3), whereas only the combination KLK10+KLK11 displayed a significant association with OS (see Table 2).
Association of KLK8, 10 and 11 mRNA expression levels with disease-free survival (DFS) and overall survival (OS) in multivariable Cox regression analysis
The independence of KLKs as prognostic factors in TNBC was studied in multivariable Cox regression analysis (Table 3). We established a base model containing age, lymph node status, and tumor size. Here, age was the only clinical parameter displaying predictive power for DFS (HR, 2.12; p = 0.025) and OS (HR, 3.18; p = 0.002). Among the tumor biological factors (added separately to the base model), KLK8 mRNA expression lost its prognostic significance for DFS, while KLK10 mRNA expression significantly contributed to the base model for DFS (HR: 2.19, p = 0.019). Elevated KLK11 mRNA expression turned out to represent an independent unfavorable predictor of OS (HR, 2.06; p = 0.044).
Table 3 Multivariate Cox regression analysis of clinical outcome in triple-negative breast cancer for tumor biological factors.
Clinicopathological parameters
|
DFS
|
OS
|
Noa
|
HR (95% CI)b
|
p
|
Noa
|
HR (95% CI)b
|
p
|
Age
|
|
|
0.025
|
|
|
0.002
|
≤ 60 years
|
54
|
1
|
|
54
|
1
|
|
> 60 years
|
41
|
2.12 (1.10-4.09)
|
|
43
|
3.18 (1.53-6.59)
|
|
Lymph node status
|
|
|
0.236
|
|
|
0.286
|
N0
|
51
|
1
|
|
53
|
1
|
|
N1/N2/N3
|
44
|
1.49 (0.77-2.87)
|
|
44
|
1.46 (0.73-2.94)
|
|
Tumor Size
|
|
|
0.305
|
|
|
0.423
|
≤20 mm
|
28
|
1
|
|
29
|
1
|
|
>20 mm
|
67
|
1.51 (0.69-3.33)
|
|
68
|
1.41 (0.61-3.29)
|
|
KLK8 mRNA c
|
|
|
0.062
|
|
|
0.390
|
low
|
46
|
1
|
|
48
|
1
|
|
high
|
49
|
1.90 (0.97-3.72)
|
|
49
|
1.36 (0.68-2.71)
|
|
KLK10 mRNA c
|
|
|
0.019
|
|
|
0.058
|
low
|
63
|
1
|
|
65
|
1
|
|
high
|
32
|
2.19 (1.14-4.20)
|
|
32
|
1.95 (0.98-3.91)
|
|
KLK11 mRNA c
|
|
|
0.113
|
|
|
0.044
|
low
|
51
|
1
|
|
53
|
1
|
|
high
|
44
|
1.70 (0.88-3.26)
|
|
44
|
2.06 (1.02-4.14)
|
|
KLK8+KLK10 c
|
|
|
0.013
|
|
|
0.163
|
low
|
39
|
1
|
|
41
|
1
|
|
high
|
56
|
2.62 (1.23-5.58)
|
|
56
|
1.70 (0.81-3.59)
|
|
KLK10+KLK11 c
|
|
|
0.049
|
|
|
0.054
|
low
|
45
|
1
|
|
47
|
1
|
|
high
|
50
|
1.97 (1.00-3.86)
|
|
50
|
2.03 (0.99-4.17)
|
|
KLK8+KLK11 c
|
|
|
0.048
|
|
|
0.171
|
low
|
31
|
1
|
|
33
|
1
|
|
high
|
64
|
2.15 (1.01-4.58)
|
|
64
|
1.72 (0.79-3.71)
|
|
KLK8+KLK0+KLK11 c
|
|
|
0.019
|
|
|
0.124
|
low
|
27
|
1
|
|
29
|
1
|
|
high
|
68
|
2.87 (1.19-6.91)
|
|
68
|
1.93 (0.83-4.48)
|
|
Chi-square test, significant p-values (p < 0.05) are indicated in bold, trends towards significance (p < 0.08) in italics.
a Number of patients.
b HR: hazard ratio (CI: confidence interval) of univariate Cox regression analysis.
c Dichotomization into a low versus high expressing group was performed as described in Table 2.
All analyzed combinations significantly contributed to the base model for DFS and thus represented independent predictive markers. Strikingly, both the low-groups of the combinations KLK8+KLK10 and KLK8+KLK10+KLK11, respectively, had an over 2.5-fold reduced risk for disease recurrence when compared to the high-groups (HR, 2.62; p = 0.013; HR: 2.87, p = 0.019; respectively).