Epilepsy surgery outcome and functional network alterations in longitudinal MEG: A minimum spanning tree analysis
Introduction
Epilepsy is common in patients with circumscribed brain abnormalities, such as primary brain tumors and mesiotemporal sclerosis. In a substantial number of patients, anti-epileptic drug treatment is ineffective (Berg, 2008, Duffau et al., 2002, Hildebrand et al., 2005, Picot et al., 2008). Many patients with medically intractable non-tumoral lesional epilepsy are referred to epilepsy surgery programs. The aim of these programs is to identify patients in whom it is possible to localize and remove the epileptogenic zone (EZ), i.e. the brain regions that need to be resected to achieve seizure freedom. This strategy succeeds in only 27–67% of these patients, depending on the specific histopathology of the lesion (Tellez-Zenteno et al., 2005). Although the primary aim of surgery in patients with brain tumors is the removal of the tumor, seizure reduction often is an important secondary aim (Chang et al., 2008). For both patient groups, seizure freedom is extremely relevant, as epilepsy is an important limiting factor for quality of life and cognitive functioning (Klein et al., 2003, Markand et al., 2000, Tellez-Zenteno et al., 2007).
The high prevalence of persistent seizures after epilepsy surgery demonstrates that the EZ is insufficiently identified and removed in these patients. The EZ is increasingly seen as an epileptogenic network instead of a localized cortical area, and removal of key regions in this network may increase success rates of epilepsy surgery (Kramer and Cash, 2012, Stam and van Straaten, 2012). Apart from the EZ, the functional organization of the brain network as a whole is disturbed in lesional epilepsy (Kramer and Cash, 2012). Overall, functional connectivity is increased particularly in the delta and theta frequency ranges (0.5–8 Hz) in MEG and EEG recordings, which is a hallmark of (tumor-related) epilepsy (Bettus et al., 2008, Douw et al., 2010a, Horstmann et al., 2010, van Dellen et al., 2012). However, network disturbances are not limited to this pathological increase in slow wave synchrony. The spatial organization or topology of functional neural networks determines to what extent the network facilitates synchronization and the spreading of seizures (Dyhrfjeld-Johnsen et al., 2007, Percha et al., 2005, Varotto et al., 2012). The healthy functional brain network is characterized by a small-world topology, which is thought to be an optimal network topology that combines global integration with local specialization of highly interconnected areas (Bullmore and Sporns, 2009). This small-world topology is lost in lesional epilepsy and primary brain tumor patients (Bartolomei et al., 2006, Bosma et al., 2009, Horstmann et al., 2010), and in particular the loss of connections that integrate different functional regions or clusters correlates with seizure frequency (Douw et al., 2010b, van Dellen et al., 2012). Interestingly, functional connectivity and network characteristics are altered by neurosurgical resections, and such changes are related to postoperative cognitive performance (Douw et al., 2008, van Dellen et al., 2013a). However, so far, surgery-induced alterations in functional network characteristics have not been correlated to postoperative seizure status, which might provide valuable insights for more accurate identification of the target areas for epilepsy surgery.
For this study several choices were made for an approach that builds upon previous work. Firstly, most measures used to characterize networks of lesional epilepsy patients depend on the number of nodes and connection density of the network, hence a normalization strategy is needed. Unfortunately, the normalization approaches that are commonly used introduce their own biases (van Wijk et al., 2010). The use of minimum spanning tree (MST) network analysis provides a principled way to construct unique networks from neurophysiological data, which uses a fixed number of nodes and edges, and is independent of average coupling strength (Boersma et al., 2013). The MST is the subset of strongest connections in the network such that all network nodes are connected, without forming loops. It might represent a critical backbone of information flow in weighted networks (i.e. it contains with high probability all the shortest paths in the network; see Fig. 1 (Wang et al., 2009)).
The use of MST analysis also solves a second methodological concern. Functional connectivity analysis and the small-world model have proven to be a fruitful starting point for studies on brain network topology in lesional epilepsy, but they provide an incomplete representation of the networks, as they do not capture other network features such as the existence of so-called hubs (Stam and van Straaten, 2012). Several studies on electrocorticography (ECoG) recordings show that hub regions, which are highly connected regions with a central role in the network, seem to characterize the EZ (Ortega et al., 2008a, Varotto et al., 2012, Wilke et al., 2010). MST analysis allows for characterization of global integration of information and hubs in one holistic model (Fig. 2), and has been used to characterize ECoG and EEG recordings of epilepsy patients (Lee et al., 2006, Ortega et al., 2008b). These studies showed that networks of left-sided and right-sided mesiotemporal sclerosis patients are dissimilar, and suggested that the EZ is characterized by nodes with a high centrality in the MST.
Finally, previous studies were either based on spatially confined intracranial recordings (ECoG studies), or, in the case of most MEG/EEG studies, were performed in signal-space instead of source-space. The use of a recently developed MEG beamformer approach in combination with proper correction for the effects of volume conduction and field spread makes it possible to study functional connectivity and networks for anatomically defined regions (Hillebrand et al., 2012). This approach makes it possible to study functional network characteristics of lesioned regions, and relate these regional characteristics to epilepsy burden (Douw et al., 2013).
In the present study, we used source-space MEG data to analyze the relation between seizure alterations and functional network alterations in lesional epilepsy patients. Firstly, a cross-sectional analysis was performed to investigate the associations between preoperative seizure frequency and functional connectivity or MST characteristics. We hypothesized a correlation between preoperative seizure frequency and average preoperative functional connectivity, especially in the lower frequency bands (0.5–8 Hz), and between preoperative seizure frequency and functional network structure characterized by pathological hubs (i.e. regions with a high centrality). Secondly, in a longitudinal analysis, the impact of surgery on these measures was compared between patients who were seizure free after surgery (SF), and patients with postoperative seizures (POS). We expected pathologically increased connectivity and pathological hubs to disappear in patients who were seizure free after surgery.
Section snippets
Patients
We included consecutive patients who, during the period 03/2010–08/2011, were referred by the Neurosurgical Center Amsterdam for MEG recordings at the VU University Medical Center. Data were collected as part of the LESION study, which is a prospective longitudinal observational study of patients eligible for lesional epilepsy or tumor surgery. Inclusion criteria were (1) adult (≥ 18 years) patients who (2) had resective surgery for (3A) medically intractable lesional epilepsy, or (3B)
Results
Twenty patients were included in this study (16 male, mean age = 36 years (SD 10)). Group characteristics are shown in Table 1 (see Inline Supplementary Table S1, Inline Supplementary Table S2 for individual patient characteristics). Localization of the lesions is shown in Fig. 4A. Nine patients were included in a previously described study (van Dellen et al., 2013a) (see Inline Supplementary Table S1). Fifteen patients were seizure free at T1, and thirteen patients were seizure free at T2. Some
Discussion
The main finding of this study is that successful epilepsy surgery is reflected in changes in functional network organization that were detected by longitudinal MEG analysis. Firstly, seizure freedom at approximately one year after resective surgery correlated with a shift towards a more integrated network organization (Fig. 5), as measured by an increase in MST leaf fraction. Secondly, MST eccentricity and betweenness centrality decreased significantly in SF patients, also in (but not
Funding
E. van Dellen was supported by the Dutch Epilepsy Foundation (NEF) grant 09-09. L. Douw was supported by the Rubicon grant of the Netherlands Organization for Scientific Research (NWO). Scanning costs were in part funded by the Amsterdam Brain Imaging Platform (ABIP), Amsterdam, The Netherlands.
Author contributions
Study design: EvD LD AH PCW JCB JJH JCR CJS; data collection: EvD LD JCB PCW; provided analysis tools: CJS AH PCW; data analysis: EvD LD AH PCW CJS; written manuscript: EvD LD AH PCW JCB JJH JCR CJS.
Acknowledgments
The authors would like to thank J.M. Meier for useful adjustments to the methods section of this manuscript; G. Engels and W. Cleijne for their assistance with data collection, and the MEG technicians of the Department of Clinical Neurophysiology, N. Sijsma, P.J. Ris and K. Plugge, for the technical assistance.
References (62)
- et al.
Unified segmentation
Neuroimage
(2005) - et al.
Disturbed functional connectivity in brain tumour patients: evaluation by graph analysis of synchronization matrices
Clin. Neurophysiol.
(2006) - et al.
Enhanced EEG functional connectivity in mesial temporal lobe epilepsy
Epilepsy Res.
(2008) - et al.
Treatment-related changes in functional connectivity in brain tumor patients: a magnetoencephalography study
Exp. Neurol.
(2008) - et al.
Cognition is related to resting-state small-world network topology: an magnetoencephalographic study
Neuroscience
(2011) - et al.
Local MEG networks: the missing link between protein expression and epilepsy in glioma patients?
Neuroimage
(2013) - et al.
A quantitative assessment of the sensitivity of whole-head MEG to activity in the adult human cortex
Neuroimage
(2002) - et al.
Frequency-dependent functional connectivity within resting-state networks: an atlas-based MEG beamformer solution
Neuroimage
(2012) - et al.
State dependent properties of epileptic brain networks: comparative graph-theoretical analyses of simultaneously recorded EEG and MEG
Clin. Neurophysiol.
(2010) - et al.
Fine tuning the correlation limit of spatio-temporal signal space separation for magnetoencephalography
J. Neurosci. Methods
(2009)
Age-related changes in modular organization of human brain functional networks
Neuroimage
Validation of head movement correction and spatiotemporal signal space separation in magnetoencephalography
Clin. Neurophysiol.
Complex network analysis of human ECoG data
Neurosci. Lett.
The organization of physiological brain networks
Clin. Neurophysiol.
Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain
Neuroimage
Connectivity in MEG resting-state networks increases after resective surgery for low-grade glioma and correlates with improved cognitive performance
Neuroimage
Local polymorphic delta activity in cortical lesions causes global decreases in functional connectivity
Neuroimage
Epileptogenic networks of type II focal cortical dysplasia: a stereo-EEG study
Neuroimage
Unified segmentation based correction of R1 brain maps for RF transmit field inhomogeneities (UNICORT)
Neuroimage
Modeling the impact of lesions in the human brain
PLoS Comput. Biol.
Realistic spatial sampling for MEG beamformer images
Hum. Brain Mapp.
The natural history of mesial temporal lobe epilepsy
Curr. Opin. Neurol.
Growing trees in child brains: graph theoretical analysis of EEG derived minimum spanning tree in 5 and 7 year old children reflects brain maturation
Brain Connect.
Disturbed functional brain networks and neurocognitive function in low-grade glioma patients: a graph theoretical analysis of resting-state MEG
Nonlinear Biomed. Phys.
Complex brain networks: graph theoretical analysis of structural and functional systems
Nat. Rev. Neurosci.
Seizure characteristics and control following resection in 332 patients with low-grade gliomas
J. Neurosurg.
Resection probability maps for quality assessment of glioma surgery without brain location bias
PLoS One
‘Functional connectivity’ is a sensitive predictor of epilepsy diagnosis after the first seizure
PLoS One
Epilepsy is related to theta band brain connectivity and network topology in brain tumor patients
BMC Neurosci.
Medically intractable epilepsy from insular low-grade gliomas: improvement after an extended lesionectomy
Acta Neurochir. (Wien)
Topological determinants of epileptogenesis in large-scale structural and functional models of the dentate gyrus derived from experimental data
J. Neurophysiol.
Cited by (100)
Network nodes in the brain
2023, Connectome Analysis: Characterization, Methods, and AnalysisAltered brain network topology during speech tracking in developmental dyslexia
2022, NeuroImageCitation Excerpt :Such increase in leaf fraction may be indicative of the presence of highly connected hub nodes and shortened path length (i.e., more efficient information flow) in the network (Stam, 2014; Tewarie et al., 2015). With reduced eccentricity, the right frontal nodes of dyslexic readers exhibited shorter distance (i.e., functionally closer) with other nodes, and hence played a more central or hub-like role in the theta-band network for word perception (Stam et al., 2014; van Dellen et al., 2014). Bilateral frontal cortex has been identified in a multiple-demand system subserving various domain-general functions, including working memory, mental programming and executive control (Diachek, Blank, Siegelman, Affourtit, & Fedorenko, 2020; Duncan, 2010; Duncan & Owen, 2000; Fedorenko, Duncan, & Kanwisher, 2013).
Potential long consequences from internal and external ecology: Loss of gut microbiota antifragility in children from an industrialized population compared with an indigenous rural lifestyle
2023, Journal of Developmental Origins of Health and Disease