Single‐neuron axonal reconstruction: The search for a wiring diagram of the brain

Microtubule-binding protein MAP1B regulates interstitial axon branching of cortical neurons via the tubulin tyrosination cycle

Regulation of directed axon guidance and branching during development is essential for the generation of neuronal networks. However, the molecular mechanisms that underlie interstitial (or collateral) axon branching in the mammalian brain remain unresolved. Here, we investigate interstitial axon branching in vivo using an approach for precise labeling of layer 2/3 callosal projection neurons (CPNs). This method allows for quantitative analysis of axonal morphology at high acuity and also manipulation of gene expression in well-defined temporal windows. We find that the GSK3β serine/threonine kinase promotes interstitial axon branching in layer 2/3 CPNs by releasing MAP1B-mediated inhibition of axon branching. Further, we find that the tubulin tyrosination cycle is a key downstream component of GSK3β/MAP1B signaling. These data suggest a cell-autonomous molecular regulation of cortical neuron axon morphology, in which GSK3β can release a MAP1B-mediated brake on interstitial axon branching upstream of the posttranslational tubulin code.

Whole-brain Optical Imaging: A Powerful Tool for Precise Brain Mapping at the Mesoscopic Level

The mammalian brain is a highly complex network that consists of millions to billions of densely-interconnected neurons. Precise dissection of neural circuits at the mesoscopic level can provide important structural information for understanding the brain. Optical approaches can achieve submicron lateral resolution and achieve "optical sectioning" by a variety of means, which has the natural advantage of allowing the observation of neural circuits at the mesoscopic level. Automated whole-brain optical imaging methods based on tissue clearing or histological sectioning surpass the limitation of optical imaging depth in biological tissues and can provide delicate structural information in a large volume of tissues. Combined with various fluorescent labeling techniques, whole-brain optical imaging methods have shown great potential in the brain-wide quantitative profiling of cells, circuits, and blood vessels. In this review, we summarize the principles and implementations of various whole-brain optical imaging methods and provide some concepts regarding their future development.

Multi-scale model of axonal and dendritic polarization by transcranial direct current stimulation in realistic head geometry

BACKGROUND

Transcranial direct current stimulation (tDCS) is a non-invasive brain stimulation modality that can alter cortical excitability. However, it remains unclear how the subcellular elements of different neuron types are polarized by specific electric field (E-field) distributions.

OBJECTIVE

To quantify neuronal polarization generated by tDCS in a multi-scale computational model.

METHODS

We embedded layer-specific, morphologically-realistic cortical neuron models in a finite element model of the E-field in a human head and simulated steady-state polarization generated by conventional primary-motor-cortex-supraorbital (M1-SO) and 4 × 1 high-definition (HD) tDCS. We quantified somatic, axonal, and dendritic polarization of excitatory pyramidal cells in layers 2/3, 5, and 6, as well as inhibitory interneurons in layers 1 and 4 of the hand knob.

RESULTS

Axonal and dendritic terminals were polarized more than the soma in all neurons, with peak axonal and dendritic polarization of 0.92 mV and 0.21 mV, respectively, compared to peak somatic polarization of 0.07 mV for 1.8 mA M1-SO stimulation. Both montages generated regions of depolarization and hyperpolarization beneath the M1 anode; M1-SO produced slightly stronger, more diffuse polarization peaking in the central sulcus, while 4 × 1 HD produced higher peak polarization in the gyral crown. The E-field component normal to the cortical surface correlated strongly with pyramidal neuron somatic polarization (R2>0.9), but exhibited weaker correlations with peak pyramidal axonal and dendritic polarization (R2:0.5-0.9) and peak polarization in all subcellular regions of interneurons (R2:0.3-0.6). Simulating polarization by uniform local E-field extracted at the soma approximated the spatial distribution of tDCS polarization but produced large errors in some regions (median absolute percent error: 7.9 %).

CONCLUSIONS

Polarization of pre- and postsynaptic compartments of excitatory and inhibitory cortical neurons may play a significant role in tDCS neuromodulation. These effects cannot be predicted from the E-field distribution alone but rather require calculation of the neuronal response.

MAP1B Regulates Cortical Neuron Interstitial Axon Branching Through the Tubulin Tyrosination Cycle

Regulation of directed axon guidance and branching during development is essential for the generation of neuronal networks. However, the molecular mechanisms that underlie interstitial axon branching in the mammalian brain remain unresolved. Here, we investigate interstitial axon branching in vivo using an approach for precise labeling of layer 2/3 callosal projection neurons (CPNs), allowing for quantitative analysis of axonal morphology at high acuity and also manipulation of gene expression in well-defined temporal windows. We find that the GSK3β serine/threonine kinase promotes interstitial axon branching in layer 2/3 CPNs by releasing MAP1B-mediated inhibition of axon branching. Further, we find that the tubulin tyrosination cycle is a key downstream component of GSK3β/MAP1B signaling. We propose that MAP1B functions as a brake on axon branching that can be released by GSK3β activation, regulating the tubulin code and thereby playing an integral role in sculpting cortical neuron axon morphology.

Multi-scale model of axonal and dendritic polarization by transcranial direct current stimulation in realistic head geometry

Background

Transcranial direct current stimulation (tDCS) is a non-invasive brain stimulation modality that can alter cortical excitability. However, it remains unclear how the subcellular elements of different neuron types are polarized by specific electric field (E-field) distributions.

Objective

To quantify neuronal polarization generated by tDCS in a multi-scale computational model.

Methods

We embedded layer-specific, morphologically-realistic cortical neuron models in a finite element model of the E-field in a human head and simulated steady-state polarization generated by conventional primary-motor-cortex-supraorbital (M1-SO) and 4×1 high-definition (HD) tDCS. We quantified somatic, axonal, and dendritic polarization of excitatory pyramidal cells in layers 2/3, 5, and 6, as well as inhibitory interneurons in layers 1 and 4 of the hand knob.

Results

Axonal and dendritic terminals were polarized more than the soma in all neurons, with peak axonal and dendritic polarization of 0.92 mV and 0.21 mV, respectively, compared to peak somatic polarization of 0.07 mV for 1.8 mA M1-SO stimulation. Both montages generated regions of depolarization and hyperpolarization beneath the M1 anode; M1-SO produced slightly stronger, more diffuse polarization peaking in the central sulcus, while 4×1 HD produced higher peak polarization in the gyral crown. Simulating polarization by uniform local E-field approximated the spatial distribution of tDCS polarization but produced large errors in some regions.

Conclusions

Polarization of pre- and postsynaptic compartments of excitatory and inhibitory cortical neurons may play a significant role in tDCS neuromodulation. These effects cannot be predicted from the E-field distribution alone but rather require calculation of the neuronal response.

Atlas-based data integration for mapping the connections and architecture of the brain

Detailed knowledge about the neural connections among regions of the brain is key for advancing our understanding of normal brain function and changes that occur with aging and disease. Researchers use a range of experimental techniques to map connections at different levels of granularity in rodent animal models, but the results are often challenging to compare and integrate. Three-dimensional reference atlases of the brain provide new opportunities for cumulating, integrating, and reinterpreting research findings across studies. Here, we review approaches for integrating data describing neural connections and other modalities in rodent brain atlases and discuss how atlas-based workflows can facilitate brainwide analyses of neural network organization in relation to other facets of neuroarchitecture.

Seeing the Forest and Its Trees Together: Implementing 3D Light Microscopy Pipelines for Cell Type Mapping in the Mouse Brain

The brain is composed of diverse neuronal and non-neuronal cell types with complex regional connectivity patterns that create the anatomical infrastructure underlying cognition. Remarkable advances in neuroscience techniques enable labeling and imaging of these individual cell types and their interactions throughout intact mammalian brains at a cellular resolution allowing neuroscientists to examine microscopic details in macroscopic brain circuits. Nevertheless, implementing these tools is fraught with many technical and analytical challenges with a need for high-level data analysis. Here we review key technical considerations for implementing a brain mapping pipeline using the mouse brain as a primary model system. Specifically, we provide practical details for choosing methods including cell type specific labeling, sample preparation (e.g., tissue clearing), microscopy modalities, image processing, and data analysis (e.g., image registration to standard atlases). We also highlight the need to develop better 3D atlases with standardized anatomical labels and nomenclature across species and developmental time points to extend the mapping to other species including humans and to facilitate data sharing, confederation, and integrative analysis. In summary, this review provides key elements and currently available resources to consider while developing and implementing high-resolution mapping methods.

Whole-Brain Reconstruction of Neurons in the Ventral Pallidum Reveals Diverse Projection Patterns

The ventral pallidum (VP) integrates reward signals to regulate cognitive, emotional, and motor processes associated with motivational salience. Previous studies have revealed that the VP projects axons to many cortical and subcortical structures. However, descriptions of the neuronal morphologies and projection patterns of the VP neurons at the single neuron level are lacking, thus hindering the understanding of the wiring diagram of the VP. In this study, we used recently developed progress in robust sparse labeling and fluorescence micro-optical sectioning tomography imaging system (fMOST) to label mediodorsal thalamus-projecting neurons in the VP and obtain high-resolution whole-brain imaging data. Based on these data, we reconstructed VP neurons and classified them into three types according to their fiber projection patterns. We systematically compared the axonal density in various downstream centers and analyzed the soma distribution and dendritic morphologies of the various subtypes at the single neuron level. Our study thus provides a detailed characterization of the morphological features of VP neurons, laying a foundation for exploring the neural circuit organization underlying the important behavioral functions of VP.

Microfluidic platforms for single neuron analysis

Single-neuron actions are the basis of brain function, as clinical sequelae, neuronal dysfunction or failure for most of the central nervous system (CNS) diseases and injuries can be identified via tracing single-neurons. The bulk analysis methods tend to miscue critical information by assessing the population-averaged outcomes. However, its primary requisite in neuroscience to analyze single-neurons and to understand dynamic interplay of neurons and their environment. Microfluidic systems enable precise control over nano-to femto-liter volumes via adjusting device geometry, surface characteristics, and flow-dynamics, thus facilitating a well-defined micro-environment with spatio-temporal control for single-neuron analysis. The microfluidic platform not only offers a comprehensive landscape to study brain cell diversity at the level of transcriptome, genome, and/or epigenome of individual cells but also has a substantial role in deciphering complex dynamics of brain development and brain-related disorders. In this review, we highlight recent advances of microfluidic devices for single-neuron analysis, i.e., single-neuron trapping, single-neuron dynamics, single-neuron proteomics, single-neuron transcriptomics, drug delivery at the single-neuron level, single axon guidance, and single-neuron differentiation. Moreover, we also emphasize limitations and future challenges of single-neuron analysis by focusing on key performances of throughput and multiparametric activity analysis on microfluidic platforms.

TISSUE CLEARING

Tissue clearing of gross anatomical samples was first described over a century ago and has only recently found widespread use in the field of microscopy. This renaissance has been driven by the application of modern knowledge of optical physics and chemical engineering to the development of robust and reproducible clearing techniques, the arrival of new microscopes that can image large samples at cellular resolution and computing infrastructure able to store and analyze large data volumes. Many biological relationships between structure and function require investigation in three dimensions and tissue clearing therefore has the potential to enable broad discoveries in the biological sciences. Unfortunately, the current literature is complex and could confuse researchers looking to begin a clearing project. The goal of this Primer is to outline a modular approach to tissue clearing that allows a novice researcher to develop a customized clearing pipeline tailored to their tissue of interest. Further, the Primer outlines the required imaging and computational infrastructure needed to perform tissue clearing at scale, gives an overview of current applications, discusses limitations and provides an outlook on future advances in the field.

Predicting Synaptic Connectivity for Large-Scale Microcircuit Simulations Using Snudda

Simulation of large-scale networks of neurons is an important approach to understanding and interpreting experimental data from healthy and diseased brains. Owing to the rapid development of simulation software and the accumulation of quantitative data of different neuronal types, it is possible to predict both computational and dynamical properties of local microcircuits in a ‘bottom-up’ manner. Simulated data from these models can be compared with experiments and ‘top-down’ modelling approaches, successively bridging the scales. Here we describe an open source pipeline, using the software Snudda, for predicting microcircuit connectivity and for setting up simulations using the NEURON simulation environment in a reproducible way. We also illustrate how to further ‘curate’ data on single neuron morphologies acquired from public databases. This model building pipeline was used to set up a first version of a full-scale cellular level model of mouse dorsal striatum. Model components from that work are here used to illustrate the different steps that are needed when modelling subcortical nuclei, such as the basal ganglia.

Uncovering Statistical Links Between Gene Expression and Structural Connectivity Patterns in the Mouse Brain

Finding links between genes and structural connectivity is of the utmost importance for unravelling the underlying mechanism of the brain connectome. In this study we identify links between the gene expression and the axonal projection density in the mouse brain, by applying a modified version of the Linked ICA method to volumetric data from the Allen Institute for Brain Science for identifying independent sources of information that link both modalities at the voxel level. We performed separate analyses on sets of projections from the visual cortex, the caudoputamen and the midbrain reticular nucleus, and we determined those brain areas, injections and genes that were most involved in independent components that link both gene expression and projection density data, while we validated their biological context through enrichment analysis. We identified representative and literature-validated cortico-midbrain and cortico-striatal projections, whose gene subsets were enriched with annotations for neuronal and synaptic function and related developmental and metabolic processes. The results were highly reproducible when including all available projections, as well as consistent with factorisations obtained using the Dictionary Learning and Sparse Coding technique. Hence, Linked ICA yielded reproducible independent components that were preserved under increasing data variance. Taken together, we have developed and validated a novel paradigm for linking gene expression and structural projection patterns in the mouse mesoconnectome, which can power future studies aiming to relate genes to brain function.

Current Status of Tissue Clearing and the Path Forward in Neuroscience

Due to the complexity and limited availability of human brain tissues, for decades, pathologists have sought to maximize information gained from individual samples, based on which (patho)physiological processes could be inferred. Recently, new understandings of chemical and physical properties of biological tissues and multiple chemical profiling have given rise to the development of scalable tissue clearing methods allowing superior optical clearing of across-the-scale samples. In the past decade, tissue clearing techniques, molecular labeling methods, advanced laser scanning microscopes, and data visualization and analysis have become commonplace. Combined, they have made 3D visualization of brain tissues with unprecedented resolution and depth widely accessible. To facilitate further advancements and applications, here we provide a critical appraisal of these techniques. We propose a classification system of current tissue clearing and expansion methods that allows users to judge the applicability of individual ones to their questions, followed by a review of the current progress in molecular labeling, optical imaging, and data processing to demonstrate the whole 3D imaging pipeline based on tissue clearing and downstream techniques for visualizing the brain. We also raise the path forward of tissue-clearing-based imaging technology, that is, integrating with state-of-the-art techniques, such as multiplexing protein imaging, in situ signal amplification, RNA detection and sequencing, super-resolution imaging techniques, multiomics studies, and deep learning, for drawing the complete atlas of the human brain and building a 3D pathology platform for central nervous system disorders.

Prediction of a Cell-Class-Specific Mouse Mesoconnectome Using Gene Expression Data

Reconstructing brain connectivity at sufficient resolution for computational models designed to study the biophysical mechanisms underlying cognitive processes is extremely challenging. For such a purpose, a mesoconnectome that includes laminar and cell-class specificity would be a major step forward. We analyzed the ability of gene expression patterns to predict cell-class and layer-specific projection patterns and assessed the functional annotations of the most predictive groups of genes. To achieve our goal we used publicly available volumetric gene expression and connectivity data and we trained computational models to learn and predict cell-class and layer-specific axonal projections using gene expression data. Predictions were done in two ways, namely predicting projection strengths using the expression of individual genes and using the co-expression of genes organized in spatial modules, as well as predicting binary forms of projection. For predicting the strength of projections, we found that ridge (L2-regularized) regression had the highest cross-validated accuracy with a median r² score of 0.54 which corresponded for binarized predictions to a median area under the ROC value of 0.89. Next, we identified 200 spatial gene modules using a dictionary learning and sparse coding approach. We found that these modules yielded predictions of comparable accuracy, with a median r² score of 0.51. Finally, a gene ontology enrichment analysis of the most predictive gene groups resulted in significant annotations related to postsynaptic function. Taken together, we have demonstrated a prediction workflow that can be used to perform multimodal data integration to improve the accuracy of the predicted mesoconnectome and support other neuroscience use cases.

Axonal branch patterning and neuronal shape diversity: roles in developmental circuit assembly: Axonal branch patterning and neuronal shape diversity in developmental circuit assembly

Recent progress in human genetics and single cell sequencing rapidly expands the list of molecular factors that offer important new contributions to our understanding of brain wiring. Yet many new molecular factors are being discovered that have never been studied in the context of neuronal circuit development. This is clearly asking for increased efforts to better understand the developmental mechanisms of circuit assembly [1]. Moreover, recent studies characterizing the developmental causes of some psychiatric diseases show impressive progress in reaching cellular resolution in their analysis. They provide concrete support emphasizing the importance of axonal branching and synapse formation as a hotspot for potential defects. Inspired by these new studies we will discuss progress but also challenges in understanding how neurite branching and neuronal shape diversity itself impacts on specificity of neuronal circuit assembly. We discuss the idea that neuronal shape acquisition itself is a key specificity factor in neuronal circuit assembly.

Skeleton optimization of neuronal morphology based on three-dimensional shape restrictions

Background

Neurons are the basic structural unit of the brain, and their morphology is a key determinant of their classification. The morphology of a neuronal circuit is a fundamental component in neuron modeling. Recently, single-neuron morphologies of the whole brain have been used in many studies. The correctness and completeness of semimanually traced neuronal morphology are credible. However, there are some inaccuracies in semimanual tracing results. The distance between consecutive nodes marked by humans is very long, spanning multiple voxels. On the other hand, the nodes are marked around the centerline of the neuronal fiber, not on the centerline. Although these inaccuracies do not seriously affect the projection patterns that these studies focus on, they reduce the accuracy of the traced neuronal skeletons. These small inaccuracies will introduce deviations into subsequent studies that are based on neuronal morphology files.

Results

We propose a neuronal digital skeleton optimization method to evaluate and make fine adjustments to a digital skeleton after neuron tracing. Provided that the neuronal fiber shape is smooth and continuous, we describe its physical properties according to two shape restrictions. One restriction is designed based on the grayscale image, and the other is designed based on geometry. These two restrictions are designed to finely adjust the digital skeleton points to the neuronal fiber centerline. With this method, we design the three-dimensional shape restriction workflow of neuronal skeleton adjustment computation. The performance of the proposed method has been quantitatively evaluated using synthetic and real neuronal image data. The results show that our method can reduce the difference between the traced neuronal skeleton and the centerline of the neuronal fiber. Furthermore, morphology metrics such as the neuronal fiber length and radius become more precise.

Conclusions

This method can improve the accuracy of a neuronal digital skeleton based on traced results. The greater the accuracy of the digital skeletons that are acquired, the more precise the neuronal morphologies that are analyzed will be.

A Single-Neuron: Current Trends and Future Prospects

The brain is an intricate network with complex organizational principles facilitating a concerted communication between single-neurons, distinct neuron populations, and remote brain areas. The communication, technically referred to as connectivity, between single-neurons, is the center of many investigations aimed at elucidating pathophysiology, anatomical differences, and structural and functional features. In comparison with bulk analysis, single-neuron analysis can provide precise information about neurons or even sub-neuron level electrophysiology, anatomical differences, pathophysiology, structural and functional features, in addition to their communications with other neurons, and can promote essential information to understand the brain and its activity. This review highlights various single-neuron models and their behaviors, followed by different analysis methods. Again, to elucidate cellular dynamics in terms of electrophysiology at the single-neuron level, we emphasize in detail the role of single-neuron mapping and electrophysiological recording. We also elaborate on the recent development of single-neuron isolation, manipulation, and therapeutic progress using advanced micro/nanofluidic devices, as well as microinjection, electroporation, microelectrode array, optical transfection, optogenetic techniques. Further, the development in the field of artificial intelligence in relation to single-neurons is highlighted. The review concludes with between limitations and future prospects of single-neuron analyses.

Whole-Brain Profiling of Cells and Circuits in Mammals by Tissue Clearing and Light-Sheet Microscopy

Tissue clearing and light-sheet microscopy have a 100-year-plus history, yet these fields have been combined only recently to facilitate novel experiments and measurements in neuroscience. Since tissue-clearing methods were first combined with modernized light-sheet microscopy a decade ago, the performance of both technologies has rapidly improved, broadening their applications. Here, we review the state of the art of tissue-clearing methods and light-sheet microscopy and discuss applications of these techniques in profiling cells and circuits in mice. We examine outstanding challenges and future opportunities for expanding these techniques to achieve brain-wide profiling of cells and circuits in primates and humans. Such integration will help provide a systems-level understanding of the physiology and pathology of our central nervous system.

Differential Striatal Axonal Arborizations of the Intratelencephalic and Pyramidal-Tract Neurons: Analysis of the Data in the MouseLight Database

There exist two major types of striatum-targeting neocortical neurons, specifically, intratelencephalic (IT) neurons and pyramidal-tract (PT) neurons. Regarding their striatal projections, it was once suggested that IT axons are extended whereas PT axons are primarily focal. However, subsequent study with an increased number of well-stained extended axons concluded that such an apparent distinction was spurious due to limited sample size. Recent work using genetically labeled neurons reintroduced the differential spatial extent of the striatal projections of IT and PT neurons through population-level analyses, complemented by observations of single axons. However, quantitative IT vs. PT comparison of a large number of axons remained to be conducted. We analyzed the data of axonal end-points of 161 IT neurons and 33 PT neurons in the MouseLight database (http://ml-neuronbrowser.janelia.org/). The number of axonal end-points in the ipsilateral striatum exhibits roughly monotonically decreasing distributions in both neuron types. Excluding neurons with no ipsilateral end-point, the distributions of the logarithm of the number of ipsilateral end-points are considerably overlapped between IT and PT neurons, although the proportion of neurons having more than 50 ipsilateral end-points is somewhat larger in IT neurons than in PT neurons. Looking at more details, among IT subpopulations in the secondary motor area (MOs), layer 5 neurons and bilateral striatum-targeting layer 2/3 neurons, but not contralateral striatum-non-targeting layer 2/3 neurons, have a larger number of ipsilateral end-points than MOs PT neurons. We also found that IT ipsilateral striatal axonal end-points are on average more widely distributed than PT end-points, especially in the medial-lateral direction. These results indicate that IT and PT striatal axons differ in the frequencies and spatial extent of end-points while there are wide varieties within each neuron type.

Focal neocortical lesions impair distant neuronal information processing

KEY POINTS

Parts of the fields of neuroscience and neurology consider the neocortex to be a functionally parcelled structure. Viewed through such a conceptual filter, there are multiple clinical observations after localized stroke lesions that seem paradoxical. We tested the effect that localized stroke-like lesions have on neuronal information processing in a part of the neocortex that is distant to the lesion using animal experiments. We find that the distant lesion degrades the quality of neuronal information processing of tactile input patterns in primary somatosensory cortex. The findings suggest that even the processing of primary sensory information depends on an intact neocortical network, with the implication that all neocortical processing may rely on widespread interactions across large parts of the cortex.

ABSTRACT

Recent clinical studies report a surprisingly weak relationship between the location of cortical brain lesions and the resulting functional deficits. From a neuroscience point of view, such findings raise questions as to what extent functional localization applies in the neocortex and to what extent the functions of different regions depend on the integrity of others. Here we provide an in-depth analysis of the changes in the function of the neocortical neuronal networks after distant focal stroke-like lesions in the anaesthetized rat. Using a recently introduced high resolution analysis of neuronal information processing, consisting of pre-set spatiotemporal patterns of tactile afferent activation against which the neuronal decoding performance can be quantified, we found that stroke-like lesions in distant parts of the cortex significantly degraded the decoding performance of individual neocortical neurons in the primary somatosensory cortex (decoding performance decreased from 30.9% to 24.2% for n = 22 neurons, Wilcoxon signed rank test, P = 0.028). This degrading effect was not due to changes in the firing frequency of the neuron (Wilcoxon signed rank test, P = 0.499) and was stronger the higher the decoding performance of the neuron, indicating a specific impact on the information processing capacity in the cortex. These findings suggest that even primary sensory processing depends on widely distributed cortical networks and could explain observations of focal stroke lesions affecting a large range of functions.