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

Morphine-Driven m6A Epitranscriptomic Neuroadaptations in Primary Cortical Cultures

Opioid overdose is the leading cause of accidental death in the United States and remains a major public health concern, despite significant resources aimed at combating opioid misuse. Neurobiological research to elucidate molecular and cellular consequences of opioid exposure is required to define avenues to explore for reversal of opioid-induced neuroadaptations. Opioids impart well-documented regulation of the transcriptome and epigenetic modifications in the brain, but opioid-induced epitranscriptomic posttranscriptional regulation of RNA is vastly understudied. N6-methyladenosine (m6A) RNA methylation is significantly enriched in the brain and involved in learning, memory, and reward. m6A modifications have not been studied in opioid use disorder, despite being the most common RNA modification. We detected significant regulation of m6A-modifying enzymes in rat primary cortical cultures following morphine treatment, including AlkB Homolog 5 (Alkbh5). The m6a demethylase ALKBH5 functions as an m6A eraser, removing m6A modifications from mRNA. We hypothesized that chronic opioid treatment regulates m6A modifications through modulation of Alkbh5 and profiled m6A modifications in primary cortical cultures following chronic morphine treatment and Alkbh5 knock-down. We observed differential regulation of m6A modifications for a common set of transcripts following morphine or Alkbh5 knock-down, and the two treatments elicited concordant m6A epitranscriptomic profiles, suggesting that a subset of morphine-driven m6A modifications may be mediated through downregulation of Alkbh5 in cortical cultures. Gene Ontology terms of commonly regulated transcripts included serotonin secretion, synapse disassembly, neuron remodeling, and immune response. Thus, we conclude that morphine can drive epitranscriptomic changes, a subset of which may occur in an Alkbh5-dependent manner.

Mapping the epigenomic landscape of post-traumatic stress disorder in human cortical neurons

The study conducted a comprehensive genome-wide analysis of differential 5mC and 5hmC modifications at both CpG and non-CpG sites in postmortem orbitofrontal neurons from 25 PTSD cases and 13 healthy controls. It was observed that PTSD patients exhibit a greater number of differential 5hmC sites compared to 5mC sites. Specifically, individuals with PTSD tend to show hyper-5mC/5hmC at CpG sites, particularly within CpG islands and promoter regions, and hypo-5mC/5hmC at non-CpG sites, especially within intragenic regions. Functional enrichment analysis indicated distinct yet interconnected roles for 5mC and 5hmC in PTSD. The 5mC marks primarily regulate cell-cell adhesion processes, whereas 5hmC marks are involved in embryonic morphogenesis and cell fate commitment. By integrating published PTSD findings from central and peripheral tissues through multi-omics approaches, several biological mechanisms were prioritized, including developmental processes, HPA axis regulation, and immune responses. Based on the consistent enrichment in developmental processes, we hypothesize that if epigenetic changes occur during early developmental stages, they may increase the risk of developing PTSD following trauma exposure. Conversely, if these epigenetic changes occur in adulthood, they may influence neuronal apoptosis and survival mechanisms.

Clinical potential of epigenetic and microRNA biomarkers in PTSD

Molecular studies identifying alterations associated with PTSD have predominantly focused on candidate genes or conducted genome-wide analyses, often encountering issues with replicability. This review aims to identify robust bi-directional epigenetic and microRNA (miRNA) regulators focusing on their functional impacts on post-traumatic stress disorder (PTSD) and their utility in clinical diagnosis, whilst examining knowledge gaps in the existing research. A systematic search was conducted across multiple databases, including Web of Science, Scopus, Global Health (CABI), and PubMed, augmented by grey literature, yielding 3465 potential articles. Ultimately, 92 studies met the inclusion criteria and were analysed to pinpoint significant epigenetic changes with clinically relevant potential in PTSD. The selected studies explored histone modifications, CpG sites, single nucleotide polymorphisms (SNPs), and miRNA biomarkers. Specifically, nine studies examined epigenetic markers, detailing the influence of methylation on chromatin accessibility at histone positions H3K4, H3K9, and H3K36 within a PTSD context. Seventy-three studies investigated DNA methylation, identifying 20 hypermethylated and five hypomethylated CpG islands consistently observed in PTSD participants. Nineteen studies linked 88 SNPs to PTSD, with only one SNP replicated within these studies. Furthermore, sixteen studies focused on miRNAs, with findings indicating 194 downregulated and 24 upregulated miRNAs were associated with PTSD. Although there are epigenetic mechanisms that are significantly affected by PTSD, a granular deconstruction of these mechanisms elucidates the need to incorporate more nuanced approaches to identifying the factors that contribute to PTSD. Technological advances in diagnostic tools are driving the need to integrate detailed participant characteristics, trauma type, genetic susceptibilities, and best practices for robust reporting. This comprehensive approach will be crucial for enhancing the translational potential of PTSD research for clinical application.

Methylome-wide association studies of psychiatric disorders

Neuroepigenetics underscores the significance of the methylome in both normal and pathological brain function, influencing neurobiological processes and psychiatric well-being. This study systematically examines methylome-wide association studies (MWAS) across various psychiatric disorders that utilize array-based and sequencing approaches on blood or brain tissue samples. The findings provide valuable insights into the early stages of neuropathogenesis in psychiatric disorders, revealing altered epigenetic mechanisms. The methylome emerges as a pivotal factor in the development and treatment of psychiatric conditions, potentially offering avenues for identifying therapeutic targets and informing treatment strategies.

Dynamic alterations in gene co-expression networks and gene-transcript associations characterize co-morbidities in cocaine use disorder

Background-

Individuals with cocaine use disorder (CUD) who attempt abstinence experience craving and relapse, which poses challenges in treatment. Longitudinal studies linking behavioral manifestations in CUD to the blood transcriptome in living individuals are limited. Therefore, we investigated the connection between drug use behaviors during abstinence with blood transcriptomics.

Methods-

We conducted a comprehensive longitudinal study involving 12 subjects (9 males, 3 females) with CUD and RNA sequencing on blood collected at a drug-free baseline, and 3, 6 & 9 months thereafter. We categorized subjects into 2 responder groups (high-low) based on scores of drug use variables, and 3 responder groups (low-intermediate-high) on days of abstinence. We investigated differential expression and gene-transcript associations across responder groups at each time point. Lastly, we examined genes that are both co-expressed and showed dynamic expression with time.

Results-

Genes with significant transcript associations between high and. intermediate days of abstinence at 9 months were notably enriched for cannabis use disorder, drinks weekly, and coronary artery disease risk genes. Time-specific gene co-expression analysis prioritized transcripts related to immune processes, cell cycle, RNA-protein synthesis, and second messenger signaling for days of abstinence.

Conclusion-

We demonstrate that abstinence reflects robust changes in drug use behaviors and the blood transcriptome in CUD. We also highlight the importance of longitudinal studies to capture complex biological processes during abstinence in CUD.

Neuronal-specific methylome and hydroxymethylome analysis reveal significant loci associated with alcohol use disorder

Background: Alcohol use disorder (AUD) is a complex condition associated with adverse health consequences that affect millions of individuals worldwide. Epigenetic modifications, including DNA methylation (5 mC), have been associated with AUD and other alcohol-related traits. Epigenome-wide association studies (EWAS) have identified differentially methylated genes associated with AUD in human peripheral and brain tissue. More recently, epigenetic studies of AUD have also evaluated DNA hydroxymethylation (5 hmC) in the human brain. However, most of the epigenetic work in postmortem brain tissue has examined bulk tissue. In this study, we investigated neuronal-specific 5 mC and 5 hmC alterations at CpG sites associated with AUD in the human orbitofrontal cortex (OFC). Methods: Neuronal nuclei from the OFC were evaluated in 34 human postmortem brain samples (10 AUD, 24 non-AUD). Reduced representation oxidative bisulfite sequencing was used to assess 5 mC and 5 hmC at the genome-wide level. Differential 5 mC and 5 hmC were evaluated using the methylKit R package and significance was set at false discovery rate < 0.05 and differential methylation > 2. Functional enrichment analyses were performed, and gene-level convergence was evaluated in an independent dataset that assessed 5 mC and 5 hmC of AUD in bulk cortical tissue. Results: We identified 417 5 mC and 363 5hmC significant differential CpG sites associated with AUD, with 59% in gene promoters. Some of the identified genes have been previously implicated in alcohol consumption, including SYK, DNMT3A for 5 mC, GAD1, DLX1, DLX2, for 5 hmC and GATA4 in both. Convergence with a previous AUD 5 mC and 5 hmC study was observed for 28 genes. We also identified 5 and 35 differential regions for 5 mC and 5 hmC, respectively. Lastly, GWAS enrichment analysis showed an association with AUD for differential 5 mC genes. Discussion: This study reveals neuronal-specific methylome and hydroxymethylome dysregulation associated with AUD, identifying both previously reported and potentially novel gene associations with AUD. Our findings provide new insights into the epigenomic dysregulation of AUD in the human brain.

Single nuclei transcriptomics in human and non-human primate striatum in opioid use disorder

In brain, the striatum is a heterogenous region involved in reward and goal-directed behaviors. Striatal dysfunction is linked to psychiatric disorders, including opioid use disorder (OUD). Striatal subregions are divided based on neuroanatomy, each with unique roles in OUD. In OUD, the dorsal striatum is involved in altered reward processing, formation of habits, and development of negative affect during withdrawal. Using single nuclei RNA-sequencing, we identified both canonical (e.g., dopamine receptor subtype) and less abundant cell populations (e.g., interneurons) in human dorsal striatum. Pathways related to neurodegeneration, interferon response, and DNA damage were significantly enriched in striatal neurons of individuals with OUD. DNA damage markers were also elevated in striatal neurons of opioid-exposed rhesus macaques. Sex-specific molecular differences in glial cell subtypes associated with chronic stress were found in OUD, particularly female individuals. Together, we describe different cell types in human dorsal striatum and identify cell type-specific alterations in OUD.

Neuronal-specific methylome and hydroxymethylome analysis reveal replicated and novel loci associated with alcohol use disorder

Alcohol use disorder (AUD) is a complex condition associated with adverse health consequences that affect millions of individuals worldwide. Epigenetic modifications, including DNA methylation (5mC), have been associated with AUD and other alcohol-related traits. Epigenome-wide association studies (EWAS) have identified differentially methylated genes associated with AUD in human peripheral and brain tissue. More recently, epigenetic studies of AUD have also evaluated DNA hydroxymethylation (5hmC) in the human brain. However, most of the epigenetic work in postmortem brain tissue has examined bulk tissue. In this study, we investigated neuronal-specific 5mC and 5hmC alterations at CpG sites associated with AUD in the human orbitofrontal cortex (OFC). Neuronal nuclei from the OFC were evaluated in 34 human postmortem brain samples (10 AUD, 24 non-AUD). Reduced representation oxidative bisulfite sequencing was used to assess 5mC and 5hmC at the genome-wide level. Differential 5mC and 5hmC were evaluated using the methylKit R package and significance was set at false discovery rate <0.05 and differential methylation >2. Functional enrichment analyses were performed and replication was evaluated replication in an independent dataset that assessed 5mC and 5hmC of AUD in bulk cortical tissue. We identified 417 5mC and 363 5hmC genome-wide significant differential CpG sites associated with AUD, with 59% in gene promoters. We also identified genes previously implicated in alcohol consumption, such as SYK, CHRM2, DNMT3A, and GATA4, for 5mC and GATA4, and GAD1, GATA4, DLX1 for 5hmC. Replication was observed for 28 CpG sites from a previous AUD 5mC and 5hmC study, including FOXP1. Lastly, GWAS enrichment analysis showed an association with AUD for differential 5mC genes. This study reveals neuronal-specific methylome and hydroxymethylome dysregulation associated with AUD. We replicated previous findings and identified novel associations with AUD for both 5mC and 5hmC marks within the OFC. Our findings provide new insights into the epigenomic dysregulation of AUD in the human brain.

Functional genomic mechanisms of opioid action and opioid use disorder: a systematic review of animal models and human studies

In the past two decades, over-prescription of opioids for pain management has driven a steep increase in opioid use disorder (OUD) and death by overdose, exerting a dramatic toll on western countries. OUD is a chronic relapsing disease associated with a lifetime struggle to control drug consumption, suggesting that opioids trigger long-lasting brain adaptations, notably through functional genomic and epigenomic mechanisms. Current understanding of these processes, however, remain scarce, and have not been previously reviewed systematically. To do so, the goal of the present work was to synthesize current knowledge on genome-wide transcriptomic and epigenetic mechanisms of opioid action, in primate and rodent species. Using a prospectively registered methodology, comprehensive literature searches were completed in PubMed, Embase, and Web of Science. Of the 2709 articles identified, 73 met our inclusion criteria and were considered for qualitative analysis. Focusing on the 5 most studied nervous system structures (nucleus accumbens, frontal cortex, whole striatum, dorsal striatum, spinal cord; 44 articles), we also conducted a quantitative analysis of differentially expressed genes, in an effort to identify a putative core transcriptional signature of opioids. Only one gene, Cdkn1a, was consistently identified in eleven studies, and globally, our results unveil surprisingly low consistency across published work, even when considering most recent single-cell approaches. Analysis of sources of variability detected significant contributions from species, brain structure, duration of opioid exposure, strain, time-point of analysis, and batch effects, but not type of opioid. To go beyond those limitations, we leveraged threshold-free methods to illustrate how genome-wide comparisons may generate new findings and hypotheses. Finally, we discuss current methodological development in the field, and their implication for future research and, ultimately, better care.

Circadian rhythm disruptions associated with opioid use disorder in synaptic proteomes of human dorsolateral prefrontal cortex and nucleus accumbens

Opioid craving and relapse vulnerability is associated with severe and persistent sleep and circadian rhythm disruptions. Understanding the neurobiological underpinnings of circadian rhythms and opioid use disorder (OUD) may prove valuable for developing new treatments for opioid addiction. Previous work indicated molecular rhythm disruptions in the human brain associated with OUD, highlighting synaptic alterations in the dorsolateral prefrontal cortex (DLPFC) and nucleus accumbens (NAc)-key brain regions involved in cognition and reward, and heavily implicated in the pathophysiology of OUD. To provide further insights into the synaptic alterations in OUD, we used mass-spectrometry based proteomics to deeply profile protein expression alterations in bulk tissue and synaptosome preparations from DLPFC and NAc of unaffected and OUD subjects. We identified 55 differentially expressed (DE) proteins in DLPFC homogenates, and 44 DE proteins in NAc homogenates, between unaffected and OUD subjects. In synaptosomes, we identified 161 and 56 DE proteins in DLPFC and NAc, respectively, of OUD subjects. By comparing homogenate and synaptosome protein expression, we identified proteins enriched specifically in synapses that were significantly altered in both DLPFC and NAc of OUD subjects. Across brain regions, synaptic protein alterations in OUD subjects were primarily identified in glutamate, GABA, and circadian rhythm signaling. Using time-of-death (TOD) analyses, where the TOD of each subject is used as a time-point across a 24-h cycle, we were able to map circadian-related changes associated with OUD in synaptic proteomes associated with vesicle-mediated transport and membrane trafficking in the NAc and platelet-derived growth factor receptor beta signaling in DLPFC. Collectively, our findings lend further support for molecular rhythm disruptions in synaptic signaling in the human brain as a key factor in opioid addiction.

Association between methylation in the promoter region of the GAD2 gene and opioid use disorder

DNA methylation is one of the epigenetic mechanisms involved in opioid use disorder. GAD2 is a key catalyticase in gamma amino butyric acid (GABA) synthesis from glutamate, that is implicated in opioid-induced rewarding effect. To reveal the relationship and the underlying mechanism between GAD2 gene methylation and opioid use disorder, we first examined and compared the methylation levels in the promoter region of the GAD2 gene in peripheral blood between 120 patients with opioid use disorder and 110 healthy controls by using a targeted approach. A diagnostic model with methylation biomarkers was established to distinguish opioid use disorder and healthy control groups. Correlations between methylation levels in the promoter region of the GAD2 gene and the duration and dosage of opioid use were then determined. Finally, the transcription factors that potentially bind to the target sequences including the detected CpG sites were predicted with the JASPAR database. Our results demonstrated that hypermethylation in the promoter region of the GAD2 gene was associated with opioid use disorder. A diagnostic model based on 10 methylation biomarkers could distinguish the opioid use disorder and healthy control groups. Several correlations between methylation levels in the GAD2 gene promoter and the duration and dosage of opioid use were observed. Transcription factors TFAP2A, Arnt and Runx1 were predicted to bind to the target sequences including several CpG sites detected in the present study in the GAD2 gene promoter. Our findings highlight and extend the role of DNA methylation in the GAD2 gene in opioid use disorder.