We present a method, open-source software, and experiments which embed arbitrary deformation vector fields produced by any method (e.g., ANTs or VoxelMorph) in the Large Deformation Diffeomorphic Metric Mapping (LDDMM) framework. This decouples formal diffeomorphic shape analysis from image registration, which has many practical benefits. Shape analysis can be added to study designs without modification to already chosen image registration methods and existing databases of deformation fields can be reanalyzed within the LDDMM framework without repeating image registrations. Pairwise time series studies can be extended to full time series regression with minimal added computing. The diffeomorphic rigor of image registration methods can be compared by embedding deformation fields and comparing projection distances. Finally, the added value of formal diffeomorphic shape analysis can be more fairly evaluated when it is derived from and compared to a baseline set of deformation fields. In brief, the method is a straightforward use of geodesic shooting in diffeomorphisms with a deformation field as the target, rather than an image. This is simpler than the image registration case which leads to a faster implementation that requires fewer user derived parameters.

Messenger RNA (mRNA) transfection enables rapid, transient protein expression without nuclear entry, providing a powerful alternative to DNA or viral delivery in post-mitotic and otherwise difficult-to-transfect cells. Although in vitro transcribed (IVT) mRNAs have revolutionized therapeutic applications, their adoption in experimental biology remains limited by challenges in synthesis, variability across cell types, and concerns about cytotoxicity. Here, we define design principles that maximize IVT mRNA performance across diverse cellular and organismal systems. Through systematic comparison of capping strategies and base modifications, including N1-methyl-pseudouridine, 5-methylcytidine, and 5-methoxyuridine, we identify modifications that enhance translation while minimizing activation of cellular stress responses. Optimized transcripts drive robust protein expression within four hours, persist for up to one week, and support multiplexed expression of structurally and functionally distinct proteins in mammalian cells, including cancer cell lines, iPSC-derived systems, primary cells, and organoids, as well as in vivo in zebrafish embryos and in less genetically tractable models such as Danionella cerebrum and sea urchin embryos. To further expand accessibility for community use, we developed mRNAbow, a platform for generating low-toxicity mRNAs encoding organelle-targeted fluorescent proteins and biosensors for multiplex imaging, with corresponding plasmids made publicly available. Together, these advances establish a generalizable framework for IVT mRNA design and expand experimental access to synthetic mRNA technologies for dissecting cellular architecture and dynamics.

Understanding how neural signals control muscle activity during behavior is a key challenge in motor neuroscience. To this end, recent advances in intramuscular multielectrode arrays have enabled high-quality multichannel recordings of many motor unit action potentials (MUAPs) in freely moving subjects. However, identifying individual MUAP events within multichannel recordings is a significant challenge for existing spike sorting methods, which are typically optimized for identifying action potentials from neurons in the brain. To overcome this challenge, we developed the Enhanced Motor Unit sorter (EMUsort), an extension of Kilosort4 (KS4) that achieves high-performance MUAP spike sorting. We applied EMUsort to high-resolution intramuscular recordings from rat forelimb during locomotion and monkey forelimb during a reaching task. EMUsort improves upon prior methods by addressing key challenges encountered with MUAP datasets, including: 1) long time delays across electrodes due to propagation along muscle fibers, 2) more complex waveform shapes compared to neuronal action potentials, and 3) a high degree of MUAP overlap due to cumulative motor unit recruitment. We compared EMUsort to existing spike sorting methods quantitatively using simulated datasets that closely emulated the rat and monkey datasets we recorded. EMUsort provided median error rate reductions of 67.5% and 49.9% during periods of high motor unit activation for the rat and monkey datasets, respectively. In sum, EMUsort provides a substantial improvement to MUAP spike sorter accuracy, especially during regions of high MUAP overlap, in an easy-to-use software package.

Single-walled carbon nanotubes (SWCNTs) functionalized with single-stranded DNAs can function as near-infrared nanosensors for molecular analytes. However, predicting which analytes elicit strong optical responses for specific nanosensors remains challenging. We developed machine learning (ML) models to predict analyte-induced fluorescence changes in a DNA–SWCNT dopamine nanosensor. Using a data set of 63 small molecules sampling chemical space around dopamine, we encoded analytes with RDKit fingerprints, with or without HOMO and LUMO energies, and applied principal component analysis to identify structural motifs associated with optical response strength. We trained support vector regression and classification models using two strategies: ensembles of 200 models and cross-validation. Regression models achieved mean R2 values of 0.2–0.4, with cross-validation outperforming ensembles, while classifiers reached mean F1 scores of ∼0.8. Cross-validation performed best for predictions on a blind set of 21 molecules. These findings show that ML can capture structure–response patterns in modest data sets and guide in silico DNA–SWCNT nanosensor design.

Dendrites integrate synaptic inputs to trigger action potentials, and dendrites carry back-propagating action potentials (bAPs) to synapses where these signals contribute to plasticity. Despite strong evidence for a rich repertoire of nonlinear dendritic excitations, the in vivo roles of these excitations in dendritic integration and back-propagation remain uncertain. Here, we used high-speed voltage imaging through a chronically implanted microprism to map membrane potential dynamics from basal to apical dendrites of CA1 neurons in mice navigating in a virtual reality environment. Despite complex dendritic branch morphology, the dynamics were largely captured by 2 or 3 electrical compartments: basal, soma, and apical. Fast dendritic spikes almost always started from bAPs, indicating that dendritic spikes are primarily a consequence rather than a cause of somatic spiking. These fast spikes sometimes triggered slower apical dendritic plateau depolarizations, which drove complex spikes at the soma. We found that the biophysics of dendritic excitability determined the distribution of simple and complex spikes across a place field. Our results show how CA1 pyramidal neurons convert synaptic inputs to spiking outputs and suggest a primary role of dendritic nonlinearities in mediating activity-dependent plasticity.

Online monitoring and quantification of neural signals has tremendous value both for neurofeedback experiments and for brain-computer interfaces. Unfortunately, established methods of online monitoring primarily involve the use of thresholded neural activity rather than sorted single-neuron spikes. The recent introduction of large-scale, high-density electrophysiology has enabled the recording of activity from hundreds of neurons simultaneously in both model organisms and human participants. This development highlights the need for a robust and easily implementable system for sorting spikes during data collection for ‘live’ analyses of neuronal signals. Here, we describe a system for live sorting of neuronal activity (LSS) based on the widely used Kilosort platform. The LSS workflow utilizes an initial period of recorded neural data to identify waveform templates using Kilosort 4. LSS then interfaces with the SpikeGLX API to retrieve small batches (e.g. 50 ms) of data and for processing online. We measured the similarity of single-neuron activity sorted live by LSS to that sorted offline in neurophysiological recordings from macaque visual cortex using Neuropixels probes. We show that LSS closely replicates the post-stimulus time histograms and visual response tuning curves of single-neurons obtained using offline sorting. Furthermore, we show that decoding neural signals online with LSS consistently outperforms online decoding of thresholded activity, and that LSS can achieve the same performance as that obtained with offline sorting.

Neural representations of information are shaped by long-range input and local network interactions. Previous studies linking neural coding and cortical connectivity have focused on input-driven activity in the sensory cortex. Here we studied neural activity in the motor cortex while mice gathered rewards with multidirectional tongue reaching. This behaviour does not require training, allowing us to probe neural coding and connectivity before activity is shaped by extended learning. Motor cortex neurons were tuned to target location and reward outcome, and typically responded during and after movements. We studied the underlying network interactions in vivo by estimating causal neural connections using an all-optical method. Mapping connectivity between more than 20,000,000 excitatory neuron pairs showed a multi-scale columnar architecture in layer 2/3 of the motor cortex. Neurons displayed local (less than 100 µm) like-to-like excitatory connectivity according to target-location tuning, and inhibition over longer spatial scales. Connectivity patterns comprised a continuum, with abundant sparsely connected neurons and rare densely connected neurons that function as network hubs. Hub neurons were weakly tuned to target location and reward outcome but influenced more neighbouring neurons. This network of neurons, encoding location and outcome of movements to different motor goals, may be a general substrate for rapid learning of complex, goal-directed behaviours.

Cells exhibit a mysterious form of selective heritable short-term memory, influencing outcomes as diverse as cell fate decisions in embryos and environmental responses in cancer cells and bacteria. Here, we present a simple theoretical framework explaining how this selective memory can arise from the reactions regulating molecular levels in cells. Our key insight is that related cells retain more similar molecular concentrations relative to random cells when a greater variance of possible concentration states is created during a single cell generation than is created by cell division across a population. This persistence of molecular similarity down a lineage constitutes a form of heritable short-term memory. We identify the biochemical networks that produce, modify, and degrade molecules as an underexplored source of these additional molecular concentration states. Using experimentally informed simulations, we find that the strength and duration of molecular similarity down a lineage depend on tunable network properties, explaining why some cellular traits persist only briefly while others last generations. These contributions to molecular concentration variance from biochemical reaction networks act in concert with gene expression and other regulatory processes to shape the protein composition of cells. Our framework yields clear, testable predictions for determining how biochemical network architectures drive non-genetic cellular inheritance.

Visual systems across species transform photoreceptor inputs into diverse perceptual representations through hierarchical networks that extract features via parallel pathways. In Drosophila, the optic lobes are layered, retinotopic visual processing centers that contain two-thirds of the brain’s neurons and support diverse visually guided behaviors. Although this architecture has long suggested hierarchical and parallel organization, a system-wide account of how behaviorally relevant visual features are routed and integrated across a complete visual system—in any animal—has remained elusive. The new male fly connectome now provides the synapse-level wiring needed to trace visual information from photoreceptors through the optic lobes and across the central brain. Applying a network-based analysis of information flow, we reveal a multi-layered architecture organized into distinct, functionally interpretable pathways. Using this framework to propagate signals through these pathways predicts receptive-field structure and feature selectivity consistent with physiological data, enabling large-scale functional annotation of thousands of neuron types. We find that distinct visual input channels are broadly distributed throughout the brain, yet converge in focal regions of feature specificity and acute spatial vision. Together, these analyses provide a neuron-level, connectome-based view of how a brain organizes and transforms visual input.

The early formation of sensorimotor circuits is essential for survival. While the development and function of exteroceptive circuits and their associated motor pathways are well characterized, far less is known about the circuits that convey viscerosensory inputs to the brain and transmit visceromotor commands from the central nervous system to internal organs. Technical limitations, such as the in utero development of viscerosensory and visceromotor circuits and the invasiveness of procedures required to access them, have hindered studies of their functional development in mammals. Using larval zebrafish—which are genetically accessible and optically transparent—we tracked, in vivo, how cardiosensory and cardiomotor neural circuits assemble and begin to function. We uncovered a staged program. First, a minimal efferent circuit suffices for heart-rate control: direct brain-to-heart vagal motor innervation is required, intracardiac neurons are not, and heart rate is governed exclusively by the motor vagus nerve. Within the hindbrain, we functionally localize a vagal premotor population that drives this early efferent control. Second, sympathetic innervation arrives and enhances the dynamics and amplitude of cardiac responses, as neurons in the most anterior sympathetic ganglia acquire the ability to drive cardiac acceleration. These neurons exhibit proportional, integral, and derivative–like relationships to heart rate, consistent with controller motifs that shape gain and dynamics. Third, vagal sensory neurons innervate the heart. Distinct subsets increase activity when heart rate falls or rises, and across spontaneous fluctuations, responses to aversive stimuli, and optogenetically evoked cardiac perturbations, their dynamics are captured by a single canonical temporal kernel with neuron-specific phase offsets, supporting a population code for heart rate. This temporally segregated maturation isolates three experimentally tractable regimes—unidirectional brain-to-heart communication, dual efferent control, and closed-loop control after sensory feedback engages—providing a framework for mechanistic dissection of organism-wide heart–brain circuits.