Using the raw scores derived from analysis of the immunohistochemistry (Fig. 1), our data showed that the expression of ATM was negatively correlated with the expression of STAT1 (rs=-0.227, P=0.036), STAT3 (rs=-0.239, P=0.026), PD-L1 (rs=-0.213, P=0.049), TILs (rs=-0.279, P=0.009) and CD8+ cells (rs=-0.239, P=0.027), but no correlation was found with PD-1+ cells and the clinicopathological features (Table. 1). The positive or negative of “rs” means positive correlation or negative correlation respectively.
STAT1 but not STAT3 up-regulates expression of PD-L1, TILs and some specific markers in TNBC
Again, using the raw scores of immunohistochemistry in Fig. 1 (Table. 2), the data showed that the expression of STAT1 was positively correlated with the expression of PD-L1 (rs=0.422, P<0.001), TILs (rs=0.312, P=0.003), FOXP3+ (rs=0.312, P=0.004), CD8+ cells (rs=0.413, P<0.001) and STAT3 (rs=0.457, P<0.001). However, no correlation between STAT3 expression was found with that of PD-L1 (rs=0.023, P=0.832), TILs (rs=0.057, P=0.603), FOXP3+ (rs=0.058, P=0.595) and CD8+ cells (rs=0.103, P=0.345). The positive or negative of “rs” means positive correlation or negative correlation respectively.
Age, TILs, Ki67 index, TNM stage, p-c-srcY419, STAT1, STAT3, TILs, PD-1, FOXP3+ and CD8+ cells were included in multivariate analysis for PD-L1 status in TNBC. Univariate analysis showed that high STAT1 {OR (95% CI): 4.798 (1.824–12.616); P = 0.001}, TILs {OR (95% CI): 7.500 (2.313–24.320); P = 0.001} and CD8+ {OR (95% CI): 2.500 (1.024–6.106); P = 0.044} were associated with high PD-L1 level, while high ATM {OR (95% CI): 0.253 (0.101– 0.634); P = 0.003} was associated with lower PD-L1 level (Fig. 2). Multivariate analysis showed that high STAT1 {OR (95% CI): 5.895 (1.820–19.091); P = 0.003}, high TILs {OR (95% CI): 11.432 (2.862–45.662); P = 0.001} and higher TNM stage (Ⅲc VS Ⅲa+Ⅲb)} {OR (95% CI): 3.073 (1.025–9.216); P = 0.045, were independent predictors for high PD-L1 in TNBC (Fig. 2).
Only in tumors with ATM low expression that PD-L1 are independent negative and STAT3 are positive prognostic factors, respectively
Sixty-five of the 86 pTNM stage III TNBC cases have complete follow-up information from 12 to 63 months, with a median follow-up time of 44.5 months. There were 28 patients relapsed and 14 patients died (Supplementary Materials Table S2).
STAT3, but not PD-L1, had significant association with DFS (P=0.011 and P=0.122, respectively) and OS (P=0.054 and P=0.163, respectively) in all 65 patients (Supplementary Data Fig. S1). In TNBC with low ATM, increased PD-L1 expression significantly reduced DFS of the patients (P=0.026, Fig. 3a); however this impact of PD-L1 on DFS was not found in TNBC with high ATM expression (P=0.525, Fig. 3b). The increased PD-L1 expression had no significant impact on the OS of patients with either low (P=0.761, Fig. 3c) or high ATM expression (P=0.332, Fig. 3d). In the low ATM TNBC, positive expression of STAT3 was associated with prolonged DFS (P=0.005, Fig. 3e) and OS (P=0.021, Fig. 3g), However these impacts were not found in ATM high group (P=0.296, Fig. 3f for DFS and P=0.652, Fig. 3h for OS).
Univariate analysis showed that STAT3, p-c-srcY419 and PD-L1 were significantly associated with DFS of patients with low expression of ATM. However, in the high expression group, no factors analyzed in the study had significantly impact on prognosis (Table 3). Univariate analysis showed that only STAT3 affected OS in patient with ATM low expression, and again all parameters analyzed did not show significant association with OS in ATM high group. Multivariate analysis proved that expression of PD-L1 was an independent prognostic factor for reduced DFS in patients with TNBCs showing low ATM expression {HR (95% CI): 4.577 (1.250 – 16.764); P = 0.022, Table 3} but not in those with high ATM expression. STAT3 expression, but not high p-c-srcY419 expression, was an independent prognostic factor for prolonged DFS in TNBC with low ATM expression {HR (95% CI): 0.205 (0.067– 0.624); P = 0.005, Table 3}, but not in the high ATM expression group. STAT3 positive expression was the only independent prognostic factor for prolonged OS in patient with low ATM expression {HR (95% CI): 0.186 (0.037– 0.925); P = 0.040, Table 4}, but not in the high ATM expression group.
Table 1
Correlation between ATM and STAT1, STAT3, PDL-1, tumor infiltrating lymphocytes (TILs) and their some specific markers, as well as other pathological features in all 86 cases of pTNM stage Ⅲ TNBC.
|
|
ATM
|
Indicators
|
M(QR) f
|
rs
|
P
|
Age
|
53(14)
|
0.113
|
0.301
|
Ki67 index (%)
|
60(45)
|
-0.009
|
0.936
|
pTNM(Ⅲc VS Ⅲa+Ⅲb)
|
-
|
-0.103
|
0.344
|
STAT1a
|
20(100)
|
-0.227
|
0.036
|
STAT3 a
|
0(50)
|
-0.239
|
0.026
|
p-c-srcY419 (H-score) a
|
20(40)
|
0.006
|
0.956
|
PD-1b
|
-
|
-0.182
|
0.093
|
PDL-1c
|
0(1)
|
-0.213
|
0.049
|
TILs d
|
10(25)
|
-0.279
|
0.009
|
FOXP3+ e
|
10.5(15)
|
-0.076
|
0.487
|
CD8+ e
|
50(80)
|
-0.239
|
0.027
|
a: Caculated by H-score (Staining intensity×Percentage of positive cells) |
b: PD-1 staining was considered positive in cytomembrane of TILs, while the other cases of our study were considered negative. |
c: Caculated by IC score: the PD-L1 expression on stromal tumor-infiltrating immune cells occupying the proportion (%) of the tumor area. |
d: Scored using the 2014 international TILs working group standard |
e: Average number of positive cells in 10 fields at high magnification(×400) |
f: Expressed by median (interquartile range) |
Table 2
the relationship between STAT1 & STAT3 and PDL-1, tumor infiltrating lymphocytes and their some specific markers, as well as other pathological features in all 86 cases of pTNM stage Ⅲ TNBC.
|
|
STAT1
|
STAT3
|
Indicators
|
M(QR) f
|
rs
|
P
|
rs
|
P
|
Age
|
53(14)
|
-0.093
|
0.393
|
-0.039
|
0.721
|
Ki67 index (%)
|
60(45)
|
0.125
|
0.252
|
-0.019
|
0.865
|
pTNM(Ⅲc VS Ⅲa+Ⅲb)
|
-
|
-0.077
|
0.481
|
-0.185
|
0.088
|
STAT1a
|
20(100)
|
-
|
-
|
0.457
|
<0.001
|
STAT3a
|
0(50)
|
0.457
|
<0.001
|
-
|
-
|
p-c-srcY419 (H-score) a
|
20(40)
|
0.239
|
0.027
|
0.234
|
0.030
|
PD-1b
|
-
|
0.189
|
0.081
|
0.121
|
0.268
|
PDL-1c
|
0(1)
|
0.422
|
<0.001
|
0.023
|
0.832
|
TILs d
|
10(25)
|
0.312
|
0.003
|
0.057
|
0.603
|
FOXP3+ e
|
10.5(15)
|
0.312
|
0.004
|
0.058
|
0.595
|
CD8+ e
|
50(80)
|
0.413
|
<0.001
|
0.103
|
0.345
|
a: Caculated by H-score (Staining intensity×Percentage of positive cells) |
b: PD-1 staining was considered positive in cytomembrane of TILs, while the other cases of our study were considered negative. |
c: Caculated by IC score: the PD-L1 expression on stromal tumor-infiltrating immune cells occupying the proportion (%) of the tumor area. |
d: Scored using the 2014 international TILs working group standard |
e: Average number of positive cells in 10 fields at high magnification(×400) |
f: Expressed by median (interquartile range) |
Table 3
The effect of STAT1, STAT3, PDL-1, TILs and their some specific markers as well as other clinicopathological features, on the invasive disease-free survival (iDFS) of ATM-low TNBC patients analyzed by univariate and multivariate analysis
Indicators
|
Univariate analysis
|
Multivariate Analysis
|
HR
|
95%CI
|
P
|
HR
|
95%CI
|
P
|
Age {≥53(median) VS <53}
|
0.923
|
0.334-2.550
|
0.877
|
|
|
|
p-c-srcY419 (H score ≥20 VS <20)
|
0.307
|
0.097-0.974
|
0.045
|
|
|
|
STAT1 (H score ≥20 VS <20)
|
0.739
|
0.268-2.041
|
0.560
|
|
|
|
STAT3 (Positive VS Negative)
|
0.244
|
0.082-0.722
|
0.011
|
0.205
|
0.067-0.624
|
0.005
|
Ki67 index (≥60% VS < 60%)
|
1.476
|
0.524-4.156
|
0.461
|
|
|
|
pTNM(Ⅲc VS Ⅲa+Ⅲb)
|
1.442
|
0.521-3.988
|
0.481
|
|
|
|
TILs (≥10% VS < 10%)
|
2.071
|
0.466-9.199
|
0.339
|
|
|
|
PD-1(Positive VS Negative)
|
1.007
|
0.343-2.953
|
0.990
|
|
|
|
PDL-1(Positive VS Negative)
|
3.769
|
1.057-13.441
|
0.041
|
4.577
|
1.250-16.764
|
0.022
|
FOXP3+ a(≥10.5 VS < 10.5)
|
1.227
|
0.344-3.388
|
0.693
|
|
|
|
CD8+ a(≥50 VS < 50)
|
0.864
|
0.307-2.429
|
0.781
|
|
|
|
a: Average number of positive cells in 10 fields at high magnification(×400) |
Table 4
The effect of STAT1, STAT3, PDL-1, TILs and their some specific markers as well as other clinicopathological features, on the overall survival (OS) of ATM-low TNBC patients analyzed by univariate and multivariate analysis
Indicators
|
Univariate analysis
|
Multivariate Analysis
|
HR
|
95%CI
|
P
|
HR
|
95%CI
|
P
|
Age {≥53(median) VS <53}
|
0.729
|
0.182-2.916
|
0.655
|
|
|
|
p-c-srcY419 (H score ≥20 VS <20)
|
0.542
|
0.129-2.272
|
0.402
|
|
|
|
STAT1 (H score ≥20 VS <20)
|
0.350
|
0.083-1.468
|
0.151
|
|
|
|
STAT3 (Positive VS Negative)
|
0.186
|
0.037-0.925
|
0.040
|
0.186
|
0.037-0.925
|
0.040
|
Ki67 index (≥60% VS < 60%)
|
3.949
|
0.791-19.726
|
0.094
|
|
|
|
pTNM(Ⅲc VS Ⅲa+Ⅲb)
|
1.334
|
0.332-5.354
|
0.685
|
|
|
|
TILs (≥10% VS < 10%)
|
2.376
|
0.292-19.343
|
0.419
|
|
|
|
PD-1(Positive VS Negative)
|
0.891
|
0.212-3.739
|
0.874
|
|
|
|
PDL-1(Positive VS Negative)
|
1.248
|
0.298-5.229
|
0.762
|
|
|
|
FOXP3+ a(≥10.5 VS < 10.5)
|
1.813
|
0.432-7.608
|
0.416
|
|
|
|
CD8+ a(≥50 VS < 50)
|
0.433
|
0.108-1.735
|
0.237
|
|
|
|
a: Average number of positive cells in 10 fields at high magnification(×400) |
Low expression of ATM and phosphorylation of tyrosine at position 419 of c-Src induces overexpression of STAT3, but not STAT1, and improves PFS of patients. To explore the regulatory mechanism of STAT3 and PD-L1, and to analyze why STAT3 and PD-L1 only impact on prognosis in TNBC with low expression of ATM, we analyzed the relationship between p-c-src Y419 and STAT1, STAT3, PD-L1, TILs and some specific tumor immune cell markers at different ATM expression levels. We found in TNBC with low expression of ATM, p-c-src Y419 was only positively associated with STAT3 (rs=0.356, P=0.022, Table 5), but not significantly with other factors including STAT1 (rs=0.064, P=0.693, Table 5). However in tumors with high expression of ATM, p-c-src Y419 was positive correlated with STAT1 (rs=0.426, P=0.004, Table 5) and TILs (rs=0.335, P=0.025, Table 5), but not significantly with other factors including STAT3 (rs=0.141, P=0.357, Table 5)
Table 5
The relationship between p-c-src Y419 and STAT1, STAT3, PDL-1, TILs and some specific tumor immune cell markers at different ATM expression levels of TNM stage Ⅲ TNBC
|
p-c-srcY419 a
|
|
ATM low (N=41)
|
ATM high (N=45)
|
Indicators
|
M(QR) f
|
rs
|
P
|
M(QR) f
|
rs
|
P
|
Ki67 index (%)
|
60(32.5)
|
0.012
|
0.941
|
60(50)
|
-0.133
|
0.383
|
STAT1a
|
30(110)
|
0.064
|
0.693
|
10(67.5)
|
0.426
|
0.004
|
STAT3a
|
5(70)
|
0.356
|
0.022
|
0(30)
|
0.141
|
0.357
|
PDL-1c
|
1(1)
|
-0.261
|
0.099
|
0(5)
|
0.056
|
0.716
|
TILs d
|
20(20)
|
0.145
|
0.365
|
10(15)
|
0.335
|
0.025
|
FOXP3+ e
|
11(13.5)
|
0.232
|
0.145
|
10(15)
|
0.087
|
0.568
|
CD8+ e
|
64(119.5)
|
0.206
|
0.196
|
40(50.5)
|
0.266
|
0.077
|
a: Caculated by H-score (Staining intensity×Percentage of positive cells) |
b: PD-1 staining was considered positive in cytomembrane of TILs, while the other cases of our study were considered negative. |
c: Caculated by IC score: the PD-L1 expression on stromal tumor-infiltrating immune cells occupying the proportion (%) of the tumor area. |
d: Scored using the 2014 international TILs working group standard |
e: Average number of positive cells in 10 fields at high magnification(×400) |
f: Expressed by median (interquartile range) |
In TNBC with low ATM expression, high expression of p-c-srcY419 was significantly associated with a better DFS than that of low expression (P=0.032, Fig. 4a), although it is no proven to be independent at multivariate analysis (Table 3) and did not show value in the impact on OS (P=0.395, Fig. 4c). In ATM high expression group, p-c-srcY419 expression was not associated with improved DFS or OS at univariate (Fig. 4b & 4d) and multivariate analysis (Table 3 & 4).