Quantification of enlarged deep medullary vein volumes in Sturge-Weber syndrome.

Background: Enlarged deep medullary veins (EDMVs) in patients with Sturge-Weber syndrome (SWS) may channel venous blood from the surface to the deep vein system in brain regions affected by the leptomeningeal venous malformation. Thus, the quantification of EDMV volume may provide an objective imaging marker for this vascular compensatory process. The present study proposes a novel analytical method to quantify enlarged EDMV volumes in the affected hemisphere of patients with unilateral SWS. Methods: Twenty young subjects, including 10 patients with unilateral SWS and 10 healthy siblings (age 14.5±6.7 and 16.0±7.0 years, respectively) underwent 3T brain MRI scanning using susceptibility-weighted imaging (SWI) and volumetric T1-weighted sequences. The proposed image analytic steps segmented EDMVs in white matter regions, defined on the volumetric T1-weighted images, by statistically associating the likelihood of intensity, location, and tubular shape on SWI. The volumes of the segmented EDMVs, calculated in each hemisphere, were compared between affected and unaffected hemispheres. EDMV volumes were also correlated with visually assessed EDMV scores, hemispheric white matter volumes, and cortical surface areas. Parametric tests including Pearson's correlation, unpaired and paired t-tests, were used. A P value <0.05 was considered statistically significant. Results: It was found that EDMVs were identified well in SWS-affected hemispheres while calcified regions were excluded. Mean EDMV volumes in the SWS-affected hemispheres were 10-12-fold greater than in the unaffected or healthy control hemispheres; while white matter volumes and cortical surface areas were lower. EDMV volumes in the SWS-affected hemispheres showed a strong positive correlation with the visual EDMV scores (r=0.88, P=0.001) and an inverse correlation with cortical surface area ratios (r=-0.65, P=0.04) but no correlation with white matter volume ratios. Conclusions: EDMVs were detected in the SWS-affected atrophic hemispheres reliably while avoiding calcified regions. The approach can be used to quantify enlarged deep cerebral veins in the human brain, which may provide a potential marker of cerebral venous remodeling.

Dynamic cortical and tractography atlases of proactive and reactive alpha and high-gamma activities.

Alpha waves-posterior dominant rhythms at 8-12 Hz reactive to eye opening and closure-are among the most fundamental EEG findings in clinical practice and research since Hans Berger first documented them in the early 20th century. Yet, the exact network dynamics of alpha waves in regard to eye movements remains unknown. High-gamma activity at 70-110 Hz is also reactive to eye movements and a summary measure of local cortical activation supporting sensorimotor or cognitive function. We aimed to build the first-ever brain atlases directly visualizing the network dynamics of eye movement-related alpha and high-gamma modulations, at cortical and white matter levels. We studied 28 patients (age: 5-20 years) who underwent intracranial EEG and electro-oculography recordings. We measured alpha and high-gamma modulations at 2167 electrode sites outside the seizure onset zone, interictal spike-generating areas and MRI-visible structural lesions. Dynamic tractography animated white matter streamlines modulated significantly and simultaneously beyond chance, on a millisecond scale. Before eye-closure onset, significant alpha augmentation occurred at the occipital and frontal cortices. After eye-closure onset, alpha-based functional connectivity was strengthened, while high gamma-based connectivity was weakened extensively in both intra-hemispheric and inter-hemispheric pathways involving the central visual areas. The inferior fronto-occipital fasciculus supported the strengthened alpha co-augmentation-based functional connectivity between occipital and frontal lobe regions, whereas the posterior corpus callosum supported the inter-hemispheric functional connectivity between the occipital lobes. After eye-opening offset, significant high-gamma augmentation and alpha attenuation occurred at occipital, fusiform and inferior parietal cortices. High gamma co-augmentation-based functional connectivity was strengthened, whereas alpha-based connectivity was weakened in the posterior inter-hemispheric and intra-hemispheric white matter pathways involving central and peripheral visual areas. Our results do not support the notion that eye closure-related alpha augmentation uniformly reflects feedforward or feedback rhythms propagating from lower to higher order visual cortex, or vice versa. Rather, proactive and reactive alpha waves involve extensive, distinct white matter networks that include the frontal lobe cortices, along with low- and high-order visual areas. High-gamma co-attenuation coupled to alpha co-augmentation in shared brain circuitry after eye closure supports the notion of an idling role for alpha waves during eye closure. These normative dynamic tractography atlases may improve understanding of the significance of EEG alpha waves in assessing the functional integrity of brain networks in clinical practice; they also may help elucidate the effects of eye movements on task-related brain network measures observed in cognitive neuroscience research.

Sevoflurane-induced high-frequency oscillations, effective connectivity and intraoperative classification of epileptic brain areas.

OBJECTIVE: To determine how sevoflurane anesthesia modulates intraoperative epilepsy biomarkers on electrocorticography, including high-frequency oscillation (HFO) effective connectivity (EC), and to investigate their relation to epileptogenicity and anatomical white matter. METHODS: We studied eight pediatric drug-resistant focal epilepsy patients who achieved seizure control after invasive monitoring and resective surgery. We visualized spatial distributions of the electrocorticography biomarkers at an oxygen baseline, three time-points while sevoflurane was increasing, and at a plateau of 2 minimum alveolar concentration (MAC) sevoflurane. HFO EC was combined with diffusion-weighted imaging, in dynamic tractography. RESULTS: Intraoperative HFO EC diffusely increased as a function of sevoflurane concentration, although most in epileptogenic sites (defined as those included in the resection); their ability to classify epileptogenicity was optimized at sevoflurane 2 MAC. HFO EC could be visualized on major white matter tracts, as a function of sevoflurane level. CONCLUSIONS: The results strengthened the hypothesis that sevoflurane-activated HFO biomarkers may help intraoperatively localize the epileptogenic zone. SIGNIFICANCE: Our results help characterize how HFOs at non-epileptogenic and epileptogenic networks respond to sevoflurane. It may be warranted to establish a normative HFO atlas incorporating the modifying effects of sevoflurane and major white matter pathways, as critical reference in epilepsy presurgical evaluation.

Deep Venous Remodeling in Unilateral Sturge-Weber Syndrome: Robust Hemispheric Differences and Clinical Correlates.

BACKGROUND: Enlarged deep medullary veins (EDMVs) in patients with Sturge-Weber syndrome (SWS) may provide compensatory venous drainage for brain regions affected by the leptomeningeal venous malformation (LVM). We evaluated the prevalence, extent, hemispheric differences, and clinical correlates of EDMVs in SWS. METHODS: Fifty children (median age: 4.5 years) with unilateral SWS underwent brain magnetic resonance imaging prospectively including susceptibility-weighted imaging (SWI); children aged 2.5 years or older also had a formal neurocognitive evaluation. The extent of EDMVs was assessed on SWI by using an EDMV hemispheric score, which was compared between patients with right and left SWS and correlated with clinical variables. RESULTS: EDMVs were present in 89% (24 of 27) of right and 78% (18 of 23) of left SWS brains. Extensive EDMVs (score >6) were more frequent in right (33%) than in left SWS (9%; P = 0.046) and commonly occurred in young children with right SWS. Patients with EDMV scores >4 had rare (less than monthly) seizures, whereas 35% (11 of 31) of patients with EDMV scores ≤4 had monthly or more frequent seizures (P = 0.003). In patients with right SWS and at least two LVM-affected lobes, higher EDMV scores were associated with higher intelligence quotient (P < 0.05). CONCLUSIONS: Enlarged deep medullary veins are common in unilateral SWS, but extensive EDMVs appear to develop more commonly and earlier in right hemispheric SWS. Deep venous remodeling may be a compensatory mechanism contributing to better clinical outcomes in some patients with SWS.

Multi-Scale Deep Learning of Clinically Acquired Multi-Modal MRI Improves the Localization of Seizure Onset Zone in Children With Drug-Resistant Epilepsy.

The present study investigates the effectiveness of a deep learning neural network for non-invasively localizing the seizure onset zone (SOZ) using multi-modal MRI data that are clinically acquired from children with drug-resistant epilepsy. A cortical parcellation was applied to localize the SOZ in cortical nodes of the epileptogenic hemisphere. At each node, the laminar surface analysis was followed to sample 1) the relative intensity of gray matter and white matter in multi-modal MRI and 2) the neighboring white matter connectivity using diffusion tractography edge strengths. A cross-validation was employed to train and test all layers of a multi-scale residual neural network (msResNet) that can classify SOZ node in an end-to-end fashion. A prediction probability of a given node belonging to the SOZ class was proposed as a non-invasive MRI marker of seizure onset likelihood. In an independent validation cohort, the proposed MRI marker provided a very large effect size of Cohen's d = 1.21 between SOZ and non-SOZ, and classified SOZ with a balanced accuracy of 0.75 in lesional and 0.67 in non-lesional MRI groups. The subsequent multi-variate logistic regression found the incorporation of the proposed MRI marker into interictal intracranial EEG (iEEG) markers further improves the differentiation between the epileptogenic focus (defined as SOZ resected during surgery) and non-epileptogenic sites (i.e., non-SOZ sites preserved during surgery) up to 15 % in non-lesional MRI group, suggesting that the proposed MRI marker could improve the localization of epileptogenic foci for successful pediatric epilepsy surgery.

Developmental organization of neural dynamics supporting auditory perception.

PURPOSE: A prominent view of language acquisition involves learning to ignore irrelevant auditory signals through functional reorganization, enabling more efficient processing of relevant information. Yet, few studies have characterized the neural spatiotemporal dynamics supporting rapid detection and subsequent disregard of irrelevant auditory information, in the developing brain. To address this unknown, the present study modeled the developmental acquisition of cost-efficient neural dynamics for auditory processing, using intracranial electrocorticographic responses measured in individuals receiving standard-of-care treatment for drug-resistant, focal epilepsy. We also provided evidence demonstrating the maturation of an anterior-to-posterior functional division within the superior-temporal gyrus (STG), which is known to exist in the adult STG. METHODS: We studied 32 patients undergoing extraoperative electrocorticography (age range: eight months to 28 years) and analyzed 2,039 intracranial electrode sites outside the seizure onset zone, interictal spike-generating areas, and MRI lesions. Patients were given forward (normal) speech sounds, backward-played speech sounds, and signal-correlated noises during a task-free condition. We then quantified sound processing-related neural costs at given time windows using high-gamma amplitude at 70-110 Hz and animated the group-level high-gamma dynamics on a spatially normalized three-dimensional brain surface. Finally, we determined if age independently contributed to high-gamma dynamics across brain regions and time windows. RESULTS: Group-level analysis of noise-related neural costs in the STG revealed developmental enhancement of early high-gamma augmentation and diminution of delayed augmentation. Analysis of speech-related high-gamma activity demonstrated an anterior-to-posterior functional parcellation in the STG. The left anterior STG showed sustained augmentation throughout stimulus presentation, whereas the left posterior STG showed transient augmentation after stimulus onset. We found a double dissociation between the locations and developmental changes in speech sound-related high-gamma dynamics. Early left anterior STG high-gamma augmentation (i.e., within 200 ms post-stimulus onset) showed developmental enhancement, whereas delayed left posterior STG high-gamma augmentation declined with development. CONCLUSIONS: Our observations support the model that, with age, the human STG refines neural dynamics to rapidly detect and subsequently disregard uninformative acoustic noises. Our study also supports the notion that the anterior-to-posterior functional division within the left STG is gradually strengthened for efficient speech-sound perception after birth.

Diffusion tractography predicts propagated high-frequency activity during epileptic spasms.

OBJECTIVE: To determine the structural networks that constrain propagation of ictal oscillations during epileptic spasm events, and compare the observed propagation patterns across patients with successful or unsuccessful surgical outcomes. METHODS: Subdural electrode recordings of 18 young patients (age 1-11 years) were analyzed during epileptic spasm events to determine ictal networks and quantify the amplitude and onset time of ictal oscillations across the cortical surface. Corresponding structural networks were generated with diffusion magnetic resonance imaging (MRI) tractography by seeding the cortical region associated with the earliest average oscillation onset time, and white matter pathways connecting active electrode regions within the ictal network were isolated. Properties of this structural network were used to predict oscillation onset times and amplitudes, and this relationship was compared across patients who did and did not achieve seizure freedom following resective surgery. RESULTS: Onset propagation patterns were relatively consistent across each patient's spasm events. An electrode's average ictal oscillation onset latency was most significantly associated with the length of direct corticocortical tracts connecting to the area with the earliest average oscillation onset (p < .001, model R2  = .54). Moreover, patients demonstrating a faster propagation of ictal oscillation signals within the corticocortical network were more likely to have seizure recurrence following resective surgery (p = .039). In addition, ictal oscillation amplitude was associated with connecting tractography length and weighted fractional anisotropy (FA) measures along these pathways (p = .002/.030, model R2  = .31/.25). Characteristics of analogous corticothalamic pathways did not show significant associations with ictal oscillation onset latency or amplitude. SIGNIFICANCE: Spatiotemporal propagation patterns of high-frequency activity in epileptic spasms align with length and FA measures from onset-originating corticocortical pathways. Considering the data in this individualized framework may help inform surgical decision-making and expectations of surgical outcomes.

Temporally and functionally distinct large-scale brain network dynamics supporting task switching.

OBJECTIVE: Our daily activities require frequent switches among competing responses at the millisecond time scale. We determined the spatiotemporal characteristics and functional significance of rapid, large-scale brain network dynamics during task switching. METHODS: This cross-sectional study investigated patients with drug-resistant focal epilepsy who played a Lumosity cognitive flexibility training game during intracranial electroencephalography (iEEG) recording. According to a given task rule, unpredictably switching across trials, participants had to swipe the screen in the direction the stimulus was pointing or moving. Using this data, we described the spatiotemporal characteristics of iEEG high-gamma augmentation occurring more intensely during switch than repeat trials, unattributable to the effect of task rule (pointing or moving), within-stimulus congruence (the direction of stimulus pointing and moving was same or different in a given trial), or accuracy of an immediately preceding response. Diffusion-weighted imaging (DWI) tractography determined whether distant cortical regions showing enhanced activation during task switch trials were directly connected by white matter tracts. Trial-by-trial iEEG analysis deduced whether the intensity of task switch-related high-gamma augmentation was altered through practice and whether high-gamma amplitude predicted the accuracy of an upcoming response among switch trials. RESULTS: The average number of completed trials during five-minute gameplay was 221.4 per patient (range: 171-285). Task switch trials increased the response times, whereas later trials reduced them. Analysis of iEEG signals sampled from 860 brain sites effectively elucidated the distinct spatiotemporal characteristics of task switch, task rule, and post-error-specific high-gamma modulations. Post-cue, task switch-related high-gamma augmentation was initiated in the right calcarine cortex after 260 ms, right precuneus after 330 ms, right entorhinal after 420 ms, and bilateral anterior middle-frontal gyri after 450 ms. DWI tractography successfully showed the presence of direct white matter tracts connecting the right visual areas to the precuneus and anterior middle-frontal regions but not between the right precuneus and anterior middle-frontal regions. Task-related high-gamma amplitudes in later trials were reduced in the calcarine, entorhinal and anterior middle-frontal regions, but increased in the precuneus. Functionally, enhanced post-cue precuneus high-gamma augmentation improved the accuracy of subsequent responses among switch trials. CONCLUSIONS: Our multimodal analysis uncovered two temporally and functionally distinct network dynamics supporting task switching. High-gamma augmentation in the visual-precuneus pathway may reflect the neural process facilitating an attentional shift to a given updated task rule. High-gamma activity in the visual-dorsolateral prefrontal pathway, rapidly reduced through practice, may reflect the cost of executing appropriate stimulus-response translation.

Naming-related spectral responses predict neuropsychological outcome after epilepsy surgery.

This prospective study determined the use of intracranially recorded spectral responses during naming tasks in predicting neuropsychological performance following epilepsy surgery. We recruited 65 patients with drug-resistant focal epilepsy who underwent preoperative neuropsychological assessment and intracranial EEG recording. The Clinical Evaluation of Language Fundamentals evaluated the baseline and postoperative language function. During extra-operative intracranial EEG recording, we assigned patients to undergo auditory and picture naming tasks. Time-frequency analysis determined the spatiotemporal characteristics of naming-related amplitude modulations, including high gamma augmentation at 70-110 Hz. We surgically removed the presumed epileptogenic zone based on the intracranial EEG and MRI abnormalities while maximally preserving the eloquent areas defined by electrical stimulation mapping. The multivariate regression model incorporating auditory naming-related high gamma augmentation predicted the postoperative changes in Core Language Score with r2 of 0.37 and in Expressive Language Index with r2 of 0.32. Independently of the effects of epilepsy and neuroimaging profiles, higher high gamma augmentation at the resected language-dominant hemispheric area predicted a more severe postoperative decline in Core Language Score and Expressive Language Index. Conversely, the model incorporating picture naming-related high gamma augmentation predicted the change in Receptive Language Index with an r2 of 0.50. Higher high gamma augmentation independently predicted a more severe postoperative decline in Receptive Language Index. Ancillary regression analysis indicated that naming-related low gamma augmentation and alpha/beta attenuation likewise independently predicted a more severe Core Language Score decline. The machine learning-based prediction model suggested that naming-related high gamma augmentation, among all spectral responses used as predictors, most strongly contributed to the improved prediction of patients showing a >5-point Core Language Score decline (reflecting the lower 25th percentile among patients). We generated the model-based atlas visualizing sites, which, if resected, would lead to such a language decline. With a 5-fold cross-validation procedure, the auditory naming-based model predicted patients who had such a postoperative language decline with an accuracy of 0.80. The model indicated that virtual resection of an electrical stimulation mapping-defined language site would have increased the relative risk of the Core Language Score decline by 5.28 (95% confidence interval: 3.47-8.02). Especially, that of an electrical stimulation mapping-defined receptive language site would have maximized it to 15.90 (95% confidence interval: 9.59-26.33). In summary, naming-related spectral responses predict neuropsychological outcomes after epilepsy surgery. We have provided our prediction model as an open-source material, which will indicate the postoperative language function of future patients and facilitate external validation at tertiary epilepsy centres.

Six-dimensional dynamic tractography atlas of language connectivity in the developing brain.

During a verbal conversation, our brain moves through a series of complex linguistic processing stages: sound decoding, semantic comprehension, retrieval of semantically coherent words, and overt production of speech outputs. Each process is thought to be supported by a network consisting of local and long-range connections bridging between major cortical areas. Both temporal and extratemporal lobe regions have functional compartments responsible for distinct language domains, including the perception and production of phonological and semantic components. This study provides quantitative evidence of how directly connected inter-lobar neocortical networks support distinct stages of linguistic processing across brain development. Novel six-dimensional tractography was used to intuitively visualize the strength and temporal dynamics of direct inter-lobar effective connectivity between cortical areas activated during each linguistic processing stage. We analysed 3401 non-epileptic intracranial electrode sites from 37 children with focal epilepsy (aged 5-20 years) who underwent extra-operative electrocorticography recording. Principal component analysis of auditory naming-related high-gamma modulations determined the relative involvement of each cortical area during each linguistic processing stage. To quantify direct effective connectivity, we delivered single-pulse electrical stimulation to 488 temporal and 1581 extratemporal lobe sites and measured the early cortico-cortical spectral responses at distant electrodes. Mixed model analyses determined the effects of naming-related high-gamma co-augmentation between connecting regions, age, and cerebral hemisphere on the strength of effective connectivity independent of epilepsy-related factors. Direct effective connectivity was strongest between extratemporal and temporal lobe site pairs, which were simultaneously activated between sentence offset and verbal response onset (i.e. response preparation period); this connectivity was approximately twice more robust than that with temporal lobe sites activated during stimulus listening or overt response. Conversely, extratemporal lobe sites activated during overt response were equally connected with temporal lobe language sites. Older age was associated with increased strength of inter-lobar effective connectivity especially between those activated during response preparation. The arcuate fasciculus supported approximately two-thirds of the direct effective connectivity pathways from temporal to extratemporal auditory language-related areas but only up to half of those in the opposite direction. The uncinate fasciculus consisted of <2% of those in the temporal-to-extratemporal direction and up to 6% of those in the opposite direction. We, for the first time, provided an atlas which quantifies and animates the strength, dynamics, and direction specificity of inter-lobar neural communications between language areas via the white matter pathways. Language-related effective connectivity may be strengthened in an age-dependent manner even after the age of 5.

Dynamic tractography-based localization of spike sources and animation of spike propagations.

OBJECTIVE: This study was undertaken to build and validate a novel dynamic tractography-based model for localizing interictal spike sources and visualizing monosynaptic spike propagations through the white matter. METHODS: This cross-sectional study investigated 1900 spike events recorded in 19 patients with drug-resistant temporal lobe epilepsy (TLE) who underwent extraoperative intracranial electroencephalography (iEEG) and resective surgery. Twelve patients had mesial TLE (mTLE) without a magnetic resonance imaging-visible mass lesion. The remaining seven had a mass lesion in the temporal lobe neocortex. We identified the leading and lagging sites, defined as those initially and subsequently (but within ≤50 ms) showing spike-related augmentation of broadband iEEG activity. In each patient, we estimated the sources of 100 spike discharges using the latencies at given electrode sites and diffusion-weighted imaging-based streamline length measures. We determined whether the spatial relationship between the estimated spike sources and resection was associated with postoperative seizure outcomes. We generated videos presenting the spatiotemporal change of spike-related fiber activation sites by estimating the propagation velocity using the streamline length and spike latency measures. RESULTS: The spike propagation velocity from the source was 1.03 mm/ms on average (95% confidence interval = .91-1.15) across 133 tracts noted in the 19 patients. The estimated spike sources in mTLE patients with International League Against Epilepsy Class 1 outcome were more likely to be in the resected area (83.9% vs. 72.3%, φ = .137, p < .001) and in the medial temporal lobe region (80.5% vs. 72.5%, φ = .090, p = .002) than those associated with the Class ≥2 outcomes. The resulting video successfully animated spike propagations, which were confined within the temporal lobe in mTLE but involved extratemporal lobe areas in lesional TLE. SIGNIFICANCE: We have, for the first time, provided dynamic tractography visualizing the spatiotemporal profiles of rapid propagations of interictal spikes through the white matter. Dynamic tractography has the potential to serve as a unique epilepsy biomarker.

Objective interictal electrophysiology biomarkers optimize prediction of epilepsy surgery outcome.

Researchers have looked for rapidly- and objectively-measurable electrophysiology biomarkers that accurately localize the epileptogenic zone. Promising candidates include interictal high-frequency oscillation and phase-amplitude coupling. Investigators have independently created the toolboxes that compute the high-frequency oscillation rate and the severity of phase-amplitude coupling. This study of 135 patients determined what toolboxes and analytic approaches would optimally classify patients achieving post-operative seizure control. Four different detector toolboxes computed the rate of high-frequency oscillation at ≥80 Hz at intracranial EEG channels. Another toolbox calculated the modulation index reflecting the strength of phase-amplitude coupling between high-frequency oscillation and slow-wave at 3-4 Hz. We defined the completeness of resection of interictally-abnormal regions as the subtraction of high-frequency oscillation rate (or modulation index) averaged across all preserved sites from that averaged across all resected sites. We computed the outcome classification accuracy of the logistic regression-based standard model considering clinical, ictal intracranial EEG and neuroimaging variables alone. We then determined how well the incorporation of high-frequency oscillation/modulation index would improve the standard model mentioned above. To assess the anatomical variability across non-epileptic sites, we generated the normative atlas of detector-specific high-frequency oscillation and modulation index. Each atlas allowed us to compute the statistical deviation of high-frequency oscillation/modulation index from the non-epileptic mean. We determined whether the model accuracy would be improved by incorporating absolute or normalized high-frequency oscillation/modulation index as a biomarker assessing interictally-abnormal regions. We finally determined whether the model accuracy would be improved by selectively incorporating high-frequency oscillation verified to have high-frequency oscillatory components unattributable to a high-pass filtering effect. Ninety-five patients achieved successful seizure control, defined as International League against Epilepsy class 1 outcome. Multivariate logistic regression analysis demonstrated that complete resection of interictally-abnormal regions additively increased the chance of success. The model accuracy was further improved by incorporating z-score normalized high-frequency oscillation/modulation index or selective incorporation of verified high-frequency oscillation. The standard model had a classification accuracy of 0.75. Incorporation of normalized high-frequency oscillation/modulation index or verified high-frequency oscillation improved the classification accuracy up to 0.82. These outcome prediction models survived the cross-validation process and demonstrated an agreement between the model-based likelihood of success and the observed success on an individual basis. Interictal high-frequency oscillation and modulation index had a comparably additive utility in epilepsy presurgical evaluation. Our empirical data support the theoretical notion that the prediction of post-operative seizure outcomes can be optimized with the consideration of both interictal and ictal abnormalities.

Four-dimensional tractography animates propagations of neural activation via distinct interhemispheric pathways.

OBJECTIVE: To visualize and validate the dynamics of interhemispheric neural propagations induced by single-pulse electrical stimulation (SPES). METHODS: This methodological study included three patients with drug-resistant focal epilepsy who underwent measurement of cortico-cortical spectral responses (CCSRs) during bilateral stereo-electroencephalography recording. We delivered SPES to 83 electrode pairs and analyzed CCSRs recorded at 268 nonepileptic electrode sites. Diffusion-weighted imaging (DWI) tractography localized the interhemispheric white matter pathways as streamlines directly connecting two electrode sites. We localized and visualized the putative SPES-related fiber activation, at each 1-ms time window, based on the propagation velocity defined as the DWI-based streamline length divided by the early CCSR peak latency. RESULTS: The resulting movie, herein referred to as four-dimensional tractography, delineated the spatiotemporal dynamics of fiber activation via the corpus callosum and anterior commissure. Longer streamline length was associated with delayed peak latency and smaller amplitude of CCSRs. The cortical regions adjacent to each fiber activation site indeed exhibited CCSRs at the same time window. CONCLUSIONS: Our four-dimensional tractography successfully animated neural propagations via distinct interhemispheric pathways. SIGNIFICANCE: Our novel animation method has the potential to help investigators in addressing the mechanistic significance of the interhemispheric network dynamics supporting physiological function.

Deep Relational Reasoning for the Prediction of Language Impairment and Postoperative Seizure Outcome Using Preoperative DWI Connectome Data of Children With Focal Epilepsy.

Prolonged seizures in children with focal epilepsy (FE) may impair language functions and often reoccur after surgical intervention. This study is aimed at developing a novel deep relational reasoning network to investigate whether conventional diffusion-weighted imaging connectome analysis can be improved when predicting expressive and receptive scores of preoperative language impairments and classifying postoperative seizure outcomes (seizure freedom or recurrence) in individual FE children. To deeply reason the dependencies of axonal connections that are sparsely distributed in the whole brain, this study proposes the "dilated CNN + RN", a dilated convolutional neural network (CNN) combined with a relation network (RN). The performance of the dilated CNN + RN was evaluated using whole brain connectome data from 51 FE children. It was found that when compared with other state-of-the-art algorithms, the dilated CNN + RN led to an average improvement of 90.2% and 97.3% in predicting expressive and receptive language scores, and 2.2% and 4% improvement in classifying seizure freedom and seizure recurrence, respectively. These improvements were independent of the prefixed connectome densities. Also, the dilated CNN + RN could provide an explainable artificial intelligence (AI) model by computing gradient-based regression/classification activation maps. This mapping analysis revealed left superior-medial frontal cortex, bilateral hippocampi, and cerebellum as crucial hubs, facilitating important connections that were most predictive of language function and seizure refractoriness after surgery.

Dynamic tractography: Integrating cortico-cortical evoked potentials and diffusion imaging.

INTRODUCTION: Cortico-cortical evoked potentials (CCEPs) are utilized to identify effective networks in the human brain. Following single-pulse electrical stimulation of cortical electrodes, evoked responses are recorded from distant cortical areas. A negative deflection (N1) which occurs 10-50 â€‹ms post-stimulus is considered to be a marker for direct cortico-cortical connectivity. However, with CCEPs alone it is not possible to observe the white matter pathways that conduct the signal or accurately predict N1 amplitude and latency at downstream recoding sites. Here, we develop a new approach, termed "dynamic tractography," which integrates CCEP data with diffusion-weighted imaging (DWI) data collected from the same patients. This innovative method allows greater insights into cortico-cortical networks than provided by each method alone and may improve the understanding of large-scale networks that support cognitive functions. The dynamic tractography model produces several fundamental hypotheses which we investigate: 1) DWI-based pathlength predicts N1 latency; 2) DWI-based pathlength negatively predicts N1 voltage; and 3) fractional anisotropy (FA) along the white matter path predicts N1 propagation velocity. METHODS: Twenty-three neurosurgical patients with drug-resistant epilepsy underwent both extraoperative CCEP recordings and preoperative DWI scans. Subdural grids of 3 â€‹mm diameter electrodes were used for stimulation and recording, with 98-128 eligible electrodes per patient. CCEPs were elicited by trains of 1 â€‹Hz stimuli with an intensity of 5 â€‹mA and recorded at a sample rate of 1 â€‹kHz. N1 peak and latency were defined as the maximum of a negative deflection within 10-50 â€‹ms post-stimulus with a z-score > 5 relative to baseline. Electrodes and DWI were coregistered to construct electrode connectomes for white matter quantification. RESULTS: Clinical variables (age, sex, number of anti-epileptic drugs, handedness, and stimulated hemisphere) did not correlate with the key outcome measures (N1 peak amplitude, latency, velocity, or DWI pathlength). All subjects and electrodes were therefore pooled into a group-level analysis to determine overall patterns. As hypothesized, DWI path length positively predicted N1 latency (R2 â€‹= â€‹0.81, ߠ​= â€‹1.51, p â€‹= â€‹4.76e-16) and negatively predicted N1 voltage (R2 â€‹= â€‹0.79, ߠ​= â€‹-0.094, p â€‹= â€‹9.30e-15), while FA predicted N1 propagation velocity (R2 â€‹= â€‹0.35, ߠ​= â€‹48.0, p â€‹= â€‹0.001). CONCLUSION: We have demonstrated that the strength and timing of the CCEP N1 is dependent on the properties of the underlying white matter network. Integrated CCEP and DWI visualization allows robust localization of intact axonal pathways which effectively interconnect eloquent cortex.

Novel Deep Learning Network Analysis of Electrical Stimulation Mapping-Driven Diffusion MRI Tractography to Improve Preoperative Evaluation of Pediatric Epilepsy.

OBJECTIVE: To investigate the clinical utility of deep convolutional neural network (DCNN) tract classification as a new imaging tool in the preoperative evaluation of children with focal epilepsy (FE). METHODS: A DCNN tract classification deeply learned spatial trajectories of DWI white matter pathways linking electrical stimulation mapping (ESM) findings from 89 children with FE, and then automatically identified white matter pathways associated with eloquent functions (i.e., primary motor, language, and vision). Clinical utility was examined by 1) measuring the nearest distance between DCNN-determined pathways and ESM, 2) evaluating the effectiveness of DCNN-determined pathways to optimize surgical margins via Kalman filter analysis, and 3) evaluating how accurately changes in DCNN-determined language pathway volume can predict changes in language ability via canonical correlation analysis. RESULTS: DCNN tract classification outperformed other existing methods, achieving an excellent accuracy of 98 % while non-invasively detecting eloquent areas within the spatial resolution of ESM (i.e., 1 cm). The Kalman filter analysis found that the preservation of brain areas within a surgical margin determined by DCNN tract classification predicted lack of postoperative deficit with a high accuracy of 92 %. Postoperative change of DCNN-determined language pathway volume showed a significant correlation with postoperative changes in language ability (R = 0.7, p   0.001). CONCLUSION: Our findings demonstrate that postoperative functional deficits substantially differ according to the extent of resected white matter, and that DCNN tract classification may offer key translational information by identifying these pathways in pediatric epilepsy surgery. SIGNIFICANCE: DCNN tract classification may be an effective tool to improve surgical outcome of children with FE.

Four-dimensional map of direct effective connectivity from posterior visual areas.

Lower- and higher-order visual cortices in the posterior brain, ranging from the medial- and lateral-occipital to fusiform regions, are suggested to support visual object recognition, whereas the frontal eye field (FEF) plays a role in saccadic eye movements which optimize visual processing. Previous studies using electrophysiology and functional MRI techniques have reported that tasks requiring visual object recognition elicited cortical activation sequentially in the aforementioned posterior visual regions and FEFs. The present study aims to provide unique evidence of direct effective connectivity outgoing from the posterior visual regions by measuring the early component (10-50 â€‹ms) of cortico-cortical spectral responses (CCSRs) elicited by weak single-pulse direct cortical electrical stimulation. We studied 22 patients who underwent extraoperative intracranial EEG recording for clinical localization of seizure foci and functionally-important brain regions. We used animations to visualize the spatiotemporal dynamics of gamma band CCSRs elicited by stimulation of three different posterior visual regions. We quantified the strength of CCSR-defined effective connectivity between the lower- and higher-order posterior visual regions as well as from the posterior visual regions to the FEFs. We found that effective connectivity within the posterior visual regions was larger in the feedforward (i.e., lower-to higher-order) direction compared to the opposite direction. Specifically, connectivity from the medial-occipital region was largest to the lateral-occipital region, whereas that from the lateral-occipital region was largest to the fusiform region. Among the posterior visual regions, connectivity to the FEF was largest from the lateral-occipital region and the mean peak latency of CCSR propagation from the lateral-occipital region to FEF was 26 â€‹ms. Our invasive study of the human brain using a stimulation-based intervention supports the model that the posterior visual regions have direct cortico-cortical connectivity pathways in which neural activity is transferred preferentially from the lower-to higher-order areas. The human brain has direct cortico-cortical connectivity allowing a rapid transfer of neural activity from the lateral-occipital region to the FEF.

Quantitative analysis of intracranial electrocorticography signals using the concept of statistical parametric mapping.

Statistical parametric mapping (SPM) is a technique with which one can delineate brain activity statistically deviated from the normative mean, and has been commonly employed in noninvasive neuroimaging and EEG studies. Using the concept of SPM, we developed a novel technique for quantification of the statistical deviation of an intracranial electrocorticography (ECoG) measure from the nonepileptic mean. We validated this technique using data previously collected from 123 patients with drug-resistant epilepsy who underwent resective epilepsy surgery. We determined how the measurement of statistical deviation of modulation index (MI) from the non-epileptic mean (rated by z-score) improved the performance of seizure outcome classification model solely based on conventional clinical, seizure onset zone (SOZ), and neuroimaging variables. Here, MI is a summary measure quantifying the strength of in-situ coupling between high-frequency activity at >150 Hz and slow wave at 3-4 Hz. We initially generated a normative MI atlas showing the mean and standard deviation of slow-wave sleep MI of neighboring non-epileptic channels of 47 patients, whose ECoG sampling involved all four lobes. We then calculated 'MI z-score' at each electrode site. SOZ had a greater 'MI z-score' compared to non-SOZ in the remaining 76 patients. Subsequent multivariate logistic regression analysis and receiver operating characteristic analysis to the combined data of all patients revealed that the full regression model incorporating all predictor variables, including SOZ and 'MI z-score', best classified the seizure outcome with sensitivity/specificity of 0.86/0.76. The model excluding 'MI z-score' worsened its sensitivity/specificity to 0.86/0.48. Furthermore, the leave-one-out analysis successfully cross-validated the full regression model. Measurement of statistical deviation of MI from the non-epileptic mean on invasive recording is technically feasible. Our analytical technique can be used to evaluate the utility of ECoG biomarkers in epilepsy presurgical evaluation.

Novel diffusion tractography methodology using Kalman filter prediction to improve preoperative benefit-risk analysis in pediatric epilepsy surgery.

OBJECTIVE: In this study the authors investigated the clinical reliability of diffusion weighted imaging maximum a posteriori probability (DWI-MAP) analysis with Kalman filter prediction in pediatric epilepsy surgery. This approach can yield a suggested resection margin as a dynamic variable based on preoperative DWI-MAP pathways. The authors sought to determine how well the suggested margin would have maximized occurrence of postoperative seizure freedom (benefit) and minimized occurrence of postoperative neurological deficits (risk). METHODS: The study included 77 pediatric patients with drug-resistant focal epilepsy (age 10.0 ± 4.9 years) who underwent resection of their presumed epileptogenic zone. In preoperative DWI tractography from the resected hemisphere, 9 axonal pathways, Ci=1-9, were identified using DWI-MAP as follows: C1-3 supporting face, hand, and leg motor areas; C4 connecting Broca's and Wernicke's areas; C5-8 connecting Broca's, Wernicke's, parietal, and premotor areas; and C9 connecting the occipital lobe and lateral geniculate nucleus. For each Ci, the resection margin, di, was measured by the minimal Euclidean distance between the voxels of Ci and the resection boundary determined by spatially coregistered postoperative MRI. If Ci was resected, di was assumed to be negative (calculated as -1 × average Euclidean distance between every voxel inside the resected Ci volume, ri). Kalman filter prediction was then used to estimate an optimal resection margin, d*i, to balance benefit and risk by approximating the relationship between di and ri. Finally, the authors defined the preservation zone of Ci that can balance the probability of benefit and risk by expanding the cortical area of Ci up to d*i on the 3D cortical surface. RESULTS: In the whole group (n = 77), nonresection of the preoperative preservation zone (i.e., actual resection margin d*i greater than the Kalman filter-defined d*i) accurately predicted the absence of postoperative motor (d*1-3: 0.93 at seizure-free probability of 0.80), language (d*4-8: 0.91 at seizure-free probability of 0.81), and visual deficits (d*9: 0.90 at seizure-free probability of 0.75), suggesting that the preservation of preoperative Ci within d*i supports a balance between postoperative functional deficit and seizure freedom. The subsequent subgroup analyses found that preservation of preoperative Ci =1-4,9 within d*i =1-4,9 may provide accurate deficit predictions independent of age and seizure frequency, suggesting that the DWI-based surgical margin can be effective for surgical planning even in young children and across a range of epilepsy severity. CONCLUSIONS: Integrating DWI-MAP analysis with Kalman filter prediction may help guide epilepsy surgery by visualizing the margins of the eloquent white matter pathways to be preserved.

Objective Detection of Eloquent Axonal Pathways to Minimize Postoperative Deficits in Pediatric Epilepsy Surgery using Diffusion Tractography and Convolutional Neural Networks.

Convolutional neural networks (CNNs) have recently been used in biomedical imaging applications with great success. In this paper, we investigated the classi?cation performance of CNN models on diffusion weighted imaging (DWI) streamlines de?ned by functional MRI (fMRI) and electrical stimulation mapping (ESM). To learn a set of discriminative and interpretable features from the extremely unbalanced dataset, we evaluated different CNN architectures with multiple loss functions (e.g., focal loss and center loss) and a soft attention mechanism, and compared our models with current state-ofthe-art methods. Through extensive experiments on streamlines collected from 70 healthy children and 70 children with focal epilepsy, we demonstrated that our deep CNN model with focal and central losses and soft attention outperforms all existing models in the literature and provides clinically acceptable accuracy (73 -100%) for the objective detection of functionally-important white matter pathways including ESM determined eloquent areas such as primary motor, aphasia, speech arrest, auditory, and visual functions. The ?ndings of this study encourage further investigations to determine if DWICNN analysis can serve as a noninvasive diagnostic tool during pediatric presurgical planning by estimating not only the location of essential cortices at the gyral level, but also the underlying ?bers connecting these cortical areas, to minimize or predict postsurgical functional de?cits. This study translates an advanced CNN model to clinical practice in the pediatric population where currently available approaches (e.g., ESM, fMRI) are suboptimal. The implementation will be released at https://github. com/HaotianMXu/Brain-?ber-classi?cation-using-CNNs.