A generative model of whole-brain effective connectivity

Atypical connectome topography and signal flow in temporal lobe epilepsy

Temporal lobe epilepsy (TLE) is the most common pharmaco-resistant epilepsy in adults. While primarily associated with mesiotemporal pathology, recent evidence suggests that brain alterations in TLE extend beyond the paralimbic epicenter and impact macroscale function and cognitive functions, particularly memory. Using connectome-wide manifold learning and generative models of effective connectivity, we examined functional topography and directional signal flow patterns between large-scale neural circuits in TLE at rest. Studying a multisite cohort of 95 patients with TLE and 95 healthy controls, we observed atypical functional topographies in the former group, characterized by reduced differentiation between sensory and transmodal association cortices, with most marked effects in bilateral temporo-limbic and ventromedial prefrontal cortices. These findings were consistent across all study sites, present in left and right lateralized patients, and validated in a subgroup of patients with histopathological validation of mesiotemporal sclerosis and post-surgical seizure freedom. Moreover, they were replicated in an independent cohort of 30 TLE patients and 40 healthy controls. Further analyses demonstrated that reduced differentiation related to decreased functional signal flow into and out of temporolimbic cortical systems and other brain networks. Parallel analyses of structural and diffusion-weighted MRI data revealed that topographic alterations were independent of TLE-related cortical thinning but partially mediated by white matter microstructural changes that radiated away from paralimbic circuits. Finally, we found a strong association between the degree of functional alterations and behavioral markers of memory dysfunction. Our work illustrates the complex landscape of macroscale functional imbalances in TLE, which can serve as intermediate markers bridging microstructural changes and cognitive impairment.

Linking fast and slow: The case for generative models

A pervasive challenge in neuroscience is testing whether neuronal connectivity changes over time due to specific causes, such as stimuli, events, or clinical interventions. Recent hardware innovations and falling data storage costs enable longer, more naturalistic neuronal recordings. The implicit opportunity for understanding the self-organised brain calls for new analysis methods that link temporal scales: from the order of milliseconds over which neuronal dynamics evolve, to the order of minutes, days, or even years over which experimental observations unfold. This review article demonstrates how hierarchical generative models and Bayesian inference help to characterise neuronal activity across different time scales. Crucially, these methods go beyond describing statistical associations among observations and enable inference about underlying mechanisms. We offer an overview of fundamental concepts in state-space modeling and suggest a taxonomy for these methods. Additionally, we introduce key mathematical principles that underscore a separation of temporal scales, such as the slaving principle, and review Bayesian methods that are being used to test hypotheses about the brain with multiscale data. We hope that this review will serve as a useful primer for experimental and computational neuroscientists on the state of the art and current directions of travel in the complex systems modelling literature.

A recurrent Hopfield network for estimating meso-scale effective connectivity in MEG

The architecture of communication within the brain, represented by the human connectome, has gained a paramount role in the neuroscience community. Several features of this communication, e.g., the frequency content, spatial topology, and temporal dynamics are currently well established. However, identifying generative models providing the underlying patterns of inhibition/excitation is very challenging. To address this issue, we present a novel generative model to estimate large-scale effective connectivity from MEG. The dynamic evolution of this model is determined by a recurrent Hopfield neural network with asymmetric connections, and thus denoted Recurrent Hopfield Mass Model (RHoMM). Since RHoMM must be applied to binary neurons, it is suitable for analyzing Band Limited Power (BLP) dynamics following a binarization process. We trained RHoMM to predict the MEG dynamics through a gradient descent minimization and we validated it in two steps. First, we showed a significant agreement between the similarity of the effective connectivity patterns and that of the interregional BLP correlation, demonstrating RHoMM's ability to capture individual variability of BLP dynamics. Second, we showed that the simulated BLP correlation connectomes, obtained from RHoMM evolutions of BLP, preserved some important topological features, e.g, the centrality of the real data, assuring the reliability of RHoMM. Compared to other biophysical models, RHoMM is based on recurrent Hopfield neural networks, thus, it has the advantage of being data-driven, less demanding in terms of hyperparameters and scalable to encompass large-scale system interactions. These features are promising for investigating the dynamics of inhibition/excitation at different spatial scales.

Controlling target brain regions by optimal selection of input nodes

The network control theory framework holds great potential to inform neurostimulation experiments aimed at inducing desired activity states in the brain. However, the current applicability of the framework is limited by inappropriate modeling of brain dynamics, and an overly ambitious focus on whole-brain activity control. In this work, we leverage recent progress in linear modeling of brain dynamics (effective connectivity) and we exploit the concept of target controllability to focus on the control of a single region or a small subnetwork of nodes. We discuss when control may be possible with a reasonably low energy cost and few stimulation loci, and give general predictions on where to stimulate depending on the subset of regions one wishes to control. Importantly, using the robustly asymmetric effective connectome instead of the symmetric structural connectome (as in previous research), we highlight the fundamentally different roles in- and out-hubs have in the control problem, and the relevance of inhibitory connections. The large degree of inter-individual variation in the effective connectome implies that the control problem is best formulated at the individual level, but we discuss to what extent group results may still prove useful.

Aberrant Effective Connectivity During Eye Gaze Processing Is Linked to Social Functioning and Symptoms in Schizophrenia

BACKGROUND

Patients with schizophrenia show abnormal gaze processing, which is associated with social dysfunction. These abnormalities are related to aberrant connectivity among brain regions that are associated with visual processing, social cognition, and cognitive control. In this study, we investigated 1) how effective connectivity during gaze processing is disrupted in schizophrenia and 2) how this may contribute to social dysfunction and clinical symptoms.

METHODS

Thirty-nine patients with schizophrenia/schizoaffective disorder (SZ) and 33 healthy control participants completed an eye gaze processing task during functional magnetic resonance imaging. Participants viewed faces with different gaze angles and performed explicit and implicit gaze processing. Four brain regions-the secondary visual cortex, posterior superior temporal sulcus, inferior parietal lobule, and posterior medial frontal cortex-were identified as nodes for dynamic causal modeling analysis.

RESULTS

Both the SZ and healthy control groups showed similar model structures for general gaze processing. Explicit gaze discrimination led to changes in effective connectivity, including stronger excitatory, bottom-up connections from the secondary visual cortex to the posterior superior temporal sulcus and inferior parietal lobule and inhibitory, top-down connections from the posterior medial frontal cortex to the secondary visual cortex. Group differences in top-down modulation from the posterior medial frontal cortex to the posterior superior temporal sulcus and inferior parietal lobule were noted, such that these inhibitory connections were attenuated in the healthy control group but further strengthened in the SZ group. Connectivity was associated with social dysfunction and symptom severity.

CONCLUSIONS

The SZ group showed notably stronger top-down inhibition during explicit gaze discrimination, which was associated with more social dysfunction but less severe symptoms among patients. These findings help pinpoint neural mechanisms of aberrant gaze processing and may serve as future targets for interventions that combine neuromodulation with social cognitive training.

The regional variation of laminar thickness in the human isocortex is related to cortical hierarchy and interregional connectivity

The human isocortex consists of tangentially organized layers with unique cytoarchitectural properties. These layers show spatial variations in thickness and cytoarchitecture across the neocortex, which is thought to support function through enabling targeted corticocortical connections. Here, leveraging maps of the 6 cortical layers based on 3D human brain histology, we aimed to quantitatively characterize the systematic covariation of laminar structure in the cortex and its functional consequences. After correcting for the effect of cortical curvature, we identified a spatial pattern of changes in laminar thickness covariance from lateral frontal to posterior occipital regions, which differentiated the dominance of infra- versus supragranular layer thickness. Corresponding to the laminar regularities of cortical connections along cortical hierarchy, the infragranular-dominant pattern of laminar thickness was associated with higher hierarchical positions of regions, mapped based on resting-state effective connectivity in humans and tract-tracing of structural connections in macaques. Moreover, we show that regions with similar laminar thickness patterns have a higher likelihood of structural connections and strength of functional connections. In sum, here we characterize the organization of laminar thickness in the human isocortex and its association with cortico-cortical connectivity, illustrating how laminar organization may provide a foundational principle of cortical function.

Effective Connectivity of Thalamocortical Interactions Following d-Amphetamine, LSD, and MDMA Administration

BACKGROUND

While the exploration of serotonergic psychedelics as psychiatric medicines deepens, so does the pressure to better understand how these compounds act on the brain.

METHODS

We used a double-blind, placebo-controlled, crossover design and administered lysergic acid diethylamide (LSD), 3,4-methylenedioxymethamphetamine (MDMA), and d-amphetamine in 25 healthy participants. By using spectral dynamic causal modeling, we mapped substance-induced changes in effective connectivity between the thalamus and different cortex types (unimodal vs. transmodal) derived from a previous study with resting-state functional magnetic resonance imaging data. Due to the distinct pharmacological modes of action of the 3 substances, we were able to investigate specific effects mainly driven by different neurotransmitter systems on thalamocortical and corticothalamic interactions.

RESULTS

Compared with placebo, all 3 substances increased the effective connectivity from the thalamus to specific unimodal cortices, whereas the influence of these cortices on the thalamus was reduced. These results indicate increased bottom-up and decreased top-down information flow between the thalamus and some unimodal cortices. However, for the amphetamines, we found the opposite effects when examining the effective connectivity with transmodal cortices, including parts of the salience network. Intriguingly, LSD increased the effective connectivity from the thalamus to both unimodal and transmodal cortices, indicating a breach in the hierarchical organization of ongoing brain activity.

CONCLUSIONS

The results advance our knowledge about the action of psychedelics on the brain and refine current models aiming to explain the underlying neurobiological processes.

Atypical connectome topography and signal flow in temporal lobe epilepsy

Temporal lobe epilepsy (TLE) is one of the most common pharmaco-resistant epilepsies in adults. While hippocampal pathology is the hallmark of this condition, emerging evidence indicates that brain alterations extend beyond the mesiotemporal epicenter and affect macroscale brain function and cognition. We studied macroscale functional reorganization in TLE, explored structural substrates, and examined cognitive associations. We investigated a multisite cohort of 95 patients with pharmaco-resistant TLE and 95 healthy controls using state-of-the-art multimodal 3T magnetic resonance imaging (MRI). We quantified macroscale functional topographic organization using connectome dimensionality reduction techniques and estimated directional functional flow using generative models of effective connectivity. We observed atypical functional topographies in patients with TLE relative to controls, manifesting as reduced functional differentiation between sensory/motor networks and transmodal systems such as the default mode network, with peak alterations in bilateral temporal and ventromedial prefrontal cortices. TLE-related topographic changes were consistent in all three included sites and reflected reductions in hierarchical flow patterns between cortical systems. Integration of parallel multimodal MRI data indicated that these findings were independent of TLE-related cortical grey matter atrophy, but mediated by microstructural alterations in the superficial white matter immediately beneath the cortex. The magnitude of functional perturbations was robustly associated with behavioral markers of memory function. Overall, this work provides converging evidence for macroscale functional imbalances, contributing microstructural alterations, and their associations with cognitive dysfunction in TLE.

The effect of lysergic acid diethylamide (LSD) on whole-brain functional and effective connectivity

Psychedelics have emerged as promising candidate treatments for various psychiatric conditions, and given their clinical potential, there is a need to identify biomarkers that underlie their effects. Here, we investigate the neural mechanisms of lysergic acid diethylamide (LSD) using regression dynamic causal modelling (rDCM), a novel technique that assesses whole-brain effective connectivity (EC) during resting-state functional magnetic resonance imaging (fMRI). We modelled data from two randomised, placebo-controlled, double-blind, cross-over trials, in which 45 participants were administered 100 μg LSD and placebo in two resting-state fMRI sessions. We compared EC against whole-brain functional connectivity (FC) using classical statistics and machine learning methods. Multivariate analyses of EC parameters revealed predominantly stronger interregional connectivity and reduced self-inhibition under LSD compared to placebo, with the notable exception of weakened interregional connectivity and increased self-inhibition in occipital brain regions as well as subcortical regions. Together, these findings suggests that LSD perturbs the Excitation/Inhibition balance of the brain. Notably, whole-brain EC did not only provide additional mechanistic insight into the effects of LSD on the Excitation/Inhibition balance of the brain, but EC also correlated with global subjective effects of LSD and discriminated experimental conditions in a machine learning-based analysis with high accuracy (91.11%), highlighting the potential of using whole-brain EC to decode or predict subjective effects of LSD in the future.

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.

Predicting Future Depressive Episodes from Resting-State fMRI with Generative Embedding

After a first episode of major depressive disorder (MDD), there is substantial risk for a long-term remitting-relapsing course. Prevention and early interventions are thus critically important. Various studies have examined the feasibility of detecting at-risk individuals based on out-of-sample predictions about the future occurrence of depression. However, functional magnetic resonance imaging (MRI) has received very little attention for this purpose so far. Here, we explored the utility of generative models (i.e. different dynamic causal models, DCMs) as well as functional connectivity (FC) for predicting future episodes of depression in never-depressed adults, using a large dataset (N=906) of task-free ("resting state") fMRI data from the UK Biobank. Connectivity analyses were conducted using timeseries from pre-computed spatially independent components of different dimensionalities. Over a three year period, 50% of participants showed indications of at least one depressive episode, while the other 50% did not. Using nested cross-validation for training and a held-out test set (80/20 split), we systematically examined the combination of 8 connectivity feature sets and 17 classifiers. We found that a generative embedding procedure based on combining regression DCM (rDCM) with a support vector machine (SVM) enabled the best predictions, both on the training set (0.63 accuracy, 0.66 area under the curve, AUC) and the test set (0.62 accuracy, 0.64 AUC; p<0.001). However, on the test set, rDCM was only slightly superior to predictions based on FC (0.59 accuracy, 0.61 AUC). Interpreting model predictions based on SHAP (SHapley Additive exPlanations) values suggested that the most predictive connections were widely distributed and not confined to specific networks. Overall, our analyses suggest (i) ways of improving future fMRI-based generative embedding approaches for the early detection of individuals at-risk for depression and that (ii) achieving accuracies of clinical utility may require combination of fMRI with other data modalities.

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.

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.

Fully Bayesian estimation of virtual brain parameters with self-tuning Hamiltonian Monte Carlo

Virtual brain models are data-driven patient-specific brain models integrating individual brain imaging data with neural mass modeling in a single computational framework, capable of autonomously generating brain activity and its associated brain imaging signals. Along the example of epilepsy, we develop an efficient and accurate Bayesian methodology estimating the parameters linked to the extent of the epileptogenic zone. State-of-the-art advances in Bayesian inference using Hamiltonian Monte Carlo (HMC) algorithms have remained elusive for large-scale differential-equations based models due to their slow convergence. We propose appropriate priors and a novel reparameterization to facilitate efficient exploration of the posterior distribution in terms of computational time and convergence diagnostics. The methodology is illustrated for in-silico dataset and then, applied to infer the personalized model parameters based on the empirical stereotactic electroencephalography (SEEG) recordings of retrospective patients. This improved methodology may pave the way to render HMC methods sufficiently easy and efficient to use, thus applicable in personalized medicine.

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.

A strategy of model space search for dynamic causal modeling in task fMRI data exploratory analysis

Dynamic causal modeling (DCM) is a tool used for effective connectivity (EC) estimation in neuroimage analysis. But it is a model-driven analysis method, and the structure of the EC network needs to be determined in advance based on a large amount of prior knowledge. This characteristic makes it difficult to apply DCM to the exploratory brain network analysis. The exploratory analysis of DCM can be realized from two perspectives: one is to reduce the computational cost of the model; the other is to reduce the model space. From the perspective of model space reduction, a model space exploration strategy is proposed, including two algorithms. One algorithm, named GreedyEC, starts with reducing EC from full model, and the other, named GreedyROI, start with adding EC from one node model. Then the two algorithms were applied to the task state functional magnetic resonance imaging (fMRI) data of visual object recognition and selected the best DCM model from the perspective of model comparison based on Bayesian model compare method. Results show that combining the results of the two algorithms can further improve the effect of DCM exploratory analysis. For convenience in application, the algorithms were encapsulated into MATLAB function based on SPM to help neuroscience researchers to analyze the brain causal information flow network. The strategy provides a model space exploration tool that may obtain the best model from the perspective of model comparison and lower the threshold of DCM analysis.

Forecasting brain activity based on models of spatiotemporal brain dynamics: A comparison of graph neural network architectures

Comprehending the interplay between spatial and temporal characteristics of neural dynamics can contribute to our understanding of information processing in the human brain. Graph neural networks (GNNs) provide a new possibility to interpret graph-structured signals like those observed in complex brain networks. In our study we compare different spatiotemporal GNN architectures and study their ability to model neural activity distributions obtained in functional MRI (fMRI) studies. We evaluate the performance of the GNN models on a variety of scenarios in MRI studies and also compare it to a VAR model, which is currently often used for directed functional connectivity analysis. We show that by learning localized functional interactions on the anatomical substrate, GNN-based approaches are able to robustly scale to large network studies, even when available data are scarce. By including anatomical connectivity as the physical substrate for information propagation, such GNNs also provide a multimodal perspective on directed connectivity analysis, offering a novel possibility to investigate the spatiotemporal dynamics in brain networks.

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.

Test-retest reliability of regression dynamic causal modeling

Regression dynamic causal modeling (rDCM) is a novel and computationally highly efficient method for inferring effective connectivity at the whole-brain level. While face and construct validity of rDCM have already been demonstrated, here we assessed its test-retest reliability—a test-theoretical property of particular importance for clinical applications—together with group-level consistency of connection-specific estimates and consistency of whole-brain connectivity patterns over sessions. Using the Human Connectome Project dataset for eight different paradigms (tasks and rest) and two different parcellation schemes, we found that rDCM provided highly consistent connectivity estimates at the group level across sessions. Second, while test-retest reliability was limited when averaging over all connections (range of mean intraclass correlation coefficient 0.24–0.42 over tasks), reliability increased with connection strength, with stronger connections showing good to excellent test-retest reliability. Third, whole-brain connectivity patterns by rDCM allowed for identifying individual participants with high (and in some cases perfect) accuracy. Comparing the test-retest reliability of rDCM connectivity estimates with measures of functional connectivity, rDCM performed favorably—particularly when focusing on strong connections. Generally, for all methods and metrics, task-based connectivity estimates showed greater reliability than those from the resting state. Our results underscore the potential of rDCM for human connectomics and clinical applications.

Optogenetic activation of striatal D1R and D2R cells differentially engages downstream connected areas beyond the basal ganglia

The basal ganglia (BG) are a group of subcortical nuclei responsible for motor and executive function. Central to BG function are striatal cells expressing D1 (D1R) and D2 (D2R) dopamine receptors. D1R and D2R cells are considered functional antagonists that facilitate voluntary movements and inhibit competing motor patterns, respectively. However, whether they maintain a uniform function across the striatum and what influence they exert outside the BG is unclear. Here, we address these questions by combining optogenetic activation of D1R and D2R cells in the mouse ventrolateral caudoputamen with fMRI. Striatal D1R/D2R stimulation evokes distinct activity within the BG-thalamocortical network and differentially engages cerebellar and prefrontal regions. Computational modeling of effective connectivity confirms that changes in D1R/D2R output drive functional relationships between these regions. Our results suggest a complex functional organization of striatal D1R/D2R cells and hint toward an interconnected fronto-BG-cerebellar network modulated by striatal D1R and D2R cells.

The dynamic connectome of speech control

Speech production relies on the orchestrated control of multiple brain regions. The specific, directional influences within these networks remain poorly understood. We used regression dynamic causal modelling to infer the whole-brain directed (effective) connectivity from functional magnetic resonance imaging data of 36 healthy individuals during the production of meaningful English sentences and meaningless syllables. We identified that the two dynamic connectomes have distinct architectures that are dependent on the complexity of task production. The speech was regulated by a dynamic neural network, the most influential nodes of which were centred around superior and inferior parietal areas and influenced the whole-brain network activity via long-ranging coupling with primary sensorimotor, prefrontal, temporal and insular regions. By contrast, syllable production was controlled by a more compressed, cost-efficient network structure, involving sensorimotor cortico-subcortical integration via superior parietal and cerebellar network hubs. These data demonstrate the mechanisms by which the neural network reorganizes the connectivity of its influential regions, from supporting the fundamental aspects of simple syllabic vocal motor output to multimodal information processing of speech motor output. This article is part of the theme issue 'Vocal learning in animals and humans'.

Learning Brain Dynamics With Coupled Low-Dimensional Nonlinear Oscillators and Deep Recurrent Networks

Many natural systems, especially biological ones, exhibit complex multivariate nonlinear dynamical behaviors that can be hard to capture by linear autoregressive models. On the other hand, generic nonlinear models such as deep recurrent neural networks often require large amounts of training data, not always available in domains such as brain imaging; also, they often lack interpretability. Domain knowledge about the types of dynamics typically observed in such systems, such as a certain type of dynamical systems models, could complement purely data-driven techniques by providing a good prior. In this work, we consider a class of ordinary differential equation (ODE) models known as van der Pol (VDP) oscil lators and evaluate their ability to capture a low-dimensional representation of neural activity measured by different brain imaging modalities, such as calcium imaging (CaI) and fMRI, in different living organisms: larval zebrafish, rat, and human. We develop a novel and efficient approach to the nontrivial problem of parameters estimation for a network of coupled dynamical systems from multivariate data and demonstrate that the resulting VDP models are both accurate and interpretable, as VDP's coupling matrix reveals anatomically meaningful excitatory and inhibitory interactions across different brain subsystems. VDP outperforms linear autoregressive models (VAR) in terms of both the data fit accuracy and the quality of insight provided by the coupling matrices and often tends to generalize better to unseen data when predicting future brain activity, being comparable to and sometimes better than the recurrent neural networks (LSTMs). Finally, we demonstrate that our (generative) VDP model can also serve as a data-augmentation tool leading to marked improvements in predictive accuracy of recurrent neural networks. Thus, our work contributes to both basic and applied dimensions of neuroimaging: gaining scientific insights and improving brain-based predictive models, an area of potentially high practical importance in clinical diagnosis and neurotechnology.

Aberrant Temporal Variability in Brain Regions during Risk Decision Making in Patients with Bipolar I Disorder: A Dynamic Effective Connectivity Study

Bipolar I disorder (BD-I) is associated with high-risk behaviors, such as suicide attempts and addictive substance abuse. Understanding brain activity exposure to risk decision making provides evidence for the treatment of BD-I patients. This study aimed to investigate the temporal dynamics of brain connectivity underlying risk decision making in patients with BD-I. A total of 101 subjects (48 BD-I patients and 53 age- and gender-matched healthy controls (HCs)) were included in this research. We analyzed the fMRI data acquired during Balloon Analog Risk Task (BART) performance. Voxel-wise dynamic effective connectivity (dEC) was employed to measure the activities in 264 brain regions. The coefficient of variation (CV) was calculated as temporal dynamics of brain connectivity. Finally, we used structural equation modeling (SEM) to determine the relationships of dEC in brain regions with clinical symptoms, behavior performances in patients. Results showed that BD-I patients exhibited increased dynamics in four lobes and exhibited decreased in three frontal regions. Besides, SEM results showed that the impulsive symptoms of patients were affected by the dEC during both resting and task states. Moreover, the dEC of left supramarginal gyrus (SMG) influenced those of left orbital frontal and right cuneus (CUN), as well as the affective symptoms and BART behaviors in patients with BD-I. Our results suggested that the altered temporal dynamics of brain connectivity might contribute to the impulsivity of BD-I during resting and task states. More importantly, the left SMG might be a therapeutic target to reduce the risk behavior in BD-I patients.

TAPAS: An Open-Source Software Package for Translational Neuromodeling and Computational Psychiatry

Psychiatry faces fundamental challenges with regard to mechanistically guided differential diagnosis, as well as prediction of clinical trajectories and treatment response of individual patients. This has motivated the genesis of two closely intertwined fields: (i) Translational Neuromodeling (TN), which develops "computational assays" for inferring patient-specific disease processes from neuroimaging, electrophysiological, and behavioral data; and (ii) Computational Psychiatry (CP), with the goal of incorporating computational assays into clinical decision making in everyday practice. In order to serve as objective and reliable tools for clinical routine, computational assays require end-to-end pipelines from raw data (input) to clinically useful information (output). While these are yet to be established in clinical practice, individual components of this general end-to-end pipeline are being developed and made openly available for community use. In this paper, we present the Translational Algorithms for Psychiatry-Advancing Science (TAPAS) software package, an open-source collection of building blocks for computational assays in psychiatry. Collectively, the tools in TAPAS presently cover several important aspects of the desired end-to-end pipeline, including: (i) tailored experimental designs and optimization of measurement strategy prior to data acquisition, (ii) quality control during data acquisition, and (iii) artifact correction, statistical inference, and clinical application after data acquisition. Here, we review the different tools within TAPAS and illustrate how these may help provide a deeper understanding of neural and cognitive mechanisms of disease, with the ultimate goal of establishing automatized pipelines for predictions about individual patients. We hope that the openly available tools in TAPAS will contribute to the further development of TN/CP and facilitate the translation of advances in computational neuroscience into clinically relevant computational assays.

Emergence of canonical functional networks from the structural connectome

How do functional brain networks emerge from the underlying wiring of the brain? We examine how resting-state functional activation patterns emerge from the underlying connectivity and length of white matter fibers that constitute its “structural connectome”. By introducing realistic signal transmission delays along fiber projections, we obtain a complex-valued graph Laplacian matrix that depends on two parameters: coupling strength and oscillation frequency. This complex Laplacian admits a complex-valued eigen-basis in the frequency domain that is highly tunable and capable of reproducing the spatial patterns of canonical functional networks without requiring any detailed neural activity modeling. Specific canonical functional networks can be predicted using linear superposition of small subsets of complex eigenmodes. Using a novel parameter inference procedure we show that the complex Laplacian outperforms the real-valued Laplacian in predicting functional networks. The complex Laplacian eigenmodes therefore constitute a tunable yet parsimonious substrate on which a rich repertoire of realistic functional patterns can emerge. Although brain activity is governed by highly complex nonlinear processes and dense connections, our work suggests that simple extensions of linear models to the complex domain effectively approximate rich macroscopic spatial patterns observable on BOLD fMRI.

GABAergic cortical network physiology in frontotemporal lobar degeneration

The clinical syndromes caused by frontotemporal lobar degeneration are heterogeneous, including the behavioural variant frontotemporal dementia (bvFTD) and progressive supranuclear palsy. Although pathologically distinct, they share many behavioural, cognitive and physiological features, which may in part arise from common deficits of major neurotransmitters such as γ-aminobutyric acid (GABA). Here, we quantify the GABAergic impairment and its restoration with dynamic causal modelling of a double-blind placebo-controlled crossover pharmaco-magnetoencephalography study. We analysed 17 patients with bvFTD, 15 patients with progressive supranuclear palsy, and 20 healthy age- and gender-matched controls. In addition to neuropsychological assessment and structural MRI, participants undertook two magnetoencephalography sessions using a roving auditory oddball paradigm: once on placebo and once on 10 mg of the oral GABA reuptake inhibitor tiagabine. A subgroup underwent ultrahigh-field magnetic resonance spectroscopy measurement of GABA concentration, which was reduced among patients. We identified deficits in frontotemporal processing using conductance-based biophysical models of local and global neuronal networks. The clinical relevance of this physiological deficit is indicated by the correlation between top-down connectivity from frontal to temporal cortex and clinical measures of cognitive and behavioural change. A critical validation of the biophysical modelling approach was evidence from parametric empirical Bayes analysis that GABA levels in patients, measured by spectroscopy, were related to posterior estimates of patients’ GABAergic synaptic connectivity. Further evidence for the role of GABA in frontotemporal lobar degeneration came from confirmation that the effects of tiagabine on local circuits depended not only on participant group, but also on individual baseline GABA levels. Specifically, the phasic inhibition of deep cortico-cortical pyramidal neurons following tiagabine, but not placebo, was a function of GABA concentration. The study provides proof-of-concept for the potential of dynamic causal modelling to elucidate mechanisms of human neurodegenerative disease, and explains the variation in response to candidate therapies among patients. The laminar- and neurotransmitter-specific features of the modelling framework, can be used to study other treatment approaches and disorders. In the context of frontotemporal lobar degeneration, we suggest that neurophysiological restoration in selected patients, by targeting neurotransmitter deficits, could be used to bridge between clinical and preclinical models of disease, and inform the personalized selection of drugs and stratification of patients for future clinical trials.

Dynamic causal modeling of eye gaze processing in schizophrenia

BACKGROUND

Abnormal eye gaze perception is related to symptoms and social functioning in schizophrenia. However, little is known about the brain network mechanisms underlying these abnormalities. Here, we employed dynamic causal modeling (DCM) of fMRI data to discover aberrant effective connectivity within networks associated with eye gaze processing in schizophrenia.

METHODS

Twenty-seven patients (schizophrenia/schizoaffective disorder, SZ) and 22 healthy controls (HC) completed an eye gaze processing task during fMRI. Participants viewed faces with different gaze angles and performed explicit gaze discrimination (Gaze: "Looking at you?" yes/no) or implicit gaze processing (Gender: "male or female?"). Four brain regions, the secondary visual cortex (Vis), posterior superior temporal sulcus (pSTS), inferior parietal lobule (IPL), and posterior medial frontal cortex (pMFC) were identified as nodes for subsequent DCM analysis.

RESULTS

SZ and HC showed similar generative model structure, but SZ showed altered connectivity for specific self-connections, inter-regional connections during all gaze processing (reduced excitatory bottom-up and enhanced inhibitory top-down connections), and modulation by explicit gaze discrimination (increased frontal inhibition of visual cortex). Altered effective connectivity was significantly associated with poorer social cognition and functioning.

CONCLUSIONS

General gaze processing in SZ is associated with distributed cortical dysfunctions and bidirectional connectivity between regions, while explicit gaze discrimination involves predominantly top-down abnormalities in the visual system. These results suggest plausible neural mechanisms underpinning gaze processing deficits and may serve as bio-markers for intervention.

Regression dynamic causal modeling for resting-state fMRI

“Resting‐state” functional magnetic resonance imaging (rs‐fMRI) is widely used to study brain connectivity. So far, researchers have been restricted to measures of functional connectivity that are computationally efficient but undirected, or to effective connectivity estimates that are directed but limited to small networks. Here, we show that a method recently developed for task‐fMRI—regression dynamic causal modeling (rDCM)—extends to rs‐fMRI and offers both directional estimates and scalability to whole‐brain networks. First, simulations demonstrate that rDCM faithfully recovers parameter values over a wide range of signal‐to‐noise ratios and repetition times. Second, we test construct validity of rDCM in relation to an established model of effective connectivity, spectral DCM. Using rs‐fMRI data from nearly 200 healthy participants, rDCM produces biologically plausible results consistent with estimates by spectral DCM. Importantly, rDCM is computationally highly efficient, reconstructing whole‐brain networks (>200 areas) within minutes on standard hardware. This opens promising new avenues for connectomics.

Advances in the computational understanding of mental illness

Computational psychiatry is a rapidly growing field attempting to translate advances in computational neuroscience and machine learning into improved outcomes for patients suffering from mental illness. It encompasses both data-driven and theory-driven efforts. Here, recent advances in theory-driven work are reviewed. We argue that the brain is a computational organ. As such, an understanding of the illnesses arising from it will require a computational framework. The review divides work up into three theoretical approaches that have deep mathematical connections: dynamical systems, Bayesian inference and reinforcement learning. We discuss both general and specific challenges for the field, and suggest ways forward.

Dynamic Causal Modeling for fMRI With Wilson-Cowan-Based Neuronal Equations

Dynamic causal modeling (DCM) is an analysis technique that has been successfully used to infer about directed connectivity between brain regions based on imaging data such as functional magnetic resonance imaging (fMRI). Most variants of DCM for fMRI rely on a simple bilinear differential equation for neural activation, making it difficult to interpret the results in terms of local neural dynamics. In this work, we introduce a modification to DCM for fMRI by replacing the bilinear equation with a non-linear Wilson-Cowan based equation and use Bayesian Model Comparison (BMC) to show that this modification improves the model evidences. Improved model evidence of the non-linear model is shown for our empirical data (imitation of facial expressions) and validated by synthetic data as well as an empirical test dataset (attention to visual motion) used in previous foundational papers. For our empirical data, we conduct the analysis for a group of 42 healthy participants who performed an imitation task, activating regions putatively containing the human mirror neuron system (MNS). In this regard, we build 540 models as one family for comparing the standard bilinear with the modified Wilson-Cowan models on the family-level. Using this modification, we can interpret the sigmoid transfer function as an averaged f-I curve of many neurons in a single region with a sigmoidal format. In this way, we can make a direct inference from the macroscopic model to detailed microscopic models. The new DCM variant shows superior model evidence on all tested data sets.

T-Hex: Tilted hexagonal grids for rapid 3D imaging

PURPOSE

The purpose of this work is to devise and demonstrate an encoding strategy for 3D MRI that reconciles high speed with flexible segmentation, uniform k-space density, and benign T 2 ∗ effects.

METHODS

Fast sampling of a 3D k-space is typically accomplished by 2D readouts per shot using EPI trains or spiral readouts. Tilted hexagonal (T-Hex) sampling is a way of acquiring more k-space volume per excitation while maintaining uniform sampling density and a smooth T 2 ∗ filter. The k-space volume covered per shot is controlled by the tilting angle. Image reconstruction is performed with a 3D extension of the iterative SENSE approach, incorporating actual field dynamics and static off-resonance. T-Hex imaging is compared with established 3D schemes in terms of speed and noise performance.

RESULTS

Tilted hexagonal acquisition is found to achieve greater imaging speed than known alternatives, particularly in combination with spiral trajectories. The interplay of the proposed 3D trajectories, array detection, and off-resonance is successfully addressed by iterative inversion of the full signal model. Enhanced coverage per shot is of greatest utility for high speed in an intermediate resolution regime of 1 to 4 mm. T-Hex EPI combines the benefits of extended coverage per shot with increased robustness against off-resonance effects.

CONCLUSION

Sampling of tilted hexagonal grids is a feasible means of gaining 3D imaging speed with near-optimal SNR efficiency and benign depiction properties. It is a particularly promising technique for time-resolved applications such as fMRI.

Whole-brain estimates of directed connectivity for human connectomics

Connectomics is essential for understanding large-scale brain networks but requires that individual connection estimates are neurobiologically interpretable. In particular, a principle of brain organization is that reciprocal connections between cortical areas are functionally asymmetric. This is a challenge for fMRI-based connectomics in humans where only undirected functional connectivity estimates are routinely available. By contrast, whole-brain estimates of effective (directed) connectivity are computationally challenging, and emerging methods require empirical validation. Here, using a motor task at 7T, we demonstrate that a novel generative model can infer known connectivity features in a whole-brain network (>200 regions, >40,000 connections) highly efficiently. Furthermore, graph-theoretical analyses of directed connectivity estimates identify functional roles of motor areas more accurately than undirected functional connectivity estimates. These results, which can be achieved in an entirely unsupervised manner, demonstrate the feasibility of inferring directed connections in whole-brain networks and open new avenues for human connectomics.

Asymmetric high-order anatomical brain connectivity sculpts effective connectivity

Bridging the gap between symmetric, direct white matter brain connectivity and neural dynamics that are often asymmetric and polysynaptic may offer insights into brain architecture, but this remains an unresolved challenge in neuroscience. Here, we used the graph Laplacian matrix to simulate symmetric and asymmetric high-order diffusion processes akin to particles spreading through white matter pathways. The simulated indirect structural connectivity outperformed direct as well as absent anatomical information in sculpting effective connectivity, a measure of causal and directed brain dynamics. Crucially, an asymmetric diffusion process determined by the sensitivity of the network nodes to their afferents best predicted effective connectivity. The outcome is consistent with brain regions adapting to maintain their sensitivity to inputs within a dynamic range. Asymmetric network communication models offer a promising perspective for understanding the relationship between structural and functional brain connectomes, both in normalcy and neuropsychiatric conditions.

Scalable surrogate deconvolution for identification of partially-observable systems and brain modeling

OBJECTIVE

For many biophysical systems, direct measurement of all state-variables, in - vivo is not feasible. Thus, a key challenge in biological modeling and signal processing is to reconstruct the activity and structure of interesting biological systems from indirect measurements. These measurements are often generated by approximately linear time-invariant dynamical interactions with the hidden system and may therefore be described as a convolution of hidden state-variables with an unknown kernel.

APPROACH

In the current work, we present an approach termed surrogate deconvolution, to directly identify such coupled systems (i.e. parameterize models). Surrogate deconvolution reframes certain non linear partially-observable identification problems, which are common in neuroscience/biology, as analytical objectives that are compatible with almost any user-chosen optimization procedure.

MAIN RESULTS

We show that the proposed technique is highly scalable, low in computational complexity, and performs competitively with the current gold-standard in partially-observable system estimation: the joint Kalman Filters (Unscented and Extended). We show the benefits of surrogate deconvolution for model identification when applied to simulations of the Local Field Potential and blood oxygen level dependent (BOLD) signal. Lastly, we demonstrate the empirical stability of Hemodynamic Response Function (HRF) kernel estimates for Mesoscale Individualized NeuroDynamic (MINDy) models of individual human brains. The recovered HRF parameters demonstrate reliable individual variation as well as a stereotyped spatial distribution, on average.

SIGNIFICANCE

These results demonstrate that surrogate deconvolution promises to enhance brain-modeling approaches by simultaneously and rapidly fitting large-scale models of brain networks and the physiological processes which generate neuroscientific measurements (e.g. hemodynamics for BOLD fMRI).

Spectral graph theory of brain oscillations

The relationship between the brain's structural wiring and the functional patterns of neural activity is of fundamental interest in computational neuroscience. We examine a hierarchical, linear graph spectral model of brain activity at mesoscopic and macroscopic scales. The model formulation yields an elegant closed-form solution for the structure-function problem, specified by the graph spectrum of the structural connectome's Laplacian, with simple, universal rules of dynamics specified by a minimal set of global parameters. The resulting parsimonious and analytical solution stands in contrast to complex numerical simulations of high dimensional coupled nonlinear neural field models. This spectral graph model accurately predicts spatial and spectral features of neural oscillatory activity across the brain and was successful in simultaneously reproducing empirically observed spatial and spectral patterns of alpha-band (8-12 Hz) and beta-band (15-30 Hz) activity estimated from source localized magnetoencephalography (MEG). This spectral graph model demonstrates that certain brain oscillations are emergent properties of the graph structure of the structural connectome and provides important insights towards understanding the fundamental relationship between network topology and macroscopic whole-brain dynamics. .

The Bayesian Virtual Epileptic Patient: A probabilistic framework designed to infer the spatial map of epileptogenicity in a personalized large-scale brain model of epilepsy spread

Despite the importance and frequent use of Bayesian frameworks in brain network modeling for parameter inference and model prediction, the advanced sampling algorithms implemented in probabilistic programming languages to overcome the inference difficulties have received relatively little attention in this context. In this technical note, we propose a probabilistic framework, namely the Bayesian Virtual Epileptic Patient (BVEP), which relies on the fusion of structural data of individuals to infer the spatial map of epileptogenicity in a personalized large-scale brain model of epilepsy spread. To invert the individualized whole-brain model employed in this study, we use the recently developed algorithms known as No-U-Turn Sampler (NUTS) as well as Automatic Differentiation Variational Inference (ADVI). Our results indicate that NUTS and ADVI accurately estimate the degree of epileptogenicity of brain regions, therefore, the hypothetical brain areas responsible for the seizure initiation and propagation, while the convergence diagnostics and posterior behavior analysis validate the reliability of the estimations. Moreover, we illustrate the efficiency of the transformed non-centered parameters in comparison to centered form of parameterization. The Bayesian framework used in this work proposes an appropriate patient-specific strategy for estimating the epileptogenicity of the brain regions to improve outcome after epilepsy surgery.

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

Neuroimaging techniques are now widely used to study human cognition. The functional associations between brain areas have become a standard proxy to describe how cognitive processes are distributed across the brain network. Among the many analysis tools available, dynamic models of brain activity have been developed to overcome the limitations of original connectivity measures such as functional connectivity. This goes in line with the many efforts devoted to the assessment of directional interactions between brain areas from the observed neuroimaging activity. This opinion article provides an overview of our model-based whole-brain effective connectivity to analyze fMRI data, while discussing the pros and cons of our approach with respect to other established approaches. Our framework relies on the multivariate Ornstein-Uhlenbeck (MOU) process and is thus referred to as MOU-EC. Once tuned, the model provides a directed connectivity estimate that reflects the dynamical state of BOLD activity, which can be used to explore cognition. We illustrate this approach using two applications on task-evoked fMRI data. First, as a connectivity measure, MOU-EC can be used to extract biomarkers for task-specific brain coordination, understood as the patterns of areas exchanging information. The multivariate nature of connectivity measures raises several challenges for whole-brain analysis, for which machine-learning tools present some advantages over statistical testing. Second, we show how to interpret changes in MOU-EC connections in a collective and model-based manner, bridging with network analysis. Our framework provides a comprehensive set of tools that open exciting perspectives to study distributed cognition, as well as neuropathologies.

Multi-Level Clustering of Dynamic Directional Brain Network Patterns and Their Behavioral Relevance

Dynamic functional connectivity (DFC) obtained from resting state functional magnetic resonance imaging (fMRI) data has been shown to provide novel insights into brain function which may be obscured by static functional connectivity (SFC). Further, DFC, and by implication how different brain regions may engage or disengage with each other over time, has been shown to be behaviorally relevant and more predictive than SFC of behavioral performance and/or diagnostic status. DFC is not a directional entity and may capture neural synchronization. However, directional interactions between different brain regions is another putative mechanism by which neural populations communicate. Accordingly, static effective connectivity (SEC) has been explored as a means of characterizing such directional interactions. But investigation of its dynamic counterpart, i.e., dynamic effective connectivity (DEC), is still in its infancy. Of particular note are methodological insufficiencies in identifying DEC configurations that are reproducible across time and subjects as well as a lack of understanding of the behavioral relevance of DEC obtained from resting state fMRI. In order to address these issues, we employed a dynamic multivariate autoregressive (MVAR) model to estimate DEC. The method was first validated using simulations and then applied to resting state fMRI data obtained in-house (N = 21), wherein we performed dynamic clustering of DEC matrices across multiple levels [using adaptive evolutionary clustering (AEC)] – spatial location, time, and subjects. We observed a small number of directional brain network configurations alternating between each other over time in a quasi-stable manner akin to brain microstates. The dominant and consistent DEC network patterns involved several regions including inferior and mid temporal cortex, motor and parietal cortex, occipital cortex, as well as part of frontal cortex. The functional relevance of these DEC states were determined using meta-analyses and pertained mainly to memory and emotion, but also involved execution and language. Finally, a larger cohort of resting-state fMRI and behavioral data from the Human Connectome Project (HCP) (N = 232, Q1–Q3 release) was used to demonstrate that metrics derived from DEC can explain larger variance in 70 behaviors across different domains (alertness, cognition, emotion, and personality traits) compared to SEC in healthy individuals.

Predicting individual clinical trajectories of depression with generative embedding

•

Patients with major depressive disorder (MDD) show variable clinical trajectories.

•

Generative embedding (GE) is used to predict clinical trajectories in MDD patients.

•

GE classifies patients with chronic depression vs. fast remission with 79% accuracy.

•

GE provides mechanistic interpretability and outperforms conventional measures.

•

Proof-of-concept that illustrates the potential of GE for clinical prediction.

Patients with major depressive disorder (MDD) show heterogeneous treatment response and highly variable clinical trajectories: while some patients experience swift recovery, others show relapsing-remitting or chronic courses. Predicting individual clinical trajectories at an early stage is a key challenge for psychiatry and might facilitate individually tailored interventions. So far, however, reliable predictors at the single-patient level are absent. Here, we evaluated the utility of a machine learning strategy – generative embedding (GE) – which combines interpretable generative models with discriminative classifiers. Specifically, we used functional magnetic resonance imaging (fMRI) data of emotional face perception in 85 MDD patients from the NEtherlands Study of Depression and Anxiety (NESDA) who had been followed up over two years and classified into three subgroups with distinct clinical trajectories. Combining a generative model of effective (directed) connectivity with support vector machines (SVMs), we could predict whether a given patient would experience chronic depression vs. fast remission with a balanced accuracy of 79%. Gradual improvement vs. fast remission could still be predicted above-chance, but less convincingly, with a balanced accuracy of 61%. Generative embedding outperformed classification based on conventional (descriptive) features, such as functional connectivity or local activation estimates, which were obtained from the same data and did not allow for above-chance classification accuracy. Furthermore, predictive performance of GE could be assigned to a specific network property: the trial-by-trial modulation of connections by emotional content. Given the limited sample size of our study, the present results are preliminary but may serve as proof-of-concept, illustrating the potential of GE for obtaining clinical predictions that are interpretable in terms of network mechanisms. Our findings suggest that abnormal dynamic changes of connections involved in emotional face processing might be associated with higher risk of developing a less favorable clinical course.

Intrinsic connectome organization across temporal scales: New insights from cross-modal approaches

The discovery of a stable, whole-brain functional connectivity organization that is largely independent of external events has drastically extended our view of human brain function. However, this discovery has been primarily based on functional magnetic resonance imaging (fMRI). The role of this whole-brain organization in fast oscillation-based connectivity as measured, for example, by electroencephalography (EEG) and magnetoencephalography (MEG) is only beginning to emerge. Here, we review studies of intrinsic connectivity and its whole-brain organization in EEG, MEG, and intracranial electrophysiology with a particular focus on direct comparisons to connectome studies in fMRI. Synthesizing this literature, we conclude that irrespective of temporal scale over four orders of magnitude, intrinsic neurophysiological connectivity shows spatial similarity to the connectivity organization commonly observed in fMRI. A shared structural connectivity basis and cross-frequency coupling are possible mechanisms contributing to this similarity. Acknowledging that a stable whole-brain organization governs long-range coupling across all timescales of neural processing motivates researchers to take "baseline" intrinsic connectivity into account when investigating brain-behavior associations, and further encourages more widespread exploration of functional connectomics approaches beyond fMRI by using EEG and MEG modalities.

A Computational Theory of Mindfulness Based Cognitive Therapy from the "Bayesian Brain" Perspective

Mindfulness Based Cognitive Therapy (MBCT) was developed to combine methods from cognitive behavioral therapy and meditative techniques, with the specific goal of preventing relapse in recurrent depression. While supported by empirical evidence from multiple clinical trials, the cognitive mechanisms behind the effectiveness of MBCT are not well understood in computational (information processing) or biological terms. This article introduces a testable theory about the computational mechanisms behind MBCT that is grounded in "Bayesian brain" concepts of perception from cognitive neuroscience, such as predictive coding. These concepts regard the brain as embodying a model of its environment (including the external world and the body) which predicts future sensory inputs and is updated by prediction errors, depending on how precise these error signals are. This article offers a concrete proposal how core concepts of MBCT-(i) the being mode (accepting whatever sensations arise, without judging or changing them), (ii) decentering (experiencing thoughts and percepts simply as events in the mind that arise and pass), and (iii) cognitive reactivity (changes in mood reactivate negative beliefs)-could be understood in terms of perceptual and metacognitive processes that draw on specific computational mechanisms of the "Bayesian brain." Importantly, the proposed theory can be tested experimentally, using a combination of behavioral paradigms, computational modelling, and neuroimaging. The novel theoretical perspective on MBCT described in this paper may offer opportunities for finessing the conceptual and practical aspects of MBCT.

Modeling functional resting-state brain networks through neural message passing on the human connectome

In this work, we propose a natural model for information flow in the brain through a neural message-passing dynamics on a structural network of macroscopic regions, such as the human connectome (HC). In our model, each brain region is assumed to have a binary behavior (active or not), the strengths of interactions among them are encoded in the anatomical connectivity matrix defined by the HC, and the dynamics of the system is defined by the Belief Propagation (BP) algorithm, working near the critical point of the network. We show that in the absence of direct external stimuli the BP algorithm converges to a spatial map of activations that is similar to the Default Mode Network (DMN) of the brain, which has been defined from the analysis of functional MRI data. Moreover, we use Susceptibility Propagation (SP) to compute the matrix of long-range correlations between the different regions and show that the modules defined by a clustering of this matrix resemble several Resting State Networks (RSN) determined experimentally. Both results suggest that the functional DMN and RSNs can be seen as simple consequences of the anatomical structure of the brain and a neural message-passing dynamics between macroscopic regions. With the new model, we explore predictions on how functional maps change when the anatomical brain network suffers structural alterations, like in Alzheimer's disease and in lesions of the Corpus Callosum. The implications and novel interpretations suggested by the model, as well as the role of criticality, are discussed.

Advancing functional connectivity research from association to causation

Cognition and behavior emerge from brain network interactions, such that investigating causal interactions should be central to the study of brain function. Approaches that characterize statistical associations among neural time series-functional connectivity (FC) methods-are likely a good starting point for estimating brain network interactions. Yet only a subset of FC methods ('effective connectivity') is explicitly designed to infer causal interactions from statistical associations. Here we incorporate best practices from diverse areas of FC research to illustrate how FC methods can be refined to improve inferences about neural mechanisms, with properties of causal neural interactions as a common ontology to facilitate cumulative progress across FC approaches. We further demonstrate how the most common FC measures (correlation and coherence) reduce the set of likely causal models, facilitating causal inferences despite major limitations. Alternative FC measures are suggested to immediately start improving causal inferences beyond these common FC measures.

fMRI in Non-human Primate: A Review on Factors That Can Affect Interpretation and Dynamic Causal Modeling Application

Dynamic causal modeling (DCM)-a framework for inferring hidden neuronal states from brain activity measurements (e. g., fMRI) and their context-dependent modulation-was developed for human neuroimaging, and has not been optimized for non-human primate (NHP) studies, which are usually done under anesthesia. Animal neuroimaging studies offer the potential to improve effective connectivity modeling using DCM through combining functional imaging with invasive procedures such as in vivo optogenetic or electrical stimulation. Employing a Bayesian approach, model parameters are estimated based on prior knowledge of conditions that might be related to neural and BOLD dynamics (e.g., requires empirical knowledge about the range of plausible parameter values). As such, we address the following questions in this review: What factors need to be considered when applying DCM to NHP data? What differences in functional networks, cerebrovascular architecture and physiology exist between human and NHPs that are relevant for DCM application? How do anesthetics affect vascular physiology, BOLD contrast, and neural dynamics-particularly, effective communication within, and between networks? Considering the factors that are relevant for DCM application to NHP neuroimaging, we propose a strategy for modeling effective connectivity under anesthesia using an integrated physiologic-stochastic DCM (IPS-DCM).

A new computational approach to estimate whole-brain effective connectivity from functional and structural MRI, applied to language development

Recently introduced effective connectivity methods allow for the in-vivo investigation of large-scale functional interactions between brain regions. However, dynamic causal modeling, the most widely used technique to date, typically captures only a few predefined regions of interest. In this study, we present an alternative computational approach to infer effective connectivity within the entire connectome and show its performance on a developmental cohort with emerging language capacities. The novel approach provides new opportunities to quantify effective connectivity changes in the human brain.

Pathophysiological and cognitive mechanisms of fatigue in multiple sclerosis

Fatigue is one of the most common symptoms in multiple sclerosis (MS), with a major impact on patients' quality of life. Currently, treatment proceeds by trial and error with limited success, probably due to the presence of multiple different underlying mechanisms. Recent neuroscientific advances offer the potential to develop tools for differentiating these mechanisms in individual patients and ultimately provide a principled basis for treatment selection. However, development of these tools for differential diagnosis will require guidance by pathophysiological and cognitive theories that propose mechanisms which can be assessed in individual patients. This article provides an overview of contemporary pathophysiological theories of fatigue in MS and discusses how the mechanisms they propose may become measurable with emerging technologies and thus lay a foundation for future personalised treatments.

A Functional Data Method for Causal Dynamic Network Modeling of Task-Related fMRI

Functional MRI (fMRI) is a popular approach to investigate brain connections and activations when human subjects perform tasks. Because fMRI measures the indirect and convoluted signals of brain activities at a lower temporal resolution, complex differential equation modeling methods (e.g., Dynamic Causal Modeling) are usually employed to infer the neuronal processes and to fit the resulting fMRI signals. However, this modeling strategy is computationally expensive and remains to be mostly a confirmatory or hypothesis-driven approach. One major statistical challenge here is to infer, in a data-driven fashion, the underlying differential equation models from fMRI data. In this paper, we propose a causal dynamic network (CDN) method to estimate brain activations and connections simultaneously. Our method links the observed fMRI data with the latent neuronal states modeled by an ordinary differential equation (ODE) model. Using the basis function expansion approach in functional data analysis, we develop an optimization-based criterion that combines data-fitting errors and ODE fitting errors. We also develop and implement a block coordinate-descent algorithm to compute the ODE parameters efficiently. We illustrate the numerical advantages of our approach using data from realistic simulations and two task-related fMRI experiments. Compared with various effective connectivity methods, our method achieves higher estimation accuracy while improving the computational speed by from tens to thousands of times. Though our method is developed for task-related fMRI, we also demonstrate the potential applicability of our method (with a simple modification) to resting-state fMRI, by analyzing both simulated and real data from medium-sized networks.

Disentangling causal webs in the brain using functional magnetic resonance imaging: A review of current approaches

In the past two decades, functional Magnetic Resonance Imaging (fMRI) has been used to relate neuronal network activity to cognitive processing and behavior. Recently this approach has been augmented by algorithms that allow us to infer causal links between component populations of neuronal networks. Multiple inference procedures have been proposed to approach this research question but so far, each method has limitations when it comes to establishing whole-brain connectivity patterns. In this paper, we discuss eight ways to infer causality in fMRI research: Bayesian Nets, Dynamical Causal Modelling, Granger Causality, Likelihood Ratios, Linear Non-Gaussian Acyclic Models, Patel's Tau, Structural Equation Modelling, and Transfer Entropy. We finish with formulating some recommendations for the future directions in this area.