A long-term storable gel-laden chip composite built in a multi-well plate enablingin situcell encapsulation for high-throughput liver model.

Hydrogels are widely used as scaffold materials for constructingin vitrothree-dimensional microphysiological systems. However, their high sensitivity to various external cues hinders the development of hydrogel-laden, microscale, and high-throughput chips. Here, we have developed a long-term storable gel-laden chip composite built in a multi-well plate, which enablesin situcell encapsulation and facilitates high-throughput analysis. Through optimized chemical crosslinking and freeze-drying method (C/FD), we have achieved a high-quality of gel-laden chip composite with excellent transparency, uniform porosity, and appropriate swelling and mechanical characteristics. Besides collagen, decellularized extracellular matrix with tissue-specific biochemical compound has been applied as chip composite. As a ready-to-use platform,in situcell encapsulation within the gel has been achieved through capillary force generated during gel reswelling. The liver-mimetic chip composite, comprising HepG2 cells or primary hepatocytes, has demonstrated favorable hepatic functionality and high sensitivity in drug testing. The developed fabrication process with improved stability of gels and storability allows chip composites to be stored at a wide range of temperatures for up to 28 d without any deformation, demonstrating off-the-shelf products. Consequently, this provides an exceptionally simple and long-term storable platform that can be utilized for an efficient tissue-specific modeling and various biomedical applications.

Association Between Taurine Level in the Hippocampus and Major Depressive Disorder in Young Women: A Proton Magnetic Resonance Spectroscopy Study at 7T.

BACKGROUND: Major depressive disorder (MDD) is characterized by depressed mood or loss of interest or pleasure. Generally, women are twice as likely as men to have depression. Taurine, a type of amino acid, plays critical roles in neuronal generation, differentiation, arborization, and formation of synaptic connections. Importantly, it enhances proliferation and synaptogenesis in the hippocampus. When injected into animals, taurine has an antidepressant effect. However, there is no in vivo evidence to show an association between taurine concentration in the human brain and the development of MDD. METHODS: Forty-one unmedicated young women with MDD (ages 18-29) and 43 healthy control participants matched for gender and age were recruited in South Korea. Taurine concentration was measured in the hippocampus, anterior cingulate cortex, and occipital cortex of the MDD and healthy control groups using proton magnetic resonance spectroscopy at 7T. Analysis of covariance was used to examine differences in taurine concentration, adjusting for age as a covariate. RESULTS: Taurine concentration in the hippocampus was lower (F1,75 = 5.729, p = .019, Δη2 = 0.073) for the MDD group (mean [SEM] = 0.91 [0.06] mM) than for the healthy control group (1.13 [0.06] mM). There was no significant difference in taurine concentration in the anterior cingulate cortex or occipital cortex between the two groups. CONCLUSIONS: This study demonstrates that a lower level of taurine concentration in the hippocampus may be a novel characteristic of MDD.

Task-residual effective connectivity of motor network in transient ischemic attack.

Transient ischemic attack (TIA) is a temporary episode of neurological dysfunction that results from focal brain ischemia. Although TIA symptoms are quickly resolved, patients with TIA have a high risk of stroke and persistent impairments in multiple domains of cognitive and motor functions. In this study, using spectral dynamic causal modeling, we investigate the changes in task-residual effective connectivity of patients with TIA during fist-closing movements. 28 healthy participants and 15 age-matched patients with TIA undergo functional magnetic resonance imaging at 7T. Here we show that during visually cued motor movement, patients with TIA have significantly higher effective connectivity toward the ipsilateral primary motor cortex and lower connectivity to the supplementary motor area than healthy controls. Our results imply that TIA patients have aberrant connections among motor regions, and these changes may reflect the decreased efficiency of primary motor function and disrupted control of voluntary movement in patients with TIA.

Kinetic evidence in favor of glyoxalase III and against deglycase activity of DJ-1.

DJ-1, a protein encoded by PARK7 plays a protective role against neurodegeneration. Since its glyoxalase III activity catalyzing methylglyoxal (MG) to lactate was discovered, DJ-1 has been re-established as a deglycase decomposing the MG-intermediates with amino acids and nucleotides (hemithioacetal and hemiaminal) rather than MG itself, but it is still debatable. Here, we have clarified that human DJ-1 directly recognizes MG, and not MG-intermediates, by monitoring the detailed catalytic processes and enantiomeric lactate products. The hemithioacetal intermediate between C106 of 15 N-labeled DJ-1 (15N DJ-1) and MG was also monitored by NMR. TRIS molecule formed stable diastereotopic complexes with MG (Kd , 1.57 ± 0.27 mM) by utilizing its three OH groups, which likely disturbed the assay of deglycase activity. The low kcat of DJ-1 for MG and its MG-induced structural perturbation may suggest that DJ-1 has a regulatory function as an in vivo sensor of reactive carbonyl stress.

A high-throughput biomimetic bone-on-a-chip platform with artificial intelligence-assisted image analysis for osteoporosis drug testing.

Although numerous organ-on-a-chips have been developed, bone-on-a-chip platforms have rarely been reported because of the high complexity of the bone microenvironment. With an increase in the elderly population, a high-risk group for bone-related diseases such as osteoporosis, it is essential to develop a precise bone-mimicking model for efficient drug screening and accurate evaluation in preclinical studies. Here, we developed a high-throughput biomimetic bone-on-a-chip platform combined with an artificial intelligence (AI)-based image analysis system. To recapitulate the key aspects of natural bone microenvironment, mouse osteocytes (IDG-SW3) and osteoblasts (MC3T3-E1) were cocultured within the osteoblast-derived decellularized extracellular matrix (OB-dECM) built in a well plate-based three-dimensional gel unit. This platform spatiotemporally and configurationally mimics the characteristics of the structural bone unit, known as the osteon. Combinations of native and bioactive ingredients obtained from the OB-dECM and coculture of two types of bone cells synergistically enhanced osteogenic functions such as osteocyte differentiation and osteoblast maturation. This platform provides a uniform and transparent imaging window that facilitates the observation of cell-cell interactions and features high-throughput bone units in a well plate that is compatible with a high-content screening system, enabling fast and easy drug tests. The drug efficacy of anti-SOST antibody, which is a newly developed osteoporosis drug for bone formation, was tested via ß-catenin translocation analysis, and the performance of the platform was evaluated using AI-based deep learning analysis. This platform could be a cutting-edge translational tool for bone-related diseases and an efficient alternative to bone models for the development of promising drugs.

Prediction of anticancer drug resistance using a 3D microfluidic bladder cancer model combined with convolutional neural network-based image analysis.

Bladder cancer is the most common urological malignancy worldwide, and its high recurrence rate leads to poor survival outcomes. The effect of anticancer drug treatment varies significantly depending on individual patients and the extent of drug resistance. In this study, we developed a validation system based on an organ-on-a-chip integrated with artificial intelligence technologies to predict resistance to anticancer drugs in bladder cancer. As a proof-of-concept, we utilized the gemcitabine-resistant bladder cancer cell line T24 with four distinct levels of drug resistance (parental, early, intermediate, and late). These cells were co-cultured with endothelial cells in a 3D microfluidic chip. A dataset comprising 2,674 cell images from the chips was analyzed using a convolutional neural network (CNN) to distinguish the extent of drug resistance among the four cell groups. The CNN achieved 95.2% accuracy upon employing data augmentation and a step decay learning rate with an initial value of 0.001. The average diagnostic sensitivity and specificity were 90.5% and 96.8%, respectively, and all area under the curve (AUC) values were over 0.988. Our proposed method demonstrated excellent performance in accurately identifying the extent of drug resistance, which can assist in the prediction of drug responses and in determining the appropriate treatment for bladder cancer patients.

A deep convolutional neural network for estimating hemodynamic response function with reduction of motion artifacts in fNIRS.

Objective. Functional near-infrared spectroscopy (fNIRS) is a neuroimaging technique for monitoring hemoglobin concentration changes in a non-invasive manner. However, subject movements are often significant sources of artifacts. While several methods have been developed for suppressing this confounding noise, the conventional techniques have limitations on optimal selections of model parameters across participants or brain regions. To address this shortcoming, we aim to propose a method based on a deep convolutional neural network (CNN).Approach. The U-net is employed as a CNN architecture. Specifically, large-scale training and testing data are generated by combining variants of hemodynamic response function (HRF) with experimental measurements of motion noises. The neural network is then trained to reconstruct hemodynamic response coupled to neuronal activity with a reduction of motion artifacts.Main results. Using extensive analysis, we show that the proposed method estimates the task-related HRF more accurately than the existing methods of wavelet decomposition and autoregressive models. Specifically, the mean squared error and variance of HRF estimates, based on the CNN, are the smallest among all methods considered in this study. These results are more prominent when the semi-simulated data contain variants of shapes and amplitudes of HRF.Significance. The proposed CNN method allows for accurately estimating amplitude and shape of HRF with significant reduction of motion artifacts. This method may have a great potential for monitoring HRF changes in real-life settings that involve excessive motion artifacts.

Errata: Best practices for fNIRS publications.

[This corrects the article DOI: 10.1117/1.NPh.8.1.012101.].

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.

Altered Effective Connectivity within the Fronto-Limbic Circuitry in Response to Negative Emotional Task in Female Patients with Major Depressive Disorder.

Background: Major depressive disorder (MDD) is a mood disorder associated with disruptions in emotional control. Previous studies have investigated abnormal regional activity and connectivity within the fronto-limbic circuit. However, condition-specific connectivity changes and their association with the pathophysiology of MDD remain unexplored. This study investigated effective connectivity in the fronto-limbic circuit induced by negative emotional processing from patients with MDD. Methods: Thirty-four unmedicated female patients with MDD and 28 healthy participants underwent event-related functional magnetic resonance imaging at 7T while viewing emotionally negative and neutral images. Brain regions whose dynamics are driven by experimental conditions were identified by using statistical parametric mapping. Effective connectivity among regions of interest was then estimated by using dynamic causal modeling. Results: Patients with MDD had lower activation of the orbitofrontal cortex (OFC) and higher activation of the parahippocampal gyrus (PHG) than healthy controls (HC). In association with these regional changes, we found that patients with MDD did not have significant modulatory connections from the primary visual cortex (V1) to OFC, whereas those connections of HC were significantly positively modulated during negative emotional processing. Regarding the PHG activity, patients with MDD had greater modulatory connection from the V1, but reduced negative modulatory connection from the OFC, compared with healthy participants. Conclusions: These results imply that disrupted effective connectivity among regions of the OFC, PHG, and V1 may be closely associated with the impaired regulation of negative emotional processing in the female patients with MDD.

Effects of Autonomous Sensory Meridian Response on the Functional Connectivity as Measured by Functional Magnetic Resonance Imaging.

Autonomous sensory meridian response (ASMR) is a sensory phenomenon in which audio-visual stimuli evoke a tingling sensation and is accompanied by a feeling of calm and relaxation. Therefore, there has been an increasing interest in using stimuli that elicit ASMR in cognitive and clinical neuroscience studies. However, neurophysiological basis of sensory-emotional experiences evoked by ASMR remain largely unexplored. In this study, we investigated how functional connectivity is changed while watching ASMR video, compared to resting state, and assessed its potential association with affective state induced by ASMR. 28 subjects participated in fMRI experiment consisting of 2 sessions (resting-state and task of viewing ASMR-eliciting video). Using a seed-based correlation analysis, we found that functional connections between the posterior cingulate cortex, and superior/middle temporal gyri, cuneus, and lingual gyrus were significantly increased during ASMR compared to resting state. In addition, we found that with the pregenual anterior cingulate cortex seed region, functional connectivity of the medial prefrontal cortex was increased during ASMR condition, relative to resting state. These results imply that ASMR can be elicited and maintained by ongoing interaction between regional activity that are mainly involved in the mentalizing and self-referential processing. We also found that ASMR-induced affective state changes (high activation negative and high activation positive state) were negatively correlated with functional connectivity involved in visual information processing, suggesting that visual information processing in response to high arousal states can be weakened by ASMR-eliciting stimuli.

Atrophic Changes and Diffusion Abnormalities of Affected Trigeminal Nerves in Trigeminal Neuralgia Using 7-T MRI.

BACKGROUND: Magnetic resonance imaging (MRI) has been widely used for visualizing trigeminal nerves in trigeminal neuralgia. OBJECTIVE: To assess atrophy and diffusion abnormalities of affected trigeminal nerves in trigeminal neuralgia with 7-T MRI. METHODS: In this prospective study, 14 patients (mean age 49 years; range 31-64 years) with trigeminal neuralgia underwent 7-T MRI. We measured trigeminal nerve volumes along their course through the pontocerebellar cistern. We also evaluated fractional anisotropy (FA) and quantitative anisotropy (QA) values within cisternal segment and pontine nuclei of the affected-side and unaffected-side trigeminal nerves, using diffusion tensor imaging (DTI). Associations between DTI metrics and Barrow Neurological Institute (BNI) pain scores were examined. RESULTS: The volumes were significantly smaller for the affected trigeminal nerves (33.83 ± 23.12 mm3) than for the unaffected ones (47.76 ± 32.48 mm3; p = 0.008). Cisternal segment FA and QA values were significantly lower in affected trigeminal nerves than in unaffected ones. However, DTI measurements in the pontine nuclei revealed no significant differences between affected-side and unaffected-side trigeminal nerves. No DTI metrics significantly correlated with BNI pain scores. CONCLUSION: Our results suggest that 7-T MRI allows identifications of atrophy and diffusion abnormalities of trigeminal nerves in trigeminal neuralgia.

Deep Residual Learning for Accelerated MRI Using Magnitude and Phase Networks.

OBJECTIVE: Accelerated magnetic resonance (MR) image acquisition with compressed sensing (CS) and parallel imaging is a powerful method to reduce MR imaging scan time. However, many reconstruction algorithms have high computational costs. To address this, we investigate deep residual learning networks to remove aliasing artifacts from artifact corrupted images. METHODS: The deep residual learning networks are composed of magnitude and phase networks that are separately trained. If both phase and magnitude information are available, the proposed algorithm can work as an iterative k-space interpolation algorithm using framelet representation. When only magnitude data are available, the proposed approach works as an image domain postprocessing algorithm. RESULTS: Even with strong coherent aliasing artifacts, the proposed network successfully learned and removed the aliasing artifacts, whereas current parallel and CS reconstruction methods were unable to remove these artifacts. CONCLUSION: Comparisons using single and multiple coil acquisition show that the proposed residual network provides good reconstruction results with orders of magnitude faster computational time than existing CS methods. SIGNIFICANCE: The proposed deep learning framework may have a great potential for accelerated MR reconstruction by generating accurate results immediately.

Dynamic causal modelling on infant fNIRS data: A validation study on a simultaneously recorded fNIRS-fMRI dataset.

Tracking the connectivity of the developing brain from infancy through childhood is an area of increasing research interest, and fNIRS provides an ideal method for studying the infant brain as it is compact, safe and robust to motion. However, data analysis methods for fNIRS are still underdeveloped compared to those available for fMRI. Dynamic causal modelling (DCM) is an advanced connectivity technique developed for fMRI data, that aims to estimate the coupling between brain regions and how this might be modulated by changes in experimental conditions. DCM has recently been applied to adult fNIRS, but not to infants. The present paper provides a proof-of-principle for the application of this method to infant fNIRS data and a demonstration of the robustness of this method using a simultaneously recorded fMRI-fNIRS single case study, thereby allowing the use of this technique in future infant studies. fMRI and fNIRS were simultaneously recorded from a 6-month-old sleeping infant, who was presented with auditory stimuli in a block design. Both fMRI and fNIRS data were preprocessed using SPM, and analysed using a general linear model approach. The main challenges that adapting DCM for fNIRS infant data posed included: (i) the import of the structural image of the participant for spatial pre-processing, (ii) the spatial registration of the optodes on the structural image of the infant, (iii) calculation of an accurate 3-layer segmentation of the structural image, (iv) creation of a high-density mesh as well as (v) the estimation of the NIRS optical sensitivity functions. To assess our results, we compared the values obtained for variational Free Energy (F), Bayesian Model Selection (BMS) and Bayesian Model Average (BMA) with the same set of possible models applied to both the fMRI and fNIRS datasets. We found high correspondence in F, BMS, and BMA between fMRI and fNIRS data, therefore showing for the first time high reliability of DCM applied to infant fNIRS data. This work opens new avenues for future research on effective connectivity in infancy by contributing a data analysis pipeline and guidance for applying DCM to infant fNIRS data.

Effective connectivity during working memory and resting states: A DCM study.

Although the relationship between resting-state functional connectivity and task-related activity has been addressed, the relationship between task and resting-state directed or effective connectivity - and its behavioral concomitants - remains elusive. We evaluated effective connectivity under an N-back working memory task in 24 participants using stochastic dynamic causal modelling (DCM) of 7 T fMRI data. We repeated the analysis using resting-state data, from the same subjects, to model connectivity among the same brain regions engaged by the N-back task. This allowed us to: (i) examine the relationship between intrinsic (task-independent) effective connectivity during resting (Arest) and task states (Atask), (ii) cluster phenotypes of task-related changes in effective connectivity (Btask) across participants, (iii) identify edges (Btask) showing high inter-individual effective connectivity differences and (iv) associate reaction times with the similarity between Btask and Arest in these edges. We found a strong correlation between Arest and Atask over subjects but a marked difference between Btask and Arest. We further observed a strong clustering of individuals in terms of Btask, which was not apparent in Arest. The task-related effective connectivity Btask varied highly in the edges from the parietal to the frontal lobes across individuals, so the three groups were clustered mainly by the effective connectivity within these networks. The similarity between Btask and Arest at the edges from the parietal to the frontal lobes was positively correlated with 2-back reaction times. This result implies that a greater change in context-sensitive coupling - from resting-state connectivity - is associated with faster reaction times. In summary, task-dependent connectivity endows resting-state connectivity with a context sensitivity, which predicts the speed of information processing during the N-back task.

Sparse SPM: Group Sparse-dictionary learning in SPM framework for resting-state functional connectivity MRI analysis.

Recent studies of functional connectivity MR imaging have revealed that the default-mode network activity is disrupted in diseases such as Alzheimer's disease (AD). However, there is not yet a consensus on the preferred method for resting-state analysis. Because the brain is reported to have complex interconnected networks according to graph theoretical analysis, the independency assumption, as in the popular independent component analysis (ICA) approach, often does not hold. Here, rather than using the independency assumption, we present a new statistical parameter mapping (SPM)-type analysis method based on a sparse graph model where temporal dynamics at each voxel position are described as a sparse combination of global brain dynamics. In particular, a new concept of a spatially adaptive design matrix has been proposed to represent local connectivity that shares the same temporal dynamics. If we further assume that local network structures within a group are similar, the estimation problem of global and local dynamics can be solved using sparse dictionary learning for the concatenated temporal data across subjects. Moreover, under the homoscedasticity variance assumption across subjects and groups that is often used in SPM analysis, the aforementioned individual and group analyses using sparse dictionary learning can be accurately modeled by a mixed-effect model, which also facilitates a standard SPM-type group-level inference using summary statistics. Using an extensive resting fMRI data set obtained from normal, mild cognitive impairment (MCI), and Alzheimer's disease patient groups, we demonstrated that the changes in the default mode network extracted by the proposed method are more closely correlated with the progression of Alzheimer's disease.

Associations of resting-state fMRI functional connectivity with flow-BOLD coupling and regional vasculature.

There has been tremendous interest in applying functional magnetic resonance imaging-based resting-state functional connectivity (rs-fcMRI) measurements to the study of brain function. However, a lack of understanding of the physiological mechanisms of rs-fcMRI limits their ability to interpret rs-fcMRI findings. In this work, the authors examine the regional associations between rs-fcMRI estimates and dynamic coupling between the blood oxygenation level-dependent (BOLD) and cerebral blood flow (CBF), as well as resting macrovascular volume. Resting-state BOLD and CBF data were simultaneously acquired using a dual-echo pseudocontinuous arterial spin labeling (pCASL) technique, whereas macrovascular volume fraction was estimated using time-of-flight MR angiography. Functional connectivity within well-known functional networks­including the default mode, frontoparietal, and primary sensory-motor networks­was calculated using a conventional seed-based correlation approach. They found the functional connectivity strength to be significantly correlated with the regional increase in CBF-BOLD coupling strength and inversely proportional to macrovascular volume fraction. These relationships were consistently observed within all functional networks considered. Their findings suggest that highly connected networks observed using rs-fcMRI are not likely to be mediated by common vascular drainage linking distal cortical areas. Instead, high BOLD functional connectivity is more likely to reflect tighter neurovascular connections, attributable to neuronal pathways.

Dynamic and static contributions of the cerebrovasculature to the resting-state BOLD signal.

Functional magnetic resonance imaging (fMRI) in the resting state, particularly fMRI based on the blood-oxygenation level-dependent (BOLD) signal, has been extensively used to measure functional connectivity in the brain. However, the mechanisms of vascular regulation that underlie the BOLD fluctuations during rest are still poorly understood. In this work, using dual-echo pseudo-continuous arterial spin labeling and MR angiography (MRA), we assess the spatio-temporal contribution of cerebral blood flow (CBF) to the resting-state BOLD signals and explore how the coupling of these signals is associated with regional vasculature. Using a general linear model analysis, we found that statistically significant coupling between resting-state BOLD and CBF fluctuations is highly variable across the brain, but the coupling is strongest within the major nodes of established resting-state networks, including the default-mode, visual, and task-positive networks. Moreover, by exploiting MRA-derived large vessel (macrovascular) volume fraction, we found that the degree of BOLD-CBF coupling significantly decreased as the ratio of large vessels to tissue volume increased. These findings suggest that the portion of resting-state BOLD fluctuations at the sites of medium-to-small vessels (more proximal to local neuronal activity) is more closely regulated by dynamic regulations in CBF, and that this CBF regulation decreases closer to large veins, which are more distal to neuronal activity.

Statistical analysis of fNIRS data: a comprehensive review.

Functional near-infrared spectroscopy (fNIRS) is a non-invasive method to measure brain activities using the changes of optical absorption in the brain through the intact skull. fNIRS has many advantages over other neuroimaging modalities such as positron emission tomography (PET), functional magnetic resonance imaging (fMRI), or magnetoencephalography (MEG), since it can directly measure blood oxygenation level changes related to neural activation with high temporal resolution. However, fNIRS signals are highly corrupted by measurement noises and physiology-based systemic interference. Careful statistical analyses are therefore required to extract neuronal activity-related signals from fNIRS data. In this paper, we provide an extensive review of historical developments of statistical analyses of fNIRS signal, which include motion artifact correction, short source-detector separation correction, principal component analysis (PCA)/independent component analysis (ICA), false discovery rate (FDR), serially-correlated errors, as well as inference techniques such as the standard t-test, F-test, analysis of variance (ANOVA), and statistical parameter mapping (SPM) framework. In addition, to provide a unified view of various existing inference techniques, we explain a linear mixed effect model with restricted maximum likelihood (ReML) variance estimation, and show that most of the existing inference methods for fNIRS analysis can be derived as special cases. Some of the open issues in statistical analysis are also described.

Lipschitz-Killing curvature based expected Euler characteristics for p-value correction in fNIRS.

Functional near-infrared spectroscopy (fNIRS) is a non-invasive imaging approach for measuring brain activities based on changes in the cerebral concentrations of hemoglobin. Recently, statistical analysis based on a general linear model (GLM) has become popular. Here, to impose statistical significance on the activation detected by fNIRS, family-wise error (FWE) rate control is important. However, unlike fMRI, in which measurements are densely sampled on a regular lattice and Gaussian smoothing makes the resulting random field homogeneous, the random fields from fNIRS are inhomogeneous due to the interpolation from sparsely and irregularly distributed optode locations. Thus, tube formula based correction has been proposed to address this issue. However, Sun's tube formula cannot be used for general random fields such as F-statistics. To overcome these difficulties, we employ the expected Euler characteristic approach based on Lipschitz-Killing curvature (LKC) to control the family-wise error rate. We compared this correction method with Sun's tube formula for t-statistics to confirm the existing method. Based on this comparison, we show that covariance estimation should be modified to consider channel-wise least-square residual correlation. These new results supplement the existing tool of statistical parameter mapping for fNIRS.