Functional connectivity (FC) patterns in functional MRI exhibit dynamic behavior on the scale of seconds, with rich spatiotemporal structure and limited sets of whole-brain, quasi-stable FC ...configurations (FC states) recurring across time and subjects. Based on previous evidence linking various aspects of cognition to group-level, minute-to-minute FC changes in localized connections, we hypothesized that whole-brain FC states may reflect the global, orchestrated dynamics of cognitive processing on the scale of seconds. To test this hypothesis, subjects were continuously scanned as they engaged in and transitioned between mental states dictated by tasks. FC states computed within windows as short as 22.5 s permitted robust tracking of cognition in single subjects with near perfect accuracy. Accuracy dropped markedly for subjects with the lowest task performance. Spatially restricting FC information decreased accuracy at short time scales, emphasizing the distributed nature of whole-brain FC dynamics, beyond univariate magnitude changes, as valuable markers of cognition.
Layer-dependent fMRI allows measurements of information flow in cortical circuits, as afferent and efferent connections terminate in different cortical layers. However, it is unknown to what level ...human fMRI is specific and sensitive enough to reveal directional functional activity across layers. To answer this question, we developed acquisition and analysis methods for blood-oxygen-level-dependent (BOLD) and cerebral-blood-volume (CBV)-based laminar fMRI and used these to discriminate four different tasks in the human motor cortex (M1). In agreement with anatomical data from animal studies, we found evidence for somatosensory and premotor input in superficial layers of M1 and for cortico-spinal motor output in deep layers. Laminar resting-state fMRI showed directional functional connectivity of M1 with somatosensory and premotor areas. Our findings demonstrate that CBV-fMRI can be used to investigate cortical activity in humans with unprecedented detail, allowing investigations of information flow between brain regions and outperforming conventional BOLD results that are often buried under vascular biases.
•A CBV-fMRI method is developed for sampling brain activity across cortical layers•We can separately measure input and output activity of primary motor cortex•CBV-fMRI can measure laminar activity with higher accuracy than BOLD-fMRI•Laminar resting-state fMRI can reveal directional connectivity between brain areas
Huber et al. demonstrate an MRI method to measure brain activity changes at the spatial resolution of cortical layers in humans. This allows investigations of directional functional connectivity, paving the way for non-invasive studies investigating information flow between brain regions.
A central challenge in the fMRI based study of functional connectivity is distinguishing neuronally related signal fluctuations from the effects of motion, physiology, and other nuisance sources. ...Conventional techniques for removing nuisance effects include modeling of noise time courses based on external measurements followed by temporal filtering. These techniques have limited effectiveness. Previous studies have shown using multi-echo fMRI that neuronally related fluctuations are Blood Oxygen Level Dependent (BOLD) signals that can be characterized in terms of changes in R2* and initial signal intensity (S0) based on the analysis of echo-time (TE) dependence. We hypothesized that if TE-dependence could be used to differentiate BOLD and non-BOLD signals, non-BOLD signal could be removed to denoise data without conventional noise modeling. To test this hypothesis, whole brain multi-echo data were acquired at 3 TEs and decomposed with Independent Components Analysis (ICA) after spatially concatenating data across space and TE. Components were analyzed for the degree to which their signal changes fit models for R2* and S0 change, and summary scores were developed to characterize each component as BOLD-like or not BOLD-like. These scores clearly differentiated BOLD-like “functional network” components from non BOLD-like components related to motion, pulsatility, and other nuisance effects. Using non BOLD-like component time courses as noise regressors dramatically improved seed-based correlation mapping by reducing the effects of high and low frequency non-BOLD fluctuations. A comparison with seed-based correlation mapping using conventional noise regressors demonstrated the superiority of the proposed technique for both individual and group level seed-based connectivity analysis, especially in mapping subcortical–cortical connectivity. The differentiation of BOLD and non-BOLD components based on TE-dependence was highly robust, which allowed for the identification of BOLD-like components and the removal of non BOLD-like components to be implemented as a fully automated procedure.
Display omitted
► ICA components were analyzed for TE-dependence. ► Component-level TE-dependence scores identified BOLD components. ► BOLD networks were identified without spatial templates or frequency thresholds. ► Removing non-BOLD components dramatically denoised data. ► Robust subcortical-cortical connectivity was revealed after denoising.
Findings from single-cell recording studies suggest that a comparison of the outputs of different pools of selectively tuned lower-level sensory neurons may be a general mechanism by which ...higher-level brain regions compute perceptual decisions. For example, when monkeys must decide whether a noisy field of dots is moving upward or downward, a decision can be formed by computing the difference in responses between lower-level neurons sensitive to upward motion and those sensitive to downward motion. Here we use functional magnetic resonance imaging and a categorization task in which subjects decide whether an image presented is a face or a house to test whether a similar mechanism is also at work for more complex decisions in the human brain and, if so, where in the brain this computation might be performed. Activity within the left dorsolateral prefrontal cortex is greater during easy decisions than during difficult decisions, covaries with the difference signal between face- and house-selective regions in the ventral temporal cortex, and predicts behavioural performance in the categorization task. These findings show that even for complex object categories, the comparison of the outputs of different pools of selectively tuned neurons could be a general mechanism by which the human brain computes perceptual decisions.
Head motion can be estimated at any point of fMRI image processing. Processing steps involving temporal interpolation (e.g., slice time correction or outlier replacement) often precede motion ...estimation in the literature. From first principles it can be anticipated that temporal interpolation will alter head motion in a scan. Here we demonstrate this effect and its consequences in five large fMRI datasets. Estimated head motion was reduced by 10-50% or more following temporal interpolation, and reductions were often visible to the naked eye. Such reductions make the data seem to be of improved quality. Such reductions also degrade the sensitivity of analyses aimed at detecting motion-related artifact and can cause a dataset with artifact to falsely appear artifact-free. These reduced motion estimates will be particularly problematic for studies needing estimates of motion in time, such as studies of dynamics. Based on these findings, it is sensible to obtain motion estimates prior to any image processing (regardless of subsequent processing steps and the actual timing of motion correction procedures, which need not be changed). We also find that outlier replacement procedures change signals almost entirely during times of motion and therefore have notable similarities to motion-targeting censoring strategies (which withhold or replace signals entirely during times of motion).
Brain functional connectivity (FC) changes have been measured across seconds using fMRI. This is true for both rest and task scenarios. Moreover, it is well accepted that task engagement alters FC, ...and that dynamic estimates of FC during and before task events can help predict their nature and performance. Yet, when it comes to dynamic FC (dFC) during rest, there is no consensus about its origin or significance. Some argue that rest dFC reflects fluctuations in on-going cognition, or is a manifestation of intrinsic brain maintenance mechanisms, which could have predictive clinical value. Conversely, others have concluded that rest dFC is mostly the result of sampling variability, head motion or fluctuating sleep states. Here, we present novel analyses suggesting that rest dFC is influenced by short periods of spontaneous cognitive-task-like processes, and that the cognitive nature of such mental processes can be inferred blindly from the data. As such, several different behaviorally relevant whole-brain FC configurations may occur during a single rest scan even when subjects were continuously awake and displayed minimal motion. In addition, using low dimensional embeddings as visualization aids, we show how FC states—commonly used to summarize and interpret resting dFC—can accurately and robustly reveal periods of externally imposed tasks; however, they may be less effective in capturing periods of distinct cognition during rest.
•Covert self-driven cognition is a key contributor to dynamic FC at rest.•Several behaviorally-relevant FC configurations occur during a single rest scan.•Propose a way to map dynamic FC onto cognitive states in an open-ended fashion.•Inferred mental states agree with prior reports of common resting mental activities.
Inferior temporal (IT) object representations have been intensively studied in monkeys and humans, but representations of the same particular objects have never been compared between the species. ...Moreover, IT's role in categorization is not well understood. Here, we presented monkeys and humans with the same images of real-world objects and measured the IT response pattern elicited by each image. In order to relate the representations between the species and to computational models, we compare response-pattern dissimilarity matrices. IT response patterns form category clusters, which match between man and monkey. The clusters correspond to animate and inanimate objects; within the animate objects, faces and bodies form subclusters. Within each category, IT distinguishes individual exemplars, and the within-category exemplar similarities also match between the species. Our findings suggest that primate IT across species may host a common code, which combines a categorical and a continuous representation of objects.
Recent advances in MRI receiver and coil technologies have significantly improved image signal-to-noise ratios (SNR) and thus temporal SNR (TSNR). These gains in SNR and TSNR have allowed the ...detection of fMRI signal changes at higher spatial resolution and therefore have increased the potential to localize small brain structures such as cortical layers and columns. The majority of current fMRI processing strategies employ multi-subject averaging and therefore require spatial smoothing and normalization, effectively negating these gains in spatial resolution higher than about 10 mm
3. Reliable detection of activation in single subjects at high resolution is becoming a more common desire among fMRI researchers who are interested in comparing individuals rather than populations. Since TSNR decreases with voxel volume, detection of activation at higher resolutions requires longer scan durations. The relationship between TSNR, voxel volume and detectability is highly non-linear. In this study, the relationship between TSNR and the necessary fMRI scan duration required to obtain significant results at varying
P values is determined both experimentally and theoretically. The results demonstrate that, with a TSNR of 50, detection of activation of above 2% requires at most 350 scan volumes (when steps are taken to remove the influence of physiological noise from the data). Importantly, these results also demonstrate that, for activation magnitude on the order of 1%, the scan duration required is more sensitive to the TSNR level than at 2%. This study showed that with voxel volumes of ∼
10 mm
3 at 3 T, and a corresponding TSNR of ∼
50, the required number of time points that guarantees detection of signal changes of 1% is about 860, but if TSNR increases by only 20%, the time for detection decreases by more than 30%. More than just being an exercise in numbers, these results imply that imaging of columnar resolution (effect size
=
1% and assuming a TR of 1 s) at 3 T will require either 10 min for a TSNR of 60 or 40 min for a TSNR of 30. The implication is that at these resolutions, TSNR is likely to be critical for determining success or failure of an experiment.
Functional connectivity analysis of resting state blood oxygen level–dependent (BOLD) functional MRI is widely used for noninvasively studying brain functional networks. Recent findings have ...indicated, however, that even small (≤1 mm) amounts of head movement during scanning can disproportionately bias connectivity estimates, despite various preprocessing efforts. Further complications for interregional connectivity estimation from time domain signals include the unaccounted reduction in BOLD degrees of freedom related to sensitivity losses from high subject motion. To address these issues, we describe an integrated strategy for data acquisition, denoising, and connectivity estimation. This strategy builds on our previously published technique combining data acquisition with multiecho (ME) echo planar imaging and analysis with spatial independent component analysis (ICA), called ME-ICA, which distinguishes BOLD (neuronal) and non-BOLD (artifactual) components based on linear echo-time dependence of signals—a characteristic property of BOLD $T_{2}^{*}$ signal changes. Here we show for 32 control subjects that this method provides a physically principled and nearly operator-independent way of removing complex artifacts such as motion from resting state data. We then describe a robust estimator of functional connectivity based on interregional correlation of BOLD-independent component coefficients. This estimator, called independent components regression, considerably simplifies statistical inference for functional connectivity because degrees of freedom equals the number of independent coefficients. Compared with traditional connectivity estimation methods, the proposed strategy results in fourfold improvements in signal-to-noise ratio, functional connectivity analysis with improved specificity, and valid statistical inference with nominal control of type 1 error in contrasts of connectivity between groups with different levels of subject motion.
Subtle changes in a subject's breathing rate or depth, which occur naturally during rest at low frequencies (<0.1 Hz), have been shown to be significantly correlated with fMRI signal changes ...throughout gray matter and near large vessels. The goal of this study was to investigate the impact of these low-frequency respiration variations on both task activation fMRI studies and resting-state functional connectivity analysis. Unlike MR signal changes correlated with the breathing motion (∼0.3 Hz), BOLD signal changes correlated with across-breath variations in respiratory volume (∼0.03 Hz) appear localized to blood vessels and regions with high blood volume, such as gray matter, similar to changes seen in response to a breath-hold challenge. In addition, the respiration-variation-induced signal changes were found to coincide with many of the areas identified as part of the ‘default mode’ network, a set of brain regions hypothesized to be more active at rest. Regions could therefore be classified as being part of a resting network based on their similar respiration-induced changes rather than their synchronized neuronal activity. Monitoring and removing these respiration variations led to a significant improvement in the identification of task-related activation and deactivation and only slight differences in regions correlated with the posterior cingulate at rest. Regressing out global signal changes or cueing the subject to breathe at a constant rate and depth resulted in an improved spatial overlap between deactivations and resting-state correlations among areas that showed deactivation.