INTRODUCTION

Difficulties in exerting control over behaviors such as drug seeking are a hallmark of drug addiction. Loss of control leads to persistence in drug taking with larger amounts consumed for longer than intended, despite a desire to quit and despite adverse personal and social consequences (American Psychiatric Association, 2000). Impairments in inhibitory control have been proposed as central to various stages of stimulant dependence, such as increasing susceptibility to initial use, transition to dependence, maintenance, as well as contributing to relapse and difficulties in maintaining abstinence (Goldstein and Volkow, 2011). Response inhibition is particularly informative in investigating control functions as in addition to face validity and objective behavioral indices, its underlying neural substrates likely mediate other self-regulation functions implicated in drug abuse (Tabibnia et al, 2011). Moreover, animal models of impaired response inhibition have proven instructive for addiction (Dalley et al, 2011). Converging evidence from patient lesion, imaging, transcranial magnetic stimulation and electrocorticography provide support for involvement of the pre-supplementary motor area (pre-SMA) in dorsomedial prefrontal cortex (dmPFC) and right ventrolateral PFC (vlPFC), encompassing the anterior insula and inferior frontal gyrus, in response inhibition (Aron et al, 2003; Floden and Stuss, 2006; Levy and Wagner, 2011; Swann et al, 2012; Swick et al, 2011).

One key challenge is identifying the relative contribution of deficient inhibitory function to different stages of stimulant dependence. Chronic stimulant users have difficulties suppressing prepotent responses, and exhibit slowed stopping compared with controls (Ersche et al, 2011; Fernandez-Serrano et al, 2012; Fillmore and Rush, 2002; Kaufman et al, 2003; Monterosso et al, 2005). However, pre-existing genetic and environment vulnerabilities likely interact with short- and long-term effects of drug use on behavior and the brain to produce these deficits. Investigating the underlying traits or endophenotypes that are intermediate between predisposing genes and clinical symptoms may prove useful (Gottesman and Gould, 2003). Endophenotypes are quantifiable traits that are (1) associated with the disorder; (2) heritable; (3) largely state independent; (4) co-segregate with the disorder within families; and (5) overrepresented in non-affected family members relative to the general population. As stimulant dependence is highly heritable (Merikangas and McClair, 2012), this approach can elucidate potential behavioral and associated brain abnormalities independent of drug taking. The role of pre-existing vulnerabilities is poorly understood, but evidence suggests impaired response regulation and control in those with a family history of drug or alcohol dependence (Acheson et al, 2011; Dawes et al, 1997; Tarter et al, 2003). Moreover, weak inhibitory control predicted later illicit drug use in individuals with and without a family history of alcoholism, but more strongly in the former (Nigg et al, 2006). Comparing drug-dependent individuals to their siblings (SIBs), who have not become drug dependent, and both to healthy control participants can further address these issues, as shared impairments would not result from drug exposure. Identifying familial vulnerability markers for drug dependence may further lead to the development of effective, preventive, and therapeutic strategies for individuals at risk.

In accordance with this rationale, inhibitory performance using the stop-signal task (Logan, 1994) has recently been explored as part of a neurocognitive endophenotype for stimulant drug dependence (Ersche et al, 2012b, 2012c). This task requires suppressing already initiated motor responses, enlisting greater inhibitory demands, and less action selection than other action inhibition tasks (Eagle et al, 2008). Stimulant-dependent individuals (SDIs) and their unaffected SIBs demonstrated slowed stopping latencies, and shared white matter abnormalities in the vlPFC in both groups were associated with this slowing (Ersche et al, 2012b, 2012c). Although some gray matter abnormalities were noted in both groups, no association was noted between gray matter integrity and response inhibition difficulties. These results suggest impaired response inhibition as part of a biomarker for familial risk of drug dependence.

Another important aspect of such a neurocognitive biomarker is neural functioning, which provides an important link between brain-based biological underpinnings of the disorder and behavior. Chronic drug abuse is associated with substantial metabolic PFC changes (Goldstein and Volkow, 2011). In go/no-go tasks, non-abstinent cocaine users demonstrated hypoactivity in no-go trials in the insula and pre-SMA (Hester and Garavan, 2004; Kaufman et al, 2003). Reduced activation in SDI was also reported in the anterior cingulate (ACC; Hester and Garavan, 2004; Kaufman et al, 2003; Li et al, 2008), in accordance with medial PFC hypo-responsiveness during errors and high-conflict situations in addiction, believed to reflect diminished error monitoring (Hester and Garavan, 2004; Nestor et al, 2011). Such results suggest executive control-related hypoactivity is not pervasive throughout the PFC with some psychological and neocortical specificity to drug-related inhibitory dysfunction (Kaufman et al, 2003).

We compared functional activation in SDI, their non-affected SIBs, and healthy CTs while performing a stop-signal task. Brain function could reveal additional manifestations of inhibitory vulnerability relating to the SDI endophenotype, but as function is highly adaptable potential protective factors may also be evident. Similar functional abnormalities between SDI and SIB reflecting predisposing vulnerability factors may be of particular utility for prevention and treatment as these could be amenable to change (Vocci, 2008). We anticipated the possibility of shared impairments between the SDI and SIB, as heritability influences on functional activation relating to executive functions have been found in the vlPFC and ACC (Koten et al, 2009; Matthews et al, 2007). Based on previously published findings (Ersche et al, 2012b, 2012c), we hypothesized impaired response inhibition in the sibling-pairs using a modified version of the stop-signal adapted for the MR environment. We further hypothesized hypoactivations in brain regions subserving response inhibition in SDI. Critically, if the SIB demonstrated similar hypoactivations then it would suggest functional difficulties in the motor inhibition system despite no exposure to chronic drug use, thus further characterizing the endophenotype. Alternatively, should the SIB demonstrate hyperactivations or other differences in the relevant neuro-circuitry, this may be indicative of compensatory responses possibly relevant to resilience.

MATERIALS AND METHODS

Participants

Recruitment and screening procedures have been described in detail elsewhere (Ersche et al, 2012b, 2012c). Briefly, all participants were aged 18–55 years, with no history of psychotic or neurodevelopmental disorder (including Attention Deficit/Hyperactivity Disorder), neurological illness, or traumatic head injury, and were fluent in English. Sibling pairs shared the same biological parents with one sibling satisfying Diagnostic and Statistical Manual, Fourth Edition, Text Revision (DSM-IV-TR; American Psychiatric Association, 2000) criteria for cocaine or amphetamine dependence, and the other with no history of substance dependence except nicotine. All participants were psychiatrically evaluated using the Structured Clinical Interview for DSM-IV (First et al, 2002) augmented with a semi-structured interview to ascertain history of drug use, mental, and physical health, including signs of acute intoxication and withdrawal. Participants were recruited by advertisements and from treatment services in East Anglia, UK. Healthy CTs had no personal or family history of drug or alcohol dependence, and were recruited from the community. Verbal IQ was assessed by the National Adult Reading Test (NART; Nelson, 1982). Drug-taking experiences in CT and SIB were assessed with the Drug Abuse Screening Test (DAST-20; Gavin et al, 1989) and Alcohol Use Disorders Identification Test (AUDIT; Saunders et al, 1993). Compulsive drug taking was assessed in the SDI with the Obsessive-Compulsive Drug Use Score (OCDUS; Franken et al, 2002). The SDI had been using stimulants for several years and reported moderate-to-high compulsive symptoms (see Table 1). Of the SDI, 17 met criteria for dependence on opiates, 9 for alcohol and 1 for cannabis. Urine screen results were positive for all, but one, SDI and negative for all SIB and CT. The study was approved by the Cambridge Research Ethics Committee (REC08/H0308/310; PI KD Ersche) and before participation, volunteers provided written informed consent. Data from these individuals as part of a larger sample have been published previously (Ersche et al, 2012b, 2012c).

Table 1 Demographic and Clinical Measures

Stop-Signal Task

Participants viewed the task stimuli via a mirror as they lay in the scanner. On go trials, participants pressed left and right buttons, in response to visual go stimuli (left and right pointing white arrows, 1000 ms). On stop trials, the go stimulus was followed by a stop signal (orange arrow pointing upwards, 300 ms) and participants had to withhold responding. Following initial pilot studies (see Supplementary Information), the MR task was modified to ensure adequate performance in SDI. There were 48 stop trials and 240 go trials, presented intermixed and counterbalanced with left and right, in a single block. The delay between go and stop stimuli, initially set to 250 ms, was adjusted individually by a tracking algorithm in 50 ms steps to allow 50% successful stopping (Logan et al, 1997). If a response was recorded before stop signal onset, it did not appear and the trial was repeated. Inter-trial-intervals were randomly jittered between 900 and 1100 ms (Whelan et al, 2012).

Scanning Acquisition

Functional imaging data were collected in one run on a Siemens TIM Trio 3-Tesla scanner (Erlangen, Germany) using whole-brain echo planer images (EPIs) with the following parameters: repetition time (TR)=2000 ms; echo time (TE)=30 ms; flip angle=78°; 32 slices with slice thickness 3 mm plus 0.75 mm gap; matrix=64 × 64; field of view (FOV)=192 × 192 mm yielding 3 × 3 mm in-plane resolution, and the number of volumes ranged from 278 to 305. T1-weighted scans were acquired for registration purposes (176 slices of 1 mm thickness, TR=2300 ms; TE =2.98 ms, TI=900 ms, flip angle=9°, FOV=240 × 256 mm).

Data Analysis

Mean go reaction time (RT) and stop-signal RT (SSRT) were calculated for each individual. In addition, percent unsuccessful stopping was computed and unsuccessful stop RT was compared with go RT. SSRT was estimated by subtracting mean stop-signal delay from correct go RT in accordance with the race model (Logan, 1994). It became apparent that a subset of participants did not follow instructions, leading to violation of the race model. Full details regarding inclusion criteria adopted to uphold the race model can be found in the Supplementary Information. Following exclusion of three participants due to poor scanning quality, analyses included 41 CT, 32 SDI, and 39 SIB from 143 participants. Behavioral and demographic data were analyzed using χ2 and analysis of variance (ANOVA) tests, significant group differences are followed by Cohen d’s effect size.

Imaging data were processed and analyzed using Statistical Parametric Mapping 8 (SPM, http://www.fil.ion.ucl.ac.uk/spm/). The first five volumes were discarded due to T1-equilibrium effects. Images were realigned and mean EPI image was co-registered to the T1-weighted image, which was segmented and warped to the Montreal Neurological Institute template using new segment and the deformations applied to the EPI volumes, which were resampled to 2 × 2 × 2 mm3. Finally, images were smoothed with a Gaussian kernel of a 8-mm FWHM.

First level analyses were performed using the general linear model as implemented in SPM8. Individual design matrixes modeled two events of interest: successful stops and unsuccessful stops by convolving onset times with a canonical hemodynamic response function with temporal and dispersion derivatives (see Supplementary Information for further details). First level contrasts were computed from stops relative to baseline and failed vs successful stops. Go trials occurred frequently, comprising the general linear model baseline. Successful and unsuccessful stops were entered into second-level random-effect analyses, comprising one-sample t-tests for each group separately and a between-groups ANOVA. Gender was included as a covariate in all second-level analyses. Regions associated with stopping comprised a search area used in the ANOVA (see Supplementary Information). Three anatomical ROIs further examined group differences associated with successful stopping, including the right and left pre-SMA and right anterior Insula/frontal Operculum (aIfO; see Supplementary Information) (Tzourio-Mazoyer et al, 2002). For unsuccessful vs successful stops, the right and left ACC comprised ROIs. All results reported were significant at P<0.05 corrected for family-wise error with small volume correction, and peak voxels in Montreal Neurological Institute coordinates. Where significant between-group differences were found, MarsBaR (Brett et al, 2002) was used to compute 8-mm spheres around peak coordinates for each individual, which were then correlated with SST performance measures, and in the SDI group, drug use measures. Results were further examined while covarying for tobacco and alcohol consumption.

RESULTS

Demographic and Clinical Measures

As seen in Table 1, age (F<1, P>0.30) and verbal IQ (F=1.1.8, P>0.30) did not differ between groups, but gender distribution differed between SDI and SIB with the latter being predominantly male (χ2(2)=19.57, P<0.001). Importantly, SIB and CT both had low drug-taking experiences and alcohol use as reflected by their low DAST-20 and AUDIT scores.

Behavioral Measures

As seen in Table 2, there were no significant group differences in go RT (F(2,109)=1.63, P=0.21). There was a group difference in SSRT (F(2,109)=3.81, P=.025) with SDI slower than CT (P<0.025, d=0.63; see Table 2). SIB were intermediate, and did not significantly differ from either CT (P=0.573, d=0.22) or SDI (P=0.185, d=0.44). The results remained when covarying for gender with significantly slowed SSRT in SDI compared with CT and intermediate performance in SIB. With regard to the race model, the percent of unsuccessful stop trials did not differ between groups (F(2,109)=0.18, P=0.834) and indicated the tracking algorithm was successful. All groups demonstrated faster unsuccessful stop than go latencies in keeping with the race model (CT: F(1, 40)=4.57, P=.039; SDI: F(1, 31)=41.197, P<0.001; SIB: F(1, 38)=13.579, P<0.001). In addition, there were no significant group differences in go RT standard deviation (F(2,109)=2.33, P=0.10) or slowing on go trials following unsuccessful stops vs successful stops (F(2,109)=0.71, P=0.50).

Table 2 Stop-signal Task Performance Measures

Neuroimaging

Each group demonstrated activation in regions commonly activated in this task, including the pre-SMA, the vlPFC encompassing the insula, and IFG in addition to the inferior parietal cortex and occipital cortex (see Figure 1).

Figure 1
figure 1

Significant brain activation maps associated with stopping in each group (P<0.05, Family-wise error). Note. Brain activation maps per group: row (a) stimulant-dependent individuals; row (b) unaffected siblings of stimulant-dependent individuals; row (c) healthy controls. Axial brain slices demonstrate main activation clusters per group at family-wise error P<0.05. The z-coordinate slices are: −4, 8, 24 and 48. Right side of each slice corresponds to the right side of the brain.

PowerPoint slide

In the ANOVA, there was a single significant cluster in the vlPFC (P=0.038 (30, 20, −10), cluster extent (KE)=165, Z=4.28). Direct comparisons indicated CT had greater activation than SDI (P=0.014 (32, 20, −10), KE=203, Z=4.48). Similarly, SIB had greater activation than SDI in a single cluster at close proximity (P=0.040 (32, 22, −6), KE=271, Z=4.22). At this stringent level, there were no differences between CT and SIB. Whole-brain pairwise comparisons at P<0.001 uncorrected in the Supplementary Information demonstrate the specificity of these results.

The ROI analyses replicated the above finding (see Figure 2). In the right aIfO, a main effect of group revealed a single cluster (P=0.025 (34, 20, −10), KE=38, Z=4.00), with CT having greater activation than SDI (P=0.004 (34, 20, −10), KE=60, Z=4.46) and similarly SIB had greater activation than SDI (P=0.045 (34, 24, −6), KE=51, Z=3.78). Additionally, SIB had greater activation than CT in the left pre-SMA (P=0.033 (−12, 8, 46), KE=7, Z=3.47).

Figure 2
figure 2

Group mean percent signal change of 8-mm spheres associated with stopping (a and b) and error monitoring (c). (a) Anterio-insula/frontal operculum, centred around coordinates: x=34, y=20, z=−10; (b) pre-SMA, centred around coordinates: x=−12, y=8, z=46; (c) anterior cingulate, centred around coordinates: x=16, y=32, z=18. Error bars index standard error of the mean. CT, controls; SDI, stimulant dependent individuals; SIB, unaffected siblings.

PowerPoint slide

Comparing unsuccessful to successful stops revealed a main effect of group in right ACC (P=0.001 (16, 32, 18), KE=25, Z=4.53), with both CT and SIB having greater activation than SDI (P<0.001 (16, 32, 18), KE=126, Z=5.03 and P=0.050 (12, 28, 18), KE=11, Z=3.37, respectively). The findings were almost identical when using all stop trials. In none of these comparisons did SDI demonstrate increased activation compared with CT and SIB. Nor did the CT show increased activations compared with SIB. There was a small correlation between aIfO activation and SSRT (t(110)=2.64, P=.009, r=−0.244 at (30, 20, −10)), but no significant correlations with go RT. Finally, there was a correlation between SSRT and duration of use in SDI (t(30)=2.20, P=.036, r=0.373). In sum, the SDI demonstrated vlPFC hypoactivation during stopping and ACC hypoactivation during unsuccessful stopping, compared with the other two groups. The SIB had increased activation in the pre-SMA compared with CT. With the exception of the error-related ACC hypoactivation, the results remained covarying with alcohol and tobacco consumption (see Supplementary Information).

Discussion

We used functional magnetic resonance imaging (fMRI) to elucidate functional correlates of response inhibition associated with family vulnerability for stimulant dependence. In all groups, stopping was associated with similar widespread activations throughout prefrontal, parietal, and occipital regions, although with smaller activation volumes in SDI (Kaufman et al, 2003). Group differences were not widespread, and were restricted to PFC regions subserving response inhibition (Levy and Wagner, 2011; Swick et al, 2011). Stopping in SDI was impaired and associated with relative hypoactivation in several PFC regions. The SIBs exhibited stopping performance that was intermediate between SDI and CTs together with increased dmPFC activation, suggesting increased effort to achieve comparable, if not worse, performance. The diverging results between SDI and SIBs suggest fMRI activation measures associated with response inhibition are unlikely to contribute to an endophenotype. Potential compensatory hyperactivations within the neurocircuitry-mediating response inhibition in the SIBs point to the possibility that underlying vulnerabilities in executive functions may be overcome in individuals at risk of stimulant drug abuse under certain conditions.

Response Inhibition as an Endophenotype for Stimulant Dependence

A putative endophenotype would reflect shared abnormalities in family members rendering them more vulnerable to drug dependence. Despite intermediate performance levels, often reported for potential endophenotypes, abnormal functional correlates differed considerably between SIBs and SDI. The SIBs, manifesting increased familial risk in the absence of drug abuse, did not demonstrate any hypoactivations. Rather they recruited key prefrontal areas successfully and their hemodynamic response resembled CTs rather than SDI. These SIBs have previously shown increased self-report impulsivity and compulsivity compared with CTs and share some abnormalities in brain structure with their stimulant-dependent SIBs (Ersche et al, 2012b, 2010). In fact, shared IFG white matter abnormalities between SDI and SIBs were associated with impaired stopping, although no gray matter abnormalities were detected in regions implicated in the current study. It is possible that measures other than BOLD signal could reveal functional abnormalities in the vlPFC. In addition, tasks sensitive to function in non-PFC areas showing gray matter shared familial abnormalities may reveal potential functional endophenotypes in future. Intriguingly, both the unaffected SIBs and recreational cocaine users have increased gray matter volume in the cerebellum (Ersche et al, 2013), which may prove to be protective. As the present study did not reveal abnormal cerebellar activation, such putative protective factors may involve processes other than response inhibition. Although task-related PFC functional activations have shown good group stability (Blokland et al, 2011) and hereditary influences (Koten et al, 2009; Matthews et al, 2007), these are typically smaller in magnitude than those reported for brain structure (Peper et al, 2007). Taken together, brain structure associated with inhibitory difficulties can be considered a potential endophenotype for stimulant drug dependence, whereas brain function appears to reveal considerably greater levels of plasticity and variability.

Potential for Compensatory Mechanisms in Unaffected SIBs

Compared with CTs, SIBs demonstrated localized dmPFC hyperactivation, accompanying slightly (albeit not significantly) worse performance. As the dmPFC, particularly the pre-SMA, is considered a key neural substrate of response inhibition (Floden and Stuss, 2006), it is plausible that the hyperactivations are indicative of compensatory mechanisms. Intriguingly, right vlPFC hyperactivations have been found with illicit substance use in adolescents (Whelan et al, 2012) and with long-term abstinence in cocaine users (Connolly et al, 2012). Presumably such functional compensations could develop during initial stages of drug use, disappear with chronic use, but reemerge following protracted abstinence. These groups did not show dmPFC hyperactivations as seen in the unaffected SIBs. This then may be a marker of neural vulnerability seen in the absence of drug use and its resulting toxic effects. Alternatively, it may reflect protective factors enabling resilience, as the SIBs avoided drug dependence despite initial vulnerability.

Consistent with this possibility, the SIBs exhibited less SSRT slowing than previously reported (Ersche et al, 2012b). The present task was adapted to enable SDI to perform in the challenging MR environment. That the unaffected SIBs exhibit only mild impairments in attention and executive functions (Ersche et al, 2012c), may have allowed them to benefit from a more salient and less unpredictable stop signal. Accordingly, better performance was observed in the current sample in the MR task than the non-MR version. Compared with the others, this was most pronounced in the SIB, who still had slower SSRTs than CTs in the non-MR version (see Supplementary Information).

The SIB performed better than anticipated compared with SDI and CTs, demonstrating that under present task demands, any underlying abnormalities in structure and predisposing vulnerability factors do not necessarily translate to abnormal functional activation. That individuals may compensate for underlying vulnerability factors, suggests avenues for exploration of protective factors against addiction such as training in volitional control. Functional activations, together with a wider array of cognitive and personality measures (Acheson et al, 2011; Ersche et al, 2012c), may also offer a means of distinguishing between endophenotypes of behavioral disinhibition across different addictions.

Drug-Related Effects on Response Inhibition

The SDI exhibited longer SSRTs and altered regional brain activations that could be attributed to a combination of increased impulsivity, compulsivity, and drug-use effects. Poor inhibitory control was accompanied by reduced activation in the anterior insula extending into the IFG (Kaufman et al, 2003; Li et al, 2008), a region consistently engaged, albeit not exclusively, in inhibitory tasks (Levy and Wagner, 2011; Swick et al, 2011). From present data, we cannot conclude unequivocally that the hypoactivation resulted from inhibitory impairments per se rather than related executive functions such as performance monitoring (Li et al, 2006), specific task characteristics (Aichert et al, 2012), or increased response control demands (Dodds et al, 2011). Nevertheless, the hypoactivations are unlikely to reflect general task difficulties as behavioral group differences were specific to SSRT. With a refined understanding of vlPFC sub-regions still emerging, future studies should examine its role in drug use more closely.

Efficient executive control requires monitoring for errors or conflicting response plans (Botvinick et al, 2004), and is particularly impaired in drug abuse (Goldstein and Volkow, 2011). The ACC hypoactivation associated with failed stopping in SDI is consistent with previous findings in response inhibition tasks (Kaufman et al, 2003; Li et al, 2008) and error/conflict monitoring tasks (Bolla et al, 2004). ACC hypoactivity has been reported in users of other drugs such as opiates (Forman et al, 2004) and cannabis (Eldreth et al, 2004). As the SIB did not show ACC hypoactivation, this may not be a useful endophenotype for stimulant drug use. Similarly, as there were no group differences in post-error slowing, this marker may be less sensitive in SDI (see also Hester et al (2007)) than in other impulsivity-related disorders such as ADHD (Schachar et al, 2004). Deficient error and conflict monitoring associated with ACC hypoactivation could have a role in drug abuse development and maintenance as well as abstinence failure (Connolly et al, 2012). Thus, inhibitory and monitoring dysfunctions likely interact synergistically contributing to disorder severity, relapse, and abstinence maintenance difficulties.

Enhancing cognitive control over drug use remains a key challenge for effective pharmacological and behavioral treatments. A temporal profile of inhibitory deficits as the disorder progresses may highlight a time window, where treatment to improve control could improve clinical outcome, including pharmaceutical treatment (eg, methylphenidate or atomoxetine) and behavioral techniques (Castells et al, 2007). Prevention may be further facilitated by cognitive training, aimed at strengthening the neural networks mediating inhibitory control and self-regulation (Vocci, 2008).

Strengths and Limitations

This study employed a considerably larger sample than previous fMRI studies of response inhibition in SDI, which together with rigorous type-I error control allowed us to highlight specific findings. The stop-signal task has well-defined neural correlates and is of theoretical relevance to drug dependence. The SDI and their SIBs have been well characterized (Ersche et al, 2012c) and this study adds an important dimension to the inhibitory endophenotype in addiction. The SDI did not refrain from taking drugs, but were not under acute administration and completed the task successfully. In short-term abstinence, the findings may be more pronounced and accompanied by additional behavioral and hemodynamic abnormalities (Fillmore et al, 2006). Accordingly, improved performance and increased task-related activation in the PFC are observed with cocaine, methylphenidate as well as with stimulant self-administration (Garavan et al, 2008; Goldstein and Volkow, 2011; Li et al, 2008). Better characterization in acute dosing, non-abstinence, short-, and long-term abstinence could clarify the magnitude of inhibitory dysfunction, although this may be affected by age and duration of use (Ersche et al, 2012a).

Many SDI, as in most previous studies, were polydrug users, with higher levels of tobacco and alcohol consumption. Critically, this is part of the drug dependence phenotype as individuals rarely restrict themselves to a single substance. In fact, such a sample would be unrepresentative and of limited generalizability. Importantly, SDI still exhibited response inhibition deficits when covarying for cigarette and alcohol use (see Supplementary Information). The vlPFC hypoactivation remained, whereas ACC hypoactivation no longer reached significance. Similar analyses in the SIBs still showed dmPFC hyperactivation. This together with evidence that smokers exhibit lower dmPFC activation (de Ruiter et al, 2012), suggests increased smoking in SIB unlikely to account for the differences observed.

In sum, we provide evidence for functional divergence within the neurocircuitry underlying response inhibition and control between SDI and their SIBs with circumscribed PFC dysfunction in SDI, and abnormal overactivations in their SIBs.