Noninvasive measurement of neuronal activity with near-infrared optical imaging

Diffuse optical tomography of the brain: effects of inaccurate baseline optical parameters and refinements using learned post-processing

Diffuse optical tomography (DOT) uses near-infrared light to image spatially varying optical parameters in biological tissues. In functional brain imaging, DOT uses a perturbation model to estimate the changes in optical parameters, corresponding to changes in measured data due to brain activity. The perturbation model typically uses approximate baseline optical parameters of the different brain compartments, since the actual baseline optical parameters are unknown. We simulated the effects of these approximate baseline optical parameters using parameter variations earlier reported in literature, and brain atlases from four adult subjects. We report the errors in estimated activation contrast, localization, and area when incorrect baseline values were used. Further, we developed a post-processing technique based on deep learning methods that can reduce the effects due to inaccurate baseline optical parameters. The method improved imaging of brain activation changes in the presence of such errors.

Low frequency oscillations reflect neurovascular coupling and disappear after cerebral death

Spectrum power analysis in the low frequency oscillations (LFO) region of functional near infrared spectroscopy (fNIRS) is a promising method to deliver information about brain activation and therefore might be used for prognostication in patients with disorders of consciousness in the neurocritical care unit alongside with established methods. In this study, we measure the cortical hemodynamic response measured by fNIRS in the LFO region following auditory and somatosensory stimulation in healthy subjects. The significant hemodynamic reaction in the contralateral hemisphere correlation with the physiologic electric response suggests neurovascular coupling. In addition, we investigate power spectrum changes in steady state measurements of cerebral death patients and healthy subjects in the LFO region, the frequency of the heartbeat and respiration. The spectral power within the LFO region was lower in the patients with cerebral death compared to the healthy subjects, whereas there were no differences in spectral power for physiological activities such as heartbeat and respiration rate. This finding indicates the cerebral origin of our low frequency measurements. Therefore, LFO measurements are a potential method to detect brain activation in patients with disorders of consciousness and cerebral death. However, further studies in patients are needed to investigate its potential clinical use.

Fast Optical Signals for Real-Time Retinotopy and Brain Computer Interface

A brain-computer interface (BCI) allows users to control external devices through brain activity. Portable neuroimaging techniques, such as near-infrared (NIR) imaging, are suitable for this goal. NIR imaging has been used to measure rapid changes in brain optical properties associated with neuronal activation, namely fast optical signals (FOS) with good spatiotemporal resolution. However, FOS have a low signal-to-noise ratio, limiting their BCI application. Here FOS were acquired with a frequency-domain optical system from the visual cortex during visual stimulation consisting of a rotating checkerboard wedge, flickering at 5 Hz. We used measures of photon count (Direct Current, DC light intensity) and time of flight (phase) at two NIR wavelengths (690 nm and 830 nm) combined with a machine learning approach for fast estimation of visual-field quadrant stimulation. The input features of a cross-validated support vector machine classifier were computed as the average modulus of the wavelet coherence between each channel and the average response among all channels in 512 ms time windows. An above chance performance was obtained when differentiating visual stimulation quadrants (left vs. right or top vs. bottom) with the best classification accuracy of ~63% (information transfer rate of ~6 bits/min) when classifying the superior and inferior stimulation quadrants using DC at 830 nm. The method is the first attempt to provide generalizable retinotopy classification relying on FOS, paving the way for the use of FOS in real-time BCI.

Multimodal assessment of the spatial correspondence between fNIRS and fMRI hemodynamic responses in motor tasks

Functional near-infrared spectroscopy (fNIRS) provides a cost-efficient and portable alternative to functional magnetic resonance imaging (fMRI) for assessing cortical activity changes based on hemodynamic signals. The spatial and temporal underpinnings of the fMRI blood-oxygen-level-dependent (BOLD) signal and corresponding fNIRS concentration of oxygenated (HbO), deoxygenated (HbR), and total hemoglobin (HbT) measurements are still not completely clear. We aim to analyze the spatial correspondence between these hemodynamic signals, in motor-network regions. To this end, we acquired asynchronous fMRI and fNIRS recordings from 9 healthy participants while performing motor imagery and execution. Using this multimodal approach, we investigated the ability to identify motor-related activation clusters in fMRI data using subject-specific fNIRS-based cortical signals as predictors of interest. Group-level activation was found in fMRI data modeled from corresponding fNIRS measurements, with significant peak activation found overlapping the individually-defined primary and premotor motor cortices, for all chromophores. No statistically significant differences were observed in multimodal spatial correspondence between HbO, HbR, and HbT, for both tasks. This suggests the possibility of translating neuronal information from fMRI into an fNIRS motor-coverage setup with high spatial correspondence using both oxy and deoxyhemoglobin data, with the inherent benefits of translating fMRI paradigms to fNIRS in cognitive and clinical neuroscience.

Comparison of cerebral metabolic rate of oxygen, blood flow, and bispectral index under general anesthesia

Significance

The optical measurement of cerebral oxygen metabolism was evaluated.

Aim

Compare optically derived cerebral signals to the electroencephalographic bispectral index (BIS) sensors to monitor propofol-induced anesthesia during surgery.

Approach

Relative cerebral metabolic rate of oxygen (rCMRO 2 ) and blood flow (rCBF) were measured by time-resolved and diffuse correlation spectroscopies. Changes were tested against the relative BIS (rBIS) ones. The synchronism in the changes was also assessed by the R-Pearson correlation.

Results

In 23 measurements, optically derived signals showed significant changes in agreement with rBIS: during propofol induction, rBIS decreased by 67% [interquartile ranges (IQR) 62% to 71%],rCMRO 2 by 33% (IQR 18% to 46%), and rCBF by 28% (IQR 10% to 37%). During recovery, a significant increase was observed for rBIS (48%, IQR 38% to 55%),rCMRO 2 (29%, IQR 17% to 39%), and rCBF (30%, IQR 10% to 44%). The significance and direction of the changes subject-by-subject were tested: the coupling between the rBIS,rCMRO 2 , and rCBF was witnessed in the majority of the cases (14/18 and 12/18 for rCBF and 19/21 and 13/18 forrCMRO 2 in the initial and final part, respectively). These changes were also correlated in time ( R > 0.69 to R = 1 , p - values < 0.05 ).

Conclusions

Optics can reliably monitorrCMRO 2 in such conditions.

Frontal lobe hemodynamics detected by functional near-infrared spectroscopy during head-up tilt table tests in patients with electrical burns

Significance

Electrical burns can cause severe damage to the nervous system, resulting in autonomic dysfunction with reduced cerebral perfusion. However, few studies have investigated these consequences.

Aim

To elucidate changes in prefrontal cerebral hemodynamics using functional near-infrared spectroscopy (fNIRS) during the head-up tilt table test (HUT) for patients with electrical burns.

Approach

We recruited 17 patients with acute electrical burns within 1 week after their accidents and 10 healthy volunteers. The NIRS parameters acquired using an fNIRS device attached to the forehead were analyzed in five distinct HUT phases.

Results

Based on their HUT response patterns, patients with electrical burns were classified into the group with abnormal HUT results (APG, n = 4) or normal HUT results (NPG, n = 13) and compared with the healthy control (HC, n = 10) participants. We found trends in hemodynamic changes during the HUT that distinguished HC, NPG, and APG. Reduced cerebral perfusion and decreased blood oxygenation during the HUT were found in both the NPG and APG groups. Patients with electrical burns had autonomic dysfunction compared to the HC participants.

Conclusions

Using fNIRS, we observed that acute-stage electrical burn injuries could affect cerebral perfusion.

Prefrontal fNIRS-based clinical data analysis of brain functions in individuals abusing different types of drugs

Background

The activation degree of the orbitofrontal cortex (OFC) functional area in drug abusers is directly related to the craving for drugs and the tolerance to punishment. Currently, among the clinical research on drug rehabilitation, there has been little analysis of the OFC activation in individuals abusing different types of drugs, including heroin, methamphetamine, and mixed drugs. Therefore, it becomes urgently necessary to clinically investigate the abuse of different drugs, so as to explore the effects of different types of drugs on the human brain.

Methods

Based on prefrontal high-density functional near-infrared spectroscopy (fNIRS), this research designs an experiment that includes resting and drug addiction induction. Hemoglobin concentrations of 30 drug users (10 on methamphetamine, 10 on heroin, and 10 on mixed drugs) were collected using fNIRS and analyzed by combining algorithm and statistics.

Results

Linear discriminant analysis (LDA), Support vector machine (SVM) and Machine-learning algorithm was implemented to classify different drug abusers. Oxygenated hemoglobin (HbO2) activations in the OFC of different drug abusers were statistically analyzed, and the differences were confirmed. Innovative findings: in both the Right-OFC and Left-OFC areas, methamphetamine abusers had the highest degree of OFC activation, followed by those abusing mixed drugs, and heroin abusers had the lowest. The same result was obtained when OFC activation was investigated without distinguishing the left and right hemispheres.

Conclusions

The findings confirmed the significant differences among different drug abusers and the patterns of OFC activations, providing a theoretical basis for personalized clinical treatment of drug rehabilitation in the future.

Measuring neuronal activity with diffuse correlation spectroscopy: a theoretical investigation

Significance: Diffuse correlation spectroscopy (DCS) measures cerebral blood flow non-invasively. Variations in blood flow can be used to detect neuronal activities, but its peak has a latency of a few seconds, which is slow for real-time monitoring. Neuronal cells also deform during activation, which, in principle, can be utilized to detect neuronal activity on fast timescales (within 100 ms) using DCS. Aims: We aim to characterize DCS signal variation quantified as the change of the decay time of the speckle intensity autocorrelation function during neuronal activation on both fast (within 100 ms) and slow (100 ms to seconds) timescales. Approach: We extensively modeled the variations in the DCS signal that are expected to arise from neuronal activation using Monte Carlo simulations, including the impacts of neuronal cell motion, vessel wall dilation, and blood flow changes. Results: We found that neuronal cell motion induces a DCS signal variation of ∼ 10 - 5 . We also estimated the contrast and number of channels required to detect hemodynamic signals at different time delays. Conclusions: From this extensive analysis, we do not expect to detect neuronal cell motion using DCS in the near future based on current technology trends. However, multi-channel DCS will be able to detect hemodynamic response with sub-second latency, which is interesting for brain-computer interfaces.

Global motion detection and censoring in high-density diffuse optical tomography

Motion-induced artifacts can significantly corrupt optical neuroimaging, as in most neuroimaging modalities. For high-density diffuse optical tomography (HD-DOT) with hundreds to thousands of source-detector pair measurements, motion detection methods are underdeveloped relative to both functional magnetic resonance imaging (fMRI) and standard functional near-infrared spectroscopy (fNIRS). This limitation restricts the application of HD-DOT in many challenging imaging situations and subject populations (e.g., bedside monitoring and children). Here, we evaluated a new motion detection method for multi-channel optical imaging systems that leverages spatial patterns across measurement channels. Specifically, we introduced a global variance of temporal derivatives (GVTD) metric as a motion detection index. We showed that GVTD strongly correlates with external measures of motion and has high sensitivity and specificity to instructed motion-with an area under the receiver operator characteristic curve of 0.88, calculated based on five different types of instructed motion. Additionally, we showed that applying GVTD-based motion censoring on both hearing words task and resting state HD-DOT data with natural head motion results in an improved spatial similarity to fMRI mapping. We then compared the GVTD similarity scores with several commonly used motion correction methods described in the fNIRS literature, including correlation-based signal improvement (CBSI), temporal derivative distribution repair (TDDR), wavelet filtering, and targeted principal component analysis (tPCA). We find that GVTD motion censoring on HD-DOT data outperforms other methods and results in spatial maps more similar to those of matched fMRI data.

Frequency-Domain Techniques for Cerebral and Functional Near-Infrared Spectroscopy

This article reviews the basic principles of frequency-domain near-infrared spectroscopy (FD-NIRS), which relies on intensity-modulated light sources and phase-sensitive optical detection, and its non-invasive applications to the brain. The simpler instrumentation and more straightforward data analysis of continuous-wave NIRS (CW-NIRS) accounts for the fact that almost all the current commercial instruments for cerebral NIRS have embraced the CW technique. However, FD-NIRS provides data with richer information content, which complements or exceeds the capabilities of CW-NIRS. One example is the ability of FD-NIRS to measure the absolute optical properties (absorption and reduced scattering coefficients) of tissue, and thus the absolute concentrations of oxyhemoglobin and deoxyhemoglobin in brain tissue. This article reviews the measured values of such optical properties and hemoglobin concentrations reported in the literature for animal models and for the human brain in newborns, infants, children, and adults. We also review the application of FD-NIRS to functional brain studies that focused on slower hemodynamic responses to brain activity (time scale of seconds) and faster optical signals that have been linked to neuronal activation (time scale of 100 ms). Another example of the power of FD-NIRS data is related to the different regions of sensitivity featured by intensity and phase data. We report recent developments that take advantage of this feature to maximize the sensitivity of non-invasive optical signals to brain tissue relative to more superficial extracerebral tissue (scalp, skull, etc.). We contend that this latter capability is a highly appealing quality of FD-NIRS, which complements absolute optical measurements and may result in significant advances in the field of non-invasive optical sensing of the brain.

Human brain function and brain-computer interfaces

Human brain function research has evolved dramatically in the last decades. In this chapter the role of modern methods of recording brain activity in understanding human brain function is explained. Current knowledge of brain function relevant to brain-computer interface (BCI) research is detailed, with an emphasis on the motor system which provides an exceptional level of detail to decoding of intended or attempted movements in paralyzed beneficiaries of BCI technology and translation to computer-mediated actions. BCI technologies that stand to benefit the most of the detailed organization of the human cortex are, and for the foreseeable future are likely to be, reliant on intracranial electrodes. These evolving technologies are expected to enable severely paralyzed people to regain the faculty of movement and speech in the coming decades.

Non-sedated functional imaging based on deep synchronization of PROPELLER MRI and NIRS

BACKGROUND AND OBJECTIVE

Periodically rotated overlapping parallel lines with enhanced reconstruction-echo planar imaging (PROPELLER-EPI) is a promising technique for non-sedated functional imaging due to its unique advantage of motion correction. However, its multiple-blades sampling blood-oxygen-level dependent (BOLD) signal leads to low sampling rate and aliasing of higher frequency physiological signal components such as the cardiac pulsation.

METHODS

In this study, we use near infrared spectroscopy (NIRS) synchronized with pulse sequences of PROPELLER-EPI, utilizing the fact that the optical sensing speed is inherently high. NIRS measures changes of oxyhemoglobin and deoxyhemoglobin to identify the transient states of on-BOLD and off-BOLD, and then labels each blade by temporal co-registration. The labeled blades from multiple epochs of a functional experiment are then used for the k-space data combination and subsequent image reconstruction. An eigenfunction model is proposed for temporal co-registration and to quantify the temporal resolution of the hemodynamic response.

RESULT

The experiment of NIRS labeled PROPELLER-EPI was carried out with the optical sampling rate of 10 Hz and the magnetic pulses repetition time of 1000 ms, and the temporal resolution is 20 times better than that of the state-of-the-art sliding-window PROPELLER-EPI. We compared the functional imaging results against the conventional magnetic resonance echo planar imaging-measured activity and achieved an accuracy of 0.9.

CONCLUSIONS

Using the synchronization of NIRS, the proposed imaging scheme provides an effective way to implement PROPELLER-EPI, which features motion free, high SNR, and enhanced spatial-temporal resolution.

Back-propagation neural network-based reconstruction algorithm for diffuse optical tomography

Diffuse optical tomography (DOT) is a promising noninvasive imaging modality and is capable of providing functional characteristics of biological tissue by quantifying optical parameters. The DOT image reconstruction is ill-posed and ill-conditioned, due to the highly diffusive nature of light propagation in biological tissues and limited boundary measurements. The widely used regularization technique for DOT image reconstruction is Tikhonov regularization, which tends to yield oversmoothed and low-quality images containing severe artifacts. It is necessary to accurately choose a regularization parameter for Tikhonov regularization. To overcome these limitations, we develop a noniterative reconstruction method, whereby optical properties are recovered based on a back-propagation neural network (BPNN). We train the parameters of BPNN before DOT image reconstruction based on a set of training data. DOT image reconstruction is achieved by implementing a single evaluation of the trained network. To demonstrate the performance of the proposed algorithm, we compare with the conventional Tikhonov regularization-based reconstruction method. The experimental results demonstrate that image quality and quantitative accuracy of reconstructed optical properties are significantly improved with the proposed algorithm.

Weighting function effects in a direct regularization method for image-guided near-infrared spectral tomography of breast cancer

Structural image-guided near-infrared spectral tomography (NIRST) has been developed as a way to use diffuse NIR spectroscopy within the context of image-guided quantification of tissue spectral features. A direct regularization imaging (DRI) method for NIRST has the value of not requiring any image segmentation. Here, we present a comprehensive investigational study to analyze the impact of the weighting function implied when weighting the recovery of optical coefficients in DRI based NIRST. This was done using simulations, phantom and clinical patient exam data. Simulations where the true object is known indicate that changes to this weighting function can vary the contrast by 10%, the contrast to noise ratio by 20% and the full width half maximum (FWHM) by 30%. The results from phantoms and human images show that a linear inverse distance weighting function appears optimal, and that incorporation of this function can generally improve the recovered total hemoglobin contrast of the tumor to the normal surrounding tissue by more than 15% in human cases.

Quantifying fast optical signal and event-related potential relationships during a visual oddball task

Event-related potentials (ERPs) have previously been used to confirm the existence of the fast optical signal (FOS) but validation methods have mainly been limited to exploring the temporal correspondence of FOS peaks to those of ERPs. The purpose of this study was to systematically quantify the relationship between FOS and ERP responses to a visual oddball task in both time and frequency domains. Near-infrared spectroscopy (NIRS) and electroencephalography (EEG) sensors were co-located over the prefrontal cortex while participants performed a visual oddball task. Fifteen participants completed 2 data collection sessions each, where they were instructed to keep a mental count of oddball images. The oddball condition produced a positive ERP at 200 ms followed by a negativity 300-500 ms after image onset in the frontal electrodes. In contrast to previous FOS studies, a FOS response was identified only in DC intensity signals and not in phase delay signals. A decrease in DC intensity was found 150-250 ms after oddball image onset with a 400-trial average in 10 of 15 participants. The latency of the positive 200 ms ERP and the FOS DC intensity decrease were significantly correlated for only 6 (out of 15) participants due to the low signal-to-noise ratio of the FOS response. Coherence values between the FOS and ERP oddball responses were found to be significant in the 3-5 Hz frequency band for 10 participants. A significant Granger causal influence of the ERP on the FOS oddball response was uncovered in the 2-6 Hz frequency band for 7 participants. Collectively, our findings suggest that, for a majority of participants, the ERP and the DC intensity signal of the FOS are spectrally coherent, specifically in narrow frequency bands previously associated with event-related oscillations in the prefrontal cortex. However, these electro-optical relationships were only found in a subset of participants. Further research on enhancing the quality of the event-related FOS signal is required before it can be practically exploited in applications such as brain-computer interfacing.

Cortical light scattering during interictal epileptic spikes in frontal lobe epilepsy in children: A fast optical signal and electroencephalographic study

OBJECTIVE

Interictal epileptic spikes (IES) represent a signature of the transient synchronous and excessive discharge of a large ensemble of cortical heterogeneous neurons. Epilepsy cannot be reduced to a hypersynchronous activation of neurons whose functioning is impaired, resulting on electroencephalogram (EEG) in epileptic seizures or IES. The complex pathophysiological mechanisms require a global approach to the interactions between neural synaptic and nonsynaptic, vascular, and metabolic systems. In the present study, we focused on the interaction between synaptic and nonsynaptic mechanisms through the simultaneous noninvasive multimodal multiscale recording of high-density EEG (HD-EEG; synaptic) and fast optical signal (FOS; nonsynaptic), which evaluate rapid changes in light scattering related to changes in membrane configuration occurring during neuronal activation of IES.

METHODS

To evaluate changes in light scattering occurring around IES, three children with frontal IES were simultaneously recorded with HD-EEG and FOS. To evaluate change in synchronization, time-frequency representation analysis of the HD-EEG was performed simultaneously around the IES. To independently evaluate our multimodal method, a control experiment with somatosensory stimuli was designed and applied to five healthy volunteers.

RESULTS

Alternating increase-decrease-increase in optical signals occurred 200 ms before to 180 ms after the IES peak. These changes started before any changes in EEG signal. In addition, time-frequency domain EEG analysis revealed alternating decrease-increase-decrease in the EEG spectral power concomitantly with changes in the optical signal during IES. These results suggest a relationship between (de)synchronization and neuronal volume changes in frontal lobe epilepsy during IES.

SIGNIFICANCE

These changes in the neuronal environment around IES in frontal lobe epilepsy observed in children, as they have been in rats, raise new questions about the synaptic/nonsynaptic mechanisms that propel the neurons to hypersynchronization, as occurs during IES. We further demonstrate that this noninvasive multiscale multimodal approach is suitable for studying the pathophysiology of the IES in patients.

Process-specific analysis in episodic memory retrieval using fast optical signals and hemodynamic signals in the right prefrontal cortex

OBJECTIVE

Memory is formed by the interaction of various brain functions at the item and task level. Revealing individual and combined effects of item- and task-related processes on retrieving episodic memory is an unsolved problem because of limitations in existing neuroimaging techniques. To investigate these issues, we analyze fast and slow optical signals measured from a custom-built continuous wave functional near-infrared spectroscopy (CW-fNIRS) system.

APPROACH

In our work, we visually encode the words to the subjects and let them recall the words after a short rest. The hemodynamic responses evoked by the episodic memory are compared with those evoked by the semantic memory in retrieval blocks. In the fast optical signal, we compare the effects of old and new items (previously seen and not seen) to investigate the item-related process in episodic memory. The Kalman filter is simultaneously applied to slow and fast optical signals in different time windows.

MAIN RESULTS

A significant task-related HbR decrease was observed in the episodic memory retrieval blocks. Mean amplitude and peak latency of a fast optical signal are dependent upon item types and reaction time, respectively. Moreover, task-related hemodynamic and item-related fast optical responses are correlated in the right prefrontal cortex.

SIGNIFICANCE

We demonstrate that episodic memory is retrieved from the right frontal area by a functional connectivity between the maintained mental state through retrieval and item-related transient activity. To the best of our knowledge, this demonstration of functional NIRS research is the first to examine the relationship between item- and task-related memory processes in the prefrontal area using single modality.

Signal processing of functional NIRS data acquired during overt speaking

Functional near-infrared spectroscopy (fNIRS) offers an advantage over traditional functional imaging methods [such as functional magnetic resonance imaging (fMRI)] by allowing participants to move and speak relatively freely. However, neuroimaging while actively speaking has proven to be particularly challenging due to the systemic artifacts that tend to be located in the critical brain areas. To overcome these limitations and enhance the utility of fNIRS, we describe methods for investigating cortical activity during spoken language tasks through refinement of deoxyhemoglobin (deoxyHb) signals with principal component analysis (PCA) spatial filtering to remove global components. We studied overt picture naming and compared oxyhemoglobin (oxyHb) and deoxyHb signals with and without global component removal using general linear model approaches. Activity in Broca's region and supplementary motor cortex was observed only when the filter was applied to the deoxyHb signal and was shown to be spatially comparable to fMRI data acquired using a similar task and to meta-analysis data. oxyHb signals did not yield expected activity in Broca's region with or without global component removal. This study demonstrates the utility of a PCA spatial filter on the deoxyHb signal in revealing neural activity related to a spoken language task and extends applications of fNIRS to natural and ecologically valid conditions.

From brain to blood vessels and back: a noninvasive optical imaging approach

The seminal work of Grinvald et al. has paved the way for the use of intrinsic optical signals measured with reflection methods for the analysis of brain function. Although this work has focused on the absorption signal associated with deoxygenation, due to its detailed mapping ability and good signal-to-noise ratio, Grinvald's group has also described other intrinsic signals related to increased blood flow, scattering effects directly related to neural activation, and pulsation effects related to arterial function. These intrinsic optical signals can also be measured using noninvasive diffuse optical topographic and tomographic imaging (DOT) methods that can be applied to humans. Here we compare the reflection and DOT methods and the evidence for each type of intrinsic signal in these two domains, with particular attention to work that has been conducted in our laboratory. This work reveals the refined two-way relationship that exists between vascular and neural phenomena in the brain: arterial health is related to normal brain structure and function, both across individuals and across brain regions within an individual, and neural function influences blood flow to specific cortical regions. DOT methods can provide quantitative tools for investigating these relationships in normal human subjects.

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.

Hemodynamic Response to Interictal Epileptiform Discharges Addressed by Personalized EEG-fNIRS Recordings

Objective: We aimed at studying the hemodynamic response (HR) to Interictal Epileptic Discharges (IEDs) using patient-specific and prolonged simultaneous ElectroEncephaloGraphy (EEG) and functional Near InfraRed Spectroscopy (fNIRS) recordings.

Methods: The epileptic generator was localized using Magnetoencephalography source imaging. fNIRS montage was tailored for each patient, using an algorithm to optimize the sensitivity to the epileptic generator. Optodes were glued using collodion to achieve prolonged acquisition with high quality signal. fNIRS data analysis was handled with no a priori constraint on HR time course, averaging fNIRS signals to similar IEDs. Cluster-permutation analysis was performed on 3D reconstructed fNIRS data to identify significant spatio-temporal HR clusters. Standard (GLM with fixed HRF) and cluster-permutation EEG-fMRI analyses were performed for comparison purposes.

Results: fNIRS detected HR to IEDs for 8/9 patients. It mainly consisted oxy-hemoglobin increases (seven patients), followed by oxy-hemoglobin decreases (six patients). HR was lateralized in six patients and lasted from 8.5 to 30 s. Standard EEG-fMRI analysis detected an HR in 4/9 patients (4/9 without enough IEDs, 1/9 unreliable result). The cluster-permutation EEG-fMRI analysis restricted to the region investigated by fNIRS showed additional strong and non-canonical BOLD responses starting earlier than the IEDs and lasting up to 30 s.

Conclusions: (i) EEG-fNIRS is suitable to detect the HR to IEDs and can outperform EEG-fMRI because of prolonged recordings and greater chance to detect IEDs; (ii) cluster-permutation analysis unveils additional HR features underestimated when imposing a canonical HR function (iii) the HR is often bilateral and lasts up to 30 s.

Combining energy and Laplacian regularization to accurately retrieve the depth of brain activity of diffuse optical tomographic data

Diffuse optical tomography (DOT) provides data about brain function using surface recordings. Despite recent advancements, an unbiased method for estimating the depth of absorption changes and for providing an accurate three-dimensional (3-D) reconstruction remains elusive. DOT involves solving an ill-posed inverse problem, requiring additional criteria for finding unique solutions. The most commonly used criterion is energy minimization (energy constraint). However, as measurements are taken from only one side of the medium (the scalp) and sensitivity is greater at shallow depths, the energy constraint leads to solutions that tend to be small and superficial. To correct for this bias, we combine the energy constraint with another criterion, minimization of spatial derivatives (Laplacian constraint, also used in low resolution electromagnetic tomography, LORETA). Used in isolation, the Laplacian constraint leads to solutions that tend to be large and deep. Using simulated, phantom, and actual brain activation data, we show that combining these two criteria results in accurate (error <2 mm) absorption depth estimates, while maintaining a two-point spatial resolution of <24 mm up to a depth of 30 mm. This indicates that accurate 3-D reconstruction of brain activity up to 30 mm from the scalp can be obtained with DOT.

Optical stimulation of primary motor cortex with 980nm infrared neural stimulation

To explore the penetration depth with short-wavelength infrared light, 980 nm pulse infrared light was used to stimulate the primary motor cortex of rat. The heating model was created to simulate the temperature distribution for 1875 nm and 980 nm infrared neural stimulation. Post-stimulus time histogram was used to observe the neural response induced by Infrared neural stimulation on primary motor cortex. The model predicted the penetration depth of 980 nm was deep into 1.2 mm. Cortical neural located between 500 μm to 1000 μm were successfully activated by 980 nm INS. The preliminary results suggested that, 980 nm pulse INS could serve as a candidate for deep tissue stimulation.

fNIRS in the developmental sciences

With the introduction of functional near-infrared spectroscopy (fNIRS) into the experimental setting, developmental scientists have, for the first time, the capacity to investigate the functional activation of the infant brain in awake, engaged participants. The advantages of fNIRS clearly outweigh the limitations, and a description of how this technology is implemented in infant populations is provided. Most fNIRS research falls into one of three content domains: object processing, processing of biologically and socially relevant information, and language development. Within these domains, there are ongoing debates about the origins and development of human knowledge, making early neuroimaging particularly advantageous. The use of fNIRS has allowed investigators to begin to identify the localization of early object, social, and linguistic knowledge in the immature brain and the ways in which this changes with time and experience. In addition, there is a small but growing body of research that provides insight into the neural mechanisms that support and facilitate learning during the first year of life. At the same time, as with any emerging field, there are limitations to the conclusions that can be drawn on the basis of current findings. We offer suggestions as to how to optimize the use of this technology to answer questions of theoretical and practical importance to developmental scientists.

Determination of detection depth of optical probe in pedicle screw measurement device

BACKGROUND

There is a high probability of accidental perforation of the vertebral pedicle wall in pedicle screw insertion surgery. A pedicle screw (PS) measurement device with an optical probe has been reported to send out a warning signal before the PS tip breaking the vertebral pedicle wall.

METHODS

In this study, we explored the detection depth of optical probe in this measurement device, which was closely related to the effective alarm distance. In the boundary, the vertebrae tissues could be treated as 2-layer models including spongy bones and compact bones. The Monte Carlo simulation and phantom models were performed to analyse and define the detection depth. Then the porcine vertebrae models were performed to obtain optical spectrum and reduced scattering coefficient, based on which the detection depths were deduced. Moreover, a comparison was made to explore the most significant pattern factor from the experiment results.

RESULTS

According to the pattern factor, an alarm threshold was successfully deduced to define the alarm distance during pedicle screw monitoring.

CONCLUSIONS

Thus, the proposed alarm standard based on detection depth provides a potential for guiding pedicle screw in surgery.

Dynamics of alpha control: preparatory suppression of posterior alpha oscillations by frontal modulators revealed with combined EEG and event-related optical signal

We investigated the dynamics of brain processes facilitating conscious experience of external stimuli. Previously, we proposed that alpha (8-12 Hz) oscillations, which fluctuate with both sustained and directed attention, represent a pulsed inhibition of ongoing sensory brain activity. Here we tested the prediction that inhibitory alpha oscillations in visual cortex are modulated by top-down signals from frontoparietal attention networks. We measured modulations in phase-coherent alpha oscillations from superficial frontal, parietal, and occipital cortices using the event-related optical signal (EROS), a measure of neuronal activity affording high spatiotemporal resolution, along with concurrently recorded EEG, while participants performed a visual target detection task. The pretarget alpha oscillations measured with EEG and EROS from posterior areas were larger for subsequently undetected targets, supporting alpha's inhibitory role. Using EROS, we localized brain correlates of these awareness-related alpha oscillations measured at the scalp to the cuneus and precuneus. Crucially, EROS alpha suppression correlated with posterior EEG alpha power across participants. Sorting the EROS data based on EEG alpha power quartiles to investigate alpha modulators revealed that suppression of posterior alpha was preceded by increased activity in regions of the dorsal attention network and decreased activity in regions of the cingulo-opercular network. Cross-correlations revealed the temporal dynamics of activity within these preparatory networks before posterior alpha modulation. The novel combination of EEG and EROS afforded localization of the sources and correlates of alpha oscillations and their temporal relationships, supporting our proposal that top-down control from attention networks modulates both posterior alpha and awareness of visual stimuli.

Neurophotonics: non-invasive optical techniques for monitoring brain functions

The aim of this review is to present the state of the art of neurophotonics, a recently founded discipline lying at the interface between optics and neuroscience. While neurophotonics also includes invasive techniques for animal studies, in this review we focus only on the non-invasive methods that use near infrared light to probe functional activity in the brain, namely the fast optical signal, diffuse correlation spectroscopy, and functional near infrared spectroscopy methods. We also present an overview of the physical principles of light propagation in biological tissues, and of the main physiological sources of signal. Finally, we discuss the open issues in models, instrumentation, data analysis and clinical approaches.

Subsurface fluorescence molecular tomography with prior information

Subsurface fluorescence molecular tomography (FMT) is an emerging technique determining fluorescence distribution by tomographic means in reflectance geometry. However, due to the highly diffusive nature of the photon propagation in biological tissues and the influence of nearer source-detector separations, stand-alone subsurface FMT could not accurately reflect the fluorophore distributions. To overcome this drawback, we propose a method to improve the performance of fluorescence imaging by coupling x-ray computed tomography (XCT) and subsurface FMT modalities. A Laplacian-type regularization matrix generated with tissue prior information obtained from XCT images is used to guide the reconstruction of fluorophore distribution. Reconstruction results of both simulation and phantom studies showed that significant improvements in localization and demarcation of fluorescent targets can be obtained with the proposed method compared to the reconstruction method without structural prior information.

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.

Intraoperative cortico–subcortical stimulations in surgery of low-grade gliomas

In order to increase the impact of surgery on the natural history of low-grade glioma, resection should be of maximum importance. Nevertheless, since low-grade gliomas are frequently located in eloquent structures, function needs to be preserved. Therefore, studying the functional organization of the brain is mandatory for each patient due to the inter-individual anatomofunctional variability, increased in tumors due to cerebral plasticity. This strategy enables performance of a resection according to functional boundaries. However, preoperative neurofunctional imaging only allows the study of the gray matter. Consequently, since low-grade glioma invades cortical and subcortical structures and shows an infiltrative progression along the fibers, the goal of this review is to focus on the techniques able to map both cortical and subcortical regions. In addition to diffusion tensor imaging, which gives only anatomical information and still needs to be validated, intraoperative direct cortico-subcortical electrostimulation is the sole current method allowing a reliable study of the individual anatomofunctional connectivity, concerning sensorimotor, language and other cognitive functions. Its actual contribution is detailed, both in clinical issues, especially the improvement of the benefit/risk ratio of low-grade glioma resection, and in fundamental applications--namely, a new door to the connectionism and cerebral plasticity.

Detection of optical neuronal signals in the visual cortex using continuous wave near-infrared spectroscopy

Near-infrared spectroscopy (NIRS) measures slow hemodynamic signals noninvasively to indirectly infer the neuronal activity in the brain. However, it remains a controversy on whether this optical measurement technique can detect the optical neuronal signal, which reflects the optical changes directly associated with neuronal activity, within the visual cortex of human and non-human primates. By carefully reviewing the important factors in the detection of optical neuronal signals, we aim to investigate the feasibility of performing NIRS measurements of optical neuronal signals within the visual cortex in humans. To ensure a strong optical neuronal response, a full-field circular black and white reversing checkerboard stimulus was presented, and the reversal frequency was carefully chosen. We used a homemade continuous wave (CW) NIRS system with high detection sensitivity (of the order of 0.1 pW) to record a large area of the visual cortex (approximately 6 × 14 cm(2)). EEG was simultaneously acquired with the optical signal. Based on the mathematical morphology, we adapted the filter proposed by Gratton et al. to remove the influence of arterial pulsation and facilitate the detection and elimination of unknown artifacts from the data. We obtained reliable optical neuronal signals in 77% of the participants (10 out of 13). The amplitudes (latencies) of the obtained optical neuronal signals corresponding to the 785 and 850 nm wavelengths were 0.017 ± 0.003% (94.7 ± 8.4 ms) and 0.025 ± 0.006% (99.0 ± 7.7 ms), respectively. There were no significant differences between the latencies of the N75 component of the visual evoked potential (VEP) and optical neuronal signals at either wavelength. This is the first study to report optical neuronal signals within the visual cortex in the intact human brain using a CW NIRS system. These results indicate the feasibility of measuring noninvasive optical neuronal signals using a CW NIRS system with high detection sensitivity.

Fast optical signals in the sensorimotor cortex: General Linear Convolution Model applied to multiple source-detector distance-based data

In this study, we applied the General Linear Convolution Model to detect fast optical signals (FOS) in the somatosensory cortex, and to study their dependence on the source-detector separation distance (2.0 to 3.5 cm) and irradiated light wavelength (690 and 830 nm). We modeled the impulse response function as a rectangular function that lasted 30 ms, with variable time delay with respect to the stimulus onset. The model was tested in a cohort of 20 healthy volunteers who underwent supra-motor threshold electrical stimulation of the median nerve. The impulse response function quantified the time delay for the maximal response at 70 ms to 110 ms after stimulus onset, in agreement with classical somatosensory-evoked potentials in the literature, previous optical imaging studies based on a grand-average approach, and grand-average based processing. Phase signals at longer wavelength were used to identify FOS for all the source-detector separation distances, but the shortest one. Intensity signals only detected FOS at the greatest distance; i.e., for the largest channel depth. There was no activation for the shorter wavelength light. Correlational analysis between the phase and intensity of FOS further confirmed diffusive rather than optical absorption changes associated with neuronal activity in the activated cortical volume. Our study demonstrates the reliability of our method based on the General Linear Convolution Model for the detection of fast cortical activation through FOS.

Pulsed 808-nm infrared laser stimulation of the auditory nerve in guinea pig cochlea

Pulsed near-infrared radiation has been proposed as an alternative stimulus for auditory nerve stimulation and could be potentially used in the design of cochlear implant. Although the infrared with high absorption coefficient of water (i.e., wavelength ranged from 1.8 to 2.2 μm) has been widely investigated, the lymph in the cochlea absorbs most of the infrared energies, and only a small part can arrive at the target auditory nerves. The present study is aimed to test whether the short-wavelength near-infrared irradiation with lower absorption coefficients can penetrate the lymph fluid to stimulate the auditory nerves. An 808-nm near-infrared laser was chosen to stimulate the auditory nerve in the guinea pig cochlea. The infrared pulse was delivered by an optical fiber that was surgically inserted near the round window membrane and oriented toward the spiral ganglion cells in the basal turn of the cochlea. The 2-Hz infrared pulses were used to stimulate the cochlea before and after the deafness with different pulse durations (100-1,000 μs). Optically evoked compound action potentials (oCAPs) were recorded during the infrared radiation. We successfully recorded oCAPs from both normal hearing animals and deafened animals. The oCAP amplitude increased with the infrared radiation energy. The preliminary experiment suggests that the near-infrared with lower absorption coefficients can effectively pass through the lymph filled in the cochlea and stimulate the auditory nerve. Further studies will optimize the deafness animal model and determine the optimal stimulation parameters.

Brain structure and spatial sensitivity profile assessing by near-infrared spectroscopy modeling based on 3D MRI data

The goal of this study is to prove that the light propagation in the head by used the 3-D optical model from in vivo MRI data set can also provide significant characteristics on the spatial sensitivity of cerebral cortex folding geometry based on Monte Carlo simulation. Thus, we proposed a MRI based approach for 3-D brain modeling of near-infrared spectroscopy (NIRS). In the results, the spatial sensitivity profile of the cerebral cortex folding geometry and the arrangement of source-detector separation have being necessarily considered for applications of functional NIRS. The optimal choice of source-detector separation is suggested within 3-3.5 cm by the received intensity with different source-detector separations and the ratio of received light from the gray and white matter layer is greater than 50%. Additionally, this study has demonstrated the capability of NIRS in not only assessing the functional but also detecting the structural change of the brain by taking advantage of the low scattering and absorption coefficients observed in CSF of sagittal view.

Quantitative evaluation of deep and shallow tissue layers' contribution to fNIRS signal using multi-distance optodes and independent component analysis

To quantify the effect of absorption changes in the deep tissue (cerebral) and shallow tissue (scalp, skin) layers on functional near-infrared spectroscopy (fNIRS) signals, a method using multi-distance (MD) optodes and independent component analysis (ICA), referred to as the MD-ICA method, is proposed. In previous studies, when the signal from the shallow tissue layer (shallow signal) needs to be eliminated, it was often assumed that the shallow signal had no correlation with the signal from the deep tissue layer (deep signal). In this study, no relationship between the waveforms of deep and shallow signals is assumed, and instead, it is assumed that both signals are linear combinations of multiple signal sources, which allows the inclusion of a "shared component" (such as systemic signals) that is contained in both layers. The method also assumes that the partial optical path length of the shallow layer does not change, whereas that of the deep layer linearly increases along with the increase of the source-detector (S-D) distance. Deep- and shallow-layer contribution ratios of each independent component (IC) are calculated using the dependence of the weight of each IC on the S-D distance. Reconstruction of deep- and shallow-layer signals are performed by the sum of ICs weighted by the deep and shallow contribution ratio. Experimental validation of the principle of this technique was conducted using a dynamic phantom with two absorbing layers. Results showed that our method is effective for evaluating deep-layer contributions even if there are high correlations between deep and shallow signals. Next, we applied the method to fNIRS signals obtained on a human head with 5-, 15-, and 30-mm S-D distances during a verbal fluency task, a verbal working memory task (prefrontal area), a finger tapping task (motor area), and a tetrametric visual checker-board task (occipital area) and then estimated the deep-layer contribution ratio. To evaluate the signal separation performance of our method, we used the correlation coefficients of a laser-Doppler flowmetry (LDF) signal and a nearest 5-mm S-D distance channel signal with the shallow signal. We demonstrated that the shallow signals have a higher temporal correlation with the LDF signals and with the 5-mm S-D distance channel than the deep signals. These results show the MD-ICA method can discriminate between deep and shallow signals.

Label-free imaging of membrane potential using membrane electromotility

Electrical activity may cause observable changes in a cell's structure in the absence of exogenous reporter molecules. In this work, we report a low-coherence interferometric microscopy technique that can detect an optical signal correlated with the membrane potential changes in individual mammalian cells without exogenous labels. By measuring milliradian-scale phase shifts in the transmitted light, we can detect changes in the cells' membrane potential. We find that the observed optical signals are due to membrane electromotility, which causes the cells to deform in response to the membrane potential changes. We demonstrate wide-field imaging of the propagation of electrical stimuli in gap-junction-coupled cell networks. Membrane electromotility-induced cell deformation may be useful as a reporter of electrical activity.

A near-infrared spectroscopy computational model for cerebral hemodynamics

Near infrared spectroscopy (NIRS) is a technique used to detect and measure changes in the concentrations of oxygenated hemoglobin, deoxygenated hemoglobin, and water in tissues based on the differential absorption, scattering, and refraction of the near infrared light. In this imaging technique, the optical properties of tissues are reconstructed from the measurements obtained from the sensors located on the boundary. A computational method for the rapid noninvasive detection ∕ quantification of cerebral hemorrhage is described using the above procedure. CFD Research Corporation's finite volume computational biology code was used to numerically mimic the NIRS procedure by (i) noninvasively 'numerically penetrating' the brain tissues and (ii) reconstructing the optical properties the presence of water, oxygenated, and deoxygenated blood. These numerical noninvasive measurements are then used to predict the extent and severity of the brain hemorrhage. The paper also discusses ideas to obtain the location and the severity of a localized injury. Two-dimensional and three-dimensional simulations are performed as a proof of concept for the numerical formulation being feasible for the above mentioned detection/quantification. The results demonstrate that this numerical NIRS formulation can be used as a noninvasive technique for both qualitative and quantitative evaluation of cerebral hemodynamics.

Fast optical signal in visual cortex: Improving detection by General Linear Convolution Model

In this study we applied the General Linear Convolution Model to fast optical signals (FOS). We modeled the Impulse Response Function (IRF) as a rectangular function lasting 30ms, with variable time delay with respect to the stimulus onset. Simulated data confirmed the feasibility of this approach and its capability of detecting simulated activations in case of very unfavorable Signal to Noise Ratio (SNR), providing better results than the grand average method. The model was tested in a cohort of 10 healthy volunteers who underwent to hemi-field visual stimulation. Experimental data quantified the IRF time delay at 80-100ms after the stimulus onset, in agreement with classical visual evoked potential literature and previous optical imaging studies based on grand average approach and a larger number of trails. FOS confirmed the expected contralateral activation in the occipital region. Correlational analysis between hemodynamic intensity signal, phase and intensity FOS supports diffusive rather than optical absorption changes associated with neuronal activity in the activated cortical volume. Our study provides a feasible method for detecting fast cortical activations by means of FOS.

A fast neuronal signal-sensitive continuous-wave near-infrared imaging system

We have developed a continuous-wave near-infrared imaging system to measure fast neuronal signals. We used a simultaneous sampling method with a separate high-speed analog-to-digital converter for each input channel, which provides a much larger point sample in a digital lock-in algorithm, higher temporal resolution, and lower crosstalk among detected channels. Without any analog filter, digital lock-in detection with a large point sample suppresses noise excellently, making the system less complex and offering better flexibility. In addition, using a custom-made collimator, more photons can reach the brain tissue due to the smaller divergence angle. Performance analysis shows high detection sensitivity (on the order of 0.1 pW) and high temporal resolution (~50 Hz, 48 channels). Simulation experiments show that intensity changes on the order of 0.01% can be resolved by our instrument when averaging over approximately 500 stimuli. In vivo experiments over the motor cortex show that our instrument can detect fast neuronal signals in the human brain.

Multimodal recording of brain activity in term newborns during photic stimulation by near-infrared spectroscopy and electroencephalography

In this study 14 healthy term newborns (postnatal mean age 2.1 days) underwent photic stimulation during sleep on two different days. Near-infrared spectroscopy (NIRS) and electroencephalography (EEG) was acquired simultaneously. The aims of the study were: to determine (i) the sensitivity and (ii) the repeatability of NIRS to detect the hemodynamic response, (iii) the sensitivity and (iv) the repeatability of EEG to detect a visual evoked potential (VEP), (v) to analyze optical data for the optical neuronal signal, and (vi) to test whether inadequate stimulation could be reason for absent hemodynamic responses. The results of the study were as follows. (i) Sensitivity of NIRS was 61.5% to detect hemodynamic responses; (ii) their reproducibility was 41.7%. A VEP was detected (iii) in 96.3% of all subjects with (iv) a reproducibility of 92.3%. (v) In two measurements data met the criteria for an optical neuronal signal. The noise level was 9.6·10-5% change in optical density. (vi) Insufficient stimulation was excluded as reason for absent hemodynamic responses. We conclude that NIRS is an promising tool to study cognitive activation and development of the brain. For clinical application, however, the sensitivity and reproducibility on an individual level needs to be improved.