Skip to main content
  • Research article
  • Open access
  • Published:

Molecular characterisation of isogenic taxane resistant cell lines identify novel drivers of drug resistance



Taxanes such as paclitaxel and docetaxel are used successfully to treat breast cancer, usually in combination with other agents. They interfere with microtubules causing cell cycle arrest; however, the mechanisms underlying the clinical effects of taxanes are yet to be fully elucidated.


Isogenic paclitaxel resistant (PACR) MDA‒MB‒231, paclitaxel resistant ZR75‒1 and docetaxel resistant (DOCR) ZR75‒1 cell lines were generated by incrementally increasing taxane dose in native cell lines in vitro. We used aCGH analysis to identify mechanisms driving taxane resistance.


Taxane resistant cell lines exhibited an 18-170 fold increased resistance to taxanes, with the ZR75-1 resistant cell lines also demonstrating cross resistance to anthracyclines. Paclitaxel treatment of native cells resulted in a G2/M block and a decrease in the G1 phase of the cell cycle. However, in the resistant cell lines, minimal changes were present. Functional network analysis revealed that the mitotic prometaphase was lost in the resistant cell lines.


This study established a model system for examining taxane resistance and demonstrated that both MDR and mitosis represent common mechanism of taxane resistance.

Peer Review reports


Breast tumours exhibit a wide degree of heterogeneity and diversity at both the cellular and molecular level. The taxanes, paclitaxel and docetaxel, are used successfully to treat breast cancer, alone or in combination with other agents [1]. Taxanes act by interfering with the spindle microtubule dynamics of the cell causing cell cycle arrest followed by cell death [2]. A significant proportion of patients progress despite treatment with taxane containing chemotherapy and there is a pressing need for both novel therapeutic options for patients failing taxane therapy and predictive biomarkers to select patients likely to benefit. Overexpression of P-glycoprotein (PgP/MDR1) is one of the most recognised mechanisms causing taxane resistance [3, 4]. However, several other candidate predictive biomarkers have been proposed in recent studies (AKT/HER2/TLE3) [57], but to date no robust, predictive diagnostic assay for taxane benefit or resistance has emerged. Whilst data suggests some patients are intrinsically resistant to taxanes and others acquire resistance to taxanes as treatment advances there is insufficient understanding of the clinical mechanisms underlying taxane resistance to develop either rational novel therapeutic or diagnostic approaches to target taxane based chemotherapy.

Progress in “targeting” conventional therapeutics such as anthracyclines and taxanes has been slow and has been hampered, in part, by a lack of focus and understanding of the key molecular events that lead to drug response or resistance in the clinical setting. Without significant progress in identifying the key molecular pathways driving drug resistance in vivo, we run the risk of continuing to seek to identify novel drugs and molecular diagnostics in a stochastic and largely unfocused manner.

Genome wide profiling of breast tumours is a powerful tool that can be used to correlate tumour characteristics to clinical outcome in patients. Many extensive studies have proposed novel and molecular subtypes of breast cancer which may have clinical relevance [812]. However few, if any, have proven effective as a basis for either targeting existing treatments or identifying novel therapeutic approaches in the context of drug resistance.

The overall aim of this study was to generate isogenic taxane-resistant breast cancer cell lines and elucidate the mechanisms that are driving resistance to taxanes in a pre-clinical model system. The studies summarised here characterise taxane resistant cell lines derived by the incremental increase of paclitaxel or docetaxel dose. The results presented demonstrate the ZR75-1 resistant cell line harbour cross-resistance to anthracyclines. An aCGH profile demonstrated a loss of mitotic pathways in the resistant cell lines indicating a potential theranostic pathway.


Cell culture and reagents

The breast cancer cell lines MDA-MB-231 and ZR75-1 (ATCC, Cedarlane Laboratories Ltd, Burlington, Canada) were cultured as monolayer in DMEM supplemented with 10% foetal calf serum, 10 mM glutamine and penicillin and streptomycin. Paclitaxel (Sigma, Oakville, Canada), docetaxel (Sigma, Oakville, Canada), epirubicin (Sigma, Oakville, Canada), doxorubicin (Sigma, Oakville, Canada) and carboplatin (Sigma, Oakville, Canada) were dissolved in dimethyl sulphoxide (DMSO) (Sigma, Oakville, Canada). Concentrated stock solutions were stored at -20°C. Drug resistant isogenic daughter cell lines were derived by incremental increases in drug concentrations over time until a stable taxane resistant phenotype was acquired. Cells were in each concentration of drug for two passages and until confluent, this ranged between 1-4weeks dependent on the dose. The following isogenic sub-lines were selected for further characterisation alongside each parent line: MDA-MB-231 25nM and 50nM paclitaxel resistant (MDA-MB-231 25PACR and MDA-MB-231 50PACR), ZR75-1 25nM and 50nM paclitaxel resistant (ZR75-1 25PACR, ZR75-1 50PACR) and 25nM and 50nM docetaxel resistant (ZR75-1 25DOCR, ZR75-1 50DOCR).

IC50 and proliferation rates of parental and isogenic drug resistant lines

Dose response curves were set up by treating cells with increasing doses of the appropriate taxane: 0, 0.3, 1, 3, 10, 30, 100, 300, 1000 or 3000nM of either paclitaxel or docetaxel. Cross resistance to epirubicin, doxorubicin and carboplatin was assessed in a similar manner. Cell suspensions (100μl) were seeded in triplicate at a density of 30,000 cells/ml in 96 well plates and grown for 24 hours, washed and treated with drug for 72 hours. After 72 hours 100μl of growth media containing 10μl of CCK8 (Promega, Madison, USA) was added to each well for 3 hours at 37°C. The plates were then shaken for 10 minutes and optical density (OD) recorded at 450nm. IC50s were calculated using GraphPad Prism 5 (San Diego, USA). Stability of taxane resistance in MDA-MB-231 25PACR was assessed by maintaining the cells for 6 months with or without paclitaxel added to the growth medium. MDA-MB-231 parental cells were maintained without paclitaxel for an equivalent period for comparison.

Flow cytometry

For cell cycle and DNA content analyses, native and resistant cells were plated in equal numbers into 6-well plates and synchronized by serum starvation overnight. Cells were then incubated with the appropriate concentration of taxane (25nM or 50nM of either docetaxel or paclitaxel), DMSO control or media alone control. The cells were collected after 24 and 48 hours, fixed with 80% ethanol and incubated with 2mg/ml RNase A (Sigma, Oakville, Canada) and 0.1mg/ml propidium iodide (Sigma, Oakville, Canada) for 30 minutes prior to analysis by flow cytometry. Data was collected by FACS Canto II and FACS Diva (both from BD Biosciences, Mississauga, Canada), and analyzed by FlowJo (Treesta, San Carlos, USA).

DNA extraction and sample preparation for array Comparative Genomic Hybridisation

DNA was extracted from cells using the Qiagen Blood and Cell Culture Maxi kit (Qiagen, Toronto, Canada). DNA was stored in TE buffer pH 8.0 at 4°C.

Microarray CGH

Cell line DNA was analysed on the Breakthrough Breast Cancer human CGH 4.6K 1.12 arrays as previously described [13]. Briefly, 1 μg of test and normal female genomic DNA, from pooled donor samples, was directly labelled with Cy3-dCTP or Cy5-dCTP (Amersham BioSciences, Amersham, UK) using a Bioprime labelling kit (Invitrogen, Paisley, UK) according to the manufacturer's protocol modified to incorporate 1.0 mM Cy dye, 0.6 mM dCTP, and 1.2 mM dATP, dGTP and dTTP. Unincorporated nucleotides were removed with MinElute purification columns (Qiagen, Crawley, UK). The labelled DNA was co-precipitated with 100 μg of Cot1 (Invitrogen, Paisley, UK), resuspended in hybridization buffer [50% formamide, 10% dextran sulphate, 2× SSC, 2% SDS, 2 mg of yeast tRNA (Invitrogen, Paisley, UK)], denatured at 75°C for 5 min, and pre-annealed for 30 min at 37°C. Slides were blocked in 10% BSA–50% formamide solution at 42°C for 45 min. The probe was subsequently applied to the slide and hybridized overnight at 42°C. Slides were washed in 2× SSC, 0.1% SDS for 15 min at 45°C; 2× SSC, 50% formamide for 15 min at 45°C; 2× SSC, 0.1% SDS for a subsequent 30 min at 45°C; and finally two 15-min washes of 0.2× SSC at room temperature. Slides were centrifuged at 1200 rpm for 2 min to dry. Each experiment was performed in duplicate as a dye swap to eliminate any labelling bias.

Image acquisition and data analysis

Slides were scanned using an Axon 4000B scanner (Axon Instruments, Burlingame, CA, USA) and images were analysed using Genepix Pro 4.1 software (Axon Instruments). The median localized background slide signal for each clone was subtracted and each clone Cy5/Cy3 ratio subjected to print-tip loess normalization [14]. Dye swap experiments were collated, bacterial artificial chromosome (BAC) clone replicate spots averaged, and clones with poor reproducibility between replicates excluded (standard deviation >0.2).

Network-based analysis

To examine whether genes showing common copy number gains or copy number losses across all three cell lines belong to a specific pathway, we conducted functional analysis of the common genes using Cytoscape Reactome Functional Interaction (FI) plugin in Cytoscape 3.0.2 (2013 FI network version). Genes were loaded using the gene set format with FI annotations and linker genes. Spectral clustering was performed to identify distinct network modules and subsequent pathway enrichment was calculated. Symbols were loaded as a gene set and interactions from the FI network 2012 version, including FI annotations and linker genes. Network modules were identified using spectral clustering and Pathway Enrichment computed for each module using the Reactome FI plugin functions. Reactome pathways exhibiting FDR values < 0.01 were considered enriched.

MDR Resistance: RNAi Transfection of ZR75-1 resistant cells

A total of 2.6 × 105 ZR75-1 25PACR cells were transfected with Lipofectamine RNAiMAX (Invitrogen, Paisley, UK) and siRNAs (each 30nM, Dharmacon, Waltman, USA) targeting MDR1, according to manufacturer’s instructions. As controls, transfection reagents without siRNAs were added (mock transfection) and cells were transfected with siRNA targeting GAPDH. After 48h cells were lysed for RNA analysis and 72h cells were lysed for protein analysis. The differences in IC50 were analysed and calculated as described above.

Western blot analysis

Total protein lysates (20 μg) were separated by SDS-PAGE according to standard protocols [15] and immunoblotting was carried out using antibodies directed against PgP-specific MDR1 (G-1) (Santa Cruz Biotechnology, Santa Cruz, CA, USA) diluted 1:1000, GAPDH (14C10) (Cell Signalling, Whitby, Canada) diluted 1:5000 and β-actin (Calbiochem, La Jolla, USA) diluted 1:10000. Horseradish peroxidase–conjugated secondary antibodies were detected by ECL chemiluminescence (Amersham Biosciences, Plc.).


Taxane resistant cell lines IC50s and cross resistance

The taxane resistant cell lines exhibited 18-170 fold increased resistance to taxanes, when IC50s were compared to those from parental cell lines, with cross resistance to both forms of taxane observed in all cell lines (Table 1). All ZR75-1 PACR and DOCR cell lines exhibited cross resistance to anthracyclines (epirubicin and doxorubicin); however, no cross-resistance was observed with carboplatin. MDA-MB-231 PACR cells were not cross-resistant to either anthracyclines or carboplatin (Table 1).Following long term (6 months) culture of MDA-MB-231 25PACR cells in the absence of drug, cells were re-challenged with taxanes and the responses compared to parental and resistant cells cultured in the presence of taxanes (Figure 1). A two way Anova analysis of the proliferation data between the native and resistant cells with or without paclitaxel was performed in a pairwise fashion. When the two resistant cell lines were compared there was no significant difference between the two lines (p = 0.09728), indicating that they exhibited a very similar paclitaxel resistant phenotype.

Table 1 IC 50 values (nM) for paclitaxel, docetaxel, epirubicin, doxorubicin and carboplatin in isogenic MDA-MB-231 and ZR75-1 cell lines
Figure 1
figure 1

MDA-MB-231 25PACR cells maintained resistance to paclitaxel after prolonged culture without exposure to taxane. MDA-MB-231 25PACR cells were separated into two groups: one maintained and passaged, as normal in the presence of paclitaxel (white bar), the other was maintained and passaged in the absence of drug (grey bar)for a period of six months. The native cells are represented by the black bar. Cells were incubated with varying concentrations of paclitaxel and cell viability determined by CCK-8 assay. The X axis shows the increasing paclitaxel concentration measured in nM. The Y axis represents the percentage of cells with untreated cells being used as a baseline of 100%.

Cell cycle specific effects of taxanes

Paclitaxel treatment of native MDA-MB-231 and ZR75-1 cells resulted in a G2/M block, and a failure to return to the G0/G1 phase (Figure 2). The G2/M population of the MDA-MB-231 native cells increased significantly from 24% to 44% upon paclitaxel exposure compared with a minimal change of 24% to 19% in the MDA-MB-231 25PACR cells. The increase of cell population at the G2/M phase was accompanied by a decrease of cell population in the G1 phase of the cell cycle for the native cells; however the resistance cell lines exhibited no change in the percentage of cells in the G1 phase. Paclitaxel treatment of native ZR75-1 cells resulted in a significant increase in the G2/M population from 18% to 72% and a decrease in the G1 population from 48% to 11%. While in the ZR75-1 25PACR cells there were minimal changes in the G2/M population from 20% to 30%, there was a decrease in the G1 population of cells from 59% to 28%. Treatment of the ZR75-1 50PACR cells with paclitaxel caused a slight decrease in the G2/M population of cells from 14% to 8% and a slight change in the G1 population of cells from 61% to 71% (Figure 2).

Figure 2
figure 2

Cell cycle analysis of native and resistant MDA-MB-231 and ZR75-1 by flow cytometry after synchronisation. The native and respective resistant cell lines were treated with 25nM or 50nM of paclitaxel. The DNA content was measured by flow cytometry to determine the distribution of cell in each phase. The histograms demonstrate the cell cycle distribution within the cell population. A MDA-MB-231 native and MDA-MB231 25PACR cells with or without 25nM paclitaxel. B MDA-MB-231 native and MDA-MB231 25PACR cells with or without 50nM paclitaxel. C ZR75-1 native and ZR75-1 25PACR cells with or without 25nM paclitaxel. D ZR75-1 native and ZR75-1 25PACR cells with or without 50nM paclitaxel. E ZR75-1 native and ZR75-1 25DOCR cells with or without 25nM docetaxel. F ZR75-1 native and ZR75-1 50DOCR cells with or without 50nM docetaxel. Standard deviation of three experiments are shown in brackets.

Array comparative genomic Hybridisation


The PACR cell lines were analysed and compared to the parental controls (Figure 3A and B). Both the 25PACR and 50PACR cells demonstrated marked gains and losses (Table 2 and Figure 3A and B). There are three common areas of genomic loss in the MDA-MB-231 cell lines that extend with increasing paclitaxel resistance in chromosome 1p, 6p and 17p. Common areas of gain include 8q and 15p.

Figure 3
figure 3

aCGH of taxane resistant cell lines. The plots show Log2Ratios of test to reference signal intensity from BAC clines in an aCGH experiment using DNA from native cells as a reference samples and DNA from resistant cells as a test sample. Navy dots represent BAC clones which remain unchanged, the green dots represent the BAC clones in which there is an area of gain on the genome, and the red dots represent the BAC clones in which there is an area of loss of the genome. The Log2ratio is measured on the Y axis and on the X axis the genome runs in chromosome order from 1 to the sex chromosomes. The p or short arm on each chromosome is followed by the q or long arm. The dotted lines represent the position of the centromere. The cbs algorithm recursively split chromosomes into segments based on the maximum t statistic estimated by each permutation (re A. MDA-MB-231 Natives vs MDA-MB-231 25PACR. B. MDA-MB-231 Natives vs MDA-MB-231 50PACR. C. ZR75-1 native cells vs ZR75-1 25PACR. D. ZR75-1 native cells vs ZR75-1 50PACR. E. ZR75-1 native cells vs ZR75-1 25DOCR. F. ZR75-1 native cells vs ZR75-1 50DOCR.

Table 2 Common areas of loss, gain, deletion and amplification identified by aCGH in MDA-MB-231 PACR, ZR75-1 PACR and ZR75-1 DOCR at the two resistance levels 25nM and 50nM when compared to the native cell line


In the ZR75-1 cell lines there were fewer genomic changes that occurred once cell becomes resistant in contrast to the MDA-MB-231 (Table 2, Figure 3C-F). There were common areas of gains and losses in the 25PACR and 50PACR cells; losses were observed in 3p, 7q, 10p, 12p and 15p. Interestingly within region 7p22.3-q11.21 the gene ABCB5, a member of the p-glycoprotein family, is present and appears to be gained. There were no common areas of gain in the ZR75-1 PACR cell lines.

25DOCR and 50DOCR cells compared with native cells show area of loss in 7q, 12p and 16q again there were no common areas of gain with the DOCR cell lines.

When comparing the data obtained from the PACR and DOCR ZR75-1 cells the sole areas of common genomic alterations were losses at 7q and 12p.

Combined analysis

When all areas of gain or loss across the 25nM resistant cell lines were combined, 295 known genes were identified as lost and 306 genes gained (Figure 4A and 4B). Following network analysis, eight modules were identified that contained significantly enriched pathways with a False Discovery Rate (FDR) <0.01. Each module contained clusters of connected genes. Module II contained 6 genes involved in the mitotic prometaphase. Interestingly, all six genes were deleted in taxane resistance cells and directly interconnected without linker genes (Figure 4C). These findings would suggest that loss of mitotic prometaphase regulatory genes is a common event associated with taxane resistance in breast cancer cells.

Figure 4
figure 4

Network-based analysis of MDA-MB-231 25PACR, ZR75-1 25PACR and ZR75-1 25DOCR taxane resistant cell lines. A. Venn diagram of genes within significant areas of gain. B. Venn diagram of genes within significant areas of loss in 3 cell lines. C. Mitotic prometaphase module identified from functional interaction network analysis. D. Bar graph shows relative expression levels of AHCTH1, CENPF, PPP2R5A, NUP133, MLP1P and NSL1 in the MDA-MB-231 native, 25PACR, 50 PACR ZR75-1 25PACR, 50PACR and ZR75-1 25DOCR and 50DOCR taxane resistant cell lines. Error bars show standard deviation of three experiments. Asterisks indicate statistical difference * = p < 0.05, ** = p < 0.001.

qRT-PCR validation of aCGH

qRT-PCR analysis was performed on the six deleted genes present in the resistant cell lines compared to parental controls. As shown in Figure 4C in the MDA-MB-231 resistant cell lines all six of the genes were downregulated compared to the parental control cells. Within the ZR75 cell lines downregulation of all resistant cell lines was demonstrated with AHCTH1 and NUP133. MLP1IP showed a decrease in expression in the both DOCR and 50PACR cells compared to the natives while the 25PACR cells showed an increase in expression.

MDR1 is a driver of taxane resistance in ZR75-1 cells only

No MDR‒1 protein expression was identified by western blotting in the MDA‒MB‒231 native, MDA-MB-231 25PACR or MDA-MB-231 50PACR cell lines (Figure 5A). There was a large increase in MDR1 protein expression in all four taxane resistant ZR75‒1 cell lines while no expression of the protein was observed in the ZR75‒1 native line. Western blot and cell proliferation assays were performed after down-regulation of MDR1 using siRNA. Western blot analysis demonstrated a reduction in MDR1 expression following transfection with siRNA (Figure 5B). In the proliferation assay MDR1 knock-down exhibited a 14- and 34-fold reduction in the IC50 concentration of paclitaxel in both the ZR75-1 25PACR and ZR75-1 50PACR cells respectively (Figure 5C and D). This corroborates previous western blot analysis suggesting MDR1 is the driver of taxane resistance in the ZR75-1 cell lines. There were no differences in α/β tubulin expression in either the ZR75 or MDA-MB-231 resistant cell lines compared to the native parental lines (Additional file 1: Figure S1). Taken together this would suggest taxane resistance in these cell lines is at least partially driven by MDR1 expression

Figure 5
figure 5

Effects of MDR-1 knockdown in ZR75-1 paclitaxel resistant cells A Western blot analysis of proteins extracted from the cell lines and probed with MDR1.Actin was used as a loading control. B. Western blot analysis of proteins extracted from ZR75-1 25PACR cells after transfections with MDR1 siRNA. Control cells were untreated, or with transfections reagents only or with siRNA targeting GAPDH; actin was used as a loading control. C. Graph shows the average IC50 of ZR75-1 25PACR cells after transfections with MDR1 siRNA. D. Graph shows the average IC50 of ZR75-1 50PACR cells after transfections with MDR1 siRNA. Asterisks indicate statistical difference * = p < 0.05, ** = p < 0.001.


The taxanes are a useful and effective group of chemotherapeutic agents that can be used as front line therapy to treat many types of cancer including breast, ovarian and prostate. Unfortunately, taxane resistance is a considerable clinical problem, and overcoming this is a key step to improving breast cancer patient survival. One way of combating this is to identify potential molecular drivers of taxane resistance so that they can be targeted with combination therapies to down-regulate the resistant phenotype.

In this study, a panel of isogenic paclitaxel resistant cell lines were generated by exposing parental cells to increasing concentrations of the appropriate taxane in vitro. We successfully generated daughter cell lines with markedly increased IC50s for taxanes; demonstrating clear resistance to these agents. Our cell cycle analysis demonstrated, in the native/parental cells, treatment with either docetaxel or paclitaxel resulted in a G2/M block. However, the drug resistant cell lines were able to overcome this G2/M block and progress through the cell cycle. Interestingly, we demonstrated that both the MDA‒MB‒231 and the ZR75‒1 native cell lines were more sensitive to docetaxel than paclitaxel which concurs with previous studies using other cell lines [16]. Other studies have suggested that docetaxel and paclitaxel affect different stages of the cell cycle with paclitaxel only targeting G2/M whilst docetaxel targets both S phase and G2/M [17]. In our hands treatment with both docetaxel and paclitaxel resulted in a G2/M block. Once reason that we may not see an S phase block is due to the concentration of drug that we were using. One study by Hernández‒Vargas used synchronized cells and then subjected them to low (2‒4nM) or high (100nM) concentrations of docetaxel. The low dose treatment caused a transient arrest and the high dose cause a prolonged arrest in mitosis. The short arrest leads to an aberrant mitosis and aneuploidy whereas the long arrest leads to mitotic slippage and tetraploidy [18]. A dual mechanism of cell cycle response has also been seen with paclitaxel treatment [19]. Low doses of paclitaxel have been shown to inhibit or retard the progression of mitosis and as a consequence alter microtubule dynamics rather than actually increasing polymer mass [20, 21]. At higher concentrations of paclitaxel cells become blocked in G2/M phase so that they cannot progress through mitosis.

Both ZR75-1 resistant cell lines showed cross-resistance to anthracyclines. This is consistent with a clinical trial of first-line treatment with anthracyclines followed by a crossover to taxanes which showed reduced response to taxanes [17], suggesting that anthracycline treatment may induce taxane cross-resistance. At a protein level both resistant ZR75-1 cell lines exhibited up-regulation of MDR1 suggesting that the resistance observed in these cells is mediated by the MDR family. The MDR family of p-glycoproteins are a common resistance mechanism observed in numerous in vitro studies. These proteins bind non-specifically to multiple chemotherapy drugs and actively export them across the cellular membrane [22, 23]. However, the clinical relevance of MDR genes remains to be elucidated. No cross-resistance was detected in the MDA-MB-231 PACR cells, consistent with previous studies performed in MCF-7 paclitaxel resistant cell lines [24] and suggests that a mechanism other than PgP glycoprotein is driving this resistance.

We performed a genomic analysis of our cell lines using aCGH. We compared each of the cell lines to DNA from pooled female blood and then compared each of the resistant lines with their respective native lines. Within the MDA-MB-231 cell, areas on chromosome 1p, 6p and 17p were lost and gains in chromosome 8q and 15p were observed. When comparing the native MDA-MB-231 cell line with the 50PACR cell lines, amplification at chromosome 1q was observed. Chromosome 1 aberrations are the most frequently described in a variety of cancers [25]. In breast cancer 1q gain is commonly observed across all subtypes, however the functional driver in this region has yet to be elucidated. There are many candidate genes; CENF, KIF14, DTL, NEK2, CKS1B, ASPM and EXO1 each of which are significantly associated with poor clinical outcome in breast cancer patients [26]. Interestingly, when functional network analysis was performed incorporating all three paclitaxel resistant cell lines, a signalling module which included genes controlling the mitotic prometaphase was identified. Five of the genes, PP2R5A, NUP133, AHCTF1, CENPF and NSL1, within this module are located on chromosome 1. Previously studies have demonstrated CENPF as both a prognostic and predictive gene in breast cancer [27]. One study showed CENPF to be associated with poor prognosis [28]. Paclitaxel enhances the stability of microtubules and mitosis is blocked at the metaphase-anaphase transition with prolonged blocking resulting in cell death. However, taxane resistant cells drug appear to have lost the ability to control this process. Drug resistant cells, when treated with taxane, progress through the cell cycle without arresting in G2/M suggesting they are bypassing a critical cell cycle checkpoint. Dysregulation of the mitotic metaphase check point is linked to chromosome instability (CIN). CIN cells become aneuploid and are associated with aggressive tumours and poor prognosis. CIN has been previously linked to taxane resistance in ovarian and colorectal cancer [29, 30]. Therefore, it would suggest that once these in vitro cell lines become taxane resistant they also become genomically unstable and therefore may be at greater risk of progression.


In conclusion, our study has established a new model system to examine mechanisms of taxane resistance in breast cancer with genomic analysis showing a mitotic prometaphase as a predictor of resistance.


  1. Effects of chemotherapy and hormonal therapy for early breast cancer on recurrence and 15-year survival: an overview of the randomised trials. Lancet. 2005, 365: 1687-1717.

  2. 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. 10.1016/j.ctrv.2012.02.011.

    Article  CAS  PubMed  Google Scholar 

  3. Brooks TA, Minderman H, O'Loughlin KL, Pera P, Ojima I, Baer MR, Bernacki RJ: Taxane-based reversal agents modulate drug resistance mediated by P-glycoprotein, multidrug resistance protein, and breast cancer resistance protein. Mol Cancer Ther. 2003, 2 (11): 1195-1205.

    CAS  PubMed  Google Scholar 

  4. O'Brien C, Cavet G, Pandita A, Hu X, Haydu L, Mohan S, Toy K, Rivers CS, Modrusan Z, Amler LC, Lackner MR: Functional genomics identifies ABCC3 as a mediator of taxane resistance in HER2-amplified breast cancer. Cancer Res. 2008, 68 (13): 5380-5389. 10.1158/0008-5472.CAN-08-0234.

    Article  PubMed  Google Scholar 

  5. Hayes DF, Thor AD, Dressler LG, Weaver D, Edgerton S, Cowan D, Broadwater G, Goldstein LJ, Martino S, Ingle JN, Henderson IC, Norton L, Winer EP, Hudis CA, Ellis MJ, Berry DA: HER2 and response to paclitaxel in node-positive breast cance. Adv Physiol Educ. 2007, 357: 1496-1506.

    CAS  Google Scholar 

  6. Bartlett JM, A'Hern R, Piper T, Ellis IO, Dowsett M, Mallon EA, Cameron DA, Johnston S, Bliss JM, Ellis P, Barrett-Lee PJ: Phosphorylation of AKT pathway proteins is not predictive of benefit of taxane therapy in early breast cancer. Breast Cancer Res Treat. 2013, 138 (3): 773-781. 10.1007/s10549-013-2489-y.

    Article  CAS  PubMed  Google Scholar 

  7. Kulkarni SA, Hicks DG, Watroba NL, Murekeyisoni C, Hwang H, Khoury T, Beck RA, Ring BZ, Estopinal NC, Schreeder MT, Seitz RS, Ross DT: TLE3 as a candidate biomarker of response to taxane therapy. Breast Cancer Res. 2009, 11 (2): R17-10.1186/bcr2241.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Curtis C, Shah SP, Chin SF, Turashvili G, Rueda OM, Dunning MJ, Speed D, Lynch AG, Samarajiwa S, Yuan Y, Graf S, Ha G, Haffari G, Bashashati A, Russell R, McKinney S, Langerod A, Green A, Provenzano E, Wishart G, Pinder S, Watson P, Markowetz F, Murphy L, Ellis I, Purushotham A, Borresen-Dale AL, Brenton JD, et al: The genomic and transcriptomic architecture of 2,000 breast tumours reveals novel subgroups. Nature. 2012, 486 (7403): 346-352.

    CAS  PubMed  PubMed Central  Google Scholar 

  9. Shah SP, Roth A, Goya R, Oloumi A, Ha G, Zhao Y, Turashvili G, Ding J, Tse K, Haffari G, Bashashati A, Prentice LM, Khattra J, Burleigh A, Yap D, Bernard V, McPherson A, Shumansky K, Crisan A, Giuliany R, Heravi-Moussavi A, Rosner J, Lai D, Birol I, Varhol R, Tam A, Dhalla N, Zeng T, Ma K, et al: The clonal and mutational evolution spectrum of primary triple-negative breast cancers. Nature. 2012, 486 (7403): 395-399.

    CAS  PubMed  Google Scholar 

  10. Perou CM, Sorlie T, Eisen MB, van de Rijn M, Jeffrey SS, Rees CA, Pollack JR, Ross DT, Johnsen H, Akslen LA, Fluge O, Pergamenschikov A, Williams C, Zhu SX, Lonning PE, Borresen-Dale AL, Brown PO, Botstein D: Molecular portraits of human breast tumours. Adv Physiol Educ. 2000, 406: 747-752.

    CAS  Google Scholar 

  11. Rouzier R, Perou CM, Symmans WF, Ibrahim N, Cristofanilli M, Anderson K, Hess KR, Stec J, Ayers M, Wagner P, Morandi P, Fan C, Rabiul I, Ross JS, Hortobagyi GN, Pusztai L: Breast cancer molecular subtypes respond differently to preoperative chemotherapy. Adv Physiol Educ. 2005, 11: 5678-5685.

    CAS  Google Scholar 

  12. Sorlie T, Perou CM, Tibshirani R, Aas T, Geisler S, Johnsen H, Hastie T, Eisen MB, van de Rijn M, Jeffrey SS, Thorsen T, Quist H, Matese JC, Brown PO, Botstein D, Lonning PE, Borresen-Dale AL: Gene expression patterns of breast carcinomas distinguish tumor subclasses with clinical implications. Proc Natl Acad Sci U S A. 2001, 98 (19): 10869-10874. 10.1073/pnas.191367098.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Natrajan R, Little SE, Sodha N, Reis-Filho JS, Mackay A, Fenwick K, Ashworth A, Perlman EJ, Dome JS, Grundy PE, Pritchard-Jones K, Jones C: Analysis by array CGH of genomic changes associated with the progression or relapse of Wilms' tumour. J Pathol. 2007, 211 (1): 52-59. 10.1002/path.2087.

    Article  CAS  PubMed  Google Scholar 

  14. Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J, Speed TP: Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Res. 2002, 30 (4): e15-10.1093/nar/30.4.e15.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Spears M, Cunningham CA, Taylor KJ, Mallon EA, Thomas JS, Kerr GR, Jack WJ, Kunkler IH, Cameron DA, Chetty U, Bartlett JM: Proximity ligation assays for isoform-specific Akt activation in breast cancer identify activated Akt1 as a driver of progression. J Pathol. 2012, 227 (4): 481-489. 10.1002/path.4022.

    Article  CAS  PubMed  Google Scholar 

  16. Hennequin C, Giocanti N, Favaudon V: S-phase specificity of cell killing by docetaxel (Taxotere) in synchronised HeLa cells. Br J Cancer. 1995, 71 (6): 1194-1198. 10.1038/bjc.1995.232.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Gligorov J, Lotz JP: Preclinical pharmacology of the taxanes: implications of the differences. Oncologist. 2004, 9 (Suppl 2): 3-8.

    Article  CAS  PubMed  Google Scholar 

  18. Hernandez-Vargas H, Palacios J, Moreno-Bueno G: Telling cells how to die: docetaxel therapy in cancer cell lines. Cell Cycle. 2007, 6 (7): 780-783. 10.4161/cc.6.7.4050.

    Article  CAS  PubMed  Google Scholar 

  19. Torres K, Horwitz SB: Mechanisms of Taxol-induced cell death are concentration dependent. Cancer Res. 1998, 58 (16): 3620-3626.

    CAS  PubMed  Google Scholar 

  20. Jordan MA, Toso RJ, Thrower D, Wilson L: Mechanism of mitotic block and inhibition of cell proliferation by taxol at low concentrations. Proc Natl Acad Sci U S A. 1993, 90 (20): 9552-9556. 10.1073/pnas.90.20.9552.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Jordan MA, Wendell K, Gardiner S, Derry WB, Copp H, Wilson L: Mitotic block induced in HeLa cells by low concentrations of paclitaxel (Taxol) results in abnormal mitotic exit and apoptotic cell death. Cancer Res. 1996, 56 (4): 816-825.

    CAS  PubMed  Google Scholar 

  22. Schinkel AH, Roelofs EM, Borst P: Characterization of the human MDR3 P-glycoprotein and its recognition by P-glycoprotein-specific monoclonal antibodies. Cancer Res. 1991, 51 (10): 2628-2635.

    CAS  PubMed  Google Scholar 

  23. Van der Bliek AM, Baas F, Van der V, Biedler JL, Meyers MB, Ozols RF, Hamilton TC, Joenje H, Borst P: Genes amplified and overexpressed in human multidrug-resistant cell lines. Cancer Res. 1988, 48 (21): 5927-5932.

    CAS  PubMed  Google Scholar 

  24. Villeneuve DJ, Hembruff SL, Veitch Z, Cecchetto M, Dew WA, Parissenti AM: cDNA microarray analysis of isogenic paclitaxel- and doxorubicin-resistant breast tumor cell lines reveals distinct drug-specific genetic signatures of resistance. Breast Cancer Res Treat. 2006, 96: 17-39. 10.1007/s10549-005-9026-6.

    Article  CAS  PubMed  Google Scholar 

  25. Kytola S, Rummukainen J, Nordgren A, Karhu R, Farnebo F, Isola J, Larsson C: Chromosomal alterations in 15 breast cancer cell lines by comparative genomic hybridization and spectral karyotyping. Genes Chromosom Cancer. 2000, 28 (3): 308-317. 10.1002/1098-2264(200007)28:3<308::AID-GCC9>3.0.CO;2-B.

    Article  CAS  PubMed  Google Scholar 

  26. Muthuswami M, Ramesh V, Banerjee S, Viveka Thangaraj S, Periasamy J, Bhaskar Rao D, Barnabas GD, Raghavan S, Ganesan K: Breast tumors with elevated expression of 1q candidate genes confer poor clinical outcome and sensitivity to Ras/PI3K inhibition. PLoS One. 2013, 8 (10): e77553-10.1371/journal.pone.0077553.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. van 't Veer LJ, Dai H, van de Vijver MJ, He YD, Hart AA, Mao M, Peterse HL, van der Kooy K, Marton MJ, Witteveen AT, Schreiber GJ, Kerkhoven RM, Roberts C, Linsley PS, Bernards R, Friend SH: Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002, 415 (6871): 530-536. 10.1038/415530a.

    Article  PubMed  Google Scholar 

  28. O'Brien SL, Fagan A, Fox EJ, Millikan RC, Culhane AC, Brennan DJ, McCann AH, Hegarty S, Moyna S, Duffy MJ, Higgins DG, Jirstrom K, Landberg G, Gallagher WM: CENP-F expression is associated with poor prognosis and chromosomal instability in patients with primary breast cancer. Int J Cancer. 2007, 120 (7): 1434-1443. 10.1002/ijc.22413.

    Article  PubMed  Google Scholar 

  29. Birkbak NJ, Eklund AC, Li Q, McClelland SE, Endesfelder D, Tan P, Tan IB, Richardson AL, Szallasi Z, Swanton C: Paradoxical relationship between chromosomal instability and survival outcome in cancer. Cancer Res. 2011, 71 (10): 3447-3452. 10.1158/0008-5472.CAN-10-3667.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Swanton C, Nicke B, Schuett M, Eklund AC, Ng C, Li Q, Hardcastle T, Lee A, Roy R, East P, Kschischo M, Endesfelder D, Wylie P, Kim SN, Chen JG, Howell M, Ried T, Habermann JK, Auer G, Brenton JD, Szallasi Z, Downward J: Chromosomal instability determines taxane response. Proc Natl Acad Sci U S A. 2009, 106 (21): 8671-8676. 10.1073/pnas.0811835106.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Pre-publication history

Download references


JK was funded by the UK Medical Research Council. ML, AM and JRF were funded through the Institute of Cancer Research. KJT and CC were funded through University of Edinburgh. MS, LL, CY and JMSB were supported by the funding from OICR. We thank the government of Ontario for funding, which is provided through the Ontario Ministry of Research and Innovation.

Author information

Authors and Affiliations


Corresponding author

Correspondence to John MS Bartlett.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

JK: design of the study and performed some of the experiments; MS experimental design, data analysis and drafted the manuscript, NL performed experimental procedures and data analysis, KJT performed experimental procedures and data analysis, LL performed experimental procedures and data analysis, CC performed experimental procedures, ML participated in aCGH experiments, AM data analysis, CY data analysis JRF participated in experimental design. JMSB experimental design, coordination and writing of the manuscript. All authors read and approved the final manuscript.

Juliet Kenicer, Melanie Spears contributed equally to this work.

Electronic supplementary material


Additional file 1: Figure S1: Western blot analysis of proteins extracted from the cell lines and probed with α/β Tubulin. GAPDH was used as a loading control. (DOCX 173 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit

The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kenicer, J., Spears, M., Lyttle, N. et al. Molecular characterisation of isogenic taxane resistant cell lines identify novel drivers of drug resistance. BMC Cancer 14, 762 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: