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.

Oligodendrocyte precursor cells ingest axons in the mouse neocortex.

Neurons in the developing brain undergo extensive structural refinement as nascent circuits adopt their mature form. This physical transformation of neurons is facilitated by the engulfment and degradation of axonal branches and synapses by surrounding glial cells, including microglia and astrocytes. However, the small size of phagocytic organelles and the complex, highly ramified morphology of glia have made it difficult to define the contribution of these and other glial cell types to this crucial process. Here, we used large-scale, serial section transmission electron microscopy (TEM) with computational volume segmentation to reconstruct the complete 3D morphologies of distinct glial types in the mouse visual cortex, providing unprecedented resolution of their morphology and composition. Unexpectedly, we discovered that the fine processes of oligodendrocyte precursor cells (OPCs), a population of abundant, highly dynamic glial progenitors, frequently surrounded small branches of axons. Numerous phagosomes and phagolysosomes (PLs) containing fragments of axons and vesicular structures were present inside their processes, suggesting that OPCs engage in axon pruning. Single-nucleus RNA sequencing from the developing mouse cortex revealed that OPCs express key phagocytic genes at this stage, as well as neuronal transcripts, consistent with active axon engulfment. Although microglia are thought to be responsible for the majority of synaptic pruning and structural refinement, PLs were ten times more abundant in OPCs than in microglia at this stage, and these structures were markedly less abundant in newly generated oligodendrocytes, suggesting that OPCs contribute substantially to the refinement of neuronal circuits during cortical development.

Mapping synapses by conjugate light-electron array tomography.

Synapses of the mammalian CNS are diverse in size, structure, molecular composition, and function. Synapses in their myriad variations are fundamental to neural circuit development, homeostasis, plasticity, and memory storage. Unfortunately, quantitative analysis and mapping of the brain's heterogeneous synapse populations has been limited by the lack of adequate single-synapse measurement methods. Electron microscopy (EM) is the definitive means to recognize and measure individual synaptic contacts, but EM has only limited abilities to measure the molecular composition of synapses. This report describes conjugate array tomography (AT), a volumetric imaging method that integrates immunofluorescence and EM imaging modalities in voxel-conjugate fashion. We illustrate the use of conjugate AT to advance the proteometric measurement of EM-validated single-synapse analysis in a study of mouse cortex.

Ultrastructural differences impact cilia shape and external exposure across cell classes in the visual cortex.

A primary cilium is a membrane-bound extension from the cell surface that contains receptors for perceiving and transmitting signals that modulate cell state and activity. Primary cilia in the brain are less accessible than cilia on cultured cells or epithelial tissues because in the brain they protrude into a deep, dense network of glial and neuronal processes. Here, we investigated cilia frequency, internal structure, shape, and position in large, high-resolution transmission electron microscopy volumes of mouse primary visual cortex. Cilia extended from the cell bodies of nearly all excitatory and inhibitory neurons, astrocytes, and oligodendrocyte precursor cells (OPCs) but were absent from oligodendrocytes and microglia. Ultrastructural comparisons revealed that the base of the cilium and the microtubule organization differed between neurons and glia. Investigating cilia-proximal features revealed that many cilia were directly adjacent to synapses, suggesting that cilia are poised to encounter locally released signaling molecules. Our analysis indicated that synapse proximity is likely due to random encounters in the neuropil, with no evidence that cilia modulate synapse activity as would be expected in tetrapartite synapses. The observed cell class differences in proximity to synapses were largely due to differences in external cilia length. Many key structural features that differed between neuronal and glial cilia influenced both cilium placement and shape and, thus, exposure to processes and synapses outside the cilium. Together, the ultrastructure both within and around neuronal and glial cilia suggest differences in cilia formation and function across cell types in the brain.

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.

Emerging roles of oligodendrocyte precursor cells in neural circuit development and remodeling.

Oligodendrocyte precursor cells (OPCs) are non-neuronal brain cells that give rise to oligodendrocytes, glia that myelinate the axons of neurons in the brain. Classically known for their contributions to myelination via oligodendrogenesis, OPCs are increasingly appreciated to play diverse roles in the nervous system, ranging from blood vessel formation to antigen presentation. Here, we review emerging literature suggesting that OPCs may be essential for the establishment and remodeling of neural circuits in the developing and adult brain via mechanisms that are distinct from the production of oligodendrocytes. We discuss the specialized features of OPCs that position these cells to integrate activity-dependent and molecular cues to shape brain wiring. Finally, we place OPCs within the context of a growing field focused on understanding the importance of communication between neurons and glia in the contexts of both health and disease.

Nanometer-scale views of visual cortex reveal anatomical features of primary cilia poised to detect synaptic spillover.

A primary cilium is a thin membrane-bound extension off a cell surface that contains receptors for perceiving and transmitting signals that modulate cell state and activity. While many cell types have a primary cilium, little is known about primary cilia in the brain, where they are less accessible than cilia on cultured cells or epithelial tissues and protrude from cell bodies into a deep, dense network of glial and neuronal processes. Here, we investigated cilia frequency, internal structure, shape, and position in large, high-resolution transmission electron microscopy volumes of mouse primary visual cortex. Cilia extended from the cell bodies of nearly all excitatory and inhibitory neurons, astrocytes, and oligodendrocyte precursor cells (OPCs), but were absent from oligodendrocytes and microglia. Structural comparisons revealed that the membrane structure at the base of the cilium and the microtubule organization differed between neurons and glia. OPC cilia were distinct in that they were the shortest and contained pervasive internal vesicles only occasionally observed in neuron and astrocyte cilia. Investigating cilia-proximal features revealed that many cilia were directly adjacent to synapses, suggesting cilia are well poised to encounter locally released signaling molecules. Cilia proximity to synapses was random, not enriched, in the synapse-rich neuropil. The internal anatomy, including microtubule changes and centriole location, defined key structural features including cilium placement and shape. Together, the anatomical insights both within and around neuron and glia cilia provide new insights into cilia formation and function across cell types in the brain.

NEURD: automated proofreading and feature extraction for connectomics.

We are now in the era of millimeter-scale electron microscopy (EM) volumes collected at nanometer resolution (Shapson-Coe et al., 2021; Consortium et al., 2021). Dense reconstruction of cellular compartments in these EM volumes has been enabled by recent advances in Machine Learning (ML) (Lee et al., 2017; Wu et al., 2021; Lu et al., 2021; Macrina et al., 2021). Automated segmentation methods can now yield exceptionally accurate reconstructions of cells, but despite this accuracy, laborious post-hoc proofreading is still required to generate large connectomes free of merge and split errors. The elaborate 3-D meshes of neurons produced by these segmentations contain detailed morphological information, from the diameter, shape, and branching patterns of axons and dendrites, down to the fine-scale structure of dendritic spines. However, extracting information about these features can require substantial effort to piece together existing tools into custom workflows. Building on existing open-source software for mesh manipulation, here we present "NEURD", a software package that decomposes each meshed neuron into a compact and extensively-annotated graph representation. With these feature-rich graphs, we implement workflows for state of the art automated post-hoc proofreading of merge errors, cell classification, spine detection, axon-dendritic proximities, and other features that can enable many downstream analyses of neural morphology and connectivity. NEURD can make these new massive and complex datasets more accessible to neuroscience researchers focused on a variety of scientific questions.

Functional connectomics reveals general wiring rule in mouse visual cortex.

To understand how the brain computes, it is important to unravel the relationship between circuit connectivity and function. Previous research has shown that excitatory neurons in layer 2/3 of the primary visual cortex of mice with similar response properties are more likely to form connections. However, technical challenges of combining synaptic connectivity and functional measurements have limited these studies to few, highly local connections. Utilizing the millimeter scale and nanometer resolution of the MICrONS dataset, we studied the connectivity-function relationship in excitatory neurons of the mouse visual cortex across interlaminar and interarea projections, assessing connection selectivity at the coarse axon trajectory and fine synaptic formation levels. A digital twin model of this mouse, that accurately predicted responses to arbitrary video stimuli, enabled a comprehensive characterization of the function of neurons. We found that neurons with highly correlated responses to natural videos tended to be connected with each other, not only within the same cortical area but also across multiple layers and visual areas, including feedforward and feedback connections, whereas we did not find that orientation preference predicted connectivity. The digital twin model separated each neuron's tuning into a feature component (what the neuron responds to) and a spatial component (where the neuron's receptive field is located). We show that the feature, but not the spatial component, predicted which neurons were connected at the fine synaptic scale. Together, our results demonstrate the "like-to-like" connectivity rule generalizes to multiple connection types, and the rich MICrONS dataset is suitable to further refine a mechanistic understanding of circuit structure and function.

CAVE: Connectome Annotation Versioning Engine.

Advances in Electron Microscopy, image segmentation and computational infrastructure have given rise to large-scale and richly annotated connectomic datasets which are increasingly shared across communities. To enable collaboration, users need to be able to concurrently create new annotations and correct errors in the automated segmentation by proofreading. In large datasets, every proofreading edit relabels cell identities of millions of voxels and thousands of annotations like synapses. For analysis, users require immediate and reproducible access to this constantly changing and expanding data landscape. Here, we present the Connectome Annotation Versioning Engine (CAVE), a computational infrastructure for immediate and reproducible connectome analysis in up-to petascale datasets (~1mm3) while proofreading and annotating is ongoing. For segmentation, CAVE provides a distributed proofreading infrastructure for continuous versioning of large reconstructions. Annotations in CAVE are defined by locations such that they can be quickly assigned to the underlying segment which enables fast analysis queries of CAVE's data for arbitrary time points. CAVE supports schematized, extensible annotations, so that researchers can readily design novel annotation types. CAVE is already used for many connectomics datasets, including the largest datasets available to date.

A scalable and modular automated pipeline for stitching of large electron microscopy datasets.

Serial-section electron microscopy (ssEM) is the method of choice for studying macroscopic biological samples at extremely high resolution in three dimensions. In the nervous system, nanometer-scale images are necessary to reconstruct dense neural wiring diagrams in the brain, so -called connectomes. The data that can comprise of up to 108 individual EM images must be assembled into a volume, requiring seamless 2D registration from physical section followed by 3D alignment of the stitched sections. The high throughput of ssEM necessitates 2D stitching to be done at the pace of imaging, which currently produces tens of terabytes per day. To achieve this, we present a modular volume assembly software pipeline ASAP (Assembly Stitching and Alignment Pipeline) that is scalable to datasets containing petabytes of data and parallelized to work in a distributed computational environment. The pipeline is built on top of the Render Trautman and Saalfeld (2019) services used in the volume assembly of the brain of adult Drosophila melanogaster (Zheng et al. 2018). It achieves high throughput by operating only on image meta-data and transformations. ASAP is modular, allowing for easy incorporation of new algorithms without significant changes in the workflow. The entire software pipeline includes a complete set of tools for stitching, automated quality control, 3D section alignment, and final rendering of the assembled volume to disk. ASAP has been deployed for continuous stitching of several large-scale datasets of the mouse visual cortex and human brain samples including one cubic millimeter of mouse visual cortex (Yin et al. 2020); Microns Consortium et al. (2021) at speeds that exceed imaging. The pipeline also has multi-channel processing capabilities and can be applied to fluorescence and multi-modal datasets like array tomography.

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.

Ca2+ store depletion causes STIM1 to accumulate in ER regions closely associated with the plasma membrane.

Stromal interacting molecule 1 (STIM1), reported to be an endoplasmic reticulum (ER) Ca(2+) sensor controlling store-operated Ca(2+) entry, redistributes from a diffuse ER localization into puncta at the cell periphery after store depletion. STIM1 redistribution is proposed to be necessary for Ca(2+) release-activated Ca(2+) (CRAC) channel activation, but it is unclear whether redistribution is rapid enough to play a causal role. Furthermore, the location of STIM1 puncta is uncertain, with recent reports supporting retention in the ER as well as insertion into the plasma membrane (PM). Using total internal reflection fluorescence (TIRF) microscopy and patch-clamp recording from single Jurkat cells, we show that STIM1 puncta form several seconds before CRAC channels open, supporting a causal role in channel activation. Fluorescence quenching and electron microscopy analysis reveal that puncta correspond to STIM1 accumulation in discrete subregions of junctional ER located 10-25 nm from the PM, without detectable insertion of STIM1 into the PM. Roughly one third of these ER-PM contacts form in response to store depletion. These studies identify an ER structure underlying store-operated Ca(2+) entry, whose extreme proximity to the PM may enable STIM1 to interact with CRAC channels or associated proteins.

The elementary unit of store-operated Ca2+ entry: local activation of CRAC channels by STIM1 at ER-plasma membrane junctions.

The activation of store-operated Ca(2+) entry by Ca(2+) store depletion has long been hypothesized to occur via local interactions of the endoplasmic reticulum (ER) and plasma membrane, but the structure involved has never been identified. Store depletion causes the ER Ca(2+) sensor stromal interacting molecule 1 (STIM1) to form puncta by accumulating in junctional ER located 10-25 nm from the plasma membrane (see Wu et al. on p. 803 of this issue). We have combined total internal reflection fluorescence (TIRF) microscopy and patch-clamp recording to localize STIM1 and sites of Ca(2+) influx through open Ca(2+) release-activated Ca(2+) (CRAC) channels in Jurkat T cells after store depletion. CRAC channels open only in the immediate vicinity of STIM1 puncta, restricting Ca(2+) entry to discrete sites comprising a small fraction of the cell surface. Orai1, an essential component of the CRAC channel, colocalizes with STIM1 after store depletion, providing a physical basis for the local activation of Ca(2+) influx. These studies reveal for the first time that STIM1 and Orai1 move in a coordinated fashion to form closely apposed clusters in the ER and plasma membranes, thereby creating the elementary unit of store-operated Ca(2+) entry.

Structure and function of axo-axonic inhibition.

Inhibitory neurons in mammalian cortex exhibit diverse physiological, morphological, molecular, and connectivity signatures. While considerable work has measured the average connectivity of several interneuron classes, there remains a fundamental lack of understanding of the connectivity distribution of distinct inhibitory cell types with synaptic resolution, how it relates to properties of target cells, and how it affects function. Here, we used large-scale electron microscopy and functional imaging to address these questions for chandelier cells in layer 2/3 of the mouse visual cortex. With dense reconstructions from electron microscopy, we mapped the complete chandelier input onto 153 pyramidal neurons. We found that synapse number is highly variable across the population and is correlated with several structural features of the target neuron. This variability in the number of axo-axonic ChC synapses is higher than the variability seen in perisomatic inhibition. Biophysical simulations show that the observed pattern of axo-axonic inhibition is particularly effective in controlling excitatory output when excitation and inhibition are co-active. Finally, we measured chandelier cell activity in awake animals using a cell-type-specific calcium imaging approach and saw highly correlated activity across chandelier cells. In the same experiments, in vivo chandelier population activity correlated with pupil dilation, a proxy for arousal. Together, these results suggest that chandelier cells provide a circuit-wide signal whose strength is adjusted relative to the properties of target neurons.

The three-dimensional micro- and nanostructure of the aortic medial lamellar unit measured using 3D confocal and electron microscopy imaging.

Changes in arterial wall composition and function underlie all forms of vascular disease. The fundamental structural and functional unit of the aortic wall is the medial lamellar unit (MLU). While the basic composition and organization of the MLU is known, three-dimensional (3D) microstructural details are tenuous, due (in part) to lack of three-dimensional data at micro- and nano-scales. We applied novel electron and confocal microscopy techniques to obtain 3D volumetric information of aortic medial microstructure at micro- and nano-scales with all constituents present. For the rat abdominal aorta, we show that medial elastin has three primary forms: with approximately 71% of total elastin as thick, continuous lamellar sheets, 27% as thin, protruding interlamellar elastin fibers (IEFs), and 2% as thick radial struts. Elastin pores are not simply holes in lamellar sheets, but are indented and gusseted openings in lamellae. Smooth muscle cells (SMCs) weave throughout the interlamellar elastin framework, with cytoplasmic extensions abutting IEFs, resulting in approximately 20 degrees radial tilt (relative to the lumen surface) of elliptical SMC nuclei. Collagen fibers are organized as large, parallel bundles tightly enveloping SMC nuclei. Quantification of the orientation of collagen bundles, SMC nuclei, and IEFs reveal that all three primary medial constituents have predominantly circumferential orientation, correlating with reported circumferentially dominant values of physiological stress, collagen fiber recruitment, and tissue stiffness. This high resolution three-dimensional view of the aortic media reveals MLU microstructure details that suggest a highly complex and integrated mural organization that correlates with aortic mechanical properties.

Retrograde regulation of synaptic vesicle endocytosis and recycling.

Sustained release of neurotransmitter depends upon the recycling of synaptic vesicles. Until now, it has been assumed that vesicle recycling is regulated by signals from the presynaptic bouton alone, but results from rat hippocampal neurons reported here indicate that this need not be the case. Fluorescence imaging and pharmacological analysis show that a nitric oxide (NO) signal generated postsynaptically can regulate endocytosis and at least one later step in synaptic vesicle recycling. The proposed retrograde pathway involves an NMDA receptor (NMDAR)-dependent postsynaptic production of NO, diffusion of NO to a presynaptic site, and a cGMP-dependent increase in presynaptic phosphatidylinositol 4,5-biphosphate (PIP2). These results indicate that the regulation of synaptic vesicle recycling may integrate a much broader range of neural activity signals than previously recognized, including postsynaptic depolarization and the activation of NMDARs at both immediate and nearby postsynaptic active zones.

A large fraction of neocortical myelin ensheathes axons of local inhibitory neurons.

Myelin is best known for its role in increasing the conduction velocity and metabolic efficiency of long-range excitatory axons. Accordingly, the myelin observed in neocortical gray matter is thought to mostly ensheath excitatory axons connecting to subcortical regions and distant cortical areas. Using independent analyses of light and electron microscopy data from mouse neocortex, we show that a surprisingly large fraction of cortical myelin (half the myelin in layer 2/3 and a quarter in layer 4) ensheathes axons of inhibitory neurons, specifically of parvalbumin-positive basket cells. This myelin differs significantly from that of excitatory axons in distribution and protein composition. Myelin on inhibitory axons is unlikely to meaningfully hasten the arrival of spikes at their pre-synaptic terminals, due to the patchy distribution and short path-lengths observed. Our results thus highlight the need for exploring alternative roles for myelin in neocortical circuits.

Late recruitment of synapsin to nascent synapses is regulated by Cdk5.

Synapse formation is a complex process that involves the recruitment and assembly of a myriad of pre- and postsynaptic proteins. Despite being present at every synapse in the vertebrate CNS, little is known about the transport, recruitment, and stabilization of synapsin at nascent synapses during development. We examined the transport and recruitment of synapsin to nascent presynaptic terminals in vivo in the developing zebrafish spinal cord. Synapsin was transported in a transport packet independently of two other presynaptic organelles: synaptic vesicle (SV) protein transport vesicles (STVs) and Piccolo-containing active zone precursor transport vesicles (PTVs). During presynaptic assembly, recruitment of all three transport packets occurred in an ordered sequence: STVs preceded PTVs, which in turn preceded synapsin. Importantly, cyclin-dependent kinase 5 (Cdk5) specifically regulated the late recruitment of synapsin transport packets at synapses. These results point to additional layers of complexity in the established mechanisms of synaptogenesis.

High-contrast en bloc staining of neuronal tissue for field emission scanning electron microscopy.

Conventional heavy metal poststaining methods on thin sections lend contrast but often cause contamination. To avoid this problem, we tested several en bloc staining techniques to contrast tissue in serial sections mounted on solid substrates for examination by field emission scanning electron microscopy (FESEM). Because FESEM section imaging requires that specimens have higher contrast and greater electrical conductivity than transmission electron microscopy (TEM) samples, our technique uses osmium impregnation (OTO) to make the samples conductive while heavily staining membranes for segmentation studies. Combining this step with other classic heavy metal en bloc stains, including uranyl acetate (UA), lead aspartate, copper sulfate and lead citrate, produced clean, highly contrasted TEM and scanning electron microscopy (SEM) samples of insect, fish and mammalian nervous systems. This protocol takes 7-15 d to prepare resin-embedded tissue, cut sections and produce serial section images.