We identified 24 out of 39 modules that were treatment-responsive in the amygdala. The majority of genes clustered into co-expressed modules with only 1,352 out of 19,783 genes in the grey category (i.e., expression profiles that failed to cluster). Six modules were responsive to both working treatments. The up-regulated treatment-responsive modules (greenyellow, orange and darkred) were enriched with known neuronal and GABAergic genes (Table S11), corroborating the cell-type enrichment analysis of the tesa- and feno-regulated genesets from DE. About half of the co-expressed GABAergic genes were also up-regulated by feno and tesa; for example, Slc32a1, the GABA vesicular transporter was up-regulated 29% by feno and 24% by tesa and Abat, the enzyme responsible for catabolism of GABA into succinic semialdehyde was up-regulated 16% for feno and 13% for tesa. We tested the up-regulated modules for overlap with a list of genes known to affect alcohol consumption (based on studies of mutant mice), and found that two of these modules contained more alcohol consumption-related genes than expected by chance (hypergeometric p < 0.01). Over half of the co-expressed alcohol consumption-related genes
