Equilibrium formally can be represented as an ensemble of uncoupled systems undergoing unbiased dynamics in which detailed balance is maintained. Many nonequilibrium processes can be described by suitable subsets of the equilibrium ensemble. Here, we employ the “weighted ensemble” (WE) simulation protocol [Huber and Kim, Biophys. J.1996, 70, 97–110] to generate equilibrium trajectory ensembles and extract nonequilibrium subsets for computing kinetic quantities. States do not need to be chosen in advance. The procedure formally allows estimation of kinetic rates between arbitrary states chosen after the simulation, along with their equilibrium populations. We also describe a related history-dependent matrix procedure for estimating equilibrium and nonequilibrium observables when phase space has been divided into arbitrary non-Markovian regions, whether in WE or ordinary simulation. In this proof-of-principle study, these methods are successfully applied and validated on two molecular systems: explicitly solvated methane association and the implicitly solvated Ala4 peptide. We comment on challenges remaining in WE calculations.

Spike sorting is the computational process of extracting the firing times of single neurons from recordings of local electrical fields. This is an important but hard problem in neuroscience, complicated by the non-stationarity of the recordings and the dense overlap in electrical fields between nearby neurons. To solve the spike sorting problem, we have continuously developed over the past eight years a framework known as Kilosort. This paper describes the various algorithmic steps introduced in different versions of Kilosort. We also report the development of Kilosort4, a new version with substantially improved performance due to new clustering algorithms inspired by graph-based approaches. To test the performance of Kilosort, we developed a realistic simulation framework which uses densely sampled electrical fields from real experiments to generate non-stationary spike waveforms and realistic noise. We find that nearly all versions of Kilosort outperform other algorithms on a variety of simulated conditions, and Kilosort4 performs best in all cases, correctly identifying even neurons with low amplitudes and small spatial extents in high drift conditions.

The human cerebral cortex, pivotal for advanced cognitive functions, is composed of six distinct layers and dozens of functionally specialized areas. The layers and areas are distinguished both molecularly, by diverse neuronal and glial cell subtypes, and structurally, through intricate spatial organization3,4. While single-cell transcriptomics studies have advanced molecular characterization of human cortical development, a critical gap exists due to the loss of spatial context during cell dissociation. Here, we utilized multiplexed error-robust fluorescence in situ hybridization (MERFISH)9, augmented with deep-learning-based cell segmentation, to examine the molecular, cellular, and cytoarchitectural development of human fetal cortex with spatially resolved single-cell resolution. Our extensive spatial atlas, encompassing 16 million single cells, spans eight cortical areas across four time points in the second and third trimesters. We uncovered an early establishment of the six-layer structure, identifiable in the laminar distribution of excitatory neuronal subtypes by mid-gestation, long before the emergence of cytoarchitectural layers. Notably, while anterior-posterior gradients of neuronal subtypes were generally observed in most cortical areas, a striking exception was the sharp molecular border between primary (V1) and secondary visual cortices (V2) at gestational week 20. Here we discovered an abrupt binary shift in neuronal subtype specification at the earliest stages, challenging the notion that continuous morphogen gradients dictate mid-gestation cortical arealization. Moreover, integrating single-nuclei RNA-sequencing and in situ whole transcriptomics revealed an early upregulation of synaptogenesis in V1-specific Layer 4 neurons, suggesting a role of synaptogenesis in this discrete border formation. Collectively, our findings underscore the crucial role of spatial relationships in determining the molecular specification of cortical layers and areas. This work not only provides a valuable resource for the field, but also establishes a spatially resolved single-cell analysis paradigm that paves the way for a comprehensive developmental atlas of the human brain.

Spike sorting is the computational process of extracting the firing times of single neurons from recordings of local electrical fields. This is an important but hard problem in neuroscience, made complicated by the nonstationarity of the recordings and the dense overlap in electrical fields between nearby neurons. To address the spike-sorting problem, we have been openly developing the Kilosort framework. Here we describe the various algorithmic steps introduced in different versions of Kilosort. We also report the development of Kilosort4, a version with substantially improved performance due to clustering algorithms inspired by graph-based approaches. To test the performance of Kilosort, we developed a realistic simulation framework that uses densely sampled electrical fields from real experiments to generate nonstationary spike waveforms and realistic noise. We found that nearly all versions of Kilosort outperformed other algorithms on a variety of simulated conditions and that Kilosort4 performed best in all cases, correctly identifying even neurons with low amplitudes and small spatial extents in high drift conditions.

Sensory cortices are active in the absence of external sensory stimuli. To understand the nature of this ongoing activity, we used two-photon calcium imaging to record from over 10,000 neurons in the visual cortex of mice awake in darkness while monitoring their behavior videographically. Ongoing population activity was multidimensional, exhibiting at least 100 significant dimensions, some of which were related to the spontaneous behaviors of the mice. The largest single dimension was correlated with the running speed and pupil area, while a 16-dimensional summary of orofacial behaviors could predict ~45% of the explainable neural variance. Electrophysiological recordings with 8 simultaneous Neuropixels probes revealed a similar encoding of high-dimensional orofacial behaviors across multiple forebrain regions. Representation of motor variables continued uninterrupted during visual stimulus presentation, occupying dimensions nearly orthogonal to the stimulus responses. Our results show that a multidimensional representation of motor state is encoded across the forebrain, and is integrated with visual input by neuronal populations in primary visual cortex.

A surprising finding of recent studies in mouse is the dominance of widespread movement-related activity throughout the brain, including in early sensory areas. In awake subjects, failing to account for movement risks misattributing movement-related activity to other (e.g., sensory or cognitive) processes. In this article, we 1) review task designs for separating task-related and movement-related activity, 2) review three 'case studies' in which not considering movement would have resulted in critically different interpretations of neuronal function, and 3) discuss functional couplings that may prevent us from ever fully isolating sensory, motor, and cognitive-related activity. Our main thesis is that neural signals related to movement are ubiquitous, and therefore ought to be considered first and foremost when attempting to correlate neuronal activity with task-related processes.

In a recent publication, Ma et al [1] claim that a transformer-based cellular segmentation method called Mediar [2] — which won a Neurips challenge — outperforms Cellpose [3] (0.897 vs 0.543 median F1 score). Here we show that this result was obtained by artificially impairing Cellpose in multiple ways. When we removed these impairments, Cellpose outperformed Mediar (0.861 vs 0.826 median F1 score on the updated test set). To further investigate the performance of transformers for cellular segmentation, we replaced the Cellpose backbone with a transformer. The transformer-Cellpose model also did not outperform the standard Cellpose (0.848 median F1 test score). Our results suggest that transformers do not advance the state-of-the-art in cellular segmentation.

Visceral sensory pathways mediate homeostatic reflexes, the dysfunction of which leads to many neurological disorders. The Bezold-Jarisch reflex (BJR), first described in 1867, is a cardioinhibitory reflex that is speculated to be mediated by vagal sensory neurons (VSNs) that also triggers syncope. However, the molecular identity, anatomical organization, physiological characteristics and behavioural influence of cardiac VSNs remain mostly unknown. Here we leveraged single-cell RNA-sequencing data and HYBRiD tissue clearing to show that VSNs that express neuropeptide Y receptor Y2 (NPY2R) predominately connect the heart ventricular wall to the area postrema. Optogenetic activation of NPY2R VSNs elicits the classic triad of BJR responses-hypotension, bradycardia and suppressed respiration-and causes an animal to faint. Photostimulation during high-resolution echocardiography and laser Doppler flowmetry with behavioural observation revealed a range of phenotypes reflected in clinical syncope, including reduced cardiac output, cerebral hypoperfusion, pupil dilation and eye-roll. Large-scale Neuropixels brain recordings and machine-learning-based modelling showed that this manipulation causes the suppression of activity across a large distributed neuronal population that is not explained by changes in spontaneous behavioural movements. Additionally, bidirectional manipulation of the periventricular zone had a push-pull effect, with inhibition leading to longer syncope periods and activation inducing arousal. Finally, ablating NPY2R VSNs specifically abolished the BJR. Combined, these results demonstrate a genetically defined cardiac reflex that recapitulates characteristics of human syncope at physiological, behavioural and neural network levels.