Multi-Omics Signatures of Alcohol Use Disorder in the Dorsal and Ventral Striatum
=================================================================================

* Lea Zillich
* Eric Poisel
* Josef Frank
* Jerome C. Foo
* Marion M. Friske
* Fabian Streit
* Lea Sirignano
* Stefanie Heilmann-Heimbach
* André Heimbach
* Per Hoffmann
* Franziska Degenhardt
* Anita C. Hansson
* Georgy Bakalkin
* Markus M. Nöthen
* Marcella Rietschel
* Rainer Spanagel
* Stephanie H. Witt

## Abstract

Alcohol Use Disorder (AUD) is a major contributor to global mortality and morbidity. Postmortem human brain tissue enables the investigation of molecular mechanisms of AUD in the neurocircuitry of addiction. We aimed to identify differentially expressed (DE) genes in the ventral and dorsal striatum between individuals with AUD and controls, and to integrate the results with findings from genome- and epigenome-wide association studies (GWAS/EWAS) to identify functionally relevant molecular mechanisms of AUD. DNA-methylation and gene expression (RNA-seq) data was generated from postmortem brain samples of 48 individuals with AUD and 51 controls from the ventral striatum (VS) and the dorsal striatal regions caudate nucleus (CN) and putamen (PUT). We identified DE genes using DESeq2, performed gene-set enrichment analysis (GSEA), and tested enrichment of DE genes in results of GWASs using MAGMA. Weighted correlation network analysis (WGCNA) was performed for DNA-methylation and gene expression data and gene overlap was tested. Differential gene expression was observed in the dorsal (FDR<0.05), but not the ventral striatum of AUD cases. In the VS, DE genes at FDR<0.25 were overrepresented in a recent GWAS of problematic alcohol use. The *ARHGEF15* gene was upregulated in all three brain regions. GSEA in CN and VS pointed towards cell-structure associated GO-terms and in PUT towards immune pathways. The WGCNA modules most strongly associated with AUD showed strong enrichment for immune response and inflammation pathways. Our integrated analysis of multi-omics data sets provides further evidence for the importance of immune-and inflammation-related processes in AUD.

## Introduction

Alcohol Use Disorder (AUD) is a major contributor to the global disease burden, with a prevalence of ∼17% among 12-month alcohol users in the US1, 2 and an estimated heritability of 49%3. Knowledge about the molecular mechanisms can foster understanding of causes and promote prevention. Recent genome-wide association studies (GWASs) have identified 29 genetic loci associated with Problematic Alcohol Use (PAU), a proxy of AUD4. While GWASs identify increasing numbers of disease-associated loci, the functional interpretation of many of these findings remains inconclusive. Analyzing the transcriptome can extend the understanding of the molecular mechanisms underlying AUD, by identifying associated gene expression patterns. Findings can in turn be integrated with results from GWASs and epigenome-wide association studies (EWASs) to identify the pathomechanisms underlying disease.

Processes in the central nervous system are considered to play a major role in the etiology of addiction, and the transition from chronic alcohol consumption to AUD5. Therefore, it is of particular interest to examine molecular changes associated with addiction in brain tissue. So far, only few studies have been conducted in postmortem human brain tissue to identify transcriptional changes associated with AUD6-8. These studies mainly focused on the prefrontal cortex (PFC) one important part of the neurocircuitry of addiction9, 10. The first transcriptome-wide study in the PFC found DE genes implicated in neuronal processes, such as myelination, neurogenesis, and neural diseases, as well as cellular processes, such as cell adhesion and apoptosis11. In Brodmann Area 9 downregulation of calcium signaling pathways has been observed in individuals with AUD compared to controls12. In the same study, a weighted gene co-expression analysis (WGCNA) pointed towards modules associated with AUD case/control status, which were enriched for nicotine and opioid signaling, as well as immune processes. Another study in the PFC (Brodmann Area 8) showed that co-expression networks associated with lifetime alcohol consumption were enriched for GWAS signals of alcohol dependence6.

Despite the importance of striatal regions in addiction processes, genome-wide human omics studies of these brain regions are still missing. The striatum is divided into the ventral striatum (VS), consisting of the nucleus accumbens and olfactory tubercle; and the dorsal striatum, comprising the caudate nucleus (CN) and putamen (PUT)13. The nucleus accumbens is involved in mediating motivational processes such as aversion and reward, which play a significant role in the development and maintenance of substance use disorders (SUD)13. In addition to regulating motor function, the CN and PUT are involved in cognitive processes relevant for addiction, such as executive functioning and cognitive control, reinforcement learning and habit formation14. Analyses of omics data from striatal regions could complement the knowledge on global molecular changes in the neurocircuitry of addiction in AUD.

In a recent EWAS of AUD in postmortem brain tissue, we identified differentially methylated CpG-sites and regions in the ventral and dorsal striatum15. Previous studies have shown the utility of integrating epigenetic and transcriptomic data in postmortem brain tissue of SUDs using weighted correlation network analysis (WGCNA)16. WGCNA clusters genes or CpG-sites into co-expressed or co-methylated modules based on correlation matrices. By relating modules to each other, WGCNA can be used for data integration, providing more insights than descriptive overlap. For example, whereas a descriptive comparison of histone H3 lysine 4 trimethylation (H3K4me3) and mRNA expression in individuals with AUD and cocaine use disorder revealed no consistent overlap between H3K4me3 trimethylation and gene expression17, a network analysis identified overlapping modules pointing towards co-expressed genes associated with H3K4me3 trimethylation6. Modules associated with AUD were enriched for CNS functions, such as synaptic transmission and regulation of neurogenesis6. WGCNA has also been used for integrating epigenetic and transcriptomic data and investigating their association with opioid use disorder (OUD) in postmortem human brain, identifying immune-related transcriptional regulation to be enriched in co-expressed and co-methylated modules18.

The aim of the present study was to investigate differential gene expression associated with AUD status in the ventral and dorsal striatum, relate these to GWAS findings, and to integrate the findings with DNA-methylation data using a network approach (WGCNA) in order to identify functionally relevant molecular mechanisms in AUD.

## Materials and Methods

### Samples

Postmortem human brain tissue from CN, PUT and VS of a total of 48 individuals with AUD and 51 control individuals (68% male) was obtained from the New South Wales Tissue Resource Centre at the University of Sydney. The Ethics Committee II of the University of Heidelberg approved the study (reference number 2021-681). After quality control (QC), the total sample sizes for each brain region were NCN = 71, NPUT = 77 and NVS = 63. Phenotypic information was assessed by next-of-kin interviews. Inclusion criteria for this study were: age>18 years, Western European Ancestry, no history of severe psychiatric or neurodevelopmental disorders, or SUDs other than AUD and nicotine use disorder or smoking. AUD was defined as meeting DSM-IV criteria for alcohol dependence and consuming 80g of alcohol a day or more (control group: <20g/day). Descriptive information can be found in Table 1 and Supplementary Table S1.

View this table:
[Table 1.](http://medrxiv.org/content/early/2022/04/05/2021.10.04.21264523/T1)

Table 1. 
Descriptive statistics of demographic data.

### RNA extraction and -sequencing

RNA was extracted from frozen tissue according to the manufacturer’s protocol using the Qiagen RNeasy microKit (Qiagen, Hilden, Germany). The RNA Integrity Number (RIN) of all samples was determined using a TapeStation 4200 (Agilent, Santa Clara, CA). RIN values of 273 samples were larger than 5.5, for which libraries were prepared using the TruSeq Stranded Total RNA Library Prep Kit (Illumina, San Diego, CA). RNA sequencing was performed on the NovaSeq 6000 (Illumina) at the Life & Brain Center in Bonn, Germany with read lengths of 2×100bp and a sequencing depth of 62.5 M read pairs per sample on average. Technical replicates were sequenced for all but four samples.

### DNA extraction and Methylation Profiling

DNA extraction, methylation profiling, and QC was performed as described in Zillich et al.15. In brief, DNA was extracted using the DNeasy extraction kit (Qiagen, Hilden, Germany); the Illumina HumanMethylation EPIC BeadChip and the Illumina HiScan array scanning system (Illumina, San Diego, CA) were used to determine DNA-methylation levels. We used an updated and customized version of the CPACOR pipeline to extract beta values from raw intensities19. Criteria for the removal of samples and probes can be found in Zillich et al.15. In the present analyses, DNA methylation data was included from all subjects from whom gene expression data was available after QC.

### Statistical Analyses

All analyses apart from QC and read mapping were performed using R version 3.6.120. An overview of the analysis workflow can be found in Figure 1. The Benjamini-Hochberg (FDR)21 procedure was used to correct for multiple testing. Differentially expressed genes were considered statistically significant at FDR<0.05. All downstream analyses were performed using genes significantly differentially expressed at FDR<0.25.

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/04/05/2021.10.04.21264523/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2022/04/05/2021.10.04.21264523/F1)

Figure 1. 
Analysis workflow of the present study.

#### Mapping and Quantification

Sequencing quality was determined using FastQC22 and 24 samples (11 cases and 13 controls) were excluded due to insufficient sequencing quality (e.g. strong overrepresentation of sequences, GC distribution). Raw reads were mapped to the human genome (hg38) using HISAT2 (v.2.1.0)23. Quantification was performed with the featureCounts function of the Rsubread package (v.2.0.1)24, with hg38 annotation.

#### Differential Gene Expression Analysis

Differential gene expression was determined using DESeq2 (v.1.26.0)25. Minimal pre-filtering was applied, removing genes with normalized counts <10 for more than two samples. Technical replicates were merged prior to differential expression analysis using the collapseReplicates function as implemented in DESeq2. For the differential gene expression analysis, we included age, sex, RIN, pH-value of the brain, and postmortem interval (PMI) as covariates, because of their known influence on gene expression26-28. To assess residual bias after adjustment for covariates, we generated Q-Q plots and calculated genomic inflation factors (Supplementary Figure 1). We further conducted a variance partition analysis using the variancePartition() function of the corresponding R package29, which confirmed the covariates. Results of this analysis can be found in Supplementary Figure 2. Results were filtered for differentially expressed (DE) genes with an absolute log2 fold change larger than 0.02. Volcano plots displaying up- and downregulation of genes for each brain region are shown in Supplementary Figure 3.

#### Gene-Set Enrichment Analysis

Gene-set enrichment analysis was performed using the R package fgsea (v.1.12.0)30, for which DE genes were ranked according to p-value. Enrichment analysis was performed for Gene-Ontology (GO) terms31 and Hallmark gene sets32 and the results were adjusted using FDR correction.

#### Cell Type Enrichment Analysis

To identify cell type specific expression signatures, we performed cell type enrichment analysis using DE genes (FDR<0.25) from the three brain regions. As a reference gene set for brain cell types, we used the “top ranked cell type-enriched genes based on human data” as provided by McKenzie, Wang 33. These contain the 1000 most enriched genes in a cell type and cover astrocytes, endothelial cells, microglia, neurons, and oligodendrocytes. Using the R package GeneOverlap (v.1.22.0)34, we assessed the overlap of AUD-status associated DE genes with markers from the different cell populations. Results were adjusted for multiple testing using the Benjamini-Hochberg method as implemented in GeneOverlap.

#### Differential Methylation Analysis

Effect sizes and p values for CpG sites were used from the EWAS results as presented in the original publication15. In brief, the EWAS model was based on methylation M-values as the dependent variable and AUD status as the predictor. As covariates, sex, age, postmortem interval (PMI), pH-value, estimated smoking, standardized neuronal cell count, and the first ten principal components of the EPIC array internal control probes were included.

#### WGCNA

Weighted correlation network analyses (WGCNA, v.1.70-3)16 were performed to identify modules of co-expressed genes and co-methylated CpG-sites. We assessed the relationship of these modules with AUD case/control status and tested the overlap between associated modules. WGCNA clusters the input matrix according to a dynamic tree cutting algorithm, using a soft power threshold that approximates the criterion of scale-free topology (Rsigned2>0.80). Resulting soft power thresholds for expression networks were 6 for CN, 5 for PUT, and 14 for VS; for methylation networks, all power thresholds were 2.

To identify methylation networks associated with gene expression, beta values from normalized intensities of all samples from which gene expression data were available were filtered for promoter-associated CpG-sites based on the manufacturer’s manifest (Illumina, San Diego, CA). The resulting 105 796 CpG-sites were used as input.

For the RNA-seq data, count matrices were normalized using the DESeq2 function normalizeCounts and variance stability transformation was applied.

Networks were constructed using following settings: minimum module size=30, mergeCutHeight=0.25, maxBlockSize=36 000. In WGCNA, modules are labeled using colors. In the results section modules are labeled according to type of data, brain region, and color assigned in the analysis, e.g. “e-VS-pink” for module “pink” from the WGCNA analysis of gene expression data in the ventral striatum. For each module, its eigengene was calculated and correlated with AUD status. Association of modules with AUD status and covariates is shown in Supplementary Figure 4. For modules associated with AUD status, we performed enrichment analysis using the GOenrichmentAnalysis function implemented in the WGCNA package for expression data and the R package missMethyl (v.1.20.4)35 for methylation modules. Further, we extracted hub genes of AUD associated WGCNA expression modules by calculating the product of module membership and gene significance for each gene of a module. Based on this score, the 10% of highest ranking genes were defined as module hub genes. To investigate the biological relevance of hub genes, protein-protein interaction networks were generated using the *Search Tool for the Retrieval of Interacting Genes/Proteins* (STRING, v.11.5)36. Graphical representation of gene networks was restricted to high confidence interactions (interaction score threshold 0.7).

#### Expression and Methylation Data Integration

To identify genes both DE and differentially methylated, we analyzed the overlap of DE genes (FDR<0.25) with the results of an EWAS (p<0.001) in the same sample15. We prioritized CpG-sites based on their functional relevance in gene expression regulation. Thus, promoter-associated CpG-sites were used in the analysis.

At the module level, gene-set overlap tests were performed using the R package GeneOverlap (v.1.22.0)34. Here, Fisher’s exact test is used to identify significant overlap. For each brain region, the overlap of the AUD-associated co-expression and co-methylation modules was tested.

#### GWAS Enrichment Analysis

We analyzed enrichment of DE genes with an FDR<0.25, and genes in AUD-associated WGCNA modules in GWAS summary statistics using Multi-marker Analysis of GenoMic Annotation (MAGMA, v.1.08b)37. We performed GWAS enrichment analysis for several SUDs, such as alcohol use disorder and problematic alcohol use4, cannabis use disorder38, and a recent GWAS comparing individuals with opioid use disorder with unexposed controls39. Bonferroni correction (n=4 tests per gene set) of p values was applied to adjust for multiple testing.

## Results

### Differential Gene Expression

Gene expression analysis of postmortem brain tissue from AUD cases and controls revealed DE genes at FDR<0.05 in both dorsal striatal regions. In the caudate nucleus, 49 DE genes were identified at FDR<0.05 (39 up- and 10 downregulated). Tubulin Tyrosine Ligase Like 4 (*TTLL4*, log2FC=0.11, p=2.3*10−8) and GATA Binding Protein 2 (*GATA2*, log2FC=-0.27, p=8.6*10−7) were the most significantly upregulated and downregulated genes, respectively. Top up- and downregulated genes in the putamen were found to be Transcription Elongation Factor A Like 2 (*TCEAL2*, log2FC=0.09, p=5.8*10−5) and Desmin (*DES*, log2FC=-0.86, p=2.6*10−6), the latter being the only significant gene after correction for multiple testing. Nine genes were downregulated in both dorsal striatal regions, with *HLA-DOB* having the highest log2FC in both regions. In the ventral striatum, no DE genes were detected at FDR<0.05. The most significant differential gene expression in the ventral striatum was observed for Ankyrin Repeat And Ubiquitin Domain Containing 1 (*ANKUB1*) which was upregulated in AUD cases (log2FC=1.35, p=5.8*10−5). In the VS of AUD cases Caseinolytic Mitochondrial Matrix Peptidase Chaperone Subunit B (*CLPB*, log2FC=-0.11, p=5.2*10−6) was the most significantly downregulated gene.

None of the DE genes at FDR<0.05 overlapped between multiple brain regions. Therefore, the less conservative significance threshold of FDR<0.25, which was also used for downstream analyses, was applied to compare the overlap of DE genes. At FDR<0.25 the cardiomyopathy associated 5 (*CMYA5*) gene showed an upregulation in both caudate nucleus and putamen. *ARHGEF15* (Rho Guanine Nucleotide Exchange Factor 15) was upregulated in all three brain regions at FDR<0.25. The Top 5 DE genes from each brain region are listed in Table 2; complete summary statistics are listed in Supplementary Table S2 (CN), S3 (PUT), and S4 (VS). Overlap between DE genes in the different brain regions is shown in Figure 2A.

View this table:
[Table 2.](http://medrxiv.org/content/early/2022/04/05/2021.10.04.21264523/T2)

Table 2. 

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/04/05/2021.10.04.21264523/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2022/04/05/2021.10.04.21264523/F2)

Figure 2. 
Venn Diagrams of gene overlap of A) DE genes at FDR < 0.25 in caudate nucleus (CN), putamen (PUT), and ventral striatum (VS), B) genes forming WGCNA expression-modules showing the strongest association with AUD status for CN, PUT, and VS, C) genes forming WGCNA expression-module “e-CN-magenta” and those forming the methylation-modules “m-CN-red” and “m-CN-midnightblue”. Network plots depicting the WGCNA gene expression modules showing the strongest association with AUD: D) module magenta from caudate nucleus, E) module black from putamen and F) module pink from ventral striatum.

### Gene-set Enrichment Analysis

Pathway analysis using a pre-ranked enrichment analysis revealed significant enrichment of dorsal striatum DE genes for several GO terms and Hallmark gene-sets. Genes in the CN were found to be related to cilia- and microtubule-associated GO-terms, while none of the Hallmark gene-sets was significantly enriched. GO-term and Hallmark gene-set analysis in PUT samples showed enrichment for immune processes, such as “acute inflammatory response to antigenic stimuli” (pFDR=0.006) and “adaptive immune response” (pFDR =0.006). In the VS the most significantly enriched GO-terms were also related to cilia and microtubules, as well as antigen processing. All GO-terms and Hallmark gene-set with FDR <0.10 are listed in Supplementary Tables S5 (CN), S6 (PUT), and S7 (VS).

### Cell Type Enrichment Analysis

In the CN, upregulated DE genes were significantly enriched for astrocytic markers (pFDR=7*10−6), whereas an enrichment for endothelial cell marker genes was detected among downregulated genes (pFDR =2*10−7). No significant cell-type enrichment of DE genes was found in the putamen and the ventral striatum. GeneOverlap heatmap visualizations for the three brain regions are displayed in Supplementary Figure 5.

### WGCNA

#### Expression

In the CN, 21 modules with a median size of 352 genes (range: 64-7 259) were identified. Module “e-CN-magenta”, consisting of 328 genes, showed the strongest positive association with AUD status (r=0.42, p=2.89*10−4). In the PUT, of the 25 modules (median size 249 genes, range: 33-5 381) identified, module “e-PUT-black” was most strongly correlated with AUD with a positive direction of effect (r=0.41, p=2.31*10−4). For expression data from the ventral striatum, 16 modules with a median size of 429 genes (range: 35-9 708) were identified; module “e-VS-pink” had the strongest positive association with AUD (r=0.41, p=0.009). Interestingly, in a GO-term analysis the three AUD associated modules were all enriched for immune processes, such as “defense response” and “inflammation response”. Gene network representation of hub genes in modules “e-CN-magenta”, “e-VS-pink”, and “e-PUT-black” revealed the signal transducer and activator of transcription 3 (*STAT3*) gene as a conserved hub node in all three brain regions (Figure 2D-2F). There was also a wide overlap of the genes in the three modules: 174 (22.54%) were partially shared between all three modules corresponding to the three brain regions, while another 21.76% were shared between at least two modules (Figure 2B). A gene network analysis of the 174 shared genes between regions identified *STAT3, TP53, ICAM1, MYC*, and *NFKBIA* as the top 5 hub nodes of the network. A visualization of the network is depicted in Figure 3.

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/04/05/2021.10.04.21264523/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2022/04/05/2021.10.04.21264523/F3)

Figure 3. 
Network plot of genes, co-expressed between the WGCNA modules e-CN-magenta, e-PUT-black and e-VS-pink.

#### Methylation

In the CN, WGCNA resulted in 36 modules with a median size of 346 CpG-sites (range: 66-41 423). Module “m-CN-red”, consisting of 2 117 CpG-sites, showed the strongest association with AUD case control status (r=-0.27, p=0.021). This module was most highly enriched for the biological processes “cell activation” (p=1.52*10−5) and “leukocyte activation” (p=2.09*10−5). In PUT 177 modules were identified (median size=57 CpG-sites, range: 30-42 248). Module “m-PUT-plum” consisted of 70 CpG-sites and was significantly associated with AUD case/control status (r=-0.29, p=0.023) and enriched for the biological processes “positive regulation of I-κB kinase/NF-κB signaling” (p=0.002) and “regulation of I-κB kinase/NF-κB signaling” (p=0.005). WGCNA in the VS methylation data resulted in 85 modules (median size=178 CpG-sites, range: 35-30 370). The module with the strongest association with AUD was “m-VS-lavender” (r=-0.29, p=0.023), which consisted of 117 CpG-sites and was enriched for the molecular function “natural killer cell lectin-like receptor binding” (p=3.43*10−4) and the biological process “susceptibility to natural killer cell mediated cytotoxicity” (p=3.65*10−4). The top 10 enriched GO-terms for all AUD-associated modules can be found in Supplementary Tables S8-S10.

### Expression and Methylation Data Integration

In the CN, 12 genes showed both differential methylation and differential gene expression. DE statistics, EWAS summary statistics and functional annotation for these genes are provided in Supplementary Table S11. No overlap was observed in the VS and PUT. At the module-level, co-expression module “e-CN-magenta” showed significant overlap with the methylation modules “m-CN-red” (p=0.003) and “m-CN-midnightblue” (p=0.014) (Figure 2C), while expression module “e-CN-purple” did not show significant overlap with the methylation modules in CN. Of the 3 AUD-associated expression modules in the VS, only “e-VS-salmon” showed significant overlap with the methylation module “m-VS-turquoise” (p=0.003), but not “m-VS-lavender”. No overlap was observed for gene expression and DNA-methylation in PUT.

### GWAS Enrichment Analysis of DE Genes and WGCNA modules

In the VS, but not in the dorsal striatum, we observed enrichment of DE genes in the GWAS signal of PAU (p=0.045). In the putamen DE genes were enriched for GWAS signal from a study comparing individuals with OUD to unexposed controls (p=0.025). None of the DE genes in any of the brain regions showed enrichment for signals from a GWAS of Cannabis Use Disorder or AUD.

From the WGCNA modules showing the strongest association with AUD, only module e-VS-pink showed significant enrichment for GWAS signals of CUD (p=0.043). None of the findings remained statistically significant after multiple testing correction. Results from the respective analyses are depicted in Supplementary Figure 6 and enrichment p values as well as the number of overlapping genes are displayed in Supplementary Table S12.

## Discussion

In the present study, we identified DE genes, co-expression networks, and pathways associated with AUD in the dorsal and ventral striatum. The results were integrated with DNA-methylation data and results from GWASs of SUDs.

We discovered that one gene (*ARHGEF15*) was consistently upregulated in all investigated brain regions of AUD cases compared to controls. *ARHGEF15* encodes a specific guanine nucleotide exchange factor for the activation of Ras homolog family member A (RhoA), a GTPase, which has been linked to higher blood pressure and hypertension over the Rho/ROCK signaling cascade40. It is postulated that the Rho Guanine Nucleotide Exchange Factor 15 negatively regulates excitatory synapse development by suppressing the synapse-promoting activity of EPHB241. EPHB2 deficiency has been linked to depression-like behaviors and memory impairments in animal studies42. In line with this, genetic variation within *ARHGEF15* has been associated with hematocrit, red blood cell count, and hemoglobin concentration43, but also with psychiatric traits, such as neuroticism and worries44 as well as bipolar disorder45.

Among the genes that were downregulated in both dorsal striatal regions, *HLA-DOB* displayed the highest fold change. HLAs of the Major Histocompatibility Class II are an essential part of the acquired immune system presenting antigens to T-lymphocytes (for review: Howell, Carter 46). The most significantly downregulated gene in the VS is *CLPB*, a mitochondrial chaperone, which has been associated with progressive brain atrophy47 and with the cellular response to alcohol-induced stress48. In a recent GWAS, *CLPB* was associated with the amount of alcohol consumed on a typical day (p=9.67*10−5, N=116 163)49.

DE genes in the ventral striatum were enriched for GWAS signals of PAU, but not AUD. This could be a result of the larger sample size of the PAU GWAS, but also point towards differences in genetic variation as responsible for differential expression.

Our results from the pathway and network analyses further underline immune-related effects of chronic alcohol exposure; the pathway and network modules most strongly associated with AUD case-control status were also enriched for immune system and inflammation processes. This was observed for all three brain regions, and both in expression and methylation data, providing further evidence for the important role of immune processes in AUD.

Gene networks derived from WGCNA hub genes similarly revealed genes related to inflammatory processes as strongly connected network nodes. Here, *STAT3* represents a conserved network hub node in all three brain regions. *STAT3* is a member of the JAK/STAT pathway and acts as a transcription factor upon activation by cytokines, hormones and growth factors50. Interestingly, a recent study assessing expression signatures of alcohol withdrawal in rats discovered a very similar gene network in the hippocampus with *STAT3* as a hub node surrounded by a network of downstream targets51. The authors also discovered increased levels of STAT3 and its neuroinflammation-related target genes in postmortem brain tissue of subjects with AUD. Activation of the *STAT3* gene network was found to be primarily restricted to astrocytes. This supports the results of the cell type enrichment analyses, where enrichment of astrocytic expression signatures was detected for upregulated DE genes in the CN.

These results strongly reflect the well-described effect of chronic alcohol exposure on different aspects of the innate and acquired immune systems52. Chronic alcohol exposure accelerates the inflammatory response and reduces anti-inflammatory cytokines52. An activated immune response in response to chronic alcohol exposure has been shown on the cell level53, as well as on the transcription53, and protein levels54, 55. In a previous EWAS, we found strong enrichment of immune processes in differentially methylated CpG-sites associated with alcohol withdrawal56. Neuroinflammation has been repeatedly associated with AUD and both the glutamate excitotoxicity and the production of acetaldehyde, key processes in AUD metabolism, have been suggested to produce an inflammatory response in the brain57. On a phenotypic level, there is also widespread overlap between symptoms of inflammation and of SUDs, such as anhedonia, depression, and decreased cognitive functioning58. In addition, in candidate gene studies in postmortem human PFC, hippocampus, and orbitofrontal cortex, increased mRNA levels of *HMGB1*, which encodes a proinflammatory cytokine and toll-like receptor genes have been associated with alcohol consumption in AUD cases, providing evidence for chronic neuroinflammation in response to alcohol59-61. Notably, there is an overlap of findings not only on the single-gene level but also on the level of pathways and networks/modules. This overlap underlines that alcohol consumption has common biological effects in different brain regions, i.e., most prominently, effects on immune and inflammation processes.

Several limitations apply to our study. First, we cannot distinguish between effects being a consequence of chronic alcohol consumption or addiction. Second, although we corrected for PMI, which can influence tissue quality as a confounding factor, it cannot be ruled out that other characteristics not easily accounted for, such as cause of death, or blood alcohol level for which the majority of individuals have missing data, influenced gene expression. Third, although the sample size is comparatively large for postmortem brain studies in the addiction field, the small number of differentially expressed genes is likely attributable to limited power. Lastly, analyzing bulk tissue does not adequately reflect the diversity of cell types across different brain regions and future studies on single-cell level are needed to investigate cell-type specific transcriptional changes associated with AUD.

It has to be noted that besides DNA methylation, epigenetic mechanisms such as histone and chromatin modifications, or microRNA expression profiles can influence gene expression and are especially important in addiction research62. Future studies should therefore expand the epigenetic profiling of AUD to include these mechanisms.

In summary, the present study provides further evidence from multi-omics data sets for the importance of immune-and inflammation-related processes in AUD. Notably, drugs that reduce neuroinflammation to reduce drinking, such as phosphodiesterases, may be promising approaches for novel treatment options for AUD. Recently published randomized controlled trials suggest that a phosphodiesterase inhibitor reduces heavy drinking whereas an antibiotic compound was not effective63, 64. A deeper understanding of the underlying mechanisms will enhance the discovery of drug targets and drive forward the development of precision medicine within in this field.

## Supporting information

Supplementary Figure 1 [[supplements/264523_file05.pdf]](pending:yes)

Supplementary Figure 2 [[supplements/264523_file06.pdf]](pending:yes)

Supplementary Figure 3 [[supplements/264523_file07.pdf]](pending:yes)

Supplementary Figure 4 [[supplements/264523_file08.pdf]](pending:yes)

Supplementary Figure 5 [[supplements/264523_file09.pdf]](pending:yes)

Supplementary Figure 6 [[supplements/264523_file10.pdf]](pending:yes)

Supplementary Figure Legends [[supplements/264523_file11.docx]](pending:yes)

Supplementary Tables [[supplements/264523_file12.xlsx]](pending:yes)

## Data Availability

Raw data are available on request. Summary statistics are presented in the Supplementary Materials.

## Conflict of Interest

The authors have nothing to disclose.

## Supplementary Information

Supplementary information is available at MP’s website.

## Data and Code Availability

Raw data and analysis code are available from the corresponding author upon reasonable request.

## Acknowledgments

The study was supported by the German Federal Ministry of Education and Research (BMBF), “A systems-medicine approach towards distinct and shared resilience and pathological mechanisms of substance use disorders” (01ZX01909), “Towards Targeted Oxytocin Treatment in Alcohol Addiction” (031L0190). ERA-NET: Psi-Alc (01EW1908), and the Deutsche Forschungsgemeinschaft (DFG) – Project-ID 402170461 – TRR 26565.

We thank Elisabeth Röbel and Claudia Schäfer-Arnold for technical assistance.

Tissues were received from the New South Wales Brain Tissue Resource Centre at the University of Sydney which is supported by the University of Sydney. Research reported in this publication was supported by the National Institute of Alcohol Abuse and Alcoholism of the National Institutes of Health under Award Number R28AA012725. The content is solely the responsibility of the authors and does not represent the official views of the National Institutes of Health.

## Footnotes

*   We have now increased the visualization of the results. We included volcano plots to exemplify the differential expression in the brain regions, as well as heatmaps of module trait relationships and networks plots to illustrate the WGCNA results. We further performed cell-type enrichment analyses and provided more information on the results in the supplementary material.

*   Received October 4, 2021.
*   Revision received April 4, 2022.
*   Accepted April 5, 2022.


*   © 2022, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## References

1.  1.World Health Organization. Global status report on alcohol and health 2018. World Health Organization 2019.
    
    

2.  2.Grant BF, Chou SP, Saha TD, Pickering RP, Kerridge BT, Ruan WJ et al. Prevalence of 12-Month Alcohol Use, High-Risk Drinking, and DSM-IV Alcohol Use Disorder in the United States, 2001-2002 to 2012-2013: Results From the National Epidemiologic Survey on Alcohol and Related Conditions. JAMA Psychiatry 2017; 74(9): 911–923.
    
    

3.  3.Verhulst B, Neale MC, Kendler KS. The heritability of alcohol use disorders: a meta-analysis of twin and adoption studies. Psychological medicine 2015; 45(5): 1061.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0033291714002165&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25171596&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2021.10.04.21264523.atom) 

4.  4.Zhou H, Sealock JM, Sanchez-Roige S, Clarke T-K, Levey DF, Cheng Z et al. Genome-wide meta-analysis of problematic alcohol use in 435,563 individuals yields insights into biology and relationships with other traits. Nature Neuroscience 2020; 23(7): 809–818.
    
    

5.  5.Volkow ND, Koob GF, McLellan AT. Neurobiologic Advances from the Brain Disease Model of Addiction. New England Journal of Medicine 2016; 374(4): 363–371.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMra1511480&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26816013&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2021.10.04.21264523.atom) 

6.  6.Farris SP, Arasappan D, Hunicke-Smith S, Harris RA, Mayfield RD. Transcriptome organization for chronic alcohol abuse in human brain. Mol Psychiatry 2015; 20(11): 1438–1447.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/mp.2014.159&link_type=DOI) 

7.  7.Ponomarev I, Wang S, Zhang L, Harris RA, Mayfield RD. Gene coexpression networks in human brain identify epigenetic modifications in alcohol dependence. J Neurosci 2012; 32(5): 1884–1897.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Njoiam5ldXJvIjtzOjU6InJlc2lkIjtzOjk6IjMyLzUvMTg4NCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzA0LzA1LzIwMjEuMTAuMDQuMjEyNjQ1MjMuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

8.  8.Liu J, Lewohl JM, Dodd PR, Randall PK, Harris RA, Mayfield RD. Gene expression profiling of individual cases reveals consistent transcriptional changes in alcoholic human brain. J Neurochem 2004; 90(5): 1050–1058.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15312160&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2021.10.04.21264523.atom) 

9.  9.Koob GF, Volkow ND. Neurocircuitry of Addiction. Neuropsychopharmacology 2010; 35(1): 217–238.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/npp.2009.110&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19710631&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2021.10.04.21264523.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000272784600012&link_type=ISI) 

10. 10.Noori HR, Spanagel R, Hansson AC. Neurocircuitry for modeling drug effects. Addict Biol 2012; 17(5): 827–864.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1369-1600.2012.00485.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22978651&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2021.10.04.21264523.atom) 

11. 11.Liu J, Lewohl JM, Harris RA, Iyer VR, Dodd PR, Randall PK et al. Patterns of gene expression in the frontal cortex discriminate alcoholic from nonalcoholic individuals. Neuropsychopharmacology 2006; 31(7): 1574–1582.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/sj.npp.1300947&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16292326&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2021.10.04.21264523.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000238418700023&link_type=ISI) 

12. 12.Kapoor M, Wang JC, Farris SP, Liu Y, McClintick J, Gupta I et al. Analysis of whole genome-transcriptomic organization in brain to identify genes associated with alcoholism. Transl Psychiatry 2019; 9(1): 89.
    
    

13. 13.Volkow ND, Morales M. The Brain on Drugs: From Reward to Addiction. Cell 2015; 162(4): 712–725.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2015.07.046&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26276628&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2021.10.04.21264523.atom) 

14. 14.Galandra C, Basso G, Cappa S, Canessa N. The alcoholic brain: neural bases of impaired reward-based decision-making in alcohol use disorders. Neurol Sci 2018; 39(3): 423–435.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10072-017-3205-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29188399&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2021.10.04.21264523.atom) 

15. 15.Zillich L, Frank J, Streit F, Friske MM, Foo JC, Sirignano L et al. Epigenome-wide association study of alcohol use disorder in five brain regions. Neuropsychopharmacology 2022; 47(4): 832–839.
    
    

16. 16.Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC bioinformatics 2008; 9(1): 559.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/1471-2105-9-559&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19114008&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2021.10.04.21264523.atom) 

17. 17.Zhou Z, Yuan Q, Mash DC, Goldman D. Substance-specific and shared transcription and epigenetic changes in the human hippocampus chronically exposed to cocaine and alcohol. Proc Natl Acad Sci U S A 2011; 108(16): 6626–6631.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTA4LzE2LzY2MjYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wNC8wNS8yMDIxLjEwLjA0LjIxMjY0NTIzLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

18. 18.Liu A, Dai Y, Mendez EF, Hu R, Fries GR, Najera KE et al. Genome-Wide Correlation of DNA Methylation and Gene Expression in Postmortem Brain Tissues of Opioid Use Disorder Patients. International Journal of Neuropsychopharmacology 2021.
    
    

19. 19.Lehne B, Drong AW, Loh M, Zhang W, Scott WR, Tan ST et al. A coherent approach for analysis of the Illumina HumanMethylation450 BeadChip improves data quality and performance in epigenome-wide association studies. Genome Biol 2015; 16(1): 37.
    
    

20. 20.Team RC. R: A language and environment for statistical computing. 2013.
    
    

21. 21.Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society: series B 1995; 57(1): 289–300.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0166-4328(01)00297-2&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1995QE45300017&link_type=ISI) 

22. 22.Andrews S. FastQC: A quality control tool for high throughput sequencing data. [https://www.bioinformatics.babraham.ac.uk/projects/fastqc/](https://www.bioinformatics.babraham.ac.uk/projects/fastqc/), 2010.
    
    

23. 23.Kim D, Paggi JM, Park C, Bennett C, Salzberg SL. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nature Biotechnology 2019; 37(8): 907–915.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41587-019-0201-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31375807&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2021.10.04.21264523.atom) 

24. 24.Liao Y, Smyth GK, Shi W. The R package Rsubread is easier, faster, cheaper and better for alignment and quantification of RNA sequencing reads. Nucleic Acids Research 2019; 47(8): e47–e47.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/NAR/GKZ114&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2021.10.04.21264523.atom) 

25. 25.Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biology 2014; 15(12): 550.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13059-014-0550-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25516281&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2021.10.04.21264523.atom) 

26. 26.Durrenberger PF, Fernando S, Kashefi SN, Ferrer I, Hauw J-J, Seilhean D et al. Effects of Antemortem and Postmortem Variables on Human Brain mRNA Quality: A BrainNet Europe Study. Journal of Neuropathology & Experimental Neurology 2010; 69(1): 70–81.
    
    

27. 27.Trabzuni D, Ryten M, Walker R, Smith C, Imran S, Ramasamy A et al. Quality control parameters on a large dataset of regionally dissected human control brains for whole genome expression studies. Journal of Neurochemistry 2011; 119(2): 275–282.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1471-4159.2011.07432.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21848658&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2021.10.04.21264523.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000295973900003&link_type=ISI) 

28. 28.Birdsill AC, Walker DG, Lue L, Sue LI, Beach TG. Postmortem interval effect on RNA and gene expression in human brain tissue. Cell and Tissue Banking 2011; 12(4): 311–318.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s10561-010-9210-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20703815&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2021.10.04.21264523.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000296636100007&link_type=ISI) 

29. 29.Hoffman GE, Schadt EE. variancePartition: interpreting drivers of variation in complex gene expression studies. BMC Bioinformatics 2016; 17(1): 483.
    
    

30. 30.Sergushichev AA. An algorithm for fast preranked gene set enrichment analysis using cumulative statistic calculation. bioRxiv 2016: 060012.
    
    

31. 31.Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM et al. Gene Ontology: tool for the unification of biology. Nature Genetics 2000; 25(1): 25–29.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/75556&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10802651&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2021.10.04.21264523.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000086884000011&link_type=ISI) 

32. 32.Liberzon A, Birger C, Thorvaldsdóttir H, Ghandi M, Mesirov Jill P, Tamayo P. The Molecular Signatures Database Hallmark Gene Set Collection. Cell Systems 2015; 1(6): 417–425.
    
    

33. 33.McKenzie AT, Wang M, Hauberg ME, Fullard JF, Kozlenkov A, Keenan A et al. Brain Cell Type Specific Gene Expression and Co-expression Network Architectures. Scientific reports 2018; 8(1): 8868.
    
    

34. 34.Shen L. GeneOverlap: Test and visualize gene overlaps. R package version 1.28.0, 2021.
    
    

35. 35.Phipson B, Maksimovic J, Oshlack A. missMethyl: an R package for analyzing data from Illumina’s HumanMethylation450 platform. Bioinformatics 2016; 32(2): 286–288.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btv560&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26424855&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2021.10.04.21264523.atom) 

36. 36.Szklarczyk D, Gable AL, Nastou KC, Lyon D, Kirsch R, Pyysalo S et al. The STRING database in 2021: customizable protein–protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Research 2021; 49(D1): D605–D612.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkaa1074&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33237311&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2021.10.04.21264523.atom) 

37. 37.de Leeuw CA, Mooij JM, Heskes T, Posthuma D. MAGMA: Generalized Gene-Set Analysis of GWAS Data. PLOS Computational Biology 2015; 11(4): e1004219.
    
    

38. 38.Johnson EC, Demontis D, Thorgeirsson TE, Walters RK, Polimanti R, Hatoum AS et al. A large-scale genome-wide association study meta-analysis of cannabis use disorder. The Lancet Psychiatry 2020; 7(12): 1032–1045.
    
    

39. 39.Polimanti R, Walters RK, Johnson EC, McClintick JN, Adkins AE, Adkins DE et al. Leveraging genome-wide data to investigate differences between opioid use vs. opioid dependence in 41,176 individuals from the Psychiatric Genomics Consortium. Mol Psychiatry 2020; 25(8): 1673–1687.
    
    

40. 40.Wirth A. Rho kinase and hypertension. Biochimica et Biophysica Acta (BBA) - Molecular Basis of Disease 2010; 1802(12): 1276–1284.
    
    

41. 41.Margolis SS, Salogiannis J, Lipton DM, Mandel-Brehm C, Wills ZP, Mardinly AR et al. EphB-Mediated Degradation of the RhoA GEF Ephexin5 Relieves a Developmental Brake on Excitatory Synapse Formation. Cell 2010; 143(3): 442–455.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2010.09.038&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21029865&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2021.10.04.21264523.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000283603900017&link_type=ISI) 

42. 42.Zhen L, Shao T, Luria V, Li G, Li Z, Xu Y et al. EphB2 Deficiency Induces Depression-Like Behaviors and Memory Impairment: Involvement of NMDA 2B Receptor Dependent Signaling. Frontiers in Pharmacology 2018; 9(862).
    
    

43. 43.Astle WJ, Elding H, Jiang T, Allen D, Ruklisa D, Mann AL et al. The Allelic Landscape of Human Blood Cell Trait Variation and Links to Common Complex Disease. Cell 2016; 167(5): 1415–1429.e1419.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2016.10.042&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27863252&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2021.10.04.21264523.atom) 

44. 44.Nagel M, Jansen PR, Stringer S, Watanabe K, de Leeuw CA, Bryois J et al. Meta-analysis of genome-wide association studies for neuroticism in 449,484 individuals identifies novel genetic loci and pathways. Nature Genetics 2018; 50(7): 920–927.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2021.10.04.21264523.atom) 

45. 45.Ikeda M, Takahashi A, Kamatani Y, Okahisa Y, Kunugi H, Mori N et al. A genome-wide association study identifies two novel susceptibility loci and trans population polygenicity associated with bipolar disorder. Mol Psychiatry 2018; 23(3): 639–647.
    
    

46. 46.Howell WM, Carter V, Clark B. The HLA system: immunobiology, HLA typing, antibody screening and crossmatching techniques. Journal of clinical pathology 2010; 63(5): 387–390.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6OToiamNsaW5wYXRoIjtzOjU6InJlc2lkIjtzOjg6IjYzLzUvMzg3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDQvMDUvMjAyMS4xMC4wNC4yMTI2NDUyMy5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

47. 47.Saunders C, Smith L, Wibrand F, Ravn K, Bross P, Thiffault I et al. CLPB Variants Associated with Autosomal-Recessive Mitochondrial Disorder with Cataract, Neutropenia, Epilepsy, and Methylglutaconic Aciduria. The American Journal of Human Genetics 2015; 96(2): 258–265.
    
    

48. 48.Tóth ME, Vígh L, Sántha M. Alcohol stress, membranes, and chaperones. Cell Stress and Chaperones 2014; 19(3): 299–309.
    
    

49. 49.Watanabe K, Stringer S, Frei O, Umicevic Mirkov M, de Leeuw C, Polderman TJC et al. A global overview of pleiotropy and genetic architecture in complex traits. Nat Genet 2019; 51(9): 1339–1348.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-019-0481-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2021.10.04.21264523.atom) 

50. 50.Nicolas CS, Amici M, Bortolotto ZA, Doherty A, Csaba Z, Fafouri A et al. The role of JAK-STAT signaling within the CNS. JAKSTAT 2013; 2(1): e22925–e22925.
    
    

51. 51.Chen W-Y, Chen H, Hamada K, Gatta E, Chen Y, Zhang H et al. Transcriptomics identifies STAT3 as a key regulator of hippocampal gene expression and anhedonia during withdrawal from chronic alcohol exposure. Translational Psychiatry 2021; 11(1): 298.
    
    

52. 52.Szabo G, Saha B. Alcohol’s Effect on Host Defense. Alcohol research : current reviews 2015; 37(2): 159–170.
    
    

53. 53.McClintick JN, Tischfield JA, Deng L, Kapoor M, Xuei X, Edenberg HJ. Ethanol activates immune response in lymphoblastoid cells. Alcohol 2019; 79: 81–91.
    
    

54. 54.Donnadieu-Rigole H, Mura T, Portales P, Duroux-Richard I, Bouthier M, Eliaou JF et al. Effects of alcohol withdrawal on monocyte subset defects in chronic alcohol users. Journal of leukocyte biology 2016; 100(5): 1191–1199.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1189/jlb.5A0216-060RR&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27256567&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2021.10.04.21264523.atom) 

55. 55.Yen CH, Ho PS, Yeh YW, Liang CS, Kuo SC, Huang CC et al. Differential cytokine levels between early withdrawal and remission states in patients with alcohol dependence. Psychoneuroendocrinology 2017; 76: 183–191.
    
    

56. 56.Witt SH, Frank J, Frischknecht U, Treutlein J, Streit F, Foo JC et al. Acute alcohol withdrawal and recovery in men lead to profound changes in DNA methylation profiles: a longitudinal clinical study. Addiction (Abingdon, England) 2020; 115(11): 2034–2044.
    
    

57. 57.de Timary P, Stärkel P, Delzenne NM, Leclercq S. A role for the peripheral immune system in the development of alcohol use disorders? Neuropharmacology 2017; 122: 148–160.
    
    

58. 58.Crews FT, Lawrimore CJ, Walter TJ, Coleman LG. The role of neuroimmune signaling in alcoholism. Neuropharmacology 2017; 122: 56–73.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuropharm.2017.01.031&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28159648&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2021.10.04.21264523.atom) 

59. 59.Crews FT, Qin L, Sheedy D, Vetreno RP, Zou J. High mobility group box 1/Toll-like receptor danger signaling increases brain neuroimmune activation in alcohol dependence. Biol Psychiatry 2013; 73(7): 602–612.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2012.09.030&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23206318&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2021.10.04.21264523.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000316685400005&link_type=ISI) 

60. 60.Coleman LG, Jr.., Zou J, Crews FT. Microglial-derived miRNA let-7 and HMGB1 contribute to ethanol-induced neurotoxicity via TLR7. J Neuroinflammation 2017; 14(1): 22.
    
    

61. 61.Vetreno RP, Qin L, Coleman Jr LG, Crews FT. Increased Toll-like Receptor-MyD88-NFκB-Proinflammatory neuroimmune signaling in the orbitofrontal cortex of human alcohol use disorder. Alcoholism: Clinical and Experimental Research 2021; 00: 1–15.
    
    

62. 62.Robison AJ, Nestler EJ. Transcriptional and epigenetic mechanisms of addiction. Nature Reviews Neuroscience 2011; 12(11): 623–637.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrn3111&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21989194&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F04%2F05%2F2021.10.04.21264523.atom) 

63. 63.Grodin EN, Bujarski S, Towns B, Burnette E, Nieto S, Lim A et al. Ibudilast, a neuroimmune modulator, reduces heavy drinking and alcohol cue-elicited neural activation: a randomized trial. Translational Psychiatry 2021; 11(1): 355.
    
    

64. 64.Petrakis IL, Ralevski E, Gueorguieva R, Sloan ME, Devine L, Yoon G et al. Targeting neuroinflammation with minocycline in heavy drinkers. Psychopharmacology (Berl) 2019; 236(10): 3013–3021.
    
    

65. 65.Heinz A, Kiefer F, Smolka MN, Endrass T, Beste C, Beck A et al. Addiction Research Consortium: Losing and regaining control over drug intake (ReCoDe)—From trajectories to mechanisms and interventions. Addict Biol 2020; 25(2): e12866.