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.

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.

Profiling neuronal methylome and hydroxymethylome of opioid use disorder in the human orbitofrontal cortex.

Opioid use disorder (OUD) is influenced by genetic and environmental factors. While recent research suggests epigenetic disturbances in OUD, this is mostly limited to DNA methylation (5mC). DNA hydroxymethylation (5hmC) has been widely understudied. We conducted a multi-omics profiling of OUD in a male cohort, integrating neuronal-specific 5mC and 5hmC as well as gene expression profiles from human postmortem orbitofrontal cortex (OUD = 12; non-OUD = 26). Single locus methylomic analysis and co-methylation analysis showed a higher number of OUD-associated genes and gene networks for 5hmC compared to 5mC; these were enriched for GPCR, Wnt, neurogenesis, and opioid signaling. 5hmC marks also showed a higher correlation with gene expression patterns and enriched for GWAS of psychiatric traits. Drug interaction analysis revealed interactions with opioid-related drugs, some used as OUD treatments. Our multi-omics findings suggest an important role of 5hmC and reveal loci epigenetically dysregulated in OFC neurons of individuals with OUD.

Decoding the role of transcriptomic clocks in the human prefrontal cortex.

Aging is a complex process with interindividual variability, which can be measured by aging biological clocks. Aging clocks are machine-learning algorithms guided by biological information and associated with mortality risk and a wide range of health outcomes. One of these aging clocks are transcriptomic clocks, which uses gene expression data to predict biological age; however, their functional role is unknown. Here, we profiled two transcriptomic clocks (RNAAgeCalc and knowledge-based deep neural network clock) in a large dataset of human postmortem prefrontal cortex (PFC) samples. We identified that deep-learning transcriptomic clock outperforms RNAAgeCalc to predict transcriptomic age in the human PFC. We identified associations of transcriptomic clocks with psychiatric-related traits. Further, we applied system biology algorithms to identify common gene networks among both clocks and performed pathways enrichment analyses to assess its functionality and prioritize genes involved in the aging processes. Identified gene networks showed enrichment for diseases of signal transduction by growth factor receptors and second messenger pathways. We also observed enrichment of genome-wide signals of mental and physical health outcomes and identified genes previously associated with human brain aging. Our findings suggest a link between transcriptomic aging and health disorders, including psychiatric traits. Further, it reveals functional genes within the human PFC that may play an important role in aging and health risk.

Central and Peripheral Immune Dysregulation in Posttraumatic Stress Disorder: Convergent Multi-Omics Evidence.

Posttraumatic stress disorder (PTSD) is a chronic and multifactorial disorder with a prevalence ranging between 6-10% in the general population and ~35% in individuals with high lifetime trauma exposure. Growing evidence indicates that the immune system may contribute to the etiology of PTSD, suggesting the inflammatory dysregulation as a hallmark feature of PTSD. However, the potential interplay between the central and peripheral immune system, as well as the biological mechanisms underlying this dysregulation remain poorly understood. The activation of the HPA axis after trauma exposure and the subsequent activation of the inflammatory system mediated by glucocorticoids is the most common mechanism that orchestrates an exacerbated immunological response in PTSD. Recent high-throughput analyses in peripheral and brain tissue from both humans with and animal models of PTSD have found that changes in gene regulation via epigenetic alterations may participate in the impaired inflammatory signaling in PTSD. The goal of this review is to assess the role of the inflammatory system in PTSD across tissue and species, with a particular focus on the genomics, transcriptomics, epigenomics, and proteomics domains. We conducted an integrative multi-omics approach identifying TNF (Tumor Necrosis Factor) signaling, interleukins, chemokines, Toll-like receptors and glucocorticoids among the common dysregulated pathways in both central and peripheral immune systems in PTSD and propose potential novel drug targets for PTSD treatment.

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.