Direct estimation of evoked hemoglobin changes by multimodality fusion imaging

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

Significance

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

Aim

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

Approach

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

Results

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

Conclusions

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

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

SIGNIFICANCE

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

AIMS

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

APPROACH

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

RESULTS

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

CONCLUSIONS

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

Increased cortical activation and enhanced functional connectivity in the prefrontal cortex ensure dynamic postural balance during dual-task obstacle negotiation in the older adults: A fNIRS study

OBJECTIVE

By analyzing the cortical activation and functional connectivity of the prefrontal cortex (PFC) during dual-task obstacle negotiation in the older adults, cognitive resources allocation and neural regulatory mechanisms of aging brain were shed light on in complex walking conditions.

METHODS

Twenty-eight healthy right-handed subjects participated in the study, including 15 men and 13 women (age: 68.6 ± 4.1 years, height: 162.96 ± 6.05 cm, weight: 63.63 ± 9.64 kg). There were four tasks: Normal Walk (NW), Obstacle Negotiation during Normal Walk (NW + ON), Walk while performing Cognitive Task (WCT), and Obstacle Negotiation during Walk while performing Cognitive Task (WCT + ON). Participants wore functional near-infrared spectroscopy (fNIRS) to collect hemodynamic signals from various regions of interest (ROIs) in the PFC, while the three-dimensional motion capture system was used to test the gait velocity. Cognitive task data was recorded.

RESULTS

In WCT + ON, the HbO₂ concentration change value (△HbO₂) of the left dorsolateral prefrontal cortex was significantly greater than that in the other three tasks (p < 0.05), and the△HbO₂ of the right dorsolateral prefrontal cortex was significantly greater than that in NW + ON (p < 0.05). The gait velocities in the four tasks were significantly different (p < 0.05) (NW > WCT > NW + ON > WCT + ON). There was no significant difference in cognitive performance between in the WCT and WCT + ON (p > 0.05). In WCT + ON, the left and right dorsolateral prefrontal areas had strong functional connectivity and the left frontal pole was most widely connected to the other ROIs. Compared to that in NW, the functional connectivity of the left prefrontal lobe was significantly enhanced in WCT + ON (p < 0.05).

CONCLUSIONS

As walking difficulty increased, the PFC activation in the older adults changed from right-sided to bilateral activation, indicating that the left PFC cognitive resources compensated for the right PFC in dual-task obstacle negotiation. The cognitive resources recruitment in dual-task obstacle negotiation might be achieved by synchronization and coordination of associated brain areas in the PFC, primarily to maintain dynamic postural balance when walking.

Best practices for fNIRS publications

The application of functional near-infrared spectroscopy (fNIRS) in the neurosciences has been expanding over the last 40 years. Today, it is addressing a wide range of applications within different populations and utilizes a great variety of experimental paradigms. With the rapid growth and the diversification of research methods, some inconsistencies are appearing in the way in which methods are presented, which can make the interpretation and replication of studies unnecessarily challenging. The Society for Functional Near-Infrared Spectroscopy has thus been motivated to organize a representative (but not exhaustive) group of leaders in the field to build a consensus on the best practices for describing the methods utilized in fNIRS studies. Our paper has been designed to provide guidelines to help enhance the reliability, repeatability, and traceability of reported fNIRS studies and encourage best practices throughout the community. A checklist is provided to guide authors in the preparation of their manuscripts and to assist reviewers when evaluating fNIRS papers.

Short-channel regression in functional near-infrared spectroscopy is more effective when considering heterogeneous scalp hemodynamics

Significance: The reliability of functional near-infrared spectroscopy (fNIRS) measurements is reduced by systemic physiology. Short-channel regression algorithms aim at removing systemic "noise" by subtracting the signal measured at a short source-detector separation (mainly scalp hemodynamics) from the one of a long separation (brain and scalp hemodynamics). In literature, incongruent approaches on the selection of the optimal regressor signal are reported based on different assumptions on scalp hemodynamics properties. Aim: We investigated the spatial and temporal distribution of scalp hemodynamics over the sensorimotor cortex and evaluated its influence on the effectiveness of short-channel regressions. Approach: We performed hand-grasping and resting-state experiments with five subjects, measuring with 16 optodes over sensorimotor areas, including eight 8-mm channels. We performed detailed correlation analyses of scalp hemodynamics and evaluated 180 hand-grasping and 270 simulated (overlaid on resting-state measurements) trials. Five short-channel regressor combinations were implemented with general linear models. Three were chosen according to literature, and two were proposed based on additional physiological assumptions [considering multiple short channels and their Mayer wave (MW) oscillations]. Results: We found heterogeneous hemodynamics in the scalp, coming on top of a global close-to-homogeneous behavior (correlation 0.69 to 0.92). The results further demonstrate that short-channel regression always improves brain activity estimates but that better results are obtained when heterogeneity is assumed. In particular, we highlight that short-channel regression is more effective when combining multiple scalp regressors and when MWs are additionally included. Conclusion: We shed light on the selection of optimal regressor signals for improving the removal of systemic physiological artifacts in fNIRS. We conclude that short-channel regression is most effective when assuming heterogeneous hemodynamics, in particular when combining spatial- and frequency-specific information. A better understanding of scalp hemodynamics and more effective short-channel regression will promote more accurate assessments of functional brain activity in clinical and research settings.

Investigation of brain functional connectivity in patients with mild cognitive impairment: A functional near-infrared spectroscopy (fNIRS) study

This study examines brain functional connectivity in both cognitively normal seniors and patients with mild cognitive impairment (MCI) to elucidate prospective markers of MCI. A homemade four-channel functional near-infrared spectroscopy (fNIRS) system was employed to measure hemodynamic responses in the subjects' prefrontal cortex during a resting state, an oddball task, a 1-back task, and a verbal fluency task. Brain functional connectivity was calculated as the Pearson correlation coefficients between fNIRS channels. The results show that during the verbal fluency task, while the healthy control (HC) group presents a significantly stronger inter-hemispheric connectivity compared to intra-hemispheric connectivity, there is no difference between the inter- and intra-hemispheric connectivity in the MCI group. In addition, a comparison between the MCI and HC connectivity reveals that the MCI group has a statistically higher right and inter-hemispheric connectivity during the resting state, but a significantly lower left and inter-hemispheric connectivity during the verbal fluency test. These findings demonstrate the potential of fNIRS to study brain functional connectivity in neurodegenerative diseases.

Brain activation during neurocognitive testing using functional near-infrared spectroscopy in patients following concussion compared to healthy controls

There is no accepted clinical imaging modality for concussion, and current imaging modalities including fMRI, DTI, and PET are expensive and inaccessible to most clinics/patients. Functional near-infrared spectroscopy (fNIRS) is a non-invasive, portable, and low-cost imaging modality that can measure brain activity. The purpose of this study was to compare brain activity as measured by fNIRS in concussed and age-matched controls during the performance of cognitive tasks from a computerized neurocognitive test battery. Participants included nine currently symptomatic patients aged 18-45 years with a recent (15-45 days) sport-related concussion and five age-matched healthy controls. The participants completed a computerized neurocognitive test battery while wearing the fNIRS unit. Our results demonstrated reduced brain activation in the concussed subject group during word memory, (spatial) design memory, digit-symbol substitution (symbol match), and working memory (X's and O's) tasks. Behavioral performance (percent-correct and reaction time respectively) was lower for concussed participants on the word memory, design memory, and symbol match tasks than controls. The results of this preliminary study suggest that fNIRS could be a useful, portable assessment tool to assess reduced brain activation and augment current approaches to assessment and management of patients following concussion.

Validation of the hypercapnic calibrated fMRI method using DOT-fMRI fusion imaging

Calibrated functional magnetic resonance imaging (fMRI) is a widely used method to investigate brain function in terms of physiological quantities such as the cerebral metabolic rate of oxygen (CMRO2). The first and one of the most common methods of fMRI calibration is hypercapnic calibration. This is achieved via simultaneous measures of the blood-oxygenation-level dependent (BOLD) and the arterial spin labeling (ASL) signals during a functional task that evokes regional changes in CMRO2. A subsequent acquisition is then required during which the subject inhales carbon dioxide for short periods of time. A calibration constant, typically labeled M, is then estimated from the hypercapnic data and is subsequently used together with the BOLD-ASL recordings to compute evoked changes in CMRO2 during the functional task. The computation of M assumes a constant CMRO2 during the CO2 inhalation, an assumption that has been questioned since the origin of calibrated fMRI. In this study we used diffuse optical tomography (DOT) together with BOLD and ASL--an alternative calibration method that does not require any gas manipulation and therefore no constant CMRO2 assumption--to cross-validate the estimation of M obtained from a traditional hypercapnic calibration. We found a high correlation between the M values (R=0.87, p<0.01) estimated using these two approaches. The findings serve to validate the hypercapnic fMRI calibration technique and suggest that the inter-subject variability routinely obtained for M is reproducible with an alternative method and might therefore reflect inter-subject physiological variability.

Environmental and Genetic Influences on Neurocognitive Development: The Importance of Multiple Methodologies and Time-Dependent Intervention

Genetic mutations and environmental factors dynamically influence gene expression and developmental trajectories at the neural, cognitive, and behavioral levels. The examples in this article cover different periods of neurocognitive development-early childhood, adolescence, and adulthood-and focus on studies in which researchers have used a variety of methodologies to illustrate the early effects of socioeconomic status and stress on brain function, as well as how allelic differences explain why some individuals respond to intervention and others do not. These studies highlight how similar behaviors can be driven by different underlying neural processes and show how a neurocomputational model of early development can account for neurodevelopmental syndromes, such as autism spectrum disorders, with novel implications for intervention. Finally, these studies illustrate the importance of the timing of environmental and genetic factors on development, consistent with our view that phenotypes are emergent, not predetermined.

Calibrating the BOLD signal during a motor task using an extended fusion model incorporating DOT, BOLD and ASL data

Multimodal imaging improves the accuracy of the localization and the quantification of brain activation when measuring different manifestations of the hemodynamic response associated with cerebral activity. In this study, we incorporated cerebral blood flow (CBF) changes measured with arterial spin labeling (ASL), Diffuse Optical Tomography (DOT) and blood oxygen level-dependent (BOLD) recordings to reconstruct changes in oxy- (ΔHbO(2)) and deoxyhemoglobin (ΔHbR). Using the Grubb relation between relative changes in CBF and cerebral blood volume (CBV), we incorporated the ASL measurement as a prior to the total hemoglobin concentration change (ΔHbT). We applied this ASL fusion model to both synthetic data and experimental multimodal recordings during a 2-s finger-tapping task. Our results show that the new approach is very powerful in estimating ΔHbO(2) and ΔHbR with high spatial and quantitative accuracy. Moreover, our approach allows the computation of baseline total hemoglobin concentration (HbT(0)) as well as of the BOLD calibration factor M on a single subject basis. We obtained an average HbT(0) of 71 μM, an average M value of 0.18 and an average increase of 13% in cerebral metabolic rate of oxygen (CMRO(2)), all of which are in agreement with values previously reported in the literature. Our method yields an independent measurement of M, which provides an alternative measurement to validate the hypercapnic calibration of the BOLD signal.

Molecular brain imaging in the multimodality era

Multimodality molecular brain imaging encompasses in vivo visualization, evaluation, and measurement of cellular/molecular processes. Instrumentation and software developments over the past 30 years have fueled advancements in multimodality imaging platforms that enable acquisition of multiple complementary imaging outcomes by either combined sequential or simultaneous acquisition. This article provides a general overview of multimodality neuroimaging in the context of positron emission tomography as a molecular imaging tool and magnetic resonance imaging as a structural and functional imaging tool. Several image examples are provided and general challenges are discussed to exemplify complementary features of the modalities, as well as important strengths and weaknesses of combined assessments. Alzheimer's disease is highlighted, as this clinical area has been strongly impacted by multimodality neuroimaging findings that have improved understanding of the natural history of disease progression, early disease detection, and informed therapy evaluation.

The developmental trajectory of brain-scalp distance from birth through childhood: implications for functional neuroimaging

Measurements of human brain function in children are of increasing interest in cognitive neuroscience. Many techniques for brain mapping used in children, including functional near-infrared spectroscopy (fNIRS), electroencephalography (EEG), magnetoencephalography (MEG) and transcranial magnetic stimulation (TMS), use probes placed on or near the scalp. The distance between the scalp and the brain is a key variable for these techniques because optical, electrical and magnetic signals are attenuated by distance. However, little is known about how scalp-brain distance differs between different cortical regions in children or how it changes with development. We investigated scalp-brain distance in 71 children, from newborn to age 12 years, using structural T1-weighted MRI scans of the whole head. Three-dimensional reconstructions were created from the scalp surface to allow for accurate calculation of brain-scalp distance. Nine brain landmarks in different cortical regions were manually selected in each subject based on the published fNIRS literature. Significant effects were found for age, cortical region and hemisphere. Brain-scalp distances were lowest in young children, and increased with age to up to double the newborn distance. There were also dramatic differences between brain regions, with up to 50% differences between landmarks. In frontal and temporal regions, scalp-brain distances were significantly greater in the right hemisphere than in the left hemisphere. The largest contributors to developmental changes in brain-scalp distance were increases in the corticospinal fluid (CSF) and inner table of the cranium. These results have important implications for functional imaging studies of children: age and brain-region related differences in fNIRS signals could be due to the confounding factor of brain-scalp distance and not true differences in brain activity.

Fluorescence molecular tomography: principles and potential for pharmaceutical research

Fluorescence microscopic imaging is widely used in biomedical research to study molecular and cellular processes in cell culture or tissue samples. This is motivated by the high inherent sensitivity of fluorescence techniques, the spatial resolution that compares favorably with cellular dimensions, the stability of the fluorescent labels used and the sophisticated labeling strategies that have been developed for selectively labeling target molecules. More recently, two and three-dimensional optical imaging methods have also been applied to monitor biological processes in intact biological organisms such as animals or even humans. These whole body optical imaging approaches have to cope with the fact that biological tissue is a highly scattering and absorbing medium. As a consequence, light propagation in tissue is well described by a diffusion approximation and accurate reconstruction of spatial information is demanding. While in vivo optical imaging is a highly sensitive method, the signal is strongly surface weighted, i.e., the signal detected from the same light source will become weaker the deeper it is embedded in tissue, and strongly depends on the optical properties of the surrounding tissue. Derivation of quantitative information, therefore, requires tomographic techniques such as fluorescence molecular tomography (FMT), which maps the three-dimensional distribution of a fluorescent probe or protein concentration. The combination of FMT with a structural imaging method such as X-ray computed tomography (CT) or Magnetic Resonance Imaging (MRI) will allow mapping molecular information on a high definition anatomical reference and enable the use of prior information on tissue's optical properties to enhance both resolution and sensitivity. Today many of the fluorescent assays originally developed for studies in cellular systems have been successfully translated for experimental studies in animals. The opportunity of monitoring molecular processes non-invasively in the intact organism is highly attractive from a diagnostic point of view but even more so for the drug developer, who can use the techniques for proof-of-mechanism and proof-of-efficacy studies. This review shall elucidate the current status and potential of fluorescence tomography including recent advances in multimodality imaging approaches for preclinical and clinical drug development.

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

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

A quantitative comparison of NIRS and fMRI across multiple cognitive tasks

Near infrared spectroscopy (NIRS) is an increasingly popular technology for studying brain function. NIRS presents several advantages relative to functional magnetic resonance imaging (fMRI), such as measurement of concentration changes in both oxygenated and deoxygenated hemoglobin, finer temporal resolution, and ease of administration, as well as disadvantages, most prominently inferior spatial resolution and decreased signal-to-noise ratio (SNR). While fMRI has become the gold standard for in vivo imaging of the human brain, in practice NIRS is a more convenient and less expensive technology than fMRI. It is therefore of interest to many researchers how NIRS compares to fMRI in studies of brain function. In the present study we scanned participants with simultaneous NIRS and fMRI on a battery of cognitive tasks, placing NIRS probes over both frontal and parietal brain regions. We performed detailed comparisons of the signals in both temporal and spatial domains. We found that NIRS signals have significantly weaker SNR, but are nonetheless often highly correlated with fMRI measurements. Both SNR and the distance between the scalp and the brain contributed to variability in the NIRS/fMRI correlations. In the spatial domain, we found that a photon path forming an ellipse between the NIRS emitter and detector correlated most strongly with the BOLD response. Taken together these findings suggest that, while NIRS can be an appropriate substitute for fMRI for studying brain activity related to cognitive tasks, care should be taken when designing studies with NIRS to ensure that: 1) the spatial resolution is adequate for answering the question of interest and 2) the design accounts for weaker SNR, especially in brain regions more distal from the scalp.

A random-effects model for group-level analysis of diffuse optical brain imaging

Diffuse optical imaging is a non-invasive technique for measuring changes in blood oxygenation in the brain. This technique is based on the temporally and spatially resolved recording of optical absorption in tissue within the near-infrared range of light. Optical imaging can be used to study functional brain activity similar to functional MRI. However, group level comparisons of brain activity from diffuse optical data are difficult due to registration of optical sensors between subjects. In addition, optical signals are sensitive to inter-subject differences in cranial anatomy and the specific arrangement of optical sensors relative to the underlying functional region. These factors can give rise to partial volume errors and loss of sensitivity and therefore must be accounted for in combining data from multiple subjects. In this work, we describe an image reconstruction approach using a parametric Bayesian model that simultaneously reconstructs group-level images of brain activity in the context of a random-effects analysis. Using this model, we demonstrate that localization accuracy and the statistical effects size of group-level reconstructions can be improved when compared to individualized reconstructions. In this model, we use the Restricted Maximum Likelihood (ReML) method to optimize a Bayesian random-effects model.

A quantitative comparison of NIRS and fMRI across multiple cognitive tasks

Near infrared spectroscopy (NIRS) is an increasingly popular technology for studying brain function. NIRS presents several advantages relative to functional magnetic resonance imaging (fMRI), such as measurement of concentration changes in both oxygenated and deoxygenated hemoglobin, finer temporal resolution, and ease of administration, as well as disadvantages, most prominently inferior spatial resolution and decreased signal-to-noise ratio (SNR). While fMRI has become the gold standard for in vivo imaging of the human brain, in practice NIRS is a more convenient and less expensive technology than fMRI. It is therefore of interest to many researchers how NIRS compares to fMRI in studies of brain function. In the present study we scanned participants with simultaneous NIRS and fMRI on a battery of cognitive tasks, placing NIRS probes over both frontal and parietal brain regions. We performed detailed comparisons of the signals in both temporal and spatial domains. We found that NIRS signals have significantly weaker SNR, but are nonetheless often highly correlated with fMRI measurements. Both SNR and the distance between the scalp and the brain contributed to variability in the NIRS/fMRI correlations. In the spatial domain, we found that a photon path forming an ellipse between the NIRS emitter and detector correlated most strongly with the BOLD response. Taken together these findings suggest that, while NIRS can be an appropriate substitute for fMRI for studying brain activity related to cognitive tasks, care should be taken when designing studies with NIRS to ensure that: 1) the spatial resolution is adequate for answering the question of interest and 2) the design accounts for weaker SNR, especially in brain regions more distal from the scalp.

Computer-aided interpretation approach for optical tomographic images

A computer-aided interpretation approach is proposed to detect rheumatic arthritis (RA) in human finger joints using optical tomographic images. The image interpretation method employs a classification algorithm that makes use of a so-called self-organizing mapping scheme to classify fingers as either affected or unaffected by RA. Unlike in previous studies, this allows for combining multiple image features, such as minimum and maximum values of the absorption coefficient for identifying affected and not affected joints. Classification performances obtained by the proposed method were evaluated in terms of sensitivity, specificity, Youden index, and mutual information. Different methods (i.e., clinical diagnostics, ultrasound imaging, magnet resonance imaging, and inspection of optical tomographic images), were used to produce ground truth benchmarks to determine the performance of image interpretations. Using data from 100 finger joints, findings suggest that some parameter combinations lead to higher sensitivities, while others to higher specificities when compared to single parameter classifications employed in previous studies. Maximum performances are reached when combining the minimum/maximum ratio of the absorption coefficient and image variance. In this case, sensitivities and specificities over 0.9 can be achieved. These values are much higher than values obtained when only single parameter classifications were used, where sensitivities and specificities remained well below 0.8.

Hierarchical Bayesian regularization of reconstructions for diffuse optical tomography using multiple priors

Diffuse optical tomography (DOT) is a non-invasive brain imaging technique that uses low-levels of near-infrared light to measure optical absorption changes due to regional blood flow and blood oxygen saturation in the brain. By arranging light sources and detectors in a grid over the surface of the scalp, DOT studies attempt to spatially localize changes in oxy- and deoxy-hemoglobin in the brain that result from evoked brain activity during functional experiments. However, the reconstruction of accurate spatial images of hemoglobin changes from DOT data is an ill-posed linearized inverse problem, which requires model regularization to yield appropriate solutions. In this work, we describe and demonstrate the application of a parametric restricted maximum likelihood method (ReML) to incorporate multiple statistical priors into the recovery of optical images. This work is based on similar methods that have been applied to the inverse problem for magnetoencephalography (MEG). Herein, we discuss the adaptation of this model to DOT and demonstrate that this approach provides a means to objectively incorporate reconstruction constraints and demonstrate this approach through a series of simulated numerical examples.

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

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

Neuroimaging of the developing brain: taking "developing" seriously

With a few notable exceptions, many studies, be they behavioral, neuroimaging, or genetic, are snapshots that compare one child group to one adult group, which capture only two points in time and tell the scientist nothing about the mechanisms underlying neural trajectories over developmental time. Thus, a distinction needs to be drawn between child neuroimaging and developmental neuroimaging, the latter approach being relevant not just to children, but to adults and the ageing brain.

Quantitative evaluation of high-density diffuse optical tomography: in vivo resolution and mapping performance

Despite the unique brain imaging capabilities and advantages of functional near-infrared spectroscopy (fNIRS), including portability and comprehensive hemodynamic measurement, widespread acceptance in the neuroimaging community has been hampered by low spatial resolution and image localization errors. While recent technical developments such as high-density diffuse optical tomography (HD-DOT) have, in principle, been shown to have superior in silico image quality, the majority of optical imaging studies are still conducted with sparse fNIRS arrays, perhaps partially because the performance increases of HD-DOT appear incremental. Without a quantitative comparative analysis between HD-DOT and fNIRS, using both simulation and in vivo neuroimaging, the implications of the new HD-DOT technology have been difficult to judge. We present a quantitative comparison of HD-DOT and two commonly used fNIRS geometries using (1) standard metrics of image quality, (2) simulated brain mapping tasks, and (3) in vivo visual cortex mapping results in adult humans. The results show that better resolution and lower positional errors are achieved with HD-DOT and that these improvements provide a substantial advancement in neuroimaging capability. In particular, we demonstrate that HD-DOT enables detailed phase-encoded retinotopic mapping, while sparse arrays are limited to imaging individual block-design visual stimuli.

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

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

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

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