Myomatrix arrays for high-definition muscle recording
Chung et al., eLife (2023)

All data and code described below can be found here: https://doi.org/10.5061/dryad.66t1g1k70

---

Overview of data structures

This data archive includes two electromyography (EMG) datasets recorded with Myomatrix arrays from behaving animals. These are the mouse triceps data shown in Figure 1(b,d) and Supplemental Figure 2(a,b,e,f), and the rhesus macaque biceps data shown in Figure 3 and Supplemental Figure 2g. Please see the manuscript for descriptions of how each dataset was recorded.

Data files are provided in both .h5 (Hierarchical Data Format 5) and .m (MATLAB) format:

Approximately 2 minutes of data recorded from a mouse running on a treadmill (file size roughly 0.4 GB)

Myomatrix_eLife_2023_EMG_DATA_mouse_triceps [.mat,.h5]

Approximately 100 minutes of data from the macaque performing planar reaching movements (file size roughly 3.8 GB for .mat, 11.2 GB for .h5)

Myomatrix_eLife_2023_EMG_DATA_macaque_biceps [.mat,.h5]

Spatial channel maps of the injectable and chronically-implanted devices used in the macaque and mouse dataset, respectively:

Myomatrix_eLife_2023_CHANNEL_MAP_mouse_triceps [.mat,.h5]

Myomatrix_eLife_2023_CHANNEL_MAP_macaque_biceps [.mat,.h5]

---

Guide to variables (EMG data files):

Each dataset (Myomatrix_eLife_2023_EMG_DATA_macaque_biceps and Myomatrix_eLife_2023_EMG_DATA_mouse_triceps) includes the following:

Variable name: EMG_data
Description: Myomatrix data (in units of microvolts) in an NxM matrix, where N is the number of time samples, and M is the number of recorded channels. M=32 in the macaque data (unipolar recordings from 32 electrode contacts), and M=16 in the mouse data (bipolar signals from 32 electrode contacts; see Note 2 below).

Notes:

1. The macaque dataset was high-pass filtered with a cutoff of 250Hz. The mouse dataset was bandpass filtered between 300-7500 Hz (see Methods).
2. The mouse data was initially recorded in unipolar mode; we subsequently subtracted the unipolar signals from adjacent channels to approximate bipolar recordings.
3. To reduce file size, EMG_data in the macaque dataset is saved as a 16-bit integer (all other variables are saved as double-precision values).
4. In the macaque recording, three of the 32 electrode contacts on the injectable Myomatrix device had impedances >1 MOhm, indicating broken or otherwise nonfunctional electrode contacts. The voltage values from these (channels 1, 4, and 17, provided in columns [1,4,17] in EMG_data) have therefore been set to zero.

Variable name: sampling_rate
Description: Sampling frequency of EMG data (in Hz).

Variable name: sample_ID_vec
Description: Time (in units of sample number, divide by sampling_rate to covert to units of time) of each spike waveform isolated by the modified version of Kilosort described in the manuscript. Note that by design, Kilosort does not necessarily label every waveform, only those that can be confidently assigned to waveform clusters.

Variable name: unit_ID_vec
Description: Motor unit identity of each sorted spike waveform. Units are numbered in decreasing order of the total number of spikes fired in the dataset, i.e. the unit that fired the most spikes is defined as unit 1, the unit that fired the second-most spikes is defined as unit 2, etc. As described the text, the mouse dataset includes a total of 6 identified motor units and the macaque dataset includes 13 sorted units.

---

Guide to variables (channel maps):
Each channel map file (Myomatrix_eLife_2023_CHANNEL_MAP_mouse_triceps and Myomatrix_eLife_2023_CHANNEL_MAP_macaque_biceps) includes the following, which can be used by Kilosort or other spike-sorting approaches to define the geometry of the electrode array:

Variable name: xcoords, ycoords

Description: These variables describe the spatial location of the recording channels, can be used by Kilosort to constrain the spatial extent over which the same unit may appear (see Methods). For the macaque dataset, xcoords and ycoords define the coordinates, in units of microns, of each of the 32 electrode contacts (which were recorded in unipolar mode, see manuscript). For the mouse dataset, $xcoords$ and $ycoords$ define the locations at the midpoint between each of the 16 pairs of electrode contacts that are subtracted to produce the bipolar signal.

Variable name: kcoords
Description: In Kilosort, kcoords provides an index defining on which electrode shank each channel is located.

Variable name: chanMap, chanMap0ind
Description: Index of channel numbers; either 1-indexed (e.g. channel 1-32) or 0-indexed (e.g. channel 0-31).

Variable name: connected

Description: In Kilosort, working channels indexed with connected=1, broken channels indexed with connected=0. As described above, in the macaque dataset, channels 1, 4, and 17 were broken, and the corresponding columns of EMG_data are set to zero. Note however that in our modified version of Kilosort, we set all connected values to 1 to maintain the size and shape of the channel map.

---

Included MATLAB code:

We provide MAKEFIG_myomatrix_EMG_and_sorting.m, which plots subsets of time (and channels) to allow users to view the EMG data along with our sorting After loading either

Myomatrix_eLife_2023_EMG_DATA_mouse_triceps.mat

or

Myomatrix_eLife_2023_EMG_DATA_macaque_biceps.mat

into the MATLAB workspace, run the following:

MAKEFIG_myomatrix_EMG_and_sorting(EMG_data,unit_ID_vec,sample_ID_vec,sampling_rate)

This code will create two figures. The first figure will show a subset of EMG channels during 20 seconds of recording, as well as the outcome of our spike sorting (see manuscript for spike sorting details. The second figure will show a summary of the average firing rate of each motor unit across the duration of the recording session. Modifying this code will allow users to view other time periods or subsets of channels.