Adaptive common average reference for in vivo multichannel local field potentials

A fully automatic multichannel neural spike sorting algorithm with spike reduction and positional feature

Objective: The sorting of neural spike data recorded by multichannel and high channel neural probes such as Neuropixels, especially in real-time, remains a significant technical challenge. Most neural spike sorting algorithms focus on sorting neural spikes post-hoc for high sorting accuracy-but reducing the processing delay for fast sorting, potentially even live sorting, is generally not possible with these algorithms.Approach: Here we report our Graph nEtwork Multichannel sorting (GEMsort) algorithm, which is largely based on graph network, to allow rapid neural spike sorting for multiple neural recording channels. This was accomplished by two innovations: In GEMsort, duplicated neural spikes recorded from multiple channels were eliminated from duplicate channels by only selecting the highest amplitude neural spike in any channel for subsequent processing. In addition, the channel from which the representative neural spike was recorded was used as an additional feature to differentiate between neural spikes recorded from different neurons having similar temporal features.Main results: Synthetic and experimentally recorded multichannel neural recordings were used to evaluate the sorting performance of GEMsort. The sorting results of GEMsort were also compared with two other state-of-the-art sorting algorithms (Kilosort and Mountainsort) in sorting time and sorting agreements.Significance: GEMsort allows rapidly sort neural spikes and is highly suitable to be implemented with digital circuitry for high processing speed and channel scalability.

Dynamic temporal neural patterns based on multichannel LFPs Identify different brain states during anesthesia in pigeons: comparison of three anesthetics

Anesthetic-induced brain activity study is crucial in avian cognitive-, consciousness-, and sleep-related research. However, the neurobiological mechanisms underlying the generation of brain rhythms and specific connectivity of birds during anesthesia are poorly understood. Although different kinds of anesthetics can be used to induce an anesthesia state, a comparison study of these drugs focusing on the neural pattern evolution during anesthesia is lacking. Here, we recorded local field potentials (LFPs) using a multi-channel micro-electrode array inserted into the nidopallium caudolateral (NCL) of adult pigeons (Columba livia) anesthetized with chloral hydrate, pelltobarbitalum natricum or urethane. Power spectral density (PSD) and functional connectivity analyses were used to measure the dynamic temporal neural patterns in NCL during anesthesia. Neural decoding analysis was adopted to calculate the probability of the pigeon's brain state and the kind of injected anesthetic. In the NCL during anesthesia, we found elevated power activity and functional connectivity at low-frequency bands and depressed power activity and connectivity at high-frequency bands. Decoding results based on the spectral and functional connectivity features indicated that the pigeon's brain states during anesthesia and the injected anesthetics can be effectively decoded. These findings provide an important foundation for future investigations on how different anesthetics induce the generation of specific neural patterns.

The Hippocampus in Pigeons Contributes to the Model-Based Valuation and the Relationship between Temporal Context States

Model-based decision-making guides organism behavior by the representation of the relationships between different states. Previous studies have shown that the mammalian hippocampus (Hp) plays a key role in learning the structure of relationships among experiences. However, the hippocampal neural mechanisms of birds for model-based learning have rarely been reported. Here, we trained six pigeons to perform a two-step task and explore whether their Hp contributes to model-based learning. Behavioral performance and hippocampal multi-channel local field potentials (LFPs) were recorded during the task. We estimated the subjective values using a reinforcement learning model dynamically fitted to the pigeon's choice of behavior. The results show that the model-based learner can capture the behavioral choices of pigeons well throughout the learning process. Neural analysis indicated that high-frequency (12-100 Hz) power in Hp represented the temporal context states. Moreover, dynamic correlation and decoding results provided further support for the high-frequency dependence of model-based valuations. In addition, we observed a significant increase in hippocampal neural similarity at the low-frequency band (1-12 Hz) for common temporal context states after learning. Overall, our findings suggest that pigeons use model-based inferences to learn multi-step tasks, and multiple LFP frequency bands collaboratively contribute to model-based learning. Specifically, the high-frequency (12-100 Hz) oscillations represent model-based valuations, while the low-frequency (1-12 Hz) neural similarity is influenced by the relationship between temporal context states. These results contribute to our understanding of the neural mechanisms underlying model-based learning and broaden the scope of hippocampal contributions to avian behavior.

Performance Baseline of Phase Transfer Entropy Methods for Detecting Animal Brain Area Interactions

Objective: Phase transfer entropy (TEθ) methods perform well in animal sensory-spatial associative learning. However, their advantages and disadvantages remain unclear, constraining their usage. Method: This paper proposes the performance baseline of the TEθ methods. Specifically, four TEθ methods are applied to the simulated signals generated by a neural mass model and the actual neural data from ferrets with known interaction properties to investigate the accuracy, stability, and computational complexity of the TEθ methods in identifying the directional coupling. Then, the most suitable method is selected based on the performance baseline and used on the local field potential recorded from pigeons to detect the interaction between the hippocampus (Hp) and nidopallium caudolaterale (NCL) in visual-spatial associative learning. Results: (1) This paper obtains a performance baseline table that contains the most suitable method for different scenarios. (2) The TEθ method identifies an information flow preferentially from Hp to NCL of pigeons at the θ band (4-12 Hz) in visual-spatial associative learning. Significance: These outcomes provide a reference for the TEθ methods in detecting the interactions between brain areas.

Goal given moment modulates the time period of gamma oscillations in nidopallium caudolaterale during the goal-directed behavior of pigeon

The cognitive processes of goal-directed navigation are believed to be organized around and serve the identification and selection of goals. Differences in LFP signals in avian nidopallium caudolaterale (NCL) under different goal location/distance information in the goal-directed behavior have been studied. However, for goals that are multifarious constructs that include various information, the modulation of goal time information on the LFP of NCL during goal-directed behavior remains unclear. In this study, we recorded the LFP activity from the NCL of eight pigeons as they performed two goal-directed decision-making tasks in a plus-maze. During the two tasks with different goal time information, spectral analysis revealed significant LFP power selectively enhanced in the slow gamma band (40-60 Hz), while the slow gamma band of LFP which could effectively decode the behavioral goal of the pigeons existed in different time periods. These findings suggest that the LFP activity in the gamma band correlates with the goal-time information, and help to shed light on the contribution of the gamma rhythm recorded from the NCL in goal-directed behavior.

Different coding characteristics between flight and freezing in dorsal periaqueductal gray of mice during exposure to innate threats

BACKGROUND

Flight and freezing are two vital defensive behaviors that mice display to avoid natural enemies. When they are exposed to innate threats, visual cues are processed and transmitted by the visual system into the emotional nuclei and finally transmitted to the periaqueductal gray (PAG) to induce defensive behaviors. However, how the dorsal PAG (dPAG) encodes the two defensive behaviors is unclear.

METHODS

Multi-array electrodes were implanted in the dPAG nuclei of C57BL/6 mice. Two kinds of visual stimuli (looming and sweeping) were used to induce defensive behaviors in mice. Neural signals under different defense behaviors were recorded, and the encoding characteristics of the two behaviors were extracted and analyzed from spike firing and frequency oscillations. Finally, synchronization of neural activity during the defense process was analyzed.

RESULTS

The neural activity between flight and freezing behaviors showed different firing patterns, and the differences in the inter-spike interval distribution were mainly reflected in the 2-10 ms period. The frequency band activities under both defensive behaviors were concentrated in the theta band; the active frequency of flight was ~8 to 10 Hz, whereas that of freezing behavior was ~6 to 8 Hz. The network connection density under both defense behaviors was significantly higher than the period before and after defensive behavior occurred, indicating that there was a high synchronization of neural activity during the defense process.

CONCLUSIONS

The dPAG nuclei of mice have different coding features between flight and freezing behaviors; during strong looming stimulation, fast neuro-instinctive decision making is required while encountering weak sweeping stimulation, and computable planning late behavior is predicted in the early stage. The frequency band activities under both defensive behaviors were concentrated in the theta band. There was a high synchronization of neural activity during the defense process, which may be a key factor triggering different defensive behaviors.

A Fast and Efficient Ensemble Transfer Entropy and Applications in Neural Signals

The ensemble transfer entropy (TEensemble) refers to the transfer entropy estimated from an ensemble of realizations. Due to its time-resolved analysis, it is adapted to analyze the dynamic interaction between brain regions. However, in the traditional TEensemble, multiple sets of surrogate data should be used to construct the null hypothesis distribution, which dramatically increases the computational complexity. To reduce the computational cost, a fast, efficient TEensemble with a simple statistical test method is proposed here, in which just one set of surrogate data is involved. To validate the improved efficiency, the simulated neural signals are used to compare the characteristics of the novel TEensemble with those of the traditional TEensemble. The results show that the time consumption is reduced by two or three magnitudes in the novel TEensemble. Importantly, the proposed TEensemble could accurately track the dynamic interaction process and detect the strength and the direction of interaction robustly even in the presence of moderate noises. The novel TEensemble reaches its steady state with the increased samples, which is slower than the traditional method. Furthermore, the effectiveness of the novel TEensemble was verified in the actual neural signals. Accordingly, the TEensemble proposed in this work may provide a suitable way to investigate the dynamic interactions between brain regions.

The superior colliculus/lateral posterior thalamic nuclei in mice rapidly transmit fear visual information through the theta frequency band

Animals perceive threat information mainly from vision, and the subcortical visual pathway plays a critical role in the rapid processing of fear visual information. The superior colliculus (SC) and lateral posterior (LP) nuclei of the thalamus are key components of the subcortical visual pathway; however, how animals encode and transmit fear visual information is unclear. To evaluate the response characteristics of neurons in SC and LP thalamic nuclei under fear visual stimuli, extracellular action potentials (spikes) and local field potential signals were recorded under looming and dimming visual stimuli. The results showed that both SC and LP thalamic nuclei were strongly responsive to looming visual stimuli but not sensitive to dimming visual stimuli. Under the looming visual stimulus, the theta (θ) frequency bands of both nuclei showed obvious oscillations, which markedly enhanced the synchronization between neurons. The functional network characteristics also indicated that the network connection density and information transmission efficiency were higher under fear visual stimuli. These findings suggest that both SC and LP thalamic nuclei can effectively identify threatening fear visual information and rapidly transmit it between nuclei through the θ frequency band. This discovery can provide a basis for subsequent coding and decoding studies in the subcortical visual pathways.

Adaptive virtual referencing for the extraction of extracellularly recorded action potentials in noisy environments

OBJECTIVE

Removal of common mode noise and artifacts from extracellularly measured action potentials, herein referred to as spikes, recorded with multi-electrode arrays (MEAs) which included severe noise and artifacts generated by an ultrahigh field (UHF) 16.4 Tesla magnetic resonance imaging (MRI) scanner.

APPROACH

An adaptive virtual referencing (AVR) algorithm is used to remove artifacts and thus enable extraction of neural spike signals from extracellular recordings in anesthetized rat brains. A 16-channel MEA with 150-micron inter-site spacing is used, and a virtual reference is created by spatially averaging the 16 signal channels which results in a reference signal without extracellular spiking activity while preserving common mode noise and artifacts. This virtual reference signal is then used as the input to an adaptive FIR filter which optimally scales and time-shifts the reference to each specific electrode site to remove the artifacts and noise.

MAIN RESULTS

By removing artifacts and reducing noise, the neural spikes at each electrode site can be well extracted, even from data originally recorded with a high noise floor due to electromagnetic interference and artifacts generated by a 16.4T MRI scanner. The AVR method enables many more spikes to be detected than would otherwise be possible. Further, the filtered spike waveforms can be well separated from each other using PCA feature extraction and semi-supervised k-means clustering. While data in a 16.4T MRI scanner contains significantly more noise and artifacts, the developed AVR method enables similar data quality to be extracted as recorded on benchtop experiments outside the MRI scanner.

SIGNIFICANCE

AVR of extracellular spike signals recorded with MEAs has not been previously reported and fills a technical need by enabling low-noise extracellular spike extraction in noisy and challenging environments such as UHF MRI that will enable further study of neuro-vascular coupling at UHF.

Hippocampus-nidopallium caudolaterale interactions exist in the goal-directed behavior of pigeon

Avian hippocampus (Hp) and nidopallium caudolaterale (NCL) are believed to play key roles in goal-directed behavior. However, it is still unclear whether there are interactions between the two brain regions in the goal-directed behavior of pigeons. To investigate the interactions between the Hp and the NCL in the goal-directed behavior, we recorded local field potential (LFP) signals from the two regions simultaneously when the pigeons performed a goal-directed decision-making task. Amplitude-amplitude coupling analysis revealed that the coupling value between the LFP recorded from the Hp and that from the NCL increased significantly (P < 0.05) in slow gamma-band (40-60 Hz) during the turning area. In addition, the LFP functional network analysis demonstrated the LFP functional connections between the Hp and the NCL increased significantly (P < 0.05) in the turning area. The result of partial directed coherence (PDC) analysis showed that the predominant direction of information flow is thought to be from the Hp to the NCL. These findings suggest that there are causal functional interactions between the Hp and the NCL by which information is transmitted between the two regions relevant to goal-directed behavior.

[Comparison of decoding performance between spike and local field potential signals during goal-directed decision-making task of pigeons]

Both spike and local field potential (LFP) signals are two of the most important candidate signals for neural decoding. At present there are numerous studies on their decoding performance in mammals, but the decoding performance in birds is still not clear. We analyzed the decoding performance of both signals recorded from nidopallium caudolaterale area in six pigeons during the goal-directed decision-making task using the decoding algorithm combining leave-one-out and k-nearest neighbor (LOO- kNN). And the influence of the parameters, include the number of channels, the position and size of decoding window, and the nearest neighbor k value, on the decoding performance was also studied. The results in this study have shown that the two signals can effectively decode the movement intention of pigeons during the this task, but in contrast, the decoding performance of LFP signal is higher than that of spike signal and it is less affected by the number of channels. The best decoding window is in the second half of the goal-directed decision-making process, and the optimal decoding window size of LFP signal (0.3 s) is shorter than that of spike signal (1 s). For the LOO- kNN algorithm, the accuracy is inversely proportional to the k value. The smaller the k value is, the larger the accuracy of decoding is. The results in this study will help to parse the neural information processing mechanism of brain and also have reference value for brain-computer interface.

Decoding Pigeon Behavior Outcomes Using Functional Connections among Local Field Potentials

Recent studies indicate that the local field potential (LFP) carries information about an animal's behavior, but issues regarding whether there are any relationships between the LFP functional networks and behavior tasks as well as whether it is possible to employ LFP network features to decode the behavioral outcome in a single trial remain unresolved. In this study, we developed a network-based method to decode the behavioral outcomes in pigeons by using the functional connectivity strength values among LFPs recorded from the nidopallium caudolaterale (NCL). In our method, the functional connectivity strengths were first computed based on the synchronization likelihood. Second, the strength values were unwrapped into row vectors and their dimensions were then reduced by principal component analysis. Finally, the behavioral outcomes in single trials were decoded using leave-one-out combined with the k-nearest neighbor method. The results showed that the LFP functional network based on the gamma-band was related to the goal-directed behavior of pigeons. Moreover, the accuracy of the network features (74 ± 8%) was significantly higher than that of the power features (61 ± 12%). The proposed method provides a powerful tool for decoding animal behavior outcomes using a neural functional network.