Cerebral vasomotion: a 0.1-Hz oscillation in reflected light imaging of neural activity

Alterations in neurovascular coupling following acute traumatic brain injury

Following acute traumatic brain injury (TBI), timely transport to a hospital can significantly improve the prognosis for recovery. There is, however, a dearth of quantitative biomarkers for brain injury that can be rapidly acquired and interpreted in active, field environments in which TBIs are frequently incurred. We explored potential functional indicators for TBI that can be noninvasively obtained through portable detection modalities, namely optical and electrophysiological approaches. By combining diffuse correlation spectroscopy with colocalized electrophysiological measurements in a mouse model of TBI, we observed concomitant alterations in sensory-evoked cerebral blood flow (CBF) and electrical potentials following controlled cortical impact. Injury acutely reduced the peak amplitude of both electrophysiological and CBF responses, which mostly recovered to baseline values within 30 min, and intertrial variability for these parameters was also acutely altered. Notably, the postinjury dynamics of the CBF overshoot and undershoot amplitudes differed significantly; whereas the amplitude of the initial peak of stimulus-evoked CBF recovered relatively rapidly, the ensuing undershoot did not appear to recover within 30 min of injury. Additionally, acute injury induced apparent low-frequency oscillatory behavior in CBF ([Formula: see text]). Histological assessment indicated that these physiological alterations were not associated with any major, persisting anatomical changes. Several time-domain features of the blood flow and electrophysiological responses showed strong correlations in recovery kinetics. Overall, our results reveal an array of stereotyped, injury-induced alterations in electrophysiological and hemodynamic responses that can be rapidly obtained using a combination of portable detection techniques.

Effect of confounding variables on hemodynamic response function estimation using averaging and deconvolution analysis: An event-related NIRS study

Slow and rapid event-related designs are used in fMRI and functional near-infrared spectroscopy (fNIRS) experiments to temporally characterize the brain hemodynamic response to discrete events. Conventional averaging (CA) and the deconvolution method (DM) are the two techniques commonly used to estimate the Hemodynamic Response Function (HRF) profile in event-related designs. In this study, we conducted a series of simulations using synthetic and real NIRS data to examine the effect of the main confounding factors, including event sequence timing parameters, different types of noise, signal-to-noise ratio (SNR), temporal autocorrelation and temporal filtering on the performance of these techniques in slow and rapid event-related designs. We also compared systematic errors in the estimates of the fitted HRF amplitude, latency and duration for both techniques. We further compared the performance of deconvolution methods based on Finite Impulse Response (FIR) basis functions and gamma basis sets. Our results demonstrate that DM was much less sensitive to confounding factors than CA. Event timing was the main parameter largely affecting the accuracy of CA. In slow event-related designs, deconvolution methods provided similar results to those obtained by CA. In rapid event-related designs, our results showed that DM outperformed CA for all SNR, especially above -5 dB regardless of the event sequence timing and the dynamics of background NIRS activity. Our results also show that periodic low-frequency systemic hemodynamic fluctuations as well as phase-locked noise can markedly obscure hemodynamic evoked responses. Temporal autocorrelation also affected the performance of both techniques by inducing distortions in the time profile of the estimated hemodynamic response with inflated t-statistics, especially at low SNRs. We also found that high-pass temporal filtering could substantially affect the performance of both techniques by removing the low-frequency components of HRF profiles. Our results emphasize the importance of characterization of event timing, background noise and SNR when estimating HRF profiles using CA and DM in event-related designs.

Cerebral functional connectivity and Mayer waves in mice: Phenomena and separability

Resting-state functional connectivity is a growing neuroimaging approach that analyses the spatiotemporal structure of spontaneous brain activity, often using low-frequency (<0.08 Hz) hemodynamics. In addition to these fluctuations, there are two other low-frequency hemodynamic oscillations in a nearby spectral region (0.1-0.4 Hz) that have been reported in the brain: vasomotion and Mayer waves. Despite how close in frequency these phenomena exist, there is little research on how vasomotion and Mayer waves are related to or affect resting-state functional connectivity. In this study, we analyze spontaneous hemodynamic fluctuations over the mouse cortex using optical intrinsic signal imaging. We found spontaneous occurrence of oscillatory hemodynamics ∼0.2 Hz consistent with the properties of Mayer waves reported in the literature. Across a group of mice (n = 19), there was a large variability in the magnitude of Mayer waves. However, regardless of the magnitude of Mayer waves, functional connectivity patterns could be recovered from hemodynamic signals when filtered to the lower frequency band, 0.01-0.08 Hz. Our results demonstrate that both Mayer waves and resting-state functional connectivity patterns can co-exist simultaneously, and that they can be separated by applying bandpass filters.

Impacting the effect of fMRI noise through hardware and acquisition choices - Implications for controlling false positive rates

We review the components of time-series noise in fMRI experiments and the effect of image acquisition parameters on the noise. In addition to helping determine the total amount of signal and noise (and thus temporal SNR), the acquisition parameters have been shown to be critical in determining the ratio of thermal to physiological induced noise components in the time series. Although limited attention has been given to this latter metric, we show that it determines the degree of spatial correlations seen in the time-series noise. The spatially correlations of the physiological noise component are well known, but recent studies have shown that they can lead to a higher than expected false-positive rate in cluster-wise inference based on parametric statistical methods used by many researchers. Based on understanding the effect of acquisition parameters on the noise mixture, we propose several acquisition strategies that might be helpful reducing this elevated false-positive rate, such as moving to high spatial resolution or using highly-accelerated acquisitions where thermal sources dominate. We suggest that the spatial noise correlations at the root of the inflated false-positive rate problem can be limited with these strategies, and the well-behaved spatial auto-correlation functions (ACFs) assumed by the conventional statistical methods are retained if the high resolution data is smoothed to conventional resolutions.

Noise and non-neuronal contributions to the BOLD signal: applications to and insights from animal studies

The BOLD signal reflects hemodynamic events within the brain, which in turn are driven by metabolic changes and neural activity. However, the link between BOLD changes and neural activity is indirect and can be influenced by a number of non-neuronal processes. Motion and physiological cycles have long been known to affect the BOLD signal and are present in both humans and animal models. Differences in physiological baseline can also contribute to intra- and inter-subject variability. The use of anesthesia, common in animal studies, alters neural activity, vascular tone, and neurovascular coupling. Most intriguing, perhaps, are the contributions from other processes that do not appear to be neural in origin but which may provide information about other aspects of neurophysiology. This review discusses different types of noise and non-neuronal contributors to the BOLD signal, sources of variability for animal studies, and insights to be gained from animal models.

Brain-heart communication: Evidence for "central pacemaker" oscillations with a dominant frequency at 0.1Hz in the cingulum

OBJECTIVES

In the brain and heart, oscillations at about 0.1Hz are conspicuous. It is therefore worthwhile to study the interaction between intrinsic BOLD oscillations (0.1Hz) and slow oscillations in heart rate interval (RRI) signals and differentiate between their neural and vascular origin.

METHODS

We studied the phase-coupling with a 3T scanner with high scanning rate between BOLD signals in 22 regions and simultaneously recorded RRI oscillations in 23 individuals in two resting states.

RESULTS

By applying a hierarchical cluster analysis, it was possible to separate two clusters of phase-coupling between slow BOLD and RRI oscillations in the midcingulum, one representative for neural and the other for vascular BOLD oscillations. About half of the participants revealed positive time delays characteristic for neural BOLD oscillations and neurally-driven RRI oscillations.

CONCLUSIONS

The results suggest that slow vascular and neural BOLD oscillations can be differentiated and that intrinsic oscillations (0.1Hz) originate in the cingulum or its close vicinity and contribute to heart rate variability (HRV).

SIGNIFICANCE

The study provides new insights into the dynamics of resting state activities, helps to explain HRV, and offers the possibility to investigate slow rhythmic neural activity changes in different brain regions without EEG recording.

Functional Connectivity of Resting Hemodynamic Signals in Submillimeter Orientation Columns of the Visual Cortex

Resting-state functional magnetic resonance imaging has been increasingly used for examining connectivity across brain regions. The spatial scale by which hemodynamic imaging can resolve functional connections at rest remains unknown. To examine this issue, deoxyhemoglobin-weighted intrinsic optical imaging data were acquired from the visual cortex of lightly anesthetized ferrets. The neural activity of orientation domains, which span a distance of 0.7-0.8 mm, has been shown to be correlated during evoked activity and at rest. We performed separate analyses to assess the degree to which the spatial and temporal characteristics of spontaneous hemodynamic signals depend on the known functional organization of orientation columns. As a control, artificial orientation column maps were generated. Spatially, resting hemodynamic patterns showed a higher spatial resemblance to iso-orientation maps than artificially generated maps. Temporally, a correlation analysis was used to establish whether iso-orientation domains are more correlated than orthogonal orientation domains. After accounting for a significant decrease in correlation as a function of distance, a small but significant temporal correlation between iso-orientation domains was found, which decreased with increasing difference in orientation preference. This dependence was abolished when using artificially synthetized orientation maps. Finally, the temporal correlation coefficient as a function of orientation difference at rest showed a correspondence with that calculated during visual stimulation suggesting that the strength of resting connectivity is related to the strength of the visual stimulation response. Our results suggest that temporal coherence of hemodynamic signals measured by optical imaging of intrinsic signals exists at a submillimeter columnar scale in resting state.

Is Vasomotion in Cerebral Arteries Impaired in Alzheimer's Disease?

A substantial body of evidence supports the hypothesis of a vascular component in the pathogenesis of Alzheimer’s disease (AD). Cerebral hypoperfusion and blood-brain barrier dysfunction have been indicated as key elements of this pathway. Cerebral amyloid angiopathy (CAA) is a cerebrovascular disorder, frequent in AD, characterized by the accumulation of amyloid-β (Aβ) peptide in cerebral blood vessel walls. CAA is associated with loss of vascular integrity, resulting in impaired regulation of cerebral circulation, and increased susceptibility to cerebral ischemia, microhemorrhages, and white matter damage. Vasomotion— the spontaneous rhythmic modulation of arterial diameter, typically observed in arteries/arterioles in various vascular beds including the brain— is thought to participate in tissue perfusion and oxygen delivery regulation. Vasomotion is impaired in adverse conditions such as hypoperfusion and hypoxia. The perivascular and glymphatic pathways of Aβ clearance are thought to be driven by the systolic pulse. Vasomotion produces diameter changes of comparable amplitude, however at lower rates, and could contribute to these mechanisms of Aβ clearance. In spite of potential clinical interest, studies addressing cerebral vasomotion in the context of AD/CAA are limited. This study reviews the current literature on vasomotion, and hypothesizes potential paths implicating impaired cerebral vasomotion in AD/CAA. Aβ and oxidative stress cause vascular tone dysregulation through direct effects on vascular cells, and indirect effects mediated by impaired neurovascular coupling. Vascular tone dysregulation is further aggravated by cholinergic deficit and results in depressed cerebrovascular reactivity and (possibly) impaired vasomotion, aggravating regional hypoperfusion and promoting further Aβ and oxidative stress accumulation.

Intraoperative optical mapping of epileptogenic cortices during non-ictal periods in pediatric patients

Complete removal of epileptogenic cortex while preserving eloquent areas is crucial in patients undergoing epilepsy surgery. In this manuscript, the feasibility was explored of developing a new methodology based on dynamic intrinsic optical signal imaging (DIOSI) to intraoperatively detect and differentiate epileptogenic from eloquent cortices in pediatric patients with focal epilepsy. From 11 pediatric patients undergoing epilepsy surgery, negatively-correlated hemodynamic low-frequency oscillations (LFOs, ~ 0.02-0.1 Hz) were observed from the exposed epileptogenic and eloquent cortical areas, as defined by electrocorticography (ECoG), using a DIOSI system. These LFOs were classified into multiple groups in accordance with their unique temporal profiles. Causal relationships within these groups were investigated using the Granger causality method, and 83% of the ECoG-defined epileptogenic cortical areas were found to have a directed influence on one or more cortical areas showing LFOs within the field of view of the imaging system. To understand the physiological origins of LFOs, blood vessel density was compared between epileptogenic and normal cortical areas and a statistically-significant difference (p < 0.05) was detected. The differences in blood-volume and blood-oxygenation dynamics between eloquent and epileptogenic cortices were also uncovered using a stochastic modeling approach. This, in turn, yielded a means by which to separate epileptogenic from eloquent cortex using hemodynamic LFOs. The proposed methodology detects epileptogenic cortices by exploiting the effective connectivity that exists within cortical regions displaying LFOs and the biophysical features contributed by the altered vessel networks within the epileptogenic cortex. It could be used in conjunction with existing technologies for epileptogenic/eloquent cortex localization and thereby facilitate clinical decision-making.

Spontaneous Low-Frequency Cerebral Hemodynamics Oscillations in Restless Legs Syndrome with Periodic Limb Movements During Sleep: A Near-Infrared Spectroscopy Study

Background and Purpose

Periodic limb movements (PLM) during sleep (PLMS) are associated with cortical and cardiovascular activation. Changes in cerebral hemodynamics caused by cortical activity can be measured using near-infrared spectroscopy (NIRS). We investigated oscillatory components of cerebral hemodynamics during PLM and different sleep stages in restless legs syndrome (RLS) patients with PLMS.

Methods

Four female RLS patients with PLMS, and four age- and sex-matched normal controls were included. PLM and sleep stages were scored using polysomnography, while the spontaneous cerebral hemodynamics was measured by NIRS. The phase and amplitude of the cerebral oxyhemoglobin concentration [HbO] and the deoxyhemoglobin concentration [Hb] low-frequency oscillations (LFOs) were evaluated during each sleep stage [waking, light sleep (LS; stages N1 and N2), slow-wave sleep (stage N3), and rapid eye movement (REM) sleep]. In RLS patients with PLMS, the cerebral hemodynamics during LS was divided into LS with and without PLM.

Results

The cerebral hemodynamics activity varied among the different sleep stages. There were changes in phase differences between [HbO] and [Hb] LFOs during the different sleep stages in the normal controls but not in the RLS patients with PLMS. The [HbO] and [Hb] LFO amplitudes were higher in the patient group than in controls during both LS with PLM and REM sleep.

Conclusions

The present study has demonstrated the presence of cerebral hemodynamics disturbances in RLS patients with PLMS, which may contribute to an increased risk of cerebrovascular events.

Vascular autorescaling of fMRI (VasA fMRI) improves sensitivity of population studies: A pilot study

The blood oxygenation level-dependent (BOLD) signal is widely used for functional magnetic resonance imaging (fMRI) of brain function in health and disease. The statistical power of fMRI group studies is significantly hampered by high inter-subject variance due to differences in baseline vascular physiology. Several methods have been proposed to account for physiological vascularization differences between subjects and hence improve the sensitivity in group studies. However, these methods require the acquisition of additional reference scans (such as a full resting-state fMRI session or ASL-based calibrated BOLD). We present a vascular autorescaling (VasA) method, which does not require any additional reference scans. VasA is based on the observation that slow oscillations (<0.1Hz) in arterial blood CO2 levels occur naturally due to changes in respiration patterns. These oscillations yield fMRI signal changes whose amplitudes reflect the blood oxygenation levels and underlying local vascularization and vascular responsivity. VasA estimates proxies of the amplitude of these CO2-driven oscillations directly from the residuals of task-related fMRI data without the need for reference scans. The estimates are used to scale the amplitude of task-related fMRI responses, to account for vascular differences. The VasA maps compared well to cerebrovascular reactivity (CVR) maps and cerebral blood volume maps based on vascular space occupancy (VASO) measurements in four volunteers, speaking to the physiological vascular basis of VasA. VasA was validated in a wide variety of tasks in 138 volunteers. VasA increased t-scores by up to 30% in specific brain areas such as the visual cortex. The number of activated voxels was increased by up to 200% in brain areas such as the orbital frontal cortex while still controlling the nominal false-positive rate. VasA fMRI outperformed previously proposed rescaling approaches based on resting-state fMRI data and can be readily applied to any task-related fMRI data set, even retrospectively.

Robust and fragile aspects of cortical blood flow in relation to the underlying angioarchitecture

We review the organizational principles of the cortical vasculature and the underlying patterns of blood flow under normal conditions and in response to occlusion of single vessels. The cortex is sourced by a two-dimensional network of pial arterioles that feeds a three-dimensional network of subsurface microvessels in close proximity to neurons and glia. Blood flow within the surface and subsurface networks is largely insensitive to occlusion of a single vessel within either network. However, the penetrating arterioles that connect the pial network to the subsurface network are bottlenecks to flow; occlusion of even a single penetrating arteriole results in the death of a 500 μm diameter cylinder of cortical tissue despite the potential for collateral flow through microvessels. This pattern of flow is consistent with that calculated from a full reconstruction of the angioarchitecture. Conceptually, collateral flow is insufficient to compensate for the occlusion of a penetrating arteriole because penetrating venules act as shunts of blood that flows through collaterals. Future directions that stem from the analysis of the angioarchitecture concern cellular-level issues, in particular the regulation of blood flow within the subsurface microvascular network, and system-level issues, in particular the role of penetrating arteriole occlusions in human cognitive impairment.

Time-frequency dynamics of the sum of intra- and extracerebral hemodynamic functional connectivity during resting-state and respiratory challenges assessed by multimodal functional near-infrared spectroscopy

Monitoring respiratory processes is important for evaluating neuroimaging data, given their influence on time-frequency dynamics of intra- and extracerebral hemodynamics. Here we investigated the time-frequency dynamics of the sum of intra- and extracerebral hemodynamic functional connectivity states during hypo- and hypercapnia by using three different respiratory challenge tasks (i.e., hyperventilation, breath-holding, and rebreathing) compared to resting-state. The sum of intra- and extracerebral hemodynamic responses were assessed using functional near-infrared spectroscopy (fNIRS) within two regions of interest (i.e., the dorsolateral and the medial prefrontal cortex). Time-frequency fNIRS analysis was performed based on wavelet transform coherence to quantify functional connectivity in terms of positive and negative phase-coupling within each region of interest. Physiological measures were assessed in the form of partial end-tidal carbon dioxide, heart rate, arterial tissue oxygen saturation, and respiration rate. We found that the three respiration challenges modulated time-frequency dynamics differently with respect to resting-state: 1) Hyperventilation and breath-holding exhibited inverse patterns of positive and negative phase-coupling. 2) In contrast, rebreathing had no significant effect. 3) Low-frequency oscillations contributed to a greater extent to time-frequency dynamics compared to high-frequency oscillations. The results highlight that there exist distinct differences in time-frequency dynamics of the sum of intra- and extracerebral functional connectivity not only between hypo- (hyperventilation) and hypercapnia but also between different states of hypercapnia (breath-holding versus rebreathing). This suggests that a multimodal assessment of intra-/extracerebral and systemic physiological changes during respiratory challenges compared to resting-state may have potential use in the differentiation between physiological and pathological respiratory behavior accompanied by the psycho-physiological state of a human.

Neuronal and physiological correlation to hemodynamic resting-state fluctuations in health and disease

Low-frequency, spatially coherent fluctuations present in functional magnetic resonance imaging time series have had a tremendous impact on brain connectomics. This work aims to explore the degree with which hemodynamic connectivity is associated with neuronal, metabolic, and vascular connectivity measures. For this purpose, GCaMP and nontransgenic mice were used to image neuronal activity and oxidative metabolism activity, respectively, along with blood-oxygenation- and cerebral blood volume (CBV)-sensitive hemodynamic changes from the same animals. Although network clusters calculated using either GCaMP (neuronal activity) or optical imaging of intrinsic signal (OIS)-BOLD (blood oxygenation) data did not exhibit strong spatial similarity, the strengths of node-to-node connectivity measured with these modalities were strongly correlated with one another. This finding suggests that hemodynamic connectivity as measured by blood oxygenation measurements, such as functional connectivity magnetic resonance imaging, is a valuable surrogate for the underlying neuronal connectivity. In nontransgenic animals, greater connectivity correlation was observed between tissue oxidative metabolism (flavoprotein autofluorescence imaging [FAI]) and blood oxygenation measurements, suggesting that metabolic contributions to hemodynamic signals are likely responsible for its significant correlation with neuronal connectivity. Lastly, a mouse model of Alzheimer's disease was used to explore the source of decreases in connectivity reported in these mice, a finding that is thought to be associated with amyloid load-driven metabolic decline. The intercluster connectivity measured by metabolic-sensitive measurements (FAI and OIS-BOLD) was maintained while vascular-only signals (OIS-CBV) provided negligible correlation. Therefore, metabolism-sensitive measurements as used in this work are better positioned to capture changes in neuronal connectivity, such that decreases in hemodynamic connectivity likely reflect decreases in oxidative metabolic function.

Different dynamic resting state fMRI patterns are linked to different frequencies of neural activity

Resting state functional magnetic resonance imaging (rsfMRI) results have indicated that network mapping can contribute to understanding behavior and disease, but it has been difficult to translate the maps created with rsfMRI to neuroelectrical states in the brain. Recently, dynamic analyses have revealed multiple patterns in the rsfMRI signal that are strongly associated with particular bands of neural activity. To further investigate these findings, simultaneously recorded invasive electrophysiology and rsfMRI from rats were used to examine two types of electrical activity (directly measured low-frequency/infraslow activity and band-limited power of higher frequencies) and two types of dynamic rsfMRI (quasi-periodic patterns or QPP, and sliding window correlation or SWC). The relationship between neural activity and dynamic rsfMRI was tested under three anesthetic states in rats: dexmedetomidine and high and low doses of isoflurane. Under dexmedetomidine, the lightest anesthetic, infraslow electrophysiology correlated with QPP but not SWC, whereas band-limited power in higher frequencies correlated with SWC but not QPP. Results were similar under isoflurane; however, the QPP was also correlated to band-limited power, possibly due to the burst-suppression state induced by the anesthetic agent. The results provide additional support for the hypothesis that the two types of dynamic rsfMRI are linked to different frequencies of neural activity, but isoflurane anesthesia may make this relationship more complicated. Understanding which neural frequency bands appear as particular dynamic patterns in rsfMRI may ultimately help isolate components of the rsfMRI signal that are of interest to disorders such as schizophrenia and attention deficit disorder.

Dependence of BOLD signal fluctuation on arterial blood CO2 and O2: Implication for resting-state functional connectivity

Blood oxygenation level dependent (BOLD) functional MRI signal is known to be modulated by the CO2 level. Typically only end-tidal CO2, rather than the arterial partial pressure of CO2 (paCO2), was measured while the arterial partial pressure of O2 (paO2) level was not controlled due to free breathing, making their contribution not separable. Especially, the influences of paO2 and paCO2 on resting-state functional connectivity are not well studied. In this study, we investigated the relationship between paCO2 and resting as well as stimulus-evoked BOLD signals under hyperoxic and hypercapnic manipulation with tight control of arterial paO2. Rats under isoflurane anesthesia were subjected to six inspired gas conditions: 47% O2 in air (Normal), adding 1%, 2% or 5% CO2, carbogen (95% O2/5% CO2), and 100% O2. Somatosensory BOLD activation was significantly increased under 100% O2, while reduced with increased paCO2 levels. However, while resting BOLD connectivity pattern expanded and bilateral correlation increased under 100% O2, the correlation coefficient between the left and right somatosensory cortex was generally not dependent on paCO2 or paO2. Interestingly, the correlation in 0.04-0.07Hz range significantly increased with CO2 levels. Intracortical electrophysiological recordings showed a similar trend as the BOLD but the neurovascular coupling varied. The results suggest that paO2 and paCO2 together rather than paCO2 alone alter the BOLD signal. The response is not purely vascular in nature but has strong neuronal origins. This should be taken into consideration when designing calibrated BOLD experiment and interpreting functional connectivity data especially in aging, under drug, or neurological disorders.

Greater Activity in the Frontal Cortex on Left Curves: A Vector-Based fNIRS Study of Left and Right Curve Driving

OBJECTIVES

In the brain, the mechanisms of attention to the left and the right are known to be different. It is possible that brain activity when driving also differs with different horizontal road alignments (left or right curves), but little is known about this. We found driver brain activity to be different when driving on left and right curves, in an experiment using a large-scale driving simulator and functional near-infrared spectroscopy (fNIRS).

RESEARCH DESIGN AND METHODS

The participants were fifteen healthy adults. We created a course simulating an expressway, comprising straight line driving and gentle left and right curves, and monitored the participants under driving conditions, in which they drove at a constant speed of 100 km/h, and under non-driving conditions, in which they simply watched the screen (visual task). Changes in hemoglobin concentrations were monitored at 48 channels including the prefrontal cortex, the premotor cortex, the primary motor cortex and the parietal cortex. From orthogonal vectors of changes in deoxyhemoglobin and changes in oxyhemoglobin, we calculated changes in cerebral oxygen exchange, reflecting neural activity, and statistically compared the resulting values from the right and left curve sections.

RESULTS

Under driving conditions, there were no sites where cerebral oxygen exchange increased significantly more during right curves than during left curves (p > 0.05), but cerebral oxygen exchange increased significantly more during left curves (p < 0.05) in the right premotor cortex, the right frontal eye field and the bilateral prefrontal cortex. Under non-driving conditions, increases were significantly greater during left curves (p < 0.05) only in the right frontal eye field.

CONCLUSIONS

Left curve driving was thus found to require more brain activity at multiple sites, suggesting that left curve driving may require more visual attention than right curve driving. The right frontal eye field was activated under both driving and non-driving conditions.

Single-trial estimation of the cerebral metabolic rate of oxygen with imaging photoplethysmography and laser speckle contrast imaging

Cortical cerebral metabolic rate of oxygen (CMRO(2)) could conventionally be measured by combining laser Doppler flowmetry and multispectral reflectance imaging across multiple trials of stimulation, which compromises the real-time capacity. Monitoring transient change of CMRO(2) has been challenging. In this Letter, imaging photoplethysmography (iPPG) and laser speckle contrast imaging were combined into a multi-modal optical imaging system for single-trial estimation of CMRO(2). In a physiologically stable experiment, the iPPG-based method showed a less than 4% variance in comparison with the conventional method over 20 trials, and its temporal stability could be comparable to that by conventional method over 6 trials. While the oxygen supply was decreased deliberately, the new method was able to detect the transient changes of CMRO(2) in real time, which could not be revealed by the conventional method.

High frequency functional brain networks in neonates revealed by rapid acquisition resting state fMRI

Understanding how spatially remote brain regions interact to form functional brain networks, and how these develop during the neonatal period, provides fundamental insights into normal brain development, and how mechanisms of brain disorder and recovery may function in the immature brain. A key imaging tool in characterising functional brain networks is examination of T2*-weighted fMRI signal during rest (resting state fMRI, rs-fMRI). The majority of rs-fMRI studies have concentrated on slow signal fluctuations occurring at <0.1 Hz, even though neuronal rhythms, and haemodynamic responses to these fluctuate more rapidly, and there is emerging evidence for crucial information about functional brain connectivity occurring more rapidly than these limits. The characterisation of higher frequency components has been limited by the sampling frequency achievable with standard T2* echoplanar imaging (EPI) sequences. We describe patterns of neonatal functional brain network connectivity derived using accelerated T2*-weighted EPI MRI. We acquired whole brain rs-fMRI data, at subsecond sampling frequency, from preterm infants at term equivalent age and compared this to rs-fMRI data acquired with standard EPI acquisition protocol. We provide the first evidence that rapid rs-fMRI acquisition in neonates, and adoption of an extended frequency range for analysis, allows identification of a substantial proportion of signal power residing above 0.2 Hz. We thereby describe changes in brain connectivity associated with increasing maturity which are not evident using standard rs-fMRI protocols. Development of optimised neonatal fMRI protocols, including use of high speed acquisition sequences, is crucial for understanding the physiology and pathophysiology of the developing brain.

Stimulus-related neuroimaging in task-engaged subjects is best predicted by concurrent spiking

The implicit goal of functional magnetic resonance imaging is to infer local neural activity. There is considerable debate, however, as to whether imaging correlates most linearly with local spiking or some local field potential (LFP) measurement. Through simultaneous neuroimaging (intrinsic-signal optical imaging) and electrode recordings from alert, task-engaged macaque monkeys, we showed previously that local electrophysiology correlates with only a specific stimulus-related imaging component. Here we show that this stimulus-related component--obtained by subtracting a substantial task-related component--is particularly linear with local spiking over a comprehensive range of response strengths. Matches to concurrent LFP measurements are, to varying degrees, poorer. As a control, we also tried matching the full imaging signal to local electrophysiology without subtracting task-related components. These control matches were consistently worse; they were, however, slightly better for gamma LFP than spiking, potentially resolving discrepancies between our findings and earlier reports favoring LFP.

Neurovascular coupling and decoupling in the cortex during voluntary locomotion

Hemodynamic signals are widely used to infer neural activity in the brain. We tested the hypothesis that hemodynamic signals faithfully report neural activity during voluntary behaviors by measuring cerebral blood volume (CBV) and neural activity in the somatosensory cortex and frontal cortex of head-fixed mice during locomotion. Locomotion induced a large and robust increase in firing rate and gamma-band (40-100 Hz) power in the local field potential in the limb representations in somatosensory cortex, and was accompanied by increases in CBV, demonstrating that hemodynamic signals are coupled with neural activity in this region. However, in the frontal cortex, CBV did not change during locomotion, but firing rate and gamma-band power both increased, indicating a decoupling of neural activity from the hemodynamic signal. These results show that hemodynamic signals are not faithful indicators of the mean neural activity in the frontal cortex during locomotion; thus, the results from fMRI and other hemodynamic imaging methodologies for studying neural processes must be interpreted with caution.

Quantitative separation of arterial and venous cerebral blood volume increases during voluntary locomotion

Voluntary locomotion is accompanied by large increases in cortical activity and localized increases in cerebral blood volume (CBV). We sought to quantitatively determine the spatial and temporal dynamics of voluntary locomotion-evoked cerebral hemodynamic changes. We measured single vessel dilations using two-photon microscopy and cortex-wide changes in CBV-related signal using intrinsic optical signal (IOS) imaging in head-fixed mice freely locomoting on a spherical treadmill. During bouts of locomotion, arteries dilated rapidly, while veins distended slightly and recovered slowly. The dynamics of diameter changes of both vessel types could be captured using a simple linear convolution model. Using these single vessel measurements, we developed a novel analysis approach to separate out spatially and temporally distinct arterial and venous components of the location-specific hemodynamic response functions (HRF) for IOS. The HRF of each pixel of was well fit by a sum of a fast arterial and a slow venous component. The HRFs of pixels in the limb representations of somatosensory cortex had a large arterial contribution, while in the frontal cortex the arterial contribution to the HRF was negligible. The venous contribution was much less localized, and was substantial in the frontal cortex. The spatial pattern and amplitude of these HRFs in response to locomotion in the cortex were robust across imaging sessions. Separating the more localized arterial component from the diffuse venous signals will be useful for dealing with the dynamic signals generated by naturalistic stimuli.

Regional variations in vascular density correlate with resting-state and task-evoked blood oxygen level-dependent signal amplitude

Functional magnetic resonance imaging (fMRI) has become one of the primary tools used for noninvasively measuring brain activity in humans. For the most part, the blood oxygen level-dependent (BOLD) contrast is used, which reflects the changes in hemodynamics associated with active brain tissue. The main advantage of the BOLD signal is that it is relatively easy to measure and thus is often used as a proxy for comparing brain function across population groups (i.e., control vs. patient). However, it is particularly weighted toward veins whose structural architecture is known to vary considerably across the brain. This makes it difficult to interpret whether differences in BOLD between cortical areas reflect true differences in neural activity or vascular structure. We therefore investigated how regional variations of vascular density (VAD) relate to the amplitude of resting-state and task-evoked BOLD signals. To address this issue, we first developed an automated method for segmenting veins in images acquired with susceptibility-weighted imaging, allowing us to visualize the venous vascular tree across the brain. In 19 healthy subjects, we then applied voxel-based morphometry (VBM) to T1-weighted images and computed regional measures of gray matter density (GMD). We found that, independent of spatial scale, regional variations in resting-state and task-evoked fMRI amplitudes were better correlated to VAD compared to GMD. Using a general linear model (GLM), it was observed that the bulk of regional variance in resting-state activity could be modeled by VAD. Cortical areas whose resting-state activity was most suppressed by VAD correction included Cuneus, Precuneus, Culmen, and BA 9, 10, and 47. Taken together, our results suggest that resting-state BOLD signals are significantly related to the underlying structure of the brain vascular system. Calibrating resting BOLD activity by venous structure may result in a more accurate interpretation of differences observed between cortical areas and/or individuals.

Entrainment of spontaneous cerebral hemodynamic oscillations to behavioral responses

Entrainment in physiological systems can be manifest in cases where phase-coupling (synchronization) between slow intrinsic oscillations and periodic motor responses, or vice versa, takes place. To test whether voluntary movement has something in common with entrainment of slow hemodynamic oscillations to motor responses, we studied blood pressure (BP), heart rate beat-to-beat intervals (RRI) and prefrontal (de)oxyhemoglobin (Hb/HbO2) during 5min of rest, 10min of self-paced, voluntary movements and 10min of stimulus-paced movements at 10s intervals in 9 subjects. Subjects were divided into 2 groups according to the timing of voluntary finger movements. It appeared that these movements occurred at relatively regular intervals of approximately 10s in 5 subjects (group A); while 4 subjects showed random or very short inter-movement intervals (group B). Two remarkable results were obtained: first, the phase coupling (COH(2)) between BP and RRI showed a significant (p=0.0061) interaction between activity (rest vs. movement) and group (A vs. B), with an increased (p=0.0003) coupling in group A. Second, the COH(2) between BP and Hb oscillations showed a significant (p=0.034) interaction between activity and group, with a decreased (p=0.079) coupling in group B. These results suggest that subjects able to initiate self-paced, voluntary movements at relatively regular intervals of ∼10s show an entrainment potential between physiological oscillations and motor responses. This also provides the first evidence that not only physiological oscillations can be entrained to motor responses, but also motor responses (voluntary movements) can be entrained to slow intrinsic oscillations.

Persistence of Functional Sensory Maps in the Absence of Cortical Layers in the Somsatosensory Cortex of Reeler Mice

In rodents, layer IV of the primary somatosensory cortex contains the barrel field, where individual, large facial whiskers are represented as a dense cluster of cells. In the reeler mouse, a model of disturbed cortical development characterized by a loss of cortical lamination, the barrel field exists in a distorted manner. Little is known about the consequences of such a highly disturbed lamination on cortical function in this model. We used in vivo intrinsic signal optical imaging together with piezo-controlled whisker stimulation to explore sensory map organization and stimulus representation in the barrel field. We found that the loss of cortical layers in reeler mice had surprisingly little incidence on these properties. The overall topological order of whisker representations is highly preserved and the functional activation of individual whisker representations is similar in size and strength to wild-type controls. Because intrinsic imaging measures hemodynamic signals, we furthermore investigated the cortical blood vessel pattern of both genotypes, where we also did not detect major differences. In summary, the loss of the reelin protein results in a widespread disturbance of cortical development which compromises neither the establishment nor the function of an ordered, somatotopic map of the facial whiskers.

Spontaneous BOLD signal fluctuations in young healthy subjects and elderly patients with chronic kidney disease

Spontaneous fluctuations in blood oxygenation level-dependent (BOLD) images are the basis of resting-state fMRI and frequently used for functional connectivity studies. However, there may be intrinsic information in the amplitudes of these fluctuations. We investigated the possibility of using the amplitude of spontaneous BOLD signal fluctuations as a biomarker for cerebral vasomotor reactivity.

We compared the coefficient of variation (CV) of the time series (defined as the temporal standard deviation of the time series divided by the mean signal intensity) in two populations: 1) Ten young healthy adults and 2) Ten hypertensive elderly subjects with chronic kidney disease (CKD).

We found a statistically significant increase (P<0.01) in the CV values for the CKD patients compared with the young healthy adults in both gray matter (GM) and white matter (WM). The difference was independent of the exact segmentation method, became more significant after correcting for physiological signals using RETROICOR, and mainly arose from very low frequency components of the BOLD signal fluctuation (f<0.025 Hz). Furthermore, there was a strong relationship between WM and GM signal fluctuation CV's (R² = 0.87) in individuals, with a ratio of about 1∶3.

These results suggest that amplitude of the spontaneous BOLD signal fluctuations may be used to assess the cerebrovascular reactivity mechanisms and provide valuable information about variations with age and different disease states.

Spectral methods for functional brain imaging

Dynamic functional imaging experiments typically generate large, multivariate data sets that contain considerable spatial and temporal complexity. The goal of this introduction is to present signal-processing techniques that allow the underlying spatiotemporal structure to be readily distilled and that also enable signal versus noise contributions to be separated.

Direct, intraoperative observation of ~0.1 Hz hemodynamic oscillations in awake human cortex: implications for fMRI

An almost sinusoidal, large amplitude ~0.1 Hz oscillation in cortical hemodynamics has been repeatedly observed in species ranging from mice to humans. However, the occurrence of 'slow sinusoidal hemodynamic oscillations' (SSHOs) in human functional magnetic resonance imaging (fMRI) studies is rarely noted or considered. As a result, little investigation into the cause of SSHOs has been undertaken, and their potential to confound fMRI analysis, as well as their possible value as a functional biomarker has been largely overlooked. Here, we report direct observation of large-amplitude, sinusoidal ~0.1 Hz hemodynamic oscillations in the cortex of an awake human undergoing surgical resection of a brain tumor. Intraoperative multispectral optical intrinsic signal imaging (MS-OISI) revealed that SSHOs were spatially localized to distinct regions of the cortex, exhibited wave-like propagation, and involved oscillations in the diameter of specific pial arterioles, indicating that the effect was not the result of systemic blood pressure oscillations. fMRI data collected from the same subject 4 days prior to surgery demonstrates that ~0.1 Hz oscillations in the BOLD signal can be detected around the same region. Intraoperative optical imaging data from a patient undergoing epilepsy surgery, in whom sinusoidal oscillations were not observed, is shown for comparison. This direct observation of the '0.1 Hz wave' in the awake human brain, using both intraoperative imaging and pre-operative fMRI, confirms that SSHOs occur in the human brain, and can be detected by fMRI. We discuss the possible physiological basis of this oscillation and its potential link to brain pathologies, highlighting its relevance to resting-state fMRI and its potential as a novel target for functional diagnosis and delineation of neurological disease.

The cortical activation pattern by a rehabilitation robotic hand: a functional NIRS study

Introduction: Clarification of the relationship between external stimuli and brain response has been an important topic in neuroscience and brain rehabilitation. In the current study, using functional near infrared spectroscopy (fNIRS), we attempted to investigate cortical activation patterns generated during execution of a rehabilitation robotic hand.

Methods: Ten normal subjects were recruited for this study. Passive movements of the right fingers were performed using a rehabilitation robotic hand at a frequency of 0.5 Hz. We measured values of oxy-hemoglobin (HbO), deoxy-hemoglobin (HbR) and total-hemoglobin (HbT) in five regions of interest: the primary sensory-motor cortex (SM1), hand somatotopy of the contralateral SM1, supplementary motor area (SMA), premotor cortex (PMC), and prefrontal cortex (PFC).

Results: HbO and HbT values indicated significant activation in the left SM1, left SMA, left PMC, and left PFC during execution of the rehabilitation robotic hand (uncorrected, p < 0.01). By contrast, HbR value indicated significant activation only in the hand somatotopic area of the left SM1 (uncorrected, p < 0.01).

Conclusions: Our results appear to indicate that execution of the rehabilitation robotic hand could induce cortical activation.

Identifying and quantifying main components of physiological noise in functional near infrared spectroscopy on the prefrontal cortex

Functional Near-Infrared Spectroscopy (fNIRS) is a promising method to study functional organization of the prefrontal cortex. However, in order to realize the high potential of fNIRS, effective discrimination between physiological noise originating from forehead skin haemodynamic and cerebral signals is required. Main sources of physiological noise are global and local blood flow regulation processes on multiple time scales. The goal of the present study was to identify the main physiological noise contributions in fNIRS forehead signals and to develop a method for physiological de-noising of fNIRS data. To achieve this goal we combined concurrent time-domain fNIRS and peripheral physiology recordings with wavelet coherence analysis (WCA). Depth selectivity was achieved by analyzing moments of photon time-of-flight distributions provided by time-domain fNIRS. Simultaneously, mean arterial blood pressure (MAP), heart rate (HR), and skin blood flow (SBF) on the forehead were recorded. WCA was employed to quantify the impact of physiological processes on fNIRS signals separately for different time scales. We identified three main processes contributing to physiological noise in fNIRS signals on the forehead. The first process with the period of about 3 s is induced by respiration. The second process is highly correlated with time lagged MAP and HR fluctuations with a period of about 10 s often referred as Mayer waves. The third process is local regulation of the facial SBF time locked to the task-evoked fNIRS signals. All processes affect oxygenated haemoglobin concentration more strongly than that of deoxygenated haemoglobin. Based on these results we developed a set of physiological regressors, which were used for physiological de-noising of fNIRS signals. Our results demonstrate that proposed de-noising method can significantly improve the sensitivity of fNIRS to cerebral signals.

All-optical osteotomy to create windows for transcranial imaging in mice

Surgical procedures as a prelude to optical imaging are a rate-limiting step in experimental neuroscience. Towards automation of these procedures, we describe the use of nonlinear optical techniques to create a thinned skull window for transcranial imaging. Metrology by second harmonic generation was used to map the surfaces of the skull and define a cutting path. Plasma-mediated laser ablation was utilized to cut bone. Mice prepared with these techniques were used to image subsurface cortical vasculature and blood flow. The viability of the brain tissue was confirmed via histological analysis and supports the utility of solely optical techniques for osteotomy and potentially other surgical procedures.

Hemodynamic and electrophysiological spontaneous low-frequency oscillations in the cortex: directional influences revealed by Granger causality

We used a combined electrophysiological/hemodynamic system to examine low-frequency oscillations (LFOs) in spontaneous neuronal activities (spike trains and local field potentials) and hemodynamic signals (cerebral blood flow) recorded from the anesthetized rat somatosensory and visual cortices. The laser Doppler flowmetry (LDF) probe was tilted slightly to approach the area in which a microelectrode array (MEA) was implanted for simultaneous recordings. Spike trains (STs) were converted into continuous-time rate functions (CRFs) using the ST instantaneous firing rates. LFOs were detected for all three of the components using the multi-taper method (MTM). The frequencies of these LFOs ranged from 0.052 to 0.167 Hz (mean±SD, 0.10±0.026 Hz) for cerebral blood flow (CBF), from 0.027 to 0.26 Hz (mean±SD, 0.12±0.041 Hz) for the CRFs of the STs and from 0.04 to 0.19 Hz (mean±SD, 0.11±0.035 Hz) for local field potentials (LFPs). We evaluated the Granger causal relationships of spontaneous LFOs among CBF, LFPs and CRFs using Granger causality (GC) analysis. Significant Granger causal relationships were observed from LFPs to CBF, from STs to CBF and from LFPs to STs at approximately 0.1 Hz. The present results indicate that spontaneous LFOs exist not only in hemodynamic components but also in neuronal activities of the rat cortex. To the best of our knowledge, the present study is the first to identify Granger causal influences among CBF, LFPs and STs and show that spontaneous LFOs carry important Granger causal influences from neural activities to hemodynamic signals.

PHYCAA+: an optimized, adaptive procedure for measuring and controlling physiological noise in BOLD fMRI

The presence of physiological noise in functional MRI can greatly limit the sensitivity and accuracy of BOLD signal measurements, and produce significant false positives. There are two main types of physiological confounds: (1) high-variance signal in non-neuronal tissues of the brain including vascular tracts, sinuses and ventricles, and (2) physiological noise components which extend into gray matter tissue. These physiological effects may also be partially coupled with stimuli (and thus the BOLD response). To address these issues, we have developed PHYCAA+, a significantly improved version of the PHYCAA algorithm (Churchill et al., 2011) that (1) down-weights the variance of voxels in probable non-neuronal tissue, and (2) identifies the multivariate physiological noise subspace in gray matter that is linked to non-neuronal tissue. This model estimates physiological noise directly from EPI data, without requiring external measures of heartbeat and respiration, or manual selection of physiological components. The PHYCAA+ model significantly improves the prediction accuracy and reproducibility of single-subject analyses, compared to PHYCAA and a number of commonly-used physiological correction algorithms. Individual subject denoising with PHYCAA+ is independently validated by showing that it consistently increased between-subject activation overlap, and minimized false-positive signal in non gray-matter loci. The results are demonstrated for both block and fast single-event task designs, applied to standard univariate and adaptive multivariate analysis models.

Spontaneous oscillations in intrinsic signals reveal the structure of cerebral vasculature

Functional imaging of intrinsic signals allows minimally invasive spatiotemporal mapping of stimulus representations in the cortex, but representations are often corrupted by stimulus-independent spatial artifacts, especially those originating from the blood vessels. In this paper, we present novel algorithms for unsupervised identification of cerebral vascularization, allowing blind separation of stimulus representations from noise. These algorithms commonly take advantage of the temporal fluctuations in global reflectance to extract anatomic information. More specifically, the phase of low-frequency oscillations relative to global fluctuations reveals local vascular identity. Arterioles can be reconstructed using their characteristically high power in those frequencies corresponding to respiration, heartbeat, and vasomotion signals. By treating the vasculature as a dynamic flow network, we finally demonstrate that direction of blood perfusion can be quantitatively visualized. Application of these methods for removal of stimulus-independent changes in reflectance permits isolation of stimulus-evoked representations even if the representation spatially overlaps with blood vessels. The algorithms can be expanded further to extract temporal information on blood flow, monitor revascularization following a focal stroke, and distinguish arterioles from venules and parenchyma.

Inter-trial variability in sensory-evoked cortical hemodynamic responses: the role of the magnitude of pre-stimulus fluctuations

Brain imaging techniques utilize hemodynamic changes that accompany brain activation. However, stimulus-evoked hemodynamic responses display considerable inter-trial variability and the sources of this variability are poorly understood. One of the sources of this response variation could be ongoing spontaneous hemodynamic fluctuations. We recently investigated this issue by measuring cortical hemodynamics in response to sensory stimuli in anesthetized rodents using 2-dimensional optical imaging spectroscopy. We suggested that sensory-evoked cortical hemodynamics displayed distinctive response characteristics and magnitudes depending on the phase of ongoing fluctuations at stimulus onset due to a linear superposition of evoked and ongoing hemodynamics (Saka et al., 2010). However, the previous analysis neglected to examine the possible influence of variability of the size of ongoing fluctuations. Consequently, data were further analyzed to examine whether the size of pre-stimulus hemodynamic fluctuations also influenced the magnitude of subsequent stimulus-evoked responses. Indeed, in the case of all individual trials, a moderate correlation between the size of the pre-stimulus fluctuations and the magnitudes of the subsequent sensory-evoked responses were observed. However, different correlations between the size of the pre-stimulus fluctuations and magnitudes of the subsequent sensory-evoked cortical hemodynamic responses could be observed depending on their phase at stimulus onset. These analyses suggest that both the size and phase of pre-stimulus fluctuations in cortical hemodynamics contribute to inter-trial variability in sensory-evoked responses.

Surrounding inhibition in rat somatosensory cortex during noxious electrical stimulation of the sciatic nerve

Surrounding inhibition is a physiologic mechanism to focus neuronal activity. Here we applied optical imaging of intrinsic signal to observing the temporal-spatial characteristic of rat primary somatosensory cortex during graded electrical stimulation of the sciatic nerve (5 Hz, duration of 2 s, 0.5 ms pulse, 1 and 10 muscle twitching threshold). We found that the magnitude and change duration (time course) of the optical signal were larger and longer with the intensity raising. The spatial extent was also wider under noxious electrical stimulus. Meanwhile, we found the inverted optical signal changes in the regions surround the activated primary somatosensory cortex. This phenomenon was similar to the inhibition surrounding focal itcal events observed by optical imaging of intrinsic signal. It suggests the surrounding inhibition under noxious electrical stimulus was probably induced by the excess discharge of excited neurons or it may provide finer discrimination during the noxious stimulus and support the view that the role of somatosensory cortex in pain localization is to finely discriminate the stimulus site.

Oscillations in cerebral haemodynamics in patients with falciparum malaria

Spontaneous oscillations in cerebral haemodynamics studied with near-infrared spectroscopy (NIRS), become impaired in several pathological conditions. We assessed the spectral characteristics of these oscillations in 20 patients with falciparum malaria admitted to Ispat General Hospital, Rourkela, India. Monitoring included continuous frontal lobe NIRS recordings within 24 h of admission (Day 0), together with single measurements of a number of clinical and chemical markers recorded on admission. Seven patients returned for follow-up measurements on recovery (FU). A 2,048 sampling-point segment of oxygenated haemoglobin concentration ([ΔHbO(2)]) data was subjected to Fourier analysis per patient, and power spectral density was derived over the very low frequency (VLF: 0.02-0.04 Hz), low frequency (LF: 0.04-0.15 Hz) and high frequency (HF: 0.15-0.4 Hz) bands. At Day 0, VLF spectral power was 21.1 ± 16.4, LF power 7.2 ± 4.6 and HF power 2.6 ± 5.0, with VLF power being statistically significantly higher than LF and HF (P < 0.005). VLF power tended to decrease in the severely ill patients and correlated negatively with heart rate (r = 0.57, P < 0.01), while LF power correlated positively with aural body temperature (r = 0.49, P < 0.05). In all but one of the patients who returned for FU measurements, VLF power increased after recovery. This may be related to autonomic dysfunction in severe malaria, a topic of little research to date. The present study demonstrated that application of NIRS in a resource-poor setting is feasible and has potential as a research tool.

Phase-amplitude investigation of spontaneous low-frequency oscillations of cerebral hemodynamics with near-infrared spectroscopy: a sleep study in human subjects

We have investigated the amplitude and phase of spontaneous low-frequency oscillations (LFOs) of the cerebral deoxy- and oxy-hemoglobin concentrations ([Hb] and [HbO]) in a human sleep study using near-infrared spectroscopy (NIRS). Amplitude and phase analysis was based on the analytic signal method, and phasor algebra was used to decompose measured [Hb] and [HbO] oscillations into cerebral blood volume (CBV) and flow velocity (CBFV) oscillations. We have found a greater phase lead of [Hb] vs. [HbO] LFOs during non-REM sleep with respect to the awake and REM sleep states (maximum increase in [Hb] phase lead: ~π/2). Furthermore, during non-REM sleep, the amplitudes of [Hb] and [HbO] LFOs are suppressed with respect to the awake and REM sleep states (maximum amplitude decrease: 87%). The associated cerebral blood volume and flow velocity oscillations are found to maintain their relative phase difference during sleep, whereas their amplitudes are attenuated during non-REM sleep. These results show the potential of phase-amplitude analysis of [Hb] and [HbO] oscillations measured by NIRS in the investigation of hemodynamics associated with cerebral physiology, activation, and pathological conditions.

Two-photon microscopy as a tool to study blood flow and neurovascular coupling in the rodent brain

The cerebral vascular system services the constant demand for energy during neuronal activity in the brain. Attempts to delineate the logic of neurovascular coupling have been greatly aided by the advent of two-photon laser scanning microscopy to image both blood flow and the activity of individual cells below the surface of the brain. Here we provide a technical guide to imaging cerebral blood flow in rodents. We describe in detail the surgical procedures required to generate cranial windows for optical access to the cortex of both rats and mice and the use of two-photon microscopy to accurately measure blood flow in individual cortical vessels concurrent with local cellular activity. We further provide examples on how these techniques can be applied to the study of local blood flow regulation and vascular pathologies such as small-scale stroke.

Color blobs in cortical areas V1 and V2 of the new world monkey Callithrix jacchus, revealed by non-differential optical imaging

Color vision is reserved to only few mammals, such as Old World monkeys and humans. Most Old World monkeys are trichromats. Among them, macaques were shown to exhibit functional domains of color-selectivity, in areas V1 and V2 of the visual cortex. Such color domains have not yet been shown in New World monkeys. In marmosets a sex-linked dichotomy results in dichromatic and trichromatic genotypes, rendering most male marmosets color-blind. Here we used trichromatic female marmosets to examine the intrinsic signal response in V1 and V2 to chromatic and achromatic stimuli, using optical imaging. To activate the subsystems individually, we used spatially homogeneous isoluminant color opponent (red/green, blue/yellow) and hue versus achromatic flicker (red/gray, green/gray, blue/gray, yellow/gray), as well as achromatic luminance flicker. In contrast to previous optical imaging studies in marmosets, we find clearly segregated color domains, similar to those seen in macaques. Red/green and red/gray flicker were found to be the appropriate stimulus for revealing color domains in single-condition maps. Blue/gray and blue/yellow flicker stimuli resulted in faint patch-patterns. A recently described multimodal vessel mapping approach allowed for an accurate alignment of the functional and anatomical datasets. Color domains were tightly colocalized with cytochrome oxidase blobs in V1 and with thin stripes in V2. Thus, our findings are in accord with 2-Deoxy-D-glucose studies performed in V1 of macaques and studies on color representation in V2. Our results suggest a similar organization of early cortical color processing in trichromats of both Old World and New World monkeys.

Spatial summation of neurometabolic coupling in the central visual pathway

Noninvasive neural imaging has become an important tool in both applied and theoretical applications. The hemodynamic properties that are measured in functional magnetic resonance imaging (fMRI), for example, are generally used to infer neuronal characteristics. In an attempt to provide empirical data to connect the hemodynamic measurements with neural function, we have conducted previous studies in which neural activity and tissue oxygen metabolic functions are determined together in co-localized regions of the central visual pathway. A basic question in this procedure is whether oxygen responses are coupled linearly in space and time with neural activity. We have previously examined temporal factors, and in the current study, spatial characteristics are addressed. We have recorded from neurons in the lateral geniculate nucleus (LGN) and striate cortex in anesthetized cats. In both structures, there is a classical receptive field (CRF) within which a neuron can be activated. There is also a region outside the CRF from which stimulation cannot activate the cell directly but can influence the response elicited from the CRF. In this investigation we have used several specific spatial stimulus patterns presented to either the CRF or the surrounding region or to both areas together in order to determine spatial response patterns. Within the CRF, we find that neural and metabolic responses sum in a nonlinear fashion but changes in these two measurements are closely coupled. For stimuli that extend beyond the CRF, neural activity is generally reduced while oxygen response exhibits uncoupled changes.

Individual variability in the shape and amplitude of the BOLD-HRF correlates with endogenous GABAergic inhibition

It has previously been demonstrated that there is a negative correlation between the amplitude of the BOLD response and resting γ amino-butyric acid (GABA) concentration in visual cortex. The work here is the first to empirically characterize individual variability in the haemodynamic response functions (HRFs) in response to a simple visual stimulus and baseline GABA concentration in a population of young adult males (n = 15, aged 20-28 years). The results demonstrate that GABA concentration is negatively correlated with BOLD response amplitude (r = -0.64, P < 0.02) and positively correlated with HRF width (r = 0.67, P < 0.002), that is, individuals with higher resting GABA concentration tend to exhibit smaller and wider HRFs. No correlations were observed with resting cerebral blood flow and GABA concentration and similarly, no correlations were observed between GABA and the proportional tissue content of the MRS voxel. We argue that correlation of the height of the HRF is supportive of the view that the previously observed correlations between BOLD amplitudes and GABA are reflective of differences in neuronal activity. However, the changes in HRF shape in individuals with higher baseline GABA levels are suggestive that differing vascular response characteristics may also make a significant contribution. Our results reinforce the view that variability in endogenous factors, such as neurotransmitter concentration, can have a profound effect on the vascular haemodynamic response. This has important implications for between-cohort fMRI studies in which variation in parameters such as GABA concentration may lead to group differences in the BOLD signal.

Vasomotion and neurovascular coupling in the visual thalamus in vivo

Spontaneous contraction and relaxation of arteries (and in some instances venules) has been termed vasomotion and has been observed in an extensive variety of tissues and species. However, its functions and underlying mechanisms are still under discussion. We demonstrate that in vivo spectrophotometry, measured simultaneously with extracellular recordings at the same locations in the visual thalamus of the cat, reveals vasomotion, measured as an oscillation (0.14hz) in the recorded oxyhemoglobin (OxyHb) signal, which appears spontaneously in the microcirculation and can last for periods of hours. During some non-oscillatory periods, maintained sensory stimulation evokes vasomotion lasting ∼30s, resembling an adaptive vascular phenomenon. This oscillation in the oxyhaemoblobin signal is sensitive to pharmacological manipulation: it is inducible by chloralose anaesthesia and it can be temporarily blocked by systemic administration of adrenaline or acetylcholine (ACh). During these oscillatory periods, neurovascular coupling (i.e. the relationship between local neural activity and the rate of blood supply to that location) appears significantly altered. This raises important questions with regard to the interpretation of results from studies currently dependent upon a linear relationship between neural activity and blood flow, such as neuroimaging.