Whole-brain, time-locked activation with simple tasks revealed using massive averaging and model-free analysis

Mental state decoders: game-changers or wishful thinking?

Decoding mental and perceptual states using fMRI has become increasingly popular over the past two decades, with numerous highly-cited studies published in high-profile journals. Nevertheless, what have we learned from these decoders? In this opinion, we argue that fMRI-based decoders are not neurophysiologically informative and are not, and likely cannot be, applicable to real-world decision-making. The former point stems from the fact that decoding models cannot disentangle neural mechanisms from their epiphenomena. The latter point stems from both logical and ethical constraints. Constructing decoders requires precious time and resources that should instead be directed toward scientific endeavors more likely to yield meaningful scientific progress.

Cortical Networks Relating to Arousal Are Differentially Coupled to Neural Activity and Hemodynamics

Even in the absence of specific sensory input or a behavioral task, the brain produces structured patterns of activity. This organized activity is modulated by changes in arousal. Here, we use wide-field voltage imaging to establish how arousal relates to cortical network voltage and hemodynamic activity in spontaneously behaving head-fixed male and female mice expressing the voltage-sensitive fluorescent FRET sensor Butterfly 1.2. We find that global voltage and hemodynamic signals are both positively correlated with changes in arousal with a maximum correlation of 0.5 and 0.25, respectively, at a time lag of 0 s. We next show that arousal influences distinct cortical regions for both voltage and hemodynamic signals. These include a broad positive correlation across most sensory-motor cortices extending posteriorly to the primary visual cortex observed in both signals. In contrast, activity in the prefrontal cortex is positively correlated to changes in arousal for the voltage signal while it is a slight net negative correlation observed in the hemodynamic signal. Additionally, we show that coherence between voltage and hemodynamic signals relative to arousal is strongest for slow frequencies below 0.15 Hz and is near zero for frequencies >1 Hz. We finally show that coupling patterns are dependent on the behavioral state of the animal with correlations being driven by periods of increased orofacial movement. Our results indicate that while hemodynamic signals show strong relations to behavior and arousal, these relations are distinct from those observed by voltage activity.

The Primate Cortical LFP Exhibits Multiple Spectral and Temporal Gradients and Widespread Task-Dependence During Visual Short-Term Memory

Although cognitive functions are hypothesized to be mediated by synchronous neuronal interactions in multiple frequency bands among widely distributed cortical areas, we still lack a basic understanding of the distribution and task dependence of oscillatory activity across the cortical map. Here, we ask how the spectral and temporal properties of the local field potential (LFP) vary across the primate cerebral cortex, and how they are modulated during visual short-term memory. We measured the LFP from 55 cortical areas in two macaque monkeys while they performed a visual delayed match to sample task. Analysis of peak frequencies in the LFP power spectra reveals multiple discrete frequency bands between 3-80 Hz that differ between the two monkeys. The LFP power in each band, as well as the Sample Entropy, a measure of signal complexity, display distinct spatial gradients across the cortex, some of which correlate with reported spine counts in layer 3 pyramidal neurons. Cortical areas can be robustly decoded using a small number of spectral and temporal parameters, and significant task dependent increases and decreases in spectral power occur in all cortical areas. These findings reveal pronounced, widespread and spatially organized gradients in the spectral and temporal activity of cortical areas. Task-dependent changes in cortical activity are globally distributed, even for a simple cognitive task.

State-dependent memory retrieval: insights from neural dynamics and behavioral perspectives

Memory retrieval is strikingly susceptible to external states (environment) and internal states (mood states and alcohol), yet we know little about the underlying mechanisms. We examined how internally generated states influence successful memory retrieval using the functional magnetic resonance imaging (fMRI) of laboratory mice during memory retrieval. Mice exhibited a strong tendency to perform memory retrieval correctly only in the reinstated mammillary body-inhibited state, in which mice were trained to discriminate auditory stimuli in go/no-go tasks. fMRI revealed that distinct auditory cues engaged differential brain regions, which were primed by internal state. Specifically, a cue associated with a reward activated the lateral amygdala, while a cue signaling no reward predominantly activated the postsubiculum. Modifying these internal states significantly altered the neural activity balance between these regions. Optogenetic inhibition of those regions in the precue period blocked the retrieval of type-specific memories. Our findings suggest that memory retrieval is under the control of two interrelated neural circuits underlying the neural basis of state-dependent memory retrieval.

fMRI with whole-brain coverage, 75-ms temporal resolution and high SNR by combining HiHi reshuffling and multiband imaging

Increasing the temporal resolution of the blood‑oxygen level-dependent (BOLD) response is usually accompanied by a decrease in repetition time and therefore also a reduction of the magnetic resonance (MR) signal due to incomplete T1 relaxation and thus a loss of signal-to-noise ratio (SNR). A previous data reordering method can achieve higher temporal sampling rate without the loss of SNR but at the cost of increased scan time. In this proof-of-principle work, we show that combining HiHi reshuffling with multiband acceleration allows us to measure the in vivo BOLD response with a 75-ms sampling rate that is decoupled from the acquisition repetition time (here 1.5 s and hence higher SNR) while covering the entire forebrain with 60 2-mm slices in a ~ 35-min scan. We provide single-voxel time-courses of the BOLD responses in the primary visual and primary motor cortices in three fMRI experiments on a 7 T scanner - 1 male (scanned twice on different days for test-retest reproducibility) and 1 female participant.

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.

An evidence-based critical review of the mind-brain identity theory

In the philosophy of mind, neuroscience, and psychology, the causal relationship between phenomenal consciousness, mentation, and brain states has always been a matter of debate. On the one hand, material monism posits consciousness and mind as pure brain epiphenomena. One of its most stringent lines of reasoning relies on a 'loss-of-function lesion premise,' according to which, since brain lesions and neurochemical modifications lead to cognitive impairment and/or altered states of consciousness, there is no reason to doubt the mind-brain identity. On the other hand, dualism or idealism (in one form or another) regard consciousness and mind as something other than the sole product of cerebral activity pointing at the ineffable, undefinable, and seemingly unphysical nature of our subjective qualitative experiences and its related mental dimension. Here, several neuroscientific findings are reviewed that question the idea that posits phenomenal experience as an emergent property of brain activity, and argue that the premise of material monism is based on a logical correlation-causation fallacy. While these (mostly ignored) findings, if considered separately from each other, could, in principle, be recast into a physicalist paradigm, once viewed from an integral perspective, they substantiate equally well an ontology that posits mind and consciousness as a primal phenomenon.

Whole-brain, gray, and white matter time-locked functional signal changes with simple tasks and model-free analysis

Recent studies have revealed the production of time-locked blood oxygenation level-dependent (BOLD) functional MRI (fMRI) signals throughout the entire brain in response to tasks, challenging the existence of sparse and localized brain functions and highlighting the pervasiveness of potential false negative fMRI findings. "Whole-brain" actually refers to gray matter, the only tissue traditionally studied with fMRI. However, several reports have demonstrated reliable detection of BOLD signals in white matter, which have previously been largely ignored. Using simple tasks and analyses, we demonstrate BOLD signal changes across the whole brain, in both white and gray matters, in similar manner to previous reports of whole brain studies. We investigated whether white matter displays time-locked BOLD signals across multiple structural pathways in response to a stimulus in a similar manner to the cortex. We find that both white and gray matter show time-locked activations across the whole brain, with a majority of both tissue types showing statistically significant signal changes for all task stimuli investigated. We observed a wide range of signal responses to tasks, with different regions showing different BOLD signal changes to the same task. Moreover, we find that each region may display different BOLD responses to different stimuli. Overall, we present compelling evidence that, just like all gray matter, essentially all white matter in the brain shows time-locked BOLD signal changes in response to multiple stimuli, challenging the idea of sparse functional localization and the prevailing wisdom of treating white matter BOLD signals as artifacts to be removed.

A functional account of stimulation-based aerobic glycolysis and its role in interpreting BOLD signal intensity increases in neuroimaging experiments

In aerobic glycolysis, oxygen is abundant, and yet cells metabolize glucose without using it, decreasing their ATP per glucose yield by 15-fold. During task-based stimulation, aerobic glycolysis occurs in localized brain regions, presenting a puzzle: why produce ATP inefficiently when, all else being equal, evolution should favor the efficient use of metabolic resources? The answer is that all else is not equal. We propose that a tradeoff exists between efficient ATP production and the efficiency with which ATP is spent to transmit information. Aerobic glycolysis, despite yielding little ATP per glucose, may support neuronal signaling in thin (< 0.5 µm), information-efficient axons. We call this the efficiency tradeoff hypothesis. This tradeoff has potential implications for interpretations of task-related BOLD "activation" observed in fMRI. We hypothesize that BOLD "activation" may index local increases in aerobic glycolysis, which support signaling in thin axons carrying "bottom-up" information, or "prediction error"-i.e., the BIAPEM (BOLD increases approximate prediction error metabolism) hypothesis. Finally, we explore implications of our hypotheses for human brain evolution, social behavior, and mental disorders.

Quantifying brain connectivity signatures by means of polyconnectomic scoring

A broad range of neuropsychiatric disorders are associated with alterations in macroscale brain circuitry and connectivity. Identifying consistent brain patterns underlying these disorders by means of structural and functional MRI has proven challenging, partly due to the vast number of tests required to examine the entire brain, which can lead to an increase in missed findings. In this study, we propose polyconnectomic score (PCS) as a metric designed to quantify the presence of disease-related brain connectivity signatures in connectomes. PCS summarizes evidence of brain patterns related to a phenotype across the entire landscape of brain connectivity into a subject-level score. We evaluated PCS across four brain disorders (autism spectrum disorder, schizophrenia, attention deficit hyperactivity disorder, and Alzheimer's disease) and 14 studies encompassing ~35,000 individuals. Our findings consistently show that patients exhibit significantly higher PCS compared to controls, with effect sizes that go beyond other single MRI metrics ([min, max]: Cohen's d = [0.30, 0.87], AUC = [0.58, 0.73]). We further demonstrate that PCS serves as a valuable tool for stratifying individuals, for example within the psychosis continuum, distinguishing patients with schizophrenia from their first-degree relatives (d = 0.42, p = 4 × 10-3, FDR-corrected), and first-degree relatives from healthy controls (d = 0.34, p = 0.034, FDR-corrected). We also show that PCS is useful to uncover associations between brain connectivity patterns related to neuropsychiatric disorders and mental health, psychosocial factors, and body measurements.

BOLD Response is more than just magnitude: Improving detection sensitivity through capturing hemodynamic profiles

Typical fMRI analyses often assume a canonical hemodynamic response function (HRF) that primarily focuses on the peak height of the overshoot, neglecting other morphological aspects. Consequently, reported analyses often reduce the overall response curve to a single scalar value. In this study, we take a data-driven approach to HRF estimation at the whole-brain voxel level, without assuming a response profile at the individual level. We then employ a roughness penalty at the population level to estimate the response curve, aiming to enhance predictive accuracy, inferential efficiency, and cross-study reproducibility. By examining a fast event-related FMRI dataset, we demonstrate the shortcomings and information loss associated with adopting the canonical approach. Furthermore, we address the following key questions: 1) To what extent does the HRF shape vary across different regions, conditions, and participant groups? 2) Does the data-driven approach improve detection sensitivity compared to the canonical approach? 3) Can analyzing the HRF shape help validate the presence of an effect in conjunction with statistical evidence? 4) Does analyzing the HRF shape offer evidence for whole-brain response during a simple task?

Physiological modeling of the BOLD signal and implications for effective connectivity: A primer

In this primer, I provide an overview of the physiological processes that contribute to the observed BOLD signal (i.e., the generative biophysical model), including their time course properties within the framework of the physiologically-informed dynamic causal modeling (P-DCM). The BOLD signal is primarily determined by the change in paramagnetic deoxygenated hemoglobin, which results from combination of changes in oxygen metabolism, and cerebral blood flow and volume. Specifically, the physiological origin of the so-called BOLD signal "transients" will be discussed, including the initial overshoot, steady-state activation and the post-stimulus undershoot. I argue that incorrect physiological assumptions in the generative model of the BOLD signal can lead to incorrect inferences pertaining to both local neuronal activity and effective connectivity between brain regions. In addition, I introduce the recent laminar BOLD signal model, which extends P-DCM to cortical depths-resolved BOLD signals, allowing for laminar neuronal activity to be determined using high-resolution fMRI data.

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 interaction between random and systematic visual stimulation and infraslow quasiperiodic spatiotemporal patterns of whole brain activity

One prominent feature of the infraslow BOLD signal during rest or task is quasi-periodic spatiotemporal pattern (QPP) of signal changes that involves an alternation of activity in key functional networks and propagation of activity across brain areas, and that is known to tie to the infraslow neural activity involved in attention and arousal fluctuations. This ongoing whole-brain pattern of activity might potentially modify the response to incoming stimuli or be modified itself by the induced neural activity. To investigate this, we presented checkerboard sequences flashing at 6Hz to subjects. This is a salient visual stimulus that is known to produce a strong response in visual processing regions. Two different visual stimulation sequences were employed, a systematic stimulation sequence in which the visual stimulus appeared every 20.3 secs and a random stimulation sequence in which the visual stimulus occurred randomly every 14~62.3 secs. Three central observations emerged. First, the two different stimulation conditions affect the QPP waveform in different aspects, i.e., systematic stimulation has greater effects on its phase and random stimulation has greater effects on its magnitude. Second, the QPP was more frequent in the systematic condition with significantly shorter intervals between consecutive QPPs compared to the random condition. Third, the BOLD signal response to the visual stimulus across both conditions was swamped by the QPP at the stimulus onset. These results provide novel insights into the relationship between intrinsic patterns and stimulated brain activity.

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.

Highlight results, don't hide them: Enhance interpretation, reduce biases and improve reproducibility

Most neuroimaging studies display results that represent only a tiny fraction of the collected data. While it is conventional to present "only the significant results" to the reader, here we suggest that this practice has several negative consequences for both reproducibility and understanding. This practice hides away most of the results of the dataset and leads to problems of selection bias and irreproducibility, both of which have been recognized as major issues in neuroimaging studies recently. Opaque, all-or-nothing thresholding, even if well-intentioned, places undue influence on arbitrary filter values, hinders clear communication of scientific results, wastes data, is antithetical to good scientific practice, and leads to conceptual inconsistencies. It is also inconsistent with the properties of the acquired data and the underlying biology being studied. Instead of presenting only a few statistically significant locations and hiding away the remaining results, studies should "highlight" the former while also showing as much as possible of the rest. This is distinct from but complementary to utilizing data sharing repositories: the initial presentation of results has an enormous impact on the interpretation of a study. We present practical examples and extensions of this approach for voxelwise, regionwise and cross-study analyses using publicly available data that was analyzed previously by 70 teams (NARPS; Botvinik-Nezer, et al., 2020), showing that it is possible to balance the goals of displaying a full set of results with providing the reader reasonably concise and "digestible" findings. In particular, the highlighting approach sheds useful light on the kind of variability present among the NARPS teams' results, which is primarily a varied strength of agreement rather than disagreement. Using a meta-analysis built on the informative "highlighting" approach shows this relative agreement, while one using the standard "hiding" approach does not. We describe how this simple but powerful change in practice-focusing on highlighting results, rather than hiding all but the strongest ones-can help address many large concerns within the field, or at least to provide more complete information about them. We include a list of practical suggestions for results reporting to improve reproducibility, cross-study comparisons and meta-analyses.

Strategic filtering of high-energy visible light expands neural correlates of functional vision particularly in older participants

In this study we assessed the neural correlates of functional vision while varying patterns of light filtration. Four filter conditions used relatively flat filtering across the visible spectrum while one filter was a step filter that selectively absorbed violet light (wavelengths below about 415 nm). Neural effects were quantified by measuring the BOLD response ((T2*-based fMRI) while subjects performed a challenging visual task (judging gap direction in Landolt Cs that randomly varied in size). In general (based on p < 0.01 directional criterion not corrected for aggregated error), as filtering increased (less interference by bright light), brain activity associated with the task also increased. This effect, even using the most conservative statistics, was most evident when using the violet filter (especially for the older subjects) despite only reducing the very highest energy portion of the visible spectrum. This finding suggests that filtering can increase neural activity associated with functional vision; such effects might be achievable through filtering just the highest visible energy (violet).

The interaction between random and systematic visual stimulation and infraslow quasiperiodic spatiotemporal patterns of whole brain activity

One prominent feature of the infraslow BOLD signal during rest or task is quasi-periodic spatiotemporal pattern (QPP) of signal changes that involves an alternation of activity in key functional networks and propagation of activity across brain areas, and that is known to tie to the infraslow neural activity involved in attention and arousal fluctuations. This ongoing whole-brain pattern of activity might potentially modify the response to incoming stimuli or be modified itself by the induced neural activity. To investigate this, we presented checkerboard sequences flashing at 6Hz to subjects. This is a salient visual stimulus that is known to produce a strong response in visual processing regions. Two different visual stimulation sequences were employed, a systematic stimulation sequence in which the visual stimulus appeared every 20.3 secs and a random stimulation sequence in which the visual stimulus occurred randomly every 14∼62.3 secs. Three central observations emerged. First, the two different stimulation conditions affect the QPP waveform in different aspects, i.e., systematic stimulation has greater effects on its phase and random stimulation has greater effects on its magnitude. Second, the QPP was more frequent in the systematic condition with significantly shorter intervals between consecutive QPPs compared to the random condition. Third, the BOLD signal response to the visual stimulus across both conditions was swamped by the QPP at the stimulus onset. These results provide novel insights into the relationship between intrinsic patterns and stimulated brain activity.

Evaluating increases in sensitivity from NORDIC for diverse fMRI acquisition strategies

As the neuroimaging field moves towards detecting smaller effects at higher spatial resolutions, and faster sampling rates, there is increased attention given to the deleterious contribution of unstructured, thermal noise. Here, we critically evaluate the performance of a recently developed reconstruction method, termed NORDIC, for suppressing thermal noise using datasets acquired with various field strengths, voxel sizes, sampling rates, and task designs. Following minimal preprocessing, statistical activation (t-values) of NORDIC processed data was compared to the results obtained with alternative denoising methods. Additionally, we examined the consistency of the estimates of task responses at the single-voxel, single run level, using a finite impulse response (FIR) model. To examine the potential impact on effective image resolution, the overall smoothness of the data processed with different methods was estimated. Finally, to determine if NORDIC alters or removes temporal information important for modeling responses, we employed an exhaustive leave-p-out cross validation approach, using FIR task responses to predict held out timeseries, quantified using R². After NORDIC, the t-values are increased, an improvement comparable to what could be achieved by 1.5 voxels smoothing, and task events are clearly visible and have less cross-run error. These advantages are achieved with smoothness estimates increasing by less than 4%, while 1.5 voxel smoothing is associated with increases of over 140%. Cross-validated R²s based on the FIR models show that NORDIC is not measurably distorting the temporal structure of the data under this approach and is the best predictor of non-denoised time courses. The results demonstrate that analyzing 1 run of data after NORDIC produces results equivalent to using 2 to 3 original runs and that NORDIC performs equally well across a diverse array of functional imaging protocols. Significance Statement: For functional neuroimaging, the increasing availability of higher field strengths and ever higher spatiotemporal resolutions has led to concomitant increase in concerns about the deleterious effects of thermal noise. Historically this noise source was suppressed using methods that reduce spatial precision such as image blurring or averaging over a large number of trials or sessions, which necessitates large data collection efforts. Here, we critically evaluate the performance of a recently developed reconstruction method, termed NORDIC, which suppresses thermal noise. Across datasets varying in field strength, voxel sizes, sampling rates, and task designs, NORDIC produces substantial gains in data quality. Both conventional t-statistics derived from general linear models and coefficients of determination for predicting unseen data are improved. These gains match or even exceed those associated with 1 voxel Full Width Half Max image smoothing, however, even such small amounts of smoothing are associated with a 52% reduction in estimates of spatial precision, whereas the measurable difference in spatial precision is less than 4% following NORDIC.

Improving the study of brain-behavior relationships by revisiting basic assumptions

Neuroimaging research has been at the forefront of concerns regarding the failure of experimental findings to replicate. In the study of brain-behavior relationships, past failures to find replicable and robust effects have been attributed to methodological shortcomings. Methodological rigor is important, but there are other overlooked possibilities: most published studies share three foundational assumptions, often implicitly, that may be faulty. In this paper, we consider the empirical evidence from human brain imaging and the study of non-human animals that calls each foundational assumption into question. We then consider the opportunities for a robust science of brain-behavior relationships that await if scientists ground their research efforts in revised assumptions supported by current empirical evidence.

Vigilance associates with the low-dimensional structure of fMRI data

The human brain exhibits rich dynamics that reflect ongoing functional states. Patterns in fMRI data, detected in a data-driven manner, have uncovered recurring configurations that relate to individual and group differences in behavioral, cognitive, and clinical traits. However, resolving the neural and physiological processes that underlie such measurements is challenging, particularly without external measurements of brain state. A growing body of work points to underlying changes in vigilance as one driver of time-windowed fMRI connectivity states, calculated on the order of tens of seconds. Here we examine the degree to which the low-dimensional spatial structure of instantaneous fMRI activity is associated with vigilance levels, by testing whether vigilance-state detection can be carried out in an unsupervised manner based on individual BOLD time frames. To investigate this question, we first reduce the spatial dimensionality of fMRI data, and apply Gaussian Mixture Modeling to cluster the resulting low-dimensional data without any a priori vigilance information. Our analysis includes long-duration task and resting-state scans that are conducive to shifts in vigilance. We observe a close alignment between low-dimensional fMRI states (data-driven clusters) and measurements of vigilance derived from concurrent electroencephalography (EEG) and behavior. Whole-brain coactivation analysis revealed cortical anti-correlation patterns that resided primarily during higher behavioral- and EEG-defined levels of vigilance, while cortical activity was more often spatially uniform in states corresponding to lower vigilance. Overall, these findings indicate that vigilance states may be detected in the low-dimensional structure of fMRI data, even within individual time frames.

Whole-brain, gray and white matter time-locked functional signal changes with simple tasks and model-free analysis

Recent studies have revealed the production of time-locked blood oxygenation-level dependent (BOLD) functional MRI (fMRI) signals throughout the entire brain in response to a task, challenging the idea of sparse and localized brain functions, and highlighting the pervasiveness of potential false negative fMRI findings. In these studies, 'whole-brain' refers to gray matter regions only, which is the only tissue traditionally studied with fMRI. However, recent reports have also demonstrated reliable detection and analyses of BOLD signals in white matter which have been largely ignored in previous reports. Here, using model-free analysis and simple tasks, we investigate BOLD signal changes in both white and gray matters. We aimed to evaluate whether white matter also displays time-locked BOLD signals across all structural pathways in response to a stimulus. We find that both white and gray matter show time-locked activations across the whole-brain, with a majority of both tissue types showing statistically significant signal changes for all task stimuli investigated. We observed a wide range of signal responses to tasks, with different regions showing very different BOLD signal changes to the same task. Moreover, we find that each region may display different BOLD responses to different stimuli. Overall, we present compelling evidence that the whole brain, including both white and gray matter, show time-locked activation to multiple stimuli, not only challenging the idea of sparse functional localization, but also the prevailing wisdom of treating white matter BOLD signals as artefacts to be removed.

Intrinsic macroscale oscillatory modes driving long range functional connectivity in female rat brains detected by ultrafast fMRI

Spontaneous fluctuations in functional magnetic resonance imaging (fMRI) signals correlate across distant brain areas, shaping functionally relevant intrinsic networks. However, the generative mechanism of fMRI signal correlations, and in particular the link with locally-detected ultra-slow oscillations, are not fully understood. To investigate this link, we record ultrafast ultrahigh field fMRI signals (9.4 Tesla, temporal resolution = 38 milliseconds) from female rats across three anesthesia conditions. Power at frequencies extending up to 0.3 Hz is detected consistently across rat brains and is modulated by anesthesia level. Principal component analysis reveals a repertoire of modes, in which transient oscillations organize with fixed phase relationships across distinct cortical and subcortical structures. Oscillatory modes are found to vary between conditions, resonating at faster frequencies under medetomidine sedation and reducing both in number, frequency, and duration with the addition of isoflurane. Peaking in power within clear anatomical boundaries, these oscillatory modes point to an emergent systemic property. This work provides additional insight into the origin of oscillations detected in fMRI and the organizing principles underpinning spontaneous long-range functional connectivity.

A complex systems perspective on psychedelic brain action

Recent findings suggesting the potential transdiagnostic efficacy of psychedelic-assisted therapy have fostered the need to deepen our understanding of psychedelic brain action. Functional neuroimaging investigations have found that psychedelics reduce the functional segregation of large-scale brain networks. However, beyond this general trend, findings have been largely inconsistent. We argue here that a perspective based on complexity science that foregrounds the distributed, interactional, and dynamic nature of brain function may render these inconsistencies intelligible. We propose that psychedelics induce a mode of brain function that is more dynamically flexible, diverse, integrated, and tuned for information sharing, consistent with greater criticality. This 'meta' perspective has the potential to unify past findings and guide intuitions toward compelling mechanistic models.

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.

Whole-brain tracking of cocaine and sugar rewards processing

Natural rewards, such as food, and sex are appetitive stimuli available for animals in their natural environment. Similarly, addictive rewards such as drugs of abuse possess strong, positive valence, but their action relies on their pharmacological properties. Nevertheless, it is believed that both of these kinds of rewards activate similar brain circuitry. The present study aimed to discover which parts of the brain process the experience of natural and addictive rewards. To holistically address this question, we used a single-cell whole-brain imaging approach to find patterns of activation for acute and prolonged sucrose and cocaine exposure. We analyzed almost 400 brain structures and created a brain-wide map of specific, c-Fos-positive neurons engaged by these rewards. Acute but not prolonged sucrose exposure triggered a massive c-Fos expression throughout the brain. Cocaine exposure on the other hand potentiated c-Fos expression with prolonged use, engaging more structures than sucrose treatment. The functional connectivity analysis unraveled an increase in brain modularity after the initial exposure to both types of rewards. This modularity was increased after repeated cocaine, but not sucrose, intake. To check whether discrepancies between the processing of both types of rewards can be found on a cellular level, we further studied the nucleus accumbens, one of the most strongly activated brain structures by both sucrose and cocaine experience. We found a high overlap between natural and addictive rewards on the level of c-Fos expression. Electrophysiological measurements of cellular correlates of synaptic plasticity revealed that natural and addictive rewards alike induce the accumulation of silent synapses. These results strengthen the hypothesis that in the nucleus accumbens drugs of abuse cause maladaptive neuronal plasticity in the circuitry that typically processes natural rewards.

Towards functional spin-echo BOLD line-scanning in humans at 7T

OBJECTIVE

Neurons cluster into sub-millimeter spatial structures and neural activity occurs at millisecond resolutions; hence, ultimately, high spatial and high temporal resolutions are required for functional MRI. In this work, we implemented a spin-echo line-scanning (SELINE) sequence to use in high spatial and temporal resolution fMRI.

MATERIALS AND METHODS

A line is formed by simply rotating the spin-echo refocusing gradient to a plane perpendicular to the excited slice and by removing the phase-encoding gradient. This technique promises a combination of high spatial and temporal resolution (250 μm, 500 ms) and microvascular specificity of functional responses. We compared SELINE data to a corresponding gradient-echo version (GELINE).

RESULTS

We demonstrate that SELINE showed much-improved line selection (i.e. a sharper line profile) compared to GELINE, albeit at the cost of a significant drop in functional sensitivity.

DISCUSSION

This low functional sensitivity needs to be addressed before SELINE can be applied for neuroscientific purposes.

Beyond the Beta Rebound: Post-Task Responses in Oscillatory Activity follow Cessation of Working Memory Processes

Post-task responses (PTRs) are transitionary responses occurring for several seconds between the end of a stimulus/task and a period of rest. The most well-studied of these are beta band (13 - 30 Hz) PTRs in motor networks following movement, often called post-movement beta rebounds, which have been shown to differ in patients with schizophrenia and autism. Previous studies have proposed that beta PTRs reflect inhibition of task-positive networks to enable a return to resting brain activity, scaling with cognitive demand and reflecting cortical self-regulation. It is unknown whether PTRs are a phenomenon of the motor system, or whether they are a more general self-modulatory property of cortex that occur following cessation of higher cognitive processes as well as movement. To test this, we recorded magnetoencephalography (MEG) responses in 20 healthy participants to a working-memory task, known to recruit cortical networks associated with higher cognition. Our results revealed PTRs in the theta, alpha and beta bands across many regions of the brain, including the dorsal attention network (DAN) and lateral visual regions. These PTRs increased significantly (p < 0.05) in magnitude with working-memory load, an effect which is independent of oscillatory modulations occurring over the task period as well as those following individual stimuli. Furthermore, we showed that PTRs are functionally related to reaction times in left lateral visual (p < 0.05) and left parietal (p < 0.1) regions, while the oscillatory responses measured during the task period are not. Importantly, motor PTRs following button presses did not modulate with task condition, suggesting that PTRs in different networks are driven by different aspects of cognition. Our findings show that PTRs are not limited to motor networks but are widespread in regions which are recruited during the task. We provide evidence that PTRs have unique properties, scaling with cognitive load and correlating significantly with behaviour. Based on the evidence, we suggest that PTRs inhibit task-positive network activity to enable a transition to rest, however, further investigation is required to uncover their role in neuroscience and pathology.

Neural resources shift under Methylphenidate: A computational approach to examine anxiety-cognition interplay

The reciprocal interplay between anxiety and cognition is well documented. Anxiety negatively impacts cognition, while cognitive engagement can down-regulate anxiety. The brain mechanisms and dynamics underlying such interplay are not fully understood. To study this question, we experimentally and orthogonally manipulated anxiety (using a threat of shock paradigm) and cognition (using methylphenidate; MPH). The effects of these manipulations on the brain and behavior were evaluated in 50 healthy participants (25 MPH, 25 placebo), using an n-back working memory fMRI task (with low and high load conditions). Behaviorally, improved response accuracy was observed as a main effect of the drug across all conditions. We employed two approaches to understand the neural mechanisms underlying MPH-based cognitive enhancement in safe and threat conditions. First, we performed a hypothesis-driven computational analysis using a mathematical framework to examine how MPH putatively affects cognitive enhancement in the face of induced anxiety across two levels of cognitive load. Second, we performed an exploratory data analysis using Topological Data Analysis (TDA)-based Mapper to examine changes in spatiotemporal brain activity across the entire cortex. Both approaches provided converging evidence that MPH facilitated greater differential engagement of neural resources (brain activity) across low and high working memory load conditions. Furthermore, load-based differential management of neural resources reflects enhanced efficiency that is most powerful during higher load and induced anxiety conditions. Overall, our results provide novel insights regarding brain mechanisms that facilitate cognitive enhancement under MPH and, in future research, may be used to help mitigate anxiety-related cognitive underperformance.

Human visual consciousness involves large scale cortical and subcortical networks independent of task report and eye movement activity

The full neural circuits of conscious perception remain unknown. Using a visual perception task, we directly recorded a subcortical thalamic awareness potential (TAP). We also developed a unique paradigm to classify perceived versus not perceived stimuli using eye measurements to remove confounding signals related to reporting on conscious experiences. Using fMRI, we discovered three major brain networks driving conscious visual perception independent of report: first, increases in signal detection regions in visual, fusiform cortex, and frontal eye fields; and in arousal/salience networks involving midbrain, thalamus, nucleus accumbens, anterior cingulate, and anterior insula; second, increases in frontoparietal attention and executive control networks and in the cerebellum; finally, decreases in the default mode network. These results were largely maintained after excluding eye movement-based fMRI changes. Our findings provide evidence that the neurophysiology of consciousness is complex even without overt report, involving multiple cortical and subcortical networks overlapping in space and time.

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.

The global neuronal workspace as a broadcasting network

A new strategy for moving forward in the characterization of the global neuronal workspace (GNW) is proposed. According to Dehaene, Changeux, and colleagues (Dehaene, 2014, pp. 304, 312; Dehaene & Changeux, 2004, 2005), broadcasting is the main function of the GNW. However, the dynamic network properties described by recent graph theoretic GNW models are consistent with many large-scale communication processes that are different from broadcasting. We propose to apply a different graph theoretic approach, originally developed for optimizing information dissemination in communication networks, which can be used to identify the pattern of frequency and phase-specific directed functional connections that the GNW would exhibit only if it were a broadcasting network.

Allostasis as a core feature of hierarchical gradients in the human brain

This paper integrates emerging evidence from two broad streams of scientific literature into one common framework: (a) hierarchical gradients of functional connectivity that reflect the brain's large-scale structural architecture (e.g., a lamination gradient in the cerebral cortex); and (b) approaches to predictive processing and one of its specific instantiations called allostasis (i.e., the predictive regulation of energetic resources in the service of coordinating the body's internal systems). This synthesis begins to sketch a coherent, neurobiologically inspired framework suggesting that predictive energy regulation is at the core of human brain function, and by extension, psychological and behavioral phenomena, providing a shared vocabulary for theory building and knowledge accumulation.

Detection of functional activity in brain white matter using fiber architecture informed synchrony mapping

A general linear model is widely used for analyzing fMRI data, in which the blood oxygenation-level dependent (BOLD) signals in gray matter (GM) evoked in response to neural stimulation are modeled by convolving the time course of the expected neural activity with a canonical hemodynamic response function (HRF) obtained a priori. The maps of brain activity produced reflect the magnitude of local BOLD responses. However, detecting BOLD signals in white matter (WM) is more challenging as the BOLD signals are weaker and the HRF is different, and may vary more across the brain. Here we propose a model-free approach to detect changes in BOLD signals in WM by measuring task-evoked increases of BOLD signal synchrony in WM fibers. The proposed approach relies on a simple assumption that, in response to a functional task, BOLD signals in relevant fibers are modulated by stimulus-evoked neural activity and thereby show greater synchrony than when measured in a resting state, even if their magnitudes do not change substantially. This approach is implemented in two technical stages. First, for each voxel a fiber-architecture-informed spatial window is created with orientation distribution functions constructed from diffusion imaging data. This provides the basis for defining neighborhoods in WM that share similar local fiber architectures. Second, a modified principal component analysis (PCA) is used to estimate the synchrony of BOLD signals in each spatial window. The proposed approach is validated using a 3T fMRI dataset from the Human Connectome Project (HCP) at a group level. The results demonstrate that neural activity can be reliably detected as increases in fMRI signal synchrony within WM fibers that are engaged in a task with high sensitivities and reproducibility.

Improving power in functional magnetic resonance imaging by moving beyond cluster-level inference

Inference in neuroimaging typically occurs at the level of focal brain areas or circuits. Yet, increasingly, well-powered studies paint a much richer picture of broad-scale effects distributed throughout the brain, suggesting that many focal reports may only reflect the tip of the iceberg of underlying effects. How focal versus broad-scale perspectives influence the inferences we make has not yet been comprehensively evaluated using real data. Here, we compare sensitivity and specificity across procedures representing multiple levels of inference using an empirical benchmarking procedure that resamples task-based connectomes from the Human Connectome Project dataset (∼1,000 subjects, 7 tasks, 3 resampling group sizes, 7 inferential procedures). Only broad-scale (network and whole brain) procedures obtained the traditional 80% statistical power level to detect an average effect, reflecting >20% more statistical power than focal (edge and cluster) procedures. Power also increased substantially for false discovery rate- compared with familywise error rate-controlling procedures. The downsides are fairly limited; the loss in specificity for broad-scale and FDR procedures was relatively modest compared to the gains in power. Furthermore, the broad-scale methods we introduce are simple, fast, and easy to use, providing a straightforward starting point for researchers. This also points to the promise of more sophisticated broad-scale methods for not only functional connectivity but also related fields, including task-based activation. Altogether, this work demonstrates that shifting the scale of inference and choosing FDR control are both immediately attainable and can help remedy the issues with statistical power plaguing typical studies in the field.

Reappraisal-related downregulation of amygdala BOLD activation occurs only during the late trial window

During cognitive reappraisal, an individual reinterprets the meaning of an emotional stimulus to regulate the intensity of their emotional response. Prefrontal cortex activity has been found to support reappraisal and is putatively thought to downregulate the amygdala response to these stimuli. The timing of these regulation-related responses during the course of a trial, however, remains poorly understood. In the current fMRI study, participants were instructed to view or reappraise negative images and then rate how negative they felt following each image. The hemodynamic response function was estimated in 11 regions of interest for the entire time course of the trial including image viewing and rating. Notably, within the amygdala there was no evidence of downregulation in the early (picture viewing) window of the trial, only in the late (rating) window, which also correlated with a behavioral measure of reappraisal success. With respect to the prefrontal regions, some (e.g., inferior frontal gyrus) showed reappraisal-related activation in the early window, whereas others (e.g., middle frontal gyrus) showed increased activation primarily in the late window. These results highlight the temporal dynamics of different brain regions during emotion regulation and suggest that the amygdala response to negative images need not be immediately dampened to achieve successful cognitive reappraisal.

Evidence for non-selective response inhibition in uncertain contexts revealed by combined meta-analysis and Bayesian analysis of fMRI data

Response inhibition is typically considered a brain mechanism selectively triggered by particular “inhibitory” stimuli or events. Based on recent research, an alternative non-selective mechanism was proposed by several authors. Presumably, the inhibitory brain activity may be triggered not only by the presentation of “inhibitory” stimuli but also by any imperative stimuli, including Go stimuli, when the context is uncertain. Earlier support for this notion was mainly based on the absence of a significant difference between neural activity evoked by equiprobable Go and NoGo stimuli. Equiprobable Go/NoGo design with a simple response time task limits potential confounds between response inhibition and accompanying cognitive processes while not preventing prepotent automaticity. However, previous neuroimaging studies used classical null hypothesis significance testing, making it impossible to accept the null hypothesis. Therefore, the current research aimed to provide evidence for the practical equivalence of neuronal activity in the Go and NoGo trials using Bayesian analysis of functional magnetic resonance imaging (fMRI) data. Thirty-four healthy participants performed a cued Go/NoGo task with an equiprobable presentation of Go and NoGo stimuli. To independently localize brain areas associated with response inhibition in similar experimental conditions, we performed a meta-analysis of fMRI studies using equal-probability Go/NoGo tasks. As a result, we observed overlap between response inhibition areas and areas that demonstrate the practical equivalence of neuronal activity located in the right dorsolateral prefrontal cortex, parietal cortex, premotor cortex, and left inferior frontal gyrus. Thus, obtained results favour the existence of non-selective response inhibition, which can act in settings of contextual uncertainty induced by the equal probability of Go and NoGo stimuli.

Video-evoked fMRI BOLD responses are highly consistent across different data acquisition sites

Naturalistic imaging paradigms, in which participants view complex videos in the scanner, are increasingly used in human cognitive neuroscience. Videos evoke temporally synchronized brain responses that are similar across subjects as well as within subjects, but the reproducibility of these brain responses across different data acquisition sites has not yet been quantified. Here, we characterize the consistency of brain responses across independent samples of participants viewing the same videos in functional magnetic resonance imaging (fMRI) scanners at different sites (Indiana University and Caltech). We compared brain responses collected at these different sites for two carefully matched datasets with identical scanner models, acquisition, and preprocessing details, along with a third unmatched dataset in which these details varied. Our overall conclusion is that for matched and unmatched datasets alike, video-evoked brain responses have high consistency across these different sites, both when compared across groups and across pairs of individuals. As one might expect, differences between sites were larger for unmatched datasets than matched datasets. Residual differences between datasets could in part reflect participant-level variability rather than scanner- or data- related effects. Altogether our results indicate promise for the development and, critically, generalization of video fMRI studies of individual differences in healthy and clinical populations alike.

A massive 7T fMRI dataset to bridge cognitive neuroscience and artificial intelligence

Extensive sampling of neural activity during rich cognitive phenomena is critical for robust understanding of brain function. Here we present the Natural Scenes Dataset (NSD), in which high-resolution functional magnetic resonance imaging responses to tens of thousands of richly annotated natural scenes were measured while participants performed a continuous recognition task. To optimize data quality, we developed and applied novel estimation and denoising techniques. Simple visual inspections of the NSD data reveal clear representational transformations along the ventral visual pathway. Further exemplifying the inferential power of the dataset, we used NSD to build and train deep neural network models that predict brain activity more accurately than state-of-the-art models from computer vision. NSD also includes substantial resting-state and diffusion data, enabling network neuroscience perspectives to constrain and enhance models of perception and memory. Given its unprecedented scale, quality and breadth, NSD opens new avenues of inquiry in cognitive neuroscience and artificial intelligence.

Providing Evidence for the Null Hypothesis in Functional Magnetic Resonance Imaging Using Group-Level Bayesian Inference

Classical null hypothesis significance testing is limited to the rejection of the point-null hypothesis; it does not allow the interpretation of non-significant results. This leads to a bias against the null hypothesis. Herein, we discuss statistical approaches to ‘null effect’ assessment focusing on the Bayesian parameter inference (BPI). Although Bayesian methods have been theoretically elaborated and implemented in common neuroimaging software packages, they are not widely used for ‘null effect’ assessment. BPI considers the posterior probability of finding the effect within or outside the region of practical equivalence to the null value. It can be used to find both ‘activated/deactivated’ and ‘not activated’ voxels or to indicate that the obtained data are not sufficient using a single decision rule. It also allows to evaluate the data as the sample size increases and decide to stop the experiment if the obtained data are sufficient to make a confident inference. To demonstrate the advantages of using BPI for fMRI data group analysis, we compare it with classical null hypothesis significance testing on empirical data. We also use simulated data to show how BPI performs under different effect sizes, noise levels, noise distributions and sample sizes. Finally, we consider the problem of defining the region of practical equivalence for BPI and discuss possible applications of BPI in fMRI studies. To facilitate ‘null effect’ assessment for fMRI practitioners, we provide Statistical Parametric Mapping 12 based toolbox for Bayesian inference.

Hyperbolic trade-off: The importance of balancing trial and subject sample sizes in neuroimaging

Here we investigate the crucial role of trials in task-based neuroimaging from the perspectives of statistical efficiency and condition-level generalizability. Big data initiatives have gained popularity for leveraging a large sample of subjects to study a wide range of effect magnitudes in the brain. On the other hand, most task-based FMRI designs feature a relatively small number of subjects, so that resulting parameter estimates may be associated with compromised precision. Nevertheless, little attention has been given to another important dimension of experimental design, which can equally boost a study's statistical efficiency: the trial sample size. The common practice of condition-level modeling implicitly assumes no cross-trial variability. Here, we systematically explore the different factors that impact effect uncertainty, drawing on evidence from hierarchical modeling, simulations and an FMRI dataset of 42 subjects who completed a large number of trials of cognitive control task. We find that, due to an approximately symmetic hyperbola-relationship between trial and subject sample sizes in the presence of relatively large cross-trial variability, 1) trial sample size has nearly the same impact as subject sample size on statistical efficiency; 2) increasing both the number of trials and subjects improves statistical efficiency more effectively than focusing on subjects alone; 3) trial sample size can be leveraged alongside subject sample size to improve the cost-effectiveness of an experimental design; 4) for small trial sample sizes, trial-level modeling, rather than condition-level modeling through summary statistics, may be necessary to accurately assess the standard error of an effect estimate. We close by making practical suggestions for improving experimental designs across neuroimaging and behavioral studies.

Is it time to put rest to rest?

The so-called resting state, in which participants lie quietly with no particular inputs or outputs, represented a paradigm shift from conventional task-based studies in human neuroimaging. Our foray into rest was fruitful from both a scientific and methodological perspective, but at this point, how much more can we learn from rest on its own? While rest still dominates in many subfields, data from tasks have empirically demonstrated benefits, as well as the potential to provide insights about the mind in addition to the brain. I argue that we can accelerate progress in human neuroscience by de-emphasizing rest in favor of more grounded experiments, including promising integrated designs that respect the prominence of self-generated activity while offering enhanced control and interpretability.

A Quantitative Data-Driven Analysis Framework for Resting-State Functional Magnetic Resonance Imaging: A Study of the Impact of Adult Age

The objective of this study is to introduce a new quantitative data-driven analysis (QDA) framework for the analysis of resting-state fMRI (R-fMRI) and use it to investigate the effect of adult age on resting-state functional connectivity (RFC). Whole-brain R-fMRI measurements were conducted on a 3T clinical MRI scanner in 227 healthy adult volunteers (N = 227, aged 18-76 years old, male/female = 99/128). With the proposed QDA framework we derived two types of voxel-wise RFC metrics: the connectivity strength index and connectivity density index utilizing the convolutions of the cross-correlation histogram with different kernels. Furthermore, we assessed the negative and positive portions of these metrics separately. With the QDA framework we found age-related declines of RFC metrics in the superior and middle frontal gyri, posterior cingulate cortex (PCC), right insula and inferior parietal lobule of the default mode network (DMN), which resembles previously reported results using other types of RFC data processing methods. Importantly, our new findings complement previously undocumented results in the following aspects: (1) the PCC and right insula are anti-correlated and tend to manifest simultaneously declines of both the negative and positive connectivity strength with subjects' age; (2) separate assessment of the negative and positive RFC metrics provides enhanced sensitivity to the aging effect; and (3) the sensorimotor network depicts enhanced negative connectivity strength with the adult age. The proposed QDA framework can produce threshold-free and voxel-wise RFC metrics from R-fMRI data. The detected adult age effect is largely consistent with previously reported studies using different R-fMRI analysis approaches. Moreover, the separate assessment of the negative and positive contributions to the RFC metrics can enhance the RFC sensitivity and clarify some of the mixed results in the literature regarding to the DMN and sensorimotor network involvement in adult aging.

Frontal and parietal background connectivity and their dynamic changes account for individual differences in the multisensory representation of peripersonal space

Functional connectivity (FC) of brain networks dynamically fluctuates during both rest and task execution. Individual differences in dynamic FC have been associated with several cognitive and behavioral traits. However, whether dynamic FC also contributes to sensorimotor representations guiding body-environment interactions, such as the representation of peripersonal space (PPS), is currently unknown. PPS is the space immediately surrounding the body and acts as a multisensory interface between the individual and the environment. We used an audio-tactile task with approaching sounds to map the individual PPS extension, and fMRI to estimate the background FC. Specifically, we analyzed FC values for each stimulus type (near and far space) and its across-trial variability. FC was evaluated between task-relevant nodes of two fronto-parietal networks (the Dorsal Attention Network, DAN, and the Fronto-Parietal Network, FPN) and a key PPS region in the premotor cortex (PM). PM was significantly connected to specific task-relevant nodes of the DAN and the FPN during the audio-tactile task, and FC was stronger while processing near space, as compared to far space. At the individual level, less PPS extension was associated with stronger premotor-parietal FC during processing of near space, while the across-trial variability of premotor-parietal and premotor-frontal FC was higher during the processing of far space. Notably, only across-trial FC variability captured the near-far modulation of space processing. Our findings indicate that PM connectivity with task-relevant frontal and parietal regions and its dynamic changes participate in the mechanisms that enable PPS representation, in agreement with the idea that neural variability plays a crucial role in plastic and dynamic sensorimotor representations.

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.

Voxel-Wise Linearity Analysis of Increments and Decrements in BOLD Responses in Human Visual Cortex Using a Contrast Adaptation Paradigm

The linearity of BOLD responses is a fundamental presumption in most analysis procedures for BOLD fMRI studies. Previous studies have examined the linearity of BOLD signal increments, but less is known about the linearity of BOLD signal decrements. The present study assessed the linearity of both BOLD signal increments and decrements in the human primary visual cortex using a contrast adaptation paradigm. Results showed that both BOLD signal increments and decrements kept linearity to long stimuli (e.g., 3 s, 6 s), yet, deviated from linearity to transient stimuli (e.g., 1 s). Furthermore, a voxel-wise analysis showed that the deviation patterns were different for BOLD signal increments and decrements: while the BOLD signal increments demonstrated a consistent overestimation pattern, the patterns for BOLD signal decrements varied from overestimation to underestimation. Our results suggested that corrections to deviations from linearity of transient responses should consider the different effects of BOLD signal increments and decrements.

Imaging faster neural dynamics with fast fMRI: A need for updated models of the hemodynamic response

Fast fMRI enables the detection of neural dynamics over timescales of hundreds of milliseconds, suggesting it may provide a new avenue for studying subsecond neural processes in the human brain. The magnitudes of these fast fMRI dynamics are far greater than predicted by canonical models of the hemodynamic response. Several studies have established nonlinear properties of the hemodynamic response that have significant implications for fast fMRI. We first review nonlinear properties of the hemodynamic response function that may underlie fast fMRI signals. We then illustrate the breakdown of canonical hemodynamic response models in the context of fast neural dynamics. We will then argue that the canonical hemodynamic response function is not likely to reflect the BOLD response to neuronal activity driven by sparse or naturalistic stimuli or perhaps to spontaneous neuronal fluctuations in the resting state. These properties suggest that fast fMRI is capable of tracking surprisingly fast neuronal dynamics, and we discuss the neuroscientific questions that could be addressed using this approach.