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.

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.].

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.

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 blood oxygenation level-dependent (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 qualitatively 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 preserve information about stimuli but are then gated off downstream in visual cortex.NEW & NOTEWORTHY When we close our eyes coarse luminance information is still accessible by the visual system. Using functional magnetic resonance imaging, we examined whether eyelid closure plays a unique role in visual processing. We discovered that while the LGN and V1 show equivalent responses when the eyes are open or closed, extrastriate cortex exhibited attenuated responses with eye closure. This suggests that when the eyes are closed, downstream visual processing is blind to this information.

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.

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.

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.

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.

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.

Model-based physiological noise removal in fast fMRI.

Recent improvements in the speed and sensitivity of fMRI acquisition techniques suggest that fast fMRI can be used to detect and precisely localize sub-second neural dynamics. This enhanced temporal resolution has enormous potential for neuroscientists. However, physiological noise poses a major challenge for the analysis of fast fMRI data. Physiological noise scales with sensitivity, and its autocorrelation structure is altered in rapidly sampled data, suggesting that new approaches are needed for physiological noise removal in fast fMRI. Existing strategies either rely on external physiological recordings, which can be noisy or difficult to collect, or employ data-driven approaches which make assumptions that may not hold true in fast fMRI. We created a statistical model of harmonic regression with autoregressive noise (HRAN) to estimate and remove cardiac and respiratory noise from the fMRI signal directly. This technique exploits the fact that cardiac and respiratory noise signals are fully sampled (rather than aliasing) when imaging at fast rates, allowing us to track and model physiology over time without requiring external physiological measurements. We then created a joint model of neural hemodynamics, and physiological and autocorrelated noise to more accurately remove noise. We first verified that HRAN accurately estimates cardiac and respiratory dynamics and that our model demonstrates goodness-of-fit in fast fMRI data. In task-driven data, we then demonstrated that HRAN is able to remove physiological noise while leaving the neural signal intact, thereby increasing detection of task-driven voxels. Finally, we established that in both simulations and fast fMRI data HRAN is able to improve statistical inferences as compared with gold-standard physiological noise removal techniques. In conclusion, we created a tool that harnesses the novel information in fast fMRI to remove physiological noise, enabling broader use of the technology to study human brain function.

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.

Fast fMRI can detect oscillatory neural activity in humans.

Oscillatory neural dynamics play an important role in the coordination of large-scale brain networks. High-level cognitive processes depend on dynamics evolving over hundreds of milliseconds, so measuring neural activity in this frequency range is important for cognitive neuroscience. However, current noninvasive neuroimaging methods are not able to precisely localize oscillatory neural activity above 0.2 Hz. Electroencephalography and magnetoencephalography have limited spatial resolution, whereas fMRI has limited temporal resolution because it measures vascular responses rather than directly recording neural activity. We hypothesized that the recent development of fast fMRI techniques, combined with the extra sensitivity afforded by ultra-high-field systems, could enable precise localization of neural oscillations. We tested whether fMRI can detect neural oscillations using human visual cortex as a model system. We detected small oscillatory fMRI signals in response to stimuli oscillating at up to 0.75 Hz within single scan sessions, and these responses were an order of magnitude larger than predicted by canonical linear models. Simultaneous EEG-fMRI and simulations based on a biophysical model of the hemodynamic response to neuronal activity suggested that the blood oxygen level-dependent response becomes faster for rapidly varying stimuli, enabling the detection of higher frequencies than expected. Accounting for phase delays across voxels further improved detection, demonstrating that identifying vascular delays will be of increasing importance with higher-frequency activity. These results challenge the assumption that the hemodynamic response is slow, and demonstrate that fMRI has the potential to map neural oscillations directly throughout the brain.

Stimulus-dependent hemodynamic response timing across the human subcortical-cortical visual pathway identified through high spatiotemporal resolution 7T fMRI.

Recent developments in fMRI acquisition techniques now enable fast sampling with whole-brain coverage, suggesting fMRI can be used to track changes in neural activity at increasingly rapid timescales. When images are acquired at fast rates, the limiting factor for fMRI temporal resolution is the speed of the hemodynamic response. Given that HRFs may vary substantially in subcortical structures, characterizing the speed of subcortical hemodynamic responses, and how the hemodynamic response shape changes with stimulus duration (i.e. the hemodynamic nonlinearity), is needed for designing and interpreting fast fMRI studies of these regions. We studied the temporal properties and nonlinearities of the hemodynamic response function (HRF) across the human subcortical visual system, imaging superior colliculus (SC), lateral geniculate nucleus of the thalamus (LGN) and primary visual cortex (V1) with high spatiotemporal resolution 7 Tesla fMRI. By presenting stimuli of varying durations, we mapped the timing and nonlinearity of hemodynamic responses in these structures at high spatiotemporal resolution. We found that the hemodynamic response is consistently faster and narrower in subcortical structures than in cortex. However, the nonlinearity in LGN is similar to that in cortex, with shorter duration stimuli eliciting larger and faster responses than would have been predicted by a linear model. Using oscillatory visual stimuli, we tested the frequency response in LGN and found that its BOLD response tracked high-frequency (0.5 Hz) oscillations. The LGN response magnitudes were comparable to V1, allowing oscillatory BOLD signals to be detected in LGN despite the small size of this structure. These results suggest that the increase in the speed and amplitude of the hemodynamic response when neural activity is brief may be the key physiological driver of fast fMRI signals, enabling detection of high-frequency oscillations with fMRI. We conclude that subcortical visual structures exhibit fast and nonlinear hemodynamic responses, and that these dynamics enable detection of fast BOLD signals even within small deep brain structures when imaging is performed at ultra-high field.

SAP102 regulates synaptic AMPAR function through a CNIH-2-dependent mechanism.

The postsynaptic density (PSD)-95-like, disk-large (DLG) membrane-associated guanylate kinase (PSD/DLG-MAGUK) family of proteins scaffold α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptor (AMPAR) complexes to the postsynaptic compartment and are postulated to orchestrate activity-dependent modulation of synaptic AMPAR functions. SAP102 is a key member of this family, present from early development, before PSD-95 and PSD-93, and throughout life. Here we investigate the role of SAP102 in synaptic transmission using a cell-restricted molecular replacement strategy, where SAP102 is expressed against the background of acute knockdown of endogenous PSD-95. We show that SAP102 rescues the decrease of AMPAR-mediated evoked excitatory postsynaptic currents (AMPAR eEPSCs) and AMPAR miniature EPSC (AMPAR mEPSC) frequency caused by acute knockdown of PSD-95. Further analysis of the mini events revealed that PSD-95-to-SAP102 replacement but not direct manipulation of PSD-95 increases the AMPAR mEPSC decay time. SAP102-mediated rescue of AMPAR eEPSCs requires AMPAR auxiliary subunit cornichon-2, whereas cornichon-2 knockdown did not affect PSD-95-mediated regulation of AMPAR eEPSC. Combining these observations, our data elucidate that PSD-95 and SAP102 differentially influence basic synaptic properties and synaptic current kinetics potentially via different AMPAR auxiliary subunits. NEW & NOTEWORTHY Synaptic scaffold proteins postsynaptic density (PSD)-95-like, disk-large (DLG) membrane-associated guanylate kinase (PSD-MAGUKs) regulate synaptic α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptor (AMPAR) function. However, the functional diversity among different PSD-MAGUKs remains to be categorized. We show that distinct from PSD-95, SAP102 increase the AMPAR synaptic current decay time, and the effect of SAP102 on synaptic AMPAR function requires the AMPAR auxiliary subunit cornichon-2. Our data suggest that PSD-MAGUKs target and modulate different AMPAR complexes to exert specific experience-dependent modification of the excitatory circuit.

Rapid fragmentation of neuronal networks at the onset of propofol-induced unconsciousness.

The neurophysiological mechanisms by which anesthetic drugs cause loss of consciousness are poorly understood. Anesthetic actions at the molecular, cellular, and systems levels have been studied in detail at steady states of deep general anesthesia. However, little is known about how anesthetics alter neural activity during the transition into unconsciousness. We recorded simultaneous multiscale neural activity from human cortex, including ensembles of single neurons, local field potentials, and intracranial electrocorticograms, during induction of general anesthesia. We analyzed local and global neuronal network changes that occurred simultaneously with loss of consciousness. We show that propofol-induced unconsciousness occurs within seconds of the abrupt onset of a slow (<1 Hz) oscillation in the local field potential. This oscillation marks a state in which cortical neurons maintain local patterns of network activity, but this activity is fragmented across both time and space. Local (<4 mm) neuronal populations maintain the millisecond-scale connectivity patterns observed in the awake state, and spike rates fluctuate and can reach baseline levels. However, neuronal spiking occurs only within a limited slow oscillation-phase window and is silent otherwise, fragmenting the time course of neural activity. Unexpectedly, we found that these slow oscillations occur asynchronously across cortex, disrupting functional connectivity between cortical areas. We conclude that the onset of slow oscillations is a neural correlate of propofol-induced loss of consciousness, marking a shift to cortical dynamics in which local neuronal networks remain intact but become functionally isolated in time and space.

Differential requirement for NMDAR activity in SAP97ß-mediated regulation of the number and strength of glutamatergic AMPAR-containing synapses.

PSD-95-like, disc-large (DLG) family membrane-associated guanylate kinase proteins (PSD/DLG-MAGUKs) are essential for regulating synaptic AMPA receptor (AMPAR) function and activity-dependent trafficking of AMPARs. Using a molecular replacement strategy to replace endogenous PSD-95 with SAP97ß, we show that the prototypic ß-isoform of the PSD-MAGUKs, SAP97ß, has distinct NMDA receptor (NMDAR)-dependent roles in regulating basic properties of AMPAR-containing synapses. SAP97ß enhances the number of AMPAR-containing synapses in an NMDAR-dependent manner, whereas its effect on the size of unitary synaptic response is not fully dependent on NMDAR activity. These effects contrast with those of PSD-95α, which increases both the number of AMPAR-containing synapses and the size of unitary synaptic responses, with or without NMDAR activity. Our results suggest that SAP97ß regulates synaptic AMPAR content by increasing surface expression of GluA1-containing AMPARs, whereas PSD-95α enhances synaptic AMPAR content presumably by increasing the synaptic scaffold capacity for synaptic AMPARs. Our approach delineates discrete effects of different PSD-MAGUKs on principal properties of glutamatergic synaptic transmission. Our results suggest that the molecular diversity of PSD-MAGUKs can provide rich molecular substrates for differential regulation of glutamatergic synapses in the brain.

Local cortical dynamics of burst suppression in the anaesthetized brain.

Burst suppression is an electroencephalogram pattern that consists of a quasi-periodic alternation between isoelectric 'suppressions' lasting seconds or minutes, and high-voltage 'bursts'. It is characteristic of a profoundly inactivated brain, occurring in conditions including hypothermia, deep general anaesthesia, infant encephalopathy and coma. It is also used in neurology as an electrophysiological endpoint in pharmacologically induced coma for brain protection after traumatic injury and during status epilepticus. Classically, burst suppression has been regarded as a 'global' state with synchronous activity throughout cortex. This assumption has influenced the clinical use of burst suppression as a way to broadly reduce neural activity. However, the extent of spatial homogeneity has not been fully explored due to the challenges in recording from multiple cortical sites simultaneously. The neurophysiological dynamics of large-scale cortical circuits during burst suppression are therefore not well understood. To address this question, we recorded intracranial electrocorticograms from patients who entered burst suppression while receiving propofol general anaesthesia. The electrodes were broadly distributed across cortex, enabling us to examine both the dynamics of burst suppression within local cortical regions and larger-scale network interactions. We found that in contrast to previous characterizations, bursts could be substantially asynchronous across the cortex. Furthermore, the state of burst suppression itself could occur in a limited cortical region while other areas exhibited ongoing continuous activity. In addition, we found a complex temporal structure within bursts, which recapitulated the spectral dynamics of the state preceding burst suppression, and evolved throughout the course of a single burst. Our observations imply that local cortical dynamics are not homogeneous, even during significant brain inactivation. Instead, cortical and, implicitly, subcortical circuits express seemingly different sensitivities to high doses of anaesthetics that suggest a hierarchy governing how the brain enters burst suppression, and emphasize the role of local dynamics in what has previously been regarded as a global state. These findings suggest a conceptual shift in how neurologists could assess the brain function of patients undergoing burst suppression. First, analysing spatial variation in burst suppression could provide insight into the circuit dysfunction underlying a given pathology, and could improve monitoring of medically-induced coma. Second, analysing the temporal dynamics within a burst could help assess the underlying brain state. This approach could be explored as a prognostic tool for recovery from coma, and for guiding treatment of status epilepticus. Overall, these results suggest new research directions and methods that could improve patient monitoring in clinical practice.

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.

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.