Opioid use disorder (OUD) looms as one of the most severe medical crises currently facing society. More effective therapeutics for OUD requires in-depth understanding of molecular changes supporting drug-taking and relapse. Recent efforts have helped advance these aims, but studies have been limited in number and scope. Here, we develop a brain reward circuit-wide atlas of opioid-induced transcriptional regulation by combining RNA sequencing (RNAseq) and heroin self-administration in male mice modeling multiple OUD-relevant conditions: acute heroin exposure, chronic heroin intake, context-induced drug-seeking following prolonged abstinence, and heroin-primed drug-seeking (i.e., “relapse”). Bioinformatics analysis of this rich dataset identified numerous patterns of molecular changes, transcriptional regulation, brain-region-specific involvement in various aspects of OUD, and both region-specific and pan-circuit biological domains affected by heroin. Integrating RNAseq data with behavioral outcomes using factor analysis to generate an “addiction index” uncovered novel roles for particular brain regions in promoting addiction-relevant behavior, and implicated multi-regional changes in affected genes and biological processes. Comparisons with RNAseq and genome-wide association studies from humans with OUD reveal convergent molecular regulation that are implicated in drug-taking and relapse, and point to novel gene candidates with high therapeutic potential for OUD. These results outline broad molecular reprogramming that may directly promote the development and maintenance of OUD, and provide a foundational resource to the field for future research into OUD mechanisms and treatment strategies.