Genome-Wide Pharmacogenomic Study on Methadone Maintenance Treatment Identifies SNP rs17180299 and Multiple Haplotypes on CYP2B6, SPON1, and GSG1L Associated with Plasma Concentrations of Methadone R- and S-enantiomers in Heroin-Dependent Patients

Abstract

Methadone maintenance treatment (MMT) is commonly used for controlling opioid dependence, preventing withdrawal symptoms, and improving the quality of life of heroin-dependent patients. A steady-state plasma concentration of methadone enantiomers, a measure of methadone metabolism, is an index of treatment response and efficacy of MMT. Although the methadone metabolism pathway has been partially revealed, no genome-wide pharmacogenomic study has been performed to identify genetic determinants and characterize genetic mechanisms for the plasma concentrations of methadone R- and S-enantiomers. This study was the first genome-wide pharmacogenomic study to identify genes associated with the plasma concentrations of methadone R- and S-enantiomers and their respective metabolites in a methadone maintenance cohort. After data quality control was ensured, a dataset of 344 heroin-dependent patients in the Han Chinese population of Taiwan who underwent MMT was analyzed. Genome-wide single-locus and haplotype-based association tests were performed to analyze four quantitative traits: the plasma concentrations of methadone R- and S-enantiomers and their respective metabolites. A significant single nucleotide polymorphism (SNP), rs17180299 (raw p = 2.24 × 10 ⁻⁸), was identified, accounting for 9.541% of the variation in the plasma concentration of the methadone R-enantiomer. In addition, 17 haplotypes were identified on SPON1, GSG1L, and CYP450 genes associated with the plasma concentration of methadone S-enantiomer. These haplotypes accounted for approximately one-fourth of the variation of the overall S-methadone plasma concentration. The association between the S-methadone plasma concentration and CYP2B6, SPON1, and GSG1L were replicated in another independent study. A gene expression experiment revealed that CYP2B6, SPON1, and GSG1L can be activated concomitantly through a constitutive androstane receptor (CAR) activation pathway. In conclusion, this study revealed new genes associated with the plasma concentration of methadone, providing insight into the genetic foundation of methadone metabolism. The results can be applied to predict treatment responses and methadone-related deaths for individualized MMTs.

Author Summary

Methadone maintenance treatment (MMT), among the most effective therapies for heroin-dependent patients, reduces craving and withdrawal symptoms, increases treatment compliance, and improves the quality of life of patients. The plasma concentration of methadone is a primary index for quantifying and determining therapy responses to MMT. This study was the first whole-genome pharmacogenomic study on MMT to locate genomic regions associated with the plasma concentration of methadone. The analysis identified a single nucleotide polymorphism (SNP) marker (rs17180299) and 17 haplotypes on the SPON1, GSG1L, and CYP450 genes, including CYP2B6 significantly associated with the plasma concentrations of methadone enantiomers. The identified genetic variations accounted for approximately 10% and 25% of the variations in plasma concentrations of methadone R- and S-enantiomers, respectively. The identified genetic variations have afforded insight into the genetic mechanism of the metabolism of MMT, and have potential to pave the way towards individualized MMTs for heroin-dependent patients.