Scientific Reports

Sci Rep. 10: 1010

Pyramidal tract and alternate motor fibers complementarily mediate motor compensation in patients after hemispherotomy

Jennifer Gaubatz¹, Leon Ernst¹, Conrad C. Prillwitz¹, Bastian David¹, Guido Lüchters², Johannes Schramm³, Bernd Weber⁴, Rainer Surges¹, Elke Hattingen⁵, Gottfried Schlaug⁶, http://orcid.org/0000-0002-2531-6701Christian E. Elger¹, http://orcid.org/0000-0002-6180-7671Theodor Rüber

theodor.rueber@ukbonn.de

¹⁷⁸

1. 0000 0000 8786 803Xgrid.15090.3dDepartment of Epileptology, University of Bonn Medical Center, Bonn, 53127 Germany

2. 0000 0001 2240 3300grid.10388.32Center for Development Research, University of Bonn, Bonn, 53113 Germany

3. 0000 0000 8786 803Xgrid.15090.3dDepartment of Neurosurgery, University of Bonn Medical Center, Bonn, 53127 Germany

4. 0000 0000 8786 803Xgrid.15090.3dInstitute for Experimental Epileptology and Cognition Research, University of Bonn Medical Center, Bonn, 53127 Germany

5. 0000 0004 1936 9721grid.7839.5Department of Neuroradiology, Goethe University Frankfurt, Frankfurt am Main, 60528 Germany

6. Stroke Recovery Laboratory, Beth Israel Deaconess Medical Center/Harvard Medical School, Boston, 02215 MA USA

7. 0000 0004 1936 9721grid.7839.5Epilepsy Center Frankfurt Rhine-Main, Department of Neurology, Goethe University Frankfurt, Frankfurt am Main, 60528 Germany

8. 0000 0004 1936 9721grid.7839.5Center for Personalized Translational Epilepsy Research (CePTER), Goethe-University Frankfurt, Frankfurt am Main, 60528 Germany

© The Author(s) 2020

DOI: 10.1038/s41598-020-57504-x

Published online: 23 January 2020

Abstract

Motor function after hemispheric lesions has been associated with the structural integrity of either the pyramidal tract (PT) or alternate motor fibers (aMF). In this study, we aimed to differentially characterize the roles of PT and aMF in motor compensation by relating diffusion-tensor-imaging-derived parameters of white matter microstructure to measures of proximal and distal motor function in patients after hemispherotomy. Twenty-five patients (13 women; mean age: 21.1 years) after hemispherotomy (at mean age: 12.4 years) underwent Diffusion Tensor Imaging and evaluation of motor function using the Fugl-Meyer Assessment and the index finger tapping test. Regression analyses revealed that fractional anisotropy of the PT explained (p = 0.050) distal motor function including finger tapping rate (p = 0.027), whereas fractional anisotropy of aMF originating in the contralesional cortex and crossing to the ipsilesional hemisphere in the pons explained proximal motor function (p = 0.001). Age at surgery was found to be the only clinical variable to explain motor function (p < 0.001). Our results are indicative of complementary roles of the PT and of aMF in motor compensation of hemispherotomy mediating distal and proximal motor compensation of the upper limb, respectively.

Introduction

Neuroimaging has substantially advanced our understanding of the neuronal mechanisms of functional motor recovery after brain lesions, thereby, enhancing prediction of motor recovery¹. Diffusion Tensor Imaging (DTI) and tractography have proven useful for the in vivo delineation and assessment of white matter pathways². Different DTI-studies have associated motor function after hemispheric lesions with the microstructural integrity of the pyramidal tract (PT) or alternate motor fibers (aMF)³–⁶. aMF are believed to constitute the imaging correlate of cortico-rubro-spinal or of cortico-reticulo-spinal pathways. They may be reconstructed by means of tractography as they descend from the precentral gyrus through the posterior limb of the internal capsule and the tegmentum pontis⁷–¹⁰. Our understanding of aMF and their role as a compensatory corticospinal system was only recently investigated in patients and has been derived as a comparative neuroanatomical approach from numerous animal studies¹¹–¹³. Some studies have suggested that aMF could compensate for the damaged PT⁶–⁸,¹⁴, while other studies concluded that the portion of intact fibers of the affected PT determines the degree of motor recovery³,⁴,¹⁵. A third possibility is that both systems contribute to motor recovery. Schulz and colleagues found no interaction between diffusivity parameters indexing the microstructural status of PT and aMF in patients after stroke, and thus concluded, that the manner in which they function is “synergistic, but independent”⁷. Our hypothesis states that PT and aMF may operate synergistically by mediating distal and proximal motor functions in the lesioned brain, respectively. Dexterity, the ability to fractionate movement of individual fingers, is believed to be enabled solely by monosynaptic connections of the PT¹⁶,¹⁷. The polysynaptic cortico-rubro-spinal and cortico-reticulo-spinal pathways have been associated with synergistic and proximal movements¹¹,¹⁸,¹⁹. We used DTI and tractography for the reconstruction and microstructural assessment of PT and bilateral aMF (crossing and unilateral/uncrossed aMF) in patients after hemispherotomy and matched healthy controls. Hemispherotomy is a hemispheric disconnection procedure, indicated as a neurosurgical ultima ratio option in the treatment of certain patients with therapy-refractory epilepsy²⁰. Hemispherotomy patients represent an especially interesting case for the investigation of motor recovery: Their lesion occurs very early in the course of development, when the brain is thought to be most plastic, and ipsilesional influences on post-operative motor recovery may be excluded. With this study, we aim to characterize the roles of PT, unilateral/uncrossed aMF (IIaMF; originating in the contralesional hemisphere) and crossing aMF (xaMF; originating in the contralesional hemisphere and crossing in the pons) in facilitating motor recovery by relating measures of motor function to tract-specific diffusivity parameters, namely fractional anisotropy (FA), and to clinical data of patients after hemispherotomy.

Results

Clinical explanatory variables of motor function

Twenty-five patients who underwent a hemispherotomy using a transsylvian approach²¹ were included in this study. Patients presented with different etiologies: twelve patients suffered from perinatal strokes or perinatal intracranial hemorrhages resulting in porencephaly; six patients showed neurodevelopmental disorders such as hemimegalencephaly, Sturge-Weber syndrome or polymicrogyria; and seven patients were diagnosed with Rasmussen- or other encephalitides between the ages of 3 and 7. Anatomical T1-weighted sequence of three exemplary patients are presented in Fig. 1. Patients with development disorders underwent surgery earliest (at 39.5 months on average), whereas patients with progressive disorders underwent surgery latest (at 212.86 months on average). The characteristics of the patients are provided in Table 1.

Figure 1View larger versionFigure 1. T1-weighted sequence of three exemplary patients after hemispherotomy.

See full tableTable 1. Overview: subjects.

A stepwise regression analysis among clinical explanatory variables revealed age at surgery as the only statistically significant (p < 0.001) explanatory variable for motor function of the upper extremity (ueMF). All other explanatory variables were non-significant [affected hemisphere (p = 0.848), etiology (p = 0.972), gender (p = 0.798), time after surgery (p = 0.724), age at scan (p = 0.158)], and thus removed from the model. Bootstrapping confirmed high reproducibility (age at surgery, p < 0.001) and suggested the inclusion of age at surgery for further regression analyses. A consecutive simple regression analysis also confirmed age at surgery as an explanatory variable of ueMF (p < 0.001), with an R² of 0.42 (see Fig. 2): the younger the subject was at surgery, the more likely it showed good motor function.

Figure 2View larger versionFigure 2. Motor function and age at surgery. Regressing motor function of the affected upper extremity with age at surgery.

Imaging explanatory variables of motor function

Unilateral/contralesional PT, unilateral/contralesional IIaMF, and xaMF originating in the contralesional hemisphere and crossing the midline to the ipsilesional hemisphere in the pons could be properly reconstructed in all subjects and tract-specific FA could be extracted (see Fig. 3 for canonical tracts of patients).

Figure 3View larger versionFigure 3. Overview of reconstructed canonical tracts. Trajectory of reconstructed canonical unilateral pyramidal tract (PT), crossed alternate motor fibers (xaMF) and unilateral/uncrossed alternate motor fibers (IIaMF) as they descend from the internal capsule to the basis pontis (PT) and the tegmentum pontis (aMF). Building of canonical tracts is described in the methods section. The operated hemisphere is depicted in bronze. z indicates axial level in MNI space.

Tract-wise FA showed statistically significant differences between the PT and aMF-tracts and between patients and controls: FA values of the PT were higher than FA values of aMF-tracts and FA-values of patients were higher than those of controls (all p < 0.001; see Fig. 4). No statistically significant FA differences between patient subgroups were found (all p > 0.5).

Figure 4View larger versionFigure 4. Boxplot of tract-wise FA values of patients and controls. Patients show higher FA values than controls and PT shows higher FA values than aMF-tracts (all p < 0.001). No statistically significant FA differences between patient subgroups (vascular | developmental | progressive disorders) were found (all p > 0.5). FA: fractional anisotropy.

As age at surgery was found to be a relevant explanatory variable of motor function, it was included as an additional independent variable in all four multiple regression analyses computed to explain global (ueMF, first regression analysis), proximal (proxMF, second regression analysis), distal motor function (distMF, third regression analysis) and maximal finger tapping rate (FTR as an additional measure of distal motor function, fourth regression analysis) of the upper extremity based on tract-specific FA of PT, of IIaMF, and of xaMF. All independent variables were found to be statistically significant explanatory variables of ueMF. This model with four explanatory variables accounted for 65% of the variance (R²) of ueMF and showed positive effects of PT, xaMF, and negative effects of IIaMF and age at surgery. For proxMF, only xaMF (p = 0.001) and age at surgery (p = 0.004) were significant independent explanatory variables showing positive and negative effects, respectively. For distMF, only PT (p = 0.050) and age at surgery (p = 0.049) were significant independent explanatory variables also showing positive effects and negative effects, respectively. The former model accounted for 60% of variance in proxMF and the latter for 46% of variance in distMF. The only independent explanatory variable of FTR was PT (p = 0.027, R² = 0.32) showing positive effects. With a variance inflation factor (VIF) < 2.5 for all variables in all regressions, we had no reason to assume multicollinearity. All results, VIF for all explanatory variables, and bivariate correlation coefficients are displayed in Table 2.

See full tableTable 2. Results of multiple linear regression analyses after bootstrapping, complemented by measures of multicollinearity and bivariate correlation coefficients.

Discussion

In summary, higher FA in the PT was associated to better distal motor function and faster finger tapping; and higher FA in xaMF explained better proximal motor function. In contrast, higher FA in IIaMF explained lower overall motor functions; and higher age at surgery predicted lower overall, proximal, and distal motor functions.

The hemispherotomy case

When compared to those patients who have had a stroke and are more commonly subject to studies on motor recovery, hemispherotomy patients present two particularly interesting characteristics as study subjects: first, brain lesions leading to hemispherotomy were typically acquired pre-, perinatally or in early childhood, at a time when the motor system has not yet fully matured and the brain is thought to be most plastic. Second, the hemispherotomy itself results in the contralesional hemisphere as the only possible mediator of post-surgical motor compensation. Hemispherotomy patients are particularly intriguing because they face a maximally invasive trigger of neuroplasticity at a time when the brain is thought to be maximally plastic.

Age at surgery

Among time after surgery, age at surgery, gender, affected hemisphere (left | right), and etiology (vascular | developmental | progressive disorders), only age at surgery could explain ueMF, proxMF or distMF after hemispherotomy. The non-significance of gender and affected hemisphere does not require interpretation. Time after surgery most likely is not significant as the mean time between age at surgery and age at scan was 8.63 years and patients are thought to have reached a static phase of recovery by then. However, the non-significance of etiology is somewhat unexpected, considering accruing evidence for better motor function of brain lesions occurring earlier in the life-span (i.e., in-born as compared to acquired pathologies). The reason most likely lies in the broad three-tier classification (vascular | developmental | progressive disorders) and the considerable variance of age of lesion within the classes. Not surprisingly, age at surgery differs between patient subgroups (developmental < vascular < progressive disorders). It is, thus, open to discussion, whether age at surgery represents a surrogate marker for age of lesion (as earlier lesion leads to earlier surgery) or as to what extent hemispherotomy itself is a stimulus for neuroplastic reorganization (presuming that earlier hemispherotomy occurs to a more plastic brain) in addition to the underlying lesion. Regarding the primarily inhibitory interhemispheric influences produced via transcallosal pathways, it has been suggested that the surgical deafferentation may result in a disinhibition of the contralesional hemisphere allowing it to undergo substantial neuroplastic reorganization²².

Interpretation of diffusivity parameters

The interpretation of our DTI results is based on the idea that FA informs the microstructural status of white matter tracts. Diffusivity parameters have been equally used to describe white matter degeneration and to pinpoint beneficial microstructural white matter alterations, which fall in the realm of neuroplasticity. Generally, reduced diffusion anisotropy is interpreted as a sign of lesion-induced or age-related degeneration after the maturational peak²³, whereas increased diffusion anisotropy is thought to mirror compensatory and training-induced plastic white matter alterations or white matter maturation²⁴,²⁵. In our study, patients show higher FA values as compared to controls, most likely indicating the result of preceding neuroplastic reorganization. The observed structure–function relations are indicative of the functional relevance of the diffusivity alterations seen.

PT

FA of the PT explained ueMF, distMF and FTR. FTR, as distMF, is seen as an indicator of dexterous hand function. The most likely structural correlate of reconstructed PT is the uncrossed portion of the contralesional and uncrossed PT²⁶. The PT is bilaterally organized at birth, but uncrossed connections are pruned around the age of one and a half years in normal development²⁷. Early lesions to the motor system lead to the preservation of these normally transient ipsilateral connections²⁸. Several electrophysiological and imaging studies have implicated this pathway as mediator of motor recovery after hemispheric lesions¹⁵,²⁹. The pyramidal system is believed to be an evolutionary prerequisite for dexterous hand function including grasping¹⁶, which is why their role as a mediator of distal arm motor function is particularly plausible.

aMF

FA of xaMF was positively associated with proxMF and ueMF. FA of IIaMF was negatively associated with ueMF. The cortico-rubro-spinal and the cortico-reticulo-spinal systems have both been described as anatomical substrates of aMF. In previous DTI-studies, only IIaMF has been investigated and tract-specific FA was found to be inversely related to motor function⁶,⁷. In this study, both IIaMF and xaMF were investigated. xaMF is of special interest as evidence for post-lesional neuroplastic reorganization of crossing rubro-spinal fibers is manifold: A myriad of animal studies investigating subcortical reinnervation has examined projections originating in the primary motor cortex, which after hemidecortation, pyramidotomy or cortical infarcts can sprout fibers crossing the midline and targeting the ipsilesional (!) red nucleus¹³,³⁰,³¹. As fibers from the magnocellular part of the red nucleus also cross the midline (decussation of Forel), cortico-rubro-spinal fibers from the contralesional hemisphere, sprouting to the ipsilesional red nucleus, cross twice on their trajectory from the cortex to the spine and may, thus, exert control over contralesional/paretic side. Results of previous studies and of the current study together lend support to the idea that xaMF corresponds to crossing cortico-rubral (!) fibers. The crossing of rubro-spinal (!) fibers from the magnocellular part of the red nucleus, however, cannot be visualized as the caudal extension of anatomy covered by DTI-volumes is limited. A more strictly unilateral route is discussed for the bilaterally organized cortico-reticulo-spinal fibers¹⁸,²⁹, hence, possibly represented by IIaMF. Whereas the positive association found between tract-specific FA of xaMF and proxMF/ueMF puts additional weight on our interpretation of the anatomical correlates of aMF, the inverse correlation between tract-specific FA of IIaMF and ueMF is more difficult to conceive. Notably, this negative association has been found in earlier studies⁶,⁷ and has been interpreted as stronger compensatory (but lastly ineffective) effort in more severely affected patients. It should be noted that FA of IIaMF was only found to explain ueMF (and not proxMF or distMF).

Limitations

The scope of the current study is set by the familiar limitations of DTI and of tractography². Challenges lie in the verification of virtual fiber bundles as well as in the interpretation of diffusivity measures in terms of their underlying anatomical substrates. The anatomical substrates of xaMF and IIaMF and their functional role in motor recovery in particular should be further elucidated. It should be noted, that the inference on the underlying anatomical correlate is made solely on the basis of the anatomical trajectory of the reconstructed tracts and on their functional profile as derived from the observed relationship between diffusivity parameters and motor function. A further limitation of this study is due to its cross-sectional design. This design does not allow to rule out the alternative (yet unlikely) hypothesis that structural alterations observed predate the hemispherotomy or the underlying lesion and, thus, represent a marker of successful post-operative or post-lesional rehabilitation.

Conclusion and outlook

Our results suggest complementary and quite possibly synergistic roles of unilateral PT and xaMF in mediating motor functions of the affected limb after large hemispheric lesions; neither unilateral PT¹⁵ nor aMF⁶ alone mediate motor function after hemispheric lesion. Judging by our results, motor recovery is driven by exploiting the functional diversity of the motor network with different systems mediating different functions. PT most likely corresponds to the uncrossed part of the pyramidal tract and mediates distMF. The neuroanatomical substrate of aMF most likely is cortico-rubro-spinal- and cortico-reticulo-spinal fibers mediating proxMF. One should be wary of inferences on neuroanatomical underpinnings of xaMF and IIaMF, however, based on previous animal and patient studies, it is suggested that xaMF corresponds to crossing cortico-rubral fibers. Our interpretation of the results is schematically summarized in Fig. 5.

Figure 5View larger versionFigure 5. Schematic of descending pathways mediating motor recovery after hemispherotomy as suggested by this study. The operated hemisphere is marked red.

Investigating motor function after hemispherotomy may result in predictive models helpful in the preoperative evaluation of post-surgery motor outcome. More importantly, our nascent understanding of motor recovery may inspire the development of rehabilitative therapies. Williams and colleagues have speculated that the “cocorticospinal phase of the red nucleus/the rubrospinal system […] may provide a protracted window opportunity to intervene with the therapeutic treatments for motor impairments following perinatal injuries”³². First evidence for the responsiveness of aMF for transcranial direct current stimulation has been found in patients after stroke³³. Also, antibodies may be used to counteract signals inhibiting axonal growth. After administering antibodies directed against the myelin protein NogoA to rodent models of stroke, enhanced sprouting of corticofugal fibers originating in the contralesional hemisphere and targeting the ipsilesional red nucleus has been observed and was associated with improved motor function³¹. Forthcoming advances in our understanding of post-lesional neuroplastic reorganization will provide means to manipulate naturally occurring processes, quite possibly aiming to optimize the collaboration of the pyramidal system and cortico-rubro- or cortico-reticulo-spinal fibers in order to amplify the endogenous potential for motor recovery.

Methods

Subjects

Twenty-five patients (13 women; mean age ± SD: 21.1 ± 7.6 years) who underwent a hemispherotomy using a transsylvian approach²¹ (at mean age ± SD: 12.4 ± 8.1 years) were included in this study. All patients were treated at the Department of Epileptology at the University of Bonn Medical Center and hemispherotomy was performed at the Department of Neurosurgery at the University of Bonn Medical Center between 1992 and 2012. Inclusion criteria were: (1) treatment with hemispherotomy, and (2) the ability to undergo a two-hour MRI scan, in addition to several hours of language and motor testing over two consecutive days. For analyses, patients were grouped according to their pathologies into “vascular”, “developmental”, and “progressive disorders”. Additionally, twenty-six healthy age- and gender-matched controls (16 women; mean age ± SD: 23.31 ± 7.28 years) without neurological or psychiatric diseases underwent the same scanning protocol. This study was approved by the Institutional Review Board of the medical faculty of the University of Bonn. All methods were performed in accordance with the guidelines and regulations of this ethics board and in accordance with the Declaration of Helsinki. Informed consent was obtained from all participants and/or their legal guardians.

Functional motor assessment

The motor function of each patient’s affected upper extremity (ueMF) was quantified using Fugl-Meyer Assessment³⁴. For differential analyses of proximal and distal motor skills, two new sub-scores were built. Items assessing proxMF such as shoulder retraction, elevation, abduction, external and internal rotation, flexion, as well as elbow extension and flexion were grouped into one subscore (maximal value: 15). Items representing distMF, such as mass extension and flexion, flexion in proximal and distal interphalangeal joints, extension in metacarpophalangeal joints, thumb adduction and opposition as well as cylinder, spherical grip, and all subitems assessing distal motor function as part of arm movements were grouped into a Fugl-Meyer Hand subscore (maximal value: 6). Additionally, an index finger tapping test was performed³⁵. FTR of the contralateral index finger was assessed over 20 seconds, averaged over two trials, and then used for further analysis.

Image acquisition

Using a 3 T scanner (Magnetom Trio®, Siemens Healthineers), a 3D T1-weighted sequence (voxel size = 1.0 × 1.0 × 1.0 mm³, TR = 1570 ms, TE = 3.42 ms, flip angle = 15°), a 3D T2-weighted sequence (voxel size = 1.0 × 1.0 × 1.0 mm³, TR = 3200 ms, TE = 455 ms, flip angle = 120°), and a diffusion tensor imaging single shot, dual echo, spin echo, echo planar imaging sequence (voxel size = 1.72 × 1.72 × 1.7 mm³, TR = 12000 ms, TE = 100 ms, flip angle = 90°) with 60 directions and a b-value of 1000 s/mm² as well as six volumes with a b-value of 0 s/mm² were acquired for all subjects in addition to other MR sequences, which are not the focus of the current study. An eight-channel head coil was used for signal reception.

Preprocessing

Data was preprocessed using FMRIB’s Software Library 5.0³⁶ and the Tolerable Obsessive Registration and Tensor Optimization Indolent Software Ensemble (TORTOISE)³⁷. For T1- and T2-weighted volumes, brain extraction was performed using FMRIB’s Brain Extraction Tool followed by bias-field correction of the data using FMRIB’s Automated Segmentation Tool³⁸. T1-weighted volumes were normalized to the Montreal Neurological Institute template (MNI 152, 1 mm³) using a 12 degrees of freedom affine registration, performed with FMRIB’s Linear Image Registration Tool and a non-linear registration of the T1-weighted volume to MNI space, performed with FMRIB’s Non-linear Image Registration Tool³⁹. In patients, lesion-masks were defined manually to exclude the lesion and only include voxels representing healthy brain tissue for the computation of the normalization warp-field.

For DTI volumes, susceptibility-induced geometric distortions were corrected by means of constrained registration⁴⁰ together with motion and eddy current correction using TORTOISE after non-brain tissue was removed by FMRIB’s Brain Extraction Tool⁴¹. FMRIB’s Diffusion Toolbox⁴²,⁴³ was used to calculate the fractional anisotropy (FA) value and the probability distribution of fiber directions for each voxel. The mean of the b0 images was used as reference image in diffusion space and linearly registered to the respective T1 volume. If not indicated otherwise, default parameters were used.

Tractography

Probabilistic tractography was performed for reconstruction of the PT, IIaMF and xaMF using FMRIB’s probtrackx2⁴⁴. Regions of interest (ROIs) for tractography were manually defined on the individual FA maps of all subjects on axial slices in native space under consistent criteria by the same operator. ROI delineation was visually guided by intensity differences on the FA map as well as by color-coded diffusion directions of the overlaid principal vector map. Three different ROIs were defined on axial slices: in the lower pons at the level of the inferior cerebellar peduncles, in the internal capsule, and in the juxtacortical white matter of the frontal lobe. For the PT, ROIs were delineated in the anterior pons (basis pontis) and for aMF, in the posterior pons (tegmentum pontis) (z in MNI space ≈ 38). The posterior limb of the internal capsule (IC) was defined as additional ROI (z in MNI space ≈ 75). A rim of white matter adjacent to the posterior bank of the precentral gyrus (M1) and extending from the deepest point of the central sulcus to the lateral crest of the precentral gyrus constituted the most superior ROI (z in MNI space ≈ 126). For the reconstruction of IIaMF and PT only unilateral/contralesional ROIs were applied. For the reconstruction of xaMF, tractography was operated with the ipsilesional pontine and the contralesional IC and M1 ROI. Exclusion masks were defined in the corpus callosum. Additionally, pontine ROIs not belonging to the respective tracts were set as exclusion masks. Tractography was run in two directions: using the precentral gyrus ROI as a seed region, the IC ROI as a waypoint mask, the pontine ROIs as waypoint/termination mask and vice-versa. After reconstruction, tracts of both directions were added to a single tract and modified with a three percent threshold. This left us with three tracts (see Fig. 3): Unilateral/contralesional PT, unilateral/contralesional IIaMF, and xaMF originating in the contralesional hemisphere and crossing the midline to the ipsilesional hemisphere in the pons. Finally, tracts were cropped at the level of the internal capsule to avoid the inclusion of voxels belonging to several tracts in regions where tracts overlap⁷,³³. FA values of all voxels belonging to the respective tract between the internal capsule and the pontine ROIs were averaged and taken for further analysis as tract-specific FA values.

The canonical tracts displayed in Fig. 3 were generated by a four-step routine: (1) Voxel intensity values of non-thresholded (!) native tracts were (numerically) normalized by the number of voxels of the seed ROI, (2) tracts were (spatially) normalized to MNI-space using computed transformation matrices, (3) voxel intensity values were additionally (numerically) normalized by the highest voxel intensity value of the respective tract, and, lastly, (4) the tracts were modified using a three percent threshold and displayed in three-dimensional space (Fig. 3).

Statistical analysis

Statistical analyses were computed with Stata/IC 14.2 for Mac (StataCorp). Regression models were run to explain motor function based on patient clinical data and tract-specific FA values. Firstly, the most relevant clinical explanatory variables of ueMF were identified using a stepwise regression analysis among the following possible influential factors: time after surgery, age at surgery, gender (female | male), affected hemisphere (left | right) and etiology (vascular | developmental | progressive disorders). Statistical significance was determined as p ≤ 0.05. Secondly, four multiple regression analyses were computed to explain ueMF, proxMF, distMF and FTR, respectively, based on imaging explanatory variables such as tract-specific FA of PT, IIaMF, and xaMF. Estimators’ robustness was confirmed using bootstrapping, (10,000 replications) initiated by a reproducible randomly set seed. Multiple regressions including tract-specific FA values of more than one tract were checked for multicollinearity with variance inflation factor analysis (VIF). Group- or tract-wise FA-differences were tested using two-tailed unpaired t-tests.

Footnotes

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

These authors contributed equally: Jennifer Gaubatz and Leon Ernst.

Acknowledgements

J.G. and L.E. received support from the BonnNi Promotionskolleg Neuroimmunology of the University of Bonn and the Else-Kröner-Fresenius Stiftung. C.P. was funded by a promotion scholarship of the BONFOR research commission of the medical faculty of the University of Bonn. C.E. reports personal fees from UCB, personal fees from Desitin, personal fees from BIAL, personal fees from Eisai, grants from DFG (Deutsche Forschungsgemeinschaft), grants from BMBF (Bundesministerium für Bildung und Forschung), and grants from Marga and Walter Boll Stiftung. TR is supported by the BONFOR research commission of the medical faculty of the University of Bonn. The authors are grateful for the kind support provided the Verein zur Förderung der Epilepsieforschung e.V. They thank Ms. Laura Schinabeck for technical assistance and Nurida Boddenberg for comments on an earlier version of the manuscript.

Author contributions

Study concept and design: J.G., L.E., J.S., B.W., R.S., E.H., G.S., C.E., T.R. Acquisition, analysis or interpretation of the data: J.G., L.E., C.P., B.D., G.L., B.W., E.H., G.S., C.E., T.R. Drafting of the manuscript: J.G., L.E., T.R. Critical revision of the manuscript for intellectual content: J.G., L.E., C.P., B.D., G.L., J.S., B.W., R.S., E.H., G.S., C.E., T.R. Statistical analysis: J.G., L.E., G.L., T.R. Administrative, technical or material support: J.S., B.W., E.H., G.S., C.E., T.R. Study supervision: B.W., E.H., C.E., T.R. All authors reviewed the manuscript.

Data availability

The data that support the findings of this study are available on request from the corresponding author. The data are not publicly available as they contain information that could compromise the privacy of research participants.

Competing interests

The authors declare no competing interests.

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group

PMC Copyright Notice

The articles available from the PMC site are protected by copyright, even though access is free. Copyright is held by the respective authors or publishers who provide these articles to PMC. Users of PMC are responsible for complying with the terms and conditions defined by the copyright holder.

Users should assume that standard copyright protection applies to articles in PMC, unless an article contains an explicit license statement that gives a user additional reuse or redistribution rights. PMC does not allow automated/bulk downloading of articles that have standard copyright protection.

See the copyright notice on the PMC site, https://www.ncbi.nlm.nih.gov/pmc/about/copyright/, for further details and specific exceptions.

References

Stinear CM, author. Prediction of motor recovery after stroke: advances in biomarkers. Lancet Neurol. 2017;16:826–836. [CrossRef] [PubMed]

Johansen-Berg H, Rushworth MFS, authors. Using diffusion imaging to study human connectional anatomy. Annu. Rev. Neurosci. 2009;32:75–94. [CrossRef] [PubMed]

Schaechter JD, et al., authors. Microstructural status of ipsilesional and contralesional corticospinal tract correlates with motor skill in chronic stroke patients. Hum. Brain Mapp. 2009;30:3461–3474. [CrossRef] [PubMed]

Nelles M, et al., authors. Functional hemispherectomy: postoperative motor state and correlation to preoperative DTI. Neuroradiology. 2015;57:1093–1102. [CrossRef] [PubMed]

Koch P, Schulz R, Hummel FC, authors. Structural connectivity analyses in motor recovery research after stroke. Ann. Clin. Transl. Neurol. 2016;3:233–244. [CrossRef] [PubMed]

Ruber T, Schlaug G, Lindenberg R, authors. Compensatory role of the cortico-rubro-spinal tract in motor recovery after stroke. Neurology. 2012;79:515–522. [CrossRef] [PubMed]

Schulz R, et al., authors. Synergistic but independent: The role of corticospinal and alternate motor fibers for residual motor output after stroke. NeuroImage Clin. 2017;15:118–124. [CrossRef] [PubMed]

Lindenberg R, et al., authors. Structural integrity of corticospinal motor fibers predicts motor impairment in chronic stroke. Neurology. 2010;74:280–287. [CrossRef] [PubMed]

Habas, C. & Cabanis, E. A. Cortical projection to the human red nucleus: Complementary results with probabilistic tractography at 3 T. Neuroradiology (2007).

Yeo SS, Chang MC, Kwon YH, Jung YJ, Jang SH, authors. Corticoreticular pathway in the human brain: Diffusion tensor tractography study. Neurosci. Lett. 2011;508:1–4. [PubMed]

Lawrence, D. G. & Kuypers, H. G. J. M. The functional organization of the motor system in the monkey: II. The effects of lesions of the descending brain-stem pathways. Brain (1968).

Belhaj-Saïf A, Cheney PD, authors. Plasticity in the distribution of the red nucleus output to forearm muscles after unilateral lesions of the pyramidal tract. J. Neurophysiol. 2000;83:3147–53. [CrossRef] [PubMed]

Z’Graggen WJ, et al., authors. Compensatory sprouting and impulse rerouting after unilateral pyramidal tract lesion in neonatal rats. J. Neurosci. 2000;20:6561–9. [CrossRef] [PubMed]

Takenobu Y, et al., authors. Motor recovery and microstructural change in rubro-spinal tract in subcortical stroke. YNICL. 2010;4:201–208

Küpper H, et al., authors. Predicting hand function after hemidisconnection. Brain. 2016;139:2456–2468. [CrossRef] [PubMed]

Phillips, C. G. & Porter, R. Corticospinal neurones. Their role in movement. Monogr. Physiol. Soc. v–xii–1–450 (1977).

Schieber MH, author. Motor Control: Basic Units of Cortical Output? Curr. Biol. 2004;14:R353–R354. [CrossRef] [PubMed]

Baker SN, author. The primate reticulospinal tract, hand function and functional recovery. J. Physiol. 2011;589:5603–5612. [CrossRef] [PubMed]

Cahill-Rowley K, Rose J, authors. Etiology of impaired selective motor control: emerging evidence and its implications for research and treatment in cerebral palsy. Dev. Med. Child Neurol. 2013;56:522–528. [CrossRef] [PubMed]

Schramm J, author. Hemispherectomy techniques. Neurosurgery Clinics of North America. 2002;13:113–134. [CrossRef] [PubMed]

Schramm J, Kral T, Clusmann H, authors. Transsylvian keyhole functional hemispherectomy. Neurosurgery. 2001;49:891–901. [PubMed]

Shimizu T, et al., authors. Motor cortical disinhibition in the unaffected hemisphere after unilateral cortical stroke. Brain. 2002;125:1896–1907. [CrossRef] [PubMed]

Werring DJ, et al., authors. Diffusion tensor imaging can detect and quantify corticospinal tract degeneration after stroke. J. Neurol. Neurosurg. Psychiatry. 2000;69:269–272. [CrossRef] [PubMed]

Kochunov P, et al., authors. Fractional anisotropy of water diffusion in cerebral white matter across the lifespan. Neurobiol. Aging. 2012;33:9–20. [CrossRef] [PubMed]

Rüber Theodor, Lindenberg Robert, Schlaug Gottfried, authors. Differential Adaptation of Descending Motor Tracts in Musicians. Cerebral Cortex. 2013;25(6):1490–1498. [CrossRef] [PubMed]

Jankowska E, author. How Can Corticospinal Tract Neurons Contribute to Ipsilateral Movements? A Question With Implications for Recovery of Motor Functions. Neurosci. 2006;12:67–79

Eyre JA, Taylor JP, Villagra F, Smith M, Miller S, authors. Evidence of activity-dependent withdrawal of corticospinal projections during human development. Neurology. 2001;57:1543–1554. [CrossRef] [PubMed]

Staudt M, et al., authors. Right-hemispheric organization of language following early left-sided brain lesions: Functional MRI topography. Neuroimage. 2002;16:954–967. [CrossRef] [PubMed]

Benecke R, Meyer BU, Freund HJ, authors. Reorganisation of descending motor pathways in patients after hemispherectomy and severe hemispheric lesions demonstrated by magnetic brain stimulation. Exp. Brain Res. 1991;83:419–426. [CrossRef] [PubMed]

Jones, T. A. & Adkins, D. L. Motor System Reorganization After Stroke: Stimulating and Training Toward Perfection. Physiology (2015).

Papadopoulos CM, et al., authors. Functional recovery and neuroanatomical plasticity following middle cerebral artery occlusion and IN-1 antibody treatment in the adult rat. Ann. Neurol. 2002;51:433–41. [CrossRef] [PubMed]

Williams PTJA, Kim S, Martin JH, authors. Postnatal Maturation of the Red Nucleus Motor Map Depends on Rubrospinal Connections with Forelimb Motor Pools. J. Neurosci. 2014;34:4432–4441. [CrossRef] [PubMed]

Zheng, X. & Schlaug, G. Structural white matter changes in descending motor tracts correlate with improvements in motor impairment after undergoing a treatment course of tDCS and physical therapy. Front. Hum. Neurosci. 9 (2015).

Fugl-Meyer AR, Jääskö L, Leyman I, Olsson S, Steglind S, authors. The post-stroke hemiplegic patient. 1. a method for evaluation of physical performance. Scand. J. Rehabil. Med. 1975;7:13–31. [PubMed]

Jäncke L, Schlaug G, Steinmetz H, authors. Hand Skill Asymmetry in Professional Musicians. Brain Cogn. 1997;34:424–432. [CrossRef] [PubMed]

Smith SM, et al., authors. Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage. 2004;23:S208–S219. [CrossRef] [PubMed]

Pierpaoli, C., Walker, L. & Irfanoglu, M. O. TORTOISE: an integrated software package for processing of diffusion MRI data. ISMRM18thannual meeting (2010).

Zhang Y, Brady M, Smith S, authors. Segmentation of brain MR images through a hidden Markov random field model and the expectation-maximization algorithm. IEEE Trans. Med. Imaging. 2001;20:45–57. [CrossRef] [PubMed]

Jenkinson M, Bannister P, Brady M, Smith S, authors. Improved Optimization for the Robust and Accurate Linear Registration and Motion Correction of Brain Images Improved Optimization for the Robust and Accurate Linear Registration and Motion Correction of Brain Images. Neuroimage. 2002;17:825–841. [CrossRef] [PubMed]

Bhushan, C., Haldar, J. P., Joshi, A. A. & Leahy, R. M. Correcting susceptibility-induced distortion in diffusion-weighted MRI using constrained nonrigid registration. In Proceedings of The 2012 Asia Pacific Signal and Information Processing Association Annual Summit and Conference 1–9 (2012).

Smith SM, author. Fast robust automated brain extraction. Hum. Brain Mapp. 2002;17:143–155. [CrossRef] [PubMed]

Behrens TEJJ, et al., authors. Characterization and Propagation of Uncertainty in Diffusion-Weighted MR Imaging. Magn. Reson. Med. 2003;50:1077–1088. [CrossRef] [PubMed]

Behrens, T. E. J., Johansen Berg, H., Jbabdi, S., Rushworth, M. F. S. & Woolrich, M. W. Probabilistic diffusion tractography with multiple fibre orientations: What can we gain? (2006).

Behrens TEJ, Berg HJ, Jbabdi S, Rushworth MFS, Woolrich MW, authors. Probabilistic diffusion tractography with multiple fibre orientations: What can we gain? Neuroimage. 2007;34:144–155. [CrossRef] [PubMed]

[Back]

Figure 1.

T1-weighted sequence of three exemplary patients after hemispherotomy.

[Back]

[Back]

Table 1.

Overview: subjects.

Patients
n = 25

Etiology
ID
Lesion Side
Gender
Etiology
Age at surgery (months)
Age at scan (years)
FA
Motor Test Scores

PT
IIaMF
xaMF
ue MF
prox MF
dist MF
FTR

vascular
1
L
M
porecenphaly
130
16
0.70
0.61
0.58
36
12
2
53.5

2
R
M
199
19
0.61
0.51
0.47
19
6
2
32.5

3
R
M
183
17
0.67
0.47
0.51
29
8
4
48.5

4
L
F
162
21
0.57
0.46
0.48
28
10
2
59

5
L
F
223
20
0.59
0.45
0.54
30
13
2
40

6
L
F
118
20
0.69
0.57
0.55
36
10
4
60.5

7
L
F
142
22
0.61
0.59
0.56
29
10
5
55.5

8
R
F
131
21
0.58
0.47
0.47
26
9
2
56.5

9
L
M
126
18
0.69
0.60
0.57
39
10
4
/

10
R
F
303
36
0.70
0.67
0.51
17
6
0
32.5

11
R
F
194
27
0.69
0.52
0.45
18
5
2
51

12
L
F
110
19
0.66
0.50
0.54
31
10
3
58

mean ± SD
7 left
8 females

168.42 ± 55.81
21.33 ± 5.40
0.65 ± 0.05
0.53 ± 0.07
0.52 ± 0.04
28.17 ± 7.20
9.08 ± 2.43
2.67 ± 1.37
49.77 ± 10.27

developmental
13
L
F
HMEG
72
25
0.62
0.46
0.48
40
11
3
55.5

14
L
F
HMEG
4
18
0.59
0.50
0.49
29
7
3
51

15
R
F
HMEG
9
11
0.71
0.52
0.48
50
10
6
62

16
R
M
PLMI
125
23
0.67
0.52
0.50
22
8
2
54

17
R
F
HMEG
8
16
0.62
0.55
0.49
28
13
2
40.5

18
R
M
SWS
19
20
0.68
0.59
0.58
42
11
4
50

mean ± SD
2 left
4 females

39.5 ± 48.93
18.83 ± 5.04
0.65 ± 0.05
0.52 ± 0.04
0.51 ± 0.04
35.17 ± 10.52
10 ± 2.19
3.33 ± 1.51
52.17 ± 7.12

progressive disorders
19
L
M
encephalitis
228
20
0.64
0.58
0.57
24
8
2
51.5

20
L
F
221
20
0.56
0.48
0.48
19
8
1
39.5

21
R
F
370
46
0.62
0.58
0.52
20
6
3
44.5

22
L
F
89
11
0.61
0.51
0.50
21
10
3
23

23
L
M
150
18
0.66
0.48
0.47
24
5
5
60

24
R
F
83
12
0.67
0.50
0.52
51
15
4
58.5

25
R
M
349
32
0.68
0.57
0.52
18
5
2
66

mean ± SD
4 left
4 females

212.86 ± 115.24
22.71 ± 12.37
0.63 ± 0.04
0.53 ± 0.05
0.51 ± 0.03
25.29 ± 11.57
8.14 ± 3.53
2.86 ± 1.35
49 ± 14.68

total
mean ± SD
13 left
16 females

149.92 ± 98.09
21.12 ± 7.68
0.64 ± 0.05
0.53 ± 0.06
0.51 ± 0.04
29.04 ± 9.70
9.04 ± 2.70
2.88 ± 1.36
48.14 ± 14.51

Controls

n = 25

total
mean ± SD

16 females

23.31

7.28

0.59

0.04

0.44

0.03

0.45

0.04

Abbreviations: distMF: distal motor function of the affected upper extremity; F: female; FTR: maximal finger tapping rate; HMEG: Hemimegalencephaly; L: left; M: male; PLMI: Polymicrogyria; proxMF: proximal motor function of the affected upper extremity; R: right; SD: standard deviation; SWS: Sturge-Weber syndrome; ueMF: motor function of the affected upper extremity; /: missing data.

[Back]

[Back]

Figure 2.

Motor function and age at surgery. Regressing motor function of the affected upper extremity with age at surgery.

[Back]

[Back]

Figure 3.

Overview of reconstructed canonical tracts. Trajectory of reconstructed canonical unilateral pyramidal tract (PT), crossed alternate motor fibers (xaMF) and unilateral/uncrossed alternate motor fibers (IIaMF) as they descend from the internal capsule to the basis pontis (PT) and the tegmentum pontis (aMF). Building of canonical tracts is described in the methods section. The operated hemisphere is depicted in bronze. z indicates axial level in MNI space.

[Back]

[Back]

Figure 4.

Boxplot of tract-wise FA values of patients and controls. Patients show higher FA values than controls and PT shows higher FA values than aMF-tracts (all p < 0.001). No statistically significant FA differences between patient subgroups (vascular | developmental | progressive disorders) were found (all p > 0.5). FA: fractional anisotropy.

[Back]

[Back]

Table 2.

Results of multiple linear regression analyses after bootstrapping, complemented by measures of multicollinearity and bivariate correlation coefficients.

Model and Data | Coefficient | SE | Z | p | +95% CI | −95% CI |  | R² | seed
ue MF: main effects |  |  |  |  |  |  |  | 0.651 | 421670
PT | 77.580 | 37.512 | 2.07 | 0.039 | 4.057 | 151.102 | 0.201 |  |
X aMF | 114.941 | 38.725 | 2.97 | 0.003 | 39.040 | 190.841 | 0.283 |  |
II aMF | −66.612 | 33.406 | −1.99 | 0.046 | −132.087 | −1.137 | 0.166 |  |
age at surgery | −0.055 | 0.014 | −3.83 | 0.000 | −0.083 | −0.027 | 0.437 |  |
prox MF: main effects |  |  |  |  |  |  |  | 0.602 | 905414
PT | −6.915 | 12.539 | −0.55 | 0.581 | −31.491 | 17.661 | 0.015 |  |
X aMF | 54.181 | 16.657 | 3.25 | 0.001 | 21.534 | 86.829 | 0.402 |  |
IIaMF | −18.567 | 15.058 | −1.23 | 0.218 | −48.080 | 10.947 | 0.106 |  |
age at surgery | −0.017 | 0.006 | −2.89 | 0.004 | −0.028 | −0.005 | 0.363 |  |
dist MF: main effects |  |  |  |  |  |  |  | 0.464 | 13587
PT | 43.873 | 22.400 | 1.96 | 0.050 | −0.029 | 87.776 | 0.172 |  |
XaMF | 36.054 | 28.083 | 1.28 | 0.199 | −18.988 | 91.095 | 0.091 |  |
IIaMF | −25.565 | 24.176 | −1.06 | 0.290 | −72.949 | 21.819 | 0.070 |  |
age at surgery | −0.023 | 0.012 | −1.97 | 0.049 | −0.046 | 0.000 | 0.258 |  |
FTR: main effects |  |  |  |  |  |  |  | 0.322 | 725623
PT | 123.875 | 55.998 | 2.21 | 0.027 | 14.121 | 233.629 | 0.204 |  |
XaMF | 76.688 | 69.772 | 1.10 | 0.272 | −60.062 | 213.439 | 0.063 |  |
IIaMF | −101.202 | 66.122 | −1.53 | 0.126 | −230.798 | 28.395 | 0.156 |  |
age at surgery | −0.001 | 0.025 | −0.04 | 0.967 | −0.050 | 0.048 | <0.001 |  |
VIF | variables | Pearson’s r (p-value) |  |  |  |  |  |
1.51 | PT | −0.028 (1) | 0.554 (0.024) | 0.321 (0.710) |  |  |  |  |
1.55 | X aMF | 0.021 (1) | 0.581 (0.014) |  |  |  |  |  |
2.21 | II aMF | 0.242 (1) |  |  |  |  |  |  |
1.13 | age at surgery |  |  |  |  |  |  |  |
|  | age at surgery | II aMF | X aMF | PT |  |  |  |

Bold font indicates p ≤ 0.05. Abbreviations: IIaMF: fractional anisotropy of unilateral/uncrossed alternate motor fibers; XaMF: fractional anisotropy of alternate motor fibers originating in the contralesional hemisphere and crossing in the pons; CI: confidence interval; distMF: distal motor function of the affected upper extremity; FTR: maximal finger tapping rate; proxMF: proximal motor function of the affected upper extremity; SE: standard error; ueMF: motor function of the affected upper extremity; PT: fractional anisotropy of the pyramidal tract; VIF: variance inflation factor; η²p: partial eta squared, η²p can be interpreted as the proportion of effect + error variance that is attributable to the effect.

[Back]

[Back]

Figure 5.

Schematic of descending pathways mediating motor recovery after hemispherotomy as suggested by this study. The operated hemisphere is marked red.

[Back]