Sensory neurons of the vagus nerve monitor distention and stretch in the gastrointestinal tract. We used genetically guided anatomical tracing, optogenetics and electrophysiology to identify and characterize two vagal sensory neuronal subtypes expressing Prox2 and Runx3. We show that these neuronal subtypes innervate the esophagus where they display regionalized innervation patterns. Electrophysiological analyses showed that they are both low threshold mechanoreceptors but possess different adaptation properties. Lastly, genetic ablation of Prox2 and Runx3 neurons demonstrated their essential roles for esophageal peristalsis and swallowing in freely behaving animals. Our work reveals the identity and function of the vagal neurons that provide mechanosensory feedback from the esophagus to the brain and could lead to better understanding and treatment of esophageal motility disorders.

Recently, the field of computational neuroscience has seen an explosion of the use of trained recurrent network models (RNNs) to model patterns of neural activity. These RNN models are typically characterized by tuned recurrent interactions between rate 'units' whose dynamics are governed by smooth, continuous differential equations. However, the response of biological single neurons is better described by all-or-none events - spikes - that are triggered in response to the processing of their synaptic input by the complex dynamics of their membrane. One line of research has attempted to resolve this discrepancy by linking the average firing probability of a population of simplified spiking neuron models to rate dynamics similar to those used for RNN units. However, challenges remain to account for complex temporal dependencies in the biological single neuron response and for the heterogeneity of synaptic input across the population. Here, we make progress by showing how to derive dynamic rate equations for a population of spiking neurons with multi-timescale adaptation properties - as this was shown to accurately model the response of biological neurons - while they receive independent time-varying inputs, leading to plausible asynchronous activity in the network. The resulting rate equations yield an insightful segregation of the population's response into dynamics that are driven by the mean signal received by the neural population, and dynamics driven by the variance of the input across neurons, with respective timescales that are in agreement with slice experiments. Further, these equations explain how input variability can shape log-normal instantaneous rate distributions across neurons, as observed in vivo. Our results help interpret properties of the neural population response and open the way to investigating whether the more biologically plausible and dynamically complex rate model we derive could provide useful inductive biases if used in an RNN to solve specific tasks.