Low-frequency calcium oscillations accompany deoxyhemoglobin oscillations in rat somatosensory cortex

Neurovascular effects of cocaine: relevance to addiction

Cocaine is a highly addictive drug, and its use is associated with adverse medical consequences such as cerebrovascular accidents that result in debilitating neurological complications. Indeed, brain imaging studies have reported severe reductions in cerebral blood flow (CBF) in cocaine misusers when compared to the brains of healthy non-drug using controls. Such CBF deficits are likely to disrupt neuro-vascular interaction and contribute to changes in brain function. This review aims to provide an overview of cocaine-induced CBF changes and its implication to brain function and to cocaine addiction, including its effects on tissue metabolism and neuronal activity. Finally, we discuss implications for future research, including targeted pharmacological interventions and neuromodulation to limit cocaine use and mitigate the negative impacts.

Spatiotemporal dedifferentiation of the global brain signal topography along the adult lifespan

Age-related variations in many regions and/or networks of the human brain have been uncovered using resting-state functional magnetic resonance imaging. However, these findings did not account for the dynamical effect the brain's global activity (global signal [GS]) causes on local characteristics, which is measured by GS topography. To address this gap, we tested GS topography including its correlation with age using a large-scale cross-sectional adult lifespan dataset (n = 492). Both GS topography and its variation with age showed frequency-specific patterns, reflecting the spatiotemporal characteristics of the dynamic change of GS topography with age. A general trend toward dedifferentiation of GS topography with age was observed in both spatial (i.e., less differences of GS between different regions) and temporal (i.e., less differences of GS between different frequencies) dimensions. Further, methodological control analyses suggested that although most age-related dedifferentiation effects remained across different preprocessing strategies, some were triggered by neuro-vascular coupling and physiological noises. Together, these results provide the first evidence for age-related effects on global brain activity and its topographic-dynamic representation in terms of spatiotemporal dedifferentiation.

The decoupling between hemodynamic parameters and neural activity implies a complex origin of spontaneous brain oscillations

Introduction

Spontaneous low-frequency oscillations play a key role in brain activity. However, the underlying mechanism and origin of low-frequency oscillations remain under debate.

Methods

Optical imaging and an electrophysiological recording system were combined to investigate spontaneous oscillations in the hemodynamic parameters and neuronal activity of awake and anesthetized mice after Nω-nitro-L-arginine methyl ester (L-NAME) administration.

Results

The spectrum of local field potential (LFP) signals was significantly changed by L-NAME, which was further corroborated by the increase in energy and spatial synchronization. The important finding was that L-NAME triggered regular oscillations in both LFP signals and hemodynamic signals. Notably, the frequency peak of hemodynamic signals can be different from that of LFP oscillations in awake mice.

Discussion

A model of the neurovascular system was proposed to interpret this mismatch of peak frequencies, supporting the view that spontaneous low-frequency oscillations arise from multiple sources.

Spatiotemporal analysis of blood plasma and blood cell flow fluctuations of cerebral microcirculation in anesthetized rats

Cerebral hemodynamics fluctuates spontaneously over broad frequency ranges. However, its spatiotemporal coherence of flow oscillations in cerebral microcirculation remains incompletely understood. The objective of this study was to characterize the spatiotemporal fluctuations of red blood cells (RBCs) and plasma flow in the rat cerebral microcirculation by simultaneously imaging their dynamic behaviors. Comparisons of changes in cross-section diameters between RBC and plasma flow showed dissociations in penetrating arterioles. The results indicate that vasomotion has the least effect on the lateral movement of circulating RBCs, resulting in variable changes in plasma layer thickness. Parenchymal capillaries exhibited slow fluctuations in RBC velocity (0.1 to 0.3 Hz), regardless of capillary diameter fluctuations (<0.1 Hz). Temporal fluctuations and the velocity of RBCs decreased significantly at divergent capillary bifurcations. The results indicate that a transit of RBCs generates flow resistance in the capillaries and that slow velocity fluctuations of the RBCs are subject to a number of bifurcations. In conclusion, the high-frequency oscillation of the blood flow is filtered at the bifurcation through the capillary networks. Therefore, a number of bifurcations in the cerebral microcirculation may contribute to the power of low-frequency oscillations.

Frequency-dependent genetic modulation of neuronal oscillations: a combined transcriptome and resting-state functional MRI study

Neuronal oscillations within certain frequency bands are assumed to associate with specific neural processes and cognitive functions. To examine this hypothesis, transcriptome-neuroimaging spatial correlation analysis was applied to resting-state functional magnetic resonance imaging data from 793 healthy individuals and gene expression data from the Allen Human Brain Atlas. We found that expression measures of 336 genes were correlated with fractional amplitude of low-frequency fluctuations (fALFF) in the slow-4 band (0.027-0.073 Hz), whereas there were no expression-fALFF correlations for the other frequency bands. Furthermore, functional enrichment analyses showed that these slow-4 fALFF-related genes were mainly enriched for ion channel, synaptic function, and neuronal system as well as many neuropsychiatric disorders. Specific expression analyses demonstrated that these genes were specifically expressed in brain tissue, in neurons, and during the late stage of cortical development. Concurrently, the fALFF-related genes were linked to multiple behavioral domains, including dementia, attention, and emotion. In addition, these genes could construct a protein-protein interaction network supported by 30 hub genes. Our findings of a frequency-dependent genetic modulation of spontaneous neuronal activity may support the concept that neuronal oscillations within different frequency bands capture distinct neurobiological processes from the perspective of underlying molecular mechanisms.

Time-frequency cross-coupling between cortical low-frequency neuronal calcium oscillations and blood oxygen metabolism evoked by ultrasound stimulation

Low-intensity transcranial ultrasound stimulation (TUS) can modulate the coupling of high-frequency (160-200 Hz) neural oscillations and cerebral blood oxygen metabolism (BOM); however, the correlation of low-frequency (0-2 Hz) neural oscillations with BOM in temporal and frequency domains under TUS remains unclear. To address this, we monitored the TUS-evoked neuronal calcium oscillations and BOM simultaneously in the mouse visual cortex by using multimodal optical imaging with a high spatiotemporal resolution. We demonstrated that TUS can significantly increase the intensity of the neuronal calcium oscillations and BOM; the peak value, peak time, and duration of calcium oscillations are functionally related to stimulation duration; TUS does not significantly increase the neurovascular coupling strength between calcium oscillations and BOM in the temporal domain; the time differences of the energy peaks between TUS-induced calcium oscillations and BOM depend on their spectral ranges; the frequency differences of the energy peaks between TUS-induced calcium oscillations and BOM depend on their time ranges; and TUS can significantly change the phase of calcium oscillations and BOM from uniform distribution to a more concentrated region. In conclusion, ultrasound stimulation can evoke the time-frequency cross-coupling between the cortical low-frequency neuronal calcium oscillations and BOM in mouse.

Cocaine's cerebrovascular vasoconstriction is associated with astrocytic Ca2+ increase in mice

Human and animal studies have reported widespread reductions in cerebral blood flow associated with chronic cocaine exposures. However, the molecular and cellular mechanisms underlying cerebral blood flow reductions are not well understood. Here, by combining a multimodal imaging platform with a genetically encoded calcium indicator, we simultaneously measured the effects of acute cocaine on neuronal and astrocytic activity, tissue oxygenation, hemodynamics and vascular diameter changes in the mouse cerebral cortex. Our results showed that cocaine constricted blood vessels (measured by vessel diameter Φ changes), decreasing cerebral total blood volume (HbT) and temporally reducing tissue oxygenation. Cellular imaging showed that the mean astrocytic Ca2+ dependent fluorescence [Formula: see text] increase in response to cocaine was weaker but longer lasting than the mean neuronal Ca2+ dependent fluorescence [Formula: see text] changes. Interestingly, while cocaine-induced [Formula: see text] increase was temporally correlated with tissue oxygenation change, the [Formula: see text] elevation after cocaine was in temporal correspondence with the long-lasting decrease in arterial blood volumes. To determine whether the temporal association between astrocytic activation and cocaine induced vasoconstriction reflected a causal association we inhibited astrocytic Ca2+ using GFAP-DREADD(Gi). Inhibition of astrocytes attenuated the vasoconstriction resulting from cocaine, providing evidence that astrocytes play a critical role in cocaine's vasoconstrictive effects in the brain. These results indicate that neurons and astrocytes play different roles in mediating neurovascular coupling in response to cocaine. Our findings implicate neuronal activation as the main driver of the short-lasting reduction in tissue oxygenation and astrocyte long-lasting activation as the driver of the persistent vasoconstriction with cocaine. Understanding the cellular and vascular interaction induced by cocaine will be helpful for future putative treatments to reduce cerebrovascular pathology from cocaine use.

Optical imaging of stimulation-evoked cortical activity using GCaMP6f and jRGECO1a

BACKGROUND

Genetically encoded calcium indicators (GECIs), especially the GCaMP-based green fluorescence GECIs have been widely used for in vivo detection of neuronal activity in rodents by measuring intracellular neuronal Ca²⁺ changes. More recently, jRGECO1a, a red shifted GECI, has been reported to detect neuronal Ca²⁺ activation. This opens the possibility of using dual-color GECIs for simultaneous interrogation of different cell populations. However, there has been no report to compare the functional difference between these two GECIs for in vivo imaging. Here, a comparative study is reported on neuronal responses to sensory stimulation using GCaMP6f and jRGECO1a that were virally delivered into the neurons in the somatosensory cortex of two different groups of animals, respectively.

METHODS

GCaMP6f and jRGECO1a GECI were virally delivered to sensory cortex. After 3-4 weeks, the animals were imaged to capture the spatiotemporal changes of neuronal Ca²⁺ and the hemodynamic responses to forepaw electrical stimulation (0.3 mA, 0.3 ms/pulse, 0.03 Hz). The stimulation-evoked neuronal Ca²⁺ transients expressed with GCaMP6f or jRGECO1a were recorded during the baseline period and after an acute cocaine administration (1 mg/kg, i.v.).

RESULTS

Histology confirmed that the efficiency of jRGECO1a and GCaMP6f expression into the cortical neurons was similar, i.e., 34%±3% and 32.7%±1.6%, respectively. Our imaging in vivo showed that the hemodynamic responses to the stimulation were the same between jRGECO1a and GCaMP6f expressed groups. Although the stimulation-evoked fluorescence change (∆F/F) and the time-to-peak of the neuronal Ca²⁺ transients were not significantly different between these two indicators, the full-width-half-maximum (FWHM) duration of the ∆F/F rise in the jRGECO1a-expressed group (0.16±0.02 s) was ~50 ms or 46% longer than that of the GCaMP6f group (0.11±0.003 s), indicating a longer recovery time in jRGECO1a than in GCaMP6f transients (P<0.01). This is likely due to the longer off rate of jRGECO1a than that of GCaMP6f. After cocaine, the time-to-peak of Ca²⁺ transients was delayed and their FWHM duration was prolonged for both expression groups, indicating that these are cocaine's effects on neuronal Ca²⁺ signaling and not artifacts due to the property differences of the GCEIs.

CONCLUSIONS

This study shows that both jRGECO1a and GCaMP6f have sufficient sensitivity for tracking single-stimulation-evoked Ca²⁺ transients to detect neuronal activities from the brain. Since these GECIs are emitted at the different wavelengths, it will be possible to use them together to characterize the activity of different cell types (e.g., neurons and astrocytes) to study brain activation and brain functional changes in normal or diseased brains.

The rostro-caudal gradient in the prefrontal cortex and its modulation by subthalamic deep brain stimulation in Parkinson's disease

Deep brain stimulation of the subthalamic nucleus (STN-DBS) alleviates motor symptoms in Parkinson’s disease (PD) but also affects the prefrontal cortex (PFC), potentially leading to cognitive side effects. The present study tested alterations within the rostro-caudal hierarchy of neural processing in the PFC induced by STN-DBS in PD. Granger-causality analyses of fast functional near-infrared spectroscopy (fNIRS) measurements were used to infer directed functional connectivity from intrinsic PFC activity in 24 PD patients treated with STN-DBS. Functional connectivity was assessed ON stimulation, in steady-state OFF stimulation and immediately after the stimulator was switched ON again. Results revealed that STN-DBS significantly enhanced the rostro-caudal hierarchical organization of the PFC in patients who had undergone implantation early in the course of the disease, whereas it attenuated the rostro-caudal hierarchy in late-implanted patients. Most crucially, this systematic network effect of STN-DBS was reproducible in the second ON stimulation measurement. Supplemental analyses demonstrated the significance of prefrontal networks for cognitive functions in patients and matched healthy controls. These findings show that the modulation of prefrontal functional networks by STN-DBS is dependent on the disease duration before DBS implantation and suggest a neurophysiological mechanism underlying the side effects on prefrontally-guided cognitive functions observed under STN-DBS.

Electrical Propagation of Vasodilatory Signals in Capillary Networks

We have developed a computational model to study electrical propagation of vasodilatory signals and arteriolar regulation of blood flow depending on the oxygen tension and agonist distribution in the capillary network. The involving key parameters of endothelial cell-to-cell electrical conductivity and plasma membrane area per unit volume were calibrated with the experimental data on an isolated endothelial tube of mouse skeletal feeding arteries. We have estimated the oxygen saturation parameters in terms of erythrocyte ATP release from the data of a left anterior descending coronary blood perfusion of dog. Regarding the acetylcholine-induced upstream conduction, our model shows that spatially uniform superfusion of acetylcholine attenuates the electrical signal propagation, and blocking calcium-activated potassium channels suppresses that attenuation. On the other hand, a local infusion of acetylcholine induces enhanced electrical propagation that corresponds to physiological relevance. Integrating the electrophysiology of endothelial tube and the electrophysiology/mechanics of a lumped arteriole, we show mechanistically that endothelial purinergic oxygen sensing of ATP released from erythrocytes and local infusion of acetylcholine are individually effective to induce vasodilatory signals to regulate blood flow in arterioles. We have recapitulated the upstream vasomotion in arterioles from the elevated oxygen tension in the downstream capillary domain. This study is a foundation for characterizing effective pharmaceutical strategies for ascending vasodilation and oxygenation.

Predicting the fMRI Signal Fluctuation with Recurrent Neural Networks Trained on Vascular Network Dynamics

Resting-state functional MRI (rs-fMRI) studies have revealed specific low-frequency hemodynamic signal fluctuations (<0.1 Hz) in the brain, which could be related to neuronal oscillations through the neurovascular coupling mechanism. Given the vascular origin of the fMRI signal, it remains challenging to separate the neural correlates of global rs-fMRI signal fluctuations from other confounding sources. However, the slow-oscillation detected from individual vessels by single-vessel fMRI presents strong correlation to neural oscillations. Here, we use recurrent neural networks (RNNs) to predict the future temporal evolution of the rs-fMRI slow oscillation from both rodent and human brains. The RNNs trained with vessel-specific rs-fMRI signals encode the unique brain oscillatory dynamic feature, presenting more effective prediction than the conventional autoregressive model. This RNN-based predictive modeling of rs-fMRI datasets from the Human Connectome Project (HCP) reveals brain state-specific characteristics, demonstrating an inverse relationship between the global rs-fMRI signal fluctuation with the internal default-mode network (DMN) correlation. The RNN prediction method presents a unique data-driven encoding scheme to specify potential brain state differences based on the global fMRI signal fluctuation, but not solely dependent on the global variance.

Ultra-slow Oscillations in fMRI and Resting-State Connectivity: Neuronal and Vascular Contributions and Technical Confounds

Ultra-slow, ∼0.1-Hz variations in the oxygenation level of brain blood are widely used as an fMRI-based surrogate of "resting-state" neuronal activity. The temporal correlations among these fluctuations across the brain are interpreted as "functional connections" for maps and neurological diagnostics. Ultra-slow variations in oxygenation follow a cascade. First, they closely track changes in arteriole diameter. Second, interpretable functional connections arise when the ultra-slow changes in amplitude of γ-band neuronal oscillations, which are shared across even far-flung but synaptically connected brain regions, entrain the ∼0.1-Hz vasomotor oscillation in diameter of local arterioles. Significant confounds to estimates of functional connectivity arise from residual vasomotor activity as well as arteriole dynamics driven by self-generated movements and subcortical common modulatory inputs. Last, methodological limitations of fMRI can lead to spurious functional connections. The neuronal generator of ultra-slow variations in γ-band amplitude, including that associated with self-generated movements, remains an open issue.

Association Between Brain Activation and Functional Connectivity

The origin of the "resting-state" brain activity recorded with functional magnetic resonance imaging (fMRI) is still uncertain. Here we provide evidence for the neurovascular origins of the amplitude of the low-frequency fluctuations (ALFF) and the local functional connectivity density (lFCD) by comparing them with task-induced blood-oxygen level dependent (BOLD) responses, which are considered a proxy for neuronal activation. Using fMRI data for 2 different tasks (Relational and Social) collected by the Human Connectome Project in 426 healthy adults, we show that ALFF and lFCD have linear associations with the BOLD response. This association was significantly attenuated by a novel task signal regression (TSR) procedure, indicating that task performance enhances lFCD and ALFF in activated regions. We also show that lFCD predicts BOLD activation patterns, as was recently shown for other functional connectivity metrics, which corroborates that resting functional connectivity architecture impacts brain activation responses. Thus, our findings indicate a common source for BOLD responses, ALFF and lFCD, which is consistent with the neurovascular origin of local hemodynamic synchrony presumably reflecting coordinated fluctuations in neuronal activity. This study also supports the development of task-evoked functional connectivity density mapping.

Shift-Invariant Canonical Polyadic Decomposition of Complex-Valued Multi-Subject fMRI Data With a Phase Sparsity Constraint

Canonical polyadic decomposition (CPD) of multi-subject complex-valued fMRI data can be used to provide spatially and temporally shared components among groups with both magnitude and phase information. However, the CPD model is not well formulated due to the large subject variability in the spatial and temporal modalities, as well as the high noise level in complex-valued fMRI data. Considering that the shift-invariant CPD can model temporal variability across subjects, we propose to further impose a phase sparsity constraint on the shared spatial maps to denoise the complex-valued components and to model the inter-subject spatial variability as well. More precisely, subject-specific time delays are first estimated for the complex-valued shared time courses in the framework of real-valued shift-invariant CPD. Source phase sparsity is then imposed on the complex-valued shared spatial maps. A smoothed l₀ norm is specifically used to reduce voxels with large phase values after phase de-ambiguity based on the small phase characteristic of BOLD-related voxels. The results from both the simulated and experimental fMRI data demonstrate improvements of the proposed method over three complex-valued algorithms, namely, tensor-based spatial ICA, shift-invariant CPD and CPD without spatiotemporal constraints. When comparing with a real-valued algorithm combining shift-invariant CPD and ICA, the proposed method detects 178.7% more contiguous task-related activations.

Interactions between stimuli-evoked cortical activity and spontaneous low frequency oscillations measured with neuronal calcium

Spontaneous brain activity has been widely used to map brain connectivity. The interactions between task-evoked brain responses and the spontaneous cortical oscillations, especially within the low frequency range of ~0.1 ​Hz, are not fully understood. Trial-to-trial variabilities in brain's response to sensory stimuli and the ability for brain to detect under noisy conditions suggest an appreciable impact of the brain state. Using a multimodality imaging platform, we simultaneously imaged neuronal Ca²⁺ and cerebral hemodynamics at baseline and in response to single-pulse forepaw stimuli in rat's somatosensory cortex. The high sensitivity of this system enables detection of responses to very weak and strong stimuli and real time determination of low frequency oscillations without averaging. Results show that the ongoing neuronal oscillations inversely modulate Ca²⁺ transients evoked by sensory stimuli. High intensity stimuli reset the spontaneous neuronal oscillations to an unpreferable excitability following the stimulus. Cerebral hemodynamic responses also inversely interact with the spontaneous hemodynamic oscillations, correlating with the neuronal Ca²⁺ transient changes. The results reveal competing interactions between spontaneous oscillations and stimulation-evoked brain activities in somatosensory cortex and the resultant hemodynamics.

Synchronized Astrocytic Ca2+ Responses in Neurovascular Coupling during Somatosensory Stimulation and for the Resting State

The role of astrocytes in neurovascular coupling (NVC) is unclear. Here, we applied a multimodality imaging approach to concomitantly measure synchronized neuronal or astrocytic Ca²⁺ and hemodynamic changes in the mouse somatosensory cortex at rest and during sensory electrical stimulation. Strikingly, we found that low-frequency stimulation (0.3-1 Hz), which consistently evokes fast neuronal Ca²⁺ transients (6.0 ± 2.7 ms latency) that always precede vascular responses, does not always elicit astrocytic Ca²⁺ transients (313 ± 65 ms latency). However, the magnitude of the hemodynamic response is increased when astrocytic transients occur, suggesting a facilitatory role of astrocytes in NVC. High-frequency stimulation (5-10 Hz) consistently evokes a large, delayed astrocytic Ca²⁺ accumulation (3.48 ± 0.09 s latency) that is temporarily associated with vasoconstriction, suggesting a role for astrocytes in resetting NVC. At rest, neuronal, but not astrocytic, Ca²⁺ fluctuations correlate with hemodynamic low-frequency oscillations. Taken together, these results support a role for astrocytes in modulating, but not triggering, NVC.

Cocaine Decreases Spontaneous Neuronal Activity and Increases Low-Frequency Neuronal and Hemodynamic Cortical Oscillations

Low-frequency oscillations (LFOs) in hemodynamics assessed by fMRI reflect synchronized neuronal activities and are the basis for mapping brain function and its disruption by drugs and disease. Here we assess if cocaine disrupts coupling between neuronal and vascular LFOs by simultaneously measuring cortical field potentials (FP) and cerebral blood flow (CBF) regarding their LFOs (0-1 Hz) spectral bandwidths in the somatosensory cortex of naïve and chronic cocaine-exposed rats at baseline and during cocaine intoxication. While across all conditions the dominant oscillation frequencies for FP and CBF LFOs were ~0.1 Hz, the bandwidth of FP LFOs was about 4.8 ± 0.67 times broader than that of CBF LFOs. Acute cocaine depressed high-frequency FP events but increased the relative intensity of neuronal and hemodynamic LFOs, an effect that was markedly accentuated in magnitude and duration in chronic cocaine-exposed animals. Neuronal LFOs were correlated with CBF LFOs in control animals but not in chronically cocaine-exposed animals, which suggests neurovascular uncoupling. The marked increases in neuronal LFOs with chronic cocaine, which we interpret to reflect increases in neuronal synchronization in the LFOs, and the uncoupling of hemodynamics with resting neuronal activities could contribute to brain dysfunction in cocaine abusers and confound the interpretation of fMRI studies.

Calcium channel blockade with nimodipine reverses MRI evidence of cerebral oedema following acute hypoxia

Acute cerebral hypoxia causes rapid calcium shifts leading to neuronal damage and death. Calcium channel antagonists improve outcomes in some clinical conditions, but mechanisms remain unclear. In 18 healthy participants we: (i) quantified with multiparametric MRI the effect of hypoxia on the thalamus, a region particularly sensitive to hypoxia, and on the whole brain in general; (ii) investigated how calcium channel antagonism with the drug nimodipine affects the brain response to hypoxia. Hypoxia resulted in a significant decrease in apparent diffusion coefficient (ADC), a measure particularly sensitive to cell swelling, in a widespread network of regions across the brain, and the thalamus in particular. In hypoxia, nimodipine significantly increased ADC in the same brain regions, normalizing ADC towards normoxia baseline. There was positive correlation between blood nimodipine levels and ADC change. In the thalamus, there was a significant decrease in the amplitude of low frequency fluctuations (ALFF) in resting state functional MRI and an apparent increase of grey matter volume in hypoxia, with the ALFF partially normalized towards normoxia baseline with nimodipine. This study provides further evidence that the brain response to acute hypoxia is mediated by calcium, and importantly that manipulation of intracellular calcium flux following hypoxia may reduce cerebral cytotoxic oedema.

An Autonomic Network: Synchrony Between Slow Rhythms of Pulse and Brain Resting State Is Associated with Personality and Emotions

The sympathetic system's role in modulating vasculature and its influence on emotions and personality led us to test the hypothesis that interactions between brain resting-state networks (RSNs) and pulse amplitude (indexing sympathetic activity) would be associated with emotions and personality. In 203 participants, we characterized RSN spatiotemporal characteristics, and phase-amplitude associations of RSN fluctuations with pulse and respiratory recordings. We found that RSNs are spatially reproducible within participants and were temporally associated with low frequencies (LFs < 0.1 Hz) in physiological signals. LF fluctuations in pulse amplitude were not related to cardiac electrical activity and preceded LF fluctuations in RSNs, while LF respiratory amplitude fluctuations followed LF fluctuations in RSNs. LF phase dispersion (PD) (lack of synchrony) between RSNs and pulse (PDpulse) (not respiratory) correlated with the common variability in measures of personality and emotions, with more synchrony being associated with more positive temperamental characteristics. Voxel-level PDpulse mapping revealed an "autonomic brain network," including sensory cortices and dorsal attention stream, with significant interactions with peripheral signals. Here, we uncover associations between pulse signal amplitude (presumably of sympathetic origin) and brain resting state, suggesting that interactions between central and autonomic nervous systems are important for characterizing personality and emotions.

Effect of combined naltrexone and bupropion therapy on the brain's functional connectivity

Background

The control of food intake in environments with easy access to highly rewarding foods is challenging to most modern societies. The combination of sustained release (SR) naltrexone and SR bupropion (NB32) has been used in weight-loss and obesity management. However, the effects of NB32 on the brain circuits implicated in the regulation of food intake are unknown. Here we used functional connectivity density (FCD) mapping to evaluate the effects of NB32 on resting brain FC.

Methods

Thirty-six healthy women underwent magnetic resonance imaging (MRI) before and after 4-week treatment with NB32 (n = 16) or with placebo (n = 20). In each imaging visit, a 5-min resting-state functional MRI scan was conducted after 15 h of fasting. The FC of brain regions showing significant group effects on FCD were subsequently assessed using seed-voxel correlation analyses. We characterized the associations between FCD measures and craving control scores in the Control of Eating Questionnaire.

Results

After NB32 treatment, the group showed lower local and global FCD than the placebo group in the right superior parietal cortex and lower local FCD in the left middle frontal gyrus. Seed-voxel correlation analysis for the right superior parietal cortex seed demonstrated higher positive FC with the dorsal anterior cingulate gyrus (ACC), bilateral insula, and left superior parietal gyrus and stronger negative FC with right inferior frontal gyrus and right superior parietal cortices for the NB32 than the placebo group. Further, the NB32 group showed a significant correlation between local FCD change after treatment in left middle frontal gyrus and craving control scores (r = 0.519, p = 0.039).

Conclusions

NB32 treatment decreased local and global FCD in superior parietal cortex and increased its connectivity with ACC (involved with saliency attribution), insula (interoception), and decreased local FCD in the medial prefrontal cortex (craving), which might underlie NB32 improved control over eating behaviors. ClinicalTrails.gov: NCT00711.

Ultra-Slow Single-Vessel BOLD and CBV-Based fMRI Spatiotemporal Dynamics and Their Correlation with Neuronal Intracellular Calcium Signals

Functional MRI has been used to map brain activity and functional connectivity based on the strength and temporal coherence of neurovascular-coupled hemodynamic signals. Here, single-vessel fMRI reveals vessel-specific correlation patterns in both rodents and humans. In anesthetized rats, fluctuations in the vessel-specific fMRI signal are correlated with the intracellular calcium signal measured in neighboring neurons. Further, the blood-oxygen-level-dependent (BOLD) signal from individual venules and the cerebral-blood-volume signal from individual arterioles show correlations at ultra-slow (<0.1 Hz), anesthetic-modulated rhythms. These data support a model that links neuronal activity to intrinsic oscillations in the cerebral vasculature, with a spatial correlation length of ∼2 mm for arterioles. In complementary data from awake human subjects, the BOLD signal is spatially correlated among sulcus veins and specified intracortical veins of the visual cortex at similar ultra-slow rhythms. These data support the use of fMRI to resolve functional connectivity at the level of single vessels.

Macroscale variation in resting-state neuronal activity and connectivity assessed by simultaneous calcium imaging, hemodynamic imaging and electrophysiology

Functional imaging of spontaneous activity continues to play an important role in the field of connectomics. The most common imaging signal used for these experiments is the blood-oxygen-level-dependent (BOLD) functional MRI (fMRI) signal, but how this signal relates to spontaneous neuronal activity remains incompletely understood. Genetically encoded calcium indicators represent a promising tool to study this problem, as they can provide brain-wide measurements of neuronal activity compared to point measurements afforded by electrophysiological recordings. However, the relationship between the calcium signal and neurophysiological parameters at the mesoscopic scale requires further systematic characterization. Therefore, we collected simultaneous resting-state measurements of electrophysiology, along with calcium and hemodynamic imaging, in lightly anesthetized mice to investigate two aims. First, we examined the relationship between each imaging signal and the simultaneously recorded electrophysiological signal in a single brain region, finding that both signals are better correlated with multi-unit activity compared to local field potentials, with the calcium signal possessing greater signal-to-noise ratio and regional specificity. Second, we used the resting-state imaging data to model the relationship between the calcium and hemodynamic signals across the brain. We found that this relationship varied across brain regions in a way that is consistent across animals, with delays increasing by600 ms towards posterior cortical regions. Furthermore, while overall functional connectivity (FC) measured by the hemodynamic signal is significantly correlated with FC measured by calcium, the two estimates were found to be significantly different. We hypothesize that these differences arise at least in part from the observed regional variation in the hemodynamic response. In total, this work highlights some of the caveats needed in interpreting hemodynamic-based measurements of FC, as well as the need for improved modeling methods to reduce this potential source of bias.

Entrainment of Arteriole Vasomotor Fluctuations by Neural Activity Is a Basis of Blood-Oxygenation-Level-Dependent "Resting-State" Connectivity

Resting-state signals in blood-oxygenation-level-dependent (BOLD) imaging are used to parcellate brain regions and define "functional connections" between regions. Yet a physiological link between fluctuations in blood oxygenation with those in neuronal signaling pathways is missing. We present evidence from studies on mouse cortex that modulation of vasomotion, i.e., intrinsic ultra-slow (0.1 Hz) fluctuations in arteriole diameter, provides this link. First, ultra-slow fluctuations in neuronal signaling, which occur as an envelope over γ-band activity, entrains vasomotion. Second, optogenetic manipulations confirm that entrainment is unidirectional. Third, co-fluctuations in the diameter of pairs of arterioles within the same hemisphere diminish to chance for separations >1.4 mm. Yet the diameters of arterioles in distant (>5 mm), mirrored transhemispheric sites strongly co-fluctuate; these correlations are diminished in acallosal mice. Fourth, fluctuations in arteriole diameter coherently drive fluctuations in blood oxygenation. Thus, entrainment of vasomotion links neuronal pathways to functional connections.

Functional connectivity structure of cortical calcium dynamics in anesthetized and awake mice

The interplay between hemodynamic-based markers of cortical activity (e.g. fMRI and optical intrinsic signal imaging), which are an indirect and relatively slow report of neural activity, and underlying synaptic electrical and metabolic activity through neurovascular coupling is a topic of ongoing research and debate. As application of resting state functional connectivity measures is extended further into topics such as brain development, aging and disease, the importance of understanding the fundamental physiological basis for functional connectivity will grow. Here we extend functional connectivity analysis from hemodynamic- to calcium-based imaging. Transgenic mice (n = 7) expressing a fluorescent calcium indicator (GCaMP6) driven by the Thy1 promoter in glutamatergic neurons were imaged transcranially in both anesthetized (using ketamine/xylazine) and awake states. Sequential LED illumination (λ = 454, 523, 595, 640nm) enabled concurrent imaging of both GCaMP6 fluorescence emission (corrected for hemoglobin absorption) and hemodynamics. Functional connectivity network maps were constructed for infraslow (0.009–0.08Hz), intermediate (0.08–0.4Hz), and high (0.4–4.0Hz) frequency bands. At infraslow and intermediate frequencies, commonly used in BOLD fMRI and fcOIS studies of functional connectivity and implicated in neurovascular coupling mechanisms, GCaMP6 and HbO₂ functional connectivity structures were in high agreement, both qualitatively and also quantitatively through a measure of spatial similarity. The spontaneous dynamics of both contrasts had the highest correlation when the GCaMP6 signal was delayed with a ~0.6–1.5s temporal offset. Within the higher-frequency delta band, sensitive to slow wave sleep oscillations in non-REM sleep and anesthesia, we evaluate the speed with which the connectivity analysis stabilized and found that the functional connectivity maps captured putative network structure within time window lengths as short as 30 seconds. Homotopic GCaMP6 functional connectivity maps at 0.4–4.0Hz in the anesthetized states show a striking correlated and anti-correlated structure along the anterior to posterior axis. This structure is potentially explained in part by observed propagation of delta-band activity from frontal somatomotor regions to visuoparietal areas. During awake imaging, this spatio-temporal quality is altered, and a more complex and detailed functional connectivity structure is observed. The combined calcium/hemoglobin imaging technique described here will enable the dissociation of changes in ionic and hemodynamic functional structure and neurovascular coupling and provide a framework for subsequent studies of neurological disease such as stroke.

Long-term optical imaging of neurovascular coupling in mouse cortex using GCaMP6f and intrinsic hemodynamic signals

Cerebral hemodynamics are modulated in response to changes in neuronal activity, a process termed neurovascular coupling (NVC), which can be disrupted by neuropsychiatric diseases (e.g., stroke, Alzheimer's disease). Thus, there is growing interest to image long-term NVC dynamics with high spatiotemporal resolutions. Here, by combining the use of a genetically-encoded calcium indicator with optical techniques, we develop a longitudinal multimodal optical imaging platform (MIP) that enabled time-lapse tracking of NVC over a relatively large field of view in the mouse somatosensory cortex at single cell and single vessel resolutions. Specifically, GCaMP6f was used as marker of neuronal activity, which along with MIP allowed us to simultaneously measure the changes in neuronal [Ca²⁺]_i fluorescence, cerebral blood flow velocity (CBFv) and hemodynamics longitudinally for more than eight weeks. We show that [Ca²⁺]_i fluorescence was detectable one week post viral injection and the damage to local microvasculature and perfusion recovered two weeks after injection. By three weeks post viral injection, maximal neuronal and CBFv responses to hindpaw stimulations were observed. Moreover, single neuronal activation in response to hindpaw stimulation was consistently recorded, followed by ∼2 s delayed dilation of contiguous microvessels. Additionally, resting-state spontaneous neuronal and hemodynamic oscillations were detectable throughout the eight weeks of study. Our results demonstrate the capability of MIP for longitudinal investigation of the organization and plasticity of the neurovascular network during resting state and during stimulation-evoked neuronal activation at high spatiotemporal resolutions.

Neural and metabolic basis of dynamic resting state fMRI

Resting state fMRI (rsfMRI) as a technique showed much initial promise for use in psychiatric and neurological diseases where diagnosis and treatment were difficult. To realize this promise, many groups have moved towards examining "dynamic rsfMRI," which relies on the assumption that rsfMRI measurements on short time scales remain relevant to the underlying neural and metabolic activity. Many dynamic rsfMRI studies have demonstrated differences between clinical or behavioral groups beyond what static rsfMRI measured, suggesting a neurometabolic basis. Correlative studies combining dynamic rsfMRI and other physiological measurements have supported this. However, they also indicate multiple mechanisms and, if using correlation alone, it is difficult to separate cause and effect. Hypothesis-driven studies are needed, a few of which have begun to illuminate the underlying neurometabolic mechanisms that shape observed differences in dynamic rsfMRI. While the number of potential noise sources, potential actual neurometabolic sources, and methodological considerations can seem overwhelming, dynamic rsfMRI provides a rich opportunity in systems neuroscience. Even an incrementally better understanding of the neurometabolic basis of dynamic rsfMRI would expand rsfMRI's research and clinical utility, and the studies described herein take the first steps on that path forward.

Alcohol affects brain functional connectivity and its coupling with behavior: greater effects in male heavy drinkers

Acute and chronic alcohol exposure significantly affect behavior but the underlying neurobiological mechanisms are still poorly understood. Here, we used functional connectivity density (FCD) mapping to study alcohol-related changes in resting brain activity and their association with behavior. Heavy drinkers (HD, N=16, 16 males) and normal controls (NM, N=24, 14 males) were tested after placebo and after acute alcohol administration. Group comparisons showed that NM had higher FCD in visual and prefrontal cortices, default mode network regions and thalamus, while HD had higher FCD in cerebellum. Acute alcohol significantly increased FCD within the thalamus, impaired cognitive and motor functions, and affected self-reports of mood/drug effects in both groups. Partial least squares regression showed that alcohol-induced changes in mood/drug effects were associated with changes in thalamic FCD in both groups. Disruptions in motor function were associated with increases in cerebellar FCD in NM and thalamus FCD in HD. Alcohol-induced declines in cognitive performance were associated with connectivity increases in visual cortex and thalamus in NM, but in HD, increases in precuneus FCD were associated with improved cognitive performance. Acute alcohol reduced 'neurocognitive coupling', the association between behavioral performance and FCD (indexing brain activity), an effect that was accentuated in HD compared with NM. Findings suggest that reduced cortical connectivity in HD contribute to decline in cognitive abilities associated with heavy alcohol consumption, whereas increased cerebellar connectivity in HD may have compensatory effects on behavioral performance. The results reveal how drinking history alters the association between brain FCD and individual differences in behavioral performance.

Statistical parametric mapping of stimuli evoked changes in total blood flow velocity in the mouse cortex obtained with extended-focus optical coherence microscopy

Functional magnetic resonance (fMRI) imaging is the current gold-standard in neuroimaging. fMRI exploits local changes in blood oxygenation to map neuronal activity over the entire brain. However, its spatial resolution is currently limited to a few hundreds of microns. Here we use extended-focus optical coherence microscopy (xfOCM) to quantitatively measure changes in blood flow velocity during functional hyperaemia at high spatio-temporal resolution in the somatosensory cortex of mice. As optical coherence microscopy acquires hundreds of depth slices simultaneously, blood flow velocity measurements can be performed over several vessels in parallel. We present the proof-of-principle of an optimised statistical parametric mapping framework to analyse quantitative blood flow timetraces acquired with xfOCM using the general linear model. We demonstrate the feasibility of generating maps of cortical hemodynamic reactivity at the capillary level with optical coherence microscopy. To validate our method, we exploited 3 stimulation paradigms, covering different temporal dynamics and stimulated limbs, and demonstrated its repeatability over 2 trials, separated by a week.

Resting-state functional MRI as a tool for evaluating brain hemodynamic responsiveness to external stimuli in rats

PURPOSE

Anesthesia is a major confounding factor in functional MRI (fMRI) experiments attributed to its effects on brain function. Recent evidence suggests that parameters obtained with resting-state fMRI (rs-fMRI) are coupled with anesthetic depth. Therefore, we investigated whether parameters obtained with rs-fMRI, such as functional connectivity (FC), are also directly related to blood-oxygen-level-dependent (BOLD) responses.

METHODS

A simple rs-fMRI protocol was implemented in a pharmacological fMRI study to evaluate the coupling between hemodynamic responses and FC under five anesthetics (α-chloralose, isoflurane, medetomidine, thiobutabarbital, and urethane). Temporal change in the FC was evaluated at 1-hour interval. Supplementary forepaw stimulation experiments were also conducted.

RESULTS

Under thiobutabarbital anesthesia, FC was clearly coupled with nicotine-induced BOLD responses. Good correlation values were also obtained under isoflurane and medetomidine anesthesia. The observations in the thiobutabarbital group were supported by forepaw stimulation experiments. Additionally, the rs-fMRI protocol revealed significant temporal changes in the FC in the α-chloralose, thiobutabarbital, and urethane groups.

CONCLUSION

Our results suggest that FC can be used to estimate brain hemodynamic responsiveness to stimuli and evaluate the level and temporal changes of anesthesia. Therefore, analysis of the fMRI baseline signal may be highly valuable tool for controlling the outcome of preclinical fMRI experiments. Magn Reson Med 78:1136-1146, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Cocaine attenuates blood flow but not neuronal responses to stimulation while preserving neurovascular coupling for resting brain activity

Cocaine affects neuronal activity and constricts cerebral blood vessels, making it difficult to determine whether cocaine-induced changes in cerebral blood flow (CBF) reflect neuronal activation or its vasoactive effects. Here we assessed the effects of acute cocaine on both resting-state and stimulation responses to investigate cocaine's effects on neurovascular coupling and to differentiate its effects on neuronal activity from its vasoactive actions. We concurrently measured cortical field potentials via thinned-skull electroencephalography recordings and CBF with laser Doppler flowmetry in the rat's somatosensory cortex for both resting state and forepaw stimulation before and following cocaine administration (1 mg kg(-1), intravenously). Results show both resting-state field potentials and CBF were depressed after cocaine administration (19.8±4.7% and 52.1±13.4%, respectively) and these changes were strongly correlated with each other (r=0.81, P<0.001), indicating that cocaine did not affect neurovascular coupling at rest and that the reduction in resting CBF reflected reduction in synchronized spontaneous neuronal activity rather than vasoconstriction. In contrast, the forepaw stimulation-evoked neuronal activity was not changed by cocaine (P=0.244), whereas the CBF to the stimulation was reduced 49.9±2.6% (P=0.028) gradually recovering ∼20 min after cocaine injection, indicating that neurovascular coupling during stimulation was temporarily disrupted by cocaine. Neurovascular uncoupling by cocaine during stimulation but not during rest indicates that distinct processes might underlie neurovascular regulation for both stimulation and spontaneous activity. The greater reductions by cocaine to the stimulation-induced CBF increases than to the background CBF should be considered when interpreting functional MRI studies comparing activation responses between controls and cocaine abusers. Neurovascular uncoupling could contribute to cocaine's neurotoxicity, particularly for stimulation conditions when CBF might be insufficient to cover for the energetic demands of neuronal tissue.

Comparison of peripheral near-infrared spectroscopy low-frequency oscillations to other denoising methods in resting state functional MRI with ultrahigh temporal resolution

PURPOSE

Functional MRI (fMRI) blood-oxygen level-dependent (BOLD) signals result not only from neuronal activation, but also from nonneuronal physiological processes. These changes, especially in the low-frequency domain (0.01-0.2 Hz), can significantly confound inferences about neuronal processes. It is crucial to effectively identify these nuisance low-frequency oscillations (LFOs).

METHOD

A high temporal resolution (repetition time, ∼0.5 s) fMRI resting state study was conducted with simultaneous physiological measurements to compare LFOs measured directly by near-infrared spectroscopy (NIRS) in the periphery and three methods that model LFOs from the respiration or cardiac signal: 1) the respiration volume per time (RVT), 2) the respiratory variation (RVRRF), and 3) the cardiac variation method (HRCRF). The LFO noise regressors from these methods were compared temporally and spatially as well as in their denoising efficiency.

RESULTS

Methods were not highly correlated with one another, temporally or spatially. The set of two NIRS LFOs combined explained over 13% of BOLD signal variance and explained equal or more variance than HRCRF and RVRRF or RVT combined (in 14 of 16 participants).

CONCLUSION

LFOs collected using NIRS in the periphery contain distinct temporal and spatial information about the LFOs in BOLD fMRI that is not contained in current low-frequency denoising methods derived from respiration and cardiac pulsation. Magn Reson Med 76:1697-1707, 2016. © 2016 International Society for Magnetic Resonance in Medicine.

Fully distributed absolute blood flow velocity measurement for middle cerebral arteries using Doppler optical coherence tomography

Doppler optical coherence tomography (DOCT) is considered one of the most promising functional imaging modalities for neuro biology research and has demonstrated the ability to quantify cerebral blood flow velocity at a high accuracy. However, the measurement of total absolute blood flow velocity (BFV) of major cerebral arteries is still a difficult problem since it is related to vessel geometry. In this paper, we present a volumetric vessel reconstruction approach that is capable of measuring the absolute BFV distributed along the entire middle cerebral artery (MCA) within a large field-of-view. The Doppler angle at each point of the MCA, representing the vessel geometry, is derived analytically by localizing the artery from pure DOCT images through vessel segmentation and skeletonization. Our approach could achieve automatic quantification of the fully distributed absolute BFV across different vessel branches. Experiments on rodents using swept-source optical coherence tomography showed that our approach was able to reveal the consequences of permanent MCA occlusion with absolute BFV measurement.

Neuroimaging of Parkinson's disease: Expanding views

Advances in molecular and structural and functional neuroimaging are rapidly expanding the complexity of neurobiological understanding of Parkinson's disease (PD). This review article begins with an introduction to PD neurobiology as a foundation for interpreting neuroimaging findings that may further lead to more integrated and comprehensive understanding of PD. Diverse areas of PD neuroimaging are then reviewed and summarized, including positron emission tomography, single photon emission computed tomography, magnetic resonance spectroscopy and imaging, transcranial sonography, magnetoencephalography, and multimodal imaging, with focus on human studies published over the last five years. These included studies on differential diagnosis, co-morbidity, genetic and prodromal PD, and treatments from L-DOPA to brain stimulation approaches, transplantation and gene therapies. Overall, neuroimaging has shown that PD is a neurodegenerative disorder involving many neurotransmitters, brain regions, structural and functional connections, and neurocognitive systems. A broad neurobiological understanding of PD will be essential for translational efforts to develop better treatments and preventive strategies. Many questions remain and we conclude with some suggestions for future directions of neuroimaging of PD.

Patterns of Spontaneous Local Network Activity in Developing Cerebral Cortex: Relationship to Adult Cognitive Function

Detecting neurodevelopμental disorders of cognition at the earliest possible stages could assist in understanding them mechanistically and ultimately in treating them. Finding early physiological predictors that could be visualized with functional neuroimaging would represent an important advance in this regard. We hypothesized that one potential source of physiological predictors is the spontaneous local network activity prominent during specific periods in development. To test this we used calcium imaging in brain slices and analyzed variations in the frequency and intensity of this early activity in one area, the entorhinal cortex (EC), in order to correlate early activity with level of cognitive function later in life. We focused on EC because of its known role in different types of cognitive processes and because it is an area where spontaneous activity is prominent during early postnatal development in rodent models of cortical development. Using rat strains (Long-Evans, Wistar, Sprague-Dawley and Brattleboro) known to differ in cognitive performance in adulthood we asked whether neonatal animals exhibit corresponding strain-related differences in EC spontaneous activity. Our results show significant differences in this activity between strains: compared to a high cognitive-performing strain, we consistently found an increase in frequency and decrease in intensity in neonates from three lower performing strains. Activity was most different in one strain considered a model of schizophrenia-like psychopathology. While we cannot necessarily infer a causal relationship between early activity and adult cognition our findings suggest that the pattern of spontaneous activity in development could be an early predictor of a developmental trajectory advancing toward sub-optimal cognitive performance in adulthood. Our results further suggest that the strength of dopaminergic signaling, by setting the balance between excitation and inhibition, is a potential underlying mechanism that could explain the observed differences in early spontaneous activity patterns.

Quantitative microvascular hemoglobin mapping using visible light spectroscopic Optical Coherence Tomography

Quantification of chromophore concentrations in reflectance mode remains a major challenge for biomedical optics. Spectroscopic Optical Coherence Tomography (SOCT) provides depth-resolved spectroscopic information necessary for quantitative analysis of chromophores, like hemoglobin, but conventional SOCT analysis methods are applicable only to well-defined specular reflections, which may be absent in highly scattering biological tissue. Here, by fitting of the dynamic scattering signal spectrum in the OCT angiogram using a forward model of light propagation, we quantitatively determine hemoglobin concentrations directly. Importantly, this methodology enables mapping of both oxygen saturation and total hemoglobin concentration, or alternatively, oxyhemoglobin and deoxyhemoglobin concentration, simultaneously. Quantification was verified by ex vivo blood measurements at various pO2 and hematocrit levels. Imaging results from the rodent brain and retina are presented. Confounds including noise and scattering, as well as potential clinical applications, are discussed.