A systematic review of GWAS identified SNPs associated with outcomes of medications for opioid use disorder

A genome-wide association, polygenic risk score and sex study on opioid use disorder treatment outcomes

Opioid use disorder continues to be a health concern with a high rate of opioid related deaths occurring worldwide. Medication Assisted Treatments (MAT) have been shown to reduce opioid withdrawal, cravings and opioid use, however variability exists in individual's treatment outcomes. Sex-specific differences have been reported in opioid use patterns, polysubstance use and health and social functioning. Candidate gene studies investigating methadone dose as an outcome have identified several candidate genes and only five genome-wide associations studies have been conducted for MAT outcomes. This study aimed to identify genetic variants associated with MAT outcomes through genome-wide association study (GWAS) and test the association between genetic variants previously associated with methadone dose through a polygenic risk score (PRS). Study outcomes include: continued opioid use, relapse, methadone dose and opioid overdose. No genome-wide significance SNPs or sex-specific results were identified. The PRS identified statistically significant results (p < 0.05) for the outcome of methadone dose (R2 = 3.45 × 10-3). No other PRS was statistically significant. This study provides evidence for association between a PRS and methadone dose. More research on the PRS to increase the variance explained is needed before it can be used as a tool to help identify a suitable methadone dose within this population.

Methylation and expression quantitative trait loci rs1799971 in the OPRM1 gene and rs4654327 in the OPRD1 gene are associated with opioid use disorder

Opioid use disorder (OUD) is a chronic and relapsing brain disease that results in significant mortality worldwide. Genetic factors are estimated to contribute to 40%-60% of the liability, with polymorphisms of opioid receptor genes implicated in this disorder. However, the mechanisms underlying these associations are not yet fully understood. In the present study, we first examined the methylation levels in the promoter region of the OPRM1, OPRD1, and OPRK1 genes in 111 healthy controls (HCs) and 120 patients with OUD, and genotyped three tag SNPs in these genes. Correlations between these SNPs and the methylation levels of the CpG sites and expression levels of the genes were analyzed. After identifying the mQTLs and eQTLs, we determined the associations between the mQTLs/eQTLs and susceptibility to and characteristics of OUD in 930 HCs and 801 patients with OUD. Our results demonstrated that SNPs rs1799971 in the OPRM1 gene and rs4654327 in the OPRD1 gene were both mQTLs and eQTLs. We observed unique correlations between mQTLs and methylation levels of several CpG sites in the OUD group compared to the HC group. Interestingly, both the two mQTLs and eQTLs were associated with the susceptibility to OUD. In conclusion, we suppose that mQTLs and eQTLs in genes may underlie the associations between certain risk genetic polymorphisms and OUD. These mQTLs and eQTLs could potentially serve as promising biomarkers for better management of opioid misuse.

Association between rs1799971 in the mu opioid receptor gene and methadone maintenance treatment response

OBJECTIVE

Genetic variations can affect individual response to methadone maintenance treatment (MMT) for heroin addiction. The A118G variant (rs1799971) in the mu opioid receptor gene (OPRM1) is a potential candidate single nucleotide polymorphism (SNP) for personalized MMT. This study determined whether rs1799971 is related to MMT response or dose.

METHODS

We recruited 286 MMT patients from a Han Chinese population. The rs1799971 genotype was determined via TaqMan genotyping assay. The genetic effect of this SNP on MMT response or dose was evaluated using logistic regression. A meta-analysis was performed to merge all available data to evaluate the role of rs1799971 in MMT using RevMan 5.3 software.

RESULTS

No statistical significance was observed in the association between the OPRM1 rs1799971 and MMT response or dose in our Chinese cohort. Meta-analysis indicated that the OPRM1 A118G variation was not significantly associated with MMT response or dose requirement.

CONCLUSION

The results suggest that rs1799971 in OPRM1 might not play a critical role alone in influencing MMT response or dose.

CpH methylome analysis in human cortical neurons identifies novel gene pathways and drug targets for opioid use disorder

INTRODUCTION

DNA methylation (DNAm), an epigenetic mechanism, has been associated with opioid use disorder (OUD) in preclinical and human studies. However, most of the studies have focused on DNAm at CpG sites. DNAm at non-CpG sites (mCpHs, where H indicates A, T, or C) has been recently shown to have a role in gene regulation and to be highly abundant in neurons. However, its role in OUD is unknown. This work aims to evaluate mCpHs in the human postmortem orbital frontal cortex (OFC) in the context of OUD.

METHODS

A total of 38 Postmortem OFC samples were obtained from the VA Brain Bank (OUD = 12; Control = 26). mCpHs were assessed using reduced representation oxidative bisulfite sequencing in neuronal nuclei. Differential analysis was performed using the "methylkit" R package. Age, ancestry, postmortem interval, PTSD, and smoking status were included as covariates. Significant mCpHs were set at q-value < 0.05. Gene Ontology (GO) and KEGG enrichment analyses were performed for the annotated genes of all differential mCpH loci using String, ShinyGO, and amiGO software. Further, all annotated genes were analyzed using the Drug gene interaction database (DGIdb).

RESULTS

A total of 2,352 differentially methylated genome-wide significant mCpHs were identified in OUD, mapping to 2,081 genes. GO analysis of genes with differential mCpH loci showed enrichment for nervous system development (p-value = 2.32E-19). KEGG enrichment analysis identified axon guidance and glutamatergic synapse (FDR 9E-4-2.1E-2). Drug interaction analysis found 3,420 interactions between the annotated genes and drugs, identifying interactions with 15 opioid-related drugs, including lofexidine and tizanidine, both previously used for the treatment of OUD-related symptoms.

CONCLUSION

Our findings suggest a role of mCpHs for OUD in cortical neurons and reveal important biological pathways and drug targets associated with the disorder.