Layer-dependent BOLD and CBV-weighted fMRI responses in the rat olfactory bulb

Functional MRI of murine olfactory bulbs at 15.2T reveals characteristic activation patters when stimulated by different odors

Thanks to its increased sensitivity, single-shot ultrahigh field functional MRI (UHF fMRI) could lead to valuable insight about subtle brain functions such as olfaction. However, UHF fMRI experiments targeting small organs next to air voids, such as the olfactory bulb, are severely affected by field inhomogeneity problems. Spatiotemporal Encoding (SPEN) is an emerging single-shot MRI technique that could provide a route for bypassing these complications. This is here explored with single-shot fMRI studies on the olfactory bulbs of male and female mice performed at 15.2T. SPEN images collected on these organs at a 108 µm in-plane resolution yielded remarkably large and well-defined responses to olfactory cues. Under suitable T2* weightings these activation-driven changes exceeded 5% of the overall signal intensity, becoming clearly visible in the images without statistical treatment. The nature of the SPEN signal intensity changes in such experiments was unambiguously linked to olfaction, via single-nostril experiments. These experiments highlighted specific activation regions in the external plexiform region and in glomeruli in the lateral part of the bulb, when stimulated by aversive or appetitive odors, respectively. These strong signal activations were non-linear with concentration, and shed light on how chemosensory signals reaching the olfactory epithelium react in response to different cues. Second-level analyses highlighted clear differences among the appetitive, aversive and neutral odor maps; no such differences were evident upon comparing male against female olfactory activation regions.

Odor-evoked layer-specific fMRI activities in the awake mouse olfactory bulb

Awake rodent fMRI is increasingly common over the use of anesthesia since it permits behavioral paradigms and does not confound normal brain function or neurovascular coupling. It is well established that adequate acclimation to the loud fMRI environment and head fixation reduces stress in the rodents and allows for whole brain imaging with little contamination from motion. However, it is unknown whether high-resolution fMRI with increased susceptibility to motion and lower sensitivity can measure small, but spatially discrete, activations in awake mice. To examine this, we used contrast-enhanced cerebral blood volume-weighted (CBVw) fMRI in the mouse olfactory bulb for its enhanced sensitivity and neural specificity. We determined that activation patterns in the glomerular layer to four different odors were spatially distinct and were consistent with previously established histological patterns. In addition, odor-evoked laminar activations were greatest in superficial layers that decreased with laminar depth, similar to previous observations. Interestingly, the fMRI response strengths in the granule cell layer were greater in awake mice than our previous anesthetized rat studies, suggesting that feedback neural activities were intact with wakefulness. We finally determined that fMRI signal changes to repeated odor exposure (i.e., olfactory adaptation) attenuated relatively more in the feedback granule cell layer compared to the input glomerular layer, which is consistent with prior observations. We, therefore, conclude that high-resolution CBVw fMRI can measure odor-specific activation patterns and distinguish changes in laminar activity of head and body restrained awake mice.

Modelling the depth-dependent VASO and BOLD responses in human primary visual cortex

Functional magnetic resonance imaging (fMRI) using a blood-oxygenation-level-dependent (BOLD) contrast is a common method for studying human brain function noninvasively. Gradient-echo (GRE) BOLD is highly sensitive to the blood oxygenation change in blood vessels; however, the spatial signal specificity can be degraded due to signal leakage from activated lower layers to superficial layers in depth-dependent (also called laminar or layer-specific) fMRI. Alternatively, physiological variables such as cerebral blood volume using the VAscular-Space-Occupancy (VASO) contrast have shown higher spatial specificity compared to BOLD. To better understand the physiological mechanisms such as blood volume and oxygenation changes and to interpret the measured depth-dependent responses, models are needed which reflect vascular properties at this scale. For this purpose, we extended and modified the "cortical vascular model" previously developed to predict layer-specific BOLD signal changes in human primary visual cortex to also predict a layer-specific VASO response. To evaluate the model, we compared the predictions with experimental results of simultaneous VASO and BOLD measurements in a group of healthy participants. Fitting the model to our experimental data provided an estimate of CBV change in different vascular compartments upon neural activity. We found that stimulus-evoked CBV change mainly occurs in small arterioles, capillaries, and intracortical arteries and that the contribution from venules and ICVs is smaller. Our results confirm that VASO is less susceptible to large vessel effects compared to BOLD, as blood volume changes in intracortical arteries did not substantially affect the resulting depth-dependent VASO profiles, whereas depth-dependent BOLD profiles showed a bias towards signal contributions from intracortical veins.

A comparative neuroimaging perspective of olfaction and higher-order olfactory processing: on health and disease

Olfactory dysfunction is often the earliest indicator of disease in a range of neurological and psychiatric disorders. One tempting working hypothesis is that pathological changes in the peripheral olfactory system where the body is exposed to many adverse environmental stressors may have a causal role for the brain alteration. Whether and how the peripheral pathology spreads to more central brain regions may be effectively studied in rodent models, and there is successful precedence in experimental models for Parkinson's disease. It is of interest to study whether a similar mechanism may underlie the pathology of psychiatric illnesses, such as schizophrenia. However, direct comparison between rodent models and humans includes challenges under light of comparative neuroanatomy and experimental methodologies used in these two distinct species. We believe that neuroimaging modality that has been the main methodology of human brain studies may be a useful viewpoint to address and fill the knowledge gap between rodents and humans in this scientific question. Accordingly, in the present review article, we focus on brain imaging studies associated with olfaction in healthy humans and patients with neurological and psychiatric disorders, and if available those in rodents. We organize this review article at three levels: 1) olfactory bulb (OB) and peripheral structures of the olfactory system, 2) primary olfactory cortical and subcortical regions, and 3) associated higher-order cortical regions. This research area is still underdeveloped, and we acknowledge that further validation with independent cohorts may be needed for many studies presented here, in particular those with human subjects. Nevertheless, whether and how peripheral olfactory disturbance impacts brain function is becoming even a hotter topic in the ongoing COVID-19 pandemic, given the risk of long-term changes of mental status associated with olfactory infection of SARS-CoV-2. Together, in this review article, we introduce this underdeveloped but important research area focusing on its implications in neurological and psychiatric disorders, with several pioneered publications.

Functional Activities Detected in the Olfactory Bulb and Associated Olfactory Regions in the Human Brain Using T2-Prepared BOLD Functional MRI at 7T

Olfaction is a fundamental sense that plays a vital role in daily life in humans, and can be altered in neuropsychiatric and neurodegenerative diseases. Blood oxygenation level-dependent (BOLD) functional magnetic resonance imaging (fMRI) using conventional echo-planar-imaging (EPI) based sequences can be challenging in brain regions important for olfactory processing, such as the olfactory bulb (OB) and orbitofrontal cortex, mainly due to the signal dropout and distortion artifacts caused by large susceptibility effects from the sinonasal cavity and temporal bone. To date, few studies have demonstrated successful fMRI in the OB in humans. T2-prepared (T2prep) BOLD fMRI is an alternative approach developed especially for performing fMRI in regions affected by large susceptibility artifacts. The purpose of this technical study is to evaluate T2prep BOLD fMRI for olfactory functional experiments in humans. Olfactory fMRI scans were performed on 7T in 14 healthy participants. T2prep BOLD showed greater sensitivity than GRE EPI BOLD in the OB, orbitofrontal cortex and the temporal pole. Functional activation was detected using T2prep BOLD in the OB and associated olfactory regions. Habituation effects and a bi-phasic pattern of fMRI signal changes during olfactory stimulation were observed in all regions. Both positively and negatively activated regions were observed during olfactory stimulation. These signal characteristics are generally consistent with literature and showed a good intra-subject reproducibility comparable to previous human BOLD fMRI studies. In conclusion, the methodology demonstrated in this study holds promise for future olfactory fMRI studies in the OB and other brain regions that suffer from large susceptibility artifacts.

Olfactory Information Storage Engages Subcortical and Cortical Brain Regions That Support Valence Determination

The olfactory bulb (OB) delivers sensory information to the piriform cortex (PC) and other components of the olfactory system. OB-PC synapses have been reported to express short-lasting forms of synaptic plasticity, whereas long-term potentiation (LTP) of the anterior PC (aPC) occurs predominantly by activating inputs from the prefrontal cortex. This suggests that brain regions outside the olfactory system may contribute to olfactory information processing and storage. Here, we compared functional magnetic resonance imaging BOLD responses triggered during 20 or 100 Hz stimulation of the OB. We detected BOLD signal increases in the anterior olfactory nucleus (AON), PC and entorhinal cortex, nucleus accumbens, dorsal striatum, ventral diagonal band of Broca, prelimbic–infralimbic cortex (PrL-IL), dorsal medial prefrontal cortex, and basolateral amygdala. Significantly stronger BOLD responses occurred in the PrL-IL, PC, and AON during 100 Hz compared with 20 Hz OB stimulation. LTP in the aPC was concomitantly induced by 100 Hz stimulation. Furthermore, 100 Hz stimulation triggered significant nuclear immediate early gene expression in aPC, AON, and PrL-IL. The involvement of the PrL-IL in this process is consistent with its putative involvement in modulating behavioral responses to odor experience. Furthermore, these results indicate that OB-mediated information storage by the aPC is embedded in a connectome that supports valence evaluation.

Postsynaptic activity of inhibitory neurons evokes hemodynamic fMRI responses

Functional MRI responses are localized to the synaptic sites of evoked inhibitory neurons, but it is unknown whether, or by what mechanisms, these neurons initiate functional hyperemia. Here, the neuronal origins of these hemodynamic responses were investigated by fMRI or local field potential and blood flow measurements during topical application of pharmacological agents when GABAergic granule cells in the rat olfactory bulb were synaptically targeted. First, to examine if postsynaptic activation of these inhibitory neurons was required for neurovascular coupling, we applied an NMDA receptor antagonist during cerebral blood volume-weighted fMRI acquisition and found that responses below the drug application site (up to ~1.5 mm) significantly decreased within ~30 min. Similarly, large decreases in granule cell postsynaptic activities and blood flow responses were observed when AMPA or NMDA receptor antagonists were applied. Second, inhibition of nitric oxide synthase preferentially decreased the initial, fast component of the blood flow response, while inhibitors of astrocyte-specific glutamate transporters and vasoactive intestinal peptide receptors did not decrease blood flow responses. Third, inhibition of GABA release with a presynaptic GABAB receptor agonist caused less reduction of neuronal and blood flow responses compared to the postsynaptic glutamate receptor antagonists. In conclusion, local hyperemia by synaptically-evoked inhibitory neurons was primarily driven by their postsynaptic activities, possibly through NMDA receptor-dependent calcium signaling that was not wholly dependent on nitric oxide.

fMRI study of olfactory processing in mice under three anesthesia protocols: Insight into the effect of ketamine on olfactory processing

Functional magnetic resonance imaging (fMRI) is a valuable tool for studying neural activations in the central nervous system of animals due to its wide spatial coverage and non-invasive nature. However, the advantages of fMRI have not been fully realized in functional studies in mice, especially in the olfactory system, possibly due to the lack of suitable anesthesia protocols with spontaneous breathing. Since mice are widely used in biomedical research, it is desirable to evaluate different anesthesia protocols for olfactory fMRI studies in mice. Dexmedetomidine (DEX) as a sedative/anesthetic has been introduced to fMRI studies in mice, but it has a limited anesthesia duration. To extend the anesthesia duration, DEX has been combined with a low dose of isoflurane (ISO) or ketamine (KET) in previous functional studies in mice. In this report, olfactory fMRI studies were performed under three anesthesia protocols (DEX alone, DEX/ISO, and DEX/KET) in three different groups of mice. Isoamyl-acetate was used as an odorant, and the odorant-induced neural activations were measured by blood oxygenation-level dependent (BOLD) fMRI. BOLD fMRI responses were observed in the olfactory bulb (OB), anterior olfactory nuclei (AON), and piriform cortex (Pir). Interestingly, BOLD fMRI activations were also observed in the prefrontal cortical region (PFC), which are most likely caused by the draining vein effect. The response in the OB showed no adaptation to either repeated odor stimulations or continuous odor exposure, but the response in the Pir showed adaptation during the continuous odor exposure. The data also shows that ISO suppresses the olfactory response in the OB and AON, while KET enhances the olfactory response in the Pir. Thus, DEX/KET should be an attractive anesthesia for olfactory fMRI in mice.

A methodological investigation of a flexible surface MRI coil to obtain functional signals from the human olfactory bulb

BACKGROUND

Mammalian olfaction begins with transduction in olfactory receptors, continues with extensive processing in the olfactory bulb, and culminates in cortical representation. Most rodent studies on the functional neuroanatomy of olfaction have concentrated on the olfactory bulb, yet whether this structure is tuned only to basic chemical features of odorants or also to higher-order perceptual features is unclear.

NEW METHOD

Whereas studies of the human brain can typically uncover involvement of higher-order feature extraction, this has not been possible in the case of the olfactory bulb, inaccessible to fMRI. The present study examined whether a novel method of acquisition using a facial coil could overcome this limitation.

RESULTS

A series of experiments provided preliminary evidence of odor-driven responses in the human olfactory bulb, and found that these responses differed between individuals.

COMPARISON WITH EXISTING METHODS AND CONCLUSIONS

The present preliminary technical achievement renders possible to design novel human odor fMRI studies by considering the olfactory system from the olfactory bulb to associative areas.

Functional MRI of the mouse olfactory system

Although olfactory dysfunction is an early warning sign of Alzheimer's and Parkinson's diseases, and is commonly present in a range of other neurodegenerative disorders, the mechanisms for its pathogenesis are not yet clear. Since fMRI allows the mapping of spatial and temporal patterns of activity in multiple brain regions simultaneously, it serves as a powerful tool to study olfactory dysfunction in animal models of neurodegenerative diseases. Nonetheless, there have been no reports to date of mapping odor-induced activation patterns beyond the olfactory bulb to the extended networks of olfactory and limbic archicortex, likely due to the small size of the mouse brain. Therefore, using an 11.7 T magnet and a blood volume-weighted fMRI technique, we mapped the functional neuroanatomy of the mouse olfactory system. Consistent with reports on imaging of the much larger human brain, we mapped activity in regions of the olfactory bulb, as well as olfactory and limbic archicortex. By using two distinct odorants, we further demonstrated odorant-specific activation patterns. Our work thus provides a methodological framework for fMRI studies of olfactory dysfunction in mouse models of neurodegeneration.

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.

Mesoscopic and microscopic imaging of sensory responses in the same animal

Imaging based on blood flow dynamics is widely used to study sensory processing. Here we investigated the extent to which local neuronal and capillary responses (two-photon microscopy) are correlated to mesoscopic responses detected with fast ultrasound (fUS) and BOLD-fMRI. Using a specialized chronic olfactory bulb preparation, we report that sequential imaging of the same mouse allows quantitative comparison of odour responses, imaged at both microscopic and mesoscopic scales. Under these conditions, functional hyperaemia occurred at the threshold of neuronal activation and fUS-CBV signals could be detected at the level of single voxels with activation maps varying according to blood velocity. Both neuronal and vascular responses increase non-linearly as a function of odour concentration, whereas both microscopic and mesoscopic vascular responses are linearly correlated to local neuronal calcium. These data establish strengths and limits of mesoscopic imaging techniques to report neural activity.

Neuronal activity leads to a local increase in blood flow and volume, a process termed hyperaemia. Here, the authors employ multiple imaging approaches of neuronal and vascular activity at varying resolution to delineate the spatiotemporal dynamics of neurovascular coupling evoked by odours in the olfactory bulb.

Effects of nicotine on odor-induced increases in regional blood flow in the olfactory bulb in rats

The present study examined the effects of nicotinic acetylcholine receptor activation on the odor-induced blood flow response in the olfactory bulb. In urethane-anesthetized rats, odor stimulation (5% amyl acetate, 30 s) produced an increase in olfactory bulb blood flow (reaching 107% ± 3% of the pre-stimulus basal values), without changes in frontal cortical blood flow or mean arterial pressure. Intravenous injection of nicotine (30 μg/kg), a nicotinic acetylcholine receptor agonist, significantly augmented the odor-induced increase response of olfactory bulb blood flow, without changes in the basal blood flow level. The nicotine-induced augmentation of the olfactory bulb blood flow response to odor was negated by dihydro-β-erythroidine, an α4β2-preferring nicotinic acetylcholine receptor antagonist. Our results suggest that the activation of α4β2-like neuronal nicotinic acetylcholine receptors in the brain potentiates an odor-induced blood flow response in the olfactory bulb.

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.

Brain Iron Distribution after Multiple Doses of Ultra-small Superparamagnetic Iron Oxide Particles in Rats

The purpose of this study is to determine the effects of high cumulative doses of ultra-small paramagnetic iron oxide (USPIO) used in neuroimaging studies. We intravenously administered 8 mg/kg of 2 USPIO compounds daily for 4 wk to male Sprague-Dawley rats (Crl:SD). Multiecho gradient-echo MRI, serum iron levels, and histology were performed at the end of dosing and after a 7-d washout period. R2* maps and quantitative susceptibility maps (QSM) were generated from multiecho gradient-echo data. R2* maps and QSM showed iron accumulation in brain ventricles on MR images acquired at the 4- and 5-wk time points. Estimates from QSM data showed ventricular iron concentration was equal to or higher than serum iron concentration. Histologic analysis revealed choroid plexus hemosiderosis and midbrain vacuolation, without iron deposition in brain parenchyma. Serum iron levels increased with administration of both compounds, and a 7-d washout period effectively reduced serum iron levels of one but not both of the compounds. High cumulative doses from multiple, frequent administrations of USPIO can lead to iron deposition in brain ventricles, resulting in persistent signal loss on T2*-weighted images. Techniques such as QSM are helpful in quantifying iron biodistribution in this situation.

Activity Patterns Elicited by Airflow in the Olfactory Bulb and Their Possible Functions

Olfactory sensory neurons (OSNs) can sense both odorants and airflows. In the olfactory bulb (OB), the coding of odor information has been well studied, but the coding of mechanical stimulation is rarely investigated. Unlike odor-sensing functions of OSNs, the airflow-sensing functions of OSNs are also largely unknown. Here, the activity patterns elicited by mechanical airflow in male rat OBs were mapped using fMRI and correlated with local field potential recordings. In an attempt to reveal possible functions of airflow sensing, the relationship between airflow patterns and physiological parameters was also examined. We found the following: (1) the activity pattern in the OB evoked by airflow in the nasal cavity was more broadly distributed than patterns evoked by odors; (2) the pattern intensity increases with total airflow, while the pattern topography with total airflow remains almost unchanged; and (3) the heart rate, spontaneous respiratory rate, and electroencephalograph power in the β band decreased with regular mechanical airflow in the nasal cavity. The mapping results provide evidence that the signals elicited by mechanical airflow in OSNs are transmitted to the OB, and that the OB has the potential to code and process mechanical information. Our functional data indicate that airflow rhythm in the olfactory system can regulate the physiological and brain states, providing an explanation for the effects of breath control in meditation, yoga, and Taoism practices.SIGNIFICANCE STATEMENT Presentation of odor information in the olfactory bulb has been well studied, but studies about breathing features are rare. Here, using blood oxygen level-dependent functional MRI for the first time in such an investigation, we explored the global activity patterns in the rat olfactory bulb elicited by airflow in the nasal cavity. We found that the activity pattern elicited by airflow is broadly distributed, with increasing pattern intensity and similar topography under increasing total airflow. Further, heart rate, spontaneous respiratory rate in the lung, and electroencephalograph power in the β band decreased with regular airflow in the nasal cavity. Our study provides further understanding of the airflow map in the olfactory bulb in vivo, and evidence for the possible mechanosensitivity functions of olfactory sensory neurons.

Olfactory function in an excitotoxic model for secondary neuronal degeneration: Role of dopaminergic interneurons

Secondary neuronal degeneration (SND) occurring in Traumatic brain injury (TBI) consists in downstream destructive events affecting cells that were not or only marginally affected by the initial wound, further increasing the effects of the primary injury. Glutamate excitotoxicity is hypothesized to play an important role in SND. TBI is a common cause of olfactory dysfunction that may be spontaneous and partially recovered. The role of the glutamate excitotoxicity in the TBI-induced olfactory dysfunction is still unknown. We investigated the effects of excitotoxicity induced by bilateral N-Methyl-D-Aspartate (NMDA) OB administration in the olfactory function, OB volumes, and subventricular zone (SVZ) and OB neurogenesis in rats. NMDA OB administration induced a decrease in the number of correct choices in the olfactory discrimination tests one week after lesions (p<0.01), and a spontaneous recovery of the olfactory deficit two weeks after lesions (p<0.05). A lack of correlation between OB volumes and olfactory function was observed. An increase in SVZ neurogenesis (Ki67+ cells, PSANCAM+ cells (p<0.01) associated with an increase in OB glomerular dopaminergic immunostaining (p<0.05) were related to olfactory function recovery. The present results show that changes in OB volumes cannot explain the recovery of the olfactory function and suggest a relevant role for dopaminergic OB interneurons in the pathophysiology of recovery of loss of smell in TBI.

Dominance of layer-specific microvessel dilation in contrast-enhanced high-resolution fMRI: Comparison between hemodynamic spread and vascular architecture with CLARITY

Contrast-enhanced cerebral blood volume-weighted (CBVw) fMRI response peaks are specific to the layer of evoked synaptic activity (Poplawsky et al., 2015), but the spatial resolution limit of CBVw fMRI is unknown. In this study, we measured the laminar spread of the CBVw fMRI evoked response in the external plexiform layer (EPL, 265 ± 65 μm anatomical thickness, mean ± SD, n = 30 locations from 5 rats) of the rat olfactory bulb during electrical stimulation of the lateral olfactory tract and examined its potential vascular source. First, we obtained the evoked CBVw fMRI responses with a 55 × 55 μm² in-plane resolution and a 500-μm thickness at 9.4 T, and found that the fMRI signal peaked predominantly in the inner half of EPL (136 ± 54 μm anatomical thickness). The mean full-width at half-maximum of these fMRI peaks was 347 ± 102 μm and the functional spread was approximately 100 or 200 μm when the effects of the laminar thicknesses of EPL or inner EPL were removed, respectively. Second, we visualized the vascular architecture of EPL from a different rat using a Clear Lipid-exchanged Anatomically Rigid Imaging/immunostaining-compatible Tissue hYdrogel (CLARITY)-based tissue preparation method and confocal microscopy. Microvascular segments with an outer diameter of <11 μm accounted for 64.3% of the total vascular volume within EPL and had a mean segment length of 55 ± 40 μm (n = 472). Additionally, vessels that crossed the EPL border had a mean segment length outside of EPL equal to 73 ± 61 μm (n = 28), which is comparable to half of the functional spread (50-100 μm). Therefore, we conclude that dilation of these microvessels, including capillaries, likely dominate the CBVw fMRI response and that the biological limit of the fMRI spatial resolution is approximately the average length of 1-2 microvessel segments, which may be sufficient for examining sublaminar circuits.

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.

Light controls cerebral blood flow in naive animals

Optogenetics is increasingly used to map brain activation using techniques that rely on functional hyperaemia, such as opto-fMRI. Here we test whether light stimulation protocols similar to those commonly used in opto-fMRI or to study neurovascular coupling modulate blood flow in mice that do not express light sensitive proteins. Combining two-photon laser scanning microscopy and ultrafast functional ultrasound imaging, we report that in the naive mouse brain, light per se causes a calcium decrease in arteriolar smooth muscle cells, leading to pronounced vasodilation, without excitation of neurons and astrocytes. This photodilation is reversible, reproducible and energy-dependent, appearing at about 0.5 mJ. These results impose careful consideration on the use of photo-activation in studies involving blood flow regulation, as well as in studies requiring prolonged and repetitive stimulations to correct cellular defects in pathological models. They also suggest that light could be used to locally increase blood flow in a controlled fashion.

Combination of optogenetics and BOLD fMRI is routinely used to map neuronal activity upon photostimulation. Here the authors show that light, shone at intensities used in optogenetic studies, dilates vessels and increases blood flow independently of exogenous light-sensitive proteins in the mouse brain.

Improved spatial accuracy of functional maps in the rat olfactory bulb using supervised machine learning approach

Functional MRI (fMRI) is a popular and important tool for noninvasive mapping of neural activity. As fMRI measures the hemodynamic response, the resulting activation maps do not perfectly reflect the underlying neural activity. The purpose of this work was to design a data-driven model to improve the spatial accuracy of fMRI maps in the rat olfactory bulb. This system is an ideal choice for this investigation since the bulb circuit is well characterized, allowing for an accurate definition of activity patterns in order to train the model. We generated models for both cerebral blood volume weighted (CBVw) and blood oxygen level dependent (BOLD) fMRI data. The results indicate that the spatial accuracy of the activation maps is either significantly improved or at worst not significantly different when using the learned models compared to a conventional general linear model approach, particularly for BOLD images and activity patterns involving deep layers of the bulb. Furthermore, the activation maps computed by CBVw and BOLD data show increased agreement when using the learned models, lending more confidence to their accuracy. The models presented here could have an immediate impact on studies of the olfactory bulb, but perhaps more importantly, demonstrate the potential for similar flexible, data-driven models to improve the quality of activation maps calculated using fMRI data.

Exploitation of temporal redundancy in compressed sensing reconstruction of fMRI studies with a prior-based algorithm (PICCS)

PURPOSE

Compressed sensing is a technique used to accelerate magnetic resonance imaging (MRI) acquisition without compromising image quality. While it has proven particularly useful in dynamic imaging procedures such as cardiac cine, very few authors have applied it to functional magnetic resonance imaging (fMRI). The purpose of the present study was to check whether the prior image constrained compressed sensing (PICCS) algorithm, which is based on an available prior image, can improve the statistical maps in fMRI better than other strategies that also exploit temporal redundancy.

METHODS

PICCS was compared to spatiotemporal total variation (TTV) and k-t FASTER, since they have already demonstrated high performance and robustness in other MRI applications, such as cardiac cine MRI and resting state fMRI, respectively. The prior image for PICCS was the average of all undersampled data. Both PICCS and TTV were solved using the split Bregman formulation. K-t FASTER algorithm relies on matrix completion to reconstruct the undersampled k-spaces. The three algorithms were evaluated using two datasets with high and low signal-to-noise ratio (SNR)-BOLD contrast-acquired in a 7 T preclinical MRI scanner and retrospectively undersampled at various rates (i.e., acceleration factors). The authors evaluated their performance in terms of the sensitivity/specificity of BOLD detection through receiver operating characteristic curves and by visual inspection of the statistical maps.

RESULTS

With high SNR studies, PICCS performed similarly to the state-of-the-art algorithms TTV and k-t FASTER and provided consistent BOLD signal at the ROI. In scenarios with low SNR and high acceleration factors, PICCS still provided consistent maps and higher sensitivity/specificity than TTV, whereas k-t FASTER failed to provide significant maps.

CONCLUSIONS

The authors performed a comparison between three reconstructions (PICCS, TTV, and k-t FASTER) that exploit temporal redundancy in fMRI. The prior-based algorithm, PICCS, preserved BOLD activation and sensitivity/specificity better than TTV and k-t FASTER in noisy scenarios. The PICCS algorithm can potentially reach an acceleration factor of ×8 and still provide BOLD contrast in the ROI with an area under the curve over 0.99.

Layer-Specific fMRI Responses to Excitatory and Inhibitory Neuronal Activities in the Olfactory Bulb

High-resolution functional magnetic resonance imaging (fMRI) detects localized neuronal activity via the hemodynamic response, but it is unclear whether it accurately identifies neuronal activity specific to individual layers. To address this issue, we preferentially evoked neuronal activity in superficial, middle, and deep layers of the rat olfactory bulb: the glomerular layer by odor (5% amyl acetate), the external plexiform layer by electrical stimulation of the lateral olfactory tract (LOT), and the granule cell layer by electrical stimulation of the anterior commissure (AC), respectively. Electrophysiology, laser-Doppler flowmetry of cerebral blood flow (CBF), and blood oxygenation level-dependent (BOLD) and cerebral blood volume-weighted (CBV) fMRI at 9.4 T were performed independently. We found that excitation of inhibitory granule cells by stimulating LOT and AC decreased the spontaneous multi-unit activities of excitatory mitral cells and subsequently increased CBF, CBV, and BOLD signals. Odor stimulation also increased the hemodynamic responses. Furthermore, the greatest CBV fMRI responses were discretely separated into the same layers as the evoked neuronal activities for all three stimuli, whereas BOLD was poorly localized with some exception to the poststimulus undershoot. In addition, the temporal dynamics of the fMRI responses varied depending on the stimulation pathway, even within the same layer. These results indicate that the vasculature is regulated within individual layers and CBV fMRI has a higher fidelity to the evoked neuronal activity compared with BOLD. Our findings are significant for understanding the neuronal origin and spatial specificity of hemodynamic responses, especially for the interpretation of laminar-resolution fMRI.

SIGNIFICANCE STATEMENT

Functional magnetic resonance imaging (fMRI) is a noninvasive, in vivo technique widely used to map function of the entire brain, including deep structures, in animals and humans. However, it measures neuronal activity indirectly by way of the vascular response. It is currently unclear how finely the hemodynamic response is regulated within single cortical layers and whether increased inhibitory neuronal activities affect fMRI signal changes. Both laminar specificity and the neural origins of fMRI are important to interpret functional maps properly, which we investigated by activating discrete rat olfactory bulb circuits.

Advances of Molecular Imaging for Monitoring the Anatomical and Functional Architecture of the Olfactory System

The olfactory system of organisms serves as a genetically and anatomically model for studying how sensory input can be translated into behavior output. Some neurologic diseases are considered to be related to olfactory disturbance, especially Alzheimer's disease, Parkinson's disease, multiple sclerosis, and so forth. However, it is still unclear how the olfactory system affects disease generation processes and olfaction delivery processes. Molecular imaging, a modern multidisciplinary technology, can provide valid tools for the early detection and characterization of diseases, evaluation of treatment, and study of biological processes in living subjects, since molecular imaging applies specific molecular probes as a novel approach to produce special data to study biological processes in cellular and subcellular levels. Recently, molecular imaging plays a key role in studying the activation of olfactory system, thus it could help to prevent or delay some diseases. Herein, we present a comprehensive review on the research progress of the imaging probes for visualizing olfactory system, which is classified on different imaging modalities, including PET, MRI, and optical imaging. Additionally, the probes' design, sensing mechanism, and biological application are discussed. Finally, we provide an outlook for future studies in this field.

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

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

fMRI study of olfaction in the olfactory bulb and high olfactory structures of rats: Insight into their roles in habituation

Cerebral blood volume (CBV) fMRI with ultrasmall superparamagnetic iron oxide particles (USPIO) as a contrast agent was used to investigate olfactory processing in rats. fMRI data were acquired in sixteen 0.75-mm coronal slices covering the olfactory bulb (OB) and higher olfactory regions (HOR), including the anterior olfactory nucleus and piriform cortex. For each animal, multiple consecutive fMRI measurements were made during a 3-h experiment session, with each measurement consisting of a baseline period, an odorant stimulation period, and a recovery period. Two different stimulation paradigms with a stimulation period of 40s or 80s, respectively, were used to study olfactory processing. Odorant-induced CBV increases were robustly observed in the OB and HOR of each individual animal. Olfactory adaptation, which is characterized by an attenuation of responses to continuous exposure or repeated stimulations, has different characteristics in the OB and HOR. For adaptation to repeated stimuli, while it was observed in both the OB and HOR, CBV responses in the HOR were attenuated more significantly than responses in the OB. In contrast, within each continuous 40-s or 80-s odor exposure, CBV responses in the OB were stable and did not show adaptation, but the CBV responses in the HOR were state dependent, with no adaptation during initial exposures, but significant adaptation during following exposures. These results support previous reports that HOR plays a more significant role than OB in olfactory habituation. The technical approach presented in this study should enable more extensive fMRI studies of olfactory processing in rats.

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

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

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

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

Compressed sensing fMRI using gradient-recalled echo and EPI sequences

Compressed sensing (CS) may be useful for accelerating data acquisitions in high-resolution fMRI. However, due to the inherent slow temporal dynamics of the hemodynamic signals and concerns of potential statistical power loss, the CS approach for fMRI (CS-fMRI) has not been extensively investigated. To evaluate the utility of CS in fMRI application, we systematically investigated the properties of CS-fMRI using computer simulations and in vivo experiments of rat forepaw sensory and odor stimulations with gradient-recalled echo (GRE) and echo planar imaging (EPI) sequences. Various undersampling patterns along the phase-encoding direction were studied and k-t FOCUSS was used as the CS reconstruction algorithm, which exploits the temporal redundancy of images. Functional sensitivity, specificity, and time courses were compared between fully-sampled and CS-fMRI with reduction factors of 2 and 4. CS-fMRI with GRE, but not with EPI, improves the statistical sensitivity for activation detection over the fully sampled data when the ratio of the fMRI signal change to noise is low. CS improves the temporal resolution and reduces temporal noise correlations. While CS reduces the functional response amplitudes, the noise variance is also reduced to make the overall activation detection more sensitive. Consequently, CS is a valuable fMRI acceleration approach, especially for GRE fMRI studies.