- Research article
- Open Access
Cisplatin-resistant triple-negative breast cancer subtypes: multiple mechanisms of resistance
BMC Cancer volume 19, Article number: 1039 (2019)
Understanding mechanisms underlying specific chemotherapeutic responses in subtypes of cancer may improve identification of treatment strategies most likely to benefit particular patients. For example, triple-negative breast cancer (TNBC) patients have variable response to the chemotherapeutic agent cisplatin. Understanding the basis of treatment response in cancer subtypes will lead to more informed decisions about selection of treatment strategies.
In this study we used an integrative functional genomics approach to investigate the molecular mechanisms underlying known cisplatin-response differences among subtypes of TNBC. To identify changes in gene expression that could explain mechanisms of resistance, we examined 102 evolutionarily conserved cisplatin-associated genes, evaluating their differential expression in the cisplatin-sensitive, basal-like 1 (BL1) and basal-like 2 (BL2) subtypes, and the two cisplatin-resistant, luminal androgen receptor (LAR) and mesenchymal (M) subtypes of TNBC.
We found 20 genes that were differentially expressed in at least one subtype. Fifteen of the 20 genes are associated with cell death and are distributed among all TNBC subtypes. The less cisplatin-responsive LAR and M TNBC subtypes show different regulation of 13 genes compared to the more sensitive BL1 and BL2 subtypes. These 13 genes identify a variety of cisplatin-resistance mechanisms including increased transport and detoxification of cisplatin, and mis-regulation of the epithelial to mesenchymal transition.
We identified gene signatures in resistant TNBC subtypes indicative of mechanisms of cisplatin. Our results indicate that response to cisplatin in TNBC has a complex foundation based on impact of treatment on distinct cellular pathways. We find that examination of expression data in the context of heterogeneous data such as drug-gene interactions leads to a better understanding of mechanisms at work in cancer therapy response.
A major goal of improved classification of cancer subtypes is to stratify patient populations and to more rapidly identify effective treatment strategies. Advances in molecular characterization of tumors not only improve classification, but also point directly to molecular mechanisms that lead to different therapeutic responses. By integrating heterogenous functional genomic data on tumor subtype characteristics, with known mechanisms and pathways and molecular response to drugs, it is possible to match drug response to tumor characteristics, thus refining treatment options.
Subtypes of TNBC
Classification of cancer subtypes relies on many criteria including histological typing, mutation status, genomic structural variations and expression profiling [1,2,3,4,5]. Breast cancers are often classified by the presence or absence of three receptors: estrogen receptor (ESR1), progesterone receptor (PGR), and the HER2 epidermal growth factor receptor (ERBB2) [6, 7]. Tumors that lack expression of all three receptors are called triple-negative breast cancer (TNBC). As many available therapies in breast cancer target one of these receptors, TNBC status limits treatment options. TNBC is particularly aggressive with higher rates of recurrence, metastasis, and mortality than other breast cancers [8, 9].
Additionally, breast cancers are typically classified as luminal, basal/myoepithelial or ERBB2- subtypes based on relation to cell types found in the normal breast . Although most TNBC cancers are characterized as basal-like, about 20% of TNBC tumors are classified as non-basal .
Two recent studies have classified TNBCs based on clustering genes that are up and down-regulated resulting in six and four molecularly defined subtypes, respectively [4, 5]. Lehmann et al. initially described and tested chemotherapy response in six TNBC subtypes: basal-like 1 (BL1), basal-like 2 (BL2), immunomodulatory (IM), mesenchymal (M), mesenchymal stem-like (MSL) and luminal androgen receptor (LAR) . In another study, Burstein et al. also used gene-expression profiling to subclassify TNBC into four subtypes: mesenchymal (MES), luminal AR (LAR), basal-like immune suppressed (BLIS) and basal-like immune activated (BLIA) . Burstein et al. compared their classifications with the Lehmann classifications and showed that there was some concordance with the LAR/LAR, MSL/MES and M/BLIS type tumors from both groups, but little discrimination of the BL1, BL2 and IM subtypes . For our analysis, we used sets from four of the subtypes described by Lehmann et al: BL1, BL2, M, and LAR  (more details below).
Treatment of TNBC
There are no targeted treatments for TNBC . Standard treatment for TNBC patients includes chemotherapy and surgery and patients often become refractory to the treatment [14, 15]. Patients that achieve a complete response during neoadjuvant therapy generally have better outcomes . Recent strategies for the treatment of TNBC define different treatments depending on BRCA gene status and CD274 (PD-L1) expression status . Treatments addressed include chemotherapy, immunotherapy, and PARP inhibitor therapy. First-line chemotherapeutic agents include taxane and anthracycline, which can be used singly or in combination, but these agents can be augmented with other treatments in cases of progression or contraindications .
TNBC and Cisplatin
Although not currently considered standard of care for TNBC, there is renewed interest in cisplatin use to treat TNBC . Cisplatin has been in use for over 40 years to treat multiple types of cancer. Substatial data correlating cisplatin sensitivity with respect to TNBC subtypes and curated data associating cisplatin with interacting genes provides a robust data collection for integrated analysis. This provides a unique opportunity to study the genetic mechanisms that underlie TNBC subtypes and their relation to cisplatin.
Currently, 22 clinical trials are exploring the use of cisplatin to treat TNBC either as a single agent or in combination with other therapies  (Search criteria were: not yet recruiting, recruiting, enrolling by invitation, and active, not recruiting accessed 01/22/2019). In particular, use of cisplatin therapy has been suggested for TNBC harboring a BRCA mutation . Cisplatin is a DNA-intercalating agent that cross-links DNA resulting in interference with RNA transcription and DNA replication activities. If the DNA lesions are not repaired, DNA-damage induced cell-cycle arrest and apoptosis are triggered [20, 21]. Cells can become resistant to cisplatin by several mechanisms including change in the accumulation of the drug in cells either by inhibited uptake or enhanced efflux, detoxification of the drug by redox mechanisms, repair of the DNA by excision repair mechanisms, or negative regulation of apoptotic mechanisms [22,23,24,25].
New insights into the biological processes associated with cisplatin in different molecular subtypes of TNBC may lead to  a better understanding of the mechanisms underlying treatment response differences,  strategies for identifying those patients that are more likely to respond robustly to chemotherapy, and  the identification of new treatment strategies.
Our approach is to integrate and analyze curated information from pathways and mechanisms obtained in multiple species with empirical data collected in tumor profiling and mechanistic experiments. This allows us to focus, in a ‘sea’ of differentially expressed genes, on genes related to specific areas of interest--in our case genes related to the biology of cisplatin. In this study, we used the GeneWeaver (GW) gene set analysis platform  to identify specific biological processes that could explain the observation that of the TNBC subtypes, BL1 and BL2 are more sensitive to cisplatin than M and LAR . We focus on these four subtypes because the MSL and IM subtypes were later shown to contain stromal cells and infiltrating lymphocytes respectively . GW comprises a database of gene sets from multiple functional genomics data resources, curated publications and user submisisons. These data resources are provided with a suite of combinatorial and statistical tools for performing set operations on user selected gene lists. This provided a platform for the comparison of genomic profiles of multiple TNBC subtypes and gene products with a chemotherapeutic drug. To create the gene sets for our study we first identified evolutionarily conserved genes that were associated with cellular or physiological responses to cisplatin. We then identified which of the genes in the conserved cisplatin-associated set were found among genes shown previously to be differentially expressed in TNBC molecular subtypes. Finally, we analyzed the differentially expressed, cisplatin-associated genes with respect to biological processes and to pathways associated with sensitivity or resistance to cisplatin (Fig. 1).
To investigate these genes in the context of TNBC, we expanded the gene-set collection in GW by constructing genes sets for the differentially regulated genes described by Lehmann et al., , thereby making gene sets for identified up- and down-regulated genes for each of the six molecular subtypes of TNBC. For our analysis, we used sets from four of the subtypes that were subsequently shown not to contain infiltrating cells: BL1, BL2, M, and LAR .
For all gene sets, we used Human Genome Nomenclature Committee (HGNC)-approved identifiers. Genes that we could not unambiguously assign to an HGNC identifier were not included. Details of the source and methods of curation are reported for each of the gene set descriptions as part of the GW record. For ontology-tagging, TNBC-gene sets were annotated with the Disease Ontology term ‘triple-receptor negative breast cancer’ (DOID:0060081), and the Human Phenotype Ontology term ‘Breast carcinoma’ (HP:0003002) ([27, 28], respectively). Gene sets with known response to cisplatin were tagged with the Chemicals of Biological Interest (ChEBI) term ‘cisplatin’ (CHEBI:27899) .
To create a set of human genes associated with cisplatin that are evolutionarily conserved, we identified gene sets associated with studies of cisplatin in GW’s database and applied combinatorial tools to selected sets as outlined below  (Fig. 2).
Using existing gene sets in GW we identified 34 cisplatin-associated gene sets that included sets obtained from GWAS studies (22 sets), MESH terms (2 sets) and the Comparative Toxicogenomics Database (CTD) (10 sets) respectively. CTD curates many aspects of gene-chemical interactions including regulatory, physical interaction, responses, and interactions that are reported as a result of interactions of cisplatin combined with other treatments . The provenance of chemical-gene associations is fully traceable back to the original source. For eample the association of the gene RAD51 with cisplatin can be traced back to three separate publications and three different species using the CTD resource (Query performed on Sept. 3, 2019).
We selected three large data sets from CTD for further analysis, one each from human, mouse and rat. The selected sets consisted of 2386 (GS125959), 883 (GS257116) and 616 (GS263765) genes from human, mouse and rat respectively. We chose these sets as ‘high-confidence’ sets because CTD data includes a large corpus of gene-chemical associations curated from published literature .
To identify genes associated with biological processes that are also evolutionarily conserved, and that therefore could be considered central to the action of cisplatin, we identified orthologous genes that share an association with cisplatin in CTD.
To examine the orthologous gene overlap of these species-specific sets, we used the GW Hierarchical Similarity (HiSim) Graph tool . This tool creates a graph in which leaves represent individual gene sets in the selection, and parent nodes represent sets of genes in the intersection of all child nodes (analysis date 9/2/19). Gene-overlap between mouse-human, rat-human and mouse-rat sets were 378, 219 and 151 genes respectively. We used the genes in the intersection of all three cisplatin-response sets to generate a new gene set of the 96 human cisplatin-associated genes whose homologs are conserved among the three species (GS271882) (Fig. 3). To supplement the data from the human CTD gene set, we performed the same analysis with an additional publically available gene set in GW, GS237976: [MeSH] Cisplatin:D002945. This analysis resulted in the identification of six more conserved genes: GJA1, CCN1, H2AX, IL10, WRN, HSP90AA1. Of these six genes only one, GJA1, was differentially expressed in the TNBC subtypes. We included these additional genes in our analysis. Gene Sets used for this study are listed in Table 1, for completeness we include sets for MSL and IM in this table but they were not used for further analysis.
Weaver gene-set identifier and the second column is the number of genes in the set and the third column is the gene-set title.
Gene set analysis
Gene sets were analyzed using the suite of tools from the GeneWeaver resource . As described above, we used the ‘HiSim Graph’ tool to enumerate and visualize intersections among the gene sets from human, mouse and rat, and the ‘Boolean Algebra’ tool to create a set of conserved genes representing the intersection of the homologs of the three sets. We used the ‘Jaccard Similarity’ tool to statistically evaluate and identify genes in the gene-set overlap between the set associated with cisplatin treatment, and sets of over- and under-expressed genes in the TNBC subtypes. We used the default parameters for all analysis tools, details of which can be found at the GeneWeaver.org website .
Gene function analysis
To identify processes enriched in gene sets and represent them in a graphical format we used the Visual Annotation Display (VLAD) tool for Gene Ontology enrichment analysis [34, 35]. First, to examine the 102 genes in the cisplatin-associated set we performed VLAD analysis to determine if those genes were enriched for processes known to represent cisplatin biology. We also tested the 20 cisplatin-associated genes that were differentially regulated in TNBC subtypes to see if their enrichment was different from the parental set, which would have indicated that those genes were enriched for a subset of processes that are involved in cisplatin biology. In all analyses, we used default parameters for VLAD enrichment analysis, and the set of UniProt-GOA human annotations as a background set . The analysis was run on September 2, 2019. The UniProt-GOA gene annotation data used was dated from February 26, 2018. Since GO annotations represent processes that occur in normal cells and we are ultimately interested in the effects these genes have with respect to cisplatin treatment, we extended the functional characterization of the cisplatin-associated genes that are differentially regulated in resistant TNBC subtypes by manually searching for evidence describing how they might contribute to cisplatin reistance or sensitivity.
An additional functional analysis was performed with the 102 genes in the cisplatin-associated set using the KEGG Mapper Search Pathway tool to interrogate Pathways and Diseases . Gene symbols were used with default parameters in the Organism-specific search mode (hsa). The analysis was performed on Sept 6, 2019.
We also ran an analysis using ‘String’, a network analysis tool that uses interaction data to functionally interrogate gene sets . The analysis was performed on Sept 8, 2019. Genes were entered using gene symbols, analysis in human was selected and all default parameters were used. GO and KEGG catagories are reported from the ‘Functional Analysis’ results.
Gene sets of differentially expressed genes in TNBC subtypes
To investigate sets of differentially regulated genes in TNBC subtypes, we created gene sets in GW for the six subtypes described by Lehmann et al . We chose these subtypes because the Lehmann analysis includes a measure of relative sensitivity to cisplatin treatment. Using the information from the supplemental data in Lehmann et al, we associated their gene symbols with unique HGNC identifiers to create 12 gene sets: i.e., an up and down-expressed set for each of the six TNBC subtypes (Table 1) . The gene sets ranged in size from 127 genes for which expression was down in the BL2 subtype, to 829 genes where expression is up in the MSL subtype. The 12 sets of up- and down-expressed genes represents 2161 unique human genes. Thirty-five genes were represented in 6 sets, and 101 genes were contained in 5 sets. One gene, KRT17 (HGNC:6427), was listed in both the up- and down-expressed sets of MSL. For further analysis, we focused on the four TNBC subtypes that represent subtypes that only contain tumor-derived cells .
Cisplatin-associated genes are enriched for processes that are consistent with the cytotoxic action and response to cisplatin
We hypothesized that by creating a gene set of evolutionarily-conserved cisplatin-interacting genes, we would select for genes that function in the fundamental actions of cisplatin. To test this, we used GO enrichment analysis to determine which biological processes were enriched in our 102 gene set. Our results confirm the validity of our strategy: we identified a set of genes that are involved in core cancer processes that are also known to be associated with action of cisplatin. Specifically, VLAD analysis showed that the 102 conserved cisplatin-associated genes were enriched for the GO biological processes: ‘aging’, ‘negative regulation of apoptotic process’, ‘apoptotic signaling pathway’, ‘response to ionizing radiation’, ‘cellular response to oxidative stress’, and ‘response to reactive oxygen species’ [Additional file 1: Table S1]. The 102 conserved genes were also enriched for the GO cellular component terms ‘chromosome, telomeric region’, ‘mitochondrion’, ‘cytosol’, ‘extracellular space’ and ‘membrane raft’ [Additional file 1: Table S1]. These results are consistent with the known mechanism of cisplatin action in which cisplatin causes oxidative stress, interacts with DNA and triggers a response that culminates in apoptosis .
We extended our GO results by interrogating the KEGG Pathway and KEGG Disease resources with the 102 conserved genes . The KEGG Disease analysis showed that our genes were most represented in a variety of different cancer types with esophageal cancer associated with the most genes  [Additional file 2: Table S2]. DNA excision repair was associated with four genes and breast cancer was associated with two. The top scorerer for the KEGG Pathway mapping analysis was ‘cancer pathways’ (36 genes) [Additional file 3: Table S3]. KEGG pathway analysis was also consistent with, and confirmed the GO enrichment analysis: apoptosis (27 genes), cellular senescence (21 genes) and stress response pathways like the P53 pathway (20 genes). The KEGG analysis also identified several viral pathways as well as the platinum drug resistance class (22 genes) [Additional file 3: Table S3].
The set was interrogated using the String Network analysis tool . Functional groupings from String were consistent with the VLAD and KEGG analysis results reported above [Additional file 4: Table S4].
A subset of cisplatin-associated differentially-expressed genes provide a signature for the resistant subtypes
Of the 102 evolutionarily conserved cisplatin-assocated genes, 20 are differentially expressed in TNBC subtypes (Table 2). Using the Jaccard Similarity Tool in GW, we compared the conserved set of cisplatin-responsive genes with the differentially expressed genes. Table 2 shows the summary of these data. Our results indicated that of the 102 cisplatin-associated genes conserved in human, mouse and rat, 16 genes were up-regulated in at least one of the four subtypes and 11 were down-regulated in at least one subtype.
Our results show that of the differentially expressed genes in each subtype, only a small proportion are associated with the set of cisplatin-interacting genes: BL1 (2:215 up and 5:251 down), BL2(4:154 up and 0:127 down), M(5:247 up and 3:446 down), and LAR (8:805 up and 3:382 down). If we examine only the set of genes that show different expression behavior in the resistant LAR and M subtypes when compared to the sensitive BL1 and BL2 subtypes, a signature of 13 genes is identified, shown in column 6 of Table 2. These results show that the differential expression of cisplatin-associated genes in breast cancer subtypes involves only a small percentage, 20 genes, of the overall genes used to characterize the subtypes and there is a set of 13 cisplatin-associated genes whose differential expression is characteristic of the two resistant subtypes.
The results of GO term enrichment analysis on the 20 differentially regulated genes for biological process are shown in Additional file 5: Table S5 [Additional file 5: Table S5]. Consistent with the conserved set of 102 cisplatin-associated genes, the 20 genes overlapping with the TNBC differentially regulated sets were also enriched for stress-response genes, aging, and genes that are involved in regulating programmed cell death. In addition terms representing the ‘CD95-death inducing complex’ and focal adhesion complexes were enriched, consistent with potential mechanisms of regulation of apoptosis and epithelial-to-mesenchymal transition mitochondrial outer membrane (p = 3.56e-05). Unlike the conserved set of genes, these 20 genes are not as significantly enriched for genes associated with telomeres (p = 1.1e-01) or nucleoplasm (p = 5.99e-02). This result shows that the subset of genes regulated in the TNBC subtypes are enriched for similar processes as the parental sets and are not biased towards other processes.
Genes that are differentially regulated in cisplatin-resistant TNBC subtypes identify a variety of mechanisms to escape cisplatin toxicity
To try to understand whether the differential regulation of the 13 cisplatin-associated genes in the LAR and M subtypes could explain the subtype’s resistance, we examined each gene individually to determine if there was evidence that the over- or under-expression of these genes correlated with resistance to cisplatin. The results of our analysis are shown in Table 3, where the LAR and M subtypes are shown to vary in their signature of cisplatin genes that are differentially regulated. Seven of the genes are exclusively differentially expressed in the LAR subtype, three in the M subtype and three are differentially expressed in both subtypes. Interestingly, the direction of the differential expression for the three common genes is in opposite directions. Examining how these genes might influence cisplatin reistance shows that, while some of the genes influence apoptosis directly, others identify different upstream mechanisms of resistance. Since cisplatin is not a first-line treatment for TNBC, most studies correlating these genes with resistance or sensitivity to cisplatin are from other cancer types. Our results suggest that these genes may also influence cisplatin sensitivity in TNBC, and may help further elucidate the mechanisms of cisplatin action in TNBC and suggest more refined strategies for cisplatin treatment.
We applied an integrated gene set analysis to identify potential biological mechanisms underlying cisplatin sensitivity in four different molecular subtypes of TNBC. We defined a set of 102 cisplatin-associated genes conserved across human, mouse, and rat, and we used knowledge about those genes to evaluate how those genes could be involved in the therapeutic response. Overall, our results show that many cisplatin-responsive genes are involved with the end stage of the effects of cisplatin treatment: cell death. Cell death is also the most globally differentially regulated process identified by cisplatin-responsive genes in all subtypes of TNBC. These results imply that agents that up-regulate apoptotic signaling, such as Trail sensitizers, should be investigated as effective global co-therapies for cisplatin treatment .
Response to Cisplatin
To specifically investigate the differences in cisplatin response with respect to each of the subtypes, we examined the genes that were uniquely differentially expressed in the resistant LAR and M subtypes. Response to cisplatin can be modulated by a number of different mechanisms: decreased cellular import or increased cellular efflux of cisplatin, detoxification of cisplatin, defective DNA repair or resistance to cell cycle arrest or cell death [25, 67, 68].
As noted previously, Lehmann et al showed that in cell lines, the BL1 and BL2 subtypes often contained mutations in one of the BRCA genes. They hypothesized that the DNA repair defect explained why BL1 and BL2 are more sensitive to cisplatin than the M or LAR subtypes. It has recently been suggested that platins or PARP inhibitors are potential treatment options for TNBC with BRCA mutations . A recent study by Zhao et al showed that other factors such as homologous recombination status may also influence cisplatin response in breast cancer . Our work suggests that in addition to BRCA mutation status, other factors may contribute to differential sensitivity of these subtypes. As described above, our results show that cisplatin-associated genes involved in cell death are differentially expressed in all TNBC subtypes, but the LAR and M subtypes have a unique signature of genes that are not differentially regulated in the same way in the BL1 or BL2 subtypes.
In particular, we find that the genes ABCC2, AKT1, BCL2L1, CASP8, CLU, GSR, NQO1 are up-regulated in the LAR subtype and MSH2 is downregulated. With the exception of CASP8, the regulation of all of these genes is consistent with reported resistance to cisplatin (Table 3). ABCC2 and GSR, specifically, represent a transporter and a glutathione metabolic enzyme, respectively, that lie in a pathway that detoxifies and transports cisplatin out of the cell [42, 54]. The increase in ABCC2 and GSR, and their actions upstream of the cell death related genes, provides us with a testable hypothesis for an additional mechanism that contributes to the relative cisplatin resistance of the LAR subtype compared to the other subtypes. That is to say, inhibition of either or both of these proteins could make LAR cells more sensitive to cisplatin treatment (Fig. 4). AKT1, CLU and NQO1 encode proteins that respond to stress, including oxidative stress, which is one of the mechanisms of cisplatin action . These three genes would contribute to cisplatin resitance in pathways downstream of GSR or ABCC2 [43, 49, 59,60,61]. BCL2L1 and CASP8 are both proteins integral to the apoptotic program. BCL2L1 is an inhibitor of apoptosis whose overexpression has been correlated with cisplatin-resistance, consistent with its upregulation in the resistant LAR subtype. The only down-regulated gene, MSH2, is a protein involved in DNA repair, although it has been shown to be necessary for the apoptotic action of cisplatin [56, 57]. The up-regulation of CASP8 is counter-indicative of cisplatin resistance, since its overexpression has been shown to make cells more sensitive to cisplatin . However, it is interesting to note that CASP8 would lie the most downstream of all of the other genes that are differentially regulated in the LAR subtype and therefore may be epistatically masked by upstream changes.
The LAR subtype also shows differential regulation of some genes also differentially regulated in the M subtype, but neither of the basal subtypes. VIM and TUBA1 are downregulated in the LAR subtype. Both VIM and TUBA1 have previously been associated with cisplatin resistance, but the causal effect remains to be determined [62, 64, 65]. In ovarian cancer cells down-regulation of VIM expression resulted in resistance to cisplatin by potentially down-regulating its import and up-regulating its export, indicating that it might also be contributing to cisplatin resistance in the LAR subtype . However, the factors controlling VIM expression and its exact role in cisplatin resistance in different cancer types are still not well understood. Some studies, including some breast-cancer studies show increased VIM expression correlates with cisplatin resistance [71,72,73]. One interesting question that arises from our analysis is whether or not the LAR subtype represents a heterogeneous population that can be further subdivided with respect to mechanisms of resistance and if so, what is the nature of the heterogeneity. Can some LAR tumors escape cisplatin by upregulating its transport out of the cell while others escape by different mechanisms such as upregulating GSR, or does a single tumor tend to accumulate multiple mechanisms of resistance? Because our analysis is retrospective and used aggregate data from previous studies, these types of questions require further investigation.
In the M subtype, some genes differentially regulated and potentially involved in cisplatin resistance differ from those identified in the LAR subtype. To fully understand the biology of cisplatin resistance in the M subtype, one area to further pursue is the epithelial-to-mesenchymal transition that results in increased VIM expression, which is downregulated in the LAR subtype.
The M subtype also shows differential up-regulation of VIM, NOX4 and TUBA1A. VCAM1 is downregulated in the M subtype. VCAM1 has also been associated with an increase in epithelial-to-mesenchymal transition and has been correlated with resistance to cisplatin [63, 64]. Overexpression of VCAM1 has been shown to confer cisplatin resistance in breast cancer cells . The downregulation of VCAM1 in the M subtype is counterintuitive to it being causative in this subtype’s lower sensitivity to cisplatin. As noted above, the expression of VIM is less well understood. Although overexpression correlates with cisplatin resistance in some contexts, it is still not well characterized mechanistically. At least two studies have shown that genes controlling the epithelial-to-mesenchymal transition, ITGB1 and TET1, confer cisplatin resistance, and those genes also increase the expression of VIM [65, 73]. The gene sets of TNBC differentially expressed genes did not include ITGB1 or TET1. NOX4 is an NADPH oxidase that generates reactive oxygen species and can make the effects of cisplatin treatment more severe. However, overexpression of NOX4 has been shown to result in normal breast cells being resistant to apoptosis . Like VCAM1, the higher differential expression of NOX4 is counterindictive for cisplatin resistance. CASP8 is also downregulated in the M subtype. In contrast to LAR, downregulation of CASP8 in the M subtype would lead to a defect in the apoptotic mechanism resulting in cisplatin resistance regardless of upstream triggers.
We have used a gene-set comparative approach to study potential mechnisms of cisplatin resitance in TNBC subtypes. Out results show that in the resistant LAR subtype a small number of genes that are differentially expressed identify a variety of potential mechanisms that can be used to escape cisplatin toxicity; transport, detoxification, and direct and indirect involvement in programmed cell death. We hypothesize that the differential expression of these genes impacts how tumors of a given subtype will respond to the agent. In the resistant M subtype, a small number of genes is also differentially regulated, but the interpretation of their contribution to resistance is less clear. The differentially regulated genes in the M subtype identify the epithelial-to-mesenchymal transition and the control of reactive oxygen species as potential regulators of cisplatin response.
By focusing on genes known to be associated with cisplatin, our method identifies (or excludes) genes that can serve as a signature in the differential response of TNBC subtypes to cisplatin treatment. This gives an advantage over global gene expression classification systems in that we can pinpoint specific gene signatures that classify with respect to a targeted area of interest, in this case with cisplatin association. Our results suggest that additional therapies to enhance the apoptotic mechanism might be globally beneficial for the treatment of all types of TNBC, while the LAR subtype might benefit from a combination treatment of cisplatin and glutathione-modulator treatment agents . For TNBC the analysis could be extended to investigate the molecular basis of the differences in response to other primary therapeutic agents such as taxane and anthracycline. One limitation to this extension is availability of data for analysis. These types of studies require existing experimental data with respect to response status and gene expression patterns for analysis and require high quality gene-chemical association data. In our study, we used existing data reported for TNBC subtypes and from the CTD resource to seed our analysis. As mentioned earlier, a limitation to this type of aggregate data is that it does not allow us to ask questions with respect to whether or not individual tumors or individual cells express different subsets of genes that confer resistance. These types of questions can be addressed in future studies in which wet-bench studies of expression from tumor samples or individual tumor cells are correlated with drug resitance or sensitivity and are analyzed in the context of high quality curated data about gene-chemical interactions. Ideally, a prospective strategy using markers such as BRCA status or PD-L1 to predict response-type would be most useful in deciding treatment options . Our results identify genes that can be further studied as useful biomarkers.
Availability of data and materials
The datasets generated and/or analysed during the current study are available in the GW repository .
Basal-like 1 subtype of Triple Negative Breast Cancer
Basal-like 2 subtype of Triple Negative Breast Cancer
Chemicals of Biological Interest
Immunomodulatory subtype of Triple Negative Breast Cancer
Luminal androgen receptor subtype of Triple Negative Breast Cancer
Mesenchymal subtype of Triple negative Breast Cancer
Mesenchymal stem-like subtype of Triple Negative Breast Cancer
Triple Negative Breast Cancer
Visual Annotation Display
Saraiva DP, Guadalupe Cabral M, Jacinto A, Braga S. How many diseases is triple negative breast cancer: the protagonism of the immune microenvironment. ESMO Open. 2017;2(4):e000208.
Alexandrov LB, Nik-Zainal S, Wedge DC, Aparicio SA, Behjati S, Biankin AV, et al. Signatures of mutational processes in human cancer. Nature. 2013;500(7463):415–21.
Menghi F, Inaki K, Woo X, Kumar PA, Grzeda KR, Malhotra A, et al. The tandem duplicator phenotype as a distinct genomic configuration in cancer. Proc Natl Acad Sci U S A. 2016;113(17):E2373–82.
Lehmann BD, Bauer JA, Chen X, Sanders ME, Chakravarthy AB, Shyr Y, et al. Identification of human triple-negative breast cancer subtypes and preclinical models for selection of targeted therapies. J Clin Invest. 2011;121(7):2750–67.
Burstein MD, Tsimelzon A, Poage GM, Covington KR, Contreras A, Fuqua SA, et al. Comprehensive genomic analysis identifies novel subtypes and targets of triple-negative breast cancer. Clin Cancer Res. 2015;21(7):1688–98.
Jhan JR, Andrechek ER. Triple-negative breast cancer and the potential for targeted therapy. Pharmacogenomics. 2017;18(17):1595–609.
Kibbe WA, Arze C, Felix V, Mitraka E, Bolton E, Fu G, et al. Disease ontology 2015 update: an expanded and updated database of human diseases for linking biomedical knowledge through disease data. Nucleic Acids Res. 2015;43(Database issue):D1071–8.
Dent R, Trudeau M, Pritchard KI, Hanna WM, Kahn HK, Sawka CA, et al. Triple-negative breast cancer: clinical features and patterns of recurrence. Clin Cancer Res. 2007;13(15 Pt 1):4429–34.
Liedtke C, Mazouni C, Hess KR, Andre F, Tordai A, Mejia JA, et al. Response to neoadjuvant therapy and long-term survival in patients with triple-negative breast cancer. J Clin Oncol. 2008;26(8):1275–81.
Perou CM, Sorlie T, Eisen MB, van de Rijn M, Jeffrey SS, Rees CA, et al. Molecular portraits of human breast tumours. Nature. 2000;406(6797):747–52.
Prat A, Adamo B, Cheang MC, Anders CK, Carey LA, Perou CM. Molecular characterization of basal-like and non-basal-like triple-negative breast cancer. Oncologist. 2013;18(2):123–33.
Lehmann BD, Jovanovic B, Chen X, Estrada MV, Johnson KN, Shyr Y, et al. Refinement of triple-negative breast Cancer molecular subtypes: implications for Neoadjuvant chemotherapy selection. PLoS One. 2016;11(6):e0157368.
Chacon RD, Costanzo MV. Triple-negative breast cancer. Breast Cancer Res. 2010;12 Suppl 2:S3.
Salgia R, Kulkarni P. The genetic/non-genetic duality of drug 'Resistance' in Cancer. Trends Cancer. 2018;4(2):110–8.
Wahba HA, El-Hadaad HA. Current approaches in treatment of triple-negative breast cancer. Cancer Biol Med. 2015;12(2):106–16.
Omarini C, Guaitoli G, Pipitone S, Moscetti L, Cortesi L, Cascinu S, et al. Neoadjuvant treatments in triple-negative breast cancer patients: where we are now and where we are going. Cancer Manag Res. 2018;10:91–103.
Caparica R, Lambertini M, de Azambuja E. How I treat metastatic triple-negative breast cancer. ESMO Open. 2019;4(Suppl 2):e000504.
Petrelli F, Barni S, Bregni G, de Braud F, Di Cosimo S. Platinum salts in advanced breast cancer: a systematic review and meta-analysis of randomized clinical trials. Breast Cancer Res Treat. 2016;160(3):425–37.
Rabik CA, Dolan ME. Molecular mechanisms of resistance and toxicity associated with platinating agents. Cancer Treat Rev. 2007;33(1):9–23.
Stordal B, Davey M. Understanding cisplatin resistance using cellular models. IUBMB Life. 2007;59(11):696–9.
Rabik CA, Maryon EB, Kasza K, Shafer JT, Bartnik CM, Dolan ME. Role of copper transporters in resistance to platinating agents. Cancer Chemother Pharmacol. 2009;64(1):133–42.
Zhu Z, Du S, Du Y, Ren J, Ying G, Yan Z. Glutathione reductase mediates drug resistance in glioblastoma cells by regulating redox homeostasis. J Neurochem. 2018;144(1):93–104.
Zamble DB, Mu D, Reardon JT, Sancar A, Lippard SJ. Repair of cisplatin--DNA adducts by the mammalian excision nuclease. Biochemistry. 1996;35(31):10004–13.
Galluzzi L, Senovilla L, Vitale I, Michels J, Martins I, Kepp O, et al. Molecular mechanisms of cisplatin resistance. Oncogene. 2012;31(15):1869–83.
Baker E, Bubier JA, Reynolds T, Langston MA, Chesler EJ. GeneWeaver: data driven alignment of cross-species genomics in biology and disease. Nucleic Acids Res. 2016;44(D1):D555–9.
Schriml LM, Mitraka E. The disease ontology: fostering interoperability between biological and clinical human disease-related data. Mamm Genome. 2015;26(9–10):584–9.
Groza T, Kohler S, Moldenhauer D, Vasilevsky N, Baynam G, Zemojtel T, et al. The human phenotype ontology: semantic unification of common and rare disease. Am J Hum Genet. 2015;97(1):111–24.
Hastings J, Owen G, Dekker A, Ennis M, Kale N, Muthukrishnan V, et al. ChEBI in 2016: improved services and an expanding collection of metabolites. Nucleic Acids Res. 2016;44(D1):D1214–9.
Baker EJ, Jay JJ, Bubier JA, Langston MA, Chesler EJ. GeneWeaver: a web-based system for integrative functional genomics. Nucleic Acids Res. 2012;40(Database issue):D1067–76.
Davis AP, Grondin CJ, Johnson RJ, Sciaky D, McMorran R, Wiegers J, et al. The comparative Toxicogenomics database: update 2019. Nucleic Acids Res. 2019;47(D1):D948–D54.
Davis AP, Grondin CJ, Johnson RJ, Sciaky D, King BL, McMorran R, et al. The comparative Toxicogenomics database: update 2017. Nucleic Acids Res. 2017;45(D1):D972–D8.
Richardson JE, Bult CJ. Visual annotation display (VLAD): a tool for finding functional themes in lists of genes. Mamm Genome. 2015;26(9–10):567–73.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. Gene Ontology Consortium Nat Genet. 2000;25(1):25–9.
Huntley RP, Sawford T, Mutowo-Meullenet P, Shypitsyna A, Bonilla C, Martin MJ, et al. The GOA database: gene ontology annotation updates for 2015. Nucleic Acids Res. 2015;43(Database issue):D1057–63.
Meng XM, Ren GL, Gao L, Yang Q, Li HD, Wu WF, et al. NADPH oxidase 4 promotes cisplatin-induced acute kidney injury via ROS-mediated programmed cell death and inflammation. Lab Investig. 2018;98(1):63–78.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–D13.
Yates B, Braschi B, Gray KA, Seal RL, Tweedie S, Bruford EA. Genenames.org: the HGNC and VGNC resources in 2017. Nucleic Acids Res. 2017;45(D1):D619–D25.
Dasari S, Tchounwou PB. Cisplatin in cancer therapy: molecular mechanisms of action. Eur J Pharmacol. 2014;740:364–78.
Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K. KEGG: new perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 2017;45(D1):D353–D61.
Tian J, Xu YY, Li L, Hao Q. MiR-490-3p sensitizes ovarian cancer cells to cisplatin by directly targeting ABCC2. Am J Transl Res. 2017;9(3):1127–38.
Liu LZ, Zhou XD, Qian G, Shi X, Fang J, Jiang BH. AKT1 amplification regulates cisplatin resistance in human lung cancer cells through the mammalian target of rapamycin/p70S6K1 pathway. Cancer Res. 2007;67(13):6325–32.
Fu X, Tian J, Zhang L, Chen Y, Hao Q. Involvement of microRNA-93, a new regulator of PTEN/Akt signaling pathway, in regulation of chemotherapeutic drug cisplatin chemosensitivity in ovarian cancer cells. FEBS Lett. 2012;586(9):1279–86.
Williams J, Lucas PC, Griffith KA, Choi M, Fogoros S, Hu YY, et al. Expression of Bcl-xL in ovarian carcinoma is associated with chemoresistance and recurrent disease. Gynecol Oncol. 2005;96(2):287–95.
Keitel U, Scheel A, Thomale J, Halpape R, Kaulfuss S, Scheel C, et al. Bcl-xL mediates therapeutic resistance of a mesenchymal breast cancer cell subpopulation. Oncotarget. 2014;5(23):11778–91.
Paul I, Chacko AD, Stasik I, Busacca S, Crawford N, McCoy F, et al. Acquired differential regulation of caspase-8 in cisplatin-resistant non-small-cell lung cancer. Cell Death Dis. 2012;3:e449.
Blanc C, Deveraux QL, Krajewski S, Janicke RU, Porter AG, Reed JC, et al. Caspase-3 is essential for procaspase-9 processing and cisplatin-induced apoptosis of MCF-7 breast cancer cells. Cancer Res. 2000;60(16):4386–90.
Ma G, Cai H, Gao L, Wang M, Wang H. sCLU regulates cisplatin chemosensitivity of lung cancer cells in vivo. World J Surg Oncol. 2015;13:80.
Zhang B, Zhang K, Liu Z, Hao F, Wang M, Li X, et al. Secreted clusterin gene silencing enhances chemosensitivity of a549 cells to cisplatin through AKT and ERK1/2 pathways in vitro. Cell Physiol Biochem. 2014;33(4):1162–75.
Wu W, Wang HD, Guo W, Yang K, Zhao YP, Jiang YG, et al. Up-regulation of Fas reverses cisplatin resistance of human small cell lung cancer cells. J Exp Clin Cancer Res. 2010;29:49.
Yang F, Long W, Xuechuan H, Xueqin L, Hongyun M, Yonghui D. Upregulation of Fas in epithelial ovarian cancer reverses the development of resistance to cisplatin. BMB Rep. 2015;48(1):30–5.
Shiu LY, Chang LC, Liang CH, Huang YS, Sheu HM, Kuo KW. Solamargine induces apoptosis and sensitizes breast cancer cells to cisplatin. Food Chem Toxicol. 2007;45(11):2155–64.
Chen HH, Kuo MT. Role of glutathione in the regulation of Cisplatin resistance in cancer chemotherapy. Met Based Drugs. 2010;2010.
Lien EC, Lyssiotis CA, Juvekar A, Hu H, Asara JM, Cantley LC, et al. Glutathione biosynthesis is a metabolic vulnerability in PI(3)K/Akt-driven breast cancer. Nat Cell Biol. 2016;18(5):572–8.
Topping RP, Wilkinson JC, Scarpinato KD. Mismatch repair protein deficiency compromises cisplatin-induced apoptotic signaling. J Biol Chem. 2009;284(21):14029–39.
Sawant A, Kothandapani A, Zhitkovich A, Sobol RW, Patrick SM. Role of mismatch repair proteins in the processing of cisplatin interstrand cross-links. DNA Repair (Amst). 2015;35:126–36.
Yu Y, Wang Y, Ren X, Tsuyada A, Li A, Liu LJ, et al. Context-dependent bidirectional regulation of the MutS homolog 2 by transforming growth factor beta contributes to chemoresistance in breast cancer cells. Mol Cancer Res. 2010;8(12):1633–42.
Kim HJ, Oh GS, Shen A, Lee SB, Choe SK, Kwon KB, et al. Augmentation of NAD(+) by NQO1 attenuates cisplatin-mediated hearing impairment. Cell Death Dis. 2014;5:e1292.
Gang GT, Kim YH, Noh JR, Kim KS, Jung JY, Shong M, et al. Protective role of NAD(P)H:quinone oxidoreductase 1 (NQO1) in cisplatin-induced nephrotoxicity. Toxicol Lett. 2013;221(3):165–75.
Syu JP, Chi JT, Kung HN. Nrf2 is the key to chemotherapy resistance in MCF7 breast cancer cells under hypoxia. Oncotarget. 2016;7(12):14659–72.
Toshimitsu H, Hashimoto K, Tangoku A, Iizuka N, Yamamoto K, Kawauchi S, et al. Molecular signature linked to acquired resistance to cisplatin in esophageal cancer cells. Cancer Lett. 2004;211(1):69–78.
Wang PC, Weng CC, Hou YS, Jian SF, Fang KT, Hou MF, et al. Activation of VCAM-1 and its associated molecule CD44 leads to increased malignant potential of breast cancer cells. Int J Mol Sci. 2014;15(3):3560–79.
Huo Y, Zheng Z, Chen Y, Wang Q, Zhang Z, Deng H. Downregulation of vimentin expression increased drug resistance in ovarian cancer cells. Oncotarget. 2016;7(29):45876–88.
Han X, Zhou Y, You Y, Lu J, Wang L, Hou H, et al. TET1 promotes cisplatin-resistance via demethylating the vimentin promoter in ovarian cancer. Cell Biol Int. 2017;41(4):405–14.
Taylor DJ, Parsons CE, Han H, Jayaraman A, Rege K. Parallel screening of FDA-approved antineoplastic drugs for identifying sensitizers of TRAIL-induced apoptosis in cancer cells. BMC Cancer. 2011;11:470.
Galluzzi L, Vitale I, Michels J, Brenner C, Szabadkai G, Harel-Bellan A, et al. Systems biology of cisplatin resistance: past, present and future. Cell Death Dis. 2014;5:e1257.
Siddik ZH. Cisplatin: mode of cytotoxic action and molecular basis of resistance. Oncogene. 2003;22(47):7265–79.
Zhao EY, Shen Y, Pleasance E, Kasaian K, Leelakumari S, Jones M, et al. Homologous recombination deficiency and platinum-based therapy outcomes in advanced breast Cancer. Clin Cancer Res. 2017;23(24):7521–30.
Yu W, Chen Y, Dubrulle J, Stossi F, Putluri V, Sreekumar A, et al. Cisplatin generates oxidative stress which is accompanied by rapid shifts in central carbon metabolism. Sci Rep. 2018;8(1):4306.
Voon DC, Huang RY, Jackson RA, Thiery JP. The EMT spectrum and therapeutic opportunities. Mol Oncol. 2017;11(7):878–91.
Haslehurst AM, Koti M, Dharsee M, Nuin P, Evans K, Geraci J, et al. EMT transcription factors snail and slug directly contribute to cisplatin resistance in ovarian cancer. BMC Cancer. 2012;12:91.
Yin HL, Wu CC, Lin CH, Chai CY, Hou MF, Chang SJ, et al. beta1 Integrin as a Prognostic and Predictive Marker in Triple-Negative Breast Cancer. Int J Mol Sci. 2016;17(9).
Graham KA, Kulawiec M, Owens KM, Li X, Desouki MM, Chandra D, et al. NADPH oxidase 4 is an oncoprotein localized to mitochondria. Cancer Biol Ther. 2010;10(3):223–31.
Jia Y, Zhang C, Zhou L, Xu H, Shi Y, Tong Z. Micheliolide overcomes KLF4-mediated cisplatin resistance in breast cancer cells by downregulating glutathione. Onco Targets Ther. 2015;8:2319–27.
This work was supported by NIH NHGRI grant R25 HG007053 Diversity action Plan for Mouse Genome Database (C. Bult, PI); NIH NCI grant P30CA034196 Jackson Laboratory Cancer Center (E. Liu, PI); NIH NHGRI U41 HG000330 Mouse Genome Informatics (C. Bult and J. Blake PIs); and NIH AA018776 (E. Chesler, PI).
The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
DPH and AH performed all gene-set analyses. DPH, AH, JM, MSM, SMM, TR, SP, and JAB contributed to curation of gene sets into the GW resource. All authors contributed to the analysis and interpretation of the results. All authors contributed to the writing and/or reviewing of the manuscript. All authors read and approved the final manuscript
Ethics approval and consent to participate
Consent for publication
The authors declare they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Gene Ontology Terms enriched in the 102 cisplatin-associated genes Description of data: The VLAD graphical output for GO Biological Process and GO Cellular Component was examined and reported in tabular format. The five most specific terms and their respective p-values are listed. The analysis was run on September 2, 2019. The UniProt-GOA gene annotation data used was dated from February 26, 2018.
KEGG-Disease analysis of 102 cisplatin associated genes Description of data: A list of the disease categories in KEGG that were associated with the 102 cisplatin-associated genes. Gene symbols were used with default parameters in the Organism-specific search mode (hsa). The analysis was performed on Sept 6, 2019. Categories with at least 2 genes are shown.
KEGG-Pathways analysis of 102 cisplatin associated genes Description of data: A list of the pathway categories in KEGG that were associated with the 102 cisplatin-associated genes. Gene symbols were used with default parameters in the Organism-specific search mode (hsa). The analysis was performed on Sept 6, 2019. Categories with > 10 genes are shown.
String network analysis of 102 cisplatin associated genes Description of data: A list of GO terms and KEGG categories identified using the String network analysis tool with respective False-discovery rates. The analysis was performed on Sept 8, 2019.
Gene Ontology Terms enriched in the 20 cisplatin-associated genes Description of data: The VLAD graphical output for GO Biological Process and GO Cellular Component was examined and reported in tabular format. The five most specific terms and their respective p-values are listed. The analysis was run on September 2, 2019. The UniProt-GOA gene annotation data used was dated from February 26, 2018.
About this article
Cite this article
Hill, D.P., Harper, A., Malcolm, J. et al. Cisplatin-resistant triple-negative breast cancer subtypes: multiple mechanisms of resistance. BMC Cancer 19, 1039 (2019). https://doi.org/10.1186/s12885-019-6278-9
- Triple-negative breast cancer
- Cisplatin sensitivity
- Cancer subtypes
- Gene expression
- Cancer genomics
- Drug response
- Functional genomics
- Data mining