Coupling between somatosensory evoked potentials and hemodynamic response in the rat

Effects of focal cortical cooling on somatosensory evoked potentials in rats

Although the focal brain cooling technique is widely used to examine brain function, the effects of cortical temperature at various levels on sensory information processing and neural mechanisms remain underexplored. To elucidate the mechanisms of temperature modulation in somatosensory processing, this study aimed to examine how P1 and N1 deflections of somatosensory evoked potentials (SEPs) depend on cortical temperature and how excitatory and inhibitory inputs contribute to this temperature dependency. SEPs were generated through electrical stimulation of the contralateral forepaw in anesthetized rats. The SEPs were recorded while cortical temperatures were altered between 17-38 °C either without any antagonists, with a gamma-aminobutyric acid type A (GABAA) receptor antagonist (gabazine), with aminomethylphosphonic acid (AMPA) receptor antagonist (NBQX), or with N-Methyl-D-aspartic acid (NMDA) receptor antagonist ([R]-CPP). The effects of different gabazine concentrations (0, 1, and 10 µM) were examined in the 35-38 °C range. The P1/N1 amplitudes and their peak-to-peak differences plotted against cortical temperature showed an inverted U relationship with a maximum at approximately 27.5 °C when no antagonists were administered. The negative correlation between these amplitudes and temperatures of ≥ 27.5 °C plateaued after gabazine administration, which occurred progressively as the gabazine concentration increased. In contrast, the correlation remained negative after the administration of NBQX and (R)-CPP. These results suggest that GABAergic inhibitory inputs contribute to the negative correlation between SEP amplitude and cortical temperature around the physiological cortical temperature.

Identifying Intraoperative Spinal Cord Injury Location from Somatosensory Evoked Potentials' Time-Frequency Components

Excessive distraction in corrective spine surgery can lead to iatrogenic distraction spinal cord injury. Diagnosis of the location of the spinal cord injury helps in early removal of the injury source. The time-frequency components of the somatosensory evoked potential have been reported to provide information on the location of spinal cord injury, but most studies have focused on contusion injuries of the cervical spine. In this study, we established 19 rat models of distraction spinal cord injury at different levels and collected the somatosensory evoked potentials of the hindlimb and extracted their time-frequency components. Subsequently, we used k-medoid clustering and naive Bayes to classify spinal cord injury at the C5 and C6 level, as well as spinal cord injury at the cervical, thoracic, and lumbar spine, respectively. The results showed that there was a significant delay in the latency of the time-frequency components distributed between 15 and 30 ms and 50 and 150 Hz in all spinal cord injury groups. The overall classification accuracy was 88.28% and 84.87%. The results demonstrate that the k-medoid clustering and naive Bayes methods are capable of extracting the time-frequency component information depending on the spinal cord injury location and suggest that the somatosensory evoked potential has the potential to diagnose the location of a spinal cord injury.

Different Time-frequency Distribution Patterns of Somatosensory Evoked Potentials in Dual- and Single-level Spinal Cord Compression

Among patients with cervical myelopathy, the most common level of stenosis at spinal cord of all ages was reported to be between cervical levels C5-6. Previous studies found that time-frequency components (TFCs) of somatosensory evoked potentials (SEPs) possess location information of spinal cord injury (SCI) in single-level deficits in the spinal cord. However, the clinical reality is that there are multiple compressions at multiple spinal cord segments. This study proposed a new algorithm to differentiate distribution patterns of SEP TFCs between the dual-level compression and the corresponding single-level compression, which is potentials in providing precise diagnosis of cervical myelopathy. In the present animal study, a group of rats with dual-level compressive (C5+6) injury to cervical spinal cord was investigated. SEPs were collected at 2 weeks after surgery, while SEP TFCs were calculated. The SEP TFCs under dual-level compression were compared to an existent dataset with one sham control group and three single level compression groups at C4, C5, C6. Behavioral evaluation showed very similar scale of injury severity between individual rats, while histology evaluation confirmed the precise location of injury. According to time-frequency distribution patterns, it showed that the middle-energy components of dual-level showed similar patterns as that of each single-level group. In addition, the low-energy components of the dual-level C5+6 group had the highest correlation with C5 (R=0.3423, p<0.01) and C6 (R=0.4000, p<0.01) groups, but much lower with C4 group (R=0.1071, p=0.012). These results indicated that SEP TFCs components possess information regarding the location of neurological lesion after spinal cord compression. It preliminarily demonstrated that SEP TFCs are likely a useful measure to provide location information of neurological lesions after compression SCI.

Multichannel near-infrared spectroscopy brain imaging system for small animals in mobile conditions

Significance: We propose a customized animal-specific head cap and an near-infrared spectroscopy (NIRS) system to obtain NIRS signals in mobile small animals. NIRS studies in mobile small animals provide a feasible solution for comprehensive brain function studies. Aim: We aim to develop and validate a multichannel NIRS system capable of performing functional brain imaging along with a closed-box stimulation kit for small animals in mobile conditions. Approach: The customized NIRS system uses light-weight long optical fibers, along with a customized light-weight head cap to securely attach the optical fibers to the mouse. A customized stimulation box was designed to perform various stimuli in a controlled environment. The system performance was tested in a visual stimulation task on eight anesthetized mice and eight freely moving mice. Results: Following the visual stimulation task, we observed a significant stimulation-related oxyhemoglobin (HbO) increase in the visual cortex of freely moving mice during the task. In contrast, HbO concentration did not change significantly in the visual cortex of anesthetized mice. Conclusions: We demonstrate the feasibility of a wearable, multichannel NIRS system for small animals in a less confined experimental design.

Electrical Source Imaging in Freely Moving Rats: Evaluation of a 12-Electrode Cortical Electroencephalography System

This work presents and evaluates a 12-electrode intracranial electroencephalography system developed at the National Institute of Mental Health (Klecany, Czech Republic) in terms of an electrical source imaging (ESI) technique in rats. The electrode system was originally designed for translational research purposes. This study demonstrates that it is also possible to use this well-established system for ESI, and estimates its precision, accuracy, and limitations. Furthermore, this paper sets a methodological basis for future implants. Source localization quality is evaluated using three approaches based on surrogate data, physical phantom measurements, and in vivo experiments. The forward model for source localization is obtained from the FieldTrip-SimBio pipeline using the finite-element method. Rat brain tissue extracted from a magnetic resonance imaging template is approximated by a single-compartment homogeneous tetrahedral head model. Four inverse solvers were tested: standardized low-resolution brain electromagnetic tomography, exact low-resolution brain electromagnetic tomography (eLORETA), linear constrained minimum variance (LCMV), and dynamic imaging of coherent sources. Based on surrogate data, this paper evaluates the accuracy and precision of all solvers within the brain volume using error distance and reliability maps. The mean error distance over the whole brain was found to be the lowest in the eLORETA solution through signal to noise ratios (SNRs) (0.2 mm for 25 dB SNR). The LCMV outperformed eLORETA under higher SNR conditions, and exhibiting higher spatial precision. Both of these inverse solvers provided accurate results in a phantom experiment (1.6 mm mean error distance across shallow and 2.6 mm across subcortical testing dipoles). Utilizing the developed technique in freely moving rats, an auditory steady-state response experiment provided results in line with previously reported findings. The obtained results support the idea of utilizing a 12-electrode system for ESI and using it as a solid basis for the development of future ESI dedicated implants.

Exploration of Cervical Myelopathy Location From Somatosensory Evoked Potentials Using Random Forests Classification

Studies using time-frequency analysis have reported that somatosensory evoked potentials provide information regarding the location of spinal cord injury. However, a better understanding of the time-frequency components derived from somatosensory evoked potentials is essential for developing more reliable algorithms that can diagnosis level (location) of cervical injury. In the present study, we proposed a random forests machine learning approach, for separating somatosensory evoked potentials depending on spinal cord state. For data acquisition, we established rat models of compression spinal cord injury at the C4, C5, and C6 levels to induce cervical myelopathy. After making the compression injury, we collected somatosensory evoked potentials and extracted their time-frequency components. We then used the random forests classification system to analyze the evoked potential dataset that was obtained from the three groups of model rats. Evaluation of the classifier performance revealed an overall classification accuracy of 84.72%, confirming that the random forests method was able to separate the time-frequency components of somatosensory evoked potentials from rats under different conditions. Features of the time-frequency components contained information that could identify the location of the cervical spinal cord injury, demonstrating the potential benefits of using time-frequency components of somatosensory evoked potentials to diagnose the level of cervical injury in cervical myelopathy.

Validating Non-invasive EEG Source Imaging Using Optimal Electrode Configurations on a Representative Rat Head Model

The curtain of technical limitations impeding rat multichannel non-invasive electroencephalography (EEG) has risen. Given the importance of this preclinical model, development and validation of EEG source imaging (ESI) is essential. We investigate the validity of well-known human ESI methodologies in rats which individual tissue geometries have been approximated by those extracted from an MRI template, leading also to imprecision in electrode localizations. With the half and fifth sensitivity volumes we determine both the theoretical minimum electrode separation for non-redundant scalp EEG measurements and the electrode sensitivity resolution, which vary over the scalp because of the head geometry. According to our results, electrodes should be at least ~3 to 3.5 mm apart for an optimal configuration. The sensitivity resolution is generally worse for electrodes at the boundaries of the scalp measured region, though, by analogy with human montages, concentrates the sensitivity enough to localize sources. Cramér-Rao lower bounds of source localization errors indicate it is theoretically possible to achieve ESI accuracy at the level of anatomical structures, such as the stimulus-specific somatosensory areas, using the template. More validation for this approximation is provided through the comparison between the template and the individual lead field matrices, for several rats. Finally, using well-accepted inverse methods, we demonstrate that somatosensory ESI is not only expected but also allows exploring unknown phenomena related to global sensory integration. Inheriting the advantages and pitfalls of human ESI, rat ESI will boost the understanding of brain pathophysiological mechanisms and the evaluation of ESI methodologies, new pharmacological treatments and ESI-based biomarkers.

Application of Functional Near-Infrared Spectroscopy to the Study of Brain Function in Humans and Animal Models

Functional near-infrared spectroscopy (fNIRS) is a noninvasive optical imaging technique that indirectly assesses neuronal activity by measuring changes in oxygenated and deoxygenated hemoglobin in tissues using near-infrared light. fNIRS has been used not only to investigate cortical activity in healthy human subjects and animals but also to reveal abnormalities in brain function in patients suffering from neurological and psychiatric disorders and in animals that exhibit disease conditions. Because of its safety, quietness, resistance to motion artifacts, and portability, fNIRS has become a tool to complement conventional imaging techniques in measuring hemodynamic responses while a subject performs diverse cognitive and behavioral tasks in test settings that are more ecologically relevant and involve social interaction. In this review, we introduce the basic principles of fNIRS and discuss the application of this technique in human and animal studies.

Component analysis of somatosensory evoked potentials for identifying spinal cord injury location

This study aims to determine whether the time-frequency components (TFCs) of somatosensory evoked potentials (SEPs) can be used to identify the specific location of a compressive spinal cord injury using a classification technique. Waveforms of SEPs after compressive injuries at various locations (C4, C5 and C6) in rat spinal cords were decomposed into a series of TFCs using a high-resolution time-frequency analysis method. A classification method based on support vector machine (SVM) was applied to the distributions of these TFCs among different pathological locations. The difference among injury locations manifests itself in different categories of SEP TFCs. High-energy TFCs of normal-state SEPs have significantly higher power and frequency than those of injury-state SEPs. The location of C5 is characterized by a unique distribution pattern of middle-energy TFCs. The difference between C4 and C6 is evidenced by the distribution pattern of low-energy TFCs. The proposed classification method based on SEP TFCs offers a discrimination accuracy of 80.2%. In this study, meaningful information contained in various SEP components was investigated and used to propose a new application of SEPs for identification of the location of pathological changes in the cervical spinal cord.

Electrophysiological and hemodynamic mismatch responses in rats listening to human speech syllables

Speech is a complex auditory stimulus which is processed according to several time-scales. Whereas consonant discrimination is required to resolve rapid acoustic events, voice perception relies on slower cues. Humans, right from preterm ages, are particularly efficient to encode temporal cues. To compare the capacities of preterms to those observed in other mammals, we tested anesthetized adult rats by using exactly the same paradigm as that used in preterm neonates. We simultaneously recorded neural (using ECoG) and hemodynamic responses (using fNIRS) to series of human speech syllables and investigated the brain response to a change of consonant (ba vs. ga) and to a change of voice (male vs. female). Both methods revealed concordant results, although ECoG measures were more sensitive than fNIRS. Responses to syllables were bilateral, but with marked right-hemispheric lateralization. Responses to voice changes were observed with both methods, while only ECoG was sensitive to consonant changes. These results suggest that rats more effectively processed the speech envelope than fine temporal cues in contrast with human preterm neonates, in whom the opposite effects were observed. Cross-species comparisons constitute a very valuable tool to define the singularities of the human brain and species-specific bias that may help human infants to learn their native language.

Shedding light on interictal epileptic spikes: An in vivo study using fast optical signal and electrocorticography

OBJECTIVE

Interictal epileptic spikes (IESs), apart from being a key marker of epileptic neuronal networks, constitute a nice model of the widespread endogenous phenomenon of neuronal hypersynchronization. Many questions concerning the mechanisms that drive neurons to hypersynchronize remain unresolved, but synaptic as well as nonsynaptic events are likely to be involved. In this study, changes in optical properties of neural tissues were observed in rats with penicillin-induced IES using fast optical signal (FOS) concomitantly with electrocorticography (ECoG).

METHODS

In this study, near-infrared optical imaging was used with ECoG to investigate variations in the optical properties of cortical tissue directly associated with neuronal activity in 15 rats. FOS changes correspond to variations of scattered light from neuronal tissue when neurons are activated. To independently evaluate our method, a control experiment on somatosensory was designed and applied to seven different rats. Time-frequency analysis was also used to track variations of (de)synchronization concomitantly with changes in optical signals during IES.

RESULTS

FOS responses revealed that changes in optical signals occurred 320 msec before to 370 msec after the IES peak. These changes started before any changes in ECoG signal. In addition, time-frequency domain electrocorticography revealed an alternating decrease-increase-decrease in the ECoG spectral power (pointing to desynchronization-synchronization-desynchronization), which occurred concomitantly with an increase-decrease-increase in relative optical signal during the IES. These results suggest a relationship between (de)synchronization and optical changes.

SIGNIFICANCE

These changes in the neuronal environment around IESs raise new questions about the mechanisms that induce changes in optical properties of neural tissues before the IES, which may provide suitable conditions for neuronal synchronization during IESs. FOS-ECoG constitutes a multimodal approach and opens new avenues to study the mechanisms of neuronal synchronization in the pathologic brain, which has clinical implications, at least in epilepsy.

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

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

The neurogenesis of P1 and N1: A concurrent EEG/LFP study

It is generally recognised that event related potentials (ERPs) of electroencephalogram (EEG) primarily reflect summed post-synaptic activity of the local pyramidal neural population(s). However, it is still not understood how the positive and negative deflections (e.g. P1, N1 etc) observed in ERP recordings are related to the underlying excitatory and inhibitory post-synaptic activity. We investigated the neurogenesis of P1 and N1 in ERPs by pharmacologically manipulating inhibitory post-synaptic activity in the somatosensory cortex of rodent, and concurrently recording EEG and local field potentials (LFPs). We found that the P1 wave in the ERP and LFP of the supragranular layers is determined solely by the excitatory post-synaptic activity of the local pyramidal neural population, as is the initial segment of the N1 wave across cortical depth. The later part of the N1 wave was modulated by inhibitory post-synaptic activity, with its peak and the pulse width increasing as inhibition was reduced. These findings suggest that the temporal delay of inhibition with respect to excitation observed in intracellular recordings is also reflected in extracellular field potentials (FPs), resulting in a temporal window during which only excitatory post-synaptic activity and leak channel activity are recorded in the ERP and evoked LFP time series. Based on these findings, we provide clarification on the interpretation of P1 and N1 in terms of the excitatory and inhibitory post-synaptic activities of the local pyramidal neural population(s).

Simultaneous epidural functional near-infrared spectroscopy and cortical electrophysiology as a tool for studying local neurovascular coupling in primates

Simultaneous measurements of intra-cortical electrophysiology and hemodynamic signals in primates are essential for relating human neuroimaging studies with intra-cortical electrophysiology in monkeys. Previously, technically challenging and resourcefully demanding techniques such as fMRI and intrinsic-signal optical imaging have been used for such studies. Functional near-infrared spectroscopy is a relatively less cumbersome neuroimaging method that uses near-infrared light to detect small changes in concentrations of oxy-hemoglobin (HbO), deoxy-hemoglobin (HbR) and total hemoglobin (HbT) in a volume of tissue with high specificity and temporal resolution. FNIRS is thus a good candidate for hemodynamic measurements in primates to acquire local hemodynamic signals during electrophysiological recordings. To test the feasibility of using epidural fNIRS with concomitant extracellular electrophysiology, we recorded neuronal and hemodynamic activity from the primary visual cortex of two anesthetized monkeys during visual stimulation. We recorded fNIRS epidurally, using one emitter and two detectors. We performed simultaneous cortical electrophysiology using tetrodes placed between the fNIRS sensors. We observed robust and reliable responses to the visual stimulation in both [HbO] and [HbR] signals, and quantified the signal-to-noise ratio of the epidurally measured signals. We also observed a positive correlation between stimulus-induced modulation of [HbO] and [HbR] signals and strength of neural modulation. Briefly, our results show that epidural fNIRS detects single-trial responses to visual stimuli on a trial-by-trial basis, and when coupled with cortical electrophysiology, is a promising tool for studying local hemodynamic signals and neurovascular coupling.

Long-latency reductions in gamma power predict hemodynamic changes that underlie the negative BOLD signal

Studies that use prolonged periods of sensory stimulation report associations between regional reductions in neural activity and negative blood oxygenation level-dependent (BOLD) signaling. However, the neural generators of the negative BOLD response remain to be characterized. Here, we use single-impulse electrical stimulation of the whisker pad in the anesthetized rat to identify components of the neural response that are related to "negative" hemodynamic changes in the brain. Laminar multiunit activity and local field potential recordings of neural activity were performed concurrently with two-dimensional optical imaging spectroscopy measuring hemodynamic changes. Repeated measurements over multiple stimulation trials revealed significant variations in neural responses across session and animal datasets. Within this variation, we found robust long-latency decreases (300 and 2000 ms after stimulus presentation) in gamma-band power (30-80 Hz) in the middle-superficial cortical layers in regions surrounding the activated whisker barrel cortex. This reduction in gamma frequency activity was associated with corresponding decreases in the hemodynamic responses that drive the negative BOLD signal. These findings suggest a close relationship between BOLD responses and neural events that operate over time scales that outlast the initiating sensory stimulus, and provide important insights into the neurophysiological basis of negative neuroimaging signals.

Contributions and complexities from the use of in vivo animal models to improve understanding of human neuroimaging signals

Many of the major advances in our understanding of how functional brain imaging signals relate to neuronal activity over the previous two decades have arisen from physiological research studies involving experimental animal models. This approach has been successful partly because it provides opportunities to measure both the hemodynamic changes that underpin many human functional brain imaging techniques and the neuronal activity about which we wish to make inferences. Although research into the coupling of neuronal and hemodynamic responses using animal models has provided a general validation of the correspondence of neuroimaging signals to specific types of neuronal activity, it is also highlighting the key complexities and uncertainties in estimating neural signals from hemodynamic markers. This review will detail how research in animal models is contributing to our rapidly evolving understanding of what human neuroimaging techniques tell us about neuronal activity. It will highlight emerging issues in the interpretation of neuroimaging data that arise from in vivo research studies, for example spatial and temporal constraints to neuroimaging signal interpretation, or the effects of disease and modulatory neurotransmitters upon neurovascular coupling. We will also give critical consideration to the limitations and possible complexities of translating data acquired in the typical animals models used in this area to the arena of human fMRI. These include the commonplace use of anesthesia in animal research studies and the fact that many neuropsychological questions that are being actively explored in humans have limited homologs within current animal models for neuroimaging research. Finally we will highlighting approaches, both in experimental animals models (e.g. imaging in conscious, behaving animals) and human studies (e.g. combined fMRI-EEG), that mitigate against these challenges.

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

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

GABAA Receptor-Mediated Bidirectional Control of Synaptic Activity, Intracellular Ca2+, Cerebral Blood Flow, and Oxygen Consumption in Mouse Somatosensory Cortex In Vivo

Neural activity regulates local increases in cerebral blood flow (ΔCBF) and the cortical metabolic rate of oxygen (ΔCMRO2) that constitutes the basis of BOLD functional neuroimaging signals. Glutamate signaling plays a key role in brain vascular and metabolic control; however, the modulatory effect of GABA is incompletely understood. Here we performed in vivo studies in mice to investigate how THIP (which tonically activates extrasynaptic GABAARs) and Zolpidem (a positive allosteric modulator of synaptic GABAARs) impact stimulation-induced ΔCBF, ΔCMRO2, local field potentials (LFPs), and fluorescent cytosolic Ca(2+) transients in neurons and astrocytes. Low concentrations of THIP increased ΔCBF and ΔCMRO2 at low stimulation frequencies. These responses were coupled to increased synaptic activity as indicated by LFP responses, and to Ca(2+) activities in neurons and astrocytes. Intermediate and high concentrations of THIP suppressed ΔCBF and ΔCMRO2 at high stimulation frequencies. Zolpidem had similar but less-pronounced effects, with similar dependence on drug concentration and stimulation frequency. Our present findings suggest that slight increases in both synaptic and extrasynaptic GABAAR activity might selectively gate and amplify transient low-frequency somatosensory inputs, filter out high-frequency inputs, and enhance vascular and metabolic responses that are likely to be reflected in BOLD functional neuroimaging signals.

Neurovascular coupling in normal aging: a combined optical, ERP and fMRI study

Brain aging is characterized by changes in both hemodynamic and neuronal responses, which may be influenced by the cardiorespiratory fitness of the individual. To investigate the relationship between neuronal and hemodynamic changes, we studied the brain activity elicited by visual stimulation (checkerboard reversals at different frequencies) in younger adults and in older adults varying in physical fitness. Four functional brain measures were used to compare neuronal and hemodynamic responses obtained from BA17: two reflecting neuronal activity (the event-related optical signal, EROS, and the C1 response of the ERP), and two reflecting functional hemodynamic changes (functional magnetic resonance imaging, fMRI, and near-infrared spectroscopy, NIRS). The results indicated that both younger and older adults exhibited a quadratic relationship between neuronal and hemodynamic effects, with reduced increases of the hemodynamic response at high levels of neuronal activity. Although older adults showed reduced activation, similar neurovascular coupling functions were observed in the two age groups when fMRI and deoxy-hemoglobin measures were used. However, the coupling between oxy- and deoxy-hemoglobin changes decreased with age and increased with increasing fitness. These data indicate that departures from linearity in neurovascular coupling may be present when using hemodynamic measures to study neuronal function.

Functional electrical stimulation post-spinal cord injury improves locomotion and increases afferent input into the central nervous system in rats

BACKGROUND

Functional electrical stimulation (FES) has been found to be effective in restoring voluntary functions after spinal cord injury (SCI) and stroke. However, the central nervous system (CNS) changes that occur in as a result of this therapy are largely unknown.

OBJECTIVE

To examine the effects of FES on the restoration of voluntary locomotor function of the CNS in a SCI rat model.

METHODS

SCI rats were instrumented with chronic FES electrodes in the hindlimb muscles and were divided into two groups: (a) FES therapy and (b) sedentary. At day 7 post-SCI, the animals were assessed for locomotion performance by using a Basso, Beattie and Bresnahan (BBB) scale. They were then anesthetized for a terminal in vivo experiment. The lumbar spinal cord and somatosensory cortex were exposed and the instrumented muscles were stimulated electrically. Associated neurovascular responses in the CNS were recorded with an intrinsic optical imaging system.

RESULTS

FES greatly improved locomotion recovery by day 7 post-SCI, as measured by BBB scores (P < 0.05): (a) FES 10 ± 2 and (b) controls 3 ± 1. Furthermore, the FES group showed a significant increase (P < 0.05) of neurovascular activation in the spinal cord and somatosensory cortex when the muscles were stimulated between 1 and 3 motor threshold (MT).

CONCLUSION

Hind limb rehabilitation with FES is an effective strategy to improve locomotion during the acute phase post-SCI. The results of this study indicate that after FES, the CNS preserves/acquires the capacity to respond to peripheral electrical stimulation.

Neurovascular signaling in the brain and the pathological consequences of hypertension

The execution and maintenance of all brain functions are dependent on a continuous flow of blood to meet the metabolic needs of the tissue. To ensure the delivery of resources required for neural processing and the maintenance of neural homeostasis, the cerebral vasculature is elaborately and extensively regulated by signaling from neurons, glia, interneurons, and perivascular nerves. Hypertension is associated with impaired neurovascular regulation of the cerebral circulation and culminates in neurodegeneration and cognitive dysfunction. Here, we review the physiological processes of neurovascular signaling in the brain and discuss mechanisms of hypertensive neurovascular dysfunction.

Neurovascular coupling: in vivo optical techniques for functional brain imaging

Optical imaging techniques reflect different biochemical processes in the brain, which is closely related with neural activity. Scientists and clinicians employ a variety of optical imaging technologies to visualize and study the relationship between neurons, glial cells and blood vessels. In this paper, we present an overview of the current optical approaches used for the in vivo imaging of neurovascular coupling events in small animal models. These techniques include 2-photon microscopy, laser speckle contrast imaging (LSCI), voltage-sensitive dye imaging (VSDi), functional photoacoustic microscopy (fPAM), functional near-infrared spectroscopy imaging (fNIRS) and multimodal imaging techniques. The basic principles of each technique are described in detail, followed by examples of current applications from cutting-edge studies of cerebral neurovascular coupling functions and metabolic. Moreover, we provide a glimpse of the possible ways in which these techniques might be translated to human studies for clinical investigations of pathophysiology and disease. In vivo optical imaging techniques continue to expand and evolve, allowing us to discover fundamental basis of neurovascular coupling roles in cerebral physiology and pathophysiology.

Compliant head probe for positioning electroencephalography electrodes and near-infrared spectroscopy optodes

A noninvasive head probe that combines near-infrared spectroscopy (NIRS) and electroencephalography (EEG) for simultaneous measurement of neural dynamics and hemodynamics in the brain is presented. It is composed of a compliant expandable mechanism that accommodates a wide range of head size variation and an elastomeric web that maintains uniform sensor contact pressure on the scalp as the mechanism expands and contracts. The design is intended to help maximize optical and electrical coupling and to maintain stability during head movement. Positioning electrodes at the inion, nasion, central, and preauricular fiducial locations mechanically shapes the probe to place 64 NIRS optodes and 65 EEG electrodes following the 10-5 scalp coordinates. The placement accuracy, precision, and scalp pressure uniformity of the sensors are evaluated. A root-mean-squared (RMS) positional precision of 0.89 ± 0.23 mm, percent arc subdivision RMS accuracy of 0.19 ± 0.15%, and mean normal force on the scalp of 2.28 ± 0.88 N at 5 mm displacement were found. Geometric measurements indicate that the probe will accommodate the full range of adult head sizes. The placement accuracy, precision, and uniformity of sensor contact pressure of the proposed head probe are important determinants of data quality in noninvasive brain monitoring with simultaneous NIRS-EEG.

Neurovascular saturation thresholds under high intensity auditory stimulation during wake

Coupling between neural activity and hemodynamic responses is important in understanding brain function, interpreting brain-imaging signals, and assessing pathological conditions. Tissue state is a major factor in neurovascular coupling and may alter the relationship between neural and hemodynamic activity. However, most neurovascular-coupling studies are performed under anesthetized or sedated states which may have severe consequences on coupling mechanisms. Our previous studies showed that following prolonged periods of sleep deprivation, evoked hemodynamic responses were muted despite consistent electrical responses, suggesting that sustained neural activity may decrease vascular compliance and limit blood perfusion. To investigate potential perfusion limitations during natural waking conditions, we simultaneously measured evoked response potentials (ERPs) and evoked hemodynamic responses using optical-imaging techniques to increase intensity auditory stimulation. The relationship between evoked hemodynamic responses and integrated ERPs followed a sigmoid relationship where the hemodynamic response approached saturation at lower stimulus intensities than the ERP. If limits in blood perfusion are caused by stretching of the vessel wall, then these results suggest there may be decreased vascular compliance due to sustained neural activity during wake, which could limit vascular responsiveness and local blood perfusion. Conditions that stress cerebral vasculature, such as sleep deprivation and some pathologies (e.g., epilepsy), may further decrease vascular compliance, limit metabolic delivery, and cause tissue trauma. While ERPs and evoked hemodynamic responses provide an indication of the correlated neural activity and metabolic demand, the relationship between these two responses is complex and the different measurement techniques are not directly correlated. Future studies are required to verify these findings and further explore neurovascular coupling during wake by assessing local field potentials, vascular expansion, hemodynamic response localization.

Neurovascular uncoupling under mild hypoxic hypoxia: an EEG-fMRI study in rats

The effects of oxygen availability on neurovascular coupling were investigated using simultaneous electroencephalography (EEG) and functional magnetic resonance imaging (fMRI), in addition to the monitoring of physiological parameters, in 16 α-chloralose-anesthetized rats. Mild hypoxic hypoxia (oxygen saturation=83.6±12.1%) induced significant reductions in fMRI responses (P<0.05) to electrical stimulation in the forepaw, but EEG responses remained unchanged. In addition, the changes in oxygen saturation were linearly correlated with the changes in the fMRI responses. These data further emphasize the importance of oxygen availability, which may regulate neurovascular coupling via the oxygen-dependent enzymatic synthesis of messenger molecules.

Anesthesia and the quantitative evaluation of neurovascular coupling

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

Usefulness of simultaneous EEG-NIRS recording in language studies

One of the most challenging tasks in neuroscience in language studies, is investigation of the brain's ability to integrate and process information. This task can only be successfully addressed by applying various assessment techniques integrated into a multimodal approach. Each of these techniques has its advantages and disadvantages, but help to elucidate certain aspects of the capacity of neural networks to process information. These methods provide information about changes in electrical, hemodynamic and metabolic activities. Ideally, they should be noninvasive in order to facilitate their use particularly in children. In the present review, we describe the advantages of simultaneous electroencephalographic (EEG) acquisition with near infrared spectroscopy (NIRS) to provide a better understanding of the mechanisms involved in cerebral activation. This coregistration is also useful to avoid misleading interpretation of NIRS, notably during the various phases of sleep. Development and implementation of the various tools required and assessment strategies are also discussed.

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

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

Evoked electrical and cerebral vascular responses during sleep and following sleep deprivation

Neuronal activity elicits vascular dilation, delivering additional blood and metabolites to the activated region. With increasing neural activity, vessels stretch and may become less compliant. Most functional imaging studies assume that limits to vascular expansion are not normally reached except under pathological conditions, with the possibility that metabolism could outpace supply. However, we previously demonstrated that evoked hemodynamic responses were larger during quiet sleep when compared to both waking and rapid eye movement (REM) sleep, suggesting that high basal activity during wake may elicit blunted evoked hemodynamic responses due to vascular expansion limits. We hypothesized that extended brain activity through sleep deprivation will further dilate blood vessels and exacerbate the blunted evoked hemodynamic responses observed during wake, and dampen responses in subsequent sleep. We measured evoked electrical and hemodynamic responses from rats using auditory clicks (0.5s, 10 Hz, 2-13s random ISIs) for 1h following 2, 4, or 6h of sleep deprivation. Time-of-day matched controls were recorded continuously for 7h. Within quiet sleep periods following deprivation, evoked response potential (ERP) amplitude did not differ; however, the evoked vascular response was smaller with longer sleep deprivation periods. These results suggest that prolonged neural activity periods through sleep deprivation may diminish vascular compliance as indicated by the blunted vascular response. Subsequent sleep may allow vessels to relax, restoring their ability to deliver blood. These results also suggest that severe sleep deprivation or chronic sleep disturbances could push the vasculature to critical limits, leading to metabolic deficit and the potential for tissue trauma.

Modeling habituation in rat EEG-evoked responses via a neural mass model with feedback

Habituation is a generic property of the neural response to repeated stimuli. Its strength often increases as inter-stimuli relaxation periods decrease. We propose a simple, broadly applicable control structure that enables a neural mass model of the evoked EEG response to exhibit habituated behavior. A key motivation for this investigation is the ongoing effort to develop model-based reconstruction of multi-modal functional neuroimaging data. The control structure proposed here is illustrated and validated in the context of a biophysical neural mass model, developed by Riera et al. (Hum Brain Mapp 27(11):896-914, 2006; 28(4):335-354, 2007), and of simplifications thereof, using data from rat EEG response to medial nerve stimuli presented at frequencies from 1 to 8 Hz. Performance was tested by predictions of both the response to the next stimulus based on the current one, and also of continued stimuli trains over 4-s time intervals based on the first stimulus in the interval, with similar success statistics. These tests demonstrate the ability of simple generative models to capture key features of the evoked response, including habituation.

Biophysical model estimation of neurovascular parameters in a rat model of healthy aging

Neuronal, vascular and metabolic factors result in a deterioration of the cerebral hemodynamic response with age. The interpretation of neuroimaging studies in the context of aging is rendered difficult due to the challenge in untangling the composite effect of these modifications. In this work we integrate multimodal optical imaging in biophysical models to investigate vascular and metabolic changes occurring in aging. Multispectral intrinsic optical imaging of an animal model of healthy aging, the LOU/c rat, is used in combination with somatosensory stimulation to study the modifications of the hemodynamic response with increasing age. Results are fitted with three macroscopic biophysical models to extract parameters, providing a phenomenological description of vascular and metabolic changes. Our results show that 1) biophysical parameters are estimable from multimodal data and 2) parameter estimates in this population change with aging.

Cerebral monitoring during carotid endarterectomy using near-infrared diffuse optical spectroscopies and electroencephalogram

Intraoperative monitoring of cerebral hemodynamics during carotid endarterectomy (CEA) provides essential information for detecting cerebral hypoperfusion induced by temporary internal carotid artery (ICA) clamping and post-CEA hyperperfusion syndrome. This study tests the feasibility and sensitivity of a novel dual-wavelength near-infrared diffuse correlation spectroscopy technique in detecting cerebral blood flow (CBF) and cerebral oxygenation in patients undergoing CEA. Two fiber-optic probes were taped on both sides of the forehead for cerebral hemodynamic measurements, and the instantaneous decreases in CBF and electroencephalogram (EEG) alpha-band power during ICA clamping were compared to test the measurement sensitivities of the two techniques. The ICA clamps resulted in significant CBF decreases (-24.7 ± 7.3%) accompanied with cerebral deoxygenation at the surgical sides (n = 12). The post-CEA CBF were significantly higher (+43.2 ± 16.9%) than the pre-CEA CBF. The CBF responses to ICA clamping were significantly faster, larger and more sensitive than EEG responses. Simultaneous monitoring of CBF, cerebral oxygenation and EEG power provides a comprehensive evaluation of cerebral physiological status, thus showing potential for the adoption of acute interventions (e.g., shunting, medications) during CEA to reduce the risks of severe cerebral ischemia and cerebral hyperperfusion syndrome.

Real-time in vivo simultaneous measurements of nitric oxide and oxygen using an amperometric dual microsensor

This paper reports a real-time study of the codynamical changes in the release of endogenous nitric oxide (NO) and oxygen (O(2)) consumption in a rat neocortex in vivo upon electrical stimulation using an amperometric NO/O(2) dual microsensor. Electrical stimulation induced transient cerebral hypoxia due to the increased metabolic demands that were not met by the blood volume inside the stimulated cortical region. A NO/O(2) dual microsensor was successfully used to monitor the pair of real-time dynamic changes in the tissue NO and O(2) contents. At the onset of electrical stimulation, there was an immediate decrease in the cortical tissue O(2) followed by a subsequent increase in the cortical tissue NO content. The averages of the maximum normalized concentration changes induced by the stimulation were a 0.41 (±0.04)-fold decrease in the O(2) and a 3.6 (±0.9)-fold increase in the NO concentrations when compared with the corresponding normalized basal levels. The peak increase in NO was always preceded by the peak decrease in O(2) in all animals (n = 11). The delay between the maximum decrease in O(2) and the maximum increase in NO varied from 3.1 to 54.8 s. This rather wide variation in the temporal associations was presumably attributed to the sparse distribution of NOS-containing neurons and the individual animal's differences in brain vasculatures, which suggests that a sensor with fine spatial resolution is needed to measure the location-specific real-time NO and O(2) contents. In summary, the developed NO/O(2) dual microsensor is effective for measuring the NO and O(2) contents in vivo. This study provides direct support for the dynamic role of NO in regulating the cerebral hemodynamics, particularly related to the tissue oxygenation.

Study of neurovascular coupling by modulating neuronal activity with GABA

Fundamental to the interpretation of neurovascular coupling is determining the neuronal activity that accounts for functional hyperemia. Recently, synaptic and not spiking activity has been found to be responsible for the hemodynamic response. Using pharmacological manipulation in rats, we want to further determine whether the cortical synaptic activity generated by the thalamic input or the subsequent synaptic activity related to secondary cortical processing is driving the hemodynamic response. In this study, we topically applied γ-aminobutyric acid (GABA) in the somatosensory cortex and used electrical forepaw stimulation to evoke neural and vascular activity. In a group of 8 animals, using laminar electrophysiology, we verified that topical application of GABA for 20min does not affect layer IV synaptic activity but reduces subsequent activity in the supragranular and infragranular layers. In another group of 8 animals, we simultaneously measured the electrical and vascular responses with scalp electroencephalography (EEG) and diffuse optical imaging (DOI), respectively. We decomposed somatosensory evoked potentials (SEP) into three major components: P1, N1, and P2, where P1 represents the thalamic input activity originating in layer IV and N1 and P2 represent the subsequent cortical transmissions. We verified that GABA infusion in the somatosensory cortex does not significantly reduce the P1 SEP component but strongly reduces the N1 and P2 components. We found that GABA also elicits a large reduction in the hemodynamic responses, which correlate with the reduction in N1 and P2 components. These results suggest that the hemodynamic response is predominantly driven by cortico-cortical interactions and not by the initial thalamocortical activity in layer IV.

Does neural input or processing play a greater role in the magnitude of neuroimaging signals?

An important constraint on how hemodynamic neuroimaging signals such as fMRI can be interpreted in terms of the underlying evoked activity is an understanding of neurovascular coupling mechanisms that actually generate hemodynamic responses. The predominant view at present is that the hemodynamic response is most correlated with synaptic input and subsequent neural processing rather than spiking output. It is still not clear whether input or processing is more important in the generation of hemodynamics responses. In order to investigate this we measured the hemodynamic and neural responses to electrical whisker pad stimuli in rat whisker barrel somatosensory cortex both before and after the local cortical injections of the GABA(A) agonist muscimol. Muscimol would not be expected to affect the thalamocortical input into the cortex but would inhibit subsequent intra-cortical processing. Pre-muscimol infusion whisker stimuli elicited the expected neural and accompanying hemodynamic responses to that reported previously. Following infusion of muscimol, although the temporal profile of neural responses to each pulse of the stimulus train was similar, the average response was reduced in magnitude by approximately 79% compared to that elicited pre-infusion. The whisker-evoked hemodynamic responses were reduced by a commensurate magnitude suggesting that, although the neurovascular coupling relationships were similar for synaptic input as well as for cortical processing, the magnitude of the overall response is dominated by processing rather than from that produced from the thalamocortical input alone.

Effects of blood glutamate scavenging on cortical evoked potentials

It is well known that traumatic or ischemic brain injury is followed by acute excitotoxicity caused by the presence of abnormally high glutamate (Glu) in brain fluids. It has recently been demonstrated that excess Glu can be eliminated from brain into blood following the intravenous administration of oxaloacetate (OxAc), which, by scavenging blood Glu, induces an enhanced and neuroprotective brain-to-blood Glu efflux. In this study, we subjected rats to intravenous OxAc administration (i.v., 12.5, 25, and 50 mg/kg, respectively), and studied its effects on somatosensory evoked cortical potentials (EPs). Against our expectation, the amplitudes of EPs did not decrease but increased in a dose- and time-dependent manner after OxAc administration. Similar effects were observed when blood Glu scavenging was enhanced by combining OxAc (12.5 mg/kgbw) with recombinant glutamate-oxaloacetate transaminase (GOT, 0.14 nmol/100 g rat). On the basis of these results, we suggest that the changes of amplitudes of the EPs involve not only a glutamatergic but also the weakening of a GABAergic component. We cannot rule out the possibility that OxAc penetrates into the brain and improves mitochondrial functions.

Revisiting the role of neurons in neurovascular coupling

In this article, we will review molecular, anatomical, physiological and pharmacological data in an attempt to better understand how excitatory and inhibitory neurons recruited by distinct afferent inputs to the cerebral cortex contribute to the coupled hemodynamic response, and how astrocytes can act as intermediaries to these neuronal populations. We aim at providing the pros and cons to the following statements that, depending on the nature of the afferent input to the neocortex, (i) different neuronal or astroglial messengers, likely acting in sequence, mediate the hemodynamic changes, (ii) some recruited neurons release messengers that directly alter blood vessel tone, (iii) others act by modulating neuronal and astroglial activity, and (iv) astrocytes act as intermediaries for both excitatory and inhibitory neurotransmitters. We will stress that a given afferent signal activates a precise neuronal circuitry that determines the mediators of the hemodynamic response as well as the level of interaction with surrounding astrocytes.

The effect of different anesthetics on neurovascular coupling

To date, the majority of neurovascular coupling studies focused on the thalamic afferents' activity in layer IV and the corresponding large spiking activity as responsible for functional hyperemia. This paper highlights the role of the secondary and late cortico-cortical transmission in neurovascular coupling. Simultaneous scalp electroencephalography (EEG) and diffuse optical imaging (DOI) measurements were obtained during multiple conditions of event-related electrical forepaw stimulation in 33 male Sprague-Dawley rats divided into 6 groups depending on the maintaining anesthetic - alpha-chloralose, pentobarbital, ketamine-xylazine, fentanyl-droperidol, isoflurane, or propofol. The somatosensory evoked potentials (SEP) were decomposed into four components and the question of which best predicts the hemodynamic responses was investigated. Results of the linear regression analysis show that the hemodynamic response is best correlated with the secondary and late cortico-cortical transmissions and not with the initial thalamic input activity in layer IV. Baseline cerebral blood flow (CBF) interacts with neural activity and influences the evoked hemodynamic responses. Finally, neurovascular coupling appears to be the same across all anesthetics used.

BOLD responses in human auditory cortex are more closely related to transient MEG responses than to sustained ones

Blood oxygen level dependent-functional magnetic resonance imaging (BOLD-fMRI) and magnetoencephalographic (MEG) signals are both coupled to postsynaptic potentials, although their relationship is incompletely understood. Here, the wide range of BOLD-fMRI and MEG responses produced by auditory cortex was exploited to better understand the BOLD-fMRI/MEG relationship. Measurements of BOLD and MEG responses were made in the same subjects using the same stimuli for both modalities. The stimuli, 24-s sequences of click trains, had duty cycles of 2.5, 25, 72, and 100%. For the 2.5% sequence, the BOLD response was elevated throughout the sequence, whereas for 100%, it peaked after sequence onset and offset and showed a diminished elevation in between. On the finer timescale of MEG, responses at 2.5% consisted of a complex of transients, including N(1)m, to each click train of the sequence, whereas for 100% the only transients occurred at sequence onset and offset between which there was a sustained elevation in the MEG signal (a sustained field). A model that separately estimated the contributions of transient and sustained MEG signals to the BOLD response best fit BOLD measurements when the transient contribution was weighted 8- to 10-fold more than the sustained one. The findings suggest that BOLD responses in the auditory cortex are tightly coupled to the neural activity underlying transient, not sustained, MEG signals.

Linear and nonlinear relationships between visual stimuli, EEG and BOLD fMRI signals

In the present study, the cascaded interactions between stimuli and neural and hemodynamic responses were modeled using linear systems. These models provided the theoretical hypotheses that were tested against the electroencephalography (EEG) and blood oxygen level dependent (BOLD) functional magnetic resonance imaging (fMRI) data recorded from human subjects during prolonged periods of repeated visual stimuli with a variable setting of the inter-stimulus interval (ISI) and visual contrast. Our results suggest that (1) neural response is nonlinear only when ISI<0.2 s, (2) BOLD response is nonlinear with an exclusively vascular origin when 0.25<ISI<4.2 s, (3) vascular response nonlinearity reflects a refractory effect, rather than a ceiling effect, and (4) there is a strong linear relationship between the BOLD effect size and the integrated power of event-related synaptic current activity, after modeling and taking into account the vascular refractory effect. These conclusions offer important insights into the origins of BOLD nonlinearity and the nature of neurovascular coupling, and suggest an effective means to quantitatively interpret the BOLD signal in terms of neural activity. The validated cross-modal relationship between fMRI and EEG may provide a theoretical basis for the integration of these two modalities.

Predicting stimulus-rate sensitivity of human somatosensory fMRI signals with MEG

With increasing stimulus rate (SR), cortical EEG and MEG responses typically decrease in amplitude whereas BOLD fMRI signals increase. To address this discrepancy, we predicted BOLD responses with squared MEG waveforms using a recently proposed energy-density model. Tactile stimuli were delivered to finger tips at SRs of 1, 4, or 10 Hz in successive 25-s blocks, and brain signals were detected from area 3b of the primary somatosensory cortex of nine healthy adults using a 306-channel whole-scalp neuromagnetometer and a 3-T fMRI magnet. The main MEG deflections decreased in amplitude as a function of SR, whereas the BOLD signals increased from 1- to 4-Hz SR, with no further change at 10 Hz. MEG energy densities, obtained over the whole stimulus train and convolved with different hemodynamic response functions, predicted both the shape and amplitude of the BOLD signals well, and incorporation of nonlinear terms into the model did not offer any further advantage. Thus, squared MEG waveforms obtained over the entire stimulus train provided an appropriate estimate of area 3b neuronal activity associated with the BOLD signal.