Distributed fMRI dynamics predict distinct EEG rhythms across sleep and wakefulness.

The brain exhibits rich oscillatory dynamics that vary across tasks and states, such as the EEG oscillations that define sleep. These oscillations play critical roles in cognition and arousal, but the brainwide mechanisms underlying them are not yet described. Using simultaneous EEG and fast fMRI in subjects drifting between sleep and wakefulness, we developed a machine learning approach to investigate which brainwide fMRI dynamics predict alpha (8-12 Hz) and delta (1-4 Hz) rhythms. We predicted moment-by-moment EEG power from fMRI activity in held-out subjects, and found that information about alpha power was represented by a remarkably small set of regions, segregated in two distinct networks linked to arousal and visual systems. Conversely, delta rhythms were diffusely represented on a large spatial scale across the cortex. These results identify distributed networks that predict delta and alpha rhythms, and establish a computational framework for investigating fMRI brainwide dynamics underlying EEG oscillations.

The temporal specificity of BOLD fMRI is systematically related to anatomical and vascular features of the human brain.

The ability to detect fast responses with functional MRI depends on the speed of hemodynamic responses to neural activity, because hemodynamic responses act as a temporal low-pass filter smoothing out rapid changes. However, hemodynamic responses (their shape and timing) are highly variable across the brain and across stimuli. This heterogeneity of responses implies that the temporal specificity of fMRI signals, or the ability of fMRI to preserve fast information, should also vary substantially across the cortex. In this work we investigated how local differences in hemodynamic response timing impact the temporal specificity of fMRI. We conducted our research using ultra-high field (7T) fMRI at high spatiotemporal resolution, using the primary visual cortex (V1) as a model area for investigation. We used visual stimuli oscillating at slow and fast frequencies to probe the temporal specificity of individual voxels. As expected, we identified substantial variability in temporal specificity, with some voxels preserving their responses to fast neural activity more effectively than others. We investigated which voxels had the highest temporal specificity and related those to anatomical and vascular features of V1. We found that low temporal specificity is only weakly explained by the presence of large veins or cerebral cortical depth. Notably, however, temporal specificity depended strongly on a voxel's position along the anterior-posterior anatomical axis of V1, with voxels within the calcarine sulcus being capable of preserving close to 25% of their amplitude as the frequency of stimulation increased from 0.05-Hz to 0.20-Hz, and voxels nearest to the occipital pole preserving less than 18%. These results indicate that detection biases in high-resolution fMRI will depend on the anatomical and vascular features of the area being imaged, and that these biases will differ depending on the timing of the underlying neuronal activity. Importantly, this spatial heterogeneity of temporal specificity suggests that it could be exploited to achieve higher specificity in some locations, and that tailored data analysis strategies may help improve the detection and interpretation of fast fMRI responses.

Current Understanding of the Anatomy, Physiology, and Magnetic Resonance Imaging of Neurofluids: Update From the 2022 "ISMRM Imaging Neurofluids Study group" Workshop in Rome.

Neurofluids is a term introduced to define all fluids in the brain and spine such as blood, cerebrospinal fluid, and interstitial fluid. Neuroscientists in the past millennium have steadily identified the several different fluid environments in the brain and spine that interact in a synchronized harmonious manner to assure a healthy microenvironment required for optimal neuroglial function. Neuroanatomists and biochemists have provided an incredible wealth of evidence revealing the anatomy of perivascular spaces, meninges and glia and their role in drainage of neuronal waste products. Human studies have been limited due to the restricted availability of noninvasive imaging modalities that can provide a high spatiotemporal depiction of the brain neurofluids. Therefore, animal studies have been key in advancing our knowledge of the temporal and spatial dynamics of fluids, for example, by injecting tracers with different molecular weights. Such studies have sparked interest to identify possible disruptions to neurofluids dynamics in human diseases such as small vessel disease, cerebral amyloid angiopathy, and dementia. However, key differences between rodent and human physiology should be considered when extrapolating these findings to understand the human brain. An increasing armamentarium of noninvasive MRI techniques is being built to identify markers of altered drainage pathways. During the three-day workshop organized by the International Society of Magnetic Resonance in Medicine that was held in Rome in September 2022, several of these concepts were discussed by a distinguished international faculty to lay the basis of what is known and where we still lack evidence. We envision that in the next decade, MRI will allow imaging of the physiology of neurofluid dynamics and drainage pathways in the human brain to identify true pathological processes underlying disease and to discover new avenues for early diagnoses and treatments including drug delivery. Evidence level: 1 Technical Efficacy: Stage 3.

Differential cortical and subcortical visual processing with eyes shut.

Closing our eyes largely shuts down our ability to see. That said, our eyelids still pass some light, allowing our visual system to coarsely process information about visual scenes, such as changes in luminance. However, the specific impact of eye closure on processing within the early visual system remains largely unknown. To understand how visual processing is modulated when eyes are shut, we used functional magnetic resonance imaging (fMRI) to measure responses to a flickering visual stimulus at high (100%) and low (10%) temporal contrasts, while participants viewed the stimuli with their eyes open or closed. Interestingly, we discovered that eye closure produced a qualitatively distinct pattern of effects across the visual thalamus and visual cortex. We found that with eyes open, low temporal contrast stimuli produced smaller responses, across the lateral geniculate nucleus (LGN), primary (V1) and extrastriate visual cortex (V2). However, with eyes closed, we discovered that the LGN and V1 maintained similar BOLD responses as the eyes open condition, despite the suppressed visual input through the eyelid. In contrast, V2 and V3 had strongly attenuated BOLD response when eyes were closed, regardless of temporal contrast. Our findings reveal a qualitative distinct pattern of visual processing when the eyes are closed - one that is not simply an overall attenuation, but rather reflects distinct responses across visual thalamocortical networks, wherein the earliest stages of processing preserves information about stimuli but is then gated off downstream in visual cortex.

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

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

Functional magnetic resonance imaging (fMRI) is a tool that can be used to non-invasively measure the activity of the human brain. Active parts of the brain require more oxygen, which increases blood flow to these areas. fMRI can detect these changes, and its signal reflects the coupling between brain activity and changes in blood flow. The mechanism that couples brain activity to blood flow is known as the 'hemodynamic response', and its timing varies across the brain. Therefore, to interpret fMRI signals correctly and use them to measure underlying brain activity, it is necessary to understand how the response changes across the brain. Current methods for probing hemodynamic response variation are either limited to specific brain regions or require patients to hold their breath ­ something not all groups of patients can do. To solve this problem, Bailes et al. investigated whether resting-state fMRI signals contain information about how the hemodynamic response changes across the brain. This information could then be used to better infer brain activity from fMRI measurements. The experiments showed that resting-state fMRI signals can be used to characterize and predict the timing of the hemodynamic response. Specifically, the frequencies in resting-state fMRI signals are impacted by changes in the hemodynamic response and can therefore be used to predict hemodynamic timing. Additionally, Bailes et al. showed that these predictions are better than those obtained in experiments requiring patients to hold their breath, which is the current gold standard. The findings also demonstrate that the information from the frequencies of resting-state fMRI signals should be interpreted carefully, as differences in these frequencies can have a non-neural origin. Bailes et al. propose a highly generalizable approach for mapping and predicting variations of the hemodynamic response across the whole brain. These findings provide insights into the time-related properties of fMRI signals that are crucial for accurate analyses. This will be of particular importance as the field moves towards fMRI studies focused on rapid neural dynamics and higher-level cognition.

The impact of the human thalamus on brain-wide information processing.

The thalamus is a small, bilateral structure in the diencephalon that integrates signals from many areas of the CNS. This critical anatomical position allows the thalamus to influence whole-brain activity and adaptive behaviour. However, traditional research paradigms have struggled to attribute specific functions to the thalamus, and it has remained understudied in the human neuroimaging literature. Recent advances in analytical techniques and increased accessibility to large, high-quality data sets have brought forth a series of studies and findings that (re-)establish the thalamus as a core region of interest in human cognitive neuroscience, a field that otherwise remains cortico-centric. In this Perspective, we argue that using whole-brain neuroimaging approaches to investigate the thalamus and its interaction with the rest of the brain is key for understanding systems-level control of information processing. To this end, we highlight the role of the thalamus in shaping a range of functional signatures, including evoked activity, interregional connectivity, network topology and neuronal variability, both at rest and during the performance of cognitive tasks.

Beyond the symptom: the biology of fatigue.

A workshop titled "Beyond the Symptom: The Biology of Fatigue" was held virtually September 27-28, 2021. It was jointly organized by the Sleep Research Society and the Neurobiology of Fatigue Working Group of the NIH Blueprint Neuroscience Research Program. For access to the presentations and video recordings, see: https://neuroscienceblueprint.nih.gov/about/event/beyond-symptom-biology-fatigue. The goals of this workshop were to bring together clinicians and scientists who use a variety of research approaches to understand fatigue in multiple conditions and to identify key gaps in our understanding of the biology of fatigue. This workshop summary distills key issues discussed in this workshop and provides a list of promising directions for future research on this topic. We do not attempt to provide a comprehensive review of the state of our understanding of fatigue, nor to provide a comprehensive reprise of the many excellent presentations. Rather, our goal is to highlight key advances and to focus on questions and future approaches to answering them.

The MotoNet: A 3 Tesla MRI-Conditional EEG Net with Embedded Motion Sensors.

We introduce a new electroencephalogram (EEG) net, which will allow clinicians to monitor EEG while tracking head motion. Motion during MRI limits patient scans, especially of children with epilepsy. EEG is also severely affected by motion-induced noise, predominantly ballistocardiogram (BCG) noise due to the heartbeat. METHODS: The MotoNet was built using polymer thick film (PTF) EEG leads and motion sensors on opposite sides in the same flex circuit. EEG/motion measurements were made with a standard commercial EEG acquisition system in a 3 Tesla (T) MRI. A Kalman filtering-based BCG correction tool was used to clean the EEG in healthy volunteers. RESULTS: MRI safety studies in 3 T confirmed the maximum heating below 1 °C. Using an MRI sequence with spatial localization gradients only, the position of the head was linearly correlated with the average motion sensor output. Kalman filtering was shown to reduce the BCG noise and recover artifact-clean EEG. CONCLUSIONS: The MotoNet is an innovative EEG net design that co-locates 32 EEG electrodes with 32 motion sensors to improve both EEG and MRI signal quality. In combination with custom gradients, the position of the net can, in principle, be determined. In addition, the motion sensors can help reduce BCG noise.

EEG-LLAMAS: A low-latency neurofeedback platform for artifact reduction in EEG-fMRI.

Simultaneous EEG-fMRI is a powerful multimodal technique for imaging the brain, but its use in neurofeedback experiments has been limited by EEG noise caused by the MRI environment. Neurofeedback studies typically require analysis of EEG in real time, but EEG acquired inside the scanner is heavily contaminated with ballistocardiogram (BCG) artifact, a high-amplitude artifact locked to the cardiac cycle. Although techniques for removing BCG artifacts do exist, they are either not suited to real-time, low-latency applications, such as neurofeedback, or have limited efficacy. We propose and validate a new open-source artifact removal software called EEG-LLAMAS (Low Latency Artifact Mitigation Acquisition Software), which adapts and advances existing artifact removal techniques for low-latency experiments. We first used simulations to validate LLAMAS in data with known ground truth. We found that LLAMAS performed better than the best publicly-available real-time BCG removal technique, optimal basis sets (OBS), in terms of its ability to recover EEG waveforms, power spectra, and slow wave phase. To determine whether LLAMAS would be effective in practice, we then used it to conduct real-time EEG-fMRI recordings in healthy adults, using a steady state visual evoked potential (SSVEP) task. We found that LLAMAS was able to recover the SSVEP in real time, and recovered the power spectra collected outside the scanner better than OBS. We also measured the latency of LLAMAS during live recordings, and found that it introduced a lag of less than 50 ms on average. The low latency of LLAMAS, coupled with its improved artifact reduction, can thus be effectively used for EEG-fMRI neurofeedback. A limitation of the method is its use of a reference layer, a piece of EEG equipment which is not commercially available, but can be assembled in-house. This platform enables closed-loop experiments which previously would have been prohibitively difficult, such as those that target short-duration EEG events, and is shared openly with the neuroscience community.

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

[This corrects the article DOI: 10.1371/journal.pbio.3002035.].

Characterizing brain dynamics during ketamine-induced dissociation and subsequent interactions with propofol using human intracranial neurophysiology.

Ketamine produces antidepressant effects in patients with treatment-resistant depression, but its usefulness is limited by its psychotropic side effects. Ketamine is thought to act via NMDA receptors and HCN1 channels to produce brain oscillations that are related to these effects. Using human intracranial recordings, we found that ketamine produces gamma oscillations in prefrontal cortex and hippocampus, structures previously implicated in ketamine's antidepressant effects, and a 3 Hz oscillation in posteromedial cortex, previously proposed as a mechanism for its dissociative effects. We analyzed oscillatory changes after subsequent propofol administration, whose GABAergic activity antagonizes ketamine's NMDA-mediated disinhibition, alongside a shared HCN1 inhibitory effect, to identify dynamics attributable to NMDA-mediated disinhibition versus HCN1 inhibition. Our results suggest that ketamine engages different neural circuits in distinct frequency-dependent patterns of activity to produce its antidepressant and dissociative sensory effects. These insights may help guide the development of brain dynamic biomarkers and novel therapeutics for depression.

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.

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

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

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.

Disruption of Brainstem Structural Connectivity in REM Sleep Behavior Disorder Using 7 Tesla Magnetic Resonance Imaging.

BACKGROUND: Isolated rapid eye movement (REM) sleep behavior disorder (iRBD) is one of the earliest manifestations of α synucleinopathies. Brainstem pathophysiology underlying REM sleep behavior disorder has been described in animal models, yet it is understudied in living humans because of the lack of an in vivo brainstem nuclei atlas and to the limited magnetic resonance imaging (MRI) sensitivity. OBJECTIVE: To investigate brainstem structural connectivity changes in iRBD patients by using an in vivo probabilistic brainstem nuclei atlas and 7 Tesla MRI. METHODS: Structural connectivity of 12 iRBD patients and 12 controls was evaluated by probabilistic tractography. Two-sided Wilcoxon rank-sum test was used to compare the structural connectivity indices across groups. RESULTS: In iRBD, we found impaired (Z = 2.6, P < 0.01) structural connectivity in 14 brainstem nuclei, including the connectivity between REM-on (eg, subcoeruleus [SubC]) and REM sleep muscle atonia (eg, medullary reticular formation) areas. CONCLUSIONS: The brainstem nuclei diagram of impaired connectivity in human iRBD expands animal models and is a promising tool to study and possibly assess prodromal synucleinopathy stages. © 2021 International Parkinson and Movement Disorder Society.

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.

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

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

The interconnected causes and consequences of sleep in the brain.

Sleep is essential for brain function in a surprisingly diverse set of ways. In the short term, lack of sleep leads to impaired memory and attention; in the longer term, it produces neurological dysfunction or even death. I discuss recent advances in understanding how sleep maintains the physiological health of the brain through interconnected systems of neuronal activity and fluid flow. The neural dynamics that appear during sleep are intrinsically coupled to its consequences for blood flow, cerebrospinal fluid dynamics, and waste clearance. Recognizing these linked causes and consequences of sleep has shed new light on why sleep is important for such disparate aspects of brain function.

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

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