If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
Address reprint requests to Michelle Hampson, Ph.D., Yale University School of Medicine, Department of Diagnostic Radiology, MRRC, TAC-N121, PO Box 208043, New Haven, CT 06520-8043
Chronic tic disorders are characterized by motor tics that are often preceded by premonitory urges to tic. Functional neuroimaging studies have documented brain activity patterns prior to and during tics, but these studies have not examined whether the activation patterns differ from those seen in normal control subjects performing similar acts.
Methods
A novel method was used to compare brain patterns during tics and intentional movements. First, the part of motor cortex specific to each patient's tic movement was identified. The brain areas activating prior to, during, and after that part of motor cortex during tics were then identified by temporally cross-correlating the time course of the localized motor region with activity in other brain areas. Given that motor cortex was active during tic execution, this yielded information regarding the brain areas active prior to, during, and after the movements. The spatiotemporal pattern of coactivation with motor cortex during tics was contrasted with that seen in healthy control subjects during intentional tic-like movements.
Results
Data from 16 adult subjects with tic disorders and 16 matched control subjects, who performed intentional movements similar to the patients' tics, revealed nearly identical patterns of cross-correlation to motor cortex throughout the brain in the two groups. However, the supplementary motor area showed a significantly broader profile of cross-correlation to motor cortex during tics than during intentional movements.
Conclusions
These findings highlight the importance of the supplementary motor area in tic generation and may point toward novel intervention strategies for individuals suffering with severe tics.
The neurophysiological basis of chronic tic disorders has been investigated using a variety of imaging techniques, including positron emission tomography (PET) (
The functional neuroanatomy of Tourette's syndrome: An FDG-PET study I. Regional changes in cerebral glucose metabolism differentiating patients and controls.
). These studies provide converging evidence that Tourette syndrome (TS) involves aberrant activity in the basal ganglia, midbrain, prefrontal cortex, and paralimbic regions and increased activity in motor cortical regions, including supplementary motor area (SMA). A radioactive water ([150]H20) PET study identified regions activating during tics, including the SMA, premotor and primary motor cortices, anterior cingulate, prefrontal cortex, insula, claustrum, putamen, and caudate nucleus (
), we hypothesized that tics are the result of a chain of neural events that may occur on a time scale slow enough to be examined using current imaging techniques. However, little is known regarding the temporal pattern of brain activity associated with this phenomenon in contrast to voluntary actions. A recent event-related fMRI study investigated brain activity patterns prior to and during tics (
). That study reported activation in several areas prior to tics, including the SMA, insula, and anterior cingulate, that later spread into primary motor cortex. These findings provide our first view of the temporal pattern of brain activity associated with tics, but they fail to address whether this sequence differs from that involved in intentional movements. The aim of this study was to examine the time course of brain activity that unfolds during motor tics and to compare that time course with the pattern of brain activity seen during intentional movements.
Both SMA and the primary somatosensory cortex (S1) are regions that have been implicated in TS and are candidate loci of the urges patients report preceding, and motivating, their tics. Electrical stimulation of the SMA in humans (using subdural electrodes in patients with epilepsy) produces both movements and urges to move reminiscent of the premonitory urges of TS. The nature of the movements (or corresponding urges) range from simple motor acts to complex movements, paralleling the range of simple to complex tics experienced in TS (
). This suggests the SMA may be involved in premonitory urges. Given that premonitory urges often have a sensory aspect, together with recent magnetoencephalographic data implicating S1 in TS (
), S1 is another promising candidate region for these urges. Therefore, we hypothesized that SMA or S1, due to involvement in premonitory urges, would show more activity prior to tics than prior to intentional movements.
Methods and Materials
Human Subjects
Sixteen individuals (13 men) with active tic disorders between 18 and 56 years of age were recruited. Only individuals who reported they could execute their tics without moving their heads and whose repertoire of tics was dominated by one or two specific movements were included. Diagnoses, tic locations, and pertinent medication profiles are presented in Table 1. Thirteen of the 16 patients experienced a bilateral movement as their most frequent tic, thus limiting our ability to investigate laterality issues. Data were obtained regarding the presence of premonitory urges for 14 of the 16 patients, all of whom reported experiencing this phenomenon prior to their tics.
A within-subject design requires accurate measurement of the temporal patterns of brain activity associated with imitated tics in the patients. During any time period in which patients are imitating tics, real tics will also be occurring that are difficult to differentiate from the imitated tics. A previous study that attempted to use a within-subject design for an event-related study of tics was unable to collect viable control data for this reason and because of problematic levels of movement in the patients during tic imitation (
). They also raised concerns regarding the normalcy of voluntary movements in the patient group and recommended that control data be collected in healthy control subjects instead (
For each patient, a healthy subject of the same gender, handedness, and sex and of similar age was recruited. Socioeconomic status (SES) was estimated using the Hollingshead Four Factor Index of Social Status and balanced across groups. This matching procedure resulted in populations with very similar age and SES distributions: patient age: mean 30 (standard deviation 9.7), healthy control subject age: 30 (9.0); patient SES: 52 (11.4), healthy control subject SES: 55 (12.0). All subjects gave informed consent to a protocol reviewed and approved by the Human Investigations Committee of the Yale School of Medicine.
Experimental Paradigm
A standard approach for obtaining time courses of brain activity associated with events requires a record of the times of occurrence of the events (in this case, tics) during an imaging session. As accurate records of tic occurrences during a scanning session can be difficult to obtain, a standard event-related protocol was not adopted in this study. A novel paradigm was utilized that does not require knowledge of the time of occurrence of tics during the imaging session. Instead, it depends upon a priori knowledge concerning a specific brain region that is activated in synchrony with tics. Correlations in activity between that (reference) region and other regions of the brain can then be computed. Assuming the tic events are responsible for a large portion of the variance in the activity level of the reference region, a set of regions that are also activated during the tics will be identified. Correlating to the activity pattern in a brain area known to be involved in events has been shown to produce maps similar to those obtained with a traditional event-related analysis (
The region assumed by this study to be active in synchrony with tics was the portion of motor cortex controlling the musculature specific for the selected tic. This region was localized functionally in each subject and then used as a reference region for correlations during scans in which the subject allowed their tics to occur naturally.
Therefore, two different types of functional scans were collected: localizers to identify the motor cortical reference region of interest (ROI) and scans to be used in the correlation analysis in which patients allowed their tics to occur at a natural rate throughout the scan. The localizers were of block design: during “on” blocks, patients imitated their most frequent tic at a much higher rate than the tics naturally occurred, and during the “off” blocks, patients allowed tics to occur at their natural rate. Although this localizer could be used to identify a set of brain areas active during intentional movements in the patient population, it could not provide a temporal sequence of neurophysiological events associated with intentional movements. As discussed in the Human Subjects subsection above, it is difficult to get a time course of brain activity associated specifically with intentional movements in tic patients. For this reason and because of concerns regarding the normalcy of intentional movements in tic patients (
), the spatiotemporal pattern of brain activity associated with intentional movements was identified separately in a group of healthy control subjects.
For control subjects, there were also two types of functional scans. There were localizing scans of block design in which the on blocks involved imitating the tic movement of the patient to whom they were matched at a high frequency and the off blocks were rest periods. For the remaining scans (used for correlation analyses), the control subjects were instructed to imitate the tic movement of the patient to whom they were matched throughout the scan. They were told they would need to produce the movement in clusters, so that some movements occurred close together in time and others further apart, but such that their total frequency over the scanning session should approximate that of the patient. Both the movement and frequency were practiced prior to the scanning session by the control subject until they reported feeling that they could make the movement easily and with the desired temporal pattern.
Data Acquisition
Subjects were scanned in a Siemens 3T Trio scanner (Siemens Medical Systems, Erlangen, Germany). Each functional scan involved the acquisition of 222 volumes, each composed of 22 4.5-mm axial slices (aligned with the anterior commissure-posterior commissure [AC-PC] line). A T2*-sensitive gradient-recalled, single-shot, echo-planar pulse sequence was used (repetition time [TR] = 1.5 sec, echo time [TE] = 25 msec, flip angle = 80, field of view [FOV] = 202 cm, 642 acquisition matrix).
After the functional data were collected, 22 T1-weighted slices were collected (prescribed in the same locations) for anatomical identification. Immediately prior to the collection of anatomical data, subjects were given a button-box and told to press the box each time they experienced a tic (if they were a patient) or imitated a tic (if they were a control subject). These button-press records were used to estimate movement frequencies during the functional scans. There was no button pressing or reference to that aspect of the experiment during or prior to the functional scans. The estimated frequency of tics in patients had a mean of 7.3 tics/min (standard deviation: 4.7). This was well matched by the control population, which had a mean of 7.7 movements/min (standard deviation: 7.7).
Data Analyses
Aside from motion correction, all other analyses were done using in-house software.
Preprocessing
Data were adjusted for slice acquisition time, motion corrected using the Statistical Parametric Mapping (SPM) algorithm (http://www.fil.ion.ucl.ac.uk/spm/), and spatially smoothed. Pixels with a median value below one fifth of the maximum median pixel value were set to zero.
Data used in correlational analyses were also low-pass filtered (.2 Hz). Prior to correlation, the mean signal from each slice was removed from the time course of activity in each pixel in that slice and from the reference time course. This is similar to the removal of global mean time course but is better at reducing noise caused by transient physiologic fluctuations that occur during part of the volume collection and that can lead to elevated correlations between slices collected next to one another in time.
), images were assigned to “tic forcing” or “resting” blocks, and a pixelwise, two-sample t test was performed. The reference region was defined for each subject as the 30 most activated pixels in primary motor cortex, with no requirement of spatial contiguity.
Correlational Analyses
The reference time course was computed by averaging across pixels in the reference region. The reference time course was correlated with the time course of every pixel within each of the scans used in the correlation analysis. The resulting correlations were transformed to a standard normal distribution (
). These z-transformed data were then averaged across scans to yield a map representing the strength of correlations to motor cortex in terms of standardized z-values.
The reference time course was then shifted one TR earlier and similar computations were performed to produce a map of areas that tended to activate 1.5 seconds before motor cortex. This process was repeated across a range of temporal shifts, resulting in nine maps representing areas activated from 6 sec prior to motor cortex to 6 sec after. The nine z-maps of the correlations from each patient were transformed to Talairach coordinates using a piecewise linear transform.
Group Composites
At each Talairach pixel, the mean across the 16 patients/control subjects was computed and compared with a surrogate distribution of 12,000,000 values (each the mean of 16 random numbers from a standard normal distribution). This surrogate distribution represented the distribution of values expected by chance when averaging random data (i.e., with no effect) across 16 subjects. The fraction of surrogates with a greater magnitude than the computed mean value across patients/control subjects was the p value for that pixel. Composite maps are shown at a Bonferroni-corrected cutoff (adjusted for total number of nonzero pixels) of p < .05.
Primary ROI Analyses
Two cortical ROIs (S1 and SMA) were defined bilaterally in the composite anatomical brain, formed by averaging the Talairach transformed anatomical data from all 32 subjects. These ROIs are shown in Supplement 1.
Values from the Talairach transformed z-maps of correlation for each subject were averaged across all pixels in each ROI and entered into an analysis of variance (ANOVA) with three repeated factors: time point (nine levels), hemisphere (two levels), and group (two levels). Subjects were treated in a paired manner, given the one-on-one matching of age, gender, handedness, movement type, and movement frequency. A time point by group interaction was expected for any region differentially involved in tics as compared with intentional movements. All results were corrected for violations of sphericity using the Greenhouse-Geisser correction. For the SMA, which showed a significant time point by group interaction, post hoc paired t tests were performed at each time point to identify time points that differed significantly between groups.
Exploratory ROI Analyses
Five subcortical/midbrain ROIs (caudate, lenticular nucleus, thalamus, ventral basal ganglia including substantia nigra and subthalamic nucleus, and the midbrain tegmentum including the periaqueductal gray) that have been repeatedly implicated in TS (
The functional neuroanatomy of Tourette's syndrome: An FDG-PET study I. Regional changes in cerebral glucose metabolism differentiating patients and controls.
) were defined bilaterally in the composite brain and analyzed in the same manner as the primary ROIs. These ROIs are also shown in Supplement 1.
Secondary Analyses of the SMA
The difference between groups in the SMA was re-evaluated for the following four subsets of our participants: 1) the 14 patients known to have premonitory urges and their matched control subjects; 2) the 11 patients with no comorbid conditions and their matched control subjects; 3) the 13 patients with a diagnosis of TS (not chronic tic disorder) and their matched control subjects; and 4) the 8 patients who were not taking neuropsychiatric medication and their matched control subjects.
In addition, we conducted an ANOVA on the patient data with three repeated factors: medication status (two levels), time point (nine levels), and hemisphere (two levels). A medication status by time point interaction would provide evidence that the time course of SMA activity in the eight patients taking neuropsychiatric medication differed from the time course of the SMA in the eight patients not taking medication.
Motion Comparisons
For each of the six SPM realignment parameters, the absolute value of the difference in the parameter between each volume and the next was added across all volumes in a given scan. For each subject, these numbers were averaged across scans and a two-sample t test was performed to determine whether the total movement along that dimension differed in patients and control subjects at a p < .05 level.
Results
The maps of cross-correlations to motor cortex during tics revealed activity throughout the motor circuit (Supplement 2). Positive correlations with motor cortex were apparent in primary and secondary motor cortex, SMA, somatosensory cortex, superior temporal gyrus, dorsal anterior cingulate cortex, inferior parietal lobule, cerebellum, orbitofrontal cortex, insula, claustrum, substantia nigra, lenticular nuclei, thalami, and red nuclei. Negative correlations were present in ventral anterior cingulate cortex, superior regions of the parietal lobe, middle and superior frontal gyri, parahippocampal regions, and caudate nuclei. Most regions peaked around the same time as motor cortex with a few exceptions. For example, the loci in orbitofrontal cortex (Brodmann area 47 in the most ventral slice) and the lenticular nuclei appeared to peak earlier.
The maps of cross-correlations to motor cortex during intentional movements were very similar to those seen during the tics (Supplement 3). Group contrast maps were computed at each time point but nothing was significant after Bonferroni correction.
ROI Analyses
The interaction between time point and group was significant in the SMA after a Bonferroni correction for the two primary ROIs examined (p = .003), indicating that this region had a different temporal profile during tics than during intentional movements. The temporal cross-correlations between motor cortex and this area are shown in Figure 1. The temporal cross-correlations for the other regions (with no significant group by time point interactions) are shown in Supplement 4.
Figure 1Temporal cross-correlation between the functionally defined motor cortical reference region and the SMA during tics (solid line) and during intentional movements in control subjects (dashed line). Asterisks indicate differences between patients and control subjects at specific time points as identified in the post hoc t tests. SMA is more active both prior to and after the motor cortical reference region during tics than during intentional movements. Mrr, motor cortical reference region; SMA, supplementary motor area.
A significant group by time point effect (p = .004) was found using the 14 patients known to experience premonitory urges and their matched control subjects.
2
A significant group by time point effect was found (p = .01) using the 11 patients who did not have comorbid conditions and their matched control subjects.
3
A significant group by time point effect was found in the SMA at an uncorrected level that approached significance at a corrected level (p = .028) using the 13 patients with a diagnosis of Tourette syndrome and their matched control subjects.
4
No significant group by time point effect was found when limiting analysis to the 8 patients who were not taking neuropsychiatric medication and their matched control subjects. However, this may be due to the large reduction in power associated with halving our sample size rather than to drug effects. To further investigate this issue, an ANOVA was conducted on the patients' SMA data with medication status as a factor. This analysis revealed no significant time point by medication status interaction, suggesting that the time course of activity in the SMA was similar across patients, regardless of whether they were taking neuropsychiatric drugs or not. Furthermore, patients both on and off drugs showed similar patterns of broadening in the time course of the SMA when compared with their matched control subjects, as illustrated in Supplement 5. Therefore, the differences we have identified in the SMA are unlikely to be due to a medication effect.
Motion Comparisons
There was no significant difference between groups in motion along any of the six dimensions.
Discussion
Cross-correlations with the time course of motor cortex during tics revealed a progression of activity through motor circuitry. The pattern of neural activity during intentional movements in control subjects was very similar to that obtained during tics. This underscores the importance of collecting control data to differentiate brain activity specific to tics from brain activity generally associated with movement.
Region of interest analyses revealed differences between patients and control subjects in the temporal cross-correlation of the motor cortex and the SMA. The SMA region had a broader profile of cross-correlation to motor cortex in the patients than in the control subjects. This suggests that the SMA has abnormally elevated activity in the seconds preceding and following tic execution. These differences are unlikely to be due to head motion, as the groups did not differ significantly in amount of motion. It is possible that greater clustering over time of motor activity in the patients could yield broader patterns of cross-correlation to motor cortex throughout the brain. However, our control subjects were explicitly instructed to produce their movements in a temporally clustered manner. In addition, our data do not suggest that there was global broadening effect in cross-correlation to motor cortex: most areas did not show such broadening. Therefore, the localized pattern of broadening identified in the SMA is unlikely to be due to a global effect such as that arising from different temporal patterns of movement in the two groups.
One possibility is that the difference found in the SMA across groups is related to premonitory urges. The SMA is a promising candidate locus for these urges given that electrical stimulation of the SMA in healthy subjects can induce urges to move that are remarkably similar to the premonitory urges of TS patients (
). If aberrant activity in the SMA does give rise to premonitory urges, there may be a relationship between the intensity of premonitory urges experienced by patients and the activity level in their SMA prior to tics. Unfortunately, we did not collect clinical data relevant to this issue. Further investigation of this issue would be very interesting but should consider the fact that urges increase when subjects inhibit their tics and thus are context dependent as well as subject specific.
In contrast to the SMA, no differences were found between patients and control subjects in the temporal pattern of their S1 activity. However, individual subjects are likely to be activating discrete portions of this large region (depending on which part of the sensory homunculus their tic/movement was associated). Therefore, the spatial averaging inherent in the ROI analyses of this study likely limited its power to identify interesting effects in S1. Thus, we cannot rule out the possibility that this region is implicated in the premonitory urges of TS.
For serious tic disorders that do not respond well to behavioral or pharmacologic intervention, more invasive treatments, such as repetitive transcranial magnetic stimulation (rTMS) and deep brain stimulation (DBS), may be needed. The most effective locus for rTMS is an open question. Trials of rTMS applied to primary motor, premotor, and prefrontal areas have improved tics only modestly or not at all (
). Given that SMA activity appears to be abnormally elevated prior to tics, this region may provide a more effective target. Two recent reports that rTMS of the SMA reduced TS symptoms (
This study adopted a novel method for extracting event-related brain patterns that involved correlating to a component region in the circuit of interest. This method does not require knowledge regarding the time of occurrence of events. In addition, it does not assume that events are instantaneous in nature or that the hemodynamic responses to closely spaced events sum linearly, as standard event-related analyses do. Therefore, the method is well-suited for event-related analyses of temporally extended events that occur in clusters, such as tics.
However, validity of this method depends upon the selection of a seed region that is involved in the events of interest. The seed region used in this study was the part of primary motor cortex controlling tic musculature. There are several reasons to believe that tics arise via neural activity in the voluntary motor pathway and involve activation of motor cortex. First, in contrast to other movement disorders, many patients with TS report the subjective perception that their tics are intentionally produced, often in response to an almost undeniable urge (
), suggesting this region behaves in an aberrant manner in patients with TS. Finally, evidence that motor cortex is involved in the generation of motor tics comes from neuroimaging studies implicating this region in TS (
The functional neuroanatomy of Tourette's syndrome: An FDG-PET study I. Regional changes in cerebral glucose metabolism differentiating patients and controls.
). Thus, the assumption of this study that motor cortex is involved in tics is consistent with the bulk of literature on the neurobiology of TS.
Interpretation of results using this method must always consider the possibility that the events of interest do not sufficiently affect the variance in the reference region to allow just the identification of event-related brain patterns. In our view, it is likely that the underlying functional connectivity to the reference region is also revealed by this technique. We maximized the degree to which our results were influenced by event-related activity patterns by correlating with primary motor cortex, an area that is known to activate robustly during motor acts, during time periods in which such motor acts were occurring regularly. However, our results may reflect differences in the functional connectivity of motor cortex between patients and control subjects, rather than just reflecting different activity patterns during tics and intentional movements. In either case, our findings highlight the importance of the SMA in tic disorders and suggest that it is influencing primary motor cortex differentially during tics as compared with intentional movements.
There are several limitations to this study. First, clinical data was not collected on all relevant patient characteristics. For example, intensity of premonitory urges, severity of illness, and duration of illness are not characterized in our patient sample. Thus, we are unable to investigate the relationship between SMA activity patterns and these variables. Second, this study used a between-group design. Therefore, we cannot rule out the possibility that differences in SMA activity patterns may be due more to differences between groups than to differences between tics and intentional movements. However, it should be noted that a prior study that attempted to use a within-subject design to compare tics and intentional movements found it very difficult to collect control data (
). The manuscript reporting that previous study also raised issues regarding the normalcy of intentional movement in tic patients and recommended that future studies employ a between-group design, as we have done. Finally, many of the patients in our study were using neuropsychiatric drugs. However, patients on and off medication had no significant differences in their SMA time courses and showed similar patterns of broadening in the time course of the SMA when compared with their matched control subjects. Therefore, our findings are unlikely to be due to a medication effect.
This work was supported by the Tourette Syndrome Association (MH) and National Institute of Mental Health K05MH076273 (JFL).
We thank V. Rekkas, E. Hanrahan, R. Fulbright, and R. Hoffman for their time and expertise; K. Martin and H. Sarofin for their technical assistance; and the subjects in the study for their participation.
All authors report no biomedical financial interests or potential conflicts of interest.