QUE inhibited the viability of breast cancer cell lines
In this study, we used the MTT assay, a commonly used cell viability assay[30] to measure the IC50 of QUE in TNBC cell line viability. Cells were treated with various concentrations of QUE for 12 hours prior to the assay. The study included five triple-negative breast cancer (TNBC) cell lines: MDA-MB-231, MDA-MB-468, MDA-MB-436, BT-20, and BT-549. Additionally, two non-TNBC cell lines, MDA-MB-435 (HER-enriched breast cancer) and MCF-7 (ER-positive breast cancer), as well as two non-tumorigenic epithelial breast cell lines, MCF-10A and MCF-12A, were included. Results demonstrated that QUE significantly inhibited the viability of all TNBC cell lines, as well as the HER-enriched breast cancer cell line MDA-MB-435 and the ER-positive breast cancer cell line MCF-7. The IC50 values for TNBC cell lines ranged from 15.3 µM to 55.2 µM, whereas the IC50 values for the two normal breast cell lines exceeded 100 µM. Among the TNBC cell lines, MDA-MB-231 and BT-20 exhibited the highest sensitivity to QUE, with IC50 values of 15.3 µM and 20.1 µM, respectively, while BT-549 showed the lowest sensitivity, with an IC50 of 55.2 µM. To further confirm the insensitivity of normal breast cells at low concentrations, cells were exposed to QUE for 36 hours and 72 hours, revealing IC50 values still exceeding 100 µM for the normal breast cells. Thus, it was concluded that normal breast cells were not significantly affected by QUE at concentrations below 64 µM. These findings suggested that QUE exhibited relative specificity in inhibiting breast cancer cells, although the differential effect between TNBC and non-TNBC cells may not be substantial. Subsequent experiments utilized MDA-MB-231 and BT-20 (the TNBC cell lines with the highest IC50 values) as well as MCF-10A and MCF-12A (the normal control cell lines), with QUE administered at a concentration of 64 µM.
ORM2 was associated with TNBC and worse survival in TNBC patients
ORM2 is one of the potential targets of QUE provided by the HERB database[31], a high-throughput experiment- and reference-guided database of traditional Chinese medicine. We queried “ORM2” for the drug target and the database search results came out with three drugs, including meletin, plumb, and quercitin. We then search for “quercitin” in the PubChem database and found that it is the same as “quercetin” (PubChem id: 5280343).
According to TCGA data, ORM2 was found to be significantly upregulated in breast cancer compared to normal breast tissues. Furthermore, within breast cancer subtypes, TNBC exhibited higher ORM2 expression compared to non-TNBC (Fig. 2A). We further investigated the association between ORM2 levels and different pN stages in TNBC and non-TNBC. The results revealed that N2 TNBC had higher ORM2 expression compared to N1, and N1 TNBC showed higher ORM2 expression compared to N0. However, in non-TNBC, there were no significant differences observed among different pN stages (Fig. 2B). These findings suggest that ORM2 may be associated with TNBC metastasis or cell survival during migration, but it may not impact the migration of non-TNBC cells.
Moreover, we analyzed the overall survival of TNBC and non-TNBC patients using the TCGA cohort. The analysis demonstrated that high ORM2 expression was significantly associated with worse overall survival in TNBC but not in non-TNBC. The hazard ratio for the high ORM2 group was 5.26 in TNBC (Fig. 2C). These data further support the notion that ORM2 affects TNBC but does not exert the same influence on non-TNBC. Consequently, this study focused exclusively on investigating TNBC.
To validate the overexpression of ORM2 in TNBC, we collected tumor samples and paired normal breast adjacent tissues from TNBC patients. Protein expression analysis was performed using western blotting, and protein levels were observed through immunohistochemistry staining. The results demonstrated that, overall, ORM2 protein was significantly upregulated in tumor samples compared to normal breast tissue samples. Among the 30 patients, 4 patients (13.3%) exhibited lower ORM2 protein levels in cancer tissue compared to normal tissue, while 26 patients (86.7%) showed higher ORM2 protein levels in cancer tissue (Fig. 2D-E). Immunohistochemistry staining of the slides also confirmed stronger ORM2 staining in tumor samples compared to normal breast tissue samples (Fig. 2F). The observed ORM2 protein overexpression was consistent with the ORM2 mRNA expression levels obtained from data mining analysis, suggesting that ORM2 protein levels are dependent on ORM2 mRNA levels. It is worth noting that the mRNA levels of ORM2 may be epigenetically regulated by transcriptional factors, requiring further exploration in future studies. These findings also indicate that TNBC may exhibit distinct ORM2 regulation compared to other breast cancer subtypes.
ORM2 positively regulated the viability of TNBC cells but not normal breast cells
To further investigate the role of ORM2 in breast cancer, we examined ORM2 expression in seven breast cancer cell lines and two normal breast cell lines. Western blotting analysis revealed significantly higher ORM2 protein levels in the breast cancer cell lines compared to the normal breast cell lines (Fig. 3A-B). Based on their high ORM2 expression, MDA-MB-231, BT-20, MCF-10A, and MCF-12A were selected for further experiments on ORM2 knockdown and overexpression. It is worth noting that both MDA-MB-231 and BT-20 are TNBC cell lines.
In this study, we successfully knocked down ORM2 expression in MDA-MB-231 and BT-20 cells, as confirmed by western blotting (Fig. 3C1-2 and D1-2). MTT assay results demonstrated that ORM2 knockdown led to reduced cell viability in both MDA-MB-231 and BT-20 cells (Fig. 3C3 and D3). Cell counting assays further confirmed a decrease in the cell number of MDA-MB-231 and BT-20 cells following ORM2 knockdown (Fig. 3C4 and D4). However, the apoptosis levels in MDA-MB-231 and BT-20 cells showed inconsistent results. ORM2 knockdown increased apoptosis in MDA-MB-231 cells but decreased apoptosis in BT-20 cells (Fig. 3C5 and D5). Considering the overall increase in cell number over time in both MDA-MB-231 and BT-20 cells, we hypothesized that the observed changes in apoptosis were not the primary factor driving the alterations in cell viability upon ORM2 knockdown. Thus, apoptosis may not play a critical role in the effect of ORM2.
Additionally, we conducted similar experiments in the two non-tumorigenic epithelial breast cell lines, MCF-10A and MCF-12A. The overexpression of ORM2 in MCF-10A and MCF-12A cells was confirmed by western blotting following vector transfection (Fig. 3E1-2 and F1-2). MTT assay results indicated that ORM2 overexpression did not affect the viability of MCF-10A and MCF-12A cells (Fig. 3E3 and F3). Cell counting assays further confirmed that ORM2 overexpression did not alter the cell number of MCF-10A and MCF-12A cells (Fig. 3E4 and F4). Furthermore, ORM2 overexpression did not significantly affect apoptosis in MCF-10A and MCF-12A cells (Fig. 3E5 and E5). These findings suggested that ORM2 positively regulated the viability of TNBC cells but did not impact normal breast cells.
QUE inhibited the expression of ORM2 and viability in TNBC cells but not in normal breast cells
We observed the down-regulation of ORM2 protein by QUE treatment in MDA-MB-231 and BT-20 cells, as evidenced by western blotting analysis (Fig. 4A1-2 and B1-2). Cell counting assays revealed a decrease in the cell number of MDA-MB-231 and BT-20 cells following QUE treatment (Fig. 4A3 and B3), while apoptosis remained unaffected (Fig. 4A4 and B4). Conversely, QUE treatment did not affect ORM2 protein levels in MCF-10A and MCF-12A cells (Fig. 4C1-2 and D1-2). Cell counting assays demonstrated that QUE exposure did not impact the cell number of MCF-10A and MCF-12A cells (Fig. 4C4 and D4). Furthermore, QUE treatment did not significantly affect apoptosis in MCF-10A and MCF-12A cells (Fig. 4C4 and D4). These findings suggest that ORM2 positively regulates the viability of TNBC cells, while it does not exert the same effect on normal breast cells.
ORM2 mediated the inhibition of QUE toward the viability of TNBC cells but not normal breast cells.
Correlation analysis revealed a negative correlation between ORM2 protein expression and the IC50 of QUE in breast cancer cell lines (Fig. 5A). This correlation suggests that ORM2 may mediate the effect of QUE on cell viability. To validate this hypothesis in TNBC, we performed ORM2 knockdown experiments in MDA-MB-231 and BT-20 cells using two distinct siRNA constructs (KD1 and KD2). We then measured the IC50 of QUE inhibition on the viability of these cell lines. The vehicle was used as a control to account for any potential impacts of plasmid transfection on cell behavior. As depicted in Fig. 5B, the IC50 of the knockdown control in MDA-MB-231 cells was 13.2 µM, whereas the IC50 of both ORM2-knockdown MDA-MB-231 cells exceeded 100 µM. Similarly, in the case of BT-20 cells, the IC50 of the knockdown control was 23.1 µM, while the IC50 of both ORM2-knockdown BT-20 cells also surpassed 100 µM. These results support the notion that ORM2 mediates the inhibitory effect of QUE on the viability of TNBC cells. Additionally, we conducted ORM2 overexpression experiments in MCF-10A and MCF-12A cells and determined the IC50 of QUE. As shown in Fig. 5B, the IC50 of the knockdown control in MCF-10A cells exceeded 100 µM, and the IC50 of ORM2-knockdown MCF-10A cells was also over 100 µM. Similarly, both the IC50 of the knockdown control and ORM2-knockdown MCF-12A cells exceeded 100 µM. Even after 72 hours of QUE exposure, no significant effect on cell viability was observed in MCF-10A and MCF-12A cells with ORM2 overexpression.
QUE increased survival and decreased expression of ORM2 in tumors in TNBC mice
To evaluate the therapeutic effect of QUE on TNBC, we utilized an orthotopic implantation TNBC mouse model. Initially, we conducted a survival analysis to determine the optimal treatment dose. The QUE-L groups showed minimal effect on survival compared to the control group. However, both QUE-M and QUE-H groups significantly improved the survival of TNBC mice, with the QUE-H group demonstrating the most favorable outcomes (Fig. 6A). Thus, we selected the QUE-H concentration for subsequent studies. To minimize survival bias, we chose day 30 as the endpoint when the majority of mice were still alive. We monitored the body weight of the mice every three days and measured tumor volume at the endpoint. The results indicated that the body weights of the mice did not exhibit significant changes. Although the QUE-treated groups showed slightly lower body weights compared to the control group, the difference was not statistically significant (Fig. 6B). However, the tumor volume in the QUE-treated groups was significantly reduced compared to the control group (Fig. 6C). Notably, the expression of ORM2 protein in the tumors was markedly decreased in the QUE-treated groups compared to the control group (Fig. 6D-E). Representative images of immunohistochemistry staining from tumors in the mice were presented in Fig. 6F.
To investigate whether attenuating ORM2 activity influences the therapeutic effect of QUE in TNBC treatment in vivo, we employed the ORM2 knockdown cell line in the in vivo study. We compared the overall treatment effect (survival) of QUE between ORM2 knockdown cell animals and ORM2 wild-type animals. The results revealed a striking difference in treatment outcomes between the treatment and vehicle groups in the wild-type cell model. However, in the ORM2 knockdown group, the difference between the treatment and vehicle groups was less pronounced. It is important to note that this experiment was limited by the transfection duration of the cell lines. While the cells developed tumors, the ORM2 knockdown may have been rapidly recovered to the same level, thereby potentially accounting for the difference in treatment effect observed in the initial few days and the final results. Notably, the vehicle group in the ORM2 knockdown model exhibited similar survival rates as the QUE treatment group in the ORM2 wild-type model, particularly before day 50. These findings suggest that ORM2 knockdown has a comparable effect to QUE treatment, consistent with previous results indicating that QUE reduces the level of ORM2.
Potential mechanism of QUE on ORM2
To gain insights into the potential mechanisms underlying the effect of QUE on ORM2 function in TNBC, we conducted an enrichment analysis of the differentially expressed genes associated with high and low ORM2 expression in TNBC, based on the TCGA cohort. The analysis revealed 244 up-regulated genes and 200 down-regulated genes that were associated with ORM2 expression (Fig. 7A). The up-regulated genes were found to be enriched in functions such as transmembrane receptor protein kinase activity, extracellular matrix structural constituent, and growth factor binding. Conversely, the down-regulated genes were associated with functions including interleukin-15 receptor activity, CXCR3 chemokine receptor binding, TAP binding, deoxycytidine deaminase activity, and peptide antigen binding (Fig. 7B).
Based on these differentially expressed genes, we constructed protein-protein interaction (PPI) networks. The PPI network for the up-regulated genes consisted of 242 nodes and 270 edges, with an average node degree of 2.23 and an average local clustering coefficient of 0.325. The PPI network for the down-regulated genes comprised 197 nodes and 22 edges, with an average node degree of 0.223 and an average local clustering coefficient of 0.129. The PPI enrichment analysis indicated a highly significant p-value of < 1.0e-16 for the up-regulated genes network and a significant p-value of < 9.22e-06 for the down-regulated genes network (Fig. 7C-D). These findings provide a comprehensive view of the potential interactions and functional associations of genes differentially expressed in relation to ORM2 in TNBC.