Skip to main content

Drug-resistant profiles of extracellular vesicles predict therapeutic response in TNBC patients receiving neoadjuvant chemotherapy

Abstract

Background

Predicting tumor responses to neoadjuvant chemotherapy (NAC) is critical for evaluating prognosis and designing treatment strategies for patients with breast cancer; however, there are no reliable biomarkers that can effectively assess tumor responses. Therefore, we aimed to evaluate the clinical feasibility of using extracellular vesicles (EVs) to predict tumor response after NAC.

Methods

Drug-resistant triple-negative breast cancer (TNBC) cell lines were successfully established, which developed specific morphologies and rapidly growing features. To detect resistance to chemotherapeutic drugs, EVs were isolated from cultured cells and plasma samples collected post-NAC from 36 patients with breast cancer.

Results

Among the differentially expressed gene profiles between parental and drug-resistant cell lines, drug efflux transporters such as MDR1, MRP1, and BCRP were highly expressed in resistant cell lines. Drug efflux transporters have been identified not only in cell lines but also in EVs released from parental cells using immunoaffinity-based EV isolation. The expression of drug resistance markers in EVs was relatively high in patients with residual disease compared to those with a pathological complete response.

Conclusions

The optimal combination of drug-resistant EV markers was significantly efficient in predicting resistance to NAC with 81.82% sensitivity and 92.86% specificity.

Peer Review reports

Background

Recent advancements in neoadjuvant chemotherapy (NAC) have broadened its applications beyond large aggressive tumors in breast cancer treatment. Candidates for preoperative systemic therapy are inoperable breast cancer, and operable breast cancer in selected patients who desire breast conservation with large primary tumor relative to breast size, who have HER2-positive disease and triple-negative breast cancer greater than clinical T2 or clinical N1. Current NAC strategies are highly personalized, considering the cancer subtype, stage, and molecular characteristics, to optimize treatment efficacy. Preferred regimens of NAC for HER2-negative breast cancer are dose-dense adriamycin and cytoxan followed by paclitaxel or weekly paclitaxel, and those for HER2-positive breast cancer are paclitaxel + trastuzumab, paclitaxel + carboplatin + trastuzumab (TCH), and TCH + pertuzumab. Moreover, these approaches might be integrated with both conventional therapies and immunotherapies to align with each patient’s unique cancer profile [1, 2]. The response to NAC remains a key prognostic indicator in several studies. Moreover, these developments in NAC protocols have notably increased the feasibility of breast-conserving surgeries, marking a significant shift in the breast cancer management landscape [3,4,5].

Clinically established therapeutics are highly effective in treating breast cancer; however, they present significant challenges. One major challenge is the limited therapeutic options for triple-negative breast cancer (TNBC), which is often resistant to chemotherapy and cross-resistant to other antitumor agents, suggesting multidrug resistance. In addition, the development of drug resistance to anthracycline + cyclophosphamide + taxane (ACT) combination chemotherapies undermines their therapeutic potential [6]. Moreover, some tumors relapse rapidly after NAC and surgery. Thus, a better understanding of the molecular mechanisms underlying drug resistance is urgently required to effectively treat TNBC.

Extracellular vesicles (EVs) are mediators of cell-to-cell communication, are surrounded by lipid bilayers, and are released from living cells. EVs carry real-time molecular information about their cell of origin, such as nucleic acids, proteins, and lipids. Among body fluid-derived EVs, tumor-derived EVs comprise vesicles released by highly heterogeneous breast cancer cells [7]. Therefore, the more accurately tumor-derived EVs can be isolated, the more accurately they can reflect the pathophysiological characteristics and behaviors of tumor cells [8]. An analysis of tumor-derived EVs collected from tumor tissues, especially at the early stages of drug resistance development, may offer insights into the possibility of screening and monitoring cancers, including breast cancer, and provide appropriate treatment options for patients in terms of precision medicine [9].

In this study, we developed several stable drug-resistant TNBC cell lines in vitro using a stepwise treatment strategy with chemotherapeutic agents for 28 weeks. We aimed to identify highly expressed biomarkers within EVs released from these cell lines and evaluate the clinical feasibility of EV-based assessments for predicting drug response in breast cancer patients receiving NAC.

Methods

Cell lines and anticancer drugs

The human TNBC cell lines HCC1395 (CRL-2324), MDA-MB-231 (HTB-26), and MDA-MB-468 (HTB-132) were purchased from the American Type Culture Collection (ATCC, Manassas, VA, USA). All cell lines were grown in RPMI-1640 (22400-089, Gibco, Carlsbad, CA, USA) supplemented with 10% heat-inactivated fetal bovine serum (FBS, 12483-020, Gibco), and 1% penicillin-streptomycin (15140-122, Gibco), and maintained in a humidified atmosphere of 5% CO2 in air at 37 °C.

Doxorubicin hydrochloride (Adriamycin, or Anthracycline chemotherapy drug, D4000) and docetaxel (Taxotere, or Taxane chemotherapy drug, D1000) were purchased from LC Laboratories (Woburn, MA, USA), and cyclophosphamide monohydrate (Cytoxan, NSC-26271) was purchased from Selleck Chemicals (Houston, TX, USA).

Induction of chemotherapy resistance in breast cancer cells

Drug-resistant sublines of each TNBC cell line were derived from each original parental cell line by continuous exposure of low to high doses of anticancer drugs (1/120 IC50, 1/90 IC50, 1/60 IC50, 1/30 IC50, 1/10 IC50, and IC50) for over 6 months. Each parental cell line was treated with anthracycline + cyclophosphamide (AC, 1:10 molar ratio) for 72 h, repeated four times. Subsequent treatments with Taxotere (T) proceeded in the same manner, and this process was defined as one cycle. The cells were maintained for 6 months while the drug concentration was gradually increased and allowed to recover for an additional month.

Evaluation of drug-resistant activity and growth rate

Sensitivity to chemotherapeutic drugs was measured with a colorimetric assay using MTT (M2003, Sigma-Aldrich, St. Louis, MO, USA). Following the treatment of cells with serial dilutions of AC or T for 72 h, MTT was added to each well and incubated for 4 h at 37 °C. Methanol was then added to each well and mixed for 30 min on an orbital shaker. The absorbance was recorded at 570 nm with a correction wavelength of 690 nm using a NanoDrop 3000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). IC50 values were calculated using Prism v9.0.0 (GraphPad, San Diego, CA, USA).

To measure the growth rate of the derived cell lines, three groups of drug-resistant sublines and their parental cell line counterparts (1 × 105) were seeded in 6-well plates and allowed to attach to the well surface. After attachment, the cells were counted every 24 h for 96 h using a cell counter (TC20, Bio-Rad, Hercules, CA, USA). The trypan blue exclusion assay was used to determine the number of viable cells in the cell suspension and evaluate the doubling time of the cells.

RNA sequence analysis

The concentration and quality of total RNA were checked using a Qubit 2.0 fluorometer (Thermo Fisher Scientific). Total RNA (10 ng) was used to prepare strand-specific barcoded RNA libraries using the Ion AmpliSeq Transcriptome Human Gene Expression Kit (Thermo Fisher Scientific) following the manufacturer’s protocol. The Ion AmpliSeq Transcriptome Human Gene Expression Kit was designed for simultaneous targeted amplification of over 20,000 human genes using a single primer pool. A short amplicon (approximately 110 bp) was obtained from each target gene. AmpliSeq sequencing data were obtained using the Torrent Mapping Alignment Program optimized for Ion Torrent sequencing data to align raw sequencing reads against a custom reference sequence set containing all transcripts targeted by the AmpliSeq kit [10, 11].

Transcriptome analysis

Differentially expressed genes (DEGs) were identified using DEGSeq (version 1.48.0) with p-values of < 0.05 and|fold changes|>2.5 threshold [12]. The function of each DEG was annotated based on the Biological Process Gene Ontology gene set (MSigDB collections, C5, BROAD Institute, Cambridge, MA, USA) and the Gene Reference into Function (GeneRIF, The National Center for Biotechnology Information) database. Volcano plots, bar plots, Venn diagrams, and heat maps were generated using ggplot2 (version 3.3.5) and the Complex Heatmap (version 2.10.0) software. All statistical analyses and visualizations were performed using R (version 4.1.2) and R Studio environment (release 077589bc).

Flow cytometry for EV surface profiling

The surface profiles of cells and EVs were analyzed using flow cytometry (FACS LSR Fortessa system, Becton Dickinson, Franklin Lakes, NJ, USA) and specific antibodies against surface proteins. Cell samples were incubated for 30 min at 4 °C in the dark with one test dose of drug resistance detection antibodies (anti-MDR1, anti-MRP1, and anti-BCRP) and rinsed twice with FACS buffer to prevent excessive reactions; fluorescein isothiocyanate (FITC)-labeled secondary fluorescent antibodies were used to detect fluorescence signals.

Breast cancer-derived EVs bound to 3-µm microbeads (SPHERO Streptavidin Coated Particles, SVP 30 − 5, Spherotech Inc., Lake Forest, IL, USA) that were conjugated with biotinylated breast cancer-targeting antibodies (anti-EpCAM, anti-ITGA2, and anti-ITGAV), during incubation to enable the isolation of breast cancer-derived EVs, following protocols in a previous study [13]. After isolation from tumor tissue, the drug resistance detection antibody was then immobilized with a fluorescent detection antibody. Detailed schematics are presented in Additional File 1: Figure S1a and S1b.

Confocal microscopy

To stain actin filaments, Alexa Fluor 594-conjugated phalloidin (A12381, Thermo Fisher Scientific) was used according to the manufacturer’s protocol, and DAPI (Vectashield, H-1200, Vector Laboratories, Inc., Newark, CA, USA) was used to stain cell nuclei. To confirm the presence of EVs, the captured EVs were detected using a primary PE-Cy7-labeled antibody against the general EV marker, CD63. EVs were also immobilized with a primary drug resistance detection antibody against MDR1 and a secondary FITC-labeled fluorescent antibody. Fluorescence images were obtained using a Zeiss LSM 700 confocal microscope (Carl Zeiss). The target and detection antibodies used in this study are listed in Additional File 1: Table S1.

Transmission electron microscopy (TEM)

A drop of the EV sample was fixed in 2% glutaraldehyde-2% paraformaldehyde and placed on a Formvar carbon-coated grid for 15 s. Droplets were removed using filter paper, and a drop of 1% uranyl acetate was added for 15 s for negative staining, removed using filter paper, and washed with a drop of distilled water. The dried grids were observed using a transmission electron microscope (JEM-1011, JEOL, Tokyo, Japan) at an acceleration voltage of 80 kV equipped with a MegaView III CCD camera.

Clinical characteristics of the participants

Clinical samples were obtained from subjects who visited Severance Hospital in South Korea in accordance with the guidelines of the independent Ethics Committee at the College of Medicine Yonsei University (IRB No. 4-2020-1292). Informed consent for the use of blood samples for research purposes was obtained from all patients. Pre-operative plasma samples were collected from the same patient before anesthesia. The criteria for subject eligibility in the analysis included (1) a confirmed pathological diagnosis of breast cancer, (2) collection of blood samples post-NAC and during the pre-operative period, and (3) hemolysis assessed before the isolation of EV to evaluate plasma sample quality. Details of the 36 individuals are shown in Additional File 1: Table S2.

Receiver operating characteristic (ROC) analysis

ROC analysis of drug-resistant EV markers was performed on data from 20 patients with TNBC using MedCalc version 20.014 (MedCalc Software Ltd., Ostend, Belgium). We used univariate ROC analysis for each marker to obtain the area under the ROC curve (AUC) and evaluate the diagnostic power of drug-resistant EV marker combinations. Optimal criterion values were calculated by considering not only sensitivity and specificity, but also disease prevalence and costs of various decisions [14]. After performing a univariate ROC analysis on each combination of drug-resistant EV markers, we chose the “Best” combination with the highest AUC and the lowest standard error of AUC. Based on these calculations, we developed a combinatorial predictive score composed of MDR1, MRP1, and BCRP (combi-3) that was utilized to predict drug response in patients with breast cancer. Statistical analyses were performed using a one-way analysis of variance or Welch’s t-test.

Results

Establishment of TNBC cell lines with a drug-resistant phenotype

To generate chemotherapy-resistant variants of each TNBC cell line, the parental cell line was exposed to gradually increasing concentrations of chemotherapeutic drugs for 4 weeks per treatment cycle (Fig. 1a). Approximately 28 weeks and 56 passages were required for three cell lines (HCC1395, MDA-MB-231, and MDA-MB-468) to acquire stable drug resistance. The generation of drug resistance in TNBC cell lines was validated by morphological observations, dose-response analyses, and gene expression patterns at the transcriptional level. We generated dose-response curves and calculated IC50 values for both parental and drug-resistant (derived) sublines to determine whether the derived sublines acquired resistance to chemotherapy (Fig. 1b). As expected, the IC50 concentrations for drug-resistant sublines in response to AC or T treatment were significantly higher than those for the wild-type cell lines. Furthermore, doubling times were shorter in resistant cells than in wild-type cells, implying that the difference in proliferative capacity might be associated with the survival advantage of drug-resistant cells under chemotherapy. In summary, long-term exposure of TNBC cell lines to chemotherapeutic drugs was accompanied by distinct morphological changes, from 2.9- to 29.7-fold higher IC50 values and from 0.8- to 0.91-fold shorter doubling times in resistant TNBC sublines, thereby confirming the successful acquisition of drug-resistant phenotypes.

Fig. 1
figure 1

(a) Schematic illustration of the drug treatment process for the generation of drug-resistant cell lines. (b) Dose-response curves and IC50 values. Cell growth curves and their associated doubling times were measured. Data are expressed as Mean ± SEM (n = 5)

During drug treatment, we observed a gradual change in the shape of the TNBC cell lines. HCC1395, MDA-MB-231, and MDA-MB-468 cells spread more extensively in adherent monolayer cultures after acquiring resistance to anticancer drugs. The significance of the morphological changes was confirmed by the visualization of nuclei and actin filaments. We observed that the morphology of drug-resistant cells was distinct from that of parental cells of the same cellular age (Fig. 2a). Next, when we measured the area of TNBC cell lines, all resistant sublines exhibited enlarged nucleus/cytoplasm in monolayer proliferation compared to the parental cells (Fig. 2b).

Fig. 2
figure 2

(a) Morphology of drug-resistant cell lines. Cell nuclei were stained with DAPI, and actin filaments were stained with phalloidin-Alexa Fluor 594. (b) Quantification of the area of TNBC cell lines. The cell areas were measured from three distinct cells in the image. Data are expressed as mean ± SD (n = 3). Significant differences between groups were determined using the Welch’s t-test (*p < 0.05 and **p < 0.01)

DEGs in resistant TNBC cell lines

To identify the differential abundance of transcripts in resistant TNBC, we performed RNA sequencing and transcriptomic analysis of the wild-type cell lines (hereafter referred to as HCC1395/WT, MDA-MB-231/WT, and MDA-MB-468/WT) and drug-resistant sublines (hereafter referred to as HCC1395/R, MDA-MB-231/R, and MDA-MB-468/R). A comparison of drug-resistant sublines with wild-type cell lines revealed a large number of significant differences in the volcano plot (Fig. 3a). Based on criteria [adjusted p-value ≤ 0.05, and fold change (FC) ≥ 2.5], 111, 395, and 151 upregulated DEGs, and 373, 613, and 378 downregulated DEGs were identified in HCC1395, MDA-MB-231, and MDA-MB-468, respectively. As the volcano plots indicate, in total, 657 significantly upregulated DEGs and 1,364 significantly downregulated DEGs were identified in the comparison between the wild-type and resistant groups; however, considering the limited expression level of downregulated DEGs not suitable for use as a diagnostic biomarker, they lost their significance in clinical analysis and were excluded from the candidate biomarker.

Fig. 3
figure 3

(a) Volcano plots [-log10(P-value) vs. log2(fold change)] for each comparison of wild-type and resistant HCC1395 (left), MDA-MB-231 (middle), and MDA-MB-468 (right) cells after long-term ACT treatment. The significant downregulated (green) and upregulated (red) transcripts are shown in each volcano plot. (b) A list of the 50 most significant genes in each upregulated DEG following long-term ACT treatment. The function of genes was annotated based on the GeneRIF database. (c) Heat map showing the relative abundance of selected drug resistance-related genes. The predicted locations of proteins encoded by drug resistance-related genes were categorized based on the Human Protein Atlas

To further explore the biological activities of the top 50 upregulated DEGs, we cross-checked drug resistance- and proliferation-related genes found in GeneRIF, a functional annotation database. We identified 14 drug resistance-related DEGs and 34 proliferation-related DEGs, suggesting that these genes play important roles in the acquisition of chemotherapy resistance in TNBC (Fig. 3b). In addition to analysis using GeneRIF, Gene Ontology enrichment analysis for biological processes (GOBP) revealed that all resistant sublines had a high gene set variation analysis enrichment score in at least one GOBP term among drug transmembrane transport, drug transport, and cellular response to drugs. Notably, HCC1395 cells had positively enriched scores for all three GOBP terms, suggesting that these cells firmly acquired drug resistance, followed by MDA-MB-231 and MDA-MB-468 (Additional File 1: Figure S2). We also identified several gene sets potentially involved in drug resistance that contributed to the elevated enrichment scores in HCC1395 cells using gene set enrichment analysis and gene network analysis (Additional File 1: Figure S3). Importantly, many of the annotated genes found in these analyses encode membranous proteins, including drug efflux transporters (e.g., ABCB1 and ABCG2), chemokine receptors (e.g., CXCR4), receptor tyrosine kinases (e.g., EGFR), and matrix metalloproteinases (e.g., MMP1), and so on (Fig. 3c). Among the membrane proteins that play vital roles in multidrug resistance, we focused on drug efflux transporters, the most frequently discovered molecules in chemoresistance [15].

Induction of drug efflux transporters by chemotherapy agent

Drug efflux transporters, also called ATP-binding cassette (ABC) transporter, are composed by 49 ABC genes arranged in eight human subfamilies [16]. Among those subfamilies, the three major types involved in cancer drug resistance are ABCB1, ABCC1, and ABCG2, which represent MDR1, MRP1, and BCRP proteins, respectively. Overexpression of these three ABC transporters can increase the efflux of drugs from cancer cells, thereby reducing the intracellular drug concentration [17]. We first investigated the expression levels of the three main drug efflux transporters in wild-type and resistant TNBC cells. Compared with those in the wild-type sublines, mean fluorescence intensity (MFI) values of MDR1, MRP1, and BCRP in the drug-resistant sublines were significantly higher (Fig. 4a, red vs. green column). Overall, stronger signals by drug efflux transporters were observed in the HCC1395, MDA-MB-231, and MDA-MB-468 drug-resistant sublines. Moreover, significant increases in drug efflux transporter expression were observed in the drug-resistant sublines after 48 h of AC or T treatment. Expression of these transporters was higher for HCC1395 and MDA-MB-231 cells and to a smaller extent for MDA-MB-468 cells in drug-resistant sublines. Long-term drug exposure in TNBC increased the immediate response capacity to drugs (Fig. 4a; green vs. orange or blue columns). These results indicate that transient and long-term exposure to chemotherapeutic agents induces the expression of cell surface drug transporters, which are primarily related to drug resistance.

Fig. 4
figure 4

(a) Detection of drug-resistant proteins on the surface of HCC1395, MDA-MB-231, and MDA-MB-468 cells under different conditions. MFI for each target is shown on the y-axis: untreated control cells (black), wild-type cells (red), stable drug-resistant cells (green), drug-resistant cells after 48 h of AC treatment (orange), and drug-resistant cells after 48 h of T treatment (blue). Data are shown as mean ± SD of three independent experiments. Significant differences between groups were determined using one-way ANOVA (*p < 0.05; **p < 0.01; and ***p < 0.001). (b) TEM images of EVs extracted from each cell. Black scale bars represent 100 nm. (c) A representative confocal image of enriched EVs by immuno-beads with a red box highlighting a SEM image of EVs attached to immuno-beads. White scale bars represent 10 μm, except for the SEM image which represents 200 nm. (d) Flow cytometry gating strategies for drug-resistant EV markers (d) Comparison between wild-type and drug-resistant clones in cell lines and EVs. Significant differences between groups were determined using one-way ANOVA (ns, not significant; *p < 0.05; **p < 0.01; and ***p < 0.001)

Detection of drug efflux transporters on EVs

Next, we hypothesized that EVs transfer their molecular cargo to adjacent cells which potentially influence the drug resistance profiles of tumor cells through the exchange of highly expressed drug efflux transporters. We designed an immunoaffinity-based breast cancer EV isolation method for cancer-derived EV surface profiling to evaluate whether experimentally accessible amounts of drug transporters could be detected in EVs. Before quantitative measurements, we verified that all EVs released from drug-resistant sublines had similar physicochemical characteristics. TEM and nanoparticle tracking analysis (NTA) confirmed the size, morphology, and concentration of EVs based on standard guidelines [18]. According to TEM observations, isolated EVs were approximately 100–200 nm in diameter and mostly spherical (Fig. 4b). When checking NTA results, the average particle sizes of EVs from drug-resistant sublines presented slightly larger particle sizes and higher EV concentrations than did EVs from wild-type cell lines (Additional File 1: Figure S4). Using confocal microscopy, we confirmed that a strong red fluorescent signal of the CD63 EV marker and a green fluorescent signal of the MDR1 drug transporter were detected on the EV attached to the microbeads, which indicates continuous activation at the surface of EVs (Fig. 4c). Scanning electron microscopy (SEM) was used to confirm the presence of tumor cell-derived EVs, along with their morphologies and sizes (Fig. 4c, red box).

Flow cytometric analyses identified three drug efflux transporters (MDR1, MRP1, and BCRP) in both the cell bodies and EVs released from parental cells. The fluorescence intensity of the control beads was used to set gating strategies for MDR1, MRP1, and BCRP. (Fig. 4d). Overall, the resistant clones were characterized by enhanced expression of drug efflux transporters. For MDR1, MRP1, and BCRP, the differences in MFI values between the wild-type and resistant clones significantly increased in EVs compared to the differences in the cell lines. The expression levels of drug efflux transporters in resistant EVs exhibited approximately 2.2 to 3.0-fold higher rates (Fig. 4e). These results suggest that the developed breast cancer-derived EV isolation method allows detection of membrane protein markers on the surface of EVs. Furthermore, EVs reflect the characteristics of the cells they are derived from and may be strongly relevant to inducing drug resistance after chemotherapy.

Correlation analysis between drug-resistant EV markers with tumor response in TNBC patients treated with NAC

In this study, we evaluated the clinical utility of EVs as biomarkers for predicting pathological complete response (pCR) after NAC using drug-resistant EV marker selection and in vitro analyses. Among the 36 subjects who participated, pCR was defined as the absence of residual invasive cancer in the breast and lymph nodes 3–6 months after NAC, as assessed by pathological examination of tumor tissues post-surgery. The participants were divided into the pCR and non-pCR groups. Furthermore, breast cancer subtypes were determined via immunohistochemistry (IHC) of tumor biopsies following the guidelines established by Severance Hospital, which are based on the NCCN guidelines. Luminal A was defined as HR-positive, HER2-negative, with low Ki-67 levels (< 15%), whereas Luminal B was HR-positive, but with HER2-positive or higher Ki-67 levels (≥ 15%). The HER2 subtype was identified as HR-negative but HER2-positive, whereas the TNBC subtype was characterized as negative for both HR and HER2.

We examined the expression levels of MDR1, MRP1, and BCRP on the surface of EVs from the plasma of 20 patients with TNBC who underwent NAC. Among the 20 patients, 10 showed residual tumors on surgical specimens and 10 achieved pCR after NAC. The expression of drug resistance markers in EVs in both groups was quantified objectively using the MFI values obtained by flow cytometry (Fig. 5a). The arbitrary cut-off values for the MFI of MDR1, MRP1, and BCRP were 506.8, 634.9, and 1175.8, respectively. At least one drug-resistant EV marker was highly expressed in all patients except two in the no pCR group, whereas most patients in the pCR group exhibited low expression levels of all drug-resistant EV markers (Additional File 1: Table S3). Notably, patients presenting with recurrence during follow-up or death showed elevated expression levels of all drug-resistant EV markers, implying a clinical correlation between drug-resistant EV markers and patient prognosis, as well as therapeutic response. The MFI values and distributions of all the EV markers differed significantly between the pCR and the no pCR groups (Fig. 5b). To assess the clinical utility of possible combinations of drug-resistant EV markers, we performed a ROC curve analysis on the same patient group and set the AUC cut-off for a good diagnostic value at > 0.8. The AUC values for MDR1, MRP1, and BCRP were 0.82 (p = 0.002), 0.91 (p < 0.001), and 0.93 (p < 0.001), respectively (Fig. 5c).

Fig. 5
figure 5

(a) Distribution of clinical variables and drug-resistant EV markers in patients with TNBC who underwent NAC. (b) Comparison of the MFI values for drug-resistant EV markers between patients with TNBC who responded to NAC and patients with TNBC who were resistant to NAC. Data are expressed as MFI in arbitrary units (a.u.) and represents the mean ± SEM of 20 individuals in each group. Statistical analysis was performed using the Welch’s t-test. (c) ROC curves for the three drug-resistant EV markers

Diagnostic potential of drug-resistant EV markers in all subtypes of patients with breast cancer treated with NAC

To better discriminate the therapeutic response and achieve the desired sensitivity and/or specificity of the markers, each biomarker was combined (Additional File 1: Figure S5). Using logistic regression, we identified the best combination of predictive biomarkers for drug response [19]. The patients with residual tumors after NAC had higher combi-3 expression, which was the optimal EV marker combination (p < 0.001; Fig. 6a, upper). Moreover, there were distinct differences in patients with breast cancer regardless of the subtype (n = 36), suggesting that the application of drug-resistant EV markers is not confined to TNBC but could be applied to all subtypes (Fig. 6a, lower). Next, we evaluated the performance of potential marker combinations in predicting drug responses based on the index using an optimal cut-off value of 0.631. The AUC values for MDR1, MRP1, and BCRP, when applied to all subtypes, were 0.78 (p = 0.002), 0.74 (p = 0.004), and 0.90 (p < 0.001), respectively. The AUC value of combi-3 was 0.92 (95% confidence interval, 0.830–1.000), with an 81.82% sensitivity and a 92.86% specificity (Fig. 6b). Therefore, we propose a more powerful model that combines EV markers to achieve the best diagnostic and prognostic accuracy for predicting treatment response in patients with breast cancer.

Fig. 6
figure 6

(a) Comparison of the combi-3 for drug-resistant EV markers between patients with breast cancer who responded to NAC and patients with breast cancer who were resistant to NAC. The top graph represents patients with TNBC (n = 20) and the lower graph represents patients with all subtypes (n = 36). Statistical analysis was performed using the Welch’s t-test. (b) ROC curves for the three drug-resistant EV markers and the combination marker

Discussion

Within the past few decades, the standard method of evaluating the effectiveness of NAC has been to identify any residual tumor in surgical specimens; however, this is not possible before surgery. Approximately 10% of breast tumors do not respond to NAC, and resistant breast tumors, particularly TNBC, lead to disease progression and poor prognosis. Accordingly, to reduce disease progression, it is necessary to predict which patients will not respond to standard treatments [20]. Circulating EVs in the blood and intracellular communication vesicles with a lipid bilayer ranging in size from 50 to 300 nm are essential for tumorigenesis, development, progression, and metastasis [21]. Multiple technologies have been developed for tumor-derived EV detection and analysis (e.g., immunoaffinity-based capture) and have greatly advanced our understanding of tumor characteristics through liquid biopsy, despite the absence of tumor tissues. EVs can shuttle bioactive molecules such as proteins and a wide variety of genetic materials from one cell to another, leading to molecular transformations in recipient cells [22,23,24,25]. There may be an interplay or synergy between tumor cells in the acquisition of drug resistance via EV exchange. Therefore, we suggest a method that focuses on the role of EVs in multidrug resistance for early therapeutic response prediction and therapy monitoring.

Our study revealed the dynamics of the epigenetic changes that lead to drug resistance after chemotherapy in a cell line model. Various molecular mechanisms are involved in chemoresistance; among these, we found significant changes in cancer stemness in drug-resistant TNBC models. For example, dysregulated TGF-β and Wnt signaling pathways may affect the overall progression to malignancy (Additional File 1: Figure S6). They are also known to play a positive role in promoting drug-resistant properties in the cancer stem cell (CSC) population [26, 27]. Furthermore, the cytotoxicity-induced morphological changes can be related to drug resistance. Cytoskeletal reconstruction induces biological changes in cancer cells with drug resistance [28]. After long-term treatment with chemotherapeutic drugs, actin stress fibers change, with a distinct feature showing the migratory dynamics of cell spreading [29]. All these phenotypic features can be predominantly attributed to their tumorigenic potential and multidrug resistance. Moreover, we focused mainly on the overexpression of ABC transporters, which results in drug resistance, a characteristic feature of CSC. Accumulating evidence from numerous studies, including ours, indicates that high expression of transmembrane proteins of this superfamily, such as MDR1, MRP1, and BCRP, is found in breast cancer, particularly in breast CSCs [30].

Despite the clear relevance of ABC transporters, which play a critical role in the development of multidrug resistance, clinical approaches for assessing these proteins in the development of drug resistance have not been successful and are yet to be elucidated [31, 32]. However, with the availability of the latest technology for EV analysis, we recommend re-evaluating the role of ABC transporters within EVs, not focusing on those in tumor tissues. Two notable singularities of our study are worth discussing. First, we found that a certain percentage of cells were stably resistant to NAC when a chemotherapeutic drug was administered continuously (Fig. 4a). A superficial explanation for this is that drug resistance may be induced by more than simple genetic alterations; dysregulation of major epigenetic factors in breast cancer may play a more important role. This hypothesis is consistent with several reports on the epigenetic control of CSCs and their influence on tumorigenesis, development, and responsiveness to therapy [33, 34]. Second, by using the enrichment of tumor-derived EVs, a higher expression of drug efflux transporters was observed in drug-resistant EVs than in wild-type EVs (Fig. 4e) [35,36,37]. This has been suggested earlier and may have clinical significance, as it indicates that resistant clones may induce bursts of ABC transport in EVs to enhance drug resistance.

Although a few studies have utilized miRNAs, EV concentrations, and cancer antigens to monitor and predict response to treatments [38,39,40], our study is a novel attempt to use tumor-derived EVs to assess ABC transporters isolated from the plasma of patients with breast cancer who are treated with NAC. In a retrospective study using plasma from 36 patients, considerable differences in the expression of MDR1, MRP1, and BCRP were detected between patients with pCR and those with residual tumors. The combination of these three parameters offered acceptable sensitivity, specificity, and accuracy in predicting the effectiveness of NAC through cumulative ROC analysis. Apart from the lack of understanding regarding the role of EVs in cancer stem cells and drug resistance, we suggest that this method is a better strategy for repeatedly screening molecular information to predict therapeutic responses. However, three minor limitations of this study merit further discussion. First, we performed NGS analysis for established drug-resistant cell models compared to the wild-type to identify drug resistance genes and validated their expression in EVs isolated from breast cancer patients undergoing NAC. This was due to our concerns regarding the challenges of accurately separating tumor-derived EVs and ensuring sufficient nucleic acid material for reliable NGS analysis. Further research is required to precisely separate tumor-derived EVs and analyze their contents. Second, there was a shortage of sufficient sample sizes to conduct significant cohort studies with high statistical validity. Another issue is the lack of longitudinal studies analyzing changes in EV markers over time in patients. Future studies will advance to larger patient groups and address how EVs crosstalk with other tumor cells of different phenotypes to obtain drug-resistant characteristics, in addition to MDR1, MRP1, and BCRP.

Conclusions

We established drug-resistant TNBC cell lines and investigated the genes involved in drug transport, such as MDR1, MRP1, and BCRP, which were highly expressed in resistant cell lines compared to that in their wild-type counterparts. The expression of MDR1, MRP1, and BCRP in breast cancer-derived EVs increased after in vitro chemotherapeutic treatment, particularly in drug-resistant cell lines. To investigate the clinical significance of drug-resistant EV markers in patients with TNBC receiving NAC, patients with residual tumors were found to have higher expression levels of MDR1, MRP1, and BCRP in EVs than in patients with a pathological complete response. Integrated analysis of MDR1, MRP1, and BCRP expression showed significant differences according to tumor response, not only in TNBC but also in other subtypes. In conclusion, our findings demonstrate that this EV marker combination could be a useful predictor for discriminating breast cancer patients with residual disease from those with no residual disease after NAC, especially in TNBC.

Data availability

The RNA-seq data generated here is available in the Gene Expression Omnibus database under the accession code GSE230273. All processed data generated in this study are included in the Supplementary Information.

Abbreviations

AC:

Anthracycline + Cyclophosphamide

ACT:

Anthracycline + Cyclophosphamide + Taxane

ATCC:

The American Type Culture Collection

AUC:

Area Under The Curve

ABC:

ATP-binding Cassette

BCRP:

Breast Cancer Resistance Protein

CSC:

Cancer Stem Cell

Combi-3:

Combinatorial Predictive Score Composed of MDR1,MRP1,and BCRP

DEGs:

Differentially Expressed Genes

EVs:

Extracellular Vesicles

FBS:

Fetal Bovine Serum

FC:

Fold Change

GOBP:

Gene Ontology Enrichment Analysis for Biological Processes

GeneRIF:

Gene Reference into Function

IC50:

Half-maximal Inhibitory Concentration

IHC:

Immunohistochemistry

MFI:

Mean Fluorescence Intensity

MRP1:

MDR-associated Protein 1

MDR1:

Multidrug Resistance Protein 1

NTA:

Nanoparticle Tracking Analysis

NAC:

Neoadjuvant Chemotherapy

pCR:

Pathological Complete Response

ROC:

Receiver Operating Characteristic

SEM:

Scanning Electron Microscopy

T:

Taxotere

TCH:

Paclitaxel + Carboplatin + Trastuzumab

TEM:

Transmission Electron Microscope

TNBC:

Triple-negative Breast Cancer

References

  1. Gradishar WJ, Moran MS, Abraham J, Abramson V, Aft R, Agnese D, Allison KH, Anderson B, Burstein HJ, Chew H, et al. NCCN Guidelines® insights: breast Cancer, Version 4.2023. J Natl Compr Canc Netw. 2023;21(6):594–608.

    Article  PubMed  Google Scholar 

  2. Zaheed M, Wilcken N, Willson ML, O’Connell DL, Goodwin A. Sequencing of anthracyclines and taxanes in neoadjuvant and adjuvant therapy for early breast cancer. Cochrane Database Syst Rev. 2019;2:CD012873.

    PubMed  Google Scholar 

  3. Colomer R, Saura C, Sanchez-Rovira P, Pascual T, Rubio IT, Burgues O, Marcos L, Rodriguez CA, Martin M, Lluch A. Neoadjuvant Management of early breast Cancer: a clinical and investigational position Statement. Oncologist. 2019;24(5):603–11.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Spring LM, Fell G, Arfe A, Sharma C, Greenup R, Reynolds KL, Smith BL, Alexander B, Moy B, Isakoff SJ, et al. Pathologic Complete Response after neoadjuvant chemotherapy and impact on breast Cancer recurrence and survival: a Comprehensive Meta-analysis. Clin Cancer Res. 2020;26(12):2838–48.

    Article  PubMed  PubMed Central  Google Scholar 

  5. An SJ, Duchesneau ED, Strassle PD, Reeder-Hayes K, Gallagher KK, Ollila DW, Downs-Canner SM, Spanheimer PM. Pathologic complete response and survival after neoadjuvant chemotherapy in cT1-T2/N0 HER2 + breast cancer. NPJ Breast Cancer. 2022;8(1):65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Luque-Bolivar A, Perez-Mora E, Villegas VE, Rondon-Lagos M. Resistance and overcoming resistance in breast Cancer. Breast Cancer (Dove Med Press). 2020;12:211–29.

    CAS  PubMed  Google Scholar 

  7. Willms E, Cabanas C, Mager I, Wood MJA, Vader P. Extracellular vesicle heterogeneity: subpopulations, isolation techniques, and diverse functions in Cancer Progression. Front Immunol. 2018;9:738.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Li SR, Man QW, Gao X, Lin H, Wang J, Su FC, Wang HQ, Bu LL, Liu B, Chen G. Tissue-derived extracellular vesicles in cancers and non-cancer diseases: Present and future. J Extracell Vesicles. 2021;10(14):e12175.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Xavier CPR, Belisario DC, Rebelo R, Assaraf YG, Giovannetti E, Kopecka J, Vasconcelos MH. The role of extracellular vesicles in the transfer of drug resistance competences to cancer cells. Drug Resist Updat. 2022;62:100833.

    Article  CAS  PubMed  Google Scholar 

  10. Papp AC, Azad AK, Pietrzak M, Williams A, Handelman SK, Igo RP Jr., Stein CM, Hartmann K, Schlesinger LS, Sadee W. AmpliSeq transcriptome analysis of human alveolar and monocyte-derived macrophages over time in response to Mycobacterium tuberculosis infection. PLoS ONE. 2018;13(5):e0198221.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Li W, Turner A, Aggarwal P, Matter A, Storvick E, Arnett DK, Broeckel U. Comprehensive evaluation of AmpliSeq transcriptome, a novel targeted whole transcriptome RNA sequencing methodology for global gene expression analysis. BMC Genomics. 2015;16:1069.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Wang L, Feng Z, Wang X, Wang X, Zhang X. DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010;26(1):136–8.

    Article  PubMed  Google Scholar 

  13. Kim MW, Park S, Lee H, Gwak H, Hyun KA, Kim JY, Jung HI, Kim SI. Multi-miRNA panel of tumor-derived extracellular vesicles as promising diagnostic biomarkers of early-stage breast cancer. Cancer Sci. 2021;112(12):5078–87.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Zweig MH, Campbell G. Receiver-operating characteristic (ROC) plots: a fundamental evaluation tool in clinical medicine. Clin Chem. 1993;39(4):561–77.

    Article  CAS  PubMed  Google Scholar 

  15. Famta P, Shah S, Chatterjee E, Singh H, Dey B, Guru SK, Singh SB, Srivastava S. Exploring new Horizons in overcoming P-glycoprotein-mediated multidrug-resistant breast cancer via nanoscale drug delivery platforms. Curr Res Pharmacol Drug Discov. 2021;2:100054.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Vasiliou V, Vasiliou K, Nebert DW. Human ATP-binding cassette (ABC) transporter family. Hum Genomics. 2009;3(3):281–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Giddings EL, Champagne DP, Wu MH, Laffin JM, Thornton TM, Valenca-Pereira F, Culp-Hill R, Fortner KA, Romero N, East J, et al. Mitochondrial ATP fuels ABC transporter-mediated drug efflux in cancer chemoresistance. Nat Commun. 2021;12(1):2804.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  18. Thery C, Witwer KW, Aikawa E, Alcaraz MJ, Anderson JD, Andriantsitohaina R, Antoniou A, Arab T, Archer F, Atkin-Smith GK, et al. Minimal information for studies of extracellular vesicles 2018 (MISEV2018): a position statement of the International Society for Extracellular Vesicles and update of the MISEV2014 guidelines. J Extracell Vesicles. 2018;7(1):1535750.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Mazzara S, Rossi RL, Grifantini R, Donizetti S, Abrignani S, Bombaci M. CombiROC: an interactive web tool for selecting accurate marker combinations of omics data. Sci Rep. 2017;7:45477.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  20. Huang M, O’Shaughnessy J, Zhao J, Haiderali A, Cortes J, Ramsey S, Briggs A, Karantza V, Aktan G, Qi CZ, et al. Evaluation of Pathologic Complete Response as a Surrogate for Long-Term Survival outcomes in Triple-negative breast Cancer. J Natl Compr Canc Netw. 2020;18(8):1096–104.

    Article  PubMed  Google Scholar 

  21. Kalluri R, LeBleu VS. The biology, function, and biomedical applications of exosomes. Science 2020;367(6478).

  22. Yu DD, Wu Y, Shen HY, Lv MM, Chen WX, Zhang XH, Zhong SL, Tang JH, Zhao JH. Exosomes in development, metastasis and drug resistance of breast cancer. Cancer Sci. 2015;106(8):959–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Lv MM, Zhu XY, Chen WX, Zhong SL, Hu Q, Ma TF, Zhang J, Chen L, Tang JH, Zhao JH. Exosomes mediate drug resistance transfer in MCF-7 breast cancer cells and a probable mechanism is delivery of P-glycoprotein. Tumour Biol. 2014;35(11):10773–9.

    Article  CAS  PubMed  Google Scholar 

  24. Chen WX, Zhong SL, Ji MH, Pan M, Hu Q, Lv MM, Luo Z, Zhao JH, Tang JH. MicroRNAs delivered by extracellular vesicles: an emerging resistance mechanism for breast cancer. Tumour Biol. 2014;35(4):2883–92.

    Article  CAS  PubMed  Google Scholar 

  25. Kovalev AA, Tsvetaeva DA, Grudinskaja TV. Role of ABC-cassette transporters (MDR1, MRP1, BCRP) in the development of primary and acquired multiple drug resistance in patients with early and metastatic breast cancer. Exp Oncol. 2013;35(4):287–90.

    CAS  PubMed  Google Scholar 

  26. Bellomo C, Caja L, Moustakas A. Transforming growth factor beta as regulator of cancer stemness and metastasis. Brit J Cancer. 2016;115(7):761–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Xu X, Zhang M, Xu F, Jiang S. Wnt signaling in breast cancer: biological mechanisms, challenges and opportunities. Mol Cancer. 2020;19(1):165.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Yao T, Cao R, Xiao W, Pan F, Li X. An optical study of drug resistance detection in endometrial cancer cells by dynamic and quantitative phase imaging. J Biophotonics. 2019;12(7):e201800443.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Kim MH, Kim J, Hong H, Lee SH, Lee JK, Jung E, Kim J. Actin remodeling confers BRAF inhibitor resistance to melanoma cells through YAP/TAZ activation. Embo J. 2016;35(5):462–78.

    Article  CAS  PubMed  Google Scholar 

  30. Begicevic RR, Falasca M. ABC transporters in Cancer Stem cells: beyond Chemoresistance. Int J Mol Sci 2017;18(11).

  31. Kim B, Stephen SL, Hanby AM, Horgan K, Perry SL, Richardson J, Roundhill EA, Valleley EM, Verghese ET, Williams BJ, et al. Chemotherapy induces Notch1-dependent MRP1 up-regulation, inhibition of which sensitizes breast cancer cells to chemotherapy. BMC Cancer. 2015;15:634.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Kim B, Fatayer H, Hanby AM, Horgan K, Perry SL, Valleley EM, Verghese ET, Williams BJ, Thorne JL, Hughes TA. Neoadjuvant chemotherapy induces expression levels of breast cancer resistance protein that predict disease-free survival in breast cancer. PLoS ONE. 2013;8(5):e62766.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  33. Keyvani-Ghamsari S, Khorsandi K, Rasul A, Zaman MK. Current understanding of epigenetics mechanism as a novel target in reducing cancer stem cells resistance. Clin Epigenetics. 2021;13(1):120.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Wainwright EN, Scaffidi P. Epigenetics and Cancer Stem cells: unleashing, hijacking, and restricting Cellular plasticity. Trends Cancer. 2017;3(5):372–86.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Modi A, Roy D, Sharma S, Vishnoi JR, Pareek P, Elhence P, Sharma P, Purohit P. ABC transporters in breast cancer: their roles in multidrug resistance and beyond. J Drug Target 2022:1–21.

  36. Zhang YS, Yang C, Han L, Liu L, Liu YJ. Expression of BCRP/ABCG2 protein in invasive breast Cancer and response to Neoadjuvant Chemotherapy. Oncol Res Treat. 2022;45(3):94–101.

    Article  CAS  PubMed  Google Scholar 

  37. He J, Fortunati E, Liu DX, Li Y. Pleiotropic roles of ABC transporters in breast Cancer. Int J Mol Sci 2021;22(6).

  38. Zhou E, Li Y, Wu F, Guo M, Xu J, Wang S, Tan Q, Ma P, Song S, Jin Y. Circulating extracellular vesicles are effective biomarkers for predicting response to cancer therapy. EBioMedicine. 2021;67:103365.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Konig L, Kasimir-Bauer S, Bittner AK, Hoffmann O, Wagner B, Santos Manvailer LF, Kimmig R, Horn PA, Rebmann V. Elevated levels of extracellular vesicles are associated with therapy failure and disease progression in breast cancer patients undergoing neoadjuvant chemotherapy. Oncoimmunology. 2017;7(1):e1376153.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Tian F, Zhang S, Liu C, Han Z, Liu Y, Deng J, Li Y, Wu X, Cai L, Qin L, et al. Protein analysis of extracellular vesicles to monitor and predict therapeutic response in metastatic breast cancer. Nat Commun. 2021;12(1):2536.

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We would like to thank MID (Medical Illustration & Design), a part of the Medical Research Support Services of Yonsei University College of Medicine, for providing excellent support with medical illustrations and Editage (www.editage.com) for English language editing.

Funding

This study was supported by the Severance Hospital Research Fund for Clinical Excellence (C-2023-0006) and the National Research Foundation of Korea Grants (2021R1I1A1A01051594 and 2022R1F1A1074605).

Author information

Authors and Affiliations

Authors

Contributions

MWK, SIK, and JYK* conceptualized the study. JYK* supervised the study. MWK, HL, SM, SL, YK, and JYK curated the data. MWK and YK analyzed the data. HL, SM, SL, and JYK interpreted the results. MWK, SIK, and JYK* drafted the manuscript, and MWK and HL substantially revised it. All authors have reviewed and approved the submitted version of the manuscript. JYK represents Joon Ye Kim and JYK* represents Jee Ye Kim.

Corresponding authors

Correspondence to Seung Il Kim or Jee Ye Kim.

Ethics declarations

Ethics approval and consent to participate

This study was approved by the Institutional Review Board of Yonsei University Health System (IRB No. 4-2020-1292; approved on January 4, 2021). Written Informed consent for the use of blood samples for research purposes was obtained from all patients, and the study was performed in accordance with the Declaration of Helsinki.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kim, M.W., Lee, H., Lee, S. et al. Drug-resistant profiles of extracellular vesicles predict therapeutic response in TNBC patients receiving neoadjuvant chemotherapy. BMC Cancer 24, 185 (2024). https://doi.org/10.1186/s12885-024-11822-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-024-11822-9

Keywords