Abstract
Substance use disorders (SUDs) induce widespread molecular dysregulation in nucleus accumbens (NAc), a brain region pivotal for coordinating motivation and reward, which is linked to neural and behavioral disturbances promoting addiction. Despite the overlapping symptomatology of SUDs, different drug classes exert partly unique influences on neural circuits, cell types, physiology, and gene expression. To better understand common and divergent molecular mechanisms governing SUD pathology, we characterized the cell-type-specific restructuring of the NAc transcriptional landscape after psychostimulant or opioid exposure. We combined fluorescence-activated nuclei sorting and deep RNA sequencing to profile NAc D1 and D2 medium spiny neurons (MSNs) across cocaine and morphine exposure paradigms, including initial exposure, prolonged withdrawal after repeated exposure, and re-exposure post-withdrawal. Our analyses reveal that D1 MSNs display many convergent transcriptional responses between the two drug classes, whereas D2 MSNs manifest highly divergent responses, with morphine causing more adaptations in this cell type. Utilizing multiscale embedded gene co-expression network analysis (MEGENA), we discerned transcriptional regulatory networks subserving biological functions altered by cocaine vs. morphine. We observed largely integrative engagement of overlapping gene networks across drug classes in D1 MSNs, but opposite regulation of key D2 networks, highlighting potential therapeutic gene network targets within MSNs. Analysis of gene regulatory systems at the level of enhancers revealed that morphine engages a unique enhancer landscape in D2 MSNs compared to cocaine. Our findings, and future work leveraging this dataset, will open avenues for the development of targeted therapeutic interventions, addressing the urgent need for more effective treatments for SUDs.