Search

A prominent trend in single-cell transcriptomics is providing spatial context alongside a characterization of each cell’s molecular state. This typically requires targeting ana prioriselection of genes, often covering less than 1% of the genome, and a key question is how to optimally determine the small gene panel. Reference data from these methods covering the whole genome is unavailable, and using single-cell RNA sequencing (scRNA-seq) datasets as a surrogate can result in suboptimal gene panels due to the fundamentally different data distributions across technologies. We address these challenges by introducing a flexible deep learning framework, PERSIST, to identify informative gene targets for spatial transcriptomics studies by leveraging existing scRNA-seq data. Using datasets spanning different brain regions, species, and scRNA-seq technologies, we show that PERSIST reliably identifies gene panels that provide more accurate prediction of the genome-wide expression profile, thereby capturing more information with fewer genes. Furthermore, PERSIST can be adapted to meet specific biological goals, such as classifying cell types or discerning neuronal electrical properties. Finally, via a simulation study based on a recentin situhybridization-based dataset, we demonstrate that PERSIST’s binarization of gene expression levels enables models trained on scRNA-seq data to generalize with input data obtained using spatial transcriptomics, despite the complex domain shift between these technologies.

Neuronal anatomy is central to the organization and function of brain cell types. However, anatomical variability within apparently homogeneous populations of cells can obscure such insights. Here, we report large-scale automation of arbor reconstruction on a dataset of 802 inhibitory neurons characterized using the Patch-seq method, which enables measurement of multiple properties from individual neurons, including local morphology and transcriptional signature. We demonstrate that these automated reconstructions can be used in the same manner as manual reconstructions to understand the relationship between many cellular properties used to define cell types. We uncover molecular correlates of laminar innervation on multiple molecularly defined neuronal subclasses and types. In particular, our results reveal molecular correlates of the variability in Layer 1 (L1) innervation even in a transcriptomically defined subpopulation of Martinotti cells in the adult mouse neocortex.

Significance Synaptic connectivity provides the foundation for our present understanding of neuronal network function, but static connectivity cannot explain learning and memory. We propose a computational role for the diversity of cortical neuronal types and their associated cell-type–specific neuromodulators in improving the efficiency of synaptic weight adjustments for task learning in neuronal networks., Brains learn tasks via experience-driven differential adjustment of their myriad individual synaptic connections, but the mechanisms that target appropriate adjustment to particular connections remain deeply enigmatic. While Hebbian synaptic plasticity, synaptic eligibility traces, and top-down feedback signals surely contribute to solving this synaptic credit-assignment problem, alone, they appear to be insufficient. Inspired by new genetic perspectives on neuronal signaling architectures, here, we present a normative theory for synaptic learning, where we predict that neurons communicate their contribution to the learning outcome to nearby neurons via cell-type–specific local neuromodulation. Computational tests suggest that neuron-type diversity and neuron-type–specific local neuromodulation may be critical pieces of the biological credit-assignment puzzle. They also suggest algorithms for improved artificial neural network learning efficiency.

Animals learn and form memories by jointly adjusting the efficacy of their synapses. How they efficiently solve the underlying temporal credit assignment problem remains elusive. Here, we re-analyze the mathematical basis of gradient descent learning in recurrent spiking neural networks (RSNNs) in light of the recent single-cell transcriptomic evidence for cell-type-specific local neuropeptide signaling in the cortex. Our normative theory posits an important role for the notion of neuronal cell types and local diffusive communication by enabling biologically plausible and efficient weight update. While obeying fundamental biological constraints, including separating excitatory vs inhibitory cell types and observing connection sparsity, we trained RSNNs for temporal credit assignment tasks spanning seconds and observed that the inclusion of local modulatory signaling improved learning efficiency. Our learning rule puts forth a novel form of interaction between modulatory signals and synaptic transmission. Moreover, it suggests a computationally efficient learning method for bio-inspired artificial intelligence.

Consistent identification of neurons in different experimental modalities is a key problem in neuroscience. While methods to perform multimodal measurements in the same set of single neurons have become available, parsing complex relationships across different modalities to uncover neuronal identity is a growing challenge. Here, we present an optimization framework to learn coordinated representations of multimodal data, and apply it to a large multimodal dataset profiling mouse cortical interneurons. Our approach reveals strong alignment between transcriptomic and electrophysiological characterizations, enables accurate cross-modal data prediction, and identifies cell types that are consistent across modalities.HighlightsCoupled autoencoders for multimodal assignment, Analysis of Patch-seq data consisting of more than 3000 cells

The balance between excitatory and inhibitory (E and I) synapses is thought to be critical for information processing in neural circuits. However, little is known about the spatial principles of E and I synapses organization across the entire dendritic tree of mammalian neurons. We developed a new, open-source, reconstruction platform for mapping the size and spatial distribution of E and I synapses received by individual, genetically-labeled, layer 2/3 cortical pyramidal neurons (PNs) in vivo. We mapped over 90,000 E and I synapses across twelve L2/3 PNs and uncovered structured organization of E and I synapses across dendritic domains as well as within individual dendritic segments in these neurons. Despite significant, domain-specific, variations in the absolute density of E and I synapses, their ratio is strikingly balanced locally across dendritic segments. Computational modeling indicates that this spatially-precise E/I balance dampens dendritic voltage fluctuations and strongly impacts neuronal firing output.

Neurons are frequently classified into distinct groups or cell types on the basis of structural, physiological, or genetic attributes. To better constrain the definition of neuronal cell types, we characterized the transcriptomes and intrinsic physiological properties of over 3,700 GABAergic mouse visual cortical neurons and reconstructed the local morphologies of 350 of those neurons. We found that most transcriptomic types (t-types) occupy specific laminar positions within mouse visual cortex, and many of those t-types exhibit consistent electrophysiological and morphological features. We observed that these properties could vary continuously between t-types, which limited the ability to predict specific t-types from other data modalities. Despite that, the data support the presence of at least 20 interneuron met-types that have congruent morphological, electrophysiological, and transcriptomic properties.HighlightsPatch-seq data obtained from >3,700 GABAergic cortical interneuronsComprehensive characterization of morpho-electric features of transcriptomic types20 interneuron met-types that have congruent properties across data modalitiesDifferent Sst met-types preferentially innervate different cortical layers

Neurons are frequently classified into distinct groups or cell types on the basis of structural, physiological, or genetic attributes. To better constrain the definition of neuronal cell types, we characterized the transcriptomes and intrinsic physiological properties of over 3,700 GABAergic mouse visual cortical neurons and reconstructed the local morphologies of 350 of those neurons. We found that most transcriptomic types (t-types) occupy specific laminar positions within mouse visual cortex, and many of those t-types exhibit consistent electrophysiological and morphological features. We observed that these properties could vary continuously between t- types, which limited the ability to predict specific t-types from other data modalities. Despite that, the data support the presence of at least 20 interneuron met-types that have congruent morphological, electrophysiological, and transcriptomic properties.

Neuropeptides, members of a large and evolutionarily ancient family of proteinaceous cell-cell signaling molecules, are widely recognized as extremely potent regulators of brain function and behavior. At the cellular level, neuropeptides are known to act mainly via modulation of ion channel and synapse function, but functional impacts emerging at the level of complex cortical synaptic networks have resisted mechanistic analysis. New findings from single-cell RNA-seq transcriptomics now illuminate intricate patterns of cortical neuropeptide signaling gene expression and new tools now offer powerful molecular access to cortical neuropeptide signaling. Here we highlight some of these new findings and tools, focusing especially on prospects for experimental and theoretical exploration of peptidergic and synaptic networks interactions underlying cortical function and plasticity., Comment: 22 pages, 4 figures, 1 table, to be published in Current Opinion in Neurobiology

Neurons are frequently classified into distinct types on the basis of structural, physiological, or genetic attributes. To better constrain the definition of neuronal cell types, we characterized the transcriptomes and intrinsic physiological properties of over 4,200 mouse visual cortical GABAergic interneurons and reconstructed the local morphologies of 517 of those neurons. We find that most transcriptomic types (t-types) occupy specific laminar positions within visual cortex, and, for most types, the cells mapping to a t-type exhibit consistent electrophysiological and morphological properties. These properties display both discrete and continuous variation among t-types. Through multimodal integrated analysis, we define 28 met-types that have congruent morphological, electrophysiological, and transcriptomic properties and robust mutual predictability. We identify layer-specific axon innervation pattern as a defining feature distinguishing different met-types. These met-types represent a unified definition of cortical GABAergic interneuron types, providing a systematic framework to capture existing knowledge and bridge future analyses across different modalities.

This paper studiesmeasurement linkage. An example from cell biology helps explain the problem: imagine for a given cell we can either sequence the cell’s RNA or we can examine its morphology, but not both. Given a cell’s morphology, what do we expect to see in its RNA? Given a cell’s RNA, what do we expect in its morphology? More broadly, given a measurement of one type, can we predict measurements of the other type? This measurement linkage problem arises in many scientific and technological fields. To solve this problem, we develop a nonparametric approach we dub the “Markov link method” (MLM). The MLM makes a conditional independence assumption that holds in many multi-measurement contexts and provides a way to estimate thelink, the conditional probability of one type of measurement given the other. We derive conditions under which the MLM estimator is consistent and we use simulated data to show that it provides accurate measures of uncertainty. We evaluate the MLM on real data generated by a pair of single-cell RNA sequencing techniques. The MLM characterizes the link between them and helps connect the two notions of cell type derived from each technique. Further, the MLM reveals that some aspects of the link cannot be determined from the available data, and suggests new experiments that would allow for better estimates.Significance StatementNovel experimental techniques are developing quickly, and each technique gives new perspectives. Ideally we would build theories that account for many perspectives at once. This is not easy. One challenge is that many experiments use measurement techniques that alter or destroy the subject, making it impossible to measure the same subject with both techniques and difficult to combine data from different experiments. In this paper we develop the Markov Link Method, a new tool that overcomes this challenge.

SUMMARYThe balance between excitatory and inhibitory (E and I) synaptic inputs is thought to be critical for information processing in neural circuits. However, little is known about the principles of spatial organization of E and I synapses across the entire dendritic tree of mammalian neurons. We developed a new, open-source, reconstruction platform for mapping the size and spatial distribution of E and I synapses received by individual, genetically-labeled, layer 2/3 cortical pyramidal neurons (PNs) in vivo. We mapped over 90,000 E and I synapses across twelve L2/3 PNs and uncovered structured organization of E and I synapses across dendritic domains as well as within individual dendritic segments in these cells. Despite significant, domain-specific, variations in the absolute density of E and I synapses, their ratio is strikingly balanced locally across dendritic segments. Computational modeling indicates that this spatially-precise E/I balance dampens dendritic voltage fluctuations and strongly impacts neuronal firing output.

We describe recent progress toward defining neuronal cell types in the mouse retina and attempt to extract lessons that may be generally useful in the mammalian brain. Achieving a comprehensive catalog of retinal cell types now appears within reach, because researchers have achieved consensus concerning two fundamental challenges. The first is accuracy—defining pure cell types rather than settling for neuronal classes that are mixtures of types. The second is completeness—developing methods guaranteed to eventually identify all cell types, as well as criteria for determining when all types have been found. Case studies illustrate how these two challenges are handled by combining state-of-the-art molecular, anatomical, and physiological techniques. Progress is also being made in observing and modeling connectivity between cell types. Scaling up to larger brain regions, such as the cortex, will require not only technical advances but also careful consideration of the challenges of accuracy and completeness.

A practical acceleration algorithm for real-time magnetic resonance imaging (MRI) is presented. Neither separate training scans nor embedded training samples are used. The Kalman filter based algorithm provides a fast and causal reconstruction of dynamic MRI acquisitions with arbitrary readout trajectories. The algorithm is tested against abrupt changes in the imaging conditions and offline reconstructions of in vivo cardiac MRI experiments are presented.

Cataloged from PDF version of article. Periodicity and discreteness are Fourier duals in the same sense as operators such as coordinate multiplication and differentiation, and translation and phase shift. The fractional Fourier transform allows interpolation between such operators which gradually evolve from one member of the dual pair to the other as the fractional order goes from zero to one. Here, we similarly discuss the interpolation between the dual properties of periodicity and discreteness, showing how one evolves into the other as the order goes from zero to one. We also discuss the concepts of partial discreteness and partial periodicity and relate them to fractional discreteness and periodicity.

Nature Communications 5: Article number: 3512 (2014); Published: 24 March 2014; Updated: 18 September 2014. There were errors associated with the Supplementary Data in the version of this Article originally published, which were introduced while the HTML was being prepared. The titles of Supplementary Data 8, 9, 10, 11, 12 and 13 were incorrectly linked to Supplementary Data 7, 8, 9, 10, 11 and 12, respectively, and the file for Supplementary Data 13 was omitted.

The importance of cell types in understanding brain function is widely appreciated but only a tiny fraction of neuronal diversity has been catalogued. Here, we exploit recent progress in genetic definition of cell types in an objective structural approach to neuronal classification. The approach is based on highly accurate quantification of dendritic arbor position relative to neurites of other cells. We test the method on a population of 363 mouse retinal ganglion cells. For each cell, we determine the spatial distribution of the dendritic arbors, or “arbor density” with reference to arbors of an abundant, well-defined interneuronal type. The arbor densities are sorted into a number of clusters that is set by comparison with several molecularly defined cell types. The algorithm reproduces the genetic classes that are pure types, and detects six newly clustered cell types that await genetic definition.

Chemotopic odour representations in the olfactory bulb are transferred to multiple forebrain areas and translated into appropriate output responses. However, a comprehensive projection map of bulbar output neurons at single-axon resolution is lacking in vertebrates. Here we unravel a projectome of the zebrafish olfactory bulb through genetic single-neuron tracing and image registration. We show that five major target regions receive distinct modes of projections from olfactory bulb glomeruli. The central portion of posterior telencephalon receives non-selective, interspersed inputs from all glomeruli, whereas the ventral telencephalon is diffusely innervated by axons from particular glomerular clusters. The right habenula and posterior tuberculum (diencephalic nuclei) receive convergent inputs from restricted and all glomerular clusters, respectively. The bulbar recurrent projections are coarsely topographic. Thus, the primary chemotopic organization is transformed into distinct sensory representations in higher olfactory centres. These findings provide a framework to understand general principles as well as species-specific features in decoding of odour information.

Time series of in vivo magnetic resonance images exhibit high levels of temporal correlation. Higher temporal resolution reconstructions are obtained by acquiring data at a fraction of the Nyquist rate and resolving the resulting aliasing using the correlation information. The dynamic imaging experiment is modeled as a linear dynamical system. A Kalman filter based unaliasing reconstruction is described for accelerated dynamic magnetic resonance imaging (MRI). The algorithm handles arbitrary readout trajectories naturally. The reconstruction is causal and very fast, making it applicable to real-time imaging. In vivo results are presented for cardiac MRI of healthy volunteers.

Continuum extensions of common dual pairs of operators are presented and consolidated, based on the fractional Fourier transform. In particular, the fractional chirp multiplication, fractional chirp convolution, and fractional scaling operators are defined and expressed in terms of their common nonfractional special cases, revealing precisely how they are interpolations of their conventional counterparts. Optical realizations of these operators are possible with use of common physical components. These three operators can be interpreted as fractional lenses, fractional free space, and fractional imaging systems, respectively. Any optical system consisting of an arbitrary concatenation of sections of free space and thin lenses can be interpreted as a fractional imaging system with spherical reference surfaces. As a special case, a system departing from the classical single-lens imaging condition can be interpreted as a fractional imaging system. © 2003 Optical Society of America.

Identifying the cell types constituting brain circuits is a fundamental question in neuroscience and motivates the generation of taxonomies based on electrophysiological, morphological and molecular single cell properties. Establishing the correspondence across data modalities and understanding the underlying principles has proven challenging. Bio-realistic computational models offer the ability to probe cause-and-effect and have historically been used to explore phenomena at the single-neuron level. Here we introduce a computational optimization workflow used for the generation and evaluation of more than 130 million single neuron models with active conductances. These models were based on 230 in vitro electrophysiological experiments followed by morphological reconstruction from the mouse visual cortex. We show that distinct ion channel conductance vectors exist that distinguish between major cortical classes with passive and h-channel conductances emerging as particularly important for classification. Next, using models of genetically defined classes, we show that differences in specific conductances predicted from the models reflect differences in gene expression in excitatory and inhibitory cell types as experimentally validated by single-cell RNA-sequencing. The differences in these conductances, in turn, explain many of the electrophysiological differences observed between cell types. Finally, we show the robustness of the herein generated single-cell models as representations and realizations of specific cell types in face of biological variability and optimization complexity. Our computational effort generated models that reconcile major single-cell data modalities that define cell types allowing for causal relationships to be examined.HighlightsGeneration and evaluation of more than 130 million single-cell models with active conductances along the reconstructed morphology faithfully recapitulate the electrophysiology of 230 in vitro experiments.Optimized ion channel conductances along the cellular morphology (‘all-active’) are characteristic of model complexity and offer enhanced biophysical realism.Ion channel conductance vectors of all-active models classify transcriptomically defined cell-types.Cell type differences in ion channel conductances predicted by the models correlate with experimentally measured single-cell gene expression differences in inhibitory (Pvalb, Sst, Htr3a) and excitatory (Nr5a1, Rbp4) classes.A set of ion channel conductances identified by comparing between cell type model populations explain electrophysiology differences between these types in simulations and brain slice experiments.All-active models recapitulate multimodal properties of excitatory and inhibitory cell types offering a systematic and causal way of linking differences between them.