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

Disuse-driven plasticity in the human thalamus and putamen

Motor adaptation in cortico-striato-thalamo-cortical loops has been studied mainly in animals using invasive electrophysiology. Here, we leverage functional neuroimaging in humans to study motor circuit plasticity in the human subcortex. We employed an experimental paradigm that combined two weeks of upper-extremity immobilization with daily resting-state and motor task fMRI before, during, and after the casting period. We previously showed that limb disuse leads to decreased functional connectivity (FC) of the contralateral somatomotor cortex (SM1) with the ipsilateral somatomotor cortex, increased FC with the cingulo-opercular network (CON) as well as the emergence of high amplitude, fMRI signal pulses localized in the contralateral SM1, supplementary motor area and the cerebellum. From our prior observations, it remains unclear whether the disuse plasticity affects the thalamus and striatum. We extended our analysis to include these subcortical regions and found that both exhibit strengthened cortical FC and spontaneous fMRI signal pulses induced by limb disuse. The dorsal posterior putamen and the central thalamus, mainly CM, VLP and VIM nuclei, showed disuse pulses and FC changes that lined up with fmri task activations from the Human connectome project motor system localizer, acquired before casting for each participant. Our findings provide a novel understanding of the role of the cortico-striato-thalamo-cortical loops in human motor plasticity and a potential link with the physiology of sleep regulation. Additionally, similarities with FC observation from Parkinson Disease (PD) questions a pathophysiological link with limb disuse.

Topological Data Analysis Captures Task-Driven fMRI Profiles in Individual Participants: A Classification Pipeline Based on Persistence

BOLD-based fMRI is the most widely used method for studying brain function. The BOLD signal while valuable, is beset with unique vulnerabilities. The most notable of these is the modest signal to noise ratio, and the relatively low temporal and spatial resolution. However, the high dimensional complexity of the BOLD signal also presents unique opportunities for functional discovery. Topological Data Analyses (TDA), a branch of mathematics optimized to search for specific classes of structure within high dimensional data may provide particularly valuable applications. In this investigation, we acquired fMRI data in the anterior cingulate cortex (ACC) using a basic motor control paradigm. Then, for each participant and each of three task conditions, fMRI signals in the ACC were summarized using two methods: a) TDA based methods of persistent homology and persistence landscapes and b) non-TDA based methods using a standard vectorization scheme. Finally, using machine learning (with support vector classifiers), classification accuracy of TDA and non-TDA vectorized data was tested across participants. In each participant, TDA-based classification out-performed the non-TDA based counterpart, suggesting that our TDA analytic pipeline better characterized task- and condition-induced structure in fMRI data in the ACC. Our results emphasize the value of TDA in characterizing task- and condition-induced structure in regional fMRI signals. In addition to providing our analytical tools for other users to emulate, we also discuss the unique role that TDA-based methods can play in the study of individual differences in the structure of functional brain signals in the healthy and the clinical brain.

A Survey on Brain Effective Connectivity Network Learning

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

Variability in Resting-State Functional Magnetic Resonance Imaging: The Effect of Body Mass, Blood Pressure, Hematocrit, and Glycated Hemoglobin on Hemodynamic and Neuronal Parameters

Introduction: Replicability has become an increasing focus within the scientific communities with the ongoing "replication crisis." One area that appears to struggle with unreliable results is resting-state functional magnetic resonance imaging (rs-fMRI). Therefore, the current study aimed at improving the knowledge of endogenous factors that contribute to inter-individual variability. Methods: Arterial blood pressure (BP), body mass, hematocrit, and glycated hemoglobin were investigated as potential sources of between-subject variability in rs-fMRI, in healthy individuals. Whether changes in resting-state networks (rs-networks) could be attributed to variability in the blood-oxygen-level-dependent (BOLD)-signal, changes in neuronal activity, or both was of special interest. Within-subject parameters were estimated by utilizing dynamic-causal modeling, as it allows to make inferences on the estimated hemodynamic (BOLD-signal dynamics) and neuronal parameters (effective connectivity) separately. Results: The results of the analyses imply that BP and body mass can cause between-subject and between-group variability in the BOLD-signal and that all the included factors can affect the underlying connectivity. Discussion: Given the results of the current and previous studies, rs-fMRI results appear to be susceptible to a range of factors, which is likely to contribute to the low degree of replicability of these studies. Interestingly, the highest degree of variability seems to appear within the much-studied default mode network and its connections to other networks. Impact statement We believe that thanks to the evidence that we have collected by analyzing the well-controlled data of the Human Connectome Project with dynamic-causal modeling (DCM) and by focusing not only on the effective connectivity, which is the typical way of using DCM, but also by analyzing the underlying hemodynamic parameters, we were able to explore the underlying vascular dependencies in a much broader perspective. Our results challenge the premise for studying changes in the default mode network as a clinical marker of disease, and we add to the growing list of factors that contribute to resting-state network variability.

Improved localization of language areas using single voxel signal analysis of unprocessed fMRI data

Activated brain regions can be visualized and localized with the use of fMRI (functional magnetic imaging). This is based on changes in the blood flow in activated regions, or more precisely on the hemodynamic response function (HRF) and the Blood-Oxygen-Level-Dependent (BOLD) effect. This study used a task-based fMRI examination with language paradigms in order to stimulate the language areas. The measured fMRI data are frequently altered by different preprocessing steps for the analysis and the display of activations. These changes can lead to discrepancies between the displayed and the truly measured location of the activations. Simple t-maps were created with unprocessed fMRI data, to provide a more realistic representation of the language areas. HRF-dependent single-voxel fMRI signal analysis was performed to improve the analyzability of these activation maps.

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.

The anterior insula channels prefrontal expectancy signals during affective processing

Expectancy shapes our perception of impending events. Although such an interplay between cognitive and affective processes is often impaired in mental disorders, it is not well understood how top-down expectancy signals modulate future affect. We therefore track the information flow in the brain during cognitive and affective processing segregated in time using task-specific cross-correlations. Participants in two independent fMRI studies (N₁ = 37 & N₂ = 55) were instructed to imagine a situation with affective content as indicated by a cue, which was then followed by an emotional picture congruent with expectancy. To correct for intrinsic covariance of brain function, we calculate resting-state cross-correlations analogous to the task. First, using factorial modeling of delta cross-correlations (task-rest) of the first study, we find that the magnitude of expectancy signals in the anterior insula cortex (AIC) modulates the BOLD response to emotional pictures in the anterior cingulate and dorsomedial prefrontal cortex in opposite directions. Second, using hierarchical linear modeling of lagged connectivity, we demonstrate that expectancy signals in the AIC indeed foreshadow this opposing pattern in the prefrontal cortex. Third, we replicate the results in the second study using a higher temporal resolution, showing that our task-specific cross-correlation approach robustly uncovers the dynamics of information flow. We conclude that the AIC arbitrates the recruitment of distinct prefrontal networks during cued picture processing according to triggered expectations. Taken together, our study provides new insights into neuronal pathways channeling cognition and affect within well-defined brain networks. Better understanding of such dynamics could lead to new applications tracking aberrant information processing in mental disorders.

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.

Disentangling Multispectral Functional Connectivity With Wavelets

The field of brain connectomics develops our understanding of the brain's intrinsic organization by characterizing trends in spontaneous brain activity. Linear correlations in spontaneous blood-oxygen level dependent functional magnetic resonance imaging (BOLD-fMRI) fluctuations are often used as measures of functional connectivity (FC), that is, as a quantity describing how similarly two brain regions behave over time. Given the natural spectral scaling of BOLD-fMRI signals, it may be useful to represent BOLD-fMRI as multiple processes occurring over multiple scales. The wavelet domain presents a transform space well suited to the examination of multiscale systems as the wavelet basis set is constructed from a self-similar rescaling of a time and frequency delimited kernel. In the present study, we utilize wavelet transforms to examine fluctuations in whole-brain BOLD-fMRI connectivity as a function of wavelet spectral scale in a sample (N = 31) of resting healthy human volunteers. Information theoretic criteria measure relatedness between spectrally-delimited FC graphs. Voxelwise comparisons of between-spectra graph structures illustrate the development of preferential functional networks across spectral bands.

Thresholding functional connectomes by means of mixture modeling

Functional connectivity has been shown to be a very promising tool for studying the large-scale functional architecture of the human brain. In network research in fMRI, functional connectivity is considered as a set of pair-wise interactions between the nodes of the network. These interactions are typically operationalized through the full or partial correlation between all pairs of regional time series. Estimating the structure of the latent underlying functional connectome from the set of pair-wise partial correlations remains an open research problem though. Typically, this thresholding problem is approached by proportional thresholding, or by means of parametric or non-parametric permutation testing across a cohort of subjects at each possible connection. As an alternative, we propose a data-driven thresholding approach for network matrices on the basis of mixture modeling. This approach allows for creating subject-specific sparse connectomes by modeling the full set of partial correlations as a mixture of low correlation values associated with weak or unreliable edges in the connectome and a sparse set of reliable connections. Consequently, we propose to use alternative thresholding strategy based on the model fit using pseudo-False Discovery Rates derived on the basis of the empirical null estimated as part of the mixture distribution. We evaluate the method on synthetic benchmark fMRI datasets where the underlying network structure is known, and demonstrate that it gives improved performance with respect to the alternative methods for thresholding connectomes, given the canonical thresholding levels. We also demonstrate that mixture modeling gives highly reproducible results when applied to the functional connectomes of the visual system derived from the n-back Working Memory task in the Human Connectome Project. The sparse connectomes obtained from mixture modeling are further discussed in the light of the previous knowledge of the functional architecture of the visual system in humans. We also demonstrate that with use of our method, we are able to extract similar information on the group level as can be achieved with permutation testing even though these two methods are not equivalent. We demonstrate that with both of these methods, we obtain functional decoupling between the two hemispheres in the higher order areas of the visual cortex during visual stimulation as compared to the resting state, which is in line with previous studies suggesting lateralization in the visual processing. However, as opposed to permutation testing, our approach does not require inference at the cohort level and can be used for creating sparse connectomes at the level of a single subject.

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.