Mol. Cells 2019; 42(3): 237-244
Published online February 12, 2019
https://doi.org/10.14348/molcells.2018.0413
© The Korean Society for Molecular and Cellular Biology
Correspondence to : *Correspondence: kimsy@kribb.re.kr (SYK); wkim@ewha.ac.kr (WK)
Understanding the mechanisms of cancer drug resistance is a critical challenge in cancer therapy. For many cancer drugs, various resistance mechanisms have been identified such as target alteration, alternative signaling pathways, epithelial–mesenchymal transition, and epigenetic modulation. Resistance may arise via multiple mechanisms even for a single drug, making it necessary to investigate multiple independent models for comprehensive understanding and therapeutic application. In particular, we hypothesize that different resistance processes result in distinct gene expression changes. Here, we present a web-based database, CDRgator (Cancer Drug Resistance navigator) for comparative analysis of gene expression signatures of cancer drug resistance. Resistance signatures were extracted from two different types of datasets. First, resistance signatures were extracted from transcriptomic profiles of cancer cells or patient samples and their resistance-induced counterparts for >30 cancer drugs. Second, drug resistance group signatures were also extracted from two large-scale drug sensitivity datasets representing ~1,000 cancer cell lines. All the datasets are available for download, and are conveniently accessible based on drug class and cancer type, along with analytic features such as clustering analysis, multidimensional scaling, and pathway analysis. CDRgator allows meta-analysis of independent resistance models for more comprehensive understanding of drug-resistance mechanisms that is difficult to accomplish with individual datasets alone (database URL:
Keywords cancer drug resistance, gene expression signatures, meta-analysis, microarray, RNA-seq analysis, transcriptome
Overcoming drug resistance is one of the most critical challenges in cancer therapy. Cancer drug resistance appears not only in patients who have genetic factors interfering with drug actions, but is also induced by long-term treatment even in initially responsive patients. Because such acquired resistance leads to significant survival disadvantages, it is important to understand the underlying mechanisms of drug resistance in different patients.
In the early studies on resistance mechanisms, genetic mutations that directly alter drug target proteins (i.e. T790M in EGFR and L1152R in ALK) received significant attention (Choi et al., 2010; Pao et al., 2005). These mutations at gatekeeper residues are a common mechanism of resistance to drugs targeting oncogenic kinases, and many patients develop resistance within a year (Holohan et al., 2013). Recently, mutations other than those in the target protein have also been actively studied including alteration of regulatory regions (Leucci et al., 2018; Melton et al., 2017). Although a number of studies on drug resistance have focused on mutations, it has become clear that various genetic events prior to the acquisition of resistance also play an important role (Leucci et al., 2018).
In addition to mutations in the target, other genetic or epigenetic changes including copy number variation (CNV, Bean et al., 2007), activation of bypass signaling (Niederst and Engelman, 2013), and epigenetic modulations (Hu and Baeg, 2017; Nyce et al., 1993) have been identified as primary causes of resistance. Resistance mechanisms are not limited to alterations in only a few genes, but involve overall cellular changes, e.g. epithelial-mesenchymal–transition (EMT) or transition to cancer stem cells (CSCs), which seem to be important factors (Huang et al., 2015; Phi et al., 2018; Wang et al., 2016b). Because such major changes in cell state and physiology cannot be represented by a few genetic alterations, investigation of genome-scale signatures such as transcriptome profiles is likely to be more informative. Transcriptome profiles reflect changes in both genetic and epigenetic conditions including disease progression, drug perturbation, and drug resistance. Therefore, gene expression signatures are useful markers to investigate drug resistance states, classify distinct resistance mechanisms, and predict treatment outcomes.
Typically, there are two strategies to study genetic features associated with drug sensitivity. One is to use large-scale pharmacogenomic datasets of cancer cell lines as model systems, with the cells being divided into two groups, i.e. sensitive and resistant, and genomic features differentiating the two groups being identified. The other is to culture cancer cells in the presence of a drug, leading to selection of resistant cells. By comparing the molecular profiles of the resistance-induced cells with those of the original ones, the relevant genomic features can be analyzed more explicitly than by using the group-based method because the origins of the cells are the same in the former method.
There are several publicly available pharmacogenomics databases such as Genomics of Drug Sensitivity in Cancer (GDSC,
In the present study, we present a web platform, CDRgator (Cancer Drug Resistance signature navigator), which is among the most comprehensive databases of acquired cancer drug resistance because it amasses both
Currently, CDRgator provides 603 resistance signatures for 37 cancer drugs representing more than 26 cancer types in total, and the number of signatures will grow as more data are collected. It allows users to browse resistance signatures based on on drugs or cancer types, and to analyze the similarity between resistance signatures. Additionally, CDRgator has a tool to identify resistance signature-matched gene sets from Kyoto Encyclopedia of Genes and Genomes (KEGG, (Kanehisa, 2004) or Gene ontology (GO) (Carbon et al., 2017) to characterize the biological processes involved in resistance. CDRgator also has the ability to compare the resistance signatures in a database with signatures input by users to filter out drugs expected to be resistant or ineffective. Using an illustrative analysis of EGFR inhibitors, we show the utility of CDRgator in understanding the diverse mechanisms of cancer drug resistance.
To obtain induced drug resistance signatures (Sind), we manually collected datasets containing transcriptomic profiles of resistance-induced cells in GEO and ArrayExpress (Fig. 1). The datasets were filtered based on the following criteria; (1) the presence of matched pairs of sensitive and resistant cells, (2) treatment with monotherapy but not combinational therapy for a specific period of time, (3) gene transcription quantification using RNA-seq or microarray (4) in human cells and not in mouse cells or xenografted mouse cells. The RNA-seq data (.fastq file) were mapped using STAR aligner (Dobin et al., 2013) and quantified using HTSeq (Anders et al., 2015). The list of the collected datasets is available with detailed information in Additional file 1.
For identifying the resistant group signature (Sgrp), we used drug sensitivity data and gene expression data of cancer cell lines from CTRP and CCLE, respectively. To designate the resistant cell line group for a given drug, we grouped cancer cell lines based on their drug responses. The drug responses were represented as area under curve (AUC) of cell growth at different drug concentrations, and normalized with the maximum area calculated assuming 100% response in given concentration ranges. Then, we performed Z-transformation of logged AUC values (Z-AUC) to specify resistant or sensitive cells. We defined cells as a resistant cell group if their Z-AUC was less than −0.5. For all drugs, we generated signatures only in cases with at least three cells in the resistant group.
The expression data were downloaded from GEO or ArrayExpress, and were processed using our protocols for microarray and RNAseq analyses. Background correction and normalization was performed using R. For RNAseq, raw fastq files were downloaded and aligned to the reference genome (GRCh37) using STAR aligner. For RANseq of cancer cell lines, BAM files were downloaded from Genomic Data Commons (GDC). To estimate expression, RNAseq read count was calculated using HTseq. The gene expression signature was obtained using differentially expressed gene (DEG) analysis of the resistant cell group or the resistance-induced cells. For datasets containing biological replicates, genes with adjusted p-values less than 0.05 were selected as DEGs using the limma package (Ritchie et al., 2015) for microarray, and the DESeq2 package (Love et al., 2014) for RNAseq. For datasets with no biological replicates, we calculated fold change in gene expression values and transformed them to z-score. Genes with absolute fold change in z-score of more than 2 were selected as DEGs.
CDRgator contains a total of 143 Sind for 37 anti-cancer drugs and 30 cell lines from 16 tissue origins. A set of 499 Sgrp in our study represented 19 drugs which were also included in Sind; in terms of details, 267 Sgrp were from 20 tissue-level cancer types (e.g. blood cancer) and 232 Sgrp were from 23 disease level-cancer types (e.g. acute lymphoblastic leukemia).
CDRgator provides similarity analysis which explores how similar two signatures are; an enrichment factor (
where
To annotate the biological functions of gene signatures, an over-representation test was conducted using the ‘cluster-Profiler’ in the R package (Yu et al., 2012), which implements a hypergeometric test of overrepresentation of a gene set against background. Each signature analysis was performed with pathway terms consisting of 10–500 genes among the KEGG pathways or GO biological processes. As the size of the signature set decreases, the number of significantly enriched pathways decreases. Thus, CDRgator also generates a p-value representing the significance of functional enrichment.
The CDRgator web service is freely available at
CDRgator is primarily composed of ‘Browse’ and ‘Analysis’ menus. The ‘Browse’ menu provides the ability to search and browse information regarding resistance signatures. The ‘Analysis’ menu performs signature similarity and functional analysis. Signature similarity analysis calculates how similar resistance signatures are, so that users can identify if signatures of specific cancers, drugs, or their classes show a significant tendency in original spaces or two dimensions through multidimensional scaling. A functional analysis shows a heatmap that is generated from an over-representation test to identify genetic and biological functions of signatures. Finally, significant genes can be visualized using KEGG pathways (
To demonstrate the utility of CDRgator, we performed a case study of induced EGFR inhibitor resistance signatures. There were 22 induced resistance signatures for five different EGFR inhibitors. Because these signatures were generated independent of each other, the following questions could be asked; i) is the resistance pattern dependent on the drug, the cancer type, or both; and ii) do drugs with similar mode-of-action (e.g. EGFR inhibitors) share changes in transcriptomic profiles.
Similarity clustering analysis of the 22 Sind for EGFR inhibitors showed a tight cluster, with mostly erlotinib resistance signatures from head and neck cancer (Fig. 3A). This suggests distinct and reproducible patterns of resistance mechanisms specific to a drug and cancer type, because similar signature patterns were generated from multiple independent experiments. This trend was also confirmed in 2D mapping of the signatures using multidimensional scaling (Fig. 3B). These signatures were distinguishable from those of other cancer types. Overall, the resistance patterns collected in CDRgator appear to be heterogeneous even for the same drug and cancer type (e.g. erlotinib (center green) in non-small cell lung cancinoma (NSCLC, green outer ring) in Fig. 3B).
Pathway analysis showed the heterogeneity in resistance signatures in more detail. In head and neck cancers, biological pathways previously known to be associated with drug resistance were also enriched in our signatures, such as up-regulation of NF-kB, NOD-like receptor, TNF, and FoxO signaling pathways. NF-kB signaling is well-known among the survival and resistance pathways against many anticancer drugs including EGFR-targeted drugs (Bentires-Alj et al., 2003; Hertlein and Byrd, 2010; Lagunas and Meléndez-Zajgla, 2008; Wang et al., 2018). It is known that the NF-kB activating complex drives cancer cells to become resistant against EGFR inhibitors (Blakely et al., 2015), and this signaling pathway could also lead to resistance against EGFR inhibitors designed to target mutated forms of EGFR (Galvani et al., 2015). In colon cancers, activated NF-kB induces drug resistance through the regulation of MDR1 gene expression (Bentires-Alj et al., 2003). EGFR triggers proto-oncogenic signals such as RAS and mitogen activated protein kinase (MAPK) signaling (Wee and Wang, 2017). NOD-like receptors are known to regulate innate immunity by forming inflammasomes, and also by activating MAPK and NF-kB signaling (Saxena and Yeretssian, 2014; Shaw et al., 2010). NOD1/2 downstream signaling was reported to confer drug resistance through RIP2-PAX5 interaction by activating MAPK and NF-kB signaling (Wang et al., 2016a). AKT-PI3K signaling is also known to be activated in EGFR drug-resistant cancer cells (Jacobsen et al., 2017; Ma et al., 2016). FOXO is among the downstream genes playing a role in cell survival and apoptosis. Drug sensitivity was restored via inhibition of FOXO in AKT-mediated EGFR inhibitor-resistant lung cancers (Sangodkar et al., 2012). Several of these pathways (
In NSCLC, additional pathways (i.e. PI3K-Akt, ECM-receptor interaction, and focal adhesion) were strongly enriched. These pathways were found to be associated with epithelial–mesenchymal transition (EMT), which has been recently highlighted as major resistance mechanism in both chemotherapy and targeted cancer drug therapy (Chen et al., 2013; Huang et al., 2015; Larue and Bellacosa, 2005; Le Bras et al., 2012; Wang et al., 2016b). As shown in the example with EGFR inhibitors, CDRgator provides rich information on the drug resistance processes in terms of heterogeneity, and dependence on cancer type and drug class.
The development of technology coupled with bioinformatics and system biology approaches has enabled the identification of genomic and transcriptome features to predict response and resistance to specific drugs. The wide range of molecular mechanisms involving these features has important implications for understanding and treating resistance. However, inhibiting one pathway can result in a relatively simple escape route for the tumor, and an integrated analysis approach involving various features to overcome or circumvent drug resistance is needed.
CDRgator was developed for studying the acquired resistance mechanisms against cancer drugs. Gene expression signatures were generated from manually collected data regarding cancer drug resistance-induced cells in public datasets. CDRgator allows inspection and comparison of these acquired resistance transcriptomic signatures. Pathway analysis enables functional interpretation of multiple signatures simultaneously. While the datasets available in CDRgator are still sparse and insufficient to provide comprehensive information on drug resistance mechanisms, it may provide a useful clue to investigate drug resistance mechanisms, and to develop the right therapeutic strategy. For example, an analysis of genome-scale transcript signatures provided by CDRgator and recurrent mutations in the regulatory region may provide a variety of data enabling resistance marker identification or mechanistic analysis. The complete datasets in CDRgator are accessible at
Mol. Cells 2019; 42(3): 237-244
Published online March 31, 2019 https://doi.org/10.14348/molcells.2018.0413
Copyright © The Korean Society for Molecular and Cellular Biology.
Su-Kyeong Jang1,4, Byung-Ha Yoon2,3,4, Seung Min Kang1,4, Yeo-Gha Yoon1, Seon-Young Kim2,3,*, and Wankyu Kim1,*
1Ewha Research Center for Systems Biology, Department of Life Science, Division of Molecular & Life Sciences, Ewha Womans University, Seoul 03760, Korea, 2Gene Editing Research Center, KRIBB, Daejeon 34141, Korea, 3Department of Functional Genomics, University of Science and Technology (UST), Daejeon 34113, Korea
Correspondence to:*Correspondence: kimsy@kribb.re.kr (SYK); wkim@ewha.ac.kr (WK)
Understanding the mechanisms of cancer drug resistance is a critical challenge in cancer therapy. For many cancer drugs, various resistance mechanisms have been identified such as target alteration, alternative signaling pathways, epithelial–mesenchymal transition, and epigenetic modulation. Resistance may arise via multiple mechanisms even for a single drug, making it necessary to investigate multiple independent models for comprehensive understanding and therapeutic application. In particular, we hypothesize that different resistance processes result in distinct gene expression changes. Here, we present a web-based database, CDRgator (Cancer Drug Resistance navigator) for comparative analysis of gene expression signatures of cancer drug resistance. Resistance signatures were extracted from two different types of datasets. First, resistance signatures were extracted from transcriptomic profiles of cancer cells or patient samples and their resistance-induced counterparts for >30 cancer drugs. Second, drug resistance group signatures were also extracted from two large-scale drug sensitivity datasets representing ~1,000 cancer cell lines. All the datasets are available for download, and are conveniently accessible based on drug class and cancer type, along with analytic features such as clustering analysis, multidimensional scaling, and pathway analysis. CDRgator allows meta-analysis of independent resistance models for more comprehensive understanding of drug-resistance mechanisms that is difficult to accomplish with individual datasets alone (database URL:
Keywords: cancer drug resistance, gene expression signatures, meta-analysis, microarray, RNA-seq analysis, transcriptome
Overcoming drug resistance is one of the most critical challenges in cancer therapy. Cancer drug resistance appears not only in patients who have genetic factors interfering with drug actions, but is also induced by long-term treatment even in initially responsive patients. Because such acquired resistance leads to significant survival disadvantages, it is important to understand the underlying mechanisms of drug resistance in different patients.
In the early studies on resistance mechanisms, genetic mutations that directly alter drug target proteins (i.e. T790M in EGFR and L1152R in ALK) received significant attention (Choi et al., 2010; Pao et al., 2005). These mutations at gatekeeper residues are a common mechanism of resistance to drugs targeting oncogenic kinases, and many patients develop resistance within a year (Holohan et al., 2013). Recently, mutations other than those in the target protein have also been actively studied including alteration of regulatory regions (Leucci et al., 2018; Melton et al., 2017). Although a number of studies on drug resistance have focused on mutations, it has become clear that various genetic events prior to the acquisition of resistance also play an important role (Leucci et al., 2018).
In addition to mutations in the target, other genetic or epigenetic changes including copy number variation (CNV, Bean et al., 2007), activation of bypass signaling (Niederst and Engelman, 2013), and epigenetic modulations (Hu and Baeg, 2017; Nyce et al., 1993) have been identified as primary causes of resistance. Resistance mechanisms are not limited to alterations in only a few genes, but involve overall cellular changes, e.g. epithelial-mesenchymal–transition (EMT) or transition to cancer stem cells (CSCs), which seem to be important factors (Huang et al., 2015; Phi et al., 2018; Wang et al., 2016b). Because such major changes in cell state and physiology cannot be represented by a few genetic alterations, investigation of genome-scale signatures such as transcriptome profiles is likely to be more informative. Transcriptome profiles reflect changes in both genetic and epigenetic conditions including disease progression, drug perturbation, and drug resistance. Therefore, gene expression signatures are useful markers to investigate drug resistance states, classify distinct resistance mechanisms, and predict treatment outcomes.
Typically, there are two strategies to study genetic features associated with drug sensitivity. One is to use large-scale pharmacogenomic datasets of cancer cell lines as model systems, with the cells being divided into two groups, i.e. sensitive and resistant, and genomic features differentiating the two groups being identified. The other is to culture cancer cells in the presence of a drug, leading to selection of resistant cells. By comparing the molecular profiles of the resistance-induced cells with those of the original ones, the relevant genomic features can be analyzed more explicitly than by using the group-based method because the origins of the cells are the same in the former method.
There are several publicly available pharmacogenomics databases such as Genomics of Drug Sensitivity in Cancer (GDSC,
In the present study, we present a web platform, CDRgator (Cancer Drug Resistance signature navigator), which is among the most comprehensive databases of acquired cancer drug resistance because it amasses both
Currently, CDRgator provides 603 resistance signatures for 37 cancer drugs representing more than 26 cancer types in total, and the number of signatures will grow as more data are collected. It allows users to browse resistance signatures based on on drugs or cancer types, and to analyze the similarity between resistance signatures. Additionally, CDRgator has a tool to identify resistance signature-matched gene sets from Kyoto Encyclopedia of Genes and Genomes (KEGG, (Kanehisa, 2004) or Gene ontology (GO) (Carbon et al., 2017) to characterize the biological processes involved in resistance. CDRgator also has the ability to compare the resistance signatures in a database with signatures input by users to filter out drugs expected to be resistant or ineffective. Using an illustrative analysis of EGFR inhibitors, we show the utility of CDRgator in understanding the diverse mechanisms of cancer drug resistance.
To obtain induced drug resistance signatures (Sind), we manually collected datasets containing transcriptomic profiles of resistance-induced cells in GEO and ArrayExpress (Fig. 1). The datasets were filtered based on the following criteria; (1) the presence of matched pairs of sensitive and resistant cells, (2) treatment with monotherapy but not combinational therapy for a specific period of time, (3) gene transcription quantification using RNA-seq or microarray (4) in human cells and not in mouse cells or xenografted mouse cells. The RNA-seq data (.fastq file) were mapped using STAR aligner (Dobin et al., 2013) and quantified using HTSeq (Anders et al., 2015). The list of the collected datasets is available with detailed information in Additional file 1.
For identifying the resistant group signature (Sgrp), we used drug sensitivity data and gene expression data of cancer cell lines from CTRP and CCLE, respectively. To designate the resistant cell line group for a given drug, we grouped cancer cell lines based on their drug responses. The drug responses were represented as area under curve (AUC) of cell growth at different drug concentrations, and normalized with the maximum area calculated assuming 100% response in given concentration ranges. Then, we performed Z-transformation of logged AUC values (Z-AUC) to specify resistant or sensitive cells. We defined cells as a resistant cell group if their Z-AUC was less than −0.5. For all drugs, we generated signatures only in cases with at least three cells in the resistant group.
The expression data were downloaded from GEO or ArrayExpress, and were processed using our protocols for microarray and RNAseq analyses. Background correction and normalization was performed using R. For RNAseq, raw fastq files were downloaded and aligned to the reference genome (GRCh37) using STAR aligner. For RANseq of cancer cell lines, BAM files were downloaded from Genomic Data Commons (GDC). To estimate expression, RNAseq read count was calculated using HTseq. The gene expression signature was obtained using differentially expressed gene (DEG) analysis of the resistant cell group or the resistance-induced cells. For datasets containing biological replicates, genes with adjusted p-values less than 0.05 were selected as DEGs using the limma package (Ritchie et al., 2015) for microarray, and the DESeq2 package (Love et al., 2014) for RNAseq. For datasets with no biological replicates, we calculated fold change in gene expression values and transformed them to z-score. Genes with absolute fold change in z-score of more than 2 were selected as DEGs.
CDRgator contains a total of 143 Sind for 37 anti-cancer drugs and 30 cell lines from 16 tissue origins. A set of 499 Sgrp in our study represented 19 drugs which were also included in Sind; in terms of details, 267 Sgrp were from 20 tissue-level cancer types (e.g. blood cancer) and 232 Sgrp were from 23 disease level-cancer types (e.g. acute lymphoblastic leukemia).
CDRgator provides similarity analysis which explores how similar two signatures are; an enrichment factor (
where
To annotate the biological functions of gene signatures, an over-representation test was conducted using the ‘cluster-Profiler’ in the R package (Yu et al., 2012), which implements a hypergeometric test of overrepresentation of a gene set against background. Each signature analysis was performed with pathway terms consisting of 10–500 genes among the KEGG pathways or GO biological processes. As the size of the signature set decreases, the number of significantly enriched pathways decreases. Thus, CDRgator also generates a p-value representing the significance of functional enrichment.
The CDRgator web service is freely available at
CDRgator is primarily composed of ‘Browse’ and ‘Analysis’ menus. The ‘Browse’ menu provides the ability to search and browse information regarding resistance signatures. The ‘Analysis’ menu performs signature similarity and functional analysis. Signature similarity analysis calculates how similar resistance signatures are, so that users can identify if signatures of specific cancers, drugs, or their classes show a significant tendency in original spaces or two dimensions through multidimensional scaling. A functional analysis shows a heatmap that is generated from an over-representation test to identify genetic and biological functions of signatures. Finally, significant genes can be visualized using KEGG pathways (
To demonstrate the utility of CDRgator, we performed a case study of induced EGFR inhibitor resistance signatures. There were 22 induced resistance signatures for five different EGFR inhibitors. Because these signatures were generated independent of each other, the following questions could be asked; i) is the resistance pattern dependent on the drug, the cancer type, or both; and ii) do drugs with similar mode-of-action (e.g. EGFR inhibitors) share changes in transcriptomic profiles.
Similarity clustering analysis of the 22 Sind for EGFR inhibitors showed a tight cluster, with mostly erlotinib resistance signatures from head and neck cancer (Fig. 3A). This suggests distinct and reproducible patterns of resistance mechanisms specific to a drug and cancer type, because similar signature patterns were generated from multiple independent experiments. This trend was also confirmed in 2D mapping of the signatures using multidimensional scaling (Fig. 3B). These signatures were distinguishable from those of other cancer types. Overall, the resistance patterns collected in CDRgator appear to be heterogeneous even for the same drug and cancer type (e.g. erlotinib (center green) in non-small cell lung cancinoma (NSCLC, green outer ring) in Fig. 3B).
Pathway analysis showed the heterogeneity in resistance signatures in more detail. In head and neck cancers, biological pathways previously known to be associated with drug resistance were also enriched in our signatures, such as up-regulation of NF-kB, NOD-like receptor, TNF, and FoxO signaling pathways. NF-kB signaling is well-known among the survival and resistance pathways against many anticancer drugs including EGFR-targeted drugs (Bentires-Alj et al., 2003; Hertlein and Byrd, 2010; Lagunas and Meléndez-Zajgla, 2008; Wang et al., 2018). It is known that the NF-kB activating complex drives cancer cells to become resistant against EGFR inhibitors (Blakely et al., 2015), and this signaling pathway could also lead to resistance against EGFR inhibitors designed to target mutated forms of EGFR (Galvani et al., 2015). In colon cancers, activated NF-kB induces drug resistance through the regulation of MDR1 gene expression (Bentires-Alj et al., 2003). EGFR triggers proto-oncogenic signals such as RAS and mitogen activated protein kinase (MAPK) signaling (Wee and Wang, 2017). NOD-like receptors are known to regulate innate immunity by forming inflammasomes, and also by activating MAPK and NF-kB signaling (Saxena and Yeretssian, 2014; Shaw et al., 2010). NOD1/2 downstream signaling was reported to confer drug resistance through RIP2-PAX5 interaction by activating MAPK and NF-kB signaling (Wang et al., 2016a). AKT-PI3K signaling is also known to be activated in EGFR drug-resistant cancer cells (Jacobsen et al., 2017; Ma et al., 2016). FOXO is among the downstream genes playing a role in cell survival and apoptosis. Drug sensitivity was restored via inhibition of FOXO in AKT-mediated EGFR inhibitor-resistant lung cancers (Sangodkar et al., 2012). Several of these pathways (
In NSCLC, additional pathways (i.e. PI3K-Akt, ECM-receptor interaction, and focal adhesion) were strongly enriched. These pathways were found to be associated with epithelial–mesenchymal transition (EMT), which has been recently highlighted as major resistance mechanism in both chemotherapy and targeted cancer drug therapy (Chen et al., 2013; Huang et al., 2015; Larue and Bellacosa, 2005; Le Bras et al., 2012; Wang et al., 2016b). As shown in the example with EGFR inhibitors, CDRgator provides rich information on the drug resistance processes in terms of heterogeneity, and dependence on cancer type and drug class.
The development of technology coupled with bioinformatics and system biology approaches has enabled the identification of genomic and transcriptome features to predict response and resistance to specific drugs. The wide range of molecular mechanisms involving these features has important implications for understanding and treating resistance. However, inhibiting one pathway can result in a relatively simple escape route for the tumor, and an integrated analysis approach involving various features to overcome or circumvent drug resistance is needed.
CDRgator was developed for studying the acquired resistance mechanisms against cancer drugs. Gene expression signatures were generated from manually collected data regarding cancer drug resistance-induced cells in public datasets. CDRgator allows inspection and comparison of these acquired resistance transcriptomic signatures. Pathway analysis enables functional interpretation of multiple signatures simultaneously. While the datasets available in CDRgator are still sparse and insufficient to provide comprehensive information on drug resistance mechanisms, it may provide a useful clue to investigate drug resistance mechanisms, and to develop the right therapeutic strategy. For example, an analysis of genome-scale transcript signatures provided by CDRgator and recurrent mutations in the regulatory region may provide a variety of data enabling resistance marker identification or mechanistic analysis. The complete datasets in CDRgator are accessible at
Chunggab Choi, Seung-Hun Oh, Jeong-Eun Noh, Yong-Woo Jeong, Soonhag Kim, Jung Jae Ko, Ok-Joon Kim, and Jihwan Song
Mol. Cells 2016; 39(4): 337-344 https://doi.org/10.14348/molcells.2016.2317Sang-Choon Lee, Myung-Ho Lim, Jin A Kim, Soo-In Lee, Jung Sun Kim, Mina Jin, Soo-Jin Kwon, Jeong-Hwan Mun, Yeon-Ki Kim, Hyun Uk Kim, Yoonkang Hur and Beom-Seok Park
Mol. Cells 2008; 26(6): 595-605 https://doi.org/10.14348/.2008.26.6.595Dawon Hong and Sunjoo Jeong
Mol. Cells 2023; 46(1): 48-56 https://doi.org/10.14348/molcells.2023.0003