Simultaneous blood oxygenation level-dependent and cerebral blood flow functional magnetic resonance imaging during forepaw stimulation in the rat

A systematic review of the relationship between magnetic resonance imaging based resting-state and structural networks in the rodent brain

Recent developments in rodent brain imaging have enabled translational characterization of functional and structural connectivity at the whole brain level in vivo. Nevertheless, fundamental questions about the link between structural and functional networks remain unsolved. In this review, we systematically searched for experimental studies in rodents investigating both structural and functional network measures, including studies correlating functional connectivity using resting-state functional MRI with diffusion tensor imaging or viral tracing data. We aimed to answer whether functional networks reflect the architecture of the structural connectome, how this reciprocal relationship changes throughout a disease, how structural and functional changes relate to each other, and whether changes follow the same timeline. We present the knowledge derived exclusively from studies that included in vivo imaging of functional and structural networks. The limited number of available reports makes it difficult to draw general conclusions besides finding a spatial and temporal decoupling between structural and functional networks during brain disease. Data suggest that when overcoming the currently limited evidence through future studies with combined imaging in various disease models, it will be possible to explore the interaction between both network systems as a disease or recovery biomarker.

Remote ischemic postconditioning increased cerebral blood flow and oxygenation assessed by magnetic resonance imaging in newborn piglets after hypoxia-ischemia

BACKGROUND

We have previously investigated neurological outcomes following remote ischemic postconditioning (RIPC) in a newborn piglet model of hypoxic-ischemic encephalopathy. The aim of this study was to further investigate potential mechanisms of neuroprotection by comparing newborn piglets subjected to global hypoxia-ischemia (HI) treated with and without RIPC with regards to measures of cerebral blood flow and oxygenation assessed by functional magnetic resonance imaging.

MATERIALS AND METHODS

A total of 50 piglets were subjected to 45 min global HI and randomized to either no treatment or RIPC treatment. Magnetic resonance imaging was performed 72 h after the HI insult with perfusion-weighted (arterial spin labeling, ASL) and oxygenation-weighted (blood-oxygen-level-dependent, BOLD) sequences in the whole brain, basal ganglia, thalamus, and cortex. Four sham animals received anesthesia and mechanical ventilation only.

RESULTS

Piglets treated with RIPC had higher measures of cerebral blood flow in all regions of interest and the whole brain (mean difference: 2.6 ml/100 g/min, 95% CI: 0.1; 5.2) compared with the untreated controls. They also had higher BOLD values in the basal ganglia and the whole brain (mean difference: 4.2 T2*, 95% CI: 0.4; 7.9). Measures were similar between piglets treated with RIPC and sham animals.

CONCLUSION

Piglets treated with RIPC had higher measures of cerebral blood flow and oxygenation assessed by magnetic resonance imaging in the whole brain and several regions of interest compared with untreated controls 72 h after the HI insult. Whether this reflects a potential neuroprotective mechanism of RIPC requires further study.

Contribution of animal models toward understanding resting state functional connectivity

Functional connectivity, which reflects the spatial and temporal organization of intrinsic activity throughout the brain, is one of the most studied measures in human neuroimaging research. The noninvasive acquisition of resting state functional magnetic resonance imaging (rs-fMRI) allows the characterization of features designated as functional networks, functional connectivity gradients, and time-varying activity patterns that provide insight into the intrinsic functional organization of the brain and potential alterations related to brain dysfunction. Functional connectivity, hence, captures dimensions of the brain's activity that have enormous potential for both clinical and preclinical research. However, the mechanisms underlying functional connectivity have yet to be fully characterized, hindering interpretation of rs-fMRI studies. As in other branches of neuroscience, the identification of the neurophysiological processes that contribute to functional connectivity largely depends on research conducted on laboratory animals, which provide a platform where specific, multi-dimensional investigations that involve invasive measurements can be carried out. These highly controlled experiments facilitate the interpretation of the temporal correlations observed across the brain. Indeed, information obtained from animal experimentation to date is the basis for our current understanding of the underlying basis for functional brain connectivity. This review presents a compendium of some of the most critical advances in the field based on the efforts made by the animal neuroimaging community.

More than just summed neuronal activity: how multiple cell types shape the BOLD response

Functional neuroimaging techniques are widely applied to investigations of human cognition and disease. The most commonly used among these is blood oxygen level-dependent (BOLD) functional magnetic resonance imaging. The BOLD signal occurs because neural activity induces an increase in local blood supply to support the increased metabolism that occurs during activity. This supply usually outmatches demand, resulting in an increase in oxygenated blood in an active brain region, and a corresponding decrease in deoxygenated blood, which generates the BOLD signal. Hence, the BOLD response is shaped by an integration of local oxygen use, through metabolism, and supply, in the blood. To understand what information is carried in BOLD signals, we must understand how several cell types in the brain-local excitatory neurons, inhibitory neurons, astrocytes and vascular cells (pericytes, vascular smooth muscle and endothelial cells), and their modulation by ascending projection neurons-contribute to both metabolism and haemodynamic changes. Here, we review the contributions of each cell type to the regulation of cerebral blood flow and metabolism, and discuss situations where a simplified interpretation of the BOLD response as reporting local excitatory activity may misrepresent important biological phenomena, for example with regards to arousal states, ageing and neurological disease. This article is part of the theme issue 'Key relationships between non-invasive functional neuroimaging and the underlying neuronal activity'.

An equal-TE ultrafast 3D gradient-echo imaging method with high tolerance to magnetic susceptibility artifacts: Application to BOLD functional MRI

PURPOSE

To develop an ultrafast 3D gradient echo-based MRI method with constant TE and high tolerance to B₀ inhomogeneity, dubbed ERASE (equal-TE rapid acquisition with sequential excitation), and to introduce its use in BOLD functional MRI (fMRI).

THEORY AND METHODS

Essential features of ERASE, including spin behavior, were characterized, and a comparison study was conducted with conventional EPI. To demonstrate high tolerance to B₀ inhomogeneity, in vivo imaging of the mouse brain with a fiber-optic implant was performed at 9.4 T, and human brain imaging (including the orbitofrontal cortex) was performed at 3 T and 7 T. To evaluate the performance of ERASE in BOLD-fMRI, the characteristics of SNR and temporal SNR were analyzed for in vivo rat brains at 9.4 T in comparison with multislice gradient-echo EPI. Percent signal changes and t-scores are also presented.

RESULTS

For both mouse brain and human brain imaging, ERASE exhibited a high tolerance to magnetic susceptibility artifacts, showing much lower distortion and signal dropout, especially in the regions involving large magnetic susceptibility effects. For BOLD-fMRI, ERASE provided higher temporal SNR and t-scores than EPI, but exhibited similar percent signal changes in in vivo rat brains at 9.4 T.

CONCLUSION

When compared with conventional EPI, ERASE is much less sensitive, not only to EPI-related artifacts such as Nyquist ghosting, but also to B₀ inhomogeneity including magnetic susceptibility effects. It is promising for use in BOLD-fMRI, providing higher temporal SNR and t-scores with constant TE when compared with EPI, although further optimization is needed for human fMRI.

Animal Functional Magnetic Resonance Imaging: Trends and Path Toward Standardization

Animal whole-brain functional magnetic resonance imaging (fMRI) provides a non-invasive window into brain activity. A collection of associated methods aims to replicate observations made in humans and to identify the mechanisms underlying the distributed neuronal activity in the healthy and disordered brain. Animal fMRI studies have developed rapidly over the past years, fueled by the development of resting-state fMRI connectivity and genetically encoded neuromodulatory tools. Yet, comparisons between sites remain hampered by lack of standardization. Recently, we highlighted that mouse resting-state functional connectivity converges across centers, although large discrepancies in sensitivity and specificity remained. Here, we explore past and present trends within the animal fMRI community and highlight critical aspects in study design, data acquisition, and post-processing operations, that may affect the results and influence the comparability between studies. We also suggest practices aimed to promote the adoption of standards within the community and improve between-lab reproducibility. The implementation of standardized animal neuroimaging protocols will facilitate animal population imaging efforts as well as meta-analysis and replication studies, the gold standards in evidence-based science.

A CMOS NMR needle for probing brain physiology with high spatial and temporal resolution

Magnetic resonance imaging and spectroscopy are versatile methods for probing brain physiology, but their intrinsically low sensitivity limits the achievable spatial and temporal resolution. Here, we introduce a monolithically integrated NMR-on-a-chip needle that combines an ultra-sensitive 300 µm NMR coil with a complete NMR transceiver, enabling in vivo measurements of blood oxygenation and flow in nanoliter volumes at a sampling rate of 200 Hz.

Effects of Hemodynamic Response Function Selection on Rat fMRI Statistical Analyses

The selection of the appropriate hemodynamic response function (HRF) for signal modeling in functional magnetic resonance imaging (fMRI) is important. Although the use of the boxcar-shaped hemodynamic response function (BHRF) and canonical hemodynamic response (CHRF) has gained increasing popularity in rodent fMRI studies, whether the selected HRF affects the results of rodent fMRI has not been fully elucidated. Here we investigated the signal change and t-statistic sensitivities of BHRF, CHRF, and impulse response function (IRF). The effect of HRF selection on different tasks was analyzed by using data collected from two groups of rats receiving either 3 mA whisker pad or 3 mA forepaw electrical stimulations (n = 10 for each group). Under whisker pad stimulation with large blood-oxygen-level dependent (BOLD) signal change (4.31 ± 0.42%), BHRF significantly underestimated signal changes (P < 0.001) and t-statistics (P < 0.001) compared with CHRF or IRF. CHRF and IRF did not provide significantly different t-statistics (P > 0.05). Under forepaw stimulation with small BOLD signal change (1.71 ± 0.34%), different HRFs provided insignificantly different t-statistics (P > 0.05). Therefore, the selected HRF can influence data analysis in rodent fMRI experiments with large BOLD responses but not in those with small BOLD responses.

Validations of apomorphine-induced BOLD activation correlations in hemiparkinsonian rhesus macaques

Identification of Parkinson's disease at the earliest possible stage of the disease may provide the best opportunity for the use of disease modifying treatments. However, diagnosing the disease during the pre-symptomatic period remains an unmet goal. To that end, we used pharmacological MRI (phMRI) to assess the function of the cortico-basal ganglia circuit in a non-human primate model of dopamine deficiency to determine the possible relationships between phMRI signals with behavioral, neurochemical, and histological measurements. Animals with unilateral treatments with the neurotoxin, 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP), that expressed stable, long-term hemiparkinsonism were challenged with the dopaminergic receptor agonist, apomorphine, and structure-specific phMRI blood oxygen level-dependent (BOLD) activation responses were measured. Behavioral, histopathological, and neurochemical measurements were obtained and correlated with phMRI activation of structures of the cortico-basal ganglia system. Greater phMRI activations in the basal ganglia and cortex were associated with slower movement speed, decreased daytime activity, or more pronounced parkinsonian features. Animals showed decreased stimulus-evoked dopamine release in the putamen and substantia nigra pars compacta and lower basal glutamate levels in the motor cortex on the MPTP-lesioned hemisphere compared to the contralateral hemisphere. The altered neurochemistry was significantly correlated with phMRI signals in the motor cortex and putamen. Finally, greater phMRI activations in the caudate nucleus correlated with fewer tyrosine hydroxylase-positive (TH⁺) nigral cells and decreased TH⁺ fiber density in the putamen. These results reveal the correlation of phMRI signals with the severity of the motor deficits and pathophysiological changes in the cortico-basal ganglia circuit.

•

Apomorphine in hemiparkinsonian animals can evoke changes in functional MRI signals.

•

Cortico-basal ganglia activation correlates to behavior, neurochemistry, histology

•

Pharmacological MRI has potential to be biomarker for Parkinson's disease.

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.

Gradient-echo and spin-echo blood oxygenation level-dependent functional MRI at ultrahigh fields of 9.4 and 15.2 Tesla

Purpose

Sensitivity and specificity of blood oxygenation level–dependent (BOLD) functional MRI (fMRI) is sensitive to magnetic field strength and acquisition methods. We have investigated gradient‐echo (GE)‐ and spin‐echo (SE)‐BOLD fMRI at ultrahigh fields of 9.4 and 15.2  Tesla.

Methods

BOLD fMRI experiments responding to forepaw stimulation were performed with 3 echo times (TE) at each echo type and B₀ in α‐chloralose–anesthetized rats. The contralateral forelimb somatosensory region was selected for quantitative analyses.

Results

At 9.4 T and 15.2 T, average baseline T₂^* (n = 9) was 26.6 and 17.1 msec, whereas baseline T₂ value (n = 9) was 35.7 and 24.5 msec, respectively. Averaged stimulation‐induced ΔR₂^* was –1.72 s^–1 at 9.4 T and –3.09 s^–1 at 15.2 T, whereas ΔR₂ was –1.19 s^–1 at 9.4 T and –1.97 s^–1 at 15.2 T. At the optimal TE of tissue T₂^* or T₂, BOLD percent changes were slightly higher at 15.2 T than at 9.4 T (GE: 7.4% versus 6.4% and SE: 5.7% versus 5.4%). The ΔR₂^* and ΔR₂ ratio of 15.2 T to 9.4 T was 1.8 and 1.66, respectively. The ratio of the macrovessel‐containing superficial to microvessel‐dominant parenchymal BOLD signal was 1.73 to 1.76 for GE‐BOLD versus 1.13 to 1.19 for SE‐BOLD, indicating that the SE‐BOLD contrast is less sensitive to macrovessels than GE‐BOLD.

Conclusion

SE‐BOLD fMRI improves spatial specificity to microvessels compared to GE‐BOLD at both fields. BOLD sensitivity is similar at the both fields and can be improved at ultrahigh fields only for thermal‐noise–dominant ultrahigh‐resolution fMRI.

The cerebral blood flow response dependency on stimulus pulse width is affected by stimulus current amplitude - a study of activation flow coupling

The coupling of cerebral blood flow (CBF) to neuronal activation, referred to as activation flow coupling (AFC), has been a fundamental brain physiology property. The stimulus-evoked CBF response was usually considered as a surrogate marker for neuronal activity in AFC studies. The selection of appropriate stimulation parameters, e.g., current amplitude and pulse width, is of great importance yet the effect of pulse width changes remained contradictory in previous studies. In this work, we use laser speckle contrast imaging (LSCI) to study the spatiotemporal CBF response to hindpaw somatosensory stimulation of different pulse widths (0.3 ms vs 1 ms) and current amplitudes (3 mA vs 6 mA) in a rodent experiment. The results showed that the change of pulse width significantly affected the CBF peak value at a lower current level (p<;0.05). In addition, the duration for observing significantly different average CBF response, denoted as td, at various pulse widths, was dependent on stimulus current amplitude. At a lower amplitude (3 mA), td was 6.5 s; While at a higher amplitude (6 mA), td was 2.5 s. It was indicated that the changes of pulse width had longer influence on the average CBF response at a lower current amplitude. Our findings may help to understand and explain the inconsistent AFC with different stimulation parameters in fundamental brain physiology.

Longitudinal, transcranial measurement of functional activation in the rat brain by diffuse correlation spectroscopy

Neural activity is an important biomarker for the presence of neurodegenerative diseases, cerebrovascular alterations, and brain trauma; furthermore, it is a surrogate marker for treatment effects. These pathologies may occur and evolve in a long time-period, thus, noninvasive, transcutaneous techniques are necessary to allow a longitudinal follow-up. In the present work, we have customized noninvasive, transcutaneous, diffuse correlation spectroscopy (DCS) to localize changes in cerebral blood flow (CBF) induced by neural activity. We were able to detect changes in CBF in the somatosensory cortex by using a model of electrical forepaw stimulation in rats. The suitability of DCS measurements for longitudinal monitoring was demonstrated by performing multiple sessions with the same animals at different ages (from 6 to 18 months). In addition, functional DCS has been cross-validated by comparison with functional magnetic resonance imaging (fMRI) in the same animals in a subset of the time-points. The overall results obtained with transcutaneous DCS demonstrates that it can be utilized in longitudinal studies safely and reproducibly to locate changes in CBF induced by neural activity in the small animal brain.

Phase imaging with multiple phase-cycled balanced steady-state free precession at 9.4 T

While phase imaging with a gradient echo (GRE) sequence is popular, phase imaging with balanced steady-state free precession (bSSFP) has been underexplored. The purpose of this study was to investigate anatomical and functional phase imaging with multiple phase-cycled bSSFP, in expectation of increasing spatial coverage of steep phase-change regions of bSSFP. Eight different dynamic 2D pass-band bSSFP studies at four phase-cycling (PC) angles and two T_E /T_R values were performed on rat brains at 9.4 T with electrical forepaw stimulation, in comparison with dynamic 2D GRE. Anatomical and functional phase images were obtained by averaging the dynamic phase images and mapping correlation between the dynamic images and the stimulation paradigm, and were compared with their corresponding magnitude images. Phase imaging with 3D pass-band and 3D transition-band bSSFP was also performed for comparison with 3D GRE phase imaging. Two strategies of combining the multiple phase-cycled bSSFP phase images were also proposed. Contrast between white matter and gray matter in bSSFP phase images significantly varied with PC angle and became twice as high as that of GRE phase images at a specific PC angle. With the same total scan time, the combined bSSFP phase images provided stronger phase contrast and visualized neuronal fiber-like structures more clearly than the GRE phase images. The combined phase images of both 3D pass-band and 3D transition-band bSSFP showed phase contrasts stronger than those of the GRE phase images in overall brain regions, even at a longer T_E of 20 ms. In contrast, phase functional MRI (fMRI) signals were weak overall and mostly located in draining veins for both bSSFP and GRE. Multiple phase-cycled bSSFP phase imaging is a promising anatomical imaging technique, while its usage as fMRI does not seem desirable with the current approach.

Lactate and glutamate dynamics during prolonged stimulation of the rat barrel cortex suggest adaptation of cerebral glucose and oxygen metabolism

A better understanding of BOLD responses stems from a better characterization of the brain's ability to metabolize glucose and oxygen. Non-invasive techniques such as functional magnetic resonance spectroscopy (fMRS) have thus been developed allowing for the reproducible assessment of metabolic changes during barrel cortex (S1BF) activations in rats. The present study aimed at further exploring the role of neurotransmitters on local and temporal changes in vascular and metabolic function in S1BF. fMRS and fMRI data were acquired sequentially in α-chloralose anesthetized rats during 32-min rest and trigeminal nerve stimulation periods. During stimulation, concentrations of lactate (Lac) and glutamate (Glu) increased in S1BF by 0.23±0.05 and 0.34±0.05μmol/g respectively in S1BF. Dynamic analysis of metabolite concentrations allowed estimating changes in cerebral metabolic rates of glucose (ΔCMR_Glc) and oxygen (ΔCMR_O2). Findings confirmed a prevalence of oxidative metabolism during prolonged S1BF activation. Habituation led to a significant BOLD magnitude decline as a function of time while both total ΔCMR_Glc and ΔCMR_O2 remained constant revealing adaptation of glucose and oxygen metabolisms to support ongoing trigeminal nerve stimulation.

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

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

Arterial spin labeling perfusion magnetic resonance imaging of non-human primates

Non-human primates (NHPs) resemble most aspects of humans in brain physiology and anatomy and are excellent animal models for translational research in neuroscience, biomedical research and pharmaceutical development. Cerebral blood flow (CBF) offers essential physiological information of the brain to examine the abnormal functionality in NHP models with cerebral vascular diseases and neurological disorders or dementia. Arterial spin labeling (ASL) perfusion MRI techniques allow for high temporal and spatial CBF measurement and are intensively used in studies of animals and humans. In this article, current high-resolution ASL perfusion MRI techniques for quantitative evaluation of brain physiology and function in NHPs are described and their applications and limitation are discussed as well.

Localization of cortical primary motor area of the hand using navigated transcranial magnetic stimulation, BOLD and arterial spin labeling fMRI

BACKGROUND

Although the relationship between neuronavigated transcranial magnetic stimulation (nTMS) and functional magnetic resonance imaging (fMRI) has been widely studied in motor mapping, it is unknown how the motor response type or the choice of motor task affect this relationship.

NEW METHOD

Centers of gravity (CoGs) and response maxima were measured with blood-oxygen-level dependent (BOLD) and arterial spin labeling (ASL) fMRI during motor tasks against nTMS CoGs and response maxima, which were mapped with motor evoked potentials (MEPs) and silent periods (SPs).

RESULTS

No differences in motor representations (CoGs and response maxima) were observed in lateral-medial direction (p=0.265). fMRI methods localized the motor representation more posterior than nTMS (p<0.001). This was not affected by the BOLD fMRI motor task (p>0.999) nor nTMS response type (p>0.999). ASL fMRI maxima did not differ from the nTMS nor BOLD fMRI CoGs (p≥0.070), but the ASL CoG was deeper in comparison to other methods (p≤0.042). The BOLD fMRI motor task did not influence the depth of the motor representation (p≥0.745). The median Euclidean distances between the nTMS and fMRI motor representations varied between 7.7mm and 14.5mm and did not differ between the methods (F≤1.23, p≥0.318).

COMPARISON WITH EXISTING METHODS

The relationship between fMRI and nTMS mapped excitatory (MEP) and inhibitory (SP) responses, and whether the choice of motor task affects this relationship, have not been studied before.

CONCLUSIONS

The congruence between fMRI and nTMS is good. The choice of nTMS motor response type nor BOLD fMRI motor task had no effect on this relationship.

Novel method for functional brain imaging in awake minimally restrained rats

Functional magnetic resonance imaging (fMRI) in rodents holds great promise for advancing our knowledge about human brain function. However, the use of anesthetics to immobilize rodents during fMRI experiments has restricted the type of questions that can be addressed using this technique. Here we describe an innovative procedure to train rats to be constrained without the need of any anesthesia during the whole procedure. We show that with 8-10 days of acclimation rats can be conscious and remain still during fMRI experiments under minimal stress. In addition, we provide fMRI results of conscious rodents in a variety of commonly used fMRI experimental paradigms, and we demonstrate the improved quality of these scans by comparing results when the same rodents were scanned under anesthesia. We confirm that the awake scanning procedure permits an improved evaluation of brain networks and brain response to external stimuli with minimal movement artifact. The present study further advances the field of fMRI in awake rodents, which provide more direct, forward and reverse, translational opportunities regarding brain functional correspondences between human and rodent research.

T2*-weighted fMRI time-to-peak of oxygen challenge in ischemic stroke

T2 (*)-weighted MRI of transient oxygen challenge (OC) showed exaggerated OC percent changes in the ischemic tissue at risk compared to normal tissue. One ambiguity is that regions with high vascular density also showed exaggerated OC percent changes. This study explored time-to-peak (TTP) of the OC percent changes to improve the utility of T2 (*)-weighted OC MRI. Experiments were performed longitudinally at 30 min, 150 min and 24 h after transient (60-min) stroke in rats. Ischemic core, normal, and mismatch tissue were classified pixel-by-pixel based on apparent diffusion coefficient and cerebral blood flow. Major findings were: (i) Delayed OC TTP was localized to and corresponded well with the perfusion-diffusion mismatch. (ii) By contrast, the exaggerated OC percent changes were less localized, with changes not only in the at-risk tissue but also in some areas of the contralesional hemisphere with venous vessel origins. (iii) The OC time-course of the mismatch tissue was biphasic, with a faster initial increase followed by a slower increase. (iv) At-risk tissue with delayed TTP and exaggerated OC was normal after reperfusion and the at-risk tissue was mostly (83 ± 18%) rescued by reperfusion as indicated by normal 24-h T2. OC TTP offers unique information toward better characterization of at-risk tissue in ischemic stroke.

Magnetic Resonance Imaging in Experimental Traumatic Brain Injury

Traumatic brain injury (TBI) is a leading cause of death and disability in the USA. Common causes of TBI include falls, violence, injuries from wars, and vehicular and sporting accidents. The initial direct mechanical damage in TBI is followed by progressive secondary injuries such as brain swelling, perturbed cerebral blood flow (CBF), abnormal cerebrovascular reactivity (CR), metabolic dysfunction, blood-brain-barrier disruption, inflammation, oxidative stress, and excitotoxicity, among others. Magnetic resonance imaging (MRI) offers the means to noninvasively probe many of these secondary injuries. MRI has been used to image anatomical, physiological, and functional changes associated with TBI in a longitudinal manner. This chapter describes controlled cortical impact (CCI) TBI surgical procedures, a few common MRI protocols used in TBI imaging, and, finally, image analysis pertaining to experimental TBI imaging in rats.

Quantitative β mapping for calibrated fMRI

The metabolic and hemodynamic dependencies of the blood oxygenation level-dependent (BOLD) signal form the basis for calibrated fMRI, where the focus is on oxidative energy demanded by neural activity. An important part of calibrated fMRI is the power-law relationship between the BOLD signal and the deoxyhemoglobin concentration, which in turn is related to the ratio between oxidative demand (CMRO2) and blood flow (CBF). The power-law dependence between BOLD signal and deoxyhemoglobin concentration is signified by a scaling exponent β. Until recently most studies assumed a β value of 1.5, which is based on numerical simulations of the extravascular BOLD component. Since the basal value of CMRO2 and CBF can vary from subject-to-subject and/or region-to-region, a method to independently measure β in vivo should improve the accuracy of calibrated fMRI results. We describe a new method for β mapping through characterizing R2' - the most sensitive relaxation component of BOLD signal (i.e., the reversible magnetic susceptibility component that is predominantly of extravascular origin at high magnetic field) - as a function of intravascular magnetic susceptibility induced by an FDA-approved superparamagnetic contrast agent. In α-chloralose anesthetized rat brain, at 9.4 T, we measured β values of ~0.8 uniformly across large neocortical swathes, with lower magnitude and more heterogeneity in subcortical areas. Comparison of β maps in rats anesthetized with medetomidine and α-chloralose revealed that β is independent of neural activity levels at these resting states. We anticipate that this method for β mapping can help facilitate calibrated fMRI for clinical studies.

Ipsi- and contralateral frontal cortex oxygenation during handgrip task does not follow decrease on maximal force output

The effect of fatiguing exercise on the ipsi- and contralateral frontal cortex has not been fully clarified. The purpose of this study was to investigate by near-infrared spectroscopy (NIRS) the frontal cortex oxygenation response to a prolonged fatiguing repetitive handgrip exercise performed at maximal voluntary contraction. It was found a significant oxyhemoglobin concentration ([HbO₂]) increase (p < 0.05), accompanied by a smaller and delayed deoxyhemoglobin concentration ([Hb]) decrease (p < 0.05), in both hemispheres. Then, it was indicated higher delayed oxygenation in ipsilateral oxygenation compared to contralateral oxygenation. These results provide further evidence that the complemental interaction between the ipsilateral and contralateral cortex during the fatiguing maximal exercise.

Compressed Sensing for fMRI: Feasibility Study on the Acceleration of Non-EPI fMRI at 9.4T

Conventional functional magnetic resonance imaging (fMRI) technique known as gradient-recalled echo (GRE) echo-planar imaging (EPI) is sensitive to image distortion and degradation caused by local magnetic field inhomogeneity at high magnetic fields. Non-EPI sequences such as spoiled gradient echo and balanced steady-state free precession (bSSFP) have been proposed as an alternative high-resolution fMRI technique; however, the temporal resolution of these sequences is lower than the typically used GRE-EPI fMRI. One potential approach to improve the temporal resolution is to use compressed sensing (CS). In this study, we tested the feasibility of k-t FOCUSS—one of the high performance CS algorithms for dynamic MRI—for non-EPI fMRI at 9.4T using the model of rat somatosensory stimulation. To optimize the performance of CS reconstruction, different sampling patterns and k-t FOCUSS variations were investigated. Experimental results show that an optimized k-t FOCUSS algorithm with acceleration by a factor of 4 works well for non-EPI fMRI at high field under various statistical criteria, which confirms that a combination of CS and a non-EPI sequence may be a good solution for high-resolution fMRI at high fields.

Neural and hemodynamic responses to optogenetic and sensory stimulation in the rat somatosensory cortex

Introducing optogenetics into neurovascular research can provide novel insights into the cell-specific control of the hemodynamic response. To generalize findings from molecular approaches, it is crucial to determine whether light-activated circuits have the same effect on the vasculature as sensory-activated ones. For that purpose, rats expressing channelrhodopsin (ChR2) specific to excitatory glutamatergic neurons were used to measure neural activity, blood flow, hemoglobin-based optical intrinsic signal, and blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI) during optogenetic and sensory stimulation. The magnitude of the evoked hemodynamic responses was monotonically correlated with optogenetic stimulus strength. The BOLD hemodynamic response function was consistent for optogenetic and sensory stimuli. The relationship between electrical activities and hemodynamic responses was comparable for optogenetic and sensory stimuli, and better explained by the local field potential (LFP) than the firing rate. The LFP was well correlated with cerebral blood flow, moderately with cerebral blood volume, and less with deoxyhemoglobin (dHb) level. The presynaptic firing rate had little impact on evoking vascular response. Contribution of the postsynaptic LFP to the blood flow response induced by optogenetic stimulus was further confirmed by the application of glutamate receptor antagonists. Overall, neurovascular coupling during optogenetic control of glutamatergic neurons largely conforms to that of a sensory stimulus.

Imaging of prolonged BOLD response in the somatosensory cortex of the rat

Blood oxygenation level-dependent (BOLD) functional MRI is a widely employed methodology in experimental and clinical neuroscience, although its nature is not fully understood. To gain insights into BOLD mechanisms and take advantage of the new functional methods, it is of interest to investigate prolonged paradigms of activation suitable for long experimental protocols and to observe any long-term modifications induced by these functional challenges. While different types of sustained stimulation paradigm have been explored in human studies, the BOLD response is typically limited to a few minutes in animal models, due to fatigue, anesthesia effects and physiological instability. In the present study, the rat forepaw was electrically stimulated for 2 h, which resulted in a prolonged and localized cortical BOLD response over that period. The stimulation paradigm, including an inter-stimulus interval (ISI) of 10 s, that is 25% of the total time, was applied at constant or variable frequency over 2 h. The steady-state level of the BOLD response was reached after 15-20 min of stimulation and was maintained until the end of the stimulation. On average, no substantial loss in activated volume was observed at the end of the stimulation, but less variability in the fraction of remaining activated volume and higher steady-state BOLD amplitude were observed when stimulation frequency was varied between 2 and 3 Hz every 5 min. We conclude that the combination of ISI and variable stimulus frequency reproducibly results in robust, prolonged and localized BOLD activation.

Imaging methodologies and applications for nutrition research: what can functional MRI offer?

Food intake is influenced by a complex regulatory system involving the integration of a wide variety of sensory inputs across multiple brain areas. Over the past decade, advances in neuroimaging using functional MRI (fMRI) have provided valuable insight into these pathways in the human brain. This review provides an outline of the methodology of fMRI, introducing the widely used blood oxygenation level-dependent contrast for fMRI and direct measures of cerebral blood flow using arterial spin labelling. A review of fMRI studies of the brain's response to taste, aroma and oral somatosensation, and how fat is sensed and mapped in the brain in relation to the pleasantness of food, and appetite control is given. The influence of phenotype on individual variability in cortical responses is addressed, and an overview of fMRI studies investigating hormonal influences (e.g. peptide YY, cholecystokinin and ghrelin) on appetite-related brain processes provided. Finally, recent developments in MR technology at ultra-high field (7 T) are introduced, highlighting the advances this can provide for fMRI studies to investigate the neural underpinnings in nutrition research. In conclusion, neuroimaging methods provide valuable insight into the mechanisms of flavour perception and appetite behaviour.

Physiological measurements using ultra-high field fMRI: a review

Functional MRI (fMRI) has grown to be the neuroimaging technique of choice for investigating brain function. This topical review provides an outline of fMRI methods and applications, with a particular emphasis on the recent advances provided by ultra-high field (UHF) scanners to allow functional mapping with greater sensitivity and improved spatial specificity. A short outline of the origin of the blood oxygenation level dependent (BOLD) contrast is provided, followed by a review of BOLD fMRI methods based on gradient-echo (GE) and spin-echo (SE) contrast. Phase based fMRI measures, as well as perfusion contrast obtained with the technique of arterial spin labelling (ASL), are also discussed. An overview of 7 T based functional neuroimaging is provided, outlining the potential advances to be made and technical challenges to be addressed.

Comparing like with like: the power of knowing where you are

Magnetic resonance imaging can now provide human brain images of structure, function, and connectivity with isotropic voxels smaller than one millimeter, and thus much smaller than the cortical thickness. This resolution, achievable in a scan time of less than 1 h, enables visualization of myeloarchitectural layer structure, intracortical variations in functional activity--recorded in changes in blood oxygenation level dependent signal or cerebral blood volume CBV--and intracortical axonal orientational structure via diffusion-weighted magnetic resonance imaging. While recent improvements in radiofrequency receiver coils now enable excellent image data to be obtained at 3T, scanning at the ultra-high field of 7 T offers further gains in signal-to-noise ratio and speed of image acquisition, with a structural image resolution of about 300 μm. These improvements throw into sharp question the strategies that have become conventional for the analysis of functional imaging data, especially the practice of spatial smoothing of raw functional data before further analysis. Creation of a native cortical map for each human subject that provides a reliable individual parcellation into cortical areas related to Brodmann Areas enables a strikingly different approach to functional image analysis. This proposed approach involves surface registration of the cortices of groups of subjects using maps of the longitudinal relaxation time T1 as an index of myelination, and methods for inferring statistical significance that do not entail spatial smoothing. The outcome should be a far more precise comparison of like-with-like cortical areas across subjects, with the potential to greatly increase experimental power, to discriminate activity in neighboring cortical areas, and to enable correlation of function and connectivity with specific cytoarchitecture. Such analyses should enable a far more convincing modeling of brain mechanisms than current graph-based methods that require gross over-simplification of brain activity patterns in order to be computationally tractable.

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.

Specificity of stimulus-evoked fMRI responses in the mouse: the influence of systemic physiological changes associated with innocuous stimulation under four different anesthetics

Functional magnetic resonance (fMRI) in mice has become an attractive tool for mechanistic studies, for characterizing models of human disease, and for evaluation of novel therapies. Yet, controlling the physiological state of mice is challenging, but nevertheless important as changes in cardiovascular parameters might affect the hemodynamic readout which constitutes the basics of the fMRI signal. In contrast to rats, fMRI studies in mice report less robust brain activation of rather widespread character to innocuous sensory stimulation. Anesthesia is known to influence the characteristics of the fMRI signal. To evaluate modulatory effects imposed by the anesthesia on stimulus-evoked fMRI responses, we compared blood oxygenation level dependent (BOLD) and cerebral blood volume (CBV) signal changes to electrical hindpaw stimulation using the four commonly used anesthetics isoflurane, medetomidine, propofol and urethane. fMRI measurements were complemented by assessing systemic physiological parameters throughout the experiment. Unilateral stimulation of the hindpaw elicited widespread fMRI responses in the mouse brain displaying a bilateral pattern irrespective of the anesthetic used. Analysis of magnitude and temporal profile of BOLD and CBV signals indicated anesthesia-specific modulation of cerebral hemodynamic responses and differences observed for the four anesthetics could be largely explained by their known effects on animal physiology. Strikingly, independent of the anesthetic used our results reveal that fMRI responses are influenced by stimulus-induced cardiovascular changes, which indicate an arousal response, even to innocuous stimulation. This may mask specific fMRI signal associated to the stimulus. Hence, studying the processing of peripheral input in mice using fMRI techniques constitutes a major challenge and adapted paradigms and/or alternative fMRI readouts should also be considered when studying sensory processing in mice.

Quantitative basis for neuroimaging of cortical laminae with calibrated functional MRI

Layer-specific neurophysiologic, hemodynamic, and metabolic measurements are needed to interpret high-resolution functional magnetic resonance imaging (fMRI) data in the cerebral cortex. We examined how neurovascular and neurometabolic couplings vary vertically in the rat's somatosensory cortex. During sensory stimulation we measured dynamic layer-specific responses of local field potential (LFP) and multiunit activity (MUA) as well as blood oxygenation level-dependent (BOLD) signal and cerebral blood volume (CBV) and blood flow (CBF), which in turn were used to calculate changes in oxidative metabolism (CMR(O2)) with calibrated fMRI. Both BOLD signal and CBV decreased from superficial to deep laminae, but these responses were not well correlated with either layer-specific LFP or MUA. However, CBF changes were quite stable across laminae, similar to LFP. However, changes in CMR(O2) and MUA varied across cortex in a correlated manner and both were reduced in superficial lamina. These results lay the framework for quantitative neuroimaging across cortical laminae with calibrated fMRI methods.

Neural and hemodynamic responses elicited by forelimb- and photo-stimulation in channelrhodopsin-2 mice: insights into the hemodynamic point spread function

Hemodynamic responses are commonly used to map brain activity; however, their spatial limits have remained unclear because of the lack of a well-defined and malleable spatial stimulus. To examine the properties of neural activity and hemodynamic responses, multiunit activity, local field potential, cerebral blood volume (CBV)-sensitive optical imaging, and laser Doppler flowmetry were measured from the somatosensory cortex of transgenic mice expressing Channelrhodopsin-2 in cortex Layer 5 pyramidal neurons. The magnitude and extent of neural and hemodynamic responses were modulated using different photo-stimulation parameters and compared with those induced by somatosensory stimulation. Photo-stimulation-evoked spiking activity across cortical layers was similar to forelimb stimulation, although their activity originated in different layers. Hemodynamic responses induced by forelimb- and photo-stimulation were similar in magnitude and shape, although the former were slightly larger in amplitude and wider in extent. Altogether, the neurovascular relationship differed between these 2 stimulation pathways, but photo-stimulation-evoked changes in neural and hemodynamic activities were linearly correlated. Hemodynamic point spread functions were estimated from the photo-stimulation data and its full-width at half-maximum ranged between 103 and 175 µm. Therefore, submillimeter functional structures separated by a few hundred micrometers may be resolved using hemodynamic methods, such as optical imaging and functional magnetic resonance imaging.

Pain response measured with arterial spin labeling

The majority of functional MRI studies of pain processing in the brain use the blood oxygenation level-dependent (BOLD) imaging approach. However, the BOLD signal is complex as it depends on simultaneous changes in blood flow, vascular volume and oxygen metabolism. Arterial spin labeling (ASL) perfusion imaging is another imaging approach in which the magnetically labeled arterial water is used as an endogenous tracer that allows for direct measurement of cerebral blood flow. In this study, we assessed the pain response in the brain using a pulsed-continuous arterial spin labeling (pCASL) approach and a thermal stimulation paradigm. Using pCASL, response to noxious stimulation was detected in somatosensory cortex, anterior cingulate cortex, anterior insula, hippocampus, amygdala, thalamus and precuneus, consistent with the pain response activation patterns detected using the BOLD imaging approach. We suggest that pCASL is a reliable alternative for functional MRI pain studies in conditions in which blood flow, volume or oxygen extraction are altered or compromised.

A temperature-stable cryo-system for high-temperature superconducting MR in-vivo imaging

To perform a rat experiment using a high-temperature superconducting (HTS) surface resonator, a cryostat is essential to maintain the rat's temperature. In this work, a compact temperature-stable HTS cryo-system, keeping animal rectal temperature at 37.4°C for more than 3 hours, was successfully developed. With this HTS cryo-system, a 40-mm-diameter Bi₂Sr₂Ca₂Cu₃O_x (Bi-2223) surface resonator at 77 K was demonstrated in a 3-Tesla MRI system. The proton resonant frequency (PRF) method was employed to monitor the rat's temperature. Moreover, the capacity of MR thermometry in the HTS experiments was evaluated by correlating with data from independent fiber-optic sensor temperature measurements. The PRF thermal coefficient was derived as 0.03 rad/°C and the temperature-monitoring architecture can be implemented to upgrade the quality and safety in HTS experiments. The signal-to-noise ratio (SNR) of the HTS surface resonator at 77 K was higher than that of a professionally made copper surface resonator at 300 K, which has the same geometry, by a 3.79-fold SNR gain. Furthermore, the temperature-stable HTS cryo-system we developed can obtain stable SNR gain in every scan. A temperature-stable HTS cryo-system with an external air-blowing circulation system is demonstrated.

fMRI in the awake marmoset: somatosensory-evoked responses, functional connectivity, and comparison with propofol anesthesia

Functional neuroimaging in animal models is essential for understanding the principles of neurovascular coupling and the physiological basis of fMRI signals that are widely used to study sensory and cognitive processing in the human brain. While hemodynamic responses to sensory stimuli have been characterized in humans, animal studies are able to combine very high resolution imaging with invasive measurements and pharmacological manipulation. To date, most high-resolution studies of neurovascular coupling in small animals have been carried out in anesthetized rodents. Here we report fMRI experiments in conscious, awake common marmosets (Callithrix jacchus), and compare responses to animals anesthetized with propofol. In conscious marmosets, robust BOLD fMRI responses to somatosensory stimulation of the forearm were found in contralateral and ipsilateral regions of the thalamus, primary (SI) and secondary (SII) somatosensory cortex, and the caudate nucleus. These responses were markedly stronger than those in anesthetized marmosets and showed a monotonic increase in the amplitude of the BOLD response with stimulus frequency. On the other hand, anesthesia significantly attenuated responses in thalamus, SI and SII, and abolished responses in caudate and ipsilateral SI. Moreover, anesthesia influenced several other aspects of the fMRI responses, including the shape of the hemodynamic response function and the interareal (SI-SII) spontaneous functional connectivity. Together, these findings demonstrate the value of the conscious, awake marmoset model for studying physiological responses in the somatosensory pathway, in the absence of anesthesia, so that the data can be compared most directly to fMRI in conscious humans.

Methylene blue potentiates stimulus-evoked fMRI responses and cerebral oxygen consumption during normoxia and hypoxia

Methylene blue USP (MB) at low doses has metabolic-enhancing and antioxidant properties and exhibits experimental neurotherapeutic benefits, but little is known about its in vivo effects on cerebral blood flow (CBF), functional evoked responses, and the associated changes in cerebral metabolic rate of oxygen (CMRO2). This study used magnetic resonance imaging (MRI) to evaluate the in vivo effects of a single intravenous MB therapeutic dose (0.5mg/kg) on basal CBF, blood oxygenation level-dependent (BOLD) and CBF responses to hypercapnic (5% CO2 in air) inhalation, as well as changes in BOLD, CBF, and CMRO2 during forepaw stimulation in the rat brain. MB did not have significant effects on arterial oxygen saturation, heart rate and fMRI responses to hypercapnia. However, MB significantly potentiated forepaw-evoked BOLD and CBF changes under normoxia. To further evaluate in vivo effects of MB under metabolic stress conditions, MRI measurements were also made under mild hypoxia (15% O2). Hypoxia per se increased evoked functional MRI responses. MB under hypoxia further potentiated forepaw-evoked BOLD, CBF and oxygen consumption responses relative to normoxia. These findings provide insights into MB's effects on cerebral hemodynamics in vivo and could help to optimize treatments in neurological diseases with mitochondrial dysfunction and oxidative stress.

Optimization of flow-sensitive alternating inversion recovery (FAIR) for perfusion functional MRI of rodent brain

Arterial spin labeling (ASL) MRI provides a noninvasive method to image perfusion, and has been applied to map neural activation in the brain. Although pulsed labeling methods have been widely used in humans, continuous ASL with a dedicated neck labeling coil is still the preferred method in rodent brain functional MRI (fMRI) to maximize the sensitivity and allow multislice acquisition. However, the additional hardware is not readily available and hence its application is limited. In this study, flow-sensitive alternating inversion recovery (FAIR) pulsed ASL was optimized for fMRI of rat brain. A practical challenge of FAIR is the suboptimal global inversion by the transmit coil of limited dimensions, which results in low effective labeling. By using a large volume transmit coil and proper positioning to optimize the body coverage, the perfusion signal was increased by 38.3% compared with positioning the brain at the isocenter. An additional 53.3% gain in signal was achieved using optimized repetition and inversion times compared with a long TR. Under electrical stimulation to the forepaws, a perfusion activation signal change of 63.7 ± 6.3% can be reliably detected in the primary somatosensory cortices using single slice or multislice echo planar imaging at 9.4 T. This demonstrates the potential of using pulsed ASL for multislice perfusion fMRI in functional and pharmacological applications in rat brain.

Anesthesia and the quantitative evaluation of neurovascular coupling

Anesthesia has broad actions that include changing neuronal excitability, vascular reactivity, and other baseline physiologies and eventually modifies the neurovascular coupling relationship. Here, we review the effects of anesthesia on the spatial propagation, temporal dynamics, and quantitative relationship between the neural and vascular responses to cortical stimulation. Previous studies have shown that the onset latency of evoked cerebral blood flow (CBF) changes is relatively consistent across anesthesia conditions compared with variations in the time-to-peak. This finding indicates that the mechanism of vasodilation onset is less dependent on anesthesia interference, while vasodilation dynamics are subject to this interference. The quantitative coupling relationship is largely influenced by the type and dosage of anesthesia, including the actions on neural processing, vasoactive signal transmission, and vascular reactivity. The effects of anesthesia on the spatial gap between the neural and vascular response regions are not fully understood and require further attention to elucidate the mechanism of vascular control of CBF supply to the underlying focal and surrounding neural activity. The in-depth understanding of the anesthesia actions on neurovascular elements allows for better decision-making regarding the anesthetics used in specific models for neurovascular experiments and may also help elucidate the signal source issues in hemodynamic-based neuroimaging techniques.

Investigation of the cerebral hemodynamic response function in single blood vessels by functional photoacoustic microscopy

The specificity of the hemodynamic response function (HRF) is determined spatially by the vascular architecture and temporally by the evolution of hemodynamic changes. Here, we used functional photoacoustic microscopy (fPAM) to investigate single cerebral blood vessels of rats after left forepaw stimulation. In this system, we analyzed the spatiotemporal evolution of the HRFs of the total hemoglobin concentration (HbT), cerebral blood volume (CBV), and hemoglobin oxygen saturation (SO(2)). Changes in specific cerebral vessels corresponding to various electrical stimulation intensities and durations were bilaterally imaged with 36 × 65-μm(2) spatial resolution. Stimulation intensities of 1, 2, 6, and 10 mA were applied for periods of 5 or 15 s. Our results show that the relative functional changes in HbT, CBV, and SO(2) are highly dependent not only on the intensity of the stimulation, but also on its duration. Additionally, the duration of the stimulation has a strong influence on the spatiotemporal characteristics of the HRF as shorter stimuli elicit responses only in the local vasculature (smaller arterioles), whereas longer stimuli lead to greater vascular supply and drainage. This study suggests that the current fPAM system is reliable for studying relative cerebral hemodynamic changes, as well as for offering new insights into the dynamics of functional cerebral hemodynamic changes in small animals.

Quantitative fMRI and oxidative neuroenergetics

The discovery of functional magnetic resonance imaging (fMRI) has greatly impacted neuroscience. The blood oxygenation level-dependent (BOLD) signal, using deoxyhemoglobin as an endogenous paramagnetic contrast agent, exposes regions of interest in task-based and resting-state paradigms. However the BOLD contrast is at best a partial measure of neuronal activity, because the functional maps obtained by differencing or correlations ignore the total neuronal activity in the baseline state. Here we describe how studies of brain energy metabolism at Yale, especially with (13)C magnetic resonance spectroscopy and related techniques, contributed to development of quantitative functional brain imaging with fMRI by providing a reliable measurement of baseline energy. This narrative takes us on a journey, from molecules to mind, with illuminating insights about neuronal-glial activities in relation to energy demand of synaptic activity. These results, along with key contributions from laboratories worldwide, comprise the energetic basis for quantitative interpretation of fMRI data.

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.

Contributions of dynamic venous blood volume versus oxygenation level changes to BOLD fMRI

Blood-oxygenation-level-dependent (BOLD) fMRI has contributions from venous oxygenation and venous cerebral blood volume (CBV) changes. To examine the relative contribution of venous CBV change (ΔCBV(v)) to BOLD fMRI, BOLD and arterial CBV changes (ΔCBV(a)) to a 40-s forepaw stimulation in six α-chloralose anesthetized rats were measured using a magnetization transfer-varied fMRI technique, while total CBV change (ΔCBV(t)) was measured with injection of iron oxide nanoparticles. ΔCBV(v) was obtained by subtracting ΔCBV(a) from ΔCBV(t). We observed a fast ΔCBV(a) response with a time constant of 2.9 ± 2.3s and a slower ΔCBV(v) response with a time constant of 13.5 ± 5.7s and an onset delay of 6.1 ± 3.3s. These results are consistent with earlier studies under different anesthesia and stimulus, supporting that fast CBV(a) and slow CBV(v) responses are generalizable. Assuming the observed post-stimulus BOLD undershoot is at least partly explained by the ΔCBV(v) contribution, the relative contribution of the ΔCBV(v)- and oxygenation-change-related components to the BOLD response was estimated. The relative ΔCBV(v) contribution increases with time during stimulation; whereby it is <0.14 during initial 10s and reaches a maximum possible value of ~0.45 relative to the oxygenation contribution during the 30-40s period after stimulus onset. Our data indicates that the contribution of venous oxygenation change to BOLD fMRI is dominant for short stimulations.

Using manganese-enhanced MRI to understand BOLD

The 1990s were designated "The Decade of the Brain" by U.S. Congress, perhaps in great anticipation of the impact that functional neuroimaging techniques would have on advancing our understanding of how the brain is functionally organized. While it is impossible to overestimate the impact of functional MRI in neuroscience, many aspects of the blood oxygenation level-dependent (BOLD) contrast remain poorly understood, in great part due to the complex relationship between neural activity and hemodynamic changes. To better understand such relationship, it is important to probe neural activity independently. Manganese-enhanced MRI (MEMRI), when used to monitor neural activity, is a technique that uses the divalent manganese ion, Mn(2+), as a surrogate measure of calcium influx. A major advantage of using Mn(2+) as a functional marker is that the contrast obtained is directly related to the accumulation of the ion in excitable cells in an activity dependent manner. As such, the contrast in MEMRI is more directly related to neural activity then hemodynamic-based fMRI techniques. In the present work, the early conceptualization of MEMRI is reviewed, and the comparative experiments that have helped provide a better understanding of the spatial specificity of BOLD signal changes in the cortex is discussed.

Perfusion MR imaging: evolution from initial development to functional studies

A critical indicator of tissue viability and function is blood delivery to the capillary bed (referred to as perfusion or tissue/capillary blood flow), so the measurement of this process has been pursued by many MR scientists. Perfusion MRI is currently an effective tool to non-invasively quantify cerebral blood flow (CBF) and to easily obtain its relative change due to neural activity or other stimulus. This article describes the author's experiences in perfusion MRI over the past quarter-century, including initial development of the field, development of a flow-sensitive alternating inversion recovery (FAIR) MRI technique, development of a functional oxygen consumption MRI measurement approach, validation of the FAIR technique, characterization of perfusion changes induced by neural activity, and determination of arterial blood volume.

Multimodal MRI of experimental stroke

Stroke is the fourth leading cause of death and the leading cause of long-term disability in USA. Brain imaging data from experimental stroke models and stroke patients have shown that there is often a gradual progression of potentially reversible ischemic injury toward infarction. Reestablishing tissue perfusion and/or treating with neuroprotective drugs in a timely fashion are expected to salvage some ischemic tissues. Diffusion-weighted imaging based on magnetic resonance imaging (MRI) in which contrast is based on water motion can detect ischemic injury within minutes after onsets, whereas computed tomography and other imaging modalities fail to detect stroke injury for at least a few hours. Along with quantitative perfusion imaging, the perfusion-diffusion mismatch which approximates the ischemic penumbra could be imaged noninvasively. This review describes recent progresses in the development and application of multimodal MRI and image analysis techniques to study ischemic tissue at risk in experimental stroke in rats.