Spatial specificity of the enhanced dip inherently induced by prolonged oxygen consumption in cat visual cortex: Implication for columnar resolution functional MRI

Time-dependent spatial specificity of high-resolution fMRI: insights into mesoscopic neurovascular coupling

High-resolution functional magnetic resonance imaging (fMRI) is becoming increasingly popular because of the growing availability of ultra-high magnetic fields which are capable of improving sensitivity and spatial resolution. However, it is debatable whether increased spatial resolutions for haemodynamic-based techniques, like fMRI, can accurately detect the true location of neuronal activity. We have addressed this issue in functional columns and layers of animals with haemoglobin-based optical imaging and different fMRI contrasts, such as blood oxygenation level-dependent, cerebral blood flow and cerebral blood volume fMRI. In this review, we describe empirical evidence primarily from our own studies on how well these fMRI signals are spatially specific to the neuronally active site and discuss insights into neurovascular coupling at the mesoscale. This article is part of the theme issue 'Key relationships between non-invasive functional neuroimaging and the underlying neuronal activity'.

Existence of Initial Dip for BCI: An Illusion or Reality

A tight coupling between the neuronal activity and the cerebral blood flow (CBF) is the motivation of many hemodynamic response (HR)-based neuroimaging modalities. The increase in neuronal activity causes the increase in CBF that is indirectly measured by HR modalities. Upon functional stimulation, the HR is mainly categorized in three durations: (i) initial dip, (ii) conventional HR (i.e., positive increase in HR caused by an increase in the CBF), and (iii) undershoot. The initial dip is a change in oxygenation prior to any subsequent increase in CBF and spatially more specific to the site of neuronal activity. Despite additional evidence from various HR modalities on the presence of initial dip in human and animal species (i.e., cat, rat, and monkey); the existence/occurrence of an initial dip in HR is still under debate. This article reviews the existence and elusive nature of the initial dip duration of HR in intrinsic signal optical imaging (ISOI), functional magnetic resonance imaging (fMRI), and functional near-infrared spectroscopy (fNIRS). The advent of initial dip and its elusiveness factors in ISOI and fMRI studies are briefly discussed. Furthermore, the detection of initial dip and its role in brain-computer interface using fNIRS is examined in detail. The best possible application for the initial dip utilization and its future implications using fNIRS are provided.

Biophysics of BOLD fMRI investigated with animal models

The widely-used BOLD fMRI signal depends on various anatomical, physiological, and imaging parameters. Thus, it is important to examine its biophysical and physiological source in order to optimize, model and accurately interpret fMRI. Animal models have been used to investigate these issues to take systematic measurements and combine with conventional invasive approaches. Here, we reviewed and discussed multiple issues, including the echo time-dependent intravascular contribution and extravascular contributions, gradient-echo vs. spin-echo fMRI, the physiological source of BOLD fMRI, arterial vs. venous cerebral blood volume change, cerebral oxygen consumption change, and arterial oxygen saturation change. We then discuss future directions of animal fMRI and translation to human fMRI. Systematic biophysical BOLD fMRI studies provide insight into the modeling and interpretation of BOLD fMRI in animals and humans.

Optogenetic investigation of the variable neurovascular coupling along the interhemispheric circuits

The interhemispheric circuit connecting the left and the right mammalian brain plays a key role in integration of signals from the left and the right side of the body. The information transfer is carried out by modulation of simultaneous excitation and inhibition. Hemodynamic studies of this circuit are inconsistent since little is known about neurovascular coupling of mixed excitatory and inhibitory signals. We investigated the variability in hemodynamic responses driven by the interhemispheric circuit during optogenetic and somatosensory activation. We observed differences in the neurovascular response based on the stimulation site - cell bodies versus distal projections. In half of the experiments, optogenetic stimulation of the cell bodies evoked a predominant post-synaptic inhibition in the other hemisphere, accompanied by metabolic oxygen consumption without coupled functional hyperemia. When the same transcallosal stimulation resulted in predominant post-synaptic excitation, the hemodynamic response was biphasic, consisting of metabolic dip followed by functional hyperemia. Optogenetic suppression of the postsynaptic excitation abolished the coupled functional hyperemia. In contrast, light stimulation at distal projections evoked consistently a metabolic response. Our findings suggest that functional hyperemia requires signals originating from the cell body and the hemodynamic response variability appears to reflect the balance between the post-synaptic excitation and inhibition.

Precision mapping of the vibrissa representation within murine primary somatosensory cortex

The ability to form an accurate map of sensory input to the brain is an essential aspect of interpreting functional brain signals. Here, we consider the somatotopic map of vibrissa-based touch in the primary somatosensory (vS1) cortex of mice. The vibrissae are represented by a Manhattan-like grid of columnar structures that are separated by inter-digitating septa. The development, dynamics and plasticity of this organization is widely used as a model system. Yet, the exact anatomical position of this organization within the vS1 cortex varies between individual mice. Targeting of a particular column in vivo therefore requires prior mapping of the activated cortical region, for instance by imaging the evoked intrinsic optical signal (eIOS) during vibrissa stimulation. Here, we describe a procedure for constructing a complete somatotopic map of the vibrissa representation in the vS1 cortex using eIOS. This enables precise targeting of individual cortical columns. We found, using C57BL/6 mice, that although the precise location of the columnar field varies between animals, the relative spatial arrangement of the columns is highly preserved. This finding enables us to construct a canonical somatotopic map of the vibrissae in the vS1 cortex. In particular, the position of any column, in absolute anatomical coordinates, can be established with near certainty when the functional representations in the vS1 cortex for as few as two vibrissae have been mapped with eIOS.This article is part of the themed issue 'Interpreting BOLD: a dialogue between cognitive and cellular neuroscience'.

Imaging the Neocortex Functional Architecture Using Multiple Intrinsic Signals: Implications for Hemodynamic-Based Functional Imaging

Optical imaging based on intrinsic signals has provided a new level of understanding of the principles underlying cortical development, organization, and function, providing a spatial resolution of up to 20 µm for mapping cortical columns in vivo. This introduction briefly reviews the development of this technique, the types of applications that have been pursued, and the general implications of some findings for other neuroimaging techniques based on hemodynamic responses (e.g., functional magnetic resonance imaging).

Spatiotemporal characteristics and vascular sources of neural-specific and -nonspecific fMRI signals at submillimeter columnar resolution

The neural specificity of hemodynamic-based functional magnetic resonance imaging (fMRI) signals is dependent on both the vascular regulation and the sensitivity of the applied fMRI technique to different types and sizes of blood vessels. In order to examine the specificity of MRI-detectable hemodynamic responses, submillimeter blood oxygenation level-dependent (BOLD) and cerebral blood volume (CBV) fMRI studies were performed in a well-established cat orientation column model at 9.4 T. Neural-nonspecific and -specific signals were separated by comparing the fMRI responses of orthogonal orientation stimuli. The BOLD response was dominantly neural-nonspecific, mostly originating from pial and intracortical emerging veins, and thus was highly correlated with baseline blood volume. Uneven baseline CBV may displace or distort small functional domains in high-resolution BOLD maps. The CBV response in the parenchyma exhibited dual spatiotemporal characteristics, a fast and early neural-nonspecific response (with 4.3-s time constant) and a slightly slower and delayed neural-specific response (with 9.4-s time constant). The nonspecific CBV signal originates from early-responding arteries and arterioles, while the specific CBV response, which is not correlated with baseline blood volume, arises from late-responding microvessels including small pre-capillary arterioles and capillaries. Our data indicate that although the neural specificity of CBV fMRI signals is dependent on stimulation duration, high-resolution functional maps can be obtained from steady-state CBV studies.

Biophysical and physiological origins of blood oxygenation level-dependent fMRI signals

After its discovery in 1990, blood oxygenation level-dependent (BOLD) contrast in functional magnetic resonance imaging (fMRI) has been widely used to map brain activation in humans and animals. Since fMRI relies on signal changes induced by neural activity, its signal source can be complex and is also dependent on imaging parameters and techniques. In this review, we identify and describe the origins of BOLD fMRI signals, including the topics of (1) effects of spin density, volume fraction, inflow, perfusion, and susceptibility as potential contributors to BOLD fMRI, (2) intravascular and extravascular contributions to conventional gradient-echo and spin-echo BOLD fMRI, (3) spatial specificity of hemodynamic-based fMRI related to vascular architecture and intrinsic hemodynamic responses, (4) BOLD signal contributions from functional changes in cerebral blood flow (CBF), cerebral blood volume (CBV), and cerebral metabolic rate of O(2) utilization (CMRO(2)), (5) dynamic responses of BOLD, CBF, CMRO(2), and arterial and venous CBV, (6) potential sources of initial BOLD dips, poststimulus BOLD undershoots, and prolonged negative BOLD fMRI signals, (7) dependence of stimulus-evoked BOLD signals on baseline physiology, and (8) basis of resting-state BOLD fluctuations. These discussions are highly relevant to interpreting BOLD fMRI signals as physiological means.

Image-based vessel-by-vessel analysis for red blood cell and plasma dynamics with automatic segmentation

The aim of the present study was to test the hypothesis that vascular tones of cortical surface and parenchymal blood flow can be dissociated depending on the perturbation. To this end, a novel image-based analytical method for quantitatively measuring vessel diameters and flow dynamics was developed. The algorithm relies on the spatiotemporal coherence of the pixel intensity changes induced by the transit of the fluorescent signals measured using confocal laser scanning fluorescent microscopy in the rat cerebral cortex. A cocktail of fluorescently labeled red blood cell (RBC) and plasma agents was administered to simultaneously compare RBC and plasma dynamics in the same vascular networks. The time to fluorescent signal appearance and the width of the fluorescent signal were measured in each segment and compared between sodium nitroprusside-induced global and sensory stimulation-induced local perturbation conditions. We observed that infusion of sodium nitroprusside induced significant vasodilation in the surface artery, particularly in the small arteries (1.8-fold increase). Vasodilation induced by sensory stimulation was observed to depend on vessel size, but significant changes were only detected for the small arteries and veins. Measurements of the time to venous appearance revealed that appearance time was extended by sodium nitroprusside, but shortened during forepaw stimulation, relative to the control condition. Both perturbations provoked the largest changes between the small artery and vein segments, indicating that the changes in the appearance time originate from blood passage through parenchymal microcirculation. These findings support the hypothesis that cortical surface vascular tone and parenchymal blood flow are individually coordinated.

Evolution of the dynamic changes in functional cerebral oxidative metabolism from tissue mitochondria to blood oxygen

The dynamic properties of the cerebral metabolic rate of oxygen consumption (CMR(O2)) during changes in brain activity remain unclear. Therefore, the spatial and temporal evolution of functional increases in CMR(O2) was investigated in the rat somato-sensory cortex during forelimb stimulation under a suppressed blood flow response condition. Temporally, stimulation elicited a fast increase in tissue mitochondria CMR(O2) described by a time constant of ~1 second measured using flavoprotein autofluorescence imaging. CMR(O2)-driven changes in the tissue oxygen tension measured using an oxygen electrode and blood oxygenation measured using optical imaging of intrinsic signal followed; however, these changes were slow with time constants of ~5 and ~10 seconds, respectively. This slow change in CMR(O2)-driven blood oxygenation partly explains the commonly observed post-stimulus blood oxygen level-dependent (BOLD) undershoot. Spatially, the changes in mitochondria CMR(O2) were similar to the changes in blood oxygenation. Finally, the increases in CMR(O2) were well correlated with the evoked multi-unit spiking activity. These findings show that dynamic CMR(O2) calculations made using only blood oxygenation data (e.g., BOLD functional magnetic resonance imaging (fMRI)) do not directly reflect the temporal changes in the tissue's mitochondria metabolic rate; however, the findings presented can bridge the gap between the changes in cellular oxidative rate and blood oxygenation.

The story of the initial dip in fMRI

Over the past 20 years much attention has been given to characterizing the spatial accuracy of fMRI based signals and to techniques that improve on its co-localization with neuronal activity. While the vast majority of fMRI studies have always used the conventional positive BOLD signal, alternative contrast options have demonstrated superior spatial specificity. One of these options surfaced shortly after the initial BOLD fMRI demonstrations and was motivated by optical imaging studies which revealed an early signal change that was much smaller but spatially more specific than the delayed positive response. This early signal change was attributed to oxygenation changes prior to any subsequent blood flow increases. After observation of this biphasic hemodynamic response in fMRI, because this early response resulted in a small MR signal decrease prior to the onset of the large signal increase, it became known as the "initial dip". While the initial dip in fMRI was subsequently reported by many studies, including those in humans, monkeys, and cats, there were conflicting views about the associated mechanisms and whether it could be generalized across brain regions or species, in addition to whether or not it would prove fruitful for neuroscience. These discrepancies, along with the implications that the initial dip might increase the spatial specificity of BOLD fMRI from 2 to 3mm to something more closely associated with neural activity, resulted in lot of buzz and controversy in the community for many years. In this review, the authors provide an account of the story of the initial dip in MR based functional imaging from the Minnesota perspective, where the first demonstrations, characterizations, and applications of the initial dip commenced.

Differential fMRI activation to noxious heat and tactile stimuli in parasylvian areas of new world monkeys

Emerging evidence supports an important role of posterior parasylvian areas in both pain and touch processing. Whether there are separate or shared networks for these sensations remains controversial. The present study compared spatial patterns of brain activation in response to unilateral nociceptive heat (47.5°C) or innocuous tactile stimulation (8-Hz vibration) to digits through high-resolution functional magnetic resonance imaging (fMRI) in squirrel monkeys. In addition, the temporal profile of heat-stimulus-evoked fMRI Blood Oxygenation Level Dependent (BOLD) signal changes was characterized. By examining high-resolution fMRI and histological measures at both the individual and the group levels, we found that both nociceptive heat and tactile stimuli elicited activation in bilateral secondary somatosensory and ventral parietal areas (S2/PV) and in ipsilateral ventral somatosensory areas (VS) and retroinsula (Ri). Bilateral posterior insular cortex (pIns) and area 7b responded preferentially to nociceptive heat stimulation. Single voxels within each activation cluster showed robust BOLD signal changes during each block of nociceptive stimulation. Across animals (n=11), nociceptive response magnitudes of contralateral VS and pIns and ipsilateral Ri were significantly greater than corresponding areas in the opposite hemisphere. In sum, both distinct and shared areas in regions surrounding the posterior sylvian fissure were activated in response to nociceptive and tactile inputs in nonhuman primates.

Cerebral oxygen delivery and consumption during evoked neural activity

Increases in neural activity evoke increases in the delivery and consumption of oxygen. Beyond observations of cerebral tissue and blood oxygen, the role and properties of cerebral oxygen delivery and consumption during changes in brain function are not well understood. This work overviews the current knowledge of functional oxygen delivery and consumption and introduces recent and preliminary findings to explore the mechanisms by which oxygen is delivered to tissue as well as the temporal dynamics of oxygen metabolism. Vascular oxygen tension measurements have shown that a relatively large amount of oxygen exits pial arterioles prior to capillaries. Additionally, increases in cerebral blood flow (CBF) induced by evoked neural activation are accompanied by arterial vasodilation and also by increases in arteriolar oxygenation. This increase contributes not only to the down-stream delivery of oxygen to tissue, but also to delivery of additional oxygen to extra-vascular spaces surrounding the arterioles. On the other hand, the changes in tissue oxygen tension due to functional increases in oxygen consumption have been investigated using a method to suppress the evoked CBF response. The functional decreases in tissue oxygen tension induced by increases in oxygen consumption are slow to evoked changes in CBF under control conditions. Preliminary findings obtained using flavoprotein autofluorescence imaging suggest cellular oxidative metabolism changes at a faster rate than the average changes in tissue oxygen. These issues are important in the determination of the dynamic changes in tissue oxygen metabolism from hemoglobin-based imaging techniques such as blood oxygenation-level dependent functional magnetic resonance imaging (fMRI).

Changes in cerebral arterial, tissue and venous oxygenation with evoked neural stimulation: implications for hemoglobin-based functional neuroimaging

Little is known regarding the changes in blood oxygen tension (P(O2)) with changes in brain function. This work aimed to measure the blood P(O2) in surface arteries and veins as well as tissue with evoked somato-sensory stimulation in the anesthetized rat. Electrical stimulation of the forepaw induced average increases in blood flow of 44% as well as increases in the tissue P(O2) of 28%. More importantly, increases in P(O2) throughout pial arteries (resting diameters=59 to 129 microm) and pial veins (resting diameters=62 to 361 microm) were observed. The largest increases in vascular P(O2) were observed in the small veins (from 33 to 40 mm Hg) and small arteries (from 78 to 88 mm Hg). The changes in oxygen saturation (S(O2)) were calculated and the largest increases were observed in small veins (Delta=+11%) while its increase in small arteries was small (Delta=+4%). The average diameter of arterial vessels was observed to increase by 4 to 6% while that of veins was not observed to change with evoked stimulation. These findings show that the increases in arterial P(O2) contribute to the hyper-oxygenation of tissue and, mostly likely, also to the signal changes in hemoglobin-based functional imaging methods (e.g. BOLD fMRI).

Functional MRI of postnatal visual development in normal and hypoxic-ischemic-injured superior colliculi

The superior colliculus (SC) is a laminated subcortical structure in the mammalian midbrain, whose superficial layers receive visual information from the retina and the visual cortex. To date, its functional organization and development in the visual system remain largely unknown. This study employed blood oxygenation level-dependent (BOLD) functional MRI to evaluate the visual responses of the SC in normally developing and severe neonatal hypoxic-ischemic (HI)-injured rat brains from the time of eyelid opening to adulthood. MRI was performed to the normal animals (n=7) at postnatal days (P) 14, 21, 28 and 60. In the HI-injured group (n=7), the ipsilesional primary and secondary visual cortices were completely damaged after unilateral ligation of the left common carotid artery at P7 followed by hypoxia for 2 h, and MRI was performed at P60. Upon unilateral flash illumination, the normal contralateral SC underwent a systematic increase in BOLD signal amplitude with age especially after the third postnatal week. However, no significant difference in BOLD signal increase was found between P14 and P21. These findings implied the presence of neurovascular coupling at the time of eyelid opening, and the progressive development of hemodynamic regulation in the subcortical visual system. In the HI-injured group at P60, the BOLD signal increases in both SC remained at the same level as the normal group at P28 though they were significantly lower than the normal group at P60. These observations suggested the residual visual functions on both sides of the subcortical brain, despite the damages to the entire ipsilesional visual cortex. The results of this study constitute important evidence on the progressive maturation of visual functions and hemodynamic responses in the normal subcortical brain, and its functional plasticity upon neonatal HI injury.

Oxygen transport in brain tissue

Oxygen is essential to maintaining normal brain function. A large body of evidence suggests that the partial pressure of oxygen (pO(2)) in brain tissue is physiologically maintained within a narrow range in accordance with region-specific brain activity. Since the transportation of oxygen in the brain tissue is mainly driven by a diffusion process caused by a concentration gradient of oxygen from blood to cells, the spatial organization of the vascular system, in which the oxygen content is higher than in tissue, is a key factor for maintaining effective transportation. In addition, a local mechanism that controls energy demand and blood flow supply plays a critical role in moment-to-moment adjustment of tissue pO(2) in response to dynamically varying brain activity. In this review, we discuss the spatiotemporal structures of brain tissue oxygen transport in relation to local brain activity based on recent reports of tissue pO(2) measurements with polarographic oxygen microsensors in combination with simultaneous recordings of neural activity and local cerebral blood flow in anesthetized animal models. Although a physiological mechanism of oxygen level sensing and control of oxygen transport remains largely unknown, theoretical models of oxygen transport are a powerful tool for better understanding the short-term and long-term effects of local changes in oxygen demand and supply. Finally, emerging new techniques for three-dimensional imaging of the spatiotemporal dynamics of pO(2) map may enable us to provide a whole picture of how the physiological system controls the balance between demand and supply of oxygen during both normal and pathological brain activity.

Is cortical vasculature functionally organized?

The cortical vasculature is a well-structured and organized system, but the extent to which it is organized with respect to the neuronal functional architecture is unknown. In particular, does vasculature follow the same functional organization as cortical columns? In principle, cortical columns that share tuning for stimulus features like orientation may often be active together and thus require oxygen and metabolic nutrients together. If the cortical vasculature is built to serve these needs, it may also tend to aggregate and amplify orientation specific signals and explain why they are available in fMRI data at very low resolution.

Hemodynamic surrogates for excitatory membrane potential change during interictal epileptiform events in rat neocortex

Hemodynamic changes in the brain are often used as surrogates for epileptic neuronal activity in both the laboratory and the clinic (e.g., intrinsic signal, functional magnetic resonance imaging and single-photon emission computed tomography) in spite of the fact that perfusion-based signals have been shown to overestimate the population of spiking neurons. In addition, mechanisms of neurovascular coupling that apply during normal cortical processing may not be relevant in pathological circumstances such as epilepsy. For these reasons, we investigated the spatiotemporal dynamics of epileptic neurovascular coupling using voltage-sensitive dyes (VSDs) to generate spatial maps of excitatory membrane activity and intrinsic optical spectroscopy (IOS) to measure deoxy-hemoglobin and total hemoglobin, i.e., cerebral blood volume (CBV), in vivo during interictal spikes in rat neocortex to examine their spatiotemporal correlations. We hypothesized that the IOS signal would correlate spatially with subthreshold excitatory activity, which involves a larger area of cortex than suprathreshold neuronal spiking. However, we found that both perfusion and oximetric signals spatially overshot the extent of the excitatory VSD signal by approximately 2x. Nevertheless, a high correlation could be found at specific time points in the evolution and dissolution of the hemodynamic signals. The increase in deoxy-hemoglobin reached the highest correlation with the excitatory VSD signal earlier than CBV signals although CBV signals correlated equally well at certain time points. The amplitude of the hemodynamic signals had a linear correlation with the amplitude of the VSD signals except for small nonlinearities in the very center of the focus and in the periphery of the surround, indicating a tight spatial coupling. Our data suggest that hemodynamic signals can accurately define the spatial extent of excitatory interictal epileptiform subthreshold membrane activity at specific time points in their evolution.

Brain tissue oxygen consumption and supply induced by neural activation: determined under suppressed hemodynamic response conditions in the anesthetized rat cerebral cortex

The dynamic changes in cerebral metabolic rate of oxygen (CMRO2) and oxygen supply during brain functions have not been well-characterized. To examine this issue, experiments with electrophysiology, oxygen microelectrode and laser-Doppler flowmetry were performed in the anesthetized rat somatosensory cortex. During neural activation, brain tissue partial pressure of oxygen (P(O2)) and local cerebral blood flow (CBF) were similarly increased. To separate the P(O2) changes originating from the increase in CMRO2 and the increase in oxygen supply, the same experiments were repeated under a vasodilator-induced hypotension condition in which evoked CBF change was minimal. In this condition, evoked P(O2) monotonically decreased, indicating an increase in CMRO2. Then, CMRO2 was determined at resting as well as activation periods using a dynamic oxygen exchange model. Our results indicated that the changes in CMRO2 were linearly related with the summation of evoked field potentials and further showed that the oxygen supply in the normal condition was about 2.5 times larger than the demand. However, this oxygen oversupply was not explainable by the change in CBF alone, but at least partly by the increase in oxygenation levels at pre-capillary arterioles (e.g., 82% to 90% O2 saturation level) when local neural activity was evoked.

Dynamics of oxygen delivery and consumption during evoked neural stimulation using a compartment model and CBF and tissue P(O2) measurements

The dynamics of blood oxygen delivery and tissue consumption produced by evoked stimulation of the rat somato-sensory cortex were investigated. Tissue oxygen tension (P(O2)) and laser Doppler flowmetry (LDF) measurements were recorded under two experimental conditions: normal, which represented both oxygen delivery and consumption, and suppressed CBF (achieved using a vasodilator), which only represented tissue oxygen consumption. Forepaw stimulation for 10 s produced increases of 27.7% and 48.8% in tissue P(O2) and LDF signal under normal conditions, respectively. The tissue P(O2) response peaked 9.8 s after stimulation onset and did not show any early transient decreases indicating that measurable oxygen deficits are not required to increase the delivery of oxygen by blood flow. Under suppressed CBF conditions, the LDF signal was mostly suppressed while the tissue P(O2) decreased by 11.7% and reached a minimum 10.8 s after stimulation onset. These data were analyzed using a dynamic model that described the transport of oxygen from blood to tissue. In order to explain the differences between the model prediction of the tissue P(O2) changes and the experimental data, several hypothetical scenarios were considered, such as changes in the vascular volume, permeability-surface area or arterial oxygenation. The increase in tissue P(O2) was found to probably require the recruitment of upstream oxygen from larger arteries as well as increases in the vascular volume at the oxygen exchange sites. The amplitude of the estimated tissue tension of oxygen delivered was about 2.7 x larger than the estimated consumption under normal conditions (45.7% vs. 17.1%, respectively).

Coupling between neuronal activity and microcirculation: implications for functional brain imaging

In the neocortex, neurons with similar response properties are often clustered together in column-like structures, giving rise to what has become known as functional architecture-the mapping of various stimulus feature dimensions onto the cortical sheet. At least partially, we owe this finding to the availability of several functional brain imaging techniques, both post-mortem and in-vivo, which have become available over the last two generations, revolutionizing neuroscience by yielding information about the spatial organization of active neurons in the brain. Here, we focus on how our understanding of such functional architecture is linked to the development of those functional imaging methodologies, especially to those that image neuronal activity indirectly, through metabolic or haemodynamic signals, rather than directly through measurement of electrical activity. Some of those approaches allow exploring functional architecture at higher spatial resolution than others. In particular, optical imaging of intrinsic signals reaches the striking detail of approximately 50 mum, and, together with other methodologies, it has allowed characterizing the metabolic and haemodynamic responses induced by sensory-evoked neuronal activity. Here, we review those findings about the spatio-temporal characteristics of neurovascular coupling and discuss their implications for functional brain imaging, including position emission tomography, and non-invasive neuroimaging techniques, such as funtional magnetic resonance imaging, applicable also to the human brain.

Trial-by-trial relationship between neural activity, oxygen consumption, and blood flow responses

Trial-by-trial variability in local field potential (LFP), tissue partial pressure of oxygen (PO2), cerebral blood flow (CBF), and deoxyhemoglobin-weighted optical imaging of intrinsic signals (OIS) were tested in the rat somatosensory cortex while fixed electrical forepaw stimulation (1.0-ms pulses with amplitude of 1.2 mA at a frequency of 6 Hz) was repeatedly applied. The changes in the cerebral metabolic rate of oxygen (CMRO2) were also evaluated using a hypotension condition established by our group based on the administration of a vasodilator. Under normal conditions, CBF, PO2, and OIS showed positive signal changes (48%, 32%, and 0.42%, respectively) following stimulation. Over multiple trials, the CBF responses were well correlated with the integral of the LFP amplitudes (sigmaLFP) (Rmean=0.78), whereas a lower correlation was found between PO2 and sigmaLFP (Rmean=0.60) and between OIS and sigmaLFP (Rmean=0.54). Under the hypotension condition the LFP responses were preserved, but the CBF responses were suppressed and the PO2 and OIS changes were negative (-12% and -0.28%, respectively). In this condition, the trial-by-trial variations in PO2 and OIS were well correlated with the variability in sigmaLFPs (Rmean= -0.77 and -0.76, respectively), indicating a single trial coupling between CMRO2 changes and sigmaLFP. These findings show that CBF and CMRO2 signals are more directly correlated with neural activity compared to blood oxygen-sensitive methods such as OIS and BOLD fMRI.

Patterns of Cerebral Blood Flow and Systemic Hemodynamics During Arithmetic Processing

The present study explored modulations in cerebral blood flow and systemic hemodynamics during the execution of a mental calculation task in 41 healthy subjects. Time course and lateralization of blood flow velocities in the medial cerebral arteries of both hemispheres were assessed using functional transcranial Doppler sonography. Indices of systemic hemodynamics were obtained using continuous blood pressure recordings. Doppler sonography revealed a biphasic left dominant rise in cerebral blood flow velocities during task execution. Systemic blood pressure increased, whereas heart period, heart period variability, and baroreflex sensitivity declined. Blood pressure and heart period proved predictive of the magnitude of the cerebral blood flow response, particularly of its initial component. Various physiological mechanisms may be assumed to be involved in cardiovascular adjustment to cognitive demands. While specific contributions of the sympathetic and parasympathetic systems may account for the observed pattern of systemic hemodynamics, flow metabolism coupling, fast neurogenic vasodilation, and cerebral autoregulation may be involved in mediating cerebral blood flow modulations. Furthermore, during conditions of high cardiovascular reactivity, systemic hemodynamic changes exert a marked influence on cerebral blood perfusion.

Lessons from fMRI about mapping cortical columns

Recently developed fMRI can map small functional structures noninvasively and repeatedly without any depth limitation. However, there has been a persistent concern as to whether the high-resolution fMRI signals actually mark the sites of increased neural activity. To examine this outstanding issue, the authors used iso-orientation columns of isoflurane-anesthetized cats as a biological model and confirmed the neural correlation of fMRI iso-orientation maps by comparing them with intrinsic optical imaging maps. The results suggest that highest fMRI signals indeed indicate the sites of increased neuronal activity. Now fMRI can be used to determine plastic and/or developmental change of functional columnar structure possibly on a layer-to-layer basis. In this review, the authors focus mainly on what technical aspects should be considered when mapping functional cortical columns, including imaging techniques and experimental design.

High-resolution neurometabolic coupling in the lateral geniculate nucleus

The relationships between neural and metabolic processes in activated brain regions are central to the interpretation of noninvasive imaging. To examine this relationship, we have used a specialized sensor to measure simultaneously tissue oxygen changes and neural activity in colocalized regions of the cat's lateral geniculate nucleus (LGN). Previous work with this sensor has shown that a decrease or increase in tissue oxygen can be elicited by selective control of the location and extent of neural activation in the LGN. In the current study, to evaluate the temporal integration and homogeneity of neurometabolic coupling, we have determined the relationship between multiunit extracellular neural activity and tissue oxygen responses to visual stimuli of various durations and contrasts. Our results show that the negative but not the positive oxygen response changes in an approximately linear manner with stimulus duration. The relationship between the negative oxygen response and neural activity is relatively constant with stimulus duration. Moreover, both negative and positive oxygen responses saturate at high stimulus contrast levels. Coupling between neural activity and negative oxygen responses is well described by a power law function. These results help elucidate differences between the initial negative and subsequent positive metabolic responses and may be directly relevant to questions concerning brain mapping with functional magnetic resonance imaging.

Optical brain imaging in vivo: techniques and applications from animal to man

Optical brain imaging has seen 30 years of intense development, and has grown into a rich and diverse field. In-vivo imaging using light provides unprecedented sensitivity to functional changes through intrinsic contrast, and is rapidly exploiting the growing availability of exogenous optical contrast agents. Light can be used to image microscopic structure and function in vivo in exposed animal brain, while also allowing noninvasive imaging of hemodynamics and metabolism in a clinical setting. This work presents an overview of the wide range of approaches currently being applied to in-vivo optical brain imaging, from animal to man. Techniques include multispectral optical imaging, voltage sensitive dye imaging and speckle-flow imaging of exposed cortex, in-vivo two-photon microscopy of the living brain, and the broad range of noninvasive topography and tomography approaches to near-infrared imaging of the human brain. The basic principles of each technique are described, followed by examples of current applications to cutting-edge neuroscience research. In summary, it is shown that optical brain imaging continues to grow and evolve, embracing new technologies and advancing to address ever more complex and important neuroscience questions.

Neural interpretation of blood oxygenation level-dependent fMRI maps at submillimeter columnar resolution

Whether conventional gradient-echo (GE) blood oxygenation-level-dependent (BOLD) functional magnetic resonance imaging (fMRI) is able to map submillimeter-scale functional columns remains debatable mainly because of the spatially nonspecific large vessel contribution, poor sensitivity and reproducibility, and lack of independent evaluation. Furthermore, if the results from optical imaging of intrinsic signals are directly applicable, regions with the highest BOLD signals may indicate neurally inactive domains rather than active columns when multiple columns are activated. To examine these issues, we performed BOLD fMRI at a magnetic field of 9.4 tesla to map orientation-selective columns of isoflurane-anesthetized cats. We could not convincingly map orientation columns using conventional block-design stimulation and differential analysis method because of large fluctuations of signals. However, we successfully obtained GE BOLD iso-orientation maps with high reproducibility (r = 0.74) using temporally encoded continuous cyclic orientation stimulation with Fourier data analysis, which reduces orientation-nonselective signals such as draining artifacts and is less sensitive to signal fluctuations. We further reduced large vessel contribution using the improved spin-echo (SE) BOLD method but with overall decreased sensitivity. Both GE and SE BOLD iso-orientation maps excluding large pial vascular regions were significantly correlated to maps with a known neural interpretation, which were obtained in contrast agent-aided cerebral blood volume fMRI and total hemoglobin-based optical imaging of intrinsic signals at a hemoglobin iso-sbestic point (570 nm). These results suggest that, unlike the expectation from deoxyhemoglobin-based optical imaging studies, the highest BOLD signals are localized to the sites of increased neural activity when column-nonselective signals are suppressed.

Hemodynamic responses in cortex investigated with optical imaging methods. Implications for functional brain mapping

During the last 20 years, optical imaging methods - either alone or in combination with other recording techniques - has proven a fruitful approach to explore both the physiological and the functional aspects of activity-evoked hemodynamic responses in cortex. One of the main advantages of optical imaging consists in its high spatio-temporal resolution (in the order of few microns and milliseconds), allowing not only to unambiguously distinguish between activity patterns relating to the underlying functional architecture and those originating from the activation of medium/large blood vessels, but also to investigate the various activity-evoked hemodynamic processes at very fine detail. Here, we briefly review the principal findings obtained by optical imaging about the spatio-temporal properties of the various hemodynamic responses in cortex, i.e., changes in blood-oxygenation, blood-volume, and, to some extent, blood-flow. We also discuss the implications of those findings for non-invasive high-resolution functional brain imaging, in particular for fMRI. Finally, we underscore the importance of novel approaches for high-resolution blood-flow imaging, in the context of the need to gather information at fine spatial detail about the blood-flow response, necessary to constrain the multiple free parameters of hemodynamic response models.

Mapping iso-orientation columns by contrast agent-enhanced functional magnetic resonance imaging: reproducibility, specificity, and evaluation by optical imaging of intrinsic signal

Activation resembling ocular dominance or orientation columns has been mapped with high-resolution functional magnetic resonance imaging (fMRI). However, the neuronal interpretation of these functional maps is unclear because of the poor sensitivity of fMRI, unknown point spread function (PSF), and lack of comparison with independent techniques. Here we show that cerebral blood volume (CBV)-weighted fMRI with a blood plasma contrast agent (monocrystalline iron oxide nanoparticles), in combination with continuous temporally encoded stimulation, can map columnar neuronal activity in the cat primary visual cortex with high sensitivity, selectivity, and reproducibility. We examined hemodynamic response PSF by comparing these CBV-based signals with oxygen metabolism-based negative blood oxygenation level-dependent signals. A significant positive correlation exists between CBV- and metabolism-based iso-orientation maps, suggesting that the hemodynamic PSF is narrower than intercolumn distances. We also compared CBV-based fMRI with optical intrinsic signal (OIS) imaging, a technique that identifies sites of increased neuronal activity, to investigate neuronal correlation. Iso-orientation maps obtained by fMRI and OIS were well matched, indicating that areas of the highest orientation-selective CBV signals correspond to sites of increased neural activity. Using CBV-based fMRI, we successfully mapped orientation-selective functional architecture in the medial bank of the visual cortex, an area inaccessible to OIS imaging. Thus, we conclude that contrast agent-based fMRI, in combination with continuous temporally encoded stimulation, is a highly sensitive technique capable of mapping neural activity at the resolution of functional columns without depth limitation.

Caffeine reduces the initial dip in the visual BOLD response at 3 T

Localized changes in oxygen consumption related to increased neural activity can result in a small and transient "initial dip" of the blood oxygenation level-dependent (BOLD) signal used in functional magnetic resonance imaging (fMRI). The initial dip has been of great interest to the fMRI community because it may provide a more accurate and localized measure of neural activity than the conventional BOLD signal increase. Although potentially useful as a technique for human brain mapping, the initial dip is not always detected and has been a source of some controversy. In this study, the BOLD response to a 4-s long visual stimulus was measured with a 3-T MRI system in 5 healthy volunteers both before and immediately after a 200-mg oral caffeine dose. The caffeine dose significantly (P < 0.001) reduced or eliminated the initial dip in all subjects. These findings suggest that caffeine usage may be a key factor in the detection of the initial dip in human fMRI studies.