Systematic analysis of overall survival and interactions between tumor mutations and drug treatment
Journal of Hematology & Oncology volume 9, Article number: 15 (2016)
Few exceptional responses in cancer treatment were attributed to a genetic predisposition of the tumor.
We analyzed a cohort of 3105 patients from 12 different cancer types and systematically sought the existence of a correlation between overall survival and the interaction of 21 antineoplastic treatments with 6 tumor mutations.
We identified a single significant correlation resulting in increased overall survival from temozolomide in lower-grade glioma with IDH1 R132H mutations. The trend could not be attributed to either the treatment or the mutation alone. Univariate and multivariate Cox regression demonstrated that this interaction stood as an independent prognostic predictor of survival.
Our results suggest infrequent instances of exceptional responses ascribable to tumor genomics yet corroborate the existence of an interaction of temozolomide with IDH1 mutations in lower-grade glioma.
The cancer genome can elicit sensitivity to certain drugs not specifically designed to target the underlying genetic aberrations. To this end, genomic markers of drug sensitivity have been systematically assessed in cancer cell lines [1, 2]. Ideally, these markers can identify patients who may better benefit from a certain antineoplastic drug [3, 4]. In contrast to the increasing availability of data about genomics of drug sensitivity in vitro , the association with improved patient survival is so far limited to few clinical cases, e.g., exceptional responses to everolimus in bladder cancers with TSC1 mutations .
Here, we sought to systematically assess if the chances of overall survival in patients with a certain cancer type and treated with a given antineoplastic drug correlate with the presence of a certain genetic mutation in the tumor. The examined cohort comprised 3105 patients, spanning 12 cancer types (with 81–731 samples for each cancer type). Collectively, 21 antineoplastic drugs were administered each in at least 20 patients (median 82; IQR 29–150). Six cancer-associated mutations were detected in at least 20 patients in this cohort: V600E in BRAF (n = 29), R132H in IDH1 (n = 108), G12V in KRAS (n = 49), H1047R in PIK3CA (n = 89), R175H in TP53 (n = 45), and V777 deletion in ZFHX3 (n = 22). After binning samples by cancer type, out of 1512 potential associations, 9 associations between overall survival, drug treatment, and tumor mutation had sufficient sample size for each covariate and were hereby tested. The hazard ratio (HR) for each interaction between drug treatment and tumor mutation in a cancer type was estimated in a multivariate analysis using a nested Cox proportional hazard regression model. We adopted a likelihood ratio test to test whether there is a significant effect of the interaction on overall survival on top of the tumor mutation and administered drug alone (Additional file 1: Table S1).
We observed a significant effect only in one scenario, the interaction between temozolomide (TMZ) and R132H mutations in IDH1 on the overall survival of lower-grade glioma (LGG) (likelihood ratio test p = 0.026). This test suggests that the correlation with survival is specific to the interaction between TMZ and R132H mutations in IDH1 and not associated with the drug treatment or the mutation per se, as demonstrated by the Kaplan-Meier curves generated for patients stratified upon these features (log-rank test p = 0.047, Fig. 1). The median overall survival for patients with the interaction was 95 months (95 % CI, 63—N.E.) and for patients without the interaction was 62 months (95 % CI, 49–87).
We detected a significant prognostic value for the interaction using a univariate Cox proportional hazard regression model (p = 0.016, Table 1). However, the interaction violated the proportional hazard assumption and showed a time-dependent effect. The univariate analysis was also run on validated prognostic factors in LGG  and additional clinical features (Table 1). A multivariate analysis based on significant factors from the univariate analysis and including a time-dependent effect for the interaction revealed an independent positive correlation between the interaction and overall survival (HR 0.09, 95 % CI 0.01–0.58, p = 0.012), which tends to diminish over time (Additional file 1: Figure S1).
In conclusion, we identified one genomic marker of drug sensitivity that was associated with better survival in patients, in contrast to patients treated with the same drug but with no detected mutation or vice versa. Indeed, mutations in IDH1 were previously implicated with good prognosis in brain tumors treated with TMZ [8, 9]. Our results independently validate these findings and further extend the reach of this correlation beyond some previous limitations . First and foremost, the cohort size allowed discerning that an increase in patient survival was exquisitely associated with the interaction between IDH1 mutations and TMZ, suggestive of a synergy between treatment and tumor genomics. Second, it specifically correlated with R132H mutations. Finally, we recovered a negative time-dependent effect of the interaction, which is reminiscent of emergence of drug resistance and in line with the genetic evolution of lower-grade glioma attributed to TMZ treatment .
Barretina J, Caponigro G, Stransky N, Venkatesan K, Margolin AA, Kim S, et al. The Cancer Cell Line Encyclopedia enables predictive modelling of anticancer drug sensitivity (vol 483, pg 603, 2012). Nature. 2012;492:290.
Garnett MJ, Edelman EJ, Heidorn SJ, Greenman CD, Dastur A, Lau KW, et al. (2012). Systematic identification of genomic markers of drug sensitivity in cancer cells. Nature 483, 570-U587.X
Garraway LA, Lander ES. Lessons from the cancer genome. Cell. 2013;153:17–37.
Stratton M, Garnett M, Edelman EJ, Heidorn S, Futreal PA, Haber D, et al. The genomics of drug sensitivity in cancer. European Journal of Cancer. 2012;48:S8.
Yang WJ, Soares J, Greninger P, Edelman EJ, Lightfoot H, Forbes S, et al. Genomics of drug sensitivity in cancer (GDSC): a resource for therapeutic biomarker discovery in cancer cells. Nucleic acids research. 2013;41:D955–61.
Iyer G, Hanrahan AJ, Milowsky MI, Al-Ahmadie H, Scott SN, Janakiraman M, et al. Genome sequencing identifies a basis for everolimus sensitivity. Science. 2012;338:221.
Pignatti F, van den Bent M, Curran D, Debruyne C, Sylvester R, Therasse P, et al. Prognostic factors for survival in adult patients with cerebral low-grade glioma. Journal of clinical oncology : official journal of the American Society of Clinical Oncology. 2002;20:2076–84.
Houillier C, Wang X, Kaloshi G, Mokhtari K, Guillevin R, Laffaire J, et al. IDH1 or IDH2 mutations predict longer survival and response to temozolomide in low-grade gliomas. Neurology. 2010;75:1560–6.
Kong DS, Kim HR, Choi YR, Seol HJ, Lee JI, Nam DH. Prognostic impact of molecular phenotype in patients with recurrent anaplastic glioma treated with prolonged administration of temozolomide. Journal of clinical neuroscience : official journal of the Neurosurgical Society of Australasia. 2015;22:1425–9.
Dubbink HJ, Taal W, van Marion R, Kros JM, van Heuvel I, Bromberg JE, et al. IDH1 mutations in low-grade astrocytomas predict survival but not response to temozolomide. Neurology. 2009;73:1792–5.
Johnson BE, Mazor T, Hong C, Barnes M, Aihara K, McLean CY, et al. Mutational analysis reveals the origin and therapy-driven evolution of recurrent glioma. Science. 2014;343:189–93.
The authors acknowledge Knut and Alice Wallenberg Foundation for financing this work.
The authors declare no conflict of interest in relation to the submitted work.
FG conceived and designed the study and performed the analyses; JN supervised the study. All authors approved the manuscript in its final form.
About this article
Cite this article
Gatto, F., Nielsen, J. Systematic analysis of overall survival and interactions between tumor mutations and drug treatment. J Hematol Oncol 9, 15 (2016). https://doi.org/10.1186/s13045-016-0249-2