Thalamocortical Hodology to Personalize Electrical Stimulation for Focal Epilepsy

Thalamocortical Hodology to Personalize Electrical Stimulation for Focal Epilepsy | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Article Thalamocortical Hodology to Personalize Electrical Stimulation for Focal Epilepsy Jorge Gonzalez-Martinez, Arianna Damiani, Sirisha Nouduri, Jonathan Ho, and 7 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-5507011/v1 This work is licensed under a CC BY 4.0 License Status: Published Journal Publication published 20 Oct, 2025 Read the published version in Nature Communications → Version 1 posted You are reading this latest preprint version Abstract Targeted electrical stimulation to specific thalamic regions offers a therapeutic approach for patients with refractory focal and generalized epilepsy who are not candidates for resective surgery. However, clinical outcome varies significantly, in particular for focal epilepsy, influenced by several factors, notably the precise anatomical and functional alignment between cortical regions generating epileptic discharges and the targeted thalamic stimulation sites. Here we hypothesized that targeting thalamic nuclei with precise anatomical and functional connections to epileptic cortical areas (an approach that we refer to as hodological matching) could enhance neuromodulatory effects on focal epileptic discharges. To investigate this, we examined three thalamic subnuclei (pulvinar nucleus, anterior nucleus, and ventral intermediate nucleus/ventral oral posterior nuclei) in a retrospective study involving 32 focal epilepsy patients. Specifically, we first identified hodologically organized thalamocortical fibers connecting these nuclei to individual seizure onset zones (SOZs), combining neuroimaging and electrophysiological techniques. Further, analysis of 216 spontaneous seizures revealed the critical role of matched thalamic nuclei in seizure development and termination. Importantly, electrical stimulation of hodologically-matched thalamic nuclei immediately suppressed intracortical interictal epileptiform discharges, contrasting with ineffective outcomes from stimulation of unmatched targets. Finally, we retrospectively evaluated 7 patients with a chronic hodologically-matched neurostimulation system, which led to a clinically relevant reduction in seizure frequency (median reduction 86.5%), that outstands the current clinical practice of unmatched targets (39%). Our results underscore the potential of hodological thalamic targeting to modulate epileptiform activity in specific cortical regions, highlighting the promise of precision medicine in thalamic neuromodulation for focal refractory epilepsy. Health sciences/Neurology/Neurological disorders/Epilepsy Health sciences/Medical research/Translational research Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Introduction Over 50 million people worldwide suffer from epilepsy, a neurological disorder characterized by recurrent seizures, stemming from synchronized and excessive paroxysmal electrical activity of neurons 1 , 2 . In most of the cases (> 60%), said seizures originate and develop in restricted and focal cortical regions of the brain 3 , condition clinically referred to as focal epilepsy. Conversely, in generalized epilepsy, pathological brain activity stems from widespread areas of the brain, affecting both hemispheres 4 , 5 . Independently of the focal or generalized onset, about 30% of epilepsy patients live with intractable disabling seizures 6 , which cannot be controlled with medical therapy. The current standard of care for medically refractory epilepsy (MRE) relies on the resection of the seizure focus, i.e., the brain region associated with earliest electrophysiological changes during a seizure event 7 . However, if said cortical areas are not amenable for safe removal (or disconnection) because of overlap with eloquent cortices, or extensive seizure foci, or a diagnosis of generalized epilepsy, neuromodulation approaches such as thalamic electrical stimulation might be an alternative 7 – 11 . Recent studies have propelled thalamic stimulation for epilepsy to the forefront of neuroscience, underscoring its emerging role as a pivotal therapeutic strategy in the field 12 . The thalamus is a deep central structure crucial for its multifaceted functions and intricate organization. Indeed, through extensive connections with the cortex and other subcortical and peripheral structures, it is integral to a wide array of functions including motor and sensory processing, consciousness, and sleep 13 – 16 . This diverse functionality is attributed to the thalamus's varied nuclear architecture and its intimate relationship with cortical areas 17 . Importantly, since the early 1950s the thalamus was demonstrated to synchronize widespread cortical activity, hence contributing to the generation of pathological epileptic discharges 18 . Indeed, altered thalamocortical synchrony was extensively observed during focal seizures and associated loss of consciousness 19 . Furthermore, structural connectivity variations was reported between the thalamus and the cortex in patients with generalized epilepsy 20 . Given these cortico-subcortical dynamics and epileptogenic interactions, thalamic electrical stimulation has been proposed since the early 80s, with indications both in focal 9 , 21 and generalized epilepsy 8 . In this regard, thalamic sites characterized by diffuse connectivity as the centromedian nucleus (CM), have been hypothesized to benefit generalized seizures, which originate and develop through widespread neuronal networks. Indeed, in a recent clinical trial, Cukiert and colleagues reported a significant seizure reduction (> 50%) with CM electrical stimulation in 90% of tested patients diagnosed with generalized epilepsy 22 . CM stimulation, instead, has been less effective in patients with focal MRE 23 . For this reason, thalamic nuclei with less diffuse connectivity, and mostly the anterior nucleus (ANT), have been explored as a possible target for neuromodulation in focal epilepsy. Already in the early 80s, Cooper et al. showed a significant clinical control of seizures in 4 out of 6 patients with stimulation of the ANT 24 . Since this seminal work, multiple studies have investigated the efficacy of ANT stimulation to treat focal epilepsy 9 , 21 , 25 , including the pivotal SANTE trial that achieved 75% median seizure reduction at 7 years follow up 26 . Thanks to these works, ANT stimulation for focal epilepsy received FDA approval in 2018. However, clinical outcomes are highly variable depending on the anatomical location of the seizure onset. Specifically, the most well-documented efficacy was observed in focal onset seizures occurring in the temporal and frontal lobe (78% and 86% seizure frequency reduction), whereas effects were limited for other cortical seizure onset zones (SOZs) (39%) 26 . This might be due to the ANT’s unique anatomical connectivity to temporal and frontal areas 27 , 28 . Exploration in other areas of the thalamus further suggests a possible link between thalamocortical anatomical connectivity and neuromodulation efficacy. Indeed, the pulvinar nucleus (PUL) of the thalamus, which has high connectivity towards the temporal and occipito-parietal lobes and present different neural firing patterns during ictal progression and termination in patients with temporal lobe epilepsy 29 – 32 , has been shown to be potentially effective for posterior quadrant epilepsy 33 , 34 . Similarly, the ventrolateral portion of the thalamus (ventral intermediate /ventral oral posterior nuclei, VIM/VOP) 35 and the cerebellum 36 , which present high structural connectivity to motor cortical areas, were tested, with positive outcome, for epilepsies organized in the agranular pre-central neocortical areas (rolandic epilepsy), which is particularly challenging for resective surgeries due to the high risks of motor deficits. However, these studies remain isolated investigations leading to a lack of consensus on which target is optimal for focal epilepsy and clinical indication on how to approach the heterogeneity of SOZs, with limited clinical outcome. Based on prior evidence, here we argue that, for focal epilepsy, thalamic electrical stimulation should match the hodology of the thalamocortical pathway to enhance therapeutic outcome. In other words, targeting electrical stimulation to the thalamic subnucleus with preferential anatomical and functional connectivity to the SOZ (an approach that we refer to as hodological matching ) will improve the efficacy of thalamic stimulation for the treatment of refractory focal epilepsy. To test this hypothesis, we studied the anatomical and functional properties of neural activity of 32 patients with focal MRE. We focused on the study of three thalamic nuclei that exhibit more specific and less diffuse thalamocortical projections as compared to the CM: the ANT, the PUL, and the VIM/VOP. First, we defined a hodological map of the thalamocortical fibers from said thalamic nuclei towards the cortex using imaging and electrophysiology. In this way, we identified the cortical areas that match each nucleus, i.e., receiving major projections. We then characterized the functional involvement of these thalamic nuclei to ictal events occurring both in hodologically-matching and unmatching SOZs. Furthermore, we showed that electrical stimulation of only hodologically-matched thalamic nuclei immediately suppressed pathological discharges i.e., interictal spikes, in the epileptogenic cortex. Finally, we retrospectively verified that such effects led to significant clinical improvement by reducing seizure frequency in 7 patients with a chronic hodologically-matched thalamic implant. Importantly, our hodology-based approach achieved a median reduction in seizure frequency of 86.5%, which dramatically outperforms unmatched thalamic targeting reported by literature (39%). In summary, our approach, i.e., targeting electrical stimulation to match the hodology of thalamocortical connections, holds significant potential to improve the therapeutic efficacy of thalamic electrical stimulation in the of focal refractory epilepsy. This strategy introduces innovative avenues for clinical targeting and therapeutic intervention, while challenging the prevailing paradigm that adopts a “one-nucleus-fits-all” SOZs. Results Study rationale and design Here we aimed at investigating whether specific thalamic subnuclei might provide increased efficacy in the treatment of focal refractory epilepsy based on their hodological matching with the cortical SOZ. To demonstrate our hypothesis, we deployed a large array of electrophysiological and imaging assessments in a unique cohort of 32 focal MRE patients (Fig. 1 a and Supplementary Table 1 for demographic and clinical information for each patient). Specifically, we analyzed intracranial electrophysiological data from 26 patients (S01-S26) undergoing (SEEG) exploration for epilepsy monitoring (Fig. 1 b). Prior to SEEG implantation, we acquired preoperative structural magnetic resonance imaging (MRI) and postoperative computed tomography (CT) for precise SEEG contacts localization performed by a certified neurosurgeon (JGM) in all patients. All the participants had at least one SEEG electrode located in one of the nuclei of interest, namely the PUL, ANT or VIM/VOP nuclei, and a diagnosis of focal epilepsy. Overall, 37 thalamic nuclei among those of interest were sampled in this cohort (Fig. 1 b, Supplementary Table 2 ). We categorized each patient’s SOZ as frontal (meaning originating from frontal lobe regions such as the orbitofrontal cortex), temporal (including mesial and lateral temporal lobe structures), rolandic (including pre- and post-central peri-rolandic regions) or parietal/occipital cortex (namely originating from posterior-quadrant areas) (Fig. 1 c, Supplementary Table 3 ). In a subgroup of patients (n = 8), we acquired and analyzed high-resolution diffusion MRI data to corroborate the anatomical organization of thalamocortical fibers (Fig. 1 f). Such connectivity pattern was also confirmed through electrophysiology with thalamocortical evoked potentials (n = 6). To study the interaction of the thalamus and the epileptic cortex, we then analyzed a total of 216 spontaneous seizures with concurrent sampling of the thalamic nuclei and cortical SOZ. Hence, each recorded seizure was categorized according to the SOZ anatomical location and the sampled thalamic nucleus during the event (Fig. 1 c). We then performed electrophysiological testing in a subset of patients. Specifically, we used high-frequency thalamic stimulation (n = 14) to test immediate effects of hodologically-matched thalamic stimulation on interictal epileptiform activity (Fig. 1 f). Finally, we assessed long-term efficacy of our hodology based neuromodulation approach in 7 patients (S16, S27-32) with a chronic stimulation implant (Fig. 1 d) and up to 3 years of follow-up (Fig. 1 d-e). Importantly, we also reported a single case of generalized epilepsy with simultaneous sampling of PUL, ANT, VIM/VOP as well as CM: this last site showed widespread structural and functional connectivity ( Extended Data Fig. 1 ) , replicating well known evidence from the literature 8 , 11 , 22 , 43 and supporting its applicability as a target in generalized, rather than focal, epilepsy. Overall, we designed and performed a comprehensive study to investigate the mechanisms and efficacy of hodologically-matched thalamic stimulation in focal refractory epilepsy. A specific hodological map of the thalamocortical pathway in the epileptic brain While previous studies proposed a structural connectivity map of the thalamocortical projections 17 , 44 , 45 , the majority of these studies were performed in healthy brains and were limited to imaging techniques, such as fiber tracking. Here, we confirmed the presence of a well-defined anatomical map of the thalamocortical projections of PUL, ANT and VIM/VOP in the epileptic brain with two complementary methodologies: neuroimaging and electrophysiology. First, we performed high-definition fiber tracking (HDFT) of high-resolution diffusion MRI data to confirm previously reported anatomical organization of cortical connectivity patterns of the PUL, ANT and VIM/VOP. We reconstructed all the likely axonal pathways between these nuclei and the cortical regions presenting seizures in our cohort (namely frontal, rolandic, parietal, occipital and temporal) (Fig. 2 a). We quantified the relative strength of these connections by calculating the volume of thalamocortical projections from each nucleus to each cortical region normalized by the total volume of fibers (Fig. 2 b). This analysis revealed a clear thalamocortical anatomical organization. Indeed, the PUL showed preferential connectivity towards the parietal and occipital regions, and the temporal lobe (mean volume value, frontal 0.11, rolandic: 0.04, parietal : 0.38, occipital : 0.27, temporal : 0.18); the ANT, instead, presented axonal projections going preferentially towards the frontal and temporal lobes ( frontal: 0.39 , rolandic: 0.05, parietal: 0.06, occipital: 0.16, temporal: 0.31 ); finally the VIM/VOP nuclei projected more remarkably to the rolandic cortex (frontal: 0.29, rolandic: 0.56 , parietal: 0.03, occipital: 0.06, temporal: 0.03). This constricted structure of thalamocortical fibers from PUL, ANT and VIM/VOP contrasts with the broad anatomical connectivity of the CM (frontal: 0.27, rolandic: 0.17, parietal: 0.19, occipital: 0.19, temporal: 0.15), hence corroborating its indication as a more suitable target in generalized epilepsy ( Extended Data Fig. 1 a-b). We then confirmed this anatomical organization by analyzing thalamocortical evoked potentials (EPs) in the different cortical areas. For this, we applied low frequency (1 Hz) electrical stimulation to the different thalamic nuclei at 1–3 mA while recording cortical neural response evoked by each pulse in the different cortical lobes (Fig. 2 d, g, j). We then analyzed stimulation triggered averages of EPs. Following thalamic stimulation, we observed robust cortical EPs in anatomically connected (i.e., matching), but not in unmatching, cortical regions (Fig. 2 e, h, k). Indeed, analysis of EPs peak to peak amplitudes showed significantly larger responses in those areas that matched the previously identified projections, further confirming a precise anatomical organization of distinct thalamic nuclei (Fig. 2 c, f, i, l). In summary, we confirmed with multiple techniques a specific anatomical organization of the thalamocortical projections of PUL, ANT and VIM/VOP in the studied population. We then define hodologically-matching (parietal, occipital and temporal for PUL, frontal and temporal for ANT, and rolandic for VIM/VOP) and unmatching (frontal and rolandic for PUL, rolandic, parietal and occipital for ANT, frontal, temporal, parietal and occipital for VIM/VOP) regions. Increased thalamocortical coupling during seizures in hodologically-matching SOZs Given this specific structural organization of the thalamocortical fibers, we hypothesized a preferential functional involvement in the course of spontaneous seizures of the thalamic subnucleus that anatomically matches the SOZ (Fig. 3 d). To test this hypothesis, we analyzed a total of 216 spontaneous (i.e., not induced by cortical electrical stimulation) ictal events. The anatomical location of the SOZ as well as the relevant time markers (pre-ictal, beginning and end) of each seizure were determined by a certified clinical team. For each seizure we computed the well validated non-linear h2 correlation coefficient 29 , 46 – 48 to identify time-resolved synchrony between each thalamic nucleus and each cortical SOZ. Specifically, we computed mean h2 coefficient during the ictal event (from seizure onset to termination) and the mean h2 coefficient during a 10 second baseline epoch spaced 2 minutes apart from the seizure’s onset 29 . We then compared mean h2 coefficient between the two phases separately for anatomically connected and non-connected thalamic nuclei to quantify the degree of specificity of thalamic involvement during spontaneous seizures. As hypothesized, we found that while anatomically matched nuclei showed a significant increase of synchrony between the thalamus and the SOZ, indicating functionally specific thalamic involvement during ictal events ( Figure a ), minimal or no correlation was found in unmatched nuclei (Fig. 3 b-c). This was consistent for all thalamic nuclei of interest and subjects (Fig. 3 e, mean h2 was 0.13 vs 0.20 in matched with p 0.05 in unmatched, for baseline and seizure respectively). Importantly, we also observed that thalamic nuclei had minimal or no correlation with non-SOZ regions ( Extended Data Fig. 3 a-b), demonstrating the specificity of the correlation with the anatomically-matched SOZ. Additionally, for the patient diagnosed with generalized epilepsy we found that PUL, ANT, and VIM/VOP nuclei had a lower correlation with the epileptogenic zone (0.16, 0.25, 0.21 mean h2 across 24 seizures and 13 channels for PUL, ANT, and VIM/VOP, respectively) as compared to the CM (0.3, Extended Data Fig. 1 c-d). This confirms the more diffused profile of the CM connectivity (both anatomically and functionally), entailing its suitability for generalized, rather than focal epilepsy treatment. Overall, these functional outcomes robustly validated the thalamocortical anatomical organization observed with fiber tracking and thalamocortical EPs, further supporting the proposed definition of hodologically-matching and unmatching nuclei. To further investigate possible mechanisms of thalamocortical interactions during focal seizures, we explored the coupling of the thalamus and SOZ in terms of seizure phases and frequency bands for matched thalamic nuclei. In particular, we aimed at verifying that the mechanism of thalamic involvement during seizures would not differ across the three nuclei of interest, since a different mechanism could imply a different design of electrical stimulation therapies. First, we asked whether maximal coupling of matched nuclei was specific to distinct seizure phases. To this aim, we computed the mean h2 correlation coefficient in three epochs: seizure initiation (including the 10 seconds preceding the first appearance of a tonic discharge in the SOZ contacts), middle part of the seizure (defined as the time interval separating seizure initiation and termination), and seizure termination (which includes the last 10 seconds of discharge). We found that the highest synchrony between the thalamus and SOZ occurred at seizure termination (Fig. 3 f). Importantly, this behavior was present for all three thalamic nuclei when hodologically (i.e., anatomically and functionally) matched to the SOZ ( Extended Data Fig. 3 c). These findings support and validate the notion that synchronization loops between the thalamus and cortex might contribute to the termination of seizures, as suggested in previous studies in animals 29 , 49 and humans 50 , 51 . Additionally, we sought to determine if the observed synchrony presented a specific spectral profile. Hence, we compared the mean h2 coefficient computed in different frequency bands of the intracerebral recordings (delta, theta, alpha, beta, low gamma, high gamma, ripple, fast ripple). No spectral specificity was revealed: indeed, for all the frequency bands except the fast ripple band, the correlation between hodologically-matched thalamic nuclei and SOZ always increased significantly during the ictal event, while it showed no variation for unmatched thalamic nuclei ( Extended Data Fig. 3 d). Hence, we asked whether specific frequency bands would better distinguish hodologically-matched and unmatched nuclei potentially suggesting biomarkers for the refinement of the stimulation design. For this, we compared the increase in ictal correlation respect to baseline between hodologically-matched and unmatched nuclei in different frequency bands. Interestingly, we found that theta, alpha and beta frequencies had a higher difference between the h2 correlation change during seizures between hodologically-matched and unmatched nuclei (Fig. 3 g). These results show that the theta-beta frequencies might be more relevant to differentiate thalamic nuclei based on their connectivity to the SOZ. All in all, these findings demonstrate that thalamic nuclei exhibit enhanced neural synchrony with hodologically-matched cortical SOZs, particularly at the period of seizure termination, and primarily driven by theta, alpha, and beta oscillations. Thalamocortical synchrony revers directionality during spontaneous seizures While the h2 index revealed an increase in correlation between the hodologically-matched thalamic nuclei and the SOZ throughout the seizure course, here we sought to determine the directionality of this coupling. Importantly, previous studies reported a change in the directional coupling of thalamocortical interactions in temporal lobe epilepsy throughout the seizure course 29 , 52 . Specifically, thalamic involvement was shown to occur early in a seizure, although it did not typically lead at the onset 29 . Additionally, animal studies have shown that the thalamus participates during seizures and plays a crucial role in their development, but it is less involved at seizure onset 49 . In order to explore whether a similar behavior was prevalent for all thalamic structures here investigated, we further expanded the analysis of hodologically-matched thalamocortical synchrony during spontaneous seizures by investigating the directionality of the coupling. To this purpose, we used Granger Causality (GC) 39 between the thalamic nuclei of PUL, ANT and VOP/VIM and the hodologically-matching SOZs. This methodology could indeed provide insights in the understanding of propagation patterns of epileptic activity helping the future design of stimulation patterns aimed at decreasing epileptiform activity. Since our previous results suggest a crucial role in the theta-alpha-beta oscillations, for this analysis we included only band-passed filtered (4–30 Hz) seizures that showed a clear increase in h2 coefficient in this frequency range and lasted at least 10 seconds. This resulted in a total of 98 seizures. For each seizure, we extracted two 10 s epochs of interest at seizure initiation and seizure termination and computed the GC matrix for the SOZ contacts and the hodologically-matching thalamic nucleus contacts (Fig. 4 a). Importantly, GC of unmatched nuclei was not explored because of the lack of correlation with SOZ. To test whether the role of different thalamic nuclei was changing over the course of the seizure, we compared the GC coefficient from the thalamus towards the SOZ (th◊SOZ) and from the SOZ towards the thalamus (SOZ◊th). At seizure initiation, in the majority of the seizures (70%), the GC coefficient of SOZ◊th was higher than the GC coefficients th◊SOZ (Fig. 4 b for all subjects together, Extended Data Fig. 4 for each individual) suggesting a leading role of the cortex. In contrast, more than half of the events (60%) presented a higher th◊SOZ coefficient at the termination of ictal events hence indicating a change in the leading anatomical structure. This could be due to 1) SOZ◊th decreasing while th◊SOZ remaining stable (or decreasing less) throughout the course of the seizure, 2) th◊SOZ increasing and SOZ◊th remaining stable (or increasing less). To test for this, we computed the relative change of the GC coefficient between termination and initiation and compared it between SOZ◊th and th◊SOZ. Interestingly, we noticed an increase in both the thalamocortical (th◊SOZ) and corticothalamic (SOZ◊th) coupling throughout the duration of the event (Fig. 4 c). However, the relative change was higher for the th◊SOZ consistently across all the subjects, further demonstrating a leading role of the thalamus at seizure termination (73.41 vs 345.6%, for SOZ◊th and th◊SOZ, respectively) (Fig. 4 c-d). Importantly, we verified that the same characteristics of thalamocortical interactions were not specific to a particular nucleus. As a matter of fact, in all studied nuclei, the absolute GC value in the initiation epoch and termination epoch showed a reverse trend, with a primary leading role of the cortex at first (85% for PUL, 59% for ANT, 60% for VIM/VOP), followed by a predominantly leading role of the thalamus at termination (63% for PUL, 71% for ANT, 52% for VIM/VOP) (Fig. 4 e). Additionally, PUL, ANT and VIM/VOP showed a significantly higher relative change for th◊SOZ throughout the seizure (Fig. 4 f). In summary, here we investigated the underlying mechanisms of thalamocortical and corticothalamic coupling involving the ANT, PUL and VIM/VOP with hodologically-matched SOZ during spontaneous seizures. Our findings suggest a potential trend where cortical regions exert a dominant influence at seizure initiation, while the thalamus may assume a more prominent role as the ictal events progresses. This pattern was observed across all three nuclei, suggesting that the development of nuclei-specific stimulation protocols -particularly regarding the timing of stimulation during the ictal event- may not be necessary to address the varying underlying mechanisms of thalamic action. Matched thalamic stimulation immediately suppresses epileptiform discharges Given the peculiar hodology that we previously found both in the structural connectivity of the PUL, ANT and VIM/VOP, as well as in the functional coupling observed during ictal events, we hypothesized that targeting electrical stimulation to the thalamic nucleus that hodologically-matches the cortical SOZ might be more effective in reducing epileptiform activity than unmatched nuclei. To test this hypothesis, we delivered bipolar thalamic stimulation through the SEEG electrodes in 14 patients. For clinical reasons, the stimulation testing was always performed after the patient had shown spontaneous seizures. We evaluated the effect of hodologically-matched and unmatched stimulation on interictal epileptiform discharges (IEDs), namely interictal spikes. We used an automatic validated algorithm 41 , 42 to extract IEDs from cortical recordings in the SOZ (Fig. 5 a), previously subdivided in 5 seconds epochs and denoised (see Methods ). We aimed at quantifying the difference between baseline (including the epochs in the 2 minutes preceding the first stimulation applied) and stimulation epochs (see Methods ). For this, we performed two distinct analysis: a) in n = 8 patients with continuous evident pathological spiking (Fig. 5 a), we computed the rate of IEDs per minute in baseline (Stim OFF) and stimulation (Stim ON) epochs; b) in n = 6 patients with more than half of the baseline epochs without pathological spikes, we computed the probability of suppressed IEDs in stimulation epochs and compared it to the probability of suppressed IEDs in baseline epochs, hence obtaining a percentage change in IED suppression probability ( Supplementary Table 2 ). In the first analysis, a positive effect of the stimulation (meaning reducing pathological IEDs) should result in a smaller IEDs rate during stimulation epochs than during baseline (stimulation off) epochs; whereas in the latter, a positive number indicated that the stimulation epochs are more likely to present no IEDs than the baseline epochs (without stimulation). As hypothesized, we found that electrical stimulation of the PUL, ANT and VIM/VOP significantly suppressed cortical IEDs when the SOZ was in a cortical region that hodologically-matched the thalamic nucleus stimulated (Fig. 5 b-d, g). Overall, our analysis revealed that hodologically-matched but not unmatched thalamic stimulation significantly reduced the IED rate in all tested patients (for matched stimulation, IED variation to baseline ranged from − 29.12% to -73% with p < 0.01 and p < 0.001, while for unmatched stimulation, it ranged from − 4.2% to -25%) (Fig. 5 g). Here, we highlight the IED rate values for a representative patient for each nucleus of interest: we observed a reduction of the IED rate of 132 vs 78 spikes/minute for the PUL (median rate at stimulation off vs stimulation on in S01, p < 0.001) (Fig. 5 b), 144 vs 60 spikes/minute for the ANT (in S24, p < 0.001) (Fig. 5 c), and 132 vs 48 spikes/minute in the VIM/VOP (in S13, p < 0.001) (Fig. 5 d). Similarly, when evaluating the percentage change in the probability of IED suppression, we found that hodologically-matched thalamic stimulation epochs were more likely to present no IEDs when compared to unmatched thalamic stimulation epochs. Indeed, In 5 out of 6 patients that received stimulation in the nucleus matching their SOZ, we found at least a 20% increase (up to over 200%) in the probability of epochs with no IEDs with stimulation on with respect to stimulation off (Fig. 5 h); whereas the stimulation effect was more inconsistent (-17 to + 7%) when the targeted nucleus was not matching the hodology of the thalamocortical fibers towards the SOZ (Fig. 5 h). This finding further suggests that the more effective thalamic target to suppress IEDs is the one matching the hodology of the SOZ. Overall, we showed that only targeted electrical stimulation of the nuclei hodologically-matching the SOZ is significantly reducing the rate of IEDs and increasing the probability that stimulation epochs present no IEDs (complete suppression). These findings are crucial in guiding approaches to address the heterogeneity of SOZs and in clinically selecting optimal stimulation targets for the electrical neuromodulation of focal epilepsies. Matched thalamic electrical stimulation in chronic epilepsy treatment Finally, to verify whether the observed immediate effects on IEDs could be translated in a long-term clinical improvement in epileptic patients, we evaluated 7 patients (one of which being part of the SEEG cohort) that were chronically implanted for thalamic stimulation in a hodologically-matching thalamic nucleus, as part of their clinical care. Across all patients, 6 (S16, S27-S31) presented focal MRE originating from the posterior quadrant regions and hence received a neurostimulation system in the PUL, while one (S32) was implanted in the ANT to treat a frontal onset MRE ( Supplementary Table 4 ). Overall, we observed a median reduction in seizure frequency of 86.5% (Fig. 6 a, b) at 15.86 ± 7.9 (mean ± SD) months of follow-up. As a representative example, for S16 seizures reduced from 5 per week to 1 per month, and the patient also self-reported a significant reduction in seizure severity, noting no missed workdays since the implant of the stimulation device, compared to missing 1–2 days per month prior to surgery. Additionally, the anti-seizure medication (ASM) intake for each patient decreased for 5/7 patients since the matched neurostimulation system was implanted, and did not increase for any of the participants (Fig. 6 c). Importantly, these results dramatically outperform the unmatched thalamic stimulation results reported in literature for the SANTE trial at 7 years follow-up. Indeed, in this previous study, unmatched stimulation (ANT stimulation for seizure onsets other than frontal or temporal) only achieved 39% reduction in seizure frequency, whereas outcome for matched SOZ were comparable to the ones we observed (78% and 86% for temporal and frontal onset, respectively) 26 . In conclusion, our results indicate that hodologically-matched thalamic stimulation provides significant therapeutic benefits, including notable reductions in the frequency and severity of disabling seizures. This approach not only outperforms current standard treatments but also lays the groundwork for more personalized and precise target selection in neuromodulation for medically refractory focal epilepsy, particularly in patients who are not candidates for curative resective or ablative interventions. Discussion Electrical stimulation of the thalamus is undergoing a rapid advance as a treatment opportunity for patients with MRE 12 , 53 , 54 , with indications in both generalized 22 , 53 , 55 and focal epilepsy 9 , 56 – 60 . While in generalized epilepsy the CM nucleus effectively reduce seizure events, in focal epilepsy, this expansion requires a profound understanding of the intricate anatomical and functional thalamocortical interactions with highly heterogeneous and individualized epilepsy phenotypes to improve clinical outcome. Here, we proposed a novel approach to select the optimal targets for thalamic electrical stimulation in focal epilepsy, considering the specificity of SOZs in relation to cortical and subcortical anatomy. Specifically, we demonstrated that leveraging the hodology-based interactions between certain subnuclei of the thalamus and the epileptogenic cortex results in stronger neuromodulator efficacy. This approach highlights the importance of thalamocortical hodology in individualizing stimulation targets, potentially advancing therapeutic outcomes, and paving the way for more effective neuromodulation strategies in patients with focal MRE who are not candidates for curative interventions. Since the early stages of epilepsy surgery and the first mechanistic studies about epileptogenicity, focus has increasingly turned towards deeper brain structures, particularly the thalamus 61 , 62 . Early animal studies by Penfield and Jasper in 1947 identified the thalamus as a potential generator of epileptic cortical activity, proposing the theory of "centrencephalic" epilepsy 18 , 63 . Indeed, this structure mediates reciprocal cortical-subcortical connections, functioning as an 'integrative hub' for functional brain networks 64 . Recently (challenging previous theories), the thalamus has been implicated not in the generation, but in the propagation of generalized and focal onset seizures, due to its role in the bi-hemispheric cortical spread of epileptiform activity 32 , 50 , 65 . In addition, more recent evidence indicates that the thalamus plays a critical role also in the network dynamics of focal-onset epilepsies 19 , 29 , 31 , but lacked a consistent and systematic analysis of the hodological patterns of propagation from different cortical and thalamic regions during ictal events. Our results contribute to this discussion. Indeed, we detected significant correlation with ictal activity in all recorded thalamic nuclei further suggesting the thalamus as a potential neuromodulation target to treat focal MRE. However, correlation between thalamic and cortical SOZ activation was specific to thalamocortical anatomical connection, as identified by high-density tractography as well as electrophysiological assessment, supporting the hypothesis that epileptic propagation patterns during seizures obey predictable anatomo-functional patterns and highlighting the need to carefully select thalamic nuclei for neuromodulation according to individual clinical and electrophysiological data. We applied this mechanistic understanding of the thalamocortical interactions to inform hodology-based electrical stimulation aimed at suppressing epileptiform activity. Our results demonstrated that hodologically-matched thalamic stimulation immediately suppresses interictal epileptiform discharges. Additionally, this targeted stimulation markedly enhances the likelihood of achieving complete suppression of IEDs during stimulation epochs, suggesting robust potential for targeted thalamic stimulation as a therapeutic strategy in epilepsy neuromodulatory management. Immediate effects in n = 14 subjects, paralleled with the drastic clinical improvement in seizure frequency with chronic matched thalamic stimulation (up to 95%) in seven patients, provide promising evidence that a hodology-based stimulation would improve clinical outcome for patients with focal epilepsy. Importantly, the comparative results show a substantial improvement over existing literature on unmatched thalamic neuromodulation (86.5% versus 39% reduction in seizure frequency), highlighting the potential need to reconsider current clinical practices, where a single target is often applied to all seizure onset zones. Mechanistically, our results identified a trend suggesting that the epileptic cortex is the dominant structure at the initiation of seizures, while the thalamus assumes a leading role throughout the course of ictal events, showing higher synchronization towards the middle and termination of the ictal events. Importantly, this pattern was consistent across all thalamic regions investigated. In this regard, in an animal model of limbic epilepsy induced by GABA antagonist (i.e., bicuculline), Aracri and colleagues demonstrated hippocampal activation first, progressing to hippocampus-thalamus synchronization, and then thalamic bursts leading to the termination of seizures 66 . Destexhe and colleagues also highlighted the role of cortically induced coherence of thalamic-generated oscillations 67 . They implicated temporal and spatial thalamocortical interactions, suggesting that, at the onset of seizures, thalamocortical cells are hyperpolarized due to the greater power of excitatory cortical projections to reticular neurons compared to thalamic relay neurons. Once the initial cortical fast phase is abolished, thalamic relay neurons projecting to the cortex depolarize synchronically, reinforcing the coherence of seizure activity arising from the cortex, and promoting the highly synchronous cortical spike phase observed during the middle and end of seizures. Our findings related to the directionality of epileptic activity through thalamocortical circuits align with this mechanistic theory, where information flows to the thalamus from the epileptic cortex at seizure initiation and back to cortical areas during the late seizure phase. This supports the hypothesis that the thalamus plays distinct roles in this dynamic process 68 . Furthermore, it suggests the possibility to design closed-loop stimulation paradigm to restrict the stimulation to precise moment of the seizure development with potentially stronger neuromodulatory effects. In this regard, our findings indicate that the theta to beta frequency range represents a key hallmark of thalamocortical connectivity 69 , which could serve as a biomarker for seizure detection and closed-loop monitoring. The primary limitation of this study is that the clinical assessment of matched electrical stimulation as a chronic therapy was evaluated in just seven patients as an observational retrospective report. Hence, inhomogeneous evaluation of all thalamic nuclei and SOZ occurred. However, these results are paralleled to a thorough investigation of the mechanisms of thalamocortical coupling in focal epilepsy, which support our findings. Additionally, even if limited in number, our evaluation of the effect of matched stimulation in a chronic stage of the treatment results in a drastic clinical improvement well above other literature reports. Larger clinical studies are now necessary to confirm the efficacy of hodologically-matched stimulation in the reduction of seizure frequency and improvements in patients’ quality of life in chronic settings. Another study limitation includes sampling bias associated with explorations of the human brain using intracranial electrodes, where conclusions can only be drawn from areas where electrodes were implanted. Nevertheless, strict inclusion and exclusion criteria minimized this limitation by focusing on subjects with highly localized seizures, many of whom became seizure-free after restricted focal resections. Additionally, the applied parameters of acute stimulation were restricted to few parameters (100 Hz only, 1–3 mA, continuous stimulation), mainly because the focus of the current study was related to the thalamic targeting and not to the different paradigms of electrical stimulation. However, several studies previously identified the optimal stimulation parameters in the high frequency range (> 100Hz) 70–72 , as in our chronic participants. Future studies will need to better refine stimulation parameters and identify optimal stimulation frequencies. In summary, understanding the thalamocortical hodological principles governing the organization of focal seizures and their potential clinical applications is invaluable. Personalizing therapies based on the connectivity of the seizure onset zone to thalamic sub-nuclei and individual seizure characteristics may lead to improved efficacy and more favorable outcomes in the neuromodulation field for epilepsy. Our work addresses the critical unmet clinical need to identify optimal neuromodulation target by demonstrating predictable patterns of connectivity in thalamocortical interactions and guiding optimal thalamic targeting based on location of ictogenesis. This hodology-based approach would impact the current clinical standard in which one target fits all the SOZ (ANT in focal epilepsy, CM in generalized epilepsy), and lay the foundations for highly individualized neuromodulation treatments. Methods Trial and participants information All experimental protocols were approved by the University of Pittsburgh Institutional Review Boards (IRB) (protocol STUDY20070113 and STUDY21020058). All participants included in this study were diagnosed with drug-refractory epilepsy and underwent stereoelectroencephalography (SEEG) implantation at the University of Pittsburgh Medical Center as a part of their standard clinical care. They were hospitalized in the epilepsy monitoring unit (EMU) from January 2021 to May 2024, for 5.04 ± 6.28 days (mean ± SD). All patients received a comprehensive neurological assessment, neuropsychological testing, routine MRI and CT and SEEG implantation. Informed consent was obtained from all patients or legal representatives. In order to study the anatomical and functional properties of the thalamus and the epileptic cortex via SEEG monitoring, we retrospectively included patients (S01-S26) in the present study if: 1) at least one SEEG electrode contact was in the PUL, ANT or VIM/VOP, and 2) at least one of the electrophysiological recordings analyzed in the study was performed (minimum 2 spontaneous seizures, thalamocortical evoked potentials, HDFT, thalamic stimulation). Exclusion criteria encompassed refusal in participating in the study and SEEG implantation complication (severe intracranial hemorrhage). We additionally included 7 participants (S16, S27-S32) that received a chronic thalamic stimulation system as part of their clinical care to evaluate long-term efficacy of matched thalamic stimulation. A comprehensive overview of the patients’ demographics and epilepsy categorization, summary of the experiments, and a summary of the data collected for each patient are reported in Supplementary Table 1&2 . Overall, we included 26 SEEG participants (20 males, 6 females) and 7 chronically implanted patients (S16, S27-S32), of age 34.8 ± 11.3 (mean ± SD). The age at seizure onset was 14.3 ± 13.6 years-old, while the duration of the epilepsy disease was 20.4 ± 13.1 years (mean ± SD). Additional information can be found in the Supplementary Tables 1, 2, 4 . SEEG implantation and localization. The number and location of SEEG electrodes implanted was pre-operatively planned individually for each patient as part of their clinical care and was not influenced by this study. The implant procedure was performed with robotic stereotactic guidance (ROSA, Zimmer-Biomet, Warsaw, IN, USA), applying bone fiducial registration with accuracy < 0.5 mm). Overall, we implanted a total of 5214 contacts (Microdeep® SEEG Electrodes, 12 to 18 channels, DIXI Medical, Marchaux-Chaudefontaine, France), with an average per patient of 200.53 ± 27.87. To target the thalamic nuclei of interest, we used the following stereotactic coordinates: 1) for PUL, in relation to the AC/PC reformatted planes, and having the PC (posterior commissure) as the reference point, coordinates are X: 2–10 mm lateral to the midline plane; Y: 0 to -6 mm posterior to PC; Z: 0 to -6 mm below the AC/PC defined horizontal plane; 2) for ANT, in relation to the AC/PC reformatted planes, and having the MC (mid commissural point) as the reference point, coordinates are X: 6–8 mm lateral to the midline plane; Y: 6 to 7 mm anterior to the MC; Z: 2 to 4 mm above the AC/PC defined horizontal plane; 3) for VIM/VOP, lateral (from about 5 to about 15 mm lateral to the AC/PC line); anterior/posterior (from about 2 to about 10 mm anterior to PC); and dorsal/ventral (from about + 1 to about − 2 mm from the AC/PC plane. To reconstruct the position of the SEEG electrode, we used CURRY software (Compumedics NeuroScan, Hamburg, Germany). First, we marked key reference points on each pre-implantation MRI T1 scan (namely the AC, PC and the nasion). These images were then co-registered with post-implantation CT, and a clinical team member manually identified and labeled all the SEEG contacts from the CT. All the contact locations of the subject MRI were visually inspected and confirmed by a certified neurosurgeon (JGM). For all subject-analysis (Fig. 1 b), the electrode coordinates in the subject space MRI were then translated to coordinates in MNI space for each patient. To confirm the contact location within the thalamic nuclei of interest we used The Human Motor Thalamus Atlas 37 . The assessment of SOZ contacts (as well as non-SOZ contacts used in the following analysis) was confirmed by an experienced neurosurgeon (JGM) ( Supplementary Table 3 ). SEEG recording and stimulation. SEEG data were recorded with Natus Quantum System EEG diagnostic and monitoring system (Natus, Pleasanton, A, USA), with a sampling rate of 2048 Hz (for S26 the sampling rate was 1024 Hz) during the extra-operative monitoring at the Epilepsy Monitoring Unit of the University of Pittsburgh Medical Center (Presbyterian Hospital). For all analysis, we applied SEEG bipolar montage to increase spatial selectivity and reduce noise levels. To deliver electrical stimulation in the thalamic contacts, we used the Nicolet Cortical Stimulator (Natus, Pleasanton, A, USA), which delivers biphasic pulses up to 100 Hz (amplitude was set to 1–3 mA and pulse duration to 60–300 µs for all patients) with bipolar electrode configuration. High-Definition Fiber Tracking. To estimate anatomical projections from the thalamus to various cortical areas, we first performed High-Definition Fiber Tracking (HDFT) of diffusion MRI data. The diffusion images were acquired on a SIEMENS Prisma Fit scanner using a diffusion sequence (2mm isotropic resolution, TE/TR = 99.2 ms/2490 ms, 257 diffusion sampling with maximum b-value 4010 s/mm²). Diffusion tensor estimation and tractography were performed using DSI studio ( http://dsi,studio.labsolver.org ). The accuracy of b-table orientation was examined by comparing fiber orientations with those of a population-averaged template. The tensor metrics were calculated using DWI with b-value lower than 1750 s/mm². For fiber tracking, we used a tracking threshold of 0, angular threshold of 0, and a step size of 0 mm. We utilized seed regions in thalamus to create white matter tracts to regions of interest in the cortex. Seed regions were selected in ANT, VIM/VOP, PUL, and CM based on the extended Human Connectome Project multimodal parcellation atlas. Regions of interest were selected in the frontal lobe, rolandic area, parietal lobe, occipital lobe, and temporal lobe. Tracks with lengths shorter than 30 mm or longer than 1000 mm were discarded. A total of 10,000 tracts were placed. Topology informed pruning was applied to the tractography with 2 interactions to remove false connections. We then quantified the volume of white matter tracts projecting from each thalamic nucleus to each cortical area. We normalized the volume of each white matter projection by the total volume projections from each thalamic nucleus (Fig. 2 a-b, Extended Data Fig. 1 a-b). Thalamocortical evoked potentials. In order to confirm the anatomical organization of thalamocortical fibers with electrophysiological techniques, thalamocortical evoked potentials were elicited by thalamic stimulation of PUL, ANT and VIM/VOP (frequency = 1 Hz, amplitude = 1–3 mA, pulse width = 300 µs) and recorded from all available SEEG channels with Natus recording system. SEEG intracerebral data from all cortical channels were filtered with a band-pass Butterworth 2nd order filter (1-1000 Hz) and DC offset was removed from all recordings. From each stimulation pulse, we extracted 520 ms epochs (20 ms prior and 500 ms after the stimulus) and computed stimulation triggered averages. Each epoch was baseline corrected with the pre-stimulus interval. We calculated the peak-to-peak amplitude of evoked potentials in each epoch as the difference between the maximum and minimum voltage value observed within the first 400 ms from the stimulus. We performed this analysis for all available channels that presented sufficient signal to noise ratio, that were located in the regions of interest for this study (not in thalamus) and that were not located in the white matter (Fig. 2 d-l, Extended Data Fig. 2 ). Acquisition and analysis of spontaneous seizures. Spontaneous seizures acquisition During the extra-operative monitoring, we recorded 216 spontaneous seizures from 23 patients. The remaining 3 (S14, S16, S19) patients either presented technical problems that invalidated the inclusion of the seizures or did not present any spontaneous events during the hospitalization. For our analysis, we collected recordings from 2 minutes before the seizure onset up to 30 seconds after the seizure termination. These time events (onset and termination) were manually marked for each event by certified medical personnel. We used such markers to define 5 epochs of interested, used in subsequent analysis: 1) baseline (corresponding to 20 seconds occurring two minutes before the seizure onset); 2) seizure initiation (including the 10 seconds preceding the first appearance of a tonic discharge in the SOZ contacts); 3) seizure termination (including the last 10 seconds of discharge); 4) middle of seizure (as the epoch between the previous two intervals 5) total seizure duration (defined as the interval between onset and termination). Non-linear correlation analysis For each recorded seizure, we computed the h2 non-linear correlation coefficient, as implemented in Anywave 38 for each pair of thalamus-SOZ bipolar contacts. The h2 coefficient between two signals is a time-resolved measure of their non-linear dependence. The coefficient was computed over a 2 s time window sliding by steps of 1 s. We computed the h2 across the full SEEG spectrum (1-500 Hz), as well as separate frequency bands: delta (1–4 Hz), theta (4–8 Hz), alpha (8–15 Hz), beta (15–30 Hz), gamma (30–45 Hz), gamma2 (55–90 Hz), ripple (80–250 Hz), fast ripple (250–500 Hz). For each time interval of interest (seizure duration, seizure initiation, middle of seizure, seizure termination), we computed the mean h2 (Fig. 3 b, e, f). To computer the change in h2 between matched and unmatched nuclei (Fig. 3 g), we calculated, for each matched and unmatched pair, the percentage increase of h2 in the seizure duration with respect to the baseline epoch, as: $$\:h2\:change\left(\%\right)=100*\left(mean\right(h{2}_{seizure})-mean(h{2}_{baseline}\left)\right)/mean\left(h{2}_{baseline}\right)\:\%$$ We then performed the subtraction between the percentage h2 changes calculated in matched thalamus-SOZ and unmatched thalamus-SOZ pairs. For the correlation between the thalamus and non-SOZ regions, we randomly selected two seizures for each patient. We computed the mean h2 of the thalamus with 3 distinct non-SOZ pair of contacts ( Supplementary Table 3 ). Granger Causality analysis To analyze the directionality of the interactions between the thalamus and the cortex, we used the Granger Causality method 39 . First, we selected all the seizures that lasted at least 10 s and showed an increase in the h2 correlation coefficient in the theta to beta band. We filtered the raw SEEG signals in the 4–30 Hz range and divided each seizure recording in 10 s non-overlapping windows. We computed the GC matrix of each matched thalamus-SOZ pair for every epoch. This computation results in a 4X4 matrix for each epoch, where the second diagonal contains the thalamocortical (th→SOZ) and corticothalamic (SOZ→th,) GC value. In addition to the absolute GC value, we quantified per percentage change in GC from initiation to termination (Fig. 4 c, d, f) as: $$\:\varDelta\:GC\:\left(\%\right)=100*(G{C}_{termination}-G{C}_{initiation})/G{C}_{initiation})\:\%$$ where termination and initiation are the 10 seconds epochs preceding the seizure onset and the seizure termination time stamps, respectively. To avoid misinterpretations given by the variability in the number of seizures recorded from each patient, we highlight the single-subject results in Extended Data Fig. 4 . In a total of 20 patients, 17 showed an increase in the proportion of seizures with the thalamus leading at termination (pie charts in Extended Data Fig. 4 ) and 18 showed higher \(\:\varDelta\:GC\:\) in the thalamocortical (and not corticothalamic) coefficient (bar plot in Extended Data Fig. 4 ). Such finding robustly confirms population-observed trends. Case report: CM connectivity in generalized epilepsy We had the unique opportunity to justify the exclusion of the CM in our study ( Extended Data Fig. 1 ), focused on focal epilepsy, by performing structural and functional connectivity analysis of the CM, as well as the PUL, ANT and VIM/VOP on a patient affected by generalized epilepsy. In particular, the patient was a female in her late 20s, with a severe medically refractory epilepsy, with more than 10 generalized seizures per day. The patient underwent SEEG exploration to investigate the appropriate thalamic target for neuromodulation as a part of her clinical care. We performed h2 correlation analysis of 24 seizures recorded with simultaneous sampling of CM, PUL, ANT, VIM/VOP as well as 13 bipolar electrodes sampling the following brain structures: insula, frontal operculum, superior temporal sulcus, superior temporal gyrus, middle temporal gyrus. We computed the mean h2 coefficient for each thalamic nucleus with the cortical electrodes. Unlike the focal epilepsy patients, where the SOZ only comprises a pair of cortical electrodes, here we averaged the h2 values of each thalamus-cortex pair. High-frequency acute thalamic stimulation. Stimulation testing and protocol During post-operative monitoring in the EMU, we had the unique opportunity to perform acute thalamic stimulation testing on 14 patients of our cohort. Importantly, we could stimulate in both matched and unmatched thalamic nuclei. We applied continuous electrical stimulation with bipolar configuration on the thalamic site of interest, at a frequency of 100 Hz. This value was chosen because it was the highest allowed by the stimulator used, and it is in line with commonly used parameters for thalamic electrical stimulation (≥ 100Hz). The amplitude and pulse width were adjusted individually for each patient according to their comfort levels, but never exceeded the ranges of 1–3 mA and 60–300 us respectively. For clinical safety reasons, the stimulation testing occurred only after the patients exhibited spontaneous seizures and in seizure-free days. The duration of the stimulation train varied from a minimum of 5 seconds to a maximum of 30 seconds, repeated over 1–10 times depending on the patient comfort and clinical state. Preprocessing and IED detection To analyze the immediate effect of thalamic stimulation we extracted the following data: i) for baseline (Stim OFF), we collected 2 minutes of interictal activity in the SOZ channels before any kind of stimulation was delivered to the patient; ii) for active stimulation period (Stim ON), we included the SOZ recording for the whole duration of stimulation up to 10 s after stimulation was stopped 33 . We excluded the first two seconds after the stimulation was delivered to reduce false positive detection due to transient stimulation-induced burst 40 . Unfortunately, for technical reasons, 4 patients had a shorter baseline period (average of 45 s). Prior epoching, we preprocessed the data in Anywave 38 by applying a band-pass filter (1-1000 Hz) and a notch filter (60Hz) and manually removing channels affected by stimulation artifact or other noise sources. On the remaining channels, we applied independent component analysis (ICA) and extracted n components (where n = number of channels − 1). This method is a common technique used to separate independent sources linearly mixed in several sensors. We visually inspected all the components and rejected the ones associated with stimulation artifacts. In this way, we could reliably apply an automatic spikes detection algorithm on the stim ON and stim OFF recordings. Each intracranial recording of interest was then divided in 5 s non-overlapping epochs, on which we computed the detection of epileptiform discharges. To avoid any confound due to the stimulation artifact, we considered as IED only interictal spikes, and not high-frequency oscillations. For this, we used Delphos software 41 , 42 , which detects interictal spikes from SEEG recording from their time-frequency representation. This detector was previously demonstrated to allow for almost 100% specificity and more than 80% sensitivity. To ensure correct spikes detection we 1) computed the analysis also on non-SOZ channels, to confirm that no spikes would be detected and 2) visually inspected the data. We saved the IED number for each Stim ON and Stim OFF epoch. IED analysis To quantify the impact of acute thalamic stimulation on the IED rate, we computed two complementary analyses, based on the data quality: a) if the IED were clearly detectable and frequent (occurring in at least 50% of the baseline epochs), we calculated the IED rate as the number of IED per minute, in each Stim ON and Stim OFF epoch; b) if more than half of the baseline epochs presented no IED, we calculated the probability of finding no spikes across all the Stim ON and Stim OFF epochs ( Supplementary Table 2 ). Even if consistent spiking was observed in baseline epochs, we included in analysis b also a subject (S26) that presented complete suppression of IEDs in stimulation epochs, hence precluding statistical analysis. For both analyses, we computed the percentage difference of these quantities for Stim ON vs stim OFF to characterize the stimulation effect. A beneficial effect of stimulation on epileptiform discharges would result in a smaller IED rate with stimulation ON for analysis a, and a positive change in analysis b, prompting to show that stimulation epochs would more likely show complete suppression of IED. Clinical outcome of chronic thalamic stimulation. We retrospectively evaluated clinical outcome of hodologically matched PUL and ANT stimulation in seven individuals (S16, S27-S32) who received a chronic neuromodulation system as part of their clinical care. Demographics and epilepsy information for these patients is reported in Supplementary Table 1 . All patients received a bilateral implant of the PUL or ANT, according to their SOZ location and as indicated by their clinical team following SEEG monitoring. Six patients received a Responsive Neurostimulator (RNS), while one patient received a DBS system (Medtronic Percept). Stimulation parameters were individually tailored for each patient according to their clinical need ( Supplementary Table 4 ). No clinical decisions were based upon this research. We retrospectively evaluated the percentage reduction in disabling seizures’ frequency and the change in antiseizure medication. The follow-up window ranged from 9 to 30 months. Data analysis and statistical procedures. All figures and analysis were performed using Matlab 2023b (Mathworks, California, US). For all box plots reported in this manuscript, the whiskers extend to the maximum spread not considering outliers, central, top, and bottom lines represent median, 25th, and 75th percentile, respectively. For all the analysis, we used the bootstrap method, which does not rely on distributional assumptions of the data, but rather resamples the quantities of interest to achieve empirical confidence intervals. For each comparison, we created a bootstrap sample (for n = 10000 repetitions) by drawing a sample with replacement from the actual data points (n = 10000 repetitions), and calculated the difference in means of the resampled data. We then applied two-tailed bootstrapping with significance levels of 0.05 (95% confidence interval), 0.01 (99% confidence interval), or 0.001 (99.9% confidence interval). The null hypothesis of no difference in the mean was rejected if 0 was not included in the confidence interval of the corresponding alpha value. If multiple comparisons were performed at once, we used a Bonferroni correction by dividing the alpha value by the number of pairwise comparisons being performed. Declarations Funding This research received financial support from NIH (1R01NS122927-01A1- Gonzalez-Martinez & Sridevi Sarma, PIs and 1R01NS131428-01A1 - Elvira Pirondini & Gonzalez-Martinez, PIs). Author contributions EP and JGM conceived the study. EP and JGM secured funding. AD, JGM, and EP designed the experiments. AD analyzed the data, with help from SN for correlation analysis and JH for DTI analysis. TA and JGM implemented patient recruitment, eligibility and monitoring and coordinated study management. AD, TA and JGM collected the data. AG, SS and ER helped with data management. AD, EP, JGM and TA helped with data interpretation. AD, JH, EP and JGM, wrote the paper and all authors contributed to its editing. All authors interpreted the data, contributed to the writing of the submitted manuscript, had final responsibility for the decision to submit the manuscript for publication, and vouch for the accuracy and completeness of the data and for the fidelity of the trial to the protocol. Competing interests The authors declare no conflicts of interests in relation to this work. Data availability The main data supporting the results in this study are available within the paper and its supplementary information. All data generated in this study and software will be uploaded in a public repository upon acceptance of the manuscript. Raw data will be available upon reasonable request to the corresponding author. Acknowledgements We would like to thank all the Epilepsy Monitoring Unit staff and Neurosurgery Staff of the University of Pittsburgh Medical Center for their support. We wish to thank all the participants and their families, who agreed to participate in the study and for the great insights into the real needs of people living with epilepsy. References Kwan, P. & Brodie, M. J. Definition of refractory epilepsy: defining the indefinable? Lancet Neurol. 9 , 27–29 (2010). Thijs, R. D., Surges, R., O’Brien, T. J. & Sander, J. W. Epilepsy in adults. Lancet Lond. Engl. 393 , 689–701 (2019). Gupta, S., Ryvlin, P., Faught, E., Tsong, W. & Kwan, P. Understanding the burden of focal epilepsy as a function of seizure frequency in the United States, Europe, and Brazil. Epilepsia Open 2 , 199–213 (2017). Gloor, P. & Fariello, R. G. Generalized epilepsy: some of its cellular mechanisms differ from those of focal epilepsy. Trends Neurosci. 11 , 63–68 (1988). Badawy, R. A. B., Curatolo, J. M., Newton, M., Berkovic, S. F. & Macdonell, R. A. L. Changes in cortical excitability differentiate generalized and focal epilepsy. Ann. Neurol. 61 , 324–331 (2007). French, J. A. Refractory Epilepsy: Clinical Overview. Epilepsia 48 , 3–7 (2007). Li, A. et al. Neural fragility as an EEG marker of the seizure onset zone. Nat. Neurosci. 24 , 1465–1474 (2021). Velasco, F. et al. Deep brain stimulation for treatment of the epilepsies: the centromedian thalamic target. in Operative Neuromodulation: Volume 2: Neural Networks Surgery (eds. Sakas, D. E. & Simpson, B. A.) 337–342 (Springer, Vienna, 2007). doi:10.1007/978-3-211-33081-4_38. Fisher, R. et al. Electrical stimulation of the anterior nucleus of thalamus for treatment of refractory epilepsy. Epilepsia 51 , 899–908 (2010). Fisher, R. S. Deep brain stimulation of thalamus for epilepsy. Neurobiol. Dis. 179 , 106045 (2023). Fisher, R. S. & Velasco, A. L. Electrical brain stimulation for epilepsy. Nat. Rev. Neurol. 10 , 261–270 (2014). Ryvlin, P., Rheims, S., Hirsch, L. J., Sokolov, A. & Jehi, L. Neuromodulation in epilepsy: state-of-the-art approved therapies. Lancet Neurol. 20 , 1038–1047 (2021). Shine, J. M., Lewis, L. D., Garrett, D. D. & Hwang, K. The impact of the human thalamus on brain-wide information processing. Nat. Rev. Neurosci. 24 , 416–430 (2023). Steriade, M. Synchronized activities of coupled oscillators in the cerebral cortex and thalamus at different levels of vigilance. Cereb. Cortex 7 , 583–604 (1997). Steriade, M. Sleep, epilepsy and thalamic reticular inhibitory neurons. Trends Neurosci. 28 , 317–324 (2005). Steriade, M. & Contreras, D. Spike-Wave Complexes and Fast Components of Cortically Generated Seizures. I. Role of Neocortex and Thalamus. J. Neurophysiol. 80 , 1439–1455 (1998). Behrens, T. E. J. et al. Non-invasive mapping of connections between human thalamus and cortex using diffusion imaging. Nat. Neurosci. 6 , 750–757 (2003). Penfield, W. Epileptic automatism and the centrencephalic integrating system. Res. Publ. - Assoc. Res. Nerv. Ment. Dis. 30 , 513–528 (1952). Arthuis, M. et al. Impaired consciousness during temporal lobe seizures is related to increased long-distance cortical–subcortical synchronization. Brain 132 , 2091–2101 (2009). Zhang, Y. et al. Thalamocortical structural connectivity abnormalities in drug-resistant generalized epilepsy: A diffusion tensor imaging study. Brain Res. 1727 , 146558 (2020). Salanova, V. et al. Long-term efficacy and safety of thalamic stimulation for drug-resistant partial epilepsy. Neurology 84 , 1017–1025 (2015). Cukiert, A., Cukiert, C. M., Burattini, J. A. & Mariani, P. P. Seizure outcome during bilateral, continuous, thalamic centromedian nuclei deep brain stimulation in patients with generalized epilepsy: a prospective, open-label study. Seizure 81 , 304–309 (2020). Valentín, A. et al. Deep brain stimulation of the centromedian thalamic nucleus for the treatment of generalized and frontal epilepsies. Epilepsia 54 , 1823–1833 (2013). Cooper, I. S. et al. Evoked metabolic responses in the limbic-striate system produced by stimulation of anterior thalamic nucleus in man. Int. J. Neurol. 18 , 179–187 (1984). Kaufmann, E. et al. Long-term evaluation of anterior thalamic deep brain stimulation for epilepsy in the European MORE registry. Epilepsia 65 , 2438–2458 (2024). Salanova, V. et al. The SANTÉ study at 10 years of follow-up: Effectiveness, safety, and sudden unexpected death in epilepsy. Epilepsia 62 , 1306–1317 (2021). Weitz, A. J. et al. Optogenetic fMRI reveals distinct, frequency-dependent networks recruited by dorsal and intermediate hippocampus stimulations. NeuroImage 107 , 229–241 (2015). Bhattacharyya, K. B. James Wenceslaus Papez, His Circuit, and Emotion. Ann. Indian Acad. Neurol. 20 , 207–210 (2017). Guye, M. et al. The role of corticothalamic coupling in human temporal lobe epilepsy. Brain J. Neurol. 129 , 1917–1928 (2006). Capecchi, F., Mothersill, I. & Imbach, L. L. The medial pulvinar as a subcortical relay in temporal lobe status epilepticus. Seizure 81 , 276–279 (2020). Pizzo, F. et al. The Ictal Signature of Thalamus and Basal Ganglia in Focal Epilepsy. Neurology 96 , e280–e293 (2021). Rosenberg, D. S. et al. Involvement of Medial Pulvinar Thalamic Nucleus in Human Temporal Lobe Seizures. Epilepsia 47 , 98–107 (2006). Ikegaya, N. et al. Thalamic stereoelectroencephalography for neuromodulation target selection: Proof of concept and review of literature of pulvinar direct electrical stimulation. Epilepsia 65 , e79–e86 (2024). Burdette, D., Mirro, E. A., Lawrence, M. & Patra, S. E. Brain-responsive corticothalamic stimulation in the pulvinar nucleus for the treatment of regional neocortical epilepsy: A case series. Epilepsia Open 6 , 611–617 (2021). Vilela-Filho, O. et al. Ventral intermediate nucleus deep brain stimulation for treatment-resistant focal aware motor seizures: illustrative case. J. Neurosurg. Case Lessons 5 , CASE2320 (2023). Cooper, I. S., Amin, I. & Gilman, S. The effect of chronic cerebellar stimulation upon epilepsy in man. Trans. Am. Neurol. Assoc. 98 , 192–196 (1973). Ilinsky, I. et al. Human Motor Thalamus Reconstructed in 3D from Continuous Sagittal Sections with Identified Subcortical Afferent Territories. eNeuro 5 , ENEURO.0060-18.2018 (2018). Colombet, B., Woodman, M., Badier, J. M. & Bénar, C. G. AnyWave: a cross-platform and modular software for visualizing and processing electrophysiological signals. J. Neurosci. Methods 242 , 118–126 (2015). Granger, C. W. J. Investigating Causal Relations by Econometric Models and Cross-spectral Methods. Econometrica 37 , 424–438 (1969). Milosevic, L. et al. Physiological mechanisms of thalamic ventral intermediate nucleus stimulation for tremor suppression. Brain 141 , 2142–2155 (2018). Roehri, N., Lina, J.-M., Mosher, J. C., Bartolomei, F. & Benar, C.-G. Time-Frequency Strategies for Increasing High-Frequency Oscillation Detectability in Intracerebral EEG. IEEE Trans. Biomed. Eng. 63 , 2595–2606 (2016). Roehri, N., Pizzo, F., Bartolomei, F., Wendling, F. & Bénar, C.-G. What are the assets and weaknesses of HFO detectors? A benchmark framework based on realistic simulations. PLOS ONE 12 , e0174702 (2017). Velasco, F. et al. Centromedian Nucleus and Epilepsy. J. Clin. Neurophysiol. 38 , 485 (2021). Yamada, K. et al. Somatotopic Organization of Thalamocortical Projection Fibers as Assessed with MR Tractography. Radiology 242 , 840–845 (2007). Jbabdi, S., Woolrich, M. W. & Behrens, T. E. J. Multiple-subjects connectivity-based parcellation using hierarchical Dirichlet process mixture models. NeuroImage 44 , 373–384 (2009). Aupy, J. et al. Cortico-striatal synchronization in human focal seizures. Brain 142 , 1282–1295 (2019). Lopes da Silva, F., Pijn, J. P. & Boeijinga, P. Interdependence of EEG signals: linear vs. nonlinear associations and the significance of time delays and phase shifts. Brain Topogr. 2 , 9–18 (1989). Bartolomei, F., Wendling, F., Bellanger, J.-J., Régis, J. & Chauvel, P. Neural networks involving the medial temporal structures in temporal lobe epilepsy. Clin. Neurophysiol. 112 , 1746–1760 (2001). Bertram, E. H., Mangan, P. S., Zhang, D., Scott, C. A. & Williamson, J. M. The Midline Thalamus: Alterations and a Potential Role in Limbic Epilepsy. Epilepsia 42 , 967–978 (2001). Cassidy, R. M. & Gale, K. Mediodorsal thalamus plays a critical role in the development of limbic motor seizures. J. Neurosci. Off. J. Soc. Neurosci. 18 , 9002–9009 (1998). Evangelista, E. et al. Does the Thalamo-Cortical Synchrony Play a Role in Seizure Termination? Front. Neurol. 6 , (2015). Soulier, H. et al. The anterior and pulvinar thalamic nuclei interactions in mesial temporal lobe seizure networks. Clin. Neurophysiol. 150 , 176–183 (2023). Vetkas, A. et al. Deep brain stimulation targets in epilepsy: Systematic review and meta-analysis of anterior and centromedian thalamic nuclei and hippocampus. Epilepsia 63 , 513–524 (2022). Krauss, J. K. et al. Technology of deep brain stimulation: current status and future directions. Nat. Rev. Neurol. 17 , 75–87 (2021). Agashe, S. et al. Centromedian Nucleus of the Thalamus Deep Brain Stimulation for Genetic Generalized Epilepsy: A Case Report and Review of Literature. Front. Hum. Neurosci. 16 , (2022). Graves, N. M. & Fisher, R. S. Chapter 16 Neurostimulation for Epilepsy, Including a Pilot Study of Anterior Nucleus Stimulation. Neurosurgery 52 , 127 (2005). Kim, S. H. et al. Long-term follow-up of anterior thalamic deep brain stimulation in epilepsy: A 11-year, single center experience. Seizure 52 , 154–161 (2017). Guo, W. et al. Defining the optimal target for anterior thalamic deep brain stimulation in patients with drug-refractory epilepsy. J. Neurosurg. 134 , 1054–1063 (2021). Oh, Y.-S. et al. Cognitive improvement after long-term electrical stimulation of bilateral anterior thalamic nucleus in refractory epilepsy patients. Seizure 21 , 183–187 (2012). Park, H. R. et al. The Role of Anterior Thalamic Deep Brain Stimulation as an Alternative Therapy in Patients with Previously Failed Vagus Nerve Stimulation for Refractory Epilepsy. Stereotact. Funct. Neurosurg. 97 , 176–182 (2019). Norden, A. D. & Blumenfeld, H. The role of subcortical structures in human epilepsy. Epilepsy Behav. 3 , 219–231 (2002). Avoli, M. A brief history on the oscillating roles of thalamus and cortex in absence seizures. Epilepsia 53 , 779–789 (2012). W, P. Highest level seizures. Res Publ Nerv Ment Dis 26 , 252–271 (1947). Hwang, K., Bertolero, M. A., Liu, W. B. & D’Esposito, M. The Human Thalamus Is an Integrative Hub for Functional Brain Networks. J. Neurosci. 37 , 5594–5607 (2017). Wu, T. Q. et al. Multisite thalamic recordings to characterize seizure propagation in the human brain. Brain 146 , 2792–2802 (2023). Aracri, P., de Curtis, M., Forcaia, G. & Uva, L. Enhanced thalamo-hippocampal synchronization during focal limbic seizures. Epilepsia 59 , 1774–1784 (2018). Destexhe, A., Contreras, D. & Steriade, M. Mechanisms underlying the synchronizing action of corticothalamic feedback through inhibition of thalamic relay cells. J. Neurophysiol. 79 , 999–1016 (1998). Schindler, K., Leung, H., Elger, C. E. & Lehnertz, K. Assessing seizure dynamics by analysing the correlation structure of multichannel intracranial EEG. Brain 130 , 65–77 (2007). Aiello, G. et al. Functional network dynamics between the anterior thalamus and the cortex in deep brain stimulation for epilepsy. Brain 146 , 4717–4735 (2023). Montgomery, E. B. & Gale, J. T. Mechanisms of action of deep brain stimulation(DBS). Neurosci. Biobehav. Rev. 32 , 388–407 (2008). Boëx, C., Vulliémoz, S., Spinelli, L., Pollo, C. & Seeck, M. High and low frequency electrical stimulation in non-lesional temporal lobe epilepsy. Seizure 16 , 664–669 (2007). Wyckhuys, T., Raedt, R., Vonck, K., Wadman, W. & Boon, P. Comparison of hippocampal Deep Brain Stimulation with high (130 Hz) and low frequency (5 Hz) on afterdischarges in kindled rats. Epilepsy Res. 88 , 239–246 (2010). Additional Declarations There is NO Competing Interest. Supplementary Files ExtendedDataFigures.docx SupplementaryInformation.docx Cite Share Download PDF Status: Published Journal Publication published 20 Oct, 2025 Read the published version in Nature Communications → Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-5507011","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Article","associatedPublications":[],"authors":[{"id":382797268,"identity":"35c42a16-15df-4b32-a4a4-729ff8aa581a","order_by":0,"name":"Jorge Gonzalez-Martinez","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA4klEQVRIiWNgGAWjYNCCAgYGxgbmA1BeAjFaDEBa2GBKidXCwMBjQJwWeffewx9+GNjYM7f3fJMuzLFj4GfPMcCrxfDMuQTDHoM0Zsaes9ukZ25LZpDseUNAy4wcgwQeg8NsjDNyt0nzbjvAYHCDkC1ALQf/GPznYZyR8wysxZ6QFnmJHMNmHoMDEkAtbBBbJAhoMeA5Y8wsY5BswNhzzNga6BceiTPPCvDb0t5j/PFNhZ29YXvzw9uF2+zk+NuTN+C35QDMUw0MDMxAmgevcrAtDTAGA0TLKBgFo2AUjAIMAACkD0GuAQql1QAAAABJRU5ErkJggg==","orcid":"","institution":"Cleveland Clinic","correspondingAuthor":true,"submittingAuthor":false,"prefix":"","firstName":"Jorge","middleName":"","lastName":"Gonzalez-Martinez","suffix":""},{"id":382797269,"identity":"15e19f9e-060e-42d2-9cde-e4ce10dce271","order_by":1,"name":"Arianna Damiani","email":"","orcid":"https://orcid.org/0009-0002-8521-865X","institution":"University of Pittsburgh","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Arianna","middleName":"","lastName":"Damiani","suffix":""},{"id":382797270,"identity":"a37ab8ce-20f5-4bdf-961c-5c003b625ba4","order_by":2,"name":"Sirisha Nouduri","email":"","orcid":"","institution":"University of Pittsburgh","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Sirisha","middleName":"","lastName":"Nouduri","suffix":""},{"id":382797271,"identity":"05be0b65-d23f-4e63-b39f-2386230ef204","order_by":3,"name":"Jonathan Ho","email":"","orcid":"","institution":"University of Pittsburgh","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Jonathan","middleName":"","lastName":"Ho","suffix":""},{"id":382797272,"identity":"29365bb5-6cd5-4e6d-ba9f-c03f0514846b","order_by":4,"name":"Steven Salazar","email":"","orcid":"","institution":"University of Pittsburgh","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Steven","middleName":"","lastName":"Salazar","suffix":""},{"id":382797273,"identity":"20fac83c-ee3a-4466-bd3d-5ec91cf93f54","order_by":5,"name":"Aude Jegou","email":"","orcid":"","institution":"University of Pittsburgh","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Aude","middleName":"","lastName":"Jegou","suffix":""},{"id":382797274,"identity":"9f1ff204-0a92-4544-b91e-2e16f367f72d","order_by":6,"name":"Eliza Reedy","email":"","orcid":"","institution":"University of Pittsburgh","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Eliza","middleName":"","lastName":"Reedy","suffix":""},{"id":382797275,"identity":"96f8e6d3-aecf-4280-941c-3a6c16bc26d2","order_by":7,"name":"Naoki Ikegaya","email":"","orcid":"","institution":"University of Pittsburgh","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Naoki","middleName":"","lastName":"Ikegaya","suffix":""},{"id":382797276,"identity":"5ee9e022-c02e-4670-aec7-4bb2d6063d6f","order_by":8,"name":"Sridevi Sarma","email":"","orcid":"","institution":"","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Sridevi","middleName":"","lastName":"Sarma","suffix":""},{"id":382797277,"identity":"c8062ab6-a520-43dd-b8f8-243719696b16","order_by":9,"name":"Thandar Aung","email":"","orcid":"","institution":"University of Pittsburgh","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Thandar","middleName":"","lastName":"Aung","suffix":""},{"id":382797278,"identity":"13ad450b-08c7-4ca8-a6bf-9353435d08ff","order_by":10,"name":"Elvira Pirondini","email":"","orcid":"https://orcid.org/0000-0002-8078-3951","institution":"University of Pittsburgh","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Elvira","middleName":"","lastName":"Pirondini","suffix":""}],"badges":[],"createdAt":"2024-11-22 21:35:11","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-5507011/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-5507011/v1","draftVersion":[],"editorialEvents":[{"content":"https://doi.org/10.1038/s41467-025-64922-w","type":"published","date":"2025-10-20T04:00:00+00:00"}],"editorialNote":"","failedWorkflow":false,"files":[{"id":70129641,"identity":"c3e76c12-5bf6-4ce1-aea6-d0b24d36ad63","added_by":"auto","created_at":"2024-11-28 15:50:27","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":676939,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eStudy design and patient demographics|\u003c/strong\u003e (\u003cstrong\u003ea) \u003c/strong\u003e\u003cem\u003eLeft:\u003c/em\u003e\u003cstrong\u003e \u003c/strong\u003eSchematic of hodological matching of PUL, ANT and VIM/VOP to specific cortical anatomical regions. \u003cem\u003eRight: \u003c/em\u003eanatomical localization of PUL, ANT and VIM/VOP nuclei on MNI space average T1-weighted MRI with axial view and magnification of the thalamic nuclei. Thalamic parcellation is atlas-guided\u003csup\u003e37\u003c/sup\u003e (A=anterior, P=posterior, M=medial, L=lateral). (\u003cstrong\u003eb\u003c/strong\u003e) \u003cem\u003eTop right:\u003c/em\u003e axial view of reconstruction of SEEG electrodes targeting thalamic nuclei of interest. \u003cem\u003eTop left:\u003c/em\u003e magnification of the contacts located in the PUL (blue), ANT (yellow) and VIM/VOP (red). \u003cem\u003eBottom: \u003c/em\u003eside and top view of the left and right hemisphere SEEG thalamic implantations. (\u003cstrong\u003ec\u003c/strong\u003e) Stratification of SEEG patients (S01-S26). \u003cem\u003eTop Left:\u003c/em\u003e pie chart of the lobar location of the SOZ. \u003cem\u003eTop Right: \u003c/em\u003epie chart of 37 thalamic SEEG contacts in PUL, ANT and VIM/VOP. A higher number of thalamic contacts than patients implies that some patients had more than one nucleus of interest implanted. \u003cem\u003eBottom:\u003c/em\u003e Stratification of the 216 recorded spontaneous seizures according to the lobar location of different SOZs (left) and the respective thalamic coverage (right). A different distribution from the top panel implies a different number of seizures for each patient, as reported in Supplementary Table 2. (\u003cstrong\u003ed\u003c/strong\u003e) \u003cem\u003eTop: \u003c/em\u003eStratification of 7 chronically implanted patients (S16, S27-32) according to the lobar location of the SOZ (left) and implanted thalamic nucleus (right). (\u003cstrong\u003ee\u003c/strong\u003e) \u003cem\u003e\u0026nbsp;\u003c/em\u003ehistogram of follow-up years of chronically implanted patients. (\u003cstrong\u003ef\u003c/strong\u003e) Design and experiments of the study. Patients received SEEG implantation and were stratified according to their SOZ and thalamic coverage. \u003cem\u003eBottom, left to right\u003c/em\u003e: we investigated the anatomical (imaging and electrophysiology), and functional coupling during seizures involving the thalamus, hence defining hodologically matching and unmatching nuclei. We then tested the effect of acute thalamic stimulation and finally the clinical outcome following chronic thalamic stimulation.\u0026nbsp;\u0026nbsp; \u0026nbsp;\u0026nbsp;\u003c/p\u003e","description":"","filename":"image1.png","url":"https://assets-eu.researchsquare.com/files/rs-5507011/v1/c7941ee4643e8b9f601eb5e2.png"},{"id":70129648,"identity":"c9d3c4eb-3b8a-4cc5-b1f4-8f8c053f17fc","added_by":"auto","created_at":"2024-11-28 15:50:28","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":804498,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eAnatomical connectivity of PUL, ANT and VIM/VOP via neuroimaging and electrophysiology|\u003c/strong\u003e (\u003cstrong\u003ea\u003c/strong\u003e) Representative example of high-definition fiber tracking (HDFT) from PUL, ANT and VIM/VOP nuclei (S01). (\u003cstrong\u003eb\u003c/strong\u003e) Volume of thalamocortical projections (mean ± standard error, n= 8 patients) from each nucleus to each cortical lobe normalized by the total volume of fibers projecting from each nucleus (mean ± standard error). (\u003cstrong\u003ec\u003c/strong\u003e) Group-level analysis of thalamocortical EPs. For each patient, the mean peak to peak amplitude is calculated across the contacts of a representative electrode for each cortical lobe. The normalized peak to peak amplitude of EP is color coded according to the stimulating nucleus evoking the response (blue for PUL, yellow for ANT, red for VIM/VOP). (\u003cstrong\u003ed, g, j\u003c/strong\u003e) Anatomical localization of SEEG implantation of representative patients (S01, S22, S13) that were tested for EPs with stimulation of PUL, ANT and VIM/VOP respectively. Each electrode is color-coded according to the brain lobe, and thalamic electrodes are indicated by an arrow. (\u003cstrong\u003ee, h, k\u003c/strong\u003e) Raw traces of stimulation triggered averages in different brain lobes from PUL, ANT and VIM/VOP respectively. Specific anatomical regions are: frontal operculum (pars opercularis), supplementary motor area, angular gyrus, inferior temporal sulcus (panel e, left to right); superior frontal sulcus, motor cingulate, posterior cingulate, posterior uncus (panel h, left to right); anterior cingulate, supplementary motor area, precuneus (panel k, left to right); (\u003cstrong\u003ef, i, l\u003c/strong\u003e) Violin plot of peak to peak amplitude of EPs in different brain lobes for PUL, ANT and VIM/VOP respectively. For all panels, statistical significance was assessed with two-tail bootstrapping with Bonferroni correction: p\u0026lt;0.05 (*), p\u0026lt;0.01 (**), p\u0026lt;0.001(***).\u003c/p\u003e","description":"","filename":"image2.png","url":"https://assets-eu.researchsquare.com/files/rs-5507011/v1/bb3d17f6b83dcaf151aeaa04.png"},{"id":70129642,"identity":"8fe47203-c6a2-433c-9ce9-bd353ca60d5c","added_by":"auto","created_at":"2024-11-28 15:50:28","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":570331,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eFunctional coupling of the PUL, ANT and VIM/VOP nuclei with SOZ during seizures| \u003c/strong\u003e(\u003cstrong\u003ea\u003c/strong\u003e) Representative examples of thalamocortical connectivity profiles of PUL (right, S03), ANT (middle, S24) and VIM/VOP (left, S26) with matched SOZ. For all nuclei, the top panel shows the cortical recording in the contacts of the SOZ, the middle panel shows the thalamic recording, and the bottom panel shows the h2 correlation coefficient. The gray shaded area includes the seizure duration (from onset to termination). (\u003cstrong\u003eb\u003c/strong\u003e) Boxplot of the mean h2 correlation for each seizure in matched and unmatched thalami, compared to mean h2 correlation during a baseline epoch preceding the seizure. Three representative patients are shown (right: S01, matched PUL, unmatched ANT; middle: S22, matched ANT, unmatched PUL; left: S12, matched PUL, unmatched VIM/VOP). (\u003cstrong\u003ec\u003c/strong\u003e) Raw traces of h2 correlation coefficient (smoothed for visualization) for a representative seizure for each patient of panel b. The gray shaded area includes the seizure duration (from onset to termination). (\u003cstrong\u003ed\u003c/strong\u003e) Schematic illustrating the hodologically matching (bold colored lines) and unmatching (dashed gray lines) thalamocortical connections. (\u003cstrong\u003ee\u003c/strong\u003e) Violin plot showing the mean h2 coefficient across seizures for each matched (top) and unmatched (bottom) thalamus-SOZ pair. (\u003cstrong\u003ef\u003c/strong\u003e) Box plot of mean h2 coefficient across seizures during initiation, middle and termination phase of each seizure. (\u003cstrong\u003eg\u003c/strong\u003e) Percentage difference of the change in h2 coefficient during seizures between matched and unmatched nuclei, across different frequency bands of interest. For all boxplots, the whiskers extend to the maximum spread not considering outliers, central, top, and bottom lines represent median, 25\u003csup\u003eth\u003c/sup\u003e, and 75\u003csup\u003eth\u003c/sup\u003e percentile, respectively. For all panels, statistical significance was assessed with two-tail bootstrapping with Bonferroni correction: p\u0026lt;0.05 (*), p\u0026lt;0.01 (**), p\u0026lt;0.001(***).\u003c/p\u003e","description":"","filename":"image3.png","url":"https://assets-eu.researchsquare.com/files/rs-5507011/v1/9544da4c2fde066aa5015753.png"},{"id":70130984,"identity":"b1d817bd-b37c-4a0e-b5c5-9d3db36d2183","added_by":"auto","created_at":"2024-11-28 15:58:28","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":433454,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eAnalysis of thalamocortical versus corticothalamic interactions during seizures| \u003c/strong\u003e(\u003cstrong\u003ea\u003c/strong\u003e) Representative seizure recording (S01). The red bar indicates the full seizure duration, while the shaded gray regions show the initiation and termination epoch. On such windows, we calculated the GC coefficient matrix that illustrates connectivity from and to the thalamic nuclei to a matched SOZ. (\u003cstrong\u003eb\u003c/strong\u003e) Pie charts illustrating the proportion of seizures with higher corticothalamic (SOZ→th, light gray) or thalamocortical (th→SOZ, dark gray) GC coefficient for the initiation (left) and termination (right) epoch. (\u003cstrong\u003ec\u003c/strong\u003e) (\u003cem\u003eLeft\u003c/em\u003e) Plot of the evolution of corticothalamic and thalamocortical GC coefficients for 10 s epochs starting two minutes before a representative seizure. The arrows indicate the relative change of GC from initiation to termination. (\u003cem\u003eRight\u003c/em\u003e) Quantification of the relative percentage change of GC coefficient from initiation to termination, for corticothalamic (light gray) and thalamocortical (dark gray) direction. (\u003cstrong\u003ed\u003c/strong\u003e) Boxplot of the relative percentage change of GC coefficient from initiation to termination for all subjects. Each data point represents the mean value across all the seizures for each subject. (\u003cstrong\u003ee\u003c/strong\u003e) Pie charts illustrating the proportion of seizures with higher corticothalamic (SOZ→th, light color) or thalamocortical (th→SOZ, dark color) GC coefficient for the initiation (left) and termination (right) epoch, separated from the three nuclei of interest (PUL at left, ANT at middle, VIM/VOP at right). (\u003cstrong\u003ef\u003c/strong\u003e) Boxplot of the relative percentage change of GC coefficient from initiation to termination for all the seizures separated for the matched thalamic nucleus investigated (PUL at left, ANT at middle, VIM/VOP at right). For PUL n=38 and n=36 seizures, for ANT n=14 and n=15 seizures, for VIM/VOP n=38 and n=35 seizures. For all boxplots, the whiskers extend to the maximum spread not considering outliers, central, top, and bottom lines represent median, 25th, and 75th percentile, respectively. Outliers were removed with quartile method (see Methods). For all panels, statistical significance was assessed with two-tail bootstrapping with Bonferroni correction: p\u0026lt;0.05 (*), p\u0026lt;0.01 (**), p\u0026lt;0.001(***).\u003c/p\u003e","description":"","filename":"image4.png","url":"https://assets-eu.researchsquare.com/files/rs-5507011/v1/b4d0bb39582d2929c4f51b4b.png"},{"id":70129657,"identity":"f073723d-4f84-4664-8aae-6900e843fde8","added_by":"auto","created_at":"2024-11-28 15:50:28","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":476141,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eHodologically-matched electrical stimulation of the thalamus better reduces pathological IEDs | (a\u003c/strong\u003e) Representative example of an interictal SEEG recording of a SOZ contact, with clearly visible IED (spikes) and their time-frequency representation exploited for automatic detection. (\u003cstrong\u003eb\u003c/strong\u003e) (\u003cem\u003eLeft\u003c/em\u003e) boxplot for IED rate/minute in stimulation OFF (n=18) and stimulation ON (n=8) epochs of the PUL (matched nucleus). (\u003cem\u003eRight\u003c/em\u003e) representative stim OFF (\u003cem\u003etop\u003c/em\u003e) and stim ON epoch (\u003cem\u003ebottom\u003c/em\u003e) with matched PUL stimulation shows immediate suppression of IED. (\u003cstrong\u003ec\u003c/strong\u003e) Representative example, as in panel b, for matched ANT stimulation (n=7 and n=17 for stimulation OFF and ON, respectively). (\u003cstrong\u003ed\u003c/strong\u003e) Representative example, as in panel b and c, for matched VIM/VOP stimulation (n=8 and n=21 for stimulation OFF and ON, respectively). (\u003cstrong\u003ee\u003c/strong\u003e) Boxplot of representative examples of unmatched thalamic stimulation on IED rate for ANT (left, n=18 vs n=8 for stimulation OFF and ON, respectively), and VIM/VOP (right, n=19 vs n=41 for stimulation OFF and ON, respectively). (\u003cstrong\u003ef\u003c/strong\u003e) Percentage change relative to baseline of the probability of IEDs suppression in matched and unmatched thalamic stimulation in two representative patients. (\u003cstrong\u003eg\u003c/strong\u003e) All-subject analysis of the IED rate percentage variation to baseline with matched (\u003cem\u003eleft\u003c/em\u003e) and unmatched (\u003cem\u003eright\u003c/em\u003e) stimulation. (\u003cstrong\u003eh\u003c/strong\u003e) All-subject analysis of the percentage change relative to baseline of the probability of IEDs suppression in matched (\u003cem\u003eleft\u003c/em\u003e) and unmatched (\u003cem\u003eright\u003c/em\u003e) stimulation. For all boxplots, the whiskers extend to the maximum spread not considering outliers, central, top, and bottom lines represent median, 25th, and 75th percentile, respectively. Outliers were removed with quartile method (see Methods). For all panels, statistical significance was assessed with two-tail bootstrapping with Bonferroni correction: p\u0026lt;0.05 (*), p\u0026lt;0.01 (**), p\u0026lt;0.001(***).\u003c/p\u003e","description":"","filename":"image5.png","url":"https://assets-eu.researchsquare.com/files/rs-5507011/v1/28564f802b32080ec70d0279.png"},{"id":70129672,"identity":"4bd2f7bc-f41d-4308-9a29-e610f9ba5620","added_by":"auto","created_at":"2024-11-28 15:50:29","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":238421,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eHodologically-matched thalamic stimulation as a neuromodulation therapy in chronic epilepsy treatment| \u003c/strong\u003e(\u003cstrong\u003ea-b\u003c/strong\u003e) Bar plot and radar plot of the percentage of seizure frequency reduction for S16, S27-S32, compared to the median seizure frequency reduction reported by the SANTE trial at 7 years follow up for frontal, temporal and other lobes. (\u003cstrong\u003ec\u003c/strong\u003e) Scatter plot of the number of antiseizure medication for each patient prior and after the matched thalamic neurostimulation implant. Dashed lines indicate no change in the numbers of ASM.\u003c/p\u003e","description":"","filename":"image6.png","url":"https://assets-eu.researchsquare.com/files/rs-5507011/v1/910ec38062d8e0eecfedba8d.png"},{"id":93996715,"identity":"f84baf1d-a118-4e18-b321-4ae3048e29b1","added_by":"auto","created_at":"2025-10-21 07:09:18","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":4902133,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-5507011/v1/1d880cbb-8d71-4f24-9cc4-c548e877d9f6.pdf"},{"id":70129652,"identity":"d2f4c225-b5ad-4404-ad56-dd3c438301ec","added_by":"auto","created_at":"2024-11-28 15:50:28","extension":"docx","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":2241445,"visible":true,"origin":"","legend":"","description":"","filename":"ExtendedDataFigures.docx","url":"https://assets-eu.researchsquare.com/files/rs-5507011/v1/f39fc5ec76f4c69d549f51b8.docx"},{"id":70129651,"identity":"79e2f04d-7401-4691-a87e-224b664658c1","added_by":"auto","created_at":"2024-11-28 15:50:28","extension":"docx","order_by":2,"title":"","display":"","copyAsset":false,"role":"supplement","size":3657086,"visible":true,"origin":"","legend":"","description":"","filename":"SupplementaryInformation.docx","url":"https://assets-eu.researchsquare.com/files/rs-5507011/v1/62842e54e0db49a2012c3e70.docx"}],"financialInterests":"There is \u003cb\u003eNO\u003c/b\u003e Competing Interest.","formattedTitle":"Thalamocortical Hodology to Personalize Electrical Stimulation for Focal Epilepsy","fulltext":[{"header":"Introduction","content":"\u003cp\u003eOver 50\u0026nbsp;million people worldwide suffer from epilepsy, a neurological disorder characterized by recurrent seizures, stemming from synchronized and excessive paroxysmal electrical activity of neurons\u003csup\u003e\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e,\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e\u003c/sup\u003e. In most of the cases (\u0026gt;\u0026thinsp;60%), said seizures originate and develop in restricted and focal cortical regions of the brain\u003csup\u003e\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e\u003c/sup\u003e, condition clinically referred to as focal epilepsy. Conversely, in generalized epilepsy, pathological brain activity stems from widespread areas of the brain, affecting both hemispheres\u003csup\u003e\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e,\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e\u003c/sup\u003e. Independently of the focal or generalized onset, about 30% of epilepsy patients live with intractable disabling seizures\u003csup\u003e\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e\u003c/sup\u003e, which cannot be controlled with medical therapy. The current standard of care for medically refractory epilepsy (MRE) relies on the resection of the seizure focus, i.e., the brain region associated with earliest electrophysiological changes during a seizure event\u003csup\u003e\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e\u003c/sup\u003e. However, if said cortical areas are not amenable for safe removal (or disconnection) because of overlap with eloquent cortices, or extensive seizure foci, or a diagnosis of generalized epilepsy, neuromodulation approaches such as thalamic electrical stimulation might be an alternative\u003csup\u003e\u003cspan additionalcitationids=\"CR8 CR9 CR10\" citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e\u003c/sup\u003e. Recent studies have propelled thalamic stimulation for epilepsy to the forefront of neuroscience, underscoring its emerging role as a pivotal therapeutic strategy in the field\u003csup\u003e\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eThe thalamus is a deep central structure crucial for its multifaceted functions and intricate organization. Indeed, through extensive connections with the cortex and other subcortical and peripheral structures, it is integral to a wide array of functions including motor and sensory processing, consciousness, and sleep\u003csup\u003e\u003cspan additionalcitationids=\"CR14 CR15\" citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e\u003c/sup\u003e. This diverse functionality is attributed to the thalamus's varied nuclear architecture and its intimate relationship with cortical areas\u003csup\u003e\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e\u003c/sup\u003e. Importantly, since the early 1950s the thalamus was demonstrated to synchronize widespread cortical activity, hence contributing to the generation of pathological epileptic discharges\u003csup\u003e\u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eIndeed, altered thalamocortical synchrony was extensively observed during focal seizures and associated loss of consciousness\u003csup\u003e\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e\u003c/sup\u003e. Furthermore, structural connectivity variations was reported between the thalamus and the cortex in patients with generalized epilepsy\u003csup\u003e\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e\u003c/sup\u003e .\u003c/p\u003e \u003cp\u003eGiven these cortico-subcortical dynamics and epileptogenic interactions, thalamic electrical stimulation has been proposed since the early 80s, with indications both in focal\u003csup\u003e\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e,\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e\u003c/sup\u003e and generalized epilepsy\u003csup\u003e\u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e\u003c/sup\u003e. In this regard, thalamic sites characterized by diffuse connectivity as the centromedian nucleus (CM), have been hypothesized to benefit generalized seizures, which originate and develop through widespread neuronal networks. Indeed, in a recent clinical trial, Cukiert and colleagues reported a significant seizure reduction (\u0026gt;\u0026thinsp;50%) with CM electrical stimulation in 90% of tested patients diagnosed with generalized epilepsy\u003csup\u003e\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e\u003c/sup\u003e. CM stimulation, instead, has been less effective in patients with focal MRE\u003csup\u003e\u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e23\u003c/span\u003e\u003c/sup\u003e. For this reason, thalamic nuclei with less diffuse connectivity, and mostly the anterior nucleus (ANT), have been explored as a possible target for neuromodulation in focal epilepsy. Already in the early 80s, Cooper et al. showed a significant clinical control of seizures in 4 out of 6 patients with stimulation of the ANT\u003csup\u003e\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e\u003c/sup\u003e. Since this seminal work, multiple studies have investigated the efficacy of ANT stimulation to treat focal epilepsy\u003csup\u003e\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e,\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e,\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e25\u003c/span\u003e\u003c/sup\u003e, including the pivotal SANTE trial that achieved 75% median seizure reduction at 7 years follow up\u003csup\u003e\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e\u003c/sup\u003e. Thanks to these works, ANT stimulation for focal epilepsy received FDA approval in 2018. However, clinical outcomes are highly variable depending on the anatomical location of the seizure onset. Specifically, the most well-documented efficacy was observed in focal onset seizures occurring in the temporal and frontal lobe (78% and 86% seizure frequency reduction), whereas effects were limited for other cortical seizure onset zones (SOZs) (39%)\u003csup\u003e26\u003c/sup\u003e. This might be due to the ANT\u0026rsquo;s unique anatomical connectivity to temporal and frontal areas\u003csup\u003e\u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e,\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eExploration in other areas of the thalamus further suggests a possible link between thalamocortical anatomical connectivity and neuromodulation efficacy. Indeed, the pulvinar nucleus (PUL) of the thalamus, which has high connectivity towards the temporal and occipito-parietal lobes and present different neural firing patterns during ictal progression and termination in patients with temporal lobe epilepsy\u003csup\u003e\u003cspan additionalcitationids=\"CR30 CR31\" citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e\u003c/sup\u003e, has been shown to be potentially effective for posterior quadrant epilepsy\u003csup\u003e\u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e,\u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e\u003c/sup\u003e. Similarly, the ventrolateral portion of the thalamus (ventral intermediate /ventral oral posterior nuclei, VIM/VOP)\u003csup\u003e35\u003c/sup\u003e and the cerebellum\u003csup\u003e\u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e36\u003c/span\u003e\u003c/sup\u003e, which present high structural connectivity to motor cortical areas, were tested, with positive outcome, for epilepsies organized in the agranular pre-central neocortical areas (rolandic epilepsy), which is particularly challenging for resective surgeries due to the high risks of motor deficits.\u003c/p\u003e \u003cp\u003eHowever, these studies remain isolated investigations leading to a lack of consensus on which target is optimal for focal epilepsy and clinical indication on how to approach the heterogeneity of SOZs, with limited clinical outcome. Based on prior evidence, here we argue that, for focal epilepsy, thalamic electrical stimulation should match the hodology of the thalamocortical pathway to enhance therapeutic outcome. In other words, targeting electrical stimulation to the thalamic subnucleus with preferential anatomical and functional connectivity to the SOZ (an approach that we refer to as \u003cem\u003ehodological matching\u003c/em\u003e) will improve the efficacy of thalamic stimulation for the treatment of refractory focal epilepsy.\u003c/p\u003e \u003cp\u003eTo test this hypothesis, we studied the anatomical and functional properties of neural activity of 32 patients with focal MRE. We focused on the study of three thalamic nuclei that exhibit more specific and less diffuse thalamocortical projections as compared to the CM: the ANT, the PUL, and the VIM/VOP. First, we defined a hodological map of the thalamocortical fibers from said thalamic nuclei towards the cortex using imaging and electrophysiology. In this way, we identified the cortical areas that match each nucleus, i.e., receiving major projections. We then characterized the functional involvement of these thalamic nuclei to ictal events occurring both in hodologically-matching and unmatching SOZs. Furthermore, we showed that electrical stimulation of only hodologically-matched thalamic nuclei immediately suppressed pathological discharges i.e., interictal spikes, in the epileptogenic cortex. Finally, we retrospectively verified that such effects led to significant clinical improvement by reducing seizure frequency in 7 patients with a chronic hodologically-matched thalamic implant. Importantly, our hodology-based approach achieved a median reduction in seizure frequency of 86.5%, which dramatically outperforms unmatched thalamic targeting reported by literature (39%).\u003c/p\u003e \u003cp\u003eIn summary, our approach, i.e., targeting electrical stimulation to match the hodology of thalamocortical connections, holds significant potential to improve the therapeutic efficacy of thalamic electrical stimulation in the of focal refractory epilepsy. This strategy introduces innovative avenues for clinical targeting and therapeutic intervention, while challenging the prevailing paradigm that adopts a \u0026ldquo;one-nucleus-fits-all\u0026rdquo; SOZs.\u003c/p\u003e"},{"header":"Results","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003eStudy rationale and design\u003c/h2\u003e \u003cp\u003eHere we aimed at investigating whether specific thalamic subnuclei might provide increased efficacy in the treatment of focal refractory epilepsy based on their hodological matching with the cortical SOZ. To demonstrate our hypothesis, we deployed a large array of electrophysiological and imaging assessments in a unique cohort of 32 focal MRE patients (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ea and \u003cb\u003eSupplementary Table\u0026nbsp;1\u003c/b\u003e for demographic and clinical information for each patient). Specifically, we analyzed intracranial electrophysiological data from 26 patients (S01-S26) undergoing (SEEG) exploration for epilepsy monitoring (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eb). Prior to SEEG implantation, we acquired preoperative structural magnetic resonance imaging (MRI) and postoperative computed tomography (CT) for precise SEEG contacts localization performed by a certified neurosurgeon (JGM) in all patients. All the participants had at least one SEEG electrode located in one of the nuclei of interest, namely the PUL, ANT or VIM/VOP nuclei, and a diagnosis of focal epilepsy. Overall, 37 thalamic nuclei among those of interest were sampled in this cohort (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eb, \u003cb\u003eSupplementary Table\u0026nbsp;2\u003c/b\u003e). We categorized each patient\u0026rsquo;s SOZ as frontal (meaning originating from frontal lobe regions such as the orbitofrontal cortex), temporal (including mesial and lateral temporal lobe structures), rolandic (including pre- and post-central peri-rolandic regions) or parietal/occipital cortex (namely originating from posterior-quadrant areas) (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ec, \u003cb\u003eSupplementary Table\u0026nbsp;3\u003c/b\u003e). In a subgroup of patients (n\u0026thinsp;=\u0026thinsp;8), we acquired and analyzed high-resolution diffusion MRI data to corroborate the anatomical organization of thalamocortical fibers (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ef). Such connectivity pattern was also confirmed through electrophysiology with thalamocortical evoked potentials (n\u0026thinsp;=\u0026thinsp;6). To study the interaction of the thalamus and the epileptic cortex, we then analyzed a total of 216 spontaneous seizures with concurrent sampling of the thalamic nuclei and cortical SOZ. Hence, each recorded seizure was categorized according to the SOZ anatomical location and the sampled thalamic nucleus during the event (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ec). We then performed electrophysiological testing in a subset of patients. Specifically, we used high-frequency thalamic stimulation (n\u0026thinsp;=\u0026thinsp;14) to test immediate effects of hodologically-matched thalamic stimulation on interictal epileptiform activity (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ef). Finally, we assessed long-term efficacy of our hodology based neuromodulation approach in 7 patients (S16, S27-32) with a chronic stimulation implant (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ed) and up to 3 years of follow-up (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ed-e).\u003c/p\u003e \u003cp\u003eImportantly, we also reported a single case of generalized epilepsy with simultaneous sampling of PUL, ANT, VIM/VOP as well as CM: this last site showed widespread structural and functional connectivity (\u003cb\u003eExtended Data\u003c/b\u003e Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e\u003cb\u003e)\u003c/b\u003e, replicating well known evidence from the literature\u003csup\u003e\u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e,\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e,\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e,\u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e43\u003c/span\u003e\u003c/sup\u003e and supporting its applicability as a target in generalized, rather than focal, epilepsy.\u003c/p\u003e \u003cp\u003eOverall, we designed and performed a comprehensive study to investigate the mechanisms and efficacy of hodologically-matched thalamic stimulation in focal refractory epilepsy.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e\n\u003ch3\u003eA specific hodological map of the thalamocortical pathway in the epileptic brain\u003c/h3\u003e\n\u003cp\u003eWhile previous studies proposed a structural connectivity map of the thalamocortical projections \u003csup\u003e\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e,\u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e,\u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e45\u003c/span\u003e\u003c/sup\u003e, the majority of these studies were performed in healthy brains and were limited to imaging techniques, such as fiber tracking. Here, we confirmed the presence of a well-defined anatomical map of the thalamocortical projections of PUL, ANT and VIM/VOP in the epileptic brain with two complementary methodologies: neuroimaging and electrophysiology.\u003c/p\u003e \u003cp\u003eFirst, we performed high-definition fiber tracking (HDFT) of high-resolution diffusion MRI data to confirm previously reported anatomical organization of cortical connectivity patterns of the PUL, ANT and VIM/VOP. We reconstructed all the likely axonal pathways between these nuclei and the cortical regions presenting seizures in our cohort (namely frontal, rolandic, parietal, occipital and temporal) (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ea). We quantified the relative strength of these connections by calculating the volume of thalamocortical projections from each nucleus to each cortical region normalized by the total volume of fibers (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eb). This analysis revealed a clear thalamocortical anatomical organization. Indeed, the PUL showed preferential connectivity towards the parietal and occipital regions, and the temporal lobe (mean volume value, frontal 0.11, rolandic: 0.04, \u003cem\u003eparietal\u003c/em\u003e: 0.38, \u003cem\u003eoccipital\u003c/em\u003e: 0.27, \u003cem\u003etemporal\u003c/em\u003e: 0.18); the ANT, instead, presented axonal projections going preferentially towards the frontal and temporal lobes (\u003cem\u003efrontal: 0.39\u003c/em\u003e, rolandic: 0.05, parietal: 0.06, occipital: 0.16, \u003cem\u003etemporal: 0.31\u003c/em\u003e); finally the VIM/VOP nuclei projected more remarkably to the rolandic cortex (frontal: 0.29, \u003cem\u003erolandic: 0.56\u003c/em\u003e, parietal: 0.03, occipital: 0.06, temporal: 0.03). This constricted structure of thalamocortical fibers from PUL, ANT and VIM/VOP contrasts with the broad anatomical connectivity of the CM (frontal: 0.27, rolandic: 0.17, parietal: 0.19, occipital: 0.19, temporal: 0.15), hence corroborating its indication as a more suitable target in generalized epilepsy (\u003cb\u003eExtended Data\u003c/b\u003e Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ea-b).\u003c/p\u003e \u003cp\u003eWe then confirmed this anatomical organization by analyzing thalamocortical evoked potentials (EPs) in the different cortical areas. For this, we applied low frequency (1 Hz) electrical stimulation to the different thalamic nuclei at 1\u0026ndash;3 mA while recording cortical neural response evoked by each pulse in the different cortical lobes (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ed, g, j). We then analyzed stimulation triggered averages of EPs. Following thalamic stimulation, we observed robust cortical EPs in anatomically connected (i.e., matching), but not in unmatching, cortical regions (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ee, h, k). Indeed, analysis of EPs peak to peak amplitudes showed significantly larger responses in those areas that matched the previously identified projections, further confirming a precise anatomical organization of distinct thalamic nuclei (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ec, f, i, l).\u003c/p\u003e \u003cp\u003eIn summary, we confirmed with multiple techniques a specific anatomical organization of the thalamocortical projections of PUL, ANT and VIM/VOP in the studied population. We then define hodologically-matching (parietal, occipital and temporal for PUL, frontal and temporal for ANT, and rolandic for VIM/VOP) and unmatching (frontal and rolandic for PUL, rolandic, parietal and occipital for ANT, frontal, temporal, parietal and occipital for VIM/VOP) regions.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e\n\u003ch3\u003eIncreased thalamocortical coupling during seizures in hodologically-matching SOZs\u003c/h3\u003e\n\u003cp\u003eGiven this specific structural organization of the thalamocortical fibers, we hypothesized a preferential functional involvement in the course of spontaneous seizures of the thalamic subnucleus that anatomically matches the SOZ (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003ed). To test this hypothesis, we analyzed a total of 216 spontaneous (i.e., not induced by cortical electrical stimulation) ictal events. The anatomical location of the SOZ as well as the relevant time markers (pre-ictal, beginning and end) of each seizure were determined by a certified clinical team. For each seizure we computed the well validated non-linear h2 correlation coefficient\u003csup\u003e\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e,\u003cspan additionalcitationids=\"CR47\" citationid=\"CR46\" class=\"CitationRef\"\u003e46\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e48\u003c/span\u003e\u003c/sup\u003e to identify time-resolved synchrony between each thalamic nucleus and each cortical SOZ. Specifically, we computed mean h2 coefficient during the ictal event (from seizure onset to termination) and the mean h2 coefficient during a 10 second baseline epoch spaced 2 minutes apart from the seizure\u0026rsquo;s onset\u003csup\u003e\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e\u003c/sup\u003e. We then compared mean h2 coefficient between the two phases separately for anatomically connected and non-connected thalamic nuclei to quantify the degree of specificity of thalamic involvement during spontaneous seizures.\u003c/p\u003e \u003cp\u003eAs hypothesized, we found that while anatomically matched nuclei showed a significant increase of synchrony between the thalamus and the SOZ, indicating functionally specific thalamic involvement during ictal events (\u003cb\u003eFigure a\u003c/b\u003e), minimal or no correlation was found in unmatched nuclei (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eb-c). This was consistent for all thalamic nuclei of interest and subjects (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003ee, mean h2 was 0.13 vs 0.20 in matched with p\u0026thinsp;\u0026lt;\u0026thinsp;0.001, and 0.12 vs 0.11 with p\u0026thinsp;\u0026gt;\u0026thinsp;0.05 in unmatched, for baseline and seizure respectively). Importantly, we also observed that thalamic nuclei had minimal or no correlation with non-SOZ regions (\u003cb\u003eExtended Data\u003c/b\u003e Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003ea-b), demonstrating the specificity of the correlation with the anatomically-matched SOZ. Additionally, for the patient diagnosed with generalized epilepsy we found that PUL, ANT, and VIM/VOP nuclei had a lower correlation with the epileptogenic zone (0.16, 0.25, 0.21 mean h2 across 24 seizures and 13 channels for PUL, ANT, and VIM/VOP, respectively) as compared to the CM (0.3, \u003cb\u003eExtended Data\u003c/b\u003e Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ec-d). This confirms the more diffused profile of the CM connectivity (both anatomically and functionally), entailing its suitability for generalized, rather than focal epilepsy treatment. Overall, these functional outcomes robustly validated the thalamocortical anatomical organization observed with fiber tracking and thalamocortical EPs, further supporting the proposed definition of hodologically-matching and unmatching nuclei.\u003c/p\u003e \u003cp\u003eTo further investigate possible mechanisms of thalamocortical interactions during focal seizures, we explored the coupling of the thalamus and SOZ in terms of seizure phases and frequency bands for matched thalamic nuclei. In particular, we aimed at verifying that the mechanism of thalamic involvement during seizures would not differ across the three nuclei of interest, since a different mechanism could imply a different design of electrical stimulation therapies. First, we asked whether maximal coupling of matched nuclei was specific to distinct seizure phases. To this aim, we computed the mean h2 correlation coefficient in three epochs: seizure initiation (including the 10 seconds preceding the first appearance of a tonic discharge in the SOZ contacts), middle part of the seizure (defined as the time interval separating seizure initiation and termination), and seizure termination (which includes the last 10 seconds of discharge). We found that the highest synchrony between the thalamus and SOZ occurred at seizure termination (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003ef). Importantly, this behavior was present for all three thalamic nuclei when hodologically (i.e., anatomically and functionally) matched to the SOZ (\u003cb\u003eExtended Data\u003c/b\u003e Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003ec). These findings support and validate the notion that synchronization loops between the thalamus and cortex might contribute to the termination of seizures, as suggested in previous studies in animals\u003csup\u003e\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e,\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e\u003c/sup\u003e and humans\u003csup\u003e\u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e50\u003c/span\u003e,\u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e51\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eAdditionally, we sought to determine if the observed synchrony presented a specific spectral profile. Hence, we compared the mean h2 coefficient computed in different frequency bands of the intracerebral recordings (delta, theta, alpha, beta, low gamma, high gamma, ripple, fast ripple). No spectral specificity was revealed: indeed, for all the frequency bands except the fast ripple band, the correlation between hodologically-matched thalamic nuclei and SOZ always increased significantly during the ictal event, while it showed no variation for unmatched thalamic nuclei (\u003cb\u003eExtended Data\u003c/b\u003e Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003ed). Hence, we asked whether specific frequency bands would better distinguish hodologically-matched and unmatched nuclei potentially suggesting biomarkers for the refinement of the stimulation design. For this, we compared the increase in ictal correlation respect to baseline between hodologically-matched and unmatched nuclei in different frequency bands. Interestingly, we found that theta, alpha and beta frequencies had a higher difference between the h2 correlation change during seizures between hodologically-matched and unmatched nuclei (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eg). These results show that the theta-beta frequencies might be more relevant to differentiate thalamic nuclei based on their connectivity to the SOZ.\u003c/p\u003e \u003cp\u003eAll in all, these findings demonstrate that thalamic nuclei exhibit enhanced neural synchrony with hodologically-matched cortical SOZs, particularly at the period of seizure termination, and primarily driven by theta, alpha, and beta oscillations.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e\n\u003ch3\u003eThalamocortical synchrony revers directionality during spontaneous seizures\u003c/h3\u003e\n\u003cp\u003eWhile the h2 index revealed an increase in correlation between the hodologically-matched thalamic nuclei and the SOZ throughout the seizure course, here we sought to determine the directionality of this coupling. Importantly, previous studies reported a change in the directional coupling of thalamocortical interactions in temporal lobe epilepsy throughout the seizure course\u003csup\u003e\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e,\u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e52\u003c/span\u003e\u003c/sup\u003e. Specifically, thalamic involvement was shown to occur early in a seizure, although it did not typically lead at the onset\u003csup\u003e\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e\u003c/sup\u003e. Additionally, animal studies have shown that the thalamus participates during seizures and plays a crucial role in their development, but it is less involved at seizure onset\u003csup\u003e\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e\u003c/sup\u003e. In order to explore whether a similar behavior was prevalent for all thalamic structures here investigated, we further expanded the analysis of hodologically-matched thalamocortical synchrony during spontaneous seizures by investigating the directionality of the coupling. To this purpose, we used Granger Causality (GC)\u003csup\u003e\u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e39\u003c/span\u003e\u003c/sup\u003e between the thalamic nuclei of PUL, ANT and VOP/VIM and the hodologically-matching SOZs. This methodology could indeed provide insights in the understanding of propagation patterns of epileptic activity helping the future design of stimulation patterns aimed at decreasing epileptiform activity. Since our previous results suggest a crucial role in the theta-alpha-beta oscillations, for this analysis we included only band-passed filtered (4\u0026ndash;30 Hz) seizures that showed a clear increase in h2 coefficient in this frequency range and lasted at least 10 seconds. This resulted in a total of 98 seizures. For each seizure, we extracted two 10 s epochs of interest at seizure initiation and seizure termination and computed the GC matrix for the SOZ contacts and the hodologically-matching thalamic nucleus contacts (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ea). Importantly, GC of unmatched nuclei was not explored because of the lack of correlation with SOZ. To test whether the role of different thalamic nuclei was changing over the course of the seizure, we compared the GC coefficient from the thalamus towards the SOZ (th\u0026loz;SOZ) and from the SOZ towards the thalamus (SOZ\u0026loz;th).\u003c/p\u003e \u003cp\u003eAt seizure initiation, in the majority of the seizures (70%), the GC coefficient of SOZ\u0026loz;th was higher than the GC coefficients th\u0026loz;SOZ (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eb for all subjects together, \u003cb\u003eExtended Data\u003c/b\u003e Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e for each individual) suggesting a leading role of the cortex. In contrast, more than half of the events (60%) presented a higher th\u0026loz;SOZ coefficient at the termination of ictal events hence indicating a change in the leading anatomical structure. This could be due to 1) SOZ\u0026loz;th decreasing while th\u0026loz;SOZ remaining stable (or decreasing less) throughout the course of the seizure, 2) th\u0026loz;SOZ increasing and SOZ\u0026loz;th remaining stable (or increasing less). To test for this, we computed the relative change of the GC coefficient between termination and initiation and compared it between SOZ\u0026loz;th and th\u0026loz;SOZ. Interestingly, we noticed an increase in both the thalamocortical (th\u0026loz;SOZ) and corticothalamic (SOZ\u0026loz;th) coupling throughout the duration of the event (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ec). However, the relative change was higher for the th\u0026loz;SOZ consistently across all the subjects, further demonstrating a leading role of the thalamus at seizure termination (73.41 vs 345.6%, for SOZ\u0026loz;th and th\u0026loz;SOZ, respectively) (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ec-d).\u003c/p\u003e \u003cp\u003eImportantly, we verified that the same characteristics of thalamocortical interactions were not specific to a particular nucleus. As a matter of fact, in all studied nuclei, the absolute GC value in the initiation epoch and termination epoch showed a reverse trend, with a primary leading role of the cortex at first (85% for PUL, 59% for ANT, 60% for VIM/VOP), followed by a predominantly leading role of the thalamus at termination (63% for PUL, 71% for ANT, 52% for VIM/VOP) (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ee). Additionally, PUL, ANT and VIM/VOP showed a significantly higher relative change for th\u0026loz;SOZ throughout the seizure (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ef).\u003c/p\u003e \u003cp\u003eIn summary, here we investigated the underlying mechanisms of thalamocortical and corticothalamic coupling involving the ANT, PUL and VIM/VOP with hodologically-matched SOZ during spontaneous seizures. Our findings suggest a potential trend where cortical regions exert a dominant influence at seizure initiation, while the thalamus may assume a more prominent role as the ictal events progresses. This pattern was observed across all three nuclei, suggesting that the development of nuclei-specific stimulation protocols -particularly regarding the timing of stimulation during the ictal event- may not be necessary to address the varying underlying mechanisms of thalamic action.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e\n\u003ch3\u003eMatched thalamic stimulation immediately suppresses epileptiform discharges\u003c/h3\u003e\n\u003cp\u003eGiven the peculiar hodology that we previously found both in the structural connectivity of the PUL, ANT and VIM/VOP, as well as in the functional coupling observed during ictal events, we hypothesized that targeting electrical stimulation to the thalamic nucleus that hodologically-matches the cortical SOZ might be more effective in reducing epileptiform activity than unmatched nuclei. To test this hypothesis, we delivered bipolar thalamic stimulation through the SEEG electrodes in 14 patients. For clinical reasons, the stimulation testing was always performed after the patient had shown spontaneous seizures. We evaluated the effect of hodologically-matched and unmatched stimulation on interictal epileptiform discharges (IEDs), namely interictal spikes. We used an automatic validated algorithm\u003csup\u003e\u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e41\u003c/span\u003e,\u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e42\u003c/span\u003e\u003c/sup\u003e to extract IEDs from cortical recordings in the SOZ (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003ea), previously subdivided in 5 seconds epochs and denoised (see \u003cb\u003eMethods\u003c/b\u003e). We aimed at quantifying the difference between baseline (including the epochs in the 2 minutes preceding the first stimulation applied) and stimulation epochs (see \u003cb\u003eMethods\u003c/b\u003e). For this, we performed two distinct analysis: a) in n\u0026thinsp;=\u0026thinsp;8 patients with continuous evident pathological spiking (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003ea), we computed the rate of IEDs per minute in baseline (Stim OFF) and stimulation (Stim ON) epochs; b) in n\u0026thinsp;=\u0026thinsp;6 patients with more than half of the baseline epochs without pathological spikes, we computed the probability of suppressed IEDs in stimulation epochs and compared it to the probability of suppressed IEDs in baseline epochs, hence obtaining a percentage change in IED suppression probability (\u003cb\u003eSupplementary Table\u0026nbsp;2\u003c/b\u003e). In the first analysis, a positive effect of the stimulation (meaning reducing pathological IEDs) should result in a smaller IEDs rate during stimulation epochs than during baseline (stimulation off) epochs; whereas in the latter, a positive number indicated that the stimulation epochs are more likely to present no IEDs than the baseline epochs (without stimulation).\u003c/p\u003e \u003cp\u003eAs hypothesized, we found that electrical stimulation of the PUL, ANT and VIM/VOP significantly suppressed cortical IEDs when the SOZ was in a cortical region that hodologically-matched the thalamic nucleus stimulated (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eb-d, g). Overall, our analysis revealed that hodologically-matched but not unmatched thalamic stimulation significantly reduced the IED rate in all tested patients (for matched stimulation, IED variation to baseline ranged from \u0026minus;\u0026thinsp;29.12% to -73% with p\u0026thinsp;\u0026lt;\u0026thinsp;0.01 and p\u0026thinsp;\u0026lt;\u0026thinsp;0.001, while for unmatched stimulation, it ranged from \u0026minus;\u0026thinsp;4.2% to -25%) (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eg). Here, we highlight the IED rate values for a representative patient for each nucleus of interest: we observed a reduction of the IED rate of 132 vs 78 spikes/minute for the PUL (median rate at stimulation off vs stimulation on in S01, p\u0026thinsp;\u0026lt;\u0026thinsp;0.001) (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eb), 144 vs 60 spikes/minute for the ANT (in S24, p\u0026thinsp;\u0026lt;\u0026thinsp;0.001) (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003ec), and 132 vs 48 spikes/minute in the VIM/VOP (in S13, p\u0026thinsp;\u0026lt;\u0026thinsp;0.001) (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003ed). Similarly, when evaluating the percentage change in the probability of IED suppression, we found that hodologically-matched thalamic stimulation epochs were more likely to present no IEDs when compared to unmatched thalamic stimulation epochs. Indeed, In 5 out of 6 patients that received stimulation in the nucleus matching their SOZ, we found at least a 20% increase (up to over 200%) in the probability of epochs with no IEDs with stimulation on with respect to stimulation off (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eh); whereas the stimulation effect was more inconsistent (-17 to +\u0026thinsp;7%) when the targeted nucleus was not matching the hodology of the thalamocortical fibers towards the SOZ (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eh). This finding further suggests that the more effective thalamic target to suppress IEDs is the one matching the hodology of the SOZ.\u003c/p\u003e \u003cp\u003eOverall, we showed that only targeted electrical stimulation of the nuclei hodologically-matching the SOZ is significantly reducing the rate of IEDs and increasing the probability that stimulation epochs present no IEDs (complete suppression). These findings are crucial in guiding approaches to address the heterogeneity of SOZs and in clinically selecting optimal stimulation targets for the electrical neuromodulation of focal epilepsies.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cdiv id=\"Sec8\" class=\"Section2\"\u003e \u003ch2\u003eMatched thalamic electrical stimulation in chronic epilepsy treatment\u003c/h2\u003e \u003cp\u003eFinally, to verify whether the observed immediate effects on IEDs could be translated in a long-term clinical improvement in epileptic patients, we evaluated 7 patients (one of which being part of the SEEG cohort) that were chronically implanted for thalamic stimulation in a hodologically-matching thalamic nucleus, as part of their clinical care. Across all patients, 6 (S16, S27-S31) presented focal MRE originating from the posterior quadrant regions and hence received a neurostimulation system in the PUL, while one (S32) was implanted in the ANT to treat a frontal onset MRE (\u003cb\u003eSupplementary Table\u0026nbsp;4\u003c/b\u003e). Overall, we observed a median reduction in seizure frequency of 86.5% (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003ea, b) at 15.86\u0026thinsp;\u0026plusmn;\u0026thinsp;7.9 (mean\u0026thinsp;\u0026plusmn;\u0026thinsp;SD) months of follow-up. As a representative example, for S16 seizures reduced from 5 per week to 1 per month, and the patient also self-reported a significant reduction in seizure severity, noting no missed workdays since the implant of the stimulation device, compared to missing 1\u0026ndash;2 days per month prior to surgery. Additionally, the anti-seizure medication (ASM) intake for each patient decreased for 5/7 patients since the matched neurostimulation system was implanted, and did not increase for any of the participants (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003ec). Importantly, these results dramatically outperform the unmatched thalamic stimulation results reported in literature for the SANTE trial at 7 years follow-up. Indeed, in this previous study, unmatched stimulation (ANT stimulation for seizure onsets other than frontal or temporal) only achieved 39% reduction in seizure frequency, whereas outcome for matched SOZ were comparable to the ones we observed (78% and 86% for temporal and frontal onset, respectively)\u003csup\u003e\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e\u003c/sup\u003e.\u003c/p\u003e \u003cp\u003eIn conclusion, our results indicate that hodologically-matched thalamic stimulation provides significant therapeutic benefits, including notable reductions in the frequency and severity of disabling seizures. This approach not only outperforms current standard treatments but also lays the groundwork for more personalized and precise target selection in neuromodulation for medically refractory focal epilepsy, particularly in patients who are not candidates for curative resective or ablative interventions.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e"},{"header":"Discussion","content":"\u003cp\u003eElectrical stimulation of the thalamus is undergoing a rapid advance as a treatment opportunity for patients with MRE\u003csup\u003e\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e,\u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e53\u003c/span\u003e,\u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e54\u003c/span\u003e\u003c/sup\u003e, with indications in both generalized\u003csup\u003e\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e,\u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e53\u003c/span\u003e,\u003cspan citationid=\"CR55\" class=\"CitationRef\"\u003e55\u003c/span\u003e\u003c/sup\u003e and focal epilepsy\u003csup\u003e\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e,\u003cspan additionalcitationids=\"CR57 CR58 CR59\" citationid=\"CR56\" class=\"CitationRef\"\u003e56\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR60\" class=\"CitationRef\"\u003e60\u003c/span\u003e\u003c/sup\u003e. While in generalized epilepsy the CM nucleus effectively reduce seizure events, in focal epilepsy, this expansion requires a profound understanding of the intricate anatomical and functional thalamocortical interactions with highly heterogeneous and individualized epilepsy phenotypes to improve clinical outcome. Here, we proposed a novel approach to select the optimal targets for thalamic electrical stimulation in focal epilepsy, considering the specificity of SOZs in relation to cortical and subcortical anatomy. Specifically, we demonstrated that leveraging the hodology-based interactions between certain subnuclei of the thalamus and the epileptogenic cortex results in stronger neuromodulator efficacy. This approach highlights the importance of thalamocortical hodology in individualizing stimulation targets, potentially advancing therapeutic outcomes, and paving the way for more effective neuromodulation strategies in patients with focal MRE who are not candidates for curative interventions.\u003c/p\u003e \u003cp\u003eSince the early stages of epilepsy surgery and the first mechanistic studies about epileptogenicity, focus has increasingly turned towards deeper brain structures, particularly the thalamus\u003csup\u003e\u003cspan citationid=\"CR61\" class=\"CitationRef\"\u003e61\u003c/span\u003e,\u003cspan citationid=\"CR62\" class=\"CitationRef\"\u003e62\u003c/span\u003e\u003c/sup\u003e. Early animal studies by Penfield and Jasper in 1947 identified the thalamus as a potential generator of epileptic cortical activity, proposing the theory of \"centrencephalic\" epilepsy\u003csup\u003e\u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e,\u003cspan citationid=\"CR63\" class=\"CitationRef\"\u003e63\u003c/span\u003e\u003c/sup\u003e. Indeed, this structure mediates reciprocal cortical-subcortical connections, functioning as an 'integrative hub' for functional brain networks\u003csup\u003e\u003cspan citationid=\"CR64\" class=\"CitationRef\"\u003e64\u003c/span\u003e\u003c/sup\u003e. Recently (challenging previous theories), the thalamus has been implicated not in the generation, but in the propagation of generalized and focal onset seizures, due to its role in the bi-hemispheric cortical spread of epileptiform activity\u003csup\u003e\u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e,\u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e50\u003c/span\u003e,\u003cspan citationid=\"CR65\" class=\"CitationRef\"\u003e65\u003c/span\u003e\u003c/sup\u003e. In addition, more recent evidence indicates that the thalamus plays a critical role also in the network dynamics of focal-onset epilepsies\u003csup\u003e\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e,\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e,\u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e\u003c/sup\u003e, but lacked a consistent and systematic analysis of the hodological patterns of propagation from different cortical and thalamic regions during ictal events. Our results contribute to this discussion. Indeed, we detected significant correlation with ictal activity in all recorded thalamic nuclei further suggesting the thalamus as a potential neuromodulation target to treat focal MRE. However, correlation between thalamic and cortical SOZ activation was specific to thalamocortical anatomical connection, as identified by high-density tractography as well as electrophysiological assessment, supporting the hypothesis that epileptic propagation patterns during seizures obey predictable anatomo-functional patterns and highlighting the need to carefully select thalamic nuclei for neuromodulation according to individual clinical and electrophysiological data.\u003c/p\u003e \u003cp\u003eWe applied this mechanistic understanding of the thalamocortical interactions to inform hodology-based electrical stimulation aimed at suppressing epileptiform activity. Our results demonstrated that hodologically-matched thalamic stimulation immediately suppresses interictal epileptiform discharges. Additionally, this targeted stimulation markedly enhances the likelihood of achieving complete suppression of IEDs during stimulation epochs, suggesting robust potential for targeted thalamic stimulation as a therapeutic strategy in epilepsy neuromodulatory management. Immediate effects in n\u0026thinsp;=\u0026thinsp;14 subjects, paralleled with the drastic clinical improvement in seizure frequency with chronic matched thalamic stimulation (up to 95%) in seven patients, provide promising evidence that a hodology-based stimulation would improve clinical outcome for patients with focal epilepsy. Importantly, the comparative results show a substantial improvement over existing literature on unmatched thalamic neuromodulation (86.5% versus 39% reduction in seizure frequency), highlighting the potential need to reconsider current clinical practices, where a single target is often applied to all seizure onset zones.\u003c/p\u003e \u003cp\u003eMechanistically, our results identified a trend suggesting that the epileptic cortex is the dominant structure at the initiation of seizures, while the thalamus assumes a leading role throughout the course of ictal events, showing higher synchronization towards the middle and termination of the ictal events. Importantly, this pattern was consistent across all thalamic regions investigated. In this regard, in an animal model of limbic epilepsy induced by GABA antagonist (i.e., bicuculline), Aracri and colleagues demonstrated hippocampal activation first, progressing to hippocampus-thalamus synchronization, and then thalamic bursts leading to the termination of seizures\u003csup\u003e\u003cspan citationid=\"CR66\" class=\"CitationRef\"\u003e66\u003c/span\u003e\u003c/sup\u003e. Destexhe and colleagues also highlighted the role of cortically induced coherence of thalamic-generated oscillations\u003csup\u003e\u003cspan citationid=\"CR67\" class=\"CitationRef\"\u003e67\u003c/span\u003e\u003c/sup\u003e. They implicated temporal and spatial thalamocortical interactions, suggesting that, at the onset of seizures, thalamocortical cells are hyperpolarized due to the greater power of excitatory cortical projections to reticular neurons compared to thalamic relay neurons. Once the initial cortical fast phase is abolished, thalamic relay neurons projecting to the cortex depolarize synchronically, reinforcing the coherence of seizure activity arising from the cortex, and promoting the highly synchronous cortical spike phase observed during the middle and end of seizures. Our findings related to the directionality of epileptic activity through thalamocortical circuits align with this mechanistic theory, where information flows to the thalamus from the epileptic cortex at seizure initiation and back to cortical areas during the late seizure phase. This supports the hypothesis that the thalamus plays distinct roles in this dynamic process\u003csup\u003e\u003cspan citationid=\"CR68\" class=\"CitationRef\"\u003e68\u003c/span\u003e\u003c/sup\u003e. Furthermore, it suggests the possibility to design closed-loop stimulation paradigm to restrict the stimulation to precise moment of the seizure development with potentially stronger neuromodulatory effects. In this regard, our findings indicate that the theta to beta frequency range represents a key hallmark of thalamocortical connectivity\u003csup\u003e\u003cspan citationid=\"CR69\" class=\"CitationRef\"\u003e69\u003c/span\u003e\u003c/sup\u003e, which could serve as a biomarker for seizure detection and closed-loop monitoring.\u003c/p\u003e \u003cp\u003eThe primary limitation of this study is that the clinical assessment of matched electrical stimulation as a chronic therapy was evaluated in just seven patients as an observational retrospective report. Hence, inhomogeneous evaluation of all thalamic nuclei and SOZ occurred. However, these results are paralleled to a thorough investigation of the mechanisms of thalamocortical coupling in focal epilepsy, which support our findings. Additionally, even if limited in number, our evaluation of the effect of matched stimulation in a chronic stage of the treatment results in a drastic clinical improvement well above other literature reports. Larger clinical studies are now necessary to confirm the efficacy of hodologically-matched stimulation in the reduction of seizure frequency and improvements in patients\u0026rsquo; quality of life in chronic settings. Another study limitation includes sampling bias associated with explorations of the human brain using intracranial electrodes, where conclusions can only be drawn from areas where electrodes were implanted. Nevertheless, strict inclusion and exclusion criteria minimized this limitation by focusing on subjects with highly localized seizures, many of whom became seizure-free after restricted focal resections. Additionally, the applied parameters of acute stimulation were restricted to few parameters (100 Hz only, 1\u0026ndash;3 mA, continuous stimulation), mainly because the focus of the current study was related to the thalamic targeting and not to the different paradigms of electrical stimulation. However, several studies previously identified the optimal stimulation parameters in the high frequency range (\u0026gt;\u0026thinsp;100Hz)\u003csup\u003e70\u0026ndash;72\u003c/sup\u003e, as in our chronic participants. Future studies will need to better refine stimulation parameters and identify optimal stimulation frequencies.\u003c/p\u003e \u003cp\u003eIn summary, understanding the thalamocortical hodological principles governing the organization of focal seizures and their potential clinical applications is invaluable. Personalizing therapies based on the connectivity of the seizure onset zone to thalamic sub-nuclei and individual seizure characteristics may lead to improved efficacy and more favorable outcomes in the neuromodulation field for epilepsy. Our work addresses the critical unmet clinical need to identify optimal neuromodulation target by demonstrating predictable patterns of connectivity in thalamocortical interactions and guiding optimal thalamic targeting based on location of ictogenesis. This hodology-based approach would impact the current clinical standard in which one target fits all the SOZ (ANT in focal epilepsy, CM in generalized epilepsy), and lay the foundations for highly individualized neuromodulation treatments.\u003c/p\u003e"},{"header":"Methods","content":"\u003cdiv id=\"Sec11\" class=\"Section2\"\u003e \u003ch2\u003eTrial and participants information\u003c/h2\u003e \u003cp\u003e All experimental protocols were approved by the University of Pittsburgh Institutional Review Boards (IRB) (protocol STUDY20070113 and STUDY21020058). All participants included in this study were diagnosed with drug-refractory epilepsy and underwent stereoelectroencephalography (SEEG) implantation at the University of Pittsburgh Medical Center as a part of their standard clinical care. They were hospitalized in the epilepsy monitoring unit (EMU) from January 2021 to May 2024, for 5.04\u0026thinsp;\u0026plusmn;\u0026thinsp;6.28 days (mean\u0026thinsp;\u0026plusmn;\u0026thinsp;SD). All patients received a comprehensive neurological assessment, neuropsychological testing, routine MRI and CT and SEEG implantation. Informed consent was obtained from all patients or legal representatives. In order to study the anatomical and functional properties of the thalamus and the epileptic cortex via SEEG monitoring, we retrospectively included patients (S01-S26) in the present study if: 1) at least one SEEG electrode contact was in the PUL, ANT or VIM/VOP, and 2) at least one of the electrophysiological recordings analyzed in the study was performed (minimum 2 spontaneous seizures, thalamocortical evoked potentials, HDFT, thalamic stimulation). Exclusion criteria encompassed refusal in participating in the study and SEEG implantation complication (severe intracranial hemorrhage). We additionally included 7 participants (S16, S27-S32) that received a chronic thalamic stimulation system as part of their clinical care to evaluate long-term efficacy of matched thalamic stimulation.\u003c/p\u003e \u003cp\u003eA comprehensive overview of the patients\u0026rsquo; demographics and epilepsy categorization, summary of the experiments, and a summary of the data collected for each patient are reported in \u003cb\u003eSupplementary Table\u0026nbsp;1\u0026amp;2\u003c/b\u003e. Overall, we included 26 SEEG participants (20 males, 6 females) and 7 chronically implanted patients (S16, S27-S32), of age 34.8\u0026thinsp;\u0026plusmn;\u0026thinsp;11.3 (mean\u0026thinsp;\u0026plusmn;\u0026thinsp;SD). The age at seizure onset was 14.3\u0026thinsp;\u0026plusmn;\u0026thinsp;13.6 years-old, while the duration of the epilepsy disease was 20.4\u0026thinsp;\u0026plusmn;\u0026thinsp;13.1 years (mean\u0026thinsp;\u0026plusmn;\u0026thinsp;SD). Additional information can be found in the \u003cb\u003eSupplementary Tables\u0026nbsp;1, 2, 4\u003c/b\u003e.\u003c/p\u003e \u003cp\u003e \u003cb\u003eSEEG implantation and localization.\u003c/b\u003e \u003c/p\u003e \u003cp\u003eThe number and location of SEEG electrodes implanted was pre-operatively planned individually for each patient as part of their clinical care and was not influenced by this study. The implant procedure was performed with robotic stereotactic guidance (ROSA, Zimmer-Biomet, Warsaw, IN, USA), applying bone fiducial registration with accuracy\u0026thinsp;\u0026lt;\u0026thinsp;0.5 mm). Overall, we implanted a total of 5214 contacts (Microdeep\u0026reg; SEEG Electrodes, 12 to 18 channels, DIXI Medical, Marchaux-Chaudefontaine, France), with an average per patient of 200.53\u0026thinsp;\u0026plusmn;\u0026thinsp;27.87. To target the thalamic nuclei of interest, we used the following stereotactic coordinates: 1) for PUL, in relation to the AC/PC reformatted planes, and having the PC (posterior commissure) as the reference point, coordinates are X: 2\u0026ndash;10 mm lateral to the midline plane; Y: 0 to -6 mm posterior to PC; Z: 0 to -6 mm below the AC/PC defined horizontal plane; 2) for ANT, in relation to the AC/PC reformatted planes, and having the MC (mid commissural point) as the reference point, coordinates are X: 6\u0026ndash;8 mm lateral to the midline plane; Y: 6 to 7 mm anterior to the MC; Z: 2 to 4 mm above the AC/PC defined horizontal plane; 3) for VIM/VOP, lateral (from about 5 to about 15 mm lateral to the AC/PC line); anterior/posterior (from about 2 to about 10 mm anterior to PC); and dorsal/ventral (from about\u0026thinsp;+\u0026thinsp;1 to about \u0026minus;\u0026thinsp;2 mm from the AC/PC plane.\u003c/p\u003e \u003cp\u003eTo reconstruct the position of the SEEG electrode, we used CURRY software (Compumedics NeuroScan, Hamburg, Germany). First, we marked key reference points on each pre-implantation MRI T1 scan (namely the AC, PC and the nasion). These images were then co-registered with post-implantation CT, and a clinical team member manually identified and labeled all the SEEG contacts from the CT. All the contact locations of the subject MRI were visually inspected and confirmed by a certified neurosurgeon (JGM).\u003c/p\u003e \u003cp\u003eFor all subject-analysis (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eb), the electrode coordinates in the subject space MRI were then translated to coordinates in MNI space for each patient. To confirm the contact location within the thalamic nuclei of interest we used The Human Motor Thalamus Atlas\u003csup\u003e\u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e37\u003c/span\u003e\u003c/sup\u003e. The assessment of SOZ contacts (as well as non-SOZ contacts used in the following analysis) was confirmed by an experienced neurosurgeon (JGM) (\u003cb\u003eSupplementary Table\u0026nbsp;3\u003c/b\u003e).\u003c/p\u003e \u003cp\u003e \u003cb\u003eSEEG recording and stimulation.\u003c/b\u003e \u003c/p\u003e \u003cp\u003eSEEG data were recorded with Natus Quantum System EEG diagnostic and monitoring system (Natus, Pleasanton, A, USA), with a sampling rate of 2048 Hz (for S26 the sampling rate was 1024 Hz) during the extra-operative monitoring at the Epilepsy Monitoring Unit of the University of Pittsburgh Medical Center (Presbyterian Hospital). For all analysis, we applied SEEG bipolar montage to increase spatial selectivity and reduce noise levels. To deliver electrical stimulation in the thalamic contacts, we used the Nicolet Cortical Stimulator (Natus, Pleasanton, A, USA), which delivers biphasic pulses up to 100 Hz (amplitude was set to 1\u0026ndash;3 mA and pulse duration to 60\u0026ndash;300 \u0026micro;s for all patients) with bipolar electrode configuration.\u003c/p\u003e \u003cp\u003e \u003cb\u003eHigh-Definition Fiber Tracking.\u003c/b\u003e \u003c/p\u003e \u003cp\u003eTo estimate anatomical projections from the thalamus to various cortical areas, we first performed High-Definition Fiber Tracking (HDFT) of diffusion MRI data. The diffusion images were acquired on a SIEMENS Prisma Fit scanner using a diffusion sequence (2mm isotropic resolution, TE/TR\u0026thinsp;=\u0026thinsp;99.2 ms/2490 ms, 257 diffusion sampling with maximum b-value 4010 s/mm\u0026sup2;). Diffusion tensor estimation and tractography were performed using DSI studio (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttp://dsi,studio.labsolver.org\u003c/span\u003e\u003cspan address=\"http://dsi,studio.labsolver.org\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e). The accuracy of b-table orientation was examined by comparing fiber orientations with those of a population-averaged template. The tensor metrics were calculated using DWI with b-value lower than 1750 s/mm\u0026sup2;. For fiber tracking, we used a tracking threshold of 0, angular threshold of 0, and a step size of 0 mm. We utilized seed regions in thalamus to create white matter tracts to regions of interest in the cortex. Seed regions were selected in ANT, VIM/VOP, PUL, and CM based on the extended Human Connectome Project multimodal parcellation atlas. Regions of interest were selected in the frontal lobe, rolandic area, parietal lobe, occipital lobe, and temporal lobe. Tracks with lengths shorter than 30 mm or longer than 1000 mm were discarded. A total of 10,000 tracts were placed. Topology informed pruning was applied to the tractography with 2 interactions to remove false connections. We then quantified the volume of white matter tracts projecting from each thalamic nucleus to each cortical area. We normalized the volume of each white matter projection by the total volume projections from each thalamic nucleus (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ea-b, \u003cb\u003eExtended Data\u003c/b\u003e Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ea-b).\u003c/p\u003e \u003cp\u003e \u003cb\u003eThalamocortical evoked potentials.\u003c/b\u003e \u003c/p\u003e \u003cp\u003eIn order to confirm the anatomical organization of thalamocortical fibers with electrophysiological techniques, thalamocortical evoked potentials were elicited by thalamic stimulation of PUL, ANT and VIM/VOP (frequency\u0026thinsp;=\u0026thinsp;1 Hz, amplitude\u0026thinsp;=\u0026thinsp;1\u0026ndash;3 mA, pulse width\u0026thinsp;=\u0026thinsp;300 \u0026micro;s) and recorded from all available SEEG channels with Natus recording system. SEEG intracerebral data from all cortical channels were filtered with a band-pass Butterworth 2nd order filter (1-1000 Hz) and DC offset was removed from all recordings. From each stimulation pulse, we extracted 520 ms epochs (20 ms prior and 500 ms after the stimulus) and computed stimulation triggered averages. Each epoch was baseline corrected with the pre-stimulus interval. We calculated the peak-to-peak amplitude of evoked potentials in each epoch as the difference between the maximum and minimum voltage value observed within the first 400 ms from the stimulus. We performed this analysis for all available channels that presented sufficient signal to noise ratio, that were located in the regions of interest for this study (not in thalamus) and that were not located in the white matter (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ed-l, \u003cb\u003eExtended Data\u003c/b\u003e Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003cb\u003eAcquisition and analysis of spontaneous seizures.\u003c/b\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec12\" class=\"Section2\"\u003e \u003ch2\u003eSpontaneous seizures acquisition\u003c/h2\u003e \u003cp\u003eDuring the extra-operative monitoring, we recorded 216 spontaneous seizures from 23 patients. The remaining 3 (S14, S16, S19) patients either presented technical problems that invalidated the inclusion of the seizures or did not present any spontaneous events during the hospitalization. For our analysis, we collected recordings from 2 minutes before the seizure onset up to 30 seconds after the seizure termination. These time events (onset and termination) were manually marked for each event by certified medical personnel. We used such markers to define 5 epochs of interested, used in subsequent analysis: 1) baseline (corresponding to 20 seconds occurring two minutes before the seizure onset); 2) seizure initiation (including the 10 seconds preceding the first appearance of a tonic discharge in the SOZ contacts); 3) seizure termination (including the last 10 seconds of discharge); 4) middle of seizure (as the epoch between the previous two intervals 5) total seizure duration (defined as the interval between onset and termination).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec13\" class=\"Section2\"\u003e \u003ch2\u003eNon-linear correlation analysis\u003c/h2\u003e \u003cp\u003eFor each recorded seizure, we computed the h2 non-linear correlation coefficient, as implemented in Anywave\u003csup\u003e\u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e38\u003c/span\u003e\u003c/sup\u003e for each pair of thalamus-SOZ bipolar contacts. The h2 coefficient between two signals is a time-resolved measure of their non-linear dependence. The coefficient was computed over a 2 s time window sliding by steps of 1 s. We computed the h2 across the full SEEG spectrum (1-500 Hz), as well as separate frequency bands: delta (1\u0026ndash;4 Hz), theta (4\u0026ndash;8 Hz), alpha (8\u0026ndash;15 Hz), beta (15\u0026ndash;30 Hz), gamma (30\u0026ndash;45 Hz), gamma2 (55\u0026ndash;90 Hz), ripple (80\u0026ndash;250 Hz), fast ripple (250\u0026ndash;500 Hz). For each time interval of interest (seizure duration, seizure initiation, middle of seizure, seizure termination), we computed the mean h2 (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eb, e, f).\u003c/p\u003e \u003cp\u003eTo computer the change in h2 between matched and unmatched nuclei (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eg), we calculated, for each matched and unmatched pair, the percentage increase of h2 in the seizure duration with respect to the baseline epoch, as:\u003cdiv id=\"Equa\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equa\" name=\"EquationSource\"\u003e\n$$\\:h2\\:change\\left(\\%\\right)=100*\\left(mean\\right(h{2}_{seizure})-mean(h{2}_{baseline}\\left)\\right)/mean\\left(h{2}_{baseline}\\right)\\:\\%$$\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003eWe then performed the subtraction between the percentage h2 changes calculated in matched thalamus-SOZ and unmatched thalamus-SOZ pairs.\u003c/p\u003e \u003cp\u003eFor the correlation between the thalamus and non-SOZ regions, we randomly selected two seizures for each patient. We computed the mean h2 of the thalamus with 3 distinct non-SOZ pair of contacts (\u003cb\u003eSupplementary Table\u0026nbsp;3\u003c/b\u003e).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec14\" class=\"Section2\"\u003e \u003ch2\u003eGranger Causality analysis\u003c/h2\u003e \u003cp\u003eTo analyze the directionality of the interactions between the thalamus and the cortex, we used the Granger Causality method\u003csup\u003e\u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e39\u003c/span\u003e\u003c/sup\u003e. First, we selected all the seizures that lasted at least 10 s and showed an increase in the h2 correlation coefficient in the theta to beta band. We filtered the raw SEEG signals in the 4\u0026ndash;30 Hz range and divided each seizure recording in 10 s non-overlapping windows. We computed the GC matrix of each matched thalamus-SOZ pair for every epoch. This computation results in a 4X4 matrix for each epoch, where the second diagonal contains the thalamocortical (th\u0026rarr;SOZ) and corticothalamic (SOZ\u0026rarr;th,) GC value. In addition to the absolute GC value, we quantified per percentage change in GC from initiation to termination (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003ec, d, f) as:\u003cdiv id=\"Equb\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equb\" name=\"EquationSource\"\u003e\n$$\\:\\varDelta\\:GC\\:\\left(\\%\\right)=100*(G{C}_{termination}-G{C}_{initiation})/G{C}_{initiation})\\:\\%$$\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003ewhere termination and initiation are the 10 seconds epochs preceding the seizure onset and the seizure termination time stamps, respectively. To avoid misinterpretations given by the variability in the number of seizures recorded from each patient, we highlight the single-subject results in \u003cb\u003eExtended Data\u003c/b\u003e Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e. In a total of 20 patients, 17 showed an increase in the proportion of seizures with the thalamus leading at termination (pie charts in \u003cb\u003eExtended Data\u003c/b\u003e Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e) and 18 showed higher \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:\\varDelta\\:GC\\:\\)\u003c/span\u003e\u003c/span\u003ein the thalamocortical (and not corticothalamic) coefficient (bar plot in \u003cb\u003eExtended Data\u003c/b\u003e Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e). Such finding robustly confirms population-observed trends.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec15\" class=\"Section2\"\u003e \u003ch2\u003eCase report: CM connectivity in generalized epilepsy\u003c/h2\u003e \u003cp\u003eWe had the unique opportunity to justify the exclusion of the CM in our study (\u003cb\u003eExtended Data\u003c/b\u003e Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e), focused on focal epilepsy, by performing structural and functional connectivity analysis of the CM, as well as the PUL, ANT and VIM/VOP on a patient affected by generalized epilepsy. In particular, the patient was a female in her late 20s, with a severe medically refractory epilepsy, with more than 10 generalized seizures per day. The patient underwent SEEG exploration to investigate the appropriate thalamic target for neuromodulation as a part of her clinical care.\u003c/p\u003e \u003cp\u003eWe performed h2 correlation analysis of 24 seizures recorded with simultaneous sampling of CM, PUL, ANT, VIM/VOP as well as 13 bipolar electrodes sampling the following brain structures: insula, frontal operculum, superior temporal sulcus, superior temporal gyrus, middle temporal gyrus. We computed the mean h2 coefficient for each thalamic nucleus with the cortical electrodes. Unlike the focal epilepsy patients, where the SOZ only comprises a pair of cortical electrodes, here we averaged the h2 values of each thalamus-cortex pair.\u003c/p\u003e \u003cp\u003e \u003cb\u003eHigh-frequency acute thalamic stimulation.\u003c/b\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec16\" class=\"Section2\"\u003e \u003ch2\u003eStimulation testing and protocol\u003c/h2\u003e \u003cp\u003eDuring post-operative monitoring in the EMU, we had the unique opportunity to perform acute thalamic stimulation testing on 14 patients of our cohort. Importantly, we could stimulate in both matched and unmatched thalamic nuclei. We applied continuous electrical stimulation with bipolar configuration on the thalamic site of interest, at a frequency of 100 Hz. This value was chosen because it was the highest allowed by the stimulator used, and it is in line with commonly used parameters for thalamic electrical stimulation (\u0026ge;\u0026thinsp;100Hz). The amplitude and pulse width were adjusted individually for each patient according to their comfort levels, but never exceeded the ranges of 1\u0026ndash;3 mA and 60\u0026ndash;300 us respectively.\u003c/p\u003e \u003cp\u003eFor clinical safety reasons, the stimulation testing occurred only after the patients exhibited spontaneous seizures and in seizure-free days. The duration of the stimulation train varied from a minimum of 5 seconds to a maximum of 30 seconds, repeated over 1\u0026ndash;10 times depending on the patient comfort and clinical state.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec17\" class=\"Section2\"\u003e \u003ch2\u003ePreprocessing and IED detection\u003c/h2\u003e \u003cp\u003eTo analyze the immediate effect of thalamic stimulation we extracted the following data: i) for baseline (Stim OFF), we collected 2 minutes of interictal activity in the SOZ channels before any kind of stimulation was delivered to the patient; ii) for active stimulation period (Stim ON), we included the SOZ recording for the whole duration of stimulation up to 10 s after stimulation was stopped\u003csup\u003e\u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e\u003c/sup\u003e. We excluded the first two seconds after the stimulation was delivered to reduce false positive detection due to transient stimulation-induced burst\u003csup\u003e\u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e40\u003c/span\u003e\u003c/sup\u003e. Unfortunately, for technical reasons, 4 patients had a shorter baseline period (average of 45 s). Prior epoching, we preprocessed the data in Anywave\u003csup\u003e\u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e38\u003c/span\u003e\u003c/sup\u003e by applying a band-pass filter (1-1000 Hz) and a notch filter (60Hz) and manually removing channels affected by stimulation artifact or other noise sources. On the remaining channels, we applied independent component analysis (ICA) and extracted n components (where n\u0026thinsp;=\u0026thinsp;number of channels \u0026minus;\u0026thinsp;1). This method is a common technique used to separate independent sources linearly mixed in several sensors. We visually inspected all the components and rejected the ones associated with stimulation artifacts. In this way, we could reliably apply an automatic spikes detection algorithm on the stim ON and stim OFF recordings.\u003c/p\u003e \u003cp\u003eEach intracranial recording of interest was then divided in 5 s non-overlapping epochs, on which we computed the detection of epileptiform discharges. To avoid any confound due to the stimulation artifact, we considered as IED only interictal spikes, and not high-frequency oscillations. For this, we used Delphos software\u003csup\u003e\u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e41\u003c/span\u003e,\u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e42\u003c/span\u003e\u003c/sup\u003e, which detects interictal spikes from SEEG recording from their time-frequency representation. This detector was previously demonstrated to allow for almost 100% specificity and more than 80% sensitivity. To ensure correct spikes detection we 1) computed the analysis also on non-SOZ channels, to confirm that no spikes would be detected and 2) visually inspected the data. We saved the IED number for each Stim ON and Stim OFF epoch.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec18\" class=\"Section2\"\u003e \u003ch2\u003eIED analysis\u003c/h2\u003e \u003cp\u003eTo quantify the impact of acute thalamic stimulation on the IED rate, we computed two complementary analyses, based on the data quality: a) if the IED were clearly detectable and frequent (occurring in at least 50% of the baseline epochs), we calculated the IED rate as the number of IED per minute, in each Stim ON and Stim OFF epoch; b) if more than half of the baseline epochs presented no IED, we calculated the probability of finding no spikes across all the Stim ON and Stim OFF epochs (\u003cb\u003eSupplementary Table\u0026nbsp;2\u003c/b\u003e). Even if consistent spiking was observed in baseline epochs, we included in analysis b also a subject (S26) that presented complete suppression of IEDs in stimulation epochs, hence precluding statistical analysis. For both analyses, we computed the percentage difference of these quantities for Stim ON vs stim OFF to characterize the stimulation effect. A beneficial effect of stimulation on epileptiform discharges would result in a smaller IED rate with stimulation ON for analysis a, and a positive change in analysis b, prompting to show that stimulation epochs would more likely show complete suppression of IED.\u003c/p\u003e \u003cp\u003e \u003cb\u003eClinical outcome of chronic thalamic stimulation.\u003c/b\u003e \u003c/p\u003e \u003cp\u003eWe retrospectively evaluated clinical outcome of hodologically matched PUL and ANT stimulation in seven individuals (S16, S27-S32) who received a chronic neuromodulation system as part of their clinical care.\u003c/p\u003e \u003cp\u003eDemographics and epilepsy information for these patients is reported in \u003cb\u003eSupplementary Table\u0026nbsp;1\u003c/b\u003e. All patients received a bilateral implant of the PUL or ANT, according to their SOZ location and as indicated by their clinical team following SEEG monitoring. Six patients received a Responsive Neurostimulator (RNS), while one patient received a DBS system (Medtronic Percept). Stimulation parameters were individually tailored for each patient according to their clinical need (\u003cb\u003eSupplementary Table\u0026nbsp;4\u003c/b\u003e). No clinical decisions were based upon this research.\u003c/p\u003e \u003cp\u003eWe retrospectively evaluated the percentage reduction in disabling seizures\u0026rsquo; frequency and the change in antiseizure medication. The follow-up window ranged from 9 to 30 months.\u003c/p\u003e \u003cp\u003e \u003cb\u003eData analysis and statistical procedures.\u003c/b\u003e \u003c/p\u003e \u003cp\u003eAll figures and analysis were performed using Matlab 2023b (Mathworks, California, US). For all box plots reported in this manuscript, the whiskers extend to the maximum spread not considering outliers, central, top, and bottom lines represent median, 25th, and 75th percentile, respectively.\u003c/p\u003e \u003cp\u003eFor all the analysis, we used the bootstrap method, which does not rely on distributional assumptions of the data, but rather resamples the quantities of interest to achieve empirical confidence intervals. For each comparison, we created a bootstrap sample (for n\u0026thinsp;=\u0026thinsp;10000 repetitions) by drawing a sample with replacement from the actual data points (n\u0026thinsp;=\u0026thinsp;10000 repetitions), and calculated the difference in means of the resampled data. We then applied two-tailed bootstrapping with significance levels of 0.05 (95% confidence interval), 0.01 (99% confidence interval), or 0.001 (99.9% confidence interval). The null hypothesis of no difference in the mean was rejected if 0 was not included in the confidence interval of the corresponding alpha value. If multiple comparisons were performed at once, we used a Bonferroni correction by dividing the alpha value by the number of pairwise comparisons being performed.\u003c/p\u003e \u003c/div\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eFunding\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis research received financial support from NIH (1R01NS122927-01A1- Gonzalez-Martinez \u0026amp; Sridevi Sarma, PIs and 1R01NS131428-01A1 - Elvira Pirondini \u0026amp; Gonzalez-Martinez, PIs).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthor contributions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eEP and JGM conceived the study. EP and JGM secured funding. AD, JGM, and EP designed the experiments. AD analyzed the data, with help from SN for correlation analysis and JH for DTI analysis. TA and JGM implemented patient recruitment, eligibility and monitoring and coordinated study management. AD, TA and JGM collected the data. AG, SS and ER helped with data management. AD, EP, JGM and TA helped with data interpretation. AD, JH, EP and JGM, wrote the paper and all authors contributed to its editing. All authors interpreted the data, contributed to the writing of the submitted manuscript, had final responsibility for the decision to submit the manuscript for publication, and vouch for the accuracy and completeness of the data and for the fidelity of the trial to the protocol.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting interests\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors declare no conflicts of interests in relation to this work.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eData availability\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe main data supporting the results in this study are available within the paper and its supplementary information. All data generated in this study and software will be uploaded in a public repository upon acceptance of \u0026nbsp;the manuscript. Raw data will be available upon reasonable request to the corresponding author.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAcknowledgements\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eWe would like to thank all the Epilepsy Monitoring Unit staff and Neurosurgery Staff of the University of Pittsburgh Medical Center for their support. We wish to thank all the participants and their families, who agreed to participate in the study and for the great insights into the real needs of people living with epilepsy.\u0026nbsp;\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eKwan, P. \u0026amp; Brodie, M. J. Definition of refractory epilepsy: defining the indefinable? \u003cem\u003eLancet Neurol.\u003c/em\u003e \u003cstrong\u003e9\u003c/strong\u003e, 27\u0026ndash;29 (2010).\u003c/li\u003e\n\u003cli\u003eThijs, R. D., Surges, R., O\u0026rsquo;Brien, T. J. \u0026amp; Sander, J. W. Epilepsy in adults. \u003cem\u003eLancet Lond. Engl.\u003c/em\u003e \u003cstrong\u003e393\u003c/strong\u003e, 689\u0026ndash;701 (2019).\u003c/li\u003e\n\u003cli\u003eGupta, S., Ryvlin, P., Faught, E., Tsong, W. \u0026amp; Kwan, P. Understanding the burden of focal epilepsy as a function of seizure frequency in the United States, Europe, and Brazil. \u003cem\u003eEpilepsia Open\u003c/em\u003e \u003cstrong\u003e2\u003c/strong\u003e, 199\u0026ndash;213 (2017).\u003c/li\u003e\n\u003cli\u003eGloor, P. \u0026amp; Fariello, R. G. Generalized epilepsy: some of its cellular mechanisms differ from those of focal epilepsy. \u003cem\u003eTrends Neurosci.\u003c/em\u003e \u003cstrong\u003e11\u003c/strong\u003e, 63\u0026ndash;68 (1988).\u003c/li\u003e\n\u003cli\u003eBadawy, R. A. B., Curatolo, J. M., Newton, M., Berkovic, S. F. \u0026amp; Macdonell, R. A. L. Changes in cortical excitability differentiate generalized and focal epilepsy. \u003cem\u003eAnn. Neurol.\u003c/em\u003e \u003cstrong\u003e61\u003c/strong\u003e, 324\u0026ndash;331 (2007).\u003c/li\u003e\n\u003cli\u003eFrench, J. A. Refractory Epilepsy: Clinical Overview. \u003cem\u003eEpilepsia\u003c/em\u003e \u003cstrong\u003e48\u003c/strong\u003e, 3\u0026ndash;7 (2007).\u003c/li\u003e\n\u003cli\u003eLi, A. \u003cem\u003eet al.\u003c/em\u003e Neural fragility as an EEG marker of the seizure onset zone. \u003cem\u003eNat. Neurosci.\u003c/em\u003e \u003cstrong\u003e24\u003c/strong\u003e, 1465\u0026ndash;1474 (2021).\u003c/li\u003e\n\u003cli\u003eVelasco, F. \u003cem\u003eet al.\u003c/em\u003e Deep brain stimulation for treatment of the epilepsies: the centromedian thalamic target. in \u003cem\u003eOperative Neuromodulation: Volume 2: Neural Networks Surgery\u003c/em\u003e (eds. Sakas, D. E. \u0026amp; Simpson, B. A.) 337\u0026ndash;342 (Springer, Vienna, 2007). doi:10.1007/978-3-211-33081-4_38.\u003c/li\u003e\n\u003cli\u003eFisher, R. \u003cem\u003eet al.\u003c/em\u003e Electrical stimulation of the anterior nucleus of thalamus for treatment of refractory epilepsy. \u003cem\u003eEpilepsia\u003c/em\u003e \u003cstrong\u003e51\u003c/strong\u003e, 899\u0026ndash;908 (2010).\u003c/li\u003e\n\u003cli\u003eFisher, R. S. Deep brain stimulation of thalamus for epilepsy. \u003cem\u003eNeurobiol. Dis.\u003c/em\u003e \u003cstrong\u003e179\u003c/strong\u003e, 106045 (2023).\u003c/li\u003e\n\u003cli\u003eFisher, R. S. \u0026amp; Velasco, A. L. Electrical brain stimulation for epilepsy. \u003cem\u003eNat. Rev. Neurol.\u003c/em\u003e \u003cstrong\u003e10\u003c/strong\u003e, 261\u0026ndash;270 (2014).\u003c/li\u003e\n\u003cli\u003eRyvlin, P., Rheims, S., Hirsch, L. J., Sokolov, A. \u0026amp; Jehi, L. Neuromodulation in epilepsy: state-of-the-art approved therapies. \u003cem\u003eLancet Neurol.\u003c/em\u003e \u003cstrong\u003e20\u003c/strong\u003e, 1038\u0026ndash;1047 (2021).\u003c/li\u003e\n\u003cli\u003eShine, J. M., Lewis, L. D., Garrett, D. D. \u0026amp; Hwang, K. The impact of the human thalamus on brain-wide information processing. \u003cem\u003eNat. Rev. Neurosci.\u003c/em\u003e \u003cstrong\u003e24\u003c/strong\u003e, 416\u0026ndash;430 (2023).\u003c/li\u003e\n\u003cli\u003eSteriade, M. Synchronized activities of coupled oscillators in the cerebral cortex and thalamus at different levels of vigilance. \u003cem\u003eCereb. Cortex\u003c/em\u003e \u003cstrong\u003e7\u003c/strong\u003e, 583\u0026ndash;604 (1997).\u003c/li\u003e\n\u003cli\u003eSteriade, M. Sleep, epilepsy and thalamic reticular inhibitory neurons. \u003cem\u003eTrends Neurosci.\u003c/em\u003e \u003cstrong\u003e28\u003c/strong\u003e, 317\u0026ndash;324 (2005).\u003c/li\u003e\n\u003cli\u003eSteriade, M. \u0026amp; Contreras, D. Spike-Wave Complexes and Fast Components of Cortically Generated Seizures. I. Role of Neocortex and Thalamus. \u003cem\u003eJ. Neurophysiol.\u003c/em\u003e \u003cstrong\u003e80\u003c/strong\u003e, 1439\u0026ndash;1455 (1998).\u003c/li\u003e\n\u003cli\u003eBehrens, T. E. J. \u003cem\u003eet al.\u003c/em\u003e Non-invasive mapping of connections between human thalamus and cortex using diffusion imaging. \u003cem\u003eNat. Neurosci.\u003c/em\u003e \u003cstrong\u003e6\u003c/strong\u003e, 750\u0026ndash;757 (2003).\u003c/li\u003e\n\u003cli\u003ePenfield, W. Epileptic automatism and the centrencephalic integrating system. \u003cem\u003eRes. Publ. - Assoc. Res. Nerv. Ment. Dis.\u003c/em\u003e \u003cstrong\u003e30\u003c/strong\u003e, 513\u0026ndash;528 (1952).\u003c/li\u003e\n\u003cli\u003eArthuis, M. \u003cem\u003eet al.\u003c/em\u003e Impaired consciousness during temporal lobe seizures is related to increased long-distance cortical\u0026ndash;subcortical synchronization. \u003cem\u003eBrain\u003c/em\u003e \u003cstrong\u003e132\u003c/strong\u003e, 2091\u0026ndash;2101 (2009).\u003c/li\u003e\n\u003cli\u003eZhang, Y. \u003cem\u003eet al.\u003c/em\u003e Thalamocortical structural connectivity abnormalities in drug-resistant generalized epilepsy: A diffusion tensor imaging study. \u003cem\u003eBrain Res.\u003c/em\u003e \u003cstrong\u003e1727\u003c/strong\u003e, 146558 (2020).\u003c/li\u003e\n\u003cli\u003eSalanova, V. \u003cem\u003eet al.\u003c/em\u003e Long-term efficacy and safety of thalamic stimulation for drug-resistant partial epilepsy. \u003cem\u003eNeurology\u003c/em\u003e \u003cstrong\u003e84\u003c/strong\u003e, 1017\u0026ndash;1025 (2015).\u003c/li\u003e\n\u003cli\u003eCukiert, A., Cukiert, C. M., Burattini, J. A. \u0026amp; Mariani, P. P. Seizure outcome during bilateral, continuous, thalamic centromedian nuclei deep brain stimulation in patients with generalized epilepsy: a prospective, open-label study. \u003cem\u003eSeizure\u003c/em\u003e \u003cstrong\u003e81\u003c/strong\u003e, 304\u0026ndash;309 (2020).\u003c/li\u003e\n\u003cli\u003eValent\u0026iacute;n, A. \u003cem\u003eet al.\u003c/em\u003e Deep brain stimulation of the centromedian thalamic nucleus for the treatment of generalized and frontal epilepsies. \u003cem\u003eEpilepsia\u003c/em\u003e \u003cstrong\u003e54\u003c/strong\u003e, 1823\u0026ndash;1833 (2013).\u003c/li\u003e\n\u003cli\u003eCooper, I. S. \u003cem\u003eet al.\u003c/em\u003e Evoked metabolic responses in the limbic-striate system produced by stimulation of anterior thalamic nucleus in man. \u003cem\u003eInt. J. Neurol.\u003c/em\u003e \u003cstrong\u003e18\u003c/strong\u003e, 179\u0026ndash;187 (1984).\u003c/li\u003e\n\u003cli\u003eKaufmann, E. \u003cem\u003eet al.\u003c/em\u003e Long-term evaluation of anterior thalamic deep brain stimulation for epilepsy in the European MORE registry. \u003cem\u003eEpilepsia\u003c/em\u003e \u003cstrong\u003e65\u003c/strong\u003e, 2438\u0026ndash;2458 (2024).\u003c/li\u003e\n\u003cli\u003eSalanova, V. \u003cem\u003eet al.\u003c/em\u003e The SANT\u0026Eacute; study at 10 years of follow-up: Effectiveness, safety, and sudden unexpected death in epilepsy. \u003cem\u003eEpilepsia\u003c/em\u003e \u003cstrong\u003e62\u003c/strong\u003e, 1306\u0026ndash;1317 (2021).\u003c/li\u003e\n\u003cli\u003eWeitz, A. J. \u003cem\u003eet al.\u003c/em\u003e Optogenetic fMRI reveals distinct, frequency-dependent networks recruited by dorsal and intermediate hippocampus stimulations. \u003cem\u003eNeuroImage\u003c/em\u003e \u003cstrong\u003e107\u003c/strong\u003e, 229\u0026ndash;241 (2015).\u003c/li\u003e\n\u003cli\u003eBhattacharyya, K. B. James Wenceslaus Papez, His Circuit, and Emotion. \u003cem\u003eAnn. Indian Acad. Neurol.\u003c/em\u003e \u003cstrong\u003e20\u003c/strong\u003e, 207\u0026ndash;210 (2017).\u003c/li\u003e\n\u003cli\u003eGuye, M. \u003cem\u003eet al.\u003c/em\u003e The role of corticothalamic coupling in human temporal lobe epilepsy. \u003cem\u003eBrain J. Neurol.\u003c/em\u003e \u003cstrong\u003e129\u003c/strong\u003e, 1917\u0026ndash;1928 (2006).\u003c/li\u003e\n\u003cli\u003eCapecchi, F., Mothersill, I. \u0026amp; Imbach, L. L. The medial pulvinar as a subcortical relay in temporal lobe status epilepticus. \u003cem\u003eSeizure\u003c/em\u003e \u003cstrong\u003e81\u003c/strong\u003e, 276\u0026ndash;279 (2020).\u003c/li\u003e\n\u003cli\u003ePizzo, F. \u003cem\u003eet al.\u003c/em\u003e The Ictal Signature of Thalamus and Basal Ganglia in Focal Epilepsy. \u003cem\u003eNeurology\u003c/em\u003e \u003cstrong\u003e96\u003c/strong\u003e, e280\u0026ndash;e293 (2021).\u003c/li\u003e\n\u003cli\u003eRosenberg, D. S. \u003cem\u003eet al.\u003c/em\u003e Involvement of Medial Pulvinar Thalamic Nucleus in Human Temporal Lobe Seizures. \u003cem\u003eEpilepsia\u003c/em\u003e \u003cstrong\u003e47\u003c/strong\u003e, 98\u0026ndash;107 (2006).\u003c/li\u003e\n\u003cli\u003eIkegaya, N. \u003cem\u003eet al.\u003c/em\u003e Thalamic stereoelectroencephalography for neuromodulation target selection: Proof of concept and review of literature of pulvinar direct electrical stimulation. \u003cem\u003eEpilepsia\u003c/em\u003e \u003cstrong\u003e65\u003c/strong\u003e, e79\u0026ndash;e86 (2024).\u003c/li\u003e\n\u003cli\u003eBurdette, D., Mirro, E. A., Lawrence, M. \u0026amp; Patra, S. E. Brain-responsive corticothalamic stimulation in the pulvinar nucleus for the treatment of regional neocortical epilepsy: A case series. \u003cem\u003eEpilepsia Open\u003c/em\u003e \u003cstrong\u003e6\u003c/strong\u003e, 611\u0026ndash;617 (2021).\u003c/li\u003e\n\u003cli\u003eVilela-Filho, O. \u003cem\u003eet al.\u003c/em\u003e Ventral intermediate nucleus deep brain stimulation for treatment-resistant focal aware motor seizures: illustrative case. \u003cem\u003eJ. Neurosurg. Case Lessons\u003c/em\u003e \u003cstrong\u003e5\u003c/strong\u003e, CASE2320 (2023).\u003c/li\u003e\n\u003cli\u003eCooper, I. S., Amin, I. \u0026amp; Gilman, S. The effect of chronic cerebellar stimulation upon epilepsy in man. \u003cem\u003eTrans. Am. Neurol. Assoc.\u003c/em\u003e \u003cstrong\u003e98\u003c/strong\u003e, 192\u0026ndash;196 (1973).\u003c/li\u003e\n\u003cli\u003eIlinsky, I. \u003cem\u003eet al.\u003c/em\u003e Human Motor Thalamus Reconstructed in 3D from Continuous Sagittal Sections with Identified Subcortical Afferent Territories. \u003cem\u003eeNeuro\u003c/em\u003e \u003cstrong\u003e5\u003c/strong\u003e, ENEURO.0060-18.2018 (2018).\u003c/li\u003e\n\u003cli\u003eColombet, B., Woodman, M., Badier, J. M. \u0026amp; B\u0026eacute;nar, C. G. AnyWave: a cross-platform and modular software for visualizing and processing electrophysiological signals. \u003cem\u003eJ. Neurosci. Methods\u003c/em\u003e \u003cstrong\u003e242\u003c/strong\u003e, 118\u0026ndash;126 (2015).\u003c/li\u003e\n\u003cli\u003eGranger, C. W. J. Investigating Causal Relations by Econometric Models and Cross-spectral Methods. \u003cem\u003eEconometrica\u003c/em\u003e \u003cstrong\u003e37\u003c/strong\u003e, 424\u0026ndash;438 (1969).\u003c/li\u003e\n\u003cli\u003eMilosevic, L. \u003cem\u003eet al.\u003c/em\u003e Physiological mechanisms of thalamic ventral intermediate nucleus stimulation for tremor suppression. \u003cem\u003eBrain\u003c/em\u003e \u003cstrong\u003e141\u003c/strong\u003e, 2142\u0026ndash;2155 (2018).\u003c/li\u003e\n\u003cli\u003eRoehri, N., Lina, J.-M., Mosher, J. C., Bartolomei, F. \u0026amp; Benar, C.-G. Time-Frequency Strategies for Increasing High-Frequency Oscillation Detectability in Intracerebral EEG. \u003cem\u003eIEEE Trans. Biomed. Eng.\u003c/em\u003e \u003cstrong\u003e63\u003c/strong\u003e, 2595\u0026ndash;2606 (2016).\u003c/li\u003e\n\u003cli\u003eRoehri, N., Pizzo, F., Bartolomei, F., Wendling, F. \u0026amp; B\u0026eacute;nar, C.-G. What are the assets and weaknesses of HFO detectors? A benchmark framework based on realistic simulations. \u003cem\u003ePLOS ONE\u003c/em\u003e \u003cstrong\u003e12\u003c/strong\u003e, e0174702 (2017).\u003c/li\u003e\n\u003cli\u003eVelasco, F. \u003cem\u003eet al.\u003c/em\u003e Centromedian Nucleus and Epilepsy. \u003cem\u003eJ. Clin. Neurophysiol.\u003c/em\u003e \u003cstrong\u003e38\u003c/strong\u003e, 485 (2021).\u003c/li\u003e\n\u003cli\u003eYamada, K. \u003cem\u003eet al.\u003c/em\u003e Somatotopic Organization of Thalamocortical Projection Fibers as Assessed with MR Tractography. \u003cem\u003eRadiology\u003c/em\u003e \u003cstrong\u003e242\u003c/strong\u003e, 840\u0026ndash;845 (2007).\u003c/li\u003e\n\u003cli\u003eJbabdi, S., Woolrich, M. W. \u0026amp; Behrens, T. E. J. Multiple-subjects connectivity-based parcellation using hierarchical Dirichlet process mixture models. \u003cem\u003eNeuroImage\u003c/em\u003e \u003cstrong\u003e44\u003c/strong\u003e, 373\u0026ndash;384 (2009).\u003c/li\u003e\n\u003cli\u003eAupy, J. \u003cem\u003eet al.\u003c/em\u003e Cortico-striatal synchronization in human focal seizures. \u003cem\u003eBrain\u003c/em\u003e \u003cstrong\u003e142\u003c/strong\u003e, 1282\u0026ndash;1295 (2019).\u003c/li\u003e\n\u003cli\u003eLopes da Silva, F., Pijn, J. P. \u0026amp; Boeijinga, P. Interdependence of EEG signals: linear vs. nonlinear associations and the significance of time delays and phase shifts. \u003cem\u003eBrain Topogr.\u003c/em\u003e \u003cstrong\u003e2\u003c/strong\u003e, 9\u0026ndash;18 (1989).\u003c/li\u003e\n\u003cli\u003eBartolomei, F., Wendling, F., Bellanger, J.-J., R\u0026eacute;gis, J. \u0026amp; Chauvel, P. Neural networks involving the medial temporal structures in temporal lobe epilepsy. \u003cem\u003eClin. Neurophysiol.\u003c/em\u003e \u003cstrong\u003e112\u003c/strong\u003e, 1746\u0026ndash;1760 (2001).\u003c/li\u003e\n\u003cli\u003eBertram, E. H., Mangan, P. S., Zhang, D., Scott, C. A. \u0026amp; Williamson, J. M. The Midline Thalamus: Alterations and a Potential Role in Limbic Epilepsy. \u003cem\u003eEpilepsia\u003c/em\u003e \u003cstrong\u003e42\u003c/strong\u003e, 967\u0026ndash;978 (2001).\u003c/li\u003e\n\u003cli\u003eCassidy, R. M. \u0026amp; Gale, K. Mediodorsal thalamus plays a critical role in the development of limbic motor seizures. \u003cem\u003eJ. Neurosci. Off. J. Soc. Neurosci.\u003c/em\u003e \u003cstrong\u003e18\u003c/strong\u003e, 9002\u0026ndash;9009 (1998).\u003c/li\u003e\n\u003cli\u003eEvangelista, E. \u003cem\u003eet al.\u003c/em\u003e Does the Thalamo-Cortical Synchrony Play a Role in Seizure Termination? \u003cem\u003eFront. Neurol.\u003c/em\u003e \u003cstrong\u003e6\u003c/strong\u003e, (2015).\u003c/li\u003e\n\u003cli\u003eSoulier, H. \u003cem\u003eet al.\u003c/em\u003e The anterior and pulvinar thalamic nuclei interactions in mesial temporal lobe seizure networks. \u003cem\u003eClin. Neurophysiol.\u003c/em\u003e \u003cstrong\u003e150\u003c/strong\u003e, 176\u0026ndash;183 (2023).\u003c/li\u003e\n\u003cli\u003eVetkas, A. \u003cem\u003eet al.\u003c/em\u003e Deep brain stimulation targets in epilepsy: Systematic review and meta-analysis of anterior and centromedian thalamic nuclei and hippocampus. \u003cem\u003eEpilepsia\u003c/em\u003e \u003cstrong\u003e63\u003c/strong\u003e, 513\u0026ndash;524 (2022).\u003c/li\u003e\n\u003cli\u003eKrauss, J. K. \u003cem\u003eet al.\u003c/em\u003e Technology of deep brain stimulation: current status and future directions. \u003cem\u003eNat. Rev. Neurol.\u003c/em\u003e \u003cstrong\u003e17\u003c/strong\u003e, 75\u0026ndash;87 (2021).\u003c/li\u003e\n\u003cli\u003eAgashe, S. \u003cem\u003eet al.\u003c/em\u003e Centromedian Nucleus of the Thalamus Deep Brain Stimulation for Genetic Generalized Epilepsy: A Case Report and Review of Literature. \u003cem\u003eFront. Hum. Neurosci.\u003c/em\u003e \u003cstrong\u003e16\u003c/strong\u003e, (2022).\u003c/li\u003e\n\u003cli\u003eGraves, N. M. \u0026amp; Fisher, R. S. Chapter 16 Neurostimulation for Epilepsy, Including a Pilot Study of Anterior Nucleus Stimulation. \u003cem\u003eNeurosurgery\u003c/em\u003e \u003cstrong\u003e52\u003c/strong\u003e, 127 (2005).\u003c/li\u003e\n\u003cli\u003eKim, S. H. \u003cem\u003eet al.\u003c/em\u003e Long-term follow-up of anterior thalamic deep brain stimulation in epilepsy: A 11-year, single center experience. \u003cem\u003eSeizure\u003c/em\u003e \u003cstrong\u003e52\u003c/strong\u003e, 154\u0026ndash;161 (2017).\u003c/li\u003e\n\u003cli\u003eGuo, W. \u003cem\u003eet al.\u003c/em\u003e Defining the optimal target for anterior thalamic deep brain stimulation in patients with drug-refractory epilepsy. \u003cem\u003eJ. Neurosurg.\u003c/em\u003e \u003cstrong\u003e134\u003c/strong\u003e, 1054\u0026ndash;1063 (2021).\u003c/li\u003e\n\u003cli\u003eOh, Y.-S. \u003cem\u003eet al.\u003c/em\u003e Cognitive improvement after long-term electrical stimulation of bilateral anterior thalamic nucleus in refractory epilepsy patients. \u003cem\u003eSeizure\u003c/em\u003e \u003cstrong\u003e21\u003c/strong\u003e, 183\u0026ndash;187 (2012).\u003c/li\u003e\n\u003cli\u003ePark, H. R. \u003cem\u003eet al.\u003c/em\u003e The Role of Anterior Thalamic Deep Brain Stimulation as an Alternative Therapy in Patients with Previously Failed Vagus Nerve Stimulation for Refractory Epilepsy. \u003cem\u003eStereotact. Funct. Neurosurg.\u003c/em\u003e \u003cstrong\u003e97\u003c/strong\u003e, 176\u0026ndash;182 (2019).\u003c/li\u003e\n\u003cli\u003eNorden, A. D. \u0026amp; Blumenfeld, H. The role of subcortical structures in human epilepsy. \u003cem\u003eEpilepsy Behav.\u003c/em\u003e \u003cstrong\u003e3\u003c/strong\u003e, 219\u0026ndash;231 (2002).\u003c/li\u003e\n\u003cli\u003eAvoli, M. A brief history on the oscillating roles of thalamus and cortex in absence seizures. \u003cem\u003eEpilepsia\u003c/em\u003e \u003cstrong\u003e53\u003c/strong\u003e, 779\u0026ndash;789 (2012).\u003c/li\u003e\n\u003cli\u003eW, P. Highest level seizures. \u003cem\u003eRes Publ Nerv Ment Dis\u003c/em\u003e \u003cstrong\u003e26\u003c/strong\u003e, 252\u0026ndash;271 (1947).\u003c/li\u003e\n\u003cli\u003eHwang, K., Bertolero, M. A., Liu, W. B. \u0026amp; D\u0026rsquo;Esposito, M. The Human Thalamus Is an Integrative Hub for Functional Brain Networks. \u003cem\u003eJ. Neurosci.\u003c/em\u003e \u003cstrong\u003e37\u003c/strong\u003e, 5594\u0026ndash;5607 (2017).\u003c/li\u003e\n\u003cli\u003eWu, T. Q. \u003cem\u003eet al.\u003c/em\u003e Multisite thalamic recordings to characterize seizure propagation in the human brain. \u003cem\u003eBrain\u003c/em\u003e \u003cstrong\u003e146\u003c/strong\u003e, 2792\u0026ndash;2802 (2023).\u003c/li\u003e\n\u003cli\u003eAracri, P., de Curtis, M., Forcaia, G. \u0026amp; Uva, L. Enhanced thalamo-hippocampal synchronization during focal limbic seizures. \u003cem\u003eEpilepsia\u003c/em\u003e \u003cstrong\u003e59\u003c/strong\u003e, 1774\u0026ndash;1784 (2018).\u003c/li\u003e\n\u003cli\u003eDestexhe, A., Contreras, D. \u0026amp; Steriade, M. Mechanisms underlying the synchronizing action of corticothalamic feedback through inhibition of thalamic relay cells. \u003cem\u003eJ. Neurophysiol.\u003c/em\u003e \u003cstrong\u003e79\u003c/strong\u003e, 999\u0026ndash;1016 (1998).\u003c/li\u003e\n\u003cli\u003eSchindler, K., Leung, H., Elger, C. E. \u0026amp; Lehnertz, K. Assessing seizure dynamics by analysing the correlation structure of multichannel intracranial EEG. \u003cem\u003eBrain\u003c/em\u003e \u003cstrong\u003e130\u003c/strong\u003e, 65\u0026ndash;77 (2007).\u003c/li\u003e\n\u003cli\u003eAiello, G. \u003cem\u003eet al.\u003c/em\u003e Functional network dynamics between the anterior thalamus and the cortex in deep brain stimulation for epilepsy. \u003cem\u003eBrain\u003c/em\u003e \u003cstrong\u003e146\u003c/strong\u003e, 4717\u0026ndash;4735 (2023).\u003c/li\u003e\n\u003cli\u003eMontgomery, E. B. \u0026amp; Gale, J. T. Mechanisms of action of deep brain stimulation(DBS). \u003cem\u003eNeurosci. Biobehav. Rev.\u003c/em\u003e \u003cstrong\u003e32\u003c/strong\u003e, 388\u0026ndash;407 (2008).\u003c/li\u003e\n\u003cli\u003eBo\u0026euml;x, C., Vulli\u0026eacute;moz, S., Spinelli, L., Pollo, C. \u0026amp; Seeck, M. High and low frequency electrical stimulation in non-lesional temporal lobe epilepsy. \u003cem\u003eSeizure\u003c/em\u003e \u003cstrong\u003e16\u003c/strong\u003e, 664\u0026ndash;669 (2007).\u003c/li\u003e\n\u003cli\u003eWyckhuys, T., Raedt, R., Vonck, K., Wadman, W. \u0026amp; Boon, P. Comparison of hippocampal Deep Brain Stimulation with high (130 Hz) and low frequency (5 Hz) on afterdischarges in kindled rats. \u003cem\u003eEpilepsy Res.\u003c/em\u003e \u003cstrong\u003e88\u003c/strong\u003e, 239\u0026ndash;246 (2010).\u003cstrong\u003e\u003c/strong\u003e\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":true,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"nature-portfolio","isNatureJournal":true,"hasQc":false,"allowDirectSubmit":false,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"","title":"Nature Portfolio","twitterHandle":"","acdcEnabled":false,"dfaEnabled":false,"editorialSystem":"ejp","reportingPortfolio":"","inReviewEnabled":true,"inReviewRevisionsEnabled":false},"keywords":"","lastPublishedDoi":"10.21203/rs.3.rs-5507011/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-5507011/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eTargeted electrical stimulation to specific thalamic regions offers a therapeutic approach for patients with refractory focal and generalized epilepsy who are not candidates for resective surgery. However, clinical outcome varies significantly, in particular for focal epilepsy, influenced by several factors, notably the precise anatomical and functional alignment between cortical regions generating epileptic discharges and the targeted thalamic stimulation sites.\u003c/p\u003e \u003cp\u003eHere we hypothesized that targeting thalamic nuclei with precise anatomical and functional connections to epileptic cortical areas (an approach that we refer to as hodological matching) could enhance neuromodulatory effects on focal epileptic discharges. To investigate this, we examined three thalamic subnuclei (pulvinar nucleus, anterior nucleus, and ventral intermediate nucleus/ventral oral posterior nuclei) in a retrospective study involving 32 focal epilepsy patients.\u003c/p\u003e \u003cp\u003eSpecifically, we first identified hodologically organized thalamocortical fibers connecting these nuclei to individual seizure onset zones (SOZs), combining neuroimaging and electrophysiological techniques. Further, analysis of 216 spontaneous seizures revealed the critical role of matched thalamic nuclei in seizure development and termination. Importantly, electrical stimulation of hodologically-matched thalamic nuclei immediately suppressed intracortical interictal epileptiform discharges, contrasting with ineffective outcomes from stimulation of unmatched targets. Finally, we retrospectively evaluated 7 patients with a chronic hodologically-matched neurostimulation system, which led to a clinically relevant reduction in seizure frequency (median reduction 86.5%), that outstands the current clinical practice of unmatched targets (39%).\u003c/p\u003e \u003cp\u003eOur results underscore the potential of hodological thalamic targeting to modulate epileptiform activity in specific cortical regions, highlighting the promise of precision medicine in thalamic neuromodulation for focal refractory epilepsy.\u003c/p\u003e","manuscriptTitle":"Thalamocortical Hodology to Personalize Electrical Stimulation for Focal Epilepsy","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2024-11-28 15:50:22","doi":"10.21203/rs.3.rs-5507011/v1","editorialEvents":[],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"nature-communications","isNatureJournal":true,"hasQc":false,"allowDirectSubmit":false,"externalIdentity":"NCOMMS","sideBox":"Learn more about [Nature Communications](http://www.nature.com/ncomms/)","snPcode":"","submissionUrl":"https://mts-ncomms.nature.com/","title":"Nature Communications","twitterHandle":"","acdcEnabled":true,"dfaEnabled":true,"editorialSystem":"ejp","reportingPortfolio":"Nature Communications","inReviewEnabled":true,"inReviewRevisionsEnabled":false}}],"origin":"","ownerIdentity":"748b0546-cae0-4014-a1ef-c0aaeefb13ff","owner":[],"postedDate":"November 28th, 2024","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"published-in-journal","subjectAreas":[{"id":40769125,"name":"Health sciences/Neurology/Neurological disorders/Epilepsy"},{"id":40769126,"name":"Health sciences/Medical research/Translational research"}],"tags":[],"updatedAt":"2025-10-21T07:09:06+00:00","versionOfRecord":{"articleIdentity":"rs-5507011","link":"https://doi.org/10.1038/s41467-025-64922-w","journal":{"identity":"nature-communications","isVorOnly":false,"title":"Nature Communications"},"publishedOn":"2025-10-20 04:00:00","publishedOnDateReadable":"October 20th, 2025"},"versionCreatedAt":"2024-11-28 15:50:22","video":"","vorDoi":"10.1038/s41467-025-64922-w","vorDoiUrl":"https://doi.org/10.1038/s41467-025-64922-w","workflowStages":[]},"version":"v1","identity":"rs-5507011","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-5507011","identity":"rs-5507011","version":["v1"]},"buildId":"7rjqhiLT3MXkJMwkYKINL","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}