Multivariate dynamical systems models for estimating causal interactions in fMRI

The subcortical brain regions influence the cortical areas during resting-state: an fMRI study

Introduction

Numerous modes or patterns of neural activity can be seen in the brain of individuals during the resting state. However, those functions do not persist long, and they are continuously altering in the brain. We have hypothesized that the brain activations during the resting state should themselves be responsible for this alteration of the activities.

Methods

Using the resting-state fMRI data of 63 healthy young individuals, we estimated the causality effects of each resting-state activation map on all other networks. The resting-state networks were identified, their causality effects on the other components were extracted, the networks with the top 20% of the causality were chosen, and the networks which were under the influence of those causal networks were also identified.

Results

Our results showed that the influence of each activation component over other components is different. The brain areas which showed the highest causality coefficients were subcortical regions, such as the brain stem, thalamus, and amygdala. On the other hand, nearly all the areas which were mostly under the causal effects were cortical regions.

Discussion

In summary, our results suggest that subcortical brain areas exert a higher influence on cortical regions during the resting state, which could help in a better understanding the dynamic nature of brain functions.

Robust and replicable functional brain signatures of 22q11.2 deletion syndrome and associated psychosis: a deep neural network-based multi-cohort study

A major genetic risk factor for psychosis is 22q11.2 deletion (22q11.2DS). However, robust and replicable functional brain signatures of 22q11.2DS and 22q11.2DS-associated psychosis remain elusive due to small sample sizes and a focus on small single-site cohorts. Here, we identify functional brain signatures of 22q11.2DS and 22q11.2DS-associated psychosis, and their links with idiopathic early psychosis, using one of the largest multi-cohort data to date. We obtained multi-cohort clinical phenotypic and task-free fMRI data from 856 participants (101 22q11.2DS, 120 idiopathic early psychosis, 101 idiopathic autism, 123 idiopathic ADHD, and 411 healthy controls) in a case-control design. A novel spatiotemporal deep neural network (stDNN)-based analysis was applied to the multi-cohort data to identify functional brain signatures of 22q11.2DS and 22q11.2DS-associated psychosis. Next, stDNN was used to test the hypothesis that the functional brain signatures of 22q11.2DS-associated psychosis overlap with idiopathic early psychosis but not with autism and ADHD. stDNN-derived brain signatures distinguished 22q11.2DS from controls, and 22q11.2DS-associated psychosis with very high accuracies (86-94%) in the primary cohort and two fully independent cohorts without additional training. Robust distinguishing features of 22q11.2DS-associated psychosis emerged in the anterior insula node of the salience network and the striatum node of the dopaminergic reward pathway. These features also distinguished individuals with idiopathic early psychosis from controls, but not idiopathic autism or ADHD. Our results reveal that individuals with 22q11.2DS exhibit a highly distinct functional brain organization compared to controls. Additionally, the brain signatures of 22q11.2DS-associated psychosis overlap with those of idiopathic early psychosis in the salience network and dopaminergic reward pathway, providing substantial empirical support for the theoretical aberrant salience-based model of psychosis. Collectively, our findings, replicated across multiple independent cohorts, advance the understanding of 22q11.2DS and associated psychosis, underscoring the value of 22q11.2DS as a genetic model for probing the neurobiological underpinnings of psychosis and its progression.

Deep learning models reveal replicable, generalizable, and behaviorally relevant sex differences in human functional brain organization

Sex plays a crucial role in human brain development, aging, and the manifestation of psychiatric and neurological disorders. However, our understanding of sex differences in human functional brain organization and their behavioral consequences has been hindered by inconsistent findings and a lack of replication. Here, we address these challenges using a spatiotemporal deep neural network (stDNN) model to uncover latent functional brain dynamics that distinguish male and female brains. Our stDNN model accurately differentiated male and female brains, demonstrating consistently high cross-validation accuracy (>90%), replicability, and generalizability across multisession data from the same individuals and three independent cohorts (N ~ 1,500 young adults aged 20 to 35). Explainable AI (XAI) analysis revealed that brain features associated with the default mode network, striatum, and limbic network consistently exhibited significant sex differences (effect sizes > 1.5) across sessions and independent cohorts. Furthermore, XAI-derived brain features accurately predicted sex-specific cognitive profiles, a finding that was also independently replicated. Our results demonstrate that sex differences in functional brain dynamics are not only highly replicable and generalizable but also behaviorally relevant, challenging the notion of a continuum in male-female brain organization. Our findings underscore the crucial role of sex as a biological determinant in human brain organization, have significant implications for developing personalized sex-specific biomarkers in psychiatric and neurological disorders, and provide innovative AI-based computational tools for future research.

Causalized convergent cross-mapping and its approximate equivalence with directed information in causality analysis

Convergent cross-mapping (CCM) has attracted increased attention recently due to its capability to detect causality in nonseparable systems under deterministic settings, which may not be covered by the traditional Granger causality. From an information-theoretic perspective, causality is often characterized as the directed information (DI) flowing from one side to the other. As information is essentially nondeterministic, a natural question is: does CCM measure DI flow? Here, we first causalize CCM so that it aligns with the presumption in causality analysis-the future values of one process cannot influence the past of the other, and then establish and validate the approximate equivalence of causalized CCM (cCCM) and DI under Gaussian variables through both theoretical derivations and fMRI-based brain network causality analysis. Our simulation result indicates that, in general, cCCM tends to be more robust than DI in causality detection. The underlying argument is that DI relies heavily on probability estimation, which is sensitive to data size as well as digitization procedures; cCCM, on the other hand, gets around this problem through geometric cross-mapping between the manifolds involved. Overall, our analysis demonstrates that cross-mapping provides an alternative way to evaluate DI and is potentially an effective technique for identifying both linear and nonlinear causal coupling in brain neural networks and other settings, either random or deterministic, or both.

The confound of hemodynamic response function variability in human resting-state functional MRI studies

Functional magnetic resonance imaging (fMRI) is an indirect measure of neural activity with the hemodynamic response function (HRF) coupling it with unmeasured neural activity. The HRF, modulated by several non-neural factors, is variable across brain regions, individuals and populations. Yet, a majority of human resting-state fMRI connectivity studies continue to assume a non-variable HRF. In this article, with supportive prior evidence, we argue that HRF variability cannot be ignored as it substantially confounds within-subject connectivity estimates and between-subjects connectivity group differences. We also discuss its clinical relevance with connectivity impairments confounded by HRF aberrations in several disorders. We present limited data on HRF differences between women and men, which resulted in a 15.4% median error in functional connectivity estimates in a group-level comparison. We also discuss the implications of HRF variability for fMRI studies in the spinal cord. There is a need for more dialogue within the community on the HRF confound, and we hope that our article is a catalyst in the process.

Neuronal dynamics of the default mode network and anterior insular cortex: Intrinsic properties and modulation by salient stimuli

The default mode network (DMN) is critical for self-referential mental processes, and its dysfunction is implicated in many neuropsychiatric disorders. However, the neurophysiological properties and task-based functional organization of the rodent DMN are poorly understood, limiting its translational utility. Here, we combine fiber photometry with functional magnetic resonance imaging (fMRI) and computational modeling to characterize dynamics of putative rat DMN nodes and their interactions with the anterior insular cortex (AI) of the salience network. Our analysis revealed neuronal activity changes in AI and DMN nodes preceding fMRI-derived DMN activations and cyclical transitions between brain network states. Furthermore, we demonstrate that salient oddball stimuli suppress the DMN and enhance AI neuronal activity and that the AI causally inhibits the retrosplenial cortex, a prominent DMN node. These findings elucidate the neurophysiological foundations of the rodent DMN, its spatiotemporal dynamical properties, and modulation by salient stimuli, paving the way for future translational studies.

Dopaminergic medication normalizes aberrant cognitive control circuit signalling in Parkinson's disease

Dopaminergic medication is widely used to alleviate motor symptoms of Parkinson's disease, but these medications also impact cognition with significant variability across patients. It is hypothesized that dopaminergic medication impacts cognition and working memory in Parkinson's disease by modulating frontoparietal-basal ganglia cognitive control circuits, but little is known about the underlying causal signalling mechanisms and their relation to individual differences in response to dopaminergic medication. Here we use a novel state-space computational model with ultra-fast (490 ms resolution) functional MRI to investigate dynamic causal signalling in frontoparietal-basal ganglia circuits associated with working memory in 44 Parkinson's disease patients ON and OFF dopaminergic medication, as well as matched 36 healthy controls. Our analysis revealed aberrant causal signalling in frontoparietal-basal ganglia circuits in Parkinson's disease patients OFF medication. Importantly, aberrant signalling was normalized by dopaminergic medication and a novel quantitative distance measure predicted individual differences in cognitive change associated with medication in Parkinson's disease patients. These findings were specific to causal signalling measures, as no such effects were detected with conventional non-causal connectivity measures. Our analysis also identified a specific frontoparietal causal signalling pathway from right middle frontal gyrus to right posterior parietal cortex that is impaired in Parkinson's disease. Unlike in healthy controls, the strength of causal interactions in this pathway did not increase with working memory load and the strength of load-dependent causal weights was not related to individual differences in working memory task performance in Parkinson's disease patients OFF medication. However, dopaminergic medication in Parkinson's disease patients reinstated the relation with working memory performance. Our findings provide new insights into aberrant causal brain circuit dynamics during working memory and identify mechanisms by which dopaminergic medication normalizes cognitive control circuits.

Developmental maturation of causal signaling hubs in voluntary control of saccades and their functional controllability

The ability to adaptively respond to behaviorally relevant cues in the environment, including voluntary control of automatic but inappropriate responses and deployment of a goal-relevant alternative response, undergoes significant maturation from childhood to adulthood. Importantly, the maturation of voluntary control processes influences the developmental trajectories of several key cognitive domains, including executive function and emotion regulation. Understanding the maturation of voluntary control is therefore of fundamental importance, but little is known about the underlying causal functional circuit mechanisms. Here, we use state-space and control-theoretic modeling to investigate the maturation of causal signaling mechanisms underlying voluntary control over saccades. We demonstrate that directed causal interactions in a canonical saccade network undergo significant maturation between childhood and adulthood. Crucially, we show that the frontal eye field (FEF) is an immature causal signaling hub in children during control over saccades. Using control-theoretic analysis, we then demonstrate that the saccade network is less controllable in children and that greater energy is required to drive FEF dynamics in children compared to adults. Our findings provide novel evidence that strengthening of causal signaling hubs and controllability of FEF are key mechanisms underlying age-related improvements in the ability to plan and execute voluntary control over saccades.

Robust, Generalizable, and Interpretable Artificial Intelligence-Derived Brain Fingerprints of Autism and Social Communication Symptom Severity

BACKGROUND

Autism spectrum disorder (ASD) is among the most pervasive neurodevelopmental disorders, yet the neurobiology of ASD is still poorly understood because inconsistent findings from underpowered individual studies preclude the identification of robust and interpretable neurobiological markers and predictors of clinical symptoms.

METHODS

We leverage multiple brain imaging cohorts and exciting recent advances in explainable artificial intelligence to develop a novel spatiotemporal deep neural network (stDNN) model, which identifies robust and interpretable dynamic brain markers that distinguish ASD from neurotypical control subjects and predict clinical symptom severity.

RESULTS

stDNN achieved consistently high classification accuracies in cross-validation analysis of data from the multisite ABIDE (Autism Brain Imaging Data Exchange) cohort (n = 834). Crucially, stDNN also accurately classified data from independent Stanford (n = 202) and GENDAAR (Gender Exploration of Neurogenetics and Development to Advanced Autism Research) (n = 90) cohorts without additional training. stDNN could not distinguish attention-deficit/hyperactivity disorder from neurotypical control subjects, highlighting the model's specificity. Explainable artificial intelligence revealed that brain features associated with the posterior cingulate cortex and precuneus, dorsolateral and ventrolateral prefrontal cortex, and superior temporal sulcus, which anchor the default mode network, cognitive control, and human voice processing systems, respectively, most clearly distinguished ASD from neurotypical control subjects in the three cohorts. Furthermore, features associated with the posterior cingulate cortex and precuneus nodes of the default mode network emerged as robust predictors of the severity of core social and communication deficits but not restricted/repetitive behaviors in ASD.

CONCLUSIONS

Our findings, replicated across independent cohorts, reveal robust individualized functional brain fingerprints of ASD psychopathology, which could lead to more objective and precise phenotypic characterization and targeted treatments.

Dynamic neural reconfiguration for distinct strategies during competitive social interactions

Information exchange between brain regions is key to understanding information processing for social decision-making, but most analyses ignore its dynamic nature. New insights on this dynamic might help us to uncover the neural correlates of social cognition in the healthy population and also to understand the malfunctioning neural computations underlying dysfunctional social behavior in patients with mental disorders. In this work, we used a multi-round bargaining game to detect switches between distinct bargaining strategies in a cohort of 76 healthy participants. These switches were uncovered by dynamic behavioral modeling using the hidden Markov model. Proposing a novel model of dynamic effective connectivity to estimate the information flow between key brain regions, we found a stronger interaction between the right temporoparietal junction (rTPJ) and the right dorsolateral prefrontal cortex (rDLPFC) for the strategic deception compared with the social heuristic strategies. The level of deception was associated with the information flow from the Brodmann area 10 to the rTPJ, and this association was modulated by the rTPJ-to-rDLPFC information flow. These findings suggest that dynamic bargaining strategy is supported by dynamic reconfiguration of the rDLPFC-and-rTPJ interaction during competitive social interactions.

Neural Algorithms and Circuits for Motor Planning

The brain plans and executes volitional movements. The underlying patterns of neural population activity have been explored in the context of movements of the eyes, limbs, tongue, and head in nonhuman primates and rodents. How do networks of neurons produce the slow neural dynamics that prepare specific movements and the fast dynamics that ultimately initiate these movements? Recent work exploits rapid and calibrated perturbations of neural activity to test specific dynamical systems models that are capable of producing the observed neural activity. These joint experimental and computational studies show that cortical dynamics during motor planning reflect fixed points of neural activity (attractors). Subcortical control signals reshape and move attractors over multiple timescales, causing commitment to specific actions and rapid transitions to movement execution. Experiments in rodents are beginning to reveal how these algorithms are implemented at the level of brain-wide neural circuits.

Neurocircuit dynamics of arbitration between decision-making strategies across obsessive-compulsive and related disorders

•

Obsessive-compulsive and related disorders (OCRD) include OCD and BDD.

•

Neural differences in decision-making arbitration may underlie OCRD symptoms.

•

Resting-state effective connectivity was used to assess arbitration circuitry.

•

Greater left putamen inhibition via left ventrolateral prefrontal cortex in OCRD.

•

Stronger left putamen inhibition was correlated with less severe symptoms.

Obsessions and compulsions are central components of obsessive–compulsive disorder (OCD) and obsessive–compulsive related disorders such as body dysmorphic disorder (BDD). Compulsive behaviours may result from an imbalance of habitual and goal-directed decision-making strategies. The relationship between these symptoms and the neural circuitry underlying habitual and goal-directed decision-making, and the arbitration between these strategies, remains unknown. This study examined resting state effective connectivity between nodes of these systems in two cohorts with obsessions and compulsions, each compared with their own corresponding healthy controls: OCD (n_OCD = 43; n_healthy = 24) and BDD (n_BDD = 21; n_healthy = 16). In individuals with OCD, the left ventrolateral prefrontal cortex, a node of the arbitration system, exhibited more inhibitory causal influence over the left posterolateral putamen, a node of the habitual system, compared with controls. Inhibitory causal influence in this connection showed a trend for a similar pattern in individuals with BDD compared with controls. Those with stronger negative connectivity had lower obsession and compulsion severity in both those with OCD and those with BDD. These relationships were not evident within the habitual or goal-directed circuits, nor were they associated with depressive or anxious symptomatology. These results suggest that abnormalities in the arbitration system may represent a shared neural phenotype across these two related disorders that is specific to obsessive–compulsive symptoms. In addition to nosological implications, these results identify potential targets for novel, circuit-specific treatments.

Neuroimaging approaches to cognition in Parkinson's disease

While direct visualization of Lewy body accumulation within the brain is not yet possible in living Parkinson's disease patients, brain imaging studies offer insights into how the buildup of Lewy body pathology impacts different regions of the brain. Unlike biological biomarkers and purely behavioral research, these brain imaging studies therefore offer a unique opportunity to relate brain localization to cognitive function and dysfunction in living patients. Magnetic resonance imaging studies can reveal physical changes in brain structure as they relate to different cognitive domains and task specific impairments. Functional imaging studies use a combination of task and resting state magnetic resonance imaging, as well as positron emission tomography and single photon emission computed tomography, and can be used to determine changes in blood flow, neuronal activation and neurochemical changes in the brain associated with PD cognition and cognitive impairments. Other unique advantages to brain imaging studies are the ability to monitor changes in brain structure and function longitudinally as patients progress and the ability to study changes in brain function when patients are exposed to different pharmacological manipulations. This is particularly true when assessing the effects of dopaminergic replacement therapy on cognitive function in Parkinson's disease patients. Together, this chapter will describe imaging studies that have helped identify structural and functional brain changes associated with cognition, cognitive impairment, and dementia in Parkinson's disease.

Dynamic causal brain circuits during working memory and their functional controllability

Control processes associated with working memory play a central role in human cognition, but their underlying dynamic brain circuit mechanisms are poorly understood. Here we use system identification, network science, stability analysis, and control theory to probe functional circuit dynamics during working memory task performance. Our results show that dynamic signaling between distributed brain areas encompassing the salience (SN), fronto-parietal (FPN), and default mode networks can distinguish between working memory load and predict performance. Network analysis of directed causal influences suggests the anterior insula node of the SN and dorsolateral prefrontal cortex node of the FPN are causal outflow and inflow hubs, respectively. Network controllability decreases with working memory load and SN nodes show the highest functional controllability. Our findings reveal dissociable roles of the SN and FPN in systems control and provide novel insights into dynamic circuit mechanisms by which cognitive control circuits operate asymmetrically during cognition.

Large-scale network dynamics in neural response to emotionally negative stimuli linked to serotonin 1A binding in major depressive disorder

Serotonergic dysfunction is implicated in major depressive disorder (MDD), but the mechanisms of this relationship remain elusive. Serotonin 1A (5-HT_1A) autoreceptors regulate brain-wide serotonin neuron firing and are positioned to assert large-scale effects on negative emotion. Here we investigated the relationship between raphe 5-HT_1A binding and brain-wide network dynamics of negative emotion. 22 healthy-volunteers (HV) and 27 medication-free participants with MDD underwent positron emission tomography (PET) using [¹¹C]CUMI-101 (CUMI) to quantify 5-HT_1A binding in midbrain raphe nuclei and functional magnetic resonance imaging (fMRI) scanning during emotionally negative picture viewing. Causal connectivity across regions responsive to negative emotion was estimated in the fMRI data using a multivariate dynamical systems model. During negative picture viewing, MDD subjects demonstrated significant hippocampal inhibition of amygdala, basal-ganglia, thalamus, orbital frontal cortex, inferior frontal gyrus and dorsomedial prefrontal cortex (IFG, dmPFC). MDD-related connectivity was not associated with raphe 5-HT_1A binding. However, greater hippocampal inhibition of amygdala, thalamus, IFG and dmPFC correlated with hippocampal 5-HT_1A binding. Correlation between hippocampal 5-HT_1A binding and the hippocampal inhibition network was specific to MDD but not HV. MDD and HV groups also differed with respect to the correlation between raphe and hippocampal 5-HT_1A binding which was more pronounced in HV. These findings suggest that increased hippocampal network inhibition in MDD is linked to hippocampal serotonergic dysfunction which may in turn arise from disrupted linkage in raphe to hippocampus serotonergic circuitry.

FMRI hemodynamic response function (HRF) as a novel marker of brain function: applications for understanding obsessive-compulsive disorder pathology and treatment response

The hemodynamic response function (HRF) represents the transfer function linking neural activity with the functional MRI (fMRI) signal, modeling neurovascular coupling. Since HRF is influenced by non-neural factors, to date it has largely been considered as a confound or has been ignored in many analyses. However, underlying biophysics suggests that the HRF may contain meaningful correlates of neural activity, which might be unavailable through conventional fMRI metrics. Here, we estimated the HRF by performing deconvolution on resting-state fMRI data from a longitudinal sample of 25 healthy controls scanned twice and 44 adults with obsessive-compulsive disorder (OCD) before and after 4-weeks of intensive cognitive-behavioral therapy (CBT). HRF response height, time-to-peak and full-width at half-maximum (FWHM) in OCD were abnormal before treatment and normalized after treatment in regions including the caudate. Pre-treatment HRF predicted treatment outcome (OCD symptom reduction) with 86.4% accuracy, using machine learning. Pre-treatment HRF response height in the caudate head and time-to-peak in the caudate tail were top-predictors of treatment response. Time-to-peak in the caudate tail, a region not typically identified in OCD studies using conventional fMRI activation or connectivity measures, may carry novel importance. Additionally, pre-treatment response height in caudate head predicted post-treatment OCD severity (R = -0.48, P = 0.001), and was associated with treatment-related OCD severity changes (R = -0.44, P = 0.0028), underscoring its relevance. With HRF being a reliable marker sensitive to brain function, OCD pathology, and intervention-related changes, these results could guide future studies towards novel discoveries not possible through conventional fMRI approaches like standard BOLD activation or connectivity.

Approaches to inferring multi-regional interactions from simultaneous population recordings: Inferring multi-regional interactions from simultaneous population recordings

Most past studies of neural representations and dynamics have focused on recordings from single brain areas. However, growing evidence of brain-wide, parallel representations of cognitive variables suggests that analyzing neural representations and dynamics in individual brain areas can benefit from understanding the context of multi-regional interactions that support them. Moreover, perturbation experiments revealed that the manner in which these parallel representations interact with each other can differ dramatically across different pairs of brain areas. Recent advances in recording technology offer a potentially powerful substrate to study how multi-regional interactions coordinate neural representations in individual brain areas and dictate behavior on a single-trial basis through simultaneous recordings of multiple brain areas. We review pragmatic approaches to studying multi-regional interactions and illustrate them in the concrete context of a rodent delayed response task paradigm.

Estimation and validation of individualized dynamic brain models with resting state fMRI

A key challenge for neuroscience is to develop generative, causal models of the human nervous system in an individualized, data-driven manner. Previous initiatives have either constructed biologically-plausible models that are not constrained by individual-level human brain activity or used data-driven statistical characterizations of individuals that are not mechanistic. We aim to bridge this gap through the development of a new modeling approach termed Mesoscale Individualized Neurodynamic (MINDy) modeling, wherein we fit nonlinear dynamical systems models directly to human brain imaging data. The MINDy framework is able to produce these data-driven network models for hundreds to thousands of interacting brain regions in just 1-3 ​min per subject. We demonstrate that the models are valid, reliable, and robust. We show that MINDy models are predictive of individualized patterns of resting-state brain dynamical activity. Furthermore, MINDy is better able to uncover the mechanisms underlying individual differences in resting state activity than functional connectivity methods.

Anxiety and Stress Alter Decision-Making Dynamics and Causal Amygdala-Dorsolateral Prefrontal Cortex Circuits During Emotion Regulation in Children

BACKGROUND

Anxiety and stress reactivity are risk factors for the development of affective disorders. However, the behavioral and neurocircuit mechanisms that potentiate maladaptive emotion regulation are poorly understood. Neuroimaging studies have implicated the amygdala and dorsolateral prefrontal cortex (DLPFC) in emotion regulation, but how anxiety and stress alter their context-specific causal circuit interactions is not known. Here, we use computational modeling to inform affective pathophysiology, etiology, and neurocircuit targets for early intervention.

METHODS

Forty-five children (10-11 years of age; 25 boys) reappraised aversive stimuli during functional magnetic resonance imaging scanning. Clinical measures of anxiety and stress were acquired for each child. Drift-diffusion modeling of behavioral data and causal circuit analysis of functional magnetic resonance imaging data, with a National Institute of Mental Health Research Domain Criteria approach, were used to characterize latent behavioral and neurocircuit decision-making dynamics driving emotion regulation.

RESULTS

Children successfully reappraised negative responses to aversive stimuli. Drift-diffusion modeling revealed that emotion regulation was characterized by increased initial bias toward positive reactivity during viewing of aversive stimuli and increased drift rate, which captured evidence accumulation during emotion evaluation. Crucially, anxiety and stress reactivity impaired latent behavioral dynamics associated with reappraisal and decision making. Anxiety and stress increased dynamic casual influences from the right amygdala to DLPFC. In contrast, DLPFC, but not amygdala, reactivity was correlated with evidence accumulation and decision making during emotion reappraisal.

CONCLUSIONS

Our findings provide new insights into how anxiety and stress in children impact decision making and amygdala-DLPFC signaling during emotion regulation, and uncover latent behavioral and neurocircuit mechanisms of early risk for psychopathology.

Tilting Together: An Information-Theoretic Characterization of Behavioral Roles in Rhythmic Dyadic Interaction

Every joint collaborative physical activity performed by a group of people, e.g., carrying a table, typically leads to the emergence of spatiotemporal coordination of individual motor behavior. Such interpersonal coordination can arise solely based on the observation of the partners' and/or the object's movements, without the presence of verbal communication. In this paper, we investigate how the social coupling between two individuals in a collaborative task translates into measured objective and subjective performance indicators recorded in two different studies. We analyse the trends in the dyadic interrelationship based on the information-theoretic measure of transfer entropy and identify emerging leader-follower roles. In our experimental paradigm, the actions of the pair of subjects are continuously and seamlessly fused, resulting in a joint control of an object simulated on a tablet computer. Subjects need to synchronize their movements with a 90° phase difference in order to keep the object (a ball) rotating precisely on a predefined circular or elliptic trajectory on a tablet device. Results demonstrate how the identification of causal dependencies in this social interaction task could reveal specific trends in human behavior and provide insights into the emergence of social sensorimotor contingencies.

Decoding Task-Specific Cognitive States with Slow, Directed Functional Networks in the Human Brain

Flexible functional interactions among brain regions mediate critical cognitive functions. Such interactions can be measured using functional magnetic resonance imaging (fMRI) data either with instantaneous (zero-lag) or lag-based (time-lagged) functional connectivity. Because the fMRI hemodynamic response is slow, and is sampled at a timescale (seconds) several orders of magnitude slower than the underlying neural dynamics (milliseconds), simulation studies have shown that lag-based fMRI functional connectivity, measured with approaches like Granger–Geweke causality (GC), provides spurious and unreliable estimates of underlying neural interactions. Experimental verification of this claim is challenging because neural ground truth connectivity is often unavailable concurrently with fMRI recordings. Here we demonstrate that, despite these widely held caveats, GC networks estimated from fMRI recordings contain useful information for classifying task-specific cognitive states. We estimated instantaneous and lag-based GC functional connectivity networks using fMRI data from 1000 participants (Human Connectome Project database). A linear classifier, trained on either instantaneous or lag-based GC, reliably discriminated among seven different task and resting brain states, with >80% cross-validation accuracy. With network simulations, we demonstrate that instantaneous and lag-based GC exploited interactions at fast and slow timescales, respectively, to achieve robust classification. With human fMRI data, instantaneous and lag-based GC identified complementary, task–core networks. Finally, variations in GC connectivity explained inter-individual variations in a variety of cognitive scores. Our findings show that instantaneous and lag-based methods reveal complementary aspects of functional connectivity in the brain, and suggest that slow, directed functional interactions, estimated with fMRI, may provide useful markers of behaviorally relevant cognitive states.

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.

Delayed correlations improve the reconstruction of the brain connectome

The brain works as a large-scale complex network, known as the connectome. The strength of the connections between two brain regions in the connectome is commonly estimated by calculating the correlations between their patterns of activation. This approach relies on the assumption that the activation of connected regions occurs together and at the same time. However, there are delays between the activation of connected regions due to excitatory and inhibitory connections. Here, we propose a method to harvest this additional information and reconstruct the structural brain connectome using delayed correlations. This delayed-correlation method correctly identifies 70% to 80% of connections of simulated brain networks, compared to only 5% to 25% of connections detected by the standard methods; this result is robust against changes in the network parameters (small-worldness, excitatory vs. inhibitory connection ratio, weight distribution) and network activation dynamics. The delayed-correlation method predicts more accurately both the global network properties (characteristic path length, global efficiency, clustering coefficient, transitivity) and the nodal network properties (nodal degree, nodal clustering, nodal global efficiency), particularly at lower network densities. We obtain similar results in networks derived from animal and human data. These results suggest that the use of delayed correlations improves the reconstruction of the structural brain connectome and open new possibilities for the analysis of the brain connectome, as well as for other types of networks.

Increasing robustness of pairwise methods for effective connectivity in magnetic resonance imaging by using fractional moment series of BOLD signal distributions

Estimating causal interactions in the brain from functional magnetic resonance imaging (fMRI) data remains a challenging task. Multiple studies have demonstrated that all current approaches to determine direction of connectivity perform poorly when applied to synthetic fMRI datasets. Recent advances in this field include methods for pairwise inference, which involve creating a sparse connectome in the first step, and then using a classifier in order to determine the directionality of connection between every pair of nodes in the second step. In this work, we introduce an advance to the second step of this procedure, by building a classifier based on fractional moments of the BOLD distribution combined into cumulants. The classifier is trained on datasets generated under the dynamic causal modeling (DCM) generative model. The directionality is inferred based on statistical dependencies between the two-node time series, for example, by assigning a causal link from time series of low variance to time series of high variance. Our approach outperforms or performs as well as other methods for effective connectivity when applied to the benchmark datasets. Crucially, it is also more resilient to confounding effects such as differential noise level across different areas of the connectome.

Dynamics of Segregation and Integration in Directional Brain Networks: Illustration in Soldiers With PTSD and Neurotrauma

Brain functioning relies on various segregated/specialized neural regions functioning as an integrated-interconnected network (i.e., metastability). Various psychiatric and neurologic disorders are associated with aberrant functioning of these brain networks. In this study, we present a novel framework integrating the strength and temporal variability of metastability in brain networks. We demonstrate that this approach provides novel mechanistic insights which enables better imaging-based predictions. Using whole-brain resting-state fMRI and a graph-theoretic framework, we integrated strength and temporal-variability of complex-network properties derived from effective connectivity networks, obtained from 87 U.S. Army soldiers consisting of healthy combat controls (n = 28), posttraumatic stress disorder (PTSD; n = 17), and PTSD with comorbid mild-traumatic brain injury (mTBI; n = 42). We identified prefrontal dysregulation of key subcortical and visual regions in PTSD/mTBI, with all network properties exhibiting lower variability over time, indicative of poorer flexibility. Larger impairment in the prefrontal-subcortical pathway but not prefrontal-visual pathway differentiated comorbid PTSD/mTBI from the PTSD group. Network properties of the prefrontal-subcortical pathway also had significant association (R ² = 0.56) with symptom severity and neurocognitive performance; and were also found to possess high predictive ability (81.4% accuracy in classifying the disorders, explaining 66-72% variance in symptoms), identified through machine learning. Our framework explained 13% more variance in behaviors compared to the conventional framework. These novel insights and better predictions were made possible by our novel framework using static and time-varying network properties in our three-group scenario, advancing the mechanistic understanding of PTSD and comorbid mTBI. Our contribution has wide-ranging applications for network-level characterization of healthy brains as well as mental disorders.

Identifying nonlinear dynamical systems via generative recurrent neural networks with applications to fMRI

A major tenet in theoretical neuroscience is that cognitive and behavioral processes are ultimately implemented in terms of the neural system dynamics. Accordingly, a major aim for the analysis of neurophysiological measurements should lie in the identification of the computational dynamics underlying task processing. Here we advance a state space model (SSM) based on generative piecewise-linear recurrent neural networks (PLRNN) to assess dynamics from neuroimaging data. In contrast to many other nonlinear time series models which have been proposed for reconstructing latent dynamics, our model is easily interpretable in neural terms, amenable to systematic dynamical systems analysis of the resulting set of equations, and can straightforwardly be transformed into an equivalent continuous-time dynamical system. The major contributions of this paper are the introduction of a new observation model suitable for functional magnetic resonance imaging (fMRI) coupled to the latent PLRNN, an efficient stepwise training procedure that forces the latent model to capture the 'true' underlying dynamics rather than just fitting (or predicting) the observations, and of an empirical measure based on the Kullback-Leibler divergence to evaluate from empirical time series how well this goal of approximating the underlying dynamics has been achieved. We validate and illustrate the power of our approach on simulated 'ground-truth' dynamical systems as well as on experimental fMRI time series, and demonstrate that the learnt dynamics harbors task-related nonlinear structure that a linear dynamical model fails to capture. Given that fMRI is one of the most common techniques for measuring brain activity non-invasively in human subjects, this approach may provide a novel step toward analyzing aberrant (nonlinear) dynamics for clinical assessment or neuroscientific research.

Variational Bayesian Approach for Causality and Contemporaneous Correlation Features Inference in Industrial Process Data

In this paper, a hybrid model is proposed to simultaneously mine causal connections and features responsible for contemporaneous correlations in a multivariate process. The model is developed by combining the vector auto-regressive exogenous model and the factor analysis model. The parameters of the resulting model are regularized using the hierarchical prior distributions for pruning insignificant/irrelevant ones from the model. It is then estimated under the variational Bayesian expectation maximization framework. The estimation is initiated with a complex model which is then systematically reduced to a simpler model that retains only the parameters corresponding to significant causal connections and contemporaneous correlations. Model reduction is carried out through a series of deterministic jumps from complex models to simpler models using a relevance criterion. The approach is illustrated with a number of simulated examples and an industrial case study.

Mapping distinct timescales of functional interactions among brain networks

Brain processes occur at various timescales, ranging from milliseconds (neurons) to minutes and hours (behavior). Characterizing functional coupling among brain regions at these diverse timescales is key to understanding how the brain produces behavior. Here, we apply instantaneous and lag-based measures of conditional linear dependence, based on Granger-Geweke causality (GC), to infer network connections at distinct timescales from functional magnetic resonance imaging (fMRI) data. Due to the slow sampling rate of fMRI, it is widely held that GC produces spurious and unreliable estimates of functional connectivity when applied to fMRI data. We challenge this claim with simulations and a novel machine learning approach. First, we show, with simulated fMRI data, that instantaneous and lag-based GC identify distinct timescales and complementary patterns of functional connectivity. Next, we analyze fMRI scans from 500 subjects and show that a linear classifier trained on either instantaneous or lag-based GC connectivity reliably distinguishes task versus rest brain states, with ~80-85% cross-validation accuracy. Importantly, instantaneous and lag-based GC exploit markedly different spatial and temporal patterns of connectivity to achieve robust classification. Our approach enables identifying functionally connected networks that operate at distinct timescales in the brain.

Advances in techniques for imposing reciprocity in brain-behavior relations

To better understand human behavior, the emerging field of model-based cognitive neuroscience seeks to anchor psychological theory to the biological substrate from which behavior originates: the brain. Despite complex dynamics, many researchers in this field have demonstrated that fluctuations in brain activity can be related to fluctuations in components of cognitive models, which instantiate psychological theories. In this review, we discuss a number of approaches for relating brain activity to cognitive models, and expand on a framework for imposing reciprocity in the inference of mental operations from the combination of brain and behavioral data.

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.

Abnormal functional integration across core brain networks in migraine without aura

Background

As a complex subjective experience, pain processing may be related to functional integration among intrinsic connectivity networks of migraine patients without aura. However, few study focused on the pattern alterations in the intrinsic connectivity networks of migraine patients without aura.

Results

Thirty-one migraine patients without aura and 31 age- and education-matched healthy controls participated in this study. After identifying the default mode network, central executive network and salience network as core intrinsic connectivity networks by using independent component analysis, functional connectivity, and effective connectivity during the resting state were used to investigate the abnormalities in intrinsic connectivity network interactions. Migraine patients without aura showed decreased functional connectivity among intrinsic connectivity networks compared with healthy controls. The strength of causal influences from the right frontoinsular cortex to the right anterior cingulate cortex became weaker, and the right frontoinsular cortex to the right medial prefrontal cortex became stronger in migraine patients without aura.

Conclusions

These changes suggested that the salience network may play a major role in the pathophysiological features of migraine patients without aura and helped us to synthesize previous findings into an aberrant network dynamical framework.

Directed functional connectivity using dynamic graphical models

There are a growing number of neuroimaging methods that model spatio-temporal patterns of brain activity to allow more meaningful characterizations of brain networks. This paper proposes dynamic graphical models (DGMs) for dynamic, directed functional connectivity. DGMs are a multivariate graphical model with time-varying coefficients that describe instantaneous directed relationships between nodes. A further benefit of DGMs is that networks may contain loops and that large networks can be estimated. We use network simulations and human resting-state fMRI (N = 500) to investigate the validity and reliability of the estimated networks. We simulate systematic lags of the hemodynamic response at different brain regions to investigate how these lags potentially bias directionality estimates. In the presence of such lag confounds (0.4-0.8 s offset between connected nodes), our method has a sensitivity of 72%-77% to detect the true direction. Stronger lag confounds have reduced sensitivity, but do not increase false positives (i.e., directionality estimates of the opposite direction). In human resting-state fMRI, the default mode network has consistent influence on the cerebellar, the limbic and the auditory/temporal networks. We also show a consistent reciprocal relationship between the visual medial and visual lateral network. Finally, we apply the method in a small mouse fMRI sample and discover a highly plausible relationship between areas in the hippocampus feeding into the cingulate cortex. We provide a computationally efficient implementation of DGM as a free software package for R.

Sparse Estimation of Resting-State Effective Connectivity From fMRI Cross-Spectra

In functional magnetic resonance imaging (fMRI), functional connectivity is conventionally characterized by correlations between fMRI time series, which are intrinsically undirected measures of connectivity. Yet, some information about the directionality of network connections can nevertheless be extracted from the matrix of pairwise temporal correlations between all considered time series, when expressed in the frequency-domain as a cross-spectral density matrix. Using a sparsity prior, it then becomes possible to determine a unique directed network topology that best explains the observed undirected correlations, without having to rely on temporal precedence relationships that may not be valid in fMRI. Applying this method on simulated data with 100 nodes yielded excellent retrieval of the underlying directed networks under a wide variety of conditions. Importantly, the method did not depend on temporal precedence to establish directionality, thus reducing susceptibility to hemodynamic variability. The computational efficiency of the algorithm was sufficient to enable whole-brain estimations, thus circumventing the problem of missing nodes that otherwise occurs in partial-brain analyses. Applying the method to real resting-state fMRI data acquired with a high temporal resolution, the inferred networks showed good consistency with structural connectivity obtained from diffusion tractography in the same subjects. Interestingly, this agreement could also be seen when considering high-frequency rather than low-frequency connectivity (average correlation: r = 0.26 for f < 0.3 Hz, r = 0.43 for 0.3 < f < 5 Hz). Moreover, this concordance was significantly better (p < 0.05) than for networks obtained with conventional functional connectivity based on correlations (average correlation r = 0.18). The presented methodology thus appears to be well-suited for fMRI, particularly given its lack of explicit dependence on temporal lag structure, and is readily applicable to whole-brain effective connectivity estimation.

Determining Excitatory and Inhibitory Neuronal Activity from Multimodal fMRI Data Using a Generative Hemodynamic Model

Hemodynamic responses, in general, and the blood oxygenation level-dependent (BOLD) fMRI signal, in particular, provide an indirect measure of neuronal activity. There is strong evidence that the BOLD response correlates well with post-synaptic changes, induced by changes in the excitatory and inhibitory (E-I) balance between active neuronal populations. Typical BOLD responses exhibit transients, such as the early-overshoot and post-stimulus undershoot, that can be linked to transients in neuronal activity, but they can also result from vascular uncoupling between cerebral blood flow (CBF) and venous cerebral blood volume (venous CBV). Recently, we have proposed a novel generative hemodynamic model of the BOLD signal within the dynamic causal modeling framework, inspired by physiological observations, called P-DCM (Havlicek et al., 2015). We demonstrated the generative model's ability to more accurately model commonly observed neuronal and vascular transients in single regions but also effective connectivity between multiple brain areas (Havlicek et al., 2017b). In this paper, we additionally demonstrate the versatility of the generative model to jointly explain dynamic relationships between neuronal and hemodynamic physiological variables underlying the BOLD signal using multi-modal data. For this purpose, we utilized three distinct data-sets of experimentally induced responses in the primary visual areas measured in human, cat, and monkey brain, respectively: (1) CBF and BOLD responses; (2) CBF, total CBV, and BOLD responses (Jin and Kim, 2008); and (3) positive and negative neuronal and BOLD responses (Shmuel et al., 2006). By fitting the generative model to the three multi-modal experimental data-sets, we showed that the presence or absence of dynamic features in the BOLD signal is not an unambiguous indication of presence or absence of those features on the neuronal level. Nevertheless, the generative model that takes into account the dynamics of the physiological mechanisms underlying the BOLD response allowed dissociating neuronal from vascular transients and deducing excitatory and inhibitory neuronal activity time-courses from BOLD data alone and from multi-modal data.

Network Configurations in the Human Brain Reflect Choice Bias during Rapid Face Processing

Network interactions are likely to be instrumental in processes underlying rapid perception and cognition. Specifically, high-level and perceptual regions must interact to balance pre-existing models of the environment with new incoming stimuli. Simultaneous electroencephalography (EEG) and fMRI (EEG/fMRI) enables temporal characterization of brain-network interactions combined with improved anatomical localization of regional activity. In this paper, we use simultaneous EEG/fMRI and multivariate dynamical systems (MDS) analysis to characterize network relationships between constitute brain areas that reflect a subject's choice for a face versus nonface categorization task. Our simultaneous EEG and fMRI analysis on 21 human subjects (12 males, 9 females) identifies early perceptual and late frontal subsystems that are selective to the categorical choice of faces versus nonfaces. We analyze the interactions between these subsystems using an MDS in the space of the BOLD signal. Our main findings show that differences between face-choice and house-choice networks are seen in the network interactions between the early and late subsystems, and that the magnitude of the difference in network interaction positively correlates with the behavioral false-positive rate of face choices. We interpret this to reflect the role of saliency and expectations likely encoded in frontal "late" regions on perceptual processes occurring in "early" perceptual regions.SIGNIFICANCE STATEMENT Our choices are affected by our biases. In visual perception and cognition such biases can be commonplace and quite curious-e.g., we see a human face when staring up at a cloud formation or down at a piece of toast at the breakfast table. Here we use multimodal neuroimaging and dynamical systems analysis to measure whole-brain spatiotemporal dynamics while subjects make decisions regarding the type of object they see in rapidly flashed images. We find that the degree of interaction in these networks accounts for a substantial fraction of our bias to see faces. In general, our findings illustrate how the properties of spatiotemporal networks yield insight into the mechanisms of how we form decisions.

Identifying disease foci from static and dynamic effective connectivity networks: Illustration in soldiers with trauma

Brain connectivity studies report group differences in pairwise connection strengths. While informative, such results are difficult to interpret since our understanding of the brain relies on region-based properties, rather than on connection information. Given that large disruptions in the brain are often caused by a few pivotal sources, we propose a novel framework to identify the sources of functional disruption from effective connectivity networks. Our approach integrates static and time-varying effective connectivity modeling in a probabilistic framework, to identify aberrant foci and the corresponding aberrant connectomics network. Using resting-state fMRI, we illustrate the utility of this novel approach in U.S. Army soldiers (N = 87) with posttraumatic stress disorder (PTSD), mild traumatic brain injury (mTBI) and combat controls. Additionally, we employed machine-learning classification to identify those significant connectivity features that possessed high predictive ability. We identified three disrupted foci (middle frontal gyrus [MFG], insula, hippocampus), and an aberrant prefrontal-subcortical-parietal network of information flow. We found the MFG to be the pivotal focus of network disruption, with aberrant strength and temporal-variability of effective connectivity to the insula, amygdala and hippocampus. These connectivities also possessed high predictive ability (giving a classification accuracy of 81%); and they exhibited significant associations with symptom severity and neurocognitive functioning. In summary, dysregulation originating in the MFG caused elevated and temporally less-variable connectivity in subcortical regions, followed by a similar effect on parietal memory-related regions. This mechanism likely contributes to the reduced control over traumatic memories leading to re-experiencing, hyperarousal and flashbacks observed in soldiers with PTSD and mTBI. Hum Brain Mapp 39:264-287, 2018. © 2017 Wiley Periodicals, Inc.

The impact of hemodynamic variability and signal mixing on the identifiability of effective connectivity structures in BOLD fMRI

PURPOSE

Multiple computational studies have demonstrated that essentially all current analytical approaches to determine effective connectivity perform poorly when applied to synthetic functional Magnetic Resonance Imaging (fMRI) datasets. In this study, we take a theoretical approach to investigate the potential factors facilitating and hindering effective connectivity research in fMRI.

MATERIALS AND METHODS

In this work, we perform a simulation study with use of Dynamic Causal Modeling generative model in order to gain new insights on the influence of factors such as the slow hemodynamic response, mixed signals in the network and short time series, on the effective connectivity estimation in fMRI studies.

RESULTS

First, we perform a Linear Discriminant Analysis study and find that not the hemodynamics itself but mixed signals in the neuronal networks are detrimental to the signatures of distinct connectivity patterns. This result suggests that for statistical methods (which do not involve lagged signals), deconvolving the BOLD responses is not necessary, but at the same time, functional parcellation into Regions of Interest (ROIs) is essential. Second, we study the impact of hemodynamic variability on the inference with use of lagged methods. We find that the local hemodynamic variability provide with an upper bound on the success rate of the lagged methods. Furthermore, we demonstrate that upsampling the data to TRs lower than the TRs in state-of-the-art datasets does not influence the performance of the lagged methods.

CONCLUSIONS

Factors such as background scale-free noise and hemodynamic variability have a major impact on the performance of methods for effective connectivity research in functional Magnetic Resonance Imaging.

Exploring connectivity with large-scale Granger causality on resting-state functional MRI

BACKGROUND

Large-scale Granger causality (lsGC) is a recently developed, resting-state functional MRI (fMRI) connectivity analysis approach that estimates multivariate voxel-resolution connectivity. Unlike most commonly used multivariate approaches, which establish coarse-resolution connectivity by aggregating voxel time-series avoiding an underdetermined problem, lsGC estimates voxel-resolution, fine-grained connectivity by incorporating an embedded dimension reduction.

NEW METHOD

We investigate application of lsGC on realistic fMRI simulations, modeling smoothing of neuronal activity by the hemodynamic response function and repetition time (TR), and empirical resting-state fMRI data. Subsequently, functional subnetworks are extracted from lsGC connectivity measures for both datasets and validated quantitatively. We also provide guidelines to select lsGC free parameters.

RESULTS

Results indicate that lsGC reliably recovers underlying network structure with area under receiver operator characteristic curve (AUC) of 0.93 at TR=1.5s for a 10-min session of fMRI simulations. Furthermore, subnetworks of closely interacting modules are recovered from the aforementioned lsGC networks. Results on empirical resting-state fMRI data demonstrate recovery of visual and motor cortex in close agreement with spatial maps obtained from (i) visuo-motor fMRI stimulation task-sequence (Accuracy=0.76) and (ii) independent component analysis (ICA) of resting-state fMRI (Accuracy=0.86).

COMPARISON WITH EXISTING METHOD(S)

Compared with conventional Granger causality approach (AUC=0.75), lsGC produces better network recovery on fMRI simulations. Furthermore, it cannot recover functional subnetworks from empirical fMRI data, since quantifying voxel-resolution connectivity is not possible as consequence of encountering an underdetermined problem.

CONCLUSIONS

Functional network recovery from fMRI data suggests that lsGC gives useful insight into connectivity patterns from resting-state fMRI at a multivariate voxel-resolution.

Extraversion and neuroticism related to the resting-state effective connectivity of amygdala

The amygdala plays a key role in emotion processing. Its functional connectivity with other brain regions has been extensively demonstrated to be associated with extraversion and neuroticism. However, how the amygdala affects other regions and is affected by others within these connectivity patterns associated with extraversion and neuroticism remains unclear. To address this issue, we investigated the effective connectivity of the amygdala using Granger causality analysis on the resting-state functional magnetic resonance imaging data of 70 participants. Results showed that extraversion was positively correlated with the influence from the right inferior occipital gyrus (IOG) to the left amygdala, and from the bilateral IOG to the right amygdala; such result may represent the neural correlates of social interactions in extraverts. Conversely, neuroticism was associated with an increased influence from right amygdala to right middle frontal gyrus and a decreased influence from right precuneus to right amygdala. This influence might affect the modulations of cognitive regulation function and self-referential processes in neurotic individuals. These findings highlight the importance of the causal influences of amygdala in explaining the individual differences in extraversion and neuroticism, and offer further insights into the specific neural networks underlying personality.

Understanding the Impact of Stroke on Brain Motor Function: A Hierarchical Bayesian Approach

Stroke is a disturbance in blood supply to the brain resulting in the loss of brain functions, particularly motor function. A study was conducted by the UCI Neurorehabilitation Lab to investigate the impact of stroke on motor-related brain regions. Functional MRI (fMRI) data were collected from stroke patients and healthy controls while the subjects performed a simple motor task. In addition to affecting local neuronal activation strength, stroke might also alter communications (i.e., connectivity) between brain regions. We develop a hierarchical Bayesian modeling approach for the analysis of multi-subject fMRI data that allows us to explore brain changes due to stroke. Our approach simultaneously estimates activation and condition-specific connectivity at the group level, and provides estimates for region/subject-specific hemodynamic response functions. Moreover, our model uses spike and slab priors to allow for direct posterior inference on the connectivity network. Our results indicate that motor-control regions show greater activation in the unaffected hemisphere and the midline surface in stroke patients than those same regions in healthy controls during the simple motor task. We also note increased connectivity within secondary motor regions in stroke subjects. These findings provide insight into altered neural correlates of movement in subjects who suffered a stroke.

Experimental Validation of Dynamic Granger Causality for Inferring Stimulus-Evoked Sub-100 ms Timing Differences from fMRI

Decoding the sequential flow of events in the human brain non-invasively is critical for gaining a mechanistic understanding of brain function. In this study, we propose a method based on dynamic Granger causality analysis to measure timing differences in brain responses from fMRI. We experimentally validate this method by detecting sub-100 ms timing differences in fMRI responses obtained from bilateral visual cortex using fast sampling, ultra-high field and an event-related visual hemifield paradigm with known timing difference between the hemifields. Classical Granger causality was previously shown to be able to detect sub-100 ms timing differences in the visual cortex. Since classical Granger causality does not differentiate between spontaneous and stimulus-evoked responses, dynamic Granger causality has been proposed as an alternative, thereby necessitating its experimental validation. In addition to detecting timing differences as low as 28 ms using dynamic Granger causality, the significance of the inference from our method increased with increasing delay both in simulations and experimental data. Therefore, it provides a methodology for understanding mental chronometry from fMRI in a data-driven way.

Joint Data Filtering and Labeling Using Gaussian Processes and Alternating Direction Method of Multipliers

Sequence labeling aims at assigning a label to every sample of a signal (or pixel of an image) while considering the sequentiality (or vicinity) of the samples. To perform this task, many works in the literature first filter and then label the data. Unfortunately, the filtering, which is performed independently from the labeling, is far from optimal and frequently makes the latter task harder. In this paper, a novel approach that trains a Gaussian process classifier and estimates the coefficients of an optimal filter jointly is presented. The new approach, based on Bayesian modeling and alternating direction method of multipliers (ADMMs) optimization, performs both tasks simultaneously. All unknowns are treated as stochastic variables, which are estimated using variational inference and filtering and labeling are linked with the use of ADMM. In the experimental section, synthetic and real experiments are presented to compare the proposed method with other existing approaches.

Atypical effective connectivity of thalamo-cortical circuits in autism spectrum disorder

Autism spectrum disorder (ASD) is a neurodevelopment disorder characterized by atypical connectivity within and across multiple brain systems. We aimed to explore information transmission from the sensory periphery to information processing centers of the brain across thalamo-cortical circuits in ASD. A large multicenter dataset from the autism brain imaging data exchange was utilized. A thalamus template derived from the Automatic Anatomic Labeling atlas was subdivided into six subregions corresponding to six cortical regions using a "winner-takes-all" strategy. Granger causality analysis (GCA) was then applied to calculate effective connectivity from subregions of the thalamus to the corresponding cortical regions. Results demonstrate reduced effective connectivity from the thalamus to left prefrontal cortex (P = 0.023), right posterior parietal cortex (P = 0.03), and bilateral temporal cortex (left: P = 0.014; right: P = 0.015) in ASD compared with healthy control (HC) participants. The GCA values of the thalamus-bilateral temporal cortex connections were significantly negatively correlated with communication scores as assessed by the autism diagnostic observation schedule in the ASD group (left: P = 0.037; right: P = 0.007). Age-related analyses showed that the strengths of the thalamus-bilateral temporal cortex connections were significantly positively correlated with age in the HC group (left: P = 0.013; right: P = 0.016), but not in the ASD group (left: P = 0.506; right: P = 0.219). These results demonstrate impaired thalamo-cortical information transmission in ASD and suggest that atypical development of thalamus-temporal cortex connections may relate to communication deficits in the disorder. Autism Res 2016, 9: 1183-1190. © 2016 International Society for Autism Research, Wiley Periodicals, Inc.

Multivariate dynamical systems-based estimation of causal brain interactions in fMRI: Group-level validation using benchmark data, neurophysiological models and human connectome project data

BACKGROUND

Causal estimation methods are increasingly being used to investigate functional brain networks in fMRI, but there are continuing concerns about the validity of these methods.

NEW METHOD

Multivariate dynamical systems (MDS) is a state-space method for estimating dynamic causal interactions in fMRI data. Here we validate MDS using benchmark simulations as well as simulations from a more realistic stochastic neurophysiological model. Finally, we applied MDS to investigate dynamic casual interactions in a fronto-cingulate-parietal control network using human connectome project (HCP) data acquired during performance of a working memory task. Crucially, since the ground truth in experimental data is unknown, we conducted novel stability analysis to determine robust causal interactions within this network.

RESULTS

MDS accurately recovered dynamic causal interactions with an area under receiver operating characteristic (AUC) above 0.7 for benchmark datasets and AUC above 0.9 for datasets generated using the neurophysiological model. In experimental fMRI data, bootstrap procedures revealed a stable pattern of causal influences from the anterior insula to other nodes of the fronto-cingulate-parietal network.

COMPARISON WITH EXISTING METHODS

MDS is effective in estimating dynamic causal interactions in both the benchmark and neurophysiological model based datasets in terms of AUC, sensitivity and false positive rates.

CONCLUSIONS

Our findings demonstrate that MDS can accurately estimate causal interactions in fMRI data. Neurophysiological models and stability analysis provide a general framework for validating computational methods designed to estimate causal interactions in fMRI. The right anterior insula functions as a causal hub during working memory.

Combining optogenetic stimulation and fMRI to validate a multivariate dynamical systems model for estimating causal brain interactions

State-space multivariate dynamical systems (MDS) (Ryali et al. 2011) and other causal estimation models are being increasingly used to identify directed functional interactions between brain regions. However, the validity and accuracy of such methods are poorly understood. Performance evaluation based on computer simulations of small artificial causal networks can address this problem to some extent, but they often involve simplifying assumptions that reduce biological validity of the resulting data. Here, we use a novel approach taking advantage of recently developed optogenetic fMRI (ofMRI) techniques to selectively stimulate brain regions while simultaneously recording high-resolution whole-brain fMRI data. ofMRI allows for a more direct investigation of causal influences from the stimulated site to brain regions activated downstream and is therefore ideal for evaluating causal estimation methods in vivo. We used ofMRI to investigate whether MDS models for fMRI can accurately estimate causal functional interactions between brain regions. Two cohorts of ofMRI data were acquired, one at Stanford University and the University of California Los Angeles (Cohort 1) and the other at the University of North Carolina Chapel Hill (Cohort 2). In each cohort, optical stimulation was delivered to the right primary motor cortex (M1). General linear model analysis revealed prominent downstream thalamic activation in Cohort 1, and caudate-putamen (CPu) activation in Cohort 2. MDS accurately estimated causal interactions from M1 to thalamus and from M1 to CPu in Cohort 1 and Cohort 2, respectively. As predicted, no causal influences were found in the reverse direction. Additional control analyses demonstrated the specificity of causal interactions between stimulated and target sites. Our findings suggest that MDS state-space models can accurately and reliably estimate causal interactions in ofMRI data and further validate their use for estimating causal interactions in fMRI. More generally, our study demonstrates that the combined use of optogenetics and fMRI provides a powerful new tool for evaluating computational methods designed to estimate causal interactions between distributed brain regions.

Patterns of effective connectivity during memory encoding and retrieval differ between patients with mild cognitive impairment and healthy older adults

Previous research has shown that there is considerable overlap in the neural networks mediating successful memory encoding and retrieval. However, little is known about how the relevant human brain regions interact during these distinct phases of memory or how such interactions are affected by memory deficits that characterize mild cognitive impairment (MCI), a condition that often precedes dementia due to Alzheimer's disease. Here we employed multivariate Granger causality analysis using autoregressive modeling of inferred neuronal time series obtained by deconvolving the hemodynamic response function from measured blood oxygenation level-dependent (BOLD) time series data, in order to examine the effective connectivity between brain regions during successful encoding and/or retrieval of object location associations in MCI patients and comparable healthy older adults. During encoding, healthy older adults demonstrated a left hemisphere dominant pattern where the inferior frontal junction, anterior intraparietal sulcus (likely involving the parietal eye fields), and posterior cingulate cortex drove activation in most left hemisphere regions and virtually every right hemisphere region tested. These regions are part of a frontoparietal network that mediates top-down cognitive control and is implicated in successful memory formation. In contrast, in the MCI patients, the right frontal eye field drove activation in every left hemisphere region examined, suggesting reliance on more basic visual search processes. Retrieval in the healthy older adults was primarily driven by the right hippocampus with lesser contributions of the right anterior thalamic nuclei and right inferior frontal sulcus, consistent with theoretical models holding the hippocampus as critical for the successful retrieval of memories. The pattern differed in MCI patients, in whom the right inferior frontal junction and right anterior thalamus drove successful memory retrieval, reflecting the characteristic hippocampal dysfunction of these patients. These findings demonstrate that neural network interactions differ markedly between MCI patients and healthy older adults. Future efforts will investigate the impact of cognitive rehabilitation of memory on these connectivity patterns.