The role of epidemic spreading in seizure dynamics and epilepsy surgery.

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.

Tau protein spreads through functionally connected neurons in Alzheimer's disease: a combined MEG/PET study.

Recent studies on Alzheimer's disease (AD) suggest that tau proteins spread through the brain following neuronal connections. Several mechanisms could be involved in this process: spreading between brain regions that interact strongly (functional connectivity); through the pattern of anatomical connections (structural connectivity); or simple diffusion. Using magnetoencephalography (MEG), we investigated which spreading pathways influence tau protein spreading by modelling the tau propagation process using an epidemic spreading model. We compared the modelled tau depositions with 18F-flortaucipir PET binding potentials at several stages of the AD continuum. In this cross-sectional study, we analysed source-reconstructed MEG data and dynamic 100-min 18F-flortaucipir PET from 57 subjects positive for amyloid-ß pathology [preclinical AD (n = 16), mild cognitive impairment (MCI) due to AD (n = 16) and AD dementia (n = 25)]. Cognitively healthy subjects without amyloid-ß pathology were included as controls (n = 25). Tau propagation was modelled as an epidemic process (susceptible-infected model) on MEG-based functional networks [in alpha (8-13 Hz) and beta (13-30 Hz) bands], a structural or diffusion network, starting from the middle and inferior temporal lobe. The group-level network of the control group was used as input for the model to predict tau deposition in three stages of the AD continuum. To assess performance, model output was compared to the group-specific tau deposition patterns as measured with 18F-flortaucipir PET. We repeated the analysis by using networks of the preceding disease stage and/or using regions with most observed tau deposition during the preceding stage as seeds. In the preclinical AD stage, the functional networks predicted most of the modelled tau-PET binding potential, with best correlations between model and tau-PET [corrected amplitude envelope correlation (AEC-c) alpha C = 0.584; AEC-c beta C = 0.569], followed by the structural network (C = 0.451) and simple diffusion (C = 0.451). Prediction accuracy declined for the MCI and AD dementia stages, although the correlation between modelled tau and tau-PET binding remained highest for the functional networks (C = 0.384; C = 0.376). Replacing the control-network with the network from the preceding disease stage and/or alternative seeds improved prediction accuracy in MCI but not in the dementia stage. These results suggest that in addition to structural connections, functional connections play an important role in tau spread, and highlight that neuronal dynamics play a key role in promoting this pathological process. Aberrant neuronal communication patterns should be taken into account when identifying targets for future therapy. Our results also suggest that this process is more important in earlier disease stages (preclinical AD/MCI); possibly, in later stages, other processes may be influential.

The dual action of glioma-derived exosomes on neuronal activity: synchronization and disruption of synchrony.

Seizures represent a frequent symptom in gliomas and significantly impact patient morbidity and quality of life. Although the pathogenesis of tumor-related seizures is not fully understood, accumulating evidence indicates a key role of the peritumoral microenvironment. Brain cancer cells interact with neurons by forming synapses with them and by releasing exosomes, cytokines, and other small molecules. Strong interactions among neurons often lead to the synchronization of their activity. In this paper, we used an in vitro model to investigate the role of exosomes released by glioma cell lines and by patient-derived glioma stem cells (GSCs). The addition of exosomes released by U87 glioma cells to neuronal cultures at day in vitro (DIV) 4, when neurons are not yet synchronous, induces synchronization. At DIV 7-12 neurons become highly synchronous, and the addition of the same exosomes disrupts synchrony. By combining Ca2+ imaging, electrical recordings from single neurons with patch-clamp electrodes, substrate-integrated microelectrode arrays, and immunohistochemistry, we show that synchronization and de-synchronization are caused by the combined effect of (i) the formation of new neuronal branches, associated with a higher expression of Arp3, (ii) the modification of synaptic efficiency, and (iii) a direct action of exosomes on the electrical properties of neurons, more evident at DIV 7-12 when the threshold for spike initiation is significantly reduced. At DIV 7-12 exosomes also selectively boost glutamatergic signaling by increasing the number of excitatory synapses. Remarkably, de-synchronization was also observed with exosomes released by glioma-associated stem cells (GASCs) from patients with low-grade glioma but not from patients with high-grade glioma, where a more variable outcome was observed. These results show that exosomes released from glioma modify the electrical properties of neuronal networks and that de-synchronization caused by exosomes from low-grade glioma can contribute to the neurological pathologies of patients with brain cancers.

Epidemic models characterize seizure propagation and the effects of epilepsy surgery in individualized brain networks based on MEG and invasive EEG recordings.

Epilepsy surgery is the treatment of choice for drug-resistant epilepsy patients. However, seizure-freedom is currently achieved in only 2/3 of the patients after surgery. In this study we have developed an individualized computational model based on MEG brain networks to explore seizure propagation and the efficacy of different virtual resections. Eventually, the goal is to obtain individualized models to optimize resection strategy and outcome. We have modelled seizure propagation as an epidemic process using the susceptible-infected (SI) model on individual brain networks derived from presurgical MEG. We included 10 patients who had received epilepsy surgery and for whom the surgery outcome at least one year after surgery was known. The model parameters were tuned in in order to reproduce the patient-specific seizure propagation patterns as recorded with invasive EEG. We defined a personalized search algorithm that combined structural and dynamical information to find resections that maximally decreased seizure propagation for a given resection size. The optimal resection for each patient was defined as the smallest resection leading to at least a 90% reduction in seizure propagation. The individualized model reproduced the basic aspects of seizure propagation for 9 out of 10 patients when using the resection area as the origin of epidemic spreading, and for 10 out of 10 patients with an alternative definition of the seed region. We found that, for 7 patients, the optimal resection was smaller than the resection area, and for 4 patients we also found that a resection smaller than the resection area could lead to a 100% decrease in propagation. Moreover, for two cases these alternative resections included nodes outside the resection area. Epidemic spreading models fitted with patient specific data can capture the fundamental aspects of clinically observed seizure propagation, and can be used to test virtual resections in silico. Combined with optimization algorithms, smaller or alternative resection strategies, that are individually targeted for each patient, can be determined with the ultimate goal to improve surgery outcome. MEG-based networks can provide a good approximation of structural connectivity for computational models of seizure propagation, and facilitate their clinical use.

The shape of memory in temporal networks.

How to best define, detect and characterize network memory, i.e. the dependence of a network's structure on its past, is currently a matter of debate. Here we show that the memory of a temporal network is inherently multidimensional, and we introduce a mathematical framework for defining and efficiently estimating the microscopic shape of memory, which characterises how the activity of each link intertwines with the activities of all other links. We validate our methodology on a range of synthetic models, and we then study the memory shape of real-world temporal networks spanning social, technological and biological systems, finding that these networks display heterogeneous memory shapes. In particular, online and offline social networks are markedly different, with the latter showing richer memory and memory scales. Our theory also elucidates the phenomenon of emergent virtual loops and provides a novel methodology for exploring the dynamically rich structure of complex systems.

Local topological moves determine global diffusion properties of hyperbolic higher-order networks.

From social interactions to the human brain, higher-order networks are key to describe the underlying network geometry and topology of many complex systems. While it is well known that network structure strongly affects its function, the role that network topology and geometry has on the emerging dynamical properties of higher-order networks is yet to be clarified. In this perspective, the spectral dimension plays a key role since it determines the effective dimension for diffusion processes on a network. Despite its relevance, a theoretical understanding of which mechanisms lead to a finite spectral dimension, and how this can be controlled, still represents a challenge and is the object of intense research. Here, we introduce two nonequilibrium models of hyperbolic higher-order networks and we characterize their network topology and geometry by investigating the intertwined appearance of small-world behavior, δ-hyperbolicity, and community structure. We show that different topological moves, determining the nonequilibrium growth of the higher-order hyperbolic network models, induce tuneable values of the spectral dimension, showing a rich phenomenology which is not displayed in random graph ensembles. In particular, we observe that, if the topological moves used to construct the higher-order network increase the area/volume ratio, then the spectral dimension continuously decreases, while the opposite effect is observed if the topological moves decrease the area/volume ratio. Our work reveals a new link between the geometry of a network and its diffusion properties, contributing to a better understanding of the complex interplay between network structure and dynamics.

Optimization of epilepsy surgery through virtual resections on individual structural brain networks.

The success of epilepsy surgery in patients with refractory epilepsy depends upon correct identification of the epileptogenic zone (EZ) and an optimal choice of the resection area. In this study we developed individualized computational models based upon structural brain networks to explore the impact of different virtual resections on the propagation of seizures. The propagation of seizures was modelled as an epidemic process [susceptible-infected-recovered (SIR) model] on individual structural networks derived from presurgical diffusion tensor imaging in 19 patients. The candidate connections for the virtual resection were all connections from the clinically hypothesized EZ, from which the seizures were modelled to start, to other brain areas. As a computationally feasible surrogate for the SIR model, we also removed the connections that maximally reduced the eigenvector centrality (EC) (large values indicate network hubs) of the hypothesized EZ, with a large reduction meaning a large effect. The optimal combination of connections to be removed for a maximal effect were found using simulated annealing. For comparison, the same number of connections were removed randomly, or based on measures that quantify the importance of a node or connection within the network. We found that 90% of the effect (defined as reduction of EC of the hypothesized EZ) could already be obtained by removing substantially less than 90% of the connections. Thus, a smaller, optimized, virtual resection achieved almost the same effect as the actual surgery yet at a considerably smaller cost, sparing on average 27.49% (standard deviation: 4.65%) of the connections. Furthermore, the maximally effective connections linked the hypothesized EZ to hubs. Finally, the optimized resection was equally or more effective than removal based on structural network characteristics both regarding reducing the EC of the hypothesized EZ and seizure spreading. The approach of using reduced EC as a surrogate for simulating seizure propagation can suggest more restrictive resection strategies, whilst obtaining an almost optimal effect on reducing seizure propagation, by taking into account the unique topology of individual structural brain networks of patients.

Growth strategy determines the memory and structural properties of brain networks.

The interplay between structure and function affects the emerging properties of many natural systems. Here we use an adaptive neural network model that couples activity and topological dynamics and reproduces the experimental temporal profiles of synaptic density observed in the brain. We prove that the existence of a transient period of relatively high synaptic connectivity is critical for the development of the system under noise circumstances, such that the resulting network can recover stored memories. Moreover, we show that intermediate synaptic densities provide optimal developmental paths with minimum energy consumption, and that ultimately it is the transient heterogeneity in the network that determines its evolution. These results could explain why the pruning curves observed in actual brain areas present their characteristic temporal profiles and they also suggest new design strategies to build biologically inspired neural networks with particular information processing capabilities.

Explosive Higher-Order Kuramoto Dynamics on Simplicial Complexes.

The higher-order interactions of complex systems, such as the brain, are captured by their simplicial complex structure and have a significant effect on dynamics. However, the existing dynamical models defined on simplicial complexes make the strong assumption that the dynamics resides exclusively on the nodes. Here we formulate the higher-order Kuramoto model which describes the interactions between oscillators placed not only on nodes but also on links, triangles, and so on. We show that higher-order Kuramoto dynamics can lead to an explosive synchronization transition by using an adaptive coupling dependent on the solenoidal and the irrotational component of the dynamics.

How Memory Conforms to Brain Development.

Nature exhibits countless examples of adaptive networks, whose topology evolves constantly coupled with the activity due to its function. The brain is an illustrative example of a system in which a dynamic complex network develops by the generation and pruning of synaptic contacts between neurons while memories are acquired and consolidated. Here, we consider a recently proposed brain developing model to study how mechanisms responsible for the evolution of brain structure affect and are affected by memory storage processes. Following recent experimental observations, we assume that the basic rules for adding and removing synapses depend on local synaptic currents at the respective neurons in addition to global mechanisms depending on the mean connectivity. In this way a feedback loop between "form" and "function" spontaneously emerges that influences the ability of the system to optimally store and retrieve sensory information in patterns of brain activity or memories. In particular, we report here that, as a consequence of such a feedback-loop, oscillations in the activity of the system among the memorized patterns can occur, depending on parameters, reminding mind dynamical processes. Such oscillations have their origin in the destabilization of memory attractors due to the pruning dynamics, which induces a kind of structural disorder or noise in the system at a long-term scale. This constantly modifies the synaptic disorder induced by the interference among the many patterns of activity memorized in the system. Such new intriguing oscillatory behavior is to be associated only to long-term synaptic mechanisms during the network evolution dynamics, and it does not depend on short-term synaptic processes, as assumed in other studies, that are not present in our model.

Synchronization in network geometries with finite spectral dimension.

Recently there is a surge of interest in network geometry and topology. Here we show that the spectral dimension plays a fundamental role in establishing a clear relation between the topological and geometrical properties of a network and its dynamics. Specifically we explore the role of the spectral dimension in determining the synchronization properties of the Kuramoto model. We show that the synchronized phase can only be thermodynamically stable for spectral dimensions above four and that phase entrainment of the oscillators can only be found for spectral dimensions greater than two. We numerically test our analytical predictions on the recently introduced model of network geometry called complex network manifolds, which displays a tunable spectral dimension.

Complex Network Geometry and Frustrated Synchronization.

The dynamics of networks of neuronal cultures has been recently shown to be strongly dependent on the network geometry and in particular on their dimensionality. However, this phenomenon has been so far mostly unexplored from the theoretical point of view. Here we reveal the rich interplay between network geometry and synchronization of coupled oscillators in the context of a simplicial complex model of manifolds called Complex Network Manifold. The networks generated by this model combine small world properties (infinite Hausdorff dimension) and a high modular structure with finite and tunable spectral dimension. We show that the networks display frustrated synchronization for a wide range of the coupling strength of the oscillators, and that the synchronization properties are directly affected by the spectral dimension of the network.

Concurrence of form and function in developing networks and its role in synaptic pruning.

A fundamental question in neuroscience is how structure and function of neural systems are related. We study this interplay by combining a familiar auto-associative neural network with an evolving mechanism for the birth and death of synapses. A feedback loop then arises leading to two qualitatively different types of behaviour. In one, the network structure becomes heterogeneous and dissasortative, and the system displays good memory performance; furthermore, the structure is optimised for the particular memory patterns stored during the process. In the other, the structure remains homogeneous and incapable of pattern retrieval. These findings provide an inspiring picture of brain structure and dynamics that is compatible with experimental results on early brain development, and may help to explain synaptic pruning. Other evolving networks-such as those of protein interactions-might share the basic ingredients for this feedback loop and other questions, and indeed many of their structural features are as predicted by our model.