PT  - JOURNAL ARTICLE
AU  - DV Vavoulis
AU  - AT Pagnamenta
AU  - SJL Knight
AU  - MM Pentony
AU  - M Armstrong
AU  - EC Galizia
AU  - S Balestrini
AU  - SM Sisodiya
AU  - JC Taylor
TI  - Whole genome sequencing identifies putative associations between genomic polymorphisms and clinical response to the antiepileptic drug <em>levetiracetam</em>
AID  - 10.1101/19007674
DP  - 2019 Jan 01
TA  - medRxiv
PG  - 19007674
4099  - http://medrxiv.org/content/early/2019/09/27/19007674.short
4100  - http://medrxiv.org/content/early/2019/09/27/19007674.full
AB  - In the context of pharmacogenomics, whole genome sequencing provides a powerful approach for identifying correlations between response variability to specific drugs and genomic polymorphisms in a population, in an unbiased manner. In this study, we employed whole genome sequencing of DNA samples from patients showing extreme response (n=72) and non-response (n=27) to the antiepileptic drug levetiracetam, in order to identify genomic variants that underlie response to the drug. Although no common SNP (MAF>5%) crossed the conventional genome-wide significance threshold of 5×10−8, we found common polymorphisms in genes SPNS3, HDC, MDGA2, NSG1 and RASGEF1C, which collectively predict clinical response to levetiracetam in our cohort with ∼91% predictive accuracy (∼94% positive predictive value, ∼85% negative predictive value). Among these genes, HDC, NSG1, MDGA2 and RASGEF1C are potentially implicated in synaptic neurotransmission, while SPNS3 is an atypical solute carrier transporter homologous to SV2A, the known molecular target of levetiracetam. Furthermore, we performed gene- and pathway-based statistical analysis on sets of rare and low-frequency variants (MAF<5%) and we identified associations between genes or pathways and response to levetiracetam. Our findings include a) the genes PRKCB and DLG2, which are involved in glutamatergic neurotransmission, a known target of anticonvulsants, including levetiracetam; b) the genes FILIP1 and SEMA6D, which are involved in axon guidance and modelling of neural connections; and c) pathways with a role in synaptic neurotransmission, such as WNT5A-dependent internalization of FZD4 and disinhibition of SNARE formation. Targeted analysis of genes involved in neurotransmitter release and transport further supports the possibility of association between drug response and genes NSG1 and DLG2. In summary, our approach to utilise whole genome sequencing on subjects with extreme response phenotypes is a feasible route to generate plausible hypotheses for investigating the genetic factors underlying drug response variability in cases of pharmaco-resistant epilepsy.AUTHOR SUMMARY Levetiracetam (LEV) is a prominent antiepileptic drug prescribed for the treatment of both focal and generalised epilepsy. The molecular mechanism mediating its action is not well understood, but it involves the modulation of synaptic neurotransmition through binding to the synaptic vesicle glycoprotein SV2A. Identifying genomic polymorphisms that predict response to the drug is important, because it can help clinicians prescribe the most appropriate treatment in a patient-specific manner. In this study, we employed whole genome sequencing (WGS) of DNA samples from extreme responders or non-responders to LEV and we identified a small group of common variants, which successfully predict response to the drug in our cohort. These variants are mostly located in genes implicated in synaptic function. Furthermore, we identified significant associations between clinical response to LEV and low-frequency variants in genes and pathways involved in excitatory neurotransmission or in the moulding of neural networks in the brain. Our approach to utilise WGS on subjects with extreme response phenotypes is a feasible route to generate plausible hypotheses on the genomic basis of pharmaco-resistant epilepsy. We expect that the rapidly decreasing cost of WGS will allow conducting similar studies on a larger scale in the near future.Competing Interest StatementThe authors have declared no competing interest.Funding StatementThe research was supported by the National Institute for Health Research (NIHR) Oxford Biomedical Research Centre based at Oxford University Hospitals NHS Trust and University of Oxford. This work was supported by Epilepsy Society, UK. UCB Pharma (Brussels, Belgium) funded the sequencing and bioinformatics work and provided input into the analytical approaches used. Part of this work was undertaken at University College London Hospitals, which received a proportion of funding from the NIHR Biomedical Research Centres funding scheme. S Balestrini was supported by the Muir Maxwell Trust. The views expressed are those of the author(s) and not necessarily those of the NHS, the NIHR or the Department of Health. This publication presents independent research commissioned by the Health Innovation Challenge Fund (R6-388 / WT 100127), a parallel funding partnership between the Wellcome Trust and the Department of Health. The views expressed in this publication are those of the authors and not necessarily those of the Wellcome Trust or the Department of Health.Author DeclarationsAll relevant ethical guidelines have been followed and any necessary IRB and/or ethics committee approvals have been obtained.YesAll necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived.YesAny clinical trials involved have been registered with an ICMJE-approved registry such as ClinicalTrials.gov and the trial ID is included in the manuscript.Not ApplicableI have followed all appropriate research reporting guidelines and uploaded the relevant Equator, ICMJE or other checklist(s) as supplementary files, if applicable.Not ApplicableAll data supporting the results presented are included in the manuscript and supporting material