Activation detection in diffuse optical imaging by means of the general linear model

Ipsilateral stimulation shows somatotopy of thumb and shoulder auricular points on the left primary somatosensory cortex using high-density fNIRS

Significance

Auriculotherapy is a technique based on stimulation applied to specific ear points. Its mechanism of active and clinical efficacy remain to be established. This study aims to assess the role that primary somatosensory cortex may play to validate auriculotherapy mechanisms.

Aim

This study examined whether tactile stimulation at specific auricular points is correlated with distinct cortical activation in the primary somatosensory cortex.

Approach

Seventeen healthy adults participated in the study. Tactile stimuli were delivered to the thumb, shoulder, and skin master points on the ear using von Frey filaments. Functional near-infrared spectroscopy was used to measure and spatially map cortical responses.

Results

This study revealed distinct hemodynamic activity patterns in response to ear point stimulation, consistent with the classic homunculus model of somatotopic organization. Ipsilateral stimulation showed specific cortical activations for the thumb and shoulder points, while contralateral stimulation showed less significant activity. Functional near-infrared spectroscopy effectively captured localized cortical responses to ear tactile stimuli, supporting the somatotopic mapping hypothesis.

Conclusion

These findings enhance the understanding of sensory processing with auricular stimulation and supports the concepts of auricular cartography that underpins some schools of auriculotherapy practice. Future research should explore bilateral cortical mapping and the integration of other neuroimaging techniques.

Investigating cortical activity during cybersickness by fNIRS

This study investigated brain responses during cybersickness in healthy adults using functional near-infrared spectroscopy (fNIRS). Thirty participants wore a head-mounted display and observed a virtual roller coaster scene that induced cybersickness. Cortical activation during the virtual roller coaster task was measured using fNIRS. Cybersickness symptoms were evaluated using a Simulator Sickness Questionnaire (SSQ) administered after the virtual rollercoaster. Pearson correlations were performed for cybersickness symptoms and the beta coefficients of hemodynamic responses. The group analysis of oxyhemoglobin (HbO) and total hemoglobin (HbT) levels revealed deactivation in the bilateral angular gyrus during cybersickness. In the Pearson correlation analyses, the HbO and HbT beta coefficients in the bilateral angular gyrus had a significant positive correlation with the total SSQ and disorientation. These results indicated that the angular gyrus was associated with cybersickness. These findings suggest that the hemodynamic response in the angular gyrus could be a biomarker for evaluating cybersickness symptoms.

Neuroimaging, wearable sensors, and blood-based biomarkers reveal hyperacute changes in the brain after sub-concussive impacts

Impacts in mixed martial arts (MMA) have been studied mainly in regard to the long-term effects of concussions. However, repetitive sub-concussive head impacts at the hyperacute phase (minutes after impact), are not understood. The head experiences rapid acceleration similar to a concussion, but without clinical symptoms. We utilize portable neuroimaging technology - transcranial Doppler (TCD) ultrasound and functional near infrared spectroscopy (fNIRS) - to estimate the extent of pre- and post-differences following contact and non-contact sparring sessions in nine MMA athletes. In addition, the extent of changes in neurofilament light (NfL) protein biomarker concentrations, and neurocognitive/balance parameters were determined following impacts. Athletes were instrumented with sensor-based mouth guards to record head kinematics. TCD and fNIRS results demonstrated significantly increased blood flow velocity (p = 0.01) as well as prefrontal (p = 0.01) and motor cortex (p = 0.04) oxygenation, only following the contact sparring sessions. This increase after contact was correlated with the cumulative angular acceleration experienced during impacts (p = 0.01). In addition, the NfL biomarker demonstrated positive correlations with angular acceleration (p = 0.03), and maximum principal and fiber strain (p = 0.01). On average athletes experienced 23.9 ± 2.9 g peak linear acceleration, 10.29 ± 1.1 rad/s peak angular velocity, and 1,502.3 ± 532.3 rad/s2 angular acceleration. Balance parameters were significantly increased following contact sparring for medial-lateral (ML) center of mass (COM) sway, and ML ankle angle (p = 0.01), illustrating worsened balance. These combined results reveal significant changes in brain hemodynamics and neurophysiological parameters that occur immediately after sub-concussive impacts and suggest that the physical impact to the head plays an important role in these changes.

Hyper-acute effects of sub-concussive soccer headers on brain function and hemodynamics

Introduction

Sub-concussive head impacts in soccer are drawing increasing research attention regarding their acute and long-term effects as players may experience thousands of headers in a single season. During these impacts, the head experiences rapid acceleration similar to what occurs during a concussion, but without the clinical implications. The physical mechanism and response to repetitive impacts are not completely understood. The objective of this work was to examine the immediate functional outcomes of sub-concussive level impacts from soccer heading in a natural, non-laboratory environment.

Methods

Twenty university level soccer athletes were instrumented with sensor-mounted bite bars to record impacts from 10 consecutive soccer headers. Pre- and post-header measurements were collected to determine hyper-acute changes, i.e., within minutes after exposure. This included measuring blood flow velocity using transcranial Doppler (TCD) ultrasound, oxyhemoglobin concentration using functional near infrared spectroscopy imaging (fNIRS), and upper extremity dual-task (UEF) neurocognitive testing.

Results

On average, the athletes experienced 30.7 ± 8.9 g peak linear acceleration and 7.2 ± 3.1 rad/s peak angular velocity, respectively. Results from fNIRS measurements showed an increase in the brain oxygenation for the left prefrontal cortex (PC) (p = 0.002), and the left motor cortex (MC) (p = 0.007) following the soccer headers. Additional analysis of the fNIRS time series demonstrates increased sample entropy of the signal after the headers in the right PC (p = 0.02), right MC (p = 0.004), and left MC (p = 0.04).

Discussion

These combined results reveal some variations in brain oxygenation immediately detected after repetitive headers. Significant changes in balance and neurocognitive function were not observed in this study, indicating a mild level of head impacts. This is the first study to observe hemodynamic changes immediately after sub-concussive impacts using non-invasive portable imaging technology. In combination with head kinematic measurements, this information can give new insights and a framework for immediate monitoring of sub-concussive impacts on the head.

Sensorimotor adaptation in spatial orientation task: a fNIRS study

In sensorimotor conflicts, the brain forms and updates a new sensorimotor relationship through sensorimotor integration. As humans adapt to new sensorimotor mapping, goal-directed movements become increasingly precise. Using functional near-infrared spectroscopy, we investigated the changes in cortical activity during sensorimotor adaptation in a spatial orientation task with sensorimotor conflict. Individuals performed a reversed spatial orientation training in which the visual feedback guiding hand movements was reversed. We measured cortical activity and spatial orientation performance, including the response time, completion number, error, and accuracy. The results revealed the continuous activation in the left SMG during sensorimotor adaptation and decreased activation in the right SAC, AG and SMG after sensorimotor adaptation. These findings indicated the contribution of the left SMG to sensorimotor adaptation and the improved efficiency of cortical activity after sensorimotor adaptation, respectively. Our studies suggest the neural mechanisms related to sensorimotor adaptation to a reversed spatial orientation task.

Decoupling brain activations of muscle-caused activations and mental intention-cause activations using the general linear model: A functional near-infrared spectroscopy study

Integrating a brain-computer interface into a lower-limb medical rehabilitation assistive device can enhance rehabilitation efficiency. The latest research in the field focuses on the decoding performance of different motions. However, the difference between muscle-caused primitive activation and mental intention-caused activation has not been fully investigated. Thus, our study tried to decouple these two kinds of cerebral activation using a general linear model (GLM). Nine healthy and right-handed subjects were recruited for a two-section experiment. They were asked to extend or flex their knees while seated in the first section of the experiment or standing in the second section. Functional near-infrared spectroscopy (fNIR) was adopted to monitor their hemodynamic changes. Two groups of paradigms (one for circle-wise analysis, the other for full-section analysis) were constructed from the experimental paradigm. Each group consisted of three (the first intention, the second intention, and the muscle activation). The constructed paradigms were fed to the Balloon model for six desired hemodynamic responses (dHRFs). The regressor of GLM consisted of three dHRFs and the corresponding motion artifacts and drifts. The simulated physiological noises were included in the structured background matrix. The results showed that all subjects had similar cerebral activation patterns for the intention to extend or flex knees. The activation during musclecaused activation was less intense than that caused by both intentions. This finding can help further research on more efficient motion intention detection and the possibility of multiple motions decoding.

The neurovascular couplings between electrophysiological and hemodynamic activities in anticipatory selective attention

Selective attention is thought to involve target enhancement and distractor inhibition processes. Here, we recorded simultaneous electroencephalographic (EEG) and functional near-infrared spectroscopy (fNIRS) data from human adults when they were pre-cued by the visual field of coming target, distractor, or both of them. From the EEG data, we found alpha power relatively decreased contralaterally to the to-be-attended target, as reflected by the positive-going alpha modulation index. Late alpha power relatively increased contralaterally to the to-be-suppressed distractor, as reflected by the negative-going alpha modulation index. From the fNIRS data, we found enhancements of hemodynamic activity over the contralateral hemisphere in response to both the target and the distractor anticipation but within nonoverlapping posterior brain regions. More importantly, we described the specific neurovascular modulation between alpha power and oxygenated hemoglobin signal, which showed a positive coupling effect during target anticipation and a negative coupling effect during distractor anticipation. Such flexible neurovascular couplings between EEG oscillation and hemodynamic activity seem to play an essential role in the final behavioral outcomes. These results provide unique neurovascular evidence for the dissociation of the mechanisms of target enhancement and distractor inhibition. Individual behavioral differences can be related to individual differences in neurovascular coupling.

Cognitive deficit in adults with ADHD lies in the cognitive state disorder rather than the working memory deficit: A functional near-infrared spectroscopy study

This study tested whether cognitive deficit in patients with adult attention deficit hyperactivity disorder (ADHD) is a working memory deficit or cognitive state disorder during the N-back task. Twenty-two adults with ADHD and twenty-four healthy controls participated in the N-back task. The functional near-infrared spectroscopy (fNIRS) was combined with three perspectives from behavioral and spatial and temporal activation characteristics of blood oxygen levels in the prefrontal cortex to examine the psychological and neuroprocessing characteristics of adult ADHD. Data were acquired using a block design during an N-back task with three memory loads. Visual stimuli were presented on a computer monitor. Behaviorally, response time and accuracy showed no significant differences between the two groups. Spatially, in the left orbitofrontal area and the left frontopolar area (Channels 4 and 11), adult ADHD had significantly higher activation levels of oxyHb in the 2-back task and lower activation levels of deoxyHb in the 3-back task than healthy controls (corrected p < 0.05). Therefore, Channel 4 in the 2-back condition and Channel 11 in the 3-back condition were used as the regions of interest (ROI). Temporally, adults with ADHD peaked earlier in the ROIs than healthy controls. Furthermore, working memory deficit was not found directly from the behavioral performance in adult ADHD. However, adult ADHD can be affected by memory load, task duration, and novelty stimulus. Our findings suggest that patients with adult ADHD have cognitive state disorder instead of working memory deficit.

Deconvolution of hemodynamic responses along the cortical surface using personalized functional near infrared spectroscopy

In functional near infrared spectroscopy (fNIRS), deconvolution analysis of oxy and deoxy-hemoglobin concentration changes allows estimating specific hemodynamic response functions (HRF) elicited by neuronal activity, taking advantage of the fNIRS excellent temporal resolution. Diffuse optical tomography (DOT) is also becoming the new standard reconstruction procedure as it is more accurate than the modified Beer Lambert law approach at the sensor level. The objective of this study was to assess the relevance of HRF deconvolution after DOT constrained along the cortical surface. We used local personalized fNIRS montages which consists in optimizing the position of fNIRS optodes to ensure maximal sensitivity to subject specific target brain regions. We carefully evaluated the accuracy of deconvolution when applied after DOT, using realistic simulations involving several HRF models at different signal to noise ratio (SNR) levels and on real data related to motor and visual tasks in healthy subjects and from spontaneous pathological activity in one patient with epilepsy. We demonstrated that DOT followed by deconvolution was able to accurately recover a large variability of HRFs over a large range of SNRs. We found good performances of deconvolution analysis for SNR levels usually encountered in our applications and we were able to reconstruct accurately the temporal dynamics of HRFs in real conditions.

Joint direct estimation of hemodynamic response function and activation level in brain functional high density diffuse optical tomography

High density diffuse optical tomography has become increasingly important to detect underlying neuronal activities. Conventional methods first estimate the time courses of the changes in the absorption coefficients for all the voxels, and then estimate the hemodynamic response function (HRF). Activation-level maps are extracted at last based on this HRF. However, the error propagation among the successive processes degrades and even misleads the final results. Besides, the computation burden is heavy. To address the above problems, a direct method is proposed in this paper to simultaneously estimate the HRF and the activation-level maps from the boundary fluxes. It is assumed that all the voxels in the same activated brain region share the same HRF but differ in the activation levels, and no prior information is imposed on the specific shape of the HRF. The dynamic simulation and phantom experiments demonstrate that the proposed method outperforms the conventional one in terms of the estimation accuracy and computation speed.

Using the General Linear Model to Improve Performance in fNIRS Single Trial Analysis and Classification: A Perspective

Within a decade, single trial analysis of functional Near Infrared Spectroscopy (fNIRS) signals has gained significant momentum, and fNIRS joined the set of modalities frequently used for active and passive Brain Computer Interfaces (BCI). A great variety of methods for feature extraction and classification have been explored using state-of-the-art Machine Learning methods. In contrast, signal preprocessing and cleaning pipelines for fNIRS often follow simple recipes and so far rarely incorporate the available state-of-the-art in adjacent fields. In neuroscience, where fMRI and fNIRS are established neuroimaging tools, evoked hemodynamic brain activity is typically estimated across multiple trials using a General Linear Model (GLM). With the help of the GLM, subject, channel, and task specific evoked hemodynamic responses are estimated, and the evoked brain activity is more robustly separated from systemic physiological interference using independent measures of nuisance regressors, such as short-separation fNIRS measurements. When correctly applied in single trial analysis, e.g., in BCI, this approach can significantly enhance contrast to noise ratio of the brain signal, improve feature separability and ultimately lead to better classification accuracy. In this manuscript, we provide a brief introduction into the GLM and show how to incorporate it into a typical BCI preprocessing pipeline and cross-validation. Using a resting state fNIRS data set augmented with synthetic hemodynamic responses that provide ground truth brain activity, we compare the quality of commonly used fNIRS features for BCI that are extracted from (1) conventionally preprocessed signals, and (2) signals preprocessed with the GLM and physiological nuisance regressors. We show that the GLM-based approach can provide better single trial estimates of brain activity as well as a new feature type, i.e., the weight of the individual and channel-specific hemodynamic response function (HRF) regressor. The improved estimates yield features with higher separability, that significantly enhance accuracy in a binary classification task when compared to conventional preprocessing—on average +7.4% across subjects and feature types. We propose to adapt this well-established approach from neuroscience to the domain of single-trial analysis and preprocessing wherever the classification of evoked brain activity is of concern, for instance in BCI.

Improved physiological noise regression in fNIRS: A multimodal extension of the General Linear Model using temporally embedded Canonical Correlation Analysis

For the robust estimation of evoked brain activity from functional Near-Infrared Spectroscopy (fNIRS) signals, it is crucial to reduce nuisance signals from systemic physiology and motion. The current best practice incorporates short-separation (SS) fNIRS measurements as regressors in a General Linear Model (GLM). However, several challenging signal characteristics such as non-instantaneous and non-constant coupling are not yet addressed by this approach and additional auxiliary signals are not optimally exploited. We have recently introduced a new methodological framework for the unsupervised multivariate analysis of fNIRS signals using Blind Source Separation (BSS) methods. Building onto the framework, in this manuscript we show how to incorporate the advantages of regularized temporally embedded Canonical Correlation Analysis (tCCA) into the supervised GLM. This approach allows flexible integration of any number of auxiliary modalities and signals. We provide guidance for the selection of optimal parameters and auxiliary signals for the proposed GLM extension. Its performance in the recovery of evoked HRFs is then evaluated using both simulated ground truth data and real experimental data and compared with the GLM with short-separation regression. Our results show that the GLM with tCCA significantly improves upon the current best practice, yielding significantly better results across all applied metrics: Correlation (HbO max. +45%), Root Mean Squared Error (HbO max. −55%), F-Score (HbO up to 3.25-fold) and p-value as well as power spectral density of the noise floor. The proposed method can be incorporated into the GLM in an easily applicable way that flexibly combines any available auxiliary signals into optimal nuisance regressors. This work has potential significance both for conventional neuroscientific fNIRS experiments as well as for emerging applications of fNIRS in everyday environments, medicine and BCI, where high Contrast to Noise Ratio is of importance for single trial analysis.

A Kalman-based tomographic scheme for directly reconstructing activation levels of brain function

In functional near-infrared spectroscopy (fNIRS), the conventional indirect approaches first separately recover the spatial distribution of the changes in the optical properties at every time point, and then extract the activation levels by a time-course analysis process at every site. In the tomographic implementation of fNIRS, i.e., diffuse optical tomography (DOT), these approaches not only suffer from the ill-posedness of the optical inversions and error propagation between the two successive steps, but also fail to achieve satisfactory temporal resolution due to the requirement for a complete data set. To cope with the above adversities of the indirect approaches, we propose herein a direct approach to tomographically reconstructing the activation levels by incorporating a Kalman scheme. Dynamic simulative and phantom experiments were conducted for the performance validation of the proposed approach, demonstrating its potentials to improve the calculated images and to relax the speed limitation of the instruments.

A miniature multi-contrast microscope for functional imaging in freely behaving animals

Neurovascular coupling, cerebrovascular remodeling and hemodynamic changes are critical to brain function, and dysregulated in neuropathologies such as brain tumors. Interrogating these phenomena in freely behaving animals requires a portable microscope with multiple optical contrast mechanisms. Therefore, we developed a miniaturized microscope with: a fluorescence (FL) channel for imaging neural activity (e.g., GCaMP) or fluorescent cancer cells (e.g., 9L-GFP); an intrinsic optical signal (IOS) channel for imaging hemoglobin absorption (i.e., cerebral blood volume); and a laser speckle contrast (LSC) channel for imaging perfusion (i.e., cerebral blood flow). Following extensive validation, we demonstrate the microscope’s capabilities via experiments in unanesthetized murine brains that include: (i) multi-contrast imaging of neurovascular changes following auditory stimulation; (ii) wide-area tonotopic mapping; (iii) EEG-synchronized imaging during anesthesia recovery; and (iv) microvascular connectivity mapping over the life-cycle of a brain tumor. This affordable, flexible, plug-and-play microscope heralds a new era in functional imaging of freely behaving animals.

Measuring multiple neurophysiologic variables usually requires bulky benchtop optical systems and working with anesthetized animals. Here the authors present a miniature portable microscope for neurovascular imaging in awake rodents, combining fluorescence, intrinsic optical signals and laser speckle contrast.

Multichannel continuous electroencephalography-functional near-infrared spectroscopy recording of focal seizures and interictal epileptiform discharges in human epilepsy: a review

Functional near-infrared spectroscopy (fNIRS) has emerged as a promising neuroimaging technique as it allows noninvasive and long-term monitoring of cortical hemodynamics. Recent work by our group and others has revealed the potential of fNIRS, combined with electroencephalography (EEG), in the context of human epilepsy. Hemodynamic brain responses attributed to epileptic events, such as seizures and interictal epileptiform discharges (IEDs), are routinely observed with a good degree of statistical significance and in concordance with clinical presentation. Recording done with over 100 channels allows sufficiently large coverage of the epileptic focus and other areas. Three types of seizures have been documented: frontal lobe seizures, temporal lobe seizures, and posterior seizures. Increased oxygenation was observed in the epileptic focus in most cases, while rapid but similar hemodynamic variations were identified in the contralateral homologous region. While investigating IEDs, it was shown that their hemodynamic effect is observable with fNIRS, that their response is associated with significant (inhibitive) nonlinearities, and that the sensitivity and specificity of fNIRS to localize the epileptic focus can be estimated in a sample of 40 patients. This paper first reviews recent EEG-fNIRS developments in epilepsy research and then describes applications to the study of focal seizures and IEDs.

Modeling Neurovascular Coupling from Clustered Parameter Sets for Multimodal EEG-NIRS

Despite significant improvements in neuroimaging technologies and analysis methods, the fundamental relationship between local changes in cerebral hemodynamics and the underlying neural activity remains largely unknown. In this study, a data driven approach is proposed for modeling this neurovascular coupling relationship from simultaneously acquired electroencephalographic (EEG) and near-infrared spectroscopic (NIRS) data. The approach uses gamma transfer functions to map EEG spectral envelopes that reflect time-varying power variations in neural rhythms to hemodynamics measured with NIRS during median nerve stimulation. The approach is evaluated first with simulated EEG-NIRS data and then by applying the method to experimental EEG-NIRS data measured from 3 human subjects. Results from the experimental data indicate that the neurovascular coupling relationship can be modeled using multiple sets of gamma transfer functions. By applying cluster analysis, statistically significant parameter sets were found to predict NIRS hemodynamics from EEG spectral envelopes. All subjects were found to have significant clustered parameters (P < 0.05) for EEG-NIRS data fitted using gamma transfer functions. These results suggest that the use of gamma transfer functions followed by cluster analysis of the resulting parameter sets may provide insights into neurovascular coupling in human neuroimaging data.

An integrated framework for joint HRF and drift estimation and HbO/HbR signal improvement in fNIRS data

Nonparametric hemodynamic response function (HRF) estimation in functional near-infrared spectroscopy (fNIRS) data plays an important role when investigating the temporal dynamics of a brain region response during activations. Assuming the drift arising from both physical and physiological effects in fNIRS data is Lipschitz continuous; a novel algorithm for joint HRF and drift estimation is derived in this paper. The proposed algorithm estimates the HRF by applying a first-order differencing to the fNIRS time series samples in order to remove the drift effect. An estimate of the drift is then obtained using a wavelet thresholding technique applied to the residuals generated by removing the estimated induced activation response from the fNIRS time-series. It is shown that the proposed HRF estimator is √N consistent whereas the estimator of the drift is asymptotically optimal. The de-drifted fNIRS oxygenated (HbO) and deoxygenated (HbR) hemoglobin responses are then obtained by removing the corresponding estimated drifts from the fNIRS time-series. Its performance is assessed using both simulated and real fNIRS data sets. The application results reveal that the proposed joint HRF and drift estimation method is efficient both computationally and in terms of accuracy. In comparison to traditional model based methods used for HRF estimation, the proposed novel method avoids the selection of a model to remove the drift component. As a result, the proposed method finds an optimal estimate of the fNIRS drift and offers a model-free approach to de-drift the HbO/HbR responses.

Optimizing the general linear model for functional near-infrared spectroscopy: an adaptive hemodynamic response function approach

An increasing number of functional near-infrared spectroscopy (fNIRS) studies utilize a general linear model (GLM) approach, which serves as a standard statistical method for functional magnetic resonance imaging (fMRI) data analysis. While fMRI solely measures the blood oxygen level dependent (BOLD) signal, fNIRS measures the changes of oxy-hemoglobin (oxy-Hb) and deoxy-hemoglobin (deoxy-Hb) signals at a temporal resolution severalfold higher. This suggests the necessity of adjusting the temporal parameters of a GLM for fNIRS signals. Thus, we devised a GLM-based method utilizing an adaptive hemodynamic response function (HRF). We sought the optimum temporal parameters to best explain the observed time series data during verbal fluency and naming tasks. The peak delay of the HRF was systematically changed to achieve the best-fit model for the observed oxy- and deoxy-Hb time series data. The optimized peak delay showed different values for each Hb signal and task. When the optimized peak delays were adopted, the deoxy-Hb data yielded comparable activations with similar statistical power and spatial patterns to oxy-Hb data. The adaptive HRF method could suitably explain the behaviors of both Hb parameters during tasks with the different cognitive loads during a time course, and thus would serve as an objective method to fully utilize the temporal structures of all fNIRS data.

State-space models of impulse hemodynamic responses over motor, somatosensory, and visual cortices

THE PAPER PRESENTS STATE SPACE MODELS OF THE HEMODYNAMIC RESPONSE (HR) OF FNIRS TO AN IMPULSE STIMULUS IN THREE BRAIN REGIONS: motor cortex (MC), somatosensory cortex (SC), and visual cortex (VC). Nineteen healthy subjects were examined. For each cortex, three impulse HRs experimentally obtained were averaged. The averaged signal was converted to a state space equation by using the subspace method. The activation peak and the undershoot peak of the oxy-hemoglobin (HbO) in MC are noticeably higher than those in SC and VC. The time-to-peaks of the HbO in three brain regions are almost the same (about 6.76 76 ± 0.2 s). The time to undershoot peak in VC is the largest among three. The HbO decreases in the early stage (~0.46 s) in MC and VC, but it is not so in SC. These findings were well described with the developed state space equations. Another advantage of the proposed method is its easy applicability in generating the expected HR to arbitrary stimuli in an online (or real-time) imaging. Experimental results are demonstrated.

Hemodynamic changes during posterior epilepsies: an EEG-fNIRS study

Posterior epilepsies are mainly characterized clinically by visual symptoms. Functional near-infrared spectroscopy (fNIRS) is an emerging non-invasive imaging technique that has the potential to monitor hemodynamic changes during epileptic activity. Combined with electroencephalography (EEG), 9 patients with posterior epilepsies were recorded using EEG-fNIRS with large sampling (19 EEG electrodes and over 100 fNIRS channels). Spikes and seizures were carefully marked on EEG traces, and convolved with a standard hemodynamic response function for general linear model (GLM) analysis. GLM results for seizures (in 3 patients) and spikes (7 patients) were broadly sensitive to the epileptic focus in 7/9 patients, and specific in 5/9 patients with fNIRS deoxyhemoglobin responses lateralized to the correct lobe, and to plausible locations within the occipital or parietal lobes. This work provides evidence that EEG-fNIRS is a sensitive technique for monitoring posterior epileptic activity.

Developmental changes in brain activation and functional connectivity during response inhibition in the early childhood brain

Response inhibition is an attention function which develops relatively early during childhood. Behavioral data suggest that by the age of 3, children master the basic task requirements for the assessment of response inhibition but performance improves substantially until the age of 7. The neuronal mechanisms underlying these developmental processes, however, are not well understood. In this study, we examined brain activation patterns and behavioral performance of children aged between 4 and 6 years compared to adults by applying a go/no-go paradigm during near-infrared spectroscopy (NIRS) brain imaging. We furthermore applied task-independent functional connectivity measures to the imaging data to identify maturation of intrinsic neural functional networks. We found a significant group×condition related interaction in terms of inhibition-related reduced right fronto-parietal activation in children compared to adults. In contrast, motor-related activation did not differ between age groups. Functional connectivity analysis revealed that in the children's group, short-range coherence within frontal areas was stronger, and long-range coherence between frontal and parietal areas was weaker, compared to adults. Our findings show that in children aged from 4 to 6 years fronto-parietal brain maturation plays a crucial part in the cognitive development of response inhibition.

Statistical analysis of high density diffuse optical tomography

High density diffuse optical tomography (HD-DOT) is a noninvasive neuroimaging modality with moderate spatial resolution and localization accuracy. Due to portability and wear-ability advantages, HD-DOT has the potential to be used in populations that are not amenable to functional magnetic resonance imaging (fMRI), such as hospitalized patients and young children. However, whereas the use of event-related stimuli designs, general linear model (GLM) analysis, and imaging statistics are standardized and routine with fMRI, such tools are not yet common practice in HD-DOT. In this paper we adapt and optimize fundamental elements of fMRI analysis for application to HD-DOT. We show the use of event-related protocols and GLM de-convolution analysis in un-mixing multi-stimuli event-related HD-DOT data. Statistical parametric mapping (SPM) in the framework of a general linear model is developed considering the temporal and spatial characteristics of HD-DOT data. The statistical analysis utilizes a random field noise model that incorporates estimates of the local temporal and spatial correlations of the GLM residuals. The multiple-comparison problem is addressed using a cluster analysis based on non-stationary Gaussian random field theory. These analysis tools provide access to a wide range of experimental designs necessary for the study of the complex brain functions. In addition, they provide a foundation for understanding and interpreting HD-DOT results with quantitative estimates for the statistical significance of detected activation foci.

A fNIRS investigation of switching and inhibition during the modified Stroop task in younger and older adults

Brain imaging studies have reported age-related differences in brain activation for attentional control functions, such as inhibition and task-switching. However, age-related differences in brain activation patterns in more than one attentional control task have rarely been studied in the same group of participants. In this study, younger and older adults completed a modified Stroop task with interference and switching conditions, using functional near infra-red spectroscopy. While interference did not reveal any significant activation of the prefrontal cortex in younger adults, switching produced an increased activation bilaterally in both the anterior dorsolateral prefrontal cortex (DLPFC) and the anterior ventrolateral prefrontal cortex (VLPFC). In older adults, an isolated right and left anterior DLPFC activation was observed even in the non-executive conditions of the Stroop task (color denomination) and the interference condition revealed activation mostly in the posterior left DLPFC and bilateral VLPFC with a small right anterior DLPFC component. Specific to older adults, switching induced an increased activation spread out bilaterally over the prefrontal cortex in the bilateral anterior DLPFC, the posterior left DLPFC and bilateral VLPFC. These results suggest that for both older and younger adults, inhibition and switching are associated with distinct patterns of prefrontal activation and that age-related differences exist in these patterns such that prefrontal activation seems to be more spread out at different sites in older adults.

Load-dependent brain activation assessed by time-domain functional near-infrared spectroscopy during a working memory task with graded levels of difficulty

We evaluated frontal brain activation during a mixed attentional/working memory task with graded levels of difficulty in a group of 19 healthy subjects, by means of time-domain functional near-infrared spectroscopy (fNIRS). Brain activation was assessed, and load-related oxy- and deoxy-hemoglobin changes were studied. Generalized linear model (GLM) was applied to the data to explore the metabolic processes occurring during the mental effort and, possibly, their involvement in short-term memorization. GLM was applied to the data twice: for modeling the task as a whole and for specifically investigating brain activation at each cognitive load. This twofold employment of GLM allowed (1) the extraction and isolation of different information from the same signals, obtained through the modeling of different cognitive categories (sustained attention and working memory), and (2) the evaluation of model fitness, by inspection and comparison of residuals (i.e., unmodeled part of the signal) obtained in the two different cases. Results attest to the presence of a persistent attentional-related metabolic activity, superimposed to a task-related mnemonic contribution. Some hemispherical differences have also been highlighted frontally: deoxy-hemoglobin changes manifested a strong right lateralization, whereas modifications in oxy- and total hemoglobin showed a medial localization. The present work successfully explored the capability of fNIRS to detect the two neurophysiological categories under investigation and distinguish their activation patterns.

Evaluation of receptive field size from higher harmonics in visuotopic mapping using continuous stimulation optical imaging

The extraction and detection of specific responses from a large amount of background noise has been the subject of a considerable body of research in brain functional imaging, and more specifically in optical intrinsic signal imaging. Recent work by Kalatsky and Stryker (2003) showed that by combining different conditions and using periodic stimuli, recording times can be reduced. Spectral decomposition is then used to provide amplitude and phase information locked to the stimulus. A drawback of the above method is that by focusing only on a single harmonic, response information is limited. The shape of the hemodynamic response function (HRF) and the temporal variations in the neural responses cannot be assessed. In this work it is argued that additional information about neural responses can be gathered by using higher harmonics. Moving bars were used to generate visuotopic maps on large portions of the cat visual cortex. Up to four simultaneously bars moving repetitively across the visual field at different frequencies were used to sample the HRF in the Fourier domain. The HRF profile obtained with continuous stimulation was spatially homogeneous throughout the cortex and similar to the HRF profile obtained using episodic stimulation. Furthermore, by modeling the optical response as a convolution between HRF and neuronal responses, the ratio of the second harmonic to the first provided an estimation of the receptive field size. This was further validated by measuring spatial frequency selectivity. Therefore, the use of higher harmonics opens new avenues to estimate receptive field size from temporal signals.

Intact neurovascular coupling during executive function in migraine without aura: Interictal near-infrared spectroscopy study

An altered neurovascular coupling has been proposed in migraine. We aimed to investigate neurovascular coupling during a mental task interictally in patients with migraine without aura (MO) by near-infrared spectroscopy (NIRS). Twelve migraineurs and 12 healthy controls were included. Using NIRS, we recorded the magnitude and latency of cortical changes in oxyhaemoglobin (HbO2) and deoxyhaemoglobin (Hb) during the colour-word matching Stroop test via 16 channels covering the forehead. We found no differences in the magnitude of responses between migraineurs and healthy subjects in the incongruent Stroop task subtracted by the neutral Stroop task on either side of the frontal cortex for HbO2 (left, P = 0.984; right, P = 0.406) or Hb (left, P = 0.689; right, P = 0.406) values. No differences in error rate ( P = 0.611) or reaction time ( P = 0.936) were found between healthy subjects and MO patients for incongruent tasks. The present study suggests that vascular reactivity and oxygen supply during a mental task in patients with MO are intact interictally.

Real-time imaging of human brain function by near-infrared spectroscopy using an adaptive general linear model

Near-infrared spectroscopy is a non-invasive neuroimaging method which uses light to measure changes in cerebral blood oxygenation associated with brain activity. In this work, we demonstrate the ability to record and analyze images of brain activity in real-time using a 16-channel continuous wave optical NIRS system. We propose a novel real-time analysis framework using an adaptive Kalman filter and a state-space model based on a canonical general linear model of brain activity. We show that our adaptive model has the ability to estimate single-trial brain activity events as we apply this method to track and classify experimental data acquired during an alternating bilateral self-paced finger tapping task.

HomER: a review of time-series analysis methods for near-infrared spectroscopy of the brain

Near-infrared spectroscopy (NIRS) is a noninvasive neuroimaging tool for studying evoked hemodynamic changes within the brain. By this technique, changes in the optical absorption of light are recorded over time and are used to estimate the functionally evoked changes in cerebral oxyhemoglobin and deoxyhemoglobin concentrations that result from local cerebral vascular and oxygen metabolic effects during brain activity. Over the past three decades this technology has continued to grow, and today NIRS studies have found many niche applications in the fields of psychology, physiology, and cerebral pathology. The growing popularity of this technique is in part associated with a lower cost and increased portability of NIRS equipment when compared with other imaging modalities, such as functional magnetic resonance imaging and positron emission tomography. With this increasing number of applications, new techniques for the processing, analysis, and interpretation of NIRS data are continually being developed. We review some of the time-series and functional analysis techniques that are currently used in NIRS studies, we describe the practical implementation of various signal processing techniques for removing physiological, instrumental, and motion-artifact noise from optical data, and we discuss the unique aspects of NIRS analysis in comparison with other brain imaging modalities. These methods are described within the context of the MATLAB-based graphical user interface program, HomER, which we have developed and distributed to facilitate the processing of optical functional brain data.

1/f noise in diffuse optical imaging and wavelet-based response estimation

In diffuse optical imaging (DOI) data analysis, the functional response is contaminated with physiological noise as in functional magnetic resonance imaging (fMRI). In this work we extend a previously proposed method for fMRI to estimate the parameters of a linear model of DOI time series. The regression is performed in the wavelet domain to infer drift coefficients at different scales and to estimate the strength of the hemodynamic response function (HRF). This multiresolution approach benefits from the whitening property of the discrete wavelet transform (DWT), which approximately decorrelates long-memory noise processes. We also show that a more accurate estimation is obtained by removing some regressors correlating with the protocol. Moreover, we observe that this improvement is related to a quantitative measure of 1/f noise. The performances of the method are first evaluated against a standard spline-cosine drift approach with simulated HRF and real background physiology. Lastly, the technique is applied to experimental event-related data acquired by near-infrared spectroscopy (NIRS).

Direct estimation of evoked hemoglobin changes by multimodality fusion imaging

In the last two decades, both diffuse optical tomography (DOT) and blood oxygen level dependent (BOLD)-based functional magnetic resonance imaging (fMRI) methods have been developed as noninvasive tools for imaging evoked cerebral hemodynamic changes in studies of brain activity. Although these two technologies measure functional contrast from similar physiological sources, i.e., changes in hemoglobin levels, these two modalities are based on distinct physical and biophysical principles leading to both limitations and strengths to each method. In this work, we describe a unified linear model to combine the complimentary spatial, temporal, and spectroscopic resolutions of concurrently measured optical tomography and fMRI signals. Using numerical simulations, we demonstrate that concurrent optical and BOLD measurements can be used to create cross-calibrated estimates of absolute micromolar deoxyhemoglobin changes. We apply this new analysis tool to experimental data acquired simultaneously with both DOT and BOLD imaging during a motor task, demonstrate the ability to more robustly estimate hemoglobin changes in comparison to DOT alone, and show how this approach can provide cross-calibrated estimates of hemoglobin changes. Using this multimodal method, we estimate the calibration of the 3 tesla BOLD signal to be -0.55%+/-0.40% signal change per micromolar change of deoxyhemoglobin.