Robustness of mutual information to inter-subject variability for automatic artefact removal from EEG

The externally recorded electroencephalogram (EEG) is contaminated with signals that do not originate from the brain, collectively known as artefacts. Thus, EEG signals must be cleaned prior to any further analysis. In particular, if the EEG is to be used in online applications such as Brain-Computer Interfaces (BCIs) the removal of artefacts must be performed in an automatic manner. This paper investigates the robustness of Mutual Information based features to inter-subject variability for use in an automatic artefact removal system. The system is based on the separation of EEG recordings into independent components using a temporal ICA method, RADICAL, and the utilisation of a Support Vector Machine for classification of the components into EEG and artefact signals. High accuracy and robustness to inter-subject variability is achieved.

Novel single trial movement classification based on temporal dynamics of EEG

Various complex oscillatory processes are involved in the generation of the motor command. The temporal dynamics of these processes were studied for movement detection from single trial electroencephalogram (EEG). Autocorrelation analysis was performed on the EEG signals to find robust markers of movement detection. The evolution of the autocorrelation function was characterised via the relaxation time of the autocorrelation by exponential curve fitting. It was observed that the decay constant of the exponential curve increased during movement, indicating that the autocorrelation function decays slowly during motor execution. Significant differences were observed between movement and no moment tasks. Additionally, a linear discriminant analysis (LDA) classifier was used to identify movement trials with a peak accuracy of 74%.

Towards natural human computer interaction in BCI

BCI systems require correct classification of signals interpreted from the brain for useful operation. To this end this paper investigates a method proposed in [1] to correctly classify a series of images presented to a group of subjects in [2]. We show that it is possible to use the proposed methods to correctly recognise the original stimuli presented to a subject from analysis of their EEG. Additionally we use a verification set to show that the trained classification method can be applied to a different set of data. We go on to investigate the issue of invariance in EEG signals. That is, the brain representation of similar stimuli is recognisable across different subjects. Finally we consider the usefulness of the methods investigated towards an improved BCI system and discuss how it could potentially lead to great improvements in the ease of use for the end user by offering an alternative, more intuitive control based mode of operation.

Toward construction of an inexpensive brain computer interface for goal oriented applications

The paper describes the implementation of an offline, low-cost Brain Computer Interface (BCI) alternative to more expensive commercial models. Using inexpensive general purpose clinical EEG acquisition hardware (Truscan32, Deymed Diagnostic) as the base unit, a synchronisation module was constructed to allow the EEG hardware to be operated precisely in time to allow for recording of automatically time stamped EEG signals. The synchronising module allows the EEG recordings to be aligned in stimulus time locked fashion for further processing by the classifier to establish the class of the stimulus, sample by sample. This allows for the acquisition of signals from the subject’s brain for the goal oriented BCI application based on the oddball paradigm. An appropriate graphical user interface (GUI) was constructed and implemented as the method to elicit the required responses (in this case Event Related Potentials or ERPs) from the subject.

Modulation of event-related desynchronization during kinematic and kinetic hand movements

Background Event-related desynchronization/synchronization (ERD/ERS) is a relative power decrease/increase of electroencephalogram (EEG) in a specific frequency band during physical motor execution and mental motor imagery, thus it is widely used for the brain-computer interface (BCI) purpose. However what the ERD really reflects and its frequency band specific role have not been agreed and are under investigation. Understanding the underlying mechanism which causes a significant ERD would be crucial to improve the reliability of the ERD-based BCI. We systematically investigated the relationship between conditions of actual repetitive hand movements and resulting ERD. Methods Eleven healthy young participants were asked to close/open their right hand repetitively at three different speeds (Hold, 1/3 Hz, and 1 Hz) and four distinct motor loads (0, 2, 10, and 15 kgf). In each condition, participants repeated 20 experimental trials, each of which consisted of rest (8–10 s), preparation (1 s) and task (6 s) periods. Under the Hold condition, participants were instructed to keep clenching their hand (i.e., isometric contraction) during the task period. Throughout the experiment, EEG signals were recorded from left and right motor areas for offline data analysis. We obtained time courses of EEG power spectrum to discuss the modulation of mu and beta-ERD/ERS due to the task conditions. Results We confirmed salient mu-ERD (8–13 Hz) and slightly weak beta-ERD (14–30 Hz) on both hemispheres during repetitive hand grasping movements. According to a 3 × 4 ANOVA (speed × motor load), both mu and beta-ERD during the task period were significantly weakened under the Hold condition, whereas no significant difference in the kinetics levels and interaction effect was observed. Conclusions This study investigates the effect of changes in kinematics and kinetics on resulting ERD during repetitive hand grasping movements. The experimental results suggest that the strength of ERD may reflect the time differentiation of hand postures in motor planning process or the variation of proprioception resulting from hand movements, rather than the motor command generated in the down stream, which recruits a group of motor neurons.

Rejoinder to 'Performance of different synchronisation measures in real data: A case study on electroencephalograhic signals

We agree with Duckrow and Albano [Phys. Rev. E 67, 063901 (2003)] and Quian Quiroga et al. [Phys. Rev. E 67, 063902 (2003)] that mutual information (MI) is a useful measure of dependence for electroencephalogram (EEG) data, but we show that the improvement seen in the performance of MI on extracting dependence trends from EEG is more dependent on the type of MI estimator rather than any embedding technique used. In an independent study we conducted in search for an optimal MI estimator, and in particular for EEG applications, we examined the performance of a number of MI estimators on the data set used by Quian Quiroga et al. in their original study, where the performance of different dependence measures on real data was investigated [Phys. Rev. E 65, 041903 (2002)]. We show that for EEG applications the best performance among the investigated estimators is achieved by k-nearest neighbors, which supports the conjecture by Quian Quiroga et al. in Phys. Rev. E 67, 063902 (2003) that the nearest neighbor estimator is the most precise method for estimating MI.

Comment on "Performance of different synchronization measures in real data: A case study on electroencephalographic signals"

We agree with Duckrow and Albano [Phys. Rev. E 67, 063901 (2003)] and Quian Quiroga [Phys. Rev. E 67, 063902 (2003)] that mutual information (MI) is a useful measure of dependence for electroencephalogram (EEG) data, but we show that the improvement seen in the performance of MI on extracting dependence trends from EEG is more dependent on the type of MI estimator rather than any embedding technique used. In an independent study we conducted in search for an optimal MI estimator, and in particular for EEG applications, we examined the performance of a number of MI estimators on the data set used by Quian Quiroga in their original study, where the performance of different dependence measures on real data was investigated [Phys. Rev. E 65, 041903 (2002)]. We show that for EEG applications the best performance among the investigated estimators is achieved by k-nearest neighbors, which supports the conjecture by Quian Quiroga in Phys. Rev. E 67, 063902 (2003) that the nearest neighbor estimator is the most precise method for estimating MI.

A novel approach to the detection of synchronisation in EEG based on empirical mode decomposition

Transient neural assemblies mediated by synchrony in particular frequency ranges are thought to underlie cognition. We propose a new approach to their detection, using empirical mode decomposition (EMD), a data-driven approach removing the need for arbitrary bandpass filter cut-offs. Phase locking is sought between modes. We explore the features of EMD, including making a quantitative assessment of its ability to preserve phase content of signals, and proceed to develop a statistical framework with which to assess synchrony episodes. Furthermore, we propose a new approach to ensure signal decomposition using EMD. We adapt the Hilbert spectrum to a time-frequency representation of phase locking and are able to locate synchrony successfully in time and frequency between synthetic signals reminiscent of EEG. We compare our approach, which we call EMD phase locking analysis (EMDPL) with existing methods and show it to offer improved time-frequency localisation of synchrony.

Perception and modification of plant flavonoid signals by rhizosphere microorganisms

Flavonoids are a diverse class of polyphenolic compounds that are produced as a result of plant secondary metabolism. They are known to play a multifunctional role in rhizospheric plant-microbe and plant-plant communication. Most familiar is their function as a signal in initiation of the legume-rhizobia symbiosis, but, flavonoids may also be signals in the establishment of arbuscular mycorrhizal symbiosis and are known agents in plant defence and in allelopathic interactions. Flavonoid perception by, and impact on, their microbial targets (e.g. rhizobia, plant pathogens) is relatively well characterized. However, potential impacts on 'non-target' rhizosphere inhabitants ('non-target' is used to distinguish those microorganisms not conventionally known as targets) have not been thoroughly investigated. Thus, this review first summarizes the conventional roles of flavonoids as nod gene inducers, phytoalexins and allelochemicals before exploring questions concerning 'non-target' impacts. We hypothesize that flavonoids act to shape rhizosphere microbial community structure because they represent a potential source of carbon and toxicity and that they impact on rhizosphere function, for example, by accelerating the biodegradation of xenobiotics. We also examine the reverse question, 'how do rhizosphere microbial communities impact on flavonoid signals?' The presence of microorganisms undoubtedly influences the quality and quantity of flavonoids present in the rhizosphere, both through modification of root exudation patterns and microbial catabolism of exudates. Microbial alteration and attenuation of flavonoid signals may have ecological consequences for below-ground plant-microbe and plant-plant interaction. We have a lack of knowledge concerning the composition, concentration and bioavailability of flavonoids actually experienced by microbes in an intact rhizosphere, but this may be addressed through advances in microspectroscopic and biosensor techniques. Through the use of plant mutants defective in flavonoid biosynthesis, we may also start to address the question of the significance of flavonoids in shaping rhizosphere community structure and function.

Life cycle of the QBO-modulated 11-year solar cycle signals in the Northern Hemispheric winter

This paper provides some insights on the quasi-biennial oscillation (QBO) modulated 11-year solar cycle (11-yr SC) signals in Northern Hemisphere (NH) winter temperature and zonal wind. Daily ERA-40 Reanalysis and ECMWF Operational data for the period of 1958-2006 were used to examine the seasonal evolution of the QBO-solar cycle relationship at various pressure levels up to the stratopause. The results show that the solar signals in the NH winter extratropics are indeed QBO-phase dependent, moving poleward and downward as winter progresses with a faster descent rate under westerly QBO than under easterly QBO. In the stratosphere, the signals are highly significant in late January to early March and have a life span of 30-50 days. Under westerly QBO, the stratospheric solar signals clearly lead and connect to those in the troposphere in late March and early April where they have a life span of 10 days. As the structure changes considerably from the upper stratosphere to the lower troposphere, the exact month when the maximum solar signals occur depends largely on the altitude chosen. For the low-latitude stratosphere, our analysis supports a vertical double-peaked structure of positive signature of the 11-yr SC in temperature, and demonstrates that this structure is further modulated by the QBO. These solar signals have a longer life span (3-4 months) in comparison to those in the extratropics. The solar signals in the lower stratosphere are stronger in early winter but weaker in late winter, while the reverse holds in the upper stratosphere.

Viral nucleolar localisation signals determine dynamic trafficking within the nucleolus

Localisation of both viral and cellular proteins to the nucleolus is determined by a variety of factors including nucleolar localisation signals (NoLSs), but how these signals operate is not clearly understood. The nucleolar trafficking of wild type viral proteins and chimeric proteins, which contain altered NoLSs, were compared to investigate the role of NoLSs in dynamic nucleolar trafficking. Three viral proteins from diverse viruses were selected which localised to the nucleolus; the coronavirus infectious bronchitis virus nucleocapsid (N) protein, the herpesvirus saimiri ORF57 protein and the HIV-1 Rev protein. The chimeric proteins were N protein and ORF57 protein which had their own NoLS replaced with those from ORF57 and Rev proteins, respectively. By analysing the sub-cellular localisation and trafficking of these viral proteins and their chimeras within and between nucleoli using confocal microscopy and photo-bleaching we show that NoLSs are responsible for different nucleolar localisations and trafficking rates.

Design and analysis of climate model experiments for the efficient estimation of anthropogenic signals

Presented herein is an experimental design that allows the effects of several radiative forcing factors on climate to be estimated as precisely as possible from a limited suite of atmosphere-only general circulation model (GCM) integrations. The forcings include the combined effect of observed changes in sea surface temperatures, sea ice extent, stratospheric (volcanic) aerosols, and solar output, plus the individual effects of several anthropogenic forcings. A single linear statistical model is used to estimate the forcing effects, each of which is represented by its global mean radiative forcing. The strong colinearity in time between the various anthropogenic forcings provides a technical problem that is overcome through the design of the experiment. This design uses every combination of anthropogenic forcing rather than having a few highly replicated ensembles, which is more commonly used in climate studies. Not only is this design highly efficient for a given number of integrations, but it also allows the estimation of (nonadditive) interactions between pairs of anthropogenic forcings. The simulated land surface air temperature changes since 1871 have been analyzed. The changes in natural and oceanic forcing, which itself contains some forcing from anthropogenic and natural influences, have the most influence. For the global mean, increasing greenhouse gases and the indirect aerosol effect had the largest anthropogenic effects. It was also found that an interaction between these two anthropogenic effects in the atmosphere-only GCM exists. This interaction is similar in magnitude to the individual effects of changing tropospheric and stratospheric ozone concentrations or to the direct (sulfate) aerosol effect. Various diagnostics are used to evaluate the fit of the statistical model. For the global mean, this shows that the land temperature response is proportional to the global mean radiative forcing, reinforcing the use of radiative forcing as a measure of climate change. The diagnostic tests also show that the linear model was suitable for analyses of land surface air temperature at each GCM grid point. Therefore, the linear model provides precise estimates of the space time signals for all forcing factors under consideration. For simulated 50-hPa temperatures, results show that tropospheric ozone increases have contributed to stratospheric cooling over the twentieth century almost as much as changes in well-mixed greenhouse gases.