Acute TMS/fMRI response explains offline TMS network effects - An interleaved TMS-fMRI study.

BACKGROUND: Transcranial magnetic stimulation (TMS) is an FDA-approved therapeutic option for treatment resistant depression. However, exact mechanisms-of-action are not fully understood and individual responses are variable. Moreover, although previously suggested, the exact network effects underlying TMS' efficacy are poorly understood as of today. Although, it is supposed that DLPFC stimulation indirectly modulates the sgACC, recent evidence is sparse. METHODS: Here, we used concurrent interleaved TMS/fMRI and state-of-the-science purpose-designed MRI head coils to delineate networks and downstream regions activated by DLPFC-TMS. RESULTS: We show that regions of increased acute BOLD signal activation during TMS resemble a resting-state brain network previously shown to be modulated by offline TMS. There was a topographical overlap in wide spread cortical and sub-cortical areas within this specific RSN#17 derived from the 1000 functional connectomes project. CONCLUSION: These data imply a causal relation between DLPFC-TMS and activation of the ACC and a broader network that has been implicated in MDD. In the broader context of our recent work, these data imply a direct relation between initial changes in BOLD activity mediated by connectivity to the DLPFC target site, and later consolidation of connectivity between these regions. These insights advance our understanding of the mechanistic targets of DLPFC-TMS and may provide novel opportunities to characterize and optimize TMS therapy in other neurological and psychiatric disorders.

When differences matter: rTMS/fMRI reveals how differences in dispositional empathy translate to distinct neural underpinnings of self-other distinction in empathy.

Self-other distinction is crucial for empathy, since it prevents the confusion of self-experienced emotions with those of others. We aimed to extend our understanding of the neurocognitive mechanisms of self-other distinction. Thirty-one female participants underwent continuous theta burst transcranial magnetic stimulation (cTBS) targeting the right supramarginal gyrus (rSMG), a sub-region of the temporoparietal junction previously shown to be involved in self-other distinction, and the vertex, a cortical control site. Right after stimulation they completed a visuo-tactile empathy task in an MRI scanner. Self-other distinction was assessed by differences in emotion judgments, and brain activity between conditions differing in the requirement for self-other distinction. Effects of brain stimulation on self-other distinction depended on individual differences in dispositional empathic understanding: cTBS of rSMG, compared to vertex, enhanced self-other distinction in participants with lower dispositional empathic understanding, but diminished it in participants with higher empathic understanding. On the neural level, this inverse relationship between empathic disposition and self-other distinction performance was linked to a reduction of cTBS-induced rSMG activity in persons with lower dispositional empathy, and an increase in those with higher dispositional empathy. These two opposite impacts of cTBS were associated with two anatomically and functionally distinct networks. These findings open up novel perspectives on the causal role of rSMG in self-other distinction and empathy. They also suggest that considering individual differences may yield novel insights into how brain stimulation affects higher-level affect and cognition, and its neural correlates.

Artificial scotoma estimation based on population receptive field mapping.

Population receptive field (pRF) mapping based on functional magnetic resonance imaging (fMRI) is an ideal method for obtaining detailed retinotopic information. One particularly promising application of pRF mapping is the estimation and quantification of visual field effects, for example scotomata in patients suffering from macular dysfunction or degeneration (MD) or hemianopic defects in patients with intracranial dysfunction. However, pRF mapping performance is influenced by a number of factors including spatial and temporal resolution, distribution of dural venous sinuses and patient performance. This study addresses the ability of current pRF methodology to assess the size of simulated scotomata in healthy individuals. The data demonstrate that central scotomata down to a radius of 2.35° (4.7° diameter) visual angle can be reliably estimated in single subjects using high spatial resolution protocols and multi-channel receive array coils.

Default mode network deactivation during emotion processing predicts early antidepressant response.

Several previous functional magnetic resonance imaging (fMRI) studies have demonstrated the predictive value of brain activity during emotion processing for antidepressant response, with a focus on clinical outcome after 6-8 weeks. However, longitudinal studies emphasize the paramount importance of early symptom improvement for the course of disease in major depressive disorder (MDD). We therefore aimed to assess whether neural activity during the emotion discrimination task (EDT) predicts early antidepressant effects, and how these predictive measures relate to more sustained response. Twenty-three MDD patients were investigated once with ultrahigh-field 7T fMRI and the EDT. Following fMRI, patients received Escitalopram in a flexible dose schema and were assessed with the Hamilton Depression Rating Scale (HAMD) before, and after 2 and 4 weeks of treatment. Deactivation of the precuneus and posterior cingulate cortex (PCC) during the EDT predicted change in HAMD scores after 2 weeks of treatment. Baseline EDT activity was not predictive of HAMD change after 4 weeks of treatment. The precuneus and PCC are integral components of the default mode network (DMN). We show that patients who exhibit stronger DMN suppression during emotion processing are more likely to show antidepressant response after 2 weeks. This is, to our knowledge, the first study to show that DMN activity predicts early antidepressant effects. However, DMN deactivation did not predict response at 4 weeks, suggesting that our finding is representative of early, likely treatment-related, yet unspecific symptom improvement. Regardless, early effects may be harnessed for optimization of treatment regimens and patient care.

Increased neural responses to empathy for pain might explain how acute stress increases prosociality.

Recent behavioral investigations suggest that acute stress can increase prosocial behavior. Here, we investigated whether increased empathy represents a potential mechanism for this finding. Using functional magnetic resonance imaging, we assessed the effects of acute stress on neural responses related to automatic and regulatory components of empathy for pain as well as subsequent prosocial behavior. Stress increased activation in brain areas associated with the automatic sharing of others' pain, such as the anterior insula, the anterior midcingulate cortex, and the primary somatosensory cortex. In addition, we found increased prosocial behavior under stress. Furthermore, activation in the anterior midcingulate cortex mediated the effects of stress on prosocial behavior. However, stressed participants also displayed stronger and inappropriate other-related responses in situations which required them to take the perspective of another person, and to regulate their automatic affective responses. Thus, while acute stress may increase prosocial behavior by intensifying the sharing of others' emotions, this comes at the cost of reduced cognitive appraisal abilities. Depending on the contextual constraints, stress may therefore affect empathy in ways that are either beneficial or detrimental.

Imaging the neuroplastic effects of ketamine with VBM and the necessity of placebo control.

In the last years a plethora of studies have investigated morphological changes induced by behavioural or pharmacological interventions using structural T1-weighted MRI and voxel-based morphometry (VBM). Ketamine is thought to exert its antidepressant action by restoring neuroplasticity. In order to test for acute impact of a single ketamine infusion on grey matter volume we performed a placebo-controlled, double-blind investigation in healthy volunteers using VBM. 28 healthy individuals underwent two MRI sessions within a timeframe of 2 weeks, each consisting of two structural T1-weighted MRIs within a single session, one before and one 45min after infusion of S-ketamine (bolus of 0.11mg/kg, followed by an maintenance infusion of 0.12mg/kg) or placebo (0.9% NaCl infusion) using a crossover design. In the repeated-measures ANOVA with time (post-infusion/pre-infusion) and medication (placebo/ketamine) as factors, no significant effect of interaction and no effect of medication was found (FWE-corrected). Importantly, further post-hoc t-tests revealed a strong "decrease" of grey matter both in the placebo and the ketamine condition over time. This effect was evident mainly in frontal and temporal regions bilaterally with t-values ranging from 4.95 to 5.31 (FWE-corrected at p<0.05 voxel level). The vulnerabilities of VBM have been repeatedly demonstrated, with reports of influence of blood flow, tissue water and direct effects of pharmacological compounds on the MRI signal. Here again, we highlight that the relationship between intervention and VBM results is apparently subject to a number of physiological influences, which are partly unknown. Future studies focusing on the effects of ketamine on grey matter should try to integrate known influential factors such as blood flow into analysis. Furthermore, the results of this study highlight the importance of a carefully performed placebo condition in pharmacological fMRI studies.

Eyetracker-based gaze correction for robust mapping of population receptive fields.

Functional MRI enables the acquisition of a retinotopic map that relates regions of the visual field to neural populations in the visual cortex. During such a "population receptive field" (PRF) experiment, stable gaze fixation is of utmost importance in order to correctly link the presented stimulus patterns to stimulated retinal regions and the resulting Blood Oxygen Level Dependent (BOLD) response of the appropriate region within the visual cortex. A method is described that compensates for unstable gaze fixation by recording gaze position via an eyetracker and subsequently modifies the input stimulus underlying the PRF analysis according to the eyetracking measures. Here we show that PRF maps greatly improve when the method is applied to data acquired with either saccadic or smooth eye movements. We conclude that the technique presented herein is useful for studies involving subjects with unstable gaze fixation, particularly elderly patient populations.

Gender transition affects neural correlates of empathy: A resting state functional connectivity study with ultra high-field 7T MR imaging.

Sex-steroid hormones have repeatedly been shown to influence empathy, which is in turn reflected in resting state functional connectivity (rsFC). Cross-sex hormone treatment in transgender individuals provides the opportunity to examine changes to rsFC over gender transition. We aimed to investigate whether sex-steroid hormones influence rsFC patterns related to unique aspects of empathy, namely emotion recognition and description as well as emotional contagion. RsFC data was acquired with 7Tesla magnetic resonance imaging in 24 male-to-female (MtF) and 33 female-to-male (FtM) transgender individuals before treatment, in addition to 33 male- and 44 female controls. Of the transgender participants, 15 MtF and 20 FtM were additionally assessed after 4 weeks and 4 months of treatment. Empathy scores were acquired at the same time-points. MtF differed at baseline from all other groups and assimilated over the course of gender transition in a rsFC network around the supramarginal gyrus, a region central to interpersonal emotion processing. While changes to sex-steroid hormones did not correlate with rsFC in this network, a sex hormone independent association between empathy scores and rsFC was found. Our results underline that 1) MtF transgender persons demonstrate unique rsFC patterns in a network related to empathy and 2) changes within this network over gender transition are likely related to changes in emotion recognition, -description, and -contagion, and are sex-steroid hormone independent.

Model-free fMRI group analysis using FENICA.

Exploratory analysis of functional MRI data allows activation to be detected even if the time course differs from that which is expected. Independent Component Analysis (ICA) has emerged as a powerful approach, but current extensions to the analysis of group studies suffer from a number of drawbacks: they can be computationally demanding, results are dominated by technical and motion artefacts, and some methods require that time courses be the same for all subjects or that templates be defined to identify common components. We have developed a group ICA (gICA) method which is based on single-subject ICA decompositions and the assumption that the spatial distribution of signal changes in components which reflect activation is similar between subjects. This approach, which we have called Fully Exploratory Network Independent Component Analysis (FENICA), identifies group activation in two stages. ICA is performed on the single-subject level, then consistent components are identified via spatial correlation. Group activation maps are generated in a second-level GLM analysis. FENICA is applied to data from three studies employing a wide range of stimulus and presentation designs. These are an event-related motor task, a block-design cognition task and an event-related chemosensory experiment. In all cases, the group maps identified by FENICA as being the most consistent over subjects correspond to task activation. There is good agreement between FENICA results and regions identified in prior GLM-based studies. In the chemosensory task, additional regions are identified by FENICA and temporal concatenation ICA that we show is related to the stimulus, but exhibit a delayed response. FENICA is a fully exploratory method that allows activation to be identified without assumptions about temporal evolution, and isolates activation from other sources of signal fluctuation in fMRI. It has the advantage over other gICA methods that it is computationally undemanding, spotlights components relating to activation rather than artefacts, allows the use of familiar statistical thresholding through deployment of a higher level GLM analysis and can be applied to studies where the paradigm is different for all subjects.

Long-term efficacy and respective potencies of botulinum toxin A and B: a randomized, double-blind study.

BACKGROUND: Mouse units (mU) are used for quantification of the biological activity of botulinum A and B toxin preparations. However, in human tissue, mU values between preparations are not equivalent and lack of clarity concerning efficacy and safety remains with regard to their respective potencies, duration of drug effect and diffusion qualities. OBJECTIVES: To compare short-term and long-term effects of Botox(®) (BOT; Allergan Inc., Irvine, CA, U.S.A.) and Neurobloc(®)/Myobloc(®) (NBC; Solstice Neurosciences Inc., Malvern, PA, U.S.A.) in different doses and dilutions in a human skin model. METHODS: In this prospective randomized double-blind study, 18 healthy volunteers (eight women and 10 men; mean ± SD age 28·4 ± 5·7 years) were injected intradermally with pure saline, BOT and NBC at 10 points in the abdomen in random order, using the BOT/NBC conversion ratio 1 : 75 and different dilution schemes. For an objective outcome, the ninhydrin sweat test was used to compare the anhidrotic areas (action halos). Ten measurements were taken during a time period of 54 weeks. RESULTS: Both preparations showed a peak effect at week 3, with significantly larger anhidrotic areas for NBC. Thereafter, however, the rate of decline was lower in BOT and after week 24, mean BOT areas were larger. The effect of dilution was higher in NBC and the mean dose equivalence conversion ratio (BOT/NBC) was 1 : 29 (area under the curve). Gender effects were seen in both products, with smaller action halos in women. CONCLUSIONS: These results have important implications in clinical routine, especially for autonomic indications.

Fully exploratory network ICA (FENICA) on resting-state fMRI data.

Independent component analysis (ICA) is one of the most valuable explorative methods for analyzing resting-state networks (RSNs) in fMRI, representing a data-driven approach that enables decomposition of high-dimensional data into discrete components. Extensions to a group-level suffer from the drawback of evaluating single-subject resting-state components of interest either using a predefined spatial template or via visual inspection. FENICA introduced in the context of group ICA methods is based solely on spatially consistency across subjects directly reflecting similar networks. Therefore, group data can be processed without further visual inspection of the single-subject components or the definition of a template (Schöpf et al., 2009). In this study FENICA was applied to fMRI resting-state data from 28 healthy subjects resulting in eight group RSNs. These RSNs resemble the spatial patterns of the following previously described networks: (1) visual network, (2) default mode network, (3) sensorimotor network, (4) dorsolateral prefrontal network, (5) temporal prefrontal network, (6) basal ganglia network, (7) auditory processing network, and (8) working memory network. This novel analysis approach for identifying spatially consistent networks across a group of subjects does not require manual or template-based selection of single-subject components and, therefore, offers a truly explorative procedure of assessing RSNs.

Modulation of hypothalamus and amygdalar activation levels with stimulus valence.

In spite of long-standing evidence showing that the hypothalamus is instrumental in generating behaviors associated with positive and negative emotions, little is known about the role of the hypothalamus in normal human emotional processing. Recent findings have suggested that the hypothalamus plays a role beyond mere control of HPA-axis function; this is also supported by the existence of rich anatomical connections between the hypothalamus and the amygdala, a region known for its important role in emotional processing. However, evidence of emotion-induced hypothalamic activity from neuroimaging studies has been inconsistent, possibly due to methodological limitations (e.g., low spatial resolution). Taking advantage of recent improvements in fMRI technology we set out to explore a possible valence-dependent modulation of hypothalamic activity. Using second order parametric analysis of high-resolution BOLD fMRI, we assessed hypothalamic activation patterns during passive viewing of visual stimuli of varying valence, and compared the results with the activity pattern in the amygdalae, i.e. nuclei with known valence-dependent activity profiles. We show that both hypothalamic and amygdalar activation is modulated by the second-order stimulus valence term, i.e., there is increased neural activity following the processing of both positive and negative stimuli. Our results suggest that the hypothalamus may serve a role in generating emotions broader than generally assumed.

[Functional magnetic resonance imaging with ultra-high fields]. / Funktionelle Magnetresonanztomographie bei ultrahohen Feldern.

Functional magnetic resonance imaging (fMRI) is currently the primary method for non-invasive functional localization in the brain. With the emergence of MR systems with field strengths of 4 Tesla and above, neuronal activation may be studied with unprecedented accuracy. In this article we present different approaches to use the improved sensitivity and specificity for expanding current fMRT resolution limits in space and time based on several 7 Tesla studies. In addition to the challenges that arise with ultra-high magnetic fields possible solutions will be discussed.

Optimized 3 T EPI of the amygdalae.

The optimum parameters for single-shot gradient-recalled (GR) EPI-based fMRI studies of the limbic region are systematically established at 3 T via their ability to mitigate intravoxel dephasing-measured via SNR and T2* in the amygdalae-and their implications for temporal resolution (or brain coverage). Conventional imaging parameters (64 x 64 matrix size and 4-6 mm thick slices) are confirmed to be inadequate for functional studies at 3 T. Measurements of main magnetic field variations across the amygdalae suggest that such variations are equal in the craniocaudal and anterior-posterior directions, and slightly lower in the mediolateral direction, with this and other considerations leading us to conclude an oblique axial orientation to be most suitable. In-plane resolution of approximately 1.7 mm was sufficient to recover signal in the area of the amygdalae. SNR was found to peak at a slice thickness of between 2.0 and 2.5 mm, dependent on the subject. T2* time in the amygdalae was measured with a standard EPI protocol to be 22 +/- 3 ms. Using the optimized (high resolution) EPI protocol proposed here, the measured T2* time increased to 48 +/- 2 ms (compared with 43 +/- 3 ms for a reference FLASH scan), only slightly lower than the cortex (49 +/- 2 ms measured with optimized EPI and 52 +/- 2 ms with FLASH). The FLASH measurement of 43 ms is taken to be a suitable effective echo time (TE(eff)) to achieve maximum BOLD sensitivity in the amygdalae. Time series data acquired with these parameters showed a 60% increase in SNR in the amygdala over that obtained with a standard low-resolution protocol and suggest sufficient SNR and BOLD sensitivity to make functional studies feasible. Arteries, but no substantial draining veins, were found in high-resolution BOLD venograms of the region. Our results indicate that EPI protocols need to be carefully optimized for structures of interest if reliable results from single subjects are to be established in this brain region.

The preparation and execution of self-initiated and externally-triggered movement: a study of event-related fMRI.

Studies of functional brain imaging in humans and single cell recordings in monkeys have generally shown preferential involvement of the medially located supplementary motor area (SMA) in self-initiated movement and the lateral premotor cortex in externally cued movement. Studies of event-related cortical potentials recorded during movement preparation, however, generally show increased cortical activity prior to self-initiated movements but little activity at early stages prior to movements that are externally cued at unpredictable times. In this study, the spatial location and relative timing of activation for self-initiated and externally triggered movements were examined using rapid event-related functional MRI. Twelve healthy right-handed subjects were imaged while performing a brief finger sequence movement (three rapid alternating button presses: index-middle-index finger) made either in response to an unpredictably timed auditory cue (between 8 to 24 s after the previous movement) or at self-paced irregular intervals. Both movement conditions involved similar strong activation of medial motor areas including the pre-SMA, SMA proper, and rostral cingulate cortex, as well as activation within contralateral primary motor, superior parietal, and insula cortex. Activation within the basal ganglia was found for self-initiated movements only, while externally triggered movements involved additional bilateral activation of primary auditory cortex. Although the level of SMA and cingulate cortex activation did not differ significantly between movement conditions, the timing of the hemodynamic response within the pre-SMA was significantly earlier for self-initiated compared with externally triggered movements. This clearly reflects involvement of the pre-SMA in early processes associated with the preparation for voluntary movement.

Characterization of BOLD activation in multi-echo fMRI data using fuzzy cluster analysis and a comparison with quantitative modeling.

A combination of multiple gradient-echo imaging and exploratory data analysis (EDA), i.e. fuzzy cluster analysis (FCA), is proposed for separation and characterization of BOLD activation in single-shot spiral functional magnetic resonance imaging (fMRI) experiments at 3 T. Differentiation of functional activation using FCA is performed by clustering pixel signal changes (DeltaS) as a function of echo time (TE). Further vascular classification is supported by the localization of activation and the comparison with a single-exponential decay model. In some subjects, an additional indication for large vessels within a voxel was found as oscillation of the fMRI signal difference vs echo time (TE). Such large vessels may be separated from small vessel activation and, therefore, our proposed procedure might prove useful if a more specific functional localization is desired in fMRI. In addition to the signal change DeltaS, DeltaT(2)*/T(2)* is significantly different between activated regions. Averaged over all eight subjects DeltaT(2)* is 1.7 +/- 0.2 ms in ROIs with the highest signal change characterized as containing large vessels, whereas in ROIs corresponding to microvascular environment average DeltaT(2)* values are 0.8 +/- 0.1 ms.

Dissociation of supplementary motor area and primary motor cortex in human subjects when comparing index and little finger movements with functional magnetic resonance imaging.

This study provides the first investigation of supplementary motor area (SMA) and primary motor cortex (MI) activation with similar movements differing only in subjective difficulty of motor control. Brain activation with simple tapping of the right index finger (well trained during daily life and easy to perform) was compared with tapping of the little finger (less trained and difficult to perform) using functional magnetic resonance imaging at 3 Tesla. Due to optimised movement standardisation, extrinsic influences on activation levels such as movement complexity, amplitude and frequency were minimised. Fifth finger tapping significantly increased the number of activated SMA voxels by 450% whereas MI activation showed no significant difference between fingers. We conclude that with similar movements the degree of subjective difficulty specifically modifies SMA but not MI activation.

Evidence for premotor cortex activity during dynamic visuospatial imagery from single-trial functional magnetic resonance imaging and event-related slow cortical potentials.

A strong correspondence has been repeatedly observed between actually performed and mentally imagined object rotation. This suggests an overlap in the brain regions involved in these processes. Functional neuroimaging studies have consistently revealed parietal and occipital cortex activity during dynamic visuospatial imagery. However, results concerning the involvement of higher-order cortical motor areas have been less consistent. We investigated if and when premotor structures are active during processing of a three-dimensional cube comparison task that requires dynamic visuospatial imagery. In order to achieve a good temporal and spatial resolution, single-trial functional magnetic resonance imaging (fMRI) and scalp-recorded event-related slow cortical potentials (SCPs) were recorded from the same subjects in two separate measurement sessions. In order to reduce inter-subject variability in brain activity due to individual differences, only male subjects (n = 13) with high task-specific ability were investigated. Functional MRI revealed consistent bilateral activity in the occipital (Brodmann area BA18/19) and parietal cortex (BA7), in lateral and medial premotor areas (BA6), the dorsolateral prefrontal cortex (BA9), and the anterior insular cortex. The time-course of SCPs indicated that task-related activity in these areas commenced approximately 550-650 ms after stimulus presentation and persisted until task completion. These results provide strong and consistent evidence that the human premotor cortex is involved in dynamic visuospatial imagery.

Spatial distribution of prostate cancers undetected on initial needle biopsies.

OBJECTIVES: The spatial distribution of cancer foci of prostate carcinomas with negative initial biopsies was compared to that of prostate carcinomas with positive initial biopsies to detect areas in which carcinomas were more frequently located when the initial biopsy was negative. METHODS: Twenty patients with prostate cancer and a negative initial biopsy trial were detected among 218 patients with preceding systematic biopsies (9.2%) in our hospital. Analysis of the prostatectomy specimens regarding cancer distribution, multifocality, tumour size, Gleason score, and stage was performed using pathohistological techniques and three-dimensional computer reconstruction. RESULTS: Prostatectomy specimens with negative initial biopsies showed more frequently cancer foci in apical (p<0.0001) and dorsal (p<0.02) prostatic compartments, higher incidence of multifocality (p<0.01), and smaller size of carcinoma foci (p<0.00003) compared to carcinomas in 81 stage-matched prostatectomy specimens with positive initial biopsies. Comparing both groups, no significant differences were noted in Gleason score of preoperative biopsies and prostatectomies, prostate weight, prostate-specific antigen (PSA) level, digital rectal examination, and patients age. CONCLUSIONS: Missing the cancer in clinically significant prostate carcinomas by current systematic biopsy techniques may also be due to an apico-dorsal cancer location, particularly in combination with multifocality and small size of carcinoma foci in large prostates. In case of reasonable clinical suspicion of prostate cancer and negative initial biopsy, an early repeat biopsy with special emphasis on the apico-dorsal peripheral zone should be envisaged.

Finger somatotopy in human motor cortex.

Although qualitative reports about somatotopic representation of fingers in the human motor cortex exist, up to now no study could provide clear statistical evidence. The goal of the present study was to reinvestigate finger motor somatotopy by means of a thorough investigation of standardized movements of the index and little finger of the right hand. Using high resolution fMRI at 3 Tesla, blood oxygenation level-dependent (BOLD) responses in a group of 26 subjects were repeatedly measured to achieve reliable statistical results. The center of mass of all activated voxels within the primary motor cortex was calculated for each finger and each run. Results of all runs were averaged to yield an individual index and little finger representation for each subject. The mean center of mass localizations for all subjects were then submitted to a paired t test. Results show a highly significant though small scale somatotopy of fingerspecific activation patterns in the order indicated by Penfields motor homunculus. In addition, considerable overlap of finger specific BOLD responses was found. Comparing various methods of analysis, the mean center of mass distance for the two fingers was 2--3 mm with overlapping voxels included and 4--5 mm with overlapping voxels excluded. Our data may be best understood in the context of the work of Schieber (1999) who recently described overlapping somatotopic gradients in lesion studies with humans.