Repeated stress triggers seeking of a starvation-like state in anxiety-prone female mice.

Elevated anxiety often precedes anorexia nervosa and persists after weight restoration. Patients with anorexia nervosa often describe self-starvation as pleasant, potentially because food restriction can be anxiolytic. Here, we tested whether repeated stress can cause animals to prefer a starvation-like state. We developed a virtual reality place preference paradigm in which head-fixed mice can voluntarily seek a starvation-like state induced by optogenetic stimulation of hypothalamic agouti-related peptide (AgRP) neurons. Prior to stress exposure, males but not females showed a mild aversion to AgRP stimulation. Strikingly, following multiple days of stress, a subset of females developed a strong preference for AgRP stimulation that was predicted by high baseline anxiety. Such stress-induced changes in preference were reflected in changes in facial expressions during AgRP stimulation. Our study suggests that stress may cause females predisposed to anxiety to seek a starvation state and provides a powerful experimental framework for investigating the underlying neural mechanisms.

A spatially-resolved transcriptional atlas of the murine dorsal pons at single-cell resolution.

The "dorsal pons", or "dorsal pontine tegmentum" (dPnTg), is part of the brainstem. It is a complex, densely packed region whose nuclei are involved in regulating many vital functions. Notable among them are the parabrachial nucleus, the Kölliker Fuse, the Barrington nucleus, the locus coeruleus, and the dorsal, laterodorsal, and ventral tegmental nuclei. In this study, we applied single-nucleus RNA-seq (snRNA-seq) to resolve neuronal subtypes based on their unique transcriptional profiles and then used multiplexed error robust fluorescence in situ hybridization (MERFISH) to map them spatially. We sampled ~1 million cells across the dPnTg and defined the spatial distribution of over 120 neuronal subtypes. Our analysis identified an unpredicted high transcriptional diversity in this region and pinpointed the unique marker genes of many neuronal subtypes. We also demonstrated that many neuronal subtypes are transcriptionally similar between humans and mice, enhancing this study's translational value. Finally, we developed a freely accessible, GPU and CPU-powered dashboard ( http://harvard.heavy.ai:6273/ ) that combines interactive visual analytics and hardware-accelerated SQL into a data science framework to allow the scientific community to query and gain insights into the data.

Cortical reactivations predict future sensory responses.

Many theories of offline memory consolidation posit that the pattern of neurons activated during a salient sensory experience will be faithfully reactivated, thereby stabilizing the pattern1,2. However, sensory-evoked patterns are not stable but, instead, drift across repeated experiences3-6. Here, to investigate the relationship between reactivations and the drift of sensory representations, we imaged the calcium activity of thousands of excitatory neurons in the mouse lateral visual cortex. During the minute after a visual stimulus, we observed transient, stimulus-specific reactivations, often coupled with hippocampal sharp-wave ripples. Stimulus-specific reactivations were abolished by local cortical silencing during the preceding stimulus. Reactivations early in a session systematically differed from the pattern evoked by the previous stimulus-they were more similar to future stimulus response patterns, thereby predicting both within-day and across-day representational drift. In particular, neurons that participated proportionally more or less in early stimulus reactivations than in stimulus response patterns gradually increased or decreased their future stimulus responses, respectively. Indeed, we could accurately predict future changes in stimulus responses and the separation of responses to distinct stimuli using only the rate and content of reactivations. Thus, reactivations may contribute to a gradual drift and separation in sensory cortical response patterns, thereby enhancing sensory discrimination7.

A spatially-resolved transcriptional atlas of the murine dorsal pons at single-cell resolution.

The "dorsal pons", or "dorsal pontine tegmentum" (dPnTg), is part of the brainstem. It is a complex, densely packed region whose nuclei are involved in regulating many vital functions. Notable among them are the parabrachial nucleus, the Kölliker Fuse, the Barrington nucleus, the locus coeruleus, and the dorsal, laterodorsal, and ventral tegmental nuclei. In this study, we applied single-nucleus RNA-seq (snRNA-seq) to resolve neuronal subtypes based on their unique transcriptional profiles and then used multiplexed error robust fluorescence in situ hybridization (MERFISH) to map them spatially. We sampled ~1 million cells across the dPnTg and defined the spatial distribution of over 120 neuronal subtypes. Our analysis identified an unpredicted high transcriptional diversity in this region and pinpointed many neuronal subtypes' unique marker genes. We also demonstrated that many neuronal subtypes are transcriptionally similar between humans and mice, enhancing this study's translational value. Finally, we developed a freely accessible, GPU and CPU-powered dashboard (http://harvard.heavy.ai:6273/) that combines interactive visual analytics and hardware-accelerated SQL into a data science framework to allow the scientific community to query and gain insights into the data.

A spatial human thymus cell atlas mapped to a continuous tissue axis.

T cells develop from circulating precursors, which enter the thymus and migrate throughout specialised sub-compartments to support maturation and selection. This process starts already in early fetal development and is highly active until the involution of the thymus in adolescence. To map the micro-anatomical underpinnings of this process in pre- vs. post-natal states, we undertook a spatially resolved analysis and established a new quantitative morphological framework for the thymus, the Cortico-Medullary Axis. Using this axis in conjunction with the curation of a multimodal single-cell, spatial transcriptomics and high-resolution multiplex imaging atlas, we show that canonical thymocyte trajectories and thymic epithelial cells are highly organised and fully established by post-conception week 12, pinpoint TEC progenitor states, find that TEC subsets and peripheral tissue genes are associated with Hassall's Corpuscles and uncover divergence in the pace and drivers of medullary entry between CD4 vs. CD8 T cell lineages. These findings are complemented with a holistic toolkit for spatial analysis and annotation, providing a basis for a detailed understanding of T lymphocyte development.

Chronic stress triggers seeking of a starvation-like state in anxiety-prone female mice.

Elevated anxiety often precedes anorexia nervosa and persists after weight restoration. Patients with anorexia nervosa often describe hunger as pleasant, potentially because food restriction can be anxiolytic. Here, we tested whether chronic stress can cause animals to prefer a starvation-like state. We developed a virtual reality place preference paradigm in which head-fixed mice can voluntarily seek a starvation-like state induced by optogenetic stimulation of hypothalamic agouti-related peptide (AgRP) neurons. Prior to stress induction, male but not female mice showed mild aversion to AgRP stimulation. Strikingly, following chronic stress, a subset of females developed a strong preference for AgRP stimulation that was predicted by high baseline anxiety. Such stress-induced changes in preference were reflected in changes in facial expressions during AgRP stimulation. Our study suggests that stress may cause females predisposed to anxiety to seek a starvation state, and provides a powerful experimental framework for investigating the underlying neural mechanisms.

Thalamic control of sensory processing and spindles in a biophysical somatosensory thalamoreticular circuit model of wakefulness and sleep.

Thalamoreticular circuitry plays a key role in arousal, attention, cognition, and sleep spindles, and is linked to several brain disorders. A detailed computational model of mouse somatosensory thalamus and thalamic reticular nucleus has been developed to capture the properties of over 14,000 neurons connected by 6 million synapses. The model recreates the biological connectivity of these neurons, and simulations of the model reproduce multiple experimental findings in different brain states. The model shows that inhibitory rebound produces frequency-selective enhancement of thalamic responses during wakefulness. We find that thalamic interactions are responsible for the characteristic waxing and waning of spindle oscillations. In addition, we find that changes in thalamic excitability control spindle frequency and their incidence. The model is made openly available to provide a new tool for studying the function and dysfunction of the thalamoreticular circuitry in various brain states.

High-throughput morphometric and transcriptomic profiling uncovers composition of naïve and sensory-deprived cortical cholinergic VIP/CHAT neurons.

Cortical neuronal networks control cognitive output, but their composition and modulation remain elusive. Here, we studied the morphological and transcriptional diversity of cortical cholinergic VIP/ChAT interneurons (VChIs), a sparse population with a largely unknown function. We focused on VChIs from the whole barrel cortex and developed a high-throughput automated reconstruction framework, termed PopRec, to characterize hundreds of VChIs from each mouse in an unbiased manner, while preserving 3D cortical coordinates in multiple cleared mouse brains, accumulating thousands of cells. We identified two fundamentally distinct morphological types of VChIs, bipolar and multipolar that differ in their cortical distribution and general morphological features. Following mild unilateral whisker deprivation on postnatal day seven, we found after three weeks both ipsi- and contralateral dendritic arborization differences and modified cortical depth and distribution patterns in the barrel fields alone. To seek the transcriptomic drivers, we developed NuNeX, a method for isolating nuclei from fixed tissues, to explore sorted VChIs. This highlighted differentially expressed neuronal structural transcripts, altered exitatory innervation pathways and established Elmo1 as a key regulator of morphology following deprivation.

A calcium-based plasticity model for predicting long-term potentiation and depression in the neocortex.

Pyramidal cells (PCs) form the backbone of the layered structure of the neocortex, and plasticity of their synapses is thought to underlie learning in the brain. However, such long-term synaptic changes have been experimentally characterized between only a few types of PCs, posing a significant barrier for studying neocortical learning mechanisms. Here we introduce a model of synaptic plasticity based on data-constrained postsynaptic calcium dynamics, and show in a neocortical microcircuit model that a single parameter set is sufficient to unify the available experimental findings on long-term potentiation (LTP) and long-term depression (LTD) of PC connections. In particular, we find that the diverse plasticity outcomes across the different PC types can be explained by cell-type-specific synaptic physiology, cell morphology and innervation patterns, without requiring type-specific plasticity. Generalizing the model to in vivo extracellular calcium concentrations, we predict qualitatively different plasticity dynamics from those observed in vitro. This work provides a first comprehensive null model for LTP/LTD between neocortical PC types in vivo, and an open framework for further developing models of cortical synaptic plasticity.

Statistical Emulation of Neural Simulators: Application to Neocortical L2/3 Large Basket Cells.

Many scientific systems are studied using computer codes that simulate the phenomena of interest. Computer simulation enables scientists to study a broad range of possible conditions, generating large quantities of data at a faster rate than the laboratory. Computer models are widespread in neuroscience, where they are used to mimic brain function at different levels. These models offer a variety of new possibilities for the neuroscientist, but also numerous challenges, such as: where to sample the input space for the simulator, how to make sense of the data that is generated, and how to estimate unknown parameters in the model. Statistical emulation can be a valuable complement to simulator-based research. Emulators are able to mimic the simulator, often with a much smaller computational burden and they are especially valuable for parameter estimation, which may require many simulator evaluations. This work compares different statistical models that address these challenges, and applies them to simulations of neocortical L2/3 large basket cells, created and run with the NEURON simulator in the context of the European Human Brain Project. The novelty of our approach is the use of fast empirical emulators, which have the ability to accelerate the optimization process for the simulator and to identify which inputs (in this case, different membrane ion channels) are most influential in affecting simulated features. These contributions are complementary, as knowledge of the important features can further improve the optimization process. Subsequent research, conducted after the process is completed, will gain efficiency by focusing on these inputs.

Visual association cortex links cues with conjunctions of reward and locomotor contexts.

Postrhinal cortex (POR) and neighboring lateral visual association areas are necessary for identifying objects and interpreting them in specific contexts, but how POR neurons encode the same object across contexts remains unclear. Here, we imaged excitatory neurons in mouse POR across tens of days prior to and throughout initial cue-reward learning and reversal learning. We assessed responses to the same cue when it was rewarded or unrewarded, during both locomotor and stationary contexts. Surprisingly, a large class of POR neurons were minimally cue-driven prior to learning. After learning, distinct clusters within this class responded selectively to a given cue when presented in a specific conjunction of reward and locomotion contexts. In addition, another class contained clusters of neurons whose cue responses were more transient, insensitive to reward learning, and adapted over thousands of presentations. These two classes of POR neurons may support context-dependent interpretation and context-independent identification of sensory cues.

Realistic retinal modeling unravels the differential role of excitation and inhibition to starburst amacrine cells in direction selectivity.

Retinal direction-selectivity originates in starburst amacrine cells (SACs), which display a centrifugal preference, responding with greater depolarization to a stimulus expanding from soma to dendrites than to a collapsing stimulus. Various mechanisms were hypothesized to underlie SAC centrifugal preference, but dissociating them is experimentally challenging and the mechanisms remain debatable. To address this issue, we developed the Retinal Stimulation Modeling Environment (RSME), a multifaceted data-driven retinal model that encompasses detailed neuronal morphology and biophysical properties, retina-tailored connectivity scheme and visual input. Using a genetic algorithm, we demonstrated that spatiotemporally diverse excitatory inputs-sustained in the proximal and transient in the distal processes-are sufficient to generate experimentally validated centrifugal preference in a single SAC. Reversing these input kinetics did not produce any centrifugal-preferring SAC. We then explored the contribution of SAC-SAC inhibitory connections in establishing the centrifugal preference. SAC inhibitory network enhanced the centrifugal preference, but failed to generate it in its absence. Embedding a direction selective ganglion cell (DSGC) in a SAC network showed that the known SAC-DSGC asymmetric connectivity by itself produces direction selectivity. Still, this selectivity is sharpened in a centrifugal-preferring SAC network. Finally, we use RSME to demonstrate the contribution of SAC-SAC inhibitory connections in mediating direction selectivity and recapitulate recent experimental findings. Thus, using RSME, we obtained a mechanistic understanding of SACs' centrifugal preference and its contribution to direction selectivity.

The Role of Hub Neurons in Modulating Cortical Dynamics.

Many neurodegenerative diseases are associated with the death of specific neuron types in particular brain regions. What makes the death of specific neuron types particularly harmful for the integrity and dynamics of the respective network is not well understood. To start addressing this question we used the most up-to-date biologically realistic dense neocortical microcircuit (NMC) of the rodent, which has reconstructed a volume of 0.3 mm3 and containing 31,000 neurons, ∼37 million synapses, and 55 morphological cell types arranged in six cortical layers. Using modern network science tools, we identified hub neurons in the NMC, that are connected synaptically to a large number of their neighbors and systematically examined the impact of abolishing these cells. In general, the structural integrity of the network is robust to cells' attack; yet, attacking hub neurons strongly impacted the small-world topology of the network, whereas similar attacks on random neurons have a negligible effect. Such hub-specific attacks are also impactful on the network dynamics, both when the network is at its spontaneous synchronous state and when it was presented with synchronized thalamo-cortical visual-like input. We found that attacking layer 5 hub neurons is most harmful to the structural and functional integrity of the NMC. The significance of our results for understanding the role of specific neuron types and cortical layers for disease manifestation is discussed.

Using subthreshold events to characterize the functional architecture of the electrically coupled inferior olive network.

The electrical connectivity in the inferior olive (IO) nucleus plays an important role in generating well-timed spiking activity. Here we combined electrophysiological and computational approaches to assess the functional organization of the IO nucleus in mice. Spontaneous fast and slow subthreshold events were commonly encountered during in vitro recordings. We show that whereas the fast events represent intrinsic regenerative activity, the slow events reflect the electrical connectivity between neurons ('spikelets'). Recordings from cell pairs revealed the synchronized occurrence of distinct groups of spikelets; their rate and distribution enabled an accurate estimation of the number of connected cells and is suggestive of a clustered organization. This study thus provides a new perspective on the functional and structural organization of the olivary nucleus and a novel experimental and theoretical approach to study electrically coupled networks.

An efficient analytical reduction of detailed nonlinear neuron models.

Detailed conductance-based nonlinear neuron models consisting of thousands of synapses are key for understanding of the computational properties of single neurons and large neuronal networks, and for interpreting experimental results. Simulations of these models are computationally expensive, considerably curtailing their utility. Neuron_Reduce is a new analytical approach to reduce the morphological complexity and computational time of nonlinear neuron models. Synapses and active membrane channels are mapped to the reduced model preserving their transfer impedance to the soma; synapses with identical transfer impedance are merged into one NEURON process still retaining their individual activation times. Neuron_Reduce accelerates the simulations by 40-250 folds for a variety of cell types and realistic number (10,000-100,000) of synapses while closely replicating voltage dynamics and specific dendritic computations. The reduced neuron-models will enable realistic simulations of neural networks at unprecedented scale, including networks emerging from micro-connectomics efforts and biologically-inspired "deep networks". Neuron_Reduce is publicly available and is straightforward to implement.

Intensify3D: Normalizing signal intensity in large heterogenic image stacks.

Three-dimensional structures in biological systems are routinely evaluated using large image stacks acquired from fluorescence microscopy; however, analysis of such data is muddled by variability in the signal across and between samples. Here, we present Intensify3D: a user-guided normalization algorithm tailored for overcoming common heterogeneities in large image stacks. We demonstrate the use of Intensify3D for analyzing cholinergic interneurons of adult murine brains in 2-Photon and Light-Sheet fluorescence microscopy, as well as of mammary gland and heart tissues. Beyond enhancement in 3D visualization in all samples tested, in 2-Photon in vivo images, this tool corrected errors in feature extraction of cortical interneurons; and in Light-Sheet microscopy, it enabled identification of individual cortical barrel fields and quantification of somata in cleared adult brains. Furthermore, Intensify3D enhanced the ability to separate signal from noise. Overall, the universal applicability of our method can facilitate detection and quantification of 3D structures and may add value to a wide range of imaging experiments.

From Neuron Biophysics to Orientation Selectivity in Electrically Coupled Networks of Neocortical L2/3 Large Basket Cells.

In the neocortex, inhibitory interneurons of the same subtype are electrically coupled with each other via dendritic gap junctions (GJs). The impact of multiple GJs on the biophysical properties of interneurons and thus on their input processing is unclear. The present experimentally based theoretical study examined GJs in L2/3 large basket cells (L2/3 LBCs) with 3 goals in mind: (1) To evaluate the errors due to GJs in estimating the cable properties of individual L2/3 LBCs and suggest ways to correct these errors when modeling these cells and the networks they form; (2) to bracket the GJ conductance value (0.05-0.25 nS) and membrane resistivity (10 000-40 000 Ω cm(2)) of L2/3 LBCs; these estimates are tightly constrained by in vitro input resistance (131 ± 18.5 MΩ) and the coupling coefficient (1-3.5%) of these cells; and (3) to explore the functional implications of GJs, and show that GJs: (i) dynamically modulate the effective time window for synaptic integration; (ii) improve the axon's capability to encode rapid changes in synaptic inputs; and (iii) reduce the orientation selectivity, linearity index, and phase difference of L2/3 LBCs. Our study provides new insights into the role of GJs and calls for caution when using in vitro measurements for modeling electrically coupled neuronal networks.