Brainstorm's Suite for Multi-unit Electrophysiology

Author: Konstantinos Nasiotis

ephys_intro.png

Introduction

We feature a series of data analysis tools specifically dedicated to Invasive electrophysiology research (e-phys). We see this suite of expanded Brainstorm features as an opportunity for enabling a methodological continuum between (multi)cellular and systems level (with EEG and MEG) data analytics and facilitate the scientific dialogue and interactions between subdisciplines in the field.

Invasive electrophysiological recordings capture the spiking activity of individual or small groups of neurons. Since spiking events are transient phenomena (the typical duration of spikes is 1ms), the sampling rate of e-phys recording systems is typically much faster than for scalp EEG and MEG recordings. State-of-the-art acquisition systems can accommodate simultaneous recordings from hundreds of channels over several multi-electrode arrays positioned across the preparation, ultimately producing vast amounts of data.

We foresee that Brainstorm will contribute to the efficient data management and processing when there is a need for efficient, open and reproducible data management and utilization of large quantities of data. Users that have Matlab’s parallel processing toolbox, will greatly benefit from the utilization of its properties. Most functions of the electrophysiology toolbox have a checkbox where users can enable parallelization of the processes.

Additionally, it is highly recommended for users that need to analyze large files, to select a large number for memory utilization on any processes that allow users to adjust the memory used, since the speed of the processes (especially the spike sorting procedures that use demultiplexing of the raw signal) can be greatly improved by the extra memory used.

Note that the operations used here are not detailed, the goal of this tutorial is not to introduce Brainstorm to new users. For in-depth explanations of the interface and theoretical foundations, please refer to the introduction tutorials.

How to cite

Researchers publishing results obtained with this electrophysiology toolbox are kindly asked to cite:

Tutorial dataset

Description

This tutorial uses a dataset collected by Nardin Nakhla in Dr. Christopher Pack’s lab, from a macaque monkey with a penetration chamber implanted in areas MT and MST. The recording probe consisted of 32 channels and 1 reference electrode.

The task for the animal was to maintain fixation on a target while a superflow stimulus (dots that rotate, expand or translate in a partially coherent mode) was presented in one of 9 possible visual locations (as shown in the figure below). The events from the stimulation system are saved along with pre-spike-sorted events.

Download the dataset

Files



ephys_import.gif

Create protocol

Import anatomy (Optional)

Users can now import the anatomy of the monkey if needed. However, the pipeline described in this tutorial does not require an anatomical MRI, you may skip this section to save some time.

floyd_mri.gif floyd_cortex.gif

Import recordings

Importing of this file might last a few minutes. Plexon .plx files cannot be accessed efficiently. In order to avoid waiting times while scrolling through the file, the .plx files are first converted into .bst files (Brainstorm binary format). This conversion causes the initial delay, but improves the user experience when reviewing the file. This is not typically the case for files from other acquisition systems.

Double-click on the Link to raw file to open the recordings viewer. This is a very easy way to inspect the recordings for artifacts or stimulation effects since only the selected part is loaded in memory on-the-fly, without importing anything in the database. Any events captured from the acquisition system are automatically loaded from the importer.

displayRaw_new.gif

Existing spike sorter results

It is important to note that some acquisition systems might have enabled their own automatic spike sorter. These events are automatically labeled with the convention that is used throughout the e-phys toolbox:

“Spikes Channel ChannelLabel”, where ChannelLabel is the Label of the electrode assigned by the acquisition system. If the acquisition system assigned the spikes’ waveforms into more than one cluster/neuron during this automatic clustering, the neuronal events would be (e.g. for an electrode that picked up spikes from 3 neurons): “Spikes Channel ChannelLabel |1|”, “Spikes Channel ChannelLabel |2|”, “Spikes Channel ChannelLabel |3|”.

For more information, users can visit this page.

The events are comprised of the signals that the stimulation system has sent (Stim On, Start, Stop etc.) and the pre-spike-sorted events that the Plexon acquisition system can perform (Spikes Channel AD01, Spike Channel AD02 etc.).

Slow interface

The Brainstorm interface tends to become slower when the number of events increase. Opening a file for visualization or opening the Pipeline editor window may freeze Matlab for several seconds when handling files with hundreds of thousands of markers. If this becomes too heavy, don't hesitate to delete the event groups you have no use for.

Import events

In this tutorial dataset, the stimulus presentation events were coded from a software that was not triggering events on the acquisition system, so they need to be manually imported.

At this stage, users can start performing their analysis based on the experimental events.

Import electrode positions (Optional)

If you have access to the 3D positions of the electrodes with which the data were recorded, you add this information to the database. This would enable additional display and processing options. In this dataset, the electrodes positions are available in "world" coordinates.

floyd_elec.gif








Feedback: Comments, bug reports, suggestions, questions
Email address (if you expect an answer):


e-phys/Introduction (last edited 2022-05-03 11:54:46 by FrancoisTadel)