3.1. The expression levels of DYNLT1 are significantly upregulated in Luminal B breast cancer
The RNA seq transcriptomic data from the Breast Cancer Gene-Expression Miner v5.0 (bc-GenExMiner v5.0) database [20] was analyzed to determine the DYNLT1 expression in breast cancer. Notably, we found that DYNLT1 presented a remarkable upward trend in healthy, tumor-adjacent, and tumor tissues, revealing a significant connection between DYNLT1 expression and tumorigenesis of breast cancer (Fig. 1A). In comparison to the basal-like, HER2-enriched, Luminal A, and normal breast-like subtypes, the Luminal B subtype had considerably upregulated levels of DYNLT1 (Fig. 1B). The expression of DYNLT1 was further confirmed through IHC analysis of 62 paraffin-embedded breast cancer specimens using ERα, PR, HER2, and Ki67 statuses for tumor classification [11]. Consistently, DYNLT1 was substantially upregulated in Luminal B breast cancer tissues when compared to Luminal A, HER2-enriched, and triple-negative breast cancer (TNBC) subtypes (Fig. 1C-E). Collectively, our results demonstrated an upregulation of DYNLT1 levels specifically in Luminal B breast cancer.
3.2. The high expression of DYNLT1 is indicative of an unfavorable clinical recurrence-free survival (RFS) outcome in Luminal B breast cancer patients undergoing tamoxifen therapy
To determine how the DYNLT1 expression levels were related to the survival of breast cancer patients, the online tool, namely the Kaplan–Meier plotter (KM) [12], was adopted to evaluate the clinical significance of DYNLT1 in breast cancer. Notably, it was revealed that DYNLT1 upregulation was correlated with dismal RFS in all breast cancer cases (Hazard Ratio (HR) = 1.24, P = 0.00023) (Supplemental Fig. 1A). Next, we explored the RFS in breast cancer patients stratified with distinct molecular subtypes. The KM analysis indicated that DYNLT1 expression upregulation was significantly related to unfavorable RFS in Luminal B (HR = 1.53, P = 0.00014), Luminal A (HR = 1.3, P = 0.002), and HER2-enriched (HR = 1.56, P = 0.014) breast cancer subtypes, but favorable in the basal-like subtype (HR = 0.76, P = 0.027) (Supplemental Fig. 1B).
Luminal breast cancer is treated with tamoxifen as a first-line drug, but the drug resistance is often inescapable. We further explored whether the expression of DYNLT1 is related to the effectiveness of endocrine therapy in Luminal breast cancer. Analysis of the KM plotter database indicated that higher DYNLT1 expression significantly linked to unfavorable RFS in all breast cancer cases with tamoxifen-only (HR = 1.47, P = 0.0085) or endocrine-include (HR = 1.29, P = 0.047) therapy, but not in endocrine-exclude cases (HR = 0.87, P = 0.083) (Fig. 2A). Interestingly and importantly, when we performed subtype-specific survival analysis, the DYNLT1 expression upregulation was significantly linked to dismal RFS in Luminal B with tamoxifen-only (HR = 2.23, P = 0.00084) or endocrine-include therapy (HR = 1.6, P = 0.0094), and Luminal A with tamoxifen-only (HR = 1.51, P = 0.048), but not in Luminal A with endocrine-include (HR = 1.3, P = 0.12) and other endocrine-exclude breast cancer cases (Luminal B, HR = 1.35, P = 0.064; Luminal A, HR = 0.84, P = 0.29) (Fig. 2B and C). Therefore, we conclude that DYNLT1 expression upregulation is correlated with the relapse of Luminal B breast cancer treated with tamoxifen, resulting in poorer clinical outcomes.
3.3. Elevated DYNLT1 expression induces resistance to tamoxifen both in vitro and in vivo
The DYNLT1 expression in the luminal breast cancer cell lines was analyzed to ascertain the possible involvement of DYNLT1 in tamoxifen resistance. Notably, DYNLT1 showed significantly upregulation in Luminal B breast cancer cell lines (ZR-75-30, BT-474, and MDA-MB-361) relative to Luminal A cells (MCF7, T47D, and HCC1428) (Fig. 3A). We searched for the IC50 value of tamoxifen in these cell lines from Genomics of Drug Sensitivity in Cancer (GDSC) [45] to investigate how DYNLT1 expression was correlated with drug sensitivity in breast cancer. Surprisingly, elevated DYNLT1 expression was associated with drug resistance to tamoxifen (R = 0.855, P = 0.030) (Fig. 3B). Subsequently, we analyzed RNA seq data from the Gene Expression Omnibus (GEO) dataset (GSE106681) [5], comparing the expression levels of DYNLT1 in MCF7 and its tamoxifen-sensitive but aggressive subline (WS8) as well as tamoxifen-resistant MCF7-WS8 subline (TAMR). Remarkably, our results revealed significantly elevated DYNLT1 levels in TAMR cells relative to those in WS8 or parental MCF7 cells, regardless of tamoxifen treatment (Fig. 3C).
Our tamoxifen-resistant MCF7-TMR cell line, which originated from the tamoxifen-sensitive Luminal A breast cancer cell line MCF7, was tested for DYNLT1 expression to provide additional confirmation of the aforementioned findings [46]. In comparison to parental MCF7 cells, MCF7-TMR cells exhibited a markedly elevated DYNLT1 expression level (Fig. 3D). We established ZR-75-30 and MCF7-TMR cell lines, stably transduced with a DYNLT1 shRNA, as well as MCF7, stably expressing DYNLT1, to study the biological function of DYNLT1 and ascertain its involvement in tamoxifen resistance progression (Fig. 3E). Then, tamoxifen’s cytotoxic activity was explored on breast cancer cell lines. DYNLT1 silencing markedly increased the sensitivity of tamoxifen in ZR-75-30 cells making IC50 shift from 168.655 to 29.376 or 31.989 µM and MCF7-TMR (IC50, 41.591 to 14.588 or 16.982) cells, while DYNLT1 overexpression decreased the sensitivity of tamoxifen in MCF7 (IC50, 11.324 to 35.237) (Fig. 3F). The results show that the DYNLT1 overexpression induces a tamoxifen-resistant phenotype in luminal breast cancer cells.
Thereafter, the impact of DYNLT1 overexpression on tamoxifen resistance in breast cancer was assessed in vivo. One week following the implantation of E2 pellets, female mice had their mammary fat pads injected with either control cells or MCF7/DYNLT1 cells, and the first dosage of tamoxifen was administered after the tumor volume was about 200 mm3. The anti-tumor impact of tamoxifen therapy was considerably diminished in the MCF7 xenografts when DYNLT1 was overexpressed, as illustrated in Fig. 4A-C. Additionally, the apoptotic rates were lower and Ki67 levels were greater in tamoxifen-treated MCF7/DYNLT1 xenografts in comparison to control tumors (Fig. 4D and E), which supports the hypothesis that upregulation of DYNLT1 confers tamoxifen-resistant tumorigenicity to Luminal B breast cancer cells. Furthermore, Mice were administered injections with either control cells or MCF7-TMR/DYNLT1-Ri cells into their mammary fat pads, and treatment with tamoxifen commenced when the lesions had attained an approximate measurement of 200 mm3. Notably, Fig. 4F and G demonstrates that DYNLT1 knockdown in the MCF7-TMR xenografts without tamoxifen treatment significantly reduced the tumor size and weight, while co-therapy with tamoxifen and DYNLT1 shRNA more effectively blocked the tumorigenic potential of MCF7-TMR cells. In comparison to control tumors, tamoxifen-treated MCF7-TMR/DYNLT1-Ri xenografts consistently exhibited a more pronounced decrease in Ki67 levels and elevated apoptotic rates (Fig. 4H and I). These results suggest that silencing DYNLT1 in tamoxifen-resistant breast cancer cells makes them more responsive to the drug.
3.4. Overexpression of DYNLT1 facilitates ligand-independent activation of the ER signaling pathway
We first used GPS-pro to determine proteins potentially interacting with DYNLT1 to probe the molecular mechanisms via which DYNLT1 leads to tamoxifen resistance (Fig. 5A). The ingenuity pathway analysis g:Profiler [22] was utilized to analyze DYNLT1-interacting proteins enriched pathways. The results revealed that DYNLT1-interacting proteins exhibited significant enrichment in Gene Ontology Molecular Functions, including binding to dynein intermediate/heavy chain, identical protein binding, and microtubule motor activity (Fig. 5B). In terms of Biological Processes, enrichment was observed in microtubule-based processes, cell division, and the cell cycle. Additionally, these proteins were found to be localized in Cellular Components such as the cytoplasmic dynein complex, microtubules, and the nucleoplasm. These biological functions are recognized to have an essential function in the development of resistance to cancer therapies. Subsequently, Gene Set Enrichment Analysis (GSEA) was executed utilizing the online Webgestalt tool [26] to analyze genes coexpressed with DYNLT1 from cBioPortal [10], using data from The Cancer Genome Atlas (TCGA) Breast invasive carcinoma (BRCA). We discovered that hormone-mediated signaling pathway activity correlates positively with DYNLT1 expression (Fig. 5C), indicating that DYNLT1 upregulation may have a key function in ER signaling pathway activation.
We designed an estrogen-responsive element (ERE)-luciferase reporter to more thoroughly investigate whether DYNLT1 controls ERα transcriptional activity. Notably, overexpression of DYNLT1 in MCF7 cells ligand-dependently enhanced ERα transcriptional activity, whereas silencing DYNLT1 restored MCF7-TMR cell sensitivity to estrogen and tamoxifen (Fig. 5D). Additionally, the expression of trefoil factor 1 (TFF1) and progesterone receptor (PGR), known as classical ERα target genes, exhibited a significant increase in DYNLT1-overexpressing MCF7 cells, which were reduced in DYNLT1-silenced MCF7-TMR cells without estrogen or tamoxifen treatment (Fig. 5E). Although DYNLT1-overexpressed MCF7 cells were unaffected by tamoxifen treatment, DYNLT1-silenced MCF7-TMR cells showed a stronger inhibition of gene expression (Fig. 5E). The findings confirm that the ER-dependent transcriptional network is functionally regulated by DYNLT1.
3.5 ELAVL1-mediated m6A modification contributes to the upregulation of DYNLT1 expression and facilitates the acquisition of tamoxifen resistance phenotype
The most prevalent type of RNA epigenetic modification in higher eukaryotic mRNAs, m6A, is demonstrated to be critically involved in several types of biological processes [7, 13, 15, 49]. Firstly, SRAMP [50] was employed to predict the occurrence of m6A modification sites on the DYNLT1 RNA sequences. Two m6A sites with high confidence and one m6A site with moderate confidence across the mRNA transcript of DYNLT1 were identified (Fig. 6A). According to the m6A RIP test, the m6A level was higher in the ZR-75-30 and MCF7-TMR cells in comparison to MCF7, indicating that m6A may be involved in DYNLT1 upregulation (Fig. 6B). To find out which m6A regulators are implicated in regulating DYNLT1 expression, the link between the expression levels of writers, erasers, and readers (WERs) of m6A and DYNLT1 was further assessed using cBioPortal in the breast cancer tissue. The results revealed that the ELAVL1, WTAP, METTL5, and HNRNPC expression was positively linked with the DYNLT1 expression (P < 0.001, Spearman R > 0.3) (Fig. 6C, and Supplemental Fig. 2A). We then used ENCORI [24] to predict the interaction between DYNLT1 mRNA and the corresponding RNA binding proteins (RBPs), showing that DYNLT1 mRNA had a possible interaction with ELAVL1 and HNRNPC (Supplemental Fig. 2B). Further, we analyzed the iCLIP-seq data conducted in MCF7 cells from the ENCORI database and identified that the ELAVL1 binding site located at chr6:158636709–158636748[-] on Exon5/5 of DYNLT1 (NM_006519), which precisely corresponded to the #3 m6A site predicted by SRAMP (Fig. 6D). As expected, RIP analysis demonstrated that DYNLT1 mRNA was highly enriched in ELAVL1-immunoprecipitated complexes, compared to the IgG control group (Fig. 6E).
Reports indicate that ELAVL1, a critical component of the N6-methyltransferase complex, improves mRNA stability in breast cancer cells that are insensitive to tamoxifen [39]. Next, we used siRNA to target ELAVL1 in breast cancer cells. In both the ZR-75-30 and MCF7-TMR cell lines, we observed a significant drop in m6A and DYNLT1 mRNA levels after ELAVL1 silencing (Fig. 6F and G). After that, we investigated how ELAVL1 siRNA affected the stability of DYNLT1 mRNA and discovered that breast cancer cells treated with this siRNA had a much shorter half-life of DYNLT1 mRNA (Fig. 6H and I), supporting the notion that ELAVL1 controlled DYNLT1 expression by influencing the stability of DYNLT1 mRNA. Overall, the above evidence confirms that silencing of ELAVL1 attenuates the m6A modification of DYNLT1 mRNA and subsequently reduces the expression of DYNLT1.
Furthermore, KM analysis indicated that ELAVL1 expression upregulation significantly related to worse RFS in all or Luminal B breast cancer cases with endocrine-include therapy (all, HR = 1.92, P = 3.2e-09; Luminal B, HR = 1.75, P = 0.00083), but not in endocrine-exclude cases (all, HR = 0.86, P = 0.081; Luminal B, HR = 1.24, P = 0.15) (Supplemental Fig. 3A and B). In Luminal B breast cancer with tamoxifen-only therapy, although the trend of the survival curve is consistent with the above conclusions, the sample size was not sufficient to attain statistical significance (HR = 1.61, P = 0.065) (Supplemental Fig. 3B). Thus, high levels of ELAVL1 predict poor clinical RFS in endocrine-treated breast cancer patients.
Moreover, we investigated how ELAVL1 affected tamoxifen resistance in parental MCF7 and MCF7-TMR breast cancer cell lines. Notably, ELAVL1 expression upregulation significantly attenuated the anti-tumor effect of tamoxifen treatment in MCF7 cells, while DYNLT1 silencing resensitized ELAVL1 overexpressing-MCF7 cells to tamoxifen treatment (Fig. 6J). Conversely, the knockdown of ELAVL1 restored sensitivity to tamoxifen treatment in MCF7-TMR cells, whereas DYNLT1 expression upregulation significantly abrogated the anti-tumor properties of ELAVL1 silencing on MCF7-TMR cells under tamoxifen therapy (Fig. 6K). Collectively, these data suggest that upregulation of ELAVL1 facilitates m6A modification in DYNLT1 mRNA and subsequently augments DYNLT1 expression, thereby promoting tamoxifen resistance in Luminal B breast cancer.