# Neuronal Avalanches to Study the Coordination of Large-Scale Brain Activity: Application to Rett Syndrome

Rosaria Rucco, Pia Bernardo, Anna Lardone, Fabio Baselice, Matteo Pesoli, Arianna Polverino, Carmela Bravaccio, Carmine Granata, Laura Mandolesi, Giuseppe Sorrentino, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Rosaria Rucco, Pia Bernardo, Anna Lardone, Fabio Baselice, Matteo Pesoli, et al. "Neuronal Avalanches to Study the Coordination of Large-Scale Brain Activity: Application to Rett Syndrome." *Frontiers in Psychology*, 2020, 11, 10.3389/fpsyg.2020.550749. [hal-03599050](https://hal.science/hal-03599050)

HAL Id: hal-03599050

https://hal.science/hal-03599050

Submitted on 4 Jul 2023

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Check for
updates

# Neuronal Avalanches to Study the Coordination of Large-Scale Brain Activity: Application to Rett Syndrome

Rosaria Rucco<fnref n="1" />,²†, Pia Bernardo³,⁴†, Anna Lardone<fnref n="1" />, Fabio Baselice⁵, Matteo Pesoli<fnref n="1" />, Arianna Polverino⁶, Carmela Bravaccio³, Carmine Granata², Laura Mandolesi⁷, Giuseppe Sorrentino<fnref n="1" />,²,⁶* and Pierpaolo Sorrentino²,⁵,⁸

OPEN ACCESS

**Edited by:**
Elisa Pedroli,
Istituto Auxologico Italiano (IRCCS),
Italy

**Reviewed by:**
Jean-Philippe Thivierge,
University of Ottawa, Canada
Ekhard K.H. Salje,
University of Cambridge,
United Kingdom

*Correspondence:
Giuseppe Sorrentino
giuseppe.sorrentino@uniparthenope.it

**Specialty section:**
This article was submitted to
Quantitative Psychology
and Measurement,
a section of the journal
Frontiers in Psychology

Received: 10 April 2020

Accepted: 28 September 2020

Published: 27 October 2020

**Citation:**
Rucco R, Bernardo P, Lardone A, Baselice F, Pesoli M, Polverino A, Bravaccio C, Granata C, Mandolesi L, Sorrentino G and Sorrentino P (2020) Neuronal Avalanches to Study the Coordination of Large-Scale Brain Activity: Application to Rett Syndrome. *Front. Psychol.* 11:550749.
doi: 10.3389/fpsyg.2020.550749

Many complex systems, such as the brain, display large-scale coordinated interactions that create ordered patterns. Classically, such patterns have been studied using the framework of criticality, i.e., at a transition point between two qualitatively distinct patterns. This kind of system is generally characterized by a scale-invariant organization, in space and time, optimally described by a power-law distribution whose slope is quantified by an exponent α. The dynamics of these systems is characterized by alternating periods of activations, called avalanches, with quiescent periods. To maximize its efficiency, the system must find a trade-off between its stability and ease of propagation of activation, which is achieved by a branching process. It is quantified by a branching parameter σ defined as the average ratio between the number of activations in consecutive time bins. The brain is itself a complex system and its activity can be described as a series of neuronal avalanches. It is known that critical aspects of brain dynamics are modeled with a branching parameter σ = , and the neuronal avalanches distribution fits well with a power law distribution exponent α = -3/2. The aim of our work was to study a self-organized criticality system in which there was a change in neuronal circuits due to genetic causes. To this end, we have compared the characteristics of neuronal avalanches in a group of 10 patients affected by Rett syndrome, during an open-eye resting-state condition estimated using magnetoencephalography, with respect to 10 healthy subjects. The analysis was performed both in broadband and in the five canonical frequency bands. We found, for both groups, a branching parameter close to 1. In this critical condition, Rett patients show a lower distribution parameter α in the delta and broadband. These results suggest that the large-scale coordination of activity occurs to a lesser extent in RTT patients.

**Keywords:** self-organized criticality, critical state, neuronal avalanches, branching process, magnetoencephalography, Rett syndrome

<footnote n="1">¹ Department of Motor Sciences and Wellness, University of Naples “Parthenope,” Naples, Italy, ² Institute of Applied Sciences and Intelligent Systems, National Research Council (CNR), Pozzuoli, Italy, ³ Department of Medical and Translational Science, Child Neuropsychiatry Unit, University of Naples “Federico II,” Naples, Italy, ⁴ Department of Neuroscience, Pediatric Psychiatry and Neurology, Santobono-Pausilipon Children’s Hospital, Naples, Italy, ⁵ Department of Engineering, University of Naples “Parthenope,” Naples, Italy, ⁶ Hermitage Capodimonte Hospital, Naples, Italy, ⁷ Department of Humanistic Studies, University of Naples “Federico II,” Naples, Italy, ⁸ Institut de Neurosciences des Systèmes, Aix-Marseille Université, Marseille, France</footnote>

<footnote>† These authors have contributed
equally to this work</footnote>
## INTRODUCTION

Rett syndrome (RTT) is a severe neurodevelopmental disorder, characterized by clinically normal development for the first 12–18 months of life, when an overall arrest and regression of the psychomotor development begins. The patients lose the acquired verbal skills, the normal motor hand function is replaced with stereotyped movements, and often ataxia is present, which seriously compromises motor skills. Intellectual disability, epilepsy, autonomic dysfunction, breathing abnormalities, anxiety, and orthopedic problems are the other most important clinical findings in RTT (Neul et al., 2010). In 1999, Amir et al. (1999) identified in a mutation in the X-linked MECP2 gene, encoding the methyl CpG-binding protein 2 (MeCP2), the most common cause of RTT.

The dysfunction of *MeCP2* provokes multiple effects, such as impaired neuronal maturation, altered GABAergic signaling, and, more importantly, a local imbalance between neuronal excitation and inhibition at the circuit level (Pohodich and Zoghbi, 2015). Importantly, all these molecular and circuital evidence are scarcely linked to the clinically evident deficits in higher cognitive functions. In fact, effective coordination of large-scale activity is considered necessary for the emergence of high-level cognitive abilities. All these features make RTT a potential model for studying how an alteration of the normal development process interferes with the activity of brain networks.

Consequently, we borrowed from physics some tools that allow the description of large-scale activity. It has been shown that the healthy brain has a higher number of configurations compared to the unhealthy brain (Sorrentino et al., 2018). Importantly, it has been argued that flexibility is possible given that the brain sets itself to operate near a critical state (Beggs and Plenz, 2003). Bak et al. (1987) first described the idea of “self-organized criticality” (SOC) defined as a system autonomously evolving to a critical state regardless of the initial condition. As he defines, “the critical state is an attractor for the dynamics of a system.” In the critical state, the activity can spread across the system in a controlled way, avoiding that activations either spread uncontrollably or quickly die out (Zapperi et al., 1995;

Harris, 2002). Such state is classically quantified by a branching parameter $\sigma$ defined as the average ratio between the number of activations in consecutive time bins (i.e., the sample units) (De Carvalho and Prado, 2000). More specifically, when $\sigma < 1$, the system is in a subcritical condition in which the activations quickly fade out. On the other hand, when $\sigma > 1$, the opposite happens, and the system is in a supercritical state whereby the activations spread across the entire system, leading to unstable runaway activations. Finally, when $\sigma = 1$, the system is in a critical state with balanced, fine-tuned propagation of the activity (Bak et al., 1987; Harris, 2002).

SOC systems are generally characterized by a scale-invariant organization in space and time. In other words, the fluctuations that occur are similar in all scales of time and space, and they can be optimally described by fractal metrics, typically a power law distribution whose slope quantifies the amplitude distribution of such fluctuations (Paczuski et al., 1996; Newman, 2007; Clauset et al., 2009).

Importantly, power law distributions do not prove criticality (Clauset et al., 2009). However, they do provide a statistical measure of the occurrence of large-scale coordinated activity, perhaps providing a useful tool to probe hypothesis about the nature of high-level cognitive deficits, such as those observed in Rett syndrome.

Generally speaking, SOC can be found in very different types of phenomena, like in many natural systems such as earthquakes (Gutenberg and Richter, 1956) or forest fires (Malamud et al., 1998), in which when a unit exceeds a threshold, other units follow it, thus giving rise to a cascade (or avalanche) that spreads across the entire system. The presence of power law-distributed bursts of activations, called avalanches, alternating with quiescent periods, characterizes the dynamics of SOC systems (Lombardi et al., 2016) likely maximizing the efficiency of communication among its elements. However, it is to be considered that the analysis of avalanches, with a power law distribution, varies in heterogeneous/disordered systems and that to get information on the non-equilibrium dynamics, it would be useful to also analyze the crackling noise (Salje and Dahmen, 2014).

In this article, we hypothesize that the genetic defect causing RTT impairs the whole-brain dynamic and relates to clinical symptoms. To test our hypothesis, we compared the distributions of the size of the neuronal avalanches, in order to capture the critical features of brain activity, in 10 patients with RTT and

**Abbreviations:** HS, healthy subjects; MECP2, methyl CpG-binding protein 2; MEG, magnetoencephalography; ROIs, regions of interests; RTT, Rett syndrome; SOC, self-organized criticality.

<table><caption>TABLE 1 | The mean value, for Rett syndrome (RTT) patients and health subject (HS) group, of branching parameter $\sigma$, for all frequency bands in all time bins.</caption><thead><tr><th rowspan="2"></th><th colspan="2">Δ t = 1</th><th colspan="2">Δ t = 2</th><th colspan="2">Δ t = 3</th><th colspan="2">Δ t = 4</th><th colspan="2">Δ t = 5</th></tr><tr><th>RTT patients</th><th>HS</th><th>RTT patients</th><th>HS</th><th>RTT patients</th><th>HS</th><th>RTT patients</th><th>HS</th><th>RTT patients</th><th>HS</th></tr></thead><tbody><tr><td>Delta band</td><td>0.900</td><td>0.901</td><td>0.932</td><td>0.911</td><td>0.943</td><td>0.919</td><td>0.955</td><td>0.971</td><td><strong>1.031</strong></td><td><strong>1.016</strong></td></tr><tr><td>Theta band</td><td><strong>0.911</strong></td><td><strong>0.964</strong></td><td>1.137</td><td>1.179</td><td>1.139</td><td>1.226</td><td>1.189</td><td>1.286</td><td>1.270</td><td>1.329</td></tr><tr><td>Alpha band</td><td><strong>0.974</strong></td><td><strong>0.964</strong></td><td>1.084</td><td>1.095</td><td>1.128</td><td>1.164</td><td>1.135</td><td>1.171</td><td>1.219</td><td>1.231</td></tr><tr><td>Beta band</td><td>0.835</td><td>0.848</td><td><strong>0.981</strong></td><td><strong>0.985</strong></td><td>1.071</td><td>1.104</td><td>1.106</td><td>1.172</td><td>1.098</td><td>1.179</td></tr><tr><td>Gamma band</td><td>0.848</td><td>0.865</td><td><strong>1.003</strong></td><td><strong>1.041</strong></td><td>1.091</td><td>1.169</td><td>1.159</td><td>1.244</td><td>1.167</td><td>1.244</td></tr><tr><td>Broadband</td><td>0.925</td><td>0.942</td><td><strong>1.021</strong></td><td><strong>1.043</strong></td><td>1.095</td><td>1.097</td><td>1.251</td><td>1.192</td><td>1.227</td><td>1.261</td></tr></tbody></table>

In bold are reported the time bin chosen for each frequency band for the further analysis.
<table><caption>**TABLE 2** | To appreciate the robustness of the results presented, we reported the variance values around the branching parameters, for RTT patients and HS group, when dataset is split in two segments.</caption><thead><tr><th rowspan="2"></th><th colspan="2">Segment 1/2</th><th colspan="2">Segment 2/2</th></tr><tr><th>RTT patients</th><th>HS</th><th>RTT patients</th><th>HS</th></tr></thead><tbody><tr><td>Delta band<br>(Δt = 5)</td><td>1.034 (0.007)</td><td>1.021 (0.003)</td><td>1.028 (0.002)</td><td>1.011 (0.004)</td></tr><tr><td>Alpha band<br>(Δt = 1)</td><td>0.973 (0.001)</td><td>0.966 (0.001)</td><td>0.975 (0.001)</td><td>0.962 (0.000)</td></tr><tr><td>Broadband<br>(Δt = 2)</td><td>1.018 (0.010)</td><td>1.045 (0.001)</td><td>1.024 (0.001)</td><td>1.041 (0.001)</td></tr></tbody></table>

For brevity, we reported only delta, alpha, and broadbands.

10n matched healthy subjects (HS). We expect that large-scale activation would be occurring less in patients showing marked deficits in higher cognitive functions. We based our analyses on source-reconstructed magnetoencephalographic (MEG) data acquired during resting state.

# MATERIALS AND METHODS

## Participants and Clinical Assessment

Participants were recruited from the Department of Translational Medical Sciences, Child Neuropsychiatry in “Federico II” University of Naples, Italy. We studied 10 female patients with clinical diagnosis of RTT (age 24.30 ± 8 years) based on the Neul revised criteria (Neul et al., 2010) and confirmed by mutation in the *MECP2* gene, and 10 healthy female individuals (HS) (age 26.10 ± 6.84 years). This study complied with the Declaration of Helsinki and was approved by the local ethics committee. Written informed consent has been granted by all participants (or their legal guardians).

## Acquisition

The data were acquired using a MEG system equipped by 163 magnetometers SQUID (Superconducting Quantum Interference Device) (Sorrentino et al., 2018), placed in a magnetically shielded room (AtB Biomag, Ulm, Germany), in order to

<figure><img src="image_4.png" /><figcaption>**FIGURE 1** | Data analysis pipeline. (A) Raw magnetoencephalography (MEG) signals recorded by 154 sensors (a subset displayed here). (B–D) Respectively noisy channel, cardiac artifact, blinking artifact, removed during preprocessing phase. (E) MEG signals after artifact removal and noise cleaning. (F) Coregistration between MEG signals and MRI template. (G) Source reconstruction (beamforming). (H) A sequence of four neuronal avalanches. (I) How to calculate the dimension of a neuronal avalanche, in particular, the first (in red square) of the (H) image.</figcaption></figure>
<table><caption>**TABLE 3** | To appreciate the robustness of the results presented, we reported the variance values around the branching parameters, for RTT patients and HS group, when dataset is split in three segments.</caption><thead><tr><th rowspan="2"></th><th colspan="2">Segment 1/3</th><th colspan="2">Segment 2/3</th><th colspan="2">Segment 3/3</th></tr><tr><th>RTT patients</th><th>HS</th><th>RTT patients</th><th>HS</th><th>RTT patients</th><th>HS</th></tr></thead><tbody><tr><td>Delta band (Δt = 5)</td><td>1.030 (0.010)</td><td>1.020 (0.005)</td><td>1.033 (0.004)</td><td>1.016 (0.009)</td><td>1.029 (0.005)</td><td>1.011 (0.002)</td></tr><tr><td>Alpha band (Δt = 1)</td><td>0.971 (0.002)</td><td>0.966 (0.001)</td><td>0.974 (0.001)</td><td>0.961 (0.000)</td><td>0.976 (0.001)</td><td>0.965 (0.001)</td></tr><tr><td>Broadband (Δt = 2)</td><td>1.018 (0.001)</td><td>1.040 (0.002)</td><td>1.020 (0.002)</td><td>1.046 (0.001)</td><td>1.024 (0.001)</td><td>1.043 (0.002)</td></tr></tbody></table>

For brevity, we reported only delta, alpha, and broadbands.

<figure><img src="image_5.png" /></figure>

reduce the external noise. Of all squids, 154 are positioned to be as close as possible to the head of the subject, while the remaining ones are more distant so as to measure environmental

noise (reference magnetometers). All the subjects of each group underwent a 7-min resting-state MEG acquisition with open eyes, divided in two segments. To evaluate the right
Avalanches Duration - SD=3

<figure><img src="image_6.png" /></figure>

**FIGURE 3** | Avalanche duration distributions for each RTT patient (first line) and for each HS (second line) for Delta (Δt = 5), Alpha (Δt = 1) and Broad (Δt = 2) bands (the threshold chosen is SD = 3). The blue lines represent the avalanche duration distributions for each subject (RTT patients in the top row, Healthy Subjects in the bottom row). The red lines represent the sum of the avalanche duration distributions across all subjects.

Size vs Duration - SD=3

<figure><img src="image_7.png" /></figure>

**FIGURE 4** | Size versus duration for each RTT patient (first line) and for each HS (second line) for delta (Δt = 5), alpha (Δt = 1), and broad (Δt = 2) bands (the threshold chosen is SD = 3).

position of the head under the helmet, we used Fastrak (Polhemus®) to acquire the position of four coils (attached to the head) and of four anatomical landmarks (nasion, right and left pre-auricular points, and vertex of the head). Before each acquisition segment, the head position in the helmet was obtained. During the acquisition, two electrodes for the electrocardiogram and two for the electrooculogram (Gross et al., 2013) were also acquired.

Preprocessing

As previously described (Jacini et al., 2018; Rucco et al., 2019), the MEG signals, after an anti-aliasing filter, were acquired with a sampling frequency of 1,024 Hz. A fourth-order Butterworth IIR band-pass filter in the 0.5- to 100-Hz band was subsequently

applied to the acquired signals. Environmental noise, measured by reference magnetometers, was removed by using the principal component analysis. MEG data were cleaned of physiological artifacts, such as eye blinking and heart activity, by means of independent component analysis (Sorriso et al., 2019). Visual inspection was used for identification of noisy channels. For all the preprocessing steps, we used the FieldTrip toolbox (Oostenveld et al., 2011).

Source Reconstruction

Time series of neuronal activity were reconstructed in 116 regions of interests (ROIs) based on the automated anatomical labeling (AAL) atlas (Tzourio-Mazoyer et al., 2002) using a linearly constrained minimum variance beam-former
Size Distributions

<figure><img src="image_8.png" /><figcaption>**FIGURE 5** | Average size distributions for both groups, for Delta (Δt = 5), Alpha (Δt = 1) and Broad (Δt = 2) bands, as a function of the threshold (from 1 to 3 SD). The blue lines represent the average of the avalanche size distributions across either the RTT patients or the HS (as indicated in the y-label). The red lines represent the fitted power law.</figcaption></figure>

algorithm (Van Veen et al., 1997; Nolte, 2003) based on
MRI template and then filtered both in broadband (0.5–
48 Hz) and in the five classical frequency bands [delta

(0.5–4.0 Hz), theta (4.0–8.0 Hz), alpha (8.0–13.0 Hz), beta
(13.0–30.0 Hz), and gamma (30.0–48.0 Hz)]. From the
116 ROIs of the AAL Atlas, we have excluded 26 ROIs
Avalanches Duration

<figure><img src="image_9.png" /><figcaption>FIGURE 6 | Avalanches duration distributions for both groups for Delta (Δt = 5), Alpha (Δt = 1) and Broad (Δt = 2) bands as a function of the threshold (from 1 to 3 SD). The blue lines represent the avalanche duration distributions across either the RTT patients or the HS (as indicated in the y-label). The red lines represent the sum of the avalanche duration distributions across all subjects.</figcaption></figure>

corresponding to the cerebellum because of their low reliability in MEG (Lardone et al., 2018). Hence, we considered a total of 90 ROIs. We have resampled the source-space time series at 512 Hz.

Signal Discretization

For each ROI, the time series were discretized with five different time bin durations Δt, each one multiple of Δtmin (19 ms), corresponding to the sampling period.

An event was identified by a positive or a negative excursion of an area, in a bin, beyond a threshold, defined as ±3 SD of the signal amplitude, which was a tradeoff between a lower threshold (leading to the detection of more spurious noise events in addition to real ones) and a higher threshold (missing real events). Indeed, MEG systems rely on squid sensors, which measure the magnetic field generated by brain activity. However, such sensors cannot identify the field direction. For this reason, both positive and negative excursions were picked up without distinction. As a result, the event was identified when the absolute value of this excursion exceeded the chosen threshold.

Branching Parameter

As previously described (Beggs and Plenz, 2003; Shriki et al., 2013), an avalanche is defined as a sequence of contiguous time bins starting when at least one ROI is active and ending when all regions are inactive. The number of events in all ROIs in an avalanche corresponds to its size.

For each subject and for each time bin size, the branching parameter σ was estimated by calculating, for each avalanche, the averaged (over all the time bins) ratio of the number of events between the subsequent time bin (descendants) and that in the current time bin (ancestors), and then averaging it over all the cascades (Bak et al., 1987). More specifically:

$$ \sigma = \frac{1}{N_{aval}} \sum_{i=1}^{N_{aval}} \sigma_i \qquad (1) $$

$$ \sigma_i = \frac{1}{N_{bin} - 1} \sum_{j=1}^{N_{bin}-1} \frac{n_{events}(j+1)}{n_{events}(j)} \qquad (2) $$
Size Duration

<figure><img src="image_10.png" /><figcaption>FIGURE 7 | Size versus duration for both groups for delta (Δt = 5), alpha (Δt = 1), and broad (Δt = 2) bands varying the threshold (1 to 3 SD).</figcaption></figure>

where σᵢ is the branching parameter of the *i*-th avalanche in the dataset, Nₜₜₜₜ is the total amount of bins in the *i*-th avalanche, nᵣᵣᵣᵣ is the total number of events active in the *j*-th bin, and Nₜₜₜₜ is the total number of avalanche in the dataset.

Neuronal Avalanches and Power Law
Parameter

The statistical distribution of the avalanche size x was assumed to
be described by a power law function:

$$
P_{\alpha} = \begin{cases} C_{\alpha}x^{\alpha} & x_{min} \le x \le x_{max} \\ 0 & \text{otherwise} \end{cases} \tag{3}
$$

where C_α is a normalization factor. The parameters x_min was
set to 1 (the minimal avalanche size), and x_max was set to 1.5
times the total number of ROIs included in the analysis. We
used truncated power law to take finite size effect into account.
The distribution parameter α was estimated from the data by
implementing a maximum likelihood algorithm (Kay, 1993).
When the branching parameter σ is equal to 1, the system
is in a critical state. In other words, there is a long distance
propagation of neuronal activity without runaway excitation. For
this reason, we calculated at which time bin, for each frequency
band, the branching parameter (as average distance) was the
closest to 1 for both groups. Such choice was made in order to
Duration vs Pause - SD=3

<figure><img src="image_11.png" /><figcaption>FIGURE 8 | Duration versus pause for both groups for delta (Δt = 5), alpha (Δt = 1), and broad (Δt = 2) bands varying the threshold (1 to 3 SD).</figcaption></figure>

ECDF of Hellinger distance

<figure><img src="image_12.png" /><figcaption>FIGURE 9 | The empirical cumulative distribution functions of the Hellinger distance for RTT patients and HS separately, in all time bins, for delta, alpha, and broadbands. The chosen threshold is SD = 3.</figcaption></figure>

have the two groups in the same condition, in order to allow a meaningful comparison of the α parameter. Furthermore, to evaluate the robustness of our data, we split our dataset both in two and in three segments, and for each of them, we calculated the sigma parameter and his variance. **Figure 1** shows the data analysis pipeline.

Goodness of Fit

We performed the goodness of fit (GOF) via the evaluation of Hellinger distance (HD) (Oosterhoff and van Zwet, 2012)

between the distribution of real data and the distribution of fitted data, in order to quantify the similarity between two probability distributions, as follows.

$$HD(P, Q) = \frac{1}{\sqrt{2}} \sqrt{\sum_{i=x_0}^{x_m} (\sqrt{P_i} - \sqrt{Q_i})^2} \quad (4)$$

where P is the probability mass function of real data, and Q is the probability mass function of fitted data.
<figure><img src="image_13.png" /><figcaption>**FIGURE 10** | Comparison of $\alpha$ parameter in broadband for time bin $\Delta t = 2$ and in delta band for time bin $\Delta t = 5$ between RTT patients (blue dots) and HS (red dots).</figcaption></figure>

The HD is in between 0 and 1; and the lower the HD, the better the fit to real data (ideally, HD = 0). Furthermore, to assess the GOF overall RTT patients and HS, we provided the empirical cumulative distribution function (ECDF) of HD.

## Statistical Analysis

We compared the distribution parameter $\alpha$ in patients and controls using permutation testing (Nichols and Holmes, 2002), for each frequency band. Specifically, at each iteration, each subject is randomly assigned to one of the two groups and then the difference between the averages of the two groups was computed. The group assignment is permuted 104 times, obtaining the null distribution of group differences, which was used to define the statistical significance of the observed difference between patients and controls. All the analyses were performed at a significance level of 0.05, using Matlab R2017a (MathWorks®) environment.

## RESULTS

As described in the Methods section, we discretized the ROI signals as a train of events, each of them representing the signal exceeding the threshold of $\pm 3$ SD.

First, we checked the critical condition in both groups, in all frequency bands, for all time bins.

We found, for both groups, a branching parameter $\sigma \approx 1$, in delta band for $\Delta t = 5$, in theta and alpha band for $\Delta t = 1$, in beta-, gamma-, and broadband for $\Delta t = 2$, as reported in **Table 1**. Moreover, to demonstrate the robustness of the estimation of the branching parameter, we estimated the variance around sigma both when we split our dataset into two segments (**Table 2**) and when we divided it into three segments (**Table 3**).

To demonstrate the strength of our data, we reported, for the chosen threshold SD = 3, the size distribution fitted with a power

law (Figure 2), duration distributions (Figure 3), and average size versus duration distributions (Figure 4) for each subject, in all time bins, for both groups separately (for brevity, we reported only delta, alpha, and broadbands). Again, we provided the effect of different thresholds (1 to 3 SD) on our data in size distributions (Figure 5), duration distribution (Figure 6), and average size versus duration distributions (Figure 7).

Furthermore, to test that our results might not be biased by spatial subsampling, we demonstrate the presence of scale time separation (Levina and Priesemann, 2017), as reported in **Figure 8**.

Last, to quantify the robustness of power law fitting, we evaluated the ECDF of HD for both groups, for all time bins in delta, alpha, and broadbands, and as reported in **Figure 9**, the highest HD is less than 0.3.

Finally, a comparison of the $\alpha$ parameter between both groups was carried out, highlighting significant differences in the delta band ($p = 0.0423$) and in broadband ($p = 0.0174$), as shown in **Figure 10**. In particular, RTT patients invariably showed a lower distribution parameter $\alpha$ in both frequency bands, for all time bins. No significant differences were found in other frequency bands.

## DISCUSSION

In this work, we set out to test the hypothesis that genetic mutations related to RTT would induce a reorganization of the whole-brain activity.

In particular, we characterized the distribution of neuronal avalanches in a group of 10 patients with Rett syndrome and 10 matched controls, using MEG, in an open-eye resting-state condition. The analysis was performed both in broadband and in the five canonical frequency bands. The study originated from recent evidence showing that the dynamics of resting-state brain activity, measured using MEG, produces scale-invariant neuronal avalanches, suggesting that the critical state is a physiological condition that is (presumably) optimal to the brain functioning. Consequently, the deviations from this optimally tuned configuration might convey the effects of the disease on large-scale coordinated activity (Meisel et al., 2012).

Our data are in line with recent findings showing that, in humans, the dynamics of resting-state brain activity is well described by a branching process where the sigma parameter is close to 1 (Meisel et al., 2013; Shriki et al., 2013).

Besides these confirmatory results, the main goal of this work is to study the frequency-specific differences in functional dynamics in RTT patients with respect to a control group. Our data showed that, around the critical state, the distribution size of neuronal avalanches, for both groups, obeys a power law. Interestingly, RTT patients show a lower value of the distribution parameter $\alpha$ both in delta and broadband, compared to controls. Furthermore, the exponents of the power law distributions vary greatly per each frequency band. This result is in line with Thompson et al., which suggests that resting-state connectivity is a frequency-dependent phenomenon (Thompson and Fransson, 2015).
The aforementioned results collectively suggest a global rearrangement in the brain functional dynamics of Rett patients. In fact, big avalanches capture the presence of widespread coordinated activations, with many different brain areas activating in a complex, meaningful sequence. The lack of large avalanches in patients might be a sign of the inability of the brain to coordinate a sufficient number of areas in sequence. One might speculate that the lack of large patterns of structured activity might be capturing the effects of the circuital changes induced by the genetic mutations present in RTT. This hypothesis might be compatible with the work by Armstrong et al., which claims that whole-brain structural changes are related to small, low-density neurons that have less dendritic branching and spine density (Armstrong et al., 1995). Furthermore, in the RTT syndrome, the dendrites of pyramidal neurons in the motor and frontal cortices are considerably shorter than in HS (Armstrong et al., 1995).

While our results are not conclusive about the brain being operating at criticality, they appear to be compatible with such hypothesis. This might be relevant given that it is known that critical dynamics facilitates optimization of processes such as information transfer and storage capacity. In particular, SOC involves a number advantages in information processing: (i) optimal information transmission because there is a perfect balance between the propagation of the signal over long distances and the resistance to saturation (Beggs and Plenz, 2003; Bertschinger and Natschläger, 2004; Kinouchi and Copelli, 2006); (ii) maximum mnemonic repertoire size, as the number of significant repetitive avalanche patterns is maximized (Haldeman and Beggs, 2005; Shew et al., 2011; Fagerholm et al., 2016); and (iii) maximized computational performance, as near the critical point, networks have better performances in different computational tasks, compared to networks with subcritical or supercritical dynamics (Bertschinger and Natschläger, 2004).

## CONCLUSION

Concluding, in our study we analyzed brain dynamics during eye-open resting state, measuring the size distribution of neuronal avalanches. We found a lower distribution parameter $\alpha$ in RTT patients in the delta frequency band and in the broadband, compared to controls. Such finding means that patients display a larger number of small avalanches and a lower number of big avalanches. These results can be explained in terms of reduced long-range coordination of neuronal activity across the brain in the pathological group. Such altered brain dynamics is likely to be suboptimal, and might underpin some of the symptomatology observed in Rett syndrome. Additionally, our results suggest that neuronal avalanches could be an innovative and advanced method to study the dynamics of brain activity. In fact, despite the rising interest in what constitutes the normal cortical dynamics in

healthy humans, the nature of the alterations induced by pathological processes remains a major question in neuroscience. Interestingly, the observation that cortical networks optimize information processing when they are in a critical state might be exploited to identify subtle, preclinical brain dysfunction. Furthermore, the study of the critical dynamics of the human brain could be enhanced by comparing normal resting avalanches with evoked states or pharmacologically modified states and assessing the sensitivity (or performing a profile analysis) of the system with respect to the branching parameter values. Finally, such approaches are inherently multivariate, taking into account the activity of all brain areas at once, hence providing solid theoretical ground to the study of emerging, holistic properties of the brain.

## DATA AVAILABILITY STATEMENT

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

## ETHICS STATEMENT

The studies involving human participants were reviewed and approved by “COMITATO ETICO UNIVERSITA' FEDERICO II” n° protocol: 392/19. The patients/participants provided their written informed consent to participate in this study.

## AUTHOR CONTRIBUTIONS

RR collected the sample, performed the MEG recordings, preprocessed the MEG data, contributed to MEG data analysis and statistical analysis, wrote the manuscript, and prepared the figures. PB collected the sample, contributed to clinical data analysis, and wrote the manuscript. AL performed the MEG recordings and preprocessed the MEG data. FB collaborated with the MEG data analysis and with the statistical analysis. MP performed the MEG recordings and preprocessed the MEG data. AP performed the MEG recordings and preprocessed the MEG data. CB collected the sample. CG, LM, and GS supervised the study. PS contributed to interpreting the results and critically revised the article. All authors read and approved the final version of the manuscript.

## FUNDING

This work has been partly supported by the University of Napoli Parthenope under the 'Bando di ricerca competitiva di ateneo per il triennio 2016-2018 progetto: implementazione di studi di connettività cerebrale mediante magnetoencefalografia in ambito multidisciplinare'.
REFERENCES

Amir, R. E., Van den Veyver, I. B., Wan, M., Tran, C. Q., Francke, U., and Zoghbi, H. Y. (1999). Rett syndrome is caused by mutations in X-linked MECP2, encoding methyl-CpG-binding protein 2. *Nat. Genet.* 23, 185–188. doi: 10.1038/13810

Armstrong, D., Dunn, J. K., Antalfy, B., and Trivedi, R. (1995). Selective dendritic alterations in the cortex of Rett syndrome. *J. Neuropathol. Exp. Neurol.* 54, 195–201. doi: 10.1097/00005072-199503000-00006

Bak, P., Tang, C., and Wiesenfeld, K. (1987). Self-organized criticality: an explanation of the 1/f noise. *Phys. Rev. Lett.* 59:381. doi: 10.1103/PhysRevLett.59.381

Beggs, J. M., and Plenz, D. (2003). Neuronal avalanches in neocortical circuits. *J. Neurosci.* 23, 11167–11177. doi: 10.1523/JNEUROSCI.23-35-11167.2003 doi: 10.1111/j.2517-6161.1995.tb02031.x

Bertschinger, N., and Natschläger, T. (2004). Real-time computation at the edge of chaos in recurrent neural networks. *Neural Comput.* 16, 1413–1436. doi: 10.1162/089976604323057443

Clauset, A., Shalizi, C. R., and Newman, M. E. J. (2009). Power-law distributions in empirical data. *SIAM Rev.* 51, 661–703. doi: 10.1137/070710111

De Carvalho, J. X., and Prado, C. P. C. (2000). Self-organized criticality in the Olami-Feder-Christensen model. *Phys. Rev. Lett.* 84, 4006–4009. doi: 10.1103/PhysRevLett.84.4006

Fagerholm, E. D., Scott, G., Shew, W. L., Song, C., Leech, R., Knöpfel, T., et al. (2016). Cortical entropy, mutual information and scale-free dynamics in waking mice. *Cereb. Cortex* 26, 3945–3952. doi: 10.1093/cercor/bhw200

Gross, J., Baillet, S., Barnes, G. R., Henson, R. N., Hillebrand, A., Jensen, O., et al. (2013). Good practice for conducting and reporting MEG research. *Neuroimage* 65, 349–363. doi: 10.1016/j.neuroimage.2012.10.001

Gutenberg, B., and Richter, C. F. (1956). Earthquake magnitude, intensity, energy, and acceleration. *Bull. Seismol. Soc. Am.* 46, 105–145.

Haldeman, C., and Beggs, J. M. (2005). Critical branching captures activity in living neural networks and maximizes the number of metastable states. *Phys. Rev. Lett.* 94:58101. doi: 10.1103/PhysRevLett.94.058101

Harris, T. E. (2002). *The Theory of Branching Processes*. Chelmsford, MA: Courier Corporation.

Jacini, F., Sorrentino, P., Lardone, A., Rucco, R., Baselice, F., Cavaliere, C., et al. (2018). Amnestic mild cognitive impairment is associated with frequency-specific brain network alterations in temporal poles. *Front. Aging Neurosci.* 10:400. doi: 10.3389/fnagi.2018.00400

Kay, S. M. (1993). *Fundamentals of Statistical Signal Processing*. Upper Saddle River, NJ: Prentice Hall.

Kinouchi, O., and Copelli, M. (2006). Optimal dynamical range of excitable networks at criticality. *Nat. Phys.* 2, 348–351. doi: 10.1038/nphys289

Lardone, A., Liparoti, M., Sorrentino, P., Rucco, R., Jacini, F., Polverino, A., et al. (2018). Mindfulness meditation is related to long-lasting changes in hippocampal functional topology during resting state: a magnetoencephalography study. *Neural Plast.* 2018:5340717. doi: 10.1155/2018/5340717

Levina, A., and Prieseemann, V. (2017). Subsampling scaling. *Nat. Commun.* 8:15140. doi: 10.1038/ncomms15140

Lombardi, F., Herrmann, H. J., Plenz, D., and de Arcangelis, L. (2016). Temporal correlations in neuronal avalanche occurrence. *Sci. Rep.* 6:24690. doi: 10.1038/srep24690

Malamud, B. D., Morein, G., and Turcotte, D. L. (1998). Forest fires: an example of self-organized critical behavior. *Science* 281, 1840–1842. doi: 10.1126/science. 281.5384.1840

Meisel, C., Olbrich, E., Shriki, O., and Achermann, P. (2013). Fading signatures of critical brain dynamics during sustained wakefulness in humans. *J. Neurosci.* 33, 17363–17372. doi: 10.1523/JNEUROSCI.1516-13.2013

Meisel, C., Storch, A., Hallmeyer-elgner, S., Bullmore, E., and Gross, T. (2012). Failure of adaptive self-organized criticality during epileptic seizure attacks. *PLoS Comput. Biol.* 8:e1002312. doi: 10.1371/journal.pcbi.1002312

Neul, J. L., Kaufmann, W. E., Glaze, D. G., Christodoulou, J., Clarke, A. J., Bahi-Buisson, N., et al. (2010). Rett syndrome: revised diagnostic criteria and nomenclature. *Ann. Neurol.* 68, 944–950. doi: 10.1002/ana.22124

Newman, M. E. J. (2007). Power laws, Pareto distributions and Zipf's law. *Contemp. Phys.* 46, 323–351. doi: 10.1080/00107510500052444

Nichols, T. E., and Holmes, A. P. (2002). Nonparametric permutation tests for functional neuroimaging: a primer with examples. *Hum. Brain Mapp.* 15, 1–25. doi: 10.1002/hbm.1058

Nolte, G. (2003). The magnetic lead field theorem in the quasi-static approximation and its use for magnetoencephalography forward calculation in realistic volume conductors. *Phys. Med. Biol.* 48, 3637–3652. doi: 10.1088/0031-9155/48/22/002

Oostenveld, R., Fries, P., Maris, E., and Schoffelen, J.-M. (2011). FieldTrip: open source software for advanced analysis of MEG, EEG, and invasive electrophysiological data. *Comput. Intell. Neurosci.* 2011:156869. doi: 10.1155/ 2011/156869

Oosterhoff, J., and van Zwet, W. R. (2012). "A note on contiguity and Hellinger distance," in *Selected Works of Willem van Zwet. Selected Works in Probability and Statistics*, eds S. van de Geer and M. Wegkamp (New York, NY: Springer), 63–72. doi: 10.1007/978-1-4614-1314-1_6

Paczuski, M., Maslov, S., and Bak, P. (1996). Avalanche dynamics in evolution, growth, and depinning models. *Phys. Rev. E Stat. Phys. Plasmas Fluids Relat. Interdiscip. Topics* 53, 414–443. doi: 10.1103/PhysRevE.53.414

Pohodich, A. E., and Zoghbi, H. Y. (2015). Rett syndrome: disruption of epigenetic control of postnatal neurological functions. *Hum. Mol. Genet.* 24, R10–R16. doi: 10.1093/hmg/ddv217

Rucco, R., Liparoti, M., Jacini, F., Baselice, F., Antenora, A., De Michele, G., et al. (2019). Mutations in the SPAST gene causing hereditary spastic paraplegia are related to global topological alterations in brain functional networks. *Neurol. Sci.* 40, 979–984. doi: 10.1007/s10072-019-3725-y

Salje, E. K. H., and Dahmen, K. A. (2014). Crackling noise in disordered materials. *Annu. Rev. Condens. Matter Phys.* 5, 233–254. doi: 10.1146/annurev- conmatphys-031113-133838

Shew, W. L., Yang, H., Yu, S., Roy, R., and Plenz, D. (2011). Information capacity and transmission are maximized in balanced cortical networks with neuronal avalanches. *J. Neurosci.* 31, 55–63. doi: 10.1523/JNEUROSCI.4637- 10.2011

Shriki, O., Alstott, J., Carver, F., Holroyd, T., Henson, R. N. A., Smith, M. L., et al. (2013). Neuronal avalanches in the resting MEG of the human brain. *J. Neurosci.* 33, 7079–7090. doi: 10.1523/JNEUROSCI.4286-12.2013

Sorrentino, P., Rucco, R., Jacini, F., Trojsi, F., Lardone, A., Fabio, B., et al. (2018). Brain functional networks become more connected as amyotrophic lateral sclerosis progresses: a source level magnetoencephalographic study. *Neuroimage Clin.* 20, 564–571. doi: 10.1016/j.nicl.2018.08.001

Sorriso, A., Sorrentino, P., Rucco, R., Mandolesi, L., Ferraioli, G., Franceschini, S., et al. (2019). An automated magnetoencephalographic data cleaning algorithm. *Comput. Methods Biomech. Biomed. Engin.* 22, 1116–1125. doi: 10.1080/ 10255842.2019.1634695

Thompson, W. H., and Fransson, P. (2015). The frequency dimension of fMRI dynamic connectivity: network connectivity, functional hubs and integration in the resting brain. *Neuroimage* 121, 227–242. doi: 10.1016/j.neuroimage.2015. 07.022

Tzourio-Mazoyer, N., Landeau, B., Papathanassiou, D., Crivello, F., Etard, O., Delcroix, N., et al. (2002). Automated anatomical labeling of activations in SPM using a macroscopic anatomical parcellation of the MNI MRI single-subject brain. *Neuroimage* 15, 273–289. doi: 10.1006/nimg.2001.0978

Van Veen, B. D., Van Drongelen, W., Yuchtman, M., and Suzuki, A. (1997). Localization of brain electrical activity via linearly constrained minimum variance spatial filtering. *IEEE Trans. Biomed. Eng.* 44, 867–880. doi: 10.1109/ 10.623056

Zapperi, S., Lauritsen, K. B., and Stanley, H. E. (1995). Self-organized branching processes: mean-field theory for avalanches. *Phys. Rev. Lett.* 75:4071. doi: 10. 1103/PhysRevLett.75.4071

**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.

Copyright © 2020 Rucco, Bernardo, Lardone, Baselice, Pesoli, Polverino, Bravaccio, Granata, Mandolesi, Sorrentino and Sorrentino. 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.