Dataset of cortical activity recorded with high spatial resolution from anesthetized ratshdl:21.15109/ARP/0QRUQKARP2021-07-193Csaba, Horváth; Lili Fanni, Tóth; István, Ulbert; Fiáth, Richárd, 2021, "Dataset of cortical activity recorded with high spatial resolution from anesthetized rats", https://hdl.handle.net/21.15109/ARP/0QRUQK, ARP, V3, UNF:6:9YijdOOa6/UIjD4UOf8bhQ== [fileUNF]Dataset of cortical activity recorded with high spatial resolution from anesthetized ratshdl:21.15109/ARP/0QRUQKhttps://doi.gin.g-node.org/10.12751/g-node.arf7ol/Csaba, HorváthLili Fanni, TóthIstván, UlbertFiáth, Richárd2017-1.2.1-NKP-2017-00002PD124175PD134196TUDFO/51757-1/2019-ITMARPFiáth, RichárdFiáth, Richárd2021-07-08Medicine, Health and Life SciencesSilicon probeHigh density recordingHigh spatial resolutionElectrophysiologyNeurophysiologyNeural recordingExtracellular action potentialSingle unit activitySpike sortingSlow oscillationSlow wave activityKetamine/xylazine anesthesiaNeocortexfakeDoiMigrationDoneHere, we present an electrophysiological dataset (~0.9 TB overall size) recorded from the neocortex of twenty rats anesthetized with ketamine/xylazine (<a href="https://doi.org/10.1038/s41597-021-00970-3">Horváth et al., 2021</a>). The wideband, spontaneous recordings (n = 109) were acquired with a single-shank silicon-based probe having 128 densely packed recording sites arranged in a 32x4 array (<a href="https://doi.org/10.1016/j.bios.2018.01.060">Fiáth et al., 2018</a>). The dataset contains the activity of a total of 7126 sorted single units extracted from all layers of the cortex. Here, we share raw neural recordings, as well as spike times, extracellular spike waveforms and several properties of units packaged in a standardized electrophysiological data format. For technical validation of our dataset, we provide the distributions of derived single unit properties along with various spike sorting quality metrics.high-density electrophysiological recordingsexperimental data<h3>Repository structure</h3> 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 identifier of the animal (from Rat01 to Rat20), probe insertion sequence (Insertion1 or Insertion2) and cortical depth (from Depth1 to Depth3). The filename of the NWB file (identifier) was constructed by concatenating the above information (e.g., Rat01_Insertion2_Depth3). 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, the average signal-to-noise ratio of single units, the degree of power line (50 Hz) noise contamination (given by the power spectral density value measured at 50 Hz, and also by the ratio of the power spectral density measured at 50 Hz to the power spectral density measured at 49 Hz), the RMS noise and RMS signal levels. <h3>NWB file structure</h3> Each NWB file contains several main groups which are similar to directories. The acquisition group contains the continuous wideband 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 general 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 and laminar location 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 units 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’). For an easier visualization of the multichannel spike waveform in two dimensions, we have also added an array which contains the mean spike waveform in the 32 x 4 shape of the electrode array (‘mean_waveform_all_channels_filt_32x4’; ’mean_waveform_all_channels_raw_32x4’). 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 units 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 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. <h3>Interacting with the NWB files</h3> 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. Loading a short segment (20.000 samples corresponding to 1 second of data) of the raw wideband recording on all (128) channels: <code> 1. nwb = nwbRead('Rat01_Insertion1_Depth1.nwb'); </code> <code> 2. dataChunk=nwb.acquisition.get('wideband_multichannel_recording').data.load([1, 1], [128, 20000]); </code> 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. Loading and plotting the mean spike waveform of a specific single unit on the peak waveform channel: <code> 3. peakChannels = nwb.units.vectordata.get('peak_waveform_channel').data.load(); </code> <code> 4. meanWaveforms = nwb.units.vectordata.get('mean_waveform_all_channels_filt').data.load(); </code> <code> 5. mySingleUnit = 11; </code> <code> 6. singleUnitWaveform = meanWaveforms(peakChannels(mySingleUnit), :, mySingleUnit); </code> <code> 7. plot(singleUnitWaveform); </code> Loading the isolation distance quality metric of all units found in a single NWB file: <code> 8. IDvalues = nwb.units.vectordata.get('isolation_distance').data.load(); </code> 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>). We can load the spikes times of a specific single unit (in seconds) the following way: <code> 9. allSpikeTimes = nwb.units.spike_times.data.load(); </code> <code> 10. spikeTimesIndex = nwb.units.spike_times_index.data.load(); </code> <code> 11. spikesOfSingleUnit2 = allSpikeTimes(spikeTimesIndex(1)+1 : spikeTimesIndex(2)); </code> SpikeInterface can also be used to load the wideband data and single unit properties (in Python, works only with version 0.13): <code> 1. import spikeextractors as se </code> <code> 2. nwbPath = 'Rat01_Insertion1_Depth1.nwb' </code> <code> 3. recording = se.NwbRecordingExtractor(nwbPath) </code> <code> 4. sorting = se.NwbSortingExtractor(nwbPath) </code> <code> 5. mySingleUnit = 2 </code> <code> 6. sorting.get_unit_property(mySingleUnit,'isolation_distance') </code> The dataset is also available at the G-Node GIN repository under the following link and DOI: <a href="https://gin.g-node.org/UlbertLab/High_Resolution_Cortical_Spikes">https://gin.g-node.org/UlbertLab/High_Resolution_Cortical_Spikes</a> <a href="https://doi.gin.g-node.org/10.12751/g-node.arf7ol/">https://doi.gin.g-node.org/10.12751/g-node.arf7ol/</a> 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.Fiáth R, Raducanu BC, Musa S, Andrei A, Lopez CM, van Hoof C, Ruther P, Aarts A, Horváth D, Ulbert I. A silicon-based neural probe with densely-packed low-impedance titanium nitride microelectrodes for ultrahigh-resolution in vivo recordings. (2018) BIOSENSORS & BIOELECTRONICS 106: 86-92. <a href="https://doi.org/10.1016/j.bios.2018.01.060">doi: 10.1016/j.bios.2018.01.060</a> Horváth C, Tóth LF, Ulbert I, Fiáth R. Dataset of cortical activity recorded with high spatial resolution from anesthetized rat (2021) SCIENTIFIC DATA 8:180. <a href="https://doi.org/10.1038/s41597-021-00970-3">doi: 10.1038/s41597-021-00970-3</a>10.1038/s41597-021-00970-3Fiáth R, Raducanu BC, Musa S, Andrei A, Lopez CM, van Hoof C, Ruther P, Aarts A, Horváth D, Ulbert I. A silicon-based neural probe with densely-packed low-impedance titanium nitride microelectrodes for ultrahigh-resolution in vivo recordings. (2018) BIOSENSORS & BIOELECTRONICS 106: 86-92. <a href="https://doi.org/10.1016/j.bios.2018.01.060">doi: 10.1016/j.bios.2018.01.060</a> Horváth C, Tóth LF, Ulbert I, Fiáth R. Dataset of cortical activity recorded with high spatial resolution from anesthetized rat (2021) SCIENTIFIC DATA 8:180. <a href="https://doi.org/10.1038/s41597-021-00970-3">doi: 10.1038/s41597-021-00970-3</a>Animal_characteristics_and_targeted_cortical_areas.tab209text/tab-separated-valuesUNF:6:xZZKbMKOGQoqMsQ0rAwr5A==Recording_characteristics.tab10910text/tab-separated-valuesUNF:6:TdSCBxj7yCmd36nx24Uy2Q==Animal IDUNF:6:B9UEOIlDI7dKqzOjS+3Kkw==Weight (g)59.60174847311639200.0315.00.020.0.315.5400.0UNF:6:3Wy26x1d95hwxdYqS7SqAw==SexUNF:6:vfJDWT9dDDWH+0bQ13C9pg==Cortical Area (Insertion1)UNF:6:RCpD1Vjx55675cmKfwPXmQ==AP coordinate (Insertion1)-3.620.0-2.871.-2.84999999999999960.00.4822087340124552-1.95UNF:6:s+2qS7UwmmwnKu3SxkWXoQ== ML coordinate (Insertion1)2.37820.02.961.95.0.25769422925140110.02.385UNF:6:J/RVbGkABi+gKGRLaiTNug== Cortical Area (Insertion2)UNF:6:vb3fMGJ7GZ2zitKw5rAqlw== AP coordinate (Insertion2)UNF:6:3z9ipommmAuisnAO/KaLgQ== ML coordinate (Insertion2)UNF:6:qU5I3cdo+bXI54+LelEV5Q==NWB file nameUNF:6:n27AEH6FbWee1lOCWNmcwA==NWB file size (GB)109.013.450.02.4111343932227113.148.135779816513761.6.99UNF:6:RHhSY5IKUftEh/673CfRPQ==Recording duration (min)37.55522935779816.61.9415.130.0109.011.2879017172541631.78UNF:6:PbsJYJvnNsQsPZ/CWFBzEA==Cortical areaUNF:6:tlikml6icW1Bq7w+zkD2TA==Single unit yield65.37614678899082.109.056.0159.02.040.9843667103069860.0UNF:6:4zDjOqSr/8mQdRo207P0Bw==Average signal-to-noise ratio of single units28.5813.72109.00.06.4613.9554128440366973.4722587194819865.UNF:6:o+xKFHmP/HluoQW4mDUCRw==Power spectral density at 50 Hz (dB/Hz)22.6532.34495412844036431.250.0109.0.47.395.1343297416531835UNF:6:RnwdjQdMRgMoowrxnNjJ/w==50 Hz/49 Hz ratio1.010.01.960.20883314154074398109.01.251.2971559633027523.UNF:6:6D0D/CPjThYw0u+60OArmA==RMS noise level (300 - 6000 Hz)9.9.4.133.410.0109.01.27909519020537444.4722935779816515UNF:6:rHLfFBVXa4WawwV8VobW7w==RMS signal level (300 - 6000 Hz)109.05.03324275897200710.8512.11733944954128431.145.17.0.0UNF:6:xy1wh+XEFQxa4PphqF0TvQ==allSingleUnits.nwbapplication/octet-streamNWB_tutorial_scipt.mtext/plainRat01_Insertion1_Depth1.nwbapplication/octet-streamRat01_Insertion1_Depth2.nwbapplication/octet-streamRat01_Insertion1_Depth3.nwbapplication/octet-streamRat01_Insertion2_Depth1.nwbapplication/octet-streamRat01_Insertion2_Depth2.nwbapplication/octet-streamRat01_Insertion2_Depth3.nwbapplication/octet-streamRat02_Insertion1_Depth1.nwbapplication/octet-streamRat02_Insertion1_Depth2.nwbapplication/octet-streamRat02_Insertion2_Depth1.nwbapplication/octet-streamRat02_Insertion2_Depth2.nwbapplication/octet-streamRat03_Insertion1_Depth1.nwbapplication/octet-streamRat03_Insertion1_Depth2.nwbapplication/octet-streamRat03_Insertion1_Depth3.nwbapplication/octet-streamRat03_Insertion2_Depth1.nwbapplication/octet-streamRat03_Insertion2_Depth2.nwbapplication/octet-streamRat04_Insertion1_Depth1.nwbapplication/octet-streamRat04_Insertion1_Depth2.nwbapplication/octet-streamRat04_Insertion1_Depth3.nwbapplication/octet-streamRat04_Insertion2_Depth1.nwbapplication/octet-streamRat04_Insertion2_Depth2.nwbapplication/octet-streamRat05_Insertion1_Depth1.nwbapplication/octet-streamRat05_Insertion1_Depth2.nwbapplication/octet-streamRat05_Insertion1_Depth3.nwbapplication/octet-streamRat05_Insertion2_Depth1.nwbapplication/octet-streamRat05_Insertion2_Depth2.nwbapplication/octet-streamRat06_Insertion1_Depth1.nwbapplication/octet-streamRat06_Insertion1_Depth2.nwbapplication/octet-streamRat06_Insertion1_Depth3.nwbapplication/octet-streamRat06_Insertion2_Depth1.nwbapplication/octet-streamRat06_Insertion2_Depth2.nwbapplication/octet-streamRat06_Insertion2_Depth3.nwbapplication/octet-streamRat07_Insertion1_Depth1.nwbapplication/octet-streamRat07_Insertion1_Depth2.nwbapplication/octet-streamRat07_Insertion1_Depth3.nwbapplication/octet-streamRat08_Insertion1_Depth1.nwbapplication/octet-streamRat08_Insertion1_Depth2.nwbapplication/octet-streamRat08_Insertion1_Depth3.nwbapplication/octet-streamRat08_Insertion2_Depth1.nwbapplication/octet-streamRat08_Insertion2_Depth2.nwbapplication/octet-streamRat08_Insertion2_Depth3.nwbapplication/octet-streamRat09_Insertion1_Depth1.nwbapplication/octet-streamRat09_Insertion1_Depth2.nwbapplication/octet-streamRat09_Insertion1_Depth3.nwbapplication/octet-streamRat09_Insertion2_Depth1.nwbapplication/octet-streamRat09_Insertion2_Depth2.nwbapplication/octet-streamRat10_Insertion1_Depth1.nwbapplication/octet-streamRat10_Insertion1_Depth2.nwbapplication/octet-streamRat10_Insertion1_Depth3.nwbapplication/octet-streamRat10_Insertion2_Depth1.nwbapplication/octet-streamRat10_Insertion2_Depth2.nwbapplication/octet-streamRat10_Insertion2_Depth3.nwbapplication/octet-streamRat11_Insertion1_Depth1.nwbapplication/octet-streamRat11_Insertion1_Depth2.nwbapplication/octet-streamRat11_Insertion1_Depth3.nwbapplication/octet-streamRat11_Insertion2_Depth1.nwbapplication/octet-streamRat11_Insertion2_Depth2.nwbapplication/octet-streamRat12_Insertion1_Depth1.nwbapplication/octet-streamRat12_Insertion1_Depth2.nwbapplication/octet-streamRat12_Insertion2_Depth1.nwbapplication/octet-streamRat12_Insertion2_Depth2.nwbapplication/octet-streamRat13_Insertion1_Depth1.nwbapplication/octet-streamRat13_Insertion1_Depth2.nwbapplication/octet-streamRat13_Insertion1_Depth3.nwbapplication/octet-streamRat13_Insertion2_Depth1.nwbapplication/octet-streamRat13_Insertion2_Depth2.nwbapplication/octet-streamRat13_Insertion2_Depth3.nwbapplication/octet-streamRat14_Insertion1_Depth1.nwbapplication/octet-streamRat14_Insertion1_Depth2.nwbapplication/octet-streamRat14_Insertion1_Depth3.nwbapplication/octet-streamRat14_Insertion2_Depth1.nwbapplication/octet-streamRat14_Insertion2_Depth2.nwbapplication/octet-streamRat14_Insertion2_Depth3.nwbapplication/octet-streamRat15_Insertion1_Depth1.nwbapplication/octet-streamRat15_Insertion1_Depth2.nwbapplication/octet-streamRat15_Insertion1_Depth3.nwbapplication/octet-streamRat15_Insertion2_Depth1.nwbapplication/octet-streamRat15_Insertion2_Depth2.nwbapplication/octet-streamRat15_Insertion2_Depth3.nwbapplication/octet-streamRat16_Insertion1_Depth1.nwbapplication/octet-streamRat16_Insertion1_Depth2.nwbapplication/octet-streamRat16_Insertion1_Depth3.nwbapplication/octet-streamRat16_Insertion1_Depth4.nwbapplication/octet-streamRat16_Insertion2_Depth1.nwbapplication/octet-streamRat16_Insertion2_Depth2.nwbapplication/octet-streamRat16_Insertion2_Depth3.nwbapplication/octet-streamRat17_Insertion1_Depth1.nwbapplication/octet-streamRat17_Insertion1_Depth2.nwbapplication/octet-streamRat17_Insertion1_Depth3.nwbapplication/octet-streamRat17_Insertion2_Depth1.nwbapplication/octet-streamRat17_Insertion2_Depth2.nwbapplication/octet-streamRat17_Insertion2_Depth3.nwbapplication/octet-streamRat18_Insertion1_Depth1.nwbapplication/octet-streamRat18_Insertion1_Depth2.nwbapplication/octet-streamRat18_Insertion1_Depth3.nwbapplication/octet-streamRat18_Insertion2_Depth1.nwbapplication/octet-streamRat18_Insertion2_Depth2.nwbapplication/octet-streamRat18_Insertion2_Depth3.nwbapplication/octet-streamRat19_Insertion1_Depth1.nwbapplication/octet-streamRat19_Insertion1_Depth2.nwbapplication/octet-streamRat19_Insertion1_Depth3.nwbapplication/octet-streamRat19_Insertion2_Depth1.nwbapplication/octet-streamRat19_Insertion2_Depth2.nwbapplication/octet-streamRat19_Insertion2_Depth3.nwbapplication/octet-streamRat20_Insertion1_Depth1.nwbapplication/octet-streamRat20_Insertion1_Depth2.nwbapplication/octet-streamRat20_Insertion1_Depth3.nwbapplication/octet-streamRat20_Insertion2_Depth1.nwbapplication/octet-streamRat20_Insertion2_Depth2.nwbapplication/octet-streamRat20_Insertion2_Depth3.nwbapplication/octet-stream