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 use [4], cannabis use disorder [38], and a recent GWAS comparing individuals with opioid use disorder with unexposed controls [39]. Bonferroni correction (n = 4 tests per gene set) of p values was applied to adjust for multiple testing.