Spatiotemporal evolution of the functional magnetic resonance imaging response to ultrashort stimuli

High Spatiotemporal Resolution Radial Encoding Single-Vessel fMRI

High-field preclinical functional MRI (fMRI) is enabled the high spatial resolution mapping of vessel-specific hemodynamic responses, that is single-vessel fMRI. In contrast to investigating the neuronal sources of the fMRI signal, single-vessel fMRI focuses on elucidating its vascular origin, which can be readily implemented to identify vascular changes relevant to vascular dementia or cognitive impairment. However, the limited spatial and temporal resolution of fMRI is hindered hemodynamic mapping of intracortical microvessels. Here, the radial encoding MRI scheme is implemented to measure BOLD signals of individual vessels penetrating the rat somatosensory cortex. Radial encoding MRI is employed to map cortical activation with a focal field of view (FOV), allowing vessel-specific functional mapping with 50 × 50 µm2 in-plane resolution at a 1 to 2 Hz sampling rate. Besides detecting refined hemodynamic responses of intracortical micro-venules, the radial encoding-based single-vessel fMRI enables the distinction of fMRI signals from vessel and peri-vessel voxels due to the different contribution of intravascular and extravascular effects.

Functional ultrasound imaging of stroke in awake rats

Anesthesia is a major confounding factor in preclinical stroke research as stroke rarely occurs in sedated patients. Moreover, anesthesia affects both brain functions and the stroke outcome acting as neurotoxic or protective agents. So far, no approaches were well suited to induce stroke while imaging hemodynamics along with simultaneous large-scale recording of brain functions in awake animals. For this reason, the first critical hours following the stroke insult and associated functional alteration remain poorly understood. Here, we present a strategy to investigate both stroke hemodynamics and stroke-induced functional alterations without the confounding effect of anesthesia, i.e., under awake condition. Functional ultrasound (fUS) imaging was used to continuously monitor variations in cerebral blood volume (CBV) in +65 brain regions/hemispheres for up to 3 hr after stroke onset. The focal cortical ischemia was induced using a chemo-thrombotic agent suited for permanent middle cerebral artery occlusion in awake rats and followed by ipsi- and contralesional whiskers stimulation to investigate on the dynamic of the thalamocortical functions. Early (0-3 hr) and delayed (day 5) fUS recording enabled to characterize the features of the ischemia (location, CBV loss), spreading depolarizations (occurrence, amplitude) and functional alteration of the somatosensory thalamocortical circuits. Post-stroke thalamocortical functions were affected at both early and later time points (0-3 hr and 5 days) after stroke. Overall, our procedure facilitates early, continuous, and chronic assessments of hemodynamics and cerebral functions. When integrated with stroke studies or other pathological analyses, this approach seeks to enhance our comprehension of physiopathologies towards the development of pertinent therapeutic interventions.

Layer-fMRI VASO with short stimuli and event-related designs at 7 T

Layers and columns are the dominant processing units in the human (neo)cortex at the mesoscopic scale. While the blood oxygenation dependent (BOLD) signal has a high detection sensitivity, it is biased towards unwanted signals from large draining veins at the cortical surface. The additional fMRI contrast of vascular space occupancy (VASO) has the potential to augment the neuroscientific interpretability of layer-fMRI results by means of capturing complementary information of locally specific changes in cerebral blood volume (CBV). Specifically, VASO is not subject to unwanted sensitivity amplifications of large draining veins. Because of constrained sampling efficiency, it has been mainly applied in combination with efficient block task designs and long trial durations. However, to study cognitive processes in neuroscientific contexts, or probe vascular reactivity, short stimulation periods are often necessary. Here, we developed a VASO acquisition procedure with a short acquisition period and sub-millimeter resolution. During visual event-related stimulation, we show reliable responses in visual cortices within a reasonable number of trials (∼20). Furthermore, the short TR and high spatial specificity of our VASO implementation enabled us to show differences in laminar reactivity and onset times. Finally, we explore the generalizability to a different stimulus modality (somatosensation). With this, we showed that CBV-sensitive VASO provides the means to capture layer-specific haemodynamic responses with high spatio-temporal resolution and is able to be used with event-related paradigms.

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.

The Current State of Functional MR Imaging for Trauma Prognostication

In this review, we discuss the basics of functional MRI (fMRI) techniques including task-based and resting state fMRI, and overview the major findings in patients with traumatic brain injury. We summarize the studies that have longitudinally evaluated the changes in brain connectivity and task-related activation in trauma patients during different phases of trauma. We discuss how these data may potentially be used for prognostication, treatment planning, or monitoring and management of trauma patients.

MP-PCA denoising of fMRI time-series data can lead to artificial activation "spreading"

MP-PCA denoising has become the method of choice for denoising MRI data since it provides an objective threshold to separate the signal components from unwanted thermal noise components. In rodents, thermal noise in the coils is an important source of noise that can reduce the accuracy of activation mapping in fMRI. Further confounding this problem, vendor data often contains zero-filling and other post-processing steps that may violate MP-PCA assumptions. Here, we develop an approach to denoise vendor data and assess activation "spreading" caused by MP-PCA denoising in rodent task-based fMRI data. Data was obtained from N = 3 mice using conventional multislice and ultrafast acquisitions (1 s and 50 ms temporal resolution, respectively), during visual stimulation. MP-PCA denoising produced SNR gains of 64% and 39% and Fourier Spectral Amplitude (FSA) increases in BOLD maps of 9% and 7% for multislice and ultrafast data, respectively, when using a small [2 2] denoising window. Larger windows provided higher SNR and FSA gains with increased spatial extent of activation that may or may not represent real activation. Simulations showed that MP-PCA denoising can incur activation "spreading" with increased false positive rate and smoother functional maps due to local "bleeding" of principal components, and that the optimal denoising window for improved specificity of functional mapping, based on Dice score calculations, depends on the data's tSNR and functional CNR. This "spreading" effect applies also to another recently proposed low-rank denoising method (NORDIC), although to a lesser degree. Our results bode well for enhancing spatial and/or temporal resolution in future fMRI work, while taking into account the sensitivity/specificity trade-offs of low-rank denoising methods.

Baseline oxygen consumption decreases with cortical depth

The cerebral cortex is organized in cortical layers that differ in their cellular density, composition, and wiring. Cortical laminar architecture is also readily revealed by staining for cytochrome oxidase-the last enzyme in the respiratory electron transport chain located in the inner mitochondrial membrane. It has been hypothesized that a high-density band of cytochrome oxidase in cortical layer IV reflects higher oxygen consumption under baseline (unstimulated) conditions. Here, we tested the above hypothesis using direct measurements of the partial pressure of O2 (pO2) in cortical tissue by means of 2-photon phosphorescence lifetime microscopy (2PLM). We revisited our previously developed method for extraction of the cerebral metabolic rate of O2 (CMRO2) based on 2-photon pO2 measurements around diving arterioles and applied this method to estimate baseline CMRO2 in awake mice across cortical layers. To our surprise, our results revealed a decrease in baseline CMRO2 from layer I to layer IV. This decrease of CMRO2 with cortical depth was paralleled by an increase in tissue oxygenation. Higher baseline oxygenation and cytochrome density in layer IV may serve as an O2 reserve during surges of neuronal activity or certain metabolically active brain states rather than reflecting baseline energy needs. Our study provides to our knowledge the first quantification of microscopically resolved CMRO2 across cortical layers as a step towards better understanding of brain energy metabolism.

Separating neuronal activity and systemic low-frequency oscillation related BOLD responses at nodes of the default mode network during resting-state fMRI with multiband excitation echo-planar imaging

Functional magnetic resonance imaging (fMRI) evaluates brain activity using blood oxygenation level-dependent (BOLD) contrast. Resting-state fMRI (rsfMRI) examines spontaneous brain function using BOLD in the absence of a task, and the default mode network (DMN) has been identified from that. The DMN is a set of nodes within the brain that appear to be active and in communication when the subject is in an awake resting state. In addition to signal changes related to neural activity, it is thought that the BOLD signal may be affected by systemic low-frequency oscillations (SysLFOs) that are non-neuronal in source and likely propagate throughout the brain to arrive at different regions at different times. However, it may be difficult to distinguish between the response due to neuronal activity and the arrival of a SysLFO in specific regions. Conventional single-shot EPI (Conv) acquisition requires a longish repetition time, but faster image acquisition has recently become possible with multiband excitation EPI (MB). In this study, we evaluated the time-lag between nodes of the DMN using both Conv and MB protocols to determine whether it is possible to distinguish between neuronal activity and SysLFO related responses during rsfMRI. While the Conv protocol data suggested that SysLFOs substantially influence the apparent time-lag of neuronal activity, the MB protocol data implied that the effects of SysLFOs and neuronal activity on the BOLD response may be separated. Using a higher time-resolution acquisition for rsfMRI might help to distinguish neuronal activity induced changes to the BOLD response from those induced by non-neuronal sources.

Temporal stability of the hemodynamic response function across the majority of human cerebral cortex

The hemodynamic response function (HRF) measured with functional magnetic resonance imaging is generated by vascular and metabolic responses evoked by brief (<4 s) stimuli. It is known that the human HRF varies across cortex, between subjects, with stimulus paradigms, and even between different measurements in the same cortical location. However, our results demonstrate that strong HRFs are remarkably repeatable across sessions separated by time intervals up to 3 months. In this study, a multisensory stimulus was used to activate and measure the HRF across the majority of cortex (>70%, with lesser reliability observed in some areas of prefrontal cortex). HRFs were measured with high spatial resolution (2‐mm voxels) in central gray matter to minimize variations caused by partial‐volume effects. HRF amplitudes and temporal dynamics were highly repeatable across four sessions in 20 subjects. Positive and negative HRFs were consistently observed across sessions and subjects. Negative HRFs were generally weaker and, thus, more variable than positive HRFs. Statistical measurements showed that across‐session variability is highly correlated to the variability across events within a session; these measurements also indicated a normal distribution of variability across cortex. The overall repeatability of the HRFs over long time scales generally supports the long‐term use of event‐related functional magnetic resonance imaging protocols.

Functional Ultrasound Neuroimaging

Functional ultrasound (fUS) is a neuroimaging method that uses ultrasound to track changes in cerebral blood volume as an indirect readout of neuronal activity at high spatiotemporal resolution. fUS is capable of imaging head-fixed or freely behaving rodents and of producing volumetric images of the entire mouse brain. It has been applied to many species, including primates and humans. Now that fUS is reaching maturity, it is being adopted by the neuroscience community. However, the nature of the fUS signal and the different implementations of fUS are not necessarily accessible to nonspecialists. This review aims to introduce these ultrasound concepts to all neuroscientists. We explain the physical basis of the fUS signal and the principles of the method, present the state of the art of its hardware implementation, and give concrete examples of current applications in neuroscience. Finally, we suggest areas for improvement during the next few years.

Auditory Target Detection Enhances Visual Processing and Hippocampal Functional Connectivity

Though dividing one's attention between two input streams typically impairs performance, detecting a behaviorally relevant stimulus can sometimes enhance the encoding of unrelated information presented at the same time. Previous research has shown that selection of this kind boosts visual cortical activity and memory for concurrent items. An important unanswered question is whether such effects are reflected in processing quality and functional connectivity in visual regions and in the hippocampus. In this fMRI study, participants were asked to memorize a stream of naturalistic images and press a button only when they heard a predefined target tone (400 or 1,200 Hz, counterbalanced). Images could be presented with a target tone, with a distractor tone, or without a tone. Auditory target detection increased activity throughout the ventral visual cortex but lowered it in the hippocampus. Enhancements in functional connectivity between the ventral visual cortex and the hippocampus were also observed following auditory targets. Multi-voxel pattern classification of image category was more accurate on target tone trials than on distractor and no tone trials in the fusiform gyrus and parahippocampal gyrus. This effect was stronger in visual cortical clusters whose activity was more correlated with the hippocampus on target tone than on distractor tone trials. In agreement with accounts suggesting that subcortical noradrenergic influences play a role in the attentional boost effect, auditory target detection also caused an increase in locus coeruleus activity and phasic pupil responses. These findings outline a network of cortical and subcortical regions that are involved in the selection and processing of information presented at behaviorally relevant moments.

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.

A line through the brain: implementation of human line-scanning at 7T for ultra-high spatiotemporal resolution fMRI

Functional magnetic resonance imaging (fMRI) is a widely used tool in neuroscience to detect neurally evoked responses, e.g. the blood oxygenation level-dependent (BOLD) signal. Typically, BOLD fMRI has millimeter spatial resolution and temporal resolution of one to few seconds. To study the sub-millimeter structures and activity of the cortical gray matter, the field needs an fMRI method with high spatial and temporal resolution. Line-scanning fMRI achieves very high spatial resolution and high sampling rate, at the cost of a sacrifice in volume coverage. Here, we present a human line-scanning implementation on a 7T MRI system. First, we investigate the quality of the saturation pulses that suppress MR signal outside the line. Second, we established the best coil combination for reconstruction. Finally, we applied the line-scanning method in the occipital lobe during a visual stimulation task, showing BOLD responses along cortical depth, every 250 µm with a 200 ms repetition time (TR). We found a good correspondence of t-statistics values with 2D gradient-echo echo planar imaging (GE-EPI) BOLD fMRI data with the same temporal resolution and voxel volume (R = 0.6 ± 0.2). In summary, we demonstrate the feasibility of line-scanning in humans and this opens line-scanning fMRI for applications in cognitive and clinical neuroscience.

Early fMRI responses to somatosensory and optogenetic stimulation reflect neural information flow

fMRI has revolutionized how neuroscientists investigate human brain functions and networks. To further advance understanding of brain functions, identifying the direction of information flow, such as thalamocortical versus corticothalamic projections, is critical. Because the early hemodynamic response at microvessels near active neurons can be detected by ultrahigh field fMRI, we propose using the onset times of fMRI responses to discern the information flow. This approach was confirmed by observing the ultrahigh spatiotemporal resolution BOLD fMRI responses to bottom-up somatosensory stimulation and top-down optogenetic stimulation of the primary motor cortex in anesthetized mice. Because ultrahigh field MRI is increasingly available, ultrahigh spatiotemporal fMRI will significantly facilitate the investigation of functional circuits in humans.

Blood oxygenation level–dependent (BOLD) functional magnetic resonance imaging (fMRI) has been widely used to localize brain functions. To further advance understanding of brain functions, it is critical to understand the direction of information flow, such as thalamocortical versus corticothalamic projections. For this work, we performed ultrahigh spatiotemporal resolution fMRI at 15.2 T of the mouse somatosensory network during forepaw somatosensory stimulation and optogenetic stimulation of the primary motor cortex (M1). Somatosensory stimulation induced the earliest BOLD response in the ventral posterolateral nucleus (VPL), followed by the primary somatosensory cortex (S1) and then M1 and posterior thalamic nucleus. Optogenetic stimulation of excitatory neurons in M1 induced the earliest BOLD response in M1, followed by S1 and then VPL. Within S1, the middle cortical layers responded to somatosensory stimulation earlier than the upper or lower layers, whereas the upper cortical layers responded earlier than the other two layers to optogenetic stimulation in M1. The order of early BOLD responses was consistent with the canonical understanding of somatosensory network connections and cannot be explained by regional variabilities in the hemodynamic response functions measured using hypercapnic stimulation. Our data demonstrate that early BOLD responses reflect the information flow in the mouse somatosensory network, suggesting that high-field fMRI can be used for systems-level network analyses.

Neuroplasticity-driven timing modulations revealed by ultrafast functional magnetic resonance imaging

Detecting neuroplasticity in global brain circuits in vivo is key for understanding myriad processes such as memory, learning, and recovery from injury. Functional Magnetic Resonance Imaging (fMRI) is instrumental for such in vivo mappings, yet it typically relies on mapping changes in spatial extent of activation or via signal amplitude modulations, whose interpretation can be highly ambiguous. Importantly, a central aspect of neuroplasticity involves modulation of neural activity timing properties. We thus hypothesized that this temporal dimension could serve as a new marker for neuroplasticity. To detect fMRI signals more associated with the underlying neural dynamics, we developed an ultrafast fMRI (ufMRI) approach facilitating high spatiotemporal sensitivity and resolution in distributed neural pathways. When neuroplasticity was induced in the mouse visual pathway via dark rearing, ufMRI indeed mapped temporal modulations in the entire visual pathway. Our findings therefore suggest a new dimension for exploring neuroplasticity in vivo.

Dynamics of the cerebral blood flow response to brief neural activity in human visual cortex

The blood oxygen-level dependent (BOLD) functional magnetic resonance imaging (fMRI) signal depends on an interplay of cerebral blood flow (CBF), oxygen metabolism, and cerebral blood volume. Despite wide usage of BOLD fMRI, it is not clear how these physiological components create the BOLD signal. Here, baseline CBF and its dynamics evoked by a brief stimulus (2 s) in human visual cortex were measured at 3T. We found a stereotypical CBF response: immediate increase, rising to a peak a few second after the stimulus, followed by a significant undershoot. The BOLD hemodynamic response function (HRF) was also measured in the same session. Strong correlations between HRF and CBF peak responses indicate that the flow responses evoked by neural activation in nearby gray matter drive the early HRF. Remarkably, peak CBF and HRF were also strongly modulated by baseline perfusion. The CBF undershoot was reliable and significantly correlated with the HRF undershoot. However, late-time dynamics of the HRF and CBF suggest that oxygen metabolism can also contribute to the HRF undershoot. Combined measurement of the CBF and HRF for brief neural activation is a useful tool to understand the temporal dynamics of neurovascular and neurometabolic coupling.

Spatial and temporal patterns of nitric oxide diffusion and degradation drive emergent cerebrovascular dynamics

Nitric oxide (NO) is a gaseous signaling molecule that plays an important role in neurovascular coupling. NO produced by neurons diffuses into the smooth muscle surrounding cerebral arterioles, driving vasodilation. However, the rate of NO degradation in hemoglobin is orders of magnitude higher than in brain tissue, though how this might impact NO signaling dynamics is not completely understood. We used simulations to investigate how the spatial and temporal patterns of NO generation and degradation impacted dilation of a penetrating arteriole in cortex. We found that the spatial location of NO production and the size of the vessel both played an important role in determining its responsiveness to NO. The much higher rate of NO degradation and scavenging of NO in the blood relative to the tissue drove emergent vascular dynamics. Large vasodilation events could be followed by post-stimulus constrictions driven by the increased degradation of NO by the blood, and vasomotion-like 0.1-0.3 Hz oscillations could also be generated. We found that these dynamics could be enhanced by elevation of free hemoglobin in the plasma, which occurs in diseases such as malaria and sickle cell anemia, or following blood transfusions. Finally, we show that changes in blood flow during hypoxia or hyperoxia could be explained by altered NO degradation in the parenchyma. Our simulations suggest that many common vascular dynamics may be emergent phenomena generated by NO degradation by the blood or parenchyma.

Transfer functions linking neural calcium to single voxel functional ultrasound signal

Functional ultrasound imaging (fUS) is an emerging technique that detects changes of cerebral blood volume triggered by brain activation. Here, we investigate the extent to which fUS faithfully reports local neuronal activation by combining fUS and two-photon microscopy (2PM) in a co-registered single voxel brain volume. Using a machine-learning approach, we compute and validate transfer functions between dendritic calcium signals of specific neurons and vascular signals measured at both microscopic (2PM) and mesoscopic (fUS) levels. We find that transfer functions are robust across a wide range of stimulation paradigms and animals, and reveal a second vascular component of neurovascular coupling upon very strong stimulation. We propose that transfer functions can be considered as reliable quantitative reporters to follow neurovascular coupling dynamics.

Neurovascular coupling refers to changes in cerebral blood flow in response to neuronal stimulation, but to what extent this change can report neuronal activation is not known. Here the authors develop transfer functions between neural calcium signals and functional ultrasound changes in blood volume in co-registered single voxel brain volumes.

Effects of route of administration on oxytocin-induced changes in regional cerebral blood flow in humans

Could nose-to-brain pathways mediate the effects of peptides such as oxytocin (OT) on brain physiology when delivered intranasally? We address this question by contrasting two methods of intranasal administration (a standard nasal spray, and a nebulizer expected to improve OT deposition in nasal areas putatively involved in direct nose-to-brain transport) to intravenous administration in terms of effects on regional cerebral blood flow during two hours post-dosing. We demonstrate that OT-induced decreases in amygdala perfusion, a key hub of the OT central circuitry, are explained entirely by OT increases in systemic circulation following both intranasal and intravenous OT administration. Yet we also provide robust evidence confirming the validity of the intranasal route to target specific brain regions. Our work has important translational implications and demonstrates the need to carefully consider the method of administration in our efforts to engage specific central oxytocinergic targets for the treatment of neuropsychiatric disorders.

Cerebral oxygenation during locomotion is modulated by respiration

In the brain, increased neural activity is correlated with increases of cerebral blood flow and tissue oxygenation. However, how cerebral oxygen dynamics are controlled in the behaving animal remains unclear. We investigated to what extent cerebral oxygenation varies during locomotion. We measured oxygen levels in the cortex of awake, head-fixed mice during locomotion using polarography, spectroscopy, and two-photon phosphorescence lifetime measurements of oxygen sensors. We find that locomotion significantly and globally increases cerebral oxygenation, specifically in areas involved in locomotion, as well as in the frontal cortex and the olfactory bulb. The oxygenation increase persists when neural activity and functional hyperemia are blocked, occurred both in the tissue and in arteries feeding the brain, and is tightly correlated with respiration rate and the phase of respiration cycle. Thus, breathing rate is a key modulator of cerebral oxygenation and should be monitored during hemodynamic imaging, such as in BOLD fMRI.

BOLD hemodynamic response function changes significantly with healthy aging

Functional magnetic resonance imaging (fMRI) has been used to infer age-differences in neural activity from the hemodynamic response function (HRF) that characterizes the blood-oxygen-level-dependent (BOLD) signal over time. BOLD literature in healthy aging lacks consensus in age-related HRF changes, the nature of those changes, and their implications for measurement of age differences in brain function. Between-study discrepancies could be due to small sample sizes, analysis techniques, and/or physiologic mechanisms. We hypothesize that, with large sample sizes and minimal analysis assumptions, age-related changes in HRF parameters could reflect alterations in one or more components of the neural-vascular coupling system. To assess HRF changes in healthy aging, we analyzed the large population-derived dataset from the Cambridge Center for Aging and Neuroscience (CamCAN) study (Shafto et al., 2014). During scanning, 74 younger (18-30 years of age) and 173 older participants (54-74 years of age) viewed two checkerboards to the left and right of a central fixation point, simultaneously heard a binaural tone, and responded via right index finger button-press. To assess differences in the shape of the HRF between younger and older groups, HRFs were estimated using FMRIB's Linear Optimal Basis Sets (FLOBS) to minimize a priori shape assumptions. Group mean HRFs were different between younger and older groups in auditory, visual, and motor cortices. Specifically, we observed increased time-to-peak and decreased peak amplitude in older compared to younger adults in auditory, visual, and motor cortices. Changes in the shape and timing of the HRF in healthy aging, in the absence of performance differences, support our hypothesis of age-related changes in the neural-vascular coupling system beyond neural activity alone. More precise interpretations of HRF age-differences can be formulated once these physiologic factors are disentangled and measured separately.

Increased Resting Hippocampal and Basal Ganglia Perfusion in People at Ultra High Risk for Psychosis: Replication in a Second Cohort

We recently reported that resting hippocampal, basal ganglia and midbrain perfusion is elevated in people at ultra high risk (UHR) for psychosis. The present study sought to replicate our previous finding in an independent UHR cohort, and examined the relationship between resting perfusion in these regions, psychosis and depression symptoms, and traumatic experiences in childhood. Pseudo-Continuous Arterial Spin Labelling (p-CASL) imaging was used to measure resting cerebral blood flow (rCBF) in 77 UHR for psychosis individuals and 25 healthy volunteers in a case-control design. UHR participants were recruited from clinical early detection services at 3 sites in the South of England. Symptoms levels were assessed using the Comprehensive Assessment of At Risk Mental States (CAARMS), the Hamilton Depression Scale (HAM-D), and childhood trauma was assessed retrospectively using the Childhood Trauma Questionnaire (CTQ). Right hippocampal and basal ganglia rCBF were significantly increased in UHR subjects compared to controls, partially replicating our previous finding in an independent cohort. In UHR participants, positive symptoms were positively correlated with rCBF in the right pallidum. CTQ scores were positively correlated with rCBF values in the bilateral hippocampus and negatively associated with rCBF in the left prefrontal cortex. Elevated resting hippocampal and basal ganglia activity appears to be a consistent finding in individuals at high risk for psychosis, consistent with data from preclinical models of the disorder. The association with childhood trauma suggests that its influence on the risk of psychosis may be mediated through an effect on hippocampal function.

The brain's hemodynamic response function rapidly changes under acute psychosocial stress in association with genetic and endocrine stress response markers

Ample evidence links dysregulation of the stress response to the risk for psychiatric disorders. However, we lack an integrated understanding of mechanisms that are adaptive during the acute stress response but potentially pathogenic when dysregulated. One mechanistic link emerging from rodent studies is the interaction between stress effectors and neurovascular coupling, a process that adjusts cerebral blood flow according to local metabolic demands. Here, using task-related fMRI, we show that acute psychosocial stress rapidly impacts the peak latency of the hemodynamic response function (HRF-PL) in temporal, insular, and prefrontal regions in two independent cohorts of healthy humans. These latency effects occurred in the absence of amplitude effects and were moderated by regulatory genetic variants of KCNJ2, a known mediator of the effect of stress on vascular responsivity. Further, hippocampal HRF-PL correlated with both cortisol response and genetic variants that influence the transcriptional response to stress hormones and are associated with risk for major depression. We conclude that acute stress modulates hemodynamic response properties as part of the physiological stress response and suggest that HRF indices could serve as endophenotype of stress-related disorders.

Effective Connectivity Measured Using Optogenetically Evoked Hemodynamic Signals Exhibits Topography Distinct from Resting State Functional Connectivity in the Mouse

Brain connectomics has expanded from histological assessment of axonal projection connectivity (APC) to encompass resting state functional connectivity (RS-FC). RS-FC analyses are efficient for whole-brain mapping, but attempts to explain aspects of RS-FC (e.g., interhemispheric RS-FC) based on APC have been only partially successful. Neuroimaging with hemoglobin alone lacks specificity for determining how activity in a population of cells contributes to RS-FC. Wide-field mapping of optogenetically defined connectivity could provide insights into the brain's structure-function relationship. We combined optogenetics with optical intrinsic signal imaging to create an efficient, optogenetic effective connectivity (Opto-EC) mapping assay. We examined EC patterns of excitatory neurons in awake, Thy1-ChR2 transgenic mice. These Thy1-based EC (Thy1-EC) patterns were evaluated against RS-FC over the cortex. Compared to RS-FC, Thy1-EC exhibited increased spatial specificity, reduced interhemispheric connectivity in regions with strong RS-FC, and appreciable connection strength asymmetry. Comparing the topography of Thy1-EC and RS-FC patterns to maps of APC revealed that Thy1-EC more closely resembled APC than did RS-FC. The more general method of Opto-EC mapping with hemoglobin can be determined for 100 sites in single animals in under an hour, and is amenable to other neuroimaging modalities. Opto-EC mapping represents a powerful strategy for examining evolving connectivity-related circuit plasticity.

Impulse response timing differences in BOLD and CBV weighted fMRI

Recent advances in BOLD fMRI scan techniques have substantially improved spatial and temporal resolution, currently reaching to sub-millimeter and sub-second levels respectively. Unfortunately, there remain physiological barriers that prevent achieving this resolution in practice. BOLD contrast relies on the hemodynamic response to neuronal activity, whose associated cerebral blood oxygenation (CBO) changes may spread over several millimeters and last several seconds. Recent reports have suggested that significant improvements may be possible with cerebral blood volume (CBV)-weighted fMRI, which highlights the CBV changes rather than the BOLD changes associated with the hemodynamic response. Nevertheless, quantitative comparisons between CBV and BOLD are sparse, in particular regarding their temporal characteristics in human brain. To address this, we studied a cohort of subjects that received injection of ferumoxytol, an intravascular iron-oxide based contrast agent that introduces strong CBV contrast. An event-related visual stimulus paradigm was used to compare the impulse response (IR) for CBV and BOLD contrast, obtained with and without ferumoxytol, respectively. Experiments performed at 7 T (n = 5) at 1.2-1.5 mm spatial and 1 s temporal resolution showed that the onset time and time-to-peak of the CBV IR averaged 0.8 and 3.5 s respectively, both 0.6 s shorter than the BOLD IR. While significant, these improvements are relatively small and not expected to lead to practical advantages for the extraction of temporal information about neural activity. Nonlinearities in the observed IR were also compared and found to be similar between the CBV and BOLD, indicating that these are likely not caused by a 'ceiling' effect in the CBO response, but rather support a previously proposed model of vascular compliance, in which changes in vascular tone elicited by a preceding stimulus affect the IR.

Investigation of the BOLD and CBV fMRI responses to somatosensory stimulation in awake marmosets (Callithrix jacchus)

Understanding the spatiotemporal features of the hemodynamic response function (HRF) to brain stimulation is essential for the correct application of neuroimaging methods to study brain function. Here, we investigated the spatiotemporal evolution of the blood oxygen level-dependent (BOLD) and cerebral blood volume (CBV) HRF in conscious, awake marmosets (Callithrix jacchus), a New World non-human primate with a lissencephalic brain and with growing use in biomedical research. The marmosets were acclimatized to head fixation and placed in a 7-T magnetic resonance imaging (MRI) scanner. Somatosensory stimulation (333-μs pulses; amplitude, 2 mA; 64 Hz) was delivered bilaterally via pairs of contact electrodes. A block design paradigm was used in which the stimulus duration increased in pseudo-random order from a single pulse up to 256 electrical pulses (4 s). For CBV measurements, 30 mg/kg of ultrasmall superparamagnetic ironoxide particles (USPIO) injected intravenously, were used. Robust BOLD and CBV HRFs were obtained in the primary somatosensory cortex (S1), secondary somatosensory cortex (S2) and caudate at all stimulus conditions. In particular, BOLD and CBV responses to a single 333-μs-long stimulus were reliably measured, and the CBV HRF presented shorter onset time and time to peak than the BOLD HRF. Both the size of the regions of activation and the peak amplitude of the HRFs grew quickly with increasing stimulus duration, and saturated for stimulus durations greater than 1 s. Onset times in S1 and S2 were faster than in caudate. Finally, the fine spatiotemporal features of the HRF in awake marmosets were similar to those obtained in humans, indicating that the continued refinement of awake non-human primate models is essential to maximize the applicability of animal functional MRI studies to the investigation of human brain function.

Quantitative relations between BOLD responses, cortical energetics, and impulse firing

The blood oxygen level-dependent (BOLD) functional magnetic resonance imaging signal arises as a consequence of changes in blood flow and oxygen usage that in turn are modulated by changes in neural activity. Much attention has been given to both theoretical and experimental aspects of the energetics but not to the neural activity. Here we identify the best energetic theory for the steady-state BOLD signal on the basis of correct predictions of experimental observations. This theory is then used, together with the recently determined relationship between energetics and neural activity, to predict how the BOLD signal changes with activity. Unlike existing treatments, this new theory incorporates a nonzero baseline activity in a completely consistent way and is thus able to account for both sustained positive and negative BOLD signals. We also show that the increase in BOLD signal for a given increase in activity is significantly smaller the larger the baseline activity, as is experimentally observed. Furthermore, the decline of the positive BOLD signal arising from deeper cortical laminae in response to an increase in neural firing is shown to arise as a consequence of the larger baseline activity in deeper laminae. Finally, we provide quantitative relations integrating BOLD responses, energetics, and impulse firing, which among other predictions give the same results as existing theories when the baseline activity is zero. NEW & NOTEWORTHY We use a recently established relation between energetics and neural activity to give a quantitative account of BOLD dependence on neural activity. The incorporation of a nonzero baseline neural activity accounts for positive and negative BOLD signals, shows that changes in neural activity give BOLD changes that are smaller the larger the baseline, and provides a basis for the observed inverse relation between BOLD responses and the depth of cortical laminae giving rise to them.

Investigating the spatiotemporal characteristics of the deoxyhemoglobin-related and deoxyhemoglobin-unrelated functional hemodynamic response across cortical layers in awake marmosets

Blood oxygenation level dependent (BOLD) functional magnetic resonance imaging (fMRI) has become a major tool to map neural activity. However, the spatiotemporal characteristics of the BOLD functional hemodynamic response across the cortical layers remain poorly understood. While human fMRI studies suffer from low spatiotemporal resolution, the use of anesthesia in animal models introduces confounding factors. Additionally, inflow contributions to the fMRI signal become non-negligible when short repetition times (TRs) are used. In the present work, we mapped the BOLD fMRI response to somatosensory stimulation in awake marmosets. To address the above technical concerns, we used a dual-echo gradient-recalled echo planar imaging (GR-EPI) sequence to separate the deoxyhemoglobin-related response (absolute T2* differences) from the deoxyhemoglobin-unrelated response (relative S0 changes). We employed a spatial saturation pulse to saturate incoming arterial spins and reduce inflow effects. Functional GR-EPI images were obtained from a single coronal slice with two different echo times (13.5 and 40.5ms) and TR=0.2s. BOLD, T2*, and S0 images were calculated and their functional responses were detected in both hemispheres of primary somatosensory cortex, from which five laminar regions (L1+2, L3, L4, L5, and L6) were derived. The spatiotemporal distribution of the BOLD response across the cortical layers was heterogeneous, with the middle layers having the highest BOLD amplitudes and shortest onset times. ΔT2* also showed a similar trend. However, functional S0 changes were detected only in L1+2, with a fast onset time. Because inflow effects were minimized, the source of S0 functional changes in L1+2 could be attributed to a reduction of cerebrospinal fluid volume fraction due to the functional increase in cerebral blood volume and to unmodeled T2* changes in the extra- and intra-venous compartments. Caution should be exercised when interpreting laminar BOLD fMRI changes in superficial layers as surrogates of underlying neural activity.

Evaluation of prolonged administration of isoflurane on cerebral blood flow and default mode network in macaque monkeys anesthetized with different maintenance doses

OBJECT

Isoflurane is a commonly used volatile anesthetic agent in clinical anesthesia and biomedical research. Prior study suggested the cerebral blood flow (CBF) and default mode network (DMN) could be changed after prolonged administration of isoflurane. The normal maintenance doses of isoflurane may vary from light (∼0.75%) to deep (∼1.5 or 2%) anesthesia. However, it is not clear how the duration effects are affected by the altered doses. The present study is aimed to examine if the duration effects are affected when isoflurane concentration is altered within normal maintenance doses.

MATERIALS AND METHODS

Adult rhesus monkeys (n=5, 8-12 years old, 8-10kg) were anesthetized and maintained at isoflurane levels 0.89±0.03%, 1.05±0.12%, or 1.19±0.08%. CBF and DMN of monkeys were examined using arterial spin-labeling perfusion and resting state functional MRI techniques.

RESULTS

the functional connectivity (FC) in the dominant DMN (posterior cingulate cortex (PCC) to anterior cingulated cortex (ACC) or media prefrontal cortex (MPFC)) decreased substantially and similarly during 4-h administration of isoflurane at any given maintenance dosage. CBF changes varied with isoflurane dosage. At the low dose (∼0.89%), CBF decreased in most brain regions. In contrast, no obvious changes was seen in those regions (except for the subcortex) when higher doses of isoflurane were applied.

CONCLUSION

FC in DMN was reduced substantially during prolonged administration of isoflurane. The FC reduction was not varying significantly with maintenance doses of isoflurane but the duration effect on CBF was dose-dependent. Such duration effects of isoflurane administration on DMN and CBF should be considered in the interpretation of the outcome in related neuroimaging studies of anesthetized subjects.

Determining Excitatory and Inhibitory Neuronal Activity from Multimodal fMRI Data Using a Generative Hemodynamic Model

Hemodynamic responses, in general, and the blood oxygenation level-dependent (BOLD) fMRI signal, in particular, provide an indirect measure of neuronal activity. There is strong evidence that the BOLD response correlates well with post-synaptic changes, induced by changes in the excitatory and inhibitory (E-I) balance between active neuronal populations. Typical BOLD responses exhibit transients, such as the early-overshoot and post-stimulus undershoot, that can be linked to transients in neuronal activity, but they can also result from vascular uncoupling between cerebral blood flow (CBF) and venous cerebral blood volume (venous CBV). Recently, we have proposed a novel generative hemodynamic model of the BOLD signal within the dynamic causal modeling framework, inspired by physiological observations, called P-DCM (Havlicek et al., 2015). We demonstrated the generative model's ability to more accurately model commonly observed neuronal and vascular transients in single regions but also effective connectivity between multiple brain areas (Havlicek et al., 2017b). In this paper, we additionally demonstrate the versatility of the generative model to jointly explain dynamic relationships between neuronal and hemodynamic physiological variables underlying the BOLD signal using multi-modal data. For this purpose, we utilized three distinct data-sets of experimentally induced responses in the primary visual areas measured in human, cat, and monkey brain, respectively: (1) CBF and BOLD responses; (2) CBF, total CBV, and BOLD responses (Jin and Kim, 2008); and (3) positive and negative neuronal and BOLD responses (Shmuel et al., 2006). By fitting the generative model to the three multi-modal experimental data-sets, we showed that the presence or absence of dynamic features in the BOLD signal is not an unambiguous indication of presence or absence of those features on the neuronal level. Nevertheless, the generative model that takes into account the dynamics of the physiological mechanisms underlying the BOLD response allowed dissociating neuronal from vascular transients and deducing excitatory and inhibitory neuronal activity time-courses from BOLD data alone and from multi-modal data.

Weak correlations between hemodynamic signals and ongoing neural activity during the resting state

Spontaneous fluctuations in hemodynamic signals in the absence of a task or overt stimulation are used to infer neural activity. We tested this coupling by simultaneously measuring neural activity and changes in cerebral blood volume (CBV) in the somatosensory cortex of awake, head-fixed mice during periods of true rest, and during whisker stimulation and volitional whisking. Here we show that neurovascular coupling was similar across states, and large spontaneous CBV changes in the absence of sensory input were driven by volitional whisker and body movements. Hemodynamic signals during periods of rest were weakly correlated with neural activity. Spontaneous fluctuations in CBV and vessel diameter persisted when local neural spiking and glutamatergic input was blocked, and during blockade of noradrenergic receptors, suggesting a non-neuronal origin for spontaneous CBV fluctuations. Spontaneous hemodynamic signals reflect a combination of behavior, local neural activity, and putatively non-neural processes.

Laminar fMRI: What can the time domain tell us?

The rapid developments in functional MRI (fMRI) acquisition methods and hardware technologies in recent years, particularly at high field (≥7 T), have enabled unparalleled visualization of functional detail at a laminar or columnar level, bringing fMRI close to the intrinsic resolution of brain function. These advances highlight the potential of high resolution fMRI to be a valuable tool to study the fundamental processing performed in cortical micro-circuits, and their interactions such as feedforward and feedback processes. Notably, because fMRI measures neuronal activity via hemodynamics, the ultimate resolution it affords depends on the spatial specificity of hemodynamics to neuronal activity at a detailed spatial scale, and by the evolution of this specificity over time. Several laminar (≤1 mm spatial resolution) fMRI studies have examined spatial characteristics of the measured hemodynamic signals across cortical depth, in light of understanding or improving the spatial specificity of laminar fMRI. Few studies have examined temporal features of the hemodynamic response across cortical depth. Temporal features of the hemodynamic response offer an additional means to improve the specificity of fMRI, and could help target neuronal processes and neurovascular coupling relationships across laminae, for example by differences in the onset times of the response across cortical depth. In this review, we discuss factors that affect the timing of neuronal and hemodynamic responses across laminae, touching on the neuronal laminar organization, and focusing on the laminar vascular organization. We provide an overview of hemodynamics across the cortical vascular tree based on optical imaging studies, and review temporal aspects of hemodynamics that have been examined across cortical depth in high spatiotemporal resolution fMRI studies. Last, we discuss the limits and potential of high spatiotemporal resolution fMRI to study laminar neurovascular coupling and neuronal processes.

The hemodynamic signal as a first-order low-pass temporal filter: Evidence and implications for neuroimaging studies

Neuronal activation triggers local changes in blood flow and hemoglobin oxygenation. These hemodynamic signals can be recorded through functional magnetic resonance imaging or intrinsic optical imaging, and allows inferring neural activity in response to stimuli. These techniques are widely used to uncover functional brain architectures. However, their accuracy suffers from distortions inherent to hemodynamic responses and noise. The analysis of these signals currently relies on models of impulse hemodynamic responses to brief stimuli. Here, in order to infer precise functional architectures, we focused on integrated signals associated to the dynamic response of functional maps. To this end, we recorded orientation and direction maps in cat primary visual cortex and compared two protocols: the conventional episodic stimulation technique and a continuous, periodic stimulation paradigm. Conventional methods show that the dynamics of activation and deactivation of the functional maps follows a linear first-order differential equation representing a low-pass filter. Comparison with the periodic stimulation methods confirmed this observation: the phase shifts and magnitude attenuations extracted at various frequencies were consistent with a low-pass filter with a 5s time constant. This dynamics presumably reflects the variations in deoxyhemoglobin mediated by arterial dilations. This dynamics open new avenues in the analysis of neuroimaging data that differs from common methods based on the hemodynamic response function. In particular, we demonstrate that inverting this first-order low-pass filter minimized the distortions of the signal and enabled a much faster and accurate reconstruction of functional maps.

Vascular density and distribution in neocortex

An amazingly wide range of complex behavior emerges from the cerebral cortex. Much of the information processing that leads to these behaviors is performed in neocortical circuits that span throughout the six layers of the cortex. Maintaining this circuit activity requires substantial quantities of oxygen and energy substrates, which are delivered by the complex yet well-organized and tightly-regulated vascular system. In this review, we provide a detailed characterization of the most relevant anatomical and functional features of the cortical vasculature. This includes a compilation of the available data on laminar variation of vascular density and the topological aspects of the microvascular system. We also review the spatio-temporal dynamics of cortical blood flow regulation and oxygenation, many aspects of which remain poorly understood. Finally, we discuss some of the important implications of vascular density, distribution, oxygenation and blood flow regulation for (laminar) fMRI.

Foundations of layer-specific fMRI and investigations of neurophysiological activity in the laminarized neocortex and olfactory bulb of animal models

Laminar organization of neuronal circuits is a recurring feature of how the brain processes information. For instance, different layers compartmentalize different cell types, synaptic activities, and have unique intrinsic and extrinsic connections that serve as units for specialized signal processing. Functional MRI is an invaluable tool to investigate laminar processing in the in vivo human brain, but it measures neuronal activity indirectly by way of the hemodynamic response. Therefore, the accuracy of high-resolution laminar fMRI depends on how precisely it can measure localized microvascular changes nearest to the site of evoked activity. To determine the specificity of fMRI responses to the true neurophysiological responses across layers, the flexibility to invasive procedures in animal models has been necessary. In this review, we will examine different fMRI contrasts and their appropriate uses for layer-specific fMRI, and how localized laminar processing was examined in the neocortex and olfactory bulb. Through collective efforts, it was determined that microvessels, including capillaries, are regulated within single layers and that several endogenous and contrast-enhanced fMRI contrast mechanisms can separate these neural-specific vascular changes from the nonspecific, especially cerebral blood volume-weighted fMRI with intravenous contrast agent injection. We will also propose some open questions that are relevant for the successful implementation of layer-specific fMRI and its potential future directions to study laminar processing when combined with optogenetics.

Time to wake up: Studying neurovascular coupling and brain-wide circuit function in the un-anesthetized animal

Functional magnetic resonance imaging (fMRI) has allowed the noninvasive study of task-based and resting-state brain dynamics in humans by inferring neural activity from blood-oxygenation-level dependent (BOLD) signal changes. An accurate interpretation of the hemodynamic changes that underlie fMRI signals depends on the understanding of the quantitative relationship between changes in neural activity and changes in cerebral blood flow, oxygenation and volume. While there has been extensive study of neurovascular coupling in anesthetized animal models, anesthesia causes large disruptions of brain metabolism, neural responsiveness and cardiovascular function. Here, we review work showing that neurovascular coupling and brain circuit function in the awake animal are profoundly different from those in the anesthetized state. We argue that the time is right to study neurovascular coupling and brain circuit function in the awake animal to bridge the physiological mechanisms that underlie animal and human neuroimaging signals, and to interpret them in light of underlying neural mechanisms. Lastly, we discuss recent experimental innovations that have enabled the study of neurovascular coupling and brain-wide circuit function in un-anesthetized and behaving animal models.

Assessing sensory versus optogenetic network activation by combining (o)fMRI with optical Ca2+ recordings

Encoding of sensory inputs in the cortex is characterized by sparse neuronal network activation. Optogenetic stimulation has previously been combined with fMRI (ofMRI) to probe functional networks. However, for a quantitative optogenetic probing of sensory-driven sparse network activation, the level of similarity between sensory and optogenetic network activation needs to be explored. Here, we complement ofMRI with optic fiber-based population Ca²⁺ recordings for a region-specific readout of neuronal spiking activity in rat brain. Comparing Ca²⁺ responses to the blood oxygenation level-dependent signal upon sensory stimulation with increasing frequencies showed adaptation of Ca²⁺ transients contrasted by an increase of blood oxygenation level-dependent responses, indicating that the optical recordings convey complementary information on neuronal network activity to the corresponding hemodynamic response. To study the similarity of optogenetic and sensory activation, we quantified the density of cells expressing channelrhodopsin-2 and modeled light propagation in the tissue. We estimated the effectively illuminated volume and numbers of optogenetically stimulated neurons, being indicative of sparse activation. At the functional level, upon either sensory or optogenetic stimulation we detected single-peak short-latency primary Ca²⁺ responses with similar amplitudes and found that blood oxygenation level-dependent responses showed similar time courses. These data suggest that ofMRI can serve as a representative model for functional brain mapping.

The Comparisons of Cerebral Hemodynamics Induced by Obstructive Sleep Apnea with Arousal and Periodic Limb Movement with Arousal: A Pilot NIRS Study

Obstructive sleep apnea syndrome (OSA) and restless legs syndrome (RLS) with periodic limb movement during sleep (PLMS) are two sleep disorders characterized by repetitive respiratory or movement events associated with cortical arousals. We compared the cerebral hemodynamic changes linked to periodic apneas/hypopneas with arousals (AHA) in four OSA-patients with periodic limb movements (PLMA) with arousals in four patients with RLS-PLMS using near-infrared spectroscopy (NIRS). AHA induced homogenous pattern of periodic fluctuations in oxygenated (HbO2) and deoxygenated (HHb) hemoglobin, i.e., the decrease of HbO2 was accompanied by an increase of HHb during the respiratory event and resolved to reverse pattern when cortical arousal started. Blood volume (BV) showed the same pattern as HHb but with relative smaller amplitude in most of the AHA events.These changing patterns were significant as Wilcoxon signed-rank tests gave p < 0.001 when comparing the area under the curve of these hemodynamic parameters to zero. By contrast, in PLMA limb movements induced periodic increments in HbO2 and BV (Wilcoxon signed-rank tests, p < 0.001), but HHb changed more heterogeneously even during the events coming from the same patient. Heart rate (HR) also showed different patterns between AHA and PLMA. It significantly decreased during the respiratory event (Wilcoxon signed-rank test, p < 0.001) and then increased after the occurrence of cortical arousal (Wilcoxon signed-rank test, p < 0.001); while in PLMA HR first increased preceding the occurrence of cortical arousal (Wilcoxon signed-rank test, p < 0.001) and then decreased. The results of this preliminary study show that both AHA and PLMA induce changes in cerebral hemodynamics. The occurrence of cortical arousal is accompanied by increased HR in both events, but by different BV changes (i.e., decreased/increased BV in AHA/PLMA, respectively). HR changes may partially account for the increased cerebral hemodynamics during PLMA; whereas in AHA probable vasodilatation mediated by hypoxia/hypercapnia is more crucial for the post-arousal hemodynamics. The differences between changes of cerebral hemodynamics and HR may indicate different pathological mechanisms behind these two sleep disorder events.

Effects of Voluntary Locomotion and Calcitonin Gene-Related Peptide on the Dynamics of Single Dural Vessels in Awake Mice

The dura mater is a vascularized membrane surrounding the brain and is heavily innervated by sensory nerves. Our knowledge of the dural vasculature has been limited to pathological conditions, such as headaches, but little is known about the dural blood flow regulation during behavior. To better understand the dynamics of dural vessels during behavior, we used two-photon laser scanning microscopy (2PLSM) to measure the diameter changes of single dural and pial vessels in the awake mouse during voluntary locomotion. Surprisingly, we found that voluntary locomotion drove the constriction of dural vessels, and the dynamics of these constrictions could be captured with a linear convolution model. Dural vessel constrictions did not mirror the large increases in intracranial pressure (ICP) during locomotion, indicating that dural vessel constriction was not caused passively by compression. To study how behaviorally driven dynamics of dural vessels might be altered in pathological states, we injected the vasodilator calcitonin gene-related peptide (CGRP), which induces headache in humans. CGRP dilated dural, but not pial, vessels and significantly reduced spontaneous locomotion but did not block locomotion-induced constrictions in dural vessels. Sumatriptan, a drug commonly used to treat headaches, blocked the vascular and behavioral the effects of CGRP. These findings suggest that, in the awake animal, the diameters of dural vessels are regulated dynamically during behavior and during drug-induced pathological states.

A Spatiotemporal Profile of In Vivo Cerebral Blood Flow Changes Following Intranasal Oxytocin in Humans

BACKGROUND

Animal and human studies highlight the role of oxytocin in social cognition and behavior and the potential of intranasal oxytocin (IN-OT) to treat social impairment in individuals with neuropsychiatric disorders such as autism. However, extensive efforts to evaluate the central actions and therapeutic efficacy of IN-OT may be marred by the absence of data regarding its temporal dynamics and sites of action in the living human brain.

METHODS

In a placebo-controlled study, we used arterial spin labeling to measure IN-OT-induced changes in resting regional cerebral blood flow (rCBF) in 32 healthy men. Volunteers were blinded regarding the nature of the compound they received. The rCBF data were acquired 15 min before and up to 78 min after onset of treatment onset (40 IU of IN-OT or placebo). The data were analyzed using mass univariate and multivariate pattern recognition techniques.

RESULTS

We obtained robust evidence delineating an oxytocinergic network comprising regions expected to express oxytocin receptors, based on histologic evidence, and including core regions of the brain circuitry underpinning social cognition and emotion processing. Pattern recognition on rCBF maps indicated that IN-OT-induced changes were sustained over the entire posttreatment observation interval (25-78 min) and consistent with a pharmacodynamic profile showing a peak response at 39-51 min.

CONCLUSIONS

Our study provides the first visualization and quantification of IN-OT-induced changes in rCBF in the living human brain unaffected by cognitive, affective, or social manipulations. Our findings can inform theoretical and mechanistic models regarding IN-OT effects on typical and atypical social behavior and guide future experiments (e.g., regarding the timing of experimental manipulations).

Optogenetic fMRI in the mouse hippocampus: Hemodynamic response to brief glutamatergic stimuli

The combination of optogenetics with functional magnetic resonance imaging is a promising tool to study the causal relationship between specific neuronal populations and global brain activity. We employed this technique to study the brain response to recruitment of glutamatergic neurons in the mouse hippocampus. The light-sensitive protein channelrhodopsin-2 was expressed in α-CamKII-positive glutamatergic neurons in the left hippocampus (N = 10). Functional magnetic resonance imaging was performed during local laser stimulation, with stimulus duration of 1 second. The hemodynamic response to these stimuli was analyzed on a whole-brain level. In a secondary analysis, we examined the impact of the stimulation locus on the dorso-ventral axis within the hippocampal formation. The hemodynamic response in the mouse hippocampus had an earlier peak and a shorter duration compared to those observed in humans. Photostimulation was associated with significantly increased blood oxygen level-dependent signal in group statistics: bilaterally in the hippocampus, frontal lobe and septum, ipsilaterally in the nucleus accumbens and contralaterally in the striatum. More dorsal position of the laser fiber was associated with a stronger activation in projection regions (insular cortex and striatum). The characterization of brain-region-specific hemodynamic response functions may enable more precise interpretation of future functional magnetic resonance imaging experiments.

Comparison of glomerular activity patterns by fMRI and wide-field calcium imaging: Implications for principles underlying odor mapping

Functional imaging signals arise from distinct metabolic and hemodynamic events at the neuropil, but how these processes are influenced by pre- and post-synaptic activities need to be understood for quantitative interpretation of stimulus-evoked mapping data. The olfactory bulb (OB) glomeruli, spherical neuropil regions with well-defined neuronal circuitry, can provide insights into this issue. Optical calcium-sensitive fluorescent dye imaging (OICa(2+)) reflects dynamics of pre-synaptic input to glomeruli, whereas high-resolution functional magnetic resonance imaging (fMRI) using deoxyhemoglobin contrast reveals neuropil function within the glomerular layer where both pre- and post-synaptic activities contribute. We imaged odor-specific activity patterns of the dorsal OB in the same anesthetized rats with fMRI and OICa(2+) and then co-registered the respective maps to compare patterns in the same space. Maps by each modality were very reproducible as trial-to-trial patterns for a given odor, overlapping by ~80%. Maps evoked by ethyl butyrate and methyl valerate for a given modality overlapped by ~80%, suggesting activation of similar dorsal glomerular networks by these odors. Comparison of maps generated by both methods for a given odor showed ~70% overlap, indicating similar odor-specific maps by each method. These results suggest that odor-specific glomerular patterns by high-resolution fMRI primarily tracks pre-synaptic input to the OB. Thus combining OICa(2+) and fMRI lays the framework for studies of OB processing over a range of spatiotemporal scales, where OICa(2+) can feature the fast dynamics of dorsal glomerular clusters and fMRI can map the entire glomerular sheet in the OB.

Transcranial functional ultrasound imaging of the brain using microbubble-enhanced ultrasensitive Doppler

Functional ultrasound (fUS) is a novel neuroimaging technique, based on high-sensitivity ultrafast Doppler imaging of cerebral blood volume, capable of measuring brain activation and connectivity in rodents with high spatiotemporal resolution (100 μm, 1 ms). However, the skull attenuates acoustic waves, so fUS in rats currently requires craniotomy or a thinned-skull window. Here we propose a non-invasive approach by enhancing the fUS signal with a contrast agent, inert gas microbubbles. Plane-wave illumination of the brain at high frame rate (500 Hz compounded sequence with three tilted plane waves, PRF = 1500Hz with a 128 element 15 MHz linear transducer), yields highly-resolved neurovascular maps. We compared fUS imaging performance through the intact skull bone (transcranial fUS) versus a thinned-skull window in the same animal. First, we show that the vascular network of the adult rat brain can be imaged transcranially only after a bolus intravenous injection of microbubbles, which leads to a 9 dB gain in the contrast-to-tissue ratio. Next, we demonstrate that functional increase in the blood volume of the primary sensory cortex after targeted electrical-evoked stimulations of the sciatic nerve is observable transcranially in presence of contrast agents, with high reproducibility (Pearson's coefficient ρ = 0.7 ± 0.1, p = 0.85). Our work demonstrates that the combination of ultrafast Doppler imaging and injection of contrast agent allows non-invasive functional brain imaging through the intact skull bone in rats. These results should ease non-invasive longitudinal studies in rodents and open a promising perspective for the adoption of highly resolved fUS approaches for the adult human brain.

•

We combined ultrafast sensitive Doppler with contrast-enhanced ultrasound imaging.

•

We retrieved highly-resolved neurovascular transcranial maps with contrast agents.

•

The presence of microbubbles compensates for the attenuation from the skull.

•

fUS is sensitive to the local hyperemia in the rat brain through the skull with microbubbles.

•

Transcranial fUS imaging allows non-invasive functional brain studies in rodents.

Mechanical restriction of intracortical vessel dilation by brain tissue sculpts the hemodynamic response

Understanding the spatial dynamics of dilation in the cerebral vasculature is essential for deciphering the vascular basis of hemodynamic signals in the brain. We used two-photon microscopy to image neural activity and vascular dynamics in the somatosensory cortex of awake behaving mice during voluntary locomotion. Arterial dilations within the histologically-defined forelimb/hindlimb (FL/HL) representation were larger than arterial dilations in the somatosensory cortex immediately outside the FL/HL representation, demonstrating that the vascular response during natural behaviors was spatially localized. Surprisingly, we found that locomotion drove dilations in surface vessels that were nearly three times the amplitude of intracortical vessel dilations. The smaller dilations of the intracortical arterioles were not due to saturation of dilation. Anatomical imaging revealed that, unlike surface vessels, intracortical vessels were tightly enclosed by brain tissue. A mathematical model showed that mechanical restriction by the brain tissue surrounding intracortical vessels could account for the reduced amplitude of intracortical vessel dilation relative to surface vessels. Thus, under normal conditions, the mechanical properties of the brain may play an important role in sculpting the laminar differences of hemodynamic responses.

Long-latency reductions in gamma power predict hemodynamic changes that underlie the negative BOLD signal

Studies that use prolonged periods of sensory stimulation report associations between regional reductions in neural activity and negative blood oxygenation level-dependent (BOLD) signaling. However, the neural generators of the negative BOLD response remain to be characterized. Here, we use single-impulse electrical stimulation of the whisker pad in the anesthetized rat to identify components of the neural response that are related to "negative" hemodynamic changes in the brain. Laminar multiunit activity and local field potential recordings of neural activity were performed concurrently with two-dimensional optical imaging spectroscopy measuring hemodynamic changes. Repeated measurements over multiple stimulation trials revealed significant variations in neural responses across session and animal datasets. Within this variation, we found robust long-latency decreases (300 and 2000 ms after stimulus presentation) in gamma-band power (30-80 Hz) in the middle-superficial cortical layers in regions surrounding the activated whisker barrel cortex. This reduction in gamma frequency activity was associated with corresponding decreases in the hemodynamic responses that drive the negative BOLD signal. These findings suggest a close relationship between BOLD responses and neural events that operate over time scales that outlast the initiating sensory stimulus, and provide important insights into the neurophysiological basis of negative neuroimaging signals.