[To determine the value of iMAX, a new electrodiagnostic method, and its application in the evaluation of peripheral motor nerve excitability].

Objective: To test the new method of iMAX (the minimum stimulus current that elicits the maximum compound muscle action potential amplitude) electrodiagnosis, verify the feasibility of this method in evaluating the excitability of peripheral motor axons, and preliminarily explore the clinical application value. Methods: This study was a cross-sectional study. A total of 50 healthy subjects were recruited from the outpatient department of Peking University Third Hospital from June 2022 to March 2023, including 25 males and 25 females, aged 25-68 (48±8) years. Eleven patients with Charcot-Marie-Pain-1A (CMT1A), 7 males and 4 females, aged 19-55 (41±13) years and 21 patients with diabetic peripheral neuropathy (DPN), 10 males and 11 females, aged 28-79 (53±16) years were enrolled in this study. iMAX of bilateral median nerves, ulnar nerves and peroneal nerves were detected in all patients. Repeatable motor responses with minimum motor threshold and amplitude of at least 0.1 mV and the minimum stimulus current intensity, at which the maximum compound muscle action potential amplitude is elicited, were measured respectively [1 mA increment is called (iUP) and, 0.1 mA adjustment is called (iMAX)].Comparison of the parameters: the parameters of threshold, iUP and iMAX were compared among different age groups, genders and sides, body mass index(BMI) values and detection time , as well as between CMT1A patients, DPN patients and healthy subjects. Results: In healthy subjects, the threshold, iUP value and iMAX value were (1.8±0.7) mA, (4.4±1.2) mA, and (4.2±1.3) mA respectively; ulnar nerve (3.1±1.6) mA, (6.8±3.2) mA, (6.4±3.2) mA; peroneal nerve (3.7±2.0) mA, (7.8±2.8) mA, (7.4±2.9) mA. There were statistically significant differences in threshold, iUP value and iMAX value among different age groups (all P<0.001).With the increase of age, there was a trend of increasing threshold, iUP, and iMAX values in different nerves, and the differences are statistically significant (all P<0.001). There were no significant differences in gender, side and detection time threshold, iUP value and iMAX value (all P>0.05). The parameters of healthy subjects with high BMI value were higher than those of healthy subjects with low BMI value(all P<0.05). Compared with the healthy subjects, the parameters of 11 CMT1A patients were significantly increased (all P<0.05), and the parameters of 21 DPN patients were slightly increased (P<0.05). Conclusion: The new iMAX method reflects the excitability of motor axons and early axonal dysfunction, which is an important supplement to the traditional nerve conduction, and can be used to monitor motor axon excitability disorders.

Neural Drive and Motor Unit Characteristics of the Serratus Anterior in Individuals With Scapular Dyskinesis.

OBJECTIVE: Scapular dyskinesis is one of the causes of shoulder disorders and involves muscle weakness in the serratus anterior. This study investigated whether motor unit (MU) recruitment and firing property, which are important for muscle exertion, have altered in serratus anterior of the individuals with scapular dyskinesis. METHODS: Asymptomatic adults with (SD) and without (control) scapular dyskinesis were analyzed. Surface electromyography (sEMG) waveforms were collected at submaximal voluntary contraction of the serratus anterior. The sEMG waveform was decomposed into MU action potential amplitude (MUAPAMP), mean firing rate (MFR), and recruitment threshold. MUs were divided into low, moderate, and high thresholds, and MU recruitment and firing properties of the groups were compared. RESULTS: High-threshold MUAPAMP was significantly smaller in the SD group than in the control group. The control group also exhibited recruitment properties that reflected the size principle, however, the SD group did not. Furthermore, the SD group had a lower MFR than the control group. CONCLUSIONS: Individuals with scapular dyskinesis exhibit altered MU recruitment properties and lower firing rates of the serratus anterior; this may be detrimental to muscle performance. Thus, it may be necessary to improve the neural drive of the serratus anterior when correcting scapular dyskinesis.

Months-long tracking of neuronal ensembles spanning multiple brain areas with Ultra-Flexible Tentacle Electrodes.

We introduce Ultra-Flexible Tentacle Electrodes (UFTEs), packing many independent fibers with the smallest possible footprint without limitation in recording depth using a combination of mechanical and chemical tethering for insertion. We demonstrate a scheme to implant UFTEs simultaneously into many brain areas at arbitrary locations without angle-of-insertion limitations, and a 512-channel wireless logger. Immunostaining reveals no detectable chronic tissue damage even after several months. Mean spike signal-to-noise ratios are 1.5-3x compared to the state-of-the-art, while the highest signal-to-noise ratios reach 89, and average cortical unit yields are ~1.75/channel. UFTEs can track the same neurons across sessions for at least 10 months (longest duration tested). We tracked inter- and intra-areal neuronal ensembles (neurons repeatedly co-activated within 25 ms) simultaneously from hippocampus, retrosplenial cortex, and medial prefrontal cortex in freely moving rodents. Average ensemble lifetimes were shorter than the durations over which we can track individual neurons. We identify two distinct classes of ensembles. Those tuned to sharp-wave ripples display the shortest lifetimes, and the ensemble members are mostly hippocampal. Yet, inter-areal ensembles with members from both hippocampus and cortex have weak tuning to sharp wave ripples, and some have unusual months-long lifetimes. Such inter-areal ensembles occasionally remain inactive for weeks before re-emerging.

Two Distinct Neuronal Populations in the Rat Parafascicular Nucleus Oppositely Encode the Engagement in Stimulus-driven Reward-seeking.

BACKGROUND: The thalamus is a phylogenetically well-preserved structure. Known to densely contact cortical regions, its role in the transmission of sensory information to the striatal complex has been widely reconsidered in recent years. METHODS: The parafascicular nucleus of the thalamus (Pf) has been implicated in the orientation of attention toward salient sensory stimuli. In a stimulus-driven reward-seeking task, we sought to characterize the electrophysiological activity of Pf neurons in rats. RESULTS: We observed a predominance of excitatory over inhibitory responses for all events in the task. Neurons responded more strongly to the stimulus compared to lever-pressing and reward collecting, confirming the strong involvement of the Pf in sensory information processing. The use of long sessions allowed us to compare neuronal responses to stimuli between trials when animals were engaged in action and those when they were not. We distinguished two populations of neurons with opposite responses: MOTIV+ neurons responded more intensely to stimuli followed by a behavioral response than those that were not. Conversely, MOTIV- neurons responded more strongly when the animal did not respond to the stimulus. In addition, the latency of excitation of MOTIV- neurons was shorter than that of MOTIV+ neurons. CONCLUSION: Through this encoding, the Pf could perform an early selection of environmental stimuli transmitted to the striatum according to motivational level.

The impact of spike timing precision and spike emission reliability on decoding accuracy.

Precisely timed and reliably emitted spikes are hypothesized to serve multiple functions, including improving the accuracy and reproducibility of encoding stimuli, memories, or behaviours across trials. When these spikes occur as a repeating sequence, they can be used to encode and decode a potential time series. Here, we show both analytically and in simulations that the error incurred in approximating a time series with precisely timed and reliably emitted spikes decreases linearly with the number of neurons or spikes used in the decoding. This was verified numerically with synthetically generated patterns of spikes. Further, we found that if spikes were imprecise in their timing, or unreliable in their emission, the error incurred in decoding with these spikes would be sub-linear. However, if the spike precision or spike reliability increased with network size, the error incurred in decoding a time-series with sequences of spikes would maintain a linear decrease with network size. The spike precision had to increase linearly with network size, while the probability of spike failure had to decrease with the square-root of the network size. Finally, we identified a candidate circuit to test this scaling relationship: the repeating sequences of spikes with sub-millisecond precision in area HVC (proper name) of the zebra finch. This scaling relationship can be tested using both neural data and song-spectrogram-based recordings while taking advantage of the natural fluctuation in HVC network size due to neurogenesis.

Anti-Hebbian plasticity drives sequence learning in striatum.

Spatio-temporal activity patterns have been observed in a variety of brain areas in spontaneous activity, prior to or during action, or in response to stimuli. Biological mechanisms endowing neurons with the ability to distinguish between different sequences remain largely unknown. Learning sequences of spikes raises multiple challenges, such as maintaining in memory spike history and discriminating partially overlapping sequences. Here, we show that anti-Hebbian spike-timing dependent plasticity (STDP), as observed at cortico-striatal synapses, can naturally lead to learning spike sequences. We design a spiking model of the striatal output neuron receiving spike patterns defined as sequential input from a fixed set of cortical neurons. We use a simple synaptic plasticity rule that combines anti-Hebbian STDP and non-associative potentiation for a subset of the presented patterns called rewarded patterns. We study the ability of striatal output neurons to discriminate rewarded from non-rewarded patterns by firing only after the presentation of a rewarded pattern. In particular, we show that two biological properties of striatal networks, spiking latency and collateral inhibition, contribute to an increase in accuracy, by allowing a better discrimination of partially overlapping sequences. These results suggest that anti-Hebbian STDP may serve as a biological substrate for learning sequences of spikes.

An artificial visual neuron with multiplexed rate and time-to-first-spike coding.

Human visual neurons rely on event-driven, energy-efficient spikes for communication, while silicon image sensors do not. The energy-budget mismatch between biological systems and machine vision technology has inspired the development of artificial visual neurons for use in spiking neural network (SNN). However, the lack of multiplexed data coding schemes reduces the ability of artificial visual neurons in SNN to emulate the visual perception ability of biological systems. Here, we present an artificial visual spiking neuron that enables rate and temporal fusion (RTF) coding of external visual information. The artificial neuron can code visual information at different spiking frequencies (rate coding) and enables precise and energy-efficient time-to-first-spike (TTFS) coding. This multiplexed sensory coding scheme could improve the computing capability and efficacy of artificial visual neurons. A hardware-based SNN with the RTF coding scheme exhibits good consistency with real-world ground truth data and achieves highly accurate steering and speed predictions for self-driving vehicles in complex conditions. The multiplexed RTF coding scheme demonstrates the feasibility of developing highly efficient spike-based neuromorphic hardware.

Drift of neural ensembles driven by slow fluctuations of intrinsic excitability.

Representational drift refers to the dynamic nature of neural representations in the brain despite the behavior being seemingly stable. Although drift has been observed in many different brain regions, the mechanisms underlying it are not known. Since intrinsic neural excitability is suggested to play a key role in regulating memory allocation, fluctuations of excitability could bias the reactivation of previously stored memory ensembles and therefore act as a motor for drift. Here, we propose a rate-based plastic recurrent neural network with slow fluctuations of intrinsic excitability. We first show that subsequent reactivations of a neural ensemble can lead to drift of this ensemble. The model predicts that drift is induced by co-activation of previously active neurons along with neurons with high excitability which leads to remodeling of the recurrent weights. Consistent with previous experimental works, the drifting ensemble is informative about its temporal history. Crucially, we show that the gradual nature of the drift is necessary for decoding temporal information from the activity of the ensemble. Finally, we show that the memory is preserved and can be decoded by an output neuron having plastic synapses with the main region.

Spiking Neural Membrane Systems with Adaptive Synaptic Time Delay.

Spiking neural membrane systems (or spiking neural P systems, SNP systems) are a new type of computation model which have attracted the attention of plentiful scholars for parallelism, time encoding, interpretability and extensibility. The original SNP systems only consider the time delay caused by the execution of rules within neurons, but not caused by the transmission of spikes via synapses between neurons and its adaptive adjustment. In view of the importance of time delay for SNP systems, which are a time encoding computation model, this study proposes SNP systems with adaptive synaptic time delay (ADSNP systems) based on the dynamic regulation mechanism of synaptic transmission delay in neural systems. In ADSNP systems, besides neurons, astrocytes that can generate adenosine triphosphate (ATP) are introduced. After receiving spikes, astrocytes convert spikes into ATP and send ATP to the synapses controlled by them to change the synaptic time delays. The Turing universality of ADSNP systems in number generating and accepting modes is proved. In addition, a small universal ADSNP system using 93 neurons and astrocytes is given. The superiority of the ADSNP system is demonstrated by comparison with the six variants. Finally, an ADSNP system is constructed for credit card fraud detection, which verifies the feasibility of the ADSNP system for solving real-world problems. By considering the adaptive synaptic delay, ADSNP systems better restore the process of information transmission in biological neural networks, and enhance the adaptability of SNP systems, making the control of time more accurate.

Sparse Spiking Neural-Like Membrane Systems on Graphics Processing Units.

The parallel simulation of Spiking Neural P systems is mainly based on a matrix representation, where the graph inherent to the neural model is encoded in an adjacency matrix. The simulation algorithm is based on a matrix-vector multiplication, which is an operation efficiently implemented on parallel devices. However, when the graph of a Spiking Neural P system is not fully connected, the adjacency matrix is sparse and hence, lots of computing resources are wasted in both time and memory domains. For this reason, two compression methods for the matrix representation were proposed in a previous work, but they were not implemented nor parallelized on a simulator. In this paper, they are implemented and parallelized on GPUs as part of a new Spiking Neural P system with delays simulator. Extensive experiments are conducted on high-end GPUs (RTX2080 and A100 80GB), and it is concluded that they outperform other solutions based on state-of-the-art GPU libraries when simulating Spiking Neural P systems.

Statistically inferred neuronal connections in subsampled neural networks strongly correlate with spike train covariances.

Statistically inferred neuronal connections from observed spike train data are often skewed from ground truth by factors such as model mismatch, unobserved neurons, and limited data. Spike train covariances, sometimes referred to as "functional connections," are often used as a proxy for the connections between pairs of neurons, but reflect statistical relationships between neurons, not anatomical connections. Moreover, covariances are not causal: spiking activity is correlated in both the past and the future, whereas neurons respond only to synaptic inputs in the past. Connections inferred by maximum likelihood inference, however, can be constrained to be causal. However, we show in this work that the inferred connections in spontaneously active networks modeled by stochastic leaky integrate-and-fire networks strongly correlate with the covariances between neurons, and may reflect noncausal relationships, when many neurons are unobserved or when neurons are weakly coupled. This phenomenon occurs across different network structures, including random networks and balanced excitatory-inhibitory networks. We use a combination of simulations and a mean-field analysis with fluctuation corrections to elucidate the relationships between spike train covariances, inferred synaptic filters, and ground-truth connections in partially observed networks.

Optimal input reverberation and homeostatic self-organization toward the edge of synchronization.

Transient or partial synchronization can be used to do computations, although a fully synchronized network is sometimes related to the onset of epileptic seizures. Here, we propose a homeostatic mechanism that is capable of maintaining a neuronal network at the edge of a synchronization transition, thereby avoiding the harmful consequences of a fully synchronized network. We model neurons by maps since they are dynamically richer than integrate-and-fire models and more computationally efficient than conductance-based approaches. We first describe the synchronization phase transition of a dense network of neurons with different tonic spiking frequencies coupled by gap junctions. We show that at the transition critical point, inputs optimally reverberate through the network activity through transient synchronization. Then, we introduce a local homeostatic dynamic in the synaptic coupling and show that it produces a robust self-organization toward the edge of this phase transition. We discuss the potential biological consequences of this self-organization process, such as its relation to the Brain Criticality hypothesis, its input processing capacity, and how its malfunction could lead to pathological synchronization and the onset of seizure-like activity.

Chronic pain enhances excitability of corticotropin-releasing factor-expressing neurons in the oval part of the bed nucleus of the stria terminalis.

We previously reported that enhanced corticotropin-releasing factor (CRF) signaling in the bed nucleus of the stria terminalis (BNST) caused the aversive responses during acute pain and suppressed the brain reward system during chronic pain. However, it remains to be examined whether chronic pain alters the excitability of CRF neurons in the BNST. In this study we investigated the chronic pain-induced changes in excitability of CRF-expressing neurons in the oval part of the BNST (ovBNSTCRF neurons) by whole-cell patch-clamp electrophysiology. CRF-Cre; Ai14 mice were used to visualize CRF neurons by tdTomato. Electrophysiological recordings from brain slices prepared from a mouse model of neuropathic pain revealed that rheobase and firing threshold were significantly decreased in the chronic pain group compared with the sham-operated control group. Firing rate of the chronic pain group was higher than that of the control group. These data indicate that chronic pain elevated neuronal excitability of ovBNSTCRF neurons.

Spiking activity in the visual thalamus is coupled to pupil dynamics across temporal scales.

The processing of sensory information, even at early stages, is influenced by the internal state of the animal. Internal states, such as arousal, are often characterized by relating neural activity to a single "level" of arousal, defined by a behavioral indicator such as pupil size. In this study, we expand the understanding of arousal-related modulations in sensory systems by uncovering multiple timescales of pupil dynamics and their relationship to neural activity. Specifically, we observed a robust coupling between spiking activity in the mouse dorsolateral geniculate nucleus (dLGN) of the thalamus and pupil dynamics across timescales spanning a few seconds to several minutes. Throughout all these timescales, 2 distinct spiking modes-individual tonic spikes and tightly clustered bursts of spikes-preferred opposite phases of pupil dynamics. This multi-scale coupling reveals modulations distinct from those captured by pupil size per se, locomotion, and eye movements. Furthermore, coupling persisted even during viewing of a naturalistic movie, where it contributed to differences in the encoding of visual information. We conclude that dLGN spiking activity is under the simultaneous influence of multiple arousal-related processes associated with pupil dynamics occurring over a broad range of timescales.

A longitudinal electrophysiological and behavior dataset for PD rat in response to deep brain stimulation.

Here we presented an electrophysiological dataset collected from layer V of the primary motor cortex (M1) and the corresponding behavior dataset from normal and hemi-parkinson rats over 5 consecutive weeks. The electrophysiological dataset was constituted by the raw wideband signal, neuronal spikes, and local field potential (LFP) signal. The open-field test was done and recorded to evaluate the behavior variation of rats among the entire experimental cycle. We conducted technical validation of this dataset through sorting the spike data to form action potential waveforms and analyzing the spectral power of LFP data, then based on these findings a closed-loop DBS protocol was developed by the oscillation activity response of M1 LFP signal. Additionally, this protocol was applied to the hemi-parkinson rat for five consecutive days while simultaneously recording the electrophysiological data. This dataset is currently the only publicly available dataset that includes longitudinal closed-loop DBS recordings, which can be utilized to investigate variations of neuronal activity within the M1 following long-term closed-loop DBS, and explore additional reliable biomarkers.

Relationship between ion currents and membrane capacitance in canine ventricular myocytes.

Current density, the membrane current value divided by membrane capacitance (Cm), is widely used in cellular electrophysiology. Comparing current densities obtained in different cell populations assume that Cm and ion current magnitudes are linearly related, however data is scarce about this in cardiomyocytes. Therefore, we statistically analyzed the distributions, and the relationship between parameters of canine cardiac ion currents and Cm, and tested if dividing original parameters with Cm had any effect. Under conventional voltage clamp conditions, correlations were high for IK1, moderate for IKr and ICa,L, while negligible for IKs. Correlation between Ito1 peak amplitude and Cm was negligible when analyzing all cells together, however, the analysis showed high correlations when cells of subepicardial, subendocardial or midmyocardial origin were analyzed separately. In action potential voltage clamp experiments IK1, IKr and ICa,L parameters showed high correlations with Cm. For INCX, INa,late and IKs there were low-to-moderate correlations between Cm and these current parameters. Dividing the original current parameters with Cm reduced both the coefficient of variation, and the deviation from normal distribution. The level of correlation between ion currents and Cm varies depending on the ion current studied. This must be considered when evaluating ion current densities in cardiac cells.

Learning with filopodia and spines: Complementary strong and weak competition lead to specialized, graded, and protected receptive fields.

Filopodia are thin synaptic protrusions that have been long known to play an important role in early development. Recently, they have been found to be more abundant in the adult cortex than previously thought, and more plastic than spines (button-shaped mature synapses). Inspired by these findings, we introduce a new model of synaptic plasticity that jointly describes learning of filopodia and spines. The model assumes that filopodia exhibit strongly competitive learning dynamics -similarly to additive spike-timing-dependent plasticity (STDP). At the same time it proposes that, if filopodia undergo sufficient potentiation, they consolidate into spines. Spines follow weakly competitive learning, classically associated with multiplicative, soft-bounded models of STDP. This makes spines more stable and sensitive to the fine structure of input correlations. We show that our learning rule has a selectivity comparable to additive STDP and captures input correlations as well as multiplicative models of STDP. We also show how it can protect previously formed memories and perform synaptic consolidation. Overall, our results can be seen as a phenomenological description of how filopodia and spines could cooperate to overcome the individual difficulties faced by strong and weak competition mechanisms.

The tuning of tuning: How adaptation influences single cell information transfer.

Sensory neurons reconstruct the world from action potentials (spikes) impinging on them. To effectively transfer information about the stimulus to the next processing level, a neuron needs to be able to adapt its working range to the properties of the stimulus. Here, we focus on the intrinsic neural properties that influence information transfer in cortical neurons and how tightly their properties need to be tuned to the stimulus statistics for them to be effective. We start by measuring the intrinsic information encoding properties of putative excitatory and inhibitory neurons in L2/3 of the mouse barrel cortex. Excitatory neurons show high thresholds and strong adaptation, making them fire sparsely and resulting in a strong compression of information, whereas inhibitory neurons that favour fast spiking transfer more information. Next, we turn to computational modelling and ask how two properties influence information transfer: 1) spike-frequency adaptation and 2) the shape of the IV-curve. We find that a subthreshold (but not threshold) adaptation, the 'h-current', and a properly tuned leak conductance can increase the information transfer of a neuron, whereas threshold adaptation can increase its working range. Finally, we verify the effect of the IV-curve slope in our experimental recordings and show that excitatory neurons form a more heterogeneous population than inhibitory neurons. These relationships between intrinsic neural features and neural coding that had not been quantified before will aid computational, theoretical and systems neuroscientists in understanding how neuronal populations can alter their coding properties, such as through the impact of neuromodulators. Why the variability of intrinsic properties of excitatory neurons is larger than that of inhibitory ones is an exciting question, for which future research is needed.

Biophysics-inspired spike rate adaptation for computationally efficient phenomenological nerve modeling.

This study introduces and evaluates the PHAST+ model, part of a computational framework designed to simulate the behavior of auditory nerve fibers in response to the electrical stimulation from a cochlear implant. PHAST+ incorporates a highly efficient method for calculating accommodation and adaptation, making it particularly suited for simulations over extended stimulus durations. The proposed method uses a leaky integrator inspired by classic biophysical nerve models. Through evaluation against single-fiber animal data, our findings demonstrate the model's effectiveness across various stimuli, including short pulse trains with variable amplitudes and rates. Notably, the PHAST+ model performs better than its predecessor, PHAST (a phenomenological model by van Gendt et al.), particularly in simulations of prolonged neural responses. While PHAST+ is optimized primarily on spike rate decay, it shows good behavior on several other neural measures, such as vector strength and degree of adaptation. The future implications of this research are promising. PHAST+ drastically reduces the computational burden to allow the real-time simulation of neural behavior over extended periods, opening the door to future simulations of psychophysical experiments and multi-electrode stimuli for evaluating novel speech-coding strategies for cochlear implants.

Asymmetric pulses delivered by a cochlear implant allow a reduction in evoked firing rate and in spatial activation in the guinea pig auditory cortex.

Despite that fact that the cochlear implant (CI) is one of the most successful neuro-prosthetic devices which allows hearing restoration, several aspects still need to be improved. Interactions between stimulating electrodes through current spread occurring within the cochlea drastically limit the number of discriminable frequency channels and thus can ultimately result in poor speech perception. One potential solution relies on the use of new pulse shapes, such as asymmetric pulses, which can potentially reduce the current spread within the cochlea. The present study characterized the impact of changing electrical pulse shapes from the standard biphasic symmetric to the asymmetrical shape by quantifying the evoked firing rate and the spatial activation in the guinea pig primary auditory cortex (A1). At a fixed charge, the firing rate and the spatial activation in A1 decreased by 15 to 25 % when asymmetric pulses were used to activate the auditory nerve fibers, suggesting a potential reduction of the spread of excitation inside the cochlea. A strong "polarity-order" effect was found as the reduction was more pronounced when the first phase of the pulse was cathodic with high amplitude. These results suggest that the use of asymmetrical pulse shapes in clinical settings can potentially reduce the channel interactions in CI users.