Edge-centric analysis of stroke patients: An alternative approach for biomarkers of lesion recovery

Higher-order functional connectivity analysis of resting-state functional magnetic resonance imaging data using multivariate cumulants

Blood-level oxygenation-dependent (BOLD) functional magnetic resonance imaging (fMRI) is the most common modality to study functional connectivity in the human brain. Most research to date has focused on connectivity between pairs of brain regions. However, attention has recently turned towards connectivity involving more than two regions, that is, higher-order connectivity. It is not yet clear how higher-order connectivity can best be quantified. The measures that are currently in use cannot distinguish between pairwise (i.e., second-order) and higher-order connectivity. We show that genuine higher-order connectivity can be quantified by using multivariate cumulants. We explore the use of multivariate cumulants for quantifying higher-order connectivity and the performance of block bootstrapping for statistical inference. In particular, we formulate a generative model for fMRI signals exhibiting higher-order connectivity and use it to assess bias, standard errors, and detection probabilities. Application to resting-state fMRI data from the Human Connectome Project demonstrates that spontaneous fMRI signals are organized into higher-order networks that are distinct from second-order resting-state networks. Application to a clinical cohort of patients with multiple sclerosis further demonstrates that cumulants can be used to classify disease groups and explain behavioral variability. Hence, we present a novel framework to reliably estimate genuine higher-order connectivity in fMRI data which can be used for constructing hyperedges, and finally, which can readily be applied to fMRI data from populations with neuropsychiatric disease or cognitive neuroscientific experiments.

Edge time series components of functional connectivity and cognitive function in Alzheimer's disease

Understanding the interrelationships of brain function as measured by resting-state magnetic resonance imaging and neuropsychological/behavioral measures in Alzheimer's disease is key for advancement of neuroimaging analysis methods in clinical research. The edge time-series framework recently developed in the field of network neuroscience, in combination with other network science methods, allows for investigations of brain-behavior relationships that are not possible with conventional functional connectivity methods. Data from the Indiana Alzheimer's Disease Research Center sample (53 cognitively normal control, 47 subjective cognitive decline, 32 mild cognitive impairment, and 20 Alzheimer's disease participants) were used to investigate relationships between functional connectivity components, each derived from a subset of time points based on co-fluctuation of regional signals, and measures of domain-specific neuropsychological functions. Multiple relationships were identified with the component approach that were not found with conventional functional connectivity. These involved attentional, limbic, frontoparietal, and default mode systems and their interactions, which were shown to couple with cognitive, executive, language, and attention neuropsychological domains. Additionally, overlapping results were obtained with two different statistical strategies (network contingency correlation analysis and network-based statistics correlation). Results demonstrate that connectivity components derived from edge time-series based on co-fluctuation reveal disease-relevant relationships not observed with conventional static functional connectivity.

Brain connectivity networks underlying resting heart rate variability in acute ischemic stroke

Acute strokes can affect heart rate variability (HRV), the mechanisms how are not well understood. We included 42 acute stroke patients (2-7 days after ischemic stroke, mean age 66 years, 16 women). For analysis of HRV, 20 matched controls (mean age 60.7, 10 women) were recruited. HRV was assessed at rest, in a supine position and individual breathing rhythmus for 5 min. The coefficient of variation (VC), the root mean square of successive differences (RMSSD), the powers of low (LF, 0.04-0.14 Hz) and high (HF, 0.15-0.50 Hz) frequency bands were extracted. HRV parameters were z-transformed related to age- and sex-matched normal subjects. Z-values < -1 indicate reduced HRV. Acute stroke lesions were marked on diffusion-weighted images employing MRIcroN and co-registered to a T1-weighted structural volume-dataset. Using independent component analysis (ICA), stroke lesions were related to HRV. Subsequently, we used the ICA-derived lesion pattern as a seed and estimated the connectivity between these brain regions and seven common functional networks, which were obtained from 50 age-matched healthy subjects (mean age 68.9, 27 women). Especially, LF and VC were frequently reduced in patients. ICA revealed one covarying lesion pattern for LF and one similar for VC, predominantly affecting the right hemisphere. Activity in brain areas corresponding to these lesions mainly impact on limbic (r = 0.55 ± 0.08) and salience ventral attention networks (0.61 ± 0.10) in the group with reduced LF power (z-score < -1), but on control and default mode networks in the group with physiological LF power (z-score > -1). No different connectivity could be found for the respective VC groups. Our results suggest that HRV alteration after acute stroke might be due to affecting resting-state brain networks.

Edge Time Series Components of Functional Connectivity and Cognitive Function in Alzheimer's Disease

Understanding the interrelationships of brain function as measured by resting-state magnetic resonance imaging and neuropsychological/behavioral measures in Alzheimer's disease is key for advancement of neuroimaging analysis methods in clinical research. The edge time-series framework recently developed in the field of network neuroscience, in combination with other network science methods, allows for investigations of brain-behavior relationships that are not possible with conventional functional connectivity methods. Data from the Indiana Alzheimer's Disease Research Center sample (53 cognitively normal control, 47 subjective cognitive decline, 32 mild cognitive impairment, and 20 Alzheimer's disease participants) were used to investigate relationships between functional connectivity components, each derived from a subset of time points based on co-fluctuation of regional signals, and measures of domain-specific neuropsychological functions. Multiple relationships were identified with the component approach that were not found with conventional functional connectivity. These involved attentional, limbic, frontoparietal, and default mode systems and their interactions, which were shown to couple with cognitive, executive, language, and attention neuropsychological domains. Additionally, overlapping results were obtained with two different statistical strategies (network contingency correlation analysis and network-based statistics correlation). Results demonstrate that connectivity components derived from edge time-series based on co-fluctuation reveal disease-relevant relationships not observed with conventional static functional connectivity.

CAMBA framework: Unveiling the brain asymmetry alterations and longitudinal changes after stroke using resting-state EEG

Hemispheric asymmetry or lateralization is a fundamental principle of brain organization. However, it is poorly understood to what extent the brain asymmetries across different levels of functional organizations are evident in health or altered in brain diseases. Here, we propose a framework that integrates three degrees of brain interactions (isolated nodes, node-node, and edge-edge) into a unified analysis pipeline to capture the sliding window-based asymmetry dynamics at both the node and hemisphere levels. We apply this framework to resting-state EEG in healthy and stroke populations and investigate the stroke-induced abnormal alterations in brain asymmetries and longitudinal asymmetry changes during poststroke rehabilitation. We observe that the mean asymmetry in patients was abnormally enhanced across different frequency bands and levels of brain interactions, with these abnormal patterns strongly associated with the side of the stroke lesion. Compared to healthy controls, patients displayed significant alterations in asymmetry fluctuations, disrupting and reconfiguring the balance of inter-hemispheric integration and segregation. Additionally, analyses reveal that specific abnormal asymmetry metrics in patients tend to move towards those observed in healthy controls after short-term brain-computer interface rehabilitation. Furthermore, preliminary evidence suggests that baseline clinical and asymmetry features can predict poststroke improvements in the Fugl-Meyer assessment of the lower extremity (mean absolute error of about 2). Overall, these findings advance our understanding of hemispheric asymmetry. Our framework offers new insights into the mechanisms underlying brain alterations and recovery after a brain lesion, may help identify prognostic biomarkers, and can be easily extended to different functional modalities.

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.

A low dimensional embedding of brain dynamics enhances diagnostic accuracy and behavioral prediction in stroke

Large-scale brain networks reveal structural connections as well as functional synchronization between distinct regions of the brain. The latter, referred to as functional connectivity (FC), can be derived from neuroimaging techniques such as functional magnetic resonance imaging (fMRI). FC studies have shown that brain networks are severely disrupted by stroke. However, since FC data are usually large and high-dimensional, extracting clinically useful information from this vast amount of data is still a great challenge, and our understanding of the functional consequences of stroke remains limited. Here, we propose a dimensionality reduction approach to simplify the analysis of this complex neural data. By using autoencoders, we find a low-dimensional representation encoding the fMRI data which preserves the typical FC anomalies known to be present in stroke patients. By employing the latent representations emerging from the autoencoders, we enhanced patients' diagnostics and severity classification. Furthermore, we showed how low-dimensional representation increased the accuracy of recovery prediction.

Multiscale effects of excitatory-inhibitory homeostasis in lesioned cortical networks: A computational study

Stroke-related disruptions in functional connectivity (FC) often spread beyond lesioned areas and, given the localized nature of lesions, it is unclear how the recovery of FC is orchestrated on a global scale. Since recovery is accompanied by long-term changes in excitability, we propose excitatory-inhibitory (E-I) homeostasis as a driving mechanism. We present a large-scale model of the neocortex, with synaptic scaling of local inhibition, showing how E-I homeostasis can drive the post-lesion restoration of FC and linking it to changes in excitability. We show that functional networks could reorganize to recover disrupted modularity and small-worldness, but not network dynamics, suggesting the need to consider forms of plasticity beyond synaptic scaling of inhibition. On average, we observed widespread increases in excitability, with the emergence of complex lesion-dependent patterns related to biomarkers of relevant side effects of stroke, such as epilepsy, depression and chronic pain. In summary, our results show that the effects of E-I homeostasis extend beyond local E-I balance, driving the restoration of global properties of FC, and relating to post-stroke symptomatology. Therefore, we suggest the framework of E-I homeostasis as a relevant theoretical foundation for the study of stroke recovery and for understanding the emergence of meaningful features of FC from local dynamics.

Edge-centric functional network predicts risk propensity in economic decision-making: evidence from a resting-state fMRI study

Despite node-centric studies revealing an association between resting-state functional connectivity and individual risk propensity, the prediction of future risk decisions remains undetermined. Herein, we applied a recently emerging edge-centric method, the edge community similarity network (ECSN), to alternatively describe the community structure of resting-state brain activity and to probe its contribution to predicting risk propensity during gambling. Results demonstrated that inter-individual variability of risk decisions correlates with the inter-subnetwork couplings spanning the visual network (VN) and default mode network (DMN), cingulo-opercular task control network, and sensory/somatomotor hand network (SSHN). Particularly, participants who have higher community similarity of these subnetworks during the resting state tend to choose riskier and higher yielding bets. And in contrast to low-risk propensity participants, those who behave high-risky show stronger couplings spanning the VN and SSHN/DMN. Eventually, based on the resting-state ECSN properties, the risk rate during the gambling task is effectively predicted by the multivariable linear regression model at the individual level. These findings provide new insights into the neural substrates of the inter-individual variability in risk propensity and new neuroimaging metrics to predict individual risk decisions in advance.