A comprehensive macaque fMRI pipeline and hierarchical atlas

A semi-automated pipeline for finite element modeling of electric field induced in nonhuman primates by transcranial magnetic stimulation

BACKGROUND

Transcranial magnetic stimulation (TMS) is used to treat a range of brain disorders by inducing an electric field (E-field) in the brain. However, the precise neural effects of TMS are not well understood. Nonhuman primates (NHPs) are used to model the impact of TMS on neural activity, but a systematic method of quantifying the induced E-field in the cortex of NHPs has not been developed.

NEW METHOD

The pipeline uses statistical parametric mapping (SPM) to automatically segment a structural MRI image of a rhesus macaque into five tissue compartments. Manual corrections are necessary around implants. The segmented tissues are tessellated into 3D meshes used in finite element method (FEM) software to compute the TMS induced E-field in the brain. The gray matter can be further segmented into cortical laminae using a volume preserving method for defining layers.

RESULTS

Models of three NHPs were generated with TMS coils placed over the precentral gyrus. Two coil configurations, active and sham, were simulated and compared. The results demonstrated a large difference in E-fields at the target. Additionally, the simulations were calculated using two different E-field solvers and were found to not significantly differ.

COMPARISON WITH EXISTING METHODS

Current methods segment NHP tissues manually or use automated methods for only the brain tissue. Existing methods also do not stratify the gray matter into layers.

CONCLUSION

The pipeline calculates the induced E-field in NHP models by TMS and can be used to plan implant surgeries and determine approximate E-field values around neuron recording sites.

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.

A Set of FMRI Quality Control Tools in AFNI: Systematic, in-depth and interactive QC with afni_proc.py and more

Quality control (QC) assessment is a vital part of FMRI processing and analysis, and a typically under-discussed aspect of reproducibility. This includes checking datasets at their very earliest stages (acquisition and conversion) through their processing steps (e.g., alignment and motion correction) to regression modeling (correct stimuli, no collinearity, valid fits, enough degrees of freedom, etc.) for each subject. There are a wide variety of features to verify throughout any single subject processing pipeline, both quantitatively and qualitatively. We present several FMRI preprocessing QC features available in the AFNI toolbox, many of which are automatically generated by the pipeline-creation tool, afni_proc.py. These items include: a modular HTML document that covers full single subject processing from the raw data through statistical modeling; several review scripts in the results directory of processed data; and command line tools for identifying subjects with one or more quantitative properties across a group (such as triaging warnings, making exclusion criteria or creating informational tables). The HTML itself contains several buttons that efficiently facilitate interactive investigations into the data, when deeper checks are needed beyond the systematic images. The pages are linkable, so that users can evaluate individual items across a group, for increased sensitivity to differences (e.g., in alignment or regression modeling images). Finally, the QC document contains rating buttons for each "QC block", as well as comment fields for each, to facilitate both saving and sharing the evaluations. This increases the specificity of QC, as well as its shareability, as these files can be shared with others and potentially uploaded into repositories, promoting transparency and open science. We describe the features and applications of these QC tools for FMRI.

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.

Using camera-guided electrode microdrive navigation for precise 3D targeting of macaque brain sites

Spatial accuracy in electrophysiological investigations is paramount, as precise localization and reliable access to specific brain regions help the advancement of our understanding of the brain's complex neural activity. Here, we introduce a novel, multi camera-based, frameless neuronavigation technique for precise, 3-dimensional electrode positioning in awake monkeys. The investigation of neural functions in awake primates often requires stable access to the brain with thin and delicate recording electrodes. This is usually realized by implanting a chronic recording chamber onto the skull of the animal that allows direct access to the dura. Most recording and positioning techniques utilize this implanted recording chamber as a holder of the microdrive or to hold a grid. This in turn reduces the degrees of freedom in positioning. To solve this problem, we require innovative, flexible, but precise tools for neuronal recordings. We instead mount the electrode microdrive above the animal on an arch, equipped with a series of translational and rotational micromanipulators, allowing movements in all axes. Here, the positioning is controlled by infrared cameras tracking the location of the microdrive and the monkey, allowing precise and flexible trajectories. To verify the accuracy of this technique, we created iron deposits in the tissue that could be detected by MRI. Our results demonstrate a remarkable precision with the confirmed physical location of these deposits averaging less than 0.5 mm from their planned position. Pilot electrophysiological recordings additionally demonstrate the accuracy and flexibility of this method. Our innovative approach could significantly enhance the accuracy and flexibility of neural recordings, potentially catalyzing further advancements in neuroscientific research.

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.

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.

The Subcortical Atlas of the Marmoset ("SAM") monkey based on high-resolution MRI and histology

A comprehensive three-dimensional digital brain atlas of cortical and subcortical regions based on MRI and histology has a broad array of applications in anatomical, functional, and clinical studies. We first generated a Subcortical Atlas of the Marmoset, called the "SAM," from 251 delineated subcortical regions (e.g. thalamic subregions, etc.) derived from high-resolution Mean Apparent Propagator-MRI, T2W, and magnetization transfer ratio images ex vivo. We then confirmed the location and borders of these segmented regions in the MRI data using matched histological sections with multiple stains obtained from the same specimen. Finally, we estimated and confirmed the atlas-based areal boundaries of subcortical regions by registering this ex vivo atlas template to in vivo T1- or T2W MRI datasets of different age groups (single vs. multisubject population-based marmoset control adults) using a novel pipeline developed within Analysis of Functional NeuroImages software. Tracing and validating these important deep brain structures in 3D will improve neurosurgical planning, anatomical tract tracer injections, navigation of deep brain stimulation probes, functional MRI and brain connectivity studies, and our understanding of brain structure-function relationships. This new ex vivo template and atlas are available as volumes in standard NIFTI and GIFTI file formats and are intended for use as a reference standard for marmoset brain research.

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.

Prefrontal-Limbic Circuitry Is Associated With Reward Sensitivity in Nonhuman Primates

BACKGROUND

Abnormal reward sensitivity is a risk factor for psychiatric disorders, including eating disorders such as overeating and binge-eating disorder, but the brain structural mechanisms that underlie it are not completely understood. Here, we sought to investigate the relationship between multimodal whole-brain structural features and reward sensitivity in nonhuman primates.

METHODS

Reward sensitivity was evaluated through behavioral economic analysis in which monkeys (adult rhesus macaques; 7 female, 5 male) responded for sweetened condensed milk (10%, 30%, 56%), Gatorade, or water using an operant procedure in which the response requirement increased incrementally across sessions (i.e., fixed ratio 1, 3, 10). Animals were divided into high (n = 6) or low (n = 6) reward sensitivity groups based on essential value for 30% milk. Multimodal magnetic resonance imaging was used to measure gray matter volume and white matter microstructure. Brain structural features were compared between groups, and their correlations with reward sensitivity for various stimuli was investigated.

RESULTS

Animals in the high sensitivity group had greater dorsolateral prefrontal cortex, centromedial amygdaloid complex, and middle cingulate cortex volumes than animals in the low sensitivity group. Furthermore, compared with monkeys in the low sensitivity group, high sensitivity monkeys had lower fractional anisotropy in the left dorsal cingulate bundle connecting the centromedial amygdaloid complex and middle cingulate cortex to the dorsolateral prefrontal cortex, and in the left superior longitudinal fasciculus 1 connecting the middle cingulate cortex to the dorsolateral prefrontal cortex.

CONCLUSIONS

These results suggest that neuroanatomical variation in prefrontal-limbic circuitry is associated with reward sensitivity. These brain structural features may serve as predictive biomarkers for vulnerability to food-based and other reward-related disorders.

Spontaneous Behavioural Recovery Following Stroke Relates to the Integrity of Parietal and Temporal Regions

Stroke is a devastating disease that results in neurological deficits and represents a leading cause of death and disability worldwide. Following a stroke, there is a degree of spontaneous recovery of function, the neural basis of which is of great interest among clinicians in their efforts to reduce disability following stroke and enhance rehabilitation. Conventionally, work on spontaneous recovery has tended to focus on the neural reorganization of motor cortical regions, with comparably little attention being paid to changes in non-motor regions and how these relate to recovery. Here we show, using structural neuroimaging in a macaque stroke model (N = 31) and by exploiting individual differences in spontaneous behavioural recovery, that the preservation of regions in the parietal and temporal cortices predict animal recovery. To characterize recovery, we performed a clustering analysis using Non-Human Primate Stroke Scale (NHPSS) scores and identified a good versus poor recovery group. By comparing the preservation of brain volumes in the two groups, we found that brain areas in integrity of brain areas in parietal, temporal and somatosensory cortex were associated with better recovery. In addition, a decoding approach performed across all subjects revealed that the preservation of specific brain regions in the parietal, somatosensory and medial frontal cortex predicted recovery. Together, these findings highlight the importance of parietal and temporal regions in spontaneous behavioural recovery.

The neural origin for asymmetric coding of surface color in the primate visual cortex

The coding privilege of end-spectral hues (red and blue) in the early visual cortex has been reported in primates. However, the origin of such bias remains unclear. Here, we provide a complete picture of the end-spectral bias in visual system by measuring fMRI signals and spiking activities in macaques. The correlated end-spectral biases between the LGN and V1 suggest a subcortical source for asymmetric coding. Along the ventral pathway from V1 to V4, red bias against green peaked in V1 and then declined, whereas blue bias against yellow showed an increasing trend. The feedforward and recurrent modifications of end-spectral bias were further revealed by dynamic causal modeling analysis. Moreover, we found that the strongest end-spectral bias in V1 was in layer 4C[Formula: see text]. Our results suggest that end-spectral bias already exists in the LGN and is transmitted to V1 mainly through the parvocellular pathway, then embellished by cortical processing.

The Subcortical Atlas of the Marmoset ("SAM") monkey based on high-resolution MRI and histology

A comprehensive three-dimensional digital brain atlas of cortical and subcortical regions based on MRI and histology has a broad array of applications for anatomical, functional, and clinical studies. We first generated a Subcortical Atlas of the Marmoset, called the "SAM," from 251 delineated subcortical regions (e.g., thalamic subregions, etc.) derived from the high-resolution MAP-MRI, T2W, and MTR images ex vivo. We then confirmed the location and borders of these segmented regions in MRI data using matched histological sections with multiple stains obtained from the same specimen. Finally, we estimated and confirmed the atlas-based areal boundaries of subcortical regions by registering this ex vivo atlas template to in vivo T1- or T2W MRI datasets of different age groups (single vs. multisubject population-based marmoset control adults) using a novel pipeline developed within AFNI. Tracing and validating these important deep brain structures in 3D improves neurosurgical planning, anatomical tract tracer injections, navigation of deep brain stimulation probes, fMRI and brain connectivity studies, and our understanding of brain structure-function relationships. This new ex vivo template and atlas are available as volumes in standard NIFTI and GIFTI file formats and are intended for use as a reference standard for marmoset brain research.

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.

[Whole-brain parcellation for macaque brain magnetic resonance images based on attention mechanism and multi-modality feature fusion]

OBJECTIVE

We propose a novel deep learning algorithm based on attention mechanism and multimodality feature fusion (DDAM) to achieve whole-brain parcellation of macaque brain magnetic resonance (MR) images.

METHODS

We collected multimodality brain MR images of 68 macaques (aged 13 to 36 months) with corresponding ground truth labels. To address the complexity and complementary nature of the multimodality data, we employed a multi- encoder structure to adapt to different modalities and performed feature extraction. In the decoder, an attention mechanism was introduced to construct the Attention-based Multimodality Feature Fusion module (AMFF). Leveraging the rich and complementary information between modalities, AMFF effectively integrated multimodality features of varying scales and complexities to enhance the performance of parcellation. We conducted ablation experiments and statistically analyzed the results.

RESULTS

The incorporation of the multi- encoder structure and attention mechanism significantly improved the performance of the network for integrating the multimodality features, and achieved an average DSC of 0.904 and an ASD as low as 0.131 for macaque whole-brain parcellation (P < 0.05). The ablation experiments validated the effectiveness of each component of the DDAM method.

CONCLUSION

The proposed DDAM method, with its enhanced ability to extract and integrate multimodality features, can significantly improve the accuracy of whole-brain MR image segmentation.

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.

Cortical origin of theta error signals

A multi-scale approach elucidated the origin of the error-related-negativity (ERN), with its associated theta-rhythm, and the post-error-positivity (Pe) in macaque supplementary eye field (SEF). Using biophysical modeling, synaptic inputs to a subpopulation of layer-3 (L3) and layer-5 (L5) pyramidal cells (PCs) were optimized to reproduce error-related spiking modulation and inter-spike intervals. The intrinsic dynamics of dendrites in L5 but not L3 error PCs generate theta rhythmicity with random phases. Saccades synchronized the phases of the theta-rhythm, which was magnified on errors. Contributions from error PCs to the laminar current source density (CSD) observed in SEF were negligible and could not explain the observed association between error-related spiking modulation in L3 PCs and scalp-EEG. CSD from recorded laminar field potentials in SEF was comprised of multipolar components, with monopoles indicating strong electro-diffusion, dendritic/axonal electrotonic current leakage outside SEF, or violations of the model assumptions. Our results also demonstrate the involvement of secondary cortical regions, in addition to SEF, particularly for the later Pe component. The dipolar component from the observed CSD paralleled the ERN dynamics, while the quadrupolar component paralleled the Pe. These results provide the most advanced explanation to date of the cellular mechanisms generating the ERN.

General anaesthesia reduces the uniqueness of brain connectivity across individuals and across species

The human brain is characterised by idiosyncratic patterns of spontaneous thought, rendering each brain uniquely identifiable from its neural activity. However, deep general anaesthesia suppresses subjective experience. Does it also suppress what makes each brain unique? Here we used functional MRI under the effects of the general anaesthetics sevoflurane and propofol to determine whether anaesthetic-induced unconsciousness diminishes the uniqueness of the human brain: both with respect to the brains of other individuals, and the brains of another species. We report that under anaesthesia individual brains become less self-similar and less distinguishable from each other. Loss of distinctiveness is highly organised: it co-localises with the archetypal sensory-association axis, correlating with genetic and morphometric markers of phylogenetic differences between humans and other primates. This effect is more evident at greater anaesthetic depths, reproducible across sevoflurane and propofol, and reversed upon recovery. Providing convergent evidence, we show that under anaesthesia the functional connectivity of the human brain becomes more similar to the macaque brain. Finally, anaesthesia diminishes the match between spontaneous brain activity and meta-analytic brain patterns aggregated from the NeuroSynth engine. Collectively, the present results reveal that anaesthetised human brains are not only less distinguishable from each other, but also less distinguishable from the brains of other primates, with specifically human-expanded regions being the most affected by anaesthesia.

Human-to-monkey transfer learning identifies the frontal white matter as a key determinant for predicting monkey brain age

The application of artificial intelligence (AI) to summarize a whole-brain magnetic resonance image (MRI) into an effective "brain age" metric can provide a holistic, individualized, and objective view of how the brain interacts with various factors (e.g., genetics and lifestyle) during aging. Brain age predictions using deep learning (DL) have been widely used to quantify the developmental status of human brains, but their wider application to serve biomedical purposes is under criticism for requiring large samples and complicated interpretability. Animal models, i.e., rhesus monkeys, have offered a unique lens to understand the human brain - being a species in which aging patterns are similar, for which environmental and lifestyle factors are more readily controlled. However, applying DL methods in animal models suffers from data insufficiency as the availability of animal brain MRIs is limited compared to many thousands of human MRIs. We showed that transfer learning can mitigate the sample size problem, where transferring the pre-trained AI models from 8,859 human brain MRIs improved monkey brain age estimation accuracy and stability. The highest accuracy and stability occurred when transferring the 3D ResNet [mean absolute error (MAE) = 1.83 years] and the 2D global-local transformer (MAE = 1.92 years) models. Our models identified the frontal white matter as the most important feature for monkey brain age predictions, which is consistent with previous histological findings. This first DL-based, anatomically interpretable, and adaptive brain age estimator could broaden the application of AI techniques to various animal or disease samples and widen opportunities for research in non-human primate brains across the lifespan.

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.

Small volume blood-brain barrier opening in macaques with a 1 MHz ultrasound phased array

The use of focused ultrasound to open the blood-brain barrier (BBB) has the potential to deliver drugs to specific regions of the brain. The size of the BBB opening and ability to localize the opening determines the spatial extent and is a limiting factor in many applications of BBB opening where targeting a small brain region is desired. Here we evaluate the performance of a system designed for small opening volumes and highlight the unique challenges associated with pushing the spatial precision of this technique. To achieve small volume openings in cortical regions of the macaque brain, we tested a custom 1 MHz array transducer integrated into a magnetic resonance image-guided focused ultrasound system. Using real-time cavitation monitoring, we demonstrated twelve instances of single sonication, small volume BBB opening with average volumes of 59 ± 37 mm3 and 184 ± 2 mm3 in cortical and subcortical targets, respectively. We found high correlation between subject-specific acoustic simulations and observed openings when incorporating grey matter segmentation (R2 = 0.8577), and the threshold for BBB opening based on simulations was 0.53 MPa. Analysis of MRI-based safety assessment and cavitation signals indicate a safe pressure range for 1 MHz BBB opening and suggest that our system can be used to deliver drugs and gene therapy to small brain regions.

High-dimensional topographic organization of visual features in the primate temporal lobe

The inferotemporal cortex supports our supreme object recognition ability. Numerous studies have been conducted to elucidate the functional organization of this brain area, but there are still important questions that remain unanswered, including how this organization differs between humans and non-human primates. Here, we use deep neural networks trained on object categorization to construct a 25-dimensional space of visual features, and systematically measure the spatial organization of feature preference in both male monkey brains and human brains using fMRI. These feature maps allow us to predict the selectivity of a previously unknown region in monkey brains, which is corroborated by additional fMRI and electrophysiology experiments. These maps also enable quantitative analyses of the topographic organization of the temporal lobe, demonstrating the existence of a pair of orthogonal gradients that differ in spatial scale and revealing significant differences in the functional organization of high-level visual areas between monkey and human brains.

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.

Disrupting nociceptive information processing flow through transcranial focused ultrasound neuromodulation of thalamic nuclei

BACKGROUND

MRI-guided transcranial focused ultrasound (MRgFUS) as a next-generation neuromodulation tool can precisely target and stimulate deep brain regions with high spatial selectivity. Combined with MR-ARFI (acoustic radiation force imaging) and using fMRI BOLD signal as functional readouts, our previous studies have shown that low-intensity FUS can excite or suppress neural activity in the somatosensory cortex.

OBJECTIVE

To investigate whether low-intensity FUS can suppress nociceptive heat stimulation-induced responses in thalamic nuclei during hand stimulation, and to determine how this suppression influences the information processing flow within nociception networks.

FINDINGS

BOLD fMRI activations evoked by 47.5 °C heat stimulation of hand were detected in 24 cortical regions, which belong to sensory, affective, and cognitive nociceptive networks. Concurrent delivery of low-intensity FUS pulses (650 kHz, 550 kPa) to the predefined heat nociceptive stimulus-responsive thalamic centromedial_parafascicular (CM_para), mediodorsal (MD), ventral_lateral (VL_ and ventral_lateral_posteroventral (VLpv) nuclei suppressed their heat responses. Off-target cortical areas exhibited reduced, enhanced, or no significant fMRI signal changes, depending on the specific areas. Differentiable thalamocortical information flow during the processing of nociceptive heat input was observed, as indicated by the time to reach 10% or 30% of the heat-evoked BOLD signal peak. Suppression of thalamic heat responses significantly altered nociceptive processing flow and direction between the thalamus and cortical areas. Modulation of contralateral versus ipsilateral areas by unilateral thalamic activity differed. Signals detected in high-order cortical areas, such as dorsal frontal (DFC) and ventrolateral prefrontal (vlPFC) cortices, exhibited faster response latencies than sensory areas.

CONCLUSIONS

The concurrent delivery of FUS suppressed nociceptive heat response in thalamic nuclei and disrupted the nociceptive network. This study offers new insights into the causal functional connections within the thalamocortical networks and demonstrates the modulatory effects of low-intensity FUS on nociceptive information processing.

versaFlow: a versatile pipeline for resolution adapted diffusion MRI processing and its application to studying the variability of the PRIME-DE database

The lack of "gold standards" in Diffusion Weighted Imaging (DWI) makes validation cumbersome. To tackle this task, studies use translational analysis where results in humans are benchmarked against findings in other species. Non-Human Primates (NHP) are particularly interesting for this, as their cytoarchitecture is closely related to humans. However, tools used for processing and analysis must be adapted and finely tuned to work well on NHP images. Here, we propose versaFlow, a modular pipeline implemented in Nextflow, designed for robustness and scalability. The pipeline is tailored to in vivo NHP DWI at any spatial resolution; it allows for maintainability and customization. Processes and workflows are implemented using cutting-edge and state-of-the-art Magnetic Resonance Imaging (MRI) processing technologies and diffusion modeling algorithms, namely Diffusion Tensor Imaging (DTI), Constrained Spherical Deconvolution (CSD), and DIstribution of Anisotropic MicrOstructural eNvironments in Diffusion-compartment imaging (DIAMOND). Using versaFlow, we provide an in-depth study of the variability of diffusion metrics computed on 32 subjects from 3 sites of the Primate Data Exchange (PRIME-DE), which contains anatomical T1-weighted (T1w) and T2-weighted (T2w) images, functional MRI (fMRI), and DWI of NHP brains. This dataset includes images acquired over a range of resolutions, using single and multi-shell gradient samplings, on multiple scanner vendors. We perform a reproducibility study of the processing of versaFlow using the Aix-Marseilles site's data, to ensure that our implementation has minimal impact on the variability observed in subsequent analyses. We report very high reproducibility for the majority of metrics; only gamma distribution parameters of DIAMOND display less reproducible behaviors, due to the absence of a mechanism to enforce a random number seed in the software we used. This should be taken into consideration when future applications are performed. We show that the PRIME-DE diffusion data exhibits a great level of variability, similar or greater than results obtained in human studies. Its usage should be done carefully to prevent instilling uncertainty in statistical analyses. This hints at a need for sufficient harmonization in acquisition protocols and for the development of robust algorithms capable of managing the variability induced in imaging due to differences in scanner models and/or vendors.

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.

Chronic stability of a neuroprosthesis comprising multiple adjacent Utah arrays in monkeys

Objective. Electrical stimulation of visual cortex via a neuroprosthesis induces the perception of dots of light ('phosphenes'), potentially allowing recognition of simple shapes even after decades of blindness. However, restoration of functional vision requires large numbers of electrodes, and chronic, clinical implantation of intracortical electrodes in the visual cortex has only been achieved using devices of up to 96 channels. We evaluated the efficacy and stability of a 1024-channel neuroprosthesis system in non-human primates (NHPs) over more than 3 years to assess its suitability for long-term vision restoration.Approach.We implanted 16 microelectrode arrays (Utah arrays) consisting of 8 × 8 electrodes with iridium oxide tips in the primary visual cortex (V1) and visual area 4 (V4) of two sighted macaques. We monitored the animals' health and measured electrode impedances and neuronal signal quality by calculating signal-to-noise ratios of visually driven neuronal activity, peak-to-peak voltages of the waveforms of action potentials, and the number of channels with high-amplitude signals. We delivered cortical microstimulation and determined the minimum current that could be perceived, monitoring the number of channels that successfully yielded phosphenes. We also examined the influence of the implant on a visual task after 2-3 years of implantation and determined the integrity of the brain tissue with a histological analysis 3-3.5 years post-implantation.Main results. The monkeys remained healthy throughout the implantation period and the device retained its mechanical integrity and electrical conductivity. However, we observed decreasing signal quality with time, declining numbers of phosphene-evoking electrodes, decreases in electrode impedances, and impaired performance on a visual task at visual field locations corresponding to implanted cortical regions. Current thresholds increased with time in one of the two animals. The histological analysis revealed encapsulation of arrays and cortical degeneration. Scanning electron microscopy on one array revealed degradation of IrOxcoating and higher impedances for electrodes with broken tips.Significance. Long-term implantation of a high-channel-count device in NHP visual cortex was accompanied by deformation of cortical tissue and decreased stimulation efficacy and signal quality over time. We conclude that improvements in device biocompatibility and/or refinement of implantation techniques are needed before future clinical use is feasible.

A revised perspective on the evolution of the lateral frontal cortex in primates

Detailed neuroscientific data from macaque monkeys have been essential in advancing understanding of human frontal cortex function, particularly for regions of frontal cortex without homologs in other model species. However, precise transfer of this knowledge for direct use in human applications requires an understanding of monkey to hominid homologies, particularly whether and how sulci and cytoarchitectonic regions in the frontal cortex of macaques relate to those in hominids. We combine sulcal pattern analysis with resting-state functional magnetic resonance imaging and cytoarchitectonic analysis to show that old-world monkey brains have the same principles of organization as hominid brains, with the notable exception of sulci in the frontopolar cortex. This essential comparative framework provides insights into primate brain evolution and a key tool to drive translation from invasive research in monkeys to human applications.

Monkey Dorsolateral Prefrontal Cortex Represents Abstract Visual Sequences during a No-Report Task

Monitoring sequential information is an essential component of our daily lives. Many of these sequences are abstract, in that they do not depend on the individual stimuli, but do depend on an ordered set of rules (e.g., chop then stir when cooking). Despite the ubiquity and utility of abstract sequential monitoring, little is known about its neural mechanisms. Human rostrolateral prefrontal cortex (RLPFC) exhibits specific increases in neural activity (i.e., "ramping") during abstract sequences. Monkey dorsolateral prefrontal cortex (DLPFC) has been shown to represent sequential information in motor (not abstract) sequence tasks, and contains a subregion, area 46, with homologous functional connectivity to human RLPFC. To test the prediction that area 46 may represent abstract sequence information, and do so with parallel dynamics to those found in humans, we conducted functional magnetic resonance imaging (fMRI) in three male monkeys. When monkeys performed no-report abstract sequence viewing, we found that left and right area 46 responded to abstract sequential changes. Interestingly, responses to rule and number changes overlapped in right area 46 and left area 46 exhibited responses to abstract sequence rules with changes in ramping activation, similar to that observed in humans. Together, these results indicate that monkey DLPFC monitors abstract visual sequential information, potentially with a preference for different dynamics in the two hemispheres. More generally, these results show that abstract sequences are represented in functionally homologous regions across monkeys and humans.SIGNIFICANCE STATEMENT Daily, we complete sequences that are "abstract" because they depend on an ordered set of rules (e.g., chop then stir when cooking) rather than the identity of individual items. Little is known about how the brain tracks, or monitors, this abstract sequential information. Based on previous human work showing abstract sequence related dynamics in an analogous area, we tested whether monkey dorsolateral prefrontal cortex (DLPFC), specifically area 46, represents abstract sequential information using awake monkey functional magnetic resonance imaging (fMRI). We found that area 46 responded to abstract sequence changes, with a preference for more general responses on the right and dynamics similar to humans on the left. These results suggest that abstract sequences are represented in functionally homologous regions across monkeys and humans.

A novel methodology for exact targeting of human and non-human primate brain structures and skull implants using atlas-based 3D reconstruction

BACKGROUND

Accurate targeting of brain areas for stimulation and/or electrophysiological recording is key in many therapeutic applications and basic neuroscience research. Nevertheless, there are currently no end-to-end packages that accommodate all steps required for exact localization, visualization and targeting regions of interest (ROIs) using standard atlases and for designing skull implants.

NEW METHOD

We have implemented a new processing pipeline that addresses this issue in macaques and humans including various preprocessing, registration, warping procedures and 3D reconstructions and provide a noncommercial open-source graphical software which we refer to as the MATLAB-based reconstruction for recording and stimulation (MATres).

RESULTS

The results of skull stripping were shown to work seamlessly in humans and monkeys. Linear and nonlinear warping of the standard atlas to the native space outperformed state-of-the-art using AFNI with improvements being more prominent in humans which had a more complex gyration geometry. The skull surface extracted by MATres using MRI images had more than 90% match with CT ground truth and could be used to design skull implants that conformed well to the skull's local curvature.

COMPARISON WITH EXISTING METHOD(S)

The accuracy of the various steps including skull stripping, standard atlas registration and skull reconstruction in MATres was compared with and shown to outperform the AFNI. The localization accuracy of the recording chambers designed with MATres and implanted in two macaque monkeys was further confirmed using MRI imaging.

CONCLUSIONS

Precise localization of ROIs offered by MATres can be used for planning electrode penetrations for recording and for shallow or deep brain stimulation (DBS).

Brain-wide and cell-specific transcriptomic insights into MRI-derived cortical morphology in macaque monkeys

Integrative analyses of transcriptomic and neuroimaging data have generated a wealth of information about biological pathways underlying regional variability in imaging-derived brain phenotypes in humans, but rarely in nonhuman primates due to the lack of a comprehensive anatomically-defined atlas of brain transcriptomics. Here we generate complementary bulk RNA-sequencing dataset of 819 samples from 110 brain regions and single-nucleus RNA-sequencing dataset, and neuroimaging data from 162 cynomolgus macaques, to examine the link between brain-wide gene expression and regional variation in morphometry. We not only observe global/regional expression profiles of macaque brain comparable to human but unravel a dorsolateral-ventromedial gradient of gene assemblies within the primate frontal lobe. Furthermore, we identify a set of 971 protein-coding and 34 non-coding genes consistently associated with cortical thickness, specially enriched for neurons and oligodendrocytes. These data provide a unique resource to investigate nonhuman primate models of human diseases and probe cross-species evolutionary mechanisms.

Small volume blood-brain barrier opening in macaques with a 1 MHz ultrasound phased array

Focused ultrasound blood-brain barrier (BBB) opening is a promising tool for targeted delivery of therapeutic agents into the brain. The volume of opening determines the extent of therapeutic administration and sets a lower bound on the size of targets which can be selectively treated. We tested a custom 1 MHz array transducer optimized for cortical regions in the macaque brain with the goal of achieving small volume openings. We integrated this device into a magnetic resonance image guided focused ultrasound system and demonstrated twelve instances of small volume BBB opening with average opening volumes of 59 ± 37 mm ³ and 184 ± 2 mm ³ in cortical and subcortical targets, respectively. We developed real-time cavitation monitoring using a passive cavitation detector embedded in the array and characterized its performance on a bench-top flow phantom mimicking transcranial BBB opening procedures. We monitored cavitation during in-vivo procedures and compared cavitation metrics against opening volumes and safety outcomes measured with FLAIR and susceptibility weighted MR imaging. Our findings show small BBB opening at cortical targets in macaques and characterize the safe pressure range for 1 MHz BBB opening. Additionally, we used subject-specific simulations to investigate variance in measured opening volumes and found high correlation (R ² = 0.8577) between simulation predictions and observed measurements. Simulations suggest the threshold for 1 MHz BBB opening was 0.53 MPa. This system enables BBB opening for drug delivery and gene therapy to be targeted to more specific brain regions.

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.

HC-Net: A hybrid convolutional network for non-human primate brain extraction

Brain extraction (skull stripping) is an essential step in the magnetic resonance imaging (MRI) analysis of brain sciences. However, most of the current brain extraction methods that achieve satisfactory results for human brains are often challenged by non-human primate brains. Due to the small sample characteristics and the nature of thick-slice scanning of macaque MRI data, traditional deep convolutional neural networks (DCNNs) are unable to obtain excellent results. To overcome this challenge, this study proposed a symmetrical end-to-end trainable hybrid convolutional neural network (HC-Net). It makes full use of the spatial information between adjacent slices of the MRI image sequence and combines three consecutive slices from three axes for 3D convolutions, which reduces the calculation consumption and promotes accuracy. The HC-Net consists of encoding and decoding structures of 3D convolutions and 2D convolutions in series. The effective use of 2D convolutions and 3D convolutions relieves the underfitting of 2D convolutions to spatial features and the overfitting of 3D convolutions to small samples. After evaluating macaque brain data from different sites, the results showed that HC-Net performed better in inference time (approximately 13 s per volume) and accuracy (mean Dice coefficient reached 95.46%). The HC-Net model also had good generalization ability and stability in different modes of brain extraction tasks.

Utilizing PET and MALDI Imaging for Discovery of a Targeted Probe for Brain Endocannabinoid α/β-Hydrolase Domain 6 (ABHD6)

Multimodal imaging provides rich biological information, which can be exploited to study drug activity, disease associated phenotypes, and pharmacological responses. Here we show discovery and validation of a new probe targeting the endocannabinoid α/β-hydrolase domain 6 (ABHD6) enzyme by utilizing positron emission tomography (PET) and matrix-assisted laser desorption/ionization (MALDI) imaging. [¹⁸F]JZP-MA-11 as the first PET ligand for in vivo imaging of the ABHD6 is reported and specific uptake in ABHD6-rich peripheral tissues and major brain regions was demonstrated using PET. A proof-of-concept study in nonhuman primate confirmed brain uptake. In vivo pharmacological response upon ABHD6 inhibition was observed by MALDI imaging. These synergistic imaging efforts used to identify biological information cannot be obtained by a single imaging modality and hold promise for improving the understanding of ABHD6-mediated endocannabinoid metabolism in peripheral and central nervous system disorders.

Mapping blood traits to structural organization of the brain in rhesus monkeys

Hematological and biochemical blood traits have been linked to brain structural characteristics in humans. However, the relationship between these two domains has not been systematically explored in nonhuman primates, which are crucial animal models for understanding the mechanisms of brain function and developing therapeutics for various disorders. Here we investigated the associations between hematological/biochemical parameters and the brain's gray matter volume and white matter integrity derived from T1-weighted and diffusion magnetic resonance imaging in 36 healthy macaques. We found that intersubject variations in basophil count and hemoglobin levels correlated with gray matter volumes in the anterior cingulum, prefrontal cortex, and putamen. Through interactions between these key elements, the blood parameters' covariation network was linked with that of the brain structures, forming overarching networks connecting blood traits with structural brain features. These networks exhibited hierarchical small-world architecture, indicating highly effective interactions between their constituent elements. In addition, different subnetworks of the brain areas or fiber tracts tended to correlate with unique groups of blood indices, revealing previously unknown brain structural organization. These results provide a quantitative characterization of the interactions between blood parameters and brain structures in macaques and may increase the understanding of the body-brain relationship and the pathogenesis of relevant disorders.

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.

Brain-wide functional connectivity of face patch neurons during rest

In resting state functional magnetic resonance imaging (fMRI), areas showing coherent hemodynamic fluctuations across the brain are operationally defined to be functionally connected. However, it is unknown how the activity of single units residing within a voxel contributes to this network structure. Here we demonstrate a shared but restricted pattern of functional connectivity among neighboring neurons residing in functionally defined face patches. Unexpectedly, such neurons also exhibited a prominent inverse correlation with thalamic structures and brainstem neuromodulatory centers. Single unit maps differed from analogous maps obtained with local field potentials and seed-based fMRI. These findings suggest that during rest, individual cortical neurons have a restricted set of functional connections, which is governed in part by anatomical projections and in part by neuromodulation.

The brain is a highly organized, dynamic system whose network architecture is often assessed through resting functional magnetic resonance imaging (fMRI) functional connectivity. The functional interactions between brain areas, including those observed during rest, are assumed to stem from the collective influence of action potentials carried by long-range neural projections. However, the contribution of individual neurons to brain-wide functional connectivity has not been systematically assessed. Here we developed a method to concurrently measure and compare the spiking activity of local neurons with fMRI signals measured across the brain during rest. We recorded spontaneous activity from neural populations in cortical face patches in the macaque during fMRI scanning sessions. Individual cells exhibited prominent, bilateral coupling with fMRI fluctuations in a restricted set of cortical areas inside and outside the face patch network, partially matching the pattern of known anatomical projections. Within each face patch population, a subset of neurons was positively coupled with the face patch network and another was negatively coupled. The same cells showed inverse correlations with distinct subcortical structures, most notably the lateral geniculate nucleus and brainstem neuromodulatory centers. Corresponding connectivity maps derived from fMRI seeds and local field potentials differed from the single unit maps, particularly in subcortical areas. Together, the results demonstrate that the spiking fluctuations of neurons are selectively coupled with discrete brain regions, with the coupling governed in part by anatomical network connections and in part by indirect neuromodulatory pathways.

Resolving the mesoscopic missing link: Biophysical modeling of EEG from cortical columns in primates

Event-related potentials (ERP) are among the most widely measured indices for studying human cognition. While their timing and magnitude provide valuable insights, their usefulness is limited by our understanding of their neural generators at the circuit level. Inverse source localization offers insights into such generators, but their solutions are not unique. To address this problem, scientists have assumed the source space generating such signals comprises a set of discrete equivalent current dipoles, representing the activity of small cortical regions. Based on this notion, theoretical studies have employed forward modeling of scalp potentials to understand how changes in circuit-level dynamics translate into macroscopic ERPs. However, experimental validation is lacking because it requires in vivo measurements of intracranial brain sources. Laminar local field potentials (LFP) offer a mechanism for estimating intracranial current sources. Yet, a theoretical link between LFPs and intracranial brain sources is missing. Here, we present a forward modeling approach for estimating mesoscopic intracranial brain sources from LFPs and predict their contribution to macroscopic ERPs. We evaluate the accuracy of this LFP-based representation of brain sources utilizing synthetic laminar neurophysiological measurements and then demonstrate the power of the approach in vivo to clarify the source of a representative cognitive ERP component. To that end, LFP was measured across the cortical layers of visual area V4 in macaque monkeys performing an attention demanding task. We show that area V4 generates dipoles through layer-specific transsynaptic currents that biophysically recapitulate the ERP component through the detailed forward modeling. The constraints imposed on EEG production by this method also revealed an important dissociation between computational and biophysical contributors. As such, this approach represents an important bridge between laminar microcircuitry, through the mesoscopic activity of cortical columns to the patterns of EEG we measure at the scalp.

Resting-State fMRI-Based Screening of Deschloroclozapine in Rhesus Macaques Predicts Dosage-Dependent Behavioral Effects

Chemogenetic techniques, such as designer receptors exclusively activated by designer drugs (DREADDs), enable transient, reversible, and minimally invasive manipulation of neural activity in vivo Their development in nonhuman primates is essential for uncovering neural circuits contributing to cognitive functions and their translation to humans. One key issue that has delayed the development of chemogenetic techniques in primates is the lack of an accessible drug-screening method. Here, we use resting-state fMRI, a noninvasive neuroimaging tool, to assess the impact of deschloroclozapine (DCZ) on brainwide resting-state functional connectivity in 7 rhesus macaques (6 males and 1 female) without DREADDs. We found that systemic administration of 0.1 mg/kg DCZ did not alter the resting-state functional connectivity. Conversely, 0.3 mg/kg of DCZ was associated with a prominent increase in functional connectivity that was mainly confined to the connections of frontal regions. Additional behavioral tests confirmed a negligible impact of 0.1 mg/kg DCZ on socio-emotional behaviors as well as on reaction time in a probabilistic learning task; 0.3 mg/kg DCZ did, however, slow responses in the probabilistic learning task, suggesting attentional or motivational deficits associated with hyperconnectivity in fronto-temporo-parietal networks. Our study highlights both the excellent selectivity of DCZ as a DREADD actuator, and the side effects of its excess dosage. The results demonstrate the translational value of resting-state fMRI as a drug-screening tool to accelerate the development of chemogenetics in primates.SIGNIFICANCE STATEMENT Chemogenetics, such as designer receptors exclusively activated by designer drugs (DREADDs), can afford control over neural activity with unprecedented spatiotemporal resolution. Accelerating the translation of chemogenetic neuromodulation from rodents to primates requires an approach to screen novel DREADD actuators in vivo Here, we assessed brainwide activity in response to a DREADD actuator deschloroclozapine (DCZ) using resting-state fMRI in macaque monkeys. We demonstrated that low-dose DCZ (0.1 mg/kg) did not change whole-brain functional connectivity or affective behaviors, while a higher dose (0.3 mg/kg) altered frontal functional connectivity and slowed response in a learning task. Our study highlights the excellent selectivity of DCZ at proper dosing, and demonstrates the utility of resting-state fMRI to screen novel chemogenetic actuators in primates.

A hierarchical anatomical framework and workflow for organizing stereotactic encephalography in epilepsy

Stereotactic electroencephalography (SEEG) is an increasingly utilized method for invasive monitoring in patients with medically intractable epilepsy. Yet, the lack of standardization for labeling electrodes hinders communication among clinicians. A rational clustering of contacts based on anatomy rather than arbitrary physical leads may help clinical neurophysiologists interpret seizure networks. We identified SEEG electrodes on post‐implant CTs and registered them to preoperative MRIs segmented according to an anatomical atlas. Individual contacts were automatically assigned to anatomical areas independent of lead. These contacts were then organized using a hierarchical anatomical schema for display and interpretation. Bipolar‐referenced signal cross‐correlations were used to compare the similarity of grouped signals within a conventional montage versus this anatomical montage. As a result, we developed a hierarchical organization for SEEG contacts using well‐accepted, free software that is based solely on their post‐implant anatomical location. When applied to three example SEEG cases for epilepsy, clusters of contacts that were anatomically related collapsed into standardized groups. Qualitatively, seizure events organized using this framework were better visually clustered compared to conventional schemes. Quantitatively, signals grouped by anatomical region were more similar to each other than electrode‐based groups as measured by Pearson correlation. Further, we uploaded visualizations of SEEG reconstructions into the electronic medical record, rendering them durably useful given the interpretable electrode labels. In conclusion, we demonstrate a standardized, anatomically grounded approach to the organization of SEEG neuroimaging and electrophysiology data that may enable improved communication among and across surgical epilepsy teams and promote a clearer view of individual seizure networks.

Selective Prefrontal-Amygdala Circuit Interactions Underlie Social and Nonsocial Valuation in Rhesus Macaques

Lesion studies in macaques suggest dissociable functions of the orbitofrontal cortex (OFC) and medial frontal cortex (MFC), with OFC being essential for goal-directed decision-making and MFC supporting social cognition. Bilateral amygdala damage results in impairments in both of these domains. There are extensive reciprocal connections between these prefrontal areas and the amygdala; however, it is not known whether the dissociable roles of OFC and MFC depend on functional interactions with the amygdala. To test this possibility, we compared the performance of male rhesus macaques (Macaca mulatta) with crossed surgical disconnection of the amygdala and either MFC (MFC × AMY, n = 4) or OFC (OFC × AMY, n = 4) to a group of unoperated controls (CON, n = 5). All monkeys were assessed for their performance on two tasks to measure the following: (1) food-retrieval latencies while viewing videos of social and nonsocial stimuli in a test of social interest and (2) object choices based on current food value using reinforcer devaluation in a test of goal-directed decision-making. Compared with the CON group, the MFC × AMY group, but not the OFC × AMY group, showed significantly reduced food-retrieval latencies while viewing videos of conspecifics, indicating reduced social valuation and/or interest. By contrast, on the devaluation task, group OFC × AMY, but not group MFC × AMY, displayed deficits on object choices following changes in food value. These data indicate that the MFC and OFC must functionally interact with the amygdala to support normative social and nonsocial valuation, respectively.SIGNIFICANCE STATEMENT Ascribing value to conspecifics (social) versus objects (nonsocial) may be supported by distinct but overlapping brain networks. Here, we test whether two nonoverlapping regions of the prefrontal cortex, the medial frontal cortex and the orbitofrontal cortex, must causally interact with the amygdala to sustain social valuation and goal-directed decision-making, respectively. We found that these prefrontal-amygdala circuits are functionally dissociable, lending support for the idea that medial frontal and orbital frontal cortex make independent contributions to cognitive appraisals of the environment. These data provide a neural framework for distinct value assignment processes and may enhance our understanding of the cognitive deficits observed following brain injury or in the development of mental health disorders.

Multilevel atlas comparisons reveal divergent evolution of the primate brain

Whether the size of the prefrontal cortex (PFC) in humans is disproportionate when compared to other species is a persistent debate in evolutionary neuroscience. This question has left the study of over/under-expansion in other structures relatively unexplored. We therefore sought to address this gap by adapting anatomical areas from the digital atlases of 18 mammalian species, to create a common interspecies classification. Our approach used data-driven analysis based on phylogenetic generalized least squares to evaluate anatomical expansion covering the whole brain. Our main finding suggests a divergence in primate evolution, orienting the stereotypical mammalian cerebral proportion toward a frontal and parietal lobe expansion in catarrhini (primate parvorder comprising old world monkeys, apes, and humans). Cerebral lobe volumes slopes plotted for catarrhini species were ranked as parietal∼frontal > temporal > occipital, contrasting with the ranking of other mammalian species (occipital > temporal > frontal∼parietal). Frontal and parietal slopes were statistically different in catarrhini when compared to other species through bootstrap analysis. Within the catarrhini's frontal lobe, the prefrontal cortex was the principal driver of frontal expansion. Across all species, expansion of the frontal lobe appeared to be systematically linked to the parietal lobe. Our findings suggest that the human frontal and parietal lobes are not disproportionately enlarged when compared to other catarrhini. Nevertheless, humans remain unique in carrying the most relatively enlarged frontal and parietal lobes in an infraorder exhibiting a disproportionate expansion of these areas.

Cerebro-cerebellar interactions in non-human 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.

The development of anesthesia was a significant advance in medicine. It allows individuals to undergo surgery without feeling pain or remembering the experience. But scientists still do not know how anesthesia works. During anesthesia, scientists have measured brain activity using electroencephalograms (EEG) and found that the brain appears to turn on and off. Comatose patients also have similar switches between bursts of electrical activity and periods of silence. This burst-suppression pattern may be related to unconsciousness.

But scientists still have many questions about how anesthesia causes burst-suppression. One challenge is that while an EEG can tell scientists when the brain turns on and off, it does not show exactly where this occurs. Another imaging method called functional Magnetic Resonance Imaging (fMRI) may fill this gap by allowing scientists to map where the brain activity occurs.

Sirmpilatze et al. have created detailed maps of burst-suppression in humans, primates, and rats under anesthesia by analyzing brain scans using fMRI. In rats, the entire outer layer or cortex of the brain underwent a synchronized pattern of burst-suppression. In humans and primates, areas of the brain like those responsible for eyesight did not follow the rest of the cortex in switching on and off.

The experiments reveal crucial differences in how rats and humans and other primates respond to anesthesia. The fMRI mapping technique Sirmpilatze et al. created may help scientists learn more about these differences and why some parts of human brains do not undergo burst-suppression. This may help scientists learn more about unconsciousness and help improve anesthesia or the care of comatose patients.

Early Alzheimer's disease-like reductions in gray matter and cognitive function with aging in nonhuman primates

Introduction

Age-related neuropathology associated with sporadic Alzheimer's disease (AD) often develops well before the onset of symptoms. Given AD's long preclinical period, translational models are needed to identify early signatures of pathological decline.

Methods

Using structural magnetic resonance imaging and cognitive assessments, we examined the relationships among age, cognitive performance, and neuroanatomy in 48 vervet monkeys (Chlorocebus aethiops sabaeus) ranging from young adults to very old.

Results

We found negative associations of age with cortical gray matter volume (P = .003) and the temporal-parietal cortical thickness meta-region of interest (P = .001). Additionally, cortical gray matter volumes predicted working memory at approximately 1-year follow-up (correct trials at the 20s delay [P = .008]; correct responses after longer delays [P = .004]).

Discussion

Cortical gray matter diminishes with age in vervets in regions relevant to AD, which may increase risk of cognitive impairment. This study lays the groundwork for future investigations to test therapeutics to delay or slow pathological decline.

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 fMRI performed at 10.5 Tesla 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 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 functional connectivity 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.