Cyclic structure with cellular precision in a vertebrate sensorimotor neural circuit.

Neuronal wiring diagrams reconstructed by electron microscopy1,2,3,4,5 pose new questions about the organization of nervous systems following the time-honored tradition of cross-species comparisons.6,7 The C. elegans connectome has been conceptualized as a sensorimotor circuit that is approximately feedforward,8,9,10,11 starting from sensory neurons proceeding to interneurons and ending with motor neurons. Overrepresentation of a 3-cell motif often known as the "feedforward loop" has provided further evidence for feedforwardness.10,12 Here, we contrast with another sensorimotor wiring diagram that was recently reconstructed from a larval zebrafish brainstem.13 We show that the 3-cycle, another 3-cell motif, is highly overrepresented in the oculomotor module of this wiring diagram. This is a first for any neuronal wiring diagram reconstructed by electron microscopy, whether invertebrate12,14 or mammalian.15,16,17 The 3-cycle of cells is "aligned" with a 3-cycle of neuronal groups in a stochastic block model (SBM)18 of the oculomotor module. However, the cellular cycles exhibit more specificity than can be explained by the group cycles-recurrence to the same neuron is surprisingly common. Cyclic structure could be relevant for theories of oculomotor function that depend on recurrent connectivity. The cyclic structure coexists with the classic vestibulo-ocular reflex arc for horizontal eye movements,19 and could be relevant for recurrent network models of temporal integration by the oculomotor system.20,21.

RealNeuralNetworks.jl: An Integrated Julia Package for Skeletonization, Morphological Analysis, and Synaptic Connectivity Analysis of Terabyte-Scale 3D Neural Segmentations.

Benefiting from the rapid development of electron microscopy imaging and deep learning technologies, an increasing number of brain image datasets with segmentation and synapse detection are published. Most of the automated segmentation methods label voxels rather than producing neuron skeletons directly. A further skeletonization step is necessary for quantitative morphological analysis. Currently, several tools are published for skeletonization as well as morphological and synaptic connectivity analysis using different computer languages and environments. Recently the Julia programming language, notable for elegant syntax and high performance, has gained rapid adoption in the scientific computing community. Here, we present a Julia package, called RealNeuralNetworks.jl, for efficient sparse skeletonization, morphological analysis, and synaptic connectivity analysis. Based on a large-scale Zebrafish segmentation dataset, we illustrate the software features by performing distributed skeletonization in Google Cloud, clustering the neurons using the NBLAST algorithm, combining morphological similarity and synaptic connectivity to study their relationship. We demonstrate that RealNeuralNetworks.jl is suitable for use in terabyte-scale electron microscopy image segmentation datasets.

Methods for Mapping Neuronal Activity to Synaptic Connectivity: Lessons From Larval Zebrafish.

For a mechanistic understanding of neuronal circuits in the brain, a detailed description of information flow is necessary. Thereby it is crucial to link neuron function to the underlying circuit structure. Multiphoton calcium imaging is the standard technique to record the activity of hundreds of neurons simultaneously. Similarly, recent advances in high-throughput electron microscopy techniques allow for the reconstruction of synaptic resolution wiring diagrams. These two methods can be combined to study both function and structure in the same specimen. Due to its small size and optical transparency, the larval zebrafish brain is one of the very few vertebrate systems where both, activity and connectivity of all neurons from entire, anatomically defined brain regions, can be analyzed. Here, we describe different methods and the tools required for combining multiphoton microscopy with dense circuit reconstruction from electron microscopy stacks of entire brain regions in the larval zebrafish.

Electron Microscopic Reconstruction of Functionally Identified Cells in a Neural Integrator.

Neural integrators are involved in a variety of sensorimotor and cognitive behaviors. The oculomotor system contains a simple example, a hindbrain neural circuit that takes velocity signals as inputs and temporally integrates them to control eye position. Here we investigated the structural underpinnings of temporal integration in the larval zebrafish by first identifying integrator neurons using two-photon calcium imaging and then reconstructing the same neurons through serial electron microscopic analysis. Integrator neurons were identified as those neurons with activities highly correlated with eye position during spontaneous eye movements. Three morphological classes of neurons were observed: ipsilaterally projecting neurons located medially, contralaterally projecting neurons located more laterally, and a population at the extreme lateral edge of the hindbrain for which we were not able to identify axons. Based on their somatic locations, we inferred that neurons with only ipsilaterally projecting axons are glutamatergic, whereas neurons with only contralaterally projecting axons are largely GABAergic. Dendritic and synaptic organization of the ipsilaterally projecting neurons suggests a broad sampling from inputs on the ipsilateral side. We also observed the first conclusive evidence of synapses between integrator neurons, which have long been hypothesized by recurrent network models of integration via positive feedback.

Automated computation of arbor densities: a step toward identifying neuronal cell types.

The shape and position of a neuron convey information regarding its molecular and functional identity. The identification of cell types from structure, a classic method, relies on the time-consuming step of arbor tracing. However, as genetic tools and imaging methods make data-driven approaches to neuronal circuit analysis feasible, the need for automated processing increases. Here, we first establish that mouse retinal ganglion cell types can be as precise about distributing their arbor volumes across the inner plexiform layer as they are about distributing the skeletons of the arbors. Then, we describe an automated approach to computing the spatial distribution of the dendritic arbors, or arbor density, with respect to a global depth coordinate based on this observation. Our method involves three-dimensional reconstruction of neuronal arbors by a supervised machine learning algorithm, post-processing of the enhanced stacks to remove somata and isolate the neuron of interest, and registration of neurons to each other using automatically detected arbors of the starburst amacrine interneurons as fiducial markers. In principle, this method could be generalizable to other structures of the CNS, provided that they allow sparse labeling of the cells and contain a reliable axis of spatial reference.

Ring-shaped neuronal networks: a platform to study persistent activity.

Persistent activity in the brain is involved in working memory and motor planning. The ability of the brain to hold information 'online' long after an initiating stimulus is a hallmark of brain areas such as the prefrontal cortex. Recurrent network loops such as the thalamocortical loop and reciprocal loops in the cortex are potential substrates that can support such activity. However, native brain circuitry makes it difficult to study mechanisms underlying such persistent activity. Here we propose a platform to study synaptic mechanisms of such persistent activity by constraining neuronal networks to a recurrent loop like geometry. Using a polymer stamping technique, adhesive proteins are transferred onto glass substrates in a precise ring shape. Primary rat hippocampal cultures were capable of forming ring-shaped networks containing 40-60 neurons. Calcium imaging of these networks show evoked persistent activity in an all-or-none manner. Blocking inhibition with bicuculline methaiodide (BMI) leads to an increase in the duration of persistent activity. These persistent phases were abolished by blockade of asynchronous neurotransmitter release by ethylene glycol tetraacetic acid (EGTA-AM).