Sympathetic activity contributes to the fMRI signal

Arousal as a universal embedding for spatiotemporal brain dynamics

Neural activity in awake organisms shows widespread and spatiotemporally diverse correlations with behavioral and physiological measurements. We propose that this covariation reflects in part the dynamics of a unified, arousal-related process that regulates brain-wide physiology on the timescale of seconds. Taken together with theoretical foundations in dynamical systems, this interpretation leads us to a surprising prediction: that a single, scalar measurement of arousal (e.g., pupil diameter) should suffice to reconstruct the continuous evolution of multimodal, spatiotemporal measurements of large-scale brain physiology. To test this hypothesis, we perform multimodal, cortex-wide optical imaging and behavioral monitoring in awake mice. We demonstrate that spatiotemporal measurements of neuronal calcium, metabolism, and blood-oxygen can be accurately and parsimoniously modeled from a low-dimensional state-space reconstructed from the time history of pupil diameter. Extending this framework to behavioral and electrophysiological measurements from the Allen Brain Observatory, we demonstrate the ability to integrate diverse experimental data into a unified generative model via mappings from an intrinsic arousal manifold. Our results support the hypothesis that spontaneous, spatially structured fluctuations in brain-wide physiology-widely interpreted to reflect regionally-specific neural communication-are in large part reflections of an arousal-related process. This enriched view of arousal dynamics has broad implications for interpreting observations of brain, body, and behavior as measured across modalities, contexts, and scales.

Effect of sleep deprivation and NREM sleep stage on physiological brain pulsations

Introduction

Sleep increases brain fluid transport and the power of pulsations driving the fluids. We investigated how sleep deprivation or electrophysiologically different stages of non-rapid-eye-movement (NREM) sleep affect the human brain pulsations.

Methods

Fast functional magnetic resonance imaging (fMRI) was performed in healthy subjects (n = 23) with synchronous electroencephalography (EEG), that was used to verify arousal states (awake, N1 and N2 sleep). Cardiorespiratory rates were verified with physiological monitoring. Spectral power analysis assessed the strength, and spectral entropy assessed the stability of the pulsations.

Results

In N1 sleep, the power of vasomotor (VLF < 0.1 Hz), but not cardiorespiratory pulsations, intensified after sleep deprived vs. non-sleep deprived subjects. The power of all three pulsations increased as a function of arousal state (N2 > N1 > awake) encompassing brain tissue in both sleep stages, but extra-axial CSF spaces only in N2 sleep. Spectral entropy of full band and respiratory pulsations decreased most in N2 sleep stage, while cardiac spectral entropy increased in ventricles.

Discussion

In summary, the sleep deprivation and sleep depth, both increase the power and harmonize the spectral content of human brain pulsations.

Resting-state global brain activity affects early β-amyloid accumulation in default mode network

It remains unclear why β-amyloid (Aβ) plaque, a hallmark pathology of Alzheimer's disease (AD), first accumulates cortically in the default mode network (DMN), years before AD diagnosis. Resting-state low-frequency ( < 0.1 Hz) global brain activity recently was linked to AD, presumably due to its role in glymphatic clearance. Here we show that the preferential Aβ accumulation in the DMN at the early stage of Aβ pathology was associated with the preferential reduction of global brain activity in the same regions. This can be partly explained by its failure to reach these regions as propagating waves. Together, these findings highlight the important role of resting-state global brain activity in early preferential Aβ deposition in the DMN.

Decoupling Between Brain Activity and Cerebrospinal Fluid Movement in Neurological Disorders

Recent research has identified a link between the global mean signal of resting-state functional MRI (fMRI) and macro-scale cerebrospinal fluid movement, indicating the potential link between this resting-state dynamic and brain waste clearance. Consistent with this notion, the strength of this coupling has been associated with multiple neurodegenerative disease pathologies, especially the build-up of toxic proteins. This article aimed to review the latest advancements in this research area, emphasizing studies on spontaneous global brain activity that is tightly linked to the global mean resting-state fMRI signal, and aimed to discuss potential mechanisms through which this activity and associated physiological modulations might affect brain waste clearance. The available evidence supports the presence of a highly organized global brain activity that is linked to arousal and memory systems. This global brain dynamic, along with its associated physiological modulations, has the potential to influence brain waste clearance through multiple pathways through multiple pathways. LEVEL OF EVIDENCE: 2 TECHNICAL EFFICACY: Stage 3.

An autonomic mode of brain activity

The relevance of interactions between autonomic and central nervous systems remains unclear for human brain function and health, particularly when both systems are challenged under sleep deprivation (SD). We measured brain activity (with fMRI), pulse and respiratory signals, and baseline brain amyloid beta burden (with PET) in healthy participants. We found that SD relative to rested wakefulness (RW) resulted in a significant increase in synchronized low frequency (LF, < 0.1 Hz) activity in an autonomically-related network (AN), including dorsal attention, visual, and sensorimotor regions, which we previously found to have consistent temporal coupling with LF pulse signal changes (regulated by sympathetic tone). SD resulted in a significant phase coherence between the LF component of the pulse signal and a medial network with peak effects in the midbrain reticular formation, and between LF component of the respiratory variations (regulated by respiratory motor output) and a cerebellar network. The LF power of AN during SD was significantly and independently correlated with pulse-medial network and respiratory-cerebellar network phase coherences (total adjusted R2 = 0.78). Higher LF power of AN during SD (but not RW) was associated with lower amyloid beta burden (Cohen's d = 0.8). In sum, SD triggered an autonomic mode of synchronized brain activity that was associated with distinct autonomic-central interactions. Findings highlight the direct relevance of global cortical synchronization to brain clearance mechanisms.

Global brain activity and its coupling with cerebrospinal fluid flow is related to tau pathology

Amyloid-β (Aβ) and tau deposition constitute Alzheimer's disease (AD) neuropathology. Cortical tau deposits first in the entorhinal cortex and hippocampus and then propagates to neocortex in an Aβ-dependent manner. Tau also tends to accumulate earlier in higher-order association cortex than in lower-order primary sensory-motor cortex. While previous research has examined the production and spread of tau, little attention has been paid to its clearance. Low-frequency (<0.1 Hz) global brain activity during the resting state is coupled with cerebrospinal fluid (CSF) flow and potentially reflects glymphatic clearance. Here we report that tau deposition in subjects with evaluated Aβ, accompanied by cortical thinning and cognitive decline, is strongly associated with decreased coupling between CSF flow and global brain activity. Substantial modulation of global brain activity is also manifested as propagating waves of brain activation between higher- and lower-order regions, resembling tau spreading. Together, the findings suggest an important role of resting-state global brain activity in AD tau pathology.

Altered brain dynamic in major depressive disorder: state and trait features

Temporal neural synchrony disruption can be linked to a variety of symptoms of major depressive disorder (MDD), including mood rigidity and the inability to break the cycle of negative emotion or attention biases. This might imply that altered dynamic neural synchrony may play a role in the persistence and exacerbation of MDD symptoms. Our study aimed to investigate the changes in whole-brain dynamic patterns of the brain functional connectivity and activity related to depression using the hidden Markov model (HMM) on resting-state functional magnetic resonance imaging (rs-fMRI) data. We compared the patterns of brain functional dynamics in a large sample of 314 patients with MDD (65.9% female; age (mean ± standard deviation): 35.9 ± 13.4) and 498 healthy controls (59.4% female; age: 34.0 ± 12.8). The HMM model was used to explain variations in rs-fMRI functional connectivity and averaged functional activity across the whole-brain by using a set of six unique recurring states. This study compared the proportion of time spent in each state and the average duration of visits to each state to assess stability between different groups. Compared to healthy controls, patients with MDD showed significantly higher proportional time spent and temporal stability in a state characterized by weak functional connectivity within and between all brain networks and relatively strong averaged functional activity of regions located in the somatosensory motor (SMN), salience (SN), and dorsal attention (DAN) networks. Both proportional time spent and temporal stability of this brain state was significantly associated with depression severity. Healthy controls, in contrast to the MDD group, showed proportional time spent and temporal stability in a state with relatively strong functional connectivity within and between all brain networks but weak averaged functional activity across the whole brain. These findings suggest that disrupted brain functional synchrony across time is present in MDD and associated with current depression severity.

Identifying the distinct spectral dynamics of laminar-specific interhemispheric connectivity with bilateral line-scanning fMRI

Despite extensive efforts to identify interhemispheric functional connectivity (FC) with resting-state (rs-) fMRI, correlated low-frequency rs-fMRI signal fluctuation across homotopic cortices originates from multiple sources. It remains challenging to differentiate circuit-specific FC from global regulation. Here, we developed a bilateral line-scanning fMRI method to detect laminar-specific rs-fMRI signals from homologous forepaw somatosensory cortices with high spatial and temporal resolution in rat brains. Based on spectral coherence analysis, two distinct bilateral fluctuation spectral features were identified: ultra-slow fluctuation (<0.04 Hz) across all cortical laminae versus Layer (L) 2/3-specific evoked BOLD at 0.05 Hz based on 4 s on/16 s off block design and resting-state fluctuations at 0.08-0.1 Hz. Based on the measurements of evoked BOLD signal at corpus callosum (CC), this L2/3-specific 0.05 Hz signal is likely associated with neuronal circuit-specific activity driven by the callosal projection, which dampened ultra-slow oscillation less than 0.04 Hz. Also, the rs-fMRI power variability clustering analysis showed that the appearance of L2/3-specific 0.08-0.1 Hz signal fluctuation is independent of the ultra-slow oscillation across different trials. Thus, distinct laminar-specific bilateral FC patterns at different frequency ranges can be identified by the bilateral line-scanning fMRI method.

Cognition is entangled with metabolism: relevance for resting-state EEG-fMRI

The brain is a living organ with distinct metabolic constraints. However, these constraints are typically considered as secondary or supportive of information processing which is primarily performed by neurons. The default operational definition of neural information processing is that (1) it is ultimately encoded as a change in individual neuronal firing rate as this correlates with the presentation of a peripheral stimulus, motor action or cognitive task. Two additional assumptions are associated with this default interpretation: (2) that the incessant background firing activity against which changes in activity are measured plays no role in assigning significance to the extrinsically evoked change in neural firing, and (3) that the metabolic energy that sustains this background activity and which correlates with differences in neuronal firing rate is merely a response to an evoked change in neuronal activity. These assumptions underlie the design, implementation, and interpretation of neuroimaging studies, particularly fMRI, which relies on changes in blood oxygen as an indirect measure of neural activity. In this article we reconsider all three of these assumptions in light of recent evidence. We suggest that by combining EEG with fMRI, new experimental work can reconcile emerging controversies in neurovascular coupling and the significance of ongoing, background activity during resting-state paradigms. A new conceptual framework for neuroimaging paradigms is developed to investigate how ongoing neural activity is "entangled" with metabolism. That is, in addition to being recruited to support locally evoked neuronal activity (the traditional hemodynamic response), changes in metabolic support may be independently "invoked" by non-local brain regions, yielding flexible neurovascular coupling dynamics that inform the cognitive context. This framework demonstrates how multimodal neuroimaging is necessary to probe the neurometabolic foundations of cognition, with implications for the study of neuropsychiatric disorders.

Sex-specific age-related changes in glymphatic function assessed by resting-state functional magnetic resonance imaging

The glymphatic system that clears out brain wastes, such as amyloid-β (Aβ) and tau, through cerebrospinal fluid (CSF) flow may play an important role in aging and dementias. However, a lack of non-invasive tools to assess the glymphatic function in humans hindered the understanding of the glymphatic changes in healthy aging. The global infra-slow (<0.1 Hz) brain activity measured by the global mean resting-state fMRI signal (gBOLD) was recently found to be coupled by large CSF movements. This coupling has been used to measure the glymphatic process and found to correlate with various pathologies of Alzheimer's disease (AD), including Aβ pathology. Using resting-state fMRI data from a large group of 719 healthy aging participants, we examined the sex-specific changes of the gBOLD-CSF coupling, as a measure of glymphatic function, over a wide age range between 36-100 years old. We found that this coupling index remains stable before around age 55 and then starts to decline afterward, particularly in females. Menopause may contribute to the accelerated decline in females.

Neural activity induced by sensory stimulation can drive large-scale cerebrospinal fluid flow during wakefulness in humans

Cerebrospinal fluid (CSF) flow maintains healthy brain homeostasis, facilitating solute transport and the exchange of brain waste products. CSF flow is thus important for brain health, but the mechanisms that control its large-scale movement through the ventricles are not well understood. While it is well established that CSF flow is modulated by respiratory and cardiovascular dynamics, recent work has also demonstrated that neural activity is coupled to large waves of CSF flow in the ventricles during sleep. To test whether the temporal coupling between neural activity and CSF flow is in part due to a causal relationship, we investigated whether CSF flow could be induced by driving neural activity with intense visual stimulation. We manipulated neural activity with a flickering checkerboard visual stimulus and found that we could drive macroscopic CSF flow in the human brain. The timing and amplitude of CSF flow was matched to the visually evoked hemodynamic responses, suggesting neural activity can modulate CSF flow via neurovascular coupling. These results demonstrate that neural activity can contribute to driving CSF flow in the human brain and that the temporal dynamics of neurovascular coupling can explain this effect.

Psychedelics and fNIRS neuroimaging: exploring new opportunities

In this Outlook paper, we explain to the optical neuroimaging community as well as the psychedelic research community the great potential of using optical neuroimaging with functional near-infrared spectroscopy (fNIRS) to further explore the changes in brain activity induced by psychedelics. We explain why we believe now is the time to exploit the momentum of the current resurgence of research on the effects of psychedelics and the momentum of the increasing progress and popularity of the fNIRS technique to establish fNIRS in psychedelic research. With this article, we hope to contribute to this development.

Brain and cardiovascular-related changes are associated with aging, hypertension, and atrial fibrillation

PURPOSE

The neural pathways in which the brain regulates the cardiovascular system is via sympathetic and parasympathetic control of the heart and sympathetic control of the systemic vasculature. Various cortical and sub-cortical sites are involved, but how these critical brain regions for cardiovascular control are altered in healthy aging and other risk conditions that may contribute to cardiovascular disease is uncertain.

METHODS

Here we review the functional and structural brain changes in healthy aging, hypertension, and atrial fibrillation - noting their potential influence on the autonomic nervous system and hence on cardiovascular control.

RESULTS

Evidence suggests that aging, hypertension, and atrial fibrillation are each associated with functional and structural changes in specific areas of the central nervous system involved in autonomic control. Increased muscle sympathetic nerve activity (MSNA) and significant alterations in the brain regions involved in the default mode network are commonly reported in aging, hypertension, and atrial fibrillation.

CONCLUSIONS

Further studies using functional and structural magnetic resonance imaging (MRI) coupled with autonomic nerve activity in healthy aging, hypertension, and atrial fibrillation promise to reveal the underlying brain circuitry modulating the abnormal sympathetic nerve activity in these conditions. This understanding will guide future therapies to rectify dysregulation of autonomic and cardiovascular control by the brain.

An orderly sequence of autonomic and neural events at transient arousal changes

Resting-state functional magnetic resonance imaging (rsfMRI) allows the study of functional brain connectivity based on spatially structured variations in neuronal activity. Proper evaluation of connectivity requires removal of non-neural contributions to the fMRI signal, in particular hemodynamic changes associated with autonomic variability. Regression analysis based on autonomic indicator signals has been used for this purpose, but may be inadequate if neuronal and autonomic activities covary. To investigate this potential co-variation, we performed rsfMRI experiments while concurrently acquiring electroencephalography (EEG) and autonomic indicator signals, including heart rate, respiratory depth, and peripheral vascular tone. We identified a recurrent and systematic spatiotemporal pattern of fMRI (named as fMRI cascade), which features brief signal reductions in salience and default-mode networks and the thalamus, followed by a biphasic global change with a sensory-motor dominance. This fMRI cascade, which was mostly observed during eyes-closed condition, was accompanied by large EEG and autonomic changes indicative of arousal modulations. Importantly, the removal of the fMRI cascade dynamics from rsfMRI diminished its correlations with various signals. These results suggest that the rsfMRI correlations with various physiological and neural signals are not independent but arise, at least partly, from the fMRI cascades and associated neural and physiological changes at arousal modulations.

A temporal sequence of thalamic activity unfolds at transitions in behavioral arousal state

Awakening from sleep reflects a profound transformation in neural activity and behavior. The thalamus is a key controller of arousal state, but whether its diverse nuclei exhibit coordinated or distinct activity at transitions in behavioral arousal state is unknown. Using fast fMRI at ultra-high field (7 Tesla), we measured sub-second activity across thalamocortical networks and within nine thalamic nuclei to delineate these dynamics during spontaneous transitions in behavioral arousal state. We discovered a stereotyped sequence of activity across thalamic nuclei and cingulate cortex that preceded behavioral arousal after a period of inactivity, followed by widespread deactivation. These thalamic dynamics were linked to whether participants subsequently fell back into unresponsiveness, with unified thalamic activation reflecting maintenance of behavior. These results provide an outline of the complex interactions across thalamocortical circuits that orchestrate behavioral arousal state transitions, and additionally, demonstrate that fast fMRI can resolve sub-second subcortical dynamics in the human brain.

Cerebrovascular activity is a major factor in the cerebrospinal fluid flow dynamics

Cerebrospinal fluid (CSF) provides physical protection to the central nervous system as well as an essential homeostatic environment for the normal functioning of neurons. Additionally, it has been proposed that the pulsatile movement of CSF may assist in glymphatic clearance of brain metabolic waste products implicated in neurodegeneration. In awake humans, CSF flow dynamics are thought to be driven primarily by cerebral blood volume fluctuations resulting from a number of mechanisms, including a passive vascular response to blood pressure variations associated with cardiac and respiratory cycles. Recent research has shown that mechanisms that rely on the action of vascular smooth muscle cells ("cerebrovascular activity") such as neuronal activity, changes in intravascular CO2, and autonomic activation from the brainstem, may lead to CSF pulsations as well. Nevertheless, the relative contribution of these mechanisms to CSF flow remains unclear. To investigate this further, we developed an MRI approach capable of disentangling and quantifying CSF flow components of different time scales associated with these mechanisms. This approach was evaluated on human control subjects (n = 12) performing intermittent voluntary deep inspirations, by determining peak flow velocities and displaced volumes between these mechanisms in the fourth ventricle. We found that peak flow velocities were similar between the different mechanisms, while displaced volumes per cycle were about a magnitude larger for deep inspirations. CSF flow velocity peaked at around 10.4 s (range 7.1-14.8 s, n = 12) following deep inspiration, consistent with known cerebrovascular activation delays for this autonomic challenge. These findings point to an important role of cerebrovascular activity in the genesis of CSF pulsations. Other regulatory triggers for cerebral blood flow such as autonomic arousal and orthostatic challenges may create major CSF pulsatile movement as well. Future quantitative comparison of these and possibly additional types of CSF pulsations with the proposed approach may help clarify the conditions that affect CSF flow dynamics.

A parsimonious description of global functional brain organization in three spatiotemporal patterns

Resting-state functional magnetic resonance imaging (MRI) has yielded seemingly disparate insights into large-scale organization of the human brain. The brain's large-scale organization can be divided into two broad categories: zero-lag representations of functional connectivity structure and time-lag representations of traveling wave or propagation structure. In this study, we sought to unify observed phenomena across these two categories in the form of three low-frequency spatiotemporal patterns composed of a mixture of standing and traveling wave dynamics. We showed that a range of empirical phenomena, including functional connectivity gradients, the task-positive/task-negative anti-correlation pattern, the global signal, time-lag propagation patterns, the quasiperiodic pattern and the functional connectome network structure, are manifestations of these three spatiotemporal patterns. These patterns account for much of the global spatial structure that underlies functional connectivity analyses and unifies phenomena in resting-state functional MRI previously thought distinct.

The Spatiotemporal Dynamics of Cerebral Autoregulation in Functional Magnetic Resonance Imaging

The thigh-cuff release (TCR) maneuver is a physiological challenge that is widely used to assess dynamic cerebral autoregulation (dCA). It is often applied in conjunction with Transcranial Doppler ultrasound (TCD), which provides temporal information of the global flow response in the brain. This established method can only yield very limited insights into the regional variability of dCA, whereas functional MRI (fMRI) has the ability to reveal the spatial distribution of flow responses in the brain with high spatial resolution. The aim of this study was to use whole-brain blood-oxygenation-level-dependent (BOLD) fMRI to characterize the spatiotemporal dynamics of the flow response to the TCR challenge, and thus pave the way toward mapping dCA in the brain. We used a data driven approach to derive a novel basis set that was then used to provide a voxel-wise estimate of the TCR associated haemodynamic response function (HRF TCR ). We found that the HRF TCR evolves with a specific spatiotemporal pattern, with gray and white matter showing an asynchronous response, which likely reflects the anatomical structure of cerebral blood supply. Thus, we propose that TCR challenge fMRI is a promising method for mapping spatial variability in dCA, which will likely prove to be clinically advantageous.

Systemic physiology augmented functional near-infrared spectroscopy: a powerful approach to study the embodied human brain

In this Outlook paper, we explain why an accurate physiological interpretation of functional near-infrared spectroscopy (fNIRS) neuroimaging signals is facilitated when systemic physiological activity (e.g., cardiorespiratory and autonomic activity) is measured simultaneously by employing systemic physiology augmented functional near-infrared spectroscopy (SPA-fNIRS). The rationale for SPA-fNIRS is twofold: (i) SPA-fNIRS enables a more complete interpretation and understanding of the fNIRS signals measured at the head since they contain components originating from neurovascular coupling and from systemic physiological sources. The systemic physiology signals measured with SPA-fNIRS can be used for regressing out physiological confounding components in fNIRS signals. Misinterpretations can thus be minimized. (ii) SPA-fNIRS enables to study the embodied brain by linking the brain with the physiological state of the entire body, allowing novel insights into their complex interplay. We envisage the SPA-fNIRS approach will become increasingly important in the future.

Task-related hemodynamic responses in human early visual cortex are modulated by task difficulty and behavioral performance

Early visual cortex exhibits widespread hemodynamic responses in the absence of visual stimulation, which are entrained to the timing of a task and not predicted by local spiking or local field potential (LFP). Such task-related responses ('TRRs') covary with reward magnitude and physiological signatures of arousal. It is unknown, however, if TRRs change on a trial-to-trial basis according to behavioral performance and task difficulty. If so, this would suggest that TRRs reflect arousal on a trial-to-trial timescale and covary with critical task and behavioral variables. We measured fMRI-BOLD responses in the early visual cortex of human observers performing an orientation discrimination task consisting of separate easy and hard runs of trials. Stimuli were presented in a small portion of one hemifield, but the fMRI response was measured in the ipsilateral hemisphere, far from the stimulus representation and focus of spatial attention. TRRs scaled in amplitude with task difficulty, behavioral accuracy, reaction time, and lapses across trials. These modulations were not explained by the influence of respiration, cardiac activity, or head movement on the fMRI signal. Similar modulations with task difficulty and behavior were observed in pupil size. These results suggest that TRRs reflect arousal and behavior on the timescale of individual trials.

Laminar-specific functional connectivity mapping with multi-slice line-scanning fMRI

Despite extensive studies detecting laminar functional magnetic resonance imaging (fMRI) signals to illustrate the canonical microcircuit, the spatiotemporal characteristics of laminar-specific information flow across cortical regions remain to be fully investigated in both evoked and resting conditions at different brain states. Here, we developed a multislice line-scanning fMRI (MS-LS) method to detect laminar fMRI signals in adjacent cortical regions with high spatial (50 μm) and temporal resolution (100 ms) in anesthetized rats. Across different trials, we detected either laminar-specific positive or negative blood-oxygen-level-dependent (BOLD) responses in the surrounding cortical region adjacent to the most activated cortex under the evoked condition. Specifically, in contrast to typical Layer (L) 4 correlation across different regions due to the thalamocortical projections for trials with positive BOLD, a strong correlation pattern specific in L2/3 was detected for trials with negative BOLD in adjacent regions, which indicated brain state-dependent laminar-fMRI responses based on corticocortical interaction. Also, in resting-state (rs-) fMRI study, robust lag time differences in L2/3, 4, and 5 across multiple cortices represented the low-frequency rs-fMRI signal propagation from caudal to rostral slices. In summary, our study provided a unique laminar fMRI mapping scheme to better characterize trial-specific intra- and inter-laminar functional connectivity in evoked and resting-state MS-LS.

Rapid fluctuations in brain oxygenation during glucose-drinking behavior in trained rats

Proper inflow of oxygen into brain tissue is essential for maintaining normal neural functions. Although oxygen levels in the brain's extracellular space depend upon a balance between its delivery from arterial blood and its metabolic consumption, the use of high-speed electrochemical detection revealed rapid increases in brain oxygen levels elicited by various salient sensory stimuli. These stimuli also increase intrabrain heat production, an index of metabolic neural activation, but these changes are slower and more prolonged than changes in oxygen levels. Therefore, under physiological conditions, the oxygen inflow into brain tissue exceeds its loss due to consumption, thus preventing any metabolic deficit. Here, we used oxygen sensors coupled with amperometry to examine the pattern of real-time oxygen fluctuations in the nucleus accumbens during glucose-drinking behavior in trained rats. Following the exposure to a glucose-containing cup, oxygen levels rapidly increased, peaked when the rat initiated drinking, and relatively decreased during consumption. Similar oxygen changes but more episodic drinking occurred when Stevia, a calorie-free sweet substance, was substituted for glucose. When water was substituted for glucose, rats tested the water but refused to consume all of it. Although the basic pattern of oxygen changes during this water test was similar to that with glucose drinking, the increases were larger. Finally, oxygen increases were significantly larger when rats were exposed to concealed glucose and made multiple unsuccessful attempts to obtain and consume it. Based on these data, we discuss the mechanisms underlying behavior-related brain oxygen fluctuations and their functional significance.NEW & NOTEWORTHY Oxygen sensors coupled with high-speed amperometry were used to examine brain oxygen fluctuations during glucose-drinking behavior in trained rats. Oxygen levels rapidly increased following presentation of a glucose-contained cup, peaking at the initiation of glucose drinking, and relatively decreasing during drinking. Oxygen increases were larger when rats were exposed to concealed glucose and made multiple attempts to obtain it. We discuss the mechanisms underlying behavior-related brain oxygen fluctuations and their functional significance.

Autonomic arousals contribute to brain fluid pulsations during sleep

During sleep, slow waves of neuro-electrical activity engulf the human brain and aid in the consolidation of memories. Recent research suggests that these slow waves may also promote brain health by facilitating the removal of metabolic waste, possibly by orchestrating the pulsatile flow of cerebro-spinal fluid (CSF) through local neural control over vascular tone. To investigate the role of slow waves in the generation of CSF pulsations, we analyzed functional MRI data obtained across the full sleep-wake cycle and during a respiratory task during wakefulness. This revealed a novel generating mechanism that relies on the autonomic regulation of cerebral vascular tone without requiring slow electrocortical activity or even sleep. Therefore, the role of CSF pulsations in brain waste clearance may, in part, depend on proper autoregulatory control of cerebral blood flow.

State-related neural influences on fMRI connectivity estimation

The spatiotemporal structure of functional magnetic resonance imaging (fMRI) signals has provided a valuable window into the network underpinnings of human brain function and dysfunction. Although some cross-regional temporal correlation patterns (functional connectivity; FC) exhibit a high degree of stability across individuals and species, there is growing acknowledgment that measures of FC can exhibit marked changes over a range of temporal scales. Further, FC can covary with experimental task demands and ongoing neural processes linked to arousal, consciousness and perception, cognitive and affective state, and brain-body interactions. The increased recognition that such interrelated neural processes modulate FC measurements has raised both challenges and new opportunities in using FC to investigate brain function. Here, we review recent advances in the quantification of neural effects that shape fMRI FC and discuss the broad implications of these findings in the design and analysis of fMRI studies. We also discuss how a more complete understanding of the neural factors that shape FC measurements can resolve apparent inconsistencies in the literature and lead to more interpretable conclusions from fMRI studies.

Single-neuron firing cascades underlie global spontaneous brain events

The resting brain consumes enormous energy and shows highly organized spontaneous activity. To investigate how this activity is manifest among single neurons, we analyzed spiking discharges of ∼10,000 isolated cells recorded from multiple cortical and subcortical regions of the mouse brain during immobile rest. We found that firing of a significant proportion (∼70%) of neurons conformed to a ubiquitous, temporally sequenced cascade of spiking that was synchronized with global events and elapsed over timescales of 5 to 10 s. Across the brain, two intermixed populations of neurons supported orthogonal cascades. The relative phases of these cascades determined, at each moment, the response magnitude evoked by an external visual stimulus. Furthermore, the spiking of individual neurons embedded in these cascades was time locked to physiological indicators of arousal, including local field potential power, pupil diameter, and hippocampal ripples. These findings demonstrate that the large-scale coordination of low-frequency spontaneous activity, which is commonly observed in brain imaging and linked to arousal, sensory processing, and memory, is underpinned by sequential, large-scale temporal cascades of neuronal spiking across the brain.

Imaging the temporal dynamics of brain states with highly sampled fMRI

The spontaneous dynamics of the brain modulate its function from moment to moment, shaping neural computation and cognition. Functional MRI (fMRI), while classically used as a tool for spatial localization, is increasingly being used to identify the temporal dynamics of brain activity. fMRI analyses focused on the temporal domain have revealed important new information about the dynamics underlying states such as arousal, attention, and sleep. Dense temporal sampling – either by using fast fMRI acquisition, or multiple repeated scan sessions within individuals – can further enrich the information present in these studies. This review focuses on recent developments in using fMRI to identify dynamics across brain states, particularly vigilance and sleep states, and the potential for highly temporally sampled fMRI to answer these questions.

Representation learning of resting state fMRI with variational autoencoder

Resting state functional magnetic resonance imaging (rsfMRI) data exhibits complex but structured patterns. However, the underlying origins are unclear and entangled in rsfMRI data. Here we establish a variational auto-encoder, as a generative model trainable with unsupervised learning, to disentangle the unknown sources of rsfMRI activity. After being trained with large data from the Human Connectome Project, the model has learned to represent and generate patterns of cortical activity and connectivity using latent variables. The latent representation and its trajectory represent the spatiotemporal characteristics of rsfMRI activity. The latent variables reflect the principal gradients of the latent trajectory and drive activity changes in cortical networks. Representational geometry captured as covariance or correlation between latent variables, rather than cortical connectivity, can be used as a more reliable feature to accurately identify subjects from a large group, even if only a short period of data is available in each subject. Our results demonstrate that VAE is a valuable addition to existing tools, particularly suited for unsupervised representation learning of resting state fMRI activity.

High γ Activity in Cortex and Hippocampus Is Correlated with Autonomic Tone during Sleep

Studies in animals have demonstrated a strong relationship between cortical and hippocampal activity, and autonomic tone. However, the extent, distribution, and nature of this relationship have not been investigated with intracranial recordings in humans during sleep. Cortical and hippocampal population neuronal firing was estimated from high γ band activity (HG) from 70 to 110 Hz in local field potentials (LFPs) recorded from 15 subjects (nine females) during nonrapid eye movement (NREM) sleep. Autonomic tone was estimated from heart rate variability (HRV). HG and HRV were significantly correlated in the hippocampus and multiple cortical sites in NREM stages N1-N3. The average correlation between HG and HRV could be positive or negative across patients given anatomic location and sleep stage and was most profound in lateral temporal lobe in N3, suggestive of greater cortical activity associated with sympathetic tone. Patient-wide correlation was related to δ band activity (1-4 Hz), which is known to be correlated with high γ activity during sleep. The percentage of statistically correlated channels was weaker in N1 and N2 as compared with N3, and was strongest in regions that have previously been associated with autonomic processes, such as anterior hippocampus and insula. The anatomic distribution of HRV-HG correlations during sleep did not reproduce those usually observed with positron emission tomography (PET) or functional magnetic resonance imaging (fMRI) during waking. This study aims to characterize the relationship between autonomic tone and neuronal firing rate during sleep and further studies are needed to investigate finer temporal resolutions, denser coverages, and different frequency bands in both waking and sleep.

The interconnected causes and consequences of sleep in the brain

[Figure: see text].

Contribution of animal models toward understanding resting state functional connectivity

Functional connectivity, which reflects the spatial and temporal organization of intrinsic activity throughout the brain, is one of the most studied measures in human neuroimaging research. The noninvasive acquisition of resting state functional magnetic resonance imaging (rs-fMRI) allows the characterization of features designated as functional networks, functional connectivity gradients, and time-varying activity patterns that provide insight into the intrinsic functional organization of the brain and potential alterations related to brain dysfunction. Functional connectivity, hence, captures dimensions of the brain's activity that have enormous potential for both clinical and preclinical research. However, the mechanisms underlying functional connectivity have yet to be fully characterized, hindering interpretation of rs-fMRI studies. As in other branches of neuroscience, the identification of the neurophysiological processes that contribute to functional connectivity largely depends on research conducted on laboratory animals, which provide a platform where specific, multi-dimensional investigations that involve invasive measurements can be carried out. These highly controlled experiments facilitate the interpretation of the temporal correlations observed across the brain. Indeed, information obtained from animal experimentation to date is the basis for our current understanding of the underlying basis for functional brain connectivity. This review presents a compendium of some of the most critical advances in the field based on the efforts made by the animal neuroimaging community.

Internal State: Dynamic, Interconnected Communication Loops Distributed Across Body, Brain, and Time

Internal state profoundly alters perception and behavior. For example, a starved fly may approach and consume foods that it would otherwise find undesirable. A socially engaged newt may remain engaged in the presence of a predator, whereas a solitary newt would otherwise attempt to escape. Yet, the definition of internal state is fluid and ill-defined. As an interdisciplinary group of scholars spanning five career stages (from undergraduate to full professor) and six academic institutions, we came together in an attempt to provide an operational definition of internal state that could be useful in understanding the behavior and the function of nervous systems, at timescales relevant to the individual. In this perspective, we propose to define internal state through an integrative framework centered on dynamic and interconnected communication loops within and between the body and the brain. This framework is informed by a synthesis of historical and contemporary paradigms used by neurobiologists, ethologists, physiologists, and endocrinologists. We view internal state as composed of both spatially distributed networks (body-brain communication loops), and temporally distributed mechanisms that weave together neural circuits, physiology, and behavior. Given the wide spatial and temporal scales at which internal state operates-and therefore the broad range of scales at which it could be defined-we choose to anchor our definition in the body. Here we focus on studies that highlight body-to-brain signaling; body represented in endocrine signaling, and brain represented in sensory signaling. This integrative framework of internal state potentially unites the disparate paradigms often used by scientists grappling with body-brain interactions. We invite others to join us as we examine approaches and question assumptions to study the underlying mechanisms and temporal dynamics of internal state.

Naturalistic Stimuli: A Paradigm for Multi-Scale Functional Characterization of the Human Brain

Movies, audio stories, and virtual reality are increasingly used as stimuli for functional brain imaging. Such naturalistic paradigms are in sharp contrast to the tradition of experimental reductionism in neuroscience research. Being complex, dynamic, and diverse, naturalistic stimuli set up a more ecologically relevant condition and induce highly reproducible brain responses across a wide range of spatiotemporal scales. Here, we review recent technical advances and scientific findings on imaging the brain under naturalistic stimuli. Then we elaborate on the premise of using naturalistic paradigms for multi-scale, multi-modal, and high-throughput functional characterization of the human brain. We further highlight the growing potential of using deep learning models to infer neural information processing from brain responses to naturalistic stimuli. Lastly, we advocate large-scale collaborations to combine brain imaging and recording data across experiments, subjects, and labs that use the same set of naturalistic stimuli.

Revisiting the neurovascular unit

The brain is supplied by an elaborate vascular network that originates extracranially and reaches deep into the brain. The concept of the neurovascular unit provides a useful framework to investigate how neuronal signals regulate nearby microvessels to support the metabolic needs of the brain, but it does not consider the role of larger cerebral arteries and systemic vasoactive signals. Furthermore, the recently emerged molecular heterogeneity of cerebrovascular cells indicates that there is no prototypical neurovascular unit replicated at all levels of the vascular network. Here, we examine the cellular and molecular diversity of the cerebrovascular tree and the relative contribution of systemic and brain-intrinsic factors to neurovascular function. Evidence supports the concept of a 'neurovascular complex' composed of segmentally diverse functional modules that implement coordinated vascular responses to central and peripheral signals to maintain homeostasis of the brain. This concept has major implications for neurovascular regulation in health and disease and for brain imaging.

Decoding the brain state-dependent relationship between pupil dynamics and resting state fMRI signal fluctuation

Pupil dynamics serve as a physiological indicator of cognitive processes and arousal states of the brain across a diverse range of behavioral experiments. Pupil diameter changes reflect brain state fluctuations driven by neuromodulatory systems. Resting-state fMRI (rs-fMRI) has been used to identify global patterns of neuronal correlation with pupil diameter changes; however, the linkage between distinct brain state-dependent activation patterns of neuromodulatory nuclei with pupil dynamics remains to be explored. Here, we identified four clusters of trials with unique activity patterns related to pupil diameter changes in anesthetized rat brains. Going beyond the typical rs-fMRI correlation analysis with pupil dynamics, we decomposed spatiotemporal patterns of rs-fMRI with principal component analysis (PCA) and characterized the cluster-specific pupil-fMRI relationships by optimizing the PCA component weighting via decoding methods. This work shows that pupil dynamics are tightly coupled with different neuromodulatory centers in different trials, presenting a novel PCA-based decoding method to study the brain state-dependent pupil-fMRI relationship.

Physiological noise modeling in fMRI based on the pulsatile component of photoplethysmograph

The blood oxygenation level-dependent (BOLD) contrast mechanism allows the noninvasive monitoring of changes in deoxyhemoglobin content. As such, it is commonly used in functional magnetic resonance imaging (fMRI) to study brain activity since levels of deoxyhemoglobin are indirectly related to local neuronal activity through neurovascular coupling mechanisms. However, the BOLD signal is severely affected by physiological processes as well as motion. Due to this, several noise correction techniques have been developed to correct for the associated confounds. The present study focuses on cardiac pulsatility fMRI confounds, aiming to refine model-based techniques that utilize the photoplethysmograph (PPG) signal. Specifically, we propose a new technique based on convolution filtering, termed cardiac pulsatility model (CPM) and compare its performance with the cardiac-related RETROICOR (Card-RETROICOR), which is a technique commonly used to model fMRI fluctuations due to cardiac pulsatility. Further, we investigate whether variations in the amplitude of the PPG pulses (PPG-Amp) covary with variations in amplitude of pulse-related fMRI fluctuations, as well as with the systemic low frequency oscillations (SLFOs) component of the fMRI global signal (GS - defined as the mean signal across all gray matter voxels). Capitalizing on 3T fMRI data from the Human Connectome Project, CPM was found to explain a significantly larger fraction of the fMRI signal variance compared to Card-RETROICOR, particularly for subjects with larger heart rate variability during the scan. The amplitude of the fMRI pulse-related fluctuations did not covary with PPG-Amp; however, PPG-Amp explained significant variance in the GS that was not attributed to variations in heart rate or breathing patterns. Our results suggest that the proposed approach can model high-frequency fluctuations due to pulsation as well as low-frequency physiological fluctuations more accurately compared to model-based techniques commonly employed in fMRI studies.

Physiological and motion signatures in static and time-varying functional connectivity and their subject identifiability

Human brain connectivity yields significant potential as a noninvasive biomarker. Several studies have used fMRI-based connectivity fingerprinting to characterize individual patterns of brain activity. However, it is not clear whether these patterns mainly reflect neural activity or the effect of physiological and motion processes. To answer this question, we capitalize on a large data sample from the Human Connectome Project and rigorously investigate the contribution of the aforementioned processes on functional connectivity (FC) and time-varying FC, as well as their contribution to subject identifiability. We find that head motion, as well as heart rate and breathing fluctuations, induce artifactual connectivity within distinct resting-state networks and that they correlate with recurrent patterns in time-varying FC. Even though the spatiotemporal signatures of these processes yield above-chance levels in subject identifiability, removing their effects at the preprocessing stage improves identifiability, suggesting a neural component underpinning the inter-individual differences in connectivity.

The effect of body posture on resting-state functional connectivity

INTRODUCTION

An important but under-investigated confound of functional MRI (fMRI) is body posture. Although it is well-established that body position changes cerebral blood flow, the amount of cerebrospinal fluid in the brain, intracranial pressure, and even the firing rate of certain cell types, there is currently no study that directly examines its effect on fMRI measurements. Moreover, fMRI is typically done in a supine position, which often does not correspond to how these processes are performed in everyday settings.

METHODS

In this study, 20 healthy adults underwent resting-state fMRI under three body positions: supine, right lateral decubitus (RLD), and left lateral decubitus (LLD). We first investigated whether there were differences in overall organization of whole-brain connectivity between conditions using graph theory. Second, we examined whether functional connectivity of two most studied default mode network (DMN) seeds to the rest of the brain was altered as a function of body position.

RESULTS

Nonparametric statistical analyses revealed that global network measures differed among conditions, with the supine and LLD showing identical results compared to the RLD. There was decreased connectivity for DMN seeds in the RLD condition compared to the supine and LLD, but there were no significant differences between the latter two conditions.

DISCUSSION

Potential mechanisms underlying these alterations include gravity, changes in physiology, and body anatomy. Our results suggest that, compared to supine and LLD, the RLD position leads to changes in whole-brain and DMN connectivity. Finally, depending on the research question, combining imaging modalities that allow for more naturalistic settings can provide a better understanding of certain phenomena.

Reduced coupling between cerebrospinal fluid flow and global brain activity is linked to Alzheimer disease-related pathology

The glymphatic system plays an important role in clearing the amyloid-β (Aβ) and tau proteins that are closely linked to Alzheimer disease (AD) pathology. Glymphatic clearance, as well as Aβ accumulation, is highly dependent on sleep, but the sleep-dependent driving forces behind cerebrospinal fluid (CSF) movements essential to the glymphatic flux remain largely unclear. Recent studies have reported that widespread, high-amplitude spontaneous brain activations in the drowsy state and during sleep, which are shown as large global signal peaks in resting-state functional magnetic resonance imaging (rsfMRI), are coupled with CSF movements, suggesting their potential link to glymphatic flux and metabolite clearance. By analyzing multimodal data from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) project, here we showed that the coupling between the global fMRI signal and CSF influx is correlated with AD-related pathology, including various risk factors for AD, the severity of AD-related diseases, the cortical Aβ level, and cognitive decline over a 2-year follow-up. These results provide critical initial evidence for involvement of sleep-dependent global brain activity, as well as the associated physiological modulations, in the clearance of AD-related brain waste.

This study reveals strong coupling between the global fMRI signal and cerebrospinal fluid influx, finding that this is correlated with Alzheimer’s disease-related pathology, disease severity, and cognitive decline. This supports a link between spontaneous low-frequency brain dynamics and Alzheimer’s disease pathology, presumably due to their role in glymphatic clearance.

The Effect of CO2 on Resting-State Functional Connectivity: Isocapnia vs. Poikilocapnia

The normal variability in breath size and frequency results in breath-to-breath variability of end-tidal PCO₂ (P_ETCO₂), the measured variable, and arterial partial pressure of carbon dioxide (PaCO₂), the independent variable affecting cerebral blood flow (CBF). This study examines the effect of variability in PaCO₂ on the pattern of resting-state functional MRI (rs-fMRI) connectivity. A region of interest (ROI)-to-ROI and Seed-to-Voxel first-level bivariate correlation, hemodynamic response function (hrf)-weighted analysis for measuring rs-fMRI connectivity was performed during two resting-state conditions: (a) normal breathing associated with breath-to-breath variation in PaCO₂ (poikilocapnia), and (b) normal breathing with breath-to-breath variability of P_ETCO₂ dampened using sequential rebreathing (isocapnia). End-tidal PCO₂ (P_ETCO₂) was used as a measurable surrogate for fluctuations of PaCO₂. During poikilocapnia, enhanced functional connections were found between the cerebellum and inferior frontal and supramarginal gyrus (SG), visual cortex and occipital fusiform gyrus; and between the primary visual network (PVN) and the hippocampal formation. During isocapnia, these associations were not seen, rather enhanced functional connections were identified in the corticostriatal pathway between the putamen and intracalacarine cortex, supracalcarine cortex (SCC), and precuneus cortex. We conclude that vascular responses to variations in P_ETCO₂, account for at least some of the observed resting state synchronization of blood oxygenation level-dependent (BOLD) signals.

fMRI-based detection of alertness predicts behavioral response variability

Levels of alertness are closely linked with human behavior and cognition. However, while functional magnetic resonance imaging (fMRI) allows for investigating whole-brain dynamics during behavior and task engagement, concurrent measures of alertness (such as EEG or pupillometry) are often unavailable. Here, we extract a continuous, time-resolved marker of alertness from fMRI data alone. We demonstrate that this fMRI alertness marker, calculated in a short pre-stimulus interval, captures trial-to-trial behavioral responses to incoming sensory stimuli. In addition, we find that the prediction of both EEG and behavioral responses during the task may be accomplished using only a small fraction of fMRI voxels. Furthermore, we observe that accounting for alertness appears to increase the statistical detection of task-activated brain areas. These findings have broad implications for augmenting a large body of existing datasets with information about ongoing arousal states, enriching fMRI studies of neural variability in health and disease.

Cortical and subcortical hemodynamic changes during sleep slow waves in human light sleep

EEG slow waves, the hallmarks of NREM sleep are thought to be crucial for the regulation of several important processes, including learning, sensory disconnection and the removal of brain metabolic wastes. Animal research indicates that slow waves may involve complex interactions within and between cortical and subcortical structures. Conventional EEG in humans, however, has a low spatial resolution and is unable to accurately describe changes in the activity of subcortical and deep cortical structures. To overcome these limitations, here we took advantage of simultaneous EEG-fMRI recordings to map cortical and subcortical hemodynamic (BOLD) fluctuations time-locked to slow waves of light sleep. Recordings were performed in twenty healthy adults during an afternoon nap. Slow waves were associated with BOLD-signal increases in the posterior brainstem and in portions of thalamus and cerebellum characterized by preferential functional connectivity with limbic and somatomotor areas, respectively. At the cortical level, significant BOLD-signal decreases were instead found in several areas, including insula and somatomotor cortex. Specifically, a slow signal increase preceded slow-wave onset and was followed by a delayed, stronger signal decrease. Similar hemodynamic changes were found to occur at different delays across most cortical brain areas, mirroring the propagation of electrophysiological slow waves, from centro-frontal to inferior temporo-occipital cortices. Finally, we found that the amplitude of electrophysiological slow waves was positively related to the magnitude and inversely related to the delay of cortical and subcortical BOLD-signal changes. These regional patterns of brain activity are consistent with theoretical accounts of the functions of sleep slow waves.

Relationships between oxygen changes in the brain and periphery following physiological activation and the actions of heroin and cocaine

Using two-sensor electrochemical recordings in freely moving rats, we examined the relationship between physiological and drug-induced oxygen fluctuations in the brain and periphery. Animals chronically implanted with oxygen sensors in the nucleus accumbens (NAc) and subcutaneous (SC) space were subjected to several mildly arousing stimuli (sound, tail-pinch and social interaction) and intravenous injections of cocaine and heroin. Arousing stimuli induced rapid increases in NAc oxygen levels followed by and correlated with oxygen decreases in the SC space. Therefore, cerebral vasodilation that increases cerebral blood flow and oxygen entry into brain tissue results from both direct neuronal activation and peripheral vasoconstriction, which redistributes arterial blood from periphery to the brain. The latter factor could also explain a similar pattern of oxygen responses found in the substantia nigra reticulata, suggesting hyperoxia as a global phenomenon with minor structural differences during early time intervals following the stimulus onset. While arousing stimuli and cocaine induced similar oxygen responses in the brain and SC space, heroin induced a biphasic down-up brain oxygen fluctuation associated with a monophasic oxygen decrease in the SC space. Oxygen decreases occurred more rapidly and stronger in the SC space, reflecting a drop in blood oxygen levels due to respiratory depression.

Multimodal EEG-fMRI: advancing insight into large-scale human brain dynamics

Advances in the acquisition and analysis of functional magnetic resonance imaging (fMRI) data are revealing increasingly rich spatiotemporal structure across the human brain. Nonetheless, uncertainty surrounding the origins of fMRI hemodynamic signals, and in the link between large-scale fMRI patterns and ongoing functional states, presently limits the neurobiological conclusions one can draw from fMRI alone. Electroencephalography (EEG) provides complementary information about neural electrical activity and state change, and simultaneously acquiring EEG together with fMRI presents unique opportunities for studying large-scale brain activity and gaining more information from fMRI itself. Here, we discuss recent progress in the use of concurrent EEG-fMRI to enrich the investigation of neural and physiological states and clarify the origins of fMRI hemodynamic signals. Throughout, we outline perspectives on future directions and open challenges.

Best practices for fNIRS publications

The application of functional near-infrared spectroscopy (fNIRS) in the neurosciences has been expanding over the last 40 years. Today, it is addressing a wide range of applications within different populations and utilizes a great variety of experimental paradigms. With the rapid growth and the diversification of research methods, some inconsistencies are appearing in the way in which methods are presented, which can make the interpretation and replication of studies unnecessarily challenging. The Society for Functional Near-Infrared Spectroscopy has thus been motivated to organize a representative (but not exhaustive) group of leaders in the field to build a consensus on the best practices for describing the methods utilized in fNIRS studies. Our paper has been designed to provide guidelines to help enhance the reliability, repeatability, and traceability of reported fNIRS studies and encourage best practices throughout the community. A checklist is provided to guide authors in the preparation of their manuscripts and to assist reviewers when evaluating fNIRS papers.

Rude mechanicals in brain haemodynamics: non-neural actors that influence blood flow

Fluctuations in blood oxygenation and flow are widely used to infer brain activity during resting-state functional magnetic resonance imaging (fMRI). However, there are strong systemic and vascular contributions to resting-state signals that are unrelated to ongoing neural activity. Importantly, these non-neural contributions to haemodynamic signals (or 'rude mechanicals') can be as large as or larger than the neurally evoked components. Here, we review the two broad classes of drivers of these signals. One is systemic and is tied to fluctuations in external drivers such as heart rate and breathing, and the robust autoregulatory mechanisms that try to maintain a constant milieu in the brain. The other class comprises local, active fluctuations that appear to be intrinsic to vascular tissue and are likely similar to active local fluctuations seen in vasculature all over the body. In this review, we describe these non-neural fluctuations and some of the tools developed to correct for them when interpreting fMRI recordings. However, we also emphasize the links between these vascular fluctuations and brain physiology and point to ways in which fMRI measurements can be used to exploit such links to gain valuable information about neurovascular health and about internal brain states. This article is part of the theme issue 'Key relationships between non-invasive functional neuroimaging and the underlying neuronal activity'.

Reconstruction of respiratory variation signals from fMRI data

Functional MRI signals can be heavily influenced by systemic physiological processes in addition to local neural activity. For example, widespread hemodynamic fluctuations across the brain have been found to correlate with natural, low-frequency variations in the depth and rate of breathing over time. Acquiring peripheral measures of respiration during fMRI scanning not only allows for modeling such effects in fMRI analysis, but also provides valuable information for interrogating brain-body physiology. However, physiological recordings are frequently unavailable or have insufficient quality. Here, we propose a computational technique for reconstructing continuous low-frequency respiration volume (RV) fluctuations from fMRI data alone. We evaluate the performance of this approach across different fMRI preprocessing strategies. Further, we demonstrate that the predicted RV signals can account for similar patterns of temporal variation in resting-state fMRI data compared to measured RV fluctuations. These findings indicate that fluctuations in respiration volume can be extracted from fMRI alone, in the common scenario of missing or corrupted respiration recordings. The results have implications for enriching a large volume of existing fMRI datasets through retrospective addition of respiratory variations information.

Short-channel regression in functional near-infrared spectroscopy is more effective when considering heterogeneous scalp hemodynamics

Significance: The reliability of functional near-infrared spectroscopy (fNIRS) measurements is reduced by systemic physiology. Short-channel regression algorithms aim at removing systemic "noise" by subtracting the signal measured at a short source-detector separation (mainly scalp hemodynamics) from the one of a long separation (brain and scalp hemodynamics). In literature, incongruent approaches on the selection of the optimal regressor signal are reported based on different assumptions on scalp hemodynamics properties. Aim: We investigated the spatial and temporal distribution of scalp hemodynamics over the sensorimotor cortex and evaluated its influence on the effectiveness of short-channel regressions. Approach: We performed hand-grasping and resting-state experiments with five subjects, measuring with 16 optodes over sensorimotor areas, including eight 8-mm channels. We performed detailed correlation analyses of scalp hemodynamics and evaluated 180 hand-grasping and 270 simulated (overlaid on resting-state measurements) trials. Five short-channel regressor combinations were implemented with general linear models. Three were chosen according to literature, and two were proposed based on additional physiological assumptions [considering multiple short channels and their Mayer wave (MW) oscillations]. Results: We found heterogeneous hemodynamics in the scalp, coming on top of a global close-to-homogeneous behavior (correlation 0.69 to 0.92). The results further demonstrate that short-channel regression always improves brain activity estimates but that better results are obtained when heterogeneity is assumed. In particular, we highlight that short-channel regression is more effective when combining multiple scalp regressors and when MWs are additionally included. Conclusion: We shed light on the selection of optimal regressor signals for improving the removal of systemic physiological artifacts in fNIRS. We conclude that short-channel regression is most effective when assuming heterogeneous hemodynamics, in particular when combining spatial- and frequency-specific information. A better understanding of scalp hemodynamics and more effective short-channel regression will promote more accurate assessments of functional brain activity in clinical and research settings.

Vigilance Effects in Resting-State fMRI

Measures of resting-state functional magnetic resonance imaging (rsfMRI) activity have been shown to be sensitive to cognitive function and disease state. However, there is growing evidence that variations in vigilance can lead to pronounced and spatially widespread differences in resting-state brain activity. Unless properly accounted for, differences in vigilance can give rise to changes in resting-state activity that can be misinterpreted as primary cognitive or disease-related effects. In this paper, we examine in detail the link between vigilance and rsfMRI measures, such as signal variance and functional connectivity. We consider how state changes due to factors such as caffeine and sleep deprivation affect both vigilance and rsfMRI measures and review emerging approaches and methodological challenges for the estimation and interpretation of vigilance effects.

The Contribution of Functional Magnetic Resonance Imaging to the Understanding of the Effects of Acute Physical Exercise on Cognition

The fact that a single bout of acute physical exercise has a positive impact on cognition is well-established in the literature, but the neural correlates that underlie these cognitive improvements are not well understood. Here, the use of neuroimaging techniques, such as functional magnetic resonance imaging (fMRI), offers great potential, which is just starting to be recognized. This review aims at providing an overview of those studies that used fMRI to investigate the effects of acute physical exercises on cerebral hemodynamics and cognition. To this end, a systematic literature survey was conducted by two independent reviewers across five electronic databases. The search returned 668 studies, of which 14 studies met the inclusion criteria and were analyzed in this systematic review. Although the findings of the reviewed studies suggest that acute physical exercise (e.g., cycling) leads to profound changes in functional brain activation, the small number of available studies and the great variability in the study protocols limits the conclusions that can be drawn with certainty. In order to overcome these limitations, new, more well-designed trials are needed that (i) use a more rigorous study design, (ii) apply more sophisticated filter methods in fMRI data analysis, (iii) describe the applied processing steps of fMRI data analysis in more detail, and (iv) provide a more precise exercise prescription.

Resting-state "physiological networks"

Slow changes in systemic brain physiology can elicit large fluctuations in fMRI time series, which manifest as structured spatial patterns of temporal correlations between distant brain regions. Here, we investigated whether such "physiological networks"-sets of segregated brain regions that exhibit similar responses following slow changes in systemic physiology-resemble patterns associated with large-scale networks typically attributed to remotely synchronized neuronal activity. By analyzing a large group of subjects from the 3T Human Connectome Project (HCP) database, we demonstrate brain-wide and noticeably heterogenous dynamics tightly coupled to either respiratory variation or heart rate changes. We show, using synthesized data generated from physiological recordings across subjects, that these physiologically-coupled fluctuations alone can produce networks that strongly resemble previously reported resting-state networks, suggesting that, in some cases, the "physiological networks" seem to mimic the neuronal networks. Further, we show that such physiologically-relevant connectivity estimates appear to dominate the overall connectivity observations in multiple HCP subjects, and that this apparent "physiological connectivity" cannot be removed by the use of a single nuisance regressor for the entire brain (such as global signal regression) due to the clear regional heterogeneity of the physiologically-coupled responses. Our results challenge previous notions that physiological confounds are either localized to large veins or globally coherent across the cortex, therefore emphasizing the necessity to consider potential physiological contributions in fMRI-based functional connectivity studies. The rich spatiotemporal patterns carried by such "physiological" dynamics also suggest great potential for clinical biomarkers that are complementary to large-scale neuronal networks.