Abstract
Background/Aim: THEMIS2, a protein that plays a role in immune regulation, has recently been identified as a key factor in breast cancer progression. This study aimed to identify novel oncogenic regulators in circulating tumor cells (CTCs) that contribute to breast cancer metastasis and chemoresistance, elucidating the regulatory axis involving THEMIS2, VEGFR2 and microRNAs (miRNAs), specifically miR-125b-5p.
Materials and Methods: Using RNA-sequencing of breast cancer CTCs, THEMIS2 was identified as a significantly upregulated gene. miRNA target predictions and expression profiling in primary tumors and normal breast specimens were conducted to evaluate candidate regulators. Functional assays including sphere formation, migration, luciferase reporter activity, VEGFR phosphorylation analysis and in vivo tumor models were employed to investigate the roles of THEMIS2 and miR-125b-5p in metastasis and drug sensitivity.
Results: THEMIS2 was found upregulated in CTCs and breast cancer cell lines with metastatic potential. miR-125b-5p, but not miR-125a-5p, was significantly down-regulated in metastatic models and identified as a direct post-transcriptional suppressor of THEMIS2. THEMIS2 enhanced VEGFR2 phosphorylation and promoted metastatic traits and resistance to docetaxel and bevacizumab. miR-125b-5p suppressed metastasis and restored chemosensitivity by targeting 3′UTR of THEMIS2, an effect reversible with anti-miR-125b or ectopic THEMIS2 expression.
Conclusion: We unraveled a novel miR-125b-5p–THEMIS2–VEGFR2 signaling axis as a key modulator of breast cancer metastasis and chemoresistance. These findings provide mechanistic insight and suggest that miR-125b-5p or THEMIS2 may serve as therapeutic targets or prognostic biomarkers in aggressive breast cancers.
- Circulating tumor cells (CTCs)
- THEMIS2
- miR-125b-5p
- VEGFR2 signaling
- breast cancer metastasis
- chemoresistance
Introduction
Metastasis remains the primary cause of mortality in patients with breast cancer that stresses the critical need for enhanced methods to understand the molecular mechanisms driving cancer spread and progression. Circulating tumor cells (CTCs) are rare malignant cells detached from the primary tumor that circulate within the bloodstream and have recently emerged as valuable liquid biopsy tools for providing insights into the metastatic cascade and tumor heterogeneity (1). Unlike traditional biopsy approaches that are invasive and often limited in capturing metastatic dynamics, analysis of CTCs offers an invasive revenue into the evolving genetic landscape associated with tumor metastasis (2, 3). Recent technological advancements have considerably improved the isolation, identification and characterization of CTCs. Single-cell RNA sequencing technologies have unveiled unique mutational profiles and gene expression patterns within CTCs that frequently diverge from those observed in the primary tumor, identifying novel oncogenic pathways and potential therapeutic targets involved in metastatic progression (2, 4).
Research on CTCs suggested epithelial-to-mesenchymal transition (EMT) and stem cell-like properties as potent underlying mechanisms breast cancer metastasis, heterogeneity and plasticity. For instance, stemness markers such as ALDH1 and EMT-associated transcription factors, including TWIST, AKT2, and PI3Kα, have both been correlated with a highly tumorigenic and chemoresistant subpopulation of CTCs, emphasizing their critical role in metastatic competency (5, 6). Vascular endothelial growth factors (VEGFs) and corresponding receptors (VEGFRs) that are known for their roles in angiogenesis are recently implicated for playing a pivotal role in the metastatic cascade by facilitating tumor growth and the establishment of distant metastatic sites (7, 8). CTCs have been demonstrated to modulate angiogenic milieu through secretion of VEGFs and other proangiogenic factors, enhancing vascular permeability and promoting metastatic colonization (9, 10). Importantly, targeted inhibition of VEGFR signaling has demonstrated potential therapeutic efficacy in preclinical models, underscoring the importance of unraveling the interactions between CTCs and angiogenesis-related pathways (11, 12).
Our study utilized RNA-seq analysis of CTCs to systematically identify whether the novel oncogene THEMIS2 is potentially associated with breast cancer metastatic potency.
Materials and Methods
Cell culture, cell proliferation and sphere formation assay. Human breast cancer cell lines including MDA-MB-231 (ATCC number is HTB-26), MDA-MB-231-IV2, MDA-MB-468 (ATCC number is HTB-132), Hs578T (ATCC number is HTB-126) and Hs578T-FPI were cultured and used for assays as previously described (13, 14). Briefly, the cells were cultured in DMEM containing 10 % FBS for no longer than 3 months and were routinely examined to be mycoplasma-free and authenticated by the assessment of short tandem repeat loci profiling using the Promega StemElite ID System (Madison, WI, USA). Cell proliferation rates were determined by staining with 0.4% Trypan blue solution using the LUNA Automated Cell Counter (Logos BioSystem, Miami, FL, USA) at the indicated time points. For colony formation, cells (1×104 cells/well) were seeded and cultured for 7 days before the colonies were stained with crystal violet. Experiments were performed in triplicates. Cancer stem cells derived from these cell lines were cultured in a stem cell selective condition to allow sphere formation as described previously (14). The experiments were repeated at least three times, and the results are presented as means±SEM. Invasive cell lines such as MDA-MB-231-IV2 and Hs578T-FPI cells were selected by injection of the MDA-MB-231 or Hs578T cells into tail vein of CB17-SCID mice and isolating lung metastasis to derive the highly metastatic subline named IV2 or FPI as previously described (15).
Ethics approval (animals). All animal experiments in this study adhered to the ARRIVE guidelines and conformed to the 3Rs principles (Replacement, Reduction, and Refinement). Our animal protocol number is CMUIACUC-2021-140-2.
The following describes our practices concerning animal care and use: 1) Housing environment: Experimental mice were housed in an Individually Ventilated Cage (IVC) system, ensuring controlled air quality and providing adequate space in accordance with IACUC guidelines. The animal facility maintained a constant temperature of 20-24°C and humidity of 40-70%, following a 12-hour light/12-hour dark cycle. Bedding and nesting material: Cages were provided with dust-free, absorbent bedding (e.g., sterilized wood shavings or paper pulp bedding), which was changed regularly (at least twice per week or as required by IACUC guidelines). For environmental enrichment, clean paper nesting material or species-appropriate shelters (e.g., PVC tubes, small cardboard boxes) were placed in the cages to meet the mice’s natural nesting and hiding behaviors. Diet: Mice had ad libitum access to autoclaved commercial rodent chow, ensuring a balanced and consistent nutrient source. Water: Sterilized drinking water was provided via water bottles or an automated watering system, with daily checks to ensure an uninterrupted supply. 2) During the experiment, the physical condition of all mice was monitored at least once daily (including weekends and public holidays) by trained personnel. Monitoring frequency was increased when necessary (e.g., in cases of abnormalities or after experimental procedures). 3) No animals became severely ill or died at any time prior to the experimental endpoint. As no animals died prior to the experimental endpoint and no animals became ill. 4) We strictly adhered to the IACUC-approved protocol for humane endpoints for experimental animals. This protocol ensures that when an animal exhibits irreversible pain or illness, euthanasia is performed promptly and decisively to minimize its suffering. Clinical signs used to determine when to euthanize animals (humane endpoints) included, but were not limited to, any one or more of the following: Significant Weight Loss: Body weight decreased by more than 15-20% of the starting research weight (or continuous decline for three consecutive days). Significant Deterioration of Body Condition Score (BCS): For example, a BCS score of 1 (emaciated, bony prominences visible). Persistent Signs of Pain or Distress: Manifested as hunching, piloerection (ruffled fur), semi-closed eyes, hunched posture, withdrawal, or slow/no response to stimuli. Severe Impairment of Mobility: Inability to walk, stand, or access food and water normally, or exhibiting lethargy and reluctance to move. Tumor-Related Symptoms: Extensive ulceration, necrosis, or severe infection of the tumor that cannot be effectively treated. Tumor size reaching the IACUC-stipulated limit (1.5 cm, 15 mm, in this study) or significantly affecting animal function even if below the limit. Tumor causing obvious compressive symptoms (e.g., dyspnea, paralysis). Inability to Eat or Drink Normally: Continuous lack of food or water intake for more than 24 hours. Respiratory Distress: Rapid, labored breathing, open-mouth breathing, or cyanosis. Other severe and irreversible disease states as assessed by a veterinarian. Once any of the above humane endpoint criteria were observed, immediate consultation with a veterinarian and euthanasia were performed. 5) Euthanasia for all animals in this study was performed in accordance with IACUC-approved standard operating procedures (SOPs). The primary method of euthanasia was inhalation of high-concentration carbon dioxide (CO2) gas, following these steps: Mice were placed individually or in small groups (to avoid overcrowding) into a CO2-pre-filled euthanasia chamber. CO2 gas was introduced gradually, or directly into a pre-filled chamber, to ensure the animals did not experience distress from sudden gas exposure. After confirming complete unconsciousness, exposure to CO2 was continued for at least 5 minutes to ensure complete death (cessation of heartbeat and respiration). To ensure complete death, a secondary method of confirmation was used after CO2 exposure, such as cervical dislocation, pneumothorax/thoracotomy and organ exsanguination. All euthanasia procedures were performed by personnel trained in the method.
Cell chemotactic migration and invasion assay. Control, specific plasmid, siRNA or miRNA transfected cells (1×105) were harvested from monolayer culture, resuspended in serum-free medium with 0.1% bovine serum albumin (Sigma-Aldrich, St. Louis, MO, USA), and then plated onto a Transwell Boyden chamber, with or without matrigel coating on the upper chamber, for invasion or migration assays, respectively. The chambers were incubated in DMEM for 24 h with the addition of 10% FBS to the lower chamber. Cells that did not move were removed using cotton swabs, and the chambers were stained with crystal violet. Photomicrographs of three random regions were captured from duplicated assay chambers. The cell numbers were counted and normalized to the control. All experiments were repeated independently at least three times.
Immunoblot analysis. Cell lysates (30 μg per lane) were separated using 10% SDS-PAGE and transblotted onto a polyvinylidene difluoride (PVDF) membrane (Millipore, Bedford, MA, USA), which was incubated at room temperature in a phosphate-buffered saline (PBS) blocking solution 5 % powdered skim milk and 0.02% sodium azide for 1 h. This was followed by further 2-h incubation with primary antibody including anti-THEMIS2 (A11109, Boster Biological Technology, Pleasanton, CA, USA), anti-p-VEGFR1 (AF4170, R&D Systems, Minneapolis, MN, USA), anti-VEGFR2 (55B11, Rabbit mAb#2479, Cell Signaling, Boston, MA, USA) anti-p-VEGFR2 (Tyr1175, 19A10, Rabbit mAb#2478, Cell Signaling), anti-p-VEGFR3 (CY1115, Cell Applications Inc., San Diego, CA, USA) or anti-β-actin (MAB1501, Chemicon, Temecula, CA, USA) at appropriate dilution (1:500 to 1:2,000). The membrane was washed three times and incubated with horseradish peroxidase conjugated secondary anti-mouse or anti-rabbit immunoglobulin (IgG) (1:5,000; Sigma-Aldrich) for 1 hour prior to signal development by chemiluminescence reagent (Amersham Pharmacia Biotech, Little Chalfont, UK).
Transfection and RNA inference. Non-specific and human THEMIS2-specific SMART-pool-targeting small interfering RNA molecules were obtained from Dharmacon (Lafayette, CO, USA), and the sequences were identical to previous study (16). Briefly, siRNAs were transfected into selected cells for 24 hours using a specialized transfection reagent (Catalog number sc-29528, Santa Cruz Biotechnology, Dallas, TX, USA) according to the manufacturer’s protocol. The pCMV6-puro plasmid containing full-length THEMIS2 cDNA was transfected into MDA-MB-231, MDA-MB-231-IV2, Hs578T or Hs578T-FPI cell lines, and stable THEMIS2-expressing cancer cell lines were selected by using puromycin (Sigma-Aldrich). miRNAs including miR-125a-5p, miR-125b-5p, anti-mi125a-5p and anti-miR125b-5p were transfected into cells using Lipofectamine 2000 (Invitrogen, Carlsbad, CA, USA). Western blot was used to determine the effects of gene silencing or overexpression.
In vivo tumor growth and metastasis assay. For the tumor growth assay, 1×106 cells of control, miR-125b-5p-transfected or miR-125b-5p-transfected/THEMIS2-overexpressing Hs578T cells were orthotopically implanted into the mammary fat pads of CB17-SCID mice. After tumors reached 0.8 to 1 cm3, tumor pieces were surgically implanted under the ovarian bursa of additional mice (n=5 for each group). In accordance with the Institutional Animal Care and Use Committee (IACUC) guidelines, the orthotopic tumors in this study’s mice are generally not to exceed 1.5 cm (15 mm). The tumor volumes were monitored twice per week using the IVIS imaging system for bioluminescence detection over a 3-week period. In vivo metastasis assay was performed, in which Hs578T-FPI cells (5×105) transfected with both of miR-125b-5p and Anti-125b-5p, miR-125b-5p alone, or control miRNA were suspended in PBS and injected individually into the tail vein of 6 to 8-week-old female CB17-SCID mice, which were then monitored for 30 to 60 days before sacrifice. Lung tissues were removed, fixed, paraffin-embedded, serially sectioned, and subjected to hematoxylin and eosin (H/E) staining. Tumor volumes were measured once a week, and tumor growth was evaluated until 4 months after implantation. All procedures were performed in accordance with the Guide for Care and Use of Laboratory Animals issued by the Institutional Animal Care and Use Committee of China Medical University, Taiwan (CMUIACUC-2021-140-2).
Kaplan-Meier plotter analysis. Kaplan-Meier plots were generated through the automatic selection of the optimal cutoff values between lower and upper quartiles, which formed the high and low expression groups. This analysis was performed on patients with breast cancer who had either received effective chemotherapy treatment or no treatment (https://kmplot.com/analysis/). Clustering of the transcriptomic profiles of 78 triple-negative breast cancer (TNBC) cases was obtained from the ONCOMINE database.
Transcriptomic sequencing data analysis. Two transcriptomic RNA-seq datasets were obtained from Gene Expression Omnibus (GEO) hosted by National Center of Biotechnology Information (NCBI, https://www.ncbi.nlm.nih.gov/geo/). The first dataset (GSE268201) was generated using single circulating tumor and normal cells that were isolated using density gradient centrifugation and CD45-based immunomagnetic depletion from liquid biopsies of breast cancer patients. The second dataset (GSE7842) that contains miRNA expression data from 93 primary human breast tumors and five normal breast samples was utilized to assess expressions of miRNA-125a-5p and miRNA-125b-5p within normal breast samples versus primary breast tumors. Briefly, the read alignment was conducted using aligners HISAT2 (16), transcript abundance was quantified by counting mapped reads using Kalisto (17), followed by differential gene expression analysis and normalization on edgeR as previously described (15).
Statistical analysis. All of the in vitro experiments were reproducible and repeated independently at least three times, and the results are presented as means±SEM. Except where otherwise noted, data presented in figures are shown as mean±standard deviation (SD). Independent samples t-tests were conducted using SPSS Statistics for Windows, version 18.0 (SPSS Inc., Chicago, IL, USA) to determine differences in the frequencies and means of the categorical data, p-values less than 0.05 were considered statistically significant.
Results
Differential THEMIS2-targeting miRNA expression unveils prognostic relevance in breast cancer metastatic potency. Detectable levels of EMT markers such as TWIST, AKT2 and PI3Kalpha in addition to ALDH1 have been shown to define a highly tumorigenic subset of CTCs and advocate the clinical value of CTC in predicting drug resistance and patient prognosis (18). Nonetheless, the molecular mechanisms underlying the role of CTC in drug resistance and prognosis remains elusive. In this study, RNA-seq data analysis using liquid biopsies from breast cancer and normal patient panels demonstrated a cluster of genes including THEMIS2 that was significantly up-regulated in the CTCs (Figure 1a). Subsequently, three potential miRNAs, including hsa-miR-125a-5p, hsa-miR-125b-5p, and hsa-miR-4319, were predicted to interact with the predicted 8-mer target site within the 3′ untranslated region (3′UTR) of the THEMIS2 gene (Figure 1b). These three potent miRNAs that specifically target THEMIS2 were interrogated by comparative analysis of miRNA expression between normal breast specimens and primary breast tumors. The heatmap analysis elicited significant differential expression of miR-125b between the normal breast samples (3.6-fold change) and breast tumors (2.8-fold change), while only slight variation for miR-125a between normal (2.9-fold change) and tumor samples (2.87-fold change) was observed (Figure 1c). Of note, the role of miR-125a-5p and miR-125b in breast cancer prognosis was extrapolated by Kaplan-Meier plotter analysis. The overall survival (OS) results showed that miR-125a-5p expression did not significantly correlate with survival (log-rank p=0.072, HR=1.23), whereas higher miR-125b-5p expression was significantly associated with improved OS (log-rank p=0.0012, HR=0.71) (Figure 1d). The potency of miR-125a and miR125b in driving breast cancer metastasis was next investigated using two pairs of metastatic breast cancer cell lines generated from xenograft models in vivo. In line with OS analysis, the quantitative PCR (qPCR) data demonstrated that mildly reduced the expression of miR-125a-5p was insignificant between the metastatic (MDA-MB-231-IV2 and Hs578T-FPI) and parental (MDA-MB-231 and Hs578T) cell lines, while miR-125b-5p expression was significantly mitigated in the metastatic cells when compared to the controls (Figure 1e).
Differential expression and prognostic relevance of miRNAs in breast cancer. (a) Heatmap illustrating DEGs between single circulating tumor cells and normal cells in liquid biopsies from breast cancer and normal patients, showing THEMIS2 as one of the top differentially up-regulated genes in CTCs. The RNA-seq data from GEO268201 were analyzed as described in Materials and Methods. (b) Potential miRNA pairing sites in the 3′UTR of THEMIS2 were predicted using TargetScan, which uncovered consequential pairing target regions for THEMIS2, identifying miRNAs including miR-125a-5p, miR-125b-5p and miR-4319. (c) Comparative heatmap analysis of miRNA expression between normal breast tissue and primary breast tumor samples providing further insights into miRNA expression for miR-125a and miR-125b. (d) Kaplan-Meier survival curves correlating high and low expression levels of miR-125a-5p and miR-125b with patient OS, indicating prognostic significance specifically for miR-125b. (e) qPCR analysis was conducted using two metastasis pairs of breast cancer cell lines (MDA-MB-231-IV and Hs578T-FIPI) by normalizing against RNU6B. By comparison to parental cell lines (MDA-MB-231 and Hs578T), significant reduction of miR-125b in metastatic breast cancer cells was observed. Histograms in (e) represent means±SD from 3 independent experiments (**p<0.01 compared to parental control cells.) DEGs: Differentially expressed genes; CTCs: circulating tumor cells; OS: overall survival.
miR-125b-5p but not miR-125a-5p impacts on metastatic potential and chemosensitivity. Given miR-125a is identified as a tumor suppressor miRNA that plays a crucial role in regulating cell proliferation, apoptosis and migration in various cancers (19, 20), the insignificant changes in miR-125a expression in breast tumors and metastatic breast cancer cell lines (Figure 1c, e) prompted us to investigate whether and how miR-125a could impact on oncogenic properties of breast cancer, by comparison to miR-125b. Colony formation assays performed on metastatic breast cancer cell lines (MDA-MB-231-IV2 and Hs578T-FPI), derived from xenograft models through selection of lung metastatic populations, revealed that only transfection with miR-125b-5p significantly (p<0.01) suppressed colony formation compared to control groups, whereas miR-125a-5p overexpression appeared to mildly enhance colony formation in both metastatic cell lines (p>0.05) (Figure 2a). Consistently, sphere formation assays utilizing Hs578T-FPI cells under stem cell-selective culture conditions showed that miR-125b-5p but not miR-125a-5p significantly impaired the cellular ability to form spheres, indicative of reduced stemness and metastatic capacity when miR-125b-5p was overexpressed (p<0.01) (Figure 2b). Further analysis through angiogenic tube formation assays using human umbilical vein endothelial cells (HUVECs) cultured in conditioned media derived from Hs578T-FPI cells demonstrated significantly impaired angiogenic potential by miR-125b-5p (p<0.01) and insignificant changes in the formation of vascular-like structures by miR-125a-5p (Figure 2c). Moreover, to evaluate the influence of these miRNAs on chemotherapy sensitivity, cell viability assays were carried out on Hs578T and MDA-MB-231-IV2 cells treated with escalating doses of docetaxel. Our data showed that miR-125b-5p overexpression significantly enhanced docetaxel sensitivity in both cell lines in a dose-dependent manner (p<0.01), whereas miR-125a-5p transfection showed no statistically significant effect on chemosensitivity (p=0.817 or 0.501) (Figure 2d). Further, the role of miR-125b-5p in mediating chemoresistance was confirmed by anti-miRNA oligonucleotides (AMOs) anti-miR-125a-5p specifically designed to target miR-125b-5p, which significantly inhibited the chemoresistance-promoting activity of miR-125b-5p (p<0.05) (Figure 2e).
Functional characterization of miR-125a-5p and miR-125b-5p on breast cancer cell metastatic potential and chemosensitivity. (a) Colony formation assays were performed by seeding metastatic breast cancer cell lines (MDA-MB-231-IV2 and Hs578T-FPI) at 1×104 cells/well and cultured for 7 days before staining with crystal violet. Statistically significant reductions in colony numbers following miR-125a-5p and miR-125b-5p transfections compared to control were observed (p<0.01). (b) Sphere formation assays involved culturing Hs578T-FPI cells and were carried out under stem cell-selective conditions. Representative microscopy images illustrating inhibitory effects of miR-125b-5p (p<0.01), but not miR-125a-5p are shown. (c) Angiogenic potential was evaluated by culturing HUVECs in conditioned media from Hs578T-FPI cells transfected with miRNAs; fluorescence microscopy depicted significantly impaired tube formation in miR-125b-5p conditioned medium (p<0.01). (d) Chemosensitivity was assessed using cell viability assays, wherein breast cancer cells (Hs578T and MDA-MB-231-IV2) transfected with miRNAs were exposed to varying concentrations of docetaxel for 24 h. (e) AMOs were specifically designed to evaluate whether inhibitory effect of miR-125b-5p on cellular chemosensitivity could be observed in the same cell lines used in (d). All the in vitro experiments were performed in triplicates and repeated three times (*p<0.05; **p<0.01.). Hs578T-FPI: Hs578T-fat pad injection; HUVECs: human umbilical vein endothelial cell; AMOs: anti-miRNA oligonucleotides.
miR-125b-5p directly targets THEMIS2 to suppress breast cancer metastasis. The molecular mechanisms by which miR-125b-5p regulates THEMIS2 expression and influences breast cancer metastatic behaviors were further investigated. Comparative sequence analysis across species revealed a conserved miR-125b-5p binding site within the 3′UTR of THEMIS2 (Figure 3a). This conserved miR-125b-5p element alongside its mutant were subcloned into a luciferase reporter system. Results from luciferase reporter assays demonstrated that miR-125b-5p significantly impaired luciferase activity in cells transfected with wild-type THEMIS2 3′UTR in both Hs578T and IV2 cells (both p<0.01), whereas mutant THEMIS2 3′UTR abolished this inhibitory effect, confirming the specificity of interaction within the evolutionarily conserved 3′UTR (Figure 3b). Consistently, qPCR and immunoblot analyses showed pronounced down-regulation of THEMIS2 expression upon miR-125b-5p transfection, which was specifically reversed when the cells were transfected with AMOs anti-miR-125b-5p (Figure 3c, d). Importantly, functional assays that depict the oncogenic consequences of THEMIS2 modulation by miR-125b-5p demonstrated significant suppression of invasive capabilities by miR-125b-5p was reversible by anti-miR-125b-5p in transwell assays (both p<0.01), whereas significantly mitigated sphere formation in stem cell-selective conditions by miR-125b-5p was similarly restored upon transfection with the AMOs (both p<0.01) (Figure 3e, f). Notably, significantly reduced metastatic potential in vivo was consistently reversible by anti-miR-125b-5p, as shown by hematoxylin and eosin staining of lung metastasis sections from xenografted mice (both p<0.01) (Figure 3g).
miR-125b-5p specifically targets THEMIS2 to inhibit metastatic potential in breast cancer cells. (a) Comparative sequence analysis identified conserved miR-125b-5p binding sites across species in THEMIS2 3′UTR at 111 to 118 nt, from which a 3′UTR mutant was designed and subcloned into a luciferase reporter vector (THEMIS2 3′UTR-mt-luc). (b) Luciferase reporter assays using wild-type (THEMIS2 3′UTR-wt-luc) and the mutant vectors in Hs578T and IV2 cells were conducted to determine the bona fide 3′UTR region of THEMIS2. (c-d) qPCR and immunoblotting were carried out to confirm THEMIS2 down-regulation by miR-125b-5p transfection and whether it was reversible by AMOs anti-miR-125b-5p. (e-f) Transwell invasion assays, sphere formation assays in vitro and (g) lung metastasis (×1.6) evaluations in xenografted mice were conducted to further investigate the functional impacts of miR-125b-5p and anti-miR-125b-5p on cellular invasiveness, sphere formation and metastatic capacities. H/E staining was performed as described in Materials and Methods and representative images were shown for each group. All the in vitro experiments were performed in triplicates and repeated three times (**p<0.01.). AMOs: Anti-miRNA oligonucleotides.
miR-125b-5p mediates suppression of breast cancer metastatic potential via specific targeting of THEMIS2-VEGFR signaling. Since angiogenesis is known to take part in cancer metastasis, and vascular endothelial growth factor (VEGF) is an important protein involved in this process. The ability of miR-125b-5p in ablating angiogenic properties from Hs578T-FPI condition medium (Figure 2c) led us to elucidate the underlying signaling mechanisms by assessing expression levels of phosphorylated VEGFR family proteins. Firstly, three VEGFRs (VEGFR1, VEGFR2 and VEGFR3) were subjected to immunoblot analyses, which indicated that miR-125b-5p markedly reduced phosphorylation of VEGFR1, VEGFR2 as well as VEGFR3 as compared to miRNA control in Hs578T cells (Figure 4a). Knockdown of THEMIS2 was subsequently employed via THEMIS2-specific siRNA that similarly resulted in impaired VEGFR2 phosphorylation in both Hs578T and IV2 cells, supporting the involvement of THEMIS2 in VEGFR-mediated signaling pathways. Conversely, ectopic expression of THEMIS2 not only led to restored VEGFR2 phosphorylation (Figure 4b and 4c) but also counteracted the inhibitory effects of miR-125b-5p on VEGFR2 phosphorylation (Figure 4d). Furthermore, the functional significance of VEGFR2 signaling modulation was confirmed by assessing cellular chemosensitivity to VEGF inhibitor bevacizumab using cell viability assays. The results demonstrated significantly increased sensitivity of miR-125b-5p-transfected cells to the bevacizumab that was considerably reversed when THEMIS2 was overexpressed (Figure 4e). Consistently, xenograft models further validated these findings in vivo, with significantly mitigated tumor weights observed following miR-125b-5p transfection (p <0.05), an effect that was significantly reversible by ectopic THEMIS2 expression (p<0.01) (Figure 4f).
miR-125b-5p suppresses breast cancer metastasis through modulation of THEMIS2-mediated VEGFR signaling pathways. (a) Immunoblotting analyses were first employed to assess phosphorylation levels of VEGFR1, VEGFR2 and VEGFR3 upon miR-125b-5p transfection; (b) THEMIS2-silencing by THEMIS2-specific siRNA or (c) THEMIS2 overexpression in Hs578T-FPI and IV2 cells was performed to determine the direct impact of THEMIS2 on VEGFR2 phosphorylation. (d) miR-125b-5p-transfected Hs578T-FPI cells that were overexpressed with THEMIS2 were examined for phosphorylation level of VEGFR2. (e) Another chemosensitivity assay was conducted by measuring cell viability of Hs578T-FPI cells that were treated with escalating doses of bevacizumab (0 to 1 mM) in miRNA control, miR-125b-5p- or miR-125b-5p/THEMIS2-transfected groups. (f) The same experimental designs were employed to in vivo tumor growth assays, wherein xenografted CB17-SCID mice with miR-125b-5p or miR-125b-5p/THEMIS2 overexpression were monitored for tumor growth as described in Materials and Methods. All the in vitro experiments were performed in triplicates and repeated three times (*p<0.05; **p<0.01.).
Discussion
The identification of novel oncogenic pathways and therapeutic targets is essential for improving the prognosis and treatment strategies in metastatic breast cancer. The present study unraveled a previously uncharacterized signaling axis involving THEMIS2, miR-125b-5p and VEGFR2, identified through RNA-seq profiling and molecular analysis of CTCs. Our findings provided first evidence that THEMIS2 directly influenced the metastatic and chemoresistant characteristics of breast cancer cells via the VEGFR2 signaling pathway. THEMIS2 knockdown experiments notably decreased phosphorylation levels of VEGFR2, confirming the interaction between THEMIS2 expression and angiogenic signaling. This relationship was corroborated by ectopic THEMIS2 expression, which not only restored VEGFR2 phosphorylation but also counteracted the suppressive effects of miR-125b-5p. In addition, in vitro assays consistently demonstrated enhanced chemosensitivity to VEGF inhibitors such as bevacizumab in miR-125b-5p-transfected cells that could be reversible upon THEMIS2 overexpression. These in vitro findings were further validated by in vivo experiments, which reiterated THEMIS2-mediated VEGFR2 signaling activation, tumor growth and metastatic potential. miRNAs are well-established as critical regulators of cancer progression and therapeutic response. Of the miRNAs examined, miR-125b-5p displayed pronounced anti-metastatic and chemosensitizing effects by directly targeting 3′UTR region of THEMIS2. In contrast to this observation, miR-125a-5p, despite its sequence similarity, elicited negligible impact on THEMIS2 modulation, metastatic capability, or chemoresistance. This differential functional outcome between miR-125b-5p and miR-125a-5p is critical, given the complexity of miRNA regulatory networks in cancer. To date, serum amyloid A-activating factor 1 (SAF-1) is the only molecule identified as a downstream target of miR-125b and elevates the activity of VEGF (21). The ectopic expression of miR-125b is shown to inhibit SAF-1 expression and down-regulates the proliferative and invasive properties of breast cancer cells.
Conversely, miR-125a-5p, which has been proposed as a tumor suppressor miRNA regulating proliferation and apoptosis in multiple cancers (22, 23), displayed only minor expression changes and no significant biological impact in metastatic breast cancer contexts. For instance, miR-125a-5p is reported to target VEGFA and mediate tumor suppression in colorectal and gastric cancers (24, 25), our findings suggest only miR-125b-5p, but not miR-125a-5p, was capable of targeting THEMIS2 to modulate phosphorylation of VEGFRs and impact on metastasis and chemosensitivity of breast cancer. Interestingly, miR-125a-5p is also demonstrated to modulate tumor suppressive properties in breast cancer albeit via targeting of other molecules such as BAP1, HDAC4 and HDAC5 (20, 26, 27). While these studies show the ability of miR-125a-5p to induce apoptosis in breast cancer, involvement of miR-125a-5p in metastatic pathways remained ambiguous. emphasizing the specificity and selective regulatory roles that miRNAs can exhibit in cancer biology. Our comprehensive assays indicated miR-125a-5p might mildly support colony formation without significantly altering metastatic or chemoresistant phenotypes, sharply contrasting with the pronounced inhibitory effects of miR-125b-5p. The distinct differences observed between miR-125b-5p and miR-125a-5p in our study are consistent with emerging evidence that highlight differential miRNA specificity even within closely related families (28, 29).
Furthermore, previous research investigating THEMIS2 in breast cancer did not implicate the VEGFR signaling pathway or identify any regulatory miRNAs. Instead, THEMIS2 was found to enhance cancer stemness and chemoresistance through MET signaling by suppressing PTP1B-mediated dephosphorylation, thereby promoting MET activation (15). Although THEMIS2 was proposed as a scaffold protein regulating EMT-driven cancer stemness in ovarian cancer and a prognostic biomarker in glioblastoma (30, 31), no studies have yet associated it with VEGFR2 activation or angiogenesis, nor identified any miRNA including miR-125b-5p as a post-transcriptional regulator of THEMIS2. The only non-coding RNA relevant to THEMIS2 in the literature is a long non-coding RNA THEMIS2-211 that elicits oncogenic properties by promoting proliferation, migration, invasion and EMT of HCC (32).
Conclusion
Collectively, this study unveils the sophisticated regulatory dynamics involving THEMIS2, miR-125b-5p and VEGFR2 as a central oncogenic pathway in breast cancer metastasis and resistance to chemotherapy. The functional disparity observed between miR-125b-5p and miR-125a-5p underlines the complexity and specificity inherent in miRNA regulation, necessitating precise therapeutic targeting strategies. The elucidation of this novel pathway provides a promising framework for developing targeted interventions aimed at curtailing metastatic spread and overcoming drug resistance, ultimately improving outcomes for patients suffering from metastatic breast cancer.
Acknowledgements
We thank the staff of the Kaplan–Meier plotter (https://kmplot.com/analysis/) for valuable contributions. This work were supported by the Ministry of Science and Technology, Taiwan (MOST 109-2320-B-039–067, MOST 110-2320-B-039–061), and China Medical University, Taiwan (CMU111-N-02, CMU112-S-48 and CMU113-S-53), Ton-Yen General Hospital, Hsinchu County, Taiwan (TYH-RP-112002) and Show Chwan Memorial Hospital (SRD-110051 to WTL).
Footnotes
Authors’ Contributions
W-C H, H-H C, W-T L and S-LT contributed to the experimental design and prepared the manuscript. W-CH, C-J C, I-S L, C-C C and Y-W T carried out the works of experiment. C-J C and H-H C provided and assisted in the collection and analysis of clinical samples. W-C H directed and supervised the research works, and revised the manuscript. All Authors reviewed and edited the manuscript. All Authors read and approved the final manuscript.
Conflicts of Interest
The Authors declare no competing financial or personal interests that could have influenced the work presented in this paper.
Artificial Intelligence (AI) Disclosure
No artificial intelligence (AI) tools, including large language models or machine learning software, were used in the preparation, analysis, or presentation of this manuscript.
- Received September 10, 2025.
- Revision received October 3, 2025.
- Accepted October 21, 2025.
- Copyright © 2026 The Author(s). Published by the International Institute of Anticancer Research.
This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.














