Abstract
Background:
Epilepsy is a widespread neurologic disorder and almost one-third of patients suffer from drug-resistant epilepsy (DRE). Neuromodulation targeting the centromediannucleus of the thalamus (CM) has been showing promising results for patients with generalized DRE who are not surgical candidates. Recently, the effect of CM- deep brain stimulation (DBS) in DRE patients was investigated in the Electrical Stimulation of Thalamus for Epilepsy of Lennox–Gastaut phenotype (ESTEL) trial, a monocentric randomized-controlled study. The same authors described a ‘cold-spot’ and a ‘sweet-spot’, which are defined as the volume of stimulation in the thalamus yielding the least and the best clinical response, respectively. However, it remains unclear which structural connections may contribute to the anti-seizure effect of the stimulation.
Objective:
We investigated the differences in structural connectivity among CM, the sweet-spot and the cold-spot. Furthermore, we tried to validate our results in a cohort of DRE patients who underwent CM-DBS or CM-RNS (responsive neurostimulation). We hypothesized that the sweet-spot would share similar structural connectivity with responder patients.
Methods:
By using the software FMRIB Software Library (FSL), probabilistic tractography was performed on 100 subjects from the Human Connectome Project to calculate the probability of connectivity of the whole CM, the sweet-spot and the cold-spot to 45 cortical and subcortical areas. Results among the three seeds were compared with multivariate analysis of variance (MANOVA). Similarly, the structural connectivity of volumes of tissue activated (VTAs) from eight DRE patients was investigated. Patients were divided into responders and non-responders based on the degree of reduction in seizure frequency, and the mean probabilities of connectivity were similarly compared between the two groups.
Results:
The sweet-spot demonstrated a significantly higher probability of connectivity (p < 0.001) with the precentral gyrus, superior frontal gyrus, and the cerebellum than the whole CM and the cold-spot. Responder patients displayed a higher probability of connectivity with both ipsilateral (p = 0.011) and contralateral cerebellum (p = 0.04) than the non-responders.
Conclusion:
Cerebellar connections seem to contribute to the beneficial effects of CM-neuromodulation in patients with drug-resistant generalized epilepsy.
Introduction
Epilepsy is a common chronic neurological disorder with an estimated pooled incidence rate of 61.4 per 100,000 people per year. 1 Antiepileptic drugs are the first line of therapy, but up to 30% of patients fail to respond to multiple or combined pharmacological therapies and are diagnosed with drug-resistant epilepsy (DRE). 2 Moreover, approximately half of DRE patients are not suitable candidates for resective surgery, which is known to have a high seizure freedom rate. 3 For this cohort of patients, neuromodulation targeting the thalamus has been investigated as a promising surgical treatment. 4 In recent years, deep brain stimulation (DBS) of the anterior nucleus of the thalamus was approved in Europe and the United States as adjunctive therapy for the treatment of patients suffering from focal seizures with or without secondary generalization and refractory to three or more antiepileptic medications.5,6 Conversely, all the neuromodulatory techniques for the treatment of generalized epilepsy remain off-label.7 –9 In this context, DBS targeting the centromedian nucleus (CM) of the thalamus seems to be particularly promising in patients suffering from generalized DRE epilepsy. The first case series of CM-DBS for the treatment of generalized DRE dates back to 1987 by Velasco et al., 10 who reported a significant reduction of clinical seizures after 3 months of cyclic stimulation in five patients suffering from generalized DRE. Following this first report, a recent meta-analysis 11 demonstrated a mean seizure reduction rate of 73.4% after a mean follow-up of 28.7 months. Similarly, responsive neurostimulation (RNS) of CM has been proposed as an alternative approach to control generalized seizures by some USA centers and showed encouraging results.12 –14
The CM is a small thalamic region measuring less than 10 mm in every dimension with an approximately 310 mm3 volume. 15 It belongs to the caudal intralaminar group of thalamic nuclei and is situated within the medial thalamus bordered superiorly by the mediodorsal nucleus, medially by the parafascicular nucleus (Pf), and posteriorly by the pulvinar. 16 The first tracing studies of CM were performed on rats, in which this nucleus is indistinguishable from Pf. 17 Based on this methodology, CM was described as the principal target of projections from the mesencephalic reticular formation and the main source of extra-cortical afferents to the basal ganglia. Therefore, murine CM was considered a pivotal regulator of motor task performance in the cortical–basal ganglia–thalamo–cortical loop. Similar anatomical studies in primates 18 showed that CM was well distinguishable from Pf and demonstrated predominant connections to frontal, prefrontal, and parietal cortices. On the other hand, the primate Pf displayed a more limbic pattern of cortical and subcortical connections. Given its widespread cortical and subcortical connections and the electrophysiologic evidence of diffuse EEG desynchronization produced by active CM stimulation, 19 Velasco’s group proposed that neuromodulation of the centromedian nucleus would be useful for the management of different types of generalized epileptic syndromes.
Recently, Warren et al. 20 performed a volume of tissue activated (VTA)-based probabilistic map analysis on a cohort of 20 Lennox–Gastaut syndrome (LGS) patients from the Electrical Stimulation of Thalamus for Epilepsy of Lennox–Gastaut phenotype (ESTEL) tiral. 21 In this first monocentric randomized-control study, the authors evaluated the clinical efficacy of CM-DBS in patients with LGS, which is a specific subtype of generalized epilepsy. They also generated a sweet-spot and a cold-spot, defined as thalamic subregions where the stimulation yielded above-mean and below-mean seizure frequency reduction, respectively. Interestingly, the sweet-spot covered the anterolateral portion of CM, but it also extended to the posterior subdivision of the ventrolateral nucleus based on the Krauth et al./Morel histological atlas. 22 In the current study, we conducted a probabilistic tractography analysis to elucidate the structural connectivity of CM. Then, we compared these connections with those of the recently described sweet- and cold-spots to find possible differences in connectivity that would explain the anti-seizure effect from CM-neuromodulation. Finally, we tried to validate our results in a cohort of eight patients suffering from chronic generalized DRE, who underwent CM-DBS or CM-RNS. Our hypothesis is that a similar connectivity pattern may be shared by the sweet-spot and the responder patients.
Material and methods
Study subjects
We used ‘minimally pre-processed’ 100 unrelated subjects (58 females, 42 males; mean age = 27) from the Human Connectome Project (HCP) database 23 (http://humanconnectome.org), whose acquisition and pre-processing are described elsewhere. 24 The results from HCP subjects were then validated in a cohort of eight patients suffering from chronic generalized DRE who underwent neuromodulation surgery. Our cohort was composed of two different populations: three adult patients (⩾18 years of age) diagnosed with DRE according to the International League Against Epilepsy (ILAE) guidelines 25 by the attending neurologists at UCLA Ronald Regan Medical Center and five pediatric patients (<18 years of age), whose diagnosis was made by attending pediatric neurologists at UCLA Mattel Children’s Hospital. Of all patients, six underwent CM-RNS and 2 CM-DBS. The choice between DBS and RNS was guided by the expertise and the preference of the attending neurologist caring for the patient. Clinical and imaging data were retrospectively collected from the clinical database of our institution.
Patients’ MRI data acquisition and pre-processing
Before surgery, each patient underwent 3 T magnetic resonance imaging (MRI) on a Siemens Prisma scanner, including T1-weighted anatomical images using the MP-RAGE sequence (TE = 2.98 ms, TR = 2.5 s, matrix = 256 × 256, isotropic 1 mm voxels, and flip angle = 15°) and single-shot spin-echo echo planar imaging for diffusion tensor imaging (TE = 91 ms, TR = 7.6 s, matrix = 128 × 128, voxel size = 1.7 mm × 1.7 mm or 2 mm × 2 mm, slice thickness = 2, b-value = 1000, and 64 directions). MRI scans of pediatric patients were acquired under general anesthesia.
For each patient, T1 images were skull stripped (BET) 26 and nonlinearly co-registered to the MNI152 space by using FNIRT tool in FMRIB Software Library (FSL) (ver. 5.0; http://fsl.fmrib.ox.ac.uk/fsl). 27 Eddy current correction was then used to apply affine registrations to each volume in the diffusion dataset and register it with the initial reference B0 volume. Before performing tractography, all diffusion data were processed using Bayesian Estimation of Diffusion Parameters Obtained Using Sampling Techniques (BEDPOSTX),28,29 which runs Markov Chain Monte Carlo sampling to build up distributions of diffusion parameters at each voxel. Then, a multi-fiber diffusion model was fitted on HCP data by using the FMRIB’s Diffusion (FDT) toolbox in FSL. 29 This model uses Bayesian techniques to estimate a probability distribution function (PDF) on the principal fiber direction at each voxel, accounting for the possibility of crossing fibers within each voxel. Three fibers were modeled per voxel, a multiplicative factor of 1 for the prior on the additional modeled fibers and 1000 iterations before sampling. 28 Finally, DTI data were linearly co-registered to the structural T1 images by using FLIRT tool and nonlinearly co-registered to the standard MNI space as described above.
Surgical technique
All patients or their legal representatives gave written informed consent for the surgical procedure. For adult patients, the CM was initially targeted indirectly based on stereotactic coordinates as reported in the literature30,31 (x = 8–10 mm lateral to the midcommissural point, y = 1 mm anterior to the posterior commissure, z = 0–1 mm superior to the intercommissural plane). The stereotactic target was then adjusted according to the patient’s anatomy on the T1-weighted sequence. For pediatric patients, the CM was manually segmented by an attending neuroradiologist (NS) based on the T1-weighted image. For all patients, the lead trajectory was chosen on the gadolinium-enhanced T1-weighted sequence. The remaining steps of the surgical procedure were similar in the two patient populations, although surgeries were performed by two different attending neurosurgeons (AAB for the adult patients and AF for the pediatric patients). Under general anesthesia, a Leksell stereotactic frame (Elekta, Stockholm, Sweden) was applied to the patient’s head and a thin-cut stereotactic CT was obtained. In a dedicated workstation (Elements, BrainLab, Kapellenstrat, Germany) the CT scan was co-registered with MRI sequences containing the stereotactic plan. SenSight directional leads (Sensight Model B33005, Medtronic Inc., Minneapolis, MN, USA) and four-contacts depth leads (DL-330-3.5, NeuroPace, Mountain View, CA, USA) were bilaterally implanted for CM-DBS and CM-RNS patients, respectively. For CM-RNS patients, the RNS Neurostimulator (model RNS-300M, NeuroPace, Mountain View, CA, USA) was implanted into the skull during the same surgical procedure. On the other hand, CM-DBS patients were discharged the following day and returned as outpatients after 1 week to implant the extension leads and a Percept PC IPG (Medtronic Inc., Minneapolis, MN, USA). Afterward, stimulation parameters were programed during the following neurologic outpatient clinic.
VTAs reconstruction
For each patient, left VTAs were computed using the stimulation parameters (active contacts and amplitude in mA) as recorded at the last follow-up visit. VTAs were estimated using the finite element method (FEM)-based FieldTrip SimBio pipeline in Lead-DBS (ver. 2.6). 32 As reported by previous studies,33,34 VTAs were modeled assuming a homogenous tissue conductivity of 0.1 S/m and thresholded to binary three-dimensional volumes by applying an electric field cutoff of 0.2 V/mm to approximate the gold standard of axon cable models. 35
Tractography seeds and targets
For HCP subjects, three seeds were chosen for the tractography analysis: a left CM mask from the THOMAS atlas 36 and two 3-mm-radius spherical ROI generated in MNI125 space with the Gaussian Kernel function based on the MNI coordinates of the cold- and sweet-spots as reported by Warren et al. 20 (Figure 1). We chose to analyze only the connections of the left CM to be consistent with the findings from the study by Warren et al. 20 In fact, these authors mirrored all VTAs from their patients to the left and conducted the whole analysis for the left hemisphere. Therefore, Warren et al. reported only the left coordinates for both the sweet- and cold-spot in the original paper. 20 Moreover, we did not expect significant differences in the connectivity between the left and right CM because the beneficial effect of unilateral CM stimulation has been known since the first reports by Velasco et al.’s group 31 and even other authors employed unilateral CM-DBS with good seizure control. 37 All three masks were co-registered to the subject diffusion space. When analyzing the connectivity of the eight DRE patients, previously generated VTAs were used as seeds after co-registration with patients’ diffusion space. For both HCP subjects and DRE patients, FreeSurfer (v6.0; http://surfer.nmr.mgh.harvard.edu) 38 was used for cortical surface reconstruction and volumetric segmentation to automatically generate target areas for the tractography analysis.34 –36 Forty-five target areas were chosen and divided into five groups: infratentorial, basal ganglia, sensorimotor, associative, and limbic. The complete list of target areas for each group is presented in Supplemental Table S1.

Seeds employed for tractography analysis in HCP subjects.
Probabilistic tractography
For HCP subjects, probabilistic tractography was performed to delineate the structural connectivity between the CM mask, and the cold- and sweet-spot with the target areas. Using the PDFs previously calculated and the PROBTRACKX function,24,25 we determined the probability of connectivity between each seed and the target areas. From each voxel in the three seed masks, we used the following parameters: 5000 streamlines generated, a 0.2 curvature threshold, and a loop check termination. All the target masks were used as waypoint masks to discard tracts not passing through the target, termination masks to terminate the pathways as soon as they reach the target, and classification masks to quantify connectivity values between the seed and target mask. 29 To measure the strength of connectivity between regions, we used the probability of connectivity, which is defined as the probability that a particular voxel lies on an existing tract. For each seed, the probability of connectivity to each target area was calculated by dividing the number of streamlines connecting the seed to the target area by the total number of streamlines sent out (5000 streamlines per seed voxel),39,40 according to the following formula:
To compensate for the distance-dependent effect, all the calculated probabilities of connectivity were corrected by the average distance of the seed to the corresponding target area.41,42 Furthermore, connections with null or sub-threshold (p > 0.01) probability of connectivity were regarded as false-positive results 43 and not included in the following analysis. Similarly, we calculated the corrected probability of connectivity of patients’ left VTA with the target areas resulted statistically significantly different in the previous analysis. Then, the patients were divided into responders and non-responders using as response criterion a > 50% seizure frequency reduction as evaluated by the attending neurologists at the last follow-up visit. Finally, the VTAs probability of connectivity with the target areas was compared between responders and non-responders.
Statistical analysis
Descriptive statistics, frequencies, and percentages were used to report demographic characteristics. A Shapiro–Wilk normality test was used to assess normality across selected variables along with a close inspection of the data plotted on histograms and Q–Q graphs. Continuous variables were reported as mean plus standard deviation. Categorical variables were reported as absolute numbers and percentages. In HCP subjects, MANOVA was employed to perform multivariate comparisons for the mean probability of connectivity among the three seeds (CM, sweet-spot, cold-spot). In this context, the Pillai–Bartlett trace (V) multivariate test was used as a test for overall statistical significance among all variables and ANOVA univariate F-statistic as a follow-up test in case of multivariate significance. Finally, two post hoc tests were used to uncover specific differences in the statistically significant connections among the three seeds: Bonferroni if Levene’s test of equality of error variances was not violated and Games–Howell in the case of Levene’s test violation. Continuous variables in the patients’ cohort were non-normally distributed, and the Quade’s non-parametric ANCOVA was used to compare means between responders and non-responders while correcting for patients’ age. For all traditional hypotheses, p values <0.05 were considered statistically significant. Statistical analysis was performed using SPSS (version 28.0, IBM, Armonk, NY, USA).
Results
Based on the probability of connectivity to target areas, 15 reliable connections were identified for CM (Supplemental Table S2; Figure 2). Overall, the CM had higher connectivity with the brainstem, pre- and postcentral gyri, superior frontal gyrus, caudate nucleus, and cerebellum. Interestingly, the CM displayed bilateral connections with the superior frontal gyrus and the cerebellum, although ipsilateral connections were predictably higher. When considering each target group separately, the CM had a higher probability of connectivity with the brainstem for the infratentorial group, the caudate for the basal ganglia group, the amygdala for the limbic group, the precentral gyrus for the sensorimotor group, and the superior frontal gyrus for the associative group.

Mean probability of connectivity of left CM seed with the selected target areas.
Comparison among the three seeds in HCP subjects
The multivariate test was statistically significant [V = 0.791, F(32,566) = 11.565, p < 0.001], and the follow-up univariate analysis confirmed significant differences among the three seeds for the connections with the following structures (Table 1): the brainstem and the contralateral cerebellum; the ipsilateral globus pallidus and the putamen; the ipsilateral precentral and postcentral gyri; the ipsilateral and contralateral superior frontal gyri, the ipsilateral superior parietal gyrus, and the inferior parietal gyri; the amygdala.
Mean probability of connectivity of the three seeds with the selected target areas for HCP subjects.
p < 0.005 was considered statistically significant.
F-statistics and p-value reported from comparison among probabilities of connectivity of the three seeds by using ANOVA.
The post hoc test (Supplemental Table S3; Figure 3) revealed that the CM proper had a statistically significantly higher probability of connectivity with the brainstem than the sweet-spot (p < 0.001) and the cold-spot (p = 0.004) but a lower one with the ipsilateral superior frontal gyrus than the other two seeds (p < 0.001). When directly comparing the sweet-spot with the other two seeds, a higher probability of connectivity was demonstrated with the contralateral cerebellum (p < 0.001 for both comparisons), the ipsilateral precentral gyrus (p < 0.001 for both comparisons), and the contralateral superior frontal gyrus (comparison with CM: p = 0.01; comparison with the cold-spot: p < 0.001), but lower connectivity was found with the globus pallidum (comparison with CM: p = 0.005; comparison with the cold-spot: p < 0.001) and the postcentral gyrus (p < 0.001 for both comparisons). Finally, the cold-spot had a higher probability of connectivity with the putamen than the two other seeds (comparison with CM: p < 0.001; comparison with the sweet-spot: p = 0.005) and with the amygdala (p = 0.027) than the sweet-spot, but cold-spot connectivity with the brainstem was significantly lower (p = 0.036) than the sweet-spot.

Main analysis for the HCP subjects.
Patient cohort
The patient cohort was composed of eight subjects diagnosed with DRE. (Table 2). Six patients (75%) underwent CM-RNS and two (25%) received CM-DBS. Most were female (62.5%), whereas the mean age at epilepsy onset and at surgery was 4.59 ± 4.322 years and 19.8 ± 17.09 years, respectively. Five patients had vagus nerve stimulator (VNS) implanted before undergoing CM-neuromodulation, but none had previous resective surgery since they were not deemed suitable candidates. After a mean post-implantation follow-up of 9.5 ± 3.84 months, the mean percentage reduction of seizure frequency was 60.5 ± 34.8%. However, four patients were considered responders as they achieved >50% seizure frequency reduction at the last follow-up visit with a mean percentage of seizure frequency reduction of 88.5 ± 7.8% against 32.5 ± 11% in non-responders (p = 0.006). Three out of four responders (75%) were LGS patients, whereas all the non-responders were non-LGS.
Clinical characteristics of the patients with chronic generalized epilepsy in the present study.
Patient’s sex (M = male, F = female); bpatient’s age at seizures onset expressed in years; cpatient’s age at implantation expressed in years; depilepsy etiology was deemed to be ‘undefined’ after a thorough clinical workout, including careful anamnesis and physical examination, different neuroimaging studies (MRI and PET), video-EEG and stereo-EEG in selected cases.
JAE, Juvenile Absence Epilepsy; LGS, Lennox–Gastaut syndrome; RNS, responsive neurostimulation.
When comparing the connectivity of VTAs between responders and non-responders with the previously identified target areas (Table 3), responder patients demonstrated a significantly higher probability of connectivity with the ipsilateral (p = 0.011) and contralateral cerebellum (p = 0.04). Furthermore, higher connectivity was found with the precentral gyrus, the ipsilateral and contralateral superior frontal gyri in responder patients and with the amygdala, the globus pallidum, the putamen, the postcentral gyrus, and the superior parietal gyrus in non-responder patients. Finally, no differences resulted in the connections with the brainstem (p = 0.89) and no connections at all (p = 0) with the inferior parietal gyrus in all the patients.
Mean probability of connectivity for patient’s left VTA with the selected target areas.
Response to stimulation is defined as >50% seizure frequency reduction as evaluated by the attending neurologists at the last follow-up visit.
p < 0.005 was considered statistically significant.
Discussion
Our probabilistic tractography study demonstrated that the CM had the strongest structural connectivity with the brainstem, basal ganglia, cerebellum, primary motor and sensory cortices, and the prefrontal cortex, confirming non-human tract tracing studies.16,18 This widespread connectivity pattern supports the involvement of the primate CM in brain functions other than motor regulation (such as attention, arousal, 44 and nociception 45 ). In fact, CM-neuromodulation has been investigated for the treatment of different neuropsychiatric pathologies such as Tourette’s syndrome, 46 neuropathic pain, 47 coma, and other vegetative states. 48 On the other hand, weaker structural connections of CM with limbic structures resulted, except for the amygdala and the hippocampus. When investigating possible differences in structural connectivity among the CM proper and the cold- and sweet-spots as described by Warren et al., 20 we demonstrated that the sweet-spot had a significantly higher probability of connectivity with the ipsilateral precentral gyrus (primary motor cortex), the brainstem, and the contralateral superior frontal gyrus and cerebellum, but lower connectivity with the postcentral gyrus (primary sensory cortex) and the superior and inferior parietal gyri.
We then confirmed these findings in a cohort of 8 DRE patients who underwent CM-neuromodulation, since VTAs from responders had a pattern of connections similar to the sweet-spot. Our results are in accordance with the analysis by Warren et al.,20,49 who showed that probabilistic CM maps associated with increased seizure reduction were more strongly connected with the cerebellum and frontal areas than with parietal cortices. Interestingly, histological tracing studies indicated that most fibers connecting the CM with frontal cortices and the cerebellum were located in its inferolateral part, which is histologically known as the ‘parvocellular’ component (pvCM) and has been indicated as the targeting site associated with the best clinical outcome by previous studies.31,50 Therefore, it is not surprising that the sweet-spot described by Warren et al. 20 was located on the border between pvCM and the posterior subdivision of the ventral lateral nucleus, which also receives afferents from the cerebellum and projects to the motor and premotor cortices. 51 In fact, the current literature agrees in assigning the frontal cortex an important pathophysiological role in generalized epilepsy.52,53 In particular, it has been shown that the inhibition threshold of neurons within the motor cortex is higher in patients with generalized epilepsy than in both healthy subjects and patients with focal epilepsy. 54 Moreover, non-invasive motor cortex recording in epileptic patients who underwent ANT-DBS demonstrated increased intracortical inhibition 55 and encouraging clinical results were reported in closed-loop stimulation targeting the motor cortex in some cases of medically and surgical refractive epilepsy.56,57 On the other hand, structural and functional alterations involving the superior frontal gyrus have been related to ictal activity in various epileptic syndromes.58,59 More recently, an anatomic parcellation of the superior frontal gyrus has been proposed, encompassing three subregions43,60: an anteromedial component connected to the anterior cingulate cortex and other limbic areas; a dorsolateral component connected with the middle and inferior frontal gyri; and a posterior component connected to the motor cortex and the basal ganglia. Based on this classification, the superior frontal gyrus can access many cortical and subcortical sites and might be a relevant hub for seizure initiation or spreading in generalized epilepsy.
Among CM connections, the cerebellar ones seem to have a particularly promising role in seizure control. In fact, we found stronger cerebellar connections with the sweet-spot and this finding was confirmed in responder patients, who displayed a statistically significant higher probability of connectivity with the cerebellum bilaterally than non-responder patients. The cerebellum has widespread connections61,62 and is now known to take part in brain functions other than sensorimotor regulation, such as language, cognition, and emotion. 63 The role of the cerebellum in epilepsy had initially been suggested by anecdotal reports of cerebellar lesions causing seizures reversible with surgery 64 and post-mortem anatomical abnormalities found in the cerebellum of epilepsy patients with varying etiologies.65,66 More recently, neuroimaging studies have reported increased cerebellar activity during and after seizures 67 and identified altered cerebellar functional connectivity as a predictor of drug resistance in epileptic patients. 68 Stemming from these findings, stimulation of the cerebellar cortex was one of the first neuromodulatory approaches used to treat DRE, 69 but it yielded disappointing results.70 –72 A possible explanation is that the complex foliation of the cerebellar cortex and its compartmentalization into functional zones (i.e. the zebrin bands) may cause a nonspecific spread of electrical stimulation. 51 As a result, the attention of researchers shifted toward the deep cerebellar nuclei, which receive converging afferents from cortical Purkinje cells and are directly connected with many brainstem and thalamic nuclei. 73 In fact, recent murine models of epilepsy showed that direct optogenetic activation of deep cerebellar nuclei can effectively stop both primary generalized 74 and hippocampal-onset seizures. 75 The tract connecting the CM with the cerebellar white matter in our analysis runs through the superior cerebellar peduncle (i.e. brachium conjunctivum), which is the output of the cerebellar dentate nucleus. Interestingly, functional connectivity between this nucleus and the primary motor cortex was altered in a recent study on patients with generalized epilepsy with unknown etiology. 76 Taken together, this evidence strengthens the idea that activation of cerebellar nuclei may stop seizure activity, perhaps by desynchronizing a wide array of subcortical and cortical structures via their widespread polysynaptic connections. 77 These conclusions assume particular importance if considering that our patient cohort was composed of both adult and pediatric subjects with different epileptic syndromes not limited to LGS. Consequently, CM connections with the cerebellum might positively impact seizure control irrespective of epilepsy etiology.
On the other hand, the role of the connections with the brainstem and the basal ganglia was not completely elucidated by our analysis. In fact, the CM displayed the highest probability of connectivity with the brainstem, but its projections seem to involve mainly the mesencephalon and the pons. The relevance of the brainstem in seizure spreading has been known since the first electrophysiologic studies in animals, 78 and its importance in the pathophysiology of epilepsy disorders other than LGS (such as generalized absence seizures 79 ) is acknowledged. Although many brainstem nuclei have been involved in epilepsy pathophysiology (such as the subthalamic nucleus, 80 the superior colliculus, 81 and the cholinergic reticular nuclei 82 ), the exact role of each nucleus in seizure generation or spreading is still unclear. Since the sweet-spot had stronger connections with the brainstem than the cold-spot, our work seemed to agree with previous findings20,49,83 in correlating CM connectivity to the brainstem with positive clinical outcomes. However, the probability of connectivity with the brainstem was almost identical between responder and non-responder patients. On the contrary, the results from HCP subjects would suggest a negative role for the connections with the basal ganglia, since the cold-spot in our analysis demonstrated higher connectivity with the globus pallidum and the putamen than the sweet-spot. By performing a VTA-based analysis similar to Warren et al.20,49 on a cohort formed by patients with different generalized epileptic syndromes, Torrez-Diaz et al. 83 demonstrated a negative correlation between seizure reduction after CM-DBS and structural connectivity with the basal ganglia, whereas Warren et al.20,49 found the opposite. These findings strengthened the emerging idea that different generalized epileptic syndromes would benefit from the stimulation of different CM subcomponents, since a diverse connectivity pattern may underlie each seizure pathophysiology. 84 In fact, the same authors from the ESTEL trial formulated a cortically driven theory for LGS pathophysiology,59,85 centered on early generated paroxysmal fast activity in prefrontal cortices and their subsequent propagation via the extrapyramidal corticoreticular pathway to the brainstem and basal ganglia. Therefore, it is reasonable that LGS patients would benefit the most from stimulation of pvCM, which is mainly connected with the frontal and prefrontal areas. 86 Even in our cohort, all LGS patients were responders. However, the differences in basal ganglia connectivity between responders and non-responders in our patients’ cohort were not statistically significant, and the role of basal ganglia connectivity remains unclear.
Limitations
The main limitation of probabilistic tractography analysis is the rate of false-positive tracts, which is common to the algorithms so far employed and particularly relevant when a whole-brain approach is selected. 87 Consequently, we chose a targeted-based approach by using an MRI-validated CM mask as seed. Nonetheless, the targeted tractography strategy is inherently biased by a priori anatomical knowledge of white matter tracts. 88 To overcome this limitation, we selected as target areas a wide array of anatomical structures spanning beyond the current evidence in the literature. Moreover, the DRE patient cohort we analyzed to validate our connectivity results was composed of a relatively small number of subjects, who were also heterogenous based on age and epilepsy etiology, programming settings, and CM targeting strategies. Although cerebellar connections were stronger in both the sweet-spot and responder patients, further studies with a larger cohort are needed to clarify the role of the connections with other important structures, such as the brainstem and basal ganglia. Finally, our work accounted only for structural connections of CM, whereas animal models of CM-DBS implementing DTI analysis with functional connections from fMRI data described a wider array of brain connectivity for this thalamic nucleus. 89 In fact, structural and functional connectivity profiles of the same anatomic area may not completely overlap, and fMRI activation may result from indirect structural connections between different brain regions. 90 Interestingly, a recent study demonstrated that epileptogenic post-stroke lesions were included in a brain network defined by functional connectivity with the cerebellum and the basal ganglia, and that VTAs from ANT-DBS patients were associated with the same functional network. 91 Therefore, further studies employing both structural and connectivity analysis are warranted to corroborate the results of our work.
Conclusion
Our study showed that the CM nucleus had stronger structural connectivity with the brainstem, the basal ganglia, the cerebellum and the primary sensorimotor and premotor areas, but weaker connectivity with structures traditionally classified as limbic. Since the sweet-spot demonstrated a higher probability of connectivity with frontal and prefrontal cortices, our analysis corroborated the beneficial effect of CM-neuromodulation for LGS patients, according to the cortically centered pathophysiology of this particular disease. Nevertheless, our results also suggested a promising role for the cerebellum, which might contribute to controlling generalized seizures by desynchronizing the whole brain via its widespread structural connections. That being said, generalized epilepsy is a complex disease and abnormal electrical activity in different brain regions may underline a similar epileptic semeiology. The fact that stronger CM connections with the cerebellum in the HCP cohort were supported by the results of the same analysis in the patients’ cohort led us to hypothesize that cerebellar activity via CM stimulation might halt seizures in generalized epilepsy. However, it is more realistic that cerebellar connections are one piece of a more complicated puzzle explaining the anti-seizure effect of CM-neuromodulation, and other connections may also exert a function for seizure control. In fact, several studies sustained that a personalized pattern of connections may benefit different patients according to their epilepsy subtypes and the choice of the surgical target for neuromodulation in generalized DRE should be tailored accordingly.69,92 In this context, our study demonstrated that probabilistic tractography can be employed to identify the most significant CM connections, which can be specifically targeted during preoperative surgical planning as previously illustrated by DBS for movement disorders93,94 and psychiatric conditions.95,96 Moreover, understanding the connectivity of surgical targets used in neuromodulation procedures for the DRE treatment assumes further importance if considered in the context of closed-loop stimulation, which has been shown to increase the response to neurostimulation in Parkinson’s disease 97 and depressed 98 patients and has also been proposed for epilepsy. 99 Therefore, our results contribute to laying the groundwork for refining the targeting of the CM nucleus in order to improve the clinical outcome after neuromodulation surgery in patients suffering from generalized DRE.
Supplemental Material
sj-docx-1-tan-10.1177_17562864231202064 – Supplemental material for Structural connections of the centromedian nucleus of thalamus and their relevance for neuromodulation in generalized drug-resistant epilepsy: insight from a tractography study
Supplemental material, sj-docx-1-tan-10.1177_17562864231202064 for Structural connections of the centromedian nucleus of thalamus and their relevance for neuromodulation in generalized drug-resistant epilepsy: insight from a tractography study by Luigi G. Remore, Ziad Rifi, Hiroki Nariai, Dawn S. Eliashiv, Aria Fallah, Benjamin D. Edmonds, Joyce H. Matsumoto, Noriko Salamon, Meskerem Tolossa, Wexin Wei, Marco Locatelli, Evangelia C. Tsolaki and Ausaf A. Bari in Therapeutic Advances in Neurological Disorders
Footnotes
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
