Clinical translation of cell-based pharmacogenomic discovery

High-throughput screening and genome-wide analyses of 44 anticancer drugs in the 1000 Genomes cell lines reveals an association of the NQO1 gene with the response of multiple anticancer drugs

Cancer patients exhibit a broad range of inter-individual variability in response and toxicity to widely used anticancer drugs, and genetic variation is a major contributor to this variability. To identify new genes that influence the response of 44 FDA-approved anticancer drug treatments widely used to treat various types of cancer, we conducted high-throughput screening and genome-wide association mapping using 680 lymphoblastoid cell lines from the 1000 Genomes Project. The drug treatments considered in this study represent nine drug classes widely used in the treatment of cancer in addition to the paclitaxel + epirubicin combination therapy commonly used for breast cancer patients. Our genome-wide association study (GWAS) found several significant and suggestive associations. We prioritized consistent associations for functional follow-up using gene-expression analyses. The NAD(P)H quinone dehydrogenase 1 (NQO1) gene was found to be associated with the dose-response of arsenic trioxide, erlotinib, trametinib, and a combination treatment of paclitaxel + epirubicin. NQO1 has previously been shown as a biomarker of epirubicin response, but our results reveal novel associations with these additional treatments. Baseline gene expression of NQO1 was positively correlated with response for 43 of the 44 treatments surveyed. By interrogating the functional mechanisms of this association, the results demonstrate differences in both baseline and drug-exposed induction.

In the burgeoning field of personalized medicine, genetic variation is recognized as a major contributor to patients’ differential responses to drugs. Lymphoblastoid cell lines (LCLs) are a consistent and convenient representation of cells used for in vitro research. Human genome sequencing with LCLs can identify new genes that influence individuals’ drug responses, including the dose-response relationship, which describes the relationship between physiological response and the amount of exposure to a substance. In this work, we conduct high-throughput screening and genome-wide association mapping using 680 LCLs from the 1000 Genomes Project to identify new genes that influence individual response to 44 widely used anticancer drugs. We found the NQO1 gene to be associated with the dose-response of several drugs, namely arsenic trioxide, erlotinib, trametinib, and the paclitaxel + epirubicin combination, and performed follow-up analyses to better understand its functional role in drug response. Our results indicate NQO1 expression is correlated with increased drug resistance and provide some evidence that SNP rs1800566 influences drug response by altering protein activity for these four treatments. With further research, NQO1 has potential use as a therapeutic target, for example, suppressing NQO1 expression to increase sensitivity to particular drugs.

Genetic Variants Contributing to Colistin Cytotoxicity: Identification of TGIF1 and HOXD10 Using a Population Genomics Approach

Colistin sulfate (polymixin E) is an antibiotic prescribed with increasing frequency for severe Gram-negative bacterial infections. As nephrotoxicity is a common side effect, the discovery of pharmacogenomic markers associated with toxicity would benefit the utility of this drug. Our objective was to identify genetic markers of colistin cytotoxicity that were also associated with expression of key proteins using an unbiased, whole genome approach and further evaluate the functional significance in renal cell lines. To this end, we employed International HapMap lymphoblastoid cell lines (LCLs) of Yoruban ancestry with known genetic information to perform a genome-wide association study (GWAS) with cellular sensitivity to colistin. Further association studies revealed that single nucleotide polymorphisms (SNPs) associated with gene expression and protein expression were significantly enriched in SNPs associated with cytotoxicity (p ≤ 0.001 for gene and p = 0.015 for protein expression). The most highly associated SNP, chr18:3417240 (p = 6.49 × 10⁻⁸), was nominally a cis-expression quantitative trait locus (eQTL) of the gene TGIF1 (transforming growth factor β (TGFβ)-induced factor-1; p = 0.021) and was associated with expression of the protein HOXD10 (homeobox protein D10; p = 7.17 × 10⁻⁵). To demonstrate functional relevance in a murine colistin nephrotoxicity model, HOXD10 immunohistochemistry revealed upregulated protein expression independent of mRNA expression in response to colistin administration. Knockdown of TGIF1 resulted in decreased protein expression of HOXD10 and increased resistance to colistin cytotoxicity. Furthermore, knockdown of HOXD10 in renal cells also resulted in increased resistance to colistin cytotoxicity, supporting the physiological relevance of the initial genomic associations.

van der Vaart A, Web BT, Bacanu SA, Kalsi G, Kendler KS, Miles MF, Dick D, Riley BP, Dumur C, Vladimirov VI. Integrating mRNA and miRNA Weighted Gene Co-Expression Networks with eQTLs in the Nucleus Accumbens of Subjects with Alcohol Dependence

Alcohol consumption is known to lead to gene expression changes in the brain. After performing weighted gene co-expression network analyses (WGCNA) on genome-wide mRNA and microRNA (miRNA) expression in Nucleus Accumbens (NAc) of subjects with alcohol dependence (AD; N = 18) and of matched controls (N = 18), six mRNA and three miRNA modules significantly correlated with AD were identified (Bonferoni-adj. p≤ 0.05). Cell-type-specific transcriptome analyses revealed two of the mRNA modules to be enriched for neuronal specific marker genes and downregulated in AD, whereas the remaining four mRNA modules were enriched for astrocyte and microglial specific marker genes and upregulated in AD. Gene set enrichment analysis demonstrated that neuronal specific modules were enriched for genes involved in oxidative phosphorylation, mitochondrial dysfunction and MAPK signaling. Glial-specific modules were predominantly enriched for genes involved in processes related to immune functions, i.e. cytokine signaling (all adj. p≤ 0.05). In mRNA and miRNA modules, 461 and 25 candidate hub genes were identified, respectively. In contrast to the expected biological functions of miRNAs, correlation analyses between mRNA and miRNA hub genes revealed a higher number of positive than negative correlations (χ² test p≤ 0.0001). Integration of hub gene expression with genome-wide genotypic data resulted in 591 mRNA cis-eQTLs and 62 miRNA cis-eQTLs. mRNA cis-eQTLs were significantly enriched for AD diagnosis and AD symptom counts (adj. p = 0.014 and p = 0.024, respectively) in AD GWAS signals in a large, independent genetic sample from the Collaborative Study on Genetics of Alcohol (COGA). In conclusion, our study identified putative gene network hubs coordinating mRNA and miRNA co-expression changes in the NAc of AD subjects, and our genetic (cis-eQTL) analysis provides novel insights into the etiological mechanisms of AD.

Combining genetic and nongenetic biomarkers to realize the promise of pharmacogenomics for inflammatory diseases

Many drugs used to treat inflammatory diseases are ineffective in a substantial proportion of patients. Identifying patients that are likely to respond to specific therapies would facilitate personalized treatment strategies that could improve outcomes while reducing costs and risks of adverse events. Despite these clear benefits, there are limited examples of predictive biomarkers of drug efficacy currently implemented into clinical practice for inflammatory diseases. We review efforts to identify genetic and nongenetic biomarkers of drug response in these diseases and consider potential benefits from combining multiple sources of biological data into multifeature predictive models.

Pharmacoethnicity in Paclitaxel-Induced Sensory Peripheral Neuropathy

PURPOSE

Paclitaxel is used worldwide in the treatment of breast, lung, ovarian, and other cancers. Sensory peripheral neuropathy is an associated adverse effect that cannot be predicted, prevented, or mitigated. To better understand the contribution of germline genetic variation to paclitaxel-induced peripheral neuropathy, we undertook an integrative approach that combines genome-wide association study (GWAS) data generated from HapMap lymphoblastoid cell lines (LCL) and Asian patients.

METHODS

GWAS was performed with paclitaxel-induced cytotoxicity generated in 363 LCLs and with paclitaxel-induced neuropathy from 145 Asian patients. A gene-based approach was used to identify overlapping genes and compare with a European clinical cohort of paclitaxel-induced neuropathy. Neurons derived from human-induced pluripotent stem cells were used for functional validation of candidate genes.

RESULTS

SNPs near AIPL1 were significantly associated with paclitaxel-induced cytotoxicity in Asian LCLs (P < 10(-6)). Decreased expression of AIPL1 resulted in decreased sensitivity of neurons to paclitaxel by inducing neurite morphologic changes as measured by increased relative total outgrowth, number of processes and mean process length. Using a gene-based analysis, there were 32 genes that overlapped between Asian LCL cytotoxicity and Asian patient neuropathy (P < 0.05), including BCR. Upon BCR knockdown, there was an increase in neuronal sensitivity to paclitaxel as measured by neurite morphologic characteristics.

CONCLUSIONS

We identified genetic variants associated with Asian paclitaxel-induced cytotoxicity and functionally validated the AIPL1 and BCR in a neuronal cell model. Furthermore, the integrative pharmacogenomics approach of LCL/patient GWAS may help prioritize target genes associated with chemotherapeutic-induced peripheral neuropathy.

Identification of genetic variants associated with capecitabine-induced hand-foot syndrome through integration of patient and cell line genomic analyses

OBJECTIVE

A primary challenge in identifying replicable pharmacogenomic markers from clinical genomewide association study (GWAS) trials in oncology is the difficulty in performing a second large clinical trial with the same drugs and dosage regimen. We sought to overcome this challenge by incorporating GWAS results from cell-based studies using the same chemotherapy as a clinical cohort.

METHODS

In this study, we test whether the overlap between genetic variants identified in a preclinical study and a clinical study on capecitabine is more than expected by chance. A GWAS of capecitabine-induced cytotoxicity was performed in 164 lymphoblastoid cell lines derived from the CEU HapMap population and compared with a GWAS of hand-foot syndrome (HFS), the most frequent capecitabine-induced adverse drug reaction, in Spanish breast and colorectal cancer patients (n=160) treated with capecitabine.

RESULTS

We observed an overlap of 16 single nucleotide polymorphisms associated with capecitabine-induced cytotoxicity (P<0.001) in lymphoblastoid cell lines and HFS (P<0.05) in patients, which is a greater overlap than expected by chance (genotype-phenotype permutation empirical P=0.015). Ten tag single nucleotide polymorphisms, which cover the overlap loci, were genotyped in a second patient cohort (n=85) and one of them, rs9936750, was associated with capecitabine-induced HFS (P=0.0076).

CONCLUSION

The enrichment results imply that cellular models of capecitabine-induced cytotoxicity may capture components of the underlying polygenic architecture of related toxicities in patients.

Institutional Profile: University of Chicago Center for Personalized Therapeutics: research, education and implementation science

Pharmacogenomics is aimed at advancing our knowledge of the genetic basis of variable drug response. The Center for Personalized Therapeutics within the University of Chicago comprises basic, translational and clinical research as well as education including undergraduate, graduate, medical students, clinical/postdoctoral fellows and faculty. The Committee on Clinical Pharmacology and Pharmacogenomics is the educational arm of the Center aimed at training clinical and postdoctoral fellows in translational pharmacology and pharmacogenomics. Research runs the gamut from basic discovery and functional studies to pharmacogenomic implementation studies to evaluate physician adoption of genetic medicine. The mission of the Center is to facilitate research, education and implementation of pharmacogenomics to realize the true potential of personalized medicine and improve the lives of patients.

Chemotherapy-induced peripheral neurotoxicity and ototoxicity: new paradigms for translational genomics

In view of advances in early detection and treatment, the 5-year relative survival rate for all cancer patients combined is now approximately 66%. As a result, there are more than 13.7 million cancer survivors in the United States, with this number increasing by 2% annually. For many patients, improvements in survival have been countered by therapy-associated adverse effects that may seriously impair long-term functional status, workplace productivity, and quality of life. Approximately 20% to 40% of cancer patients given neurotoxic chemotherapy develop chemotherapy-induced peripheral neurotoxicity (CIPN), which represents one of the most common and potentially permanent nonhematologic side effects of chemotherapy. Permanent bilateral hearing loss and/or tinnitus can result from several ototoxic therapies, including cisplatin- or carboplatin-based chemotherapy. CIPN and ototoxicity represent important challenges because of the lack of means for effective prevention, mitigation, or a priori identification of high-risk patients, and few studies have applied modern genomic approaches to understand underlying mechanisms/pathways. Translational genomics, including cell-based models, now offer opportunities to make inroads for the first time to develop preventive and interventional strategies for CIPN, ototoxicity, and other treatment-related complications. This commentary provides current perspective on a successful research strategy, with a focus on cisplatin, developed by an experienced, transdisciplinary group of researchers and clinicians, representing pharmacogenomics, statistical genetics, neurology, hearing science, medical oncology, epidemiology, and cancer survivorship. Principles outlined herein are applicable to the construction of research programs in translational genomics with strong clinical relevance and highlight unprecedented opportunities to understand, prevent, and treat long-term treatment-related morbidities.

Integrating cell-based and clinical genome-wide studies to identify genetic variants contributing to treatment failure in neuroblastoma patients

High-risk neuroblastoma is an aggressive malignancy with high rates of treatment failure. We evaluated genetic variants associated with in vitro sensitivity to two derivatives of cyclophosphamide for association with clinical response in a separate replication cohort of neuroblastoma patients (n=2,709). Lymphoblastoid cell lines (LCLs) were exposed to increasing concentrations of 4-hydroperoxycyclophosphamide [4HC n=422] and phosphoramide mustard [PM n=428] to determine sensitivity. Genome-wide association studies (GWAS) were performed to identify single nucleotide polymorphisms (SNPs) associated with 4HC and PM sensitivity. SNPs consistently associated with LCL sensitivity were analyzed for associations with event-free survival in patients. Two linked SNPs, rs9908694 and rs1453560, were found to be associated with PM sensitivity in LCLs across populations and were associated with event-free survival in all patients (P=0.01) and within the high-risk subset (P=0.05). Our study highlights the value of cell-based models to identify candidate variants that may predict response to treatment in patients with cancer.

Lymphoblastoid cell lines models of drug response: successes and lessons from this pharmacogenomic model

A new standard for medicine is emerging that aims to improve individual drug responses through studying associations with genetic variations. This field, pharmacogenomics, is undergoing a rapid expansion due to a variety of technological advancements that are enabling higher throughput with reductions in cost. Here we review the advantages, limitations, and opportunities for using lymphoblastoid cell lines (LCL) as a model system for human pharmacogenomic studies. There are a wide range of publicly available resources with genome-wide data available for LCLs from both related and unrelated populations, removing the cost of genotyping the data for drug response studies. Furthermore, in contrast to human clinical trials or in vivo model systems, with high-throughput in vitro screening technologies, pharmacogenomics studies can easily be scaled to accommodate large sample sizes. An important component to leveraging genome-wide data in LCL models is association mapping. Several methods are discussed herein, and include multivariate concentration response modeling, issues with multiple testing, and successful examples of the 'triangle model' to identify candidate variants. Once candidate gene variants have been determined, their biological roles can be elucidated using pathway analyses and functionally confirmed using siRNA knockdown experiments. The wealth of genomics data being produced using related and unrelated populations is creating many exciting opportunities leading to new insights into the genetic contribution and heritability of drug response.

EPS8 inhibition increases cisplatin sensitivity in lung cancer cells

Cisplatin, a commonly used chemotherapeutic, is associated with ototoxicity, renal toxicity and neurotoxicity, thus identifying means to increase the therapeutic index of cisplatin may allow for improved outcomes. A SNP (rs4343077) within EPS8, discovered through a genome wide association study of cisplatin-induced cytotoxicity and apoptosis in lymphoblastoid cell lines (LCLs), provided impetus to further study this gene. The purpose of this work was to evaluate the role of EPS8 in cellular susceptibility to cisplatin in cancerous and non-cancerous cells. We used EPS8 RNA interference to determine the effect of decreased EPS8 expression on LCL and A549 lung cancer cell sensitivity to cisplatin. EPS8 knockdown in LCLs resulted in a 7.9% increase in cisplatin-induced survival (P = 1.98 × 10(-7)) and an 8.7% decrease in apoptosis (P = 0.004) compared to control. In contrast, reduced EPS8 expression in lung cancer cells resulted in a 20.6% decrease in cisplatin-induced survival (P = 5.08 × 10(-5)). We then investigated an EPS8 inhibitor, mithramycin A, as a potential agent to increase the therapeutic index of cisplatin. Mithramycin A decreased EPS8 expression in LCLs resulting in decreased cellular sensitivity to cisplatin as evidenced by lower caspase 3/7 activation following cisplatin treatment (42.7% ± 6.8% relative to control P = 0.0002). In 5 non-small-cell lung carcinoma (NSCLC) cell lines, mithramycin A also resulted in decreased EPS8 expression. Adding mithramycin to 4 NSCLC cell lines and a bladder cancer cell line, resulted in increased sensitivity to cisplatin that was significantly more pronounced in tumor cell lines than in LCL lines (p<0.0001). An EGFR mutant NSCLC cell line (H1975) showed no significant change in sensitivity to cisplatin with the addition of mithramycin treatment. Therefore, an inhibitor of EPS8, such as mithramycin A, could improve cisplatin treatment by increasing sensitivity of tumor relative to normal cells.

Genetic and epigenetic variants contributing to clofarabine cytotoxicity

2-chloro-2-fluoro-deoxy-9-D-arabinofuranosyladenine (Clofarabine), a purine nucleoside analog, is used in the treatment of hematologic malignancies and as induction therapy for stem cell transplantation. The discovery of pharmacogenomic markers associated with chemotherapeutic efficacy and toxicity would greatly benefit the utility of this drug. Our objective was to identify genetic and epigenetic variants associated with clofarabine toxicity using an unbiased, whole genome approach. To this end, we employed International HapMap lymphoblastoid cell lines (190 LCLs) of European (CEU) or African (YRI) ancestry with known genetic information to evaluate cellular sensitivity to clofarabine. We measured modified cytosine levels to ascertain the contribution of genetic and epigenetic factors influencing clofarabine-mediated cytotoxicity. Association studies revealed 182 single nucleotide polymorphisms (SNPs) and 143 modified cytosines associated with cytotoxicity in both populations at the threshold P ≤ 0.0001. Correlation between cytotoxicity and baseline gene expression revealed 234 genes at P ≤ 3.98 × 10(-6). Six genes were implicated as: (i) their expression was directly correlated to cytotoxicity, (ii) they had a targeting SNP associated with cytotoxicity, and (iii) they had local modified cytosines associated with gene expression and cytotoxicity. We identified a set of three SNPs and three CpG sites targeting these six genes explaining 43.1% of the observed variation in phenotype. siRNA knockdown of the top three genes (SETBP1, BAG3, KLHL6) in LCLs revealed altered susceptibility to clofarabine, confirming relevance. As clofarabine's toxicity profile includes acute kidney injury, we examined the effect of siRNA knockdown in HEK293 cells. siSETBP1 led to a significant change in HEK293 cell susceptibility to clofarabine.

Cancer pharmacogenomics: early promise, but concerted effort needed

The past decade has brought together substantial advances in human genome analysis and a maturation of understanding of tumor biology. Although there is much progress still to be made, there are now several prominent examples in which tumor-associated somatic mutations have been used to identify cellular signaling pathways in tumors. This in turn has led to the development of targeted therapies, with somatic mutations serving as genomic predictors of tumor response and providing new leads for drug development. There is also a realization that germline DNA variants can help optimize cancer drug dosing and predict the susceptibility of patients to the adverse side effects of these drugs-knowledge that ultimately can be used to improve the benefit:risk ratio of cancer treatment for individual patients.

Integration of cell line and clinical trial genome-wide analyses supports a polygenic architecture of Paclitaxel-induced sensory peripheral neuropathy

PURPOSE

We sought to show the relevance of a lymphoblastoid cell line (LCL) model in the discovery of clinically relevant genetic variants affecting chemotherapeutic response by comparing LCL genome-wide association study (GWAS) results to clinical GWAS results.

EXPERIMENTAL DESIGN

A GWAS of paclitaxel-induced cytotoxicity was conducted in 247 LCLs from the HapMap Project and compared with a GWAS of sensory peripheral neuropathy in patients with breast cancer (n = 855) treated with paclitaxel in the Cancer and Leukemia Group B (CALGB) 40101 trial. Significant enrichment was assessed by permutation resampling analysis.

RESULTS

We observed an enrichment of LCL cytotoxicity-associated single-nucleotide polymorphisms (SNP) in the sensory peripheral neuropathy-associated SNPs from the clinical trial with concordant allelic directions of effect (empirical P = 0.007). Of the 24 SNPs that overlap between the clinical trial (P < 0.05) and the preclinical cytotoxicity study (P < 0.001), 19 of them are expression quantitative trait loci (eQTL), which is a significant enrichment of this functional class (empirical P = 0.0447). One of these eQTLs is located in RFX2, which encodes a member of the DNA-binding regulatory factor X family. Decreased expression of this gene by siRNA resulted in increased sensitivity of Neuroscreen-1(NS-1; rat pheochromocytoma) cells to paclitaxel as measured by reduced neurite outgrowth and increased cytotoxicity, functionally validating the involvement of RFX2 in nerve cell response to paclitaxel.

CONCLUSIONS

The enrichment results and functional example imply that cellular models of chemotherapeutic toxicity may capture components of the underlying polygenic architecture of related traits in patients.

Pharmacogenomics of chemotherapeutic susceptibility and toxicity

The goal of personalized medicine is to tailor a patient's treatment strategy on the basis of his or her unique genetic make-up. The field of oncology is beginning to incorporate many of the strategies of personalized medicine, especially within the realm of pharmacogenomics, which is the study of how inter-individual genetic variation determines drug response or toxicity. A main objective of pharmacogenomics is to facilitate physician decision-making regarding optimal drug selection, dose and treatment duration on a patient-by-patient basis. Recent advances in genome-wide genotyping and sequencing technologies have supported the discoveries of a number of pharmacogenetic markers that predict response to chemotherapy. However, effectively implementing these pharmacogenetic markers in the clinic remains a major challenge. This review focuses on the contribution of germline genetic variation to chemotherapeutic toxicity and response, and discusses the utility of genome-wide association studies and use of lymphoblastoid cell lines (LCLs) in pharmacogenomic studies. Furthermore, we highlight several recent examples of genetic variants associated with chemotherapeutic toxicity or response in both patient cohorts and LCLs, and discuss the challenges and future directions of pharmacogenomic discovery for cancer treatment.