Dynamical interactions reconfigure the gradient of cortical timescales

A method for estimating dynamic functional network connectivity gradients (dFNG) from ICA captures smooth inter-network modulation

Dynamic functional network connectivity (dFNC) analysis is a widely used approach for studying brain function and offering insight into how brain networks evolve over time. Typically, dFNC studies utilized fixed spatial maps and evaluate transient changes in coupling among time courses estimated from independent component analysis (ICA). This manuscript presents a complementary approach that relaxes this assumption by spatially reordering the components dynamically at each timepoint to optimize for a smooth gradient in the FNC (i.e., a smooth gradient among ICA connectivity values). Several methods are presented to summarize dynamic FNC gradients (dFNGs) over time, starting with static FNC gradients (sFNGs), then exploring the reordering properties as well as the dynamics of the gradients themselves. We then apply this approach to a dataset of schizophrenia (SZ) patients and healthy controls (HC). Functional dysconnectivity between different brain regions has been reported in schizophrenia, yet the neural mechanisms behind it remain elusive. Using resting state fMRI and ICA on a dataset consisting of 151 schizophrenia patients and 160 age and gender-matched healthy controls, we extracted 53 intrinsic connectivity networks (ICNs) for each subject using a fully automated spatially constrained ICA approach. We develop several summaries of our functional network connectivity gradient analysis, both in a static sense, computed as the Pearson correlation coefficient between full time series, and a dynamic sense, computed using a sliding window approach followed by reordering based on the computed gradient, and evaluate group differences. Static connectivity analysis revealed significantly stronger connectivity between subcortical (SC), auditory (AUD) and visual (VIS) networks in patients, as well as hypoconnectivity in sensorimotor (SM) network relative to controls. sFNG analysis highlighted distinctive clustering patterns in patients and HCs along cognitive control (CC)/ default mode network (DMN), SC/ AUD/ SM/ cerebellar (CB), and VIS gradients. Furthermore, we observed significant differences in the sFNGs between groups in SC and CB domains. dFNG analysis suggested that SZ patients spend significantly more time in a SC/ CB state based on the first gradient, while HCs favor the DMN state. For the second gradient, however, patients exhibited significantly higher activity in CB/ VIS domains, contrasting with HCs' DMN engagement. The gradient synchrony analysis conveyed more shifts between SM/ SC networks and transmodal CC/ DMN networks in patients. In addition, the dFNG coupling revealed distinct connectivity patterns between SC, SM and CB centroids in SZ patients compared to HCs. To recap, our results advance our understanding of brain network modulation by examining smooth connectivity trajectories. This provides a more complete spatiotemporal summary of the data, contributing to the growing body of current literature regarding the functional dysconnectivity in schizophrenia patients. By employing dFNG, we highlight a new perspective to capture large scale fluctuations across the brain while maintaining the convenience of brain networks and low dimensional summary measures.

Altered spread of waves of activities at large scale is influenced by cortical thickness organization in temporal lobe epilepsy: a magnetic resonance imaging-high-density electroencephalography study

Temporal lobe epilepsy is a brain network disorder characterized by alterations at both the structural and the functional levels. It remains unclear how structure and function are related and whether this has any clinical relevance. In the present work, we adopted a novel methodological approach investigating how network structural features influence the large-scale dynamics. The functional network was defined by the spatio-temporal spreading of aperiodic bursts of activations (neuronal avalanches), as observed utilizing high-density electroencephalography in patients with temporal lobe epilepsy. The structural network was modelled as the region-based thickness covariance. Loosely speaking, we quantified the similarity of the cortical thickness of any two brain regions, both across groups and at the individual level, the latter utilizing a novel approach to define the subject-wise structural covariance network. In order to compare the structural and functional networks (at the nodal level), we studied the correlation between the probability that a wave of activity would propagate from a source to a target region and the similarity of the source region thickness as compared with other target brain regions. Building on the recent evidence that large-waves of activities pathologically spread through the epileptogenic network in temporal lobe epilepsy, also during resting state, we hypothesize that the structural cortical organization might influence such altered spatio-temporal dynamics. We observed a stable cluster of structure-function correlation in the bilateral limbic areas across subjects, highlighting group-specific features for left, right and bilateral temporal epilepsy. The involvement of contralateral areas was observed in unilateral temporal lobe epilepsy. We showed that in temporal lobe epilepsy, alterations of structural and functional networks pair in the regions where seizures propagate and are linked to disease severity. In this study, we leveraged on a well-defined model of neurological disease and pushed forward personalization approaches potentially useful in clinical practice. Finally, the methods developed here could be exploited to investigate the relationship between structure-function networks at subject level in other neurological conditions.

Topological changes of fast large-scale brain dynamics in mild cognitive impairment predict early memory impairment: a resting-state, source reconstructed, magnetoencephalography study

Functional connectivity has been used as a framework to investigate widespread brain interactions underlying cognitive deficits in mild cognitive impairment (MCI). However, many functional connectivity metrics focus on the average of the periodic activities, disregarding the aperiodic bursts of activity (i.e., the neuronal avalanches) characterizing the large-scale dynamic activities of the brain. Here, we apply the recently described avalanche transition matrix framework to source-reconstructed magnetoencephalography signals in a cohort of 32 MCI patients and 32 healthy controls to describe the spatio-temporal features of neuronal avalanches and explore their topological properties. Our results showed that MCI patients showed a more centralized network (as assessed by higher values of the degree divergence and leaf fraction) as compared to healthy controls. Furthermore, we found that the degree divergence (in the theta band) was predictive of hippocampal memory impairment. These findings highlight the role of the changes of aperiodic bursts in clinical conditions and may contribute to a more thorough phenotypical assessment of patients.

Living on the edge: network neuroscience beyond nodes

Network neuroscience has emphasized the connectional properties of neural elements - cells, populations, and regions. This has come at the expense of the anatomical and functional connections that link these elements to one another. A new perspective - namely one that emphasizes 'edges' - may prove fruitful in addressing outstanding questions in network neuroscience. We highlight one recently proposed 'edge-centric' method and review its current applications, merits, and limitations. We also seek to establish conceptual and mathematical links between this method and previously proposed approaches in the network science and neuroimaging literature. We conclude by presenting several avenues for future work to extend and refine existing edge-centric analysis.

Brain fingerprint is based on the aperiodic, scale-free, neuronal activity

Subject differentiation bears the possibility to individualize brain analyses. However, the nature of the processes generating subject-specific features remains unknown. Most of the current literature uses techniques that assume stationarity (e.g., Pearson's correlation), which might fail to capture the non-linear nature of brain activity. We hypothesize that non-linear perturbations (defined as neuronal avalanches in the context of critical dynamics) spread across the brain and carry subject-specific information, contributing the most to differentiability. To test this hypothesis, we compute the avalanche transition matrix (ATM) from source-reconstructed magnetoencephalographic data, as to characterize subject-specific fast dynamics. We perform differentiability analysis based on the ATMs, and compare the performance to that obtained using Pearson's correlation (which assumes stationarity). We demonstrate that selecting the moments and places where neuronal avalanches spread improves differentiation (P < 0.0001, permutation testing), despite the fact that most of the data (i.e., the linear part) are discarded. Our results show that the non-linear part of the brain signals carries most of the subject-specific information, thereby clarifying the nature of the processes that underlie individual differentiation. Borrowing from statistical mechanics, we provide a principled way to link emergent large-scale personalized activations to non-observable, microscopic processes.