Evaluating the expression of heat shock protein 27 and topoisomerase II α in a retrospective cohort of patients diagnosed with locally advanced breast cancer and treated with neoadjuvant anthracycline-based chemotherapies

Background: Neoadjuvant anthracycline-based chemotherapy (NAC) is a major regimen for the treatment of local advanced breast cancer (LABC), while resistance to NAC remains a paramount clinical obstacle. To investigate the role of heat shock protein 27 (Hsp27) and/or topoisomerase IIα (TopoIIα) in LABC patients treated with NAC, we performed this retrospective study.

Methods: Associations of Hsp27 transcripts with clinic-pathological characteristics, survival and drug response were investigated in public databases. Hsp27-related genes were identified, followed by functional enrichment analyses. Besides, two protein-protein interaction networks were built. Then, tumors from 103 patients who were diagnosed with LABC and received NAC were collected, and Hsp27 and TopoIIα were examined by Immunohistochemistry (IHC). Chi-square or Fisher’s exact tests were performed, as well as survival analyses.

Results: Either at the transcriptional level in public databases or at the protein level tested by IHC, a high level of Hsp27 was associated with aggressive tumor characteristics such as lymph node invasion and chemotherapy resistance. Hsp27-related genes mostly involved in the metabolic pathway and the gamete generation biological process. An elevated Hsp27 indicated a poor prognosis in patients with breast cancer (log-rank test P = 0.002 and 0.004 for disease-free survival [DFS] and overall survival [OS], respectively), while it might not be an independent predictor. Of note, tumors with high TopoIIα expression (TopoIIα+) was less likely to express Hsp27 (Hsp27+), in contrast to those with TopoIIα negativity (31.1% vs. 86.2%, P<0.001), and survival analyses revealed that patients with Hsp27+ and TopoIIα- tumors had a significantly lower DFS and OS (log-rank test P < 0.001 and 0.001, respectively), in contrast to the other three groups.

Conclusions: Hsp27 was associated with aggressive breast cancers and more predictable for the prognosis of LABC patients treated with NAC when concomitantly considering TopoIIα expression.

Introduction

Locally advanced breast cancer (LABC) composed around 30% to 40% of all breast cancers in developing countries (1). To reduce the tumor burden for later surgical adjuvant therapy, neoadjuvant chemotherapy (NAC) is commonly used for LABC, while multidrug resistance weakens the therapeutic efficacy and remains the most critical challenge (2, 3) Therefore, identifying patients who may respond to specific treatment regimens before or during the early stage of therapy is essential for avoiding ineffective treatments and achieving better outcome.

Heat shock protein 27 (Hsp27) was reported to overexpress in various cancers, and its overexpression is associated with tumorigenesis, metastasis, and invasiveness (4). Of note, Hsp27 is known for being involved in chemotherapy resistance; related medications include doxorubicin, gemcitabine, erlotinib, 5-fluorouracil, and temozolomide, etc (5). As for breast cancer, Hsp27 was demonstrated to be overexpressed in the aggressive forms, e.g., higher levels in tumors with lymph node metastasis compared to tumors without lymph node metastasis (6), while its value as a prognostic indicator remains arguable (7–9). One early study indicated nuclear Hsp27 was increased after induction chemotherapy (10), that was confirmed later in another study (11). However, this study failed to detect the prognostic value of Hsp27 in LABC patients treated with doxorubicin- or epirubicin-based monochemotherapy (11). Considered of rare evidence specific to LABC and critical choice of NAC to later surgery and outcomes (12), it is worth to further clarity its role in NAC with and without other biomarkers.

DNA topoisomerase IIα (TopoIIα) is one molecular target for the anthracycline class of chemotherapeutic drugs. Overexpression of TopoIIα predicts the responsiveness of breast cancer to anthracycline-based therapy (13). Ectopic overexpression of Hsp27 reduces the expression of TopoIIα in breast cancer cells (14). However, the relationship between Hsp27 and TopoIIα in breast cancer, especially the prognostic value of the combination of Hsp27 and TopoIIα in LABC has not been defined.

In this study, we first worked on public databases and investigated the relation between Hsp27 transcripts and clinic-pathological features, prognosis and chemotherapy resistance. Hsp27-related genes were identified followed by functional enrichment analyses, and two protein-protein interaction networks were built. Then, we detected the expression of Hsp27 and TopoIIα in pretreatment biopsies of LABC patients treated with NAC, and evaluated the association between the expression of Hsp27 and the efficacy of neoadjuvant chemotherapy, and the other clinicopathological parameters. Moreover, we assessed the association between Hsp27 expression and prognosis of LABC, and verified the predictive value of Hsp27 and TopoIIα in combination for the prognostication of LABC patients.

Materials and methods

Public databases

Associations of Hsp27 and clinic-pathological characteristics in breast cancer

To examine the difference of Hsp27 expression between normal tissues and breast cancers, and the association between Hsp27 and clinic-pathological characteristics namely TNM stage, nodal metastasis and molecular subtypes in TCGA, UALCAN (The University of ALabama at Birmingham Cancer data analysis Portal, https://ualcan.path.uab.edu/index.html, accessed on May 11st, 2023), a database providing easy access to publicly available cancer OMICS data, was utilized (15).

Impact of Hsp27 on survival in patients with breast cancer

To identify the prognostic value of Hsp27 in breast cancer, KM Plotter database (http://kmplot.com, accessed on May 11st, 2023), that contained gene expression and survival information from 4929 patients with breast cancer, was employed. According to the median Hsp27 expression at the transcriptional level, patients were divided into two groups of high and low expression, and survival curves related to overall survival (OS), replase-free survival (RFS) and distant metastasis-free survival (DMFS) were created with log-rank test (16). Additionally, PrognoScan, a database collecting publicly available cancer microarray datasets and employing the minimum P-value approach for grouping patients for survival analysis (http://dna00.bio.kyutech.ac.jp/PrognoScan/index.html, accessed on May 13^th, 2023), was utilized to analyze the correlation between Hsp27 expression and survival in several individual datasets (17).

Association between Hsp27 and drug resistence in patients with breast cancer

To examine the relation between Hsp27 at transcriptional level and cancer drug response, CTR-DB (http://ctrdb.ncpsb.org.cn/home, accessed on May 11st, 2023), a database covering 28 tumor types, 123 drugs and 5139 tissue samples and providing drug response information, was employed (18).

Identification of Hsp27-related genes and functional enrichment

To identify Hsp27-related genes, Pearson’s correlation coefficients were calculated, and a volcano plot and heat maps were created in LinkedOmics database (http://www.linkedomics.org/login.php, accessed on May 11st, 2023), an online platform for analyzing 32 TCGA cancer-associated multi-dimensional datasets (19). Hsp27-related genes with correlation coefficients greater than 0.3 or less than -0.3 were extracted. Then, ClueGO and EnrichmentMap in NetworkAnalyst 3.0 tool, (accessed on May 11st, 2023) were unilized to perform functional enrichment analyses based on these potentially co-expressed genes (20). Otherwise, the GeneMANIA database (http://www.genemania.org, accessed on May 11st, 2023) was applied to build an interactive network related to Hsp27. A protein-protein interaction (PPI) network was constructed using the STRING online database (https://string-db.org/, accessed on May 11st, 2023), which aims to integrate all known and predicted associations between proteins, including both physical interactions and functional associations.

Patients and tissue samples

A total of 103 patients diagnosed with LABC who received anthracycline (i.e., doxorubicin or epirubicin)-based NAC (median age: 51 years; range: 28-80 years) in the Cancer Hospital of Shantou University Medical College between December 2008 and September 2013 were included, and we collected the formalin-fixed, paraffin-embedded specimens of both pre-surgery needle core biopsies and post-surgery mastectomized tumors, and clinic-pathological information including age, menopausal state, tumor size (T stage), lymph node status (N stage), clinical TNM stage, histological type, expression of ER, progesterone receptor (PgR), human epidermal growth factor receptor-2 (HER-2) and Ki67. The clinical tumor stage (TNM stage) was defined in accordance with the American Joint Committee on Cancer (AJCC) 6^th edition of the Cancer Staging Manual (2002). The molecular classification was grouped in accordance with the 13^th St. Gallen International Breast Cancer Conference Expert Panel (2013) (21). ER, PgR, HER-2 and Ki67 scores were evaluated by immunohistochemical (IHC) staining in pre-surgery samples, and Ki67 scores were also detected in post-surgery samples. Immunostaining and scoring ER, PgR, HER-2 and Ki67 were conducted by experienced pathologists during the clinical practice, that followed the standard operating procedure in the pathology department. The HER-2 and Ki67 antibodies were purchased from MXB (KIT-0043 and MAB-0672, Fuzhou, China); ER and PgR were purchased from Roche (790-4325 and 790-4296, Basel, Switzerland). The cut-off value for ER and PgR positivity was set at ≥10% (22). HER-2 was identified positive if IHC staining showed a score of 3, or if fluorescence in situ hybridization (FISH) showed >2.2-fold amplification of the HER-2 gene in tumors, and the FISH test was only performed for tumors that were scored at 2 in IHC staining (23). A cut-off of Ki67 to distinguish luminal A from luminal B was 14% (24). Prognosis was determined by annual phone interview and medical records. Informed consent for the use of their samples was obtained from all patients. This study was approved by the medical ethics committee of the Cancer Hospital of Shantou University Medical College (No. 201733).

Immunohistochemistry assay

IHC for Hsp27 and TopoIIα were carried out using a standard Envision complex method (25). Briefly, needle core biopsies were fixed in 10% buffered formalin and embedded in paraffin. After deparaffinization and rehydration, endogenous peroxidase activity was blocked with 3% hydrogen peroxide for 10 min. Then antigen retrieval was performed by heating the slides in a microwave oven for 20 min in an EDTA buffer (pH 8.0). The dilution of Hsp27 and TopoIIα primary antibodies (ab2790, ab52934, Abcam, Cambridge, UK) was 1:400 and 1:200, respectively. IHC staining was carried out by a GTVision antibody complex (anti-mouse/rabbit) method using a GTVision™ Detection Kit (GK5006, Gene Tech, Shanghai, China) and 3,3’-diaminobenzidine (DAB) as the chromogen substrate. A negative control was obtained by replacing the primary antibody with mouse IgG1-kappa monoclonal antibody (ab18443, Abcam, Cambridge, UK).

IHC staining for Hsp27 was scored, as described by a combination of intensity (0, no staining; 1, weak staining; 2, moderate staining; 3, strong staining) and proportion (0, < 5% of tumor cells stained; 1, 5-25% positive cells; 2, 26-50% positive cells; 3, 51-75% positive cells; 4, more than 76% positive cells) scores (26). If the multiplication product between staining intensity and the proportion of positive cells was >4 (the median multiplication value), expression was defined as Hsp27-positive; otherwise, samples were regarded as Hsp27-negative. TopoIIα status was scored for positivity in four random fields (×400 magnification) to get the total percent-score. Then, the median values of TopoIIα were chosen as cut-offs for high or low expression (27). Two pathologists independently assessed the cellular location and immunostaining intensity, and a third observer determined if any discrepancies happened.

Statistical analysis

All analyses were performed in SPSS (version 13.0). Associations of Hsp27 expression with TopoIIα and clinicopathological features were analyzed using the chi-square test or Fisher’s exact test. Disease-free survival (DFS) in months was defined as the period from surgery to the date of relapse or metastasis or last contact, which came first. OS in months was defined as the period from surgery to the date of last contact or death from any cause, which came first. Kaplan-Meier curves were created according to Hsp27 expressions and tested by log-rank test. Univariate and Multivariate Cox regressions were used to study the effects of Hsp27 expression on DFS and OS. P < 0.05 was considered to indicate statistical significance, and all resulting P-values were two-tailed.

Results

Associations between Hsp27 at the transcriptional level and clinic-pathological features in breast cancer

According to UALCAN database, there was a higher Hsp27 expression in breast cancers compared with normal group (Figure 1A), and an elevation of Hsp27 was in tumors at Stage 3 compared to those Stage 1 and 2 (Figure 1B). A fail of statistical significance at Stage 4 might be related to a small sample size. In term of nodal metastasis status, an increase of Hsp27 expression was seen along with increased numbers of nodes with metastasis (Figure 1C). Regarding molecular subtypes, Hsp27 expression was the highest in luminal breast cancer, followed by that in Her2-positive types and TNBC (Figure 1D). Moreover, according to CTR-DB, a high expression level of Hsp27 significantly indicated resistance to TFAC (Cyclophosphamide, Doxorubicin, fluorouracil, paclitaxel) and epirubicin, and the area under the AUC curve was greater than 0.7 (Supplementary Figure 1). In short, these findings suggest that Hsp27 mRNA levels were related to an aggressive character and drug resistance in patients with breast cancer.

FIGURE 1

Figure 1 Box plots show Hsp27 expression among different clinical parameters using the UALCAN database. Difference of Hsp27 expression was separately evaluated for (A) sample types, (B) cancer stages, (C) nodal metastasis status and (D) breast cancer subclasses. N0: no regional lymph node metastasis; N1: metastases in 1 to 3 axillary lymph nodes; N2: metastases in 4 to 9 axillary lymph nodes; N3: metastases in 10 or more axillary lymph nodes. **P<0.01, ***P<0.001.

Prognostic values of Hsp27 at the transcriptional level in breast cancer

According to Kaplan-Meier plotter database, elevation of Hsp27 at the transcriptional level was associated significantly with a poor prognosis in breast cancer (HRs, 95% CIs for OS: 1.41, 1.16-1.71, P=5.6e-04; RFS: 1.36, 1.23-1.5, P=3.9e-09; DMFS: 1.26, 1.08-1.47, P=0.0031. Figures 2A–C). Similar findings were seen when patients were separately restricted to those who did not receive chemotherapy (Figure 2D) and those who received chemotherapy (Figure 2E). Regarding individual datasets, consistent findings remained using PrognoScan database (Figures 2F–H; Supplementary Figure 2). All these results demonstrate that Hsp27 expression might be correlated with a poor outcome in patients with breast cancer.

FIGURE 2

Figure 2 Survival curves were created to evaluate the prognostic value of Hsp27 in breast cancer. (A–C) OS, RFS and DMFS curves were created despite chemotherapy in the Kaplan Meier plotter database. (D) OS curves were created for patients without chemotherapy in the Kaplan Meier plotter database. (E) OS curves were created for patients receiving chemotherapy in the Kaplan Meier plotter database. (F–H) OS, RFS and DMFS curves were created using the PrognoScan database (dataset: GSE1456). All dotted lines represented confidence bands.

Hsp27-related genes and functions in breast cancer

According to LinkedOmics database, a total of 3079 genes and 3072 genes were positively and negatively correlated with Hsp27, respectively (false positive rate, FDR<0.05, Figure 3A). Ranking by absolute values of pearson correlation coefficients, the top 50 genes positively and negatively related to Hsp27 were separately shown in Figures 3B, C. An overall description of co-expressed genes is detailed in Supplementary Table 1. Enrichment analyses of KEGG pathway based on potentially co-expressed genes indicated related pathways such as “cAMP signaling pathway”, “metabolic pathways” and “Oxytocin signaling pathway” (Supplementary Table 2; Figure 3D). Enrichment analyses based on GO database showed that the related biological processes were gamete generation, anion transport and sexual reproduction (Supplementary Table 3), and the related cellular components were cytosol, vesicle and cytoplasmic vesicle, and the related molecular functions were purine nucleotide binding, purine ribonucleotide binding (Supplementary Tables 4, 5). Further, we used GeneMANIA to construct a gene-gene interaction network, that showed twenty genes with the highest altered frequencies related to Hsp27, including MAPKAPK2, MAPKAPK3, MAPKAPK5 (Figure 3E). Consistently, these genes were mostly observed in a protein-protein interaction (PPI) network generated on the STRING database (Figure 3F). These findings suggest that Hsp27 was related to biological metabolism, and mitogen-activated protein kinases.

FIGURE 3

Figure 3 Hsp27-relevant genes and functions were estimated in breast cancer. (A) The global Hsp27 highly correlated genes were identified by Pearson test in the LinkedOmics database (dataset: TCGA). Heat maps separately showed the top-ranked 50 genes positively (B) and negatively (C) correlated with Hsp27 in breast cancer. (D) Functional enrichment analyses based on KEGG pathway were performed for Hsp27 related genes. (E) The gene-gene interaction network for Hsp27 (i.e., HSPB1) was constructed using GeneMANIA. (F) The PPI network of Hsp27 was generated using STRING.

Associations of Hsp27 proteins with clinicopathological features, TopoIIα and therapeutic efficacy in patients with LABC

Hsp27 and TopoIIα proteins were detected in tumors by IHC assays (Figure 4). As shown in Table 1, out of 103 LABC patients, 48 (46.6%) were Hsp27-positive, Hsp27 (+), while 55 (53.4%) patients were Hsp27-negative, Hsp27 (-). The percentage of Hsp27 (+) tumors comprised 22.2% of the N1 group (4/18) and 51.8% of the N2-N3 group (44/85), leading to a significant association between Hsp27 and lymph node status (P = 0.022). Besides, Hsp27 was related to histological grade (P< 0.001); there was only 17.1% of well or moderately differentiated tumors (7/48), but 66.1% of poorly differentiated tumors (41/62) being detected Hsp27 (+). Additionally, Hsp27 (-) tumors showed a percentage of 60.3% (47/78) to have NAC-reduced Ki67, but only a percentage of 32.0% to have NAC-elevated Ki67 (P=0.014). These findings again confirm the aggressive role of Hsp27 in breast cancer and its role involving in chemotherapeutic resistance. No significant differences were found for age, menopausal status, TNM stage, ER status, PgR status, HER-2 status and molecular subtypes.

FIGURE 4

Figure 4 Hsp27 and TopoIIα proteins were detected by IHC in breast cancers. Tumor tissues were separately stained (A) negative, (B) week, (C) moderate, (D) strong for Hsp27, and (E) negative, (F) week, (G) moderate, (H) strong for TopoIIα. The magnification was 400×.

TABLE 1

Table 1 Relationship of Hsp27 expression with clinicopathological features,TopoIIα and therapeutic effect of anthracycline-based NAC in LABC patients.

As shown in Table 1, 74 out of 103 tumors (71.8%) displayed TopoIIα-positive, among which 23 (23/74, 31.1%) were Hsp27 (+) tumors. However, 25 out of 29 TopoIIα-negative tumors displayed Hsp27 (+) (25/29, 86.2%), demonstrating a negative correlation between TopoIIα and Hsp27. Additionally, compared with patients who showed complete response (CR) and partial response (PR) to chemotherapy, patients with stable or progressive disease had a higher Hsp27 expression (32/79, 40.5% vs. 16/24, 66.7%, P = 0.024). All patients (11/11, 100%, P = 0.001) in the Hsp27 (-) group showed pathologic complete response (pCR) following NAC treatment.

Effect of Hsp27 expression on DFS and OS

In terms of DFS, the follow-up period ranged from 3 to 65 months (median: 31 months). Patients with Hsp27(+) tumors tended to have a shorter DFS than those with Hsp27 (-) tumors (P = 0.001, Figure 5A), consistent with the result of the univariate Cox regression (HR: 3.028, 95%CI: 1.506-6.086, P= 0.002, Table 2). Conversely, a shorter DFS was associated with TopoIIα negativity (HR: 0.270; 95%CI: 0.139~0.521; P< 0.001, Table 2). In the multivariate analysis (Table 2, right), only TopoIIα (HR: 0.316; 95%CI: 0.162~0.617; P = 0.001) remained while Hsp27 failed to be an independent prognostic factor for DFS.

FIGURE 5

Figure 5 Kaplan-Meier curves were created to examine the associations of Hsp27 and TopoIIα proteins with DFS or OS. According to the expression of Hsp27, survival curves were separately created to compare (A) disease-free survival rate and (B) overall survival rate. According the expression of Hsp27 and TopoIIa in combination, survival curves were separately created to compare (C) disease-free survival rate and (D) overall survival rate.

TABLE 2

Table 2 Cox proportional hazard regression analysis of DFS in LABC.

Regarding OS, the follow-up period ranged from 10 to 68 months (median: 34 months). Hsp27 (+) patients showed a lower OS than those Hsp27 (-) patients (P= 0.001, Figure 5B), in accordance with Cox regression analysis (HR: 6.213; 95% CI: 1.765~21.869; P = 0.004; Table 3). Additionally, TopoIIα (HR: 0.147; 95% CI: 0.051~0.425; P< 0.001) was a beneficial predictor. Furthermore, as shown in Table 3, the multivariate analysis showed that only TopoIIα (HR: 0.214; 95%CI: 0.074~0.619; P = 0.004) was an independent predictor for OS in all 103 LABC patients.

TABLE 3

Table 3 Cox proportional hazard regression analysis of OS in LABC.

Association between expression of TopoIIα and Hsp27 in combination and resistence to anthracycline-based NAC in LABC patients

Considering the prognostic values of Hsp27 and TopoIIα, their prognostic value in combination on survival was further examined. There were 4 patients with negative expression of both TopoIIα and Hsp27 [TopoIIα (-) Hsp27 (-), 4/103], 25 patients with negative expression of TopoIIα and positive expression of Hsp27 [TopoIIα (-) Hsp27 (+), 25/103, 24.3%], 51 patients with positive expression of TopoIIα and negative expression of Hsp27 [TopoIIα (+) Hsp27 (-), 51/103, 49.5%] and 23 patients with positive expression of both Hsp27 and TopoIIα [TopoIIα (+) Hsp27 (+), 23/103, 22.3%] (Table 4). Kaplan-Meier analysis indicated that patients characterized by TopoIIα (-) Hsp27 (+) had significantly lower DFS and OS (P< 0.001 and 0.001, respectively) than those in the other three groups (Figures 5C, D).

TABLE 4

Table 4 Kaplan-Meier analysis for patients with different expressions of TopoIIα and Hsp27.

Discussion

The NAC administration represents a pivotal therapeutic approach for patients with LABC. Nevertheless, the development of resistance to NAC remains a paramount clinical obstacle. Thus, identification of biomarkers capable of accurately predicting the responsiveness of patients to NAC is an urgent and justifiable pursuit.

The overexpression of Hsp27 has been linked to a range of biological processes such as aggressiveness, metastasis, apoptotic cell death inhibition immunity, drug resistance, and tumor prognosis, etc (28, 29). Although several cell line studies have indicated a potential association between Hsp27 and the responsiveness of breast cancer to anthracyclines (30, 31), less clinical investigations have been conducted to confirm the prognostic role of Hsp27 in breast cancers. In patients with node-negative breast cancer, an early study failed to see significant correlations between Hsp27 expression and histopathological features or Ki67 or DFS, but Hsp27 was associated with a short OS (7), that could not even be seen in patients with node-positive breast carcinomas (32). As Ki67 has been clinically utilized to define highly proliferative breast cancers (33), we previously reported that Ki67 decrease during neoadjuvant chemotherapy predicts favorable relapse-free survival in patients with LABC (34). The present study, otherwise, proposes that Ki67 in Hsp27-positive LABC patients is significantly elevated by NAC, suggesting that Hsp27 may predict poor prognosis of LABC patients treated with NAC. The current study additionally found a significant higher expression of Hsp27 in poorly differentiated LABC, which is consistent with our previous report showing Hsp27 overexpression in breast cancer stem-like cells (35).

One previous study similarly included LABC patients (n = 60) who were treated with doxorubicin- or epirubicin-based neoadjuvant monochemotherapy, indicated nuclear Hsp27 expression was significantly increased in invasive cells after doxorubicin/epirubicin administration, but did not find correlations between Hsp27 and the clinical and pathological response to neoadjuvant therapy, and the prognostic value of Hsp27 (11). However, this current study argued an existing association of Hsp27 with aggressive features, response to neoadjuvant therapy, and survival (using KM curves), albeit Hsp27 might not be an independent prognostic indicator for DFS or OS. The inconsistent findings may be due to the difference of Hsp27 antibody, evaluation method, genetic background, and treatments etc.

Previous studies have suggested that high expression of Hsp27 may induce drug resistance of tumor cells (14). Consistently, we found that patients with stable and progressive disease had a higher percentage of Hsp27 expression, and all patients in the Hsp27 negative group showed pCR. In addition, our study suggests that the LABC patients with high levels of Hsp27 tend to have a shorter DFS and OS than those with negative Hsp27 expression, but might not be an independent prognostic indicator, that could be closely related to TopoIIα.

A previous study demonstrated that Hsp27 overexpression may inhibit TopoIIα expression in breast cancer cells (14). Our results further pointed out that patients with Hsp27-positive and TopoIIα-negative tumors benefited from a considerably lower DFS and OS. However, the underlying molecular mechanisms remain unclear and warrant further study. In summary, our results suggest that Hsp27 expression is associated with the NAC outcome in LABC, though it is not an independent prognostic factor for DFS and OS. Further prospective investigations of a larger number of LABC patients treated with NAC will be necessary to determine whether Hsp27 could predict efficacy of treatment. More importantly, we co-evaluate the expression status of Hsp27 and TopoIIα to identify that higher Hsp27/lower TopoIIα could prognose worse outcome for NAC resistance in breast cancer.

Conclusions

Although either Hsp27 or TopoIIα might be able to predict LABC patient responsiveness to anthracycline-based neoadjuvant treatments, combination of these two biomarkers could more precisely prognose the outcome of these patients.

Data availability statement

The publicly available datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material. The dataset from the Cancer Hospital of Shantou UniversityMedical College is available on reasonable request.

Ethics statement

The studies involving human participants were reviewed and approved by Cancer Hospital of Shantou University Medical College. The patients/participants provided their written informed consent to participate in this study.

Author contributions

(I) Conception and design: YZ, YC. (II) Administrative support: YC, WH. (III) Provision of study materials or patients: YX, CH, WH, LH. (IV) Collection and assembly of data: FZ. (V) Funding Acquisition: YZ, CH. (VI) Manuscript writing: All authors. (VII) Final approval of manuscript: All authors. All authors contributed to the article and approved the submitted version.

Funding

This research was supported by (1) the Medical Scientific Research Foundation of Guangdong Province, China (Grant No. B2018049) and the Youth Research Grant from Shantou University Medical College Cancer Hospital (Grant No. 2018A002, 2018-36), and the Shantou Municipal Science-Technology Bureau (Grant No. 2022-122) to YZ; (2) the Shantou Municipal Science-Technology Bureau (Grant No. 2017-166; 2022-169) and the Project supported by the Natural Science Foundation of Guangdong Province, China (Grant No. 2018A0303130279) to Chaoqun Hong; (3) the Guangdong Basic and Applied Basic Research Foundation (Grant No. 2020A1515010094) and the Special Fund for Science and Technology of Guangdong Province in 2020 (Grant No. 200110115891683) to FZ; (4) the Special Fund for Science and Technology of Guangdong Province in 2021 (Grant No. 210729156901814) and the National Natural Science Foundation of China (Grant No. 82272670) to YC.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2023.1067179/full#supplementary-material

References

1. Chen W, Sun K, Zheng R, Zeng H, Zhang S, Xia C, et al. Cancer incidence and mortality in China, 2014. Chin J Cancer Res (2018) 30(1):1–12. doi: 10.21147/j.issn.1000-9604.2018.01.01

CrossRef Full Text | Google Scholar

2. Schlam I, Tarantino P, Tolaney SM. Overcoming resistance to HER2-directed therapies in breast cancer. Cancers (Basel) (2022) 14(16):3996. doi: 10.3390/cancers14163996

CrossRef Full Text | Google Scholar

3. Hurvitz SA, Martin M, Symmans WF, Jung KH, Huang CS, Thompson AM, et al. Neoadjuvant trastuzumab, pertuzumab, and chemotherapy versus trastuzumab emtansine plus pertuzumab in patients with HER2-positive breast cancer (KRISTINE): a randomised, open-label, multicentre, phase 3 trial. Lancet Oncol (2018) 19(1):115–26. doi: 10.1016/S1470-2045(17)30716-7

CrossRef Full Text | Google Scholar

4. Yun CW, Kim HJ, Lim JH, Lee SH. Heat shock proteins: agents of cancer development and therapeutic targets in anti-cancer therapy. Cells (2019) 9(1):60. doi: 10.3390/cells9010060

CrossRef Full Text | Google Scholar

5. Lampros M, Vlachos N, Voulgaris S, Alexiou GA. The Role of Hsp27 in Chemotherapy Resistance. Biomedicines (2022) 10(4). doi: 10.3390/biomedicines10040897

CrossRef Full Text | Google Scholar

6. Alberti G, Vergilio G, Paladino L, Barone R, Cappello F, Conway de Macario E, et al. The chaperone system in breast cancer: roles and therapeutic prospects of the molecular chaperones hsp27, hsp60, hsp70, and hsp90. Int J Mol Sci (2022) 23(14):7792. doi: 10.3390/ijms23147792

CrossRef Full Text | Google Scholar

7. Thanner F, Sutterlin MW, Kapp M, Rieger L, Morr AK, Kristen P, et al. Heat shock protein 27 is associated with decreased survival in node-negative breast cancer patients. Anticancer Res (2005) 25(3A):1649–53.

Google Scholar

8. Love S, King RJ. A 27 kDa heat shock protein that has anomalous prognostic powers in early and advanced breast cancer. Br J Cancer (1994) 69(4):743–8. doi: 10.1038/bjc.1994.140

CrossRef Full Text | Google Scholar

9. Oesterreich S, Hilsenbeck SG, Ciocca DR, Allred DC, Clark GM, Chamness GC, et al. The small heat shock protein HSP27 is not an independent prognostic marker in axillary lymph node-negative breast cancer patients. Clin Cancer Res (1996) 2(7):1199–206.

Google Scholar

10. Vargas-Roig LM, Gago FE, Tello O, Aznar JC, Ciocca DR. Heat shock protein expression and drug resistance in breast cancer patients treated with induction chemotherapy. Int J Cancer (1998) 79(5):468–75. doi: 10.1002/(SICI)1097-0215(19981023)79:5<468::AID-IJC4>3.0.CO;2-Z

CrossRef Full Text | Google Scholar

11. Nadin SB, Sottile ML, Montt-Guevara MM, Gauna GV, Daguerre P, Leuzzi M, et al. Prognostic implication of HSPA (HSP70) in breast cancer patients treated with neoadjuvant anthracycline-based chemotherapy. Cell Stress Chaperones (2014) 19(4):493–505. doi: 10.1007/s12192-013-0475-2

CrossRef Full Text | Google Scholar

12. Aebi S, Karlsson P, Wapnir IL. Locally advanced breast cancer. Breast (2022) 62 Suppl 1(Suppl 1):S58–62. doi: 10.1016/j.breast.2021.12.011

CrossRef Full Text | Google Scholar

13. Wang J, Xu B, Yuan P, Zhang P, Li Q, Ma F, et al. TOP2A amplification in breast cancer is a predictive marker of anthracycline-based neoadjuvant chemotherapy efficacy. Breast Cancer Res Treat (2012) 135(2):531–7. doi: 10.1007/s10549-012-2167-5

CrossRef Full Text | Google Scholar

14. Hansen RK, Parra I, Lemieux P, Oesterreich S, Hilsenbeck SG, Fuqua SA, et al. Hsp27 overexpression inhibits doxorubicin-induced apoptosis in human breast cancer cells. Breast Cancer Res Treat (1999) 56(2):187–96. doi: 10.1023/A:1006207009260

CrossRef Full Text | Google Scholar

15. Chandrashekar DS, Karthikeyan SK, Korla PK, Patel H, Shovon AR, Athar M, et al. UALCAN: An update to the integrated cancer data analysis platform. Neoplasia (2022) 25:18–27. doi: 10.1016/j.neo.2022.01.001

CrossRef Full Text | Google Scholar

16. Gyorffy B. Survival analysis across the entire transcriptome identifies biomarkers with the highest prognostic power in breast cancer. Comput Struct Biotechnol J (2021) 19:4101–9. doi: 10.1016/j.csbj.2021.07.014

CrossRef Full Text | Google Scholar

17. Mizuno H, Kitada K, Nakai K, Sarai A. PrognoScan: a new database for meta-analysis of the prognostic value of genes. BMC Med Genomics (2009) 2:18. doi: 10.1186/1755-8794-2-18

CrossRef Full Text | Google Scholar

18. Liu Z, Liu J, Liu X, Wang X, Xie Q, Zhang X, et al. CTR-DB, an omnibus for patient-derived gene expression signatures correlated with cancer drug response. Nucleic Acids Res (2022) 50(D1):D1184–99. doi: 10.1093/nar/gkab860

CrossRef Full Text | Google Scholar

19. Vasaikar SV, Straub P, Wang J, Zhang B. LinkedOmics: analyzing multi-omics data within and across 32 cancer types. Nucleic Acids Res (2018) 46(D1):D956–63. doi: 10.1093/nar/gkx1090

CrossRef Full Text | Google Scholar

20. Zhou G, Soufan O, Ewald J, Hancock REW, Basu N, Xia J. NetworkAnalyst 3.0: a visual analytics platform for comprehensive gene expression profiling and meta-analysis. Nucleic Acids Res (2019) 47(W1):W234–41. doi: 10.1093/nar/gkz240

CrossRef Full Text | Google Scholar

21. Goldhirsch A, Winer EP, Coates AS, Gelber RD, Piccart-Gebhart M, Thurlimann B, et al. Personalizing the treatment of women with early breast cancer: highlights of the St Gallen International Expert Consensus on the Primary Therapy of Early Breast Cancer 2013. Ann Oncol (2013) 24(9):2206–23. doi: 10.1093/annonc/mdt303

CrossRef Full Text | Google Scholar

22. Prabhu JS, Korlimarla A, Desai K, Alexander A, Raghavan R, Anupama C, et al. A Majority of Low (1-10%) ER Positive Breast Cancers Behave Like Hormone Receptor Negative Tumors. J Cancer (2014) 5(2):156–65. doi: 10.7150/jca.7668

CrossRef Full Text | Google Scholar

23. Wolff AC, Hammond ME, Schwartz JN, Hagerty KL, Allred DC, Cote RJ, et al. American Society of Clinical Oncology/College of American Pathologists guideline recommendations for human epidermal growth factor receptor 2 testing in breast cancer. Arch Pathol Lab Med (2007) 131(1):18–43. doi: 10.5858/2007-131-18-ASOCCO

CrossRef Full Text | Google Scholar

24. Cheang MC, Chia SK, Voduc D, Gao D, Leung S, Snider J, et al. Ki67 index, HER2 status, and prognosis of patients with luminal B breast cancer. J Natl Cancer Inst (2009) 101(10):736–50. doi: 10.1093/jnci/djp082

CrossRef Full Text | Google Scholar

25. Patel RM, Folpe AL. Immunohistochemistry for human telomerase reverse transcriptase catalytic subunit (hTERT): a study of 143 benign and Malignant soft tissue and bone tumours. Pathology (2009) 41(6):527–32. doi: 10.1080/00313020903121614

CrossRef Full Text | Google Scholar

26. Ghosh A, Lai C, McDonald S, Suraweera N, Sengupta N, Propper D, et al. HSP27 expression in primary colorectal cancers is dependent on mutation of KRAS and PI3K/AKT activation status and is independent of TP53. Exp Mol Pathol (2013) 94(1):103–8. doi: 10.1016/j.yexmp.2012.09.001

CrossRef Full Text | Google Scholar

27. Mukherjee A, Shehata M, Moseley P, Rakha E, Ellis I, Chan S, et al. Topo2alpha protein expression predicts response to anthracycline combination neo-adjuvant chemotherapy in locally advanced primary breast cancer. Br J Cancer (2010) 103(12):1794–800. doi: 10.1038/sj.bjc.6605960

CrossRef Full Text | Google Scholar

28. Banerjee S, Lin CF, Skinner KA, Schiffhauer LM, Peacock J, Hicks DG, et al. Heat shock protein 27 differentiates tolerogenic macrophages that may support human breast cancer progression. Cancer Res (2011) 71(2):318–27. doi: 10.1158/0008-5472.CAN-10-1778

CrossRef Full Text | Google Scholar

29. Wang X, Chen M, Zhou J, Zhang X. HSP27, 70 and 90, anti-apoptotic proteins, in clinical cancer therapy (Review). Int J Oncol (2014) 45(1):18–30. doi: 10.3892/ijo.2014.2399

CrossRef Full Text | Google Scholar

30. Diaz-Chavez J, Fonseca-Sanchez MA, Arechaga-Ocampo E, Flores-Perez A, Palacios-Rodriguez Y, Dominguez-Gomez G, et al. Proteomic profiling reveals that resveratrol inhibits HSP27 expression and sensitizes breast cancer cells to doxorubicin therapy. PloS One (2013) 8(5):e64378. doi: 10.1371/journal.pone.0064378

CrossRef Full Text | Google Scholar

31. Sims JT, Ganguly SS, Bennett H, Friend JW, Tepe J, Plattner R, et al. Imatinib reverses doxorubicin resistance by affecting activation of STAT3-dependent NF-kappaB and HSP27/p38/AKT pathways and by inhibiting ABCB1. PloS One (2013) 8(1):e55509. doi: 10.1371/journal.pone.0055509

CrossRef Full Text | Google Scholar

32. Tetu B, Brisson J, Landry J, Huot J. Prognostic significance of heat-shock protein-27 in node-positive breast carcinoma: an immunohistochemical study. Breast Cancer Res Treat (1995) 36(1):93–7. doi: 10.1007/BF00690189

CrossRef Full Text | Google Scholar

33. Nielsen TO, Leung SCY, Rimm DL, Dodson A, Acs B, Badve S, et al. Assessment of Ki67 in Breast Cancer: Updated Recommendations From the International Ki67 in Breast Cancer Working Group. J Natl Cancer Inst (2021) 113(7):808–19. doi: 10.1093/jnci/djaa201

CrossRef Full Text | Google Scholar

34. Chen C, Zhang Y, Huang Z, Wu J, Huang W, Zhang G. Decrease in the Ki67 index during neoadjuvant chemotherapy predicts favorable relapse-free survival in patients with locally advanced breast cancer. Cancer Biol Med (2019) 16(3):575–86. doi: 10.20892/j.issn.2095-3941.2018.0423

CrossRef Full Text | Google Scholar

35. Li G, Zhao F, Cui Y. Proteomics using mammospheres as a model system to identify proteins deregulated in breast cancer stem cells. Curr Mol Med (2013) 13(3):459–63.

Google Scholar