Augmenting interictal mapping with neurovascular coupling biomarkers by structured factorization of epileptic EEG and fMRI data

Studying the Alzheimer's disease continuum using EEG and fMRI in single-modality and multi-modality settings

Alzheimer's disease (AD) is a biological, clinical continuum that covers the preclinical, prodromal, and clinical phases of the disease. Early diagnosis and identification of the stages of Alzheimer's disease (AD) are crucial in clinical practice. Ideally, biomarkers should reflect the underlying process (pathological or otherwise), be reproducible and non-invasive, and allow repeated measurements over time. However, the currently known biomarkers for AD are not suitable for differentiating the stages and predicting the trajectory of disease progression. Some objective parameters extracted using electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) are widely applied to diagnose the stages of the AD continuum. While electroencephalography (EEG) has a high temporal resolution, fMRI has a high spatial resolution. Combined EEG and fMRI (EEG-fMRI) can overcome single-modality drawbacks and obtain multi-dimensional information simultaneously, and it can help explore the hemodynamic changes associated with the neural oscillations that occur during information processing. This technique has been used in the cognitive field in recent years. This review focuses on the different techniques available for studying the AD continuum, including EEG and fMRI in single-modality and multi-modality settings, and the possible future directions of AD diagnosis using EEG-fMRI.

Mesoscopic mapping of hemodynamic responses and neuronal activity during pharmacologically induced interictal spikes in awake and anesthetized mice

Imaging hemodynamic responses to interictal spikes holds promise for presurgical epilepsy evaluations. Understanding the hemodynamic response function is crucial for accurate interpretation. Prior interictal neurovascular coupling data primarily come from anesthetized animals, impacting reliability. We simultaneously monitored calcium fluctuations in excitatory neurons, hemodynamics, and local field potentials (LFP) during bicuculline-induced interictal events in both isoflurane-anesthetized and awake mice. Isoflurane significantly affected LFP amplitude but had little impact on the amplitude and area of the calcium signal. Anesthesia also dramatically blunted the amplitude and latency of the hemodynamic response, although not its area of spread. Cerebral blood volume change provided the best spatial estimation of excitatory neuronal activity in both states. Targeted silencing of the thalamus in awake mice failed to recapitulate the impact of anesthesia on hemodynamic responses suggesting that isoflurane's interruption of the thalamocortical loop did not contribute either to the dissociation between the LFP and the calcium signal nor to the alterations in interictal neurovascular coupling. The blood volume increase associated with interictal spikes represents a promising mapping signal in both the awake and anesthetized states.

Low-intensity ultrasound stimulation modulates cortical neurovascular coupling in an attention deficit hyperactivity disorder rat model

Attention deficit hyperactivity disorder is accompanied by changes in cranial nerve function and cerebral blood flow (CBF). Low-intensity ultrasound stimulation can modulate brain neural activity in attention deficit hyperactivity disorder. However, to date, the modulatory effects of low-intensity ultrasound stimulation on CBF and neurovascular coupling in attention deficit hyperactivity disorder have not been reported. To address this question, Sprague-Dawley, Wistar-Kyoto, and spontaneously hypertensive (attention deficit hyperactivity disorder (ADHD) rat model) rats were divided into the control and low-intensity ultrasound stimulation (LIUS) groups. Cortical electrical stimulation was used to induce cortical excitability in different types of rats, and a penetrable laser speckle contrast imaging (LSCI) system and electrodes were used to evaluate the electrical stimulation-induced CBF, cortical excitability, and neurovascular coupling in free-moving rats. The CBF, cortical excitability, and neurovascular coupling (NVC) under cortical electrical stimulation in the attention deficit hyperactivity disorder rats were significantly different from those in the Sprague-Dawley and Wistar-Kyoto rats. We also found that low-intensity ultrasound stimulation significantly interfered with the cortical excitability and neurovascular coupling induced by cortical electrical stimulation in rats with attention deficit hyperactivity disorder. Our findings suggest that neurovascular coupling is a potential biomarker for attention deficit hyperactivity disorder. Furthermore, low-intensity ultrasound stimulation can improve abnormal brain function in attention deficit hyperactivity disorder and lay a research foundation for its application in the clinical treatment of attention deficit hyperactivity disorder.

Exploring Neural Mechanisms of Reward Processing Using Coupled Matrix Tensor Factorization: A Simultaneous EEG-fMRI Investigation

BACKGROUND

It is crucial to understand the neural feedback mechanisms and the cognitive decision-making of the brain during the processing of rewards. Here, we report the first attempt for a simultaneous electroencephalography (EEG)-functional magnetic resonance imaging (fMRI) study in a gambling task by utilizing tensor decomposition.

METHODS

First, the single-subject EEG data are represented as a third-order spectrogram tensor to extract frequency features. Next, the EEG and fMRI data are jointly decomposed into a superposition of multiple sources characterized by space-time-frequency profiles using coupled matrix tensor factorization (CMTF). Finally, graph-structured clustering is used to select the most appropriate model according to four quantitative indices.

RESULTS

The results clearly show that not only are the regions of interest (ROIs) found in other literature activated, but also the olfactory cortex and fusiform gyrus which are usually ignored. It is found that regions including the orbitofrontal cortex and insula are activated for both winning and losing stimuli. Meanwhile, regions such as the superior orbital frontal gyrus and anterior cingulate cortex are activated upon winning stimuli, whereas the inferior frontal gyrus, cingulate cortex, and medial superior frontal gyrus are activated upon losing stimuli.

CONCLUSION

This work sheds light on the reward-processing progress, provides a deeper understanding of brain function, and opens a new avenue in the investigation of neurovascular coupling via CMTF.

Deconvolution of the Functional Ultrasound Response in the Mouse Visual Pathway Using Block-Term Decomposition

Functional ultrasound (fUS) indirectly measures brain activity by detecting changes in cerebral blood volume following neural activation. Conventional approaches model such functional neuroimaging data as the convolution between an impulse response, known as the hemodynamic response function (HRF), and a binarized representation of the input signal based on the stimulus onsets, the so-called experimental paradigm (EP). However, the EP may not characterize the whole complexity of the activity-inducing signals that evoke the hemodynamic changes. Furthermore, the HRF is known to vary across brain areas and stimuli. To achieve an adaptable framework that can capture such dynamics of the brain function, we model the multivariate fUS time-series as convolutive mixtures and apply block-term decomposition on a set of lagged fUS autocorrelation matrices, revealing both the region-specific HRFs and the source signals that induce the hemodynamic responses. We test our approach on two mouse-based fUS experiments. In the first experiment, we present a single type of visual stimulus to the mouse, and deconvolve the fUS signal measured within the mouse brain's lateral geniculate nucleus, superior colliculus and visual cortex. We show that the proposed method is able to recover back the time instants at which the stimulus was displayed, and we validate the estimated region-specific HRFs based on prior studies. In the second experiment, we alter the location of the visual stimulus displayed to the mouse, and aim at differentiating the various stimulus locations over time by identifying them as separate sources.

Advanced neuroimaging techniques in epilepsy

PURPOSE OF REVIEW

We review significant advances in epilepsy imaging in recent years.

RECENT FINDINGS

Structural MRI at 7T with optimization of acquisition and postacquisition image processing increases the diagnostic yield but artefactual findings remain a challenge. MRI analysis from multiple sites indicates different atrophy patterns and white matter diffusion abnormalities in temporal lobe and generalized epilepsies, with greater abnormalities close to the presumed seizure source. Structural and functional connectivity relate to seizure spread and generalization; longitudinal studies are needed to clarify the causal relationship of these associations. Diffusion MRI may help predict surgical outcome and network abnormalities extending beyond the epileptogenic zone. Three-dimensional multimodal imaging can increase the precision of epilepsy surgery, improve seizure outcome and reduce complications. Language and memory fMRI are useful predictors of postoperative deficits, and lead to risk minimization. FDG PET is useful for clinical studies and specific ligands probe the pathophysiology of neurochemical fluxes and receptor abnormalities.

SUMMARY

Improved structural MRI increases detection of abnormalities that may underlie epilepsy. Diffusion, structural and functional MRI indicate the widespread associations of epilepsy syndromes. These can assist stratification of surgical outcome and minimize risk. PET has continued utility clinically and for research into the pathophysiology of epilepsies.

Early soft and flexible fusion of electroencephalography and functional magnetic resonance imaging via double coupled matrix tensor factorization for multisubject group analysis

Data fusion refers to the joint analysis of multiple datasets that provide different (e.g., complementary) views of the same task. In general, it can extract more information than separate analyses can. Jointly analyzing electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) measurements has been proved to be highly beneficial to the study of the brain function, mainly because these neuroimaging modalities have complementary spatiotemporal resolution: EEG offers good temporal resolution while fMRI is better in its spatial resolution. The EEG–fMRI fusion methods that have been reported so far ignore the underlying multiway nature of the data in at least one of the modalities and/or rely on very strong assumptions concerning the relation of the respective datasets. For example, in multisubject analysis, it is commonly assumed that the hemodynamic response function is a priori known for all subjects and/or the coupling across corresponding modes is assumed to be exact (hard). In this article, these two limitations are overcome by adopting tensor models for both modalities and by following soft and flexible coupling approaches to implement the multimodal fusion. The obtained results are compared against those of parallel independent component analysis and hard coupling alternatives, with both synthetic and real data (epilepsy and visual oddball paradigm). Our results demonstrate the clear advantage of using soft and flexible coupled tensor decompositions in scenarios that do not conform with the hard coupling assumption.

Three-way parallel group independent component analysis: Fusion of spatial and spatiotemporal magnetic resonance imaging data

Advances in imaging acquisition techniques allow multiple imaging modalities to be collected from the same subject. Each individual modality offers limited yet unique views of the functional, structural, or dynamic temporal features of the brain. Multimodal fusion provides effective ways to leverage these complementary perspectives from multiple modalities. However, the majority of current multimodal fusion approaches involving functional magnetic resonance imaging (fMRI) are limited to 3D feature summaries that do not incorporate its rich temporal information. Thus, we propose a novel three‐way parallel group independent component analysis (pGICA) fusion method that incorporates the first‐level 4D fMRI data (temporal information included) by parallelizing group ICA into parallel ICA via a unified optimization framework. A new variability matrix was defined to capture subject‐wise functional variability and then link it to the mixing matrices of the other two modalities. Simulation results show that the three‐way pGICA provides highly accurate cross‐modality linkage estimation under both weakly and strongly correlated conditions, as well as comparable source estimation under different noise levels. Results using real brain imaging data identified one linked functional–structural–diffusion component associated to differences between schizophrenia and controls. This was replicated in an independent cohort, and the identified components were also correlated with major cognitive domains. Functional network connectivity revealed visual–subcortical and default mode‐cerebellum pairs that discriminate between schizophrenia and controls. Overall, both simulation and real data results support the use of three‐way pGICA to identify multimodal spatiotemporal links and to pursue the study of brain disorders under a single unifying multimodal framework.

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.

Blind Visualization of Task-Related Networks From Visual Oddball Simultaneous EEG-fMRI Data: Spectral or Spatiospectral Model?

Various disease conditions can alter EEG event-related responses and fMRI-BOLD signals. We hypothesized that event-related responses and their clinical alterations are imprinted in the EEG spectral domain as event-related (spatio)spectral patterns (ERSPat). We tested four EEG-fMRI fusion models utilizing EEG power spectra fluctuations (i.e., absolute spectral model - ASM; relative spectral model - RSM; absolute spatiospectral model - ASSM; and relative spatiospectral model - RSSM) for fully automated and blind visualization of task-related neural networks. Two (spatio)spectral patterns (high δ 4 band and low β 1 band) demonstrated significant negative linear relationship (p FWE < 0.05) to the frequent stimulus and three patterns (two low δ 2 and δ 3 bands, and narrow θ 1 band) demonstrated significant positive relationship (p < 0.05) to the target stimulus. These patterns were identified as ERSPats. EEG-fMRI F-map of each δ 4 model showed strong engagement of insula, cuneus, precuneus, basal ganglia, sensory-motor, motor and dorsal part of fronto-parietal control (FPCN) networks with fast HRF peak and noticeable trough. ASM and RSSM emphasized spatial statistics, and the relative power amplified the relationship to the frequent stimulus. For the δ 4 model, we detected a reduced HRF peak amplitude and a magnified HRF trough amplitude in the frontal part of the FPCN, default mode network (DMN) and in the frontal white matter. The frequent-related β 1 patterns visualized less significant and distinct suprathreshold spatial associations. Each θ 1 model showed strong involvement of lateralized left-sided sensory-motor and motor networks with simultaneous basal ganglia co-activations and reduced HRF peak and amplified HRF trough in the frontal part of the FPCN and DMN. The ASM θ 1 model preserved target-related EEG-fMRI associations in the dorsal part of the FPCN. For δ 4, β 1, and θ 1 bands, all models provided high local F-statistics in expected regions. The most robust EEG-fMRI associations were observed for ASM and RSSM.