Abstract

Corticocortical interactions between the primary motor cortex, the ventral premotor cortex and posterior parietal motor areas, such as the anterior and caudal intraparietal sulcus, are relevant for skilled voluntary hand function. It remains unclear to what extent these brain regions and their interactions also contribute to basic motor functions after stroke. We hypothesized that white matter integrity of the underlying parietofrontal motor pathways between these brain regions might relate to residual motor function after stroke. Twenty-five chronic stroke patients were recruited (aged 64 ± 8.8 years, range 46–75, 17 males, one left-handed) and evaluated 34 months after stroke (range 12–169 months) by means of grip force, pinch force and the Fugl-Meyer assessment of the upper extremity. Based on these measures, motor function was estimated applying a factor analysis with principal component extraction. Using diffusion tensor imaging and probabilistic tractography we reconstructed probable intrahemispheric trajectories between the primary motor cortex, the ventral premotor cortex and the anterior and caudal intraparietal sulcus in each patient. White matter integrity was estimated for each individual tract by means of fractional anisotropy. Generalized linear modelling was used to relate tract-related fractional anisotropy to the motor function. We found that the white matter integrity of the fibre tracts connecting the ventral premotor cortex and the primary motor cortex (P < 0.001) and the anterior intraparietal sulcus and the ventral premotor cortex (P < 0.01) positively correlated with motor function. The other tracts investigated did not show a similar structure-behaviour association. Providing first structural connectivity data for parietofrontal connections in chronic stroke patients, the present results indicate that both the ventral premotor cortex and the posterior parietal cortex might play a relevant role in generating basic residual motor output after stroke.

Interactions between posterior parietal motor areas, ventral premotor cortex and primary motor cortex support skilled hand movements. Using tract-related structure-behaviour analysis, Schulz et al. reveal that the integrity of connections between primary motor cortex and both ventral premotor cortex and the anterior intraparietal sulcus affects residual motor output after stroke.

Interactions between posterior parietal motor areas, ventral premotor cortex and primary motor cortex support skilled hand movements. Using tract-related structure-behaviour analysis, Schulz et al. reveal that the integrity of connections between primary motor cortex and both ventral premotor cortex and the anterior intraparietal sulcus affects residual motor output after stroke.

Introduction

A stroke causing a motor deficit leads to widespread time-dependent alterations of motor-related neural activity (Rehme et al., 2012) which is paralleled by changes of functional (Wang et al., 2010) and effective interregional connectivity patterns (Grefkes et al., 2008; Rehme et al., 2011). Numerous studies have contributed to the understanding of the primary motor cortex (M1), the premotor cortices and the supplementary motor area as important nodes for motor recovery and functioning after stroke (Rehme and Grefkes, 2013). For instance, on the ipsilesional hemisphere, a positive prefrontal–premotor coupling was found to increase after stroke (Sharma et al., 2009) while premotor–premotor (Sharma et al., 2009) and premotor–motor couplings were reduced (Grefkes et al., 2008, 2010; Rehme et al., 2011). Such alterations in brain connectivity have been found functionally relevant as the reinstatement of normal coupling patterns within the ipsilesional hemisphere over time has been related to the amount of recovery (Rehme et al., 2011).

Compared to other secondary motor areas, less is known about the functional relevance of the ventral premotor cortex and the posterior parietal cortex after stroke. Studies in healthy participants have confirmed that the posterior parietal cortex and its pathways linking different secondary motor areas along the intraparietal sulcus (IPS) with ventral premotor cortex and ipsilateral M1 (Rizzolatti et al., 1998; Grefkes and Fink, 2005) mediate skilled voluntary movements, such as reaching and grasping and using objects and tools (Vingerhoets, 2014). These features account for dexterous hand function which is often critically affected after stroke. As an important node in a dorsolateral parietofrontal circuit, the anterior portion of the lateral bank of the IPS has been found to contribute critically to grasping (Davare et al., 2010, 2011), precision grip (Ehrsson et al., 2001) and online-adjustment of grip force (Dafotakis et al., 2008) during object manipulation and exploration (Binkofski et al., 1999). Anterior IPS provides ventral premotor cortex with visual information about object properties and influences the interactions between ventral premotor cortex and M1 to generate task-specific motor commands. Also more posterior regions along the intraparietal sulcus, homologous to the caudal IPS in monkeys, showed grasp-related activation (Faillenot et al., 1997; Culham et al., 2003). Animal (Tsutsui et al., 2003) and human studies (Shikata et al., 2001; Filimon et al., 2007; Koch et al., 2010; Mruczek et al., 2013) have documented that caudal IPS might contribute to the integration of spatial object features relevant for further movement planning and execution.

Structural connections between M1 and ventral premotor cortex as well as ventral premotor cortex and both anterior and caudal IPS have been found in animals (Cavada and Goldman-Rakic, 1989; Luppino et al., 1999; Rozzi et al., 2006; Schmahmann et al., 2007; Borra et al., 2008; Gharbawie et al., 2011) and humans (Makris et al., 2005; Tomassini et al., 2007; Koch et al., 2010). Sparse pathways linking also the anterior IPS-associated supramarginal gyrus and caudal IPS-related angularis gyrus with M1 have been recently described in humans (Koch et al., 2010). The structural substrate of the parietofrontal connections is the superior longitudinal fascicle (SLF). Its major component SLF II predominantly connects caudal inferior-parietal brain regions corresponding to the angularis gyrus with dorsolateral frontal cortices. Its ventral component (SLF III) mainly derives from the more anterior supramarginal gyrus and targets predominantly ventral premotor and prefrontal cortices while a dorsal component (SLF I) connects mediodorsal parts of the parietal cortex with the mediodorsal frontal lobe (Makris et al., 2005).

There is only limited knowledge about the relevance of the parietofrontal motor pathways after stroke. A recent meta-analysis has described a consistent activation of anterior IPS during various movements of the affected limb (Rehme et al., 2012). Resting-state functional MRI has shown an increased connectivity between the contralesional superior parietal lobule and the ipsilesional dorsal premotor cortex that positively related to motor function after subcortical stroke (Wang et al., 2010). Effective connectivity analyses have revealed a decreased coupling of the ipsilesional superior parietal cortex on M1 and the supplementary motor area (Inman et al., 2012).

Given the importance of the parietofrontal interactions for skilled hand function in healthy participants and growing evidence for stroke-related alterations of neural activity in the posterior parietal cortex, we questioned whether the microstructural integrity of specific parietofrontal white matter pathways might contribute also to basic residual motor function after stroke. Indeed, previous imaging studies have clearly shown that the structural integrity of corticocortical pathways is an important basis for neuronal information throughput and relevant for behaviour (Schulz et al., 2014). Using diffusion-tensor imaging we aimed at (i) reconstructing probable intrahemispheric connections between M1, ventral premotor cortex, anterior IPS and caudal IPS; and (ii) examining the extent to which tract-related microstructural integrity correlates with preserved motor function in the chronic stage of recovery after stroke.

Materials and methods

Subjects

Twenty-five patients (aged 64 ± 8.8 years, range 46–75, 17 male, one left-handed) with a broad spectrum of first-ever cortical and subcortical ischaemic strokes (10 in the dominant hemisphere, see Fig. 1 for lesion location and Table 1 for clinical data) were recruited. Patients were evaluated 34 months after stroke (range 12–169 months) by means of grip force, pinch force, the Fugl-Meyer assessment of the upper extremity (UEFM) and the modified Rankin Scale. Based on grip and pinch force and the UEFM score, one composite motor function score was calculated for the affected hand using a factor analysis with principal component analysis. Explaining 73.4% of the variance, the first eigenvariate was used for further correlative analysis. To account for general inter-subject variability in motor performance, a motor function score was also calculated for the unaffected hand based on grip and pinch force variability. A group of 20 healthy age- (67 ± 2.4 years) and gender-matched (10 male) participants was also recruited. Participants gave written informed consent according to the Declaration of Helsinki. The study was approved by the local ethics committee.

Stroke lesions. All masks of stroke lesions were brought to the right hemisphere (affected hemisphere, AH) and overlaid on T1 template in MNI standard space. Axial slices were chosen with given z-values. The colour bar indicates the number of subjects in which voxels lay within a stroke lesion. UH = unaffected hemisphere.
Figure 1

Stroke lesions. All masks of stroke lesions were brought to the right hemisphere (affected hemisphere, AH) and overlaid on T1 template in MNI standard space. Axial slices were chosen with given z-values. The colour bar indicates the number of subjects in which voxels lay within a stroke lesion. UH = unaffected hemisphere.

Table 1

Clinical data

IDAgeGenderDominant hemisphereStroke lesion locationAHTime after stroke (months)Modified Rankin ScaleGrip force (kg)
Pinch force (kg)
UEFMMF
AUAUAU
0169MLTCR97225.050.07.510.047−0.31.2
0259MLPONR15123.745.04.56.058−0.50.0
0358FLPONR20125.029.015.516.0641.91.5
0471MLCRR12137.337.76.36.0660.9−0.5
0574MLTCL33218.038.04.56.046−1.3−0.4
0649FLCRR71119.028.07.06.056−0.3−1.0
0769MRCR, PLICL12147.753.39.711.0662.11.7
0851MLCIR80313.554.06.012.537−1.72.1
0971MLCIR13034.740.08.57.2661.2−0.1
1058MLCR, BG, CIR69210.029.04.010.032−2.40.0
1157MLPLICL20039.052.07.79.7601.01.3
1271MLCRR27230.046.010.012.5490.51.6
1366MLCI, CRR14128.737.48.17.0640.7−0.2
1474FLCR, CIR12126.728.55.86.5660.2−0.9
1572FLCR, CIR36212.024.54.56.039−1.9−1.3
1669FLMOR12023.221.95.46.9660.0−1.2
1763MLPREL28040.044.07.67.5641.20.3
1875MLMEDR15214.038.07.09.046−1.00.3
1963MLPLICL12233.442.48.410.0661.10.8
2059MLPREL26140.639.07.59.2591.00.4
2173FLCR, BGL169118.317.65.04.560−0.6−2.0
2252MLCI, BGL19142.347.312.010.0662.31.1
2346FLMEDL1617.16.73.12.961−1.3−3.1
2466MLCR, PLICL12133.533.87.86.9661.0−0.5
2551FLCI, TCL1239.426.42.16.516−3.5−1.0
IDAgeGenderDominant hemisphereStroke lesion locationAHTime after stroke (months)Modified Rankin ScaleGrip force (kg)
Pinch force (kg)
UEFMMF
AUAUAU
0169MLTCR97225.050.07.510.047−0.31.2
0259MLPONR15123.745.04.56.058−0.50.0
0358FLPONR20125.029.015.516.0641.91.5
0471MLCRR12137.337.76.36.0660.9−0.5
0574MLTCL33218.038.04.56.046−1.3−0.4
0649FLCRR71119.028.07.06.056−0.3−1.0
0769MRCR, PLICL12147.753.39.711.0662.11.7
0851MLCIR80313.554.06.012.537−1.72.1
0971MLCIR13034.740.08.57.2661.2−0.1
1058MLCR, BG, CIR69210.029.04.010.032−2.40.0
1157MLPLICL20039.052.07.79.7601.01.3
1271MLCRR27230.046.010.012.5490.51.6
1366MLCI, CRR14128.737.48.17.0640.7−0.2
1474FLCR, CIR12126.728.55.86.5660.2−0.9
1572FLCR, CIR36212.024.54.56.039−1.9−1.3
1669FLMOR12023.221.95.46.9660.0−1.2
1763MLPREL28040.044.07.67.5641.20.3
1875MLMEDR15214.038.07.09.046−1.00.3
1963MLPLICL12233.442.48.410.0661.10.8
2059MLPREL26140.639.07.59.2591.00.4
2173FLCR, BGL169118.317.65.04.560−0.6−2.0
2252MLCI, BGL19142.347.312.010.0662.31.1
2346FLMEDL1617.16.73.12.961−1.3−3.1
2466MLCR, PLICL12133.533.87.86.9661.0−0.5
2551FLCI, TCL1239.426.42.16.516−3.5−1.0

Age (in years) and gender (M = male; F = female), dominant hemisphere (L = left; R = right), stroke location (MED = media infarct; POS = posterior infarct; PLIC = posterior limb of the internal capsule; CR = corona radiate; CI = capsula interna; TC = thalamocapsular; LC = lenticocapsular; MO = medulla oblongata; BG = basal ganglia; PRE = precentral gyrus; PON = pons), affected hemisphere (AH), time in months after stroke. Absolute grip and pinch force values (in kg), and the UEFM are merged to one composite motor function score (MF) for the affected hand (A). Grip and pinch force values of the unaffected hand (U) were also merged to one composite motor function score.

Table 1

Clinical data

IDAgeGenderDominant hemisphereStroke lesion locationAHTime after stroke (months)Modified Rankin ScaleGrip force (kg)
Pinch force (kg)
UEFMMF
AUAUAU
0169MLTCR97225.050.07.510.047−0.31.2
0259MLPONR15123.745.04.56.058−0.50.0
0358FLPONR20125.029.015.516.0641.91.5
0471MLCRR12137.337.76.36.0660.9−0.5
0574MLTCL33218.038.04.56.046−1.3−0.4
0649FLCRR71119.028.07.06.056−0.3−1.0
0769MRCR, PLICL12147.753.39.711.0662.11.7
0851MLCIR80313.554.06.012.537−1.72.1
0971MLCIR13034.740.08.57.2661.2−0.1
1058MLCR, BG, CIR69210.029.04.010.032−2.40.0
1157MLPLICL20039.052.07.79.7601.01.3
1271MLCRR27230.046.010.012.5490.51.6
1366MLCI, CRR14128.737.48.17.0640.7−0.2
1474FLCR, CIR12126.728.55.86.5660.2−0.9
1572FLCR, CIR36212.024.54.56.039−1.9−1.3
1669FLMOR12023.221.95.46.9660.0−1.2
1763MLPREL28040.044.07.67.5641.20.3
1875MLMEDR15214.038.07.09.046−1.00.3
1963MLPLICL12233.442.48.410.0661.10.8
2059MLPREL26140.639.07.59.2591.00.4
2173FLCR, BGL169118.317.65.04.560−0.6−2.0
2252MLCI, BGL19142.347.312.010.0662.31.1
2346FLMEDL1617.16.73.12.961−1.3−3.1
2466MLCR, PLICL12133.533.87.86.9661.0−0.5
2551FLCI, TCL1239.426.42.16.516−3.5−1.0
IDAgeGenderDominant hemisphereStroke lesion locationAHTime after stroke (months)Modified Rankin ScaleGrip force (kg)
Pinch force (kg)
UEFMMF
AUAUAU
0169MLTCR97225.050.07.510.047−0.31.2
0259MLPONR15123.745.04.56.058−0.50.0
0358FLPONR20125.029.015.516.0641.91.5
0471MLCRR12137.337.76.36.0660.9−0.5
0574MLTCL33218.038.04.56.046−1.3−0.4
0649FLCRR71119.028.07.06.056−0.3−1.0
0769MRCR, PLICL12147.753.39.711.0662.11.7
0851MLCIR80313.554.06.012.537−1.72.1
0971MLCIR13034.740.08.57.2661.2−0.1
1058MLCR, BG, CIR69210.029.04.010.032−2.40.0
1157MLPLICL20039.052.07.79.7601.01.3
1271MLCRR27230.046.010.012.5490.51.6
1366MLCI, CRR14128.737.48.17.0640.7−0.2
1474FLCR, CIR12126.728.55.86.5660.2−0.9
1572FLCR, CIR36212.024.54.56.039−1.9−1.3
1669FLMOR12023.221.95.46.9660.0−1.2
1763MLPREL28040.044.07.67.5641.20.3
1875MLMEDR15214.038.07.09.046−1.00.3
1963MLPLICL12233.442.48.410.0661.10.8
2059MLPREL26140.639.07.59.2591.00.4
2173FLCR, BGL169118.317.65.04.560−0.6−2.0
2252MLCI, BGL19142.347.312.010.0662.31.1
2346FLMEDL1617.16.73.12.961−1.3−3.1
2466MLCR, PLICL12133.533.87.86.9661.0−0.5
2551FLCI, TCL1239.426.42.16.516−3.5−1.0

Age (in years) and gender (M = male; F = female), dominant hemisphere (L = left; R = right), stroke location (MED = media infarct; POS = posterior infarct; PLIC = posterior limb of the internal capsule; CR = corona radiate; CI = capsula interna; TC = thalamocapsular; LC = lenticocapsular; MO = medulla oblongata; BG = basal ganglia; PRE = precentral gyrus; PON = pons), affected hemisphere (AH), time in months after stroke. Absolute grip and pinch force values (in kg), and the UEFM are merged to one composite motor function score (MF) for the affected hand (A). Grip and pinch force values of the unaffected hand (U) were also merged to one composite motor function score.

Brain imaging

A 3 T Siemens Skyra MRI scanner was used to acquire both diffusion-weighted and high-resolution T1-weighted anatomical images. For diffusion-weighted imaging, 75 axial slices were obtained covering the whole brain with gradients (b = 1500 mm2/s) applied along 64 non-collinear directions with the sequence parameters: repetition time = 10 000 ms, echo time = 82 ms, field of view = 256 × 204, slice thickness = 2 mm, in-plane resolution = 2 × 2 mm. The complete data set consisted of 2 × 64 b1500 images and additionally one b0 image at the beginning and one after the first 64 images. For anatomical imaging, a 3D magnetization-prepared, rapid acquisition gradient-echo sequence (MPRAGE) was used with the following parameters: repetition time = 2500 ms, echo time = 2.12 ms, field of view = 240 × 192 mm, 256 axial slices, slice thickness = 0.94 mm, in-plane resolution = 0.94 × 0.94 mm

Preprocessing

The analysis of the diffusion-weighted and anatomical images was conducted using the FSL software package 5.1 (http://www.fmrib.ox.ac.uk/fsl). After correcting for eddy currents and head motion and brain extraction, fractional anisotropy maps were calculated fitting the diffusion tensor model at each voxel. Individual fractional anisotropy maps were then registered non-linearly to the Montreal Neurological Institute (MNI) standard space. The non-linear transformation was also performed with the anatomical images after non-linear co-registration to the individual fractional anisotropy map.

Creation of cortical seed masks

To calculate cortical seed and target masks for the tractography, the T1 structural image was used for brain segmentation into white and grey matter. The Freesurfer image analysis suite (http://surfer.nmr.mgh.harvard.edu/) was used for subsequent automatic cortical parcellation. Surface-based cortical grey matter masks were then transferred to FSL, refined in an automated fashion and used to calculate individual masks related to the boundary between white and grey matter including the whole extent of M1, ventral premotor cortex and the intraparietal sulcus. To bias the tract reconstruction towards the hand representations within each motor area of interest, functional imaging and electrophysiological data were incorporated into the mask creation as recently introduced for the reconstruction of corticospinal (Schulz et al., 2012) and corticocortical connections in young and older participants (Schulz et al., 2014): 500 voxels adjacent to peak coordinates of interest (Supplementary Table 1) within each region were automatically selected using an in-house Matlab script (Matlab R2010 b, Mathworks) (Supplementary material) providing small motor masks within each anatomical region related to hand function and standardized in size and relation to the cortical grey matter/white matter boundary.

Probabilistic tractography of corticocortical and corticospinal trajectories

Probabilistic tractography was used to reconstruct the connections between (i) M1 and ventral premotor cortex; (ii) M1 and anterior and caudal IPS; (iii) ventral premotor cortex and these two parietal brain regions; and also (iv) between anterior and caudal IPS. Therefore, 100 000 streamlines were sent bi-directionally from both the seed and target regions applying exclusion masks to guide this first tractography. The thresholded and refined combination of both resulting output distributions (Hughes et al., 2012) was then used to estimate a tract-specific exclusion mask in which the second tractography was conducted to avoid erroneous trajectories. This procedure was found to allow a reliable reconstruction of the trajectories with—compared to others—even very small structural connectivity probabilities. The probabilistic tractography distribution was finally analysed applying a spectrum of different thresholds from 1% to 20% to estimate the tract-related fractional anisotropy—a measure of white matter integrity for both the affected and unaffected hemisphere. Details on the tractography and the fractional anisotropy calculation can be found in the online Supplementary material. In agreement with previous studies, proportional fractional anisotropy values (ratio affected/unaffected hemisphere) were calculated to assess fractional anisotropy asymmetry and used for the correlative analyses (Schulz et al., 2012). To assess white matter disintegrity in the stroke patients in comparison to healthy controls, 20 healthy age- and gender matched controls were also recruited and underwent brain imaging. To account for the distribution of dominant and non-dominant hemispheres affected, the right and left hemispheres of the control subjects were randomly assigned to the ‘affected’ (right) and ‘unaffected’ (left) hemisphere, respectively. As the corticospinal tract from M1 critically involves residual motor function after stroke (Schulz et al., 2012), we aimed at addressing structure-function relationships for the corticocortical connections between parietal brain region and M1 and premotor cortex while accounting for the integrity of the corticospinal tract. The corticospinal tract reconstruction was conducted in the controls as recently suggested (Schulz et al., 2012), seeding from the hand representation in M1 downwards via the posterior limb of the internal capsule, the pontomedullar junction and the ventral medulla oblongata. After thresholding each output distribution and binarization, the common group average of the corticospinal tract was defined by those voxels that were found in at least 65% of the number of subjects contributing to the tract (Schulz et al., 2012). For the stroke patients, the binarized corticospinal tracts were used to estimate individual fractional anisotropy values at the level from the mesencephalon to the cerebral peduncle (MNI coordinates z = −25 to z = −20). Finally, tract-related fractional anisotropy values were reported as proportional values (ratio affected/unaffected hemisphere). Details on the corticospinal tract reconstruction can be found in the Supplementary material.

Statistics

For within- and between-group comparison of behavioural outcome measures, tract topography and proportional fractional anisotropy values we used Student’s t-tests. False discovery rate (FDR) correction was applied to correct for multiple comparisons. The relationship between the asymmetry of the tract-related fractional anisotropy values and motor function was inferred using generalized linear modelling. Including the proportional fractional anisotropy value of each tract of interest individually, we estimated six models. Corticospinal tract integrity (proportional fractional anisotropy) and motor function of the unaffected hand were included as covariates. Moreover, own observations and previous reports (Harris and Eng, 2006) have suggested that patients with lesions to the dominant hemisphere demonstrate better motor function than those with the non-dominant hand affected. Thus, this information was also included as a categorical factor. Model assumptions were met after z-transformation of the dependent and continuous independent variables. Statistical significance was assumed at P-values ≤ 0.05. All results are given as mean ± standard error of the mean (SEM). Statistical analysis was conducted using SPSS 21 software (IBM Corp.).

Results

Probabilistic tractography of intrahemispheric parietofrontal motor pathways

Probable trajectories connecting M1 with premotor cortex, M1 with anterior and caudal IPS as well as tracts connecting premotor cortex with anterior IPS and caudal IPS were successfully obtained in stroke patients and healthy controls (Fig. 2). More specifically, the fibre tracts between M1 and both premotor cortex and anterior IPS as well as between premotor cortex and anterior IPS were continuously reproducible in all stroke patients. By contrast, connections targeting caudal IPS were much more difficult to reconstruct: for caudal IPS–premotor cortex, three stroke patients had to be excluded from further analysis, for caudal IPS-M1 11 patients, for caudal IPS-anterior IPS two patients due to invalid tract reconstructions. Figure 2 provides an overview of the 3D spatial courses of the different tracts of interest. It shows that for M1 and premotor cortex, there is significant convergence of fibres towards the hand knob area in M1 and the anterior bank of the caudal precentral gyrus for ventral premotor cortex. Contrarily, the spatial variability is increased in deeper white matter with a more distributed course along the fibre tracts. For M1–premotor cortex, the present results are in good agreement with previous diffusion tensor imaging data on corticocortical connectivity in young and older participants (Schulz et al., 2014). In contrast to M1 or premotor cortex, the 3D spatial variability of the trajectories targeting the parietal cortex was higher. This is in line with the considerable inter-subject variability in local (i.e. gyral) brain anatomy of the parietal cortex. The trajectory variability map clearly shows that the fibres towards caudal IPS target an area on the medial wall of the caudal intraparietal sulcus whereas the seed/target region for anterior IPS was located—as hypothesized and calculated—on the anterior lateral aspect of the intraparietal sulcus.

Trajectory variability map for connections of interest in chronic stroke patients. Probable corticocortical connections for the tracts of interest are overlaid on T1 template in MNI standard space on axial slices with z-values given. Colour bar indicates the number of subjects in which voxels are considered part of the tracts of interest. Note that individual subject-specific binarized tracts of interest were used to calculate tract-related fractional anisotropy and not the group average (for illustration only). Tracts in subjects with left-sided lesions were flipped over the mid-sagital plane. UH = unaffected hemisphere; AH = affected hemisphere; aIPS = anterior IPS; cIPS = caudal IPS; PMv = premotor cortex.
Figure 2

Trajectory variability map for connections of interest in chronic stroke patients. Probable corticocortical connections for the tracts of interest are overlaid on T1 template in MNI standard space on axial slices with z-values given. Colour bar indicates the number of subjects in which voxels are considered part of the tracts of interest. Note that individual subject-specific binarized tracts of interest were used to calculate tract-related fractional anisotropy and not the group average (for illustration only). Tracts in subjects with left-sided lesions were flipped over the mid-sagital plane. UH = unaffected hemisphere; AH = affected hemisphere; aIPS = anterior IPS; cIPS = caudal IPS; PMv = premotor cortex.

In addition to the 3D spatial overlap of binarized, individual trajectories, a centre of gravity analysis allowed a more detailed assessment of topographic aspects of the long-range parietofrontal motor pathways contributing to the SLF.

As illustrated in Fig. 3, the four tracts connecting the parietal brain regions with M1 and premotor cortex exhibited spatially distinct courses within SLF II and III. At the mid-level of the parietofrontal connections (y = −34, MNI, Fig. 3), the tracts connecting anterior IPS and premotor cortex were located in a lateral and caudal location, likely within SLF III whereas fibres targeting M1 were located in a more medial and cranial position. Connections from caudal IPS to M1 and caudal IPS to premotor cortex seem to continue this distinct lateromedial and caudocranial topography and belong to SFL II (Fig. 3). Corroborating the 3D results from the trajectory variability maps in Fig. 2, the centre of gravity analysis also showed significant convergence with decreasing distance to the seed/target regions. Tractography was also conducted in healthy participants to compare the distributions between stroke patients and healthy controls. On visual inspection, there were similar parietofrontal connections in both groups (Supplementary Fig. 2). More specifically, stroke patients and controls showed comparable tract-related centres of gravity at the mid-level of the parietofrontal connections (Table 2 and Supplementary Fig. 3 for trajectories of the healthy controls). To correct the influence of the microstructural parameters of the parietofrontal connections to the integrity of the corticospinal tract, we also performed a probabilistic tractography of the corticofugal fibres from M1. These results are presented in the Supplementary Fig. 1.

Centre of gravity analysis of parietofrontal connections in chronic stroke patients. Tract-related centres of gravity were calculated for each subject and each tract and threshold from y = −90 to y = 10 in steps of 2 mm. For the final coordinate, centres of gravity were averaged across the thresholds. Notably, only those y-values were presented in which more than two thresholds contributed to the final coordinate. The centre of gravity distribution is shown on sagital slices individually for the affected (AH) and unaffected (UH) hemisphere, on one horizontal slice and one coronar slice at y = −34 in MNI standard space. Table 2 provides statistics on the centre of gravity analysis at the coronar level. Please note that the caudal course from M1 from y = −20 towards higher y-values is caused by individual variability of fibres trajectories. aIPS = anterior IPS; cIPS = caudal IPS; PMv = premotor cortex.
Figure 3

Centre of gravity analysis of parietofrontal connections in chronic stroke patients. Tract-related centres of gravity were calculated for each subject and each tract and threshold from y = −90 to y = 10 in steps of 2 mm. For the final coordinate, centres of gravity were averaged across the thresholds. Notably, only those y-values were presented in which more than two thresholds contributed to the final coordinate. The centre of gravity distribution is shown on sagital slices individually for the affected (AH) and unaffected (UH) hemisphere, on one horizontal slice and one coronar slice at y = −34 in MNI standard space. Table 2 provides statistics on the centre of gravity analysis at the coronar level. Please note that the caudal course from M1 from y = −20 towards higher y-values is caused by individual variability of fibres trajectories. aIPS = anterior IPS; cIPS = caudal IPS; PMv = premotor cortex.

Table 2

Tract-related centres of gravity at y = −34 (MNI) in stroke patients and healthy controls

Anterior IPS-premotor cortex
Caudal IPS-premotor cortex
Anterior IPS-M1
Caudal IPS-M1
AHUHAHUHAHUHAHUH
xStroke38.4 (± 0.55)−38.5 (± 0.73)27.5 (± 1.06)−25.7 (±0.83)32.6 (±0.83)−30.8 (±0.5)26.4 (±0.5)−25.8 (±0.67)
Control37.5 (±0.74)−37.4 (± 0.87)25.5 (±1.01)−25.1 (±0.9)31.0 (±0.86)−31.3 (±0.69)25.9 (±0.62)−24.6 (±0.71)
P-value0.650.650.650.860.650.760.760.65
zStroke35.4 (± 0.31)35.9 (± 0.45)41.4 (±0.65)43.5 (±0.8)39.2 (±0.62)39.6 (±0.53)42.7 (±1.03)44.2 (±1.17)
Control35.6 (± 0.25)35.8 (±0.6)43.0 (±1.26)44.3 (±0.8)40.2 (±0.8)39.4 (±0.66)45.0 (±0.98)46.0 (±0.84)
P-value0.860.980.650.720.650.860.650.65
Anterior IPS-premotor cortex
Caudal IPS-premotor cortex
Anterior IPS-M1
Caudal IPS-M1
AHUHAHUHAHUHAHUH
xStroke38.4 (± 0.55)−38.5 (± 0.73)27.5 (± 1.06)−25.7 (±0.83)32.6 (±0.83)−30.8 (±0.5)26.4 (±0.5)−25.8 (±0.67)
Control37.5 (±0.74)−37.4 (± 0.87)25.5 (±1.01)−25.1 (±0.9)31.0 (±0.86)−31.3 (±0.69)25.9 (±0.62)−24.6 (±0.71)
P-value0.650.650.650.860.650.760.760.65
zStroke35.4 (± 0.31)35.9 (± 0.45)41.4 (±0.65)43.5 (±0.8)39.2 (±0.62)39.6 (±0.53)42.7 (±1.03)44.2 (±1.17)
Control35.6 (± 0.25)35.8 (±0.6)43.0 (±1.26)44.3 (±0.8)40.2 (±0.8)39.4 (±0.66)45.0 (±0.98)46.0 (±0.84)
P-value0.860.980.650.720.650.860.650.65

Centres of gravity of parietofrontal connections targeting M1 and premotor cortex at the mid-level of the connections (y = −34) for stroke patients and controls. x and z values are presented individually for the affected (AH) and unaffected (UH) hemisphere. Results are given as mean ± SEM. Individual Student’s t-tests were used to compare tract-related coordinates between stroke patients and controls. P-values were fully FDR-corrected for multiple comparisons.

Table 2

Tract-related centres of gravity at y = −34 (MNI) in stroke patients and healthy controls

Anterior IPS-premotor cortex
Caudal IPS-premotor cortex
Anterior IPS-M1
Caudal IPS-M1
AHUHAHUHAHUHAHUH
xStroke38.4 (± 0.55)−38.5 (± 0.73)27.5 (± 1.06)−25.7 (±0.83)32.6 (±0.83)−30.8 (±0.5)26.4 (±0.5)−25.8 (±0.67)
Control37.5 (±0.74)−37.4 (± 0.87)25.5 (±1.01)−25.1 (±0.9)31.0 (±0.86)−31.3 (±0.69)25.9 (±0.62)−24.6 (±0.71)
P-value0.650.650.650.860.650.760.760.65
zStroke35.4 (± 0.31)35.9 (± 0.45)41.4 (±0.65)43.5 (±0.8)39.2 (±0.62)39.6 (±0.53)42.7 (±1.03)44.2 (±1.17)
Control35.6 (± 0.25)35.8 (±0.6)43.0 (±1.26)44.3 (±0.8)40.2 (±0.8)39.4 (±0.66)45.0 (±0.98)46.0 (±0.84)
P-value0.860.980.650.720.650.860.650.65
Anterior IPS-premotor cortex
Caudal IPS-premotor cortex
Anterior IPS-M1
Caudal IPS-M1
AHUHAHUHAHUHAHUH
xStroke38.4 (± 0.55)−38.5 (± 0.73)27.5 (± 1.06)−25.7 (±0.83)32.6 (±0.83)−30.8 (±0.5)26.4 (±0.5)−25.8 (±0.67)
Control37.5 (±0.74)−37.4 (± 0.87)25.5 (±1.01)−25.1 (±0.9)31.0 (±0.86)−31.3 (±0.69)25.9 (±0.62)−24.6 (±0.71)
P-value0.650.650.650.860.650.760.760.65
zStroke35.4 (± 0.31)35.9 (± 0.45)41.4 (±0.65)43.5 (±0.8)39.2 (±0.62)39.6 (±0.53)42.7 (±1.03)44.2 (±1.17)
Control35.6 (± 0.25)35.8 (±0.6)43.0 (±1.26)44.3 (±0.8)40.2 (±0.8)39.4 (±0.66)45.0 (±0.98)46.0 (±0.84)
P-value0.860.980.650.720.650.860.650.65

Centres of gravity of parietofrontal connections targeting M1 and premotor cortex at the mid-level of the connections (y = −34) for stroke patients and controls. x and z values are presented individually for the affected (AH) and unaffected (UH) hemisphere. Results are given as mean ± SEM. Individual Student’s t-tests were used to compare tract-related coordinates between stroke patients and controls. P-values were fully FDR-corrected for multiple comparisons.

Tract-related white matter integrity and motor function after stroke

Table 3 summarizes the tract-related absolute fractional anisotropy values for the affected and unaffected hemispheres of the stroke patients and the data for the healthy controls. Individual Student’s t-tests did not reveal significant differences between the groups (Supplementary Fig. 4A). This was also shown by an additional repeated-measures ANOVA including the subset of 12 stroke patients and 19 controls in which all tracts were successfully reconstructed (Supplementary material). The proportional fractional anisotropy values were also comparable between stroke patients and controls (Supplementary Fig. 4B).

Table 3

Tract-related white matter integrity in stroke patients and healthy controls

Premotor cortex-M1
Anterior IPS-premotor cortex
Caudal IPS-premotor cortex
Anterior IPS-M1
Caudal IPS-M1
Caudal IPS-anterior IPS
AHUHAHUHAHUHAHUHAHUHAHUH
Fractional anisotropyStroke0.35 (±0.01)0.38 (±0.005)0.35 (±0.013)0.38 (±0.006)0.38 (±0.011)0.4 (±0.007)0.37 (±0.01)0.39 (±0.006)0.39 (±0.009)0.4 (±0.009)0.34 (±0.01)0.36 (±0.006)
Control0.37 (±0.007)0.38 (±0.007)0.38 (±0.009)0.38 (±0.01)0.4 (±0.01)0.39 (±0.009)0.39 (±0.007)0.39 (±0.008)0.4 (±0.006)0.4 (±0.008)0.36 (±0.009)0.35 (±0.009)
P-value0.290.960.290.960.290.770.290.770.470.960.330.96
Premotor cortex-M1
Anterior IPS-premotor cortex
Caudal IPS-premotor cortex
Anterior IPS-M1
Caudal IPS-M1
Caudal IPS-anterior IPS
AHUHAHUHAHUHAHUHAHUHAHUH
Fractional anisotropyStroke0.35 (±0.01)0.38 (±0.005)0.35 (±0.013)0.38 (±0.006)0.38 (±0.011)0.4 (±0.007)0.37 (±0.01)0.39 (±0.006)0.39 (±0.009)0.4 (±0.009)0.34 (±0.01)0.36 (±0.006)
Control0.37 (±0.007)0.38 (±0.007)0.38 (±0.009)0.38 (±0.01)0.4 (±0.01)0.39 (±0.009)0.39 (±0.007)0.39 (±0.008)0.4 (±0.006)0.4 (±0.008)0.36 (±0.009)0.35 (±0.009)
P-value0.290.960.290.960.290.770.290.770.470.960.330.96

We estimated fractional anisotropy in individual white matter tracts of these six connections in the affected (AH) and unaffected (UH) hemisphere of stroke patients and healthy control subjects. Results are given as mean ± SEM. Comparing stroke patients with healthy controls we performed individual Student’s t-tests. The P-values were fully FDR corrected for multiple testing.

Table 3

Tract-related white matter integrity in stroke patients and healthy controls

Premotor cortex-M1
Anterior IPS-premotor cortex
Caudal IPS-premotor cortex
Anterior IPS-M1
Caudal IPS-M1
Caudal IPS-anterior IPS
AHUHAHUHAHUHAHUHAHUHAHUH
Fractional anisotropyStroke0.35 (±0.01)0.38 (±0.005)0.35 (±0.013)0.38 (±0.006)0.38 (±0.011)0.4 (±0.007)0.37 (±0.01)0.39 (±0.006)0.39 (±0.009)0.4 (±0.009)0.34 (±0.01)0.36 (±0.006)
Control0.37 (±0.007)0.38 (±0.007)0.38 (±0.009)0.38 (±0.01)0.4 (±0.01)0.39 (±0.009)0.39 (±0.007)0.39 (±0.008)0.4 (±0.006)0.4 (±0.008)0.36 (±0.009)0.35 (±0.009)
P-value0.290.960.290.960.290.770.290.770.470.960.330.96
Premotor cortex-M1
Anterior IPS-premotor cortex
Caudal IPS-premotor cortex
Anterior IPS-M1
Caudal IPS-M1
Caudal IPS-anterior IPS
AHUHAHUHAHUHAHUHAHUHAHUH
Fractional anisotropyStroke0.35 (±0.01)0.38 (±0.005)0.35 (±0.013)0.38 (±0.006)0.38 (±0.011)0.4 (±0.007)0.37 (±0.01)0.39 (±0.006)0.39 (±0.009)0.4 (±0.009)0.34 (±0.01)0.36 (±0.006)
Control0.37 (±0.007)0.38 (±0.007)0.38 (±0.009)0.38 (±0.01)0.4 (±0.01)0.39 (±0.009)0.39 (±0.007)0.39 (±0.008)0.4 (±0.006)0.4 (±0.008)0.36 (±0.009)0.35 (±0.009)
P-value0.290.960.290.960.290.770.290.770.470.960.330.96

We estimated fractional anisotropy in individual white matter tracts of these six connections in the affected (AH) and unaffected (UH) hemisphere of stroke patients and healthy control subjects. Results are given as mean ± SEM. Comparing stroke patients with healthy controls we performed individual Student’s t-tests. The P-values were fully FDR corrected for multiple testing.

For each tract of interest, one generalized linear model was estimated to infer the association between white matter integrity and motor function whose predictors are summarized by forest plots in Fig. 4 (see Supplementary Table 2 for statistical details). We found significant positive contributions of the microstructural integrity of the fibre tracts connecting ventral premotor cortex with M1 (B = 0.352; P < 0.001) and also anterior IPS with premotor cortex (B = 0.298; P < 0.01) for the motor function of the affected hand. By contrast, the microstructural integrity of the connections between anterior IPS and M1 (B = 0.213; P = 0.07) as well caudal IPS and premotor cortex (B = 0.287; P = 0.07), M1 (B = −0.068; P = 0.82) and anterior IPS (B = 0.003; P = 0.99) did not show significant associations with the preserved motor function after stroke (Fig. 4). Notably, across all parietofrontal connections both the integrity of the corticospinal tracts and the motor function of the unaffected hand were relevant predictors of residual motor function of the affected hand after stroke.

Tract-related white matter integrity of parietofrontal connections and motor function after stroke. Forest plots of the Beta values (± 95% confidence intervals) illustrate the statistical influence of proportional fractional anisotropy values for each tract of interest for residual motor function after stroke. Applying individual generalized linear models, the influence of each tract was adjusted for three covariates: white matter integrity (proportional fractional anisotropy values) of the corticospinal tract (CST), the affected side (SIDE, dominant versus non-dominant hemisphere) and the motor function (MF) of the unaffected hand. aIPS = anterior IPS; cIPS = caudal IPS; PMv = premotor cortex. ***P < 0.001; **P < 0.01; *P < 0.05.
Figure 4

Tract-related white matter integrity of parietofrontal connections and motor function after stroke. Forest plots of the Beta values (± 95% confidence intervals) illustrate the statistical influence of proportional fractional anisotropy values for each tract of interest for residual motor function after stroke. Applying individual generalized linear models, the influence of each tract was adjusted for three covariates: white matter integrity (proportional fractional anisotropy values) of the corticospinal tract (CST), the affected side (SIDE, dominant versus non-dominant hemisphere) and the motor function (MF) of the unaffected hand. aIPS = anterior IPS; cIPS = caudal IPS; PMv = premotor cortex. ***P < 0.001; **P < 0.01; *P < 0.05.

Discussion

White-matter integrity of parietofrontal motor pathways and motor function after stroke

As a main finding of the present study, tract-related proportional fractional anisotropy as a marker for microstructural white matter integrity showed a significant association with basic residual motor function of the affected hand for M1–ventral premotor cortex and anterior IPS–ventral premotor cortex connections. Thus, the present analysis provides first probabilistic structural connectivity data in chronic stroke patients which indicate that, aside from interactions with M1, i.e. between ventral premotor cortex and M1, also interactions between secondary parietal and frontal motor areas (anterior IPS and ventral premotor cortex) influence residual motor function after stroke. These data extend prior knowledge about motor-related parietofrontal interactions, both in stroke patients and healthy participants.

So far, some studies have contributed to the understanding of M1, ventral premotor cortex and the posterior parietal cortex as important areas for motor recovery and residual motor function after stroke. For instance, animal studies showed that ventral premotor cortex was able to increase its connections to the somatosensory cortex (Dancause et al., 2005) or enlarge its hand representation after a lesion to M1 (Frost et al., 2003). In chronic stroke patients, a meta-analysis has recently evidenced an increase in task-related premotor cortex activation both on the ipsilesional and contralesional hemisphere compared to controls (Rehme et al., 2012). Early after stroke, a positive coupling onto the ipsilesional M1 was found to be significantly reduced in the stroke patients, which could be related to the residual motor function (Rehme et al., 2011). Providing structural connectivity data for M1-premotor cortex, the present study strengthens these functional data by showing that the integrity of the underlying white matter tracts connecting M1 with ventral premotor cortex also influences the residual motor function.

In terms of the functional role of parietal brain regions for recovery and motor function after stroke, the available data are still limited. We found that the integrity of the anterior IPS–premotor cortex connection appears to be relevant for hand motor function after stroke. This importance of anterior IPS was similarly suggested by a meta-analysis on functional imaging data in chronic stroke patients, which found a consistent task-related activation of anterior IPS during various movements of the affected limb (Rehme et al., 2012). Also structural analyses on cortical strokes have demonstrated that deficits in parietal sensorimotor integration areas significantly impact on residual motor function and recovery after stroke: for instance, one longitudinal study has found that parietal cortex lesions impaired the overall outcome while lesions to M1 influenced rather the dynamics of motor recovery after stroke (Abela et al., 2012). Another study reported that grey matter atrophy in a cluster on the angularis gyrus was associated with less favourable recovery after cortical sensorimotor stroke (Abela et al., 2014). Providing structural data, the present results corroborate previous functional connectivity data: for instance it was shown that increased functional connectivity between ipsilesional M1 and ipsilesional inferior parietal lobe was associated with improved motor recovery whereas a decreased parietofrontal connectivity was related to bad outcome (Yin et al., 2012). Similar findings were also obtained for functional connectivity between the superior parietal cortex and M1 suggesting that motor control deficits following stroke might be partly caused by a disconnection between parietal higher-order motor guidance systems and the primary motor network (Inman et al., 2012). Another study (Park et al., 2011) has addressed recovery-related changes of altered parietofrontal functional connectivity. The authors found a reduced connectivity between contralesional posterior parietal cortex and M1 shortly after stroke. In the early subacute stage 1 and 6 months after the ischaemic event the connectivity between M1 and ipsilesional posterior parietal brain regions was increased while the reduction of connectivity with the contralesional posterior parietal cortex was persistent. However, an association with behaviour was not detected. Adding structural connectivity data for the connection between anterior IPS and premotor cortex, the present study strengthens the view that posterior parietal areas might contribute to motor function after stroke via corticocortical interactions with frontal brain regions such as premotor cortex.

Previous studies in healthy participants have commonly related M1, ventral premotor cortex and anterior IPS-related brain activity and their interactions to rather skilled hand functions, such as visually guided precision grip (Ehrsson et al., 2001), power grip (Ehrsson et al., 2000) and reaching and grasping (Ehrsson et al., 2001; Witt et al., 2008; Davare et al., 2010, 2011). Specifically, ventral premotor cortex activity concerns the integration of information processed by parietal secondary motor areas such as anterior IPS into more M1-related executive information processing. In this context, anterior IPS was found to contribute also to object manipulation and exploration (Binkofski et al., 1999) by continuously converting spatial specification of target location into motor plans (Olivier et al., 2007; Davare et al., 2010, 2011). Thereby, anterior IPS provides premotor cortex with visual information about object properties and influences grasp-related and muscle-specific premotor cortex-M1 interactions (Tunik et al., 2005; Olivier et al., 2007; Koch et al., 2008; Davare et al., 2009, 2010; Reichenbach et al., 2011; Karabanov et al., 2013). The present data suggest that, different from healthy participants, in chronic stroke patients, this structure-behaviour relationship is not only observed in such complex grasp- or reach-related functions but already in basic hand and arm motor functions.

The question arises why we did not find similar associations between tract-related white matter integrity and motor function for the other four pathways of interest. For the direct connections between caudal IPS and M1 the tract reconstruction was difficult, though not caused by a specific co-localization of the tract courses and the stroke lesions (Supplementary material and Supplementary Table 3). Koch et al. (2010) found only sparse connections between the angularis gyrus and M1 in healthy participants. Given the positive findings for functional connectivity between parietal brain regions and M1 (Koch et al., 2007, 2008; Karabanov et al., 2013; Chao et al., 2015) this might indicate that the connectivity from the parietal areas of interest with M1 is probably mediated indirectly via premotor cortex rather than directly. This was suggested by a recent study which combined transcranial magnetic stimulation and diffusion-tensor imaging in healthy participants (Koch et al., 2010) and might explain why anterior IPS–M1 and caudal IPS–M1 did not show significant relationships with motor function. For caudal IPS-related connections, the absent association with motor function could also be explained by its predominant role in integrating visual information in grasping and reaching movements and tool use (Faillenot et al., 1997; Shikata et al., 2001; Culham et al., 2003; Tsutsui et al., 2003; Filimon et al., 2007; Koch et al., 2010; Mruczek et al., 2013). Therefore the present structural data suggest that caudal portions of the intraparietal sulcus do not relate to rather basic hand and arm function after stroke. In this context future studies are needed to further assess possible associations between parietofrontal connections and more complex movements in stroke patients such as reaching and grasping in space.

Generalized linear models were used to infer tract-related structure-behaviour relationships of corticocortical parietofrontal connection in chronic stroke patients while adjusting the associations for the integrity of the corticospinal tract. In agreement with previous work, this measure was found to critically contribute to motor function after stroke (Schulz et al., 2012). Aside from the corticospinal integrity, we found that the integrity between M1 and premotor cortex as well as anterior IPS and premotor cortex explained additional variance in motor function. The positive influence of the third covariate, whether the dominant or non-dominant hand was affected by the stroke, was also in agreement with one previous report (Harris and Eng, 2006), likely caused by more intensive training and use of the dominant hand.

Topography of intrahemispheric parietofrontal motor pathways

Probable intrahemispheric trajectories connecting hand representations of anterior IPS, caudal IPS, ventral premotor cortex and M1 were successfully reconstructed in chronic stroke patients. Similar results with comparable individual trajectories were also obtained in a group of age- and gender-matched healthy controls suggesting valid tract reconstruction in patients, even with cortical stroke locations. A centre of gravity analysis revealed detailed insights into a distinct topography of intrahemispheric long-range fibres connecting anterior and caudal IPS with both premotor cortex and M1. We found a lateromedial and caudocranial distribution of the fibres in the order anterior IPS–premotor cortex, anterior IPS–M1, caudal IPS–M1 and caudal IPS–premotor cortex localizing with distinct components of the SLF. Thereby, anterior IPS–premotor cortex and anterior IPS–M1 were most likely localized in the SLF III component at the laterocaudal extent of the distribution. This localization is in good agreement with previous imaging data in healthy humans (Croxson et al., 2005; Makris et al., 2005; Koch et al., 2010) and anatomical animal work suggesting that SLF III predominantly originates from the supramarginal gyrus and connects prominently to ventral premotor cortex and prefrontal brain regions (Petrides and Pandya, 1984; Schmahmann et al., 2007). In addition we also found a connection to M1. By contrast, a recent probabilistic tractography analysis between the supramarginal gyrus and both premotor cortex and M1 showed only very sparse connections targeting M1 (Koch et al., 2010). Moreover, in animals, connections to M1 were observed even only for the superior parietal lobe (Johnson et al., 1996). Tract reconstructions for caudal IPS-related connections were found to localize in SLF II (Croxson et al., 2005; Makris et al., 2005), the major component of the SLF. SLF II originates predominantly from areas around the angularis gyrus and targets not only dorsolateral frontal brain regions (Petrides and Pandya, 1984; Luppino et al., 1999; Makris et al., 2005; Rozzi et al., 2006; Schmahmann et al., 2007) but also ventral premotor cortex (Croxson et al., 2005). However, compared to anterior IPS, the tract reconstruction between caudal IPS and M1 was less reliable, as previosuly reported (Koch et al., 2010).

In summary, the present analysis adds further structural probabilistic data on parietofrontal connectivity in the human brain, both in healthy participants and stroke patients. Previous imaging studies on parietofrontal connectivity and SLF topography have either used visual image inspection on relative location and orientation of diffusion properties of the tissue for each SLF component (Makris et al., 2005) or probabilistic tractography from large seed and target masks (Koch et al., 2010). Here we performed individual probabilistic tractography from and to rather circumscribed hand representations. This approach of using functional imaging or electrophysiology-based hand representations within functional brain regions for probabilistic tractography has been recently introduced for probabilistic tractography of corticospinal (Schulz et al., 2012) and corticocortical connections (Schulz et al., 2014). As similar trajectories and measures of integrity were obtained in 25 stroke patients and 20 healthy controls, we conclude that the actual trajectories are of reasonable quality to be used for correlative analysis.

Limitations of the study

Some limitations of the present study have to be mentioned. First, we found significant associations between a composite motor function score and proportional fractional anisotropy values within the tracts of interest. The question arises whether this relationship is stroke-specific or a general age-related phenomenon also found in healthy controls. Unfortunately, in the present control sample we did not have behavioural data. Hence we were not able to answer this question. However, a recent analysis of multiple corticocortical connections including M1–premotor cortex in healthy aged participants did not find a significant association between tract-related white matter integrity and motor performance (Schulz et al., 2014). These data favour the view that the present findings might be rather stroke-specific than a general phenomenon. Future studies will have to include correlative analyses in controls to investigate the stroke specificity of the present structural findings. Second, we used probabilistic tractography to reconstruct individual fibre tracts between M1, ventral premotor cortex, anterior IPS and caudal IPS. Seed and target masks were individually tailored according to the subject’s brain anatomy. Exclusion and waypoint masks were used to guide the individual tract reconstruction. Still, the courses of the tracts remain probabilistic and might differ in the individual brain. Third, previous functional imaging and electrophysiological data published were used to bias the tract-reconstruction to probable hand representations. To what extent cortical reorganization after stroke might change gross cortical representations and therefore influence the courses of the corticocortical connections remains unclear. The combination of functional and diffusion tensor imaging in the individual participant might help to improve the estimation of task-related seed and target regions and the reconstruction of the tracts. Fourth, the aim of the present work was to investigate the association between white matter integrity of corticocortical parietofrontal pathways and residual motor function after stroke, adjusting to the integrity of the corticospinal tract. It has to be stated that, despite the amount of variance explained, the corticocortical network analysed here was limited to areas of interest, M1, ventral premotor cortex, anterior IPS and caudal IPS, and thus not complete. More complex modelling on larger sample sizes is needed to assess the involvement of parietofrontal structural connectivity in the context of the structural properties of extended neuronal subcortico-cortico-cortical circuits, including also prefrontal, subcortical and cerebellar structures.

Conclusion

In conclusion, the present results provide structural connectivity data in chronic stroke patients which indicate that the integrity of the connection (i) between ventral premotor cortex and M1; and (ii) between anterior IPS and ventral premotor cortex contribute to motor function after stroke. This supports the view that not only frontal secondary motor areas but also parietal brain regions along the IPS might play a relevant role in generating basic residual motor output after stroke.

Funding

This research was supported by the German Research Foundation (SFB 936-C4 to F.C.H., SFB 936-C1 to C.G. and SFB936-C2 to G.T.) and the Werner-Otto-Foundation (to F.C.H).

Supplementary material

Supplementary material is available at Brain online.

Abbreviations

    Abbreviations
     
  • IPS

    intraparietal sulcus

  •  
  • SLF

    superior longitudinal fascicle

  •  
  • UEFM

    Fugl-Meyer assessment of the upper extremity

References

Abela
E
Missimer
J
Wiest
R
Federspiel
A
Hess
C
Sturzenegger
M
, et al. 
Lesions to primary sensory and posterior parietal cortices impair recovery from hand paresis after stroke
PLoS One
2012
, vol. 
7
 pg. 
e31275
 
Abela
E
Seiler
A
Missimer
JH
Federspiel
A
Hess
CW
Sturzenegger
M
, et al. 
Grey matter volumetric changes related to recovery from hand paresis after cortical sensorimotor stroke
Brain Struct Funct
2014
 
Advance Access published on June 7, 2014
Binkofski
F
Buccino
G
Posse
S
Seitz
RJ
Rizzolatti
G
Freund
H
A fronto-parietal circuit for object manipulation in man: evidence from an fMRI-study
Eur J Neurosci
1999
, vol. 
11
 (pg. 
3276
-
86
)
Borra
E
Belmalih
A
Calzavara
R
Gerbella
M
Murata
A
Rozzi
S
, et al. 
Cortical connections of the macaque anterior intraparietal (AIP) area
Cereb Cortex
2008
, vol. 
18
 (pg. 
1094
-
111
)
Cavada
C
Goldman-Rakic
PS
Posterior parietal cortex in rhesus monkey: II. Evidence for segregated corticocortical networks linking sensory and limbic areas with the frontal lobe
J Comp Neurol
1989
, vol. 
287
 (pg. 
422
-
45
)
Chao
C-C
Karabanov
AN
Paine
R
Carolina de Campos
A
Kukke
SN
Wu
T
, et al. 
Induction of motor associative plasticity in the posterior parietal cortex-primary motor network
Cereb Cortex
2015
, vol. 
25
 (pg. 
365
-
73
)
Croxson
PL
Johansen-Berg
H
Behrens
TEJ
Robson
MD
Pinsk
MA
Gross
CG
, et al. 
Quantitative investigation of connections of the prefrontal cortex in the human and macaque using probabilistic diffusion tractography
J Neurosci
2005
, vol. 
25
 (pg. 
8854
-
66
)
Culham
JC
Danckert
SL
DeSouza
JFX
Gati
JS
Menon
RS
Goodale
MA
Visually guided grasping produces fMRI activation in dorsal but not ventral stream brain areas
Exp Brain Res
2003
, vol. 
153
 (pg. 
180
-
9
)
Dafotakis
M
Sparing
R
Eickhoff
SB
Fink
GR
Nowak
DA
On the role of the ventral premotor cortex and anterior intraparietal area for predictive and reactive scaling of grip force
Brain Res
2008
, vol. 
1228
 (pg. 
73
-
80
)
Dancause
N
Barbay
S
Frost
SB
Plautz
EJ
Chen
D
Zoubina
EV
, et al. 
Extensive cortical rewiring after brain injury
J Neurosci
2005
, vol. 
25
 (pg. 
10167
-
79
)
Davare
M
Kraskov
A
Rothwell
JC
Lemon
RN
Interactions between areas of the cortical grasping network
Curr Opin Neurobiol
2011
, vol. 
21
 (pg. 
565
-
70
)
Davare
M
Montague
K
Olivier
E
Rothwell
JC
Lemon
RN
Ventral premotor to primary motor cortical interactions during object-driven grasp in humans
Cortex
2009
, vol. 
45
 (pg. 
1050
-
7
)
Davare
M
Rothwell
JC
Lemon
RN
Causal connectivity between the human anterior intraparietal area and premotor cortex during grasp
Curr Biol
2010
, vol. 
20
 (pg. 
176
-
81
)
Ehrsson
HH
Fagergren
A
Jonsson
T
Westling
G
Johansson
RS
Forssberg
H
Cortical activity in precision- versus power-grip tasks: an fMRI study
J Neurophysiol
2000
, vol. 
83
 (pg. 
528
-
36
)
Ehrsson
HH
Fagergren
E
Forssberg
H
Differential fronto-parietal activation depending on force used in a precision grip task: an fMRI study
J Neurophysiol
2001
, vol. 
85
 (pg. 
2613
-
23
)
Faillenot
I
Toni
I
Decety
J
Grégoire
MC
Jeannerod
M
Visual pathways for object-oriented action and object recognition: functional anatomy with PET
Cereb Cortex
1997
, vol. 
7
 (pg. 
77
-
85
)
Filimon
F
Nelson
JD
Hagler
DJ
Sereno
MI
Human cortical representations for reaching: mirror neurons for execution, observation, and imagery
Neuroimage
2007
, vol. 
37
 (pg. 
1315
-
28
)
Frost
SB
Barbay
S
Friel
KM
Plautz
EJ
Nudo
RJ
Reorganization of remote cortical regions after ischemic brain injury: a potential substrate for stroke recovery
J Neurophysiol
2003
, vol. 
89
 (pg. 
3205
-
14
)
Gharbawie
OA
Stepniewska
I
Qi
H
Kaas
JH
Multiple parietal-frontal pathways mediate grasping in macaque monkeys
J Neurosci
2011
, vol. 
31
 (pg. 
11660
-
77
)
Grefkes
C
Fink
GR
The functional organization of the intraparietal sulcus in humans and monkeys
J Anat
2005
, vol. 
207
 (pg. 
3
-
17
)
Grefkes
C
Nowak
DA
Eickhoff
SB
Dafotakis
M
Küst
J
Karbe
H
, et al. 
Cortical connectivity after subcortical stroke assessed with functional magnetic resonance imaging
Ann Neurol
2008
, vol. 
63
 (pg. 
236
-
46
)
Grefkes
C
Nowak
DA
Wang
LE
Dafotakis
M
Eickhoff
SB
Fink
GR
Modulating cortical connectivity in stroke patients by rTMS assessed with fMRI and dynamic causal modeling
Neuroimage
2010
, vol. 
50
 (pg. 
233
-
42
)
Harris
JE
Eng
JJ
Individuals with the dominant hand affected following stroke demonstrate less impairment than those with the nondominant hand affected
Neurorehabil Neural Repair
2006
, vol. 
20
 (pg. 
380
-
9
)
Hughes
EJ
Bond
J
Svrckova
P
Makropoulos
A
Ball
G
Sharp
DJ
, et al. 
Regional changes in thalamic shape and volume with increasing age
Neuroimage
2012
, vol. 
63
 (pg. 
1134
-
42
)
Inman
CS
James
GA
Hamann
S
Rajendra
JK
Pagnoni
G
Butler
AJ
Altered resting-state effective connectivity of fronto-parietal motor control systems on the primary motor network following stroke
Neuroimage
2012
, vol. 
59
 (pg. 
227
-
37
)
Johnson
PB
Ferraina
S
Bianchi
L
Caminiti
R
Cortical networks for visual reaching: physiological and anatomical organization of frontal and parietal lobe arm regions
Cereb Cortex
1996
, vol. 
6
 (pg. 
102
-
19
)
Karabanov
AN
Chao
C-C
Paine
R
Hallett
M
Mapping different intra-hemispheric parietal-motor networks using twin Coil TMS
Brain Stimul
2013
, vol. 
6
 (pg. 
384
-
9
)
Koch
G
Cercignani
M
Pecchioli
C
Versace
V
Oliveri
M
Caltagirone
C
, et al. 
In vivo definition of parieto-motor connections involved in planning of grasping movements
Neuroimage
2010
, vol. 
51
 (pg. 
300
-
12
)
Koch
G
Fernandez Del Olmo
M
Cheeran
B
Ruge
D
Schippling
S
Caltagirone
C
, et al. 
Focal stimulation of the posterior parietal cortex increases the excitability of the ipsilateral motor cortex
J Neurosci
2007
, vol. 
27
 (pg. 
6815
-
22
)
Koch
G
Fernandez Del Olmo
M
Cheeran
B
Schippling
S
Caltagirone
C
Driver
J
, et al. 
Functional interplay between posterior parietal and ipsilateral motor cortex revealed by twin-coil transcranial magnetic stimulation during reach planning toward contralateral space
J Neurosci
2008
, vol. 
28
 (pg. 
5944
-
53
)
Luppino
G
Murata
A
Govoni
P
Matelli
M
Largely segregated parietofrontal connections linking rostral intraparietal cortex (areas AIP and VIP) and the ventral premotor cortex (areas F5 and F4)
Exp Brain Res
1999
, vol. 
128
 (pg. 
181
-
7
)
Makris
N
Kennedy
DN
McInerney
S
Sorensen
AG
Wang
R
Caviness
VS
, et al. 
Segmentation of subcomponents within the superior longitudinal fascicle in humans: a quantitative, in vivo, DT-MRI study
Cereb Cortex
2005
, vol. 
15
 (pg. 
854
-
69
)
Mruczek
REB
von Loga
IS
Kastner
S
The representation of tool and non-tool object information in the human intraparietal sulcus
J Neurophysiol
2013
, vol. 
109
 (pg. 
2883
-
96
)
Olivier
E
Davare
M
Andres
M
Fadiga
L
Precision grasping in humans: from motor control to cognition
Curr Opin Neurobiol
2007
, vol. 
17
 (pg. 
644
-
8
)
Park
C
Chang
WH
Ohn
SH
Kim
ST
Bang
OY
Pascual-Leone
A
, et al. 
Longitudinal changes of resting-state functional connectivity during motor recovery after stroke
Stroke
2011
, vol. 
42
 (pg. 
1357
-
62
)
Petrides
M
Pandya
DN
Projections to the frontal cortex from the posterior parietal region in the rhesus monkey
J Comp Neurol
1984
, vol. 
228
 (pg. 
105
-
16
)
Rehme
AK
Eickhoff
SB
Rottschy
C
Fink
GR
Grefkes
C
Activation likelihood estimation meta-analysis of motor-related neural activity after stroke
Neuroimage
2012
, vol. 
59
 (pg. 
2771
-
82
)
Rehme
AK
Eickhoff
SB
Wang
LE
Fink
GR
Grefkes
C
Dynamic causal modeling of cortical activity from the acute to the chronic stage after stroke
Neuroimage
2011
, vol. 
55
 (pg. 
1147
-
58
)
Rehme
AK
Grefkes
C
Cerebral network disorders after stroke: evidence from imaging-based connectivity analyses of active and resting brain states in humans
J Physiol
2013
, vol. 
591
 (pg. 
17
-
31
)
Reichenbach
A
Bresciani
J-P
Peer
A
Bülthoff
HH
Thielscher
A
Contributions of the PPC to online control of visually guided reaching movements assessed with fMRI-guided TMS
Cereb Cortex
2011
, vol. 
21
 (pg. 
1602
-
12
)
Rizzolatti
G
Luppino
G
Matelli
M
The organization of the cortical motor system: new concepts. Electroencephalogr
Clin Neurophysiol
1998
, vol. 
106
 (pg. 
283
-
296
)
Rozzi
S
Calzavara
R
Belmalih
A
Borra
E
Gregoriou
GG
Matelli
M
, et al. 
Cortical connections of the inferior parietal cortical convexity of the macaque monkey
Cereb Cortex
2006
, vol. 
16
 (pg. 
1389
-
417
)
Schmahmann
JD
Pandya
DN
Wang
R
Dai
G
D’Arceuil
HE
de Crespigny
AJ
, et al. 
Association fibre pathways of the brain: parallel observations from diffusion spectrum imaging and autoradiography
Brain
2007
, vol. 
130
 (pg. 
630
-
53
)
Schulz
R
Park
C-H
Boudrias
M-H
Gerloff
C
Hummel
FC
Ward
NS
Assessing the integrity of corticospinal pathways from primary and secondary cortical motor areas after stroke
Stroke
2012
, vol. 
43
 (pg. 
2248
-
51
)
Schulz
R
Zimerman
M
Timmermann
JE
Wessel
MJ
Gerloff
C
Hummel
FC
White matter integrity of motor connections related to training gains in healthy aging
Neurobiol Aging
2014
, vol. 
35
 (pg. 
1404
-
1411
)
Sharma
N
Baron
J-C
Rowe
JB
Motor imagery after stroke: relating outcome to motor network connectivity
Ann Neurol
2009
, vol. 
66
 (pg. 
604
-
16
)
Shikata
E
Hamzei
F
Glauche
V
Knab
R
Dettmers
C
Weiller
C
, et al. 
Surface orientation discrimination activates caudal and anterior intraparietal sulcus in humans: an event-related fMRI study
J Neurophysiol
2001
, vol. 
85
 (pg. 
1309
-
14
)
Tomassini
V
Jbabdi
S
Klein
JC
Behrens
TEJ
Pozzilli
C
Matthews
PM
, et al. 
Diffusion-weighted imaging tractography-based parcellation of the human lateral premotor cortex identifies dorsal and ventral subregions with anatomical and functional specializations
J Neurosci
2007
, vol. 
27
 (pg. 
10259
-
69
)
Tsutsui
K-I
Jiang
M
Sakata
H
Taira
M
Short-term memory and perceptual decision for three-dimensional visual features in the caudal intraparietal sulcus (Area CIP)
J Neurosci
2003
, vol. 
23
 (pg. 
5486
-
95
)
Tunik
E
Frey
SH
Grafton
ST
Virtual lesions of the anterior intraparietal area disrupt goal-dependent on-line adjustments of grasp
Nat Neurosci
2005
, vol. 
8
 (pg. 
505
-
11
)
Vingerhoets
G
Contribution of the posterior parietal cortex in reaching, grasping, and using objects and tools
Front Psychol
2014
, vol. 
5
 pg. 
151
 
Wang
L
Yu
C
Chen
H
Qin
W
He
Y
Fan
F
, et al. 
Dynamic functional reorganization of the motor execution network after stroke
Brain
2010
, vol. 
133
 (pg. 
1224
-
38
)
Witt
ST
Laird
AR
Meyerand
ME
Functional neuroimaging correlates of finger-tapping task variations: an ALE meta-analysis
Neuroimage
2008
, vol. 
42
 (pg. 
343
-
56
)
Yin
D
Song
F
Xu
D
Peterson
BS
Sun
L
Men
W
, et al. 
Patterns in cortical connectivity for determining outcomes in hand function after subcortical stroke
PLoS One
2012
, vol. 
7
 pg. 
e52727
 

Author notes

*These authors contributed equally to this work.

Supplementary data