Sustained deep-tissue voltage recording using a fast indicator evolved for two-photon microscopy.

Genetically encoded voltage indicators are emerging tools for monitoring voltage dynamics with cell-type specificity. However, current indicators enable a narrow range of applications due to poor performance under two-photon microscopy, a method of choice for deep-tissue recording. To improve indicators, we developed a multiparameter high-throughput platform to optimize voltage indicators for two-photon microscopy. Using this system, we identified JEDI-2P, an indicator that is faster, brighter, and more sensitive and photostable than its predecessors. We demonstrate that JEDI-2P can report light-evoked responses in axonal termini of Drosophila interneurons and the dendrites and somata of amacrine cells of isolated mouse retina. JEDI-2P can also optically record the voltage dynamics of individual cortical neurons in awake behaving mice for more than 30 min using both resonant-scanning and ULoVE random-access microscopy. Finally, ULoVE recording of JEDI-2P can robustly detect spikes at depths exceeding 400 µm and report voltage correlations in pairs of neurons.

Reconstruction of neocortex: Organelles, compartments, cells, circuits, and activity.

We assembled a semi-automated reconstruction of L2/3 mouse primary visual cortex from ∼250 × 140 × 90 µm3 of electron microscopic images, including pyramidal and non-pyramidal neurons, astrocytes, microglia, oligodendrocytes and precursors, pericytes, vasculature, nuclei, mitochondria, and synapses. Visual responses of a subset of pyramidal cells are included. The data are publicly available, along with tools for programmatic and three-dimensional interactive access. Brief vignettes illustrate the breadth of potential applications relating structure to function in cortical circuits and neuronal cell biology. Mitochondria and synapse organization are characterized as a function of path length from the soma. Pyramidal connectivity motif frequencies are predicted accurately using a configuration model of random graphs. Pyramidal cells receiving more connections from nearby cells exhibit stronger and more reliable visual responses. Sample code shows data access and analysis.

State-dependent pupil dilation rapidly shifts visual feature selectivity.

To increase computational flexibility, the processing of sensory inputs changes with behavioural context. In the visual system, active behavioural states characterized by motor activity and pupil dilation1,2 enhance sensory responses, but typically leave the preferred stimuli of neurons unchanged2-9. Here we find that behavioural state also modulates stimulus selectivity in the mouse visual cortex in the context of coloured natural scenes. Using population imaging in behaving mice, pharmacology and deep neural network modelling, we identified a rapid shift in colour selectivity towards ultraviolet stimuli during an active behavioural state. This was exclusively caused by state-dependent pupil dilation, which resulted in a dynamic switch from rod to cone photoreceptors, thereby extending their role beyond night and day vision. The change in tuning facilitated the decoding of ethological stimuli, such as aerial predators against the twilight sky10. For decades, studies in neuroscience and cognitive science have used pupil dilation as an indirect measure of brain state. Our data suggest that, in addition, state-dependent pupil dilation itself tunes visual representations to behavioural demands by differentially recruiting rods and cones on fast timescales.

Dense and Persistent Odor Representations in the Olfactory Bulb of Awake Mice.

Recording and analysis of neural activity are often biased toward detecting sparse subsets of highly active neurons, masking important signals carried in low-magnitude and variable responses. To investigate the contribution of seemingly noisy activity to odor encoding, we used mesoscale calcium imaging from mice of both sexes to record odor responses from the dorsal surface of bilateral olfactory bulbs (OBs). The outer layer of the mouse OB is comprised of dendrites organized into discrete "glomeruli," which are defined by odor receptor-specific sensory neuron input. We extracted activity from a large population of glomeruli and used logistic regression to classify odors from individual trials with high accuracy. We then used add-in and dropout analyses to determine subsets of glomeruli necessary and sufficient for odor classification. Classifiers successfully predicted odor identity even after excluding sparse, highly active glomeruli, indicating that odor information is redundantly represented across a large population of glomeruli. Additionally, we found that random forest (RF) feature selection informed by Gini inequality (RF Gini impurity, RFGI) reliably ranked glomeruli by their contribution to overall odor classification. RFGI provided a measure of "feature importance" for each glomerulus that correlated with intuitive features like response magnitude. Finally, in agreement with previous work, we found that odor information persists in glomerular activity after the odor offset. Together, our findings support a model of OB odor coding where sparse activity is sufficient for odor identification, but information is widely, redundantly available across a large population of glomeruli, with each glomerulus representing information about more than one odor.

Endobronchial Ultrasound for the Screening of Pulmonary Embolism in Patients with Suspected Lung Cancer: A Prospective Cohort Study.

BACKGROUND: Patients with lung cancer exhibit increased risk of pulmonary embolism (PE). While the contrast phase of computed tomography of the chest in the diagnostic work-up of suspected chest malignancy does not allow reliable detection of PE, it may be feasible to screen for present PE during endobronchial ultrasound (EBUS) examination. OBJECTIVES: The aim of this study was to establish if screening during EBUS for PE in patients with suspected lung cancer is feasible and if positive findings are predictive of PE. METHODS: Patients undergoing EBUS due to suspicion of malignancy of the chest were prospectively enrolled. The pulmonary arteries were assessed during EBUS using a standardized protocol. Patients in whom PE suspicion was raised were referred to confirmatory imaging. RESULTS: From December 2020 to August 2021, 100 patients were included. Median time for vascular assessment during EBUS was 2 min (Q1-Q3: 1-3 min). EBUS identified two suspected PEs (2%), and the number needed to scan was 50. The positive predictive value of EBUS for PE was 100%. CONCLUSION: EBUS for PE screening seems feasible and with limited time use. The PPV of positive findings for the diagnosis of PE is high, but the utility is somewhat limited by a high number needed to scan even in a high-risk population. Based on our findings, we believe that EBUS assessment of the pulmonary vasculature may have a role as a routine screening tool for PE. The assessment for PE should be implemented in EBUS training programmes, as operators should be able to recognize incidental PEs.

In vivo three-photon imaging of activity of GCaMP6-labeled neurons deep in intact mouse brain.

High-resolution optical imaging is critical to understanding brain function. We demonstrate that three-photon microscopy at 1,300-nm excitation enables functional imaging of GCaMP6s-labeled neurons beyond the depth limit of two-photon microscopy. We record spontaneous activity from up to 150 neurons in the hippocampal stratum pyramidale at ∼1-mm depth within an intact mouse brain. Our method creates opportunities for noninvasive recording of neuronal activity with high spatial and temporal resolution deep within scattering brain tissues.

Follicle-stimulating hormone and luteinizing hormone increase Ca2+ in the granulosa cells of mouse ovarian follicles†.

In mammalian ovarian follicles, follicle stimulating hormone (FSH) and luteinizing hormone (LH) signal primarily through the G-protein Gs to elevate cAMP, but both of these hormones can also elevate Ca2+ under some conditions. Here, we investigate FSH- and LH-induced Ca2+ signaling in intact follicles of mice expressing genetically encoded Ca2+ sensors, Twitch-2B and GCaMP6s. At a physiological concentration (1 nM), FSH elevates Ca2+ within the granulosa cells of preantral and antral follicles. The Ca2+ rise begins several minutes after FSH application, peaks at ∼10 min, remains above baseline for another ∼10 min, and depends on extracellular Ca2+. However, suppression of the FSH-induced Ca2+ increase by reducing extracellular Ca2+ does not inhibit FSH-induced phosphorylation of MAP kinase, estradiol production, or the acquisition of LH responsiveness. Like FSH, LH also increases Ca2+, when applied to preovulatory follicles. At a physiological concentration (10 nM), LH elicits Ca2+ oscillations in a subset of cells in the outer mural granulosa layer. These oscillations continue for at least 6 h and depend on the activity of Gq family G-proteins. Suppression of the oscillations by Gq inhibition does not inhibit meiotic resumption, but does delay the time to 50% ovulation by about 3 h. In summary, both FSH and LH increase Ca2+ in the granulosa cells of intact follicles, but the functions of these Ca2+ rises are only starting to be identified.

Community-based benchmarking improves spike rate inference from two-photon calcium imaging data.

In recent years, two-photon calcium imaging has become a standard tool to probe the function of neural circuits and to study computations in neuronal populations. However, the acquired signal is only an indirect measurement of neural activity due to the comparatively slow dynamics of fluorescent calcium indicators. Different algorithms for estimating spike rates from noisy calcium measurements have been proposed in the past, but it is an open question how far performance can be improved. Here, we report the results of the spikefinder challenge, launched to catalyze the development of new spike rate inference algorithms through crowd-sourcing. We present ten of the submitted algorithms which show improved performance compared to previously evaluated methods. Interestingly, the top-performing algorithms are based on a wide range of principles from deep neural networks to generative models, yet provide highly correlated estimates of the neural activity. The competition shows that benchmark challenges can drive algorithmic developments in neuroscience.

Optogenetic manipulation of medullary neurons in the locust optic lobe.

The locust is a widely used animal model for studying sensory processing and its relation to behavior. Due to the lack of genomic information, genetic tools to manipulate neural circuits in locusts are not yet available. We examined whether Semliki Forest virus is suitable to mediate exogenous gene expression in neurons of the locust optic lobe. We subcloned a channelrhodopsin variant and the yellow fluorescent protein Venus into a Semliki Forest virus vector and injected the virus into the optic lobe of locusts ( Schistocerca americana). Fluorescence was observed in all injected optic lobes. Most neurons that expressed the recombinant proteins were located in the first two neuropils of the optic lobe, the lamina and medulla. Extracellular recordings demonstrated that laser illumination increased the firing rate of medullary neurons expressing channelrhodopsin. The optogenetic activation of the medullary neurons also triggered excitatory postsynaptic potentials and firing of a postsynaptic, looming-sensitive neuron, the lobula giant movement detector. These results indicate that Semliki Forest virus is efficient at mediating transient exogenous gene expression and provides a tool to manipulate neural circuits in the locust nervous system and likely other insects. NEW & NOTEWORTHY Using Semliki Forest virus, we efficiently delivered channelrhodopsin into neurons of the locust optic lobe. We demonstrate that laser illumination increases the firing of the medullary neurons expressing channelrhodopsin and elicits excitatory postsynaptic potentials and spiking in an identified postsynaptic target neuron, the lobula giant movement detector neuron. This technique allows the manipulation of neuronal activity in locust neural circuits using optogenetics.

Cortical acetylcholine dynamics are predicted by cholinergic axon activity and behavior state.

Even under spontaneous conditions and in the absence of changing environmental demands, awake animals alternate between increased or decreased periods of alertness. These changes in brain state can occur rapidly, on a timescale of seconds, and neuromodulators such as acetylcholine (ACh) are thought to play an important role in driving these spontaneous state transitions. Here, we perform the first simultaneous imaging of ACh sensors and GCaMP-expressing axons in vivo, to examine the spatiotemporal properties of cortical ACh activity and release during spontaneous changes in behavioral state. We observed a high correlation between simultaneously recorded basal forebrain axon activity and neuromodulator sensor fluorescence around periods of locomotion and pupil dilation. Consistent with volume transmission of ACh, increases in axon activity were accompanied by increases in local ACh levels that fell off with the distance from the nearest axon. GRAB-ACh fluorescence could be accurately predicted from axonal activity alone, providing the first validation that neuromodulator axon activity is a reliable proxy for nearby neuromodulator levels. Deconvolution of fluorescence traces allowed us to account for the kinetics of the GRAB-ACh sensor and emphasized the rapid clearance of ACh for smaller transients outside of running periods. Finally, we trained a predictive model of ACh fluctuations from the combination of pupil size and running speed; this model performed better than using either variable alone, and generalized well to unseen data. Overall, these results contribute to a growing understanding of the precise timing and spatial characteristics of cortical ACh during fast brain state transitions.

Cortical acetylcholine dynamics are predicted by cholinergic axon activity and behavior state.

Acetylcholine (ACh) is thought to play a role in driving the rapid, spontaneous brain-state transitions that occur during wakefulness; however, the spatiotemporal properties of cortical ACh activity during these state changes are still unclear. We perform simultaneous imaging of GRAB-ACh sensors, GCaMP-expressing basal forebrain axons, and behavior to address this question. We observed a high correlation between axon and GRAB-ACh activity around periods of locomotion and pupil dilation. GRAB-ACh fluorescence could be accurately predicted from axonal activity alone, and local ACh activity decreased at farther distances from an axon. Deconvolution of GRAB-ACh traces allowed us to account for sensor kinetics and emphasized rapid clearance of small ACh transients. We trained a model to predict ACh from pupil size and running speed, which generalized well to unseen data. These results contribute to a growing understanding of the precise timing and spatial characteristics of cortical ACh during fast brain-state transitions.

Cell-type-specific inhibitory circuitry from a connectomic census of mouse visual cortex.

Mammalian cortex features a vast diversity of neuronal cell types, each with characteristic anatomical, molecular and functional properties. Synaptic connectivity powerfully shapes how each cell type participates in the cortical circuit, but mapping connectivity rules at the resolution of distinct cell types remains difficult. Here, we used millimeter-scale volumetric electron microscopy1 to investigate the connectivity of all inhibitory neurons across a densely-segmented neuronal population of 1352 cells spanning all layers of mouse visual cortex, producing a wiring diagram of inhibitory connections with more than 70,000 synapses. Taking a data-driven approach inspired by classical neuroanatomy, we classified inhibitory neurons based on the relative targeting of dendritic compartments and other inhibitory cells and developed a novel classification of excitatory neurons based on the morphological and synaptic input properties. The synaptic connectivity between inhibitory cells revealed a novel class of disinhibitory specialist targeting basket cells, in addition to familiar subclasses. Analysis of the inhibitory connectivity onto excitatory neurons found widespread specificity, with many interneurons exhibiting differential targeting of certain subpopulations spatially intermingled with other potential targets. Inhibitory targeting was organized into "motif groups," diverse sets of cells that collectively target both perisomatic and dendritic compartments of the same excitatory targets. Collectively, our analysis identified new organizing principles for cortical inhibition and will serve as a foundation for linking modern multimodal neuronal atlases with the cortical wiring diagram.

Functional connectivity and tuning curves in populations of simultaneously recorded neurons.

How interactions between neurons relate to tuned neural responses is a longstanding question in systems neuroscience. Here we use statistical modeling and simultaneous multi-electrode recordings to explore the relationship between these interactions and tuning curves in six different brain areas. We find that, in most cases, functional interactions between neurons provide an explanation of spiking that complements and, in some cases, surpasses the influence of canonical tuning curves. Modeling functional interactions improves both encoding and decoding accuracy by accounting for noise correlations and features of the external world that tuning curves fail to capture. In cortex, modeling coupling alone allows spikes to be predicted more accurately than tuning curve models based on external variables. These results suggest that statistical models of functional interactions between even relatively small numbers of neurons may provide a useful framework for examining neural coding.

Conserved neural dynamics and computations across species in olfaction.

Interpreting chemical information and translating it into ethologically relevant output is a common challenge of olfactory systems across species. Are computations performed by olfactory circuits conserved across species to overcome these common challenges? To understand this, we compared odor responses in the locust antennal lobe (AL) and mouse olfactory bulb (OB). We found that odors activated nearly mutually exclusive neural ensembles during stimulus presentation ('ON response') and after stimulus termination ('OFF response'). Strikingly, ON and OFF responses evoked by a single odor were anticorrelated with each other. 'Inverted' OFF responses enhanced contrast between odors experienced close together in time. Notably, OFF responses persisted long after odor termination in both AL and OB networks, indicating a form of short-term memory. Taken together, our results reveal key neurodynamic features underlying olfactory computations that are conserved across insect and mammalian olfactory systems.

Bipartite invariance in mouse primary visual cortex.

A defining characteristic of intelligent systems, whether natural or artificial, is the ability to generalize and infer behaviorally relevant latent causes from high-dimensional sensory input, despite significant variations in the environment. To understand how brains achieve generalization, it is crucial to identify the features to which neurons respond selectively and invariantly. However, the high-dimensional nature of visual inputs, the non-linearity of information processing in the brain, and limited experimental time make it challenging to systematically characterize neuronal tuning and invariances, especially for natural stimuli. Here, we extended "inception loops" - a paradigm that iterates between large-scale recordings, neural predictive models, and in silico experiments followed by in vivo verification - to systematically characterize single neuron invariances in the mouse primary visual cortex. Using the predictive model we synthesized Diverse Exciting Inputs (DEIs), a set of inputs that differ substantially from each other while each driving a target neuron strongly, and verified these DEIs' efficacy in vivo. We discovered a novel bipartite invariance: one portion of the receptive field encoded phase-invariant texture-like patterns, while the other portion encoded a fixed spatial pattern. Our analysis revealed that the division between the fixed and invariant portions of the receptive fields aligns with object boundaries defined by spatial frequency differences present in highly activating natural images. These findings suggest that bipartite invariance might play a role in segmentation by detecting texture-defined object boundaries, independent of the phase of the texture. We also replicated these bipartite DEIs in the functional connectomics MICrONs data set, which opens the way towards a circuit-level mechanistic understanding of this novel type of invariance. Our study demonstrates the power of using a data-driven deep learning approach to systematically characterize neuronal invariances. By applying this method across the visual hierarchy, cell types, and sensory modalities, we can decipher how latent variables are robustly extracted from natural scenes, leading to a deeper understanding of generalization.

Pattern completion and disruption characterize contextual modulation in mouse visual cortex.

A key role of sensory processing is integrating information across space. Neuronal responses in the visual system are influenced by both local features in the receptive field center and contextual information from the surround. While center-surround interactions have been extensively studied using simple stimuli like gratings, investigating these interactions with more complex, ecologically-relevant stimuli is challenging due to the high dimensionality of the stimulus space. We used large-scale neuronal recordings in mouse primary visual cortex to train convolutional neural network (CNN) models that accurately predicted center-surround interactions for natural stimuli. These models enabled us to synthesize surround stimuli that strongly suppressed or enhanced neuronal responses to the optimal center stimulus, as confirmed by in vivo experiments. In contrast to the common notion that congruent center and surround stimuli are suppressive, we found that excitatory surrounds appeared to complete spatial patterns in the center, while inhibitory surrounds disrupted them. We quantified this effect by demonstrating that CNN-optimized excitatory surround images have strong similarity in neuronal response space with surround images generated by extrapolating the statistical properties of the center, and with patches of natural scenes, which are known to exhibit high spatial correlations. Our findings cannot be explained by theories like redundancy reduction or predictive coding previously linked to contextual modulation in visual cortex. Instead, we demonstrated that a hierarchical probabilistic model incorporating Bayesian inference, and modulating neuronal responses based on prior knowledge of natural scene statistics, can explain our empirical results. We replicated these center-surround effects in the multi-area functional connectomics MICrONS dataset using natural movies as visual stimuli, which opens the way towards understanding circuit level mechanism, such as the contributions of lateral and feedback recurrent connections. Our data-driven modeling approach provides a new understanding of the role of contextual interactions in sensory processing and can be adapted across brain areas, sensory modalities, and species.

Towards a Foundation Model of the Mouse Visual Cortex.

Understanding the brain's perception algorithm is a highly intricate problem, as the inherent complexity of sensory inputs and the brain's nonlinear processing make characterizing sensory representations difficult. Recent studies have shown that functional models-capable of predicting large-scale neuronal activity in response to arbitrary sensory input-can be powerful tools for characterizing neuronal representations by enabling high-throughput in silico experiments. However, accurately modeling responses to dynamic and ecologically relevant inputs like videos remains challenging, particularly when generalizing to new stimulus domains outside the training distribution. Inspired by recent breakthroughs in artificial intelligence, where foundation models-trained on vast quantities of data-have demonstrated remarkable capabilities and generalization, we developed a "foundation model" of the mouse visual cortex: a deep neural network trained on large amounts of neuronal responses to ecological videos from multiple visual cortical areas and mice. The model accurately predicted neuronal responses not only to natural videos but also to various new stimulus domains, such as coherent moving dots and noise patterns, underscoring its generalization abilities. The foundation model could also be adapted to new mice with minimal natural movie training data. We applied the foundation model to the MICrONS dataset: a study of the brain that integrates structure with function at unprecedented scale, containing nanometer-scale morphology, connectivity with >500,000,000 synapses, and function of >70,000 neurons within a ~1mm3 volume spanning multiple areas of the mouse visual cortex. This accurate functional model of the MICrONS data opens the possibility for a systematic characterization of the relationship between circuit structure and function. By precisely capturing the response properties of the visual cortex and generalizing to new stimulus domains and mice, foundation models can pave the way for a deeper understanding of visual computation.

Distinct organization of two cortico-cortical feedback pathways.

Neocortical feedback is critical for attention, prediction, and learning. To mechanically understand its function requires deciphering its cell-type wiring. Recent studies revealed that feedback between primary motor to primary somatosensory areas in mice is disinhibitory, targeting vasoactive intestinal peptide-expressing interneurons, in addition to pyramidal cells. It is unknown whether this circuit motif represents a general cortico-cortical feedback organizing principle. Here we show that in contrast to this wiring rule, feedback between higher-order lateromedial visual area to primary visual cortex preferentially activates somatostatin-expressing interneurons. Functionally, both feedback circuits temporally sharpen feed-forward excitation eliciting a transient increase-followed by a prolonged decrease-in pyramidal cell activity under sustained feed-forward input. However, under feed-forward transient input, the primary motor to primary somatosensory cortex feedback facilitates bursting while lateromedial area to primary visual cortex feedback increases time precision. Our findings argue for multiple cortico-cortical feedback motifs implementing different dynamic non-linear operations.

Binary and analog variation of synapses between cortical pyramidal neurons.

Learning from experience depends at least in part on changes in neuronal connections. We present the largest map of connectivity to date between cortical neurons of a defined type (layer 2/3 [L2/3] pyramidal cells in mouse primary visual cortex), which was enabled by automated analysis of serial section electron microscopy images with improved handling of image defects (250 × 140 × 90 µm3 volume). We used the map to identify constraints on the learning algorithms employed by the cortex. Previous cortical studies modeled a continuum of synapse sizes by a log-normal distribution. A continuum is consistent with most neural network models of learning, in which synaptic strength is a continuously graded analog variable. Here, we show that synapse size, when restricted to synapses between L2/3 pyramidal cells, is well modeled by the sum of a binary variable and an analog variable drawn from a log-normal distribution. Two synapses sharing the same presynaptic and postsynaptic cells are known to be correlated in size. We show that the binary variables of the two synapses are highly correlated, while the analog variables are not. Binary variation could be the outcome of a Hebbian or other synaptic plasticity rule depending on activity signals that are relatively uniform across neuronal arbors, while analog variation may be dominated by other influences such as spontaneous dynamical fluctuations. We discuss the implications for the longstanding hypothesis that activity-dependent plasticity switches synapses between bistable states.

Periodicity and evoked responses in motor cortex.

Spiking in primary motor cortex (MI) exhibits a characteristic beta frequency periodicity, but the functional relevance of this rhythmic firing is controversial. We simultaneously recorded multiple single units and local field potentials in MI in two monkeys (Macaca mulatta) during continuous, self-paced movements to serially presented targets. We find that the appearance of each new target evokes precisely timed spiking in MI at a characteristic latency but that the exact timing of this response varies depending on its relationship to the phase of the ongoing beta range oscillation. As a result of this interaction between evoked spiking and endogenous beta periodicity, we find that the amount of information about target location encoded in the spiking of MI neurons is not simply a function of elapsed time but depends also on oscillatory phase. Our results suggest that periodicity may be an important feature of the early stages of sensorimotor processing in the cortical motor system.