oncoPredict: an R package for predicting in vivo or cancer patient drug response and biomarkers from cell line screening data.

Cell line drug screening datasets can be utilized for a range of different drug discovery applications from drug biomarker discovery to building translational models of drug response. Previously, we described three separate methodologies to (1) correct for general levels of drug sensitivity to enable drug-specific biomarker discovery, (2) predict clinical drug response in patients and (3) associate these predictions with clinical features to perform in vivo drug biomarker discovery. Here, we unite and update these methodologies into one R package (oncoPredict) to facilitate the development and adoption of these tools. This new OncoPredict R package can be applied to various in vitro and in vivo contexts for drug and biomarker discovery.

Institutional Profile: Pharmacogenomic research in R Stephanie Huang Laboratory.

The Huang Lab was established in 2009 at the University of Chicago and has since been active in conducting pharmacogenomic research. Our laboratory's main research focus is translational pharmacogenomics with a particular interest in the pharmacogenomics of anticancer agents. By systematically evaluating the human genome and its relationships to drug response and toxicity, our goal is to develop clinically useful models that predict risk for adverse drug reactions and nonresponse prior to administration of chemotherapy. Specifically, the theme of our research evolved around the idea of cell-based pharmacogenomics, which utilizes in vitro models for biomarker discovery and prediction-model construction, followed by in vivo validation. We routinely use cell lines (derived from healthy and diseased individuals as well as commercially available cancer cell lines) and clinical samples to discover and functionally characterize genetic variation and gene, miRNA, and long noncoding RNA expression for their roles in drug sensitivity.

Discovering drugs to overcome chemoresistance in ovarian cancers based on the cancer genome atlas tumor transcriptome profile.

Ovarian cancer accounts for the highest mortality among gynecologic cancers, mainly due to intrinsic or acquired chemoresistance. While mechanistic-based methods have been used to identify compounds that can overcome chemoresistance, an effective comprehensive drug screening has yet to be developed. We applied a transcriptome based drug sensitivity prediction method, to the Cancer Genome Atlas (TCGA) ovarian cancer dataset to impute patient tumor response to over 100 different drugs. By stratifying patients based on their predicted response to standard of care (SOC) chemotherapy, we identified drugs that are likely more sensitive in SOC resistant ovarian tumors. Five drugs (ABT-888, BIBW2992, gefitinib, AZD6244 and lenalidomide) exhibit higher efficacy in SOC resistant ovarian tumors when multi-platform of transcriptome profiling methods were employed. Additional in vitro and clinical sample validations were carried out and verified the effectiveness of these agents. Our candidate drugs hold great potential to improve clinical outcome of chemoresistant ovarian cancer.

Integrative analyses of genetic variation, epigenetic regulation, and the transcriptome to elucidate the biology of platinum sensitivity.

BACKGROUND: Using genome-wide genetic, gene expression, and microRNA expression (miRNA) data, we developed an integrative approach to investigate the genetic and epigenetic basis of chemotherapeutic sensitivity. RESULTS: Through a sequential multi-stage framework, we identified genes and miRNAs whose expression correlated with platinum sensitivity, mapped these to genomic loci as quantitative trait loci (QTLs), and evaluated the associations between these QTLs and platinum sensitivity. A permutation analysis showed that top findings from our approach have a much lower false discovery rate compared to those from a traditional GWAS of drug sensitivity. Our approach identified five SNPs associated with 10 miRNAs and the expression level of 15 genes, all of which were associated with carboplatin sensitivity. Of particular interest was one SNP (rs11138019), which was associated with the expression of both miR-30d and the gene ABCD2, which were themselves correlated with both carboplatin and cisplatin drug-specific phenotype in the HapMap samples. Functional study found that knocking down ABCD2 in vitro led to increased apoptosis in ovarian cancer cell line SKOV3 after cisplatin treatment. Over-expression of miR-30d in vitro caused a decrease in ABCD2 expression, suggesting a functional relationship between the two. CONCLUSIONS: We developed an integrative approach to the investigation of the genetic and epigenetic basis of human complex traits. Our approach outperformed standard GWAS and provided hints at potential biological function. The relationships between ABCD2 and miR-30d, and ABCD2 and platin sensitivity were experimentally validated, suggesting a functional role of ABCD2 and miR-30d in sensitivity to platinating agents.

Apolipoprotein A5 gene variants and the risk of coronary heart disease: a case­control study and meta­analysis.

Previous studies have shown that apolipoprotein A5 (APOA5) gene variants are genetic determinants of the concentration of triglycerides, which are a known risk factor for coronary heart disease (CHD). Using the standardized coronary angiography method, 290 CHD patients and 198 non­CHD controls were recruited from Ningbo Lihuili Hospital. In addition, 331 unrelated healthy volunteers were recruited as healthy controls from Ningbo Ximen Community residents. Three variants of the APOA5 gene, S19W, ­1131T>C and 553G>T, were analyzed for their association with CHD. Under a dominant inheritance model, ­1131CT>C was shown to be a CHD risk factor (P=0.030; OR, 1.422; 95% CI, 1.036­1.952). The single nucleotide polymorphism, 553G>T, was found to correlate with the severity of CHD in males (P=0.032). Meta­analysis showed that ­1131T>C was significantly associated with CHD (P<0.0001). By contrast, negative correlations with CHD were observed for S19W and 553G>T. In the present case­control study, APOA5 gene variants were not found to correlate with the risk of CHD in the populations studied; however, ­1131CT>C was shown to be a CHD risk factor under a dominant inheritance model. Meta­analysis showed a significant contribution of ­1131T>C to the risk of CHD, implying an ethnic difference in APOA5 gene variants.

Genetic variation that predicts platinum sensitivity reveals the role of miR-193b* in chemotherapeutic susceptibility.

Platinum agents are the backbone of cancer chemotherapy. Recently, we identified and replicated the role of a single nucleotide polymorphism (SNP, rs1649942) in predicting platinum sensitivity both in vitro and in vivo. Using the CEU samples from the International HapMap Project, we found the same SNP to be a master regulator of multiple gene expression phenotypes, prompting us to investigate whether rs1649942-mediated regulation of miRNAs may in part contribute to variation in platinum sensitivity. To these ends, 60 unrelated HapMap CEU I/II samples were used for our discovery-phase study using high-throughput genome-wide miRNA and gene expression profiling. Examining the relationships among rs1649942, its gene expression targets, genome-wide miRNA expression, and cellular sensitivity to carboplatin and cisplatin, we identified 2 platinum-associated miRNAs (miR-193b* and miR-320) that inhibit the expression of 5 platinum-associated genes (CRIM1, IFIT2, OAS1, KCNMA1, and GRAMD1B). We further replicated the relationship between the expression of miR-193b*, CRIM1, IFIT2, KCNMA1, and GRAMD1B, and platinum sensitivity in a separate HapMap CEU III dataset. We then showed that overexpression of miR-193b* in a randomly selected HapMap cell line results in resistance to both carboplatin and cisplatin. This relationship was also found in 7 ovarian cancer cell lines from NCI60 dataset and confirmed in an OVCAR-3 that overexpression of miR-193b* leads to increased resistance to carboplatin. Our findings highlight a potential mechanism of action for a previously observed genotype-survival outcome association. Further examination of miR-193b* in platinum sensitivity in ovarian cancer is warranted.

Chemotherapeutic-induced apoptosis: a phenotype for pharmacogenomics studies.

AIM: To determine whether cellular apoptosis is a suitable phenotypic trait for pharmacogenomics studies by evaluating caspase 3/7-mediated activity in lymphoblastoid cell lines after treatment with six chemotherapeutic agents: 5'-deoxyfluorouridine, pemetrexed, cytarabine, paclitaxel, carboplatin, and cisplatin. MATERIALS AND METHODS: Using monozygotic twin pair and sibling pair lymphoblastoid cell lines, we identified conditions for measurement of caspase 3/7 activity in lymphoblastoid cell lines. Genome-wide association studies were performed with over 2 million single nucleotide polymorphisms (SNPs) and cisplatin-induced apoptosis in HapMap CEU cell lines (n=77). RESULTS: Although treatment with 5'-deoxyfluorouridine and pemetrexed for up to 24 h resulted in low levels of apoptosis or interindividual variation in caspase-dependent cell death; paclitaxel, cisplatin, carboplatin, and cytarabine treatment for 24 h resulted in 9.4-fold, 9.1-fold, 7.0-fold, and 6.0-fold increases in apoptosis relative to control, respectively. There was a weak correlation between caspase activity and cytotoxicity (r(2)=0.03-0.29) demonstrating that cytotoxicity and apoptosis are two distinct phenotypes that may produce independent genetic associations. Estimated heritability (h(2)) for apoptosis was 0.57 and 0.29 for cytarabine (5 and 40 µmol/l, respectively), 0.22 for paclitaxel (12.5 nmol/l), and 0.34 for cisplatin (5 µmol/l). In the genome-wide association study using the HapMap CEU panel, we identified a significant enrichment of cisplatin-induced cytotoxicity SNPs within the significant cisplatin-induced apoptosis SNPs and an enrichment of expression quantitative trait loci (eQTL). Among these eQTLs, we identified several eQTLs with known function related to apoptosis and/or cytotoxicity. CONCLUSION: Our study identifies apoptosis as a phenotype for pharmacogenomic studies in lymphoblastoid cell lines after treatment with paclitaxel, cisplatin, carboplatin, and cytarabine that may have utility for discovering biomarkers to predict response to certain chemotherapeutics.

Comprehensive evaluation of the contribution of X chromosome genes to platinum sensitivity.

Using a genome-wide gene expression data set generated from Affymetrix GeneChip Human Exon 1.0ST array, we comprehensively surveyed the role of 322 X chromosome gene expression traits on cellular sensitivity to cisplatin and carboplatin. We identified 31 and 17 X chromosome genes whose expression levels are significantly correlated (after multiple testing correction) with sensitivity to carboplatin and cisplatin, respectively, in the combined HapMap CEU (Utah residents with ancestry from northern and western Europe) and YRI (Yoruba in Ibahan, Nigeria) populations (false discovery rate, FDR < 0.05). Of those, 14 overlap for both cisplatin and carboplatin. Using an independent gene expression quantification method, the Illumina Sentrix Human-6 Expression BeadChip, measured on the same HapMap cell lines, we found that 4 and 2 of these genes are significantly associated with carboplatin and cisplatin sensitivity, respectively, in both analyses. Two genes, CTPS2 and DLG3, were identified by both genome-wide gene expression analyses as correlated with cellular sensitivity to both platinating agents. The expression of DLG3 gene was also found to correlate with cellular sensitivity to platinating agents in NCI-60 cancer cell lines. In addition, we evaluated whether the expression of X chromosome genes contributed to the observed differences in sensitivity to the platinums between CEU and YRI-derived cell lines. Of the 34 distinct genes significantly correlated with either carboplatin or cisplatin sensitivity, 14 are differentially expressed (defined as P < 0.05) between CEU and YRI. Thus, sex chromosome genes play a role in cellular sensitivity to platinating agents and differences in the expression level of these genes are an important source of variation that should be included in comprehensive pharmacogenomic studies.

PACdb: a database for cell-based pharmacogenomics.

We have developed Pharmacogenomics And Cell database (PACdb), a results database that makes available relationships between single nucleotide polymorphisms, gene expression, and cellular sensitivity to various drugs in cell-based models to help determine genetic variants associated with drug response. The current version also supports summary analysis on differentially expressed genes between the HapMap samples of European and African ancestry, as well as queries for summary information of correlations between gene expression and pharmacological phenotypes. At present, data generated on the following anticancer agents are included: carboplatin, cisplatin, etoposide, daunorubicin, and cytarabine (Ara-C). The database is also available to assist in the investigation of the effects of potential confounding variables (e.g. cell proliferation rate) in lymphoblastoid cell lines. PACdb will be regularly updated to include more drugs and new datasets (e.g. baseline microRNA levels). PACdb will be linked into PharmGKB to benefit the next wave of pharmacogenetic and pharmacogenomic discovery.

Identification of genetic variants and gene expression relationships associated with pharmacogenes in humans.

OBJECTIVES: The very important pharmacogenes (VIPs) were selected by Pharmacogenetic Research Network (National Institutes of Health-PGRN) owing to their significant effects on drug treatment both at the pharmacokinetic and pharmacodynamic levels. Our objective was to identify single nucleotide polymorphisms (SNPs) that potentially affected the expression of these genes or potential SNP-gene interactions involved to improve our understanding of genetic effects on drug therapy. BASIC METHODS: Gene expression was evaluated in 176 International HapMap lymphoblastoid cell lines derived from CEU (CEPH, Utah residents with ancestry from northern and western Europe; n=87) and YRI (Yoruba in Ibadan, Nigeria; n=89) using Affymetrix GeneChip Human Exon 1.0 ST arrays (Affymetrix Laboratory, Affymetrix Inc., Santa Clara, California, USA) with interrogation of greater than 17,000 human genes. Genome-wide association was performed between over two million publicly available HapMap SNPs and gene expression. MAIN RESULTS: The expression of two PGRN-VIPs (GSTT1 and GSTM1) are significantly associated with SNPs within 2.5 Mb of the genes; whereas the expression of three and ten PGRN-VIPs are significantly associated with distant-acting SNPs in CEU and YRI, respectively. In addition, three and four PGRN-VIPs harbor SNPs that are distantly associated with other gene expressions in CEU and YRI, respectively. PRINCIPAL CONCLUSION: Using this information, one may identify genetic variants that are significantly associated with the expression of any set of genes of interest; or evaluate potential gene-gene interaction through SNP expression relationships.

Mapping genes that contribute to daunorubicin-induced cytotoxicity.

Daunorubicin is an anthracycline antibiotic agent used in the treatment of hematopoietic malignancies. Toxicities associated with this agent include myelosuppression and cardiotoxicity; however, the genes or genetic determinants that contribute to these toxicities are unknown. We present an unbiased genome-wide approach that incorporates heritability, whole-genome linkage analysis, and linkage-directed association to uncover genetic variants contributing to the sensitivity to daunorubicin-induced cytotoxicity. Cell growth inhibition in 324 Centre d' Etude du Polymorphisme Humain lymphoblastoid cell lines (24 pedigrees) was evaluated following treatment with daunorubicin for 72 h. Heritability analysis showed a significant genetic component contributing to the cytotoxic phenotypes (h2 = 0.18-0.63 at 0.0125, 0.025, 0.05, 0.1, 0.2, and 1.0 mumol/L daunorubicin and at the IC50, the dose required to inhibit 50% cell growth). Whole-genome linkage scans at all drug concentrations and IC50 uncovered 11 regions with moderate peak LOD scores (> 1.5), including 4q28.2 to 4q32.3 with a maximum LOD score of 3.18. The quantitative transmission disequilibrium tests were done using 31,312 high-frequency single-nucleotide polymorphisms (SNP) located in the 1 LOD confidence interval of these 11 regions. Thirty genes were identified as significantly associated with daunorubicin-induced cytotoxicity (P < or = 2.0 x 10(-4), false discovery rate < or = 0.1). Pathway and functional gene ontology analysis showed that these genes were overrepresented in the phosphatidylinositol signaling system, axon guidance pathway, and GPI-anchored proteins family. Our findings suggest that a proportion of susceptibility to daunorubicin-induced cytotoxicity may be controlled by genetic determinants and that analysis using linkage-directed association studies with dense SNP markers can be used to identify the genetic variants contributing to cytotoxicity.

Effect of population and gender on chemotherapeutic agent-induced cytotoxicity.

Large interindividual variance is observed in both response and toxicity associated with chemotherapy. Our goal is to identify factors that contribute to chemotherapy-induced toxicity. To this end, we used EBV-transformed B-lymphoblastoid HapMap cell lines derived from 30 Yoruban trios (African descent) and 30 Centre d' Etude du Polymorphisme Humain (CEPH) trios (European descent) to evaluate population- and gender-specific differences in cytotoxicity of carboplatin, cisplatin, daunorubicin, and etoposide using a high-throughput, short-term cytotoxicity assay. The IC(50) was compared for population- and gender-specific differences for the four drugs. We observed large interindividual variance in IC(50) values for carboplatin, cisplatin, daunorubicin, and etoposide for both Yoruban and CEPH populations (range from 8- to 433-fold). Statistically significant differences in carboplatin and daunorubicin IC(50) were shown when comparing Yoruban cell lines (n = 89) to CEPH cell lines (n = 87; P = 0.002 and P = 0.029, respectively). This population difference in treatment induced cytotoxicity was not seen for either cisplatin or etoposide. In the Yoruban population, cell lines derived from females were less sensitive to platinating agents than males [median carboplatin IC(50), 29.1 versus 24.6 micromol/L (P = 0.012); median cisplatin IC(50), 7.0 versus 6.0 micromol/L (P = 0.020) in female and male, respectively]. This difference was not observed in the CEPH population. These results show that population and gender may affect risk for toxicities associated with certain chemotherapeutic agents.