Differential chromatin accessibility analysis between individuals with and without AUD was performed for each cell type using DESeq2 with the default parameters. Sex, age, and ethnic origin were included as covariates in the models. Genes with p values of less than 0.2 (corrected for multiple-hypothesis testing using the Benjamini-Hochberg method) were deemed significant. Regions residing in promoter regions of genes was determined using R package ChIPSeeker.101 Namely, the function annotatePeak() was used, with parameters: TxDb = TxDb.Hsapiens.UCSC.hg38.knownGene, annoDb="org.Hs.eg.db", and tssRegion = c(−1000, 1000).
