AUTHOR=Nagamatsu Sheila T. , Rompala Gregory , Hurd Yasmin L. , Núñez-Rios Diana L. , Montalvo-Ortiz Janitza L. , Traumatic Stress Brain Research Group , Alvarez Victor E. , Benedek David , Che Alicia , Cruz Dianne A. , Davis David A. , Girgenti Matthew J. , Hoffman Ellen , Holtzheimer Paul E. , Huber Bertrand R. , Kaye Alfred , Krystal John H. , Labadorf Adam T. , Keane Terence M. , Logue Mark W. , McKee Ann , Marx Brian , Mash Deborah , Miller Mark W. , Noller Crystal , JM-O , Scott William K. , Schnurr Paula , Stein Thor , Ursano Robert , Williamson Douglas E. , Wolf Erika J. , Young Keith A. TITLE=CpH methylome analysis in human cortical neurons identifies novel gene pathways and drug targets for opioid use disorder JOURNAL=Frontiers in Psychiatry VOLUME=13 YEAR=2023 URL=https://www.frontiersin.org/journals/psychiatry/articles/10.3389/fpsyt.2022.1078894 DOI=10.3389/fpsyt.2022.1078894 ISSN=1664-0640 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.