A Systematic Evaluation of Interneuron Morphology Representations for Cell Type Discrimination

NeoCoMM: A neocortical neuroinspired computational model for the reconstruction and simulation of epileptiform events

BACKGROUND

Understanding the pathophysiological dynamics that underline Interictal Epileptiform Events (IEEs) such as epileptic spikes, spike-and-waves or High-Frequency Oscillations (HFOs) is of major importance in the context of neocortical refractory epilepsy, as it paves the way for the development of novel therapies. Typically, these events are detected in Local Field Potential (LFP) recordings obtained through depth electrodes during pre-surgical investigations. Although essential, the underlying pathophysiological mechanisms for the generation of these epileptic neuromarkers remain unclear. The aim of this paper is to propose a novel neurophysiologically relevant reconstruction of the neocortical microcircuitry in the context of epilepsy. This reconstruction intends to facilitate the analysis of a comprehensive set of parameters encompassing physiological, morphological, and biophysical aspects that directly impact the generation and recording of different IEEs.

METHOD

a novel microscale computational model of an epileptic neocortical column was introduced. This model incorporates the intricate multilayered structure of the cortex and allows for the simulation of realistic interictal epileptic signals. The proposed model was validated through comparisons with real IEEs recorded using intracranial stereo-electroencephalography (SEEG) signals from both humans and animals. Using the model, the user can recreate epileptiform patterns observed in different species (human, rodent, and mouse) and study the intracellular activity associated with these patterns.

RESULTS

Our model allowed us to unravel the relationship between glutamatergic and GABAergic synaptic transmission of the epileptic neural network and the type of generated IEE. Moreover, sensitivity analyses allowed for the exploration of the pathophysiological parameters responsible for the transitions between these events. Finally, the presented modeling framework also provides an Electrode Tissue Model (ETI) that adds realism to the simulated signals and offers the possibility of studying their sensitivity to the electrode characteristics.

CONCLUSION

The model (NeoCoMM) presented in this work can be of great use in different applications since it offers an in silico framework for sensitivity analysis and hypothesis testing. It can also be used as a starting point for more complex studies.

High-throughput analysis of dendrite and axonal arbors reveals transcriptomic correlates of neuroanatomy

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 neuronal morphology reconstruction and analysis on a dataset of 813 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 some, but not all, cellular properties used to define cell types. We uncover gene expression correlates of laminar innervation on multiple transcriptomically defined neuronal subclasses and types. In particular, our results reveal correlates of the variability in Layer 1 (L1) axonal innervation in a transcriptomically defined subpopulation of Martinotti cells in the adult mouse neocortex.

Loss or gain of function? Effects of ion channel mutations on neuronal firing depend on the neuron type

Introduction

Clinically relevant mutations to voltage-gated ion channels, called channelopathies, alter ion channel function, properties of ionic currents, and neuronal firing. The effects of ion channel mutations are routinely assessed and characterized as loss of function (LOF) or gain of function (GOF) at the level of ionic currents. However, emerging personalized medicine approaches based on LOF/GOF characterization have limited therapeutic success. Potential reasons are among others that the translation from this binary characterization to neuronal firing is currently not well-understood-especially when considering different neuronal cell types. In this study, we investigate the impact of neuronal cell type on the firing outcome of ion channel mutations.

Methods

To this end, we simulated a diverse collection of single-compartment, conductance-based neuron models that differed in their composition of ionic currents. We systematically analyzed the effects of changes in ion current properties on firing in different neuronal types. Additionally, we simulated the effects of known mutations in KCNA1 gene encoding the K_V1.1 potassium channel subtype associated with episodic ataxia type 1 (EA1).

Results

These simulations revealed that the outcome of a given change in ion channel properties on neuronal excitability depends on neuron type, i.e., the properties and expression levels of the unaffected ionic currents.

Discussion

Consequently, neuron-type specific effects are vital to a full understanding of the effects of channelopathies on neuronal excitability and are an important step toward improving the efficacy and precision of personalized medicine approaches.

Single-Cell and Subcellular Analysis Using Ultrahigh Resolution 21 T MALDI FTICR Mass Spectrometry

The mammalian brain contains ∼20,000 distinct lipid species that contribute to its structural organization and function. The lipid profiles of cells change in response to a variety of cellular signals and environmental conditions that result in modulation of cell function through alteration of phenotype. The limited sample material combined with the vast chemical diversity of lipids makes comprehensive lipid profiling of individual cells challenging. Here, we leverage the resolving power of a 21 T Fourier-transform ion cyclotron resonance (FTICR) mass spectrometer for chemical characterization of individual hippocampal cells at ultrahigh mass resolution. The accuracy of the acquired data allowed differentiation of freshly isolated and cultured hippocampal cell populations, as well as finding differences in lipids between the soma and neuronal processes of the same cell. Differences in lipids include TG 42:2 observed solely in the cell bodies and SM 34:1;O2 found only in the cellular processes. The work represents the first mammalian single cells analyzed at ultrahigh resolution and is an advance in the performance of mass spectrometry (MS) for single-cell research.

Semi-Automated Quantitative Evaluation of Neuron Developmental Morphology In Vitro Using the Change-Point Test

Neuron morphology gives rise to distinct axons and dendrites and plays an essential role in neuronal functionality and circuit dynamics. In rat hippocampal neurons, morphological development occurs over roughly one week in vitro. This development has been qualitatively described as occurring in 5 stages. Still, there is a need to quantify cell growth to monitor cell culture health, understand cell responses to sensory cues, and compare experimental results and computational growth model predictions. To address this need, embryonic rat hippocampal neurons were observed in vitro over six days, and their processes were quantified using both standard morphometrics (degree, number of neurites, total length, and tortuosity) and new metrics (distance between change points, relative turning angle, and the number of change points) based on the Change-Point Test to track changes in path trajectories. Of the standard morphometrics, the total length of neurites per cell and the number of endpoints were significantly different between 0.5, 1.5, and 4 days in vitro, which are typically associated with Stages 2-4. Using the Change-Point Test, the number of change points and the average distance between change points per cell were also significantly different between those key time points. This work highlights key quantitative characteristics, both among common and novel morphometrics, that can describe neuron development in vitro and provides a foundation for analyzing directional changes in neurite growth for future studies.

Polymer Physics-Based Classification of Neurons

Recognizing that diverse morphologies of neurons are reminiscent of structures of branched polymers, we put forward a principled and systematic way of classifying neurons that employs the ideas of polymer physics. In particular, we use 3D coordinates of individual neurons, which are accessible in recent neuron reconstruction datasets from electron microscope images. We numerically calculate the form factor, F(q), a Fourier transform of the distance distribution of particles comprising an object of interest, which is routinely measured in scattering experiments to quantitatively characterize the structure of materials. For a polymer-like object consisting of n monomers spanning over a length scale of r, F(q) scales with the wavenumber [Formula: see text] as [Formula: see text] at an intermediate range of q, where [Formula: see text] is the fractal dimension or the inverse scaling exponent ([Formula: see text]) characterizing the geometrical feature ([Formula: see text]) of the object. F(q) can be used to describe a neuron morphology in terms of its size ([Formula: see text]) and the extent of branching quantified by [Formula: see text]. By defining the distance between F(q)s as a measure of similarity between two neuronal morphologies, we tackle the neuron classification problem. In comparison with other existing classification methods for neuronal morphologies, our F(q)-based classification rests solely on 3D coordinates of neurons with no prior knowledge of morphological features. When applied to publicly available neuron datasets from three different organisms, our method not only complements other methods but also offers a physical picture of how the dendritic and axonal branches of an individual neuron fill the space of dense neural networks inside the brain.

Large scale similarity search across digital reconstructions of neural morphology

Most functions of the nervous system depend on neuronal and glial morphology. Continuous advances in microscopic imaging and tracing software have provided an increasingly abundant availability of 3D reconstructions of arborizing dendrites, axons, and processes, allowing their detailed study. However, efficient, large-scale methods to rank neural morphologies by similarity to an archetype are still lacking. Using the NeuroMorpho.Org database, we present a similarity search software enabling fast morphological comparison of hundreds of thousands of neural reconstructions from any species, brain regions, cell types, and preparation protocols. We compared the performance of different morphological measurements: 1) summary morphometrics calculated by L-Measure, 2) persistence vectors, a vectorized descriptor of branching structure, 3) the combination of the two. In all cases, we also investigated the impact of applying dimensionality reduction using principal component analysis (PCA). We assessed qualitative performance by gauging the ability to rank neurons in order of visual similarity. Moreover, we quantified information content by examining explained variance and benchmarked the ability to identify occasional duplicate reconstructions of the same specimen. We also compared two different methods for selecting the number of principal components using this benchmark. The results indicate that combining summary morphometrics and persistence vectors with applied PCA using maximum likelihood based automatic dimensionality selection provides an information rich characterization that enables efficient and precise comparison of neural morphology. We have deployed the similarity search as open-source online software both through a user-friendly graphical interface and as an API for programmatic access.

nGauge: Integrated and Extensible Neuron Morphology Analysis in Python

The study of neuron morphology requires robust and comprehensive methods to quantify the differences between neurons of different subtypes and animal species. Several software packages have been developed for the analysis of neuron tracing results stored in the standard SWC format. The packages, however, provide relatively simple quantifications and their non-extendable architecture prohibit their use for advanced data analysis and visualization. We developed nGauge, a Python toolkit to support the parsing and analysis of neuron morphology data. As an application programming interface (API), nGauge can be referenced by other popular open-source software to create custom informatics analysis pipelines and advanced visualizations. nGauge defines an extendable data structure that handles volumetric constructions (e.g. soma), in addition to the SWC linear reconstructions, while remaining lightweight. This greatly extends nGauge's data compatibility.

Topological Sholl descriptors for neuronal clustering and classification

Neuronal morphology is a fundamental factor influencing information processing within neurons and networks. Dendritic morphology in particular can widely vary among cell classes, brain regions, and animal species. Thus, accurate quantitative descriptions allowing classification of large sets of neurons is essential for their structural and functional characterization. Current robust and unbiased computational methods that characterize groups of neurons are scarce. In this work, we introduce a novel technique to study dendritic morphology, complementing and advancing many of the existing techniques. Our approach is to conceptualize the notion of a Sholl descriptor and associate, for each morphological feature, and to each neuron, a function of the radial distance from the soma, taking values in a metric space. Functional distances give rise to pseudo-metrics on sets of neurons which are then used to perform the two distinct tasks of clustering and classification. To illustrate the use of Sholl descriptors, four datasets were retrieved from the large public repository https://neuromorpho.org/ comprising neuronal reconstructions from different species and brain regions. Sholl descriptors were subsequently computed, and standard clustering methods enhanced with detection and metric learning algorithms were then used to objectively cluster and classify each dataset. Importantly, our descriptors outperformed conventional morphometric techniques (L-Measure metrics) in several of the tested datasets. Therefore, we offer a novel and effective approach to the analysis of diverse neuronal cell types, and provide a toolkit for researchers to cluster and classify neurons.

Quantifying neuronal structural changes over time using dynamic morphometrics

Brain circuit development involves tremendous structural formation and rearrangement of dendrites, axons, and the synaptic connections between them. Direct studies of neuronal morphogenesis are now possible through recent developments in multiple technologies, including single-neuron labeling, time-lapse imaging in intact tissues, and 4D rendering software capable of tracking neural growth over periods spanning minutes to days. These methods allow detailed quantification of structural changes of neurons over time, called dynamic morphometrics, providing new insights into fundamental growth patterns, underlying molecular mechanisms, and the intertwined influences of external factors, including neural activity, and intrinsic genetic programs. Here, we review the methods of dynamic morphometrics sampling and analyses, focusing on their applications to studies of activity-driven dendritogenesis in vertebrate systems.

Phenotypic variation of transcriptomic cell types in mouse motor cortex

Cortical neurons exhibit extreme diversity in gene expression as well as in morphological and electrophysiological properties1,2. Most existing neural taxonomies are based on either transcriptomic3,4 or morpho-electric5,6 criteria, as it has been technically challenging to study both aspects of neuronal diversity in the same set of cells7. Here we used Patch-seq8 to combine patch-clamp recording, biocytin staining, and single-cell RNA sequencing of more than 1,300 neurons in adult mouse primary motor cortex, providing a morpho-electric annotation of almost all transcriptomically defined neural cell types. We found that, although broad families of transcriptomic types (those expressing Vip, Pvalb, Sst and so on) had distinct and essentially non-overlapping morpho-electric phenotypes, individual transcriptomic types within the same family were not well separated in the morpho-electric space. Instead, there was a continuum of variability in morphology and electrophysiology, with neighbouring transcriptomic cell types showing similar morpho-electric features, often without clear boundaries between them. Our results suggest that neuronal types in the neocortex do not always form discrete entities. Instead, neurons form a hierarchy that consists of distinct non-overlapping branches at the level of families, but can form continuous and correlated transcriptomic and morpho-electrical landscapes within families.