View: |
Part 1: Document Description
|
Citation |
|
---|---|
Title: |
Dataset for the study "Slow insertion of silicon probes improves the quality of acute neuronal recordings" |
Identification Number: |
hdl:21.15109/ARP/GBWSNE |
Distributor: |
ARP |
Date of Distribution: |
2021-11-11 |
Version: |
4 |
Bibliographic Citation: |
Fiáth, Richárd; Márton, Adrienn Lilla; Mátyás, Ferenc; Pinke, Domonkos; Márton, Gergely; Tóth, Kinga; Ulbert, István, 2021, "Dataset for the study "Slow insertion of silicon probes improves the quality of acute neuronal recordings"", https://hdl.handle.net/21.15109/ARP/GBWSNE, ARP, V4, UNF:6:/P2P9rrxrjZGxPQsAPYkNQ== [fileUNF] |
Citation |
|
Title: |
Dataset for the study "Slow insertion of silicon probes improves the quality of acute neuronal recordings" |
Identification Number: |
hdl:21.15109/ARP/GBWSNE |
Authoring Entity: |
Fiáth, Richárd (ELKH Természettudományi Kutatóközpont) |
Márton, Adrienn Lilla (ELKH Természettudományi Kutatóközpont) |
|
Mátyás, Ferenc (ELKH Természettudományi Kutatóközpont) |
|
Pinke, Domonkos (ELKH Természettudományi Kutatóközpont) |
|
Márton, Gergely (ELKH Természettudományi Kutatóközpont) |
|
Tóth, Kinga (ELKH Természettudományi Kutatóközpont) |
|
Ulbert, István (ELKH Természettudományi Kutatóközpont) |
|
Other identifications and acknowledgements: |
Horváth, Csaba |
Grant Number: |
2017-1.2.1-NKP-2017-00002; KTIA_NAP_13-2-2015-0010 |
Grant Number: |
VEKOP-2.3.2-16-2017-00013 |
Grant Number: |
PD124175, PD121015, PD121123, FK124434 and KKP126998 |
Grant Number: |
Pilot project to archive research data |
Distributor: |
ARP |
Access Authority: |
Fiáth, Richárd |
Depositor: |
Fiáth, Richárd |
Date of Deposit: |
2021-07-12 |
Holdings Information: |
https://hdl.handle.net/21.15109/ARP/GBWSNE |
Study Scope |
|
Keywords: |
Medicine, Health and Life Sciences, Silicon probe, High-density recordings, Neocortex, Slow oscillation, Electrophysiological recording, Neural implant, Signal-to-noise ratio, fakeDoiMigrationDone |
Topic Classification: |
Neural implants |
Abstract: |
Here we provide the dataset for our study "Slow insertion of silicon probes improves the quality of acute neuronal recordings". Abstract of the study: Neural probes designed for extracellular recording of brain electrical activity are traditionally implanted with an insertion speed between 1 µm/s and 1 mm/s into the brain tissue. Although the physical effects of insertion speed on the tissue are well studied, there is a lack of research investigating how the quality of the acquired electrophysiological signal depends on the speed of probe insertion. In this study, we used four different insertion speeds (0.002 mm/s, 0.02 mm/s, 0.1 mm/s, 1 mm/s) to implant high-density silicon probes into deep layers of the somatosensory cortex of ketamine/xylazine anesthetized rats. After implantation, various qualitative and quantitative properties of the recorded cortical activity were compared across different speeds in an acute manner. Our results demonstrate that after the slowest insertion both the signal-to-noise ratio and the number of separable single units were significantly higher compared with those measured after inserting probes at faster speeds. Furthermore, the amplitude of recorded spikes as well as the quality of single unit clusters showed similar speed-dependent differences. Post hoc quantification of the neuronal density around the probe track showed a significantly higher number of NeuN-labelled cells after the slowest insertion compared with the fastest insertion. Our findings suggest that advancing rigid probes slowly (~1 µm/s) into the brain tissue might result in less tissue damage, and thus in neuronal recordings of improved quality compared with measurements obtained after inserting probes with higher speeds. |
Notes: |
<h3>Repository structure</h3> <p>Each recording and corresponding metadata, single unit properties and quality metrics were packaged in the <a href="https://www.nwb.org/">Neurodata Without Borders: Neurophysiology version 2.0 (NWB:N 2.0) data format</a> using the <a href="https://github.com/NeurodataWithoutBorders/matnwb/">MatNWB API</a>. A single NWB file was created from each recording. NWB files were placed in folders based on the used silicon probe (32-channel or 128-channel probe), insertion speed (e.g, from IS1 to IS4) and the identifier of the animal (e.g., from Rat01 to Rat10). The filename of the NWB file (identifier) was constructed by concatenating the probe type, insertion speed and the animal ID (e.g., 128ch_IS1_Rat01). The insertion speeds are the followng:<br><br> <b>32channel_silicon_probe_recordings:</b><br> Insertion speed 1 (IS1): 2 μm/s (0.002 mm/s)<br> Insertion speed 2 (IS2): 1 mm/s<br><br> <b>128channel_silicon_probe_recordings:</b><br> Insertion speed 1 (IS1): 2 μm/s (0.002 mm/s)<br> Insertion speed 2 (IS2): 20 μm/s (0.02 mm/s)<br> Insertion speed 3 (IS3): 100 μm/s (0.1 mm/s)<br> Insertion speed 4 (IS4): 1 mm/s<br><br> <b>128channel_silicon_probe_recordings_with_histology:</b><br> Insertion speed 1 (IS1): 2 μm/s (0.002 mm/s)<br> Insertion speed 2 (IS2): 1 mm/s<br><br> </p> In the case of the 128-channel recordings located in the '128channel_silicon_probe_recordings' folder, spike sorting was carried out again using a newer version of the spike sorting software (Kilosort2). All other recordings contain the original spike sorting results done with Kilosort1 and published in the paper.<br><br> <p>The CSV file named "Animal_characteristics_and_targeted_cortical_areas" contains information about the weight and sex of animals, as well as about the cortical area and stereotaxic coordinates corresponding to each probe insertion. The "Recording_characteristics" CSV file lists several useful properties for each NWB file including the file size, the duration of the recording, the cortical location, the single unit yield and the average signal-to-noise ratio of single units.</p> <br> <h3>NWB file structure</h3> <p>Each NWB file contains several main groups which are similar to directories. The <i>acquisition</i> group contains the continuous wideband 32- and 128-channel data (‘wideband_multichannel_recording’) in a compressed form, as well as several parameters related to the raw data such as the measurement unit or the data conversion number. The <i>general</i> group contains metadata about the experiments and consists of several subgroups, related to the recording probe (‘general/devices’; ‘general/extracellular_ephys’) or the subjects of the experiments (‘general/subject’). Former subgroups carry information about the probe location (brain area and stereotaxic coordinates) and the relative positions of recordings sites, while the latter contains metadata about the animal (e.g., sex, species, subject ID, or weight). Information about spike sorting and single units and corresponding data are available in the <i>units</i> group. For each unit, we included here the mean and standard deviation of their spike waveform on all channels, calculated both from the filtered (‘mean_waveform_all_channels_filt’; ’waveform_sd_all_channels_filt’) and the wideband data (‘mean_waveform_all_channels_raw’; ’ waveform_sd_all_channels_raw’). Furthermore, the spike waveform recorded on the channel with the largest spike (i.e., peak waveform channel) was saved separately (‘mean_waveform_peak_channel_filt’, ‘mean_waveform_peak_channel_raw’). Several single unit properties and cluster quality metrics, as well as the spike times and spike count of each unit were saved in the <i>units</i> group. Furthermore, to aid users in selecting and analyzing a subset of this dataset appropriate for their research goals, we also created an NWB file (‘allSingleUnits.nwb’) which contains all single units collected for a particular insertion speed, with all the properties listed above, along with the identifier of the recording (‘units/session_id’) and the cortical area (‘units/cortical_area’) they originate from. The structure of NWB files can be explored using the <a href="https://www.hdfgroup.org/downloads/hdfview/">HDFView</a> software.<p> <br> <h3>Interacting with the NWB files</h3> <p>Users can import data from NWB files using the <a href="https://github.com/NeurodataWithoutBorders/pynwb/">PyNWB</a> and <a href="https://github.com/NeurodataWithoutBorders/matnwb/">MatNWB</a> APIs, or using <a href="https://github.com/SpikeInterface">SpikeInterface</a>. Loaded samples of the raw data have to be multiplied by a conversion number (0.195) to get the amplitudes in microvolts. Here we provide some examples how users can import data from NWB files using the MATLAB-based MatNWB API.</p> <br> Loading a short segment (20.000 samples corresponding to 1 second of data) of the raw wideband recording on all (128) channels: <br> <br> <code> 1. nwb = nwbRead('128ch_IS1_Rat01.nwb'); </code> <br> <code> 2. dataChunk=nwb.acquisition.get('wideband_multichannel_recording').data.load([1, 1], [128, 20000]); </code> <br> <br> To load data from a 32-channel recording just use 32 instead of 128. The recording site-channel mapping can be found in the published paper (see above at the field "Related publication"). <br><br> It is important to note that TimeSeries data types in NWB files are stored with time in the first dimension and channels in the second, but dimensions are reversed in MatNWB. <br> Loading and plotting the mean spike waveform of a specific single unit on the peak waveform channel: <br> <br> <code> 3. peakChannels = nwb.units.vectordata.get('peak_waveform_channel').data.load(); </code> <br> <code> 4. meanWaveforms = nwb.units.vectordata.get('mean_waveform_all_channels_filt').data.load(); </code> <br> <code> 5. mySingleUnit = 11; </code> <br> <code> 6. singleUnitWaveform = meanWaveforms(peakChannels(mySingleUnit), :, mySingleUnit); </code> <br> <code> 7. plot(singleUnitWaveform); </code> <br> <br> Loading the isolation distance quality metric of all units found in a single NWB file: <br> <br> <code> 8. IDvalues = nwb.units.vectordata.get('isolation_distance').data.load(); </code> <br> <br> Spike times are stored in a special structure called ragged arrays consisting of two vectors. The spike_times vector contains all spike times of all single units concatenated one after the other, while the spike_times_index vector stores where the spike times of individual single units are located in the spike_index vector (see also <a href="https://neurodatawithoutborders.github.io/matnwb/tutorials/html/ecephys.html#H_97F533F8">https://neurodatawithoutborders.github.io/matnwb/tutorials/html/ecephys.html#H_97F533F8</a>). <br> We can load the spikes times of a specific single unit (in seconds) the following way: <br> <br> <code> 9. allSpikeTimes = nwb.units.spike_times.data.load(); </code> <br> <code> 10. spikeTimesIndex = nwb.units.spike_times_index.data.load(); </code> <br> <code> 11. spikesOfSingleUnit2 = allSpikeTimes(spikeTimesIndex(1)+1 : spikeTimesIndex(2)); </code> <br> <br> SpikeInterface can also be used to load the wideband data and single unit properties (in Python; works only with version 0.13): <br> <br> <code> 1. import spikeextractors as se </code> <br> <code> 2. nwbPath = '128ch_IS1_Rat01.nwb' </code> <br> <code> 3. recording = se.NwbRecordingExtractor(nwbPath) </code> <br> <code> 4. sorting = se.NwbSortingExtractor(nwbPath) </code> <br> <code> 5. mySingleUnit = 2 </code> <br> <code> 6. sorting.get_unit_property(mySingleUnit,'isolation_distance') </code> <br> <br> We also provide a Matlab script ("NWB_tutorial_script.m") which can be used to load, visualize and preprocess (e.g., filter) files in the dataset using MatNWB. <br> |
Methodology and Processing |
|
Sources Statement |
|
Data Access |
|
Other Study Description Materials |
|
Related Publications |
|
Citation |
|
Title: |
Fiáth R, Márton AL, Mátyás F, Pinke D, Márton G, Tóth K, Ulbert I. Slow insertion of silicon probes improves the quality of acute neuronal recordings. (2019) Scientific Reports 9: 111 |
Identification Number: |
10.1038/s41598-018-36816-z |
Bibliographic Citation: |
Fiáth R, Márton AL, Mátyás F, Pinke D, Márton G, Tóth K, Ulbert I. Slow insertion of silicon probes improves the quality of acute neuronal recordings. (2019) Scientific Reports 9: 111 |
File Description--f6070 |
|
File: Animal_characteristics_and_targeted_cortical_areas.tab |
|
|
|
Notes: |
UNF:6:nmcbHuRlz9YXzyZYoA5NdQ== |
File Description--f6071 |
|
File: Recording_characteristics.tab |
|
|
|
Notes: |
UNF:6:sr0ZjD7dYJzglxm491Rp1Q== |
List of Variables: |
|
Variables |
|
f6070 Location: |
Variable Format: character Notes: UNF:6:jxTn5a6/XARCQccxfczaEw== |
f6070 Location: |
Summary Statistics: Min. 190.0; Mean 287.5438596491228; StDev 50.22380986837056; Valid 57.0; Max. 390.0 Variable Format: numeric Notes: UNF:6:2SLF4bmSAirvlJkyIN+Ngw== |
f6070 Location: |
Variable Format: character Notes: UNF:6:9FnfaZp7jWqe/m6uF5363g== |
f6070 Location: |
Variable Format: character Notes: UNF:6:IJLUqR3snAZD2O/yKc7E3w== |
f6070 Location: |
Summary Statistics: StDev 0.7273582766811976; Min. -5.2; Mean -2.807719298245614; Valid 57.0; Max. -1.0; Variable Format: numeric Notes: UNF:6:Z4H25d3zbDjuJxmaQqwGlg== |
f6070 Location: |
Summary Statistics: Max. 3.55; Valid 57.0; Mean 2.6228070175438596; Min. 1.97; StDev 0.40182954091240575 Variable Format: numeric Notes: UNF:6:G/dpGCXKOXnqe9f7fwcOZA== |
f6070 Location: |
Variable Format: character Notes: UNF:6:B8pP0BNCVIpY5GG6kB2cIA== |
f6070 Location: |
Summary Statistics: StDev 0.16329931618554513; Valid 10.0; Min. -3.0; Mean -2.8; Max. -2.5; Variable Format: numeric Notes: UNF:6:t3cOojXKCUymztcE4imgWQ== |
f6070 Location: |
Summary Statistics: Max. 3.1; Valid 10.0; Mean 2.75; Min. 2.29; StDev 0.3416300793418389; Variable Format: numeric Notes: UNF:6:FFFio3bMgBTJpfZr25NB8w== |
f6071 Location: |
Variable Format: character Notes: UNF:6:WKl32paipUewJITisCjZqQ== |
f6071 Location: |
Summary Statistics: StDev 3.360416970584419; Mean 7.57089552238806; Max. 10.4; Min. 2.39; Valid 67.0; Variable Format: numeric Notes: UNF:6:R0hD8KkT8fFU+IetYJt8+w== |
f6071 Location: |
Summary Statistics: Max. 50.8; StDev 0.8074732073275169; Mean 45.3444776119403; Valid 67.0; Min. 45.0 Variable Format: numeric Notes: UNF:6:9kFBEBKdWQPpZU0qESXPlQ== |
f6071 Location: |
Variable Format: character Notes: UNF:6:vGqtiKPBHuZigdMxvveuLw== |
f6071 Location: |
Summary Statistics: Max. 183.0; StDev 37.12359715060317; Mean 48.537313432835816; Min. 5.0; Valid 67.0 Variable Format: numeric Notes: UNF:6:v7W109fOcRoHPzPcYCxQrw== |
f6071 Location: |
Summary Statistics: Valid 67.0; Max. 26.64; Mean 11.974029850746268; Min. 4.09; StDev 3.913053159541595 Variable Format: numeric Notes: UNF:6:Gj95SS3UtMsuKt95czd07w== |
Label: |
NWB_tutorial_scipt.m |
Notes: |
text/plain |
Label: |
128ch_IS1_Rat01.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS1_Rat02.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS1_Rat03.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS1_Rat04.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS1_Rat05.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS1_Rat06.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS1_Rat07.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS1_Rat08.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS1_Rat09.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS1_Rat10.nwb |
Notes: |
application/octet-stream |
Label: |
allSingleUnits.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS2_Rat01.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS2_Rat02.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS2_Rat03.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS2_Rat04.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS2_Rat05.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS2_Rat06.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS2_Rat07.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS2_Rat08.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS2_Rat09.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS2_Rat10.nwb |
Notes: |
application/octet-stream |
Label: |
allSingleUnits.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS3_Rat01.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS3_Rat02.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS3_Rat03.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS3_Rat04.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS3_Rat05.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS3_Rat06.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS3_Rat07.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS3_Rat08.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS3_Rat09.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS3_Rat10.nwb |
Notes: |
application/octet-stream |
Label: |
allSingleUnits.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS4_Rat01.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS4_Rat02.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS4_Rat03.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS4_Rat04.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS4_Rat05.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS4_Rat06.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS4_Rat08.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS4_Rat09.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS4_Rat10.nwb |
Notes: |
application/octet-stream |
Label: |
allSingleUnits.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS1_Rat01_H.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS1_Rat02_H.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS1_Rat03_H.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS1_Rat04_H.nwb |
Notes: |
application/octet-stream |
Label: |
allSingleUnits.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS2_Rat01_H.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS2_Rat02_H.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS2_Rat03_H.nwb |
Notes: |
application/octet-stream |
Label: |
128ch_IS2_Rat04_H.nwb |
Notes: |
application/octet-stream |
Label: |
allSingleUnits.nwb |
Notes: |
application/octet-stream |
Label: |
32ch_IS1_Rat01_Insertion1.nwb |
Notes: |
application/octet-stream |
Label: |
32ch_IS1_Rat01_Insertion2.nwb |
Notes: |
application/octet-stream |
Label: |
32ch_IS1_Rat02_Insertion1.nwb |
Notes: |
application/octet-stream |
Label: |
32ch_IS1_Rat02_Insertion2.nwb |
Notes: |
application/octet-stream |
Label: |
32ch_IS1_Rat03_Insertion1.nwb |
Notes: |
application/octet-stream |
Label: |
32ch_IS1_Rat03_Insertion2.nwb |
Notes: |
application/octet-stream |
Label: |
32ch_IS1_Rat04_Insertion1.nwb |
Notes: |
application/octet-stream |
Label: |
32ch_IS1_Rat04_Insertion2.nwb |
Notes: |
application/octet-stream |
Label: |
32ch_IS1_Rat05_Insertion1.nwb |
Notes: |
application/octet-stream |
Label: |
32ch_IS1_Rat05_Insertion2.nwb |
Notes: |
application/octet-stream |
Label: |
allSingleUnits.nwb |
Notes: |
application/octet-stream |
Label: |
32ch_IS2_Rat01_Insertion1.nwb |
Notes: |
application/octet-stream |
Label: |
32ch_IS2_Rat01_Insertion2.nwb |
Notes: |
application/octet-stream |
Label: |
32ch_IS2_Rat02_Insertion1.nwb |
Notes: |
application/octet-stream |
Label: |
32ch_IS2_Rat02_Insertion2.nwb |
Notes: |
application/octet-stream |
Label: |
32ch_IS2_Rat03_Insertion1.nwb |
Notes: |
application/octet-stream |
Label: |
32ch_IS2_Rat03_Insertion2.nwb |
Notes: |
application/octet-stream |
Label: |
32ch_IS2_Rat04_Insertion1.nwb |
Notes: |
application/octet-stream |
Label: |
32ch_IS2_Rat04_Insertion2.nwb |
Notes: |
application/octet-stream |
Label: |
32ch_IS2_Rat05_Insertion1.nwb |
Notes: |
application/octet-stream |
Label: |
32ch_IS2_Rat05_Insertion2.nwb |
Notes: |
application/octet-stream |
Label: |
allSingleUnits.nwb |
Notes: |
application/octet-stream |