Cohort description and patients’ distribution into subcohorts
In total, data and samples from 131 patients with early ER+/HER2- BC were analysed. The main characteristics of the patients and tumours, their surgical management and the pathological changes after NET are summarised in Table 1. The study population presented a mean age at diagnosis of 70 (47-93) and the mean NET duration before surgery was 9 months (2-40). Letrozole was the endocrine treatment administered in most cases (94.7%) and breast-conserving surgery was performed in 79% of the patients. Of note, one patient was diagnosed with bilateral disease, and her two tumours were independently considered in our analyses. Moreover, two patients achieved a pCR after NET (1.5% of patients) and, consequently, their post-NET surgical samples were not available for biological evaluation. Apart from these mentioned exceptions, matched samples from each patient were analysed (a core diagnostic biopsy or pre-NET sample, and the surgical specimen or post-NET sample). Suppl Mat Fig S1 shows the distribution of patients in 3 subcohorts and the different analyses performed in each of them. The distribution of patients into each subcohort took into account the sample availability (e.g. location) and sample quality required for each analysis. Our subcohorts are representative of the general cohort as observed in Suppl Mat Tab 1.
PAM50 analysis (intrinsic subtype and ROR score) in core biopsies can predict NET response
The 50 gene expression-based PAM50 assay allows the classification of BC into 5 intrinsic biological subtypes and generates risk of recurrence (ROR) scores for each sample [39]. However, the utility of PAM50 to predict treatment benefit requires further investigation in the specific context of NET. In this regard, we selected a representative subcohort of ER+/HER2- BC patients treated with NET (subcohort #1 in Suppl Mat Fig S1) and we obtained the PAM50 data (intrinsic subtype and ROR scores) on their diagnostic biopsies (pre-NET) and surgery specimens (post-NET) (Table 2). As expected, due to the characteristics of our cohort, low number of patients presented HER2-E or basal subtypes.
First, to check if our subcohort behaved similarly to others of this kind, we validated previous findings regarding the predictive value of Luminal A (LumA) and Luminal B (LumB) baseline PAM50 categorization in ER+ BC treated with NET [33]. As expected, our results showed that LumA presented lower Ki67 after treatment, ΔKi67 and mPEPI score than LumB, indicating a better biological response to NET and better prognosis (Figure 1A-B, Suppl Mat Fig S2A). Moreover, to add further evidence of the predictive value of PAM50 evaluation at baseline (pre-NET), we analysed ROR score data. We observed that tumours with high ROR scores had worse response to NET evaluated as higher Ki67 % at surgery and a lower change in Ki67 (measured by ΔKi67) together with a worse prognosis with a higher mPEPI score (Figure 1C-D, Suppl Mat Fig S2B and Table 2). Moreover, tumour cellularity size (TCS), a novel biomarker for NET response described for the first time by Lopez-Velazco et al. [25] correlates with ROR-S at baseline (pre-NET) (Figure 1E). No differences were found for ROR-S and ROR-P (Suppl Mat Fig S2C-F). These results support that PAM50 analysis in the diagnostic biopsy could help to personalize the use of NET, conjointly with other clinical parameters.
Conversion from Luminal A to Normal-like intrinsic subtype after NET is associated with better response
As mentioned above, NET is a therapeutic strategy that allows the characterisation of tumour response by comparing diagnostic biopsies (pre-NET) and surgical specimens (post-NET) [17,31]. Nowadays, limited data exist regarding intrinsic subtype changes in paired pre- and post- NET samples of ER+/HER2- BC patients. That is why we characterised the changes in the intrinsic subtype upon NET in matched samples and analysed their relationship to treatment response.
The distribution of PAM50-based intrinsic subtypes among pre- and post-NET samples and the value of the main NET biomarkers of each change group are described in Table 2, Figure 2 and Suppl Mat Fig S3A-C. Our results showed that the two most represented profiles of intrinsic subtype after NET were: the group of tumours with a persistent LumA status (36% of our patients) and the group of tumours with a conversion of LumA to Normal-like status upon therapy (26% of patients). We wondered whether there could be a difference in response/prognosis to NET between these two groups which share the LumA pre-NET status. Our results showed that tumours changing from LumA to Normal-like after NET presented: lower Ki67 levels at surgery, a significant decrease in Ki67 (ΔKi67), lower mPEPI score and lower TCS (Figure 2C-E, Suppl Mat Fig S3D). Importantly, Lopez-Velazco et al. [25] proposed a cutoff for TCS at 2.5mm. We found that 92% of Persistent LumA tumours presented a high TCS (>2.5mm), while LumA to normal subgroup was enriched in low TCS (≤ 2.5 mm) tumours (Figure 2F). Later, we assessed whether these results could be conditioned by the tumour size or epithelial cellularity content. Our analyses showed that our two populations of interest did not statistically differ in these parameters although a tendency is observed regarding to cellularity (Figure 2G-H).
Next, we analysed if the change in ROR may also help us to evaluate NET response and prognosis. We observed that, in our cohort, ROR values are significantly lower after NET compared to the ones obtained pre-NET (Table 2, Figure 3A, Suppl Mat Fig S3E). Interestingly, Figure 3A shows that none of tumours changing from Luminal A (LumA) to normal-like (normal) intrinsic subtype increases its ROR after NET. The potential of ROR change (∆ROR) between pre- and post-NET samples as a marker of NET response and prognosis has not been largely explored in previous studies. Our results indicate that ∆ROR value should be taken into consideration, at least in the case of ROR-S. First, we plotted ∆ROR-S values from all patients (Figure 3B) and found 3 subpopulations: i) ∆ROR-S ≥ 0 (including patients with no change or increase of ROR-S after NET; 14 out of 59); ii) 0>ΔROR-S>-2.5 (including patients with a slight decrease in ROR after NET; 31 out of 59); and iii) ΔROR-S≤-2.5 (including patients with an accentuated change in ROR after NET; 14 out of 59). ΔROR-S values are significantly different between ii) and iii) according to statistical analyses, indicating that a cut-off at ΔROR-S = -2.5 could stratify two subpopulations of patients. At this point, we wondered whether these two subpopulations might show differences in their NET response and prognosis. Statistically significant differences in Ki67, ΔKi67 and mPEPI between these novel subgroups of patients subdivided taking into account their ΔROR-S values (Figure 3C-E). Interestingly, subgroups ii) and iii) present clear differences, being subgroup ii) more similar to i) than to iii). TCS results showed a trend to lower TCS in the iii) group, although no statistically significant was reached (Figure 4F) likely due to insufficient number of patients per group. Thus, this novel stratification of patients into 3 subgroups according to their changes in ROR-S (ΔROR-S, establishing cut-off at ΔROR-S=0 and ΔROR-S=-2.5) has revealed relevant differences in NET response and prognosis between patients depending on the magnitude of their ROR decrease after NET.
In summary, our data indicate that analysing the change in the intrinsic subtype or ROR value upon NET can provide information about response to NET and patient prognosis that could help in the clinical practice to, for example, determine adjuvant therapy.
Quantification of p53 levels before and after NET may contribute to asess NET response
Previous studies indicate that biological markers (e.g. senescence or apoptosis markers) may be useful to determine the efficacy of a neoadjuvant therapy (mainly, chemotherapy) in BC patients [50–53]. However, few data are published in the NET setting in ER+/HER2- BC. Here, we study the protein levels of different biomarkers of apoptosis (Bcl2), cancer stem cell phenotype (Sox2) and senescence (p53, p21 and p16) by IHC in matched pre- and post-NET samples as well as the value of their change after NET. For that, we selected a novel subcohort (subcohort #2 in Suppl Mat Fig S1), representative of our global cohort of ER+/HER2- BC patients treated with NET, composed of 61 patients.
The quantifications of the mentioned markers in each sample type are shown in Figure 4A and Suppl Mat Fig S4A-D. The percentage of positive cells for p53 and p21, two senescence markers, showed a clear decrease after NET (FC≥-2), while the percentage of positive cells for Bcl2 (apoptosis marker) presented also a significant decrease after NET but the magnitude change was lower (FC=-1.1). To determine the predictive/prognostic value of these markers, we compared them (pre-NET, post-NET or changes) with the validated predictive NET biomarkers (Ki67 expression at surgery, ∆Ki67 and with the prognostic mPEPI score). These analyses were performed with all the markers, although Sox2, Bcl2 and p16 did not show any association (data not shown). First, our results showed that, in the pre-NET samples, only p53 might have predictive capacity since a higher expression was associated with higher Ki67 and mPEPI grading, indicating poor response and prognosis after NET (Figure 4B-C). Next, we studied the protein levels of the mentioned markers in our post-NET samples. We observed that the expression of p53 and p21 was inversely associated with response to NET (Figure 4D-E and Suppl Mat Fig 5A-B). Finally, we decided to study if the change in the mentioned markers before and after NET may be related to NET response. Importantly, we found that the change in levels of p53 (∆p53) after NET may be associated to response/prognosis due to their correlation with mPEPI score and TCS (cutoff at 2.5mm) (Figure 5F-G). The change in levels of p21 (∆p21) after NET may be also associated to NET response according to its correlation with Ki67 (Suppl Mat Fig S5C). However, further validation is needed to establish their value as long-term prognostic biomarkers.
In summary, in the search for new biomarkers of response with long-term prognostic value upon NET, our results suggest that, among another biological biomarkers evaluated, the percentage of positive cells for p53 (assessed in pre-NET and post-NET samples and its change) showed promising results in terms of its ability to characterise residual cell populations with prognostic interest, deserving this further validation.
A change in HER2 status is observed after NET, although it is not associated to tumour response
HER2-low status has gained huge attention during the last years in the field due to the fact that some studies point out that it should be considered as an independent biologic subtype distinct from HER2-0 BC, although its clinic implications are not well-elucidated [49,54]. Here, we explored the relationship of HER2-low status with NET response and prognosis. For that, we defined a novel subcohort (subcohort #3 in Suppl Mat Fig S1), representative of our global cohort of ER+/HER2- BC patients treated with NET, composed of 99 patients. We collected retrospectively the patients’ data related to HER2 in their diagnostic biopsies (pre-NET) and surgical specimens (post-NET) (Table 3 and Figure 5A-B).
First, we studied if subdividing pre-NET samples into HER2-0 and HER2-low may predict NET response. We observed that, at diagnosis, pathologists’ analyses indicated that most of our patients could be classified as HER2-low (HER2-low: n=66; HER2-0: n=33), but we did not find any relationship between HER2 status and NET predictive/prognostic markers (Ki67, mPEPI and ∆Ki67; Suppl Mat Fig 6A-C). Similarly, when we analysed if the stratification of post-NET samples into HER2-0 and HER2-low groups may be related to Ki67 levels and PEPI score, we did not obtain any positive result (HER2-low: n=32; HER2-0: n=62) (Suppl Mat Fig 6D-F).
Importantly, by comparing HER2 status between pre- and post-NET samples (Figure 5A-B), we observed a statistically significant decrease in HER2 scoring after NET and, consequently, enrichment in HER2-0 samples (Suppl Mat Fig 7A). Around half of the samples presented a discordant pre- vs post-NET HER2 status (40 of 94 patients (43%) when comparing HER2-0 vs HER2-low; Table 3). In detail, 36 patients (38%) changed their HER2 status from HER2-low (pre-NET) to HER2-0 (post NET) (Suppl Mat Fig 7A and Table 3), being this the most prevalent change observed in our cohort. Thus, we checked whether this change may be used as a marker of response to NET (compared to HER2-0 and HER2-low persistent patients), but we did not observe any relationship with the aforementioned NET prognostic/predictive markers (Suppl Mat Fig 7B-D). We explored other analyses, comparing different groups, but we did not obtain any positive relationship between HER2 status change after NET and NET response (data not shown).
Taken together, the results of our ER+/HER2- cohort indicate that, after NET, a generalized decrease in HER2 status can be observed but it is not associated to NET response. Analysing pre-NET and post-NET data individually neither supports the interpretation of HER2-low status as a distinct biologic entity in the context of NET response in ER+/HER2- patients.