ABC transporter polymorphisms are associated with irinotecan pharmacokinetics and neutropenia.

Neutropenia is a common dose-limiting toxicity associated with irinotecan treatment. Although UGT1A1 variants have been associated with neutropenia, a fraction of neutropenia risk remains unaccounted for. To identify additional genetic markers contributing to variability in irinotecan pharmacokinetics and neutropenia, a regression analysis was performed in 78 irinotecan-treated patients to analyze comprehensively three hepatic efflux transporter genes (ABCB1, ABCC1 and ABCG2). rs6498588 (ABCC1) and rs12720066 (ABCB1) were associated with increased SN-38 exposure, and rs17501331 (ABCC1) and rs12720066 were associated with lower absolute neutrophil count nadir. rs6498588 and a variant in high linkage disequilibrium are located in transcriptionally active regions or are predicted to alter transcription factor binding sites. While enhancer activity was not evident in vitro for genomic regions containing these single-nucleotide polymorphisms, rs6498588 was significantly associated with ABCC1 expression in human liver. These results suggest that genetic variation in ABCC1 and ABCB1 may contribute to irinotecan-induced neutropenia by altering expression of transporters involved in irinotecan metabolite disposition.

Transcriptomic variation of pharmacogenes in multiple human tissues and lymphoblastoid cell lines.

Variation in the expression level and activity of genes involved in drug disposition and action ('pharmacogenes') can affect drug response and toxicity, especially when in tissues of pharmacological importance. Previous studies have relied primarily on microarrays to understand gene expression differences, or have focused on a single tissue or small number of samples. The goal of this study was to use RNA-sequencing (RNA-seq) to determine the expression levels and alternative splicing of 389 Pharmacogenomics Research Network pharmacogenes across four tissues (liver, kidney, heart and adipose) and lymphoblastoid cell lines, which are used widely in pharmacogenomics studies. Analysis of RNA-seq data from 139 different individuals across the 5 tissues (20-45 individuals per tissue type) revealed substantial variation in both expression levels and splicing across samples and tissue types. Comparison with GTEx data yielded a consistent picture. This in-depth exploration also revealed 183 splicing events in pharmacogenes that were previously not annotated. Overall, this study serves as a rich resource for the research community to inform biomarker and drug discovery and use.

Genetic advances uncover mechanisms of chemotherapy-induced peripheral neuropathy.

Chemotherapy-induced peripheral neuropathy (CIPN) is a common dose-limiting toxicity experienced in 30-40% of patients undergoing treatment with various chemotherapeutics, including taxanes, vinca alkaloids, epothilones, proteasome inhibitors, and thalidomide. Importantly, CIPN significantly affects a patient's quality of life. Recent genetic association studies are enhancing our understanding of CIPN pathophysiology and serve as a foundation for identification of genetic biomarkers to predict toxicity risk and for the development of novel strategies for prevention and treatment.

Metabolomic and Genome-wide Association Studies Reveal Potential Endogenous Biomarkers for OATP1B1.

Transporter-mediated drug-drug interactions (DDIs) are a major cause of drug toxicities. Using published genome-wide association studies (GWAS) of the human metabolome, we identified 20 metabolites associated with genetic variants in organic anion transporter, OATP1B1 (P < 5 × 10-8 ). Of these, 12 metabolites were significantly higher in plasma samples from volunteers dosed with the OATP1B1 inhibitor, cyclosporine (CSA) vs. placebo (q-value < 0.2). Conjugated bile acids and fatty acid dicarboxylates were among the metabolites discovered using both GWAS and CSA administration. In vitro studies confirmed tetradecanedioate (TDA) and hexadecanedioate (HDA) were novel substrates of OATP1B1 as well as OAT1 and OAT3. This study highlights the use of multiple datasets for the discovery of endogenous metabolites that represent potential in vivo biomarkers for transporter-mediated DDIs. Future studies are needed to determine whether these metabolites can serve as qualified biomarkers for organic anion transporters. Quantitative relationships between metabolite levels and modulation of transporters should be established.

Clinical Pharmacogenetics Implementation Consortium Guidelines for HLA-B Genotype and Abacavir Dosing: 2014 update.

The Clinical Pharmacogenetics Implementation Consortium (CPIC) Guidelines for HLA-B Genotype and Abacavir Dosing were originally published in April 2012. We reviewed recent literature and concluded that none of the evidence would change the therapeutic recommendations in the original guideline; therefore, the original publication remains clinically current. However, we have updated the Supplementary Material online and included additional resources for applying CPIC guidelines to the electronic health record. Up-to-date information can be found at PharmGKB (http://www.pharmgkb.org).

Polygenic inheritance of paclitaxel-induced sensory peripheral neuropathy driven by axon outgrowth gene sets in CALGB 40101 (Alliance).

Peripheral neuropathy is a common dose-limiting toxicity for patients treated with paclitaxel. For most individuals, there are no known risk factors that predispose patients to the adverse event, and pathogenesis for paclitaxel-induced peripheral neuropathy is unknown. Determining whether there is a heritable component to paclitaxel-induced peripheral neuropathy would be valuable in guiding clinical decisions and may provide insight into treatment of and mechanisms for the toxicity. Using genotype and patient information from the paclitaxel arm of CALGB 40101 (Alliance), a phase III clinical trial evaluating adjuvant therapies for breast cancer in women, we estimated the variance in maximum grade and dose at first instance of sensory peripheral neuropathy. Our results suggest that paclitaxel-induced neuropathy has a heritable component, driven in part by genes involved in axon outgrowth. Disruption of axon outgrowth may be one of the mechanisms by which paclitaxel treatment results in sensory peripheral neuropathy in susceptible patients.

Association of CYP2C9*2 with bosentan-induced liver injury.

Bosentan (Tracleer) is an endothelin receptor antagonist prescribed for the treatment of pulmonary arterial hypertension (PAH). Its use is limited by drug-induced liver injury (DILI). To identify genetic markers of DILI, association analyses were performed on 56 Caucasian PAH patients receiving bosentan. Twelve functional polymorphisms in five genes (ABCB11, ABCC2, CYP2C9, SLCO1B1, and SLCO1B3) implicated in bosentan pharmacokinetics were tested for associations with alanine aminotransferase (ALT), aspartate aminotransferase (AST), and DILI. After adjusting for body mass index, CYP2C9*2 was the only polymorphism associated with ALT, AST, and DILI (ß = 2.16, P = 0.024; ß = 1.92, P = 0.016; odds ratio 95% CI = 2.29-∞, P = 0.003, respectively). Bosentan metabolism by CYP2C9*2 in vitro was significantly reduced compared with CYP2C9*1 and was comparable to that by CYP2C9*3. These results suggest that CYP2C9*2 is a potential genetic marker for prediction of bosentan-induced liver injury and warrants investigation for the optimization of bosentan treatment.

Functional characterization of ABCC2 promoter polymorphisms and allele-specific expression.

Multidrug resistance protein 2 (MRP2, ABCC2) is an efflux membrane transporter highly expressed in liver, kidney and intestine with important physiological and pharmacological roles. The goal of this study was to investigate the functional significance of promoter region polymorphisms in ABCC2 and potential allele-specific expression. Twelve polymorphisms in the 1.6 kb region upstream of the translation start site were identified by resequencing 247 DNA samples from ethnically diverse individuals. Luciferase reporter gene assays showed that ABCC2 -24C>T both alone and as part of a common haplotype (-24C>T/-1019A>G/-1549G>A) increased promoter function 35% compared with the reference sequence (P<0.0001). No other common variants or haplotypes affected ABCC2 promoter activity. Allele-specific expression was also investigated as a mechanism to explain reported associations of the synonymous ABCC2 3972C>T variant with pharmacokinetic phenotypes. In Caucasian liver samples (n=41) heterozygous for the 3972C>T polymorphism, the 3972C allele was preferentially transcribed relative to the 3972T allele (P<0.0001). This allelic imbalance was particularly apparent in samples with haplotypes containing two or three promoter/untranslated region variants (-1549G>A, -1019A>G and -24C>T). The observed allelic imbalance was not associated with hepatic or renal ABCC2 mRNA expression. Additional mechanisms will need to be explored to account for the interindividual variation in ABCC2 expression and MRP2 function.

Clinical pharmacogenetics implementation consortium guidelines for HLA-B genotype and abacavir dosing.

Human leukocyte antigen B (HLA-B) is responsible for presenting peptides to immune cells and plays a critical role in normal immune recognition of pathogens. A variant allele, HLA-B*57:01, is associated with increased risk of a hypersensitivity reaction to the anti-HIV drug abacavir. In the absence of genetic prescreening, hypersensitivity affects ~6% of patients and can be life-threatening with repeated dosing. We provide recommendations (updated periodically at http://www.pharmkgb.org) for the use of abacavir based on HLA-B genotype.

Discovery of regulatory elements in human ATP-binding cassette transporters through expression quantitative trait mapping.

ATP-binding cassette (ABC) membrane transporters determine the disposition of many drugs, metabolites and endogenous compounds. Coding region variation in ABC transporters is the cause of many genetic disorders, but much less is known about the genetic basis and functional outcome of ABC transporter expression level variation. We used genotype and mRNA transcript level data from human lymphoblastoid cell lines to assess population and gender differences in ABC transporter expression, and to guide the discovery of genomic regions involved in transcriptional regulation. Nineteen of 49 ABC genes were differentially expressed between individuals of African, Asian and European descent, suggesting an important influence of race on expression level of ABC transporters. Twenty-four significant associations were found between transporter transcript levels and proximally located genetic variants. Several of the associations were experimentally validated in reporter assays. Through influencing ABC expression levels, these single-nucleotide polymorphisms may affect disease susceptibility and response to drugs.

Functional characterization of liver enhancers that regulate drug-associated transporters.

Little is known about how genetic variations in enhancers influence drug response. In this study, we investigated whether nucleotide variations in enhancers that regulate drug transporters can alter their expression levels. Using comparative genomics and liver-specific transcription factor binding site (TFBS) analyses, we identified evolutionary conserved regions (ECRs) surrounding nine liver membrane transporters that interact with commonly used pharmaceuticals. The top 50 ECRs were screened for enhancer activity in vivo, of which five--located around ABCB11, SLC10A1, SLCO1B1, SLCO1A2, and SLC47A1--exhibited significant enhancer activity. Common variants identified in a large ethnically diverse cohort (n = 272) were assayed for differential enhancer activity, and three variants were found to have significant effects on reporter activity as compared with the reference allele. In addition, one variant was associated with reduced SLCO1A2 mRNA expression levels in human liver tissues, and another was associated with increased methotrexate (MTX) clearance in patients. This work provides a general model for the rapid characterization of liver enhancers and identifies associations between enhancer variants and drug response.

The pharmacogenomics of membrane transporters project: research at the interface of genomics and transporter pharmacology.

Since the cloning of the first membrane transporter, our understanding of the role of transporters in clinical drug disposition and response has grown enormously. In parallel, large-scale genome-wide variation studies and the emerging field of pharmacogenomics have ushered in a new understanding of variations in drug response. At the crossroads of pharmacogenomics and transporter biology is the National Institutes of Health-funded Pharmacogenomics of Membrane Transporters (PMT) project, centered at the University of California, San Francisco.

The pharmacogenetics research network: from SNP discovery to clinical drug response.

The NIH Pharmacogenetics Research Network (PGRN) is a collaborative group of investigators with a wide range of research interests, but all attempting to correlate drug response with genetic variation. Several research groups concentrate on drugs used to treat specific medical disorders (asthma, depression, cardiovascular disease, addiction of nicotine, and cancer), whereas others are focused on specific groups of proteins that interact with drugs (membrane transporters and phase II drug-metabolizing enzymes). The diverse scientific information is stored and annotated in a publicly accessible knowledge base, the Pharmacogenetics and Pharmacogenomics Knowledge base (PharmGKB). This report highlights selected achievements and scientific approaches as well as hypotheses about future directions of each of the groups within the PGRN. Seven major topics are included: informatics (PharmGKB), cardiovascular, pulmonary, addiction, cancer, transport, and metabolism.

Amodiaquine metabolism is impaired by common polymorphisms in CYP2C8: implications for malaria treatment in Africa.

Metabolism of the antimalarial drug amodiaquine (AQ) into its primary metabolite, N-desethylamodiaquine, is mediated by CYP2C8. We studied the frequency of CYP2C8 variants in 275 malaria-infected patients in Burkina Faso, the metabolism of AQ by CYP2C8 variants, and the impact of other drugs on AQ metabolism. The allele frequencies of CYP2C8*2 and CYP2C8*3 were 0.155 and 0.003, respectively. No evidence was seen for influence of CYP2C8 genotype on AQ efficacy or toxicity, but sample size limited these assessments. The variant most common in Africans, CYP2C8(*)2, showed defective metabolism of AQ (threefold higher K(m) and sixfold lower intrinsic clearance), and CYP2C8(*)3 had markedly decreased activity. Considering drugs likely to be coadministered with AQ, the antiretroviral drugs efavirenz, saquinavir, lopinavir, and tipranavir were potent CYP2C8 inhibitors at clinically relevant concentrations. Variable CYP2C8 activity owing to genetic variation and drug interactions may have important clinical implications for the efficacy and toxicity of AQ.

ABCB1 pharmacogenetics: progress, pitfalls, and promise.

In 1976, Juliano and Ling(1) reported expression of a 170 kDa protein in colchicine-resistant Chinese hamster ovary (CHO) cells that was absent in drug-sensitive cells. Because this protein altered cellular permeability to colchicine, the authors named it P-glycoprotein (P-gp).(1) P-gp overexpression was described in tumor samples and leukemic cells.(2) High homology with bacterial transporters suggested that P-gp was an efflux transporter, modulating intracellular xenobiotic concentrations.(3) In 1986, the gene encoding P-gp was discovered and designated MDR1 (HUGO name ABCB1).(4) Immunohistochemical studies demonstrated P-gp expression in tissues with secretory or excretory functions (liver, kidney, and gastrointestinal tract) and at blood-tissue barrier sites, such as the blood-brain barrier.(5) This pattern of expression indicated that P-gp may influence xenobiotic response and toxicity, either through pharmacokinetic or pharmacodynamic effects.(6)

Irinotecan-induced diarrhea: functional significance of the polymorphic ABCC2 transporter protein.

Interindividual pharmacokinetic variability of the anticancer agent irinotecan is high. Life-threatening diarrhea is observed in up to 25% of patients receiving irinotecan and has been related with irinotecan pharmacokinetics and UGT1A1 genotype status. Here, we explore the association of ABCC2 (MRP2) polymorphisms and haplotypes with irinotecan disposition and diarrhea. A cohort of 167 Caucasian cancer patients who were previously assessed for irinotecan pharmacokinetics (90-min infusion given every 21 days), toxicity, and UGT1A1*28 genotype were genotyped for polymorphisms in ABCC2 using Pyrosequencing. Fifteen ABCC2 haplotypes were identified in the studied patients. The haplotype ABCC2*2 was associated with lower irinotecan clearance (28.3 versus 31.6 l/h; P=0.020). In patients who did not carry a UGT1A1*28 allele, a significant reduction of severe diarrhea was noted in patients with the ABCC2*2 haplotype (10 versus 44%; odds ratio, 0.15; 95% confidence interval, 0.04-0.61; P=0.005). This effect was not observed in patients with at least one UGT1A1*28 allele (32 versus 20%; odds ratio, 1.87; 95% confidence interval, 0.49-7.05; P=0.354). This study suggests that the presence of the ABCC2*2 haplotype is associated with less irinotecan-related diarrhea, maybe as a consequence of reduced hepatobiliary secretion of irinotecan. As the association was seen in patients not genetically predisposed at risk for diarrhea due to UGT1A1*28, confirmatory studies of the relationships of ABCC2 genotypes and irinotecan disposition and toxicity are warranted.

Differential renal gene expression in prehypertensive and hypertensive spontaneously hypertensive rats.

Development of hypertension stems from both environmental and genetic factors wherein the kidney plays a central role. Spontaneously hypertensive rats (SHR) and the nonhypertensive Wistar-Kyoto (WKY) controls are widely used as a model for studying hypertension. The present study examined the renal gene expression profiles between SHR and WKY at a prehypertensive stage (3 wk of age) and hypertensive stage (9 wk of age). Additionally, age-related changes in gene expression patterns were examined from 3 to 9 wk in both WKY and SHR. Five to six individual kidney samples of the same experimental group were pooled together, and quadruplicate hybridizations were performed using the National Institute of Environmental Health Sciences Rat version 2.0 Chip, which contains approximately 6,700 genes. Twenty two genes were found to be differentially expressed between SHR and WKY at 3 wk of age, and 104 genes were differentially expressed at 9 wk of age. Soluble epoxide hydrolase (Ephx2) was found to be significantly upregulated in SHR at both time points and was the predominant outlier. Conversely, elastase 1 (Ela1) was found to be the predominant gene downregulated in SHR at both time points. Analysis of profiles at 3 vs. 9 wk of age identified 508 differentially expressed genes in WKY rats. In contrast, only 211 genes were found to be differentially expressed during this time period in SHR. The altered gene expression patterns observed in the age-related analysis suggested significant differences in the vascular extracellular matrix system between SHR and WKY kidney. Together, our data highlight the complexity of hypertension and the numerous genes involved in and affected by this condition.