Where Do Core Thalamocortical Axons Terminate in Mammalian Neocortex When There Is No Cytoarchitecturally Distinct Layer 4?

Although the mammalian cerebral cortex is most often described as a hexalaminar structure, there are cortical areas (primary motor cortex) and species (elephants, cetaceans, and hippopotami), where a cytoarchitecturally indistinct, or absent, layer 4 is noted. Thalamocortical projections from the core, or first order, thalamic system terminate primarily in layers 4/inner 3. We explored the termination sites of core thalamocortical projections in cortical areas and in species where there is no cytoarchitecturally distinct layer 4 using the immunolocalization of vesicular glutamate transporter 2, a known marker of core thalamocortical axon terminals, in 31 mammal species spanning the eutherian radiation. Several variations from the canonical cortical column outline of layer 4 and core thalamocortical inputs were noted. In shrews/microchiropterans, layer 4 was present, but many core thalamocortical projections terminated in layer 1 in addition to layers 4 and inner 3. In primate primary visual cortex, the sublaminated layer 4 was associated with a specialized core thalamocortical projection pattern. In primate primary motor cortex, no cytoarchitecturally distinct layer 4 was evident and the core thalamocortical projections terminated throughout layer 3. In the African elephant, cetaceans, and river hippopotamus, no cytoarchitecturally distinct layer 4 was observed and core thalamocortical projections terminated primarily in inner layer 3 and less densely in outer layer 3. These findings are contextualized in terms of cortical processing, perception, and the evolutionary trajectory leading to an indistinct or absent cortical layer 4.

Comparative basolateral amygdala connectomics reveals dissociable single-neuron projection patterns to frontal cortex in macaques and mice.

Basolateral amygdala (BLA) is a key hub for affect in the brain,1,2,3 and dysfunction within this area contributes to a host of psychiatric disorders.4,5 BLA is extensively and reciprocally interconnected with frontal cortex,6,7,8,9 and some aspects of its function are evolutionarily conserved across rodents, anthropoid primates, and humans.10 Neuron density in BLA is substantially lower in primates compared to murine rodents,11 and frontal cortex (FC) is dramatically expanded in primates, particularly the more anterior granular and dysgranular areas.12,13,14 Yet, how these anatomical differences influence the projection patterns of single BLA neurons to frontal cortex across rodents and primates is unknown. Using a barcoded connectomic approach, we assessed the single BLA neuron connections to frontal cortex in mice and macaques. We found that BLA neurons are more likely to project to multiple distinct parts of FC in mice than in macaques. Further, while single BLA neuron projections to nucleus accumbens were similarly organized in mice and macaques, BLA-FC connections differed substantially. Notably, BLA connections to subcallosal anterior cingulate cortex (scACC) in macaques were least likely to branch to other medial frontal cortex areas compared to perigenual ACC (pgACC). This pattern of connections was reversed in the mouse homologues of these areas, infralimbic and prelimbic cortex (IL and PL), mirroring functional differences between rodents and non-human primates. Taken together, these results indicate that BLA connections to FC are not linearly scaled from mice to macaques and instead the organization of single-neuron BLA connections is distinct between these species.

Proteomic features of gray matter layers and superficial white matter of the rhesus monkey neocortex: comparison of prefrontal area 46 and occipital area 17.

In this novel large-scale multiplexed immunofluorescence study we comprehensively characterized and compared layer-specific proteomic features within regions of interest of the widely divergent dorsolateral prefrontal cortex (A46) and primary visual cortex (A17) of adult rhesus monkeys. Twenty-eight markers were imaged in rounds of sequential staining, and their spatial distribution precisely quantified within gray matter layers and superficial white matter. Cells were classified as neurons, astrocytes, oligodendrocytes, microglia, or endothelial cells. The distribution of fibers and blood vessels were assessed by quantification of staining intensity across regions of interest. This method revealed multivariate similarities and differences between layers and areas. Protein expression in neurons was the strongest determinant of both laminar and regional differences, whereas protein expression in glia was more important for intra-areal laminar distinctions. Among specific results, we observed a lower glia-to-neuron ratio in A17 than in A46 and the pan-neuronal markers HuD and NeuN were differentially distributed in both brain areas with a lower intensity of NeuN in layers 4 and 5 of A17 compared to A46 and other A17 layers. Astrocytes and oligodendrocytes exhibited distinct marker-specific laminar distributions that differed between regions; notably, there was a high proportion of ALDH1L1-expressing astrocytes and of oligodendrocyte markers in layer 4 of A17. The many nuanced differences in protein expression between layers and regions observed here highlight the need for direct assessment of proteins, in addition to RNA expression, and set the stage for future protein-focused studies of these and other brain regions in normal and pathological conditions.

Deep learning-based localization algorithms on fluorescence human brain 3D reconstruction: a comparative study using stereology as a reference.

3D reconstruction of human brain volumes at high resolution is now possible thanks to advancements in tissue clearing methods and fluorescence microscopy techniques. Analyzing the massive data produced with these approaches requires automatic methods able to perform fast and accurate cell counting and localization. Recent advances in deep learning have enabled the development of various tools for cell segmentation. However, accurate quantification of neurons in the human brain presents specific challenges, such as high pixel intensity variability, autofluorescence, non-specific fluorescence and very large size of data. In this paper, we provide a thorough empirical evaluation of three techniques based on deep learning (StarDist, CellPose and BCFind-v2, an updated version of BCFind) using a recently introduced three-dimensional stereological design as a reference for large-scale insights. As a representative problem in human brain analysis, we focus on a 4 -cm 3 portion of the Broca's area. We aim at helping users in selecting appropriate techniques depending on their research objectives. To this end, we compare methods along various dimensions of analysis, including correctness of the predicted density and localization, computational efficiency, and human annotation effort. Our results suggest that deep learning approaches are very effective, have a high throughput providing each cell 3D location, and obtain results comparable to the estimates of the adopted stereological design.

Integrated platform for multiscale molecular imaging and phenotyping of the human brain.

Understanding cellular architectures and their connectivity is essential for interrogating system function and dysfunction. However, we lack technologies for mapping the multiscale details of individual cells and their connectivity in the human organ-scale system. We developed a platform that simultaneously extracts spatial, molecular, morphological, and connectivity information of individual cells from the same human brain. The platform includes three core elements: a vibrating microtome for ultraprecision slicing of large-scale tissues without losing cellular connectivity (MEGAtome), a polymer hydrogel-based tissue processing technology for multiplexed multiscale imaging of human organ-scale tissues (mELAST), and a computational pipeline for reconstructing three-dimensional connectivity across multiple brain slabs (UNSLICE). We applied this platform for analyzing human Alzheimer's disease pathology at multiple scales and demonstrating scalable neural connectivity mapping in the human brain.

Comparative anatomy of the caudate nucleus in canids and felids: Associations with brain size, curvature, cross-sectional properties, and behavioral ecology.

The evolutionary history of canids and felids is marked by a deep time separation that has uniquely shaped their behavior and phenotype toward refined predatory abilities. The caudate nucleus is a subcortical brain structure associated with both motor control and cognitive, emotional, and executive functions. We used a combination of three-dimensional imaging, allometric scaling, and structural analyses to compare the size and shape characteristics of the caudate nucleus. The sample consisted of MRI scan data obtained from six canid species (Canis lupus lupus, Canis latrans, Chrysocyon brachyurus, Lycaon pictus, Vulpes vulpes, Vulpes zerda), two canid subspecies (Canis lupus familiaris, Canis lupus dingo), as well as three felids (Panthera tigris, Panthera uncia, Felis silvestris catus). Results revealed marked conservation in the scaling and shape attributes of the caudate nucleus across species, with only slight deviations. We hypothesize that observed differences in caudate nucleus size and structure for the domestic canids are reflective of enhanced cognitive and emotional pathways that possibly emerged during domestication.

Increased NLRP1 mRNA and Protein Expression Suggests Inflammasome Activation in the Dorsolateral Prefrontal and Medial Orbitofrontal Cortex in Schizophrenia.

Schizophrenia is a complex mental condition, with key symptoms marked for diagnosis including delusions, hallucinations, disorganized thinking, reduced emotional expression, and social dysfunction. In the context of major developmental hypotheses of schizophrenia, notably those concerning maternal immune activation and neuroinflammation, we studied NLRP1 expression and content in the postmortem brain tissue of 10 schizophrenia and 10 control subjects. In the medial orbitofrontal cortex (Brodmann's area 11/12) and dorsolateral prefrontal cortex (area 46) from both hemispheres of six schizophrenia subjects, the NLRP1 mRNA expression was significantly higher than in six control brains (p < 0.05). As the expression difference was highest for the medial orbitofrontal cortex in the right hemisphere, we assessed NLRP1-immunoreactive pyramidal neurons in layers III, V, and VI in the medial orbitofrontal cortex in the right hemisphere of seven schizophrenia and five control brains. Compared to controls, we quantified a significantly higher number of NLRP1-positive pyramidal neurons in the schizophrenia brains (p < 0.01), suggesting NLRP1 inflammasome activation in schizophrenia subjects. Layer III pyramidal neuron dysfunction aligns with working memory deficits, while impairments of pyramidal neurons in layers V and VI likely disrupt predictive processing. We propose NLRP1 inflammasome as a potential biomarker and therapeutic target in schizophrenia.

Glycine is a transmitter in the human and chimpanzee cochlear nuclei.

Introduction: Auditory information is relayed from the cochlea via the eighth cranial nerve to the dorsal and ventral cochlear nuclei (DCN, VCN). The organization, neurochemistry and circuitry of the cochlear nuclei (CN) have been studied in many species. It is well-established that glycine is an inhibitory transmitter in the CN of rodents and cats, with glycinergic cells in the DCN and VCN. There are, however, major differences in the laminar and cellular organization of the DCN between humans (and other primates) and rodents and cats. We therefore asked whether there might also be differences in glycinergic neurotransmission in the CN. Methods: We studied brainstem sections from humans, chimpanzees, and cats. We used antibodies to glycine receptors (GLYR) to identify neurons receiving glycinergic input, and antibodies to the neuronal glycine transporter (GLYT2) to immunolabel glycinergic axons and terminals. We also examined archival sections immunostained for calretinin (CR) and nonphosphorylated neurofilament protein (NPNFP) to try to locate the octopus cell area (OCA), a region in the VCN that rodents has minimal glycinergic input. Results: In humans and chimpanzees we found widespread immunolabel for glycine receptors in DCN and in the posterior (PVCN) and anterior (AVCN) divisions of the VCN. We found a parallel distribution of GLYT2-immunolabeled fibers and puncta. The data also suggest that, as in rodents, a region containing octopus cells in cats, humans and chimpanzees has little glycinergic input. Discussion: Our results show that glycine is a major transmitter in the human and chimpanzee CN, despite the species differences in DCN organization. The sources of the glycinergic input to the CN in humans and chimpanzees are not known.

Tempo and mode of gene expression evolution in the brain across primates.

Primate evolution has led to a remarkable diversity of behavioral specializations and pronounced brain size variation among species (Barton, 2012; DeCasien and Higham, 2019; Powell et al., 2017). Gene expression provides a promising opportunity for studying the molecular basis of brain evolution, but it has been explored in very few primate species to date (e.g. Khaitovich et al., 2005; Khrameeva et al., 2020; Ma et al., 2022; Somel et al., 2009). To understand the landscape of gene expression evolution across the primate lineage, we generated and analyzed RNA-seq data from four brain regions in an unprecedented eighteen species. Here, we show a remarkable level of variation in gene expression among hominid species, including humans and chimpanzees, despite their relatively recent divergence time from other primates. We found that individual genes display a wide range of expression dynamics across evolutionary time reflective of the diverse selection pressures acting on genes within primate brain tissue. Using our samples that represent a 190-fold difference in primate brain size, we identified genes with variation in expression most correlated with brain size. Our study extensively broadens the phylogenetic context of what is known about the molecular evolution of the brain across primates and identifies novel candidate genes for the study of genetic regulation of brain evolution.

The Neurovascular Unit as a Locus of Injury in Low-Level Blast-Induced Neurotrauma.

Blast-induced neurotrauma has received much attention over the past decade. Vascular injury occurs early following blast exposure. Indeed, in animal models that approximate human mild traumatic brain injury or subclinical blast exposure, vascular pathology can occur in the presence of a normal neuropil, suggesting that the vasculature is particularly vulnerable. Brain endothelial cells and their supporting glial and neuronal elements constitute a neurovascular unit (NVU). Blast injury disrupts gliovascular and neurovascular connections in addition to damaging endothelial cells, basal laminae, smooth muscle cells, and pericytes as well as causing extracellular matrix reorganization. Perivascular pathology becomes associated with phospho-tau accumulation and chronic perivascular inflammation. Disruption of the NVU should impact activity-dependent regulation of cerebral blood flow, blood-brain barrier permeability, and glymphatic flow. Here, we review work in an animal model of low-level blast injury that we have been studying for over a decade. We review work supporting the NVU as a locus of low-level blast injury. We integrate our findings with those from other laboratories studying similar models that collectively suggest that damage to astrocytes and other perivascular cells as well as chronic immune activation play a role in the persistent neurobehavioral changes that follow blast injury.

Betz cells of the primary motor cortex.

Betz cells, named in honor of Volodymyr Betz (1834-1894), who described them as "giant pyramids" in the primary motor cortex of primates and other mammalian species, are layer V extratelencephalic projection (ETP) neurons that directly innervate α-motoneurons of the brainstem and spinal cord. Despite their large volume and circumferential dendritic architecture, to date, no single molecular criterion has been established that unequivocally distinguishes adult Betz cells from other layer V ETP neurons. In primates, transcriptional signatures suggest the presence of at least two ETP neuron clusters that contain mature Betz cells; these are characterized by an abundance of axon guidance and oxidative phosphorylation transcripts. How neurodevelopmental programs drive the distinct positional and morphological features of Betz cells in humans remains unknown. Betz cells display a distinct biphasic firing pattern involving early cessation of firing followed by delayed sustained acceleration in spike frequency and magnitude. Few cell type-specific transcripts and electrophysiological characteristics are conserved between rodent layer V ETP neurons of the motor cortex and primate Betz cells. This has implications for the modeling of disorders that affect the motor cortex in humans, such as amyotrophic lateral sclerosis (ALS). Perhaps vulnerability to ALS is linked to the evolution of neural networks for fine motor control reflected in the distinct morphomolecular architecture of the human motor cortex, including Betz cells. Here, we discuss histological, molecular, and functional data concerning the position of Betz cells in the emerging taxonomy of neurons across diverse species and their role in neurological disorders.

Metabotropic Glutamate Receptor 2 Expression Is Chronically Elevated in Male Rats With Post-Traumatic Stress Disorder Related Behavioral Traits Following Repetitive Low-Level Blast Exposure.

Many military veterans who experienced blast-related traumatic brain injuries in the conflicts in Iraq and Afghanistan currently suffer from chronic cognitive and mental health problems that include depression and post-traumatic stress disorder (PTSD). Male rats exposed to repetitive low-level blast develop cognitive and PTSD-related behavioral traits that are present for more than 1 year after exposure. We previously reported that a group II metabotropic receptor (mGluR2/3) antagonist reversed blast-induced behavioral traits. In this report, we explored mGluR2/3 expression following blast exposure in male rats. Western blotting revealed that mGluR2 protein (but not mGluR3) was increased in all brain regions studied (anterior cortex, hippocampus, and amygdala) at 43 or 52 weeks after blast exposure but not at 2 weeks or 6 weeks. mGluR2 RNA was elevated at 52 weeks while mGluR3 was not. Immunohistochemical staining revealed no changes in the principally presynaptic localization of mGluR2 by blast exposure. Administering the mGluR2/3 antagonist LY341495 after behavioral traits had emerged rapidly reversed blast-induced effects on novel object recognition and cued fear responses 10 months following blast exposure. These studies support alterations in mGluR2 receptors as a key pathophysiological event following blast exposure and provide further support for group II metabotropic receptors as therapeutic targets in the neurobehavioral effects that follow blast injury.

Human-specific features and developmental dynamics of the brain N-glycome.

Comparative "omics" studies have revealed unique aspects of human neurobiology, yet an evolutionary perspective of the brain N-glycome is lacking. We performed multiregional characterization of rat, macaque, chimpanzee, and human brain N-glycomes using chromatography and mass spectrometry and then integrated these data with complementary glycotranscriptomic data. We found that, in primates, the brain N-glycome has diverged more rapidly than the underlying transcriptomic framework, providing a means for rapidly generating additional interspecies diversity. Our data suggest that brain N-glycome evolution in hominids has been characterized by an overall increase in complexity coupled with a shift toward increased usage of α(2-6)-linked N-acetylneuraminic acid. Moreover, interspecies differences in the cell type expression pattern of key glycogenes were identified, including some human-specific differences, which may underpin this evolutionary divergence. Last, by comparing the prenatal and adult human brain N-glycomes, we uncovered region-specific neurodevelopmental pathways that lead to distinct spatial N-glycosylation profiles in the mature brain.

Segmentation of supragranular and infragranular layers in ultra-high resolution 7T ex vivo MRI of the human cerebral cortex.

Accurate labeling of specific layers in the human cerebral cortex is crucial for advancing our understanding of neurodevelopmental and neurodegenerative disorders. Leveraging recent advancements in ultra-high resolution ex vivo MRI, we present a novel semi-supervised segmentation model capable of identifying supragranular and infragranular layers in ex vivo MRI with unprecedented precision. On a dataset consisting of 17 whole-hemisphere ex vivo scans at 120 µm, we propose a multi-resolution U-Nets framework (MUS) that integrates global and local structural information, achieving reliable segmentation maps of the entire hemisphere, with Dice scores over 0.8 for supra- and infragranular layers. This enables surface modeling, atlas construction, anomaly detection in disease states, and cross-modality validation, while also paving the way for finer layer segmentation. Our approach offers a powerful tool for comprehensive neuroanatomical investigations and holds promise for advancing our mechanistic understanding of progression of neurodegenerative diseases.

Neuronal TIMP2 regulates hippocampus-dependent plasticity and extracellular matrix complexity.

Functional output of the hippocampus, a brain region subserving memory function, depends on highly orchestrated cellular and molecular processes that regulate synaptic plasticity throughout life. The structural requirements of such plasticity and molecular events involved in this regulation are poorly understood. Specific molecules, including tissue inhibitor of metalloproteinases-2 (TIMP2) have been implicated in plasticity processes in the hippocampus, a role that decreases with brain aging as expression is lost. Here, we report that TIMP2 is highly expressed by neurons within the hippocampus and its loss drives changes in cellular programs related to adult neurogenesis and dendritic spine turnover with corresponding impairments in hippocampus-dependent memory. Consistent with the accumulation of extracellular matrix (ECM) in the hippocampus we observe with aging, we find that TIMP2 acts to reduce accumulation of ECM around synapses in the hippocampus. Moreover, its deletion results in hindrance of newborn neuron migration through a denser ECM network. A novel conditional TIMP2 knockout (KO) model reveals that neuronal TIMP2 regulates adult neurogenesis, accumulation of ECM, and ultimately hippocampus-dependent memory. Our results define a mechanism whereby hippocampus-dependent function is regulated by TIMP2 and its interactions with the ECM to regulate diverse processes associated with synaptic plasticity.

Evolutionary scaling and cognitive correlates of primate frontal cortex microstructure.

Investigating evolutionary changes in frontal cortex microstructure is crucial to understanding how modifications of neuron and axon distributions contribute to phylogenetic variation in cognition. In the present study, we characterized microstructural components of dorsolateral prefrontal cortex, orbitofrontal cortex, and primary motor cortex from 14 primate species using measurements of neuropil fraction and immunohistochemical markers for fast-spiking inhibitory interneurons, large pyramidal projection neuron subtypes, serotonergic innervation, and dopaminergic innervation. Results revealed that the rate of evolutionary change was similar across these microstructural variables, except for neuropil fraction, which evolves more slowly and displays the strongest correlation with brain size. We also found that neuropil fraction in orbitofrontal cortex layers V-VI was associated with cross-species variation in performance on experimental tasks that measure self-control. These findings provide insight into the evolutionary reorganization of the primate frontal cortex in relation to brain size scaling and its association with cognitive processes.

Enhancing the brain MRI at ultra-high field systems using a meta-array structure.

BACKGROUND: The main advantage of ultra-high field (UHF) magnetic resonance neuroimaging is theincreased signal-to-noise ratio (SNR) compared with lower field strength imaging. However, the wavelength effect associated with UHF MRI results in radiofrequency (RF) inhomogeneity, compromising whole brain coverage for many commercial coils. Approaches to resolving this issue of transmit field inhomogeneity include the design of parallel transmit systems (PTx), RF pulse design, and applying passive RF shimming such as high dielectric materials. However, these methods have some drawbacks such as unstable material parameters of dielectric pads, high-cost, and complexity of PTx systems. Metasurfaces are artificial structures with a unique platform that can control the propagation of the electromagnetic (EM) waves, and they are very promising for engineering EM device. Implementation of meta-arrays enhancing MRI has been explored previously in several studies. PURPOSE: The aim of this study was to assess the effect of new meta-array technology on enhancing the brain MRI at 7T. A meta-array based on a hybrid structure consisting of an array of broadside-coupled split-ring resonators and high-permittivity materials was designed to work at the Larmor frequency of a 7 Tesla (7T) MRI scanner. When placed behind the head and neck, this construct improves the SNR in the region of the cerebellum,brainstem and the inferior aspect of the temporal lobes. METHODS: Numerical electromagnetic simulations were performed to optimize the meta-array design parameters and determine the RF circuit configuration. The resultant transmit-efficiency and signal sensitivity improvements were experimentally analyzed in phantoms followed by healthy volunteers using a 7T whole-body MRI scanner equipped with a standard one-channel transmit, 32-channel receive head coil. Efficacy was evaluated through acquisition with and without the meta-array using two basic sequences: gradient-recalled-echo (GRE) and turbo-spin-echo (TSE). RESULTS: Experimental phantom analysis confirmed two-fold improvement in the transmit efficiency and 1.4-fold improvement in the signal sensitivity in the target region. In vivo GRE and TSE images with the meta-array in place showed enhanced visualization in inferior regions of the brain, especially of the cerebellum, brainstem, and cervical spinal cord. CONCLUSION: Addition of the meta-array to commonly used MRI coils can enhance SNR to extend the anatomical coverage of the coil and improve overall MRI coil performance. This enhancement in SNR can be leveraged to obtain a higher resolution image over the same time slot or faster acquisition can be achieved with same resolution. Using this technique could improve the performance of existing commercial coils at 7T for whole brain and other applications.

A cellular resolution atlas of Broca's area.

Brain cells are arranged in laminar, nuclear, or columnar structures, spanning a range of scales. Here, we construct a reliable cell census in the frontal lobe of human cerebral cortex at micrometer resolution in a magnetic resonance imaging (MRI)-referenced system using innovative imaging and analysis methodologies. MRI establishes a macroscopic reference coordinate system of laminar and cytoarchitectural boundaries. Cell counting is obtained with a digital stereological approach on the 3D reconstruction at cellular resolution from a custom-made inverted confocal light-sheet fluorescence microscope (LSFM). Mesoscale optical coherence tomography enables the registration of the distorted histological cell typing obtained with LSFM to the MRI-based atlas coordinate system. The outcome is an integrated high-resolution cellular census of Broca's area in a human postmortem specimen, within a whole-brain reference space atlas.

Transcriptomic cytoarchitecture reveals principles of human neocortex organization.

Variation in cytoarchitecture is the basis for the histological definition of cortical areas. We used single cell transcriptomics and performed cellular characterization of the human cortex to better understand cortical areal specialization. Single-nucleus RNA-sequencing of 8 areas spanning cortical structural variation showed a highly consistent cellular makeup for 24 cell subclasses. However, proportions of excitatory neuron subclasses varied substantially, likely reflecting differences in connectivity across primary sensorimotor and association cortices. Laminar organization of astrocytes and oligodendrocytes also differed across areas. Primary visual cortex showed characteristic organization with major changes in the excitatory to inhibitory neuron ratio, expansion of layer 4 excitatory neurons, and specialized inhibitory neurons. These results lay the groundwork for a refined cellular and molecular characterization of human cortical cytoarchitecture and areal specialization.

Single basolateral amygdala neurons in macaques exhibit distinct connectional motifs with frontal cortex.

Basolateral amygdala (BLA) projects widely across the macaque frontal cortex, and amygdalo-frontal projections are critical for appropriate emotional responding and decision making. While it is appreciated that single BLA neurons branch and project to multiple areas in frontal cortex, the organization and frequency of this branching has yet to be fully characterized. Here, we determined the projection patterns of more than 3,000 macaque BLA neurons. We found that one-third of BLA neurons had two or more distinct projection targets in frontal cortex and subcortical structures. The patterns of single BLA neuron projections to multiple areas were organized into repeating motifs that targeted distinct sets of areas in medial and ventral frontal cortex, indicative of separable BLA networks. Our findings begin to reveal the rich structure of single-neuron connections in the non-human primate brain, providing a neuroanatomical basis for the role of BLA in coordinating brain-wide responses to valent stimuli.