Epilepsy surgery is the treatment of choice for drug-resistant epilepsy patients, but only leads to seizure freedom for roughly two in three patients. To address this problem, we designed a patient-specific epilepsy surgery model combining large-scale magnetoencephalography (MEG) brain networks with an epidemic spreading model. This simple model was enough to reproduce the stereo-tactical electroencephalography (SEEG) seizure propagation patterns of all patients (N = 15), when considering the resection areas (RA) as the epidemic seed. Moreover, the goodness of fit of the model predicted surgical outcome. Once adapted for each patient, the model can generate alternative hypothesis of the seizure onset zone and test different resection strategies in silico. Overall, our findings indicate that spreading models based on patient-specific MEG connectivity can be used to predict surgical outcomes, with better fit results and greater reduction on seizure propagation linked to higher likelihood of seizure freedom after surgery. Finally, we introduced a population model that can be individualized by considering only the patient-specific MEG network, and showed that it not only conserves but improves the group classification. Thus, it may pave the way to generalize this framework to patients without SEEG recordings, reduce the risk of overfitting and improve the stability of the analyses.

How temporal modulations in functional interactions are shaped by the underlying anatomical connections remains an open question. Here, we analyse the role of structural eigenmodes, in the formation and dissolution of temporally evolving functional brain networks using resting-state magnetoencephalography and diffusion magnetic resonance imaging data at the individual subject level. Our results show that even at short timescales, phase and amplitude connectivity can partly be expressed by structural eigenmodes, but hardly by direct structural connections. Albeit a stronger relationship was found between structural eigenmodes and time-resolved amplitude connectivity. Time-resolved connectivity for both phase and amplitude was mostly characterised by a stationary process, superimposed with very brief periods that showed deviations from this stationary process. For these brief periods, dynamic network states were extracted that showed different expressions of eigenmodes. Furthermore, the eigenmode expression was related to overall cognitive performance and co-occurred with fluctuations in community structure of functional networks. These results implicate that ongoing time-resolved resting-state networks, even at short timescales, can to some extent be understood in terms of activation and deactivation of structural eigenmodes and that these eigenmodes play a role in the dynamic integration and segregation of information across the cortex, subserving cognitive functions.

Large-scale neurophysiological networks are often reconstructed from band-pass filtered time series derived from magnetoencephalography (MEG) data. Common practice is to reconstruct these networks separately for different frequency bands and to treat them independently. Recent evidence suggests that this separation may be inadequate, as there can be significant coupling between frequency bands (interlayer connectivity). A multilayer network approach offers a solution to analyze frequency-specific networks in one framework. We propose to use a recently developed network reconstruction method in conjunction with phase oscillator models to estimate interlayer connectivity that optimally fits the empirical data. This approach determines interlayer connectivity based on observed frequency-specific time series of the phase and a connectome derived from diffusion weighted imaging. The performance of this interlayer reconstruction method was evaluated in-silico. Our reconstruction of the underlying interlayer connectivity agreed to very high degree with the ground truth. Subsequently, we applied our method to empirical resting-state MEG data obtained from healthy subjects and reconstructed two-layered networks consisting of either alpha-to-beta or theta-to-gamma band connectivity. Our analysis revealed that interlayer connectivity is dominated by a multiplex structure, i.e. by one-to-one interactions for both alpha-to-beta band and theta-to-gamma band networks. For theta-gamma band networks, we also found a plenitude of interlayer connections between distant nodes, though weaker connectivity relative to the one-to-one connections. Our work is an stepping stone towards the identification of interdependencies across frequency-specific networks. Our results lay the ground for the use of the promising multilayer framework in this field with more-informed and justified interlayer connections.