Kinase shRNA screening reveals that TAOK3 enhances microtubule-targeted drug resistance of breast cancer cells via the NF-κB signaling pathway
Cell Communication and Signaling volume 18, Article number: 164 (2020)
Chemotherapy is currently one of the most effective treatments for advanced breast cancer. Anti-microtubule agents, including taxanes, eribulin and vinca-alkaloids are one of the primary major anti-breast cancer chemotherapies; however, chemoresistance remains a problem that is difficult to solve. We aimed to discover novel candidate protein targets to combat chemoresistance in breast cancer.
A lentiviral shRNA-based high-throughput screening platform was designed and developed to screen the global kinome to find new therapeutic targets in paclitaxel-resistant breast cancer cells. The phenotypes were confirmed with alternative expression in vitro and in vivo. Molecular mechanisms were investigated using global phosphoprotein arrays and expression microarrays. Global microarray analysis was performed to determine TAOK3 and genes that induced paclitaxel resistance.
A serine/threonine kinase gene, TAOK3, was identified from 724 screened kinase genes. TAOK3 shRNA exhibited the most significant reduction in IC50 values in response to paclitaxel treatment. Ectopic downregulation of TAOK3 resulted in paclitaxel-resistant breast cancer cells sensitize to paclitaxel treatment in vitro and in vivo. The expression of TAOK3 also was correlated to sensitivity to two other anti-microtubule drugs, eribulin and vinorelbine. Our TAOK3-modulated microarray analysis indicated that NF-κB signaling played a major upstream regulation role. TAOK3 inhibitor, CP43, and shRNA of NF-κB both reduced the paclitaxel resistance in TAOK3 overexpressed cells. In clinical microarray databases, high TAOK3 expressed breast cancer patients had poorer prognoses after adjuvant chemotherapy.
Here we identified TAOK3 overexpression increased anti-microtubule drug resistance through upregulation of NF-κB signaling, which reduced cell death in breast cancer. Therefore, inhibition of the interaction between TAOK3 and NF-κB signaling may have therapeutic implications for breast cancer patients treated with anti-microtubule drugs.
Breast cancer is currently the most common cancer among women worldwide as well as the third leading cause of cancer deaths in the United States . For locally advanced breast cancer, either paclitaxel- or docetaxel-containing adjuvant regimens are commonly used in clinical practice, with significantly better prognoses than other cytotoxic agents . A new generation of anti-microtubule inhibitors, such as eribulin, has also shown survival benefits in refractory metastatic breast cancer, which indicates the important role of microtubule-targeted drugs in preventing breast cancer recurrence and controlling progression [3, 4]. However, taxane resistance eventually developes in approximately 90% of patients despite excellent initial therapeutic efficacy . Several mechanisms of taxane resistance has previously been described: overexpression of multidrug resistance protein (MDR) genes, class III β-tubulin, and epithelial-mesenchymal transition [6,7,8,9,10]. The development of new agents to overcome resistance to taxane or other microtubule-targeting drugs is, therefore, indispensable to the advancement of disease treatment .
A systemic biology approach to screening of drug sensitivity has been used for the past 10 years in cancer treatment [12,13,14,15]. In order to discover novel therapeutic agents, shRNA is particularly valuable in identifying the mechanism of action of a compound with new anticancer indications and identifying potential targets. Furthermore, a kinome-based shRNA library for screening kinase protein inhibitors is manageable; there are 724 kinase proteins, and most of them have been implicated in primary cancer processes such as activating metastasis, sustaining proliferation, resisting cell death, etc [16,17,18] Currently, more than 10,000 candidate compounds have been patented to inhibit kinase activation . Some kinase inhibitors have been approved for anti-cancer therapies, for example, BCR–ABL inhibitors have been approved for chronic myeloid leukemia , BRAF inhibitors have been approved for melanoma  and HER2 (also known as ERBB2) inhibitors have been approved breast cancer . Other kinase inhibitors, such as Janus kinase (JAK) inhibitor tofacitinib, have been approved by FDA for the treatment of rheumatoid arthritis and could potentially be repurposed as novel anticancer agents.
In the present study, we utilized 724 shRNAs to generate a kinase library that we used to screen for novel genes responsible for taxane resistance in breast cancer. Knock down of a top ranking gene, TAOK3, overcome taxane resistance in breast cancer both in vitro and in vivo. Further pathway analysis revealed that TAOK3 may activate NF-κB signaling and blocking TAOK3-NF-κB signaling reduced the paclitaxel-resistance. Furthermore, clinical data were correlated with poor prognosis in breast cancer patients with high TAOK3 expression who accepted adjuvant therapy. This study reveals a new potential target for anticancer kinase inhibitors and future breast cancer treatments.
Cell lines and media
The 293 T, BT-20, BT-483, Hs578T, and MB 157 cell lines were cultured in Dulbecco’s modified Eagle’s medium (DMEM, Gibco, USA) supplemented with 10 mM L-glutamine and 10% fetal bovine serum. The AU565, Hcc38, Hcc70, Hcc1143, Hcc1937, Hcc1806, and T-47D cell lines were cultured in RPMI-1640 (Gibco, USA) medium supplemented with 10 mM L-glutamine and 10% fetal bovine serum. SKBR3 cells were grown in McCoy’s 5A medium (Gibco, USA) supplemented with 10 mM L-glutamine and 10% fetal bovine serum. MCF-7 cells were grown in MEM with 10% fetal bovine serum, 0.01 mg/mL human recombinant insulin and 10 mM L-glutamine. The cells were maintained in an incubator at 37 °C in 5% CO2. MDA-MB-231, MDA-MB-453, and MDA-MB-468 cells were cultured in L-15 medium supplemented with 10 mM L-glutamine and 10% fetal bovine serum. The cells were incubated in a humidified 37 °C incubator without extra CO2 supply. All cell lines were purchased from American Type Culture Collection (ATCC).
shRNA lentivirus production on solid-phase transfection platform
pGIPZ-shRNA glycerol stocks were purchased from Open Biosystems. DNA extraction was performed in a 96-well plate using a semi-automated Biomek-FX liquid handler and a ChargeSwitch nucleic acid purification kit (Invitrogen, USA). The negative control, which contains the pGIPZ vector backbone with a portion of scrambled shRNA (non-silencing gene), was included in each coated plate. A total of 150 ng of pGIPZ-shRNA plasmid and lentiviral packaging plasmids were used for plate coating on 96-well plates. Arrest-In (600 ng) (Open Biosystems, USA) was used as a transfection reagent, and the DNA-arrest-in complex was plated before the addition of 25% gelatin (Sigma, USA) . The plates with DNA complexes were dried at 65 °C for 90 min before storage at − 80 °C.
Cell viability assay
Paclitaxel (Sigma Aldrich, USA) was dissolved in DMSO and diluted in the growth medium. Eribulin (Halaven®, NerPharMa S.r.l., Italy) and vinorelbine (Navelbine®, Pierre Fabre, France) were diluted in the growth medium. High throughput cell viability assays were performed on 384-well white plates using Cell Titer Glo (Promega, USA). For general cytotoxicity assays, 2000 cells were seeded into each well 24 h before drug treatment. AlamarBlue (Invitrogen, USA) was added at 1/10 the total volume in a well and incubated for 2–4 h. Fluorescence intensity was recorded using an ELISA reader (PerkinElmer, USA). A cell imaging device, IncuCyte (IncuCyte, USA), was used to determine the cell confluence on 96-well plates over time.
Transfection, transduction, and cDNA cloning
All transfections for high-throughput screening were performed in DNA-coated 96-well plates with 50,000,293 T cells per well. The viral supernatants were collected at 48 and 72 h post transfection and stored at − 80 °C. Breast cancer cell line Au565 was transduced using the virus collected from transfection in 293 T cells. Transduction for screening was performed in 384-well plates with 8 μg/mL of polybrene (Sigma Aldrich, USA). Stable breast cancer cell lines with targeted shRNA and cDNA were generated with 10 μg/mL puromycin and blasticidin, respectively (InvivoGen, Hong Kong). TAOK3 cDNA was cloned from an ORF clone and sub-cloned into pLenti6.3 Gateway vector using Gateway cloning systems according to the manufacturer’s protocol (Invitrogen, USA).
RNA extraction and real-time quantitative PCR
Total RNA was extracted using Tri-reagent (Invitrogen, USA) and chloroform. The cDNA was synthesized by reverse transcriptase (Stratagene, USA) at 42 °C. Real-time PCR was performed using SyBr Green (Fermentas, Canada) and specific TAOK3 primers (5’gtgggcacaccttactggat3’ and 5’aacgttggggagtcattctg3’). Real-time PCR was performed in a BioRad 96-well real-time PCR detection system.
Total RNA was extracted with the RNeasy Mini kit (Qiagen, USA) and qualified with a Bioanalyzer (Agilent Technologies, USA). All samples were analyzed using Affymetrix GeneChip Human Genome U133 plus 2.0 arrays according to the manufacturer’s instructions. The data were normalized and analyzed by GeneSpring software (Agilent Tech., USA). Genes that changed more than threshold (1.5- and 2-fold) were sorted and further submitted to a computational simulation using Ingenuity Pathway Analysis (IPA, QIAGEN, USA) online tools to predict potential upstream regulators and the canonical pathways (pathways that represent common properties of a particular signaling module).
Protein extraction and Western blotting
Protein was extracted using RIPA buffer (20 mM Tris-HCl at pH 7.4, 150 mM NaCl, 0.5% Nonidet P-40, 1 mM EDTA, 50 μg/mL leupeptin, 30 μg/mL aprotinin, and 1 mM phenylmethylsulfonyl fluoride) containing proteinase inhibitors. Protein concentration was determined with the BCA kit (Thermo Scientific, Rockford, USA) using BSA as the standard. Approximately 20–100 μg of protein was loaded in an SDS-PAGE (TRIS-based), and blotting was performed on a nitrocellulose membrane (Amersham, Arlington Heights, IL, USA). Antibodies against TAOK3 (1:1000, #10158–2-AP, Proteintech, USA), phospho-p38 (1:1000, #4511, Cell signaling Tech.), p38 (1:2000, #9212, Cell Signaling Tech.), phospho-p65 (1:2000, #3033, Cell Signaling Tech.), p65 (1/2000, #4764, Cell Signaling Tech.), phospho-p53 (1:1000, #2521, Cell Signaling Tech.), p53 (1:500, #sc-126, Santa Cruz), caspase-3, (1:1000, #9662, Cell Signaling Tech.), PARP, (1:1000, #9542, Cell Signaling Tech.), β-actin (1:10000, Sigma) and α-tubulin (1:10000, Sigma) were diluted in blocking buffer. A secondary anti-mouse or anti-rabbit antibodies conjugated with HRP (Jackson ImmunoResearch Lab., USA) was used with 1:5000 dilution in blocking buffer. Visualization of the western blots was performed using the ECL Pro set (PerkinElmer) and X-ray radiography.
Caspase assays were performed on white 96-well plates according to the manufacturer’s protocol using caspase-3 Glo (Promega, USA). Approximately 20,000 cells were seeded onto the 96-well plate, and paclitaxel was added to the cells at 24 h before the caspase assay. The luciferase activity was measured using a Victor3 photometer, and the relative caspase activity was normalized with the corresponding AlamarBlue values.
Tissue slides were dewaxed felled by the detection of DNA cleavages using fluorescein-dUTP labeling with the enzyme terminal deoxynucleotidyl transferase (TdT) (In situ cell death detection kit, fluorescein, Roche, Switzerland). After labeling, slides were covered with mounting medium and sealed with nail polish. The images were scanned by the fluorescent slide scanner (Scanscope FL, Aperio, USA).
pGL4-NF-κB response vectors (Promega, USA) were transfected into cells for 48 h. The promoter activity was determined with 0.15 mg/mL luciferin and recorded by an IVIS Spectrum (PerkinElmer, USA). The activity values were normalized to the fluorescence intensity of AlamarBlue (Invitrogen, USA).
Xenograft mouse model
All animal experiments were approved by the Academia Sinica Institutional Animal Care and Unitization Committee (14–03-665). Immunodeficient (NOD-SCID) female mice (6–8 weeks old) were used. A total of 5 × 106 cells were re-suspended in 0.1 mL of PBS and injected subcutaneously. Tumor sizes were measured weekly and volume was calculated by 1/2 ab2 mm3. Taxol (paclitaxel, Bristol-Myers Squibb) in 0.2 mL PBS was administered via tail vein injection twice weekly. The Hs578T group was treated at approximately 6 mg/kg, the Hcc1806 high dose group was treated at approximately 3 mg/kg and the low dose group was treated at approximately 1 mg/kg.
All statistical calculations were performed in Excel, and the results are shown as the means ± standard deviations. The non-silencing control was used as the internal reference on each plate. Statistical significance was performed using a 2-samples t-test and a one-way ANOVA with a p-value less than 0.05.
Identification of paclitaxel-resistant gene targets in breast cancer cells
After evaluating the IC50 of paclitaxel in 15 breast cancer cell lines, we found that Au565 was the most paclitaxel-resistant cell line to paclitaxel (Fig. 1a, Table S1). For the screening method, we used a solid phase transfection system to produce shRNA lentivirus in 96 well plate. First, we determined the optimal conditions to achieve the highest transfection efficiency (Fig. 1b). We also evaluated the stability of the transfected, coating plates at different temperatures and time points. We found that 3 day of plate coating was optimal and that the efficiency was maintained for more than 30 days in the refrigerator (Fig. 1c). Next, we used 293 T cells for high-throughput lentiviral screening. In the screening, we prepared the shRNA lentivirus in 96-well plates containing 724 kinase genes; altogether, we created 1646 shRNA clones (Fig. 1d). We compare the results of the “before” and “after” paclitaxel treatment groups (the altered ratio distributions are shown in Fig. 2a). We selected the top 50 candidate genes whose knockdown induced paclitaxel sensitivities (listed in Table S2). Many of the identified kinases belonged to the MAPK, PI3K-AKT, or NF-κB signaling pathway. We selected some novel kinase genes for second round validation, and our results showed that the knockdown of TAOK3 created the highest sensitivity to paclitaxel (Fig. 2b). We found the TAOK3 protein expression level was positive correlated to paclitaxel IC50 values in breast cancer cell lines (Fig. 2c and d). Taken together, these facts indicate that TAOK3 may be associated with paclitaxel resistance.
Knockdown of TAOK3 increased chemosensitivity in breast cancer cells
First, we determined the endogenous expression levels of TAOK3 in 12 breast cancer cell lines. The expression of TAOK3 in normal breast epithelial cells (H184B5F5/M10 and MCF-10A) was significantly lower than in breast cancer cells (Fig. 2c). To determine the effect of TAOK3 on paclitaxel sensitivity, we knocked down the TAOK3 expression in the high expression cell line, AU565, Hcc1806, and, SKBR3, with two shRNA clones (Figure S1A). Cells with TAOK3 shRNA treatment showed enhanced sensitivity to paclitaxel (Fig. 3a and b). SKBR3 cells, which had moderate TAOK3 expression, had an increased enhancement of paclitaxel sensitivity (Figure S1B). In addition, we overexpressed TAOK3 in the low endogenous TAOK3 cells, Hs578T and MB157. The overexpression of TAOK3 conferred the sensitive breast cancer cells with an increased resistance to paclitaxel (Fig. 3c and d). The overexpression of TAOK3 increased the original IC50 of paclitaxel by 8.9-fold and 11.1-fold in Hs578T and MB157 cells, respectively. Next, we treated cells with the anti-microtubule drugs eribulin and vinorelbine . Downregulated TAOK3 expression in breast cancer cells enhanced the response to the drugs in non-silenced cells (Fig. 3e and g). In TAOK3 overexpression cells, cells exhibited higher drug resistance than the control group (Fig. 3f and h). However, these changes in drug response were not observed with the drugs cisplatin and doxorubicin, which interact with DNA (Figure S2). The chemoresistance effect of TAOK3 was specific to microtubule-targeted drugs. Notably, the TAOK3-related growth inhibition or enhancement was paclitaxel-dependent, as no significant difference in cell growth were observed in cells with either TAOK3 shRNA or cDNA when compared to the control cells (Figure S3).
TAOK3 knockdown increased taxane-induced cell apoptosis
To understand the mechanism of TAOK3-dependent drug resistance, we measured the cell death by detecting caspase-3 activity, sub-G1 percentage, and PARP cleavage after treatment with paclitaxel in TAOK3-manipulated cells. The results showed that caspase-3/7 was prominently activated upon treatment with an IC50 dose of paclitaxel in AU565 and MB157 cells after the depletion of TAOK3 (Fig. 4a). There were no significant differences in caspase-3/7 activation observed at 24 h in control cells treated with paclitaxel or in untreated cells. Conversely, a significant difference was observed in shTAOK3-modified cell lines treated with paclitaxel. The overexpression of TAOK3 conferred MB157 and BT483 cells with more resistance to paclitaxel treatment, and no significant difference in caspase-3/7 activity was observed in the TAOK3-overexpressed cells (Fig. 4b). For the cell cycle analysis, we harvested MB157 cells after treating with 0, 0.33, 1, 3.33, 10 and 33.33 nM paclitaxel for 24 h. At 3.33 nM treatment, the highest sub-G1 percentage was observed in the TAOK3 knockdown groups. In contrast, the sub-G1 percentage in the cell population with TAOK3 overexpressed did not dramatically increased after treatment with paclitaxel (Fig. 4c). However, if the dose of paclitaxel was too high, the sub-G1 percentage decreased the cells became arrested at the G2 phase. We also determined the cleavage of caspase-3 and PARP in apoptotic Hs578T cells with TAOK3 overexpression compared to the control with 0.1, 1 and 10 μM paclitaxel treatment after 24 h incubation. The results showed that TAOK3 overexpression reduced the cleavage of both caspase-3 and PARP (Fig. 4d).
Overexpression of TAOK3 increased chemoresistance in breast cancer cells in vivo
To access the effect of paclitaxel treatment on different TAOK3 expression levels in vivo, we determined the paclitaxel response effects in Hcc1806 and Hs578T using an in vivo subcutaneous xenograft tumor model. The experimental mice were injected with the cells at two sides: one side was injected with control cells and the other side was injected with the TAKO3-overexpressed clone. When the tumor sizes reached approximately 0.5 cm in diameter, paclitaxel (TAXOL® reagent) was intravenously injected into mice. The experimental mice received the drug twice a week until the tumors were bigger than 1500 mm3. Without paclitaxel treatment, there was no significant difference between the control and TAOK3-overexpressed groups. However, the paclitaxel treatment significantly inhibited the tumor growth of control cells (73%), compared to the TAOK3-overexpressed group (17%) at week 7 (Fig. 5a). A comparison of the tumor weight between the 4 groups revealed significant inhibition after paclitaxel treatment in the control group (Fig. 5b). In addition, tumor sections of similar size were selected for TUNEL staining to evaluate cell death after paclitaxel treatment. The results showed that the paclitaxel-treated group had a smaller percentage of positive staining in non-necrosis regions of Hs578T-TAOK3 tumor than Hs578T-VC tumor (Figs. 5c, S4). In the Hcc1806-shTAOK3 groups, the data showed dose-dependent enhanced paclitaxel sensitivity (Fig. 5d). High doses of paclitaxel decreased tumor weights more significantly than low dose (shown in Fig. 5e). The TUNEL staining also revealed that cell death was higher in the Hcc1806-shTAOK3 tumors than in the Hcc1806-NS tumors. (Figs. 5f, S4).
TAOK3 regulated paclitaxel resistance via the NF-κB/PTGS2 signaling pathway
We performed RNA microarray analysis to analyze 522 genes from Hs578T cells with TAOK3 overexpression and Hcc1806 cells with TAOK3 knockdown (Fig. 6a). After the IPA upstream analysis, both sets of cells showed enrichment of SP1, TWIST1, CHUK and NF-κB signaling, indicating that these pathways are involved in upstream regulatory functions (Fig. 6b). We also stained a phosphoprotein array to determine the changes of various phosphokinase (Figure S5A). After semi-quantitation, phosphorylation of p53 showed the largest enhanced in overexpression condition (Figure S5B). However, NF-κB-associated kinases were not presented in the phosphoprotein array. Genes associated with RELA, NF-κB and p53 pathways and genes that changed more than 2-fold in both cell sets were showed in Fig. 6c (Figure S4C and D). In this signaling pathway, we found that PLA2G4A, PTGS2, and PDE4B (shown as red diamonds) exhibited at least 10 times increase with TAOK3 overexpression. We confirmed NF-κB activity in these cells by using a NF-κB promoter assay. Our result showed that the expression of TAOK3 and NF-κB activity were positively correlated (shown in Fig. 6d). The RNA expression levels of PTGS2, PLA2G4A, and PDE4B were determined with real-time PCR (Fig. 6e). PTGS2 is known as an NF-κB-directed target. These phenotypes were validated with an analysis of endogenous protein expression (Fig. 6e). Phospho-p65 was found to be upregulated in TAOK3-overexpressed cells; however, it was undetectable in Hcc1806-NS and Hcc1806-shTAOK3 cells. The results of the PTGS2 analysis similar to the p-p65 results. Phospho-p53 was found to be upregulated and downregulated in TAOK3-modified cells. Hcc1806 had two-base-pair insertion at codon 256; thus, the protein was the wrong molecular weight, but phospho-p53 was found with a very low expression level. In addition, a mutant p53 protein (V157F) was also found in Hs578T cells. Furthermore, other cells such as AU545 (wildtype p53), SKBR3 (L175H), and MB157 (p53, del 26 bp at 261) exhibited responses to paclitaxel. These results indicate that the interaction with TAOK3 and p53 may not be critical to the cell death induced by paclitaxel. Taken together, we conclude that TAOK3 reduced paclitaxel cytotoxicity through the activity of the NF-κB signaling pathway.
Reduce paclitaxel resistance with a TAOK3 inhibitor and NF-κB shRNA in TAOK3-overexpressed cells
To overcome the TAOK3 effects on paclitaxel resistance, we tested a specific inhibitor of the TAOK3 family, CP43 , and used NF-κB shRNA on a TAOK3 overexpression clone. In Hs578T cells, we found that CP43 increased the percentage of mitotic cells (histone H3-pS10-positive) in all clones; however, in TAOK3-overexpressed cells, the percentage increase was dose-dependent (Fig. 7a). When we knocked down TAOK3, the response to CP43 were reduced (Fig. 7b). We also found that when we combined CP43 and paclitaxel, CP43 reduced paclitaxel resistance in TAOK3-overexpressed cells but only slightly reduced paclitaxel resistance in vector control (Fig. 7c). We further knocked down NF-κB in Hs578T cells, and we did not observe a dramatic change in mitotic cells (Figure S6A). We also did not observe the synergistic effects of NF-κB shRNA and paclitaxel in the vector control group (Figure S6B). Conversely, we did notice a reduction in paclitaxel resistance after NF-κB shRNA knockdown in the TOAK3-overexpressed cells (Fig. 7d, clone T3-shN2). Such results may indicate that inhibition of TAOK3-NF-κB signaling is a potential treatment on reducing paclitaxel-resistance in breast cancer cells.
TAOK3 overexpression correlated to poor prognosis of breast cancer patients who accepted adjuvant chemotherapy
Using the bioinformatics tool “Kaplan-Meier Plotter (https://kmplot.com/analysis/)” , which contains both publicly available array profiles and clinical data, we evaluated whether TAOK3 was associated with the prognoses for breast cancer patients. Because taxane-based adjuvant chemotherapy is used in most cases of high-risk breast cancer, we focused on patients with subsequent adjuvant chemotherapy (except endocrine therapy); we found that higher TAOK3 expression was significantly correlated to poor recurrence-free survival (HR = 1.7(1.2–2.41), p = 0.0024) (Fig. 8a) rather than post progression survival (Fig. 8b). This difference was not observed in patients who accepted endocrine-only adjuvant treatments (Fig. 8c). We further investigated possible correlation of TAOK3 with other chemotherapy drugs. Using the GSE16446 database, which contains data about breast cancer patients who accepted only epirubicin as neoadjuvant therapy, we found no significant difference in recurrence-free survival between high vs. low TAOK3-expressed patients (HR = 0.61(0.27–1.38), p = 0.23) (Fig. 8d). The above results provide more evidence that TAOK3 is associated with breast cancer prognosis specifically after taxane treatment.
In this study, we established a simple and stable solid-phase transfection method to systematically produce many different shRNA lentiviruses. We prepared a kinase library contained 724 genes to screen for paclitaxel-resistance genes in the most resistant breast cancer cell line. The top 50 candidate genes included many genes related to paclitaxel resistance, such as AKT1, SRC, etc., but to date, many have not been mentioned in the literature [27, 28]. After a similarity score analysis (score > 0.5) using the DAVID bioinformatics database , we found two major subsets in seven related clusters. One group included ABL2, YES1, BMX, and LCK, which have been associated with paclitaxel resistance [30,31,32,33]. The other group included MINK1, TSSK2, MARK1, MARK2, WNK3, VRK3, and TAOK3; however, this group has rarely been mentioned. The candidates that are showed in Fig. 2b are the genes that have rarely been reported to be involved in taxane resistance. TAOK3 shRNA initiated the most significant decline of paclitaxel IC50 values. We also found that the expression of TAOK3 in breast cancer cell lines was generally higher than that in normal cells (Fig. 2c). The expression of TAOK3 was positive correlative to the IC50 of paclitaxel in breast cancer cell lines (Fig. 2d). This phenomenon was also observed when different sensitivity to paclitaxel were caused by adjusting the expression of TAOK3 in different cell lines (Fig. 3a-d). In addition, we also found that TAOK3 expression was related to other anti-microtubule agents but relatively unlinked to the action of DNA damage agents (Figs. 3e-h, S2). The overexpression of TAOK3 rendered cells less likely to die due to paclitaxel, and in vivo experiments showed a similar pattern. In addition, we also found that the necrosis area was decreased when tumor with TAOK3 overexpression but without paclitaxel treatment (Figure S7). In vitro TAOK3 overexpression reduced the cleavage of caspase-3 and decreased DNA breakdown induced by paclitaxel treatment in vivo.
TAOK3 (thousand and one kinase 3) belongs to the STE20 family of serine/threonine kinase. Unlike other STE20-like kinase that indirectly enhanced the activity of JNK, ERK, and p38MAPK signaling, TAOK3 inhibited the activity of JNK but does not affect other MAPK signaling [34,35,36,37]. There are three known TAO kinases: TAOK1 ~ 3, which are activated by stress; for example, TAOK2 was found to render cells resistant to irradiation by enhancing the capability of initiating DNA damage-induced G2/M arrest . TAOK3 was also found to directly phosphorylated the LATS1/2 and MST1/2 and act as a negative regulator to the effectors YAP/TAZ in the Hippo pathway . Paclitaxel sensitivity has been highly correlated to ERK and p38 MAPK signaling [40, 41]. However, the modulation of TOAK3 expression did not affect JNK, ERK and p38 MAPK signaling in our TAOK3 overexpression kinase protein array. The microarray analysis indicated complex downstream signaling interactions of RELA and NFKB1 that induced the phosphorylation of MAPK, NF-κB and PI3K-AKT signaling pathways, which are involved in the paclitaxel response [42,43,44]. These signaling pathways enhance cell survival ability and regulates the anti-apoptosis proteins. Our data shown that NF-κB signaling was the most upregulated. The activation of NF-κB signaling has been reported to play a critical role in paclitaxel resistance in ovarian, non-small cell lung cancer, and breast cancer cells [45,46,47]. In our study, the overexpression of TAOK3 increased phosphorylation of NF-κB. The NF-κB signaling has been found to participate in crosstalk with other signaling pathways including AKT and p53 [48,49,50]. One study on the interaction between YAP and NF-κB has suggested verteporfin as a potential compound to reverse paclitaxel resistance in lung cancer . Currently, there is no evidence that TAOK3 directly binds and activates the NF-κB protein; therefore, indirect modulation is more plausible. The links between TAOK3 and NF-κB still need to be explored.
Mitotic slippage is one of the primary mechanisms of action to paclitaxel . Initially, mitotic arrest protects cells from chromosome segregation and generation of aneuploid cells. Normally, cells that slip out of mitosis stop dividing, become senescent or die . However, some resistant cells can withstand mitotic catastrophe and survive. In a previous study, CP43 was found to inhibit TAOK1 and TAOK2, which caused an increase in the percentage of cells found in the mitotic phase via disrupting of the spindle assembly in nucleus . In cell cycle analysis, this effect of CP43 is similar to the effect of paclitaxel. However, the effect of CP43 on TAOK3 is unknown. In our microarray analysis. TAOK3 modulation had no effect on the expression of TAOK1 or TAOK2 (Table S3). However, we observed that the number of cells in mitotic changed after CP43 treatment. Thus, TAOK3 may encourage mitotic slippage to reduce cell death or senescence; therefore, it would require higher CP43 to reach the threshold. More cells died when paclitaxel and CP43 were combined together (Fig. 7c). In another experiment, we used shRNA NF-κB to block TAOK3-NF-κB signaling in TAOK3-overexpressed cells. The effects were different than those of CP43. The number of cells in the mitotic phase was slightly reduced when NF-κB expression was knockdown. The synergistic effect with paclitaxel was shown in TAOK3-overexpressed cells but not in controls. This could be because NF-κB signaling was not critical to mitotic slippage but still played a role in paclitaxel resistance in the TAOK3-overexpressed cells.
Interestingly, we also found that TAOK3 is associated with chemosensitivity to other anti-microtubule agents such as eribulin and vinorelbine (Fig. 3e to h). Prior clinical trials have shown that eribulin is effective after prior exposure to taxanes [3, 54]. Moreover, vinorelbine is commonly used in chemotherapy after taxane exposure [55, 56]. Our findings warrant further investigation into taxane-exposed tumors and any correlation with the effects of subsequent treatment with eribulin or vinorebine. Our research suggests that TAOK3 may play a role in determining the chemo-sensitivity of tumors treated with anti-microtubule agents. However, more mechanistic studies are needed.
In summary, we found that TAOK3 expression enhanced the paclitaxel resistance of breast cancer cells via the NF-κB signaling pathway. However, the mechanism of TAOK3-NF-κB-PTGS2 pathway remains unclear. Although there are no known specific inhibitors of TAOK3, there are several commercial NF-κB inhibitors, and the inhibition of NF-κB signaling could provide a putative resolution for TAOK3-associated anti-microtubule drug resistance. In the future, TAOK3 as a molecular target in cancer treatment should be evaluated.
In this study, we screened paclitaxel response-associated kinases and provided evidence that TAOK3 overexpression reduced the sensitivity to anti-microtubule drug in breast cancer cells and was correlated with poor outcomes in patients. TAOK3 to NF-κB signaling could be a new target for drug development and therapeutic strategies for breast cancer.
Availability of data and materials
Clinical sample analysis was from the Kaplan-Meier Plotter database. Please refer the caption Fig. 8.
Short hairpin RNA
Small interfering RNA
Nuclear factor kappa-light-chain-enhancer of activated B cells
Multidrug resistance protein
Breakpoint cluster region
Abelson murine leukemia viral oncogene
B-Raf murine sarcoma viral oncogene homolog B
Human epidermal growth factor receptor 2
Erb-b2 receptor tyrosine kinase 2
Enzyme-linked immunosorbent assay
Polymerase chain reaction
Ingenuity pathway analysis
Poly (ADP-ribose) polymerase
Mitogen-activated protein kinase
The half maximal inhibitory concentration
Terminal deoxynucleotidyl transferase dUTP nick end labeling
Specificity protein 1
Twist-related protein 1
Conserved helix-loop-helix ubiquitous kinase
Nuclear factor NF-kappa-B p65 subunit
Cytosolic phospholipase A2
Prostaglandin-endoperoxide synthase 2
cAMP-specific 3′,5′-cyclic phosphodiesterase 4B
Extracellular signal–regulated kinase
Mitogen-activated protein kinase kinase
c-Jun N-terminal kinase
The cancer genome atlas
Siegel RL, Miller KD, Jemal A. Cancer statistics, 2018. CA Cancer J Clin. 2018;68(1):7–30.
Gradishar WJ, Anderson BO, Balassanian R, et al. NCCN guidelines insights breast Cancer, version 1.2016. J Natl Compr Cancer Netw. 2015;13(12):1475–85.
Cortes J, O'Shaughnessy J, Loesch D, et al. Eribulin monotherapy versus treatment of physician's choice in patients with metastatic breast cancer (EMBRACE): a phase 3 open-label randomised study. Lancet. 2011;377(9769):914–23.
Kaufman PA, Awada A, Twelves C, et al. Phase III open-label randomized study of eribulin mesylate versus capecitabine in patients with locally advanced or metastatic breast cancer previously treated with an anthracycline and a taxane. J Clin Oncol. 2015;33(6):594–601.
Murray S, Briasoulis E, Linardou H, Bafaloukos D, Papadimitriou C. Taxane resistance in breast cancer: mechanisms, predictive biomarkers and circumvention strategies. Cancer Treat Rev. 2012;38(7):890–903.
Zelnak A. Overcoming taxane and anthracycline resistance. Breast J. 2010;16(3):309–12.
Yang Q, Huang J, Wu Q, et al. Acquisition of epithelial-mesenchymal transition is associated with Skp2 expression in paclitaxel-resistant breast cancer cells. Br J Cancer. 2014;110(8):1958–67.
Dai M, Zhang C, Ali A, et al. CDK4 regulates cancer stemness and is a novel therapeutic target for triple-negative breast cancer. Sci Rep. 2016;6:35383.
Asaduzzaman M, Constantinou S, Min H, et al. Tumour suppressor EP300, a modulator of paclitaxel resistance and stemness, is downregulated in metaplastic breast cancer. Breast Cancer Res Treat. 2017;163(3):461–74.
Stengel C, Newman SP, Leese MP, Potter BV, Reed MJ, Purohit A. Class III beta-tubulin expression and in vitro resistance to microtubule targeting agents. Br J Cancer. 2010;102(2):316–24.
Rivera E, Gomez H. Chemotherapy resistance in metastatic breast cancer: the evolving role of ixabepilone. Breast Cancer Res. 2010;12(Suppl 2):S2.
Swanton C, Szallasi Z, Brenton JD, Downward J. Functional genomic analysis of drug sensitivity pathways to guide adjuvant strategies in breast cancer. Breast Cancer Res. 2008;10(5):214.
van de Vijver MJ, He YD, van't Veer LJ, et al. A gene-expression signature as a predictor of survival in breast cancer. N Engl J Med. 2002;347(25):1999–2009.
van ’t Veer LJ, Dai H, van de Vijver MJ, et al. Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002;415(6871):530–6.
Etemadmoghadam D, deFazio A, Beroukhim R, et al. Integrated genome-wide DNA copy number and expression analysis identifies distinct mechanisms of primary chemoresistance in ovarian carcinomas. Clin Cancer Res. 2009;15(4):1417–27.
Miller SM, Goulet DR, Johnson GL. Targeting the breast Cancer Kinome. J Cell Physiol. 2017;232(1):53–60.
Donnella HJ, Webber JT, Levin RS, et al. Kinome rewiring reveals AURKA limits PI3K-pathway inhibitor efficacy in breast cancer. Nat Chem Biol. 2018;14(8):768–77.
Murray HC, Dun MD, Verrills NM. Harnessing the power of proteomics for identification of oncogenic, druggable signalling pathways in cancer. Expert Opin Drug Discov. 2017;12(5):431–47.
Akritopoulou-Zanze I, Hajduk PJ. Kinase-targeted libraries: the design and synthesis of novel, potent, and selective kinase inhibitors. Drug Discov Today. 2009;14(5–6):291–7.
Hunter T. Treatment for chronic myelogenous leukemia: the long road to imatinib. J Clin Invest. 2007;117(8):2036–43.
Young K, Minchom A, Larkin J. BRIM-1, −2 and −3 trials: improved survival with vemurafenib in metastatic melanoma patients with a BRAF(V600E) mutation. Future Oncol. 2012;8(5):499–507.
U.S. Food and Drug Administration. Breast cancer drug approved for new indication. Womens Health (Lond). 2010;6(2):173.
Chang FH, Lee CH, Chen MT, et al. Surfection: a new platform for transfected cell arrays. Nucleic Acids Res. 2004;32(3):e33.
Jordan MA, Wilson L. Microtubules as a target for anticancer drugs. Nat Rev Cancer. 2004;4(4):253–65.
Koo CY, Giacomini C, Reyes-Corral M, et al. Targeting TAO kinases using a new inhibitor compound delays mitosis and induces mitotic cell death in centrosome amplified breast Cancer cells. Mol Cancer Ther. 2017;16(11):2410–21.
Gyorffy B, Lanczky A, Eklund AC, et al. An online survival analysis tool to rapidly assess the effect of 22,277 genes on breast cancer prognosis using microarray data of 1,809 patients. Breast Cancer Res Treat. 2010;123(3):725–31.
Hawthorne VS, Huang WC, Neal CL, Tseng LM, Hung MC, Yu D. ErbB2-mediated Src and signal transducer and activator of transcription 3 activation leads to transcriptional up-regulation of p21Cip1 and chemoresistance in breast cancer cells. Mol Cancer Res. 2009;7(4):592–600.
Swanton C, Marani M, Pardo O, et al. Regulators of mitotic arrest and ceramide metabolism are determinants of sensitivity to paclitaxel and other chemotherapeutic drugs. Cancer Cell. 2007;11(6):498–512.
Huang da W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.
Chen L, Cao H, Feng Y. MiR-199a suppresses prostate cancer paclitaxel resistance by targeting YES1. World J Urol. 2018;36(3):357–65.
Uckun FM, Dibirdik I, Qazi S, et al. Anti-breast cancer activity of LFM-A13, a potent inhibitor of polo-like kinase (PLK). Bioorg Med Chem. 2007;15(2):800–14.
Duan Z, Zhang J, Ye S, et al. A-770041 reverses paclitaxel and doxorubicin resistance in osteosarcoma cells. BMC Cancer. 2014;14:681.
Katayama MLH, Vieira R, Andrade VP, et al. Stromal Cell Signature Associated with Response to Neoadjuvant Chemotherapy in Locally Advanced Breast Cancer. Cells. 2019;8:1566.
Tassi E, Biesova Z, Di Fiore PP, Gutkind JS, Wong WT. Human JIK, a novel member of the STE20 kinase family that inhibits JNK and is negatively regulated by epidermal growth factor. J Biol Chem. 1999;274(47):33287–95.
Chen Z, Raman M, Chen L, Lee SF, Gilman AG, Cobb MH. TAO (thousand-and-one amino acid) protein kinases mediate signaling from carbachol to p38 mitogen-activated protein kinase and ternary complex factors. J Biol Chem. 2003;278(25):22278–83.
Hutchison M, Berman KS, Cobb MH. Isolation of TAO1, a protein kinase that activates MEKs in stress-activated protein kinase cascades. J Biol Chem. 1998;273(44):28625–32.
Yustein JT, Xia L, Kahlenburg JM, Robinson D, Templeton D, Kung HJ. Comparative studies of a new subfamily of human Ste20-like kinases: homodimerization, subcellular localization, and selective activation of MKK3 and p38. Oncogene. 2003;22(40):6129–41.
Raman M, Earnest S, Zhang K, Zhao Y, Cobb MH. TAO kinases mediate activation of p38 in response to DNA damage. EMBO J. 2007;26(8):2005–14.
Plouffe SW, Meng Z, Lin KC, et al. Characterization of hippo pathway components by gene inactivation. Mol Cell. 2016;64(5):993–1008.
Bacus SS, Gudkov AV, Lowe M, et al. Taxol-induced apoptosis depends on MAP kinase pathways (ERK and p38) and is independent of p53. Oncogene. 2001;20(2):147–55.
McCubrey JA, Steelman LS, Abrams SL, et al. Roles of the RAF/MEK/ERK and PI3K/PTEN/AKT pathways in malignant transformation and drug resistance. Adv Enzym Regul. 2006;46:249–79.
Okano J, Rustgi AK. Paclitaxel induces prolonged activation of the Ras/MEK/ERK pathway independently of activating the programmed cell death machinery. J Biol Chem. 2001;276(22):19555–64.
Rajput S, Volk-Draper LD, Ran S. TLR4 is a novel determinant of the response to paclitaxel in breast cancer. Mol Cancer Ther. 2013;12(8):1676–87.
Huang S, Wang D, Zhang S, et al. Tunicamycin potentiates paclitaxel-induced apoptosis through inhibition of PI3K/AKT and MAPK pathways in breast cancer. Cancer Chemother Pharmacol. 2017;80:685–96.
Sun NK, Huang SL, Chang TC, Chao CC. TLR4 and NFkappaB signaling is critical for taxol resistance in ovarian carcinoma cells. J Cell Physiol. 2018;233(3):2489–501.
Xiang F, Ni Z, Zhan Y, et al. Increased expression of MyD88 and association with paclitaxel resistance in breast cancer. Tumour Biol. 2016;37(5):6017–25.
Chen X, Sun X, Guan J, et al. Rsf-1 influences the sensitivity of non-small cell lung Cancer to paclitaxel by regulating NF-kappaB pathway and its downstream proteins. Cell Physiol Biochem. 2017;44(6):2322–36.
Hoesel B, Schmid JA. The complexity of NF-kappaB signaling in inflammation and cancer. Mol Cancer. 2013;12:86.
Yun M, Lee D, Park MN, et al. Cinnamaldehyde derivative (CB-PIC) sensitizes chemo-resistant cancer cells to drug-induced apoptosis via suppression of MDR1 and its upstream STAT3 and AKT signalling. Cell Physiol Biochem. 2015;35(5):1821–30.
D'Antona L, Amato R, Talarico C, et al. SI113, a specific inhibitor of the Sgk1 kinase activity that counteracts cancer cell proliferation. Cell Physiol Biochem. 2015;35(5):2006–18.
Pan W, Wang Q, Zhang Y, et al. Verteporfin can reverse the paclitaxel resistance induced by YAP over-expression in HCT-8/T cells without Photoactivation through inhibiting YAP expression. Cell Physiol Biochem. 2016;39(2):481–90.
Visconti R, Grieco D. Fighting tubulin-targeting anticancer drug toxicity and resistance. Endocr Relat Cancer. 2017;24(9):T107–T17.
Gascoigne KE, Taylor SS. How do anti-mitotic drugs kill cancer cells? J Cell Sci. 2009;122(Pt 15):2579–85.
Lorusso V, Cinieri S, Latorre A, et al. Efficacy and safety of eribulin in taxane-refractory patients in the 'real world'. Future Oncol. 2017;13(11):971–8.
Stravodimou A, Voutsadakis IA. A systematic review and meta-analysis of the combination of Vinorelbine and Lapatinib in patients with Her2-positive metastatic breast Cancer. Anticancer Res. 2019;39(7):3295–301.
Blancas I, Aguirre E, Morales S, et al. Real-world data on the efficacy and safety of weekly oral vinorelbine in breast cancer patients previously treated with anthracycline or taxane-based regimens. Clin Transl Oncol. 2019;21(4):459–66.
We like to thank Miss Yi-Cian Lin for her assistance.
This research was supported by Academia Sinica and Ministry of Science and Technology [AS-SUMMIT-108], [AS-SUMMIT-109], [MOST-108-3114-Y-001-002], and [ASKPQ-109-BioMed] to Michael Hsiao.
Ethics approval and consent to participate
Clinical samples were written informed consent in this study.
All animal experiment protocol (14–03-665) was approved by the Academia Sinica Institutional Animal Care and Unitization Committee.
Consent for publication
The Authors claim no competing interests in this study.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
IC50 of paclitaxel among breast cancer cell lines. Supplement Table 2. Top 50 candidate list from kenome shRNA screening. Supplement Table 3. The fold change of TAOK family in microarray with TAOK3 modulation. Supplement Figure 1. Protein expression TAOK3 and cell viability of paclitaxel in SKBR3 (A) The effects of TAOK3 shRNA in Au565, Hcc1806 and SKBR3 and overexpressed TAOK3 in MB157 and Hs578t. (B) Paclitaxel sensitivity changes of shTAOK3 SKBR3 cells with paclitaxel treatment. Supplement Figure 2. The effect of cisplatin and doxorubicin with alternative TAOK3 expression. (A) Cell viability assay of cisplatin among Hcc1806-NS, Hcc1806-shTAOK3–1 and Hcc1806-shTAOK3–2. (B) Cell viability assay of cisplatin in Hs578t-VC and Hs578t-TAOK3. (C) Cell viability assay of doxorubicin among Hcc1806-NS, Hcc1806-shTAOK3–1 and Hcc1806-shTAOK3–2. (D) Cell viability assay of doxorubicin in Hs578t-VC and Hs578t-TAOK3. Supplement Figure 3. Growth effects of TAOK3 alternation. (A) Growth curve of shTAOK3 MB157 cells (B) Growth curve of TAOK3 overexpression MB157 cells. Supplement Figure 4. TUNEL staining in a similar size subcutaneous xenograft tumor of Hs578t-VC, Hs578t-TAOK3, Hcc1806-NS and Hcc1806-shTAOK3 with paclitaxel (Hs578t: 6 mg/kg and Hcc1806: 3 mg/kg) treatment for 24 h. Supplement Figure 5. Phosphokinase array and microarray analysis of TAOK3 affection. (A) The bot blot image of phosphoprotein array between Hs578t-VC and Hs578t-TAOK3. (B) Bar chart of top 10 increasing phosphorylated proteins. The semi-quantitation was measured with ImageJ.The network of intersection genes based on upstream analysis in (C) TAOK3 overexpression and (D) shRNA knockdown cells. The number showed the fold change of probe from microarray data. Supplement Figure 6. The effects of NF-κB shRNAs in Hs578T with TAOK3 modulation cells. (A) The mitotic percentage changes of NF-B shRNAs and control in Hs578T overexpressed and control cells. (B) The cytotoxicity of paclitaxel of NF-κB shRNAs and control in Hs578T control cells. Supplement Figure 7. IHC staining of TAOK3 in xenograft tumor. Cross-sections of alternative TAOK3 expression xenograft tumor without paclitaxel treatment with TAOK3 IHC staining.
About this article
Cite this article
Lai, TC., Fang, CY., Jan, YH. et al. Kinase shRNA screening reveals that TAOK3 enhances microtubule-targeted drug resistance of breast cancer cells via the NF-κB signaling pathway. Cell Commun Signal 18, 164 (2020). https://doi.org/10.1186/s12964-020-00600-2