The structure of axonal arbors controls how signals from individual neurons are routed within the mammalian brain. However, the arbors of very few long-range projection neurons have been reconstructed in their entirety, as axons with diameters as small as 100 nm arborize in target regions dispersed over many millimeters of tissue. We introduce a platform for high-resolution, three-dimensional fluorescence imaging of complete tissue volumes that enables the visualization and reconstruction of long-range axonal arbors. This platform relies on a high-speed two-photon microscope integrated with a tissue vibratome and a suite of computational tools for large-scale image data. We demonstrate the power of this approach by reconstructing the axonal arbors of multiple neurons in the motor cortex across a single mouse brain.

The Neurodata Without Borders (NWB) initiative promotes data standardization in neuroscience to increase research reproducibility and opportunities. In the first NWB pilot project, neurophysiologists and software developers produced a common data format for recordings and metadata of cellular electrophysiology and optical imaging experiments. The format specification, application programming interfaces, and sample datasets have been released.

Neural activity maintains representations that bridge past and future events, often over many seconds. Network models can produce persistent and ramping activity, but the positive feedback that is critical for these slow dynamics can cause sensitivity to perturbations. Here we use electrophysiology and optogenetic perturbations in the mouse premotor cortex to probe the robustness of persistent neural representations during motor planning. We show that preparatory activity is remarkably robust to large-scale unilateral silencing: detailed neural dynamics that drive specific future movements were quickly and selectively restored by the network. Selectivity did not recover after bilateral silencing of the premotor cortex. Perturbations to one hemisphere are thus corrected by information from the other hemisphere. Corpus callosum bisections demonstrated that premotor cortex hemispheres can maintain preparatory activity independently. Redundancy across selectively coupled modules, as we observed in the premotor cortex, is a hallmark of robust control systems. Network models incorporating these principles show robustness that is consistent with data.

The anterolateral motor cortex (ALM) and ventral medial (VM) thalamus are functionally linked to support persistent activity during motor planning. We analyzed the underlying synaptic interconnections using optogenetics and electrophysiology in mice (female/male). In cortex, thalamocortical (TC) axons from VM thalamus excited VM-projecting pyramidal tract (PT) neurons in layer 5B of ALM. These axons also strongly excited layer 2/3 neurons (which strongly excite PT neurons, as previously shown) but not VM-projecting corticothalamic (CT) neurons in layer 6. The strongest connections in the VM → PT circuit were localized to apical tuft dendrites of PT neurons, in layer 1. These tuft inputs were selectively augmented after blocking hyperpolarization-activated cyclic nucleotide-gated (HCN) channels. In thalamus, axons from ALM PT neurons excited ALM-projecting VM neurons, located medially in VM. These axons provided weak input to neurons in mediodorsal nucleus, and little or no input either to neurons in the GABAergic reticular thalamic nucleus or to neurons in VM projecting to primary motor cortex (M1). Conversely, M1 PT axons excited M1- but not ALM-projecting VM neurons. Our findings indicate, first, a set of cell type-specific connections forming an excitatory thalamo-cortico-thalamic loop for ALM ↔ VM communication and a circuit-level substrate for supporting reverberant activity in this system. Second, a key feature of this loop is the prominent involvement of layer 1 synapses onto apical dendrites, a subcellular compartment with distinct signaling properties, including HCN-mediated gain control. Third, the segregation of the ALM ↔ VM loop from M1-related circuits of VM adds cellular-level support for the concept of parallel pathway organization in the motor system.SIGNIFICANCE STATEMENT Anterolateral motor cortex (ALM), a higher-order motor area in the mouse, and ventromedial (VM) thalamus are anatomically and functionally linked, but their synaptic interconnections at the cellular level are unknown. Our results show that ALM pyramidal tract neurons monosynaptically excite ALM-projecting thalamocortical neurons in a medial subdivision of VM thalamus, and vice versa. The thalamo-cortico-thalamic loop formed by these recurrent connections constitutes a circuit-level substrate for supporting reverberant activity in this system.

Most functional plasticity studies in the cortex have focused on layers (L) II/III and IV, whereas relatively little is known of LV. Structural measurements of dendritic spines in vivo suggest some specialization among LV cell subtypes. We therefore studied experience-dependent plasticity in the barrel cortex using intracellular recordings to distinguish regular spiking (RS) and intrinsic bursting (IB) subtypes. Postsynaptic potentials and suprathreshold responses in vivo revealed a remarkable dichotomy in RS and IB cell plasticity; spared whisker potentiation occurred in IB but not RS cells while deprived whisker depression occurred in RS but not IB cells. Similar RS/IB differences were found in the LII/III to V connections in brain slices. Modeling studies showed that subthreshold changes predicted the suprathreshold changes. These studies demonstrate the major functional partition of plasticity within a single cortical layer and reveal the LII/III to LV connection as a major excitatory locus of cortical plasticity.

The primary auditory cortex (A1) is organized tonotopically, with neurons sensitive to high and low frequencies arranged in a rostro-caudal gradient. We used laser scanning photostimulation in acute slices to study the organization of local excitatory connections onto layers 2 and 3 (L2/3) of the mouse A1. Consistent with the organization of other cortical regions, synaptic inputs along the isofrequency axis (orthogonal to the tonotopic axis) arose predominantly within a column. By contrast, we found that local connections along the tonotopic axis differed from those along the isofrequency axis: some input pathways to L3 (but not L2) arose predominantly out-of-column. In vivo cell-attached recordings revealed differences between the sound-responsiveness of neurons in L2 and L3. Our results are consistent with the hypothesis that auditory cortical microcircuitry is specialized to the one-dimensional representation of frequency in the auditory cortex.

Spine growth and retraction with synapse formation and elimination plays an important role in shaping brain circuits during development and in the adult brain, yet the temporal relationship between spine morphogenesis and the formation of functional synapses remains poorly defined. We imaged hippocampal pyramidal neurons to identify spines of different ages. We then used two-photon glutamate uncaging, whole-cell recording, and Ca(2+) imaging to analyze the properties of nascent spines and their older neighbors. New spines expressed glutamate-sensitive currents that were indistinguishable from mature spines of comparable volumes. Some spines exhibited negligible AMPA receptor-mediated responses, but the occurrence of these "silent" spines was uncorrelated with spine age. In contrast, NMDA receptor-mediated Ca(2+) accumulations were significantly lower in new spines. New spines reconstructed using electron microscopy made synapses. Our data support a model in which outgrowth and enlargement of nascent spines is tightly coupled to formation and maturation of glutamatergic synapses.

Cortical maps, consisting of orderly arrangements of functional columns, are a hallmark of the organization of the cerebral cortex. However, the microorganization of cortical maps at the level of single neurons is not known, mainly because of the limitations of available mapping techniques. Here, we used bulk loading of Ca(2+) indicators combined with two-photon microscopy to image the activity of multiple single neurons in layer (L) 2/3 of the mouse barrel cortex in vivo. We developed methods that reliably detect single action potentials in approximately half of the imaged neurons in L2/3. This allowed us to measure the spiking probability following whisker deflection and thus map the whisker selectivity for multiple neurons with known spatial relationships. At the level of neuronal populations, the whisker map varied smoothly across the surface of the cortex, within and between the barrels. However, the whisker selectivity of individual neurons recorded simultaneously differed greatly, even for nearest neighbors. Trial-to-trial correlations between pairs of neurons were high over distances spanning multiple cortical columns. Our data suggest that the response properties of individual neurons are shaped by highly specific subcolumnar circuits and the momentary intrinsic state of the neocortex.

Primary sensory cortical areas receive information through multiple thalamic channels. In the rodent whisker system, lemniscal and paralemniscal thalamocortical projections, from the ventral posteromedial nucleus (VPM) and posterior medial nucleus (POm) respectively, carry distinct types of sensory information to cortex. Little is known about how these separate streams of activity are parsed and integrated within the neocortical microcircuit. We used quantitative laser scanning photostimulation to probe the organization of functional thalamocortical and ascending intracortical projections in the mouse barrel cortex. To map the thalamocortical projections, we recorded from neocortical excitatory neurons while stimulating VPM or POm. Neurons in layers (L)4, L5, and L6A received dense input from thalamus (L4, L5B, and L6A from VPM; and L5A from POm), whereas L2/3 neurons rarely received thalamic input. We further mapped the lemniscal and paralemniscal circuits from L4 and L5A to L2/3. Lemniscal L4 neurons targeted L3 within a column. Paralemniscal L5A neurons targeted a superficial band (thickness, 60 mum) of neurons immediately below L1, defining a functionally distinct L2 in the mouse barrel cortex. L2 neurons received input from lemniscal L3 cells and paralemniscal L5A cells spread over multiple columns. Our data indicate that lemniscal and paralemniscal information is segregated into interdigitated cortical layers.

Ca(2+) influx through NMDA receptors (NMDA-Rs) triggers synaptic plasticity, gene transcription, and cytotoxicity, but little is known about the regulation of NMDA-Rs themselves. We used two-photon glutamate uncaging to activate NMDA-Rs on individual dendritic spines in rat CA1 neurons while we measured NMDA-R currents at the soma and [Ca(2+)] changes in spines. Low-frequency uncaging trains induced Ca(2+)-dependent long-term depression of NMDA-R-mediated synaptic currents. Additionally, uncaging trains caused a reduction in the Ca(2+) accumulation per unit of NMDA-R current in spines due to a reduction in the fraction of the NMDA-R current carried by Ca(2+). Induction of depression of NMDA-R-mediated Ca(2+) influx required activation of NR2B-containing receptors. Receptors in single spines depressed rapidly in an all-or-none manner. These adaptive changes in NMDA-R function likely play a critical role in metaplasticity and in stabilizing activity levels in neuronal networks with Hebbian synapses.

During active somatosensation, neural signals expected from movement of the sensors are suppressed in the cortex, whereas information related to touch is enhanced. This tactile suppression underlies low-noise encoding of relevant tactile features and the brain's ability to make fine tactile discriminations. Layer (L) 4 excitatory neurons in the barrel cortex, the major target of the somatosensory thalamus (VPM), respond to touch, but have low spike rates and low sensitivity to the movement of whiskers. Most neurons in VPM respond to touch and also show an increase in spike rate with whisker movement. Therefore, signals related to self-movement are suppressed in L4. Fast-spiking (FS) interneurons in L4 show similar dynamics to VPM neurons. Stimulation of halorhodopsin in FS interneurons causes a reduction in FS neuron activity and an increase in L4 excitatory neuron activity. This decrease of activity of L4 FS neurons contradicts the "paradoxical effect" predicted in networks stabilized by inhibition and in strongly-coupled networks. To explain these observations, we constructed a model of the L4 circuit, with connectivity constrained by in vitro measurements. The model explores the various synaptic conductance strengths for which L4 FS neurons actively suppress baseline and movement-related activity in layer 4 excitatory neurons. Feedforward inhibition, in concert with recurrent intracortical circuitry, produces tactile suppression. Synaptic delays in feedforward inhibition allow transmission of temporally brief volleys of activity associated with touch. Our model provides a mechanistic explanation of a behavior-related computation implemented by the thalamocortical circuit.

The neocortex contains a multitude of cell types that are segregated into layers and functionally distinct areas. To investigate the diversity of cell types across the mouse neocortex, here we analysed 23,822 cells from two areas at distant poles of the mouse neocortex: the primary visual cortex and the anterior lateral motor cortex. We define 133 transcriptomic cell types by deep, single-cell RNA sequencing. Nearly all types of GABA (γ-aminobutyric acid)-containing neurons are shared across both areas, whereas most types of glutamatergic neurons were found in one of the two areas. By combining single-cell RNA sequencing and retrograde labelling, we match transcriptomic types of glutamatergic neurons to their long-range projection specificity. Our study establishes a combined transcriptomic and projectional taxonomy of cortical cell types from functionally distinct areas of the adult mouse cortex.

Point-scanning two-photon microscopy enables high-resolution imaging within scattering specimens such as the mammalian brain, but sequential acquisition of voxels fundamentally limits its speed. We developed a two-photon imaging technique that scans lines of excitation across a focal plane at multiple angles and computationally recovers high-resolution images, attaining voxel rates of over 1 billion Hz in structured samples. Using a static image as a prior for recording neural activity, we imaged visually evoked and spontaneous glutamate release across hundreds of dendritic spines in mice at depths over 250 µm and frame rates over 1 kHz. Dendritic glutamate transients in anesthetized mice are synchronized within spatially contiguous domains spanning tens of micrometers at frequencies ranging from 1-100 Hz. We demonstrate millisecond-resolved recordings of acetylcholine and voltage indicators, three-dimensional single-particle tracking and imaging in densely labeled cortex. Our method surpasses limits on the speed of raster-scanned imaging imposed by fluorescence lifetime.

Task-related neural activity is widespread across populations of neurons during goal-directed behaviors. However, little is known about the synaptic reorganization and circuit mechanisms that lead to broad activity changes. Here we trained a subset of neurons in a spiking network with strong synaptic interactions to reproduce the activity of neurons in the motor cortex during a decision-making task. Task-related activity, resembling the neural data, emerged across the network, even in the untrained neurons. Analysis of trained networks showed that strong untrained synapses, which were independent of the task and determined the dynamical state of the network, mediated the spread of task-related activity. Optogenetic perturbations suggest that the motor cortex is strongly-coupled, supporting the applicability of the mechanism to cortical networks. Our results reveal a cortical mechanism that facilitates distributed representations of task-variables by spreading the activity from a subset of plastic neurons to the entire network through task-independent strong synapses.

Neuronal cell types are the nodes of neural circuits that determine the flow of information within the brain. Neuronal morphology, especially the shape of the axonal arbor, provides an essential descriptor of cell type and reveals how individual neurons route their output across the brain. Despite the importance of morphology, few projection neurons in the mouse brain have been reconstructed in their entirety. Here we present a robust and efficient platform for imaging and reconstructing complete neuronal morphologies, including axonal arbors that span substantial portions of the brain. We used this platform to reconstruct more than 1,000 projection neurons in the motor cortex, thalamus, subiculum, and hypothalamus. Together, the reconstructed neurons constitute more than 85 meters of axonal length and are available in a searchable online database. Axonal shapes revealed previously unknown subtypes of projection neurons and suggest organizational principles of long-range connectivity.

Genetically encoded voltage indicators (GEVIs) enable monitoring of neuronal activity at high spatial and temporal resolution. However, the utility of existing GEVIs has been limited by the brightness and photostability of fluorescent proteins and rhodopsins. We engineered a GEVI, called Voltron, that uses bright and photostable synthetic dyes instead of protein-based fluorophores, thereby extending the number of neurons imaged simultaneously in vivo by a factor of 10 and enabling imaging for significantly longer durations relative to existing GEVIs. We used Voltron for in vivo voltage imaging in mice, zebrafish, and fruit flies. In the mouse cortex, Voltron allowed single-trial recording of spikes and subthreshold voltage signals from dozens of neurons simultaneously over a 15-minute period of continuous imaging. In larval zebrafish, Voltron enabled the precise correlation of spike timing with behavior.

The subcellular locations of synapses on pyramidal neurons strongly influences dendritic integration and synaptic plasticity. Despite this, there is little quantitative data on spatial distributions of specific types of synaptic input. Here we use array tomography (AT), a high-resolution optical microscopy method, to examine thalamocortical (TC) input onto layer 5 pyramidal neurons. We first verified the ability of AT to identify synapses using parallel electron microscopic analysis of TC synapses in layer 4. We then use large-scale array tomography (LSAT) to measure TC synapse distribution on L5 pyramidal neurons in a 1.00 × 0.83 × 0.21 mm(3) volume of mouse somatosensory cortex. We found that TC synapses primarily target basal dendrites in layer 5, but also make a considerable input to proximal apical dendrites in L4, consistent with previous work. Our analysis further suggests that TC inputs are biased toward certain branches and, within branches, synapses show significant clustering with an excess of TC synapse nearest neighbors within 5-15 μm compared to a random distribution. Thus, we show that AT is a sensitive and quantitative method to map specific types of synaptic input on the dendrites of entire neurons. We anticipate that this technique will be of wide utility for mapping functionally-relevant anatomical connectivity in neural circuits.

Fluorescent protein-based sensors for detecting neuronal activity have been developed largely based on non-neuronal screening systems. However, the dynamics of neuronal state variables (e.g., voltage, calcium, etc.) are typically very rapid compared to those of non-excitable cells. We developed an electrical stimulation and fluorescence imaging platform based on dissociated rat primary neuronal cultures. We describe its use in testing genetically-encoded calcium indicators (GECIs). Efficient neuronal GECI expression was achieved using lentiviruses containing a neuronal-selective gene promoter. Action potentials (APs) and thus neuronal calcium levels were quantitatively controlled by electrical field stimulation, and fluorescence images were recorded. Images were segmented to extract fluorescence signals corresponding to individual GECI-expressing neurons, which improved sensitivity over full-field measurements. We demonstrate the superiority of screening GECIs in neurons compared with solution measurements. Neuronal screening was useful for efficient identification of variants with both improved response kinetics and high signal amplitudes. This platform can be used to screen many types of sensors with cellular resolution under realistic conditions where neuronal state variables are in relevant ranges with respect to timing and amplitude.

The mouse is an increasingly prominent model for the analysis of mammalian neuronal circuits. Neural circuits ultimately have to be probed during behaviors that engage the circuits. Linking circuit dynamics to behavior requires precise control of sensory stimuli and measurement of body movements. Head-fixation has been used for behavioral research, particularly in non-human primates, to facilitate precise stimulus control, behavioral monitoring and neural recording. However, choice-based, perceptual decision tasks by head-fixed mice have only recently been introduced. Training mice relies on motivating mice using water restriction. Here we describe procedures for head-fixation, water restriction and behavioral training for head-fixed mice, with a focus on active, whisker-based tactile behaviors. In these experiments mice had restricted access to water (typically 1 ml/day). After ten days of water restriction, body weight stabilized at approximately 80% of initial weight. At that point mice were trained to discriminate sensory stimuli using operant conditioning. Head-fixed mice reported stimuli by licking in go/no-go tasks and also using a forced choice paradigm using a dual lickport. In some cases mice learned to discriminate sensory stimuli in a few trials within the first behavioral session. Delay epochs lasting a second or more were used to separate sensation (e.g. tactile exploration) and action (i.e. licking). Mice performed a variety of perceptual decision tasks with high performance for hundreds of trials per behavioral session. Up to four months of continuous water restriction showed no adverse health effects. Behavioral performance correlated with the degree of water restriction, supporting the importance of controlling access to water. These behavioral paradigms can be combined with cellular resolution imaging, random access photostimulation, and whole cell recordings.

Activity in the mouse anterior lateral motor cortex (ALM) instructs directional movements, often seconds before movement initiation. It is unknown whether this preparatory activity is localized to ALM or widely distributed within motor cortex. Here we imaged activity across motor cortex while mice performed a whisker-based object localization task with a delayed, directional licking response. During tactile sensation and the delay epoch, object location was represented in motor cortex areas that are medial and posterior relative to ALM, including vibrissal motor cortex. Preparatory activity appeared first in deep layers of ALM, seconds before the behavioral response, and remained localized to ALM until the behavioral response. Later, widely distributed neurons represented the outcome of the trial. Cortical area was more predictive of neuronal selectivity than laminar location or axonal projection target. Motor cortex therefore represents sensory, motor, and outcome information in a spatially organized manner.