Systematic Analysis of Genetic and Pathway Determinants of Eribulin Sensitivity across 100 Human Cancer Cell Lines from the Cancer Cell Line Encyclopedia (CCLE)

SIMPLE SUMMARY: Despite decades of clinical use and detailed understandings of their mechanisms of action, clinically useful predictive biomarkers for approved microtubule targeting agents such as eribulin, paclitaxel and vinorelbine have eluded development. Our results now provide the basis for gene expression-based, drug-specific predictive biomarkers for eribulin and vinorelbine, as well as deeper understandings of the molecular pathways associated with eribulin and vinorelbine response. ABSTRACT: Eribulin, a natural product-based microtubule targeting agent with cytotoxic and noncytotoxic mechanisms, is FDA approved for certain patients with advanced breast cancer and liposarcoma. To investigate the feasibility of developing drug-specific predictive biomarkers, we quantified antiproliferative activities of eribulin versus paclitaxel and vinorelbine against 100 human cancer cell lines from the Cancer Cell Line Encyclopedia, and correlated results with publicly available databases to identify genes and pathways associated with eribulin response, either uniquely or shared with paclitaxel or vinorelbine. Mean expression ratios of 11,985 genes between the most and least sensitive cell line quartiles were sorted by p-values and drug overlaps, yielding 52, 29 and 80 genes uniquely associated with eribulin, paclitaxel and vinorelbine, respectively. Further restriction to minimum 2-fold ratios followed by reintroducing data from the middle two quartiles identified 9 and 13 drug-specific unique fingerprint genes for eribulin and vinorelbine, respectively; surprisingly, no gene met all criteria for paclitaxel. Interactome and Reactome pathway analyses showed that unique fingerprint genes of both drugs were primarily associated with cellular signaling, not microtubule-related pathways, although considerable differences existed in individual pathways identified. Finally, four-gene (C5ORF38, DAAM1, IRX2, CD70) and five-gene (EPHA2, NGEF, SEPTIN10, TRIP10, VSIG10) multivariate regression models for eribulin and vinorelbine showed high statistical correlation with drug-specific responses across the 100 cell lines and accurately calculated predicted mean IC50s for the most and least sensitive cell line quartiles as surrogates for responders and nonresponders, respectively. Collectively, these results provide a foundation for developing drug-specific predictive biomarkers for eribulin and vinorelbine.