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.

Ultra-slow fMRI fluctuations in the fourth ventricle as a marker of drowsiness.

Wakefulness levels modulate estimates of functional connectivity (FC), and, if unaccounted for, can become a substantial confound in resting-state fMRI. Unfortunately, wakefulness is rarely monitored due to the need for additional concurrent recordings (e.g., eye tracking, EEG). Recent work has shown that strong fluctuations around 0.05Hz, hypothesized to be CSF inflow, appear in the fourth ventricle (FV) when subjects fall asleep, and that they correlate significantly with the global signal. The analysis of these fluctuations could provide an easy way to evaluate wakefulness in fMRI-only data and improve our understanding of FC during sleep. Here we evaluate this possibility using the 7T resting-state sample from the Human Connectome Project (HCP). Our results replicate the observation that fourth ventricle ultra-slow fluctuations (∼0.05Hz) with inflow-like characteristics (decreasing in intensity for successive slices) are present in scans during which subjects did not comply with instructions to keep their eyes open (i.e., drowsy scans). This is true despite the HCP data not being optimized for the detection of inflow-like effects. In addition, time-locked BOLD fluctuations of the same frequency could be detected in large portions of grey matter with a wide range of temporal delays and contribute in significant ways to our understanding of how FC changes during sleep. First, these ultra-slow fluctuations explain half of the increase in global signal that occurs during descent into sleep. Similarly, global shifts in FC between awake and sleep states are driven by changes in this slow frequency band. Second, they can influence estimates of inter-regional FC. For example, disconnection between frontal and posterior components of the Defulat Mode Network (DMN) typically reported during sleep were only detectable after regression of these ultra-slow fluctuations. Finally, we report that the temporal evolution of the power spectrum of these ultra-slow FV fluctuations can help us reproduce sample-level sleep patterns (e.g., a substantial number of subjects descending into sleep 3 minutes following scanning onset), partially rank scans according to overall drowsiness levels, and predict individual segments of elevated drowsiness (at 60 seconds resolution) with 71% accuracy.

Layer-specific activation in human primary somatosensory cortex during tactile temporal prediction error processing.

The human brain continuously generates predictions of incoming sensory input and calculates corresponding prediction errors from the perceived inputs to update internal predictions. In human primary somatosensory cortex (area 3b), different cortical layers are involved in receiving the sensory input and generation of error signals. It remains unknown, however, how the layers in the human area 3b contribute to the temporal prediction error processing. To investigate prediction error representation in the area 3b across layers, we acquired layer-specific functional magnetic resonance imaging (fMRI) data at 7T from human area 3b during a task of index finger poking with no-delay, short-delay and long-delay touching sequences. We demonstrate that all three tasks increased activity in both superficial and deep layers of area 3b compared to the random sensory input. The fMRI signal was differentially modulated solely in the deep layers rather than the superficial layers of area 3b by the delay time. Compared with the no-delay stimuli, activity was greater in the deep layers of area 3b during the short-delay stimuli but lower during the long-delay stimuli. This difference activity features in the superficial and deep layers suggest distinct functional contributions of area 3b layers to tactile temporal prediction error processing. The functional segregation in area 3b across layers may reflect that the excitatory and inhibitory interplay in the sensory cortex contributions to flexible communication between cortical layers or between cortical areas.

Different activation signatures in the primary sensorimotor and higher-level regions for haptic three-dimensional curved surface exploration.

Haptic object perception begins with continuous exploratory contact, and the human brain needs to accumulate sensory information continuously over time. However, it is still unclear how the primary sensorimotor cortex (PSC) interacts with these higher-level regions during haptic exploration over time. This functional magnetic resonance imaging (fMRI) study investigates time-dependent haptic object processing by examining brain activity during haptic 3D curve and roughness estimations. For this experiment, we designed sixteen haptic stimuli (4 kinds of curves × 4 varieties of roughness) for the haptic curve and roughness estimation tasks. Twenty participants were asked to move their right index and middle fingers along the surface twice and to estimate one of the two features-roughness or curvature-depending on the task instruction. We found that the brain activity in several higher-level regions (e.g., the bilateral posterior parietal cortex) linearly increased as the number of curves increased during the haptic exploration phase. Surprisingly, we found that the contralateral PSC was parametrically modulated by the number of curves only during the late exploration phase but not during the early exploration phase. In contrast, we found no similar parametric modulation activity patterns during the haptic roughness estimation task in either the contralateral PSC or in higher-level regions. Thus, our findings suggest that haptic 3D object perception is processed across the cortical hierarchy, whereas the contralateral PSC interacts with other higher-level regions across time in a manner that is dependent upon the features of the object.

Idiosynchrony: From shared responses to individual differences during naturalistic neuroimaging.

Two ongoing movements in human cognitive neuroscience have researchers shifting focus from group-level inferences to characterizing single subjects, and complementing tightly controlled tasks with rich, dynamic paradigms such as movies and stories. Yet relatively little work combines these two, perhaps because traditional analysis approaches for naturalistic imaging data are geared toward detecting shared responses rather than between-subject variability. Here, we review recent work using naturalistic stimuli to study individual differences, and advance a framework for detecting structure in idiosyncratic patterns of brain activity, or "idiosynchrony". Specifically, we outline the emerging technique of inter-subject representational similarity analysis (IS-RSA), including its theoretical motivation and an empirical demonstration of how it recovers brain-behavior relationships during movie watching using data from the Human Connectome Project. We also consider how stimulus choice may affect the individual signal and discuss areas for future research. We argue that naturalistic neuroimaging paradigms have the potential to reveal meaningful individual differences above and beyond those observed during traditional tasks or at rest.

Sub-millimeter fMRI reveals multiple topographical digit representations that form action maps in human motor cortex.

The human brain coordinates a wide variety of motor activities. On a large scale, the cortical motor system is topographically organized such that neighboring body parts are represented by neighboring brain areas. This homunculus-like somatotopic organization along the central sulcus has been observed using neuroimaging for large body parts such as the face, hands and feet. However, on a finer scale, invasive electrical stimulation studies show deviations from this somatotopic organization that suggest an organizing principle based on motor actions rather than body part moved. It has not been clear how the action-map organization principle of the motor cortex in the mesoscopic (sub-millimeter) regime integrates into a body map organization principle on a macroscopic scale (cm). Here we developed and applied advanced mesoscopic (sub-millimeter) fMRI and analysis methodology to non-invasively investigate the functional organization topography across columnar and laminar structures in humans. Compared to previous methods, in this study, we could capture locally specific blood volume changes across entire brain regions along the cortical curvature. We find that individual fingers have multiple mirrored representations in the primary motor cortex depending on the movements they are involved in. We find that individual digits have cortical representations up to 3 â€‹mm apart from each other arranged in a column-like fashion. These representations are differentially engaged depending on whether the digits' muscles are used for different motor actions such as flexion movements, like grasping a ball or retraction movements like releasing a ball. This research provides a starting point for non-invasive investigation of mesoscale topography across layers and columns of the human cortex and bridges the gap between invasive electrophysiological investigations and large coverage non-invasive neuroimaging.

Imaging the spontaneous flow of thought: Distinct periods of cognition contribute to dynamic functional connectivity during rest.

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.

Visual temporal frequency preference shows a distinct cortical architecture using fMRI.

Studies of visual temporal frequency preference typically examine frequencies under 20 Hz and measure local activity to evaluate the sensitivity of different cortical areas to variations in temporal frequencies. Most of these studies have not attempted to map preferred temporal frequency within and across visual areas, nor have they explored in detail, stimuli at gamma frequency, which recent research suggests may have potential clinical utility. In this study, we address this gap by using functional magnetic resonance imaging (fMRI) to measure response to flickering visual stimuli varying in frequency from 1 to 40 Hz. We apply stimulation in both a block design to examine task response and a steady-state design to examine functional connectivity. We observed distinct activation patterns between 1 Hz and 40 Hz stimuli. We also found that the correlation between medial thalamus and visual cortex was modulated by the temporal frequency. The modulation functions and tuned frequencies are different for the visual activity and thalamo-visual correlations. Using both fMRI activity and connectivity measurements, we show evidence for a temporal frequency specific organization across the human visual system.

Efficacy of different dynamic functional connectivity methods to capture cognitively relevant information.

Given the dynamic nature of the human brain, there has been an increasing interest in investigating short-term temporal changes in functional connectivity, also known as dynamic functional connectivity (dFC), i.e., the time-varying inter-regional statistical dependence of blood oxygenation level-dependent (BOLD) signal within the constraints of a single scan. Numerous methodologies have been proposed to characterize dFC during rest and task, but few studies have compared them in terms of their efficacy to capture behavioral and clinically relevant dynamics. This is mostly due to lack of a well-defined ground truth, especially for rest scans. In this study, with a multitask dataset (rest, memory, video, and math) serving as ground truth, we investigated the efficacy of several dFC estimation techniques at capturing cognitively relevant dFC modulation induced by external tasks. We evaluated two framewise methods (dFC estimates for a single time point): dynamic conditional correlation (DCC) and jackknife correlation (JC); and five window-based methods: sliding window correlation (SWC), sliding window correlation with L1-regularization (SWC_L1), a combination of DCC and SWC called moving average DCC (DCC_MA), multiplication of temporal derivatives (MTD), and a variant of jackknife correlation called delete-d jackknife correlation (dJC). The efficacy is defined as each dFC metric's ability to successfully subdivide multitask scans into cognitively homogenous segments (even if those segments are not temporally continuous). We found that all window-based dFC methods performed well for commonly used window lengths (WL ≥ 30sec), with sliding window methods (SWC, SWC_L1) as well as the hybrid DCC_MA approach performing slightly better. For shorter window lengths (WL ≤ 15sec), DCC_MA and dJC produced the best results. Neither framewise method (i.e., DCC and JC) led to dFC estimates with high accuracy.

Techniques for blood volume fMRI with VASO: From low-resolution mapping towards sub-millimeter layer-dependent applications.

Quantitative cerebral blood volume (CBV) fMRI has the potential to overcome several specific limitations of BOLD fMRI. It provides direct physiological interpretability and promises superior localization specificity in applications of sub-millimeter resolution fMRI applications at ultra-high magnetic fields (7T and higher). Non-invasive CBV fMRI using VASO (vascular space occupancy), however, is inherently limited with respect to its data acquisition efficiency, restricting its imaging coverage and achievable spatial and temporal resolution. This limitation may be reduced with recent advanced acceleration and reconstruction strategies that allow two-dimensional acceleration, such as in simultaneous multi-slice (SMS) 2D-EPI or 3D-EPI in combination with CAIPIRINHA field-of-view shifting. In this study, we sought to determine the functional sensitivity and specificity of these readout strategies with VASO over a broad range of spatial resolutions; spanning from low spatial resolution (3mm) whole-cortex to sub-millimeter (0.75mm) slab-of-cortex (for cortical layer-dependent applications). In the thermal-noise-dominated regime of sub-millimeter resolutions, 3D-EPI-VASO provides higher temporal stability and sensitivity to detect changes in CBV compared to 2D-EPI-VASO. In this regime, 3D-EPI-VASO unveils task activation located in the cortical laminae with little contamination from surface veins, in contrast to the cortical surface weighting of GE-BOLD fMRI. In the physiological-noise-dominated regime of lower resolutions, however, 2D-SMS-VASO shows superior performance compared to 3D-EPI-VASO. Due to its superior sensitivity at a layer-dependent level, 3D-EPI VASO promises to play an important role in future neuroscientific applications of layer-dependent fMRI.

A functional connectivity-based neuromarker of sustained attention generalizes to predict recall in a reading task.

Sustaining attention to the task at hand is a crucial part of everyday life, from following a lecture at school to maintaining focus while driving. Lapses in sustained attention are frequent and often problematic, with conditions such as attention deficit hyperactivity disorder affecting millions of people worldwide. Recent work has had some success in finding signatures of sustained attention in whole-brain functional connectivity (FC) measures during basic tasks, but since FC can be dynamic and task-dependent, it remains unclear how fully these signatures would generalize to a more complex and naturalistic scenario. To this end, we used a previously defined whole-brain FC network - a marker of attention that was derived from a sustained attention task - to predict the ability of participants to recall material during a free-viewing reading task. Though the predictive network was trained on a different task and set of participants, the strength of FC in the sustained attention network predicted reading recall significantly better than permutation tests where behavior was scrambled to simulate chance performance. To test the generalization of the method used to derive the sustained attention network, we applied the same method to our reading task data to find a new FC network whose strength specifically predicts reading recall. Even though the sustained attention network provided significant prediction of recall, the reading network was more predictive of recall accuracy. The new reading network's spatial distribution indicates that reading recall is highest when temporal pole regions have higher FC with left occipital regions and lower FC with bilateral supramarginal gyrus. Right cerebellar to right frontal connectivity is also indicative of poor reading recall. We examine these and other differences between the two predictive FC networks, providing new insight into the task-dependent nature of FC-based performance metrics.

Tracking ongoing cognition in individuals using brief, whole-brain functional connectivity patterns.

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.

Evaluation of multi-echo ICA denoising for task based fMRI studies: Block designs, rapid event-related designs, and cardiac-gated fMRI.

Multi-echo fMRI, particularly the multi-echo independent component analysis (ME-ICA) algorithm, has previously proven useful for increasing the sensitivity and reducing false positives for functional MRI (fMRI) based resting state connectivity studies. Less is known about its efficacy for task-based fMRI, especially at the single subject level. This work, which focuses exclusively on individual subject results, compares ME-ICA to single-echo fMRI and a voxel-wise T2(⁎) weighted combination of multi-echo data for task-based fMRI under the following scenarios: cardiac-gated block designs, constant repetition time (TR) block designs, and constant TR rapid event-related designs. Performance is evaluated primarily in terms of sensitivity (i.e., activation extent, activation magnitude, percent detected trials and effect size estimates) using five different tasks expected to evoke neuronal activity in a distributed set of regions. The ME-ICA algorithm significantly outperformed all other evaluated processing alternatives in all scenarios. Largest improvements were observed for the cardiac-gated dataset, where ME-ICA was able to reliably detect and remove non-neural T1 signal fluctuations caused by non-constant repetition times. Although ME-ICA also outperformed the other options in terms of percent detection of individual trials for rapid event-related experiments, only 46% of all events were detected after ME-ICA; suggesting additional improvements in sensitivity are required to reliably detect individual short event occurrences. We conclude the manuscript with a detailed evaluation of ME-ICA outcomes and a discussion of how the ME-ICA algorithm could be further improved. Overall, our results suggest that ME-ICA constitutes a versatile, powerful approach for advanced denoising of task-based fMRI, not just resting-state data.

Task Dependence, Tissue Specificity, and Spatial Distribution of Widespread Activations in Large Single-Subject Functional MRI Datasets at 7T.

It was recently shown that when large amounts of task-based blood oxygen level-dependent (BOLD) data are combined to increase contrast- and temporal signal-to-noise ratios, the majority of the brain shows significant hemodynamic responses time-locked with the experimental paradigm. Here, we investigate the biological significance of such widespread activations. First, the relationship between activation extent and task demands was investigated by varying cognitive load across participants. Second, the tissue specificity of responses was probed using the better BOLD signal localization capabilities of a 7T scanner. Finally, the spatial distribution of 3 primary response types--namely positively sustained (pSUS), negatively sustained (nSUS), and transient--was evaluated using a newly defined voxel-wise waveshape index that permits separation of responses based on their temporal signature. About 86% of gray matter (GM) became significantly active when all data entered the analysis for the most complex task. Activation extent scaled with task load and largely followed the GM contour. The most common response type was nSUS BOLD, irrespective of the task. Our results suggest that widespread activations associated with extremely large single-subject functional magnetic resonance imaging datasets can provide valuable information about the functional organization of the brain that goes undetected in smaller sample sizes.

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

The brain is the body's largest energy consumer, even in the absence of demanding tasks. Electrophysiologists report on-going neuronal firing during stimulation or task in regions beyond those of primary relationship to the perturbation. Although the biological origin of consciousness remains elusive, it is argued that it emerges from complex, continuous whole-brain neuronal collaboration. Despite converging evidence suggesting the whole brain is continuously working and adapting to anticipate and actuate in response to the environment, over the last 20 y, task-based functional MRI (fMRI) have emphasized a localizationist view of brain function, with fMRI showing only a handful of activated regions in response to task/stimulation. Here, we challenge that view with evidence that under optimal noise conditions, fMRI activations extend well beyond areas of primary relationship to the task; and blood-oxygen level-dependent signal changes correlated with task-timing appear in over 95% of the brain for a simple visual stimulation plus attention control task. Moreover, we show that response shape varies substantially across regions, and that whole-brain parcellations based on those differences produce distributed clusters that are anatomically and functionally meaningful, symmetrical across hemispheres, and reproducible across subjects. These findings highlight the exquisite detail lying in fMRI signals beyond what is normally examined, and emphasize both the pervasiveness of false negatives, and how the sparseness of fMRI maps is not a result of localized brain function, but a consequence of high noise and overly strict predictive response models.

Connectivity trajectory across lifespan differentiates the precuneus from the default network.

The default network of the human brain has drawn much attention due to its relevance to various brain disorders, cognition, and behavior. However, its functional components and boundaries have not been precisely defined. There is no consensus as to whether the precuneus, a hub in the functional connectome, acts as part of the default network. This discrepancy is more critical for brain development and aging studies: it is not clear whether age has a stronger impact on the default network or precuneus, or both. We used Generalized Ranking and Averaging Independent Component Analysis by Reproducibility (gRAICAR) to investigate the lifespan trajectories of intrinsic functional networks. By estimating individual-specific spatial components and aligning them across subjects, gRAICAR measures the spatial variation of component maps across a population without constraining the same components to appear in every subject. In a cross-lifespan fMRI dataset (N=126, 7-85years old), we observed stronger age dependence in the spatial pattern of a precuneus-dorsal posterior cingulate cortex network compared to the default network, despite the fact that the two networks exhibit considerable spatial overlap and temporal correlation. These results remained even when analyses were restricted to a subpopulation with very similar head motion across age. Our analyses further showed that the two networks tend to merge with increasing age. Post-hoc analyses of functional connectivity confirmed the distinguishable cross-lifespan trajectories between the two networks. Based on these observations, we proposed a dynamic model of cross-lifespan functional segregation and integration between the two networks, suggesting that the precuneus network may have a different functional role than the default network, which declines with age. These findings have implications for understanding the functional roles of the default network, gaining insight into its dynamics throughout life, and guiding interpretation of alterations in brain disorders.

Neuroanatomical changes observed over the course of a human pregnancy.

Pregnancy is a period of profound hormonal and physiological change experienced by millions of women annually, yet the neural changes unfolding in the maternal brain throughout gestation have not been studied in humans. Leveraging precision imaging, we mapped neuroanatomical changes in an individual from preconception through two years postpartum. Pronounced decreases in gray matter volume and cortical thickness were evident across the brain, which stand in contrast to increases in white matter microstructural integrity, ventricle volume, and cerebrospinal fluid, with few regions untouched by the transition to motherhood. This dataset serves as the first comprehensive map of the human brain across gestation, providing an open-access resource for the brain imaging community to stimulate further exploration and discovery.

Visual imagery vividness correlates with afterimage conscious perception.

Afterimages are illusory, visual conscious perceptions. A widely accepted theory is that afterimages are caused by retinal signaling that continues after the physical disappearance of a light stimulus. However, afterimages have been reported without preceding visual, sensory stimulation (e.g., conditioned afterimages and afterimages induced by illusory vision). These observations suggest the role of top-down, brain mechanisms in afterimage conscious perception. Therefore, some afterimages may share perceptual features with sensory-independent conscious perceptions (e.g., imagery, hallucinations, and dreams) that occur without bottom-up, sensory input. In the current investigation, we tested for a link between the vividness of visual imagery and afterimage conscious perception. Participants reported their vividness of visual imagery and perceived sharpness, contrast, and duration of negative afterimages. The afterimage perceptual features were acquired using perception matching paradigms that were validated on image stimuli. Relating these perceptual reports revealed that the vividness of visual imagery positively correlated with afterimage contrast and sharpness. These behavioral results support shared neural mechanisms between visual imagery and afterimages. This study encourages future research combining neurophysiology recording methods and afterimage paradigms to directly examine the neural mechanisms of afterimage conscious perception.

The Past, Present, and Future of the Brain Imaging Data Structure (BIDS).

The Brain Imaging Data Structure (BIDS) is a community-driven standard for the organization of data and metadata from a growing range of neuroscience modalities. This paper is meant as a history of how the standard has developed and grown over time. We outline the principles behind the project, the mechanisms by which it has been extended, and some of the challenges being addressed as it evolves. We also discuss the lessons learned through the project, with the aim of enabling researchers in other domains to learn from the success of BIDS.

Dynamic functional connectivity: promise, issues, and interpretations.

The brain must dynamically integrate, coordinate, and respond to internal and external stimuli across multiple time scales. Non-invasive measurements of brain activity with fMRI have greatly advanced our understanding of the large-scale functional organization supporting these fundamental features of brain function. Conclusions from previous resting-state fMRI investigations were based upon static descriptions of functional connectivity (FC), and only recently studies have begun to capitalize on the wealth of information contained within the temporal features of spontaneous BOLD FC. Emerging evidence suggests that dynamic FC metrics may index changes in macroscopic neural activity patterns underlying critical aspects of cognition and behavior, though limitations with regard to analysis and interpretation remain. Here, we review recent findings, methodological considerations, neural and behavioral correlates, and future directions in the emerging field of dynamic FC investigations.