Symptom-provocation fMRI analysis Functional imaging data were analyzed using statistical parametric mapping (SPM12; http://www.fil.ion.ucl.ac.uk/spm) employing an epoch-related model. The main symptom(s) for each patient was defined and used to construct a contrast of responses to provocation vs. neutral pictures. The ensuing contrasts were masked with an inclusive prefrontal mask and the statistical parametric map thresholded at P < 0.001 uncorrected. Cluster extent was defined using AlphaSIM (threshold P < 0.05). Note that due to weak responses during symptom activation in patients 4 and 6, the SPMs for these patients were first thresholded at P < 0.01 and 0.005 uncorrected, respectively. After the clinical protocol had ended, the projections to the cortical activated areas on the MOCSS for each patient in their main symptomatic dimension from the striatum were calculated using probabilistic tractography. The connectivity between each of the contacts and the fMRI activations was calculated as the number of fibers of the projection of each contact (sphere of diameter 4 mm around the contact) to the cortical activated areas.