Multi-subject and multi-task experimental validation of the hierarchical Bayesian diffuse optical tomography algorithm

Surface-based integration approach for fNIRS-fMRI reliability assessment

INTRODUCTION

Studies integrating functional near-infrared spectroscopy (fNIRS) with functional MRI (fMRI) employ heterogeneous methods in defining common regions of interest in which similarities are assessed. Therefore, spatial agreement and temporal correlation may not be reproducible across studies. In the present work, we address this issue by proposing a novel method for integration and analysis of fNIRS and fMRI over the cortical surface.

MATERIALS AND METHODS

Eighteen healthy volunteers (age mean±SD 30.55 ± 4.7, 7 males) performed a motor task during non-simultaneous fMRI and fNIRS acquisitions. First, fNIRS and fMRI data were integrated by projecting subject- and group-level source maps over the cortical surface mesh to define anatomically constrained functional ROIs (acfROI). Next, spatial agreement and temporal correlation were quantified as Dice Coefficient (DC) and Pearson's correlation coefficient between fNIRS-fMRI in the acfROIs.

RESULTS

Subject-level results revealed moderate to substantial spatial agreement (DC range 0.43 - 0.64), confirmed at the group-level only for blood oxygenation level-dependent (BOLD) signal vs. HbO2 (0.44 - 0.69), while lack of agreement was found for BOLD vs. HbR in some instances (0.05 - 0.49). Subject-level temporal correlation was moderate to strong (0.79 - 0.85 for BOLD vs. HbO2 and -0.62 to -0.72 for BOLD vs. HbR), while an overall strong correlation was found for group-level results (0.95 - 0.98 for BOLD vs. HbO2 and -0.91 to -0.94 for BOLD vs. HbR).

CONCLUSION

The proposed method directly compares fNIRS and fMRI by projecting individual source maps to the cortical surface. Our results indicate spatial and temporal correspondence between fNIRS and fMRI, and promotes the use of fNIRS when more ecological acquision settings are required, such as longitudinal monitoring of brain activity before and after rehabilitation.

Brain space image reconstruction of functional near-infrared spectroscopy using a Bayesian adaptive fused sparse overlapping group lasso model

Significance

Functional near-infrared spectroscopy (fNIRS) is a noninvasive technology that uses low levels of nonionizing light in the range of red and near-infrared to record changes in the optical absorption and scattering of the underlying tissue that can be used to infer blood flow and oxygen changes during brain activity. The challenges and difficulties of reconstructing spatial images of hemoglobin changes from fNIRS data are mainly caused by the illposed nature of the optical inverse model.

Aim

We describe a Bayesian approach combining several lasso-based regularizations to apply anatomy-prior information to solving the inverse model.

Approach

We built a Bayesian hierarchical model to solve the Bayesian adaptive fused sparse overlapping group lasso (Ba-FSOGL) model. The method is evaluated and validated using simulation and experimental datasets.

Results

We apply this approach to the simulation and experimental datasets to reconstruct a known brain activity. The reconstructed images and statistical plots are shown.

Conclusion

We discuss the adaptation of this method to fNIRS data and demonstrate that this approach provides accurate image reconstruction with a low false-positive rate, through numerical simulations and application to experimental data collected during motor and sensory tasks.

Diffuse optical reconstructions of functional near infrared spectroscopy data using maximum entropy on the mean

Functional near-infrared spectroscopy (fNIRS) measures the hemoglobin concentration changes associated with neuronal activity. Diffuse optical tomography (DOT) consists of reconstructing the optical density changes measured from scalp channels to the oxy-/deoxy-hemoglobin concentration changes within the cortical regions. In the present study, we adapted a nonlinear source localization method developed and validated in the context of Electro- and Magneto-Encephalography (EEG/MEG): the Maximum Entropy on the Mean (MEM), to solve the inverse problem of DOT reconstruction. We first introduced depth weighting strategy within the MEM framework for DOT reconstruction to avoid biasing the reconstruction results of DOT towards superficial regions. We also proposed a new initialization of the MEM model improving the temporal accuracy of the original MEM framework. To evaluate MEM performance and compare with widely used depth weighted Minimum Norm Estimate (MNE) inverse solution, we applied a realistic simulation scheme which contained 4000 simulations generated by 250 different seeds at different locations and 4 spatial extents ranging from 3 to 40[Formula: see text] along the cortical surface. Our results showed that overall MEM provided more accurate DOT reconstructions than MNE. Moreover, we found that MEM was remained particularly robust in low signal-to-noise ratio (SNR) conditions. The proposed method was further illustrated by comparing to functional Magnetic Resonance Imaging (fMRI) activation maps, on real data involving finger tapping tasks with two different montages. The results showed that MEM provided more accurate HbO and HbR reconstructions in spatial agreement with the main fMRI cluster, when compared to MNE.

Diffuse Optical Tomography Using fNIRS Signals Measured from the Skull Surface of the Macaque Monkey

Diffuse optical tomography (DOT), as a functional near-infrared spectroscopy (fNIRS) technique, can estimate three-dimensional (3D) images of the functional hemodynamic response in brain volume from measured optical signals. In this study, we applied DOT algorithms to the fNIRS data recorded from the surface of macaque monkeys’ skulls when the animals performed food retrieval tasks using either the left- or right-hand under head-free conditions. The hemodynamic response images, reconstructed by DOT with a high sampling rate and fine voxel size, demonstrated significant activations at the upper limb regions of the primary motor area in the central sulcus and premotor, and parietal areas contralateral to the hands used in the tasks. The results were also reliable in terms of consistency across different recording dates. Time-series analyses of each brain area revealed preceding activity of premotor area to primary motor area consistent with previous physiological studies. Therefore, the fNIRS–DOT protocol demonstrated in this study provides reliable 3D functional brain images over a period of days under head-free conditions for region-of-interest–based time-series analysis.

Evaluation of a personalized functional near infra-red optical tomography workflow using maximum entropy on the mean

In the present study, we proposed and evaluated a workflow of personalized near infra-red optical tomography (NIROT) using functional near-infrared spectroscopy (fNIRS) for spatiotemporal imaging of cortical hemodynamic fluctuations. The proposed workflow from fNIRS data acquisition to local 3D reconstruction consists of: (a) the personalized optimal montage maximizing fNIRS channel sensitivity to a predefined targeted brain region; (b) the optimized fNIRS data acquisition involving installation of optodes and digitalization of their positions using a neuronavigation system; and (c) the 3D local reconstruction using maximum entropy on the mean (MEM) to accurately estimate the location and spatial extent of fNIRS hemodynamic fluctuations along the cortical surface. The workflow was evaluated on finger-tapping fNIRS data acquired from 10 healthy subjects for whom we estimated the reconstructed NIROT spatiotemporal images and compared with functional magnetic resonance imaging (fMRI) results from the same individuals. Using the fMRI activation maps as our reference, we quantitatively compared the performance of two NIROT approaches, the MEM framework and the conventional minimum norm estimation (MNE) method. Quantitative comparisons were performed at both single subject and group-level. Overall, our results suggested that MEM provided better spatial accuracy than MNE, while both methods offered similar temporal accuracy when reconstructing oxygenated (HbO) and deoxygenated hemoglobin (HbR) concentration changes evoked by finger-tapping. Our proposed complete workflow was made available in the brainstorm fNIRS processing plugin-NIRSTORM, thus providing the opportunity for other researchers to further apply it to other tasks and on larger populations.

Diffuse optical tomography to measure functional changes during motor tasks: a motor imagery study

The present work shows the spatial reliability of the diffuse optical tomography (DOT) system in a group of healthy subjects during a motor imagery task. Prior to imagery task performance, the subjects executed a motor task based on the finger to thumb opposition for motor training, and to corroborate the DOT spatial localization during the motor execution. DOT technology and data treatment allows us to distinguish oxy- and deoxyhemoglobin at the cerebral gyri level unlike the cerebral activations provided by fMRI series that were processed using different approaches. Here we show the DOT reliability showing functional activations at the cerebral gyri level during motor execution and motor imagery, which provide subtler cerebral activations than the motor execution. These results will allow the use of the DOT system as a monitoring device in a brain computer interface.

Functional near-infrared-spectroscopy-based measurement of changes in cortical activity in macaques during post-infarct recovery of manual dexterity

Because compensatory changes in brain activity underlie functional recovery after brain damage, monitoring of these changes will help to improve rehabilitation effectiveness. Functional near-infrared spectroscopy (fNIRS) has the potential to measure brain activity in freely moving subjects. We recently established a macaque model of internal capsule infarcts and an fNIRS system for use in the monkey brain. Here, we used these systems to study motor recovery in two macaques, for which focal infarcts of different sizes were induced in the posterior limb of the internal capsule. Immediately after the injection, flaccid paralysis was observed in the hand contralateral to the injected hemisphere. Thereafter, dexterous hand movements gradually recovered over months. After movement recovery, task-evoked hemodynamic responses increased in the ventral premotor cortex (PMv). The response in the PMv of the infarcted (i.e., ipsilesional) hemisphere increased in the monkey that had received less damage. In contrast, the PMv of the non-infarcted (contralesional) hemisphere was recruited in the monkey with more damage. A pharmacological inactivation experiment with muscimol suggested the involvement of these areas in dexterous hand movements during recovery. These results indicate that fNIRS can be used to evaluate brain activity changes crucial for functional recovery after brain damage.

Resting-State Functional Connectivity Estimated With Hierarchical Bayesian Diffuse Optical Tomography

Resting-state functional connectivity (RSFC) has been generally assessed with functional magnetic resonance imaging (fMRI) thanks to its high spatial resolution. However, fMRI has several disadvantages such as high cost and low portability. In addition, fMRI may not be appropriate for people with metal or electronic implants in their bodies, with claustrophobia and who are pregnant. Diffuse optical tomography (DOT), a method of neuroimaging using functional near-infrared spectroscopy (fNIRS) to reconstruct three-dimensional brain activity images, offers a non-invasive alternative, because fNIRS as well as fMRI measures changes in deoxygenated hemoglobin concentrations and, in addition, fNIRS is free of above disadvantages. We recently proposed a hierarchical Bayesian (HB) DOT algorithm and verified its performance in terms of task-related brain responses. In this study, we attempted to evaluate the HB DOT in terms of estimating RSFC. In 20 healthy males (21-38 years old), 10 min of resting-state data was acquired with 3T MRI scanner or high-density NIRS on different days. The NIRS channels consisted of 96 long (29-mm) source-detector (SD) channels and 56 short (13-mm) SD channels, which covered bilateral frontal and parietal areas. There were one and two resting-state runs in the fMRI and fNIRS experiments, respectively. The reconstruction performances of our algorithm and the two currently prevailing algorithms for DOT were evaluated using fMRI signals as a reference. Compared with the currently prevailing algorithms, our HB algorithm showed better performances in both the similarity to fMRI data and inter-run reproducibility, in terms of estimating the RSFC.

Non-contact acquisition of brain function using a time-extracted compact camera

A revolution in functional brain imaging techniques is in progress in the field of neurosciences. Optical imaging techniques, such as high-density diffuse optical tomography (HD-DOT), in which source-detector pairs of probes are placed on subjects' heads, provide better portability than conventional functional magnetic resonance imaging (fMRI) equipment. However, these techniques remain costly and can only acquire images at up to a few measurements per square centimetre, even when multiple detector probes are employed. In this study, we demonstrate functional brain imaging using a compact and affordable setup that employs nanosecond-order pulsed ordinary laser diodes and a time-extracted image sensor with superimposition capture of scattered components. Our technique can simply and easily attain a high density of measurement points without requiring probes to be attached, and can directly capture two-dimensional functional brain images. We have demonstrated brain activity imaging using a phantom that mimics the optical properties of an adult human head, and with a human subject, have measured cognitive brain activation while the subject is solving simple arithmetical tasks.

Development of multi-directional functional near-infrared spectroscopy system for human neuroimaging studies

Multi-directional measurement using multi-directional light sources and multi-directional photodetectors drastically increases the amount of observation data without increasing the number of optical probes. In this study, we developed a novel multi-directional functional near-infrared spectroscopy (fNIRS) system for human neuroimaging studies. We tested our system by measuring the cortical hemodynamic changes of a single subject during a motor task and compared them with the same subject's functional magnetic resonance imaging (fMRI) data. We detected the direction-dependent fNIRS signals that originate from the cortical hemodynamic changes that are consistent with the fMRI data.

Is it possible to measure hemodynamic changes in the prefrontal cortex through the frontal sinus using continuous wave DOT systems?

The present work shows the capability of near infrared (NIR) light to reach the cerebral cortex through the frontal sinus using continuous-wave techniques (CW-DOT) in a dual study. On the one hand, changes in time during the tracking of a blood dye in the prefrontal cortex were monitored. On the other hand, hemodynamic changes induced by low frequency of transcranial magnetic stimulation applied on the prefrontal cortex were recorded. The results show how NIR light projected through the frontal sinus reaches the cerebral cortex target, providing enough information to have a reliable measurement of cortical hemodynamic changes using CW-DOT.

Comparison of source localization techniques in diffuse optical tomography for fNIRS application using a realistic head model

Functional near-infrared spectroscopy (fNIRS) is a non-invasive imaging technique that elicits growing interest for research and clinical applications. In the last decade, efforts have been made to develop a mathematical framework in order to image the effective sources of hemoglobin variations in brain tissues. Different approaches can be used to impose additional information or constraints when reconstructing the cerebral images of an ill-posed problem. The goal of this study is to compare the performance and limitations of several source localization techniques in the context of fNIRS tomography using individual anatomical magnetic resonance imaging (MRI) to model light propagation. The forward problem is solved using a Monte Carlo simulation of light propagation in the tissues. The inverse problem has been linearized using the Rytov approximation. Then, Tikhonov regularization applied to least squares, truncated singular value decomposition, back-projection, L1-norm regularization, minimum norm estimates, low resolution electromagnetic tomography and Bayesian model averaging techniques are compared using a receiver operating characteristic analysis, blurring and localization error measures. Using realistic simulations (n = 450) and data acquired from a human participant, this study depicts how these source localization techniques behave in a human head fNIRS tomography. When compared to other methods, Bayesian model averaging is proposed as a promising method in DOT and shows great potential to improve specificity, accuracy, as well as to reduce blurring and localization error even in presence of noise and deep sources. Classical reconstruction methods, such as regularized least squares, offer better sensitivity but higher blurring; while more novel L1-based method provides sparse solutions with small blurring and high specificity but lower sensitivity. The application of these methods is also demonstrated experimentally using visual fNIRS experiment with adult participant.

Comparing diffuse optical tomography and functional magnetic resonance imaging signals during a cognitive task: pilot study

Diffuse optical tomography (DOT) measures concentration changes in both oxy- and deoxyhemoglobin providing three-dimensional images of local brain activations. A pilot study, which compares both DOT and functional magnetic resonance imaging (fMRI) volumes through t-maps given by canonical statistical parametric mapping (SPM) processing for both data modalities, is presented. The DOT series were processed using a method that is based on a Bayesian filter application on raw DOT data to remove physiological changes and minimum description length application index to select a number of singular values, which reduce the data dimensionality during image reconstruction and adaptation of DOT volume series to normalized standard space. Therefore, statistical analysis is performed with canonical SPM software in the same way as fMRI analysis is done, accepting DOT volumes as if they were fMRI volumes. The results show the reproducibility and ruggedness of the method to process DOT series on group analysis using cognitive paradigms on the prefrontal cortex. Difficulties such as the fact that scalp-brain distances vary between subjects or cerebral activations are difficult to reproduce due to strategies used by the subjects to solve arithmetic problems are considered. T-images given by fMRI and DOT volume series analyzed in SPM show that at the functional level, both DOT and fMRI measures detect the same areas, although DOT provides complementary information to fMRI signals about cerebral activity.

Diffuse optical tomography using multi-directional sources and detectors

Diffuse optical tomography (DOT) is an advanced imaging method used to visualize the internal state of biological tissues as 3D images. However, current continuous-wave DOT requires high-density probe arrays for measurement (less than 15-mm interval) to gather enough information for 3D image reconstruction, which makes the experiment time-consuming. In this paper, we propose a novel DOT measurement system using multi-directional light sources and multi-directional photodetectors instead of high-density probe arrays. We evaluated this system's multi-directional DOT through computer simulation and a phantom experiment. From the results, we achieved DOT with less than 5-mm localization error up to a 15-mm depth with low-density probe arrays (30-mm interval), indicating that the multi-directional measurement approach allows DOT without requiring high-density measurement.