Fig. 15.1
Definition of the neuronal avalanche. (a) Brain slice with the recording site. (b) Comparison of seizure-like activity of ictal events, followed by tonic and clonic bursts with spontaneous activity. (c) Example of the collective time step, which is framed in 4 ms time bins from eight channels. The definition of an avalanche is separated by blank activity at the beginning and end of the events. The activated electrodes are counted as the avalanche size, and each event’s lifetime is the summation of the total time frames. (d) Distribution of different avalanche sizes plotted on a log-log scale. The neuronal avalanche could follow the power-law distribution, and its slope could be calculated as the α value. The original data were shuffled in order to disturb their spatiotemporal arrangement. Scale bar = 1 mm in A (Adapted with permission from Wu et al. 2014 BMC Neuroscience)
Seizure-like activity consists of ictal and tonic bursts (Fig. 15.1b, red line, upper panel), followed by clonic bursts (Fig. 15.1b, gray line, upper panel). The time-points at which each nLFPs exceeded the specific threshold is marked as a unit raster in the lower panel of Fig. 15.1b. The time units were summed together in a timescale-binned plot to calculate the avalanche size and lifetime. Figure 15.1c shows the cumulated time units framed by a time window (Δt, gray regions) from eight-channel recording. According to Wu et al. 2014, the avalanche is defined as “a series of activity separated by a blank at the beginning and end of the events”. The avalanche size was calculated as the total number of electrodes with active units, while the lifetime was calculated as the summation of the total time frame in each avalanche event. The distribution of the avalanche size with its probability P(s) were plotted on a log-log scale. A neuronal avalanche that has a fitted straight-line slope of α value indicates a power-law relationship and the event’s dependence. The 4AP-Bic group (red solid line) showed a power-law distribution with the α value around −1.5 (Fig. 15.1d). An important test was performed to assess whether event dependence is essential for the power-law relationship. Thus, the event dependence was disturbed, by randomly shuffling the data with regard to the order of temporal sequence and spatial arrangement of the events. Both the shuffling data (dashed red line) and spontaneous activity data (black line) showed an exponential distribution, which is a type of Poisson distribution in which the events occur independently (Fig. 15.1d).
15.2.2 Power Law Distributions
Power law distributions of event sizes are often seen in complex phenomena including semiconductor devices (Levinshtein et al. 2005), forest fires, earthquakes, phase transitions, financial market fluctuations, snow avalanches and many other instances (Bak 1996). For example, earthquake models integrate local rules in which forces at one site are distributed to nearest neighbors without dissipation. This conservation of forces is similar to the conservation of synaptic strengths (Royer and Pare 2003) and it could be a mechanism responsible for maintaining a network near the critical point. Computer simulations indicate that networks can be kept nearly critical levels when the total sum of synaptic strengths soars near a constant value (Hsu and Beggs 2006). This could be accomplished through a mechanism like synaptic scaling (Turrigiano and Nelson 2000), which has been observed experimentally. Finally, recently “burned” areas in forest fire models are refractory, while unburned areas are more likely to ignite. This balance of refractoriness and excitability combines to maintain the system near the critical point. Recent models of neuronal avalanches (Levina et al. 2005) have suggested that short-term synaptic depression and facilitation may also serve to drive neuronal networks toward the critical point where avalanches occur. Thus, an understanding of power laws in diverse complex systems can suggest mechanisms that might underlie criticality in neuronal networks. Computational models have also explored the potential relationship between neuronal avalanches and epilepsy (Beggs 2008; Hsu and Beggs 2006; Hsu et al. 2007, 2008).
Although the normal brain firing activity can be moderately-synchronous, in epileptic seizures, a group of neurons begin firing in an abnormal, excessive, and synchronized manner (McPhee and Hammer 2010; National Institute for Health and Clinical Excellence 2012; Plenz 2012; Yang et al. 2012). This excessive firing may occur due to structural or functional anomalies within the epileptic brain. Such abnormal activity could be produced either by hyper excited neurons acting independently or it could involve abnormal interactions among many neurons. Many forms of collective activity including waves, spirals, oscillations, synchrony, and neuronal avalanches have been identified in abnormal epileptic firing. All these emergent activity patterns have been hypothesized to show pathologic signatures associated with epilepsy (Hobbs et al. 2010). From this perspective, epileptic activity would occur when regulatory mechanisms failed and the network entered a super critical regimen. Operating at the critical point depends on the appropriate balance between excitation and inhibition implies a structured activity that is far from random. In addition, the activity in one neuron would, on average, lead to activity in other neurons, amplifying activity excessively and possibly leading to seizures (Hobbs et al. 2010). Therefore, when network activity is randomly shuffled, it no longer follows a power law distribution characteristic of avalanches (Beggs 2008).
Hobbs et al. (2010) examined neural activity from human and rat cortical tissue in local cortical networks using 60 channel multielectrode arrays to record local field potentials in brain slices removed from the most active epileptogenic area (as identified by intraoperative electrocorticography). In the human cortical tissue they found periods of pronounced hyperexcitability and lack of a clear power law in avalanche size distributions. Analysis showed that during these periods, there was a significant positive correlation between the branching parameter and the firing rate, suggesting a positive feedback loop. This aspect was not present in the activity examined in rat tissue. These results indicate that cortical tissue removed from pediatric epilepsy patients produces aberrant neuronal avalanches (Hobbs et al. 2010).
15.2.3 Self-Organized Criticality
Self-organized criticality (SOC) represents a property of complex dynamic systems that evolve to a critical state, capable of producing scale-free energy fluctuations. A characteristic feature of dynamical systems exhibiting SOC is the power-law probability distributions that describe the dynamics of energy release. Worrell and colleagues in their study investigated the probability distribution of in vivo pathological energy fluctuations in human epileptic hippocampus by analyzing data from seven consecutive patients with temporal lobe epilepsy who required depth electrode iEEG monitoring during evaluation for epilepsy surgery (Worrell et al. 2002). Contacts that recorded the earliest clear seizure onset on iEEG delineated the seizure onset zone that was determined by visual inspection as located within the mesial (middle) temporal lobe in each patient. Typical waveforms of interest were epileptiform spikes, sharp waves, and sharp and slow wave complexes. They concluded that the probability densities of interictal epileptic energy fluctuations and the quiescent time between successive fluctuations exhibit power-law scaling, which provides evidence for SOC in human epileptic hippocampus. They hypothesized that interictal epileptiform discharges are a mechanism for energy release within epileptic brain, and that these events may provide a method for identifying the network involved in seizure generation and they even may assure a physiological mechanism for preventing seizures (Worrell et al. 2002).
15.3 Epileptic Seizures
15.3.1 Definitions, Cause, Symptoms and Features of Epileptic Seizures
Epilepsy has been defined as a chronic, complex neurological disorder associated with abnormal electrical activity in the brain, marked by sudden recurrent episodes of sensory/motor disturbance, changes in behavior, with or without loss of consciousness, and/or convulsions (Nunes et al. 2012). This excessive activity can be produced by hyperexcited neurons acting independently or involve abnormal interactions among many neurons (Hobbs et al. 2010). Epilepsy is one of the most common chronic neurological affections, with more than 50 million patients worldwide and approximate 2 million new cases each year. Epilepsy responds to treatment in about 70 % of cases. Not all cases of epilepsy are lifelong, a substantial number of people improve to the point that medication is no longer needed. One seizure does not signal epilepsy (up to 10 % of people worldwide have one seizure during their lifetimes) (Chang and Lowenstein 2003). Epilepsy has been defined since 2005 by two or more unprovoked seizures, >24 h apart (Fisher et al. 2005). The definition of epilepsy has been recently revised in order to consider the diagnosis even after occurrence of a single unprovoked seizure, providing that we can demonstrate an enduring predisposition for recurrence, similar to the general recurrence risk after two seizures (more than 60 %) over the next 10 years. For instance a patient can have a single unprovoked seizure after trauma, a stroke, or central nervous system infection. A patient with such brain disorders has a high risk of developing epilepsy after a single unprovoked seizure.
The tendency to respond to particular stimuli with seizures also meets the conceptual definition of epilepsy. Reflex epilepsies are also associated with an enduring abnormal predisposition to have recurrent seizures.
The document also mentions as a general agreement that epilepsy should no longer be considered a disorder but a disease, the term implying a more substantial, long lasting derangement of neuronal functionality (Fisher et al. 2014).
15.3.1.1 Classification of Seizures and Epilepsies
Generalized epileptic seizures are defined as originating within bilaterally spread networks that, at some point, are rapidly involved. These networks do not necessary include the entire cortex and could involve both cortical and subcortical structures. Generalized seizures can be asymmetric in clinical appearance. Focal epileptic seizures are defined as originating within networks restricted to one hemisphere. They may be localized in a small area or have a wider spread. Depending on the seizure type, ictal onset (“area of cortex that initiates clinical seizures” (Rosenow and Luders 2001) is consistent from one seizure to another), with distinctive propagation patterns that can sometimes imply the contralateral hemisphere. Focal seizures might provoke alteration of consciousness or awareness and could evolve into secondary tonic-clonic generalization. Auras are purely subjective clinical manifestation usually occurring at the very onset of a focal seizure as a warning event. The epileptic nature is sometimes difficult to prove if the aura symptoms do not evolve into a more objective clinical pattern (Luders et al. 1998). Focal seizures might develop abruptly without any warning symptoms.
15.3.1.2 Cause Types (Etiology)
The International League Against Epilepsy (ILAE) Commission on Classification and Terminology has simplified the classification of seizures and they proposed the following three concepts:
1.
Genetic epilepsy represents the result of a known or supposed genetic error in which seizures are the main symptom of the disorder. In recent years, some types of epilepsy have been correlated to mutations in genes, mostly involving ion channels, assumptions derived from specific molecular genetic studies or from suitable adapted family studies. For example a mutation in the GABA1 gene has been detected in some members of a family with juvenile myoclonic epilepsy (Cossette et al. 2002). No definitive conclusions can be made because there is no knowledge regarding specific surrounding influences as causes of or factors that contribute to these forms of epilepsy.
2.
Structural or metabolic: This type of epilepsy refers to a condition or disease with a metabolic or structural background associated with a high risk of developing epilepsy. Structural lesions are frequently associated with focal seizures. They include cerebral changes resulting from, trauma, stroke, perinatal brain damage, malformations, brain tumors, and infections. These epilepsies are characterized by the presence in epileptogenic foci of residual neurons with no afferents and less dendritic spines, destroyed probably by infections, trauma, stroke, or other lesion. Structural lesions of the cerebral cortex can disrupt, inhibitory GABAergic interneurons, thereby minimizing the inhibition that controls large pyramidal cells (Menzler et al. 2011).
3.
“Unknown cause”: There are several cases wherein the underlying pathophysiology can’t be identified. These are usually encountered in children and young adults but can occur at any age in persons who have a family history of epilepsy or seizures. This epilepsy may have a genetic defect at its origin or it may be the result of a distinct and yet unknown disorder (Berg et al. 2010).
15.3.2 Power Law Distributions in Epilepsy
Diseases of central nervous system are often associated with altered brain dynamics (Expert et al. 2010). It has been hypothesized that the dynamical properties characterizing a critical state may be considered as an important marker of brain well-being in both health and disease (Plenz 2012). During epileptic seizures the distribution of phase-locking intervals (PLI) is providing additional evidence for the criticality hypothesis. Furthermore, deriving the distribution of PLI from electrocorticogram (ECoG) data as an indicator of critical brain dynamics has shown that the system deviates from scale-free behavior during seizures. All scales closely follow a power-law probability distribution during pre-ictal time intervals with the exponent between 22 and 23.5 (see Fig. 15.2). The apparent robustness of the power-law against exact conditions (different anatomical regions with varying number of channels) strengthens the hypothesis of the relevance of a critical state in human brain dynamics. While the PLI distribution followed a power-law in time intervals preceding the seizure onset, a deviation from power-law behavior was observed in intervals containing the seizure attack.
Fig. 15.2
The distribution of phase-locking intervals deviates from a power-law during epileptic seizures. Top: The electrocorticogram (ECoG) recording shows the onset of a focal epileptic seizure attack around 300 s time. Bottom: Cumulative distributions of phase-locking intervals (PLI) are obtained during three time intervals of 150 s: preictal (left), ictal (middle) and postictal (right). Dashed lines indicate a power-law with exponent 23.1. While the distribution appears to follow a power-law during the pre-ictal period, intervals of increased phase-locking disturb this characteristic distribution with the onset of seizure activity. Data shown are from a patient scale 3, corresponding to the frequency band 25–12.5 Hz (Adapted with permission from Meisel et al. 2012, PLOS Computational Biology)
Figure 15.2 shows the distribution of PLI derived from a pre-ictal, an ictal and a postictal time interval. The probability to find longer PLI increased during attacks thereby destroying the scale-free property of the original distribution. After the seizure this distribution slowly relaxed back to a power-law. In Fig. 15.2 this relaxation is not yet complete in the postictal time interval as there is still some residual seizure dynamics in the ECoG recording.
15.3.3 Self-Organized Criticality in the Brain
Critical systems have been defined as systems that are close to a critical point, near the boundary of an order-disorder phase transition. At criticality, these systems can avoid being trapped in one of two extreme cases: a disordered state (when interactions are too weak and the system is dominated by noise) or a globally ordered state in which all elements are locked (when interactions are too strong and the system is completely static). A dualism is essential for a complex system, like the brain, to function: it must maintain some order to ensure coherent functioning (i.e., generate a reproducible behavior in response to a certain stimulus) while allowing for a certain degree of disorder to enable flexibility (i.e., adapt to varying external conditions). Such dualism is possible at criticality.
While many degrees of order/disorder are possible, the subtle balance between order and disorder at criticality manifests itself in certain general statistical properties: critical systems exhibit spatial and temporal correlations that are long range (i.e., on scales that are larger than those on which mutual interactions take effect) and follow power-law distributions. Recent research has shown that brain networks, by being in the critical state, optimize their response to inputs and maximize their information processing ability (Shew and Plentz 2013).
Haimovici and colleagues have presented a simple brain model that, if tuned to criticality, explains the broad range of experimental observations of human brain activity, in particular, reproducing key findings obtained with functional magnetic resonance imaging (fMRI) (Haimovici et al. 2013). fMRI research has been able to deliver an important observation: the human brain at rest exhibits a large-scale spatiotemporal organization into distinct functional networks—so-called resting state networks (RSN) (Fox and Raichle 2007). RSNs are areas of the resting brain—measured in subjects performing any cognitive, or motor tasks—in which fluctuations of neural activity are correlated, as revealed by the fact that BOLD signal fluctuations within the same network are synchronous. Each RSN can be related to a specific set of cortical areas associated with certain functions: cognitive, sensory (visual, auditory), and motor RSNs, for instance, have been identified (see Fig. 15.3).
Fig. 15.3
Functional magnetic resonance imaging (fMRI) experiments have revealed that the brain at rest is organized into several areas in which fluctuations of brain activities are correlated, so-called resting state networks (RSN). From top to bottom: medial visual (VisM), lateral visual (VisL), auditory (Aud), and sensory-motor (SM) RSNs. (Right columns) Results from the work of Haimovici et al. 2013 show that a simple model can reproduce the statistical properties of RSNs only if the model is tuned to criticality (at TC) (Adapted with permission from Haimovici et al. 2013, Phys Rev Letters)
What Haimovici and colleagues found is that such a model can lead to activity clusters similar to those found for RSN activity (Haimovici et al. 2013). But for the model to match the experimental data, the activation threshold had to be set exactly at the level at which their model becomes critical, as illustrated in Fig. 15.3. At criticality, their model predicts a number of statistical properties that are consistent with experiments: the brain forms activity clusters whose size follows a power law with slope of −3/2, the hallmark of neuronal avalanches (Beggs and Plenz 2003), corresponding to a peak in the size of the second-largest cluster, as found in percolation models (Margolina et al. 1982); the correlation length (the distance at which two points in the system behave independently) and its fluctuations diverge and match those seen in human brain data.
15.4 Modular Signatures in Epileptic Seizures
15.4.1 Functional vs. Structural Modularity of the Epileptic Brain
The human neocortex consists of a large number of minicolumns in parallel vertical arrays (Mountcastle 1957, 1997; Buxhoeveden and Casanova 2002; Casanova et al. 2007; Shepherd and Grillner 2010). Minicolumns are the first step in a nested series of nodes or echelons of increasing complexity (Mountcastle 1997; Szentagothai 1975). Within minicolumns, cortical neurons are aggregated into five horizontal layers (or laminae), namely two supra-granular, one granular and two infra-granular layers. Other levels of modular organization include multiple minicolumns, macrocolumns, and large-scale networks of macrocolumns that are interconnected with the entire brain (Buxhoeveden and Casanova 2002; Opris and Casanova 2014).
In contrast to the sparse but organized connectivity between the modules of control subjects, brain connectivity of epileptic patients shows a configuration where nodes in a functional module are more connected to different functional modules. Recently, Vaessen et al. (2014) examined connectivity of whole brains of children with frontal lobe epilepsy (FLE) and compared their structural and functional connectivity with the same in healthy controls. Their measurements of the functional connectivity was derived from the dynamic fluctuations of the fMRI, while the structural connectivity was determined from fiber tractograms of diffusion weighted MRI. The whole brain network patterns of connectivity were characterized with graph theoretical metrics and further decomposed into modules. Then, the graph metrics with the extracted connectivity within and between modules were related to cognitive performance.
As shown in Fig. 15.4a, the modularity algorithm extracted four modules from the averaged functional connectivity matrix over all subjects. Spatial organization of module #1 was considered the “default mode network” with network nodes distributed in the frontal, temporal and parietal lobes (Vaessen et al. 2014). Module #2 was distributed over the frontal and subcortical regions. Module #3 was localized in the occipital lobe. Module #4 nodes were distributed over frontal, temporal and occipital regions. All four modules seemed to be symmetric to the inter-hemispheric fissure. As shown in Figs. 15.4a, b when the structural connections were organized according to the modularity index of the frontal cortex, the structural modularity revealed bilateral structural sub networks. For the structural connectivity matrix the modularity algorithm determined only two modules, which were highly symmetric over the two hemispheres. After functional connections were organized according to these structural connectivity modules, the sub organization vanished (Fig. 15.4c, d). An interactive view of functional connectivity and structural connectivity was displayed as between-module, within-module averaged over all modules and individual within-module connections (Fig. 15.4e, f). It was further shown that functional “disturbances” of epileptic children were related to increased clustering and stronger modularity compared to healthy controls, which was accompanied by stronger within- and weaker between-module functional connectivity. While structural modularity increased with stronger cognitive impairment, it was concluded that decreased coupling between large-scale functional network modules may represent a hallmark for impaired cognition in childhood FLE.
Fig. 15.4
Functional vs. structural brain modularity. (a) Functional connectivity showing the average connection matrices sorted by module. Colored rectangles indicate the modules. High within-module connectivity is illustrated by the higher values (more hot colors), while between-module connectivity is shown by more sparse (more cold colors). (b) Structural connectivity was sorted by functional modules. The functional modules are organized bilaterally, while the structural connectivity has strong inter-hemispheric connectivity and low intra-hemispheric connectivity clearly visible in the block patterns. (c) Functional connectivity matrix was sorted by the modular organization derived from the structural connectivity. The two structural connectivity modules represent the left and right hemisphere. Functional connectivity shows that strong inter-hemispheric connections are present within the two modules. (d) The structural connectivity was sorted by structural connectivity modularity. Strong intra-hemispheric connections are obvious, while inter-hemispheric connections (and thus between-module connections) are weaker. (e) Illustration of modular organization of functional connectivity. Within-module connections are colored as in panel A. (f) An alternative presentation of the modular organization. The nodes of each separate module are depicted spatially segregated. The gray lines indicate the between-module connections (Adapted with permission from Vaessen et al. 2014, PLOS one)
15.4.2 Spatial Scale of Epileptogenicity Biomarkers Matches Minicolumns Size
High-frequency oscillations (HFO) in the 80(100)–500 Hz range are considered an important biomarker of cortical epileptogenicity. Although HFO may be present in non-epilepogenic brain structures, it has been shown that the ripple (<250 Hz) and particularly fast ripples (>250 Hz) subbands have are highly specific to the seizure onset zone (Zijlmans et al. 2012). Using dense 2D microelectrode arrays (MEA) having a spacing of 0.4 mm between electrodes, Schevon and colleagues have shown that about ~90 % of the recorded HFOs were limited to a single channel (Schevon et al. 2009). This is consistent with the idea that most of the spontaneous HFO events are confined within a minicolumn, and spreading the activity requires a recruitment process of several minicolumns, possibly through an avalanche process following a power-law, as suggested by several other studies using intracranial recordings (Worrell et al. 2002; Wu et al. 2014).