Neural systems language: a formal modeling language for the systematic description, unambiguous communication, and automated digital curation of neural connectivity

Network architecture of intrinsic connectivity in a mammalian spinal cord (the central nervous system's caudal sector)

The vertebrate spinal cord (SP) is the long, thin extension of the brain forming the central nervous system's caudal sector. Functionally, the SP directly mediates motor and somatic sensory interactions with most parts of the body except the face, and it is the preferred model for analyzing relatively simple reflex behaviors. Here, we analyze the organization of axonal connections between the 50 gray matter regions forming the bilaterally symmetric rat SP. The assembled dataset suggests that there are about 385 of a possible 2,450 connections between the 50 regions for a connection density of 15.7%. Multiresolution consensus cluster analysis reveals a hierarchy of structure-function subsystems in this neural network, with 4 subsystems at the top level and 12 at the bottom-level. The top-level subsystems include a) a bilateral subsystem related most clearly to somatic and autonomic motor functions and centered in the ventral horn and intermediate zone; b) a bilateral subsystem associated with general somatosensory functions and centered in the base, neck, and head of the dorsal horn; and c) a pair of unilateral, bilaterally symmetric subsystems associated with nociceptive information processing and occupying the apex of the dorsal horn. The intrinsic SP network displayed no hubs, rich club, or small-world attributes, which are common measures of global functionality. Advantages and limitations of our methodology are discussed in some detail. The present work is part of a comprehensive project to assemble and analyze the neurome of a mammalian nervous system and its interactions with the body.

Intrinsic circuitry of the rhombicbrain (central nervous system's intermediate sector) in a mammal

The rhombicbrain (rhombencephalon or intermediate sector) is the vertebrate central nervous system part between the forebrain-midbrain (rostral sector) and spinal cord (caudal sector), and it has three main divisions: pons, cerebellum, and medulla. Using a data-driven approach, here we examine intrinsic rhombicbrain (intrarhombicbrain) network architecture that in rat consists of 52,670 possible axonal connections between 230 gray matter regions (115 bilaterally symmetrical pairs). Our analysis indicates that only 8,089 (15.4%) of these connections exist. Multiresolution consensus cluster analysis yields a nested hierarchy model of rhombicbrain subsystems that at the top level are associated with 1) the cerebellum and vestibular nuclei, 2) orofacial-pharyngeal-visceral integration, and 3) auditory connections; the bottom level has 68 clusters, ranging in size from 2 to 11 regions. The model provides a basis for functional hypothesis development and interrogation. More granular network analyses performed on the intrinsic connectivity of individual and combined main rhombicbrain divisions (pons, cerebellum, medulla, pons + cerebellum, and pons + medulla) demonstrate the mutability of network architecture in response to the addition or subtraction of connections. Clear differences between the structure-function network architecture of the rhombicbrain and forebrain-midbrain are discussed, with a stark comparison provided by the subsystem and small-world organization of the cerebellar cortex and cerebral cortex. Future analysis of the connections within and between the forebrain-midbrain and rhombicbrain will provide a model of brain neural network architecture in a mammal.

Structure-function subsystem models of female and male forebrain networks integrating cognition, affect, behavior, and bodily functions

The forebrain is the first of three primary vertebrate brain subdivisions. Macrolevel network analysis in a mammal (rat) revealed that the 466 gray matter regions composing the right and left sides of the forebrain are interconnected by 35,738 axonal connections forming a large set of overlapping, hierarchically arranged subsystems. This hierarchy is bilaterally symmetrical and sexually dimorphic, and it was used to create a structure-function conceptual model of intraforebrain network organization. Two mirror image top-level subsystems are presumably the most fundamental ontogenetically and phylogenetically. They essentially form the right and left forebrain halves and are relatively weakly interconnected. Each top-level subsystem in turn has two second-level subsystems. A ventromedial subsystem includes the medial forebrain bundle, functionally coordinating instinctive survival behaviors with appropriate physiological responses and affect. This subsystem has 26/24 (female/male) lowest-level subsystems, all using a combination of glutamate and GABA as neurotransmitters. In contrast, a dorsolateral subsystem includes the lateral forebrain bundle, functionally mediating voluntary behavior and cognition. This subsystem has 20 lowest-level subsystems, and all but 4 use glutamate exclusively for their macroconnections; no forebrain subsystems are exclusively GABAergic. Bottom-up subsystem analysis is a powerful engine for generating testable hypotheses about mechanistic explanations of brain function, behavior, and mind based on underlying circuit organization. Targeted computational (virtual) lesioning of specific regions of interest associated with Alzheimer's disease, clinical depression, and other disorders may begin to clarify how the effects spread through the entire forebrain network model.

Cognitive and Non-Cognitive Outcomes Associated with Student Engagement in a Novel Brain Chemoarchitecture Mapping Course-Based Undergraduate Research Experience

Course-based undergraduate research experiences (CUREs) engage emerging scholars in the authentic process of scientific discovery, and foster their development of content knowledge, motivation, and persistence in the science, technology, engineering, and mathematics (STEM) disciplines. Importantly, authentic research courses simultaneously offer investigators unique access to an extended population of students who receive education and mentoring in conducting scientifically relevant investigations and who are thus able to contribute effort toward big-data projects. While this paradigm benefits fields in neuroscience, such as atlas-based brain mapping of nerve cells at the tissue level, there are few documented cases of such laboratory courses offered in the domain. Here, we describe a curriculum designed to address this deficit, evaluate the scientific merit of novel student-produced brain atlas maps of immunohistochemically-identified nerve cell populations for the rat brain, and assess shifts in science identity, attitudes, and science communication skills of students engaged in the introductory-level Brain Mapping and Connectomics (BM&C) CURE. BM&C students reported gains in research and science process skills following participation in the course. Furthermore, BM&C students experienced a greater sense of science identity, including a greater likelihood to discuss course activities with non-class members compared to their non-CURE counterparts. Importantly, evaluation of student-generated brain atlas maps indicated that the course enabled students to produce scientifically valid products and make new discoveries to advance the field of neuroanatomy. Together, these findings support the efficacy of the BM&C course in addressing the relatively esoteric demands of chemoarchitectural brain mapping.

The network architecture of rat intrinsic interbrain (diencephalic) macroconnections

The endbrain and interbrain form 2 great vertebrate forebrain divisions, and the interbrain is subdivided into the hypothalamus ventrally and thalamus dorsally. General organizing principles of intrainterbrain axonal circuitry were examined here at the level of gray matter regions using network analysis tools in a mammal with the most complete available dataset-before examining interbrain input-output relationships with other nervous system parts. The dataset was curated expertly from the neuroanatomical literature using experimental axonal pathway-tracing methods, and evidence from 74,242 connection reports indicates the existence of 10,836 macroconnections of the possible 49,062 macroconnections between the 222 gray matter regions forming the right and left halves of the interbrain. Two identical sets of 6 putative hubs were identified in the intrainterbrain network and form a continuous tissue mass in a part of the right and left medial hypothalamus associated functionally with physiological mechanisms controlling bodily functions. The intrainterbrain network shows only weak evidence of small-world attributes, rich club organization is absent, and multiresolution consensus cluster analysis indicates a solution with only 3 top-level subsystems or modules. In contrast, a previous analysis employing the same methodology to the significantly denser 244-node intraendbrain network revealed 2 identical sets of 13 hubs, small-world and rich club attributes, and 4 top-level subsystems. These differences in intrinsic network architecture across subdivisions suggest that intrinsic connections shape regional functional specialization to a varying extent, in part driven by differences in density and centrality, with extrinsic input-output connectivity playing a greater role in subdivisions that are sparser and less centralized.

The network organization of rat intrathalamic macroconnections and a comparison with other forebrain divisions

The thalamus is 1 of 4 major divisions of the forebrain and is usually subdivided into epithalamus, dorsal thalamus, and ventral thalamus. The 39 gray matter regions comprising the large dorsal thalamus project topographically to the cerebral cortex, whereas the much smaller epithalamus (2 regions) and ventral thalamus (5 regions) characteristically project subcortically. Before analyzing extrinsic inputs and outputs of the thalamus, here, the intrinsic connections among all 46 gray matter regions of the rat thalamus on each side of the brain were expertly collated and subjected to network analysis. Experimental axonal pathway-tracing evidence was found in the neuroanatomical literature for the presence or absence of 99% of 2,070 possible ipsilateral connections and 97% of 2,116 possible contralateral connections; the connection density of ipsilateral connections was 17%, and that of contralateral connections 5%. One hub, the reticular thalamic nucleus (of the ventral thalamus), was found in this network, whereas no high-degree rich club or clear small-world features were detected. The reticular thalamic nucleus was found to be primarily responsible for conferring the property of complete connectedness to the intrathalamic network in the sense that there is, at least, one path of finite length between any 2 regions or nodes in the network. Direct comparison with previous investigations using the same methodology shows that each division of the forebrain (cerebral cortex, cerebral nuclei, thalamus, hypothalamus) has distinct intrinsic network topological organization. A future goal is to analyze the network organization of connections within and among these 4 divisions of the forebrain.

Subsystem organization of axonal connections within and between the right and left cerebral cortex and cerebral nuclei (endbrain)

The endbrain (telencephalon) is at the rostral end of the central nervous system and is primarily responsible for supporting cognition and affect. Structurally, it consists of right and left cerebral hemispheres, each parceled into multiple cortical and nuclear gray matter regions. The global network organization of axonal macroconnections between the 244 regions forming the endbrain was analyzed with a multiresolution consensus clustering (MRCC) method that provides a hierarchical description of community clustering (modules or subsystems) within the network. Experimental evidence was collated from the neuroanatomical literature for the existence of 10,002 of a possible 59,292 connections within the network, and they cluster into four top-level subsystems and 60 bottom-level subsystems arranged in a 50-level hierarchy. Two top-level subsystems are bihemispheric: One deals with auditory and visual information, and the other corresponds broadly to the default mode network. The other two top-level subsystems are bilaterally symmetrical, and each deals broadly with somatic and visceral information. Because the entire endbrain connection matrix was assembled from multiple subconnectomes, it was easy to show that the status of a region as a connectivity hub is not absolute but, instead, depends on the size and coverage of its anatomical neighborhood. It was also shown numerically that creating an ultradense connection matrix by converting all "absent" connections to a "very weak" connection weight has virtually no effect on the clustering hierarchy. The next logical step in this project is to complete the forebrain connectome by adding the thalamus and hypothalamus (together, the interbrain) to the endbrain analysis.

Organizing principles for the cerebral cortex network of commissural and association connections

Cognition is supported by a network of axonal connections between gray matter regions within and between right and left cerebral cortex. Global organizing principles of this circuitry were examined with network analysis tools applied to monosynaptic association (within one side) and commissural (between sides) connections between all 77 cortical gray matter regions in each hemisphere of the rat brain. The analysis used 32,350 connection reports expertly collated from published pathway tracing experiments, and 5,394 connections of a possible 23,562 were identified, for a connection density of 23%-of which 20% (1,084) were commissural. Network community detection yielded a stable bihemispheric six-module solution, with an identical set in each hemisphere of three modules topographically forming a lateral core and medial shell arrangement of cortical regions. Functional correlations suggest the lateral module deals preferentially with environmental sensory-motor interactions and the ventromedial module deals preferentially with visceral control, affect, and short-term memory, whereas the dorsomedial module resembles the default mode network. Analysis of commissural connections revealed a set of unexpected rules to help generate hypotheses. Most notably, there is an order of magnitude more heterotopic than homotopic projections; all cortical regions send more association than commissural connections, and for each region, the latter are always a subset of the former; the number of association connections from each cortical region strongly correlates with the number of its commissural connections; and the module (dorsomedial) lying closest to the corpus callosum has the most complete set of commissural connections-and apparently the most complex function.

Network architecture of the cerebral nuclei (basal ganglia) association and commissural connectome

The cerebral nuclei form the ventral division of the cerebral hemisphere and are thought to play an important role in neural systems controlling somatic movement and motivation. Network analysis was used to define global architectural features of intrinsic cerebral nuclei circuitry in one hemisphere (association connections) and between hemispheres (commissural connections). The analysis was based on more than 4,000 reports of histologically defined axonal connections involving all 45 gray matter regions of the rat cerebral nuclei and revealed the existence of four asymmetrically interconnected modules. The modules form four topographically distinct longitudinal columns that only partly correspond to previous interpretations of cerebral nuclei structure-function organization. The network of connections within and between modules in one hemisphere or the other is quite dense (about 40% of all possible connections), whereas the network of connections between hemispheres is weak and sparse (only about 5% of all possible connections). Particularly highly interconnected regions (rich club and hubs within it) form a topologically continuous band extending through two of the modules. Connection path lengths among numerous pairs of regions, and among some of the network's modules, are relatively long, thus accounting for low global efficiency in network communication. These results provide a starting point for reexamining the connectional organization of the cerebral hemispheres as a whole (right and left cerebral cortex and cerebral nuclei together) and their relation to the rest of the nervous system.

Golgi: Interactive Online Brain Mapping

Golgi (http://www.usegolgi.com) is a prototype interactive brain map of the rat brain that helps researchers intuitively interact with neuroanatomy, connectomics, and cellular and chemical architecture. The flood of “-omic” data urges new ways to help researchers connect discrete findings to the larger context of the nervous system. Here we explore Golgi’s underlying reasoning and techniques and how our design decisions balance the constraints of building both a scientifically useful and usable tool. We demonstrate how Golgi can enhance connectomic literature searches with a case study investigating a thalamocortical circuit involving the Nucleus Accumbens and we explore Golgi’s potential and future directions for growth in systems neuroscience and connectomics.

Connections of the juxtaventromedial region of the lateral hypothalamic area in the male rat

Evolutionary conservation of the hypothalamus attests to its critical role in the control of fundamental behaviors. However, our knowledge of hypothalamic connections is incomplete, particularly for the lateral hypothalamic area (LHA). Here we present the results of neuronal pathway-tracing experiments to investigate connections of the LHA juxtaventromedial region, which is parceled into dorsal (LHAjvd) and ventral (LHAjvv) zones. Phaseolus vulgaris leucoagglutinin (PHAL, for outputs) and cholera toxin B subunit (CTB, for inputs) coinjections were targeted stereotaxically to the LHAjvd/v.

Results: LHAjvd/v connections overlapped highly but not uniformly. Major joint outputs included: Bed nuc. stria terminalis (BST), interfascicular nuc. (BSTif) and BST anteromedial area, rostral lateral septal (LSr)- and ventromedial hypothalamic (VMH) nuc., and periaqueductal gray. Prominent joint LHAjvd/v input sources included: BSTif, BST principal nuc., LSr, VMH, anterior hypothalamic-, ventral premammillary-, and medial amygdalar nuc., and hippocampal formation (HPF) field CA1. However, LHAjvd HPF retrograde labeling was markedly more abundant than from the LHAjvv; in the LSr this was reversed. Furthermore, robust LHAjvv (but not LHAjvd) targets included posterior- and basomedial amygdalar nuc., whereas the midbrain reticular nuc. received a dense input from the LHAjvd alone. Our analyses indicate the existence of about 500 LHAjvd and LHAjvv connections with about 200 distinct regions of the cerebral cortex, cerebral nuclei, and cerebrospinal trunk. Several highly LHAjvd/v-connected regions have a prominent role in reproductive behavior. These findings contrast with those from our previous pathway-tracing studies of other LHA medial and perifornical tier regions, with different connectional behavioral relations. The emerging picture is of a highly differentiated LHA with extensive and far-reaching connections that point to a role as a central coordinator of behavioral control.

Architecture of the cerebral cortical association connectome underlying cognition

Cognition presumably emerges from neural activity in the network of association connections between cortical regions that is modulated by inputs from sensory and state systems and directs voluntary behavior by outputs to the motor system. To reveal global architectural features of the cortical association connectome, network analysis was performed on >16,000 reports of histologically defined axonal connections between cortical regions in rat. The network analysis reveals an organization into four asymmetrically interconnected modules involving the entire cortex in a topographic and topologic core-shell arrangement. There is also a topographically continuous U-shaped band of cortical areas that are highly connected with each other as well as with the rest of the cortex extending through all four modules, with the temporal pole of this band (entorhinal area) having the most cortical association connections of all. These results provide a starting point for compiling a mammalian nervous system connectome that could ultimately reveal novel correlations between genome-wide association studies and connectome-wide association studies, leading to new insights into the cellular architecture supporting cognition.

Controlling feeding behavior by chemical or gene-directed targeting in the brain: what's so spatial about our methods?

Intracranial chemical injection (ICI) methods have been used to identify the locations in the brain where feeding behavior can be controlled acutely. Scientists conducting ICI studies often document their injection site locations, thereby leaving kernels of valuable location data for others to use to further characterize feeding control circuits. Unfortunately, this rich dataset has not yet been formally contextualized with other published neuroanatomical data. In particular, axonal tracing studies have delineated several neural circuits originating in the same areas where ICI injection feeding-control sites have been documented, but it remains unclear whether these circuits participate in feeding control. Comparing injection sites with other types of location data would require careful anatomical registration between the datasets. Here, a conceptual framework is presented for how such anatomical registration efforts can be performed. For example, by using a simple atlas alignment tool, a hypothalamic locus sensitive to the orexigenic effects of neuropeptide Y (NPY) can be aligned accurately with the locations of neurons labeled by anterograde tracers or those known to express NPY receptors or feeding-related peptides. This approach can also be applied to those intracranial "gene-directed" injection (IGI) methods (e.g., site-specific recombinase methods, RNA expression or interference, optogenetics, and pharmacosynthetics) that involve viral injections to targeted neuronal populations. Spatial alignment efforts can be accelerated if location data from ICI/IGI methods are mapped to stereotaxic brain atlases to allow powerful neuroinformatics tools to overlay different types of data in the same reference space. Atlas-based mapping will be critical for community-based sharing of location data for feeding control circuits, and will accelerate our understanding of structure-function relationships in the brain for mammalian models of obesity and metabolic disorders.