|Home | About | Journals | Submit | Contact Us | Français|
The serotonergic neurotoxin, 3,4-methylenedioxymethamphetamine (MDMA/Ecstasy), is a highly popular recreational drug. Human recreational MDMA users have neurocognitive and neuropsychiatric impairments, and human neuroimaging data are consistent with animal reports of serotonin neurotoxicity. However, functional neuroimaging studies have not found consistent effects of MDMA on brain neurophysiology in human users. Several lines of evidence suggest that studying MDMA effects in visual system might reveal the general cortical and subcortical neurophysiological consequences of MDMA use. We used 3T functional magnetic resonance imaging during visual stimulation to compare visual system lateral geniculate nucleus (LGN) and Brodmann Area (BA) 17 and BA 18 activation in 20 long abstinent (479.95±580.65 days) MDMA users and 20 non-MDMA user controls. Lifetime quantity of MDMA use was strongly positively correlated with blood oxygenation level-dependent (BOLD) signal intensity in bilateral LGN (rs=0.59; p=0.007), BA 17 (rs=0.50; p=0.027), and BA 18 (rs=0.48; p=0.031), and with the spatial extent of activation in BA 17 (rs=0.059; p=0.007) and BA 18 (rs=0.55; p=0.013). There were no between-group differences in brain activation in any region, but the heaviest MDMA users showed a significantly greater spatial extent of activation than controls in BA 17 (p=0.031) and BA 18 (p=0.049). These results suggest that human recreational MDMA use may be associated with a long-lasting increase in cortical excitability, possibly through loss of serotonin input to cortical and subcortical regions. When considered in the context of previous results, cortical hyper-excitability may be a biomarker for MDMA-induced serotonin neurotoxicity.
The substituted amphetamine, 3,4-methylenedioxymethamphetamine (MDMA/Ecstasy), is widely used as a recreational drug owing to its acute stimulant, hallucinogenic, and pro-social effects that are produced by MDMA's acute actions on brain serotonin (5-HT), dopamine, and norepinephrine (Bankson and Cunningham, 2001). Despite evidence that MDMA causes long-lasting serotonergic neurotoxicity, MDMA use remains high in the United States and European Union (NSDUH, 2009; EMCDDA, 2010). Because of MDMA's effects on neurotransmitter physiology (Yamamoto et al, 2010) and case reports suggesting that it may facilitate insight and engagement (Liester et al, 1992), MDMA is in clinical research trials as an adjunct to psychotherapy for anxiety disorders (NIH, 2008a, 2008b), opening the possibility that MDMA may become an approved pharmacotherapeutic agent.
In this context of widespread recreational exposure and potential therapeutic administration, much remains unknown with regard to the long-lasting effects of MDMA on the human brain. In animals, long-lasting neurotoxicity of MDMA is manifested as a dose-related loss of 5-HT axons or as sustained reductions in 5-HT neuronal signaling (Capela et al, 2009; Wilson et al, 1989; Fischer et al, 1995; Hatzidimitriou et al, 1999; Green et al, 2003; Lyles and Cadet, 2003), although evidence for axon loss is not found in all studies (Baumann et al, 2007; Wang et al, 2007; Fantegrossi et al, 2004). Human MDMA users have 5-HT alterations consistent with the animal neurotoxicity data that include reduced 5-HT transporter binding (McCann et al, 2008; Kish et al, 2010), increased 5-HT2A receptors (Reneman et al, 2002), and reduced 5-HT metabolite levels in cerebrospinal fluid (Ricaurte et al, 1990). Recreational MDMA users also have chronic neurocognitive (Kalechstein et al, 2007) and neuropsychiatric (Karlsen et al, 2008) impairments.
Despite this evidence for long-lasting brain effects, functional neuroimaging studies of human MDMA users have not found a consistent pattern of altered brain neurophysiology in association with MDMA exposure (Cowan, 2007; Cowan et al, 2008a; Daumann et al, 2003, 2004, 2005; Jager et al, 2008a, 2008b; Jacobsen et al, 2004; Moeller et al, 2004; Roberts et al, 2009; Raj et al, 2009). Several factors, including dose–response effects, polydrug effects, and task design, may account for the absence of a consistent neurophysiological effect of MDMA exposure in extant functional neuroimaging studies.
Several lines of evidence suggest that the visual system may be ideal as a paradigm system for revealing chronic MDMA effects. First, 5-HT axons densely innervate structures in the visual pathway, including the lateral geniculate nucleus (LGN) of the thalamus and the occipital cortex (Waterhouse et al, 1990), which contains the primary and secondary visual cortical regions (Tork, 1990). Second, MDMA administration produces long-lasting or permanent loss of 5-HT axon labeling in primate occipital cortex (Hatzidimitriou et al, 1999), which is prominent in distal occipital regions. Third, 5-HT influences visual cortical excitatory and inhibitory tone (Moreau et al, 2010), suggesting that MDMA-induced loss or reduction in 5-HT signaling would be expected to alter visual cortical excitability. Fourth, human MDMA use is associated with reduced occipital cortex 5-HT reuptake transporter (eg, McCann et al, 1998, 2005; Semple et al, 1999; Kish et al, 2010), increased 5-HT2A receptor binding (Reneman et al, 2002; reviewed by Cowan, 2007), and altered visual activation (Cowan et al, 2006). Fifth, MDMA users have altered visual perception (Brown et al, 2007; Oliveri and Calvo, 2003).
Given the lack of clarity in the human literature with regard to the persistent neurophysiological effects of MDMA on brain function, we designed a study to address unresolved questions with regard to MDMA's effects in human recreational users. We used visual stimulation in conjunction with the functional magnetic resonance imaging (fMRI) blood oxygenation level-dependent (BOLD) method to examine visual stimulus-evoked activation in the LGN and visual cortical subregions in Brodmann Area (BA) 17 and BA 18 in abstinent MDMA users and non-MDMA user control subjects. On the basis of previous MDMA research and the known neurophysiological role of 5-HT in visual regions, our primary hypothesis was that greater lifetime MDMA exposure would predict greater visual system stimulus-evoked signal intensity and spatial extent of activation that is not due to other drug effects and that does not resolve with abstinence. If these hypotheses are supported, the results would suggest that MDMA-associated toxicity leads to an increase in cortical excitability, likely due in part to altered 5-HT function.
This study used a mixed, within-subject and between-groups design to assess the association of MDMA use with brain activation in the LGN of the thalamus and BA 17 and BA 18 of cerebral occipital cortex during visual stimulation in abstinent MDMA users and non-MDMA (control) subjects.
Forty-five subjects took part in the visual stimulation paradigm: data from five subjects were excluded for technical failure. Research subjects were men or women, 18–35 years old, and reported previous MDMA use (MDMA users group; N=20) or other recreational drug or alcohol use (non-MDMA user controls; N=20). Subjects provided written informed consent and the protocol was approved by the Vanderbilt University Institutional Review Board and conformed to principles of the World Medical Association's Declaration of Helsinki.
Subjects enrolled in this study took part in a multimodal imaging study of MDMA effects, but did not overlap with those in our previous visual system study (Cowan et al, 2006). Some participated in our study of motor function (Karageorgiou et al, 2009). Subjects were screened with: clinician version of the Structured Clinical Interview for DSM-IV Axis I Disorders (First et al, 1997) or Mini-International Neuropsychiatric Interview (Sheehan et al, 1998) urine drug testing (Triage Drugs of Abuse Panel, Biosite Diagnostics, San Diego, CA) and, if female, pregnancy screening (Sure-Vue Urine hCG, Fisher HealthCare, Houston, TX). Subjects were re-tested for drugs and pregnancy, and also passed an alcohol breathalyzer test (Alco Sensor III, Intoximeters, St Louis, MO) to confirm recent abstinence (see Abstinence Criteria) on each study day.
Subjects reported abstinence of at least 14 days for drug use and 48h for alcohol use. Exclusions were as follows: general medical conditions, contraindications to MRI scanning, pregnancy, lifetime history of Axis-I psychiatric diagnoses, except for substance use/abuse/dependence-related disorders, current substance dependence (other than nicotine and caffeine), positive urine drug or alcohol screen, use of psychoactive or vasoactive medications within 6 weeks of enrollment, and head injury with loss of consciousness greater than 30min.
We assessed the extent of self-reported substance use with a drug use history questionnaire, utilizing a time-line follow-back method (Fals-Stewart et al, 2000). The questionnaire contains items to indicate when a specific drug was last used and assays recent and remote exposures as number of episodes and approximate amount of the drug used (Cowan et al, 2006, 2007). An episode of drug use was defined as use of a form of a drug (eg, snorting, injecting) within a single 24-h period. We relied on lifetime measures of drug use for our analyses because: (1) our primary aim was to assess for chronic cumulative MDMA effects; (2) the average duration of MDMA abstinence was over 1 year (therefore many subjects did not have past year or past month use); and (3) our previous studies have found that lifetime MDMA use is most strongly associated with brain activation outcomes.
We used a modified two-color visual stimulation paradigm that has proven sensitive in our previous studies to the effects of gender (Cowan et al, 2000), MDMA use (Cowan et al, 2006), and amphetamine administration (Cowan et al, 2008b) on occipital cortical BOLD fMRI signal. Subjects viewed red and blue light stimuli at a flash rate of 8Hz at three intensity levels (low, medium, and high, based on the dynamic range of the computer display) with within-run randomization of color and intensity conditions. We chose a block design to maximize our ability to detect activation differences (Birn et al, 2002). We used three intensity levels to avoid floor or ceiling effects on activation. Each run contained four repetitions of each color and intensity. The duration of each light block was 10s and each baseline period (black screen) was 14s. Visual stimuli were generated using Matlab (Mathworks, Natick, MA) on a Macintosh computer. Stimuli were displayed using either liquid crystal diode (Resonance Technologies) goggles (for studies performed on the GE scanner) or onto a screen behind the subject's head using an Avotec projector and were viewed using a mirror placed within the bore of the scanner (for studies performed on the Philips scanner).
Thirty-one subjects were scanned with a Philips 3T Intera Achieva MRI Scanner (Philips Medical Systems, Andover, MA) and nine subjects with a 3T GE Signa whole body scanner (GE Healthcare, Milwaukee, WI). The following parameters were common to both scanners: 288 collected dynamics within a functional run between 582 and 592s, 19 or 28 slices either 4.5 or 5mm thick with a gap from 0.00 to 0.45mm; TR=2s, TE=35ms; 70° or 79° flip angle; and FOV=240 × 240, matrix=128 × 128 pixels. Functional runs were co-registered to T1-weighted 3D anatomical images, acquired under the following parameters: 128 or 170 slices, 1.00 or 1.20mm slice thickness, gap=0mm; TE=3.7–4.6ms, TR=7.9–8.9ms; 2°–8° flip angle; and matrix=256 × 256 pixels, FOV=240 × 170mm, or 256 × 170mm. We assessed for potential differential scanner/delivery method effects using multiple approaches (see Online Supplement, Supplementary Table S1, and Supplementary Figure S1).
Images were preprocessed and analyzed using in-house software and SPM5 (Wellcome Department of Cognitive Neuroscience, London, UK). Functional and anatomical data were co-registered and images were normalized to the SPM T1 template image, and then smoothed with a full-width half maximum 8-mm Gaussian kernel.
We used a mass univariate approach as implemented in SPM5 to create individual voxel-level t-maps for stimulus-induced activation for each stimulus condition (red or blue color at low, medium, or high intensity) for each subject. T-maps generated in the first-level analysis were used for second-level analyses. In addition, the stimulus-evoked signal intensity data was extracted for additional analyses outside SPM5.
On the basis of the known anatomy of the visual system and on findings from our earlier studies, we examined stimulus-evoked activation in six a priori ROIs (right and left LGN; right and left BA 17; and right and left BA 18). We used the Wake Forest University Pickatlas tool (v.2.4) in SPM (Maldjian et al, 2003; Tzourio-Mazoyer et al, 2002) to dilate the LGN mask by one voxel in all dimensions to account for potential individual variations in this small structure. We transposed these six bilateral ROIs into cluster images using the MarsBaR ‘build' function. We used Monte Carlo simulation as implemented in AlphaSim (AFNI, NIMH) to determine a corrected statistical threshold for activation to visual stimulation in the a priori ROIs. To avoid partial volume effects arising from averaging physiologically non-relevant brain regions within the full anatomical BA 17 and BA 18 ROIs with those activated by the visual stimulus, we defined the stimulus-relevant portion of the ROI for BA 17 and BA 18 as those voxels activated during combined red and blue high stimulus intensity across both MDMA users and controls at a voxel-level p<0.05, with an extent threshold of 90 voxels (as determined using AlphaSim) for family-wise corrected p<0.05. We used these physiological ROI masks for bilateral BA 17 and BA 18 for subsequent analyses of signal intensity and spatial extent of activated voxels by extracting this information and exporting individual results to SPSS for further analysis. We used a threshold of p<0.001 uncorrected for multiple comparisons in the small volumes to determine the spatial extent of activated voxels in the physiological ROIs within BA 17 and BA 18. Because there were no activated voxels detected in the LGN using the corrected thresholds that were applied to measure BA 17 and BA 18 activation, we did not analyze the spatial extent of activation within the LGN and extracted only the signal intensity measure. We used an in-house MATLAB script to extract ROI data on a voxel level (to determine the proportion of activated voxels to generate the spatial extent measure) and MarsBaR to extract the signal intensity data.
Descriptive statistics were used to summarize and initially inspect the distributions of demographic and study measures. Because MDMA and other drug use were non-normally distributed, we employed non-parametric approaches for within- and between-group comparisons for activation outcome measures.
Within-group analysis. Because our earlier fMRI studies (Cowan et al, 2006; Karageorgiou et al, 2009; Raj et al, 2009) revealed that lifetime MDMA exposure correlated with task-evoked brain activation across cohorts, tasks, and brain regions, our primary hypothesis was that lifetime MDMA exposure within the MDMA user group would predict visual system activation. Therefore, our primary analytic approach was the assessment of the relationships of stimulus-evoked signal changes with lifetime MDMA use (Spearman's rank correlations in SPSS; Spearman's rho; rs). We examined activation effects of the visual stimulation paradigm by measuring stimulus-evoked BOLD signal intensity and the spatial extent of stimulus-induced activation.
Between-group analysis. We compared between-group activation in the six a priori anatomical ROIs (right and left LGN, BA 17, and BA 18) within SPM using the contrast of red and blue high-intensity stimulus vs baseline. In addition, we used in-house Matlab scripts to extract the ROI data for stimulus-evoked BOLD signal change and compared mean signal intensity data between groups using Mann–Whitney tests in SPSS.
We examined patterns of activation from the second-level analyses for each color and stimulus intensity. Because the overall intensity of the average signal change was small and there were no clearly observed color or hemisphere-specific differences in the activation patterns, we chose initially to aggregate activation outcome data from the right and left hemispheres and for red and blue color. An analysis of the relationship of stimulus intensity from the aggregated data (low, medium, or high) with BOLD signal intensity in LGN, BA 17, and BA 18 revealed that BOLD signal intensity increased with increasing stimulus intensity with no evidence for a ceiling effect (Figure 1). The second-level analyses revealed also that the most extensive activation was produced by the combination of red and blue high-intensity stimuli. Therefore, all subsequent analyses were conducted using activation data (as BOLD signal intensity or spatial extent of activation) averaging signal obtained in red and blue high stimulus intensity conditions and from right and left hemispheres. Figure 2 shows the resulting pattern of activation in LGN, BA 17, and BA 18 following high-intensity visual stimulation. Activated voxels were not detected in the LGN using corrected statistical thresholds, but activation was robust in BA 17 and BA 18. Therefore (as noted above), signal intensity was used as the outcome measure for analyses in the LGN, whereas both signal intensity and spatial extent of activation were used as outcomes for BA 17 and BA 18.
Sample characteristics are summarized in Table 1. No statistically significant differences in age or gender distributions were found between the MDMA and non-MDMA user groups (p>0.05, Mann–Whitney and χ2 tests). No statistically significant associations of age or sex with any activation measure within any region were found (p>0.05, Spearman's rank correlations). Additional sample characteristics are provided in Online Supplement in Section ‘Additional Sample Characteristics'.
Effect of MDMA. Activation produced by the high-intensity stimulation is shown in Figure 2. The major finding from the within-subject correlation analysis was that greater lifetime MDMA use was statistically significantly associated with greater stimulus-evoked visual activation for signal intensity in LGN, BA 17, and BA 18 (Table 2; Figure 3), and for spatial extent of activation in BA 17 and BA 18 (Table 2) (p>0.05, Spearman's Rank correlations). There was no statistically significant association of the duration of MDMA abstinence with any outcome measure for any ROI.
Effect of combined MDMA and methamphetamine. Of the 20 MDMA subjects, seven also reported use of methamphetamine (METH). Because METH produces serotonergic toxicity similar to that caused by MDMA (Peat et al, 1983; Woolverton et al, 1989), we reasoned that both drugs might produce similar and therefore additive effects on brain activation, especially if the observed findings are due to altered 5-HT neurophysiology. To explore this possibility, we created a measure of substituted amphetamine exposure that combined MDMA and METH exposures. As we did not have a priori data to suggest the mg equivalent 5-HT toxicity produced by each drug in human beings, we simply summed the lifetime quantity of MDMA and METH use (in mg), shown as MDMA+METH in Table 2. Overall, the composite measure accounted for 32% of the variance (as the square of rs) for LGN signal intensity vs 35% for MDMA alone. For BA 17 signal intensity, the composite measure accounted for 35, vs 25% of the variance for MDMA alone, with similar results found for BA 17 spatial extent, and for both signal intensity and spatial extent in BA 18. As for MDMA, there was no statistically significant association of the duration of METH abstinence with any activation measure in the METH-exposed subgroup.
Other drug effects. Although our primary analysis revealed that greater lifetime MDMA use, and greater lifetime use of MDMA+METH was associated with greater activation in all visual regions, MDMA users in our cohort (consistent with worldwide use patterns) reported exposure to multiple recreational drugs other than MDMA. We therefore examined the data for the influence of other drugs. Because partial correlations require data for all three variables in a specific analysis (activation measure, MDMA use, specific other drug use) and to maintain a sample of size of at least 10, these analyses were conducted only for those drugs used by at least 50% of the MDMA cohort (Table 1; alcohol, cannabis, cocaine, codeine, LSD, opium, psilocybin, and sedatives). There were no statistically significant associations of lifetime quantity of drug use with activation outcome variables (Spearman's rank correlation; p>0.05).
As noted in the Materials and Methods, we acquired data on two different scanners, each employing a different stimulus delivery method (goggle vs projection screen). It is theoretically possible that the use of the different scanner/stimulus delivery methods might have affected our results. We conducted an analysis of the potential effects of these differences on our findings.
Lifetime quantity of MDMA use by scanner/stimulus method. MDMA use levels showed mean differences and correlations with scanner type at a trend (p<0.10 level). Mean MDMA use was slightly greater for MDMA users studied on the GE scanner (p=0.076; Mann–Whitney) (Supplementary Table S1). Scanner type was nonsignificantly correlated with lifetime quantity of MDMA use in mg (rs=0.41, p=0.074; Spearman's rank correlation).
Relationship of scanner to correlation of lifetime MDMA use with signal intensity. Supplementary Figure S1 depicts scatter plots of data from GE and Philips scanners to show the degree of overlap and pattern of data from the two scanners. The MDMA users studied on the GE scanner had higher minimum use levels than the Philips group, but there was good overlap for use levels in the higher use ranges.
Partial correlation analysis controlling for scanner effects. Although the graphic analysis suggests no scanner/stimulus-specific differences in slope of the relationship, we conducted a statistical analysis of the data adjusting for scanner/stimulus condition. We anticipated that as MDMA use was correlated at a trend level with scanner/stimulus method and as the additional scanner covariate would reduce the available degrees of freedom, a partial correlation analysis adjusted for scanner/stimulus would yield a weaker relationship between MDMA use and regional signal intensity. However, despite these factors, the partial correlation of lifetime quantity of MDMA use with regional visual stimulus-evoked signal intensity remained positive, but significantly so only for LGN (LGN, r=0.60, p=0.006; BA 17, r=0.35, p=0.092; and BA 18, r=0.351, p=0.141).
As predicted by data from our earlier study, there was general overlap between control group and MDMA user group activation measures, and there were no statistically significant differences between the MDMA and control groups in average signal intensity within the LGN, BA 17, or BA 18 (Table 3; Figure 4). In addition, there were no statistically significant between-group differences in the spatial extent of activation for BA 17 and BA 18 (Table 3; Figure 5). As noted above, spatial extent analysis was not performed for the LGN.
Given the relationship between lifetime MDMA exposure and activation shown above for the within-group analysis, we anticipated that the lack of between-group differences may have been secondary to the fact that the MDMA user group was not homogenous—that is, individuals with greater MDMA exposure would have higher mean activation measures than those with lower exposure. To address this issue, we conducted a subgroup analysis using a median split procedure to divide the MDMA user group into high- and low-exposure groups. To aid in controlling for Type I error, we used a non-parametric ANOVA (Kruskal–Wallis), followed by a post hoc Mann–Whitney test. This analysis revealed that the high-exposure MDMA group had significantly greater spatial extent of activation when compared with the control group, but there were no significant between-group differences in signal intensity (Table 3). In contrast to the greater spatial extent of activation seen for the comparison of high-exposure MDMA users to the control group, the only significant difference between low-exposure MDMA users and control subjects was found in the LGN where the low-exposure MDMA subgroup had lower signal intensity than the controls (Table 3). Only the signal intensity difference in LGN met the adjusted α-level (p<0.017) after accounting for the number of comparisons.
The primary finding of this study is that greater lifetime MDMA exposure was associated with greater stimulus-evoked visual system activation and this relationship was not affected by the duration of MDMA abstinence. Greater cumulative lifetime MDMA use was associated with greater stimulus-evoked visual system BOLD signal intensity in the LGN, BA 17, and BA 18, and with greater spatial extent of activated voxels in BA 17 and BA 18. Although there were no between-group differences in activation outcome measures when the entire MDMA group was compared with the non-MDMA-exposed control group, subgroup analyses suggested that high-exposure MDMA users tended to have greater activation than controls, most notably for the spatial extent of activation. Additive effects of MDMA and METH on both signal intensity and spatial extent of activation in BA 17 and BA 18 were also observed, supporting the possibility that the 5-HT neurotoxicity common to both drugs may account in part for the observed results. However, this finding needs to be interpreted cautiously until replicated. The long average duration of self-reported MDMA abstinence and the positive exposure–activation relationship, despite a long average duration of abstinence, suggests that MDMA use may cause graded, long-lasting, and potentially permanent alterations in human visual system neurophysiology. The detected exposure-related increases in task- or stimulus-evoked activations may be a neurophysiological biomarker for MDMA neurotoxicity.
As reviewed previously (Cowan, 2007; Cowan et al, 2008a), previous studies employing more complex cognitive paradigms found a range of activation patterns following chronic MDMA exposure (Daumann et al, 2003, 2004, 2005; Jager et al, 2008a, 2008b; Jacobsen et al, 2004; Moeller et al, 2004). Subsequent studies have revealed a similarly mixed picture when using fMRI to examine brain activation during cognitive paradigms (Roberts et al, 2009). Our own recent report of semantic processing showed a largely inverse association of lifetime MDMA use and signal intensity, but this inverse association was found after a complex subtraction analysis of word and non-word tasks in regions weakly activated by the employed paradigm (Raj et al, 2009). Several behavioral and pharmacological aspects of recreational MDMA use likely account for the lack of consistent findings with regard to MDMA effects on brain neurophysiology in previous reports. First, MDMA toxicity in animals is dose-dependent (Green et al, 2003), suggesting that human MDMA imaging research study design might need to account for individual variation in drug exposure. Second, MDMA users worldwide have high levels of exposure to other drugs (Wish et al, 2006), which may thwart attempts to specifically associate MDMA use with altered brain function. These polydrug effects may be accounted for in part by comparison with control groups having other drug (but no MDMA) exposure and several investigators have attempted to include drug-matched controls in their study design. However, matching average use across groups does not account fully for individual patterns in drug exposure and also controls only for the most commonly used substances (e.g., cannabis). Third, most studies investigating MDMA effects on human neurophysiology have employed paradigms that involve the performance of complex cognitive tasks during functional neuroimaging (reviewed by Cowan, 2007). However, interpreting functional neuroimaging findings arising from complex cognitive tasks in terms of underlying brain neurophysiology is challenging. Cognitive processing is more widely distributed in the brain than simple sensory or motor function, and brain activation results from cognitive paradigms are often reported as contrasts of tasks, potentially obscuring relationships to more basic neurophysiology (reviewed by Cowan, 2007). Our experimental approach was specifically designed to attempt to account for these factors.
In this study and in our earlier report of visual system function in MDMA users (Cowan et al, 2006), we did not find between-group differences in average stimulus-evoked activation when MDMA users were compared with non-MDMA-exposed control subjects. However, in both our initial and current reports, we found that the MDMA user subgroup having the greatest lifetime MDMA exposure showed a greater spatial extent of activation than non-MDMA-exposed controls. This suggests that studies restricting enrollment to individuals with heavier MDMA use would be more likely to find greater activation in the MDMA user group, whereas studies enrolling MDMA users with lower MDMA exposure or a broad range of MDMA exposure might find opposite or no between-group differences. It is also possible that the effect observed here for the visual system is not present in all cortical regions, even for simple sensory tasks. However, an evoked potential study of 5-HT-related auditory processing found, similarly to our fMRI findings, that lifetime MDMA use was associated with greater amplitude of cortical response to auditory stimuli (Croft et al, 2001).
The absence of strong between-group differences in activation does not suggest that MDMA is safe or that there are no differences in brain function in MDMA users and controls. The increased spatial extent of activation in the higher-use MDMA group suggests that between-group differences emerge with greater MDMA exposure levels, which is consistent with the strong exposure–response effect shown in the within-group analysis for MDMA. Two additional factors may account for the lack of strong differences in the between-group comparison. First, MDMA has multiple mechanisms of potential neurotoxicity. Each of these may manifest differently in terms of neurophysiological effects as seen in the fMRI signal, with effects of 5-HT axon loss dominating at higher MDMA dosages. Second, other drugs may also be associated with reductions in the fMRI signal (Raj et al, 2009) and the effects of 5-HT axon loss may therefore be obscured in MDMA users with very low MDMA exposures. This is a pattern consistent with the clear exposure–response relationship shown within the MDMA user group.
Our study used the fMRI BOLD method to assay for MDMA effects on brain physiology. The fMRI BOLD method indirectly assays brain neurophysiology because increased focal neuronal synaptic or metabolic activity is associated with increased local blood flow (reviewed by Ekstrom, 2010). Increased task- or stimulus-evoked regional blood flow above a defined baseline is considered ‘activation' of a brain region and activation can be reported as the magnitude of BOLD signal intensity increase above baseline (a continuous measure of the signal change) or as the number, volume, or spatial extent of activated voxels (a discrete measure of voxels in which signal change exceeds a defined statistical threshold). Using these metrics, our initial report (Cowan et al, 2006), in a cohort not overlapping with the current group, suggested that greater MDMA use was associated with greater spatial extent of occipital activation. However, this early study was not designed to account for cumulative MDMA dosing or to control for the duration of MDMA abstinence or for polydrug effects, and methodological limitations of the scanner technology prevented the isolation of effects in specific anatomical regions. These results replicate the initial finding with regard to spatial activation, but not with regard to signal intensity. We suspected that the failure to detect an association of MDMA use and signal intensity in the original study was likely due to partial volume effects of averaging a small area of physiologically relevant task-evoked signal across a broader anatomically relevant region of interest. We extend that finding here by: (1) revealing that greater MDMA use is associated with not only greater spatial extent of activation, but also with greater signal intensity, and (2) by showing alterations in visual thalamus (LGN) and specific occipital cortical subregions, BA 17 and BA 18. These visual pathway findings, when coupled with evidence from our recent fMRI report of MDMA effects during motor task performance (Karageorgiou et al, 2009), show that simple sensory or motor fMRI paradigms are sensitive assays of potential serotonergic neurotoxicity associated with MDMA use.
We found that an aggregate measure combining lifetime MDMA exposure with lifetime METH exposure accounted for a greater proportion of the variance in activation in visual cortical regions (BA 17 and BA 18) than for MDMA exposure alone. Although METH effects on a simple visual fMRI paradigm have not been previously reported to our knowledge, several studies of brain activation during cognitive paradigms in METH-dependent individuals have found a mixed pattern of activation, similar to that observed in MDMA users, when compared with control subjects during a facial affect matching task (Payer et al, 2008) and a delay discounting task (Monterosso et al, 2007; Hoffman et al, 2008).
Associations of MDMA use with altered activation measures in these data do not prove that MDMA exposure altered 5-HT function or caused 5-HT axon loss. However, several indirect lines of evidence suggest that serotonergic toxicity likely explains in part the observed findings. First, animal models have repeatedly revealed that MDMA exerts a dose-related toxicity to serotonergic axons (reviewed by Green et al, 2003) and our data show a similar dose–response effect on activation measures. Second, although MDMA produces neuroplastic changes as altered spines or dendrite shape in medial prefrontal cortex and nucleus accumbens in rats (Ball et al, 2009), and may induce scattered neuronal cell death (Schmued, 2003), we are not aware of evidence to suggest that these effects are prominent in visual system regions. Third, as discussed in the section that follows, 5-HT plays a major role in excitation/inhibition balance in the cortex (Moreau et al, 2010), and specifically with respect to cortical pyramidal response to thalamic inputs in visual system where serotonin limits the lateral spread of incoming excitation. Finally, MDMA and METH produce similar effects on 5-HT axon markers (eg, Peat et al, 1983; Woolverton et al, 1989; reviewed by Quinton and Yamamoto, 2006), and the aggregate MDMA+METH measure explained a greater portion of the variance in activation in both visual cortical regions.
Several lines of evidence converge to suggest that the most prominent effect of MDMA-induced 5-HT axon loss may be cortical hyper-excitability. We have previously speculated that (Karageorgiou et al, 2009), as have others (Jacobson et al, 2005; Vannini et al, 2007), increased activation in the presence of normal task performance is consistent with reduced cortical efficiency. If interpreted in terms of the amount of signal change or the spatial extent in cortical or subcortical regions activated by the visual stimulus, these and earlier findings (Cowan et al, 2006; Karageorgiou et al, 2009) would be consistent with less efficient brain functioning in association with MDMA use, in that more neural work (as indirectly indexed by increased signal and therefore potential metabolic demand) and more neural volume (as indexed by the increased spatial extent of spread) are utilized for sensory response. In line with our initial visual study and this report, we also recently reported that a motor task produced greater activation in cortico-striato-thalamic loops as lifetime MDMA use increased (Karageorgiou et al, 2009). In aggregate, these studies suggest that greater task- or stimulus-evoked activation with greater MDMA exposure may be a general neurophysiological and dose-related consequence of recreational MDMA use in human beings, especially when assayed with relatively simple visual and motor activation paradigms. We now propose that the neural basis for the observed reduction in cortical efficiency (measured as increased activation) in association with MDMA use is due in part to an increase in cortical excitability that is secondary to MDMA-induced reduction in serotonin signaling. As discussed below, evidence consistent with MDMA-associated increases in cortical excitability has been reported from electrophysiological studies in mice and in human visual cortex. These results provide further neurophysiological support for MDMA-induced cortical, and possibly thalamic, hyper-excitability.
How might greater MDMA exposure lead to task- or stimulus-evoked cortical hyper-excitability? The origins of the fMRI BOLD signal are not fully understood, but the bulk of what is measured as ‘activation' in fMRI studies is thought to reflect indirectly graded synaptic inputs to a brain region and their consequences (measured as the local field potential) and local neuronal processing, determined in part by glutamate release (Logothetis, 2008; Bennett et al, 2008). However, the precise neural and metabolic events contributing to the BOLD signal are complex and these likely vary with local neural circuitry (reviewed by Ekstrom, 2010). Factoring in visual system neural circuitry and current knowledge of the BOLD signal in visual cortex, we have developed a tentative, highly simplified, neural model (Figure 6) to help explain the positive correlation of MDMA exposure and activation measures in visual system and as reflective of the net effect of MDMA exposure on LGN and cortical output (pyramidal) neurons. 5-HT influences retino-geniculate synaptic transmission in a frequency-dependent manner through multiple receptors and through pre- and post-synaptic mechanisms (Seeburg et al, 2004). Dorsal raphe (the source for the fine diameter 5-HT axons most vulnerable to substituted amphetamines) stimulation in rats leads to the inhibition of LGN response to retinal stimulation and decreased thalamo-cortical excitation (Yoshida et al, 1984). Therefore, loss or reduction of dorsal raphe 5-HT input to LGN may result in increased stimulus-evoked thalamo-cortical excitation. In neocortex, 5-HT receptors of multiple subtypes are present on thalamo-cortical afferents, excitatory, and inhibitory neurons. In pyramidal neurons, 5-HT is both excitatory (via 5-HT2A receptors on dendrites and spines) and inhibitory (via 5-HT1A receptors on the axon hillock) (Rauch et al, 2008). 5-HT1A agonists inhibit the lateral spread of cortical excitation, suggesting that loss of agonist signaling at these receptors might increase spatial spread of excitation in the cortex. This model reflects only the potential net effects of MDMA exposure and does not address cellular or receptor level alterations. MDMA has been associated with the loss of retinal ganglion and other retinal cells in animals (Miranda et al, 2007), but we do not include this information in the present model because the functional implications are unclear.
Directly relevant to this study, Oliveri and Calvo (2003) used transcranial magnetic stimulation to elicit visual phosphenes (which they defined as a conscious report of light sensation in the absence of a visual stimulus) in MDMA users and controls and found that the stimulus threshold for generating phosphenes was lower in MDMA users. Giorgi et al (2005) used a mouse model of MDMA administration and found that MDMA increased hippocampus and entorhinal cortex excitability as measured using electroencephalography, a kainate-induced seizure threshold determination paradigm, and glucose uptake. More globally, 5-HT plays a role in setting the gain for inputs to pyramidal neurons (Higgs et al, 2006) and blocks or reduces seizure activity in multiple seizure models (reviewed by Bagdy et al, 2007), further supporting an overall role for this neurotransmitter of reducing cortical excitability. If the BOLD signal largely reflects the consequences of glutamate signaling (Ku et al, 2008; Bennett et al, 2008), our data thus far are consistent with thalamo-cortical, cortical-cortical, and cortico-thalamic increases in glutamate release reflected in increased BOLD signal. As summarized in Figure 6, loss or reduction of dorsal raphe 5-HT input to LGN, BA 17, and BA 18 could result in increased LGN excitability to retinal ganglion cell input, as well as increased local excitatory (via recurrent excitation; Cowan and Wilson, 1994) and feedforward and feedback projections between LGN, BA 17, and BA 18. However, given the complexity of 5-HT innervation and action and the uncertainty of the origin of the BOLD signal, this interpretation is necessarily speculative and is presented mainly to aid understanding and as a framework for future studies aimed at refining this idea. Notably, factors other than MDMA-induced 5-HT neurotoxicity, including neuroplastic changes (Ball et al, 2009), focal neuronal death (Schmued, 2003), altered regional metabolite concentrations or brain gray matter density (Cowan et al, 2003, 2007, 2008a), or other neurotransmitter or neurotoxic effects (Yamamoto et al, 2010), could account in part or in total for the observed findings.
Several factors limit our ability to draw definitive conclusions with regard to the role of MDMA in producing altered brain activation. A cross-sectional study design as employed here cannot assess causation. User self-report of previous MDMA exposure is subject to inaccuracies through recall bias, deception, forgetting, and impurity in Ecstasy pills sold as MDMA. Drug abstinence was verified only for the time period in which drugs or metabolites are detected by urine screens (2–3 days for most drugs), with the result that some subjects could have under-reported recent drug use. However, recent exposure to drugs other than MDMA seems unlikely to have produced the observed exposure–response association with lifetime MDMA use. We grouped data acquired using two scanning/stimulation methods; however, our analysis of the effect of scanner/stimulus on outcomes did not suggest that methodological factors influenced the results. In addition, these findings are consistent with our earlier report using a non-overlapping cohort of MDMA users studied using a 1.5T scanner and different stimulus methods (Cowan et al, 2006), as well as with our recent motor system study performed in a partially overlapping cohort, but on a single scanner (Karageorgiou et al, 2009).
Because our hypotheses predicted that MDMA use would be associated with altered spatial extent of activation and to avoid partial volume effects, we devised an experimental approach that permitted us to isolate the physiologically relevant (task-relevant) sub-region in larger a priori anatomically defined areas (LGN, BA 17, and BA 18). An alternative approach has been used in retinotopic mapping methods to precisely isolate primary (V1) and secondary (V2) visual regions. However, the interpretation of retinotopic mapping, to our knowledge, is ambiguous in the presence of neurophysiological alterations that influence potentially both the probability of declaring a voxel as activated (via effects on increased signal intensity, and thus increased signal-to-noise ratio) and the spatial extent of activation (via effects on the spread of recurrent excitation and inhibition). And, it is not clear that knowing MDMA effects specifically in V1 and V2 would add to our understanding of the general effects of MDMA on cortical neurophysiology.
It is possible that pre-existing brain differences are responsible for the observed results, but we are not aware of a mechanism through which altered visual system excitability would predispose to higher levels of MDMA use. The observed fMRI findings cannot be assumed to imply altered conscious visual perception in this cohort and future studies using biophysical measures in conjunction with fMRI studies are necessary to address this issue. Taffe et al (2001, 2002, 2003) reported that MDMA administration regimens in monkeys that deplete cortical 5-HT by greater than 75% (measured in post-mortem assays) produced detectable behavioral deficits only during tryptophan depletion or during administration of drugs active at 5-HT receptors. These findings suggest that, despite chronic MDMA-induced reductions in brain 5-HT content, compensatory brain mechanisms may sustain brain functions in non-human primates.
However, we think that the significance of this report lies not in detecting visual system dysfunction per se, but more broadly in showing a potential chronic and widespread shift in brain excitability that may be a neural signature of MDMA toxicity. Functional MRI is an indirect method that is influenced both by neuronal and vascular function, and it is therefore possible that the observed effects may not be due to altered neuron to neuron communication. However, we believe this is unlikely because human MDMA use has not been associated with chronically altered cerebral blood flow (Chang et al, 2000). McCann and Ricaurte (2007) and McCann et al (2009a, 2009b) have found that MDMA users have altered sleep, including higher rates of sleep apnea and increased vulnerability of cognitive performance to the effects of sleep deprivation. We did not assay sleep parameters in this cohort, and therefore we cannot rule out the possibility that sleep deprivation may have influenced our results.
These results suggest that recreational MDMA exposure leads to chronic alterations in visual system neurophysiology. These results, and those of others using different techniques, indicate that the observed association of MDMA use with brain activation may be secondary to cortical hyper-excitability upon exposure to a sensory stimulus. Additional studies are needed to replicate this finding, to determine which additional cortical regions show this effect, whether stimulus intensity or task demands influence the observed pattern, and perhaps to show prospectively that MDMA exposure causes altered cortical neurophysiology. If the observed associations of MDMA use with increased activation result from shifts in cortical excitatory thresholds, the implications of this finding and the underlying neural basis for this phenomena need to be more precisely determined. Notably, our study excluded individuals who reported co-morbid psychiatric conditions such as major depression or anxiety disorders. As these disorders are linked to 5-HT function and are reported at higher rates in MDMA users, our study may have excluded a subgroup having greater vulnerability to MDMA toxicity. Additional studies specifically examining MDMA effects in subjects with co-morbid psychiatric conditions appear urgently warranted. Given the many millions of individuals worldwide who have used or continue to use MDMA and ongoing clinical trials investigating MDMA's therapeutic potential, large prospective studies of brain function in abstinent users seem essential to determine whether MDMA-associated effects remain stable, improve, or worsen over time.
This work was supported by NSF (IGERT DGE-0801634 to ALB), NIDA/NIH (R01 DA01537 and R21 DA020149 to RLC, and T32DA021123 to EJC), NIMH/NIH (K01 MH083052 to JUB), and NCRR/NIH (Vanderbilt CTSA Grant UL1 RR024975).
ALB, JUB, EJC JGL, CJC, NDW, AC, TW, CR D I, CC, and RMS report no potential conflicts of interest related to this report. MSD performed consultation or received compensation from the American Economic Association, Belmont University College of Health Sciences, and the University of Oklahoma Health Sciences Center. RLC received publication royalties from Lippincott Williams and Wilkins and consultant income from the Southwest Michigan First Life Science Fund, the University of West Alabama, Novo Nordisk, and Shire Pharmaceuticals.
Supplementary Information accompanies the paper on the Neuropsychopharmacology website (http://www.nature.com/npp)