The relative brain signal variability increases in the behavioral variant of frontotemporal dementia and Alzheimer's disease but not in schizophrenia

Overlapping symptoms between Alzheimer's disease (AD), behavioral variant of frontotemporal dementia (bvFTD), and schizophrenia (SZ) can lead to misdiagnosis and delays in appropriate treatment, especially in cases of early-onset dementia. To determine the potential of brain signal variability as a diagnostic tool, we assessed the coefficient of variation of the BOLD signal (CVBOLD) in 234 participants spanning bvFTD (n = 53), AD (n = 17), SZ (n = 23), and controls (n = 141). All underwent functional and structural MRI scans. Data unveiled a notable increase in CVBOLD in bvFTD patients across both datasets (local and international, p < 0.05), revealing an association with clinical scores (CDR and MMSE, r = 0.46 and r = -0.48, p < 0.0001). While SZ and control group demonstrated no significant differences, a comparative analysis between AD and bvFTD patients spotlighted elevated CVBOLD in the frontopolar cortices for the latter (p < 0.05). Furthermore, CVBOLD not only presented excellent diagnostic accuracy for bvFTD (AUC 0.78-0.95) but also showcased longitudinal repeatability. During a one-year follow-up, the CVBOLD levels increased by an average of 35% in the bvFTD group, compared to a 2% increase in the control group (p < 0.05). Our findings suggest that CVBOLD holds promise as a biomarker for bvFTD, offering potential for monitoring disease progression and differentiating bvFTD from AD and SZ.

Resting state neurophysiology of agonist-antagonist myoneural interface in persons with transtibial amputation

The agonist-antagonist myoneural interface (AMI) is an amputation surgery that preserves sensorimotor signaling mechanisms of the central-peripheral nervous systems. Our first neuroimaging study investigating AMI subjects conducted by Srinivasan et al. (2020) focused on task-based neural signatures, and showed evidence of proprioceptive feedback to the central nervous system. The study of resting state neural activity helps non-invasively characterize the neural patterns that prime task response. In this study on resting state functional magnetic resonance imaging in AMI subjects, we compared functional connectivity in patients with transtibial AMI (n = 12) and traditional (n = 7) amputations (TA). To test our hypothesis that we would find significant neurophysiological differences between AMI and TA subjects, we performed a whole-brain exploratory analysis to identify a seed region; namely, we conducted ANOVA, followed by t-test statistics to locate a seed in the salience network. Then, we implemented a seed-based connectivity analysis to gather cluster-level inferences contrasting our subject groups. We show evidence supporting our hypothesis that the AMI surgery induces functional network reorganization resulting in a neural configuration that significantly differs from the neural configuration after TA surgery. AMI subjects show significantly less coupling with regions functionally dedicated to selecting where to focus attention when it comes to salient stimuli. Our findings provide researchers and clinicians with a critical mechanistic understanding of the effect of AMI amputation on brain networks at rest, which has promising implications for improved neurorehabilitation and prosthetic control.

Hemodynamic response function description in patients with glioma

INTRODUCTION

Functional magnetic resonance imaging is a powerful tool that has provided many insights into cognitive sciences. Yet, as its analysis is mostly based on the knowledge of an a priori canonical hemodynamic response function (HRF), its reliability in patients' applications has been questioned. There have been reports of neurovascular uncoupling in patients with glioma, but no specific description of the Hemodynamic Response Function (HRF) in glioma has been reported so far. The aim of this work is to describe the HRF in patients with glioma.

METHODS

Forty patients were included. MR images were acquired on a 1.5T scanner. Activated clusters were identified using a fuzzy general linear model; HRFs were adjusted with a double-gamma function. Analyses were undertaken considering the tumor grade, age, sex, tumor location, and activated location.

RESULTS

Differences are found in the occipital, limbic, insular, and sub-lobar areas, but not in the frontal, temporal, and parietal lobes. The presence of a glioma slows the time-to-peak and onset times by 5.2 and 3.8 % respectively; high-grade gliomas present 8.1 % smaller HRF widths than low-grade gliomas.

DISCUSSION AND CONCLUSION

There is significant HRF variation due to the presence of glioma, but the magnitudes of the observed differences are small. Most processing pipelines should be robust enough for this magnitude of variation and little if any impact should be visible on functional maps. The differences that have been observed in the literature between functional mapping obtained with magnetic resonance vs. that obtained with direct electrostimulation during awake surgery are more probably due to the intrinsic difference in the mapping process: fMRI mapping detects all recruited areas while intra-surgical mapping indicates only the areas indispensable for the realization of a certain task. Surgical mapping might not be the gold standard to use when trying to validate the fMRI mapping process.

Blink-related arousal network surges are shaped by cortical vigilance states

The vigilance state and the excitability of cortical networks impose wide-range effects on brain dynamics that arousal surges could promptly modify. We previously reported an association between spontaneous eye-blinks and BOLD activation in the brain arousal ascending network (AAN) and in thalamic nuclei based on 3T MR resting state brain images. Here we aimed to replicate our analyses using 7T MR images in a larger cohort of participants collected from the Human Connectome Project (HCP), which also contained simultaneous eye-tracking recordings, and to assess the interaction between the blink-associated arousal surges and the vigilance states. For this purpose, we compared blink associated BOLD activity under a vigilant versus a drowsy state, a classification made based on the pupillary data obtained during the fMRI scans. We conducted two main analyses: i) Cross-correlation analysis between the BOLD signal and blink events (eye blink time-series were convolved with the canonical and also with the temporal derivative of the Hemodynamic Response Function, HRF) within preselected regions of interests (ROIs) (i.e., brainstem AAN, thalamic and cerebellar nuclei) together with an exploratory voxel-wise analyses to assess the whole-brain, and ii) blink-event analysis of the BOLD signals to reveal the signal changes onset to the blinks in the preselected ROIs. Consistent with our prior findings on 3T MRI, we showed significant positive cross correlations between BOLD peaks in brainstem and thalamic nuclei that preceded or were overlapping with blink moments and that sharply decreased post-blink. Whole brain analysis revealed blink-related activation that was strongest in cerebellum, insula, lateral geniculate nucleus (LGN) and visual cortex. Drowsiness impacted HRF BOLD (enhancing it), time-to-peak (delaying it) and post-blink BOLD activity (accentuating decreases). Responses in the drowsy state could be related to the differences in the excitability of cortical, subcortical and cerebellar tissue, such that cerebellar and thalamic regions involved in visual attention processing were more responsive for the vigilant state, but AAN ROIs, as well as cerebellar and thalamic ROIs connected to pre-motor, frontal, temporal and DMN regions were less responsive. Such qualitative and quantitative differences in the blink related BOLD signal changes could reflect delayed cortical processing and the ineffectiveness of arousal surges during states of drowsiness. Future studies that manipulate arousal are needed to corroborate a mechanistic interaction of arousal surges with vigilance states and cortical excitability.

Spectral dynamic causal modeling: A didactic introduction and its relationship with functional connectivity

We present a didactic introduction to spectral dynamic causal modeling (DCM), a Bayesian state-space modeling approach used to infer effective connectivity from noninvasive neuroimaging data. Spectral DCM is currently the most widely applied DCM variant for resting-state functional MRI analysis. Our aim is to explain its technical foundations to an audience with limited expertise in state-space modeling and spectral data analysis. Particular attention will be paid to cross-spectral density, which is the most distinctive feature of spectral DCM and is closely related to functional connectivity, as measured by (zero-lag) Pearson correlations. In fact, the model parameters estimated by spectral DCM are those that best reproduce the cross-correlations between all measurements-at all time lags-including the zero-lag correlations that are usually interpreted as functional connectivity. We derive the functional connectivity matrix from the model equations and show how changing a single effective connectivity parameter can affect all pairwise correlations. To complicate matters, the pairs of brain regions showing the largest changes in functional connectivity do not necessarily coincide with those presenting the largest changes in effective connectivity. We discuss the implications and conclude with a comprehensive summary of the assumptions and limitations of spectral DCM.

Traces of EEG-fMRI coupling reveals neurovascular dynamics on sleep inertia

Upon emergence from sleep, individuals experience temporary hypo-vigilance and grogginess known as sleep inertia. During the transient period of vigilance recovery from prior nocturnal sleep, the neurovascular coupling (NVC) may not be static and constant as assumed by previous neuroimaging studies. Stemming from this viewpoint of sleep inertia, this study aims to probe the NVC changes as awakening time prolongs using simultaneous EEG-fMRI. The time-lagged coupling between EEG features of vigilance and BOLD-fMRI signals, in selected regions of interest, was calculated with one pre-sleep and three consecutive post-awakening resting-state measures. We found marginal changes in EEG theta/beta ratio and spectral slope across post-awakening sessions, demonstrating alterations of vigilance during sleep inertia. Time-varying EEG-fMRI coupling as awakening prolonged was evidenced by the changing time lags of the peak correlation between EEG alpha-vigilance and fMRI-thalamus, as well as EEG spectral slope and fMRI-anterior cingulate cortex. This study provides the first evidence of potential dynamicity of NVC occurred in sleep inertia and opens new avenues for non-invasive neuroimaging investigations into the neurophysiological mechanisms underlying brain state transitions.

Macroscopic resting-state brain dynamics are best described by linear models

It is typically assumed that large networks of neurons exhibit a large repertoire of nonlinear behaviours. Here we challenge this assumption by leveraging mathematical models derived from measurements of local field potentials via intracranial electroencephalography and of whole-brain blood-oxygen-level-dependent brain activity via functional magnetic resonance imaging. We used state-of-the-art linear and nonlinear families of models to describe spontaneous resting-state activity of 700 participants in the Human Connectome Project and 122 participants in the Restoring Active Memory project. We found that linear autoregressive models provide the best fit across both data types and three performance metrics: predictive power, computational complexity and the extent of the residual dynamics unexplained by the model. To explain this observation, we show that microscopic nonlinear dynamics can be counteracted or masked by four factors associated with macroscopic dynamics: averaging over space and over time, which are inherent to aggregated macroscopic brain activity, and observation noise and limited data samples, which stem from technological limitations. We therefore argue that easier-to-interpret linear models can faithfully describe macroscopic brain dynamics during resting-state conditions.

Static and dynamic fMRI-derived functional connectomes represent largely similar information

Functional connectivity (FC) of blood oxygen level-dependent (BOLD) fMRI time series can be estimated using methods that differ in sensitivity to the temporal order of time points (static vs. dynamic) and the number of regions considered in estimating a single edge (bivariate vs. multivariate). Previous research suggests that dynamic FC explains variability in FC fluctuations and behavior beyond static FC. Our aim was to systematically compare methods on both dimensions. We compared five FC methods: Pearson's/full correlation (static, bivariate), lagged correlation (dynamic, bivariate), partial correlation (static, multivariate), and multivariate AR model with and without self-connections (dynamic, multivariate). We compared these methods by (i) assessing similarities between FC matrices, (ii) by comparing node centrality measures, and (iii) by comparing the patterns of brain-behavior associations. Although FC estimates did not differ as a function of sensitivity to temporal order, we observed differences between the multivariate and bivariate FC methods. The dynamic FC estimates were highly correlated with the static FC estimates, especially when comparing group-level FC matrices. Similarly, there were high correlations between the patterns of brain-behavior associations obtained using the dynamic and static FC methods. We conclude that the dynamic FC estimates represent information largely similar to that of the static FC.

The virtual aging brain: Causal inference supports interhemispheric dedifferentiation in healthy aging

The mechanisms of cognitive decline and its variability during healthy aging are not fully understood, but have been associated with reorganization of white matter tracts and functional brain networks. Here, we built a brain network modeling framework to infer the causal link between structural connectivity and functional architecture and the consequent cognitive decline in aging. By applying in-silico interhemispheric degradation of structural connectivity, we reproduced the process of functional dedifferentiation during aging. Thereby, we found the global modulation of brain dynamics by structural connectivity to increase with age, which was steeper in older adults with poor cognitive performance. We validated our causal hypothesis via a deep-learning Bayesian approach. Our results might be the first mechanistic demonstration of dedifferentiation during aging leading to cognitive decline.

Comparison of hemodynamic response functions obtained from resting-state functional MRI and invasive electrophysiological recordings in rats

Resting-state functional MRI (rs-fMRI) is a popular technology that has enriched our understanding of brain and spinal cord functioning, including how different regions communicate (connectivity). But fMRI is an indirect measure of neural activity capturing blood hemodynamics. The hemodynamic response function (HRF) interfaces between the unmeasured neural activity and measured fMRI time series. The HRF is variable across brain regions and individuals, and is modulated by non-neural factors. Ignoring this HRF variability causes errors in FC estimates. Hence, it is crucial to reliably estimate the HRF from rs-fMRI data. Robust techniques have emerged to estimate the HRF from fMRI time series. Although such techniques have been validated non-invasively using simulated and empirical fMRI data, thorough invasive validation using simultaneous electrophysiological recordings, the gold standard, has been elusive. This report addresses this gap in the literature by comparing HRFs derived from invasive intracranial electroencephalogram recordings with HRFs estimated from simultaneously acquired fMRI data in six epileptic rats. We found that the HRF shape parameters (HRF amplitude, latency and width) were not significantly different (p>0.05) between ground truth and estimated HRFs. In the single pathological region, the HRF width was marginally significantly different (p=0.03). Our study provides preliminary invasive validation for the efficacy of the HRF estimation technique in reliably estimating the HRF non-invasively from rs-fMRI data directly. This has a notable impact on rs-fMRI connectivity studies, and we recommend that HRF deconvolution be performed to minimize HRF variability and improve connectivity estimates.

Association of Intrinsic Functional Connectivity between the Locus Coeruleus and Salience Network with Attentional Ability

The locus coeruleus (LC) is a brainstem region associated with broad neural arousal because of norepinephrine production, but it has increasingly been associated with specific cognitive processes. These include sustained attention, with deficits associated with various neuropsychological disorders. Neural models of attention deficits have focused on interrupted dynamics between the salience network (SAL) with the frontoparietal network, which has been associated with task-switching and processing of external stimuli, respectively. Conflicting findings for these regions suggest the possibility of upstream signaling leading to attention dysfunction, and recent research suggests LC involvement. In this study, resting-state functional connectivity and behavioral performance on an attention task was examined within 584 individuals. Analysis revealed significant clusters connected to LC activity in the SAL. Given previous findings that attention deficits may be caused by SAL network switching dysfunctions, findings here further suggest that dysfunction in LC-SAL connectivity may impair attention.

Aging affects the phase coherence between spontaneous oscillations in brain oxygenation and neural activity

The risk of neurodegenerative disorders increases with age, due to reduced vascular nutrition and impaired neural function. However, the interactions between cardiovascular dynamics and neural activity, and how these interactions evolve in healthy aging, are not well understood. Here, the interactions are studied by assessment of the phase coherence between spontaneous oscillations in cerebral oxygenation measured by fNIRS, the electrical activity of the brain measured by EEG, and cardiovascular functions extracted from ECG and respiration effort, all simultaneously recorded. Signals measured at rest in 21 younger participants (31.1 ± 6.9 years) and 24 older participants (64.9 ± 6.9 years) were analysed by wavelet transform, wavelet phase coherence and ridge extraction for frequencies between 0.007 and 4 Hz. Coherence between the neural and oxygenation oscillations at ∼ 0.1 Hz is significantly reduced in the older adults in 46/176 fNIRS-EEG probe combinations. This reduction in coherence cannot be accounted for in terms of reduced power, thus indicating that neurovascular interactions change with age. The approach presented promises a noninvasive means of evaluating the efficiency of the neurovascular unit in aging and disease.

Resting-state fMRI signals contain spectral signatures of local hemodynamic response timing

Functional magnetic resonance imaging (fMRI) has proven to be a powerful tool for noninvasively measuring human brain activity; yet, thus far, fMRI has been relatively limited in its temporal resolution. A key challenge is understanding the relationship between neural activity and the blood-oxygenation-level-dependent (BOLD) signal obtained from fMRI, generally modeled by the hemodynamic response function (HRF). The timing of the HRF varies across the brain and individuals, confounding our ability to make inferences about the timing of the underlying neural processes. Here, we show that resting-state fMRI signals contain information about HRF temporal dynamics that can be leveraged to understand and characterize variations in HRF timing across both cortical and subcortical regions. We found that the frequency spectrum of resting-state fMRI signals significantly differs between voxels with fast versus slow HRFs in human visual cortex. These spectral differences extended to subcortex as well, revealing significantly faster hemodynamic timing in the lateral geniculate nucleus of the thalamus. Ultimately, our results demonstrate that the temporal properties of the HRF impact the spectral content of resting-state fMRI signals and enable voxel-wise characterization of relative hemodynamic response timing. Furthermore, our results show that caution should be used in studies of resting-state fMRI spectral properties, because differences in fMRI frequency content can arise from purely vascular origins. This finding provides new insight into the temporal properties of fMRI signals across voxels, which is crucial for accurate fMRI analyses, and enhances the ability of fast fMRI to identify and track fast neural dynamics.

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.

Pupil size reflects activation of subcortical ascending arousal system nuclei during rest

Neuromodulatory nuclei that are part of the ascending arousal system (AAS) play a crucial role in regulating cortical state and optimizing task performance. Pupil diameter, under constant luminance conditions, is increasingly used as an index of activity of these AAS nuclei. Indeed, task-based functional imaging studies in humans have begun to provide evidence of stimulus-driven pupil-AAS coupling. However, whether there is such a tight pupil-AAS coupling during rest is not clear. To address this question, we examined simultaneously acquired resting-state fMRI and pupil-size data from 74 participants, focusing on six AAS nuclei: the locus coeruleus, ventral tegmental area, substantia nigra, dorsal and median raphe nuclei, and cholinergic basal forebrain. Activation in all six AAS nuclei was optimally correlated with pupil size at 0-2 s lags, suggesting that spontaneous pupil changes were almost immediately followed by corresponding BOLD-signal changes in the AAS. These results suggest that spontaneous changes in pupil size that occur during states of rest can be used as a noninvasive general index of activity in AAS nuclei. Importantly, the nature of pupil-AAS coupling during rest appears to be vastly different from the relatively slow canonical hemodynamic response function that has been used to characterize task-related pupil-AAS coupling.

Static and dynamic functional connectomes represent largely similar information

Functional connectivity (FC) of blood-oxygen-level-dependent (BOLD) fMRI time series can be estimated using methods that differ in sensitivity to the temporal order of time points (static vs. dynamic) and the number of regions considered in estimating a single edge (bivariate vs. multivariate). Previous research suggests that dynamic FC explains variability in FC fluctuations and behavior beyond static FC. Our aim was to systematically compare methods on both dimensions. We compared five FC methods: Pearson's/full correlation (static, bivariate), lagged correlation (dynamic, bivariate), partial correlation (static, multivariate) and multivariate AR model with and without self-connections (dynamic, multivariate). We compared these methods by (i) assessing similarities between FC matrices, (ii) by comparing node centrality measures, and (iii) by comparing the patterns of brain-behavior associations. Although FC estimates did not differ as a function of sensitivity to temporal order, we observed differences between the multivariate and bivariate FC methods. The dynamic FC estimates were highly correlated with the static FC estimates, especially when comparing group-level FC matrices. Similarly, there were high correlations between the patterns of brain-behavior associations obtained using the dynamic and static FC methods. We conclude that the dynamic FC estimates represent information largely similar to that of the static FC.

Heritability of human "directed" functional connectome

INTRODUCTION

The functional connectivity patterns in the brain are highly heritable; however, it is unclear how genetic factors influence the directionality of such "information flows." Studying the "directionality" of the brain functional connectivity and assessing how heritability modulates it can improve our understanding of the human connectome.

METHODS

Here, we investigated the heritability of "directed" functional connections using a state-space formulation of Granger causality (GC), in conjunction with blind deconvolution methods accounting for local variability in the hemodynamic response function. Such GC implementation is ideal to explore the directionality of functional interactions across a large number of networks. Resting-state functional magnetic resonance imaging data were drawn from the Human Connectome Project (total n = 898 participants). To add robustness to our findings, the dataset was randomly split into a "discovery" and a "replication" sample (each with n = 449 participants). The two cohorts were carefully matched in terms of demographic variables and other confounding factors (e.g., education). The effect of shared environment was also modeled.

RESULTS

The parieto- and prefronto-cerebellar, parieto-prefrontal, and posterior-cingulate to hippocampus connections showed the highest and most replicable heritability effects with little influence by shared environment. In contrast, shared environmental factors significantly affected the visuo-parietal and sensory-motor directed connectivity.

CONCLUSION

We suggest a robust role of heritability in influencing the directed connectivity of some cortico-subcortical circuits implicated in cognition. Further studies, for example using task-based fMRI and GC, are warranted to confirm the asymmetric effects of genetic factors on the functional connectivity within cognitive networks and their role in supporting executive functions and learning.

Resting state neurophysiology of agonist-antagonist myoneural interface in persons with transtibial amputation

The agonist-antagonist myoneural interface (AMI) is a novel amputation surgery that preserves sensorimotor signaling mechanisms of the central-peripheral nervous systems. Our first neuroimaging study investigating AMI subjects (Srinivasan et al., Sci. Transl. Med. 2020) focused on task-based neural signatures, and showed evidence of proprioceptive feedback to the central nervous system. The study of resting state neural activity helps non-invasively characterize the neural patterns that prime task response. In this first study on resting state fMRI in AMI subjects, we compared resting state functional connectivity in patients with transtibial AMI (n=12) and traditional (n=7) amputations, as well as biologically intact control subjects (n=10). We hypothesized that the AMI surgery will induce functional network reorganization that significantly differs from the traditional amputation surgery and also more closely resembles the neural configuration of controls. We found AMI subjects to have lower connectivity with salience and motor seed regions compared to traditional amputees. Additionally, with connections affected in traditional amputees, AMI subjects exhibited a connectivity pattern more closely resembling controls. Lastly, sensorimotor connectivity in amputee cohorts was significantly associated with phantom sensation (R2=0.7, p=0.0008). These findings provide researchers and clinicians with a critical mechanistic understanding of the effects of the AMI surgery on the brain at rest, spearheading future research towards improved prosthetic control and embodiment.

Resting-state fMRI signals contain spectral signatures of local hemodynamic response timing

Functional magnetic resonance imaging (fMRI) has proven to be a powerful tool for noninvasively measuring human brain activity; yet, thus far, fMRI has been relatively limited in its temporal resolution. A key challenge is understanding the relationship between neural activity and the blood-oxygenation-level-dependent (BOLD) signal obtained from fMRI, generally modeled by the hemodynamic response function (HRF). The timing of the HRF varies across the brain and individuals, confounding our ability to make inferences about the timing of the underlying neural processes. Here we show that resting-state fMRI signals contain information about HRF temporal dynamics that can be leveraged to understand and characterize variations in HRF timing across both cortical and subcortical regions. We found that the frequency spectrum of resting-state fMRI signals significantly differs between voxels with fast versus slow HRFs in human visual cortex. These spectral differences extended to subcortex as well, revealing significantly faster hemodynamic timing in the lateral geniculate nucleus of the thalamus. Ultimately, our results demonstrate that the temporal properties of the HRF impact the spectral content of resting-state fMRI signals and enable voxel-wise characterization of relative hemodynamic response timing. Furthermore, our results show that caution should be used in studies of resting-state fMRI spectral properties, as differences can arise from purely vascular origins. This finding provides new insight into the temporal properties of fMRI signals across voxels, which is crucial for accurate fMRI analyses, and enhances the ability of fast fMRI to identify and track fast neural dynamics.

Investigating dynamic causal network with unified Granger causality analysis

BACKGROUND

Dynamic coupling phenomena characterize a widespread fundamental mechanism for the functional brain, which involves large-scale interactions at a multi-level. The Granger causality analysis (GCA) provides a data-driven procedure to investigate causal connections and has the potential to be a powerful dynamic capturing tool.

NEW METHOD

In this paper, distinct from the conventional two-stage scheme of most GCA methods, we suggest a unified GCA (uGCA) method incorporating a sliding window to further capture dynamic connections. And the uGCA method integrates all related procedures into the same space by a single mathematical theory, which involves a description length guided framework.

RESULTS

Through synthetic data experiments and real fMRI data experiments, we illustrated the effectiveness and priority of the proposed uGCA method.

COMPARISON WITH EXISTING METHODS

By varying the data length, we have demonstrated its superiority to conventional GCA in synthetic data experiments. We further illustrated the outstanding capability of their dynamic causal investigation in the fMRI data, involving serial mental arithmetic tasks under visual and auditory stimuli, respectively, one can evaluate the performance of different methods by accessing their network similarities among different stimuli. When varying windows size and step size of the sliding window, respectively, compared with conventional GCA, the uGCA identified higher network similarities while ensuring more robust performance.

CONCLUSIONS

The stability and effectiveness of uGCA will show it an advantage in the further research of multi-level dynamic coupling and characterizing.

autohrf-an R package for generating data-informed event models for general linear modeling of task-based fMRI data

The analysis of task-related fMRI data at the level of individual participants is commonly based on general linear modeling (GLM), which allows us to estimate the extent to which the BOLD signal can be explained by the task response predictors specified in the event model. The predictors are constructed by convolving the hypothesized time course of neural activity with an assumed hemodynamic response function (HRF). However, our assumptions about the components of brain activity, including their onset and duration, may be incorrect. Their timing may also differ across brain regions or from person to person, leading to inappropriate or suboptimal models, poor fit of the model to actual data, and invalid estimates of brain activity. Here, we present an approach that uses theoretically driven models of task response to define constraints on which the final model is computationally derived using actual fMRI data. Specifically, we developed autohrf-an R package that enables the evaluation and data-driven estimation of event models for GLM analysis. The highlight of the package is the automated parameter search that uses genetic algorithms to find the onset and duration of task predictors that result in the highest fitness of GLM based on the fMRI signal under predefined constraints. We evaluated the usefulness of the autohrf package on two original datasets of task-related fMRI activity, a slow event-related spatial working memory study and a mixed state-item study using the flanker task, and on a simulated slow event-related working memory data. Our results suggest that autohrf can be used to efficiently construct and evaluate better task-related brain activity models to gain a deeper understanding of BOLD task response and improve the validity of model estimates. Our study also highlights the sensitivity of fMRI analysis with GLM to precise event model specification and the need for model evaluation, especially in complex and overlapping event designs.

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.

Pre- and post-surgery brain tumor multimodal magnetic resonance imaging data optimized for large scale computational modelling

We present a dataset of magnetic resonance imaging (MRI) data (T1, diffusion, BOLD) acquired in 25 brain tumor patients before the tumor resection surgery, and six months after the surgery, together with the tumor masks, and in 11 controls (recruited among the patients’ caregivers). The dataset also contains behavioral and emotional scores obtained with standardized questionnaires. To simulate personalized computational models of the brain, we also provide structural connectivity matrices, necessary to perform whole-brain modelling with tools such as The Virtual Brain. In addition, we provide blood-oxygen-level-dependent imaging time series averaged across regions of interest for comparison with simulation results. An average resting state hemodynamic response function for each region of interest, as well as shape maps for each voxel, are also contributed.

Measurement(s)	BOLD signal • Diffusion Anisotropy
Technology Type(s)	Functional Magnetic Resonance Imaging • Diffusion Weighted Imaging
Factor Type(s)	Surgery
Sample Characteristic - Organism	Homo sapiens

Visualization of the Dynamic Brain Activation Pattern during a Decision-Making Task

Decision making is a complex process involving various parts of the brain which are active during different times. It is challenging to measure externally the exact instant when any given region becomes active during the decision-making process. Here, we propose the development and validation of an algorithm to extract and visualize the dynamic functional brain activation information from the observed fMRI data. We propose the use of a regularized deconvolution model to simultaneously map various activation regions within the brain and track how different activation regions changes with time, thus providing both spatial and temporal brain activation information. The proposed technique was validated using simulated data and then applied to a simple decision-making task for identification of various brain regions involved in different stages of decision making. Using the results of the dynamic activation for the decision-making task, we were able to identify key brain regions involved in some of the phases of decision making. The visualization aspect of the algorithm allows us to actually see the flow of activation (and deactivation) in the form of a motion picture. The dynamic estimate may aid in understanding the causality of activation between various brain regions in a better way in future fMRI brain studies.

Anomalous and heterogeneous characteristics of the BOLD hemodynamic response function in white matter

Detailed knowledge of the BOLD hemodynamic response function (HRF) is crucial for accurate analyses and interpretation of functional MRI data. Considerable efforts have been made to characterize the HRF in gray matter (GM), but much less attention has been paid to BOLD effects in white matter (WM). However, several recent reports have demonstrated reliable detection and analyses of WM BOLD signals both after stimulation and in a resting state. WM and GM differ in composition, energy requirements, and blood flow, so their neurovascular couplings also may well be different. We aimed to derive a comprehensive characterization of the HRF in WM across a population, including accurate measurements of its shape and its variation along and between WM pathways, using resting-state fMRI acquisitions. Our results show that the HRF is significantly different between WM and GM. Features of the HRF, such as a prominent initial dip, show strong relationships with features of the tissue microstructure derived from diffusion imaging, and these relationships differ between WM and GM, consistent with BOLD signal fluctuations reflecting different energy demands and neurovascular couplings in tissues of different composition and function. We also show that the HRF varies in shape significantly along WM pathways and is different between different WM pathways, suggesting the temporal evolution of BOLD signals after an event vary in different parts of the WM. These features of the HRF in WM are especially relevant for interpretation of the biophysical basis of BOLD effects in WM.

Temporal stability of the hemodynamic response function across the majority of human cerebral cortex

The hemodynamic response function (HRF) measured with functional magnetic resonance imaging is generated by vascular and metabolic responses evoked by brief (<4 s) stimuli. It is known that the human HRF varies across cortex, between subjects, with stimulus paradigms, and even between different measurements in the same cortical location. However, our results demonstrate that strong HRFs are remarkably repeatable across sessions separated by time intervals up to 3 months. In this study, a multisensory stimulus was used to activate and measure the HRF across the majority of cortex (>70%, with lesser reliability observed in some areas of prefrontal cortex). HRFs were measured with high spatial resolution (2‐mm voxels) in central gray matter to minimize variations caused by partial‐volume effects. HRF amplitudes and temporal dynamics were highly repeatable across four sessions in 20 subjects. Positive and negative HRFs were consistently observed across sessions and subjects. Negative HRFs were generally weaker and, thus, more variable than positive HRFs. Statistical measurements showed that across‐session variability is highly correlated to the variability across events within a session; these measurements also indicated a normal distribution of variability across cortex. The overall repeatability of the HRFs over long time scales generally supports the long‐term use of event‐related functional magnetic resonance imaging protocols.

The Atypical Effective Connectivity of Right Temporoparietal Junction in Autism Spectrum Disorder: A Multi-Site Study

Social function impairment is the core deficit of autism spectrum disorder (ASD). Although many studies have investigated ASD through a variety of neuroimaging tools, its brain mechanism of social function remains unclear due to its complex and heterogeneous symptoms. The present study aimed to use resting-state functional magnetic imaging data to explore effective connectivity between the right temporoparietal junction (RTPJ), one of the key brain regions associated with social impairment of individuals with ASD, and the whole brain to further deepen our understanding of the neuropathological mechanism of ASD. This study involved 1,454 participants from 23 sites from the Autism Brain Imaging Data Exchange (ABIDE) public dataset, which included 618 individuals with ASD and 836 with typical development (TD). First, a voxel-wise Granger causality analysis (GCA) was conducted with the RTPJ selected as the region of interest (ROI) to investigate the differences in effective connectivity between the ASD and TD groups in every site. Next, to obtain further accurate and representative results, an image-based meta-analysis was implemented to further analyze the GCA results of each site. Our results demonstrated abnormal causal connectivity between the RTPJ and the widely distributed brain regions and that the connectivity has been associated with social impairment in individuals with ASD. The current study could help to further elucidate the pathological mechanisms of ASD and provides a new perspective for future research.

Effect of Intranasal Oxytocin on Resting-state Effective Connectivity in Schizophrenia

OBJECTIVES

Evidence from several lines of research suggests the critical role of neuropeptide oxytocin in social cognition and social behavior. Though a few studies have examined the effect of oxytocin on clinical symptoms of schizophrenia, the underlying neurobiological changes are underexamined. Hence, in this study, we examined the effect of oxytocin on the brain's effective connectivity in schizophrenia.

METHODS

31 male patients with schizophrenia (SCZ) and 21 healthy male volunteers (HV) underwent resting functional magnetic resonance imaging scans with intra-nasal oxytocin (24 IU) and placebo administered in counterbalanced order. We conducted a whole-brain effective connectivity analysis using a multivariate vector autoregressive granger causality model. We performed a conjunction analysis to control for spurious changes and canonical correlation analysis between changes in connectivity and clinical and demographic variables.

RESULTS

Three connections, sourced from the left caudate survived the FDR correction threshold with the conjunction analysis; connections to the left supplementary motor area, left precentral gyrus, and left frontal inferior triangular gyrus. At baseline, SCZ patients had significantly weaker connectivity from caudate to these three regions. Oxytocin, but not placebo, significantly increased the strength of connectivity in these connections. Better cognitive insight and lower negative symptoms were associated with a greater increase in connectivity with oxytocin.

CONCLUSIONS

These findings provide a preliminary mechanistic understanding of the effect of oxytocin on brain connectivity in schizophrenia. The study findings provide the rationale to examine the potential utility of oxytocin for social cognitive deficits in schizophrenia.

Estrogen, Cognitive Performance, and Functional Imaging Studies: What Are We Missing About Neuroprotection?

Menopause transition can be interpreted as a vulnerable state characterized by estrogen deficiency with detrimental systemic effects as the low-grade chronic inflammation that appears with aging and partly explains age-related disorders as cancer, diabetes mellitus and increased risk of cognitive impairment. Over the course of a lifetime, estrogen produces several beneficial effects in healthy neurological tissues as well as cardioprotective effects, and anti-inflammatory effects. However, clinical evidence on the efficacy of hormone treatment in menopausal women has failed to confirm the benefit reported in observational studies. Unambiguously, enhanced verbal memory is the most robust finding from longitudinal and cross-sectional studies, what merits consideration for future studies aiming to determine estrogen neuroprotective efficacy. Estrogen related brain activity and functional connectivity remain, however, unexplored. In this context, the resting state paradigm may provide valuable information about reproductive aging and hormonal treatment effects, and their relationship with brain imaging of functional connectivity may be key to understand and anticipate estrogen cognitive protective effects. To go in-depth into the molecular and cellular mechanisms underlying rapid-to-long lasting protective effects of estrogen, we will provide a comprehensive review of cognitive tasks used in animal studies to evaluate the effect of hormone treatment on cognitive performance and discuss about the tasks best suited to the demonstration of clinically significant differences in cognitive performance to be applied in human studies. Eventually, we will focus on studies evaluating the DMN activity and responsiveness to pharmacological stimulation in humans.

Characterization of Hemodynamic Alterations in Schizophrenia and Bipolar Disorder and Their Effect on Resting-State fMRI Functional Connectivity

Common and distinct neural bases of Schizophrenia (SZ) and bipolar disorder (BP) have been explored using resting-state fMRI (rs-fMRI) functional connectivity (FC). However, fMRI is an indirect measure of neural activity, which is a convolution of the hemodynamic response function (HRF) and latent neural activity. The HRF, which models neurovascular coupling, varies across the brain within and across individuals, and is altered in many psychiatric disorders. Given this background, this study had three aims: quantifying HRF aberrations in SZ and BP, measuring the impact of such HRF aberrations on FC group differences, and exploring the genetic basis of HRF aberrations. We estimated voxel-level HRFs by deconvolving rs-fMRI data obtained from SZ (N = 38), BP (N = 19), and matched healthy controls (N = 35). We identified HRF group differences (P < .05, FDR corrected) in many regions previously implicated in SZ/BP, with mediodorsal, habenular, and central lateral nuclei of the thalamus exhibiting HRF differences in all pairwise group comparisons. Thalamus seed-based FC analysis revealed that ignoring HRF variability results in false-positive and false-negative FC group differences, especially in insula, superior frontal, and lingual gyri. HRF was associated with DRD2 gene expression (P < .05, 1.62 < |Z| < 2.0), as well as with medication dose (P < .05, 1.75 < |Z| < 3.25). In this first study to report HRF aberrations in SZ and BP, we report the possible modulatory effect of dopaminergic signalling on HRF, and the impact that HRF variability can have on FC studies in clinical samples. To mitigate the impact of HRF variability on FC group differences, we suggest deconvolution during data preprocessing.

Modeling the Hemodynamic Response Function Using EEG-fMRI Data During Eyes-Open Resting-State Conditions and Motor Task Execution

Being able to accurately quantify the hemodynamic response function (HRF) that links the blood oxygen level dependent functional magnetic resonance imaging (BOLD-fMRI) signal to the underlying neural activity is important both for elucidating neurovascular coupling mechanisms and improving the accuracy of fMRI-based functional connectivity analyses. In particular, HRF estimation using BOLD-fMRI is challenging particularly in the case of resting-state data, due to the absence of information about the underlying neuronal dynamics. To this end, using simultaneously recorded electroencephalography (EEG) and fMRI data is a promising approach, as EEG provides a more direct measure of neural activations. In the present work, we employ simultaneous EEG-fMRI to investigate the regional characteristics of the HRF using measurements acquired during resting conditions. We propose a novel methodological approach based on combining distributed EEG source space reconstruction, which improves the spatial resolution of HRF estimation and using block-structured linear and nonlinear models, which enables us to simultaneously obtain HRF estimates and the contribution of different EEG frequency bands. Our results suggest that the dynamics of the resting-state BOLD signal can be sufficiently described using linear models and that the contribution of each band is region specific. Specifically, it was found that sensory-motor cortices exhibit positive HRF shapes, whereas the lateral occipital cortex and areas in the parietal cortex, such as the inferior and superior parietal lobule exhibit negative HRF shapes. To validate the proposed method, we repeated the analysis using simultaneous EEG-fMRI measurements acquired during execution of a unimanual hand-grip task. Our results reveal significant associations between BOLD signal variations and electrophysiological power fluctuations in the ipsilateral primary motor cortex, particularly for the EEG beta band, in agreement with previous studies in the literature.

Electronic Cigarette Vaping Did Not Enhance the Neural Process of Working Memory for Regular Cigarette Smokers

Background

Electronic cigarettes (e-cigs) as substitute devices for regular tobacco cigarettes (r-cigs) have been increasing in recent times. We investigated neuronal substrates of vaping e-cigs and smoking r-cigs from r-cig smokers.

Methods

Twenty-two r-cig smokers made two visits following overnight smoking cessation. Functional magnetic resonance imaging (fMRI) data were acquired while participants watched smoking images. Participants were then allowed to smoke either an e-cig or r-cig until satiated and fMRI data were acquired. Their craving levels and performance on the Montreal Imaging Stress Task and a 3-back alphabet/digit recognition task were obtained and analyzed using two-way repeated-measures analysis of variance. Regions-of-interest (ROIs) were identified by comparing the abstained and satiated conditions. Neuronal activation within ROIs was regressed on the craving and behavioral data separately.

Results

Craving was more substantially reduced by smoking r-cigs than by vaping e-cigs. The response time (RT) for the 3-back task was significantly shorter following smoking r-cigs than following vaping e-cigs (interaction: F (1, 17) = 5.3, p = 0.035). Neuronal activations of the right vermis (r = 0.43, p = 0.037, CI = [-0.05, 0.74]), right caudate (r = 0.51, p = 0.015, CI = [0.05, 0.79]), and right superior frontal gyrus (r = −0.70, p = 0.001, CI = [−0.88, −0.34]) were significantly correlated with the RT for the 3-back task only for smoking r-cigs.

Conclusion

Our findings suggest that insufficient satiety from vaping e-cigs for r-cigs smokers may be insignificant effect on working memory function.

Changes in Cortical Activity in Stroke Survivors Undergoing Botulinum Neurotoxin Therapy for Treatment of Focal Spasticity

Background: Botulinum NeuroToxin-A (BoNT-A) relieves muscle spasticity and increases range of motion necessary for stroke rehabilitation. Determining the effects of BoNT-A therapy on brain neuroplasticity could help physicians customize its use and predict its outcome.

Objective: The purpose of this study was to investigate the effects of Botulinum Toxin-A therapy for treatment of focal spasticity on brain activation and functional connectivity.

Design: We used functional Magnetic Resonance Imaging (fMRI) to track changes in blood oxygen-level dependent (BOLD) activation and functional connectivity associated with BoNT-A therapy in nine chronic stroke participants, and eight age-matched controls. Scans were acquired before BoNT-A injections (W0) and 6 weeks after the injections (W6). The task fMRI scan consisted of a block design of alternating mass finger flexion and extension. The voxel-level changes in BOLD activation, and pairwise changes in functional connectivity were analyzed for BoNT-A treatment (stroke W0 vs. W6).

Results: BoNT-A injection therapy resulted in significant increases in brain activation in the contralesional premotor cortex, cingulate gyrus, thalamus, superior cerebellum, and in the ipsilesional sensory integration area. Lastly, cerebellar connectivity correlated with the Fugl-Meyer assessment of motor impairment before injection, while premotor connectivity correlated with the Fugl-Meyer score after injection.

Conclusion: BoNT-A therapy for treatment of focal spasticity resulted in increased brain activation in areas associated with motor control, and cerebellar connectivity correlated with motor impairment before injection. These results suggest that neuroplastic effects might take place in response to improvements in focal spasticity.

Functional Connectivity-Based Prediction of Autism on Site Harmonized ABIDE Dataset

OBJECTIVE

The larger sample sizes available from multi-site publicly available neuroimaging data repositories makes machine-learning based diagnostic classification of mental disorders more feasible by alleviating the curse of dimensionality. However, since multi-site data are aggregated post-hoc, i.e. they were acquired from different scanners with different acquisition parameters, non-neural inter-site variability may mask inter-group differences that are at least in part neural in origin. Hence, the advantages gained by the larger sample size in the context of machine-learning based diagnostic classification may not be realized.

METHODS

We address this issue using harmonization of multi-site neuroimaging data using the ComBat technique, which is based on an empirical Bayes formulation to remove inter-site differences in data distributions, to improve diagnostic classification accuracy. Specifically, we demonstrate this using ABIDE (Autism Brain Imaging Data Exchange) multi-site data for classifying individuals with Autism from healthy controls using resting state fMRI-based functional connectivity data.

RESULTS

Our results show that higher classification accuracies across multiple classification models can be obtained (especially for models based on artificial neural networks) from multi-site data post harmonization with the ComBat technique as compared to without harmonization, outperforming earlier results from existing studies using ABIDE. Furthermore, our network ablation analysis facilitated important insights into autism spectrum disorder pathology and the connectivity in networks shown to be important for classification covaried with verbal communication impairments in Autism.

CONCLUSION

Multi-site data harmonization using ComBat improves neuroimaging-based diagnostic classification of mental disorders.

SIGNIFICANCE

ComBat has the potential to make AI-based clinical decision-support systems more feasible in psychiatry.

rsHRF: A toolbox for resting-state HRF estimation and deconvolution

The hemodynamic response function (HRF) greatly influences the intra- and inter-subject variability of brain activation and connectivity, and might confound the estimation of temporal precedence in connectivity analyses, making its estimation necessary for a correct interpretation of neuroimaging studies. Additionally, the HRF shape itself is a useful local measure. However, most algorithms for HRF estimation are specific for task-related fMRI data, and only a few can be directly applied to resting-state protocols. Here we introduce rsHRF, a Matlab and Python toolbox that implements HRF estimation and deconvolution from the resting-state BOLD signal. We first provide an overview of the main algorithm, practical implementations, and then demonstrate the feasibility and usefulness of rsHRF by validation experiments with a publicly available resting-state fMRI dataset. We also provide tools for statistical analyses and visualization. We believe that this toolbox may significantly contribute to a better analysis and understanding of the components and variability of BOLD signals.

Revisiting Nonlinear Functional Brain Co-activations: Directed, Dynamic, and Delayed

The center stage of neuro-imaging is currently occupied by studies of functional correlations between brain regions. These correlations define the brain functional networks, which are the most frequently used framework to represent and interpret a variety of experimental findings. In the previous study, we first demonstrated that the relatively stronger blood oxygenated level dependent (BOLD) activations contain most of the information relevant to understand functional connectivity, and subsequent work confirmed that a large compression of the original signals can be obtained without significant loss of information. In this study, we revisit the correlation properties of these epochs to define a measure of nonlinear dynamic directed functional connectivity (nldFC) across regions of interest. We show that the proposed metric provides at once, without extensive numerical complications, directed information of the functional correlations, as well as a measure of temporal lags across regions, overall offering a different and complementary perspective in the analysis of brain co-activation patterns. In this study, we provide further details for the computations of these measures and for a proof of concept based on replicating existing results from an Autistic Syndrome database, and discuss the main features and advantages of the proposed strategy for the study of brain functional correlations.

Investigating mechanisms of fast BOLD responses: The effects of stimulus intensity and of spatial heterogeneity of hemodynamics

Recent studies have demonstrated that fast fMRI can track neural activity well above the temporal limit predicted by the canonical hemodynamic response model. While these findings are promising, the biophysical mechanisms underlying these fast fMRI phenomena remain underexplored. In this study, we discuss two aspects of the hemodynamic response, complementary to several existing hypotheses, that can accommodate faster fMRI dynamics beyond those predicted by the canonical model. First, we demonstrate, using both visual and somatosensory paradigms, that the timing and shape of hemodynamic response functions (HRFs) vary across graded levels of stimulus intensity-with lower-intensity stimulation eliciting faster and narrower HRFs. Second, we show that as the spatial resolution of fMRI increases, voxel-wise HRFs begin to deviate from the canonical model, with a considerable portion of voxels exhibiting faster temporal dynamics than predicted by the canonical HRF. Collectively, both stimulus/task intensity and image resolution can affect the sensitivity of fMRI to fast brain activity, which may partly explain recent observations of fast fMRI signals. It is further noteworthy that, while the present investigations focus on fast neural responses, our findings suggest that a revised hemodynamic model may benefit the many fMRI studies using paradigms with wide ranges of contrast levels (e.g., resting or naturalistic conditions) or with modern, high-resolution MR acquisitions.

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.

Structurally constrained effective brain connectivity

The relationship between structure and function is of interest in many research fields involving the study of complex biological processes. In neuroscience in particular, the fusion of structural and functional data can help to understand the underlying principles of the operational networks in the brain. To address this issue, this paper proposes a constrained autoregressive model leading to a representation of effective connectivity that can be used to better understand how the structure modulates the function. Or simply, it can be used to find novel biomarkers characterizing groups of subjects. In practice, an initial structural connectivity representation is re-weighted to explain the functional co-activations. This is obtained by minimizing the reconstruction error of an autoregressive model constrained by the structural connectivity prior. The model has been designed to also include indirect connections, allowing to split direct and indirect components in the functional connectivity, and it can be used with raw and deconvoluted BOLD signal. The derived representation of dependencies was compared to the well known dynamic causal model, giving results closer to known ground-truth. Further evaluation of the proposed effective network was performed on two typical tasks. In a first experiment the direct functional dependencies were tested on a community detection problem, where the brain was partitioned using the effective networks across multiple subjects. In a second experiment the model was validated in a case-control task, which aimed at differentiating healthy subjects from individuals with autism spectrum disorder. Results showed that using effective connectivity leads to clusters better describing the functional interactions in the community detection task, while maintaining the original structural organization, and obtaining a better discrimination in the case-control classification task.

Hemodynamic matrix factorization for functional magnetic resonance imaging

The General Linear Model (GLM) used in task-fMRI relates activated brain areas to extrinsic task conditions. The translation of resulting neural activation into a hemodynamic response is commonly approximated with a linear convolution model using a hemodynamic response function (HRF). There are two major limitations in GLM analysis. Firstly, the GLM assumes that neural activation is either on or off and matches the exact stimulus duration in the corresponding task timings. Secondly, brain networks observed in resting-state fMRI experiments present also during task experiments, but the GLM approach models these task-unrelated brain activity as noise. A novel kernel matrix factorization approach, called hemodynamic matrix factorization (HMF), is therefore proposed that addresses both limitations by assuming that task-related and task-unrelated brain activity can be modeled with the same convolution model as in GLM analysis. By contrast to the GLM, the proposed HMF is a blind source separation (BSS) technique, which decomposes fMRI data into modes. Each mode comprises of a neural activation time course and a spatial mapping. Two versions of HMF are proposed in which the neural activation time course of each mode is convolved with either the canonical HRF or predetermined subject-specific HRFs. Firstly, HMF with the canonical HRF is applied to two open-source cohorts. These cohorts comprise of several task experiments including motor, incidental memory, spatial coherence discrimination, verbal discrimination task and a very short localization task, engaging multiple parts of the eloquent cortex. HMF modes were obtained whose neural activation time course followed original task timings and whose corresponding spatial map matched cortical areas known to be involved in the respective task processing. Secondly, the alignment of these neural activation time courses to task timings were further improved by replacing the canonical HRF with subject-specific HRFs during HMF mode computation. In addition to task-related modes, HMF also produced seemingly task-unrelated modes whose spatial maps matched known resting-state networks. The validity of a fMRI task experiment relies on the assumption that the exposure to a stimulus for a given time causes an imminent increase in neural activation of equal duration. The proposed HMF is an attempt to falsify this assumption and allows to identify subject task participation that does not comply with the experiment instructions.

Connective profiles and antagonism between dynamic and static connectivity underlying generalized epilepsy

This study aims to characterize the connective profiles and the coupling relationship between dynamic and static functional connectivity (dFC and sFC) in large-scale brain networks in patients with generalized epilepsy (GE). Functional, structural and diffuse MRI data were collected from 83 patients with GE and 106 matched healthy controls (HC). Resting-state BOLD time course was deconvolved to neural time course using a blind hemodynamic deconvolution method. Then, five connective profiles, including the structural connectivity (SC) and BOLD/neural time course-derived sFC/dFC networks, were constructed based on the proposed whole brain atlas. Network-level weighted correlation probability (NWCP) were proposed to evaluate the association between dFC and sFC. Both the BOLD signal and neural time course showed highly concordant findings and the present study emphasized the consistent findings between two functional approaches. The patients with GE showed hypervariability and enhancement of FC, and notably decreased SC in the subcortical network. Besides, increased dFC, weaker anatomic links, and complex alterations of sFC were observed in the default mode network of GE. Moreover, significantly increased SC and predominantly increased sFC were found in the frontoparietal network. Remarkably, antagonism between dFC and sFC was observed in large-scale networks in HC, while patients with GE showed significantly decreased antagonism in core epileptic networks. In sum, our study revealed distinct connective profiles in different epileptic networks and provided new clues to the brain network mechanism of epilepsy from the perspective of antagonism between dynamic and static functional connectivity.

Reconfiguration of Directed Functional Connectivity Among Neurocognitive Networks with Aging: Considering the Role of Thalamo-Cortical Interactions

A complete picture of how subcortical nodes, such as the thalamus, exert directional influence on large-scale brain network interactions across age remains elusive. Using directed functional connectivity and weighted net causal outflow on resting-state fMRI data, we provide evidence of a comprehensive reorganization within and between neurocognitive networks (default mode: DMN, salience: SN, and central executive: CEN) associated with age and thalamocortical interactions. We hypothesize that thalamus subserves both modality-specific and integrative hub role in organizing causal weighted outflow among large-scale neurocognitive networks. To this end, we observe that within-network directed functional connectivity is driven by thalamus and progressively weakens with age. Secondly, we find that age-associated increase in between CEN- and DMN-directed functional connectivity is driven by both the SN and the thalamus. Furthermore, left and right thalami act as a causal integrative hub exhibiting substantial interactions with neurocognitive networks with aging and play a crucial role in reconfiguring network outflow. Notably, these results were largely replicated on an independent dataset of matched young and old individuals. Our findings strengthen the hypothesis that the thalamus is a key causal hub balancing both within- and between-network connectivity associated with age and maintenance of cognitive functioning with aging.

Distinction of High- and Low-Frequency Repetitive Transcranial Magnetic Stimulation on the Functional Reorganization of the Motor Network in Stroke Patients

Objective

To investigate the functional reorganization of the motor network after repetitive transcranial magnetic stimulation (rTMS) in stroke patients with motor dysfunction and the distinction between high-frequency rTMS (HF-rTMS) and low-frequency rTMS (LF-rTMS).

Methods

Thirty-three subcortical stroke patients were enrolled and assigned to the HF-rTMS group, LF-rTMS group, and sham group. Each patient of rTMS groups received either 10.0 Hz rTMS over the ipsilesional primary motor cortex (M1) or 1.0 Hz rTMS over the contralesional M1 for 10 consecutive days. A resting-state functional magnetic resonance imaging (fMRI) scan and neurological examinations were performed at baseline and after rTMS. The motor network and functional connectivities intramotor network with the core brain regions including the bilateral M1, premotor area (PMA), and supplementary motor area (SMA) were calculated. Comparisons of functional connectivities and Pearson correlation analysis between functional connectivity changes and behavioral improvement were calculated.

Results

Significant motor improvement was found after rTMS in all groups which was larger in two rTMS groups than in the sham group. The functional connectivities of the motor network were significantly increased in bilateral M1, SMA, and contralesional PMA after real rTMS. These changes were only detected in the regions of the ipsilesional hemisphere in the HF-rTMS group and in the regions of the contralesional hemisphere in the LF-rTMS group. Significantly changed functional connectivities of the intramotor network were found between the ipsilesional M1 and SMA and contralesional PMA, between contralesional M1 and contralesional SMA, between contralesional SMA and ipsilesional SMA and contralesional PMA in the HF-rTMS group in which the changed connectivity between ipsilesional M1 and contralesional PMA was obviously correlated with the motor improvement. In addition, the functional connectivity of the intramotor network between ipsilesional M1 and contralesional PMA was significantly higher in the HF-rTMS group than in the LF-rTMS group.

Conclusion

Both HF-rTMS and LF-rTMS have a positive effect on motor recovery in patients with subcortical stroke and could promote the reorganization of the motor network. HF-rTMS may contribute more to the functional connectivity reorganization of the ipsilesional motor network and realize greater benefit to the motor recovery.

Dynamic time warping outperforms Pearson correlation in detecting atypical functional connectivity in autism spectrum disorders

Resting state fMRI (rsfMRI) is frequently used to study brain function, including in clinical populations. Similarity of blood-oxygen-level-dependent (BOLD) fluctuations during rsfMRI between brain regions is thought to reflect intrinsic functional connectivity (FC), potentially due to history of coactivation. To quantify similarity, studies have almost exclusively relied on Pearson correlation, which assumes linearity and can therefore underestimate FC if the hemodynamic response function differs regionally or there is BOLD signal lag between timeseries. Here we show in three cohorts of children, adolescents and adults, with and without autism spectrum disorders (ASDs), that measuring the similarity of BOLD signal fluctuations using non-linear dynamic time warping (DTW) is more robust to global signal regression (GSR), has higher test-retest reliability and is more sensitive to task-related changes in FC. Additionally, when comparing FC between individuals with ASDs and typical controls, more group differences are detected using DTW. DTW estimates are also more related to ASD symptom severity and executive function, while Pearson correlation estimates of FC are more strongly associated with respiration during rsfMRI. Together these findings suggest that non-linear methods such as DTW improve estimation of resting state FC, particularly when studying clinical populations whose hemodynamics or neurovascular coupling may be altered compared to typical controls.

Dynamic brain-to-brain concordance and behavioral mirroring as a mechanism of the patient-clinician interaction

The patient-clinician interaction can powerfully shape treatment outcomes such as pain but is often considered an intangible "art of medicine" and has largely eluded scientific inquiry. Although brain correlates of social processes such as empathy and theory of mind have been studied using single-subject designs, specific behavioral and neural mechanisms underpinning the patient-clinician interaction are unknown. Using a two-person interactive design, we simultaneously recorded functional magnetic resonance imaging (hyperscanning) in patient-clinician dyads, who interacted via live video, while clinicians treated evoked pain in patients with chronic pain. Our results show that patient analgesia is mediated by patient-clinician nonverbal behavioral mirroring and brain-to-brain concordance in circuitry implicated in theory of mind and social mirroring. Dyad-based analyses showed extensive dynamic coupling of these brain nodes with the partners' brain activity, yet only in dyads with pre-established clinical rapport. These findings introduce a putatively key brain-behavioral mechanism for therapeutic alliance and psychosocial analgesia.

Functional connectome organization predicts conversion to psychosis in clinical high-risk youth from the SHARP program

The emergence of prodromal symptoms of schizophrenia and their evolution into overt psychosis may stem from an aberrant functional reorganization of the brain during adolescence. To examine whether abnormalities in connectome organization precede psychosis onset, we performed a functional connectome analysis in a large cohort of medication-naive youth at risk for psychosis from the Shanghai At Risk for Psychosis (SHARP) study. The SHARP program is a longitudinal study of adolescents and young adults at Clinical High Risk (CHR) for psychosis, conducted at the Shanghai Mental Health Center in collaboration with neuroimaging laboratories at Harvard and MIT. Our study involved a total of 251 subjects, including 158 CHRs and 93 age-, sex-, and education-matched healthy controls. During 1-year follow-up, 23 CHRs developed psychosis. CHRs who would go on to develop psychosis were found to show abnormal modular connectome organization at baseline, while CHR non-converters did not. In all CHRs, abnormal modular connectome organization at baseline was associated with a threefold conversion rate. A region-specific analysis showed that brain regions implicated in early-course schizophrenia, including superior temporal gyrus and anterior cingulate cortex, were most abnormal in terms of modular assignment. Our results show that functional changes in brain network organization precede the onset of psychosis and may drive psychosis development in at-risk youth.

Quality and denoising in real-time functional magnetic resonance imaging neurofeedback: A methods review

Neurofeedback training using real-time functional magnetic resonance imaging (rtfMRI-NF) allows subjects voluntary control of localised and distributed brain activity. It has sparked increased interest as a promising non-invasive treatment option in neuropsychiatric and neurocognitive disorders, although its efficacy and clinical significance are yet to be determined. In this work, we present the first extensive review of acquisition, processing and quality control methods available to improve the quality of the neurofeedback signal. Furthermore, we investigate the state of denoising and quality control practices in 128 recently published rtfMRI-NF studies. We found: (a) that less than a third of the studies reported implementing standard real-time fMRI denoising steps, (b) significant room for improvement with regards to methods reporting and (c) the need for methodological studies quantifying and comparing the contribution of denoising steps to the neurofeedback signal quality. Advances in rtfMRI-NF research depend on reproducibility of methods and results. Notably, a systematic effort is needed to build up evidence that disentangles the various mechanisms influencing neurofeedback effects. To this end, we recommend that future rtfMRI-NF studies: (a) report implementation of a set of standard real-time fMRI denoising steps according to a proposed COBIDAS-style checklist (https://osf.io/kjwhf/), (b) ensure the quality of the neurofeedback signal by calculating and reporting community-informed quality metrics and applying offline control checks and (c) strive to adopt transparent principles in the form of methods and data sharing and support of open-source rtfMRI-NF software. Code and data for reproducibility, as well as an interactive environment to explore the study data, can be accessed at https://github.com/jsheunis/quality-and-denoising-in-rtfmri-nf.