Retinotopic mapping of adult human visual cortex with high-density diffuse optical tomography

Distinct Time-Resolved Brain-Wide Coactivations in Oxygenated and Deoxygenated Hemoglobin

OBJECTIVE

Human resting-state networks (RSNs) estimated from oxygenated (HbO) and deoxygenated hemoglobin (HbR) data exhibit strong similarities, while task-based studies show different dynamics in HbR and HbO responses. Such a discrepancy might be explained due to time-averaged estimations of RSNs. Our study investigated differences between HbO and HbR on time-resolved brain-wide coactivation patterns (CAPs).

METHODS

Diffuse optical tomography was reconstructed from resting-state whole-head functional near-infrared spectroscopy data of HbR and HbO in individual healthy participants. Time-averaged RSNs were obtained using the group-level independent component analysis. Time-resolved CAPs were estimated using a clustering approach on the time courses of all obtained RSNs. Characteristics of the RSNs and CAPs from HbR and HbO were compared.

RESULTS

Spatial patterns of HbR and HbO RSNs exhibited significant similarities. Meanwhile, HbR CAPs revealed much more organized spatial and dynamic characteristics than HbO CAPs. The entire set of HbR CAPs suggests a superstructure resulted from brain-wide neuronal dynamics, which is less evident in the set of HbO CAPs. These differences between HbO and HbR CAPs were consistently replicated in individual session data.

CONCLUSION

Our results suggest that human resting brain-wide neuronal activations are preserved better in time-resolved brain-wide patterns, i.e., CAPs, from HbR than those from HbO, while such a difference is lost between time-averaged HbR and HbO RSNs.

SIGNIFICANCE

Our results reveal, for the first time, HbR concentration fluctuations are more directly coupled with resting dynamics of brain-wide neuronal activations in human brains.

Identifying Naturalistic Movies from Human Brain Activity with High-Density Diffuse Optical Tomography

Modern neuroimaging modalities, particularly functional MRI (fMRI), can decode detailed human experiences. Thousands of viewed images can be identified or classified, and sentences can be reconstructed. Decoding paradigms often leverage encoding models that reduce the stimulus space into a smaller yet generalizable feature set. However, the neuroimaging devices used for detailed decoding are non-portable, like fMRI, or invasive, like electrocorticography, excluding application in naturalistic use. Wearable, non-invasive, but lower-resolution devices such as electroencephalography and functional near-infrared spectroscopy (fNIRS) have been limited to decoding between stimuli used during training. Herein we develop and evaluate model-based decoding with high-density diffuse optical tomography (HD-DOT), a higher-resolution expansion of fNIRS with demonstrated promise as a surrogate for fMRI. Using a motion energy model of visual content, we decoded the identities of novel movie clips outside the training set with accuracy far above chance for single-trial decoding. Decoding was robust to modulations of testing time window, different training and test imaging sessions, hemodynamic contrast, and optode array density. Our results suggest that HD-DOT can translate detailed decoding into naturalistic use.

Mapping brain function in adults and young children during naturalistic viewing with high-density diffuse optical tomography

Human studies of early brain development have been limited by extant neuroimaging methods. MRI scanners present logistical challenges for imaging young children, while alternative modalities like functional near-infrared spectroscopy have traditionally been limited by image quality due to sparse sampling. In addition, conventional tasks for brain mapping elicit low task engagement, high head motion, and considerable participant attrition in pediatric populations. As a result, typical and atypical developmental trajectories of processes such as language acquisition remain understudied during sensitive periods over the first years of life. We evaluate high-density diffuse optical tomography (HD-DOT) imaging combined with movie stimuli for high resolution optical neuroimaging in awake children ranging from 1 to 7 years of age. We built an HD-DOT system with design features geared towards enhancing both image quality and child comfort. Furthermore, we characterized a library of animated movie clips as a stimulus set for brain mapping and we optimized associated data analysis pipelines. Together, these tools could map cortical responses to movies and contained features such as speech in both adults and awake young children. This study lays the groundwork for future research to investigate response variability in larger pediatric samples and atypical trajectories of early brain development in clinical populations.

Anatomical Modeling and Optimization of Speckle Contrast Optical Tomography

Traditional methods for mapping cerebral blood flow (CBF), such as positron emission tomography and magnetic resonance imaging, offer only isolated snapshots of CBF due to scanner logistics. Speckle contrast optical tomography (SCOT) is a promising optical technique for mapping CBF. However, while SCOT has been established in mice, the method has not yet been demonstrated in humans - partly due to a lack of anatomical reconstruction methods and uncertainty over the optimal design parameters. Herein we develop SCOT reconstruction methods that leverage MRI-based anatomical head models and finite-element modeling of the SCOT forward problem (NIRFASTer). We then simulate SCOT for CBF perturbations to evaluate sensitivity of imaging performance to exposure time and SD-distances. We find image resolution comparable to intensity-based diffuse optical tomography at superficial cortical tissue depth (~1.5 cm). Localization errors can be reduced by including longer SD-measurements. With longer exposure times speckle contrast decreases, however, noise decreases faster, resulting in a net increase in SNR. Specifically, extending exposure time from 10μs to 10ms increased SCOT SNR by 1000X. Overall, our modeling methods provide anatomically-based image reconstructions that can be used to evaluate a broad range of tissue conditions, measurement parameters, and noise sources and inform SCOT system design.

Ultra-high density imaging arrays for diffuse optical tomography of human brain improve resolution, signal-to-noise, and information decoding

Functional magnetic resonance imaging (fMRI) has dramatically advanced non-invasive human brain mapping and decoding. Functional near-infrared spectroscopy (fNIRS) and high-density diffuse optical tomography (HD-DOT) non-invasively measure blood oxygen fluctuations related to brain activity, like fMRI, at the brain surface, using more-lightweight equipment that circumvents ergonomic and logistical limitations of fMRI. HD-DOT grids have smaller inter-optode spacing (∼13 mm) than sparse fNIRS (∼30 mm) and therefore provide higher image quality, with spatial resolution ∼1/2 that of fMRI. Herein, simulations indicated reducing inter-optode spacing to 6.5 mm would further improve image quality and noise-resolution tradeoff, with diminishing returns below 6.5 mm. We then constructed an ultra-high-density DOT system (6.5-mm spacing) with 140 dB dynamic range that imaged stimulus-evoked activations with 30-50% higher spatial resolution and repeatable multi-focal activity with excellent agreement with participant-matched fMRI. Further, this system decoded visual stimulus position with 19-35% lower error than previous HD-DOT, throughout occipital cortex.

Diffuse Optical Tomography Spatial Prior for EEG Source Localization in Human Visual Cortex

Electroencephalography (EEG) and diffuse optical tomography (DOT) are imaging methods which are widely used for neuroimaging. While the temporal resolution of EEG is high, the spatial resolution is typically limited. DOT, on the other hand, has high spatial resolution, but the temporal resolution is inherently limited by the slow hemodynamics it measures. In our previous work, we showed using computer simulations that when using the results of DOT reconstruction as the spatial prior for EEG source reconstruction, high spatio-temporal resolution could be achieved. In this work, we experimentally validate the algorithm by alternatingly flashing two visual stimuli at a speed that is faster than the temporal resolution of DOT. We show that the joint reconstruction using both EEG and DOT clearly resolves the two stimuli temporally, and the spatial confinement is drastically improved in comparison to reconstruction using EEG alone.

Scalable, modular continuous wave functional near-infrared spectroscopy system (Spotlight)

Significance

We present a fiberless, portable, and modular continuous wave-functional near-infrared spectroscopy system, Spotlight, consisting of multiple palm-sized modules-each containing high-density light-emitting diode and silicon photomultiplier detector arrays embedded in a flexible membrane that facilitates optode coupling to scalp curvature.

Aim

Spotlight's goal is to be a more portable, accessible, and powerful functional near-infrared spectroscopy (fNIRS) device for neuroscience and brain-computer interface (BCI) applications. We hope that the Spotlight designs we share here can spur more advances in fNIRS technology and better enable future non-invasive neuroscience and BCI research.

Approach

We report sensor characteristics in system validation on phantoms and motor cortical hemodynamic responses in a human finger-tapping experiment, where subjects wore custom 3D-printed caps with two sensor modules.

Results

The task conditions can be decoded offline with a median accuracy of 69.6%, reaching 94.7% for the best subject, and at a comparable accuracy in real time for a subset of subjects. We quantified how well the custom caps fitted to each subject and observed that better fit leads to more observed task-dependent hemodynamic response and better decoding accuracy.

Conclusions

The advances presented here should serve to make fNIRS more accessible for BCI applications.

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.

Multi-modulated frequency domain high density diffuse optical tomography

Frequency domain (FD) high density diffuse optical tomography (HD-DOT) utilising varying or combined modulation frequencies (mFD) has shown to theoretically improve the imaging accuracy as compared to conventional continuous wave (CW) measurements. Using intensity and phase data from a solid inhomogeneous phantom (NEUROPT) with three insertable rods containing different contrast anomalies, at modulation frequencies of 78 MHz, 141 MHz and 203 MHz, HD-DOT is applied and quantitatively evaluated, showing that mFD outperforms FD and CW for both absolute (iterative) and temporal (linear) tomographic imaging. The localization error (LOCA), full width half maximum (FWHM) and effective resolution (ERES) were evaluated. Across all rods, the LOCA of mFD was 61.3% better than FD and 106.1% better than CW. For FWHM, CW was 6.0% better than FD and mFD and for ERES, mFD was 1.20% better than FD and 9.83% better than CW. Using mFD data is shown to minimize the effect of inherently noisier FD phase data whilst maximising its strengths through improved contrast.

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.

Correction of global physiology in resting-state functional near-infrared spectroscopy

Significance: Resting-state functional connectivity (RSFC) analyses of functional near-infrared spectroscopy (fNIRS) data reveal cortical connections and networks across the brain. Motion artifacts and systemic physiology in evoked fNIRS signals present unique analytical challenges, and methods that control for systemic physiological noise have been explored. Whether these same methods require modification when applied to resting-state fNIRS (RS-fNIRS) data remains unclear. Aim: We systematically examined the sensitivity and specificity of several RSFC analysis pipelines to identify the best methods for correcting global systemic physiological signals in RS-fNIRS data. Approach: Using numerically simulated RS-fNIRS data, we compared the rates of true and false positives for several connectivity analysis pipelines. Their performance was scored using receiver operating characteristic analysis. Pipelines included partial correlation and multivariate Granger causality, with and without short-separation measurements, and a modified multivariate causality model that included a non-traditional zeroth-lag cross term. We also examined the effects of pre-whitening and robust statistical estimators on performance. Results: Consistent with previous work on bivariate correlation models, our results demonstrate that robust statistics and pre-whitening are effective methods to correct for motion artifacts and autocorrelation in the fNIRS time series. Moreover, we found that pre-filtering using principal components extracted from short-separation fNIRS channels as part of a partial correlation model was most effective in reducing spurious correlations due to shared systemic physiology when the two signals of interest fluctuated synchronously. However, when there was a temporal lag between the signals, a multivariate Granger causality test incorporating the short-separation channels was better. Since it is unknown if such a lag exists in experimental data, we propose a modified version of Granger causality that includes the non-traditional zeroth-lag term as a compromising solution. Conclusions: A combination of pre-whitening, robust statistical methods, and partial correlation in the processing pipeline to reduce autocorrelation, motion artifacts, and global physiology are suggested for obtaining statistically valid connectivity metrics with RS-fNIRS. Further studies should validate the effectiveness of these methods using human data.

Prefrontal cortex supports speech perception in listeners with cochlear implants

Cochlear implants are neuroprosthetic devices that can restore hearing in people with severe to profound hearing loss by electrically stimulating the auditory nerve. Because of physical limitations on the precision of this stimulation, the acoustic information delivered by a cochlear implant does not convey the same level of acoustic detail as that conveyed by normal hearing. As a result, speech understanding in listeners with cochlear implants is typically poorer and more effortful than in listeners with normal hearing. The brain networks supporting speech understanding in listeners with cochlear implants are not well understood, partly due to difficulties obtaining functional neuroimaging data in this population. In the current study, we assessed the brain regions supporting spoken word understanding in adult listeners with right unilateral cochlear implants (n=20) and matched controls (n=18) using high-density diffuse optical tomography (HD-DOT), a quiet and non-invasive imaging modality with spatial resolution comparable to that of functional MRI. We found that while listening to spoken words in quiet, listeners with cochlear implants showed greater activity in the left prefrontal cortex than listeners with normal hearing, specifically in a region engaged in a separate spatial working memory task. These results suggest that listeners with cochlear implants require greater cognitive processing during speech understanding than listeners with normal hearing, supported by compensatory recruitment of the left prefrontal cortex.

Different brain activation patterns in the prefrontal area between self-paced and high-speed driving tasks

The purpose of this study was to investigate the effects on prefrontal cortex brain activity when participants attempted to stop a car accurately at a stop line when driving at different speeds using functional near-infrared spectroscopy (fNIRS). Twenty healthy subjects with driving experience drove their own cars for a distance of 60 m five times each at their own pace or as fast as possible. The variation in the distance between the stop line and the car was not significantly different between the self-paced and high-speed tasks. However, oxygenated hemoglobin concentration in the prefrontal cortex was significantly higher in the high-speed task than in the self-paced task. These findings suggest that driving at high speed requires more divided attention than driving at self-paced speed, even though the participants were able to stop the car at the same distance from the target. This study shows the advantages and usefulness of fNIRS .

Photoacoustic imaging of squirrel monkey cortical and subcortical brain regions during peripheral electrical stimulation

The investigation of neuronal activity in non-human primate models is of critical importance due to their genetic similarity to human brains. In this study, we tested the feasibility of using photoacoustic imaging for the detection of cortical and subcortical responses due to peripheral electrical stimulation in a squirrel monkey model. Photoacoustic computed tomography and photoacoustic microscopy were applied on squirrel monkeys for real-time deep subcortical imaging and optical-resolution cortical imaging, respectively. The electrically evoked hemodynamic changes in primary somatosensory cortex, premotor cortices, primary motor cortex, and underlying subcortical areas were measured. Hemodynamic responses were observed in both cortical and subcortical brain areas at the cortices during external stimulation, demonstrating the feasibility of photoacoustic technique for functional imaging of non-human primate brain.

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\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\text {cm}^2$$\end{document}cm2 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.

Imaging affective and non-affective touch processing in two-year-old children

Touch is an important component of early parent-child interaction and plays a critical role in the socio-emotional development of children. However, there are limited studies on touch processing amongst children in the age range from one to three years. The present study used frequency-domain diffuse optical tomography (DOT) to investigate the processing of affective and non-affective touch over left frontotemporal brain areas contralateral to the stimulated forearm in two-year-old children. Affective touch was administered by a single stroke with a soft brush over the child's right dorsal forearm at 3 cm/s, while non-affective touch was provided by multiple brush strokes at 30 cm/s. We found that in the insula, the total haemoglobin (HbT) response to slow brushing was significantly greater than the response to fast brushing (slow > fast). Additionally, a region in the postcentral gyrus, Rolandic operculum and superior temporal gyrus exhibited greater response to fast brushing than slow brushing (fast > slow). These findings confirm that an adult-like pattern of haemodynamic responses to affective and non-affective touch can be recorded in two-year-old subjects using DOT. To improve the accuracy of modelling light transport in the two-year-old subjects, we used a published age-appropriate atlas and deformed it to match the exterior shape of each subject's head. We estimated the combined scalp and skull, and grey matter (GM) optical properties by fitting simulated data to calibrated and coupling error corrected phase and amplitude measurements. By utilizing a two-compartment cerebrospinal fluid (CSF) model, the accuracy of estimation of GM optical properties and the localization of activation in the insula was improved. The techniques presented in this paper can be used to study neural development of children at different ages and illustrate that the technology is well-tolerated by most two-year-old children and not excessively sensitive to subject movement. The study points the way towards exciting possibilities in functional imaging of deeper functional areas near sulci in small children.

Label-free photoacoustic computed tomography of mouse cortical responses to retinal photostimulation using a pair-wise correlation map

The lack of a non-invasive or minimally invasive imaging technique has long been a challenge to investigating brain activities in mice. Functional magnetic resonance imaging and the more recently developed diffuse optical imaging both suffer from limited spatial resolution. Photoacoustic (PA) imaging combines the sensitivity of optical excitation to hemodynamic changes and ultrasound detection's relatively high spatial resolution. In this study, we evaluated the feasibility of using a label-free, real-time PA computed tomography (PACT) system to measure visually evoked hemodynamic responses within the primary visual cortex (V1) in mice. Photostimulation of the retinas evoked significantly faster and stronger V1 responses in wild-type mice than in age-matched rod/cone-degenerate mice, consistent with known differences between rod/cone- vs. melanopsin-mediated photoreception. In conclusion, the PACT system in this study has sufficient sensitivity and spatial resolution to resolve visual cortical hemodynamics during retinal photostimulation, and PACT is a potential tool for investigating visually evoked brain activities in mouse models of retinal diseases.

A method for mapping retinal images in early visual cortical areas

The visual cortex has been a heavily studied region in neuroscience due to many factors, not the least of which is its well-defined retinotopic organization. This organization makes it possible to predict the general location of cortical regions that stimuli will activate during visual tasks. However, the precise and accurate mapping of these regions in human patients takes time, effort, and participant compliance that can be difficult in many patient populations. In humans, this retino-cortical mapping has typically been done using functional localizers which maximally activate the area of interest, and then the activation profile is thresholded and converted to a binary mask region of interest (ROI). An alternative method involves performing population receptive field (pRF) mapping of the whole visual field and choosing vertices whose pRF centers fall within the stimulus. This method ignores the spatial extent of the pRF which changes dramatically between central and peripheral vision. Both methods require a dedicated functional scan and depend on participants' stable fixation. The aim of this project was to develop a user-friendly method that can transform a retinal object of interest (for example, an image, a retinal lesion, or a preferred locus for fixation) from retinal space to its expected representation on the cortical surface without a functional scan. We modeled the retinal representation of each cortical vertex as a 2D Gaussian with a location and spatial extent given by a previously published retinotopic atlas. To identify how affected any cortical vertex would be by a given retinal object, we took the product of the retinal object with the Gaussian pRF of that cortical vertex. Normalizing this value gives the expected response of a given vertex to the retinal object. This method was validated using BOLD data obtained using a localizer with discrete visual stimuli, and showed good agreement to predicted values. Cortical localization of a visual stimulus or retinal defect can be obtained using our publicly available software, without a functional scan. Our software may benefit research with disease populations who have trouble maintaining stable fixation.

Ultrafast and Ultrahigh-Resolution Diffuse Optical Tomography for Brain Imaging with Sensitivity Equation based Noniterative Sparse Optical Reconstruction (SENSOR)

We introduce a novel image reconstruction method for time-resolved diffuse optical tomography (DOT) that yields submillimeter resolution in less than a second. This opens the door to high-resolution real-time DOT in imaging of the brain activity. We call this approach the sensitivity equation based noniterative sparse optical reconstruction (SENSOR) method. The high spatial resolution is achieved by implementing an asymptotic l ₀-norm operator that guarantees to obtain sparsest representation of reconstructed targets. The high computational speed is achieved by employing the nontruncated sensitivity equation based noniterative inverse formulation combined with reduced sensing matrix and parallel computing. We tested the new method with numerical and experimental data. The results demonstrate that the SENSOR algorithm can achieve 1 mm³ spatial-resolution optical tomographic imaging at depth of ∼60 mean free paths (MFPs) in 20∼30 milliseconds on an Intel Core i9 processor.

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.

Investigation of effect of modulation frequency on high-density diffuse optical tomography image quality

Significance: By incorporating multiple overlapping functional near-infrared spectroscopy (fNIRS) measurements, high-density diffuse optical tomography (HD-DOT) images human brain function with fidelity comparable to functional magnetic resonance imaging (fMRI). Previous work has shown that frequency domain high-density diffuse optical tomography (FD-HD-DOT) may further improve image quality over more traditional continuous wave (CW) HD-DOT. Aim: The effects of modulation frequency on image quality as obtainable with FD-HD-DOT is investigated through simulations with a realistic noise model of functional activations in human head models, arising from 11 source modulation frequencies between CW and 1000 MHz. Approach: Simulations were performed using five representative head models with an HD regular grid of 158 light sources and 166 detectors and an empirically derived noise model. Functional reconstructions were quantitatively assessed with multiple image quality metrics including the localization error (LE), success rate, full width at half maximum, and full volume at half maximum (FVHM). All metrics were evaluated against CW-based models. Results: Compared to CW, localization accuracy is improved by >40% throughout brain depths of 13 to 25 mm below the surface with 300 to 500 MHz modulation frequencies. Additionally, the reliable field of view in brain tissue is enlarged by 35% to 48% within an optimal frequency of 300 MHz after considering realistic noise, depending on the dynamic range of the system. Conclusions: These results point to the tremendous opportunities in further development of high bandwidth FD-HD-DOT system hardware for applications in human brain mapping.

Quantitative evaluation of frequency domain measurements in high density diffuse optical tomography

SIGNIFICANCE

High density diffuse optical tomography (HD-DOT) as applied in functional near-infrared spectroscopy (fNIRS) is largely limited to continuous wave (CW) data. Using a single modulation frequency, frequency domain (FD) HD-DOT has recently demonstrated better localization of focal activation as compared to CW data. We show that combining CW and FD measurements and multiple modulation frequencies increases imaging performance in fNIRS.

AIM

We evaluate the benefits of multiple modulation frequencies, combining different frequencies as well as CW data in fNIRS HD-DOT.

APPROACH

A layered model was used, with activation occurring within a cortex layer. CW and FD measurements were simulated at 78, 141, and 203 MHz with and without noise. The localization error, full width half maximum, and effective resolution were evaluated.

RESULTS

Across the average of the three metrics, at 141 MHz, FD performed 8.4% better than CW, and the combination of CW and FD was 21.7% better than CW. FD measurements at 203 MHz performed 5% better than 78 MHz. Moreover, the three combined modulation frequencies of FD and CW performed up to 3.92% better than 141 MHz alone.

CONCLUSIONS

We show that combining CW and FD measurements offers better performance than FD alone, with higher modulation frequencies increasing accuracy. Combining CW and FD measurements at multiple modulation frequencies yields the best overall performance.

Transcranial brain atlas-based optimization for functional near-infrared spectroscopy optode arrangement: Theory, algorithm, and application

The quality of optode arrangement is crucial for group imaging studies when using functional near-infrared spectroscopy (fNIRS). Previous studies have demonstrated the promising effectiveness of using transcranial brain atlases (TBAs), in a manual and intuition-based way, to guide optode arrangement when individual structural MRI data are unavailable. However, the theoretical basis of using TBA to optimize optode arrangement remains unclear, which leads to manual and subjective application. In this study, we first describe the theoretical basis of TBA-based optimization of optode arrangement using a mathematical framework. Second, based on the theoretical basis, an algorithm is proposed for automatically arranging optodes on a virtual scalp. The resultant montage is placed onto the head of each participant guided by a low-cost and portable navigation system. We compared our method with the widely used 10/20-system-assisted optode arrangement procedure, using finger-tapping and working memory tasks as examples of both low- and high-level cognitive systems. Performance, including optode montage designs, locations on each participant's scalp, brain activation, as well as ground truth indices derived from individual MRI data were evaluated. The results give convergent support for our method's ability to provide more accurate, consistent and efficient optode arrangements for fNIRS group imaging than the 10/20 method.

Neuroimaging of depression with diffuse optical tomography during repetitive transcranial magnetic stimulation

Repetitive transcranial magnetic stimulation (rTMS) is an effective and safe treatment for depression; however, its potential has likely been hindered due to non-optimized targeting, unclear ideal stimulation parameters, and lack of information regarding how the brain is physiologically responding during and after stimulation. While neuroimaging is ideal for obtaining such critical information, existing modalities have been limited due to poor resolutions, along with significant noise interference from the electromagnetic spectrum. In this study, we used a novel diffuse optical tomography (DOT) device in order to advance our understanding of the neurophysiological effects of rTMS in depression. Healthy and depressed subjects aged 18–70 were recruited. Treatment parameters were standardized with targeting of the left dorsolateral prefrontal cortex with a magnetic field intensity of 100% of motor threshold, pulse frequency of 10 per second, a 4 s stimulation time and a 26 s rest time. DOT imaging was simultaneously acquired from the contralateral dorsolateral prefrontal cortex. Six healthy and seven depressed subjects were included for final analysis. Hemoglobin changes and volumetric three-dimensional activation patterns were successfully captured. Depressed subjects were observed to have a delayed and less robust response to rTMS with a decreased volume of activation compared to healthy subjects. In this first-in-human study, we demonstrated the ability of DOT to safely and reliably capture and compare cortical response patterns to rTMS in depressed and healthy subjects. We introduced this emerging optical functional imaging modality as a novel approach to investigating targeting, new treatment parameters, and physiological effects of rTMS in depression.

Evaluating a new generation of wearable high-density diffuse optical tomography technology via retinotopic mapping of the adult visual cortex

Significance: High-density diffuse optical tomography (HD-DOT) has been shown to approach the resolution and localization accuracy of blood oxygen level dependent-functional magnetic resonance imaging in the adult brain by exploiting densely spaced, overlapping samples of the probed tissue volume, but the technique has to date required large and cumbersome optical fiber arrays. Aim: To evaluate a wearable HD-DOT system that provides a comparable sampling density to large, fiber-based HD-DOT systems, but with vastly improved ergonomics. Approach: We investigated the performance of this system by replicating a series of classic visual stimulation paradigms, carried out in one highly sampled participant during 15 sessions to assess imaging performance and repeatability. Results: Hemodynamic response functions and cortical activation maps replicate the results obtained with larger fiber-based systems. Our results demonstrate focal activations in both oxyhemoglobin and deoxyhemoglobin with a high degree of repeatability observed across all sessions. A comparison with a simulated low-density array explicitly demonstrates the improvements in spatial localization, resolution, repeatability, and image contrast that can be obtained with this high-density technology. Conclusions: The system offers the possibility for minimally constrained, spatially resolved functional imaging of the human brain in almost any environment and holds particular promise in enabling neuroscience applications outside of the laboratory setting. It also opens up new opportunities to investigate populations unsuited to traditional imaging technologies.

Processing pipeline for image reconstructed fNIRS analysis using both MRI templates and individual anatomy

Significance: Image reconstruction of fNIRS data is a useful technique for transforming channel-based fNIRS into a volumetric representation and managing spatial variance based on optode location. We present an innovative integrated pipeline for image reconstruction of fNIRS data using either MRI templates or individual anatomy. Aim: We demonstrate a pipeline with accompanying code to allow users to clean and prepare optode location information, prepare and standardize individual anatomical images, create the light model, run the 3D image reconstruction, and analyze data in group space. Approach: We synthesize a combination of new and existing software packages to create a complete pipeline, from raw data to analysis. Results: This pipeline has been tested using both templates and individual anatomy, and on data from different fNIRS data collection systems. We show high temporal correlations between channel-based and image-based fNIRS data. In addition, we demonstrate the reliability of this pipeline with a sample dataset that included 74 children as part of a longitudinal study taking place in Scotland. We demonstrate good correspondence between data in channel space and image reconstructed data. Conclusions: The pipeline presented here makes a unique contribution by integrating multiple tools to assemble a complete pipeline for image reconstruction in fNIRS. We highlight further issues that may be of interest to future software developers in the field.

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.

Signal regression in frequency-domain diffuse optical tomography to remove superficial signal contamination

Significance: Signal contamination is a major hurdle in functional near-infrared spectroscopy (fNIRS) of the human head as the NIR signal is contaminated with the changes corresponding to superficial tissue, therefore occluding the functional information originating from the cerebral region. For continuous wave, this is generally handled through linear regression of the shortest source-detector (SD) distance intensity measurement from all of the signals. Although phase measurements utilizing frequency domain (FD) provide deeper tissue sampling, the use of the shortest SD distance phase measurement for regression of superficial signal contamination can lead to misleading results, therefore suppressing cortical signals. Aim: An approach for FD fNIRS that utilizes a short-separation intensity signal directly to regress both intensity and phase measurements, providing a better regression of superficial signal contamination from both data-types, is proposed. Approach: Simulated data from realistic models of the human head are used, and signal regression using both intensity and phase-based components of the FD fNIRS is evaluated. Results: Intensity-based phase regression achieves a suppression of superficial signal contamination by 68% whereas phase-based phase regression is only by 13%. Phase-based phase regression is also shown to generate false-positive signals from the cortex, which are not desirable. Conclusions: Intensity-based phase regression provides a better methodology for minimizing superficial signal contamination in FD fNIRS.

Effect of Prolonged Pressure on Hemodynamics of Sacral Tissues Assessed by Diffuse Optical Imaging: A Pilot Study

Pressure injuries (PIs) are wounds resulting from prolonged pressure exerting on the skin and underlying tissues over bony prominences (e.g., lower back, heels, shoulders) in bed-bound patients and wheelchair users. Minimizing pressure has long been considered the most effective preventative method, and current guidelines require visual skin inspection and repositioning every two hours. However, these strategies are often applied deficiently and do not adequately prevent PIs from becoming penetrating wounds. Recent studies attribute the development of PIs to cell deformation, inflammatory, and ischemic damages that cumulatively propagate from the microscale (death of few cells) to the macroscale (tissue necrosis) within one to several hours. Although the nature of the PI pathogenesis is complex and multifactorial, measuring tissue alterations in real-time may elucidate the origination mechanism and ultimately allow detecting PIs at the earliest stage. In this pilot study, we evaluated the ability of diffuse optical imaging (DOI) to assess hemodynamic changes resulting from prolonged pressure on the sacral tissues in five healthy volunteers laying immobile in a supine position for 2 hours. A thin, body-conforming optical imaging probe encompassing 256 optodes arranged in a regularly spaced grid over a 160 × 160 mm area was used to construct DOI volumetric images representing changes of oxyhemoglobin (HbO₂) and deoxyhemoglobin (HHb) concentration from a zeroed baseline. After 2 hours of continuous body weight pressure, hemodynamic images in all subjects were substantially dissimilar from their individual baseline. We also found that hemodynamic similarity computed pairwise across subjects exhibited a high value and limited variability around the mean, thus denoting a consistent level of image similarity across subjects. These preliminary results indicate that prolonged pressure causes distinctive hemodynamic patterns that can be effectively investigated with DOI and that monitoring functional changes over time holds potential for clarifying the development mechanisms of PIs.

Amplitude of fNIRS Resting-State Global Signal Is Related to EEG Vigilance Measures: A Simultaneous fNIRS and EEG Study

Recently, functional near-infrared spectroscopy (fNIRS) has been utilized to image the hemodynamic activities and connectivity in the human brain. With the advantage of economic efficiency, portability, and fewer physical constraints, fNIRS enables studying of the human brain at versatile environment and various body positions, including at bed side and during exercise, which complements the use of functional magnetic resonance imaging (fMRI). However, like fMRI, fNIRS imaging can be influenced by the presence of a strong global component. Yet, the nature of the global signal in fNIRS has not been established. In this study, we investigated the relationship between fNIRS global signal and electroencephalogram (EEG) vigilance using simultaneous recordings in resting healthy subjects in high-density and whole-head montage. In Experiment 1, data were acquired at supine, sitting, and standing positions. Results found that the factor of body positions significantly affected the amplitude of the resting-state fNIRS global signal, prominently in the frequency range of 0.05-0.1 Hz but not in the very low frequency range of less than 0.05 Hz. As a control, the task-induced fNIRS or EEG responses to auditory stimuli did not differ across body positions. However, EEG vigilance plays a modulatory role in the fNIRS signals in the frequency range of less than 0.05 Hz: resting-state sessions of low EEG vigilance measures are associated with high amplitudes of fNIRS global signals. Moreover, in Experiment 2, we further examined the epoch-to-epoch fluctuations in concurrent fNIRS and EEG data acquired from a separate group of subjects and found a negative temporal correlation between EEG vigilance measures and fNIRS global signal amplitudes. Our study for the first time revealed that vigilance as a neurophysiological factor modulates the resting-state dynamics of fNIRS, which have important implications for understanding and processing the noises in fNIRS signals.

Cap-Based Transcranial Optical Tomography in an Awake Infant

Although Blood Oxygenation Level Dependent (BOLD) functional MRI (fMRI) is widely used to examine brain function in adults, the need for general anesthesia limits its practical utility in infants and small children. Functional Near-Infrared Spectroscopy - Diffuse Optical Tomography (fNIRS-DOT) imaging promises to be an alternative brain network imaging technique. Yet current versions of continuous-wave fNIRS-DOT systems are restricted to the cortical surface measurements and do not probe deep structures that are frequently injured especially in premature infants. Herein we report a transcranial near infrared optical imaging system, called Cap-based Transcranial Optical Tomography (CTOT) able to image whole brain hemodynamic activity with 3 seconds of data acquisition time. We show the system is capable of whole brain oxygenation mapping in an awake child, and that tomographically reconstructed static CTOT-derived oxy- and deoxygenated blood volumes are spatially correlated with the time-averaged BOLD fMRI volumes. By removing time bottlenecks in the current system, dynamic CTOT mapping should be possible, which would then enable evaluation of functional connectivity in awake infants.

Decoding visual information from high-density diffuse optical tomography neuroimaging data

Background:

Neural decoding could be useful in many ways, from serving as a neuroscience research tool to providing a means of augmented communication for patients with neurological conditions. However, applications of decoding are currently constrained by the limitations of traditional neuroimaging modalities. Electrocorticography requires invasive neurosurgery, magnetic resonance imaging (MRI) is too cumbersome for uses like daily communication, and alternatives like functional near-infrared spectroscopy (fNIRS) offer poor image quality. High-density diffuse optical tomography (HD-DOT) is an emerging modality that uses denser optode arrays than fNIRS to combine logistical advantages of optical neuroimaging with enhanced image quality. Despite the resulting promise of HD-DOT for facilitating field applications of neuroimaging, decoding of brain activity as measured by HD-DOT has yet to be evaluated.

Objective:

To assess the feasibility and performance of decoding with HD-DOT in visual cortex.

Methods and Results:

To establish the feasibility of decoding at the single-trial level with HD-DOT, a template matching strategy was used to decode visual stimulus position. A receiver operating characteristic (ROC) analysis was used to quantify the sensitivity, specificity, and reproducibility of binary visual decoding. Mean areas under the curve (AUCs) greater than 0.97 across 10 imaging sessions in a highly sampled participant were observed. ROC analyses of decoding across 5 participants established both reproducibility in multiple individuals and the feasibility of inter-individual decoding (mean AUCs > 0.7), although decoding performance varied between individuals. Phase-encoded checkerboard stimuli were used to assess more complex, non-binary decoding with HD-DOT. Across 3 highly sampled participants, the phase of a 60° wide checkerboard wedge rotating 10° per second through 360° was decoded with a within-participant error of 25.8±24.7°. Decoding between participants was also feasible based on permutation-based significance testing.

Conclusions:

Visual stimulus information can be decoded accurately, reproducibly, and across a range of detail (for both binary and non-binary outcomes) at the single-trial level (without needing to block-average test data) using HD-DOT data. These results lay the foundation for future studies of more complex decoding with HD-DOT and applications in clinical populations.

A wide field-of-view, modular, high-density diffuse optical tomography system for minimally constrained three-dimensional functional neuroimaging

The ability to produce high-quality images of human brain function in any environment and during unconstrained movement of the subject has long been a goal of neuroimaging research. Diffuse optical tomography, which uses the intensity of back-scattered near-infrared light from multiple source-detector pairs to image changes in haemoglobin concentrations in the brain, is uniquely placed to achieve this goal. Here, we describe a new generation of modular, fibre-less, high-density diffuse optical tomography technology that provides exceptional sensitivity, a large dynamic range, a field-of-view sufficient to cover approximately one-third of the adult scalp, and also incorporates dedicated motion sensing into each module. Using in-vivo measures, we demonstrate a noise-equivalent power of 318 fW, and an effective dynamic range of 142 dB. We describe the application of this system to a novel somatomotor neuroimaging paradigm that involves subjects walking and texting on a smartphone. Our results demonstrate that wearable high-density diffuse optical tomography permits three-dimensional imaging of the human brain function during overt movement of the subject; images of somatomotor cortical activation can be obtained while subjects move in a relatively unconstrained manner, and these images are in good agreement with those obtained while the subjects remain stationary. The scalable nature of the technology we described here paves the way for the routine acquisition of high-quality, three-dimensional, whole-cortex diffuse optical tomography images of cerebral haemodynamics, both inside and outside of the laboratory environment, which has profound implications for neuroscience.

Fixation-related NIRS indexes retinotopic occipital processing of parafoveal preview during natural reading

While word frequency and predictability effects have been examined extensively, any evidence on interactive effects as well as parafoveal influences during whole sentence reading remains inconsistent and elusive. Novel neuroimaging methods utilize eye movement data to account for the hemodynamic responses of very short events such as fixations during natural reading. In this study, we used the rapid sampling frequency of near-infrared spectroscopy (NIRS) to investigate neural responses in the occipital and orbitofrontal cortex to word frequency and predictability. We observed increased activation in the right ventral occipital cortex when the fixated word N was of low frequency, which we attribute to an enhanced cost during saccade planning. Importantly, unpredictable (in contrast to predictable) low frequency words increased the activity in the left dorsal occipital cortex at the fixation of the preceding word N-1, presumably due to an upcoming breach of top-down modulated expectation. Opposite to studies that utilized a serial presentation of words (e.g. Hofmann et al., 2014), we did not find such an interaction in the orbitofrontal cortex, implying that top-down timing of cognitive subprocesses is not required during natural reading. We discuss the implications of an interactive parafoveal-on-foveal effect for current models of eye movements.

Advantages of double density alignment of fNIRS optodes to evaluate cortical activities related to phonological short-term memory using NIRS-SPM

Phonological short-term memory (STM) and working memory (WM) capacity, which is temporal storage of phonological information, facilitates language development. This capacity is often impaired in cochlear implant (CI) users who have congenital hearing loss. Functional near-infrared spectroscopy (fNIRS) is potentially useful to reveal underlying mechanisms of this impairment due to its tolerance to magnetic and electrical artifacts generated by CIs. The spatial resolution of the standard fNIRS, however, seems inadequate to evaluate cortical activity associated with the maintenance of phonological information in the STM/WM tasks. In the present study recruiting 14 normal hearing adults, we applied a double density alignment of fNIRS optodes to improve spatial resolution in the generalized linear model (GLM)-based statistical analysis using NIRS-SPM, in which cortical activities were estimated during each of three stages (encoding, maintenance, and retrieval) in pseudoword STM tasks with auditory or visual presentation. Since the double density alignment of fNIRS optodes contains two sets of standard density arrays, in the off-line analysis the measured cortical hemodynamic responses can be analyzed as data from two independent standard density arrays as well as those from one double density array (DD-array) which has two times higher density of channels. The two standard arrays demonstrated a similar pattern of cortical activation at each stage of the auditory and visual tasks, which proved the reliability of our fNIRS analysis, but failed to detect significant cortical activation in 2 of 12 conditions including the maintenance stage in the visual task. On the other hand, DD-array revealed significant cortical activation in all conditions. These differences were observed when estimated cortical activation was localized in small regions, which suggests higher spatial resolution in DD-array than the standard arrays. In our knowledge, this is the first clinical study supporting the previous experimental phantom study which demonstrated improvement of spatial resolution in the double density arrangement of fNIRS optodes. These findings imply that the double density alignment of fNIRS optodes improves reliability and spatial resolution in fNIRS-based estimation of cortical activity in the STM/WM studies, although further studies are required to determine usefulness of this method in the CI population.

Using anatomically defined regions-of-interest to adjust for head-size and probe alignment in functional near-infrared spectroscopy

Significance: Functional near-infrared spectroscopy (fNIRS) uses surface-placed light sources and detectors to record underlying changes in the brain due to fluctuations in hemoglobin levels and oxygenation. Since these measurements are recorded from the surface of the scalp, the mapping from underlying regions-of-interest (ROIs) in the brain space to the fNIRS channel space measurements depends on the registration of the sensors, the anatomy of the head/brain, and the sensitivity of these diffuse measurements through the tissue. However, small displacements in the probe position can change the distribution of recorded brain activity across the fNIRS measurements. Aim: We propose an approach using either individual or atlas-based brain-space anatomical information to define ROI-based statistical hypotheses to test the null involvement of specific regions, which allows us to test the analogous ROI across subjects while adjusting for fNIRS probe placement and sensitivity differences due to head size variations without a localizer task. Approach: We use the optical forward model to project the underlying brain-space ROI into a tapered contrast vector, which defines the relative weighting of the fNIRS channels contributing to the ROI and allows us to test the null hypothesis of no brain activity in this region during a functional task. We demonstrate this method through simulation and compare the sensitivity-specificity of this approach to other conventional methods. Results: We examine the performance of this method in the scenario where head size and probe registration are both an accurately known parameters and where this is subject to unknown experimental errors. This method is compared with the performance of the conventional method using 364 different simulation parameter combinations. Conclusion: The proposed method is always recommended in ROI-based analysis, since it significantly improves the analysis performance without a localizer task, wherever the fNIRS probe registration is known or unknown.

4D Functional Imaging of the Rat Brain Using a Large Aperture Row-Column Array

Functional ultrasound imaging (fUS) recently emerged as a promising neuroimaging modality to image and monitor brain activity based on cerebral blood volume response (CBV) and neurovascular coupling. fUS offers very good spatial and temporal resolutions compared to fMRI gold standard as well as simplicity and portability. It was recently extended to 4D fUS imaging in preclinical settings although this approach remains limited and complex. Indeed 4D fUS requires a 2D matrix probe and specific hardware able to drive the N² elements of the probe with thousands of electronic channels. Several under-sampling approaches are currently investigated to limit the channel count and spread ultrasound 4D modalities. Among them, the Row Column Addressing (RCA) approach combined with ultrafast imaging is a compelling alternative using only N + N channels. We present a large field of view RCA probe prototype of 128 + 128 channels and 15 MHz central frequency adapted for preclinical imaging. Based on the Orthogonal Plane Wave compounding scheme, we were able to perform 4D vascular brain acquisitions at high volume rate. Doppler volumes of the whole rat brain were obtained in vivo at high rates (23 dB CNR at 156 Hz and 19 dB CNR at 313 Hz). Visual and whiskers stimulations were performed and the corresponding CBV increases were reconstructed in 3D with successful functional activation detected in the superior colliculus and somato-sensorial cortex respectively. This proof of concept study demonstrates for the first time the use of a low-channel count RCA array for in vivo 4D fUS imaging in the whole rat brain.

Retinal Prosthetic Approaches to Enhance Visual Perception for Blind Patients

Retinal prostheses are implantable devices that aim to restore the vision of blind patients suffering from retinal degeneration, mainly by artificially stimulating the remaining retinal neurons. Some retinal prostheses have successfully reached the stage of clinical trials; however, these devices can only restore vision partially and remain insufficient to enable patients to conduct everyday life independently. The visual acuity of the artificial vision is limited by various factors from both engineering and physiological perspectives. To overcome those issues and further enhance the visual resolution of retinal prostheses, a variety of retinal prosthetic approaches have been proposed, based on optimization of the geometries of electrode arrays and stimulation pulse parameters. Other retinal stimulation modalities such as optics, ultrasound, and magnetics have also been utilized to address the limitations in conventional electrical stimulation. Although none of these approaches have been clinically proven to fully restore the function of a degenerated retina, the extensive efforts made in this field have demonstrated a series of encouraging findings for the next generation of retinal prostheses, and these could potentially enhance the visual acuity of retinal prostheses. In this article, a comprehensive and up-to-date overview of retinal prosthetic strategies is provided, with a specific focus on a quantitative assessment of visual acuity results from various retinal stimulation technologies. The aim is to highlight future directions toward high-resolution retinal prostheses.

The present and future use of functional near-infrared spectroscopy (fNIRS) for cognitive neuroscience

The past few decades have seen a rapid increase in the use of functional near-infrared spectroscopy (fNIRS) in cognitive neuroscience. This fast growth is due to the several advances that fNIRS offers over the other neuroimaging modalities such as functional magnetic resonance imaging and electroencephalography/magnetoencephalography. In particular, fNIRS is harmless, tolerant to bodily movements, and highly portable, being suitable for all possible participant populations, from newborns to the elderly and experimental settings, both inside and outside the laboratory. In this review we aim to provide a comprehensive and state-of-the-art review of fNIRS basics, technical developments, and applications. In particular, we discuss some of the open challenges and the potential of fNIRS for cognitive neuroscience research, with a particular focus on neuroimaging in naturalistic environments and social cognitive neuroscience.

Imaging System Based on Silicon Photomultipliers and Light Emitting Diodes for Functional Near-Infrared Spectroscopy

We built a fiber-less prototype of an optical system with 156 channels each one consisting of an optode made of a silicon photomultiplier (SiPM) and a pair of light emitting diodes (LEDs) operating at 700 nm and 830 nm. The system uses functional near-infrared spectroscopy (fNIRS) and diffuse optical tomography (DOT) imaging of the cortical activity of the human brain at frequencies above 1 Hz. In this paper, we discuss testing and system optimization performed through measurements on a multi-layered optical phantom with mechanically movable parts that simulate near-infrared light scattering inhomogeneities. The baseline optical characteristics of the phantom are carefully characterized and compared to those of human tissues. Here we discuss several technical aspects of the system development, such as LED light output drift and its possible compensation, SiPM linearity, corrections of channel signal differences, and signal-to-noise ratio (SNR). We implement an imaging algorithm that investigates large phantom regions. Thanks to the use of SiPMs, very large source-to-detector distances are acquired with a high SNR and 2 Hz time resolution. The overall results demonstrate the high potentialities of a system based on SiPMs for fNIRS/DOT human brain imaging applications.

Relationship between maternal pregnancy-related anxiety and infant brain responses to emotional speech - a pilot study

BACKGROUND

Maternal pregnancy-related anxiety (PRA) is reportedly related to neurodevelopmental outcomes of infants. However, the relationship between maternal PRA and the processing of emotions in the infant brain has not been extensively studied with neuroimaging. The objective of the present pilot study is to investigate the relationship between maternal PRA and infant hemodynamic responses to emotional speech at two months of age.

METHODS

The study sample included 19 mother-infant dyads from a general sample of a population of Caucasian mothers. Self-reported Pregnancy-Related Anxiety Questionnaire (PRAQ-R2) data was collected from mothers during pregnancy at gestational weeks (gwks) 24 (N = 19) and 34 (N = 18). When their infants were two months old, the infants' brains functional responses to emotional speech in the left fronto-temporoparietal cortex were recorded using diffuse optical tomography (DOT).

RESULTS

Maternal PRAQ-R2 scores at gwk 24 correlated negatively with the total hemoglobin (HbT) responses to sad speech on both sides of the temporoparietal junction (Spearman's rank correlation coefficient ρ = -0.87). The correlation was significantly greater at gwk 24 than gwk 34 (ρ = -0.42).

LIMITATIONS

The field of view of the measurement did not include the right hemisphere or parts of the frontal cortex. The sample size is moderate and the mothers were relatively highly educated, thus there may be some differences between the study sample and the general population.

CONCLUSIONS

Maternal pregnancy-related anxiety may affect child brain emotion processing development. Further research is needed to understand the functional and developmental significance of the findings.

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.

Portable, field-based neuroimaging using high-density diffuse optical tomography

Behavioral and cognitive tests in individuals who were malnourished as children have revealed malnutrition-related deficits that persist throughout the lifespan. These findings have motivated recent neuroimaging investigations that use highly portable functional near-infrared spectroscopy (fNIRS) instruments to meet the demands of brain imaging experiments in low-resource environments and enable longitudinal investigations of brain function in the context of long-term malnutrition. However, recent studies in healthy subjects have demonstrated that high-density diffuse optical tomography (HD-DOT) can significantly improve image quality over that obtained with sparse fNIRS imaging arrays. In studies of both task activations and resting state functional connectivity, HD-DOT is beginning to approach the data quality of fMRI for superficial cortical regions. In this work, we developed a customized HD-DOT system for use in malnutrition studies in Cali, Colombia. Our results evaluate the performance of the HD-DOT instrument for assessing brain function in a cohort of malnourished children. In addition to demonstrating portability and wearability, we show the HD-DOT instrument's sensitivity to distributed brain responses using a sensory processing task and measurements of homotopic functional connectivity. Task-evoked responses to the passive word listening task produce activations localized to bilateral superior temporal gyrus, replicating previously published work using this paradigm. Evaluating this localization performance across sparse and dense reconstruction schemes indicates that greater localization consistency is associated with a dense array of overlapping optical measurements. These results provide a foundation for additional avenues of investigation, including identifying and characterizing a child's individual malnutrition burden and eventually contributing to intervention development.

Influence of contrast-reversing frequency on the amplitude and spatial distribution of visual cortex hemodynamic responses

Visual stimulation is one of the most commonly used paradigms for cerebral cortex function investigation. Experiments typically involve presenting to a volunteer a black-and-white checkerboard with contrast-reversing at a frequency of 4 to 16 Hz. The aim of the present study was to investigate the influence of the flickering frequency on the amplitude of changes in the concentration of oxygenated and deoxygenated hemoglobin. The hemoglobin concentrations were assessed with the use of a high resolution diffuse optical tomography method. Spatial distributions of changes in hemoglobin concentrations overlaying the visual cortex are shown for various stimuli frequencies. Moreover, the hemoglobin concentration changes obtained for different source-detector separations (from 1.5 to 5.4 cm) are presented. Our results demonstrate that the flickering frequency had a statistically significant effect on the induced oxyhemoglobin changes (p < 0,001). The amplitude of oxy hemoglobin concentration changes at a frequency of 8 Hz was higher in comparison with that measured at 4 Hz :[median(25th-75thpercentiles) 1.24 (0.94-1.71) vs. 0.92(0.73-1.28)µM, p < 0.001]; 12 Hz:[1.24 (0.94-1.71) vs. 1.04 (0.78-1.32) µM, p < 0.001]; and 16 Hz:[1.24 (0.94-1.71) vs. 1.15(0.87-1.48) µM, p < 0.001]. No significant differences were observed between the size of an area of activation for various frequencies. The demonstrated superiority of 8 Hz over other frequencies can advance understanding of visual stimulations and help guide future fNIRS protocols.

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.

Functional near infrared spectroscopy using spatially resolved data to account for tissue scattering: A numerical study and arm-cuff experiment

Functional Near-Infrared Spectroscopy (fNIRS) aims to recover changes in tissue optical parameters relating to tissue hemodynamics, to infer functional information in biological tissue. A widely-used application of fNIRS relies on continuous wave (CW) methodology that utilizes multiple distance measurements on human head for study of brain health. The typical method used is spatially resolved spectroscopy (SRS), which is shown to recover tissue oxygenation index (TOI) based on gradient of light intensity measured between two detectors. However, this methodology does not account for tissue scattering which is often assumed. A new parameter recovery algorithm is developed, which directly recovers both the scattering parameter and scaled chromophore concentrations and hence TOI from the measured gradient of light-attenuation at multiple wavelengths. It is shown through simulations that in comparison to conventional SRS which estimates cerebral TOI values with an error of ±12.3%, the proposed method provides more accurate estimate of TOI exhibiting an error of ±5.7% without any prior assumptions of tissue scatter, and can be easily implemented within CW fNIRS systems. Using an arm-cuff experiment, the obtained TOI using the proposed method is shown to provide a higher and more realistic value as compared to utilizing any prior assumptions of tissue scatter.

Approach to optimize 3-dimensional brain functional activation image with high resolution: a study on functional near-infrared spectroscopy

In this study, 3-dimensional (3-D) enhanced brain-function-map generation and estimation methodology is presented. Optical signals were modelled in the form of numerical optimization problem to infer the best existing waveform of canonical hemodynamic response function. Inter-channel activity patterns were also estimated. The estimation of activation of inter-channel gap depends on the minimization of generalized cross-validation. 3-D brain activation maps were produced through inverse discrete cosine transform. The proposed algorithm acquired significant results for 3-D functional maps with high resolution, in comparison with that of 2-D functional t-maps. A comprehensive analysis by exhibiting images corresponding to several layers has also been appended.

Mapping brain function during naturalistic viewing using high-density diffuse optical tomography

Naturalistic stimuli, such as movies, more closely recapitulate "real life" sensory processing and behavioral demands relative to paradigms that rely on highly distilled and repetitive stimulus presentations. The rich complexity inherent in naturalistic stimuli demands an imaging system capable of measuring spatially distributed brain responses, and analysis tools optimized for unmixing responses to concurrently presented features. In this work, the combination of passive movie viewing with high-density diffuse optical tomography (HD-DOT) is developed as a platform for naturalistic brain mapping. We imaged healthy young adults during free viewing of a feature film using HD-DOT and observed reproducible, synchronized cortical responses across a majority of the field-of-view, most prominently in hierarchical cortical areas related to visual and auditory processing, both within and between individuals. In order to more precisely interpret broad patterns of cortical synchronization, we extracted visual and auditory features from the movie stimulus and mapped the cortical responses to the features. The results demonstrate the sensitivity of HD-DOT to evoked responses during naturalistic viewing, and that feature-based decomposition strategies enable functional mapping of naturalistic stimulus processing, including human-generated speech.

Hyper-spectral Recovery of Cerebral and Extra-Cerebral Tissue Properties Using Continuous Wave Near-Infrared Spectroscopic Data

Near-infrared spectroscopy (NIRS) is widely used as a non-invasive method to monitor the hemodynamics of biological tissue. A common approach of NIRS relies on continuous wave (CW) methodology, i.e. utilizing intensity-only measurements, and, in general, assumes homogeneity in the optical properties of the biological tissue. However, in monitoring cerebral hemodynamics within humans, this assumption is not valid due to complex layered structure of the head. The NIRS signal that contains information about cerebral blood hemoglobin levels is also contaminated with extra-cerebral tissue hemodynamics, and any recovery method based on such a priori homogenous approximation would lead to erroneous results. In this work, utilization of hyper-spectral intensity only measurements (i.e., CW) at multiple distances are presented and are shown to recover two-layered tissue properties along with the thickness of top layer, using an analytical solution for a two-layered semi-infinite geometry. It is demonstrated that the recovery of tissue oxygenation index (TOI) of both layers can be achieved with an error of 4.4%, with the recovered tissue thickness of 4% error. When the data is measured on a complex tissue such as the human head, it is shown that the semi-infinite recovery model can lead to uncertain results, whereas, when using an appropriate model accounting for the tissue-boundary structure, the tissue oxygenation levels are recovered with an error of 4.2%, and the extra-cerebral tissue thickness with an error of 11.8%. The algorithm is finally used together with human subject data, demonstrating robustness in application and repeatability in the recovered parameters that adhere well to expected published parameters.

High-density functional diffuse optical tomography based on frequency-domain measurements improves image quality and spatial resolution

Measurements of dynamic near-infrared (NIR) light attenuation across the human head together with model-based image reconstruction algorithms allow the recovery of three-dimensional spatial brain activation maps. Previous studies using high-density diffuse optical tomography (HD-DOT) systems have reported improved image quality over sparse arrays. These HD-DOT systems incorporated multidistance overlapping continuous wave measurements that only recover differential intensity attenuation. We investigate the potential improvement in reconstructed image quality due to the additional incorporation of phase shift measurements, which reflect the time-of-flight of the measured NIR light, within the tomographic reconstruction from high-density measurements. To evaluate image reconstruction with and without the additional phase information, we simulated point spread functions across a whole-scalp field of view in 24 subject-specific anatomical models using an experimentally derived noise model. The addition of phase information improves the image quality by reducing localization error by up to 59% and effective resolution by up to 21% as compared to using the intensity attenuation measurements alone. Furthermore, we demonstrate that the phase data enable images to be resolved at deeper brain regions where intensity data fail, which is further supported by utilizing experimental data from a single subject measurement during a retinotopic experiment.