Altered metabolic-functional coupling in the epileptogenic network could predict surgical outcomes of mesial temporal lobe epilepsy

Objective

To investigate the relationship between glucose metabolism and functional activity in the epileptogenic network of patients with mesial temporal lobe epilepsy (MTLE) and to determine whether this relationship is associated with surgical outcomes.

Methods

¹⁸F-FDG PET and resting-state functional MRI (rs-fMRI) scans were performed on a hybrid PET/MR scanner in 38 MTLE patients with hippocampal sclerosis (MR-HS), 35 MR-negative patients and 34 healthy controls (HC). Glucose metabolism was measured using ¹⁸F-FDG PET standardized uptake value ratio (SUVR) relative to cerebellum; Functional activity was obtained by fractional amplitude of low-frequency fluctuation (fALFF). The betweenness centrality (BC) of metabolic covariance network and functional network were calculated using graph theoretical analysis. Differences in SUVR, fALFF, BC and the spatial voxel-wise SUVR-fALFF couplings of the epileptogenic network, consisting of default mode network (DMN) and thalamus, were evaluated by Mann-Whitney U test (using the false discovery rate [FDR] for multiple comparison correction). The top ten SUVR-fALFF couplings were selected by Fisher score to predict surgical outcomes using logistic regression model.

Results

The results showed decreased SUVR-fALFF coupling in the bilateral middle frontal gyrus (P_FDR = 0.0230, P_FDR = 0.0296) in MR-HS patients compared to healthy controls. Coupling in the ipsilateral hippocampus was marginally increased (P_FDR = 0.0802) in MR-HS patients along with decreased BC of metabolic covariance network and functional network (P_FDR = 0.0152; P_FDR = 0.0429). With Fisher score ranking, the top ten SUVR-fALFF couplings in regions from DMN and thalamic subnuclei could predict surgical outcomes with the best performance being a combination of ten SUVR-fALFF couplings with an AUC of 0.914.

Conclusion

These findings suggest that the altered neuroenergetic coupling in the epileptogenic network is associated with surgical outcomes of MTLE patients, which may provide insight into their pathogenesis and help with preoperative evaluation.

Neurovascular coupling changes in patients with magnetic resonance imaging negative focal epilepsy

Brain neuron activity is closely related to cerebral blood flow (CBF) changes. Alterations in the regional homogeneity (ReHo) and CBF occur in patients with magnetic resonance imaging negative focal epilepsy (FEP-MRI^-). However, the coupling alterations of ReHo and CBF in FEP-MRI^- remain unclear. The study aims to explore neurovascular coupling alterations and their clinical implication in FEP-MRI^-. We collected resting-state magnetic resonance imaging (MRI) data from 31 healthy controls (HCs) and 48 patients with FEP-MRI^-,including three-dimensional (3D) T1-weighted imaging, 3D arterial spin labeling (ASL) imaging,and resting-state functional MRI (rs-fMRI). The CBF and ReHo values were calculated from the ASL and rs-fMRI data, respectively. The CBF/ReHo ratio per voxel and whole-brain CBF-ReHo coupling were compared between the two groups. Correlation analysis involved the CBF/ReHo ratio and clinical indicators in FEP-MRI^-. Patients with FEP-MRI^- showed significantly increased cross-subject CBF-ReHo and global cross-voxel CBF-ReHo coupling. The CBF/ReHo ratio was higher in the bilateral orbitofrontal gyrus, right parietal lobe, and right middle frontal gyrus of patients with FEP-MRI^-. Nevertheless, this ratio was lower in the bilateral supplementary motor areas, the left middle and posterior cingulate gyrus, and the right central sulcus cover. The CBF/ReHo ratio was markedly correlated with cognitive function, memory, intelligence, and epilepsy duration in the above abnormal brain regions. CBF/ReHo ratio may be useful as an indicator of neuropathological mechanisms. These results support the hypothesis that CBF/ReHo ratio relates to the neuropathological mechanisms of FEP-MRI^-. Furthermore, it offers new perspectives for studying the mechanisms of MRI-negative epilepsy.

Diffusion kurtosis imaging and diffusion tensor imaging parameters applied to white matter and gray matter of patients with anti-N-methyl-D-aspartate receptor encephalitis

Objectives

To compare parameters of diffusion tensor imaging (DTI) and diffusion kurtosis imaging (DKI) to evaluate which can better describe the microstructural changes of anti-N-methyl-D-aspartate receptor (NMDAR) encephalitis patients and to characterize the non-Gaussian diffusion patterns of the whole brain and their correlation with neuropsychological impairments in these patients.

Materials and methods

DTI and DKI parameters were measured in 57 patients with anti-NMDAR encephalitis and 42 healthy controls. Voxel-based analysis was used to evaluate group differences between white matter and gray matter separately. The modified Rankin Scale (mRS) was used to evaluate the severity of the neurofunctional recovery of patients, the Montreal Cognitive Assessment (MoCA) was used to assess global cognitive performance, and the Hamilton Depression Scale (HAMD) and fatigue severity scale (FSS) were used to evaluate depressive and fatigue states.

Results

Patients with anti-NMDAR encephalitis showed significantly decreased radial kurtosis (RK) in the right extranucleus in white matter (P < 0.001) and notably decreased kurtosis fractional anisotropy (KFA) in the right precuneus, the right superior parietal gyrus (SPG), the left precuneus, left middle occipital gyrus, and left superior occipital gyrus in gray matter (P < 0.001). Gray matter regions with decreased KFA overlapped with those with decreased RK in the left middle temporal gyrus, superior temporal gyrus (STG), supramarginal gyrus (SMG), postcentral gyrus (POCG), inferior parietal but supramarginal gyrus, angular gyrus (IPL) and angular gyrus (ANG) (P < 0.001). The KFA and RK in the left ANG, IPL and POCG correlated positively with MoCA scores. KFA and RK in the left ANG, IPL, POCG and SMG correlated negatively with mRS scores. KFA in the left precuneus and right SPG as well as RK in the left STG correlated negatively with mRS scores. No significant correlation between KFA and RK in the abnormal brain regions and HAMD and FSS scores was found.

Conclusion

The microstructural changes in gray matter were much more extensive than those in white matter in patients with anti-NMDAR encephalitis. The brain damage reflected by DKI parameters, which have higher sensitivity than parameters of DTI, correlated with cognitive impairment and the severity of the neurofunctional recovery.

Long-term development of dynamic changes in neurovascular coupling after acute temporal lobe epilepsy

Epilepsy is an abnormal brain state that may be induced by synchronous neuronal activation and also abnormalities in energy metabolism or the oxygen supply vascular system. Neurovascular coupling (NVC), the relationship between neuron, capillary, and penetrating artery, remains unexplored on a fine-scale with respect to the pathology process after acute temporal lobe epilepsy (TLE). Here we use two-photon microscopy (TPM) to provide high temporal-spatial resolution imaging to identify changes in NVC during spontaneous and electro-stimulated (ES) states in awake mice. Implantation of a long-term craniotomy window allowed TPM recording of the pathological development after the acute Kainic Acid temporal lobe epilepsy model. Our results provide direct evidence that the capillary and penetrating artery are not correlated to rhythmic neuronal activity during acute epilepsy. During the CSD period, NVC shows a strong correlation. We demonstrate that NVC exhibits nonlinear dynamics after status epilepticus. Furthermore, the vascular correlation to neuronal signals in spontaneous and ES states shows dynamic changes which correlate to the evolution after acute TLE. Understanding NVC in all TLE stages, from the acute through the TLE pathological development, may provide new therapeutic pathways.

Arterial blood stealing as a mechanism of negative BOLD response: From the steady-flow with nonlinear phase separation to a windkessel-based model

Blood Oxygen Level Dependent (BOLD) signal indirectly characterizes neuronal activity by measuring hemodynamic and metabolic changes in the nearby microvasculature. A deeper understanding of how localized changes in electrical, metabolic and hemodynamic factors translate into a BOLD signal is crucial for the interpretation of functional brain imaging techniques. While positive BOLD responses (PBR) are widely considered to be linked with neuronal activation, the origins of negative BOLD responses (NBR) have remained largely unknown. As NBRs are sometimes observed in close proximity of regions with PBR, a blood "stealing" effect, i.e., redirection of blood from a passive periphery to the area with high neuronal activity, has been postulated. In this study, we used the Hagen-Poiseuille equation to model hemodynamics in an idealized microvascular network that account for the particulate nature of blood and nonlinearities arising from the red blood cell (RBC) distribution (i.e., the Fåhraeus, Fåhraeus-Lindqvist and the phase separation effects). Using this detailed model, we evaluate determinants driving this "stealing" effect in a microvascular network with geometric parameters within physiological ranges. Model simulations predict that during localized cerebral blood flow (CBF) increases due to neuronal activation-hyperemic response, blood from surrounding vessels is reallocated towards the activated region. This stealing effect depended on the resistance of the microvasculature and the uneven distribution of RBCs at vessel bifurcations. A parsimonious model consisting of two-connected windkessel regions sharing a supplying artery was proposed to simulate the stealing effect with a minimum number of parameters. Comparison with the detailed model showed that the parsimonious model can reproduce the observed response for hematocrit values within the physiological range for different species. Our novel parsimonious model promise to be of use for statistical inference (top-down analysis) from direct blood flow measurements (e.g., arterial spin labeling and laser Doppler/Speckle flowmetry), and when combined with theoretical models for oxygen extraction/diffusion will help account for some types of NBRs.

Impaired neurovascular coupling and cognitive deficits in anti-N-methyl-D-aspartate receptor encephalitis

Anti-N-methyl-D-aspartate receptor (anti-NMDAR) encephalitis is a recently identified autoimmune disorder with heterogeneous neurological, psychiatric, and cognitive manifestations. The NMDAR is a key signaling node for neurovascular coupling, the mechanism by which cerebral blood perfusion is enhanced to meet local metabolic requirements from increased neuronal activity. Therefore, anti-NMDAR encephalitis may disrupt neurovascular coupling and induce cognitive deficits. This study examined neurovascular coupling and cognitive function in anti-NMDAR encephalitis patients to identify prognostic biomarkers, reveal potential pathogenic mechanisms, and provide clues to possible therapeutic strategies. In this study, twenty-three anti-NMDAR encephalitis patients and thirty healthy controls received neuropsychological testing and multimodal magnetic resonance imaging (MRI). Cerebral blood flow (CBF) was calculated from arterial spin labeling, and regional homogeneity (ReHo) was computed from functional MRI. Pearson's correlation coefficients between CBF and ReHo were calculated to obtain neurovascular coupling. At the whole gray matter level, CBF‒ReHo coupling was reduced in patients compared to healthy controls. At the regional level, CBF‒ReHo was significantly lower among patients in the precentral gyrus, frontal gyrus, insula, cuneus, inferior parietal lobe, supramarginal gyrus, angular gyrus, precuneus, temporal gyrus, and temporal pole. Reduced CBF‒ReHo in the left superior medial frontal gyrus of patients was significantly correlated with a deficit in verbal inhibition control, and the reduced CBF‒ReHo in the left insula was significantly correlated with impaired executive function. In conclusion, anti-NMDAR encephalitis is associated with both global and regional disruptions in neurovascular coupling that may in turn lead to deficits in specific cognitive domains.

Identification of Negative BOLD Responses in Epilepsy Using Windkessel Models

Alongside positive blood oxygenation level–dependent (BOLD) responses associated with interictal epileptic discharges, a variety of negative BOLD responses (NBRs) are typically found in epileptic patients. Previous studies suggest that, in general, up to four mechanisms might underlie the genesis of NBRs in the brain: (i) neuronal disruption of network activity, (ii) altered balance of neurometabolic/vascular couplings, (iii) arterial blood stealing, and (iv) enhanced cortical inhibition. Detecting and classifying these mechanisms from BOLD signals are pivotal for the improvement of the specificity of the electroencephalography–functional magnetic resonance imaging (EEG-fMRI) image modality to identify the seizure-onset zones in refractory local epilepsy. This requires models with physiological interpretation that furnish the understanding of how these mechanisms are fingerprinted by their BOLD responses. Here, we used a Windkessel model with viscoelastic compliance/inductance in combination with dynamic models of both neuronal population activity and tissue/blood O₂ to classify the hemodynamic response functions (HRFs) linked to the above mechanisms in the irritative zones of epileptic patients. First, we evaluated the most relevant imprints on the BOLD response caused by variations of key model parameters. Second, we demonstrated that a general linear model is enough to accurately represent the four different types of NBRs. Third, we tested the ability of a machine learning classifier, built from a simulated ensemble of HRFs, to predict the mechanism underlying the BOLD signal from irritative zones. Cross-validation indicates that these four mechanisms can be classified from realistic fMRI BOLD signals. To demonstrate proof of concept, we applied our methodology to EEG-fMRI data from five epileptic patients undergoing neurosurgery, suggesting the presence of some of these mechanisms. We concluded that a proper identification and interpretation of NBR mechanisms in epilepsy can be performed by combining general linear models and biophysically inspired models.

Mesoscopic Mapping of Ictal Neurovascular Coupling in Awake Behaving Mice Using Optical Spectroscopy and Genetically Encoded Calcium Indicators

Unambiguously identifying an epileptic focus with high spatial resolution is a challenge, especially when no anatomic abnormality can be detected. Neurovascular coupling (NVC)-based brain mapping techniques are often applied in the clinic despite a poor understanding of ictal NVC mechanisms, derived primarily from recordings in anesthetized animals with limited spatial sampling of the ictal core. In this study, we used simultaneous wide-field mesoscopic imaging of GCamp6f and intrinsic optical signals (IOS) to record the neuronal and hemodynamic changes during acute ictal events in awake, behaving mice. Similar signals in isoflurane-anesthetized mice were compared to highlight the unique characteristics of the awake condition. In awake animals, seizures were more focal at the onset but more likely to propagate to the contralateral hemisphere. The HbT signal, derived from an increase in cerebral blood volume (CBV), was more intense in awake mice. As a result, the “epileptic dip” in hemoglobin oxygenation became inconsistent and unreliable as a mapping signal. Our data indicate that CBV-based imaging techniques should be more accurate than blood oxygen level dependent (BOLD)-based imaging techniques for seizure mapping in awake behaving animals.

Autonomic Dysfunction Contributes to Impairment of Cerebral Autoregulation in Patients with Epilepsy

Patients with epilepsy frequently experience autonomic dysfunction and impaired cerebral autoregulation. The present study investigates autonomic function and cerebral autoregulation in patients with epilepsy to determine whether these factors contribute to impaired autoregulation. A total of 81 patients with epilepsy and 45 healthy controls were evaluated, assessing their sudomotor, cardiovagal, and adrenergic functions using a battery of autonomic nervous system (ANS) function tests, including the deep breathing, Valsalva maneuver, head-up tilting, and Q-sweat tests. Cerebral autoregulation was measured by transcranial Doppler examination during the breath-holding test, the Valsalva maneuver, and the head-up tilting test. Autonomic functions were impaired during the interictal period in patients with epilepsy compared to healthy controls. The three indices of cerebral autoregulation—the breath-holding index (BHI), an autoregulation index calculated in phase II of the Valsalva maneuver (ASI), and cerebrovascular resistance measured in the second minute during the head-up tilting test (CVR_2-min)—all decreased in patients with epilepsy. ANS dysfunction correlated significantly with impairment of cerebral autoregulation (measured by BHI, ASI, and CVR_2-min), suggesting that the increased autonomic dysfunction in patients with epilepsy may augment the dysregulation of cerebral blood flow. Long-term epilepsy, a high frequency of seizures, and refractory epilepsy, particularly temporal lobe epilepsy, may contribute to advanced autonomic dysfunction and impaired cerebral autoregulation. These results have implications for therapeutic interventions that aim to correct central autonomic dysfunction and impairment of cerebral autoregulation, particularly in patients at high risk for sudden, unexplained death in epilepsy.

Altered coupling between resting-state glucose metabolism and functional activity in epilepsy

OBJECTIVE

Altered functional activities and hypometabolism have been found in medial temporal lobe epilepsy patients with hippocampal sclerosis (mTLE-HS). Hybrid PET/MR scanners provide opportunities to explore the relationship between resting-state energy consumption and functional activities, but whether repeated seizures disturb the bioenergetic coupling and its relationship with seizure outcomes remain unknown.

METHODS

¹⁸ F-FDG PET and resting-state functional MRI (rs-fMRI) scans were performed with hybrid PET/MR in 26 patients with mTLE-HS and in healthy controls. Energy consumption was quantified by ¹⁸ F-FDG standardized uptake value ratio(SUVR) relative to cerebellum. Spontaneous neural activities were estimated using regional homogeneity (ReHo), fractional amplitude of low frequency fluctuations (fALFF) from rs-fMRI. Between-group differences in SUVR and rs-fMRI derived metrics were evaluated by two-sample t test. Voxel-wise spatial correlations were explored between SUVR and ReHo, fALFF across gray matter and compared between groups. Furthermore, the relationships between altered fALFF/SUVR and ReHo/SUVR coupling and surgical outcomes were evaluated.

RESULTS

Both the patients and healthy controls showed significant positive correlations between SUVR and rs-fMRI metrics. Spatial correlations between SUVR and fMRI-derived metrics across gray matter were significantly higher in patients with mTLE-HS compared with healthy controls (fALFF/SUVR, P < 0.001; ReHo/SUVR, P = 0.022). Higher fALFF/SUVR couplings were found in patients who had Engel class IA after surgery than all other (P = 0.025), while altered ReHo/SUVR couplings (P = 0.097) were not.

CONCLUSION

These findings demonstrated altered bioenergetic coupling across gray matter and its relationship with seizure outcomes, which may provide novel insights into pathogenesis of mTLE-HS and potential biomarkers for epilepsy surgery planning.

Neuropsychological Deficits in Patients with Electrical Status Epilepticus During Sleep: A Non-invasive Analysis of Neurovascular Coupling

To evaluate the effects of electrical status epilepticus during sleep (ESES) on cerebral blood flow (CBF) and explore the associated neuro-vascular coupling and neuropsychological deficits. 19 ESES patients were recruited to undergo real-time transcranial doppler ultrasonography (TCD) and video-EEG monitoring (vEEG). Patients were grouped based on their cognitive functions or their EEG patterns. The mean cerebral blood flow velocity (CBFV_m) of the unilateral middle cerebral artery was measured using TCD and was used to calculate various relevant parameters. The 19 patients participated in a total of 54 effective TCD-vEEG monitoring sessions. We found a significant effect of clinical severity for the following measurements: spike wave index (SWI), peak and average deep sleep stage (N3) CBFV_m, peak, average and minimum deep sleep and awake CBFV_m, and CBFV_m oscillations during deep sleep. Nevertheless, CBFV_m oscillations were not related to SWI. Furthermore, CBFV_m oscillations revealed a statistically significant difference between the near-ESES and asymmetric-ESES groups. CBFV_m oscillations may reflect the neuro-vascular coupling process associated with ESES disfunction. Understanding the relationship between CBFV_m oscillations and epileptic activity will be important for assessing the neuropsychological damage associated with ESES and for developing treatment options for this and other diseases.

waveCSD: A method for estimating transmembrane currents originated from propagating neuronal activity in the neocortex: Application to study cortical spreading depression

BACKGROUND

Recent years have witnessed an upsurge in the development of methods for estimating current source densities (CSDs) in the neocortical tissue from their recorded local field potential (LFP) reflections using microelectrode arrays. Among these, methods utilizing linear arrays work under the assumption that CSDs vary as a function of cortical depth; whereas they are constant in the tangential direction, infinitely or in a confined cylinder. This assumption is violated in the analysis of neuronal activity propagating along the neocortical sheet, e.g. propagation of alpha waves or cortical spreading depression.

NEW METHOD

Here, we developed a novel mathematical method (waveCSD) for CSD analysis of LFPs associated with a planar wave of neocortical neuronal activity propagating at a constant velocity towards a linear probe.

RESULTS

Results show that the algorithm is robust to the presence of noise in LFP data and uncertainties in knowledge of propagation velocity. Also, results show high level of accuracy of the method in a wide range of electrode resolutions. Using in vivo experimental recordings from the rat neocortex, we employed waveCSD to characterize transmembrane currents associated with cortical spreading depressions.

COMPARISON WITH EXISTING METHOD(S)

Simulation results indicate that waveCSD has a significantly higher reconstruction accuracy compared to the widely-used inverse CSD method (iCSD), and the regularized kernel CSD method (kCSD), in the analysis of CSDs originating from propagating neuronal activity.

CONCLUSIONS

The waveCSD method provides a theoretical platform for estimation of transmembrane currents from their LFPs in experimental paradigms involving wave propagation.

Histological Characterization of the Irritative Zones in Focal Cortical Dysplasia Using a Preclinical Rat Model

Current clinical practice in focal epilepsy involves brain source imaging (BSI) to localize brain areas where from interictal epileptiform discharges (IEDs) emerge. These areas, named irritative zones, have been useful to define candidate seizures-onset zones during pre-surgical workup. Since human histological data are mostly available from final resected zones, systematic studies characterizing pathophysiological mechanisms and abnormal molecular/cellular substrates in irritative zones—independent of them being epileptogenic—are challenging. Combining BSI and histological analysis from all types of irritative zones is only possible through the use of preclinical animal models. Here, we recorded 32-channel spontaneous electroencephalographic data from rats that have focal cortical dysplasia (FCD) and chronic seizures. BSI for different IED subtypes was performed using the methodology presented in Bae et al. (2015). Post-mortem brain sections containing irritative zones were stained to quantify anatomical, functional, and inflammatory biomarkers specific for epileptogenesis, and the results were compared with those obtained using the contralateral healthy brain tissue. We found abnormal anatomical structures in all irritative zones (i.e., larger neuronal processes, glioreactivity, and vascular cuffing) and larger expressions for neurotransmission (i.e., NR2B) and inflammation (i.e., ILβ1, TNFα and HMGB1). We conclude that irritative zones in this rat preclinical model of FCD comprise abnormal tissues disregarding whether they are actually involved in icto-genesis or not. We hypothesize that seizure perpetuation happens gradually; hence, our results could support the use of IED-based BSI for the early diagnosis and preventive treatment of potential epileptic foci. Further verifications in humans are yet needed.

Quantitative assessment of hemodynamic and structural characteristics of in vivo brain tissue using total diffuse reflectance spectrum measured in a non-contact fashion

Here we present a new methodology that investigates the intrinsic structural and hemodynamic characteristics of in vivo brain tissue, in a non-contact fashion, and can be easily incorporated in an intra-operative environment. Within this methodology, relative total diffuse reflectance spectra (R_TD(λ)) were acquired from targets using a hybrid spectroscopy imaging system. A spectral interpretation algorithm was subsequently applied to R_TD(λ) to retrieve optical properties related to the compositional and structural characteristics of each target. Estimation errors of the proposed methodology were computationally evaluated using a Monte Carlo simulation model for photon migration under various conditions. It was discovered that this new methodology could handle moderate noise and achieve very high accuracy, but only if the refractive index of the target is known. The accuracy of the technique was also validated using a series of tissue phantom studies, and consistent and accurate estimates of μ_s'(λ)/μ_a(λ) were obtained from all the phantoms tested. Finally, a small-scale animal study was conducted to demonstrate the clinical utility of the reported method, wherein a forepaw stimulation model was utilized to induce transient hemodynamic responses in somatosensory cortices. With this approach, significant stimulation-related changes (p < 0.001) in cortical hemodynamic and structural characteristics were successfully measured.