Depth Sensitivity and Source-Detector Separations for Near Infrared Spectroscopy Based on the Colin27 Brain Template

Understanding the spatial and depth sensitivity of non-invasive near-infrared spectroscopy (NIRS) measurements to brain tissue–i.e., near-infrared neuromonitoring (NIN) – is essential for designing experiments as well as interpreting research findings. However, a thorough characterization of such sensitivity in realistic head models has remained unavailable. In this study, we conducted 3,555 Monte Carlo (MC) simulations to densely cover the scalp of a well-characterized, adult male template brain (Colin27). We sought to evaluate: (i) the spatial sensitivity profile of NIRS to brain tissue as a function of source-detector separation, (ii) the NIRS sensitivity to brain tissue as a function of depth in this realistic and complex head model, and (iii) the effect of NIRS instrument sensitivity on detecting brain activation. We found that increasing the source-detector (SD) separation from 20 to 65 mm provides monotonic increases in sensitivity to brain tissue. For every 10 mm increase in SD separation (up to ∼45 mm), sensitivity to gray matter increased an additional 4%. Our analyses also demonstrate that sensitivity in depth (S) decreases exponentially, with a “rule-of-thumb” formula S = 0.75*0.85^depth. Thus, while the depth sensitivity of NIRS is not strictly limited, NIN signals in adult humans are strongly biased towards the outermost 10–15 mm of intracranial space. These general results, along with the detailed quantitation of sensitivity estimates around the head, can provide detailed guidance for interpreting the likely sources of NIRS signals, as well as help NIRS investigators design and plan better NIRS experiments, head probes and instruments.     

Spatiotemporal Rank Filtering Improves Image Quality Compared to Frame Averaging in 2-Photon Laser Scanning Microscopy

Live imaging of biological specimens using optical microscopy is limited by tradeoffs between spatial and temporal resolution, depth into intact samples, and phototoxicity. Two-photon laser scanning microscopy (2P-LSM), the gold standard for imaging turbid samples in vivo, has conventionally constructed images with sufficient signal-to-noise ratio (SNR) generated by sequential raster scans of the focal plane and temporal integration of the collected signals. Here, we describe spatiotemporal rank filtering, a nonlinear alternative to temporal integration, which makes more efficient use of collected photons by selectively reducing noise in 2P-LSM images during acquisition. This results in much higher SNR while preserving image edges and fine details. Practically, this allows for at least a four fold decrease in collection times, a substantial improvement for time-course imaging in biological systems.     

Modelling and detecting deep brain activity with MEG and EEG

We introduce an anatomical and electrophysiological model of deep brain structures dedicated to magnetoencephalography (MEG) and electroencephalography (EEG) source imaging. So far, most imaging inverse models considered that MEG/EEG surface signals were predominantly produced by cortical, hence superficial, neural currents. Here we question whether crucial deep brain structures such as the basal ganglia and the hippocampus may also contribute to distant, scalp MEG and EEG measurements. We first design a realistic anatomical and electrophysiological model of these structures and subsequently run Monte-Carlo experiments to evaluate the respective sensitivity of the MEG and EEG to signals from deeper origins. Results indicate that MEG/EEG may indeed localize these deeper generators, which is confirmed here from experimental MEG data reporting on the modulation of alpha (10–12 Hz) brain waves.     

Measuring Granger Causality between cortical regions from voxelwise fMRI BOLD signals with LASSO

Functional brain network studies using the Blood Oxygen-Level Dependent (BOLD) signal from functional Magnetic Resonance Imaging (fMRI) are becoming increasingly prevalent in research on the neural basis of human cognition. An important problem in functional brain network analysis is to understand directed functional interactions between brain regions during cognitive performance. This problem has important implications for understanding top-down influences from frontal and parietal control regions to visual occipital cortex in visuospatial attention, the goal motivating the present study. A common approach to measuring directed functional interactions between two brain regions is to first create nodal signals by averaging the BOLD signals of all the voxels in each region, and to then measure directed functional interactions between the nodal signals. Another approach, that avoids averaging, is to measure directed functional interactions between all pairwise combinations of voxels in the two regions. Here we employ an alternative approach that avoids the drawbacks of both averaging and pairwise voxel measures. In this approach, we first use the Least Absolute Shrinkage Selection Operator (LASSO) to pre-select voxels for analysis, then compute a Multivariate Vector AutoRegressive (MVAR) model from the time series of the selected voxels, and finally compute summary Granger Causality (GC) statistics from the model to represent directed interregional interactions. We demonstrate the effectiveness of this approach on both simulated and empirical fMRI data. We also show that averaging regional BOLD activity to create a nodal signal may lead to biased GC estimation of directed interregional interactions. The approach presented here makes it feasible to compute GC between brain regions without the need for averaging. Our results suggest that in the analysis of functional brain networks, careful consideration must be given to the way that network nodes and edges are defined because those definitions may have important implications for the validity of the analysis.     

Microelectrode arrays and application to biosensing devices

An electrochemical microanalytical system consisting of a microelectrode array, a micromachined flow-through assembly, and a multichannel potentiostat were constructed for highly sensitive biosensing. Thin-film platinum microelectrode arrays consisting of four interdigitated microelectrodes (IDAs), which are spaced in the sub-micrometer range, were fabricated using silicon technology. On top of this chip, a micromachined flow-through cell was mounted. Using a home made miniaturized multipotentiostat, amperometric measurements of the individual electrodes at different and changing potentials, respective to a single reference electrode, were performed simultaneously. The signal generation, signal processing and the analytical system were controlled by a computer (PC type) and special software. An improved sensor sensitivity was achieved by multielectrode detection and averaging of the IDA responses.

By applying both the oxidation and reduction potentials of reversible redox molecules to pairs of the interdigitated electrodes, an increased current generation can be observed. Thus the steady state current of mediators such as benzoquinone can be amplified by a factor of 30 compared with conventional electrodes. This measuring principle was applied to determine of the activity of hydrolases by detecting the enzyme generated p-aminophenol in the nanomolar range. By combining both, the averaging and the recycling procedures, the detection limit of amperometric biosensing devices may be lowered by about one and a half orders of magnitude.     

精准定位脑刺激对运动功能调控研究进展

脑刺激是神经科学研究的重要手段,传统的经颅磁刺激和经颅电刺激等脑刺激方法尽管能调控运动功能(包括减轻运动性障碍疾病的运动障碍、提高运动能力等),但存在空间分辨率低且无法刺激深部脑组织的局限性.近年来迅速发展的深部脑刺激(deep brain stimulation,DBS)、光遗传学、经颅超声刺激(transcranial ultrasound stimulation,TUS)、时间干涉(temporal interference,TI)等精准定位脑刺激方法,具有空间分辨率高、可聚焦深部脑组织等优点.本文综述了上述几种脑刺激方法的原理、特点,对运动功能调控的研究进展,以及面临的挑战和发展前景,从而为神经科学研究提供更好的研究工具,为临床实践提供更多的干预治疗手段.     

MIDA: A Multimodal Imaging-Based Detailed Anatomical Model of the Human Head and Neck

Computational modeling and simulations are increasingly being used to complement experimental testing for analysis of safety and efficacy of medical devices. Multiple voxel- and surface-based whole- and partial-body models have been proposed in the literature, typically with spatial resolution in the range of 1–2 mm and with 10–50 different tissue types resolved. We have developed a multimodal imaging-based detailed anatomical model of the human head and neck, named “MIDA”. The model was obtained by integrating three different magnetic resonance imaging (MRI) modalities, the parameters of which were tailored to enhance the signals of specific tissues: i) structural T1- and T2-weighted MRIs; a specific heavily T2-weighted MRI slab with high nerve contrast optimized to enhance the structures of the ear and eye; ii) magnetic resonance angiography (MRA) data to image the vasculature, and iii) diffusion tensor imaging (DTI) to obtain information on anisotropy and fiber orientation. The unique multimodal high-resolution approach allowed resolving 153 structures, including several distinct muscles, bones and skull layers, arteries and veins, nerves, as well as salivary glands. The model offers also a detailed characterization of eyes, ears, and deep brain structures. A special automatic atlas-based segmentation procedure was adopted to include a detailed map of the nuclei of the thalamus and midbrain into the head model. The suitability of the model to simulations involving different numerical methods, discretization approaches, as well as DTI-based tensorial electrical conductivity, was examined in a case-study, in which the electric field was generated by transcranial alternating current stimulation. The voxel- and the surface-based versions of the models are freely available to the scientific community.     

Adaptation is not required to explain the long-term response of axons to molecular gradients

It has been suggested that growth cones navigating through the developing nervous system might display adaptation, so that their response to gradient signals is conserved over wide variations in ligand concentration. Recently however, a new chemotaxis assay that allows the effect of gradient parameters on axonal trajectories to be finely varied has revealed a decline in gradient sensitivity on either side of an optimal concentration. We show that this behavior can be quantitatively reproduced with a computational model of axonal chemotaxis that does not employ explicit adaptation. Two crucial components of this model required to reproduce the observed sensitivity are spatial and temporal averaging. These can be interpreted as corresponding, respectively, to the spatial spread of signaling effects downstream from receptor binding, and to the finite time over which these signaling effects decay. For spatial averaging, the model predicts that an effective range of roughly one-third of the extent of the growth cone is optimal for detecting small gradient signals. For temporal decay, a timescale of about 3 minutes is required for the model to reproduce the experimentally observed sensitivity.     

Decoding of attentional state using local field potentials

We review recent efforts to decode visual spatial attention from different types of brain signals, such as spikes and local field potentials (LFPs). Combining signals from more electrodes improves decoding, but the pattern of improvement varies considerably depending on the signal as well as the task (for example, decoding of sensory stimulus/motor intention versus location of attention). We argue that this pattern of results conveys important information not only about the usefulness of a particular brain signal for decoding attention, but also about the spatial scale over which attention operates in the brain. The spatial scale, in turn, likely depends on the extent of underlying mechanisms such as normalization, gain control via excitation–inhibition interactions, and neuromodulatory regulation of attention.     

Simulation analysis of the ability to estimate motor unit propagation velocity non-invasively by different two-channel methods and types of multi-electrodes.

Ability to estimate motor unit propagation velocity correctly using different two-channel methods for delay estimation and different non-invasive spatial filters was analysed by simulation. It was established that longitudinal double difference electrodes could be not a better choice than simple bipolar parallel electrodes. Spatial filtration with a new multi-electrode (performing difference between signals detected by two transversal double difference electrodes positioned along the muscle fibres) promises to give the best estimate. Delay estimation between reference points is more preferable than that based on the cross-correlation technique, which is considerably sensitive to the fundamental properties of the muscle fibre extracellular fields. Preliminary averaging and approximation of the appropriate parts of the signals around chosen reference points could reduce the larger noise sensitivity and the effects of local tissue inhomogeneities as well as eliminate the sampling problem. A correct estimate of the propagation velocity could be impossible, even in the case of not very deep motor units (15 or 10 mm, depending on the spatial filter used) with relatively long (about 120 mm) muscle fibres. In the case of fibres with asymmetrical location of the end-plates in respect to the fibre ends, the propagation velocity estimates could be additionally biased above the longer semilength of the motor unit fibres.     

Adaptive Filtering Methods for Identifying Cross-Frequency Couplings in Human EEG

Oscillations have been increasingly recognized as a core property of neural responses that contribute to spontaneous, induced, and evoked activities within and between individual neurons and neural ensembles. They are considered as a prominent mechanism for information processing within and communication between brain areas. More recently, it has been proposed that interactions between periodic components at different frequencies, known as cross-frequency couplings, may support the integration of neuronal oscillations at different temporal and spatial scales. The present study details methods based on an adaptive frequency tracking approach that improve the quantification and statistical analysis of oscillatory components and cross-frequency couplings. This approach allows for time-varying instantaneous frequency, which is particularly important when measuring phase interactions between components. We compared this adaptive approach to traditional band-pass filters in their measurement of phase-amplitude and phase-phase cross-frequency couplings. Evaluations were performed with synthetic signals and EEG data recorded from healthy humans performing an illusory contour discrimination task. First, the synthetic signals in conjunction with Monte Carlo simulations highlighted two desirable features of the proposed algorithm vs. classical filter-bank approaches: resilience to broad-band noise and oscillatory interference. Second, the analyses with real EEG signals revealed statistically more robust effects (i.e. improved sensitivity) when using an adaptive frequency tracking framework, particularly when identifying phase-amplitude couplings. This was further confirmed after generating surrogate signals from the real EEG data. Adaptive frequency tracking appears to improve the measurements of cross-frequency couplings through precise extraction of neuronal oscillations.     

Few-sweep estimation of evoked potential based on a generalized subspace approach

A generalized subspace approach is proposed for single channel brain evoked potential (EP) extraction from background electroencephalogram (EEG) signal. The method realizes the optimum estimate of EP signal from the observable noisy signal. The underlying principle is to project the signal and noise into signal and noise coefficient subspace respectively by applying projection matrix at first. Secondly, coefficient weighting matrix is achieved based on the autocorrelation matrices of the noise and the noisy signal. With the coefficient weighting matrix, we can remove the noise projection coefficients and estimate the signal ones. EP signal is then obtained by averaging the signals estimated with the reconstruction matrix. Given different signal-to-noise ratio (SNR) conditions, the algorithm can estimate the EP signal with only two sweeps observable noisy signals. Our approach is shown to have excellent capability of estimating EP signal even in poor SNR conditions. The interference of spontaneous EEG has been eliminated with significantly improved SNR. The simulation results have demonstrated the effectiveness and superior performance of the proposed method.     

Vestibular-Somatosensory Interactions: Effects of Passive Whole-Body Rotation on Somatosensory Detection

Vestibular signals are strongly integrated with information from several other sensory modalities. For example, vestibular stimulation was reported to improve tactile detection. However, this improvement could reflect either a multimodal interaction or an indirect interaction driven by vestibular effects on spatial attention and orienting. Here we investigate whether natural vestibular activation induced by passive whole-body rotation influences tactile detection. In particular, we assessed the ability to detect faint tactile stimuli to the fingertips of the left and right hand during spatially congruent or incongruent rotations. We found that passive whole-body rotations significantly enhanced sensitivity to faint shocks, without affecting response bias. Critically, this enhancement of somatosensory sensitivity did not depend on the spatial congruency between the direction of rotation and the hand stimulated. Thus, our results support a multimodal interaction, likely in brain areas receiving both vestibular and somatosensory signals.     

Assessment of Subcortical Source Localization Using Deep Brain Activity Imaging Model with Minimum Norm Operators: A MEG Study

Subcortical structures are involved in many healthy and pathological brain processes. It is crucial for many studies to use magnetoencephalography (MEG) to assess the ability to detect subcortical generators. This study aims to assess the source localization accuracy and to compare the characteristics of three inverse operators in the specific case of subcortical generators. MEG has a low sensitivity to subcortical sources mainly because of their distance from sensors and their complex cyto-architecture. However, we show that using a realistic anatomical and electrophysiological model of deep brain activity (DBA), the sources make measurable contributions to MEG sensors signals. Furthermore, we study the point-spread and cross-talk functions of the wMNE, sLORETA and dSPM inverse operators to characterize distortions in cortical and subcortical regions and to study how noise-normalization methods can improve or bias accuracy. We then run Monte Carlo simulations with neocortical and subcortical activations. In the case of single hippocampus patch activations, the results indicate that MEG can indeed localize the generators in the head and the body of the hippocampus with good accuracy. We then tackle the question of simultaneous cortical and subcortical activations. wMNE can detect hippocampal activations that are embedded in cortical activations that have less than double their amplitude, but it does not completely correct the bias to more superficial sources. dSPM and sLORETA can still detect hippocampal activity above this threshold, but such detection might include the creation of ghost deeper sources. Finally, using the DBA model, we showed that the detection of weak thalamic modulations of ongoing brain activity is possible.     

Frequency sensitivity in Hodgkin–Huxley systems

 The frequency sensitivity of weak periodic signal detection has been studied via numerical simulations for both a single neuron and a neuronal network. The dependence of the critical amplitude of the signal upon its frequency and a resonance between the intrinsic oscillations of a neuron and the signal could account for the frequency sensitivity. In the presence of both a subthreshold periodic signal and noise, the signal-to-noise ratio (SNR) of the output of either a single neuron or a neuronal network present the typical characteristics of stochastic resonance. In particular, there exists a frequency-sensitive range of 30–100 Hz, and for signals with frequencies within this range the SNRs have large values. This implies that the system under consideration (a single neuron or a neuronal network) is more sensitive to the detection of periodic signals, and the frequency sensitivity may be of a functional significance to signal processing. Received: 26 October 1999 / Accepted in revised form: 25 July 2000     

Discrete-wavelet-transform-based noise removal and feature extraction for ECG signals

Nowadays, doctors use electrocardiogram (ECG) to diagnose heart diseases commonly. However, some nonideal effects are often distributed in ECG. Discrete wavelet transform (DWT) is efficient for nonstationary signal analysis. In this paper, the Symlets sym5 is chosen as the wavelet function to decompose recorded ECG signals for noise removal. Soft-thresholding method is then applied for feature detection. To detect ECG features, R peak of each heart beat is first detected, and the onset and offset of the QRS complex are then detected. Finally, the signal is reconstructed to remove high frequency interferences and applied with adaptive searching window and threshold to detect P and T waves. We use the MIT-BIH arrhythmia database for algorithm verification. For noise reduction, the SNR improvement is achieved at least 10 dB at SNR 5 dB, and most of the improvement SNR are better than other methods at least 1 dB at different SNR. When applying to the real portable ECG device, all R peaks can be detected when patients walk, run, or move at the speed below 9 km/h. The performance of delineation on database shows in our algorithm can achieve high sensitivity in detecting ECG features. The QRS detector attains a sensitivity over 99.94%, while detectors of P and T waves achieve 99.75% and 99.7%, respectively.     

Improvement of low-level light imaging performance using optical clearing method

Low-level light-emitting imaging technique often detects the light emerged at the tissue surface that is generated internally from a specific target. However, in most cases, the high scattering nature of biological tissue limits the sensitivity and spatial resolution of this imaging modality. In this paper, we report that a significant improvement of chemiluminescence (CL) imaging performance in terms of both sensitivity and spatial resolution can be achieved by use of the topical application of glycerol solution onto tissue sample, i.e. optical clearing approach. Monte Carlo (MC) simulation of internally-launched point source shows that the decrease of scattering coefficient of turbid medium, which can be achieved by optical tissue clearing approach, causes stronger peak intensity with a narrower full-width at half-maximum (FWHM). The improvement becomes more significant with the source depth increasing from 1 to 5 mm. The experimental results shows that tissue clearing with 50% glycerol solution could largely improve the brightness and the spatial resolution of CL imaging when the target is covered by biological tissue with a thickness of either 1 or 3mm. This method could have potential applications for the in vivo low-level light imaging techniques.     

NMR spectroscopy reveals unexpected structural variation at the protein–protein interface in MHC class I molecules

It is well established that non-uniform sampling (NUS) allows acquisition of multi-dimensional NMR spectra at a resolution that cannot be obtained with traditional uniform acquisition through the indirect dimensions. However, the impact of NUS on the signal-to-noise ratio (SNR) and sensitivity are less well documented. SNR and sensitivity are essential aspects of NMR experiments as they define the quality and extent of data that can be obtained. This is particularly important for spectroscopy with low concentration samples of biological macromolecules. There are different ways of defining the SNR depending on how to measure the noise, and the distinction between SNR and sensitivity is often not clear. While there are defined procedures for measuring sensitivity with high concentration NMR standards, such as sucrose, there is no clear or generally accepted definition of sensitivity when comparing different acquisition and processing methods for spectra of biological macromolecules with many weak signals close to the level of noise. Here we propose tools for estimating the SNR and sensitivity of NUS spectra with respect to sampling schedule and reconstruction method. We compare uniformly acquired spectra with NUS spectra obtained in the same total measuring time. The time saving obtained when only 1/k of the Nyquist grid points are sampled is used to measure k-fold more scans per increment. We show that judiciously chosen NUS schedules together with suitable reconstruction methods can yield a significant increase of the SNR within the same total measurement time. Furthermore, we propose to define the sensitivity as the probability to detect weak peaks and show that time-equivalent NUS can considerably increase this detection sensitivity. The sensitivity gain increases with the number of NUS indirect dimensions. Thus, well-chosen NUS schedules and reconstruction methods can significantly increase the information content of multidimensional NMR spectra of challenging biological macromolecules.     

Investigation of the source‐detector separation in near infrared spectroscopy for healthy and clinical applications

Understanding near infrared light propagation in tissue is vital for designing next generation optical brain imaging devices. Monte Carlo (MC) simulations provide a controlled mechanism to characterize and evaluate contributions of diverse near infrared spectroscopy (NIRS) sensor configurations and parameters. In this study, we developed a multilayer adult digital head model under both healthy and clinical settings and assessed light‐tissue interaction through MC simulations in terms of partial differential pathlength, mean total optical pathlength, diffuse reflectance, detector light intensity and spatial sensitivity profile of optical measurements. The model incorporated four layers: scalp, skull, cerebrospinal‐fluid and cerebral cortex with and without a customizable lesion for modeling hematoma of different sizes and depths. The effect of source‐detector separation (SDS) on optical measurements' sensitivity to brain tissue was investigated. Results from 1330 separate simulations [(4 lesion volumes × 4 lesion depths for clinical +3 healthy settings) × 7 SDS × 10 simulation = 1330)] each with 100 million photons indicated that selection of SDS is critical to acquire optimal measurements from the brain and recommended SDS to be 25 to 35 mm depending on the wavelengths to obtain optical monitoring of the adult brain function. The findings here can guide the design of future NIRS probes for functional neuroimaging and clinical diagnostic systems.      

Using Information Theory to Determine Optimum Pixel Size and Shape for Ecological Studies: Aggregating Land Surface Characteristics in Arctic Ecosystems

Quantifying vegetation structure and function is critical for modeling ecological processes, and an emerging challenge is to apply models at multiple spatial scales. Land surface heterogeneity is commonly characterized using rectangular pixels, whose length scale reflects that of remote sensing measurements or ecological models rather than the spatial scales at which vegetation structure and function varies. We investigated the ‘optimum’ pixel size and shape for averaging leaf area index (LAI) measurements in relatively large (85 m² estimates on a 600 × 600-m² grid) and small (0.04 m² measurements on a 40 × 40-m² grid) patches of sub-Arctic tundra near Abisko, Sweden. We define the optimum spatial averaging operator as that which preserves the information content (IC) of measured LAI, as quantified by the normalized Shannon entropy (E _S,n) and Kullback–Leibler divergence (D _KL), with the minimum number of pixels. Based on our criterion, networks of Voronoi polygons created from triangulated irregular networks conditioned on hydrologic and topographic indices are often superior to rectangular shapes for averaging LAI at some, frequently larger, spatial scales. In order to demonstrate the importance of information preservation when upscaling, we apply a simple, validated ecosystem carbon flux model at the landscape level before and after spatial averaging of land surface characteristics. Aggregation errors are minimal due to the approximately linear relationship between flux and LAI, but large errors of approximately 45% accrue if the normalized difference vegetation index (NDVI) is averaged without preserving IC before conversion to LAI due to the nonlinear NDVI-LAI transfer function. Author Contributions:  PS devised and undertook the analyses and wrote the paper. MW devised and implemented the measurement plan, and reviewed the analysis. LS assisted in the spatial data analysis and derivation of the terrain indices. RAB provided the macro-scale dataset. APB generated the DEM from aircraft data. JGE provided meteorological data. MvW generated the micro-scale field data with the help of Lorna Street and Sven Rasmussen. All authors contributed to the text.