Nonmotor regions encode path-related information during movements

Stereo-EEG recordings extend known distributions of canonical movement-related oscillations

Objective.Previous electrophysiological research has characterized canonical oscillatory patterns associated with movement mostly from recordings of primary sensorimotor cortex. Less work has attempted to decode movement based on electrophysiological recordings from a broader array of brain areas such as those sampled by stereoelectroencephalography (sEEG), especially in humans. We aimed to identify and characterize different movement-related oscillations across a relatively broad sampling of brain areas in humans and if they extended beyond brain areas previously associated with movement.Approach.We used a linear support vector machine to decode time-frequency spectrograms time-locked to movement, and we validated our results with cluster permutation testing and common spatial pattern decoding.Main results.We were able to accurately classify sEEG spectrograms during a keypress movement task versus the inter-trial interval. Specifically, we found these previously-described patterns: beta (13-30 Hz) desynchronization, beta synchronization (rebound), pre-movement alpha (8-15 Hz) modulation, a post-movement broadband gamma (60-90 Hz) increase and an event-related potential. These oscillatory patterns were newly observed in a wide range of brain areas accessible with sEEG that are not accessible with other electrophysiology recording methods. For example, the presence of beta desynchronization in the frontal lobe was more widespread than previously described, extending outside primary and secondary motor cortices.Significance.Our classification revealed prominent time-frequency patterns which were also observed in previous studies that used non-invasive electroencephalography and electrocorticography, but here we identified these patterns in brain regions that had not yet been associated with movement. This provides new evidence for the anatomical extent of the system of putative motor networks that exhibit each of these oscillatory patterns.

Shared and distinct voxel-based lesion-symptom mappings for spasticity and impaired movement in the hemiparetic upper limb

Hemiparesis and spasticity are common co-occurring manifestations of hemispheric stroke. The relationship between impaired precision and force in voluntary movement (hemiparesis) and the increment in muscle tone that stems from dysregulated activity of the stretch reflex (spasticity) is far from clear. Here we aimed to elucidate whether variation in lesion topography affects hemiparesis and spasticity in a similar or dis-similar manner. Voxel-based lesion-symptom mapping (VLSM) was used to assess the impact of lesion topography on (a) upper limb paresis, as reflected by the Fugl-Meyer Assessment scale for the upper limb and (b) elbow flexor spasticity, as reflected by the Tonic Stretch Reflex Threshold, in 41 patients with first-ever stroke. Hemiparesis and spasticity were affected by damage to peri-Sylvian cortical and subcortical regions and the putamen. Hemiparesis (but not spasticity) was affected by damage to the corticospinal tract at corona-radiata and capsular levels, and by damage to white-matter association tracts and additional regions in the temporal cortex and pallidum. VLSM conjunction analysis showed only a minor overlap of brain voxels where the existence of damage affected both hemiparesis and spasticity, suggesting that control of voluntary movement and regulation of muscle tone at rest involve largely separate parts of the motor network.

Non-motor Brain Regions in Non-dominant Hemisphere Are Influential in Decoding Movement Speed

Sensorimotor control studies have predominantly focused on how motor regions of the brain relay basic movement-related information such as position and velocity. However, motor control is often complex, involving the integration of sensory information, planning, visuomotor tracking, spatial mapping, retrieval and storage of memories, and may even be emotionally driven. This suggests that many more regions in the brain are involved beyond premotor and motor cortices. In this study, we exploited an experimental setup wherein activity from over 87 non-motor structures of the brain were recorded in eight human subjects executing a center-out motor task. The subjects were implanted with depth electrodes for clinical purposes. Using training data, we constructed subject-specific models that related spectral power of neural activity in six different frequency bands as well as a combined model containing the aggregation of multiple frequency bands to movement speed. We then tested the models by evaluating their ability to decode movement speed from neural activity in the test data set. The best models achieved a correlation of 0.38 ± 0.03 (mean ± standard deviation). Further, the decoded speeds matched the categorical representation of the test trials as correct or incorrect with an accuracy of 70 ± 2.75% across subjects. These models included features from regions such as the right hippocampus, left and right middle temporal gyrus, intraparietal sulcus, and left fusiform gyrus across multiple frequency bands. Perhaps more interestingly, we observed that the non-dominant hemisphere (ipsilateral to dominant hand) was most influential in decoding movement speed.

An exploratory data analysis method for identifying brain regions and frequencies of interest from large-scale neural recordings

High-resolution whole brain recordings have the potential to uncover unknown functionality but also present the challenge of how to find such associations between brain and behavior when presented with a large number of regions and spectral frequencies. In this paper, we propose an exploratory data analysis method that sorts through a massive quantity of multivariate neural recordings to quickly extract a subset of brain regions and frequencies that encode behavior. This approach combines existing tools and exploits low-rank approximation of matrices without a priori selection of regions and frequency bands for analysis. In detail, the spectral content of neural activity across all frequencies of each recording contact is computed and represented as a matrix. Then, the rank-1 approximation of the matrix is computed using singular value decomposition and the associated singular vectors are extracted. The temporal singular vector, which captures the salient features of the spectrogram, is then correlated to the trial-varying behavioral signal. The distribution of correlations for each brain region is efficiently computed and used to find a subset of regions and frequency bands of interest for further examination. As an illustration, we apply this approach to a data set of local field potentials collected using stereoelectroencephalography from a human subject performing a reaching task. Using the proposed procedure, we produced a comprehensive set of brain regions and frequencies related to our specific behavior. We demonstrate how this tool can produce preliminary results that capture neural patterns related to behavior and aid in formulating data-driven hypotheses, hence reducing the time it takes for any scientist to transition from the exploratory to the confirmatory phase.