Our study focused on human brain transcriptomes and the genetic risks of cigarettes per day (CPD) to investigate the neurogenetic mechanisms of individual variation in nicotine use severity. We constructed whole-brain and intramodular region-specific coexpression networks using BrainSpan's transcriptomes, and the genomewide association studies identified risk variants of CPD, confirmed the associations between CPD and each gene set in the region-specific subnetworks using an independent dataset, and conducted bioinformatic analyses. Eight brain-region-specific coexpression subnetworks were identified in association with CPD: amygdala, hippocampus, medial prefrontal cortex (MPFC), orbitofrontal cortex (OPFC), dorsolateral prefrontal cortex, striatum, mediodorsal nucleus of the thalamus (MDTHAL), and primary motor cortex (M1C). Each gene set in the eight subnetworks was associated with CPD. We also identified three hub proteins encoded by GRIN2A in the amygdala, PMCA2 in the hippocampus, MPFC, OPFC, striatum, and MDTHAL, and SV2B in M1C. Intriguingly, the pancreatic secretion pathway appeared in all the significant protein interaction subnetworks, suggesting pleiotropic effects between cigarette smoking and pancreatic diseases. The three hub proteins and genes are implicated in stress response, drug memory, calcium homeostasis, and inhibitory control. These findings provide novel evidence of the neurogenetic underpinnings of smoking severity.