Spatial correlations of laminar BOLD and CBV responses to rat whisker stimulation with neuronal activity localized by Fos expression

High-resolution awake mouse fMRI at 14 Tesla

High-resolution awake mouse fMRI remains challenging despite extensive efforts to address motion-induced artifacts and stress. This study introduces an implantable radiofrequency (RF) surface coil design that minimizes image distortion caused by the air/tissue interface of mouse brains while simultaneously serving as a headpost for fixation during scanning. Furthermore, this study provides a thorough acclimation method used to accustom animals to the MRI environment minimizing motion induced artifacts. Using a 14T scanner, high-resolution fMRI enabled brain-wide functional mapping of visual and vibrissa stimulation at 100×100×200μm resolution with a 2s per frame sampling rate. Besides activated ascending visual and vibrissa pathways, robust BOLD responses were detected in the anterior cingulate cortex upon visual stimulation and spread through the ventral retrosplenial area (VRA) with vibrissa air-puff stimulation, demonstrating higher-order sensory processing in association cortices of awake mice. In particular, the rapid hemodynamic responses in VRA upon vibrissa stimulation showed a strong correlation with the hippocampus, thalamus, and prefrontal cortical areas. Cross-correlation analysis with designated VRA responses revealed early positive BOLD signals at the contralateral barrel cortex (BC) occurring 2 seconds prior to the air-puff in awake mice with repetitive stimulation, which was not detected using a randomized stimulation paradigm. This early BC activation indicated a learned anticipation through the vibrissa system and association cortices in awake mice under continuous training of repetitive air-puff stimulation. This work establishes a high-resolution awake mouse fMRI platform, enabling brain-wide functional mapping of sensory signal processing in higher association cortical areas.

Alpha-180 spin-echo based line-scanning method for high resolution laminar-specific fMRI

Laminar-specific functional magnetic resonance imaging (fMRI) has been widely used to study circuit-specific neuronal activity by mapping spatiotemporal fMRI response patterns across cortical layers. Hemodynamic responses reflect indirect neuronal activity given limit of spatial and temporal resolution. Previous gradient-echo based line-scanning fMRI (GELINE) method was proposed with high temporal (50 ms) and spatial (50 µm) resolution to better characterize the fMRI onset time across cortical layers by employing 2 saturation RF pulses. However, the imperfect RF saturation performance led to poor boundary definition of the reduced region of interest (ROI) and aliasing problems outside of the ROI. Here, we propose α (alpha)-180 spin-echo-based line-scanning fMRI (SELINE) method to resolve this issue by employing a refocusing 180° RF pulse perpendicular to the excitation slice. In contrast to GELINE signals peaked at the superficial layer, we detected varied peaks of laminar-specific BOLD signals across deeper cortical layers with the SELINE method, indicating the well-defined exclusion of the large drain-vein effect with the spin-echo sequence. Furthermore, we applied the SELINE method with 200 ms TR to sample the fast hemodynamic changes across cortical layers with a less draining vein effect. In summary, this SELINE method provides a novel acquisition scheme to identify microvascular-sensitive laminar-specific BOLD responses across cortical depth.

Lactate as a determinant of neuronal excitability, neuroenergetics and beyond

Over the last decades, lactate has emerged as important energy substrate for the brain fueling of neurons. A growing body of evidence now indicates that it is also a signaling molecule modulating neuronal excitability and activity as well as brain functions. In this review, we will briefly summarize how different cell types produce and release lactate. We will further describe different signaling mechanisms allowing lactate to fine-tune neuronal excitability and activity, and will finally discuss how these mechanisms could cooperate to modulate neuroenergetics and higher order brain functions both in physiological and pathological conditions.

Capillary responses to functional and pathological activations rely on the capillary states at rest

Brain capillaries play a crucial role in maintaining cellular viability and thus preventing neurodegeneration. The aim of this study was to characterize the brain capillary morphology at rest and during neural activation based on a big data analysis from three-dimensional microangiography. Neurovascular responses were measured using a genetic calcium sensor expressed in neurons and microangiography with two-photon microscopy, while neural acivity was modulated by stimulation of contralateral whiskers or by a seizure evoked by kainic acid. For whisker stimulation, 84% of the capillary sites showed no detectable diameter change. The remaining 10% and 6% were dilated and constricted, respectively. Significant differences were observed for capillaries in the diameter at rest between the locations of dilation and constriction. Even the seizures resulted in 44% of the capillaries having no detectable change in diameter, while 56% of the capillaries dilated. The extent of dilation was dependent on the diameter at rest. In conclusion, big data analysis on brain capillary morphology has identified at least two types of capillary states: capillaries with diameters that are relatively large at rest and stable over time regardless of neural activity and capillaries whose diameters are relatively small at rest and vary according to neural activity.

The current status and trend of the functional magnetic resonance combined with stimulation in animals

As a non-radiative, non-invasive imaging technique, functional magnetic resonance imaging (fMRI) has excellent effects on studying the activation of blood oxygen levels and functional connectivity of the brain in human and animal models. Compared with resting-state fMRI, fMRI combined with stimulation could be used to assess the activation of specific brain regions and the connectivity of specific pathways and achieve better signal capture with a clear purpose and more significant results. Various fMRI methods and specific stimulation paradigms have been proposed to investigate brain activation in a specific state, such as electrical, mechanical, visual, olfactory, and direct brain stimulation. In this review, the studies on animal brain activation using fMRI combined with different stimulation methods were retrieved. The instruments, experimental parameters, anesthesia, and animal models in different stimulation conditions were summarized. The findings would provide a reference for studies on estimating specific brain activation using fMRI combined with stimulation.

Thalamic activations in rat brain by fMRI during tactile (forepaw, whisker) and non-tactile (visual, olfactory) sensory stimulations

The thalamus is a crucial subcortical hub that impacts cortical activity. Tracing experiments in animals and post-mortem humans suggest rich morphological specificity of the thalamus. Very few studies reported rodent thalamic activations by functional MRI (fMRI) as compared to cortical activations for different sensory stimuli. Here, we show different portions of the rat thalamus in response to tactile (forepaw, whisker) and non-tactile (visual, olfactory) sensory stimuli with high field fMRI (11.7T) using a custom-build quadrature surface coil to capture high sensitivity signals from superficial and deep brain regions simultaneously. Results demonstrate reproducible thalamic activations during both tactile and non-tactile stimuli. Forepaw and whisker stimuli activated broader regions within the thalamus: ventral posterior lateral (VPL), ventral posterior medial (VPM), lateral posterior mediorostral (LPMR) and posterior medial (POm) thalamic nuclei. Visual stimuli activated dorsal lateral geniculate nucleus (DLG) of the thalamus but also parts of the superior/inferior colliculus, whereas olfactory stimuli activated specifically the mediodorsal nucleus of the thalamus (MDT). BOLD activations in LGN and MDT were much stronger than in VPL, VPM, LPMR and POm. These fMRI-based thalamic activations suggest that forepaw and whisker (i.e., tactile) stimuli engage VPL, VPM, LPMR and POm whereas visual and olfactory (i.e., non-tactile) stimuli, respectively, recruit DLG and MDT exclusively.

Improvement of sensitivity and specificity for laminar BOLD fMRI with double spin-echo EPI in humans at 7 T

Mapping mesoscopic cortical functional units such as columns or laminae is increasingly pursued by ultra-high field (UHF) functional magnetic resonance imaging (fMRI). The most popular approach for high-resolution fMRI is currently gradient-echo (GE) blood oxygenation level-dependent (BOLD) fMRI. However, its spatial accuracy is reduced due to its sensitivity to draining vessels, including pial veins, whereas spin-echo (SE) BOLD signal is expected to have higher spatial accuracy, albeit with lower sensitivity than the GE-BOLD signal. Here, we introduce a new double spin-echo (dSE) echo-planar imaging (EPI) method to improve the sensitivity of SE-BOLD contrast by averaging two spin-echoes using three radiofrequency pulses. Human fMRI experiments were performed with slices perpendicular to the central sulcus between motor and sensory cortices at 7 T during fist-clenching with touching. First, we evaluated the feasibility of single-shot dSE-EPI for BOLD fMRI with 1.5 mm isotropic resolution and found that dSE-BOLD fMRI has higher signal-to-noise ratio (SNR), temporal SNR (tSNR), and higher functional sensitivity than conventional SE-BOLD fMRI. Second, to investigate the laminar specificity of dSE-BOLD fMRI, we implemented a multi-shot approach to achieve 0.8-mm isotropic resolution with sliding-window reconstruction. Unlike GE-BOLD fMRI, the cortical profile of dSE-BOLD fMRI peaked at ~ 1.0 mm from the surface of the primary motor and sensory cortices, demonstrating an improvement of laminar specificity in humans over GE-BOLD fMRI. The proposed multi-shot dSE-EPI method is viable for high spatial resolution UHF-fMRI studies in the pursuit of resolving mesoscopic functional units.

One-pot synthesis of carboxymethyl-dextran coated iron oxide nanoparticles (CION) for preclinical fMRI and MRA applications

Superparamagnetic iron-oxide nanoparticles are robust contrast agents for magnetic resonance imaging (MRI) used for sensitive structural and functional mapping of the cerebral blood volume (CBV) when administered intravenously. To date, many CBV-MRI studies are conducted with Feraheme, manufactured for the clinical treatment of iron-deficiency. Unfortunately, Feraheme is currently not available outside the United States due to commercial and regulatory constraints, making CBV-MRI methods either inaccessible or very costly to achieve. To address this barrier, we developed a simple, one-pot recipe to synthesize Carboxymethyl-dextran coated Iron Oxide Nanoparticles, namely, “CION”, suitable for preclinical CBV-MRI applications. Here we disseminate a step-by-step instruction of our one-pot synthesis protocol, which allows CION to be produced in laboratories with minimal cost. We also characterized different CION-conjugations by manipulating polymer to metal stoichiometric ratio in terms of their size, surface chemistry, and chemical composition, and shifts in MR relaxivity and pharmacokinetics. We performed several proof-of-concept experiments in vivo, demonstrating the utility of CION for functional and structural MRI applications, including hypercapnic CO₂ challenge, visual stimulation, targeted optogenetic stimulation, and microangiography. We also present evidence that CION can serve as a cross-modality research platform by showing concurrent in vivo optical and MRI measurement of CBV using fluorescent-labeled CION. The simplicity and cost-effectiveness of our one-pot synthesis method should allow researchers to reproduce CION and tailor the relaxivity and pharmacokinetics according to their imaging needs. It is our hope that this work makes CBV-MRI more openly available and affordable for a variety of research applications.

Validating layer-specific VASO across species

Cerebral blood volume (CBV) has been shown to be a robust and important physiological parameter for quantitative interpretation of functional (f)MRI, capable of delivering highly localized mapping of neural activity. Indeed, with recent advances in ultra-high-field (≥7T) MRI hardware and associated sequence libraries, it has become possible to capture non-invasive CBV weighted fMRI signals across cortical layers. One of the most widely used approaches to achieve this (in humans) is through vascular-space-occupancy (VASO) fMRI. Unfortunately, the exact contrast mechanisms of layer-dependent VASO fMRI have not been validated for human fMRI and thus interpretation of such data is confounded. Here we validate the signal source of layer-dependent SS-SI VASO fMRI using multi-modal imaging in a rat model in response to neuronal activation (somatosensory cortex) and respiratory challenge (hypercapnia). In particular VASO derived CBV measures are directly compared to concurrent measures of total haemoglobin changes from high resolution intrinsic optical imaging spectroscopy (OIS). Quantified cortical layer profiling is demonstrated to be in agreement between VASO and contrast enhanced fMRI (using monocrystalline iron oxide nanoparticles, MION). Responses show high spatial localisation to layers of cortical processing independent of confounding large draining veins which can hamper BOLD fMRI studies, (depending on slice positioning). Thus, a cross species comparison is enabled using VASO as a common measure. We find increased VASO based CBV reactivity (3.1 ± 1.2 fold increase) in humans compared to rats. Together, our findings confirm that the VASO contrast is indeed a reliable estimate of layer-specific CBV changes. This validation study increases the neuronal interpretability of human layer-dependent VASO fMRI as an appropriate method in neuroscience application studies, in which the presence of large draining intracortical and pial veins limits neuroscientific inference with BOLD fMRI.

Optical Imaging-Based Guidance of Viral Microinjections and Insertion of a Laminar Electrophysiology Probe Into a Predetermined Barrel in Mouse Area S1BF

Wide-field Optical Imaging of Intrinsic Signals (OI-IS; Grinvald et al., 1986) is a method for imaging functional brain hemodynamic responses, mainly used to image activity from the surface of the cerebral cortex. It localizes small functional modules – such as cortical columns – with great spatial resolution and spatial specificity relative to the site of increases in neuronal activity. OI-IS is capable of imaging responses either through an intact or thinned skull or following a craniotomy. Therefore, it is minimally invasive, which makes it ideal for survival experiments. Here we describe OI-IS-based methods for guiding microinjections of optogenetics viral vectors in proximity to small functional modules (S1 barrels) of the cerebral cortex and for guiding the insertion of electrodes for electrophysiological recording into such modules. We validate our proposed methods by tissue processing of the cerebral barrel field area, revealing the track of the electrode in a predetermined barrel. In addition, we demonstrate the use of optical imaging to visualize the spatial extent of the optogenetics photostimulation, making it possible to estimate one of the two variables that conjointly determine which region of the brain is stimulated. Lastly, we demonstrate the use of OI-IS at high-magnification for imaging the upper recording contacts of a laminar probe, making it possible to estimate the insertion depth of all contacts relative to the surface of the cortex. These methods support the precise positioning of microinjections and recording electrodes, thus overcoming the variability in the spatial position of fine-scale functional modules.

An in-vivo study of BOLD laminar responses as a function of echo time and static magnetic field strength

Layer specific functional MRI requires high spatial resolution data. To compensate the associated poor signal to noise ratio it is common to integrate the signal from voxels at a given cortical depth. If the region is sufficiently large then physiological noise will be the dominant noise source. In this work, activation profiles in response to the same visual stimulus are compared at 1.5 T, 3 T and 7 T using a multi-echo, gradient echo (GE) FLASH sequence, with a 0.75 mm isotropic voxel size and the cortical integration approach. The results show that after integrating over a cortical volume of 40, 60 and 100 mm³ (at 7 T, 3 T, and 1.5 T, respectively), the signal is in the physiological noise dominated regime. The activation profiles obtained are similar for equivalent echo times. BOLD-like noise is found to be the dominant source of physiological noise. Consequently, the functional contrast to noise ratio is not strongly echo-time or field-strength dependent. We conclude that laminar GE-BOLD fMRI at lower field strengths is feasible but that larger patches of cortex will need to be examined, and that the acquisition efficiency is reduced.

Validating Linear Systems Analysis for Laminar fMRI: Temporal Additivity for Stimulus Duration Manipulations

Advancements in ultra-high field (7 T and higher) magnetic resonance imaging (MRI) scanners have made it possible to investigate both the structure and function of the human brain at a sub-millimeter scale. As neuronal feedforward and feedback information arrives in different layers, sub-millimeter functional MRI has the potential to uncover information processing between cortical micro-circuits across cortical depth, i.e. laminar fMRI. For nearly all conventional fMRI analyses, the main assumption is that the relationship between local neuronal activity and the blood oxygenation level dependent (BOLD) signal adheres to the principles of linear systems theory. For laminar fMRI, however, directional blood pooling across cortical depth stemming from the anatomy of the cortical vasculature, potentially violates these linear system assumptions, thereby complicating analysis and interpretation. Here we assess whether the temporal additivity requirement of linear systems theory holds for laminar fMRI. We measured responses elicited by viewing stimuli presented for different durations and evaluated how well the responses to shorter durations predicted those elicited by longer durations. We find that BOLD response predictions are consistently good predictors for observed responses, across all cortical depths, and in all measured visual field maps (V1, V2, and V3). Our results suggest that the temporal additivity assumption for linear systems theory holds for laminar fMRI. We thus show that the temporal additivity assumption holds across cortical depth for sub-millimeter gradient-echo BOLD fMRI in early visual cortex.

Linear systems analysis for laminar fMRI: Evaluating BOLD amplitude scaling for luminance contrast manipulations

A fundamental assumption of nearly all functional magnetic resonance imaging (fMRI) analyses is that the relationship between local neuronal activity and the blood oxygenation level dependent (BOLD) signal can be described as following linear systems theory. With the advent of ultra-high field (7T and higher) MRI scanners, it has become possible to perform sub-millimeter resolution fMRI in humans. A novel and promising application of sub-millimeter fMRI is measuring responses across cortical depth, i.e. laminar imaging. However, the cortical vasculature and associated directional blood pooling towards the pial surface strongly influence the cortical depth-dependent BOLD signal, particularly for gradient-echo BOLD. This directional pooling may potentially affect BOLD linearity across cortical depth. Here we assess whether the amplitude scaling assumption for linear systems theory holds across cortical depth. For this, we use stimuli with different luminance contrasts to elicit different BOLD response amplitudes. We find that BOLD amplitude across cortical depth scales with luminance contrast, and that this scaling is identical across cortical depth. Although nonlinearities may be present for different stimulus configurations and acquisition protocols, our results suggest that the amplitude scaling assumption for linear systems theory across cortical depth holds for luminance contrast manipulations in sub-millimeter laminar BOLD fMRI.

Delta Rhythm Orchestrates the Neural Activity Underlying the Resting State BOLD Signal via Phase-amplitude Coupling

Spontaneous ongoing neuronal activity is a prominent feature of the mammalian brain. Temporal and spatial patterns of such ongoing activity have been exploited to examine large-scale brain network organization and function. However, the neurophysiological basis of this spontaneous brain activity as detected by resting-state functional Magnetic Resonance Imaging (fMRI) remains poorly understood. To this end, multi-site local field potentials (LFP) and blood oxygenation level-dependent (BOLD) fMRI were simultaneously recorded in the rat striatum along with local pharmacological manipulation of striatal activity. Results demonstrate that delta (δ) band LFP power negatively, while beta (β) and gamma (γ) band LFPs positively correlated with BOLD fluctuation. Furthermore, there was strong cross-frequency phase-amplitude coupling (PAC), with the phase of δ LFPs significantly modulating the amplitude of the high frequency signal. Enhancing dopaminergic neuronal activity significantly reduced ventral striatal functional connectivity, δ LFP-BOLD correlation, and the PAC effect. These data suggest that different frequency bands of the LFP contribute distinctively to BOLD spontaneous fluctuation and that PAC is the organizing mechanism through which low frequency LFPs orchestrate neural activity that underlies resting state functional connectivity.

Functional connectivity with cortical depth assessed by resting state fMRI of subregions of S1 in squirrel monkeys

Whereas resting state blood oxygenation-level dependent (BOLD) functional MRI has been widely used to assess functional connectivity between cortical regions, the laminar specificity of such measures is poorly understood. This study aims to determine: (a) whether the resting state functional connectivity (rsFC) between two functionally related cortical regions varies with cortical depth, (b) the relationship between layer-resolved tactile stimulus-evoked activation pattern and interlayer rsFC pattern between two functionally distinct but related somatosensory areas 3b and 1, and (c) the effects of spatial resolution on rsFC measures. We examined the interlayer rsFC between areas 3b and 1 of squirrel monkeys under anesthesia using tactile stimulus-driven and resting state BOLD acquisitions at submillimeter resolution. Consistent with previous observations in the areas 3b and 1, we detected robust stimulus-evoked BOLD activations with foci were confined mainly to the upper layers (centered at 21% of the cortical depth). By carefully placing seeds in upper, middle, and lower layers of areas 3b and 1, we observed strong rsFC between upper and middle layers of these two areas. The layer-resolved activation patterns in areas 3b and 1 agree with their interlayer rsFC patterns, and are consistent with the known anatomical connections between layers. In summary, using BOLD rsFC pattern, we identified an interlayer interareal microcircuit that shows strong intrinsic functional connections between upper and middle layer areas 3b and 1. RsFC can be used as a robust invasive tool to probe interlayer corticocortical microcircuits.

In Vivo Visualization of Active Polysynaptic Circuits With Longitudinal Manganese-Enhanced MRI (MEMRI)

Manganese-enhanced magnetic resonance imaging (MEMRI) is a powerful tool for in vivo non-invasive whole-brain mapping of neuronal activity. Mn²⁺ enters active neurons via voltage-gated calcium channels and increases local contrast in T₁-weighted images. Given the property of Mn²⁺ of axonal transport, this technique can also be used for tract tracing after local administration of the contrast agent. However, MEMRI is still not widely employed in basic research due to the lack of a complete description of the Mn²⁺ dynamics in the brain. Here, we sought to investigate how the activity state of neurons modulates interneuronal Mn²⁺ transport. To this end, we injected mice with low dose MnCl₂ 2. (i.p., 20 mg/kg; repeatedly for 8 days) followed by two MEMRI scans at an interval of 1 week without further MnCl₂ injections. We assessed changes in T₁ contrast intensity before (scan 1) and after (scan 2) partial sensory deprivation (unilateral whisker trimming), while keeping the animals in a sensory enriched environment. After correcting for the general decay in Mn²⁺ content, whole brain analysis revealed a single cluster with higher signal in scan 1 compared to scan 2: the left barrel cortex corresponding to the right untrimmed whiskers. In the inverse contrast (scan 2 > scan 1), a number of brain structures, including many efferents of the left barrel cortex were observed. These results suggest that continuous neuronal activity elicited by ongoing sensory stimulation accelerates Mn²⁺ transport from the uptake site to its projection terminals, while the blockage of sensory-input and the resulting decrease in neuronal activity attenuates Mn²⁺ transport. The description of this critical property of Mn²⁺ dynamics in the brain allows a better understanding of MEMRI functional mechanisms, which will lead to more carefully designed experiments and clearer interpretation of the results.

Cortical depth profiles of luminance contrast responses in human V1 and V2 using 7 T fMRI

Neural activity in early visual cortex is modulated by luminance contrast. Cortical depth (i.e., laminar) contrast responses have been studied in monkey early visual cortex, but not in humans. In addition to the high spatial resolution needed and the ensuing low signal-to-noise ratio, laminar studies in humans using fMRI are hampered by the strong venous vascular weighting of the fMRI signal. In this study, we measured luminance contrast responses in human V1 and V2 with high-resolution fMRI at 7 T. To account for the effect of intracortical ascending veins, we applied a novel spatial deconvolution model to the fMRI depth profiles. Before spatial deconvolution, the contrast response in V1 showed a slight local maximum at mid cortical depth, whereas V2 exhibited a monotonic signal increase toward the cortical surface. After applying the deconvolution, both V1 and V2 showed a pronounced local maximum at mid cortical depth, with an additional peak in deep grey matter, especially in V1. Moreover, we found a difference in contrast sensitivity between V1 and V2, but no evidence for variations in contrast sensitivity as a function of cortical depth. These findings are in agreement with results obtained in nonhuman primates, but further research will be needed to validate the spatial deconvolution approach.

Autobiography of James S. Hyde

The papers, book chapters, reviews, and patents by James S. Hyde in the bibliography of this document have been separated into EPR and MRI sections, and within each section by topics. Within each topic, publications are listed chronologically. A brief summary is provided for each patent listed. A few publications and patents that do not fit this schema have been omitted. This list of publications is preceded by a scientific autobiography that focuses on selected topics that are judged to have been of most scientific importance. References to many of the publications and patents in the bibliography are made in the autobiography.

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

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

Neuronal networks and mediators of cortical neurovascular coupling responses in normal and altered brain states

Brain imaging techniques that use vascular signals to map changes in neuronal activity, such as blood oxygenation level-dependent functional magnetic resonance imaging, rely on the spatial and temporal coupling between changes in neurophysiology and haemodynamics, known as 'neurovascular coupling (NVC)'. Accordingly, NVC responses, mapped by changes in brain haemodynamics, have been validated for different stimuli under physiological conditions. In the cerebral cortex, the networks of excitatory pyramidal cells and inhibitory interneurons generating the changes in neural activity and the key mediators that signal to the vascular unit have been identified for some incoming afferent pathways. The neural circuits recruited by whisker glutamatergic-, basal forebrain cholinergic- or locus coeruleus noradrenergic pathway stimulation were found to be highly specific and discriminative, particularly when comparing the two modulatory systems to the sensory response. However, it is largely unknown whether or not NVC is still reliable when brain states are altered or in disease conditions. This lack of knowledge is surprising since brain imaging is broadly used in humans and, ultimately, in conditions that deviate from baseline brain function. Using the whisker-to-barrel pathway as a model of NVC, we can interrogate the reliability of NVC under enhanced cholinergic or noradrenergic modulation of cortical circuits that alters brain states.This article is part of the themed issue 'Interpreting BOLD: a dialogue between cognitive and cellular neuroscience'.

Low- but Not High-Frequency LFP Correlates with Spontaneous BOLD Fluctuations in Rat Whisker Barrel Cortex

Resting-state magnetic resonance imaging (rsMRI) is thought to reflect ongoing spontaneous brain activity. However, the precise neurophysiological basis of rsMRI signal remains elusive. Converging evidence supports the notion that local field potential (LFP) signal in the high-frequency range correlates with fMRI response evoked by a task (e.g., visual stimulation). It remains uncertain whether this relationship extends to rsMRI. In this study, we systematically modulated LFP signal in the whisker barrel cortex (WBC) by unilateral deflection of rat whiskers. Results show that functional connectivity between bilateral WBC was significantly modulated at the 2 Hz, but not at the 4 or 6 Hz, stimulus condition. Electrophysiologically, only in the low-frequency range (<5 Hz) was the LFP power synchrony in bilateral WBC significantly modulated at 2 Hz, but not at 4- or 6-Hz whisker stimulation, thus distinguishing these 2 experimental conditions, and paralleling the findings in rsMRI. LFP power synchrony in other frequency ranges was modulated in a way that was neither unique to the specific stimulus conditions nor parallel to the fMRI results. Our results support the hypothesis that emphasizes the role of low-frequency LFP signal underlying rsMRI.

Prefrontal dopamine regulates fear reinstatement through the downregulation of extinction circuits

Prevention of relapses is a major challenge in treating anxiety disorders. Fear reinstatement can cause relapse in spite of successful fear reduction through extinction-based exposure therapy. By utilising a contextual fear-conditioning task in mice, we found that reinstatement was accompanied by decreased c-Fos expression in the infralimbic cortex (IL) with reduction of synaptic input and enhanced c-Fos expression in the medial subdivision of the central nucleus of the amygdala (CeM). Moreover, we found that IL dopamine plays a key role in reinstatement. A reinstatement-inducing reminder shock induced c-Fos expression in the IL-projecting dopaminergic neurons in the ventral tegmental area, and the blocking of IL D1 signalling prevented reduction of synaptic input, CeM c-Fos expression, and fear reinstatement. These findings demonstrate that a dopamine-dependent inactivation of extinction circuits underlies fear reinstatement and may explain the comorbidity of substance use disorders and anxiety disorders.

Reduced microvascular volume and hemispherically deficient vasoreactivity to hypercapnia in acute ischemia: MRI study using permanent middle cerebral artery occlusion rat model

Vasoreactivity to hypercapnia has been used for assessing cerebrovascular tone and control altered by ischemic stroke. Despite the high prognostic potential, traits of hypercapnia-induced hemodynamic changes have not been fully characterized in relation with baseline vascular states and brain tissue damage. To monitor cerebrovascular responses, T2- and T2*-weighted magnetic resonance imaging (MRI) images were acquired alternatively using spin- and gradient-echo echo plannar imaging (GESE EPI) sequence with 5% CO2 gas inhalation in normal (n=5) and acute stroke rats (n=10). Dynamic relative changes in cerebrovascular volume (CBV), microvascular volume (MVV), and vascular size index (VSI) were assessed from regions of interest (ROIs) delineated by the percent decrease of apparent diffusion coefficient (ADC). The baseline CBV was not affected by middle cerebral artery occlusion (MCAO) whereas the baseline MVV in ischemic areas was significantly lower than that in the rest of the brain and correlated with ADC. Vasoreactivity to hypercapnic challenge was considerably attenuated in the entire ipsilesional hemisphere including normal ADC regions, in which unsolicited, spreading depression-associated increases of CBV and MVV were observed. The lesion-dependent inhomogeneity in baseline MVV indicates the effective perfusion reserve for accurately delineating the true ischemic damage while the cascade of neuronal depolarization is probably responsible for the hemispherically lateralized changes in overall neurovascular physiology.

Functional connectivity in rat brain at 200 μm resolution

The somatosensory functional magnetic resonance imaging (fMRI) response to electrical stimulation of the middle phalange of the second digit of four rats at a spatial resolution of 200 μm cubic at 9.4 T is reported. At high threshold (p<0.002), activated voxels encompass a penetrating vein that passes across the gray matter. These voxels lie mostly in three contiguous slices perpendicular to the pial surface. This activation is assigned to the representation in the forepaw barrel subfield (FBS) of a single cortical column of this phalange. In addition, the activation of the indusium griseum (IG) is visualized robustly. Voxels revealed by fMRI were used to observe functional connectivity to other voxels of the sensorimotor cortex using fcMRI. Results of this experiment were analyzed as a function of decreasing threshold, which exhibited spreading connectivity that revealed S2, M1/M2, and contralateral S1. Noting that every cubic millimeter of tissue contains 125 voxels, connectivity patterns are complex. It is hypothesized that they reflect connections within gray matter by association fibers. S2 and IG revealed connectivities with many voxels across the sensorimotor cortex. These regions also showed subregional variation of connectivity. A 1-cm-diameter surface coil with a local low-noise RF amplifier was used in these studies. The usual region of sensitivity (ROS) of such a coil is 1 cm diameter by 0.5 cm depth. Significant connectivity was observed between time courses of voxels that were within the ROS and voxels that were outside, which extends the volume of tissue that can be observed by the methods of this article.

Quantitative activity-induced manganese-dependent MRI for characterizing cortical layers in the primary somatosensory cortex of the rat

The ability of Mn(2+) to follow Ca(2+) pathways upon stimulation transform them into remarkable surrogate markers of neuronal activity using activity-induced manganese-dependent MRI (AIM-MRI). In the present study, a precise follow-up of physiological parameters during MnCl2 and mannitol infusions improved the reproducibility of AIM-MRI allowing in-depth evaluation of the technique. Pixel-by-pixel T1 data were investigated using histogram distributions in the barrel cortex (BC) and the thalamus before and after Mn(2+) infusion, after blood brain barrier opening and after BC activation. Mean BC T1 values dropped significantly upon trigeminal nerve (TGN) stimulation (-38 %, P = 0.02) in accordance with previous literature findings. T1 histogram distributions showed that 34 % of T1s in the range 600-1500 ms after Mn(2+ )+ mannitol infusions shifted to 50-350 ms after TGN stimulation corresponding to a twofold increase of the percentage of pixels with the lowest T1s in BC. Moreover, T1 changes in response to stimulation increased significantly from superficial cortical layers (I-III) to deeper layers (V-VI). Cortical cytoarchitecture detection during a functional paradigm was performed extending the potential of AIM-MRI. Quantitative AIM-MRI could thus offer a means to interpret local neural activity across cortical layers while identification of the role of calcium dynamics in vivo during brain activation could play a key role in resolving neurovascular coupling mechanisms.

Optical Coherence Tomography angiography reveals laminar microvascular hemodynamics in the rat somatosensory cortex during activation

The BOLD (blood-oxygen-level dependent) fMRI (functional Magnetic Resonance Imaging) signal is shaped, in part, by changes in red blood cell (RBC) content and flow across vascular compartments over time. These complex dynamics have been challenging to characterize directly due to a lack of appropriate imaging modalities. In this study, making use of infrared light scattering from RBCs, depth-resolved Optical Coherence Tomography (OCT) angiography was applied to image laminar functional hyperemia in the rat somatosensory cortex. After defining and validating depth-specific metrics for changes in RBC content and speed, laminar hemodynamic responses in microvasculature up to cortical depths of >1mm were measured during a forepaw stimulus. The results provide a comprehensive picture of when and where changes in RBC content and speed occur during and immediately following cortical activation. In summary, the earliest and largest microvascular RBC content changes occurred in the middle cortical layers, while post-stimulus undershoots were most prominent superficially. These laminar variations in positive and negative responses paralleled known distributions of excitatory and inhibitory synapses, suggesting neuronal underpinnings. Additionally, the RBC speed response consistently returned to baseline more promptly than RBC content after the stimulus across cortical layers, supporting a "flow-volume mismatch" of hemodynamic origin.

Direct 3D Analyses Reveal Barrel-Specific Vascular Distribution and Cross-Barrel Branching in the Mouse Barrel Cortex

Whether vascular distribution is spatially specific among cortical columns is a fundamental yet controversial question. Here, we have obtained 1-μm resolution 3D datasets that cover the whole mouse barrel cortex by combining Nissl staining with micro-optical sectioning tomography to simultaneously visualize individual cells and blood vessels, including capillaries. Pinpointing layer IV of the posteromedial barrel subfield, direct 3D reconstruction and quantitative analysis showed that (1) penetrating vessels preferentially locate in the interbarrel septa/barrel wall (75.1%) rather than the barrel hollows, (2) the branches of 70% penetrating vessels only reach the neighboring but not always all the neighboring barrels and the other 30% extend beyond the neighboring barrels and may provide cross-barrel blood supply or drainage, (3) the branches of 59.6% penetrating vessels reach all the neighboring barrels, while the rest only reach part of them, and (4) the length density of microvessels in the interbarrel septa/barrel wall is lower than that in the barrel hollows with a ratio of 0.92. These results reveal that the penetrating vessels and microvessels exhibit a barrel-specific organization, whereas the branches of penetrating vessels do not, which suggests a much more complex vascular distribution pattern among cortical columns than previously thought.

Repeated BOLD-fMRI imaging of deep brain stimulation responses in rats

Functional magnetic resonance imaging (fMRI) provides a picture of the global spatial activation pattern of the brain. Interest is growing regarding the application of fMRI to rodent models to investigate adult brain plasticity. To date, most rodent studies used an electrical forepaw stimulation model to acquire fMRI data, with α-chloralose as the anesthetic. However, α-chloralose is harmful to animals, and not suitable for longitudinal studies. Moreover, peripheral stimulation models enable only a limited number of brain regions to be studied. Processing between peripheral regions and the brain is multisynaptic, and renders interpretation difficult and uncertain. In the present study, we combined the medetomidine-based fMRI protocol (a noninvasive rodent fMRI protocol) with chronic implantation of an MRI-compatible stimulation electrode in the ventroposterior (VP) thalamus to repetitively sample thalamocortical responses in the rat brain. Using this model, we scanned the forebrain responses evoked by the VP stimulation repeatedly of individual rats over 1 week. Cortical BOLD responses were compared between the 2 profiles obtained at day1 and day8. We discovered reproducible frequency- and amplitude-dependent BOLD responses in the ipsilateral somatosensory cortex (S1). The S1 BOLD responses during the 2 sessions were conserved in maximal response amplitude, area size (size ratio from 0.88 to 0.91), and location (overlap ratio from 0.61 to 0.67). The present study provides a long-term chronic brain stimulation protocol for studying the plasticity of specific neural circuits in the rodent brain by BOLD-fMRI.

FMRI of deep brain stimulation at the rat ventral posteromedial thalamus

BACKGROUND

Functional magnetic resonance imaging (fMRI) of deep brain stimulation (DBS) has potentials to reveal neuroanatomical connectivity of a specific brain region in vivo.

OBJECTIVE

This study aimed to demonstrate frequency and amplitude tunings of the thalamocortical tract using DBS fMRI at the rat ventral posteromedial thalamus.

METHODS

Blood oxygenation level dependent (BOLD) fMRI data were acquired in a total of twelve rats at a high-field 11.7 T MRI scanner with modulation of nine stimulus frequencies (1-40 Hz) and seven stimulus amplitudes (0.2-3.6 mA).

RESULTS

BOLD response in the barrel cortex peaked at 25 Hz. The response increased with stimulus amplitude and reached a plateau at 1 mA. Cortical spreading depolarization (CSD) was observed occasionally after DBS that carries >10% BOLD waves spanning the entire ipsilateral cortex.

CONCLUSION

fMRI is sensitive to the frequency effect of DBS and has potential to investigate the function of a particular neuroanatomical pathway.

Resting state functional connectivity: its physiological basis and application in neuropharmacology

Brain structures do not work in isolation; they work in concert to produce sensory perception, motivation and behavior. Systems-level network activity can be investigated by resting state magnetic resonance imaging (rsMRI), an emerging neuroimaging technique that assesses the synchrony of the brain's ongoing spontaneous activity. Converging evidence reveals that rsMRI is able to consistently identify distinct spatiotemporal patterns of large-scale brain networks. Dysregulation within and between these networks has been implicated in a number of neurodegenerative and neuropsychiatric disorders, including Alzheimer's disease and drug addiction. Despite wide application of this approach in systems neuroscience, the physiological basis of these fluctuations remains incompletely understood. Here we review physiological studies in electrical, metabolic and hemodynamic fluctuations that are most pertinent to the rsMRI signal. We also review recent applications to neuropharmacology - specifically drug effects on resting state fluctuations. We speculate that the mechanisms governing spontaneous fluctuations in regional oxygenation availability likely give rise to the observed rsMRI signal. We conclude by identifying several open questions surrounding this technique. This article is part of the Special Issue Section entitled 'Neuroimaging in Neuropharmacology'.

Cerebral blood volume MRI with intravascular superparamagnetic iron oxide nanoparticles

The cerebral blood volume (CBV) is a crucial physiological indicator of tissue viability and vascular reactivity. Thus, noninvasive CBV mapping has been of great interest. For this, ultrasmall superparamagnetic iron oxide (USPIO) nanoparticles, including monocrystalline iron oxide nanoparticles, can be used as long-half-life, intravascular susceptibility agents of CBV MRI measurements. Moreover, CBV-weighted functional MRI (fMRI) with USPIO nanoparticles provides enhanced sensitivity, reduced large vessel contribution and improved spatial specificity relative to conventional blood oxygenation level-dependent fMRI, and measures a single physiological parameter that is easily interpretable. We review the physiochemical and magnetic properties, and pharmacokinetics, of USPIO nanoparticles in brief. We then extensively discuss quantifications of baseline CBV, vessel size index and functional CBV change. We also provide reviews of dose-dependent sensitivity, vascular filter function, specificity, characteristics and impulse response function of CBV fMRI. Examples of CBV fMRI specificity at the laminar and columnar resolution are provided. Finally, we briefly review the application of CBV measurements to functional and pharmacological studies in animals. Overall, the use of USPIO nanoparticles can determine baseline CBV and its changes induced by functional activity and pharmacological interventions.

Ultra high-resolution fMRI and electrophysiology of the rat primary somatosensory cortex

High-resolution functional-magnetic-resonance-imaging (fMRI) has been used to study brain functions at increasingly finer scale, but whether fMRI can accurately reflect layer-specific neuronal activities is less well understood. The present study investigated layer-specific cerebral-blood-volume (CBV) fMRI and electrophysiological responses in the rat cortex. CBV fMRI at 40×40 μm in-plane resolution was performed on an 11.7-T scanner. Electrophysiology used a 32-channel electrode array that spanned the entire cortical depth. Graded electrical stimulation was used to study activations in different cortical layers, exploiting the notion that most of the sensory-specific neurons are in layers II-V and most of the nociceptive-specific neurons are in layers V-VI. CBV response was strongest in layer IV of all stimulus amplitudes. Current source density analysis showed strong sink currents at cortical layers IV and VI. Multi-unit activities mainly appeared at layers IV-VI and peaked at layer V. Although our measures showed scaled activation profiles during modulation of stimulus amplitude and failed to detect specific recruitment at layers V and VI during noxious electrical stimuli, there appears to be discordance between CBV fMRI and electrophysiological peak responses, suggesting neurovascular uncoupling at laminar resolution. The technique implemented in the present study offers a means to investigate intracortical neurovascular function in the normal and diseased animal models at laminar resolution.

Cerebral microvascular network geometry changes in response to functional stimulation

The cortical microvessels are organized in an intricate, hierarchical, three-dimensional network. Superimposed on this anatomical complexity is the highly complicated signaling that drives the focal blood flow adjustments following a rise in the activity of surrounding neurons. The microvascular response to neuronal activation remains incompletely understood. We developed a custom two photon fluorescence microscopy acquisition and analysis to obtain 3D maps of neuronal activation-induced changes in the geometry of the microvascular network of the primary somatosensory cortex of anesthetized rats. An automated, model-based tracking algorithm was employed to reconstruct the 3D microvascular topology and represent it as a graph. The changes in the geometry of this network were then tracked, over time, in the course of electrical stimulation of the contralateral forepaw. Both dilatory and constrictory responses were observed across the network. Early dilatory and late constrictory responses propagated from deeper to more superficial cortical layers while the response of the vertices that showed initial constriction followed by later dilation spread from cortical surface toward increasing cortical depths. Overall, larger caliber adjustments were observed deeper inside the cortex. This work yields the first characterization of the spatiotemporal pattern of geometric changes on the level of the cortical microvascular network as a whole and provides the basis for bottom-up modeling of the hemodynamically-weighted neuroimaging signals.

Isolation rearing significantly perturbs brain metabolism in the thalamus and hippocampus

Psychosocial neglect during childhood severely impairs both behavioral and physical health. The isolation rearing model in rodents has been employed by our group and others to study this clinical problem at a basic level. We previously showed that immediate early gene (IEG) expression in the hippocampus and medial prefrontal cortex (mPFC) is decreased in isolation-reared (IR) compared to group-reared (GR) rats. In the current study, we sought to evaluate: (1) whether these changes in IEG expression would be detected by the measurement of brain glucose metabolism using positron emission tomography (PET) with fluorodeoxyglucose (FDG) and (2) whether PET FDG could illuminate other brain regions with different glucose metabolism in IR compared to GR rats. We found that there were significant differences in FDG uptake in the hippocampus that were consistent with our findings for IEG expression (decreased mean FDG uptake in IR rats). In contrast, in the mPFC, the FDG uptake between IR and GR rats did not differ. Finally, we found decreased mean FDG uptake in the thalamus of the IR rats, a region we had not previously examined. The results suggest that PET FDG has the potential to be utilized as a biomarker of molecular changes in the hippocampus. Further, the differences found in thalamic brain FDG uptake suggest that further investigation of this region at the molecular and cellular levels may provide an important insight into the neurobiological basis of the adverse clinical outcomes found in children exposed to psychosocial deprivation.

Contrast-enhanced functional blood volume imaging (CE-fBVI): enhanced sensitivity for brain activation in humans using the ultrasmall superparamagnetic iron oxide agent ferumoxytol

Functional MRI (fMRI) brain studies performed in the presence of a steady-state or "blood pool" contrast agent yields activation maps that are weighted for cerebral blood volume (CBV). Previous animal experiments suggest significant contrast-to-noise ratio (CNR) improvements, but these studies have not yet been performed in humans due to the lack of availability of a suitable agent. Here we report the use of the USPIO ferumoxytol (AMAG Pharmaceuticals, Inc., Cambridge, MA) for functional brain activation in humans, termed contrast enhanced functional blood volume imaging (CE-fBVI). Four subjects were scanned during a unilateral finger tapping task with standard blood-oxygen level dependent (BOLD) imaging before contrast and CE-fBVI after contrast injection. The CE-fBVI response showed both a fast (5.8±1.3 s) and a slow (75.3±27.5 s) component of CBV response to stimuli. A significant CNR gain of approximately 2-3 was found for CE-fBVI compared to BOLD fMRI. Interestingly, less susceptibility-related signal dropouts were observed in the inferior frontal and temporal lobes with CE-fBVI. The combination of higher CNR and better spatial specificity, enabled by CE-fBVI using blood pool USPIO contrast agent opens the door to higher resolution brain mapping.

The history and role of long duration stimulation in fMRI

During the past 20 years, BOLD fMRI has developed towards a central and fundamental tool in neuroscience. It has been shown that the BOLD response provides an indicator of neuronal activity in the brain. Consequently, for an accurate interpretation of findings in BOLD MRI experiments and to draw meaningful conclusions about the temporal evolution of neural events, a deep understanding of the nature of the BOLD contrast has become of essential importance. Since the dynamics of the major direct determinants of the BOLD signal (CBF, CBV and CMRO(2)) range between seconds and minutes, long duration stimulation was an early key strategy needed to study and understand the BOLD characteristics. This paper summarizes and discusses the thoughts and rationales of the long duration stimulation studies.

Neuronal inhibition and excitation, and the dichotomic control of brain hemodynamic and oxygen responses

Brain's electrical activity correlates strongly to changes in cerebral blood flow (CBF) and the cerebral metabolic rate of oxygen (CMRO(2)). Subthreshold synaptic processes correlate better than the spike rates of principal neurons to CBF, CMRO(2) and positive BOLD signals. Stimulation-induced rises in CMRO(2) are controlled by the ATP turnover, which depends on the energy used to fuel the Na,K-ATPase to reestablish ionic gradients, while stimulation-induced CBF responses to a large extent are controlled by mechanisms that depend on Ca(2+) rises in neurons and astrocytes. This dichotomy of metabolic and vascular control explains the gap between the stimulation-induced rises in CMRO(2) and CBF, and in turn the BOLD signal. Activity-dependent rises in CBF and CMRO(2) vary within and between brain regions due to differences in ATP turnover and Ca(2+)-dependent mechanisms. Nerve cells produce and release vasodilators that evoke positive BOLD signals, while the mechanisms that control negative BOLD signals by activity-dependent vasoconstriction are less well understood. Activation of both excitatory and inhibitory neurons produces rises in CBF and positive BOLD signals, while negative BOLD signals under most conditions correlate to excitation of inhibitory interneurons, but there are important exceptions to that rule as described in this paper. Thus, variations in the balance between synaptic excitation and inhibition contribute dynamically to the control of metabolic and hemodynamic responses, and in turn the amplitude and polarity of the BOLD signal. Therefore, it is not possible based on a negative or positive BOLD signal alone to decide whether the underlying activity goes on in principal or inhibitory neurons.

Acute nicotine-induced tachyphylaxis is differentially manifest in the limbic system

Rapid tolerance develops to many of nicotine's behavioral and autonomic effects. A better understanding of the spatiotemporal patterns in neuronal activity as a consequence of acute nicotine tolerance (tachyphylaxis) may help explain its commonly found inverted 'U'-shaped biphasic dose-effect relationship on various behaviors. To this end, we employed high-resolution functional magnetic resonance imaging and relative cerebral blood volume (rCBV) as a marker of neuronal activity, to characterize the regional development of acute tolerance as a function of nicotine dose in naïve, anesthetized rats. A single intravenous nicotine injection at 0.1 and 0.3, but not 0.03 mg/kg, significantly increased neuronal activity in many neocortical areas. In contrast, dose-dependent increases in rCBV were most pronounced in limbic regions, such that responses seen at 0.1 mg/kg nicotine in accumbens, hippocampus, amygdala, and several other limbic areas were not seen following 0.3 mg/kg nicotine. Finally, whereas profound tolerance was observed in many cortical regions after the second of two paired nicotine injections at either 0.1 or 0.3 mg/kg, subcortical limbic structures showed only a weak trend for tolerance. Lack of rCBV changes in animals receiving nicotine methiodide, a quaternary nicotine analog that does not cross the blood-brain barrier, supports a direct neuronal effect of nicotine rather than an action on the vasculature. These data provide pharmacodynamic insight into the regional heterogeneity of nicotine tachyphylaxis development, which may be relevant to behavioral and neurobiological mechanisms associated with repeated tobacco consumption.

Pyramidal neurons are "neurogenic hubs" in the neurovascular coupling response to whisker stimulation

The whisker-to-barrel cortex is widely used to study neurovascular coupling, but the cellular basis that underlies the perfusion changes is still largely unknown. Here, we identified neurons recruited by whisker stimulation in the rat somatosensory cortex using double immunohistochemistry for c-Fos and markers of glutamatergic and GABAergic neurons, and investigated in vivo their contribution along with that of astrocytes in the evoked perfusion response. Whisker stimulation elicited cerebral blood flow (CBF) increases concomitantly with c-Fos upregulation in pyramidal cells that coexpressed cyclooxygenase-2 (COX-2) and GABA interneurons that coexpressed vasoactive intestinal polypeptide and/or choline acetyltransferase, but not somatostatin or parvalbumin. The evoked CBF response was decreased by blockade of NMDA (MK-801, -37%), group I metabotropic glutamate (MPEP+LY367385, -40%), and GABA-A (picrotoxin, -31%) receptors, but not by GABA-B, VIP, or muscarinic receptor antagonism. Picrotoxin decreased stimulus-induced somatosensory evoked potentials and CBF responses. Combined blockade of GABA-A and NMDA receptors yielded an additive decreasing effect (-61%) of the evoked CBF compared with each antagonist alone, demonstrating cooperation of both excitatory and inhibitory systems in the hyperemic response. Blockade of prostanoid synthesis by inhibiting COX-2 (indomethacin, NS-398), expressed by ∼40% of pyramidal cells but not by astrocytes, impaired the CBF response (-50%). The hyperemic response was also reduced (-40%) after inhibition of astroglial oxidative metabolism or epoxyeicosatrienoic acids synthesis. These results demonstrate that changes in pyramidal cell activity, sculpted by specific types of inhibitory GABA interneurons, drive the CBF response to whisker stimulation and, further, that metabolically active astrocytes are also required.

Frequency tuning in the rat whisker barrel cortex revealed through RBC flux maps

The rodent whisker barrel cortex is ideal for studies related to sensory processing and neural plasticity in the brain. However, its small spatial dimensions challenge optical and other imaging technologies mapping cortical hemodynamics as functional resolution (the ability to spatially and selectively discriminate signals from microvascular compartments) limit measurement accuracy. To precisely map hemodynamic activity within the rat posteriomedial barrel subfield (PMBSF), we used functional Laser Doppler Imaging (fLDI) at high spatial resolution with optimized detection and analysis. In this configuration, we demonstrate prominent whisker deflection-induced fLDI hemodynamic responses from microvascular regions indicating the technique's specificity to smaller vessel compartments. Clusters of fLDI activation were confined within the PMBSF region during deflection of either single or all whiskers. Stereotaxic co-ordinate mapping was performed over all animals leading to an average maximum activity cluster at +5.3, -3.5 from the Bregma. The maximum activity cluster during all whisker stimulation combined with the principal activation cluster during deflection of the C1 whisker were used as a reference to characterize the fLDI maps within the PMBSF. fLDI activation area increased with the frequency of whisker deflection. In a quantitative analysis, we reveal the increase in the spatial extent of fLDI activation with stimulation frequency as spatially non-uniform with a bias towards the caudal region for low and rostral region for higher stimulation frequencies.

Direct imaging of macrovascular and microvascular contributions to BOLD fMRI in layers IV-V of the rat whisker-barrel cortex

The spatiotemporal characteristics of the hemodynamic response to increased neural activity were investigated at the level of individual intracortical vessels using BOLD-fMRI in a well-established rodent model of somatosensory stimulation at 11.7 T. Functional maps of the rat barrel cortex were obtained at 150 × 150 × 500 μm spatial resolution every 200 ms. The high spatial resolution allowed separation of active voxels into those containing intracortical macro vessels, mainly vein/venules (referred to as macrovasculature), and those enriched with arteries/capillaries and small venules (referred to as microvasculature) since the macro vessel can be readily mapped due to the fast T2 decay of blood at 11.7 T. The earliest BOLD response was observed within layers IV-V by 0.8s following stimulation and encompassed mainly the voxels containing the microvasculature and some confined macrovasculature voxels. By 1.2s, the BOLD signal propagated to the macrovasculature voxels where the peak BOLD signal was 2-3 times higher than that of the microvasculature voxels. The BOLD response propagated in individual venules/veins far from neuronal sources at later times. This was also observed in layers IV-V of the barrel cortex after specific stimulation of separated whisker rows. These results directly visualized that the earliest hemodynamic changes to increased neural activity occur mainly in the microvasculature and spread toward the macrovasculature. However, at peak response, the BOLD signal is dominated by penetrating venules even at layers IV-V of the cortex.

MRI of neuronal plasticity in rodent models

Modifications in the behavior and architecture of neuronal networks are well documented to occur in association with learning and memory, as well as following injury. These plasticity mechanisms are crucial to ensure adequate processing of stimuli, and they also dictate the degree of recovery following peripheral or central nervous system injury. Nevertheless, the underlying neuronal mechanisms that determine the degree of plasticity of neuronal pathways are not fully understood. Recent developments in animal-dedicated magnetic resonance imaging (MRI) scanners and related hardware afford a high spatial and temporal resolution, making functional MRI and manganese-enhanced MRI emerging tools for studying reorganization of neuronal pathways in rodent models. Many of the observed changes in neuronal functions in rodent's brains following injury discussed here agree with clinical human fMRI findings. This demonstrates that animal model imaging can have a significant clinical impact in the neuronal plasticity and rehabilitation arenas.

Temporal dynamics and spatial specificity of arterial and venous blood volume changes during visual stimulation: implication for BOLD quantification

Determination of compartment-specific cerebral blood volume (CBV) changes is important for understanding neurovascular physiology and quantifying blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI). In isoflurane-anesthetized cats, we measured the spatiotemporal responses of arterial CBV (CBV(a)) and total CBV (CBV(t)) induced by a 40-second visual stimulation, using magnetization transfer (MT)-varied BOLD and contrast-agent fMRI techniques at 9.4 T. To determine the venous CBV (CBV(v)) change, we calculated the difference between CBV(t) and CBV(a) changes. The dynamic response of CBV(a) was an order of magnitude faster than that of CBV(v), while the magnitude of change under steady-state conditions was similar between the two. Following stimulation offset, ΔCBV(a) showed small poststimulus undershoots, while ΔCBV(v) slowly returned to baseline. The largest CBV(a) and CBV(t) response occurred after 10 seconds of simulation in cortical layer 4, which we identified as the stripe of Gennari by T(1)-weighted MRI. The CBV(v) response, however, was not specific across the cortical layers during the entire stimulation period. Our data indicate that rapid, more-specific arterial vasodilation is followed by slow, less-specific venous dilation. Our finding implies that the contribution of CBV(v) changes to BOLD signals is significant for long, but not short, stimulation periods.