Brain-computer interfaces for neurorehabilitation.

Brain-computer interfaces (BCIs) enable control of computers and other assistive devices, such as neuro-prostheses, which are used for communication, movement restoration, neuro-modulation, and muscle stimulation, by using only signals measured directly from the brain. A BCI creates a new output channel for the brain to a computer or a device. This requires retrieval of signals of interest from the brain, and its use for neuro-rehabilitation by means of interfacing the signals to a computerized device. Brain signals such as action potentials from single neurons or nerve fibers, extracellular local field potentials (LFPs), electrocorticograms, electroencephalogram and its components such as the event-related brain potentials, real-time functional magnetic resonance imaging, near-infrared spectroscopy, and magneto-encephalogram have been used. BCIs are envisaged to be useful for communication, control and self-regulation of brain function. BCIs employ neurofeedback to enable operant conditioning to allow the user to learn using it. Paralytic conditions arising from stroke or other diseases are being targeted for BCI application. Neurofeedback strategies ranging from sensory feedback to direct brain stimulation are being employed. Existing BCIs are limited in their throughput in terms of letters per minute or commands per minute, and need extensive training to use the BCI. Further, they can cause rapid fatigue due to use and have limited adaptability to changes in the patient's brain state. The challenge before BCI technology for neuro-rehabilitation today is to enable effective clinical use of BCIs with minimal effort to set up and operate.

Acousto-optic modulation by pulsed optical excitation: implications to imaging in turbid media.

We show that the transient response of acoustically modulated optical flux in a turbid medium irradiated by a pulsed point source of light is delayed in time relative to the light-alone flux obtained in the absence of acoustic modulation. The time delay is shown to result from an initial phase of flux reversal, as determined by the time point of the input pulse onset with reference to the ultrasound cycle. Both the time delay and amplitude of modulation are shown to be dependent on the effective attenuation coefficient of the medium. Application of a periodic train of excitation pulses spaced at equal intervals at, or in multiples of, the ultrasound period enables a time-locked detection of the modulated light, without the deleterious effects caused by speckle artifacts.

Anisotropic scattering of discrete particle arrays.

Far-field intensities of light scattered from a linear centro-symmetric array illuminated by a plane wave of incident light are estimated at a series of detector angles. The intensities are computed from the superposition of E-fields scattered by the individual array elements. An average scattering phase function is used to model the scattered fields of individual array elements. The nature of scattering from the array is investigated using an image (theta-phi plot) of the far-field intensities computed at a series of locations obtained by rotating the detector angle from 0 degrees to 360 degrees, corresponding to each angle of incidence in the interval [0 degrees 360 degrees]. The diffraction patterns observed from the theta-Phi plot are compared with those for isotropic scattering. In the absence of prior information on the array geometry, the intensities corresponding to theta-Phi pairs satisfying the Bragg condition are used to estimate the phase function. An algorithmic procedure is presented for this purpose and tested using synthetic data. The relative error between estimated and theoretical values of the phase function is shown to be determined by the mean spacing factor, the number of elements, and the far-field distance. An empirical relationship is presented to calculate the optimal far-field distance for a given specification of the percentage error.

Acousto-optic modulation of a point-scatterer array.

Acoustic modulation of light scattering from a linear centrosymmetric array is analyzed by considering far-field contributions due to optoelastic (OE) effect and acoustically induced translation of the array elements. The modulated light intensity is shown to vary sinusoidally at the acoustic frequency when the physical constants representative of the above effects are within ranges of their physical limits. The OE and translation components of the acousto-optic (AO) signal are shown to be in phase quadrature, each exhibiting a double-sided maxima when expressed as a function of the detector angle.

On the eigenfunctions of acousto-optic modulation in a homogeneously absorbing optical medium.

A closed form solution for time-averaged modulated fluence rate is presented for acoustically modulated diffusive light propagation in a medium. The solution assumes that the component of modulated light flux in the direction of acoustic pressure variation is zero.

New insights into image processing of cortical blood flow monitors using laser speckle imaging.

Laser speckle imaging has increasingly become a viable technique for real-time medical imaging. However, the computational intricacies and the viewing experience involved limit its usefulness for real-time monitors such as those intended for neurosurgical applications. In this paper, we propose a new technique, tLASCA, which processes statistics primarily in the temporal direction using the laser speckle contrast analysis (LASCA) equation, proposed by Briers and Webster. This technique is thoroughly compared with the existing techniques for signal processing of laser speckle images, including, the spatial-based sLASCA and the temporal-based modified laser speckle imaging (mLSI) techniques. sLASCA is an improvement of the basic LASCA technique. In sLASCA, the derived contrasts are further averaged over a predetermined number of raw speckle images. mLSI, on the other hand, is the technique in which temporal statistics are processed using the equation described by Ohtsubo and Asakura. tLASCA preserves the original image resolution similar to mLSI. tLASCA outperforms sLASCA (window size M = 5) with faster convergence of K values (5.32 versus 20.56 s), shorter per-frame processing time (0.34 versus 2.51 s), and better subjective and objective quality evaluations of contrast images. tLASCA also outperforms mLSI with faster convergence of K values (5.32 s) compared to N values (10.44 s), shorter per-frame processing time (0.34 versus 0.91 s), smaller intensity fluctuations among frames (8%-10% versus 15%-35%), and better subjective and objective quality evaluations of contrast images. As laser speckle imaging becomes an important tool for real-time monitoring of blood flows and vascular perfusion, tLASCA is proven to be the technique of choice.

Statistical mapping of speckle autocorrelation for visualization of hyperaemic responses to cortical stimulation.

Statistically mapped speckle autocorrelation images (SAR) were used to track the hemodynamically active perfusion regions in the rat cortex during and following DC current stimulation with high transverse spatial resolution (38 um). The SAR images provided a spatio-temporal information about the net activation patterns of Cerebral Blood Flow (CBF) changes over a period of time as against those changes for each frame interval estimated using spatial contrasts derived from the first order spatial statistics. Thus the information about the relative maxima of perfusion during a Transient Hyperaemic Episode (THE) across different regions in the imaging window could be identified without the need for actually having to estimate the spatial contrast maps of the imaged region for each frame contained in the time window of observation. With the application of DC stimulation, the regions with a high correlation in the temporal fluctuations were representative of the areas that underwent least changes in activation. By varying the intensity of stimulation, THEs were observed for stimulation current densities in the range 0.1-3.8 mA/mm2 using both the derived speckle contrast maps and concurrently on a Laser Doppler Flow meter, with its probe positioned 1 mm from the site of stimulation. For current densities below the lower threshold of stimulation, the SAR images revealed an unprecedented reduction in the surge amplitude at sites distal to the region of stimulation. This was accompanied by an increase in pixel areas representing minimally active regions of perfusion ("perfusion islets") with no identifiable peak in the hemodynamic responses estimated from speckle contrast variations. The SAR images can be a useful tool for visualization of slow wave perfusion dynamics during cortical stimulation.

Where the BOLD signal goes when alpha EEG leaves.

Previous studies using simultaneous EEG and fMRI recordings have yielded discrepant results regarding the topography of brain activity in relation to spontaneous power fluctuations in the alpha band of the EEG during eyes-closed rest. Here, we explore several possible explanations for this discrepancy by re-analyzing in detail our previously reported data. Using single subject analyses as a starting point, we found that alpha power decreases are associated with fMRI signal increases that mostly follow two distinct patterns: either 'visual' areas in the occipital lobe or 'attentional' areas in the frontal and parietal lobe. On examination of the EEG spectra corresponding to these two fMRI patterns, we found greater relative theta power in sessions yielding the 'visual' fMRI pattern during alpha desynchronization and greater relative beta power in sessions yielding the 'attentional' fMRI pattern. The few sessions that fell into neither pattern featured the overall lowest theta and highest beta power. We conclude that the pattern of brain activation observed during spontaneous power reduction in the alpha band depends on the general level of brain activity as indexed over a broader spectral range in the EEG. Finally, we relate these findings to the concepts of 'resting state' and 'default mode' and discuss how - as for sleep - EEG-based criteria might be used for staging brain activity during wakefulness.

Early adaptations in somatosensory cortex after focal ischemic injury to motor cortex.

In response to a lesion, intact regions of cortex in both hemispheres undergo adaptive changes in network function. For example, changes in excitability and intracortical inhibition in primary motor cortex (M1) were reported after lesioning contralateral or ipsilateral brain regions. Close interactions exist between M1 and primary somatosensory cortex (S1) within one hemisphere. Therefore, we hypothesized that lasting modifications would occur in S1 excitability after lesioning ipsilateral M1. Imaging of intrinsic optical signals (IOS, at 570 nm) was used to investigate the evolution of the somatosensory cortical response evoked by contralateral median nerve stimulation during the first hour after a photothrombotic lesion to M1 (caudal motor cortex) of the rat (n=10). Control rats (n=6) received no lesion. Perfusion was monitored by Laser speckle imaging and the extent of the resulting lesion was determined histologically. Control animals did not show evidence for reduced perfusion, infarction, or changes in IOS. M1 infarction led to a significant increase in evoked response amplitude, duration, and area of activation, and a shortening of latencies. These parameters reached a plateau around 50 min after ischemia. These results indicate S1 hyperexcitability after M1 injury. Whether these adaptations contribute to functional deficits or play a role in recovery, remains to be determined.

Imaging the development of an ischemic core following photochemically induced cortical infarction in rats using Laser Speckle Contrast Analysis (LASCA).

Laser Speckle Contrast Analysis (LASCA) has been used to image the development of an ischemic core following photochemically induced infarction in rats up to 1 h post-lesion. Using LASCA, we have been able to image a central ischemic core which had little or no perfusion surrounded by a penumbral region with reduced perfusion. In addition, we have shown the existence of a surrounding region of hyperemic tissue. A potential feature of this imaging approach is its capability to track cerebral blood flow (CBF) changes in the region within and outside the ischemic core besides demonstrating the real-time progression of the ischemic core into the penumbral region. We have demonstrated the continuous disruption of CBF to the ischemic core that eventually affected the blood supply to the surrounding regions. The penumbral flow is shown to exhibit a sudden increase post-ischemic induction followed by a slow decline to the final baseline level. Interestingly, we observed an interaction (P < 0.03) between penumbral flow peak and the time effects of increase in pixel area from the infarct region to the surrounding penumbral region. Using a paired-sample t test, we observed that the mean pixel area was larger for the infarct region than for the penumbral region (P < 0.004) during the time interval between the induction of ischemia and the time point of peak flow in the penumbral region.

Wavelet entropy for subband segmentation of EEG during injury and recovery.

In this paper, subband wavelet entropy (SWE) is used for the segmentation of electroencephalographic signals (EEG) recorded during injury and recovery following global cerebral ischemia. Wavelet analysis is used to decompose the EEG into standard clinical subbands followed by computation of the Shannon entropy. The EEG was measured from rodent brains in a controlled experimental brain injury model by hypoxic-ischemic cardiac arrest. Results show that while the relative EEG power failed to reveal the order of bursting activity associated with recovery, SWE was used to segment the EEG and delineate the initial bursting periods in each subband. Based on entropy variations obtained from a cohort of animals with graded levels of hypoxic-ischemic cardiac arrest, an intermittent pattern of bursting was observed in the high frequency bands.

Prediction of PTZ-induced seizures using wavelet-based residual entropy of cortical and subcortical field potentials.

Our proposed algorithm for seizure prediction is based on the principle that seizure build-up is always preceded by constantly changing bursting levels. We use a novel measure of residual subband wavelet entropy (RSWE) to directly estimate the entropy of bursts, which is otherwise obscured by the ongoing background activity. Our results are obtained using a slow infusion anesthetized pentylenetetrazol (PTZ) rat model in which we record field potentials (FPs) from frontal cortex and two thalamic areas (anterior and posterior nuclei). In each frequency band, except for the theta-delta frequency bands, we observed a significant build-up of RSWE from the preictal period to the first ictal event (p < or = 0.05) in cortex. Significant differences were observed between cortical and thalamic RSWE (p < or = 0.05) subsequent to seizure development. A key observation is the twofold increase in mean cortical RSWE from the preictal to interictal period. Exploiting this increase, we develop a slope change detector to discern early acceleration of entropy and predict the approaching seizure. We use multiple observations through sequential detection of slope changes to enhance the sensitivity of our prediction. Using the proposed method applied to a cohort of four rats subjected to PTZ infusion, we were able to predict the first seizure episode 28 min prior to its occurrence.

Early restitution of electrocorticogram predicts subsequent behavioral recovery from cardiac arrest.

Previous studies have shown that parameters of EEG restitution reflect the severity of global hypoxic-ischemic brain injury. Here, the hypothesis is tested that patterns of EEG restitution during the first 4 hours predict later behavioral recovery. Time course and correlations between behavior, electrocorticogram (EcoG), and neuronal injury were investigated in a rodent model of asphyctic cardiac arrest. Forty Wistar rats were subjected to 5 minutes of asphyxia and cardiopulmonary resuscitation. Behavior was assessed by repeated scoring of neurodeficits and open field activity until euthanasia at 48 hours. Electrocorticographic bursting occurred at 13.2 +/- 4 minutes after resuscitation. Bursts increased in frequency and duration until the EcoG reverted to a continuous signal. The resuscitation-continuous EcoG interval correlated with the first appearance of spontaneous movements (r = 0.80, P < 0.05). Larger intervals were associated with hyperactivity in the open field at 24 hours (r = 0.61, P < 0.05), indicating a more severe behavioral deficit. Larger intervals were also associated with worse 48-hour neurodeficit scores (P < 0.05). Neuronal damage in the hippocampus correlated with the degree of open field hyperactivity at 14 hours (P < 0.05). These findings demonstrate a close temporal and prognostic relationship between electrical and behavioral recovery after hypoxic-ischemic brain injury.