Abstract
The acceptability of antidepressant drugs partly depends on genetic factors. The list of genes involved in antidepressant response, including Adverse Drug Reactions (ADRs) is broad and contains both drug-metabolizing enzymes (pharmacogenes) and genes involved in pharmacodynamics. Variants in pharmacogenes are traditionally reported in the form of star alleles and are partially annotated with known phenotypic consequences. As it is unfeasible to analyze all genotype-phenotype pairs, computational approaches remain the practical solution. A pharmacogenetic framework to predict responses to antidepressant drug treatment would provide great benefit to patients. In this study, we present a scoring system (PharmGScore) to assess both rare and common genetic variant burden across multiple genes. The PharmGScore is constructed by normalizing and aggregating existing, well-established computational variant predictors (CADD, Fathmm-xf, PROVEAN, Mutation Assessor). We show that this score effectively distinguishes no and decreased function from normal and increased function pharmacogenetic variants reported in PharmVar (PharmGScore AUC = 0.86). PharmGScore has improved performance when compared to its component scores (AUCs: CADD = 0.79; FATHMM-XF = 0.81; PROVEAN = 0.81; Mutation Assessor = 0.75). We then apply the PharmGScore to the 200k exome sequences of the UK Biobank (UKB). We report the overrepresentation of UKB participants with high (>50) gene PharmGScore for CYP2C19 and CYP2C9 and with high (>100) compound PharmGScore from nine pharmacogenes within a group with an antidepressant toxicity diagnostic code (T43.2). We then analyze all UKB participants that received any antidepressant toxicity or ADR diagnosis (n = 602). We indicate genes for which a higher burden may be associated with antidepressant toxicity or ADRs and confirm the known roles of CYP2C19 and CYP2D6 in this process. Finally, we show that patients who experienced ADRs to antidepressants in the therapeutic process or accidental poisoning with antidepressants have a higher PharmGScore composed of nine cytochrome P450 genes. Our study proposes a novel paradigm to assess the compound genetic variant burden associated with antidepressant response from exome sequencing data. This approach can be further applied to a user-defined set of genes to investigate other pharmacological traits.
Competing Interest Statement
The authors have declared no competing interest.
Funding Statement
The research in this manuscript was funded by the National Science Center grant SONATINA UMO-2021/40/C/NZ2/00218 awarded to MB and the statutory funds of the Maj Institute of Pharmacology Polish Academy of Sciences. The computational part of this research was supported by PLGrid Infrastructure.
Author Declarations
I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.
Yes
The details of the IRB/oversight body that provided approval or exemption for the research described are given below:
The UK Biobank study was conducted under generic approval from the NHS National Research Ethics Service (approval letter dated 17 June 2011, Ref 11/NW/0382).
I confirm that all necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived, and that any patient/participant/sample identifiers included were not known to anyone (e.g., hospital staff, patients or participants themselves) outside the research group so cannot be used to identify individuals.
Yes
I understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).
Yes
I have followed all appropriate research reporting guidelines, such as any relevant EQUATOR Network research reporting checklist(s) and other pertinent material, if applicable.
Yes
Data Availability
All the code used in this project is available in a GitHub repository (https://github.com/ippas/ifpan-pharmgscore-manuscript). The star allele definitions were downloaded from the Pharmacogene Variation Consortium (PharmVar v 5.1.6) at www.pharmvar.org (Gaedigk et al. 2021). The component scores of the PharmGScore are publicly available: CADD (Rentzsch et al. 2019), FathmmXF (Rogers et al. 2018), PROVEAN (Choi and Chan 2015) and MutationAssessor (Reva, Antipin, and Sander 2011). The UK Biobank data used in this research is available to approved researchers and has been accessed through the application 62979 Impact of pharmacogenetic profiles on depression treatment outcomes. The pre-computed PharmGScore for all possible single-nucleotide variants (SNV) in WES is available for download as tsv: https://pharmgscore.labpgx.com/ (Supplementary Table A)
