Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2023 Jan 19:13:1078894.
doi: 10.3389/fpsyt.2022.1078894. eCollection 2022.

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

Sheila T Nagamatsu  1   2   3 , Gregory Rompala  4 , Yasmin L Hurd  4 , Diana L Núñez-Rios  1   2   3 , Janitza L Montalvo-Ortiz  1   2   3 ; Traumatic Stress Brain Research Group
Collaborators, Affiliations

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

Sheila T Nagamatsu et al. Front Psychiatry. .

Abstract

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.

Keywords: epigenetic; methylation; non-CpG site; opioid; orbitofrontal cortex; postmortem human brain.

PubMed Disclaimer

Conflict of interest statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Figures

FIGURE 1
FIGURE 1
Evaluation of differential mCpHs associated with OUD. (A) Genomic location of differential mCpHs. (B) Cell-type enrichment of differential mCpHs sites. (C) Venn diagram showing the comparison between differential mCpHs and mCpGs.
FIGURE 2
FIGURE 2
Differential mCpHs sites associated with OUD. (A) Circus plot showing hypermethylated and hypomethylated sites. (B) Volcano plot of the differential sites. (C) Genomic location of hypermethylated and (D) hypomethylated differential sites.
FIGURE 3
FIGURE 3
Functional enrichment. Enrichment analysis was conducted using EnrichR (A,B) and STRING software (C,D) for GO biological process, wikiPathways, KEGG, and RCTM.
FIGURE 4
FIGURE 4
Protein-protein interaction network analysis. We used STRING selecting experiments and co-expression with the highest confidence value (0.9).
FIGURE 5
FIGURE 5
Gene-drug interactions with opioids. The figure shows all genes that interact with 15 selected opioids. OPRK1 interacts with most opioid-related drugs.
See this image and copyright information in PMC

References

    1. Lyden J, Binswanger I. The united states opioid epidemic. Semin Perinatol. (2019) 43:123–31. - PMC - PubMed
    1. Lancet T. A time of crisis for the opioid epidemic in the USA. Lancet (London, England). (2021) 398:277. - PMC - PubMed
    1. Deak J, Zhou H, Galimberti M, Levey D, Wendt F, Sanchez-Roige S, et al. Genome-wide association study in individuals of European and African ancestry and multi-trait analysis of opioid use disorder identifies 19 independent genome-wide significant risk loci. Mol Psychiatry. (2022) 27:3970–9. 10.1038/s41380-022-01709-1 - DOI - PMC - PubMed
    1. Zhou H, Rentsch C, Cheng Z, Kember R, Nunez Y, Sherva R, et al. Association of OPRM1 functional coding variant with opioid use disorder: a genome-wide association study. JAMA Psychiatry. (2020) 77:1072–80. - PMC - PubMed
    1. Jang H, Shin W, Lee J, Do J. CpG and non-CpG methylation in epigenetic gene regulation and brain function. Genes. (2017) 8:6. - PMC - PubMed