Self-sustained activity in a small-world network of excitable neurons

Abnormal suppression of thermal transport by long-range interactions in networks

Heat and electricity are two fundamental forms of energy widely utilized in our daily lives. Recently, in the study of complex networks, there is growing evidence that they behave significantly different at the micro-nanoscale. Here, we use a small-world network model to investigate the effects of reconnection probability p and decay exponent α on thermal and electrical transport within the network. Our results demonstrate that the electrical transport efficiency increases by nearly one order of magnitude, while the thermal transport efficiency falls off a cliff by three to four orders of magnitude, breaking the traditional rule that shortcuts enhance energy transport in small-world networks. Furthermore, we elucidate that phonon localization is a crucial factor in the weakening of thermal transport efficiency in small-world networks by characterizing the density of states, phonon participation ratio, and nearest-neighbor spacing distribution. These insights will pave new ways for designing thermoelectric materials with high electrical conductance and low thermal conductance.

Does the brain behave like a (complex) network? I. Dynamics

Graph theory is now becoming a standard tool in system-level neuroscience. However, endowing observed brain anatomy and dynamics with a complex network structure does not entail that the brain actually works as a network. Asking whether the brain behaves as a network means asking whether network properties count. From the viewpoint of neurophysiology and, possibly, of brain physics, the most substantial issues a network structure may be instrumental in addressing relate to the influence of network properties on brain dynamics and to whether these properties ultimately explain some aspects of brain function. Here, we address the dynamical implications of complex network, examining which aspects and scales of brain activity may be understood to genuinely behave as a network. To do so, we first define the meaning of networkness, and analyse some of its implications. We then examine ways in which brain anatomy and dynamics can be endowed with a network structure and discuss possible ways in which network structure may be shown to represent a genuine organisational principle of brain activity, rather than just a convenient description of its anatomy and dynamics.

MAPT rs17649553 T allele is associated with better verbal memory and higher small-world properties in Parkinson's disease

Currently, over 90 genetic loci have been found to be associated with Parkinson's disease (PD) in genome-wide association studies, nevertheless, the effects of these genetic variants on the clinical features and brain structure of PD patients are largely unknown. This study investigated the effects of microtubule-associated protein tau (MAPT) rs17649553 (C>T), a genetic variant associated with reduced PD risk, on the clinical manifestations and brain networks of PD patients. We found MAPT rs17649553 T allele was associated with better verbal memory in PD patients. In addition, MAPT rs17649553 significantly shaped the topology of gray matter covariance network and white matter network. Both the network metrics in gray matter covariance network and white matter network were correlated with verbal memory, however, the mediation analysis showed that it was the small-world properties in white matter network that mediated the effects of MAPT rs17649553 on verbal memory. These results suggest that MAPT rs17649553 T allele is associated with higher small-world properties in structural network and better verbal memory in PD.

Reconstruction of sparse recurrent connectivity and inputs from the nonlinear dynamics of neuronal networks

Reconstructing the recurrent structural connectivity of neuronal networks is a challenge crucial to address in characterizing neuronal computations. While directly measuring the detailed connectivity structure is generally prohibitive for large networks, we develop a novel framework for reverse-engineering large-scale recurrent network connectivity matrices from neuronal dynamics by utilizing the widespread sparsity of neuronal connections. We derive a linear input-output mapping that underlies the irregular dynamics of a model network composed of both excitatory and inhibitory integrate-and-fire neurons with pulse coupling, thereby relating network inputs to evoked neuronal activity. Using this embedded mapping and experimentally feasible measurements of the firing rate as well as voltage dynamics in response to a relatively small ensemble of random input stimuli, we efficiently reconstruct the recurrent network connectivity via compressive sensing techniques. Through analogous analysis, we then recover high dimensional natural stimuli from evoked neuronal network dynamics over a short time horizon. This work provides a generalizable methodology for rapidly recovering sparse neuronal network data and underlines the natural role of sparsity in facilitating the efficient encoding of network data in neuronal dynamics.

Organic Disordered Semiconductors as Networks Embedded in Space and Energy

Organic disordered semiconductors have a growing importance because of their low cost, mechanical flexibility, and multiple applications in thermoelectric devices, biosensors, and optoelectronic devices. Carrier transport consists of variable-range hopping between localized quantum states, which are disordered in both space and energy within the Gaussian disorder model. In this paper, we model an organic disordered semiconductor system as a network embedded in both space and energy so that a node represents a localized state while a link encodes the probability (or, equivalently, the Miller-Abrahams hopping rate) for carriers to hop between nodes. The associated network Laplacian matrix allows for the study of carrier dynamics using edge-centric random walks, in which links are activated by the corresponding carrier hopping rates. Our simulation work suggests that at room temperature the network exhibits a strong propensity for small-network nature, a beneficial property that in network science is related to the ease of exchanging information, particles, or energy in many different systems. However, this is not the case at low temperature. Our analysis suggests that there could be a parallelism between the well-known dependence of carrier mobility on temperature and the potential emergence of the small-world property with increasing temperature.

Controlled generation of self-sustained oscillations in complex artificial neural networks

Spatially distinct, self-sustained oscillations in artificial neural networks are fundamental to information encoding, storage, and processing in these systems. Here, we develop a method to induce a large variety of self-sustained oscillatory patterns in artificial neural networks and a controlling strategy to switch between different patterns. The basic principle is that, given a complex network, one can find a set of nodes-the minimum feedback vertex set (mFVS), whose removal or inhibition will result in a tree-like network without any loop structure. Reintroducing a few or even a single mFVS node into the tree-like artificial neural network can recover one or a few of the loops and lead to self-sustained oscillation patterns based on these loops. Reactivating various mFVS nodes or their combinations can then generate a large number of distinct neuronal firing patterns with a broad distribution of the oscillation period. When the system is near a critical state, chaos can arise, providing a natural platform for pattern switching with remarkable flexibility. With mFVS guided control, complex networks of artificial neurons can thus be exploited as potential prototypes for local, analog type of processing paradigms.

Diffusion tensor imaging brain structural clustering patterns in major depressive disorder

Using magnetic resonance diffusion tensor imaging data from 45 patients with major depressive disorder (MDD) and 41 healthy controls (HCs), network indices based on a 246-region Brainnetcome Atlas were investigated in the two groups, and in the MDD subgroups that were subgrouped based on their duration of the disease. Correlation between the network indices and the duration of illness was also examined. Differences were observed between the MDDS subgroup (short disease duration) and the HC group, but not between the MDD and HC groups. Compared with the HCs, the clustering coefficient (CC) values of MDDS were higher in precentral gyrus, and caudal lingual gyrus; the CC of MDDL subgroup (long disease duration) was higher in postcentral gyrus and dorsal granular insula in the right hemisphere. Network resilience analyses showed that the MDDS group was higher than the HC group, representing relatively more randomized networks in the diseased brains. The correlation analyses showed that the caudal lingual gyrus in the right hemisphere and the rostral lingual gyrus in the left hemisphere were particularly correlated with disease duration. The analyses showed that duration of the illness appears to have an impact on the networking patterns. Networking abnormalities in MDD patients could be blurred or hidden by the heterogeneity of the MDD clinical subgroups. Brain plasticity may introduce a recovery effect to the abnormal network patterns seen in patients with a relative short term of the illness, as the abnormalities may disappear in MDDL .

Spatiotemporal characteristics in systems of diffusively coupled excitable slow-fast FitzHugh-Rinzel dynamical neurons

In this paper, we study an excitable, biophysical system that supports wave propagation of nerve impulses. We consider a slow-fast, FitzHugh-Rinzel neuron model where only the membrane voltage interacts diffusively, giving rise to the formation of spatiotemporal patterns. We focus on local, nonlinear excitations and diverse neural responses in an excitable one- and two-dimensional configuration of diffusively coupled FitzHugh-Rinzel neurons. The study of the emerging spatiotemporal patterns is essential in understanding the working mechanism in different brain areas. We derive analytically the coefficients of the amplitude equations in the vicinity of Hopf bifurcations and characterize various patterns, including spirals exhibiting complex geometric substructures. Furthermore, we derive analytically the condition for the development of antispirals in the neighborhood of the bifurcation point. The emergence of broken target waves can be observed to form spiral-like profiles. The spatial dynamics of the excitable system exhibits two- and multi-arm spirals for small diffusive couplings. Our results reveal a multitude of neural excitabilities and possible conditions for the emergence of spiral-wave formation. Finally, we show that the coupled excitable systems with different firing characteristics participate in a collective behavior that may contribute significantly to irregular neural dynamics.

Emergent dynamics in excitable flow systems

Flow networks can describe many natural and artificial systems. We present a model for a flow system that allows for volume accumulation, includes conduits with a nonlinear relation between current and pressure difference, and can be applied to networks of arbitrary topology. The model displays complex dynamics, including self-sustained oscillations in the absence of any dynamics in the inputs and outputs. In this work we analytically show the origin of self-sustained oscillations for the one-dimensional case. We numerically study the behavior of systems of arbitrary topology under different conditions: we discuss their excitability, the effect of different boundary conditions, and wave propagation when the network has regions of conduits with linear conductance.

Emergence of spatially periodic diffusive waves in small-world neuronal networks

It has been observed in experiment that the anatomical structure of neuronal networks in the brain possesses the feature of small-world networks. Yet how the small-world structure affects network dynamics remains to be fully clarified. Here we study the dynamics of a class of small-world networks consisting of pulse-coupled integrate-and-fire (I&F) neurons. Under stochastic Poisson drive, we find that the activity of the entire network resembles diffusive waves. To understand its underlying mechanism, we analyze the simplified regular-lattice network consisting of firing-rate-based neurons as an approximation to the original I&F small-world network. We demonstrate both analytically and numerically that, with strongly coupled connections, in the absence of noise, the activity of the firing-rate-based regular-lattice network spatially forms a static grating pattern that corresponds to the spatial distribution of the firing rate observed in the I&F small-world neuronal network. We further show that the spatial grating pattern with different phases comprise the continuous attractor of both the I&F small-world and firing-rate-based regular-lattice network dynamics. In the presence of input noise, the activity of both networks is perturbed along the continuous attractor, which gives rise to the diffusive waves. Our numerical simulations and theoretical analysis may potentially provide insights into the understanding of the generation of wave patterns observed in cortical networks.

Winfree loop sustained oscillation in two-dimensional excitable lattices: Prediction and realization

The problem of self-sustained oscillations in excitable complex networks is the central issue under investigation, among which the prediction and the realization of self-sustained oscillations in different kinds of excitable networks are the challenging tasks. In this paper, we extensively investigate the prediction and the realization of a Winfree loop sustained oscillation (WLSO) in two-dimensional (2D) excitable lattices. By analyzing the network structure, the fundamental oscillation source structure (FOSS) of WLSO in a 2D excitable lattice is exposed explicitly. For the suitable combinations of system parameters, the Winfree loop can self-organize on the FOSS to form an oscillation source sustaining the oscillation, and these suitable parameter combinations are predicted by calculating the minimum Winfree loop length and have been further confirmed in numerical simulations. However, the FOSS cannot spontaneously offer the WLSO in 2D excitable lattices in usual cases due to the coupling bidirectionality and the symmetry properties of the lattice. A targeted protection scheme of the oscillation source is proposed by overcoming these two drawbacks. Finally, the WLSO is realized in the 2D excitable lattice successfully.

Network structure and input integration in competing firing rate models for decision-making

Making a decision among numerous alternatives is a pervasive and central undertaking encountered by mammals in natural settings. While decision making for two-option tasks has been studied extensively both experimentally and theoretically, characterizing decision making in the face of a large set of alternatives remains challenging. We explore this issue by formulating a scalable mechanistic network model for decision making and analyzing the dynamics evoked given various potential network structures. In the case of a fully-connected network, we provide an analytical characterization of the model fixed points and their stability with respect to winner-take-all behavior for fair tasks. We compare several means of input integration, demonstrating a more gradual sigmoidal transfer function is likely evolutionarily advantageous relative to binary gain commonly utilized in engineered systems. We show via asymptotic analysis and numerical simulation that sigmoidal transfer functions with smaller steepness yield faster response times but depreciation in accuracy. However, in the presence of noise or degradation of connections, a sigmoidal transfer function garners significantly more robust and accurate decision-making dynamics. For fair tasks and sigmoidal gain, our model network also exhibits a stable parameter regime that produces high accuracy and persists across tasks with diverse numbers of alternatives and difficulties, satisfying physiological energetic constraints. In the case of more sparse and structured network topologies, including random, regular, and small-world connectivity, we show the high-accuracy parameter regime persists for biologically realistic connection densities. Our work shows how neural system architecture is potentially optimal in making economic, reliable, and advantageous decisions across tasks.

The impact of spike-frequency adaptation on balanced network dynamics

A dynamic balance between strong excitatory and inhibitory neuronal inputs is hypothesized to play a pivotal role in information processing in the brain. While there is evidence of the existence of a balanced operating regime in several cortical areas and idealized neuronal network models, it is important for the theory of balanced networks to be reconciled with more physiological neuronal modeling assumptions. In this work, we examine the impact of spike-frequency adaptation, observed widely across neurons in the brain, on balanced dynamics. We incorporate adaptation into binary and integrate-and-fire neuronal network models, analyzing the theoretical effect of adaptation in the large network limit and performing an extensive numerical investigation of the model adaptation parameter space. Our analysis demonstrates that balance is well preserved for moderate adaptation strength even if the entire network exhibits adaptation. In the common physiological case in which only excitatory neurons undergo adaptation, we show that the balanced operating regime in fact widens relative to the non-adaptive case. We hypothesize that spike-frequency adaptation may have been selected through evolution to robustly facilitate balanced dynamics across diverse cognitive operating states.

Effect of inhibitory spike-timing-dependent plasticity on fast sparsely synchronized rhythms in a small-world neuronal network

We consider the Watts-Strogatz small-world network (SWN) consisting of inhibitory fast spiking Izhikevich interneurons. This inhibitory neuronal population has adaptive dynamic synaptic strengths governed by the inhibitory spike-timing-dependent plasticity (iSTDP). In previous works without iSTDP, fast sparsely synchronized rhythms, associated with diverse cognitive functions, were found to appear in a range of large noise intensities for fixed strong synaptic inhibition strengths. Here, we investigate the effect of iSTDP on fast sparse synchronization (FSS) by varying the noise intensity D. We employ an asymmetric anti-Hebbian time window for the iSTDP update rule [which is in contrast to the Hebbian time window for the excitatory STDP (eSTDP)]. Depending on values of D, population-averaged values of saturated synaptic inhibition strengths are potentiated [long-term potentiation (LTP)] or depressed [long-term depression (LTD)] in comparison with the initial mean value, and dispersions from the mean values of LTP/LTD are much increased when compared with the initial dispersion, independently of D. In most cases of LTD where the effect of mean LTD is dominant in comparison with the effect of dispersion, good synchronization (with higher spiking measure) is found to get better via LTD, while bad synchronization (with lower spiking measure) is found to get worse via LTP. This kind of Matthew effect in inhibitory synaptic plasticity is in contrast to that in excitatory synaptic plasticity where good (bad) synchronization gets better (worse) via LTP (LTD). Emergences of LTD and LTP of synaptic inhibition strengths are intensively investigated via a microscopic method based on the distributions of time delays between the pre- and the post-synaptic spike times. Furthermore, we also investigate the effects of network architecture on FSS by changing the rewiring probability p of the SWN in the presence of iSTDP.

Circuit variability interacts with excitatory-inhibitory diversity of interneurons to regulate network encoding capacity

Local interneurons (LNs) in the Drosophila olfactory system exhibit neuronal diversity and variability, yet it is still unknown how these features impact information encoding capacity and reliability in a complex LN network. We employed two strategies to construct a diverse excitatory-inhibitory neural network beginning with a ring network structure and then introduced distinct types of inhibitory interneurons and circuit variability to the simulated network. The continuity of activity within the node ensemble (oscillation pattern) was used as a readout to describe the temporal dynamics of network activity. We found that inhibitory interneurons enhance the encoding capacity by protecting the network from extremely short activation periods when the network wiring complexity is very high. In addition, distinct types of interneurons have differential effects on encoding capacity and reliability. Circuit variability may enhance the encoding reliability, with or without compromising encoding capacity. Therefore, we have described how circuit variability of interneurons may interact with excitatory-inhibitory diversity to enhance the encoding capacity and distinguishability of neural networks. In this work, we evaluate the effects of different types and degrees of connection diversity on a ring model, which may simulate interneuron networks in the Drosophila olfactory system or other biological systems.

Neuronal Population Activity in Spinal Motor Circuits: Greater Than the Sum of Its Parts

The core elements of stereotypical movements such as locomotion, scratching and breathing are generated by networks in the lower brainstem and the spinal cord. Ensemble activities in spinal motor networks had until recently been merely a black box, but with the emergence of ultra-thin Silicon multi-electrode technology it was possible to reveal the spiking activity of larger parts of the network. A series of experiments revealed unexpected features of spinal networks, such as multiple spiking regimes and lognormal firing rate distributions. The lognormality renders the widespread idea of a typical firing rate ± standard deviation an ill-suited description, and therefore these findings define a new arithmetic of motor networks. Focusing on the population activity behind motor pattern generation this review summarizes this advance and discusses its implications.

High-energy brain dynamics during anesthesia-induced unconsciousness

Characterizing anesthesia-induced alterations to brain network dynamics provides a powerful framework to understand the neural mechanisms of unconsciousness. To this end, increased attention has been directed at how anesthetic drugs alter the functional connectivity between brain regions as defined through neuroimaging. However, the effects of anesthesia on temporal dynamics at functional network scales is less well understood. Here, we examine such dynamics in view of the free-energy principle, which postulates that brain dynamics tend to promote lower energy (more organized) states. We specifically engaged the hypothesis that such low-energy states play an important role in maintaining conscious awareness. To investigate this hypothesis, we analyzed resting-state BOLD fMRI data from human volunteers during wakefulness and under sevoflurane general anesthesia. Our approach, which extends an idea previously used in the characterization of neuron-scale populations, involves thresholding the BOLD time series and using a normalized Hamiltonian energy function derived from the Ising model. Our major finding is that the brain spends significantly more time in lower energy states during eyes-closed wakefulness than during general anesthesia. This effect is especially pronounced in networks thought to be critical for maintaining awareness, suggesting a crucial cognitive role for both the structure and the dynamical landscape of these networks.

We show that activity in the human brain, as captured by functional magnetic resonance imaging (fMRI), is more organized during wakefulness than during general anesthesia. This increased organization corresponds to a decrease in a statistical-physics-inspired energy measure among brain regions of shared functional specialization (resting-state networks) that have putative roles in conscious awareness and attention. Characterizing the energy distributions in this way reveals significant changes in the dynamics of brain activity in different states of consciousness, insights that are not observable in the average functional connectivity data alone. Our results indicate that the ability of brain networks to sustain stable representations, via their dynamics, may be crucial for consciousness and cognition.

Stochastic spike synchronization in a small-world neural network with spike-timing-dependent plasticity

We consider the Watts-Strogatz small-world network (SWN) consisting of subthreshold neurons which exhibit noise-induced spikings. This neuronal network has adaptive dynamic synaptic strengths governed by the spike-timing-dependent plasticity (STDP). In previous works without STDP, stochastic spike synchronization (SSS) between noise-induced spikings of subthreshold neurons was found to occur in a range of intermediate noise intensities. Here, we investigate the effect of additive STDP on the SSS by varying the noise intensity. Occurrence of a "Matthew" effect in synaptic plasticity is found due to a positive feedback process. As a result, good synchronization gets better via long-term potentiation of synaptic strengths, while bad synchronization gets worse via long-term depression. Emergences of long-term potentiation and long-term depression of synaptic strengths are intensively investigated via microscopic studies based on the pair-correlations between the pre- and the post-synaptic IISRs (instantaneous individual spike rates) as well as the distributions of time delays between the pre- and the post-synaptic spike times. Furthermore, the effects of multiplicative STDP (which depends on states) on the SSS are studied and discussed in comparison with the case of additive STDP (independent of states). These effects of STDP on the SSS in the SWN are also compared with those in the regular lattice and the random graph.

Spatiotemporal dynamics in excitable homogeneous random networks composed of periodically self-sustained oscillation

The collective behaviors of networks are often dependent on the network connections and bifurcation parameters, also the local kinetics plays an important role in contributing the consensus of coupled oscillators. In this paper, we systematically investigate the influence of network structures and system parameters on the spatiotemporal dynamics in excitable homogeneous random networks (EHRNs) composed of periodically self-sustained oscillation (PSO). By using the dominant phase-advanced driving (DPAD) method, the one-dimensional (1D) Winfree loop is exposed as the oscillation source supporting the PSO, and the accurate wave propagation pathways from the oscillation source to the whole network are uncovered. Then, an order parameter is introduced to quantitatively study the influence of network structures and system parameters on the spatiotemporal dynamics of PSO in EHRNs. Distinct results induced by the network structures and the system parameters are observed. Importantly, the corresponding mechanisms are revealed. PSO influenced by the network structures are induced not only by the change of average path length (APL) of network, but also by the invasion of 1D Winfree loop from the outside linking nodes. Moreover, PSO influenced by the system parameters are determined by the excitation threshold and the minimum 1D Winfree loop. Finally, we confirmed that the excitation threshold and the minimum 1D Winfree loop determined PSO will degenerate as the system size is expanded.

Statistical complexity is maximized in a small-world brain

In this paper, we study a network of Izhikevich neurons to explore what it means for a brain to be at the edge of chaos. To do so, we first constructed the phase diagram of a single Izhikevich excitatory neuron, and identified a small region of the parameter space where we find a large number of phase boundaries to serve as our edge of chaos. We then couple the outputs of these neurons directly to the parameters of other neurons, so that the neuron dynamics can drive transitions from one phase to another on an artificial energy landscape. Finally, we measure the statistical complexity of the parameter time series, while the network is tuned from a regular network to a random network using the Watts-Strogatz rewiring algorithm. We find that the statistical complexity of the parameter dynamics is maximized when the neuron network is most small-world-like. Our results suggest that the small-world architecture of neuron connections in brains is not accidental, but may be related to the information processing that they do.

Minimum Winfree loop determines self-sustained oscillations in excitable Erdös-Rényi random networks

The investigation of self-sustained oscillations in excitable complex networks is very important in understanding various activities in brain systems, among which the exploration of the key determinants of oscillations is a challenging task. In this paper, by investigating the influence of system parameters on self-sustained oscillations in excitable Erdös-Rényi random networks (EERRNs), the minimum Winfree loop (MWL) is revealed to be the key factor in determining the emergence of collective oscillations. Specifically, the one-to-one correspondence between the optimal connection probability (OCP) and the MWL length is exposed. Moreover, many important quantities such as the lower critical connection probability (LCCP), the OCP, and the upper critical connection probability (UCCP) are determined by the MWL. Most importantly, they can be approximately predicted by the network structure analysis, which have been verified in numerical simulations. Our results will be of great importance to help us in understanding the key factors in determining persistent activities in biological systems.

Ablation as targeted perturbation to rewire communication network of persistent atrial fibrillation

Persistent atrial fibrillation (AF) can be viewed as disintegrated patterns of information transmission by action potential across the communication network consisting of nodes linked by functional connectivity. To test the hypothesis that ablation of persistent AF is associated with improvement in both local and global connectivity within the communication networks, we analyzed multi-electrode basket catheter electrograms of 22 consecutive patients (63.5 ± 9.7 years, 78% male) during persistent AF before and after the focal impulse and rotor modulation-guided ablation. Eight patients (36%) developed recurrence within 6 months after ablation. We defined communication networks of AF by nodes (cardiac tissue adjacent to each electrode) and edges (mutual information between pairs of nodes). To evaluate patient-specific parameters of communication, thresholds of mutual information were applied to preserve 10% to 30% of the strongest edges. There was no significant difference in network parameters between both atria at baseline. Ablation effectively rewired the communication network of persistent AF to improve the overall connectivity. In addition, successful ablation improved local connectivity by increasing the average clustering coefficient, and also improved global connectivity by decreasing the characteristic path length. As a result, successful ablation improved the efficiency and robustness of the communication network by increasing the small-world index. These changes were not observed in patients with AF recurrence. Furthermore, a significant increase in the small-world index after ablation was associated with synchronization of the rhythm by acute AF termination. In conclusion, successful ablation rewires communication networks during persistent AF, making it more robust, efficient, and easier to synchronize. Quantitative analysis of communication networks provides not only a mechanistic insight that AF may be sustained by spatially localized sources and global connectivity, but also patient-specific metrics that could serve as a valid endpoint for therapeutic interventions.

Synaptic Impairment and Robustness of Excitatory Neuronal Networks with Different Topologies

Synaptic deficiencies are a known hallmark of neurodegenerative diseases, but the diagnosis of impaired synapses on the cellular level is not an easy task. Nonetheless, changes in the system-level dynamics of neuronal networks with damaged synapses can be detected using techniques that do not require high spatial resolution. This paper investigates how the structure/topology of neuronal networks influences their dynamics when they suffer from synaptic loss. We study different neuronal network structures/topologies by specifying their degree distributions. The modes of the degree distribution can be used to construct networks that consist of rich clubs and resemble small world networks, as well. We define two dynamical metrics to compare the activity of networks with different structures: persistent activity (namely, the self-sustained activity of the network upon removal of the initial stimulus) and quality of activity (namely, percentage of neurons that participate in the persistent activity of the network). Our results show that synaptic loss affects the persistent activity of networks with bimodal degree distributions less than it affects random networks. The robustness of neuronal networks enhances when the distance between the modes of the degree distribution increases, suggesting that the rich clubs of networks with distinct modes keep the whole network active. In addition, a tradeoff is observed between the quality of activity and the persistent activity. For a range of distributions, both of these dynamical metrics are considerably high for networks with bimodal degree distribution compared to random networks. We also propose three different scenarios of synaptic impairment, which may correspond to different pathological or biological conditions. Regardless of the network structure/topology, results demonstrate that synaptic loss has more severe effects on the activity of the network when impairments are correlated with the activity of the neurons.

Development of structural correlations and synchronization from adaptive rewiring in networks of Kuramoto oscillators

Synchronization of non-identical oscillators coupled through complex networks is an important example of collective behavior, and it is interesting to ask how the structural organization of network interactions influences this process. Several studies have explored and uncovered optimal topologies for synchronization by making purposeful alterations to a network. On the other hand, the connectivity patterns of many natural systems are often not static, but are rather modulated over time according to their dynamics. However, this co-evolution and the extent to which the dynamics of the individual units can shape the organization of the network itself are less well understood. Here, we study initially randomly connected but locally adaptive networks of Kuramoto oscillators. In particular, the system employs a co-evolutionary rewiring strategy that depends only on the instantaneous, pairwise phase differences of neighboring oscillators, and that conserves the total number of edges, allowing the effects of local reorganization to be isolated. We find that a simple rule-which preserves connections between more out-of-phase oscillators while rewiring connections between more in-phase oscillators-can cause initially disordered networks to organize into more structured topologies that support enhanced synchronization dynamics. We examine how this process unfolds over time, finding a dependence on the intrinsic frequencies of the oscillators, the global coupling, and the network density, in terms of how the adaptive mechanism reorganizes the network and influences the dynamics. Importantly, for large enough coupling and after sufficient adaptation, the resulting networks exhibit interesting characteristics, including degree-frequency and frequency-neighbor frequency correlations. These properties have previously been associated with optimal synchronization or explosive transitions in which the networks were constructed using global information. On the contrary, by considering a time-dependent interplay between structure and dynamics, this work offers a mechanism through which emergent phenomena and organization can arise in complex systems utilizing local rules.

Optimal Detection of a Localized Perturbation in Random Networks of Integrate-and-Fire Neurons

Experimental and theoretical studies suggest that cortical networks are chaotic and coding relies on averages over large populations. However, there is evidence that rats can respond to the short stimulation of a single cortical cell, a theoretically unexplained fact. We study effects of single-cell stimulation on a large recurrent network of integrate-and-fire neurons and propose a simple way to detect the perturbation. Detection rates obtained from simulations and analytical estimates are similar to experimental response rates if the readout is slightly biased towards specific neurons. Near-optimal detection is attained for a broad range of intermediate values of the mean coupling between neurons.

Emergence of ultrafast sparsely synchronized rhythms and their responses to external stimuli in an inhomogeneous small-world complex neuronal network

We consider an inhomogeneous small-world network (SWN) composed of inhibitory short-range (SR) and long-range (LR) interneurons, and investigate the effect of network architecture on emergence of synchronized brain rhythms by varying the fraction of LR interneurons p_long. The betweenness centralities of the LR and SR interneurons (characterizing the potentiality in controlling communication between other interneurons) are distinctly different. Hence, in view of the betweenness, SWNs we consider are inhomogeneous, unlike the "canonical" Watts-Strogatz SWN with nearly the same betweenness centralities. For small p_long, the load of communication traffic is much concentrated on a few LR interneurons. However, as p_long is increased, the number of LR connections (coming from LR interneurons) increases, and then the load of communication traffic is less concentrated on LR interneurons, which leads to better efficiency of global communication between interneurons. Sparsely synchronized rhythms are thus found to emerge when passing a small critical value p_long^(c)(≃0.16). The population frequency of the sparsely synchronized rhythm is ultrafast (higher than 100 Hz), while the mean firing rate of individual interneurons is much lower (∼30 Hz) due to stochastic and intermittent neural discharges. These dynamical behaviors in the inhomogeneous SWN are also compared with those in the homogeneous Watts-Strogatz SWN, in connection with their network topologies. Particularly, we note that the main difference between the two types of SWNs lies in the distribution of betweenness centralities. Unlike the case of the Watts-Strogatz SWN, dynamical responses to external stimuli vary depending on the type of stimulated interneurons in the inhomogeneous SWN. We consider two cases of external time-periodic stimuli applied to sub-populations of the LR and SR interneurons, respectively. Dynamical responses (such as synchronization suppression and enhancement) to these two cases of stimuli are studied and discussed in relation to the betweenness centralities of stimulated interneurons, representing the effectiveness for transfer of stimulation effect in the whole network.

A small change in neuronal network topology can induce explosive synchronization transition and activity propagation in the entire network

We here study explosive synchronization transitions and network activity propagation in networks of coupled neurons to provide a new understanding of the relationship between network topology and explosive dynamical transitions as in epileptic seizures and their propagations in the brain. We model local network motifs and configurations of coupled neurons and analyze the activity propagations between a group of active neurons to their inactive neuron neighbors in a variety of network configurations. We find that neuronal activity propagation is limited to local regions when network is highly clustered with modular structures as in the normal brain networks. When the network cluster structure is slightly changed, the activity propagates to the entire network, which is reminiscent of epileptic seizure propagation in the brain. Finally, we analyze intracranial electroencephalography (IEEG) recordings of a seizure episode from a epilepsy patient and uncover that explosive synchronization-like transition occurs around the clinically defined onset of seizure. These findings may provide a possible mechanism for the recurrence of epileptic seizures, which are known to be the results of aberrant neuronal network structure and/or function in the brain.

Autaptic Connections Shift Network Excitability and Bursting

We examine the role of structural autapses, when a neuron synapses onto itself, in driving network-wide bursting behavior. Using a simple spiking model of neuronal activity, we study how autaptic connections affect activity patterns, and evaluate if controllability significantly affects changes in bursting from autaptic connections. Adding more autaptic connections to excitatory neurons increased the number of spiking events and the number of network-wide bursts. We observed excitatory synapses contributed more to bursting behavior than inhibitory synapses. We evaluated if neurons with high average controllability, predicted to push the network into easily achievable states, affected bursting behavior differently than neurons with high modal controllability, thought to influence the network into difficult to reach states. Results show autaptic connections to excitatory neurons with high average controllability led to higher burst frequencies than adding the same number of self-looping connections to neurons with high modal controllability. The number of autapses required to induce bursting was lowered by adding autapses to high degree excitatory neurons. These results suggest a role of autaptic connections in controlling network-wide bursts in diverse cortical and subcortical regions of mammalian brain. Moreover, they open up new avenues for the study of dynamic neurophysiological correlates of structural controllability.

The shaping of intrinsic membrane potential oscillations: positive/negative feedback, ionic resonance/amplification, nonlinearities and time scales

The generation of intrinsic subthreshold (membrane potential) oscillations (STOs) in neuronal models requires the interaction between two processes: a relatively fast positive feedback that favors changes in voltage and a slower negative feedback that opposes these changes. These are provided by the so-called resonant and amplifying gating variables associated to the participating ionic currents. We investigate both the biophysical and dynamic mechanisms of generation of STOs and how their attributes (frequency and amplitude) depend on the model parameters for biophysical (conductance-based) models having qualitatively different types of resonant currents (activating and inactivating) and an amplifying current. Combinations of the same types of ionic currents (same models) in different parameter regimes give rise to different types of nonlinearities in the voltage equation: quasi-linear, parabolic-like and cubic-like. On the other hand, combinations of different types of ionic currents (different models) may give rise to the same type of nonlinearities. We examine how the attributes of the resulting STOs depend on the combined effect of these resonant and amplifying ionic processes, operating at different effective time scales, and the various types of nonlinearities. We find that, while some STO properties and attribute dependencies on the model parameters are determined by the specific combinations of ionic currents (biophysical properties), and are different for models with different such combinations, others are determined by the type of nonlinearities and are common for models with different types of ionic currents. Our results highlight the richness of STO behavior in single cells as the result of the various ways in which resonant and amplifying currents interact and affect the generation and termination of STOs as control parameters change. We make predictions that can be tested experimentally and are expected to contribute to the understanding of how rhythmic activity in neuronal networks emerge from the interplay of the intrinsic properties of the participating neurons and the network connectivity.

Pattern Selection in Network of Coupled Multi-Scroll Attractors

Multi-scroll chaotic attractor makes the oscillator become more complex in dynamic behaviors. The collective behaviors of coupled oscillators with multi-scroll attractors are investigated in the regular network in two-dimensional array, which the local kinetics is described by an improved Chua circuit. A feasible scheme of negative feedback with diversity is imposed on the network to stabilize the spatial patterns. Firstly, the Chua circuit is improved by replacing the nonlinear term with Sine function to generate infinite aquariums so that multi-scroll chaotic attractors could be generated under appropriate parameters, which could be detected by calculating the Lyapunov exponent in the parameter region. Furthermore, negative feedback with different gains (D1, D2) is imposed on the local square center area A2 and outer area A1 of the network, it is found that spiral wave, target wave could be developed in the network under appropriate feedback gain with diversity and size of controlled area. Particularly, homogeneous state could be reached after synchronization by selecting appropriate feedback gain and controlled size in the network. Finally, the distribution for statistical factors of synchronization is calculated in the two-parameter space to understand the transition of pattern region. It is found that developed spiral waves, target waves often are associated with smaller factor of synchronization. These results show that emergence of sustained spiral wave and continuous target wave could be effective for further suppression of spatiotemporal chaos in network by generating stable pacemaker completely.

Small-World Propensity and Weighted Brain Networks

Quantitative descriptions of network structure can provide fundamental insights into the function of interconnected complex systems. Small-world structure, diagnosed by high local clustering yet short average path length between any two nodes, promotes information flow in coupled systems, a key function that can differ across conditions or between groups. However, current techniques to quantify small-worldness are density dependent and neglect important features such as the strength of network connections, limiting their application in real-world systems. Here, we address both limitations with a novel metric called the Small-World Propensity (SWP). In its binary instantiation, the SWP provides an unbiased assessment of small-world structure in networks of varying densities. We extend this concept to the case of weighted brain networks by developing (i) a standardized procedure for generating weighted small-world networks, (ii) a weighted extension of the SWP, and (iii) a method for mapping observed brain network data onto the theoretical model. In applying these techniques to compare real-world brain networks, we uncover the surprising fact that the canonical biological small-world network, the C. elegans neuronal network, has strikingly low SWP. These metrics, models, and maps form a coherent toolbox for the assessment and comparison of architectural properties in brain networks.

The Fundamental Structure and the Reproduction of Spiral Wave in a Two-Dimensional Excitable Lattice

In this paper we have systematically investigated the fundamental structure and the reproduction of spiral wave in a two-dimensional excitable lattice. A periodically rotating spiral wave is introduced as the model to reproduce spiral wave artificially. Interestingly, by using the dominant phase-advanced driving analysis method, the fundamental structure containing the loop structure and the wave propagation paths has been revealed, which can expose the periodically rotating orbit of spiral tip and the charity of spiral wave clearly. Furthermore, the fundamental structure is utilized as the core for artificial spiral wave. Additionally, the appropriate parameter region, in which the artificial spiral wave can be reproduced, is studied. Finally, we discuss the robustness of artificial spiral wave to defects.

Low-rank network decomposition reveals structural characteristics of small-world networks

Small-world networks occur naturally throughout biological, technological, and social systems. With their prevalence, it is particularly important to prudently identify small-world networks and further characterize their unique connection structure with respect to network function. In this work we develop a formalism for classifying networks and identifying small-world structure using a decomposition of network connectivity matrices into low-rank and sparse components, corresponding to connections within clusters of highly connected nodes and sparse interconnections between clusters, respectively. We show that the network decomposition is independent of node indexing and define associated bounded measures of connectivity structure, which provide insight into the clustering and regularity of network connections. While many existing network characterizations rely on constructing benchmark networks for comparison or fail to describe the structural properties of relatively densely connected networks, our classification relies only on the intrinsic network structure and is quite robust with respect to changes in connection density, producing stable results across network realizations. Using this framework, we analyze several real-world networks and reveal new structural properties, which are often indiscernible by previously established characterizations of network connectivity.

Memory-induced mechanism for self-sustaining activity in networks

We study a mechanism of activity sustaining on networks inspired by a well-known model of neuronal dynamics. Our primary focus is the emergence of self-sustaining collective activity patterns, where no single node can stay active by itself, but the activity provided initially is sustained within the collective of interacting agents. In contrast to existing models of self-sustaining activity that are caused by (long) loops present in the network, here we focus on treelike structures and examine activation mechanisms that are due to temporal memory of the nodes. This approach is motivated by applications in social media, where long network loops are rare or absent. Our results suggest that under a weak behavioral noise, the nodes robustly split into several clusters, with partial synchronization of nodes within each cluster. We also study the randomly weighted version of the models where the nodes are allowed to change their connection strength (this can model attention redistribution) and show that it does facilitate the self-sustained activity.

Effect of intermodular connection on fast sparse synchronization in clustered small-world neural networks

We consider a clustered network with small-world subnetworks of inhibitory fast spiking interneurons and investigate the effect of intermodular connection on the emergence of fast sparsely synchronized rhythms by varying both the intermodular coupling strength J(inter) and the average number of intermodular links per interneuron M(syn)(inter). In contrast to the case of nonclustered networks, two kinds of sparsely synchronized states such as modular and global synchronization are found. For the case of modular sparse synchronization, the population behavior reveals the modular structure, because the intramodular dynamics of subnetworks make some mismatching. On the other hand, in the case of global sparse synchronization, the population behavior is globally identical, independently of the cluster structure, because the intramodular dynamics of subnetworks make perfect matching. We introduce a realistic cross-correlation modularity measure, representing the matching degree between the instantaneous subpopulation spike rates of the subnetworks, and examine whether the sparse synchronization is global or modular. Depending on its magnitude, the intermodular coupling strength J(inter) seems to play "dual" roles for the pacing between spikes in each subnetwork. For large J(inter), due to strong inhibition it plays a destructive role to "spoil" the pacing between spikes, while for small J(inter) it plays a constructive role to "favor" the pacing between spikes. Through competition between the constructive and the destructive roles of J(inter), there exists an intermediate optimal J(inter) at which the pacing degree between spikes becomes maximal. In contrast, the average number of intermodular links per interneuron M(syn)(inter) seems to play a role just to favor the pacing between spikes. With increasing M(syn)(inter), the pacing degree between spikes increases monotonically thanks to the increase in the degree of effectiveness of global communication between spikes. Furthermore, we employ the realistic sub- and whole-population order parameters, based on the instantaneous sub- and whole-population spike rates, to determine the threshold values for the synchronization-unsynchronization transition in the sub- and whole populations, and the degrees of global and modular sparse synchronization are also measured in terms of the realistic sub- and whole-population statistical-mechanical spiking measures defined by considering both the occupation and the pacing degrees of spikes. It is expected that our results could have implications for the role of the brain plasticity in some functional behaviors associated with population synchronization.

Broadband macroscopic cortical oscillations emerge from intrinsic neuronal response failures

Broadband spontaneous macroscopic neural oscillations are rhythmic cortical firing which were extensively examined during the last century, however, their possible origination is still controversial. In this work we show how macroscopic oscillations emerge in solely excitatory random networks and without topological constraints. We experimentally and theoretically show that these oscillations stem from the counterintuitive underlying mechanism-the intrinsic stochastic neuronal response failures (NRFs). These NRFs, which are characterized by short-term memory, lead to cooperation among neurons, resulting in sub- or several- Hertz macroscopic oscillations which coexist with high frequency gamma oscillations. A quantitative interplay between the statistical network properties and the emerging oscillations is supported by simulations of large networks based on single-neuron in-vitro experiments and a Langevin equation describing the network dynamics. Results call for the examination of these oscillations in the presence of inhibition and external drives.

Global segregation of cortical activity and metastable dynamics

Cortical activity exhibits persistent metastable dynamics. Assemblies of neurons transiently couple (integrate) and decouple (segregate) at multiple spatiotemporal scales; both integration and segregation are required to support metastability. Integration of distant brain regions can be achieved through long range excitatory projections, but the mechanism supporting long range segregation is not clear. We argue that the thalamocortical matrix connections, which project diffusely from the thalamus to the cortex and have long been thought to support cortical gain control, play an equally-important role in cortical segregation. We present a computational model of the diffuse thalamocortical loop, called the competitive cross-coupling (CXC) spiking network. Simulations of the model show how different levels of tonic input from the brainstem to the thalamus could control dynamical complexity in the cortex, directing transitions between sleep, wakefulness and high attention or vigilance. The model also explains how mutually-exclusive activity could arise across large portions of the cortex, such as between the default-mode and task-positive networks. It is robust to noise but does not require noise to autonomously generate metastability. We conclude that the long range segregation observed in brain activity and required for global metastable dynamics could be provided by the thalamocortical matrix, and is strongly modulated by brainstem input to the thalamus.

A framework for analyzing the relationship between gene expression and morphological, topological, and dynamical patterns in neuronal networks

BACKGROUND

A key point in developmental biology is to understand how gene expression influences the morphological and dynamical patterns that are observed in living beings.

NEW METHOD

In this work we propose a methodology capable of addressing this problem that is based on estimating the mutual information and Pearson correlation between the intensity of gene expression and measurements of several morphological properties of the cells. A similar approach is applied in order to identify effects of gene expression over the system dynamics. Neuronal networks were artificially grown over a lattice by considering a reference model used to generate artificial neurons. The input parameters of the artificial neurons were determined according to two distinct patterns of gene expression and the dynamical response was assessed by considering the integrate-and-fire model.

RESULTS

As far as single gene dependence is concerned, we found that the interaction between the gene expression and the network topology, as well as between the former and the dynamics response, is strongly affected by the gene expression pattern. In addition, we observed a high correlation between the gene expression and some topological measurements of the neuronal network for particular patterns of gene expression.

COMPARISON WITH EXISTING METHODS

To our best understanding, there are no similar analyses to compare with.

CONCLUSIONS

A proper understanding of gene expression influence requires jointly studying the morphology, topology, and dynamics of neurons. The proposed framework represents a first step towards predicting gene expression patterns from morphology and connectivity.

Measuring predictability of autonomous network transitions into bursting dynamics

Understanding spontaneous transitions between dynamical modes in a network is of significant importance. These transitions may separate pathological and normal functions of the brain. In this paper, we develop a set of measures that, based on spatio-temporal features of network activity, predict autonomous network transitions from asynchronous to synchronous dynamics under various conditions. These metrics quantify spike-timing distributions within a narrow time window as a function of the relative location of the active neurons. We applied these metrics to investigate the properties of these transitions in excitatory-only and excitatory-and-inhibitory networks and elucidate how network topology, noise level, and cellular heterogeneity affect both the reliability and the timeliness of the predictions. The developed measures can be calculated in real time and therefore potentially applied in clinical situations.

Noise-induced burst and spike synchronizations in an inhibitory small-world network of subthreshold bursting neurons

We are interested in noise-induced firings of subthreshold neurons which may be used for encoding environmental stimuli. Noise-induced population synchronization was previously studied only for the case of global coupling, unlike the case of subthreshold spiking neurons. Hence, we investigate the effect of complex network architecture on noise-induced synchronization in an inhibitory population of subthreshold bursting Hindmarsh-Rose neurons. For modeling complex synaptic connectivity, we consider the Watts-Strogatz small-world network which interpolates between regular lattice and random network via rewiring, and investigate the effect of small-world connectivity on emergence of noise-induced population synchronization. Thus, noise-induced burst synchronization (synchrony on the slow bursting time scale) and spike synchronization (synchrony on the fast spike time scale) are found to appear in a synchronized region of the [Formula: see text]-[Formula: see text] plane ([Formula: see text]: synaptic inhibition strength and [Formula: see text]: noise intensity). As the rewiring probability [Formula: see text] is decreased from 1 (random network) to 0 (regular lattice), the region of spike synchronization shrinks rapidly in the [Formula: see text]-[Formula: see text] plane, while the region of the burst synchronization decreases slowly. We separate the slow bursting and the fast spiking time scales via frequency filtering, and characterize the noise-induced burst and spike synchronizations by employing realistic order parameters and statistical-mechanical measures introduced in our recent work. Thus, the bursting and spiking thresholds for the burst and spike synchronization transitions are determined in terms of the bursting and spiking order parameters, respectively. Furthermore, we also measure the degrees of burst and spike synchronizations in terms of the statistical-mechanical bursting and spiking measures, respectively.

Emergence of dynamical complexity related to human heart rate variability

We apply the refined composite multiscale entropy (MSE) method to a one-dimensional directed small-world network composed of nodes whose states are binary and whose dynamics obey the majority rule. We find that the resulting fluctuating signal becomes dynamically complex. This dynamical complexity is caused (i) by the presence of both short-range connections and long-range shortcuts and (ii) by how well the system can adapt to the noisy environment. By tuning the adaptability of the environment and the long-range shortcuts we can increase or decrease the dynamical complexity, thereby modeling trends found in the MSE of a healthy human heart rate in different physiological states. When the shortcut and adaptability values increase, the complexity in the system dynamics becomes uncorrelated.

Role of long cycles in excitable dynamics on graphs

Topological cycles in excitable networks can play an important role in maintaining the network activity. When properly activated, cycles act as dynamic pacemakers, sustaining the activity of the whole network. Most previous research has focused on the contributions of short cycles to network dynamics. Here, we identify the specific cycles that are used during different runs of activation in sparse random graphs, as a basis of characterizing the contribution of cycles of any length. Both simulation and a refined mean-field approach evidence a decrease in the cycle usage when the cycle length increases, reflecting a trade-off between long time for recovery after excitation and low vulnerability to out-of-phase external excitations. In spite of this statistical observation, we find that the successful usage of long cycles, though rare, has important functional consequences for sustaining network activity: The average cycle length is the main feature of the cycle length distribution that affects the average lifetime of activity in the network. Particularly, use of long, rather than short, cycles correlates with higher lifetime, and cutting shortcuts in long cycles tends to increase the average lifetime of the activity. Our findings, thus, emphasize the essential, previously underrated role of long cycles in sustaining network activity. On a more general level, the findings underline the importance of network topology, particularly cycle structure, for self-sustained network dynamics.

Slow synaptic dynamics in a network: from exponential to power-law forgetting

We investigate a mean-field model of interacting synapses on a directed neural network. Our interest lies in the slow adaptive dynamics of synapses, which are driven by the fast dynamics of the neurons they connect. Cooperation is modeled from the usual Hebbian perspective, while competition is modeled by an original polarity-driven rule. The emergence of a critical manifold culminating in a tricritical point is crucially dependent on the presence of synaptic competition. This leads to a universal 1/t power-law relaxation of the mean synaptic strength along the critical manifold and an equally universal 1/√[t] relaxation at the tricritical point, to be contrasted with the exponential relaxation that is otherwise generic. In turn, this leads to the natural emergence of long- and short-term memory from different parts of parameter space in a synaptic network, which is the most original and important result of our present investigations.

Emergence of self-sustained oscillations in excitable Erdös-Rényi random networks

We investigate the emergence of self-sustained oscillations in excitable Erdös-Rényi random networks (EERRNs). Interestingly, periodical self-sustained oscillations have been found at a moderate connection probability P. For smaller or larger P, the system evolves into a homogeneous rest state with distinct mechanisms. One-dimensional Winfree loops are discovered as the sources to maintain the oscillations. Moreover, by analyzing these oscillation sources, we propose two criteria to explain the spatiotemporal dynamics obtained in EERRNs. Finally, the two critical connection probabilities for which self-sustained oscillations can emerge are approximately predicted based on these two criteria.

Analysis of epileptic seizures with complex network

Epilepsy is a disease of abnormal neural activities involving large area of brain networks. Until now the nature of functional brain network associated with epilepsy is still unclear. Recent researches indicate that the small world or scale-free attributes and the occurrence of highly clustered connection patterns could represent a general organizational principle in the human brain functional network. In this paper, we seek to find whether the small world or scale-free property of brain network is correlated with epilepsy seizure formation. A mass neural model was adopted to generate multiple channel EEG recordings based on regular, small world, random, and scale-free network models. Whether the connection patterns of cortical networks are directly associated with the epileptic seizures was investigated. The results showed that small world and scale-free cortical networks are highly correlated with the occurrence of epileptic seizures. In particular, the property of small world network is more significant during the epileptic seizures.