PharmacoGx: an R package for analysis of large pharmacogenomic datasets

A review of connectivity map and computational approaches in pharmacogenomics

Large-scale perturbation databases, such as Connectivity Map (CMap) or Library of Integrated Network-based Cellular Signatures (LINCS), provide enormous opportunities for computational pharmacogenomics and drug design. A reason for this is that in contrast to classical pharmacology focusing at one target at a time, the transcriptomics profiles provided by CMap and LINCS open the door for systems biology approaches on the pathway and network level. In this article, we provide a review of recent developments in computational pharmacogenomics with respect to CMap and LINCS and related applications.

Pharmacogenomics and big genomic data: from lab to clinic and back again

The field of pharmacogenomics is an area of great potential for near-term human health impacts from the big genomic data revolution. Pharmacogenomics research momentum is building with numerous hypotheses currently being investigated through the integration of molecular profiles of different cell lines and large genomic data sets containing information on cellular and human responses to therapies. Additionally, the results of previous pharmacogenetic research efforts have been formulated into clinical guidelines that are beginning to impact how healthcare is conducted on the level of the individual patient. This trend will only continue with the recent release of new datasets containing linked genotype and electronic medical record data. This review discusses key resources available for pharmacogenomics and pharmacogenetics research and highlights recent work within the field.

Personalised drug repositioning for Clear Cell Renal Cell Carcinoma using gene expression

Reversal of cancer gene expression is predictive of therapeutic potential and can be used to find new indications for existing drugs (drug repositioning). Gene expression reversal potential is currently calculated, in almost all studies, by pre-aggregating all tumour samples into a single group signature or a limited number of molecular subtype signatures. Here, we investigate whether drug repositioning based on individual tumour sample gene expression signatures outperforms the use of tumour group and subtype signatures. The tumour signatures were created using 534 tumour samples and 72 matched normal samples from 530 clear cell renal cell carcinoma (ccRCC) patients. More than 20,000 drug signatures were extracted from the CMAP and LINCS databases. We show that negative enrichment of individual tumour samples correlated (Spearman's rho = 0.15) much better with the amount of differentially expressed genes in drug signatures than with the tumour group signature (Rho = 0.08) and the 4 tumour subtype signatures (Rho 0.036-0.11). Targeted drugs used against ccRCC, such as sirolimus and temsirolimus, which could not be identified with the pre-aggregated tumour signatures could be recovered using individual sample analysis. Thus, drug repositioning can be personalized by taking into account the gene expression profile of the individual's tumour sample.

Disruption of the anaphase-promoting complex confers resistance to TTK inhibitors in triple-negative breast cancer

TTK protein kinase (TTK), also known as Monopolar spindle 1 (MPS1), is a key regulator of the spindle assembly checkpoint (SAC), which functions to maintain genomic integrity. TTK has emerged as a promising therapeutic target in human cancers, including triple-negative breast cancer (TNBC). Several TTK inhibitors (TTKis) are being evaluated in clinical trials, and an understanding of the mechanisms mediating TTKi sensitivity and resistance could inform the successful development of this class of agents. We evaluated the cellular effects of the potent clinical TTKi CFI-402257 in TNBC models. CFI-402257 induced apoptosis and potentiated aneuploidy in TNBC lines by accelerating progression through mitosis and inducing mitotic segregation errors. We used genome-wide CRISPR/Cas9 screens in multiple TNBC cell lines to identify mechanisms of resistance to CFI-402257. Our functional genomic screens identified members of the anaphase-promoting complex/cyclosome (APC/C) complex, which promotes mitotic progression following inactivation of the SAC. Several screen candidates were validated to confer resistance to CFI-402257 and other TTKis using CRISPR/Cas9 and siRNA methods. These findings extend the observation that impairment of the APC/C enables cells to tolerate genomic instability caused by SAC inactivation, and support the notion that a measure of APC/C function could predict the response to TTK inhibition. Indeed, an APC/C gene expression signature is significantly associated with CFI-402257 response in breast and lung adenocarcinoma cell line panels. This expression signature, along with somatic alterations in genes involved in mitotic progression, represent potential biomarkers that could be evaluated in ongoing clinical trials of CFI-402257 or other TTKis.

Statin-Induced Cancer Cell Death Can Be Mechanistically Uncoupled from Prenylation of RAS Family Proteins

The statin family of drugs preferentially triggers tumor cell apoptosis by depleting mevalonate pathway metabolites farnesyl pyrophosphate (FPP) and geranylgeranyl pyrophosphate (GGPP), which are used for protein prenylation, including the oncoproteins of the RAS superfamily. However, accumulating data indicate that activation of the RAS superfamily are poor biomarkers of statin sensitivity, and the mechanism of statin-induced tumor-specific apoptosis remains unclear. Here we demonstrate that cancer cell death triggered by statins can be uncoupled from prenylation of the RAS superfamily of oncoproteins. Ectopic expression of different members of the RAS superfamily did not uniformly sensitize cells to fluvastatin, indicating that increased cellular demand for protein prenylation cannot explain increased statin sensitivity. Although ectopic expression of HRAS increased statin sensitivity, expression of myristoylated HRAS did not rescue this effect. HRAS-induced epithelial-to-mesenchymal transition (EMT) through activation of zinc finger E-box binding homeobox 1 (ZEB1) sensitized tumor cells to the antiproliferative activity of statins, and induction of EMT by ZEB1 was sufficient to phenocopy the increase in fluvastatin sensitivity; knocking out ZEB1 reversed this effect. Publicly available gene expression and statin sensitivity data indicated that enrichment of EMT features was associated with increased sensitivity to statins in a large panel of cancer cell lines across multiple cancer types. These results indicate that the anticancer effect of statins is independent from prenylation of RAS family proteins and is associated with a cancer cell EMT phenotype.Significance: The use of statins to target cancer cell EMT may be useful as a therapy to block cancer progression. Cancer Res; 78(5); 1347-57. ©2017 AACR.

xCell: digitally portraying the tissue cellular heterogeneity landscape

Tissues are complex milieus consisting of numerous cell types. Several recent methods have attempted to enumerate cell subsets from transcriptomes. However, the available methods have used limited sources for training and give only a partial portrayal of the full cellular landscape. Here we present xCell, a novel gene signature-based method, and use it to infer 64 immune and stromal cell types. We harmonized 1822 pure human cell type transcriptomes from various sources and employed a curve fitting approach for linear comparison of cell types and introduced a novel spillover compensation technique for separating them. Using extensive in silico analyses and comparison to cytometry immunophenotyping, we show that xCell outperforms other methods. xCell is available at http://xCell.ucsf.edu/.

Gene isoforms as expression-based biomarkers predictive of drug response in vitro

Next-generation sequencing technologies have recently been used in pharmacogenomic studies to characterize large panels of cancer cell lines at the genomic and transcriptomic levels. Among these technologies, RNA-sequencing enable profiling of alternatively spliced transcripts. Given the high frequency of mRNA splicing in cancers, linking this feature to drug response will open new avenues of research in biomarker discovery. To identify robust transcriptomic biomarkers for drug response across studies, we develop a meta-analytical framework combining the pharmacological data from two large-scale drug screening datasets. We use an independent pan-cancer pharmacogenomic dataset to test the robustness of our candidate biomarkers across multiple cancer types. We further analyze two independent breast cancer datasets and find that specific isoforms of IGF2BP2, NECTIN4, ITGB6, and KLHDC9 are significantly associated with AZD6244, lapatinib, erlotinib, and paclitaxel, respectively. Our results support isoform expressions as a rich resource for biomarkers predictive of drug response.

GRcalculator: an online tool for calculating and mining dose-response data

BACKGROUND

Quantifying the response of cell lines to drugs or other perturbagens is the cornerstone of pre-clinical drug development and pharmacogenomics as well as a means to study factors that contribute to sensitivity and resistance. In dividing cells, traditional metrics derived from dose-response curves such as IC 50 , AUC, and E max , are confounded by the number of cell divisions taking place during the assay, which varies widely for biological and experimental reasons. Hafner et al. (Nat Meth 13:521-627, 2016) recently proposed an alternative way to quantify drug response, normalized growth rate (GR) inhibition, that is robust to such confounders. Adoption of the GR method is expected to improve the reproducibility of dose-response assays and the reliability of pharmacogenomic associations (Hafner et al. 500-502, 2017).

RESULTS

We describe here an interactive website ( www.grcalculator.org ) for calculation, analysis, and visualization of dose-response data using the GR approach and for comparison of GR and traditional metrics. Data can be user-supplied or derived from published datasets. The web tools are implemented in the form of three integrated Shiny applications (grcalculator, grbrowser, and grtutorial) deployed through a Shiny server. Intuitive graphical user interfaces (GUIs) allow for interactive analysis and visualization of data. The Shiny applications make use of two R packages (shinyLi and GRmetrics) specifically developed for this purpose. The GRmetrics R package is also available via Bioconductor and can be used for offline data analysis and visualization. Source code for the Shiny applications and associated packages (shinyLi and GRmetrics) can be accessed at www.github.com/uc-bd2k/grcalculator and www.github.com/datarail/gr_metrics .

CONCLUSIONS

GRcalculator is a powerful, user-friendly, and free tool to facilitate analysis of dose-response data. It generates publication-ready figures and provides a unified platform for investigators to analyze dose-response data across diverse cell types and perturbagens (including drugs, biological ligands, RNAi, etc.). GRcalculator also provides access to data collected by the NIH LINCS Program ( http://www.lincsproject.org /) and other public domain datasets. The GRmetrics Bioconductor package provides computationally trained users with a platform for offline analysis of dose-response data and facilitates inclusion of GR metrics calculations within existing R analysis pipelines. These tools are therefore well suited to users in academia as well as industry.

Integrative Cancer Pharmacogenomics to Infer Large-Scale Drug Taxonomy

Identification of drug targets and mechanism of action (MoA) for new and uncharacterized anticancer drugs is important for optimization of treatment efficacy. Current MoA prediction largely relies on prior information including side effects, therapeutic indication, and chemoinformatics. Such information is not transferable or applicable for newly identified, previously uncharacterized small molecules. Therefore, a shift in the paradigm of MoA predictions is necessary toward development of unbiased approaches that can elucidate drug relationships and efficiently classify new compounds with basic input data. We propose here a new integrative computational pharmacogenomic approach, referred to as Drug Network Fusion (DNF), to infer scalable drug taxonomies that rely only on basic drug characteristics toward elucidating drug-drug relationships. DNF is the first framework to integrate drug structural information, high-throughput drug perturbation, and drug sensitivity profiles, enabling drug classification of new experimental compounds with minimal prior information. DNF taxonomy succeeded in identifying pertinent and novel drug-drug relationships, making it suitable for investigating experimental drugs with potential new targets or MoA. The scalability of DNF facilitated identification of key drug relationships across different drug categories, providing a flexible tool for potential clinical applications in precision medicine. Our results support DNF as a valuable resource to the cancer research community by providing new hypotheses on compound MoA and potential insights for drug repurposing. Cancer Res; 77(11); 3057-69. ©2017 AACR.

DNA replication stress: a source of APOBEC3B expression in breast cancer

APOBEC cytidine deaminases have been implicated as major contributors to the mutation burden in many cancers on the basis of their mutational signature. A new experimental study sheds light on the inciting factors, linking APOBEC3B expression to oncogene- and drug-induced replication stress.

Revisiting inconsistency in large pharmacogenomic studies

In 2013, we published a comparative analysis of mutation and gene expression profiles and drug sensitivity measurements for 15 drugs characterized in the 471 cancer cell lines screened in the Genomics of Drug Sensitivity in Cancer (GDSC) and Cancer Cell Line Encyclopedia (CCLE). While we found good concordance in gene expression profiles, there was substantial inconsistency in the drug responses reported by the GDSC and CCLE projects. We received extensive feedback on the comparisons that we performed. This feedback, along with the release of new data, prompted us to revisit our initial analysis. We present a new analysis using these expanded data, where we address the most significant suggestions for improvements on our published analysis - that targeted therapies and broad cytotoxic drugs should have been treated differently in assessing consistency, that consistency of both molecular profiles and drug sensitivity measurements should be compared across cell lines, and that the software analysis tools provided should have been easier to run, particularly as the GDSC and CCLE released additional data. Our re-analysis supports our previous finding that gene expression data are significantly more consistent than drug sensitivity measurements. Using new statistics to assess data consistency allowed identification of two broad effect drugs and three targeted drugs with moderate to good consistency in drug sensitivity data between GDSC and CCLE. For three other targeted drugs, there were not enough sensitive cell lines to assess the consistency of the pharmacological profiles. We found evidence of inconsistencies in pharmacological phenotypes for the remaining eight drugs. Overall, our findings suggest that the drug sensitivity data in GDSC and CCLE continue to present challenges for robust biomarker discovery. This re-analysis provides additional support for the argument that experimental standardization and validation of pharmacogenomic response will be necessary to advance the broad use of large pharmacogenomic screens.

Revisiting inconsistency in large pharmacogenomic studies

In 2013, we published a comparative analysis of mutation and gene expression profiles and drug sensitivity measurements for 15 drugs characterized in the 471 cancer cell lines screened in the Genomics of Drug Sensitivity in Cancer (GDSC) and Cancer Cell Line Encyclopedia (CCLE). While we found good concordance in gene expression profiles, there was substantial inconsistency in the drug responses reported by the GDSC and CCLE projects. We received extensive feedback on the comparisons that we performed. This feedback, along with the release of new data, prompted us to revisit our initial analysis. We present a new analysis using these expanded data, where we address the most significant suggestions for improvements on our published analysis - that targeted therapies and broad cytotoxic drugs should have been treated differently in assessing consistency, that consistency of both molecular profiles and drug sensitivity measurements should be compared across cell lines, and that the software analysis tools provided should have been easier to run, particularly as the GDSC and CCLE released additional data. Our re-analysis supports our previous finding that gene expression data are significantly more consistent than drug sensitivity measurements. Using new statistics to assess data consistency allowed identification of two broad effect drugs and three targeted drugs with moderate to good consistency in drug sensitivity data between GDSC and CCLE. For three other targeted drugs, there were not enough sensitive cell lines to assess the consistency of the pharmacological profiles. We found evidence of inconsistencies in pharmacological phenotypes for the remaining eight drugs. Overall, our findings suggest that the drug sensitivity data in GDSC and CCLE continue to present challenges for robust biomarker discovery. This re-analysis provides additional support for the argument that experimental standardization and validation of pharmacogenomic response will be necessary to advance the broad use of large pharmacogenomic screens.

Revisiting inconsistency in large pharmacogenomic studies

In 2013, we published a comparative analysis mutation and gene expression profiles and drug sensitivity measurements for 15 drugs characterized in the 471 cancer cell lines screened in the Genomics of Drug Sensitivity in Cancer (GDSC) and Cancer Cell Line Encyclopedia (CCLE). While we found good concordance in gene expression profiles, there was substantial inconsistency in the drug responses reported by the GDSC and CCLE projects. We received extensive feedback on the comparisons that we performed. This feedback, along with the release of new data, prompted us to revisit our initial analysis. Here we present a new analysis using these expanded data in which we address the most significant suggestions for improvements on our published analysis - that targeted therapies and broad cytotoxic drugs should have been treated differently in assessing consistency, that consistency of both molecular profiles and drug sensitivity measurements should both be compared across cell lines, and that the software analysis tools we provided should have been easier to run, particularly as the GDSC and CCLE released additional data.             Our re-analysis supports our previous finding that gene expression data are significantly more consistent than drug sensitivity measurements. The use of new statistics to assess data consistency allowed us to identify two broad effect drugs and three targeted drugs with moderate to good consistency in drug sensitivity data between GDSC and CCLE. For three other targeted drugs, there were not enough sensitive cell lines to assess the consistency of the pharmacological profiles. We found evidence of inconsistencies in pharmacological phenotypes for the remaining eight drugs.             Overall, our findings suggest that the drug sensitivity data in GDSC and CCLE continue to present challenges for robust biomarker discovery. This re-analysis provides additional support for the argument that experimental standardization and validation of pharmacogenomic response will be necessary to advance the broad use of large pharmacogenomic screens.

Integrating heterogeneous drug sensitivity data from cancer pharmacogenomic studies

The consistency of in vitro drug sensitivity data is of key importance for cancer pharmacogenomics. Previous attempts to correlate drug sensitivities from the large pharmacogenomics databases, such as the Cancer Cell Line Encyclopedia (CCLE) and the Genomics of Drug Sensitivity in Cancer (GDSC), have produced discordant results. We developed a new drug sensitivity metric, the area under the dose response curve adjusted for the range of tested drug concentrations, which allows integration of heterogeneous drug sensitivity data from the CCLE, the GDSC, and the Cancer Therapeutics Response Portal (CTRP). We show that there is moderate to good agreement of drug sensitivity data for many targeted therapies, particularly kinase inhibitors. The results of this largest cancer cell line drug sensitivity data analysis to date are accessible through the online portal, which serves as a platform for high power pharmacogenomics analysis.

Assessment of pharmacogenomic agreement

In 2013 we published an analysis demonstrating that drug response data and gene-drug associations reported in two independent large-scale pharmacogenomic screens, Genomics of Drug Sensitivity in Cancer (GDSC) and Cancer Cell Line Encyclopedia (CCLE), were inconsistent. The GDSC and CCLE investigators recently reported that their respective studies exhibit reasonable agreement and yield similar molecular predictors of drug response, seemingly contradicting our previous findings. Reanalyzing the authors' published methods and results, we found that their analysis failed to account for variability in the genomic data and more importantly compared different drug sensitivity measures from each study, which substantially deviate from our more stringent consistency assessment. Our comparison of the most updated genomic and pharmacological data from the GDSC and CCLE confirms our published findings that the measures of drug response reported by these two groups are not consistent. We believe that a principled approach to assess the reproducibility of drug sensitivity predictors is necessary before envisioning their translation into clinical settings.