Identification of genetic variants contributing to cisplatin-induced cytotoxicity by use of a genomewide approach

Pharmacogenomic predictor of long-term residual chemotherapy-induced peripheral neuropathy in ovarian cancer survivors: A substudy of the GINECO Vivrovaire study

BACKGROUND

Chemotherapy (CT) remains a backbone treatment of epithelial ovarian cancer (EOC) inducing persistent peripheral neuropathy (CIPN). Using a dedicated patient-reported outcome tool, this study investigated persistent CIPN and its pharmacogenetic predictors in a cohort of long-term EOC survivors.

METHODS

Vivrovaire was a French multicenter cohort of patients with EOC free of disease 3 years after CT completion. Persistent CIPN was assessed using the FACT/GOG-Ntx4 self-questionnaire. The association of homozygous (hom) or heterozygous (het) single nucleotide polymorphisms (SNPs) in selected genes was evaluated.

RESULTS

130 patients were included with a median time from CT completion of 63 [35-180] months. The median CIPN score was 37 [18-44], with 35 (26.9%) patients reporting severe CIPN (<33). SNPs were identified as follows: CYP2C8 [hom, n = 32 (24.6%)/het, n = 99, (76.2%)]; CYP3A4 [hom, n = 0 (0%)/het, n = 8 (6.2%)], ERCC1 [hom, n = 21 (16.2%)/het, n = 57 (43.8%)], and XPC [hom, n = 45 (34.6%)/het, n = 66 (50.8%)]. In univariate analysis, the identification of ≥1 hom SNP was associated with a lower CIPN score (continuous variable; p = 0.045). Patients harboring hom or het CYP2C8_rs1934951 SNP reported more likely severe CIPN (threshold <33) score (OR 2.482; 95% CI [1.126-5.47], p = 0.024). In the multivariate analyses, age, interval from CT completion, type and number of CT courses were not significantly associated with CIPN score (OR 5.165, 95% CI [0.478-55.83], p = 0.176).

CONCLUSIONS

Persistent CIPN is common among ovarian cancer long-term survivors. CYP2C8_rs1934951 SNP may be associated with severe residual CIPN in EOC survivors. More studies are warranted to identify predictive factors of CIPN.

Genomic profile of Parkinson's disease in Asians

Parkinson's Disease (PD) has witnessed an alarming rise in prevalence, highlighting the suboptimal nature of early diagnostic and therapeutic strategies. To address this issue, genetic testing has emerged as a potential avenue. In this comprehensive review, we have meticulously summarized the variants associated with PD in Asian populations. Our review reveals that these variants exert their influence on diverse biological pathways, encompassing the autophagy-lysosome pathway, cholesterol metabolism, circadian rhythm regulation, immune system response, and synaptic function. Conventionally, PD has been linked to other diseases; however, our findings shed light on a shared genetic susceptibility among these conditions, implying an underlying pathophysiological mechanism that unifies them. Moreover, it is noteworthy that these PD-associated variants can significantly impact drug responses during therapeutic interventions. This review not only provides a consolidated overview of the genetic variants associated with PD in Asian populations but also contributes novel insights into the intricate relationships between PD and other diseases by elucidating shared genetic components. These findings underscore the importance of personalized approaches in diagnosing and treating PD based on individual genetic profiles to optimize patient outcomes.

The Expression of Trace Amine-Associated Receptors (TAARs) in Breast Cancer Is Coincident with the Expression of Neuroactive Ligand-Receptor Systems and Depends on Tumor Intrinsic Subtype

Currently, the contribution of trace amine-associated receptors (TAARs) to breast cancer (BC) is recognized, but their associations with various pathological characteristics are not yet understood. There is accumulated transcriptomic data for BC tumors, which are represented in publicly accessible databases. We estimated TAARs' (including TAAR1, TAAR2, TAAR5, TAAR6, TAAR8, and TAAR9) associations with BC stage, grade, and molecular subtypes in these data and identified that the expression of all TAARs was associated with more unfavorable cancer subtypes, including basal-like and HER2-positive tumors. Also, the significant upregulation of all TAARs was demonstrated in circulating tumor cells compared to the metastatic lesions. Considering that co-expressed genes are more likely to be involved in the same biologic processes, we analyzed genes that are co-expressed with TAARs in BC. These gene sets were enriched with the genes of the olfactory transduction pathway and neuroactive ligand-receptor interaction participants. TAARs are co-expressed with G-protein-coupled receptors of monoamine neurotransmitters including dopamine, norepinephrine, and serotonin as well as with other neuroactive ligand-specific receptors. Since TAAR1 is able to modulate the activity of monoamine receptors that are involved in the regulation of BC growth, TAAR1 and potentially other TAARs may be regarded as prospective therapeutic targets for breast cancer.

Uncovering a Genetic Polymorphism Located in Huntingtin Associated Protein 1 in Modulation of Central Pain Sensitization Signaling Pathways

Fibromyalgia syndrome (FMS) is characterized by widespread pain and increased sensitivity to nociceptive stimulus or tenderness. While familial aggregation could suggest a potential hereditary component in FMS development, isolation of genetic determinants has proven difficult due to the multi-factorial nature and complexity of the syndrome. Central sensitization is thought to be one of the key mechanisms leading to FMS in a subset of patients. Enhanced central pain signaling can be measured using the Nociceptive Flexion Reflex (NFR) or RIII threshold. We performed a genome-wide association study (GWAS) using an array to genotype 258,756 human genetic polymorphisms in 225 FMS patients and 77 healthy volunteers and searched for genetic variants associated with a lowered NFR threshold. We have identified a potential association between a single nucleotide polymorphism resulting in a common non-synonymous coding mutation in the Huntingtin associated protein 1 (HAP1) gene (rs4796604, MAF = 0.5) and the NFR threshold (p = 4.78E−06). The Hap1 protein is involved in trafficking and is particularly enriched in neurons. Our results suggest a possible involvement of the neuronal trafficking protein HAP1 in modulating pain signaling pathways and thus participate in the establishment of the NFR threshold.

Pattern of TAAR5 Expression in the Human Brain Based on Transcriptome Datasets Analysis

Trace amine-associated receptors (TAAR) recognize organic compounds, including primary, secondary, and tertiary amines. The TAAR5 receptor is known to be involved in the olfactory sensing of innate socially relevant odors encoded by volatile amines. However, emerging data point to the involvement of TAAR5 in brain functions, particularly in the emotional behaviors mediated by the limbic system which suggests its potential contribution to the pathogenesis of neuropsychiatric diseases. TAAR5 expression was explored in datasets available in the Gene Expression Omnibus, Allen Brain Atlas, and Human Protein Atlas databases. Transcriptomic data demonstrate ubiquitous low TAAR5 expression in the cortical and limbic brain areas, the amygdala and the hippocampus, the nucleus accumbens, the thalamus, the hypothalamus, the basal ganglia, the cerebellum, the substantia nigra, and the white matter. Altered TAAR5 expression is identified in Down syndrome, major depressive disorder, or HIV-associated encephalitis. Taken together, these data indicate that TAAR5 in humans is expressed not only in the olfactory system but also in certain brain structures, including the limbic regions receiving olfactory input and involved in critical brain functions. Thus, TAAR5 can potentially be involved in the pathogenesis of brain disorders and represents a valuable novel target for neuropsychopharmacology.

Advances in Our Understanding of the Molecular Mechanisms of Action of Cisplatin in Cancer Therapy

Cisplatin and other platinum-based chemotherapeutic drugs have been used extensively for the treatment of human cancers such as bladder, blood, breast, cervical, esophageal, head and neck, lung, ovarian, testicular cancers, and sarcoma. Cisplatin is commonly administered intravenously as a first-line chemotherapy for patients suffering from various malignancies. Upon absorption into the cancer cell, cisplatin interacts with cellular macromolecules and exerts its cytotoxic effects through a series of biochemical mechanisms by binding to Deoxyribonucleic acid (DNA) and forming intra-strand DNA adducts leading to the inhibition of DNA synthesis and cell growth. Its primary molecular mechanism of action has been associated with the induction of both intrinsic and extrinsic pathways of apoptosis resulting from the production of reactive oxygen species through lipid peroxidation, activation of various signal transduction pathways, induction of p53 signaling and cell cycle arrest, upregulation of pro-apoptotic genes/proteins, and down-regulation of proto-oncogenes and anti-apoptotic genes/proteins. Despite great clinical outcomes, many studies have reported substantial side effects associated with cisplatin monotherapy, while others have shown substantial drug resistance in some cancer patients. Hence, new formulations and several combinational therapies with other drugs have been tested for the purpose of improving the clinical utility of cisplatin. Therefore, this review provides a comprehensive understanding of its molecular mechanisms of action in cancer therapy and discusses the therapeutic approaches to overcome cisplatin resistance and side effects.

Response and Toxicity to Cytarabine Therapy in Leukemia and Lymphoma: From Dose Puzzle to Pharmacogenomic Biomarkers

Simple Summary

In this review, the authors propose a crosswise examination of cytarabine-related issues ranging from the spectrum of clinical activity and severe toxicities, through updated cellular pharmacology and drug formulations, to the genetic variants associated with drug-induced phenotypes. Cytarabine (cytosine arabinoside; Ara-C) in multiagent chemotherapy regimens is often used for leukemia or lymphoma treatments, as well as neoplastic meningitis. Chemotherapy regimens can induce a suboptimal clinical outcome in a fraction of patients. The individual variability in clinical response to Leukemia & Lymphoma treatments among patients appears to be associated with intracellular accumulation of Ara-CTP due to genetic variants related to metabolic enzymes. The review provides exhaustive information on the effects of Ara-C-based therapies, the adverse drug reaction will also be provided including bone pain, ocular toxicity (corneal pain, keratoconjunctivitis, and blurred vision), maculopapular rash, and occasional chest pain. Evidence for predicting the response to cytarabine-based treatments will be highlighted, pointing at their significant impact on the routine management of blood cancers.

Abstract

Cytarabine is a pyrimidine nucleoside analog, commonly used in multiagent chemotherapy regimens for the treatment of leukemia and lymphoma, as well as for neoplastic meningitis. Ara-C-based chemotherapy regimens can induce a suboptimal clinical outcome in a fraction of patients. Several studies suggest that the individual variability in clinical response to Leukemia & Lymphoma treatments among patients, underlying either Ara-C mechanism resistance or toxicity, appears to be associated with the intracellular accumulation and retention of Ara-CTP due to genetic variants related to metabolic enzymes. Herein, we reported (a) the latest Pharmacogenomics biomarkers associated with the response to cytarabine and (b) the new drug formulations with optimized pharmacokinetics. The purpose of this review is to provide readers with detailed and comprehensive information on the effects of Ara-C-based therapies, from biological to clinical practice, maintaining high the interest of both researcher and clinical hematologist. This review could help clinicians in predicting the response to cytarabine-based treatments.

Genetically regulated expression underlies cellular sensitivity to chemotherapy in diverse populations

Most cancer chemotherapeutic agents are ineffective in a subset of patients; thus, it is important to consider the role of genetic variation in drug response. Lymphoblastoid cell lines (LCLs) in 1000 Genomes Project populations of diverse ancestries are a useful model for determining how genetic factors impact the variation in cytotoxicity. In our study, LCLs from three 1000 Genomes Project populations of diverse ancestries were previously treated with increasing concentrations of eight chemotherapeutic drugs, and cell growth inhibition was measured at each dose with half-maximal inhibitory concentration (IC50) or area under the dose-response curve (AUC) as our phenotype for each drug. We conducted both genome-wide association studies (GWAS) and transcriptome-wide association studies (TWAS) within and across ancestral populations. We identified four unique loci in GWAS and three genes in TWAS to be significantly associated with the chemotherapy-induced cytotoxicity within and across ancestral populations. In the etoposide TWAS, increased STARD5 predicted expression associated with decreased etoposide IC50 (P = 8.5 × 10-8). Functional studies in A549, a lung cancer cell line, revealed that knockdown of STARD5 expression resulted in the decreased sensitivity to etoposide following exposure for 72 (P = 0.033) and 96 h (P = 0.0001). By identifying loci and genes associated with cytotoxicity across ancestral populations, we strive to understand the genetic factors impacting the effectiveness of chemotherapy drugs and to contribute to the development of future cancer treatment.

Progenitor cell therapy for acquired pediatric nervous system injury: Traumatic brain injury and acquired sensorineural hearing loss

While cell therapies hold remarkable promise for replacing injured cells and repairing damaged tissues, cell replacement is not the only means by which these therapies can achieve therapeutic effect. For example, recent publications show that treatment with varieties of adult, multipotent stem cells can improve outcomes in patients with neurological conditions such as traumatic brain injury and hearing loss without directly replacing damaged or lost cells. As the immune system plays a central role in injury response and tissue repair, we here suggest that multipotent stem cell therapies achieve therapeutic effect by altering the immune response to injury, thereby limiting damage due to inflammation and possibly promoting repair. These findings argue for a broader understanding of the mechanisms by which cell therapies can benefit patients.

Pathway-Based Drug-Repurposing Schemes in Cancer: The Role of Translational Bioinformatics

Cancer is a set of complex pathologies that has been recognized as a major public health problem worldwide for decades. A myriad of therapeutic strategies is indeed available. However, the wide variability in tumor physiology, response to therapy, added to multi-drug resistance poses enormous challenges in clinical oncology. The last years have witnessed a fast-paced development of novel experimental and translational approaches to therapeutics, that supplemented with computational and theoretical advances are opening promising avenues to cope with cancer defiances. At the core of these advances, there is a strong conceptual shift from gene-centric emphasis on driver mutations in specific oncogenes and tumor suppressors-let us call that the silver bullet approach to cancer therapeutics-to a systemic, semi-mechanistic approach based on pathway perturbations and global molecular and physiological regulatory patterns-we will call this the shrapnel approach. The silver bullet approach is still the best one to follow when clonal mutations in driver genes are present in the patient, and when there are targeted therapies to tackle those. Unfortunately, due to the heterogeneous nature of tumors this is not the common case. The wide molecular variability in the mutational level often is reduced to a much smaller set of pathway-based dysfunctions as evidenced by the well-known hallmarks of cancer. In such cases "shrapnel gunshots" may become more effective than "silver bullets". Here, we will briefly present both approaches and will abound on the discussion on the state of the art of pathway-based therapeutic designs from a translational bioinformatics and computational oncology perspective. Further development of these approaches depends on building collaborative, multidisciplinary teams to resort to the expertise of clinical oncologists, oncological surgeons, and molecular oncologists, but also of cancer cell biologists and pharmacologists, as well as bioinformaticians, computational biologists and data scientists. These teams will be capable of engaging on a cycle of analyzing high-throughput experiments, mining databases, researching on clinical data, validating the findings, and improving clinical outcomes for the benefits of the oncological patients.

Association of Genetic Variants at TRPC6 With Chemotherapy-Related Heart Failure

Background: Our previous GWAS identified genetic variants at six novel loci that were associated with a decline in left ventricular ejection fraction (LVEF), p < 1 × 10⁻⁵ in 1,191 early breast cancer patients from the N9831 clinical trial of chemotherapy plus trastuzumab. In this study we sought replication of these loci.

Methods: We tested the top loci from the GWAS for association with chemotherapy-related heart failure (CRHF) using 26 CRHF cases from N9831 and 984 patients from the Mayo Clinic Biobank which included CRHF cases (N = 12) and control groups of patients treated with anthracycline +/– trastuzumab without HF (N = 282) and patients with HF that were never treated with anthracycline or trastuzumab (N = 690). We further examined associated loci in the context of gene expression and rare coding variants using a TWAS approach in heart left ventricle and Sanger sequencing, respectively. Doxorubicin-induced apoptosis and cardiomyopathy was modeled in human iPSC-derived cardiomyocytes and endothelial cells and a mouse model, respectively, that were pre-treated with GsMTx-4, an inhibitor of TRPC6.

Results:TRPC6 5′ flanking variant rs57242572-T was significantly more frequent in cases compared to controls, p = 0.031, and rs61918162-T showed a trend for association, p = 0.065. The rs61918162 T-allele was associated with higher TRPC6 expression in the heart left ventricle. We identified a single TRPC6 rare missense variant (rs767086724, N338S, prevalence 0.0025% in GnomAD) in one of 38 patients (2.6%) with CRHF. Pre-treatment of cardiomyocytes and endothelial cells with GsMTx4 significantly reduced doxorubicin-induced apoptosis. Similarly, mice treated with GsMTx4 had significantly improved doxorubicin-induced cardiac dysfunction.

Conclusions: Genetic variants that are associated with increased TRPC6 expression in the heart and rare TRPC6 missense variants may be clinically useful as risk factors for CRHF. GsMTx-4 may be a cardioprotective agent in patients with TRPC6 risk variants. Replication of the genetic associations in larger well-characterized samples and functional studies are required.

Pharmacogenomics in epithelial ovarian cancer first-line treatment outcome: validation of GWAS-associated NRG3 rs1649942 and BRE rs7572644 variants in an independent cohort

The identification of predictive biomarkers for the first-line treatment of epithelial ovarian cancer (EOC) remains a challenge. Although genome-wide association studies (GWAS) have identified several genetic polymorphisms as predictors of EOC clinical outcome, the subsequent validation has not yet been performed. This study aims to validate the influence of Neuregulin 3 (NRG3) rs1649942 and Brain and reproductive organ-expressed (TNFRSF1A modulator) (BRE) rs7572644 GWAS-identified variants in an independent cohort of EOC patients from the North region of Portugal (n = 339) submitted to first-line treatment. Polymorphism genotypes were determined by real-time PCR using validated assays. Patients carrying the NRG3 rs1649942 A allele presented a significantly longer overall survival (OS) when compared to GG-genotype patients (log-rank test, P = 0.011) in the FIGO IV stage subgroup. No impact was observed for early-stage patients or considering disease-free survival (DFS) as an outcome. For FIGO I/II stage patients, BRE rs7572644 C allele carriers exhibit a decreased OS (P = 0.014) and DFS (P = 0.032) when compared to TT-homozygous patients. Furthermore, a Multivariate Cox regression analysis revealed a three-fold increase in the risk of death (HR, 3.09; P = 0.015) and recurrence (HR, 3.33; P = 0.009) for FIGO I/II C allele carriers. No significant impact was observed for late-stage patients. The BRE rs7572644 and NRG3 rs1649942 genetic variants were validated in an independent cohort of EOC Portuguese patients, particularly in specific subgroups considering FIGO staging. Further functional post-GWAS analyses are indispensable to understand the biological mechanisms underlying the observed results.

MTOR Pathway-Based Discovery of Genetic Susceptibility to L-DOPA-Induced Dyskinesia in Parkinson's Disease Patients

Dyskinesia induced by L-DOPA administration (LID) is one of the most invalidating adverse effects of the gold standard treatment restoring dopamine transmission in Parkinson's disease (PD). However, LID manifestation in parkinsonian patients is variable and heterogeneous. Here, we performed a candidate genetic pathway analysis of the mTOR signaling cascade to elucidate a potential genetic contribution to LID susceptibility, since mTOR inhibition ameliorates LID in PD animal models. We screened 64 single nucleotide polymorphisms (SNPs) mapping to 57 genes of the mTOR pathway in a retrospective cohort of 401 PD cases treated with L-DOPA (70 PD with moderate/severe LID and 331 with no/mild LID). We performed classic allelic, genotypic, and epistatic analyses to evaluate the association of individual or combinations of SNPs with LID onset and with LID severity after initiation of L-DOPA treatment. As for the time to LID onset, we found significant associations with SNP rs1043098 in the EIF4EBP2 gene and also with an epistatic interaction involving EIF4EBP2 rs1043098, RICTOR rs2043112, and PRKCA rs4790904. For LID severity, we found significant association with HRAS rs12628 and PRKN rs1801582 and also with a four-loci epistatic combination involving RPS6KB1 rs1292034, HRAS rs12628, RPS6KA2 rs6456121, and FCHSD1 rs456998. These findings indicate that the mTOR pathway contributes genetically to LID susceptibility. Our study could help to identify the most susceptible PD patients to L-DOPA in order to prevent the appearance of early and/or severe LID in a future. This information could also be used to stratify PD patients in clinical trials in a more accurate way.

Inner ear exosomes and their potential use as biomarkers

Exosomes are nanovesicles involved in intercellular communications. They are released by a variety of cell types; however, their presence in the inner ear has not been described in the literature. The aims of this study were to determine if exosomes are present in the inner ear and, if present, characterize the changes in their protein content in response to ototoxic stress. In this laboratory investigation, inner ear explants of 5-day-old Wistar rats were cultured and treated with either cisplatin or gentamicin. Hair cell damage was assessed by confocal microscopy. Exosomes were isolated using ExoQuick, serial centrifugation, and mini-column methods. Confirmation and characterization of exosomes was carried out using transmission electron microscopy (TEM), ZetaView, BCA protein analysis, and proteomics. Vesicles with a typical size distribution for exosomes were observed using TEM and ZetaView. Proteomic analysis detected typical exosome markers and markers for the organ of Corti. There was a statistically significant reduction in the exosome protein level and number of particles per cubic centimeter when the samples were exposed to ototoxic stress. Proteomic analysis also detected clear differences in protein expression when ototoxic medications were introduced. Significant changes in the proteomes of the exosomes were previously described in the context of hearing loss and ototoxic treatment. This is the first report describing exosomes derived from the inner ear. These findings may present an opportunity to conduct further studies with the hope of using exosomes as a biomarker to monitor inner ear function in the future.

DNA methylation of a novel PAK4 locus influences ototoxicity susceptibility following cisplatin and radiation therapy for pediatric embryonal tumors

Background

Ototoxicity is a common adverse side effect of platinum chemotherapy and cranial radiation therapy; however, individual susceptibility is highly variable. Therefore, our objective was to conduct an epigenome-wide association study to identify differentially methylated cytosine-phosphate-guanine (CpG) sites associated with ototoxicity susceptibility among cisplatin-treated pediatric patients with embryonal tumors.

Methods

Samples were collected for a discovery cohort (n = 62) and a replication cohort (n = 18) of medulloblastoma and primitive neuroectodermal tumor patients. Posttreatment audiograms were evaluated using the International Society of Paediatric Oncology (SIOP) Boston Ototoxicity Scale. Genome-wide associations between CpG methylation and ototoxicity were examined using multiple linear regression, controlling for demographic and treatment factors.

Results

The mean cumulative dose of cisplatin was 330 mg/m2 and the mean time from end of therapy to the last available audiogram was 6.9 years. In the discovery analysis of 435233 CpG sites, 6 sites were associated with ototoxicity grade (P < 5 × 10-5) after adjusting for confounders. Differential methylation at the top CpG site identified in the discovery cohort (cg14010619, PAK4 gene) was replicated (P = 0.029) and reached genome-wide significance (P = 2.73 × 10-8) in a combined analysis. These findings were robust to a sensitivity analysis evaluating other potential confounders.

Conclusions

We identified and replicated a novel CpG methylation loci (cg14010619) associated with ototoxicity severity. Methylation at cg14010619 may modify PAK4 activity, which has been implicated in cisplatin resistance in malignant cell lines.

Integration of genetic and functional genomics data to uncover chemotherapeutic induced cytotoxicity

Identifying genetic variants associated with chemotherapeutic induced toxicity is an important step towards personalized treatment of cancer patients. However, annotating and interpreting the associated genetic variants remains challenging because each associated variant is a surrogate for many other variants in the same region. The issue is further complicated when investigating patterns of associated variants with multiple drugs. In this study, we used biological knowledge to annotate and compare genetic variants associated with cellular sensitivity to mechanistically distinct chemotherapeutic drugs, including platinating agents (cisplatin, carboplatin), capecitabine, cytarabine, and paclitaxel. The most significantly associated SNPs from genome wide association studies of cellular sensitivity to each drug in lymphoblastoid cell lines derived from populations of European (CEU) and African (YRI) descent were analyzed for their enrichment in biological pathways and processes. We annotated genetic variants using higher-level biological annotations in efforts to group variants into more interpretable biological modules. Using the higher-level annotations, we observed distinct biological modules associated with cell line populations as well as classes of chemotherapeutic drugs. We also integrated genetic variants and gene expression variables to build predictive models for chemotherapeutic drug cytotoxicity and prioritized the network models based on the enrichment of DNA regulatory data. Several biological annotations, often encompassing different SNPs, were replicated in independent datasets. By using biological knowledge and DNA regulatory information, we propose a novel approach for jointly analyzing genetic variants associated with multiple chemotherapeutic drugs.

Pharmacogenetics of drug-induced ototoxicity caused by aminoglycosides and cisplatin

Aminoglycosides and the anticancer drug cisplatin can cause permanent hearing loss, which impacts patients' quality of life and results in considerable subsequent costs. Since patients' individual susceptibility to aminoglycoside- and cisplatin-induced ototoxicity varies considerably, strategies are needed to identify patients at risk, who may require alternative treatments or specific protection strategies. For both drugs, various genetic variants were linked to an increased or decreased risk for ototoxicity. Except for the association between the A1555G mitochondrial DNA mutation and aminoglycoside ototoxicity, their evidence is considered low because study cohorts were often small and replication studies either missing or contradictory. This review summarizes the pharmacogenetic markers linked to aminoglycoside- or cisplatin-induced ototoxicity and discusses reasons for replication failure and future perspective.

Rethinking ovarian cancer genomics: where genome-wide association studies stand?

Genome-wide association studies (GWAS) allow the finding of genetic variants associated with several traits. Regarding ovarian cancer (OC), 15 GWAS have been conducted since 2009, with the discovery of 49 SNPs associated with disease susceptibility and 46 with impact in the clinical outcome of patients (p < 5.00 × 10-2). Among them, 14 variants reached the genome-wide significance (p < 5.00 × 10−8). Despite the results obtained, they should be validated in independent sets. So far, five validation studies have been conducted which could confirm the association of 12 OC susceptibility SNPs. Consequently, post-GWAS studies are crucial unravel the biological plausibility of GWAS’ findings and the allelic spectrum of OC.

Transcriptional signature of lymphoblastoid cell lines of BRCA1, BRCA2 and non-BRCA1/2 high risk breast cancer families

Approximately 25% of hereditary breast cancer cases are associated with a strong familial history which can be explained by mutations in BRCA1 or BRCA2 and other lower penetrance genes. The remaining high-risk families could be classified as BRCAX (non-BRCA1/2) families.

Gene expression involving alternative splicing represents a well-known mechanism regulating the expression of multiple transcripts, which could be involved in cancer development. Thus using RNA-seq methodology, the analysis of transcriptome was undertaken to potentially reveal transcripts implicated in breast cancer susceptibility and development.

RNA was extracted from immortalized lymphoblastoid cell lines of 117 women (affected and unaffected) coming from BRCA1, BRCA2 and BRCAX families. Anova analysis revealed a total of 95 transcripts corresponding to 85 different genes differentially expressed (Bonferroni corrected p-value <0.01) between those groups. Hierarchical clustering allowed distinctive subgrouping of BRCA1/2 subgroups from BRCAX individuals. We found 67 transcripts, which could discriminate BRCAX from BRCA1/BRCA2 individuals while 28 transcripts discriminate affected from unaffected BRCAX individuals.

To our knowledge, this represents the first study identifying transcripts differentially expressed in lymphoblastoid cell lines from major classes of mutation-related breast cancer subgroups, namely BRCA1, BRCA2 and BRCAX. Moreover, some transcripts could discriminate affected from unaffected BRCAX individuals, which could represent potential therapeutic targets for breast cancer treatment.

Methods to analyze big data in pharmacogenomics research

The scale and scope of pharmacogenomics research continues to expand as the cost and efficiency of molecular data generation techniques advance. These new technologies give rise to enormous opportunity for the identification of important genetic and genomic factors important for drug treatment response. With this opportunity come significant challenges. Most of these can be categorized as 'big data' issues, facing not only pharmacogenomics, but other fields in the life sciences as well. In this review, we describe some of the analysis techniques and tools being implemented for genetic/genomic discovery in pharmacogenomics.

Polymorphic variation in TPMT is the principal determinant of TPMT phenotype: A meta-analysis of three genome-wide association studies

Thiopurine-related hematotoxicity in pediatric acute lymphoblastic leukemia (ALL) and inflammatory bowel diseases has been linked to genetically defined variability in thiopurine S-methyltransferase (TPMT) activity. While gene testing of TPMT is being clinically implemented, it is unclear if additional genetic variation influences TPMT activity with consequences for thiopurine-related toxicity. To examine this possibility, we performed a genome-wide association study (GWAS) of red blood cell TPMT activity in 844 Estonian individuals and 245 pediatric ALL cases. Additionally, we correlated genome-wide genotypes to human hepatic TPMT activity in 123 samples. Only genetic variants mapping to chromosome 6, including the TPMT gene region, were significantly associated with TPMT activity (P < 5.0 × 10^-8 ) in each of the three GWAS and a joint meta-analysis of 1,212 cases (top hit P = 1.2 × 10^-72 ). This finding is consistent with TPMT genotype being the primary determinant of TPMT activity, reinforcing the rationale for genetic testing of TPMT alleles in routine clinical practice to individualize mercaptopurine dosage.

Computational discovery of transcription factors associated with drug response

This study integrates gene expression, genotype and drug response data in lymphoblastoid cell lines with transcription factor (TF)-binding sites from ENCODE (Encyclopedia of Genomic Elements) in a novel methodology that elucidates regulatory contexts associated with cytotoxicity. The method, GENMi (Gene Expression iN the Middle), postulates that single-nucleotide polymorphisms within TF-binding sites putatively modulate its regulatory activity, and the resulting variation in gene expression leads to variation in drug response. Analysis of 161 TFs and 24 treatments revealed 334 significantly associated TF-treatment pairs. Investigation of 20 selected pairs yielded literature support for 13 of these associations, often from studies where perturbation of the TF expression changes drug response. Experimental validation of significant GENMi associations in taxanes and anthracyclines across two triple-negative breast cancer cell lines corroborates our findings. The method is shown to be more sensitive than an alternative, genome-wide association study-based approach that does not use gene expression. These results demonstrate the utility of GENMi in identifying TFs that influence drug response and provide a number of candidates for further testing.

Pharmacogenomics of platinum-based chemotherapy sensitivity in NSCLC: toward precision medicine

Lung cancer is one of the leading causes of cancer-related death in the world. Platinum-based chemotherapy is the first-line treatment for non-small-cell lung cancer (NSCLC), however, the therapeutic efficiency varies remarkably among individuals. A large number of pharmacogenomics studies aimed to identify genetic variations which can be used to predict platinum response. Those studies are leading NSCLC treatment to the new era of precision medicine. In the current review, we provided a comprehensive update on the main recent findings of genetic variations which can be used to predict platinum sensitivity in the NSCLC patients.

Clinical Evaluation of Cisplatin Sensitivity of Germline Polymorphisms in Neoadjuvant Chemotherapy for Urothelial Cancer

BACKGROUND

Level 1 evidence has demonstrated increased overall survival with cisplatin-based neoadjuvant chemotherapy for patients with muscle-invasive urothelial cancer. Usage remains low, however, in part because neoadjuvant chemotherapy will not be effective for every patient. To identify the patients most likely to benefit, we evaluated germline pharmacogenomic markers for association with neoadjuvant chemotherapy sensitivity in 2 large cohorts of patients with urothelial cancer.

PATIENTS AND METHODS

Patients receiving neoadjuvant cisplatin-based chemotherapy for muscle-invasive urothelial cancer were eligible. Nine germline single nucleotide polymorphisms (SNPs) potentially conferring platinum sensitivity were tested for an association with a complete pathologic response to neoadjuvant chemotherapy (pT0) or elimination of muscle-invasive cancer (<pT2).

RESULTS

The data from 205 patients were analyzed-59 patients were included in the discovery set and 146 in an independent replication cohort-from 3 institutions. The stage pT0 (26%) and < pT2 (50%) rates were consistent across the discovery and replication populations. Using a multivariate recessive genetic model, rs244898 in RARS (odds ratio, 6.8; 95% confidence interval, 1.8-28.9; P = .006) and rs7937567 in GALNTL4 (odds ratio, 4.8; 95% confidence interval, 1.1-22.6; P = .04) were associated with pT0 in the discovery set. Despite these large effects, neither were associated with achievement of pT0 in the replication set. A third SNP, rs10964552, was associated with stage < pT2 in the discovery set but also failed to replicate.

CONCLUSION

Germline SNPs previously associated with platinum sensitivity were not associated with the neoadjuvant chemotherapy response in a large replication cohort of patients with urothelial cancer. These results emphasize the need for replication when evaluating pharmacogenomic markers and demonstrate that multi-institutional efforts are feasible and will be necessary to achieve advances in urothelial cancer pharmacogenomics.

Genome-wide Analysis Identifies Novel Loci Associated with Ovarian Cancer Outcomes: Findings from the Ovarian Cancer Association Consortium

PURPOSE

Chemotherapy resistance remains a major challenge in the treatment of ovarian cancer. We hypothesize that germline polymorphisms might be associated with clinical outcome.

EXPERIMENTAL DESIGN

We analyzed approximately 2.8 million genotyped and imputed SNPs from the iCOGS experiment for progression-free survival (PFS) and overall survival (OS) in 2,901 European epithelial ovarian cancer (EOC) patients who underwent first-line treatment of cytoreductive surgery and chemotherapy regardless of regimen, and in a subset of 1,098 patients treated with ≥ 4 cycles of paclitaxel and carboplatin at standard doses. We evaluated the top SNPs in 4,434 EOC patients, including patients from The Cancer Genome Atlas. In addition, we conducted pathway analysis of all intragenic SNPs and tested their association with PFS and OS using gene set enrichment analysis.

RESULTS

Five SNPs were significantly associated (P ≤ 1.0 × 10(-5)) with poorer outcomes in at least one of the four analyses, three of which, rs4910232 (11p15.3), rs2549714 (16q23), and rs6674079 (1q22), were located in long noncoding RNAs (lncRNAs) RP11-179A10.1, RP11-314O13.1, and RP11-284F21.8, respectively (P ≤ 7.1 × 10(-6)). ENCODE ChIP-seq data at 1q22 for normal ovary show evidence of histone modification around RP11-284F21.8, and rs6674079 is perfectly correlated with another SNP within the super-enhancer MEF2D, expression levels of which were reportedly associated with prognosis in another solid tumor. YAP1- and WWTR1 (TAZ)-stimulated gene expression and high-density lipoprotein (HDL)-mediated lipid transport pathways were associated with PFS and OS, respectively, in the cohort who had standard chemotherapy (pGSEA ≤ 6 × 10(-3)).

CONCLUSIONS

We have identified SNPs in three lncRNAs that might be important targets for novel EOC therapies.

In vitro human cell line models to predict clinical response to anticancer drugs

In vitro human cell line models have been widely used for cancer pharmacogenomic studies to predict clinical response, to help generate pharmacogenomic hypothesis for further testing, and to help identify novel mechanisms associated with variation in drug response. Among cell line model systems, immortalized cell lines such as Epstein-Barr virus (EBV)-transformed lymphoblastoid cell lines (LCLs) have been used most often to test the effect of germline genetic variation on drug efficacy and toxicity. Another model, especially in cancer research, uses cancer cell lines such as the NCI-60 panel. These models have been used mainly to determine the effect of somatic alterations on response to anticancer therapy. Even though these cell line model systems are very useful for initial screening, results from integrated analyses of multiple omics data and drug response phenotypes using cell line model systems still need to be confirmed by functional validation and mechanistic studies, as well as validation studies using clinical samples. Future models might include the use of patient-specific inducible pluripotent stem cells and the incorporation of 3D culture which could further optimize in vitro cell line models to improve their predictive validity.

Genome-wide analysis correlates Ayurveda Prakriti

The practice of Ayurveda, the traditional medicine of India, is based on the concept of three major constitutional types (Vata, Pitta and Kapha) defined as “Prakriti”. To the best of our knowledge, no study has convincingly correlated genomic variations with the classification of Prakriti. In the present study, we performed genome-wide SNP (single nucleotide polymorphism) analysis (Affymetrix, 6.0) of 262 well-classified male individuals (after screening 3416 subjects) belonging to three Prakritis. We found 52 SNPs (p ≤ 1 × 10⁻⁵) were significantly different between Prakritis, without any confounding effect of stratification, after 10⁶ permutations. Principal component analysis (PCA) of these SNPs classified 262 individuals into their respective groups (Vata, Pitta and Kapha) irrespective of their ancestry, which represent its power in categorization. We further validated our finding with 297 Indian population samples with known ancestry. Subsequently, we found that PGM1 correlates with phenotype of Pitta as described in the ancient text of Caraka Samhita, suggesting that the phenotypic classification of India’s traditional medicine has a genetic basis; and its Prakriti-based practice in vogue for many centuries resonates with personalized medicine.

Common variants in ACYP2 influence susceptibility to cisplatin-induced hearing loss

Taking a genome-wide association study approach, we identified inherited genetic variations in ACYP2 associated with cisplatin-related ototoxicity (rs1872328: P = 3.9 × 10(-8), hazard ratio = 4.5) in 238 children with newly diagnosed brain tumors, with independent replication in 68 similarly treated children. The ACYP2 risk variant strongly predisposed these patients to precipitous hearing loss and was related to ototoxicity severity. These results point to new biology underlying the ototoxic effects of platinum agents.

Association of an inherited genetic variant with vincristine-related peripheral neuropathy in children with acute lymphoblastic leukemia

IMPORTANCE

With cure rates of childhood acute lymphoblastic leukemia (ALL) exceeding 85%, there is a need to mitigate treatment toxicities that can compromise quality of life, including peripheral neuropathy from vincristine treatment.

OBJECTIVE

To identify genetic germline variants associated with the occurrence or severity of vincristine-induced peripheral neuropathy in children with ALL.

DESIGN, SETTING, AND PARTICIPANTS

Genome-wide association study of patients in 1 of 2 prospective clinical trials for childhood ALL that included treatment with 36 to 39 doses of vincristine. Genome-wide single-nucleotide polymorphism (SNP) analysis and vincristine-induced peripheral neuropathy were assessed in 321 patients from whom DNA was available: 222 patients (median age, 6.0 years; range, 0.1-18.8 years) enrolled in 1994-1998 in the St Jude Children's Research Hospital protocol Total XIIIB with toxic effects follow-up through January 2001, and 99 patients (median age, 11.4 years; range, 3.0-23.8 years) enrolled in 2007-2010 in the Children's Oncology Group (COG) protocol AALL0433 with toxic effects follow-up through May 2011. Human leukemia cells and induced pluripotent stem cell neurons were used to assess the effects of lower CEP72 expression on vincristine sensitivity.

EXPOSURE

Treatment with vincristine at a dose of 1.5 or 2.0 mg/m2.

MAIN OUTCOMES AND MEASURES

Vincristine-induced peripheral neuropathy was assessed at clinic visits using National Cancer Institute criteria and prospectively graded as mild (grade 1), moderate (grade 2), serious/disabling (grade 3), or life threatening (grade 4).

RESULTS

Grade 2 to 4 vincristine-induced neuropathy during continuation therapy occurred in 28.8% of patients (64/222) in the St Jude cohort and in 22.2% (22/99) in the COG cohort. A SNP in the promoter region of the CEP72 gene, which encodes a centrosomal protein involved in microtubule formation, had a significant association with vincristine neuropathy (meta-analysis P = 6.3×10(-9)). This SNP had a minor allele frequency of 37% (235/642), with 50 of 321 patients (16%; 95% CI, 11.6%-19.5%) homozygous for the risk allele (TT at rs924607). Among patients with the high-risk CEP72 genotype (TT at rs924607), 28 of 50 (56%; 95% CI, 41.2%-70.0%) developed at least 1 episode of grade 2 to 4 neuropathy, a higher rate than in patients with the CEP72 CC or CT genotypes (58/271 patients [21.4%; 95% CI, 16.9%-26.7%]; P = 2.4×10(-6)). The severity of neuropathy was greater in patients homozygous for the TT genotype compared with patients with the CC or CT genotype (2.4-fold by Poisson regression [P<.0001] and 2.7-fold based on mean grade of neuropathy: 1.23 [95% CI, 0.74-1.72] vs 0.45 [95% CI, 0.3-0.6]; P = .004 by t test). Reducing CEP72 expression in human neurons and leukemia cells increased their sensitivity to vincristine.

CONCLUSIONS AND RELEVANCE

In this preliminary study of children with ALL, an inherited polymorphism in the promoter region of CEP72 was associated with increased risk and severity of vincristine-related peripheral neuropathy. If replicated in additional populations, this finding may provide a basis for safer dosing of this widely prescribed anticancer agent.

Measuring circadian clock function in human cells

Circadian clocks are present in most cells and are essential for maintenance of daily rhythms in physiology, mood, and cognition. Thus, not only neurons of the central circadian pacemaker but also many other peripheral tissues possess the same functional and self-sustained circadian clocks. Surprisingly, however, their properties vary widely within the human population. In recent years, this clock variance has been studied extensively both in health and in disease using robust lentivirus-based reporter technologies to probe circadian function in human peripheral cells as proxies for those in neurologically and physiologically relevant but inaccessible tissues. The same procedures can be used to investigate other conserved signal transduction cascades affecting multiple aspects of human physiology, behavior, and disease. Accessing gene expression variation within human populations via these powerful in vitro cell-based technologies could provide important insights into basic phenotypic diversity or to better interpret patterns of gene expression variation in disease.

Pharmacogenetic predictor of extrapyramidal symptoms induced by antipsychotics: multilocus interaction in the mTOR pathway

Antipsychotic (AP) treatment-emergent extrapyramidal symptoms (EPS) are acute adverse reactions of APs. The aim of the present study is to analyze gene-gene interactions in nine genes related to the mTOR pathway, in order to develop genetic predictors of the appearance of EPS. 243 subjects (78 presenting EPS: 165 not) from three cohorts participated in the present study: Cohort 1, patients treated with risperidone, (n=114); Cohort 2, patients treated with APs other than risperidone (n=102); Cohort 3, AP-naïve patients with first-episode psychosis treated with risperidone, paliperidone or amisulpride, n=27. We analyzed gene-gene interactions by multifactor dimensionality reduction assay (MDR). In Cohort 1, we identified a four-way interaction, including the rs1130214 (AKT1), rs456998 (FCHSD1), rs7211818 (Raptor) and rs1053639 (DDIT4), that correctly predicted 97 of the 114 patients (85% accuracy). We validated the predictive power of the four-way interaction in Cohort 2 and in Cohort 3 with 86% and 88% accuracy respectively. We develop and validate a powerful pharmacogenetic predictor of AP-induced EPS. For the first time, the mTOR pathway has been related to EPS susceptibility and AP response. However, validation in larger and independent populations will be necessary for optimal generalization.

Prediction of individual response to anticancer therapy: historical and future perspectives

Since the introduction of chemotherapy for cancer treatment in the early 20th century considerable efforts have been made to maximize drug efficiency and at the same time minimize side effects. As there is a great interpatient variability in response to chemotherapy, the development of predictive biomarkers is an ambitious aim for the rapidly growing research area of personalized molecular medicine. The individual prediction of response will improve treatment and thus increase survival and life quality of patients. In the past, cell cultures were used as in vitro models to predict in vivo response to chemotherapy. Several in vitro chemosensitivity assays served as tools to measure miscellaneous endpoints such as DNA damage, apoptosis and cytotoxicity or growth inhibition. Twenty years ago, the development of high-throughput technologies, e.g. cDNA microarrays enabled a more detailed analysis of drug responses. Thousands of genes were screened and expression levels were correlated to drug responses. In addition, mutation analysis became more and more important for the prediction of therapeutic success. Today, as research enters the area of -omics technologies, identification of signaling pathways is a tool to understand molecular mechanism underlying drug resistance. Combining new tissue models, e.g. 3D organoid cultures with modern technologies for biomarker discovery will offer new opportunities to identify new drug targets and in parallel predict individual responses to anticancer therapy. In this review, we present different currently used chemosensitivity assays including 2D and 3D cell culture models and several -omics approaches for the discovery of predictive biomarkers. Furthermore, we discuss the potential of these assays and biomarkers to predict the clinical outcome of individual patients and future perspectives.

Heat shock protein-mediated protection against Cisplatin-induced hair cell death

Cisplatin is a highly successful and widely used chemotherapy for the treatment of various solid malignancies in both adult and pediatric patients. Side effects of cisplatin treatment include nephrotoxicity and ototoxicity. Cisplatin ototoxicity results from damage to and death of cells in the inner ear, including sensory hair cells. We showed previously that heat shock inhibits cisplatin-induced hair cell death in whole-organ cultures of utricles from adult mice. Since heat shock protein 70 (HSP70) is the most upregulated HSP in response to heat shock, we investigated the role of HSP70 as a potential protectant against cisplatin-induced hair cell death. Our data using utricles from HSP70 (-/-) mice indicate that HSP70 is necessary for the protective effect of heat shock against cisplatin-induced hair cell death. In addition, constitutive expression of inducible HSP70 offered modest protection against cisplatin-induced hair cell death. We also examined a second heat-inducible protein, heme oxygenase-1 (HO-1, also called HSP32). HO-1 is an enzyme responsible for the catabolism of free heme. We previously showed that induction of HO-1 using cobalt protoporphyrin IX (CoPPIX) inhibits aminoglycoside-induced hair cell death. Here, we show that HO-1 also offers significant protection against cisplatin-induced hair cell death. HO-1 induction occurred primarily in resident macrophages, with no detectable expression in hair cells or supporting cells. Depletion of macrophages from utricles abolished the protective effect of HO-1 induction. Together, our data indicate that HSP induction protects against cisplatin-induced hair cell death, and they suggest that resident macrophages mediate the protective effect of HO-1 induction.

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.

Expression quantitative trait analyses to identify causal genetic variants for type 2 diabetes susceptibility

Type 2 diabetes (T2D) is a common metabolic disorder which is caused by multiple genetic perturbations affecting different biological pathways. Identifying genetic factors modulating the susceptibility of this complex heterogeneous metabolic phenotype in different ethnic and racial groups remains challenging. Despite recent success, the functional role of the T2D susceptibility variants implicated by genome-wide association studies (GWAS) remains largely unknown. Genetic dissection of transcript abundance or expression quantitative trait (eQTL) analysis unravels the genomic architecture of regulatory variants. Availability of eQTL information from tissues relevant for glucose homeostasis in humans opens a new avenue to prioritize GWAS-implicated variants that may be involved in triggering a causal chain of events leading to T2D. In this article, we review the progress made in the field of eQTL research and knowledge gained from those studies in understanding transcription regulatory mechanisms in human subjects. We highlight several novel approaches that can integrate eQTL analysis with multiple layers of biological information to identify ethnic-specific causal variants and gene-environment interactions relevant to T2D pathogenesis. Finally, we discuss how the eQTL analysis mediated search for “missing heritability” may lead us to novel biological and molecular mechanisms involved in susceptibility to T2D.

Integrative "omic" analysis for tamoxifen sensitivity through cell based models

It has long been observed that tamoxifen sensitivity varies among breast cancer patients. Further, ethnic differences of tamoxifen therapy between Caucasian and African American have also been reported. Since most studies have been focused on Caucasian people, we sought to comprehensively evaluate genetic variants related to tamoxifen therapy in African-derived samples. An integrative "omic" approach developed by our group was used to investigate relationships among endoxifen (an active metabolite of tamoxifen) sensitivity, SNP genotype, mRNA and microRNA expressions in 58 HapMap YRI lymphoblastoid cell lines. We identified 50 SNPs that associate with cellular sensitivity to endoxifen through their effects on 34 genes and 30 microRNA expression. Some of these findings are shared in both Caucasian and African samples, while others are unique in the African samples. Among gene/microRNA that were identified in both ethnic groups, the expression of TRAF1 is also correlated with tamoxifen sensitivity in a collection of 44 breast cancer cell lines. Further, knock-down TRAF1 and over-expression of hsa-let-7i confirmed the roles of hsa-let-7i and TRAF1 in increasing tamoxifen sensitivity in the ZR-75-1 breast cancer cell line. Our integrative omic analysis facilitated the discovery of pharmacogenomic biomarkers that potentially affect tamoxifen sensitivity.

Protein quantitative trait loci identify novel candidates modulating cellular response to chemotherapy

Annotating and interpreting the results of genome-wide association studies (GWAS) remains challenging. Assigning function to genetic variants as expression quantitative trait loci is an expanding and useful approach, but focuses exclusively on mRNA rather than protein levels. Many variants remain without annotation. To address this problem, we measured the steady state abundance of 441 human signaling and transcription factor proteins from 68 Yoruba HapMap lymphoblastoid cell lines to identify novel relationships between inter-individual protein levels, genetic variants, and sensitivity to chemotherapeutic agents. Proteins were measured using micro-western and reverse phase protein arrays from three independent cell line thaws to permit mixed effect modeling of protein biological replicates. We observed enrichment of protein quantitative trait loci (pQTLs) for cellular sensitivity to two commonly used chemotherapeutics: cisplatin and paclitaxel. We functionally validated the target protein of a genome-wide significant trans-pQTL for its relevance in paclitaxel-induced apoptosis. GWAS overlap results of drug-induced apoptosis and cytotoxicity for paclitaxel and cisplatin revealed unique SNPs associated with the pharmacologic traits (at p<0.001). Interestingly, GWAS SNPs from various regions of the genome implicated the same target protein (p<0.0001) that correlated with drug induced cytotoxicity or apoptosis (p≤0.05). Two genes were functionally validated for association with drug response using siRNA: SMC1A with cisplatin response and ZNF569 with paclitaxel response. This work allows pharmacogenomic discovery to progress from the transcriptome to the proteome and offers potential for identification of new therapeutic targets. This approach, linking targeted proteomic data to variation in pharmacologic response, can be generalized to other studies evaluating genotype-phenotype relationships and provide insight into chemotherapeutic mechanisms.

The central dogma of biology explains that DNA is transcribed to mRNA that is further translated into protein. Many genome-wide studies have implicated genetic variation that influences gene expression and that ultimately affect downstream complex traits including response to drugs. However, because of technical limitations, few studies have evaluated the contribution of genetic variation on protein expression and ensuing effects on downstream phenotypes. To overcome this challenge, we used a novel technology to simultaneously measure the baseline expression of 441 proteins in lymphoblastoid cell lines and compared them with publicly available genetic data. To further illustrate the utility of this approach, we compared protein-level measurements with chemotherapeutic induced apoptosis and cell-growth inhibition data. This study demonstrates the importance of using protein information to understand the functional consequences of genetic variants identified in genome-wide association studies. This protein data set will also have broad utility for understanding the relationship between other genome-wide studies of complex traits.

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.

In vitro and in vivo mouse models for pharmacogenetic studies

The identification of causative genes underlying biomedically relevant phenotypes, particularly complex multigenic traits, is of vital interest to modern medicine. Using genome-wide association analysis, many studies have successfully identified thousands of loci (called quantitative trait loci or QTL), some of these associating with drug response phenotypes. However, the determination and validation of putative genes has been much more challenging. The actions of drugs, both efficacious and deleterious, are complex phenotypes that are controlled or influenced in part by genetic mechanisms.Investigation for genetic correlates of complex traits and pharmacogenetic traits is often difficult to perform in human studies due to cost, availability of relevant sample population, and limited ability to control for environmental effects. These challenges can be circumvented with the use of mouse models for pharmacogenetic studies. In addition, the mouse can be treated at sub- and supratherapeutic doses and subjected to invasive procedures, which can facilitate measures of drug response phenotypes, making identification of pharmacogenetically relevant genes more feasible. The availability of multiple mouse genetic and phenotypic resources is an additional benefit to using the mouse for pharmacogenetic studies.Here, we describe the contribution of animal models, specifically the mouse, towards the field of pharmacogenetics. In this chapter, we describe different mouse models, including the knockout mouse, recombinant mouse inbred strains, in vitro mouse cell-based assays, as well as novel experimental approaches like the Collaborative Cross recombinant mouse inbred panel, which can be applied to preclinical pharmacogenetics research. These approaches can be used to assess drug response phenotypes that are difficult to model in humans, thereby facilitating drug discovery, development, and application.

Discovery and functional assessment of gene variants in the vascular endothelial growth factor pathway

Angiogenesis is a host-mediated mechanism in disease pathophysiology. The vascular endothelial growth factor (VEGF) pathway is a major determinant of angiogenesis, and a comprehensive annotation of the functional variation in this pathway is essential to understand the genetic basis of angiogenesis-related diseases. We assessed the allelic heterogeneity of gene expression, population specificity of cis expression quantitative trait loci (eQTLs), and eQTL function in luciferase assays in CEU and Yoruba people of Ibadan, Nigeria (YRI) HapMap lymphoblastoid cell lines in 23 resequenced genes. Among 356 cis-eQTLs, 155 and 174 were unique to CEU and YRI, respectively, and 27 were shared between CEU and YRI. Two cis-eQTLs provided mechanistic evidence for two genome-wide association study findings. Five eQTLs were tested for function in luciferase assays and the effect of two KRAS variants was concordant with the eQTL effect. Two eQTLs found in each of PRKCE, PIK3C2A, and MAP2K6 could predict 44%, 37%, and 45% of the variance in gene expression, respectively. This is the first analysis focusing on the pattern of functional genetic variation of the VEGF pathway genes in CEU and YRI populations and providing mechanistic evidence for genetic association studies of diseases for which angiogenesis plays a pathophysiologic role.

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.

Genome-wide discovery of genetic variants affecting tamoxifen sensitivity and their clinical and functional validation

BACKGROUND

Beyond estrogen receptor (ER), there are no validated predictors for tamoxifen (TAM) efficacy and toxicity. We utilized a genome-wide cell-based model to comprehensively evaluate genetic variants for their contribution to cellular sensitivity to TAM.

DESIGN

Our discovery model incorporates multidimensional datasets, including genome-wide genotype, gene expression, and endoxifen-induced cellular growth inhibition in the International HapMap lymphoblastoid cell lines (LCLs). Genome-wide findings were further evaluated in NCI60 cancer cell lines. Gene knock-down experiments were performed in four breast cancer cell lines. Genetic variants identified in the cell-based model were examined in 245 Caucasian breast cancer patients who underwent TAM treatment.

RESULTS

We identified seven novel single-nucleotide polymorphisms (SNPs) associated with endoxifen sensitivity through the expression of 10 genes using the genome-wide integrative analysis. All 10 genes identified in LCLs were associated with TAM sensitivity in NCI60 cancer cell lines, including USP7. USP7 knock-down resulted in increasing resistance to TAM in four breast cancer cell lines tested, which is consistent with the finding in LCLs and in the NCI60 cells. Furthermore, we identified SNPs that were associated with TAM-induced toxicities in breast cancer patients, after adjusting for other clinical factors.

CONCLUSION

Our work demonstrates the utility of a cell-based model in genome-wide identification of pharmacogenomic markers.

The role of inherited TPMT and COMT genetic variation in cisplatin-induced ototoxicity in children with cancer

Ototoxicity is a debilitating side effect of platinating agents with substantial interpatient variability. We sought to evaluate the association of thiopurine S-methyltransferase (TPMT) and catechol O-methyltransferase (COMT) genetic variations with cisplatin-related hearing damage in the context of frontline pediatric cancer treatment protocols. In 213 children from the St. Jude Medulloblastoma-96 and -03 protocols, hearing loss was related to younger age (P = 0.013) and craniospinal irradiation (P = 0.001), but did not differ by TPMT or COMT variants. Results were similar in an independent cohort of 41 children from solid-tumor frontline protocols. Functional hearing loss or hair cell damage was not different in TPMT knockout vs. wild-type mice following cisplatin treatment, and neither TPMT nor COMT variant was associated with cisplatin cytotoxicity in lymphoblastoid cell lines. In conclusion, our results indicated that TPMT or COMT genetic variation was not related to cisplatin ototoxicity in children with cancer and did not influence cisplatin-induced hearing damage in laboratory models.

Integrative genomics: quantifying significance of phenotype-genotype relationships from multiple sources of high-throughput data

Given recent advances in the generation of high-throughput data such as whole-genome genetic variation and transcriptome expression, it is critical to come up with novel methods to integrate these heterogeneous datasets and to assess the significance of identified phenotype-genotype relationships. Recent studies show that genome-wide association findings are likely to fall in loci with gene regulatory effects such as expression quantitative trait loci (eQTLs), demonstrating the utility of such integrative approaches. When genotype and gene expression data are available on the same individuals, we and others developed methods wherein top phenotype-associated genetic variants are prioritized if they are associated, as eQTLs, with gene expression traits that are themselves associated with the phenotype. Yet there has been no method to determine an overall p-value for the findings that arise specifically from the integrative nature of the approach. We propose a computationally feasible permutation method that accounts for the assimilative nature of the method and the correlation structure among gene expression traits and among genotypes. We apply the method to data from a study of cellular sensitivity to etoposide, one of the most widely used chemotherapeutic drugs. To our knowledge, this study is the first statistically sound quantification of the overall significance of the genotype-phenotype relationships resulting from applying an integrative approach. This method can be easily extended to cases in which gene expression data are replaced by other molecular phenotypes of interest, e.g., microRNA or proteomic data. This study has important implications for studies seeking to expand on genetic association studies by the use of omics data. Finally, we provide an R code to compute the empirical false discovery rate when p-values for the observed and simulated phenotypes are available.

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.

ATHENA: a tool for meta-dimensional analysis applied to genotypes and gene expression data to predict HDL cholesterol levels

Technology is driving the field of human genetics research with advances in techniques to generate high-throughput data that interrogate various levels of biological regulation. With this massive amount of data comes the important task of using powerful bioinformatics techniques to sift through the noise to find true signals that predict various human traits. A popular analytical method thus far has been the genome-wide association study (GWAS), which assesses the association of single nucleotide polymorphisms (SNPs) with the trait of interest. Unfortunately, GWAS has not been able to explain a substantial proportion of the estimated heritability for most complex traits. Due to the inherently complex nature of biology, this phenomenon could be a factor of the simplistic study design. A more powerful analysis may be a systems biology approach that integrates different types of data, or a meta-dimensional analysis. For this study we used the Analysis Tool for Heritable and Environmental Network Associations (ATHENA) to integrate high-throughput SNPs and gene expression variables (EVs) to predict high-density lipoprotein cholesterol (HDL-C) levels. We generated multivariable models that consisted of SNPs only, EVs only, and SNPs + EVs with testing r-squared values of 0.16, 0.11, and 0.18, respectively. Additionally, using just the SNPs and EVs from the best models, we generated a model with a testing r-squared of 0.32. A linear regression model with the same variables resulted in an adjusted r-squared of 0.23. With this systems biology approach, we were able to integrate different types of high-throughput data to generate meta-dimensional models that are predictive for the HDL-C in our data set. Additionally, our modeling method was able to capture more of the HDL-C variation than a linear regression model that included the same variables.