Frequencies of pharmacogenomic alleles across biogeographic groups in a large-scale biobank.

Pharmacogenomics (PGx) is an integral part of precision medicine and contributes to the maximization of drug efficacy and reduction of adverse drug event risk. Accurate information on PGx allele frequencies improves the implementation of PGx. Nonetheless, curating such information from published allele data is time and resource intensive. The limited number of allelic variants in most studies leads to an underestimation of certain alleles. We applied the Pharmacogenomics Clinical Annotation Tool (PharmCAT) on an integrated 200K UK Biobank genetic dataset (N = 200,044). Based on PharmCAT results, we estimated PGx frequencies (alleles, diplotypes, phenotypes, and activity scores) for 17 pharmacogenes in five biogeographic groups: European, Central/South Asian, East Asian, Afro-Caribbean, and Sub-Saharan African. PGx frequencies were distinct for each biogeographic group. Even biogeographic groups with similar proportions of phenotypes were driven by different sets of dominant PGx alleles. PharmCAT also identified "no-function" alleles that were rare or seldom tested in certain groups by previous studies, e.g., SLCO1B1∗31 in the Afro-Caribbean (3.0%) and Sub-Saharan African (3.9%) groups. Estimated PGx frequencies are disseminated via the PharmGKB (The Pharmacogenomics Knowledgebase: www.pharmgkb.org). We demonstrate that genetic biobanks such as the UK Biobank are a robust resource for estimating PGx frequencies. Improving our understanding of PGx allele and phenotype frequencies provides guidance for future PGx studies and clinical genetic test panel design, and better serves individuals from wider biogeographic backgrounds.

How to Run the Pharmacogenomics Clinical Annotation Tool (PharmCAT).

Pharmacogenomics (PGx) investigates the genetic influence on drug response and is an integral part of precision medicine. While PGx testing is becoming more common in clinical practice and may be reimbursed by Medicare/Medicaid and commercial insurance, interpreting PGx testing results for clinical decision support is still a challenge. The Pharmacogenomics Clinical Annotation Tool (PharmCAT) has been designed to tackle the need for transparent, automatic interpretations of patient genetic data. PharmCAT incorporates a patient's genotypes, annotates PGx information (allele, genotype, and phenotype), and generates a report with PGx guideline recommendations from the Clinical Pharmacogenetics Implementation Consortium (CPIC) and/or the Dutch Pharmacogenetics Working Group (DPWG). PharmCAT has introduced new features in the last 2 years, including a variant call format (VCF) Preprocessor, the inclusion of DPWG guidelines, and functionalities for PGx research. For example, researchers can use the VCF Preprocessor to prepare biobank-scale data for PharmCAT. In addition, PharmCAT enables the assessment of novel partial and combination alleles that are composed of known PGx variants and can call CYP2D6 genotypes based on single and deletions in the input VCF file. This tutorial provides materials and detailed step-by-step instructions for how to use PharmCAT in a versatile way that can be tailored to users' individual needs.

An Evidence-Based Framework for Evaluating Pharmacogenomics Knowledge for Personalized Medicine.

Clinical annotations are one of the most popular resources available on the Pharmacogenomics Knowledgebase (PharmGKB). Each clinical annotation summarizes the association between variant-drug pairs, shows relevant findings from the curated literature, and is assigned a level of evidence (LOE) to indicate the strength of support for that association. Evidence from the pharmacogenomic literature is curated into PharmGKB as variant annotations, which can be used to create new clinical annotations or added to existing clinical annotations. This means that the same clinical annotation can be worked on by multiple curators over time. As more evidence is curated into PharmGKB, the task of maintaining consistency when assessing all the available evidence and assigning an LOE becomes increasingly difficult. To remedy this, a scoring system has been developed to automate LOE assignment to clinical annotations. Variant annotations are scored according to certain attributes, including study size, reported P value, and whether the variant annotation supports or fails to find an association. Clinical guidelines or US Food and Drug Administration (FDA)-approved drug labels which give variant-specific prescribing guidance are also scored. The scores of all annotations attached to a clinical annotation are summed together to give a total score for the clinical annotation, which is used to calculate an LOE. Overall, the system increases transparency, consistency, and reproducibility in LOE assignment to clinical annotations. In combination with increased standardization of how clinical annotations are written, use of this scoring system helps to ensure that PharmGKB clinical annotations continue to be a robust source of pharmacogenomic information.

Genomewide Association Study of Platelet Reactivity and Cardiovascular Response in Patients Treated With Clopidogrel: A Study by the International Clopidogrel Pharmacogenomics Consortium.

Antiplatelet response to clopidogrel shows wide variation, and poor response is correlated with adverse clinical outcomes. CYP2C19 loss-of-function alleles play an important role in this response, but account for only a small proportion of variability in response to clopidogrel. An aim of the International Clopidogrel Pharmacogenomics Consortium (ICPC) is to identify other genetic determinants of clopidogrel pharmacodynamics and clinical response. A genomewide association study (GWAS) was performed using DNA from 2,750 European ancestry individuals, using adenosine diphosphate-induced platelet reactivity and major cardiovascular and cerebrovascular events as outcome parameters. GWAS for platelet reactivity revealed a strong signal for CYP2C19*2 (P value = 1.67e-33). After correction for CYP2C19*2 no other single-nucleotide polymorphism reached genomewide significance. GWAS for a combined clinical end point of cardiovascular death, myocardial infarction, or stroke (5.0% event rate), or a combined end point of cardiovascular death or myocardial infarction (4.7% event rate) showed no significant results, although in coronary artery disease, percutaneous coronary intervention, and acute coronary syndrome subgroups, mutations in SCOS5P1, CDC42BPA, and CTRAC1 showed genomewide significance (lowest P values: 1.07e-09, 4.53e-08, and 2.60e-10, respectively). CYP2C19*2 is the strongest genetic determinant of on-clopidogrel platelet reactivity. We identified three novel associations in clinical outcome subgroups, suggestive for each of these outcomes.

Pharmacogenomics Clinical Annotation Tool (PharmCAT).

Pharmacogenomics (PGx) decision support and return of results is an active area of precision medicine. One challenge of implementing PGx is extracting genomic variants and assigning haplotypes in order to apply prescribing recommendations and information from the Clinical Pharmacogenetics Implementation Consortium (CPIC), the US Food and Drug Administration (FDA), the Pharmacogenomics Knowledgebase (PharmGKB), etc. Pharmacogenomics Clinical Annotation Tool (PharmCAT) (i) extracts variants specified in guidelines from a genetic data set derived from sequencing or genotyping technologies, (ii) infers haplotypes and diplotypes, and (iii) generates a report containing genotype/diplotype-based annotations and guideline recommendations. We describe PharmCAT and a pilot validation project comparing results for 1000 Genomes Project sequences of Coriell samples with corresponding Genetic Testing Reference Materials Coordination Program (GeT-RM) sample characterization. PharmCAT was highly concordant with the GeT-RM data. PharmCAT is available in GitHub to evaluate, test, and report results back to the community. As precision medicine becomes more prevalent, our ability to consistently, accurately, and clearly define and report PGx annotations and prescribing recommendations is critical.

Pharmacogenomic polygenic response score predicts ischaemic events and cardiovascular mortality in clopidogrel-treated patients.

AIMS: Clopidogrel is prescribed for the prevention of atherothrombotic events. While investigations have identified genetic determinants of inter-individual variability in on-treatment platelet inhibition (e.g. CYP2C19*2), evidence that these variants have clinical utility to predict major adverse cardiovascular events (CVEs) remains controversial. METHODS AND RESULTS: We assessed the impact of 31 candidate gene polymorphisms on adenosine diphosphate (ADP)-stimulated platelet reactivity in 3391 clopidogrel-treated coronary artery disease patients of the International Clopidogrel Pharmacogenomics Consortium (ICPC). The influence of these polymorphisms on CVEs was tested in 2134 ICPC patients (N = 129 events) in whom clinical event data were available. Several variants were associated with on-treatment ADP-stimulated platelet reactivity (CYP2C19*2, P = 8.8 × 10-54; CES1 G143E, P = 1.3 × 10-16; CYP2C19*17, P = 9.5 × 10-10; CYP2B6 1294 + 53 C > T, P = 3.0 × 10-4; CYP2B6 516 G > T, P = 1.0 × 10-3; CYP2C9*2, P = 1.2 × 10-3; and CYP2C9*3, P = 1.5 × 10-3). While no individual variant was associated with CVEs, generation of a pharmacogenomic polygenic response score (PgxRS) revealed that patients who carried a greater number of alleles that associated with increased on-treatment platelet reactivity were more likely to experience CVEs (ß = 0.17, SE 0.06, P = 0.01) and cardiovascular-related death (ß = 0.43, SE 0.16, P = 0.007). Patients who carried eight or more risk alleles were significantly more likely to experience CVEs [odds ratio (OR) = 1.78, 95% confidence interval (CI) 1.14-2.76, P = 0.01] and cardiovascular death (OR = 4.39, 95% CI 1.35-14.27, P = 0.01) compared to patients who carried six or fewer of these alleles. CONCLUSION: Several polymorphisms impact clopidogrel response and PgxRS is a predictor of cardiovascular outcomes. Additional investigations that identify novel determinants of clopidogrel response and validating polygenic models may facilitate future precision medicine strategies.

PGxMine: Text mining for curation of PharmGKB.

Precision medicine tailors treatment to individuals personal data including differences in their genome. The Pharmacogenomics Knowledgebase (PharmGKB) provides highly curated information on the effect of genetic variation on drug response and side effects for a wide range of drugs. PharmGKB's scientific curators triage, review and annotate a large number of papers each year but the task is challenging. We present the PGxMine resource, a text-mined resource of pharmacogenomic associations from all accessible published literature to assist in the curation of PharmGKB. We developed a supervised machine learning pipeline to extract associations between a variant (DNA and protein changes, star alleles and dbSNP identifiers) and a chemical. PGxMine covers 452 chemicals and 2,426 variants and contains 19,930 mentions of pharmacogenomic associations across 7,170 papers. An evaluation by PharmGKB curators found that 57 of the top 100 associations not found in PharmGKB led to 83 curatable papers and a further 24 associations would likely lead to curatable papers through citations. The results can be viewed at https://pgxmine.pharmgkb.org/ and code can be downloaded at https://github.com/jakelever/pgxmine.

The association of obesity and coronary artery disease genes with response to SSRIs treatment in major depression.

Selective serotonin reuptake inhibitors (SSRIs) are first-line antidepressants for the treatment of major depressive disorder (MDD). However, treatment response during an initial therapeutic trial is often poor and is difficult to predict. Heterogeneity of response to SSRIs in depressed patients is partly driven by co-occurring somatic disorders such as coronary artery disease (CAD) and obesity. CAD and obesity may also be associated with metabolic side effects of SSRIs. In this study, we assessed the association of CAD and obesity with treatment response to SSRIs in patients with MDD using a polygenic score (PGS) approach. Additionally, we performed cross-trait meta-analyses to pinpoint genetic variants underpinnings the relationship of CAD and obesity with SSRIs treatment response. First, PGSs were calculated at different p value thresholds (PT) for obesity and CAD. Next, binary logistic regression was applied to evaluate the association of the PGSs to SSRIs treatment response in a discovery sample (ISPC, N = 865), and in a replication cohort (STAR*D, N = 1,878). Finally, a cross-trait GWAS meta-analysis was performed by combining summary statistics. We show that the PGSs for CAD and obesity were inversely associated with SSRIs treatment response. At the most significant thresholds, the PGS for CAD and body mass index accounted 1.3%, and 0.8% of the observed variability in treatment response to SSRIs, respectively. In the cross-trait meta-analyses, we identified (1) 14 genetic loci (including NEGR1, CADM2, PMAIP1, PARK2) that are associated with both obesity and SSRIs treatment response; (2) five genetic loci (LINC01412, PHACTR1, CDKN2B, ATXN2, KCNE2) with effects on CAD and SSRIs treatment response. Our findings implicate that the genetic variants of CAD and obesity are linked to SSRIs treatment response in MDD. A better SSRIs treatment response might be achieved through a stratified allocation of treatment for MDD patients with a genetic risk for obesity or CAD.

Genome-wide and candidate gene approaches of clopidogrel efficacy using pharmacodynamic and clinical end points-Rationale and design of the International Clopidogrel Pharmacogenomics Consortium (ICPC).

RATIONALE: The P2Y12 receptor inhibitor clopidogrel is widely used in patients with acute coronary syndrome, percutaneous coronary intervention, or ischemic stroke. Platelet inhibition by clopidogrel shows wide interpatient variability, and high on-treatment platelet reactivity is a risk factor for atherothrombotic events, particularly in high-risk populations. CYP2C19 polymorphism plays an important role in this variability, but heritability estimates suggest that additional genetic variants remain unidentified. The aim of the International Clopidogrel Pharmacogenomics Consortium (ICPC) is to identify genetic determinants of clopidogrel pharmacodynamics and clinical response. STUDY DESIGN: Based on the data published on www.clinicaltrials.gov, clopidogrel intervention studies containing genetic and platelet function data were identified for participation. Lead investigators were invited to share DNA samples, platelet function test results, patient characteristics, and cardiovascular outcomes to perform candidate gene and genome-wide studies. RESULTS: In total, 17 study sites from 13 countries participate in the ICPC, contributing individual patient data from 8,829 patients. Available adenosine diphosphate-stimulated platelet function tests included vasodilator-stimulated phosphoprotein assay, light transmittance aggregometry, and the VerifyNow P2Y12 assay. A proof-of-principle analysis based on genotype data provided by each group showed a strong and consistent association between CYP2C19*2 and platelet reactivity (P value=5.1 × 10-40). CONCLUSION: The ICPC aims to identify new loci influencing clopidogrel efficacy by using state-of-the-art genetic approaches in a large cohort of clopidogrel-treated patients to better understand the genetic basis of on-treatment response variability.

Association of the Polygenic Scores for Personality Traits and Response to Selective Serotonin Reuptake Inhibitors in Patients with Major Depressive Disorder.

Studies reported a strong genetic correlation between the Big Five personality traits and major depressive disorder (MDD). Moreover, personality traits are thought to be associated with response to antidepressants treatment that might partly be mediated by genetic factors. In this study, we examined whether polygenic scores (PGSs) derived from the Big Five personality traits predict treatment response and remission in patients with MDD who were prescribed selective serotonin reuptake inhibitors (SSRIs). In addition, we performed meta-analyses of genome-wide association studies (GWASs) on these traits to identify genetic variants underpinning the cross-trait polygenic association. The PGS analysis was performed using data from two cohorts: the Pharmacogenomics Research Network Antidepressant Medication Pharmacogenomic Study (PGRN-AMPS, n = 529) and the International SSRI Pharmacogenomics Consortium (ISPC, n = 865). The cross-trait GWAS meta-analyses were conducted by combining GWAS summary statistics on SSRIs treatment outcome and on the personality traits. The results showed that the PGS for openness and neuroticism were associated with SSRIs treatment outcomes at p < 0.05 across PT thresholds in both cohorts. A significant association was also found between the PGS for conscientiousness and SSRIs treatment response in the PGRN-AMPS sample. In the cross-trait GWAS meta-analyses, we identified eight loci associated with (a) SSRIs response and conscientiousness near YEATS4 gene and (b) SSRI remission and neuroticism eight loci near PRAG1, MSRA, XKR6, ELAVL2, PLXNC1, PLEKHM1, and BRUNOL4 genes. An assessment of a polygenic load for personality traits may assist in conjunction with clinical data to predict whether MDD patients might respond favorably to SSRIs.

Life beyond kinases: structure-based discovery of sorafenib as nanomolar antagonist of 5-HT receptors.

Of great interest in recent years has been computationally predicting the novel polypharmacology of drug molecules. Here, we applied an "induced-fit" protocol to improve the homology models of 5-HT(2A) receptor, and we assessed the quality of these models in retrospective virtual screening. Subsequently, we computationally screened the FDA approved drug molecules against the best induced-fit 5-HT(2A) models and chose six top scoring hits for experimental assays. Surprisingly, one well-known kinase inhibitor, sorafenib, has shown unexpected promiscuous 5-HTRs binding affinities, K(i) = 1959, 56, and 417 nM against 5-HT(2A), 5-HT(2B), and 5-HT(2C), respectively. Our preliminary SAR exploration supports the predicted binding mode and further suggests sorafenib to be a novel lead compound for 5HTR ligand discovery. Although it has been well-known that sorafenib produces anticancer effects through targeting multiple kinases, carefully designed experimental studies are desirable to fully understand whether its "off-target" 5-HTR binding activities contribute to its therapeutic efficacy or otherwise undesirable side effects.

Clinical assessment incorporating a personal genome.

BACKGROUND: The cost of genomic information has fallen steeply, but the clinical translation of genetic risk estimates remains unclear. We aimed to undertake an integrated analysis of a complete human genome in a clinical context. METHODS: We assessed a patient with a family history of vascular disease and early sudden death. Clinical assessment included analysis of this patient's full genome sequence, risk prediction for coronary artery disease, screening for causes of sudden cardiac death, and genetic counselling. Genetic analysis included the development of novel methods for the integration of whole genome and clinical risk. Disease and risk analysis focused on prediction of genetic risk of variants associated with mendelian disease, recognised drug responses, and pathogenicity for novel variants. We queried disease-specific mutation databases and pharmacogenomics databases to identify genes and mutations with known associations with disease and drug response. We estimated post-test probabilities of disease by applying likelihood ratios derived from integration of multiple common variants to age-appropriate and sex-appropriate pre-test probabilities. We also accounted for gene-environment interactions and conditionally dependent risks. FINDINGS: Analysis of 2.6 million single nucleotide polymorphisms and 752 copy number variations showed increased genetic risk for myocardial infarction, type 2 diabetes, and some cancers. We discovered rare variants in three genes that are clinically associated with sudden cardiac death-TMEM43, DSP, and MYBPC3. A variant in LPA was consistent with a family history of coronary artery disease. The patient had a heterozygous null mutation in CYP2C19 suggesting probable clopidogrel resistance, several variants associated with a positive response to lipid-lowering therapy, and variants in CYP4F2 and VKORC1 that suggest he might have a low initial dosing requirement for warfarin. Many variants of uncertain importance were reported. INTERPRETATION: Although challenges remain, our results suggest that whole-genome sequencing can yield useful and clinically relevant information for individual patients. FUNDING: National Institute of General Medical Sciences; National Heart, Lung And Blood Institute; National Human Genome Research Institute; Howard Hughes Medical Institute; National Library of Medicine, Lucile Packard Foundation for Children's Health; Hewlett Packard Foundation; Breetwor Family Foundation.

Predicting new molecular targets for known drugs.

Although drugs are intended to be selective, at least some bind to several physiological targets, explaining side effects and efficacy. Because many drug-target combinations exist, it would be useful to explore possible interactions computationally. Here we compared 3,665 US Food and Drug Administration (FDA)-approved and investigational drugs against hundreds of targets, defining each target by its ligands. Chemical similarities between drugs and ligand sets predicted thousands of unanticipated associations. Thirty were tested experimentally, including the antagonism of the beta(1) receptor by the transporter inhibitor Prozac, the inhibition of the 5-hydroxytryptamine (5-HT) transporter by the ion channel drug Vadilex, and antagonism of the histamine H(4) receptor by the enzyme inhibitor Rescriptor. Overall, 23 new drug-target associations were confirmed, five of which were potent (<100 nM). The physiological relevance of one, the drug N,N-dimethyltryptamine (DMT) on serotonergic receptors, was confirmed in a knockout mouse. The chemical similarity approach is systematic and comprehensive, and may suggest side-effects and new indications for many drugs.

CHIP deficiency decreases longevity, with accelerated aging phenotypes accompanied by altered protein quality control.

During the course of biological aging, there is a gradual accumulation of damaged proteins and a concomitant functional decline in the protein degradation system. Protein quality control is normally ensured by the coordinated actions of molecular chaperones and the protein degradation system that collectively help to maintain protein homeostasis. The carboxyl terminus of Hsp70-interacting protein (CHIP), a ubiquitin ligase/cochaperone, participates in protein quality control by targeting a broad range of chaperone substrates for proteasome degradation via the ubiquitin-proteasome system, demonstrating a broad involvement of CHIP in maintaining cytoplasmic protein quality control. In the present study, we have investigated the influence that protein quality control exerts on the aging process by using CHIP-/- mice. CHIP deficiency in mice leads to a markedly reduced life span, along with accelerated age-related pathophysiological phenotypes. These features were accompanied by indications of accelerated cellular senescence and increased indices of oxidative stress. In addition, CHIP-/- mice exhibit a deregulation of protein quality control, as indicated by elevated levels of toxic oligomer proteins and a decline in proteasome activity. Taken together, these data reveal that impaired protein quality control contributes to cellular senescence and implicates CHIP-dependent quality control mechanisms in the regulation of mammalian longevity in vivo.

The pharmacogenetics and pharmacogenomics knowledge base: accentuating the knowledge.

PharmGKB is a knowledge base that captures the relationships between drugs, diseases/phenotypes and genes involved in pharmacokinetics (PK) and pharmacodynamics (PD). This information includes literature annotations, primary data sets, PK and PD pathways, and expert-generated summaries of PK/PD relationships between drugs, diseases/phenotypes and genes. PharmGKB's website is designed to effectively disseminate knowledge to meet the needs of our users. PharmGKB currently has literature annotations documenting the relationship of over 500 drugs, 450 diseases and 600 variant genes. In order to meet the needs of whole genome studies, PharmGKB has added new functionalities, including browsing the variant display by chromosome and cytogenetic locations, allowing the user to view variants not located within a gene. We have developed new infrastructure for handling whole genome data, including increased methods for quality control and tools for comparison across other data sources, such as dbSNP, JSNP and HapMap data. PharmGKB has also added functionality to accept, store, display and query high throughput SNP array data. These changes allow us to capture more structured information on phenotypes for better cataloging and comparison of data. PharmGKB is available at www.pharmgkb.org.

CHIP activates HSF1 and confers protection against apoptosis and cellular stress.

Induction of molecular chaperones is the characteristic protective response to environmental stress, and is regulated by a transcriptional program that depends on heat shock factor 1 (HSF1), which is normally under negative regulatory control by molecular chaperones Hsp70 and Hsp90. In metazoan species, the chaperone system also provides protection against apoptosis. We demonstrate that the dual function co-chaperone/ubiquitin ligase CHIP (C-terminus of Hsp70-interacting protein) regulates activation of the stress-chaperone response through induced trimerization and transcriptional activation of HSF1, and is required for protection against stress-induced apoptosis in murine fibroblasts. The consequences of this function are demonstrated by the phenotype of mice lacking CHIP, which develop normally but are temperature-sensitive and develop apoptosis in multiple organs after environmental challenge. CHIP exerts a central and unique role in tuning the response to stress at multiple levels by regulation of protein quality control and transcriptional activation of stress response signaling.