Contralesional White Matter Alterations in Patients After Hemispherotomy

Introduction

Cerebral lesions may cause two types of structural changes in the brain: degeneration and neuroplastic reorganization (Rossini et al., 2003; Murphy and Corbett, 2009; Dodd et al., 2017). Both degeneration and neuroplastic reorganization can occur in the ipsilesional hemisphere or the contralesional, intact hemisphere (Buetefisch, 2015; Dodd et al., 2017). The concept that neuroplastic reorganization occurs in the contralesional hemisphere is particularly intriguing, as it entails the idea that one hemisphere can compensate for the losses of the other. However, this idea gets complicated by the primarily inhibitory interhemispheric influences produced via transcallosal pathways (Ferbert et al., 1992; Murase et al., 2004). It has been implied that inhibitory influences from the ipsilesional hemisphere hamper successful neuronal reorganization of the contralesional hemisphere and, thereby, prevent compensatory functional reorganization (Manganotti et al., 2008; Chieffo et al., 2013). In patients after hemispherotomy, the influence of the ipsilesional hemisphere on the contralesional hemisphere is completely removed: hemispherotomy is a hemispheric disconnection procedure indicated as a neurosurgical ultima ratio option in the treatment of selected patients with medically therapy-refractory focal epilepsy (Schramm, 2002). It has been discussed to which degree the hemispherotomy is a stimulus for neuroplastic reorganization in addition to the underlying lesion, conceivably causing structural white matter alterations (Shimizu et al., 2002; Choi et al., 2010; Nelles et al., 2015; Küpper et al., 2016). The results of preceding degeneration and neuroplastic reorganization of white matter may be reflected by Diffusion Tensor Imaging (DTI). In DTI, the MR signal is sensitized to the water diffusion rate. Directional diffusion is described as anisotropic, while directionally unrestricted diffusion is described as isotropic. DTI datasets may be used to model the underlying diffusivity characteristics in each voxel. The most common three-dimensional model mapping of anisotropic (and isotropic) diffusion is the diffusion tensor (Basser et al., 1994a,b). The most commonly used scalar DTI measure is fractional anisotropy (FA; Werring et al., 2000; Beaulieu, 2002; Puig et al., 2017). Importantly, FA may inform investigators about the microstructural white matter status. In this study, we used DTI to evaluate contralesional white matter changes after extended unilateral early brain lesions in the absence or presence of the ipsilesional hemisphere by comparing epilepsy patients who had undergone transsylvian functional hemispherotomy (hemispherotomy group) to nonsurgical patients with similar pathologies (non-surgery group). Thereby, the influence of the hemispherotomy on the plastic reorganization of the ipsilesional hemisphere could be estimated. Also, a group of healthy controls was examined. We hypothesized, that neuroplastic reorganization as indicated by higher FA could be observed in both patient groups, but was more pronounced in the hemispherotomy group, and within this group, in patients with earlier disease onset.

Materials and Methods

Subjects

Fifty-one epilepsy patients (32 females, mean age ± SD: 25.69 ± 12.77 years) with early acquired brain lesions (24 left-sided) and 43 age-matched healthy controls (26 females, mean age ± SD: 25.07 ± 11.33 years) with no history of neurologic or psychiatric disorders were included in this study. For the hemispherotomy group, inclusion criteria were the following: (1) treatment at the Department of Epileptology at the University of Bonn Medical Center between 1992 and 2014; (2) transsylvian functional hemispherotomy (Schramm et al., 2001) at the Department of Neurosurgery at the University of Bonn Medical Center; (3) onset of pathology and diagnosed epilepsy before the age of 8 years; (4) damage to the corticospinal system and impaired motor function; and (5) ability to undergo MR-scanning. For the non-surgery group, the inclusion criteria were the following: (1) treatment at the Department of Epileptology at the University of Bonn Medical Center between 1992 and 2014; (2) no neurosurgical procedure in patient history; (3) onset of pathology and diagnosed epilepsy before the age of 8 years; (4) damage to the corticospinal system and impaired motor function; and (5) ability to undergo MR-scanning. Thirty of the 51 patients met the inclusion criteria of the hemispherotomy group. They underwent transsylvian functional hemispherotomy at mean age 14.23 years as a treatment for pharmacoresistant epilepsy and were examined postoperatively (19 females, 16 left-sided, mean age ± SD: 23.03 ± 9.71 years). Twenty five patients of the hemispherotomy group were also analyzed as part of a previously published study (Gaubatz et al., 2020). Twenty-one patients met the inclusion criteria of the non-surgery group (13 females, eight left-sided, mean age ± SD: 29.43 ± 15.69 years). Twenty-seven of the whole group of 51 patients suffered from perinatal strokes or intracranial hemorrhages resulting in porencephaly, 12 patients were diagnosed with Rasmussen- and other encephalitides, nine patients presented neurodevelopmental disorders such as hemimegalencephaly, Sturge–Weber syndrome, polymicrogyria, or schizencephaly, two patients showed early tumor lesions and there was one patient with stroke in early childhood. In a second step, the hemispherotomy group was further divided into two subgroups according to the onset of their pathology. A cut-off age of 6 months was chosen to distinguish between an early and a late-onset subgroup (with all porencephaly patients being part of the early-onset group). The characteristics of all patients are specified in Table 1. The study was approved by the local Institutional Review Board and all participants and/or their legal guardians gave written informed consent.

TABLE 1

Table 1. Overview: subjects.

Image Acquisition

Magnetic resonance imaging data was acquired using a 3T Magnetom Trio (Siemens Healthineers). T1, T2 and DTI data were collected for all patients and healthy controls. As part of a scanner update, a new head coil was implemented in October 2014 leading to minimal changes in sequence parameters: all scans acquired before the scanner update were run with an eight-channel-coil and a scanning routine containing a 3D MPRAGE sequence (resolution = 1.0 × 1.0 × 1.0 mm³, TR = 1,570 ms, TE = 3.42 ms, flip angle = 15°), a 3D T2-weighted sequence (resolution = 1.0 × 1.0 × 1.0 mm³, TR = 3,200 ms, TE = 455 ms, flip angle = 120°), and a single-shot diffusion-weighted sequence (resolution = 1.72 × 1.72 × 1.7 mm³, TR = 12,000 ms, TE = 100 ms, flip angle = 90°) with 60 diffusion-encoding directions and a b-value of 1,000 s/mm² as well as six baseline volumes with a b-value of 0 s/mm². Sequences acquired after the scanner update were run with a 32-channel head-coil equally including a MPRAGE sequence (resolution = 0.8 × 0.8 × 0.8 mm³, TR = 1,660 ms, TE = 2.54 ms, flip angle = 9°), T2-weighted sequence (resolution = 0.8 × 0.8 × 0.8 mm³, TR = 3,200 ms, TE = 401 ms, flip angle = 120°), and a single-shot diffusion-weighted sequence (resolution 1.72 × 1.72 × 1.7 mm³, TR = 9,000 ms, TE = 87 ms, flip angle = 90°) with 60 diffusion-encoding directions and a b-value of 1,000 s/mm² as well as six baseline volumes with a b-value of 0 s/mm².

Preprocessing

Preprocessing was run using FMRIB’s Software Library 5.0 (Jenkinson et al., 2012) and the Tolerable Obsessive Registration and Tensor Optimization Indolent Software Ensemble (TORTOISE; Pierpaoli et al., 2010). In a first step, T1- and T2-weighted volumes were brain extracted (Smith, 2002) and corrected for intensity inhomogeneities (Zhang et al., 2001). As extensive lesions, as present in our study group, tend to bias the normalization process, individual masks covering compromised structures of the lesional hemisphere were manually drawn on the T1-weighted volumes. By inverting those masks, only healthy tissue was considered during the computation of all inter-individual transformations. Thus, T1-weighted volumes were normalized to the Montreal Neurological Institute template (MNI152 nonlinear, 6th generation, 1 mm resolution) using linear and non-linear registration (Jenkinson and Smith, 2001). For DTI volumes, susceptibility-induced geometric distortions were corrected using constrained registration (Bhushan et al., 2012) on the same step as motion and eddy current correction. Diffusion tensors were fitted to each voxel in the corrected DTI data and FA was calculated. Finally, the mean of each patients’ six b0 baseline volumes was calculated as a reference image for a boundary-based registration (BBR) to the respective T1-weighted volume.

TBSS

Here, we opted to analyze our datasets using tract-based spatial statistics (TBSS; Smith et al., 2006). TBSS is part of the FMRIB’s Software Library 5.0 and has been developed to face two of the biggest challenges of conventional voxel-wise analysis in DTI analysis: multiple comparison and alignment. In TBSS, voxel-wise analysis is constrained to the FA skeleton, which represents the centers of all white matter tracts common to the study group (Andersson et al., 2007a,b). This FA skeleton is applied to mitigate any residual misalignment after a common mapping function has been applied. In the standard processing pipeline, individual FA maps are normalized applying a nonlinear transformation. Normalization describes the process of aligning an individual brain to a template brain using a special function. In the process of normalization, every voxel is assigned specific coordinates (Snook et al., 2007). After normalization, these coordinates are supposed to correspond to certain anatomical structures. At the end of successful normalization, all brains align well with each other. In the TBSS standard routine, a mean FA map (in standard space) is generated based on all normalized brains and used to create the white matter skeleton. Next, the skeleton is projected on the individual FA map but is deformed based on a local search for maximal FA values. Lastly, voxel-wise testing may be performed within the skeleton as described above. It is not surprising that the registration of healthy brains will never work perfectly, let alone the registration of lesioned brains. The general advice to create a study-specific template (Smith et al., 2006) would mean the following for our study: Non-linear transformation of each FA map to the “most representative” one of the group, affine-alignment to the MNI152 space and transformation of every image into 1 × 1 × 1 mm MNI152 space by combining the warp-field to the target FA image with the affine transformation from that target to MNI152 space. However, since most brains of the patient group presented extensive lesions, a “most representative” image cannot be determined. Instead, each FA-map was transformed to MNI152 space by combining the above mentioned BBR-based registration of the individual mean b0 volume to the respective T1-weighted image with the warp-field aligning this T1-weighted image to MNI152 space. Volumes of patients with lesions in the left hemisphere were flipped along the x-axis, allowing voxel-wise comparisons within the healthy hemisphere of all patients. For between-group comparisons, the same proportions of healthy controls as in the patient groups, respectively, were flipped along the x-axis to exclude the effects of physiological brain asymmetry. Lesion masks were delineated manually on the individual T1 volumes in native space and controlled by a second, independent rater. Normalized lesion masks were binarized and merged to create a canonical lesion mask. This canonical lesion mask was modified with a 50% threshold and used to restrain the TBSS skeleton in the lesional hemisphere to exclude deformed tissue structures from the analysis (see Figure 1).

FIGURE 1

Figure 1. Mean fractional anisotropy (FA) skeleton and canonical lesion mask. Mean FA skeleton (red to yellow) of all subjects and canonical lesion mask (rainbow) created from all patients shown on the FMRIB58_FA template. Volumes of patients with left-hemispheric lesions are flipped along the x-axis. Coordinates are provided in MNI standard space.

Statistics

Statistical analyses were performed with FSL’s tool for nonparametric permutation inference (Winkler et al., 2014). To test whether the FA values of the hemispherotomy group, the non-surgery group, and the healthy control group would show significant differences, we ran a one-factor, three-level analysis of covariance (ANCOVA). Age at scan and a variable indicating whether subjects had been scanned before or after the scanner update were both included as covariates. Directions of effects shown in the overall F-test were determined with three post hoc unpaired (two-sample) t-tests between all groups. Group effects in all post hoc t-tests were equally adjusted for scanner update and age at scan. The anatomical location of resulting clusters was identified with the JHU ICBM-DTI-81 white-matter labels atlas (Mori et al., 2005). In a second approach, FA values were subjected to a two-sample unpaired t-test between subgroups of hemispherotomy patients with an early (21 patients, 13 females, mean age ± SD: 23.29 ± 9.30 years) and late (nine patients, five females, mean age ± SD: 25.11 ± 11.11 years) onset of pathology. Scanner update and age at surgery were distributed equally in both subgroups and, thus, had not to be adjusted for this model. Lastly, we performed a voxel-wise regression analysis to link the volume of the lesion masks to voxel-wise FA values. All analyses were run using Monte Carlo permutation tests with 10,000 permutations. As we performed voxel-wise testing, as many statistical tests as there are voxels in the white matter skeleton were carried out. To keep the conventional level of significance at 5% and to control for the number of false positives, results were corrected for family-wise error if not indicated otherwise. Threshold-free cluster enhancement was performed to avoid an arbitrary initial cluster-forming threshold for final voxel-wise inference (Smith and Nichols, 2009). Resulting clusters were thickened to aid visualization.

Results

Intergroup FA Differences

The voxel-wise ANCOVA controlling for scanner update and age at scan as covariates yielded significant effects of the group factor (p < 0.05): clusters were widely distributed in the contralesional hemisphere, whereas they only covered frontal areas and cerebral peduncle in the lesional hemisphere. Adjusting for scanner update and age at scan, FA values were subjected to two post hoc t-tests comparing the hemispherotomy group and the healthy control group (A) and the non-surgery group and the healthy control group (B). A: widespread clusters indicating higher FA values in hemispherotomized patients as compared to healthy controls were found primarily in the internal capsule, the external capsule, superior and posterior corona radiata, posterior thalamic radiations, and brainstem regions in the contralesional hemisphere (p < 0.05; see Figure 2A, green). In the ipsilesional hemisphere, clusters indicating comparatively higher FA values in the hemispherotomy group were only found in small parts of the brainstem (which had not been deafferented). Comparatively lower FA values in the hemispherotomy group were found in residual callosal pathways, in perilesional frontal and occipital regions as well as in the ipsilesional medullary pyramid (see Figure 2A, red). B: clusters indicating lower FA values in the non-surgery group as compared to healthy controls resembled those confirming lower FA values of the hemispherotomy group as compared to healthy controls in A (Figure 2B, red). In B however, clusters were more widespread as they covered larger parts of the temporal and occipital lobes. Only very small clusters proving higher FA values in non-surgery patients as compared to healthy controls emerged in the contralesional hemisphere and did not survive correction for multiple comparisons (Figure 2B, green), whereas results of all other contrasts applied to show intergroup FA differences were corrected for family-wise error.

FIGURE 2

Figure 2. Intergroup FA differences. (A) Voxel-wise post hoc tests in FA between hemispherotomy group and healthy controls. (B) Voxel-wise post hoc tests in FA between non-surgery group and healthy controls. If not indicated otherwise, tests were corrected for family-wise error and p < 0.05. z indicates axial coordinate in MNI standard space. Please note that the canonical lesions masks differ between the post hoc tests as they are built out of the individual lesion masks of the respective groups contrasted. Skeletonized results were thickened for visualization purposes.

Intragroup FA Differences Within the Hemispherotomy Group

A voxel-wise t-test (p < 0.05) was used to determine differences between subgroups of hemispherotomy patients with early and late-onset of pathology. Patients with late-onset of pathology (>0.5 years old at onset) exhibited higher FA values in residual ipsilesional brainstem regions (Figure 3, blue), while patients with early onset of pathology (<0.5 years old at onset) presented higher FA values in frontal regions and within the corona radiata of the contralesional hemisphere (Figure 3, pink). Results of the intragroup comparison did not survive correction for multiple comparisons. The opposite contrasts did not yield significant results.

FIGURE 3

Figure 3. Intragroup FA differences. Voxel-wise t-test in FA between hemispherotomy subgroups with early-onset and late-onset pathologies (p < 0.05, uncorrected for multiple comparisons). Copper voxels indicate a canonical lesion mask. z indicates axial coordinate in MNI space. Skeletonized results were thickened for visualization purposes.

Association Between FA Alterations and Lesion Size

A voxel-wise regression analysis aiming to explain FA values by lesion size did not yield statistically significant results (p > 0.45), even when omitting a correction for multiple comparisons (p > 0.10).

Discussion

In this study, we contrast contralesional FA differences of two epilepsy patient groups with juvenile brain lesions after hemispherotomy and without neurosurgery. Our hypothesis was partly confirmed by the comparison between patients after hemispherotomy and healthy controls resulting in widespread supra-threshold voxels which indicated higher FA values and most likely extensive neuroplastic reorganization in the hemispherotomy group. However, the same comparison with the non-surgery group left us with widespread supra-threshold voxels indicating comparatively lower FA values, most likely indicating more pronounced degeneration in nonsurgical patients, which is not what we hypothesized. Hemispherotomy/the removal of ipsilesional influences seems to be a prerequisite for contralesional neuroplastic remodeling.

Contralesional White Matter Alterations With and Without Ipsilesional Influences

These results alone do not allow for definite inferences to be drawn, however, they yield insights when put in the context of previous studies. It is a common notion that changes in white matter microstructure happen within a genetically determined framework (Rüber et al., 2016). This framework changes across the lifespan leaving varying room for the white matter effects of the environment. DTI may detect even subtle white matter alterations, however, it remains challenging to specifically attribute these alterations to the influence of genetic factors, environment, or age (Rüber et al., 2016). It is a common notion that illness, aging and a loss of cognitive functions is associated with decreased FA values, whereas learning and rehabilitation are accompanied by increased FA values (Bengtsson et al., 2005; Landi et al., 2011; Engvig et al., 2012; Lebel et al., 2012; Zheng and Schlaug, 2015). Previous studies have been quite controversial concerning the role of the contralesional hemisphere in functional recovery after extended unilateral brain lesions. It is easy to conceive that the contralesional hemisphere is the only possible compensator left in patients after hemispherotomy. In patients with unilateral brain lesions but without deafferentation, some imaging studies have found no signs of structural white matter affection of the contralesional hemisphere (Lindenberg et al., 2010a; Rüber et al., 2012), whereas others found DTI indices of contralesional structural reorganization partly correlating with the degree of motor recovery (Schaechter et al., 2009; Crofts et al., 2011). It should be noted that beside contralesional reorganization, also contralesional degeneration has been observed, most likely constituting the structural correlate of what has been termed diaschisis (von Monakow, 1914; Seitz et al., 1999): the post-lesional change of function constituted in brain areas remote from the anatomical site of damage.

Apart from our results, much speaks in favor of the underlying lesion being a powerful stimulus for neuroplastic reorganization: most convincingly, language organization is induced in the right hemisphere when the lesion has occurred to the left one (Staudt et al., 2002b). Furthermore, a large group of hemispherotomy patients shows unchanged motor scores directly after surgery (Choi et al., 2010; Küpper et al., 2016), while some groups reported even better postoperative motor skills after hemispherotomy (Shimizu et al., 2002; van der Kolk et al., 2013; Küpper et al., 2016). Also, patients with pre-operative mutism have anecdotally been reported to acquire speech after hemispherotomy (Vargha-Khadem et al., 1997; Bulteau et al., 2017; two non-published cases in our department). It can easily be seen how the neuroplastic reorganization of the contralesional hemisphere is hampered by the transcallosal inhibition. In an attempt to diminish interhemispheric imbalance, the ipsilesional hemisphere conceivably shows abnormal hyperexcitability. A recent transcranial magnetic stimulation (TMS) study suggests that the recovery of stroke patients is achieved by the excitability increase of the contralesional hemisphere (Stinear et al., 2015). Higher excitability of the ipsilesional hemisphere equals stronger inhibition of the contralesional hemisphere. The surgical deafferentation would, thus, result in a disinhibition of the contralesional hemisphere allowing the contralesional hemisphere to undergo substantial neuroplastic reorganization (Shimizu et al., 2002). The hemispherotomy in that sense may allow unmasking of the plastic potential of the contralesional hemisphere by taking full advantage of what has been termed “degeneracy”—the ability of structurally nonidentical biological elements (like the two hemispheres) to exert the same function (Noppeney et al., 2004). See Figure 4 for a schematic.

FIGURE 4

Figure 4. Schematic. The juvenile brain lesion creates an interhemispheric imbalance with prevailing inhibitory influences. These inhibitory influences are removed by the procedure of hemispherotomy.

Degeneracy is a biological principle of paramount importance for post-lesional recovery. Additionally, neurorehabilitation may leverage the understanding of interhemispheric influences in the lesioned brain. A burgeoning body of evidence (Hummel and Cohen, 2006; Lindenberg et al., 2010b; Grefkes and Fink, 2016; Ladenbauer et al., 2017) from studies suggests that non-invasive brain stimulation such as transcranial direct current stimulation (tDCS) or TMS may enhance recovery. By downregulating the excitability of the contralesional hemisphere, by upregulating the excitability of the ipsilesional hemisphere or by combining both recoveries after stroke may be facilitated. Non-invasive brain stimulation has not been applied to hemispherotomy patients, however, our results inspire trials aimed at “preparing” patients for surgery by downregulating the ipsilesional hemisphere and anticipating the postoperative status.

Early vs. Late Onset of Pathology

The obtained effect of age at onset on FA values is in line with the previously mentioned idea that early brain lesions may lead to the reinforcement of uncrossed pyramidal tract of the ipsilesional hemisphere: the pyramidal tract is bilaterally organized at birth, but uncrossed connections are normally pruned around the age of one and a half years (Eyre et al., 2001). In patients with early lesions, these normally transient unilateral connections may be preserved (Staudt et al., 2002a) and the contralesional hemisphere overtakes motor function. The preservation of these connections is most likely reflected by comparatively higher FA values in the contralesional frontal lobe, whereas patients with late lesions had to recruit ipsilesional regions for functional compensation as indicated by comparatively higher FA values in the residual ipsilesional brainstem. As it is too late for the preservation of the (contralesional) uncrossed pyramidal tract, the ipsilesional hemisphere itself has to compensate the lesion.

Limitations

However, one should be wary when interpreting our results concerning the fact that some variables known to influence the white matter status could not be controlled for: patients after hemispherotomy may not necessarily have larger lesions, however, it may be guessed that they have a higher seizure burden (which ultimately lead to the indication for hemispherotomy) and a higher antiepileptic drug load, both thought to lower FA values (Jehi, 2017). In the hemispherotomy group, however, comparatively higher FA values have been found. Another constraint of this study is set by the restrictive limits of DTI in the assessment of white matter microstructure. On a microstructural scale, the relation between DTI indices and their biological substrates remains to be elucidated. It is important to acknowledge that there is no “one-to-one association” between any DTI index and a single white matter component (Johansen-Berg and Rushworth, 2009). Distinct microstructural processes may lead to the same DTI observations since there is a coexistence of different biological substrates influencing each other (Cercignani and Gandini Wheeler-Kingshott, 2018). The findings of this study shall, thus, not directly be interpreted as correlates of distinct microstructural processes, but rather as indirect correlates of loss of function or recovery previously defined by previous studies. It is a limitation of this study, that no function-structure associations were shown. Furthermore, results of the subgroup-comparison were not corrected for multiple comparisons and we cannot overemphasize the fact that they should, thus, be interpreted with great caution. Lastly, the most important limitation of our study lies in the heterogenicity of the study groups. The number of patients with lesions in the left and the right hemisphere was not balanced and also there were several different etiologies present. Regarding the unbalanced lesion sides in the patient group, it should be noted that in this study, patient brains and an equal proportion of the respective control group were flipped along the x-axis as part of the preprocessing. While this step is inevitable for a joint analysis of the contralesional hemisphere, it has been shown that the white matter architecture of the brain shows a left > right asymmetry in FA (Takao et al., 2011). In a DTI twin study (Jahanshad et al., 2010), 20% to 40% of FA brain asymmetry was even found to be genetically determined. The asymmetry of white matter architecture is being accounted for by flipping the equal proportion of the respective control group of healthy subjects, however, what we do not know is whether a left hemisphere changes differently after lesion/hemispherotomy occurs to the right hemisphere as compared to a right hemisphere after lesion/hemispherotomy occurs to the left one. Regarding the different pathologies in the patient groups (perinatal strokes or intracranial hemorrhages resulting in porencephaly = 27 patients, Rasmussen- and other encephalitides = 12 patients, hemimegalencephaly, Sturge–Weber syndrome, polymicrogyria, or schizencephaly = 9 patients, early tumor lesions = 2 patients), it is likely that they will have different effects on white matter degeneration and reorganization. Whereas the lesion size does not seem to be influential on white matter diffusivity characteristics (as indicated by our regression analysis, see “Association Between FA Alterations and Lesion Size” section), the time onset of the lesion along with the lifespan surely is (as seen in the comparison between patients with early and late-onset pathology). While the concept of an all-decisive “critical phase” in early childhood with a later loss of plastic potential has been softened, partly due to imaging studies with older patients (Jäncke, 2009), the core of the so-called Kennard principle (Kennard, 1936) with a higher plastic potential in early childhood or adolescence still applies (Dennis, 2010). However, one has to keep in mind that hemispherotomy is a neurosurgical procedure extremely rarely performed. Between 100 and 160 yearly hemispherotomies/hemispherectomies have been estimated to be performed in the United States between 2000 and 2009 (Lin et al., 2015). It is, thus, extremely difficult to ascertain a homogenous group of hemispherotomy patients large enough to perform analyses with sufficient statistical power. It is also for this reason that we opted to use TBSS for the analysis of our datasets, as TBSS gains gain statistical power from the dimensionality reduction of projecting FA values on the white matter skeleton (Bach et al., 2014). While this fact explains our limitations, it does not overcome them. A future longitudinal study investigating how white matter microstructure changes before and after hemispherotomy would be one way to overcome some of the limitations of the current study.

Conclusion

Lower FA values unexpectedly found in the white matter infrastructure of the contralesional hemisphere of epilepsy patients with juvenile brain lesion when compared to healthy controls most likely reflect secondary degeneration and may constitute the structural correlate of diaschisis. This degeneration may be aggravated and neuroplastic reorganization in this hemisphere may be hindered by the inhibitory influence of the ipsilesional hemisphere produced via transcallosal pathways. The contrasting result of higher FA values in the contralesional hemisphere in a group of patients with similar pathologies who had undergone hemispherotomy as compared to healthy controls may be interpreted as a correlate of preceding neuroplastic reorganization of the contralesional hemisphere. This reorganization has most likely been enabled by the lacking inhibitory influence of the ipsilesional hemisphere resulting in the disinhibition of the contralesional one. Our study probes models of interhemispheric balance by applying it to two patient groups with early brain lesions after and without hemispherotomy. The pattern of degenerative and plastic structural alterations is most likely a result of interhemispheric inhibition and, therefore, substantiates these models. We conclude, that the procedure of hemispherotomy should be seen as a major plastic stimulus for reorganization in the contralesional hemisphere. Complete disconnection of inhibitory influences from the ipsilesional hemisphere releases the full rehabilitative potential of the contralesional hemisphere for functional recovery in patients with extended unilateral lesions.

Data Availability Statement

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.

Ethics Statement

The studies involving human participants were reviewed and approved by Ethics Committee of the Medical Faculty of the University of Bonn. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.

Author Contributions

JG, CP, EH, BW, HV, RS, CE, and TR contributed to the conception and design of the study. JG, CP, and LE organized the database. JG, LE, BD, and TR performed the statistical analysis. JG wrote the first draft of the manuscript. CH and TR wrote sections of the manuscript. EH, BW, HV, RS, CE, and TR gave administrative, technical or material support. EH, BW, HV, RS, CE, and TR supervised the study. All authors contributed to manuscript revision, read and approved the submitted version. JG and CP have contributed equally to this work.

Funding

Funded by the BONFOR research commission of the medical faculty of the University of Bonn (2015–6–08). JG and LE hold a promotion scholarship of the BonnNi graduate school funded by the Else-Kröner-Fresenius-Stiftung. CP holds a promotion scholarship of the BONFOR research commission of the medical faculty of the University of Bonn.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

Part of the results of this study was presented as a talk (“Contralesional white matter alterations with and without ipsilesional influences”) at the annual meeting of the German Neurological Society (DGN) 2018, which was held from October 30th to November 3rd in Berlin. We are grateful for the kind support provided by the Verein zur Förderung der Epilepsieforschung e.V. We thank Ms. Laura Schinabeck (University of Bonn, Bonn, Germany) for technical assistance and Ms. Karen Wörmann and Ms. Nurida Boddenberg (Department of Epileptology, University of Bonn Medical Center, Bonn, Germany) for editing the manuscript for non-intellectual content.

References

Andersson, J. L. R., Jenkinson, M., Smith, S., and Andersson, J. (2007a). Non-linear optimisation FMRIB technial report TR07JA1. Available online at: https://www.fmrib.ox.ac.uk/datasets/techrep/tr07ja1/tr07ja1.pdf.

Google Scholar

Andersson, J. L. R., Jenkinson, M., Smith, S., and Andersson, J. (2007b). Non-linear registration aka spatial normalisation FMRIB technial report TR07JA2. Available online at: https://www.fmrib.ox.ac.uk/datasets/techrep/tr07ja2/tr07ja2.pdf.

Google Scholar

Bach, M., Laun, F. B., Leemans, A., Tax, C. M. W., Biessels, G. J., Stieltjes, B., et al. (2014). Methodological considerations on tract-based spatial statistics (TBSS). NeuroImage 100, 358–369. doi: 10.1016/j.neuroimage.2014.06.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Basser, P. J., Mattiello, J., and Lebihan, D. (1994a). Estimation of the effective self-diffusion tensor from the NMR Spin Echo. J. Magn. Reson. 103, 247–254. doi: 10.1006/jmrb.1994.1037

PubMed Abstract | CrossRef Full Text | Google Scholar

Basser, P. J., Mattiello, J., and LeBihan, D. (1994b). MR diffusion tensor spectroscopy and imaging. Biophys. J. 66, 259–267. doi: 10.1017/cbo9780511544958.007

CrossRef Full Text | Google Scholar

Beaulieu, C. (2002). The basis of anisotropic water diffusion in the nervous system-a technical review. NMR Biomed. 15, 435–455. doi: 10.1002/nbm.782

PubMed Abstract | CrossRef Full Text | Google Scholar

Bengtsson, S. L., Nagy, Z., Skare, S., Forsman, L., Forssberg, H., and Ullén, F. (2005). Extensive piano practicing has regionally specific effects on white matter development. Nat. Neurosci. 8, 1148–1150. doi: 10.1038/nn1516

PubMed Abstract | CrossRef Full Text | Google Scholar

Bhushan, C., Haldar, J. P., Joshi, A. A., and Leahy, R. M. (2012). Correcting susceptibility-induced distortion in diffusion-weighted (MRI) using constrained nonrigid registration. Signal Inf. Process. Assoc. Annu. Summit. Conf. APSIPA Asia Pac. 2012, 1–9.

PubMed Abstract | Google Scholar

Buetefisch, C. M. (2015). Role of the contralesional hemisphere in post-stroke recovery of upper extremity motor function. Front. Neurol. 6:214. doi: 10.3389/fneur.2015.00214

PubMed Abstract | CrossRef Full Text | Google Scholar

Bulteau, C., Jambaqué, I., Chiron, C., Rodrigo, S., Dorfmüller, G., Dulac, O., et al. (2017). Language plasticity after hemispherotomy of the dominant hemisphere in 3 patients: implication of non-linguistic networks. Epilepsy Behav. 69, 86–94. doi: 10.1016/j.yebeh.2017.01.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Cercignani, M., and Gandini Wheeler-Kingshott, C. (2018). From micro- to macro-structures in multiple sclerosis: what is the added value of diffusion imaging. NMR Biomed. 32:e3888. doi: 10.1002/nbm.3888

PubMed Abstract | CrossRef Full Text | Google Scholar

Chieffo, R., Inuggi, A., Straffi, L., Coppi, E., Gonzalez-Rosa, J., Spagnolo, F., et al. (2013). Mapping early changes of cortical motor output after subcortical stroke: a transcranial magnetic stimulation study. Brain Stimul. 6, 322–329. doi: 10.1016/j.brs.2012.06.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Choi, J. T., Vining, E. P. G., Mori, S., and Bastian, A. J. (2010). Sensorimotor function and sensorimotor tracts after hemispherectomy. Neuropsychologia 48, 1192–1199. doi: 10.1016/j.neuropsychologia.2009.12.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Crofts, J. J., Higham, D. J., Bosnell, R., Jbabdi, S., Matthews, P. M., Behrens, T. E. J., et al. (2011). Network analysis detects changes in the contralesional hemisphere following stroke. NeuroImage 54, 161–169. doi: 10.1016/j.neuroimage.2010.08.032

PubMed Abstract | CrossRef Full Text | Google Scholar

Dennis, M. (2010). Margaret Kennard (1899–1975): not a “principle” of brain plasticity but a founding mother of developmental neuropsychology. Cortex 46, 1043–1059. doi: 10.1016/j.cortex.2009.10.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Dodd, K. C., Nair, V. A., and Prabhakaran, V. (2017). Role of the contralesional vs. ipsilesional hemisphere in stroke recovery. Front. Hum. Neurosci. 11:469. doi: 10.3389/fnhum.2017.00469

PubMed Abstract | CrossRef Full Text | Google Scholar

Engvig, A., Fjell, A. M., Westlye, L. T., Moberget, T., Sundseth, Ø., Larsen, V. A., et al. (2012). Memory training impacts short-term changes in aging white matter: a longitudinal diffusion tensor imaging study. Hum. Brain Mapp. 33, 2390–2406. doi: 10.1002/hbm.21370

PubMed Abstract | CrossRef Full Text | Google Scholar

Eyre, J. A., Taylor, J. P., Villagra, F., Smith, M., and Miller, S. (2001). Evidence of activity-dependent withdrawal of corticospinal projections during human development. Neurology 57, 1543–1554. doi: 10.1212/WNL.57.9.1543

PubMed Abstract | CrossRef Full Text | Google Scholar

Ferbert, A., Priori, A., Rothwell, J. C., Day, B. L., Colebatch, J. G., and Marsden, C. D. (1992). Interhemispheric inhibition of the human motor cortex. J. Physiol. 453, 525–546. doi: 10.1113/jphysiol.1992.sp019243

PubMed Abstract | CrossRef Full Text | Google Scholar

Gaubatz, J., Ernst, L., Prillwitz, C. C., David, B., Lüchters, G., Schramm, J., et al. (2020). Pyramidal tract and alternate motor fibers complementarily mediate motor compensation in patients after hemispherotomy. Sci. Rep. 10:1010. doi: 10.1038/s41598-020-57504-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Grefkes, C., and Fink, G. R. (2016). Noninvasive brain stimulation after stroke. Curr. Opin. Neurol. 29, 714–720. doi: 10.1097/WCO.0000000000000395

CrossRef Full Text | Google Scholar

Hummel, F. C., and Cohen, L. G. (2006). Non-invasive brain stimulation: a new strategy to improve neurorehabilitation after stroke? Lancet Neurol. 5, 708–712. doi: 10.1016/S1474-4422(06)70525-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Jahanshad, N., Lee, A. D., Barysheva, M., McMahon, K. L., de Zubicaray, G. I., Martin, N. G., et al. (2010). Genetic influences on brain asymmetry: a DTI study of 374 twins and siblings. NeuroImage 52, 455–469. doi: 10.1016/j.neuroimage.2010.04.236

PubMed Abstract | CrossRef Full Text | Google Scholar

Jäncke, L. (2009). The plastic human brain. Restor. Neurol. Neurosci. 27, 521–538. doi: 10.3233/RNN-2009-0519

PubMed Abstract | CrossRef Full Text | Google Scholar

Jehi, L. (2017). Searching for autoimmune epilepsy: why, where and when? Epilepsy Curr. 17, 363–364. doi: 10.5698/1535-7597.17.6.363

PubMed Abstract | CrossRef Full Text | Google Scholar

Jenkinson, M., Beckmann, C., Behrens, T. E., Woolrich, M. W., and Smith, S. M. (2012). FSL. NeuroImage 62, 782–790. doi: 10.1016/j.neuroimage.2011.09.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Jenkinson, M., and Smith, S. (2001). A global optimisation method for robust affine registration of brain images. Med. Image Anal. 5, 143–156. doi: 10.1016/s1361-8415(01)00036-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Johansen-Berg, H., and Rushworth, M. F. S. (2009). Using diffusion imaging to study human connectional anatomy. Annu. Rev. Neurosci. 32, 75–94. doi: 10.1146/annurev.neuro.051508.135735

PubMed Abstract | CrossRef Full Text | Google Scholar

Kennard, M. A. (1936). Age and other factors in motor recovery from precentral lesions in monkeys. Am. J. Physiol. Content. 115, 138–146. doi: 10.1152/ajplegacy.1936.115.1.138

CrossRef Full Text | Google Scholar

Küpper, H., Kudernatsch, M., Pieper, T., Groeschel, S., Tournier, J. D., Raffelt, D., et al. (2016). Predicting hand function after hemidisconnection. Brain 139, 2456–2468. doi: 10.1093/brain/aww170

PubMed Abstract | CrossRef Full Text | Google Scholar

Ladenbauer, J., Ladenbauer, J., Külzow, N., de Boor, R., Avramova, E., Grittner, U., et al. (2017). Promoting sleep oscillations and their functional coupling by transcranial stimulation enhances memory consolidation in mild cognitive impairment. J. Neurosci. 37, 7111–7124. doi: 10.1523/JNEUROSCI.0260-17.2017

PubMed Abstract | CrossRef Full Text | Google Scholar

Landi, S. M., Baguear, F., and Della-Maggiore, V. (2011). One week of motor adaptation induces structural changes in primary motor mortex that predict long-term memory one year later. J. Neurosci. 31, 11808–11813. doi: 10.1523/JNEUROSCI.2253-11.2011

PubMed Abstract | CrossRef Full Text | Google Scholar

Lebel, C., Gee, M., Camicioli, R., Wieler, M., Martin, W., and Beaulieu, C. (2012). Diffusion tensor imaging of white matter tract evolution over the lifespan. NeuroImage 60, 340–352. doi: 10.1016/j.neuroimage.2011.11.094

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, Y., Harris, D. A., Curry, D. J., and Lam, S. (2015). Trends in outcomes, complications and hospitalization costs for hemispherectomy in the United States for the years 2000–2009. Epilepsia 56, 139–146. doi: 10.1111/epi.12869

PubMed Abstract | CrossRef Full Text | Google Scholar

Lindenberg, R., Renga, V., Zhu, L. L., Betzler, F., Alsop, D., and Schlaug, G. (2010a). Structural integrity of corticospinal motor fibers predicts motor impairment in chronic stroke. Neurology 74, 280–287. doi: 10.1212/wnl.0b013e3181ccc6d9

CrossRef Full Text | Google Scholar

Lindenberg, R., Renga, V., Zhu, L. L., Nair, D., and Schlaug, G. (2010b). Bihemispheric brain stimulation facilitates motor recovery in chronic stroke patients. Neurology 75, 2176–2184. doi: 10.1212/WNL.0b013e318202013a

PubMed Abstract | CrossRef Full Text | Google Scholar

Manganotti, P., Acler, M., Zanette, G. P., Smania, N., and Fiaschi, A. (2008). Motor cortical disinhibition during early and late recovery after stroke. Neurorehabil. Neural Repair 22, 396–403. doi: 10.1177/1545968307313505

PubMed Abstract | CrossRef Full Text | Google Scholar

Mori, S., Wakana, S., Nagae-Poetscher, L. M., and Van Zijl, P. C. M. (2005). MRI Atlas of Human White Matter. Elsevier Science 2005, eBook ISBN: 9780080456164.

Google Scholar

Murase, N., Duque, J., Mazzocchio, R., and Cohen, L. G. (2004). Influence of interhemispheric interactions on motor function in chronic stroke. Ann. Neurol. 55, 400–409. doi: 10.1002/ana.10848

PubMed Abstract | CrossRef Full Text | Google Scholar

Murphy, T. H., and Corbett, D. (2009). Plasticity during stroke recovery: from synapse to behaviour. Nat. Rev. Neurosci. 10, 861–872. doi: 10.1038/nrn2735

PubMed Abstract | CrossRef Full Text | Google Scholar

Nelles, M., Urbach, H., Sassen, R., Schoene-Bake, J. C., Tschampa, H., Traeber, F., et al. (2015). Functional hemispherectomy: postoperative motor state and correlation to preoperative DTI. Neuroradiology 57, 1093–1102. doi: 10.1007/s00234-015-1564-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Noppeney, U., Friston, K. J., and Price, C. J. (2004). Degenerate neuronal systems sustaining cognitive functions. J. Anat. 205, 433–442. doi: 10.1111/j.0021-8782.2004.00343.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Pierpaoli, C., Walker, L., and Irfanoglu, M. O. (2010). “TORTOISE: an integrated software package for processing of diffusion MRI data,” in ISMRM 18th Annual Meeting, Stockholm.

Google Scholar

Puig, J., Blasco, G., Schlaug, G., Stinear, C. M., Daunis-i-Estadella, P., Biarnes, C., et al. (2017). Diffusion tensor imaging as a prognostic biomarker for motor recovery and rehabilitation after stroke. Neuroradiology 59, 343–351. doi: 10.1007/s00234-017-1816-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Rossini, P. M., Calautti, C., Pauri, F., and Baron, J. C. (2003). Post-stroke plastic reorganisation in the adult brain. Lancet Neurol. 2, 493–502. doi: 10.1016/s1474-4422(03)00485-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Rüber, T., Elger, C. E., and Weber, B. (2016). “Diffusion tensor imaging and tractography,” in Neuroeconomics an Introduction, eds M. Reuter and C. Montag (Berlin: Springer), 411–424.

Google Scholar

Rüber, T., Schlaug, G., and Lindenberg, R. (2012). Compensatory role of the cortico-rubro-spinal tract in motor recovery after stroke. Neurology 79, 515–522. doi: 10.1212/WNL.0b013e31826356e8

PubMed Abstract | CrossRef Full Text | Google Scholar

Schaechter, J. D., Fricker, Z. P., Perdue, K. L., Helmer, K. G., Vangel, M. G., Greve, D. N., et al. (2009). Microstructural status of ipsilesional and contralesional corticospinal tract correlates with motor skill in chronic stroke patients. Hum. Brain Mapp. 30, 3461–3474. doi: 10.1002/hbm.20770

PubMed Abstract | CrossRef Full Text | Google Scholar

Schramm, J. (2002). Hemispherectomy techniques. Neurosurg. Clin. N. Am. 13, 113–134. doi: 10.1016/s1042-3680(02)80012-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Schramm, J., Kral, T., and Clusmann, H. (2001). Transsylvian keyhole functional hemispherectomy. Neurosurgery 49, 891–901. doi: 10.1227/00006123-200110000-00021

PubMed Abstract | CrossRef Full Text | Google Scholar

Seitz, R. J., Azari, N. P., Knorr, U., Binkofski, F., Herzog, H., and Freund, H. J. (1999). The role of diaschisis in stroke recovery. Stroke 30, 1844–1850. doi: 10.1161/01.str.30.9.1844

PubMed Abstract | CrossRef Full Text | Google Scholar

Shimizu, T., Hosaki, A., Hino, T., Sato, M., Komori, T., Hirai, S., et al. (2002). Motor cortical disinhibition in the unaffected hemisphere after unilateral cortical stroke. Brain 125, 1896–1907. doi: 10.1093/brain/awf183

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, S. M. (2002). Fast robust automated brain extraction. Hum. Brain Mapp. 17, 143–155. doi: 10.1002/hbm.10062

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, S. M., Jenkinson, M., Johansen-Berg, H., Rueckert, D., Nichols, T. E., Mackay, C. E., et al. (2006). Tract-based spatial statistics: voxelwise analysis of multi-subject diffusion data. NeuroImage 31, 1487–1505. doi: 10.1016/j.neuroimage.2006.02.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, S. M., and Nichols, T. E. (2009). Threshold-free cluster enhancement: addressing problems of smoothing, threshold dependence and localisation in cluster inference. NeuroImage 44, 83–98. doi: 10.1016/j.neuroimage.2008.03.061

PubMed Abstract | CrossRef Full Text | Google Scholar

Snook, L., Plewes, C., and Beaulieu, C. (2007). Voxel based versus region of interest analysis in diffusion tensor imaging of neurodevelopment. NeuroImage 34, 243–252. doi: 10.1016/j.neuroimage.2006.07.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Staudt, M., Grodd, W., Gerloff, C., Erb, M., Stitz, J., and Krägeloh-Mann, I. (2002a). Two types of ipsilateral reorganization in congenital hemiparesis: a TMS and fMRI study. Brain 125, 2222–2237. doi: 10.1093/brain/awf227

PubMed Abstract | CrossRef Full Text | Google Scholar

Staudt, M., Lidzba, K., Grodd, W., Wildgruber, D., Erb, M., and Krägeloh-Mann, I. (2002b). Right-hemispheric organization of language following early left-sided brain lesions: functional MRI topography. NeuroImage 16, 954–967. doi: 10.1006/nimg.2002.1108

PubMed Abstract | CrossRef Full Text | Google Scholar

Stinear, C. M., Petoe, M. A., and Byblow, W. D. (2015). Primary motor cortex excitability during recovery after stroke: implications for neuromodulation. Brain Stimul. 8, 1183–1190. doi: 10.1016/j.brs.2015.06.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Takao, H., Abe, O., Yamasue, H., Aoki, S., Sasaki, H., Kasai, K., et al. (2011). Gray and white matter asymmetries in healthy individuals aged 21–29 years: a voxel-based morphometry and diffusion tensor imaging study. Hum. Brain Mapp. 32, 1762–1773. doi: 10.1002/hbm.21145

PubMed Abstract | CrossRef Full Text | Google Scholar

van der Kolk, N. M., Boshuisen, K., van Empelen, R., Koudijs, S. M., Staudt, M., van Rijen, P. C., et al. (2013). Etiology-specific differences in motor function after hemispherectomy. Epilepsy Res. 103, 221–230. doi: 10.1016/j.eplepsyres.2012.08.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Vargha-Khadem, F., Carr, L. J., Isaacs, E., Brett, E., Adams, C., and Mishkin, M. (1997). Onset of speech after left hemispherectomy in a nine-year-old boy. Brain 120, 159–182. doi: 10.1093/brain/120.1.159

PubMed Abstract | CrossRef Full Text | Google Scholar

von Monakow, C. (1914). Lokalisation im Gehirn und Funktionelle Störungen Induziert Durch Kortikale Läsionen. Wiesbaden, Germany: Bergmann JF.

Google Scholar

Werring, D. J., Toosy, A. T., Clark, C. A., Parker, G. J. M., Barker, G. J., Miller, D. H., et al. (2000). Diffusion tensor imaging can detect and quantify corticospinal tract degeneration after stroke. J. Neurol. Neurosurg. Psychiatry 69, 269–272. doi: 10.1136/jnnp.69.2.269

PubMed Abstract | CrossRef Full Text | Google Scholar

Winkler, A. M., Ridgway, G. R., Webster, M. A., Smith, S. M., and Nichols, T. E. (2014). Permutation inference for the general linear model. NeuroImage 92, 381–397. doi: 10.1016/j.neuroimage.2014.01.060

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Brady, M., and Smith, S. (2001). Segmentation of brain MR images through a hidden Markov random field model and the expectation-maximization algorithm. IEEE Trans. Med. Imaging 20, 45–57. doi: 10.1109/42.906424

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, X., and Schlaug, G. (2015). 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:229. doi: 10.3389/fnhum.2015.00229

PubMed Abstract | CrossRef Full Text | Google Scholar