Skip to main content

Main menu

  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE

User menu

  • Log out
  • Log in
  • My Cart

Search

  • Advanced search
Journal of Neuroscience
  • Log out
  • Log in
  • My Cart
Journal of Neuroscience

Advanced Search

Submit a Manuscript
  • HOME
  • CONTENT
    • Early Release
    • Featured
    • Current Issue
    • Issue Archive
    • Collections
    • Podcast
  • ALERTS
  • FOR AUTHORS
    • Information for Authors
    • Fees
    • Journal Clubs
    • eLetters
    • Submit
  • EDITORIAL BOARD
  • ABOUT
    • Overview
    • Advertise
    • For the Media
    • Rights and Permissions
    • Privacy Policy
    • Feedback
  • SUBSCRIBE
PreviousNext
Research Articles, Systems/Circuits

Dynamical Mechanisms of Interictal Resting-State Functional Connectivity in Epilepsy

Julie Courtiol, Maxime Guye, Fabrice Bartolomei, Spase Petkoski and Viktor K. Jirsa
Journal of Neuroscience 15 July 2020, 40 (29) 5572-5588; DOI: https://doi.org/10.1523/JNEUROSCI.0905-19.2020
Julie Courtiol
1Aix-Marseille Univ, Inserm, INS, Institut de Neurosciences des Systèmes, 13005 Marseille, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Maxime Guye
2Aix-Marseille Univ, CNRS, Centre de Résonance Magnétique et Biologique et Médicale (CRMBM), 13005 Marseille, France
3Assistance Publique-Hôpitaux de Marseille, Hôpital de La Timone, CEMEREM, Pôle d'Imagerie Médicale, CHU, 13005 Marseille, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Fabrice Bartolomei
1Aix-Marseille Univ, Inserm, INS, Institut de Neurosciences des Systèmes, 13005 Marseille, France
4Assistance Publique-Hôpitaux de Marseille, Hôpital de La Timone, Service de Neurophysiologie Clinique, CHU, 13005 Marseille, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Spase Petkoski
1Aix-Marseille Univ, Inserm, INS, Institut de Neurosciences des Systèmes, 13005 Marseille, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Spase Petkoski
Viktor K. Jirsa
1Aix-Marseille Univ, Inserm, INS, Institut de Neurosciences des Systèmes, 13005 Marseille, France
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

Drug-resistant focal epilepsy is a large-scale brain networks disorder characterized by altered spatiotemporal patterns of functional connectivity (FC), even during interictal resting state (RS). Although RS-FC-based metrics can detect these changes, results from RS functional magnetic resonance imaging (RS-fMRI) studies are unclear and difficult to interpret, and the underlying dynamical mechanisms are still largely unknown. To better capture the RS dynamics, we phenomenologically extended the neural mass model of partial seizures, the Epileptor, by including two neuron subpopulations of epileptogenic and nonepileptogenic type, making it capable of producing physiological oscillations in addition to the epileptiform activity. Using the neuroinformatics platform The Virtual Brain, we reconstructed 14 epileptic and 5 healthy human (of either sex) brain network models (BNMs), based on individual anatomical connectivity and clinically defined epileptogenic heatmaps. Through systematic parameter exploration and fitting to neuroimaging data, we demonstrated that epileptic brains during interictal RS are associated with lower global excitability induced by a shift in the working point of the model, indicating that epileptic brains operate closer to a stable equilibrium point than healthy brains. Moreover, we showed that functional networks are unaffected by interictal spikes, corroborating previous experimental findings; additionally, we observed higher excitability in epileptogenic regions, in agreement with the data. We shed light on new dynamical mechanisms responsible for altered RS-FC in epilepsy, involving the following two key factors: (1) a shift of excitability of the whole brain leading to increased stability; and (2) a locally increased excitability in the epileptogenic regions supporting the mixture of hyperconnectivity and hypoconnectivity in these areas.

SIGNIFICANCE STATEMENT Advances in functional neuroimaging provide compelling evidence for epilepsy-related brain network alterations, even during the interictal resting state (RS). However, the dynamical mechanisms underlying these changes are still elusive. To identify local and network processes behind the RS-functional connectivity (FC) spatiotemporal patterns, we systematically manipulated the local excitability and the global coupling in the virtual human epileptic patient brain network models (BNMs), complemented by the analysis of the impact of interictal spikes and fitting to the neuroimaging data. Our results suggest that a global shift of the dynamic working point of the brain model, coupled with locally hyperexcitable node dynamics of the epileptogenic networks, provides a mechanistic explanation of the epileptic processes during the interictal RS period. These, in turn, are associated with the changes in FC.

  • brain dynamics
  • brain network model
  • epilepsy
  • fMRI
  • functional connectivity
  • resting state

Introduction

Drug-resistant focal epilepsy is a large-scale brain networks disorder (Bartolomei et al., 2005, 2008, 2017). These networks are typically studied using resting state (RS)-functional connectivity (FC)—the statistical interdependencies in the spontaneous activity between brain regions (Biswal et al., 1995)—which can be performed during interictal periods (Guye et al., 2008; Tracy and Doucet, 2015). The majority of RS-FC studies use fMRI without a precise hypothesis of the organization of the epileptogenic network (Waites et al., 2006; Voets et al., 2014). Various contradictory findings were reported, such as globally increased/decreased connectivity (Pereira et al., 2010; Pittau et al., 2012; Wirsich et al., 2016), or concomitantly increased and decreased connectivity in neighboring networks (Bettus et al., 2009, 2010; Liao et al., 2010; Su et al., 2015). In contrast, most RS-FC studies using magnetoencephalographic (MEG) and electroencephalographic (EEG) recordings indicate predominantly increased connectivity between the regions found to be less connected in fMRI (Bettus et al., 2008; Bartolomei et al., 2013a,b; Schevon et al., 2007; Lagarde et al., 2018; Li Hegner et al., 2018). These discrepancies could be attributed to the different functional aspects captured by each modality or the different RS-FC estimates applied (Ridley et al., 2017). Vigilance fluctuations during the RS-fMRI scans could also be associated with (dynamic) changes in connectivity (Haimovici et al., 2017; Laumann et al., 2017). The effect of interictal epileptiform discharges (IEDs) on the abnormalities of brain networks, such as interictal spikes, is also largely neglected in these studies and might constitute another confounding factor (Bettus et al., 2008; Coito et al., 2016).

In contrast to conventional RS analyses, studies based on brain network models (BNMs) have provided important insights into the mechanisms underlying RS dynamics in health (Ghosh et al., 2008; Deco and Corbetta, 2011) and disease (Cabral et al., 2013; Demirtaş et al., 2017). Spontaneous activity is shaped by the structural connectivity (SC) and the dynamical working point, which refers to the parameters of the model maximally fitting the empirical data. An important common feature of healthy BNMs is the emergence of structured RS fluctuations when the system operates near criticality, namely at the border of a qualitative change in behavior (Deco and Jirsa, 2012; Deco et al., 2013), and shows maximal metastability (i.e., variability of synchronization; Deco et al., 2017a). In epilepsy, the presence of abnormalities even outside the seizure, suggests an alteration of this feature. Hence, we hypothesized that the resting epileptic brain deviates from the healthy critical state. In particular, epileptogenic areas tend to exhibit higher excitability compared with “healthy” regions (Beggs and Plenz, 2003; Monto et al., 2007; Arviv et al., 2016).

Here, we reconstructed 14 epileptic and 5 healthy human BNMs with individual SC and spatial epileptogenicity using The Virtual Brain (TVB; Sanz-Leon et al., 2013, 2015), a neuroinformatics platform that enables biologically realistic modeling and simulation of brain network dynamics using connectome-based approaches and directly linking them to various brain imaging modalities. We developed a population model capable of producing physiological oscillations beyond the epileptiform activity by linearly weighting healthy and epileptiform contributions as ratios of epileptogenicity using the neural mass model Epileptor, which was originally designed to reproduce seizure dynamics (Jirsa et al., 2014). Our proposed approach results from the presumed existence of different zones with varying pathophysiological activity (Lüders et al., 2006), reflecting a gradual change of epileptogenicity. In dynamical systems, discrete onsets of novel behavior are described by bifurcations, which are parametrized by continuous control parameters. The concept of a continuous degree of epileptogenicity has been previously introduced in Jirsa et al. (2014) and applied for the characterization of epileptogenic zones in Proix et al. (2017). Here, we formalized this notion in a region-specific model framework through the extended Epileptor. We then investigated the optimal working point in epileptic BNMs against health through parameter exploration, and tested the effect of in silico IEDs and the local shift in the dynamics of epileptogenic regions. This approach allowed us to identify dynamical mechanisms that can explain the empirically observed RS-FC alterations.

Materials and Methods

Patient selection and data acquisition

We selected 14 patients in whom drug-resistant epilepsy had been diagnosed (mean ± SD age, 31.6 ± 11.3 years; 7 females), previously analyzed in Proix et al. (2017). Patients have different types of partial epilepsy accounting for different epileptogenic zone localization and underwent a standard comprehensive presurgical evaluation including structural, functional, and diffusion tensor scans. The findings from this evaluation and details of the patients are summarized in Table 1

View this table:
  • View inline
  • View popup
Table 1.

Clinical characteristics of patients with focal epilepsy

MRI scans were performed on a Magnetom Verio 3T MR-scanner (Siemens) at the Center for Magnetic Resonance in Biology and Medicine (CRMBM, Marseille, France). T1-weighted anatomical images were acquired with a 3D-MPRAGE sequence (repetition time = 1900 ms; echo time = 2.19 ms; inversion time = 900 ms; voxel size = 1.0 × 1.0 × 1.0 mm3), the diffusion MRI images used a diffusion tensor imaging (DTI)-MR sequence (repetition time = 10 700 ms; echo time = 95 ms; angular gradient set of 64 directions; b-weighting of 1000 s/mm2; 60 contiguous slices; voxel size = 2.0 × 2.0 × 2.0 mm3) and the 20 min resting-state functional MRI images recorded by a BOLD-sensitized EPI T2*-weighted sequence (350 volumes; repetition time = 3600 ms; echo time = 27 ms; voxel size = 2.0 × 2.0 × 2.5 mm; 50 slices; flip angle = 90°). During the resting-state protocol, the patients were asked to stay awake and keep their eyes closed.

Additionally, five healthy control subjects, with no history of neurologic or psychiatric disease who had undergone the same MRI protocol, were also selected to evaluate the dynamical characteristics of subject-specific RS-FC compared with patients with epilepsy and were used in a personalized BNM. All participants signed an informed consent form according to the rules of the local ethics committee (Comité de Protection des Personnes Marseille 2).

Data preprocessing

Anatomical MRI data preprocessing and structural reconstruction.

The reconstruction of the subjects' individual brain network topography and connection topology within the 3D physical space was performed using an in-house pipeline for automatic processing of multimodal neuroimaging data based on publicly available neuroimaging tools and customized for TVB (https://github.com/the-virtual-brain/tvb-recon). The current version of the pipeline used in this study evolved from a version described previously in Proix et al. (2016).

In short, the pipeline proceeds as follows for each subject: First, the brain anatomy was reconstructed from the T1-weighted images using the recon-all command from the FreeSurfer package (version 6.0.0; Fischl, 2012). The T1-weighted images and the generated parcellation volume were then aligned with the diffusion-weighted images (DWIs) using the linear registration tool flirt from FSL package (FMRIB Software Library, version 6.0; Jenkinson et al., 2012). The correlation ratio cost function was used for the alignment with 12 degrees of freedom. The tractography was performed using the tools from the MRtrix package (version 0.3.15; Tournier et al., 2012). The fiber orientation distributions were estimated from the DWI using spherical deconvolution (Tournier et al., 2007) by the dwi2fod tool with the response function estimated by the dwi2response tool using Tournier's algorithm (Tournier et al., 2013). Afterward, 15 million tracts were generated using the tckgen tool by a probabilistic tractography algorithm, iFOD2, based on a second-order integration over fiber orientation distribution (Tournier et al., 2010) with the streamlines seeded randomly within the brain volume. Finally, the SC weights matrix Cij was generated by the tck2connectome tool by counting the number of tracts connecting the regions in the Desikan–Killiany parcellation generated by FreeSurfer including 70 cortical regions and 14 subcortical regions (Desikan et al., 2006). No threshold was used to prune weaker edges; however, the BNM is sensitive to connection strength, thereby effectively discarding the effect of smaller weights. All processed data were formatted to facilitate import into TVB, and each Cij was normalized to unity as Embedded Image. Note that as the patients have small lesions (in terms of volume), the preprocessing of their MRI scan does not necessitate filling the missing portions of the brain at the lesion site by the brain tissue of the nonlesioned hemisphere homologous (see Falcon et al., 2015 and Aerts et al., 2020, for a detailed description of such procedure).

Functional MRI data preprocessing.

RS-fMRI data preprocessing was performed with the FSL feat (FMRI Expert Analysis Tool, version 5.0.9; Woolrich et al., 2001) toolbox with standard parameters and not discarding any independent component analysis (ICA) components. The preprocessing steps included the following: discarding the first six images of each scanning run to allow the MRI signal to reach steady state, high-pass temporal filtering to remove slow time drifts (100 s high-pass filter), motion and slice timing correction, spatial smoothing with Gaussian kernel (FWHM = 5 mm), brain extraction, and 12 linear registration to the MNI space.

Then, functional data were registered to the subject's T1-weighted images and parcellated according to FreeSurfer segmentation. By inverting the mapping rule found by registration, anatomical segmentations were mapped to the functional space and average blood-oxygen-level dependent (BOLD) signal time series for each region were generated by computing the spatial mean over all voxel time series of each region. Finally, to limit the effects of physiological noise, the overall time series were temporally low-passed filtered, removing frequencies >0.1 Hz (Cordes et al., 2001).

Note that global regression was not performed as it shifts the distribution of the correlation values of the RS-FC and, in particular, allows the introduction of spurious negative correlations and the underestimation of the true-positive ones (Fox et al., 2009; Murphy et al., 2009). In addition, no additional physiological noise removal was applied, although it is a standard preprocessing procedure when analyzing FC in experimental RS-fMRI studies (Bettus et al., 2009, 2010; Pereira et al., 2010). However, there is also a large number of computational works where this is not the case, even from the same laboratory groups (Schirner et al., 2015, 2018; Deco et al., 2017a; Demirtaş et al., 2017). Having the focus on the novel features of the model and how it can be applied to data fitting, we are less interested in specific activation patterns that could indeed be affected by different acquisition and processing steps, but are more interested in capturing the systematic differences appearing between groups, making the issue of physiological noise of less relevance. Nevertheless, we investigated the effect of head motion and the presence of significant group differences in framewise displacement (FD) for both translational and rotational parameters as movement metric (Power et al., 2012). For the rotational parameters, degrees of arc were converted to millimeters by calculating displacement on the surface of a sphere with a radius of 50 mm. Although epileptic patients tend to move more (mean FD = 0.3605 mm) than healthy control subjects (mean FD = 0.1823 mm) during fMRI acquisition, as observed in previous studies (Lemieux et al., 2007), the difference in mean FD is not significant (two-tailed Wilcoxon rank-sum test: p = 0.15; Ridley et al., 2015). This indicates that residual motion-related signals did not make a qualitative change in the forthcoming group results.

Individual brain model of resting-state in epilepsy using TVB.

To obtain the BOLD signal during rest, we first virtualized 14 patients with epilepsy and 5 healthy control subjects using The Virtual Brain (version 1.5.6; http://www.thevirtualbrain.org) and simulated the spontaneous neural activity of the brain. TVB is a free open-source neuroinformatics tool designed to aid in the exploration of network mechanisms of brain function and associated pathologies. TVB provides the possibility to feed computational neuronal network models with information about structural and functional imaging data including population [(stereotactic, S)EEG/MEG] activity, spatially highly resolved whole-brain metabolic/vascular signals (fMRI) and global measures of neuronal connections (DTI), for intact as well as pathologically altered connectivity. TVB is model agnostic and offers a wide range of neural population models to be used as network nodes. Manipulations of network parameters within The Virtual Brain allow researchers and clinicians to test the effects of experimental paradigms, interventions (e.g., stimulation and surgery), and therapeutic strategies (e.g., pharmaceutical interventions targeting local areas). The computational environment enables the user to visualize the simulated data in 2D and 3D and to perform data analyses in the same way as commonly performed with empirical data.

Individual virtual human brains were built based on the mutual interactions, in other words, the linear summation as specified in the current subsection, of the local brain region dynamics coupled through the underlying empirical anatomical SC matrix Cij of each subject obtained from the anatomical and diffusion MRI scans. The interareal connections were weighted by the strength specified in the SC matrix and by a common scaling factor that multiplies all interareal connection strengths. To explore interictal RS network dynamics, we developed a network node model capable of expressing both regionally specific physiological and epileptiform activity. A parameter p may scale the relative contribution of each type of activity to the overall node dynamics, as follows: Embedded Image where yepileptogenic and yhealthy are the respective epileptic and healthy neurons of the population activity Y.

Previously, regional epileptiform dynamics have been described by the Epileptor model (Jirsa et al., 2014), which was initially designed to realistically reproduce the temporal dynamics of epileptic seizures and also include a slow coupling variable that is responsible for the switching between ictal and interictal states (Fig. 1B,C). Naively, the novel neural population model could be interpreted as separated into epileptic and healthy neurons, which is possible, but certainly an oversimplification. The gradual parametrization of epileptogenicity is principally a functional differentiation, which may find various mechanistic realizations resulting in the same pathophysiology. The initial Epileptor model comprises three different time scales interacting together and accounting for various electrographic patterns. The fastest and intermediate time scales are two coupled oscillators [(x1, y1) and (x2, y2)], accounting respectively for the low-voltage fast discharges (i.e., very fast oscillations; Fig. 1A, top) and spike-and-wave discharges (Fig. 1A, middle). The slowest time scale is responsible for leading the autonomous switch between interictal and ictal states and is driven by a slow-permittivity variable z (Fig. 1A, bottom). This switching is accompanied by a direct current (DC) shift (Fig. 1A, top), which has been recorded in vitro and in vivo (Ikeda et al., 1999; Vanhatalo et al., 2003; Jirsa et al., 2014).

At some distance before and after seizures (i.e., during the interictal state), the epileptic brain appears to operate “normally” and expresses its rich dynamic repertoire of diverse brain states when driven by noise; these brain states are known as resting-state networks (RSNs; Fox et al., 2005). On fast time scales of 10–500 ms, electrographic recordings identify characteristic oscillatory modes of brain activity showing transient spindle-like behaviors (i.e., fast damped subthreshold oscillations), which repeat themselves intermittently. These waves patterns are strongly dominated by α waves (8–12 Hz; da Silva et al., 1997; Buzsaki, 2006). The existence of fast subthreshold oscillations is a distinguishable feature of systems near a Hopf bifurcation (Fig. 1D; Izhikevich, 2007). In Epileptor, a Hopf bifurcation can be configured for m = –0.5 and Embedded Image (El Houssaini et al., 2020). However, such parametrization (Embedded Image) sets the second population (x2, y2) far from its bifurcation, resulting in the loss of interictal spikes when the system is destabilized by noise. To address this problem, we extended Epileptor by another fast time scale of coupled oscillators (x3, y3) accounting for transient spindle-like patterns and mathematically equivalent to the normal form of a supercritical Hopf bifurcation (Stefanescu and Jirsa, 2008; Kuznetsov, 2013). This system was already used to retrieve RSNs (Ghosh et al., 2008; Freyer et al., 2011, 2012; Spiegler et al., 2016). Thus, the extended version of Epileptor equations read as follows: Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image where Embedded Image Embedded Image Embedded Image and m = 0, Embedded Image and γ = 0.01. The characteristic frequency rate d, fixed to 0.02, sets the natural frequency of the third subsystem, ∼10 Hz, the most powerful frequency peak observed in electrographic recordings at rest (da Silva et al., 1997; Buzsaki, 2006). The degree of epileptogenicity or excitability of a brain region is represented through the values x0 and a. If x0 is greater than a critical value (i.e., in supercritical regime), Embedded Image, the brain region can trigger seizures autonomously; otherwise, it is in its equilibrium state (i.e., in subcritical regime). In a similar way, if a is greater than a critical value, acritic = 1.74, the brain region enters in a stable limit cycle; otherwise, it is in a stable fixed point. We note that in the original Epileptor, τ0 = 2857, b2 = 2, and τ2 = 10, the main differences being that IED propagation is shifted toward the interictal period as g increases more rapidly and mean spike frequency (the number of IEDs per minute) is decreased to be more realistic. A more detailed description of the original Epileptor model can be found in Jirsa et al. (2014), with an extended bifurcation analysis of its parameters in El Houssaini et al. (2015, 2020). Recent studies link the Epileptor to physiological mechanisms of extracellular potassium accumulation (Chizhov et al., 2018). We coupled the network nodes by permittivity coupling for the Epileptor subpopulation following Proix et al. (2014), and by a fast diffusive coupling for the Hopf subpopulation. The whole-brain network activity is then described by the following equations: Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image Embedded Image where Cij are the weights of the subjects-based SC matrix, and Ks and Krs are the respective large-scale scaling parameters of the connectivity weights for the Epileptor and Hopf subpopulations. Note that as the BNM of one brain region already takes into account the effect of its internal connectivity, the connection of a region to itself was set to 0 in the connectivity matrix Cij for the simulations. Also, we assumed the neural transmission via Cij as instantaneous. Although time delays due to the tract propagation can be of crucial importance for the study of synchronization (Ghosh et al., 2008; Petkoski et al., 2016, 2018; Petkoski and Jirsa, 2019), here being on the phenomenological level, we assumed their impact to be encompassed in the neural masses and we neglected them for the sake of the computational cost. The local bifurcation parameter ai and the global scaling parameter Krs are the control parameters with which we studied, by extensive search, to find the optimal dynamical working region of the brain model where the simulations maximally fitted the empirical functional data. The global coupling parameter Ks was fixed to an arbitrary value of 0.1 such that the full system can reproduce seizure (or spike) spread pattern following the clinical criteria of each patient (Table 2). The final parameter values were chosen to fit the extended Epileptor against the experimental data, where the combination of the three ensemble x-variables was matched visually against the electrographic signatures of epileptic patient recordings. We found that plotting the following: Embedded Image where the parameter pi scales the proportion of the respective activity during the different processes (i.e., ictal and interictal period), as a function of time, bore a striking resemblance with the empirical SEEG signals. From this model, three different case scenarios can then be defined. The first one is the resting state without any epileptiform activity (or spike-free resting state). In this case, the epileptogenic population is silent and p is set to 0.1 (Fig. 2A). The second and third cases are resting state with interictal spikes (or spiking resting state) and resting state with seizures (or bursting resting state), respectively. Here, the epileptogenic population is dominant and p is set closer to 1, or in this case to 0.9 for the epileptogenic zone (EZ)/primary irritative zone (IZ1) and 0.7 for the propagation zone (PZ)/secondary irritative zone (IZ2; Fig. 2B,C). See Definition of the epileptogenic brain networks subsection for the definition of the different zones. For the purpose of our study, we focused only on the two-first case scenarios.

View this table:
  • View inline
  • View popup
Table 2.

Results of EZ/IZ1 and PZ/IZ2 prediction from SEEG signals for each patient

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Dynamics of the original Epileptor and supercritical Hopf oscillator model. A, B, Time series of the first x1 (top), second x2 (middle), and z (bottom) state variables (A) as well as the original Epileptor model (B; expressed in terms of Embedded Image) are plotted. C, The seizure trajectory is approximated in a 3D physical space defined by the state variables (y1 andEmbedded Image) and by the slow permittivity variable z. D, Depending on the local bifurcation parameter a, each region has a supercritical Hopf bifurcation at a = acritic such that, for Embedded Image the region is in a stable fixed point (subcritical regime) and the system corresponds to a damped oscillatory state, whereas for Embedded Image the region enters in a stable limit cycle and the system switches to an oscillatory state (supercritical regime). The closer a node operates to the critical point, the larger and longer lasting is the oscillation (compared ψ1 and ψ2). When the critical point is reached, the node intrinsically performs a rhythm of constant magnitude (see ψ3).

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Extended Epileptor dynamics and spatial distribution of epileptogenicity. A–C, Time series of the three different case scenarios of the Extended Epileptor model, namely resting state without any epileptiform activity (with p = 0.1 in Eq. 19 and x0 = –2.5 in Eq. 14; A), resting state with interictal spikes (here with p = 0.7 and x0 = –2.07 as example; B), and resting state with seizures (here with p = 0.7 and x0 = –1.9 as example; C). D, When embedded into the network and in the spike-free resting-state scenario (A), the regions set in the EZ/IZ1 (red nodes) and PZ/IZ2 (blue nodes) have a very low excitability value x0,i (Embedded Image) and high value ai (ai > acritic), and all other regions set in the NIZ (green nodes) have the same excitability value x0,i but an elevated excitability value ai (ai = acritic); in the spiking resting-state scenario, the regions in the EZ/IZ1 have a low x0,i (Embedded Image) and high ai (ai > acritic), the regions in PZ/IZ2 have lower x0,i (Embedded Image) and high ai (ai > acritic), and all other regions in the NIZ have a very low x0,i (Embedded Image) and elevated ai (ai = acritic; B). The white links represent the anatomical structural connections. See Definition of the epileptogenic brain networks for the description of the different zones.

Numerical implementation

All the simulations were performed with TVB using a stochastic Heun's integration scheme (Mannella, 2002) with a time step of 0.1 ms, a simulation time of 20 min, and random initial conditions drawn from a normal distribution (Embedded Image). Additive white Gaussian noise was introduced in the state variables x2 and y2, with a mean of 0 and noise strength of 0.00025, as well as in x3 with a mean of 0 and noise strength of 0.02. The data corresponding to the first 20 s were always discarded from the analysis to avoid initial transient dynamics.

To test the emergence of ultra-slow fluctuations, we estimated the BOLD signal changes associated with the simulated neural activity (Eq. 19) using the Ballon–Windkessel hemodynamic model (Friston et al., 2003) implemented in TVB. The Ballon–Windkessel model describes the coupling of perfusion to the BOLD signal, with a dynamical model of the transduction of neural activity into perfusion changes. The simulated BOLD signal was downsampled at 3600 ms to match the time resolution (TR) of the empirical fMRI signals. The global mean signal was then regressed out from each region's time series, and temporal bandpass filtering was performed to retain frequencies between 0.01 and 0.1 Hz using a third-order Butterworth filter to reproduce empirical conditions.

Definition of the epileptogenic brain networks

Epileptogenic brain networks were evaluated by two different methods. The first one consisted of the visual inspection and interpretation by the expert epileptologist (F.B.) of the different measurement modalities gathered throughout the two-step procedure (noninvasive and invasive) of the comprehensive presurgical evaluation of each patient. The second method consisted in the application of signal-processing techniques on the invasive SEEG measurements, which have been used in previous studies (Bartolomei et al., 2013b, 2017). In particular, SEEG signals were used to refine the zones involved by different epileptogenic processes (Bettus et al., 2011), as follows: the EZ/IZ1, the PZ/IZ2, and the noninvolved zone (NIZ) are determined for each patient. EZ/IZ1 is defined as the subset of brain regions involved in the generation of seizures that may also exhibit IEDs. PZ/IZ2 is defined as those regions only secondarily involved in seizures and that produce interictal spikes. Finally, NIZ is defined as structures without epileptiform discharges during clinical monitoring. The identified zones for each patient are provided in Table 2.

For the simulations, the different regions set respectively in EZ/IZ1 and PZ/IZ2 were used to reconstruct the epileptogenic networks in our individual virtual epileptic brains. We defined a spatial distribution map of epileptogenicity or excitability where each brain region, i, was characterized by an excitability value, x0,i, which quantifies the ability of the Epileptor subpopulation to trigger an epileptogenic discharge or not; and an excitability value ai, which quantifies the ability of the Hopf subpopulation to generate self-sustained oscillations. The spatial map of epileptogenicity comprises the excitability values of the EZ/IZ1, PZ/IZ2, and all other NIZ regions (Fig. 2D).

Resting-state functional connectivity analysis

FC describes brain function through estimates of the covariant links between two signals (originating from different brain regions), reflecting how different brain areas coordinate their activities. Functional connections were explored from both local (i.e., at the regional level) and global (i.e., at the whole-brain level) perspectives using a set of four measures of their spatiotemporal dynamics that have been widely used to estimate FC from RS-fMRI signals (for review, see Smitha et al., 2017). An additional metric was used to explore the complexity of the fMRI signals. Figure 3 shows a summary of these metrics. The simulated and empirical data are analyzed using TVB analysis tools.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Methods for measuring and fitting RS-FC. A, sPCC is estimated as the Pearson correlation coefficients between each pair of brain regions over the whole time window. The fitting of the sPCC is measured by the Pearson correlation coefficient between corresponding elements of the upper triangular part of the matrices. B, dPCC is estimated as the series of sPCC matrices from the windowed segment. C, FCD is estimated as the Pearson correlation coefficient between each pair of sPCC matrices. D, For comparing the FCD statistics, the upper triangular elements of the matrices were collected, and the simulated and empirical distribution were compared by means of the Kolmogorov–Smirnov distance between them. E, The phase synchronization of the BOLD signals is measured by the Kuramoto order parameter R. F, The coherence and the metastability are estimated as the mean Embedded Image and the standard deviation σR of the Kuramoto order parameter across time. The phase synchrony is fitted by minimizing the absolute value of the difference between the empirical and simulated coherence and metastability. G, Schematic illustration of the coarse-graining (left) and an exemplary time series is shown to illustrate the procedure for calculating sample entropy for the case m = 2 and a given positive real value r (denoted by the height of the colored bands; right). H, The MSE profiles are fitted by minimizing their RMSD.

Local functional connectivity

Static PCC.

The first metric applied is the most classic and widely used to infer the strength of functional connections, namely the Pearson correlation coefficients (PCCs). The PCC consists of estimating the (linear) temporal correlations between each pair of brain regions over the whole time window acquisition [we talk about “static” PCC (sPCC); Bandettini et al., 1993; Biswal et al., 1995]; that is, for each pair of regions i and j, the corresponding time series xi(t) and xj(t) of size N are used to calculate the correlation coefficients rij (Fig. 3A), as follows: Embedded Image where the notation Embedded Image denotes the mathematical mean operator. This sPCC was estimated in the same way from the empirical and simulated fMRI BOLD signals. Since the matrices are symmetric, we compared the modeled and empirical sPCC matrices by using the Pearson correlation coefficient between the corresponding elements of their upper triangular part.

Dynamic PCC.

To consider the temporal dynamics of RS-FC, we computed the FC using a sliding temporal window approach (Hutchison et al., 2013; Allen et al., 2014; Hansen et al., 2015). In this approach, each full-length BOLD signals of 20 min (either empirical or simulated) was split up into a mean of 328 sliding windows of 60 s, overlapping by56.4 s (i.e., an increment of 1 TR). For each sliding window centered in time t, we calculated a separate sPCC matrix [Embedded Image]. This procedure resulted in a series of sPCC matrices that describe the time-resolved behavior of connectivity over the entire duration of the fMRI BOLD signals, namely the dynamic PCC (dPCC; Fig. 3B). Subsequently, to study the evolution of the dPCC over time with reduced dimensionality, we computed a time-versus-time matrix representing the functional connectivity dynamic (FCD; Hansen et al., 2015), where each entry (ti, tj) was defined by the Pearson correlation coefficient between the upper triangular parts of the two matrices Embedded Image and Embedded Image centered at time ti and tj, respectively (Fig. 3C).

For comparing the FCD statistics, we collected the upper triangular elements of the matrices and generated their cumulative distribution function. Then, we compared them through the Kolmogorov–Smirnov (K–S) distance, which quantifies the maximal difference between the two samples (Fig. 3D; Deco et al., 2017a).

Global functional connectivity

The global network dynamics were quantified using coherence and metastability, which are defined by the mean and the standard deviation of synchronization over time (Shanahan, 2010; Wildie and Shanahan, 2012; Cabral et al., 2014). The phase synchrony of the empirical and simulated network dynamics was evaluated using the Kuramoto order parameter R(t) (Kuramoto, 1984), which is defined as follows (Fig. 3E): Embedded Image where N is the number of brain regions in the network and Embedded Image is the instantaneous phase of each region j estimated using Hilbert transform (Glerean et al., 2012; Ponce-Alvarez et al., 2015). R measures the phase uniformity and varies between 0 for a network in a fully desynchronized—or incoherent—state and 1 for full synchronization. Thus, for calculating the coherence Embedded Image and metastability σR of the empirical and simulated BOLD signals, we computed the instantaneous phase Embedded Image of the signal of each region j using the Hilbert transform. The Hilbert transform yields the associated analytical signals s(t), representing a signal in the time domain as a rotating vector with an instantaneous phase, Embedded Image, and an instantaneous amplitude, A(t) [i.e., Embedded Image]. The phase and amplitude are given by the argument and the modulus, respectively, of the complex signal z(t), given by Embedded Image, where i is the imaginary unit and Embedded Image is the Hilbert transform of s(t).The simulated coherence and metastability were compared with the empirical values of synchronization between the different brain regions across time by using as a measure of similarity the absolute value of the difference (Wildie and Shanahan, 2012; Deco et al., 2017b; Lee and Frangou, 2017).

Complexity analysis: multiscale entropy

Introduced as a measure of complexity of physiological signals, multiscale entropy (MSE; Costa et al., 2002, 2005) is thought to evaluate the presence of long-range correlations in the behavior under scrutiny by taking into consideration the multiple time scales on which the overall dynamics operate. Technically, MSE consists in computing sample entropy (Richman and Moorman, 2000) of the successively coarse-grained signals, namely downsampled versions of the original signal at different (time) scales. For a detailed description of MSE measure and its relevance for the analysis of signal complexity see Courtiol et al. (2016).

Briefly, a time series Embedded Image of size N is divided into segments that correspond to consecutive nonoverlapping time intervals of length τSF, where τSF represents the scale factor and takes integer values ≥1. All data points in each segment are replaced by the average value of that segment, thus producing a new time series, Embedded Image, called a coarse-grained time series, the length of which is equal to that of the original, N, divided by the scale factor τSF (Fig. 3G), as follows: Embedded Image

For scale 1, y1 is simply the original time series. Then, the sample entropy (SampEn), quantifying the predictability or regularity of a time series (Richman and Moorman, 2000), is applied for each time series Embedded Image. It is defined as the negative natural logarithm of the conditional probability that two similar sequences of m consecutive data points in a time series of length Embedded Image within a given tolerance r normalized to the standard deviation of the time series, will remain similar when the next point m + 1 is also included in the sequence (Fig. 3G). Then, for each time scale τSF: Embedded Image where the quantity Embedded Image is the conditional probability described above. Note that BOLD time series usually comprise few data points, and the coarse-grained procedure in MSE with large scale factor may result in short data length and, subsequently, in unreliable SampEn estimation. Prior studies suggested that data length of 10m to 20m should be sufficient for a robust calculation (Richman and Moorman, 2000). Therefore, following several studies using MSE in fMRI recordings, we set m to 1 and r to 0.35 (Nakagawa et al., 2013; Yang et al., 2013) and computed MSE over scales 1:13 (3.6–46.8 s). The empirical and simulated MSE profiles were compared using the root mean square distance (RMSD; Fig. 3H).

Global similarity

Based on the previously described five fitting metrics and having the considerations that (1) absolute value of the difference is better as it gets closer to 0, (2) the Pearson correlation coefficient is better as it gets closer to 1, (3) K–S is better as it gets closer to 0, and (4) RMSD is better as it gets closer to 0, we developed an additional fitting metric, the global similarity (GS), to express all these conditions in a single numerical value (Kehoe et al., 2017). MSE is the only metric that portrays the complexity, hence we gave more importance to RMSD in the expression of GS. Thus, GS depends quadratically on RMSD, while it depends only linearly on the other metrics, as follows: Embedded Image where Coh and Meta represent the absolute value of the difference between the empirical and simulated coherence and metastability, respectively; corr is the Pearson correlation coefficient between the empirical and simulated sPCC matrices, and KS is the K–S distance between the empirical and simulated FCD histograms. Then, we used GS to find the optimal working region in an automated way by defining it as the values of the free model parameters, namely the local bifurcation parameter a and the global coupling strength Krs, where GS exhibited its minimal value. Note that the values of GS were normalized (between [1,2]).

Statistical analysis

When comparing two SC or RS-FC distributions from a particular measure of connectivity, we performed a two-tailed Wilcoxon rank-sum test, testing the null hypothesis that there existed no difference between the calculated measures at both distributions. The statistical significance level was set to p < 0.05, and so the z score was greater than Embedded Image.

Results

The main goal of our study was to identify network mechanisms behind the RS-FC patterns observed in epileptic patients, and to characterize how epileptogenic brain regions express themselves outside the seizure (i.e., during the interictal state). To this end, we developed a novel large-scale BNM linking the underlying anatomical SC of each patient (derived from MRI scans) with the local functional dynamics of each brain area to emulate the characteristics of spontaneous whole-brain dynamics, as observed in functional neuroimaging data.

The spatiotemporal structures of spontaneous fMRI BOLD fluctuations were characterized using the following five RS-FC estimates: (1) sPCC and (2) dPCC (or FCD; the static PCC describes the mean spatial structure of the resting-state activity, whereas the dynamic PCC captures the temporal structure of those spatial correlations); (3) coherence and (4) metastability, which quantify the level of phase synchronization between brain regions across time; and (5) MSE, which measures the level of complexity within each brain region across multiple temporal scales.

Structural reorganization in epileptic patients not captured by SC weights

The SC from the brains of epileptic patients and healthy control subjects were compared using two typical graph theory measures, namely node strength and streamline counts. These metrics can provide good predictions of the optimal working point of the brain model, as SC weights appear in the brain model equations as a linear combination with the remote brain region activity. The mean values of the streamline counts and mean structural strength of the nodes did not differ significantly between the two groups, and also between patients with MRI-positive and MRI-negative results. However, the small sample size of the healthy group limits the statistical power of this analysis.

Shift of the optimal working point of brain model toward lower excitability in patients with epilepsy

We first investigated the dynamical properties of the optimal working point of the brain model that was able to fit the characteristics of the empirical fMRI data from 14 epileptic patients and 5 healthy control subjects. fMRI BOLD activity was simulated with the extended Epileptor model in the case of absence of epileptiform discharges (i.e., in the spike-free resting-state scenario with p = 0.1 for all brain regions in Eq. 19) with the brain regions coupled through the empirical subject-specific SC matrix Cij as extracted from the DTI-based tractography (normalized to unity). Using our model, we studied whether the patients' SC led to altered RS-FC at the dynamical level, and, particularly, how this depended on the two free parameters of the model, namely, the local bifurcation parameter ai, that was homogeneously modified over all regions (i.e., ai = a, for all i), and the global coupling Krs. We performed a parameter space exploration of these two parameters and characterized the whole-brain RS-FC of each subject. For this, we used the previously introduced metrics static and dynamic PCC, coherence, metastability, and MSE for each set of parameters (a, Krs). In this paragraph, we additionally used the GS to refine the working point. These six metrics characterize computationally the bifurcation properties of the full network dynamics. Then, we computed the best fit between the simulated and empirical RS-FC-based metric values. All five parameter spaces presented in Figure 4A–C display a full picture of the spatiotemporal organization of the system.

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Whole-brain parameter space explorations and RS-FC fitting for an exemplary epileptic patient and healthy control subject. A–C, The six metrics for assessing the RS-FC fitting between the simulated and empirical data are shown color coded as a function of the global coupling Krs (x-axis) and the bifurcation parameter a (y-axis). Top row, Correlation coefficient between the empirical and simulated sPCC matrices (left) and the K–S distance between the empirical and simulated FCD histograms (right). Middle row, Absolute value of the difference between the empirical and simulated coherence (left) and metastability (right). Bottom row, RMSD between the empirical and simulated MSE matrices (left) and GS metric (right). The horizontal black dashed lines represent the critical point. B–D, Simulated sPCC, FCD, and MSE matrices for the optimal working point. For comparison, the same matrices are plotted for the empirical data (denoted by*).

Figure 4A shows the results of this analysis for an exemplary epileptic patient (P1). As can be observed in Figure 4A (top row, left panel), the correlation between the empirical and simulated sPCC matrices is sensitive to the large-scale coupling strength Krs and the bifurcation parameter a. The model shows the best agreement with empirical data for a broad range of parameters (Fig. 4A, area indicated with hot colors); indeed, a large region of model parameters is consistent with the empirical data and reaches a correlation of up to 0.40 between empirical and simulated sPCC matrices. This is within the values of similarity reported in previous studies that have simulated subject-specific brain dynamics at the large-scale level (see for example Deco et al., 2013). In Figure 4A (right), the best fitting of the spatiotemporal characteristics of the empirical RS-fMRI data can be found at the minimum of the K–S distance between the empirical and simulated FCD histograms for a smaller range of parameters (Fig. 4A, area indicated with cold colors). In Figure 4A (middle row), the phase synchronization behavior of the system obtained from the simulations is compared with the empirical fMRI BOLD data of the patient for the same range of model parameters. In the range of best agreement with FCD (Fig. 4A, cold colors), the absolute difference between the empirical and simulated coherence (Fig. 4A, left), and metastability (Fig. 4A, right), respectively, is minimized (Fig. 4A, area also indicated with cold colors). Note that the calculated coherence (or degree of phase synchrony) is moderate (Embedded Image), indicating that the system is neither fully synchronized (Embedded Image close to 1) nor incoherent (Embedded Image close to 0), and the coherence exhibits maximal variability or metastability (Embedded Image). Moreover in the last row of Figure 4A, the RMSD between the empirical and simulated MSE matrices displays the best agreement for a limited range of parameters close to the Hopf bifurcation (Fig. 4A, horizontal black dashed line). In Figure 4A (right), the GS is plotted to refine the best-fit zone between all metrics (area indicated with cold colors). All other patients exhibit RS-FC configuration patterns that were highly similar to those of patient P1 (Fig. 5A, median group results of GS).

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

Group comparison of the optimal working point in the brain models. Median group results. A, B, Level of fitting of the GS between the simulated and empirical data are shown color coded as a function of the bifurcation parameter a (y-axis) and the global coupling Krs (x-axis) for the epileptic group (EP; A) and the healthy group (HC; B). The horizontal black dashed lines represent the critical point. C, Minimal level of fitting of the GS as a function of the local bifurcation parameter values a (left) and coupling strength Krs (right), for the epileptic group (red curves) and control group (green curves). Areas of faded colors represent the standard error intervals. D, Distribution of optimal local bifurcation parameter values a (left) and coupling parameter values Krs (right), for EPs (red violin plot) and HCs (green violin plot). *p = 0.038.

Subsequently, the same analysis was performed on the simulated data obtained from the healthy control subjects. The results of this analysis are presented in Figure 4C for an exemplary healthy subject. The plots show that the data obtained with the model indicate an optimal working point into a similar range of parameters compared with patient P1. All other subjects have very similar fitting patterns (Fig. 5B, median group results of GS), although the best agreement with the data is obtained for a larger range of parameters, in particular in the supercritical regime (i.e., self-sustained oscillations).

Furthermore, we compared the dynamic working point of the two groups with regard to the GS parameter space in Figure 5C. The GS between the empirical and modeled data is more sensitive to the local bifurcation parameter a (Fig. 5C, left) than the global coupling strength Krs (Fig. 5C, right) for both groups. In particular, for the epileptic patients, the optimal working point is obtained when the brain regions operate in the vicinity of the supercritical Hopf bifurcation (i.e., close to Embedded Image; Fig. 5D, left). For the healthy subjects, the best fit is obtained in the supercritical regime (Fig. 5D, left). However, there is a quite large region with a U-shaped dependence of Krs on a, where there is relatively good fitting (Fig. 5B). This could suggest either that there are several working points for the healthy brain or that it could be a result of the small number of subjects. The dynamic working point shifts significantly (p = 0.038) to lower bifurcation parameter values in epilepsy compared with the healthy state, where the system remains in the oscillatory region. The simulated data obtained at the optimal working point for the exemplary epileptic patient P1 (i.e., at a = 1.74 and Krs = 20) and the healthy control subject (i.e., at a = 1.74 and Krs = 20) are shown in Figure 4B–D. For comparison, the same matrices are plotted for the empirical data.

Note that due to the lack of simultaneous physiological recording during the RS-fMRI acquisition, no additional physiological noise modeling has been applied and some aliasing artifacts (e.g., cardiac and respiratory rhythms) could be still present in the data. However, regressing them out also removes a meaningful part of the brain signal. Indeed, we performed an ICA-based noise removal for all of our subjects and the best-fit GS decreases by ∼54% for the epileptic group and by ∼70% for the control group, although the landscape and working region of the model remain topologically identical and quantitatively similar. Thus, this preprocessing procedure does not make a qualitative change in the results. The denoising of the RS-fMRI is still being debated in the field (Bright and Murphy, 2015) and more in-depth investigation of this issue is beyond the scope of our study here, which has a modeling focus.

IEDs do not change connectivity patterns

The impact of the presence of IEDs in brain signals on RS-FC features was systematically evaluated by simulating the spontaneous activity with the extended Epileptor in the spiking resting-state scenario (where the parameter p was set according to the relative contribution of healthy and interictal activity in Eq. 19). The local bifurcation parameter ai was fixed for all the regions at the value ai = a = 1.74, where the optimal similarity was observed between the model and the data in a large majority of patients, while the global coupling Krs was set at the patient-specific optimal value.

We first introduced one spiking region in the BNM and we investigated whether the spike frequency, modulated by the parameter b2 in Equation 15, led to changes in RS-FC at the whole-brain dynamical level and, in particular, how this depended on the structural strength of that region. Parameter space exploration was performed by successively modifying the selected spiking region [called the region of interest (ROI)] in the SC that produced IEDs, and by varying the parameter b2. For each set of parameters (ROI, b2), we calculated the changes in whole-brain RS-FC compared with a control case scenario (i.e., without interictal spikes or in spike-free resting state) using the metrics defined above. The results are shown in Figure 6, where the BNM was computed from the connectome of P1. Figure 6A shows that the RS-FC patterns are very slightly affected by the presence and the frequency of interictal spikes, which is increased for lower b2, as illustrated in Figure 6B, and by the structural node strength of the region. Indeed, all the parameter spaces are fully homogeneous. We then performed a deeper analysis by comparing the sPCC connectivity of the selected spiking region with the rest of the brain to the control scenario. The comparison was done by applying a two-tailed Wilcoxon rank-sum test (Fig. 7A), which reveals a widespread significant decrease in the connectivity (p = 0.04). This means higher connectivity in the control case, as the frequency of interictal spikes decreases for higher b2 (i.e., when the frequency of spikes becomes comparable to the sampling frequency of BOLD signals).

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6.

Impact of interictal spikes on the whole-brain RS-FC. A, C, The five metrics for assessing changes between simulated data with and without IEDs in brain signals are shown color coded as a function of spike frequencies linked to the parameter b2 (x-axis) and spiking region ROI sorted by structural node strength (A) and the size of the patient-specific IZ (y-axis; C). Top row, Correlation coefficient between the two simulated sPCC matrices (left) and the K–S distance between the two simulated FCD histograms (right). Middle row, Absolute value of the difference between the two simulated coherence (left) and metastability (right). Bottom row, RMSD between the two simulated MSE profiles. B, Left, Number of spikes (in minutes) in the mid-connected brain region and (right) its respective neural activity time series for two different values of the parameter b2 (black dashed lines in A).

Figure 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 7.

Statistical comparison of the sPCC matrices between the spiking or hyperexcitable case to the control scenario. The z-score values are color coded as a function of spike frequency b2 (top row) or bifurcation parameter a (x-axis; bottom row), and spiking or hyperexcitable region ROI sorted by structural node strength (left column) or the size of the patient-specific IZ (y-axis; middle and right columns). A–F, Comparisons are for spiking or hyperexcitable region with the rest of the brain (A–D), within IZ (B–E), and IZ with the rest of the brain (C–F).

Furthermore, we analyzed the impact of the presence of a spiking subnetwork on the RS-FC, by simulating individual epileptic brain models with an IZ. The epileptogenicity x0,i of the regions belonging respectively to IZ1, IZ2, and NIZ was set according to clinical criteria (Table 2). The spatial distribution of epileptogenicity values is thus heterogeneous across the brain network (Fig. 2D), with a small value x0 for regions in IZ1 (Embedded Image), smaller excitability for regions in IZ2 (Embedded Image), and very low excitability for the other regions (Embedded Image). The BOLD signal was then derived from the neural activity (Eq. 19) using the Ballon-Windkessel model, on which we calculated the RS-FC metrics. The combined results for all the patients are presented in Figure 6C. Although increasing the size of IZ increases the impact of spikes, changes are very small for all the metrics. Notably, higher spike frequencies seem to be completely irrelevant for low-frequency BOLD RS-FC. Subsequently, we compared sPCC connectivity of patient-specific IZ with the rest of the brain (Fig. 7C), and within IZ (Fig. 7B) to the control case scenario (i.e., no spiking subnetwork). Figure 7B and C, shows a significant decrease in connectivity in the spiking scenario (p = 0.00012 and Embedded Image, respectively), albeit increasing the size of the IZ increases the impact of spikes. Note that for the case of two regions in the IZ (Fig. 7B), statistical comparison cannot be performed.

Finally, we assessed the level of fitting between the empirical and modeled data including the patient-specific IZ (Fig. 8A, median group results). The combined results for all the patients show that the measures of similarity are almost unaffected when varying the frequency of the spikes, which is in agreement with the results of the systematic analysis described above and that are shown in Fig. 6C.

Figure 8.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 8.

RS-FC fitting for all the epileptic patients. A, B, Median group of the five observables for assessing the level of fitting between empirical and simulated data are shown as a function of spike frequency b2 (A) and local bifurcation parameter a (B). Left column, Correlation coefficient between the empirical and simulated sPCC matrices (light curve) and the K–S distance between the empirical and simulated FCD histograms (dark curve). Middle column, Absolute value of the difference between the empirical and simulated coherence (dark curve) and metastability (light curve). Right column, RMSD between the empirical and simulated MSE profiles. Areas of faded colors represent the standard error intervals.

Perturbations of local brain dynamics in epileptogenic brain regions trigger network-wide connectivity

The alteration of the individual local dynamics in the epileptogenic brain regions is a potential mechanism that could explain the observed empirical differences in RS-FC features between the different brain regions in epileptic brains. To test this hypothesis, we conducted a computational experiment similar to that described in the previous paragraph. Specifically, based on our model in the spike-free resting-state scenario (with p = 0.1 for all brain regions in Eq. 19), the local bifurcation parameter ai of the epileptogenic areas was systematically modified toward the supercritical regime (i.e., self-sustained oscillations), implying a hyperexcitability of that region, to characterize changes in global and regional connectivity due to a local shift in optimal dynamics.

One hyperexcitable region was first introduced in the BNMs, and we examined the effect of the level of excitability a on the resulting simulated RS-FC values and, in particular, how this depended on the structural strength of that region. Parameter space exploration was performed by successively modifying the selected hyperexcitable region (ROI) in SC, and by varying the local bifurcation of that region. For each set of parameters (ROI, a), we calculated the changes in whole-brain RS-FC compared with the control case scenario (i.e., ai = a = 1.74, for all i) using the metrics defined before. The results, presented in Figure 9A, where the BNM was computed from the connectome of patient P1, show that a hyperexcitable node can induce fluctuations in RS-FC patterns at the whole-brain level, in particular, if this node is at least moderately anatomically connected in the network and/or for very large divergence from criticality, that is increased for higher excitability a, as illustrated in Figure 9B. We then analyzed how the sPCC connectivity of the hyperexcitable region with the rest of the brain was impacted compared with the control case. Figure 7D reveals a widespread significant increase in connectivity in the hyperexcitable scenario (p = 0.04), as the bifurcation parameter a diverges from criticality.

Figure 9.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 9.

Impact of local hyperexcitability on the whole-brain RS-FC. A, C, The five observables for assessing changes between simulated data with and without hyperexcitable nodes in the network are shown color coded as a function of (y-axis) hyperexcitable region ROI sorted by structural node strength (A) and the size of the patient-specific hyperexcitable subnetwork IZ, and (x-axis) local bifurcation parameter a (C). Top row, Correlation coefficient between the two simulated sPCC matrices (left) and K–S distance between the two simulated FCD histograms (right). Middle row, Absolute value of the difference between the two simulated coherence (left) and metastability (right). Bottom row, RMSD between the two simulated MSE profiles. The red dashed lines represent the critical point. B, Neural activity time series for two different values of the parameter a for the mid-connected brain region (A, black dashed lines). D, Simulated sPCC, FCD, and MSE matrices for the optimal parameter a from the patient P1 (right). For comparison, the same matrices are plotted for the empirical data (left) and for the simulated control data (middle).

Furthermore, we investigated the impact that a hyperexcitable subnetwork could have on the RS-FC patterns, by simulating individual epileptic brain models with patient-specific IZ (including IZ1 and IZ2) replaced by a hyperexcitable zone. As shown in Figure 9C, changes spread faster (i.e., for smaller divergence from criticality) as the size of the network increases, although sPCC and MSE are less impacted than other metrics. Subsequently, the comparison between the sPCC connectivity of the hyperexcitable patient-specific subnetwork with other regions (Fig. 7F), and within the hyperexcitable subnetwork (Fig. 7E), to the control scenario show a significant increase in the connectivity in the hyperexcitable case (p = 0.04 for both panels).

Finally, we assessed the fitting between the empirical and simulated data for all the patients (Fig. 8B, median group results). As can be observed in Figure 8B (left), the correlation coefficient between the empirical and simulated sPCC matrices (Fig. 8B, light blue curve) is increased as the bifurcation parameter a diverges from criticality, and, in a similar way, the K–S distance (Fig. 8B, dark blue curve) decreases for higher a. In contrast, in Figure 8B (middle), the level of similarity between the empirical and simulated coherence (Fig. 8B, dark blue curve) slightly increases as a becomes higher and decreases for larger a; while it remains almost similar for the metastability (Fig. 8B, light blue curve). Regarding the complexity in Figure 8B (right), the difference between the empirical and simulated MSE profiles also slightly increases for larger parameter a. The simulated sPCC, FCD, and MSE matrices obtained at the best-fit point for patient P1 are shown in Figure 9D (depicted by *). For comparison, the same matrices are plotted for the empirical data (Fig. 9D, left), and for the control scenario (Fig. 9D, middle).

Discussion

Understanding and characterizing the dynamics of epileptogenic networks within the interictal periods are critical for improving the surgical management of drug-resistant partial epilepsy. To identify the relevant variables contributing to the differences observed in brain imaging between epileptic and healthy brains, we analyzed the organization of brain network activations. These variables are physiologically unspecific, as multiple mechanisms can lead to the same behavior, but are useful in research and clinical applications. Examples include the epileptogenicity of a brain region, which maps directly between a model parameter and clinical hypothesis, and the global coupling of the brain. The former is an excitability parameter that describes the interplay between oscillations and noisy dynamics, while the latter sets the range of coactivation between brain regions and global synchrony. Using the BNM, we demonstrated that epileptic brains operated at a different global working point induced by a shift toward lower excitability and less oscillatory activity compared with healthy control subjects. We also showed that the presence of in silico IEDs did not modify BOLD fluctuations and connectivity patterns compared with the spike-free signals. Moreover, we established that the epileptogenic regions displayed higher excitability and more oscillatory activity, enhancing the level of similarity with the empirical data.

To this end, we extended the phenomenological neural mass model of partial seizures, Epileptor (Jirsa et al., 2014), tuned to express regionally specific physiological oscillations in addition to the epileptiform discharges. Hence, reflecting RS activity that closely resembles SEEG recordings. The extension was made using linearly coupled oscillators near a supercritical Hopf bifurcation to model the spontaneous local field potential-like signal (Ghosh et al., 2008; Deco et al., 2017a). Although the parameters of the model cannot be directly interpreted in biological terms, the degree of epileptogenicity or excitability of each brain region, characterized by the continuous control parameters x0 and a, was motivated by the presumed existence of the different SEEG-defined zones that have been linked to the surgical prognosis (Proix et al., 2017).

Through systematic parameter exploration and fitting to neuroimaging data, we showed that the region where the brain model fits the data best lies in the supercritical region in health, characterized by self-sustained oscillations. In epilepsy, it shifts at the edge of the bifurcation, characterized by a mixed noisy oscillatory behavior. This shift toward lower excitability implies that the neural excitation during the interictal periods free of epileptiform discharges is also affected, but opposed to the (hyperexcitable) ictal state (Scharfman, 2007). These findings are in agreement with recent studies reporting a divergence from the optimal healthy bifurcation parameter in disease, such as in schizophrenia and Alzheimer's disease (Cabral et al., 2013; Demirtaş et al., 2017). This indicates that in pathologic brains, RS emerges from a different dynamical regime compared with healthy brains. The shift may result from dynamic compensatory processes that preserve neurophysiological functioning in the presence of hyperexcitable pathologic epileptiform network activity. Other influences should be considered as well, like the effect of antiepileptic drugs, indeed, patients remained on their daily medication during the scan (Wandschneider et al., 2014) or vigilance fluctuations, as sleep is characterized by a shift toward lower excitability compared with wake RS (Jobst et al., 2017). Yet, patients are more alert than control subjects due to the higher anxiety level of the investigation. Furthermore, we demonstrated that the optimal global coupling was similar in epileptic and healthy groups. The result suggests that either the neural transmission is preserved in the disease, or region modeling of the cortex and subcortical areas does not capture the structural reorganization in epilepsy (Besson et al., 2017). This is corroborated by the nonsignificant difference found between the SC weights of the two groups. Also, it should be noted that varying the weights has less influence on the propagation activity than the topology of the SC matrix (Proix et al., 2017), demonstrating the predictive power of individual connectome-based models. Another important finding is that the best fitting of the complexity in BOLD signals using MSE was present only in a narrow range of parameters when a is close to the bifurcation. In other words, we demonstrated a better way of constraining BNMs, rendering the fitting process more complex and complete.

The simulation of IEDs in brain signals showed that the whole-brain RS-FC was only slightly modified when comparing to spike-free signals, and the patient-specific connectivity was almost unaffected when adding the effect of IZ. These results were independent of the spike frequency and the SC node strength of the involved areas, suggesting that the FC in epileptic brain is related to neuronal activity other than IEDs. All patients having long-standing drug-resistant epilepsies, so we can hypothesize that the epileptogenic activity may evolve from IEDs related to more permanent connectivity changes as a result of plasticity triggered by the pathologic activity, causing metabolic changes in the BOLD data. This hypothesis is corroborated by some studies showing changes over time of the size of EZ (Bartolomei et al., 2008) or alterations in SC (Bonilha et al., 2006). Our results are consistent with those of previous studies investigating connectivity patterns in the presence/absence of IEDs, and using ROI-based or ICA-based FC methods, showing that the features of the brain networks were similar in fMRI or SEEG datasets with and without IEDs (Bettus et al., 2008; van Houdt et al., 2015). We extended these results by considering the nonstationarities within the brain signals using dynamic FC and MSE analyses. We showed limited to no effect of IEDs on RS-FC changes.

The alteration of the local dynamics of the epileptogenic regions showed that whole-brain RS-FC changes to local hyperexcitability can be widespread, including regions local to and distant from the presumed focus when compared with control signals. The performance of the model was also improved, implying that the pathologic areas operate in the supercritical regime. Several works revealed similar dynamical properties by quantifying local long-range temporal (auto)correlations using detrended fluctuation analysis (Parish et al., 2004; Monto et al., 2007; Tavares et al., 2015). They suggest a persistent network abnormality due to exposure of the neuronal networks to the epileptic activity. We can then speculate that, even in the absence of visible epileptic discharges, the specific alteration of local high excitability might help to localize the epileptogenic regions in focal epilepsy by separating these areas to healthy ones, and to inform the likely outcome after surgery. Note that we had no intention to select an optimal value for each region, but rather to mimic a scenario where epileptogenic brain areas exhibited higher excitability or hyperexcitability in parallel.

To test for consistency with prior RS-fMRI findings mainly based on (linear) temporal correlation (i.e., sPCC), we compared the direction of the alterations induced in the BNMs either by the presence of IEDs or hyperexcitability within the epileptogenic networks. In line with previous work (Bettus et al., 2009, 2010; Su et al., 2015), our results showed a significant hypoconnectivity within the spiking networks, and in their connections with the rest of the brain. We also observed a significant hyperconnectivity within the hyperexcitable networks, and in their connections with the rest of the brain, corroborating empirical results (Pereira et al., 2010; Pittau et al., 2012; Wirsich et al., 2016). This supports the hypothesis that long-standing epileptic activity results in a more permanent connectivity change. Moreover, the variation of RS-FC induced by nonstationarities during the interictal RS period was recently investigated via sliding window analyses reporting both increased (Laufs et al., 2014; Pedersen et al., 2017; Ridley et al., 2017) and decreased (Nedic et al., 2015) connectivity in the epileptogenic networks. It is quite probable that these discrepancies within the same modality could be explained by the different methodologies used or also by the vigilance fluctuations that affect RS-FC estimates (Tagliazucchi and Laufs, 2014).

Note that the results obtained in this study are constrained by the modeling framework that, apart from the level of epileptogenicity, keeps all the other parameters constant across brain regions. In this way, while capturing some basic biological principles, and especially those that cause interictal discharges or higher excitability of the epileptogenic regions, we are likely to miss some system features. Regional heterogeneity of the parameters of the computational models has already proven to improve the fit of the FC (Deco et al., 2019; Wang et al., 2019). However, without a clear physiological prior to the regional distribution of excitability or some other dynamical or biological parameter, there is a risk to overfit the data. Hence the current approach, which in the future could be improved using parameters obtained by Bayesian inference or other fitting techniques (Deco et al., 2017a; Jirsa et al., 2017; Hashemi et al., 2020), would go beyond epileptogenicity.

Placed in a personalized BNM framework, the novel hybrid model combines generic network dynamic features of neuronal populations at rest and during seizures with patient-specific (compared with BNM informed by generic assumptions or average brains) predictive power of measured network connectivity, that is accessible in human brain imaging. Therefore, despite the loss of direct causality of biophysical mechanisms, the approach advocated here is more explanatory due to its general nature than if it were generated using detailed mechanistic models.

Footnotes

  • The authors declare no competing financial interests.

  • This research was supported by the Fondation pour la Recherche Médicale (Grant DIC20161236442 to V.K.J.), European Union's Horizon 2020 research and innovation programme under grant agreement No. 720270 (SGA1), No. 785907 (SGA2), and No. 945539 (SGA3) Human Brain Project, and the SATT Sud-Est (TVB-Epilepsy). This work has been carried out within the Fédération Hospitalo-Universitaire EPINEXT with the support of the A*MIDEX project (Grant ANR-11-IDEX-0001-02) and the Recherche Hospitalo-Universitaire EPINOV (Grant ANR-17-RHUS-0004) funded by the “Investissements d'Avenir” French Government program managed by the French National Research Agency (ANR).

  • Correspondence should be addressed to Spase Petkoski at spase.petkoski{at}univ-amu.fr or Viktor K. Jirsa at viktor.jirsa{at}univ-amu.fr

This is an open-access article distributed under the terms of the Creative Commons Attribution License Creative Commons Attribution 4.0 International, which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.

References

  1. ↵
    1. Aerts H,
    2. Schirner M,
    3. Dhollander T,
    4. Jeurissen B,
    5. Achten E,
    6. Van Roost D,
    7. Ritter P,
    8. Marinazzo D
    (2020) Modeling brain dynamics after tumor resection using the virtual brain. Neuroimage 213:116738. doi:10.1016/j.neuroimage.2020.116738 pmid:32194282
    OpenUrlCrossRefPubMed
  2. ↵
    1. Allen EA,
    2. Damaraju E,
    3. Plis SM,
    4. Erhardt EB,
    5. Eichele T,
    6. Calhoun VD
    (2014) Tracking whole-brain connectivity dynamics in the resting state. Cereb Cortex 24:663–676. doi:10.1093/cercor/bhs352 pmid:23146964
    OpenUrlCrossRefPubMed
  3. ↵
    1. Arviv O,
    2. Medvedovsky M,
    3. Sheintuch L,
    4. Goldstein A,
    5. Shriki O
    (2016) Deviations from critical dynamics in interictal epileptiform activity. J Neurosci 36:12276–12292. doi:10.1523/JNEUROSCI.0809-16.2016 pmid:27903734
    OpenUrlAbstract/FREE Full Text
  4. ↵
    1. Bandettini PA,
    2. Jesmanowicz A,
    3. Wong EC,
    4. Hyde JS
    (1993) Processing strategies for time-course data sets in functional MRI of the human brain. Magn Reson Med 30:161–173. doi:10.1002/mrm.1910300204 pmid:8366797
    OpenUrlCrossRefPubMed
  5. ↵
    1. Bartolomei F,
    2. Chauvel P,
    3. Wendling F
    (2005) Spatio-temporal dynamics of neuronal networks in partial epilepsy. Rev Neurol (Paris) 161:767–780. doi:10.1016/s0035-3787(05)85136-7 pmid:16244559
    OpenUrlCrossRefPubMed
  6. ↵
    1. Bartolomei F,
    2. Wendling F,
    3. Chauvel P
    (2008) The concept of an epileptogenic network in human partial epilepsies. Neurochirurgie 54:174–184. doi:10.1016/j.neuchi.2008.02.013 pmid:18420227
    OpenUrlCrossRefPubMed
  7. ↵
    1. Bartolomei F,
    2. Guye M,
    3. Wendling F
    (2013a) Abnormal binding and disruption in large scale networks involved in human partial seizures. EPJ Nonlinear Biomed Phys 1:4. doi:10.1140/epjnbp11
    OpenUrlCrossRef
  8. ↵
    1. Bartolomei F,
    2. Bettus G,
    3. Stam C,
    4. Guye M
    (2013b) Interictal network properties in mesial temporal lobe epilepsy: a graph theoretical study from intracerebral recordings. Clin Neurophysiol 124:2345–2353. doi:10.1016/j.clinph.2013.06.003 pmid:23810635
    OpenUrlCrossRefPubMed
  9. ↵
    1. Bartolomei F,
    2. Lagarde S,
    3. Wendling F,
    4. McGonigal A,
    5. Jirsa V,
    6. Guye M,
    7. Bénar C
    (2017) Defining epileptogenic networks: contribution of SEEG and signal analysis. Epilepsia 58:1131–1147. doi:10.1111/epi.13791 pmid:28543030
    OpenUrlCrossRefPubMed
  10. ↵
    1. Besson P,
    2. Bandt SK,
    3. Proix T,
    4. Lagarde S,
    5. Jirsa VK,
    6. Ranjeva J-P,
    7. Bartolomei F,
    8. Guye M
    (2017) Anatomic consistencies across epilepsies: a stereotactic-EEG informed high-resolution structural connectivity study. Brain 140:2639–2652. doi:10.1093/brain/awx181 pmid:28969369
    OpenUrlCrossRefPubMed
  11. ↵
    1. Bonilha L,
    2. Rorden C,
    3. Appenzeller S,
    4. Coan AC,
    5. Cendes F,
    6. Li LM
    (2006) Gray matter atrophy associated with duration of temporal lobe epilepsy. Neuroimage 32:1070–1079. doi:10.1016/j.neuroimage.2006.05.038 pmid:16872843
    OpenUrlCrossRefPubMed
  12. ↵
    1. Bright MG,
    2. Murphy K
    (2015) Is fMRI “noise” really noise? Resting state nuisance regressors remove variance with network structure. Neuroimage 114:158–169. doi:10.1016/j.neuroimage.2015.03.070 pmid:25862264
    OpenUrlCrossRefPubMed
  13. ↵
    1. Beggs JM,
    2. Plenz D
    (2003) Neuronal avalanches in neocortical circuits. J Neurosci 23:11167–11177. pmid:14657176
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. Bettus G,
    2. Wendling F,
    3. Guye M,
    4. Valton L,
    5. Régis J,
    6. Chauvel P,
    7. Bartolomei F
    (2008) Enhanced EEG functional connectivity in mesial temporal lobe epilepsy. Epilepsy Res 81:58–68. doi:10.1016/j.eplepsyres.2008.04.020 pmid:18547787
    OpenUrlCrossRefPubMed
  15. ↵
    1. Bettus G,
    2. Guedj E,
    3. Joyeux F,
    4. Confort-Gouny S,
    5. Soulier E,
    6. Laguitton V,
    7. Cozzone PJ,
    8. Chauvel P,
    9. Ranjeva J-P,
    10. Bartolomei F,
    11. Guye M
    (2009) Decreased basal FMRI functional connectivity in epileptogenic networks and contralateral compensatory mechanisms. Hum Brain Mapp 30:1580–1591. doi:10.1002/hbm.20625 pmid:18661506
    OpenUrlCrossRefPubMed
  16. ↵
    1. Bettus G,
    2. Bartolomei F,
    3. Confort-Gouny S,
    4. Guedj E,
    5. Chauvel P,
    6. Cozzone PJ,
    7. Ranjeva JP,
    8. Guye M
    (2010) Role of resting state functional connectivity MRI in presurgical investigation of mesial temporal lobe epilepsy. J Neurol Neurosurg Psychiatry 81:1147–1154. doi:10.1136/jnnp.2009.191460
    OpenUrlAbstract/FREE Full Text
  17. ↵
    1. Bettus G,
    2. Ranjeva J-P,
    3. Wendling F,
    4. Bénar CG,
    5. Confort-Gouny S,
    6. Régis J,
    7. Chauvel P,
    8. Cozzone PJ,
    9. Lemieux L,
    10. Bartolomei F,
    11. Guye M
    (2011) Interictal functional connectivity of human epileptic networks assessed by intracerebral EEG and bold signal fluctuations. PLoS One 6:e20071. doi:10.1371/journal.pone.0020071 pmid:21625517
    OpenUrlCrossRefPubMed
  18. ↵
    1. Biswal B,
    2. Zerrin Yetkin F,
    3. Haughton VM,
    4. Hyde JS
    (1995) Functional connectivity in the motor cortex of resting human brain using echo-planar MRI. Magn Reson Med 34:537–541. doi:10.1002/mrm.1910340409
    OpenUrlCrossRefPubMed
  19. ↵
    1. Buzsaki G
    (2006) Rhythms of the brain. New York: Oxford UP.
  20. ↵
    1. Cabral J,
    2. Fernandes HM,
    3. Van Hartevelt TJ,
    4. James AC,
    5. Kringelbach ML,
    6. Deco G
    (2013) Structural connectivity in schizophrenia and its impact on the dynamics of spontaneous functional networks. Chaos 23:046111. doi:10.1063/1.4851117 pmid:24387590
    OpenUrlCrossRefPubMed
  21. ↵
    1. Cabral J,
    2. Luckhoo H,
    3. Woolrich M,
    4. Joensson M,
    5. Mohseni H,
    6. Baker A,
    7. Kringelbach ML,
    8. Deco G
    (2014) Exploring mechanisms of spontaneous functional connectivity in MEG: how delayed network interactions lead to structured amplitude envelopes of band-pass filtered oscillations. Neuroimage 90:423–435. doi:10.1016/j.neuroimage.2013.11.047
    OpenUrlCrossRefPubMed
  22. ↵
    1. Chizhov AV,
    2. Zefirov AV,
    3. Amakhin DV,
    4. Smirnova EY,
    5. Zaitsev AV
    (2018) Minimal model of interictal and ictal discharges “epileptor-2”. PLoS Comput Biol 14:e1006186. doi:10.1371/journal.pcbi.1006186
    OpenUrlCrossRefPubMed
  23. ↵
    1. Coito A,
    2. Genetti M,
    3. Pittau F,
    4. Iannotti GR,
    5. Thomschewski A,
    6. Höller Y,
    7. Trinka E,
    8. Wiest R,
    9. Seeck M,
    10. Michel CM,
    11. Plomp G,
    12. Vulliemoz S
    (2016) Altered directed functional connectivity in temporal lobe epilepsy in the absence of interictal spikes: a high density EEG study. Epilepsia 57:402–411. doi:10.1111/epi.13308 pmid:26890734
    OpenUrlCrossRefPubMed
  24. ↵
    1. Cordes D,
    2. Haughton VM,
    3. Arfanakis K,
    4. Carew JD,
    5. Turski PA,
    6. Moritz CH,
    7. Quigley MA,
    8. Meyerand ME
    (2001) Frequencies contributing to functional connectivity in the cerebral cortex in” resting-state” data. Am J Neuroradiol 22:1326–1333.
    OpenUrlAbstract/FREE Full Text
  25. ↵
    1. Costa M,
    2. Goldberger AL,
    3. Peng CK
    (2002) Multiscale entropy analysis of complex physiologic time series. Phys Rev Lett 89:068102. doi:10.1103/PhysRevLett.89.068102 pmid:12190613
    OpenUrlCrossRefPubMed
  26. ↵
    1. Costa M,
    2. Goldberger AL,
    3. Peng CK
    (2005) Multiscale entropy analysis of biological signals. Phys Rev E Stat Nonlin Soft Matter Phys 71:021906. doi:10.1103/PhysRevE.71.021906 pmid:15783351
    OpenUrlCrossRefPubMed
  27. ↵
    1. Courtiol J,
    2. Perdikis D,
    3. Petkoski S,
    4. Müller V,
    5. Huys R,
    6. Sleimen-Malkoun R,
    7. Jirsa VK
    (2016) The multiscale entropy: guidelines for use and interpretation in brain signal analysis. J Neurosci Methods 273:175–190. doi:10.1016/j.jneumeth.2016.09.004 pmid:27639660
    OpenUrlCrossRefPubMed
  28. ↵
    1. da Silva FL,
    2. Pijn J,
    3. Velis D,
    4. Nijssen P
    (1997) Alpha rhythms: noise, dynamics and models. Int J Psychophysiol 26:237–249. doi:10.1016/S0167-8760(97)00767-8 pmid:9203006
    OpenUrlCrossRefPubMed
  29. ↵
    1. Deco G,
    2. Corbetta M
    (2011) The dynamical balance of the brain at rest. Neuroscientist 17:107–123. doi:10.1177/1073858409354384 pmid:21196530
    OpenUrlCrossRefPubMed
  30. ↵
    1. Deco G,
    2. Jirsa VK
    (2012) Ongoing cortical activity at rest: criticality, multistability, and ghost attractors. J Neurosci 32:3366–3375. doi:10.1523/JNEUROSCI.2523-11.2012
    OpenUrlAbstract/FREE Full Text
  31. ↵
    1. Deco G,
    2. Ponce-Alvarez A,
    3. Mantini D,
    4. Romani GL,
    5. Hagmann P,
    6. Corbetta M
    (2013) Resting-state functional connectivity emerges from structurally and dynamically shaped slow linear fluctuations. J Neurosci 33:11239–11252. doi:10.1523/JNEUROSCI.1091-13.2013 pmid:23825427
    OpenUrlAbstract/FREE Full Text
  32. ↵
    1. Deco G,
    2. Kringelbach ML,
    3. Jirsa VK,
    4. Ritter P
    (2017a) The dynamics of resting fluctuations in the brain: metastability and its dynamical cortical core. Sci Rep 7:3095. doi:10.1038/s41598-017-03073-5 pmid:28596608
    OpenUrlCrossRefPubMed
  33. ↵
    1. Deco G,
    2. Tagliazucchi E,
    3. Laufs H,
    4. Sanjuán A,
    5. Kringelbach ML
    (2017b) Novel intrinsic ignition method measuring local-global integration characterizes wakefulness and deep sleep. eNeuro 4:ENEURO.0106-17.2017. doi:10.1523/ENEURO.0106-17.2017
    OpenUrlAbstract/FREE Full Text
  34. ↵
    1. Deco G,
    2. Cruzat J,
    3. Cabral J,
    4. Tagliazucchi E,
    5. Laufs H,
    6. Logothetis NK,
    7. Kringelbach ML
    (2019) Awakening: predicting external stimulation to force transitions between different brain states. Proc Natl Acad Sci U S A 116:18088–18097. doi:10.1073/pnas.1905534116 pmid:31427539
    OpenUrlAbstract/FREE Full Text
  35. ↵
    1. Demirtaş M,
    2. Falcon C,
    3. Tucholka A,
    4. Gispert JD,
    5. Molinuevo JL,
    6. Deco G
    (2017) A whole-brain computational modeling approach to explain the alterations in resting-state functional connectivity during progression of Alzheimer's disease. Neuroimage Clin 16:343–354. doi:10.1016/j.nicl.2017.08.006 pmid:28861336
    OpenUrlCrossRefPubMed
  36. ↵
    1. Desikan RS,
    2. Ségonne F,
    3. Fischl B,
    4. Quinn BT,
    5. Dickerson BC,
    6. Blacker D,
    7. Buckner RL,
    8. Dale AM,
    9. Maguire RP,
    10. Hyman BT,
    11. Albert MS,
    12. Killiany RJ
    (2006) An automated labeling system for subdividing the human cerebral cortex on MRI scans into gyral based regions of interest. Neuroimage 31:968–980. doi:10.1016/j.neuroimage.2006.01.021 pmid:16530430
    OpenUrlCrossRefPubMed
  37. ↵
    1. El Houssaini K,
    2. Ivanov AI,
    3. Bernard C,
    4. Jirsa VK
    (2015) Seizures, refractory status epilepticus, and depolarization block as endogenous brain activities. Phys Rev E Stat Nonlin Soft Matter Phys 91:010701. doi:10.1103/PhysRevE.91.010701 pmid:25679555
    OpenUrlCrossRefPubMed
  38. ↵
    1. El Houssaini K,
    2. Bernard C,
    3. Jirsa VK
    (2020) The epileptor model: a systematic mathematical analysis linked to the dynamics of seizures, refractory status epilepticus and depolarization block. eNeuro 7:ENEURO.0485-18.2019. doi:10.1523/ENEURO.0485-18.2019
    OpenUrlAbstract/FREE Full Text
  39. ↵
    1. Falcon MI,
    2. Riley JD,
    3. Jirsa VK,
    4. McIntosh AR,
    5. Shereen AD,
    6. Chen EE,
    7. Solodkin A
    (2015) The virtual brain: modeling biological correlates of recovery after chronic stroke. Front Neurol 6:228. doi:10.3389/fneur.2015.00228 pmid:26579071
    OpenUrlCrossRefPubMed
  40. ↵
    1. Fischl B
    (2012) Freesurfer. Neuroimage 62:774–781. doi:10.1016/j.neuroimage.2012.01.021 pmid:22248573
    OpenUrlCrossRefPubMed
  41. ↵
    1. Fox MD,
    2. Snyder AZ,
    3. Vincent JL,
    4. Corbetta M,
    5. Van Essen DC,
    6. Raichle ME
    (2005) The human brain is intrinsically organized into dynamic, anticorrelated functional networks. Proc Natl Acad Sci U S A 102:9673–9678. doi:10.1073/pnas.0504136102 pmid:15976020
    OpenUrlAbstract/FREE Full Text
  42. ↵
    1. Fox MD,
    2. Zhang D,
    3. Snyder AZ,
    4. Raichle ME
    (2009) The global signal and observed anticorrelated resting state brain networks. J Neurophysiol 101:3270–3283. doi:10.1152/jn.90777.2008 pmid:19339462
    OpenUrlCrossRefPubMed
  43. ↵
    1. Freyer F,
    2. Roberts JA,
    3. Becker R,
    4. Robinson PA,
    5. Ritter P,
    6. Breakspear M
    (2011) Biophysical mechanisms of multistability in resting-state cortical rhythms. J Neurosci 31:6353–6361. doi:10.1523/JNEUROSCI.6693-10.2011 pmid:21525275
    OpenUrlAbstract/FREE Full Text
  44. ↵
    1. Freyer F,
    2. Roberts JA,
    3. Ritter P,
    4. Breakspear M
    (2012) A canonical model of multistability and scale-invariance in biological systems. PLoS Comput Biol 8:e1002634. doi:10.1371/journal.pcbi.1002634
    OpenUrlCrossRefPubMed
  45. ↵
    1. Friston K,
    2. Harrison L,
    3. Penny W
    (2003) Dynamic causal modelling. Neuroimage 19:1273–1302. doi:10.1016/S1053-8119(03)00202-7
    OpenUrlCrossRefPubMed
  46. ↵
    1. Ghosh A,
    2. Rho Y,
    3. McIntosh AR,
    4. Kötter R,
    5. Jirsa VK
    (2008) Noise during rest enables the exploration of the brain's dynamic repertoire. PLoS Comput Biol 4:e1000196. doi:10.1371/journal.pcbi.1000196
    OpenUrlCrossRefPubMed
  47. ↵
    1. Glerean E,
    2. Salmi J,
    3. Lahnakoski JM,
    4. Jääskeläinen IP,
    5. Sams M
    (2012) Functional magnetic resonance imaging phase synchronization as a measure of dynamic functional connectivity. Brain Connect 2:91–101. doi:10.1089/brain.2011.0068 pmid:22559794
    OpenUrlCrossRefPubMed
  48. ↵
    1. Guye M,
    2. Bartolomei F,
    3. Ranjeva JP
    (2008) Imaging structural and functional connectivity: towards a unified definition of human brain organization? Curr Opin Neurol 21:393–403. doi:10.1097/WCO.0b013e3283065cfb pmid:18607198
    OpenUrlCrossRefPubMed
  49. ↵
    1. Haimovici A,
    2. Tagliazucchi E,
    3. Balenzuela P,
    4. Laufs H
    (2017) On wakefulness fluctuations as a source of bold functional connectivity dynamics. Sci Rep 7:5908. doi:10.1038/s41598-017-06389-4 pmid:28724928
    OpenUrlCrossRefPubMed
  50. ↵
    1. Hansen EC,
    2. Battaglia D,
    3. Spiegler A,
    4. Deco G,
    5. Jirsa VK
    (2015) Functional connectivity dynamics: modeling the switching behavior of the resting state. Neuroimage 105:525–535. doi:10.1016/j.neuroimage.2014.11.001 pmid:25462790
    OpenUrlCrossRefPubMed
  51. ↵
    1. Hashemi M,
    2. Vattikonda AN,
    3. Sip V,
    4. Guye M,
    5. Bartolomei F
    6. Woodman MM
    7. Jirsa VK
    (2020) The Bayesian virtual epileptic patient: a probabilistic framework designed to infer the spatial map of epileptogenicity in a personalized large-scale brain model of epilepsy spread. NeuroImage 217:116839.
    OpenUrl
  52. ↵
    1. Hutchison RM,
    2. Womelsdorf T,
    3. Allen EA,
    4. Bandettini PA,
    5. Calhoun VD,
    6. Corbetta M,
    7. Penna SD,
    8. Duyn JH,
    9. Glover GH,
    10. Gonzalez-Castillo J,
    11. Handwerker DA,
    12. Keilholz S,
    13. Kiviniemi V,
    14. Leopold DA,
    15. de Pasquale F,
    16. Sporns O,
    17. Walter M,
    18. Chang C
    (2013) Dynamic functional connectivity: promise, issues, and interpretations. Neuroimage 80:360–378. doi:10.1016/j.neuroimage.2013.05.079 pmid:23707587
    OpenUrlCrossRefPubMed
  53. ↵
    1. Ikeda A,
    2. Taki W,
    3. Kunieda T,
    4. Terada K,
    5. Mikuni N,
    6. Nagamine T,
    7. Yazawa S,
    8. Ohara S,
    9. Hori T,
    10. Kaji R,
    11. Kimura J,
    12. Shibasaki H
    (1999) Focal ictal direct current shifts in humanepilepsy as studied by subdural and scalp recording. Brain 122:827–838. doi:10.1093/brain/122.5.827
    OpenUrlCrossRefPubMed
  54. ↵
    1. Izhikevich EM
    (2007) Dynamical systems in neuroscience. Cambridge, MA: MIT.
  55. ↵
    1. Jenkinson M,
    2. Beckmann CF,
    3. Behrens TE,
    4. Woolrich MW,
    5. Smith SM
    (2012) Fsl. Neuroimage 62:782–790. doi:10.1016/j.neuroimage.2011.09.015 pmid:21979382
    OpenUrlCrossRefPubMed
  56. ↵
    1. Jirsa VK,
    2. Stacey WC,
    3. Quilichini PP,
    4. Ivanov AI,
    5. Bernard C
    (2014) On the nature of seizure dynamics. Brain 137:2210–2230. doi:10.1093/brain/awu133 pmid:24919973
    OpenUrlCrossRefPubMed
  57. ↵
    1. Jirsa VK,
    2. Proix T,
    3. Perdikis D,
    4. Woodman MM,
    5. Wang H,
    6. Gonzalez-Martinez J,
    7. Bernard C,
    8. Bénar C,
    9. Guye M,
    10. Chauvel P,
    11. Bartolomei F
    (2017) The virtual epileptic patient: individualized whole-brain models of epilepsy spread. Neuroimage 145:377–388. doi:10.1016/j.neuroimage.2016.04.049 pmid:27477535
    OpenUrlCrossRefPubMed
  58. ↵
    1. Jobst BM,
    2. Hindriks R,
    3. Laufs H,
    4. Tagliazucchi E,
    5. Hahn G,
    6. Ponce-Alvarez A,
    7. Stevner ABA,
    8. Kringelbach ML,
    9. Deco G
    (2017) Increased stability and breakdown of brain effective connectivity during slow-wave sleep: mechanistic insights from whole-rain computational modelling. Sci Rep 7:4634. doi:10.1038/s41598-017-04522-x pmid:28680119
    OpenUrlCrossRefPubMed
  59. ↵
    1. Kehoe PD,
    2. Saenger VM,
    3. Lisofsky N,
    4. Kühn S,
    5. Kringelbach ML,
    6. Schwarzbach J,
    7. Deco G
    (2017) Consistent local dynamics in the brain across sessions are revealed by whole brain modeling of resting state activity. bioRxiv. Advance online publication. Retrieved June 9, 2020. doi:10.1101/104232
    OpenUrlCrossRef
  60. ↵
    1. Kuramoto Y
    (1984) Chemical oscillations, waves and turbulence. Berlin: Springer.
  61. ↵
    1. Kuznetsov YA
    (2013) Elements of applied bifurcation theory. New York: Springer Science and Business Media.
  62. ↵
    1. Lagarde S,
    2. Roehri N,
    3. Lambert I,
    4. Trebuchon A,
    5. McGonigal A,
    6. Carron R,
    7. Scavarda D,
    8. Milh M,
    9. Pizzo F,
    10. Colombet B,
    11. Giusiano B,
    12. Medina Villalon S,
    13. Guye M,
    14. Bénar C-G,
    15. Bartolomei F
    (2018) Interictal stereotactic-EEG functional connectivity in refractory focal epilepsies. Brain 141:2966–2980. doi:10.1093/brain/awy214 pmid:30107499
    OpenUrlCrossRefPubMed
  63. ↵
    1. Laufs H,
    2. Rodionov R,
    3. Thornton R,
    4. Duncan JS,
    5. Lemieux L,
    6. Tagliazucchi E
    (2014) Altered FMRI connectivity dynamics in temporal lobe epilepsy might explain seizure semiology. Front Neurol 5:175. doi:10.3389/fneur.2014.00175 pmid:25309503
    OpenUrlCrossRefPubMed
  64. ↵
    1. Laumann TO,
    2. Snyder AZ,
    3. MA,
    4. Gordon EM,
    5. Gratton C,
    6. Adeyemo B,
    7. Gilmore AW,
    8. Nelson SM,
    9. Berg JJ,
    10. Greene DJ,
    11. McCarthy JE,
    12. Tagliazucchi E,
    13. Laufs H,
    14. Schlaggar BL,
    15. Dosenbach NUF,
    16. Petersen SE
    (2017) On the stability of BOLD fMRI correlations. Cereb Cortex 27:4719–4732.
    OpenUrlCrossRefPubMed
  65. ↵
    1. Lee WH,
    2. Frangou S
    (2017) Linking functional connectivity and dynamic properties of resting-state networks. Sci Rep 7:16610. doi:10.1038/s41598-017-16789-1 pmid:29192157
    OpenUrlCrossRefPubMed
  66. ↵
    1. Lemieux L,
    2. Salek-Haddadi A,
    3. Lund TE,
    4. Laufs H,
    5. Carmichael D
    (2007) Modelling large motion events in fMRI studies of patients with epilepsy. Magn Reson Imaging 25:894–901. doi:10.1016/j.mri.2007.03.009 pmid:17490845
    OpenUrlCrossRefPubMed
  67. ↵
    1. Li Hegner Y,
    2. Marquetand J,
    3. Elshahabi A,
    4. Klamer S,
    5. Lerche H,
    6. Braun C,
    7. Focke NK
    (2018) Increased functional MEG connectivity as a hallmark of MRI-negative focal and generalized epilepsy. Brain Topogr 31:863–874. doi:10.1007/s10548-018-0649-4 pmid:29766384
    OpenUrlCrossRefPubMed
  68. ↵
    1. Liao W,
    2. Zhang Z,
    3. Pan Z,
    4. Mantini D,
    5. Ding J,
    6. Duan X,
    7. Luo C,
    8. Lu G,
    9. Chen H
    (2010) Altered functional connectivity and small-world in mesial temporal lobe epilepsy. PLoS One 5:e8525. doi:10.1371/journal.pone.0008525 pmid:20072616
    OpenUrlCrossRefPubMed
  69. ↵
    1. Lüders HO,
    2. Najm I,
    3. Nair D,
    4. Widdess-Walsh P,
    5. Bingman W
    (2006) The epileptogenic zone: general principles. Epileptic Disord 8:S1–S9.
    OpenUrlPubMed
  70. ↵
    1. Mannella R
    (2002) Integration of stochastic differential equations on a computer. Int J Mod Phys C 13:1177–1194. doi:10.1142/S0129183102004042
    OpenUrlCrossRef
  71. ↵
    1. Monto S,
    2. Vanhatalo S,
    3. Holmes MD,
    4. Palva JM
    (2007) Epileptogenic neocortical networks are revealed by abnormal temporal dynamics in seizure-free subdural EEG. Cereb Cortex 17:1386–1393. doi:10.1093/cercor/bhl049 pmid:16908492
    OpenUrlCrossRefPubMed
  72. ↵
    1. Murphy K,
    2. Birn RM,
    3. Handwerker DA,
    4. Jones TB,
    5. Bandettini PA
    (2009) The impact of global signal regression on resting state correlations: are anti-correlated networks introduced? Neuroimage 44:893–905. doi:10.1016/j.neuroimage.2008.09.036 pmid:18976716
    OpenUrlCrossRefPubMed
  73. ↵
    1. Nakagawa TT,
    2. Jirsa VK,
    3. Spiegler A,
    4. McIntosh AR,
    5. Deco G
    (2013) Bottom up modeling of the connectome: linking structure and function in the resting brain and their changes in aging. Neuroimage 80:318–329. doi:10.1016/j.neuroimage.2013.04.055 pmid:23629050
    OpenUrlCrossRefPubMed
  74. ↵
    1. Nedic S,
    2. Stufflebeam SM,
    3. Rondinoni C,
    4. Velasco TR,
    5. dos Santos AC,
    6. Leite JP,
    7. Gargaro AC,
    8. Mujica-Parodi LR,
    9. Ide JS
    (2015) Using network dynamic FMRI for detection of epileptogenic foci. BMC Neurol 15:262. doi:10.1186/s12883-015-0514-y pmid:26689596
    OpenUrlCrossRefPubMed
  75. ↵
    1. Parish L,
    2. Worrell G,
    3. Cranstoun S,
    4. Stead S,
    5. Pennell P,
    6. Litt B
    (2004) Long-range temporal correlations in epileptogenic and non-epileptogenic human hippocampus. Neuroscience 125:1069–1076. doi:10.1016/j.neuroscience.2004.03.002 pmid:15120866
    OpenUrlCrossRefPubMed
  76. ↵
    1. Pedersen M,
    2. Omidvarnia A,
    3. Curwood EK,
    4. Walz JM,
    5. Rayner G,
    6. Jackson GD
    (2017) The dynamics of functional connectivity in neocortical focal epilepsy. Neuroimage Clin 15:209–214. doi:10.1016/j.nicl.2017.04.005 pmid:28529877
    OpenUrlCrossRefPubMed
  77. ↵
    1. Pereira FRS,
    2. Alessio A,
    3. Sercheli MS,
    4. Pedro T,
    5. Bilevicius E,
    6. Rondina JM,
    7. Ozelo HFB,
    8. Castellano G,
    9. Covolan RJM,
    10. Damasceno BP,
    11. Cendes F
    (2010) Asymmetrical hippocampal connectivity in mesial temporal lobe epilepsy: evidence from resting state FMRI. BMC Neurosci 11:66. doi:10.1186/1471-2202-11-66 pmid:20525202
    OpenUrlCrossRefPubMed
  78. ↵
    1. Petkoski S,
    2. Jirsa VK
    (2019) Transmission time delays organize the brain network synchronization. Philos Trans A Math Phys Eng Sci 377:20180132. doi:10.1098/rsta.2018.0132 pmid:31329065
    OpenUrlCrossRefPubMed
  79. ↵
    1. Petkoski S,
    2. Spiegler A,
    3. Proix T,
    4. Aram P,
    5. Temprado JJ,
    6. Jirsa VK
    (2016) Heterogeneity of time delays determines synchronization of coupled oscillators. Phys Rev E 94:012209. doi:10.1103/PhysRevE.94.012209 pmid:27575125
    OpenUrlCrossRefPubMed
  80. ↵
    1. Petkoski S,
    2. Palva JM,
    3. Jirsa VK
    (2018) Phase-lags in large scale brain synchronization: methodological considerations and in-silico analysis. PLoS Comput Biol 14:e1006160. doi:10.1371/journal.pcbi.1006160 pmid:29990339
    OpenUrlCrossRefPubMed
  81. ↵
    1. Pittau F,
    2. Grova C,
    3. Moeller F,
    4. Dubeau F,
    5. Gotman J
    (2012) Patterns of altered functional connectivity in mesial temporal lobe epilepsy. Epilepsia 53:1013–1023. doi:10.1111/j.1528-1167.2012.03464.x pmid:22578020
    OpenUrlCrossRefPubMed
  82. ↵
    1. Ponce-Alvarez A,
    2. Deco G,
    3. Hagmann P,
    4. Romani GL,
    5. Mantini D,
    6. Corbetta M
    (2015) Resting-state temporal synchronization networks emerge from connectivity topology and heterogeneity. PLoS Comput Biol 11:e1004100. doi:10.1371/journal.pcbi.1004100
    OpenUrlCrossRefPubMed
  83. ↵
    1. Power JD,
    2. Barnes KA,
    3. Snyder AZ,
    4. Schlaggar BL,
    5. Petersen SE
    (2012) Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. Neuroimage 59:2142–2154. doi:10.1016/j.neuroimage.2011.10.018 pmid:22019881
    OpenUrlCrossRefPubMed
  84. ↵
    1. Proix T,
    2. Bartolomei F,
    3. Chauvel P,
    4. Bernard C,
    5. Jirsa VK
    (2014) Permittivity coupling across brain regions determines seizure recruitment in partial epilepsy. J Neurosci 34:15009–15021. doi:10.1523/JNEUROSCI.1570-14.2014 pmid:25378166
    OpenUrlAbstract/FREE Full Text
  85. ↵
    1. Proix T,
    2. Spiegler A,
    3. Schirner M,
    4. Rothmeier S,
    5. Ritter P,
    6. Jirsa VK
    (2016) How do parcellation size and short-range connectivity affect dynamics in large-scale brain network models? Neuroimage 142:135–149. doi:10.1016/j.neuroimage.2016.06.016 pmid:27480624
    OpenUrlCrossRefPubMed
  86. ↵
    1. Proix T,
    2. Bartolomei F,
    3. Guye M,
    4. Jirsa VK
    (2017) Individual brain structure and modelling predict seizure propagation. Brain 140:641–654. doi:10.1093/brain/awx004 pmid:28364550
    OpenUrlCrossRefPubMed
  87. ↵
    1. Richman JS,
    2. Moorman JR
    (2000) Physiological time-series analysis using approximate entropy and sample entropy. Am J Physiol Heart Circ Physiol 278:H2039–H2049. doi:10.1152/ajpheart.2000.278.6.H2039 pmid:10843903
    OpenUrlCrossRefPubMed
  88. ↵
    1. Ridley B,
    2. Wirsich J,
    3. Bettus G,
    4. Rodionov R,
    5. Murta T,
    6. Chaudhary U,
    7. Carmichael D,
    8. Thornton R,
    9. Vulliemoz S,
    10. McEvoy A,
    11. Wendling F,
    12. Bartolomei F,
    13. Ranjeva J-P,
    14. Lemieux L,
    15. Guye M
    (2017) Simultaneous intracranial EEG-FMRI shows inter-modality correlation in time-resolved connectivity within normal areas but not within epileptic regions. Brain Topogr 30:639–655. doi:10.1007/s10548-017-0551-5 pmid:28194612
    OpenUrlCrossRefPubMed
  89. ↵
    1. Ridley BGY,
    2. Rousseau C,
    3. Wirsich J,
    4. Le Troter A,
    5. Soulier E,
    6. Confort-Gouny S,
    7. Bartolomei F,
    8. Ranjeva JP,
    9. Achard S,
    10. Guye M
    (2015) Nodal approach reveals differential impact of lateralized focal epilepsies on hub reorganization. Neuroimage 118:39–48. doi:10.1016/j.neuroimage.2015.05.096 pmid:26070261
    OpenUrlCrossRefPubMed
  90. ↵
    1. Sanz-Leon P,
    2. Knock SA,
    3. Woodman MM,
    4. Domide L,
    5. Mersmann J,
    6. McIntosh AR,
    7. Jirsa V
    (2013) The virtual brain: a simulator of primate brain network dynamics. Front Neuroinform 7:10. doi:10.3389/fninf.2013.00010 pmid:23781198
    OpenUrlCrossRefPubMed
  91. ↵
    1. Sanz-Leon P,
    2. Knock SA,
    3. Spiegler A,
    4. Jirsa VK
    (2015) Mathematical framework for large-scale brain network modeling in the virtual brain. Neuroimage 111:385–430. doi:10.1016/j.neuroimage.2015.01.002 pmid:25592995
    OpenUrlCrossRefPubMed
  92. ↵
    1. Scharfman HE
    (2007) The neurobiology of epilepsy. Curr Neurol Neurosci Rep 7:348–354. doi:10.1007/s11910-007-0053-z pmid:17618543
    OpenUrlCrossRefPubMed
  93. ↵
    1. Schevon CA,
    2. Cappell J,
    3. Emerson R,
    4. Isler J,
    5. Grieve P,
    6. Goodman R,
    7. Mckhann G,
    8. Weiner H,
    9. Doyle W,
    10. Kuzniecky R,
    11. Devinsky O,
    12. Gilliam F
    (2007) Cortical abnormalities in epilepsy revealed by local EEG synchrony. Neuroimage 35:140–148. doi:10.1016/j.neuroimage.2006.11.009
    OpenUrlCrossRefPubMed
  94. ↵
    1. Schirner M,
    2. Rothmeier S,
    3. Jirsa VK,
    4. McIntosh AR,
    5. Ritter P
    (2015) An automated pipeline for constructing personalized virtual brains from multimodal neuroimaging data. Neuroimage 117:343–357. doi:10.1016/j.neuroimage.2015.03.055 pmid:25837600
    OpenUrlCrossRefPubMed
  95. ↵
    1. Schirner M,
    2. McIntosh AR,
    3. Jirsa V,
    4. Deco G,
    5. Ritter P
    (2018) Inferring multi-scale neural mechanisms with brain network modelling. eLife 7:e28927. doi:10.7554/eLife.28927
    OpenUrlCrossRef
  96. ↵
    1. Shanahan M
    (2010) Metastable chimera states in community-structured oscillator networks. Chaos 20:013108. doi:10.1063/1.3305451 pmid:20370263
    OpenUrlCrossRefPubMed
  97. ↵
    1. Smitha K,
    2. Akhil Raja K,
    3. Arun K,
    4. Rajesh P,
    5. Thomas B,
    6. Kapilamoorthy T,
    7. Kesavadas C
    (2017) Resting state FMRI: a review on methods in resting state connectivity analysis and resting state networks. Neuroradiol J 30:305–317. doi:10.1177/1971400917697342 pmid:28353416
    OpenUrlCrossRefPubMed
  98. ↵
    1. Spiegler A,
    2. Hansen EC,
    3. Bernard C,
    4. McIntosh AR,
    5. Jirsa VK
    (2016) Selective activation of resting state networks following focal stimulation in a connectome-based network model of the human brain. eNeuro 3:ENEURO.0068-16.2016. doi:10.1523/ENEURO.0068-16.2016
    OpenUrlAbstract/FREE Full Text
  99. ↵
    1. Stefanescu RA,
    2. Jirsa VK
    (2008) A low dimensional description of globally coupled heterogeneous neural networks of excitatory and inhibitory neurons. PLoS Comput Biol 4:e1000219. doi:10.1371/journal.pcbi.1000219 pmid:19008942
    OpenUrlCrossRefPubMed
  100. ↵
    1. Su L,
    2. An J,
    3. Ma Q,
    4. Qiu S,
    5. Hu D
    (2015) Influence of resting-state network on lateralization of functional connectivity in mesial temporal lobe epilepsy. AJNR Am J Neuroradiol 36:1479–1487. doi:10.3174/ajnr.A4346 pmid:26021622
    OpenUrlAbstract/FREE Full Text
  101. ↵
    1. Tagliazucchi E,
    2. Laufs H
    (2014) Decoding wakefulness levels from typical FMRI resting-state data reveals reliable drifts between wakefulness and sleep. Neuron 82:695–708. doi:10.1016/j.neuron.2014.03.020 pmid:24811386
    OpenUrlCrossRefPubMed
  102. ↵
    1. Tavares V,
    2. Ribeiro AS,
    3. Capela C,
    4. Cerqueira L,
    5. Ferreira HA
    (2015) Epileptogenic focus localization: a new approach. EJNMMI Phys 2 [Suppl 1]:A81. doi:10.1186/2197-7364-2-S1-A81 pmid:26956343
    OpenUrlCrossRefPubMed
  103. ↵
    1. Tournier JD,
    2. Calamante F,
    3. Connelly A
    (2007) Robust determination of the fibre orientation distribution in diffusion MRI: non-negativity constrained super-resolved spherical deconvolution. Neuroimage 35:1459–1472. doi:10.1016/j.neuroimage.2007.02.016 pmid:17379540
    OpenUrlCrossRefPubMed
  104. ↵
    1. Tournier JD,
    2. Calamante F,
    3. Connelly A
    (2010) Improved probabilistic streamlines tractography by 2nd order integration over fibre orientation distributions. Paper presented at 2010 ISMRM, Stockholm, Sweden, May.
  105. ↵
    1. Tournier JD,
    2. Calamante F,
    3. Connelly A
    (2012) Mrtrix: diffusion tractography in crossing fiber regions. Int J Imaging Syst Technol 22:53–66. doi:10.1002/ima.22005
    OpenUrlCrossRef
  106. ↵
    1. Tournier JD,
    2. Calamante F,
    3. Connelly A
    (2013) Determination of the appropriate b value and number of gradient directions for high-angular-resolution diffusion-weighted imaging. NMR Biomed 26:1775–1786. doi:10.1002/nbm.3017 pmid:24038308
    OpenUrlCrossRefPubMed
  107. ↵
    1. Tracy JI,
    2. Doucet GE
    (2015) Resting-state functional connectivity in epilepsy: growing relevance for clinical decision making. Curr Opin Neurol 28:158–165. doi:10.1097/WCO.0000000000000178 pmid:25734954
    OpenUrlCrossRefPubMed
  108. ↵
    1. van Houdt PJ,
    2. Ossenblok PP,
    3. Colon AJ,
    4. Hermans KH,
    5. Verdaasdonk RM,
    6. Boon PA,
    7. de Munck JC
    (2015) Are epilepsy-related FMRI components dependent on the presence of interictal epileptic discharges in scalp EEG? Brain Topogr 28:606–618. doi:10.1007/s10548-014-0407-1 pmid:25315607
    OpenUrlCrossRefPubMed
  109. ↵
    1. Vanhatalo S,
    2. Holmes M,
    3. Tallgren P,
    4. Voipio J,
    5. Kaila K,
    6. Miller J
    (2003) Very slow EEG responses lateralize temporal lobe seizures an evaluation of non-invasive DC-EEG. Neurology 60:1098–1104. doi:10.1212/01.wnl.0000052993.37621.cc pmid:12682313
    OpenUrlAbstract/FREE Full Text
  110. ↵
    1. Voets NL,
    2. Zamboni G,
    3. Stokes MG,
    4. Carpenter K,
    5. Stacey R,
    6. Adcock JE
    (2014) Aberrant functional connectivity in dissociable hippocampal networks is associated with deficits in memory. J Neurosci 34:4920–4928. doi:10.1523/JNEUROSCI.4281-13.2014 pmid:24695711
    OpenUrlAbstract/FREE Full Text
  111. ↵
    1. Waites AB,
    2. Briellmann RS,
    3. Saling MM,
    4. Abbott DF,
    5. Jackson GD
    (2006) Functional connectivity networks are disrupted in left temporal lobe epilepsy. Ann Neurol 59:335–343. doi:10.1002/ana.20733 pmid:16404743
    OpenUrlCrossRefPubMed
  112. ↵
    1. Wandschneider B,
    2. Stretton J,
    3. Sidhu M,
    4. Centeno M,
    5. Kozák LR,
    6. Symms M,
    7. Thompson PJ,
    8. Duncan JS,
    9. Koepp MJ
    (2014) Levetiracetam reduces abnormal network activations in temporal lobe epilepsy. Neurology 83:1508–1512. doi:10.1212/WNL.0000000000000910
    OpenUrlCrossRefPubMed
  113. ↵
    1. Wang P,
    2. Kong R,
    3. Kong X,
    4. Liégeois R,
    5. Orban C,
    6. Deco G,
    7. van den Heuvel MP,
    8. Thomas Yeo B
    (2019) Inversion of a large-scale circuit model reveals a cortical hierarchy in the dynamic resting human brain. Sci Adv 5:eaat7854. doi:10.1126/sciadv.aat7854
    OpenUrlFREE Full Text
  114. ↵
    1. Wildie M,
    2. Shanahan M
    (2012) Metastability and chimera states in modular delay and pulse-coupled oscillator networks. Chaos 22:043131. doi:10.1063/1.4766592 pmid:23278066
    OpenUrlCrossRefPubMed
  115. ↵
    1. Wirsich J,
    2. Perry A,
    3. Ridley B,
    4. Proix T,
    5. Golos M,
    6. Bénar C,
    7. Ranjeva J-P,
    8. Bartolomei F,
    9. Breakspear M,
    10. Jirsa V,
    11. Guye M
    (2016) Whole-brain analytic measures of network communication reveal increased structure-function correlation in right temporal lobe epilepsy. Neuroimage Clin 11:707–718. doi:10.1016/j.nicl.2016.05.010 pmid:27330970
    OpenUrlCrossRefPubMed
  116. ↵
    1. Woolrich MW,
    2. Ripley BD,
    3. Brady M,
    4. Smith SM
    (2001) Temporal autocorrelation in univariate linear modeling of FMRI data. Neuroimage 14:1370–1386. doi:10.1006/nimg.2001.0931 pmid:11707093
    OpenUrlCrossRefPubMed
  117. ↵
    1. Yang AC,
    2. Huang C-C,
    3. Yeh H-L,
    4. Liu M-E,
    5. Hong C-J,
    6. Tu P-C,
    7. Chen J-F,
    8. Huang NE,
    9. Peng C-K,
    10. Lin C-P,
    11. Tsai S-J
    (2013) Complexity of spontaneous bold activity in default mode network is correlated with cognitive function in normal male elderly: a multiscale entropy analysis. Neurobiol Aging 34:428–438. doi:10.1016/j.neurobiolaging.2012.05.004 pmid:22683008
    OpenUrlCrossRefPubMed
Back to top

In this issue

The Journal of Neuroscience: 40 (29)
Journal of Neuroscience
Vol. 40, Issue 29
15 Jul 2020
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Advertising (PDF)
  • Ed Board (PDF)
Email

Thank you for sharing this Journal of Neuroscience article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Dynamical Mechanisms of Interictal Resting-State Functional Connectivity in Epilepsy
(Your Name) has forwarded a page to you from Journal of Neuroscience
(Your Name) thought you would be interested in this article in Journal of Neuroscience.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Print
View Full Page PDF
Citation Tools
Dynamical Mechanisms of Interictal Resting-State Functional Connectivity in Epilepsy
Julie Courtiol, Maxime Guye, Fabrice Bartolomei, Spase Petkoski, Viktor K. Jirsa
Journal of Neuroscience 15 July 2020, 40 (29) 5572-5588; DOI: 10.1523/JNEUROSCI.0905-19.2020

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Respond to this article
Request Permissions
Share
Dynamical Mechanisms of Interictal Resting-State Functional Connectivity in Epilepsy
Julie Courtiol, Maxime Guye, Fabrice Bartolomei, Spase Petkoski, Viktor K. Jirsa
Journal of Neuroscience 15 July 2020, 40 (29) 5572-5588; DOI: 10.1523/JNEUROSCI.0905-19.2020
Reddit logo Twitter logo Facebook logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Introduction
    • Materials and Methods
    • Results
    • Discussion
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

Keywords

  • brain dynamics
  • brain network model
  • epilepsy
  • fMRI
  • functional connectivity
  • resting state

Responses to this article

Respond to this article

Jump to comment:

No eLetters have been published for this article.

Related Articles

Cited By...

More in this TOC Section

Research Articles

  • Chronic desipramine reverses deficits in cell activity, norepinephrine innervation and anxiety-depression phenotypes in fluoxetine-resistant cF1ko mice
  • Dissociative effects of age on neural differentiation at the category and item level
  • Featural representation and internal noise underlie the eccentricity effect in contrast sensitivity
Show more Research Articles

Systems/Circuits

  • Cortical Synchrony and Information Flow during Transition from Wakefulness to Light Non-Rapid Eye Movement Sleep
  • Hypersensitivity to Distractors in Fragile X Syndrome from Loss of Modulation of Cortical VIP Interneurons
  • The Structure of Hippocampal CA1 Interactions Optimizes Spatial Coding across Experience
Show more Systems/Circuits
  • Home
  • Alerts
  • Visit Society for Neuroscience on Facebook
  • Follow Society for Neuroscience on Twitter
  • Follow Society for Neuroscience on LinkedIn
  • Visit Society for Neuroscience on Youtube
  • Follow our RSS feeds

Content

  • Early Release
  • Current Issue
  • Issue Archive
  • Collections

Information

  • For Authors
  • For Advertisers
  • For the Media
  • For Subscribers

About

  • About the Journal
  • Editorial Board
  • Privacy Policy
  • Contact
(JNeurosci logo)
(SfN logo)

Copyright © 2023 by the Society for Neuroscience.
JNeurosci Online ISSN: 1529-2401

The ideas and opinions expressed in JNeurosci do not necessarily reflect those of SfN or the JNeurosci Editorial Board. Publication of an advertisement or other product mention in JNeurosci should not be construed as an endorsement of the manufacturer’s claims. SfN does not assume any responsibility for any injury and/or damage to persons or property arising from or related to any use of any material contained in JNeurosci.