Large-scale pharmacogenomics based drug discovery for ITGB3 dependent chemoresistance in mesenchymal lung cancer
Molecular Cancer volume 17, Article number: 175 (2018)
Even when targets responsible for chemoresistance are identified, drug development is often hampered due to the poor druggability of these proteins. We systematically analyzed therapy-resistance with a large-scale cancer cell transcriptome and drug-response datasets and predicted the candidate drugs based on the gene expression profile. Our results implicated the epithelial–mesenchymal transition as a common mechanism underlying resistance to chemotherapeutic drugs. Notably, we identified ITGB3, whose expression was abundant in both drug resistance and mesenchymal status, as a promising target to overcome chemoresistance. We also confirmed that depletion of ITGB3 sensitized cancer cells to conventional chemotherapeutic drugs by modulating the NF-κB signaling pathway. Considering the poor druggability of ITGB3 and the lack of feasible drugs to directly inhibit this protein, we took an in silico screening for drugs mimicking the transcriptome-level changes caused by knockdown of ITGB3. This approach successfully identified atorvastatin as a novel candidate for drug repurposing, paving an alternative path to drug screening that is applicable to undruggable targets.
Recent studies in both in vitro cell and in vivo animal models demonstrated that the epithelial-mesenchymal transition (EMT), a major cause of metastasis, is closely associated with chemoresistance . These are consistent with the reports that cancer patients with mesenchymal gene signatures have poor prognoses or exhibit therapy resistance . However, due to the poor druggability of the EMT-associated proteins responsible for chemoresistance (e.g., ZEB1/2, SNAI2, SOX4, etc.), it is important to develop alternative strategies to make ‘undruggable but attractive targets’ druggable.
To elucidate the mechanisms underlying chemoresistance, we examined the gene expression profiles of 804 cancer cell lines, as well as their responses to anti-cancer drugs, using data from Cancer Therapeutics Response Portal (CTRP) (Additional file 1: Figure S1A). For each drug, cell lines were classified as resistant or sensitive group, and differentially expressed genes (DEGs) in each resistant group were selected (Additional file 1: Figure S1B). EMT was the most frequently up-regulated phathway in the resistant group across most chemotherapeutics (26 out of 32 drugs), and targeted drugs (15 out of 20 drugs) (Fig. 1a and Additional file 1: Figure S1C). Among the down-regulated genes, ‘immune & inflammatory response’ pathways were highly enriched (Additional file 1: Figure S1C). It would be noteworthy that upregulation of interferon signaling contributes to efficacy of chemotherapy . Given the previous studies supporting EMT as a major chemoresistance mechanism , we chose to investigate potential targets among the EMT signature genes, of which the vast majority were up-regulated by chemotherapeutics (Additional file 1: Figure S1C). Using the mesenchymal-type lung cancer cells (A549TD; hereafter, TD) from the A549 lung cancer cell line  showing a clear gene signature of ‘hallmark of EMT’ (Additional file 1: Figure S1E) and high resistance to etoposide treatment  (Additional file 1: Figure S1F), we showed the expression profiles of TD cells cluster together with those of doxorubicin-resistant cell lines, whereas the parental line A549 with the sensitive cells (Fig. 1b). As predicted, high resistance of TD cells to other conventional chemoradiotherapies such as etoposide (ETO), camptothecin (CPT) (Fig. 1c), ionizing radiation (IR) (Fig. 1d) and doxorubicin (Fig. 1e-f), all of which trigger apoptosis by inducing DNA damage.
For prediction of genes for chemoresistance, we defined the frequency of a gene in the up-regulated DEGs as its ‘chemoresistance score’; this value was strongly correlated with the degree of overexpression in TD (Spearman correlation = 0.44, P < 0.0004, Fig. 2a). Two genes that satisfied all three of the criteria; i) genes frequently represented in the up DEGs of the 32 chemotherapeutic drugs; ii) genes overexpressed in TD cells; and iii) EMT-related genes were ITGB3 and CTGF (Fig. 2a). Of note, the role of CTGF in chemoresistance and its antagonism for chemosensitization has been determined . In particular, we found that the IC50 of doxorubicin is positively correlated with the ITGB3 expression level (Additional file 1: Figure S2A). Analyzing a genome-scaling RNAi screening for 501 cancer cells, Project Achilles , 105 cells (21%) showed significant dependency on ITGB3 (Additional file 1: Figure S2B), which was among the top ~ 10% vulnerable genes (Additional file 1: Figure S2C). Moreover, ITGB3-dependency was more significant in the resistant cells for most chemotherapeutic drugs (Fig. 2b and Additional file 1: Figure S2D). Such trend was most evident in eight chemotherapeutic drugs (Fisher’s method P < 1.1 × 10− 5, Additional file 1: Figure S2E). Consistently, loss of ITGB3 in TD cells (shβ#3 and shβ#4) (Additional file 1: Figure S2F) increased sensitivity to doxorubicin treatment (Fig. 2c and Additional file 1: Figure S2G–I). Similar results were obtained with CPT, IR, and ETO (Additional file 1: Figure S2J-L). Given that expression of ITGB3 was sufficient to restore the chemoresistance of shβ#3 TD cells (Additional file 1: Figure S2M) and even increase the chemoresistance of A549, the parental cell line of TD (Fig. 2d), we conclude that ITGB3 expression is solely sufficient to induce chemoresistance.
Among the down-regulated pathways by ITGB3 depletion, NF-κB was the signaling pathways most strongly associated with patients’ survival (Fig. 2e). Given that depletion of ITGB3 down-regulated NF-κB-dependent survival factors (IL8, XIAP, PLAU, BIRC2/3, BCL2, or BCL2L1), and induced negative feedback regulators such as NFKBIA and TNFAIP3 (encoding IκBα and A20 deubiquitinase, respectively), we hypothesized that inhibition of NF-κB signaling would be a key process required for cell sensitization (Fig. 2f and Additional file 1: Figure S3A). Consistently, highly chemoresistant TD cells exhibited higher levels of NF-κB activity than A549 cells, as determined by nuclear p65 localization (Fig. 2g and Additional file 1: Figure S3B), the protein level of IκBα (Additional file 1: Figure S3C) and NF-κB reporter activity (Additional file 1: Figure S3D). Furthermore, loss of ITGB3 markedly attenuated NF-κB activity, as determined by NF-κB reporter activity (Fig. 2h) and acetylation of p65 (which is critical for its DNA-binding affinity) (Fig. 2i) as well as the level of nuclear p65 (Additional file 1: Figure S3E–F). In adverse, ectopic expression of ITGB3 restored NF-κB activity, as determined by acetylation of p65 and level of IκB (Additional file 1: Figure S3G). Together, these data indicate that ITGB3 expression is closely associated with NF-κB activity. According to this prediction that elevation of NF-κB activity by ITGB3 expression could be a primary cause of the elevated hazard ratio (Fig. 2e), we assessed the cytotoxicity of doxorubicin following abrogation of NF-κB activity. Depletion of p65 with siRNA, was sufficient to sensitize cells to doxorubicin treatment (Fig. 2j), suggesting that the increase in NF-κB activity mediated by ITGB3 expression is responsible for acquisition of chemoresistance.
To predict chemosensitizing drug candidates, we leveraged drug-induced transcriptome data from the Connectivity Map (CMap) and searched for drugs with expression signatures similar to those of ITGB3 depletion or NF-κB signaling inhibition (Fig. 3a). Among the candidate drugs, atorvastatin (ATV) was the only drug identified based on both ITGB3 and NF-κB gene signatures (Fig. 3b and Additional file 1: Figure S4A-B). As predicted, pretreatment with ATV significantly sensitized TD cells to doxorubicin (Fig. 3c). It is noteworthy that, despite clear induction of ITGB3 by treatment of doxorubicin, combined treatment of ATV and doxorubicin increased the rate of cell death (Fig. 3d), decreased cell viability (Fig. 3e) and apoptotic cell death (Fig. 3f). The increase in chemosensitivity following ATV treatment occurred in parallel with a reduction in NF-κB reporter activity (Fig. 3g). As similar as induction of NFKBIA (encoding IκBα) and suppression of BCL2L1 (encoding BCL-xL) by ITGB3 depletion (Additional file 1: Figure S3A), which may account for the decrease in NF-κB activity as well as pro-survival activity, short-course ATV treatment increased expression of IκBα and decreased the level of BCL-xL (Fig. 3h), suggesting that attenuation of the NF-κB–dependent pro-survival pathway by ATV leads to chemosensitization. Consistently, H460 cancer cells with mesenchymal gene expression (Additional file 1: Figure S5A) and high ITGB3 and IL6 expression (Additional file 1: Figure S5B), which was repressed by loss of ITGB3 (Additional file 1: Figure S5C and D) became more chemosensitive by ATV treatment (Additional file 1: Figure S5F), lowing NF-κB activity (Additional file 1: Figure S5E). In other hand, H358 cancer cells with epithelial gene expression (Additional file 1: Figure S5G) were likely to acquire chemoresistance (Additional file 1: Figure S5H) and increased NF-κB activation (Additional file 1: Figure S5I) by ITGB3 ectopic expression, which were weakened by ATV treatment (Additional file 1: Figure S5H and I). Conversely, depletion of ITGB3 promoted chemosensitivity in H358 (Additional file 1: Figure S5J).
Most targets responsible for acquired chemoresistance in cancers, identified during extensive molecular mechanistic studies, remain undrugged  due to poor druggability or possible side effects by direct inhibition. Thus, we took advantage of CMap approach based on a large-scale drug-induced transcriptome dataset and identified ATV, one of the world’s best-selling drugs for hyperlipidemia, as a candidate drug for abrogating the pro-survival and chemoresistance effect of ITGB3; specifically, we showed that the transcriptional profile of ATV-treated cells was similar to that of ITGB3 knockdown. Consistently, ATV has a radiosensitizing effect on prostate cancer cells . Although the inhibitory effect of STATINs on NF-κB is varied markedly , our predictive analysis identified ATV as a top-ranking candidate, strongly validating our data-driven approach.
By integrating pharmacogenomics and chemical genomic data, we successfully identified both a therapeutic target and a novel chemosensitizing drug to overcome resistance to multiple chemotherapeutic drugs. Our approach can be applied to a wide range of targets beyond those associated with EMT, paving an alternative path to drug discovery even for undruggable targets.
Area under fitted curve
Cancer Cell Line Encyclopedia
Connective Tissue Growth Factor
Cancer Therapeutics Response Portal
Differentially expressed gene
Genomics Data Common
Genomics of Drug Sensitivity in Cancer
Integrin beta 3
Fischer KR, Durrans A, Lee S, Sheng J, Li F, Wong ST, Choi H, El Rayes T, Ryu S, Troeger J, et al. Epithelial-to-mesenchymal transition is not required for lung metastasis but contributes to chemoresistance. Nature. 2015;527:472–6.
De Sousa E, Melo F, Wang X, Jansen M, Fessler E, Trinh A, de Rooij LPMH, de Jong JH, de Boer OJ, van Leersum R, Bijlsma MF, et al. Poor-prognosis colon cancer is defined by a molecularly distinct subtype and develops from serrated precursor lesions. Nat Med. 2013;19:614–8.
Sistigu A, Yamazaki T, Vacchelli E, Chaba K, Enot DP, Adam J, Vitale I, Goubar A, Baracco EE, Remedios C, et al. Cancer cell-autonomous contribution of type I interferon signaling to the efficacy of chemotherapy. Nat Med. 2014;20:1301–9.
Hong SK, Park JR, Kwon OS, Kim KT, Bae GY, Cha HJ. Induction of integrin beta3 by sustained ERK activity promotes the invasiveness of TGFbeta-induced mesenchymal tumor cells. Cancer Lett. 2016;376:339–46.
Kwon O-S, Hong S-K, Kwon S-J, Go Y-H, Oh E, Cha H-J. BCL2 induced by LAMTOR3/MAPK is a druggable target of chemoradioresistance in mesenchymal lung cancer. Cancer Lett. 2017;403:48–58.
Neesse A, Frese KK, Bapiro TE, Nakagawa T, Sternlicht MD, Seeley TW, Pilarsky C, Jodrell DI, Spong SM, Tuveson DA. CTGF antagonism with mAb FG-3019 enhances chemotherapy response without increasing drug delivery in murine ductal pancreas cancer. Proc Natl Acad Sci U S A. 2013;110:12325–30.
Tsherniak A, Vazquez F, Montgomery PG, Weir BA, Kryukov G, Cowley GS, Gill S, Harrington WF, Pantel S, Krill-Burger JM, et al. Defining a cancer dependency map. Cell. 2017;170:564–76 e516.
Dang CV, Reddy EP, Shokat KM, Soucek L. Drugging the ‘undruggable’ cancer targets. Nat Rev Cancer. 2017;17:502–8.
He M, Hu J, Xia Y. Large scale expressed sequence tag (EST) analysis of Metarhizium acridum infecting Locusta migratoria reveals multiple strategies for fungal adaptation to the host cuticle. Curr Genet. 2012;58:265–79.
Hilgendorff A, Muth H, Parviz B, Staubitz A, Haberbosch W, Tillmanns H, Holschermann H. Statins differ in their ability to block NF-kappaB activation in human blood monocytes. Int J Clin Pharmacol Ther. 2003;41:397–401.
We appreciate Seon-Young Kim and Dong-Uk Kim at Korea Research Institute of Bioscience and Biotechnology (KRIBB) for helpful discussions.
This work was supported by a grant from the National Research Foundation of Korea (NRF-2017M3C9A5028691 from HJC, NRF-2017M3C9A5028690 from WKK and NRF-2017R1A6A3A11030794 from HSL).
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Hong, SK., Lee, H., Kwon, OS. et al. Large-scale pharmacogenomics based drug discovery for ITGB3 dependent chemoresistance in mesenchymal lung cancer. Mol Cancer 17, 175 (2018). https://doi.org/10.1186/s12943-018-0924-8