Neuronal and physiological correlation to hemodynamic resting-state fluctuations in health and disease

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.

Relationships between correlated spikes, oxygen and LFP in the resting-state primate

Resting-state functional MRI (rsfMRI) provides a view of human brain organization based on correlation patterns of blood oxygen level dependent (BOLD) signals recorded across the whole brain. The neural basis of resting-state BOLD fluctuations and their correlation remains poorly understood. We simultaneously recorded oxygen level, spikes, and local field potential (LFP) at multiple sites in awake, resting monkeys. Following a spike, the average local oxygen and LFP voltage responses each resemble a task-driven BOLD response, with LFP preceding oxygen by 0.5 s. Between sites, features of the long-range correlation patterns of oxygen, LFP, and spikes are similar to features seen in rsfMRI. Most of the variance shared between sites lies in the infraslow frequency band (0.01-0.1 Hz) and in the infraslow envelope of higher-frequency bands (e.g. gamma LFP). While gamma LFP and infraslow LFP are both strong correlates of local oxygen, infraslow LFP explains significantly more of the variance shared between correlated oxygen signals than any other electrophysiological signal. Together these findings are consistent with a causal relationship between infraslow LFP and long-range oxygen correlations in the resting state.

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.

In vivo imaging of calcium and glutamate responses to intracortical microstimulation reveals distinct temporal responses of the neuropil and somatic compartments in layer II/III neurons

OBJECTIVE

Intracortical microelectrode implants can generate a tissue response hallmarked by glial scarring and neuron cell death within 100-150 μm of the biomaterial device. Many have proposed that any performance decline in intracortical microstimulation (ICMS) due to this foreign body tissue response could be offset by increasing the stimulation amplitude. The mechanisms of this approach are unclear, however, as there has not been consensus on how increasing amplitude affects the spatial and temporal recruitment patterns of ICMS.

APPROACH

We clarify these unknowns using in vivo two-photon imaging of mice transgenically expressing the calcium sensor GCaMP6s in Thy1 neurons or virally expressing the glutamate sensor iGluSnFr in neurons. Calcium and neurotransmitter activity are tracked in the neuronal somas and neuropil during long-train stimulation in Layer II/III of somatosensory cortex.

MAIN RESULTS

Neural calcium activity and glutamate release are dense and strongest within 20-40 μm around the electrode, falling off with distance from the electrode. Neuronal calcium increases with higher amplitude stimulations. During prolonged stimulation trains, a sub-population of somas fail to maintain calcium activity. Interestingly, neuropil calcium activity is 3-fold less correlated to somatic calcium activity for cells that drop-out during the long stimulation train compared to cells that sustain activity throughout the train. Glutamate release is apparent only within 20 μm of the electrode and is sustained for at least 10s after cessation of the 15 and 20 μA stimulation train, but not lower amplitudes.

SIGNIFICANCE

These results demonstrate that increasing amplitude can increase the radius and intensity of neural recruitment, but it also alters the temporal response of some neurons. Further, dense glutamate release is highest within the first 20 μm of the electrode site even at high amplitudes, suggesting that there may be spatial limitations to the amplitude parameter space. The glutamate elevation outlasts stimulation, suggesting that high-amplitude stimulation may affect neurotransmitter re-uptake. This ultimately suggests that increasing the amplitude of ICMS device stimulation may fundamentally alter the temporal neural response, which could have implications for using amplitude to improve the ICMS effect or "offset" the effects of glial scarring.

Multi-scale, multi-modal analysis uncovers complex relationship at the brain tissue-implant neural interface: new emphasis on the biological interface

OBJECTIVE

Implantable neural electrode devices are important tools for neuroscience research and have an increasing range of clinical applications. However, the intricacies of the biological response after implantation, and their ultimate impact on recording performance, remain challenging to elucidate. Establishing a relationship between the neurobiology and chronic recording performance is confounded by technical challenges related to traditional electrophysiological, material, and histological limitations. This can greatly impact the interpretations of results pertaining to device performance and tissue health surrounding the implant.

APPROACH

In this work, electrophysiological activity and immunohistological analysis are compared after controlling for motion artifacts, quiescent neuronal activity, and material failure of devices in order to better understand the relationship between histology and electrophysiological outcomes.

MAIN RESULTS

Even after carefully accounting for these factors, the presence of viable neurons and lack of glial scarring does not convey single unit recording performance.

SIGNIFICANCE

To better understand the biological factors influencing neural activity, detailed cellular and molecular tissue responses were examined. Decreases in neural activity and blood oxygenation in the tissue surrounding the implant, shift in expression levels of vesicular transporter proteins and ion channels, axon and myelin injury, and interrupted blood flow in nearby capillaries can impact neural activity around implanted neural interfaces. Combined, these tissue changes highlight the need for more comprehensive, basic science research to elucidate the relationship between biology and chronic electrophysiology performance in order to advance neural technologies.

Spontaneous activity forms a foundation for odor-evoked activation maps in the rat olfactory bulb

Fluctuations in spontaneous activity have been observed by many neuroimaging techniques, but because these resting-state changes are not evoked by stimuli, it is difficult to determine how they relate to task-evoked activations. We conducted multi-modal neuroimaging scans of the rat olfactory bulb, both with and without odor, to examine interaction between spontaneous and evoked activities. Independent component analysis of spontaneous fluctuations revealed resting-state networks, and odor-evoked changes revealed activation maps. We constructed simulated activation maps using resting-state networks that were highly correlated to evoked activation maps. Simulated activation maps derived by intrinsic optical signal (IOS), which covers the dorsal portion of the glomerular sheet, significantly differentiated one odor's evoked activation map from the other two. To test the hypothesis that spontaneous activity of the entire glomerular sheet is relevant for representing odor-evoked activations, we used functional magnetic resonance imaging (fMRI) to map the entire glomerular sheet. In contrast to the IOS results, the fMRI-derived simulated activation maps significantly differentiated all three odors' evoked activation maps. Importantly, no evoked activation maps could be significantly differentiated using simulated activation maps produced using phase-randomized resting-state networks. Given that some highly organized resting-state networks did not correlate with any odors' evoked activation maps, we posit that these resting-state networks may characterize evoked activation maps associated with odors not studied. These results emphasize that fluctuations in spontaneous activity form a foundation for active processing, signifying the relevance of resting-state mapping to functional neuroimaging.

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.

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.

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.

Phase-synchronization-based parcellation of resting state fMRI signals reveals topographically organized clusters in early visual cortex

Resting-state fMRI is widely used to study brain function and connectivity. However, interpreting patterns of resting state (RS) fMRI activity remains challenging as they may arise from different neuronal mechanisms than those triggered by exogenous events. Currently, this limits the use of RS-fMRI for understanding cortical function in health and disease. Here, we examine the phase synchronization (PS) properties of blood-oxygen level dependent (BOLD) signals obtained during visual field mapping (VFM) and RS with 7T fMRI. This data-driven approach exploits spatiotemporal covariations in the phase of BOLD recordings to establish the presence of clusters of synchronized activity. We find that, in both VFM and RS data, selecting the most synchronized neighboring recording sites identifies spatially localized PS clusters that follow the topographic organization of the visual cortex. However, in activity obtained during VFM, PS is spatially more extensive than in RS activity, likely reflecting stimulus-driven interactions between local responses. Nevertheless, the similarity of the PS clusters obtained for RS and stimulus-driven fMRI suggest that they share a common neuroanatomical origin. Our finding justifies and facilitates direct comparison of RS and stimulus-evoked activity.

Mesoscale Mapping of Mouse Cortex Reveals Frequency-Dependent Cycling between Distinct Macroscale Functional Modules

Connectivity mapping based on resting-state activity in mice has revealed functional motifs of correlated activity. However, the rules by which motifs organize into larger functional modules that lead to hemisphere wide spatial-temporal activity sequences is not clear. We explore cortical activity parcellation in head-fixed, quiet awake GCaMP6 mice from both sexes by using mesoscopic calcium imaging. Spectral decomposition of spontaneous cortical activity revealed the presence of two dominant frequency modes (<1 and ∼3 Hz), each of them associated with a unique spatial signature of cortical macro-parcellation not predicted by classical cytoarchitectonic definitions of cortical areas. Based on assessment of 0.1-1 Hz activity, we define two macro-organizing principles: the first being a rotating polymodal-association pinwheel structure around which activity flows sequentially from visual to barrel then to hindlimb somatosensory; the second principle is correlated activity symmetry planes that exist on many levels within a single domain such as intrahemispheric reflections of sensory and motor cortices. In contrast, higher frequency activity >1 Hz yielded two larger clusters of coactivated areas with an enlarged default mode network-like posterior region. We suggest that the apparent constrained structure for intra-areal cortical activity flow could be exploited in future efforts to normalize activity in diseases of the nervous system.SIGNIFICANCE STATEMENT Increasingly, functional connectivity mapping of spontaneous activity is being used to reveal the organization of the brain. However, because the brain operates across multiple space and time domains a more detailed understanding of this organization is necessary. We used in vivo wide-field calcium imaging of the indicator GCaMP6 in head-fixed, awake mice to characterize the organization of spontaneous cortical activity at different spatiotemporal scales. Correlation analysis defines the presence of two to three superclusters of activity that span traditionally defined functional territories and were frequency dependent. This work helps define the rules for how different cortical areas interact in time and space. We provide a framework necessary for future studies that explore functional reorganization of brain circuits in disease models.

Biomarkers, designs, and interpretations of resting-state fMRI in translational pharmacological research: A review of state-of-the-Art, challenges, and opportunities for studying brain chemistry

A decade of research and development in resting-state functional MRI (RSfMRI) has opened new translational and clinical research frontiers. This review aims to bridge between technical and clinical researchers who seek reliable neuroimaging biomarkers for studying drug interactions with the brain. About 85 pharma-RSfMRI studies using BOLD signal (75% of all) or arterial spin labeling (ASL) were surveyed to investigate the acute effects of psychoactive drugs. Experimental designs and objectives include drug fingerprinting dose-response evaluation, biomarker validation and calibration, and translational studies. Common biomarkers in these studies include functional connectivity, graph metrics, cerebral blood flow and the amplitude and spectrum of BOLD fluctuations. Overall, RSfMRI-derived biomarkers seem to be sensitive to spatiotemporal dynamics of drug interactions with the brain. However, drugs cause both central and peripheral effects, thus exacerbate difficulties related to biological confounds, structured noise from motion and physiological confounds, as well as modeling and inference testing. Currently, these issues are not well explored, and heterogeneities in experimental design, data acquisition and preprocessing make comparative or meta-analysis of existing reports impossible. A unifying collaborative framework for data-sharing and data-mining is thus necessary for investigating the commonalities and differences in biomarker sensitivity and specificity, and establishing guidelines. Multimodal datasets including sham-placebo or active control sessions and repeated measurements of various psychometric, physiological, metabolic and neuroimaging phenotypes are essential for pharmacokinetic/pharmacodynamic modeling and interpretation of the findings. We provide a list of basic minimum and advanced options that can be considered in design and analyses of future pharma-RSfMRI studies. Hum Brain Mapp 38:2276-2325, 2017. © 2017 Wiley Periodicals, Inc.

Correlation of hemodynamic and fluorescence signals under resting state conditions in mice's barrel field cortex

Both neurons and astrocytes are known to affect local vascular response in the brain following neuronal activity. In order to differentiate the contributions of each cell type to the hemodynamic response during stimulation and resting state, intrinsic optical signal (IOI) was recorded synchronized with fluorescence imaging of calcium concentration sensitive dye Oregon Green BAPTA-1 AM. By changing the stimulation parameters (frequency and duration), it was possible to individually promote neuronal and glial responses and to compare them to levels of oxy (HbO), deoxy (HbR) and total (HbT) hemoglobin concentrations. Finally, resting state recordings were done to investigate the possible correlation between hemoglobin fluctuation and calcium transients, based on different frequency bands associated either with neuronal or glial activity.

Photoelectric artefact from optogenetics and imaging on microelectrodes and bioelectronics: New Challenges and Opportunities

Bioelectronics, electronic technologies that interface with biological systems, are experiencing rapid growth in terms of technology development and applications, especially in neuroscience and neuroprosthetic research. The parallel growth with optogenetics and in vivo multi-photon microscopy has also begun to generate great enthusiasm for simultaneous applications with bioelectronic technologies. However, emerging research showing artefact contaminated data highlight the need for understanding the fundamental physical principles that critically impact experimental results and complicate their interpretation. This review covers four major topics: 1) material dependent properties of the photoelectric effect (conductor, semiconductor, organic, photoelectric work function (band gap)); 2) optic dependent properties of the photoelectric effect (single photon, multiphoton, entangled biphoton, intensity, wavelength, coherence); 3) strategies and limitations for avoiding/minimizing photoelectric effects; and 4) advantages of and applications for light-based bioelectronics (photo-bioelectronics).