Baseline characteristics according to HER2 expression
The majority of 2,099 HR positive cases were HER2 low expressed tumor (n = 1,732, 82.5%). No significant difference was found when comparing the age, menopause status, histological grade, progesterone receptor (PR status, American Joint Committee on Cancer (AJCC) stage, RS range and treatment after surgery between HER2-low and HER2 negative subgroups (Table 1). The HER2 negative subgroup had a higher proportion of specific pathologic types beyond invasive ductal and lobular cancer (16.1% vs 10.5%, p = 0.01) and tumors with high histologic grade (p = 0.02). We also divided ER IHC expression into two categories based on the median percentage of 90%. The HER2 low tumors showed a higher ratio of concurrent ER high expression than HER2 negative patients (78.2% vs 58.6%, p < 0.01).
Table 1
Basic Characteristics of 2,099 HR + breast cancer patients from SJTU-BCDB
Characteristics
|
HER2 negative
n = 367
|
HER2 low
n = 1,732
|
p-value
|
Age
≤50
>50
|
110(30.0)
257(70.0)
|
562(32.4)
1,170(67.6)
|
0.39
|
Menopause status
Premenopausal
Postmenopausal
|
123(33.5)
244(66.5)
|
607(35.0)
1,125(65.0)
|
0.61
|
Pathology
IDC
ILC
others
|
294(80.1)
14(3.8)
59(16.1)
|
1,478(85.3)
72(4.2)
182(10.5)
|
0.01
|
Grade
I
II
III
Unknown
|
36(9.8)
173(47.1)
77(21.0)
81(22.1)
|
147(8.5)
1,027(59.3)
312(18.0)
246(14.2)
|
0.02
|
ER positivity(%)
<90
≥90
|
152(41.4)
215(58.6)
|
377(21.8)
1355(78.2)
|
< 0.01
|
PR
Positive
Negative
|
313(85.3)
54(14.7)
|
1,540(88.9)
192(11.1)
|
0.06
|
Ki67(%)
<20
≥20
|
213(58.0)
154(42.0)
|
951(54.9)
781(45.1)
|
0.30
|
pT
1
2–4
|
262(71.4)
105(28.6)
|
1,215(70.2)
517(29.8)
|
0.68
|
pN
0
1–2
|
313(85.3)
54(14.7)
|
1,448(83.6)
284(16.4)
|
0.47
|
RS
<18
18–30
>30
|
83(22.6)
185(50.4)
99(27.0)
|
361(20.8)
914(52.8)
457(26.4)
|
0.67
|
Endocrine therapy
AI(± OFS)
SERM(± OFS)
others
|
246(67.0)
120(32.7)
1(0.3)
|
1,114(64.3)
616(35.6)
2(0.1)
|
0.29
|
Chemotherapy
No
Yes
|
196(53.4)
170(46.6)
|
859(49.6)
872(50.4)
|
0.19
|
Radiotherapy
No
Yes
|
197(53.7)
169(46.3)
|
903(52.2)
828(47.8)
|
0.60
|
Abbreviation: IDC, invasive ductal carcinoma; ILC, invasive lobular carcinoma; ER, estrogen receptor; PR, progesterone receptor; RS, recurrence score; AI, aromatase inhibitor; SERM, selective estrogen receptor modulator; OFS, ovarian function suppression. |
The survival outcomes of patients according to different HER2 status
A total of 146 DFS events happened during the follow-up period (33 and 123 in HER2 negative and low subgroups respectively). The Kaplan-Meier curves (Fig. 1) revealed no significant difference between HER2 negative and low subgroups in HR + patients. We further performed survival analyses according to different RS, and the subpopulation treatment effect pattern plot (STEPP) analysis22 (Fig. 3) showed that in patients with RS > 30, HER2 low subgroup showed a better DFS than HER2 negative subgroup. We compared the predictive value of RS in HR + patients with different HER2 status (Fig. 2). Surprisingly, in HER2 negative patients, the Kaplan-Meier curves showed significant different survival outcomes between patients with high risk and those with low/intermediate risk (p = 0.003). However, the recurrence score was not significantly correlated with survival outcomes in HER2 low patients (p = 0.11).
Univariate And Multivariate Analysis Of Clinicopathological Predictors For Dfs
We further investigated the impact of clinicopathological factors on DFS in subgroups divided by HER2 status respectively (Table 2). In HER2 negative population, higher pN grade (hazard ratio = 2.59, 95% CI 1.08–6.26, p = 0.03) and RS > 30 (hazard ratio = 2.39, 95% CI 1.14–5.04, p = 0.02) was associated with worse survival outcomes after adjustment. However, in patients with HER2 low tumors, higher RS and pN stage lost their prognostic value both in univariate and multivariate analysis.
Table 2
The impact of clinicopathological factors on DFS according to HER2 status
Factors
|
Univariate Analysis
|
p-value
|
Multivariate Analysis
|
p-value
|
HER2 negative
|
Age(y)
>50 vs ≤ 50
|
1.08(0.51–2.66)
|
0.85
|
|
|
Menopause
Post vs Pre
|
1.37(0.65–2.88)
|
0.42
|
|
|
Grade
III vs I-II
|
1.67(0.76–3.65)
|
0.20
|
|
|
ER positivity(%)
≥90 vs < 90
|
0.46(0.21–1.02)
|
0.06
|
0.53(0.23–1.20)
|
0.13
|
PR status
Pos vs Neg
|
0.39(0.18–0.81)
|
0.01
|
0.62(0.28–1.39)
|
0.24
|
Ki67 positivity(%)
≥20 vs < 20
|
1.68(0.82–3.43)
|
0.16
|
|
|
pT
2–4 vs 1
|
1.99(0.98–4.01)
|
< 0.01
|
1.51(0.70–3.24)
|
0.29
|
pN
≥ 1 vs 0
|
2.84(1.24–6.49)
|
0.01
|
2.59(1.08–6.26)
|
0.03
|
RS
>30 vs ≤ 30
|
3.10(1.56–6.17)
|
< 0.01
|
2.39(1.14–5.04)
|
0.02
|
Chemotherapy
No vs Yes
|
1.44(0.72–2.86)
|
0.30
|
|
|
HER2 low
|
Age(y)
>50 vs ≤ 50
|
0.90(0.62–1.31)
|
0.60
|
|
|
Menopause
Post vs Pre
|
0.91(0.63–1.30)
|
0.60
|
|
|
Grade
III vs I-II
|
1.73(1.15–0.62)
|
0.01
|
1.30(0.82–2.05)
|
0.27
|
ER positivity(%)
≥90 vs < 90
|
0.66(0.45–0.98)
|
0.04
|
0.52(0.23–1.19)
|
0.12
|
PR status
Pos vs Neg
|
0.76(0.47–1.24)
|
0.27
|
|
|
Ki67 positivity(%)
≥20 vs < 20
|
1.71(1.20–2.45)
|
< 0.01
|
1.47(0.93–2.34)
|
0.10
|
pT
2–4 vs 1
|
2.19(1.54–3.13)
|
< 0.01
|
1.39(0.91–2.11)
|
0.13
|
pN
≥ 1 vs 0
|
1.15(0.70–1.91)
|
0.58
|
1.20 (0.68–2.10)
|
0.53
|
RS
> 30 vs ≤ 30
|
1.47(1.02–2.12)
|
0.04
|
0.86(0.49–1.52)
|
0.60
|
Chemotherapy
No vs Yes
|
1.44(1.00-2.06)
|
0.05
|
1.13(0.70–1.81)
|
0.62
|
Abbreviation: ER, estrogen receptor; PR, progesterone receptor; RS, recurrence score; HER2, human epidermal growth factor receptor 2. |
The distribution of genetic risks and gene expression between HER2 negative and HER2-low subgroups
We compared the expression of individual RS genes between HER2 low and negative groups. The results showed no significant difference of gene expression except for HER2 gene. The level of HER2 gene mRNA was slightly higher in HER2 low subgroup than in HER2 negative one (Fig. 2).
The impact of RS modules on DFS according to HER2 status
We further evaluated the predictive roles of separate RS modules and genes in patients according to HER2 status and observed different patterns (Table 4).
Table 4 The impact of RS modules and genes on DFS according to HER2 status
|
|
Factors
|
HER2-negative
|
p-value
|
HER2-low
|
p-value
|
pinteraction
|
ER Module
ER
PgR
Bcl2
SCUBE2
ER Module
|
1.20(0.88-1.65)
0.84(0.67-1.06)
0.88(0.58-1.32)
0.89(0.70-1.13)
0.84(0.58-1.21)
|
0.25
0.15
0.53
0.34
0.34
|
1.01(0.85-1.19)
0.97(0.87-1.08)
0.87(0.71-1.06)
0.97(0.86-1.09)
0.92(0.76-1.11)
|
0.95
0.55
0.16
0.56
0.38
|
0.36
0.30
0.67
0.39
0.52
|
Proliferation Module
|
|
Ki67
STK15
Survivin
CCNB1
MYBL2
Prol Module
|
2.20(1.46-3.30)
1.85(1.26-2.72)
1.72(1.25-2.37)
1.48(0.90-2.44)
1.52(1.15-2.02)
7.22(2.49-20.93)
|
<0.01
<0.01
<0.01
0.12
<0.01
<0.01
|
1.02(0.88-1.19)
1.02(0.89-1.17)
1.06(0.89-1.26)
1.19(0.95-1.50)
1.08(0.90-1.28)
0.81(0.40-1.66)
|
0.80
0.78
0.53
0.14
0.41
0.57
|
<0.01
<0.01
0.03
0.64
0.12
<0.01
|
HER2 Module
|
|
|
|
|
|
GRB7
HER2
HER2 Module
|
1.39(1.04-1.87)
1.02(0.70-1.49)
1.32(0.71-2.48)
|
0.03
0.91
0.38
|
1.01(0.82-1.24)
0.99(0.81-1.23)
1.89(1.02-3.48)
|
0.92
0.98
0.04
|
0.27
0.76
0.32
|
Invasion Module
|
MMP11
CTSL2
Inv Module
|
1.11(0.83-1.48)
1.39(1.06-1.83)
1.58(1.04-2.38)
|
0.49
0.02
0.03
|
1.20(1.04-1.39)
1.08(0.90-1.29)
1.30(1.05-1.61)
|
0.02
0.42
0.02
|
0.53
0.30
0.68
|
GSTM1
|
1.15(0.82-1.62)
|
0.43
|
0.85(0.70-1.02)
|
0.08
|
0.32
|
CD68
|
1.57(1.02-2.41)
|
0.04
|
1.02(0.80-1.30)
|
0.87
|
0.20
|
BAG1
|
1.13(0.73-1.77)
|
0.58
|
1.01(0.82-1.24)
|
0.96
|
0.99
|
In HER2 negative subgroup, most of the proliferation related genes played important parts in predicting DFS, as well as the whole module. Although GRB7 in HER2 module acted as a possible predictor of poorer survival (hazard ratio = 1.39, 95%CI 1.04–1.87, p = 0.03), HER2 mRNA expression and HER2 module didn’t had significant prognostic value. The increased invasion module score was also associated with worse survival outcome. No significance was association was found between ER module and DFS.
In HER2 low subgroup, the proliferation related genes lost prognostic value. In addition, unlike in the HER2 negative subgroup, the higher score of HER2 module might predict increased risk of DFS events (hazard ratio = 1.89, 95%CI 1.02–3.48, p = 0.04). The invasion module still served as a negative prognostic factor, while the correlation between ER related genes and DFS was not obvious in HER2 low as in HER2 negative patients.
The interaction test for HER2 status and proliferation module showed that there were significant interactive effects of most proliferation related genes and HER2 status on DFS.