Abstract
The human brain contains billions of neurons that flexibly interconnect to support local and global computational spans. As neuronal activity propagates through the neural medium, it approaches a critical state hedged between ordered and disordered system regimes. Recent work demonstrates that this criticality coincides with the small-world topology, a network arrangement that accommodates both local (subcritical) and global (supercritical) system properties. On one hand, operating near criticality is thought to offer several neurocomputational advantages, e.g., high-dynamic range, efficient information capacity, and information transfer fidelity. On the other hand, aberrations from the critical state have been linked to diverse pathologies of the brain, such as post-traumatic epileptiform seizures and disorders of consciousness. Modulation of brain activity, through neuromodulation, presents an attractive mode of treatment to alleviate such neurological disorders, but a tractable neural framework is needed to facilitate clinical progress. Using a variation on the generative small-world model of Watts and Strogatz and Kuramoto's model of coupled oscillators, we show that the topological and dynamical properties of the small-world network are divided into two functional domains based on the range of connectivity, and that these domains play distinct roles in shaping the behavior of the critical state. We demonstrate that short-range network connections shape the dynamics of the system, e.g., its volatility and metastability, whereas long-range connections drive the system state, e.g., a seizure. Together, these findings lend support to combinatorial neuromodulation approaches that synergistically normalize the system dynamic while mobilizing the system state.
Introduction
The human brain is thought to contain billions of neurons that densely interconnect across short and long spatial distances (von Bartheld et al., 2016). The pattern of neuronal activity hinges on the anatomical and functional medium by which it is generated, and in which it propagates (Figure 1) (Wolfram, 1984a,b; Perc, 2007; Wang et al., 2010). In a hypothetical lattice, where nodes are highly ordered and hold no long-range shortcuts, signals tend to fizzle out locally due to the resistance that high-nodal separation exerts on global transmission (Shew and Plenz, 2012). This contrast with more disordered graphs where signals tend to overwhelm the global network span through dense interconnectivity. Intermediately, in the “small-world” network formed by integrating just a few long-range shortcuts into an otherwise ordered lattice (Watts and Strogatz, 1998), signals tend to reverberate, perched on the edge of chaos in a so-called “critical” state (Shew and Plenz, 2012; Kim and Lim, 2015). Intriguingly, it is thought that the functional topology of the brain tends to this criticality (Takagi, 2018), flexibly maneuvering it based on an immediate operational needs; by dynamically recruiting or abandoning short- and long-range functional connections, e.g., through coherence of neuroelectric oscillations (Singer, 1999; Buzsáki, 2006; Akam and Kullmann, 2014), or neuroplasticity (Dan and Poo, 2004; Shin and Kim, 2006), the brain maneuvers clustered and disordered topological phases tuned to local and global operational spans, respectively. Within this theoretical framework, the dynamics of the brain essentially reflect a dialectic on one hand pulling the brain to its topological extremes (Poil et al., 2012; Shew and Plenz, 2012; Hesse and Gross, 2014), while, on the other hand, keeping it near the critical state (Shin and Kim, 2006; Hesse and Gross, 2014; Priesemann, 2015; Takagi, 2018). Operating near criticality is thought to offer several neurocomputational advantages, e.g., high-dynamic range, efficient information capacity, and information transfer fidelity. In turn, aberrations from criticality have been theorized to underpin distinct neuropathologies, such as post-traumatic epilepsy and consciousness disorders (Colombo et al., 2016). Certainly, the malfunction of long- and short-range functional connections, by injury or otherwise, could have disastrous effects on the dynamics of the brain (Pevzner et al., 2016). In this Original Research article, we investigate specifically how long- and short-range connections affect the topological and dynamical properties of the small-world network. Our results indicate that short-range connections shape the dynamics of the system, whereas long-range connections define its state. We discuss the implications of these differential effects on clinical neuromodulation.
Figure 1
Methods
See Table 1 for model parameters.
Table 1
| Description | Notation | Notes |
|---|---|---|
| Network topology | ||
| Short-range connections | h, H | [4–50] |
| Long-range connections | g, G | [0.001–10] |
| Adjacency matrix | Amn | |
| Small-world coefficient | ω | [−1–+1] |
| Kuramoto-type simulation | ||
| Oscillatory phase | θm | [0°–360°] |
| Network synchrony | r | [0–1] |
| Natural frequency | ϵ | Gaussian: μ = 0, σ = 1 |
| Total nodes | N | 1,000 |
Model parameters.
Network Generation
To keep a manageable number of free parameters, and to reduce the artifacts of boundary conditions, we restricted our analysis to a generative ring network model based on the small-world model of Watts and Strogatz (1998). These ring networks were generated using custom Python code based on the open-source module networkx. Briefly, N = 1,000 nodes were each wired to their h nearest neighbors, thus denoted “short-range” connections (for h well below saturation, h < < N). Next, each node on average received an additional set of g random, yet unique, wires, which were denoted “long-range” connections (as wires of g did not equal those of h). Concretely, long-range connectivity g was generated by a nested loop given by g = T× u, where T is the maximum number of additional wires per node, and u is the average fraction of these actuated. To approximate biological gray-white matter ratios, while retaining a connected graph, we kept the range of short-range connectivity ~10 times that of long-range connections (Bajada et al., 2019; Mota et al., 2019).
Finally, each network was defined by its adjacency matrix, Amn, which was used for network simulation analysis (see Kuramoto's Model of Coupled Oscillators).
Small-World Coefficient
To quantify the extent to which a network resembled a small-world network, we computed the small-world coefficient ω (Telesford et al., 2011). Essentially, the small-world coefficient compares the resemblance of a network to a perfectly ordered vs. a perfectly disordered arrangement based on the extent to which the nodes of the network are clustered and the extent to which they are separated. The small-world coefficient is defined as:
where L is the average shortest path length between nodes in the network, and C is the degree of clustering (Figure 1). The disordered and ordered networks were generated based on the long-range connectivity given by g = T× u (see Network Generation). For the perfectly ordered network, no long-range connections were added, thus u = 0, and consequently g = 0. For the perfectly disordered network, the maximum number of long-range connections was introduced, thus u = 1, and consequently g = T.
The network parameters C and L were computed using common graph theory methods. Concretely, clustering C was computed as the network transitivity, such that:
where ∇ is the number of closed triplets in the network, and Tr is the maximum number of triplets. The average shortest path length L was given by:
where V is the set of nodes in the network, D(s,t) is the shortest path length from node s to t, and N is the total number of nodes. Thus, when network separation L ≈ Ldisordered, and network clustering C < < Cordered, the small-world coefficient ω ≈ +1, meaning that the network approximates a perfectly disordered graph. Similarly, for the perfectly ordered lattice, when network clustering C ≈ Cordered and network separation L >> Ldisordered, the small-world coefficient approximates ω ≈ −1. Crucially, the small-world topology is defined as the critical state possessing both qualities, namely, network clustering similar to an ordered lattice C ≈ Cordered,and network separation similar to a disordered graph L ≈ Ldisordered; thus, the small-world coefficient tends to ω ≈ 0 as the network tends to the critical small-world arrangement.
Kuramoto's Model of Coupled Oscillators
Each node of the network was modeled as a coupled Kuramoto-type oscillator (Yamamoto et al., 2018), described by the set of N-coupled differential equations (Breakspear et al., 2010):
where the nth oscillator with a natural frequency εn adjusts its phase velocity based on the pair-wise phase interactions with its coupled peers (provided by the adjacency matrix Amn, see Network Generation). The internodal coupling was K = 3, and the natural frequencies were distributed according to the Gaussian probability density with mean ε0 = 0. The state of the node (n = 1, …, N) was thus defined by its phase θ, which was calculated by the Livermore Solver for Ordinary Differential Equations (LSODA) method with a dynamic time step.
The degree of synchrony in the network was quantified by the order parameter r, given by:
where ψ is the mean phase of the set of oscillators N, and the scalar r represents the order, or phase uniformity, of the network. An open-source Python implementation of the Kuramoto oscillatory system is available online, which was used to generate the simulation data presented here (Damicelli, 2021).
Stability and Attractiveness Analysis
To compute the stability of different network states, we set the initial synchrony level of the network via the initial nodal phases θm0. Thus, for initial synchrony r0 = 0.5, on average 50% of the nodes of network had equal phases in the initial state. Then, at a predefined time-step Δt = 250, not necessarily in the steady-state, the deviation of the network from the initial synchrony level was computed, revealing the stability of the initial state. Specifically, larger deviations reflect weaker stability. Repeating this process for all combinations of initial synchrony levels and connectivity parameters produces the stability heat maps depicted in Figure 4A.
To calculate the attractiveness of different network states, we checked which synchrony levels the networks shifted to during the simulations and scored the end synchrony level based on the size of the shift. For instance, for initial synchrony r0 = 0.5, and long-range connectivity g = 0.001, one network might end up in end synchrony level r1 = 0.0. This adds a score of s = |r0 – r1| = 0.5 to the end synchrony state r1. The synchrony states holding the highest cumulative scores had the highest attractiveness. Repeating this process for all the combinations of initial synchrony levels and connectivity parameters produces the attractiveness heat maps given in Figure 4B.
Results
To examine the effects of long- and short-range connections on the topological and dynamical properties of the small-world network, we applied a variation on small-world model of Watts and Strogatz (1998) and the Kuramoto model of coupled oscillators (Kuramoto, 1984). Concretely, we generated an ordered ring lattice consisting of N = 1,000 nodes, each node connected to its h nearest neighbors. To this base, we added, on average, g unique long-range connections per node. Thus, by definition, long-range connections g were topologically distinct from their short-range correlates h, holding true for short-range connectivity well below network saturation, h < < N. Then, to quantify the extent to which a network resembled a small-world network, we computed the small-world coefficient ω (Telesford et al., 2011). Essentially, the small-world coefficient compares the resemblance of a network to a perfectly ordered vs. a perfectly disordered arrangement based on the extent to which the nodes of the network are clustered and the extent to which they are separated. More specifically, ordered, subcritical lattices tend to ω ≈ –1, having high-clustering and high-separation parameters; disordered, supercritical graphs tend to ω ≈ +1, having low-clustering and low-separation parameters; and critical small-world topologies tend to ω ≈ 0, having both the ordered and disordered tendencies balanced out (see Methods). Within this definition, we visualized the topological behavior of the network by plotting the small-world coefficient ω as a function of the long- and short-range connectivity g and h, respectively (Figure 2). Finally, we used the same topological framework to generate networks of coupled Kuramoto-type oscillators (Kuramoto, 1984).
Figure 2
Long-Range Connections Dominate the Topological State
We first examined the roles of short- and long-range connections in defining the topological state of the network, specifically by keeping one parameter static (uppercase letters H and G) while modulating the other (lowercase letters h and g) (Figure 2). We found that modulations of the long-range connectivity g offered a near full topological range despite the underlying static short-range connectivity H (Figure 2A; ~70% ± 0.05; mean ± SEM). Yet, the opposite was not the case: Invariant to the underlying long-range connectivity G, increases to the short-range connectivity h all converged to the critical state (Figures 2A,B). In general, less than half of the topological range was attainable by short-range modulation alone. Thus, short-range connections appear to be poorly suited as a modulator of the topological state.
For further examination, we computed the first derivative of the topological state to reveal the state mobility Δω of the network, i.e., how readily the network moved from one topological state to another via changes to its connectivity parameters h and g (Figure 2C). We found that across all underlying long-range connectivities G, modulation of the short-range connectivity h had near null effects on the topological state. Modulation of the long-range connectivity g of the network, however, offered potent state mobilization within the subcritical and critical regimes, but near null mobility approaching supercriticality. The dominant role of long-range connections on the topological state was confirmed by dominance analysis (R2 ~ 0.648 for long-range vs. R2 ~ 0.003 for short-range connections).
These results together indicate that the topological state of the small-world network is dominantly defined by the long-range connectivity (Watts and Strogatz, 1998) and that the topological mobility of the network is the most potent well below supercriticality (Carhart-Harris et al., 2014).
Short-Range Connections Shape the Topological Dynamics
Next, we evaluated how the underlying short-range connectivity H affects the topological behavior of the network, as reflected by the shape of the topological state curves (Figure 2). We found that as the static short-range connectivity H was reduced, the state curve steepened about the critical point, thus, contracting and “right-shifting” the critical regime to higher values of the long-range connectivity g (Figure 2A). This indicates that, to sustain the small-world criticality, poorly clustered networks (low H) must integrate long-range connections to a greater extent, yet within a narrower limit.
We then calculated the difference in state mobility of networks that had a high-static short-range connectivity (H = 100) and a low-static short-range connectivity (H = 10) (Figure 2C). In this difference plot, negative values reflect a reduction in the state mobility of the network, which essentially equates to a stabilization of the topological state (and oppositely for the positive values). Intriguingly, we found that, as the static short-range connectivity H was reduced, the stability of the topological state shifted to the subcritical regime, strongly destabilizing the small-world criticality (Figure 2D). This indicates that the short-range connectivity of the network has fundamental effects on the stability of the network across diverse topological regimes.
Network Synchronizability
To extend our topological findings, we examined the synchronization properties offered by small-world networks of varying short- and long-range connectivity parameters (Figure 3). To this end, we quantified the global network synchrony at the steady-state using Kuramoto's order parameter r, which reflects that the global phase uniformity of the network nodes (see Methods).
Figure 3
By gradually integrating long-range connections into the network structure, our simulations show that the synchronizability abruptly reaches a critical point at which the network shifts from a state of low synchrony to near-complete synchrony (Figure 3A). Such “explosive synchronization,” a critical transitioning, is characteristic for the Kuramoto-type coupled oscillators (Kuramoto, 1984; Gómez-Gardeñes et al., 2011; Boccaletti et al., 2016), and mirrors the topological criticality of small-world networks (Figure 1) (Watts and Strogatz, 1998).
Next, we modeled the synchronizability using four-parameter logistic regression (Figure 3B). Like in our topological findings, we found that the slope of the critical transitioning b steepened as the static short-range connectivity H was reduced, indicating a destabilization and narrowing of the critical regime (Figure 3B, inset). The minimal synchronizability a of the network moreover related proportionally to the static short-range connectivity, indicating baseline synchronization hinged separately on short-range interactions. Indeed, as more short-range connections are introduced, the network ultimately reaches a point of saturation where global synchronization becomes deterministic, invariant to topological modulations (Barahona and Pecora, 2002).
Finally, by calculating the divergence of the network from a predefined initial synchrony level (see Methods), we examined the stability and attractiveness of different network states (Figure 4). First, these data confirm that the stability of the critical regime narrows as short-range connectivity decreases. Second, as short-range connections are removed, we find that the network becomes increasingly attracted to subcritical synchrony states (Figure 4B), which aligns with the topological destabilization that favors subcriticality presented earlier (cf. red curves in Figures 2D, 4D, and inset).
Figure 4
Network Metastability Depends on Short-Range Connections
Our data show that the functional interactions of the network converge as the static short-range connectivity decreases (Figure 3B, note narrowing standard deviation). Accordingly, we found that the long-range connectivity g in the poorly clustered network (low H) had a very high predictive power (PPS) on the global synchrony r of the network, whereas highly clustered networks (high H) were generally poorly predictable (H = 10, PPS = 0.93; H = 100, PPS = 0.43) (Wetschoreck et al., 2020). Furthermore, as the short-range connectivity of the networks tends to saturation (H to N = 1,000), the PPS drops to 0. More specifically, we find that the PPS is linearly proportional to the short-range connectivity of the network (PPS = −0.010H + 0.985;R2 = 0.99) (Supplementary Figure 1).
The PPS can be used to assess the metastability of the network. Thought to be inherent to cognition (Alderson et al., 2020), metastability defines a dynamical regime that accommodates flexible interactions of network nodes without stagnating in the fixed positions (Hellyer et al., 2015). Thus, our results show that metastability of the network depends linearly on the underlying short-range connectivity (by dominance analysis, R2 = 0.438 for the short-range vs. R2 = 0.136 for long-range connections). These simulation data altogether mirror our topological findings by suggesting that short-range connections are pivotal for the network's system dynamics (Figure 4).
Discussion
We have investigated the effects of short- and long-range connections on the topological and dynamical properties of the small-world network. Converging with previous work (Watts and Strogatz, 1998), we demonstrate, first, that long-range connections determine the topological and functional state of the network. Second, we show that short-range connections shape the dynamics of the system, i.e., the stability of the system across diverse topological regimes (Figures 2, 4). Our findings together provide evidence that short- and long-range connections play distinct roles in shaping the behavior of the small-world network.
The topological properties of a network have fundamental effects on the activity taking place on it (Strogatz, 2001). Several works have, for instance, analyzed the spread of infectious disease in small-world networks, finding fluctuations between sporadic endemic and self-sustaining epidemic infectious cycles based on network disorder (Kuperman and Abramson, 2001; Rüdiger et al., 2020). Others have examined the synchronizability of coupled oscillators on small-world graphs (Barahona and Pecora, 2002; Nishikawa et al., 2003). Later, such simulations have been expanded to examine cortical oscillations and neuroplasticity (Maistrenko et al., 2007; Breakspear et al., 2010).
The human brain is a complex system sustained by the interactions of billions of neurons across local and global spatial scales. Previous work has shown that the functional topology of the brain tends to a small-world-like criticality that accommodates both local (subcritical) and global (supercritical) system properties (Bassett and Bullmore, 2017; Takagi, 2018). The hypothesis that the brain maintains a proximity to the critical state stems from the premise of superior computational adaptability to rapidly changing operational demands (Massobrio et al., 2015a). Contention posits, however, that signatures for criticality, e.g., power-law distributions, could be artifacts of sampling (Touboul and Destexhe, 2010; Marsili et al., 2013), multiplicative noise (Sornette, 1998) or emerge from “hidden variables” not necessarily linked to network topology (Aitchison et al., 2016; Morrell et al., 2021). While an exhaustive review is beyond the scope of this discussion (Beggs and Timme, 2012), we note that diverse data supports the relationship between critical neural dynamics and small-world topologies (Massobrio et al., 2015b; Tan and Cheong, 2017; Takagi, 2018) and the presence of critical signatures in human fMRI (Kitzbichler et al., 2009), local field potentials (Petermann et al., 2009), spike data (Friedman et al., 2012), human brain oscillations (Poil et al., 2008), and artificial neural networks (Shin and Kim, 2006). Indeed, congruent with a near-critical regime (Priesemann, 2015), the brain operates within a wide dynamic range that accommodate high-level cognition through global neural coordination (Taylor et al., 2015), and low-activity states, such as anesthesia (Brown et al., 2010), and, to some extent, sleeping (Priesemann et al., 2013; Tagliazucchi and van Someren, 2017), marked by weaker, more fragmented interactions outside the local milieu.
It is believed that neural oscillations, or “brain waves,” mediate short- and long-range neural connectivity through high- and low-frequency wavebands, respectively (Kopell et al., 2000; Buzsáki, 2006; Tiesinga and Sejnowski, 2009). In essence, the wave interference of oscillating neural populations facilitates the selective transfer of information (Singer, 1999; Buzsáki, 2006; Akam and Kullmann, 2010). Thus, it has been hypothesized that the malfunction of such neural interactions may have deleterious effects on the brain's system dynamics (Uhlhaas and Singer, 2006; Pevzner et al., 2016). In agreement with this premise, our results indicate that impairments to the network's short-range connectivity destabilize the small-world criticality in favor of extreme network regimes, i.e., sub- and supercriticality (Figures 2D, 4C). Such departure from criticality has been linked to large-scale fMRI signatures of unconsciousness (Tagliazucchi et al., 2016).
Subcritical networks tend to be states of desynchronization and clustering that perturb global network processing, e.g., cognition (Roozenbeek et al., 2013). Congruently, our simulations show that sparsely clustered networks, with poor short-range connectivity, exhibit weak metastability (Supplementary Figure 1), which has been correlated with deficits in cognitive flexibility (Hellyer et al., 2015).
Notably, our stability analysis indicates that damage to the short-range connectivity of the network could produce a “repellant peak” that effectively barricades the critical regime, trapping the network activity in a subcritical trough (Figure 4D). Such “subcritical entrapment” aligns with the behavioral heterogeneity of persistent disorders of consciousness (Giacino et al., 2014), e.g., partial retainment of cognitive processing, and lends theoretical support to the rehabilitation of the system dynamic, e.g., through short-range neural potentiation.
The supercritical network, on the other hand, tends to hypersynchrony, broadly resembling the state of seizures (Szaflarski et al., 2014; Zimmern, 2020). Indeed, researchers have argued that epileptiform seizures reflect a critical–supercritical transition (Arviv et al., 2016; Bauer et al., 2017; Freestone et al., 2017), which was recently supported by a strong electroencephalographic sign in human patients (Scheffer et al., 2009; Maturana et al., 2020). Similarly, Gerster et al. report that artificial neuronal oscillators on supercritical small-world graphs mirror electroencephalographic epileptic patterns (Gerster et al., 2020). The refractoriness of some types of epilepsy could thus reflect an underlying destabilization of the critical regime by elimination of the short-range connections, such as through cortical dysgenesis or brain trauma (Semah et al., 1998). Interestingly, recent work on the Kuramoto model has shown that generalized resource constraints seed the network to self-terminating supercritical episodes (Frolov and Hramov, 2021), consistent with epileptic recurrences.
One potential mechanism for the disruption of short-range neural connectivity may be an injury to key brain hubs that contain a high-cumulative weight of short-range connections (Gratton et al., 2012; Zhou et al., 2012; Haimovici et al., 2016; Yuan et al., 2017). Indeed, hubs, e.g., the cingulate cortex, have been shown to be instrumental for cognitive performance (Fagerholm et al., 2015; Li et al., 2019), and have profound effects on the functional connectivity of simulated networks (Aerts et al., 2016). It is interesting to note that compensation to injury could thus predictably be offered by the recruitment, or hyperactivity, of dense hub regions, which has been widely hypothesized (Hillary et al., 2011, 2015; Tang et al., 2012; Iraji et al., 2016), e.g., in components of the default mode network (Zhou et al., 2012).
Our findings altogether lend support to combinatorial neuromodulation strategies that target short- and long-range neural connectivity differentially, to normalize the system dynamic and mobilize the system state, respectively. Future work will target components of short- and long-range neural communication, e.g., through pharmacological neurostimulation via amantadine to preferentially enhance low-frequency brain oscillations (Ott et al., 2018; Ma and Zafonte, 2020), direct current stimulation of deep brain structures, e.g., hippocampal theta (Lee et al., 2013), or modulation of cerebral cortex gamma (Pink et al., 2019), e.g., using cell-type-specific optogenetic or pharmacogenetic modulation (Liu et al., 2020), or non-invasive transcranial magnetic stimulation at low frequencies (Farzan et al., 2012).
Study Limitations
There are several limitations to this study. First, while providing a useful conceptual framework, Watts and Strogatz's ring model does not reflect real brain connectivity known to contain non-random edge distributions, e.g., “rich hubs” (van den Heuvel and Sporns, 2011), and a scale-free degree distribution (Eguíluz et al., 2005). Still, reduced topologies, e.g., generative small-worlds (Netoff et al., 2004; Perc, 2007; Tekin and Tagluk, 2017), and randomized graphs (van Vreeswijk and Sompolinsky, 1996; Tsodyks et al., 2000) remain valuable to neural network analysis by offering a controlled computational environment with manageable parameters and optimized network conditions.
Second, Kuramoto's oscillatory model represents a reduction of the complex interactions of distributed neural populations (Singer, 1999; Buzsáki, 2006). It is plausible that fuller physiological models would provide deeper insights into the precise mechanisms of such neural interactions. In support of the applicability of Kuramoto's equations, however, simulations have previously been applied to macaque (Honey and Sporns, 2008), and human brain research (Kitzbichler et al., 2009; Cabral et al., 2014), showing high congruence between simulation data and resting-state activity (Cabral et al., 2014; Vuksanović and Hövel, 2014). More broadly, reduced models (Siettos and Starke, 2016), such as two-state units (van Vreeswijk and Sompolinsky, 1996), and the FitzHugh–Nagumo model (Perc, 2007; Gerster et al., 2020), have been used extensively to examine complex network behaviors, such as self-organized balanced states (van Vreeswijk and Sompolinsky, 1996). Similarly, the abstraction offered by Kuramoto's model allows tractable simulations and analyses, holding high value for the investigation of more fundamental principles of oscillatory dynamics (Breakspear et al., 2010), such as the functional division of network connectivity examined here.
Despite these limitations, this study provides important insights into the relationship between network connectivity and critical system dynamics, which are broadly consistent with empirical reports and previous work (Haimovici et al., 2016). Future research should apply brain connectomic data and fuller network simulations to extend these findings.
Funding
We acknowledge the following grants: Lundbeck Foundation (DANDRITE-R248-2016-2518; R344-2020-300; and R351-2020-1095), Novo Nordisk Foundation (NNF20OC0064395), and European Research Council Starting (638730) grants to KY.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The datatsets presented in this study can be found in the online repository: https://github.com/simonarvin/connectivity_smallworld.
Author contributions
SA and AG conceived the project. SA designed the project, performed the computations, and analyzed the data. SA, AG, and KY interpreted the data and wrote the manuscript. All authors contributed to the article and approved the submitted version.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fncom.2021.783474/full#supplementary-material
Supplementary Figure 1Predictive power of long-range connectivity g on network synchrony r. The predictive power of long-range connections on the synchrony state of the network is linearly proportional to short-range connectivity of the network. In fact, as short-range connectivity H tends to saturation N = 1,000, the predictive power tends to 0. g, long-range connectivity; r, network synchrony; H, short-range connectivity; N, total nodes; PPS, predictive power.
References
1
AertsH.FiasW.CaeyenberghsK.MarinazzoD. (2016). Brain networks under attack: robustness properties and the impact of lesions. Brain139, 3063–3083. 10.1093/brain/aww194
2
AitchisonL.CorradiN.LathamP. E. (2016). Zipf's law arises naturally when there are underlying, unobserved variables. PLoS Comput. Biol.12, e1005110. 10.1371/journal.pcbi.1005110
3
AkamT.KullmannD. M. (2010). Oscillations and filtering networks support flexible routing of information. Neuron67, 308–320. 10.1016/j.neuron.2010.06.019
4
AkamT.KullmannD. M. (2014). Oscillatory multiplexing of population codes for selective communication in the mammalian brain. Nat. Rev. Neurosci.111–122. 10.1038/nrn3668
5
AldersonT. H.BokdeA. L. W.KelsoJ. A. S.MaguireL.CoyleD. (2020). Metastable neural dynamics underlies cognitive performance across multiple behavioural paradigms. Hum. Brain Mapp.41, 3212–3234. 10.1002/hbm.25009
6
ArvivO.MedvedovskyM.SheintuchL.GoldsteinA.ShrikiO. (2016). Deviations from critical dynamics in interictal epileptiform activity. J. Neurosci.36, 12276–12292. 10.1523/JNEUROSCI.0809-16.2016
7
BajadaC. J.SchreiberJ.CaspersS. (2019). Fiber length profiling: a novel approach to structural brain organization. Neuroimage186, 164–173. 10.1016/j.neuroimage.2018.10.070
8
BarahonaM.PecoraL. M. (2002). Synchronization in small-world systems. Phys. Rev. Lett.89, 054101. 10.1103/PhysRevLett.89.054101
9
BassettD. S.BullmoreE. T. (2017). Small-world brain networks revisited. Neuroscientist23, 499–516. 10.1177/1073858416667720
10
BauerP. R.ThijsR. D.LambertsR. J.VelisD. N.VisserG. H.TolnerE. A.et al. (2017). Dynamics of convulsive seizure termination and postictal generalized EEG suppression. Brain140, 655–668. 10.1093/brain/aww322
11
BeggsJ. M.TimmeN. (2012). Being critical of criticality in the brain. Front. Physiol.3, 163. 10.3389/fphys.2012.00163
12
BoccalettiS.AlmendralJ. A.GuanS.LeyvaI.LiuZ.Sendiña-NadalI.et al. (2016). Explosive transitions in complex networks' structure and dynamics: percolation and synchronization. Phys. Rep.660, 1–94. 10.1016/j.physrep.2016.10.004
13
BreakspearM.HeitmannS.DaffertshoferA. (2010). Generative models of cortical oscillations: neurobiological implications of the Kuramoto model. Front. Hum. Neurosci.4, 190. 10.3389/fnhum.2010.00190
14
BrownE. N.LydicR.SchiffN. D. (2010). General anesthesia, sleep, and coma. N. Engl. J. Med.363, 2638–2650. 10.1056/NEJMra0808281
15
BuzsákiG. (2006). Rhythms of the Brain. Oxford, UK; New York: Oxford University Press.
16
CabralJ.KringelbachM. L.DecoG. (2014). Exploring the network dynamics underlying brain activity during rest. Prog. Neurobiol.114, 102–131. 10.1016/j.pneurobio.2013.12.005
17
Carhart-HarrisR. L.LeechR.HellyerP. J.ShanahanM.FeildingA.TagliazucchiE.et al. (2014). The entropic brain: a theory of conscious states informed by neuroimaging research with psychedelic drugs. Front. Hum. Neurosci.8, 20. 10.3389/fnhum.2014.00020
18
ColomboM. A.WeiY.RamautarJ. R.Linkenkaer-HansenK.TagliazucchiE.Van SomerenE. J. W. (2016). More severe insomnia complaints in people with stronger long-range temporal correlations in wake resting-state EEG. Front. Physiol.7, 576. 10.3389/fphys.2016.00576
19
DamicelliF. (2021). kuramoto: Python implementation of the Kuramoto model on graphs. Github. Available online at: https://github.com/fabridamicelli/kuramoto (accessed August 4, 2021).
20
DanY.PooM.-M. (2004). Spike timing-dependent plasticity of neural circuits. Neuron44, 23–30. 10.1016/j.neuron.2004.09.007
21
EguíluzV. M.ChialvoD. R.CecchiG. A.BalikiM.ApkarianA. V. (2005). Scale-free brain functional networks. Phys. Rev. Lett.94, 018102. 10.1103/PhysRevLett.94.018102
22
FagerholmE. D.HellyerP. J.ScottG.LeechR.SharpD. J. (2015). Disconnection of network hubs and cognitive impairment after traumatic brain injury. Brain138, 1696–1709. 10.1093/brain/awv075
23
FarzanF.BarrM. S.SunY.FitzgeraldP. B.DaskalakisZ. J. (2012). Transcranial magnetic stimulation on the modulation of gamma oscillations in schizophrenia. Ann. N. Y. Acad. Sci.1265, 25–35. 10.1111/j.1749-6632.2012.06543.x
24
FreestoneD. R.KarolyP. J.CookM. J. (2017). A forward-looking review of seizure prediction. Curr. Opin. Neurol.30, 167–173. 10.1097/WCO.0000000000000429
25
FriedmanN.ItoS.BrinkmanB. A. W.ShimonoM.DeVilleR. E. L.DahmenK. A.et al. (2012). Universal critical dynamics in high resolution neuronal avalanche data. Phys. Rev. Lett.108, 208102. 10.1103/PhysRevLett.108.208102
26
FrolovN.HramovA. (2021). Extreme synchronization events in a Kuramoto model: the interplay between resource constraints and explosive transitions. Chaos31, 063103. 10.1063/5.0055156
27
GersterM.BernerR.SawickiJ.ZakharovaA.ŠkochA.HlinkaJ.et al. (2020). FitzHugh-Nagumo oscillators on complex networks mimic epileptic-seizure-related synchronization phenomena. Chaos30, 123130. 10.1063/5.0021420
28
GiacinoJ. T.FinsJ. J.LaureysS.SchiffN. D. (2014). Disorders of consciousness after acquired brain injury: the state of the science. Nat. Rev. Neurol.10, 99–114. 10.1038/nrneurol.2013.279
29
Gómez-GardeñesJ.GómezS.ArenasA.MorenoY. (2011). Explosive synchronization transitions in scale-free networks. Phys. Rev. Lett.106, 128701. 10.1103/PhysRevLett.106.128701
30
GrattonC.NomuraE. M.PérezF.D'EspositoM. (2012). Focal brain lesions to critical locations cause widespread disruption of the modular organization of the brain. J. Cogn. Neurosci.24, 1275–1285. 10.1162/jocn_a_00222
31
HaimoviciA.BalenzuelaP.TagliazucchiE. (2016). Dynamical signatures of structural connectivity damage to a model of the brain posed at criticality. Brain Connect.6, 759–771. 10.1089/brain.2016.0455
32
HellyerP. J.ScottG.ShanahanM.SharpD. J.LeechR. (2015). Cognitive flexibility through metastable neural dynamics is disrupted by damage to the structural connectome. J. Neurosci.35, 9050–9063. 10.1523/JNEUROSCI.4648-14.2015
33
HesseJ.GrossT. (2014). Self-organized criticality as a fundamental property of neural systems. Front. Syst. Neurosci.8, 166. 10.3389/fnsys.2014.00166
34
HillaryF. G.RomanC. A.VenkatesanU.RajtmajerS. M.BajoR.CastellanosN. D. (2015). Hyperconnectivity is a fundamental response to neurological disruption. Neuropsychology29, 59–75. 10.1037/neu0000110
35
HillaryF. G.SlocombJ.HillsE. C.FitzpatrickN. M.MedagliaJ. D.WangJ.et al. (2011). Changes in resting connectivity during recovery from severe traumatic brain injury. Int. J. Psychophysiol.82, 115–123. 10.1016/j.ijpsycho.2011.03.011
36
HoneyC. J.SpornsO. (2008). Dynamical consequences of lesions in cortical networks. Hum. Brain Mapp.29, 802–809. 10.1002/hbm.20579
37
IrajiA.ChenH.WisemanN.WelchR. D.O'NeilB. J.HaackeE. M.et al. (2016). Compensation through functional hyperconnectivity: a longitudinal connectome assessment of mild traumatic brain injury. Neural Plast.2016, 4072402. 10.1155/2016/4072402
38
KimS.-Y.LimW. (2015). Effect of small-world connectivity on fast sparsely synchronized cortical rhythms. Phys. A Stat. Mech. Appl.421, 109–123. 10.1016/j.physa.2014.10.019
39
KitzbichlerM. G.SmithM. L.ChristensenS. R.BullmoreE. (2009). Broadband criticality of human brain network synchronization. PLoS Comput. Biol.5, e1000314. 10.1371/journal.pcbi.1000314
40
KopellN.ErmentroutG. B.WhittingtonM. A.TraubR. D. (2000). Gamma rhythms and beta rhythms have different synchronization properties. Proc. Nat. Acad. Sci. U.S.A.97, 1867. 10.1073/pnas.97.4.1867
41
KupermanM.AbramsonG. (2001). Small world effect in an epidemiological model. Phys. Rev. Lett.86, 2909–2912. 10.1103/PhysRevLett.86.2909
42
KuramotoY. (1984). Cooperative dynamics of oscillator community: a study based on lattice of rings. Progr. Theoret. Phys.79, 223–240. 10.1143/PTPS.79.223
43
LeeD. J.GurkoffG. G.IzadiA.BermanR. F.EkstromA. D.MuizelaarJ. P.et al. (2013). Medial septal nucleus theta frequency deep brain stimulation improves spatial working memory after traumatic brain injury. J. Neurotrauma30, 131–139. 10.1089/neu.2012.2646
44
LiF.LuL.ChenH.WangP.ChenY.-C.ZhangH.et al. (2019). Disrupted brain functional hub and causal connectivity in acute mild traumatic brain injury. Aging11, 10684–10696. 10.18632/aging.102484
45
LiuL.XuH.WangJ.LiJ.TianY.ZhengJ.et al. (2020). Cell type-differential modulation of prefrontal cortical GABAergic interneurons on low gamma rhythm and social interaction. Sci. Adv.6, eaay4073. 10.1126/sciadv.aay4073
46
MaH. M.ZafonteR. D. (2020). Amantadine and memantine: a comprehensive review for acquired brain injury. Brain Inj.34, 299–315. 10.1080/02699052.2020.1723697
47
MaistrenkoY. L.LysyanskyB.HauptmannC.BurylkoO.TassP. A. (2007). Multistability in the Kuramoto model with synaptic plasticity. Phys. Rev. E Stat. Nonlin. Soft Matter Phys.75, 066207. 10.1103/PhysRevE.75.066207
48
MarsiliM.MastromatteoI.RoudiY. (2013). On sampling and modeling complex systems. arXiv arXiv:1301.3622. 10.1088/1742-5468/2013/09/P09003
49
MassobrioP.de ArcangelisL.PasqualeV.JensenH. J.PlenzD. (2015a). Criticality as a signature of healthy neural systems. Front. Syst. Neurosci.9, 22. 10.3389/fnsys.2015.00022
50
MassobrioP.PasqualeV.MartinoiaS. (2015b). Self-organized criticality in cortical assemblies occurs in concurrent scale-free and small-world networks. Sci. Rep.5, 10578. 10.1038/srep10578
51
MaturanaM. I.MeiselC.DellK.KarolyP. J.D'SouzaW.GraydenD. B.et al. (2020). Critical slowing down as a biomarker for seizure susceptibility. Nat. Commun.11, 2172. 10.1038/s41467-020-15908-3
52
MorrellM. C.SederbergA. J.NemenmanI. (2021). Latent dynamical variables produce signatures of spatiotemporal criticality in large biological systems. Phys. Rev. Lett.126, 118302. 10.1103/PhysRevLett.126.118302
53
MotaB.Dos SantosS. E.Ventura-AntunesL.Jardim-MessederD.NevesK.KazuR. S.et al. (2019). White matter volume and white/gray matter ratio in mammalian species as a consequence of the universal scaling of cortical folding. Proc. Natl. Acad. Sci. USA.116, 15253–15261. 10.1073/pnas.1716956116
54
NetoffT. I.ClewleyR.ArnoS.KeckT.WhiteJ. A. (2004). Epilepsy in small-world networks. J. Neurosci.24, 8075–8083. 10.1523/JNEUROSCI.1509-04.2004
55
NishikawaT.MotterA. E.LaiY.-C.HoppensteadtF. C. (2003). Heterogeneity in oscillator networks: are smaller worlds easier to synchronize?Phys. Rev. Lett.91, 014101. 10.1103/PhysRevLett.91.014101
56
OttT.WestendorffS.NiederA. (2018). Dopamine receptors influence internally generated oscillations during rule processing in primate prefrontal cortex. J. Cogn. Neurosci.30, 770–784. 10.1162/jocn_a_01248
57
PercM. (2007). Effects of small-world connectivity on noise-induced temporal and spatial order in neural media. Chaos Solit. Fract.31, 280–291. 10.1016/j.chaos.2005.10.018
58
PetermannT.ThiagarajanT. C.LebedevM. A.NicolelisM. A. L.ChialvoD. R.PlenzD. (2009). Spontaneous cortical activity in awake monkeys composed of neuronal avalanches. Proc. Natl. Acad. Sci. U.S.A.106, 15921–15926. 10.1073/pnas.0904089106
59
PevznerA.IzadiA.LeeD. J.ShahlaieK.GurkoffG. G. (2016). Making waves in the brain: what are oscillations, and why modulating them makes sense for brain injury. Front. Syst. Neurosci.10, 30–30. 10.3389/fnsys.2016.00030
60
PinkA. E.WilliamsC.AldermanN.StoffelsM. (2019). The use of repetitive transcranial magnetic stimulation (rTMS) following traumatic brain injury (TBI): a scoping review. Neuropsychol. Rehabil. 31, 479–505. 10.1080/09602011.2019.1706585
61
PoilS.-S.HardstoneR.MansvelderH. D.Linkenkaer-HansenK. (2012). Critical-state dynamics of avalanches and oscillations jointly emerge from balanced excitation/inhibition in neuronal networks. J. Neurosci.32, 9817–9823. 10.1523/JNEUROSCI.5990-11.2012
62
PoilS.-S.van OoyenA.Linkenkaer-HansenK. (2008). Avalanche dynamics of human brain oscillations: relation to critical branching processes and temporal correlations. Hum. Brain Mapp.29, 770–777. 10.1002/hbm.20590
63
PriesemannV. (2015). Self-organization to sub-criticality. BMC Neurosci.16, O19. 10.1186/1471-2202-16-S1-O19
64
PriesemannV.ValderramaM.WibralM.Le Van QuyenM. (2013). Neuronal avalanches differ from wakefulness to deep sleep–evidence from intracranial depth recordings in humans. PLoS Comput. Biol.9, e1002985. 10.1371/journal.pcbi.1002985
65
RoozenbeekB.MaasA. I. R.MenonD. K. (2013). Changing patterns in the epidemiology of traumatic brain injury. Nat. Rev. Neurol.9, 231–236. 10.1038/nrneurol.2013.22
66
RüdigerS.PlietzschA.SaguésF.SokolovI. M.KurthsJ. (2020). Epidemics with mutating infectivity on small-world networks. Sci. Rep.10, 5919. 10.1038/s41598-020-62597-5
67
SchefferM.BascompteJ.BrockW. A.BrovkinV.CarpenterS. R.DakosV.et al. (2009). Early-warning signals for critical transitions. Nature461, 53–59. 10.1038/nature08227
68
SemahF.PicotM. C.AdamC.BroglinD.ArzimanoglouA.BazinB.et al. (1998). Is the underlying cause of epilepsy a major prognostic factor for recurrence?Neurology51, 1256–1262. 10.1212/WNL.51.5.1256
69
ShewW. L.PlenzD. (2012). The functional benefits of criticality in the cortex. Neuroscientist19, 88–100. 10.1177/1073858412445487
70
ShinC.-W.KimS. (2006). Self-organized criticality and scale-free properties in emergent functional neural networks. Phys. Rev. E Stat. Nonlin. Soft Matter Phys.74, 045101. 10.1103/PhysRevE.74.045101
71
SiettosC.StarkeJ. (2016). Multiscale modeling of brain dynamics: from single neurons and networks to mathematical tools. Wiley Interdiscip. Rev. Syst. Biol. Med.8, 438–458. 10.1002/wsbm.1348
72
SingerW. (1999). Neuronal synchrony: a versatile code for the definition of relations?Neuron24, 49–65, 111–25. 10.1016/S0896-6273(00)80821-1
73
SornetteD. (1998). Multiplicative processes and power laws. Phys. Rev. E57, 4811–4813. 10.1103/PhysRevE.57.4811
74
StrogatzS. H. (2001). Exploring complex networks. Nature410, 268–276. 10.1038/35065725
75
SzaflarskiJ. P.NazzalY.DreerL. E. (2014). Post-traumatic epilepsy: current and emerging treatment options. Neuropsychiatr. Dis. Treat.10, 1469–1477. 10.2147/NDT.S50421
76
TagliazucchiE.ChialvoD. R.SiniatchkinM.AmicoE.BrichantJ.-F.BonhommeV.et al. (2016). Large-scale signatures of unconsciousness are consistent with a departure from critical dynamics. J. R. Soc. Interface13, 20151027. 10.1098/rsif.2015.1027
77
TagliazucchiE.van SomerenE. J. W. (2017). The large-scale functional connectivity correlates of consciousness and arousal during the healthy and pathological human sleep cycle. Neuroimage160, 55–72. 10.1016/j.neuroimage.2017.06.026
78
TakagiK. (2018). Information-based principle induces small-world topology and self-organized criticality in a large scale brain network. Front. Comput. Neurosci.12, 65. 10.3389/fncom.2018.00065
79
TanT. L.CheongS. A. (2017). Statistical complexity is maximized in a small-world brain. PLoS ONE12, e0183918. 10.1371/journal.pone.0183918
80
TangC. Y.EavesE.Dams-O'ConnorK.HoL.LeungE.WongE.et al. (2012). Diffuse disconnectivity in tBi: a resting state fMri anD Dti stuDy. Transl. Neurosci.3, 9–14. 10.2478/s13380-012-0003-3
81
TaylorP.HobbsJ. N.BurroniJ.SiegelmannH. T. (2015). The global landscape of cognition: hierarchical aggregation as an organizational principle of human cortical networks and functions. Sci. Rep.5, 18112. 10.1038/srep18112
82
TekinR.TaglukM. E. (2017). Effects of small-world rewiring probability and noisy synaptic conductivity on slow waves: cortical network. Neural Comput.29, 679–715. 10.1162/NECO_a_00932
83
TelesfordQ. K.JoyceK. E.HayasakaS.BurdetteJ. H.LaurientiP. J. (2011). The ubiquity of small-world networks. Brain Connect.1, 367–375. 10.1089/brain.2011.0038
84
TiesingaP.SejnowskiT. J. (2009). Cortical enlightenment: are attentional gamma oscillations driven by ING or PING?Neuron63, 727–732. 10.1016/j.neuron.2009.09.009
85
TouboulJ.DestexheA. (2010). Can power-law scaling and neuronal avalanches arise from stochastic dynamics?PLoS ONE5, e8982. 10.1371/journal.pone.0008982
86
TsodyksM.UzielA.MarkramH. (2000). Synchrony generation in recurrent networks with frequency-dependent synapses. J. Neurosci.20, RC50. 10.1523/JNEUROSCI.20-01-j0003.2000
87
UhlhaasP. J.SingerW. (2006). Neural synchrony in brain disorders: relevance for cognitive dysfunctions and pathophysiology. Neuron52, 155–168. 10.1016/j.neuron.2006.09.020
88
van den HeuvelM. P.SpornsO. (2011). Rich-club organization of the human connectome. J. Neurosci.31, 15775. 10.1523/JNEUROSCI.3539-11.2011
89
van VreeswijkC.SompolinskyH. (1996). Chaos in neuronal networks with balanced excitatory and inhibitory activity. Science274, 1724. 10.1126/science.274.5293.1724
90
von BartheldC. S.BahneyJ.Herculano-HouzelS. (2016). The search for true numbers of neurons and glial cells in the human brain: a review of 150 years of cell counting. J. Comp. Neurol.524, 3865–3895. 10.1002/cne.24040
91
VuksanovićV.HövelP. (2014). Functional connectivity of distant cortical regions: role of remote synchronization and symmetry in interactions. Neuroimage97, 1–8. 10.1016/j.neuroimage.2014.04.039
92
WangQ.PercM.DuanZ.ChenG. (2010). Impact of delays and rewiring on the dynamics of small-world neuronal networks with two types of coupling. Phys. A Stat. Mech. Appl.389, 3299–3306. 10.1016/j.physa.2010.03.031
93
WattsD. J.StrogatzS. H. (1998). Collective dynamics of “small-world” networks. Nature393, 440–442. 10.1038/30918
94
WetschoreckF.KrabelT.KrishnamurthyS. (2020). 8080labs/ppscore: zenodo release. 10.5281/zenodo.4091345
95
WolframS. (1984a). Cellular automata as models of complexity. Nature311, 419–424. 10.1038/311419a0
96
WolframS. (1984b). Universality and complexity in cellular automata. Physica D10, 1–35. 10.1016/0167-2789(84)90245-8
97
YamamotoH.KubotaS.ShimizuF. A.Hirano-IwataA.NiwanoM. (2018). Effective subnetwork topology for synchronizing interconnected networks of coupled phase oscillators. Front. Comput. Neurosci.12, 17. 10.3389/fncom.2018.00017
98
YuanB.FangY.HanZ.SongL.HeY.BiY. (2017). Brain hubs in lesion models: Predicting functional network topology with lesion patterns in patients. Sci. Rep.7, 17908. 10.1038/s41598-017-17886-x
99
ZhouY.MilhamM. P.LuiY. W.MilesL.ReaumeJ.SodicksonD. K.et al. (2012). Default-mode network disruption in mild traumatic brain injury. Radiology265, 882–892. 10.1148/radiol.12120748
100
ZimmernV. (2020). Why brain criticality is clinically relevant: a scoping review. Front. Neural Circuits14, 54. 10.3389/fncir.2020.00054
Summary
Keywords
small-world, neuromodulation, neural oscillations, topology, simulation, network, criticality
Citation
Arvin S, Glud AN and Yonehara K (2022) Short- and Long-Range Connections Differentially Modulate the Dynamics and State of Small-World Networks. Front. Comput. Neurosci. 15:783474. doi: 10.3389/fncom.2021.783474
Received
26 September 2021
Accepted
03 December 2021
Published
25 January 2022
Volume
15 - 2021
Edited by
Ioanna Sandvig, Norwegian University of Science and Technology, Norway
Reviewed by
Matjaž Perc, University of Maribor, Slovenia; Claudia Battistin, Norwegian University of Science and Technology, Norway; Alexander E. Hramov, Innopolis University, Russia
Updates
Copyright
© 2022 Arvin, Glud and Yonehara.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Simon Arvin sarv@dandrite.au.dkKeisuke Yonehara keisuke.yonehara@dandrite.au.dk
†These authors share first authorship
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.