Investigation of Linear Coupling Between Single-Event Blood Flow Responses and Interictal Discharges in a Model of Experimental Epilepsy

Preictal neuronal and vascular activity precedes the onset of childhood absence seizure: direct current potential shifts and their correlation with hemodynamic activity

SIGNIFICANCE AIMS

The neurovascular mechanisms underlying the initiation of absence seizures and their dynamics are still not well understood. The objective of this study was to better noninvasively characterize the dynamics of the neuronal and vascular network at the transition from the interictal state to the ictal state of absence seizures and back to the interictal state using a combined electroencephalography (EEG), functional near-infrared spectroscopy (fNIRS), and diffuse correlation spectroscopy (DCS) approach. The second objective was to develop hypotheses about the neuronal and vascular mechanisms that propel the networks to the 3-Hz spikes and wave discharges (SWDs) observed during absence seizures.

APPROACHES

We evaluated the simultaneous changes in electrical (neuronal) and optical dynamics [hemodynamic, with changes in (Hb) and cerebral blood flow] of 8 pediatric patients experiencing 25 typical childhood absence seizures during the transition from the interictal state to the absence seizure by simultaneously performing EEG, fNIRS, and DCS.

RESULTS

Starting from ∼ 20    s before the onset of the SWD, we observed a transient direct current potential shift that correlated with alterations in functional fNIRS and DCS measurements of the cerebral hemodynamics detecting the preictal changes.

DISCUSSION

Our noninvasive multimodal approach highlights the dynamic interactions between the neuronal and vascular compartments that take place in the neuronal network near the time of the onset of absence seizures in a very specific cerebral hemodynamic environment. These noninvasive approaches contribute to a better understanding of the electrical hemodynamic environment prior to seizure onset. Whether this may ultimately be relevant for diagnostic and therapeutic approaches requires further evaluation.

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.

A novel dynamic network imaging analysis method reveals aging-related fragmentation of cortical networks in mouse

Network analysis of large-scale neuroimaging data is a particularly challenging computational problem. Here, we adapt a novel analytical tool, the community dynamic inference method (CommDy), for brain imaging data from young and aged mice. CommDy, which was inspired by social network theory, has been successfully used in other domains in biology; this report represents its first use in neuroscience. We used CommDy to investigate aging-related changes in network metrics in the auditory and motor cortices by using flavoprotein autofluorescence imaging in brain slices and in vivo. We observed that auditory cortical networks in slices taken from aged brains were highly fragmented compared to networks observed in young animals. CommDy network metrics were then used to build a random-forests classifier based on NMDA receptor blockade data, which successfully reproduced the aging findings, suggesting that the excitatory cortical connections may be altered during aging. A similar aging-related decline in network connectivity was also observed in spontaneous activity in the awake motor cortex, suggesting that the findings in the auditory cortex reflect general mechanisms during aging. These data suggest that CommDy provides a new dynamic network analytical tool to study the brain and that aging is associated with fragmentation of intracortical networks.

What Triggers the Interictal Epileptic Spike? A Multimodal Multiscale Analysis of the Dynamic of Synaptic and Non-synaptic Neuronal and Vascular Compartments Using Electrical and Optical Measurements

Interictal spikes (IISs) may result from a disturbance of the intimate functional balance between various neuronal (synaptic and non-synaptic), vascular, and metabolic compartments. To better characterize the complex interactions within these compartments at different scales we developed a simultaneous multimodal-multiscale approach and measure their activity around the time of the IIS. We performed such measurements in an epileptic rat model (n = 43). We thus evaluated (1) synaptic dynamics by combining electrocorticography and multiunit activity recording in the time and time-frequency domain, (2) non-synaptic dynamics by recording modifications in light scattering induced by changes in the membrane configuration related to cell activity using the fast optical signal, and (3) vascular dynamics using functional near-infrared spectroscopy and, independently but simultaneously to the electrocorticography, the changes in cerebral blood flow using diffuse correlation spectroscopy. The first observed alterations in the measured signals occurred in the hemodynamic compartments a few seconds before the peak of the IIS. These hemodynamic changes were followed by changes in coherence and then synchronization between the deep and superficial neural networks in the 1 s preceding the IIS peaks. Finally, changes in light scattering before the epileptic spikes suggest a change in membrane configuration before the IIS. Our multimodal, multiscale approach highlights the complexity of (1) interactions between the various neuronal, vascular, and extracellular compartments, (2) neural interactions between various layers, (3) the synaptic mechanisms (coherence and synchronization), and (4) non-synaptic mechanisms that take place in the neuronal network around the time of the IISs in a very specific cerebral hemodynamic environment.

Dynamic oxygen changes during status epilepticus and subsequent endogenous kindling

OBJECTIVE

Brain tissue oxygen (partial oxygen pressure [pO2 ]) levels are tightly regulated to stay within the normoxic zone, with deviations on either side resulting in impaired brain function. Whereas pathological events such as ischemic attacks and brief seizures have previously been shown to result in pO2 levels well below the normoxic zone, oxygen levels during prolonged status epilepticus (SE) and the subsequent endogenous kindling period are unknown.

METHODS

We utilized two models of acquired temporal lobe epilepsy in rats: intrahippocampal kainic acid infusion and prolonged perforant pathway stimulation. Local tissue oxygen was measured in the dorsal hippocampus using an optode during and for several weeks following SE.

RESULTS

We observed hyperoxia in the hippocampus during induced SE in both models. Following termination of SE, 88% of rats initiated focal self-generated spiking activity in the hippocampus within the first 7 days, which was associated with dynamic oxygen changes. Self-generated and recurring epileptiform activity subsequently organized into higher-frequency bursts that became progressively longer and were ultimately associated with behavioral seizures that became more severe with time and led to postictal hypoxia.

SIGNIFICANCE

Induced SE and self-generated recurrent epileptiform activity can have profound and opposing effects on brain tissue oxygenation that may serve as a biomarker for ongoing pathological activity in the brain.

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.

Effect of temperature on FAD and NADH-derived signals and neurometabolic coupling in the mouse auditory and motor cortex

Tight coupling of neuronal metabolism to synaptic activity is critical to ensure that the supply of metabolic substrates meets the demands of neuronal signaling. Given the impact of temperature on metabolism, and the wide fluctuations of brain temperature observed during clinical hypothermia, we examined the effect of temperature on neurometabolic coupling. Intrinsic fluorescence signals of the oxidized form of flavin adenine dinucleotide (FAD) and the reduced form of nicotinamide adenine dinucleotide (NADH), and their ratios, were measured to assess neural metabolic state and local field potentials were recorded to measure synaptic activity in the mouse brain. Brain slice preparations were used to remove the potential impacts of blood flow. Tight coupling between metabolic signals and local field potential amplitudes was observed at a range of temperatures below 29 °C. However, above 29 °C, the metabolic and synaptic signatures diverged such that FAD signals were diminished, but local field potentials retained their amplitude. It was also observed that the declines in the FAD signals seen at high temperatures (and hence the decoupling between synaptic and metabolic events) are driven by low FAD availability at high temperatures. These data suggest that neurometabolic coupling, thought to be critical for ensuring the metabolic health of the brain, may show temperature dependence, and is related to temperature-dependent changes in FAD supplies.

A study of the electro-haemodynamic coupling using simultaneously acquired intracranial EEG and fMRI data in humans

In current fMRI studies designed to map BOLD changes related to interictal epileptiform discharges (IED), which are recorded on simultaneous EEG, the information contained in the morphology and field extent of the EEG events is exclusively used for their classification. Usually, a BOLD predictor based on IED onset times alone is constructed, effectively treating all events as identical. We used intracranial EEG (icEEG)-fMRI data simultaneously recorded in humans to investigate the effect of including any of the features: amplitude, width (duration), slope of the rising phase, energy (area under the curve), or spatial field extent (number of contacts over which the sharp wave was observed) of the fast wave of the IED (the sharp wave), into the BOLD model, to better understand the neurophysiological origin of sharp wave-related BOLD changes, in the immediate vicinity of the recording contacts. Among the features considered, the width was the only one found to explain a significant amount of additional variance, suggesting that the amplitude of the BOLD signal depends more on the duration of the underlying field potential (reflected in the sharp wave width) than on the degree of neuronal activity synchrony (reflected in the sharp wave amplitude), and, consequently, that including inter-event variations of the sharp wave width in the BOLD signal model may increase the sensitivity of forthcoming EEG-fMRI studies of epileptic activity.

A New Computational Model for Neuro-Glio-Vascular Coupling: Astrocyte Activation Can Explain Cerebral Blood Flow Nonlinear Response to Interictal Events

Developing a clear understanding of the relationship between cerebral blood flow (CBF) response and neuronal activity is of significant importance because CBF increase is essential to the health of neurons, for instance through oxygen supply. This relationship can be investigated by analyzing multimodal (fMRI, PET, laser Doppler…) recordings. However, the important number of intermediate (non-observable) variables involved in the underlying neurovascular coupling makes the discovery of mechanisms all the more difficult from the sole multimodal data. We present a new computational model developed at the population scale (voxel) with physiologically relevant but simple equations to facilitate the interpretation of regional multimodal recordings. This model links neuronal activity to regional CBF dynamics through neuro-glio-vascular coupling. This coupling involves a population of glial cells called astrocytes via their role in neurotransmitter (glutamate and GABA) recycling and their impact on neighboring vessels. In epilepsy, neuronal networks generate epileptiform discharges, leading to variations in astrocytic and CBF dynamics. In this study, we took advantage of these large variations in neuronal activity magnitude to test the capacity of our model to reproduce experimental data. We compared simulations from our model with isolated epileptiform events, which were obtained in vivo by simultaneous local field potential and laser Doppler recordings in rats after local bicuculline injection. We showed a predominant neuronal contribution for low level discharges and a significant astrocytic contribution for higher level discharges. Besides, neuronal contribution to CBF was linear while astrocytic contribution was nonlinear. Results thus indicate that the relationship between neuronal activity and CBF magnitudes can be nonlinear for isolated events and that this nonlinearity is due to astrocytic activity, highlighting the importance of astrocytes in the interpretation of regional recordings.

Interneurons contribute to the hemodynamic/metabolic response to epileptiform discharges

Interpretation of hemodynamic responses in epilepsy is hampered by an incomplete understanding of the underlying neurovascular coupling, especially the contributions of excitation and inhibition. We made simultaneous multimodal recordings of local field potentials (LFPs), firing of individual neurons, blood flow, and oxygen level in the somatosensory cortex of anesthetized rats. Epileptiform discharges induced by bicuculline injections were used to trigger large local events. LFP and blood flow were robustly coupled, as were LFP and tissue oxygen. In a parametric linear model, LFP and the baseline activities of cerebral blood flow and tissue partial oxygen tension contributed significantly to blood flow and oxygen responses. In an analysis of recordings from 402 neurons, blood flow/tissue oxygen correlated with the discharge of putative interneurons but not of principal cells. Our results show that interneuron activity is important in the vascular and metabolic responses during epileptiform discharges.

Using patient-specific hemodynamic response function in epileptic spike analysis of human epilepsy: a study based on EEG-fNIRS

Functional near-infrared spectroscopy (fNIRS) can be combined with electroencephalography (EEG) to continuously monitor the hemodynamic signal evoked by epileptic events such as seizures or interictal epileptiform discharges (IEDs, aka spikes). As estimation methods assuming a canonical shape of the hemodynamic response function (HRF) might not be optimal, we sought to model patient-specific HRF (sHRF) with a simple deconvolution approach for IED-related analysis with EEG-fNIRS data. Furthermore, a quadratic term was added to the model to account for the nonlinearity in the response when IEDs are frequent. Prior to analyzing clinical data, simulations were carried out to show that the HRF was estimable by the proposed deconvolution methods under proper conditions. EEG-fNIRS data of five patients with refractory focal epilepsy were selected due to the presence of frequent clear IEDs and their unambiguous focus localization. For each patient, both the linear sHRF and the nonlinear sHRF were estimated at each channel. Variability of the estimated sHRFs was seen across brain regions and different patients. Compared with the SPM8 canonical HRF (cHRF), including these sHRFs in the general linear model (GLM) analysis led to hemoglobin activations with higher statistical scores as well as larger spatial extents on all five patients. In particular, for patients with frequent IEDs, nonlinear sHRFs were seen to provide higher sensitivity in activation detection than linear sHRFs. These observations support using sHRFs in the analysis of IEDs with EEG-fNIRS data.

Mapping interictal epileptic discharges using mutual information between concurrent EEG and fMRI

OBJECTIVE

The mapping of haemodynamic changes related to interictal epileptic discharges (IED) in simultaneous electroencephalography (EEG) and functional MRI (fMRI) studies is usually carried out by means of EEG-correlated fMRI analyses where the EEG information specifies the model to test on the fMRI signal. The sensitivity and specificity critically depend on the accuracy of EEG detection and the validity of the haemodynamic model. In this study we investigated whether an information theoretic analysis based on the mutual information (MI) between the presence of epileptic activity on EEG and the fMRI data can provide further insights into the haemodynamic changes related to interictal epileptic activity. The important features of MI are that: 1) both recording modalities are treated symmetrically; 2) no requirement for a-priori models for the haemodynamic response function, or assumption of a linear relationship between the spiking activity and BOLD responses, and 3) no parametric model for the type of noise or its probability distribution is necessary for the computation of MI.

METHODS

Fourteen patients with pharmaco-resistant focal epilepsy underwent EEG-fMRI and intracranial EEG and/or surgical resection with positive postoperative outcome (seizure freedom or considerable reduction in seizure frequency) was available in 7/14 patients. We used nonparametric statistical assessment of the MI maps based on a four-dimensional wavelet packet resampling method. The results of MI were compared to the statistical parametric maps obtained with two conventional General Linear Model (GLM) analyses based on the informed basis set (canonical HRF and its temporal and dispersion derivatives) and the Finite Impulse Response (FIR) models.

RESULTS

The MI results were concordant with the electro-clinically or surgically defined epileptogenic area in 8/14 patients and showed the same degree of concordance as the results obtained with the GLM-based methods in 12 patients (7 concordant and 5 discordant). In one patient, the information theoretic analysis improved the delineation of the irritative zone compared with the GLM-based methods.

DISCUSSION

Our findings suggest that an information theoretic analysis can provide clinically relevant information about the BOLD signal changes associated with the generation and propagation of interictal epileptic discharges. The concordance between the MI, GLM and FIR maps support the validity of the assumptions adopted in GLM-based analyses of interictal epileptic activity with EEG-fMRI in such a manner that they do not significantly constrain the localization of the epileptogenic zone.

Nonlinear hemodynamic responses in human epilepsy: a multimodal analysis with fNIRS-EEG and fMRI-EEG

Functional magnetic resonance imaging (fMRI) combined with electroencephalography (fMRI-EEG) is a neuroimaging technique based on the blood oxygenation level dependent (BOLD) signal which has been shown to be useful in the study of epilepsy for the localization of the epileptogenic focus. Functional near-infrared spectroscopy (fNIRS) combined with EEG (fNIRS-EEG) is another imaging technique based on the measurement of oxygenated and deoxygenated hemoglobin with complementary clinical potential in epilepsy, for continuous patient monitoring, language lateralization, and focus localization. In this work fMRI-EEG and fNIRS-EEG are used to quantify nonlinear hemodynamic responses in three cases of human refractory focal epilepsy, by using the Volterra kernel expansion up to second order. Prior to analyzing real data, extensive simulations are carried out to show that nonlinearities are estimable. The Volterra methodology is then applied to multimodal data recorded from 3 epileptic patients selected for their frequent spiking activity. Care is taken to account for variability of hemodynamic responses due to other causes than Volterra nonlinearities. Statistically significant nonlinearities are observed for all patients and all modalities. Good concordance between fNIRS and fMRI is found for both the amplitude of the Volterra responses, and, with limitations, in the localization of the epileptic focus and regions of inverted responses (negative BOLD signals). In one patient, Volterra nonlinearities allowed epileptic focus identification with fMRI, while analyses without nonlinearities failed to see it. In simulations when nonlinearities were included, analysis without Volterra nonlinearities performed poorly. These two observations suggest routinely checking for nonlinearities in functional imaging of patients presenting with frequent spikes.

Preictal changes in cerebral haemodynamics: review of findings and insights from intracerebral EEG

The possibility of recording changes in brain signals occurring before epileptic seizures is of considerable interest, both as markers for seizure anticipation and as a window into the mechanisms of seizure generation. Several studies have reported preictal changes on electrophysiological traces. More recently, observations have been made of changes occurring on haemodynamic signals before interictal events or before seizures, often without concurrent changes observed on electrophysiology. We present here a critical review of these findings, in optical imaging, SPECT and fMRI, followed by a discussion based on data from intracerebral EEG.

Tissue hypoxia correlates with intensity of interictal spikes

Interictal spikes (IISs) represent burst firing of a small focal population of hypersynchronous, hyperexcitable cells. Whether cerebral blood flow (CBF) is adequate to meet the metabolic demands of this dramatic increase in membrane excitability is unknown. Positron emission tomography, single photon emission computed tomography, and functional magnetic resonance imaging studies have shown increases in CBF and hypometabolism, thus indicating the likelihood of adequate perfusion. We measured tissue oxygenation and CBF in a rat model of IIS using oxygen electrodes and laser-Doppler flowmetry. A ∼3-second dip in tissue oxygenation was shown, followed by more prolonged tissue hyperoxygenation, in spite of a 25% increase in CBF. Increases in the number of spikes, as well as in their amplitude and spike width further amplified these responses, and a decrease in interspike interval decreased the CBF response. Altering the anesthetic did not influence our results. Taken together, these findings indicate that frequent, high-amplitude IISs may produce significant tissue hypoxia, which has implications for patients with epilepsy and noninvasive techniques of seizure localization.