Previous statistical voxelwise lesion-behavior mapping (VLBM) studies have demonstrated that spatial neglect is associated with cortical and subcortical gray matter damage. However, it has also been suggested that the disorder may result from white matter injury. Our aim was to investigate the white matter connectivity in a large sample of 140 stroke patients. We combined a VLBM approach with the histological maps of the human white matter fiber tracts provided by the Jülich probabilistic cytoarchitectonic atlas. We found that damage of right perisylvian white matter connections—the superior longitudinal fasciculus, the inferior occipitofrontal fasciculus, and the superior occipitofrontal fasciculus—is a typical finding in patients with spatial neglect. However, the analysis also revealed that the largest portion of the lesion area, namely between 89.1% and 96.6%, affected brain structures other than the perisylvian white matter fiber tracts. Predominantly, these included gray matter structures such as the superior temporal, inferior parietal, inferior frontal, and insular cortices, as well as subcortically the putamen and the caudate nucleus. Damage of gray matter structures thus appears to be a strong predictor of spatial neglect.
Spatial neglect is a common consequence of stroke and is typically observed in individuals who suffer from right hemisphere strokes (Pedersen et al. 1997; Hillis et al. 2005; Becker and Karnath 2007). Patients with spatial neglect spontaneously deviate their eyes and head toward the ipsilesional, right side (Fruhmann Berger et al. 2006) and tend to ignore information presented on their contralesional side. Clinically, spatial neglect is a strong predictor for poor recovery on a wide range of everyday tasks (Jehkonen et al. 2000). The syndrome also provides theoretical insight into the normal function of the spatial attentional system in humans.
Structural imaging can be combined with new statistical methods to correlate the location of brain injury with behavioral deficits (for review, Rorden and Karnath 2004). These techniques can offer an objective measure for lesion-behavior relationships. Indeed, a number of rigorous statistical methods have been developed for conducting voxelwise lesion-behavior mapping (VLBM), where one does not make any explicit, a priori assumptions regarding brain anatomy (Frank et al. 1997; Bates et al. 2003; Kimberg et al. 2007; Rorden et al. 2007).
Recently, a series of studies have examined the anatomical correlates of spatial neglect by using these statistical VLBM methods (Mort et al. 2003; Karnath et al. 2004; Committeri et al. 2007). The results from these studies were interpreted as implicating cortical and subcortical gray matter structures. However, none of these studies explicitly attempted to assess involvement of white matter connections. This raises the possibility that these studies might have a bias toward identifying gray matter structures rather than nearby white matter fiber bundles. It was even speculated that “…a particular form of disconnection might have greater predictive value than the localization of gray matter lesions concerning the [neglect] patients’ deficits and disabilities” (Bartolomeo et al. 2007, p. 2484).
Our aim was to address the issue of involvement of gray versus white matter tissue in spatial neglect. In particular, we were interested in investigating the relationship of white matter injury to the presence of the disorder in a large sample of 140 stroke patients. In theory, this work should help identify if white matter injury is a common predictor of spatial neglect and, if so, which of the fiber bundles typically are involved.
The first step of this approach was a reanalysis of a large imaging dataset of right hemisphere stroke lesions described by Karnath et al. (2004) consisting of a 7-year sample of 140 patients by using a new statistical VLBM protocol. We have recently advocated the use of the Liebermeister measure for binomial data in VLBM studies (Rorden et al. 2007). As this test does not assume fixed marginals, it is more sensitive than the previously used chi-square test (which is overly conservative). The second step of our present approach of white matter fiber tract analysis is the combination of the statistical lesion map resulting from the VLBM procedure with a recently developed set of cytoarchitectonic probabilistic maps of white matter fiber connections (Amunts and Zilles 2001; Zilles et al. 2002; Bürgel et al. 2006). This atlas provides stereotaxic information on the location and variability of cortical areas and white matter fiber tracts in the Montreal Neurological Institute (MNI) reference space (http://www.fz-juelich.de/inb/inb-3/index.php?index=29). In contrast to the reference brain of the Talairach and Tournoux atlas (Talairach and Tournoux 1988) or the MNI single subject or group templates (Evans et al. 1992; Collins et al. 1994; Tzourio-Mazoyer et al. 2002), the probabilistic cytoarchitectonic atlas is based on the analysis of the cytoarchitecture in a sample of 10 different postmortem human brains. By using a modified myelin staining technique, Bürgel et al. (2006) were able to distinguish and map individual tracts at microscopic resolution within the densely packed fibers of the white matter in these 10 human brains. Digitalization and normalization of results derived 3-dimensional, probabilistic maps of 10 different long white matter fiber bundles. The maps illustrate the relative frequency with which a certain fiber tract of the 10 normal human brains was present on an MNI reference brain in a voxel. They thus serve as a measure of intersubject variability of the human white matter fiber bundle anatomy for each voxel of the reference space.
Subjects and Methods
Data were drawn from 140 consecutively admitted stroke patients with circumscribed right hemisphere lesions from a well-defined recruitment area belonging to the University of Tübingen, sampled over a period of 7 years (Karnath et al. 2004). The full details regarding patient characteristics, test criteria, etc., are reported in our previous work (Karnath et al. 2004). In short, following standardized testing for spatial neglect by using the “Letter cancellation” task (Weintraub and Mesulam 1985), the “Bells test” (Gauthier et al. 1989), the “Baking tray task” (Tham and Tegnér 1996), and a copying task (Johannsen and Karnath 2004), the subjects were divided into a group of 78 patients with spatial neglect and 62 control patients who did not show the disorder. Magnetic resonance imaging (MRI) or computerized tomography (Spiral-CT) was carried out in each subject. Lesion mapping was aided by diffusion-weighted imaging for MRI occurring within the first 48 h poststroke and T2-weighted fluid-attenuated inversion recovery sequences when MRI was conducted 48 h or later after the stroke. The mean time between lesion and the MRI was 5.0 days (standard deviation [SD] 5.4). In those subjects who underwent the CT imaging protocol, the mean time since lesion and the CT was 6.7 days (SD 8.4). All lesions were drawn manually (while being blind for the diagnosis of spatial neglect) on axial slices of a T1-weighted template MRI scan from the MNI (http://www.bic.mni.mcgill.ca/cgi/icbm_view) using the MRIcro software package (Rorden and Brett 2000) with a 1 × 1 mm in-plane resolution. This template is approximately oriented to match Talairach space (Talairach and Tournoux 1988) and is distributed with MRIcro. Lesions were mapped onto the slices that correspond to z-coordinates −40, −32, −24, −16, −8, 0, 8, 16, 24, 32, 40, and 50 mm in Talairach space by using the identical or the closest matching axial slices of each individual. The full details regarding imaging protocols etc. are reported in our previous work (Karnath et al. 2004).
Analysis of White Matter Lesion Anatomy
For the present investigation, we used the new nonparametric mapping software that is distributed as part of the MRIcron toolset (Rorden et al. 2007). Specifically, we performed the nonparametric Liebermeister test to identify voxels that when injured predicted the presence of spatial neglect. The resulting statistical map was adjusted for multiple comparisons by using a Bonferroni-corrected P < 0.05 threshold.
Therefore, our initial analysis was based on this statistical map that illustrated all the voxels that predicted the presence of spatial neglect. In a first analysis, we determined the percentage of these voxels that comprised white matter compared with gray matter tissue. We segmented the individual T1-weighted template MRI scan from the MNI (http://www.bic.mni.mcgill.ca/cgi/icbm_view) by using Statistical parametric mapping routines (SPM5; http://www.fil.ion.ucl.ac.uk/spm/). The percentage of overlap of white versus gray matter tissue with the statistical lesion map was determined by using MRIcron software (Rorden et al. 2007).
In order to identify the specific fiber tracts that predict spatial neglect, we segmented the thresholded Liebermeister statistical map using the white matter fiber tracts from the human probabilistic cytoarchitectonic atlas (Bürgel et al. 2006). This atlas is in the same space as the MNI reference brain, with each atlas map illustrating the relative frequency with which a certain fiber tract of 10 normal human brains was present (e.g., a 30% value reflects that the fiber tract was present in that voxel for 3 out of 10 brains). The number of overlapping voxels between the statistical lesion map and the anatomical fiber tract map in the right hemisphere was determined by using the SPM anatomy toolbox of the Jülich atlas (http://www.fz-juelich.de/inb/inb-3//spm_anatomy_toolbox; Eickhoff et al. 2005).
We also conducted a logistic regression that included each white matter fiber tract from the Jülich probabilistic atlas (Bürgel et al. 2006). The extent of injury to each fiber tract was calculated for each individual on a voxel-by-voxel basis, weighting each lesioned voxel by the frequency of the target fiber tract in the atlas. These continuous variables were used to predict the absence or presence of spatial neglect, adding overall lesion volume as a predictor. Note that logistic regression computes the partial correlation for each factor and adding an extra variable (such as overall lesion volume) can increase (if the variable predominantly describes variance that is not accounted for by other variables) or decrease (if the variable predominantly describes variance that is redundant with another variable) statistical significance of the other variables. Finally, we conducted a restricted version of this analysis, where we excluded primary sensory and primary motor tracts (acoustic radiation [AR], corticospinal tract [CT], optic radiation [OR]). Our reasoning was that these fiber bundles are associated with lower level perception and motor control, whereas neglect behavior is often seen as distinctively different from injury that is restricted to such primary processing. On explicit testing, the symptoms of spatial neglect appear to clearly dissociate from the symptoms associated with injury to these primary fiber tracts. Specifically, many neglect patients are not physically blind or deaf, rather they fail to attend and orient toward left-sided visual or auditory information.
Figure 1 shows all voxels that were—after adjustment for multiple comparisons—significantly more frequently affected in those patients showing spatial neglect than in those without the disorder by using the new Liebermeister approach. In comparison with our previous results obtained with chi-square statistics (Karnath et al. 2004), the present study revealed a larger territory of significant difference. This reflects the overly conservative nature of the original chi-square test, which resulted in fewer voxels surviving the statistical threshold.
For a first analysis of the white matter lesion anatomy, we determined the percentage to which the statistical lesion map (Fig. 1) affected white compared with gray matter brain tissue. We found that 36.9% of the map overlapped with the white matter portion of the segmented MNI brain, whereas 63.1% of the map overlapped with the gray matter segment. White matter brain tissue thus was less prominently affected than gray matter tissue in those patients showing spatial neglect than in those without the disorder. By using the anatomical parcellation of the MNI single subject brain by Tzourio-Mazoyer et al. (2002) implemented in MRIcron, we found the most statistically significant voxels at MNI coordinates (x 66; y −16; z 8) within the right superior temporal gyrus (STG) and coordinates (x 28; y −7; z 24) within the right periventricular white matter. Regions significantly damaged also encompassed large parts of the right insula, the planum temporale, and operculum (including the opercular portion of the inferior frontal gyrus). The lesion area further extended into the pre- and postcentral gyri, the inferior parietal lobule, and subcortically into the putamen, head of caudate nucleus, as well as external and extreme capsule.
For a closer analysis of the white matter lesion anatomy, we then determined which fiber tracts were affected within the white matter portion of the lesion territory. Figure 2a illustrates the statistical lesion map obtained from the VLBM analysis (Fig. 1) in homogenous brown color. It was superimposed onto the anatomical maps of the white matter fiber tracts from the probabilistic Jülich atlas. Figure 3 gives the quantitative analysis of these data. When we included every single voxel of the anatomical atlas that is, even those voxels in which fiber bundle tissue had been found in only 1 out of the 10 postmortem brains, we revealed that 22.4% of the superior longitudinal fasciculus (SLF), 26.3% of the inferior occipitofrontal fasciculus (IOF), and 39.6% of the superior occipitofrontal fasciculus (SOF) as well as 7.5% of the uncinate fasciculus (UF), 18.4% of the CT, and 46.9% of the AR were located within the territory that was more frequently affected in the group with spatial neglect than in controls (Figs. 2 and 3). There was no overlap between this area and the optic radiation, corpus callosum, fornix, and cingulum from the Jülich atlas. Following the threshold procedure suggested for defining anatomical regions of interest on the basis of the Jülich atlas for functional imaging data (Eickhoff et al. 2005, 2006), we then focused on voxels located within the 30% isocontours of the atlas, that is, those voxels in which a fiber tract was present in at least 3 (30%) of the 10 postmortem brains, to increase the reliability of fiber tract anatomy in the reference map. Under this condition, we revealed that 7.0% of the SLF, 8.2% of the IOF, and 12.7% of the SOF as well as 0.6% of the UF, 9.5% of the CT, and 15.0% of the AR were in the territory identified by the statistically thresholded map.
The results from the logistic regression are shown in Table 1. When overall lesion volume was taken into account, the CT was the only white matter fiber tract that predicted spatial neglect, surviving adjustment for multiple comparisons by a Bonferroni-corrected P < 0.05 threshold. This model correctly categorized 55 of the 62 patients without neglect and 64 of the 78 patients exhibiting spatial neglect. A correlation matrix for all factors is presented in Table 2. It revealed that lesion volume (r = 0.548) and the CT (r = 0.517) were the strongest individual predictors of spatial neglect. We strongly suspect that the disconnection of the CT has no direct relevance for the presence of spatial neglect (but rather for the frequent additional presence of right-sided hemiparesis in these patients), and we therefore theorize that the association of the CT to spatial neglect observed with the logistic regression reflects its proximity to neighboring cortical and subcortical areas. Note also that the numerical trend was for injury to the cingulum and fornix to be negatively correlated with neglect, suggesting that these regions are distant from the core areas responsible for spatial neglect.
|Region||P value||P value (partial model)|
|Region||P value||P value (partial model)|
Note: A logistic regression was conducted with the absence or presence of spatial neglect as the dependent variable and the extent of injury to each of the white matter fiber tracts from the Jülich atlas as the independent variables. Overall lesion volume was included as a predictor. Logistic regression assesses partial correlation—measuring whether adding the variable improves the model when all the other variables are already known. Asterisks (*) highlight regions that survived a 5% Bonferroni correction (P < 0.0045 for the 11-factor comparison). We also conducted a partial model where the AR, Corpus Callosum, and OR were excluded (Bonferroni threshold P < 0.00625 for the 8-factor comparison). The results from this partial model are shown in the right column.
Note: Also included is presence of spatial neglect (Neglect) and overall lesion volume (Volume). CC, corpus callosum.
We also conducted a logistic regression that excluded the primary sensory and primary motor fiber bundles (AR, CT, OR). This model included 8 variables (7 fiber bundles and overall lesion volume). The results from this analysis are shown in Table 1. Two variables (IOF and lesion volume) exceeded the Bonferroni threshold of 0.00625 (P < 0.05 corrected for multiple comparisons), with the UF approaching significance. The reason these fiber tracts are not detected in the full model is likely due to their strong correlation with the CT (as shown in Table 2). Despite the reduced number of factors, this model correctly categorized 53 of the 62 patients without neglect and 63 of the 78 patients exhibiting spatial neglect.
Finally, we determined the overall proportion of the lesion territory that overlapped with the white matter fiber tracts identified by our analyses. The CT and the AR were not included due to their obvious involvement in primary motor and acoustic processing. Rather, in this calculation, we focused on those tracts that might be of relevance for spatial neglect in case of its disconnection, that is, the SLF, IOF, SOF, or UF. With the inclusion of every single voxel of the anatomical atlas, that is, even those voxels in which fiber bundle tissue had been found in only 1 out of the 10 postmortem brains, we revealed 10.9% of overlap. In reference to the 30% isocontours of the atlas, we found an overlap of 3.4%. This indicates that between 89.1% and 96.6% of the area, which was significantly more damaged in patients with spatial neglect than in controls, affected neural structures other than the long white matter association fibers SLF, IOF, SOF, and UF, namely predominantly the cortical and subcortical gray matter structures listed above (Fig. 1) as well as the primary motor and primary sensory fiber bundels CT and AR in the white matter.
In order to assess the involvement of gray versus white matter tissue in spatial neglect, we conducted a reanalysis of the 140 patients sample described by Karnath et al. (2004). We found that only 36.9% of the territory significantly more affected in patients with spatial neglect than in controls overlapped with white matter tissue. In contrast, 63.1% of this territory affected cortical and subcortical gray matter brain structures. By using the 30% isocontours from the probabilistic cytoarchitectonic atlas provided by Bürgel et al. (2006), that is, the regions where in at least 3 brains from the 10 postmortem brains fiber tract tissue was found histologically, 7.0% of the SLF, 8.2% of the IOF, 12.7% of the SOF, and only 0.6% of the UF overlapped with the territory significantly more affected in patients with spatial neglect than in controls.
The SLF, SOF, and IOF belong to a densely connected perisylvian network in the human right hemisphere composed of the inferior parietal, dorsolateral frontal, superior temporal, and insular cortices (Schmahmann and Pandya 2006; Catani et al. 2007; Schmahmann et al. 2007; Gharabaghi A, Kunath F, Erb M, Grodd W, Tatagiba M, Karnath H-O, unpublished data). The SLF is the major cortical association fiber pathway linking parietal and frontal cortices in both humans and monkeys. It is subdivided into different components (Petrides and Pandya 1984; Makris et al. 2005; Schmahmann et al. 2007). Subcomponent II (SLF II) links the inferior parietal lobule and intraparietal sulcus with the posterior and caudal prefrontal cortex, whereas SLF III connects the rostral inferior parietal lobule with the ventral part of premotor and prefrontal cortex. A further fiber tract stems from the caudal part of the STG, arches around the caudal end of the Sylvian fissure, and extends to the lateral prefrontal cortex along with the SLF II fibers. This latter fiber tract is termed the arcuate fasciculus (AF) (Burdach 1819–1826; Dejerine and Dejerine-Klumpke 1895; Catani et al. 2002, 2005; Makris et al. 2005). Some authors, including those providing the probabilistic white matter fiber tract maps of the Jülich atlas (Bürgel et al. 2006), regarded the AF as a fourth subcomponent of the SLF. Fiber tract “SLF” of the Jülich atlas histologically encompasses subcomponents SLF I–III and the AF (K. Amunts, personal communication).
A fiber bundle situated in close proximity to the SLF II and the AF is the SOF. The SOF was the association fiber tract that was found to be affected most prominently in the present analysis. This long association bundle bidirectionally extends from the inferior parietal lobule and dorsomedial parastriate occipital cortex to caudal, dorsal, and medial frontal lobe areas (Catani et al. 2002; Bürgel et al. 2006; Makris et al. 2007; Schmahmann et al. 2007). Some of its fibers intermingle with SLF II and AF fibers. The IOF runs through the isthmus of the temporal lobe medial to the lower insula and connects ventrolateral frontal cortex with posterior temporal, inferior parietal, and occipital cortices (Nieuwenhuys et al. 1988; Catani et al. 2002; Bürgel et al. 2006).
By using intraoperative electrical stimulation and diffusion tensor imaging (DTI) tractography, Bartolomeo et al. (2007) suggested—in single cases—disconnection of the SLF (Thiebaut de Schotten et al. 2005) and of the IOF (Urbanski et al. 2008) to be associated with spatial neglect. Further, a DTI tractography analysis by He et al. (2007) revealed damage to the SLF and the AF in 5 patients with severe spatial neglect but not in 5 patients with mild neglect. In contrast to these studies, the present investigation used a different technique to investigate the question of whether or not white matter association fibers are predictive of acute spatial neglect. We used a probabilistic cytoarchitectonic atlas based on the histological findings in 10 adult postmortem brains. Further and most importantly, the present study used a large, representative sample of 140 right hemisphere stroke patients.
The majority of voxels where injury predicted spatial neglect compromised gray matter tissue (63.1%), with white matter only accounting for 36.9% of these voxels. However, we found that damage of right perisylvian white matter connections indeed is a typical finding in patients with spatial neglect. The SOF was more prominently affected than the SLF or IOF. However, we did not find a strong predictive relationship between white matter injury and spatial neglect. Overall, the portion of involvement of these fiber bundles was not very high.
Logistic regression analysis identified damage of the CT as the only significant white matter predictor of spatial neglect. Because disconnection of the CT has no known direct relevance for the presence of spatial neglect (but rather for the frequent additional presence of right-sided hemiparesis in these patients), the association of the CT to spatial neglect presumably reflects its proximity to neighboring cortical and subcortical areas. Logistic regression only revealed the IOF as a significant white matter predictor of spatial neglect (with the UF approaching significance) when the primary sensory and primary motor fiber tracts (AR, CT, OR) were excluded from the model. Finally, the percentage of overlap between the statistical lesion map and all the perisylvian network connections (histologically defined via the Jülich atlas) turned out to be small between 3.4% and 10.9%.
It is possible that the histological Jülich atlas (Bürgel et al. 2006) may slightly underestimate the true overlap between the statistical lesion map and the perisylvian white matter network connections. The authors of the atlas cautiously mentioned the possibility that because white matter tracts were mapped on the basis of sections in the coronal plane, fibers that proceed rostrocaudally may have been underestimated. On the other hand, the microscopic mapping used for constructing the Jülich atlas results in a more precise delineation of fiber tracts and enables mapping in some portions of fiber tracts that are beyond the resolution of alternative mapping tools, for example, DTI tractography (Bürgel et al. 2006).
Our findings are in contrast with those of Doricchi and Tomaiuolo (2003), who report that white matter injury strongly predicts spatial neglect. There are clear methodological differences as their population was less acute (mean poststroke interval of 121 days). Therefore, it is possible that white matter involvement may be a marker of those patients who have persistent symptoms. This could be directly studied by conducting an investigation with a large group of patients with chronic injury. The relatively small sample size in the study of Doricchi and Tomaiuolo (2003) makes statistical inference difficult. Although they report Fisher exacts tests that implicated the corticospinal tract, it is unclear how this test was conducted. Specifically, there is no suggestion of controlling for multiple comparisons, the precise region compared (e.g., how was this tract identified) or whether this region was selected post hoc based on the maximum overlap. A further striking difference between this finding and the present observations as well as other work by our group (Karnath et al. 2002) is the high incidence of damage to the putamen in Doricchi and Tomaiuolo's chronic control group who were not experiencing spatial neglect.
The present study cannot decide whether or not spatial neglect should best be interpreted as a “disconnection syndrome” (Watson et al. 1974, 1978, Mesulam and Geschwind 1978; Bartolomeo et al. 2007). Our analyses only demonstrate that the involvement of white matter fiber tracts in the typical lesion territory of patients with spatial neglect is rather small. The overlap between the statistical lesion map and all the perisylvian white matter network connections was between 3.4% and 10.9%. This indicates that between 89.1% and 96.6% of the lesion area affected brain structures other than the perisylvian white matter fiber tracts, namely cortical and subcortical gray matter structures and the primary motor and primary sensory fiber bundles CT and AR. The gray matter structures involved were the superior temporal, inferior parietal, inferior frontal, and insular cortices as well as subcortically the putamen and the caudate nucleus. Hillis et al. (2002) have observed that white matter lesions cause spatial neglect only if the subcortical damage provokes additional malperfusion of cortical gray matter structures in the ipsilesional, right hemisphere. In combination with the present findings, we conclude that damage of gray matter structures appears to be a strong predictor of spatial neglect.
Bundesministerium für Bildung und Forschung (BMBF-Verbundprojekt “Räumliche Orientierung” 01GW0641), the Deutsche Forschungsgemeinschaft (SFB 550-A4); National Institutes of Health (R01 NS054266); and European Union (PERACT—Marie Curie Early Stage Training MEST-CT-2004-504321).
Conflict of Interest: None declared.