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

Synaptic architecture of leg and wing premotor control networks in Drosophila

Animal movement is controlled by motor neurons (MNs), which project out of the central nervous system to activate muscles1. MN activity is coordinated by complex premotor networks that facilitate the contribution of individual muscles to many different behaviours2-6. Here we use connectomics7 to analyse the wiring logic of premotor circuits controlling the Drosophila leg and wing. We find that both premotor networks cluster into modules that link MNs innervating muscles with related functions. Within most leg motor modules, the synaptic weights of each premotor neuron are proportional to the size of their target MNs, establishing a circuit basis for hierarchical MN recruitment. By contrast, wing premotor networks lack proportional synaptic connectivity, which may enable more flexible recruitment of wing steering muscles. Through comparison of the architecture of distinct motor control systems within the same animal, we identify common principles of premotor network organization and specializations that reflect the unique biomechanical constraints and evolutionary origins of leg and wing motor control.

Synaptic architecture of leg and wing premotor control networks in Drosophila

Animal movement is controlled by motor neurons (MNs), which project out of the central nervous system to activate muscles. MN activity is coordinated by complex premotor networks that allow individual muscles to contribute to many different behaviors. Here, we use connectomics to analyze the wiring logic of premotor circuits controlling the Drosophila leg and wing. We find that both premotor networks cluster into modules that link MNs innervating muscles with related functions. Within most leg motor modules, the synaptic weights of each premotor neuron are proportional to the size of their target MNs, establishing a circuit basis for hierarchical MN recruitment. In contrast, wing premotor networks lack proportional synaptic connectivity, which may allow wing steering muscles to be recruited with different relative timing. By comparing the architecture of distinct limb motor control systems within the same animal, we identify common principles of premotor network organization and specializations that reflect the unique biomechanical constraints and evolutionary origins of leg and wing motor control.

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 to automate a variety of tasks that would otherwise require extensive manual effort, such as state of the art automated post-hoc proofreading of merge errors, cell classification, spine detection, axon-dendritic proximities, and computation of other features. These features enable many downstream analyses of neural morphology and connectivity, making these new massive and complex datasets more accessible to neuroscience researchers focused on a variety of scientific questions.

Multidimensional analysis of cortical interneuron synaptic features reveals underlying synaptic heterogeneity

Cortical interneurons represent a diverse set of neuronal subtypes characterized in part by their striking degree of synaptic specificity. However, little is known about the extent of synaptic diversity because of the lack of unbiased methods to extract synaptic features among interneuron subtypes. Here, we develop an approach to aggregate image features from fluorescent confocal images of interneuron synapses and their post-synaptic targets, in order to characterize the heterogeneity of synapses at fine scale. We started by training a model that recognizes pre- and post-synaptic compartments and then determines the target of each genetically-identified interneuron synapse in vitro and in vivo. Our model extracts hundreds of spatial and intensity features from each analyzed synapse, constructing a multidimensional data set, consisting of millions of synapses, which allowed us to perform an unsupervised analysis on this dataset, uncovering novel synaptic subgroups. The subgroups were spatially distributed in a highly structured manner that revealed the local underlying topology of the postsynaptic environment. Dendrite-targeting subgroups were clustered onto subdomains of the dendrite along the proximal to distal axis. Soma-targeting subgroups were enriched onto different postsynaptic cell types. We also find that the two main subclasses of interneurons, basket cells and somatostatin interneurons, utilize distinct strategies to enact inhibitory coverage. Thus, our analysis of multidimensional synaptic features establishes a conceptual framework for studying interneuron synaptic diversity.

Neuronal "parts list" and wiring diagram for a visual system

A catalog of neuronal cell types has often been called a "parts list" of the brain, and regarded as a prerequisite for understanding brain function. In the optic lobe of Drosophila, rules of connectivity between cell types have already proven essential for understanding fly vision. Here we analyze the fly connectome to complete the list of cell types intrinsic to the optic lobe, as well as the rules governing their connectivity. We more than double the list of known types. Most new cell types contain between 10 and 100 cells, and integrate information over medium distances in the visual field. Some existing type families (Tm, Li, and LPi) at least double in number of types. We introduce a new Sm interneuron family, which contains more types than any other, and three new families of cross-neuropil types. Self-consistency of cell types is demonstrated through automatic assignment of cells to types by distance in high-dimensional feature space, and further validation is provided by algorithms that select small subsets of discriminative features. Cell types with similar connectivity patterns divide into clusters that are interpretable in terms of motion, object, and color vision. Our work showcases the advantages of connectomic cell typing: complete and unbiased sampling, a rich array of features based on connectivity, and reduction of the connectome to a drastically simpler wiring diagram of cell types, with immediate relevance for brain function and development.

Emergence of co-tuning in inhibitory neurons as a network phenomenon mediated by randomness, correlations, and homeostatic plasticity

Cortical excitatory neurons show clear tuning to stimulus features, but the tuning properties of inhibitory interneurons are ambiguous. While inhibitory neurons have been considered to be largely untuned, some studies show that some parvalbumin-expressing (PV) neurons do show feature selectivity and participate in co-tuned subnetworks with pyramidal neurons. In this study, we first use mean-field theory to demonstrate that a combination of homeostatic plasticity governing the synaptic dynamics of the connections from PV to excitatory neurons, heterogeneity in the excitatory postsynaptic potentials that impinge on PV neurons, and shared correlated input from layer 4 results in the functional and structural self-organization of PV subnetworks. Second, we show that structural and functional feature tuning of PV neurons emerges more clearly at the network level, i.e., that population-level measures identify functional and structural co-tuning of PV neurons that are not evident in pairwise individual-level measures. Finally, we show that such co-tuning can enhance network stability at the cost of reduced feature selectivity.

Rabies virus-based barcoded neuroanatomy resolved by single-cell RNA and in situ sequencing

Mapping the connectivity of diverse neuronal types provides the foundation for understanding the structure and function of neural circuits. High-throughput and low-cost neuroanatomical techniques based on RNA barcode sequencing have the potential to map circuits at cellular resolution and a brain-wide scale, but existing Sindbis virus-based techniques can only map long-range projections using anterograde tracing approaches. Rabies virus can complement anterograde tracing approaches by enabling either retrograde labeling of projection neurons or monosynaptic tracing of direct inputs to genetically targeted postsynaptic neurons. However, barcoded rabies virus has so far been only used to map non-neuronal cellular interactions in vivo and synaptic connectivity of cultured neurons. Here we combine barcoded rabies virus with single-cell and in situ sequencing to perform retrograde labeling and transsynaptic labeling in the mouse brain. We sequenced 96 retrogradely labeled cells and 295 transsynaptically labeled cells using single-cell RNA-seq, and 4130 retrogradely labeled cells and 2914 transsynaptically labeled cells in situ. We found that the transcriptomic identities of rabies virus-infected cells can be robustly identified using both single-cell RNA-seq and in situ sequencing. By associating gene expression with connectivity inferred from barcode sequencing, we distinguished long-range projecting cortical cell types from multiple cortical areas and identified cell types with converging or diverging synaptic connectivity. Combining in situ sequencing with barcoded rabies virus complements existing sequencing-based neuroanatomical techniques and provides a potential path for mapping synaptic connectivity of neuronal types at scale.

Postsynaptic cell type and synaptic distance do not determine efficiency of monosynaptic rabies virus spread measured at synaptic resolution

Retrograde monosynaptic tracing using glycoprotein-deleted rabies virus is an important component of the toolkit for investigation of neural circuit structure and connectivity. It allows for the identification of first-order presynaptic connections to cell populations of interest across both the central and peripheral nervous system, helping to decipher the complex connectivity patterns of neural networks that give rise to brain function. Despite its utility, the factors that influence the probability of transsynaptic rabies spread are not well understood. While it is well established that expression levels of rabies glycoprotein used to trans-complement G-deleted rabies can result in large changes in numbers of inputs labeled per starter cell (convergence index [CI]), it is not known how typical values of CI relate to the proportions of synaptic contacts or input neurons labeled. And it is not known whether inputs to different cell types, or synaptic contacts that are more proximal or distal to the cell body, are labeled with different probabilities. Here, we use a new rabies virus construct that allows for the simultaneous labeling of pre- and postsynaptic specializations to quantify the proportion of synaptic contacts labeled in mouse primary visual cortex. We demonstrate that with typical conditions about 40% of first-order presynaptic excitatory synapses to cortical excitatory and inhibitory neurons are labeled. We show that using matched tracing conditions there are similar proportions of labeled contacts onto L4 excitatory pyramidal, somatostatin (Sst) inhibitory, and vasoactive intestinal peptide (Vip) starter cell types. Furthermore, we find no difference in the proportions of labeled excitatory contacts onto postsynaptic sites at different subcellular locations.

Connecting single-cell transcriptomes to projectomes in mouse visual cortex

The mammalian brain is composed of diverse neuron types that play different functional roles. Recent single-cell RNA sequencing approaches have led to a whole brain taxonomy of transcriptomically-defined cell types, yet cell type definitions that include multiple cellular properties can offer additional insights into a neuron's role in brain circuits. While the Patch-seq method can investigate how transcriptomic properties relate to the local morphological and electrophysiological properties of cell types, linking transcriptomic identities to long-range projections is a major unresolved challenge. To address this, we collected coordinated Patch-seq and whole brain morphology data sets of excitatory neurons in mouse visual cortex. From the Patch-seq data, we defined 16 integrated morpho-electric-transcriptomic (MET)-types; in parallel, we reconstructed the complete morphologies of 300 neurons. We unified the two data sets with a multi-step classifier, to integrate cell type assignments and interrogate cross-modality relationships. We find that transcriptomic variations within and across MET-types correspond with morphological and electrophysiological phenotypes. In addition, this variation, along with the anatomical location of the cell, can be used to predict the projection targets of individual neurons. We also shed new light on infragranular cell types and circuits, including cell-type-specific, interhemispheric projections. With this approach, we establish a comprehensive, integrated taxonomy of excitatory neuron types in mouse visual cortex and create a system for integrated, high-dimensional cell type classification that can be extended to the whole brain and potentially across species.

The elusive brain perivascular fibroblast: a potential role in vascular stability and homeostasis

In the brain, perivascular fibroblasts (PVFs) reside within the perivascular spaces (PVSs) of arterioles and large venules, however their physiological and pathophysiological roles remain largely unknown. PVFs express numerous extracellular matrix proteins that are found in the basement membrane and PVS surrounding large diameter vessels. PVFs are sandwiched between the mural cell layer and astrocytic endfeet, where they are poised to interact with mural cells, perivascular macrophages, and astrocytes. We draw connections between the more well-studied PVF pro-fibrotic response in ischemic injury and the less understood thickening of the vascular wall and enlargement of the PVS described in dementia and neurodegenerative diseases. We postulate that PVFs may be responsible for stability and homeostasis of the brain vasculature, and may also contribute to changes within the PVS during disease.

Rabies virus-based barcoded neuroanatomy resolved by single-cell RNA and in situ sequencing

Mapping the connectivity of diverse neuronal types provides the foundation for understanding the structure and function of neural circuits. High-throughput and low-cost neuroanatomical techniques based on RNA barcode sequencing have the potential to map circuits at cellular resolution and a brain-wide scale, but existing Sindbis virus-based techniques can only map long-range projections using anterograde tracing approaches. Rabies virus can complement anterograde tracing approaches by enabling either retrograde labeling of projection neurons or monosynaptic tracing of direct inputs to genetically targeted postsynaptic neurons. However, barcoded rabies virus has so far been only used to map non-neuronal cellular interactions in vivo and synaptic connectivity of cultured neurons. Here we combine barcoded rabies virus with single-cell and in situ sequencing to perform retrograde labeling and transsynaptic labeling in the mouse brain. We sequenced 96 retrogradely labeled cells and 295 transsynaptically labeled cells using single-cell RNA-seq, and 4,130 retrogradely labeled cells and 2,914 transsynaptically labeled cells in situ. We found that the transcriptomic identities of rabies virus-infected cells can be robustly identified using both single-cell RNA-seq and in situ sequencing. By associating gene expression with connectivity inferred from barcode sequencing, we distinguished long-range projecting cortical cell types from multiple cortical areas and identified cell types with converging or diverging synaptic connectivity. Combining in situ sequencing with barcoded rabies virus complements existing sequencing-based neuroanatomical techniques and provides a potential path for mapping synaptic connectivity of neuronal types at scale.

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.

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.

Integrating EM and Patch-seq data: Synaptic connectivity and target specificity of predicted Sst transcriptomic types

Neural circuit function is shaped both by the cell types that comprise the circuit and the connections between those cell types ¹ . Neural cell types have previously been defined by morphology ^{2, 3} , electrophysiology ^{4, 5} , transcriptomic expression ^6-8 , connectivity ^9-13 , or even a combination of such modalities ^14-16 . More recently, the Patch-seq technique has enabled the characterization of morphology (M), electrophysiology (E), and transcriptomic (T) properties from individual cells ^17-20 . Using this technique, these properties were integrated to define 28, inhibitory multimodal, MET-types in mouse primary visual cortex ²¹ . It is unknown how these MET-types connect within the broader cortical circuitry however. Here we show that we can predict the MET-type identity of inhibitory cells within a large-scale electron microscopy (EM) dataset and these MET-types have distinct ultrastructural features and synapse connectivity patterns. We found that EM Martinotti cells, a well defined morphological cell type ^{22, 23} known to be Somatostatin positive (Sst+) ^{24, 25} , were successfully predicted to belong to Sst+ MET-types. Each identified MET-type had distinct axon myelination patterns and synapsed onto specific excitatory targets. Our results demonstrate that morphological features can be used to link cell type identities across imaging modalities, which enables further comparison of connectivity in relation to transcriptomic or electrophysiological properties. Furthermore, our results show that MET-types have distinct connectivity patterns, supporting the use of MET-types and connectivity to meaningfully define cell types.