Search

To decipher the organizational logic of complex brain circuits, it is important to chart long-distance pathways while preserving micron-level accuracy of local network. However, mapping the neuronal projections with individual-axon resolution in the large and complex primate brain is still challenging. Herein, we describe a highly efficient pipeline for three-dimensional mapping of the entire macaque brain with subcellular resolution. The pipeline includes a novel poly-N-acryloyl glycinamide (PNAGA)-based embedding method for long-term structure and fluorescence preservation, high-resolution and rapid whole-brain optical imaging, and image post-processing. The cytoarchitectonic information of the entire macaque brain was acquired with a voxel size of 0.32 μm × 0.32 μm × 10 μm, showing its anatomical structure with cell distribution, density, and shape. Furthermore, thanks to viral labeling, individual long-distance projection axons from the frontal cortex were for the first time reconstructed across the entire brain hemisphere with a voxel size of 0.65 μm × 0.65 μm × 3 μm. Our results show that individual cortical axons originating from the prefrontal cortex simultaneously target multiple brain regions, including the visual cortex, striatum, thalamus, and midbrain. This pipeline provides an efficient method for cellular and circuitry investigation of the whole macaque brain with individual-axon resolution, and can shed light on brain function and disorders.

An essential step toward understanding brain function is to establish a structural framework with cellular resolution on which multi-scale datasets spanning molecules, cells, circuits and systems can be integrated and interpreted1. Here, as part of the collaborative Brain Initiative Cell Census Network (BICCN), we derive a comprehensive cell type-based anatomical description of one exemplar brain structure, the mouse primary motor cortex, upper limb area (MOp-ul). Using genetic and viral labelling, barcoded anatomy resolved by sequencing, single-neuron reconstruction, whole-brain imaging and cloud-based neuroinformatics tools, we delineated the MOp-ul in 3D and refined its sublaminar organization. We defined around two dozen projection neuron types in the MOp-ul and derived an input–output wiring diagram, which will facilitate future analyses of motor control circuitry across molecular, cellular and system levels. This work provides a roadmap towards a comprehensive cellular-resolution description of mammalian brain architecture., Multi-modal analysis is used to generate a 3D atlas of the upper limb area of the mouse primary motor cortex, providing a framework for future studies of motor control circuitry.

The cortico–basal ganglia–thalamo–cortical loop is one of the fundamental network motifs in the brain. Revealing its structural and functional organization is critical to understanding cognition, sensorimotor behaviour, and the natural history of many neurological and neuropsychiatric disorders. Classically, this network is conceptualized to contain three information channels: motor, limbic and associative1–4. Yet this three-channel view cannot explain the myriad functions of the basal ganglia. We previously subdivided the dorsal striatum into 29 functional domains on the basis of the topography of inputs from the entire cortex5. Here we map the multi-synaptic output pathways of these striatal domains through the globus pallidus external part (GPe), substantia nigra reticular part (SNr), thalamic nuclei and cortex. Accordingly, we identify 14 SNr and 36 GPe domains and a direct cortico-SNr projection. The striatonigral direct pathway displays a greater convergence of striatal inputs than the more parallel striatopallidal indirect pathway, although direct and indirect pathways originating from the same striatal domain ultimately converge onto the same postsynaptic SNr neurons. Following the SNr outputs, we delineate six domains in the parafascicular and ventromedial thalamic nuclei. Subsequently, we identify six parallel cortico–basal ganglia–thalamic subnetworks that sequentially transduce specific subsets of cortical information through every elemental node of the cortico– basal ganglia–thalamic loop. Thalamic domains relay this output back to the originating corticostriatal neurons of each subnetwork in a bona fide closed loop., Mesoscale connectomic mapping of the cortico–basal ganglia–thalamic network reveals key architectural and information processing features.

The impressive functions of the brain rely on an extensive connectivity matrix between specific neurons, the architecture of which is frequently characterized by one brain nucleus/region connecting to multiple targets, either via collaterals of the same projection neuron or several, differentially specified neurons. Delineating the fine architecture of projection neuron subsets in a specific brain region could greatly facilitate its circuit, computational, and functional resolution. Here, we developed multiple fluorescent rabies viruses (RV) to delineate the fine organization of corticothalamic projection neuron subsets in the primary visual cortex (V1). By simultaneously retrograde labeling multiple distinct subsets of corticothalamic projection neurons in V1 from their target nuclei in thalamus (dLGN, LP, LD), we observed that V1-dLGN corticothalamic projection neurons were densely concentrated in layer VI, except for several sparsely scattered neurons in layer V, while V1-LP and V1-LD corticothalamic projection neurons were localized to both layers V and VI. Meanwhile, we observed a fraction of V1 corticothalamic projection neurons targeting two thalamic nuclei, which was further confirmed by fMOST whole-brain imaging. The multiple fluorescent RV tracing tools can be extensively applied to resolve the architecture of projection neuron subsets in certain brain regions, with a strong potential to delineate the computational and functional organization of these brain regions.

The medial prefrontal cortex (mPFC) contains populations of GABAergic interneurons that play different roles in cognition and emotion. Their local and long-range inputs are incompletely understood. We used monosynaptic rabies viral tracers in combination with fluorescence micro-optical sectioning tomography to generate a whole-brain atlas of direct long-range inputs to GABAergic interneurons in the mPFC of male mice. We discovered that three subtypes of GABAergic interneurons in two areas of the mPFC are innervated by same upstream areas. Input from subcortical upstream areas includes cholinergic neurons from the basal forebrain and serotonergic neurons (which co-release glutamate) from the raphe nuclei. Reconstruction of single-neuron morphology revealed novel substantia innominata-anteromedial thalamic nucleus-mPFC and striatum-anteromedial thalamic nucleus-mPFC circuits. Based on the projection logic of individual neurons, we classified cortical and hippocampal input neurons into several types. This atlas provides the anatomical foundation for understanding the functional organization of the mPFC.

Many people affected by fragile X syndrome (FXS) and autism spectrum disorders have sensory processing deficits, such as hypersensitivity to auditory, tactile, and visual stimuli. Like FXS in humans, loss of Fmr1 in rodents also cause sensory, behavioral, and cognitive deficits. However, the neural mechanisms underlying sensory impairment, especially vision impairment, remain unclear. It remains elusive whether the visual processing deficits originate from corrupted inputs, impaired perception in the primary sensory cortex, or altered integration in the higher cortex, and there is no effective treatment. In this study, we used a genetic knockout mouse model (Fmr1KO), in vivo imaging, and behavioral measurements to show that the loss of Fmr1 impaired signal processing in the primary visual cortex (V1). Specifically, Fmr1KO mice showed enhanced responses to low-intensity stimuli but normal responses to high-intensity stimuli. This abnormality was accompanied by enhancements in local network connectivity in V1 microcircuits and increased dendritic complexity of V1 neurons. These effects were ameliorated by the acute application of GABAA receptor activators, which enhanced the activity of inhibitory neurons, or by reintroducing Fmr1 gene expression in knockout V1 neurons in both juvenile and young-adult mice. Overall, V1 plays an important role in the visual abnormalities of Fmr1KO mice and it could be possible to rescue the sensory disturbances in developed FXS and autism patients.

The dorsal raphe nucleus (DR) and median raphe nucleus (MR) contain populations of glutamatergic and GABAergic neurons regulating diverse behavioral functions. Their whole-brain input-output circuits remain incompletely understood. We used viral tracing combined with fluorescence micro-optical sectioning tomography to generate a comprehensive whole-brain atlas of inputs and outputs of glutamatergic and GABAergic neurons in the DR and MR. We discovered that these neurons receive inputs from similar upstream brain regions. The glutamatergic and GABAergic neurons in the same raphe nucleus have divergent projection patterns with differences in critical brain regions. Specifically, MR glutamatergic neurons project to the lateral habenula via multiple pathways. Correlation and cluster analysis indicated that glutamatergic and GABAergic neurons in the same raphe nucleus receive inputs from heterogeneous neurons in upstream brain regions and send different collateral projections. This connectivity atlas provides insights into the cell heterogeneity, anatomical connectivity and behavioral functions of the raphe nucleus.

The dorsal raphe nucleus (DR) and median raphe nucleus (MR) contain populations of glutamatergic and GABAergic neurons that regulate diverse behavioral functions. However, their whole-brain input-output circuits remain incompletely elucidated. We used viral tracing combined with fluorescence micro-optical sectioning tomography to generate a comprehensive whole-brain atlas of inputs and outputs of glutamatergic and GABAergic neurons in the DR and MR. We found that these neurons received inputs from similar upstream brain regions. The glutamatergic and GABAergic neurons in the same raphe nucleus had divergent projection patterns with differences in critical brain regions. Specifically, MR glutamatergic neurons projected to the lateral habenula through multiple pathways. Correlation and cluster analysis revealed that glutamatergic and GABAergic neurons in the same raphe nucleus received heterogeneous inputs and sent different collateral projections. This connectivity atlas further elucidates the anatomical architecture of the raphe nuclei, which could facilitate better understanding of their behavioral functions.

Axonopathy is a pathological feature observed in both Alzheimer’s disease (AD) patients and animal models. However, identifying the temporal and regional progression of axonopathy during AD development remains elusive. Using the fluorescence micro-optical sectioning tomography system, we acquired whole-brain datasets in the early stage of 5xFAD/Thy1-GFP-M mice. We reported that among GFP labeled axons, GFP-positive axonopathy first formed in the lateral septal nucleus, subiculum, and medial mammillary nucleus. The axonopathy further increased in most brain regions during aging. However, most of the axonopathic varicosities disappeared significantly in the medial mammillary nucleus after 8 weeks old. Continuous three-dimensional datasets showed that axonopathy in the medial mammillary nucleus was mainly located on axons from hippocampal GFP-positive neurons. Using the rabies viral tracer in combination with immunohistochemistry, we found that axons in the medial mammillary nucleus from the subiculum were susceptible to lesions that prior to the occurrence of behavioral disorders. In conclusion, we created an early-stage spatiotemporal map of axonopathy in 5xFAD/Thy1-GFP-M mice and identified specific neural circuits which are vulnerable to axon lesions in an AD mouse model. These findings underline the importance of early interventions for AD, and may contribute to the understanding of its progression and its early symptom treatment.

Ever since the seminal findings of Ramon y Cajal, dendritic and axonal morphology has been recognized as a defining feature of neuronal types. Yet our knowledge concerning the diversity of neuronal morphologies, in particular distal axonal projection patterns, is extremely limited. To systematically obtain single neuron full morphology on a brain-wide scale, we established a platform with five major components: sparse labeling, whole-brain imaging, reconstruction, registration, and classification. We achieved sparse, robust and consistent fluorescent labeling of a wide range of neuronal types by combining transgenic or viral Cre delivery with novel transgenic reporter lines. We acquired high-resolution whole-brain fluorescent images from a large set of sparsely labeled brains using fluorescence micro-optical sectioning tomography (fMOST). We developed a set of software tools for efficient large-volume image data processing, registration to the Allen Mouse Brain Common Coordinate Framework (CCF), and computer-assisted morphological reconstruction. We reconstructed and analyzed the complete morphologies of 1,708 neurons from the striatum, thalamus, cortex and claustrum. Finally, we classified these cells into multiple morphological and projection types and identified a set of region-specific organizational rules of long-range axonal projections at the single cell level. Specifically, different neuron types from different regions follow highly distinct rules in convergent or divergent projection, feedforward or feedback axon termination patterns, and between-cell homogeneity or heterogeneity. Major molecularly defined classes or types of neurons have correspondingly distinct morphological and projection patterns, however, we also identify further remarkably extensive morphological and projection diversity at more fine-grained levels within the major types that cannot presently be accounted for by preexisting transcriptomic subtypes. These insights reinforce the importance of full morphological characterization of brain cell types and suggest a plethora of ways different cell types and individual neurons may contribute to the function of their respective circuits.

The cortico-basal ganglia-thalamic loop is one of the fundamental network motifs in the brain. Revealing its structural and functional organization is critical to understanding cognition, sensorimotor behavior, and the natural history of many neurological and neuropsychiatric diseases. Classically, the basal ganglia is conceptualized to contain three primary information output channels: motor, limbic, and associative. However, given the roughly 65 cortical areas and two dozen thalamic nuclei that feed into the dorsal striatum, a three-channel view is overly simplistic for explaining the myriad functions of the basal ganglia. Recent works from our lab and others have subdivided the dorsal striatum into numerous functional domains based on convergent and divergent inputs from the cortex and thalamus. To complete this work, we generated a comprehensive data pool of ∼700 injections placed across the striatum, external globus pallidus (GPe), substantia nigra pars reticulata (SNr), thalamic nuclei, and cortex. We identify 14 domains of SNr, 36 in the GPe, and 6 in the parafascicular and ventromedial thalamic nuclei. Subsequently, we identify 6 parallel cortico-basal ganglia-thalamic subnetworks that sequentially transduce specific subsets of cortical information with complex patterns of convergence and divergence through every elemental node of the entire cortico-basal ganglia loop. These experiments reveal multiple important novel features of the cortico-basal ganglia network motif. The prototypical sub-network structure is characterized by a highly interconnected nature, with cortical information processing through one or more striatal nodes, which send a convergent output to the SNr and a more parallelized output to the GPe; the GPe output then converges with the SNr. A domain of the thalamus receives the nigral output, and is interconnected with both the striatal domains and the cortical areas that filter into its nigral input source. This study provides conceptual advancement of our understanding of the structural and functional organization of the classic cortico-basal ganglia network.

An essential step toward understanding brain function is to establish a cellular-resolution structural framework upon which multi-scale and multi-modal information spanning molecules, cells, circuits and systems can be integrated and interpreted. Here, through a collaborative effort from the Brain Initiative Cell Census Network (BICCN), we derive a comprehensive cell type-based description of one brain structure - the primary motor cortex upper limb area (MOp-ul) of the mouse. Applying state-of-the-art labeling, imaging, computational, and neuroinformatics tools, we delineated the MOp-ul within the Mouse Brain 3D Common Coordinate Framework (CCF). We defined over two dozen MOp-ul projection neuron (PN) types by their anterograde targets; the spatial distribution of their somata defines 11 cortical sublayers, a significant refinement of the classic notion of cortical laminar organization. We further combine multiple complementary tracing methods (classic tract tracing, cell type-based anterograde, retrograde, and transsynaptic viral tracing, high-throughput BARseq, and complete single cell reconstruction) to systematically chart cell type-based MOp input-output streams. As PNs link distant brain regions at synapses as well as host cellular gene expression, our construction of a PN type resolution MOp-ul wiring diagram will facilitate an integrated analysis of motor control circuitry across the molecular, cellular, and systems levels. This work further provides a roadmap towards a cellular resolution description of mammalian brain architecture.

Dissection of the neural circuits of the cerebral cortex is essential for studying mechanisms underlying brain function. Herein, combining a retrograde rabies tracing system with fluorescent micro-optical sectional tomography, we investigated long-range input neurons of corticotropin-releasing hormone containing neurons in the six main cortical areas, including the prefrontal, somatosensory, motor, auditory, and visual cortices. The whole brain distribution of input neurons showed similar patterns to input neurons distributed mainly in the adjacent cortical areas, thalamus, and basal forebrain. Reconstruction of continuous three-dimensional datasets showed the anterior and middle thalamus projected mainly to the rostral cortex whereas the posterior and lateral projected to the caudal cortex. In the basal forebrain, immunohistochemical staining showed these cortical areas received afferent information from cholinergic neurons in the substantia innominata and lateral globus pallidus, whereas cholinergic neurons in the diagonal band nucleus projected strongly to the prefrontal and visual cortex. Additionally, dense neurons in the zona incerta and ventral hippocampus were found to project to the prefrontal cortex. These results showed general patterns of cortical input circuits and unique connection patterns of each individual area, allowing for valuable comparisons among the organisation of different cortical areas and new insight into cortical functions.

The segmentation of brain region contours in three dimensions is critical for the analysis of different brain structures, and advanced approaches are emerging continuously within the field of neurosciences. With the development of high-resolution micro-optical imaging, whole-brain images can be acquired at the cellular level. However, brain regions in microscopic images are aggregated by discrete neurons with blurry boundaries, the complex and variable features of brain regions make it challenging to accurately segment brain regions. Manual segmentation is a reliable method, but is unrealistic to apply on a large scale. Here, we propose an automated brain region segmentation framework, DeepBrainSeg, which is inspired by the principle of manual segmentation. DeepBrainSeg incorporates three feature levels to learn local and contextual features in different receptive fields through a dual-pathway convolutional neural network (CNN), and to provide global features of localization by image registration and domain-condition constraints. Validated on biological datasets, DeepBrainSeg can not only effectively segment brain-wide regions with high accuracy (Dice ratio > 0.9), but can also be applied to various types of datasets and to datasets with noises. It has the potential to automatically locate information in the brain space on the large scale.

Summary Ventral hippocampus (vHIP) and medial prefrontal cortex (mPFC) are both critical regions for social behaviors. However, how their interactions affect social behavior is not well understood. By viral tracing, optogenetics, chemogenetics, and fiber photometry, we demonstrated that inhibition of vHIP or direct projections from vHIP to mPFC impaired social memory expression. Via rabies retrograde tracing, we found that all three major GABAergic neurons in mPFC received direct inputs from vHIP. Activation of parvalbumin positive (PV+) neurons in mPFC but not somatostatin positive (SST+) neurons can rescue the social memory impairment caused by vHIP inhibition. Furthermore, fiber photometry results demonstrated that social behaviors preferentially recruited PV+ neurons and inhibition of hippocampal neurons disrupted the activity of PV+ neurons during social interactions. These results revealed a new mechanism of how vHIP and mPFC regulate social behavior in complementarity with the existing neural circuitry mechanism., Graphical Abstract, Highlights • Inhibition of vHIP or direct vHIP-mPFC pathway disrupts social memory expression • Social behaviors preferentially recruit PV+ neurons in mPFC • Activation of PV+ neurons in mPFC rescue the vHIP-related impairment of social memory • Inhibition of VIP+ neurons in mPFC rescue the vHIP-related impairment of social memory, Behavioral Neuroscience; Neuroscience; Systems Neuroscience

The caudal forelimb area (CFA) of the mouse cortex is essential in many forelimb movements, and diverse types of GABAergic interneuron in the CFA are distinct in the mediation of cortical inhibition in motor information processing. However, their long-range inputs remain unclear. In the present study, we combined the monosynaptic rabies virus system with Cre driver mouse lines to generate a whole-brain map of the inputs to three major inhibitory interneuron types in the CFA. We discovered that each type was innervated by the same upstream areas, but there were quantitative differences in the inputs from the cortex, thalamus, and pallidum. Comparing the locations of the interneurons in two sub-regions of the CFA, we discovered that their long-range inputs were remarkably different in distribution and proportion. This whole-brain mapping indicates the existence of parallel pathway organization in the forelimb subnetwork and provides insight into the inhibitory processes in forelimb movement to reveal the structural architecture underlying the functions of the CFA. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (10.1007/s12264-019-00458-6) contains supplementary material, which is available to authorized users.

The elegant functions of the brain are facilitated by sophisticated connections between neurons, the architecture of which is frequently characterized by one nucleus connecting to multiple targets via projection neurons. Delineating the sub-nucleus fine architecture of projection neurons in a certain nucleus could greatly facilitate its circuit, computational, and functional resolution. Here, we developed multi-fluorescent rabies virus to delineate the fine organization of corticothalamic projection neuron subsets in the primary visual cortex (V1). By simultaneously labeling multiple distinct subsets of corticothalamic projection neurons in V1 from their target nuclei in thalamus (dLGN, LP, LD), we observed that V1-dLGN corticothalamic neurons were densely concentrated in layer VI, except for several sparsely scattered neurons in layer V, while V1-LP and V1-LD corticothalamic neurons were localized to both layers V and VI. Meanwhile, we observed a fraction of V1 corticothalamic neurons targeting multiple thalamic nuclei, which was further confirmed by fMOST whole-brain imaging. We further conceptually proposed an upgraded sub-nucleus tracing system with higher throughput (21 subsets) for more complex architectural tracing. The multi-fluorescent RV tracing tool can be extensively applied to resolve architecture of projection neuron subsets, with a strong potential to delineate the computational and functional organization of these nuclei.

Mapping the cytoarchitecture of the whole brain can reveal the organizational logic of neural systems. However, this remains a significant challenge, especially for gyrencephalic brains with a large volume. Here we propose an integrated pipeline for generating a cytoarchitectonic atlas with single-cell resolution of the whole brain. To analyze a large-volume brain, we used a modified en-bloc Nissl staining protocol to achieve uniform staining of large-scale brain specimens from ferret (Mustela putorius furo). By combining whole-brain imaging and big data processing, we established strategies for parsing cytoarchitectural information at a voxel resolution of 0.33 μm × 0.33 μm × 1 μm and terabyte-scale data analysis. Using the cytoarchitectonic datasets for adult ferret brain, we identified giant pyramidal neurons in ferret brains and provide the first report of their morphological diversity, neurochemical phenotype, and distribution patterns in the whole brain in three dimensions. This pipeline will facilitate studies on the organization and development of the mammalian brains, from that of rodents to the gyrencephalic brains of ferret and even primates.

The excitatory neurons in the visual cortex are of great significance for us in understanding brain functions. However, the diverse neuron types and their morphological properties have not been fully deciphered. In this paper, we applied the brain-wide positioning system (BPS) to image the entire brain of two Thy1-eYFP H-line male mice at 0.2 μm × 0.2 μm × 1 μm voxel resolution. A total of 103 neurons were reconstructed in layers 5 and 6 of the visual cortex with single-axon-level resolution. Based on the complete topology of neurons and the inherent positioning function of the imaging method, we classified the observed neurons into six types according to their apical dendrites and somata location: star pyramidal cells in layer 5 (L5-sp), slender-tufted pyramidal cells in layer 5 (L5-st), tufted pyramidal cells in layer 5 (L5-tt), spiny stellate-like cells in layer 6 (L6-ss), star pyramidal cells in layer 6 (L6-sp), and slender-tufted pyramidal cells in layer 6 (L6-st). By examining the axonal projection patterns of individual neurons, they can be categorized into three modes: ipsilateral circuit connection neurons, callosal projection neurons and corticofugal projection neurons. Correlating the two types of classifications, we have found that there are at least two projection modes comprised in the former defined neuron types except for L5-tt. On the other hand, each projection mode may consist of four dendritic types defined in this study. The axon projection mode only partially correlates with the apical dendrite feature. This work has demonstrated a paradigm for resolving the visual cortex through single-neuron-level quantification and has shown potential to be extended to reveal the connectome of other defined sensory and motor systems.

Paraffin embedding is widely used in microscopic imaging for preparing biological specimens. However, owing to significant fluorescence quenching during the embedding process, it is not compatible with fluorescent-labeling techniques, such as transgenic and viral labeling using green fluorescent protein (GFP). Here, we investigate the quenching mechanism and optimize the embedding process to improve the preservation of fluorescence intensity. The results show that dehydration is the main reason for fluorescence quenching during paraffin embedding, caused by the full denaturation of GFP molecules in ethyl alcohol. To evaluate fluorescent and morphological preservation, we modified the embedding process using tertiary butanol (TBA) instead of ethyl alcohol. Fluorescence intensity following TBA dehydration increased 12.08-fold of that observed in the traditional method. We obtained uniform fluorescence maintenance throughout the whole mouse brain, while the continuous apical dendrites, spines, and axon terminals were shown evenly within the cortex, hippocampus, and the amygdala. Moreover, we embedded a whole rat brain labeled with AAV in the prelimbic cortex (Prl). With the axon terminals in different areas, such as the caudate putamen, thalamus, and pyramidal tract, the results showed a continuous tract of Prl neurons throughout the whole brain. This method was also suitable for tdTomota labeled samples. These findings indicate that this modified embedding method could be compatible with GFP and provides a potential turning point for applications in the fluorescent labeling of samples.

Long-range neuronal circuits play an important role in motor and sensory information processing. Determining direct synaptic inputs of excited and inhibited neurons is important for understanding the circuit mechanisms involved in regulating movement. Here, we used the monosynaptic rabies tracing technique, combined with fluorescent micro-optical sectional tomography, to characterize the brain-wide input to the motor cortex (MC). The whole brain dataset showed that the main excited and inhibited neurons in the MC received inputs from similar brain regions with a quantitative difference. With 3D reconstruction we found that the distribution of input neurons, that target the primary and secondary MC, had different patterns. In the cortex, the neurons projecting to the primary MC mainly distributed in the lateral and anterior portion, while those to the secondary MC distributed in the medial and posterior portion. The input neurons in the subcortical areas also showed the topographic shift model, as in the thalamus, the neurons distributed as outer and inner shells while the neurons in the claustrum and amygdala were in the ventral and dorsal part, respectively. These results lay the anatomical foundation to understanding the organized pattern of motor circuits and the functional differences between the primary and secondary MC.

Neural information is encoded by action potentials delivered by neurons. Which component of neural activity constitutes the basic unit carrying information is still a controversial issue. In this paper, stimulation experiments using a network of hippocampal neurons cultured on a multi-electrode array are used to investigate this issue. The results show that for a set of pulse stimuli with varying voltage amplitude, the neuronal response to the fronto-potential sequence encodes more information through the moment of fronto-potential delivery than the number of fronto-potential deliveries, and that neurons at each locus are encoded independently of each other. After the addition of bicuculline inhibited the GABAA receptors, the information capacity decreased and the temporal resolution decreased, but the neurons at each site were still encoded independently. The results suggest that the encoding of stimulus amplitude in the cultured hippocampal neuronal network is better with spike timing than with count, and the effect of timing encoding is dependent on GABAA receptors.

Summary In response to stressors, individuals adopt different behavioral styles, which are essential for survival and form the basis of differential susceptibility to stress-related disorders. Corticotropin-releasing factor (CRF) and the medial prefrontal cortex (mPFC) have predominantly been studied in behavioral response to stress, while the role of mPFC CRF neurons is poorly understood. Using morphology, electrophysiology, and calcium imaging approaches, we characterized mPFC CRF neurons as a unique subtype of GABAergic inhibitory interneurons that were directly engaged in the tail suspension challenge. Genetic ablation or chemogenetic inhibition of dorsal mPFC (dmPFC) CRF neurons increased immobility under the tail-suspension and forced-swimming challenges and induced social avoidance behavior, whereas activation had the opposite effect on the same measures. Furthermore, increasing CRF neuronal activity promoted durable resilience to repeated social defeat stress. These results uncover a critical role of mPFC CRF interneurons in bidirectionally controlling motivated behavioral style selection under stress.

Resin embedding has been widely used for precise imaging of fluorescently labeled biological samples with optical and electron microscopy. The low preservation rate of fluorescence, especially for red fluorescent proteins, has limited the application of resin embedding in multifluorescent protein-labeled samples. Here, we optimized the embedding method to retain the intensity of multiple fluorescent proteins during resin embedding. By reducing the polymerization temperature from 50 to 35°C and adding a fluorescent protein protection reagent during the embedding process, we successfully increased the fluorescence preservation rate by nearly twofold for red fluorescent proteins, including tdTomato, mCherry, and DsRed. Meanwhile, the background fluorescence decreased significantly in the optimized embedding method. This method is suitable not only for red fluorescent protein-labeled samples but also for blue (BFP) and green fluorescent protein (GFP)-labeled samples. We embedded brains labeled with BFP, DsRed, and GFP via AAV and rabies virus and acquired the distribution of input neurons to different cortical areas. With GFP/tdTomato double-labeled samples in resin, we obtained the cholinergic projectome of the pedunculopontine tegmental nucleus (PPTg) and the distribution of cholinergic neurons at single-neuron resolution in the whole brain simultaneously. Input cholinergic terminals from the PPTg were found to innervate the cholinergic soma and fiber in the neocortex, basal forebrain and brainstem, indicating that local cholinergic neurons received long-range cholinergic modulation from the midbrain. Our optimized method is useful for embedding multicolor fluorescent protein-labeled samples to acquire multidimensional structural information on neural circuits at single-neuron resolution in the whole brain.

Urination (also called micturition) is thought to be regulated by a neural network that is distributed in both subcortical and cortical regions. Previously, urination-related neurons have been identified in subcortical structures such as the pontine micturition center (also known as Barrington's nucleus). However, the origin of the descending cortical pathway and how it interfaces with this subcortical circuit to permit voluntary initiation of urination remain elusive. Here we identified a small cluster of layer 5 neurons in the primary motor cortex whose activities tightly correlate with the onset of urination in freely behaving mice and increase dramatically during territorial marking. Optogenetically activating these neurons elicits contraction of the bladder and initiates urination, through their projections to the pontine micturition center, while silencing or ablating them impairs urination and causes retention of urine. Together these results reveal a novel cortical component upstream of the pontine micturition center that is critically involved in urination.

The cholinergic system in the brain plays crucial roles in regulating sensory and motor functions as well as cognitive behaviors by modulating neuronal activity. Understanding the organization of the cholinergic system requires a complete map of cholinergic neurons and their axon arborizations throughout the entire brain at the level of single neurons. Here, we report a comprehensive whole-brain atlas of the cholinergic system originating from various cortical and subcortical regions of the mouse brain. Using genetically labeled cholinergic neurons together with whole-brain reconstruction of optical images at 2-μm resolution, we obtained quantification of the number and soma volume of cholinergic neurons in 22 brain areas. Furthermore, by reconstructing the complete axonal arbors of fluorescently labeled single neurons from a subregion of the basal forebrain at 1-μm resolution, we found that their projections to the forebrain and midbrain showed neuronal subgroups with distinct projection specificity and diverse arbor distribution within the same projection area. These results suggest the existence of distinct subtypes of cholinergic neurons that serve different regulatory functions in the brain and illustrate the usefulness of complete reconstruction of neuronal distribution and axon projections at the mesoscopic level.

Understanding amazingly complex brain functions and pathologies requires a complete cerebral vascular atlas in stereotaxic coordinates. Making a precise atlas for cerebral arteries and veins has been a century-old objective in neuroscience and neuropathology. Using micro-optical sectioning tomography (MOST) with a modified Nissl staining method, we acquired five mouse brain data sets containing arteries, veins, and microvessels. Based on the brain-wide vascular spatial structures and brain regions indicated by cytoarchitecture in one and the same mouse brain, we reconstructed and annotated the vascular system atlas of both arteries and veins of the whole mouse brain for the first time. The distributing patterns of the vascular system within the brain regions were acquired and our results show that the patterns of individual vessels are different from each other. Reconstruction and statistical analysis of the microvascular network, including derivation of quantitative vascular densities, indicate significant differences mainly in vessels with diameters less than 8 μm and large than 20 μm across different brain regions. Our precise cerebral vascular atlas provides an important resource and approach for quantitative studies of brain functions and diseases.

Corticotropin-releasing hormone (CRH), with widespread expression in the brain, plays a key role in modulating a series of behaviors, including anxiety, arousal, motor function, learning and memory. Previous studies have focused on some brain regions with densely distributed CRH neurons such as paraventricular hypothalamic nucleus (PVH) and bed nuclei of the stria terminalis (BST) and revealed some basic structural and functional knowledge of CRH neurons. However, there is no systematic analysis of brain-wide distribution of CRH neurons. Here, we performed a comprehensive study of CRH neurons in CRH-IRES-Cre;Ai3 mice via automatic imaging and stereoscopic cell counting in a whole mouse brain. We acquired four datasets of the CRH distributions with co-localized cytoarchitecture at a voxel resolution of 0.32 × 0.32 × 2 μm using brain-wide positioning system (BPS). Next, we precisely located and counted the EYFP-labeled neurons in different regions according to propidium iodide (PI) counterstained anatomical reference using Neuronal Global Position System (NeuroGPS). In particular, dense EYFP expression was found in piriform area (PIR), BST, central amygdalar nucleus (CEA), PVH, Barrington’s nucleus (B) and inferior olivary complex (IO). Considerable CRH neurons were also found in main olfactory bulb (MOB), medial preoptic nucleus (MPN), pontine gray (PG), tegmental reticular nucleus (TRN), external cuneate nucleus (ECU) and midline thalamus. We reconstructed and compared the soma morphology of CRH neurons in 11 brain regions. The results demonstrated that CRH neurons had regional diversities of both cell distribution and soma morphology. This anatomical knowledge enhances the current understanding of the functions of CRH neurons. These results also demonstrated the ability of our platform to accurately orient, reconstruct and count neuronal somas in three-dimension for type-specific neurons in the whole brain, making it feasible to answer the fundamental neuroscience question of exact numbers of various neurons in the whole brain.

Working memory (WM) plays a crucial role in human brain functions. The application of brain connectivity analysis helps to understand the brain network properties in WM. Combination of functional and effective connectivity can provide new insights for exploring network attributes. Nevertheless, few studies have combined these two modes in WM. Near-infrared spectroscopy was used to investigate the connectivity properties in the prefrontal cortex (PFC) during n-back (0-back and 2-back) tasks by combining functional and effective connectivity analysis. Our results demonstrated that the PFC network showed small-world properties in both WM tasks. The characteristic path length was significantly longer in the 2-back task than in the 0-back task, while there was no obvious difference in the clustering coefficient between two tasks. Regarding the effective connectivity, the Granger causality (GC) was higher for right PFC→left PFC than for left PFC→right PFC in the 2-back task. Compared with the 0-back task, GC of right PFC→left PFC was higher in the 2-back task. Our findings show that, along with memory load increase, long range connections in PFC are enhanced and this enhancement might be associated with the stronger information flow from right PFC to left PFC.

SummaryChronic stress could trigger maladaptive changes associated with stress-related mental disorders; however, the underlying mechanisms remain elusive. In this study, we found that exposing juvenile male rats to repeated stress significantly impaired the temporal order recognition memory, a cognitive process controlled by the prefrontal cortex (PFC). Concomitantly, significantly reduced AMPAR- and NMDAR-mediated synaptic transmission and glutamate receptor expression were found in PFC pyramidal neurons from repeatedly stressed animals. All these effects relied on activation of glucocorticoid receptors and the subsequent enhancement of ubiquitin/proteasome-mediated degradation of GluR1 and NR1 subunits, which was controlled by the E3 ubiquitin ligase Nedd4-1 and Fbx2, respectively. Inhibition of proteasomes or knockdown of Nedd4-1 and Fbx2 in PFC prevented the loss of glutamatergic responses and recognition memory in stressed animals. Our results suggest that repeated stress dampens PFC glutamatergic transmission by facilitating glutamate receptor turnover, which causes the detrimental effect on PFC-dependent cognitive processes.

Homeostatic plasticity plays a critical role in the stability of neuronal activities. Here, with high-density hippocampal networks cultured on multi-electrode arrays (MEAs), the transformation of spontaneous neuronal firing patterns induced by 1microM tetrodotoxin was clarified. Once tetrodotoxin was washed out after a 4-h treatment, spontaneous activities rose significantly with spike rate increasing approximately three times, and synchronized burst oscillations appeared throughout the network, with the cross-correlation coefficient between the active sites rising from 0.06+/-0.03 to 0.27+/-0.05. The long-term recording showed that the oscillations lasted for more than 4h before the network recovered. These results suggest that short-term treatment by tetrodotoxin may induce the homeostatically enhanced neuronal excitability, and that the spontaneous synchronized oscillations should be an indicator of homeostatic plasticity in cultured neuronal network. Furthermore, the non-invasive and long-term recording with MEAs as a novel sensing system is identified to be appropriate for pharmacological investigations of neuronal plasticity at the network level.

Cognitive dysfunction may result from abnormality of ionotropic glutamate receptors. Although various forms of synaptic plasticity in learning that rely on altering of glutamate receptors have been considered, the evidence is insufficient from an informatics view. Dynamics could reflect neuroinformatics encoding, including temporal pattern encoding, spatial pattern encoding, and energy distribution. Discovering informatics encoding is fundamental and crucial to understanding the working principle of the neural system. In this article, we analyzed the dynamic characteristics of response activities during learning training in cultured hippocampal networks under normal and abnormal conditions of ionotropic glutamate receptors, respectively. The rate, which is one of the temporal configurations, was decreased markedly by inhibition of α-amino-3-hydroxy-5-methylisoxazole-4-proprionic acid (AMPA) receptors. Moreover, the energy distribution in different characteristic frequencies was changed markedly by inhibition of AMPA receptors. Spatial configurations, including regularization, correlation, and synchrony, were changed significantly by inhibition of N-methyl-d-aspartate receptors. These results suggest that temporal pattern encoding and energy distribution of response activities in cultured hippocampal neuronal networks during learning training are modulated by AMPA receptors, whereas spatial pattern encoding of response activities is modulated by N-methyl-d-aspartate receptors.

Spontaneous synchronized bursts seem to play a key role in brain functions such as learning and memory. Still controversial is the characterization of spontaneous synchronized bursts in neuronal networks after learning training, whether depression or promotion. By taking advantages of the main features of the microelectrode array (MEA) technology (i.e. multisite recordings, stable and long-term coupling with the biological preparation), we analyzed changes of spontaneous synchronized bursts in cultured hippocampal neuronal networks after learning training. And for this purpose, a learning model at networking level on MEA system was constructed, and analysis of spontaneous synchronized burst activity modulation was presented. Preliminary results show that, the number of burst was increased by 154%, burst duration was increased by 35%, and the number of spikes per burst was increased by 124%, while interburst interval decreased by 44% with learning. In particular, correlation and synchrony of neuronal activities in networks were enhanced by 51% and 36%, respectively, with learning. In contrast, dynamic properties of neuronal networks were not changed much when the network was under "non-learning" condition. These results indicate that firing, association and synchrony of spontaneous bursts in neuronal networks were promoted by learning. Furthermore, from these observations, we are encouraged to think of a more engineered system based on in vitro hippocampal neurons, as a novel sensitive system for electrophysiological evaluations.

Spontaneous neuronal activity plays an important role in development. However, the mechanism that underlies the long-term spontaneous developmental change of cultured neuronal networks in vitro is not well understood. To investigate the contribution of inhibitory and excitatory connections to the development of neuronal networks, dissociated neurons from an embryonic rat hippocampal formation were cultured on a multi-electrode array plate and spontaneous activities were recorded by multi-channel system. These spontaneous activities were compared to bicuculline-induced firings, which were recorded by 60 electrodes simultaneously from 1 to 14 weeks in vitro (WIV). The phenomena showed that the spontaneous firing activities changed from an initial pattern of synchronized bursts to a later pattern of high frequency random spikes. The bicuculline-induced firing activities transformed from a pattern of synchronized bursts throughout all active sites in 3 WIV, to a pattern of local synchronized or random spikes appearing in the intervals of synchronized bursts after 11 WIV, while the firing rate hardly changed. Kynurenic acid, a broad-spectrum glutamate receptor antagonist, blocked all activities while CNQX inhibited only the local synchronized or random spikes. These suggest that the inhibitory connection was age-dependent degraded in vitro and the developmental spontaneous firing pattern was built by the homeostatic balance of the excitatory-inhibitory connection networks. Long-term cultures on MEA provided a useful tool to measure the relationship between spontaneous developmental change and pharmacological influence in vitro.

There are some unsolvable fundamental questions, such as cell type classification, neural circuit tracing and neurovascular coupling, though great progresses are being made in neuroscience. Because of the structural features of neurons and neural circuits, the solution of these questions needs us to break through the current technology of neuroanatomy for acquiring the exactly fine morphology of neuron and vessels and tracing long-distant circuit at axonal resolution in the whole brain of mammals. Combined with fast-developing labeling techniques, efficient whole-brain optical imaging technology emerging at the right moment presents a huge potential in the structure and function research of specific-function neuron and neural circuit. In this review, we summarize brain-wide optical tomography techniques, review the progress on visible brain neuronal/vascular networks benefit from these novel techniques, and prospect the future technical development.

We will report 3D brainwide neuron morphology, neural circuit and vascular network visualization protocol for a whole mouse brain, including developed staining methods and one micron voxel resolution using micro-optical sectioning tomography (MOST).

Self-organized criticality emerged in neural activity is one of the key concepts to describe the formation and the function of developing neuronal networks. The relationship between critical dynamics and neural development is both theoretically and experimentally appealing. However, whereas it is well-known that cortical networks exhibit a rich repertoire of activity patterns at different stages during in vitro maturation, dynamical activity patterns through the entire neural development still remains unclear. Here we show that a series of metastable network states emerged in the developing and “aging” process of hippocampal networks cultured from dissociated rat neurons. The unidirectional sequence of state transitions could be only observed in networks showing power-law scaling of distributed neuronal avalanches. Our data suggest that self-organized criticality may guide spontaneous activity into a sequential succession of homeostatically-regulated transient patterns during development, which may help to predict the tendency of neural development at early ages in the future.

Slow oscillations in the hippocampus are correlated with memory consolidation and brain diseases. The characteristic firings of the hippocampal network in vitro are still poorly understood. Here, spontaneous oscillations (~0.004 Hz) were found in high-density hippocampal networks by multi-electrode arrays after 30 days in vitro. This kind of spontaneous activity was characterized by periodic synchronized superbursts, which persisted for approximately 60 s at long intervals. Additionally, 1-Hz stimulation (duration

The spontaneous synchronous activity is a common behavior in a developing brain and plays a critical role in establishing appropriate connections and certain clinical diseases. Therefore, the investigation of the synchronous firing is important for understanding the formation of functional circuits and their implications in the network plasticity. In a limited period of time during development, the neuronal networks show synchronous activities, which occur simultaneously on a large amount of cells and varies wildly among different preparations. In this study, the spontaneous synchronous bursts are observed during the development of cultured neuron networks on multi-electrode array. The initiating site of a round of spontaneous synchronous burst, estimated from the relative delays of onsets of activities between electrodes, distributed randomly from each burst, while our statistical results confirmed that the positions of such initiating sites are stable. By calculating the cross-correlation function of the spike trains recorded from different electrodes simultaneously, the spreading mode and the spreading topography of the synchronized bursting activity were described. To access the changes in firing patterns in disinhibited cultured networks, the spontaneous activities were compared with the firings when the network exposed to bicuculline, the blocker of GABAA receptor. The results showed that the generation of synchronous bursts in cultured neuron networks is governed by the level of spontaneous activities and by the balance between excitation and inhibition circuits.

GABAergic neurons, considered as the chief inhibitory neurons in the mammalian hippocampus, comprise a myriad of diverse neuronal subtypes. However, due to difficulties in the identification of GABAergic neurons in living tissues, characteristics of GABAergic neurons involved in the balance controll of neural network remains unclear. To facilitate the study of these neurons, neurons that express enhanced green flourscent protein (GFP) from neonatal transgenic mice were isolated and purified to examine its biological characters in vitro. The multi-electrode array and confocal microscopy is ideal for repeatly recording over hours or days. Combining neural activities recorded simultaneously from 60 channels in a multi-electrode array with morphological changes imaged with the confocal microscopy technique, we have succeeded in analysing the electrophysiological and morphological characteristics of GABAergic neurons in cultured hippocampal networks. The result showed that morphological changes might relate with signal processing in neural network. These suggest that the cultured transgenic neural networks of hippocampal neurons grown on microelectrode arrays could be a useful model system for future studies on properties of neuron ensembles and roles of GABAergic neurons in living neuronal network.

Burst as a unit of information coding is widely investigated in the developing central nervous system. However the mechanism underling the bursts generate and disappear is unclear at present. Neurons cultured on the multi-electrode arrays, are spontaneously active, and show complex pattern with random spikes and bursts firing. With long-term recording, the course of bursts generation and disappearance was detected. The results showed that the firing pattern could transform from random spikes to bursts firing. In the beginning, the random spikes rate decreased, accompanied with bursts occurred once in a while. It appeared both single spikes and bursts at the same time. After that, the random spikes disappeared. Spontaneous activity displayed a regular occurrence of bursts with shorter interspike interval. During such bursts the firing rate at the active sites was increased dramatically. After several seconds, firing rate decreased, interburst interval extended, accompanied with the occurrence of random spikes, opposite to the beginning. At last, bursts disappeared and the networks just fired in random spikes. The observation showed that the complex electrophysiological activities of the cultured neuronal networks could implicate the spontaneous generation of burst firing. Understanding how bursts generate and disappear might be significant for deeply investigating the function and mechanism of bursts information coding.

In this paper, we study the modulation of low frequency closed-loop electric stimulation on spontaneous activity in cultured hippocampal neuronal networks. First, we plated monolayer cultures of hippocampal neurons from rat embryos (E18) on multi-electrode arrays and the experiments were performed in the networks from the second week to the sixth week continuously. During the experiments, we detected the spontaneous spikes of the networks firstly, and then stimulated the networks at low frequency (0.2 Hz or 1 Hz) stimulation respectively until a desired response was observed 20-80msec after a stimulus. The protocol was closed-loop. After that, we detected the spontaneous spikes of the networks. It is observed that the spontaneous activity in the developing networks is developing, which is oscillatory and periodic. Low frequency (0.2 Hz or 1 Hz) stimulation enhanced the spontaneous synchronous burst activity of the developing networks. These results implicated that activity-dependent mechanism in the modulation of plasticity of synaptic transmission in the cultured neuronal networks. Closed-loop stimulation will give a better view on the functional significance of networks activities. Besides, close-loop stimulation could set up the stimulus-reward system in the neuronal networks, which is of great benefit to the plasticity of synaptic transmission in the cultured neuronal networks.

Spontaneous firing play an important role in development of neuronal network. Activity-dependent modification of synaptic efficacy is widely recognized as a cellular basis of learning, memory, and development plasticity. Little is known of the activity-dependent modification of the synchronized spontaneous firing of the hippocampal networks. Long-term recording of spontaneous activity in cultured hippocampal neuronal networks was carried out using substrates containing multi-electrode array (MEA). Spontaneous uncorrelated firing appeared within a week and transformed progressively into synchronized pattern. During the development, these synchronized firings became into oscillation pattern and the synchronization has little change. By paired stimulation from adjacent electrodes in the network, the synchronized firing form a larger network burst. These results suggest that synchronized spontaneous spikes show the development of neuronal network and electronical stimulation could change the development.

Learning is crucial for brain functioning, especially for hippocampus. To understand the basis of physiological information processing of hippocampus, we developed the networks of purified-cultured rat hippocampus neurons on a planar 60-microelectrode array MEA. We describe how these networks are prepared, their structure, and the forms of spontaneous activity they exhibit. We stimulate the network at a low frequency until a desired predefined response is observed 50- 10 msec after a stimulus, at which point the stimulus is stopped for 5 min. Repeated cycles of this procedure ultimately lead to the desired response being directly elicited by the stimulus. These results demonstrate that there is learning in networks composed of purified-cultured hippocampus neurons.