To uncover AUD-associated regulatory programs in each cell type, we performed differential module expression analysis. We first pre-processed the pseudo-bulk gene expression count matrix by: (1) normalizing for cell depth, (2) log-transforming, and (3) z-scoring expression across all donors. We then estimated module activity in each donor as the mean expression of module genes. Using a two-tailed two-sample t-test, we identified differentially active modules between individuals with and without AUD.