Abstract
Functional polymorphisms in the dopamine transporter gene (DAT1 or SLC6A3) modulate responsiveness to salient stimuli, such that carriers of one 9R-allele of DAT1 (compared with homozygote carriers of the 10R-allele) show heightened reactivity to drug-related reinforcement in addiction. Here, using multimodal neuroimaging and behavioral dependent variables in 73 human cocaine-addicted individuals and 47 healthy controls, we hypothesized and found that cocaine-addicted carriers of a 9R-allele exhibited higher responses to drug cues, but only among individuals who had used cocaine within 72 h of the study as verified by positive cocaine urine screens (a state characterized by intense craving). Importantly, this responsiveness to drug cues was reliably preserved across multimodal imaging and behavioral probes: psychophysiological event-related potentials, self-report, simulated cocaine choice, and fMRI. Because drug cues contribute to relapse, our results identify the DAT1R 9R-allele as a vulnerability allele for relapse especially during early abstinence (e.g., detoxification).
Introduction
The corticobasal ganglia reward circuit, which comprises the ventral striatum and main cortical and subcortical input regions, including the orbitofrontal cortex (OFC), anterior cingulate cortex, and midbrain, is a central hub that mediates the execution of motivated behaviors (Haber and Knutson, 2010). This reward circuit relies heavily on dopamine (DA) neurotransmission, which has core roles in reward prediction (Schultz, 2010), incentive salience (Berridge, 2007), and arousal (Horvitz, 2000), among others. Such DA neurotransmission and consequent responsiveness to salient stimuli are further modulated by functional polymorphisms (i.e., sequence variations) in dopaminergic genes. For example, a polymorphism in the 3′-UTR of the DA transporter gene (DAT1 or SLC6A3) that produces common alleles with 9-repeats (9R) and 10-repeats (10R) has been linked to DA transporter density in the dorsal and ventral striatum (higher density in 9R-allele carriers) (Shumay et al., 2011). Furthermore, the length of the polymorphic region was associated with phasic DA release, wherein carriers of one 9R-allele have lower tonic but higher phasic DA levels than homozygous carriers of the 10R-allele (van Dyck et al., 2005).
Importantly, prior studies have suggested that this DAT1 9R-allele may confer vulnerability for sensitivity to drug cues in addiction. In cigarette smokers, for example, DAT1 9R-allele carriers (compared with smokers homozygous for the 10R-repeat allele) reported stronger craving when exposed to smoking cues (Erblich et al., 2005) or stress (Erblich et al., 2004); 9R-allele smokers also showed more smoking-induced DA release in the caudate and nucleus accumbens as measured with positron emission tomography (Brody et al., 2006), and more responsiveness to smoking-related cues in corticolimbic regions, including the OFC as measured with perfusion fMRI (Franklin et al., 2011), effects that have been correlated with subjective craving (Franklin et al., 2011) and/or attention bias to the drug cues (Wetherill et al., 2012). Together, current evidence raises the intriguing hypothesis that the lower tonic yet increased phasic DA neurotransmission conferred by the DAT1 9R-allele augments response to drug-related cues, an association perhaps further modulated by drug-relevant state variables (e.g., craving).
The current study tested whether drug cue reactivity is modulated by DAT1 in combination with another highly important drug-relevant state variable: recent cocaine use. Recent cocaine use was objectively operationalized as the presence of cocaine metabolites in urine, a variable that modulates response to reinforcement (including to cocaine cues) (Moeller et al., 2010; Parvaz et al., 2012) and predicts treatment outcome (Poling et al., 2007; Ahmadi et al., 2009; García-Fernández et al., 2011). As part of a multimodal imaging and behavioral design, individuals with cocaine use disorder (CUD) and healthy controls were exposed to cocaine- and noncocaine stimuli, during which we measured event-related potentials (ERPs), self-reported valence and arousal ratings, simulated cocaine choice, and fMRI responsiveness. We hypothesized that CUD carriers of a 9R-allele, especially those with positive cocaine urine screens, would show heightened responsiveness to cocaine stimuli.
Materials and Methods
Subjects
Seventy-three CUD (63 males) and 47 healthy controls (41 males) participated in this research, recruited through advertisements, local treatment facilities, and word of mouth; all provided written informed consent in accordance with the local Institutional Review Board. Exclusion criteria were as follows: (1) head trauma (with a loss of consciousness > 30 min); (2) any psychiatric, medical, or neurological disorder requiring hospitalization or regular monitoring (except for substance abuse disorders and highly comorbid disorders, such as depression and post-traumatic stress disorder); (3) current use (within the last 6 months) of psychoactive medications; (4) in the healthy controls, current or past history of substance use disorder (other than nicotine) or any other psychiatric illness; and (5) positive urine screens for drugs of abuse other than cocaine (avoiding possible interference of other drugs when interpreting results, such as the use of benzodiazepines to assuage withdrawal symptoms). These subjects underwent a comprehensive diagnostic interview, which consisted of the following: (1) structured clinical interview for DSM-IV axis I disorders (First et al., 1996); (2) Addiction Severity Index (McLellan et al., 1992), a semistructured interview instrument used to assess history and severity of substance-related problems in seven problem areas (medical, employment, legal, alcohol, other drug use, family-social functioning, and psychological status); (3) Cocaine Selective Severity Assessment Scale (Kampman et al., 1998), measuring cocaine abstinence/withdrawal signs and symptoms (i.e., sleep impairment, anxiety, energy levels, craving, and depressive symptoms) 24 h within the time of interview; (4) Severity of Dependence Scale (Gossop et al., 1992); and (5) Cocaine Craving Questionnaire (Tiffany et al., 1993). This interview identified the following cocaine-related diagnoses in CUD: current cocaine dependence (N = 52), cocaine dependence in partial remission (N = 11), and cocaine dependence in full remission (N = 10). Current axis I comorbidities were identified in nine CUD subjects, including marijuana use disorders (N = 3), alcohol use disorders (N = 6), ecstasy abuse (N = 1), opiate (heroin) dependence (N = 2), major depressive disorder (N = 1), and post-traumatic stress disorder (N = 3). Fifty subjects reported past comorbidities, including marijuana use disorders (N = 33), alcohol use disorders (N = 30), stimulant abuse (N = 1), opiate dependence (N = 3), phencyclidine use disorders (N = 2), major depressive disorder (N = 8), and post-traumatic stress disorder (N = 3).
Genetics screening
Using DNA extracted from the peripheral blood, all subjects were genotyped by PCR as previously described (Shumay et al., 2011) for the 3′UTR VNTR of the SLC6A3. Because cocaine blocks the DA transporter and because of the known deficits in DA functioning in CUD (Volkow et al., 2011), this population is appropriate and important for investigating the effects of DAT1 polymorphisms (although we accordingly focused on DAT1, two other dopaminergic genes were also examined as described in the Control analyses section). From these genetic analyses, and consistent with recent neuroimaging studies in addiction (Franklin et al., 2011; Wetherill et al., 2012), subjects were partitioned into those who were homozygous carriers of the 10R-allele (10R/10R genotype, N = 70) and those who were carriers of at least one 9R-allele (encompassing 9R/10R and 9R/9R genotypes, N = 50). Observed frequency of the major DAT1 genotypes were close to expected according to Hardy-Weinberg assumptions [χ2(1) = 0.38, not significant].
Cocaine urine status
The presence of cocaine metabolites in urine was ascertained with a triage urine panel for abused drugs (Biopsych) in all subjects on study day. Our prior studies have similarly used this cocaine urine status variable (positive vs negative) to subgroup CUD, with effects on neuropsychological functioning (Woicik et al., 2009), reward-associated ERPs (Parvaz et al., 2012), and simulated drug choice (Moeller et al., 2010). Thus, although the urine grouping itself was determined post hoc, the decision to group based on this variable was firmly a priori. Results of this urine screen showed that 35 CUD tested positive for cocaine in urine (CUD+), objectively confirming use within 72 h before study time; 38 CUD tested negative for cocaine in urine (CUD−). Importantly, although CUD+ had used cocaine within 72 h of the study, none was acutely intoxicated upon arrival as determined by the study staff. As expected, cocaine craving was greater in CUD+ than CUD− (across genotypes; independent t(87) = 5.57, p < 0.001). Table 1 provides comprehensive information on the current study sample.
Main study procedures
Study procedures encompassed four main components, all described in detail below.
ERPs
During passive viewing of pleasant, unpleasant, neutral, and cocaine pictures [2000 ms per picture; 30 pictures per picture category; the first three categories selected from the International Affective Picture System (Lang et al., 2008), and the cocaine pictures selected from freely available online and in-house collections (Moeller et al., 2009)], scalp-recorded ERPs [specifically, the late positive potential (LPP), which is thought to index stimulus salience (Hajcak et al., 2010)] were measured by EEG. Consistent with our prior research (Moeller et al., 2012), LPPs for each picture type were defined as the average activity that appeared 400–2000 ms after stimulus onset at the Cz, FCz, FC1, FC2, and Fz electrodes; the average activity in the 200 ms window before picture onset served as the baseline.
Psychophysiological recordings.
A 64 silver-silver chloride electrode cap was positioned according to the International 10/20 System and was used to obtain continuous EEG (Neuroscan) and electro-oculogram recordings (using a frontocentral electrode as ground). Electrodes were placed to record horizontal and vertical eye movements. The EEG was digitized at a rate of 500 Hz and amplified with a gain of 250, and a bandpass filter of 0–70 Hz. The amplifiers were calibrated before each recording. Electrode impedances did not exceed 10 kΩ for any electrodes used in the analysis.
All bioelectric signals were analyzed off-line using Statistical Parametric Mapping (SPM8) for MEG/EEG (Wellcome Department of Cognitive Neurology, London; www.fil.ion.ucl.ac.uk/spm/) and custom MATLAB code (MathWorks). Data were filtered with low and high cutoffs of 0.01 and 30 Hz, respectively, and were then rereferenced to the averaged electrical activity from all 64 scalp sites. The artifact rejection procedure identified a voltage step of >75 μV between sample points and a peak-to-peak voltage difference of 150 μV within an epoch. Additional artifacts were identified through visual inspection and subsequently rejected. Robust averaging was also used to remove artifacts (Wager et al., 2005).
Picture ratings
Immediately after passive viewing, subjects rated each of these pictures on valence (pleasantness) and arousal from 1 to 9 (higher numbers = higher pleasantness/arousal).
Picture choice
Choice to view these pictures was measured under explicit contingencies (when choice was made between two fully visible side-by-side images: explicit task) and under more probabilistic contingencies (when choice was made between pictures hidden under flipped-over cards: probabilistic task) (Moeller et al., 2009). In prior work using these choice tasks, the more CUD chose to view cocaine stimuli specifically over pleasant stimuli, the higher was the actual cocaine use as measured both concurrently (Moeller et al., 2009) and prospectively (Moeller et al., 2013).
Explicit choice task.
The explicit choice task prompted subjects to choose between two simultaneously (side-by-side) presented images of different categories to be enlarged to cover the computer screen with continuous button pressing for up to 5000 ms. If subjects failed to respond for 500 ms, the images returned from their full-screen size to the side-by-side display. The total number of button presses for each category was summed across 70 choice trials and reflected subjects' choice and motivation (through continuous button pressing) to view images of each category.
Probabilistic choice task.
The probabilistic choice task allowed subjects to demonstrate preference under less certain task contingencies than those in the explicit task. In this task, subjects were presented with four virtual and flipped-over card decks on a computer screen. Upon selection of a deck (via a single button press), an image from the selected deck covered the screen for 2000 ms of passive viewing. Each deck contained 26 pictures of a particular category (e.g., cocaine pictures), and four pictures of other categories (e.g., pleasant, unpleasant, neutral) intermixed within each deck for a total of 30 images in each deck. Once a subject selected a particular deck eight times, the next task repetition (run) began with a different arrangement of decks. Each subject completed four runs, and the total number of cards selected from each category was summed across the runs. The probabilistic arrangement of pictures reduced awareness of deck identity but still allowed subjects to establish a preference for certain categories of pictures.
fMRI drug word task
Subjects completed an fMRI drug word task that has been extensively described previously (e.g., Goldstein et al., 2009; Konova et al., 2012). In brief, the task consisted of eight 3.4 min task repetitions (four drug, four neutral), each containing two blocks of 20 drug or neutral words, interleaved with a 20 s white fixation cross overlaid on a black background. Subjects performed each word sequence under one of four counterbalanced monetary reward amounts (50¢, 25¢, 1¢, or 0¢), gained for correct performance for up to $75 of real money (effects of money were not investigated in this study). Each word trial consisted of a 500 ms fixation cross, a 2000 ms word presentation (for word reading), a 500 ms response window (response was made using a Cedrus brand Lumina model LP-400), and a 500 ms feedback slide (correct/incorrect). Word color order was pseudorandomized across all task runs.
Image acquisition.
Scanning was performed on a 4T whole-body Varian/Siemens MRI scanner. BOLD responses were measured as a function of time using a T2*-weighted single-shot gradient-echo EPI sequence (TE/TR = 20/1600 ms, 4 mm slice thickness, 1 mm gap, 33 coronal slices, 20 cm field of view, 64 × 64 matrix size, 90° flip angle, 200 kHz bandwidth with ramp sampling, 128 time points, and 4 dummy scans, discarded to avoid nonequilibrium effects in the fMRI signal). Anatomical images were collected using a T1-weighted 3D-MDEFT sequence (Lee et al., 1995) (TE/TR = 7/15 ms, 0.94 × 0.94 × 0.94 mm spatial resolution, 144 axial slices, 256 readout and 192 × 96 phase-encoding steps, 16 min scan time). A modified T2-weighted hyperecho image (TE/TR = 42/10000 ms, echo train length = 16, 256 × 256 matrix size, 30 coronal slices, 0.86 × 0.86 mm in-plane resolution, 5 mm slice thickness, no gap, 2 min scan time) was also acquired.
Image processing.
Subsequent analyses were performed with the Statistical Parametric Mapping package (SPM8; Wellcome Department of Cognitive Neurology, London) running on MATLAB version 7.7 (MathWorks). A six-parameter rigid body transformation (3 rotations, 3 translations) was used for image realignment and to correct for head motion; 2 mm displacement and 2° rotation in any of the axes in any of the task runs were used as criteria for acceptable motion. Spatial normalization to a standard EPI template (Montreal Neurological Institute) was performed using a 12-parameter affine transformation, resulting in a final voxel size of 3 × 3 × 3 mm. An 8 mm3 full-width at half maximum Gaussian kernel was used to smooth the data. A general linear model and a box-car design convolved with a canonical hemodynamic response function and high-pass filter (cutoff frequency: 1/520 s) were used to calculate individual BOLD-fMRI maps.
Main analyses
Because recent reviews have suggested that cocaine-induced DA levels increase the salience of cocaine rewards while simultaneously decreasing the salience of natural (nondrug) reinforcers (Goldstein and Volkow, 2011), and that comparison between competing reinforcers is important in differentiating addicted individuals from those who self-administer drugs because of a lack of other viable options (Ahmed, 2010), the primary comparison of interest for our dependent measures was between two salient picture categories for CUD: cocaine and pleasant images. To obtain these measures for each dependent measure, raw scores for pleasant variables were subtracted from raw scores for the cocaine variables, creating the cocaine > pleasant contrast for LPPs, valence and arousal ratings, explicit choice, and probabilistic choice. However, given the centrality of the drug > neutral comparison in the addiction literature, we also created and inspected contrasts in which raw scores for neutral variables were subtracted from raw scores for the cocaine variables. For the fMRI data, the primary contrast of interest was percentage signal change for the drug word block minus the fixation baseline for each subject (drug word > baseline fixation); indeed, because there were no pleasant, non–drug-related words in the task, the drug > fixation contrast represents the closest analog to the picture tasks. However, similarly to the other variables, a second contrast of interest was also inspected, reflecting percentage signal change for the drug word block minus the neutral word block for each subject (drug word > neutral word). For all fMRI analyses, we inspected a priori (and independent) ROIs, all created using PickAtlas: anatomical masks in the lateral OFC [encompassing Brodmann Areas (BA) 11 and 47] and medial OFC (encompassing portions of BA 11 and the gyrus rectus that were distinct from the lateral OFC ROI), and a 10 mm spherical volume of interest in the ventral striatum centered at the peak coordinates from Franklin et al. (2011) as follows: x = ±4, y = 4, z = −6. The OFC was our main ROI because this region is critical in computing subjective value, even in situations not requiring a choice (Padoa-Schioppa and Cai, 2011); importantly, this region was also shown to be modulated by DAT1 (Franklin et al., 2011). Although this latter study showed a 9R-allele > 10R/10R effect that was primarily localized to the medial OFC (Franklin et al., 2011), inspection of the cluster indicates that it also extended to more lateral portions of the OFC, justifying our approach of inspecting both OFC subregions. Also leaning on this same study (Franklin et al., 2011), we tested for a similar pattern of effects in the ventral striatum.
Statistical analyses primarily consisted of between-group 3 (diagnosis: CUD+, CUD−, control) × 2 (DAT1 genotype: 10R/10R vs 9R/10R or 9R/9R) ANOVAs, but we also performed correlation analyses between our dependent variables and select drug use variables (specifically, those that are relevant to abstinence/recent drug use: craving, withdrawal, and duration of cocaine abstinence). To analyze the fMRI data, the selected contrasts were entered into one-way ANOVAs in SPM8, where subjects were partitioned by diagnosis (CUD+, CUD−, control) and DAT1 genotype (10R/10R vs 9R/10R or 9R/9R). Then, for each subject, fMRI BOLD activity (drug > fixation or drug > neutral) in each entire ROI (lateral and medial OFC, ventral striatum) was extracted using MARSBAR and compared between the groups in SPSS. In addition to these central ROI analyses, we tested for whole-brain diagnosis × DAT1 interactions at p < 0.05 family-wise error corrected at the voxel-level. For all statistical analyses, all subjects with available data were included. Significance was set at p < 0.05 for all ANOVAs, and p < 0.01 for all correlations (the latter to minimize type I error).
Control analyses
Apart from our main analyses, we also performed control analyses that were meant to bolster our main results. These control analyses proceeded in four steps.
Control tasks
To rule out potentially confounding factors, we analyzed results from two additional neuropsychological tasks of attention and inhibitory control, respectively, for which the dependent variables were unspecific to the subjects' illness (cocaine addiction): the Attention Network Task (ANT) (Fan et al., 2002) and the Stroop Color-Word Test (Stroop, 1935). In particular, the ANT was administered to account for effects of global attention. During this task, subjects responded quickly to neutral visual cues in different directions on a computer screen. Accuracy and reaction time were used to measure three factors of attention, including alerting (response readiness), orientating (scanning/selection of information), and executive control (conflict resolution). The Stroop Color-Word Test (Stroop, 1935), a classical executive function task (measuring suppression of automatic response tendencies), was administered to account for effects of executive function. During this task, subjects were asked to name color words printed in either their congruent or incongruent colors; their contrast provides a measure of interference that was examined in the current study. By including these control tasks (i.e., that used neutral, non–drug-associated contexts), we were able to test whether our results were specific to a salient, drug-associated context or could instead be explained by more basic, but related, cognitive functions. These tasks were completed as part of a comprehensive neuropsychological battery described in detail previously (Woicik et al., 2009).
Covariate analysis
Analyses of covariance (ANCOVAs) were used to control for demographic and other relevant clinical variables that could potentially explain our findings (i.e., continuous variables that showed similar diagnosis × DAT1 interactions, or categorical variables that showed group differences between the six cell groupings resulting from splitting variables as a function of diagnosis × DAT1; Table 1).
Additional ERP component
Although we had a priori hypotheses for the LPP, we tested an additional ERP component. Specifically, we probed for potential effects of the N2, which is a negativity that occurs temporally before the LPP and that results from attention to an unexpected stimulus (such as during an oddball paradigm) (Patel and Azzam, 2005). Because our tasks were expected to tap into salience processing rather than general attention processing, we did not expect a diagnosis × DAT1 interaction for the N2.
Additional dopaminergic genes
Although we had a priori hypotheses for the DAT1 gene, we examined two additional DA-related genes that have been associated with addiction and/or inhibitory control. In particular, along with 3′UTR VNTR of the SLC6A3, all subjects were genotyped for the DRD4- and PER2-VNTRs. The DRD4 exon III VNTR has been previously associated with poorer inhibitory control (Congdon et al., 2008); it has also been associated with substance abuse through novelty seeking (Ray et al., 2009), perhaps especially in men (Laucht et al., 2005, 2007). The VNTR polymorphism in intron 3 of PER2 has been associated with cocaine addiction and striatal DA D2 receptor (DRD2) availability (Shumay et al., 2012). For DRD4, after partitioning our sample based on the presence of the 7R-allele, we performed three separate 3 (diagnosis: CUD+, CUD−, control) × 2 (DRD4 genotype: 7R-allele present vs 7R-allele absent) ANOVAs; for PER2, after partitioning our sample into three genotype groups (high/high, high/low, and low/low), we performed three separate 3 (diagnosis: CUD+, CUD−, control) × 3 (PER2 genotype: high/high, high/low, and low/low) ANOVAs. The dependent variables in these ANOVAs were the dependent variables that showed significant diagnosis × DAT1 interactions in the main results.
Results
Main dependent measures: cocaine versus pleasant
LPPs, ratings, and choice
We analyzed our main dependent measures using 3 (diagnosis: CUD+, CUD−, control) × 2 (DAT1 genotype: 10R/10R vs 9R/10R or 9R/9R) ANOVAs, with dependent variables that included cocaine > pleasant LPPs, self-reported valence and arousal, and probabilistic and explicit choice (see Table 2 for statistics). Both CUD groups showed greater cocaine > pleasant responsiveness than the healthy controls across all five variables as expected, with diagnosis main effects reaching significance for cocaine > pleasant valence, arousal, and probabilistic choice (Table 2). No group main effects of DAT1 reached significance. Importantly, there was a diagnosis × DAT1 interaction for the LPPs (p = 0.016), explained by a urine group difference (CUD+ > CUD−) in the 9R-allele carriers but not in the 10R/10R subjects (Fig. 1A). Although the omnibus diagnosis × DAT1 interactions did not reach significance for the other four variables, in a pattern similar to that of the LPP, and consistent with our hypotheses, valence ratings (Fig. 1B) and probabilistic choice (Fig. 1C) showed the expected urine group differences (CUD+ > CUD−) that emerged only in the 9R-allele subjects.
fMRI response
We also analyzed the three fMRI ROIs (lateral OFC, medial OFC, ventral striatum) with 3 (diagnosis: CUD+, CUD−, control) × 2 (DAT1 genotype: 10R/10R vs 9R/10R or 9R/9R) ANOVAs, but this time with the drug > fixation fMRI response as the dependent variables. The lateral OFC ROI showed a main effect of group, no main effect of DAT1, and a diagnosis × DAT1 interaction (p = 0.011). Consistent with results for the other modalities, this interaction was explained by a urine group difference (CUD+ > CUD−) in the 9R-allele carriers but not in the 10R/10R subjects (Fig. 1D). The other two ROIs showed no main effects or interactions (Table 2).
Combined modalities
Here we examined responsiveness to the cocaine cues across modalities. z-scores were computed for each of the eight dependent variables (LPPs, valence, arousal, explicit choice, probabilistic choice, fMRI lateral OFC, medial OFC, and ventral striatum), which were then averaged to create a composite index of reactivity. The diagnosis × DAT1 interaction was significant (p = 0.017), as expected driven by a urine group difference (CUD+ > CUD−) in the 9R-allele carriers but not in the 10R/10R subjects (Table 2). Indeed, a subsequent planned comparison of this composite score showed that the 9R-allele CUD+ exhibited higher cocaine-related responsiveness compared with the other five diagnosis-genotype group combinations (t(114) = 3.97, p < 0.001). Even after removing three potential outliers (as identified by a diagnosis × genotype boxplot on this aggregate variable), both the interaction and the central CUD+ > CUD− comparison remained significant (both p < 0.01). Thus, although the 9R-allele CUD− subgroup was relatively small (N = 10), our effects were not driven by an extreme outlier in this subgroup.
Main dependent measures: cocaine versus neutral
Although our primary hypotheses pertained to the cocaine > pleasant contrast, we also analyzed the more commonly reported cocaine > neutral contrast (see Table 3 for statistics). We again used 3 (diagnosis: CUD+, CUD−, control) × 2 (DAT1 genotype: 10R/10R vs 9R/10R or 9R/9R) ANOVAs, for each of the eight dependent measures. Results of these ANOVAs revealed diagnosis main effects for valence, arousal, probabilistic choice, and explicit choice, again such that both CUD groups showed higher cocaine > neutral responsiveness than controls (Table 3). There were again no significant main effects of DAT1. The only diagnosis × DAT1 interaction to reach significance, including a newly created cocaine > neutral composite variable, was for cocaine > neutral explicit choice (p = 0.027); and even for this variable, the CUD+ > CUD− comparisons did not reach significance (Table 3). Together, compared with results for the cocaine > pleasant contrast, results for the cocaine > neutral contrast generally exhibited similar directionality but were of weaker magnitude, further supporting our a priori decision to focus on the cocaine > pleasant contrast.
Correlations
To provide additional validity to the results using cocaine urine status, correlation analyses were conducted with measures relevant to abstinence in CUD. Specifically, for the dependent variables that showed significant diagnosis × DAT1 interactions (cocaine > pleasant LPPs, fMRI drug > fixation OFC response, and the combined score of cocaine responsiveness across modalities), we examined their correlations with (1) cocaine craving (total score) (Tiffany et al., 1993), (2) cocaine withdrawal (total score) (Kampman et al., 1998), (3) current frequency of cocaine use (days per week, last 30 d), and (4) duration of current cocaine abstinence (with current frequency of cocaine use and duration of current cocaine abstinence assessed during the clinical interview). These correlations were conducted separately as a function of DAT1 genotype to examine whether correlations differed between 9R-allele carriers and the 10R/10R group, and tests of correlation differences were computed as appropriate. We hypothesized that carriers of a 9R-allele would show higher correlations between cocaine-related responsiveness with craving, withdrawal, current use, and/or abstinence.
Supporting this hypothesis, only in the 9R-allele subjects, the combined cocaine reactivity score (averaged across all eight dependent variables) correlated at the specified p < 0.01 threshold with greater cocaine craving (9R-allele, r = 0.58, p = 0.002; 10R/10R, r = 0.34, p = 0.021; although the correlation difference was not significant, z = 1.20, p = 0.23) (Fig. 2A) and fewer days of cocaine abstinence (9R-allele, Spearman r = −0.62, p = 0.001; 10R/10R, r = −0.07, p = 0.65; correlation difference, z = 2.59, p = 0.015) (Fig. 2B). We observed a similar pattern of effects between cocaine > pleasant LPPs and higher self-reported cocaine craving (9R-allele, r = 0.62, p = 0.002; 10R/10R, r = 0.06, p = 0.68; correlation difference, z = 2.40, p = 0.016). No other correlations reached significance at p < 0.01.
Control analyses
Control analyses helped bolster conclusions from our main analyses. (1) Neutral tasks: Diagnosis × DAT1 interactions were not found when analyzing nonsalient control tasks (Stroop Color-Word Test, p > 0.4; ANT alerting, p > 0.7; and ANT conflict, p > 0.5); and importantly, on these neutral tasks, there were no urine group differences (as seen in LPPs, pleasantness ratings, probabilistic choice, fMRI, and composite score) in either 9R-allele carriers or the 10R/10R genotype for the ANT subscales (all p > 0.05) or Stroop (all p > 0.2). (2) Covariates: Although severity of dependence (Gossop et al., 1992) showed a diagnosis (urine status) × DAT1 interaction (F(2,114) = 4.47, p < 0.05) in CUD, this interaction showed a different pattern of effects than for the cocaine > pleasant reactivity variables (higher severity scores for 9R-allele CUD−); in addition, when covaried in an ANCOVA, this variable did not attenuate the diagnosis × DAT1 interaction for LPP, fMRI, or composite score (all p < 0.05). Similarly, another potential covariate, smoking history, which is a categorical variable that differed between the diagnosis and DAT1 groups (Table 1), did not attenuate the diagnosis × DAT1 interaction for LPP, fMRI, or composite score when covaried in an ANCOVA (all p < 0.05). (3) Additional ERP component: The diagnosis × DAT1 interaction for cocaine > pleasant N2 did not reach significance as expected (p > 0.2). (4) Additional dopaminergic genes: Diagnosis × genotype interactions (on cocaine > pleasant LPPs, drug > neutral lateral OFC response, combined cocaine reactivity measure) were not significant for either DRD4 or PER2 (all p > 0.1).
Discussion
The dopamine transporter (DAT) is a main regulator of extracellular DA levels in striatum (Gainetdinov et al., 1999), and preclinical studies showed that it modulated conditioning to cocaine (Medvedev et al., 2005). In the current study, we hypothesized and found that DAT1 played a prominent role in modulating drug cue reactivity in human cocaine addiction, with the most pronounced reactivity emerging for carriers of a DAT1 9R-allele who had used cocaine within 72 h of the study (evidenced by cocaine metabolites in urine). Notably, this pattern of increased reactivity in urine-positive 9R-allele carriers was preserved across multimodal dependent variables (cocaine > pleasant LPPs, drug > fixation fMRI lateral OFC response, and a composite score reflecting these two neuroimaging variables plus self-reports, simulated cocaine choice behavior, and two additional ROIs of high relevance). Correlation analyses with cocaine craving, and with days since last use, showed the association with these abstinence-related measures when continuously measured. Control analyses showed that these results were not driven by covariates and also did not emerge for other DA genes, other global attentional processing or inhibitory control tasks, or for an additional ERP component. We interpret our findings according to the modulatory role of DAT1 on phasic/tonic DA release. In particular, the 9R-allele contributes to a neural phenotype characterized by lowered tonic but increased phasic DA firing (van Dyck et al., 2005), which in turn is associated with increased reward-related responsiveness (Dreher et al., 2009; Forbes et al., 2009). In CUD, this sensitive DA profile may predispose individuals toward increased reactivity to drug cues in dopaminergically innervated brain regions inclusive of the OFC (Franklin et al., 2011).
One explanation of why cocaine responsiveness in the 9R-allele carriers was most powerfully expressed in CUD+ is that these individuals may have been in short-term cocaine withdrawal, a state of perturbed DA functioning previously found to be modulated by DAT1. For example, 9R-allele carriers exhibited greater susceptibility to withdrawal-related delirium and seizures in alcoholism (Sander et al., 1997), cocaine-induced paranoia in cocaine addiction (Gelernter et al., 1994), and psychosis (lasting for at least 1 month after discontinuation of use) in methamphetamine addiction (Ujike et al., 2003). Previous studies have indicated that CUD in both short-term and longer-term withdrawal exhibit decreased striatal DA activity (Volkow et al., 1997), including decreases in nonstimulated DA release (Martinez et al., 2009). Considering that this dysregulated DA state might be accentuated in CUD 9R-allele carriers (because of reduced tonic DA activity associated with this allele), and that low tonic DA levels are expected to enhance phasic DA cell firing in response to conditioned cues (Wanat et al., 2009), it indeed follows that these subjects showed heightened sensitivity to drug cues especially while possibly experiencing acute withdrawal (i.e., when tonic DA levels were ostensibly reduced above and beyond the influence of the DAT1 9R-allele). Similar logic was put forth to explain the ability of a DA agonist, small-dose oral methylphenidate (which elevates tonic DA levels), to attenuate phasic DA responses to conditioned drug cues in CUD (Volkow et al., 2010). In addition to withdrawal, another potential mechanism underlying our urine status effects may have involved drug priming (i.e., because CUD+, although not acutely intoxicated during the study, self-administered cocaine within 72 h of the study). Similarly to withdrawal, drug priming leads to increased craving and drug-taking (Donny et al., 2004; Mahoney et al., 2007) that is possibly undergirded by the OFC (in concert with other interconnected regions) (Schmidt et al., 2005). Although the mechanisms underlying our urine status effects therefore remain to be clarified in future studies (i.e., whether withdrawal or priming mainly drive our effects), urine status itself is important for addiction research insofar as it predicts poorer treatment outcome (Poling et al., 2007; Ahmadi et al., 2009; García-Fernández et al., 2011).
Our results have important implications for theorizing about the DAT1 polymorphism. First, our results help inform previous studies suggesting that 9R carriers may be more sensitive to an environmental challenge (in the current study, short-term abstinence). For example, 9R-allele carriers with attention deficit/hyperactivity disorder showed a stronger correlation between low childhood maternal warmth and increased conduct and emotional problems later in life than non–9R-allele carriers (Sonuga-Barke et al., 2009). The 9R-allele was also associated with an increased lifetime risk of post-traumatic stress disorder (Chang et al., 2012), indicative of intransigent sensitivity to a distressing experience. Second, although the presence of the 9R-allele has been suggested to reduce the need for novelty and reward by external stimuli (Sabol et al., 1999), our data suggest the opposite conclusion when such individuals have recently used cocaine. To further test this gene × abstinence hypothesis, future studies can expose addicted individuals to drug cues while experimentally controlling the length of current drug abstinence (e.g., before vs after a laboratory-supervised dose of cocaine in individuals then prospectively monitored for abstinence). More broadly, potentially interesting new studies can be undertaken to test the hypothesis of satiation versus deprivation even beyond drug addiction, extending into other psychopathologies with possible dopaminergic underpinnings that are similarly characterized by excessive reward seeking and impaired inhibition (e.g., obesity, pathological gambling).
Limitations of this study include the following: (1) because the same pictures were used for all picture tasks, repeated exposure to the pictures may have reduced participants' responsiveness to the images. Indeed, because the ERP recordings were administered first, we cannot rule out the possibility that null interaction effects for the other dependent measures was the result of habituation. Future studies could use different image sets for different tasks; here, the goal was consistency of stimuli across picture-viewing modalities. (2) Because the fMRI task did not include pleasant stimuli (and used words instead of pictures), future studies could integrate pleasant stimuli into an fMRI picture cue-reactivity task. This design would also allow for follow-up on the comparably high cocaine > pleasant LPP responsiveness in the urine-negative 10R/10R subjects (e.g., does it replicate with fMRI when the comparison is with pleasant stimuli?). (3) As alluded to above, the mechanisms underlying positive cocaine urine screens remain to be clarified, including dissociating withdrawal from priming. In addition, acute withdrawal per se needs to be dissociated from the stress associated with cocaine discontinuation. Future studies could assess whether carriers of the 9R-allele also show enhanced reactivity while undergoing stress that is independent from withdrawal. (4) Current sample size constraints did not allow for data partitioning by gender or ethnicity. In light of previous studies that have uncovered a significant relationship between 3′-UTR VNTR genotype and DAT density in the striatum for whites but not for blacks (Shumay et al., 2011), future studies of this kind could recruit equal numbers of individuals from both races. Another future recruitment goal should be to include more 9R-allele carriers with drug-negative urines and/or longer current abstinence periods, as the size of this subgroup was relatively small (N = 10). Although a notable strength of our approach is that our central CUD+ > CUD− contrast in 9R-allele subjects was observed across multiple methodological approaches, and although boxplots did not reveal any large outliers that could potentially drive our findings, results with these 9R-allele CUD− should be interpreted with some caution until larger sample sizes are accrued for this subgroup. It could also be interesting to recruit homozygote carriers of the 9R/9R genotype, which were scarce in our sample (N = 8 across all available subjects). Recruitment of more 9R/9R subjects would be interesting in light of other research that has observed diminished response in the 9R/9R genotype to amphetamine in healthy adults (Lott et al., 2005) and methylphenidate in children with attention deficit hyperactivity disorder (Stein et al., 2005); there is also a seeming contradiction between these studies that suggests blunted responsiveness to DA agonists in 9R/9R subjects (Lott et al., 2005; Stein et al., 2005) and another study suggesting heightened responsiveness to methylphenidate in 9R-allele carriers with binge eating disorder (Davis et al., 2007), which indeed merits further investigation.
In conclusion, results of this study help to elucidate the DAT gene's modulation of behavioral and neural responses to drug cues in addiction, suggesting that such responsiveness depends not only on the trait (gene) but also the state (abstinence) of the individual. Identifying such a specific, well-defined phenotype has the potential to illuminate mechanistic, biologically informed pathways of risk for psychopathology (Sweitzer et al., 2012). Because sensitivity to drug cues is associated with a more severe addiction phenotype (Volkow et al., 2006), cocaine-addicted 9R-allele carriers could possibly benefit from additional therapeutic intervention to help regulate reactivity to drug-associated stimuli, especially during the early stages of treatment, enabling more appropriate and efficient allocation of scarce clinical resources and improved clinical outcomes. In particular, this subgroup of addicted individuals could be targeted for specific treatments and/or medications based on their genetic profiles (an emerging field of study subsumed under the name “pharmacogenetics”) (Haile et al., 2009), which has shown initial promise in moderating the effects of naltrexone on drinking behavior (Anton et al., 2012) and cue-reactivity (Schacht et al., 2012) in alcohol dependence.
Footnotes
This work was supported by National Institute on Drug Abuse Grant 1R01DA023579 to R.Z.G., Grant 1F32DA030017–01 to S.J.M., and Grant 1F32DA033088 to M.A.P. We thank Michail Misyrlis, Thomas Maloney, Patricia A. Woicik, Dardo Tomasi, Ruiliang Wang, and Gene-Jack Wang for assistance.
The authors declare no competing financial interests.
- Correspondence should be addressed to Dr. Rita Z. Goldstein, One Gustave L. Levy Place, Box 1230, New York, NY 10029-6574. rita.goldstein{at}mssm.edu