Dynamic physiological modeling for functional diffuse optical tomography

Short-separation regression incorporated diffuse optical tomography image reconstruction modeling for high-density functional near-infrared spectroscopy

Significance

Short-separation (SS) regression and diffuse optical tomography (DOT) image reconstruction, two widely adopted methods in functional near-infrared spectroscopy (fNIRS), were demonstrated to individually facilitate the separation of brain activation and physiological signals, with further improvement using both sequentially. We hypothesized that doing both simultaneously would further improve the performance.

Aim

Motivated by the success of these two approaches, we propose a method, SS-DOT, which applies SS and DOT simultaneously.

Approach

The method, which employs spatial and temporal basis functions to represent the hemoglobin concentration changes, enables us to incorporate SS regressors into the time series DOT model. To benchmark the performance of the SS-DOT model against conventional sequential models, we use fNIRS resting state data augmented with synthetic brain response as well as data acquired during a ball squeezing task. The conventional sequential models comprise performing SS regression and DOT.

Results

The results show that the SS-DOT model improves the image quality by increasing the contrast-to-background ratio by a threefold improvement. The benefits are marginal at small brain activation.

Conclusions

The SS-DOT model improves the fNIRS image reconstruction quality.

Simultaneous multimodal fNIRS-EEG recordings reveal new insights in neural activity during motor execution, observation, and imagery

Motor execution, observation, and imagery are important skills used in motor learning and rehabilitation. The neural mechanisms underlying these cognitive-motor processes are still poorly understood. We used a simultaneous recording of functional near-infrared spectroscopy (fNIRS) and electroencephalogram (EEG) to elucidate the differences in neural activity across three conditions requiring these processes. Additionally, we used a new method called structured sparse multiset Canonical Correlation Analysis (ssmCCA) to fuse the fNIRS and EEG data and determine the brain regions of neural activity consistently detected by both modalities. Unimodal analyses revealed differentiated activation between conditions; however, the activated regions did not fully overlap across the two modalities (fNIRS: left angular gyrus, right supramarginal gyrus, as well as right superior and inferior parietal lobes; EEG: bilateral central, right frontal, and parietal). These discrepancies might be because fNIRS and EEG detect different signals. Using fused fNIRS-EEG data, we consistently found activation over the left inferior parietal lobe, superior marginal gyrus, and post-central gyrus during all three conditions, suggesting that our multimodal approach identifies a shared neural region associated with the Action Observation Network (AON). This study highlights the strengths of using the multimodal fNIRS-EEG fusion technique for studying AON. Neural researchers should consider using the multimodal approach to validate their findings.

How much do time-domain functional near-infrared spectroscopy (fNIRS) moments improve estimation of brain activity over traditional fNIRS?

SIGNIFICANCE

Advances in electronics have allowed the recent development of compact, high channel count time domain functional near-infrared spectroscopy (TD-fNIRS) systems. Temporal moment analysis has been proposed for increased brain sensitivity due to the depth selectivity of higher order temporal moments. We propose a general linear model (GLM) incorporating TD moment data and auxiliary physiological measurements, such as short separation channels, to improve the recovery of the HRF.

AIMS

We compare the performance of previously reported multi-distance TD moment techniques to commonly used techniques for continuous wave (CW) fNIRS hemodynamic response function (HRF) recovery, namely block averaging and CW GLM. Additionally, we compare the multi-distance TD moment technique to TD moment GLM.

APPROACH

We augmented resting TD-fNIRS moment data (six subjects) with known synthetic HRFs. We then employed block averaging and GLM techniques with "short-separation regression" designed both for CW and TD to recover the HRFs. We calculated the root mean square error (RMSE) and the correlation of the recovered HRF to the ground truth. We compared the performance of equivalent CW and TD techniques with paired t-tests.

RESULTS

We found that, on average, TD moment HRF recovery improves correlations by 98% and 48% for HbO and HbR respectively, over CW GLM. The improvement on the correlation for TD GLM over TD moment is 12% (HbO) and 27% (HbR). RMSE decreases 56% and 52% (HbO and HbR) for TD moment compared to CW GLM. We found no statistically significant improvement in the RMSE for TD GLM compared to TD moment.

CONCLUSIONS

Properly covariance-scaled TD moment techniques outperform their CW equivalents in both RMSE and correlation in the recovery of the synthetic HRFs. Furthermore, our proposed TD GLM based on moments outperforms regular TD moment analysis, while allowing the incorporation of auxiliary measurements of the confounding physiological signals from the scalp.

Lawrence K, Tachtsidis I, Tong Y, Torricelli A, Urner T, Wabnitz H, Wolf M, Wolf U, Xu S, Yang C, Yodh AG, Yücel MA, Zhou W. Optical imaging and spectroscopy for the study of the human brain: status report

This report is the second part of a comprehensive two-part series aimed at reviewing an extensive and diverse toolkit of novel methods to explore brain health and function. While the first report focused on neurophotonic tools mostly applicable to animal studies, here, we highlight optical spectroscopy and imaging methods relevant to noninvasive human brain studies. We outline current state-of-the-art technologies and software advances, explore the most recent impact of these technologies on neuroscience and clinical applications, identify the areas where innovation is needed, and provide an outlook for the future directions.

A Model-Based Iterative Learning Approach for Diffuse Optical Tomography

Diffuse optical tomography (DOT) utilises near-infrared light for imaging spatially distributed optical parameters, typically the absorption and scattering coefficients. The image reconstruction problem of DOT is an ill-posed inverse problem, due to the non-linear light propagation in tissues and limited boundary measurements. The ill-posedness means that the image reconstruction is sensitive to measurement and modelling errors. The Bayesian approach for the inverse problem of DOT offers the possibility of incorporating prior information about the unknowns, rendering the problem less ill-posed. It also allows marginalisation of modelling errors utilising the so-called Bayesian approximation error method. A more recent trend in image reconstruction techniques is the use of deep learning, which has shown promising results in various applications from image processing to tomographic reconstructions. In this work, we study the non-linear DOT inverse problem of estimating the (absolute) absorption and scattering coefficients utilising a 'model-based' learning approach, essentially intertwining learned components with the model equations of DOT. The proposed approach was validated with 2D simulations and 3D experimental data. We demonstrated improved absorption and scattering estimates for targets with a mix of smooth and sharp image features, implying that the proposed approach could learn image features that are difficult to model using standard Gaussian priors. Furthermore, it was shown that the approach can be utilised in compensating for modelling errors due to coarse discretisation enabling computationally efficient solutions. Overall, the approach provided improved computation times compared to a standard Gauss-Newton iteration.

Frequency-domain analysis of fNIRS fluctuations induced by rhythmic mental arithmetic

Functional near‐infrared spectroscopy (fNIRS) is an increasingly used technology for imaging neural correlates of cognitive processes. However, fNIRS signals are commonly impaired by task‐evoked and spontaneous hemodynamic oscillations of non‐cerebral origin, a major challenge in fNIRS research. In an attempt to isolate the task‐evoked cortical response, we investigated the coupling between hemodynamic changes arising from superficial and deep layers during mental effort. For this aim, we applied a rhythmic mental arithmetic task to induce cyclic hemodynamic fluctuations suitable for effective frequency‐resolved measurements. Twenty university students aged 18–25 years (eight males) underwent the task while hemodynamic changes were monitored in the forehead using a newly developed NIRS device, capable of multi‐channel and multi‐distance recordings. We found significant task‐related fluctuations for oxy‐ and deoxy‐hemoglobin, highly coherent across shallow and deep tissue layers, corroborating the strong influence of surface hemodynamics on deep fNIRS signals. Importantly, after removing such surface contamination by linear regression, we show that the frontopolar cortex response to a mental math task follows an unusual inverse oxygenation pattern. We confirm this finding by applying for the first time an alternative method to estimate the neural signal, based on transfer function analysis and phasor algebra. Altogether, our results demonstrate the feasibility of using a rhythmic mental task to impose an oscillatory state useful to separate true brain functional responses from those of non‐cerebral origin. This separation appears to be essential for a better understanding of fNIRS data and to assess more precisely the dynamics of the neuro‐visceral link.

We proposed the use of rhythmic mental arithmetic tasks to induce cyclic oscillations in multi‐distance fNIRS signals measured on the forehead, suitable for effective frequency‐domain analysis to better identify the actual neural functional response. We confirm the impairment of deep signals by task‐evoked non‐cerebral confounds, while providing evidence for an inverse oxygenation response in the frontopolar cortex.

Multivariate Kalman filter regression of confounding physiological signals for real-time classification of fNIRS data

Significance: Functional near-infrared spectroscopy (fNIRS) is a noninvasive technique for measuring hemodynamic changes in the human cortex related to neural function. Due to its potential for miniaturization and relatively low cost, fNIRS has been proposed for applications, such as brain-computer interfaces (BCIs). The relatively large magnitude of the signals produced by the extracerebral physiology compared with the ones produced by evoked neural activity makes real-time fNIRS signal interpretation challenging. Regression techniques incorporating physiologically relevant auxiliary signals such as short separation channels are typically used to separate the cerebral hemodynamic response from the confounding components in the signal. However, the coupling of the extra-cerebral signals is often noninstantaneous, and it is necessary to find the proper delay to optimize nuisance removal. Aim: We propose an implementation of the Kalman filter with time-embedded canonical correlation analysis for the real-time regression of fNIRS signals with multivariate nuisance regressors that take multiple delays into consideration. Approach: We tested our proposed method on a previously acquired finger tapping dataset with the purpose of classifying the neural responses as left or right. Results: We demonstrate computationally efficient real-time processing of 24-channel fNIRS data (400 samples per second per channel) with a two order of selective magnitude decrease in cardiac signal power and up to sixfold increase in the contrast-to-noise ratio compared with the nonregressed signals. Conclusion: The method provides a way to obtain better distinction of brain from non-brain signals in real time for BCI application with fNIRS.

Evaluation of fNIRS signal components elicited by cognitive and hypercapnic stimuli

Functional near infrared spectroscopy (fNIRS) measurements are confounded by signal components originating from multiple physiological causes, whose activities may vary temporally and spatially (across tissue layers, and regions of the cortex). Furthermore, the stimuli can induce evoked effects, which may lead to over or underestimation of the actual effect of interest. Here, we conducted a temporal, spectral, and spatial analysis of fNIRS signals collected during cognitive and hypercapnic stimuli to characterize effects of functional versus systemic responses. We utilized wavelet analysis to discriminate physiological causes and employed long and short source-detector separation (SDS) channels to differentiate tissue layers. Multi-channel measures were analyzed further to distinguish hemispheric differences. The results highlight cardiac, respiratory, myogenic, and very low frequency (VLF) activities within fNIRS signals. Regardless of stimuli, activity within the VLF band had the largest contribution to the overall signal. The systemic activities dominated the measurements from the short SDS channels during cognitive stimulus, but not hypercapnic stimulus. Importantly, results indicate that characteristics of fNIRS signals vary with type of the stimuli administered as cognitive stimulus elicited variable responses between hemispheres in VLF band and task-evoked temporal effect in VLF, myogenic and respiratory bands, while hypercapnic stimulus induced a global response across both hemispheres.

Investigating Language and Domain-General Processing in Neurotypicals and Individuals With Aphasia - A Functional Near-Infrared Spectroscopy Pilot Study

Brain reorganization patterns associated with language recovery after stroke have long been debated. Studying mechanisms of spontaneous and treatment-induced language recovery in post-stroke aphasia requires a network-based approach given the potential for recruitment of perilesional left hemisphere language regions, homologous right hemisphere language regions, and/or spared bilateral domain-general regions. Recent hardware, software, and methodological advances in functional near-infrared spectroscopy (fNIRS) make it well-suited to examine this question. fNIRS is cost-effective with minimal contraindications, making it a robust option to monitor treatment-related brain activation changes over time. Establishing clear activation patterns in neurotypical adults during language and domain-general cognitive processes via fNIRS is an important first step. Some fNIRS studies have investigated key language processes in healthy adults, yet findings are challenging to interpret in the context of methodological limitations. This pilot study used fNIRS to capture brain activation during language and domain-general processing in neurotypicals and individuals with aphasia. These findings will serve as a reference when interpreting treatment-related changes in brain activation patterns in post-stroke aphasia in the future. Twenty-four young healthy controls, seventeen older healthy controls, and six individuals with left hemisphere stroke-induced aphasia completed two language tasks (i.e., semantic feature, picture naming) and one domain-general cognitive task (i.e., arithmetic) twice during fNIRS. The probe covered bilateral frontal, parietal, and temporal lobes and included short-separation detectors for scalp signal nuisance regression. Younger and older healthy controls activated core language regions during semantic feature processing (e.g., left inferior frontal gyrus pars opercularis) and lexical retrieval (e.g., left inferior frontal gyrus pars triangularis) and domain-general regions (e.g., bilateral middle frontal gyri) during hard versus easy arithmetic as expected. Consistent with theories of post-stroke language recovery, individuals with aphasia activated areas outside the traditional networks: left superior frontal gyrus and left supramarginal gyrus during semantic feature judgment; left superior frontal gyrus and right precentral gyrus during picture naming; and left inferior frontal gyrus pars opercularis during arithmetic processing. The preliminary findings in the stroke group highlight the utility of using fNIRS to study language and domain-general processing in aphasia.

Data Processing in Functional Near-Infrared Spectroscopy (fNIRS) Motor Control Research

FNIRS pre-processing and processing methodologies are very important-how a researcher chooses to process their data can change the outcome of an experiment. The purpose of this review is to provide a guide on fNIRS pre-processing and processing techniques pertinent to the field of human motor control research. One hundred and twenty-three articles were selected from the motor control field and were examined on the basis of their fNIRS pre-processing and processing methodologies. Information was gathered about the most frequently used techniques in the field, which included frequency cutoff filters, wavelet filters, smoothing filters, and the general linear model (GLM). We discuss the methodologies of and considerations for these frequently used techniques, as well as those for some alternative techniques. Additionally, general considerations for processing are discussed.

Diffuse optical tomography in the human brain: A briefly review from the neurophysiology to its applications

The present work describes the use of noninvasive diffuse optical tomography (DOT) technology to measure hemodynamic changes, providing relevant information which helps to understand the basis of neurophysiology in the human brain. Advantages such as portability, direct measurements of hemoglobin state, temporal resolution, non‐restricted movements as occurs in magnetic resonance imaging (MRI) devices mean that DOT technology can be used in research and clinical fields. In this review we covered the neurophysiology, physical principles underlying optical imaging during tissue‐light interactions, and technology commonly used during the construction of a DOT device including the source‐detector requirements to improve the image quality. DOT provides 3D cerebral activation images due to complex mathematical models which describe the light propagation inside the tissue head. Moreover, we describe briefly the use of Bayesian methods for raw DOT data filtering as an alternative to linear filters widely used in signal processing, avoiding common problems such as the filter selection or a false interpretation of the results which is sometimes due to the interference of background physiological noise with neural activity.

Quantitative comparison of correction techniques for removing systemic physiological signal in functional near-infrared spectroscopy studies

Significance: Isolating task-evoked brain signals from background physiological noise (e.g., cardiac, respiratory, and blood pressure fluctuations) poses a major challenge for the analysis of functional near-infrared spectroscopy (fNIRS) data. Aim: The performance of several analytic methods to separate background physiological noise from brain activity including spatial and temporal filtering, regression, component analysis, and the use of short-separation (SS) measurements were quantitatively compared. Approach: Using experimentally recorded background signals (breath-hold task), receiver operating characteristics simulations were performed by adding various levels of additive synthetic "brain" responses in order to examine the sensitivity and specificity of several previously proposed analytic approaches. Results: We found that the use of SS fNIRS channels as regressors of no-interest within a linear regression model was the best performing approach examined. Furthermore, we found that the addition of all available SS data, including all recorded channels and both hemoglobin species, improved the method performance despite the additional degrees-of-freedom of the models. When SS data were not available, we found that principal component filtering using a separate baseline scan was the best alternative. Conclusions: The use of multiple SS measurements as regressors of no interest implemented in a robust, iteratively prewhitened, general linear model has the best performance of the tested existing methods.

Bibliometric evaluation of 2000-2019 publications on functional near-infrared spectroscopy

This study aimed to explore and analyze research trends and frontiers on functional near-infrared spectroscopy (fNIRS) in the past 20 years and identify collaboration networks. fNIRS-related publications from 2000 to 2019 were retrieved from the Web of Science database. A total of 1727 publications satisfied the search criteria. Bibliometric visualization analysis of active authors, journals, institutions, countries, references, and keywords were conducted. The number of annual related publications remarkably increased over the years. Fallgatter published the largest number of fNIRS-related papers (83). Neuroimage not only had the largest number of papers published in the first 10 journals (157 articles) but also had the highest impact factor (IF, 2018 ​= ​5.812). The University of Tubingen had the highest number of fNIRS-related publications in the past 20 years. The United States ranked first in terms of comprehensive influence in this field. In recent years, burst keywords (e.g., infant, social interaction, and older adult) and a series of references with citation burst provided clues on research frontiers.

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.

Enhancing fNIRS Analysis Using EEG Rhythmic Signatures: An EEG-Informed fNIRS Analysis Study

Neurovascular coupling represents the relationship between changes in neuronal activity and cerebral hemodynamics. Concurrent Electroencephalography (EEG) and functional near-infrared spectroscopy (fNIRS) recording and integration analysis has emerged as a promising multi-modal neuroimaging approach to study the neurovascular coupling as it provides complementary properties with regard to high temporal and moderate spatial resolution of brain activity. In this study we developed an EEG-informed-fNIRS analysis framework to investigate the neuro-correlate between neuronal activity and cerebral hemodynamics by identifying specific EEG rhythmic modulations which contribute to the improvement of the fNIRS-based general linear model (GLM) analysis. Specifically, frequency-specific regressors derived from EEG were used to construct design matrices to guide the GLM analysis of the fNIRS signals collected during a hand grasp task. Our results showed that the EEG-informed fNIRS GLM analysis, especially the alpha and beta band, revealed significantly higher sensitivity and specificity in localizing the task-evoked regions compared to the canonical boxcar model, demonstrating the strong correlations between hemodynamic response and EEG rhythmic modulations. Results also indicated that analysis based on the deoxygenated hemoglobin (HbR) signal slightly outperformed the oxygenated hemoglobin (HbO)-based analysis. The findings in our study not only validate the feasibility of enhancing fNIRS GLM analysis using simultaneously recorded EEG signals, but also provide a new perspective to study the neurovascular coupling of brain activity.

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.

Association of Residents' Neural Signatures With Stress Resilience During Surgery

Importance

Intraoperative stressors may compound cognitive load, prompting performance decline and threatening patient safety. However, not all surgeons cope equally well with stress, and the disparity between performance stability and decline under high cognitive demand may be characterized by differences in activation within brain areas associated with attention and concentration such as the prefrontal cortex (PFC).

Objective

To compare PFC activation between surgeons demonstrating stable performance under temporal stress with those exhibiting stress-related performance decline.

Design, Setting, and Participants

Cohort study conducted from July 2015 to September 2016 at the Imperial College Healthcare National Health Service Trust, England. One hundred two surgical residents (postgraduate year 1 and greater) were invited to participate, of which 33 agreed to partake.

Exposures

Participants performed a laparoscopic suturing task under 2 conditions: self-paced (SP; without time-per-knot restrictions), and time pressure (TP; 2-minute per knot time restriction).

Main Outcomes and Measures

A composite deterioration score was computed based on between-condition differences in task performance metrics (task progression score [arbitrary units], error score [millimeters], leak volume [milliliters], and knot tensile strength [newtons]). Based on the composite score, quartiles were computed reflecting performance stability (quartile 1 [Q1]) and decline (quartile 4 [Q4]). Changes in PFC oxygenated hemoglobin concentration (HbO2) measured at 24 different locations using functional near-infrared spectroscopy were compared between Q1 and Q4. Secondary outcomes included subjective workload (Surgical Task Load Index) and heart rate.

Results

Of the 33 participants, the median age was 33 years, the range was 29 to 56 years, and 27 were men (82%). The Q1 residents demonstrated task-induced increases in HbO2 across the bilateral ventrolateral PFC (VLPFC) and right dorsolateral PFC in the SP condition and in the VLPFC in the TP condition. In contrast, Q4 residents demonstrated decreases in HbO2 in both conditions. The magnitude of PFC activation (change in HbO2) was significantly greater in Q1 than Q4 across the bilateral VLPFC during both SP (mean [SD] left VLPFC: Q1, 0.44 [1.30] μM; Q4, -0.21 [2.05] μM; P < .001; right VLPFC: Q1, 0.46 [1.12] μM; Q4, -0.15 [2.14] μM; P < .001) and TP (mean [SD] left VLPFC: Q1, 0.44 [1.36] μM; Q4, -0.03 [1.83] μM; P = .001; right VLPFC: Q1, 0.49 [1.70] μM; Q4, -0.32 [2.00] μM; P < .001) conditions. There were no significant between-group differences in Surgical Task Load Index or heart rate in either condition.

Conclusions and Relevance

Performance stability within TP is associated with sustained prefrontal activation indicative of preserved attention and concentration, whereas performance decline is associated with prefrontal deactivation that may represent task disengagement.

More homogeneous capillary flow and oxygenation in deeper cortical layers correlate with increased oxygen extraction

Our understanding of how capillary blood flow and oxygen distribute across cortical layers to meet the local metabolic demand is incomplete. We addressed this question by using two-photon imaging of resting-state microvascular oxygen partial pressure (PO₂) and flow in the whisker barrel cortex in awake mice. Our measurements in layers I-V show that the capillary red-blood-cell flux and oxygenation heterogeneity, and the intracapillary resistance to oxygen delivery, all decrease with depth, reaching a minimum around layer IV, while the depth-dependent oxygen extraction fraction is increased in layer IV, where oxygen demand is presumably the highest. Our findings suggest that more homogeneous distribution of the physiological observables relevant to oxygen transport to tissue is an important part of the microvascular network adaptation to local brain metabolism. These results will inform the biophysical models of layer-specific cerebral oxygen delivery and consumption and improve our understanding of the diseases that affect cerebral microcirculation.

Motion artifact detection and correction in functional near-infrared spectroscopy: a new hybrid method based on spline interpolation method and Savitzky-Golay filtering

Motion artifact contamination in near-infrared spectroscopy (NIRS) data has become an important challenge in realizing the full potential of NIRS for real-life applications. Various motion correction algorithms have been used to alleviate the effect of motion artifacts on the estimation of the hemodynamic response function. While smoothing methods, such as wavelet filtering, are excellent in removing motion-induced sharp spikes, the baseline shifts in the signal remain after this type of filtering. Methods, such as spline interpolation, on the other hand, can properly correct baseline shifts; however, they leave residual high-frequency spikes. We propose a hybrid method that takes advantage of different correction algorithms. This method first identifies the baseline shifts and corrects them using a spline interpolation method or targeted principal component analysis. The remaining spikes, on the other hand, are corrected by smoothing methods: Savitzky-Golay (SG) filtering or robust locally weighted regression and smoothing. We have compared our new approach with the existing correction algorithms in terms of hemodynamic response function estimation using the following metrics: mean-squared error, peak-to-peak error ([Formula: see text]), Pearson's correlation ([Formula: see text]), and the area under the receiver operator characteristic curve. We found that spline-SG hybrid method provides reasonable improvements in all these metrics with a relatively short computational time. The dataset and the code used in this study are made available online for the use of all interested researchers.

Optical mapping of the dominant frequency of brain signal oscillations in motor systems

Recent neuroimaging studies revealed that the dominant frequency of neural oscillations is brain-region-specific and can vary with frequency-specific reorganization of brain networks during cognition. In this study, we examined the dominant frequency in low-frequency neural oscillations represented by oxygenated hemoglobin measurements after the hemodynamic response function (HRF) deconvolution. Twenty-nine healthy college subjects were recruited to perform a serial finger tapping task at the frequency of 0.2 Hz. Functional near-infrared spectroscopy (fNIRS) was applied to record the hemodynamic signals over the primary motor cortex, supplementary motor area (SMA), premotor cortex, and prefrontal area. We then explored the low frequency steady-state brain response (lfSSBR), which was evoked in the motor systems at the fundamental frequency (0.2 Hz) and its harmonics (0.4, 0.6, and 0.8 Hz). In particular, after HRF deconvolution, the lfSSBR at the frequency of 0.4 Hz in the SMA was identified as the dominant frequency. Interestingly, the domain frequency exhibited the correlation with behavior data such as reaction time, indicating that the physiological implication of lfSSBR is related to the brain anatomy, stimulus frequency and cognition. More importantly, the HRF deconvolution showed its capability for recovering signals probably reflecting neural-level events and revealing the physiological meaning of lfSSBR.

fNIRS can robustly measure brain activity during memory encoding and retrieval in healthy subjects

Early intervention in Alzheimer’s Disease (AD) requires novel biomarkers that can capture changes in brain activity at an early stage. Current AD biomarkers are expensive and/or invasive and therefore unsuitable for use as screening tools, but a non-invasive, inexpensive, easily accessible screening method could be useful in both clinical and research settings. Prior studies suggest that especially paired-associate learning tasks may be useful in detecting the earliest memory impairment in AD. Here, we investigated the utility of functional Near Infrared Spectroscopy in measuring brain activity from prefrontal, parietal and temporal cortices of healthy adults (n = 19) during memory encoding and retrieval under a face-name paired-associate learning task. Our findings demonstrate that encoding of novel face-name pairs compared to baseline as well as compared to repeated face-name pairs resulted in significant activation in left dorsolateral prefrontal cortex while recalling resulted in activation in dorsolateral prefrontal cortex bilaterally. Moreover, brain response to recalling was significantly higher than encoding in medial, superior and middle frontal cortices for novel faces. Overall, this study shows that fNIRS can reliably measure cortical brain activation during a face-name paired-associate learning task. Future work will include similar measurements in populations with progressing memory deficits.

Characterization and correction of the false-discovery rates in resting state connectivity using functional near-infrared spectroscopy

Functional near-infrared spectroscopy (fNIRS) is a noninvasive neuroimaging technique that uses low levels of red to near-infrared light to measure changes in cerebral blood oxygenation. Spontaneous (resting state) functional connectivity (sFC) has become a critical tool for cognitive neuroscience for understanding task-independent neural networks, revealing pertinent details differentiating healthy from disordered brain function, and discovering fluctuations in the synchronization of interacting individuals during hyperscanning paradigms. Two of the main challenges to sFC-NIRS analysis are (i) the slow temporal structure of both systemic physiology and the response of blood vessels, which introduces false spurious correlations, and (ii) motion-related artifacts that result from movement of the fNIRS sensors on the participants’ head and can introduce non-normal and heavy-tailed noise structures. In this work, we systematically examine the false-discovery rates of several time- and frequency-domain metrics of functional connectivity for characterizing sFC-NIRS. Specifically, we detail the modifications to the statistical models of these methods needed to avoid high levels of false-discovery related to these two sources of noise in fNIRS. We compare these analysis procedures using both simulated and experimental resting-state fNIRS data. Our proposed robust correlation method has better performance in terms of being more reliable to the noise outliers due to the motion artifacts.

Effect of confounding variables on hemodynamic response function estimation using averaging and deconvolution analysis: An event-related NIRS study

Slow and rapid event-related designs are used in fMRI and functional near-infrared spectroscopy (fNIRS) experiments to temporally characterize the brain hemodynamic response to discrete events. Conventional averaging (CA) and the deconvolution method (DM) are the two techniques commonly used to estimate the Hemodynamic Response Function (HRF) profile in event-related designs. In this study, we conducted a series of simulations using synthetic and real NIRS data to examine the effect of the main confounding factors, including event sequence timing parameters, different types of noise, signal-to-noise ratio (SNR), temporal autocorrelation and temporal filtering on the performance of these techniques in slow and rapid event-related designs. We also compared systematic errors in the estimates of the fitted HRF amplitude, latency and duration for both techniques. We further compared the performance of deconvolution methods based on Finite Impulse Response (FIR) basis functions and gamma basis sets. Our results demonstrate that DM was much less sensitive to confounding factors than CA. Event timing was the main parameter largely affecting the accuracy of CA. In slow event-related designs, deconvolution methods provided similar results to those obtained by CA. In rapid event-related designs, our results showed that DM outperformed CA for all SNR, especially above -5 dB regardless of the event sequence timing and the dynamics of background NIRS activity. Our results also show that periodic low-frequency systemic hemodynamic fluctuations as well as phase-locked noise can markedly obscure hemodynamic evoked responses. Temporal autocorrelation also affected the performance of both techniques by inducing distortions in the time profile of the estimated hemodynamic response with inflated t-statistics, especially at low SNRs. We also found that high-pass temporal filtering could substantially affect the performance of both techniques by removing the low-frequency components of HRF profiles. Our results emphasize the importance of characterization of event timing, background noise and SNR when estimating HRF profiles using CA and DM in event-related designs.

Image reconstruction of oxidized cerebral cytochrome C oxidase changes from broadband near-infrared spectroscopy data

In diffuse optical tomography (DOT), overlapping and multidistance measurements are required to reconstruct depth-resolved images of oxy- ([Formula: see text]) and deoxy- (HHb) hemoglobin concentration changes occurring in the brain. These can be considered an indirect measure of brain activity, under the assumption of intact neurovascular coupling. Broadband systems also allow changes in the redox state of cytochrome c oxidase (oxCCO) to be measured, which can be an important biomarker when neurovascular coupling is impaired. We used DOT to reconstruct images of [Formula: see text], [Formula: see text], and [Formula: see text] from data acquired with a broadband system. Four healthy volunteers were measured while performing a visual stimulation task (4-Hz inverting checkerboard). The broadband system was configured to allow multidistance and overlapping measurements of the participants' visual cortex with 32 channels. A multispectral approach was employed to reconstruct changes in concentration of the three chromophores during the visual stimulation. A clear and focused activation was reconstructed in the left occipital cortex of all participants. The difference between the residuals of the three-chromophore model and of the two-chromophore model (recovering only [Formula: see text] and [Formula: see text]) exhibits a spectrum similar to that of oxCCO. These results form a basis for further studies aimed to further optimize image reconstruction of [Formula: see text].

Design of multichannel functional near-infrared spectroscopy system with application to propofol and sevoflurane anesthesia monitoring

Monitoring the changes of cerebral hemodynamics and the state of consciousness during general anesthesia (GA) is clinically important. There is a great need for developing advanced detectors to investigate the physiological processes of the brain during GA. We developed a multichanneled, functional near-infrared spectroscopy (fNIRS) system device and applied it to GA operation monitoring. The cerebral hemodynamic data from the forehead of 11 patients undergoing propofol and sevoflurane anesthesia were analyzed. The concentration changes of oxygenated hemoglobin, deoxygenated hemoglobin, total hemoglobin, and cerebral tissue heart rate were determined from the raw optical information based on the discrete stationary wavelet transform. This custom-made device provides an easy-to-build solution for continuous wave-fNIRS system, with customized specifications. The developed device has a potential value in cerebral monitoring in clinical settings.

Modelling confounding effects from extracerebral contamination and systemic factors on functional near-infrared spectroscopy

Haemodynamics-based neuroimaging is widely used to study brain function. Regional blood flow changes characteristic of neurovascular coupling provide an important marker of neuronal activation. However, changes in systemic physiological parameters such as blood pressure and concentration of CO2 can also affect regional blood flow and may confound haemodynamics-based neuroimaging. Measurements with functional near-infrared spectroscopy (fNIRS) may additionally be confounded by blood flow and oxygenation changes in extracerebral tissue layers. Here we investigate these confounds using an extended version of an existing computational model of cerebral physiology, ‘BrainSignals’. Our results show that confounding from systemic physiological factors is able to produce misleading haemodynamic responses in both positive and negative directions. By applying the model to data from previous fNIRS studies, we demonstrate that such potentially deceptive responses can indeed occur in at least some experimental scenarios. It is therefore important to record the major potential confounders in the course of fNIRS experiments. Our model may then allow the observed behaviour to be attributed among the potential causes and hence reduce identification errors.

•

Confounding of fNIRS haemoglobin signals is simulated using a computational model.

•

Model is extended to simulate scalp haemodynamics.

•

Changes in blood pressure and CO₂ can mimic and mask functional activation.

•

Experimental recording of systemic factors is recommended to aid interpretation.

Mayer waves reduce the accuracy of estimated hemodynamic response functions in functional near-infrared spectroscopy

Analysis of cerebral hemodynamics reveals a wide spectrum of oscillations ranging from 0.0095 to 2 Hz. While most of these oscillations can be filtered out during analysis of functional near-infrared spectroscopy (fNIRS) signals when estimating stimulus evoked hemodynamic responses, oscillations around 0.1 Hz are an exception. This is due to the fact that they share a common spectral range with typical stimulus evoked hemodynamic responses from the brain. Here we investigate the effect of hemodynamic oscillations around 0.1 Hz on the estimation of hemodynamic response functions from fNIRS data. Our results show that for an expected response of ~1 µM in oxygenated hemoglobin concentration (HbO), Mayer wave oscillations with an amplitude > ~1 µM at 0.1 Hz reduce the accuracy of the estimated response as quantified by a 3 fold increase in the mean squared error and decrease in correlation (R(2) below 0.78) when compared to the true HRF. These results indicate that the amplitude of oscillations at 0.1 Hz can serve as an objective metric of the expected HRF estimation accuracy. In addition, we investigated the effect of short separation regression on the recovered HRF, and found that this improves the recovered HRF when large amplitude 0.1 Hz oscillations are present in fNIRS data. We suspect that the development of other filtering strategies may provide even further improvement.

Dynamic filtering improves attentional state prediction with fNIRS

Brain activity can predict a person's level of engagement in an attentional task. However, estimates of brain activity are often confounded by measurement artifacts and systemic physiological noise. The optimal method for filtering this noise - thereby increasing such state prediction accuracy - remains unclear. To investigate this, we asked study participants to perform an attentional task while we monitored their brain activity with functional near infrared spectroscopy (fNIRS). We observed higher state prediction accuracy when noise in the fNIRS hemoglobin [Hb] signals was filtered with a non-stationary (adaptive) model as compared to static regression (84% ± 6% versus 72% ± 15%).

Commentary on the statistical properties of noise and its implication on general linear models in functional near-infrared spectroscopy

Functional near-infrared spectroscopy (fNIRS) is a noninvasive neuroimaging technique that uses low levels of light to measure changes in cerebral blood oxygenation levels. In the majority of NIRS functional brain studies, analysis of this data is based on a statistical comparison of hemodynamic levels between a baseline and task or between multiple task conditions by means of a linear regression model: the so-called general linear model. Although these methods are similar to their implementation in other fields, particularly for functional magnetic resonance imaging, the specific application of these methods in fNIRS research differs in several key ways related to the sources of noise and artifacts unique to fNIRS. In this brief communication, we discuss the application of linear regression models in fNIRS and the modifications needed to generalize these models in order to deal with structured (colored) noise due to systemic physiology and noise heteroscedasticity due to motion artifacts. The objective of this work is to present an overview of these noise properties in the context of the linear model as it applies to fNIRS data. This work is aimed at explaining these mathematical issues to the general fNIRS experimental researcher but is not intended to be a complete mathematical treatment of these concepts.

Pressure modulation algorithm to separate cerebral hemodynamic signals from extracerebral artifacts

We introduce and validate a pressure measurement paradigm that reduces extracerebral contamination from superficial tissues in optical monitoring of cerebral blood flow with diffuse correlation spectroscopy (DCS). The scheme determines subject-specific contributions of extracerebral and cerebral tissues to the DCS signal by utilizing probe pressure modulation to induce variations in extracerebral blood flow. For analysis, the head is modeled as a two-layer medium and is probed with long and short source-detector separations. Then a combination of pressure modulation and a modified Beer-Lambert law for flow enables experimenters to linearly relate differential DCS signals to cerebral and extracerebral blood flow variation without a priori anatomical information. We demonstrate the algorithm's ability to isolate cerebral blood flow during a finger-tapping task and during graded scalp ischemia in healthy adults. Finally, we adapt the pressure modulation algorithm to ameliorate extracerebral contamination in monitoring of cerebral blood oxygenation and blood volume by near-infrared spectroscopy.

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.

Correspondence of electroencephalography and near-infrared spectroscopy sensitivities to the cerebral cortex using a high-density layout

This study investigates the correspondence of the cortical sensitivity of electroencephalography (EEG) and near-infrared spectroscopy (NIRS). EEG forward model sensitivity to the cerebral cortex was calculated for 329 EEG electrodes following the 10-5 EEG positioning system using a segmented structural magnetic resonance imaging scan of a human subject. NIRS forward model sensitivity was calculated for the same subject using 156 NIRS source-detector pairs selected from 32 source and 32 detector optodes positioned on the scalp using a subset of the 10-5 EEG positioning system. Sensitivity correlations between colocalized NIRS source-detector pair groups and EEG channels yielded R = 0.46 ± 0.08. Groups of NIRS source-detector pairs with maximum correlations to EEG electrode sensitivities are tabulated. The mean correlation between the point spread functions for EEG and NIRS regions of interest (ROI) was R = 0.43 ± 0.07. Spherical ROIs with radii of 26 mm yielded the maximum correlation between EEG and NIRS averaged across all cortical mesh nodes. These sensitivity correlations between EEG and NIRS should be taken into account when designing multimodal studies of neurovascular coupling and when using NIRS as a statistical prior for EEG source localization.

Further improvement in reducing superficial contamination in NIRS using double short separation measurements

Near-Infrared Spectroscopy (NIRS) allows the recovery of the evoked hemodynamic response to brain activation. In adult human populations, the NIRS signal is strongly contaminated by systemic interference occurring in the superficial layers of the head. An approach to overcome this difficulty is to use additional NIRS measurements with short optode separations to measure the systemic hemodynamic fluctuations occurring in the superficial layers. These measurements can then be used as regressors in the post-experiment analysis to remove the systemic contamination and isolate the brain signal. In our previous work, we showed that the systemic interference measured in NIRS is heterogeneous across the surface of the scalp. As a consequence, the short separation measurement used in the regression procedure must be located close to the standard NIRS channel from which the evoked hemodynamic response of the brain is to be recovered. Here, we demonstrate that using two short separation measurements, one at the source optode and one at the detector optode, further increases the performance of the short separation regression method compared to using a single short separation measurement. While a single short separation channel produces an average reduction in noise of 33% for HbO, using a short separation channel at both source and detector reduces noise by 59% compared to the standard method using a general linear model (GLM) without short separation. For HbR, noise reduction of 3% is achieved using a single short separation and this number goes to 47% when two short separations are used. Our work emphasizes the importance of integrating short separation measurements both at the source and at the detector optode of the standard channels from which the hemodynamic response is to be recovered. While the implementation of short separation sources presents some difficulties experimentally, the improvement in noise reduction is significant enough to justify the practical challenges.

Further improvement in reducing superficial contamination in NIRS using double short separation measurements

Near-Infrared Spectroscopy (NIRS) allows the recovery of the evoked hemodynamic response to brain activation. In adult human populations, the NIRS signal is strongly contaminated by systemic interference occurring in the superficial layers of the head. An approach to overcome this difficulty is to use additional NIRS measurements with short optode separations to measure the systemic hemodynamic fluctuations occurring in the superficial layers. These measurements can then be used as regressors in the post-experiment analysis to remove the systemic contamination and isolate the brain signal. In our previous work, we showed that the systemic interference measured in NIRS is heterogeneous across the surface of the scalp. As a consequence, the short separation measurement used in the regression procedure must be located close to the standard NIRS channel from which the evoked hemodynamic response of the brain is to be recovered. Here, we demonstrate that using two short separation measurements, one at the source optode and one at the detector optode, further increases the performance of the short separation regression method compared to using a single short separation measurement. While a single short separation channel produces an average reduction in noise of 33% for HbO, using a short separation channel at both source and detector reduces noise by 59% compared to the standard method using a general linear model (GLM) without short separation. For HbR, noise reduction of 3% is achieved using a single short separation and this number goes to 47% when two short separations are used. Our work emphasizes the importance of integrating short separation measurements both at the source and at the detector optode of the standard channels from which the hemodynamic response is to be recovered. While the implementation of short separation sources presents some difficulties experimentally, the improvement in noise reduction is significant enough to justify the practical challenges.

Autoregressive model based algorithm for correcting motion and serially correlated errors in fNIRS

Systemic physiology and motion-induced artifacts represent two major sources of confounding noise in functional near infrared spectroscopy (fNIRS) imaging that can reduce the performance of analyses and inflate false positive rates (i.e., type I errors) of detecting evoked hemodynamic responses. In this work, we demonstrated a general algorithm for solving the general linear model (GLM) for both deconvolution (finite impulse response) and canonical regression models based on designing optimal pre-whitening filters using autoregressive models and employing iteratively reweighted least squares. We evaluated the performance of the new method by performing receiver operating characteristic (ROC) analyses using synthetic data, in which serial correlations, motion artifacts, and evoked responses were controlled via simulations, as well as using experimental data from children (3-5 years old) as a source baseline physiological noise and motion artifacts. The new method outperformed ordinary least squares (OLS) with no motion correction, wavelet based motion correction, or spline interpolation based motion correction in the presence of physiological and motion related noise. In the experimental data, false positive rates were as high as 37% when the estimated p-value was 0.05 for the OLS methods. The false positive rate was reduced to 5-9% with the proposed method. Overall, the method improves control of type I errors and increases performance when motion artifacts are present.

Patient-oriented simulation based on Monte Carlo algorithm by using MRI data

Background

Although Monte Carlo simulations of light propagation in full segmented three-dimensional MRI based anatomical models of the human head have been reported in many articles. To our knowledge, there is no patient-oriented simulation for individualized calibration with NIRS measurement. Thus, we offer an approach for brain modeling based on image segmentation process with in vivo MRI T1 three-dimensional image to investigate the individualized calibration for NIRS measurement with Monte Carlo simulation.

Methods

In this study, an individualized brain is modeled based on in vivo MRI 3D image as five layers structure. The behavior of photon migration was studied for this individualized brain detections based on three-dimensional time-resolved Monte Carlo algorithm. During the Monte Carlo iteration, all photon paths were traced with various source-detector separations for characterization of brain structure to provide helpful information for individualized design of NIRS system.

Results

Our results indicate that the patient-oriented simulation can provide significant characteristics on the optimal choice of source-detector separation within 3.3 cm of individualized design in this case. Significant distortions were observed around the cerebral cortex folding. The spatial sensitivity profile penetrated deeper to the brain in the case of expanded CSF. This finding suggests that the optical method may provide not only functional signal from brain activation but also structural information of brain atrophy with the expanded CSF layer. The proposed modeling method also provides multi-wavelength for NIRS simulation to approach the practical NIRS measurement.

Conclusions

In this study, the three-dimensional time-resolved brain modeling method approaches the realistic human brain that provides useful information for NIRS systematic design and calibration for individualized case with prior MRI data.

Implicit and explicit prior information in near-infrared spectral imaging: accuracy, quantification and diagnostic value

Near-infrared spectroscopy (NIRS) of tissue provides quantification of absorbers, scattering and luminescent agents in bulk tissue through the use of measurement data and assumptions. Prior knowledge can be critical about things such as (i) the tissue shape and/or structure, (ii) spectral constituents, (iii) limits on parameters, (iv) demographic or biomarker data, and (v) biophysical models of the temporal signal shapes. A general framework of NIRS imaging with prior information is presented, showing that prior information datasets could be incorporated at any step in the NIRS process, with the general workflow being: (i) data acquisition, (ii) pre-processing, (iii) forward model, (iv) inversion/reconstruction, (v) post-processing, and (vi) interpretation/diagnosis. Most of the development in NIRS has used ad hoc or empirical implementations of prior information such as pre-measured absorber or fluorophore spectra, or tissue shapes as estimated by additional imaging tools. A comprehensive analysis would examine what prior information maximizes the accuracy in recovery and value for medical diagnosis, when implemented at separate stages of the NIRS sequence. Individual applications of prior information can show increases in accuracy or improved ability to estimate biochemical features of tissue, while other approaches may not. Most beneficial inclusion of prior information has been in the inversion/reconstruction process, because it solves the mathematical intractability. However, it is not clear that this is always the most beneficial stage.

The utility of near-infrared spectroscopy in the regression of low-frequency physiological noise from functional magnetic resonance imaging data

Near-infrared spectroscopy (NIRS) signals have been shown to correlate with resting-state BOLD-fMRI data across the whole brain volume, particularly at frequencies below 0.1Hz. While the physiological origins of this correlation remain unclear, its existence may have a practical application in minimizing the background physiological noise present in BOLD-fMRI recordings. We performed simultaneous, resting-state fMRI and 28-channel NIRS in seven adult subjects in order to assess the utility of NIRS signals in the regression of physiological noise from fMRI data. We calculated the variance of the residual error in a general linear model of the baseline fMRI signal, and the reduction of this variance achieved by including NIRS signals in the model. In addition, we introduced a sequence of simulated hemodynamic response functions (HRFs) into the resting-state fMRI data of each subject in order to quantify the effectiveness of NIRS signals in optimizing the recovery of that HRF. For comparison, these calculations were also performed using a pulse and respiration RETROICOR model. Our results show that the use of 10 or more NIRS channels can reduce variance in the residual error by as much as 36% on average across the whole cortex. However the same number of low-pass filtered white noise regressors is shown to produce a reduction of 19%. The RETROICOR model obtained a variance reduction of 6.4%. Our HRF simulation showed that the mean-squared error (MSE) between the recovered and true HRFs is reduced by 21% on average when 10 NIRS channels are applied and by introducing an optimized time lag between the NIRS and fMRI time series, a single NIRS channel can provide an average MSE reduction of 14%. The RETROICOR model did not provide a significant change in MSE. By each of the metrics calculated, NIRS recording is shown to be of significant benefit to the regression of low-frequency physiological noise from fMRI data.

Diffuse optical multipatch technique for tissue oxygenation monitoring: clinical study in intensive care unit

Diffuse optical multipatch technique is used to assess spatial variations in absorption and scattering in biological tissue, by monitoring changes in the concentration of oxyhemoglobin and deoxyhemoglobin. In our preliminary study, the temporal tracings of tissue oxygenation are measured using diffuse optical multipatch measurement and a venous occlusion test, employing normal subjects and ICU patients suffering from sepsis and heart failure. In experiments, obvious differences in tissue oxygenation signals were observed among all three groups. This paper discusses the physiological relevance of tissue oxygenation with respect to disease.

Improved recovery of the hemodynamic response in diffuse optical imaging using short optode separations and state-space modeling

Diffuse optical imaging (DOI) allows the recovery of the hemodynamic response associated with evoked brain activity. The signal is contaminated with systemic physiological interference which occurs in the superficial layers of the head as well as in the brain tissue. The back-reflection geometry of the measurement makes the DOI signal strongly contaminated by systemic interference occurring in the superficial layers. A recent development has been the use of signals from small source-detector separation (1cm) optodes as regressors. Since those additional measurements are mainly sensitive to superficial layers in adult humans, they help in removing the systemic interference present in longer separation measurements (3 cm). Encouraged by those findings, we developed a dynamic estimation procedure to remove global interference using small optode separations and to estimate simultaneously the hemodynamic response. The algorithm was tested by recovering a simulated synthetic hemodynamic response added over baseline DOI data acquired from 6 human subjects at rest. The performance of the algorithm was quantified by the Pearson R(2) coefficient and the mean square error (MSE) between the recovered and the simulated hemodynamic responses. Our dynamic estimator was also compared with a static estimator and the traditional adaptive filtering method. We observed a significant improvement (two-tailed paired t-test, p<0.05) in both HbO and HbR recovery using our Kalman filter dynamic estimator compared to the traditional adaptive filter, the static estimator and the standard GLM technique.

Kalman estimator- and general linear model-based on-line brain activation mapping by near-infrared spectroscopy

BACKGROUND

Near-infrared spectroscopy (NIRS) is a non-invasive neuroimaging technique that recently has been developed to measure the changes of cerebral blood oxygenation associated with brain activities. To date, for functional brain mapping applications, there is no standard on-line method for analysing NIRS data.

METHODS

In this paper, a novel on-line NIRS data analysis framework taking advantages of both the general linear model (GLM) and the Kalman estimator is devised. The Kalman estimator is used to update the GLM coefficients recursively, and one critical coefficient regarding brain activities is then passed to a t-statistical test. The t-statistical test result is used to update a topographic brain activation map. Meanwhile, a set of high-pass filters is plugged into the GLM to prevent very low-frequency noises, and an autoregressive (AR) model is used to prevent the temporal correlation caused by physiological noises in NIRS time series. A set of data recorded in finger tapping experiments is studied using the proposed framework.

RESULTS

The obtained results suggest that the method can effectively track the task related brain activation areas, and prevent the noise distortion in the estimation while the experiment is running. Thereby, the potential of the proposed method for real-time NIRS-based brain imaging was demonstrated.

CONCLUSIONS

This paper presents a novel on-line approach for analysing NIRS data for functional brain mapping applications. This approach demonstrates the potential of a real-time-updating topographic brain activation map.

Bayesian filtering of human brain hemodynamic activity elicited by visual short-term maintenance recorded through functional near-infrared spectroscopy (fNIRS)

Functional near-infrared spectroscopy (fNIRS) is a neuroimaging technique that measures changes in oxy-hemoglobin (ΔHbO) and deoxy-hemoglobin (ΔHbR) concentration associated with brain activity. The signal acquired with fNIRS is naturally affected by disturbances engendering from ongoing physiological activity (e.g., cardiac, respiratory, Mayer wave) and random measurement noise. Despite its several drawbacks, the so-called conventional averaging (CA) is still widely used to estimate the hemodynamic response function (HRF) from noisy signal. One such drawback is related to the number of trials necessary to derive stable HRF functions adopting the CA approach, which must be substantial (N >> 50). In this work, a pre-processing procedure to remove artifacts followed by the application of a non-parametric Bayesian approach is proposed that capitalizes on a priori available knowledge about HRF and noise. Results with the proposed Bayesian approach were compared with CA and with a straightforward band-pass filtering approach. On simulated data, a five times lower estimation error on HRF was obtained with respect to that obtained by CA, and 2.5 times lower than that obtained by band pass filtering. On real data, the improvement achieved by the present method was attested by an increase in the contrast to noise ratio (CNR) and by a reduced variability in single trial estimation. An application of the present Bayesian approach is illustrated that was optimized to monitor changes in hemodynamic activity reflecting variations in visual short-term memory load in humans, which are notoriously hard to detect using functional magnetic resonance imaging (fMRI). In particular, statistical analyses of HRFs recorded during a memory task established with high reliability the crucial role of the intraparietal sulcus and the intra-occipital sulcus in posterior areas of the human brain in visual short-term memory maintenance.

Resting state functional connectivity of the whole head with near-infrared spectroscopy

Resting state connectivity aims to identify spontaneous cerebral hemodynamic fluctuations that reflect neuronal activity at rest. In this study, we investigated the spatial-temporal correlation of hemoglobin concentration signals over the whole head during the resting state. By choosing a source-detector pair as a seed, we calculated the correlation value between its time course and the time course of all other source-detector combinations, and projected them onto a topographic map. In all subjects, we found robust spatial interactions in agreement with previous fMRI and NIRS findings. Strong correlations between the two opposite hemispheres were seen for both sensorimotor and visual cortices. Correlations in the prefrontal cortex were more heterogeneous and dependent on the hemodynamic contrast. HbT provided robust, well defined maps, suggesting that this contrast may be used to better localize functional connectivity. The effects of global systemic physiology were also investigated, particularly low frequency blood pressure oscillations which give rise to broad regions of high correlation and mislead interpretation of the results. These results confirm the feasibility of using functional connectivity with optical methods during the resting state, and validate its use to investigate cortical interactions across the whole head.

Diffuse Optics for Tissue Monitoring and Tomography

This review describes the diffusion model for light transport in tissues and the medical applications of diffuse light. Diffuse optics is particularly useful for measurement of tissue hemodynamics, wherein quantitative assessment of oxy- and deoxy-hemoglobin concentrations and blood flow are desired. The theoretical basis for near-infrared or diffuse optical spectroscopy (NIRS or DOS, respectively) is developed, and the basic elements of diffuse optical tomography (DOT) are outlined. We also discuss diffuse correlation spectroscopy (DCS), a technique whereby temporal correlation functions of diffusing light are transported through tissue and are used to measure blood flow. Essential instrumentation is described, and representative brain and breast functional imaging and monitoring results illustrate the workings of these new tissue diagnostics.

Topographic localization of brain activation in diffuse optical imaging using spherical wavelets

Diffuse optical imaging is a non-invasive technique that uses near-infrared light to measure changes in brain activity through an array of sensors placed on the surface of the head. Compared to functional MRI, optical imaging has the advantage of being portable while offering the ability to record functional changes in both oxy- and deoxy-hemoglobin within the brain at a high temporal resolution. However, the reconstruction of accurate spatial images of brain activity from optical measurements represents an ill-posed and underdetermined problem that requires regularization. These reconstructions benefit from incorporating prior information about the underlying spatial structure and function of the brain. In this work, we describe a novel image reconstruction approach which uses surface-based wavelets derived from structural MRI to incorporate high-resolution anatomical and structural prior information about the brain. This surface-based approach is used to approximate brain activation patterns through the reconstruction and presentation of topographical (two-dimensional) maps of brain activation directly onto the folded surface of the cortex. The set of wavelet coefficients is directly estimated by a truncated singular-value decomposition based pseudo-inversion of the wavelet projection of the optical forward model. We use a reconstruction metric based on Shannon entropy which quantifies the sparse loading of the wavelet coefficients and is used to determine the optimal truncation and regularization of this inverse model. In this work, examples of the performance of this model are illustrated for several cases of numerical simulation and experimental data with comparison to functional magnetic resonance imaging.

Comparison of principal and independent component analysis in removing extracerebral interference from near-infrared spectroscopy signals

Near-infrared spectroscopy (NIRS) is a method for noninvasive estimation of cerebral hemodynamic changes. Principal component analysis (PCA) and independent component analysis (ICA) can be used for decomposing a set of signals to underlying components. Our objective is to determine whether PCA or ICA is more efficient in identifying and removing scalp blood flow interference from multichannel NIRS signals. Concentration changes of oxygenated (HbO(2)) and deoxygenated (HbR) hemoglobin are measured on the forehead with multichannel NIRS during hyper- and hypocapnia. PCA and ICA are used separately to identify and remove signal contribution from extracerebral tissue, and the resulting estimates of cerebral responses are compared to the expected cerebral responses. Both methods were able to reduce extracerebral contribution to the signals, but PCA typically performs equal to or better than ICA. The improvement in 3-cm signal quality achieved with both methods is comparable to increasing the source-detector separation from 3 to 5 cm. Especially PCA appears to be well suited for use in NIRS applications where the cerebral activation is diffuse, such as monitoring of global cerebral oxygenation and hemodynamics. Performance differences between PCA and ICA could be attributed primarily to different criteria for identifying the surface effect.