Mapping hippocampal and thalamic atrophy in epilepsy: A 7-T magnetic resonance imaging study.

OBJECTIVE: Epilepsy patients are often grouped together by clinical variables. Quantitative neuroimaging metrics can provide a data-driven alternative for grouping of patients. In this work, we leverage ultra-high-field 7-T structural magnetic resonance imaging (MRI) to characterize volumetric atrophy patterns across hippocampal subfields and thalamic nuclei in drug-resistant focal epilepsy. METHODS: Forty-two drug-resistant epilepsy patients and 13 controls with 7-T structural neuroimaging were included in this study. We measured hippocampal subfield and thalamic nuclei volumetry, and applied an unsupervised machine learning algorithm, Latent Dirichlet Allocation (LDA), to estimate atrophy patterns across the hippocampal subfields and thalamic nuclei of patients. We studied the association between predefined clinical groups and the estimated atrophy patterns. Additionally, we used hierarchical clustering on the LDA factors to group patients in a data-driven approach. RESULTS: In patients with mesial temporal sclerosis (MTS), we found a significant decrease in volume across all ipsilateral hippocampal subfields (false discovery rate-corrected p [pFDR] < .01) as well as in some ipsilateral (pFDR < .05) and contralateral (pFDR < .01) thalamic nuclei. In left temporal lobe epilepsy (L-TLE) we saw ipsilateral hippocampal and some bilateral thalamic atrophy (pFDR < .05), whereas in right temporal lobe epilepsy (R-TLE) extensive bilateral hippocampal and thalamic atrophy was observed (pFDR < .05). Atrophy factors demonstrated that our MTS cohort had two atrophy phenotypes: one that affected the ipsilateral hippocampus and one that affected the ipsilateral hippocampus and bilateral anterior thalamus. Atrophy factors demonstrated posterior thalamic atrophy in R-TLE, whereas an anterior thalamic atrophy pattern was more common in L-TLE. Finally, hierarchical clustering of atrophy patterns recapitulated clusters with homogeneous clinical properties. SIGNIFICANCE: Leveraging 7-T MRI, we demonstrate widespread hippocampal and thalamic atrophy in epilepsy. Through unsupervised machine learning, we demonstrate patterns of volumetric atrophy that vary depending on disease subtype. Incorporating these atrophy patterns into clinical practice could help better stratify patients to surgical treatments and specific device implantation strategies.

iEEG-recon: A fast and scalable pipeline for accurate reconstruction of intracranial electrodes and implantable devices.

OBJECTIVE: Clinicians use intracranial electroencephalography (iEEG) in conjunction with noninvasive brain imaging to identify epileptic networks and target therapy for drug-resistant epilepsy cases. Our goal was to promote ongoing and future collaboration by automating the process of "electrode reconstruction," which involves the labeling, registration, and assignment of iEEG electrode coordinates on neuroimaging. We developed a standalone, modular pipeline that performs electrode reconstruction. We demonstrate our tool's compatibility with clinical and research workflows and its scalability on cloud platforms. METHODS: We created iEEG-recon, a scalable electrode reconstruction pipeline for semiautomatic iEEG annotation, rapid image registration, and electrode assignment on brain magnetic resonance imaging (MRI). Its modular architecture includes a clinical module for electrode labeling and localization, and a research module for automated data processing and electrode contact assignment. To ensure accessibility for users with limited programming and imaging expertise, we packaged iEEG-recon in a containerized format that allows integration into clinical workflows. We propose a cloud-based implementation of iEEG-recon and test our pipeline on data from 132 patients at two epilepsy centers using retrospective and prospective cohorts. RESULTS: We used iEEG-recon to accurately reconstruct electrodes in both electrocorticography and stereoelectroencephalography cases with a 30-min running time per case (including semiautomatic electrode labeling and reconstruction). iEEG-recon generates quality assurance reports and visualizations to support epilepsy surgery discussions. Reconstruction outputs from the clinical module were radiologically validated through pre- and postimplant T1-MRI visual inspections. We also found that our use of ANTsPyNet deep learning-based brain segmentation for electrode classification was consistent with the widely used FreeSurfer segmentations. SIGNIFICANCE: iEEG-recon is a robust pipeline for automating reconstruction of iEEG electrodes and implantable devices on brain MRI, promoting fast data analysis and integration into clinical workflows. iEEG-recon's accuracy, speed, and compatibility with cloud platforms make it a useful resource for epilepsy centers worldwide.

Mapping hippocampal glutamate in mesial temporal lobe epilepsy with glutamate weighted CEST (GluCEST) imaging.

Temporal lobe epilepsy (TLE) is one of the most common subtypes of focal epilepsy, with mesial temporal sclerosis (MTS) being a common radiological and histopathological finding. Accurate identification of MTS during presurgical evaluation confers an increased chance of good surgical outcome. Here we propose the use of glutamate-weighted chemical exchange saturation transfer (GluCEST) magnetic resonance imaging (MRI) at 7 Tesla for mapping hippocampal glutamate distribution in epilepsy, allowing to differentiate lesional from non-lesional mesial TLE. We demonstrate that a directional asymmetry index, which quantifies the relative difference between GluCEST contrast in hippocampi ipsilateral and contralateral to the seizure onset zone, can differentiate between sclerotic and non-sclerotic hippocampi, even in instances where traditional presurgical MRI assessments did not provide evidence of sclerosis. Overall, our results suggest that hippocampal glutamate mapping through GluCEST imaging is a valuable addition to the presurgical epilepsy evaluation toolbox.

Resting state functional connectivity demonstrates increased segregation in bilateral temporal lobe epilepsy.

OBJECTIVE: Temporal lobe epilepsy (TLE) is the most common type of focal epilepsy. An increasingly identified subset of patients with TLE consists of those who show bilaterally independent temporal lobe seizures. The purpose of this study was to leverage network neuroscience to better understand the interictal whole brain network of bilateral TLE (BiTLE). METHODS: In this study, using a multicenter resting state functional magnetic resonance imaging (rs-fMRI) data set, we constructed whole-brain functional networks of 19 patients with BiTLE, and compared them to those of 75 patients with unilateral TLE (UTLE). We quantified resting-state, whole-brain topological properties using metrics derived from network theory, including clustering coefficient, global efficiency, participation coefficient, and modularity. For each metric, we computed an average across all brain regions, and iterated this process across network densities. Curves of network density vs each network metric were compared between groups. Finally, we derived a combined metric, which we term the "integration-segregation axis," by combining whole-brain average clustering coefficient and global efficiency curves, and applying principal component analysis (PCA)-based dimensionality reduction. RESULTS: Compared to UTLE, BiTLE had decreased global efficiency (p = .031), and decreased whole brain average participation coefficient across a range of network densities (p = .019). Modularity maximization yielded a larger number of smaller communities in BiTLE than in UTLE (p = .020). Differences in network properties separate BiTLE and UTLE along the integration-segregation axis, with regions within the axis having a specificity of up to 0.87 for BiTLE. Along the integration-segregation axis, UTLE patients with poor surgical outcomes were distributed in the same regions as BiTLE, and network metrics confirmed similar patterns of increased segregation in both BiTLE and poor outcome UTLE. SIGNIFICANCE: Increased interictal whole-brain network segregation, as measured by rs-fMRI, is specific to BiTLE, as well as poor surgical outcome UTLE, and may assist in non-invasively identifying this patient population prior to intracranial electroencephalography or device implantation.

Quantifying trial-by-trial variability during cortico-cortical evoked potential mapping of epileptogenic tissue.

OBJECTIVE: Measuring cortico-cortical evoked potentials (CCEPs) is a promising tool for mapping epileptic networks, but it is not known how variability in brain state and stimulation technique might impact the use of CCEPs for epilepsy localization. We test the hypotheses that (1) CCEPs demonstrate systematic variability across trials and (2) CCEP amplitudes depend on the timing of stimulation with respect to endogenous, low-frequency oscillations. METHODS: We studied 11 patients who underwent CCEP mapping after stereo-electroencephalography electrode implantation for surgical evaluation of drug-resistant epilepsy. Evoked potentials were measured from all electrodes after each pulse of a 30 s, 1 Hz bipolar stimulation train. We quantified monotonic trends, phase dependence, and standard deviation (SD) of N1 (15-50 ms post-stimulation) and N2 (50-300 ms post-stimulation) amplitudes across the 30 stimulation trials for each patient. We used linear regression to quantify the relationship between measures of CCEP variability and the clinical seizure-onset zone (SOZ) or interictal spike rates. RESULTS: We found that N1 and N2 waveforms exhibited both positive and negative monotonic trends in amplitude across trials. SOZ electrodes and electrodes with high interictal spike rates had lower N1 and N2 amplitudes with higher SD across trials. Monotonic trends of N1 and N2 amplitude were more positive when stimulating from an area with higher interictal spike rate. We also found intermittent synchronization of trial-level N1 amplitude with low-frequency phase in the hippocampus, which did not localize the SOZ. SIGNIFICANCE: These findings suggest that standard approaches for CCEP mapping, which involve computing a trial-averaged response over a .2-1 Hz stimulation train, may be masking inter-trial variability that localizes to epileptogenic tissue. We also found that CCEP N1 amplitudes synchronize with ongoing low-frequency oscillations in the hippocampus. Further targeted experiments are needed to determine whether phase-locked stimulation could have a role in localizing epileptogenic tissue.

Hyperspectral wide-field-of-view imaging to study dynamic microcirculatory changes during hypoxia.

Hyperspectral imaging (HSI) provides a fast, reliable, and noninvasive way for the study of vascular microcirculation in animal models. Rapid hyperspectral imaging of large portions of the microcirculatory preparation is critical for understanding the function and regulation of vascular microcirculatory networks. This report presents the application of an off-the-shelf benchtop HSI linear scanning system to acquire larger field-of-view images of microcirculatory preparations. The HSI line detector was displaced perpendicular to the scanning direction to map larger areas, with a rate of displacement determined by the scanning rate and the exposure time. Analysis of the collected image was used to assess dynamic changes in the microcirculation. The system records dynamic changes in microvascular hemoglobin oxygen (HbO2) saturation and vascular morphology during hypoxia and reoxygenation and has similar acquisition speeds to commonly referenced spectral-scanning HSI systems. In addition, the HbO2 saturations collected via HSI closely correlate with those collected by phosphorescence quenching microscopy. The reported system enables dynamic functional imaging of the microcirculation for broad experimental and clinical applications.NEW & NOTEWORTHY This study presents a novel bench setup and algorithm to measure intravascular hemoglobin oxygen saturation in microcirculation. Wide-field hyperspectral imaging allows for rapid quantification of intravascular changes in hemoglobin saturation. The method described in this manuscript can expand the understanding of oxygen delivery to tissues in vivo.

Increased hemoglobin affinity for oxygen with GBT1118 improves hypoxia tolerance in sickle cell mice.

Therapeutic agents that increase the Hb affinity for oxygen (O2) could, in theory, lead to decreased O2 release from Hb and impose a hypoxic risk to tissues. In this study, GBT1118, an allosteric modifier of Hb affinity for O2, was used to assess the impact of increasing Hb affinity for O2 on brain tissue oxygenation, blood pressure, heart rate, O2 delivery, and tolerance to hypoxia in Townes transgenic sickle cell disease (SCD) mice. Brain oxygenation and O2 delivery were studied during normoxia and severe hypoxic challenges. Chronic treatment with GBT1118 increased Hb affinity for O2, reducing the Po2 for 50% HbO2 saturation (P50) in SCD mice from 31 mmHg to 18 mmHg. This treatment significantly reduced anemia, increasing hematocrit by 33%, improved cardiac output (CO), and O2 delivery and extraction. Chronically increasing Hb affinity for O2 with GBT1118 preserved cortical O2 tension during normoxia, improved cortical O2 tension during hypoxia, and increased tolerance to severe hypoxia in SCD mice. Independent of hematological changes induced by chronic treatment, a single dose of GBT1118 significantly improved tolerance to hypoxia, highlighting the benefits of increasing Hb affinity for O2 and consequently reducing sickling of RBCs in blood during hypoxia in SCD.NEW & NOTEWORTHY Chronic pharmacologically increased hemoglobin affinity for oxygen in sickle cell disease mice alleviated hematological consequences of sickle cell disease, increasing RBC half-life, hematocrit, and hemoglobin concentration, while also decreasing reticulocyte count. Additionally, chronically increased hemoglobin affinity for oxygen significantly improved survival as well as cortical tissue oxygenation in sickle cell disease mice during hypoxia, suggesting that oxygen delivery and utilization is improved by increased hemoglobin affinity for oxygen.

Tumor vascular status controls oxygen delivery facilitated by infused polymerized hemoglobins with varying oxygen affinity.

Oxygen (O2) delivery facilitated by hemoglobin (Hb)-based O2 carriers (HBOCs) is a promising strategy to increase the effectiveness of chemotherapeutics for treatment of solid tumors. However, the heterogeneous vascular structures present within tumors complicates evaluating the oxygenation potential of HBOCs within the tumor microenvironment. To account for spatial variations in the vasculature and tumor tissue that occur during tumor growth, we used a computational model to develop artificial tumor constructs. With these simulated tumors, we performed a polymerized human hemoglobin (hHb) (PolyhHb) enhanced oxygenation simulation accounting for differences in the physiologic characteristics of human and mouse blood. The results from this model were used to determine the potential effectiveness of different treatment options including a top load (low volume) and exchange (large volume) infusion of a tense quaternary state (T-State) PolyhHb, relaxed quaternary state (R-State) PolyhHb, and a non O2 carrying control. Principal component analysis (PCA) revealed correlations between the different regimes of effectiveness within the different simulated dosage options. In general, we found that infusion of T-State PolyhHb is more likely to decrease tissue hypoxia and modulate the metabolic rate of O2 consumption. Though the developed models are not a definitive descriptor of O2 carrier interaction in tumor capillary networks, we accounted for factors such as non-uniform vascular density and permeability that limit the applicability of O2 carriers during infusion. Finally, we have used these validated computational models to establish potential benchmarks to guide tumor treatment during translation of PolyhHb mediated therapies into clinical applications.

Red Blood Cell Metabolic Responses to Torpor and Arousal in the Hibernator Arctic Ground Squirrel.

Arctic ground squirrels provide a unique model to investigate metabolic responses to hibernation in mammals. During winter months these rodents are exposed to severe hypothermia, prolonged fasting, and hypoxemia. In the light of their role in oxygen transport/off-loading and owing to the absence of nuclei and organelles (and thus de novo protein synthesis capacity), mature red blood cells have evolved metabolic programs to counteract physiological or pathological hypoxemia. However, red blood cell metabolism in hibernation has not yet been investigated. Here we employed targeted and untargeted metabolomics approaches to investigate erythrocyte metabolism during entrance to torpor to arousal, with a high resolution of the intermediate time points. We report that torpor and arousal promote metabolism through glycolysis and pentose phosphate pathway, respectively, consistent with previous models of oxygen-dependent metabolic modulation in mature erythrocytes. Erythrocytes from hibernating squirrels showed up to 100-fold lower levels of biomarkers of reperfusion injury, such as the pro-inflammatory dicarboxylate succinate. Altered tryptophan metabolism during torpor was here correlated to the accumulation of potentially neurotoxic catabolites kynurenine, quinolinate, and picolinate. Arousal was accompanied by alterations of sulfur metabolism, including sudden spikes in a metabolite putatively identified as thiorphan (level 1 confidence)-a potent inhibitor of several metalloproteases that play a crucial role in nociception and inflammatory complication to reperfusion secondary to ischemia or hemorrhage. Preliminary studies in rats showed that intravenous injection of thiorphan prior to resuscitation mitigates metabolic and cytokine markers of reperfusion injury, etiological contributors to inflammatory complications after shock.

Image Annotation by Eye Tracking: Accuracy and Precision of Centerlines of Obstructed Small-Bowel Segments Placed Using Eye Trackers.

Small-bowel obstruction (SBO) is a common and important disease, for which machine learning tools have yet to be developed. Image annotation is a critical first step for development of such tools. This study assesses whether image annotation by eye tracking is sufficiently accurate and precise to serve as a first step in the development of machine learning tools for detection of SBO on CT. Seven subjects diagnosed with SBO by CT were included in the study. For each subject, an obstructed segment of bowel was chosen. Three observers annotated the centerline of the segment by manual fiducial placement and by visual fiducial placement using a Tobii 4c eye tracker. Each annotation was repeated three times. The distance between centerlines was calculated after alignment using dynamic time warping (DTW) and statistically compared to clinical thresholds for diagnosis of SBO. Intra-observer DTW distance between manual and visual centerlines was calculated as a measure of accuracy. These distances were 1.1 ± 0.2, 1.3 ± 0.4, and 1.8 ± 0.2 cm for the three observers and were less than 1.5 cm for two of three observers (P < 0.01). Intra- and inter-observer DTW distances between centerlines placed with each method were calculated as measures of precision. These distances were 0.6 ± 0.1 and 0.8 ± 0.2 cm for manual centerlines, 1.1 ± 0.4 and 1.9 ± 0.6 cm for visual centerlines, and were less than 3.0 cm in all cases (P < 0.01). Results suggest that eye tracking-based annotation is sufficiently accurate and precise for small-bowel centerline annotation for use in machine learning-based applications.

Artificial intelligence in epilepsy - applications and pathways to the clinic.

Artificial intelligence (AI) is rapidly transforming health care, and its applications in epilepsy have increased exponentially over the past decade. Integration of AI into epilepsy management promises to revolutionize the diagnosis and treatment of this complex disorder. However, translation of AI into neurology clinical practice has not yet been successful, emphasizing the need to consider progress to date and assess challenges and limitations of AI. In this Review, we provide an overview of AI applications that have been developed in epilepsy using a variety of data modalities: neuroimaging, electroencephalography, electronic health records, medical devices and multimodal data integration. For each, we consider potential applications, including seizure detection and prediction, seizure lateralization, localization of the seizure-onset zone and assessment for surgical or neurostimulation interventions, and review the performance of AI tools developed to date. We also discuss methodological considerations and challenges that must be addressed to successfully integrate AI into clinical practice. Our goal is to provide an overview of the current state of the field and provide guidance for leveraging AI in future to improve management of epilepsy.

Utility of intracranial EEG networks depends on re-referencing and connectivity choice.

Studies of intracranial EEG networks have been used to reveal seizure generators in patients with drug-resistant epilepsy. Intracranial EEG is implanted to capture the epileptic network, the collection of brain tissue that forms a substrate for seizures to start and spread. Interictal intracranial EEG measures brain activity at baseline, and networks computed during this state can reveal aberrant brain tissue without requiring seizure recordings. Intracranial EEG network analyses require choosing a reference and applying statistical measures of functional connectivity. Approaches to these technical choices vary widely across studies, and the impact of these technical choices on downstream analyses is poorly understood. Our objective was to examine the effects of different re-referencing and connectivity approaches on connectivity results and on the ability to lateralize the seizure onset zone in patients with drug-resistant epilepsy. We applied 48 pre-processing pipelines to a cohort of 125 patients with drug-resistant epilepsy recorded with interictal intracranial EEG across two epilepsy centres to generate intracranial EEG functional connectivity networks. Twenty-four functional connectivity measures across time and frequency domains were applied in combination with common average re-referencing or bipolar re-referencing. We applied an unsupervised clustering algorithm to identify groups of pre-processing pipelines. We subjected each pre-processing approach to three quality tests: (i) the introduction of spurious correlations; (ii) robustness to incomplete spatial sampling; and (iii) the ability to lateralize the clinician-defined seizure onset zone. Three groups of similar pre-processing pipelines emerged: common average re-referencing pipelines, bipolar re-referencing pipelines and relative entropy-based connectivity pipelines. Relative entropy and common average re-referencing networks were more robust to incomplete electrode sampling than bipolar re-referencing and other connectivity methods (Friedman test, Dunn-Sidák test P < 0.0001). Bipolar re-referencing reduced spurious correlations at non-adjacent channels better than common average re-referencing (Δ mean from machine ref = -0.36 versus -0.22) and worse in adjacent channels (Δ mean from machine ref = -0.14 versus -0.40). Relative entropy-based network measures lateralized the seizure onset hemisphere better than other measures in patients with temporal lobe epilepsy (Benjamini-Hochberg-corrected P < 0.05, Cohen's d: 0.60-0.76). Finally, we present an interface where users can rapidly evaluate intracranial EEG pre-processing choices to select the optimal pre-processing methods tailored to specific research questions. The choice of pre-processing methods affects downstream network analyses. Choosing a single method among highly correlated approaches can reduce redundancy in processing. Relative entropy outperforms other connectivity methods in multiple quality tests. We present a method and interface for researchers to optimize their pre-processing methods for deriving intracranial EEG brain networks.

Enhancing the Diagnostic Utility of ASL Imaging in Temporal Lobe Epilepsy through FlowGAN: An ASL to PET Image Translation Framework.

Background and Significance: Positron Emission Tomography (PET) using fluorodeoxyglucose (FDG-PET) is a standard imaging modality for detecting areas of hypometabolism associated with the seizure onset zone (SOZ) in temporal lobe epilepsy (TLE). However, FDG-PET is costly and involves the use of a radioactive tracer. Arterial Spin Labeling (ASL) offers an MRI-based quantification of cerebral blood flow (CBF) that could also help localize the SOZ, but its performance in doing so, relative to FDG-PET, is limited. In this study, we seek to improve ASL's diagnostic performance by developing a deep learning framework for synthesizing FDG-PET-like images from ASL and structural MRI inputs. Methods: We included 68 epilepsy patients, out of which 36 had well lateralized TLE. We compared the coupling between FDG-PET and ASL CBF values in different brain regions, as well as the asymmetry of these values across the brain. We additionally assessed each modality's ability to lateralize the SOZ across brain regions. Using our paired PET-ASL data, we developed FlowGAN, a generative adversarial neural network (GAN) that synthesizes PET-like images from ASL and T1-weighted MRI inputs. We tested our synthetic PET images against the actual PET images of subjects to assess their ability to reproduce clinically meaningful hypometabolism and asymmetries in TLE. Results: We found variable coupling between PET and ASL CBF values across brain regions. PET and ASL had high coupling in neocortical temporal and frontal brain regions (Spearman's r > 0.30, p < 0.05) but low coupling in mesial temporal structures (Spearman's r < 0.30, p > 0.05). Both whole brain PET and ASL CBF asymmetry values provided good separability between left and right TLE subjects, but PET (AUC = 0.96, 95% CI: [0.88, 1.00]) outperformed ASL (AUC = 0.81; 95% CI: [0.65, 0.96]). FlowGAN-generated images demonstrated high structural similarity to actual PET images (SSIM = 0.85). Globally, asymmetry values were better correlated between synthetic PET and original PET than between ASL CBF and original PET, with a mean correlation increase of 0.15 (95% CI: [0.07, 0.24], p<0.001, Cohen's d = 0.91). Furthermore, regions that had poor ASL-PET correlation (e.g. mesial temporal structures) showed the greatest improvement with synthetic PET images. Conclusions: FlowGAN improves ASL's diagnostic performance, generating synthetic PET images that closely mimic actual FDG-PET in depicting hypometabolism associated with TLE. This approach could improve non-invasive SOZ localization, offering a promising tool for epilepsy presurgical assessment. It potentially broadens the applicability of ASL in clinical practice and could reduce reliance on FDG-PET for epilepsy and other neurological disorders.

Disparities in seizure outcomes revealed by large language models.

OBJECTIVE: Large-language models (LLMs) can potentially revolutionize health care delivery and research, but risk propagating existing biases or introducing new ones. In epilepsy, social determinants of health are associated with disparities in care access, but their impact on seizure outcomes among those with access remains unclear. Here we (1) evaluated our validated, epilepsy-specific LLM for intrinsic bias, and (2) used LLM-extracted seizure outcomes to determine if different demographic groups have different seizure outcomes. MATERIALS AND METHODS: We tested our LLM for differences and equivalences in prediction accuracy and confidence across demographic groups defined by race, ethnicity, sex, income, and health insurance, using manually annotated notes. Next, we used LLM-classified seizure freedom at each office visit to test for demographic outcome disparities, using univariable and multivariable analyses. RESULTS: We analyzed 84 675 clinic visits from 25 612 unique patients seen at our epilepsy center. We found little evidence of bias in the prediction accuracy or confidence of outcome classifications across demographic groups. Multivariable analysis indicated worse seizure outcomes for female patients (OR 1.33, P ≤ .001), those with public insurance (OR 1.53, P ≤ .001), and those from lower-income zip codes (OR ≥1.22, P ≤ .007). Black patients had worse outcomes than White patients in univariable but not multivariable analysis (OR 1.03, P = .66). CONCLUSION: We found little evidence that our LLM was intrinsically biased against any demographic group. Seizure freedom extracted by LLM revealed disparities in seizure outcomes across several demographic groups. These findings quantify the critical need to reduce disparities in the care of people with epilepsy.

Subcortical Functional Connectivity Gradients in Temporal Lobe Epilepsy.

Background and Motivation: Functional gradients have been used to study differences in connectivity between healthy and diseased brain states, however this work has largely focused on the cortex. Because the subcortex plays a key role in seizure initiation in temporal lobe epilepsy (TLE), subcortical functional-connectivity gradients may help further elucidate differences between healthy brains and TLE, as well as differences between left (L)-TLE and right (R)-TLE. Methods: In this work, we calculated subcortical functional-connectivity gradients (SFGs) from resting-state functional MRI (rs-fMRI) by measuring the similarity in connectivity profiles of subcortical voxels to cortical gray matter voxels. We performed this analysis in 23 R-TLE patients and 32 L-TLE patients (who were otherwise matched for age, gender, disease specific characteristics, and other clinical variables), and 16 controls. To measure differences in SFGs between L-TLE and R-TLE, we quantified deviations in the average functional gradient distributions, as well as their variance, across subcortical structures. Results: We found an expansion, measured by increased variance, in the principal SFG of TLE relative to controls. When comparing the gradient across subcortical structures between L-TLE and R-TLE, we found that abnormalities in the ipsilateral hippocampal gradient distributions were significantly different between L-TLE and R-TLE. Conclusion: Our results suggest that expansion of the SFG is characteristic of TLE. Subcortical functional gradient differences exist between left and right TLE and are driven by connectivity changes in the hippocampus ipsilateral to the seizure onset zone.

Improved Seizure Onset-Zone Lateralization in Temporal Lobe Epilepsy using 7T Resting-State fMRI: A Direct Comparison with 3T.

Objective: Resting-state functional magnetic resonance imaging (rs-fMRI) at ultra high-field strengths (≥7T) is known to provide superior signal-to-noise and statistical power than comparable acquisitions at lower field strengths. In this study, we aim to provide a direct comparison of the seizure onset-zone (SOZ) lateralizing ability of 7T rs-fMRI and 3T rs-fMRI. Methods: We investigated a cohort of 70 temporal lobe epilepsy (TLE) patients. A paired cohort of 19 patients had 3T and 7T rs-fMRI acquisitions for direct comparison between the two field strengths. Forty-three patients had only 3T, and 8 patients had only 7T rs-fMRI acquisitions. We quantified the functional connectivity between the hippocampus and other nodes within the default mode network (DMN) using seed-to-voxel connectivity, and measured how hippocampo-DMN connectivity could inform SOZ lateralization at 7T and 3T field strengths. Results: Differences between hippocampo-DMN connectivity ipsilateral and contralateral to the SOZ were significantly higher at 7T (pFDR=0.008) than at 3T (pFDR=0.80) when measured in the same subjects. We found that our ability to lateralize the SOZ, by distinguishing subjects with left TLE from subjects with right TLE, was superior at 7T (AUC = 0.97) than 3T (AUC = 0.68). Our findings were reproduced in extended cohorts of subjects scanned at either 3T or 7T. Our rs-fMRI findings at 7T, but not 3T, are consistent and highly correlated (Spearman Rho=0.65) with clinical FDG-PET lateralizing hypometabolism. Significance: We show superior SOZ lateralization in TLE patients when using 7T relative to 3T rs-fMRI, supporting the adoption of high-field strength functional imaging in the epilepsy presurgical evaluation.

Subcortical functional connectivity gradients in temporal lobe epilepsy.

BACKGROUND AND MOTIVATION: Functional gradients have been used to study differences in connectivity between healthy and diseased brain states, however this work has largely focused on the cortex. Because the subcortex plays a key role in seizure initiation in temporal lobe epilepsy (TLE), subcortical functional-connectivity gradients may help further elucidate differences between healthy brains and TLE, as well as differences between left (L)-TLE and right (R)-TLE. METHODS: In this work, we calculated subcortical functional-connectivity gradients (SFGs) from resting-state functional MRI (rs-fMRI) by measuring the similarity in connectivity profiles of subcortical voxels to cortical gray matter voxels. We performed this analysis in 24 R-TLE patients and 31 L-TLE patients (who were otherwise matched for age, gender, disease specific characteristics, and other clinical variables), and 16 controls. To measure differences in SFGs between L-TLE and R-TLE, we quantified deviations in the average functional gradient distributions, as well as their variance, across subcortical structures. RESULTS: We found an expansion, measured by increased variance, in the principal SFG of TLE relative to controls. When comparing the gradient across subcortical structures between L-TLE and R-TLE, we found that abnormalities in the ipsilateral hippocampal gradient distributions were significantly different between L-TLE and R-TLE. CONCLUSION: Our results suggest that expansion of the SFG is characteristic of TLE. Subcortical functional gradient differences exist between left and right TLE and are driven by connectivity changes in the hippocampus ipsilateral to the seizure onset zone.

iEEG-recon: A Fast and Scalable Pipeline for Accurate Reconstruction of Intracranial Electrodes and Implantable Devices.

Background: Collaboration between epilepsy centers is essential to integrate multimodal data for epilepsy research. Scalable tools for rapid and reproducible data analysis facilitate multicenter data integration and harmonization. Clinicians use intracranial EEG (iEEG) in conjunction with non-invasive brain imaging to identify epileptic networks and target therapy for drug-resistant epilepsy cases. Our goal was to promote ongoing and future collaboration by automating the process of "electrode reconstruction," which involves the labeling, registration, and assignment of iEEG electrode coordinates on neuroimaging. These tasks are still performed manually in many epilepsy centers. We developed a standalone, modular pipeline that performs electrode reconstruction. We demonstrate our tool's compatibility with clinical and research workflows and its scalability on cloud platforms. Methods: We created iEEG-recon, a scalable electrode reconstruction pipeline for semi-automatic iEEG annotation, rapid image registration, and electrode assignment on brain MRIs. Its modular architecture includes three modules: a clinical module for electrode labeling and localization, and a research module for automated data processing and electrode contact assignment. To ensure accessibility for users with limited programming and imaging expertise, we packaged iEEG-recon in a containerized format that allows integration into clinical workflows. We propose a cloud-based implementation of iEEG-recon, and test our pipeline on data from 132 patients at two epilepsy centers using retrospective and prospective cohorts. Results: We used iEEG-recon to accurately reconstruct electrodes in both electrocorticography (ECoG) and stereoelectroencephalography (SEEG) cases with a 10 minute running time per case, and ~20 min for semi-automatic electrode labeling. iEEG-recon generates quality assurance reports and visualizations to support epilepsy surgery discussions. Reconstruction outputs from the clinical module were radiologically validated through pre- and post-implant T1-MRI visual inspections. Our use of ANTsPyNet deep learning approach for brain segmentation and electrode classification was consistent with the widely used Freesurfer segmentation. Discussion: iEEG-recon is a valuable tool for automating reconstruction of iEEG electrodes and implantable devices on brain MRI, promoting efficient data analysis, and integration into clinical workflows. The tool's accuracy, speed, and compatibility with cloud platforms make it a useful resource for epilepsy centers worldwide. Comprehensive documentation is available at https://ieeg-recon.readthedocs.io/en/latest/.

Quantifying interictal intracranial EEG to predict focal epilepsy.

Introduction: Intracranial EEG (IEEG) is used for 2 main purposes, to determine: (1) if epileptic networks are amenable to focal treatment and (2) where to intervene. Currently these questions are answered qualitatively and sometimes differently across centers. There is a need for objective, standardized methods to guide surgical decision making and to enable large scale data analysis across centers and prospective clinical trials. Methods: We analyzed interictal data from 101 patients with drug resistant epilepsy who underwent presurgical evaluation with IEEG. We chose interictal data because of its potential to reduce the morbidity and cost associated with ictal recording. 65 patients had unifocal seizure onset on IEEG, and 36 were non-focal or multi-focal. We quantified the spatial dispersion of implanted electrodes and interictal IEEG abnormalities for each patient. We compared these measures against the "5 Sense Score (5SS)," a pre-implant estimate of the likelihood of focal seizure onset, and assessed their ability to predict the clinicians' choice of therapeutic intervention and the patient outcome. Results: The spatial dispersion of IEEG electrodes predicted network focality with precision similar to the 5SS (AUC = 0.67), indicating that electrode placement accurately reflected pre-implant information. A cross-validated model combining the 5SS and the spatial dispersion of interictal IEEG abnormalities significantly improved this prediction (AUC = 0.79; p<0.05). The combined model predicted ultimate treatment strategy (surgery vs. device) with an AUC of 0.81 and post-surgical outcome at 2 years with an AUC of 0.70. The 5SS, interictal IEEG, and electrode placement were not correlated and provided complementary information. Conclusions: Quantitative, interictal IEEG significantly improved upon pre-implant estimates of network focality and predicted treatment with precision approaching that of clinical experts. We present this study as an important step in building standardized, quantitative tools to guide epilepsy surgery.

Resting-state background features demonstrate multidien cycles in long-term EEG device recordings.

Background: Longitudinal EEG recorded by implanted devices is critical for understanding and managing epilepsy. Recent research reports patient-specific, multi-day cycles in device-detected epileptiform events that coincide with increased likelihood of clinical seizures. Understanding these cycles could elucidate mechanisms generating seizures and advance drug and neurostimulation therapies. Objective/Hypothesis: We hypothesize that seizure-correlated cycles are present in background neural activity, independent of interictal epileptiform spikes, and that neurostimulation may disrupt these cycles. Methods: We analyzed regularly-recorded seizure-free data epochs from 20 patients implanted with a responsive neurostimulation (RNS) device for at least 1.5 years, to explore the relationship between cycles in device-detected interictal epileptiform activity (dIEA), clinician-validated interictal spikes, background EEG features, and neurostimulation. Results: Background EEG features tracked the cycle phase of dIEA in all patients (AUC: 0.63 [0.56 - 0.67]) with a greater effect size compared to clinically annotated spike rate alone (AUC: 0.55 [0.53-0.61], p < 0.01). After accounting for circadian variation and spike rate, we observed significant population trends in elevated theta and beta band power and theta and alpha connectivity features at the cycle peaks (sign test, p < 0.05). In the period directly after stimulation we observe a decreased association between cycle phase and EEG features compared to background recordings (AUC: 0.58 [0.55-0.64]). Conclusions: Our findings suggest that seizure-correlated dIEA cycles are not solely due to epileptiform discharges but are associated with background measures of brain state; and that neurostimulation may disrupt these cycles. These results may help elucidate mechanisms underlying seizure generation, provide new biomarkers for seizure risk, and facilitate monitoring, treating, and managing epilepsy with implantable devices.