fMRI functional connectivity as an indicator of interictal epileptic discharges

Frequency-specific segregation and integration of human cerebral cortex: An intrinsic functional atlas

The cognitive and behavioral functions of the human brain are supported by its frequency multiplexing mechanism. However, there is limited understanding of the dynamics of the functional network topology. This study aims to investigate the frequency-specific topology of the functional human brain using 7T rs-fMRI data. Frequency-specific parcellations were first performed, revealing frequency-dependent dynamics within the frontoparietal control, parietal memory, and visual networks. An intrinsic functional atlas containing 456 parcels was proposed and validated using stereo-EEG. Graph theory analysis suggested that, in addition to the task-positive vs. task-negative organization observed in static networks, there was a cognitive control system additionally from a frequency perspective. The reproducibility and plausibility of the identified hub sets were confirmed through 3T fMRI analysis, and their artificial removal had distinct effects on network topology. These results indicate a more intricate and subtle dynamics of the functional human brain and emphasize the significance of accurate topography.

Altered spread of waves of activities at large scale is influenced by cortical thickness organization in temporal lobe epilepsy: a magnetic resonance imaging-high-density electroencephalography study

Temporal lobe epilepsy is a brain network disorder characterized by alterations at both the structural and the functional levels. It remains unclear how structure and function are related and whether this has any clinical relevance. In the present work, we adopted a novel methodological approach investigating how network structural features influence the large-scale dynamics. The functional network was defined by the spatio-temporal spreading of aperiodic bursts of activations (neuronal avalanches), as observed utilizing high-density electroencephalography in patients with temporal lobe epilepsy. The structural network was modelled as the region-based thickness covariance. Loosely speaking, we quantified the similarity of the cortical thickness of any two brain regions, both across groups and at the individual level, the latter utilizing a novel approach to define the subject-wise structural covariance network. In order to compare the structural and functional networks (at the nodal level), we studied the correlation between the probability that a wave of activity would propagate from a source to a target region and the similarity of the source region thickness as compared with other target brain regions. Building on the recent evidence that large-waves of activities pathologically spread through the epileptogenic network in temporal lobe epilepsy, also during resting state, we hypothesize that the structural cortical organization might influence such altered spatio-temporal dynamics. We observed a stable cluster of structure-function correlation in the bilateral limbic areas across subjects, highlighting group-specific features for left, right and bilateral temporal epilepsy. The involvement of contralateral areas was observed in unilateral temporal lobe epilepsy. We showed that in temporal lobe epilepsy, alterations of structural and functional networks pair in the regions where seizures propagate and are linked to disease severity. In this study, we leveraged on a well-defined model of neurological disease and pushed forward personalization approaches potentially useful in clinical practice. Finally, the methods developed here could be exploited to investigate the relationship between structure-function networks at subject level in other neurological conditions.

Estimation of fMRI responses related to epileptic discharges using Bayesian hierarchical modeling

Simultaneous electroencephalography-functional MRI (EEG-fMRI) is a unique and noninvasive method for epilepsy presurgical evaluation. When selecting voxels by null-hypothesis tests, the conventional analysis may overestimate fMRI response amplitudes related to interictal epileptic discharges (IEDs), especially when IEDs are rare. We aimed to estimate fMRI response amplitudes represented by blood oxygen level dependent (BOLD) percentage changes related to IEDs using a hierarchical model. It involves the local and distributed hemodynamic response homogeneity to regularize estimations. Bayesian inference was applied to fit the model. Eighty-two epilepsy patients who underwent EEG-fMRI and subsequent surgery were included in this study. A conventional voxel-wise general linear model was compared to the hierarchical model on estimated fMRI response amplitudes and on the concordance between the highest response cluster and the surgical cavity. The voxel-wise model overestimated fMRI responses compared to the hierarchical model, evidenced by a practically and statistically significant difference between the estimated BOLD percentage changes. Only the hierarchical model differentiated brief and long-lasting IEDs with significantly different BOLD percentage changes. Overall, the hierarchical model outperformed the voxel-wise model on presurgical evaluation, measured by higher prediction performance. When compared with a previous study, the hierarchical model showed higher performance metric values, but the same or lower sensitivity. Our results demonstrated the capability of the hierarchical model of providing more physiologically reasonable and more accurate estimations of fMRI response amplitudes induced by IEDs. To enhance the sensitivity of EEG-fMRI for presurgical evaluation, it may be necessary to incorporate more appropriate spatial priors and bespoke decision strategies.

Integration of white matter architecture to stereo-EEG better describes epileptic spike propagation

OBJECTIVE

Stereo-electroencephalography (SEEG)-derived epilepsy networks are used to better understand a patient's epilepsy; however, a unimodal approach provides an incomplete picture. We combine tractography and SEEG to determine the relationship between spike propagation and the white matter architecture and to improve our understanding of spike propagation mechanisms.

METHODS

Probablistic tractography from diffusion imaging (dMRI) of matched subjects from the Human Connectome Project (HCP) was combined with patient-specific SEEG-derived spike propagation networks. Two regions-of-interest (ROIs) with a significant spike propagation relationship constituted a Propagation Pair.

RESULTS

In 56 of 59 patients, Propagation Pairs were more often tract-connected as compared to all ROI pairs (p < 0.01; d = -1.91). The degree of spike propagation between tract-connected ROIs was greater (39 ± 21%) compared to tract-unconnected ROIs (31 ± 18%; p < 0.0001). Within the same network, ROIs receiving propagation earlier were more often tract-connected to the source (59.7%) as compared to late receivers (25.4%; p < 0.0001).

CONCLUSIONS

Brain regions involved in spike propagation are more likely to be connected by white matter tracts. Between nodes, presence of tracts suggests a direct course of propagation, whereas the absence of tracts suggests an indirect course of propagation.

SIGNIFICANCE

We demonstrate a logical and consistent relationship between spike propagation and the white matter architecture.

Resting-State fMRI Can Detect Alterations in Seizure Onset and Spread Regions in Patients with Non-Lesional Epilepsy: A Pilot Study

Introduction

Epilepsy is defined as non-lesional (NLE) when a lesion cannot be localized via standard neuroimaging. NLE is known to have a poor response to surgery. Stereotactic electroencephalography (sEEG) can detect functional connectivity (FC) between zones of seizure onset (OZ) and early (ESZ) and late (LSZ) spread. We examined whether resting-state fMRI (rsfMRI) can detect FC alterations in NLE to see whether noninvasive imaging techniques can localize areas of seizure propagation to potentially target for intervention.

Methods

This is a retrospective study of 8 patients with refractory NLE who underwent sEEG electrode implantation and 10 controls. The OZ, ESZ, and LSZ were identified by generating regions around sEEG contacts that recorded seizure activity. Amplitude synchronization analysis was used to detect the correlation of the OZ to the ESZ. This was also done using the OZ and ESZ of each NLE patient for each control. Patients with NLE were compared to controls individually using Wilcoxon tests and as a group using Mann-Whitney tests. Amplitude of low-frequency fluctuations (ALFF), fractional ALFF (fALFF), regional homogeneity (ReHo), degree of centrality (DoC), and voxel-mirrored homotopic connectivity (VMHC) were calculated as the difference between NLE and controls and compared between the OZ and ESZ and to zero. A general linear model was used with age as a covariate with Bonferroni correction for multiple comparisons.

Results

Five out of 8 patients with NLE showed decreased correlations from the OZ to the ESZ. Group analysis showed patients with NLE had lower connectivity with the ESZ. Patients with NLE showed higher fALFF and ReHo in the OZ but not the ESZ, and higher DoC in the OZ and ESZ. Our results indicate that patients with NLE show high levels of activity but dysfunctional connections in seizure-related areas.

Discussion

rsfMRI analysis showed decreased connectivity directly between seizure-related areas, while FC metric analysis revealed increases in local and global connectivity in seizure-related areas. FC analysis of rsfMRI can detect functional disruption that may expose the pathophysiology underlying NLE.

Parcellating the human brain using resting-state dynamic functional connectivity

Brain cartography has expanded substantially over the past decade. In this regard, resting-state functional connectivity (FC) plays a key role in identifying the locations of putative functional borders. However, scant attention has been paid to the dynamic nature of functional interactions in the human brain. Indeed, FC is typically assumed to be stationary across time, which may obscure potential or subtle functional boundaries, particularly in regions with high flexibility and adaptability. In this study, we developed a dynamic FC (dFC)-based parcellation framework, established a new functional human brain atlas termed D-BFA (DFC-based Brain Functional Atlas), and verified its neurophysiological plausibility by stereo-EEG data. As the first dFC-based whole-brain atlas, the proposed D-BFA delineates finer functional boundaries that cannot be captured by static FC, and is further supported by good correspondence with cytoarchitectonic areas and task activation maps. Moreover, the D-BFA reveals the spatial distribution of dynamic variability across the brain and generates more homogenous parcels compared with most alternative parcellations. Our results demonstrate the superiority and practicability of dFC in brain parcellation, providing a new template to exploit brain topographic organization from a dynamic perspective. The D-BFA will be publicly available for download at https://github.com/sliderplm/D-BFA-618.

Clinical applications of magnetic resonance imaging based functional and structural connectivity

Advances in computational neuroimaging techniques have expanded the armamentarium of imaging tools available for clinical applications in clinical neuroscience. Non-invasive, in vivo brain MRI structural and functional network mapping has been used to identify therapeutic targets, define eloquent brain regions to preserve, and gain insight into pathological processes and treatments as well as prognostic biomarkers. These tools have the real potential to inform patient-specific treatment strategies. Nevertheless, a realistic appraisal of clinical utility is needed that balances the growing excitement and interest in the field with important limitations associated with these techniques. Quality of the raw data, minutiae of the processing methodology, and the statistical models applied can all impact on the results and their interpretation. A lack of standardization in data acquisition and processing has also resulted in issues with reproducibility. This limitation has had a direct impact on the reliability of these tools and ultimately, confidence in their clinical use. Advances in MRI technology and computational power as well as automation and standardization of processing methods, including machine learning approaches, may help address some of these issues and make these tools more reliable in clinical use. In this review, we will highlight the current clinical uses of MRI connectomics in the diagnosis and treatment of neurological disorders; balancing emerging applications and technologies with limitations of connectivity analytic approaches to present an encompassing and appropriate perspective.

Localization of Epileptic Foci Based on Simultaneous EEG-fMRI Data

Combining functional magnetic resonance imaging (fMRI) and electroencephalography (EEG) enables a non-invasive investigation of the human brain function and evaluation of the correlation of these two important modalities of brain activity. This paper explores recent reports on using advanced simultaneous EEG–fMRI methods proposed to map the regions and networks involved in focal epileptic seizure generation. One of the applications of EEG and fMRI combination as a valuable clinical approach is the pre-surgical evaluation of patients with epilepsy to map and localize the precise brain regions associated with epileptiform activity. In the process of conventional analysis using EEG–fMRI data, the interictal epileptiform discharges (IEDs) are visually extracted from the EEG data to be convolved as binary events with a predefined hemodynamic response function (HRF) to provide a model of epileptiform BOLD activity and use as a regressor for general linear model (GLM) analysis of the fMRI data. This review examines the methodologies involved in performing such studies, including techniques used for the recording of EEG inside the scanner, artifact removal, and statistical analysis of the fMRI signal. It then discusses the results reported for patients with primary generalized epilepsy and patients with different types of focal epileptic disorders. An important matter that these results have brought to light is that the brain regions affected by interictal epileptic discharges might not be limited to the ones where they have been generated. The developed methods can help reveal the regions involved in or affected by a seizure onset zone (SOZ). As confirmed by the reviewed literature, EEG–fMRI provides information that comes particularly useful when evaluating patients with refractory epilepsy for surgery.

Enhanced Regional Functional Connectivity Indicates Seizure Onset Zone

This project aims to explore if stronger functional connectivity (FC) exists in the maximal BOLD response of EEG/fMRI analysis when it is concordant with seizure-onset-zone (SOZ). Twenty-six patients with drug-resistant focal epilepsy who had an EEG/fMRI and later underwent stereo-EEG implantation were included. Different types of IEDs were labeled in scalp EEG and IED-related maximal BOLD responses were evaluated separately, each constituting one study. After evaluating concordance between maximal BOLD and SOZ, twenty-seven studies were placed in the concordant group and eight in the discordant group. We evaluated the local connectivity and ipsilaterally distant connectivity difference between the maximal BOLD and the contralateral homotopic region. Significantly stronger local FC was found for the maximal BOLD in the concordant group (p < 0.05, Bonferroni corrected). 52% of the studies in the concordant group and 13% in the discordant group had a significant difference compared to healthy subjects (p < 0.05, uncorrected). The finding suggests that, when concordant with the SOZ, the maximal BOLD is more likely to have stronger local FC compared to its contralateral counterpart. This asymmetry in functional connectivity may help to noninvasively improve the specificity of EEG/fMRI analysis.