Details of collection, processing and analysis of recordings are described in the published article. Each recording and corresponding metadata, single unit properties and quality metrics were packaged in the\nNeurodata Without Borders: Neurophysiology version 2.0 (NWB:N 2.0) data format using the NeuroConv Python package. A single NWB file was created from each recording. NWB files were placed in folders based on the brain area (Neocortex or Thalamus), the animal model (Rat, Mouse or Human) and the channel number (e.g., 256, 128, 64, 32 or 16). The filename of the NWB file (identifier) was constructed by concatenating the above information and the identifier of the animal (e.g., Rat01_Neocortex_256channel.nwb). In the case of multiple probe insertions being performed in a single animal, or recordings were carried out at multiple brain depths, this information was also incorporated in the NWB filename (e.g., Rat01_Thalamus_256ch_Insertion1.nwb). The original human Neuropixels recordings can be found in the Dryad repository.
\n\nThe CSV file named \"Animal_characteristics_\" contains information about the subjects (e.g., age, weight, sex). The \"Recording_characteristics\" CSV file lists several useful properties for each NWB file including the file size, the duration of the recording, the cortical or thalamic location and stereotaxic coordinates, the single unit yield or the average signal-to-noise ratio of single units.
\nEach NWB file contains several main groups which are similar to directories. The acquisition group contains the continuous wideband multichannel data (‘ElectricalSeriesRaw’) 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), 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 from the wideband data (‘waveform_mean’; ’waveform_sd’). Several single unit properties (e.g., half_width, unit_name, peak_channel) and cluster quality metrics (e.g., isolation_distance, amplitude_cutoff, presence ratio), as well as the spike times and spike count of each unit were saved in the units group.The structure of NWB files can be explored using the HDFView software.
\n
\n\n
After downloading, the structure of NWB files can be explored using the freely available HDFView software. Users can import data from NWB files offline using the PyNWB and MatNWB APIs, or using SpikeInterface. \nLoaded samples of the raw data have to be multiplied by a conversion number (e.g., 0.195 for multichannel recordings) 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.
\n\n1.\tnwb = nwbRead('Rat01_Neocortex_128channel.nwb');\n\n\n2.\tdataChunk=nwb.acquisition.get('ElectricalSeriesRaw').data.load([1, 1], [128, 20000]);\n\n\n3.\tpeakChannels = nwb.units.vectordata.get('peak_waveform_channel').data.load();\n\n\n4.\tmeanWaveforms = nwb.units.vectordata.get('mean_waveform_all_channels_filt').data.load();\n\n\n5.\tmySingleUnit = 11;\n\n\n6.\tsingleUnitWaveform = meanWaveforms(peakChannels(mySingleUnit), :, mySingleUnit);\n\n\n7.\tplot(singleUnitWaveform);\n\n\n8.\tIDvalues = nwb.units.vectordata.get('isolation_distance').data.load();\n\n\n9.\tallSpikeTimes = nwb.units.spike_times.data.load();\n\n\n10.\tspikeTimesIndex = nwb.units.spike_times_index.data.load();\n\n\n11.\tspikesOfSingleUnit2 = allSpikeTimes(spikeTimesIndex(1)+1 : spikeTimesIndex(2));\n\n\n1.\timport spikeextractors as se\n\n\n2.\tnwbPath = 'Rat01_Neocortex_128channel'\n\n\n3.\trecording = se.NwbRecordingExtractor(nwbPath)\n\n\n4.\tsorting = se.NwbSortingExtractor(nwbPath)\n\n\n5.\tmySingleUnit = 2\n\n\n6.\tsorting.get_unit_property(mySingleUnit,'isolation_distance')\n\n