Model-based whole-brain effective connectivity to study distributed cognition in health and disease

Structurally informed resting-state effective connectivity recapitulates cortical hierarchy

Interregional brain communication is mediated by the brain's physical wiring (i.e., structural connectivity). Yet, it remains unclear whether models describing directed, functional interactions between latent neuronal populations-effective connectivity-benefit from incorporating macroscale structural connectivity. Here, we assess a hierarchical empirical Bayes method: structural connectivity-based priors constrain the inversion of group-level resting-state effective connectivity, using subject-level posteriors as input; subsequently, group-level posteriors serve as empirical priors for re-evaluating subject-level effective connectivity. This approach permits knowledge of the brain's structure to inform inference of (multilevel) effective connectivity. In 17 resting-state brain networks, we find that a positive, monotonic relationship between structural connectivity and the prior probability of group-level effective connectivity generalizes across sessions and samples. Providing further validation, we show that inter-network differences in the coupling between structural and effective connectivity recapitulate a well-known unimodal-transmodal hierarchy. Thus, our results provide support for the use of our method over structurally uninformed alternatives.

Reduced Prefrontal-Thalamic Theta Flow During Working Memory Retrieval in APP/PS1 Mice

Background

Working memory deficits in Alzheimer's disease (AD) are linked to impairments in the retrieval of stored memory information. However, research on the mechanism of impaired working memory retrieval in Alzheimer's disease is still lacking.

Objective

The medial prefrontal cortex (mPFC) and mediodorsal thalamus (MD) are involved in memory retrieval. The purpose of this study is to investigate the functional interactions and information transmission between mPFC and MD in the AD model.

Methods

We recorded local field potentials from mPFC and MD while the mice (APP/PS1 transgenic model and control) performed a T-maze spatial working memory task. The temporal dynamics of oscillatory activity and bidirectional information flow between mPFC and MD were assessed during the task phases.

Results

We mainly found a significant decrease in theta flow from mPFC to MD in APP/PS1 mice during retrieval.

Conclusions

Our results indicate an important role of the mPFC-MD input for retrieval and the disrupted information transfer from mPFC to MD may be the underlying mechanism of working memory deficits in APP/PS1 mice.

Bridging functional and anatomical neural connectivity through cluster synchronization

The dynamics of the brain results from the complex interplay of several neural populations and is affected by both the individual dynamics of these areas and their connection structure. Hence, a fundamental challenge is to derive models of the brain that reproduce both structural and functional features measured experimentally. Our work combines neuroimaging data, such as dMRI, which provides information on the structure of the anatomical connectomes, and fMRI, which detects patterns of approximate synchronous activity between brain areas. We employ cluster synchronization as a tool to integrate the imaging data of a subject into a coherent model, which reconciles structural and dynamic information. By using data-driven and model-based approaches, we refine the structural connectivity matrix in agreement with experimentally observed clusters of brain areas that display coherent activity. The proposed approach leverages the assumption of homogeneous brain areas; we show the robustness of this approach when heterogeneity between the brain areas is introduced in the form of noise, parameter mismatches, and connection delays. As a proof of concept, we apply this approach to MRI data of a healthy adult at resting state.

Violations of the fluctuation-dissipation theorem reveal distinct nonequilibrium dynamics of brain states

The brain is a nonequilibrium system whose dynamics change in different brain states, such as wakefulness and deep sleep. Thermodynamics provides the tools for revealing these nonequilibrium dynamics. We used violations of the fluctuation-dissipation theorem to describe the hierarchy of nonequilibrium dynamics associated with different brain states. Together with a whole-brain model fitted to empirical human neuroimaging data, and deriving the appropriate analytical expressions, we were able to capture the deviation from equilibrium in different brain states that arises from asymmetric interactions and hierarchical organization.

A topological classifier to characterize brain states: When shape matters more than variance

Despite the remarkable accuracies attained by machine learning classifiers to separate complex datasets in a supervised fashion, most of their operation falls short to provide an informed intuition about the structure of data, and, what is more important, about the phenomena being characterized by the given datasets. By contrast, topological data analysis (TDA) is devoted to study the shape of data clouds by means of persistence descriptors and provides a quantitative characterization of specific topological features of the dataset under scrutiny. Here we introduce a novel TDA-based classifier that works on the principle of assessing quantifiable changes on topological metrics caused by the addition of new input to a subset of data. We used this classifier with a high-dimensional electro-encephalographic (EEG) dataset recorded from eleven participants during a previous decision-making experiment in which three motivational states were induced through a manipulation of social pressure. We calculated silhouettes from persistence diagrams associated with each motivated state with a ready-made band-pass filtered version of these signals, and classified unlabeled signals according to their impact on each reference silhouette. Our results show that in addition to providing accuracies within the range of those of a nearest neighbour classifier, the TDA classifier provides formal intuition of the structure of the dataset as well as an estimate of its intrinsic dimension. Towards this end, we incorporated variance-based dimensionality reduction methods to our dataset and found that in most cases the accuracy of our TDA classifier remains essentially invariant beyond a certain dimension.

Computational modelling in disorders of consciousness: Closing the gap towards personalised models for restoring consciousness

Disorders of consciousness are complex conditions characterised by persistent loss of responsiveness due to brain injury. They present diagnostic challenges and limited options for treatment, and highlight the urgent need for a more thorough understanding of how human consciousness arises from coordinated neural activity. The increasing availability of multimodal neuroimaging data has given rise to a wide range of clinically- and scientifically-motivated modelling efforts, seeking to improve data-driven stratification of patients, to identify causal mechanisms for patient pathophysiology and loss of consciousness more broadly, and to develop simulations as a means of testing in silico potential treatment avenues to restore consciousness. As a dedicated Working Group of clinicians and neuroscientists of the international Curing Coma Campaign, here we provide our framework and vision to understand the diverse statistical and generative computational modelling approaches that are being employed in this fast-growing field. We identify the gaps that exist between the current state-of-the-art in statistical and biophysical computational modelling in human neuroscience, and the aspirational goal of a mature field of modelling disorders of consciousness; which might drive improved treatments and outcomes in the clinic. Finally, we make several recommendations for how the field as a whole can work together to address these challenges.

Exploring neural heterogeneity in inattention and hyperactivity

Inattention and hyperactivity are cardinal symptoms of Attention Deficit Hyperactivity Disorder (ADHD). These characteristics have also been observed across a range of other neurodevelopmental conditions, such as autism and dyspraxia, suggesting that they might best be studied across diagnostic categories. Here, we evaluated the associations between inattention and hyperactivity behaviours and features of the structural brain network (connectome) in a large transdiagnostic sample of children (Centre for Attention, Learning, and Memory; n = 383). In our sample, we found that a single latent factor explains 77.6% of variance in scores across multiple questionnaires measuring inattention and hyperactivity. Partial Least-Squares (PLS) regression revealed that variability in this latent factor could not be explained by a linear component representing nodewise properties of connectomes. We then investigated the type and extent of neural heterogeneity in a subset of our sample with clinically-elevated levels of inattention and hyperactivity. Multidimensional scaling combined with k-means clustering revealed two neural subtypes in children with elevated levels of inattention and hyperactivity (n = 232), differentiated primarily by nodal communicability-a measure which demarcates the extent to which neural signals propagate through specific brain regions. These different clusters had similar behavioural profiles, which included high levels of inattention and hyperactivity. However, one of the clusters scored higher on multiple cognitive assessment measures of executive function. We conclude that inattention and hyperactivity are so common in children with neurodevelopmental difficulties because they emerge through multiple different trajectories of brain development. In our own data, we can identify two of these possible trajectories, which are reflected by measures of structural brain network topology and cognition.

From Neurotransmitters to Networks: Transcending Organisational Hierarchies with Molecular-informed Functional Imaging

The human brain exhibits complex interactions across micro, meso-, and macro-scale organisational principles. Recent synergistic multi-modal approaches have begun to link micro-scale information to systems level dynamics, transcending organisational hierarchies and offering novel perspectives into the brain's function and dysfunction. Specifically, the distribution of micro-scale properties (such as receptor density or gene expression) can be mapped onto macro-scale measures from functional MRI to provide novel neurobiological insights. Methodological approaches to enrich functional imaging analyses with molecular information are rapidly evolving, with several streams of research having developed relatively independently, each offering unique potential to explore the trans-hierarchical functioning of the brain. Here, we address the three principal streams of research - spatial correlation, molecular-enriched network, and in-silico whole brain modelling analyses - to provide a critical overview of the different sources of molecular information, how this information can be utilised within analyses of fMRI data, the merits and pitfalls of each methodology, and, through the use of key examples, highlight their promise to shed new light on key domains of neuroscientific inquiry.

A Survey on Brain Effective Connectivity Network Learning

Human brain effective connectivity characterizes the causal effects of neural activities among different brain regions. Studies of brain effective connectivity networks (ECNs) for different populations contribute significantly to the understanding of the pathological mechanism associated with neuropsychiatric diseases and facilitate finding new brain network imaging markers for the early diagnosis and evaluation for the treatment of cerebral diseases. A deeper understanding of brain ECNs also greatly promotes brain-inspired artificial intelligence (AI) research in the context of brain-like neural networks and machine learning. Thus, how to picture and grasp deeper features of brain ECNs from functional magnetic resonance imaging (fMRI) data is currently an important and active research area of the human brain connectome. In this survey, we first show some typical applications and analyze existing challenging problems in learning brain ECNs from fMRI data. Second, we give a taxonomy of ECN learning methods from the perspective of computational science and describe some representative methods in each category. Third, we summarize commonly used evaluation metrics and conduct a performance comparison of several typical algorithms both on simulated and real datasets. Finally, we present the prospects and references for researchers engaged in learning ECNs.

Effective connectivity analysis of brain networks of mathematically gifted adolescents using transfer entropy

Using functional neuroimaging, electrophysiological techniques and neural data processing techniques, neuroscientists have found that mathematically gifted adolescents exhibit unusual neurocognitive features in the activation of task-related brain regions. Hemispheric information interaction, functional reorganization of networks, and utilization of task-related brain regions are beneficial to rapid and efficient task processing. Based on Granger causality channel selection, the transfer entropy (TE) value between effective channels was computed, and the information flow patterns in the directed functional brain networks derived from electroencephalography (EEG) data during deductive reasoning tasks were explored. We evaluated the workspace configuration patterns of the brain network and the global integration characteristics of separated brain regions using node strength, motif, directed clustering coefficient and characteristic path length in the brain networks of mathematically gifted adolescents with effective connectivity. The empirical results demonstrated that a more integrated functional network at the global level and a more efficient clique at the local level support a pattern of workspace configuration in the mathematically gifted brain that is more conducive to task-related information processing.

Characterization of regional differences in resting-state fMRI with a data-driven network model of brain dynamics

Model-based data analysis of whole-brain dynamics links the observed data to model parameters in a network of neural masses. Recently, studies focused on the role of regional variance of model parameters. Such analyses however necessarily depend on the properties of preselected neural mass model. We introduce a method to infer from the functional data both the neural mass model representing the regional dynamics and the region- and subject-specific parameters while respecting the known network structure. We apply the method to human resting-state fMRI. We find that the underlying dynamics can be described as noisy fluctuations around a single fixed point. The method reliably discovers three regional parameters with clear and distinct role in the dynamics, one of which is strongly correlated with the first principal component of the gene expression spatial map. The present approach opens a novel way to the analysis of resting-state fMRI with possible applications for understanding the brain dynamics during aging or neurodegeneration.

Amortized Bayesian inference on generative dynamical network models of epilepsy using deep neural density estimators

Whole-brain modeling of epilepsy combines personalized anatomical data with dynamical models of abnormal activities to generate spatio-temporal seizure patterns as observed in brain imaging data. Such a parametric simulator is equipped with a stochastic generative process, which itself provides the basis for inference and prediction of the local and global brain dynamics affected by disorders. However, the calculation of likelihood function at whole-brain scale is often intractable. Thus, likelihood-free algorithms are required to efficiently estimate the parameters pertaining to the hypothetical areas, ideally including the uncertainty. In this study, we introduce the simulation-based inference for the virtual epileptic patient model (SBI-VEP), enabling us to amortize the approximate posterior of the generative process from a low-dimensional representation of whole-brain epileptic patterns. The state-of-the-art deep learning algorithms for conditional density estimation are used to readily retrieve the statistical relationships between parameters and observations through a sequence of invertible transformations. We show that the SBI-VEP is able to efficiently estimate the posterior distribution of parameters linked to the extent of the epileptogenic and propagation zones from sparse intracranial electroencephalography recordings. The presented Bayesian methodology can deal with non-linear latent dynamics and parameter degeneracy, paving the way for fast and reliable inference on brain disorders from neuroimaging modalities.

Entropy production of multivariate Ornstein-Uhlenbeck processes correlates with consciousness levels in the human brain

Consciousness is supported by complex patterns of brain activity which are indicative of irreversible nonequilibrium dynamics. While the framework of stochastic thermodynamics has facilitated the understanding of physical systems of this kind, its application to infer the level of consciousness from empirical data remains elusive. We faced this challenge by calculating entropy production in a multivariate Ornstein-Uhlenbeck process fitted to Functional magnetic resonance imaging brain activity recordings. To test this approach, we focused on the transition from wakefulness to deep sleep, revealing a monotonous relationship between entropy production and the level of consciousness. Our results constitute robust signatures of consciousness while also advancing our understanding of the link between consciousness and complexity from the fundamental perspective of statistical physics.

Dynamic sensitivity analysis: Defining personalised strategies to drive brain state transitions via whole brain modelling

Traditionally, in neuroimaging, model-free analyses are used to find significant differences between brain states via signal detection theory. Depending on the a priori assumptions about the underlying data, different spatio-temporal features can be analysed. Alternatively, model-based techniques infer features from the data and compare significance from model parameters. However, to assess transitions from one brain state to another remains a challenge in current paradigms. Here, we introduce a "Dynamic Sensitivity Analysis" framework that quantifies transitions between brain states in terms of stimulation ability to rebalance spatio-temporal brain activity towards a target state such as healthy brain dynamics. In practice, it means building a whole-brain model fitted to the spatio-temporal description of brain dynamics, and applying systematic stimulations in-silico to assess the optimal strategy to drive brain dynamics towards a target state. Further, we show how Dynamic Sensitivity Analysis extends to various brain stimulation paradigms, ultimately contributing to improving the efficacy of personalised clinical interventions.

Directional prefrontal-thalamic information flow is selectively required during spatial working memory retrieval

Introduction

Spatial working memory is a kind of short-term memory that allows temporarily storing and manipulating spatial information. Evidence suggests that spatial working memory is processed through three distinctive phases: Encoding, maintenance, and retrieval. Though the medial prefrontal cortex (mPFC) and mediodorsal thalamus (MD) are involved in memory retrieval, how the functional interactions and information transfer between mPFC and MD remains largely unclear.

Methods

We recorded local field potentials (LFPs) from mPFC and MD while mice performed a spatial working memory task in T-maze. The temporal dynamics of functional interactions and bidirectional information flow between mPFC and MD was quantitatively assessed by using directed transfer function.

Results

Our results showed a significantly elevated information flow from mPFC to MD, varied in time and frequency (theta in particular), accompanying successful memory retrieval.

Discussion

Elevated theta information flow, a feature that was absent on error trials, indicates an important role of the directional information transfer from mPFC to MD for memory retrieval.

Advancing brain network models to reconcile functional neuroimaging and clinical research

Functional magnetic resonance imaging (fMRI) captures information on brain function beyond the anatomical alterations that are traditionally visually examined by neuroradiologists. However, the fMRI signals are complex in addition to being noisy, so fMRI still faces limitations for clinical applications. Here we review methods that have been proposed as potential solutions so far, namely statistical, biophysical and decoding models, with their strengths and weaknesses. We especially evaluate the ability of these models to directly predict clinical variables from their parameters (predictability) and to extract clinically relevant information regarding biological mechanisms and relevant features for classification and prediction (interpretability). We then provide guidelines for useful applications and pitfalls of such fMRI-based models in a clinical research context, looking beyond the current state of the art. In particular, we argue that the clinical relevance of fMRI calls for a new generation of models for fMRI data, which combine the strengths of both biophysical and decoding models. This leads to reliable and biologically meaningful model parameters, which thus fulfills the need for simultaneous interpretability and predictability. In our view, this synergy is fundamental for the discovery of new pharmacological and interventional targets, as well as the use of models as biomarkers in neurology and psychiatry.

From descriptive connectome to mechanistic connectome: Generative modeling in functional magnetic resonance imaging analysis

As a newly emerging field, connectomics has greatly advanced our understanding of the wiring diagram and organizational features of the human brain. Generative modeling-based connectome analysis, in particular, plays a vital role in deciphering the neural mechanisms of cognitive functions in health and dysfunction in diseases. Here we review the foundation and development of major generative modeling approaches for functional magnetic resonance imaging (fMRI) and survey their applications to cognitive or clinical neuroscience problems. We argue that conventional structural and functional connectivity (FC) analysis alone is not sufficient to reveal the complex circuit interactions underlying observed neuroimaging data and should be supplemented with generative modeling-based effective connectivity and simulation, a fruitful practice that we term "mechanistic connectome." The transformation from descriptive connectome to mechanistic connectome will open up promising avenues to gain mechanistic insights into the delicate operating principles of the human brain and their potential impairments in diseases, which facilitates the development of effective personalized treatments to curb neurological and psychiatric disorders.

Driving brain state transitions in major depressive disorder through external stimulation

Major depressive disorder (MDD) as a dysfunction of neural circuits and brain networks has been established in modern neuroimaging sciences. However, the brain state transitions between MDD and health through external stimulation remain unclear, which limits translation to clinical contexts and demonstrable clinical utility. We propose a framework of the large-scale whole-brain network model for MDD linking the underlying anatomical connectivity with functional dynamics obtained from diffusion tensor imaging (DTI) and functional magnetic resonance imaging (fMRI). Then, we further explored the optimal brain regions to promote the transition of brain states between MDD and health through external stimulation of the model. Based on the whole-brain model successfully fitting the brain state space in MDD and the health, we demonstrated that the transition from MDD to health is achieved by the excitatory activation of the limbic system and from health to MDD by the inhibitory stimulation of the reward circuit. Our finding provides novel biophysical evidence for the neural mechanism of MDD and its recovery and allows the discovery of new stimulation targets for MDD recovery.

One session of fMRI-Neurofeedback training on motor imagery modulates whole-brain effective connectivity and dynamical complexity

In the past decade, several studies have shown that Neurofeedback (NFB) by functional magnetic resonance imaging can alter the functional coupling of targeted and non-targeted areas. However, the causal mechanisms underlying these changes remain uncertain. Here, we applied a whole-brain dynamical model to estimate Effective Connectivity (EC) profiles of resting-state data acquired before and immediately after a single-session NFB training for 17 participants who underwent motor imagery NFB training and 16 healthy controls who received sham feedback. Within-group and between-group classification analyses revealed that only for the NFB group it was possible to accurately discriminate between the 2 resting-state sessions. NFB training-related signatures were reflected in a support network of direct connections between areas involved in reward processing and implicit learning, together with regions belonging to the somatomotor, control, attention, and default mode networks, identified through a recursive-feature elimination procedure. By applying a data-driven approach to explore NFB-induced changes in spatiotemporal dynamics, we demonstrated that these regions also showed decreased switching between different brain states (i.e. metastability) only following real NFB training. Overall, our findings contribute to the understanding of NFB impact on the whole brain’s structure and function by shedding light on the direct connections between brain areas affected by NFB training.

Meditation-induced effects on whole-brain structural and effective connectivity

In the past decades, there has been a growing scientific interest in characterizing neural correlates of meditation training. Nonetheless, the mechanisms underlying meditation remain elusive. In the present work, we investigated meditation-related changes in functional dynamics and structural connectivity (SC). For this purpose, we scanned experienced meditators and control (naive) subjects using magnetic resonance imaging (MRI) to acquire structural and functional data during two conditions, resting-state and meditation (focused attention on breathing). In this way, we aimed to characterize and distinguish both short-term and long-term modifications in the brain’s structure and function. First, to analyze the fMRI data, we calculated whole-brain effective connectivity (EC) estimates, relying on a dynamical network model to replicate BOLD signals’ spatio-temporal structure, akin to functional connectivity (FC) with lagged correlations. We compared the estimated EC, FC, and SC links as features to train classifiers to predict behavioral conditions and group identity. Then, we performed a network-based analysis of anatomical connectivity. We demonstrated through a machine-learning approach that EC features were more informative than FC and SC solely. We showed that the most informative EC links that discriminated between meditators and controls involved several large-scale networks mainly within the left hemisphere. Moreover, we found that differences in the functional domain were reflected to a smaller extent in changes at the anatomical level as well. The network-based analysis of anatomical pathways revealed strengthened connectivity for meditators compared to controls between four areas in the left hemisphere belonging to the somatomotor, dorsal attention, subcortical and visual networks. Overall, the results of our whole-brain model-based approach revealed a mechanism underlying meditation by providing causal relationships at the structure-function level.

Directed Brain Connectivity Identifies Widespread Functional Network Abnormalities in Parkinson's Disease

Parkinson's disease (PD) is a neurodegenerative disorder characterized by topological abnormalities in large-scale functional brain networks, which are commonly analyzed using undirected correlations in the activation signals between brain regions. This approach assumes simultaneous activation of brain regions, despite previous evidence showing that brain activation entails causality, with signals being typically generated in one region and then propagated to other ones. To address this limitation, here, we developed a new method to assess whole-brain directed functional connectivity in participants with PD and healthy controls using antisymmetric delayed correlations, which capture better this underlying causality. Our results show that whole-brain directed connectivity, computed on functional magnetic resonance imaging data, identifies widespread differences in the functional networks of PD participants compared with controls, in contrast to undirected methods. These differences are characterized by increased global efficiency, clustering, and transitivity combined with lower modularity. Moreover, directed connectivity patterns in the precuneus, thalamus, and cerebellum were associated with motor, executive, and memory deficits in PD participants. Altogether, these findings suggest that directional brain connectivity is more sensitive to functional network differences occurring in PD compared with standard methods, opening new opportunities for brain connectivity analysis and development of new markers to track PD progression.

Effective connectivity extracts clinically relevant prognostic information from resting state activity in stroke

Recent resting-state functional MRI studies in stroke patients have identified two robust biomarkers of acute brain dysfunction: a reduction of inter-hemispheric functional connectivity between homotopic regions of the same network, and an abnormal increase of ipsi-lesional functional connectivity between task-negative and task-positive resting-state networks. Whole-brain computational modelling studies, at the individual subject level, using undirected effective connectivity derived from empirically measured functional connectivity, have shown a reduction of measures of integration and segregation in stroke as compared to healthy brains. Here we employ a novel method, first, to infer whole-brain directional effective connectivity from zero-lagged and lagged covariance matrices, then, to compare it to empirically measured functional connectivity for predicting stroke versus healthy status, and patient performance (zero, one, multiple deficits) across neuropsychological tests. We also investigated the accuracy of functional connectivity versus model effective connectivity in predicting the long-term outcome from acute measures. Both functional and effective connectivity predicted healthy from stroke individuals significantly better than the chance-level; however, accuracy for the effective connectivity was significantly higher than for functional connectivity at 1- to 2-week, 3-month and 1-year post-stroke. Predictive functional connections mainly included those reported in previous studies (within-network inter-hemispheric and between task-positive and -negative networks intra-hemispherically). Predictive effective connections included additional between-network links. Effective connectivity was a better predictor than functional connectivity of the number of behavioural domains in which patients suffered deficits, both at 2-week and 1-year post-onset of stroke. Interestingly, patient deficits at 1-year time-point were better predicted by effective connectivity values at 2 weeks rather than at 1-year time-point. Our results thus demonstrate that the second-order statistics of functional MRI resting-state activity at an early stage of stroke, derived from a whole-brain effective connectivity, estimated in a model fitted to reproduce the propagation of neuronal activity, has pertinent information for clinical prognosis.

Crosslink: An R Package for Network Visualization of Grouped Nodes

The demand for network visualization of relationships between nodes attributed to different categories grows in various biomedical research scenarios, such as gene regulatory networks, drug-target networks, ligand-receptor interactions and association networks of multi-omics elements. Elegantly visualizing the relationships between nodes with complex metadata of nodes and edges appended may inspire new insights. Here, we developed the crosslink R package, tailored for network visualization of grouped nodes, to provide a series of flexible functions for generating network diagrams. We first designed a CrossLink class for storage of metadata about nodes and edges and manipulation of node coordinates. Then affine transformation and function mapping transformation are implemented to perform fundamental node coordinates transformation by groups, based on which various network layouts can be defined easily. For convenience, we predefined several commonly used layouts, including row, column, arc, polygon and hive, which also can be combined in one layout. Finally, we designed a user-friendly wrapper function to draw network connections, aesthetic mappings of metadata and decoration with related annotation graphs in one interface by taking advantage of the powerful ggplot2 system. Overall, the crosslink R package is easy-to-use for achieving complex visualization of a network diagram of grouped nodes surrounded by associated annotation graphs.

Multiscale neural modeling of resting-state fMRI reveals executive-limbic malfunction as a core mechanism in major depressive disorder

Major depressive disorder (MDD) represents a grand challenge to human health and society, but the underlying pathophysiological mechanisms remain elusive. Previous neuroimaging studies have suggested that MDD is associated with abnormal interactions and dynamics in two major neural systems including the default mode - salience (DMN-SAL) network and the executive - limbic (EXE-LIM) network, but it is not clear which network plays a central role and which network plays a subordinate role in MDD pathophysiology. To address this question, we refined a newly developed Multiscale Neural Model Inversion (MNMI) framework and applied it to test whether MDD is more affected by impaired circuit interactions in the DMN-SAL network or the EXE-LIM network. The model estimates the directed connection strengths between different neural populations both within and between brain regions based on resting-state fMRI data collected from normal healthy subjects and patients with MDD. Results show that MDD is primarily characterized by abnormal circuit interactions in the EXE-LIM network rather than the DMN-SAL network. Specifically, we observe reduced frontoparietal effective connectivity that potentially contributes to hypoactivity in the dorsolateral prefrontal cortex (dlPFC), and decreased intrinsic inhibition combined with increased excitation from the superior parietal cortex (SPC) that potentially lead to amygdala hyperactivity, together resulting in activation imbalance in the PFC-amygdala circuit that pervades in MDD. Moreover, the model reveals reduced PFC-to-hippocampus excitation but decreased SPC-to-thalamus inhibition in MDD population that potentially lead to hypoactivity in the hippocampus and hyperactivity in the thalamus, consistent with previous experimental data. Overall, our findings provide strong support for the long-standing limbic-cortical dysregulation model in major depression but also offer novel insights into the multiscale pathophysiology of this debilitating disease.

Brain States and Transitions: Insights from Computational Neuroscience

Within the field of computational neuroscience there are great expectations of finding new ways to rebalance the complex dynamic system of the human brain through controlled pharmacological or electromagnetic perturbation. Yet many obstacles remain between the ability to accurately predict how and where best to perturb to force a transition from one brain state to another. The foremost challenge is a commonly agreed definition of a given brain state. Recent progress in computational neuroscience has made it possible to robustly define brain states and force transitions between them. Here, we review the state of the art and propose a framework for determining the functional hierarchical organization describing any given brain state. We describe the latest advances in creating sophisticated whole-brain computational models with interacting neuronal and neurotransmitter systems that can be studied fully in silico to predict and design novel pharmacological and electromagnetic interventions to rebalance them in disease.

Sparse coupled logistic regression to estimate co-activation and modulatory influences of brain regions

Accurate mapping of the functional interactions between remote brain areas with resting-state functional magnetic resonance imaging requires the quantification of their underlying dynamics. In conventional methodological pipelines, a spatial scale of interest is first selected, and dynamic analysis then proceeds at this hypothesised level of complexity. If large-scale functional networks or states are studied, more local regional rearrangements are then not described, potentially missing important neurobiological information. Here, we propose a novel mathematical framework that jointly estimates resting-state functional networks, and spatially more localised cross-regional modulations. To do so, the changes in activity of each brain region are modelled by a logistic regression including co-activation coefficients (reflective of network assignment, as they highlight simultaneous activations across areas) and causal interplays (denoting finer regional cross-talks, when one region active at timetmodulates thettot+1 transition likelihood of another area). A two-parameter L1 regularisation scheme is used to make these two sets of coefficients sparse: one controls overall sparsity, while the other governs the trade-off between co-activations and causal interplays, enabling to properly fit the data despite the yet unknown balance between both types of couplings. Across a range of simulation settings, we show that the framework successfully retrieves the two types of cross-regional interactions at once. Performance across noise and sample size settings was globally on par with that of other existing methods, with the potential to reveal more precise information missed by alternative approaches. Preliminary application to experimental data revealed that in the resting brain, co-activations and causal modulations co-exist with a varying balance across regions. Our methodological pipeline offers a conceptually elegant alternative for the assessment of functional brain dynamics, and can be downloaded at https://c4science.ch/source/Sparse_logistic_regression.git.

Stable biomarker identification for predicting schizophrenia in the human connectome

Schizophrenia, as a psychiatric disorder, has recognized brain alterations both at the structural and at the functional magnetic resonance imaging level. The developing field of connectomics has attracted much attention as it allows researchers to take advantage of powerful tools of network analysis in order to study structural and functional connectivity abnormalities in schizophrenia. Many methods have been proposed to identify biomarkers in schizophrenia, focusing mainly on improving the classification performance or performing statistical comparisons between groups. However, the stability of biomarkers selection has been for long overlooked in the connectomics field. In this study, we follow a machine learning approach where the identification of biomarkers is addressed as a feature selection problem for a classification task. We perform a recursive feature elimination and support vector machines (RFE-SVM) approach to identify the most meaningful biomarkers from the structural, functional, and multi-modal connectomes of healthy controls and patients. Furthermore, the stability of the retrieved biomarkers is assessed across different subsamplings of the dataset, allowing us to identify the affected core of the pathology. Considering our technique altogether, it demonstrates a principled way to achieve both accurate and stable biomarkers while highlighting the importance of multi-modal approaches to brain pathology as they tend to reveal complementary information.

Deep Temporal Organization of fMRI Phase Synchrony Modes Promotes Large-Scale Disconnection in Schizophrenia

Itinerant dynamics of the brain generates transient and recurrent spatiotemporal patterns in neuroimaging data. Characterizing metastable functional connectivity (FC) - particularly at rest and using functional magnetic resonance imaging (fMRI) - has shaped the field of dynamic functional connectivity (DFC). Mainstream DFC research relies on (sliding window) correlations to identify recurrent FC patterns. Recently, functional relevance of the instantaneous phase synchrony (IPS) of fMRI signals has been revealed using imaging studies and computational models. In the present paper, we identify the repertoire of whole-brain inter-network IPS states at rest. Moreover, we uncover a hierarchy in the temporal organization of IPS modes. We hypothesize that connectivity disorder in schizophrenia (SZ) is related to the (deep) temporal arrangement of large-scale IPS modes. Hence, we analyze resting-state fMRI data from 68 healthy controls (HC) and 51 SZ patients. Seven resting-state networks (and their sub-components) are identified using spatial independent component analysis. IPS is computed between subject-specific network time courses, using analytic signals. The resultant phase coupling patterns, across time and subjects, are clustered into eight IPS states. Statistical tests show that the relative expression and mean lifetime of certain IPS states have been altered in SZ. Namely, patients spend (45%) less time in a globally coherent state and a subcortical-centered state, and (40%) more time in states reflecting anticoupling within the cognitive control network, compared to the HC. Moreover, the transition profile (between states) reveals a deep temporal structure, shaping two metastates with distinct phase synchrony profiles. A metastate is a collection of states such that within-metastate transitions are more probable than across. Remarkably, metastate occupation balance is altered in SZ, in favor of the less synchronous metastate that promotes disconnection across networks. Furthermore, the trajectory of IPS patterns is less efficient, less smooth, and more restricted in SZ subjects, compared to the HC. Finally, a regression analysis confirms the diagnostic value of the defined IPS measures for SZ identification, highlighting the distinctive role of metastate proportion. Our results suggest that the proposed IPS features may be used for classification studies and for characterizing phase synchrony modes in other (clinical) populations.