Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Pharmacol., 09 February 2026

Sec. Pharmacology of Anti-Cancer Drugs

Volume 17 - 2026 | https://doi.org/10.3389/fphar.2026.1752328

This article is part of the Research TopicAdvancement of Chemotherapy in Breast Cancer: Predictive Markers, Resistance Mechanism and Therapeutic StrategiesView all 15 articles

Tumour microenvironment and HER2-low status dictate response and resistance to anthracycline-taxane chemotherapy in premenopausal TNBC: a retrospective multicohort study

Shuanglong Cai&#x;Shuanglong Cai1Shaohong Yu&#x;Shaohong Yu2Quan Zhou&#x;Quan Zhou3Fangyuan KuangFangyuan Kuang1Lingyan RaoLingyan Rao1Zhen FangZhen Fang1Xiaoxin Zheng,,
&#x;Xiaoxin Zheng4,5,6*Yong Shi
&#x;Yong Shi1*Jingdan Li
&#x;Jingdan Li1*
  • 1Department of Thyroid and Breast Surgery, Comprehensive Breast Health Center, The First Affiliated Hospital of Lishui University, Lishui People’s Hospital, Lishui, Zhejiang, China
  • 2Peking Union Medical College Hospital, Peking Union Medical College and Chinese Academy of Medical Sciences, Beijing, China
  • 3Department of Gynecology, Fujian Provincial Hospital, Fuzhou University Affiliated Provincial Hospital, Fuzhou, Fujian, China
  • 4Department of Cardiology, Fuzhou University Affiliated Provincial Hospital, Shengli Clinical Medical College of Fujian Medical University, Fujian Provincial Hospital, Fuzhou, Fujian, China
  • 5Fujian Provincial Key Laboratory of Cardiovascular Disease, Fujian Cardiovascular Institute, Fujian Provincial Center for Geriatrics, Fujian Clinical Medical Research Center for Cardiovascular Diseases, Fuzhou, Fujian, China
  • 6Fujian Heart Failure Center Alliance, Fuzhou, Fujian, China

Background: Premenopausal women with triple-negative breast cancer (TNBC) exhibit considerable heterogeneity in their response to standard anthracycline and taxane-based chemotherapy, yet the underlying mechanisms remain poorly understood. We aimed to investigate the combined role of the tumour immune microenvironment and HER2-low status in predicting chemosensitivity and intrinsic resistance in this specific population.

Methods: We retrospectively analysed data from 767 premenopausal patients with TNBC across two Chinese medical centres. All patients underwent primary surgery followed by adjuvant chemotherapy based on anthracyclines and taxanes. Patients were randomly assigned to training and internal validation cohorts. Independent predictors of DFS were identified using multivariable Cox proportional hazards regression, and a nomogram was constructed accordingly. The model’s discrimination was assessed using the concordance index (C-index) and time-dependent receiver operating characteristic (ROC) curve analysis, while calibration was evaluated with calibration curves.

Results: Multivariable analysis identified lower stromal tumour-infiltrating lymphocyte (sTIL) expression levels, and HER2 IHC 2+/FISH-negative status as independent factors associated with poorer DFS, besides that T3/T4 staging, higher N staging. A nomogram integrating these four variables demonstrated excellent predictive accuracy for DFS, with a C-index of 0.862 in the training set and 0.861 in the validation set. The area under the ROC curve (AUC) for predicting 3-year DFS was 0.907 and 0.908 in the training and validation sets, respectively.

Conclusion: Our findings reveal that an immune-poor tumour microenvironment and HER2-low biology are key, complementary determinants of intrinsic resistance to standard chemotherapy in premenopausal TNBC. These readily available biomarkers provide a mechanistic rationale for patient stratification, suggesting that sTIL-low tumours might benefit from immunomodulatory strategies, while HER2-low tumours represent a candidate population for novel antibody-drug conjugates. This paradigm shift from empirical to biomarker-informed therapy could help overcome chemoresistance in this high-risk group.

1 Introduction

Breast cancer remains the most commonly diagnosed malignancy and the leading cause of cancer-related mortality among women worldwide (Bray et al., 2024). Triple-negative breast cancer (TNBC), characterized by the absence of estrogen receptor, progesterone receptor, and HER2 amplification, represents a particularly aggressive subtype with limited therapeutic options and a heightened propensity for early recurrence (Bianchini et al., 2016; Foulkes et al., 2010). Within this heterogeneous disease, premenopausal patients often exhibit distinct clinicopathological and molecular features, including enriched immunogenic signatures and homologous recombination deficiency, which may influence both tumour behaviour and therapeutic response (Azim et al., 2012; Bareche et al., 2018).

Anthracycline and taxane-based chemotherapy remains the cornerstone of adjuvant treatment for early-stage TNBC (Liedtke et al., 2008). However, intrinsic and acquired chemoresistance significantly undermines its efficacy, leading to disparate clinical outcomes (Gonzalez-Angulo et al., 2011; Wahba and El-Hadaad, 2015). The mechanisms underlying chemoresistance are multifactorial, encompassing tumour-intrinsic factors such as apoptotic evasion, enhanced DNA damage repair, and cancer stem cell persistence, as well as microenvironmental influences including immune contexture and stromal composition (Holohan et al., 2013; Shibue and Weinberg, 2017; Gottesman, 2002).

Emerging evidence underscores the prognostic significance of stromal tumour-infiltrating lymphocytes (sTILs) in TNBC, where higher densities correlate with improved pathological complete response and survival following chemotherapy (M Geurts et al., 2024; Salgado et al., 2015). Concurrently, The biological and clinical significance of HER2-low breast cancer, defined as immunohistochemistry (IHC) 1+ or IHC 2+ with negative in situ hybridization (FISH), has recently come into focus. Supported by evidence from clinical trials of novel anti-HER2 antibody–drug conjugates, this subset has emerged as a therapeutically actionable entity (Saura et al., 2025; Modi et al., 2025). Growing preliminary evidence further suggests that HER2-low tumours may represent a distinct biological subtype with unique therapeutic vulnerabilities (Tarantino et al., 2020; Agostinetto et al., 2021). Recent studies, including those by Tarantino P and Schettini F et al., have begun to delineate the underlying biology of HER2-low breast cancer across subtypes, including triple-negative breast cancer (Tarantino et al., 2023; Schettini et al., 2025).

Despite these advances, validated biomarkers for stratifying chemosensitivity in premenopausal TNBC remain scarce. Most existing prognostic tools neither account for the unique biology of this demographic nor integrate readily accessible biomarkers reflective of tumour-immune interplay (Sparano et al., 2018; Parker et al., 2009). Therefore, we aimed to develop and validate a clinically applicable nomogram that incorporates key clinicopathological and microenvironmental variables-specifically sTIL density and HER2-low status-to predict disease-free survival (DFS) in premenopausal women with TNBC treated with anthracycline-taxane chemotherapy.

2 Methods

2.1 Study design and patient cohort

This retrospective, multicentre cohort study consecutively enrolled 767 female patients with TNBC from Peking Union Medical College Hospital and Lishui People’s Hospital of Zhejiang Province between 1 Jan 2016, and 31 Dec 2021 (Figure 1).

Figure 1
Flowchart illustrating patient selection criteria and study groups from Peking Union Medical College Hospital and Lishui People's Hospital database between January 1, 2016, and December 31, 2021. Patients included are premenopausal women with newly diagnosed, operable, primary unilateral invasive breast cancer and a triple-negative subtype. Those excluded had metastatic disease at diagnosis, a history of neoadjuvant chemotherapy, or incomplete data. A total of 767 patients treated with postoperative adjuvant chemotherapy were enrolled, divided into training sets (548 patients) and validation sets (219 patients).

Figure 1. Patient selection flowchart for this study.

Inclusion Criteria: (1) Histopathologically confirmed invasive breast carcinoma; (2) Negative for estrogen receptor, progesterone receptor, and HER2 (defined as ER/PR <1% positive cells, HER2 IHC 0 or 1+, or 2+ with negative FISH); (3) Premenopausal status at diagnosis; (4) Underwent curative-intent surgery (modified radical mastectomy or breast-conserving surgery) at initial diagnosis; (5) Received adjuvant chemotherapy based on anthracyclines and taxanes (sequential or combination regimens) postoperatively; (6) Had complete clinical, pathological, and follow-up data.

Exclusion Criteria: (1) Bilateral breast cancer or concurrent other malignancies; (2) History of neoadjuvant chemotherapy; (3) Presence of distant metastasis at diagnosis; (4) Incomplete follow-up data.

2.2 Variable definition and data collection

Collected variables included: age at diagnosis, family history of breast cancer, surgical approach, postoperative pathological T staging, postoperative pathological N staging, tumor grade, lymph-vascular invasion, sTIL expression levels, Her2 expression levels, Ki-67 index, CK5/6 expression, EGFR expression, and receipt of radiotherapy.

HER2 Status Assessment: HER2 protein expression was evaluated by IHC, and equivocal cases (IHC 2+) were reflexively tested by FISH according to the American Society of Clinical Oncology/College of American Pathologists (ASCO/CAP) guidelines (Wolff et al., 2023). To ensure consistency and mitigate inter-institutional variability-a known challenge in HER2-low scoring—all IHC and FISH results underwent a centralized review by two experienced breast pathologists who were blinded to patient outcomes. HER2-low status was defined according to the standard definition as IHC 1+ or IHC 2+/FISH-negative. For primary and sensitivity analyses, patients were categorized into three groups: HER2 IHC 0, HER2 IHC 1+, and HER2 IHC 2+/FISH-negative.

sTIL Assessment: sTILs were evaluated by two experienced pathologists blinded to clinical outcomes, according to international guidelines. sTILs were defined as the percentage of tumour stroma area occupied by mononuclear inflammatory cells (lymphocytes and plasma cells). sTIL density was categorised into three groups: low (≤10%), intermediate (10%–40%), and high (>40%).

2.3 Follow-up strategy

Follow-up started from the date of initial breast cancer diagnosis and ended on 1 June 2025. Follow-up methods included outpatient visits and telephone interviews. Data collected encompassed: local or regional recurrence, distant metastasis, time to first recurrence or metastasis, site of first recurrence or metastasis, and vital status at the last follow-up.

Disease-free survival (DFS) was defined as the time from surgery date to the first occurrence of any of the following events: local recurrence, regional recurrence, distant metastasis, contralateral breast cancer, or death from any cause.

2.4 Statistical analysis

Analyses were performed using R software (version 4.3.1). The patient cohort was randomly split into a training set (70%) and an internal validation set (30%). Categorical variables are presented as numbers (percentages) and compared using the χ2 test or Fisher’s exact test. Continuous variables are presented as mean ± standard deviation or median (interquartile range) based on their distribution and compared using the t-test or Mann-Whitney U test.

In the training set, univariable Cox regression analysis (P < 0.05) was used to screen prognostic factors associated with DFS. Significant variables were subsequently included in a multivariable Cox proportional hazards regression model to identify independent prognostic factors. Results are presented as hazard ratios (HRs) with 95% confidence intervals (CIs).

Based on the final independent prognostic factors, a nomogram for predicting 1-year, 3-year, and 5-year DFS probabilities was constructed using the rms package. The model’s discrimination was assessed using the C-index and time-dependent AUC analysis. Calibration curves were plotted to evaluate the agreement between predicted probabilities and observed outcomes (assessed by the Kaplan-Meier method). For key prognostic variables, Kaplan-Meier survival curves were generated, and differences between groups were compared using the log-rank test.

All tests were two-sided, and a P value <0.05 was considered statistically significant.

3 Results

3.1 Patient characteristics and follow-up

The final cohort comprised 767 premenopausal TNBC patients with a median follow-up of 76.0 months. The median age was 49.0 years. The cohort was randomly divided into training (n = 548) and internal validation (n = 219) sets. Baseline clinicopathological characteristics were largely balanced between the training and validation sets, except for family history, postoperative pathological T staging, Ki-67 index, and EGFR status (all P < 0.05, Table 1). All other variables showed no significant differences (all P > 0.05, Table 1). In the training set, 420 DFS events and 128 deaths were recorded over a median follow-up of 73.7 months; the validation set experienced 170 DFS events and 49 deaths over 75.4 months.

Table 1
www.frontiersin.org

Table 1. The clinical-pathological characteristics of TNBC patients in both the training and validation sets.

3.2 Prognostic nomogram for disease-free survival

The multivariate Cox regression revealed that postoperative pathological T staging (T3+T4,HR = 3.434,95%CI:1.975-5.97,P < 0.001), postoperative pathological N staging (N1,HR = 3.204,95%CI:2.101-4.885,P < 0.001; N2,HR = 5.284,95%CI:2.507-11.139,P < 0.001; N3,HR = 7.587,95%CI:4.381-13.14,P < 0.001),sTIL expression levels (sTIL intermediate expression, HR = 0.408, 95% CI: 0.275-0.604,P < 0.001; sTIL high expression, HR = 0.044, 95% CI: 0.021-0.091,P < 0.001), and Her2 expression levels (IHC 2+/Fish-, HR = 2.392, 95% CI: 1.513-3.782,P < 0.001) were independent predictive factors for DFS in premenopausal TNBC patients (Table 2). A nomogram integrating these variables was constructed to estimate 1-, 3-, and 5-year DFS probabilities (Figure 2). The model exhibited excellent discrimination, with C-indices of 0.862 (training) and 0.861 (validation). Time-dependent receiver operating characteristic (ROC) analysis further confirmed its robust predictive accuracy, with area under the curve (AUC) values for 3-year DFS of 0.907 and 0.908 in the training and validation sets, respectively (Figures 3, 4). Calibration curves demonstrated close alignment between predicted and observed outcomes across all time points (Figures 510). Kaplan-Meier analyses visually affirmed the prognostic stratification afforded by each variable (Figures 1114).

Table 2
www.frontiersin.org

Table 2. COX Univariate and Multivariate Analysis of Disease-free survival for TNBC patients.

Figure 2
Nomogram for predicting disease-free survival (DFS) probabilities, including scales for postoperative pathological T and N staging, sTIL and Her2 expression levels. It provides total points with corresponding 1-year, 3-year, and 5-year DFS probabilities.

Figure 2. Prognostic nomogram for predicting disease-free survival time of premenopausal TNBC patients.

Figure 3
ROC curve comparing model performance at different time intervals with three colored lines: yellow for one year (AUC 0.790), blue for three years (AUC 0.907), and red for five years (AUC 0.909). Sensitivity is plotted against one minus specificity.

Figure 3. Roc curve of the training set.

Figure 4
ROC curve plot showing sensitivity versus one minus specificity. Three lines represent different time intervals: 1-year in orange (AUC: 0.841), 3-year in blue (AUC: 0.908), and 5-year in red (AUC: 0.918). A diagonal reference line is included.

Figure 4. Roc curve of the validation set.

Figure 5
Calibration plot showing predicted probability on the x-axis and actual probability on the y-axis. Data points are near the diagonal line, indicating good model calibration. Error bars represent uncertainty.

Figure 5. Calibration curve of the training set at 1 Year.

Figure 6
Calibration plot showing predicted versus actual probabilities with a dashed diagonal reference line. Data points with error bars are shown along the solid line, indicating prediction accuracy. The actual probabilities range from zero to one.

Figure 6. Calibration curve of the training set at 3 Years.

Figure 7
Scatter plot displaying actual probability versus predicted probability with a line of best fit. Blue data points with error bars show deviation from the diagonal line, indicating prediction accuracy. Axis ranges from 0.0 to 1.0.

Figure 7. Calibration curve of the training set at 5 Years.

Figure 8
Calibration plot showing actual probability versus predicted probability. Data points are plotted along with vertical error bars indicating variability. A diagonal dotted line represents perfect calibration where predicted equals actual probability.

Figure 8. Calibration curve of the validation set at 1 Years.

Figure 9
Calibration plot showing predicted probability on the x-axis and actual probability on the y-axis. Points with error bars closely follow the diagonal line, indicating good model calibration.

Figure 9. Calibration curve of the validation set at 3 Years.

Figure 10
Calibration plot showing predicted versus actual probabilities, with data points connected by lines. Error bars are present, indicating uncertainty. The dashed diagonal line represents perfect calibration where predicted and actual probabilities match.

Figure 10. Calibration curve of the validation set at 5 Years.

Figure 11
Kaplan-Meier survival curve showing survival probability over time in months for three groups: T1 in red, T2 in blue, and T3+T4 in green. T1 and T2 have higher survival probabilities than T3+T4. The p-value is less than 0.0001.

Figure 11. KM curves stratified by T staging in the Training Set.

Figure 12
Kaplan-Meier survival curves show four groups (N0 in red, N1 in blue, N2 in green, N3 in black) over time in months. Survival probability decreases most rapidly for N3 and least for N2. A p-value is less than 0.0001, indicating statistical significance.

Figure 12. KM curves stratified by N staging in the Training Set.

Figure 13
Kaplan-Meier survival curve displays survival probabilities over time for three risk groups: low (red), intermediate (blue), and high (green). The high-risk group maintains the highest survival probability throughout. The p-value is less than 0.0001, indicating significant differences between groups.

Figure 13. KM curves stratified by sTIL expression levels in the Training Set.

Figure 14
Kaplan-Meier survival plot showing survival probability over time in months for three immunohistochemistry (IHC) groups: IHC 0 (red line), IHC 1+ (blue line), and IHC 2+/Fish- (green line). The x-axis represents time in months, and the y-axis shows survival probability. The significant p-value is 0.00064.

Figure 14. KM curves stratified by Her2 expression levels in the Training Set.

4 Discussion

In this large, multicentre retrospective study, we developed and validated a prognostic nomogram specifically tailored for premenopausal women with TNBC receiving anthracycline–taxane chemotherapy. The model-incorporating T stage, N stage, sTIL density, and HER2 immunohistochemistry (IHC) status-demonstrated high predictive accuracy and clinical utility, addressing a critical unmet need in this high-risk population.

Our findings confirm the prognostic relevance of the tumour immune microenvironment in this cohort. Low sTIL density was independently associated with worse disease-free survival, corroborating previous studies that have consistently linked lower sTIL levels with inferior outcomes in patients treated with DNA-damaging chemotherapies (M Geurts et al., 2024; Denkert et al., 2018; Adams et al., 2019). Clinically, these results align with the concept that tumours with low sTILs frequently exhibit an immunologically “cold” phenotype, which may inform future therapeutic strategies aimed at enhancing antitumour immunity in such patients (Sharma et al., 2017; Emens and Middleton, 2015).

We also observed that HER2 IHC 2+/fluorescence in situ hybridization (FISH)-negative status—but not IHC 1+—independently predicted poorer disease-free survival, challenging several earlier reports that associated HER2-zero expression with inferior outcomes in TNBC (Miglietta et al., 2025; Won and Spruck, 2020). This discrepancy may reflect variations in patient selection, treatment setting, or HER2 scoring methodology. Multivariable analysis under the standard HER2-low definition (IHC 1+ or IHC 2+/FISH-negative) indicated that the prognostic risk was predominantly driven by the IHC 2+/FISH-negative subgroup. This suggests that HER2-low TNBC may not represent a homogeneous entity and that IHC 2+/FISH-negative tumours could harbour distinct biological properties—a hypothesis that warrants further validation.

Growing preclinical and clinical evidence indicates that HER2-low breast cancers exhibit altered receptor trafficking, activation of compensatory signalling pathways (e.g., PI3K/AKT, MAPK), and specific tumour–stroma interactions, which may collectively compromise chemotherapy efficacy (Tarantino et al., 2020; Schettini et al., 2021; Modi et al., 2022). Translationaly, the recent efficacy of novel antibody–drug conjugates such as trastuzumab deruxtecan in HER2-low metastatic breast cancer underscores the clinical relevance of this subset and highlights its potential susceptibility to targeted approaches (Banerji et al., 2019; Gennari et al., 2021). Our data suggest that premenopausal patients with HER2-low TNBC—particularly those with IHC 2+/FISH-negative disease—represent a candidate population for adjuvant antibody–drug conjugate trials. Such a strategy could eventually shift the therapeutic paradigm from empirical chemotherapy toward biomarker-guided targeted treatment.

Moreover, TNBC is a molecularly heterogeneous disease, encompassing distinct subtypes such as those described by Lehmann BD et al. (basal-like, immunomodulatory, mesenchymal, etc.) (Lehmann et al., 2011). The prognostic role of sTILs and HER2-low status may vary across these subtypes, a nuance not captured in our retrospective analysis. Future prospective studies integrating transcriptomic subtyping with tumour microenvironment features are warranted to refine personalised prognostic models.

By restricting our analysis to premenopausal women, we minimised confounding from age-related biological differences, such as the higher prevalence of PIK3CA mutations in older patients (Ciriello et al., 2015; Cancer Genome Atlas Network, 2012). This enhances the biological coherence and clinical applicability of our model for this specific demographic. The nomogram’s performance (C-index >0.86) surpasses that of many conventional prognostic tools and genomic assays (Paik et al., 2004; Cardoso et al., 2019), supporting its utility in routine practice for risk-adapted therapeutic decision-making.

From a translational perspective, our model provides a rationale for personalised adjuvant strategies. High-risk patients-those with low sTILs and HER2-low tumours-may benefit from treatment intensification, including capecitabine maintenance (Wang et al., 2021; Masuda et al., 2017) or PARP inhibition in BRCA-mutated cases (J Tutt et al., 2021; Robson et al., 2017). Conversely, low-risk patients could be spared unnecessary toxicity. Furthermore, this tool may facilitate refined patient selection for clinical trials evaluating immunomodulatory agents or novel antibody–drug conjugates in the early-stage setting.

5 Limitations

Several limitations of our study warrant consideration. First, its retrospective design introduces potential selection bias, despite multicentre enrolment. Second, all participants were of Chinese ancestry; external validation in diverse ethnic and geographic populations is essential to confirm generalisability (Eniu et al., 2008). Third, although patients were randomly assigned to training and validation sets, some baseline variables (e.g., family history, postoperative pathological T staging, Ki-67 index, and EGFR status) showed uneven distributions. While our model maintained high predictive performance across cohorts, such imbalances could theoretically introduce selection bias or affect generalizability. Future validation in larger, multi-ethnic cohorts is warranted to confirm the model’s stability. Fourth, although HER2 and sTIL assessments followed international guidelines and underwent central review, HER2 scoring reproducibility remains a challenge in TNBC, and sTIL evaluation is semi-quantitative. Future iterations could benefit from digital pathology, automated image analysis, or standardized immune-based assays to improve objectivity and reproducibility (Klauschen et al., 2018; Vanguri et al., 2022). Fifth, we did not incorporate germline BRCA status, a known determinant of TNBC prognosis and treatment response (Tutt et al., 2010; Tung et al., 2020). Integrating genomic and transcriptomic data in future prospective studies may further refine predictive accuracy and biological insight (Curtis et al., 2012; Pereira et al., 2016).

6 Conclusion

In conclusion, we have identified and validated low stromal TILs and HER2-low status as key biomarkers of intrinsic resistance to anthracycline-taxane chemotherapy in premenopausal TNBC. Beyond prognostication, these factors illuminate distinct biological pathways underlying treatment failure-an immunosuppressive microenvironment and the unique biology of HER2-low tumors. This mechanistic insight provides a compelling rationale for future pharmacological strategies: leveraging immunotherapy to inflame sTIL-low tumors and deploying novel ADCs to target HER2-low subtypes. Our findings thus offer a practical framework for patient stratification and pave the way for more personalized, mechanism-driven adjuvant therapy trials in this challenging disease.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving humans were approved by Medical Ethics Committee of Peking Union Medical College Hospital and Lishui People’s Hospital of Zhejiang Province. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

SC: Writing – review and editing, Writing – original draft. SY: Formal Analysis, Writing – original draft. QZ: Writing – review and editing, Methodology. FK: Formal Analysis, Writing – review and editing. LR: Resources, Writing – review and editing. ZF: Formal Analysis, Writing – review and editing. XZ: Project administration, Funding acquisition, Supervision, Validation, Writing – review and editing. YS: Project administration, Formal Analysis, Supervision, Writing – review and editing. JL: Project administration, Formal Analysis, Supervision, Validation, Writing – review and editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This study was supported in part by grants from the Doctoral Research Initiation Fund (No. 2025bs001) from Lishui Hospital of Wenzhou Medical University, the First Affiliated Hospital of Lishui University, and Lishui People’s Hospital, Zhejiang Province, Zhejiang Provincial Medical and Health Science and Technology Plan Project (No. 2022KY1439), Green Valley Medical Talents Cultivation Project in Lishui, Zhejiang Province, and the Natural Science Foundation of Fujian Province (2025J01080).

Acknowledgements

We thank all the authors listed for their contributions to this study.

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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.

Footnotes

Abbreviations:TNBC, Triple-Negative Breast Cancer; DFS, Disease-Free Survival; sTILs, Stromal Tumor-Infiltrating Lymphocytes; IHC, Immunohistochemistry; HER2, Human Epidermal Growth Factor Receptor 2; FISH, Fluorescence In Situ Hybridization; ER, Estrogen Receptor; PR, Progesterone Receptor; ROC, Receiver Operating Characteristic; AUC, Area Under the ROC Curve; C-index, Concordance Index; HR, Hazard Ratio; CI, Confidence Interval; ADC, :Antibody-Drug Conjugate; BRCA, Breast Cancer gene; PIK3CA, Phosphatidylinositol-4,5-Bisphosphate 3-Kinase Catalytic Subunit Alpha; PARP, Poly(ADP-ribose) Polymerase.

References

Adams, S., Gatti-Mays, M. E., Kalinsky, K., Korde, L. A., Sharon, E., Amiri-Kordestani, L., et al. (2019). Current landscape of immunotherapy in breast cancer: a review. JAMA Oncol. 5 (8), 1205–1214. doi:10.1001/jamaoncol.2018.7147

PubMed Abstract | CrossRef Full Text | Google Scholar

Agostinetto, E., Rediti, M., Fimereli, D., Debien, V., Piccart, M., Aftimos, P., et al. (2021). HER2-Low breast cancer: molecular characteristics and prognosis. Cancers (Basel) 13 (11), 2824. doi:10.3390/cancers13112824

PubMed Abstract | CrossRef Full Text | Google Scholar

Azim, Jr H. A., Michiels, S., Bedard, P. L., Singhal, S. K., Criscitiello, C., Ignatiadis, M., et al. (2012). Elucidating prognosis and biology of breast cancer arising in young women using gene expression profiling. Clin. Cancer Res. 18 (5), 1341–1351. doi:10.1158/1078-0432.CCR-11-2599

PubMed Abstract | CrossRef Full Text | Google Scholar

Banerji, U., Van Herpen, C. M. L., Saura, C., Thistlethwaite, F., Lord, S., Moreno, V., et al. (2019). Trastuzumab duocarmazine in locally advanced and metastatic solid tumours and HER2-expressing breast cancer: a phase 1 dose-escalation and dose-expansion study. Lancet Oncol. 20 (8), 1124–1135. doi:10.1016/S1470-2045(19)30328-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Bareche, Y., Venet, D., Ignatiadis, M., Aftimos, P., Piccart, M., Rothe, F., et al. (2018). Unravelling triple-negative breast cancer molecular heterogeneity using an integrative multiomic analysis. Ann. Oncol. 29 (4), 895–902. doi:10.1093/annonc/mdy024

PubMed Abstract | CrossRef Full Text | Google Scholar

Bianchini, G., Balko, J. M., Mayer, I. A., Sanders, M. E., and Gianni, L. (2016). Triple-negative breast cancer: challenges and opportunities of a heterogeneous disease. Nat. Rev. Clin. Oncol. 13 (11), 674–690. doi:10.1038/nrclinonc.2016.66

PubMed Abstract | CrossRef Full Text | Google Scholar

Bray, F., Laversanne, M., Sung, H., Ferlay, J., Siegel, R. L., Soerjomataram, I., et al. (2024). Global cancer statistics 2022: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 74 (3), 229–263. doi:10.3322/caac.21834

PubMed Abstract | CrossRef Full Text | Google Scholar

Cancer Genome Atlas Network (2012). Comprehensive molecular portraits of human breast tumours. Nature 490 (7418), 61–70. doi:10.1038/nature11412

PubMed Abstract | CrossRef Full Text | Google Scholar

Cardoso, F., Kyriakides, S., Ohno, S., Penault-Llorca, F., Poortmans, P., Rubio, I. T., et al. (2019). Early breast cancer: ESMO clinical Practice Guidelines for diagnosis, treatment and follow-up. Ann. Oncol. 30 (8), 1194–1220. doi:10.1093/annonc/mdz173

PubMed Abstract | CrossRef Full Text | Google Scholar

Ciriello, G., Gatza, M. L., Beck, A. H., Wilkerson, M. D., Rhie, S. K., Pastore, A., et al. (2015). Comprehensive molecular portraits of invasive lobular breast cancer. Cell 163 (2), 506–519. doi:10.1016/j.cell.2015.09.033

PubMed Abstract | CrossRef Full Text | Google Scholar

Curtis, C., Shah, S. P., Chin, S. F., Turashvili, G., Rueda, O. M., Dunning, M. J., et al. (2012). The genomic and transcriptomic architecture of 2,000 breast tumours reveals novel subgroups. Nature 486 (7403), 346–352. doi:10.1038/nature10983

PubMed Abstract | CrossRef Full Text | Google Scholar

Denkert, C., Minckwitz, G. V., Darb-Esfahani, S., Lederer, B., Heppner, B. I., Weber, K. E., et al. (2018). Tumour-infiltrating lymphocytes and prognosis in different subtypes of breast cancer: a pooled analysis of 3771 patients treated with neoadjuvant therapy. Lancet Oncol. 19 (1), 40–50. doi:10.1016/S1470-2045(17)30904-X

PubMed Abstract | CrossRef Full Text | Google Scholar

Emens, L. A., and Middleton, G. (2015). The interplay of immunotherapy and chemotherapy: harnessing potential synergies. Cancer Immunol. Res. 3 (5), 436–443. doi:10.1158/2326-6066.CIR-15-0064

PubMed Abstract | CrossRef Full Text | Google Scholar

Eniu, A., Carlson, R. W., El Saghir, N. S., Bines, J., Bese, N. S., Vorobiof, D., et al. (2008). Guideline implementation for breast healthcare in low- and middle-income countries: treatment resource allocation. Cancer 113 (8 Suppl. l), 2269–2281. doi:10.1002/cncr.23843

PubMed Abstract | CrossRef Full Text | Google Scholar

Foulkes, W. D., Smith, I. E., and Reis-Filho, J. S. (2010). Triple-negative breast cancer. N. Engl. J. Med. 363 (20), 1938–1948. doi:10.1056/NEJMra1001389

PubMed Abstract | CrossRef Full Text | Google Scholar

Gennari, A., André, F., Barrios, C. H., Cortés, J., Azambuja, E. D., DeMichele, A., et al. (2021). ESMO Clinical Practice Guideline for the diagnosis, staging and treatment of patients with metastatic breast cancer. Ann. Oncol. 32 (12), 1475–1495. doi:10.1016/j.annonc.2021.09.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Gonzalez-Angulo, A. M., Timms, K. M., Liu, S. Y., Chen, H. Q., Litton, J. K., Potter, J., et al. (2011). Incidence and outcome of BRCA mutations in unselected patients with triple receptor-negative breast cancer. Clin. Cancer Res. 17 (5), 1082–1089. doi:10.1158/1078-0432.CCR-10-2560

PubMed Abstract | CrossRef Full Text | Google Scholar

Gottesman, M. M. (2002). Mechanisms of cancer drug resistance. Annu. Rev. Med. 53, 615–627. doi:10.1146/annurev.med.53.082901.103929

PubMed Abstract | CrossRef Full Text | Google Scholar

Holohan, C., Van Schaeybroeck, S., Longley, D. B., and Johnston, P. G. (2013). Cancer drug resistance: an evolving paradigm. Nat. Rev. Cancer 13 (10), 714–726. doi:10.1038/nrc3599

PubMed Abstract | CrossRef Full Text | Google Scholar

J Tutt, A. N., Garber, J. E., Kaufman, B., Viale, G., Fumagalli, D., Rastogi, P., et al. (2021). Adjuvant olaparib for patients with BRCA1-or BRCA2-Mutated breast cancer. N. Engl. J. Med. 384 (25), 2394–2405. doi:10.1056/NEJMoa2105215

PubMed Abstract | CrossRef Full Text | Google Scholar

Klauschen, F., Müller, K. R., Binder, A., Bockmayr, M., Hägele, M., Seegerer, P., et al. (2018). Scoring of tumor-infiltrating lymphocytes: from visual estimation to machine learning. Semin. Cancer Biol. 52 (Pt 2), 151–157. doi:10.1016/j.semcancer.2018.07.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Lehmann, B. D., Bauer, J. A., Chen, X., Sanders, M. E., Chakravarthy, A. B., Yu, S., et al. (2011). Identification of human triple-negative breast cancer subtypes and preclinical models for selection of targeted therapies. J. Clin. Invest 121 (7), 2750–2767. doi:10.1172/JCI45014

PubMed Abstract | CrossRef Full Text | Google Scholar

Liedtke, C., Mazouni, C., Hess, K. R., André, F., Tordai, A., Mejia, J. A., et al. (2008). Response to neoadjuvant therapy and long-term survival in patients with triple-negative breast cancer. J. Clin. Oncol. 26 (8), 1275–1281. doi:10.1200/JCO.2007.14.4147

PubMed Abstract | CrossRef Full Text | Google Scholar

M Geurts, V. C., Balduzzi, S., Steenbruggen, T. G., Linn, S. C., Siesling, S., Badve, S. S., et al. (2024). Tumor-Infiltrating lymphocytes in patients with stage I triple-negative breast cancer untreated with chemotherapy. JAMA Oncol. 10 (8), 1077–1086. doi:10.1001/jamaoncol.2024.1917

PubMed Abstract | CrossRef Full Text | Google Scholar

Masuda, N., Lee, S. J., Ohtani, S., Im, Y. H., Lee, E. S., Yokota, I., et al. (2017). Adjuvant capecitabine for breast cancer after preoperative chemotherapy. N. Engl. J. Med. 376 (22), 2147–2159. doi:10.1056/NEJMoa1612645

PubMed Abstract | CrossRef Full Text | Google Scholar

Miglietta, F., Collesei, A., Vernieri, C., Giarratano, T., Giorgi, C. A., Girardi, F., et al. (2025). Development of two machine learning models to predict conversion from primary HER2-0 breast cancer to HER2-low metastases: a proof-of-concept study. ESMO Open 10 (1), 104087. doi:10.1016/j.esmoop.2024.104087

PubMed Abstract | CrossRef Full Text | Google Scholar

Modi, S., Jacot, W., Yamashita, T., Sohn, J., Vidal, M., Tokunaga, E., et al. (2022). Trastuzumab deruxtecan in previously treated HER2-Low advanced breast cancer. N. Engl. J. Med. 387 (1), 9–20. doi:10.1056/NEJMoa2203690

PubMed Abstract | CrossRef Full Text | Google Scholar

Modi, S., Jacot, W., Iwata, H., Park, Y. H., Losada, M. V., Li, W., et al. (2025). Trastuzumab deruxtecan in HER2-low metastatic breast cancer: long-term survival analysis of the randomized, phase 3 DESTINY-Breast04 trial. Nat. Med. 31 (12), 4205–4213. doi:10.1038/s41591-025-03981-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Paik, S., Shak, S., Tang, G., Kim, C., Baker, J., Cronin, M., et al. (2004). A multigene assay to predict recurrence of tamoxifen-treated, node-negative breast cancer. N. Engl. J. Med. 351 (27), 2817–2826. doi:10.1056/NEJMoa041588

PubMed Abstract | CrossRef Full Text | Google Scholar

Parker, J. S., Mullins, M., U Cheang, M. C., Leung, S., Voduc, D., Vickery, T., et al. (2009). Supervised risk predictor of breast cancer based on intrinsic subtypes. J. Clin. Oncol. 27 (8), 1160–1167. doi:10.1200/JCO.2008.18.1370

PubMed Abstract | CrossRef Full Text | Google Scholar

Pereira, B., Chin, S. F., Rueda, O. M., Moen Vollan, H. K., Provenzano, E., Bardwell, H. A., et al. (2016). The somatic mutation profiles of 2,433 breast cancers refines their genomic and transcriptomic landscapes. Nat. Commun. 7 (7), 11479. doi:10.1038/ncomms11479

PubMed Abstract | CrossRef Full Text | Google Scholar

Robson, M., Im, S. A., Senkus, E., Xu, B. H., Domchek, S. M., Masuda, N., et al. (2017). Olaparib for metastatic breast cancer in patients with a germline BRCA mutation. N. Engl. J. Med. 377 (6), 523–533. doi:10.1056/NEJMoa1706450

PubMed Abstract | CrossRef Full Text | Google Scholar

Salgado, R., Denkert, C., Demaria, S., Sirtaine, N., Klauschen, F., Pruneri, G., et al. (2015). The evaluation of tumor-infiltrating lymphocytes (TILs) in breast cancer: recommendations by an International TILs working Group 2014. Ann. Oncol. 26 (2), 259–271. doi:10.1093/annonc/mdu450

PubMed Abstract | CrossRef Full Text | Google Scholar

Saura, C., Cortés, J., Modi, S., Kim, S. B., Hamilton, E., Hurvitz, S. A., et al. (2025). Pooled analysis by best confirmed response to trastuzumab deruxtecan and related biomarkers in patients with HER2-positive metastatic breast cancer from DESTINY-Breast01, DESTINY-Breast02, and DESTINY-Breast03. Ann. Oncol. 17. doi:10.1016/j.annonc.2025.11.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Schettini, F., Chic, N., Brasó-Maristany, F., Paré, L., Pascual, T., Conte, B., et al. (2021). Clinical, pathological, and PAM50 gene expression features of HER2-low breast cancer. NPJ Breast Cancer 7 (1), 1. doi:10.1038/s41523-020-00208-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Schettini, F., Nucera, S., Pascual, T., Martínez-Sáez, O., Sánchez-Bayona, R., Conte, B., et al. (2025). Efficacy and safety of antibody-drug conjugates in pretreated HER2-low metastatic breast cancer: a systematic review and network meta-analysis. Cancer Treat. Rev. 132, 102865. doi:10.1016/j.ctrv.2024.102865

PubMed Abstract | CrossRef Full Text | Google Scholar

Sharma, P., Hu-Lieskovan, S., Wargo, J. A., and Ribas, A. (2017). Primary, adaptive, and acquired resistance to cancer immunotherapy. Cell 168 (4), 707–723. doi:10.1016/j.cell.2017.01.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Shibue, T., and Weinberg, R. A. (2017). EMT, CSCs, and drug resistance: the mechanistic link and clinical implications. Nat. Rev. Clin. Oncol. 14 (10), 611–629. doi:10.1038/nrclinonc.2017.44

PubMed Abstract | CrossRef Full Text | Google Scholar

Sparano, J. A., Gray, R. J., Makower, D. F., Pritchard, K. I., Albain, K. S., Hayes, D. F., et al. (2018). Adjuvant chemotherapy guided by a 21-Gene expression assay in breast cancer. N. Engl. J. Med. 379 (2), 111–121. doi:10.1056/NEJMoa1804710

PubMed Abstract | CrossRef Full Text | Google Scholar

Tarantino, P., Hamilton, E., Tolaney, S. M., Cortes, J., Morganti, S., Ferraro, E., et al. (2020). HER2-Low breast cancer: pathological and clinical landscape. J. Clin. Oncol. 38 (17), 1951–1962. doi:10.1200/JCO.19.02488

PubMed Abstract | CrossRef Full Text | Google Scholar

Tarantino, P., Viale, G., Press, M. F., Hu, X., Penault-Llorca, F., Bardia, A., et al. (2023). ESMO expert consensus statements (ECS) on the definition, diagnosis, and management of HER2-low breast cancer. Ann. Oncol. 34 (8), 645–659. doi:10.1016/j.annonc.2023.05.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Tung, N. M., Boughey, J. C., Pierce, L. J., Robson, M. E., Bedrosian, I., Dietz, J. R., et al. (2020). Management of hereditary breast cancer: American Society of clinical oncology, American Society for radiation oncology, and Society of Surgical Oncology Guideline. J. Clin. Oncol. 38 (18), 2080–2106. doi:10.1200/JCO.20.00299

PubMed Abstract | CrossRef Full Text | Google Scholar

Tutt, A., Robson, M., Garber, J. E., Domchek, S. M., Audeh, M. W., Weitzel, J. N., et al. (2010). Oral poly(ADP-ribose) polymerase inhibitor olaparib in patients with BRCA1 or BRCA2 mutations and advanced breast cancer: a proof-of-concept trial. Lancet 376 (9737), 235–244. doi:10.1016/S0140-6736(10)60892-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Vanguri, R. S., Luo, J., Aukerman, A. T., Egger, J. V., Fong, C. J., Horvat, N., et al. (2022). Multimodal integration of radiology, pathology and genomics for prediction of response to PD-(L)1 blockade in patients with non-small cell lung cancer. Nat. Cancer 3 (10), 1151–1164. doi:10.1038/s43018-022-00416-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Wahba, H. A., and El-Hadaad, H. A. (2015). Current approaches in treatment of triple-negative breast cancer. Cancer Biol. Med. 12 (2), 106–116. doi:10.7497/j.issn.2095-3941.2015.0030

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Wang, S. S., Huang, H., Cai, L., Zhao, L., Peng, R. J., et al. (2021). Effect of capecitabine maintenance therapy using lower dosage and higher frequency vs observation on disease-free survival among patients with early-stage triple-negative breast cancer who had received standard treatment: the SYSUCC-001 randomized clinical trial. JAMA 325 (1), 50–58. doi:10.1001/jama.2020.23370

PubMed Abstract | CrossRef Full Text | Google Scholar

Wolff, A. C., Somerfield, M. R., Dowsett, M., H Hammond, M. E., Hayes, D. F., McShane, L. M., et al. (2023). Human epidermal growth factor receptor 2 testing in breast cancer: ASCO-College of American Pathologists Guideline update. J. Clin. Oncol. 41 (22), 3867–3872. doi:10.1200/JCO.22.02864

PubMed Abstract | CrossRef Full Text | Google Scholar

Won, K. A., and Spruck, C. (2020). Triple-negative breast cancer therapy: current and future perspectives. Int. J. Oncol. 57 (6), 1245–1261. doi:10.3892/ijo.2020.5135

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: anthracycline-taxanechemotherapy, biomarkers, chemoresistance, HER2-low, premenopausal triple-negative breast cancer, stromal tumor infiltrating lymphocytes (sTILs)

Citation: Cai S, Yu S, Zhou Q, Kuang F, Rao L, Fang Z, Zheng X, Shi Y and Li J (2026) Tumour microenvironment and HER2-low status dictate response and resistance to anthracycline-taxane chemotherapy in premenopausal TNBC: a retrospective multicohort study. Front. Pharmacol. 17:1752328. doi: 10.3389/fphar.2026.1752328

Received: 23 November 2025; Accepted: 13 January 2026;
Published: 09 February 2026.

Edited by:

Zhendong Shi, Tianjin Medical University Cancer Institute and Hospital, China

Reviewed by:

Tanzeel Huma, Wake Forest University, United States
Xuliren Wang, Fudan University, China

Copyright © 2026 Cai, Yu, Zhou, Kuang, Rao, Fang, Zheng, Shi and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Xiaoxin Zheng, NzI5NjA3NDI1QHFxLmNvbQ==; Yong Shi, c2hpeW9uZzExMUAxNjMuY29t; Jingdan Li, bGlqaW5nZGFuQGZveG1pbC5jb20=

These authors have contributed equally to this work and share first authorship

These authors have contributed equally to this work

Disclaimer: 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.