Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Prediction of Hand Trajectory from Electrocorticography Signals in Primary Motor Cortex

  • Chao Chen ,

    Contributed equally to this work with: Chao Chen, Duk Shin, Hidenori Watanabe

    Affiliation Department of Information Processing, Tokyo Institute of Technology, Yokohama, Japan

  • Duk Shin ,

    Contributed equally to this work with: Chao Chen, Duk Shin, Hidenori Watanabe

    shinduk@cns.pi.titech.ac.jp

    Affiliation Precision and Intelligence Laboratory, Tokyo Institute of Technology, Yokohama, Japan

  • Hidenori Watanabe ,

    Contributed equally to this work with: Chao Chen, Duk Shin, Hidenori Watanabe

    Affiliation Department of Developmental Physiology, National Institute for Physiological Sciences, National Institutes of Natural Sciences, Okazaki, Japan

  • Yasuhiko Nakanishi,

    Affiliation Precision and Intelligence Laboratory, Tokyo Institute of Technology, Yokohama, Japan

  • Hiroyuki Kambara,

    Affiliation Precision and Intelligence Laboratory, Tokyo Institute of Technology, Yokohama, Japan

  • Natsue Yoshimura,

    Affiliation Precision and Intelligence Laboratory, Tokyo Institute of Technology, Yokohama, Japan

  • Atsushi Nambu,

    Affiliation Graduate University for Advanced Studies (SOKENDAI), Hayama, Japan

  • Tadashi Isa,

    Affiliations Department of Developmental Physiology, National Institute for Physiological Sciences, National Institutes of Natural Sciences, Okazaki, Japan, Graduate University for Advanced Studies (SOKENDAI), Hayama, Japan

  • Yukio Nishimura,

    Affiliations Department of Developmental Physiology, National Institute for Physiological Sciences, National Institutes of Natural Sciences, Okazaki, Japan, Graduate University for Advanced Studies (SOKENDAI), Hayama, Japan, Precursory Research for Embryonic Science and Technology, Japan Science and Technology Agency, Tokyo, Japan

  • Yasuharu Koike

    Affiliations Department of Information Processing, Tokyo Institute of Technology, Yokohama, Japan, Precision and Intelligence Laboratory, Tokyo Institute of Technology, Yokohama, Japan, CREST, Japan Science and Technology Agency, Kawaguchi, Japan

Abstract

Due to their potential as a control modality in brain-machine interfaces, electrocorticography (ECoG) has received much focus in recent years. Studies using ECoG have come out with success in such endeavors as classification of arm movements and natural grasp types, regression of arm trajectories in two and three dimensions, estimation of muscle activity time series and so on. However, there still remains considerable work to be done before a high performance ECoG-based neural prosthetic can be realized. In this study, we proposed an algorithm to decode hand trajectory from 15 and 32 channel ECoG signals recorded from primary motor cortex (M1) in two primates. To determine the most effective areas for prediction, we applied two electrode selection methods, one based on position relative to the central sulcus (CS) and another based on the electrodes' individual prediction performance. The best coefficients of determination for decoding hand trajectory in the two monkeys were 0.4815±0.0167 and 0.7780±0.0164. Performance results from individual ECoG electrodes showed that those with higher performance were concentrated at the lateral areas and areas close to the CS. The results of prediction according with different numbers of electrodes based on proposed methods were also shown and discussed. These results also suggest that superior decoding performance can be achieved from a group of effective ECoG signals rather than an entire ECoG array.

Introduction

Over the past two decades, brain-machine interfaces (BMI) have been developed utilizing the growing understanding of brain function and the development of technology to measure brain activity. BMIs translate brain signals into commands for controlling devices such as cursors [1], spelling devices [2], robot arms, and neural prosthetics [3][6]. This new communication pathway has not only the potential to help to disabled persons but also provide insight into the motor system of the brain [7][13]. A number of methods have been developed to measure brain signals. BMIs are mainly categorized into two types, invasive and non-invasive BMIs, according to the signal source. BMI systems have been developed using modalities such as multi-neuron activity [14], [15], local field potentials [16], [17], electroencephalography [1], [7], [18], [19], and functional magnetic resonance imaging [20].

Electrocorticography (ECoG) has been in focus as a less invasive recording method for BMIs [21][36] since the first ECoG-based BMI succeeded in one-dimensional cursor control in human subjects [21]. ECoG signals have higher signal-to-noise ratio and spatiotemporal resolution than non-invasive recording methods, because ECoG electrodes are laid on the surface of the cerebral cortex. ECoG recording has also been shown to have long-term stability [22], [23], and its level of clinical risk is lower compared with invasive methods, because the electrodes do not penetrate the brain. Classifications of arm movement direction [24], [25], 3D cursor control [26], natural grasp type [27], [28], and hand posture [29], [30] have been achieved by using ECoG signals. Two-dimensional [31][33] and three-dimensional (3D) hand trajectories [22], [23] and muscle activities [34] have been decoded using epidural or subdural ECoG signals in time series. Despite these successes, however, which locations are most effective for ECoG-based hand trajectory prediction and how different numbers of effective ECoG signals affect decoding performance are still open questions.

In this study, and in investigation of these questions, we attempted to decode hand trajectory from ECoG signals. We recorded 15 and 32 ECoG signals of the primary motor cortex (Ml) and 3D hand positioning simultaneously in two Japanese monkeys while they performed reaching and grasping tasks. We predicted 3D hand trajectories using our previous signal preprocessing method [34] and a partial least squares (PLS) method. Two methods for electrode selection were proposed in order to examine the questions previously mentioned. Prediction performances with different combination of electrodes using the proposed decoding methods were compared. Both methods showed equivalent ability to predict hand trajectories. Our results indicated that 3D hand trajectories can be predicted using nine or ten ECoG signals and that ECoG electrodes with higher performance were concentrated at the lateral areas and areas close to the central sulcus (CS).

Methods

Ethics statement

All experimental procedures were performed in accordance with the Guidelines for Proper Conduct of Animal Experiments of the Science Council of Japan and approved by the Committee for Animal Experiments at the National Institutes of Natural Sciences (Approval No.: 11A157). The animals' welfare and steps taken to ameliorate suffering were in accordance with the recommendations of the Weatherall report, “The use of non-human primates in research”. The animals were monitored closely, and their welfare was assessed on a daily basis, or several times a day if necessary. This included veterinary examinations to ensure that they were not suffering, as well as the use of analgesics, antiemetics, or antibiotic therapy if necessary. The animals were housed individually on a 12-hour light/dark cycle and provided a rubber toy and ample food and water in their home cage. No animals were sacrificed in this study.

Behavioral Task

Two Japanese macaques (Monkey A: male, 8.9 kg; Monkey B: female, 4.7 kg) were trained to perform right hand reaching, grasping, pulling, and releasing tasks as shown in Figure 1A. The monkeys performed these tasks repeatedly and continuously for over 700 s. Monkey A performed a total of 134 trials, and monkey B performed 248 trials.

thumbnail
Figure 1. Behavioral task and location of ECoG electrodes used in decoding.

A) Monkeys performed right hand reaching, grasping, pulling, and releasing tasks in a 3D workspace. During the task, ECoG and hand positioning were recorded simultaneously. B) The planar-surface platinum electrode arrays were implanted on the gyrus between the central sulcus (CS) and the arcuate sulcus (AS) in the primary motor area in left hemisphere. The locations of all 15 and 32 electrodes in monkey A and monkey B are shown with defined channel numbers. Locations of 3, 6, 9, 12, and 15 or 16 electrode groups used in decoding are denoted with green solid, purple dotted, blue dotted, brown dotted, and blue solid lines, respectively. The column nearest the CS was column #1 in the rostral-caudal direction, and the row in the medial-lateral direction was row *1. Note that electrodes inside the blue line were used in both column and row decoding.

https://doi.org/10.1371/journal.pone.0083534.g001

Both monkeys underwent surgery to implant an ECoG electrode array under anesthesia after they completed behavioral training. We chronically implanted a platinum ECoG array (Unique Medical Corporation, Tokyo, Japan) over the left M1, which contained 15 (monkey A: 5×3 grid) and 32 (monkey B: 4×8 grid) channel electrodes, as shown in Figure 1B. Electrode locations were identified from anatomical views during surgery, preparation of brain and postoperative x-ray images. The midline of the brain was estimated from the sagittal suture and used to landmark electrode locations in the medial-lateral direction. The center electrode 8 in monkey A was approximately 15 mm from the midline. The medial electrodes in monkey B were located approximately 14 mm from the midline. The ECoG arrays were nearly parallel to the midline. In the rostral-caudal direction, for both monkeys A and B, centers of electrodes in column #1 (monkey A:3, 6, 9, 12, 15; monkey B: 7, 15, 23) were placed 1–2 mm rostral of the central sulcus. Electrode 31 in monkey B was placed 1–2 mm caudal of the central sulcus. Descriptions of the technical and surgical details can be found in our previous work [34].

Data recording

ECoG signals were sampled at 4 kHz using an acquisition processor system (Plexon MAP System; Plexon, Inc., Dallas, US). ECoG signals were filtered with band-pass filters through multi-channel bio-signal amplifiers (monkey A: 1.5 Hz high-pass and 1 kHz low-pass analog filters, MEG-6116, Nihon Kohden Corporation, Tokyo, Japan; monkey B: 0.7 Hz high-pass and 8 kHz low-pass analog filters, Plexon, Inc., Dallas, USA).

3D-positions of various points of the right arm were recorded using reflective markers tracked with an optical motion capture system (Eagle Digital System; Motion Analysis Corporation, Santa Rosa, CA). The system used twelve infrared cameras operating at 200 frames/s to track the positions of multiple reflective markers (4-mm-diameter spheroids). A total of fourteen markers were attached to the right arm of each monkey but we used only the wrist marker to extract hand positioning. In addition to optical data, the motion capture system also recorded analog signals from the external stimulator (SEN-8203; Nihon Kohden Corporation, Tokyo, Japan) for synchronization with the neural recordings. The neural data were down-sampled to 500 samples per second, and the motion data were up-sampled to 500 samples per second to match the neural data, similar in manner to our previous work [35].

Preprocessing and feature selection

Raw ECoG signals were re-referenced to a common average reference (CAR) to increase the signal-to-noise ratio in the preprocessing phase. The CAR method calculates the mean of all channels, and subtracts this value from the selected output channels [36], [37].

Nine specific frequency bands were selected for further analysis: δ (1.5∼4 Hz), θ (4∼8 Hz), α (8∼14 Hz), β1 (14∼20 Hz), β2 (20∼30 Hz), γ1 (30∼50 Hz), γ2 (50∼90 Hz), γ3 (90∼120 Hz), and γ4 (120∼150 Hz). These specific bands were selected due to their correlation with motor activity, as shown in previous ECoG-based BMI studies [27][34]. Band-pass filters for each of the nine frequency bands were used to transform the re-referenced ECoG signals into nine separate time series. Then, each time series was rectified and smoothed with a Gaussian filter of 0.1 s width (σ: 0.04 s). Finally, the smoothed time series at time t were z-score normalized to produce the final ECoG source signal as follows:(1)

where, i and j are the electrode channel and the frequency band, respectively. and denote the mean value and the standard deviation of over a 2 s interval before time t, respectively. These became the final ECoG feature signals for use in hand trajectory prediction. An example ECoG feature signal during a trial movement is shown in Figure 2.

thumbnail
Figure 2. Example of measured trajectory and frequency band feature data during a movement task.

Frequency band feature data were sorted into channels and frequency bands, as shown at the top. The X, Y, and Z positioning data recorded from the markers attached to the hand of monkey B, are shown at the bottom.

https://doi.org/10.1371/journal.pone.0083534.g002

Partial Least Squares Regression

Partial least squares regression (PLS) was used to decode the 3D hand positioning from ECoG. Because of its utility in variable selection and dimension reduction, PLS has been widely used in the fields of brain imaging, computational chemistry, data mining, and others [22], [23], [38][41].

The 3D hand positioning at time t, , was decoded using the ECoG feature signal over a 0.6 s interval before time t and can be described as (2)

where, p represents the predicted value of each xyz-coordinate, is 30 ms, are the weights according to the ECoG feature signal at electrode channel i, frequency band j, and time , and is the bias.

The PLS methods calculates a set of orthogonal factors called latent variables to model the relationship between two sets of data. Ten-fold cross validation was used to evaluate prediction by the model. To avoid over-fitting, the predictive error sum of squares (PRESS) was calculated to find the optimal number of latent variables in the PLS model, which can be described as(3)

where is the predicted hand position, and is the observed hand position.

Two methods for electrode selection

To investigate which electrode locations were more effective, we decreased the number of electrodes for prediction using two methods and compared their respective performance.

In the first method, electrodes were selected based on their implantation position. Previous physiological studies have shown that cortico-motoneuronal cells that encode muscle-activation patterns reflected in EMG activity are located predominantly in the anterior bank of the central sulcus (CS) [42], [43]. Our previous work [34] also showed that the area close to the CS might be key to decoding muscle activity. Therefore, we selected electrodes in groups of 3, 6, 9, 12, and 15 or 16 electrodes, expanding in distance from the CS as shown in Figure 1B. We refer to this method hereafter as location-based selection.

For the second selection method, electrodes were chosen based on prediction performance. Performance values for the PLS model using only one electrode were calculated and sorted by their coefficients of determination (R2). Then, electrodes with high performance were added in turn to train a new PLS model. To investigate the effective frequency band for prediction, performance values for the PLS model using only δ (1.5∼4 Hz), γ3 (90∼120 Hz), and γ4 (120∼150 Hz) bands were also calculated for each electrode. We refer to this method hereafter as performance-based selection.

Analysis

The entire 700 s of experiment data were divided into two parts, 500 s of training data and 200 s of test data. Ten-fold cross validation was employed to train the PLS model on the 500 s of training data. Then, the 200 s of test data were used to evaluate the PLS model.

PRESS values were calculated to find the optimal number of latent variables in the PLS model. Smaller PRESS values were associated with greater PLS model performance. Typically, the PRESS value decreases when effective latent variables are added to train the model. Then, if over-fitting occurs, the PRESS value increases. A good choice is to stop adding latent variables as soon as the PRESS value increases. In this study, however, the PRESS value decreased quickly when the number of latent variables was within 20, but then plateaued soon after, as shown in Figure 3 and Figure S1. Thus, we selected 20 as the optimal number of latent variables.

thumbnail
Figure 3. Predictive error sum of squares in model training for monkey B.

The blue line and green line show predictive error sum of squares (PRESS) and R2 values, respectively, for different numbers of latent variables used in the PLS model. The optimal number of 20 is denoted with the red dotted line.

https://doi.org/10.1371/journal.pone.0083534.g003

Weights of the prediction model were analyzed to evaluate the contribution of each the nine frequency bands used in this study. The contribution of frequency band Confb was calculated as (4)

Where, are the weights associated with the ECoG feature signal at electrode i, frequency band j, and time .

In addition, 3D hand trajectories were predicted using each of the nine frequency bands of the ECoG feature signals to investigate their individual contributions to prediction.

Result

Prediction with the location-based selection method

3D hand trajectories were first decoded using all the ECoG electrodes. For monkey A, the mean R2 value and standard deviation (STD) after 10-fold cross validation were 0.4840±0.0118, and mean R2 using the test data was 0.4806. For monkey B, the mean R2 values after 10-fold cross validation and using the test data were 0.8424±0.0032 (Mean±STD) and 0.7328, respectively.

We verified how decoding performance changes depending on the number of effective ECoG signals. Positions for the groups of 3, 6, 9, 12, and 15 or 16 ECoG electrodes selected to decode hand trajectories are shown in Figure 1B. For monkey A, R2 values for X, Y, and Z positioning were 0.4724, 0.4695, and 0.4997, respectively, obtained using all 15 electrodes. For monkey B with all 32 electrodes, R2 values for X, Y, and Z positioning were 0.7126, 0.7644 and 0.7263, respectively. One example of continuous prediction is shown in Figure 4 (see also Fig. S2 for monkey A). Actual and predicted hand trajectories in 3D space for a single trial are also shown in Figure 5 and Figure S3. The best R2 values for X, Y, and Z positioning were 0.7288, 0.7677, and 0.7526, respectively, obtained using 9 electrodes of monkey B.

thumbnail
Figure 4. Decoding results using 32 electrodes selected with the location-based method for monkey B.

Example of 3D hand trajectory prediction over 100 s of test data using 32 channel ECoG signals. R2 values between the predicted (blue) and observed (red) trajectories for X-, Y-, and Z-positions are shown.

https://doi.org/10.1371/journal.pone.0083534.g004

thumbnail
Figure 5. Decoding results for monkey B in three dimensional space.

A) Example of 3D hand trajectory prediction for one trial movement using 32 channel ECoG signals. The predicted and observed trajectories in 3D space are depicted in blue and red, respectively. The unfilled circles represent the start point of movement. The two triangles mark hand position at equivalent time points during movement. Solid circles denote the end point of movement. B, C, and D) The predicted (blue) and observed (red) trajectories shown in the X–Y, X–Z, and Y–Z planes, respectively.

https://doi.org/10.1371/journal.pone.0083534.g005

Figure 6 shows prediction results over 8 s of test data using 3, 6, 9, 12, and 16 ECoG electrodes for monkey B (Fig. S4). With the location-based selection method, 67.3% and 92.9% of the best R2 values were achieved with 3 electrodes for monkeys A and B, respectively. Best R2 percentages using 6 electrodes were 85.3% and 96.9%, and 97.9% and 100% using 9 electrodes.

thumbnail
Figure 6. Decoding results for monkey B with different electrode numbers selected using the location-based method.

Example prediction of 3D hand positioning over 8 s of test data using 3, 6, 9, 12, and 16 electrodes. The red solid line depicts actual trajectories. The green solid line, purple dotted line, light blue dotted line, brown dotted line, and blue solid line represent predicted trajectories using 3, 6, 9, 12, and 16 electrodes, respectively.

https://doi.org/10.1371/journal.pone.0083534.g006

Prediction with the performance-based selection method

For the performance-based method, prediction results for each individual electrode are shown in Figure 7A. For monkey A, R2 values ranged from 0.0903 to 0.2407. The highest R2 value was achieved with electrode 10. For monkey B, R2 values ranged from 0.3566 to 0.6269. The highest R2 value was achieved with electrode 23. Prediction results for each electrode using δ (1.5∼4 Hz), γ3 (90∼120 Hz), and γ4 (120∼150 Hz) bands are shown in Figure 7B, 7C, 7D, respectively. R2 values using the δ band ranged from 0.00 to 0.05 and from −0.06 to 0.37 for monkey A and monkey B, respectively. R2 values using the γ3 band ranged from 0.01 to 0.11 and from 0.01 to 0.50, respectively. The R2 values using γ4 ranged from 0.01 to 0.16 and from 0.01 to 0.47, respectively. Performances for the γ3 and γ4 bands were similar and generally higher than those of the δ band. For both monkeys, the most effective electrodes were concentrated at the lateral areas and areas close to the CS, especially for the γ3 and γ4 bands.

thumbnail
Figure 7. Decoding results for each individual electrode and performance of the two electrode selection methods.

A) Prediction performance for each electrode using all frequency bands. The color map in each electrode represents performance of that electrode. R2 values ranged from 0.00 to 0.24 for monkey A and 0.00 to 0.63 for monkey B. B) Prediction performance of each electrode using the δ band. For monkey A, R2 ranged from 0.00 to 0.24. For monkey B, R2 ranged from 0.00 to 0.63. C) Prediction performance of each electrode using the γ3 band. For monkey A, R2 ranged from 0.00 to 0.24. For monkey B, R2 ranged from 0.00 to 0.63. D) Prediction performance of each electrode using the γ4 band. For monkey A, R2 ranged from 0.00 to 0.24. For monkey B, R2 ranged from 0.00 to 0.63. E) The blue solid and blue dotted lines represent the decoding performance for monkey A using the performance-based method and location-based method, respectively. The red solid and red dotted lines represent the decoding performance for monkey B using the performance-based method and location-based method, respectively. The solid stars and circles denote the highest performances for the performance-based method and location-based method, respectively.

https://doi.org/10.1371/journal.pone.0083534.g007

Summary of the two electrode selection methods

Performance details of two electrode selection methods are shown in Figure 7E. For both monkeys, performance was improved quickly as the number of electrodes used increased from 1 to 9. The performance curves fluctuated only slightly when using 10 electrodes and above. The best R2 values were achieved using 13 and 10 electrodes for monkeys A and B, respectively.

For both methods, the principle is to select more effective electrodes in prediction. As shown in Figure 7A, 7B, 7C, and 7D, higher performance electrodes are concentrated at the lateral areas and near areas of CS. This result is consistent with the principle of the location-based selection method.

To confirm this principle, columns electrodes were also used to predicted hand trajectory. Prediction results in the rostral-caudal direction, and in the medial-lateral direction are shown in Table 1 and Table 2, respectively.

thumbnail
Table 1. Prediction results using location-based electrode selection in the rostral-caudal direction.

https://doi.org/10.1371/journal.pone.0083534.t001

thumbnail
Table 2. Prediction results using location-based electrode selection in the medial-lateral direction.

https://doi.org/10.1371/journal.pone.0083534.t002

The highest performance in the rostral-caudal direction was achieved using column #1 in both monkeys. For monkey A, performance using column #3 was higher than that using #2. This might have been an effect of the presence of electrode 10 (Figure 7A). For monkey B, performance using column #3 was second highest, and performance using column #2 was higher than that using column #4. Highest performance in the medial-lateral direction was achieved using row *2 in both monkeys. For monkey A, performance using row *1 was higher than that using rows *3, *4, and *5. For monkey B, performance using row *3 was second highest, and may have been due to the effect of the δ band at electrode 13 (Figure 7B). Performance using row *1 was higher than that using row *4. Generally, higher performance rows and columns are at the lateral areas and areas near the CS.

Analysis of specific frequency bands

Weights of the nine frequency bands in the prediction model were calculated and are shown in Figure 8A as percent contributions. For monkey B, γ3 (90∼120 Hz) provided the highest contribution. The contributions of δ (1.5∼4 Hz) and γ4 (120∼150 Hz) were higher than those of θ (4∼8 Hz), α (8∼14 Hz), β1 (14∼20 Hz), and β2 (14∼20 Hz).

thumbnail
Figure 8. Contribution and performance of specific frequency bands for monkey B.

A) Each bar represents the weight of each frequency band in the PLS model. From top to bottom, the graphs depict contributions to X, Y, and Z positioning, respectively. B) Prediction results using each frequency band. We performed a two-way ANOVA with effects positioning and frequency bands. Each bar represents the mean R2 value of X, Y, and Z positioning when comparing observed and predicted hand trajectories. Significant differences between mean R2 values are denoted with * (p<0.001).

https://doi.org/10.1371/journal.pone.0083534.g008

3D hand trajectories were predicted by using each frequency band of the ECoG feature signals individually (Table 3). A two-way ANOVA was employed to judge two effects (X, Y, and Z positioning, and the nine frequency bands). No significant differences in prediction performance between X, Y, and Z positioning were observed in both monkeys (monkey A: F2, 16  = 3.61, p  = 0.051; monkey B: F2, 16  = 1.96, p  = 0.173). Significant differences in prediction performance were observed between frequency bands, (monkey A: F8, 16  = 14.16, p  = 6.41×10−6; monkey B: F8, 16  = 52.39, s  = 4.99×10−10), as shown in Figure 8B. The prediction performances using δ, γ2, γ3, and γ4 bands were also significantly higher than that of other bands. Prediction performance of the θ and γ1 bands was significantly higher than that of β2.

thumbnail
Table 3. Prediction results using individual frequency bands.

https://doi.org/10.1371/journal.pone.0083534.t003

Discussion

This study decoded 3D hand trajectories from ECoG signals in Ml and showed that most effective electrodes were concentrated at the lateral areas and areas close to the CS. Comparisons between prediction results suggest that a selection of effective ECoG signals may be better choice than a whole ECoG array. Our results also suggested that ECoG signals are of ample quality and efficiency to control a high performance neural prosthetic.

Which locations are most effective for prediction?

Carmena et. al. (2003) reported that neuron activity recorded from Ml showed greater efficacy than that from dorsal premotor cortex, supplementary motor cortex, posterior parietal cortex, and primary somatosensory cortex. Previous ECoG studies have also used signals mainly from the primary motor area [27][29]. We chose M1 based on those previous results and evaluated the optimal locations in M1. As shown in Figure 7A, ECoG signals from the lateral areas and near areas of CS also showed greater efficacy in prediction, especially in the δ, γ3, and γ4 bands (Figure 7B, 7C, 7D).

How did different numbers of ECoG electrodes affect performance?

As shown in in Figure 7E, the best mean R2 values for monkeys A and B were 0.4805 and 0.7496, respectively, in the location-based selection, and 0.4815 and 0.7780 in the performance-based selection. Both methods, therefore, appear to have equivalent ability to predict hand trajectories.

For both monkeys, performance improved quickly as the number of electrodes used increased from 1 to 9. The performance curves fluctuated only slightly when using 10 electrodes and above. Best decoding performance was achieved using a relatively small number of electrodes, 13 and 10 electrodes in the performance-based selection for monkey A and monkey B, respectively. The performances curves of this study are similar to the results of a previous neuron activity-based study [14], which selected different numbers of high sensitivity neurons in decoding kinematic variables. These results suggest that best decoding performance can be achieved from a relatively small number of effective ECoG signals. However, it should also be noted that decoding performance is not simply related to the electrode number but may more closely depend on the density of electrodes within the effective areas. Still, with the potential utility of wireless transmission technology in ECoG [44], [45], a relatively smaller number of electrodes would provide the benefit of lower power consumption, extending the usage time for wireless BMIs.

Which frequency bands are most effective?

To evaluate the efficacy of specific frequency bands in trajectory decoding, we compared prediction performances of the nine physiologically-based frequency bands with 10 Hz-width fractionized frequency bands from 0 to 150 Hz. The physiologically-based method produced nearly the same or better results (R2  = 0.7328) than the fractionized frequency method (R2  = 0.6815) for monkey B. These results suggest that the usage of physiological frequency bands is more effective than non-physiological fractionized frequency bands.

Weight analysis for the PLS model and the results of decoding performance using each of the nine frequency bands showed that the δ, γ2, γ3, and γ4 bands were more effective than other bands in this study. Previous ECoG studies have shown the importance of the high γ band in motor decoding and BMI control, such as the 60–80 Hz band in prediction of 3D hand trajectories in monkeys [22], [23], 70–110 Hz in controlling a 3D cursor in humans [26], and 56–128 Hz in grasp detection in humans [28]. The importance of the δ band is also supported by our previous ECoG work [34], and is consistent with a previous study [27], which employed a low-frequency band (2–6 Hz) to classify natural grasp types.

Supporting Information

Figure S1.

Predictive error sum of squares in model training for monkey A. The blue line and green line show predictive error sum of squares (PRESS) and R2 values, respectively, for different numbers of latent variables used in the PLS model. The optimal number of 20 is denoted with the red dotted line.

https://doi.org/10.1371/journal.pone.0083534.s001

(EPS)

Figure S2.

Decoding results using 15 electrodes for monkey A. Example of prediction of 3D hand positions during 100 seconds test data by using 15 channel ECoG signals. The R2 value between the predicted (blue) and observed (red) trajectories for X-, Y-, and Z-positions are shown.

https://doi.org/10.1371/journal.pone.0083534.s002

(EPS)

Figure S3.

Decoding results for monkey A in three dimensional space. A) Example of 3D hand trajectory prediction for one trial movement using 3 electrodes. The predicted and observed trajectories in 3D space are depicted in blue and red, respectively. The unfilled circles represent the start point of movement. The two triangles mark hand position at equivalent time points during movement. Solid circles denote the end point of movement. B, C, and D) The predicted (blue) and observed (red) trajectories shown in the X–Y, X–Z, and Y–Z planes, respectively.

https://doi.org/10.1371/journal.pone.0083534.s003

(EPS)

Figure S4.

Decoding results for monkey A with different electrode numbers selected using the location-based method. Example prediction of 3D hand positioning over 8 s of test data using 3, 6, 9, 12, and 15 electrodes. The red solid line depicts actual trajectories. The green solid line, purple dotted line, light blue dotted line, brown dotted line, and blue solid line represent predicted trajectories using 3, 6, 9, 12, and 15 electrodes, respectively.

https://doi.org/10.1371/journal.pone.0083534.s004

(EPS)

Acknowledgments

We thank M. Togawa (NIPS) and Y. Yamanishi (NIPS) for their efforts in this study. We also thank C. S. DaSalla for proofreading the manuscript.

Author Contributions

Conceived and designed the experiments: DS HW AN TI Y. Nishimura YK. Performed the experiments: HW AN TI Y. Nishimura. Analyzed the data: CC DS Y. Nakanishi HK NY. Contributed reagents/materials/analysis tools: CC DS. Wrote the paper: CC DS TI.

References

  1. 1. Wolpaw JR, Mcfarland DJ, Neat GW, Forneris CA (1991) An Eeg-Based Brain-Computer Interface for Cursor Control. Electroencephalography and Clinical Neurophysiology 78: 252–259.
  2. 2. Birbaumer N, Ghanayim N, Hinterberger T, Iversen I, Kotchoubey B, et al. (1999) A spelling device for the paralysed. Nature 398: 297–298.
  3. 3. Chapin JK, Moxon KA, Markowitz RS, Nicolelis MAL (1999) Real-time control of a robot arm using simultaneously recorded neurons in the motor cortex. Nature Neuroscience 2: 664–670.
  4. 4. Wessberg J, Stambaugh CR, Kralik JD, Beck PD, Laubach M, et al. (2000) Real-time prediction of hand trajectory by ensembles of cortical neurons in primates. Nature 408: 361–365.
  5. 5. Taylor DM, Tillery SIH, Schwartz AB (2002) Direct cortical control of 3D neuroprosthetic devices. Science 296: 1829–1832.
  6. 6. Carmena JM, Lebedev MA, Crist RE, O'Doherty JE, Santucci DM, et al. (2003) Learning to control a brain-machine interface for reaching and grasping by primates. PLoS Biol 1: e42.
  7. 7. Wolpaw JR, McFarland DJ (2004) Control of a two-dimensional movement signal by a noninvasive brain-computer interface in humans. Proc Natl Acad Sci USA 101: 17849–17854.
  8. 8. Lebedev MA, Carmena JM, O'Doherty JE, Zacksenhouse M, Henriquez CS, et al. (2005) Cortical ensemble adaptation to represent velocity of an artificial actuator controlled by a brain-machine interface. Journal of Neuroscience 25: 4681–4693.
  9. 9. Hochberg LR, Serruya MD, Friehs GM, Mukand JA, Saleh M, et al. (2006) Neuronal ensemble control of prosthetic devices by a human with tetraplegia. Nature 442: 164–171.
  10. 10. Fagg AH, Hatsopoulos NG, de Lafuente V, Moxon KA, Nemati S, et al. (2007) Biomimetic brain machine interfaces for the control of movement. Journal of Neuroscience 27: 11842–11846.
  11. 11. Velliste M, Perel S, Spalding MC, Whitford AS, Schwartz AB (2008) Cortical control of a prosthetic arm for self-feeding. Nature 453: 1098–1101.
  12. 12. Hochberg LR, Bacher D, Jarosiewicz B, Masse NY, Simeral JD, et al. (2012) Reach and grasp by people with tetraplegia using a neurally controlled robotic arm. Nature 485: 372–375.
  13. 13. Hauschild M, Mulliken GH, Fineman I, Loeb GE, Andersen RA (2012) Cognitive signals for brain-machine interfaces in posterior parietal cortex include continuous 3D trajectory commands. Proc Natl Acad Sci USA 109: 17075–17080.
  14. 14. Sanchez JC, Carmena JM, Lebedev MA, Nicolelis MAL, Harris JG, et al. (2004) Ascertaining the importance of neurons to develop better brain-machine interfaces. IEEE Transactions on Biomedical Engineering 51: 943–953.
  15. 15. Tankus A, Fried I, Shoham S (2012) Sparse decoding of multiple spike trains for brain-machine interfaces. Journal of Neural Engineering 9: 054001.
  16. 16. Rickert J, de Oliveira SC, Vaadia E, Aertsen A, Rotter S, et al. (2005) Encoding of movement direction in different frequency ranges of motor cortical local field potentials. Journal of Neuroscience 25: 8815–8824.
  17. 17. Kellis S, Miller K, Thomson K, Brown R, House P, et al. (2010) Decoding spoken words using local field potentials recorded from the cortical surface. Journal of Neural Engineering 7: 056007.
  18. 18. Thulasidas M, Guan C, Wu JK (2006) Robust classification of EEG signal for brain-computer interface. IEEE Transactions on Neural Systems and Rehabilitation Engineering 14: 24–29.
  19. 19. Hadjidimitriou SK, Hadjileontiadis LJ (2012) Toward an EEG-Based Recognition of Music Liking Using Time-Frequency Analysis. IEEE Transactions on Biomedical Engineering 59: 3498–3510.
  20. 20. Shibata K, Watanabe T, Sasaki Y, Kawato M (2011) Perceptual learning incepted by decoded fMRI neurofeedback without stimulus presentation. Science 334: 1413–1415.
  21. 21. Leuthardt EC, Schalk G, Wolpaw JR, Ojemann JG, Moran DW (2004) A brain-computer interface using electrocorticographic signals in humans. Journal of Neural Engineering 1: 63–71.
  22. 22. Chao ZC, Nagasaka Y, Fujii N (2010) Long-term asynchronous decoding of arm motion using electrocorticographic signals in monkeys. Frontiers in neuroengineering 3: 3.
  23. 23. Shimoda K, Nagasaka Y, Chao ZC, Fujii N (2012) Decoding continuous three-dimensional hand trajectories from epidural electrocorticographic signals in Japanese macaques. Journal of Neural Engineering 9: 036015.
  24. 24. Wilson JA, Felton EA, Garell PC, Schalk G, Williams JC (2006) ECoG factors underlying multimodal control of a brain-computer interface. IEEE Trans Neural Syst Rehabil Eng 14: 246–250.
  25. 25. Chin CM, Popovic MR, Thrasher A, Cameron T, Lozano A, et al. (2007) Identification of arm movements using correlation of electrocorticographic spectral components and kinematic recordings. Journal of Neural Engineering 4: 146–158.
  26. 26. Wang W, Collinger JL, Degenhart AD, Tyler-Kabara EC, Schwartz AB, et al. (2013) An electrocorticographic brain interface in an individual with tetraplegia. PloS ONE 8: e55344.
  27. 27. Pistohl T, Schulze-Bonhage A, Aertsen A, Mehring C, Ball T (2012) Decoding natural grasp types from human ECoG. Neuroimage 59: 248–260.
  28. 28. Pistohl T, Schmidt TSB, Ball T, Schulze-Bonhage A, Aertsen A, et al. (2013) Grasp Detection from Human ECoG during Natural Reach-to-Grasp Movements. PloS ONE 8: e54658.
  29. 29. Yanagisawa T, Hirata M, Saitoh Y, Kishima H, Matsushita K, et al. (2012) Electrocorticographic control of a prosthetic arm in paralyzed patients. Annals of Neurology 71: 353–361.
  30. 30. Chestek CA, Gilja V, Blabe CH, Foster BL, Shenoy KV, et al. (2013) Hand posture classification using electrocorticography signals in the gamma band over human sensorimotor brain areas. Journal of Neural Engineering 10: 026002.
  31. 31. Schalk G, Kubanek J, Miller KJ, Anderson NR, Leuthardt EC, et al. (2007) Decoding two-dimensional movement trajectories using electrocorticographic signals in humans. Journal of Neural Engineering 4: 264–275.
  32. 32. Schalk G, Miller KJ, Anderson NR, Wilson JA, Smyth MD, et al. (2007) Two-dimensional movement control using electrocorticographic signals in humans. Journal of Neural Engineering 5: 75–84.
  33. 33. Pistohl T, Ball T, Schulze-Bonhage A, Aertsen A, Mehring C (2008) Prediction of arm movement trajectories from ECoG-recordings in humans. Journal of Neuroscience Methods 167: 105–114.
  34. 34. Shin D, Watanabe H, Kambara H, Nambu A, Isa T, et al. (2012) Prediction of Muscle Activities from Electrocorticograms in Primary Motor Cortex of Primates. PloS ONE 7: e47992.
  35. 35. Watanabe H, Sato MA, Suzuki T, Nambu A, Nishimura Y, et al. (2012) Reconstruction of movement-related intracortical activity from micro-electrocorticogram array signals in monkey primary motor cortex. Journal of Neural Engineering 9: 036006.
  36. 36. McFarland DJ, McCane LM, David SV, Wolpaw JR (1997) Spatial filter selection for EEG-based communication. Electroencephalography and Clinical Neurophysiology 103: 386–394.
  37. 37. Ludwig KA, Miriani RM, Langhals NB, Joseph MD, Anderson DJ, et al. (2009) Using a common average reference to improve cortical neuron recordings from Microelectrode Arrays. Journal of Neurophysiology 101: 1679–1689.
  38. 38. Wold S, Ruhe A, Wold H, Dunn WJ (1984) The collinearity problem in linear regression – the partial least squares (PLS) approach to generalized Inverses. Siam Journal on Scientific and Statistical Computing 5: 735–743.
  39. 39. Geladi P, Kowalski B (1986) Partial least-squares regression: a tutorial. Analytica Chimica Acta 185: 1–17.
  40. 40. Rosipal R, Kramer N (2006) Overview and recent advances in partial least squares. Subspace, Latent Structure and Feature Selection 3940: 34–51.
  41. 41. Krishnan A, Williams LJ, McIntosh AR, Abdi H (2011) Partial Least Squares (PLS) methods for neuroimaging: a tutorial and review. Neuroimage 56: 455–475.
  42. 42. Rathelot JA, Strick PL (2009) Subdivisions of primary motor cortex based on cortico-motoneuronal cells. Proc Natl Acad Sci USA. 106: 918–23.
  43. 43. Griffin DM, Hudson HM, Belhaj-Saïf A, McKiernan BJ, Cheney PD (2008) Do corticomotoneuronal cells predict target muscle EMG activity? Journal of Neurophysiology 99: 1169–986.
  44. 44. Bjorninen T, Muller R, Ledochowitsch P, Sydanheimo L, Ukkonen L, et al. (2012) Design of wireless links to implanted brain-machine interface microelectronic systems. IEEE Antennas and Wireless Propagation Letters 11: 1663–1666.
  45. 45. Hirata M, Matsushita K, Suzuki T, Yoshida T, Sato F, et al. (2011) A fully-implantable wireless system for human brain-machine interfaces using brain surface electrodes: W-HERBS. IEICE Transactions on Communications E94-B: 2448–2453.