The Subcortical Atlas of the Rhesus Macaque (SARM) for neuroimaging

Mapping sagittal-plane reference brain atlas of the cynomolgus macaque (Macaca fascicularis) based on consecutive cytoarchitectonic images

The brain atlas is essential for exploring the anatomical structure and function of the brain. Non-human primates, such as cynomolgus macaque, have received increasing attention due to their genetic similarity to humans. However, current macaque brain atlases only offer coarse sections with intervals along the coronal direction, failing to meet the needs of single-cell resolution studies in functional and multi-omics research of the macaque brain. To address this issue, we utilized fluorescence micro-optical sectioning tomography to obtain sub-micron resolution cytoarchitectonic images of the macaque brain at the sagittal plane. Based on the obtained 8000 image sequences, a reference brain atlas comprising 45 sagittal sections was created, delineating 270 brain regions other than the cortex. Additionally, a website was established to share the reference atlas corresponding image data. This study is expected to provide an essential dataset and tool for scientists studying the macaque brain.

Evaluation of registration-based vs. manual segmentation of rhesus macaque brain MRIs

With increasing numbers of magnetic resonance imaging (MRI) datasets becoming publicly available, researchers and clinicians alike have turned to automated methods of segmentation to enable population-level analyses of these data. Although prior research has evaluated the extent to which automated methods recapitulate "gold standard" manual segmentation methods in the human brain, such an evaluation has not yet been carried out for segmentation of MRIs of the macaque brain. Macaques offer the important opportunity to bridge gaps between microanatomical studies using invasive methods like tract tracing, neural recordings, and high-resolution histology and non-invasive macroanatomical studies using methods like MRI. As such, it is important to evaluate whether automated tools derive data of sufficient quality from macaque MRIs to bridge these gaps. We tested the relationship between automated registration-based segmentation using an open source and actively maintained NHP imaging analysis pipeline (AFNI) and gold standard manual segmentation of 4 structures (2 cortical: anterior cingulate cortex and insula; 2 subcortical: amygdala and caudate) across 37 rhesus macaques (Macaca mulatta). We identified some variability in the strength of correlation between automated and manual segmentations across neural regions and differences in relationships with demographic variables like age and sex between the two techniques.

Evolutionarily conserved fMRI network dynamics in the mouse, macaque, and human brain

Evolutionarily relevant networks have been previously described in several mammalian species using time-averaged analyses of fMRI time-series. However, fMRI network activity is highly dynamic and continually evolves over timescales of seconds. Whether the dynamic organization of resting-state fMRI network activity is conserved across mammalian species remains unclear. Using frame-wise clustering of fMRI time-series, we find that intrinsic fMRI network dynamics in awake male macaques and humans is characterized by recurrent transitions between a set of 4 dominant, neuroanatomically homologous fMRI coactivation modes (C-modes), three of which are also plausibly represented in the male rodent brain. Importantly, in all species C-modes exhibit species-invariant dynamic features, including preferred occurrence at specific phases of fMRI global signal fluctuations, and a state transition structure compatible with infraslow coupled oscillator dynamics. Moreover, dominant C-mode occurrence reconstitutes the static organization of the fMRI connectome in all species, and is predictive of ranking of corresponding fMRI connectivity gradients. These results reveal a set of species-invariant principles underlying the dynamic organization of fMRI networks in mammalian species, and offer novel opportunities to relate fMRI network findings across the phylogenetic tree.

Cross-Species Convergence of Functional Connectivity Changes in Thalamic Pain Across Human Patients and Model Macaques

Thalamic pain can be understood as a network reorganization disorder. This study aimed to investigate functional connectivity (FC) in human patients and a macaque model of thalamic pain. In humans, resting-state FC was compared between patients with thalamic pain and healthy individuals. Furthermore, resting-state FC was compared in macaques, before and after the induction of thalamic pain in the same individuals. FC between the amygdala of the unaffected hemisphere and the brainstem was significantly higher in patients with thalamic pain. More specifically, a significantly higher FC was observed between the basolateral amygdala and the ventral tegmental area, which also significantly predicted the value of a visual analog scale of pain intensity in individual patients. The macaque model of thalamic pain also exhibited a significantly higher FC between the amygdala of the unaffected hemisphere and the brainstem, particularly between the basolateral amygdala and the midbrain. Furthermore, the previously reported significantly higher FC between the amygdala and the mediodorsal nucleus of the thalamus in macaques with thalamic pain was also reproduced in the human patients. Therefore, the present results suggest that the FC changes in the regions associated with emotion, memory, motivation, and reward are part of the underlying mechanisms of thalamic pain onset present in both human patients and model macaques. This cross-species convergence provides new insights into the neurological mechanisms underlying thalamic pain, paving the way for further studies and the development of therapeutic strategies. PERSPECTIVE: This article presents that the FC changes in the regions associated with emotion, motivation, and reward are part of the underlying mechanisms of thalamic pain in humans and macaques.

Brain structure and activity predicting cognitive maturation in adolescence

Cognitive abilities of primates, including humans, continue to improve through adolescence 1,2. While a range of changes in brain structure and connectivity have been documented 3,4, how they affect neuronal activity that ultimately determines performance of cognitive functions remains unknown. Here, we conducted a multilevel longitudinal study of monkey adolescent neurocognitive development. The developmental trajectory of neural activity in the prefrontal cortex accounted remarkably well for working memory improvements. While complex aspects of activity changed progressively during adolescence, such as the rotation of stimulus representation in multidimensional neuronal space, which has been implicated in cognitive flexibility, even simpler attributes, such as the baseline firing rate in the period preceding a stimulus appearance had predictive power over behavior. Unexpectedly, decreases in brain volume and thickness, which are widely thought to underlie cognitive changes in humans 5 did not predict well the trajectory of neural activity or cognitive performance changes. Whole brain cortical volume in particular, exhibited an increase and reached a local maximum in late adolescence, at a time of rapid behavioral improvement. Maturation of long-distance white matter tracts linking the frontal lobe with areas of the association cortex and subcortical regions best predicted changes in neuronal activity and behavior. Our results provide evidence that optimization of neural activity depending on widely distributed circuitry effects cognitive development in adolescence.

State-specific Regulation of Electrical Stimulation in the Intralaminar Thalamus of Macaque Monkeys: Network and Transcriptional Insights into Arousal

Long-range thalamocortical communication is central to anesthesia-induced loss of consciousness and its reversal. However, isolating the specific neural networks connecting thalamic nuclei with various cortical regions for state-specific anesthesia regulation is challenging, with the biological underpinnings still largely unknown. Here, simultaneous electroencephalogram-fuctional magnetic resonance imaging (EEG-fMRI) and deep brain stimulation are applied to the intralaminar thalamus in macaques under finely-tuned propofol anesthesia. This approach led to the identification of an intralaminar-driven network responsible for rapid arousal during slow-wave oscillations. A network-based RNA-sequencing analysis is conducted of region-, layer-, and cell-specific gene expression data from independent transcriptomic atlases and identifies 2489 genes preferentially expressed within this arousal network, notably enriched in potassium channels and excitatory, parvalbumin-expressing neurons, and oligodendrocytes. Comparison with human RNA-sequencing data highlights conserved molecular and cellular architectures that enable the matching of homologous genes, protein interactions, and cell types across primates, providing novel insight into network-focused transcriptional signatures of arousal.

A brain subcortical segmentation tool based on anatomy attentional fusion network for developing macaques

Magnetic Resonance Imaging (MRI) plays a pivotal role in the accurate measurement of brain subcortical structures in macaques, which is crucial for unraveling the complexities of brain structure and function, thereby enhancing our understanding of neurodegenerative diseases and brain development. However, due to significant differences in brain size, structure, and imaging characteristics between humans and macaques, computational tools developed for human neuroimaging studies often encounter obstacles when applied to macaques. In this context, we propose an Anatomy Attentional Fusion Network (AAF-Net), which integrates multimodal MRI data with anatomical constraints in a multi-scale framework to address the challenges posed by the dynamic development, regional heterogeneity, and age-related size variations of the juvenile macaque brain, thus achieving precise subcortical segmentation. Specifically, we generate a Signed Distance Map (SDM) based on the initial rough segmentation of the subcortical region by a network as an anatomical constraint, providing comprehensive information on positions, structures, and morphology. Then we construct AAF-Net to fully fuse the SDM anatomical constraints and multimodal images for refined segmentation. To thoroughly evaluate the performance of our proposed tool, over 700 macaque MRIs from 19 datasets were used in this study. Specifically, we employed two manually labeled longitudinal macaque datasets to develop the tool and complete four-fold cross-validations. Furthermore, we incorporated various external datasets to demonstrate the proposed tool's generalization capabilities and promise in brain development research. We have made this tool available as an open-source resource at https://github.com/TaoZhong11/Macaque_subcortical_segmentation for direct application.

Macaque Brainnetome Atlas: A multifaceted brain map with parcellation, connection, and histology

The rhesus macaque (Macaca mulatta) is a crucial experimental animal that shares many genetic, brain organizational, and behavioral characteristics with humans. A macaque brain atlas is fundamental to biomedical and evolutionary research. However, even though connectivity is vital for understanding brain functions, a connectivity-based whole-brain atlas of the macaque has not previously been made. In this study, we created a new whole-brain map, the Macaque Brainnetome Atlas (MacBNA), based on the anatomical connectivity profiles provided by high angular and spatial resolution ex vivo diffusion MRI data. The new atlas consists of 248 cortical and 56 subcortical regions as well as their structural and functional connections. The parcellation and the diffusion-based tractography were evaluated with invasive neuronal-tracing and Nissl-stained images. As a demonstrative application, the structural connectivity divergence between macaque and human brains was mapped using the Brainnetome atlases of those two species to uncover the genetic underpinnings of the evolutionary changes in brain structure. The resulting resource includes: (1) the thoroughly delineated Macaque Brainnetome Atlas (MacBNA), (2) regional connectivity profiles, (3) the postmortem high-resolution macaque diffusion and T2-weighted MRI dataset (Brainnetome-8), and (4) multi-contrast MRI, neuronal-tracing, and histological images collected from a single macaque. MacBNA can serve as a common reference frame for mapping multifaceted features across modalities and spatial scales and for integrative investigation and characterization of brain organization and function. Therefore, it will enrich the collaborative resource platform for nonhuman primates and facilitate translational and comparative neuroscience research.

Deep brain stimulation induces white matter remodeling and functional changes to brain-wide networks

Deep brain stimulation (DBS) is an emerging therapeutic option for treatment resistant neurological and psychiatric disorders, most notably depression. Despite this, little is known about the anatomical and functional mechanisms that underlie this therapy. Here we targeted stimulation to the white matter adjacent to the subcallosal anterior cingulate cortex (SCC-DBS) in macaques, modeling the location in the brain proven effective for depression. We demonstrate that SCC-DBS has a selective effect on white matter macro- and micro-structure in the cingulum bundle distant to where stimulation was delivered. SCC-DBS also decreased functional connectivity between subcallosal and posterior cingulate cortex, two areas linked by the cingulum bundle and implicated in depression. Our data reveal that white matter remodeling as well as functional effects contribute to DBS's therapeutic efficacy.

Transient brain activity dynamics discriminate levels of consciousness during anesthesia

The awake mammalian brain is functionally organized in terms of large-scale distributed networks that are constantly interacting. Loss of consciousness might disrupt this temporal organization leaving patients unresponsive. We hypothesize that characterizing brain activity in terms of transient events may provide a signature of consciousness. For this, we analyze temporal dynamics of spatiotemporally overlapping functional networks obtained from fMRI transient activity across different anesthetics and levels of anesthesia. We first show a striking homology in spatial organization of networks between monkeys and humans, indicating cross-species similarities in resting-state fMRI structure. We then track how network organization shifts under different anesthesia conditions in macaque monkeys. While the spatial aspect of the networks is preserved, their temporal dynamics are highly affected by anesthesia. Networks express for longer durations and co-activate in an anesthetic-specific configuration. Additionally, hierarchical brain organization is disrupted with a consciousness-level-signature role of the default mode network. In conclusion, large-scale brain network temporal dynamics capture differences in anesthetic-specific consciousness-level, paving the way towards a clinical translation of these cortical signature.

Beyond faces: the contribution of the amygdala to visual processing in the macaque brain

The amygdala is present in a diverse range of vertebrate species, such as lizards, rodents, and primates; however, its structure and connectivity differs across species. The increased connections to visual sensory areas in primate species suggests that understanding the visual selectivity of the amygdala in detail is critical to revealing the principles underlying its function in primate cognition. Therefore, we designed a high-resolution, contrast-agent enhanced, event-related fMRI experiment, and scanned 3 adult rhesus macaques, while they viewed 96 naturalistic stimuli. Half of these stimuli were social (defined by the presence of a conspecific), the other half were nonsocial. We also nested manipulations of emotional valence (positive, neutral, and negative) and visual category (faces, nonfaces, animate, and inanimate) within the stimulus set. The results reveal widespread effects of emotional valence, with the amygdala responding more on average to inanimate objects and animals than faces, bodies, or social agents in this experimental context. These findings suggest that the amygdala makes a contribution to primate vision that goes beyond an auxiliary role in face or social perception. Furthermore, the results highlight the importance of stimulus selection and experimental design when probing the function of the amygdala and other visually responsive brain regions.

Macaque claustrum, pulvinar and putative dorsolateral amygdala support the cross-modal association of social audio-visual stimuli based on meaning

Social communication draws on several cognitive functions such as perception, emotion recognition and attention. The association of audio-visual information is essential to the processing of species-specific communication signals. In this study, we use functional magnetic resonance imaging in order to identify the subcortical areas involved in the cross-modal association of visual and auditory information based on their common social meaning. We identified three subcortical regions involved in audio-visual processing of species-specific communicative signals: the dorsolateral amygdala, the claustrum and the pulvinar. These regions responded to visual, auditory congruent and audio-visual stimulations. However, none of them was significantly activated when the auditory stimuli were semantically incongruent with the visual context, thus showing an influence of visual context on auditory processing. For example, positive vocalization (coos) activated the three subcortical regions when presented in the context of positive facial expression (lipsmacks) but not when presented in the context of negative facial expression (aggressive faces). In addition, the medial pulvinar and the amygdala presented multisensory integration such that audiovisual stimuli resulted in activations that were significantly higher than those observed for the highest unimodal response. Last, the pulvinar responded in a task-dependent manner, along a specific spatial sensory gradient. We propose that the dorsolateral amygdala, the claustrum and the pulvinar belong to a multisensory network that modulates the perception of visual socioemotional information and vocalizations as a function of the relevance of the stimuli in the social context. SIGNIFICANCE STATEMENT: Understanding and correctly associating socioemotional information across sensory modalities, such that happy faces predict laughter and escape scenes predict screams, is essential when living in complex social groups. With the use of functional magnetic imaging in the awake macaque, we identify three subcortical structures-dorsolateral amygdala, claustrum and pulvinar-that only respond to auditory information that matches the ongoing visual socioemotional context, such as hearing positively valenced coo calls and seeing positively valenced mutual grooming monkeys. We additionally describe task-dependent activations in the pulvinar, organizing along a specific spatial sensory gradient, supporting its role as a network regulator.

The neural basis of resting-state fMRI functional connectivity in fronto-limbic circuits revealed by chemogenetic manipulation

Measures of fMRI resting-state functional connectivity (rs-FC) are an essential tool for basic and clinical investigations of fronto-limbic circuits. Understanding the relationship between rs-FC and the underlying patterns of neural activity in these circuits is therefore vital. Here we introduced inhibitory designer receptors exclusively activated by designer drugs (DREADDs) into the amygdala of two male macaques. We evaluated the causal effect of activating the DREADD receptors on rs-FC and neural activity within circuits connecting amygdala and frontal cortex. Activating the inhibitory DREADD increased rs-FC between amygdala and ventrolateral prefrontal cortex. Neurophysiological recordings revealed that the DREADD-induced increase in fMRI rs-FC was associated with increased local field potential coherency in the alpha band (6.5-14.5 Hz) between amygdala and ventrolateral prefrontal cortex. Thus, our multi-modal approach reveals the specific signature of neuronal activity that underlies rs-FC in fronto-limbic circuits.

Evolutionarily conserved neural responses to affective touch in monkeys transcend consciousness and change with age

Affective touch-a slow, gentle, and pleasant form of touch-activates a different neural network than which is activated during discriminative touch in humans. Affective touch perception is enabled by specialized low-threshold mechanoreceptors in the skin with unmyelinated fibers called C tactile (CT) afferents. These CT afferents are conserved across mammalian species, including macaque monkeys. However, it is unknown whether the neural representation of affective touch is the same across species and whether affective touch's capacity to activate the hubs of the brain that compute socioaffective information requires conscious perception. Here, we used functional MRI to assess the preferential activation of neural hubs by slow (affective) vs. fast (discriminative) touch in anesthetized rhesus monkeys (Macaca mulatta). The insula, anterior cingulate cortex (ACC), amygdala, and secondary somatosensory cortex were all significantly more active during slow touch relative to fast touch, suggesting homologous activation of the interoceptive-allostatic network across primate species during affective touch. Further, we found that neural responses to affective vs. discriminative touch in the insula and ACC (the primary cortical hubs for interoceptive processing) changed significantly with age. Insula and ACC in younger animals differentiated between slow and fast touch, while activity was comparable between conditions for aged monkeys (equivalent to >70 y in humans). These results, together with prior studies establishing conserved peripheral nervous system mechanisms of affective touch transduction, suggest that neural responses to affective touch are evolutionarily conserved in monkeys, significantly impacted in old age, and do not necessitate conscious experience of touch.

Involvement of neurons in the non-human primate anterior striatum in proactive inhibition

Behaving as desired requires selecting the appropriate behavior and inhibiting the selection of inappropriate behavior. This inhibitory function involves multiple processes, such as reactive and proactive inhibition, instead of a single process. In this study, macaque monkeys were required to perform a task in which they had to sequentially select (accept) or refuse (reject) a choice. Neural activity was recorded from the anterior striatum, which is considered to be involved in behavioral inhibition, focusing on the distinction between proactive and reactive inhibitions. We identified neurons with significant activity changes during the rejection of bad objects. Cluster analysis revealed three distinct groups, of which one showed obviously increased activity during object rejection, suggesting its involvement in proactive inhibition. This activity pattern was consistent irrespective of the rejection method, indicating a role beyond mere saccadic suppression. Furthermore, minimal activity changes during the fixation task indicated that these neurons were not primarily involved in reactive inhibition. In conclusion, these findings suggest that the anterior striatum plays a crucial role in cognitive control and orchestrates goal-directed behavior through proactive inhibition, which may be critical in understanding the mechanisms of behavioral inhibition dysfunction that occur in patients with basal ganglia disease.

Differential functional organization of amygdala-medial prefrontal cortex networks in macaque and human

Over the course of evolution, the amygdala (AMG) and medial frontal cortex (mPFC) network, involved in behavioral adaptation, underwent structural changes in the old-world monkey and human lineages. Yet, whether and how the functional organization of this network differs remains poorly understood. Using resting-state functional magnetic resonance imagery, we show that the functional connectivity (FC) between AMG nuclei and mPFC regions differs between humans and awake macaques. In humans, the AMG-mPFC FC displays U-shaped pattern along the corpus callosum: a positive FC with the ventromedial prefrontal (vmPFC) and anterior cingulate cortex (ACC), a negative FC with the anterior mid-cingulate cortex (MCC), and a positive FC with the posterior MCC. Conversely, in macaques, the negative FC shifted more ventrally at the junction between the vmPFC and the ACC. The functional organization divergence of AMG-mPFC network between humans and macaques might help understanding behavioral adaptation abilities differences in their respective socio-ecological niches.

Magnetic Resonance Imaging-Guided Frameless Stereotactic Injections of the Bilateral Cerebellar Dentate Nuclei in Nonhuman Primates: Technical Note

BACKGROUND AND OBJECTIVES

Nonhuman primates (NHPs) are important preclinical models for evaluating therapeutics because of their anatomophysiological similarities to humans, and can be especially useful for testing new delivery targets. With the growing promise of cell and gene therapies for the treatment of neurological diseases, it is important to ensure the accurate and safe delivery of these agents to target structures in the brain. However, a standard guideline or method has not been developed for stereotactic targeting in NHPs. In this article, we describe the safe use of a magnetic resonance imaging-guided frameless stereotactic system to target bilateral cerebellar dentate nuclei for accurate, real-time delivery of viral vector in NHPs.

METHODS

Seventeen rhesus macaques (Macaca mulatta) underwent stereotactic surgery under real-time MRI guidance using the ClearPoint® system. Bilateral cerebellar dentate nuclei were targeted through a single parietal entry point with a transtentorial approach. Fifty microliters of contrast-impregnated infusate was delivered to each dentate nucleus, and adjustments were made as necessary according to real-time MRI monitoring of delivery. Perioperative clinical outcomes and postoperative volumes of distribution were recorded.

RESULTS

All macaques underwent bilateral surgery successfully. Superficial pin site infection occurred in 4/17 (23.5%) subjects, which resolved with antibiotics. Two episodes of transient neurological deficit (anisocoria and unilateral weakness) were recorded, which did not require additional postoperative treatment and resolved over time. Volume of distribution of infusate achieved satisfactory coverage of target dentate nuclei, and only 1 incidence (2.9%) of cerebrospinal fluid penetration was recorded. Mean volume of distribution was 161.22 ± 39.61 mm3 (left, 173.65 ± 48.29; right, 148.80 ± 23.98).

CONCLUSION

MRI-guided frameless stereotactic injection of bilateral cerebellar dentate nuclei in NHPs is safe and feasible. The use of this technique enables real-time modification of the surgical plan to achieve adequate target coverage and can be readily translated to clinical use.

Pharmacological modulation of dopamine D1 and D2 receptors reveals distinct neural networks related to probabilistic learning in non-human primates

The neurotransmitter dopamine (DA) has a multifaceted role in healthy and disordered brains through its action on multiple subtypes of dopaminergic receptors. How modulation of these receptors controls behavior by altering connectivity across intrinsic brain-wide networks remains elusive. Here we performed parallel behavioral and resting-state functional MRI experiments after administration of two different DA receptor antagonists in macaque monkeys. Systemic administration of SCH-23390 (D1 antagonist) disrupted probabilistic learning when subjects had to learn new stimulus-reward associations and diminished functional connectivity (FC) in cortico-cortical and fronto-striatal connections. By contrast, haloperidol (D2 antagonist) improved learning and broadly enhanced FC in cortical connections. Further comparison between the effect of SCH-23390/haloperidol on behavioral and resting-state FC revealed specific cortical and subcortical networks associated with the cognitive and motivational effects of DA, respectively. Thus, we reveal the distinct brain-wide networks that are associated with the dopaminergic control of learning and motivation via DA receptors.

Neural Integration of Audiovisual Sensory Inputs in Macaque Amygdala and Adjacent Regions

Integrating multisensory inputs to generate accurate perception and guide behavior is among the most critical functions of the brain. Subcortical regions such as the amygdala are involved in sensory processing including vision and audition, yet their roles in multisensory integration remain unclear. In this study, we systematically investigated the function of neurons in the amygdala and adjacent regions in integrating audiovisual sensory inputs using a semi-chronic multi-electrode array and multiple combinations of audiovisual stimuli. From a sample of 332 neurons, we showed the diverse response patterns to audiovisual stimuli and the neural characteristics of bimodal over unimodal modulation, which could be classified into four types with differentiated regional origins. Using the hierarchical clustering method, neurons were further clustered into five groups and associated with different integrating functions and sub-regions. Finally, regions distinguishing congruent and incongruent bimodal sensory inputs were identified. Overall, visual processing dominates audiovisual integration in the amygdala and adjacent regions. Our findings shed new light on the neural mechanisms of multisensory integration in the primate brain.

HumanBrainAtlas: an in vivo MRI dataset for detailed segmentations

We introduce HumanBrainAtlas, an initiative to construct a highly detailed, open-access atlas of the living human brain that combines high-resolution in vivo MR imaging and detailed segmentations previously possible only in histological preparations. Here, we present and evaluate the first step of this initiative: a comprehensive dataset of two healthy male volunteers reconstructed to a 0.25 mm isotropic resolution for T1w, T2w, and DWI contrasts. Multiple high-resolution acquisitions were collected for each contrast and each participant, followed by averaging using symmetric group-wise normalisation (Advanced Normalisation Tools). The resulting image quality permits structural parcellations rivalling histology-based atlases, while maintaining the advantages of in vivo MRI. For example, components of the thalamus, hypothalamus, and hippocampus are often impossible to identify using standard MRI protocols-can be identified within the present data. Our data are virtually distortion free, fully 3D, and compatible with the existing in vivo Neuroimaging analysis tools. The dataset is suitable for teaching and is publicly available via our website (hba.neura.edu.au), which also provides data processing scripts. Instead of focusing on coordinates in an averaged brain space, our approach focuses on providing an example segmentation at great detail in the high-quality individual brain. This serves as an illustration on what features contrasts and relations can be used to interpret MRI datasets, in research, clinical, and education settings.

Rhesus monkeys exhibiting spontaneous ritualistic behaviors resembling obsessive-compulsive disorder

Obsessive-compulsive disorder (OCD) is a chronic and debilitating psychiatric disorder that affects ∼2%-3% of the population globally. Studying spontaneous OCD-like behaviors in non-human primates may improve our understanding of the disorder. In large rhesus monkey colonies, we found 10 monkeys spontaneously exhibiting persistent sequential motor behaviors (SMBs) in individual-specific sequences that were repetitive, time-consuming and stable over prolonged periods. Genetic analysis revealed severely damaging mutations in genes associated with OCD risk in humans. Brain imaging showed that monkeys with SMBs had larger gray matter (GM) volumes in the left caudate nucleus and lower fractional anisotropy of the corpus callosum. The GM volume of the left caudate nucleus correlated positively with the daily duration of SMBs. Notably, exposure to a stressor (human presence) significantly increased SMBs. In addition, fluoxetine, a serotonergic medication commonly used for OCD, decreased SMBs in these monkeys. These findings provide a novel foundation for developing better understanding and treatment of OCD.

The neural basis of resting-state fMRI functional connectivity in fronto-limbic circuits revealed by chemogenetic manipulation

Measures of fMRI resting-state functional connectivity (rs-FC) are an essential tool for basic and clinical investigations of fronto-limbic circuits. Understanding the relationship between rs-FC and neural activity in these circuits is therefore vital. Here we introduced inhibitory designer receptors exclusively activated by designer drugs (DREADDs) into the macaque amygdala and activated them with a highly selective and potent DREADD agonist, deschloroclozapine. We evaluated the causal effect of activating the DREADD receptors on rs-FC and neural activity within circuits connecting amygdala and frontal cortex. Interestingly, activating the inhibitory DREADD increased rs-FC between amygdala and ventrolateral prefrontal cortex. Neurophysiological recordings revealed that the DREADD-induced increase in fMRI rs-FC was associated with increased local field potential coherency in the alpha band (6.5-14.5Hz) between amygdala and ventrolateral prefrontal cortex. Thus, our multi-disciplinary approach reveals the specific signature of neuronal activity that underlies rs-FC in fronto-limbic circuits.

Longitudinal positron emission tomography and postmortem analysis reveals widespread neuroinflammation in SARS-CoV-2 infected rhesus macaques

BACKGROUND

Coronavirus disease 2019 (COVID-19) patients initially develop respiratory symptoms, but they may also suffer from neurological symptoms. People with long-lasting effects after acute infections with severe respiratory syndrome coronavirus 2 (SARS-CoV-2), i.e., post-COVID syndrome or long COVID, may experience a variety of neurological manifestations. Although we do not fully understand how SARS-CoV-2 affects the brain, neuroinflammation likely plays a role.

METHODS

To investigate neuroinflammatory processes longitudinally after SARS-CoV-2 infection, four experimentally SARS-CoV-2 infected rhesus macaques were monitored for 7 weeks with 18-kDa translocator protein (TSPO) positron emission tomography (PET) using [¹⁸F]DPA714, together with computed tomography (CT). The baseline scan was compared to weekly PET-CTs obtained post-infection (pi). Brain tissue was collected following euthanasia (50 days pi) to correlate the PET signal with TSPO expression, and glial and endothelial cell markers. Expression of these markers was compared to brain tissue from uninfected animals of comparable age, allowing the examination of the contribution of these cells to the neuroinflammatory response following SARS-CoV-2 infection.

RESULTS

TSPO PET revealed an increased tracer uptake throughout the brain of all infected animals already from the first scan obtained post-infection (day 2), which increased to approximately twofold until day 30 pi. Postmortem immunohistochemical analysis of the hippocampus and pons showed TSPO expression in cells expressing ionized calcium-binding adaptor molecule 1 (IBA1), glial fibrillary acidic protein (GFAP), and collagen IV. In the hippocampus of SARS-CoV-2 infected animals the TSPO⁺ area and number of TSPO⁺ cells were significantly increased compared to control animals. This increase was not cell type specific, since both the number of IBA1⁺TSPO⁺ and GFAP⁺TSPO⁺ cells was increased, as well as the TSPO⁺ area within collagen IV⁺ blood vessels.

CONCLUSIONS

This study manifests [¹⁸F]DPA714 as a powerful radiotracer to visualize SARS-CoV-2 induced neuroinflammation. The increased uptake of [¹⁸F]DPA714 over time implies an active neuroinflammatory response following SARS-CoV-2 infection. This inflammatory signal coincides with an increased number of TSPO expressing cells, including glial and endothelial cells, suggesting neuroinflammation and vascular dysregulation. These results demonstrate the long-term neuroinflammatory response following a mild SARS-CoV-2 infection, which potentially precedes long-lasting neurological symptoms.

Multiple routes of communication within the amygdala-mPFC network: A comparative approach in humans and macaques

The network formed by the amygdala (AMG) and the medial Prefrontal Cortex (mPFC), at the interface between our internal and external environment, has been shown to support some important aspects of behavioral adaptation. Whether and how the anatomo-functional organization of this network evolved across primates remains unclear. Here, we compared AMG nuclei morphological characteristics and their functional connectivity with the mPFC in humans and macaques to identify potential homologies and differences between these species. Based on selected studies, we highlight two subsystems within the AMG-mPFC circuits, likely involved in distinct temporal dynamics of integration during behavioral adaptation. We also show that whereas the mPFC displays a large expansion but a preserved intrinsic anatomo-functional organization, the AMG displays a volume reduction and morphological changes related to specific nuclei. We discuss potential commonalities and differences in the dialogue between AMG nuclei and mPFC in humans and macaques based on available data.

High-throughput sequencing of macaque basolateral amygdala projections reveals dissociable connectional motifs with frontal cortex

The basolateral amygdala (BLA) projects widely across the macaque frontal cortex1-4, and amygdalo-frontal projections are critical for optimal emotional responding5 and decision-making6. Yet, little is known about the single-neuron architecture of these projections: namely, whether single BLA neurons project to multiple parts of the frontal cortex. Here, we use MAPseq7 to determine the projection patterns of over 3000 macaque BLA neurons. We found that one-third of BLA neurons have two or more distinct targets in parts of frontal cortex and of subcortical structures. Further, we reveal non-random structure within these branching patterns such that neurons with four targets are more frequently observed than those with two or three, indicative of widespread networks. Consequently, these multi-target single neurons form distinct networks within medial and ventral frontal cortex consistent with their known functions in regulating mood and decision-making. Additionally, we show that branching patterns of single neurons shape functional networks in the brain as assessed by fMRI-based functional connectivity. These results provide a neuroanatomical basis for the role of the BLA in coordinating brain-wide responses to valent stimuli8 and highlight the importance of high-resolution neuroanatomical data for understanding functional networks in the brain.

A multimodal imaging-guided software for access to primate brains

Background

Imaging-guided access to the brain has become a routine procedure for various research and clinical applications, including drug administration, neurophysiological recording, and sampling tissue. Therefore, open-source software is required to handle such datasets in these specific applications.

New methods

Here, we proposed an open-source tool utilizing different imaging modalities for automating the steps to access the brain. This tool provides means for easily calculating the coordination of the area of interest concerning a specific point of entry. The source and documentation are available at this link.

Results

We have used this software for three different applications: electrophysiological recording, drug infusion in the nonhuman primate brain, and guided biopsy procedure in the human brain. We performed a neural recording of two monkeys' prefrontal cortex and inferior temporal cortex using this software in submillimeter resolution. We also applied our procedure for infusion in the putamen and caudate nuclei in both hemispheres of another group of rhesus monkeys with histological proof in one animal. More so, we validated this software in the human subjects that underwent biopsy surgery with the commercial software used in human biopsy surgery.

Comparison with existing methods

Our software uses different imaging modalities by co-registering them. This will provide structural details of the skull and brain tissue. We can calculate each brain region's coordination at the point of entry by re-slicing the images. Atlas-based image segmentation were implemented in our software. Three mentioned applications of our software in neuroscience will be further discussed in this paper.

Conclusion

In our procedure, working with different imaging modalities provides a precise estimation of the specific region in the brain related to the location of implants or stereotaxic frames. There is no limitation to using metal implants in this procedure.

The use of chemogenetic actuator ligands in nonhuman primate DREADDs-fMRI

Designer Receptors Exclusively Activated by Designer Drugs (DREADDs) are engineered receptors that allow for genetically targeted, reversible manipulation of cellular activity via systemic drug administration. DREADD induced manipulations are initiated via the binding of an actuator ligand. Therefore, the use of DREADDs is contingent on the availability of actuator ligands. Actuator ligands low-dose clozapine (CLZ) and deschloroclozapine (DCZ) are highly selective for DREADDs, and, upon binding, induce physiological and behavioral changes in rodents and nonhuman primates (NHPs). Despite this reported specificity, both CLZ and DCZ have partial affinity for a variety of endogenous receptors and can induce dose-specific changes even in naïve animals. As such, this study aimed to examine the effects of CLZ and DCZ on resting-state functional connectivity (rs-FC) and intrinsic neural timescales (INTs) in naïve NHPs. In doing so, we evaluated whether CLZ and DCZ - in the absence of DREADDs - are inert by examining these ligands' effects on the intrinsic functional properties of the brain. Low-dose DCZ did not induce consistent changes in rs-FC or INTs prior to the expression of DREADDs; however, a high dose resulted in subject-specific changes in rs-FC and INTs. In contrast, CLZ administration induced consistent changes in rs-FC and INTs prior to DREADD expression in our subjects. Our results caution against the use of CLZ by explicitly demonstrating the impact of off-target effects that can confound experimental results. Altogether, these data endorse the use of low dose DCZ for future DREADD-based experiments.

Fusion of quantitative susceptibility maps and T1-weighted images improve brain tissue contrast in primates

Recent progress in quantitative susceptibility mapping (QSM) has enabled the accurate delineation of submillimeter-scale subcortical brain structures in humans. However, the simultaneous visualization of cortical, subcortical, and white matter structure remains challenging, utilizing QSM data solely. Here we present TQ-SILiCON, a fusion method that enhances the contrast of cortex and subcortical structures and provides an excellent white matter delineation by combining QSM and conventional T1-weighted (T1w) images. In this study, we first applied QSM in the macaque monkey to map iron-rich subcortical structures. Implementing the same QSM acquisition and analysis methods allowed a similar accurate delineation of subcortical structures in humans. However, the QSM contrast of white and cortical gray matter was not sufficient for appropriate segmentation. Applying automatic brain tissue segmentation to TQ-SILiCON images of the macaque improved the classification of subcortical brain structures as compared to the single T1 contrast by maintaining an excellent white to cortical gray matter contrast. Furthermore, we validated our dual-contrast fusion approach in humans and similarly demonstrated improvements in automated segmentation of the cortex and subcortical structures. We believe the proposed contrast will facilitate translational studies in nonhuman primates to investigate the pathophysiology of neurodegenerative diseases that affect subcortical structures such as the basal ganglia in humans.

A broadly tuned network for affective body language in the macaque brain

Body language is a powerful tool that we use to communicate how we feel, but it is unclear whether other primates also communicate in this way. Here, we use functional magnetic resonance imaging to show that the body-selective patches in macaques are activated by affective body language. Unexpectedly, we found these regions to be tolerant of naturalistic variation in posture as well as species; the bodies of macaques, humans, and domestic cats all evoked a stronger response when they conveyed fear than when they conveyed no affect. Multivariate analyses confirmed that the neural representation of fear-related body expressions was species-invariant. Collectively, these findings demonstrate that, like humans, macaques have body-selective brain regions in the ventral visual pathway for processing affective body language. These data also indicate that representations of body stimuli in these regions are built on the basis of emergent properties, such as socio-affective meaning, and not just putative image properties.

Abnormal regional homogeneity and amplitude of low frequency fluctuation in chronic kidney patients with and without dialysis

Purpose

The study characterizes regional homogeneity (ReHo) and amplitude of low frequency fluctuations (ALFF) in abnormal regions of brain in patients of chronic kidney disease (CKD).

Materials and methods

A total of 64 patients of CKD were divided into 26 cases of non-dialysis-dependent chronic kidney disease (NDD-CKD), and 38 cases of dialysis-dependent chronic kidney disease (DD-CKD). A total of 43 healthy controls (normal control, NC) were also included. All subjects underwent resting-state functional magnetic resonance imaging (rs-fMRI). ALFF and ReHo data was processed for monitoring the differences in spontaneous brain activity between the three groups. ALFF and ReHo values of extracted differential brain regions were correlated to the clinical data and cognitive scores of CKD patients.

Results

Non-dialysis-dependent group has increased ALFF levels in 13 brain regions while that of DD group in 28 brain regions as compared with NC group. ReHo values are altered in six brain regions of DD group. ALFF is correlated with urea nitrogen and ReHo with urea nitrogen and creatinine. DD group has altered ReHo in two brain regions compared with NDD group. The differences are located in basal ganglia, cerebellar, and hippocampus regions.

Conclusion

Abnormal activity in basal ganglia, cerebellar, and hippocampal regions may be involved in the cognitive decline of CKD patients. This link can provide theoretical basis for understanding the cognitive decline.

Cerebro-cerebellar interactions in nonhuman primates examined by optogenetic functional magnetic resonance imaging

Functional magnetic resonance imaging (fMRI) is a promising approach for the simultaneous and extensive scanning of whole-brain activities. Optogenetics is free from electrical and magnetic artifacts and is an ideal stimulation method for combined use with fMRI. However, the application of optogenetics in nonhuman primates (NHPs) remains limited. Recently, we developed an efficient optogenetic intracortical microstimulation method of the primary motor cortex (M1), which successfully induced forelimb movements in macaque monkeys. Here, we aimed to investigate how optogenetic M1 stimulation causes neural modulation in the local and remote brain regions in anesthetized monkeys using 7-tesla fMRI. We demonstrated that optogenetic stimulation of the M1 forelimb and hindlimb regions successfully evoked robust direct and remote fMRI activities. Prominent remote activities were detected in the anterior and posterior lobes in the contralateral cerebellum, which receive projections polysynaptically from the M1. We further demonstrated that the cerebro-cerebellar projections from these M1 regions were topographically organized, which is concordant with the somatotopic map in the cerebellar cortex previously reported in macaques and humans. The present study significantly enhances optogenetic fMRI in NHPs, resulting in profound understanding of the brain network, thereby accelerating the translation of findings from animal models to humans.

Spatial signatures of anesthesia-induced burst-suppression differ between primates and rodents

During deep anesthesia, the electroencephalographic (EEG) signal of the brain alternates between bursts of activity and periods of relative silence (suppressions). The origin of burst-suppression and its distribution across the brain remain matters of debate. In this work, we used functional magnetic resonance imaging (fMRI) to map the brain areas involved in anesthesia-induced burst-suppression across four mammalian species: humans, long-tailed macaques, common marmosets, and rats. At first, we determined the fMRI signatures of burst-suppression in human EEG-fMRI data. Applying this method to animal fMRI datasets, we found distinct burst-suppression signatures in all species. The burst-suppression maps revealed a marked inter-species difference: in rats, the entire neocortex engaged in burst-suppression, while in primates most sensory areas were excluded-predominantly the primary visual cortex. We anticipate that the identified species-specific fMRI signatures and whole-brain maps will guide future targeted studies investigating the cellular and molecular mechanisms of burst-suppression in unconscious states.

Intrinsic timescales as an organizational principle of neural processing across the whole rhesus macaque brain

Hierarchical temporal dynamics are a fundamental computational property of the brain; however, there are no whole brain, noninvasive investigations into timescales of neural processing in animal models. To that end, we used the spatial resolution and sensitivity of ultrahigh field functional magnetic resonance imaging (fMRI) performed at 10.5 T to probe timescales across the whole macaque brain. We uncovered within-species consistency between timescales estimated from fMRI and electrophysiology. Crucially, we extended existing electrophysiological hierarchies to whole-brain topographies. Our results validate the complementary use of hemodynamic and electrophysiological intrinsic timescales, establishing a basis for future translational work. Further, with these results in hand, we were able to show that one facet of the high-dimensional functional connectivity (FC) topography of any region in the brain is closely related to hierarchical temporal dynamics. We demonstrated that intrinsic timescales are organized along spatial gradients that closely match FC gradient topographies across the whole brain. We conclude that intrinsic timescales are a unifying organizational principle of neural processing across the whole brain.

Anatomical variability, multi-modal coordinate systems, and precision targeting in the marmoset brain

Localising accurate brain regions needs careful evaluation in each experimental species due to their individual variability. However, the function and connectivity of brain areas is commonly studied using a single-subject cranial landmark-based stereotactic atlas in animal neuroscience. Here, we address this issue in a small primate, the common marmoset, which is increasingly widely used in systems neuroscience. We developed a non-invasive multi-modal neuroimaging-based targeting pipeline, which accounts for intersubject anatomical variability in cranial and cortical landmarks in marmosets. This methodology allowed creation of multi-modal templates (MarmosetRIKEN20) including head CT and brain MR images, embedded in coordinate systems of anterior and posterior commissures (AC-PC) and CIFTI grayordinates. We found that the horizontal plane of the stereotactic coordinate was significantly rotated in pitch relative to the AC-PC coordinate system (10 degrees, frontal downwards), and had a significant bias and uncertainty due to positioning procedures. We also found that many common cranial and brain landmarks (e.g., bregma, intraparietal sulcus) vary in location across subjects and are substantial relative to average marmoset cortical area dimensions. Combining the neuroimaging-based targeting pipeline with robot-guided surgery enabled proof-of-concept targeting of deep brain structures with an accuracy of 0.2 mm. Altogether, our findings demonstrate substantial intersubject variability in marmoset brain and cranial landmarks, implying that subject-specific neuroimaging-based localization is needed for precision targeting in marmosets. The population-based templates and atlases in grayordinates, created for the first time in marmoset monkeys, should help bridging between macroscale and microscale analyses.

Interconnected sub-networks of the macaque monkey gustatory connectome

Macroscopic taste processing connectivity was investigated using functional magnetic resonance imaging during the presentation of sour, salty, and sweet tastants in anesthetized macaque monkeys. This examination of taste processing affords the opportunity to study the interactions between sensory regions, central integrators, and effector areas. Here, 58 brain regions associated with gustatory processing in primates were aggregated, collectively forming the gustatory connectome. Regional regression coefficients (or β-series) obtained during taste stimulation were correlated to infer functional connectivity. This connectivity was then evaluated by assessing its laterality, modularity and centrality. Our results indicate significant correlations between same region pairs across hemispheres in a bilaterally interconnected scheme for taste processing throughout the gustatory connectome. Using unbiased community detection, three bilateral sub-networks were detected within the graph of the connectome. This analysis revealed clustering of 16 medial cortical structures, 24 lateral structures, and 18 subcortical structures. Across the three sub-networks, a similar pattern was observed in the differential processing of taste qualities. In all cases, the amplitude of the response was greatest for sweet, but the network connectivity was strongest for sour and salty tastants. The importance of each region in taste processing was computed using node centrality measures within the connectome graph, showing centrality to be correlated across hemispheres and, to a smaller extent, region volume. Connectome hubs exhibited varying degrees of centrality with a prominent leftward increase in insular cortex centrality. Taken together, these criteria illustrate quantifiable characteristics of the macaque monkey gustatory connectome and its organization as a tri-modular network, which may reflect the general medial-lateral-subcortical organization of salience and interoception processing networks.

High-resolution mapping and digital atlas of subcortical regions in the macaque monkey based on matched MAP-MRI and histology

Subcortical nuclei and other deep brain structures are known to play an important role in the regulation of the central and peripheral nervous systems. It can be difficult to identify and delineate many of these nuclei and their finer subdivisions in conventional MRI due to their small size, buried location, and often subtle contrast compared to neighboring tissue. To address this problem, we applied a multi-modal approach in ex vivo non-human primate (NHP) brain that includes high-resolution mean apparent propagator (MAP)-MRI and five different histological stains imaged with high-resolution microscopy in the brain of the same subject. By registering these high-dimensional MRI data to high-resolution histology data, we can map the location, boundaries, subdivisions, and micro-architectural features of subcortical gray matter regions in the macaque monkey brain. At high spatial resolution, diffusion MRI in general, and MAP-MRI in particular, can distinguish a large number of deep brain structures, including the larger and smaller white matter fiber tracts as well as architectonic features within various nuclei. Correlation with histology from the same brain enables a thorough validation of the structures identified with MAP-MRI. Moreover, anatomical details that are evident in images of MAP-MRI parameters are not visible in conventional T1-weighted images. We also derived subcortical template "SC21" from segmented MRI slices in three-dimensions and registered this volume to a previously published anatomical template with cortical parcellation (Reveley et al., 2017; Saleem and Logothetis, 2012), thereby integrating the 3D segmentation of both cortical and subcortical regions into the same volume. This newly updated three-dimensional D99 digital brain atlas (V2.0) is intended for use as a reference standard for macaque neuroanatomical, functional, and connectional imaging studies, involving both cortical and subcortical targets. The SC21 and D99 digital templates are available as volumes and surfaces in standard NIFTI and GIFTI formats.

Visualization of iron-rich subcortical structures in non-human primates in vivo by quantitative susceptibility mapping at 3T MRI

Magnetic resonance imaging (MRI) is now an essential tool in the field of neuroscience involving non-human primates (NHP). Structural MRI scanning using T1-weighted (T1w) or T2-weighted (T2w) images provides anatomical information, particularly for experiments involving deep structures such as the basal ganglia and cerebellum. However, for certain subcortical structures, T1w and T2w image contrasts are insufficient for their detection of important anatomical details. To better visualize such structures in the macaque brain, we applied a relatively new method called quantitative susceptibility mapping (QSM), which enhances tissue contrast based on the local tissue magnetic susceptibility. The QSM significantly improved the visualization of important structures, including the ventral pallidum (VP), globus pallidus external and internal segments (GPe and GPi), substantia nigra (SN), subthalamic nucleus (STN) in the basal ganglia and the dentate nucleus (DN) in the cerebellum. We quantified this the contrast enhancement by systematically comparing of contrast-to-noise ratios (CNRs) of QSM images relative to the corresponding T1w and T2w images. In addition, QSM values of some structures were correlated to the age of the macaque subjects. These results identify the QSM method as a straightforward and useful tool for clearly visualizing details of subcortical structures that are invisible with more traditional scanning sequences.

A comprehensive macaque fMRI pipeline and hierarchical atlas

Functional neuroimaging research in the non-human primate (NHP) has been advancing at a remarkable rate. The increase in available data establishes a need for robust analysis pipelines designed for NHP neuroimaging and accompanying template spaces to standardize the localization of neuroimaging results. Our group recently developed the NIMH Macaque Template (NMT), a high-resolution population average anatomical template and associated neuroimaging resources, providing researchers with a standard space for macaque neuroimaging . Here, we release NMT v2, which includes both symmetric and asymmetric templates in stereotaxic orientation, with improvements in spatial contrast, processing efficiency, and segmentation. We also introduce the Cortical Hierarchy Atlas of the Rhesus Macaque (CHARM), a hierarchical parcellation of the macaque cerebral cortex with varying degrees of detail. These tools have been integrated into the neuroimaging analysis software AFNI to provide a comprehensive and robust pipeline for fMRI processing, visualization and analysis of NHP data. AFNI's new @animal_warper program can be used to efficiently align anatomical scans to the NMT v2 space, and afni_proc.py integrates these results with full fMRI processing using macaque-specific parameters: from motion correction through regression modeling. Taken together, the NMT v2 and AFNI represent an all-in-one package for macaque functional neuroimaging analysis, as demonstrated with available demos for both task and resting state fMRI.