|
||||||||||||||||||||
|
|
Use of resting-state fMRI in planning epilepsy surgery
Correspondence Address: Source of Support: None, Conflict of Interest: None DOI: 10.4103/neuroindia.NI_823_16
Epileptic seizures result from abnormal neuronal excitability and synchronization, affecting 0.5–1% of the population worldwide. Although anti-seizure drugs are often effective, a significant number of patients with epilepsy continue to experience refractory seizures and are candidates for surgical resection. Whereas standard presurgical evaluation has relied on intracranial electroencephalography (icEEG) and direct cortical stimulation to identify epileptogenic tissue and areas of cortex for which resection would produce clinical deficits, the invasive nature and limited spatial extent of icEEG has led to the investigation of less invasive imaging modalities as adjunctive tools in the presurgical workup. In the past few decades, functional connectivity MRI has emerged as a promising approach for presurgical mapping, leading to a surge in the number of proposed methods and biomarkers for identifying epileptogenic tissue. This review focuses on recent advances in the use of functional connectivity MRI toward its application for presurgical planning, including epilepsy localization and eloquent cortex mapping. Keywords: Eloquent cortex, epilepsy, epileptogenic zone, functional connectivity magnetic resonance imaging, presurgical planning
Epilepsy is one of the most prevalent, serious, chronic neurological disorders, affecting nearly 0.5–1% of the population worldwide.[1] For the 20–30% of patients suffering from epilepsy whose seizures are resistant to medications,[2] surgical resection of epileptogenic tissue may be an effective method of treatment. However, determining candidacy for epilepsy surgery requires sophisticated testing to precisely delineate the epileptogenic tissue and its anatomic relationship with the eloquent cortex. This is needed to achieve postoperative seizure freedom while minimizing the likelihood of postoperative focal neurologic deficits. Current standard noninvasive presurgical evaluation employs multiple imaging and electrophysiological modalities, almost always including structural magnetic resonance imaging (MRI) and surface electroencephalography (EEG). However, such noninvasive assessment is often limited by normal structural MRI and nonfocal surface EEG. As such, it is nonconfirmatory for the epileptogenic zone (EZ) in a considerable proportion of patients with pharmacoresistant epilepsy.[3],[4] Intracranial EEG (icEEG) remains the necessary test in such situations for delineation of the EZ, albeit with limitations, including the extent of coverage possible. As the invasive nature of intracranial EEG carries risks, noninvasive methods are preferable. The goal of the evaluation is the identification of the EZ, which is defined as the area necessary and sufficient for spontaneous seizure generation. Recent studies, however, have broadened this concept, suggesting that ictal activity may begin in a spatially localized EZ which spreads to, and is modulated by, a larger network of brain structures.[5] In contrast to the spatially localized EZ, an epileptogenic network is dynamic and may change with a longer epilepsy duration.[6] Strategic disruption of the epileptogenic network may allow for seizure elimination while decreasing the risk of recurrence. This approach is now being conceptualized as results from investigations into the epileptic network are being reported. If effective, the epileptic network characterization would alter the approach to epilepsy surgery. There has been an increasing interest in functional MRI (fMRI) investigations of the human connectome, and such investigations have considerable value in understanding epilepsy networks with a potential to advance surgical treatment. This use of fMRI differs from the more established approach at many epilepsy surgery centers of presurgically mapping the eloquent cortex, and sometimes includes a simultaneous EEG to localize regions generating interictal or ictal activity.[7],[8],[9],[10] For mapping the eloquent function or epileptic activity, the approach seeks to identify areas of cortex with fMRI blood-oxygen-level dependent (BOLD) signal activations associated with particular tasks. such as in eloquent cortex mapping, or in response to EEG events, as in the identification of interictal or ictal-onset zones. This use of fMRI has at times supplanted the traditional approaches for identifying the eloquent cortex, such as the Wada test for language lateralization, or cortical mapping by electrical stimulation during electrocorticography. As Wada testing may be used to lateralize but not localize language, and because of the invasive nature of direct cortical stimulation, task-based fMRI can be superior in both anatomic detail and scope, as well as safer by virtue of its noninvasiveness. Questions relating to the optimal task choice, as well as the extrapolation of task designs from healthy to patient populations, however, provide a challenge. As task-based fMRI depends on the patient's ability to perform the needed task, it may be suboptimal in children and in patients with neurological deficits, altered mental status, or under sedation. In addition, the need to perform each task separately requires long acquisition times if multiple functional sites are to be investigated. Functional connectivity MRI (fcMRI) differs from task-based fMRI by recording a resting-state. Rather than performing an explicit task, the patient is asked to lie still and not to think about anything in particular. It is during this undirected mental state that the default mode network is active.[11] Instead of associating the fMRI BOLD signal to a behavior or EEG finding, connectivity is calculated from spontaneous BOLD signal fluctuations and identifying the relationship in the fluctuations across cerebral regions. Estimated correlations between low-frequency fluctuations (<0.1 Hz) have been found to reproducibly reflect the underlying intrinsic neural architecture,[12] and are presumed to be generated by underlying fluctuations in neural activity. These correlations may be complementary in identifying seizure-generating networks as well as defining eloquent cortex or networks. Graph theory provides an additional network-level perspective by quantifying how each given voxel or region-of-interest is connected to the remainder of the network. Network “hubs,” or regions, which are central to the functioning of the network, may be quantified through measures such as degree or betweenness centrality, whereas measures of local and/or long-range connectivity may be quantified through measures such as the clustering coefficient or path length, respectively.[13] With respect to both localization of the EZ and eloquent cortex mapping, there is a need for less invasive and more reliable diagnostic tools. In this article, we review the converging evidence from fMRI connectomic analyses toward their application to presurgical planning. In particular, we focus on the clinical impact of this evidence for localization of epileptogenic tissue and eloquent cortex mapping. A summary of the evidence underlying the following discussion is provided in [Table 1] and [Table 2].
Electroencephalography–functional magnetic resonance imaging Interictal epileptiform discharges (IEDs) reflect abnormal interictal neural activity and arise from the so-called irritative zone (IZ). One approach for mapping the IZ utilizes simultaneous EEG–fMRI recordings, which combine the advantages of the specificity of EEG for epileptic abnormality and the spatial resolution of fMRI to map changes in the cerebral hemodynamic flow corresponding to IEDs. IEDs are identified in EEG data, convolved with a hemodynamic response function (HRF), and treated as regressors of interest in a General Linear Model (GLM). Areas that are activated during IEDs are then identified as the IZ [Figure 1]. While EEG has classically been used for IZ identification, the spatial resolution and tomography of fMRI is significantly higher than EEG, permitting a depiction of the IZ's hemodynamic map. Studies report that greater correspondence between the IZ delineated by EEG–fMRI and resected tissue leads to a greater probability of postoperative seizure freedom. One group found that 58% of the patients with resection areas concordant with the IZ and 19% of the patients with discordant areas attain good postsurgical seizure outcomes (Engel 1-2).[14] Similarly, 80% of the patients with concordant IZ and seizure-onset zone (SOZ) attained >50% reduction in seizure frequency one year after surgery, compared to only 17% of the patients with discordant IZ and SOZ.[15] Such studies suggest that simultaneous EEG–fMRI may be a useful adjunct in the preoperative evaluation of epileptogenic cortex.
However, the EEG–fMRI approach to localizing epileptogenic tissue suffers from a few limitations:
Functional connectivity magnetic resonance imaging fcMRI provides a promising alternative approach for identifying the IZ without necessitating assumptions of the GLM. In particular, IED-associated changes in the graph theory measures of network properties, rather than in BOLD amplitude, may facilitate IZ delineation. As graph measures are calculated based on the number of connections for each voxel, this approach obviates the need for GLM assumptions. One approach, for example, spatially maps the IZ using IED-associated changes in local degree centrality, a graph measure of local network connectivity. The clinically presumed EZ had increases in local degree centrality during IEDs compared to baseline. Furthermore, whereas IZ identified using this method was concordant with the presumed EZ in all six patients in the sample, GLM-based methods were concordant in only three of the six patients.[26] Although large sample validation is needed, the potential promise of non-task-based approaches to IZ mapping may be because of the ability to detect BOLD changes not identifiable under GLM-based assumptions, such as when IEDs do not appear on the scalp EEG or when IEDs cluster in a short period, violating assumptions of linear neurovascular coupling.[26] Various fcMRI graph theory studies have supported the potential for resting-state graph theory measures to lateralize the clinically presumed epileptogenic focus.[27],[28],[29],[30] However, fewer studies have addressed the sensitivity and specificity of these measures for focus lateralization.[31],[32],[33] Functional connectivity studies have also shown that abnormalities in the strength of functional connectivity in the absence of IED information may serve as a marker for the EZ. The direction in which functional connectivity changes near the suspected EZ has not been definitely established. While several fcMRI studies have also found evidence that local functional connectivity is increased for epileptogenic tissue,[34],[35],[36] other research has shown decreases.[33],[37] Indeed, icEEG studies have shown increases in functional connectivity in the EZ.[38],[39] Discrepant conclusions may potentially be due to patient heterogeneity, small sample sizes, varying definitions of local neighborhoods, and different local connectivity metrics and modalities. Differences between the SOZ,[35],[37] IZ,[34] and a priori defined regions of interest (ROIs)[33],[36] may also contribute. Recent evidence furthermore indicates that connectivity is dynamic both within and across scans,[40],[41],[42] such that traditional static measures of connectivity provide only an average measure of functional connectivity. Local measures of connectivity exhibit significant temporal nonstationarity,[42] potentially impacting whether increases or decreases in functional connectivity (FC) are identified. Discrepant results from fcMRI compared to icEEG may also result from neurovascular decoupling in epilepsy and differences in timescales.[33],[39] Other changes in brain connectivity in epilepsy Other changes have also been identified in brain connectivity that contribute to understanding the pathophysiology of epilepsy and may provide future targets for intervention by surgery or devices. Connectivity to spatially remote regions appears altered for the presumed EZ and may indicate diffuse epileptic networks. However, the direction of this alteration is not consistent. Some studies have investigated FC between the IZ and remote regions and found evidence that FC is decreased,[34] whereas other studies on the SOZ have found that connectivity to remote regions is increased. Among patients with left TLE, the superior left temporal pole has been found to have increased connectivity with neighboring regions such as within the left temporal and frontal cortex; however, decreased connectivity to more remote regions such as bilateral dorsal posterior cingulate cortex, contralateral inferior temporal gyrus, and contralateral thalamus occurs.[36] Another study of patients with TLE found that the hippocampus has increased connectivity to key areas of the limbic system and frontal lobes, and decreased connectivity to various sensorimotor areas.[30] These variations may reflect the fact that increases in remote functional connectivity may be less specific to epileptogenic tissue than increases in local functional connectivity.[35] In general, data from EEG–fMRI suggests that more diffuse epileptogenic networks are associated with poor postoperative outcomes.[14],[15],[43] Measures of the diffusivity of epileptogenic networks, such as laterality indices of functional connectivity computed from the presumed epileptogenic focus,[43] may be useful for the prediction of surgical success. Some studies have investigated whether combining local and remote connectivity may be useful as a marker for localization,[37],[44] based on the hypothesis that contralateral changes in connectivity may result as a compensatory mechanism for abnormalities ipsilateral to the seizure focus. As various aspects of functional connectivity appear to be affected by the progression of epilepsy, including homogenization of connectivity strengths [45] and the diminishment of functional connectivity in ipsilateral networks,[46] markers which account for both local and remote changes may be useful in providing markers less dependent on epilepsy burden. Various approaches to considering relative amounts of ipsilateral and contralateral connectivity have been considered, including computing voxel-wise difference maps between intra- and interhemispheric connectivity,[44] ROI-based approaches which compare the local connectivity in an ROI containing the SOZ to that contralateral to the clinically presumed epileptogenic focus,[37] and ROI-based approaches which compare ipsilateral versus contralateral connectivity among predefined ROIs thought to be typically involved in epileptogenic networks.[33] These studies consistently suggest that functional connectivity is decreased ipsilateral to the seizure focus relative to the contralateral hemisphere. Less chaotic BOLD dynamics within the functional deficit zone, defined as the area of cortex that is interictally dysfunctional, have also been identified in association with increased local functional connectivity, which may reflect increased density of inhibitory or excitatory inputs.[36]
Task-based fMRI has traditionally been used for the preoperative delineation of eloquent cortex. During a typical task-based fMRI session, a stimulus is administered to the patient. Finger-tapping paradigms are commonly used to localize sensorimotor cortex, whereas verb-generation or sentence-completion paradigms are commonly used to localize language cortex. Areas with BOLD signal change associated with these tasks are identified as the corresponding cortical representation area of the cortex. Resting-state fMRI (rs-fMRI) provides a noninvasive approach and allows for the simultaneous mapping of multiple functions. Several approaches have been proposed for applying rs-fMRI for preoperative localization of eloquent cortex. This would obviate the task paradigm and would be easier for the patient and potentially be more reliable. The resting-state approach has been directed toward sensorimotor and language networks in particular. A traditional approach to mapping functional networks involves seed-based cross-correlation analysis [Figure 2]a. In this approach, a seed voxel or cluster is determined based on brain anatomy or previous functional studies, and the functional network is then extracted by computing the univariate correlations between the BOLD signal of each voxel with the seed region.[47] Seed-based cross-correlation analysis has been found to yield sensorimotor maps consistent with electrical stimulation mapping (ESM), and to perform at least as well as task-based fMRI for consistency with ESM.[48],[49] An alternative data-driven technique for identifying functional networks in rs-fMRI data, which does not require predefined seeds, is spatial independent component analysis (ICA). In spatial ICA, the observed rs-fMRI data is mathematically decomposed into spatially nonoverlapping and temporally coherent regions, which correspond to various resting-state networks [Figure 2]b.[50],[51] Several studies have consistently shown that ICA can identify sensorimotor networks in healthy individuals.[52],[53],[54],[55],[56] Studies in patients with epilepsy or brain tumors have also found that ICA of rs-fMRI data produces sensorimotor maps similar to those obtained from task-based GLM mapping.[57] Comparisons of maps elicited through rs-fMRI also appear to be consistent with sensorimotor maps elicited through ESM in patients with epilepsy.[49],[58] In general, seed-based approaches have been found to be potentially problematic and generally inferior to ICA because of the sensitivity of seed placement.[59] ICA, which does not require predefined seeds, is thought to offer an advantage in this case. However, ICA has its own drawbacks, including the need to specify the number of independent components as well as for user intervention to identify which components reflect neural systems versus artifacts. Template matching methods have been proposed in which group-level ICA is used to identify appropriate networks on the individual patient level.[60] Other approaches avoid reliance on templates by combining patient-level data into a single ICA followed by back-reconstruction of individual patient ICs,[61] or other approaches such as “partner-matching,” which identifies reproducible components by clustering components within or across patients.[61] In addition, findings from ICA can also be influenced by the number of ICs if the number of ICs is too small.[59] A critical number of ICs below which ICA results are less robust was estimated empirically as nine,[59] although many studies choose a relatively higher order ICA (>60 ICs),[51],[57],[62] which reduces the loss of important information within a reasonable computation time. In the case of the sensorimotor cortex, seed-based mapping may efficiently map sensorimotor cortex as indicated by comparisons with ESM, and appears to yield similar maps to those obtained through ICA.
Identification of the language area is more difficult because of the greater locational variability across individuals. Furthermore, patients with epilepsy may have additional lesion-induced architectural rearrangements, further increasing the variability of the language network's spatial distribution.[63] Therefore, employing seed-based cross-correlation analysis is more challenging for the language network, and particularly for patients with abnormal brain anatomy. ICA or neural network methods may hold more promise than seed-based methods for mapping language cortex. One recent study on healthy controls has shown that blind-source separation of rs-fMRI data using ICA, followed by semi-automatic classification of the language network, yields good overlap with the language network obtained from task fMRI.[64] Recent studies to evaluate the feasibility of rs-fMRI for identifying language networks among patients with epilepsy have also yielded favorable results. An initial study among six patients with temporal and frontal lobe epilepsies showed fair sensitivity and specificity for rs-fMRI language mapping compared to ESM when back-propagation in artificial neural networks was used.[58] While promising, the reliance of this approach on an adequate training set may lead to difficulties, particularly when dealing with patients whose brains exhibit developmental or plastic effects. Future developments in data-sharing collaborations, such as through the Human Connectome Project [65] and the Brain Genomics Superstruct Project,[66] may be helpful in this regard. ICA provides an avenue for extracting functional networks while avoiding reliance on a training set. A fair amount of overlap has been found between language networks extracted using ICA from rs-fMRI data compared to the language network obtained from task fMRI, in a mixed group of 12 patients with brain tumors and three patients with epilepsy. Furthermore, stability of language IC maps was found to be relatively stable above 50 components.[67]
The attainment of seizure freedom while preserving functions for patients with medication-resistant epilepsy depends on the precision with which areas with the greatest likelihood to compose the epileptogenic zone can be identified. While icEEG remains the gold standard for identifying the epileptogenic zone, the spatial coverage and noninvasive nature of rs-fMRI suggests its potential as an adjunct in presurgical EZ localization and eloquent cortex mapping. Rapidly evolving research involving multimodal imaging and dynamic functional connectivity is anticipated to further the ability of rs-fMRI to identify markers for epileptogenicity. These developments may allow for better understanding of how epileptogenesis results from abnormal connectivity, as well as guide more effective surgical interventions. Acknowledgements Funding/support for this research was provided by: (1) NIH-NINDS K23 Grant NS044936 (JMS); (2) The Leff Family Foundation (JMS); (3) P01 NS02808 (JE); (4) R01 NS33310 (JE); (5) U01 NS42372 (JE); (6) P20 NS80181 (JE). Financial support and sponsorship Nil. Conflicts of interest There are no conflicts of interest
[Figure 1], [Figure 2]
[Table 1], [Table 2]
|
|
|||||