Workflow and atlas system for brain-wide mapping of axonal connectivity in rat

The efferent connections of the orbitofrontal, posterior parietal, and insular cortex of the rat brain

The orbitofrontal, posterior parietal, and insular cortices are sites of higher-order cognitive processing implicated in a wide range of behaviours, including working memory, attention guiding, decision making, and spatial navigation. To better understand how these regions contribute to such functions, we need detailed knowledge about the underlying structural connectivity. Several tract-tracing studies have investigated specific aspects of orbitofrontal, posterior parietal and insular connectivity, but a digital resource for studying the cortical and subcortical projections from these areas in detail is not available. We here present a comprehensive collection of brightfield and fluorescence microscopic images of serial coronal sections from 49 rat brain tract-tracing experiments, in which discrete injections of the anterograde tracers biotinylated dextran amine and/or Phaseolus vulgaris leucoagglutinin were placed in the orbitofrontal, parietal, or insular cortex. The images are spatially registered to the Waxholm Space Rat brain atlas. The image collection, with corresponding reference atlas maps, is suitable as a reference framework for investigating the brain-wide efferent connectivity of these cortical association areas.

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.

Myelination of Callosal Axons Is Hampered by Early and Late Forelimb Amputation in Rats

Deafferentation is an important determinant of plastic changes in the CNS, which consists of a loss of inputs from the body periphery or from the CNS itself. Although cortical reorganization has been well documented, white matter plasticity was less explored. Our goal was to investigate microstructural interhemispheric connectivity changes in early and late amputated rats. For that purpose, we employed diffusion-weighted magnetic resonance imaging, as well as Western blotting, immunohistochemistry, and electron microscopy of sections of the white matter tracts to analyze the microstructural changes in the corticospinal tract and in the corpus callosum (CC) sector that contains somatosensory fibers integrating cortical areas representing the forelimbs and compare differences in rats undergoing forelimb amputation as neonates, with those amputated as adults. Results showed that early amputation induced decreased fractional anisotropy values and reduction of total myelin amount in the cerebral peduncle contralateral to the amputation. Both early and late forelimb amputations induced decreased myelination of callosal fibers. While early amputation affected myelination of thinner axons, late amputation disrupted axons of all calibers. Since the CC provides a modulation of inhibition and excitation between the hemispheres, we suggest that the demyelination observed among callosal fibers may misbalance this modulation.

Connection strength of the macaque connectome augments topological and functional network attributes

Mammalian brains constitute complex organized networks of neural projections. On top of their binary topological organization, the strength (or weight) of these neural projections can be highly variable across connections and is thus likely of additional importance to the overall topological and functional organization of the network. Here we investigated the specific distribution pattern of connection strength in the macaque connectome. We performed weighted and binary network analysis on the cortico-cortical connectivity of the macaque provided by the unique tract-tracing dataset of Markov and colleagues (2014) and observed in both analyses a small-world, modular and rich club organization. Moreover, connectivity strength showed a distribution augmenting the architecture identified in the binary network version by enhancing both local network clustering and the central infrastructure for global topological communication and integration. Functional consequences of this topological distribution were further examined using the Kuramoto model for simulating interactions between brain regions and showed that the connectivity strength distribution across connections enhances synchronization within modules and between rich club hubs. Together, our results suggest that neural pathway strength promotes topological properties in the macaque connectome for local processing and global network integration.

Spatial registration of serial microscopic brain images to three-dimensional reference atlases with the QuickNII tool

Modern high throughput brain wide profiling techniques for cells and their morphology, connectivity, and other properties, make the use of reference atlases with 3D coordinate frameworks essential. However, anatomical location of observations made in microscopic sectional images from rodent brains is typically determined by comparison with 2D anatomical reference atlases. A major challenge in this regard is that microscopic sections often are cut with orientations deviating from the standard planes used in the reference atlases, resulting in inaccuracies and a need for tedious correction steps. Overall, efficient tools for registration of large series of section images to reference atlases are currently not widely available. Here we present QuickNII, a stand-alone software tool for semi-automated affine spatial registration of sectional image data to a 3D reference atlas coordinate framework. A key feature in the tool is the capability to generate user defined cut planes through the reference atlas, matching the orientation of the cut plane of the sectional image data. The reference atlas is transformed to match anatomical landmarks in the corresponding experimental images. In this way, the spatial relationship between experimental image and atlas is defined, without introducing distortions in the original experimental images. Following anchoring of a limited number of sections containing key landmarks, transformations are propagated across the entire series of sectional images to reduce the amount of manual steps required. By having coordinates assigned to the experimental images, further analysis of the distribution of features extracted from the images is greatly facilitated.

Navigating the Murine Brain: Toward Best Practices for Determining and Documenting Neuroanatomical Locations in Experimental Studies

In experimental neuroscientific research, anatomical location is a key attribute of experimental observations and critical for interpretation of results, replication of findings, and comparison of data across studies. With steadily rising numbers of publications reporting basic experimental results, there is an increasing need for integration and synthesis of data. Since comparison of data relies on consistently defined anatomical locations, it is a major concern that practices and precision in the reporting of location of observations from different types of experimental studies seem to vary considerably. To elucidate and possibly meet this challenge, we have evaluated and compared current practices for interpreting and documenting the anatomical location of measurements acquired from murine brains with different experimental methods. Our observations show substantial differences in approach, interpretation and reproducibility of anatomical locations among reports of different categories of experimental research, and strongly indicate that ambiguous reports of anatomical location can be attributed to missing descriptions. Based on these findings, we suggest a set of minimum requirements for documentation of anatomical location in experimental murine brain research. We furthermore demonstrate how these requirements have been applied in the EU Human Brain Project to optimize workflows for integration of heterogeneous data in common reference atlases. We propose broad adoption of some straightforward steps for improving the precision of location metadata and thereby facilitating interpretation, reuse and integration of data.

Diverse Long-Range Axonal Projections of Excitatory Layer 2/3 Neurons in Mouse Barrel Cortex

Excitatory projection neurons of the neocortex are thought to play important roles in perceptual and cognitive functions of the brain by directly connecting diverse cortical and subcortical areas. However, many aspects of the anatomical organization of these inter-areal connections are unknown. Here, we studied long-range axonal projections of excitatory layer 2/3 neurons with cell bodies located in mouse primary somatosensory barrel cortex (wS1). As a population, these neurons densely projected to secondary whisker somatosensory cortex (wS2) and primary/secondary whisker motor cortex (wM1/2), with additional axon in the dysgranular zone surrounding the barrel field, perirhinal temporal association cortex and striatum. In three-dimensional reconstructions of 6 individual wS2-projecting neurons and 9 individual wM1/2-projecting neurons, we found that both classes of neurons had extensive local axon in layers 2/3 and 5 of wS1. Neurons projecting to wS2 did not send axon to wM1/2, whereas a small subset of wM1/2-projecting neurons had relatively weak projections to wS2. A small fraction of projection neurons solely targeted wS2 or wM1/2. However, axon collaterals from wS2-projecting and wM1/2-projecting neurons were typically also found in subsets of various additional areas, including the dysgranular zone, perirhinal temporal association cortex and striatum. Our data suggest extensive diversity in the axonal targets selected by individual nearby cortical long-range projection neurons with somata located in layer 2/3 of wS1.

Brain-Wide Mapping of Axonal Connections: Workflow for Automated Detection and Spatial Analysis of Labeling in Microscopic Sections

Axonal tracing techniques are powerful tools for exploring the structural organization of neuronal connections. Tracers such as biotinylated dextran amine (BDA) and Phaseolus vulgaris leucoagglutinin (Pha-L) allow brain-wide mapping of connections through analysis of large series of histological section images. We present a workflow for efficient collection and analysis of tract-tracing datasets with a focus on newly developed modules for image processing and assignment of anatomical location to tracing data. New functionality includes automatic detection of neuronal labeling in large image series, alignment of images to a volumetric brain atlas, and analytical tools for measuring the position and extent of labeling. To evaluate the workflow, we used high-resolution microscopic images from axonal tracing experiments in which different parts of the rat primary somatosensory cortex had been injected with BDA or Pha-L. Parameters from a set of representative images were used to automate detection of labeling in image series covering the entire brain, resulting in binary maps of the distribution of labeling. For high to medium labeling densities, automatic detection was found to provide reliable results when compared to manual analysis, whereas weak labeling required manual curation for optimal detection. To identify brain regions corresponding to labeled areas, section images were aligned to the Waxholm Space (WHS) atlas of the Sprague Dawley rat brain (v2) by custom-angle slicing of the MRI template to match individual sections. Based on the alignment, WHS coordinates were obtained for labeled elements and transformed to stereotaxic coordinates. The new workflow modules increase the efficiency and reliability of labeling detection in large series of images from histological sections, and enable anchoring to anatomical atlases for further spatial analysis and comparison with other data.

Three-Dimensional Histology Volume Reconstruction of Axonal Tract Tracing Data: Exploring Topographical Organization in Subcortical Projections from Rat Barrel Cortex

Topographical organization is a hallmark of the mammalian brain, and the spatial organization of axonal connections in different brain regions provides a structural framework accommodating specific patterns of neural activity. The presence, amount, and spatial distribution of axonal connections are typically studied in tract tracing experiments in which axons or neurons are labeled and examined in histological sections. Three-dimensional (3-D) reconstruction techniques are used to achieve more complete visualization and improved understanding of complex topographical relationships. 3-D reconstruction approaches based on manually or semi-automatically recorded spatial points representing axonal labeling have been successfully applied for investigation of smaller brain regions, but are not practically feasible for whole-brain analysis of multiple regions. We here reconstruct serial histological images from four whole brains (originally acquired for conventional microscopic analysis) into volumetric images that are spatially registered to a 3-D atlas template. The aims were firstly to evaluate the quality of the 3-D reconstructions and the usefulness of the approach, and secondly to investigate axonal projection patterns and topographical organization in rat corticostriatal and corticothalamic pathways. We demonstrate that even with the limitations of the original routine histological material, the 3-D reconstructed volumetric images allow efficient visualization of tracer injection sites and axonal labeling, facilitating detection of spatial distributions and across-case comparisons. Our results further show that clusters of S1 corticostriatal and corticothalamic projections are distributed within narrow, elongated or spherical subspaces extending across the entire striatum / thalamus. We conclude that histology volume reconstructions facilitate mapping of spatial distribution patterns and topographical organization. The reconstructed image volumes are shared via the Rodent Brain Workbench (www.rbwb.org).

Confocal stereology: an efficient tool for measurement of microscopic structures

Quantitative measurements of geometric forms or counting of objects in microscopic specimens is an essential tool in studies of microstructure. Confocal stereology represents a contemporary approach to the evaluation of microscopic structures by using a combination of stereological methods and confocal microscopy. 3-D images acquired by confocal microscopy can be used for the estimation of geometrical characteristics of microscopic structures by stereological methods, based on the evaluation of optical sections within a thick slice and using computer-generated virtual test probes. Such methods can be used for estimating volume, number, surface area and length using relevant spatial probes, which are generated by specific software. The interactions of the probes with the structure under study are interactively evaluated. An overview of the methods of confocal stereology developed during the past 30 years is presented. Their advantages and pitfalls in comparison with other methods for measurement of geometrical characteristics of microscopic structures are discussed.

Anatomical landmarks for registration of experimental image data to volumetric rodent brain atlasing templates

BACKGROUND

Assignment of anatomical reference is a key step in integration of the rapidly expanding collection of rodent brain data. Landmark-based registration facilitates spatial anchoring of diverse types of data not suitable for automated methods operating on voxel-based image information.

NEW TOOL

Here we propose a standardized set of anatomical landmarks for registration of whole brain imaging datasets from the mouse and rat brain, and in particular for integration of experimental image data in Waxholm Space (WHS).

RESULTS

Sixteen internal landmarks of the C57BL/6J mouse brain have been reliably identified: by different individuals, independent of their experience in anatomy; across different MRI contrasts (T1, T2, T2(*)) and other modalities (Nissl histology and block-face anatomy); in different specimens; in different slice acquisition angles; and in different image resolutions. We present a registration example between T1-weighted MRI and the mouse WHS template using these landmarks and reaching fairly high accuracy. Landmark positions identified in the mouse WHS template are shared through the Scalable Brain Atlas, accompanied by graphical and textual guidelines for locating each landmark. We identified 14 of the 16 landmarks in the WHS template for the Sprague Dawley rat.

COMPARISON WITH EXISTING METHODS

This landmark set can withstand substantial differences in acquisition angle, imaging modality, and is less vulnerable to subjectivity.

CONCLUSIONS

This facilitates registration of multimodal 3D brain data to standard coordinate spaces for mouse and rat brain taking a step toward the creation of a common rodent reference system; raising data sharing to a qualitatively higher level.

Cyberinfrastructure for the digital brain: spatial standards for integrating rodent brain atlases

Biomedical research entails capture and analysis of massive data volumes and new discoveries arise from data-integration and mining. This is only possible if data can be mapped onto a common framework such as the genome for genomic data. In neuroscience, the framework is intrinsically spatial and based on a number of paper atlases. This cannot meet today's data-intensive analysis and integration challenges. A scalable and extensible software infrastructure that is standards based but open for novel data and resources, is required for integrating information such as signal distributions, gene-expression, neuronal connectivity, electrophysiology, anatomy, and developmental processes. Therefore, the International Neuroinformatics Coordinating Facility (INCF) initiated the development of a spatial framework for neuroscience data integration with an associated Digital Atlasing Infrastructure (DAI). A prototype implementation of this infrastructure for the rodent brain is reported here. The infrastructure is based on a collection of reference spaces to which data is mapped at the required resolution, such as the Waxholm Space (WHS), a 3D reconstruction of the brain generated using high-resolution, multi-channel microMRI. The core standards of the digital atlasing service-oriented infrastructure include Waxholm Markup Language (WaxML): XML schema expressing a uniform information model for key elements such as coordinate systems, transformations, points of interest (POI)s, labels, and annotations; and Atlas Web Services: interfaces for querying and updating atlas data. The services return WaxML-encoded documents with information about capabilities, spatial reference systems (SRSs) and structures, and execute coordinate transformations and POI-based requests. Key elements of INCF-DAI cyberinfrastructure have been prototyped for both mouse and rat brain atlas sources, including the Allen Mouse Brain Atlas, UCSD Cell-Centered Database, and Edinburgh Mouse Atlas Project.

The physiological plausibility of time-varying Granger-causal modeling: normalization and weighting by spectral power

Time-varying connectivity methods are increasingly used to study directed interactions between brain regions from electrophysiological signals. These methods often show good results in simulated data but it is unclear to what extent connectivity results obtained from real data are physiologically plausible. Here we introduce a benchmark approach using multichannel somatosensory evoked potentials (SEPs) measured across rat cortex, where the structural and functional connectivity is relatively simple and well-understood. Rat SEPs to whisker stimulation are exclusively initiated by contralateral primary sensory cortex (S1), at known latencies, and with activity spread from S1 to specific cortical regions. This allows for a comparison of time-varying connectivity measures according to fixed criteria. We thus evaluated the performance of time-varying Partial Directed Coherence (PDC) and the Directed Transfer Function (DTF), comparing row- and column-wise normalization and the effect of weighting by the power spectral density (PSD). The benchmark approach revealed clear differences between methods in terms of physiological plausibility, effect size and temporal resolution. The results provide a validation of time-varying directed connectivity methods in an animal model and suggest a driving role for ipsilateral S1 in the later part of the SEP. The benchmark SEP dataset is made freely available.

Brain-wide map of efferent projections from rat barrel cortex

The somatotopically organized whisker barrel field of the rat primary somatosensory (S1) cortex is a commonly used model system for anatomical and physiological investigations of sensory processing. The neural connections of the barrel cortex have been extensively mapped. But most investigations have focused on connections to limited regions of the brain, and overviews in the literature of the connections across the brain thus build on a range of material from different laboratories, presented in numerous publications. Furthermore, given the limitations of the conventional journal article format, analyses and interpretations are hampered by lack of access to the underlying experimental data. New opportunities for analyses have emerged with the recent release of an online resource of experimental data consisting of collections of high-resolution images from 6 experiments in which anterograde tracers were injected in S1 whisker or forelimb representations. Building on this material, we have conducted a detailed analysis of the brain wide distribution of the efferent projections of the rat barrel cortex. We compare our findings with the available literature and reports accumulated in the Brain Architecture Management System (BAMS₂) database. We report well-known and less known intracortical and subcortical projections of the barrel cortex, as well as distinct differences between S1 whisker and forelimb related projections. Our results correspond well with recently published overviews, but provide additional information about relative differences among S1 projection targets. Our approach demonstrates how collections of shared experimental image data are suitable for brain-wide analysis and interpretation of connectivity mapping data.

Blood capillary length estimation from three-dimensional microscopic data by image analysis and stereology

Studies of the capillary bed characterized by its length or length density are relevant in many biomedical studies. A reliable assessment of capillary length from two-dimensional (2D), thin histological sections is a rather difficult task as it requires physical cutting of such sections in randomized directions. This is often technically demanding, inefficient, or outright impossible. However, if 3D image data of the microscopic structure under investigation are available, methods of length estimation that do not require randomized physical cutting of sections may be applied. Two different rat brain regions were optically sliced by confocal microscopy and resulting 3D images processed by three types of capillary length estimation methods: (1) stereological methods based on a computer generation of isotropic uniform random virtual test probes in 3D, either in the form of spatial grids of virtual "slicer" planes or spherical probes; (2) automatic method employing a digital version of the Crofton relations using the Euler characteristic of planar sections of the binary image; and (3) interactive "tracer" method for length measurement based on a manual delineation in 3D of the axes of capillary segments. The presented methods were compared in terms of their practical applicability, efficiency, and precision.

Optogenetic approaches for functional mouse brain mapping

To better understand the connectivity of the brain, it is important to map both structural and functional connections between neurons and cortical regions. In recent years, a set of optogenetic tools have been developed that permit selective manipulation and investigation of neural systems. These tools have enabled the mapping of functional connections between stimulated cortical targets and other brain regions. Advantages of the approach include the ability to arbitrarily stimulate brain regions that express opsins, allowing for brain mapping independent of behavior or sensory processing. The ability of opsins to be rapidly and locally activated allows for investigation of connectivity with spatial resolution on the order of single neurons and temporal resolution on the order of milliseconds. Optogenetic methods for functional mapping have been applied in experiments ranging from in vitro investigation of microcircuits, to in vivo probing of inter-regional cortical connections, to examination of global connections within the whole brain. We review recently developed functional mapping methods that use optogenetic single-point stimulation in the rodent brain and employ cellular electrophysiology, evoked motor movements, voltage sensitive dyes (VSDs), calcium indicators, or functional magnetic resonance imaging (fMRI) to assess activity. In particular we highlight results using red-shifted organic VSDs that permit high temporal resolution imaging in a manner spectrally separated from Channelrhodopsin-2 (ChR2) activation. VSD maps stimulated by ChR2 were dependent on intracortical synaptic activity and were able to reflect circuits used for sensory processing. Although the methods reviewed are powerful, challenges remain with respect to finding approaches that permit selective high temporal resolution assessment of stimulated activity in animals that can be followed longitudinally.

CoCoMac 2.0 and the future of tract-tracing databases

The CoCoMac database contains the results of several hundred published axonal tract-tracing studies in the macaque monkey brain. The combined results are used for constructing the macaque macro-connectome. Here we discuss the redevelopment of CoCoMac and compare it to six connectome-related projects: two online resources that provide full access to raw tracing data in rodents, a connectome viewer for advanced 3D graphics, a partial but highly detailed rat connectome, a brain data management system that generates custom connectivity matrices, and a software package that covers the complete pipeline from connectivity data to large-scale brain simulations. The second edition of CoCoMac features many enhancements over the original. For example, a search wizard is provided for full access to all tables and their nested dependencies. Connectivity matrices can be computed on demand in a user-selected nomenclature. A new data entry system is available as a preview, and is to become a generic solution for community-driven data entry in manually collated databases. We conclude with the question whether neuronal tracing will remain the gold standard to uncover the wiring of brains, thereby highlighting developments in human connectome construction, tracer substances, polarized light imaging, and serial block-face scanning electron microscopy.

Combining collation and annotation efforts toward completion of the rat and mouse connectomes in BAMS

Many different independently published neuroanatomical parcellation schemes (brain maps, nomenclatures, or atlases) can exist for a particular species, although one scheme (a standard scheme) is typically chosen for mapping neuroanatomical data in a particular study. This is problematic for building connection matrices (connectomes) because the terms used to name structures in different parcellation schemes differ widely and interrelationships are seldom defined. Therefore, data sets cannot be compared across studies that have been mapped on different neuroanatomical atlases without a reliable translation method. Because resliceable 3D brain models for relating systematically and topographically different parcellation schemes are still in the first phases of development, it is necessary to rely on qualitative comparisons between regions and tracts that are either inserted directly by neuroanatomists or trained annotators, or are extracted or inferred by collators from the available literature. To address these challenges, we developed a publicly available neuroinformatics system, the Brain Architecture Knowledge Management System (BAMS; http://brancusi.usc.edu/bkms). The structure and functionality of BAMS is briefly reviewed here, as an exemplar for constructing interrelated connectomes at different levels of the mammalian central nervous system organization. Next, the latest version of BAMS rat macroconnectome is presented because it is significantly more populated with the number of inserted connectivity reports exceeding a benchmark value (50,000), and because it is based on a different classification scheme. Finally, we discuss a general methodology and strategy for producing global connection matrices, starting with rigorous mapping of data, then inserting and annotating it, and ending with online generation of large-scale connection matrices.

Probing tissue microstructure with restriction spectrum imaging: Histological and theoretical validation

Water diffusion magnetic resonance imaging (dMRI) is a powerful tool for studying biological tissue microarchitectures in vivo. Recently, there has been increased effort to develop quantitative dMRI methods to probe both length scale and orientation information in diffusion media. Diffusion spectrum imaging (DSI) is one such approach that aims to resolve such information based on the three-dimensional diffusion propagator at each voxel. However, in practice, only the orientation component of the propagator function is preserved when deriving the orientation distribution function. Here, we demonstrate how a straightforward extension of the linear spherical deconvolution (SD) model can be used to probe tissue orientation structures over a range (or "spectrum") of length scales with minimal assumptions on the underlying microarchitecture. Using high b-value Cartesian q-space data on a rat brain tissue sample, we demonstrate how this "restriction spectrum imaging" (RSI) model allows for separating the volume fraction and orientation distribution of hindered and restricted diffusion, which we argue stems primarily from diffusion in the extraneurite and intraneurite water compartment, respectively. Moreover, we demonstrate how empirical RSI estimates of the neurite orientation distribution and volume fraction capture important additional structure not afforded by traditional DSI or fixed-scale SD-like reconstructions, particularly in gray matter. We conclude that incorporating length scale information in geometric models of diffusion offers promise for advancing state-of-the-art dMRI methods beyond white matter into gray matter structures while allowing more detailed quantitative characterization of water compartmentalization and histoarchitecture of healthy and diseased tissue.