Dynamics of large-scale fMRI networks: Deconstruct brain activity to build better models of brain function

Control energy detects discrepancies in good vs. poor readers' structural-functional coupling during a rhyming task

Neuroimaging studies have identified functional and structural brain circuits that support reading. However, much less is known about how reading-related functional dynamics are constrained by white matter structure. Network control theory proposes that cortical brain dynamics are linearly determined by the white matter connectome, using control energy to evaluate the difficulty of the transition from one cognitive state to another. Here we apply this approach to linking brain dynamics with reading ability and disability in school-age children. A total of 51 children ages 8.25 -14.6 years performed an in-scanner rhyming task in visual and auditory modalities, with orthographic (spelling) and phonological (rhyming) similarity manipulated across trials. White matter structure and fMRI activation were used conjointly to compute the control energy of the reading network in each condition relative to a null fixation state. We then tested differences in control energy across trial types, finding higher control energy during non-word trials than word trials, and during incongruent trials than congruent trials. ROI analyses further showed a dissociation between control energy of the left fusiform and superior temporal gyrus depending on stimulus modality, with higher control energy for visual modalities in fusiform and higher control energy for auditory modalities in STG. Together, this study highlights that control theory can explain variations on cognitive demands in higher-level abilities such as reading, beyond what can be inferred from either functional or structural MRI measures alone.

Score-driven exponential random graphs: A new class of time-varying parameter models for temporal networks

Motivated by the increasing abundance of data describing real-world networks that exhibit dynamical features, we propose an extension of the exponential random graph models (ERGMs) that accommodates the time variation of its parameters. Inspired by the fast-growing literature on dynamic conditional score models, each parameter evolves according to an updating rule driven by the score of the ERGM distribution. We demonstrate the flexibility of score-driven ERGMs (SD-ERGMs) as data-generating processes and filters and show the advantages of the dynamic version over the static one. We discuss two applications to temporal networks from financial and political systems. First, we consider the prediction of future links in the Italian interbank credit network. Second, we show that the SD-ERGM allows discriminating between static or time-varying parameters when used to model the U.S. Congress co-voting network dynamics.

Impact of a mindfulness-based intervention on neurobehavioral functioning and its association with large-scale brain networks in preterm young adolescents

AIM

Adolescents born very preterm (VPT; <32 weeks of gestation) face an elevated risk of executive, behavioral, and socioemotional difficulties. Evidence suggests beneficial effects of mindfulness-based intervention (MBI) on these abilities. This study seeks to investigate the association between the effects of MBI on executive, behavioral, and socioemotional functioning and reliable changes in large-scale brain networks dynamics during rest in VPT young adolescents who completed an 8-week MBI program.

METHODS

Neurobehavioral assessments and resting-state functional magnetic resonance imaging were performed before and after MBI in 32 VPT young adolescents. Neurobehavioral abilities in VPT participants were compared with full-term controls. In the VPT group, dynamic functional connectivity was extracted by using the innovation-driven coactivation patterns framework. The reliable change index was used to quantify change after MBI. A multivariate data-driven approach was used to explore associations between MBI-related changes on neurobehavioral measures and temporal brain dynamics.

RESULTS

Compared with term-born controls, VPT adolescents showed reduced executive and socioemotional functioning before MBI. After MBI, a significant improvement was observed for all measures that were previously reduced in the VPT group. The increase in executive functioning, only, was associated with reliable changes in the duration of activation of large-scale brain networks, including frontolimbic, amygdala-hippocampus, dorsolateral prefrontal, and visual networks.

CONCLUSION

The improvement in executive functioning after an MBI was associated with reliable changes in large-scale brain network dynamics during rest. These changes encompassed frontolimbic, amygdala-hippocampus, dorsolateral prefrontal, and visual networks that are related to different executive processes including self-regulation, attentional control, and attentional awareness of relevant sensory stimuli.

On null models for temporal small-worldness in brain dynamics

Brain dynamics can be modeled as a temporal brain network starting from the activity of different brain regions in functional magnetic resonance imaging (fMRI) signals. When validating hypotheses about temporal networks, it is important to use an appropriate statistical null model that shares some features with the treated empirical data. The purpose of this work is to contribute to the theory of temporal null models for brain networks by introducing the random temporal hyperbolic (RTH) graph model, an extension of the random hyperbolic (RH) graph, known in the study of complex networks for its ability to reproduce crucial properties of real-world networks. We focus on temporal small-worldness which, in the static case, has been extensively studied in real-world complex networks and has been linked to the ability of brain networks to efficiently exchange information. We compare the RTH graph model with standard null models for temporal networks and show it is the null model that best reproduces the small-worldness of resting brain activity. This ability to reproduce fundamental features of real brain networks, while adding only a single parameter compared with classical models, suggests that the RTH graph model is a promising tool for validating hypotheses about temporal brain networks.

Transient brain activity dynamics discriminate levels of consciousness during anesthesia

The awake mammalian brain is functionally organized in terms of large-scale distributed networks that are constantly interacting. Loss of consciousness might disrupt this temporal organization leaving patients unresponsive. We hypothesize that characterizing brain activity in terms of transient events may provide a signature of consciousness. For this, we analyze temporal dynamics of spatiotemporally overlapping functional networks obtained from fMRI transient activity across different anesthetics and levels of anesthesia. We first show a striking homology in spatial organization of networks between monkeys and humans, indicating cross-species similarities in resting-state fMRI structure. We then track how network organization shifts under different anesthesia conditions in macaque monkeys. While the spatial aspect of the networks is preserved, their temporal dynamics are highly affected by anesthesia. Networks express for longer durations and co-activate in an anesthetic-specific configuration. Additionally, hierarchical brain organization is disrupted with a consciousness-level-signature role of the default mode network. In conclusion, large-scale brain network temporal dynamics capture differences in anesthetic-specific consciousness-level, paving the way towards a clinical translation of these cortical signature.

Task sub-type states decoding via group deep bidirectional recurrent neural network

Decoding brain states under different cognitive tasks from functional magnetic resonance imaging (fMRI) data has attracted great attention in the neuroimaging filed. However, the well-known temporal dependency in fMRI sequences has not been fully exploited in existing studies, due to the limited temporal-modeling capacity of the backbone machine learning algorithms and rigid training sample organization strategies upon which the brain decoding methods are built. To address these limitations, we propose a novel method for fine-grain brain state decoding, namely, group deep bidirectional recurrent neural network (Group-DBRNN) model. We first propose a training sample organization strategy that consists of a group-task sample generation module and a multiple-scale random fragment strategy (MRFS) module to collect training samples that contain rich task-relevant brain activity contrast (i.e., the comparison of neural activity patterns between different tasks) and maintain the temporal dependency. We then develop a novel decoding model by replacing the unidirectional RNNs that are widely used in existing brain state decoding studies with bidirectional stacked RNNs to better capture the temporal dependency, and by introducing a multi-task interaction layer (MTIL) module to effectively model the task-relevant brain activity contrast. Our experimental results on the Human Connectome Project task fMRI dataset (7 tasks consisting of 23 task sub-type states) show that the proposed model achieves an average decoding accuracy of 94.7% over the 23 fine-grain sub-type states. Meanwhile, our extensive interpretations of the intermediate features learned in the proposed model via visualizations and quantitative assessments of their discriminability and inter-subject alignment evidence that the proposed model can effectively capture the temporal dependency and task-relevant contrast.

Cerebello-basal Ganglia Networks and Cortical Network Global Efficiency

The cerebellum (CB) and basal ganglia (BG) each have topographically distinct functional subregions that are functionally and anatomically interconnected with cortical regions through discrete thalamic loops and with each other via disynaptic connections, with previous work detailing high levels of functional connectivity between these phylogenetically ancient regions. It was posited that this CB-BG network provides support for cortical systems processing, spanning cognitive, emotional, and motor domains, implying that subcortical network measures are strongly related to cortical network measures (Bostan & Strick, 2018); however, it is currently unknown how network measures within distinct CB-BG networks relate to cortical network measures. Here, 122 regions of interest comprising cognitive and motor CB-BG networks and 7 canonical cortical resting-state were used to investigate whether the integration (quantified using global efficiency, GE) of cognitive CB-BG network (CCBN) nodes and their segregation from motor CB-BG network (MCBN) nodes is related to cortical network GE and segregation in 233 non-related, right-handed participants (Human Connectome Project-1200). CCBN GE positively correlated with GE in the default mode, motor, and auditory networks and MCBN GE positively correlated with GE in all networks, except the default mode and emotional. MCBN segregation was related to motor network segregation. These findings highlight the CB-BG network's potential role in cortical networks associated with executive function, task switching, and verbal working memory. This work has implications for understanding cortical network organization and cortical-subcortical interactions in healthy adults and may help in determining biomarkers and deciphering subcortical differences seen in disease states.

Large-scale brain network dynamics in very preterm children and relationship with socio-emotional outcomes: an exploratory study

BACKGROUND

Children born very preterm (VPT; <32 weeks' gestation) are at high risk of neurodevelopmental and behavioural difficulties associated with atypical brain maturation, including socio-emotional difficulties. The analysis of large-scale brain network dynamics during rest allows us to investigate brain functional connectivity and its association with behavioural outcomes.

METHODS

Dynamic functional connectivity was extracted by using the innovation-driven co-activation patterns framework in VPT and full-term children aged 6-9 to explore changes in spatial organisation, laterality and temporal dynamics of spontaneous large-scale brain activity (VPT, n = 28; full-term, n = 12). Multivariate analysis was used to explore potential biomarkers for socio-emotional difficulties in VPT children.

RESULTS

The spatial organisation of the 13 retrieved functional networks was comparable across groups. Dynamic features and lateralisation of network brain activity were also comparable for all brain networks. Multivariate analysis unveiled group differences in associations between dynamical functional connectivity parameters with socio-emotional abilities.

CONCLUSION

In this exploratory study, the group differences observed might reflect reduced degrees of maturation of functional architecture in the VPT group in regard to socio-emotional abilities. Dynamic features of functional connectivity could represent relevant neuroimaging markers and inform on potential mechanisms through which preterm birth leads to neurodevelopmental and behavioural disorders.

IMPACT

Spatial organisation of the retrieved resting-state networks was comparable between school-aged very preterm and full-term children. Dynamic features and lateralisation of network brain activity were also comparable across groups. Multivariate pattern analysis revealed different patterns of association between dynamical functional connectivity parameters and socio-emotional abilities in the very preterm and full-term groups. Findings suggest a reduced degree of maturation of the functional architecture in the very preterm group in association with socio-emotional abilities.

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.

Using deep clustering to improve fMRI dynamic functional connectivity analysis

Dynamic functional connectivity (dFC) analysis of resting-state fMRI data is commonly performed by calculating sliding-window correlations (SWC), followed by k-means clustering in order to assign each window to a given state. Studies using synthetic data have shown that k-means performance is highly dependent on sliding window parameters and signal-to-noise ratio. Additionally, sources of heterogeneity between subjects may affect the accuracy of group-level clustering, thus affecting measurements of dFC state temporal properties such as dwell time and fractional occupancy. This may result in spurious conclusions regarding differences between groups (e.g. when comparing a clinical population to healthy controls). Therefore, is it important to quantify the ability of k-means to estimate dFC state temporal properties when applied to cohorts of multiple subjects, and to explore ways in which clustering performance can be maximised. Here, we explore the use of dimensionality reduction methods prior to clustering in order to map high-dimensional data to a lower dimensional space, providing salient features to the subsequent clustering step. We assess the use of deep autoencoders for dimensionality reduction prior to applying k-means clustering to the encoded data. We compare this deep clustering method to dimensionality reduction using principle component analysis (PCA), uniform manifold approximation and projection (UMAP), as well as applying k-means to the original feature space using either L1 or L2 distance. We provide extensive quantitative evaluation of clustering performance using synthetic datasets, representing data from multiple heterogeneous subjects. In synthetic data we find that deep clustering gives the best performance, while other approaches are often insufficient to capture temporal properties of dFC states. We then demonstrate the application of each method to real-world data from human subjects and show that the choice of dimensionality reduction method has a significant effect on group-level measurements of state temporal properties.

Post-stroke reorganization of transient brain activity characterizes deficits and recovery of cognitive functions

Functional magnetic resonance imaging (fMRI) has been widely employed to study stroke pathophysiology. In particular, analyses of fMRI signals at rest were directed at quantifying the impact of stroke on spatial features of brain networks. However, brain networks have intrinsic time features that were, so far, disregarded in these analyses. In consequence, standard fMRI analysis failed to capture temporal imbalance resulting from stroke lesions, hence restricting their ability to reveal the interdependent pathological changes in structural and temporal network features following stroke. Here, we longitudinally analyzed hemodynamic-informed transient activity in a large cohort of stroke patients (n = 103) to assess spatial and temporal changes of brain networks after stroke. Metrics extracted from the hemodynamic-informed transient activity were replicable within- and between-individuals in healthy participants, hence supporting their robustness and their clinical applicability. While large-scale spatial patterns of brain networks were preserved after stroke, their durations were altered, with stroke subjects exhibiting a varied pattern of longer and shorter network activations compared to healthy individuals. Specifically, patients showed a longer duration in the lateral precentral gyrus and anterior cingulum, and a shorter duration in the occipital lobe and in the cerebellum. These temporal alterations were associated with white matter damage in projection and association pathways. Furthermore, they were tied to deficits in specific behavioral domains as restoration of healthy brain dynamics paralleled recovery of cognitive functions (attention, language and spatial memory), but was not significantly correlated to motor recovery. These findings underscore the critical importance of network temporal properties in dissecting the pathophysiology of brain changes after stroke, thus shedding new light on the clinical potential of time-resolved methods for fMRI analysis.

Ranking of communities in multiplex spatiotemporal models of brain dynamics

As a relatively new field, network neuroscience has tended to focus on aggregate behaviours of the brain averaged over many successive experiments or over long recordings in order to construct robust brain models. These models are limited in their ability to explain dynamic state changes in the brain which occurs spontaneously as a result of normal brain function. Hidden Markov Models (HMMs) trained on neuroimaging time series data have since arisen as a method to produce dynamical models that are easy to train but can be difficult to fully parametrise or analyse. We propose an interpretation of these neural HMMs as multiplex brain state graph models we term Hidden Markov Graph Models. This interpretation allows for dynamic brain activity to be analysed using the full repertoire of network analysis techniques. Furthermore, we propose a general method for selecting HMM hyperparameters in the absence of external data, based on the principle of maximum entropy, and use this to select the number of layers in the multiplex model. We produce a new tool for determining important communities of brain regions using a spatiotemporal random walk-based procedure that takes advantage of the underlying Markov structure of the model. Our analysis of real multi-subject fMRI data provides new results that corroborate the modular processing hypothesis of the brain at rest as well as contributing new evidence of functional overlap between and within dynamic brain state communities. Our analysis pipeline provides a way to characterise dynamic network activity of the brain under novel behaviours or conditions.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s41109-022-00454-2.

Functional Connectivity of the Brain Across Rodents and Humans

Resting-state functional magnetic resonance imaging (rs-fMRI), which measures the spontaneous fluctuations in the blood oxygen level-dependent (BOLD) signal, is increasingly utilized for the investigation of the brain's physiological and pathological functional activity. Rodents, as a typical animal model in neuroscience, play an important role in the studies that examine the neuronal processes that underpin the spontaneous fluctuations in the BOLD signal and the functional connectivity that results. Translating this knowledge from rodents to humans requires a basic knowledge of the similarities and differences across species in terms of both the BOLD signal fluctuations and the resulting functional connectivity. This review begins by examining similarities and differences in anatomical features, acquisition parameters, and preprocessing techniques, as factors that contribute to functional connectivity. Homologous functional networks are compared across species, and aspects of the BOLD fluctuations such as the topography of the global signal and the relationship between structural and functional connectivity are examined. Time-varying features of functional connectivity, obtained by sliding windowed approaches, quasi-periodic patterns, and coactivation patterns, are compared across species. Applications demonstrating the use of rs-fMRI as a translational tool for cross-species analysis are discussed, with an emphasis on neurological and psychiatric disorders. Finally, open questions are presented to encapsulate the future direction of the field.

Metastable dynamics of neural circuits and networks

Cortical neurons emit seemingly erratic trains of action potentials or "spikes," and neural network dynamics emerge from the coordinated spiking activity within neural circuits. These rich dynamics manifest themselves in a variety of patterns, which emerge spontaneously or in response to incoming activity produced by sensory inputs. In this Review, we focus on neural dynamics that is best understood as a sequence of repeated activations of a number of discrete hidden states. These transiently occupied states are termed "metastable" and have been linked to important sensory and cognitive functions. In the rodent gustatory cortex, for instance, metastable dynamics have been associated with stimulus coding, with states of expectation, and with decision making. In frontal, parietal, and motor areas of macaques, metastable activity has been related to behavioral performance, choice behavior, task difficulty, and attention. In this article, we review the experimental evidence for neural metastable dynamics together with theoretical approaches to the study of metastable activity in neural circuits. These approaches include (i) a theoretical framework based on non-equilibrium statistical physics for network dynamics; (ii) statistical approaches to extract information about metastable states from a variety of neural signals; and (iii) recent neural network approaches, informed by experimental results, to model the emergence of metastable dynamics. By discussing these topics, we aim to provide a cohesive view of how transitions between different states of activity may provide the neural underpinnings for essential functions such as perception, memory, expectation, or decision making, and more generally, how the study of metastable neural activity may advance our understanding of neural circuit function in health and disease.

Avalanche criticality in individuals, fluid intelligence, and working memory

The critical brain hypothesis suggests that efficient neural computation can be achieved through critical brain dynamics. However, the relationship between human cognitive performance and scale‐free brain dynamics remains unclear. In this study, we investigated the whole‐brain avalanche activity and its individual variability in the human resting‐state functional magnetic resonance imaging (fMRI) data. We showed that though the group‐level analysis was inaccurate because of individual variability, the subject wise scale‐free avalanche activity was significantly associated with maximal synchronization entropy of their brain activity. Meanwhile, the complexity of functional connectivity, as well as structure–function coupling, is maximized in subjects with maximal synchronization entropy. We also observed order–disorder phase transitions in resting‐state brain dynamics and found that there were longer times spent in the subcritical regime. These results imply that large‐scale brain dynamics favor the slightly subcritical regime of phase transition. Finally, we showed evidence that the neural dynamics of human participants with higher fluid intelligence and working memory scores are closer to criticality. We identified brain regions whose critical dynamics showed significant positive correlations with fluid intelligence performance and found that these regions were located in the prefrontal cortex and inferior parietal cortex, which were believed to be important nodes of brain networks underlying human intelligence. Our results reveal the possible role that avalanche criticality plays in cognitive performance and provide a simple method to identify the critical point and map cortical states on a spectrum of neural dynamics, ranging from subcriticality to supercriticality.

Large-scale functional network dynamics in human callosal agenesis: Increased subcortical involvement and preserved laterality

In the human brain, the corpus callosum is the major white-matter commissural tract enabling the transmission of sensory-motor, and higher level cognitive information between homotopic regions of the two cerebral hemispheres. Despite developmental absence (i.e., agenesis) of the corpus callosum (AgCC), functional connectivity is preserved, including interhemispheric connectivity. Subcortical structures have been hypothesised to provide alternative pathways to enable this preservation. To test this hypothesis, we used functional Magnetic Resonance Imaging (fMRI) recordings in children with AgCC and typically developing children, and a time-resolved approach to retrieve temporal characteristics of whole-brain functional networks. We observed an increased engagement of the cerebellum and amygdala/hippocampus networks in children with AgCC compared to typically developing children. There was little evidence that laterality of activation networks was affected in AgCC. Our findings support the hypothesis that subcortical structures play an essential role in the functional reconfiguration of the brain in the absence of a corpus callosum.

Brain Connectivity Studies on Structure-Function Relationships: A Short Survey with an Emphasis on Machine Learning

This short survey reviews the recent literature on the relationship between the brain structure and its functional dynamics. Imaging techniques such as diffusion tensor imaging (DTI) make it possible to reconstruct axonal fiber tracks and describe the structural connectivity (SC) between brain regions. By measuring fluctuations in neuronal activity, functional magnetic resonance imaging (fMRI) provides insights into the dynamics within this structural network. One key for a better understanding of brain mechanisms is to investigate how these fast dynamics emerge on a relatively stable structural backbone. So far, computational simulations and methods from graph theory have been mainly used for modeling this relationship. Machine learning techniques have already been established in neuroimaging for identifying functionally independent brain networks and classifying pathological brain states. This survey focuses on methods from machine learning, which contribute to our understanding of functional interactions between brain regions and their relation to the underlying anatomical substrate.

Structural control energy of resting-state functional brain states reveals less cost-effective brain dynamics in psychosis vulnerability

How the brain's white-matter anatomy constrains brain activity is an open question that might give insights into the mechanisms that underlie mental disorders such as schizophrenia. Chromosome 22q11.2 deletion syndrome (22q11DS) is a neurodevelopmental disorder with an extremely high risk for psychosis providing a test case to study developmental aspects of schizophrenia. In this study, we used principles from network control theory to probe the implications of aberrant structural connectivity for the brain's functional dynamics in 22q11DS. We retrieved brain states from resting-state functional magnetic resonance images of 78 patients with 22q11DS and 85 healthy controls. Then, we compared them in terms of persistence control energy; that is, the control energy that would be required to persist in each of these states based on individual structural connectivity and a dynamic model. Persistence control energy was altered in a broad pattern of brain states including both energetically more demanding and less demanding brain states in 22q11DS. Further, we found a negative relationship between persistence control energy and resting-state activation time, which suggests that the brain reduces energy by spending less time in energetically demanding brain states. In patients with 22q11DS, this behavior was less pronounced, suggesting a deficiency in the ability to reduce energy through brain activation. In summary, our results provide initial insights into the functional implications of altered structural connectivity in 22q11DS, which might improve our understanding of the mechanisms underlying the disease.

Brain Density Clustering Analysis: A New Approach to Brain Functional Dynamics

BACKGROUND

A number of studies in recent years have explored whole-brain dynamic connectivity using pairwise approaches. There has been less focus on trying to analyze brain dynamics in higher dimensions over time.

METHODS

We introduce a new approach that analyzes time series trajectories to identify high traffic nodes in a high dimensional space. First, functional magnetic resonance imaging (fMRI) data are decomposed using spatial ICA to a set of maps and their associated time series. Next, density is calculated for each time point and high-density points are clustered to identify a small set of high traffic nodes. We validated our method using simulations and then implemented it on a real data set.

RESULTS

We present a novel approach that captures dynamics within a high dimensional space and also does not use any windowing in contrast to many existing approaches. The approach enables one to characterize and study the time series in a potentially high dimensional space, rather than looking at each component pair separately. Our results show that schizophrenia patients have a lower dynamism compared to healthy controls. In addition, we find patients spend more time in nodes associated with the default mode network and less time in components strongly correlated with auditory and sensorimotor regions. Interestingly, we also found that subjects oscillate between state pairs that show opposite spatial maps, suggesting an oscillatory pattern.

CONCLUSION

Our proposed method provides a novel approach to analyze the data in its native high dimensional space and can possibly provide new information that is undetectable using other methods.

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.

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.

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

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

Normalization of aberrant pretherapeutic dynamic functional connectivity of extrastriate visual system in patients who underwent thalamotomy with stereotactic radiosurgery for essential tremor: a resting-state functional MRI study

OBJECTIVE

The tremor circuitry has commonly been hypothesized to be driven by one or multiple pacemakers within the cerebello-thalamo-cortical pathway, including the cerebellum, contralateral motor thalamus, and primary motor cortex. However, previous studies, using multiple methodologies, have advocated that tremor could be influenced by changes within the right extrastriate cortex, at both the structural and functional level. The purpose of this work was to evaluate the role of the extrastriate cortex in tremor generation and further arrest after left unilateral stereotactic radiosurgery thalamotomy (SRS-T).

METHODS

The authors considered 12 healthy controls (HCs, group 1); 15 patients with essential tremor (ET, right-sided, drug-resistant; group 2) before left unilateral SRS-T; and the same 15 patients (group 3) 1 year after the intervention, to account for delayed effects. Blood oxygenation level-dependent functional MRI during resting state was used to characterize the dynamic interactions of the right extrastriate cortex, comparing HC subjects against patients with ET before and 1 year after SRS-T. In particular, the authors applied coactivation pattern analysis to extract recurring whole-brain spatial patterns of brain activity over time.

RESULTS

The authors found 3 different sets of coactivating regions within the right extrastriate cortex in HCs and patients with pretherapeutic ET, reminiscent of the "cerebello-visuo-motor," "thalamo-visuo-motor" (including the targeted thalamus), and "basal ganglia and extrastriate" networks. The occurrence of the first pattern was decreased in pretherapeutic ET compared to HCs, whereas the other two patterns showed increased occurrences. This suggests a misbalance between the more prominent cerebellar circuitry and the thalamo-visuo-motor and basal ganglia networks. Multiple regression analysis showed that pretherapeutic standard tremor scores negatively correlated with the increased occurrence of the thalamo-visuo-motor network, suggesting a compensatory pathophysiological trait. Clinical improvement after SRS-T was related to changes in occurrences of the basal ganglia and extrastriate cortex circuitry, which returned to HC values after the intervention, suggesting that the dynamics of the extrastriate cortex had a role in tremor generation and further arrest after the intervention.

CONCLUSIONS

The data in this study point to a broader implication of the visual system in tremor generation, and not only through visual feedback, given its connections to the dorsal visual stream pathway and the cerebello-thalamo-cortical circuitry, with which its dynamic balance seems to be a crucial feature for reduced tremor. Furthermore, SRS-T seems to bring abnormal pretherapeutic connectivity of the extrastriate cortex to levels comparable to those of HC subjects.

Temporal sequences of brain activity at rest are constrained by white matter structure and modulated by cognitive demands

A diverse set of white matter connections supports seamless transitions between cognitive states. However, it remains unclear how these connections guide the temporal progression of large-scale brain activity patterns in different cognitive states. Here, we analyze the brain's trajectories across a set of single time point activity patterns from functional magnetic resonance imaging data acquired during the resting state and an n-back working memory task. We find that specific temporal sequences of brain activity are modulated by cognitive load, associated with age, and related to task performance. Using diffusion-weighted imaging acquired from the same subjects, we apply tools from network control theory to show that linear spread of activity along white matter connections constrains the probabilities of these sequences at rest, while stimulus-driven visual inputs explain the sequences observed during the n-back task. Overall, these results elucidate the structural underpinnings of cognitively and developmentally relevant spatiotemporal brain dynamics.

Characterization of network switching in disorder of consciousness at multiple time scales

OBJECTIVE

Recent works have shown that flexible information processing is closely related to the reconfiguration of human brain networks underlying brain functions. However, the role of network switching for consciousness is poorly explored and whether such transition can indicate the behavioral performance of patients with disorders of consciousness (DOC) remains unknown. Here, we investigate the relationship between the switching of brain networks (states) over time and the consciousness levels.

APPROACH

By applying multilayer network methods, we calculated time-resolved functional connectivity from source-level EEG data in different frequency bands. At various time scales, we explored how the human brain changes its community structure and traverses across defined network states (integrated and segregated states) in subjects with different consciousness levels.

MAIN RESULTS

Network switching in the human brain is decreased with increasing time scale opposite to that in random systems. Transitions of community assignment (denoted by flexibility) are negatively correlated with the consciousness levels (particularly in the alpha band) at short time scales. At long time scales, the opposite trend is found. Compared to healthy controls, patients show a new balance between dynamic segregation and integration, with decreased proportion and mean duration of segregated state (contrary to those of integrated state) at small scales.

SIGNIFICANCE

These findings may contribute to the development of EEG-based network analysis and shed new light on the pathological mechanisms of neurological disorders like DOC.

Comparison of feature selection methods based on discrimination and reliability for fMRI decoding analysis

BACKGROUND

Feature selection is a crucial step in the machine learning methods that are currently used to assist with decoding brain states from fMRI data. This step can be based on either feature discrimination or feature reliability, but there is no clear evidence indicating which method is more suitable for fMRI data.

METHODS

We used ANOVA and Kendall's concordance coefficient as proxies for the two kinds of feature selection criteria. The performances of both methods were compared using different subject and feature numbers. The study included 987 subjects from the Human Connectome Project (HCP).

RESULTS

Classification performance suggested that features based on discrimination were more capable of distinguishing between various brain states for any number of subjects or extracted features. In addition, reliability-based features were always more stable than other features, and these properties (discernment and stability) of features, to some degree, related to the number of subjects and features. Furthermore, when the number of extracted features increased, the feature distributions also gradually extended from occipital lobe to more association regions of the brain.

CONCLUSION

The results from this study provide empirical guides for feature selection for the prediction of individual brain states.

Time-resolved effective connectivity in task fMRI: Psychophysiological interactions of Co-Activation patterns

Investigating context-dependent modulations of Functional Connectivity (FC) with functional magnetic resonance imaging is crucial to reveal the neurological underpinnings of cognitive processing. Most current analysis methods hypothesise sustained FC within the duration of a task, but this assumption has been shown too limiting by recent imaging studies. While several methods have been proposed to study functional dynamics during rest, task-based studies are yet to fully disentangle network modulations. Here, we propose a seed-based method to probe task-dependent modulations of brain activity by revealing Psychophysiological Interactions of Co-activation Patterns (PPI-CAPs). This point process-based approach temporally decomposes task-modulated connectivity into dynamic building blocks which cannot be captured by current methods, such as PPI or Dynamic Causal Modelling. Additionally, it identifies the occurrence of co-activation patterns at single frame resolution as opposed to window-based methods. In a naturalistic setting where participants watched a TV program, we retrieved several patterns of co-activation with a posterior cingulate cortex seed whose occurrence rates and polarity varied depending on the context; on the seed activity; or on an interaction between the two. Moreover, our method exposed the consistency in effective connectivity patterns across subjects and time, allowing us to uncover links between PPI-CAPs and specific stimuli contained in the video. Our study reveals that explicitly tracking connectivity pattern transients is paramount to advance our understanding of how different brain areas dynamically communicate when presented with a set of cues.

Neural circuits of idiopathic Normal Pressure Hydrocephalus: A perspective review of brain connectivity and symptoms meta-analysis

Idiopathic normal pressure hydrocephalus (iNPH) is a prevalent reversible neurological disorder characterized by impaired locomotion, cognition and urinary control with ventriculomegaly. Symptoms can be relieved with cerebrospinal fluid drainage, which makes iNPH the leading cause of reversible dementia. Because of a limited understanding of pathophysiological mechanisms, unspecific symptoms and the high prevalence of comorbidity (i.e. Alzheimer's disease), iNPH is largely underdiagnosed. For these reasons, there is an urgent need for developing noninvasive quantitative biomarkers for iNPH diagnosis and prognosis. Structural and functional changes of brain circuits in relation to symptoms and treatment response are expected to deliver major advances in this direction. We review structural and functional brain connectivity findings in iNPH and complement those findings with iNPH symptom meta-analyses in healthy populations. Our goal is to reinforce our conceptualization of iNPH as to brain network mechanisms and foster the development of new hypotheses for future research and treatment options.

Triple Network Model Dynamically Revisited: Lower Salience Network State Switching in Pre-psychosis

Emerging evidence has attributed altered network coordination between the default mode, central executive, and salience networks (DMN/CEN/SAL) to disturbances seen in schizophrenia, but little is known for at-risk psychosis stages. Moreover, pinpointing impairments in specific network-to-network interactions, although essential to resolve possibly distinct harbingers of conversion to clinically diagnosed schizophrenia, remains particularly challenging. We addressed this by a dynamic approach to functional connectivity, where right anterior insula brain interactions were examined through co-activation pattern (CAP) analysis. We utilized resting-state fMRI in 19 subjects suffering from subthreshold delusions and hallucinations (UHR), 28 at-risk for psychosis with basic symptoms describing only self-experienced subclinical disturbances (BS), and 29 healthy controls (CTR) matched for age, gender, handedness, and intelligence. We extracted the most recurring CAPs, compared their relative occurrence and average dwell time to probe their temporal expression, and quantified occurrence balance to assess the putative loss of competing relationships. Our findings substantiate the pivotal role of the right anterior insula in governing CEN-to-DMN transitions, which appear dysfunctional prior to the onset of psychosis, especially when first attenuated psychotic symptoms occur. In UHR subjects, it is longer active in concert with the DMN and there is a loss of competition between a SAL/DMN state, and a state with insula/CEN activation paralleled by DMN deactivation. These features suggest that abnormal network switching disrupts one's capacity to distinguish between the internal world and external environment, which is accompanied by inflexibility and an excessive awareness to internal processes reflected by prolonged expression of the right anterior insula-default mode co-activation pattern.

A deconvolution algorithm for multi-echo functional MRI: Multi-echo Sparse Paradigm Free Mapping

This work introduces a novel algorithm for deconvolution of the BOLD signal in multi-echo fMRI data: Multi-echo Sparse Paradigm Free Mapping (ME-SPFM). Assuming a linear dependence of the BOLD percent signal change on the echo time (TE) and using sparsity-promoting regularized least squares estimation, ME-SPFM yields voxelwise time-varying estimates of the changes in the apparent transverse relaxation (ΔR₂^⁎) without prior knowledge of the timings of individual BOLD events. Our results in multi-echo fMRI data collected during a multi-task event-related paradigm at 3 Tesla demonstrate that the maps of R₂^⁎ changes obtained with ME-SPFM at the times of the stimulus trials show high spatial and temporal concordance with the activation maps and BOLD signals obtained with standard model-based analysis. This method yields estimates of ΔR₂^⁎ having physiologically plausible values. Owing to its ability to blindly detect events, ME-SPFM also enables us to map ΔR₂^⁎ associated with spontaneous, transient BOLD responses occurring between trials. This framework is a step towards deciphering the dynamic nature of brain activity in naturalistic paradigms, resting-state or experimental paradigms with unknown timing of the BOLD events.

The potential of MR-Encephalography for BCI/Neurofeedback applications with high temporal resolution

Real-time functional magnetic resonance imaging (rt-fMRI) enables the update of various brain-activity measures during an ongoing experiment as soon as a new brain volume is acquired. However, the recorded Blood-oxygen-level dependent (BOLD) signal also contains physiological artifacts such as breathing and heartbeat, which potentially cause misleading false positive effects especially problematic in brain-computer interface (BCI) and neurofeedback (NF) setups. The low temporal resolution of echo planar imaging (EPI) sequences (which is in the range of seconds) prevents a proper separation of these artifacts from the BOLD signal. MR-Encephalography (MREG) has been shown to provide the high temporal resolution required to unalias and correct for physiological fluctuations and leads to increased specificity and sensitivity for mapping task-based activation and functional connectivity as well as for detecting dynamic changes in connectivity over time. By comparing a simultaneous multislice echo planar imaging (SMS-EPI) sequence and an MREG sequence using the same nominal spatial resolution in an offline analysis for three different experimental fMRI paradigms (perception of house and face stimuli, motor imagery, Stroop task), the potential of this novel technique for future BCI and NF applications was investigated. First, adapted general linear model pre-whitening which accounts for the high temporal resolution in MREG was implemented to calculate proper statistical results and be able to compare these with the SMS-EPI sequence. Furthermore, the respiration- and cardiac pulsation-related signals were successfully separated from the MREG signal using independent component analysis which were then included as regressors for a GLM analysis. Only the MREG sequence allowed to clearly separate cardiac pulsation and respiration components from the signal time course. It could be shown that these components highly correlate with the recorded respiration and cardiac pulsation signals using a respiratory belt and fingertip pulse plethysmograph. Temporal signal-to-noise ratios of SMS-EPI and MREG were comparable. Functional connectivity analysis using partial correlation showed a reduced standard error in MREG compared to SMS-EPI. Also, direct time course comparisons by down-sampling the MREG signal to the SMS-EPI temporal resolution showed lower variance in MREG. In general, we show that the higher temporal resolution is beneficial for fMRI time course modeling and this aspect can be exploited in offline application but also, is especially attractive, for real-time BCI and NF applications.

Music in premature infants enhances high-level cognitive brain networks

Preterm babies are cared for in neonatal intensive care units (NICU), which are busy places with a lot of mechanical noise increasingly recognized to disrupt normal brain development. NICUs therefore invest in developmental care procedures, with music for example, but neurobiological evidence for these interventions is missing. We present results from a clinical trial to study the effects of a music intervention on preterm infants’ brain development. Based on resting-state fMRI, we provide evidence that music enhanced connectivity in a brain circuitry involving the salience network with regions implicated in sensory and higher-order cognitive functions, previously found to be altered in preterm infants. To our knowledge, this study is unique in observing an impact of music on brain development in preterm newborns.

Neonatal intensive care units are willing to apply environmental enrichment via music for preterm newborns. However, no evidence of an effect of music on preterm brain development has been reported to date. Using resting-state fMRI, we characterized a circuitry of interest consisting of three network modules interconnected by the salience network that displays reduced network coupling in preterm compared with full-term newborns. Interestingly, preterm infants exposed to music in the neonatal intensive care units have significantly increased coupling between brain networks previously shown to be decreased in premature infants: the salience network with the superior frontal, auditory, and sensorimotor networks, and the salience network with the thalamus and precuneus networks. Therefore, music exposure leads to functional brain architectures that are more similar to those of full-term newborns, providing evidence for a beneficial effect of music on the preterm brain.

Robust Recovery of Temporal Overlap Between Network Activity Using Transient-Informed Spatio-Temporal Regression

Functional magnetic resonance imaging is a non-invasive tomographic imaging modality that has provided insights into system-level brain function. New analysis methods are emerging to study the dynamic behavior of brain activity. The innovation-driven co-activation pattern (iCAP) approach is one such approach that relies on the detection of timepoints with a significant transient activity to subsequently retrieve spatially and temporally overlapping large-scale brain networks. To recover temporal profiles of the iCAPs for further time-resolved analysis, spatial patterns are fitted back to the activity-inducing signals. In this crucial step, spatial dependences can hinder the recovery of temporal overlapping activity. To overcome this effect, we propose a novel back-projection method that optimally fits activity-inducing signals given a set of transient timepoints and spatial maps of iCAPs, thus taking into account both spatial and temporal constraints. Validation on simulated data shows that transient-based constraints improve the quality of fitted time courses. Further evaluation on experimental data demonstrates that overfitting and underfitting are prevented by the use of optimized spatio-temporal constraints. Spatial and temporal properties of resulting iCAPs support that brain activity is characterized by the recurrent co-activation and co-deactivation of spatially overlapping large-scale brain networks. This new approach opens new avenues to explore the brain's dynamic core.

Brain imaging of locomotion in neurological conditions

Impaired locomotion is a frequent and major source of disability in patients with neurological conditions. Different neuroimaging methods have been used to understand the brain substrates of locomotion in various neurological diseases (mainly in Parkinson's disease) during actual walking, and while resting (using mental imagery of gait, or brain-behavior correlation analyses). These studies, using structural (i.e., MRI) or functional (i.e., functional MRI or functional near infra-red spectroscopy) brain imaging, electrophysiology (i.e., EEG), non-invasive brain stimulation (i.e., transcranial magnetic stimulation, or transcranial direct current stimulation) or molecular imaging methods (i.e., PET, or SPECT) reveal extended brain networks involving both grey and white matters in key cortical (i.e., prefrontal cortex) and subcortical (basal ganglia and cerebellum) regions associated with locomotion. However, the specific roles of the various pathophysiological mechanisms encountered in each neurological condition on the phenotype of gait disorders still remains unclear. After reviewing the results of individual brain imaging techniques across the common neurological conditions, such as Parkinson's disease, dementia, stroke, or multiple sclerosis, we will discuss how the development of new imaging techniques and computational analyses that integrate multivariate correlations in "large enough datasets" might help to understand how individual pathophysiological mechanisms express clinically as an abnormal gait. Finally, we will explore how these new analytic methods could drive our rehabilitative strategies.

Spatio-temporal dynamics of resting-state brain networks improve single-subject prediction of schizophrenia diagnosis

Correlation in functional MRI activity between spatially separated brain regions can fluctuate dynamically when an individual is at rest. These dynamics are typically characterized temporally by measuring fluctuations in functional connectivity between brain regions that remain fixed in space over time. Here, dynamics in functional connectivity were characterized in both time and space. Temporal dynamics were mapped with sliding-window correlation, while spatial dynamics were characterized by enabling network regions to vary in size (shrink/grow) over time according to the functional connectivity profile of their constituent voxels. These temporal and spatial dynamics were evaluated as biomarkers to distinguish schizophrenia patients from controls, and compared to current biomarkers based on static measures of resting-state functional connectivity. Support vector machine classifiers were trained using: (a) static, (b) dynamic in time, (c) dynamic in space, and (d) dynamic in time and space characterizations of functional connectivity within canonical resting-state brain networks. Classifiers trained on functional connectivity dynamics mapped over both space and time predicted diagnostic status with accuracy exceeding 91%, whereas utilizing only spatial or temporal dynamics alone yielded lower classification accuracies. Static measures of functional connectivity yielded the lowest accuracy (79.5%). Compared to healthy comparison individuals, schizophrenia patients generally exhibited functional connectivity that was reduced in strength and more variable. Robustness was established with replication in an independent dataset. The utility of biomarkers based on temporal and spatial functional connectivity dynamics suggests that resting-state dynamics are not trivially attributable to sampling variability and head motion.

Applications of community detection techniques to brain graphs: Algorithmic considerations and implications for neural function

The human brain can be represented as a graph in which neural units such as cells or small volumes of tissue are heterogeneously connected to one another through structural or functional links. Brain graphs are parsimonious representations of neural systems that have begun to offer fundamental insights into healthy human cognition, as well as its alteration in disease. A critical open question in network neuroscience lies in how neural units cluster into densely interconnected groups that can provide the coordinated activity that is characteristic of perception, action, and adaptive behaviors. Tools that have proven particularly useful for addressing this question are community detection approaches, which can identify communities or modules: groups of neural units that are densely interconnected with other units in their own group but sparsely interconnected with units in other groups. In this paper, we describe a common community detection algorithm known as modularity maximization, and we detail its applications to brain graphs constructed from neuroimaging data. We pay particular attention to important algorithmic considerations, especially in recent extensions of these techniques to graphs that evolve in time. After recounting a few fundamental insights that these techniques have provided into brain function, we highlight potential avenues of methodological advancements for future studies seeking to better characterize the patterns of coordinated activity in the brain that accompany human behavior. This tutorial provides a naive reader with an introduction to theoretical considerations pertinent to the generation of brain graphs, an understanding of modularity maximization for community detection, a resource of statistical measures that can be used to characterize community structure, and an appreciation of the usefulness of these approaches in uncovering behaviorally-relevant network dynamics in neuroimaging data.