- Primary research
- Open Access
Integrating rare genetic variants into pharmacogenetic drug response predictions
© The Author(s). 2018
- Received: 6 March 2018
- Accepted: 8 May 2018
- Published: 25 May 2018
Variability in genes implicated in drug pharmacokinetics or drug response can modulate treatment efficacy or predispose to adverse drug reactions. Besides common genetic polymorphisms, recent sequencing projects revealed a plethora of rare genetic variants in genes encoding proteins involved in drug metabolism, transport, and response.
To understand the global importance of rare pharmacogenetic gene variants, we mapped the variability in 208 pharmacogenes by analyzing exome sequencing data from 60,706 unrelated individuals and estimated the importance of rare and common genetic variants using a computational prediction framework optimized for pharmacogenetic assessments. Our analyses reveal that rare pharmacogenetic variants were strongly enriched in mutations predicted to cause functional alterations. For more than half of the pharmacogenes, rare variants account for the entire genetic variability. Each individual harbored on average a total of 40.6 putatively functional variants, rare variants accounting for 10.8% of these. Overall, the contribution of rare variants was found to be highly gene- and drug-specific. Using warfarin, simvastatin, voriconazole, olanzapine, and irinotecan as examples, we conclude that rare genetic variants likely account for a substantial part of the unexplained inter-individual differences in drug metabolism phenotypes.
Combined, our data reveal high gene and drug specificity in the contributions of rare variants. We provide a proof-of-concept on how this information can be utilized to pinpoint genes for which sequencing-based genotyping can add important information to predict drug response, which provides useful information for the design of clinical trials in drug development and the personalization of pharmacological treatment.
- Personalized medicine
- ADME genes
- Genetic variability
- Drug response
The response of patients to medical treatment is influenced by a variety of physiological, pathological, environmental, and genetic factors . These inter-individual differences can lower treatment efficacy or manifest in adverse drug reactions (ADRs), which are estimated to cause around 6.5% of all hospital admissions . Overall, the genetic makeup of a patient accounts for 20–30% of the inter-individual variability in drug response , but for certain clinically important drugs, such as metoprolol and torsemide, twin studies suggested genetic contributions to the variability in their pharmacokinetics of up to 90% .
Genetic variability in phase I and phase II enzymes, transporters, cytochrome reductases, and nuclear receptors, hereafter jointly termed pharmacogenes, can modulate drug absorption, distribution, metabolism, and excretion (ADME), thereby shaping human drug response and the risk of ADRs. Prominent examples include associations of common TPMT variants with hematological toxicity of 6-mercaptopurines, ultrarapid metabolism of CYP2D6 with codeine toxicity, and effects of specific CYP2C19 polymorphisms on the response to clopidogrel or proton pump inhibitors . Yet, a substantial fraction of the heritable variability in drug response cannot be explained by these common variants, suggesting that other genetic factors are important contributors. In recent years, increasing capacities and decreasing costs of next-generation sequencing (NGS) platforms have facilitated large-scale studies of genetic variation and NGS assays are becoming increasingly implemented in clinical diagnostics . Importantly, NGS-based analyses revealed that over 90% of the overall genetic variability in pharmacogenes is allotted to rare genetic variants, but the impact of rare genetic variability on drug pharmacokinetics has not been systematically evaluated (Additional file 1: Table S1) [7–10].
We therefore analyzed the distribution of rare and common gene variants in the 208 clinically most relevant pharmacogenes of 60,706 unrelated individuals and leveraged these genetic variability profiles to predict the relevance of rare SNVs for the pharmacokinetics or ADR risks for several clinically important drugs of diverse therapeutic areas. Based on these analyses, we conclude that the contribution of rare genetic variants is gene- and drug-specific and can account for a substantial part of the unexplained genetic inter-individual variability in drug response. Furthermore, we highlight genes for which comprehensive NGS-based genotyping instead of candidate SNP interrogations can reveal important additional information to personalize pharmacological treatment strategies. The presented data incentivizes the consideration of rare pharmacogenetic variants for the guidance of personalized drug therapy and holds important implications for the design of clinical trials.
Human sequencing data from 60,706 unrelated individuals was obtained from the Exome Aggregation Consortium (ExAC) database , a platform that provides summary frequency information of exonic genetic variants from 17 large-scale sequencing projects. Notably, consistency of the individual data sets is assured by reprocessing of all raw data through the same bioinformatic pipelines. We complemented the obtained variants with six non-exonic variants from the 1000 Genomes Project  that define CYP1A2*1C (rs2069514), CYP1A2*1F (rs762551), CYP2C19*17 (rs12248560), CYP3A4*22 (rs35599367), CYP3A5*3 (rs776746), and UGT1A1*28 (rs8175347). Variants with less than 10,000 called high-quality alleles were not considered. Novel variants were defined relative to dbSNP release 135.
Loss-of-function (LoF) intolerance scores were provided by Lek et al. based on the expectation-maximization algorithm . In brief, low scores (< 0.1) indicate that the number of protein-truncating variants is similar to what is expected by chance, whereas high scores (> 0.9) indicate much fewer of such variants are observed than would be expected, suggesting haploinsufficiency. Aggregated functional variant frequency is defined as the sum of MAFs of all variants predicted to be deleterious. Variants with MAF ≤ 0.01 were considered as rare, and variants with MAF > 0.01 were considered as common.
Variant analyses and computational functionality predictions
Computational algorithms mostly use evolutionary conservation as a metric to predict whether a given variant likely has functional effects. Importantly, we previously evaluated 18 current functionality prediction methods and found that their predictive performance was low for poorly conserved genes, such as cytochrome P450s .
Here, we therefore used a functionality prediction method that we previously developed . In brief, using high-quality experimental data for 123 pharmacogenetic alleles, Zhou et al. tailored the parameterization of 18 different algorithms specifically for ADME genes and integrated the results of multiple prediction methods into an ADME-optimized prediction framework . Finally, the model’s performance was validated in an independent validation cohort of additional 121 experimentally characterized variants. Overall, the model achieved 92% sensitivity and 95% specificity for loss-of-function and functionally neutral variants, respectively, thereby substantially outperforming previous computational tools on pharmacogenomics data sets.
Analysis of the genetic landscape in 208 human pharmacogenes
Evolutionary constraints in transporters as well as in phase 1 and phase 2 drug-metabolizing genes were low as judged by the large numbers of loss-of-function variants identified in these genes (LoF intolerance score = 0.08 ± 0.02 SEM; see the “Methods” section for details). In contrast, nuclear receptors and other selected genes with importance for drug response were highly LoF-intolerant (LoF intolerance score = 0.53 ± 0.11 SEM), comparable to values observed across haploinsufficient genes (LoF intolerance score > 0.5; Fig. 1c) . Importantly, the vast majority of variants were rare (98.5%; MAF < 1%) or very rare (96.2%; MAF < 0.1%) and more than half (51.1%) of all variants were only detected in a single individual, highlighting the genetic diversity in human pharmacogenes (Fig. 1d).
Rare variants contribute substantially to functional variability
Notably, rare pharmacogenetic variants were strongly enriched in mutations predicted to cause functional alterations (Fig. 2b), consistent with previous reports [12, 14, 15]. In the 208 pharmacogenes combined, each individual harbored on average a total of 40.6 putatively functional variants (Fig. 2c). Rare variants accounted for 4.4 (10.8%) of these functional variants, of which 1.8, 1.7, 0.7, and 0.2 were allotted to transporters, phase 1, phase 2, and other pharmacogenes, respectively (Fig. 2c).
Prediction of the importance of rare genetic variants for drug response
Given the significant contribution of rare genetic variants to the functional variability in pharmacogenes, we considered it of importance to include rare variants into predictions of drug response. Using the genetic information as template, we analyzed the contribution of rare genetic variants for drug pharmacokinetics and/or drug response, focusing on five drugs with well-characterized pharmacology and substantial unexplained inter-individual variability. Specifically, we evaluated the relevance of rare SNVs for the anticoagulant warfarin, the HMG-CoA reductase inhibitor simvastatin, the antifungal voriconazole, the antipsychotic olanzapine, and the antineoplastic agent irinotecan (Additional file 2: Table S2). We first estimated the relative importance of different genetic factors for drug metabolism phenotypes of the specific drugs based on extensive literature analysis. Subsequently in a second step, we integrated these evaluations with our genetic variability data to derive assessments of the impact of rare genetic variants on the pharmacokinetics of or response to the given drug.
Warfarin is a racemic mixture of the R- and S-stereoisomers of which the S-form is at least five times more potent. Warfarin response is influenced by common genetic polymorphisms in CYP2C9, CYP4F2, and VKORC1, which jointly explain up to 45% of warfarin dose requirements . Yet also other genes, such as CYP3A4, CYP1A2, EPHX1, and ABCB1, have been implicated in warfarin pharmacokinetics [17–19]. However, despite this extensive knowledge of warfarin transport and metabolism, around 40% of the variability in warfarin dose requirements remains unexplained by common genetic variants and other patient-specific factors .
ADRs related to high-dose simvastatin therapy are strongly linked to the common (12.9% MAF) genetic variant rs4149056 in SLCO1B1 (encoding the transporter OATP1B1) with an odds ratio of 4.5 per copy of the risk allele . Toxicity is caused by an impaired hepatic uptake of the drug that results in elevated plasma concentrations of simvastatin acid, which have been shown to cause myotoxicity in vitro  (Fig. 3d). In our analyses, we correctly predicted the functional impact of rs4149056 and did not find additional common variants with reduced functionality. However, we identified rare deleterious variants with an aggregated frequency of 1.2%, which are estimated to jointly explain 8.7% of the genetic basis of SLCO1B1 variability (Fig. 3e). Similarly, rare variants are expected to contribute substantially to the metabolism and transport of simvastatin with an aggregated rare functional variant frequency of 1.3, 1.1, and 0.8% for ABCB1, CYP3A4, and ABCG2, respectively (Fig. 3f).
Voriconazole efficacy and ADRs
Serum olanzapine levels
The therapeutic benefits for schizophrenic or bipolar patients when treated with the antipsychotic olanzapine are limited by extensive inter-individual variability in olanzapine serum concentrations, which can result in exposure levels outside the therapeutic interval . As olanzapine serum levels are directly linked to the likelihood of therapeutic success  and the risk of ADRs [31, 32], an individualization of dosing regimens promises to increase treatment success rates.
While the metabolism of olanzapine is well characterized, the influence of genetic factors is more controversial (Fig. 4d). CYP2D6 and UGT2B10 hydroxylate or glucuronidate olanzapine, respectively, but so far, no study demonstrated clinically relevant effects of haplotypes of these genes on olanzapine exposure. In contrast, multiple clinical association studies found that genetic variants in CYP1A2, FMO3, and UGT1A4 could explain up to 50% of differences in olanzapine serum levels, whereas other studies failed to replicate such associations ( and references therein). Rare variants substantially contribute to UGT1A4 and UGT2B10 variability and are predicted to account for 100 and 47% of the genetically encoded inter-individual variability in the functionality of these genes (Fig. 4e, f). On the contrary, we estimate that rare variants only explain 6.3, 5.5, and 1.6% of the variability in CYP2D6, FMO3, and CYP1A2, respectively.
Genetic variants in UGT1A1 have been reproducibly linked to neutropenia and diarrhea toxicity in various ethnicities and dosing regimens . Furthermore, multiple polymorphisms in the transporter genes ABCB1, ABCC1, ABCC2, ABCG2, and SLCO1B1 have been implicated in irinotecan clearance and/or risk of toxicity [37–40]. While associations between CYP3A genotype and irinotecan pharmacokinetics are controversial, incorporation of CYP3A activity data into dosing calculations have resulted in reduced incidence of severe neutropenia .
Interestingly, our computational analyses of population-scale sequencing data indicate that rare genetic variants are important factors for irinotecan activation and transport (Fig. 5b). The aggregated frequency of rare deleterious alleles in CES1 and CES2 were 1.5 and 0.4%, respectively, accounting for 19.6 and 100% of the functional genetic variants in these genes. Similarly, for SLCO1B1 and ABCC1, 8.7 and 40.5% of all deleterious variants were assigned to rare variants, while for ABCB1, ABCC2, and ABCG2, no common variants with functional consequences were identified. Notably though, we do not consider here variants with functional impacts that do not result in changes of the gene product, such as the synonymous variants rs1045642 (ABCB1 I1145I) and rs1128503 (ABCB1 G412G) or the UTR variant rs717620 in ABCC2. In contrast, inter-individual variability in UGT1A1 is primarily due to common polymorphisms, such as UGT1A1*28, and only an estimated 4.5% are allotted to rare SNVs.
Combined, our evaluations indicate that rare genetic variants in pharmacogenes have the potential to explain a substantial part of the unexplained genetic variability in drug metabolism phenotypes. Examples were selected for which gene-drug interactions were well studied, and we speculate that the relative importance of rare variants is even higher for less extensively characterized drugs. Furthermore, we give indications about the extent of genetically encoded functional variability that would be missed when only considering common genetic variants, thereby providing guidance for the optimal drug-specific choice of genotyping strategy.
From a drug development perspective, an appropriate pharmacokinetic profile is of key importance to achieve the desired spatial and temporal exposure pattern of a given drug of interest. However, genetic variants in ADME genes encoding for transporters, drug-metabolizing enzymes, or nuclear receptors modulate drug pharmacokinetics and thus impact treatment efficacy and the risk of ADRs. Consequently, 190 drugs approved by the US Food and Drug Administration (FDA) and 155 drugs approved by EMA currently contain pharmacogenetic information in their labels, of which many are related to drug pharmacokinetics [42, 43]. Besides well-characterized common polymorphisms, ADME genes harbor a plethora of rare genetic variants that are not interrogated by current pharmacogenomic genotyping panels. By leveraging large-scale whole-exome sequencing data from 60,706 individuals, we present here the first analysis in which we systematically integrated information about rare genetic variability into predictions of pharmacokinetic variability (Additional file 2: Table S2).
Individual in silico functionality prediction algorithms distinguish deleterious from neutral variants with sensitivities and specificities between 60 and 90%. Furthermore, by using computational methods optimized for the evaluation of pharmacogenes with low evolutionary constraints, we were able to show that up to 92% sensitivity and 95% specificity can be achieved for loss-of-function and functionally neutral variants, respectively . Furthermore, rare copy number variations in pharmacogenes, accounting for up to 1% of all loss-of-function alleles, are an additional source of genetic variability with relevance for drug metabolism phenotypes that are commonly not considered by computational functionality prediction algorithms .
Notably, genetic variability in non-coding regions has been demonstrated to have important influence on phenotypic traits . However, while promising progress has been made regarding the prediction of the effects of those variants based on DNA sequence , no current prediction framework can reliably predict the functionality of non-coding genetic variation, such as synonymous variants or variants located in UTRs, promoters, or enhancers. With respect to CYP alleles, such mutations represent < 1% of all functionally important variant alleles described (https://www.pharmvar.org/genes). We therefore restricted our analyses to the evaluation of LOF variants and variants that directly affect the amino acid sequence of their respective gene products.
Thus, while the predictive power of current functionality prediction methods is still not sufficient to support a recommendation of these tools for genetic counseling of individual patients, our data indicate however that leveraging of NGS technology can yield significant amounts of additional information for pharmacogenomic predictions on a population scale. Accordingly, we advocate for the development of a widened perspective in which conclusions about the functionality of a gene product are not solely based on the interrogation of few common variants. Rather, we recommend that the entire spectrum of genetic variability, including rare or novel variants, should be considered and integrated into gene activity scores. This holistic perspective is especially important as rare polymorphisms are enriched in variants that alter the functionality of the gene product and are the sole genetic factors for variability of more than 50% of the pharmacogenes analyzed (Fig. 2b).
Relating the functional inventory of pharmacogenetic variability to the pharmacology of selected drugs of interest can provide important insights into predicted hotspots of unexplained inter-individual differences in drug metabolism-related phenotypes. In this work, we estimated the relative contribution of rare genetic variants to the variability in pharmacokinetics and/or ADR risk of five clinically important drugs from different therapeutic areas. Rare variants are estimated to contribute significantly to the inter-individual variability of warfarin pharmacokinetics and irinotecan toxicity accounting for 18.4% of deleterious CYP2C9 alleles and > 40% of the variability in irinotecan transport (Figs. 3a and 5). In contrast, the relative importance of rare variants is expected to be lower for the metabolism of simvastatin, voriconazole, and olanzapine for which rare variants only contribute between 1.6 and 8.7% of the key metabolic and/or transport processes. Thus, we find that the relevance of rare genetic variants is highly drug-specific, depending on the gene products involved. These findings suggest that it is likely that the inter-individual variability in pharmacokinetics and response for certain drugs is to a large extent determined by rare genetic variability, which is important to consider particularly in drug development. Integrating pharmacological information of the drug of interest with information about the distribution of rare variants in pharmacogenes can guide the design of the genotyping strategy most suitable to reveal important additional genetic factors that improve the prediction of drug metabolism phenotypes.
For many drugs, genetic variability cannot be directly translated into effectiveness or ADR susceptibility due to the interdependency of different metabolic pathways. For instance, genetic or pharmacological inhibition of the major pathway of a given drug can result in a shunt to an alternative otherwise negligible metabolic route, as observed for oxycodone . To date, data regarding the effects of such gene-gene or gene-drug interactions is sparse, which complicates predictions of drug effectiveness or safety even when rare genetic variants are incorporated into the analyses. Thus, while the consideration of pharmacogenomic information including rare genetic variants promises to improve ADME predictions, further work particularly in physiologically based pharmacokinetic (PBPK) modeling is necessary to reliably predict treatment outcomes for the individual patient.
Currently, genotyping is largely based on the interrogation of well-characterized, common polymorphisms. This strategy neglects the impacts of rare variants as experimental in vitro or in vivo data that demonstrate their functional impact is not available. Yet, due to rapidly decreasing sequencing times and costs, we suggest that current NGS technology in combination with more advanced computational prediction methods could already today facilitate the refinement of individualized predictions regarding drug efficacy and its propensity to cause ADRs and thereby to contribute to the implementation of pharmacogenetic markers into routine care [48, 49]. Furthermore, the multitude of ongoing sequencing projects on unprecedented scale, such as the 100K Genomes Project run by the British Department of Health and the 1 Million Genomes project as part of the Chinese Precision Medicine Initiative, will soon provide a wealth of information about non-transcribed, regulatory regions and improve linkage information between variants, which will allow to expand the scope of computational analyses from variants to haplotypes. For the field of pharmacogenomics, these developments hold promise to further increase the predictive power of gene-drug response predictions and to allow more accurate estimates of drug response across more narrowly stratified subpopulations.
We present results suggesting that integration of rare pharmacogenomic variability can improve predictions of drug pharmacokinetics compared to the use of candidate variants. This information is important for drug development and clinical care as well as for future preemptive pharmacogenomic advice. Furthermore, these data incentivize the design of prospective trials using NGS-based genotyping for specific medications, such as warfarin and irinotecan, to assess whether clinical outcomes can be improved.
The authors would like to thank the Exome Aggregation Consortium and the groups that provided exome variant data for comparison.
This study was supported by the European Union’s Horizon 2020 research and innovation program U-PGx (grant agreement no. 668353), by the Swedish Research Council (grant agreement numbers 2015-02760, 2016-01153, and 2016-01154), and by the Lennart Philipson Foundation.
Availability of data and materials
The data analyzed in this study is publically available under the ODC Open Database License through the server of the Exome Aggregation Consortium.
MI-S and VML designed the study. YZ collected the exome sequencing data. VML devised and conducted the data analysis. SM, MI-S, and VML wrote the manuscript. All authors discussed and agreed on the final version of the manuscript.
Ethics approval and consent to participate
In this study, only publically available, anonymized data released under the Fort Lauderdale Agreement was analyzed, and thus, separate ethical approval was not required.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Lauschke VM, Ingelman-Sundberg M. The importance of patient-specific factors for hepatic drug response and toxicity. Int J Mol Sci. 2016;17(10):1714–27. https://doi.org/10.3390/ijms17101714.View ArticlePubMed CentralGoogle Scholar
- Pirmohamed M, James S, Meakin S, et al. Adverse drug reactions as cause of admission to hospital: prospective analysis of 18 820 patients. BMJ. 2004;329(7456):15–9. https://doi.org/10.1136/bmj.329.7456.15.View ArticlePubMedPubMed CentralGoogle Scholar
- Sim SC, Kacevska M, Ingelman-Sundberg M. Pharmacogenomics of drug-metabolizing enzymes: a recent update on clinical implications and endogenous effects. Pharmacogenomics J. 2012;13(1):1–11. https://doi.org/10.1038/tpj.2012.45.View ArticlePubMedGoogle Scholar
- Matthaei J, Brockmöller J, Tzvetkov MV, et al. Heritability of metoprolol and torsemide pharmacokinetics. Clin Pharmacol Ther. 2015;98(6):611–21. https://doi.org/10.1002/cpt.258.View ArticlePubMedGoogle Scholar
- Lauschke VM, Milani L, Ingelman-Sundberg M. Pharmacogenomic biomarkers for improved drug therapy-recent progress and future developments. AAPS J. 2017;20(1):4. https://doi.org/10.1208/s12248-017-0161-x.View ArticlePubMedGoogle Scholar
- Singh RR, Luthra R, Routbort MJ, et al. Implementation of next generation sequencing in clinical molecular diagnostic laboratories: advantages, challenges and potential. Exp Rev Precision Med Drug Dev. 2016;1(1):109–20. https://doi.org/10.1080/23808993.2015.1120401.View ArticleGoogle Scholar
- Gordon AS, Tabor HK, Johnson AD, et al. Quantifying rare, deleterious variation in 12 human cytochrome P450 drug-metabolism genes in a large-scale exome dataset. Hum Mol Genet. 2014;23(8):1957–63. https://doi.org/10.1093/hmg/ddt588.View ArticlePubMedGoogle Scholar
- Fujikura K, Ingelman-Sundberg M, Lauschke VM. Genetic variation in the human cytochrome P450 supergene family. Pharmacogenet Genomics. 2015;25(12):584–94. https://doi.org/10.1097/FPC.0000000000000172.View ArticlePubMedGoogle Scholar
- Kozyra M, Ingelman-Sundberg M, Lauschke VM. Rare genetic variants in cellular transporters, metabolic enzymes, and nuclear receptors can be important determinants of interindividual differences in drug response. Genet Med. 2017;19(1):20–9. https://doi.org/10.1038/gim.2016.33.View ArticlePubMedGoogle Scholar
- Bush WS, Crosslin DR, Owusu-Obeng A, et al. Genetic variation among 82 pharmacogenes: the PGRNseq data from the eMERGE network. Clin Pharmacol Ther. 2016;100(2):160–9. https://doi.org/10.1002/cpt.350.View ArticlePubMedPubMed CentralGoogle Scholar
- Lek M, Karczewski KJ, Minikel EV, et al. Analysis of protein-coding genetic variation in 60,706 humans. Nature. 2016;536(7616):285–91. https://doi.org/10.1038/nature19057.View ArticlePubMedPubMed CentralGoogle Scholar
- Abecasis GR, Chakravarti A, Donnelly P, et al. A global reference for human genetic variation. Nature. 2015;526(7571):68–74. https://doi.org/10.1038/nature15393.View ArticlePubMedGoogle Scholar
- Zhou Y, Mkrtchian S, Kumondai M, et al. An optimized prediction framework to assess the functional impact of pharmacogenetic variants. Pharmacogenomics J. (In review)Google Scholar
- Nelson MR, Wegmann D, Ehm MG, et al. An abundance of rare functional variants in 202 drug target genes sequenced in 14,002 people. Science. 2012;337(6090):100–4. https://doi.org/10.1126/science.1217876.View ArticlePubMedPubMed CentralGoogle Scholar
- Tennessen JA, Bigham AW, O’Connor TD, et al. Evolution and functional impact of rare coding variation from deep sequencing of human exomes. Science. 2012;337(6090):64–9. https://doi.org/10.1126/science.1219240.View ArticlePubMedPubMed CentralGoogle Scholar
- Johnson JA, Cavallari LH. Warfarin pharmacogenetics. Trends Cardiovasc Med. 2015;25(1):33–41. https://doi.org/10.1016/j.tcm.2014.09.001.View ArticlePubMedGoogle Scholar
- Kaminsky LS, Zhang Z-Y. Human P450 metabolism of warfarin. Pharmacol Ther. 1997;73(1):67–74. https://doi.org/10.1016/s0163-7258(96)00140-4.View ArticlePubMedGoogle Scholar
- Wadelius M, Sörlin K, Wallerman O, et al. Warfarin sensitivity related to CYP2C9, CYP3A5, ABCB1 (MDR1) and other factors. Pharmacogenomics J. 2003;4(1):40–8. https://doi.org/10.1038/sj.tpj.6500220.View ArticleGoogle Scholar
- Pautas E, Moreau C, Gouin-Thibault I, et al. Genetic factors (VKORC1, CYP2C9, EPHX1, and CYP4F2) are predictor variables for warfarin response in very elderly, frail inpatients. Clin Pharmacol Ther. 2009;87(1):57–64. https://doi.org/10.1038/clpt.2009.178.View ArticlePubMedGoogle Scholar
- The International Warfarin Pharmacogenetics Consortium, Klein TE, Altman RB, et al. Estimation of the warfarin dose with clinical and pharmacogenetic data. N Engl J Med. 2009;360(8):753–64. https://doi.org/10.1056/NEJMoa0809329.View ArticlePubMed CentralGoogle Scholar
- Hodges LM, Markova SM, Chinn LW, et al. Very important pharmacogene summary: ABCB1 (MDR1, P-glycoprotein). Pharmacogenet Genomics. 2011;21(3):152–61. https://doi.org/10.1097/FPC.0b013e3283385a1c.View ArticlePubMedPubMed CentralGoogle Scholar
- SEARCH Collaborative Group, Link E, Parish S, et al. SLCO1B1 variants and statin-induced myopathy—a genomewide study. N Engl J Med. 2008;359(8):789–99. https://doi.org/10.1056/NEJMoa0801936.View ArticleGoogle Scholar
- Kobayashi M, Chisaki I, Narumi K, et al. Association between risk of myopathy and cholesterol-lowering effect: a comparison of all statins. Life Sci. 2008;82(17–18):969–75. https://doi.org/10.1016/j.lfs.2008.02.019.View ArticlePubMedGoogle Scholar
- Hyland R, Jones BC, Smith DA. Identification of the cytochrome P450 enzymes involved in the N-oxidation of voriconazole. Drug Metab Dispos. 2003;31(5):540–7.View ArticlePubMedGoogle Scholar
- Yanni SB, Annaert PP, Augustijns P, et al. Role of flavin-containing monooxygenase in oxidative metabolism of voriconazole by human liver microsomes. Drug Metab Dispos. 2008;36(6):1119–25. https://doi.org/10.1124/dmd.107.019646.View ArticlePubMedPubMed CentralGoogle Scholar
- Owusu Obeng A, Egelund EF, Alsultan A, et al. CYP2C19 polymorphisms and therapeutic drug monitoring of voriconazole: are we ready for clinical implementation of pharmacogenomics? Pharmacotherapy. 2014;34(7):703–18. https://doi.org/10.1002/phar.1400.View ArticlePubMedGoogle Scholar
- Niwa T, Hata T. The effect of genetic polymorphism on the inhibition of azole antifungal agents against CYP2C9-mediated metabolism. J Pharm Sci. 2016;105(3):1345–8. https://doi.org/10.1016/j.xphs.2016.01.007.View ArticlePubMedGoogle Scholar
- He HR, Sun JY, Ren XD, et al. Effects of CYP3A4 polymorphisms on the plasma concentration of voriconazole. Eur J Clin Microbiol Infect Dis. 2014;34(4):811–9. https://doi.org/10.1007/s10096-014-2294-5.View ArticlePubMedGoogle Scholar
- Patel MX, Bowskill S, Couchman L, et al. Plasma olanzapine in relation to prescribed dose and other factors. J Clin Psychopharmacol. 2011;31(4):411–7. https://doi.org/10.1097/JCP.0b013e318221b408.View ArticlePubMedGoogle Scholar
- Perry PJ, Lund BC, Sanger T, et al. Olanzapine plasma concentrations and clinical response: acute phase results of the North American Olanzapine Trial. J Clin Psychopharmacol. 2001;21(1):14–20.View ArticlePubMedGoogle Scholar
- Perry PJ, Argo TR, Carnahan RM, et al. The association of weight gain and olanzapine plasma concentrations. J Clin Psychopharmacol. 2005;25(3):250–4. https://doi.org/10.1097/01.jcp.0000162800.64378.82.View ArticlePubMedGoogle Scholar
- Kinon BJ, Volavka J, Stauffer V, et al. Standard and higher dose of olanzapine in patients with schizophrenia or schizoaffective disorder. J Clin Psychopharmacol. 2008;28(4):392–400. https://doi.org/10.1097/JCP.0b013e31817e63a5.View ArticlePubMedGoogle Scholar
- Söderberg MM, Dahl M-L. Pharmacogenetics of olanzapine metabolism. Pharmacogenomics. 2013;14(11):1319–36. https://doi.org/10.2217/pgs.13.120.View ArticlePubMedGoogle Scholar
- Fuchs CS, Moore MR, Harker G, et al. Phase III comparison of two irinotecan dosing regimens in second-line therapy of metastatic colorectal cancer. J Clin Oncol. 2003;21(5):807–14. https://doi.org/10.1200/JCO.2003.08.058.View ArticlePubMedGoogle Scholar
- Brandi G, Dabard J, Raibaud P, et al. Intestinal microflora and digestive toxicity of irinotecan in mice. Clin Cancer Res. 2006;12(4):1299–307. https://doi.org/10.1158/1078-0432.CCR-05-0750.View ArticlePubMedGoogle Scholar
- Campbell JM, Stephenson MD, Bateman E, et al. Irinotecan-induced toxicity pharmacogenetics: an umbrella review of systematic reviews and meta-analyses. Pharmacogenomics J. 2017;17(1):21–8. https://doi.org/10.1038/tpj.2016.58.View ArticlePubMedGoogle Scholar
- Han J-Y, Lim H-S, Yoo Y-K, et al. Associations of ABCB1, ABCC2, and ABCG2 polymorphisms with irinotecan-pharmacokinetics and clinical outcome in patients with advanced non-small cell lung cancer. Cancer. 2007;110(1):138–47. https://doi.org/10.1002/cncr.22760.View ArticlePubMedGoogle Scholar
- Han J-Y, Lim H-S, Park YH, et al. Integrated pharmacogenetic prediction of irinotecan pharmacokinetics and toxicity in patients with advanced non-small cell lung cancer. Lung Cancer. 2009;63(1):115–20. https://doi.org/10.1016/j.lungcan.2007.12.003.View ArticlePubMedGoogle Scholar
- Innocenti F, Kroetz DL, Schuetz E, et al. Comprehensive pharmacogenetic analysis of irinotecan neutropenia and pharmacokinetics. J Clin Oncol. 2009;27(16):2604–14. https://doi.org/10.1200/JCO.2008.20.6300.View ArticlePubMedPubMed CentralGoogle Scholar
- Li M, Seiser EL, Baldwin RM, et al. ABC transporter polymorphisms are associated with irinotecan pharmacokinetics and neutropenia. Pharmacogenomics J. 2018;18(1):35–42. https://doi.org/10.1038/tpj.2016.75.View ArticlePubMedGoogle Scholar
- van der Bol JM, Mathijssen RHJ, Creemers GJM, et al. A CYP3A4 phenotype-based dosing algorithm for individualized treatment of irinotecan. Clin Cancer Res. 2010;16(2):736–42. https://doi.org/10.1158/1078-0432.CCR-09-1526.View ArticlePubMedGoogle Scholar
- https://www.fda.gov/Drugs/ScienceResearch/ResearchAreas/Pharmacogenetics/ucm083378.htm. Accessed 20 May 2017.
- Ehmann F, Caneva L, Prasad K, et al. Pharmacogenomic information in drug labels: European Medicines Agency perspective. Pharmacogenomics J. 2015;15(3):201–10. https://doi.org/10.1038/tpj.2014.86.View ArticlePubMedGoogle Scholar
- Santos M, Niemi M, Hiratsuka M, et al. Novel copy-number variations in pharmacogenes contribute to interindividual differences in drug pharmacokinetics. Genet Med. 2017; https://doi.org/10.1038/gim.2017.156. [Epub ahead of print]
- Maurano MT, Humbert R, Rynes E, et al. Systematic localization of common disease-associated variation in regulatory DNA. Science. 2012;337(6099):1190–5. https://doi.org/10.1126/science.1222794.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee D, Gorkin DU, Baker M, et al. A method to predict the impact of regulatory variants from DNA sequence. Nat Genet. 2015;47(8):955–61. https://doi.org/10.1038/ng.3331.View ArticlePubMedPubMed CentralGoogle Scholar
- Samer CF, Daali Y, Wagner M, et al. The effects of CYP2D6 and CYP3A activities on the pharmacokinetics of immediate release oxycodone. Br J Pharmacol. 2010;160(4):907–18. https://doi.org/10.1111/j.1476-5381.2010.00673.x.View ArticlePubMedPubMed CentralGoogle Scholar
- Lauschke VM, Ingelman-Sundberg M. How to consider rare genetic variants in personalized drug therapy. Clin Pharmacol Ther. 2018;19:20. https://doi.org/10.1002/cpt.976.View ArticleGoogle Scholar
- Lauschke VM, Ingelman-Sundberg M. Precision medicine and rare genetic variants. Trends Pharmacol Sci. 2016;37(2):85–6. https://doi.org/10.1016/j.tips.2015.10.006.View ArticlePubMedGoogle Scholar
- Spitzer M, Wildenhain J, Rappsilber J, et al. BoxPlotR: a web tool for generation of box plots. Nat Methods. 2014;11(2):121–2. https://doi.org/10.1038/nmeth.2811.View ArticlePubMedPubMed CentralGoogle Scholar