A comprehensive macaque fMRI pipeline and hierarchical atlas

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

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 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.

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 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.

Frontal Cortical Functional Connectivity Is Impacted by Anaesthesia in Macaques

A critical aspect of neuroscience is to establish whether and how brain networks evolved across primates. To date, most comparative studies have used resting-state functional magnetic resonance imaging (rs-fMRI) in anaesthetized nonhuman primates and in awake humans. However, anaesthesia strongly affects rs-fMRI signals. The present study investigated the impact of the awareness state (anaesthesia vs. awake) within the same group of macaque monkeys on the rs-fMRI functional connectivity organization of a well-characterized network in the human brain, the cingulo-frontal lateral network. Results in awake macaques show that rostral seeds in the cingulate sulcus exhibited stronger correlation strength with rostral compared to caudal lateral frontal cortical areas, while more caudal seeds displayed stronger correlation strength with caudal compared to anterior lateral frontal cortical areas. Critically, this inverse rostro-caudal functional gradient was abolished under anaesthesia. This study demonstrated a similar functional connectivity (FC) organization of the cingulo-frontal cortical network in awake macaque to that previously uncovered in the human brain pointing toward a preserved FC organization from macaque to human. However, it can only be observed in awake state suggesting that this network is sensitive to anaesthesia and warranting significant caution when comparing FC patterns across species under different states.

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.

Combining brain perturbation and neuroimaging in non-human primates

Brain perturbation studies allow detailed causal inferences of behavioral and neural processes. Because the combination of brain perturbation methods and neural measurement techniques is inherently challenging, research in humans has predominantly focused on non-invasive, indirect brain perturbations, or neurological lesion studies. Non-human primates have been indispensable as a neurobiological system that is highly similar to humans while simultaneously being more experimentally tractable, allowing visualization of the functional and structural impact of systematic brain perturbation. This review considers the state of the art in non-human primate brain perturbation with a focus on approaches that can be combined with neuroimaging. We consider both non-reversible (lesions) and reversible or temporary perturbations such as electrical, pharmacological, optical, optogenetic, chemogenetic, pathway-selective, and ultrasound based interference methods. Method-specific considerations from the research and development community are offered to facilitate research in this field and support further innovations. We conclude by identifying novel avenues for further research and innovation and by highlighting the clinical translational potential of the methods.