GRK5 Is Associated with the Regulation of Methadone Dosage in Heroin Dependence

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.

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.

Common and distinguishing genetic factors for substance use behavior and disorder: an integrated analysis of genomic and transcriptomic studies from both human and animal studies

BACKGROUND AND AIMS

Genomic and transcriptomic findings greatly broaden the biological knowledge regarding substance use. However, systematic convergence and comparison evidence of genome-wide findings is lacking for substance use. Here, we combined all the genome-wide findings from both substance use behavior and disorder (SUBD) and identified common and distinguishing genetic factors for different SUBDs.

METHODS

Systemic literature search for genome-wide association (GWAS) and RNA-seq studies of alcohol/nicotine/drug use behavior (partially meets or not reported diagnostic criteria) and alcohol use behavior and disorder (AUBD), nicotine use behavior and disorder (NUBD) and drug use behavior and disorder (DUBD) was performed using PubMed and the GWAS catalog. Drug use was focused upon cannabis, opioid, cocaine and methamphetamine use. GWAS studies required case-control or case/cohort samples. RNA-seq studies were based on brain tissues. The genes which contained significant single nucleotide polymorphism (P ≤ 1 × 10^-6 ) in GWAS and reported as significant in RNA-seq studies were extracted. Pathway enrichment was performed by using Metascape. Gene interaction networks were identified by using the Protein Interaction Network Analysis database.

RESULTS

Total SUBD-related 2910 genes were extracted from 75 GWAS studies (2 773 889 participants) and 17 RNA-seq studies. By overlapping the genes and pathways of AUBD, NUBD and DUBD, four shared genes (CACNB2, GRIN2B, PLXDC2 and PKNOX2), four shared pathways [two Gene Ontology (GO) terms of 'modulation of chemical synaptic transmission', 'regulation of trans-synaptic signaling', two Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of 'dopaminergic synapse', 'cocaine addiction'] were identified (significantly higher than random, P < 1 × 10^-5 ). The top shared KEGG pathways (Benjamini-Hochberg-corrected P-value < 0.05) in the pairwise comparison of AUBD versus DUBD, NUBD versus DUBD, AUBD versus NUBD were 'Epstein-Barr virus infection', 'protein processing in endoplasmic reticulum' and 'neuroactive ligand-receptor interaction', respectively. We also identified substance-specific genetic factors: i.e. ADH1B and ALDH2 were unique for AUBD, while CHRNA3 and CHRNA4 were unique for NUBD.

CONCLUSIONS

This systematic review identifies the shared and unique genes and pathways for alcohol, nicotine and drug use behaviors and disorders at the genome-wide level and highlights critical biological processes for the common and distinguishing vulnerability of substance use behaviors and disorders.

Association of the D-amino acid oxidase gene with methadone dose in heroin dependent patients under methadone maintenance treatment

Methadone is a synthetic opioid used for the maintenance treatment (MMT) of heroin dependence. It primarily binds to the μ-opioid receptor (MOR; with its gene, namely OPRM1). Methadone is also an N-methyl-D-aspartate (NMDA) receptor antagonist. The role of NMDA receptor in the regulatory mechanisms of methadone dosage in heroin dependent patients is so far not clear. D-amino acid oxidase (DAO) is an important enzyme that indirectly activates the NMDA receptor through its effect on the D-serine level. To test the hypothesis that genetic polymorphisms in the DAO gene are associated with methadone treatment dose and responses, we selected four single nucleotide polymorphisms (SNPs) in DAO from the literature reports of the Taiwanese population. SNPs were genotyped in 344 MMT patients. In this study, we identified a functional SNP rs55944529 in the DAO gene that reveals a modest but significant association with the methadone dosage in the recessive model of analysis (P = 0.003) and plasma concentrations (P = 0.003) in MMT patients. However, it did not show association with plasma methadone concentration in multiple linear regression analysis. It is also associated with the methadone adverse reactions of dry mouth (P = 0.002), difficulty with urination (P = 0.0003) in the dominant model, and the withdrawal symptoms of yawning (P = 0.005) and gooseflesh skin (P = 0.004) in the recessive model. Our results suggest a role of the indirect regulatory mechanisms of the NMDA reporter, possibly via the DAO genetic variants, in the methadone dose and some adverse reactions in MMT patients.

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.

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

BACKGROUND

Patients with opioid use disorder (OUD) display an interindividual variability in their response to medications for opioid use disorder (MOUD). A genetic basis may explain the variability in this response. However, no consensus has been reached regarding which genetic variants significantly contribute to MOUD outcomes.

OBJECTIVES

This systematic review aims to summarize genome-wide significant findings on MOUD outcomes and critically appraise the quality of the studies involved.

METHODS

Databases searched from inception until August 21st, 2020 include: MEDLINE, Web of Science, EMBASE, CINAHL and Pre-CINAHL, GWAS Catalog and GWAS Central. The included studies had to be GWASs that assessed MOUD in an OUD population. All studies were screened in duplicate. The quality of the included studies was scored and assessed using the Q-Genie tool. Quantitative analysis, as planned in the protocol, was not feasible, so the studies were analyzed qualitatively.

RESULTS

Our search identified 7292 studies. Five studies meeting the eligibility criteria were included. However, only three studies reported results that met our significance threshold of p ≤ 1.0 × 10-7. In total, 43 genetic variants were identified. Variants corresponding to CNIH3 were reported to be associated with daily heroin injection in Europeans, OPRM1, TRIB2, and ZNF146 with methadone dose in African Americans, EYS with methadone dose in Europeans, and SPON1 and intergenic regions in chromosomes 9 and 3 with plasma concentrations of S-methadone, R-methadone, and R-EDDP, respectively, in Han Chinese.

LIMITATIONS

The limitations of this study include not being able to synthesize the data in a quantitative way and a conservative eligibility and data collection model.

CONCLUSION

The results from this systematic review will aid in highlighting significant genetic variants that can be replicated in future OUD pharmacogenetics research to ascertain their role in patient-specific MOUD outcomes. Systematic review registration number CRD42020169121.

Add-On Selective Estrogen Receptor Modulators for Methadone Maintenance Treatment

Methadone maintenance treatment (MMT) remains the cornerstone for the management of opiate abuse. However, MMT can be associated with complex factors, including complications during the tolerance phase, the inability of some patients to maintain treatment effects during the tapering or abstinence phases, and the development of methadone dependence. Previous studies have revealed a sex disparity in MMT efficacy, showing that women undergoing MMT experiencing an increase in psychological symptoms compared with men and suggesting a link between disparate responses and the effects of estrogen signaling on methadone metabolism. More specifically, estradiol levels are positively associated with MMT dosing, and the expression of a single-nucleotide polymorphism (SNP) associated with estrogen receptor (ER) regulation is also associated with MMT dosing. In addition to performing mechanistic dissections of estrogen signaling in the presence of methadone, past studies have also proposed the targeting of estrogen signaling during MMT. The present report provides an overview of the relevant literature regarding sex effects, including differences in sex hormones and their potential impacts on MMT regimens. Moreover, this article provides a pharmacological perspective on the targeting of estrogen signals through the use of selective ER modulators (SERMs) during MMT. Preliminary preclinical experiments were also performed to evaluate the potential effects of targeting estrogen signaling with tamoxifen on methadone metabolism.

G protein-coupled receptor kinases are associated with Alzheimer's disease pathology

AIM

Alzheimer's disease (AD) is characterised by extracellular deposition of amyloid-β (Aβ) in amyloid plaques and intracellular aggregation and accumulation of hyperphosphorylated tau in neurofibrillary tangles (NFTs). Although several kinases have been identified to contribute to the pathological phosphorylation of tau, kinase-targeted therapies for AD have not been successful in clinical trials. Critically, the kinases responsible for numerous identified tau phosphorylation sites remain unknown. G protein-coupled receptor (GPCR) kinases (GRKs) have recently been implicated in phosphorylation of non-GPCR substrates, for example, tubulin and α-synuclein, and in neurological disorders, including schizophrenia and Parkinson's disease. Accordingly, we investigated the involvement of GRKs in the pathophysiology of AD.

METHODS

We performed a comprehensive immunohistochemical and biochemical analysis of the ubiquitously expressed GRKs, namely, GRK2, 3, 5 and 6, in postmortem human brain tissue of control subjects and AD patients.

RESULTS

GRKs display unique cell-type-specific expression patterns in neurons, astrocytes and microglia. Levels of GRKs 2, 5 and 6 are specifically decreased in the CA1 region of the AD hippocampus. Biochemical evidence indicates that the GRKs differentially associate with total, soluble and insoluble pools of tau in the AD brain. Complementary immunohistochemical studies indicate that the GRKs differentially colocalise with total tau, phosphorylated tau and NFTs. Notably, GRKs 3 and 5 also colocalise with amyloid plaques.

CONCLUSION

These studies establish a link between GRKs and the pathological phosphorylation and accumulation of tau and amyloid pathology in AD brains and suggest a novel role for these kinases in regulation of the pathological hallmarks of AD.

Socioeconomic Deprivation Index Is Associated With Psychiatric Disorders: An Observational and Genome-wide Gene-by-Environment Interaction Analysis in the UK Biobank Cohort

BACKGROUND

Psychiatric disorders are among the largest and fastest-growing categories of the global disease burden. However, limited effort has been made to further elucidate associations between socioeconomic factors and psychiatric disorders from a genetic perspective.

METHODS

We randomly divided 501,882 participants in the UK Biobank cohort with socioeconomic Townsend deprivation index (TDI) data into a discovery cohort and a replication cohort. For both cohorts, we first conducted regression analyses to evaluate the associations between the TDI and common psychiatric disorders or traits, including anxiety, bipolar disorder, self-harm, and depression (based on self-reported depression and Patient Health Questionnaire scores). We then performed a genome-wide gene-by-environment interaction study using PLINK 2.0 with the TDI as an environmental factor to explore interaction effects.

RESULTS

In the discovery cohort, significant associations were observed between the TDI and psychiatric disorders (p < 4.00 × 10^-16), including anxiety (odds ratio [OR] = 1.08, 95% confidence interval [CI] = 1.07-1.10), bipolar disorder (OR = 1.42, 95% CI = 1.36-1.48), self-harm (OR = 1.21, 95% CI = 1.19-1.23), self-reported depression (OR = 1.22, 95% CI = 1.20-1.24), and Patient Health Questionnaire scores (β = .07, SE = 0.004). We observed similar significant associations in the replication cohort. In addition, multiple candidate loci were identified by the genome-wide gene-by-environment interaction study, including rs10886438 at 10q26.11 (GRK5) (p = 5.72 × 10^-11) for Patient Health Questionnaire scores and rs162553 at 2p22.2 (CYP1B1) (p = 2.25 × 10^-9) for self-harm.

CONCLUSIONS

Our findings suggest the relevance of the TDI to psychiatric disorders. The genome-wide gene-by-environment interaction study identified several candidate genes interacting with the TDI, providing novel clues for understanding the biological mechanism of associations between the TDI and psychiatric disorders.

Targeting GRK5 for Treating Chronic Degenerative Diseases

G protein-coupled receptors (GPCRs) are the largest family of cell-surface receptors and they are responsible for the transduction of extracellular signals, regulating almost all aspects of mammalian physiology. These receptors are specifically regulated by a family of serine/threonine kinases, called GPCR kinases (GRKs). Given the biological role of GPCRs, it is not surprising that GRKs are also involved in several pathophysiological processes. Particular importance is emerging for GRK5, which is a multifunctional protein, expressed in different cell types, and it has been found located in single or multiple subcellular compartments. For instance, when anchored to the plasma membrane, GRK5 exerts its canonical function, regulating GPCRs. However, under certain conditions (e.g., pro-hypertrophic stimuli), GRK5 translocates to the nucleus of cells where it can interact with non-GPCR-related proteins as well as DNA itself to promote “non-canonical” signaling, including gene transcription. Importantly, due to these actions, several studies have demonstrated that GRK5 has a pivotal role in the pathogenesis of chronic-degenerative disorders. This is true in the cardiac cells, tumor cells, and neurons. For this reason, in this review article, we will inform the readers of the most recent evidence that supports the importance of targeting GRK5 to prevent the development or progression of cancer, cardiovascular, and neurological diseases.

Genetic implications in quality palliative care and preventing opioid crisis in cancer-related pain management

The prevalence of cancer-related pain is 64% among patients with metastatic, advanced, or terminal cancer, 59% among patients undergoing anticancer treatment, and 33% among patients who completed curative treatment. According to the World Health Organization cancer pain relief guidelines, opioid analgesics are the mainstay analgesic therapy in addition to conventional first-step analgesics, such as non-steroidal anti-inflammatory drugs and acetaminophen. The indications for strong opioids have recently been expanded to include mild-to-moderate pain in addition to moderate-to-severe pain. The U.S. Centers for Disease Control and Prevention guidelines emphasize that realistic expectations should be weighed against potential serious harm from opioids, rather than relying on the unrealized long-term benefits of these drugs. Therefore, treatment strategies for both cancer-related chronic or acute pain have been unfortunately deviated from opioid analgesics. The barriers hindering adequate cancer-related pain management with opioid analgesics are related to the inadequate knowledge of opioid analgesics (e.g., effective dose, adverse effects, and likelihood of addiction or tolerance). To achieve adequate opioid availability, these barriers should be overcome in a clinically suitable manner. Genetic assessments could play an important role in overcoming challenges in opioid management. To balance the improvement in opioid availability and the prevention of opioid misuse and addiction, the following two considerations concerning opioids and genetic polymorphisms warrant attention: (A) pain severity, opioid sensitivity, and opioid tolerance; and (B) vulnerability to opioid dependence and addiction.

Endogenous Opiates and Behavior: 2018

This paper is the forty-first consecutive installment of the annual anthological review of research concerning the endogenous opioid system, summarizing articles published during 2018 that studied the behavioral effects of molecular, pharmacological and genetic manipulation of opioid peptides and receptors as well as effects of opioid/opiate agonists and antagonists. The review is subdivided into the following specific topics: molecular-biochemical effects and neurochemical localization studies of endogenous opioids and their receptors (2), the roles of these opioid peptides and receptors in pain and analgesia in animals (3) and humans (4), opioid-sensitive and opioid-insensitive effects of nonopioid analgesics (5), opioid peptide and receptor involvement in tolerance and dependence (6), stress and social status (7), learning and memory (8), eating and drinking (9), drug abuse and alcohol (10), sexual activity and hormones, pregnancy, development and endocrinology (11), mental illness and mood (12), seizures and neurologic disorders (13), electrical-related activity and neurophysiology (14), general activity and locomotion (15), gastrointestinal, renal and hepatic functions (16), cardiovascular responses (17), respiration and thermoregulation (18), and immunological responses (19).

Genetic variants in NECTIN4 encoding an adhesion molecule are associated with continued opioid use

Methadone is a synthetic opioid used as maintenance treatment for patients addicted to heroin. Skin irritation is one of the adverse events caused by opioid use. 344 methadone maintenance treatment (MMT) patients were recruited with records and measurements on methadone dose, plasma methadone concentrations, and treatment emergent symptom scales (TESS). 15 patients reported with skin irritation. Five SNPs located within the NECTIN4 genetic region were genotyped. The NECTIN4 gene within the adherens junction interaction pathway was associated with methadone dose in pathway-based genome wide association analyses (P = 0.0008). Three highly-linked SNPs, rs11265549, rs3820097, and rs4656978, were significantly associated with methadone dose (P = 0.0003), plasma concentrations of R,S-methadone (P = 0.0004) and TNF-α (P = 0.010) in all 344 MMT patients, and with self-report skin irritation symptom scores (P = 0.010) in the 15 MMT patients who reported with skin irritation. To identify the possible roles of plasma level of Nectin-4 in the responses to MMT and opioid use, additional age- and gender-matched 51 controls and 83 methadone-free abstinent former heroin users were recruited. Plasma level of Nectin-4 was the highest in MMT patients among the three groups. The results suggest involvement of genetic variants on NECTIN4 in methadone dose. Plasma Nectin-4 level is likely an indicator for continued use of opioids.

Genetic polymorphisms in the opioid receptor delta 1 (OPRD1) gene are associated with methadone dose in methadone maintenance treatment for heroin dependence

Delta opioid receptor (DOR) is well known to be involved in heroin dependence. This study tested the hypothesis that single nucleotide polymorphisms (SNPs) in the opioid receptor delta 1 (OPRD1) gene coding region are associated with treatment responses in a methadone maintenance therapy (MMT) cohort in Taiwan. Three hundred forty-four MMT patients were recruited. Diastolic/systolic blood pressure, heart rate, methadone dosage, and plasma concentrations of methadone were recorded. Twenty-five SNPs located within the OPRD1 genetic region were selected and genotyped from the genomic DNA of all 344 participants. After pairwise tagger analyses, tagger SNP rs204047 showed a significant association with methadone dosage (P = 0.0019), and tagger SNPs rs204047 and rs797397 were significantly associated with plasma R, S-methadone concentrations (P < 0.0006) in patients tested negative in the urine morphine test, which indicated patients with a better response to MMT. The major genotype carriers showed a higher methadone dosage and higher plasma concentrations of R, S-methadone than the minor genotype carriers. The results indicated that OPRD1 genetic variants were associated with methadone dosage and methadone plasma concentration in MMT patients with a negative morphine test result.

National Institute on Drug Abuse genomics consortium white paper: Coordinating efforts between human and animal addiction studies

The National Institute on Drug Abuse Genetics and Epigenetics Cross-Cutting Research Team convened a diverse group of researchers, clinicians, and healthcare providers on the campus of the University of California, San Diego, in June 2018. The goal was to develop strategies to integrate genetics and phenotypes across species to achieve a better understanding of substance use disorders through associations between genotypes and addictive behaviors. This conference (a) discussed progress in harmonizing large opioid genetics cohorts, (b) discussed phenotypes that are used for genetics studies in humans, (c) examined phenotypes that are used for genetics studies in animal models, (d) identified synergies and gaps in phenotypic analyses of human and animal models and (e) identified strategies to integrate genetics and genomics data with phenotypes across species. The meeting consisted of panels that focused on phenotype harmonization (Dr. Laura Bierut, Dr. Olivier George, Dr. Dan Larach and Dr. Sesh Mudumbai), translating genetic findings between species (Dr. Elissa Chesler, Dr. Gary Peltz and Dr. Abraham Palmer), interpreting and understanding allelic variations (Dr. Vanessa Troiani and Dr. Tamara Richards) and pathway conservation in animal models and human studies (Dr. Robert Hitzemann, Dr. Huda Akil and Dr. Laura Saba). There were also updates that were provided by large consortia (Dr. Susan Tapert, Dr. Danielle Dick, Dr. Howard Edenberg and Dr. Eric Johnson). Collectively, the conference was convened to discuss progress and changes in genome-wide association studies.

Relevance of CYP2B6 and CYP2D6 genotypes to methadone pharmacokinetics and response in the OPAL study

AIMS

Our study aimed to evaluate the impacts of the cytochrome P450 (CYP) 2B6-G516T and CYP2D6 genetic polymorphisms on pharmacokinetic and clinical parameters in patients receiving methadone maintenance treatment.

METHODS

Opioid PhArmacoLogy (OPAL) was a clinical survey of the sociodemographic characteristics, history and consequences of pathology associated with methadone maintenance treatment response and current addictive comorbidities. A subgroup of 72 methadone patients was genotyped.

RESULTS

When comparing the three CYP2B6 genotype groups, the methadone (R)- and (S)-methadone enantiomer concentrations/doses (concentrations relative to doses) were different (P = .029, P = .0019). The CYP2D6 phenotypes did not seem to be relevant with regard to methadone levels. On multivariate analysis, neither the CYP2B6 genotype nor the CYP2D6 phenotype explained the (R)-methadone concentration/dose values (P = .92; P = .86); the (S)-methadone concentration/dose values (P = .052; P = .95 [although there was a difference between the TT group and GT and GG groups {P = .019}]); or opiate cessation (P = .12; P = .90).

CONCLUSION

The genotyping of CYP2B6 G516T could be an interesting tool to explore methadone intervariability.